Strand-specific nucleotide composition bias in echinoderm and vertebrate mitochondrial genomes

Characterization and Phylogenetic Analysis of the Complete Mitochondrial Genome of Triplophysa microphthalma

The complete mitochondrial genome has been extensively utilized in studies related to phylogenetics, offering valuable perspectives on evolutionary relationships. The mitochondrial genome of the fine-eyed plateau loach, Triplophysa microphthalma, has not attracted much attention, although this species is endemic to China. In this study, we characterized the mitochondrial genome of T. microphthalma and reassessed the classification status of its genus. The complete mitochondrial genome of T. microphthalma was 16,591 bp and contained thirty-seven genes, including thirteen protein-coding genes (PCGs), two ribosomal RNA genes (rRNAs), and twenty-two transfer RNA genes (tRNAs). All but one of the thirteen PCGs had the regular start codon ATG; the gene cox1 started with GTG. Six PCGs had incomplete stop codons (T--). These thirteen PCGs are thought to have evolved under purifying selection, and the mitogenome shared a high degree of similarity with the genomes of species within the genus Leptobotia. All tRNA genes exhibited the standard clover-shaped structure, with the exception of the trnS1 gene, which lacked a DHU stem. A phylogenetic analysis indicated that T. microphthalma was more closely related to species within the genus Triplophysa than to those in Barbatula. The present study contributes valuable genomic information for T. microphthalma, and offers new perspectives on the phylogenetic relationships among species of Triplophysa and Barbatula. The findings also provide essential data that can inform the management and conservation strategies for T. microphthalma and other species of Triplophysa and Barbatula.

Characterization of the Complete Mitochondrial Genome of Schizothorax kozlovi (Cypriniformes, Cyprinidae, Schizothorax) and Insights into the Phylogenetic Relationships of Schizothorax

Schizothorax kozlovi is an endemic and vulnerable fish species found in the upper Yangtze River in China. Over the past few years, the population resources of S. kozlovi have been nearly completely depleted owing to multiple contributing threats. While the complete mitochondrial genomes serve as important molecular markers for phylogenetic and genetic studies, the mitochondrial genome of S. kozlovi has still received little attention. In this study, we analyzed the characterization of the mitochondrial genome of S. kozlovi and investigated the phylogenetic relationships of Schizothorax. The complete mitochondrial genome of S. kozlovi was 16,585 bp in length, which contained thirty-seven genes (thirteen protein-coding genes (PCGs), two ribosomal RNA genes (rRNAs), twenty-two transfer RNA genes (tRNAs)) and two non-coding regions for the origin of light strand (OL) and the control region (CR). There were nine overlapping regions and seventeen intergenic spacers regions in the mitochondrial genome. The genome also showed a bias towards A + T content (55.01%) and had a positive AT-skew (0.08) and a negative GC-skew (-0.20). All the PCGs employed the ATG or GTG as the start codon and TAA, TAG, or single T as the stop codon. Additionally, all of the tRNAs displayed a typical cloverleaf secondary structure, except trnS1 which lacked the D arm. The phylogenetic analysis, based on the maximum likelihood (ML) and Bayesian inference (BI) methods, revealed that the topologies of the phylogenetic tree divided the Schizothorax into four clades and did not support the classification of Schizothorax based on morphology. The phylogenetic status of S. kozlovi was closely related to that of S. chongi. The present study provides valuable genomic information for S. kozlovi and new insights in phylogenetic relationships of Schizothorax. These data could also offer fundamental references and guidelines for the management and conservation of S. kozlovi and other species of Schizothorax.

Novel mitochondrial genome rearrangements including duplications and extensive heteroplasmy could underlie temperature adaptations in Antarctic notothenioid fishes

Mitochondrial genomes are known for their compact size and conserved gene order, however, recent studies employing long-read sequencing technologies have revealed the presence of atypical mitogenomes in some species. In this study, we assembled and annotated the mitogenomes of five Antarctic notothenioids, including four icefishes (Champsocephalus gunnari, C. esox, Chaenocephalus aceratus, and Pseudochaenichthys georgianus) and the cold-specialized Trematomus borchgrevinki. Antarctic notothenioids are known to harbor some rearrangements in their mt genomes, however the extensive duplications in icefishes observed in our study have never been reported before. In the icefishes, we observed duplications of the protein coding gene ND6, two transfer RNAs, and the control region with different copy number variants present within the same individuals and with some ND6 duplications appearing to follow the canonical Duplication-Degeneration-Complementation (DDC) model in C. esox and C. gunnari. In addition, using long-read sequencing and k-mer analysis, we were able to detect extensive heteroplasmy in C. aceratus and C. esox. We also observed a large inversion in the mitogenome of T. borchgrevinki, along with the presence of tandem repeats in its control region. This study is the first in using long-read sequencing to assemble and identify structural variants and heteroplasmy in notothenioid mitogenomes and signifies the importance of long-reads in resolving complex mitochondrial architectures. Identification of such wide-ranging structural variants in the mitogenomes of these fishes could provide insight into the genetic basis of the atypical icefish mitochondrial physiology and more generally may provide insights about their potential role in cold adaptation.

Mitogenomic Codon Usage Patterns of Superfamily Certhioidea (Aves, Passeriformes): Insights into Asymmetrical Bias and Phylogenetic Implications

The superfamily Certhioidea currently comprises five families. Due to the rapid diversification, the phylogeny of Certhioidea is still controversial. The advent of next generation sequencing provides a unique opportunity for a mitogenome-wide study. Here, we first provided six new complete mitogenomes of Certhioidea (Certhia americana, C. familiaris, Salpornis spilonota, Cantorchilus leucotis, Pheugopedius coraya, and Pheugopedius genibarbis). We further paid attention to the genomic characteristics, codon usages, evolutionary rates, and phylogeny of the Certhioidea mitogenomes. All mitogenomes we analyzed displayed typical ancestral avian gene order with 13 protein-coding genes (PCGs), 22 tRNAs, 2 rRNAs, and one control region (CR). Our study indicated the strand-biased compositional asymmetry might shape codon usage preferences in mitochondrial genes. In addition, natural selection might be the main factor in shaping the codon usages of genes. Additionally, evolutionary rate analyses indicated all mitochondrial genes were under purifying selection. Moreover, MT-ATP8 and MT-CO1 were the most rapidly evolving gene and conserved genes, respectively. According to our mitophylogenetic analyses, the monophylies of Troglodytidae and Sittidae were strongly supported. Importantly, we suggest that Salpornis should be separated from Certhiidae and put into Salpornithidae to maintain the monophyly of Certhiidae. Our findings are useful for further evolutionary studies within Certhioidea.

Legacies of domestication, Neolithic diffusion and trade between Indian subcontinent and Island Southeast Asia shape maternal genetic diversity of Andaman cattle

Andaman cattle is a precious indigenous livestock species endemic to Andaman and Nicobar Islands, India. Till date, origin and genetic makeup of the breed which is warranted for breed conservation is not known. Moreover, the spread of zebu cattle from Indus valley to different parts of Island Southeast Asia (ISEA) is not properly understood. Here, we report the genetic diversity, population structure of Andaman cattle and their evolution in the context of epicentre of zebu domestication and ISEA. High genetic diversity in complete mitochondrial D-loop sequences indicated the ability of the breed to withstand impending climate change. Total 81 haplotypes were detected and all of them except three belonged to Bos indicus. The presence of taurine haplotypes in Andaman cattle indicate introgression by European-derived cattle. A poor phylogenetic signal of Andaman cattle with genetic affinities with cattle of Indian subcontinent and ISEA was observed. The poor phylogenetic structure may be due to multidirectional gene flow from Indian subcontinent and ISEA, with which Andaman shares a close cultural and trade relationship from Neolithic age. We hypothesize that Andaman cattle is the outcome of Neolithic diffusion from centre of zebu domestication along with multidirectional commercial exchange between Indian subcontinent and ISEA.

A replication-linked mutational gradient drives somatic mutation accumulation and influences germline polymorphisms and genome composition in mitochondrial DNA

Mutations in mitochondrial DNA (mtDNA) cause maternally inherited diseases, while somatic mutations are linked to common diseases of aging. Although mtDNA mutations impact health, the processes that give rise to them are under considerable debate. To investigate the mechanism by which de novo mutations arise, we analyzed the distribution of naturally occurring somatic mutations across the mouse and human mtDNA obtained by Duplex Sequencing. We observe distinct mutational gradients in G→A and T→C transitions delimited by the light-strand origin and the mitochondrial Control Region (mCR). The gradient increases unequally across the mtDNA with age and is lost in the absence of DNA polymerase γ proofreading activity. In addition, high-resolution analysis of the mCR shows that important regulatory elements exhibit considerable variability in mutation frequency, consistent with them being mutational ‘hot-spots’ or ‘cold-spots’. Collectively, these patterns support genome replication via a deamination prone asymmetric strand-displacement mechanism as the fundamental driver of mutagenesis in mammalian DNA. Moreover, the distribution of mtDNA single nucleotide polymorphisms in humans and the distribution of bases in the mtDNA across vertebrate species mirror this gradient, indicating that replication-linked mutations are likely the primary source of inherited polymorphisms that, over evolutionary timescales, influences genome composition during speciation.

Understanding the codon usage patterns of mitochondrial CO genes among Amphibians

A universal phenomenon of using synonymous codons unequally in coding sequences known as codon usage bias (CUB) is observed in all forms of life. Mutation and natural selection drive CUB in many species but the relative role of evolutionary forces varies across species, genes and genomes. We studied the CUB in mitochondrial (mt) CO genes from three orders of Amphibia using bioinformatics approach as no work was reported yet. We observed that CUB of mt CO genes of Amphibians was weak across different orders. Order Caudata had higher CUB followed by Gymnophiona and Anura for all genes and CUB also varied across genes. Nucleotide composition analysis showed that CO genes were AT-rich. The AT content in Caudata was higher than that in Gymnophiona while Anura showed the least content. Multiple investigations namely nucleotide composition, correspondence analysis, parity plot analysis showed that the interplay of mutation pressure and natural selection caused CUB in these genes. Neutrality plot suggested the involvement of natural selection was more than the mutation pressure. The contribution of natural selection was higher in Anura than Gymnophiona and the lowest in Caudata. The codons CGA, TGA, AAA were found to be highly favoured by nature across all genes and orders.

Comprehensive whole genome survey analyses of male and female brown-spotted flathead fish Platycephalus sp.1

The flathead fish Platycephalus sp.1 is an ecologically and commercially important marine fish in the northwestern Pacific with notable sexual differences in growth and development. Yet the genomic data of this species is lacking. In the present study, whole genome sequencing of two individuals (one male and one female) of Platycephalus sp.1 were conducted to provide fundamental genomic information. The genome sizes were estimated to be 674.96 Mb (male) and 684.15 Mb (female) by using k-mer analyses. The heterozygosity and repeat ratios suggested possible male heterogamety of this species. The draft genome sequences were initially assembled and genome-wide microsatellite motifs were identified. Besides, the complete mitochondrial genome sequences were assembled and the phylogenetic analyses genetically supported the validation of Platycephalus sp.1. The reported genomic data and genetic markers in this study could be useful in future comparative genomics and evolutionary biology studies.

Analysis of compositional properties and codon usage bias of mitochondrial CYB gene in anura, urodela and gymnophiona

We delineated the pattern of synonymous codon usage bias (CUB) and its determinants in mitochondrial CYB gene of respiratory chain across different amphibian groups namely orders anura, urodela and gymnophiona. We observed that CUB was low in CYB gene of amphibia. The gymnophionans had comparatively high bias followed by urodeles and anurans. The codons namely TCA, CCA, CAA, CGA, TGA, AAA and ACA were over-represented in all three orders. The codons such as GCC and TCC were over-represented in anura whereas in urodela, the over-represented codons were TTA, CTA, ATA, GTA, GAA, GGA and GCA. In gymnophiona, GCC, TTA, CTA, ATA, GTA, GAA and GGA codons were over-represented. The regression analysis between effective number of codons (ENC) and nucleobase at the 3rd position revealed that nucleobase A and C influenced CUB positively in order anura, while in urodela and gymnophiona, nucleobase A and T influenced the CUB positively. Mutation pressure and natural selection mutually illustrate the CUB of CYB gene (complex III gene) of amphibia as elucidated by correlation analysis between 3rd nucleotide in a codon and overall nucleotide content of the gene. However, neutrality plot showed that natural selection was the dominant evolutionary factor of CUB.

The First Complete Mitochondrial Genome of the Flathead Cociella crocodilus (Scorpaeniformes: Platycephalidae) and the Phylogenetic Relationships within Scorpaeniformes Based on Whole Mitogenomes

Complete mitochondrial genomes (mitogenomes) are important molecular markers for understanding the phylogenetics of various species. Although recent studies on the mitogenomes of the Scorpaeniformes species have been greatly advanced, information regarding molecular studies and the taxonomic localization of Platycephalidae is still sparse. To further analyze the phylogeny of Platycephalidae, we sequenced the complete mitogenome of Cociella crocodilus of the Platycephalidae family within Scorpaeniformes for the first time. The mitogenome was 17,314 bp in length, contained two ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), 13 protein-coding genes (PCGs), and two typical noncoding control regions (the control region (CR) and origin of the light strand (OL)). All PCGs used standard initiation codons ATG, apart from cox1. The majority of the tRNA genes could be folded into cloverleaf secondary structures, whereas the secondary structure of tRNASer (AGN) lacked a dihydrouridine (DHU) stem. The CR contained several conserved sequence blocks (CSBs) and eight tandem repeats. In addition, the phylogenetic relationship based on the concatenated nucleotides sequences of 13 PCGs indicated that the Platycephalidae species are relatively basal in the phylogenetic relationships of Scorpaeniformes. Our results may not only advance the origin and the evolution of Scorpaeniformes, but also provide information for the genetic evolution and taxonomy of the teleostean species.

The complete mitochondrial genomes of Tarsiger cyanurus and Phoenicurus auroreus: a phylogenetic analysis of Passeriformes

Passeriformes is the largest group within aves and the phylogenetic relationships between Passeriformes have caused major disagreement in ornithology. Particularly, the phylogenetic relationships between muscicapoidea and sylvioidea are complex, and their taxonomic boundaries have not been clearly defined. Our aim was to study the status of two bird species: Tarsiger cyanurus and Phoenicurus auroreus. Furthermore, we analyzed the phylogenetic relationships of Passeriformes. Complete mitochondrial DNA (mtDNA) sequences of both species were determined and the lengths were 16,803 (T. cyanurus) and 16,772 bp (P. auroreus), respectively. Thirteen protein-coding genes, 22 tRNA genes, two rRNA genes, and one control region were identified in these mtDNAs. The contents of A and T at the base compositions was significantly higher than the content of G and C, and this AT skew was positive, while the GC skew was negative. The monophyly of Passeriformes is divided into four major clades: Corvoidea, Sylvioidea, Passeroidea, and Musicicapoidea. Paridae should be separated from the superfamily Sylvioidea and placed within the superfamily Muscicapoidea. The family Muscicapidae and Corvida were paraphyly, while Carduelis and Emberiza were grouped as a sister taxon. The relationships between some species of the order passeriformes may remain difficult to resolve despite an effort to collect additional characters for phylogenetic analysis. Current research of avian phylogeny should focus on adding characters and taxa and use both effectively to obtain a better resolution for deeper and shallow nodes.

Multiple independent structural dynamic events in the evolution of snake mitochondrial genomes

BACKGROUND

Mitochondrial DNA sequences have long been used in phylogenetic studies. However, little attention has been paid to the changes in gene arrangement patterns in the snake's mitogenome. Here, we analyzed the complete mitogenome sequences and structures of 65 snake species from 14 families and examined their structural patterns, organization and evolution. Our purpose was to further investigate the evolutionary implications and possible rearrangement mechanisms of the mitogenome within snakes.

RESULTS

In total, eleven types of mitochondrial gene arrangement patterns were detected (Type I, II, III, III-A, III-B, III-B1, III-C, III-D, III-E, III-F, III-G), with mitochondrial genome rearrangements being a major trend in snakes, especially in Alethinophidia. In snake mitogenomes, the rearrangements mainly involved three processes, gene loss, translocation and duplication. Within Scolecophidia, the O_L was lost several times in Typhlopidae and Leptotyphlopidae, but persisted as a plesiomorphy in the Alethinophidia. Duplication of the control region and translocation of the tRNA^Leu gene are two visible features in Alethinophidian mitochondrial genomes. Independently and stochastically, the duplication of pseudo-Pro (P*) emerged in seven different lineages of unequal size in three families, indicating that the presence of P* was a polytopic event in the mitogenome.

CONCLUSIONS

The WANCY tRNA gene cluster and the control regions and their adjacent segments were hotspots for mitogenome rearrangement. Maintenance of duplicate control regions may be the source for snake mitogenome structural diversity.

Mitochondrial genomes of three Tetrigoidea species and phylogeny of Tetrigoidea

The mitochondrial genomes (mitogenomes) of Formosatettix qinlingensis, Coptotettix longjiangensis and Thoradonta obtusilobata (Orthoptera: Caelifera: Tetrigoidea) were sequenced in this study, and almost the entire mitogenomes of these species were determined. The mitogenome sequences obtained for the three species were 15,180, 14,495 and 14,538 bp in length, respectively, and each sequence included 13 protein-coding genes (PCGs), partial sequences of rRNA genes (rRNAs), tRNA genes (tRNAs) and a A + T-rich region. The order and orientation of the gene arrangement pattern were identical to that of most Tetrigoidea species. Some conserved spacer sequences between trnS(UCN) and nad1 were useful to identify Tetrigoidea and Acridoidea. The Ka/Ks value of atp8 between Trachytettix bufo and other four Tetrigoidea species indicated that some varied sites in this gene might be related with the evolution of T. bufo. The three Tetrigoidea species were compared with other Caelifera. At the superfamily level, conserved sequences were observed in intergenic spacers, which can be used for superfamily level identification between Tetrigoidea and Acridoidea. Furthermore, a phylogenomic analysis was conducted based on the concatenated data sets from mitogenome sequences of 24 species of Orthoptera in the superorders Caelifera and Ensifera. Both maximum likelihood and bayesian inference analyses strongly supported Acridoidea and Tetrigoidea as forming monophyletic groups. The relationships among six Tetrigoidea species were (((((Tetrix japonica, Alulatettix yunnanensis), Formosatettix qinlingensis), Coptotettix longjiangensis), Trachytettix bufo), Thoradonta obtusilobata).

Cellular and population level processes influence the rate, accumulation and observed frequency of inherited and somatic mtDNA mutations

Mitochondria are found in all animals and have the unique feature of containing multiple copies of their own small, circular DNA genome (mtDNA). The rate and pattern of mutation accumulation in the mtDNA are influenced by molecular, cellular and population level processes. We distinguish between inherited and somatic mtDNA mutations and review evidence for the often-made assumption that mutations accumulate at a higher rate in mtDNA than in nuclear DNA (nDNA). We conclude that the whole genome mutation accumulation rate is higher for mtDNA than for nDNA but include the caveat that rates overlap considerably between the individual mtDNA- and nDNA-encoded genes. Next, we discuss the postulated causal mechanisms for the high rate of mtDNA mutation accumulation in both inheritance and in somatic cells. Perhaps unexpectedly, mtDNA is resilient to many mutagens of nDNA but is prone to errors of replication. We then consider the influence of maternal inheritance, recombination and selection on the observed accumulation pattern of inherited mtDNA mutations. Finally, we discuss environmental influences of temperature and diet on the observed frequency of inherited and somatic mtDNA mutations. We conclude that it is necessary to understand the cellular processes to fully interpret the pattern of mutations and how they influence our interpretations of evolution and disease.

The mitochondrial genome sequences of the round goby and the sand goby reveal patterns of recent evolution in gobiid fish

Background

Vertebrate mitochondrial genomes are optimized for fast replication and low cost of RNA expression. Accordingly, they are devoid of introns, are transcribed as polycistrons and contain very little intergenic sequences. Usually, vertebrate mitochondrial genomes measure between 16.5 and 17 kilobases (kb).

Results

During genome sequencing projects for two novel vertebrate models, the invasive round goby and the sand goby, we found that the sand goby genome is exceptionally small (16.4 kb), while the mitochondrial genome of the round goby is much larger than expected for a vertebrate. It is 19 kb in size and is thus one of the largest fish and even vertebrate mitochondrial genomes known to date. The expansion is attributable to a sequence insertion downstream of the putative transcriptional start site. This insertion carries traces of repeats from the control region, but is mostly novel. To get more information about this phenomenon, we gathered all available mitochondrial genomes of Gobiidae and of nine gobioid species, performed phylogenetic analyses, analysed gene arrangements, and compared gobiid mitochondrial genome sizes, ecological information and other species characteristics with respect to the mitochondrial phylogeny. This allowed us amongst others to identify a unique arrangement of tRNAs among Ponto-Caspian gobies.

Conclusions

Our results indicate that the round goby mitochondrial genome may contain novel features. Since mitochondrial genome organisation is tightly linked to energy metabolism, these features may be linked to its invasion success. Also, the unique tRNA arrangement among Ponto-Caspian gobies may be helpful in studying the evolution of this highly adaptive and invasive species group. Finally, we find that the phylogeny of gobiids can be further refined by the use of longer stretches of linked DNA sequence.

Electronic supplementary material

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

Two complete mitochondrial genomes from Praticolella mexicana Perez, 2011 (Polygyridae) and gene order evolution in Helicoidea (Mollusca, Gastropoda)

Helicoidea is a diverse group of land snails with a global distribution. While much is known regarding the relationships of helicoid taxa, comparatively little is known about the evolution of the mitochondrial genome in the superfamily. We sequenced two complete mitochondrial genomes from Praticolellamexicana Perez, 2011 representing the first such data from the helicoid family Polygyridae, and used them in an evolutionary analysis of mitogenomic gene order. We found the mitochondrial genome of Praticolellamexicana to be 14,008 bp in size, possessing the typical 37 metazoan genes. Multiple alternate stop codons are used, as are incomplete stop codons. Mitogenome size and nucleotide content is consistent with other helicoid species. Our analysis of gene order suggested that Helicoidea has undergone four mitochondrial rearrangements in the past. Two rearrangements were limited to tRNA genes only, and two involved protein coding genes.

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.

Multiple Factors Drive Replicating Strand Composition Bias in Bacterial Genomes

Composition bias from Chargaff’s second parity rule (PR2) has long been found in sequenced genomes, and is believed to relate strongly with the replication process in microbial genomes. However, some disagreement on the underlying reason for strand composition bias remains. We performed an integrative analysis of various genomic features that might influence composition bias using a large-scale dataset of 1111 genomes. Our results indicate (1) the bias was stronger in obligate intracellular bacteria than in other free-living species (p-value = 0.0305); (2) Fusobacteria and Firmicutes had the highest average bias among the 24 microbial phyla analyzed; (3) the strength of selected codon usage bias and generation times were not observably related to strand composition bias (p-value = 0.3247); (4) significant negative relationships were found between GC content, genome size, rearrangement frequency, Clusters of Orthologous Groups (COG) functional subcategories A, C, I, Q, and composition bias (p-values < 1.0 × 10⁻⁸); (5) gene density and COG functional subcategories D, F, J, L, and V were positively related with composition bias (p-value < 2.2 × 10⁻¹⁶); and (6) gene density made the most important contribution to composition bias, indicating transcriptional bias was associated strongly with strand composition bias. Therefore, strand composition bias was found to be influenced by multiple factors with varying weights.

Two nearly complete mitogenomes of wheat stem borers, Cephus pygmeus (L.) and Cephus sareptanus Dovnar-Zapolskij (Hymenoptera: Cephidae): an unusual elongation of rrnS gene

Two nearly complete mitochondrial genomes (mitogenomes) of wheat stem borers, Cephus pygmeus and Cephus sareptanus (Hymenoptera: Cephidae), were sequenced, characterised and compared with the previously known mitogenome of Cephus cinctus. The gene orders are mostly conserved, except for translocation of trnM and swapped position of trnI and trnQ. An A+T bias was found, but a deviation from strand asymmetry was also detected on the J strand. All protein coding genes (PCGs) are initiated by ATN codons, except for nad1, nad2 and atp8, and all are terminated with TAA, TA- or T- as a stop codon. The predicted secondary structures of rrnS and rrnL genes are mostly consistent with reported hymenopteran species. However, an unusual elongation in rrnS, not know elsewhere in the order, was discovered in Cephus species. Three autonomous sequences detected in domains I and II are mainly responsible for the length expansions.

In vivo occupancy of mitochondrial single-stranded DNA binding protein supports the strand displacement mode of DNA replication

Mitochondrial DNA (mtDNA) encodes for proteins required for oxidative phosphorylation, and mutations affecting the genome have been linked to a number of diseases as well as the natural ageing process in mammals. Human mtDNA is replicated by a molecular machinery that is distinct from the nuclear replisome, but there is still no consensus on the exact mode of mtDNA replication. We here demonstrate that the mitochondrial single-stranded DNA binding protein (mtSSB) directs origin specific initiation of mtDNA replication. MtSSB covers the parental heavy strand, which is displaced during mtDNA replication. MtSSB blocks primer synthesis on the displaced strand and restricts initiation of light-strand mtDNA synthesis to the specific origin of light-strand DNA synthesis (OriL). The in vivo occupancy profile of mtSSB displays a distinct pattern, with the highest levels of mtSSB close to the mitochondrial control region and with a gradual decline towards OriL. The pattern correlates with the replication products expected for the strand displacement mode of mtDNA synthesis, lending strong in vivo support for this debated model for mitochondrial DNA replication.

Mitochondria are cytoplasmatic organelles that produce most of the adenosine triphosphate (ATP) used by the cell as a source of chemical energy. A subset of proteins required for ATP production is encoded by a distinct mitochondrial DNA genome (mtDNA). Proper maintenance of mtDNA is essential, since mutations or depletion of this circular molecule may lead to a number of different diseases and also contribute to normal ageing. We are interested in the molecular mechanisms that ensure correct replication and propagation of mtDNA. Even if many of the responsible enzymes have been identified, there is still a debate within our scientific field regarding the exact mode of mtDNA replication. We have here used a combination of in vitro biochemistry and in vivo protein-DNA interaction characterization to address this question. Our findings demonstrate that the mitochondrial single-stranded DNA-binding protein (mtSSB) restricts initiation of mtDNA replication to a specific origin of replication. By characterizing how mtSSB interacts with the two strands of mtDNA in vivo, we are able to directly demonstrate the relevance of one proposed mode of mitochondrial DNA replication and at the same time seriously question the validity of other, alternative modes that have been proposed over the years.

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.

Complete mitochondrial genomes of five skippers (Lepidoptera: Hesperiidae) and phylogenetic reconstruction of Lepidoptera

We sequenced mitogenomes of five skippers (family Hesperiidae, Lepidoptera) to obtain further insight into the characteristics of butterfly mitogenomes and performed phylogenetic reconstruction using all available gene sequences (PCGs, rRNAs, and tRNAs) from 85 species (20 families in eight superfamilies). The general genomic features found in the butterflies also were found in the five skippers: a high A+T composition (79.3%-80.9%), dominant usage of TAA stop codon, similar skewness pattern in both strands, consistently length intergenic spacer sequence between tRNA(Gln) and ND2 (64-87 bp), conserved ATACTAA motif between tRNA(Ser (UCN)) and ND1, and characteristic features of the A+T-rich region (the ATAGA motif, varying length of poly-T stretch, and poly-A stretch). The start codon for COI was CGA in four skippers as typical, but Lobocla bifasciatus evidently possessed canonical ATG as start codon. All species had the ancestral arrangement tRNA(Asn)/tRNA(Ser (AGN)), instead of the rearrangement tRNA(Ser (AGN))/tRNA(Asn), found in another skipper species (Erynnis). Phylogenetic analyses using all available genes (PCGs, rRNAS, and tRNAs) yielded the consensus superfamilial relationships ((((((Bombycoidea+Noctuoidea+Geometroidea)+Pyraloidea)+Papilionoidea)+Tortricoidea)+Yponomeutoidea)+Hepialoidea), confirming the validity of Macroheterocera (Bombycoidea, Noctuoidea, and Geometroidea in this study) and its sister relationship to Pyraloidea. Within Rhopalocera (butterflies and skippers) the familial relationships (Papilionidae+(Hesperiidae+(Pieridae+((Lycaenidae+Riodinidae)+Nymphalidae)))) were strongly supported in all analyses (0.98-1 by BI and 96-100 by ML methods), rendering invalid the superfamily status for Hesperioidea. On the other hand, current mitogenome-based phylogeny did not find consistent superfamilial relationships among Noctuoidea, Geometroidea, and Bombycoidea and the familial relationships within Bombycoidea between analyses, requiring further taxon sampling in future studies.

Genetic Diversity of mtDNA D-loop and Maternal Origin of Three Chinese Native Horse Breeds

In order to protect the genetic resource of native horse breeds, the genetic diversity of mitochondrial DNA (mtDNA) D-loop of three native horse breeds in western China were investigated. Forty-three 600 bp mtDNA D-loop sequences were analyzed by PCR and sequencing techniques, 33 unique haplotypes with 70 polymorphic sites were detected in these horses, which account for 11.67% of 600 bp sequence analyzed, showing the abundant genetic diversity of the three native horse breeds in western China. The Neighbour-Joining (NJ) phylogenetic tree based on 247 bp of 43 D-loop sequences demonstrated the presence of seven major lineages (A to G), indicating that the three native horse breeds in western China originated from multiple maternal origins. Consistent with the front, the NJ phylogenetic tree based on 600 bp of mtDNA D-loop sequences of 43 Chinese western native horses and 81 sequences of six horse breeds from GenBank indicated that the three horse breeds had distributed into the seven major lineages (A to G). The structure of the phylogenic tree is often blurred because the variation in a short segment of the mitochondrial genome is often accompanied by high levels of recurrent mutations. Consequently, longer D-loop sequences are helpful in achieving a higher level of molecular resolution in horses.

Structure and evolution of the mitochondrial genome of Exorista sorbillans: the Tachinidae (Diptera: Calyptratae) perspective

The first complete mitochondrial genome (mitogenome) of Tachinidae Exorista sorbillans (Diptera) is sequenced by PCR-based approach. The circular mitogenome is 14,960 bp long and has the representative mitochondrial gene (mt gene) organization and order of Diptera. All protein-coding sequences are initiated with ATN codon; however, the only exception is Cox I gene, which has a 4-bp ATCG putative start codon. Ten of the thirteen protein-coding genes have a complete termination codon (TAA), but the rest are seated on the H strand with incomplete codons. The mitogenome of E. sorbillans is biased toward A+T content at 78.4 %, and the strand-specific bias is in reflection of the third codon positions of mt genes, and their T/C ratios as strand indictor are higher on the H strand more than those on the L strand pointing at any strain of seven Diptera flies. The length of the A+T-rich region of E. sorbillans is 106 bp, including a tandem triple copies of a13-bp fragment. Compared to Haematobia irritans, E. sorbillans holds distant relationship with Drosophila. Phylogenetic topologies based on the amino acid sequences, supporting that E. sorbillans (Tachinidae) is clustered with strains of Calliphoridae and Oestridae, and superfamily Oestroidea are polyphyletic groups with Muscidae in a clade.

A three-dimensional topology of complex I inferred from evolutionary correlations

Background

The quaternary structure of eukaryotic NADH:ubiquinone oxidoreductase (complex I), the largest complex of the oxidative phosphorylation, is still mostly unresolved. Furthermore, it is unknown where transiently bound assembly factors interact with complex I. We therefore asked whether the evolution of complex I contains information about its 3D topology and the binding positions of its assembly factors. We approached these questions by correlating the evolutionary rates of eukaryotic complex I subunits using the mirror-tree method and mapping the results into a 3D representation by multidimensional scaling.

Results

More than 60% of the evolutionary correlation among the conserved seven subunits of the complex I matrix arm can be explained by the physical distance between the subunits. The three-dimensional evolutionary model of the eukaryotic conserved matrix arm has a striking similarity to the matrix arm quaternary structure in the bacterium Thermus thermophilus (rmsd=19 Å) and supports the previous finding that in eukaryotes the N-module is turned relative to the Q-module when compared to bacteria. By contrast, the evolutionary rates contained little information about the structure of the membrane arm. A large evolutionary model of 45 subunits and assembly factors allows to predict subunit positions and interactions (rmsd = 52.6 Å). The model supports an interaction of NDUFAF3, C8orf38 and C2orf56 during the assembly of the proximal matrix arm and the membrane arm. The model further suggests a tight relationship between the assembly factor NUBPL and NDUFA2, which both have been linked to iron-sulfur cluster assembly, as well as between NDUFA12 and its paralog, the assembly factor NDUFAF2.

Conclusions

The physical distance between subunits of complex I is a major correlate of the rate of protein evolution in the complex I matrix arm and is sufficient to infer parts of the complex’s structure with high accuracy. The resulting evolutionary model predicts the positions of a number of subunits and assembly factors.

Base-pairing versatility determines wobble sites in tRNA anticodons of vertebrate mitogenomes

BACKGROUND

Vertebrate mitochondrial genomes typically have one transfer RNA (tRNA) for each synonymous codon family. This limited anticodon repertoire implies that each tRNA anticodon needs to wobble (establish a non-Watson-Crick base pairing between two nucleotides in RNA molecules) to recognize one or more synonymous codons. Different hypotheses have been proposed to explain the factors that determine the nucleotide composition of wobble sites in vertebrate mitochondrial tRNA anticodons. Until now, the two major postulates--the "codon-anticodon adaptation hypothesis" and the "wobble versatility hypothesis"--have not been formally tested in vertebrate mitochondria because both make the same predictions regarding the composition of anticodon wobble sites. The same is true for the more recent "wobble cost hypothesis".

PRINCIPAL FINDINGS

In this study we have analyzed the occurrence of synonymous codons and tRNA anticodon wobble sites in 1553 complete vertebrate mitochondrial genomes, focusing on three fish species with mtDNA codon usage bias reversal (L-strand is GT-rich). These mitogenomes constitute an excellent opportunity to study the evolution of the wobble nucleotide composition of tRNA anticodons because due to the reversal the predictions for the anticodon wobble sites differ between the existing hypotheses. We observed that none of the wobble sites of tRNA anticodons in these unusual mitochondrial genomes coevolved to match the new overall codon usage bias, suggesting that nucleotides at the wobble sites of tRNA anticodons in vertebrate mitochondrial genomes are determined by wobble versatility.

CONCLUSIONS/SIGNIFICANCE

Our results suggest that, at wobble sites of tRNA anticodons in vertebrate mitogenomes, selection favors the most versatile nucleotide in terms of wobble base-pairing stability and that wobble site composition is not influenced by codon usage. These results are in agreement with the "wobble versatility hypothesis".

The relation between the neutrality index for mitochondrial genes and the distribution of mutational effects on fitness

We explore factors affecting patterns of polymorphism and divergence (as captured by the neutrality index) at mammalian mitochondrial loci. To do this, we develop a population genetic model that incorporates a fraction of neutral amino acid sites, mutational bias, and a probability distribution of selection coefficients against new nonsynonymous mutations. We confirm, by reanalyzing publicly available datasets, that the mitochondrial cyt-b gene shows a broad range of neutrality indices across mammalian taxa, and explore the biological factors that can explain this observation. We find that observed patterns of differences in the neutrality index, polymorphism, and divergence are not caused by differences in mutational bias. They can, however, be explained by a combination of a small fraction of neutral amino acid sites, weak selection acting on most amino acid mutations, and differences in effective population size among taxa.

Current Status of Echinoderm Genome Analysis - What do we Know?

Echinoderms have long served as model organisms for a variety of biological research, especially in the field of developmental biology. Although the genome of the purple sea urchin Strongylocentrotus purpuratus has been sequenced, it is the only echinoderm whose whole genome sequence has been reported. Nevertheless, data is rapidly accumulating on the chromosomes and genomic sequences of all five classes of echinoderms, including the mitochondrial genomes and Hox genes. This blossoming new data will be essential for estimating the phylogenetic relationships among echinoderms, and also to examine the underlying mechanisms by which the diverse morphologies of echinoderms have arisen.

The first mitochondrial genome for the wasp superfamily Platygastroidea: the egg parasitoid Trissolcus basalis

The nearly complete mitochondrial (mt) genome of an egg parasitoid, Trissolcus basalis (Wollaston), was sequenced using both 454 and Illumina next-generation sequencing technologies. A portion of the noncoding region remained unsequenced, possibly owing to the presence of repeats. The sequenced portion of the genome is 15,768 bp and has a high A+T content (84.2%), as is typical for hymenopteran mt genomes. A total of 36 of the 37 genes normally present in animal mt genomes were located. The one exception was trnR; a truncated version of this gene is present between trnS(1) and nd5, but it is unclear whether this gene fragment could code for the entire trnR gene. The mt gene arrangement of T. basalis is different from other Proctotrupomorpha mt genomes, with a number of trn genes in different positions. However, no shared derived gene rearrangements were identified in the present study. Bayesian analyses of mt genomes from 29 hymenopteran taxa and seven other orders of holometabolous insects support some uncontroversial evolutionary relationships, but indicate that much higher levels of taxonomic sampling are necessary for the resolution of family and superfamily relationships.

Complete mitochondrial genome sequence from an endangered Indian snake, Python molurus molurus (Serpentes, Pythonidae)

This paper reports the complete mitochondrial genome sequence of an endangered Indian snake, Python molurus molurus (Indian Rock Python). A typical snake mitochondrial (mt) genome of 17258 bp length comprising of 37 genes including the 13 protein coding genes, 22 tRNA genes, and 2 ribosomal RNA genes along with duplicate control regions is described herein. The P. molurus molurus mt. genome is relatively similar to other snake mt. genomes with respect to gene arrangement, composition, tRNA structures and skews of AT/GC bases. The nucleotide composition of the genome shows that there are more A-C % than T-G% on the positive strand as revealed by positive AT and CG skews. Comparison of individual protein coding genes, with other snake genomes suggests that ATP8 and NADH3 genes have high divergence rates. Codon usage analysis reveals a preference of NNC codons over NNG codons in the mt. genome of P. molurus. Also, the synonymous and non-synonymous substitution rates (ka/ks) suggest that most of the protein coding genes are under purifying selection pressure. The phylogenetic analyses involving the concatenated 13 protein coding genes of P. molurus molurus conformed to the previously established snake phylogeny.

Mitogenomic analysis of the genus Panthera

The complete sequences of the mitochondrial DNA genomes of Panthera tigris, Panthera pardus, and Panthera uncia were determined using the polymerase chain reaction method. The lengths of the complete mitochondrial DNA sequences of the three species were 16990, 16964, and 16773 bp, respectively. Each of the three mitochondrial DNA genomes included 13 protein-coding genes, 22 tRNA, two rRNA, one O(L)R, and one control region. The structures of the genomes were highly similar to those of Felis catus, Acinonyx jubatus, and Neofelis nebulosa. The phylogenies of the genus Panthera were inferred from two combined mitochondrial sequence data sets and the complete mitochondrial genome sequences, by MP (maximum parsimony), ML (maximum likelihood), and Bayesian analysis. The results showed that Panthera was composed of Panthera leo, P. uncia, P. pardus, Panthera onca, P. tigris, and N. nebulosa, which was included as the most basal member. The phylogeny within Panthera genus was N. nebulosa (P. tigris (P. onca (P. pardus, (P. leo, P. uncia)))). The divergence times for Panthera genus were estimated based on the ML branch lengths and four well-established calibration points. The results showed that at about 11.3 MYA, the Panthera genus separated from other felid species and then evolved into the several species of the genus. In detail, N. nebulosa was estimated to be founded about 8.66 MYA, P. tigris about 6.55 MYA, P. uncia about 4.63 MYA, and P. pardus about 4.35 MYA. All these estimated times were older than those estimated from the fossil records. The divergence event, evolutionary process, speciation, and distribution pattern of P. uncia, a species endemic to the central Asia with core habitats on the Qinghai-Tibetan Plateau and surrounding highlands, mostly correlated with the geological tectonic events and intensive climate shifts that happened at 8, 3.6, 2.5, and 1.7 MYA on the plateau during the late Cenozoic period.

[Strong strand specific composition bias-a genomic character of some obligate parasites or symbionts]

DNA replication includes a set of asymmetric mechanisms, which is a division into lagging and leading strands. The former is synthesized continuously whereas the synthesis for the latter is discontinuous. Such a asymmetric mechanism leads to distinct nucleotide composition of these two strands. Strands specific nucleotide composition bias was originally found in genomes of echinoderm and vertebrate mitochondria and then in several bacterial genomes. With the rapid growth in the number of sequenced genomes, many bacteria and even eukaryotes are found to have the consistent strand composition bias. In some bacteria, the extent of strand specific composition bias was so strong that genes on the two replicating strands could be separated according to their codon usages. Till now, 11 obligate intracellular bacteria have been found to have separate codon usages according to whether genes located on the leading or lagging strands. However, there is still not a well-accepted theory that could interpret the reason for the occurrence of separate codon usages in some special bacterial genomes and not in others. This paper reviews the related works and points out its open problems.

Mitochondrial genomes and Doubly Uniparental Inheritance: new insights from Musculista senhousia sex-linked mitochondrial DNAs (Bivalvia Mytilidae)

BACKGROUND

Doubly Uniparental Inheritance (DUI) is a fascinating exception to matrilinear inheritance of mitochondrial DNA (mtDNA). Species with DUI are characterized by two distinct mtDNAs that are inherited either through females (F-mtDNA) or through males (M-mtDNA). DUI sex-linked mitochondrial genomes share several unusual features, such as additional protein coding genes and unusual gene duplications/structures, which have been related to the functionality of DUI. Recently, new evidence for DUI was found in the mytilid bivalve Musculista senhousia. This paper describes the complete sex-linked mitochondrial genomes of this species.

RESULTS

Our analysis highlights that both M and F mtDNAs share roughly the same gene content and order, but with some remarkable differences. The Musculista sex-linked mtDNAs have differently organized putative control regions (CR), which include repeats and palindromic motifs, thought to provide sites for DNA-binding proteins involved in the transcriptional machinery. Moreover, in male mtDNA, two cox2 genes were found, one (M-cox2b) 123bp longer.

CONCLUSIONS

The complete mtDNA genome characterization of DUI bivalves is the first step to unravel the complex genetic signals allowing Doubly Uniparental Inheritance, and the evolutionary implications of such an unusual transmission route in mitochondrial genome evolution in Bivalvia. The observed redundancy of the palindromic motifs in Musculista M-mtDNA may have a role on the process by which sperm mtDNA becomes dominant or exclusive of the male germline of DUI species. Moreover, the duplicated M-COX2b gene may have a different, still unknown, function related to DUI, in accordance to what has been already proposed for other DUI species in which a similar cox2 extension has been hypothesized to be a tag for male mitochondria.

Crocodilian phylogeny inferred from twelve mitochondrial protein-coding genes, with new complete mitochondrial genomic sequences for Crocodylus acutus and Crocodylus novaeguineae

We report complete mitochondrial genomic sequences for Crocodylus acutus and Crocodylus novaeguineae, whose gene orders match those of other crocodilians. Phylogenetic analyses based on the sequences of 12 mitochondrial protein-coding genes support monophyly of two crocodilian taxonomic families, Alligatoridae (genera Alligator, Caiman, and Paleosuchus) and Crocodylidae (genera Crocodylus, Gavialis, Mecistops, Osteolaemus, and Tomistoma). Our results are consistent with monophyly of all crocodilian genera. Within Alligatoridae, genus Alligator is the sister taxon of a clade comprising Caiman and Paleosuchus. Within Crocodylidae, the basal phylogenetic split separates a clade comprising Gavialis and Tomistoma from a clade comprising Crocodylus, Mecistops, and Osteolaemus. Mecistops and Osteolaemus form the sister taxon to Crocodylus. Within Crocodylus, we sampled five Indopacific species, whose phylogenetic ordering is ((C. mindorensis, C. novaeguineae), (C. porosus, (C. siamensis, C. palustris))). The African species C. niloticus and New World species C. acutus form the sister taxon to the Indopacific species, although our sampling lacks three other New World species and an Australian species of Crocodylus.

The complete mitochondrial genome and phylogenetic analysis of the Debao pony (Equus caballus)

The Debao pony (Equus caballus) is the most important local variety of domestic horses, and is strictly protected by the Chinese government. Their average adult withers height is 94.42±3.76 cm for males and 98.35±4.55 cm for females, respectively. In the present study, the complete sequence of the Debao pony mitochondrial genome was determined (GenBank Accession No. EU939445), and was found to be similar to other equine mitochondrial genomes. However, there were 85 nucleotide substitutions in the 13 protein-coding genes; the percentage of substitution was 0.8±0.1. Polymorphisms of mtDNA control regions were analyzed with restriction fragment length polymorphism (RFLP), and 19 haplotypes were found, with a genetic diversity of 0.77. Neighbor-Joining (NJ) and Minimum Evolution (ME) trees based on complete control regions of mtDNA were constructed with the Maximum Composite Likelihood (MCL) method. The analysis indicated that the origins of the Debao pony were scattered in the various branches of the phylogenetic tree. The results from the present study suggest that the Debao pony is derived from multi-matrilineal origins of the species.

A deviant genetic code in the reduced mitochondrial genome of the picoplanktonic green alga Pycnococcus provasolii

Reduction in size of flagellated chlorophytes occurred multiple times during evolution, providing the opportunity to study the consequences of cell reduction on genome architecture. Recent investigations on the chloroplast genomes of the tiny prasinophyceans Ostreococcus tauri (Mamiellales), Micromonas sp. RCC299 (Mamiellales), and Pycnococcus provasolii (Pseudocourfieldiales) highlighted their extreme compaction and reduced gene repertoires. Genome compaction is also exemplified by the Ostreococcus and Micromonas mitochondrial DNAs (mtDNAs) although they have retained almost all of the about 65 genes presumably present in the mitochondria of ancestral prasinophyceans. In this study, the mitochondrial genome of Pycnococcus was sequenced and compared to those of previously examined chlorophytes. Our results document the first case where cellular reduction of a free-living alga was accompanied by marked reduction in gene content of both the mitochondrial and chloroplast genomes. At 24,321 bp, the intronless Pycnococcus mitochondrial genome falls within the lower size range displayed by green algal mtDNAs. The 36 conserved genes, specifying two rRNAs with conventional structures, 16 tRNAs and 18 proteins, are all encoded on the same DNA strand and represent 88% of the genome. Besides a pronounced codon bias, the protein-coding genes feature a variant genetic code characterized by the use of TGA (normally a stop codon) to code for tryptophan, and the unprecedented use of TTA and TTG (normally leucine codons) as stop codons. We conclude that substantial reduction of the mitochondrial genome occurred in at least three independent chlorophyte lineages and that this process entailed a number of convergent changes in these lineages.

Evolution of the mitochondrial genome in snakes: gene rearrangements and phylogenetic relationships

Background

Snakes as a major reptile group display a variety of morphological characteristics pertaining to their diverse behaviours. Despite abundant analyses of morphological characters, molecular studies using mitochondrial and nuclear genes are limited. As a result, the phylogeny of snakes remains controversial. Previous studies on mitochondrial genomes of snakes have demonstrated duplication of the control region and translocation of trnL to be two notable features of the alethinophidian (all serpents except blindsnakes and threadsnakes) mtDNAs. Our purpose is to further investigate the gene organizations, evolution of the snake mitochondrial genome, and phylogenetic relationships among several major snake families.

Results

The mitochondrial genomes were sequenced for four taxa representing four different families, and each had a different gene arrangement. Comparative analyses with other snake mitochondrial genomes allowed us to summarize six types of mitochondrial gene arrangement in snakes. Phylogenetic reconstruction with commonly used methods of phylogenetic inference (BI, ML, MP, NJ) arrived at a similar topology, which was used to reconstruct the evolution of mitochondrial gene arrangements in snakes.

Conclusion

The phylogenetic relationships among the major families of snakes are in accordance with the mitochondrial genomes in terms of gene arrangements. The gene arrangement in Ramphotyphlops braminus mtDNA is inferred to be ancestral for snakes. After the divergence of the early Ramphotyphlops lineage, three types of rearrangements occurred. These changes involve translocations within the IQM tRNA gene cluster and the duplication of the CR. All phylogenetic methods support the placement of Enhydris plumbea outside of the (Colubridae + Elapidae) cluster, providing mitochondrial genomic evidence for the familial rank of Homalopsidae.

Mitogenomic evaluation of the historical biogeography of cichlids toward reliable dating of teleostean divergences

BACKGROUND

Recent advances in DNA sequencing and computation offer the opportunity for reliable estimates of divergence times between organisms based on molecular data. Bayesian estimations of divergence times that do not assume the molecular clock use time constraints at multiple nodes, usually based on the fossil records, as major boundary conditions. However, the fossil records of bony fishes may not adequately provide effective time constraints at multiple nodes. We explored an alternative source of time constraints in teleostean phylogeny by evaluating a biogeographic hypothesis concerning freshwater fishes from the family Cichlidae (Perciformes: Labroidei).

RESULTS

We added new mitogenomic sequence data from six cichlid species and conducted phylogenetic analyses using a large mitogenomic data set. We found a reciprocal monophyly of African and Neotropical cichlids and their sister group relationship to some Malagasy taxa (Ptychochrominae sensu Sparks and Smith). All of these taxa clustered with a Malagasy + Indo/Sri Lankan clade (Etroplinae sensu Sparks and Smith). The results of the phylogenetic analyses and divergence time estimations between continental cichlid clades were much more congruent with Gondwanaland origin and Cretaceous vicariant divergences than with Cenozoic transmarine dispersal between major continents.

CONCLUSION

We propose to add the biogeographic assumption of cichlid divergences by continental fragmentation as effective time constraints in dating teleostean divergence times. We conducted divergence time estimations among teleosts by incorporating these additional time constraints and achieved a considerable reduction in credibility intervals in the estimated divergence times.

The complete mitochondrial genome of the Antarctic springtail Cryptopygus antarcticus (Hexapoda: Collembola)

BACKGROUND

Mitogenomics data, i.e. complete mitochondrial genome sequences, are popular molecular markers used for phylogenetic, phylogeographic and ecological studies in different animal lineages. Their comparative analysis has been used to shed light on the evolutionary history of given taxa and on the molecular processes that regulate the evolution of the mitochondrial genome. A considerable literature is available in the fields of invertebrate biochemical and ecophysiological adaptation to extreme environmental conditions, exemplified by those of the Antarctic. Nevertheless, limited molecular data are available from terrestrial Antarctic species, and this study represents the first attempt towards the description of a mitochondrial genome from one of the most widespread and common collembolan species of Antarctica.

RESULTS

In this study we describe the mitochondrial genome of the Antarctic collembolan Cryptopygus antarcticus Willem, 1901. The genome contains the standard set of 37 genes usually present in animal mtDNAs and a large non-coding fragment putatively corresponding to the region (A+T-rich) responsible for the control of replication and transcription. All genes are arranged in the gene order typical of Pancrustacea. Three additional short non-coding regions are present at gene junctions. Two of these are located in positions of abrupt shift of the coding polarity of genes oriented on opposite strands suggesting a role in the attenuation of the polycistronic mRNA transcription(s). In addition, remnants of an additional copy of trnL(uag) are present between trnS(uga) and nad1. Nucleotide composition is biased towards a high A% and T% (A+T = 70.9%), as typically found in hexapod mtDNAs. There is also a significant strand asymmetry, with the J-strand being more abundant in A and C. Within the A+T-rich region, some short sequence fragments appear to be similar (in position and primary sequence) to those involved in the origin of the N-strand replication of the Drosophila mtDNA.

CONCLUSION

The mitochondrial genome of C. antarcticus shares several features with other pancrustacean genomes, although the presence of unusual non-coding regions is also suggestive of molecular rearrangements that probably occurred before the differentiation of major collembolan families. Closer examination of gene boundaries also confirms previous observations on the presence of unusual start and stop codons, and suggests a role for tRNA secondary structures as potential cleavage signals involved in the maturation of the primary transcript. Sequences potentially involved in the regulation of replication/transcription are present both in the A+T-rich region and in other areas of the genome. Their position is similar to that observed in a limited number of insect species, suggesting unique replication/transcription mechanisms for basal and derived hexapod lineages. This initial description and characterization of the mitochondrial genome of C. antarcticus will constitute the essential foundation prerequisite for investigations of the evolutionary history of one of the most speciose collembolan genera present in Antarctica and other localities of the Southern Hemisphere.