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

The effect of mitochondrial recombination on fertilization success in blue mussels

The presence of doubly uniparental inheritance (DUI) in bivalves represents a unique mode of mitochondrial transmission, whereby paternal (male-transmitted M-type) and maternal (female-transmitted F-type) haplotypes are transmitted to offspring separately. Male embryos retain both haplotypes, but the M-type is selectively removed from females. Due to the presence of heteroplasmy in males, mtDNA can recombine resulting in a 'masculinized' haplotype referred to as Mf-type. While mtDNA recombination is usually rare, it has been recorded in multiple mussel species across the Northern Hemisphere. Given that mitochondria are the powerhouse of the cell, different mtDNA haplotypes may have different selective advantages under diverse environmental conditions. This may be particularly important for sperm fitness and fertilization success. In this study we aimed to i) determine the presence, prevalence of the Mf-type in Australian blue mussels (Mytilus sp.) and ii) investigate the effect of Mf-mtDNA on sperm performance (a fitness correlate). We found a high prevalence of recombined mtDNA (≈35 %) located within the control region of the mitochondrial genome, which occurred only in specimens that contained Southern Hemisphere mtDNA. The presence of two female mitotypes were identified in the studied mussels, one likely originating from the Northern Hemisphere, and the other either representing the endemic M. planulatus species or introduced genotypes from the Southern Hemisphere. Despite having recombination events present in a third of the studied population, analysis of sperm performance indicated no difference in fertilization success related to mitotype.

Epigenetic variation mediated by lncRNAs accounts for adaptive genomic differentiation of the endemic blue mussel Mytiluschilensis

Epigenetic variation affects gene expression without altering the underlying DNA sequence of genes controlling ecologically relevant phenotypes through different mechanisms, one of which is long non-coding RNAs (lncRNAs). This study identified and evaluated the gene expression of lncRNAs in the gill and mantle tissues of Mytilus chilensis individuals from two ecologically different sites: Cochamó (41°S) and Yaldad (43°S), southern Chile, both impacted by climatic-related conditions and by mussel farming given their use as seedbeds. Sequences identified as lncRNAs exhibited tissue-specific differences, mapping to 3.54 % of the gill transcriptome and 1.96 % of the mantle transcriptome, representing an average of 2.76 % of the whole transcriptome. Using a high fold change value (≥|100|), we identified 43 and 47 differentially expressed lncRNAs (DE-lncRNAs) in the gill and mantle tissue of individuals sampled from Cochamó and 21 and 17 in the gill and mantle tissue of individuals sampled from Yaldad. Location-specific DE-lncRNAs were also detected in Cochamó (65) and Yaldad (94) samples. Via analysis of the differential expression of neighboring protein-coding genes, we identified enriched GO terms related to metabolic, genetic, and environmental information processing and immune system functions, reflecting how the impact of local ecological conditions may influence the M. chilensis (epi)genome expression. These DE-lncRNAs represent complementary biomarkers to DNA sequence variation for maintaining adaptive differences and phenotypic plasticity to cope with natural and human-driven perturbations.

Screening of genes encoding proteins that interact with Nrf2: Probing a cDNA library from Mytilus coruscus using a yeast two-hybrid system

The Nuclear factor Erythroid 2-related factor 2 (Nrf2) is the most important endogenous antioxidant factor in organisms, and it has been demonstrated that it exerts extensive control over the immune response by interacting with crucial innate immunity components directly or indirectly. Although Nrf2 has been widely confirmed to be involved in stress resistance in mammals and some fish, its contribution to mollusks oxidative stress resistance has not frequently been documented. In this investigation, total RNA was taken from the digestive gland of M. coruscus, and a cDNA library was constructed and screened using the GATEWAY recombination technology. The Nrf2 cDNA sequence of M. coruscus was cloned into the pGBKT7 vector to prepare the bait plasmid. Using yeast two-hybrid system, after auxotrophic medium screening, sequencing, and bioinformatics analysis, 13 binding proteins that interacted with Nrf2 were finally identified. They were QM-like protein, 40S ribosomal protein S4 (RPS4), ribosomal protein S2 (RPS2), ribosomal protein L12 (RPL12), EF1-alpha mRNA for elongation factor 1 alpha (eEF1-alpha), ferritin, alpha-amylase, trypsin, vdg3, period clock protein, cyclophilin A isoform 1 (CYP A), serine protease CFSP2, histone variant H2A.Z (H2A.Z). For a better understanding the physiological function of Nrf2 in animals and as a potential target for future research on protein roles in Nrf2 interactions, it is crucial to clarify these protein interactions.

Mussel power: Scoping a nature-based solution to microplastic debris

Microplastics are a prolific environmental contaminant. Curbing microplastic pollution requires an array of globally relevant interventions, including source-reduction and curative measures. A novel, nature-based solution to microplastics is proposed, in which mussels are deployed in aquatic ecosystems to act as microplastic biofilters, removing waterborne microplastics and repackaging them into biodeposits that are subsequently captured and removed. Blue mussels (Mytilus edulis) were used to establish the feasibility of such an approach. In the laboratory, mussels were exposed to representative microplastics in a flume tank; at an initial concentration of 1000 microplastics L-1, mussels reduced waterborne microplastic concentrations at an average rate of 40,146 microplastics kg-1 h-1. Mussel faeces sank irrespective of microplastic content, with average sinking velocities of 223-266 m day-1. Modelling predicts ∼3 × 109 mussels deployed on ropes at the mouths of estuaries could remove 4% of waterborne microplastics discharged from rivers. Mussels were successfully deployed in a prototype biodeposit collection system in an urban marina, with 5.0 kg of mussels removing and repackaging 239.9 ± 145.9 microplastics and anthropogenic particles day-1 into their faeces. These results provide impetus for further development of nature-based solutions targeting plastic debris.

Chromosome-Level Genome Assembly of the Blue Mussel Mytilus chilensis Reveals Molecular Signatures Facing the Marine Environment

The blue mussel Mytilus chilensis is an endemic and key socioeconomic species inhabiting the southern coast of Chile. This bivalve species supports a booming aquaculture industry, which entirely relies on artificially collected seeds from natural beds that are translocated to diverse physical-chemical ocean farming conditions. Furthermore, mussel production is threatened by a broad range of microorganisms, pollution, and environmental stressors that eventually impact its survival and growth. Herein, understanding the genomic basis of the local adaption is pivotal to developing sustainable shellfish aquaculture. We present a high-quality reference genome of M. chilensis, which is the first chromosome-level genome for a Mytilidae member in South America. The assembled genome size was 1.93 Gb, with a contig N50 of 134 Mb. Through Hi-C proximity ligation, 11,868 contigs were clustered, ordered, and assembled into 14 chromosomes in congruence with the karyological evidence. The M. chilensis genome comprises 34,530 genes and 4795 non-coding RNAs. A total of 57% of the genome contains repetitive sequences with predominancy of LTR-retrotransposons and unknown elements. Comparative genome analysis of M. chilensis and M. coruscus was conducted, revealing genic rearrangements distributed into the whole genome. Notably, transposable Steamer-like elements associated with horizontal transmissible cancer were explored in reference genomes, suggesting putative relationships at the chromosome level in Bivalvia. Genome expression analysis was also conducted, showing putative genomic differences between two ecologically different mussel populations. The evidence suggests that local genome adaptation and physiological plasticity can be analyzed to develop sustainable mussel production. The genome of M. chilensis provides pivotal molecular knowledge for the Mytilus complex.

The complete mitochondrial genome of a fouling mussel, Xenostrobus atratus (Mollusca: Mytilidae), and its phylogenetic implication

In this study, we report the female-lineage mitochondrial genome of Xenostrobus atratus for the first time. The circular mitochondrial genome is 14,806 bp in length and contains 12 protein-coding genes, 22 transfer RNA genes, and two ribosomal RNA genes. All genes are encoded on the heavy strand. The genome composition is A + T biased (66.6%), with 25.2% A, 41.4% T, 21.7% G and 11.7% C. A Bayesian inference (BI) phylogenetic tree was constructed based on the mitochondrial genomes of X. atratus and 46 other Mytilidae species. Our results demonstrate that X. atratus and Limnoperna fortunei have distinct lineages, opposing synonymizing Xenostrobus within Limnoperna. According to this study, the validity of the subfamily Limnoperninae and genus Xenostrobus is strongly supported. However, there is still an urgent need for more mitochondrial data to decide to which subfamily X. atratus belongs.

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.

Use of limpets as alternative to mussels in metal pollution monitoring; application in the Canary Islands

Metal pollution monitoring programs make use of organisms, such as mussels, as biomonitors to evaluate and compare the metal pollution status of coastal areas worldwide. Despite the widespread distribution of mussels of the genus Mytilus, there are places where these organisms are absent or where their abundance is insufficient for biomonitoring purposes, such as in the Canary Islands (Spain). This study considers the use of limpets of the genus Patella as alternative/complementary species in metal pollution monitoring, and compares the bioaccumulation capacity of mussels and limpets collected simultaneously during several sampling campaigns at 11 sampling sites along the continental Spanish Atlantic coastline. Results show that there are some differences in bioaccumulation capacity between limpets and mussels, with limpets showing in average higher concentrations of Cd, Cr, Ni, As and Cu, and lower concentrations of Zn, while concentrations of Hg and Pb were similar for both types of organisms. Significant correlations between metal concentration in mussels and limpets were found for most metals, with the notable exception of Cd, that showed a very different bioaccumulation pattern in both types of organisms. Limpet to mussel metal concentration ratios were derived that can be used to compare the results of metal biomonitoring using either limpets or mussels, and limpets were successfully used in a monitoring survey in the Canary Islands.

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.

Biomarker-based assessment of sublethal toxicity of organic UV filters (ensulizole and octocrylene) in a sentinel marine bivalve Mytilus edulis

The global occurrence of organic UV filters in the marine environment is of increasing ecotoxicological concern. Here we assessed the toxicity of UV filters ensulizole and octocrylene in the blue mussels Mytilus edulis exposed to 10 or 100 μg l^-1 of octocrylene and ensulizole for two weeks. An integrated battery of biochemical and molecular biomarkers related to xenobiotics metabolism and cellular toxicity (including oxidative stress, DNA damage, apoptosis, autophagy and inflammation) was used to assess the toxicity of these UV filters in the mussels. Octocrylene (but not ensulizole) accumulated in the mussel tissues during the waterborne exposures. Both studied UV filters induced sublethal toxic effects in M. edulis at the investigated concentrations. These effects involved induction of oxidative stress, genotoxicity (indicated by upregulation of DNA damage sensing and repair markers), upregulation of apoptosis and inflammation, and dysregulation of the xenobiotic biotransformation system. Octocrylene induced cellular stress in a concentration-dependent manner, whereas ensulizole appeared to be more toxic at the lower (10 μg l^-1) studied concentration than at 100 μg l^-1. The different concentration-dependence of sublethal effects and distinct toxicological profiles of ensulizole and octocrylene show that the environmental toxicity is not directly related to lipophilicity and bioaccumulation potential of these UV filters and demonstrate the importance of using bioassays for toxicity assessment of emerging pollutants in coastal marine ecosystems.

Where are we? Towards an understanding of the selective accumulation of microplastics in mussels

Mussels are suggested as bioindicators of marine microplastic pollution. However, they are selective in regards to accumulation of microplastics. To make studies more targeted and comparable, ultimately helping to determine the suitability of the mussel as a bioindicator species for microplastic exposure, we review the published literature that has directly or indirectly demonstrated particle selection in mussels. The reported difference between microplastic levels in mussel tissues and environmental matrices provides evidence for their selective uptake characteristics. Both the organ-specific fate characteristics of microplastics, and the different movement patterns of microplastics in the same organ, show that selective translocation processes take place. The selective elimination is reflected in multiple aspects which include (1) the different characteristics of microplastics in excretion and mussel body; (2) the different retention time of various microplastics in mussels; and (3) the tissue-specific change in the numbers of microplastics during the depuration process. This selectivity is affected by the characteristics of the microplastics, the environmental, or laboratory exposure concentrations, feeding status, and other factors. There are still many research gaps and contradictory viewpoints in this field due to this complexity. The current methodology needs improvement and a breakthrough in standardization.

Blue mussels of the Mytilus edulis species complex from South America: The application of species delimitation models to DNA sequence variation

Smooth-shelled blue mussels, Mytilus spp., have a worldwide antitropical distribution and are ecologically and economically important. Mussels of the Mytilus edulis species complex have been the focus of numerous taxonomic and biogeographical studies, in particular in the Northern hemisphere, but the taxonomic classification of mussels from South America remains unclear. The present study analysed 348 mussels from 20 sites in Argentina, Chile, Uruguay and the Falkland Islands on the Atlantic and Pacific coasts of South America. We sequenced two mitochondrial locus, Cytochrome c Oxidase subunit I (625 bp) and 16S rDNA (443 bp), and one nuclear gene, ribosomal 18S rDNA (1770 bp). Mitochondrial and nuclear loci were analysed separately and in combination using maximum likelihood and Bayesian inference methods to identify the combination of the most informative dataset and model. Species delimitation using five different models (GMYC single, bGMYC, PTP, bPTP and BPP) revealed that the Mytilus edulis complex in South America is represented by three species: native M. chilensis, M. edulis, and introduced Northern Hemisphere M. galloprovincialis. However, all models failed to delimit the putative species Mytilus platensis. In contrast, however, broad spatial scale genetic structure in South America using Geneland software to analyse COI sequence variation revealed a group of native mussels (putatively M. platensis) in central Argentina and the Falkland Islands. We discuss the scope of species delimitation methods and the use of nuclear and mitochondrial genetic data to the recognition of species within the Mytilus edulis complex at regional and global scales.

Adaptive Differences in Gene Expression in Farm-Impacted Seedbeds of the Native Blue Mussel Mytilus chilensis

The study of adaptive population differences is relevant for evolutionary biology, as it evidences the power of selective local forces relative to gene flow in maintaining adaptive phenotypes and their underlying genetic determinants. However, human-mediated hybridization through habitat translocations, a common and recurrent aquaculture practice where hybrids could eventually replace local genotypes, risk populations' ability to cope with perturbations. The endemic marine mussel Mytilus chilensis supports a booming farming industry in the inner sea of Chiloé Island, southern Chile, which entirely relies on artificially collected seeds from natural beds that are translocated to ecologically different fattening centers. A matter of concern is how farm-impacted seedbeds will potentially cope with environmental shifts and anthropogenic perturbations. This study provides the first de novo transcriptome of M. chilensis; assembled from tissue samples of mantles and gills of individuals collected in ecologically different farm-impacted seedbeds, Cochamó (41°S) and Yaldad (43°S). Both locations and tissue samples differentially expressed transcripts (DETs) in candidate adaptive genes controlling multiple fitness traits, involved with metabolism, genetic and environmental information processing, and cellular processes. From 189,743 consensus contigs assembled: 1,716 (Bonferroni p value ≤ 0.05) were DETs detected in different tissues of samples from different locations, 210 of them (fold change ≥ | 100|) in the same tissue of samples from a different location, and 665 (fold change ≥ | 4|) regardless of the tissue in samples from a different location. Site-specific DETs in Cochamó (169) and Yaldad (150) in candidate genes controlling tolerance to temperature and salinity shifts, and biomineralization exhibit a high number of nucleotide genetic variants with regular occurrence (frequency > 99%). This novel M. chilensis transcriptome should help assessing and monitoring the impact of translocations in wild and farm-impacted mussel beds in Chiloé Island. At the same time, it would help designing effective managing practices for conservation, and translocation traceability.

Reproductive interference hampers species coexistence despite conspecific sperm precedence

Negative interspecific mating interactions, known as reproductive interference, can hamper species coexistence in a local patch and promote niche partitioning or geographical segregation of closely related species. Conspecific sperm precedence (CSP), which occurs when females that have mated with both conspecific and heterospecific males preferentially use conspecific sperm for fertilization, might contribute to species coexistence by mitigating the costs of interspecific mating and hybridization. We discussed whether two species exhibiting CSP can coexist in a local environment in the presence of reproductive interference. First, using a behaviorally explicit mathematical model, we demonstrated that two species characterized by negative mating interactions are unlikely to coexist because the costs of reproductive interference, such as loss of mating opportunity with conspecific partners, are inevitably incurred when individuals of both species are present. Second, we experimentally examined differences in mating activity and preference in two Harmonia ladybird species known to exhibit CSP. These behavioral differences may lead to local extinction of H. yedoensis because of reproductive interference by H. axyridis. This prediction is consistent with field observations that H. axyridis uses various food sources and habitats whereas H. yedoensis is confined to a less preferred prey item and a pine tree habitat. Finally, by a comparative approach, we observed that niche partitioning or parapatric distribution, but not sympatric coexistence in the same habitat, is maintained between species with CSP belonging to a wide range of taxa, including vertebrates and invertebrates living in aquatic or terrestrial environments. Taken together, it is possible that reproductive interference may destabilize local coexistence even in closely related species that exhibit CSP.

Computational modelling of wet adhesive mussel foot proteins (Bivalvia): Insights into the evolutionary convolution in diverse perspectives

Underwater adhesion in mussels (Bivalvia) is an extreme adaptation to achieve robust and firm wet adhesion in the freshwater/brackish/ocean, which biochemically shaped through millions of years. The protein-based adhesion has huge prospective in various fields like industry, medical, etc. Currently, no comprehensive records related to the systematic documentation of structural and functional properties of Mussel foot proteins (Mfps). In this study, we identified the nine species of bivalves in which the complete sequence of at least one adhesive protein is known. The insilico characterization revealed the specific physio-chemical structural and functional characters of each Mfps. The evolutionary analyses of selected bivalves are mainly based on Mfps, Mitogenome, and TimeTree. The outcome of the works has great applications for designing biomimetic materials in future.

Experimental infection of Mytilus edulis by two Vibrio splendidus-related strains: Determination of pathogenicity level of strains and influence of the origin and annual cycle of mussels on their sensitivity

This study aimed at assessing the pathogenicity of two Vibrio splendidus-related species and evaluating the influence of the origin and annual life cycle of mussels on their sensitivity during a bacterial challenge. Thus, in vivo infection assays were made with Vibrio crassostreae 7T4_12 and Vibrio splendidus 3G1_6, over, respectively, thirteen and 9 months, on adult blue mussels from five recruitment areas in France. Two bacterial concentrations were tested: one consistent with the loads of Vibrio spp. in environment and mussel tissues (~10⁵  CFU/ml) and another one much higher (~10⁸  CFU/ml). The tested environmental concentration has no pathogenic effect whatever the time of year, the strain used and the origin of mussels. However, at the highest concentration, a pathogenic effect was observed only at specific moments, and one of the origins appeared to be more resistant. The physiological state of mussels-depending on the time of year-seemed significant in mussels' sensitivity, as their recruitment origin. This study is the first to test the pathogenicity of V. splendidus-related strains at concentrations close to what is found in the wild, over the annual cycle of mussels, and considering their origin.

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.

Mantle transcriptome sequencing of Mytilus spp. and identification of putative biomineralization genes

In molluscs, the shell secreted by mantle tissue during the biomineralization process is the first barrier against predators and mechanical damage. Changing environmental conditions, such as ocean acidification, influence shell strength and thus protection of the soft body within. Mussels are marine bivalves with important commercial and ecological value worldwide. Despite this importance, the proteins involved in the biomineralization and pigmentation processes in Mytilus spp. remain unclear, as does taxonomy of Mytilus taxa, though there have been many molecular studies. To further understanding in these areas, this study aimed to characterize and compare mantle transcriptomes of four mussel taxa using next generation sequencing. Mussels representing four taxa, were collected from several localities and RNA from mantle tissue was extracted. RNA sequences obtained were assembled, annotated and potential molecular markers, including simple sequence repeats (SSRs) and single nucleotide polymorphisms (SNPs) were identified. Candidate contigs putatively related to biomineralization and pigmentation processes were then selected and several transcripts were chosen for phylogenetic analyses from the Bivalvia class. Transcriptome comparisons between Mytilus taxa, including gene ontology (GO) enrichment analysis and orthologues identification were performed. Of assembled contigs, 46.57%, 37.28% and 17.53% were annotated using NCBI NR, GO and Kyoto Encyclopedia of Genes and Genomes databases, respectively. Potential SSRs (483) and SNPs (1,497) were identified. Results presented a total of 1,292 contigs putatively involved in biomineralization and melanogenesis. Phylogenetic analyses of α-carbonic anhydrase, chitinase and tyrosinase revealed complex evolutionary history and diversity of these genes, which may be a result of duplication events or adaptation to different environments in mussels and other bivalves. Enrichment analyses revealed GO terms associated with pH and thermal response in Mytilus edulis from the North Sea and M. galloprovincialis from the Mediterranean Sea. The phylogenetic analysis within the genus Mytilus revealed M. californianus and M. coruscus to be genetically more distant from the other taxa: M. trossulus, M. edulis, M. chilensis and M. galloprovincialis. This work represents the first mantle transcriptome comparison between Mytilus taxa and provides contigs putatively involved in biomineralization.

Purification and characterization of an antimicrobial peptide mytichitin-chitin binding domain from the hard-shelled mussel, Mytilus coruscus

An antimicrobial peptide with 55 amino acid residues was purified by C₁₈ reversed-phase high-performance liquid chromatography (HPLC) from foot extract of the hard-shelled mussel, Mytilus coruscus. This peptide showed strong antimicrobial activity against Gram-positive and Gram-negative bacteria, as well as fungi. The purified peptide was determined to have a molecular mass of 6202 Da by matrix-assisted laser desorption/ionization time-of-flight mass spectrophotometry (MALDI-TOF/MS). The identified 20-amino acid sequence of the purified peak by Edman degradation shared 100% identity with the N-terminal regions of mytichitin-1, mytichitin-2, mytichitin-3, mytichitin-4, mytichitin-5, and chitinase-like protein-1, and so was named mytichitin-CBD. The cDNA of mytichitin-CBD was cloned and sequenced by rapid amplification of cDNA ends (RACE). The mRNA transcripts were mainly detected in foot tissue, and they were up-regulated and peaked at 4 h after bacterial infection. We constructed and expressed recombinant mytichitin-CBD protein which displayed antimicrobial activity against Gram-negative bacteria Gram-positive bacteria and the fungus as well as anti-parasitic activity against scuticociliates. The results of this study demonstrate that the peptide isolated from M. coruscus is related to the innate immune system of this marine invertebrate and is a possible alternative to antibiotics.

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).

Single nucleotide polymorphisms in native South American Atlantic coast populations of smooth shelled mussels: hybridization with invasive European Mytilus galloprovincialis

BACKGROUND

Throughout the world, harvesting of mussels Mytilus spp. is based on the exploitation of natural populations and aquaculture. Aquaculture activities include transfers of spat and live adult mussels between various geographic locations, which may result in large-scale changes in the world distribution of Mytilus taxa. Mytilus taxa are morphologically similar and difficult to distinguish. In spite of much research on taxonomy, evolution and geographic distribution, the native Mytilus taxa of the Southern Hemisphere are poorly understood. Recently, single nucleotide polymorphisms (SNPs) have been used to clarify the taxonomic status of populations of smooth shelled mussels from the Pacific coast of South America. In this paper, we used a set of SNPs to characterize, for the first time, populations of smooth shelled mussels Mytilus from the Atlantic coast of South America.

RESULTS

Mytilus spp. samples were collected from eastern South America. Six reference samples from the Northern Hemisphere were used: Mytilus edulis from USA and Northern Ireland, Mytilus trossulus from Canada, and Mytilus galloprovincialis from Spain and Italy. Two other reference samples from the Southern Hemisphere were included: M. galloprovincialis from New Zealand and Mytilus chilensis from Chile. Fifty-five SNPs were successfully genotyped, of which 51 were polymorphic. Population genetic analyses using the STRUCTURE program revealed the clustering of eight populations from Argentina (Mytilus platensis) and the clustering of the sample from Ushuaia with M. chilensis from Chile. All individuals in the Puerto Madryn (Argentina) sample were identified as M. platensis × M. galloprovincialis F2 (88.89%) hybrids, except one that was classified as Mediterranean M. galloprovincialis. No F1 hybrids were observed.

CONCLUSIONS

We demonstrate that M. platensis (or Mytilus edulis platensis) and M. chilensis are distinct native taxa in South America, which indicates that the evolutionary histories of Mytilus taxa along the Atlantic and Pacific coasts differ. M. platensis is endangered by hybridization with M. galloprovincialis that was introduced from Europe into the Puerto Madryn area in Argentina, presumably by accidental introduction via ship traffic. We confirm the occurrence of a native M. chilensis population in southern Argentina on the coast of Patagonia.

Multiple reversals of strand asymmetry in molluscs mitochondrial genomes, and consequences for phylogenetic inferences

Strand asymmetry in nucleotide composition is a remarkable feature of animal mitochondrial genomes. The strand-specific bias in the nucleotide composition of the mtDNA has been known to be highly problematic for phylogenetic analyses. Here, the strand asymmetry was compared across 140 mollusc species and analyzed for a mtDNA fragment including twelve protein-coding genes. The analyses show that almost all species in Gastropoda (except Heterobranchia) and all species in Bivalvia present reversals of strand bias. The skew values on individual genes for all codon positions (P₁₂₃), third codon positions (P₃), and fourfold redundant third codon positions (P_4FD) indicated that CG skews are the best indicators of strand asymmetry. The differences in the patterns of strand asymmetry significantly influenced the amino acid composition of the encoded proteins. These biases are most striking for the amino acids Valine, Cysteine, Asparagine and Threonines, which appear to have evolved asymmetrical exchanges in response to shifts in nucleotide composition. Molluscs with strong variability of genome architectures (ARs) are usually characterized by a reversal of the usual strand bias. Phylogenetic analyses show that reversals of asymmetric mutational constraints have consequences on the phylogenetic inferences, as taxa characterized by reverse strand bias (Heterobranchia and Bivalvia) tend to group together due to long-branch attraction (LBA) artifacts. Neutral Transitions Excluded (NTE) model did not overcome the problem of heterogeneous biases present in molluscs mt genomes, suggested it may not be appropriate for molluscs mt genome data. Further refinement phylogenetic models may help us better understand internal relationships among these diverse organisms.

Human mitochondrial ribosomes can switch structural tRNAs - but when and why?

High resolution cryoEM of mammalian mitoribosomes revealed the unexpected presence of mitochondrially encoded tRNA as a structural component of mitochondrial large ribosomal subunit (mt-LSU). Our previously published data identified that only mitochondrial (mt-) tRNA^Phe and mt-tRNA^Val can be incorporated into mammalian mt-LSU and within an organism there is no evidence of tissue specific variation. When mt-tRNA^Val is limiting, human mitoribosomes can integrate mt-tRNA^Phe instead to generate a translationally competent monosome. Here we discuss the possible reasons for and consequences of the observed plasticity of the structural mt-tRNA integration. We also indicate potential direction for further research that could help our understanding of the mechanistic and evolutionary aspects of this unprecedented system.

Blue mussels (Mytilus edulis spp.) as sentinel organisms in coastal pollution monitoring: A review

The blue mussel (Mytilus spp.) is widely used as a bioindicator for monitoring of coastal water pollution (mussel watch programs). Herein we provide a review of this study field with emphasis on: the suitability of Mytilus spp. as environmental sentinels; uptake and bioaccumulation patterns of key pollutant classes; the use of Mytilus spp. in mussel watch programs; recent trends in Norwegian mussel monitoring; environmental quality standards and background concentrations of key contaminants; pollutant effect biomarkers; confounding factors; particulate contaminants (microplastics, engineered nanomaterials); climate change; harmonization of monitoring procedures; and the use of deployed mussels (transplant caging) in pollution monitoring. Lastly, the overall state of the art of blue mussel pollution monitoring is discussed and some important issues for future research and development are highlighted.

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.