Structure and evolution of the avian mitochondrial control region

The structural and evolutionary characteristics of the mitochondrial control region were studied by using control region sequences of 68 avian species. The distribution of the variable nucleotide positions within the control region was found to be genus specific and not dependant on the level of divergence, as suggested before. Saturation was shown to occur at the level of divergence of 10% in pairwise comparisons of the control region sequences, as has also been reported for the third codon positions in ND2 and cytochrome b genes of mtDNA. The ratio of control region vs cytochrome b divergence in pairwise comparisons of the sequences was shown to vary from 0.13 to 21.65, indicating that the control region is not always the most variable region of the mtDNA, but also that there are differences in the rate of divergence among the lineages. Only two of the conserved sequence blocks localized earlier for other species, D box and CSB-1, were found to show a considerable amount of sequence conservation across the avian and mammalian sequences. Additionally, a novel avian-specific sequence block was found.     

Sequence and properties of the human KB cell and mouse L cell D-loop regions of mitochondrial DNA.

The sequences of the displacement-loop (D-loop) regions of mitochondrial DNA (mtDNA) from mouse L cells and human KB cells have been determined and provide physical maps to aid in the identification of sequences involved in the regulation of replication and expression of mammalian mtDNA. Both D-loop regions are bounded by the genes for tRNAPhe and tRNAPro. This region contains the most highly divergent sequences in mtDNA with the exceptions of three small conserved sequence blocks near the 5' ends of D-loop strands, a 225 nucleotide conserved sequence block in the center of the D-loop strand template region, and a short sequence associated with the 3' ends of D-loop strands. A sequence similar to that associated with the 3' termini of D-loop strands overlaps one of the conserved sequence blocks near the 5' ends of D-loop strands. The large, central conserved sequence probably does not code for a protein since no open reading frames are discretely conserved. Numerous symmetric sequences and potential secondary structures exist in these sequences, but none appear to be clearly conserved between species.     

The main regulatory region of mammalian mitochondrial DNA: Structure-function model and evolutionary pattern

Summary The evolution of the main regulatory region (D-loop) of the mammalian mitochondrial genome was analyzed by comparing the sequences of eight mammalian species: human, common chimpanzee, pygmy chimpanzee, dolphin, cow, rat, mouse, and rabbit. The best alignment of the sequences was obtained by optimization of the sequence similarities common to all these species.The two peripheral left and right D-loop domains, which contain the main regulatory elements so far discovered, evolved rapidly in a species-specific manner generating heterogeneity in both length and base composition. They are prone to the insertion and deletion of elements and to the generation of short repeats by replication slippage. However, the preservation of some sequence blocks and similar cloverleaf-like structures in these regions, indicates a basic similarity in the regulatory mechanisms of the mitochondrial genome in all mammalian species.We found, particularly in the right domain, significant similarities to the telomeric sequences of the mitochondrial (mt) and nuclear DNA ofTetrahymena thermophila. These sequences may be interpreted as relics of telomeres present in ancestral linear forms of mtDNA or may simply represent efficient templates of RNA primase-like enzymes.Due to their peculiar evolution, the two peripheral domains cannot be used to estimate in a quantitative way the genetic distances between mammalian species. On the other hand the central domain, highly conserved during evolution, behaves as a good molecular clock.Reliable estimates of the times of divergence between closely and distantly related species were obtained from the central domain using a Markov model and assuming nonhomogeneous evolution of nucleotide sites.     

Efficient site-specific cleavage by RNase MRP requires interaction with two evolutionarily conserved mitochondrial RNA sequences.

RNase MRP is a site-specific endonuclease that processes primer mitochondrial RNA from the leading-strand origin of mitochondrial DNA replication. Using deletional analysis and saturation mutagenesis, we have determined the substrate requirements for cleavage by mouse mitochondrial RNase MRP. Two regions of sequence homology among vertebrate mitochondrial RNA primers, conserved sequence blocks II and III, were found to be critical for both efficient and accurate cleavage; a third region of sequence homology, conserved sequence block I, was dispensable. Analysis of insertion and deletion mutations within conserved sequence block II demonstrated that the specificity of RNase MRP accommodates the natural sequence heterogeneity of conserved sequence block II in vivo. Heterologous assays with human RNase MRP and mutated mouse mitochondrial RNA substrates indicated that sequences essential for substrate recognition are conserved between mammalian species.     

Rates of nuclear and cytoplasmic mitochondrial DNA sequence divergence in mammals

Differential rates of nucleotide substitution among different gene segments and between distinct evolutionary lineages is well documented among mitochondrial genes and is likely a consequence of locus-specific selective constraints that delimit mutational divergence over evolutionary time. We compared sequence variation of 18 homologous loci (15 coding genes and 3 parts of the control region) among 10 mammalian mitochondrial DNA genomes which allowed us to describe different mitochondrial evolutionary patterns and to produce an estimation of the relative order of gene divergence. The relative rates of divergence of mitochondrial DNA genes in the family Felidae were estimated by comparing their divergence from homologous counterpart genes included in nuclear mitochondrial DNA (Numt, pronounced "new might"), a genomic fossil that represents an ancient transfer of 7.9 kb of mitochondrial DNA to the nuclear genome of an ancestral species of the domestic cat (Felis catus). Phylogenetic analyses of mitochondrial (mtDNA) sequences with multiple outgroup species were conducted to date the ancestral node common to the Numt and the cytoplasmic (Cymt) mtDNA genes and to calibrate the rate of sequence divergence of mitochondrial genes relative to nuclear homologous counterparts. By setting the fastest substitution rate as strictly mutational, an empirical "selective retardation index" is computed to quantify the sum of all constraints, selective and otherwise, that limit sequence divergence of mitochondrial gene sequences over time.      

Human Mitochondrial DNA Variation and Evolution: Analysis of Nucleotide Sequences from Seven Individuals

We have analyzed nucleotide sequence variation in an approximately 900-base pair region of the human mitochondrial DNA molecule encompassing the heavy strand origin of replication and the D-loop. Our analysis has focused on nucleotide sequences available from seven humans. Average nucleotide diversity among the sequences is 1.7%, several-fold higher than estimates from restriction endonuclease site variation in mtDNA from these individuals and previously reported for other humans. This disparity is consistent with the rapidly evolving nature of this noncoding region. However, several instances of convergent or parallel gain and loss of restriction sites due to multiple substitutions were observed. In addition, other results suggest that restriction site (as well as pairwise sequence) comparisons may underestimate the total number of substitutions that have occurred since the divergence of two mtDNA sequences from a common ancestral sequence, even at low levels of divergence. This emphasizes the importance of recognizing the large standard errors associated with estimates of sequence variability, particularly when constructing phylogenies among closely related sequences. Analysis of the observed number and direction of substitutions revealed several significant biases, most notably a strand dependence of substitution type and a 32-fold bias favoring transitions over transversions. The results also revealed a significantly nonrandom distribution of nucleotide substitutions and sequence length variation. Significantly more multiple substitutions were observed than expected for these closely related sequences under the assumption of uniform rates of substitution. The bias for transitions has resulted in predominantly convergent or parallel changes among the observed multiple substitutions. There is no convincing evidence that recombination has contributed to the mtDNA sequence diversity we have observed.     

Nucleotide sequence of a Xenopus laevis mitochondrial DNA fragment containing the D-loop, flanking tRNA genes and the apocytochrome b gene

Extensive corrections of the nucleotide sequence of the Xenopus laevis mitochondrial (mt) displacement (D) loop and surrounding genes [Wong et al., Nucl. Acids Res. 11 (1983) 4977-4995] are reported, including addition of two stretches of nucleotides and 60 scattered modifications. The additional sequences presented here correspond to the apocytochrome b gene, the tRNAGlu gene and part of URF6. This allows us to propose a conformational model for the X. laevis apocytochrome b protein and also permits comparisons with mammalian mtDNA. The D-loop sequence is poorly conserved except for sequences involved in the regulation of the mt genome (conserved sequence blocks and the DNA polymerase stop sequences). On the other hand, all genes show marked conservation both of their nucleotide sequence and their respective location on the mt genome. Organization of the genetic information described for mammalian mtDNA also holds for the X. laevis mtDNA. This result strongly suggests that all animal vertebrate mtDNAs have followed the same evolutionary pathway.     

Origin and evolution of homologous repeated sequences in the mitochondrial DNA control region of shrews

The complete mitochondrial DNA (mtDNA) control region was amplified and directly sequenced in two species of shrew, Crocidura russula and Sorex araneus (Insectivora, Mammalia). The general organization is similar to that found in other mammals: a central conserved region surrounded by two more variable domains. However, we have found in shrews the simultaneous presence of arrays of tandem repeats in potential locations where repeats tend to occur separately in other mammalian species. These locations correspond to regions which are associated with a possible interruption of the replication processes, either at the end of the three-stranded D-loop structure or toward the end of the heavy-strand replication. In the left domain the repeated sequences (R1 repeats) are 78 bp long, whereas in the right domain the repeats are 12 bp long in C. russula and 14 bp long in S. araneus (R2 repeats). Variation in the copy number of these repeated sequences results in mtDNA control region length differences. Southern blot analysis indicates that level of heteroplasmy (more than one mtDNA form within an individual) differs between species. A comparative study of the R2 repeats in 12 additional species representing three shrew subfamilies provides useful indications for the understanding of the origin and the evolution of these homologous tandemly repeated sequences. An asymmetry in the distribution of variants within the arrays, as well as the constant occurrence of shorter repeated sequences flanking only one side of the R2 arrays, could be related to asymmetry in the replication of each strand of the mtDNA molecule. The pattern of sequence and length variation within and between species, together with the capability of the arrays to form stable secondary structures, suggests that the dominant mechanism involved in the evolution of these arrays in unidirectional replication slippage.      

Assessing Nuclear and Mitochondrial DNA Sequence Variation Within Steinernema (Rhabditida: Steinernematidae)

DNA sequence analysis was used to characterize the nuclear ribosomal DNA ITS1 region and a portion of the COII and 16S rDNA genes of the mitochondrial genome from Steinernema entomopathogenic nematodes. Nuclear ITS1 nucleotide divergence among seven Steinernema spp. ranged from 6 to 22%, and mtDNA divergence among five species ranged from 12 to 20%. No intraspecific variation was observed among three S. feltiae strains. Phylogenetic analysis of both nuclear and mitochondrial DNA sequences confirms the existing morphological relationships of several Steinernema species. Both the rDNA ITS1 and mtDNA sequences were useful for resolving relationships among Steinernema taxa.     

Molecular analyses of mtDNA deletion mutations in microdissected skeletal muscle fibers from aged rhesus monkeys

Mitochondrial DNA (mtDNA) deletion mutations co-localize with electron transport system (ETS) abnormalities in rhesus monkey skeletal muscle fibers. Using laser capture microdissection in conjunction with PCR and DNA sequence analysis, mitochondrial genomes from single sections of ETS abnormal fibers were characterized. All ETS abnormal fibers contained mtDNA deletion mutations. Deletions were large, removing 20-78% of the genome, with some to nearly all of the functional genes lost. In one-third of the deleted genomes, the light strand origin was deleted, whereas the heavy strand origin of replication was conserved in all fibers. A majority (27/39) of the deletion mutations had direct repeat sequences at their breakpoints and most (36/39) had one breakpoint within or in close proximity to the cytochrome b gene. Several pieces of evidence support the clonality of the mtDNA deletion mutation within an ETS abnormal region of a fiber: (a) only single, smaller than wild-type, PCR products were obtained from each ETS abnormal region; (b) the amplification of mtDNA from two regions of the same ETS abnormal fiber identified identical deletion mutations, and (c) a polymorphism was observed at nucleotide position 16103 (A and G) in the wild-type mtDNA of one animal (sequence analysis of an ETS abnormal region revealed that mtDNA deletion mutations contained only A or G at this position). Species-specific differences in the regions of the genomes lost as well as the presence of direct repeat sequences at the breakpoints suggest mechanistic differences in deletion mutation formation between rodents and primates.     

Intraspecific nucleotide sequence differences in the major noncoding region of human mitochondrial DNA.

Nucleotide sequences of the major noncoding region of human mitochondrial DNA (mtDNA) from 95 human placentas have been determined. These sequences include at least a 482-bp-long region encompassing most of the D-loop-forming region. Comparisons of these sequences with those previously determined have revealed remarkable features of nucleotide substitutions and insertion/deletion events. The nucleotide diversity among the sequences is estimated as 1.45%, which is three- to fourfold higher than the corresponding value estimated from restriction-enzyme analysis of whole mtDNA genome. A hypervariable region has also been defined. In this 14-bp region, 17 different sequences were detected. More than 97% of the base changes are transitions. A significantly nonrandom distribution of nucleotide substitutions and sequence length variations were also noted. The phylogenetic analysis indicates that diversity among the negroids is much larger than that among the caucasoids or the mongoloids. In fact, part of the negroids first diverged from other humans in the phylogenetic tree. A striking finding in the phylogenetic analysis is that the mongoloids can be separated into two distinct groups. Divergence of part of the mongoloids follows the earliest divergence of part of the negroids. The remainder of the mongoloids subsequently diverged together with the caucasoids. This observation confirmed our earlier study, which clearly demonstrated, by the restriction-enzyme analysis, existence of two distinct groups in the Japanese.     

Three divergent mitochondrial genomes from California populations of the copepod Tigriopus californicus

Previous work on the harpacticoid copepod Tigriopus californicus has focused on the extensive population differentiation in three mtDNA protein coding genes (COXI, COXII, Cytb). In order to get a more complete understanding of mtDNA evolution in this species, we sequenced three complete mitochondrial genomes (one from each of three California populations) and compared them to two published mtDNA genomes from an Asian congener, Tigriopus japonicus. Several features of the mtDNA genome appear to be conserved within the genus: 1) the unique order of the protein coding genes, rRNA genes and most of the tRNA genes, 2) the genome is compact, varying between 14.3 and 14.6 kb, and 3) all genes are encoded on the same strand of the mtDNA. Within T. californicus, extremely high levels of nucleotide divergence (>20%) are observed across much of the mitochondrial genome. Inferred amino acid sequences of the proteins encoded in the mtDNAs also show high levels of divergence; at the extreme, the three ND3 variants in T. californicus showed >25% amino acid substitutions, compared with <3% amino acid divergence at the previously studied COXI locus. Unusual secondary structures make functional assignments of some tRNAs difficult. The only apparent tRNA(trp) in these genomes completely overlaps the 5' end of the 16S rRNA in all three T. californicus mtDNAs. Although not previously noted, this feature is also conserved in T. japonicus mtDNAs; whether this sequence is processed into a functional tRNA has not been determined. The putative control region contains a duplicated segment of different length (from 88 to 155 bp) in each of the T. californicus sequences. In each case, the duplicated segments are not tandem repeats; despite their different lengths, the distance between the start of the first and the start of the second repeat is conserved (520 bp). The functional significance, if any, of this repeat structure remains unknown.     

STR polymorphism of mtDNA D-loop in rhesus macaques of Bangladesh

Molecular variation of mitochondrial DNA (mtDNA) was investigated for rhesus macaques (Macaca mulatta) of Bangladesh. A partial sequence (583–599 bp) of mtDNA containing the second variable region of the D-loop was compared for 39 individuals from five localities in the country. A total of seven haplotypes were detected with substitutional or insertion/deletion mutations. They contained a unique polymorphism of pentanucleotide STRs (short tandem repeats). There were at least four different length types, from two to five repeats of the unit nucleotide. One site of substitution and one site of single nucleotide insertion/deletion were also involved in the polymorphism. The mutation hot spots of the STR polymorphism were located between the first and second conserved sequence blocks (CSB1 and CSB2), as observed previously in some other mammals. The geographical distribution of the STR polymorphism revealed local differences; the northeastern population was polymorphic with three STR haplotypes, but other local populations were simply monomorphic with a single STR haplotype. Molecular phylogenetic analysis with reported sequences from outside Bangladesh indicated a low substitution diversity of mtDNA in Bangladesh. Clustering results suggested a close relationship to India and divergence from Laos and China.     

鸮形目4种鸟类线粒体调控区全序列的测定与比较研究

利用Long-PCR和Primer Walking的方法对鸮形目的短耳鸮、长耳鸮、纵纹腹小鸮、灰林鸮4种鸟类的线粒体调控区进行了全序列测定。结果表明：短耳鸮的调控区跃度为3290bp;长耳鸮为2848bp;纵纹腹小鸮为2444bp;灰林鸮为1771bp。短耳鸮的调控区长度是4种鸮中最大的,并且是目前已知最大的鸟类线粒体调控区。这4种鸮类调控区的基本结构和其他鸟类相似,按照碱基变化速率的不同可以分为3个区：碱基变化速率较快的外围区域Ⅰ、Ⅲ和保守的中间区域Ⅱ。这4种鸟类调控区的3’端均存在大量的串联重复序列,短耳鸮为126bp单元重复7次和78bp单元重复14次;长耳鸮为127bp单元重复8次和78bp单几重复6次;纵纹腹小鸮有3个重复单元,分别为89bp单元重复3次、77bp单元重复4次和71bp单元重复6次;灰林鸮仅有1个单元的串联重复为78bp重复5次。调控区中串联重复序列可能是由链的滑动错配产生,另外这些重复序列都能形成热力学稳定的多重茎环二级结构,而且在重复序列中还发现一些保守基序,这说明重复序列可能具有一定的生理功能,影响调控区的调重控功能从而影响线粒体基因组的复制和转录。     

Sequence length and variation in the mitochondrial control region of two freshwater gobiid fishes belonging to Rhinogobius (Teleostei: Gobioidei)

The control region (D-loop) of mitochondrial DNA (mtDNA) was amplified and sequenced for eight samples of the rhinogobies Rhinogobius maculafasciatus and R. giurinus from Taiwan and southern China. The control regions of both species are of 841–842 bp; the length of these sequences being the most compact among all known sequences in teleost fishes. Three conserved sequence blocks (CSB) were observed. The full D-loop and tRNA Phe gene sequences were determined and compared with other fishes. The interspecific sequence divergence between the two species is 11.3–11.7%; and the intraspecific variation in R. guirinus 0.8–1.8%. Results suggest that the control region of Rhinogobius is informative for phylogenetic reconstruction at both intraspecific and interspecific levels in this gobiid genus.     

Mitochondrial DNA Sequence Variation in a Sea Star (Leptasteriasspp.) Species Complex

A 551-bp region of a PCR product containing the putative mitochondrial control region and flanking sequences was analyzed for sequence variation among 19 sea stars representing 10 previously described PCR–RFLP haplotypes within a cryptic species complex (Leptasteriasspp.). Most (97%) of the sequence variation was interhaplotypic rather than intrahaplotypic, which greatly reduced the utility of sequence polymorphisms in this mtDNA region as markers of intrahaplotypic population structure and gene flow. The estimated number of transition and transversion substitutions per nucleotide site, corrected for multiple hits, was 0.0364 and 0.0158, respectively. Most of the sequence variation occurred in the first half of the putative control region. Phylogenetic analysis (both maximum parsimony and maximum likelihood) revealed three well-supported clades, but the position of two PCR-RFLP haplotypes was not completely resolved. Low intraspecific mtDNA sequence divergence over large geographic distances may be a general pattern for echinoderm species.     

Sequence variation and structural conservation in the D-loop region and flanking genes of mitochondrial DNA from Japanese pond frogs

The nucleotide sequences of the D-loop region and its flanking genes of the mitochondrial DNA (mtDNA) from Japanese pond frogs were determined by the methods of PCR, cloning, and sequencing. The frogs belonged to two species, one subspecies, and one local race. The gene arrangements adjacent to the D-loop region were analyzed. The frogs shared a unique mitochondrial gene order that was found in Rana catesbeiana; i.e., cyt b--D-loop region--tRNA(Leu(CUN))--tRNA(Thr)--tRNA(Pro)--tRNA(Phe)--12S rRNA. The arrangements of the three tRNA genes of these frogs were different from those of X. laevis, a species which has the same overall structure as in mammals. Highly repetitive sequences with repeat units (16-bp or 17-bp sequence specific for each taxon) were found in the D-loop region. The length of repetitive sequences varied from 0.6 kbp to 1.2 kbp, and caused the extensive size variation in mtDNA. Several short sequence elements such as putative TAS, OH, CSB-1, and CSB-2 were found in the D-loop region of these frogs. The sequences of these short regulatory elements were conserved in R. catesbeiana, X. laevis, and also in human. The comparison of sequence divergences of the D-loop region and its adjacent genes among various taxa revealed that the rates of nucleotide substitutions depend on genes. The nucleotide sequences of the 3'-side segment of the D-loop region were the most variable among taxa, whereas those of the tRNA and 12S rRNA genes were the most conservative.