Complete mitochondrial genome of the Asian green musselPerna viridis(Bivalvia, Mytilidae)

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.

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.

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.

Population connectivity and genetic structure of Asian green mussel, Perna viridis along Indian waters assessed using mitochondrial markers

Perna viridis (Linnaeus, 1758), the Asian green mussel, belonging to the family Mytilidae is widely distributed along the Indian coast. The species is majorly found in southeastern countries and is considered an ideal candidate for aquaculture due to its high nutritional value and growth rate. Obtaining their genetic information is essential for their sustainable capture-based production. In the present study, genetic variation, population structure, and demographic processes of the populations across the distribution of this species were assessed using the mitochondrial DNA ATPase6 and cytb gene. In total, we selected 170 samples from five localities across the Indian subcontinent including Andaman Sea. Sequence analysis of partial cytb (885 bp) and ATPase6 (714 bp) genes revealed 45 and 58 haplotypes, respectively. The significant coefficient of genetic differentiation (F_ST: 0.255 for cytb and 0.252 for ATPase6) and analyses of molecular variance indicated three varieties of stocks, namely Arabian Sea, Bay of Bengal, and Andaman Sea. All the populations showed low nucleotide diversity, suggesting severe historical bottleneck events and high haplotype diversity, indicating population expansion. The genetic variation and demographic process reported in this study will form the baseline information for framing policies, which can be adopted while planning stock specific ranching and relaying programmes in the Indian subcontinent with view to enhance and manage the fishery.

A mitochondrial genome phylogeny of Mytilidae (Bivalvia: Mytilida)

The family Mytilidae is a family of bivalve mussels that are distributed worldwide in diverse marine habitats. Within the family, classification systems and phylogenetic relationships among subfamilies remain not yet fully resolved. In this study, we newly determined 9 mitochondrial genome sequences from 7 subfamilies: Bathymodiolus thermophilus (Bathymodiolinae), Modiolus nipponicus (Modiolinae), Lithophaga curta (the first representative of Lithophaginae), Brachidontes mutabilis (Brachidontinae), Mytilisepta virgata (Brachidontinae), Mytilisepta keenae (Brachidontinae), Crenomytilus grayanus (Mytilinae), Gregariella coralliophaga (Crenellinae), and Septifer bilocularis (the first representative of Septiferinae). Phylogenetic trees using maximum likelihood and Bayesian inference methods for 28 mitochondrial genomes (including 19 previously published sequences) showed two major clades with high support values: Clade 1 ((Bathymodiolinae + Modiolinae) + (Lithophaginae + Limnoperninae)) and Clade 2 (((Mytilinae + Crenellinae) + Septiferinae) + Brachidontinae). The position of the genus Lithophaga (representing Lithophaginae) differed from a previously published molecular phylogeny. Divergence time analysis with a molecular clock indicated that lineage splitting among the major subfamilies of Mytilidae (including the habitat transition from marine to freshwater environments by ancestral Limnoperninae) occurred in the Mesozoic period, coinciding with high diversification rates of marine fauna during that time. This is the first mitochondrial genome-based phylogenetic study of the Mytilidae that covers nearly all subfamily members, excluding the subfamily Dacrydiinae.

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 mitochondrial genome of the golden mussel Limnoperna fortunei and comparative mitogenomics of Mytilidae

Here we describe the mitochondrial genome of the golden mussel Limnoperna fortunei, an Asian bivalve which has become one of the most aggressive invasive species in Japan and South America. The mitochondrial genome of L. fortunei does not present conserved gene arrangement when compared to the other Mytilidae species suggesting a high degree of gene recombination in the mitochondria of this clade. In addition, the golden mussel mitogenome encodes two copies of tRNA(Lys) and presents a putative pseudogene for the atp8 gene sequence that encodes a 27 amino acid peptide containing an in-frame stop codon. The presence of this pseudogene raises the question as to whether atp8 is encoded in some bivalve mitochondrial genomes or not. The phylogenetic analysis of all complete mitochondrial genomes available from Mytilidae mussels confirmed the close evolutionary relationships among bivalves from the genus Mytilys and placed L. fortunei coming from a more ancestral branch on the family. The supermatrix phylogeny described used the concatenation of all 12 genes from the mitochondria and disputed the monophyly of the genus Perna, as Perna perna was shown to be more closely related to Brachidontes exustus than to Perna viridis. The comparative analysis of mitogenome synteny also confirmed the polyphyly of the genus Perna. The complete and annotated mitogenome has been published in GenBank under the accession number KP756905.

Linking environmental heavy metal concentrations and salinity gradients with metal accumulation and their effects: A case study in 3 mussel species of Vitória estuary and Espírito Santo bay, Southeast Brazil

The present study was conducted to link the heavy metal load in three species of mussels (Perna perna, Mytella falcata and Mytella guyanensis) from the estuaries and bays around Vitória island, south-east of Brazil, with the salinity gradient and the heavy metal levels in the abiotic environment (including water, suspended particulate matter (SPM) and sediment). Primarily based on the salinity gradient, a total of 26 sites around Vitória Island were selected for sampling of water, SPM, sediments and organisms. Besides tissue metal levels, the condition index and energy stores (glycogen, lipid and protein) were quantified as an indicator of fitness in response to metal pollution. Dissolved metals in water indicate that Cd and Mn content was higher along Espírito Santo Bay, while Al, Co, Cu, Cr and Fe were elevated in the sites with low salinity such as river mouths, estuarine and sewage canals. Likewise, suspended matter sampled from low salinity sites showed a higher heavy metal load compared to moderate and high salinity sites. Though mussels were sampled from different sites, the contamination for Cd, Cu, Fe and Mn was higher in mussels inhabiting low salinity sites (M. guyanensis and M. falcata) compared to P. perna, a high saline water inhabitant. However, a higher Zn body burden was observed for P. perna compared to Mytella species. Tissue Fe accumulation (but not Mn and Zn) correlated with heavy metal levels in suspended material for all three species, and for M. falcata this correlation also existed for Cd and Cu. Energy store and condition index in all mussels varied depending on the sampling sites and correlated with salinity gradient rather than tissue metal concentration. Overall, metal concentration in mussels did not exceed the safe levels as per the international standards for metals, and would be of no risk for human consumption.

De novo transcriptome analysis of Perna viridis highlights tissue-specific patterns for environmental studies

BACKGROUND

The tropical green-lipped mussel Perna viridis is a common biomonitor throughout the Indo-Pacific region that is used for environmental monitoring and ecotoxicological investigations. However, there is limited molecular data available regarding this species. We sought to establish a global transcriptome database from the tissues of adductor muscle, gills and the hepatopancreas of P. viridis in an effort to advance our understanding of the molecular aspects involved during specific toxicity responses in this sentinel species.

RESULTS

Illumina sequencing results yielded 544,272,542 high-quality filtered reads. After de novo assembly using Trinity, 233,257 contigs were generated with an average length of 1,264 bp and an N50 length of 2,868 bp; 192,879 assembled transcripts and 150,111 assembled unigenes were obtained after clustering. A total of 93,668 assembled transcripts (66,692 assembled genes) with putative functions for protein domains were predicted based on InterProScan analysis. Based on similarity searches, 44,713 assembled transcripts and 25,319 assembled unigenes were annotated with at least one BLAST hit. A total of 21,262 assembled transcripts (11,947 assembled genes) were annotated with at least one well-defined Gene Ontology (GO) and 5,131 assembled transcripts (3,181 assembled unigenes) were assigned to 329 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. The quantity of assembled unigenes and transcripts obtained from male and female mussels were similar but varied among the three studied tissues, with the highest numbers recorded in the gills, followed by the hepatopancreas, and then the adductor muscle. Multivariate analyses revealed strong tissue-specific patterns among the three different tissues, but not between sexes in terms of expression profiles for annotated genes in various GO terms, and genes associated with stress responses and degradation of xenobiotics. The expression profiles of certain selected genes in each tissue type were further validated using real-time quantitative polymerase chain reaction assays and a similar tissue-specific trend was seen.

CONCLUSIONS

The extensive sequence data generated from this study will provide a valuable molecular resource for facilitating environmental studies with P. viridis, and highlight the importance of tissue-specific approaches in the future.