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 growing amount of genome information and transcriptomes data available allows for a better understanding of biological processes. However, analysis of complex transcriptomic experimental designs involving different conditions, tissues, or times is relevant. This study proposes a novel approach to analyze complex data sets combining transcriptomes and miRNAs at the chromosome-level genome. Atlantic salmon smolts were transferred to seawater under two strategies: (i) fish group exposed to gradual salinity changes (GSC) and (ii) fish group exposed to a salinity shock (SS). Gills, intestine, and head kidney samples were used for total RNA extraction, followed by mRNA and small RNA illumina sequencing. Different expression patterns among the tissues and treatments were observed through a whole-genome transcriptomic approach. Chromosome regions highly expressed between experimental conditions included a great abundance of transposable elements. In addition, differential expression analysis showed a greater number of transcripts modulated in response to SS in gills and head kidney. miRNA expression analysis suggested a small number of miRNAs involved in the smoltification process. However, target analysis of these miRNAs showed a regulatory role in growth, stress response, and immunity. This study is the first to evidence the interplaying among mRNAs and miRNAs and the structural relationship at the genome level during Atlantic salmon smoltification.

Salmon aquaculture is constantly threatened by pathogens that impact fish health, welfare, and productivity, including the sea louse Caligus rogercresseyi. This marine ectoparasite is mainly controlled through delousing drug treatments that have lost efficacy. Therein, strategies such as salmon breeding selection represent a sustainable alternative to produce fish with resistance to sea lice. This study explored the whole-transcriptome changes in Atlantic salmon families with contrasting resistance phenotypes against lice infestation. In total, 121 Atlantic salmon families were challenged with 35 copepodites per fish and ranked after 14 infestation days. Skin and head kidney tissue from the top two lowest (R) and highest (S) infested families were sequenced by the Illumina platform. Genome-scale transcriptome analysis showed different expression profiles between the phenotypes. Significant differences in chromosome modulation between the R and S families were observed in skin tissue. Notably, the upregulation of genes associated with tissue repairs, such as collagen and myosin, was found in R families. Furthermore, skin tissue of resistant families showed the highest number of genes associated with molecular functions such as ion binding, transferase, and cytokine activity, compared with the susceptible. Interestingly, lncRNAs differentially modulated in the R/S families are located near genes associated with immune response, which are upregulated in the R family. Finally, SNPs variations were identified in both salmon families, where the resistant ones showed the highest number of SNPs variations. Remarkably, among the genes with SPNs, genes associated with the tissue repair process were identified. This study reported Atlantic salmon chromosome regions exclusively expressed in R or S Atlantic salmon families’ phenotypes. Furthermore, due to the presence of SNPs and high expression of tissue repair genes in the resistant families, it is possible to suggest mucosal immune activation associated with the Atlantic salmon resistance to sea louse infestation.

The increase in hypoxia events, a result of climate change in coastal and fjord ecosystems, impacts the health and survival of mussels. These organisms deploy physiological and molecular responses as an adaptive mechanism to maintain cellular homeostasis under environmental stress. However, the specific effects of hypoxia on mussels of socioeconomic interest, such as Mytilus chilensis, are unknown. Using RNA-seq, we investigated the transcriptomic profiles of the gills, digestive gland, and adductor muscle of M. chilensis under hypoxia (10 days at 2 mg L⁻¹) and reoxygenation (10 days at 6 mg L⁻¹). There were 15,056 differentially expressed transcripts identified in gills, 11,864 in the digestive gland, and 9862 in the adductor muscle. The response varied among tissues, showing chromosomal changes in Chr1, Chr9, and Chr10 during hypoxia. Hypoxia regulated signaling genes in the Toll-like, mTOR, citrate cycle, and apoptosis pathways in gills, indicating metabolic and immunological alterations. These changes suggest that hypoxia induced a metabolic shift in mussels, reducing reliance on aerobic respiration and increasing reliance on anaerobic metabolism. Furthermore, hypoxia appeared to suppress the immune response, potentially increasing disease susceptibility, with negative implications for the mussel culture industry and natural bed populations. This study provides pivotal insights into metabolic and immunological adaptations to hypoxia in M. chilensis, offering candidate genes for adaptive traits.

The host’s physiological history and environment determine the microbiome structure. In that sense, the strategy used for the salmon transfer to seawater after parr-smolt transformation may influence the Atlantic salmon’s intestinal microbiota. Therefore, this study aimed to explore the diversity and abundance of the Atlantic salmon intestinal microbiota and metagenome functional prediction during seawater transfer under three treatments. One group was exposed to gradual salinity change (GSC), the other to salinity shock (SS), and the third was fed with a functional diet (FD) before the seawater (SW) transfer. The microbial profile was assessed through full-16S rRNA gene sequencing using the Nanopore platform. In addition, metagenome functional prediction was performed using PICRUSt2. The results showed an influence of salinity changes on Atlantic salmon gut microbiota richness, diversity, and taxonomic composition. The findings reveal that GSC and the FD increased the Atlantic salmon smolt microbiota diversity, suggesting a positive association between the intestinal microbial community and fish health during seawater transfer. The reported knowledge can be applied to surveil the microbiome in smolt fish production, improving the performance of Atlantic salmon to seawater transfer.

Caligus rogercresseyi, commonly known as sea louse, is an ectoparasite copepod that impacts the salmon aquaculture in Chile, causing losses of hundreds of million dollars per year. This pathogen is mainly controlled by immersion baths with delousing drugs, which can lead to resistant traits selection in lice populations. Bioassays are commonly used to assess louse drug sensitivity, but the current procedures may mask relevant molecular responses. This study aimed to discover novel coding genes and non-coding RNAs that could evidence drug sensitivity at the genomic level. Sea lice samples from populations with contrasting sensitivity to delousing drugs were collected. Bioassays using azamethiphos, cypermethrin, and deltamethrin drugs were conducted to evaluate the sensitivity and to collect samples for RNA-sequencing. Transcriptome sequencing was conducted on samples exposed to each drug to evaluate the presence of coding and non-coding RNAs associated with the response of these compounds. The results revealed specific transcriptome patterns in lice exposed to azamethiphos, deltamethrin, and cypermethrin drugs. Enrichment analyses of Gene Ontology terms showed specific biological processes and molecular functions associated with each delousing drug analyzed. Furthermore, novel long non-coding RNAs (lncRNAs) were identified in C. rogercresseyi and tightly linked to differentially expressed coding genes. A significant correlation between gene transcription patterns and phenotypic effects was found in lice collected from different salmon farms with contrasting drug treatment efficacies. The significant correlation among gene transcription patterns with the historical background of drug sensitivity suggests novel molecular mechanisms of pharmacological resistance in lice populations.