A resourceful genome: updating the functional repertoire and evolutionary role of animal mitochondrial DNAs

Morpometric and molecular characterization of Surguli goat through CO1 gene in district Kohat

The present study was designed with the aim to study morphometric characterization as well as phylogeny and diversity of the local Surguli goat at their breeding tract district Kohat through mitochondrial DNA region, i.e., Cytochrome C Oxidase Subunit One (CO1) gene. Morphometric data and blood samples were collected from thirty (30) pure goats. Morphometric analysis showed that sex had significant effect (p < 0.05) on body weight, body length, hearth girth and horn length while no significant effect (p > 0.05) was observed for other characteristics. The results also indicated that age had significant effect (p < 0.05) on height at rump, ear length, horn length and tail length while no significant effect (p > 0.05) was observed for other characteristics. The phylogenetic analysis through CO1 nucleotide sequences within nucleotide range 1-767 showed nine polymorphic sites segregating into eight haplotypes. The mean intraspecific diversity and mean interspecific diversity were calculated as 0.23 and 2.36%, respectively. Phylogenetic tree revealed that Capra Ibex and native Surguli goat have common ancestors. The morphometric and molecular results obtained from the present study can be exploited as a selection tool for breeding and overall improvement.

Mitochondrially mediated RNA interference, a retrograde signaling system affecting nuclear gene expression

Several functional classes of short noncoding RNAs are involved in manifold regulatory processes in eukaryotes, including, among the best characterized, miRNAs. One of the most intriguing regulatory networks in the eukaryotic cell is the mito-nuclear crosstalk: recently, miRNA-like elements of mitochondrial origin, called smithRNAs, were detected in a bivalve species, Ruditapes philippinarum. These RNA molecules originate in the organelle but were shown in vivo to regulate nuclear genes. Since miRNA genes evolve easily de novo with respect to protein-coding genes, in the present work we estimate the probability with which a newly arisen smithRNA finds a suitable target in the nuclear transcriptome. Simulations with transcriptomes of 12 bivalve species suggest that this probability is high and not species specific: one in a hundred million (1 × 10-8) if five mismatches between the smithRNA and the 3' mRNA are allowed, yet many more are allowed in animals. We propose that novel smithRNAs may easily evolve as exaptation of the pre-existing mitochondrial RNAs. In turn, the ability of evolving novel smithRNAs may have played a pivotal role in mito-nuclear interactions during animal evolution, including the intriguing possibility of acting as speciation trigger.

Comparative mitogenomics of Brachiopods reveals conservatism in articulate species and unusualness in inarticulate species

BACKGROUND

Brachiopods are a phylum of marine invertebrates with over 10,000 fossil species. Today, there are fewer than 500 extant species assigned to the class Articulata or Inarticulata and for which knowledge of evolutionary genetics and genomics is still poor. Until now, complete mitogenome sequences of two inarticulate species and four articulate species were available.

METHODS AND RESULTS

The complete mitogenome of the inarticulate brachiopod species Lingula reevii (20,778 bp) was obtained by using next generation sequencing. It contains 12 protein-coding genes (the annotation of atp8 is unsure), two ribosomal RNA genes, 26 transfer RNA genes, and one supernumerary ORF that is also conserved in the inarticulate species Lingula anatina. It is hypothesized that this ORF could represent a Lingula-specific mtORFan gene (without obvious homology to other genes). Comparative mitogenomics indicate the mitochondrial gene order of L. reevii is unique among brachiopods, and that compared to articulate species, inarticulate species exhibit massive mitogenome rearrangements, deviant ATP8 protein sequences and supernumerary ORFs, possibly representing species- or lineage-specific mtORFan genes.

CONCLUSION

The results of this study enrich genetics knowledge of extant brachiopods, which may eventually help to test hypotheses about their decline.

High-altitude adaptation is accompanied by strong signatures of purifying selection in the mitochondrial genomes of three Andean waterfowl

Evidence from a variety of organisms points to convergent evolution on the mitochondria associated with a physiological response to oxygen deprivation or temperature stress, including mechanisms for high-altitude adaptation. Here, we examine whether demography and/or selection explains standing mitogenome nucleotide diversity in high-altitude adapted populations of three Andean waterfowl species: yellow-billed pintail (Anas georgica), speckled teal (Anas flavirostris), and cinnamon teal (Spatula cyanoptera). We compared a total of 60 mitogenomes from each of these three duck species (n = 20 per species) across low and high altitudes and tested whether part(s) or all of the mitogenome exhibited expected signatures of purifying selection within the high-altitude populations of these species. Historical effective population sizes (Ne) were inferred to be similar between high- and low-altitude populations of each species, suggesting that selection rather than genetic drift best explains the reduced genetic variation found in mitochondrial genes of high-altitude populations compared to low-altitude populations of the same species. Specifically, we provide evidence that establishment of these three Andean waterfowl species in the high-altitude environment, coincided at least in part with a persistent pattern of negative purifying selection acting on oxidative phosphorylation (OXPHOS) function of the mitochondria. Our results further reveal that the extent of gene-specific purifying selection has been greatest in the speckled teal, the species with the longest history of high-altitude occupancy.

Comparative Mitogenomics Reveals Cryptic Species in Sillago ingenuua McKay, 1985 (Perciformes: Sillaginidae)

It is unreliable to identify marine fishes only by external morphological features. Species misidentification brings great challenges to fishery research, resource monitoring and ecomanagement. Sillago ingenuua is an important part of commercial marine fishes, and in which, the morphological differences between different groups are not obvious. Here, we compared different geographical groups of S. ingenuua which were collected from Xiamen, Dongshan, Keelung, Songkhla and Java. The results showed that all samples of S. ingenuua were similar in external morphological characteristics and the shape of the swim bladder, but there were two distinctive lineages which were flagged as cryptic species based on DNA barcoding. The comparative mitogenomic results showed that S. ingenuua A and S. ingenuua B were identical in structural organization and gene arrangement. Their nucleotide composition and codon usage were also similar. A phylogenetic analysis was performed based on 13 concatenated PCGs from eight Sillago species. The results showed that the genetic distance between S. ingenuua A and S. ingenuua B was large (D = 0.069), and this genetic distance was large enough to reveal that S. ingenuua A and S. ingenuua B might be different species.

Comparative analysis of the mitochondrial genomes of the family Mactridae (Mollusca: Venerida) and their phylogenetic implications

Taxonomy and phylogenetic relationships within the family Mactridae have remained debatable because of the plasticity of morphological characteristics and the lack of accurate molecular data, thereby resulting in abundant synonyms and taxa rearrangements. Mitochondrial genomes (mitogenomes) have been widely used as powerful tools to reconstruct phylogenies of various groups of mollusks; however, they have not been used for studying the phylogeny of mactrids specifically. In the present study, mitogenomes of seven Mactridae species, namely Mactra chinensis, Mactra cygnus, Mactra quadrangularis, Mactra cumingii, Mactrinula dolabrata, Raeta pulchella, and Raeta sp., were sequenced by Illumina high-throughput sequencing, and a comparative mitochondrial genomic analysis was conducted. The newly sequenced mitogenomes were double-stranded circular molecules, with all functional genes encoded on the heavy strand. All the new mactrid mitogenomes had two rRNA genes (12S and 16S), 13 protein-coding genes (PCGs) (atp6, cox1, cox2, cox3, cytb, nad1, nad2, nad3, nad4, nad4l, nad5, nad6, and atp8), and 22 tRNAs. The mitogenomes showed considerable variation in AT content, GC skew, and AT skew. The results of the phylogenetic analysis confirmed monophyly of the family Mactridae and suggested that genera Mactrinula, Spisula, Rangia, and Mulinia should not be placed under subfamily Mactrinae. Our results supported that potential cryptic species existed in Mactra antiquata. We also proposed subfamily Kymatoxinae should belong to the family Mactridae rather than Anatinellidae and Mactra alta in China should be Mactra cygnus. Additionally, conservation in functional gene arrangement was found in genera Mactra, Raeta, and Lutraria. But gene orders in S. sachalinensis and S. solida were quite different, questioning their congeneric relationship. Our results further suggested that the taxonomy within the family Mactridae requires an integrative revision.

LAMP-based molecular sexing in a gonochoric marine bivalve (Macoma balthica rubra) with divergent sex-specific mitochondrial genomes

Taking advantage of the unique system of doubly uniparental inheritance (DUI) of mitochondria, we developed a reliable molecular method to sex individuals of the marine bivalve Macoma balthica rubra. In species with DUI (~100 known bivalves), both sexes transmit their mitochondria: males bear both a male- and female-type mitogenome, while females bear only the female type. Male and female mitotypes are sufficiently divergent to reliably PCR-amplify them specifically. Loop-mediated isothermal amplification (LAMP) is a precise, economical and portable alternative to PCR for molecular sexing and we demonstrate its application in this context. We used 154 individuals sampled along the Atlantic coast of France and sexed microscopically by gonad examination to test for the congruence among gamete type, PCR sexing and LAMP sexing. We show an exact match among the sexing results from these three methods using the male and female mt-cox1 genes. DUI can be disrupted in inter-specific hybrids, causing unexpected distribution of mitogenomes, such as homoplasmic males or heteroplasmic females. To our knowledge, DUI disruption at the intra-specific scale has never been tested. We applied our sexing protocol to control for unexpected heteroplasmy caused by hybridization between divergent genetic lineages and found no evidence of disruption in the mode of mitochondrial inheritance in M. balthica rubra. We propose LAMP as a useful tool to accelerate eco-evolutionary studies of DUI. It offers the opportunity to investigate the potential role of, previously unaccounted-for, sex-specific patterns such as sexual selection or sex-specific dispersal bias in the evolution of free-spawning benthic species.

Phylogenomic resolution of Imparidentia (Mollusca: Bivalvia) diversification through mitochondrial genomes

Despite significant advances in the phylogenomics of bivalves over the past decade, the higher-level phylogeny of Imparidentia (a superorder of Heterodonta) remains elusive. Here, a total of five new mitochondrial sequences (Chama asperella, Chama limbula, Chama dunkeri, Barnea manilensis and Ctena divergens) was added to provide resolution in nodes that required additional study. Although the monophyly of Lucinida remains less clear, the results revealed the overall backbone of the Imparidentia tree and the monophyly of Imparidentia. Likewise, most relationships among the five major Imparidentia lineages-Lucinida, Cardiida, Adapedonta, Myida and Venerida-were addressed with a well-supported topology. Basal relationships of Imparidentia recovered Lucinidae as the sister group to all remaining imparidentian taxa. Thyasiridae is a sister group to other imparidentian bivalves (except Lucinidae species) which is split into Cardiida, Adapedonta and the divergent clade of Neoheterodontei. Neoheterodontei was comprised of Venerida and Myida, the former of which now also contains Chamidae as the sister group to all the remaining venerid taxa. Moreover, molecular divergence times were inferred by calibrating nine nodes in the Imparidentia tree of life by extinct taxa. The origin of these major clades ranged from Ordovician to Permian with the diversification through the Palaeozoic to Mesozoic. Overall, the results obtained in this study demonstrate a better-resolved Imparidentia phylogeny based on mitochondrial genomes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-023-00178-x.

A tale of two paths: The evolution of mitochondrial recombination in bivalves with doubly uniparental inheritance

In most animals, mitochondrial DNA is strictly maternally inherited and non-recombining. One exception to this pattern is called doubly uniparental inheritance (DUI), a phenomenon involving the independent transmission of female and male mitochondrial genomes. DUI is known only from the molluskan class Bivalvia. The phylogenetic distribution of male-transmitted mitochondrial DNA (M mtDNA) in bivalves is consistent with several evolutionary scenarios, including multiple independent gains, losses, and varying degrees of recombination with female-transmitted mitochondrial DNA (F mtDNA). In this study, we use phylogenetic methods to test M mtDNA origination hypotheses and infer the prevalence of mitochondrial recombination in bivalves with DUI. Phylogenetic modeling using site concordance factors supported a single origin of M mtDNA in bivalves coupled with recombination acting over long evolutionary timescales. Ongoing mitochondrial recombination is present in Mytilida and Venerida, which results in a pattern of concerted evolution of F mtDNA and M mtDNA. Mitochondrial recombination could be favored to offset the deleterious effects of asexual inheritance and maintain mitonuclear compatibility across tissues. Cardiida and Unionida have gone without recent recombination, possibly due to an extension of the COX2 gene in male mitochondrial DNA. The loss of recombination could be connected to the role of M mtDNA in sex determination or sexual development. Our results support that recombination events may occur throughout the mitochondrial genomes of DUI species. Future investigations may reveal more complex patterns of inheritance of recombinants, which could explain the retention of signal for a single origination of M mtDNA in protein-coding genes.

Ancient endosymbiont-mediated transmission of a selfish gene provides a model for overcoming barriers to gene transfer into animal mitochondrial genomes

In contrast to bilaterian animals, non-bilaterian mitochondrial genomes contain atypical genes, often attributed to horizontal gene transfer (HGT) as an ad hoc explanation. Although prevalent in plants, HGT into animal mitochondrial genomes is rare, lacking suitable explanatory models for their occurrence. HGT of the mismatch DNA repair gene (mtMutS) from giant viruses to octocoral (soft corals and their kin) mitochondrial genomes provides a model for how barriers to HGT to animal mitochondria may be overcome. A review of the available literature suggests that this HGT was mediated by an alveolate endosymbiont infected with a lysogenic phycodnavirus that enabled insertion of the homing endonuclease containing mtMutS into octocoral mitochondrial genomes. We posit that homing endonuclease domains and similar selfish elements play a crucial role in such inter-domain gene transfers. Understanding the role of selfish genetic elements in HGT has the potential to aid development of tools for manipulating animal mitochondrial DNA.

A comparative analysis of mitochondrial ORFs provides new insights on expansion of mitochondrial genome size in Arcidae

BACKGROUND

Arcidae, comprising about 260 species of ark shells, is an ecologically and economically important lineage of bivalve mollusks. Interestingly, mitochondrial genomes of several Arcidae species are 2-3 times larger than those of most bilaterians, and are among the largest bilaterian mitochondrial genomes reported to date. The large mitochondrial genome size is mainly due to expansion of unassigned regions (regions that are functionally unassigned). Previous work on unassigned regions of Arcidae mtDNA genomes has focused on nucleotide-level analyses to observe sequence characteristics, however the origin of expansion remains unclear.

RESULTS

We assembled six new mitogenomes and sequenced six transcriptomes of Scapharca broughtonii to identify conserved functional ORFs that are transcribed in unassigned regions. Sixteen lineage-specific ORFs with different copy numbers were identified from seven Arcidae species, and 11 of 16 ORFs were expressed and likely biologically active. Unassigned regions of 32 Arcidae mitogenomes were compared to verify the presence of these novel mitochondrial ORFs and their distribution. Strikingly, multiple structural analyses and functional prediction suggested that these additional mtDNA-encoded proteins have potential functional significance. In addition, our results also revealed that the ORFs have a strong connection to the expansion of Arcidae mitochondrial genomes and their large-scale duplication play an important role in multiple expansion events. We discussed the possible origin of ORFs and hypothesized that these ORFs may originate from duplication of mitochondrial genes.

CONCLUSIONS

The presence of lineage-specific mitochondrial ORFs with transcriptional activity and potential functional significance supports novel features for Arcidae mitochondrial genomes. Given our observation and analyses, these ORFs may be products of mitochondrial gene duplication. These findings shed light on the origin and function of novel mitochondrial genes in bivalves and provide new insights into evolution of mitochondrial genome size in metazoans.

Mitochondrial genomic analyses provide new insights into the "missing" atp8 and adaptive evolution of Mytilidae

Background

Mytilidae, also known as marine mussels, are widely distributed in the oceans worldwide. Members of Mytilidae show a tremendous range of ecological adaptions, from the species distributed in freshwater to those that inhabit in deep-sea. Mitochondria play an important role in energy metabolism, which might contribute to the adaptation of Mytilidae to different environments. In addition, some bivalve species are thought to lack the mitochondrial protein-coding gene ATP synthase F0 subunit 8. Increasing studies indicated that the absence of atp8 may be caused by annotation difficulties for atp8 gene is characterized by highly divergent, variable length.

Results

In this study, the complete mitochondrial genomes of three marine mussels (Xenostrobus securis, Bathymodiolus puteoserpentis, Gigantidas vrijenhoeki) were newly assembled, with the lengths of 14,972 bp, 20,482, and 17,786 bp, respectively. We annotated atp8 in the sequences that we assembled and the sequences lacking atp8. The newly annotated atp8 sequences all have one predicted transmembrane domain, a similar hydropathy profile, as well as the C-terminal region with positively charged amino acids. Furthermore, we reconstructed the phylogenetic trees and performed positive selection analysis. The results showed that the deep-sea bathymodiolines experienced more relaxed evolutionary constraints. And signatures of positive selection were detected in nad4 of Limnoperna fortunei, which may contribute to the survival and/or thriving of this species in freshwater.

Conclusions

Our analysis supported that atp8 may not be missing in the Mytilidae. And our results provided evidence that the mitochondrial genes may contribute to the adaptation of Mytilidae to different environments.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08940-8.

The longest mitochondrial protein in metazoans is encoded by the male-transmitted mitogenome of the bivalve Scrobicularia plana

Cytochrome c oxidase subunit II (COX2) is one of the three mitochondrially encoded proteins of the complex IV of the respiratory chain that catalyses the reduction of oxygen to water. The cox2 gene spans about 690 base pairs in most animal species and produces a protein composed of approximately 230 amino acids. We discovered an extreme departure from this pattern in the male-transmitted mitogenome of the bivalve Scrobicularia plana with doubly uniparental inheritance (DUI) of mitochondrial DNA (mtDNA), which possesses an important in-frame insertion of approximately 4.8 kb in its cox2 gene. This feature—an enlarged male cox2 gene—is found in many species with DUI; the COX2 protein can be up to 420 amino acids long. Through RT-PCRs, immunoassays and comparative genetics, the evolution and functionality of this insertion in S. plana were characterized. The in-frame insertion is conserved among individuals from different populations and bears the signature of purifying selection seemingly indicating maintenance of functionality. Its transcription and translation were confirmed: this gene produces a polypeptide of 1892 amino acids, making it the largest metazoan COX2 protein known to date. We hypothesize that these extreme modifications in the COX2 protein affect the metabolism of mitochondria containing the male-transmitted mtDNA in Scrobicularia plana.

Did doubly uniparental inheritance (DUI) of mtDNA originate as a cytoplasmic male sterility (CMS) system?

Animal and plant species exhibit an astonishing diversity of sexual systems, including environmental and genetic determinants of sex, with the latter including genetic material in the mitochondrial genome. In several hermaphroditic plants for example, sex is determined by an interaction between mitochondrial cytoplasmic male sterility (CMS) genes and nuclear restorer genes. Specifically, CMS involves aberrant mitochondrial genes that prevent pollen development and specific nuclear genes that restore it, leading to a mixture of female (male-sterile) and hermaphroditic individuals in the population (gynodioecy). Such a mitochondrial-nuclear sex determination system is thought to be rare outside plants. Here, we present one possible case of CMS in animals. We hypothesize that the only exception to the strict maternal mtDNA inheritance in animals, the doubly uniparental inheritance (DUI) system in bivalves, might have originated as a mitochondrial-nuclear sex-determination system. We document and explore similarities that exist between DUI and CMS, and we propose various ways to test our hypothesis.

Characterisation of the Complete Mitochondrial Genome of Critically Endangered Mustela lutreola (Carnivora: Mustelidae) and Its Phylogenetic and Conservation Implications

In this paper, a complete mitochondrial genome of the critically endangered European mink Mustela lutreola L., 1761 is reported. The mitogenome was 16,504 bp in length and encoded the typical 13 protein-coding genes, two ribosomal RNA genes and 22 transfer RNA genes, and harboured a putative control region. The A+T content of the entire genome was 60.06% (A > T > C > G), and the AT-skew and GC-skew were 0.093 and −0.308, respectively. The encoding-strand identity of genes and their order were consistent with a collinear gene order characteristic for vertebrate mitogenomes. The start codons of all protein-coding genes were the typical ATN. In eight cases, they were ended by complete stop codons, while five had incomplete termination codons (TA or T). All tRNAs had a typical cloverleaf secondary structure, except tRNA^Ser(AGC) and tRNA^Lys, which lacked the DHU stem and had reduced DHU loop, respectively. Both rRNAs were capable of folding into complex secondary structures, containing unmatched base pairs. Eighty-one single nucleotide variants (substitutions and indels) were identified. Comparative interspecies analyses confirmed the close phylogenetic relationship of the European mink to the so-called ferret group, clustering the European polecat, the steppe polecat and the black-footed ferret. The obtained results are expected to provide useful molecular data, informing and supporting effective conservation measures to save M. lutreola.

The Mitochondrial Genome of a Freshwater Pelagic Amphipod Macrohectopus branickii Is among the Longest in Metazoa

There are more than 350 species of amphipods (Crustacea) in Lake Baikal, which have emerged predominantly through the course of endemic radiation. This group represents a remarkable model for studying various aspects of evolution, one of which is the evolution of mitochondrial (mt) genome architectures. We sequenced and assembled the mt genome of a pelagic Baikalian amphipod species Macrohectopus branickii. The mt genome is revealed to have an extraordinary length (42,256 bp), deviating significantly from the genomes of other amphipod species and the majority of animals. The mt genome of M. branickii has a unique gene order within amphipods, duplications of the four tRNA genes and Cox2, and a long non-coding region, that makes up about two thirds of the genome’s size. The extension of the mt genome was most likely caused by multiple duplications and inversions of regions harboring ribosomal RNA genes. In this study, we analyzed the patterns of mt genome length changes in amphipods and other animal phyla. Through a statistical analysis, we demonstrated that the variability in the mt genome length may be a characteristic of certain phyla and is primarily conferred by expansions of non-coding regions.

The complete mitogenome of the inarticulate brachiopod Glottidia pyramidata reveals insights into gene order variation, deviant ATP8 and mtORFans in the Brachiopoda

Brachiopods are a clade of marine organisms with a tremendously diverse and abundant fossil record but with fewer than 500 species extant today. Even if a better understanding of their biology and genetics could help to test hypotheses about their impressive decline, knowledge of genetics and evolutionary genomics in extant brachiopods is very poor. Here, we present the complete mitochondrial genome sequence of the inarticulate Glottidia pyramidata, an eastern North American extant representative of the phylum Brachiopoda. Besides the general characteristics of the sequenced mitogenome, we present its unusual features such as deviant ATP8 protein sequence and supernumerary ORFs, and also unique gene order, considering the available genome sequences of other brachiopod species.

Mitochondrial Genomic Landscape: A Portrait of the Mitochondrial Genome 40 Years after the First Complete Sequence

Notwithstanding the initial claims of general conservation, mitochondrial genomes are a largely heterogeneous set of organellar chromosomes which displays a bewildering diversity in terms of structure, architecture, gene content, and functionality. The mitochondrial genome is typically described as a single chromosome, yet many examples of multipartite genomes have been found (for example, among sponges and diplonemeans); the mitochondrial genome is typically depicted as circular, yet many linear genomes are known (for example, among jellyfish, alveolates, and apicomplexans); the chromosome is normally said to be “small”, yet there is a huge variation between the smallest and the largest known genomes (found, for example, in ctenophores and vascular plants, respectively); even the gene content is highly unconserved, ranging from the 13 oxidative phosphorylation-related enzymatic subunits encoded by animal mitochondria to the wider set of mitochondrial genes found in jakobids. In the present paper, we compile and describe a large database of 27,873 mitochondrial genomes currently available in GenBank, encompassing the whole eukaryotic domain. We discuss the major features of mitochondrial molecular diversity, with special reference to nucleotide composition and compositional biases; moreover, the database is made publicly available for future analyses on the MoZoo Lab GitHub page.

A proposed method for analyzing molecular signatures to detect hermaphroditism in freshwater mussels: a case study using the eastern floater (Pyganodon cataracta)

Many freshwater mussels (order Unionida) have an unusual system of doubly uniparental inheritance (DUI) of mitochondrial (mt) DNA. In species with DUI, males possess a female-transmitted (F-type) mt genome and a male-transmitted (M-type) mt genome. These genomes contain non-canonical open reading frame (orf) genes, referred to as f-orf and m-orf, present in F and M mt genomes, respectively. These genes have been implicated in sexual development in Unionida. When gonochoric species become hermaphroditic, which has happened several times in Unionida, they lose their M-type mt genome and f-orf genes evolve dramatically. Resulting F-ORF proteins are highly divergent in terms of primary nucleotide sequence, inferred amino acids, and hydrophobic properties; these genes (and proteins) are referred to as hermaphroditic orfs or h-orfs (and H-ORFs). We investigated patterns of hydrophobicity divergence for H-ORF proteins in hermaphrodites versus F-ORF proteins in closely related gonochoric species against cytochrome c oxidase subunit 1 (cox1) divergences. This approach was used to assess whether cryptic hermaphrodites can be detected. Although we did not detect evidence for the recent transition of any populations of eastern floaters (Pyganodon cataracta (Say, 1817)) to hermaphroditism, our analyses demonstrate that molecular signatures in mtDNA can be used to detect hermaphroditism in freshwater mussels.

The inevitability of ATP as a transmitter in the carotid body

Atmospheric oxygen concentrations rose markedly at several points in evolutionary history. Each of these increases was followed by an evolutionary leap in organismal complexity, and thus the cellular adaptions we see today have been shaped by the levels of oxygen within our atmosphere. In eukaryotic cells, oxygen is essential for the production of adenosine 5'-triphosphate (ATP) which is the 'Universal Energy Currency' of life. Aerobic organisms survived by evolving precise mechanisms for converting oxygen within the environment into energy. Higher mammals developed specialised organs for detecting and responding to changes in oxygen content to maintain gaseous homeostasis for survival. Hypoxia is sensed by the carotid bodies, the primary chemoreceptor organs which utilise multiple neurotransmitters one of which is ATP to evoke compensatory reflexes. Yet, a paradox is presented in oxygen sensing cells of the carotid body when during periods of low oxygen, ATP is seemingly released in abundance to transmit this signal although the synthesis of ATP is theoretically halted because of its dependence on oxygen. We propose potential mechanisms to maintain ATP production in hypoxia and summarise recent data revealing elevated sensitivity of purinergic signalling within the carotid body during conditions of sympathetic overactivity and hypertension. We propose the carotid body is hypoxic in numerous chronic cardiovascular and respiratory diseases and highlight the therapeutic potential for modulating purinergic transmission.

Bivalve molluscs as model systems for studying mitochondrial biology

The class Bivalvia is a highly successful and ancient taxon including ∼25,000 living species. During their long evolutionary history bivalves adapted to a wide range of physicochemical conditions, habitats, biological interactions, and feeding habits. Bivalves can have strikingly different size, and despite their apparently simple body plan, they evolved very different shell shapes, and complex anatomic structures. One of the most striking features of this class of animals is their peculiar mitochondrial biology: some bivalves have facultatively anaerobic mitochondria that allow them to survive prolonged periods of anoxia/hypoxia. Moreover, more than 100 species have now been reported showing the only known evolutionarily stable exception to the strictly maternal inheritance of mitochondria in animals, named doubly uniparental inheritance. Mitochondrial activity is fundamental to eukaryotic life, and thanks to their diversity and uncommon features, bivalves represent a great model system to expand our knowledge about mitochondrial biology, so far limited to a few species. We highlight recent works studying mitochondrial biology in bivalves at either genomic or physiological level. A link between these two approaches is still missing, and we believe that an integrated approach and collaborative relationships are the only possible ways to be successful in such endeavour.

Mitochondrial Short-Term Plastic Responses and Long-Term Evolutionary Dynamics in Animal Species

How do species respond or adapt to environmental changes? The answer to this depends partly on mitochondrial epigenetics and genetics, new players in promoting adaptation to both short- and long-term environmental changes. In this review, we explore how mitochondrial epigenetics and genetics mechanisms, such as mtDNA methylation, mtDNA-derived noncoding RNAs, micropeptides, mtDNA mutations, and adaptations, can contribute to animal plasticity and adaptation. We also briefly discuss the challenges in assessing mtDNA adaptive evolution. In sum, this review covers new advances in the field of mitochondrial genomics, many of which are still controversial, and discusses processes still somewhat obscure, and some of which are still quite speculative and require further robust experimentation.

HERMES: An improved method to test mitochondrial genome molecular synapomorphies among clades

Mitochondrial chromosomes have diversified among eukaryotes and many different architectures and features are now acknowledged for this genome. Here we present the improved HERMES index, which can measure and quantify the amount of molecular change experienced by mitochondrial genomes. We test the improved approach with ten molecular phylogenetic studies based on complete mitochondrial genomes, representing six bilaterian Phyla. In most cases, HERMES analysis spotted out clades or single species with peculiar molecular synapomorphies, allowing to identify phylogenetic and ecological patterns. The software presented herein handles linear, circular, and multi-chromosome genomes, thus widening the HERMES scope to the complete eukaryotic domain.

A Naturally Heteroplasmic Clam Provides Clues about the Effects of Genetic Bottleneck on Paternal mtDNA

Mitochondrial DNA (mtDNA) is present in multiple copies within an organism. Since these copies are not identical, a single individual carries a heterogeneous population of mtDNAs, a condition known as heteroplasmy. Several factors play a role in the dynamics of the within-organism mtDNA population: among them, genetic bottlenecks, selection, and strictly maternal inheritance are known to shape the levels of heteroplasmy across mtDNAs.

In Metazoa, the only evolutionarily stable exception to the strictly maternal inheritance of mitochondria is the doubly uniparental inheritance (DUI), reported in 100+ bivalve species. In DUI species, there are two highly divergent mtDNA lineages, one inherited through oocyte mitochondria (F-type) and the other through sperm mitochondria (M-type). Having both parents contributing to the mtDNA pool of the progeny makes DUI a unique system to study the dynamics of mtDNA populations. Since, in bivalves, the spermatozoon has few mitochondria (4–5), M-type mtDNA faces a tight bottleneck during embryo segregation, one of the narrowest mitochondrial bottlenecks investigated so far.

Here, we analyzed the F- and M-type mtDNA variability within individuals of the DUI species Ruditapes philippinarum and investigated for the first time the effects of such a narrow bottleneck affecting mtDNA populations. As a potential consequence of this narrow bottleneck, the M-type mtDNA shows a large variability in different tissues, a condition so pronounced that it leads to genotypes from different tissues of the same individual not to cluster together. We believe that such results may help understanding the effect of low population size on mtDNA bottleneck.

Characterization of the complete mitochondrial genome of Macrotocinclus affinis (Siluriformes; Loricariidae) and phylogenetic studies of Siluriformes

The vertebrate mitochondrial genome is typically circular molecules made up of 14,000 to 16,000 bp, including 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (12 s rRNA and 16 s rRNA) and a control region. Compared with nuclear DNA, mitochondrial DNA has a higher mutation rate, so it is one of the most effective and reliable molecular markers in fish phylogeny. Macrotocinclus affinis was the only species in Macrotocinclus (it was classified as Otocinclus in the past) and currently lacks genetic information. Most of the current researches are based on the mitochondrial Cytb gene and RAG1 and RAG2 nuclear genes to study the phylogenetic analysis of Siluriformes. So, the study provides the characteristic features of the Macrotocinclus affinis mitochondrial genome and this is the first time that the phylogenetic relationship of Siluriformes has been reconstructed based on COI. We aimed to sequence the entire mitochondrial genome of Macrotocinclus affinis using conventional PCR techniques and to clarify its phylogenetic status in Siluriformes by using the COI sequence of mitochondria. In this study, we sequenced the whole mitochondrial genome of this species yielding a 16,632 bp circular assembly composed of the typical vertebrate mitochondrial features. It contains 13 protein-coding genes, two rRNA genes, 22 tRNA genes, a putative control region, and one origin of replication on the light-strand. The overall base composition includes A (30.07%), T (24.43%), C (29.43%) and G (16.01%). The genome composition is A + T biased (54.5%), and exhibits AT-skew (0.1036) and GC-skew (-0.2962). Moreover, the 13 PCGs encode 3850 amino acids in total. The result of the phylogenetic tree supports Macrotocinclus affinis has a closest relationship with Otocinclus cf. hoppei far. These results will help to understand the characteristics of the mitochondrial genome of Macrotocinclus affinis and provide molecular basis for the evolutionary relationship of Loricariidae.

Genetic characterization and phylogenetic study of Lakor goat from Southwest Maluku Regency based on mitochondrial COI gene

Aim:

This study is aimed at characterizing the genetic and phylogenetic structure of Lakor goats as indigenous livestock from the Southwest Maluku Regency based on mitochondrial COI gene sequences.

Materials and Methods:

The genomes of 103 follicle samples from Lakor goats, collected from Lakor Island, were analyzed. The polymerase chain reaction was used to amplify 1548 bp of the mitochondrial COI gene using two primer pairs (COIA and COIB). Following sequencing, genetic variation and phylogenetic relationship were established using MEGA version X software.

Results:

The results of multiple COI gene alignment of the total sequences identified four polymorphic nucleotides that function as genetic markers between individual animals within the Lakor goat population. These correspond to positions 228 (A-G), 519 (G-A), 900 (C-T), and 1266 (T-C). Phylogenetic signals based on the COI gene showed that Lakor goat breed is a monophyletic group or single clade with a bootstrap value of 100% by the neighbor-joining (NJ) and maximum likelihood (ML) evolutionary models. This data indicated that evolutionarily, the Lakor goat breed has a very close kinship with three goat breeds from China: The Meigu goat (KM 244714.1), Chinese Tibet (Capra hircus) (KJ 940969.1), and C. hircus (KP 677510.1). Phylogenetic information based on the cladistics system classified the Lakor goat as a single clade (monophyletic group). The low-genetic diversity within populations indicates that there has been an inbreeding depression occurring at a very high frequency.

Conclusion:

We conclude that the Lakor goat may be divided into a single clade or monophyletic group based on the COI gene sequence. Four nucleotides were identified that can be used as genetic markers among individual animals within the Lakor goat population, as well as C. hircus and others as derived from GenBank data. The Lakor goat population has a high level of inbreeding depression as a result of geographical isolation, which supports the formation of a monophyletic group with different genetic characteristics, and does not allow the introduction of males from other breeds. Phylogenetic signals indicated that Capra aegagrus (bezoar) is the ancestor of the native goats in Indonesia, including the Lakor goats.

Variation in protein gene and intron content among land plant mitogenomes

The functional content of the mitochondrial genome (mitogenome) is highly diverse across eukaryotes. Among land plants, our understanding of the variation in mitochondrial gene and intron content is improving from concerted efforts to densely sample mitogenomes from diverse land plants. Here I review the current state of knowledge regarding the diversity in content of protein genes and introns in the mitogenomes of all major land plant lineages. Mitochondrial protein gene content is largely conserved among mosses and liverworts, but it varies substantially among and within other land plant lineages due to convergent losses of genes encoding ribosomal proteins and, to a lesser extent, genes for proteins involved in cytochrome c maturation and oxidative phosphorylation. Mitochondrial intron content is fairly stable within each major land plant lineage, but highly variable among lineages, resulting from occasional gains and many convergent losses over time. Trans-splicing has evolved dozens of times in various vascular plant lineages, particularly those with relatively higher rates of mitogenomic rearrangement. Across eukaryotes, mitochondrial protein gene and intron content has been shaped massive convergent evolution.

Mitogenomics reveals phylogenetic relationships of Arcoida (Mollusca, Bivalvia) and multiple independent expansions and contractions in mitochondrial genome size

Arcoida, comprising about 570 species of blood cockles, is an ecologically and economically important lineage of bivalve molluscs. Current classification of arcoids is largely based on morphology, which shows widespread homoplasy. Despite two recent studies employing multi-locus analyses with broad sampling of Arcoida, evolutionary relationships among major lineages remain controversial. Interestingly, mitochondrial genomes of several ark shell species are 2-3 times larger than those found in most bilaterians, and are among the largest bilaterian mitochondrial genomes reported to date. These results highlight the need of detailed phylogenetic study to explore evolutionary relationships within Arcoida so that the evolution of mitochondrial genome size can be understood. To this end, we sequenced 17 mitochondrial genomes and compared them with publicly available data, including those from other lineages of Arcoida with emphasis on the subclade Arcoidea species. Our phylogenetic analyses indicate that Noetiidae, Cucullaeidae and Glycymerididae are nested within a polyphyletic Arcidae. Moreover, we find multiple independent expansions and potential contractions of mitochondrial genome size, suggesting that the large mitochondrial genome is not a shared ancestral feature in Arcoida. We also examined tandem repeats and inverted repeats in non-coding regions and investigated the presence of such repeats with relation to genome size variation. Our results suggest that tandem repeats might facilitate intraspecific mitochondrial genome size variation, and that inverted repeats, which could be derived from transposons, might be responsible for mitochondrial genome expansions and contractions. We show that mitochondrial genome size in Arcoida is more dynamic than previously understood and provide insights into evolution of mitochondrial genome size variation in metazoans.

Clues of in vivo nuclear gene regulation by mitochondrial short non-coding RNAs

Gene expression involves multiple processes, from transcription to translation to the mature, functional peptide, and it is regulated at multiple levels. Small RNA molecules are known to bind RNA messengers affecting their fate in the cytoplasm (a process generically termed ‘RNA interference’). Such small regulatory RNAs are well-known to be originated from the nuclear genome, while the role of mitochondrial genome in RNA interference was largely overlooked. However, evidence is growing that mitochondrial DNA does provide the cell a source of interfering RNAs. Small mitochondrial highly transcribed RNAs (smithRNAs) have been proposed to be transcribed from the mitochondrion and predicted to regulate nuclear genes. Here, for the first time, we show in vivo clues of the activity of two smithRNAs in the Manila clam, Ruditapes philippinarum. Moreover, we show that smithRNAs are present and can be annotated in representatives of the three main bilaterian lineages; in some cases, they were already described and assigned to a small RNA category (e.g., piRNAs) given their biogenesis, while in other cases their biogenesis remains unclear. If mitochondria may affect nuclear gene expression through RNA interference, this opens a plethora of new possibilities for them to interact with the nucleus and makes metazoan mitochondrial DNA a much more complex genome than previously thought.

Effects of mitonuclear combination and thermal acclimation on the energetic phenotype

Activity of the oxidative phosphorylation complexes rely on intimately associated subunits encoded by the mitochondrial and nuclear genomes. Given the key role of this system in adenosine triphosphate production, genes from both genomes must coevolve. A combination of northern redbelly dace (Chrosomus eos) or finescale dace (C. neogaeus) mitochondrial genome with a C. eos nuclear genome allows for a close examination of a naturally occurring disruption of mitonuclear coevolution. We, therefore, investigated the combined effect of mitonuclear genotypes, acclimation, and temperature on the activity of enzymes linked with the energy metabolism in a sympatric population of wild type and cybrid. As expected, the activity of the nuclear-encoded citrate synthase was only influenced by temperature while the cytochrome c oxidase (composed of nuclear and mitochondrial subunits from wild type and cybrid individuals) responded differently to temperature. This study provides clear evidence of the extent by which mitonuclear coadaptation could influence aerobic metabolism.

Metabolic remodelling associated with mtDNA: insights into the adaptive value of doubly uniparental inheritance of mitochondria

Mitochondria produce energy through oxidative phosphorylation (OXPHOS), which depends on the expression of both nuclear and mitochondrial DNA (mtDNA). In metazoans, a striking exception from strictly maternal inheritance of mitochondria is doubly uniparental inheritance (DUI). This unique system involves the maintenance of two highly divergent mtDNAs (F- and M-type, 8-40% of nucleotide divergence) associated with gametes, and occasionally coexisting in somatic tissues. To address whether metabolic differences underlie this condition, we characterized the OXPHOS activity of oocytes, spermatozoa, and gills of different species through respirometry. DUI species express different gender-linked mitochondrial phenotypes in gametes and partly in somatic tissues. The M-phenotype is specific to sperm and entails (i) low coupled/uncoupled respiration rates, (ii) a limitation by the phosphorylation system, and (iii) a null excess capacity of the final oxidases, supporting a strong control over the upstream complexes. To our knowledge, this is the first example of a phenotype resulting from direct selection on sperm mitochondria. This metabolic remodelling suggests an adaptive value of mtDNA variations and we propose that bearing sex-linked mitochondria could assure the energetic requirements of different gametes, potentially linking male-energetic adaptation, mitotype preservation and inheritance, as well as resistance to both heteroplasmy and ageing.

Unorthodox features in two venerid bivalves with doubly uniparental inheritance of mitochondria

In animals, strictly maternal inheritance (SMI) of mitochondria is the rule, but one exception (doubly uniparental inheritance or DUI), marked by the transmission of sex-specific mitogenomes, has been reported in bivalves. Associated with DUI is a frequent modification of the mitochondrial cox2 gene, as well as additional sex-specific mitochondrial genes not involved in oxidative phosphorylation. With the exception of freshwater mussels (for 3 families of the order Unionida), these DUI-associated features have only been shown in few species [within Mytilidae (order Mytilida) and Veneridae (order Venerida)] because of the few complete sex-specific mitogenomes published for these orders. Here, we present the complete sex-specific mtDNAs of two recently-discovered DUI species in two families of the order Venerida, Scrobicularia plana (Semelidae) and Limecola balthica (Tellinidae). These species display the largest differences in genome size between sex-specific mitotypes in DUI species (>10 kb), as well as the highest mtDNA divergences (sometimes reaching >50%). An important in-frame insertion (>3.5 kb) in the male cox2 gene is partly responsible for the differences in genome size. The S. plana cox2 gene is the largest reported so far in the Kingdom Animalia. The mitogenomes may be carrying sex-specific genes, indicating that general mitochondrial features are shared among DUI species.

Peptides encoded by noncoding genes: challenges and perspectives

In recent years, noncoding gene (NCG) translation events have been frequently discovered. The resultant peptides, as novel findings in the life sciences, perform unexpected functions of increasingly recognized importance in many fundamental biological and pathological processes. The emergence of these novel peptides, in turn, has advanced the field of genomics while indispensably aiding living organisms. The peptides from NCGs serve as important links between extracellular stimuli and intracellular adjustment mechanisms. These peptides are also important entry points for further exploration of the mysteries of life that may trigger a new round of revolutionary biotechnological discoveries. Insights into NCG-derived peptides will assist in understanding the secrets of life and the causes of diseases, and will also open up new paths to the treatment of diseases such as cancer. Here, a critical review is presented on the action modes and biological functions of the peptides encoded by NCGs. The challenges and future trends in searching for and studying NCG peptides are also critically discussed.

Faraway, so close. The comparative method and the potential of non-model animals in mitochondrial research

Inference from model organisms has been the engine for many discoveries in life science, but indiscriminate generalization leads to oversimplifications and misconceptions. Model organisms and inductive reasoning are irreplaceable: there is no other way to tackle the complexity of living systems. At the same time, it is not advisable to infer general patterns from a restricted number of species, which are very far from being representative of the diversity of life. Not all models are equal. Some organisms are suitable to find similarities across species, other highly specialized organisms can be used to focus on differences. In this opinion piece, we discuss the dominance of the mechanistic/reductionist approach in life sciences and make a case for an enhanced application of the comparative approach to study processes in all their various forms across different organisms. We also enlist some rising animal models in mitochondrial research, to exemplify how non-model organisms can be chosen in a comparative framework. These taxa often do not possess implemented tools and dedicated methods/resources. However, because of specific features, they have the potential to address still unanswered biological questions. Finally, we discuss future perspectives and caveats of the comparative method in the age of 'big data'. This article is part of the theme issue 'Linking the mitochondrial genotype to phenotype: a complex endeavour'.

Linking the mitochondrial genotype to phenotype: a complex endeavour

Finding causal links between genotype and phenotype is a major issue in biology, even more in mitochondrial biology. First of all, mitochondria form complex networks, undergoing fission and fusion and we do not know how such dynamics influence the distribution of mtDNA variants across the mitochondrial network and how they affect the phenotype. Second, the non-Mendelian inheritance of mitochondrial genes can have sex-specific effects and the mechanism of mitochondrial inheritance is still poorly understood, so it is not clear how selection and/or drift act on mtDNA genetic variation in each generation. Third, we still do not know how mtDNA expression is regulated; there is growing evidence for a convoluted mechanism that includes RNA editing, mRNA stability/turnover, post-transcriptional and post-translational modifications. Fourth, mitochondrial activity differs across species as a result of several interacting processes such as drift, adaptation, genotype-by-environment interactions, mitonuclear coevolution and epistasis. This issue will cover several aspects of mitochondrial biology along the path from genotype to phenotype, and it is subdivided into four sections focusing on mitochondrial genetic variation, on the relationship among mitochondria, germ line and sex, on the role of mitochondria in adaptation and phenotypic plasticity, and on some future perspectives in mitochondrial research. This article is part of the theme issue 'Linking the mitochondrial genotype to phenotype: a complex endeavour'.

The conserved mitochondrial genomes of Drosophila mercatorum (Diptera: Drosophilidae) with different reproductive modes and phylogenetic implications

Fruit flies (Drosophilidae: Drosophila) are commonly found in daily life and have long been used as model organisms in biology researches. Drosophila mercatorum is one important member of the Drosophila genus and has been used to study centrosome assembly of cells. In this study, we sequenced and analyzed the mitochondrial genome (mitogenome) of D. mercatorum, finding that it contains the typical structure of 37 genes and a control region. The arrangement of mitochondrial genes is in accordance with that in other Drosophila species, which is considered the ancestral organization of insects' mitogenomes. Phylogenetic analyses were performed based on 23 species of Drosophila. Our results supported two monophyletic subgenera, Drosophila and Sophophora, except for D. willistoni which was presented as an early offshoot of Drosophila. The topology ((D. yakuba + D. erecta) + D. melanogaster) was supported. We further compared the mitogenomes of parthenogenesis and sexual reproduction strains of D. mercatorum. However, only one synonymous mutation in COI gene was identified, indicating mitogenomic evolution is not strongly correlated with the different reproductive modes of this species. Taken together, our results demonstrate that mitogenome is an effective molecular marker that can be further used in phylogenetic studies of Drosophila and other organisms.

Putative Mitochondrial Sex Determination in the Bivalvia: Insights From a Hybrid Transcriptome Assembly in Freshwater Mussels

Bivalves exhibit an astonishing diversity of sexual systems, with genetic and environmental determinants of sex, and possibly the only example of mitochondrial genes influencing sex determination pathways in animals. In contrast to all other animal species in which strict maternal inheritance (SMI) of mitochondria is the rule, bivalves possess a system known as doubly uniparental inheritance (DUI) of mitochondria in which maternal and paternal mitochondria (and their corresponding female-transmitted or F mtDNA and male-transmitted or M mtDNA genomes) are transmitted within a species. Species with DUI also possess sex-associated mtDNA-encoded proteins (in addition to the typical set of 13), which have been hypothesized to play a role in sex determination. In this study, we analyzed the sex-biased transcriptome in gonads of two closely-related freshwater mussel species with different reproductive and mitochondrial transmission modes: the gonochoric, DUI species, Utterbackia peninsularis, and the hermaphroditic, SMI species, Utterbackia imbecillis. Through comparative analysis with other DUI and non-DUI bivalve transcriptomes already available, we identify common male and female-specific genes, as well as SMI and DUI-related genes, that are probably involved in sex determination and mitochondrial inheritance in this animal group. Our results contribute to the understanding of what could be the first animal sex determination system involving the mitochondrial genome.

Comparative genomic analysis identifies small open reading frames (sORFs) with peptide-encoding features in avian 16S rDNA

The mitochondrial genome (mt-DNA) functional repertoire has recently been enriched in mammals by the identification of functional small open reading frames (sORFs) embedded in ribosomal DNAs. Through comparative genomic analyses the presence of putatively functional sORFs was investigated in birds. Alignment of available avian mt-DNA sequences revealed highly conserved regions containing four putative sORFs that presented low insertion/deletion polymorphism rate (<0.1%) and preserved in frame start/stop codons in >80% of species. Detected sORFs included avian homologs of human Humanin and Short-Humanin-Like-Peptide 6 and two new sORFs not yet described in mammals. The amino-acid sequences of the four putative encoded peptides were strongly conserved among birds, with amino-acid p-distances (5.6 to 25.4%) similar to those calculated for typical avian mt-DNA-encoded proteins (14.8%). Conservation resulted from either drastic conservation of the nucleotide sequence or negative selection pressure. These data extend to birds the possibility that mitochondrial rDNA may encode small bioactive peptides.

Computational Characterization of the mtORF of Pocilloporid Corals: Insights into Protein Structure and Function in Stylophora Lineages from Contrasting Environments

More than a decade ago, a new mitochondrial Open Reading Frame (mtORF) was discovered in corals of the family Pocilloporidae and has been used since then as an effective barcode for these corals. Recently, mtORF sequencing revealed the existence of two differentiated Stylophora lineages occurring in sympatry along the environmental gradient of the Red Sea (18.5°C to 33.9°C). In the endemic Red Sea lineage RS_LinB, the mtORF and the heat shock protein gene hsp70 uncovered similar phylogeographic patterns strongly correlated with environmental variations. This suggests that the mtORF too might be involved in thermal adaptation. Here, we used computational analyses to explore the features and putative function of this mtORF. In particular, we tested the likelihood that this gene encodes a functional protein and whether it may play a role in adaptation. Analyses of full mitogenomes showed that the mtORF originated in the common ancestor of Madracis and other pocilloporids, and that it encodes a transmembrane protein differing in length and domain architecture among genera. Homology-based annotation and the relative conservation of metal-binding sites revealed traces of an ancient hydrolase catalytic activity. Furthermore, signals of pervasive purifying selection, lack of stop codons in 1830 sequences analyzed, and a codon-usage bias similar to that of other mitochondrial genes indicate that the protein is functional, i.e., not a pseudogene. Other features, such as intrinsically disordered regions, tandem repeats, and signals of positive selection particularly in StylophoraRS_LinB populations, are consistent with a role of the mtORF in adaptive responses to environmental changes.

[The underestimated coding potential of mitochondrial DNA]

Mitochondria are ancient organelles that emerged from the endosymbiosis of free-living proto-bacteria. They still retain a semi-autonomous genetic system with a small genome. Mitochondrial DNA (mtDNA) codes for 13 essential proteins for the production of ATP, the sequences of which are relatively conserved across Metazoans. The discovery of additional mitochondria-derived peptides (MDPs) indicates an underestimated coding potential. Humanin, an anti-apoptotic peptide, is likely independently transcribed from within the 16S rRNA gene, as are recently described SHLPs. MOTS-c, discovered in silico, has been demonstrated to be involved in metabolism and insulin sensitivity. Gau, is a positionally conserved open reading frame (ORF) sequence found in the antisense strand of the COX1 gene and its corresponding peptide is strictly colocalized with mitochondrial markers. In bivalves with doubly uniparental inheritance of mtDNA, male and female mtDNAs each carry a separate additional gene possibly involved in sex determination. Other MDPs likely exist and their investigation will shed light on the underestimated functional repertoire of mitochondria.

Structure-Related Differences between Cytochrome Oxidase I Proteins in a Stable Heteroplasmic Mitochondrial System

Many bivalve species have two types of mitochondrial DNA passed independently through the female line (F genome) and male line (M genome). Here we study the cytochrome oxidase I protein in such bivalve species and provide evidence for differences between the F and M proteins in amino acid property values, particularly relating to hydrophobicity and helicity. The magnitude of these differences varies between different regions of the protein and the change from the ancestor is most marked in the M protein. The observed changes occur in parallel and in the same direction in the different species studied. Two possible causes are considered, first relaxation of purifying selection with drift and second positive selection. These may operate in different ways in different regions of the protein. Many different amino acid substitutions contribute in a small way to the observed variation, but substitutions involving alanine and serine have a quantitatively large effect. Some of these substitutions are potential targets for phosphorylation and some are close to residues of functional importance in the catalytic mechanism. We propose that the observed changes in the F and M proteins might contribute to functional differences between them relating to ATP production and mitochondrial membrane potential with implications for sperm function.

The complete mitochondrial genome of the hermaphroditic freshwater mussel Anodonta cygnea (Bivalvia: Unionidae): in silico analyses of sex-specific ORFs across order Unionoida

BACKGROUND

Doubly uniparental inheritance (DUI) of mitochondrial DNA in bivalves is a fascinating exception to strictly maternal inheritance as practiced by all other animals. Recent work on DUI suggests that there may be unique regions of the mitochondrial genomes that play a role in sex determination and/or sexual development in freshwater mussels (order Unionoida). In this study, one complete mitochondrial genome of the hermaphroditic swan mussel, Anodonta cygnea, is sequenced and compared to the complete mitochondrial genome of the gonochoric duck mussel, Anodonta anatina. An in silico assessment of novel proteins found within freshwater bivalve species (known as F-, H-, and M-open reading frames or ORFs) is conducted, with special attention to putative transmembrane domains (TMs), signal peptides (SPs), signal cleavage sites (SCS), subcellular localization, and potential control regions. Characteristics of TMs are also examined across freshwater mussel lineages.

RESULTS

In silico analyses suggests the presence of SPs and SCSs and provides some insight into possible function(s) of these novel ORFs. The assessed confidence in these structures and functions was highly variable, possibly due to the novelty of these proteins. The number and topology of putative TMs appear to be maintained among both F- and H-ORFs, however, this is not the case for M-ORFs. There does not appear to be a typical control region in H-type mitochondrial DNA, especially given the loss of tandem repeats in unassigned regions when compared to F-type mtDNA.

CONCLUSION

In silico analyses provides a useful tool to discover patterns in DUI and to navigate further in situ analyses related to DUI in freshwater mussels. In situ analysis will be necessary to further explore the intracellular localizations and possible role of these open reading frames in the process of sex determination in freshwater mussel.

Evolution as an ecosystem process: insights from genomics

Evolution is a fundamental ecosystem process. The study of genomic variation of organisms can not only improve our understanding of evolutionary processes, but also of contemporary and future ecosystem dynamics. We argue that integrative research between the fields of genomics and ecosystem ecology could generate new insights. Specifically, studies of biodiversity and ecosystem functioning, evolutionary rescue, and eco-evolutionary dynamics could all benefit from information about variation in genome structure and the genetic architecture of traits, whereas genomic studies could benefit from information about the ecological context of evolutionary dynamics. We propose new ways to help link research on functional genomic diversity with (reciprocal) interactions between phenotypic evolution and ecosystem change. Despite numerous challenges, we anticipate that the wealth of genomic data being collected on natural populations will improve our understanding of ecosystems.

Bijective codon transformations show genetic code symmetries centered on cytosine's coding properties

Homology of some RNAs with template DNA requires systematic exchanges between nucleotides. Such exchanges produce 'swinger' RNA along 23 bijective transformations (nine symmetric, X ↔ Y; and 14 asymmetric, X → Y → Z → X, for example A ↔ C and A → C → G → A, respectively). Here, analyses compare amino acids coded by swinger-transformed codons to those coded by untransformed codons, defining coding invariance after transformations. Swinger transformations cluster according to coding invariance in four groups characterized by transformations into cytosine (C = C, T → C, A → C, and G → C). C's central mutational coding role shows that swinger transformations constrained genetic code genesis. Coding invariance post-transformations correlate positively/negatively with mitochondrial swinger transcription/lepidosaurian body temperature. Presumably, low/high temperatures stabilize/revert rare swinger polymerization modes, producing long swinger sequences/point mutations, respectively. Coding invariance after swinger transformations might compensate effects of swinger polymerizations in species with low body temperatures. Hypothetically, swinger transcription increased coding potential of RNA self-replicating protolife systems under heating/cooling cycles.

The complete mitochondrial genome of parasitic nematode Camallanus cotti: extreme discontinuity in the rate of mitogenomic architecture evolution within the Chromadorea class

BACKGROUND

Complete mitochondrial genomes are much better suited for the taxonomic identification and phylogenetic studies of nematodes than morphology or traditionally-used molecular markers, but they remain unavailable for the entire Camallanidae family (Chromadorea). As the only published mitogenome in the Camallanina suborder (Dracunculoidea superfamily) exhibited a unique gene order, the other objective of this research was to study the evolution of mitochondrial architecture in the Spirurida order. Thus, we sequenced the complete mitogenome of the Camallanus cotti fish parasite and conducted structural and phylogenomic comparative analyses with all available Spirurida mitogenomes.

RESULTS

The mitogenome is exceptionally large (17,901 bp) among the Chromadorea and, with 46 (pseudo-) genes, exhibits a unique architecture among nematodes. Six protein-coding genes (PCGs) and six tRNAs are duplicated. An additional (seventh) tRNA (Trp) was probably duplicated by the remolding of tRNA-Ser2 (missing). Two pairs of these duplicated PCGs might be functional; three were incomplete and one contained stop codons. Apart from Ala and Asp, all other duplicated tRNAs are conserved and probably functional. Only 19 unique tRNAs were found. Phylogenomic analysis included Gnathostomatidae (Spirurina) in the Camallanina suborder.

CONCLUSIONS

Within the Nematoda, comparable PCG duplications were observed only in the enoplean Mermithidae family, but those result from mitochondrial recombination, whereas characteristics of the studied mitogenome suggest that likely rearrangement mechanisms are either a series of duplications, transpositions and random loss events, or duplication, fragmentation and subsequent reassembly of the mitogenome. We put forward a hypothesis that the evolution of mitogenomic architecture is extremely discontinuous, and that once a long period of stasis in gene order and content has been punctuated by a rearrangement event, such a destabilised mitogenome is much more likely to undergo subsequent rearrangement events, resulting in an exponentially accelerated evolutionary rate of mitogenomic rearrangements. Implications of this model are particularly important for the application of gene order similarity as an additive source of phylogenetic information. Chromadorean nematodes, and particularly Camallanina clade (with C. cotti as an example of extremely accelerated rate of rearrangements), might be a good model to further study this discontinuity in the dynamics of mitogenomic evolution.