This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrowReprints and Permissions
Right arrow Copyright Information
Right arrow Books from ASM Press
Right arrow MicrobeWorld
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Pettersson, B.
Right arrow Articles by Johansson, K.-E.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Pettersson, B.
Right arrow Articles by Johansson, K.-E.

 Previous Article  |  Next Article 

J Bacteriol, May 1998, p. 2350-2358, Vol. 180, No. 9
0021-9193/98/$04.00+0
Copyright © 1998, American Society for Microbiology. All rights reserved.

Molecular Evolution of Mycoplasma capricolum subsp. capripneumoniae Strains, Based on Polymorphisms in the 16S rRNA Genes

Bertil Pettersson,1 Göran Bölske,2 François Thiaucourt,3 Mathias Uhlén,1 and Karl-Erik Johansson2,4,*

Department of Biochemistry and Biotechnology, The Royal Institute of Technology, S-100 44 Stockholm,1 Department of Bacteriology, National Veterinary Institute,2 and Department of Veterinary Microbiology, Swedish University of Agricultural Sciences,4 S-750 07 Uppsala, Sweden, and Département d'élevage et de médecine véterinaire, Centre de coopération internationale en recherche agronomique pour le dévelopment, F-34032 Montpellier, France3

Received 18 November 1997/Accepted 23 February 1998

    ABSTRACT
Top
Abstract
Introduction
Materials & Methods
Results & Discussion
References

Mycoplasma capricolum subsp. capripneumoniae belongs to the so-called Mycoplasma mycoides cluster and is the causal agent of contagious caprine pleuropneumonia (CCPP). All members of the M. mycoides cluster have two rRNA operons. The sequences of the 16S rRNA genes of both rRNA operons from 20 strains of M. capricolum subsp. capripneumoniae of different geographical origins in Africa and Asia were determined. Nucleotide differences which were present in only one of the two operons (polymorphisms) were detected in 24 positions. The polymorphisms were not randomly distributed in the 16S rRNA genes, and some of them were found in regions of low evolutionary variability. Interestingly, 11 polymorphisms were found in all the M. capricolum subsp. capripneumoniae strains, thus defining a putative ancestor. A sequence length difference between the 16S rRNA genes in a poly(A) region and 12 additional polymorphisms were found in only one or some of the strains. A phylogenetic tree was constructed by comparative analysis of the polymorphisms, and this tree revealed two distinct lines of descent. The nucleotide substitution rate of strains within line II was up to 50% higher than within line I. A tree was also constructed from individual operonal 16S rRNA sequences, and the sequences of the two operons were found to form two distinct clades. The topologies of both clades were strikingly similar, which supports the use of 16S rRNA sequence data from homologous operons for phylogenetic studies. The strain-specific polymorphism patterns of the 16S rRNA genes of M. capricolum subsp. capripneumoniae may be used as epidemiological markers for CCPP.

    INTRODUCTION
Top
Abstract
Introduction
Materials & Methods
Results & Discussion
References

rRNA sequences are, in general, believed to show low variability between and within species or subspecies. Heterogeneities in 16S rRNA genes have been reported but only to a minor extent. Nevertheless, both macroheterogeneities and microheterogeneities are known to exist. Macroheterogeneities involving large insertions ranging from 50 to several hundred nucleotides have been observed, e.g., in the archaeon Pyrobaculum aerophilum (6) and in the (eu)bacteria Campylobacter helveticus (29), Desulfotomaculum australicum (36), and a spore-forming Bacillus species (43). The first two are examples of species with split-gene formations of the rRNA genes, and the insertions are defined as intervening sequences. Therefore, although the intervening sequence is present in the structural gene and is even represented in the primary transcript, it is absent in the mature 16S rRNA molecule and does not contribute to the structure of the gene product. The last two species, however, have unusually long extensions of helices 6 and 49, respectively, according to Van de Peer et al. (54), which are positioned within the hypervariable regions V1 and V5, respectively, following the nomenclature of Gray et al. (11). These idiosyncrasies will cause problems only if they are not present in all of the rrn operons and when PCR based sequencing is used for determination of the nucleotide sequence. The resulting extra characters will be removed from the alignment which is used for inferring the tree and therefore do not constitute phylogenetic insignia. Microheterogeneities are probably by far more common than macroheterogeneities, and they are likely to be reported more frequently when we start looking for them. Clayton et al. recently observed that slightly different 16S rRNA sequences were deposited into the data banks for different strains belonging to the same species (8). Sequencing errors as the only plausible explanation for this variability were ruled out by the authors. Instead, most of the differences were believed to be real and to be caused by intraspecific variations. However, microheterogeneities in the form of nucleotide differences between the rrn operons, so-called polymorphisms, within a species and the extent to which they occur are not known. Examples of species where polymorphisms have been identified are Haloarcula marismortui (33), Bacillus sporothermodurans (42), and members of the class Mollicutes (16, 39-41, 44, 45).

About 175 species have been recognized within the class Mollicutes, and discoveries of new species are constantly being reported. The trivial name "mollicutes" will be used herein to avoid confusion with members of the genus Mycoplasma. The mollicutes have a small genome with a low G+C content and lack a cell wall, but they are phylogenetically related to gram-positive bacteria with a low G+C content in their genomes. Phylogenetic analysis of the mollicutes based on 16S rRNA sequences (55) has, together with other data, resulted in a revised taxonomy of this class, which is now composed of eight genera (53). The tree of the mollicutes revealed five distinct groups, of which one was named the spiroplasma group (55). A cluster of mycoplasmas within the spiroplasma group, which is of particular importance in veterinary medicine, is the so-called Mycoplasma mycoides cluster. All the members of the M. mycoides cluster are closely related, and some of them are difficult to differentiate by conventional techniques. Analysis of rRNA sequences also showed that M. putrefaciens is related to the members of the M. mycoides cluster (55). The following six mollicutes (9) denoted as species, subspecies or strains are included in the classical M. mycoides cluster: Mycoplasma capricolum subsp. capripneumoniae, Mycoplasma capricolum subsp. capricolum, Mycoplasma mycoides subsp. capri, Mycoplasma mycoides subsp. mycoides type LC, Mycoplasma mycoides subsp. mycoides type SC, and Mycoplasma sp. bovine serogroup 7. The M. mycoides cluster can be subdivided into the M. capricolum species group and the M. capri species group (41). M. capricolum subsp. capripneumoniae, formerly Mycoplasma sp. strain F38 (26), which belongs to the M. capricolum species group, causes contagious caprine pleuropneumonia (CCPP). CCPP is a goat disease of great concern in Africa and Asia (26, 30) and is included in the B list of communicable animal diseases of the Office International des Epizooties (22). CCPP was first described at the end of the last century (20, 51) and was shown to be caused by M. capricolum subsp. capripneumoniae in 1976 by MacOwan and Minette (31). More than 30 countries have declared that they have detected CCPP, but the organism has been isolated from goats in only 11 countries (48). A diagnostic method for CCPP based on PCR of the 16S rRNA genes from M. capricolum subsp. capripneumoniae and restriction enzyme analysis of the PCR product has been developed (45). The members of the M. mycoides cluster have two rRNA operons, designated rrnA and rrnB (7), and the above diagnostic method for CCPP was based on a polymorphism (nucleotide difference between the two 16S rRNA genes in the same strain) which was found to be unique for M. capricolum subsp. capripneumoniae. This polymorphism can therefore be easily used as a molecular marker for this subspecies, because it is localized in a restriction site for PstI (45). This unique polymorphism has been shown to be present in 16 strains of M. capricolum subsp. capripneumoniae from different parts of the world, but it was not present in 39 strains representing other species or subspecies of the M. mycoides cluster (4, 45).

The sequences of the 16S rRNA genes from both rRNA operons were recently determined for all members of the M. mycoides cluster, and several polymorphisms (or microheterogeneities) were identified (41). The strains F38T and 4/2LC of M. capricolum subsp. capripneumoniae were found to be unique among the members of the M. mycoides cluster, because they had the largest number of polymorphisms (15 and 17, respectively). M. mycoides subsp. mycoides type SC, the causal agent of contagious bovine pleuropneumonia, was found to have eight polymorphisms and a sequence length difference of two adenosines, whereas the other members of the M. mycoides cluster had only one, two, or three polymorphisms and sequence length variations were not found (41). Other mollicutes from ruminants have also been shown to have only few polymorphisms (44). Furthermore, the polymorphism patterns in the 16S rRNA genes were not identical for the two strains F38T and 4/2LC, even though they represent the same subspecies. Thus, M. capricolum subsp. capripneumoniae is unique among the mollicutes, which so far have been analyzed for the presence of polymorphisms in the 16S rRNA genes, which also supports the current classification of M. capricolum subsp. capripneumoniae into at least a separate subspecies (26). In the present study, we therefore determined the sequences of both 16S rRNA genes from 20 M. capricolum subsp. capripneumoniae strains from different geographical origins to investigate if polymorphisms can be used as epidemiological markers. The variations in the polymorphism patterns were found to be surprisingly great within M. capricolum subsp. capripneumoniae and could therefore be used to study molecular evolution within this subspecies.

    MATERIALS AND METHODS
Top
Abstract
Introduction
Materials & Methods
Results & Discussion
References

M. capricolum subsp. capripneumoniae strains, growth conditions, and sample preparations. The M. capricolum subsp. capripneumoniae strains used in this work are listed in Table 1, and the geographical distribution of the strains is shown in Fig. 1. All strains were grown in Hayflick's medium supplemented with pyruvate (50). Samples for PCR and DNA sequencing were prepared as described previously (41).

                              
View this table:
[in this window]
[in a new window]
  		TABLE 1.   M. capricolum subsp. capripneumoniae strains from which the 16S rRNA genes have been sequenced earlier or in this worka


View larger version (31K):
[in this window]
[in a new window]
  		FIG. 1.   Geographical origins of the M. capricolum subsp. capripneumoniae strains which have been analyzed. The affected countries are shown with bold boundaries. When the site of isolation is known, it is indicated with a symbol (Table 1).

In vitro amplification and DNA sequencing of the 16S rRNA genes. The 16S rRNA genes from both rRNA operons were amplified by seminested PCR. The reverse primers were biotinylated for magnetic separation of the strands. One outer primer pair was complementary to the universal regions U1 and U8 (11) and was used to amplify the 16S rRNA genes of both rRNA operons (21, 41). The other outer primer pair was complementary to the flanking regions of the 16S rRNA gene of the rrnB operon and was used for specific amplification of the gene from this operon (41). Thereafter, the respective product was used in a seminested amplification as described previously (21, 41). The biotinylated amplicons from the 16S rRNA genes were immobilized onto streptavidin-coated superparamagnetic beads (Dynabeads M280; Dynal AS, Oslo, Norway), and single strands suitable for bidirectional DNA sequencing were obtained by magnetic separation (18, 40). The sequencing reactions, performed by the method of Sanger et al. (47) with bacteriophage T7 DNA polymerase, were carried out automatically (17, 38, 40, 41). Detailed protocols and descriptions of primers for solid-phase 16S rDNA sequencing have been published (21, 38, 40, 41, 44).

Evaluation of sequence data. The sequence of the 16S rRNA gene of the rrnA operon was deduced from the sequences obtained by PCR with general primers and with rrnB-specific primers (41). A secondary-structure model of the 16S rRNA molecule transcribed from the rrnB operon of M. capricolum subsp. capripneumoniae was constructed by modification of the Postscript file of the 16S rRNA molecule of M. capricolum subsp. capricolum (12) retrieved from the Ribosomal Database Project (32). The polymorphisms of the different M. capricolum subsp. capripneumoniae strains were introduced into the model as described previously (41).

Two kinds of phylogenetic trees were inferred from the rRNA sequences. One tree was based on mutations and was inferred from the consensus sequences by comparative analysis of the sequence differences compiled in Table 2. By this procedure, all of the polymorphisms in the 16S rRNA genes which were found to be common to all strains constituted the microheterogeneity pattern identical to that of the ancestor. Thereafter, the exceptionally found polymorphisms were used to extrapolate the evolutionary lineages originating from this hypothetical ancestor. One evolutionary event (i.e., one polymorphism or sequence truncation) was given a value of 1. The second tree was derived from a corrected distance matrix. The one-parameter model of Jukes and Cantor (24) was used to correct the matrix at single locations, assuming equal frequencies and identical substitution rates of all nucleotides of the individual 16S rRNA sequences. The dendrogram was computed by using the tree-building method of Saitou and Nei (46).

                              
View this table:
[in this window]
[in a new window]
  		TABLE 2.   Mutational events in the 16S rRNA genes of different strains of M. capricolum subsp. capripneumoniaea

Nucleotide sequence accession numbers. The sequences of the 16S rRNA genes from the rrnA and the rrnB operons of the M. capricolum subsp. capripneumoniae strains have been deposited in GenBank (National Center for Biotechnology Information, Bethesda, Md.) under the accession numbers listed in Table 1.

    RESULTS AND DISCUSSION
Top
Abstract
Introduction
Materials & Methods
Results & Discussion
References

Nucleotide sequence heterogeneity in 16S rRNA genes. Despite the reports of variability within 16S rRNA genes (see Introduction), the extent to which microheterogeneities occur has not been extensively studied. This is partially due to the limitations in the detection techniques commonly used for rDNA sequencing. In the present study, we found that the 16S rRNA genes can differ by >1% between the operons within an individual strain of M. capricolum subsp. capripneumoniae. We have shown for M. capricolum subsp. capripneumoniae that the interoperon variation in, for instance, strains 4/2LC and GL102 is 1.1% as calculated from data in Table 2. We have also observed a strain-to-strain variation for a specific operon of up to 0.26% between several of the strains (Table 2).

The importance of determining the intraspecific variability in 16S rRNA genes for developing reliable phylogenetic hypothesis has recently been pointed out (8, 41). The usefulness of characterizing microheterogeneities is also demonstrated in this work, emphasizing the importance of sequencing the 16S rRNA genes of all operons (if more than one is present) and investigating the strain-to-strain variations in the 16S rRNA genes. It is, of course, also extremely important to deposit sequences in GenBank under the correct species name and strain designation. Consequently, a bacterial DNA or protein sequence should preferably not be accepted in a data bank without a recognized strain designation.

Nucleotide sequences of the 16S rRNA genes of M. capricolum subsp. capripneumoniae. The 16S rRNA gene sequences of both operons from 20 strains of the species M. capricolum subsp. capripneumoniae listed in Table 1 were subjected to bidirectional solid-phase rDNA sequencing, which has previously been shown to give very accurate data with the possibility of detecting heterogeneity, e.g., in viral populations and in 16S rRNA genes (see, e.g., references 16, 21, 28, 38, 40-42, and 44). Ten different 16S rRNA polymorphism patterns were found among the 20 M. capricolum subsp. capripneumoniae strains. Furthermore, a sequence length variation of one adenosine in a poly(A) region was found between the 16S rRNA genes of the two rRNA operons in nine of the strains. The poly(A) region is situated between nucleotides 444 and 449 and is not homologous to the poly(A) region with sequence length variations found in M. mycoides subsp. mycoides SC (between positions 1264 and 1270) reported previously (41). Length variations between the 16S rRNA genes of the members of the M. mycoides cluster have so far been observed only within poly(A) regions. These sequence length variations are probably caused by a process known as replication slippage (57).

A primary isolate of a mycoplasma may represent several clonal variants of a particular strain. A sample was therefore also analyzed for clonal variants with respect to 16S rRNA gene sequences among the M. capricolum subsp. capripneumoniae isolates. The 16S rRNA sequences of both operons obtained from a sample from Uganda were compared with those of the cloned strain M79/93 originating from this sample. No sequence differences were observed, which indicates that in a particular sample, only one 16S rRNA clone variant is normally present. The fact that the two nucleotide alternatives in a polymorphism were always present in a 1:1 ratio also supports this observation.

The hypothetical ancestor of the M. capricolum subsp. capripneumoniae strains and two lines of descent. The 24 evolutionary events were compiled and compared as detailed in Table 2. Eleven polymorphic positions were found to be common to all 20 strains, thus defining a hypothetical ancestor from which the M. capricolum subsp. capripneumoniae strains have evolved. The sequences of the two 16S rRNA genes of the hypothetical ancestor can be derived from Table 2. The remaining 12 polymorphisms and the length variation which occurred in only some of the strains most probably arose due to later evolutionary events. The polymorphisms in the 16S rRNA genes of the M. capricolum subsp. capripneumoniae strains (Table 2) were converted into a cladogram (Fig. 2) by comparative analysis. Both the length polymorphism and a nucleotide polymorphism were treated as one event. The horizontal lines are proportional to the number of events. The tree shown in Fig. 2 was rooted by the 11 ancestral evolutionary events and revealed two major lines of descent.


View larger version (7K):
[in this window]
[in a new window]
  		FIG. 2.   Phylogenetic tree based on mutational events of the 16S rRNA genes of the M. capricolum subsp. capripneumoniae strains. The tree was constructed by comparative analysis and with the 11 polymorphisms of the putative ancestor used as the root. Two major lines of descent (I and II) can be seen. The positions and types of the polymorphisms are shown on the axis.

Evolutionary line I contained the 9 M. capricolum subsp. capripneumoniae strains which shared the synapomorphy (444-449)A/- (Table 2; Fig. 2), situated in a poly(A) segment of the molecule in a region which is characterized by bilaterally bulged residues (Fig. 3). It remains to be shown whether the presence or absence of a sixth adenosine in this region of the 16S rRNA molecule will affect the function of the mature small ribosomal subunit from the strains of line I. The strains of line I showed a very low interstrain variability, and only two additional polymorphic positions, 509T/C and 875G/A were found, thus bifurcating into sublines IA and IB (Table 2; Fig. 2). Strain 95043 from Niger differed from the hypothetical ancestor only by the truncation of one adenosine in the 16S rRNA gene of the rrnB operon.


View larger version (31K):
[in this window]
[in a new window]
  		FIG. 3.   Secondary-structure model of the 16S rRNA molecule from the rrnB operon from M. capricolum subsp. capripneumoniae. All identified polymorphisms are indicated with arrows. They are numbered, and their composition (Table 2) is given according to the letter code of the International Union of Biochemistry. The nucleotide residues of the rrnA and rrnB operons are written at the arrows as N/N, respectively. The segment containing the sequence length variation between rrnA and rrnB is denoted by an arc. This model has been adapted from the secondary-structure model of the 16S rRNA molecule of M. capricolum subsp. capricolum described by Gutell et al. (13, 14).

The 11 M. capricolum subsp. capripneumoniae strains of evolutionary line II shared the two polymorphisms 404C/A and 1080C/A (Table 2; Fig. 2). The earliest member of line II was M. capricolum subsp. capripneumoniae 9081 from Oman, and line II branched into the two sublines IIA and IIB with a fork of three additional branches at the node of IIB (Fig. 2). This line is by far the most heterogeneous, and the evolutionary drift seems to be more pronounced for the 16S rRNA genes in the strains of line II than in those of line I. The two lines are separated by only three mutational events but can nevertheless be recognized as true phylogenetic entities, since none of the polymorphisms present in only some of the strains were shared among M. capricolum subsp. capripneumoniae strains of the two lines (Table 2). This observation can therefore be regarded as a consistency measure for the two lines, and it also indicates that sequence homogenization due to gene conversion has not occurred as was observed for the tuf genes of Salmonella typhimurium (1).

The 16S rRNA gene sequences of the five M. capricolum subsp. capripneumoniae strains, which were used in a recent study on the DNA relatedness of strains from the M. capricolum species group (5), have also been determined in this work. Although the hybridization values were very similar for all these strains (87 to 89%), the 16S rRNA sequence data clearly showed that two of the strains (F38T and 7/1a) belonged to line II and three of the strains (G1943/80, G280/80, and G94/83) belonged to line I (Fig. 2; Table 2).

Polymorphisms as epidemiological markers and diagnostic targets. All members of the M. mycoides cluster are closely related as judged by biochemistry (10, 41), serological reactions (25, 35), DNA-DNA reassociation studies (2), and 16S rRNA sequence analysis (41). It is therefore important to identify polymorphisms in the 16S rRNA genes, because they can be used in diagnostic applications to distinguish between closely related species (45). The polymorphisms may also be useful as epidemiological markers. Two major lines of descent have been defined for M. capricolum subsp. capripneumoniae strains in this work. Representatives of line I were found only in Africa, whereas representatives of line II were found in both Africa and Asia (Fig. 1). A representative of line I, strain 95043, was isolated in Niger in 1995, and a strain of subline IB (Fig. 2; Table 2) was isolated in the neighboring country of Chad. Several strains of subline IA were isolated in Kenya and Uganda. The occurrence of the same type in neighboring countries, such as IA in Kenya and Uganda and IIB in Kenya and Sudan, might well be explained by trade and other movement of animals across these borders. An area of eastern Uganda and western Kenya is, for instance, populated by the same tribe, and the border is not well respected.

Strains of line II were found in countries bordering maritime routes (Mediterranean sea: Turkey, Tunisia; Indian Ocean: Kenya, Oman; Red Sea: Ethiopia, Sudan). Strain 9081, with a mutation pattern of line II, was isolated in Oman in 1990 and differed from the hypothetical ancestor by two polymorphisms in positions 404 and 1080. The polymorphism type pattern of subline IIA is the closest to that of line II and differs by one extra polymorphism in position 232. Strains of this type were isolated in Dubai in 1991. The polymorphism type pattern of the subline IIB differed from that of line II by the two polymorphisms in positions 94 and 1403. Strains of this type have been isolated in Kenya in 1976 (strain F38) and Sudan in 1989. Strains of types IIB1 and IIB2 differed from the strains of type IIB by one polymorphism in positions 538 and 684, respectively. A strain of type IIB1 was isolated in Ethiopia in 1992, and representatives of type IIB2 were isolated in Tunisia 1980 and Oman 1988. The Oman strain 7/1a originated from a goat which was imported from Turkey. Strains with a polymorphism pattern of type IIB3 differed from type IIB by two polymorphisms in positions 1255 and 1297. Representatives of type IIB3 were isolated in Oman.

Three types of M. capricolum subsp. capripneumoniae were found in Oman, namely, II, IIB2, and IIB3. A plausible explanation for this would be the extensive importation of goats from several countries in the Arabic peninsula for the religious feasts. These animals are normally slaughtered, but some might be kept alive and can thus infect indigenous herds. It is also noteworthy that two countries which were most likely to have CCPP-infected animals in the 19th century, Algeria (51) and South Africa (20), might have had strains of type IIB2. Algeria has a border with Tunisia, where type IIB2 was found in 1980, and South Africa imported Angora goats, which most probably brought CCPP from Turkey in 1881. The Turkish goat imported to Oman in 1988 was also found to be of type IIB2. In general, no correlation could be seen between the polymorphism pattern and the time of isolation of the strains.

Secondary structure analysis of the polymorphic positions and their implications. The positions of the 24 mutational events in the 16S rRNA molecule of M. capricolum subsp. capripneumoniae are shown in the secondary-structure model in Fig. 3. The nucleotide substitution rates according to a recently published variability map of the 16S rRNA molecule (54) are included in Table 2, showing the variability in each of the polymorphic positions. For comparison, consensus 16S rRNA sequences from an alignment of the mollicutes (32, 39) were computed by using different cutoff values for which a residue in a certain position was present. These percentages are included in Table 2, indicating that positions of high, intermediate, and low variability found in (eu)bacteria (54) also seem to hold for the variability in the 16S rRNA genes of mycoplasmas. Of the 11 polymorphisms of the ancestor, 6 (positions 180, 452, 672, 844, 1151, and 1446) were located in the high-variability groups 1 and 2 as defined by Van de Peer et al. (54). In contrast, 6 of the 12 microheterogeneities found in only some of the M. capricolum subsp. capripneumoniae strains (i.e., positions 94, 404, 509, 538, 684, and 897) belonged to group 5 or 6 of low variability. Moreover, frequently observed polymorphisms situated in stems were found predominantly in positions of relatively high variability (positions 452, 672, 1080, 1151, 1403, and 1446), while rarely found microheterogeneities often were situated in highly conserved locales in the 16S rRNA molecule (positions 538, 871, 875, 897, and 1255). It is noteworthy that all but 875 of the rarely found polymorphisms in sites of low variability occur in the starting or ending base pair of a stem (positions 538, 871, 897, and 1255). The reason for this is not known, but the base pairing at these positions might be rather flexible, at least in these mycoplasmas. This observation indicates that polymorphisms shared among all strains occur in rather variable positions while rarely found polymorphisms are present in more highly conserved positions of the 16S rRNA molecule of M. capricolum subsp. capripneumoniae.

A total of 13 polymorphisms are situated in stems (Fig. 3). A compensatory mutation in the opposite position to stabilize the actual stem in the other molecule was never observed. However, the general alternatives to noncanonical base pairing (13, 54) were followed in most cases. Therefore, a guanosine, a uridine, and an adenosine residue in one strand can have C or U, A or G, and U or C, respectively, as their counterparts. Moreover, the polymorphisms 509T/C, 875G/A, and 1255A/G indicated that cytidines can tolerate A or G as pairing nucleotides and that a C · G base pair can be substituted with a C · A pair, in certain positions. Therefore, a C · A pair might not necessarily affect the stability of a stem (54), and plausible explanations are a protonated C or a tautomeric configuration of A or C in these positions, which have been suggested in previous reports (19, 54). Interestingly, over 40 polymorphisms have been localized to stem regions of the 16S rRNA molecule of mycoplasmas (39, 41, 44, also see above) and without a compensatory mutation in the corresponding nucleotide position. This finding is contradictory to the polymorphisms in the 16S rRNA genes of the recently described B. sporothermodurans (42), where most of the microheterogeneities found in stems were associated with compensatory mutations in the complementary nucleotide position.

All polymorphisms of the M. capricolum subsp. capripneumoniae strains were checked against the corresponding positions of the 16S rRNA mutation database which provides a list of mutated positions in Escherichia coli (52). Only position 897 in M. capricolum subsp. capripneumoniae (912 in E. coli)T/C was listed in the database. A change of a C to a U in this position in E. coli has been shown to confer resistance to streptomycin (34). The actual polymorphism 897(912)T/C is present in strain GL102, a laboratory strain originating from strain Gabés, which has been passaged 102 times in the laboratory. It has not been possible to find whether streptomycin has been used during any of the passages. Nevertheless, this transition indicated that mutations may also occur in positions of low variability (Table 2) with a significant rate in the 16S rRNA genes.

We have so far characterized over 80 microheterogeneities in the 16S rRNA genes of different Mycoplasma species isolated from different hosts (16, 39, 41, 44). A comparison of these polymorphisms with those determined in this study confirmed that microheterogeneities in the 16S rRNA genes of mycoplasmas are distributed throughout the molecule but are rarely seen in universal regions (11) of the molecule. Polymorphisms with identical nucleotide compositions and present in different species have been found in only 3 of the more than 80 positions.

Evolution of the 16S rRNA genes of M. capricolum subsp. capripneumoniae: rrnA operon versus rrnB operon. The 16S rRNA sequences of the rrnA operon were 99.73 to 99.93% similar for the different M. capricolum subsp. capripneumoniae strains. Identical values were found for the rrnB operon. The percent similarity ranged from 98.77 to 99.20% between the two operons. This means that the similarity between 16S rRNA genes within a strain is lower (11 to 17 nucleotide differences) than that between 16S rRNA genes of homologous operons in different strains (0 to 4 nucleotide differences). A phylogenetic tree was constructed from 16S rRNA sequences of the individual operons (Fig. 4). Mycoplasma putrefaciens served as the outgroup, and the branches are denoted according to the corresponding polymorphism pattern (Table 1). The phylogenetic analysis revealed two distinct entities, one of which consisted of sequences from the rrnA operon and the other consisted of sequences from the rrnB operon. A slightly higher evolutionary rate was observed for the sequences of the rrnB operon as judged from the branch lengths (Fig. 4). The tree implies that the rrnB operon has evolved by one to three further nucleotide changes. Strikingly, the topologies for the respective clade are very similar, and the line I strains form deep branches of both clades. Descendants belonging to line II form more recent branches, where M. capricolum subsp. capripneumoniae 9081 of subline II is placed intermediately in the respective operonal cluster. This shows that the 16S rRNA genes from the two operons coevolved and that both are phylogenetically informative. The tree illustrates the importance of using 16S rRNA sequences from homologous operons to infer the phylogeny when polymorphisms are present and the species to be analyzed are closely related.


View larger version (13K):
[in this window]
[in a new window]
  		FIG. 4.   Phylogenetic tree based on 16S rRNA sequences of the individual operons obtained by the neighbor-joining method (46). The tree revealed two distinct clades representing sequences of the rrnA and rrnB operons, indicating that the individual operons are monophyletic and that gene conversion has not taken place.

Conclusion. The evolution within mycoplasmas has been reported to be unusually rapid (56). In a recent study, we suggested that members of the M. mycoides cluster could be used as a model system for molecular evolution by mapping the polymorphisms in the two 16S rRNA genes (41). In the present work, we have characterized microheterogeneities in strains of M. capricolum subsp. capripneumoniae, a subspecies of the M. capricolum species group (41) of the M. mycoides cluster (55), and have constructed two kinds of trees. One tree was based on mutational events from consensus sequences, and this tree formed two lines of descent without underlying synapomorphies (close parallelism as a result of common inherited genetic factors causing incomplete synapomorphy). In a tree based on individual operons, both 16S rRNA genes were found to evolve and to reflect similar phylogeny. Therefore, strains of M. capricolum subsp. capripneumoniae constitute a very useful model for studies of molecular evolution. M. capricolum subsp. capripneumoniae may also be a suitable subspecies for comparing the evolution of other genes with the evolution observed for the 16S rRNA genes.

    ACKNOWLEDGMENTS

We are grateful to Diarmaid Hughes for valuable discussions on gene conversion. We also thank Marianne Persson for skillfully performing all PCR experiments and Joseph Tully, Henning Ernø, Gareth Jones, and Hezron Wesonga for generous supply of strains.

This work has been financially supported by grants from the Göran Gustafsson Foundation and the Swedish Engineering Science Research Council to M.U. and from the Commission of the EU (DG XII) for "Contagious caprine pleuropneumonia: distribution, evaluation of vaccines, and development of new tools" (contract IC18-CT95-0007), from the Swedish Council for Forestry and Agricultural Research, and from the Swedish International Development Cooperation Agency to K.-E.J. This study also forms a part of the EU research collaboration COST 826 on Ruminants' mycoplasmoses.

    FOOTNOTES

* Corresponding author. Mailing address: Department of Bacteriology, National Veterinary Institute, P.O. Box 7073, S-750 07 Uppsala, Sweden. Phone: 46 18 67 40 00. Fax: 46 18 30 91 62. E-mail: Kaggen{at}sva.se.

    REFERENCES
Top
Abstract
Introduction
Materials & Methods
Results & Discussion
References

1. Abdulkarim, F., and D. Hughes. 1996. Homologous recombination between the tuf genes of Salmonella typhimurium. J. Mol. Biol. 260:506-522[Medline].
2. Askaa, G., H. Ernø, and M. O. Ojo. 1978. Bovine mycoplasmas: classification of groups related to Mycoplasma mycoides. Acta Vet. Scand. 19:166-178[Medline].
3. Bölske, G., K.-E. Johansson, R. Heinonen, P. A. Panvuga, and E. Twinamasiko. 1995. Contagious caprine pleuropneumonia in Uganda and isolation of Mycoplasma capricolum subspecies capripneumoniae from goats and sheep. Vet. Rec. 137:594[Medline].
4. Bölske, G., J. G. Mattsson, C. Ros Bascuñana, K. Bergström, H. Wesonga, and K.-E. Johansson. 1996. Diagnosis of contagious caprine pleuropneumonia by detection and identification of Mycoplasma capricolum subsp. capripneumoniae by PCR and restriction enzyme analysis. J. Clin. Microbiol. 34:785-791[Abstract].
5. Bonnet, F., C. Saillard, J. M. Bové, R. H. Leach, D. L. Rose, G. S. Cottew, and J. G. Tully. 1993. DNA relatedness between field isolates of mycoplasma F38 group, the agent of contagious caprine pleuropneumonia and strains of Mycoplasma capricolum. Int. J. Syst. Bacteriol. 43:597-602[Abstract/Free Full Text].
6. Burggraf, S., N. Larsen, C. R. Woese, and K. O. Stetter. 1993. An intron within the 16S ribosomal RNA gene of the archeon Pyrobaculum aerophilum. Proc. Natl. Acad. Sci. USA 90:2547-2550[Abstract/Free Full Text].
7. Christiansen, G., and H. Ernø. 1990. RFLP in rRNA genes of Mycoplasma capricolum, the caprine F38-like group and the bovine serogroup 7. Zentralbl. Bakteriol. Suppl. 20:479-488.
8. Clayton, R. A., G. Sutton, P. S. Hinkle, Jr., C. Bult, and C. Fields. 1995. Intraspecific variation in small-subunit rRNA sequences in GenBank: why single sequences may not adequately represent prokaryotic taxa. Int. J. Syst. Bacteriol. 45:595-599[Abstract/Free Full Text].
9. Cottew, G. S., A. Breard, A. J. DaMassa, H. Ernø, R. H. Leach, P. C. Lefèvre, A. W. Rodwell, and G. R. Smith. 1987. Taxonomy of the Mycoplasma mycoides cluster. Isr. J. Med. Sci. 23:632-635[Medline].
10. Ernø, H., and L. Stipkovits. 1973. Bovine mycoplasmas: cultural and biochemical studies. II. Acta Vet. Scand. 14:450-463[Medline].
11. Gray, M. W., D. Sankoff, and R. J. Cedergren. 1984. On the evolutionary descent of organisms and organelles: a global phylogeny based on a highly conserved structural core in small subunit ribosomal RNA. Nucleic Acids Res. 12:5837-5852[Abstract/Free Full Text].
12. Gutell, R. R. 1994. Collection of small subunit (16S- and 16S-like) ribosomal RNA structures: 1994. Nucleic Acids Res. 22:3502-3507[Abstract/Free Full Text].
13. Gutell, R. R., N. Larsen, and C. R. Woese. 1994. Lessons from an evolving rRNA: 16S and 23S rRNA structures from a comparative perspective. Microbiol. Rev. 58:10-26[Abstract/Free Full Text].
14. Gutell, R. R., B. Weiser, C. R. Woese, and H. F. Noller. 1985. Comparative anatomy of 16-S-like ribosomal RNA. Prog. Nucleic Acid Res. Mol. Biol. 32:155-216[Medline].
15. Harbi, M. S. M. A., M. S. El Tahir, K. J. MacOwan, and A. A. Nayil. 1981. Mycoplasma strain F38 and contagious caprine pleuropneumonia in Sudan. Vet. Rec. 108:261[Medline].
16. Heldtander, M., B. Pettersson, J. G. Tully, and K.-E. Johansson. 1998. Sequences of the 16S rRNA genes and phylogeny of the goat mycoplasmas Mycoplasma adleri, Mycoplasma auris, Mycoplasma cottewii, and Mycoplasma yeatsii. Int. J. Syst. Bacteriol. 48:263-268[Abstract/Free Full Text].
17. Hultman, T., S. Bergh, T. Moks, and M. Uhlén. 1991. Bidirectional solid-phase sequencing of in vitro-amplified plasmid DNA. BioTechniques 10:84-93. [Medline]
18. Hultman, T., S. Ståhl, E. Hornes, and M. Uhlén. 1989. Direct solid phase sequencing of genomic and plasmid DNA using magnetic beads as solid support. Nucleic Acids Res. 17:4937-4946[Abstract/Free Full Text].
19. Hunter, W. N., T. Brown, N. N. Anand, and O. Kennard. 1986. Structure of an adenine · cytosine base pair in DNA and its implications for mismatch repair. Nature (London) 320:552-555[Medline].
20. Hutcheon, D. 1881. Contagious pleuropneumonia in Angora goats. Vet. J. 13:171-180.
21. Johansson, K.-E., M. U. K. Heldtander, and B. Pettersson. 1998. Characterization of mycoplasmas by PCR and sequence analysis with universal 16S rDNA primers, p. 145-165. In R. J. Miles, and R. A. J. Nicholas (ed.), Methods in molecular medicine, vol. 104. Mycoplasma protocols. Humana Press Inc., Totowa, N.J.
22. Jones, G. E. 1992. Contagious caprine pleuropneumonia, p. 376-392. In M. Truszczynsky, J. Pearson, and S. Edwards (ed.), The OIE manual of standards for diagnostic tests and vaccines. Office International des Epizooties, Paris, France.
23. Jones, G. E., and A. R. Wood. 1988. Microbiological and serological studies on caprine pneumonias in Oman. Res. Vet. Sci. 44:125-131[Medline].
24. Jukes, T. H., and C. R. Cantor. 1969. Evolution of protein molecules, p. 21-132. In H. H. Munro (ed.), Mammalian protein metabolism. Academic Press, Inc., New York, N.Y.
25. Kanyi Kibe, M., and G. R. Smith. 1984. A study of F38-type and related mycoplasmas by mycoplasmaemia and cross-immunization tests in mice. J. Hyg. Camb. 93:465-473.
26. Leach, R. H., H. Ernø, and K. J. MacOwan. 1993. Proposal for designation of F38-type caprine mycoplasmas as Mycoplasma capricolum subsp. capripneumoniae subsp. nov. and consequent obligatory relegation of strains currently classified as M. capricolum (Tully, Barile, Edward, Theodore, and Ernø 1974) to an additional new subspecies, M. capricolum subsp. capricolum subsp. nov. Int. J. Syst. Bacteriol. 43:603-605[Abstract/Free Full Text].
27. Lefèvre, P. C., A. Breard, I. Al Faroukh, and S. Buron. 1987. Mycoplasma species F38 isolated in Chad. Vet. Rec. 121:575-576[Medline].
28. Leitner, T., E. Halapi, G. Scarlatti, P. Rossi, J. Albert, E. M. Fenyö, and M. Uhlén. 1993. Analysis of heterogeneous viral populations by direct DNA sequencing. BioTechniques 15:120-126. [Medline]
29. Linton, D., F. E. Dewhirst, J. P. Clewley, R. J. Owen, A. P. Burnens, and J. Stanley. 1994. Two types of 16S rRNA gene are found in Campylobacter helveticus: analysis, applications and characterization of the intervening sequence found in some strains. Microbiology 140:847-855[Abstract/Free Full Text].
30. MacOwan, K. J. 1984. Role of mycoplasma strain F38 in contagious caprine pleuropneumonia. Isr. J. Med. Sci. 20:979-981[Medline].
31. MacOwan, K. J., and J. E. Minette. 1976. A mycoplasma from acute contagious caprine pleuropneumonia in Kenya. Trop. Anim. Health Prod. 8:91-95[Medline].
32. Maidak, B. L., G. J. Olsen, N. Larsen, R. Overbeek, M. J. McCaughey, and C. R. Woese. 1996. The ribosomal database project (RDP). Nucleic Acids Res. 24:82-85[Abstract/Free Full Text].
33. Mevarech, M., S. Hirsch-Twizer, S. Goldman, E. Yakobson, H. Eisenberg, and P. P. Dennis. 1989. Isolation and characterization of the rRNA gene clusters of Halobacterium marismortui. J. Bacteriol. 171:3479-3485[Abstract/Free Full Text].
34. Montandon, P. E., R. Wagner, and E. Stutz. 1986. E. coli ribosomes with a C912 to U base change in the 16S rRNA are streptomycin resistant. EMBO J. 5:3705-3708[Medline].
35. Olsson, B., G. Bölske, K. Bergström, and K.-E. Johansson. 1990. Analysis of caprine mycoplasmas and mycoplasma infections in goats using two-dimensional electrophoresis and immunoblotting. Electrophoresis 11:861-869[Medline].
36. Patel, B. K. C., C. A. Love, and E. Stackebrandt. 1992. Helix 6 of the 16S rRNA of the bacterium Desulfotomaculum australicum exhibits an unusual structural idiosyncrasy. Nucleic Acids Res. 20:5483[Free Full Text].
37. Perreau, P., A. Breard, and C. Le Goff. 1984. Infection Expérimentale de la chévre par les souches de mycoplasme de type F.38 (pleuropneumonie contagieuse caprine). Ann. Microbiol. (Paris) 135A:119-124.
38. Pettersson, B. 1997. In Direct solid-phase 16S rDNA sequencing: a tool in bacterial phylogeny. Ph.D. thesis. Royal Institute of Technology, Stockholm, Sweden.
39. Pettersson, B., and K.-E. Johansson. Unpublished data.
40. Pettersson, B., K.-E. Johansson, and M. Uhlén. 1994. Sequence analysis of 16S rRNA from mycoplasmas by direct solid-phase DNA sequencing. Appl. Environ. Microbiol. 60:2456-2461[Abstract/Free Full Text].
41. Pettersson, B., T. Leitner, M. Ronaghi, G. Bölske, M. Uhlén, and K.-E. Johansson. 1996. Phylogeny of the Mycoplasma mycoides cluster as determined by sequence analysis of the 16S rRNA genes from the two rRNA operons. J. Bacteriol. 178:4131-4142[Abstract/Free Full Text].
42. Pettersson, B., F. Lembke, P. Hammer, E. Stackebrandt, and F. G. Priest. 1996. Bacillus sporothermodurans, a new species producing highly heat-resistant endospores. Int. J. Syst. Bacteriol. 46:759-764[Abstract/Free Full Text].
43. Pettersson, B., and F. G. Priest. Unpublished data.
44. Pettersson, B., M. Uhlén, and K.-E. Johansson. 1996. Phylogeny of some mycoplasmas from ruminants based on 16S rRNA sequences and definition of a new cluster within the hominis group. Int. J. Syst. Bacteriol. 46:1093-1098[Abstract/Free Full Text].
45. Ros Bascuñana, C., J. G. Mattsson, G. Bölske, and K.-E. Johansson. 1994. Characterization of the 16S rRNA genes from Mycoplasma sp. strain F38 and development of an identification system based on PCR. J. Bacteriol. 176:2577-2586[Abstract/Free Full Text].
46. Saitou, N., and M. Nei. 1987. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 4:406-425[Abstract].
47. Sanger, F., S. Nicklen, and A. R. Coulson. 1977. DNA sequencing with chain-terminating inhibitors. Proc. Natl. Acad. Sci. USA 74:5463-5467[Abstract/Free Full Text].
48. Thiaucourt, F., and G. Bölske. 1996. Contagious caprine pleuropneumonia and other pulmonary mycoplasmoses of sheep and goats. Rev. Sci. Tech. OIE 15:1397-1414.
49. Thiaucourt, F., A. Breard, P. C. Lefèvre, and G. Y. Mebratu. 1992. Contagious caprine pleuropneumonia in Ethiopia. Vet. Rec. 131:585[Medline].
50. Thiaucourt, F., C. Guérin, V. Mady, and P. C. Lefèvre. 1992. Diagnostic de la pleuropneumonie contagieuse caprine: améliorations récentes. Rev. Sci. Tech. OIE 11:859-865.
51. Thomas, P. 1893. Rapport médical sur la maladie des chèvres désgnée par les Arabes sous le nom de Bou-Frida gui a sévi epizootiquement pendant l'hiver 1872-1873 dans le cercle de Djelfa, p. 1-35. In A. Jourdan (ed.), Rapport médical sur le Bou-Frida. Gouvernment général civil de L'Algérie, Algiers, Algeria.
52. Triman, K. L. 1996. The 16S ribosomal RNA mutation database (16SMDB). Nucleic Acids Res. 24:166-168[Abstract/Free Full Text].
53. Tully, J. G., J. M. Bové, F. Laigret, and R. F. Whitcomb. 1993. Revised taxonomy of the class Mollicutes: proposed elevation of a monophyletic cluster of arthropod-associated mollicutes to ordinal rank (Entomoplasmatales ord. nov.), with provision for familial rank to separate species with nonhelical morphology (Entomoplasmataceae fam. nov.) from helical species (Spiroplasmataceae), and emended descriptions of the order Mycoplasmatales, family Mycoplasmataceae. Int. J. Syst. Bacteriol. 43:378-385[Abstract/Free Full Text].
54. Van de Peer, Y., S. Chapelle, and R. De Wachter. 1996. A quantitative map of nucleotide substitution rates in bacterial rRNA. Nucleic Acids Res. 24:3381-3391[Abstract/Free Full Text].
55. Weisburg, W. G., J. G. Tully, D. L. Rose, J. P. Petzel, H. Oyaizu, D. Yang, L. Mandelco, J. Sechrest, T. G. Lawrence, J. Van Etten, J. Maniloff, and C. R. Woese. 1989. A phylogenetic analysis of the mycoplasmas: basis for their classification. J. Bacteriol. 171:6455-6467[Abstract/Free Full Text].
56. Woese, C. R., E. Stackebrandt, and W. Ludwig. 1985. What are mycoplasmas: the relationship of tempo and mode in bacterial evolution. J. Mol. Evol. 21:305-316.
57. Wolfson, R., K. G. Higgins, and B. B. Sears. 1991. Evidence for replication slippage in the evolution of Oenothera chloroplast DNA. Mol. Biol. Evol. 8:709-720[Abstract].

J Bacteriol, May 1998, p. 2350-2358, Vol. 180, No. 9
0021-9193/98/$04.00+0
Copyright © 1998, American Society for Microbiology. All rights reserved.


This article has been cited by other articles:

  • Morandi, A., Zhaxybayeva, O., Gogarten, J. P., Graf, J. (2005). Evolutionary and Diagnostic Implications of Intragenomic Heterogeneity in the 16S rRNA Gene in Aeromonas Strains. J. Bacteriol. 187: 6561-6564 [Abstract] [Full Text]  
  • Nicolas, M. M., Stalis, I. H., Clippinger, T. L., Busch, M., Nordhausen, R., Maalouf, G., Schrenzel, M. D. (2005). Systemic Disease in Vaal Rhebok (Pelea capreolus) Caused by Mycoplasmas in the Mycoides Cluster. J. Clin. Microbiol. 43: 1330-1340 [Abstract] [Full Text]  
  • Picard, F. J., Ke, D., Boudreau, D. K., Boissinot, M., Huletsky, A., Richard, D., Ouellette, M., Roy, P. H., Bergeron, M. G. (2004). Use of tuf Sequences for Genus-Specific PCR Detection and Phylogenetic Analysis of 28 Streptococcal Species. J. Clin. Microbiol. 42: 3686-3695 [Abstract] [Full Text]  
  • Coram, N. J., Rawlings, D. E. (2002). Molecular Relationship between Two Groups of the Genus Leptospirillum and the Finding that Leptospirillum ferriphilum sp. nov. Dominates South African Commercial Biooxidation Tanks That Operate at 40{degrees}C. Appl. Environ. Microbiol. 68: 838-845 [Abstract] [Full Text]  
  • Persson, A., Pettersson, B., Bölske, G., Johansson, K.-E. (1999). Diagnosis of Contagious Bovine Pleuropneumonia by PCR-Laser- Induced Fluorescence and PCR-Restriction Endonuclease Analysis Based on the 16S rRNA Genes of Mycoplasma mycoides subsp. mycoides SC. J. Clin. Microbiol. 37: 3815-3821 [Abstract] [Full Text]  

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrowReprints and Permissions
Right arrow Copyright Information
Right arrow Books from ASM Press
Right arrow MicrobeWorld
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Pettersson, B.
Right arrow Articles by Johansson, K.-E.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Pettersson, B.
Right arrow Articles by Johansson, K.-E.

Copyright © 2009 by the American Society for Microbiology.   For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»