The complete mitochondrial DNA (mtDNA) of the donkey and mtDNA comparisons among four closely related mammalian species-pairs

The nucleotide sequence of the complete mitochondrial genome of the donkey, Equus asinus, was determined. The length of the molecule is 16,670 bp. The length, however, is not absolute due to pronounced heteroplasmy caused by variable numbers of two types of repetitive motifs in the control region. The sequence of the repeats is (a) 5'-CACACCCA and (b) 5'-TGCGCGCA, respectively. The order of (a) and (b) can be expressed as {n[2(a)+(b)]+m(a)}. In 32 different clones analyzed the number of n and m ranged from 0 to 9 and 1 to 7. The two rRNA genes, the 13 peptide-coding genes, and the 22 tRNA genes of the donkey and the horse, Equus caballus, were compared in detail. Total nucleotide difference outside the control region was 6.9%. Nucleotide difference between peptide-coding genes ranged from 6.4% to 9.4% with a mean of 8.0%. In the inferred protein sequences of the 13 peptide-coding genes the amino acid difference was 0.2-8.8%, and the mean for the 13 concatenated amino acid sequences was 1.9%. In the 22 tRNA genes, the mean difference was 3.5%, and that in the two rRNA genes was 4.1%. The mtDNA differences between the donkey and the horse suggest that the evolutionary separation of the two species occurred approximately 9 million years ago. Analyses of differences among the mtDNAs of three other species-pairs, harbor seal/grey seal, fin whale/blue whale, and Homo/common chimpanzee, showed that the relative evolutionary rate of individual peptide-coding genes varies among different species-pairs and modes of comparison. The findings show that the superimposition of sequence data of one lineage for resolving and dating evolutionary divergences of other lineages should be performed with caution unless based on comprehensive data.     

Contrasting population structure from nuclear intron sequences and mtDNA of humpback whales

Powerful analyses of population structure require information from multiple genetic loci. To help develop a molecular toolbox for obtaining this information, we have designed universal oligonucleotide primers that span conserved intron-exon junctions in a wide variety of animal phyla. We test the utility of exon-primed, intron-crossing amplifications by analyzing the variability of actin intron sequences from humpback, blue, and bowhead whales and comparing the results with mitochondrial DNA (mtDNA) haplotype data. Humpback actin introns fall into two major clades that exist in different frequencies in different oceanic populations. It is surprising that Hawaii and California populations, which are very distinct in mtDNAs, are similar in actin intron alleles. This discrepancy between mtDNA and nuclear DNA results may be due either to differences in genetic drift in mitochondrial and nuclear genes or to preferential movement of males, which do not transmit mtDNA to offspring, between separate breeding grounds. Opposing mtDNA and nuclear DNA results can help clarify otherwise hidden patterns of structure in natural populations.     

Isolation of T cell receptor gamma/delta chain constant genes from minke whale

The polymerase chain reaction using oligonucleotides based on consensus sequences from other species allowed us to amplify minke whale (Balaenoptera acutorostrata) T cell receptor (TcR) gamma and delta constant genes. Two types of Cgamma clones were obtained which only differ by one mismatch. These minke whale Cgamma clones showed 83% nucleotide and 70% amino acid sequence similarity to the corresponding region of the human TRGC genes. The minke whale Cdelta sequences were all identical, and showed 86% nucleotide and 77.8% amino acid sequence similarity to the human TRDC gene. The whale Cgamma and Cdelta clones were used as probes for Southern blot analysis. The data confirmed that there is a unique Cdelta gene in minke whale. Two different Cgamma genes were detected, one of them also cross-hybridizing with a human Cgamma probe, suggesting that the Cgamma locus in minke whale consists of at least two different constant genes.     

Evolution of fragmented mitochondrial ribosomal RNA genes in Chlamydomonas

The fragmented mitochondrial ribosomal RNAs (rRNAs) of the green algae Chlamydomonas eugametos and Chlamydomonas reinhardtii are discontinuously encoded in subgenic modules that are scrambled in order and interspersed with protein coding and tRNA genes. The mitochondrial rRNA genes of these two algae differ, however, in both the distribution and organization of rRNA coding information within their respective genomes. The objectives of this study were (1) to examine the phylogenetic relationships between the mitochondrial rRNA gene sequences of C. eugametos and C. reinhardtii and those of the conventional mitochondrial rRNA genes of the green alga, Prototheca wickerhamii, and land plants and (2) to attempt to deduce the evolutionary pathways that gave rise to the unusual mitochondrial rRNA gene structures in the genus Chlamydomonas. Although phylogenetic analysis revealed an affiliation between the mitochondrial rRNA gene sequences of the two Chlamydomonas taxa to the exclusion of all other mitochondrial rRNA gene sequences tested, no specific affiliation was noted between the Chlamydomonas sequences and P. wickerhamii or land plants. Calculations of the minimal number of transpositions required to convert hypothetical ancestral rRNA gene organizations to the arrangements observed for C. eugametos and C. reinhardtii mitochondrial rRNA genes, as well as a limited survey of the size of mitochondrial rRNAs in other members of the genus, lead us to propose that the last common ancestor of Chlamydomonas algae contained fragmented mitochondrial rRNA genes that were nearly co-linear with conventional rRNA genes.     

Molecular studies on two variant repeat types of the common cetacean DNA satellite of the sperm whale, and the relationship between Physeteridae (sperm whales) and Ziphiidae (beaked whales)

In the sperm whale (Physeter macrocephalus) two different repeat types (A and B) of the common cetacean DNA satellite were identified. The evolution of each group of repeats appears to be independent from that of the other. The sequence similarity between the two groups is less than the similarity between group A and repeats of the satellite in related whale species. The systematic relationship within and between the families Physeteridae (sperm whales) and Ziphiidae (beaked whales) was addressed by both sequence analysis of the satellite and comparisons with the families Delphinidae and Phocoenidae. The mysticete blue whale (Balaenoptera musculus) was used as an outgroup in the comparisons. The molecular phylogeny, when maximum-parsimony analysis and the neighbor-joining method were used, grouped together species of each family. At the family level the ziphiids grouped closet to the families Phocoenidae and Delphinidae. The similarities between the common cetacean satellite of the blue whale and the sperm whale were greater than those between the blue whale and the other odontocetes included, suggesting that the evolution of the satellite is slower in the sperm whale than in the other odontocetes.     

Gene rearrangements in snake mitochondrial genomes: highly concerted evolution of control-region-like sequences duplicated and inserted into a tRNA gene cluster

Mitochondrial DNA (mtDNA) regions corresponding to two major tRNA gene clusters were amplified and sequenced for the Japanese pit viper, himehabu. In one of these clusters, which in most vertebrates characterized to date contains three tightly connected genes for tRNA(Ile), and tRNA(Gln), and tRNA(Met), a sequence of approximately 1.3 kb was found to be inserted between the genes for tRNA(Ile) and tRNA(Gln). The insert consists of a control-region-like sequence possessing some conserved sequence blocks, and short flanking sequences which may be folded into tRNA(Pro), tRNA(Phe), and tRNA(Leu) genes. Several other snakes belonging to different families were also found to possess a control-region-like sequence and tRNA(Leu) gene between the tRNA(Ile)and tRNA(Gln) genes. We also sequenced a region surrounded by genes for cytochrome b and 12S rRNA, where the control region and genes for tRNA(Pro) and tRNA(Phe) are normally located in the mtDNAs of most vertebrates. In this region of three examined snakes, a control-region-like sequence exists that is almost completely identical to the one found between the tRNA(Ile) and tRNA(Gln) genes. The mtDNAs of these snakes thus possess two nearly identical control-region-like sequences which are otherwise divergent to a large extent between the species. These results suggest that the duplicate state of the control-region-like sequences has long persisted in snake mtDNAs, possibly since the original insertion of the control-region-like sequence and tRNA(Leu) gene into the tRNA gene cluster, which occurred in the early stage of the divergence of snakes. It is also suggested that the duplicated control-region-like sequences at two distant locations of mtDNA have evolved concertedly by a mechanism such as frequent gene conversion. The secondary structures of the determined tRNA genes point to the operation of simplification pressure on the T psi C arm of snake mitochondrial tRNAs.     

Physical and genetic map of the mitochondrial genome of Cryphonectria parasitica Ep155

In the chestnut-blight fungus, Cryphonectria parasitica, a cytoplasmically transmissible (infectious) form of hypovirulence is associated with mitochondrial DNA (mtDNA) mutations that cause respiratory deficiencies. To facilitate the characterization of such mutations, a restriction map including the probable location of 13 genes was constructed for a relatively well-characterized virulent strain of the fungus, Ep155. The physical map is based on the order of all fragments generated by cleavage of the mtDNA by the PstI restriction endonuclease and includes some of the cleavage sites for HindIII, EcoRI, and XbaI. It was constructed from hybridization patterns of cloned mtDNA fragments with Southern blots of mtDNA digested with the four restriction enzymes. On this map, the probable locations of genes commonly found in the mitochondrial genomes of ascomycetes were determined by low-stringency hybridization of cloned Neurospora crassa mitochondrial gene probes to Southern blots of C. parasitica mtDNA. The data indicate that the mtDNA of strain Ep155 is a circular molecule of approximately 157 kbp and ranks among the largest mitochondrial chromosomes observed so far in fungi. The mtDNAs of 11 different C. parasitica isolates range in size from 135 to 157 kbp and in relatedness from 68 to 100 percent, as estimated from restriction-fragment polymorphisms. In addition to the typical mtDNA, the mitochondria of some isolates of the fungus contain double-stranded DNA plasmids consisting of nucleotide sequences not represented in the mtDNA of Ep155.     

New Bradyrhizobium japonicum strains that possess high copy numbers of the repeated sequence RS alpha

In a survey of DNA fingerprints of indigenous Bradyrhizobium japonicum with the species-specific repeated sequences RS alpha and RS beta, 21 isolates from three field sites showed numerous RS-specific hybridization bands. The isolates were designated highly reiterated sequence-possessing (HRS) isolates, and their DNA hybridization profiles were easily distinguished from the normal patterns. Some HRS isolates from two field sites possessed extremely high numbers of RS alpha copies, ranging from 86 to 175 (average, 128), and showed shifts and duplications of nif- and hup-specific hybridization bands. The HRS isolates exhibited slower growth than normal isolates, although no difference in symbiotic properties was detected between the HRS and normal isolates. Nucleotide sequence analysis of 16S rRNA genes showed that HRS isolates were strains of B. japonicum. There was no difference in the spectra of serological and hydrogenase groupings of normal and HRS isolates. Some HRS isolates possessed a tandem repeat RS alpha dimer that is similar to the structure of (IS30)2, which was shown to cause a burst of transpositional rearrangements in Escherichia coli. The results suggest that HRS isolates are derived from normal isolates in individual fields by genome rearrangements that may be mediated by insertion sequences such as RS alpha.     

Molecular evolution of Mycoplasma capricolum subsp. capripneumoniae strains, based on polymorphisms in the 16S rRNA genes

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.     

Phylogenetic analysis of Acinetobacter strains based on the nucleotide sequences of gyrB genes and on the amino acid sequences of their products

Partial nucleotide sequences of the gyrB genes (DNA gyrase B subunit genes) of 15 Acinetobacter strains, including the type and reference strains of genomic species 1 to 12 (A. calcoaceticus [genomic species 1], A. baumannii [genomic species 2], Acinetobacter genomic species 3, A. haemolyticus [genomic species 4], A. junii [genomic species 5], Acinetobacter genomic species 6, A. johnsonii [genomic species 7], A. lwoffii [genomic species 8], Acinetobacter genomic species 9, Acinetobacter genomic species 10, Acinetobacter genomic species 11, and A. radioresistens [genomic species 12]), were determined by sequencing the PCR-amplified fragments of gyrB. The gyrB sequence homology among these Acinetobacter strains ranged from 69.6 to 99.7%. A phylogenetic analysis, using the gyrB sequences, indicates that genomic species 1, 2, and 3 formed one cluster (87.3 to 90.3% identity), while genomic species 8 and 9 formed another cluster (99.7% identity). These results are consistent with those of DNA-DNA hybridization and of biochemical systematics. On the other hand, the topology of the published phylogenetic tree based on the 16S rRNA sequences of the Acinetobacter strains was quite different from that of the gyrB-based tree. The numbers of substitution in the 16S rRNA gene sequences were not high enough to construct a reliable phylogenetic tree. The gyrB-based analysis indicates that the genus Acinetobacter is highly diverse and that a reclassification of this genus would be required.     

Phylogenetic differences between particle-associated and planktonic ammonia-oxidizing bacteria of the beta subdivision of the class Proteobacteria in the Northwestern Mediterranean Sea

The aim of this study was to determine if there were differences between the types of ammonia-oxidizing bacteria of the beta subdivision of the class Proteobacteria associated with particulate material and planktonic samples obtained from the northwestern Mediterranean Sea. A nested PCR procedure performed with ammonia oxidizer-selective primers was used to amplify 16S rRNA genes from extracted DNA. The results of partial and full-length sequence analyses of 16S rRNA genes suggested that different groups of ammonia-oxidizing bacteria were associated with the two sample types. The particle-associated sequences were predominantly related to Nitrosomonas eutropha, while the sequences obtained from the planktonic samples were related to a novel marine Nitrosospira group (cluster 1) for which there is no cultured representative yet. A number of oligonucleotide probes specific for different groups of ammonia oxidizers were used to estimate the relative abundance of sequence types in samples of clone libraries. The planktonic libraries contained lower proportions of ammonia oxidizer clones (0 to 26%) than the particulate material libraries (9 to 83%). Samples of the planktonic and particle-associated libraries showed that there were depth-related differences in the ammonia oxidizer populations, with the highest number of positive clones in the particle-associated sample occurring at a depth of 700 m. The greatest difference between planktonic and particle-associated populations occurred at a depth of 400 m, where only 4% of the clones in the planktonic library were identified as Nitrosomonas clones, while 96% of these clones were identified as clones that were related to the marine Nitrosospira species. Conversely, all ammonia oxidizer-positive clones obtained from the particle-associated library were members of the Nitrosomonas group. This is the first indication that Nitrosomonas species and Nitrosospira species may occupy at least two distinct environmental niches in marine environments. The occurrence of these groups in different niches may result from differences in physiological properties and, coupled with the different environmental conditions associated with these niches, may lead to significant differences in the nature and rates of nitrogen cycling in these environments.     

Species-specific segregation of gender-associated mitochondrial DNA types in an area where two mussel species (Mytilus edulis and M. trossulus) hybridize

In each of the mussel species Mytilus edulis and M. trossulus there exist two types of mtDNA, the F type transmitted through females and the M type transmitted through males. Because the two species produce fertile hybrids in nature, F and M types of one may introgress into the other. We present the results from a survey of a population in which extensive hybridization occurs between these two species. Among specimens classified as "pure" M. edulis or "pure" M. trossulus on the basis of allozyme analysis, we observed no animal that carried the F or the M mitotype of the other species. In most animals of mixed nuclear background, an individual's mtDNA came from the species that contributed the majority of the individual's nuclear genes. Most importantly, the two mtDNA types in post-F1 male hybrids were of the same species origin. We interpret this to mean that there are intrinsic barriers to the exchange of mtDNA between these two species. Because such barriers were not noted in other hybridizing species pairs (many being even less interfertile than M. edulis and M. trossulus), their presence in Mytilus could be another feature of the unusual mtDNA system in this genus.     

Osseointegrated endosseous implants, University of Liège concepts. Various clinical applications

The occurrence of optic neuropathy in patients with MS-like disorders who carry one of the pathogenetically significant LHON mutations as well as the higher incidence of maternal transmission in familial cases of MS support the hypothesis that mitochondrial genes may be implicated in susceptibility to MS. We sequenced the entire mtDNA of six children with MS who developed optic neuritis as early and prominent visual involvement. The analysis revealed a high degree of nucleotide variations relative to the standard mtDNA sequence. After excluding various synonymous nucleotide changes and common neutral polymorphisms, eight discrete novel missense mutations within the protein coding, tRNA or rRNA genes were detected. None of the eight polymorphic sites were found in common between the patients with MS. Of particular interest was the observation that five of six children carried a total of nine secondary LHON mutations at nucleotide positions 4216, 4917 or 13708. We conclude that variation in mtDNA is unlikely to contribute to genetic predisposition for MS. However, secondary LHON mutations may be regarded as additional risk factor for developing prominent optic nerve involvement. The association of individual sets of mtDNA variations with phenotypic presentation in certain subgroups of MS patients remains to be clarified.     

Evolutionary relationships among members of the genus Chlamydia based on 16S ribosomal DNA analysis

Nucleotide sequences from strains of the four species currently in the genus Chlamydia, C. pecorum, C. pneumoniae, C. psittaci, and C. trachomatis were investigated. In vitro-amplified RNA genes of the ribosomal small subunit from 30 strains of C. pneumoniae and C. pecorum were subjected to solid-phase DNA sequencing of both strands. The human isolates of C. pneumoniae differed in only one position in the 16S rRNA gene, indicating genetic homogeneity among these strains. Interestingly, horse isolate N16 of C. pneumoniae was found to be closely related to the human isolates of this species, with a 98.9% nucleotide similarity between their 16S rRNA sequences. The type strain and koala isolates of C. pecorum were also found to be very similar to each other, possessing two different 16S rRNA sequences with only one-nucleotide difference. Furthermore, the C. pecorum strains truncated the 16S rRNA molecule by one nucleotide compared to the molecules of the other chlamydial species. This truncation was found to result in loss of a unilaterally bulged nucleotide, an attribute present in all other eubacteria. The phylogenetic structure of the genus Chlamydia was determined by analysis of 16S rRNA sequences. All phylogenetic trees revealed a distinct line of descent of the family Chlamydiaceae built of two main clusters which we denote the C. pneumoniae cluster and the C. psittaci cluster. The clusters were verified by bootstrap analysis of the trees and signature nucleotide analysis. The former cluster contained the human isolates of C. pneumoniae and equine strain N16. The latter cluster consisted of C. psittaci, C. pecorum, and C. trachomatis. The members of the C. pneumoniae cluster showed tight clustering and strain N16 is likely to be a subspecies of C. pneumoniae since these strains also share some antigenic cross-reactivity and clustering of major outer membrane protein gene sequences. C. psittaci and strain N16 branched early out of the respective cluster, and interestingly, their inclusion bodies do not stain with iodine. Furthermore, they also share less reliable features like normal elementary body morphology and plasmid content. Therefore, the branching order presented here is very likely a true reflection of evolution, with strain N16 of the species C. pneumoniae and C. psittaci forming early branches of their respective cluster and with C. trachomatis being the more recently evolved species within the genus Chlamydia.     

Origin of a fungal symbiont of perennial ryegrass by interspecific hybridization of a mutualist with the ryegrass choke pathogen, Epichlo? typhina

Seed-borne fungal symbionts (endophytes) provide many cool-season grass species with biological protection from biotic and abiotic stresses. The endophytes are asexual, whereas closely related sexual species of genus Epichlo? (Clavicipitales) cause grass choke disease. Perennial ryegrass (Lolium perenne) is a host of two endophyte taxa, LpTG-1 (L. perenne endophyte taxonomic grouping one = Acremonium lolii) and LpTG-2, as well as the choke pathogen, Epichlo? typhina (represented by isolate E8). Relationships among these fungi and other Epichlo? species were investigated by analysis of gene sequences, DNA polymorphisms and allozymes. The results indicate that LpTG-2 is a heteroploid derived from an interspecific hybrid. The LpTG-2 isolates had two copies each of nine out of ten genes analyzed (the exception being the rRNA gene locus), and the profiles for seven of these were composites of those from E. typhina E8 and A. lolii isolate Lp5. Molecular phylogenetic analysis grouped the two beta-tubulin genes of LpTG-2 into separate clades. One (tub2-1) was related to that of E. typhina E8, and the other (tub2-2) to that of A. lolii. The mitochondrial DNA profile of LpTG-2 was similar to that of A. lolii, but its rRNA gene sequence grouped it with E. typhina E8. A proposed model for the evolution of LpTG-2 involves infection of a L. perenne-A. lolii symbiotum by E. typhina, followed by hybridization of the two fungi. Such interspecific hybridization may be a common and important mechanism for genetic variation in Epichlo? endophytes.     

rpoB sequence analysis as a novel basis for bacterial identification

Comparison of the sequences of conserved genes, most commonly those encoding 16S rRNA, is used for bacterial genotypic identification. Among some taxa, such as the Enterobacteriaceae, variation within this gene does not allow confident species identification. We investigated the usefulness of RNA polymerase beta-subunit encoding gene (rpoB) sequences as an alternative tool for universal bacterial genotypic identification. We generated a database of partial rpoB for 14 Enterobacteriaceae species and then assessed the intra- and interspecies divergence between the rpoB and the 16S rRNA genes by pairwise comparisons. We found that levels of divergence between the rpoB sequences of different strains were markedly higher than those between their 16S rRNA genes. This higher discriminatory power was further confirmed by assigning 20 blindly selected clinical isolates to the correct enteric species on the basis of rpoB sequence comparison. Comparison of rpoB sequences from Enterobacteriaceae was also used as the basis for their phylogenetic analysis and demonstrated the genus Klebsiella to be polyphyletic. The trees obtained with rpoB were more compatible with the currently accepted classification of Enterobacteriaceae than those obtained with 16S rRNA. These data indicate that rpoB is a powerful identification tool, which may be useful for universal bacterial identification.     

Evaluation of intraspecies genetic variation within the 16S rRNA gene of Mycoplasma hominis and detection by polymerase chain reaction

Mycoplasma hominis is a heterogeneous species with DNA-DNA hybridization values ranging from 51 to 100%. We report here the sequencing of the 16S rRNA gene of a strain (183) that greatly differs from the type strain (PG21) of this species. Comparison of 16S rDNA sequences from these two strains showed limited differences, indicating that the two strains belong to the same rRNA species complex. Using these nucleotide sequence data, we established a rapid method for the detection of M. hominis by using polymerase chain reaction. This method was shown to be sensitive and specific when tested with reference strains and clinical isolates.     

16S-23S and 23S-5S intergenic spacer regions of lactobacilli: nucleotide sequence, secondary structure and comparative analysis

Lactobacilli have been used as industrial starters for a long time, but in several cases their identification was, and still is, neither easy nor reliable. The aim of the present work was to examine whether the intergenic spacer regions could be of value in the identification of Lactobacillus species. For that purpose, the polymerase chain reaction (PCR) was used to amplify 16S-23S and 23S-5S spacer regions of Lactobacillus (L.) acidophilus, L. delbrueckii subsp. bulgaricus, L. casei, L. helveticus and L. curvatus. The PCR products were directly sequenced, and two forms of ribosomal RNA (rrn) operons were identified in each species studied: one with tandem tRNA(Ile)/tRNA(Ala) genes and the other one without tRNA genes. Our study revealed that the rrn operons of Lactobacillus species studied comprise the genes of 16S, 23S and 5S rRNA, in that order. Only the tRNA genes and the rRNA processing stems are highly conserved in spacer regions of lactobacilli. The divergence between the lactobacilli spacer region sequences arises from insertions and deletions of short sequences. These sequences could be interesting candidates for the development of species-specific probes. Theoretical RNA/RNA secondary structure models of the interaction between the two spacer region sequences were constructed. In conclusion, the two spacer region sequences may prove to be a useful alternative to 16S and 23S rDNA sequencing for designing species-specific probes and for establishing phylogenetic relationships between closely related species such as L. curvatus and L. casei or L. acidophilus and L. helveticus.