Comparative Genomics of Marine Mussels (Mytilus spp.) Gender Associated mtDNA: Rapidly Evolving atp8

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.

SNP discovery and genetic structure in blue mussel species using low coverage sequencing and a medium density 60 K SNP-array

Blue mussels from the genus Mytilus are an abundant component of the benthic community, found in the high latitude habitats. These foundation species are relevant to the aquaculture industry, with over 2 million tonnes produced globally each year. Mussels withstand a wide range of environmental conditions and species from the Mytilus edulis complex readily hybridize in regions where their distributions overlap. Significant effort has been made to investigate the consequences of environmental stress on mussel physiology, reproductive isolation, and local adaptation. Yet our understanding on the genomic mechanisms underlying such processes remains limited. In this study, we developed a multi species medium-density 60 K SNP-array including four species of the Mytilus genus. SNPs included in the platform were called from 138 mussels from 23 globally distributed mussel populations, sequenced using a whole-genome low coverage approach. The array contains polymorphic SNPs which capture the genetic diversity present in mussel populations thriving across a gradient of environmental conditions (~59 K SNPs) and a set of published and validated SNPs informative for species identification and for diagnosis of transmissible cancer (610 SNPs). The array will allow the consistent genotyping of individuals, facilitating the investigation of ecological and evolutionary processes in these taxa. The applications of this array extend to shellfish aquaculture, contributing to the optimization of this industry via genomic selection of blue mussels, parentage assignment, inbreeding assessment and traceability. Further applications such as genome wide association studies (GWAS) for key production traits and those related to environmental resilience are especially relevant to safeguard aquaculture production under climate change.

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.

Full Mitochondrial Genomes Reveal Species Differences between the Venerid Clams Ruditapes philippinarum and R. variegatus

In natural sea areas along the coast of China, venerid clams Ruditapes philippinarum and R. variegatus exhibit similar adult shell forms and are especially difficult to distinguish as spat and juveniles. This study used comparative mitochondrial genome analysis to reveal differences between these species. The results showed that: (1) the mitochondrial genomes of R. philippinarum and R. variegatus share a large number of similar gene clusters arranged in consistent order, yet they also display noncommon genes, with both gene rearrangements and random losses found; (2) the 13 protein-coding genes in R. philippinarum as well as two-fold and four-fold degenerate sites in R. variegatus have an evident AT bias; (3) the Ka/Ks ratio of the mitochondrial ATP8 gene was significantly higher in R. philippinarum than in R. variegatus, and an analysis of selection pressure revealed that the mitochondrial NADH dehydrogenase subunit 2 gene and NADH dehydrogenase subunit 6 gene of R. variegatus were under great selective pressure during its evolution; and finally, (4) the two species clustered into one branch on a phylogenetic tree, further affirming their phylogenetic closeness. Based on these results, we speculate that the species differences between R. variegatus and R. philippinarum are largely attributable to adaptive evolution to the environment. The present findings provide a reference for the development of germplasm identification.

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.

Confirmation of the first intronic sequence in the bivalvian mitochondrial genome of Macoma balthica (Linnaeus, 1758)

In 2020, the first male-type mitochondrial genome from the clam Macoma balthica was published. Apart from the unusual doubly uniparental inheritance of mtDNA, scientists observed a unique (over 4k bp long) extension in the middle of the cox2 gene. We have attempted to replicate these data by NGS DNA sequencing and explore further the expression of the long cox2 gene. In our study, we report an even longer cox2 gene (over 5.5 kbp) with no stop codon separating conserved cox2 domains, as well as, based on the rtPCR, a lower relative gene expression pattern of the middle part of the gene (5' = 1; mid = 0.46; 3' = 0.89). Lastly, we sequenced the cox2 gene transcript proving the excision of the intronic sequence.

Complete mitochondrial genome of Pectocera sp. (Elateridae: Dendrometrinae: Oxynopterini) and its phylogenetic implications

The Elateridae family (click beetles) represents a highly diverse lineage that possesses a specialized clicking mechanism to startle predators. At present, however, phylogenetic relationships, especially among recognized subfamilies, remain contentious. Mitochondrial genomes (mitogenomes) can help resolve previously intractable phylogenetic relationships using morphological or limited molecular data. Here, we report the complete mitogenome of Pectocera sp. (Elateridae: Dendrometrinae: Oxynopterini), which was 15,962 bp in length and showed a typical gene number and order as most beetle mitogenomes, including 13 protein-coding genes (PCGs), 22 transfer RNA genes, 2 ribosomal RNA genes, and 1 noncoding control region (AT-rich region). Comparative genomic analyses showed a high degree of feature similarity among Pectocera sp. and other click beetles. Evolutionary analysis of all PCGs based on the nonsynonymous/synonymous substitution rate ratio (ω) indicated that cox1 and atp8 exhibited the lowest and highest evolutionary rates, respectively, and that the evolutionary rates of all PCGs, except for cox3, nad2, and nad3, were lower than the average ω of click beetles. Phylogenetic analyses based on concatenated and coalescent approaches indicated that Pectocera sp. was sister to Campsosternus auratus in the same tribe (Oxynopterini) with high support. This study offers insight into the mitogenomic basis of Pectocera sp. and provides an important data resource for exploring the taxonomy, phylogeny, and evolution of click beetles.

Adaptive mitochondrial genome functioning in ecologically different farm-impacted natural seedbeds of the endemic blue mussel Mytilus chilensis

We assessed the adaptive contribution of the mitochondrial genes involved with the respiratory chain and oxidative phosphorylation of the blue mussel Mytilus chilensis, a native and heavily exploited species in the inner sea of Chiloé Island, southern Chile. The assembled mitochondrial transcriptome of individuals from two ecologically different farm-impacted natural seedbeds, Cochamó (41°S) and Yaldad (42°S), represented about 4.5% of the whole de novo transcriptome of the species and showed location and tissue (gills, mantle) specific expression differences in 13 protein-coding mitochondrial genes. The RNA-Seq analysis detected differences in the number of up-regulated mitogenes between individuals from Cochamó (7) and Yaldad (11), some being tissue-specific (ND4L and COX2). However, the analysis did not detect transcripts-per-million (TPM = 0) of ND2 and ND5 in gills and ATP6 in mantle samples from Cochamó. Likewise, for ND6 and ATP8 in any sample. Several monomorphic location-specific mitochondrial genetic variants were detected in samples from Cochamó (78) and Yaldad (207), representing standing genetic variability to optimize mitochondrial functioning under local habitats. Overall, these mitochondrial transcriptomic differences reflect the impact of environmental conditions on the mitochondrial genome functioning and offer new markers to assess the effects on mussel fitness of habitat translocations, a routine industry practice. Likewise, these mitochondrial markers should help monitor and maintain adaptive population differences in this keystone and heavily exploited native species.

No evidence of DUI in the Mediterranean alien species Brachidontes pharaonis (P. Fisher, 1870) despite mitochondrial heteroplasmy

Two genetically different mitochondrial haplogroups of Brachidontes pharaonis (p-distance 6.8%) have been identified in the Mediterranean Sea. This hinted at a possible presence of doubly uniparental inheritance in this species. To ascertain this possibility, we sequenced two complete mitogenomes of Brachidontes pharaonis mussels and performed a qPCR analysis to measure the relative mitogenome copy numbers of both mtDNAs. Despite the presence of two very similar regions composed entirely of repetitive sequences in the two haplogroups, no recombination between mitogenomes was detected. In heteroplasmic individuals, both mitogenomes were present in the generative tissues of both sexes, which argues against the presence of doubly uniparental inheritance in this species.

The First Eight Mitogenomes of Leaf-Mining Dactylispa Beetles (Coleoptera: Chrysomelidae: Cassidinae) Shed New Light on Subgenus Relationships

The taxonomic classification of Dactylispa, a large genus of leaf-mining beetles, is problematic because it is currently based on morphology alone. Here, the first eight mitochondrial genomes of Dactylispa species, which were used to construct the first molecular phylogenies of this genus, are reported. The lengths of the eight mitogenomes range from 17,189 bp to 20,363 bp. All of the mitochondrial genomes include 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs), and 1 A + T-rich region. According to the nonsynonymous/synonymous mutation ratio (Ka/Ks) of all PCGs, the highest and the lowest evolutionary rates were found for atp8 and cox1, respectively, which is a common phenomenon among animals. According to relative synonymous codon usage, UUA(L) has the highest frequency. With two Gonophorini species as the outgroup, mitogenome-based phylogenetic trees of the eight Dactylispa species were constructed using maximum likelihood (ML) and Bayesian inference (BI) methods based on the PCGs, tRNAs, and rRNAs. Two DNA-based phylogenomic inferences and one protein-based phylogenomic inference support the delimitation of the subgenera Dactylispa s. str. and Platypriella as proposed in the system of Chen et al. (1986). However, the subgenus Triplispa is not recovered as monophyletic. The placement of Triplispa species requires further verification and testing with more species. We also found that both adult body shape and host plant relationship might explain the subgeneric relationships among Dactylispa beetles to a certain degree.

Mitogenomic phylogeny and fossil-calibrated mutation rates for all F- and M-type mtDNA genes of the largest freshwater mussel family, the Unionidae (Bivalvia)

The Unionidae represent an excellent model taxon for unravelling the drivers of freshwater diversity, but, phylogeographic studies on Southeast Asian taxa are hampered by lack of a comprehensive phylogeny and mutation rates for this fauna. We present complete female- (F) and male-type (M) mitogenomes of four genera of the Southeast Asian clade Contradentini+Rectidentini. We calculate substitution rates for the mitogenome, the 13 protein-coding genes, the two ribosomal units and three commonly used fragments (co1, nd1 and 16S) of both F- and M-mtDNA, based on a fossil-calibrated, mitogenomic phylogeny of the Unionidae. Phylogenetic analyses, including an M+F concatenated dataset, consistently recovers a monophyletic Gonideinae. Subfamily-level topology is congruent with that of a previous nuclear genomic study and with patterns in mitochondrial gene order, suggesting Unionidae F-type 2 as a synapomorphy of the Gonideinae. Our phylogeny indicates that the clades Contradentini+Rectidentini and Lamprotulini+Pseudodontini+Gonideini split in the early Cretaceous (~125 Mya), and that the crown group of Contradentini+Rectidentini originated in the late Cretaceous (~79 Mya). Most gonideine tribes originated during the early Palaeogene. Substitution rates were comparable to those previously published for F-type co1 and 16S for certain Unionidae and Margaritiferidae species (pairs).

Semimytilus algosus: first known hermaphroditic mussel with doubly uniparental inheritance of mitochondrial DNA

Doubly uniparental inheritance (DUI) of mitochondrial DNA is a rare phenomenon occurring in some freshwater and marine bivalves and is usually characterized by the mitochondrial heteroplasmy of male individuals. Previous research on freshwater Unionida mussels showed that hermaphroditic species do not have DUI even if their closest gonochoristic counterparts do. No records showing DUI in a hermaphrodite have ever been reported. Here we show for the first time that the hermaphroditic mussel Semimytilus algosus (Mytilida), very likely has DUI, based on the complete sequences of both mitochondrial DNAs and the distribution of mtDNA types between male and female gonads. The two mitogenomes show considerable divergence (34.7%). The presumably paternal M type mitogenome dominated the male gonads of most studied mussels, while remaining at very low or undetectable levels in the female gonads of the same individuals. If indeed DUI can function in the context of simultaneous hermaphroditism, a change of paradigm regarding its involvement in sex determination is needed. It is apparently associated with gonadal differentiation rather than with sex determination in bivalves.

Complete Mitogenome of a Leaf-Mining Buprestid Beetle, Trachys auricollis, and Its Phylogenetic Implications

A complete mitogenome of Trachys auricollis is reported, and a mitogenome-based phylogenetic tree of Elateriformia with all protein-coding genes (PCGs), rRNAs, and tRNAs is presented for the first time. The complete mitochondrial genome of T. auricollis is 16,429 bp in size and contains 13 PCGs, two rRNA genes, 22 tRNA genes, and an A + T-rich region. The A + T content of the entire genome is approximately 71.1%, and the AT skew and GC skew are 0.10 and -0.20, respectively. According to the the nonsynonymous substitution rate to synonymous substitution rates (Ka/Ks) of all PCGs, the highest and lowest evolutionary rates were observed for atp8 and cox1, respectively, which is a common finding among animals. The start codons of all PCGs are the typical ATN. Ten PCGs have complete stop codons, but three have incomplete stop codons with T or TA. As calculated based on the relative synonymous codon usage (RSCU) values, UUA(L) is the codon with the highest frequency. Except for trnS1, all 22 tRNA genes exhibit typical cloverleaf structures. The A + T-rich region of T. auricollis is located between rrnS and the trnI-trnG-trnM gene cluster, with six 72-bp tandem repeats. Both maximum likelihood (ML) and Bayesian (BI) trees suggest that Buprestoidea is close to Byrrhoidea and that Buprestoidea and Byrrhoidea are sister groups of Elateroidea, but the position of Psephenidae is undetermined. The inclusion of tRNAs might help to resolve the phylogeny of Coleoptera.

Sperm competitive advantage of a rare mitochondrial haplogroup linked to differential expression of mitochondrial oxidative phosphorylation genes

Maternal inheritance of mitochondria creates a sex-specific selective sieve through which mitochondrial mutations harmful to males but not females accumulate and contribute to sexual differences in longevity and disease susceptibility. Because eggs and sperm are under disruptive selection, sperm are predicted to be particularly vulnerable to the genetic load generated by maternal inheritance, yet evidence for mitochondrial involvement in male fertility is limited and controversial. Here, we exploit the coexistence of two divergent mitochondrial haplogroups (A and B2) in a Neotropical arachnid to investigate the role of mitochondria in sperm competition. DNA profiling demonstrated that B2-carrying males sired more than three times as many offspring in sperm competition experiments than A males, and this B2 competitive advantage cannot be explained by female mitochondrial haplogroup or male nuclear genetic background. RNA-Seq of testicular tissues implicates differential expression of mitochondrial oxidative phosphorylation (OXPHOS) genes in the B2 competitive advantage, including a 22-fold upregulation of atp8 in B2 males. Previous comparative genomic analyses have revealed functionally significant amino acid substitutions in differentially expressed genes, indicating that the mitochondrial haplogroups differ not only in expression but also in DNA sequence and protein functioning. However, mitochondrial haplogroup had no effect on sperm number or sperm viability, and, when females were mated to a single male, neither male haplogroup, female haplogroup nor the interaction between male/female haplogroup significantly affected female reproductive success. Our findings therefore suggest that mitochondrial effects on male reproduction may often go undetected in noncompetitive contexts and may prove more important in nature than is currently appreciated.

Mitogenomics of Perumytilus purpuratus (Bivalvia: Mytilidae) and its implications for doubly uniparental inheritance of mitochondria

Animal mitochondria are usually inherited through the maternal lineage. The exceptional system allowing fathers to transmit their mitochondria to the offspring exists in some bivalves. Its taxonomic spread is poorly understood and new mitogenomic data are needed to fill the gap. Here, we present for the first time the two divergent mitogenomes from Chilean mussel Perumytilus purpuratus. The existence of these sex-specific mitogenomes confirms that this species has the doubly uniparental inheritance (DUI) of mitochondria. The genetic distance between the two mitochondrial lineages in P. purpuratus is not only much bigger than in the Mytilus edulis species complex but also greater than the distance observed in Musculista senhousia, the only other DUI-positive member of the Mytilidae family for which both complete mitochondrial genomes were published to date. One additional, long ORF (open reading frame) is present exclusively in the maternal mitogenome of P. purpuratus. This ORF evolves under purifying selection, and will likely be a target for future DUI research.

Actively transcribed and expressed atp8 gene in Mytilus edulis mussels

Background

Animal mitochondrial genomes typically encode 37 genes: 13 proteins, 22 tRNAs and two rRNAs. However, many species represent exceptions to that rule. Bivalvia along with Nematoda and Platyhelminthes are often suspected to fully or partially lack the ATP synthase subunit 8 (atp8) gene. This raises the question as to whether they are really lacking this gene or is this maybe an annotation problem? Among bivalves, Mytilus edulis has been inferred to lack an ATP8 gene since the characterization of its mitochondrial genome in 1992. Even though recent bioinformatic analyses suggested that atp8 is present in Mytilus spp., due to high divergence in predicted amino acid sequences, the existence of a functional atp8 gene in this group remains controversial.

Results

Here we demonstrate that M. edulis mitochondrial open reading frames suggested to be atp8 (in male and female mtDNAs) are actively translated proteins. We also provide evidence that both proteins are an integral part of the ATP synthase complex based on in-gel detection of ATP synthase activity and two-dimensional Blue-Native and SDS polyacrylamide electrophoresis.

Conclusion

Many organisms (e.g., Bivalvia along with Nematoda and Platyhelminthes) are considered to be lacking certain mitochondrial genes often only based on poor similarity between protein coding gene sequences in genetically closed species. In some situations, this may lead to the inference that the ATP8 gene is absent, when it is in fact present, but highly divergent. This shows how important complementary role protein-based approaches, such as those in the present study, can provide to bioinformatic, genomic studies (i.e., ability to confirm the presence of a gene).

Disruption of doubly uniparental inheritance of mitochondrial DNA associated with hybridization area of European Mytilus edulis and Mytilus trossulus in Norway

Doubly uniparental inheritance of mitochondria (DUI) is best known in the blue mussel Mytilus. Under this model, two types of mitochondrial DNA exist: female type (F), transmitted from females to offspring of both genders, and male type (M), transmitted exclusively from males to sons. The mitogenomes are usually highly divergent, but an occasional replacement of a typical M genome by a particular F genome has been postulated to explain reduction of this divergence. Disruption of the DUI model has been reported in hybridization areas. Here, we present a new case of DUI disruption in a hybrid M. trossulus/M. edulis population from the North Sea (Norway). No M haplotypes derived from M. trossulus were identified in this population. Typical M haplotypes derived from M. edulis (ME) were rare. Two F-type haplogroups were found: one derived from M. edulis (FE) and the second derived from M. trossulus (FT). Many haplotypes from the FT group were recombinants, with the male CR sequence coming from the M. trossulus genome (FT1 haplogroup) in contrast to M. edulis CR as in the Baltic. FT1 haplotypes were abundant in the studied population, including homoplasmic females. However, males significantly more often carried these haplotypes; therefore, male heteroplasmy involved the original FE and recombinant FT, indicating that the FT genome undergoes masculinization. Structural similarity of FT1 CR with previously reported, masculinized Baltic haplotypes, which were derived from FE/ME recombination, provides further evidence that CR M-F recombination is a prerequisite for masculinization, also in the context of native M. trossulus mtDNA.

First complete female mitochondrial genome in four bivalve species genus Donax and their phylogenetic relationships within the Veneroida order

BACKGROUND

Four species of the genus Donax (D. semistriatus, D. trunculus, D. variegatus and D. vittatus) are common on Iberian Peninsula coasts. Nevertheless, despite their economic importance and overexploitation, scarce genetic resources are available. In this work, we newly determined the complete mitochondrial genomes of these four representatives of the family Donacidae, with the aim of contributing to unveil phylogenetic relationships within the Veneroida order, and of developing genetic markers being useful in wedge clam identification and authentication, and aquaculture stock management.

PRINCIPAL FINDINGS

The complete female mitochondrial genomes of the four species vary in size from 17,044 to 17,365 bp, and encode 13 protein-coding genes (including the atp8 gene), 2 rRNAs and 22 tRNAs, all located on the same strand. A long non-coding region was identified in each of the four Donax species between cob and cox2 genes, presumably corresponding to the Control Region. The Bayesian and Maximum Likelihood phylogenetic analysis of the Veneroida order indicate that all four species of Donax form a single clade as a sister group of other bivalves within the Tellinoidea superfamily. However, although Tellinoidea is actually monophyletic, none of its families are monophyletic.

CONCLUSIONS

Sequencing of complete mitochondrial genomes provides highly valuable information to establish the phylogenetic relationships within the Veneroida order. Furthermore, we provide here significant genetic resources for further research and conservation of this commercially important fishing resource.

Complete female mitochondrial genome of Mytilus chilensis

The controversy surrounding the origin of antitropical distribution of Mytilus mussels and the taxonomic status of southern hemisphere populations remain unsolved, despite the efforts. One of the limiting factors remains the lack of the complete sequences of the representative mitochondrial genomes which would allow their proper comparison with the relatively well-represented northern hemisphere congeneric mussels. To fill this gap we sequenced the representative maternal (F) genome of a native Chilean mussel. The genome is 16,748bp long and structurally identical to the northern hemisphere M. edulis and M. galloprovincialis F genomes. However, the genetic distance from them (≈5%) is twice as high as the maximum distance between them (<2.5%). Thus, the notion that the name M. chilensis should be used for native Chilean Mytilus mussels, with the same rank as M. galloprovincialis and M. edulis is supported.

Sequence motifs associated with paternal transmission of mitochondrial DNA in the horse mussel, Modiolus modiolus (Bivalvia: Mytilidae)

In the majority of metazoans paternal mitochondria represent evolutionary dead-ends. In many bivalves, however, this paradigm does not hold true; both maternal and paternal mitochondria are inherited. Herein, we characterize maternal and paternal mitochondrial control regions of the horse mussel, Modiolus modiolus (Bivalvia: Mytilidae). The maternal control region is 808bp long, while the paternal control region is longer at 2.3kb. We hypothesize that the size difference is due to a combination of repeated duplications within the control region of the paternal mtDNA genome, as well as an evolutionarily ancient recombination event between two sex-associated mtDNA genomes that led to the insertion of a second control region sequence in the genome that is now transmitted via males. In a comparison to other mytilid male control regions, we identified two evolutionarily Conserved Motifs, CMA and CMB, associated with paternal transmission of mitochondrial DNA. CMA is characterized by a conserved purine/pyrimidine pattern, while CMB exhibits a specific 13bp nucleotide string within a stem and loop structure. The identification of motifs CMA and CMB in M. modiolus extends our understanding of Sperm Transmission Elements (STEs) that have recently been identified as being associated with the paternal transmission of mitochondria in marine bivalves.

The complete maternal and paternal mitochondrial genomes of Unio crassus: Mitochondrial molecular clock and the overconfidence of molecular dating

The availability of a rapidly growing number of complete mitochondrial genome sequences provokes high confidence dating approaches. However, even if the congruence between mitochondrial and nuclear markers is reasonable, the resulting topologies are frequently questionable. The unique opportunity to study the evolutionary history of two independent mitochondrial genomes in one phylogenetic context exists in the freshwater mussels family Unionidae. The two lineages function under doubly uniparental inheritance since well before the emergence of the family. Despite the relatively high number of available complete sequences of maternally inherited genomes, comparative analyses are limited by the small number of sequences of counterpart paternally inherited genomes. We have sequenced for the first time the representative set of five sequences (two maternal and three paternal) from the species Unio crassus. Comparative analysis of the phylogenies reconstructed using relevant mitogenomic data available in GenBank (13 species in total) reveal that single - genome inferences are congruent only if the relaxed clock is assumed.

Comparative Large-Scale Mitogenomics Evidences Clade-Specific Evolutionary Trends in Mitochondrial DNAs of Bivalvia

Despite the figure of complete bivalve mitochondrial genomes keeps growing, an assessment of the general features of these genomes in a phylogenetic framework is still lacking, despite the fact that bivalve mitochondrial genomes are unusual under different aspects. In this work, we constructed a dataset of one hundred mitochondrial genomes of bivalves to perform the first systematic comparative mitogenomic analysis, developing a phylogenetic background to scaffold the evolutionary history of the class' mitochondrial genomes. Highly conserved domains were identified in all protein coding genes; however, four genes (namely, atp6, nad2, nad4L, and nad6) were found to be very divergent for many respects, notwithstanding the overall purifying selection working on those genomes. Moreover, the atp8 gene was newly annotated in 20 mitochondrial genomes, where it was previously declared as lacking or only signaled. Supernumerary mitochondrial proteins were compared, but it was possible to find homologies only among strictly related species. The rearrangement rate on the molecule is too high to be used as a phylogenetic marker, but here we demonstrate for the first time in mollusks that there is correlation between rearrangement rates and evolutionary rates. We also developed a new index (HERMES) to estimate the amount of mitochondrial evolution. Many genomic features are phylogenetically congruent and this allowed us to highlight three main phases in bivalve history: the origin, the branching of palaeoheterodonts, and the second radiation leading to the present-day biodiversity.

In silico analyses of mitochondrial ORFans in freshwater mussels (Bivalvia: Unionoida) provide a framework for future studies of their origin and function

BACKGROUND

Many species of bivalves exhibit a unique system of mtDNA transmission named Doubly Uniparental Inheritance (DUI). Under this system, species have two distinct, sex-linked mitochondrial genomes: the M-type mtDNA, which is transmitted by males to male offspring and found in spermatozoa, and the F-type mtDNA, which is transmitted by females to all offspring, and found in all tissues of females and in somatic tissues of males. Bivalves with DUI also have sex-specific mitochondrial ORFan genes, (M-orf in the M mtDNA, F-orf in the F mtDNA), which are open reading frames having no detectable homology and no known function. DUI ORFan proteins have previously been characterized in silico in a taxonomically broad array of bivalves including four mytiloid, one veneroid and one unionoid species. However, the large evolutionary distance among these taxa prevented a meaningful comparison of ORFan properties among these divergent lineages. The present in silico study focuses on a suite of more closely-related Unionoid freshwater mussel species to provide more reliably interpretable information on patterns of conservation and properties of DUI ORFans. Unionoid species typically have separate sexes, but hermaphroditism also occurs, and hermaphroditic species lack the M-type mtDNA and possess a highly mutated version of the F-orf in their maternally transmitted mtDNA (named H-orf in these taxa). In this study, H-orfs and their respective proteins are analysed for the first time.

RESULTS

Despite a rapid rate of evolution, strong structural and functional similarities were found for M-ORF proteins compared among species, and among the F-ORF and H-ORF proteins across the studied species. In silico analyses suggest that M-ORFs have a role in transport and cellular processes such as signalling, cell cycle and division, and cytoskeleton organisation, and that F-ORFs may be involved in cellular traffic and transport, and in immune response. H-ORFs appear to be structural glycoproteins, which may be involved in signalling, transport and transcription. Our results also support either a viral or a mitochondrial origin for the ORFans.

CONCLUSIONS

Our findings reveal striking structural and functional similarities among proteins encoded by mitochondrial ORFans in freshwater mussels, and strongly support a role for these genes in the DUI mechanism. Our analyses also support the possibility of DUI systems with elements of different sources/origins and different mechanisms of action in the distantly-related DUI taxa. Parallel situations to the novel mitochondrially-encoded functions of freshwater mussel ORFans present in some other eukaryotes are also discussed.

Mitogenomics of southern hemisphere blue mussels (Bivalvia: Pteriomorphia): Insights into the evolutionary characteristics of the Mytilus edulis complex

Marine blue mussels (Mytilus spp.) are widespread species that exhibit an antitropical distribution with five species occurring in the Northern Hemisphere (M. trossulus, M. edulis, M. galloprovincialis, M. californianus and M. coruscus) and three in the Southern Hemisphere (M. galloprovincialis, M. chilensis and M. platensis). Species limits in this group remain controversial, in particular for those forms that live in South America. Here we investigated structural characteristics of marine mussels mitogenomes, based on published F mtDNA sequences of Northern Hemisphere species and two newly sequenced South American genomes, one from the Atlantic M. platensis and another from the Pacific M. chilensis. These mitogenomes exhibited similar architecture to those of other genomes of Mytilus, including the presence of the Atp8 gene, which is missing in most of the other bivalves. Our evolutionary analysis of mitochondrial genes indicates that purifying selection is the predominant force shaping the evolution of the coding genes. Results of our phylogenetic analyses supported the monophyly of Pteriomorphia and fully resolved the phylogenetic relationships among its five orders. Finally, the low genetic divergence of specimens assigned to M. chilensis and M. platensis suggests that these South American marine mussels represent conspecific variants rather than distinct species.

The Complete Female- and Male-Transmitted Mitochondrial Genome of Meretrix lamarckii

Bivalve mitochondrial genomes show many uncommon features, like additional genes, high rates of gene rearrangement, high A-T content. Moreover, Doubly Uniparental Inheritance (DUI) is a distinctive inheritance mechanism allowing some bivalves to maintain and transmit two separate sex-linked mitochondrial genomes. Many bivalve mitochondrial features, such as gene extensions or additional ORFs, have been proposed to be related to DUI but, up to now, this topic is far from being understood. Several species are known to show this unusual organelle inheritance but, being widespread only among Unionidae and Mytilidae, DUI distribution is unclear. We sequenced and characterized the complete female- (F) and male-transmitted (M) mitochondrial genomes of Meretrix lamarckii, which, in fact, is the second species of the family Veneridae where DUI has been demonstrated so far. The two mitochondrial genomes are comparable in length and show roughly the same gene content and order, except for three additional tRNAs found in the M one. The two sex-linked genomes show an average nucleotide divergence of 16%. A 100-aminoacid insertion in M. lamarckii M-cox2 gene was found; moreover, additional ORFs have been found in both F and M Long Unassigned Regions of M. lamarckii. Even if no direct involvement in DUI process has been demonstrated so far, the finding of cox2 insertions and supernumerary ORFs in M. lamarckii both strengthens this hypothesis and widens the taxonomical distribution of such unusual features. Finally, the analysis of inter-sex genetic variability shows that DUI species form two separate clusters, namely Unionidae and Mytilidae+Veneridae; this dichotomy is probably due to different DUI regimes acting on separate taxa.

Mitochondrial genomes of the Baltic clam Macoma balthica (Bivalvia: Tellinidae): setting the stage for studying mito-nuclear incompatibilities

BACKGROUND

Allopatric divergence across lineages can lead to post-zygotic reproductive isolation upon secondary contact and disrupt coevolution between mitochondrial and nuclear genomes, promoting emergence of genetic incompatibilities. A previous F ST scan on the transcriptome of the Baltic clam Macoma balthica highlighted several genes potentially involved in mito-nuclear incompatibilities (MNIs). As proteins involved in the mitochondrial oxidative phosphorylation (OXPHO) chain are prone to MNIs and can contribute to the maintenance of genetic barriers, the mitochondrial genomes of six Ma. balthica individuals spanning two secondary contact zones were sequenced using the Illumina MiSeq plateform.

RESULTS

The mitogenome has an approximate length of 16,806 bp and encodes 13 protein-coding genes, 2 rRNAs and 22 tRNAs, all located on the same strand. atp8, a gene long reported as rare in bivalves, was detected. It encodes 42 amino acids and is putatively expressed and functional. A large unassigned region was identified between rrnS and tRNA (Met) and could likely correspond to the Control Region. Replacement and synonymous mutations were mapped on the inferred secondary structure of all protein-coding genes of the OXPHO chain. The atp6 and atp8 genes were characterized by background levels of replacement mutations, relative to synonymous mutations. However, most nad genes (notably nad2 and nad5) were characterized by an elevated proportion of replacement mutations.

CONCLUSIONS

Six nearly complete mitochondrial genomes were successfully assembled and annotated, providing the necessary roadmap to study MNIs at OXPHO loci. Few replacement mutations were mapped on mitochondrial-encoded ATP synthase subunits, which is in contrast with previous data on nuclear-encoded subunits. Conversely, the high population divergence and the prevalence of non-synonymous mutations at nad genes are congruent with previous observations from the nuclear transcriptome. This further suggest that MNIs between subunits of Complex I of the OXPHO chain, coding for NADH dehydrogenase, may play a role in maintaining barriers to gene flow in Ma. balthica.

Mitogenomics of recombinant mitochondrial genomes of Baltic Sea Mytilus mussels

Recombination in the control region (CR) of Mytilus mitochondrial DNA (mtDNA) was originally reported based on the relatively short, sequenced fragments of mitochondrial genomes. Recombination outside the CR has been reported recently with the suggestion that such processes are common in Mytilus. We have fully sequenced a set of 11 different mitochondrial haplotypes representing the high diversity of paternally inherited mitochondrial genomes of Baltic Sea Mytilus mussels, including the haplotype close to the native Mytilus trossulus mitochondrial genome, which was thought to have been entirely eliminated from this population. Phylogenetic and comparative analysis showed that the recombination is limited to the vicinity of the CR in all sequenced genomes. Coding sequence comparison indicated that all paternally inherited genomes showed increased accumulation of nonsynonymous substitutions, including the genomes which switched their transmission route very recently. The acquisition of certain CR sequences through recombination with highly divergent paternally inherited genomes seems to precede and favor the switch, but it is not a prerequisite for this process. Interspecies hybridization in the Baltic Sea during the recent 10,000 years created conditions for both structural and evolutionary mitochondrial instability which resulted in the observed variation and dynamics of mtDNA in Baltic Sea Mytilus mussels. In conclusion, the data shows that the effects of mitochondrial recombination are limited to the CR of few phylogenetic lineages.

Complete male mitochondrial genome of Anodonta anatina (Mollusca: Unionidae)

Anodonta anatina is a freshwater mussel of the family Unionidae. These mussels have a unique mitochondria inheritance system named doubly uniparental inheritance (DUI). Under DUI males have two, potentially very divergent mitochondrial genomes: F-type inherited from mother and M-type inherited from father. F-type is present in soma whereas M-type is present in gonadal tissues and sperm. Here we report two M-type sequences of complete mitochondrial genomes from Anodonta anatina. They are 16,906 bp long and their sequences are similar (0.1% divergence). The genome organization is identical to the other Unionidean M-type genomes published to date. There are 38 genes, including the recently described M-type specific M ORF. The presence of tRNA-like repeat in one of the noncoding regions, suggests that the control region is located in this area. Nucleotide composition is quite extreme, with AT content (66.2%) higher than in any other of the six published Unionidean M genomes.

Complete male mitochondrial genomes of European Mytilus edulis mussels

Doubly uniparental inheritance (DUI) results in the existence of two gender-specific, divergent mtDNA lineages within a single species. Under DUI, the female genome (F) is transmitted from mothers to the whole offspring, and the male genome (M) is transmitted exclusively from fathers to sons. This system was first described in a marine mussels Mytilus edulis inhabiting European coastal waters, over a decade ago. Despite that, the complete sequence of the M genome from the European M. edulis mussels remained unknown. Here we announce it for the first time. The announcement is based on the two haplotypes isolated from heteroplasmic males of European M. edulis sampled at two moderately distant locations: southern North Sea and western Baltic. The two M genomes are quite similar both in length (16,631 and 16,632 bp) and in sequence (98.3%). Furthermore, both newly sequenced genomes are closely related to the genomes described from Baltic M. trossulus.

Glacial history of the European marine mussels Mytilus, inferred from distribution of mitochondrial DNA lineages

Mussels of the genus Mytilus have been used to assess the circumglacial phylogeography of the intertidal zone. These mussels are representative components of the intertidal zone and have rapidly evolving mitochondrial DNA, suitable for high resolution phylogeographic analyses. In Europe, the three Mytilus species currently share mitochondrial haplotypes, owing to the cases of extensive genetic introgression. Genetic diversity of Mytilus edulis, Mytilus trossulus and Mytilus galloprovincialis was studied using a 900-bp long part of the most variable fragment of the control region from one of their two mitochondrial genomes. To this end, 985 specimens were sampled along the European coasts, at sites ranging from the Black Sea to the White Sea. The relevant DNA fragments were amplified, sequenced and analyzed. Contrary to the earlier findings, our coalescence and nested cladistics results show that only a single M. edulis glacial refugium existed in the Atlantic. Despite that, the species survived the glaciation retaining much of its diversity. Unsurprisingly, M. galloprovincialis survived in the Mediterranean Sea. In a relatively short time period, around the climatic optimum at 10 ky ago, the species underwent rapid expansion coupled with population differentiation. Following the expansion, further contemporary gene flow between populations was limited.

mtDNA sequence, phylogeny and evolution of laboratory mice

Laboratory mice are important tools for biomedical research. Aiming to investigate the phylogeny and evolution of laboratory mice, we investigated the mtDNA sequences of classic inbred strains, classic outbred stocks and wild-derived inbred strains. Our results showed that the most classic outbred stocks and classic inbred strains are descended from a single mtDNA ancestor. The phylogenic analysis supports the topology of M. m. castaneus/M. m. domesticus as sister subspecies, and the divergence time between the two sister subspecies and M. m. musculus was 493,000 (435,000-557,000) years ago. Furthermore, the mtDNA polymorphisms accumulated in the last 100years in the laboratory mice are under a relaxed purifying selection.

Recombinant mitochondrial genome with standard transmission route from Mediterranean mussel Mytilus galloprovincialis

Several bivalve species, including marine mussels Mytilus are atypical in having two gender-specific and highly divergent mtDNA genomes. This peculiar genetic system allows not only the recombination to occur but also facilitates its detection. Previous reports associated the existence of mosaic recombinant haplotypes with the switch of their transmission route. Here we report nearly complete sequence of a mitochondrial genome isolated from a homoplasmic female individual of Mediterranean Mytilus galloprovincialis. The genome has clear phylogenetic affinity with and organization identical to the M. galloprovincialis female haplotypes, in the coding part. However, the genome is very large, approximately 20,600 bp long, exclusively due to a long and complex control region. It contains an array of repeats, some of which are degenerated. A large part of the control region is derived from the paternal genome. This finding shows that not all haplotypes with recombinant control regions must be paternally inherited in Mytilus.

Co-expressed mitochondrial genomes: recently masculinized, recombinant mitochondrial genome is co-expressed with the female-transmitted mtDNA genome in a male Mytilus trossulus mussel from the Baltic Sea

BACKGROUND

Few exceptions have been described from strict maternal inheritance of mitochondrial DNA in animals, including sea mussels (Mytilidae), clams (Donacidae, Veneridae and Solenidae) and freshwater mussels (Unionoidae) order. In these bivalves mitochondria and their DNA are transferred through two separate routes. The females inherit only the maternal mitochondrial DNA whereas the males inherit maternal as well as paternal mitochondrial DNA, which is usually present only in gonads and sperm. The mechanism controlling this phenomenon is unclear but leads to the existence of two separate mitochondrial DNA lineages in a single species. The lineages are usually well differentiated: up to 20-50% divergence in nucleotide sequence. Occasionally, a maternal mitochondrial DNA can invade the paternal transmission route, eventually replacing the diverged M-type and lowering the divergence. Such role reversal (masculinization) event has happened recently in the Mytilus population of the Baltic Sea which consists of M. edulis × M. trossulus hybrids, but the functional status of the resulting mitochondrial genome was unknown.

RESULTS

In this paper we sequenced transcripts from one specimen that was identified as male carrying both the female mitochondrial genome and a recently masculinized mitochondrial genome. Additionally, the analysis of the control region has showed that the recently masculinized, recombinant genome, not only has an M-type control region and all coding regions derived from the F-type, but also is transcriptionally active along side the maternally inherited F-type genome. In the comparative analysis, the two genomes exhibit different substitution patterns, typical for the M vs. F genome comparisons. The genetic distances and ratios of non-synonymous substitutions also suggest that one of the genomes is transitioning from the maternal to the paternal inheritance mode, consistent with its recent masculinization.

CONCLUSION

We have shown, for the first time, that the recently masculinized mitochondrial genome is active and that it accumulates excess of non-synonymous substitutions across its coding sequence. This suggests, that, under certain cytonuclear incompatibility conditions, masculinization may serve to restore the endangered functionality of the paternally inherited genome. This is also another example of a mitochondrial genome in which the recombination in the control region predated its transition from paternal to maternal transmission route.

The rRNA and tRNA transcripts of maternally and paternally inherited mitochondrial DNAs of Mytilus galloprovincialis suggest presence of a "degradosome" in mussel mitochondria and necessitate the re-annotation of the l-rRNA/CR boundary

Species of the genus Mytilus carry two mitochondrial genomes in obligatory coexistence; one transmitted though the eggs (the F type) and one through the sperm (the M type). We have studied the 3' and 5' ends of rRNA and tRNA transcripts using RT-PCR and RNA circularization techniques in both the F and M genomes of Mytilus galloprovincialis. We have found polyadenylated and non-adenylated transcripts for both ribosomal and transfer RNAs. In all these genes the 5' ends of the transcripts coincided with the first nucleotide of the annotated genes, but the 3' ends were heterogeneous. The l-rRNA 3' end is 47 or 48 nucleotides upstream from the one assigned by a previous annotation, which makes the adjacent first domain (variable domain one, VD1) of the main control region (CR) correspondingly longer. We have observed s-rRNA and l-rRNA transcripts with truncated 3' end and polyadenylated tRNA transcripts carrying the CCA trinucleotide. We have also detected polyadenylated RNA remnants carrying the sequences of the control region, which strongly suggests RNA degradation activity and thus presence of degradosomes in Mytilus mitochondria.

20 years since the introduction of DNA barcoding: from theory to application

Traditionally, taxonomic identification has relied upon morphological characters. In the last two decades, molecular tools based on DNA sequences of short standardised gene fragments, termed DNA barcodes, have been developed for species discrimination. The most common DNA barcode used in animals is a fragment of the cytochrome c oxidase (COI) mitochondrial gene, while for plants, two chloroplast gene fragments from the RuBisCo large subunit (rbcL) and maturase K (matK) genes are widely used. Information gathered from DNA barcodes can be used beyond taxonomic studies and will have far-reaching implications across many fields of biology, including ecology (rapid biodiversity assessment and food chain analysis), conservation biology (monitoring of protected species), biosecurity (early identification of invasive pest species), medicine (identification of medically important pathogens and their vectors) and pharmacology (identification of active compounds). However, it is important that the limitations of DNA barcoding are understood and techniques continually adapted and improved as this young science matures.

The complete mitochondrial genome of the land snail Cornu aspersum (Helicidae: Mollusca): intra-specific divergence of protein-coding genes and phylogenetic considerations within Euthyneura

The complete sequences of three mitochondrial genomes from the land snail Cornu aspersum were determined. The mitogenome has a length of 14050 bp, and it encodes 13 protein-coding genes, 22 transfer RNA genes and two ribosomal RNA genes. It also includes nine small intergene spacers, and a large AT-rich intergenic spacer. The intra-specific divergence analysis revealed that COX1 has the lower genetic differentiation, while the most divergent genes were NADH1, NADH3 and NADH4. With the exception of Euhadra herklotsi, the structural comparisons showed the same gene order within the family Helicidae, and nearly identical gene organization to that found in order Pulmonata. Phylogenetic reconstruction recovered Basommatophora as polyphyletic group, whereas Eupulmonata and Pulmonata as paraphyletic groups. Bayesian and Maximum Likelihood analyses showed that C. aspersum is a close relative of Cepaea nemoralis, and with the other Helicidae species form a sister group of Albinaria caerulea, supporting the monophyly of the Stylommatophora clade.

Molecular population genetics of male and female mitochondrial genomes in subarctic Mytilus trossulus

The doubly uniparental inheritance system allows for the use of two independent mitochondrial genomes for population history investigations. Under this system, two lineages of mitochondrial DNA (mtDNA) exist and males are typically heteroplasmic, having the additional, usually divergent, mitochondrial genome inherited from their male parent. This additional mtDNA typically evolves faster, potentially allowing for insight into more recent events in population history. Few studies did explore this possibility in marine mussels Mytilus showing its usefulness. Recent observations of the Mytilus trossulus mussels who have retained their native mtDNA in European waters posed the question of their origin. Are they part of a population present, but previously undetected, or is this a potentially human mediated, ongoing spread of an invasive species? To tackle this question, we amplified with species-specific primers and sequenced an approximately 1,200-bp-long fragment spanning COIII and ND2 genes from both mitochondrial genomes of mussels sampled at five locations worldwide, covering the whole M. trossulus range. The overall pattern of polymorphisms is compatible with the entirely postglacial history of the whole species, indicating a very deep bottleneck at last glacial maximum, with possible retention of the whole species in a single refugium, and the effective population size of no more than a few thousands. Both analyses of molecular variance and isolation with migration (IM) models point at the West Atlantic as the source of the European M. trossulus mussels, at least the ones who retained their native mtDNA. The hypothesis that this is an ongoing, human-mediated process was considered. To this end, comparison with the well-known case: the introduction of congeneric mussel, Mytilus galloprovincialis, from Mediterranean Sea to Asia was used. This introduction occurred within the last 100 years. The results inferred by the IM model suggest that the timing and structure of transatlantic migration of M. trossulus differs significantly from the M. galloprovincialis case: it is more than 1,000 years old and involves a much larger fraction of the ancestral population. Therefore, most likely, this invasion is not a human-mediated process.