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.

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.

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.

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.

Vibrio sp. ArtGut-C1, a polyhydroxybutyrate producer isolated from the gut of the aquaculture live diet Artemia (Crustacea)

BACKGROUND: Vibrio species display variable and plastic fitness strategies to survive and interact with multiple hosts, including marine aquaculture species that are severely affected by pathogenic Vibrios. The culturable Vibrio sp. strain ArtGut-C1, the focus of this study, provides new evidence of such phenotypic plasticity as it accumulates polyhydroxybutyrate (PHB), a biodegradable polymer with anti-pathogen activity, particularly in the marine larviculture phase. The strain was isolated from the gut of laboratory-reared Artemia individuals, the live diet and PHB carrier used in larviculture. Its main phenotypic properties, taxonomic status and genomic properties are reported based on the whole-genome sequencing. RESULTS: Vibrio sp. ArtGut-C1 yielded 72.6% PHB of cells' dry weight at 25 C. The genomic average nucleotide identity (ANI) shows it is closely related to V. diabolicus (ANI: 88.6%). Its genome contains 5,236,997- bp with 44.8% GC content, 3,710 protein-coding sequences, 96 RNA, 9 PHB genes functionally related to PHB metabolic pathways, and several genes linked to competing and colonizing abilities. CONCLUSIONS: This culturable PHB-accumulating Vibrio strain shows high genomic and phenotypic variability. It may be used as a natural pathogen biocontrol in the marine hatchery and as a potential cell factory for PHB production.

Distribution and growth of Vibrio parahaemolyticus in southern Chilean clams (Venus antiqua) and blue mussels (Mytilus chilensis).

We evaluated the distribution and growth of Vibrio parahaemolyticus in the inland sea of southern Chile, where the world's largest foodborne gastroenteritis outbreak by the pandemic strain O3:K6 occurred in 2005. Intertidal samples of Mytilus chilensis and Venus antiqua were collected around port towns between 41°28'S and 43°07'S, during April to May 2011 and January to March 2012. We used most probable number real-time polymerase chain reaction (MPN-PCR) for enumeration of the tlh, tdh, and trh genes in freshly harvested bivalves and after a controlled postharvest temperature abuse. Pathogenic markers (tdh+ or trh+) were not detected. Total V. parahaemolyticus (tlh+) in freshly harvested samples reached up to 0.38 and 3.66 log MPN/g in 2011 and 2012, respectively, with values close to or above 3 log MPN/g only near Puerto Montt (41°28'S, 72°55'W). Enrichments by temperature abuse (>2 log MPN/g) occurred mainly in the same zone, regardless of the year, suggesting that both natural or anthropogenic exposure to high temperatures were more critical. Lower salinity and higher sea surface temperature in Reloncaví Sound and Reloncaví Estuary were consistent with our observations and allowed confirmation of the existence of a high-risk zone near Puerto Montt. Based on the results, a strategy focused on risk management inside this defined hazard zone is recommended.