Structure and evolution of the mitochondrial control region of the pollen beetle Meligethes thalassophilus (Coleoptera: Nitidulidae).

The organization of the mitochondrial DNA (mtDNA) control region (CR) of the pollen beetle Meligethes thalassophilus is described. This mtDNA CR represents the longest sequenced for beetles so far, since the entire nucleotide sequence ranges from approximately 5000 to approximately 5500 bp. The CR of M. thalassophilus is organized in three distinct domains: a conserved domain near the tRNAIle gene, a variable domain flanking the 12S rRNA gene, and a relatively large central tandem array made up of a variable number of approximately 170 bp repeats that is responsible for the intraspecific length variation observed. Like other CRs found in insects, the M. thalassophilus CR contains two long homopolymeric runs that may be involved in mtDNA replication. Furthermore, conserved stem-and-loop structures in the repetitive domain were identified and their possible role in generating length variation is examined. Intraspecific comparison of the tandem repeat elements of M. thalassophilus suggests mechanisms of concerted evolution leading to homogenization of the repetitive region. The utility of such an array of tandem repeats as a genetic marker for assessing population-level variability and evolutionary relationships among populations is discussed. Finally, the technical difficulties found in isolating the mtDNA CR in beetles are remarked upon.     

Evidence for recombination of mtDNA in the marine mussel Mytilus trossulus from the Baltic

A number of studies have claimed that recombination occurs in animal mtDNA, although this evidence is controversial. Ladoukakis and Zouros (2001) provided strong evidence for mtDNA recombination in the COIII gene in gonadal tissue in the marine mussel Mytilus galloprovincialis from the Black Sea. The recombinant molecules they reported had not however become established in the population from which experimental animals were sampled. In the present study, we provide further evidence of the generality of mtDNA recombination in Mytilus by reporting recombinant mtDNA molecules in a related mussel species, Mytilus trossulus, from the Baltic. The mtDNA region studied begins in the 16S rRNA gene and terminates in the cytochrome b gene and includes a major noncoding region that may be analogous to the D-loop region observed in other animals. Many bivalve species, including some Mytilus species, are unusual in that they have two mtDNA genomes, one of which is inherited maternally (F genome) the other inherited paternally (M genome). Two recombinant variants reported in the present study have population frequencies of 5% and 36% and appear to be mosaic for F-like and M-like sequences. However, both variants have the noncoding region from the M genome, and both are transmitted to sperm like the M genome. We speculate that acquisition of the noncoding region by the recombinant molecules has conferred a paternal role on mtDNA genomes that otherwise resemble the F genome in sequence.     

Hypervariability of the D-loop region in mitochondrial DNA of Russian sturgeon <Emphasis Type="Italic">Acipenser gueldenstaedtii</Emphasis> (Acipenseriformes,Acipenseridae)

The structure of the D-loop region in mitochondrial DNA (mtDNA) of Russian sturgeon Acipenser gueldenstaedtii from the Azov Sea population was studied with the method of direct sequencing. Interindividual heteroplasmy of the length of mtDNA in the region of D-loop realized by the presence of a different number of tandem repeats (82 pairs of bases) was found. Analysis of tandem repeats in the D-loop region in mtDNA in the studied sample (28 individuals) revealed eight mitotypes differed in the pattern of nucleotide substitution and in the number of tandem repeats (2, 3, and 4 repeats). Revealed mitotypes can be considered as potential genetic markers for different biological groups, schools, or seasonal races of A. gueldenstaedtii.     

Replication Intermediates of the Linear Mitochondrial DNA of Candida parapsilosis Suggest a Common Recombination Based Mechanism for Yeast Mitochondria

Variation in the topology of mitochondrial DNA (mtDNA) in eukaryotes evokes the question if differently structured DNAs are replicated by a common mechanism. RNA-primed DNA synthesis has been established as a mechanism for replicating the circular animal/mammalian mtDNA. In yeasts, circular mtDNA molecules were assumed to be templates for rolling circle DNA-replication. We recently showed that in Candida albicans, which has circular mapping mtDNA, recombination driven replication is a major mechanism for replicating a complex branched mtDNA network. Careful analyses of C. albicans-mtDNA did not reveal detectable amounts of circular DNA molecules. In the present study we addressed the question of how the unit sized linear mtDNA of Candida parapsilosis terminating at both ends with arrays of tandem repeats (mitochondrial telomeres) is replicated. Originally, we expected to find replication intermediates diagnostic of canonical bi-directional replication initiation at the centrally located bi-directional promoter region. However, we found that the linear mtDNA of Candida parapsilosis also employs recombination for replication initiation. The most striking findings were that the mitochondrial telomeres appear to be hot spots for recombination driven replication, and that stable RNA:DNA hybrids, with a potential role in mtDNA replication, are also present in the mtDNA preparations.     

Site-specific deletions of the mitochondrial genome in the Pearson marrow-pancreas syndrome.

The Pearson marrow-pancreas syndrome is a fatal disorder involving the hematopoietic system and the exocrine pancreas in early infancy. We have previously shown that this disease results from a widespread defect of oxidative phosphorylation. Here, we describe deletions of the mitochondrial (mt) genome between repeated 8- to 13-bp sequences as consistent features of the disease. Studying a series of nine unrelated children, including the patient originally reported by H. Pearson, we found five different types of direct repeats at the boundaries of the mtDNA deletions and we provided evidence for conservation of the 3'-repeated sequence in the deletions. In addition, we found a certain degree of homology between the nucleotide composition of the direct repeats and several structures normally involved in mtDNA replication and mtRNA processing. These results are consistent either with the recognition and cleavage of a particular DNA sequence with a factor of still unknown origin or with a homologous recombination between direct-repeat mtDNA sequences in the Pearson syndrome.     

Length and sequence variation in the mitochondrialDNA control region of the Etruscan freshwater goby Padogobius nigricans (Teleostei, Gobiidae)

The primary structures of the mtDNA control region, the flanking tRNA^Pro and tRNA^Phe genes and the partial 12S rRNA gene of the freshwater Etruscan goby Padogobius nigricans (Teleostei, Gobiidae) were determined. A tandem array of variable 17-bp repeats, named R1–7, was found to be responsible for most of the observed length and sequence variation between individuals.     

Sequence analysis of the complete mitochondrial DNA molecule of the hedgehog,Erinaceus europaeus,and the phylogenetic position of the Lipotyphla

The sequence of the mitochondrial DNA (mtDNA) molecule of the European hedgehog (Erinaceus europaeus) was determined. The length of the sequence presented is 17,442 nucleotides (nt). The molecule is thus the largest eutherian mtDNA molecule so far reported. The organization of the molecule conforms with that of other eutherians, but the control region of the molecule is exceptionally long, 1,988 nt, due to the presence of repeated motifs at two different positions in the 3 part of the control region. The length of the control region is not absolute due to pronounced heteroplasmy caused by variable numbers of the motif TACGCA in one of the repetitive regions. The sequence presented includes 46 repeats of this type. The other repeated region is composed of different AT-rich repeats. This region was identical among four clones studied. Comparison of mitochondrial peptide-coding genes identified a separate position of the hedgehog among several mammalian orders. The concatenated protein sequence of the 13 peptide-coding genes was used in a phylogenetic study using the opossum as outgroup. The position of the hedgehog sequence was basal among the other eutherian sequences included: human, rat, mouse, cow, blue whale, harbor seal, and horse. The analysis did not resolve the relationship among carnivores, perissodactyls, and artiodactyls/cetaceans, suggesting a closer relationship among these orders than acknowledged by classical approaches. Correspondence to: U. Arnason     

Recombination sequences in plant mitochondrial genomes: diversity and homologies to known mitochondrial genes.

Several plant mitochondrial genomes contain repeated sequences that are postulated to be sites of homologous intragenomic recombination (1-3). In this report, we have used filter hybridizations to investigate sequence relationships between the cloned mitochondrial DNA (mtDNA) recombination repeats from turnip, spinach and maize and total mtDNA isolated from thirteen species of angiosperms. We find that strong sequence homologies exist between the spinach and turnip recombination repeats and essentially all other mitochondrial genomes tested, whereas a major maize recombination repeat does not hybridize to any other mtDNA. The sequences homologous to the turnip repeat do not appear to function in recombination in any other genome, whereas the spinach repeat hybridizes to reiterated sequences within the mitochondrial genomes of wheat and two species of pokeweed that do appear to be sites of recombination. Thus, although intragenomic recombination is a widespread phenomenon in plant mitochondria, it appears that different sequences either serve as substrates for this function in different species, or else surround a relatively short common recombination site which does not cross-hybridize under our experimental conditions. Identified gene sequences from maize mtDNA were used in heterologous hybridizations to show that the repeated sequences implicated in recombination in turnip and spinach/pokeweed/wheat mitochondria include, or are closely linked to genes for subunit II of cytochrome c oxidase and 26S rRNA, respectively. Together with previous studies indicating that the 18S rRNA gene in wheat mtDNA is contained within a recombination repeat (3), these results imply an unexpectedly frequent association between recombination repeats and plant mitochondrial genes.     

Inter-species sequence diversity in the replication initiation region of Paramecium mitochondrial DNA

Mitochondrial DNA from Paramecium aurelia is a linear molecule. Replication is initiated at a unique cross-linked molecular terminus. During replication dimer length molecules, consisting of two head-to-head monomers, are generated. We have cloned the head-to-head dimer initiation region from five different species and several stocks (or races) within species and determined its DNA sequence. For all species, this dimer initiation region consists of a central non-palindromic sequence containing almost exclusively A and T, arranged in an array of direct tandem repeats. In an intra-species comparison, the sequences of the repeat units are relatively homogeneous; inter-species comparisons, however, show diversity except for a conserved "Goldberg-Hogness box", T-A-T-A-A-A-T-A. The size of a repeat unit and the number of repeats within a molecule can vary over a wide range, even in an intra-species comparison. Because of these wide inter-species variations observed, it is likely that the function of this region imposes few constraints on the sequence other than its high A + T content and possibly a Goldberg-Hogness box. The array of direct tandem repeats may have arisen from unequal recombination or crossover within this region. Adjacent to the non-palindromic region is a transcribed sequence which is highly conserved for all species and presumably represents a gene coding region.     

Generation of VNTRs and heteroplasmy by sequence turnover in the mitochondrial control region of two elephant seal species

We describe an unusual repetitive DNA region located in the 3′ end of the light (L)-strand in the mitochondrial control region of two elephant seal species. The array of tandem repeats shows both VNTR (variable-number tandem repeat) and sequence variation and is absent from 12 compared mammalian species, except for the occurrence in the same location of a distinct repetitive region in rabbit mtDNA and a similar repeat in the harbor seal. The sequence composition and arrangement of the repeats differ considerably between the northern elephant seal (Mirounga angustirostris) and the southern species (M. leonina) despite an estimated divergence time of 1 MY (based on an mtDNA-RNA gene and the nonrepetitive control region). Analysis of repeat sequence relationships within and between species indicate that divergence in sequence and structure of repeats has involved both slippagelike and unequal crossingover processes of turnover, generating very high levels of divergence and heteroplasmy. Presented at the NATO Advanced Research Workshop onGenome Organization and Evolution, Spetsai, Greece, 16–22 September 1992     

Origin and Evolution of Tandem Repeats in the Mitochondrial DNA Control Region of Shrikes (<Emphasis Type="BoldItalic">Lanius</Emphasis> spp.)

The origin and evolution of a 128-bp tandem repeat in the mtDNA control region of shrikes (Lanius: Aves) were investigated. The tandem repeat is present in only two species, L. excubitor and L. ludovicianus. In contrast to the variation in repeat number in L. ludovicianus, all individuals of three subspecies of L. excubitor had three repeats. Comparative analysis suggests that a short direct repeat, and a secondary structure including the tandem repeat and a downstream inverted repeat, may be important in the origin of the tandem repeat by slipped-strand mispairing and its subsequent turnover. Homogenization of repeat sequences is most simply explained by expansion and contraction of the repeat array. Surprisingly, mtDNA sequences from L. excubitor were found to be paraphyletic with respect to L. ludovicianus. These results show the utility of a comparative analysis for insights into the evolutionary dynamics of mtDNA tandem repeats.[Reviewing Editor: Martin Kreitman]     

Nonneutral evolution of tandem repeats in the mitochondrial DNA control region of lagomorphs

The mitochondrial DNA of the European rabbit (Oryctolagus cuniculus) contains a tandem array of 153-bp repeats in the vicinity of the replication origin of the H-stand. Variation among molecules in the number of these repeats results in inter- and intraindividual length polymorphism (heteroplasmy). Generally, in an individual, one predominant molecular type is observed, the others representing a low percentage of the mtDNA content. At the tissue level, we observe a particular distribution of this polymorphism in the gonads compared with liver, kidneys, or brain, implying a relationship between the differentiation status of the cells and the types of new mtDNA molecules which appear and accumulate during lifetime. Similar tandem repeats were also found in the mtDNA noncoding region of European hares (Lepus europaeus), a cottontail (Sylvilagus floridanus), and a pika (Ochotona rufescens). The lengths and the sequences of these units evolve rapidly and in a concerted way, but the number of repeats is maintained in a narrow range, and an internal 20-bp segment is highly conserved. Constraints restrict the evolution of the primary sequence of these repeated units, the number of which is probably controlled by a stabilizing selection.      

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.     

One-way sequencing of multiple amplicons from tandem repetitive mitochondrial DNA control region

Repetitive DNA sequences not only exist abundantly in eukaryotic nuclear genomes, but also occur as tandem repeats in many animal mitochondrial DNA (mtDNA) control regions. Due to concerted evolution, these repetitive sequences are highly similar or even identical within a genome. When long repetitive regions are the targets of amplification for the purpose of sequencing, multiple amplicons may result if one primer has to be located inside the repeats. Here, we show that, without separating these amplicons by gel purification or cloning, directly sequencing the mitochondrial repeats with the primer outside repetitive region is feasible and efficient. We exemplify it by sequencing the mtDNA control region of the mosquito Aedes albopictus, which harbors typical large tandem DNA repeats. This one-way sequencing strategy is optimal for population surveys.     

Unique features in the mitochondrial D-loop region of the European seabass Dicentrarchus labrax

We have cloned and sequenced the displacement-loop (D-loop) region of the mitochondrial DNA (mtDNA) from the European seabass Dicentrarchus labrax (Dl). This sequencing revealed the presence of four tandemly repeated elements (R1, R2, R3 and R4); the individual variation in mtDNA total length is entirely accounted for by their variable number. The individuals examined also possessed an imperfect copy of one of the tandem repeats (ΨR2). At least one termination-associated sequence (TAS) is present in each of the repeats and in two copies 5′ upstream from the tandem array as well. The alignment of the Dl D-loop region with D-loop sequences from four other Teleosts and one Chondrosteus showed the Dl sequence to be larger than that of other fish. The extraordinary length of the D1 D-loop sequence is also due to the 5′ and 3′ regions that are flanking the tandem array, the largest ones to date analyzed in fish. In this study, we also report the unique organization and localization of putative TAS and conserved-sequence block (CSB) elements, and the presence of a conserved 218-bp sequence in the D1 D-loop region.     

The Agrocybe aegerita mitochondrial genome contains two inverted repeats of the nad4 gene arisen by duplication on both sides of a linear plasmid integration site

The Agrocybe aegerita mitochondrial genome possesses two polB genes with linear plasmid origin. The cloning and sequencing of the regions flanking Aa-polB P1 revealed two large inverted repeats (higher than 2421 nt) separated by a single copy region of 5834 nt. Both repeats contain identical copies of the nad4 gene. The single copy region contains two disrupted genes with plasmid origin Aa-polB P1 and a small ORF homologous to a small gene described in two basidiomycete linear plasmids. The phylogenetic analyses argue in favor of a same plasmid origin for both genes but, surprisingly, these genes were separated by a mitochondrial tRNA-Met. Both strands of the complete region containing the two nad4 inverted copies and the tRNA-Met appear to be transcribed on large polycistronic mRNAs. A model summarizing the events that would have occurred is proposed: (1) capture of the tRNA by the plasmid before its integration in the mtDNA or acquisition of the tRNA gene by recombination after the plasmid integration, (2) integration of the plasmid in the mtDNA, accompanied by a large duplication containing the nad4 gene and (3) erosion of the plasmid sequences by large deletions and mutations.     

The Agrocybe aegerita mitochondrial genome contains two inverted repeats of the nad4 gene arisen by duplication on both sides of a linear plasmid integration site.

The Agrocybe aegerita mitochondrial genome possesses two polB genes with linear plasmid origin. The cloning and sequencing of the regions flanking Aa-polB P1 revealed two large inverted repeats (higher than 2421 nt) separated by a single copy region of 5834 nt. Both repeats contain identical copies of the nad4 gene. The single copy region contains two disrupted genes with plasmid origin Aa-polB P1 and a small ORF homologous to a small gene described in two basidiomycete linear plasmids. The phylogenetic analyses argue in favor of a same plasmid origin for both genes but, surprisingly, these genes were separated by a mitochondrial tRNA-Met. Both strands of the complete region containing the two nad4 inverted copies and the tRNA-Met appear to be transcribed on large polycistronic mRNAs. A model summarizing the events that would have occurred is proposed: (1) capture of the tRNA by the plasmid before its integration in the mtDNA or acquisition of the tRNA gene by recombination after the plasmid integration, (2) integration of the plasmid in the mtDNA, accompanied by a large duplication containing the nad4 gene and (3) erosion of the plasmid sequences by large deletions and mutations.     

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.      

Evidence of animal mtDNA recombination between divergent populations of the potato cyst nematode Globodera pallida

Recombination is typically assumed to be absent in animal mitochondrial genomes (mtDNA). However, the maternal mode of inheritance means that recombinant products are indistinguishable from their progenitor molecules. The majority of studies of mtDNA recombination assess past recombination events, where patterns of recombination are inferred by comparing the mtDNA of different individuals. Few studies assess contemporary mtDNA recombination, where recombinant molecules are observed as direct mosaics of known progenitor molecules. Here we use the potato cyst nematode, Globodera pallida, to investigate past and contemporary recombination. Past recombination was assessed within and between populations of G. pallida, and contemporary recombination was assessed in the progeny of experimental crosses of these populations. Breeding of genetically divergent organisms may cause paternal mtDNA leakage, resulting in heteroplasmy and facilitating the detection of recombination. To assess contemporary recombination we looked for evidence of recombination between the mtDNA of the parental populations within the mtDNA of progeny. Past recombination was detected between a South American population and several UK populations of G. pallida, as well as between two South American populations. This suggests that these populations may have interbred, paternal mtDNA leakage occurred, and the mtDNA of these populations subsequently recombined. This evidence challenges two dogmas of animal mtDNA evolution; no recombination and maternal inheritance. No contemporary recombination between the parental populations was detected in the progeny of the experimental crosses. This supports current arguments that mtDNA recombination events are rare. More sensitive detection methods may be required to adequately assess contemporary mtDNA recombination in animals.     

RepeatAround: a software tool for finding and visualizing repeats in circular genomes and its application to a human mtDNA database

RepeatAround is a Windows based software tool designed to find "direct repeats", "inverted repeats", "mirror repeats" and "complementary repeats", from 3 to 64bp length, in circular genomes. It processes input files directly extracted from GenBank database, providing visualisation of the repeats location in the genomic structure, so that for instance, in most mtDNAs the user can check if the repeats are located in coding or non-coding region (and in the first case in which gene), and how far apart the repeat pair(s) are. Besides the visual tool, it provides other outputs in a spreadsheet containing information on the number and location of the repeats, facilitating graphic analyses. Several genomes can be inputed simultaneously, for phylogenetic comparison purposes. Other capabilities of the software are the generation of random circular genomes, for statistical evaluation of comparison between observed repeats distributions with their shuffled counterparts, as well as the search for specific motifs, allowing an easy confirmation of repeats flanking a newly detected rearrangement. As an example of the programme's applications we analysed the Direct Repeats distribution in a large human mtDNA database. Results showed that Direct Repeats, even the larger ones, are evenly distributed among the human mtDNA haplogroups, enabling us to state that, based only on the repetitive motifs, no haplogroup is particularly more or less prone to mtDNA macrodeletions.