A de novo Full-Length mRNA Transcriptome Generated From Hybrid-Corrected PacBio Long-Reads Improves the Transcript Annotation and Identifies Thousands of Novel Splice Variants in Atlantic Salmon

Characterization of the impact of dietary immunostimulant CpG on the expression of mRNA biomarkers involved in the immune responses in Atlantic salmon (Salmo salar)

Infectious diseases have significantly impacted Atlantic salmon aquaculture worldwide. Modulating fish immunity with immunostimulant-containing functional feeds could be an effective strategy in mitigating disease problems. Previously, we characterized the impact of polyriboinosinic polyribocytidylic acid (pIC) and formalin-killed typical Aeromonas salmonicida bacterin on miRNA expression in Atlantic salmon fed a commercial diet with and without immunostimulant CpG. A set of miRNA biomarkers of Atlantic salmon head kidney responding to pIC and/or bacterin immune stimulations was identified (Xue et al., 2019) [1]. Herein, we report a complementary qPCR study that investigated the impact of the pIC, bacterin and dietary CpG on the expression of immune-relevant mRNAs (n = 31) using the same samples as in the previous study (Xue et al., 2019) [1]. Twenty-six of these genes were predicted target transcripts of the pIC- and/or bacterin-responsive miRNAs identified in the earlier study. The current data showed that pIC and/or bacterin stimulations significantly modulated the majority of the qPCR-analyzed genes involved in various immune pathways. Some genes responded to both stimulations (e.g. tnfa, il10rb, ifng, irf9, cxcr3, campb) while others appeared to be stimulation specific [e.g. irf3, irf7a, il1r1, mxa, mapk3 (pIC only); clra (bacterin only)]. A. salmonicida bacterin stimulation produced a strong inflammatory response (e.g. higher expression of il1b, il8a and tnfa), while salmon stimulated with pIC showed robust interferon responses (both type I and II). Furthermore, the current data indicated significant down-regulation of immune-relevant transcripts (e.g. tlr9, irf5, il1r1, hsp90ab1, itgb2) by dietary immunostimulant CpG, especially among pre-injection and PBS-injected fish. Together with our prior miRNA study, the present research provided complementary information on Atlantic salmon anti-viral and anti-bacterial immune responses and on how dietary CpG may modulate these responses.

Effect of prolonged feeding of broodstock diet with increased inclusion of essential n-3 fatty acids on maturing and spawning performance in 3-year-old Atlantic salmon (Salmo salar)

Atlantic salmon (Salmo salar) broodstock recruits are normally fed a specialized diet with a higher content of essential nutrients for a limited time period prior to fasting and transfer to freshwater. Typically, this period lasts for about six months, but may vary among producers. Reduced use of marine ingredients in commercial salmon diets during the last decades has affected the content of essential nutrients, such as n-3 long chained polyunsaturated fatty acids (LC-PUFA), minerals and vitamins. Furthermore, to minimize the risk of losses and implement new breeding achievements faster, breeding companies have shortened the production cycle of broodstock from 4 to 3 years, which may affect the number of fish that are large enough to mature. In the present study, we have extended the broodstock feeding period from 6 to 15 months prior to the freshwater transfer giving a higher content of n-3 LC-PUFA (higher inclusion of marine oils) from February to December (Phase 1), and thereafter a diet with a higher energy content to ensure growth towards the spring and maturation (Phase 2). Four sea cages with approximately 80.000 salmon postsmolt, two sea cages with males and two with females, were given a control diet and an experimental diet. Samples were taken in Phase 1 at start (1.7 kg), mid (3.4 kg) and end Phase 1/start of Phase 2 (8.3 kg), and end of Phase 2 (13.4 kg). The fish were thereafter fasted, and selected fish transferred to landbased freshwater tanks where light and temperature were used to manipulate the spawning time of the fish in two groups (early or late). Due to disease in the facility, measures of egg quality and hatching were only obtained from the early group. During the trial and spawning period, biometrical measurements were recorded, and samples of liver, gonad, fillet and red blood cells (RBC) were collected for fatty acid composition and blood plasma for analysis of lipid and health-related parameters. Samples were also collected for gonadal transcriptomic analysis by microarray and qPCR (end Phase 2) and plasma steroids (end Phase 2, mid maturation and spawning). Males fed the test diet had a larger body size compared to the control group at the end of Phase 2, while no differences were observed between dietary groups for the females. Total mortality in the trial was lower in the test group compared to the control, losses were caused mainly by sea lice treatments, loser fish or cardiomyopathy syndrome (CMS). The dietary LC-PUFA levels in the test diet were reflected in the tissues particularly during Phase 1, but only different in the fillet samples and eggs at the end of Phase 2 and at spawning. Plasma sex steroids content increased at mid maturation and showed lower levels of androgens and estrogens in females fed the test diet compared to the control. At the end of Phase 2, transcriptional analysis showed upregulation of steroidogenic enzymes, although not reflected in changes in plasma steroids in Phase 2, indicating changes to come during maturation. The differences in LC-PUFA content in tissues and plasma steroids did not appear to affect fecundity, sperm quality, egg survival or hatching rate, but the test group had larger eggs compared to the control in the early spawner-group. Prolonged feeding of n-3 LC-PUFA to pre-puberty Atlantic salmon broodstock appears to be important for higher survival in challenging sea cage environments and has an effect on sex steroid production that, together with high energy diet during early maturation, cause the test group to produce larger eggs.

Catch of the Day: New Serum Amyloid A (SAA) Antibody Is a Valuable Tool to Study Fish Health in Salmonids

Serum amyloid A (SAA) proteins belong to a family of acute-phase reactants, playing an integral role in defending the organism from pathological damage. Despite a wealth of data on the regulation of SAA transcripts in teleosts, there is only limited information on these proteins' abundance in fish. The aim of this study is to characterise SAA protein levels in salmonids using a newly developed antibody specific to salmonid SAA. The salmonid SAA antibody detected SAA and accurately discriminated between stimulated and control specimens from rainbow trout macrophage cell line (RTS-11) in vitro, as well as rainbow trout challenged with Aeromonas salmonicida- or flagellin-stimulated Atlantic salmon in vivo. The presence of SAA protein was analysed in RTS-11 cell line supernatants, liver, and spleen samples using ELISA, immunoblotting, and immunohistochemistry. This study is the first to characterise SAA protein levels in salmonids in vivo and in vitro. The newly developed salmonid SAA antibody was able to discriminate between stimulated and unstimulated specimens, showing that it can be used to study the acute-phase response in salmonids with the potential to be further developed into assays to monitor and evaluate health in wild and farmed fish.

Full-length transcriptome from different life stages of cobia (Rachycentron canadum, Rachycentridae)

Cobia (Rachycentron canadum, Rachycentridae) is one of the prospective species for mariculture. The transcriptome-based study on cobia was hampered by an inadequate reference genome and a lack of full-length cDNAs. We used a long-read based sequencing technology (PacBio Sequel II Iso-Seq3 SMRT) to obtain complete transcriptome sequences from larvae, juveniles, and various tissues of adult cobia, and a single SMRTcell generated 99 gigabytes of data and 51,205,946,694 bases. A total of 8609435, 7441673 and 9140164 subreads were generated from the larval, juvenile, and adult sample pools, with mean sub-read lengths of 2109.9, 1988.2 and 1996.2 bp, respectively. All samples were combined to increase transcript recovery and clustered into 35661 high-quality reads. This is the first report on a full-length transcriptome from R. canadum. Our results illustrate a significant increase in the identified amount of cobia LncRNAs and alternatively spliced transcripts, which will help improve genome annotation. Furthermore, this information will be beneficial for nutrigenomics and functional studies on cobia and other commercially important mariculture species.

CODA: a combo-Seq data analysis workflow

The analysis of the combined mRNA and miRNA content of a biological sample can be of interest for answering several research questions, like biomarkers discovery, or mRNA-miRNA interactions. However, the process is costly and time-consuming, separate libraries need to be prepared and sequenced on different flowcells. Combo-Seq is a library prep kit that allows us to prepare combined mRNA-miRNA libraries starting from very low total RNA. To date, no dedicated bioinformatics method exists for the processing of Combo-Seq data. In this paper, we describe CODA (Combo-seq Data Analysis), a workflow specifically developed for the processing of Combo-Seq data that employs existing free-to-use tools. We compare CODA with exceRpt, the pipeline suggested by the kit manufacturer for this purpose. We also evaluate how Combo-Seq libraries analysed with CODA perform compared with conventional poly(A) and small RNA libraries prepared from the same samples. We show that using CODA more successfully trimmed reads are recovered compared with exceRpt, and the difference is more dramatic with short sequencing reads. We demonstrate how Combo-Seq identifies as many genes and fewer miRNAs compared to the standard libraries, and how miRNA validation favours conventional small RNA libraries over Combo-Seq. The CODA code is available at https://github.com/marta-nazzari/CODA.

Full-Length Transcriptome of Myotis pilosus as a Reference Resource and Mining of Auditory and Immune Related Genes

Rickett's big-footed bat, Myotis pilosus, which belongs to the family Vespertilionida, is the only known piscivorous bat in East Asia. Accurate whole genome and transcriptome annotations are essential for the study of bat biological evolution. The lack of a whole genome for M. pilosus has limited our understanding of the molecular mechanisms underlying the species' evolution, echolocation, and immune response. In the present work, we sequenced the entire transcriptome using error-corrected PacBio single-molecule real-time (SMRT) data. Then, a total of 40 GB of subreads were generated, including 29,991 full-length non-chimeric (FLNC) sequences. After correction by Illumina short reads and de-redundancy, we obtained 26,717 error-corrected isoforms with an average length of 3018.91 bp and an N50 length of 3447 bp. A total of 1528 alternative splicing (AS) events were detected by transcriptome structural analysis. Furthermore, 1032 putative transcription factors (TFs) were identified, with additional identification of several long non-coding RNAs (lncRNAs) with high confidence. Moreover, several key genes, including PRL-2, DPP4, Glul, and ND1 were also identified as being associated with metabolism, immunity, nervous system processes, and auditory perception. A multitude of pattern recognition receptors was identified, including NLR, RLR, SRCR, the antiviral molecule IRF3, and the IFN receptor subunit IFNAR1. High-quality reference genomes at the transcriptome level may be used to quantify gene or transcript expression, evaluate alternative splicing levels, identify novel transcripts, and enhance genome annotation in bats.

Expression Analysis of Moritella viscosa-Challenged Atlantic Salmon Identifies Disease-Responding Genes, MicroRNAs and Their Predicted Target Genes and Pathways

Moritella viscosa is a bacterial pathogen causing winter-ulcer disease in Atlantic salmon. The lesions on affected fish lead to increased mortality, decreased fish welfare, and inferior meat quality in farmed salmon. MicroRNAs (miRNAs) are small non-coding RNAs involved in post-transcriptional regulation by guiding the miRNA-induced silencing complex to specific mRNA transcripts (target genes). The goal of this study was to identify miRNAs responding to Moritella viscosa in salmon by investigating miRNA expression in the head-kidney and the muscle/skin from lesion sites caused by the pathogen. Protein coding gene expression was investigated by microarray analysis in the same materials. Seventeen differentially expressed guide-miRNAs (gDE-miRNAs) were identified in the head-kidney, and thirty-nine in lesion sites, while the microarray analysis reproduced the differential expression signature of several thousand genes known as infection-responsive. In silico target prediction and enrichment analysis suggested that the gDE-miRNAs were predicted to target genes involved in immune responses, hemostasis, angiogenesis, stress responses, metabolism, cell growth, and apoptosis. The majority of the conserved gDE-miRNAs (e.g., miR-125, miR-132, miR-146, miR-152, miR-155, miR-223 and miR-2188) are known as infection-responsive in other vertebrates. Collectively, the findings indicate that gDE-miRNAs are important post-transcriptional gene regulators of the host response to bacterial infection.

Differential Expression of miRNAs and Their Predicted Target Genes Indicates That Gene Expression in Atlantic Salmon Gill Is Post-Transcriptionally Regulated by miRNAs in the Parr-Smolt Transformation and Adaptation to Sea Water

Smoltification (parr-smolt transformation) is a complex developmental process consisting of developmental changes that lead to remodeling of the Atlantic salmon gill. Here, the expression changes of miRNAs and mRNAs were studied by small-RNA sequencing and microarray analysis, respectively, to identify miRNAs and their predicted targets associated with smoltification and subsequent sea water adaptation (SWA). In total, 18 guide miRNAs were identified as differentially expressed (gDE miRNAs). Hierarchical clustering analysis of expression changes divided these into one cluster of 13 gDE miRNAs with decreasing expression during smoltification and SWA that included the miRNA-146, miRNA-30 and miRNA-7132 families. Another smaller cluster that showed increasing expression consisted of miR-101a-3p, miR-193b-5p, miR-499a-5p, miR-727a-3p and miR-8159-5p. The gDE miRNAs were predicted to target 747 of the genes (DE mRNAs), showing expression changes in the microarray analysis. The predicted targets included genes encoding NKA-subunits, aquaporin-subunits, cystic fibrosis transmembrane conductance regulator and the solute carrier family. Furthermore, the predicted target genes were enriched in biological processes associated with smoltification and SWA (e.g., immune system, reactive oxygen species, stress response and extracellular matrix organization). Collectively, the results indicate that remodeling of the gill involves the post-transcriptional regulation of gene expression by the characterized gDE miRNAs.

Modern venomics-Current insights, novel methods, and future perspectives in biological and applied animal venom research

Venoms have evolved >100 times in all major animal groups, and their components, known as toxins, have been fine-tuned over millions of years into highly effective biochemical weapons. There are many outstanding questions on the evolution of toxin arsenals, such as how venom genes originate, how venom contributes to the fitness of venomous species, and which modifications at the genomic, transcriptomic, and protein level drive their evolution. These questions have received particularly little attention outside of snakes, cone snails, spiders, and scorpions. Venom compounds have further become a source of inspiration for translational research using their diverse bioactivities for various applications. We highlight here recent advances and new strategies in modern venomics and discuss how recent technological innovations and multi-omic methods dramatically improve research on venomous animals. The study of genomes and their modifications through CRISPR and knockdown technologies will increase our understanding of how toxins evolve and which functions they have in the different ontogenetic stages during the development of venomous animals. Mass spectrometry imaging combined with spatial transcriptomics, in situ hybridization techniques, and modern computer tomography gives us further insights into the spatial distribution of toxins in the venom system and the function of the venom apparatus. All these evolutionary and biological insights contribute to more efficiently identify venom compounds, which can then be synthesized or produced in adapted expression systems to test their bioactivity. Finally, we critically discuss recent agrochemical, pharmaceutical, therapeutic, and diagnostic (so-called translational) aspects of venoms from which humans benefit.

Expression Analysis in Atlantic Salmon Liver Reveals miRNAs Associated with Smoltification and Seawater Adaptation

Simple Summary

Smoltification is a developmental process that preadapts Atlantic salmon for a life in seawater. Suboptimal smoltification and poor timing of transfer to seawater is associated with increased mortality. MicroRNAs (miRNAs) are small non-coding genes. They regulate gene expression post-transcriptionally as part of the miRNA induce silencing complex (miRISC) where they guide miRISC to particular mRNAs (target genes). The aim of this study was to identify Atlantic salmon miRNAs expressed in liver that are associated with smoltification and adaptation to seawater as well as to predict their target genes. In total, 62 guide miRNAs were identified, and by their expression patterns they were clustered into three groups. Target gene predictions followed by gene enrichment analysis of the predicted targets indicated that the guide miRNAs were involved in post-transcriptional regulation of important smoltification associated biological processes. Some of these were energy metabolism, protein metabolism and transport, circadian rhythm, stress and immune response. Together, the results indicate that certain miRNAs are involved in the regulation of many of the important changes occurring in the liver during this developmental transition.

Abstract

Optimal smoltification is crucial for normal development, growth, and health of farmed Atlantic salmon in seawater. Here, we characterize miRNA expression in liver to reveal whether miRNAs regulate gene expression during this developmental transition. Expression changes of miRNAs and mRNAs was studied by small-RNA sequencing and microarray analysis, respectively. This revealed 62 differentially expressed guide miRNAs (gDE-miRNAs) that could be divided into three groups with characteristic dynamic expression patterns. Three of miRNA families are known as highly expressed in liver. A rare arm shift was observed during smoltification in the Atlantic salmon-specific novel-ssa-miR-16. The gDE-miRNAs were predicted to target 2804 of the genes revealing expression changes in the microarray analysis. Enrichment analysis revealed that targets were significantly enriched in smoltification-associated biological process groups. These included lipid and cholesterol synthesis, carbohydrate metabolism, protein metabolism and protein transport, immune system genes, circadian rhythm and stress response. The results indicate that gDE-miRNAs may regulate many of the changes associated with this developmental transition in liver. The results pave the way for validation of the predicted target genes and further study of gDE-miRNA and their targets by functional assays.

Long-read isoform sequencing reveals tissue-specific isoform expression between active and hibernating brown bears (Ursus arctos)

Understanding hibernation in brown bears (Ursus arctos) can provide insight into some human diseases. During hibernation, brown bears experience periods of insulin resistance, physical inactivity, extreme bradycardia, obesity, and the absence of urine production. These states closely mimic aspects of human diseases such as type 2 diabetes, muscle atrophy, as well as renal and heart failure. The reversibility of these states from hibernation to active season enables the identification of mediators with possible therapeutic value for humans. Recent studies have identified genes and pathways that are differentially expressed between active and hibernation seasons in bears. However, little is known about the role of differential expression of gene isoforms on hibernation physiology. To identify both distinct and novel mRNA isoforms, full-length RNA-sequencing (Iso-Seq) was performed on adipose, skeletal muscle, and liver from three individual bears sampled during both active and hibernation seasons. The existing reference genome annotation was improved by combining it with the Iso-Seq data. Short-read RNA-sequencing data from six individuals were mapped to the new reference annotation to quantify differential isoform usage (DIU) between tissues and seasons. We identified differentially expressed isoforms in all three tissues, to varying degrees. Adipose had a high level of DIU with isoform switching, regardless of whether the genes were differentially expressed. Our analyses revealed that DIU, even in the absence of differential gene expression, is an important mechanism for modulating genes during hibernation. These findings demonstrate the value of isoform expression studies and will serve as the basis for deeper exploration into hibernation biology.

MicroSalmon: A Comprehensive, Searchable Resource of Predicted MicroRNA Targets and 3'UTR Cis-Regulatory Elements in the Full-Length Sequenced Atlantic Salmon Transcriptome

Complete 3′UTRs unambiguously assigned to specific mRNA isoforms from the Atlantic salmon full-length (FL) transcriptome were collected into a 3′UTRome. miRNA response elements (MREs) and other cis-regulatory motifs were subsequently predicted and assigned to 3′UTRs of all FL-transcripts. The MicroSalmon GitHub repository provides all results. RNAHybrid and sRNAtoolbox tools predicted the MREs. UTRscan and the Teiresias algorithm predicted other 3′UTR cis-acting motifs, both known vertebrate motifs and putative novel motifs. MicroSalmon provides search programs to retrieve all FL-transcripts targeted by a miRNA (median number 1487), all miRNAs targeting an FL-transcript (median number 27), and other cis-acting motifs. As thousands of FL-transcripts may be targets of each miRNA, additional experimental strategies are necessary to reduce the likely true and relevant targets to a number that may be functionally validated. Low-complexity motifs known to affect mRNA decay in vertebrates were over-represented. Many of these were enriched in the terminal end, while purine- or pyrimidine-rich motifs with unknown functions were enriched immediately downstream of the stop codon. Furthermore, several novel complex motifs were over-represented, indicating conservation and putative function. In conclusion, MicroSalmon is an extensive and useful, searchable resource for study of Atlantic salmon transcript regulation by miRNAs and cis-acting 3′UTR motifs.

PacBio Iso-Seq Improves the Rainbow Trout Genome Annotation and Identifies Alternative Splicing Associated With Economically Important Phenotypes

Rainbow trout is an important model organism that has received concerted international efforts to study the transcriptome. For this purpose, short-read sequencing has been primarily used over the past decade. However, these sequences are too short of resolving the transcriptome complexity. This study reported a first full-length transcriptome assembly of the rainbow trout using single-molecule long-read isoform sequencing (Iso-Seq). Extensive computational approaches were used to refine and validate the reconstructed transcriptome. The study identified 10,640 high-confidence transcripts not previously annotated, in addition to 1,479 isoforms not mapped to the current Swanson reference genome. Most of the identified lncRNAs were non-coding variants of coding transcripts. The majority of genes had multiple transcript isoforms (average ∼3 isoforms/locus). Intron retention (IR) and exon skipping (ES) accounted for 56% of alternative splicing (AS) events. Iso-Seq improved the reference genome annotation, which allowed identification of characteristic AS associated with fish growth, muscle accretion, disease resistance, stress response, and fish migration. For instance, an ES in GVIN1 gene existed in fish susceptible to bacterial cold-water disease (BCWD). Besides, under five stress conditions, there was a commonly regulated exon in prolyl 4-hydroxylase subunit alpha-2 (P4HA2) gene. The reconstructed gene models and their posttranscriptional processing in rainbow trout provide invaluable resources that could be further used for future genetics and genomics studies. Additionally, the study identified characteristic transcription events associated with economically important phenotypes, which could be applied in selective breeding.