Evolution by DNA turnover in the control region of vertebrate mitochondrial DNA

Species-specific variation in mitochondrial genome tandem repeat polymorphisms in hares (Lepus spp., Lagomorpha, Leporidae) provides insight into their evolution

The non-coding regions of the mitochondrial DNAs (mtDNAs) of hares, rabbits, and pikas (Lagomorpha) contain short (∼20 bp) and long (130-160 bp) tandem repeats, absent in related mammalian orders. In the presented study, we provide in-depth analysis for mountain hare (Lepus timidus) and brown hare (L. europaeus) mtDNA non-coding regions, together with a species- and population-level analysis of tandem repeat variation. Mountain hare short tandem repeats (SRs) as well as other analyzed hare species consist of two conserved 10 bp motifs, with only brown hares exhibiting a single, more variable motif. Long tandem repeats (LRs) also differ in sequence and copy number between species. Mountain hares have four to seven LRs, median value five, while brown hares exhibit five to nine LRs, median value six. Interestingly, introgressed mountain hare mtDNA in brown hares obtained an intermediate LR length distribution, with median copy number being the same as with conspecific brown hare mtDNA. In contrast, transfer of brown hare mtDNA into cultured mtDNA-less mountain hare cells maintained the original LR number, whereas the reciprocal transfer caused copy number instability, suggesting that cellular environment rather than the nuclear genomic background plays a role in the LR maintenance. Due to their dynamic nature and separation from other known conserved sequence elements on the non-coding region of hare mitochondrial genomes, the tandem repeat elements likely to represent signatures of ancient genetic rearrangements. clarifying the nature and dynamics of these rearrangements may shed light on the possible role of NCR repeated elements in mitochondria and in species evolution.

Selected Lark Mitochondrial Genomes Provide Insights into the Evolution of Second Control Region with Tandem Repeats in Alaudidae (Aves, Passeriformes)

The control region (CR) regulates the replication and transcription of the mitochondrial genome (mitogenome). Some avian mitogenomes possess two CRs, and the second control region (CR2) may enhance replication and transcription; however, the CR2 in lark mitogenome appears to be undergoing loss and is accompanied by tandem repeats. Here, we characterized six lark mitogenomes from Alaudala cheleensis, Eremophila alpestris, Alauda razae, and Calandrella cinerea and reconstructed the phylogeny of Passerida. Through further comparative analysis among larks, we traced the evolutionary process of CR2. The mitochondrial gene orders were conserved in all published lark mitogenomes, with Cytb-trnT-CR1-trnP-ND6-trnE-remnant CR2 with tandem repeat-trnF-rrnS. Phylogenetic analysis revealed Alaudidae and Panuridae are sister groups at the base of Sylvioidea, and sporadic losses of CR2 may occur in their common ancestor. CR sequence and phylogeny analysis indicated CR2 tandem repeats were generated within CR2, originating in the ancestor of all larks, rather than inherited from CR1. The secondary structure comparison of tandem repeat units within and between species suggested slipped-strand mispairing and DNA turnover as suitable models for explaining the origin and evolution of these repeats. This study reveals the evolutionary process of the CR2 containing tandem repeat in Alaudidae, providing reference for understanding the evolutionary characteristics and dynamics of tandem repeats.

Extreme variation in patterns of tandem repeats in mitochondrial control region of yellow-browed tits (Sylviparus modestus, Paridae)

To investigate the evolutionary pattern and origins of tandem repeats in the mitochondrial control region of the yellow-browed tit (Sylviparus modestus), the control region and another four mitochondrial loci from fifteen individuals were analyzed. A 117-bp tandem repeat unit that repeated once, twice or three times in different individuals was found, and a rarely reported arrangement for this tandem repeats region that a 5' imperfect copy at its downstream and a 3' imperfect copy at its upstream was observed. The haplotype network, phylogenetic trees, and ancestral state reconstruction of the combined dataset of five loci suggested multiple origins of the same repeat number. The turnover model via slipped-strand mispairing was introduced to interpret the results, because mispairing occurred so frequently that multiple origins of certain repeat number were observed. Insertion via recombination should be a better explanation for the origin of this tandem repeat unit, considering characteristics of the combined sequence of the 3' and 5' imperfect copy, including identification of its homolog in other passerines and its predicted secondary structure.

An examination of the origin and evolution of additional tandem repeats in the mitochondrial DNA control region of Japanese sika deer (Cervus Nippon)

Tandem repeat units are only detected in the left domain of the mitochondrial DNA control region in sika deer. Previous studies showed that Japanese sika deer have more tandem repeat units than its cousins from the Asian continent and Taiwan, which often have only three repeat units. To determine the origin and evolution of these additional repeat units in Japanese sika deer, we obtained the sequence of repeat units from an expanded dataset of the control region from all sika deer lineages. The functional constraint is inferred to act on the first repeat unit because this repeat has the least sequence divergence in comparison to the other units. Based on slipped-strand mispairing mechanisms, the illegitimate elongation model could account for the addition or deletion of these additional repeat units in the Japanese sika deer population. We also report that these additional repeat units could be occurring in the internal positions of tandem repeat regions, possibly via coupling with a homogenization mechanism within and among these lineages. Moreover, the increased number of repeat units in the Japanese sika deer population could reflect a balance between mutation and selection, as well as genetic drift.

Limited phylogenetic distribution of a long tandem-repeat cluster in the mitochondrial control region in Bubo (Aves, Strigidae) and cluster variation in Blakiston's fish owl (Bubo blakistoni)

To investigate the phylogenetic position of Blakiston's fish owl (Bubo blakistoni), we sequenced the mitochondrial (mt) DNA control region and cytochrome b (cyt b) for nine Bubo species. Maximum-likelihood analyses of combined control region and cyt b sequences, and cyt b sequences alone, showed that species formerly placed in genus Ketupa comprise a monophyletic group. Unexpectedly, we discovered a long cluster of 20-25 tandem repeat units 77 or 78bp long in the third control region domain in four of the nine Bubo species for which the control region was sequenced (B. blakistoni, B. flavipes, and B. ketupu in the Ketupa clade; B. lacteus), leading to overall control region lengths of 3.0-3.8kpb estimated from agarose gel electrophoresis. The control region in B. lacteus is the longest (3.8kbp) reported to date in vertebrates. Sequencing of eight repeat units at each end of the cluster in 20 B. blakistoni individuals detected several types of repeat units 77 or 78bp long, and six patterns in the order of unit types. The occurrence of a repeat cluster in all three species examined in the Ketupa clade suggests their common ancestor also had a cluster, whereas a maximum parsimony tree showed repeat-unit types grouping by species, rather than by paralog groups, suggesting independent origins of the clusters. We reconcile these results with a turnover model, in which the range in cluster-length variation and unit types at the 5' end are hypothetically functionally constrained by the protein-binding function of the control region, but otherwise there is a continual turnover of units in evolutionary time, with new unit types arising through mutations, proliferating by duplication of single and double repeat blocks, and being lost through deletion. Estimated free energies for reconstructed secondary structures of single and especially pairs of repeat units were higher than for homologous single-unit blocks in species lacking a repeat cluster, supporting slipped-strand mispairing as the mechanism of cluster turnover.

Inferring recent historic abundance from current genetic diversity

Recent historic abundance is an elusive parameter of great importance for conserving endangered species and understanding the pre-anthropogenic state of the biosphere. The number of studies that have used population genetic theory to estimate recent historic abundance from contemporary levels of genetic diversity has grown rapidly over the last two decades. Such assessments often yield unexpectedly large estimates of historic abundance. We review the underlying theory and common practices of estimating recent historic abundance from contemporary genetic diversity, and critically evaluate the potential issues at various estimation steps. A general issue of mismatched spatio-temporal scales between the estimation itself and the objective of the estimation emerged from our assessment; genetic diversity-based estimates of recent historic abundance represent long-term averages, whereas the objective typically is an estimate of recent abundance for a specific population. Currently, the most promising approach to estimate the difference between recent historic and contemporary abundance requires that genetic data be collected from samples of similar spatial and temporal duration. Novel genome-enabled inference methods may be able to utilize additional information of dense genome-wide distributions of markers, such as of identity-by-descent tracts, to infer recent historic abundance from contemporary samples only.

Mitochondrial DNA STR analysis as a tool for studying the green sea turtle (Chelonia mydas) populations: the Mediterranean Sea case study

The Mediterranean population of the green sea turtle Chelonia mydas is critically endangered. Genetic analysis of this population using the ordinary haplotyping system, based on sequence analysis of a segment of the mitochondrial DNA (mtDNA) D-loop (control region), revealed very little variation. The most common haplotype, CM-A13, was observed in all but three individuals in hundreds of samples in previous studies. In search for a more informative marker we sequenced the 3' of the mitochondrial control region which contains an AT-rich microsatellite. We found a unique pattern that consists of four AT short tandem repeats (STRs) with varying copy numbers. This allowed us to construct a new haplotyping system composed of four different STR sizes for each mtDNA sequence. Our new mitochondrial STR (mtSTR) haplotyping approach revealed 33 different haplotypes within the nesting and stranded sea turtles along the Mediterranean Israeli seashore. The Israeli coast nesting females had 10 different haplotypes that can be used for monitoring and conservation purposes. The mtSTR haplotyping system can clearly assist in fingerprinting of individual turtles. Moreover, it can be used for estimating phylogenetic distances within populations. This case study shows that the mtSTR haplotyping is applicable for the study of global green sea turtle populations and could also be considered as markers of genetic variability in other species.

Radiation and speciation of pelagic organisms during periods of global warming: the case of the common minke whale, Balaenoptera acutorostrata

How do populations of highly mobile species inhabiting open environments become reproductively isolated and evolve into new species? We test the hypothesis that elevated ocean-surface temperatures can facilitate allopatry among pelagic populations and thus promote speciation. Oceanographic modelling has shown that increasing surface temperatures cause localization and reduction of upwelling, leading to fragmentation of feeding areas critical to pelagic species. We test our hypothesis by genetic analyses of populations of two closely related baleen whales, the Antarctic minke whale (Balaenoptera bonaerensis) and common minke whale (Balaenoptera acutorostrata) whose current distributions and migration patterns extent are largely determined by areas of consistent upwelling with high primary production. Phylogeographic and population genetic analyses of mitochondrial DNA control-region nucleotide sequences collected from 467 whales sampled in four different ocean basins were employed to infer the evolutionary relationship among populations of B. acutorostrata by rooting an intraspecific phylogeny with a population of B. bonaerensis. Our findings suggest that the two species diverged in the Southern Hemisphere less than 5 million years ago (Ma). This estimate places the speciation event during a period of extended global warming in the Pliocene. We propose that elevated ocean temperatures in the period facilitated allopatric speciation by disrupting the continuous belt of upwelling maintained by the Antarctic Circumpolar Current. Our analyses revealed that the current populations of B. acutorostrata likely diverged after the Pliocene some 1.5 Ma when global temperatures had decreased and presumably coinciding with the re-establishment of the polar-equatorial temperature gradient that ultimately drives upwelling. In most population samples, we detected genetic signatures of exponential population expansions, consistent with the notion of increasing carrying capacity after the Pliocene. Our hypothesis that prolonged periods of global warming facilitate speciation in pelagic marine species that depend on upwelling should be tested by comparative analyses in other pelagic species.

A novel mitochondrial intergenic spacer reflecting population structure of Pacific oyster

Nucleotide sequence divergence in a novel major mitochondrial DNA intergenic spacer (IGS) of Pacific oyster Crassostrea gigas was analyzed for 29 cultured individuals within the Goseong population (Korea). A total of 7 variable sites were detected within the IGS, and the relative frequency of nucleotide alteration was determined to be 1.16%;. All alterations were due to a single nucleotide substitution, and 5 transitions and 2 transversions were observed. Among 29 specimens, only 8 haplotypes could be identified, and 6 of the haplotypes were unique to particular specimens. Pairwise genetic diversity of all 8 haplotypes was calculated to be 0.412+/-0.134 from multiple sequence substitutions based on the two-parameter model. The phylogenetic tree obtained for these haplotypes according to the neighbor-joining method illustrated a single cluster of linkages, which comprised 5 haplotypes associated with 23 specimens, while the other 3 haplotypes associated with 6 specimens were scattered. The results indicate that the IGS is higher polymorphic and thus more suitable as a genetic marker for population structure analysis of Pacific oyster than the mtDNA coding regions, such as cytochrome c oxidase I and 16S ribosomal RNA genes.

Independent occurrences of multiple repeats in the control region of mitochondrial DNA of white-tailed deer

Deer in the genera Mazama and Odocoileus generally have two copies of a 75-base-pair (bp) repeat in the left domain of the control region of the mitochondrial DNA (mtDNA). Phylogenetic analyses further suggest an ancient origin for the duplication supporting a previously stated contention that this event occurred before the separation of Mazama and Odocoileus. However, white-tailed deer (Odocoileus virginianus) had three or four copies of a 75-bp repeat in the control region of their mtDNA in 7.8% of the individuals analyzed, and all of these animals were from the coastal plain of the southeastern United States. When copy 3 is present, it is very similar in sequence to copy 2, but variation suggests that copy 3 probably evolved multiple times from copy 2. The pattern of phylogenetic clustering of the haplotypes from across the coastal plain also suggests that phenotypes with three or four copies of the repeat have originated multiple times. The 44 observed haplotypes showed strong spatial subdivision across the area with subpopulations frequently showing complete shifts in haplotype frequencies from others taken from nearby areas. Many of the subpopulations right along the coast or on adjacent barrier islands have a limited number of haplotypes as would occur in populations undergoing drift because of small numbers of breeding females and limited female dispersal.

Origin and evolution of tandem repeats in the mitochondrial DNA control region of shrikes (Lanius 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.

Sequence variation in the mitochondrial DNA control region of wild African cheetahs (Acinonyx jubatus)

Five hundred and twenty-five bp of mitochondrial control region were sequenced and analysed for 20 Acinonyx jubatus and one Felis catus. These sequences were compared with published sequences from another domestic cat, 20 ocelots (Leopardus pardalus) and 11 margays (Leopardus weidii). The intraspecific population divergence in cheetahs was found to be less than in the other cats. However variation was present and distinct groups of cheetahs were discernible. The 80 bp RS2 repetitive sequence motif previously described in other felids was found in four copies in cheetah. The repeat units probably have the ability to form secondary structure and may have some function in the regulation of control region replication. The two central repeat units in cheetah show homogenization that may have arisen by convergent evolution.

Do mitochondria recombine in humans?

Until very recently, mitochondria were thought to be clonally inherited through the maternal line in most higher animals. However, three papers published in 2000 claimed population-genetic evidence of recombination in human mitochondrial DNA. Here I review the current state of the debate. I review the evidence for the two main pathways by which recombination might occur: through paternal leakage and via a mitochondrial DNA sequence in the nuclear genome. There is no strong evidence for either pathway, although paternal leakage seems a definite possibility. However, the population-genetic evidence, although not conclusive, is strongly suggestive of recombination in mitochondrial DNA. The implications of non-clonality for our understanding of human and mitochondrial evolution are discussed.

Mammalian mitochondrial D-loop region structural analysis: identification of new conserved sequences and their functional and evolutionary implications

This paper reports the first comprehensive analysis of Displacement loop (D-loop) region sequences from ten different mammalian orders. It represents a systematic evolutionary study at the molecular level on regulatory homologous regions in organisms belonging to a well defined class, mammalia, which radiated about 150 million years ago (Mya). We have aligned and analyzed 26 complete D-loop region sequences available in the literature and the fat dormouse sequence, recently determined in our laboratory. The novelty of our alignment consists of the extensive manual revision of the preliminary output obtained by computer program to optimize sequence similarity, particularly for the two peripheral domains displaying heterogeneity in length and the presence of repeated sequences. The multialignment is available at the WWW site: http://www.ba.cnr.it/dloop.html. Our comparative study has allowed us to identify new conserved sequence blocks present in all the species under consideration and events of insertion/deletion which have important implications in both functional and evolutionary aspects. In particular we have detected two blocks, about 60 bp long, extended termination associated sequences (ETAS1 and ETAS2) conserved in all the organisms considered. Evaluation against experimental work suggests a possible functional role of ETAS1 and ETAS2 in the regulation of replication and transcription and targeted experimental approaches. The analyses on conserved sequence blocks (CSBs) clearly indicate that CSB1 is the only very essential element, common to all mammalian mt genomes, while CSB2 and CSB3 could be involved in different though related functions, probably species specific, and thus more linked to nuclear mitochondrial coevolutionary processes. Our hypothesis on the different functional implications of the conserved elements, CSBs and TASs, reported so far as main regulatory signals, would explain the different conservation of these elements in evolution. Moreover the intra-order comparison of the D-loop regions highlights peculiar features useful to define the evolutionary dynamics of this region in closely related species.

SSCP analysis of pig mitochondrial DNA D-loop region polymorphism

The sequence polymorphism that occurs in the mitochondrial DNA (mtDNA) displacement (D)-loop region is useful as a cytoplasmic DNA marker. We cloned the mtDNA D-loop regions of five breeds of pig by polymerase chain reaction (PCR) and determined their sequences. The sequence diversities in D-loop regions among five breeds of pig were located in the starting area of heavy-strand replication. From these sequences, we designed primers for PCR-mediated single-strand conformation polymorphism (PCR-SSCP) analysis that amplified the most polymorphic 227 bp fragment of the D-loop region. The results of PCR-SSCP analysis clearly showed that four types of polymorphism (A to D) are found in Landrace (A), Large White (A, B), Duroc (A), Göttingen miniature pig (B) and Meishan (C, D). The same polymorphisms were also detected from each porcine embryo by this method. Our results show that PCR-SSCP analysis is useful in detecting polymorphisms in the D-loop region of pigs and pig embryos.

Polymorphic sequence in the D-loop region of equine mitochondrial DNA

The D-loop regions in equine mitochondrial DNA were cloned from three thoroughbred horses by polymerase chain reaction (PCR). The total number of bases in the D-loop region were 1114 bp, 1115 bp and 1146 bp. The equine D-loop region is A/T rich like many other mammalian D-loops. The large central conserved sequence block and small conserved sequence blocks 1, 2 and 3, that are common to other mammals, were observed. Between conserved sequence blocks 1 and 2 there were tandem repeats of an 8 bp equine-specific sequence TGTGCACC, and the number of tandem repeats differed among individual horses. The base composition in the unit of these repeats is G/C rich as are the short repeats in the D-loops of rabbit and pig. Comparing DNA sequences between horse and other mammals, the difference in the D-loop region length is mostly due to the difference in the number of DNA sequences at both extremities. The similarities of the DNA sequences are in the middle part of the D-loop. In comparison of the sequences among three thoroughbred horses, it was determined that the region between tRNA(Pro) and the large central conserved sequence block was the richest in variation. PCR primers in the D-loop region were designed and the expected maternal inheritance was confirmed by PCR-RFLP (restriction fragment length polymorphism).