Expansion/contraction of mammalian mitochondrial DNA repeats in Escherichia coli mimics the mitochondrial heteroplasmy.

Length polymorphism due to tandem repeats is a common feature in animal mitochondrial DNA. The rabbit mitochondrial genome contains a 20 bp repeat domain, which generates a general heteroplasmic state. The observed polymorphic patterns suggest a dynamic equilibrium between gain and loss of units that maintains the copy number in the range 3-19 repeat units. In the apparent absence of recombination, slipped-strand mispairing during replication appears to be the primary cause of additions and deletions. To investigate this hypothesis we have set up a plasmid assay in Escherichia coli. A variable number of repeat units was inserted into a plasmid in both orientations relative to the colE1 origin of replication. Our data show that (i) a minimum unit number (>3) is necessary to generate length polymorphs, (ii) the number of events increases with the length tract, (iii) an excess of additions over deletions is found when the copy number is less than 10 and the trend is reversed when it is over 10, (iv) the frequency of deletions-additions is dependent on the orientation, (v) the polymorphism patterns are different according to the orientation. The length polymorphic pattern generated in the bacteria, in one orientation, mimics that observed in the mitochondria, suggesting that slipped mispairing between repeated sequences during DNA replication is responsible for the mitochondrial heteroplasmic state.     

The mitochondrial control region of Cervidae: evolutionary patterns and phylogenetic content

The mitochondrial control region (CR) sequence, also known as the D- loop, has been determined for six Cervidae (Artiodactyla, Ruminantia): the red and fallow deers (subfamily Cervinae), the brocket deer and two roe deers (subfamily Odocoileinae), and the Chinese water deer (Hydropotinae). These new sequences have been aligned with available cervid and bovid orthologues. Comparative analyses indicate that the 5'- peripheral domain exhibits a 75-bp length polymorphism near sequences associated with the termination of the H-strand replication. The New World Odocoileinae possess the longest cervid CR due to the presence of an additional 47-bp tandem repeat, located in the 3'-peripheral domain, downstream of the initiation site for H-strand replication (OH) and the first conserved sequence block (CSB-1). This insertion represents a duplication spanning the OH to CSB-1 region and constitutes an exclusive synapomorphy for New World Odocoileinae. Phylogenetic analyses of the complete CR support the paraphyly of antlered deers due to the nesting of the antlerless Hydropotes within Odocoileinae. Capreolus is the closest relative of Hydropotes, and the divergence of this Old World Odocoileinae clade may have occurred between 8.7 and 10.4 MYA. The conserved central domain of CR can be aligned across ungulates and indicates the Pecora monophyly, their close association with cetaceans, and the earlier emergence of suiformes.      

Sequence, organization, and evolution of the A+T region of Drosophila melanogaster mitochondrial DNA

The long (4.6-kb) A+T region of Drosophila melanogaster mitochondrial DNA has been cloned and sequenced. The A+T region is organized in two large arrays of tandemly repeated DNA sequence elements, with nonrepetitive intervening and flanking sequences comprising only 22% of its length. The first repeat array consists of five repeats of 338-373 bp. The second consists of four intact 464-bp repeats and a fifth partial repeat of 137 bp. Three DNA sequence elements are found to be highly conserved in D. melanogaster and in several Drosophila species with short A+T regions. These include a 300-bp DNA sequence element that overlaps the DNA replication origin and two thymidylate stretches identified on opposite DNA strands. We conclude that the length heterogeneity observed in the A+T regulatory region in mitochondrial DNAs from the genus Drosophila results from the expansion (and contraction) of the number of repeated DNA sequence elements. We also propose that the 300-bp conserved DNA sequence element, in conjunction with another primary sequence determinant, perhaps the adjacent thymidylate stretch, functions in the regulation of mitochondrial DNA replication.      

Sequence variation and evolution of the mitochondrial DNA control region in the musk shrew, Suncus murinus

The complete mitochondrial DNA (mtDNA) control region was cloned and sequenced in the musk shrew, Suncus murinus, Insectivora. The general aspect was similar to that found in other mammals. We have found in two locations of this region the presence of arrays of tandem repeats like those in other shrew species. One array was located in the left domain containing the termination-associated sequences (TAS) and the length of a copy was 77 bp. The other repeats were situated upstream from the recognition site for the end of H-strand replication in the right domain and were 20 bp long. The left halves of the control region containing the former repeats were sequenced and compared in several laboratory lines and wild animals from different localities, variations in copy number of repeated sequences were found both among individuals and within an individual. A comparative study of repeated sequences provides useful indication for the origin and evolution of tandem repeated sequences. Strand slippage and mispairing during replication of mtDNA with concerted manner is currently regarded as a dominant theory to account molecular mechanism for tandemly repeated sequences, and the pattern of sequence and length variation in our study supports this theory. Our results, however, suggest that the evolution of the repeated sequences containing the TAS in the musk shrew might go through the process of two steps; at the first step one complete repeated and several incomplete repeated sequences had reproduced in common ancestor of the shrew, and the second stage step-up of complete repeated sequences occurred with concerted evolution after differentiation into continental and insular groups.     

Evolution of Repeated Sequence Arrays in the D-Loop Region of Bat Mitochondrial DNA

Analysis of mitochondrial DNA control region sequences from 41 species of bats representing 11 families revealed that repeated sequence arrays near the tRNA-Pro gene are present in all vespertilionine bats. Across 18 species tandem repeats varied in size from 78 to 85 bp and contained two to nine repeats. Heteroplasmy ranged from 15% to 63%. Fewer repeats among heteroplasmic than homoplasmic individuals in a species with up to nine repeats indicates selection may act against long arrays. A lower limit of two repeats and more repeats among heteroplasmic than homoplasmic individuals in two species with few repeats suggests length mutations are biased. Significant regressions of heteroplasmy, θ and π, on repeat number further suggest that repeat duplication rate increases with repeat number. Comparison of vespertilionine bat consensus repeats to mammal control region sequences revealed that tandem repeats of similar size, sequence and number also occur in shrews, cats and bighorn sheep. The presence of two conserved protein-binding sequences in all repeat units indicates that convergent evolution has occurred by duplication of functional units. We speculate that D-loop region tandem repeats may provide signal redundancy and a primitive repair mechanism in the event of somatic mutations to these binding sites.     

Length and Sequence Variation in Evening Bat D-Loop Mtdna

Length variation in D-loop mitochondrial DNA was observed after amplification with the polymerase chain reaction (PCR) in 28% of 195 evening bats, Nycticeius humeralis, from seven colonies. Nucleotide sequences of PCR products show that this heteroplasmy is characterized by an 81-bp region which is tandemly repeated five to eight times. Southern blots using PCR products as probes on HaeIII genomic digests confirm the presence of heteroplasmy. Furthermore, densitometry of electrophoresed PCR products from 109 mother-offspring pairs indicate that heteroplasmy is stably transmitted from mother to offspring with one exception: a heteroplasmic offspring had a homoplastic mother and sib. Nucleotide sequences from this family reveal that a repeat duplication and deletion occurred. The observed mutation rate per generation, mu, for length polymorphism is comparable to an independent estimate, mu = 10(-2), based on hierarchical diversity statistics. With the exception of the repeat nearest the proline tRNA gene, sequence similarities between repeats within bats are consistent with a model of concerted evolution due to unidirectional replication slippage. Selection is inferred to act on the first repeat because in comparison to other repeats it has the least sequence divergence among bats, the fewest transversions, and the lowest minimum free energy associated with folding.     

鸮形目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次。调控区中串联重复序列可能是由链的滑动错配产生,另外这些重复序列都能形成热力学稳定的多重茎环二级结构,而且在重复序列中还发现一些保守基序,这说明重复序列可能具有一定的生理功能,影响调控区的调重控功能从而影响线粒体基因组的复制和转录。     

Different mechanisms inferred from sequences of human mitochondrial DNA deletions in ocular myopathies.

We have sequenced the deletion borders of the muscle mitochondrial DNA from 24 patients with heteroplasmic deletions. The length of these deletions varies from 2.310 bp to 8.476 bp and spans from position 5.786 to 15.925 of the human mitochondrial genome preserving the heavy chain and light chain origins of replication. 12 cases are common deletions identical to the mutation already described by other workers and characterized by 13 bp repeats at the deletion boundaries, one of these repeats being retained during the deletion process. The other cases (10 out of 12) have shown deletions which have not been previously described. All these deletions are located in the H strand DNA region which is potentially single stranded during mitochondrial DNA replication. In two cases, the retained Adenosine from repeat closed to the heavy strand origin of replication would indicate slippage mispairing. Furthermore in one patient two mt DNA molecules have been cloned and their sequences showed the difference of four nucleotides in the breakpoint of the deletion, possibly dued to slippage mispairing. Taken together our results suggest that deletions occur either by slippage mispairing or by internal recombination at the direct repeat level. They also suggest that different mechanisms account for the deletions since similarly located deletions may display different motives at the boundaries including the absence of any direct repeat.     

Mechanism of mitochondrial DNA replication in mouse L-cells: RNA priming during the initiation of heavy-strand synthesis

The major form of mouse L-cell mitochondrial DNA contains a small displacement loop at the replication origin, created by synthesis of a 550 to 670-nucleotide portion of the heavy strand. These short heavy-strand segments remain hydrogen-bonded to the parental light strand and are collectively termed 7 S mitochondrial DNA. The unique location of these 7 S mitochondrial DNAs at the heavy-strand origin suggests that they may function as primers in the synthesis of full-length heavy strands. Ribonucleotides have been detected at the 5′-end of some of these molecules, which are most likely remnants of primer RNAs. Using 5′-end labeling in vitro, we have determined that these ribonucleotides occur at several discrete positions along the nucleotide sequence of the origin region, which suggests that there may be variability in the precise initiation point of RNA priming or in the location of the switchover from RNA priming to DNA synthesis. The length of 5′-end RNA was estimated by alkali treatment of mitochondrial DNA prior to end labeling. A range of one to ten ribonucleotides was hydrolyzed from the 5′-end of some 7 S mitochondrial DNA strands. This is the first evidence of RNA priming at a eukaryotic cell DNA replication origin.     

Variation in repeat length and heteroplasmy of the mitochondrial DNA control region along a core-edge gradient in the eastern spadefoot toad (Pelobates syriacus)

Peripheral populations are those situated at the distribution margins of a species and are often subjected to more extreme abiotic and biotic conditions than those at the core. Here, we hypothesized that shorter repeat length and fewer heteroplasmic mitochondrial DNA (mtDNA) copies, which are associated with more efficient mitochondrial function, may be related to improved survival under extreme environmental conditions. We sampled eastern spadefoot toads (mostly as tadpoles) from 43 rain pools distributed along a 300-km gradient from core to edge of the species' distribution. We show that mean pool tandem repeat length and heteroplasmy increase from edge to core, even after controlling for body size. We evaluate several alternative hypotheses and propose the Fisher hypothesis as the most likely explanation. However, additional sequential sampling and experimental studies are required to determine whether selection under extreme conditions, or alternative mechanisms, could account for the gradient in heteroplasmy and repeat length in the mtDNA control region.     

Length and restriction site heteroplasmy in the mitochondrial DNA of american shad (alosa sapidissima)

Restriction endonuclease analysis was used to assess mitochondrial DNA (mtDNA) variation in American shad (Alosa sapidissima) collected from 14 rivers ranging from Florida to Quebec. Two types of heteroplasmy were observed, one involving a major length polymorphism and the other a single restriction site. Shad mtDNA occurred in two principal size classes, 18.3 and 19.8 kb. Of 244 shad examined, 30 were heteroplasmic and carried both size classes of mtDNA in varying proportions; the remainder were homoplasmic for the smaller size class of mtDNA. The large mtDNA variant occurred most frequently at the southern end of the range, and except for two individuals from Nova Scotia, was not detected among shad from rivers north of the Delaware. In contrast, ten shad heteroplasmic for a SalI restriction site originated from rivers ranging from South Carolina to Nova Scotia. DNA mapping and hybridization experiments indicated that the length polymorphism is in the D-loop-containing region and consists of a tandemly repeated 1.5-kb DNA sequence occurring in two and three copies, respectively, in the two major size classes of shad mtDNA. Continuous length variation up to approximately 40 bp occurs among copies of the repeat both within and among individuals. Restriction site data support the conclusion that both forms of heteroplasmy in shad mtDNA have originated more than once.     

Structural characterization of site-specific discontinuities associated with replication origins of minicircle DNA from Crithidia fasciculata

The kinetoplast DNA of trypanosomes is comprised of thousands of DNA minicircles and 20-50 maxicircles catenated into a single network. Replication intermediates of minicircle DNA from the trypanosomatid species Crithidia fasciculata contain site-specific discontinuities in both heavy (H) and light (L) strands. These discontinuities map to two small regions situated 180 degrees apart on the minicircle; each region has two sites at which a discontinuity can occur, one on each strand. We have determined the position of these discontinuities on the minicircle DNA sequence and have characterized their structure. H-strand discontinuities occur within a 4-5-nucleotide sequence and consist of single nicks, only one of which appears to be a DNA-DNA junction. Characterization of the remaining H-strand nicks indicates a structure other than a typical DNA-DNA or DNA-RNA junction. Discontinuities on the L-strand can be either a nick or a short gap which overlaps a 12-nucleotide sequence universally conserved among minicircles from various trypanosome species. Up to 6 nucleotides are hydrolyzed from the 5' terminus facing the gap upon treatment with alkali, suggesting the presence of an RNA primer. Based on the structures of minicircle replication intermediates, we present a model for replication of minicircle DNA in which the site-specific discontinuities closely coincide with the origins of replication.     

Complete mitochondrial genome of the miiuy croaker Miichthys miiuy (Perciformes, Sciaenidae) with phylogenetic consideration

The complete sequence of the 16,493 nucleotide mitochondrial genome from the single species of the family Sciaenidae, the miiuy croaker, Miichthys miiuy, was determined. The nucleotide sequences of M. miiuy mitochondrial DNA have been compared with those of three other Sciaenidae fishes. The contents of the M. miiuy mitochondrial genome are 13 protein-coding genes, two ribosomal RNA genes and 22 transfer RNA genes, and two non-coding regions (L-strand replication origin and control region), the gene order of which is identical to that observed in most vertebrates. The L-strand replication origin of M. miiuy is not pyrimidine-rich compared to those of most bony fishes. Within the control region, we identified the extended termination associated sequence domain, the central conserved sequence block domain and the conserved sequence block domain, while the typical central conserved blocks CSB-D, -E and -F could not be detected in the three other Sciaenidae species. In the ML phylogenetic analyses, the monophyly of Pseudosciaeniae was not supported, which is against with the morphological results. Collichthys niveatus is most closely related to Larimichthys polyactis, and Collichthys and Larimichthys may be merged into one genus, based on the current datasets.