Evolution of duplicate control regions in the mitochondrial genomes of metazoa: a case study with Australasian Ixodes ticks

To investigate the evolution pattern and phylogenetic utility of duplicate control regions (CRs) in mitochondrial (mt) genomes, we sequenced the entire mt genomes of three Ixodes species and part of the mt genomes of another 11 species. All the species from the Australasian lineage have duplicate CRs, whereas the other species have one CR. Sequence analyses indicate that the two CRs of the Australasian Ixodes ticks have evolved in concert in each species. In addition to the Australasian Ixodes ticks, species from seven other lineages of metazoa also have mt genomes with duplicate CRs. Accumulated mtDNA sequence data from these metazoans and two recent experiments on replication of mt genomes in human cell lines with duplicate CRs allowed us to re-examine four intriguing questions about the presence of duplicate CRs in the mt genomes of metazoa: (1) Why do some mt genomes, but not others, have duplicate CRs? (2) How did mt genomes with duplicate CRs evolve? (3) How could the nucleotide sequences of duplicate CRs remain identical or very similar over evolutionary time? (4) Are duplicate CRs phylogenetic markers? It appears that mt genomes with duplicate CRs have a selective advantage in replication over mt genomes with one CR. Tandem duplication followed by deletion of genes is the most plausible mechanism for the generation of mt genomes with duplicate CRs. Once duplicate CRs occur in an mt genome, they tend to evolve in concert, probably by gene conversion. However, there are lineages where gene conversion may not always occur, and, thus, the two CRs may evolve independently in these lineages. Duplicate CRs have much potential as phylogenetic markers at low taxonomic levels, such as within genera, within families, or among families, but not at high taxonomic levels, such as among orders.     

Comparison research and phylogenetic implications of mitochondrial control regions in four soft-shelled turtles of Trionychia (Reptilia, Testudinata)

The mitochondrial control regions (CRs) and flanking sequences of Pelodiscus sinensis, Apalone ferox, Palea steindachneri and Carettochelys insculpta were obtained using Long-PCR with gene-specific primers. The CR lengths of the four species were 1843 bp, 1356 bp, 1725 bp, and 969 bp. The base composition percentages of A+T were 60.5%, 60.7%, 65.7%, 64.7%, respectively. Combined with CR sequences of other three soft-shelled turtles published in GenBank (Pelodiscus sinensis, Korea, AY962573; Dogania subplana, AF366350; Lissemys punctata, EF050073), we compared the CR structures and identified three functional domains (TAS, CD and CSB) in which conserved sequence blocks (TAS, CSB -F, CSB-1, CSB-2 and CSB-3) were also successfully identified according to their sequence similarities to those of other turtles. The variable numbers of tandem repeats (VNTRs 1) with 50–52 bp motif were identified at 5′-end of CR among the five soft-shelled turtles P. sinensis (China), P. sinensis (Korea), A. ferox, P. steindachneri, D. subplana. The copy number of the VNTRs varied from 5 to 15. VNTRs 2 with 2–11 bp motif were identified in the 3′- end of CR among all of the six soft-shelled turtles with variable number of motifs from 4 to 29. Moreover, VNTRs 3 with 6 bp motif were identified between CSB-1 and CSB-2 of CR both in P. sinensis (China) and P. sinensis (Korea), in which the number of motifs varied from 19 to 29. The types and distribution of VNTRs of the six soft-shelled turtles were also discussed. With Alligator mississippiensis as an outgroup, combined with the CR sequences (excluding VNTRs) of other five turtles which were published in GenBank, the molecular phylogenetic trees were constructed using PAUP 4.0b10 and MrBayes ver. 3.0. The results strongly supported the monophyly of Carretochelyidae and Carettochelyidae as sister group to an assemblage of Cryptodira. Our research suggested that the earliest phylogenetic tree splits into three separated basal branches; the Pelomedusidira (Pelomedusa subrufa), the Carettochelyidae (C. insculpta), and an assemblage of Cryptodira and the C. insculpta that might be a representation of distinctive suborder.     

Comparison of complete mitochondrial DNA control regions among five Asian freshwater turtle species and their phylogenetic relationships

The complete mitochondrial DNA (mtDNA) control regions (CR), cytochrome b (Cyt b), NADH dehydrogenase subunit 4 (ND4) and cytochrome coxidase subunit I (CO I) genes of four Asian freshwater turtles, Mauremys japonica, Ocadia sinensis, M. mutica, and Annamemys annamensis, were sequenced using universal PCR and long-PCR techniques. Combined with CR sequences of Chinemys reevesii, the composition and structure of CR of the five species were compared and analyzed. Three functional domains (TAS, CD and CSB) in CR and their conserved sequences (TAS, CSB-F, CSB-1, CSB-2, and CSB-3) were identified based on sequence similarity to those of other turtles. At the 3' end of CSB, six type motifs of variable number of tandem repeats (VNTRs) of five species were recognized, in which the TTATATTA motif may be the VNTR motif of the ancestral species of these five turtles. Comparison of nucleotide divergences among Cyt b, ND4, CO I, and CR of 11 turtle species using transitions + transversions and transversions-only methods supported the conclusion that CR evolved 2.6- to 5.7-fold faster than the other mtDNA genes. After excluding VNTRs of CR, molecular phylogenetic trees were constructed with maximum parsimony, maximum likelihood and Bayesian inference methods. The results supported an expanded clade of Mauremys, which included species formerly in Ocadia, Chinemys, Mauremys, and Annamemys; this was also reflected in the results of VNTR analysis.     

Mitogenomes of two neotropical bird species and the multiple independent origin of mitochondrial gene orders in Passeriformes

At least four mitogenome arrangements occur in Passeriformes and differences among them are derived from an initial tandem duplication involving a segment containing the control region (CR), followed by loss or reduction of some parts of this segment. However, it is still unclear how often duplication events have occurred in this bird order. In this study, the mitogenomes from two species of Neotropical passerines (Sicalis olivascens and Lepidocolaptes angustirostris) with different gene arrangements were first determined. We also estimated how often duplication events occurred in Passeriformes and if the two CR copies demonstrate a pattern of concerted evolution in Sylvioidea. One tissue sample for each species was used to obtain the mitogenomes as a byproduct using next generation sequencing. The evolutionary history of mitogenome rearrangements was reconstructed mapping these characters onto a mitogenome Bayesian phylogenetic tree of Passeriformes. Finally, we performed a Bayesian analysis for both CRs from some Sylvioidea species in order to evaluate the evolutionary process involving these two copies. Both mitogenomes described comprise 2 rRNAs, 22 tRNAs, 13 protein-codon genes and the CR. However, S. olivascens has 16,768 bp showing the ancestral avian arrangement, while L. angustirostris has 16,973 bp and the remnant CR2 arrangement. Both species showed the expected gene order compared to their closest relatives. The ancestral state reconstruction suggesting at least six independent duplication events followed by partial deletions or loss of one copy in some lineages. Our results also provide evidence that both CRs in some Sylvioidea species seem to be maintained in an apparently functional state, perhaps by concerted evolution, and that this mechanism may be important for the evolution of the bird mitogenome.     

Drosophila duplicate genes evolve new functions on the fly

Gene duplication is thought to play a key role in phenotypic innovation. While several processes have been hypothesized to drive the retention and functional evolution of duplicate genes, their genomic contributions have never been determined. We recently developed the first genome-wide method to classify these processes by comparing distances between expression profiles of duplicate genes and their ancestral single-copy orthologs. Application of our approach to spatial gene expression profiles in two Drosophila species revealed that a majority of young duplicate genes possess new functions, and that new functions are acquired rapidly—often within a few million years. Surprisingly, new functions tend to arise in younger copies of duplicate gene pairs. Moreover, we found that young duplicates are often specifically expressed in testes, whereas old duplicates are broadly expressed across several tissues, providing strong support for the hypothetical “out-of-testes” origin of new genes. In this Extra View, I discuss our findings in the context of theoretical predictions about gene duplication, with a particular emphasis on the importance of natural selection in the evolution of novel phenotypes.     

Duplication and Concerted Evolution of the Mitochondrial Control Region in the Parrot Genus Amazona

We report a duplication and rearrangement of the mitochondrialgenome involving the control region of parrots in the genusAmazona. This rearrangement results in a gene order of cytochromeb/tRNA^Thr/pND6/pGlu/CR1/tRNA^Pro/NADH dehydrogenase 6/tRNA^Glu/CR2/tRNA^Phe/12srRNA, where CR1 and CR2 refer to duplicate control regions,and pND6 and pGlu indicate presumed pseudogenes. In contrastto previous reports of duplications involving the control regionsof birds, neither copy of the parrot control region shows anyindications of degeneration. Rather, both copies contain manyof the conserved sequence features typically found in aviancontrol regions, including the goose hairpin, TASs, the F, C,and D boxes, conserved sequence box 1 (CSB1), and an apparenthomolog to the mammalian CSB3. We conducted a phylogenetic analysisof homologous portions of the duplicate control regions from21 individuals representing four species of Amazona (A. ochrocephala,A. autumnalis, A. farinosa, and A. amazonica) and Pionus chalcopterus.This analysis revealed that an individual's two control regioncopies (i.e., the paralogous copies) were typically more closelyrelated to one another than to corresponding segments of otherindividuals (i.e., the orthologous copies). The average sequencedivergence of the paralogous control region copies within anindividual was 1.4%, versus a mean value of 4.1% between controlregion orthologs representing nearest phylogenetic neighbors.No differences were found between the paralogous copies in eitherthe rate or the pattern in which the two copies accumulatedbase pair changes. This pattern suggests concerted evolutionof the two control regions, perhaps through occasional geneconversion events. We estimated that gene conversion eventsoccurred on average every 34,670 ± 18,400 years basedon pairwise distances between the paralogous control regionsequences of each individual. Our results add to the growingbody of work indicating that under some circumstances duplicatedmitochondrial control regions are retained through evolutionarytime rather than degenerating and being lost, presumably dueto selection for a small mitochondrial genome.     

Isolation and Characterization of 14 Polymorphic Microsatellite Loci in the Big-Headed Turtle (Platysternon megacephalum)

The big-headed turtle (Platysternon megacephalum) is critically endangered because of overharvesting, illegal trade, and habitat destruction. Assessment of genetic variability in existing populations becomes very important to the taxonomy and conservation of this species. Here we describe 14 microsatellite loci isolated from an enriched genomic library of the big-headed turtle, and the polymorphisms of these loci were assessed in 28 individuals from Huizhou, Heyuan, Zhaoqing, and Shaoguan of Guangdong, China. The range of polymorphism information content is 0.305–0.738, and no evidence of significant linkage disequilibrium was found among any pairs of loci. These 14 new polymorphic microsatellite loci can be used in population genetics, taxonomy, phylogeography, behavior ecology, and conservation efforts of Platysternon megacephalum.     

Generalized linear mixed model for segregation distortion analysis

Background

Concerted evolution refers to the pattern in which copies of multigene families show high intraspecific sequence homogeneity but high interspecific sequence diversity. Sequence homogeneity of these copies depends on relative rates of mutation and recombination, including gene conversion and unequal crossing over, between misaligned copies. The internally repetitive intergenic spacer (IGS) is located between the genes for the 28S and 18S ribosomal RNAs. To identify patterns of recombination and/or homogenization within IGS repeat arrays, and to identify regions of the IGS that are under functional constraint, we analyzed 13 complete IGS sequences from 10 individuals representing four species in the Daphnia pulex complex.

Results

Gene conversion and unequal crossing over between misaligned IGS repeats generates variation in copy number between arrays, as has been observed in previous studies. Moreover, terminal repeats are rarely involved in these events. Despite the occurrence of recombination, orthologous repeats in different species are more similar to one another than are paralogous repeats within species that diverged less than 4 million years ago. Patterns consistent with concerted evolution of these repeats were observed between species that diverged 8-10 million years ago. Sequence homogeneity varies along the IGS; the most homogeneous regions are downstream of the 28S rRNA gene and in the region containing the core promoter. The inadvertent inclusion of interspecific hybrids in our analysis uncovered evidence of both inter- and intrachromosomal recombination in the nonrepetitive regions of the IGS.

Conclusions

Our analysis of variation in ribosomal IGS from Daphnia shows that levels of homogeneity within and between species result from the interaction between rates of recombination and selective constraint. Consequently, different regions of the IGS are on substantially different evolutionary trajectories.     

Recombination and Speciation: Loci Near Centromeres Are More Differentiated Than Loci Near Telomeres Between Subspecies of the European Rabbit (Oryctolagus cuniculus)

Recent empirical and theoretical studies suggest that regions of restricted recombination play an important role in the formation of new species. To test this idea, we studied nucleotide variation in two parapatric subspecies of the European rabbit (Oryctolagus cuniculus). We surveyed five loci near centromeres, where recombination is expected to be suppressed, and five loci near telomeres, where recombination is expected to be higher. We analyzed this multilocus data set using a divergence-with-gene flow framework and we report three main findings. First, we estimated that these subspecies diverged ~1.8 MYA and maintained large effective population sizes (O. c. algirus N_e ≈ 1,600,000 and O. c. cuniculus N_e ≈ 780,000). Second, we rejected a strict allopatric model of divergence without gene flow; instead, high rates of gene flow were inferred in both directions. Third, we found different patterns between loci near centromeres and loci near telomeres. Loci near centromeres exhibited higher levels of linkage disequilibrium than loci near telomeres. In addition, while all loci near telomeres showed little differentiation between subspecies, three of five loci near centromeres showed strong differentiation. These results support a view of speciation in which regions of low recombination can facilitate species divergence in the presence of gene flow.     

Molecular and Cytological Characterization of Centromeric Retrotransposons in a Wild Relative of Rice, Oryza granulata

Centromeric retrotransposons (CRs) are important component of the functional centromeres of rice chromosomes. To track the evolution of the CR elements in genus Oryza, we sequenced the orthologous region of the rice centromere 8 (Cen8) in O. granulata and analyzed transposons in this region. A total of 12 bacterial artificial chromosomes (BACs) that span the centromeric region in O. granulata were sequenced. The O. granulate centromeric sequences are composed of as much as 85% of transposons, higher than any other reported eukaryotic centromeres. Ten novel LTR retrotransposon families were identified but a single retrotransposon, Gran3, constitutes nearly 43% of the centromeric sequences. Integration times of complete LTR retrotransposons indicate that the centromeric region had a massive insertion of LTR retrotransposons within 4.5 million year (Myr), which indicates a recent expansion of the centromere in O. granulata after the radiation of the Oryza genus. Two retrotransposon families, OGRetro7 and OGRetro9, show sequence similarity with the canonical CRs from rice and maize. Both OGRetro7 and OGRetro9 are highly concentrated in the centromeres of O. granulata chromosomes. Furthermore, strong hybridization signals were detected in all Oryza species but in O. brachyantha with the OGRetro7 and OGRetro9 probes. Characterization of the centromeric retrotransposons in O. granulata confirms the conservation of the CRs in the Oryza genus and provides a resource for comparative analysis of centromeres and centromere evolution among the Oryza genus and other genomes.     

Mitochondrial DNA-like sequences in the nuclear genome of the opossum genus Didelphis (Marsupialia: Didelphidae).

This study reports the occurrence of mtDNA-like sequences in the nuclear genome of the opossum genus Didelphis (Didelphidae, Marsupialia). A specific primer pair designed to amplify a region encompassing a 3' terminal 118 bp region of the cytochrome b gene, the Thr and Pro tRNA genes, and a 489 bp region of the D-loop of the D. virginiana mtDNA, was used in highly stringent PCR reactions. These PCR reactions resulted in several fragments per individual varying in size from 259 bp to 1 kb. The sequencing of some of these fragments showed the occurrence of paralogous mtDNA-like sequences among the PCR amplified fragments. Analyses of qualitative aspects of these sequences, their transition/transversion ratios, and phylogenetic relationships were conclusive in showing the occurrence of mtDNA-like sequences in the nuclear genome of the genus Didelphis. Comparisons and phylogenetic analysis of orthologous mtDNA from the four Didelphis species and paralogous nuclear sequences suggested that mtDNA migration to the nuclear genome occurred more than once in Didelphis evolution.     

Mitochondrial DNA Diversity of Modern,Ancient and Wild Sheep (Ovis gmelinii anatolica) from Turkey: New Insights on the Evolutionary History of Sheep

In the present study, to contribute to the understanding of the evolutionary history of sheep, the mitochondrial (mt) DNA polymorphisms occurring in modern Turkish native domestic (n = 628), modern wild (Ovis gmelinii anatolica) (n = 30) and ancient domestic sheep from Oylum Höyük in Kilis (n = 33) were examined comparatively with the accumulated data in the literature. The lengths (75 bp/76 bp) of the second and subsequent repeat units of the mtDNA control region (CR) sequences differentiated the five haplogroups (HPGs) observed in the domestic sheep into two genetic clusters as was already implied by other mtDNA markers: the first cluster being composed of HPGs A, B, D and the second cluster harboring HPGs C, E.To manifest genetic relatedness between wild Ovis gmelinii and domestic sheep haplogroups, their partial cytochrome B sequences were examined together on a median-joining network. The two parallel but wider aforementioned clusters were observed also on the network of Ovis gmelenii individuals, within which domestic haplogroups were embedded. The Ovis gmelinii wilds of the present day appeared to be distributed on two partially overlapping geographic areas parallel to the genetic clusters that they belong to (the first cluster being in the western part of the overall distribution). Thus, the analyses suggested that the domestic sheep may be the products of two maternally distinct ancestral Ovis gmelinii populations.Furthermore, Ovis gmelinii anatolica individuals exhibited a haplotype of HPG A (n = 22) and another haplotype (n = 8) from the second cluster which was not observed among the modern domestic sheep. HPG E, with the newly observed members (n = 11), showed signs of expansion. Studies of ancient and modern mtDNA suggest that HPG C frequency increased in the Southeast Anatolia from 6% to 22% some time after the beginning of the Hellenistic period, 500 years Before Common Era (BCE).     

Estimates of linkage disequilibrium and effective population size in rainbow trout

Background

Mitochondrial DNA (mtDNA) is widely used in population genetic and phylogenetic studies in animals. However, such studies can generate misleading results if the species concerned contain nuclear copies of mtDNA (Numts) as these may amplify in addition to, or even instead of, the authentic target mtDNA. The aim of this study was to determine if Numts are present in Aedes aegypti mosquitoes, to characterise any Numts detected, and to assess the utility of using mtDNA for population genetics studies in this species.

Results

BLAST searches revealed large numbers of Numts in the Ae. aegypti nuclear genome on 146 supercontigs. Although the majority are short (80% < 300 bp), some Numts are almost full length mtDNA copies. These long Numts are not due to misassembly of the nuclear genome sequence as the Numt-nuclear genome junctions could be recovered by amplification and sequencing. Numt evolution appears to be a complex process in Ae. aegypti with ongoing genomic integration, fragmentation and mutation and the secondary movement of Numts within the nuclear genome. The PCR amplification of the putative mtDNA nicotinamide adenine dinucleotide dehydrogenase subunit 4 (ND4) gene from 166 Southeast Asian Ae. aegypti mosquitoes generated a network with two highly divergent lineages (clade 1 and clade 2). Approximately 15% of the ND4 sequences were a composite of those from each clade indicating Numt amplification in addition to, or instead of, mtDNA. Clade 1 was shown to be composed at least partially of Numts by the removal of clade 1-specific bases from composite sequences following enrichment of the mtDNA. It is possible that all the clade 1 sequences in the network were Numts since the clade 2 sequences correspond to the known mitochondrial genome sequence and since all the individuals that produced clade 1 sequences were also found to contain clade 2 mtDNA-like sequences using clade 2-specific primers. However, either or both sets of clade sequences could have Numts since the BLAST searches revealed two long Numts that match clade 2 and one long Numt that matches clade 1. The substantial numbers of mutations in cloned ND4 PCR products also suggest there are both recently-derived clade 1 and clade 2 Numt sequences.

Conclusion

We conclude that Numts are prevalent in Ae. aegypti and that it is difficult to distinguish mtDNA sequences due to the presence of recently formed Numts. Given this, future population genetic or phylogenetic studies in Ae. aegypti should use nuclear, rather than mtDNA, markers.     

Two Sympatric Phylogroups of the Chinese Water Deer (Hydropotes inermis) Identified by Mitochondrial DNA Control Region and Cytochrome b Gene Analyses

Mitochondrial DNA (mtDNA) control region (927 bp) and cytochrome b gene (1,140 bp) sequences of the Chinese water deer (Hydropotes inermis) from China and Korea were obtained to examine the taxonomic status of two subspecies, H. i. inermis from China and H. i. argyropus from Korea. Two sympatric mtDNA clades (a major clade from China and Korea and a minor clade from Korea) with an average genetic distance of 2.1% in the control region and 1.3% in the cytochrome b gene were detected. These findings are not consistent with the current classification by pelage color. We propose a reconsideration of the validity of the subspecies designation by the statistical comparison of morphological characters including body color. The major common mtDNA phylogroup in the two allopatric subspecies could be explained by the contiguous distribution of the Chinese water deer from east China to Korea until recent years. The restriction in the range and number of the Chinese subspecies after the last glacier might have caused the disappearance of the minor phylogroup in China. The taxonomic status of the two groups in Korea should be clarified using nuclear DNA marker analyses as well as morphological characters including pelage color.