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

Structure and variation of the mitochondrial genome of fishes

Background

The mitochondrial (mt) genome has been used as an effective tool for phylogenetic and population genetic analyses in vertebrates. However, the structure and variability of the vertebrate mt genome are not well understood. A potential strategy for improving our understanding is to conduct a comprehensive comparative study of large mt genome data. The aim of this study was to characterize the structure and variability of the fish mt genome through comparative analysis of large datasets.

Results

An analysis of the secondary structure of proteins for 250 fish species (248 ray-finned and 2 cartilaginous fishes) illustrated that cytochrome c oxidase subunits (COI, COII, and COIII) and a cytochrome bc1 complex subunit (Cyt b) had substantial amino acid conservation. Among the four proteins, COI was the most conserved, as more than half of all amino acid sites were invariable among the 250 species. Our models identified 43 and 58 stems within 12S rRNA and 16S rRNA, respectively, with larger numbers than proposed previously for vertebrates. The models also identified 149 and 319 invariable sites in 12S rRNA and 16S rRNA, respectively, in all fishes. In particular, the present result verified that a region corresponding to the peptidyl transferase center in prokaryotic 23S rRNA, which is homologous to mt 16S rRNA, is also conserved in fish mt 16S rRNA. Concerning the gene order, we found 35 variations (in 32 families) that deviated from the common gene order in vertebrates. These gene rearrangements were mostly observed in the area spanning the ND5 gene to the control region as well as two tRNA gene cluster regions (IQM and WANCY regions). Although many of such gene rearrangements were unique to a specific taxon, some were shared polyphyletically between distantly related species.

Conclusions

Through a large-scale comparative analysis of 250 fish species mt genomes, we elucidated various structural aspects of the fish mt genome and the encoded genes. The present results will be important for understanding functions of the mt genome and developing programs for nucleotide sequence analysis. This study demonstrated the significance of extensive comparisons for understanding the structure of the mt genome.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3054-y) contains supplementary material, which is available to authorized users.

Novel duplication pattern of the mitochondrial control region in Cantor's Giant softshell turtle Pelochelys cantorii

Cantor's Giant Softshell Turtle, Pelochelys cantorii has become one of the most critically endangered species in the world. When comparative analyses of the P. cantorii complete mitochondrial genome sequences were conducted, we discovered a duplication of a segment of the control region in the mitochondrial genome of P. cantorii. The duplication is characterized by two copies of conserved sequence box 2 (CSB2) and CSB3 in a single control region. In contrast to previous reports of duplications involving the control regions of other animals, this particular pattern of duplications appears to be unique to P. cantorii. Copies of the CSB2 and CSB3 show many of the conserved sequence features typically found in mitochondrial control regions, and rare differences were found between the paralogous copies. Using the primer design principle of simple sequence repeats (SSR) and the reference sequence of the duplicated CSBs, specific primers were designed to amplify the duplicated CSBs. These primers were validated among different individuals and populations of P. cantorii. This unique duplication structure suggests the two copies of the CSB2 and CSB3 may have arisen through occasional tandem duplication and subsequent concerted evolution.

The role of mitochondrial DNA copy number, variants, and haplotypes in farm animal developmental outcome

The vast majority of cellular energy is generated through the process of oxidative phosphorylation, which takes place in the electron transport chain in the mitochondria. The electron transport chain is encoded by 2 genomes, the chromosomal and the mitochondrial genomes. Mitochondrial DNA is associated with a number of traits, which include tolerance to heat, growth and physical performance, meat and milk quality, and fertility. Mitochondrial genomes can be clustered into groups known as mtDNA haplotypes. Mitochondrial DNA haplotypes are a potential genetic source for manipulating phenotypes in farm animals. The use of assisted reproductive technologies, such as nuclear transfer, allows favorable chromosomal genetic traits to be mixed and matched with sought after mtDNA haplotype traits. As a result super breeds can be generated.

Characterization of the complete mitochondrial genome of Cynoglossus gracilis and a comparative analysis with other Cynoglossinae fishes

Mitochondrial genomes can provide basic information for phylogenetic analysis and evolutionary studies. We present here the mitochondrial genome of Cynoglossus gracilis, which is 16,565bp in length. Numerous distinct regions were identified, including 13 protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes, a light-strand replication origin, and a control region. Interestingly, we detected rearrangement of genes in C. gracilis, including a control region translocation, tRNA(Gln) gene inversion, and tRNA(Ile) gene shuffling. Additionally, a phylogenetic analysis based on the nucleotide sequences of the 13 PCGs using maximum likelihood and Bayesian inference methods reveals that C. gracilis is closely related to Cynoglossus semilaevis. This study provides important mitogenomic data for analyzing phylogenetic relationships in the Cynoglossinae.

Identification of species origin of meat and meat products on the DNA basis: a review

The adulteration/substitution of meat has always been a concern for various reasons such as public health, religious factors, wholesomeness, and unhealthy competition in meat market. Consumer should be protected from these malicious practices of meat adulterations by quick, precise, and specific identification of meat animal species. Several analytical methodologies have been employed for meat speciation based on anatomical, histological, microscopic, organoleptic, chemical, electrophoretic, chromatographic, or immunological principles. However, by virtue of their inherent limitations, most of these techniques have been replaced by the recent DNA-based molecular techniques. In the last decades, several methods based on polymerase chain reaction have been proposed as useful means for identifying the species origin in meat and meat products, due to their high specificity and sensitivity, as well as rapid processing time and low cost. This review intends to provide an updated and extensive overview on the DNA-based methods for species identification in meat and meat products.

Complete mitochondrial genome of the golden takin (Budorcas taxicolor bedfordi)

The complete mitochondrial genome sequence of Budorcas taxicolor bedfordi was 16 662-bp long, containing 13 protein-coding genes (PCGs), 22 tRNA genes, 2 rRNA genes and a control region. All of the PCGs begin with the typical ATN start codon. Stop codons TAA and AGA are present in the PCGs; exceptions are ND2, COIII, ND3 and ND4, which possess incomplete termination codon (T––). Secondary structure prediction of the 22 tRNA genes revealed the absence of the DHU arm in tRNASer(AGN).

Strictly conserved tri-nucleotide motif "CAT" is associated with TAS DNA protein-binding sites in human mitochondrial DNA control region

The mitochondrial DNA (mtDNA) control region is a highly variable segment that contains functional elements that control mtDNA transcription and replication. By analysis of the polymorphic nucleotide spectrum of that segment, we aimed to identify the most conserved sites that should be associated with these elements. For that aim, we analyzed 50 033 human mtDNA control region sequences (mtDNA positions 16 066-16 374). We identified 10 conserved tri-nucleotides, one conserved tetra-nucleotide, and one conserved penta-nucleotide, containing six repetitions of the motif CAT, and two of its complement motif ATG (p value < 2 × 10 ^- 4). Three other appearances of the tri-nucleotide CAT were almost perfectly preserved. The positions of the preserved CAT elements are associated with the location of previously identified termination-associated sequences (TAS) which are the binding locations for proteins involved in mtDNA replication. We, therefore, hypothesize that the CAT tri-nucleotide elements within the control region may be the binding sites for TAS proteins and are directly involved in mtDNA transcription and replication.

Extremely low genetic diversity of mtDNA control region and remarkable population differentiation of Ichthyophis bannanicus (Amphibia: Gymnophiona)

It is widely accepted that the mitochondrial DNA (mtDNA) control region (CR) gene evolves more quickly than protein-encoding genes, such as ND2 and Cyt b, with few exceptions. However, some species have a different evolution pattern. In the present study, we sequenced the mtDNA CR partial sequences (454 base pairs) of 142 individuals from five sampling sites of Ichthyophis bannanicus and compared the genetic diversity and structure with the information from the NADH dehydrogenase subunit 2 (ND2) and Cytochrome b (Cyt b) genes within this species. Extremely low genetic diversity was found in the mtDNA CR compared with those of the ND2 and Cyt b genes. These results showed that the relatively mean clock rate of the CR was broadly lower than those of the ND2 (about 2.55 times) and Cyt b (about 3.14 times) genes. Despite the extremely low genetic diversity of CR, the population structure analysis identified two groups, Xishuangbanna and Northern Vietnam-Yulin-Yangchun-Deqing, which indicated that the Red River systems may have acted as gene-flow barriers for I. bannanicus.

Complete mitochondrial DNA sequence of the endangered fish (Bahaba taipingensis): Mitogenome characterization and phylogenetic implications

To understand the systematic status of Bahabataipingensis within Sciaenidae, the complete mitochondrial genome (mitogenome) sequence of Chinese bahaba has recently been determined by long PCR and primer walking methods. The complete mitochondrial genome is 16500 bp in length and contains 37 mitochondrial genes (13 protein-coding genes, 2 ribosomal RNA genes and 22 transfer RNA genes) as well as a control region (CR) as other bony fishes. Within the control region, we identified the

extended termination associated sequence domain

(ETAS), the central conserved sequence block domain (CSB-D, SCB-E and CSB-F) and the conserved sequence block domain (CSB-1, CSB-2 and CSB-3). Phylogenetic analyses revealed that Bahabataipingensis is more closely related to Pseudosciaeniae than Argyrosominae and Sciaeninae. Additionally, Bahabataipingensis is the sister taxon of Miichthysmiiuy, and those two are sister to Collichthys plus Larimichthys.

Complete mitochondrial genome and phylogenic analysis of the mudskipper Scartelaos gigas (Perciformes, Gobiidae)

In this paper, the complete mitochondrial genome sequence of Scartelaos gigas was firstly determined. The circular genome (16 717 bp) comprises 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, and 1 control region. The overall base composition of S. gigas is 28.9% for C, 28.3% for A, 26.4% for T, 16.4% for G, with a slight A + T bias of 54.7%. In the control region, the termination-associated sequence and conserved sequence block domains were found, but the tandem repeat structure was not found. It has the typical vertebrate mitochondrial gene arrangement. The phylogenic analysis using the Neighbor-Joining method showed that the fishes belonging to Gobiidae, Odontoburidae, and Eleotridae formed three branches grouped with other fishes into one clade which separated from the mammals. We hope that the results from the present study will provide useful molecular information for the further studies on genetic structure and demographic history of S. gigas.

Rearrangement and evolution of mitochondrial genomes in parrots

Mitochondrial genome rearrangements that result in control region duplication have been described for a variety of birds, but the mechanisms leading to their appearance and maintenance remain unclear, and their effect on sequence evolution has not been explored. A recent survey of mitochondrial genomes in the Psittaciformes (parrots) found that control region duplications have arisen independently at least six times across the order. We analyzed complete mitochondrial genome sequences from 20 parrot species, including representatives of each lineage with control region duplications, to document the gene order changes and to examine effects of genome rearrangements on patterns of sequence evolution. The gene order previously reported for Amazona parrots was found for four of the six independently derived genome rearrangements, and a previously undescribed gene order was found in Prioniturus luconensis, representing a fifth clade with rearranged genomes; the gene order resulting from the remaining rearrangement event could not be confirmed. In all rearranged genomes, two copies of the control region are present and are very similar at the sequence level, while duplicates of the other genes involved in the rearrangement show signs of degeneration or have been lost altogether. We compared rates of sequence evolution in genomes with and without control region duplications and did not find a consistent acceleration or deceleration associated with the duplications. This could be due to the fact that most of the genome rearrangement events in parrots are ancient, and additionally, to an effect of body size on evolutionary rate that we found for mitochondrial but not nuclear sequences. Base composition analyses found that relative to other birds, parrots have unusually strong compositional asymmetry (AT- and GC-skew) in their coding sequences, especially at fourfold degenerate sites. Furthermore, we found higher AT skew in species with control region duplications. One potential cause for this compositional asymmetry is that parrots have unusually slow mtDNA replication. If this is the case, then any replicative advantage provided by having a second control region could result in selection for maintenance of both control regions once duplicated.

Regulation of DNA replication at the end of the mitochondrial D-loop involves the helicase TWINKLE and a conserved sequence element

The majority of mitochondrial DNA replication events are terminated prematurely. The nascent DNA remains stably associated with the template, forming a triple-stranded displacement loop (D-loop) structure. However, the function of the D-loop region of the mitochondrial genome remains poorly understood. Using a comparative genomics approach we here identify two closely related 15 nt sequence motifs of the D-loop, strongly conserved among vertebrates. One motif is at the D-loop 5'-end and is part of the conserved sequence block 1 (CSB1). The other motif, here denoted coreTAS, is at the D-loop 3'-end. Both these sequences may prevent transcription across the D-loop region, since light and heavy strand transcription is terminated at CSB1 and coreTAS, respectively. Interestingly, the replication of the nascent D-loop strand, occurring in a direction opposite to that of heavy strand transcription, is also terminated at coreTAS, suggesting that coreTAS is involved in termination of both transcription and replication. Finally, we demonstrate that the loading of the helicase TWINKLE at coreTAS is reversible, implying that this site is a crucial component of a switch between D-loop formation and full-length mitochondrial DNA replication.

The mitochondrial genome of the Saunders's gull Chroicocephalus saundersi (Charadriiformes: Laridae) and a higher phylogeny of shorebirds (Charadriiformes)

The complete mitogenome of Chroicocephalus saundersi was characterized and compared with the 6 published Charadriiformes mitogenomes. The mitogenome of C. saundersi is a closed circular molecule 16,739 bp in size, and contains 37 genes and a control region. The AT and GC skews are positive and negative, respectively, and in agreement with those of the other Charadriiformes mitogenomes. The mitogenome of C. saundersi contains 3 start codons (ATG, GTG, and ATT), 4 stop codons (TAA, TAG, AGG, and AGA), and an incomplete stop codon (T-) in 13 PCGs. A codon usage analysis of all available Charadriiformes mitogenomes showed that the ATG (78%) and TAA (50.5%) were the most common start codon and stop codon, respectively. An unusual start codon, ATT, is commonly found in the ND3s of Charadriiformes mitogenomes, whereas the more common start codons, ATC and ATA, are rarely found. In all the Laridae species, one extra cytosine was inserted at position 174 in ND3. The control region of C. saundersi is 1180-bp long, with a nucleotide composition of 30.2% A, 28.6% T, 27.3% C, and 14.0% G. Variable numbers of tandem repeats (VNTRs) with nine copies of the 10 bp repeat sequence (AACAACAAAC) are found within the CSB domain of the control region. The ML/BI analyses, based on the amino acids of the 13 mitochondrial PCGs, strongly support the monophyly of the order Charadriiformes, with the suborder Lari considered sister to the Scolopaci, which is in turn a sister group to the suborder Charadrii.

Complete Taiwanese Macaque (Macaca cyclopis) Mitochondrial Genome: Reference-Assisted de novo Assembly with Multiple k-mer Strategy

The Taiwanese (Formosan) macaque (Macaca cyclopis) is the only nonhuman primate endemic to Taiwan. This primate species is valuable for evolutionary studies and as subjects in medical research. However, only partial fragments of the mitochondrial genome (mitogenome) of this primate species have been sequenced, not mentioning its nuclear genome. We employed next-generation sequencing to generate 2 x 90 bp paired-end reads, followed by reference-assisted de novo assembly with multiple k-mer strategy to characterize the M. cyclopis mitogenome. We compared the assembled mitogenome with that of other macaque species for phylogenetic analysis. Our results show that, the M. cyclopis mitogenome consists of 16,563 nucleotides encoding for 13 protein-coding genes, 2 ribosomal RNAs and 22 transfer RNAs. Phylogenetic analysis indicates that M. cyclopis is most closely related to M. mulatta lasiota (Chinese rhesus macaque), supporting the notion of Asia-continental origin of M. cyclopis proposed in previous studies based on partial mitochondrial sequences. Our work presents a novel approach for assembling a mitogenome that utilizes the capabilities of de novo genome assembly with assistance of a reference genome. The availability of the complete Taiwanese macaque mitogenome will facilitate the study of primate evolution and the characterization of genetic variations for the potential usage of this species as a non-human primate model for medical research.

The complete mitochondrial genome sequence of the little egret (Egretta garzetta)

Many phylogenetic questions in the Ciconiiformes remain unresolved and complete mitogenome data are urgently needed for further molecular investigation. In this work, we determined the complete mitogenome sequence of the little egret (Egretta garzetta). The genome was 17,361 bp in length and the gene organization was typical of other avian mtDNA. In protein-coding genes (PCGs), a C insertion was found in ND3, and COIII and ND4 terminated with incomplete stop codons (T). tRNA-Val and tRNA-Ser (AGY) were unable to fold into canonical cloverleaf secondary structures because they had lost the DHU arms. Long repetitive sequences consisting of five types of tandem repeats were found at the 3' end of Domain III in the control region. A phylogenetic analysis of 11 species of Ciconiiformes was done using complete mitogenome data and 12 PCGs. The tree topologies obtained with these two strategies were identical, which strongly confirmed the monophyly of Ardeidae, Threskiorothidae and Ciconiidae. The phylogenetic analysis also revealed that Egretta was more closely related to Ardea than to Nycticorax in the Ardeidae, and Platalea was more closely related to Threskiornis than to Nipponia in the Threskiornithidae. These findings contribute to our understanding of the phylogenetic relationships of Ciconiiformes based on complete mitogenome data.

Do cryptic species exist in Hoplobatrachus rugulosus? An examination using four nuclear genes, the cyt b gene and the complete MT genome

he Chinese tiger frog Hoplobatrachus rugulosus is widely distributed in southern China, Malaysia, Myanmar, Thailand, and Vietnam. It is listed in Appendix II of CITES as the only Class II nationally-protected frog in China. The bred tiger frog known as the Thailand tiger frog, is also identified as H. rugulosus. Our analysis of the Cyt b gene showed high genetic divergence (13.8%) between wild and bred samples of tiger frog. Unexpected genetic divergence of the complete mt genome (14.0%) was also observed between wild and bred samples of tiger frog. Yet, the nuclear genes (NCX1, Rag1, Rhod, Tyr) showed little divergence between them. Despite this and their very similar morphology, the features of the mitochondrial genome including genetic divergence of other genes, different three-dimensional structures of ND5 proteins, and gene rearrangements indicate that H. rugulosus may be a cryptic species complex. Using Bayesian inference, maximum likelihood, and maximum parsimony analyses, Hoplobatrachus was resolved as a sister clade to Euphlyctis, and H. rugulosus (BT) as a sister clade to H. rugulosus (WT). We suggest that we should prevent Thailand tiger frogs (bred type) from escaping into wild environments lest they produce hybrids with Chinese tiger frogs (wild type).

Complete mitochondrial genome of Acrossocheilus parallens (Cypriniformes, Barbinae)

The omnivorous fish Acrossocheilus parallens (Cypriniformes, Barbinae) is endemic to China. In this paper, the complete mitochondrial genome sequence of A. parallens was firstly determined. The circle genome (16,592 bp) comprises 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes and 1 control region. The overall base composition of A. parallens is 31.0% for A, 24.7% for T, 28.0% for C, and 16.3% for G, with a slight A + T bias of 55.7%. The termination-associated sequence and conserved sequence block domains were found in the control region. It has the typical vertebrate mitochondrial gene arrangement.

The complete mitochondrial genome of silver fox (Caniformia: Canidae)

Silver fox is color variant of Vulpes vulpes. At present, there are few studies on phylogeny of Canidae and Caniformia. In this article, we determined and described the complete mitogenome of silver fox for the first time, which is 16,723 bp in length, containing 13 protein-coding genes, 2 rRNA genes, 22 tRNA genes, one origin of replication on the light-strand (OL) and a putative control region (CR). The overall base composition is 31.4% A, 27.9% T, 26.0% C, 14.7% G, respectively, with a AT bias (59.3%). Ten protein-coding genes use the initiation codon ATG while ND2, ND3 and ND5 use ATA. Most of them have TAA as the stop codon, except ND2 uses TAG, Cytb uses AGA, and COX3, ND3, ND4 use an incomplete stop codon TA. The information is expected to provide useful molecular data for further taxonomic and phylogenetic studies of Canidae and Caniformia.

Comparative mitochondrial genomics and phylogenetic relationships of the Crossoptilon species (Phasianidae, Galliformes)

BACKGROUND

Phasianidae is a family of Galliformes containing 38 genera and approximately 138 species, which is grouped into two tribes based on their morphological features, the Pheasants and Partridges. Several studies have attempted to reconstruct the phylogenetic relationships of the Phasianidae, but many questions still remain unaddressed, such as the taxonomic status and phylogenetic relationships among Crossoptilon species. The mitochondrial genome (mitogenome) has been extensively used to infer avian genetic diversification with reasonable resolution. Here, we sequenced the entire mitogenomes of three Crossoptilon species (C. harmani, C. mantchuricum and C. crossoptilon) to investigate their evolutionary relationship among Crossoptilon species.

RESULTS

The complete mitogenomes of C. harmani, C. mantchuricum and C. crossoptilon are 16682 bp, 16690 bp and 16680 bp in length, respectively, encoding a standard set of 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, and a putative control region. C. auritum and C. mantchuricum are more closely related genetically, whereas C. harmani is more closely related to C. crossoptilon. Crossoptilon has a closer relationship with Lophura, and the following phylogenetic relationship was reconstructed: ((Crossoptilon + Lophura) + (Phasianus + Chrysolophus)). The divergence time between the clades C. harmani-C. crossoptilon and C. mantchuricum-C. auritum is consistent with the uplift of the Tibetan Plateau during the Tertiary Pliocene. The Ka/Ks analysis showed that atp8 gene in the Crossoptilon likely experienced a strong selective pressure in adaptation to the plateau environment.

CONCLUSIONS

C. auritum with C. mantchuricum and C. harmani with C. crossoptilon form two pairs of sister groups. The genetic distance between C. harmani and C. crossoptilon is far less than the interspecific distance and is close to the intraspecific distance of Crossoptilon, indicating that C. harmani is much more closely related to C. crossoptilon. Our mito-phylogenomic analysis supports the monophyly of Crossoptilon and its closer relationship with Lophura. The uplift of Tibetan Plateau is suggested to impact the divergence between C. harmani-C. crossoptilon clade and C. mantchuricum-C. auritum clade during the Tertiary Pliocene. Atp8 gene in the Crossoptilon species might have experienced a strong selective pressure for adaptation to the plateau environment.

Characterization of the Macropodus opercularis complete mitochondrial genome and family Channidae taxonomy using Illumina-based de novo transcriptome sequencing

In this study, the complete mitochondrial genome of Macropodus opercularis was sequenced using Illumina-based de novo transcriptome technology and annotated using bioinformatic tools. The circular mitochondrial genome was 16,496bp in length and contained two ribosomal RNAs, 13 protein-coding genes, 22 transfer RNA genes, and the control region. The gene composition and order were similar to suborder Anabantoidei. Phylogenetic analyses using concatenated amino acid and nucleotide sequences of the 13 protein-coding genes with two different methods (Neighbor-joining and Bayesian analysis) both highly supported the close relationship of M. opercularis to M. ocellatus, consistent with previous classifications based on morphological and molecular studies. Furthermore, family Channidae and Parachanna insignis were clustered in the same clade. Our results supported the inclusion of family Channidae in suborder Channoidei. The complete mitochondrial genome of M. opercularis will provide genetic markers for better understanding species identification, population genetics and phylogeographics of freshwater fishes.

Complete mitochondrial genome of Cuora trifasciata (Chinese three-striped box turtle), and a comparative analysis with other box turtles

Cuora trifasciata has become one of the most critically endangered species in the world. The complete mitochondrial genome of C. trifasciata (Chinese three-striped box turtle) was determined in this study. Its mitochondrial genome is a 16,575-bp-long circular molecule that consists of 37 genes that are typically found in other vertebrates. And the basic characteristics of the C. trifasciata mitochondrial genome were also determined. Moreover, a comparison of C. trifasciata with Cuora cyclornata, Cuora pani and Cuora aurocapitata indicated that the four mitogenomics differed in length, codons, overlaps, 13 protein-coding genes (PCGs), ND3, rRNA genes, control region, and other aspects. Phylogenetic analysis with Bayesian inference and maximum likelihood based on 12 protein-coding genes of the genus Cuora indicated the phylogenetic position of C. trifasciata within Cuora. The phylogenetic analysis also showed that C. trifasciata from Vietnam and China formed separate monophyletic clades with different Cuora species. The results of nucleotide base compositions, protein-coding genes and phylogenetic analysis showed that C. trifasciata from these two countries may represent different Cuora species.

Gene rearrangements in gekkonid mitochondrial genomes with shuffling, loss, and reassignment of tRNA genes

BACKGROUND

Vertebrate mitochondrial genomes (mitogenomes) are 16-18 kbp double-stranded circular DNAs that encode a set of 37 genes. The arrangement of these genes and the major noncoding region is relatively conserved through evolution although gene rearrangements have been described for diverse lineages. The tandem duplication-random loss model has been invoked to explain the mechanisms of most mitochondrial gene rearrangements. Previously reported mitogenomic sequences for geckos rarely included gene rearrangements, which we explore in the present study.

RESULTS

We determined seven new mitogenomic sequences from Gekkonidae using a high-throughput sequencing method. The Tropiocolotes tripolitanus mitogenome involves a tandem duplication of the gene block: tRNAArg, NADH dehydrogenase subunit 4L, and NADH dehydrogenase subunit 4. One of the duplicate copies for each protein-coding gene may be pseudogenized. A duplicate copy of the tRNAArg gene appears to have been converted to a tRNAGln gene by a C to T base substitution at the second anticodon position, although this gene may not be fully functional in protein synthesis. The Stenodactylus petrii mitogenome includes several tandem duplications of tRNALeu genes, as well as a translocation of the tRNAAla gene and a putative origin of light-strand replication within a tRNA gene cluster. Finally, the Uroplatus fimbriatus and U. ebenaui mitogenomes feature the apparent loss of the tRNAGlu gene from its original position. Uroplatus fimbriatus appears to retain a translocated tRNAGlu gene adjacent to the 5' end of the major noncoding region.

CONCLUSIONS

The present study describes several new mitochondrial gene rearrangements from Gekkonidae. The loss and reassignment of tRNA genes is not very common in vertebrate mitogenomes and our findings raise new questions as to how missing tRNAs are supplied and if the reassigned tRNA gene is fully functional. These new examples of mitochondrial gene rearrangements in geckos should broaden our understanding of the evolution of mitochondrial gene arrangements.

The inversion of the Control Region in three mitogenomes provides further evidence for an asymmetric model of vertebrate mtDNA replication

Mitochondrial genomes are known to have a strong strand-specific compositional bias that is more pronounced at fourfold redundant sites of mtDNA protein-coding genes. This observation suggests that strand asymmetries, to a large extent, are caused by mutational asymmetric mechanisms. In vertebrate mitogenomes, replication and not transcription seems to play a major role in shaping compositional bias. Hence, one can better understand how mtDNA is replicated – a debated issue – through a detailed picture of mitochondrial genome evolution. Here, we analyzed the compositional bias (AT and GC skews) in protein-coding genes of almost 2,500 complete vertebrate mitogenomes. We were able to identify three fish mitogenomes with inverted AT/GC skew coupled with an inversion of the Control Region. These findings suggest that the vertebrate mitochondrial replication mechanism is asymmetric and may invert its polarity, with the leading-strand becoming the lagging-strand and vice-versa, without compromising mtDNA maintenance and expression. The inversion of the strand-specific compositional bias through the inversion of the Control Region is in agreement with the strand-displacement model but it is also compatible with the RITOLS model of mtDNA replication.

Identification of species-specific nuclear insertions of mitochondrial DNA (numts) in gorillas and their potential as population genetic markers

The first hyper-variable region (HV1) of the mitochondrial control region (MCR) has been widely used as a molecular tool in population genetics, but inadvertent amplification of nuclear translocated copies of mitochondrial DNA (numts) in gorillas has compromised the use of mitochondrial DNA in population genetic studies. At least three putative classes (I, II, III) of gorilla-specific HV1 MCR numts have been uncovered over the past decade. However, the number, size and location of numt loci in gorillas and other apes are completely unknown. Furthermore, little work to date has assessed the utility of numts as candidate population genetic markers. In the present study, we screened Bacterial Artificial Chromosome (BAC) genomic libraries in the chimpanzee and gorilla to compare patterns of mitochondrial-wide insertion in both taxa. We conducted an intensive BLAST search for numts in the gorilla genome and compared the prevalence of numt loci originating from the MCR with other great ape taxa. Additional gorilla-specific MCR numts were retrieved either through BAC library screens or using an anchored-PCR (A-PCR) amplification using genomic DNA from five unrelated gorillas. Locus-specific primers were designed to identify numt insertional polymorphisms and evaluate their potential as population genetic markers. Mitochondrial-wide surveys of chimpanzee and gorilla BACs showed that the number of numts does not differ between these two taxa. However, MCR numts are more abundant in chimpanzees than in other great apes. We identified and mapped 67 putative gorilla-specific numts, including two that contain the entire HV1 domain, cluster with sequences from two numt classes (I, IIb) and will likely co-amplify with mitochondrial sequences using most published HV1 primers. However, phylogenetic analysis coupled with post-hoc analysis of mitochondrial variation can successfully differentiate nuclear sequences. Insertional polymorphisms were evident in three out of five numts examined, indicating their potential utility as molecular markers. Taken together, these findings demonstrate the potentially powerful insight that numts could make in uncovering population history in gorillas and other mammals.

Control region length dynamics potentially drives amino acid evolution in tarsier mitochondrial genomes

Patterns and processes of molecular evolution critically influence inferences in phylogeny and phylogeography. Within primates, a shift in evolutionary rates has been identified as the rationale for contrasting findings from mitochondrial and nuclear DNA studies as to the position of Tarsius. While the latter now seems settled, we sequenced complete mitochondrial genomes of three Sulawesi tarsiers (Tarsius dentatus, T. lariang, and T. wallacei) and analyzed substitution rates among tarsiers and other primates to infer driving processes of molecular evolution. We found substantial length polymorphism of the D-loop within tarsier individuals, but little variation of predominant lengths among them, regardless of species. Length variation was due to repetitive elements in the CSB domain-minisatellite motifs of 35 bp length and microsatellite motifs of 6 bp length. Amino acid evolutionary rates were second highest among major primate taxa relative to nucleotide substitution rates. We observed many radical possibly function-altering amino acid changes that were rarely driven by positive selection and thus potentially slightly deleterious or neutral. We hypothesize that the observed pattern of an increased amino acid evolutionary rate in tarsier mitochondrial genomes may be caused by hitchhiking of slightly deleterious mutations with favored D-loop length variants selected for maximizing replication success within the cell or the mitochondrion.

Complete sequence of the mitochondrial genome of Odontamblyopus rubicundus (Perciformes: Gobiidae): genome characterization and phylogenetic analysis

Odontamblyopus rubicundus is a species of gobiid fishes, inhabits muddy-bottomed coastal waters. In this paper, the first complete mitochondrial genome sequence of O. rubicundus is reported. The complete mitochondrial genome sequence is 17119 bp in length and contains 13 protein-coding genes, two rRNA genes, 22 tRNA genes, a control region and an L-strand origin as in other teleosts. Most mitochondrial genes are encoded on H-strand except for ND6 and seven tRNA genes. Some overlaps occur in protein-coding genes and tRNAs ranging from 1 to 7 bp. The possibly nonfunctional L-strand origin folded into a typical stem-loop secondary structure and a conserved motif (5'-GCCGG-3') was found at the base of the stem within the tRNACys gene. The TAS, CSB-2 and CSB-3 could be detected in the control region. However, in contrast to most of other fishes, the central conserved sequence block domain and the CSB-1 could not be recognized in O. rubicundus, which is consistent with Acanthogobius hasta (Gobiidae). In addition, phylogenetic analyses based on different sequences of species of Gobiidae and different methods showed that the classification of O. rubicundus into Odontamblyopus due to morphology is debatable.

Heteroplasmy and ancient translocation of mitochondrial DNA to the nucleus in the Chinese Horseshoe Bat (Rhinolophus sinicus) complex

The utility and reliability of mitochondrial DNA sequences in phylogenetic and phylogeographic studies may be compromised by widespread and undetected nuclear mitochondrial copies (numts) as well as heteroplasmy within individuals. Both numts and heteroplasmy are likely to be common across diverse taxa yet few studies have characterised their frequencies and variation at the intra-specific level. Here we report the presence of both numts and heteroplasmy in the mitochondrial control region of the Chinese horseshoe bat Rhinolophus sinicus. In total we generated 123 sequences from 18 bats, which contained two different numt clades (i.e. Numt-1 and Numt-2) and one mtDNA clade. The sequence divergence between Numt-1 and Numt-2 was 16.8% and each numt type was found in all four R. sinicus taxa, suggesting either two ancient translocations of mitochondrial DNA into the nucleus from the same source taxon, or a single translocation from different source taxa that occurred before the split of R. sinicus into different lineages. Within the mtDNA clade, phylogenetic relationships among the four taxa of R. sinicus were similar to those seen in previous results. Based on PCR comparisons, heteroplasmy was inferred between almost all individuals of R. sinicus with respect to sequence variation. Consistent with introgression of mtDNA between Central sinicus and septentrionalis, individuals from these two taxa exhibited similar signatures of repeated sequences in the control region. Our study highlights the importance of testing for the presence of numts and heteroplasmy when applying mtDNA markers to phylogenetic studies.

Mitochondrial DNA variants in obesity

Heritability estimates for body mass index (BMI) variation are high. For mothers and their offspring higher BMI correlations have been described than for fathers. Variation(s) in the exclusively maternally inherited mitochondrial DNA (mtDNA) might contribute to this parental effect. Thirty-two to 40 mtDNA single nucleotide polymorphisms (SNPs) were available from genome-wide association study SNP arrays (Affymetrix 6.0). For discovery, we analyzed association in a case-control (CC) sample of 1,158 extremely obese children and adolescents and 435 lean adult controls. For independent confirmation, 7,014 population-based adults were analyzed as CC sample of n = 1,697 obese cases (BMI ≥ 30 kg/m2) and n = 2,373 normal weight and lean controls (BMI<25 kg/m2). SNPs were analyzed as single SNPs and haplogroups determined by HaploGrep. Fisher's two-sided exact test was used for association testing. Moreover, the D-loop was re-sequenced (Sanger) in 192 extremely obese children and adolescents and 192 lean adult controls. Association testing of detected variants was performed using Fisher's two-sided exact test. For discovery, nominal association with obesity was found for the frequent allele G of m.8994G/A (rs28358887, p = 0.002) located in ATP6. Haplogroup W was nominally overrepresented in the controls (p = 0.039). These findings could not be confirmed independently. For two of the 252 identified D-loop variants nominal association was detected (m.16292C/T, p = 0.007, m.16189T/C, p = 0.048). Only eight controls carried the m.16292T allele, five of whom belonged to haplogroup W that was initially enriched among these controls. m.16189T/C might create an uninterrupted poly-C tract located near a regulatory element involved in replication of mtDNA. Though follow-up of some D-loop variants still is conceivable, our hypothesis of a contribution of variation in the exclusively maternally inherited mtDNA to the observed larger correlations for BMI between mothers and their offspring could not be substantiated by the findings of the present study.

Complete mitochondrial genome of the orange-spotted spinefoot Siganus guttatus (Perciformes, Siganidae)

The complete mitochondrial genome of the orange-spotted spinefoot Siganus guttatus (Perciformes, Siganidae) is presented in this paper. The entire sequence, 16,505 bp in length, is comprised of 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes and 1 control region. Structurally, the composition and order of genes are same as those found in other rabbitfishes. Except ND6 and 8 tRNA, all other mitochondrial genes were encoded on the heavy strand. Overall base compositions of mitogenome are 29.3% of A, 29.3% of C, 25.7% of T, and 15.7% of G, showing an obvious anti-G bias which was commonly found in fishes. Three conserved sequence blocks (CSB-1, CSB-2 and CSB-3) and long series of di- and mono-nucleotide microsatellite repeats in 3' end are observed in the putative control region. The complete mitogenome sequence of S. guttatus could provide useful information for further phylogenetic and population genetic studies.

In D-loop: 40 years of mitochondrial 7S DNA

Given the tiny size of the mammalian mitochondrial genome, at only 16.5 kb, it is often surprising how little we know about some of its molecular features, and the molecular mechanisms governing its maintenance. One such conundrum is the biogenesis and function of the mitochondrial displacement loop (D-loop). The mitochondrial D-loop is a triple-stranded region found in the major non-coding region (NCR) of many mitochondrial genomes, and is formed by stable incorporation of a third, short DNA strand known as 7S DNA. In this article we review the current affairs regarding the main features of the D-loop structure, the diverse frequency of D-loops in the mtDNAs of various species and tissues, and also the mechanisms of its synthesis and turnover. This is followed by an account of the possible functions of the mitochondrial D-loop that have been proposed over the last four decades. In the last section, we discuss the potential links of the D-loop with mammalian ageing.

Phylogeography and postglacial recolonization of Europe by Rhinolophus hipposideros: evidence from multiple genetic markers

The demographic history of Rhinolophus hipposideros (lesser horseshoe bat) was reconstructed across its European, North African and Middle-Eastern distribution prior to, during and following the most recent glaciations by generating and analysing a multimarker data set. This data set consisted of an X-linked nuclear intron (Bgn; 543 bp), mitochondrial DNA (cytb-tRNA-control region; 1630 bp) and eight variable microsatellite loci for up to 373 individuals from 86 localities. Using this data set of diverse markers, it was possible to determine the species' demography at three temporal stages. Nuclear intron data revealed early colonization into Europe from the east, which pre-dates the Quaternary glaciations. The mtDNA data supported multiple glacial refugia across the Mediterranean, the largest of which were found in the Ibero-Maghreb region and an eastern location (Anatolia/Middle East)-that were used by R. hipposideros during the most recent glacial cycles. Finally, microsatellites provided the most recent information on these species' movements since the Last Glacial Maximum and suggested that lineages that had diverged into glacial refugia, such as in the Ibero-Maghreb region, have remained isolated. These findings should be used to inform future conservation management strategies for R. hipposideros and show the power of using a multimarker data set for phylogeographic studies.

Comparative analysis of mitochondrial genomes in distinct nuclear ploidy loach Misgurnus anguillicaudatus and its implications for polyploidy evolution

Misgurnus anguillicaudatus has several natural ploidy types. To investigate whether nuclear polyploidy have an impact on mitochondrial DNA (mtDNA), the complete mitochondrial genomes (mitogenomes) of five distinct ploidy M. anguillicaudatus (natural diploid, triploid, tetraploid, pentaploid and hexaploid), which were collected in central China, were sequenced and analyzed. The five mitogenomes share the same gene arrangement and have similar gene size, base composition and codon usage pattern. The most variable regions of the mitogenome were the protein-coding genes, especially the ND4L (5.39% mutation rate). Most variations occurred in tetraploids. The phylogenetic tree showed that the tetraploid M. anguillicaudatus separated early from other ploidy loaches. Meanwhile, the mitogenomes from pentaploids, and hexaploids have the closest phylogenetic relations, but far from that of tetraploids, implying that pentaploids and hexaploids could not be formed from tetraploids, possibly from the diploids and triploids. The genus Misgurnus species were divided into two divergent inter-genus clades, and the five ploidy M. anguillicaudatus were monophyletic, which support the hypotheses about the mitochondrial introgression in loach species.

Complete mitochondrial genome of the Bloch's gizzard shad Nematalosa nasus (Clupeiformes: Clupeidae)

The complete mitochondrial genome sequence of Nematalosa nasus was determined in this paper. The entire sequence was 16,674 bp in length and the overall base compositions were 26.5% of A, 28.4% of C, 25.4% of T and 19.7% of G. The mitogenome contained 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes and two non-coding regions (the control region and the light strand replication origin). Gene arrangement was similar to that observed in most other vertebrates. Except for the ND6 genes and eight tRNA genes, all other genes were encoded on the heavy strand. Within the control region, the extended termination associated sequence was identified, as well as four central conserved sequence blocks (CSB-F, CSB-E, CSB-D and CSB-C) and three conserved sequence blocks (CSB-1, CSB-2 and CSB-3).

Complete mitochondrial genome of the mudskipper Boleophthalmus boddarti (Perciformes, Gobiidae)

The Boddart's goggle-eyed mudskipper, Boleophthalmus boddarti (Perciformes, Gobiidae) is an amphibious fish, inhabiting brackish waters of estuaries and builds burrows in soft mud along the intertidal zone. In this paper, the complete mitochondrial genome sequence of B. boddarti was firstly determined. The circle genome (16,727 bp) comprises 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes and 1 control region. The overall base composition of B. boddarti is 29.1% for C, 28.9% for A, 25.9% for T, and 16.0% for G, with a slight A + T bias of 54.8%. The termination-associated sequence, conserved sequence block domains, and a 131-bp tandem repeat were found in the control region. It has the typical vertebrate mitochondrial gene arrangement.

Genetic diversity analysis of mitochondrial DNA control region in artificially propagated Chinese sucker Myxocyprinus asiaticus

The genetic diversity of the three major artificially propagated populations of Chinese sucker, an endangered freshwater fish species, was investigated using the sequences of mitochondrial DNA (mtDNA) control regions. Among the 89 individuals tested, 66 variable sites (7.26%) and 10 haplotypes were detected (Haplotype diversity Hd = 0.805, Nucleotide diversity π = 0.0287). In general, genetic diversity was lower in artificially propagated populations than in wild populations. This reduction in genetic diversity may be due to population bottlenecks, genetic drift and human selection. A stepping-stone pattern of gene flow was detected in the populations studied, showing much higher gene flow between neighbouring populations. To increase the genetic diversity, wild lineages should be introduced, and more lineages should be shared among artificially propagated populations.

The complete mitochondrial genome of Cuon alpinus lepturus

The complete mitochondrial genome of a female Chinese dhole (Cuon alpinus lepturus) was sequenced for the first time. It contains 13 protein-coding genes, 2 rRNA genes, 22 tRNA genes and 1 control region.

The control region of mitochondrial DNA shows an unusual CpG and non-CpG methylation pattern

DNA methylation is a common epigenetic modification of the mammalian genome. Conflicting data regarding the possible presence of methylated cytosines within mitochondrial DNA (mtDNA) have been reported. To clarify this point, we analysed the methylation status of mtDNA control region (D-loop) on human and murine DNA samples from blood and cultured cells by bisulphite sequencing and methylated/hydroxymethylated DNA immunoprecipitation assays. We found methylated and hydroxymethylated cytosines in the L-strand of all samples analysed. MtDNA methylation particularly occurs within non-C-phosphate-G (non-CpG) nucleotides, mainly in the promoter region of the heavy strand and in conserved sequence blocks, suggesting its involvement in regulating mtDNA replication and/or transcription. We observed DNA methyltransferases within the mitochondria, but the inactivation of Dnmt1, Dnmt3a, and Dnmt3b in mouse embryonic stem (ES) cells results in a reduction of the CpG methylation, while the non-CpG methylation shows to be not affected. This suggests that D-loop epigenetic modification is only partially established by these enzymes. Our data show that DNA methylation occurs in the mtDNA control region of mammals, not only at symmetrical CpG dinucleotides, typical of nuclear genome, but in a peculiar non-CpG pattern previously reported for plants and fungi. The molecular mechanisms responsible for this pattern remain an open question.

Mitochondrial DNA rearrangements in health and disease--a comprehensive study

Mitochondrial DNA (mtDNA) rearrangements cause a wide variety of highly debilitating and often fatal disorders and have been implicated in aging and age-associated disease. Here, we present a meta-analytical study of mtDNA deletions (n = 730) and partial duplications (n = 37) using information from more than 300 studies published over the last 30 years. We show that both classes of mtDNA rearrangements are unequally distributed among disorders and their breakpoints have different genomic locations. We also demonstrate that 100% of cases with sporadic mtDNA deletions and 97.3% with duplications have no breakpoints in the 16,071 breakage hotspot site, in contrast with deletions from healthy and aged tissues. Notably, most deletions removing a section of the D-loop are found in tumors. Deleted mtDNA molecules lacking the origin of L-strand replication (O(L)) represent only 9.5% of all reported cases, whereas extra origins of replication occur in all duplications. As previously shown for deletions, imperfect stretches of homology are common in duplication breakpoints. Finally, we provide a dedicated Website with detailed information on deleted/duplicated mtDNA regions to facilitate the design of efficient methods for identification and screening of rearranged mitochondrial genomes (available at http://www.portugene.com/mtDNArearrangements.html).

Genetic variation and phylogenetic relationship between three serow species of the genus Capricornis based on the complete mitochondrial DNA control region sequences

The molecular evidence of phylogenetic status regarding the Formosan serow (Capricornis swinhoei) is not robust and little is known about the genetic diversity of the Sumatran serow (Capricornis sumatraensis), which partly is due to the hardness in sample collection. Here we determined the sequences of the complete mitochondrial DNA control region (1,014 bp) of 19 Sumatran-serow individuals. Nine new haplotypes were defined based on 78 variable sites. Combined analysis with other 32 haplotypes downloaded from the public database, including 1 Sumatran-serow, 11 Formosan-serow and 20 Japanese-serow (Capricornis crispus) haplotypes, a relatively high level of nucleotide diversity was first observed in Sumatran serow (π = 0.0249). By comparative analysis with structural consensus sequences from other mammals, we have identified central, left and right domains and depicted the putative functional structure, including extend termination associated sequences and conserve sequence blocks, in mtDNA control region. The alignment of mtDNA control region revealed that both Sumatran and Japanese serow have two tandem repeats (TRs), but three TRs in Formosan serow. Phylogenetic analyses revealed that the Formosan serow is distinct species with the Japanese serow, but a sister group with the Sumatran serow. The divergence time estimated among three serow species revealed that Pleistocene climate changes and the uplift of Qinghai-Tibetan plateau might play an important role in the genetic differentiation of the serows. These results mainly provide the convinced evidence on the genetic relationship between three serow species.

Mitochondrial genome sequence of Sibynophis chinensis (Squamata, Colubridae)

The complete mitochondrial (mt) genome of the Chinese many-toothed snake, Sibynophis chinensis, was sequenced and found to be 17,163 bp in length. The arrangement of 13 protein-coding genes, tRNAs and rRNAs was identical to that of other common snake mt genomes. The mt protein-coding genes of S. chinensis utilized ATA, ATG, ATA and GTG as initiation codons and AGA, AGG, TAA, TAG and T as termination codons. Among three tRNA clusters (LQM, WANCY and HSL), LQM was found instead of IQM, which is common in other vertebrates. We also identified two control regions that contained several conserved elements known as conserved sequence blocks and termination-associated sequences related to mt replication and transcription.