Exploring the transcriptome of Atlantic salmon (Salmo salar) skin, a major defense organ

Transcriptome Profiling of Gene Expression in Atlantic Salmon (Salmo salar) at Early Stage of Development

For Atlantic salmon development, the most critical phase is the early development stage from egg to fry through alevin. However, the studies investigating the early development of Atlantic salmon based on RNA-seq are scarce and focus only on one stage of development. Therefore, using the RNA-seq technology, the assessment of different gene expressions of various early development stages (egg, alevin, and fry) was performed on a global scale. Over 22 GB of clean data was generated from 9 libraries with three replicates for each stage with over 90% mapping efficiency. A total of 5534 genes were differentially expressed, among which 19, 606, and 826 genes were specifically expressed in each stage, respectively. The transcriptome analysis showed that the number of differentially expressed genes (DEGs) increased as the Atlantic salmon progressed in development from egg to fry stage. In addition, gene ontology enrichment demonstrated that egg and alevin stages are characterized by upregulation of genes involved in spinal cord development, neuron projection morphogenesis, axonogenesis, and cytoplasmic translation. At the fry stage, upregulated genes were enriched in the muscle development process (muscle cell development, striated muscle cell differentiation, and muscle tissue development), immune system (defense response and canonical NF-kappaB signal transduction), as well as epidermis development. These results suggest that the early development of Atlantic salmon is characterized by a dynamic shift in gene expression and DEGs between different stages, which provided a solid foundation for the investigation of Atlantic salmon development.

Transcriptomic landscape of Atlantic salmon (Salmo salar L.) skin

In this study, we present the first spatial transcriptomic atlas of Atlantic salmon skin using the Visium Spatial Gene Expression protocol. We utilized frozen skin tissue from 4 distinct sites, namely the operculum, pectoral and caudal fins, and scaly skin at the flank of the fish close to the lateral line, obtained from 2 Atlantic salmon (150 g). High-quality frozen tissue sections were obtained by embedding tissue in optimal cutting temperature media prior to freezing and sectioning. Further, we generated libraries and spatial transcriptomic maps, achieving a minimum of 80 million reads per sample with mapping efficiencies ranging from 79.3 to 89.4%. Our analysis revealed the detection of over 80,000 transcripts and nearly 30,000 genes in each sample. Among the tissue types observed in the skin, the epithelial tissues exhibited the highest number of transcripts (unique molecular identifier counts), followed by muscle tissue, loose and fibrous connective tissue, and bone. Notably, the widest nodes in the transcriptome network were shared among the epithelial clusters, while dermal tissues showed less consistency, which is likely attributable to the presence of multiple cell types at different body locations. Additionally, we identified collagen type 1 as the most prominent gene family in the skin, while keratins were found to be abundant in the epithelial tissue. Furthermore, we successfully identified gene markers specific to epithelial tissue, bone, and mesenchyme. To validate their expression patterns, we conducted a meta-analysis of the microarray database, which confirmed high expression levels of these markers in mucosal organs, skin, gills, and the olfactory rosette.

Heritability of Immunity Traits and Resistance of Atlantic Salmon against the Sea Louse Caligus rogercresseyi

The immune response of Atlantic salmon to sea lice has been extensively studied, but we still do not know the mechanisms by which some fish become resistant and others do not. In this study, we estimated the heritabilities of three key proteins associated with the innate immunity and resistance of Salmo salar against the sea louse Caligus rogercresseyi. In particular, we quantified the abundance of 2 pro-inflammatory cytokines, Tnfα and Il-8, and an antioxidant enzyme, Nkef, in Atlantic salmon skin and gill tissue from 21 families and 268 individuals by indirect ELISA. This covers a wide parasite load range from low or resistant (mean sea lice ± SE = 8.7 ± 0.9) to high or susceptible (mean sea lice ± SE = 43.3 ± 2.0). Our results showed that susceptible fish had higher levels of Nkef and Tnfα than resistant fish in their gills and skin, although gill Il-8 was higher in resistant fish, while no significant differences were found in the skin. Furthermore, moderate to very high heritable genetic variation was estimated for Nkef (h2 skin: 0.96 ± 0.14 and gills: 0.97 ± 0.11) and Tnfα (h2 skin: 0.53 ± 0.17 and gills: 0.32 ± 0.14), but not for Il-8 (h2 skin: 0.22 ± 0.12 ns and gills: 0.09 ± 0.08 ns). This work provides evidence that Nkef and Tnfα protein expressions are highly heritable and related to resistance against sea lice in Atlantic salmon.

Integrated Analysis of lncRNA and circRNA Mediated ceRNA Regulatory Networks in Skin Reveals Innate Immunity Differences Between Wild-Type and Yellow Mutant Rainbow Trout (Oncorhynchus mykiss)

Fish skin is a vital immune organ that forms the first protective barrier preventing entry of external pathogens. Rainbow trout is an important aquaculture fish species that is farmed worldwide. However, our knowledge of innate immunity differences between wild-type (WR_S) and yellow mutant rainbow trout (YR_S) remains limited. In this study, we performed whole transcriptome analysis of skin from WR_S and YR_S cultured in a natural flowing water pond. A total of 2448 mRNAs, 1630 lncRNAs, 22 circRNAs and 50 miRNAs were found to be differentially expressed (DE). Among these DEmRNAs, numerous key immune-related genes, including ifih1, dhx58, trim25, atp6v1e1, tap1, tap2, cd209, hsp90a.1, nlrp3, nlrc3, and several other genes associated with metabolism (gstp1, nampt, naprt and cd38) were identified. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of DEmRNAs revealed that many were significantly enriched in innate immune-related GO terms and pathways, including NAD+ADP-ribosyltransferase activity, complement binding, immune response and response to bacterium GO terms, and RIG-I-like receptor signaling, NOD-like receptor signaling and phagosome KEGG pathways. Furthermore, the immune-related competing endogenous RNA networks were constructed, from which we found that lncRNAs MSTRG.11484.2, MSTRG.32014.1 and MSTRG.29012.1 regulated at least three immune-related genes (ifih1, dhx58 and irf3) through PC-5p-43254_34, PC-3p-28352_70 and bta-miR-11987_L-1R-1_1ss8TA, and tap2 was regulated by two circRNAs (circRNA5279 and circRNA5277) by oni-mir-124a-2-p5_1ss13GA. The findings expand our understanding of the innate immune system of rainbow trout, and lay the foundation for further study of immune mechanisms and disease resistance breeding.

Transcriptomic characterisation of a common skin lesion in farmed chinook salmon

Little is known about host responses of farmed Chinook salmon with skin lesions, despite the lesions being associated with increased water temperatures and elevated mortality rates. To address this shortfall, a transcriptomic approach was used to characterise the molecular landscape of spot lesions, the most commonly reported lesion type in New Zealand Chinook salmon, versus healthy appearing skin in fish with and without spot lesions. Many biological (gene ontology) pathways were enriched in lesion adjacent tissue, relative to control skin tissue, including proteolysis, fin regeneration, calcium ion binding, mitochondrial transport, actin cytoskeleton organisation, epithelium development, and tissue development. In terms of specific transcripts of interest, pro-inflammatory cytokines (interleukin 1β and tumour necrosis factor), annexin A1, mucin 2, and calreticulin were upregulated, while cathepsin H, mucin 5AC, and perforin 1 were downregulated in lesion tissue. In some instances, changes in gene expression were consistent between lesion and healthy appearing skin from the same fish relative to lesion free fish, suggesting that host responses weren't limited to the site of the lesion. Goblet cell density in skin histological sections was not different between skin sample types. Collectively, these results provide insights into the physiological changes associated with common spot lesions in farmed Chinook salmon.

Comparative analysis of the miRNA-mRNA regulation networks in turbot (Scophthalmus maximus L.) following Vibrio anguillarum infection

MicroRNAs could not only regulate posttranscriptional silencing of target genes in eukaryotic organisms, but also have positive effect on their target genes as well. These microRNAs have been reported to be involved in mucosal immune responses to pathogen infection in teleost. Therefore, we constructed the immune-related miRNA-mRNA networks in turbot intestine following Vibrio anguillarum infection. In our results, 1550 differentially expressed (DE) genes and 167 DE miRNAs were identified. 113 DE miRNAs targeting 89 DE mRNAs related to immune response were used to construct miRNA-mRNA interaction networks. Functional analysis showed that target genes were associated with synthesis and degradation of ketone bodies, mucin type O-Glycan biosynthesis, homologous recombination, biotin metabolism, and intestinal immune network for IgA production that were equivalent to the function of IgT and IgM in fish intestine. Finally, 10 DE miRNAs and 7 DE mRNAs were selected for validating the accuracy of high-throughput sequencing results by qRT-PCR. The results of this study will provide valuable information for the elucidation of the regulation mechanisms of miRNA-mRNA interactions involved in disease resistance in teleost mucosal immune system.

Elucidating bacterial adhesion to mucosal surface by an original AFM approach

Background

Fish skin represents an ancient vertebrate mucosal surface, sharing characteristics with other mucosal surfaces including those of the intestine. The skin mucosa is continuously exposed to microbes in the surrounding water and is therefore important in the first line defense against environmental pathogens by preventing bacteria from accessing the underlying surfaces. Understanding the microbe-host interactions at the fish skin mucosa is highly relevant in order to understand and control infection, commensalism, colonization, persistence, infection, and disease. Here we investigate the interactions between the pathogenic bacteria Aeromonas salmonicida (A. salmonicida) and Yersinia ruckeri (Y. ruckeri), respectively, and the skin mucosal surface of Atlantic salmon fry using AFM force spectroscopy.

Results

The results obtained revealed that when retracting probes functionalized with bacteria from surfaces coated with immobilized mucins, isolated from salmon mucosal surfaces, rupture events reflecting the disruption of adhesive interactions were observed, with rupture strengths centered around 200 pN. However, when retracting probes functionalized with bacteria from the intact mucosal surface of salmon fish fry no adhesive interactions could be detected. Furthermore, rheological measurements revealed a near fluid-like behavior for the fish fry skin mucus. Taken together, the experimental data indicate that the adhesion between the mucin molecules within the mucous layer may be significantly weaker than the interaction between the bacteria and the mucin molecules. The bacteria, immobilized on the AFM probe, do bind to individual mucins in the mucosal layer, but are released from the near fluid mucus with little resistance upon retraction of the AFM probe, to which they are immobilized.

Conclusion

The data provided in the current paper reveal that A. salmonicida and Y. ruckeri do bind to the immobilized mucins. However, when retracting the bacteria from intact mucosal surfaces, no adhesive interactions are detected. These observations suggest a mechanism underlying the protective function of the mucosal surface based on the clearing of potential threats by adhering them to loosely attached mucus that is subsequently released from the fish skin.

Survival of metazoan parasites in fish: Putting into context the protective immune responses of teleost fish

Defence mechanisms of fish can be divided into specific and non-specific that act in concert and are often interdependent. Most fish in both wild and cultured populations are vulnerable to metazoan parasites. Endoparasitic helminths include several species of digeneans, cestodes, nematodes, and acanthocephalans. Although they may occur in large numbers, helminth infections rarely result in fish mortality. Conversely, some ectoparasites cause mass mortality in farmed fish. Given the importance of fish innate immunity, this review addresses non-specific defence mechanisms of fish against metazoan parasites, with emphasis on granulocyte responses involving mast cells, neutrophils, macrophages, rodlet cells, and mucous cells. Metazoan parasites are important disease agents that affect wild and farmed fish and can induce high economic loss and, as pathogen organisms, deserve considerable attention. The paper will provide our light and transmission electron microscopy data on metazoan parasites-fish innate immune and neuroendocrine systems. Insights about the structure and functions of the cell types listed above and a brief account of the effects and harms of each metazoan taxon to specific fish apparati/organs will be presented.

Stress Impairs Skin Barrier Function and Induces α2-3 Linked N-Acetylneuraminic Acid and Core 1 O-Glycans on Skin Mucins in Atlantic Salmon, Salmo salar

The skin barrier consists of mucus, primarily comprising highly glycosylated mucins, and the epithelium. Host mucin glycosylation governs interactions with pathogens and stress is associated with impaired epithelial barrier function. We characterized Atlantic salmon skin barrier function during chronic stress (high density) and mucin O-glycosylation changes in response to acute and chronic stress. Fish held at low (LD: 14–30 kg/m³) and high densities (HD: 50-80 kg/m³) were subjected to acute stress 24 h before sampling at 17 and 21 weeks after start of the experiment. Blood parameters indicated primary and secondary stress responses at both sampling points. At the second sampling, skin barrier function towards molecules was reduced in the HD compared to the LD group (P_app mannitol; p < 0.01). Liquid chromatography–mass spectrometry revealed 81 O-glycan structures from the skin. Fish subjected to both chronic and acute stress had an increased proportion of large O-glycan structures. Overall, four of the O-glycan changes have potential as indicators of stress, especially for the combined chronic and acute stress. Stress thus impairs skin barrier function and induces glycosylation changes, which have potential to both affect interactions with pathogens and serve as stress indicators.

Skin transcriptome, tissue distribution of mucin genes and discovery of simple sequence repeats in crucian carp (Carassius auratus)

Crucian carp (Carassius auratus) is one of the major freshwater species and important food fish in China. Fish skin acts as the first line of defence against pathogens, yet its molecular and immune mechanism remains unclear. In this study, a de novo transcriptome assembly of C. auratus skin was performed with the Illumina Hiseq 2000 platform. A total of 49,154,776 clean reads were assembled, among which 60,824 (46.86%), 37,103 (28.59%), 43,269 (33.33%) unigenes were annotated against National Center for Biotechnology Information, Gene Onotology and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases, respectively. KEGG Orthology categories were significantly involved in immune system (20.50%), signal transduction (18.04%) and mucosal mucin genes (e.g., muc2, muc5AC, muc5B, muc17, muc18). The high expression of muc18 gene was observed in brain; that of muc2 in intestine; and that of muc5AC in skin, liver, spleen, intestine and muscle. Moreover, the potential 28,928 simple sequence repeats with the three most abundant dinucleotide repeat motifs (AC/GT, AG/CT, AT/AT) were detected in C. auratus. To authors' knowledge, this is the first report to describe the transcriptome analysis of C. auratus skin, and the outcome of this study contributed to the understanding of mucosal immune response of the skin and molecular markers in cyprinid species.

Effect of dietary isoleucine on skin mucus barrier and epithelial physical barrier functions of hybrid bagrid catfish Pelteobagrus vachelli × Leiocassis longirostris

The study investigated the effects of dietary isoleucine (Ile) on skin mucus barrier and epithelial physical barrier functions of hybrid bagrid catfish Pelteobagrus vachelli × Leiocassis longirostris. A total of 630 fish (33.11 ± 0.09 g) were fed semi-purified isonitrogenous diets containing 5.0 (control), 7.5, 10.0, 12.5, 15.0, 17.5, and 20.0 g Ile kg ^-1 diet for 8 weeks. The results indicated that dietary Ile increased (P < 0.05) in skin (1) mucus protein content and antimicrobial activity against three gram-negative bacteria (Aeromonas hydrophila, Escherichia coli, and Yersinia ruckeri) and two gram-positive bacteria (Streptococcus agalactiae and Staphylococcus aureus), (2) mucus lysofew information is available about the influencezyme (LZM), acid phosphatase (ACP), and alkaline phosphatase (AKP) activities, and complement 3 and 4 (C3 and C4) and immunoglobulin M (IgM) contents, (3) intelectin 1 (intl1), intelectin 2 (intl2), c-type-lysozyme (c-LZM), g-type-lysozyme (g-LZM), and β-defensin mRNA levels. Dietary Ile decreased (P < 0.05) reactive oxygen species (ROS), malondialdehyde (MDA), and protein carbonyl (PC) contents, and up-regulated (P < 0.05) CuZnSOD, GST, GPX1a, muc5ac, muc5b, zonula occludens-1 (ZO-1), zonula occludens-2 (ZO-2), occludin, and claudin 3 mRNA levels in skin. These results indicated that Ile improved skin mucus barrier function via increasing mucus protein, C3 and C4, and IgM contents and antibacterial factors activities, and promoted epithelial physical barrier function via decreasing skin antioxidant damage and improving tight junction structure in hybrid bagrid catfish.

Transcriptome analysis of the Larimichthys polyactis under heat and cold stress

The small yellow croaker (Larimichthys polyactis) is an economically important marine fish that is widely distributed in the East Sea, Yellow Sea and Bohai of China. However, the wild populations of L. polyactis are severely depleted, and there is currently a developing large-scale artificial propagation of this fish for aquaculture. However, the current variety of L. polyactis that is cultivated is not capable to coping with large fluctuations in temperature. Therefore, it is important to understand the molecular mechanisms that are activated in response to temperature stress in the L. polyactis. Here, we conducted transcriptomic analysis of the liver of L. polyactis under heat and cold stress. A total of 270,844,888, 265,727,006 and 259,666,218 clean reads were generated from high temperature group, low temperature group and control group, respectively, and comparing expression of genes in these transcriptomes, 10,878 unigenes that were differential expressed were identified. Seventeen of the differentially expressed unigenes were validated by qRT-PCR. Pathway enrichment analysis identified that the ER pathway, immune signaling pathway and metabolic response pathway were affected by temperature stress. The results of this study provide a comprehensive overview of temperature stress-induced transcriptional patterns in liver tissues of the L. polyactis. In addition, these results can guide future molecular studies of heat and cold stress response in this species for improving the stock used for aquaculture.

Transcriptome analyses reveal alterations in muscle metabolism, immune responses and reproductive behavior of Japanese mantis shrimp (Oratosquilla oratoria) at different cold temperature

Low temperature reduction is thought to cause widespread effects on the physical and behavioral traits of marine organisms, which include metabolic processes, immune responses, and reproductive behavior. Crustaceans are generally considered sensitive to temperature reduction due to the lack of efficient regulators. To better understand the molecular regulatory mechanisms of crustacean exposure to cold stress, Japanese mantis shrimp (Oratosquilla oratoria) was chosen as a representative crustacean. Transcriptomic responses in O. oratoria from five temperatures (25 °C, 22 °C, 19 °C, 16 °C, and 13 °C) were studied using RNA-seq. A total of 64.91 Gb of clean transcriptomic data were generated in 10 libraries and then spliced into 52,107 unigenes with an average length of 1089 bp and an N50 length of 1872 bp. A total of 14,841 unigenes was annotated in at least one database using Blastx alignment. Compared with the control temperature (25 °C), 7, 21, 58, and 236 unigenes were significantly differentially expressed at 22 °C, 19 °C, 16 °C, and 13 °C, respectively. GO analysis showed that 6, 20, 27, and 35 terms were significantly enriched at 22 °C, 19 °C, 16 °C, and 13 °C, respectively. In addition, 2, 5, 2, and 10 significant pathways were presented at 22 °C, 19 °C, 16 °C, and 13 °C, respectively. Combining NR, GO, and KEGG annotation information, many genes significantly differentially expressed at low temperatures may be associated with metabolic processes, immune response, and reproductive behavior. Additionally, we reconstructed the phylogenetic relationship based on 366 orthologous genes and the predicted differentiation time of O. oratoria and P. vannamei range from 212.82 to 365.30 Mya. Furthermore, 16 orthologous genes were identified as PSGs and 30 orthologous genes were identified as FEGs and these adaptive genes were associated with energy metabolism, stress response and immunity, and multiple cellular processing. These results provide fundamental information about molecular mechanisms regulating cold stress response of O. oratoria.

Transcription profiles of skin and head kidney from goldfish suffering hemorrhagic septicemia with an emphasis on the TLR signaling pathway

Hemorrhagic septicemia is an acute, highly fatal disease that affects goldfish (Carassius auratus). To gain a better understanding of related immune genes, the transcriptomes of the skin and head kidney of goldfish suffering hemorrhagic septicemia were sequenced, assembled, and characterized. Based on functional annotation, an extensive and diverse catalog of expressed genes were identified in both the skin and head kidney. As two different organs, pair-wise comparison identified 122/77 unigenes up/down-regulated (two-fold change with P<0.05) in the skin and head kidney. Most genes of the immune pathways were expressed and isolated in both skin and head kidney, including interferon (IFN) transcription factors 1-10 and Toll-like receptors (TLRs). Interferon regulatory factor 3 (IRF3), a key IFN transcription factor, was up-regulated at the transcriptional level by polyriboinosinic: polyribocytidylic acid (poly I:C) challenge and regulated the IFN response by increasing the activity of IFN-β and IFN-stimulated response element (ISRE)-containing promoter. This study will benefit the identification and understanding of novel genes that play important roles in the immunological reactions of fish suffering from hemorrhagic septicemia.

Comparative analysis of the transcriptome of the Amazonian fish species Colossoma macropomum (tambaqui) and hybrid tambacu by next generation sequencing

Background

The C. macropomum is a characiform fish from the Amazon basin that has been hybridized with other pacu species to produce commercial hybrids, such as the tambacu. However, little is known of the functional genomics of the parental species or these hybrid forms. The transcriptome of C. macropomum and tambacu were sequenced using 454 Roche platform (pyrosequencing) techniques to characterize the domains of Gene Ontology (GO) and to evaluate the levels of gene expression in the two organisms.

Results

The 8,188,945 reads were assembled into 400,845 contigs. A total of 58,322 contigs were annotated with a predominance of biological processes for both organisms, as determined by Gene Ontology (GO). Similar numbers of metabolic pathways were identified in both the C. macropomum and the tambacu, with the metabolism category presenting the largest number of transcripts. The BUSCO analysis indicated that our assembly was more than 40% complete. We identified 21,986 genes for the two fishes. The P and Log2FC values indicated significant differences in the levels of gene expression, with a total of 600 up-regulated genes.

Conclusion

In spite of the lack of a reference genome, the functional annotation was successful, and confirmed a considerable difference in the specificity and levels of gene expression between the two organisms. This report provides a comprehensive baseline for the genetic management of these commercially important fishes, in particular for the identification of specific genes that may represent markers involved in the immunity, growth, and fertility of these organisms, with potential practical applications in aquaculture management.

The Change of Teleost Skin Commensal Microbiota Is Associated With Skin Mucosal Transcriptomic Responses During Parasitic Infection by Ichthyophthirius multifillis

Teleost skin serves as the first line of defense against invading pathogens, and contain a skin-associated lymphoid tissue (SALT) that elicit gut-like immune responses against antigen stimulation. Moreover, exposed to the water environment and the pathogens therein, teleost skin is also known to be colonized by diverse microbial communities. However, little is known about the interactions between microbiota and the teleost skin mucosal immune system, especially dynamic changes about the interactions under pathogen infection. We hypothesized that dramatic changes of microbial communities and strong mucosal immune response would be present in the skin of aquatic vertebrate under parasite infection. To confirm this hypothesis, we construct an infected model with rainbow trout (Oncorhynchus mykiss), which was experimentally challenged by Ichthyophthirius multifiliis (Ich). H & E staining of trout skin indicates the successful invasion of Ich and shows the morphological changes caused by Ich infection. Critically, increased mRNA expression levels of immune-related genes were detected in trout skin from experimental groups using qRT-PCR, which were further studied by RNA-Seq analysis. Here, through transcriptomics, we detected that complement factors, pro-inflammatory cytokines, and antimicrobial genes were strikingly induced in the skin of infected fish. Moreover, high alpha diversity values of microbiota in trout skin from the experimental groups were discovered. Interestingly, we found that Ich infection led to a decreased abundance of skin commensals and increased colonization of opportunistic bacteria through 16S rRNA pyrosequencing, which were mainly characterized by lose of Proteobacteria and increased intensity of Flavobacteriaceae. To our knowledge, our results suggest for the first time that parasitic infection could inhibit symbionts and offer opportunities for other pathogens' secondary infection in teleost skin.

Biological and Ecological Roles of External Fish Mucus: A Review

Fish mucus layers are the main surface of exchange between fish and the environment, and they possess important biological and ecological functions. Fish mucus research is increasing rapidly, along with the development of high-throughput techniques, which allow the simultaneous study of numerous genes and molecules, enabling a deeper understanding of the fish mucus composition and its functions. Fish mucus plays a major role against fish infections, and research has mostly focused on the study of fish mucus bioactive molecules (e.g., antimicrobial peptides and immune-related molecules) and associated microbiota due to their potential in aquaculture and human medicine. However, external fish mucus surfaces also play important roles in social relationships between conspecifics (fish shoaling, spawning synchronisation, suitable habitat finding, or alarm signals) and in interspecific interactions such as prey-predator relationships, parasite–host interactions, and symbiosis. This article reviews the biological and ecological roles of external (gills and skin) fish mucus, discussing its importance in fish protection against pathogens and in intra and interspecific interactions. We also discuss the advances that “omics” sciences are bringing into the fish mucus research and their importance in studying the fish mucus composition and functions.

Gene expression analysis of Atlantic salmon gills reveals mucin 5 and interleukin 4/13 as key molecules during amoebic gill disease

Amoebic gill disease (AGD) is one of the main diseases affecting Atlantic salmon (Salmo salar L.) mariculture. Hallmarks of AGD are hyperplasia of the lamellar epithelium and increased production of gill mucus. This study investigated the expression of genes involved in mucus secretion, cell cycle regulation, immunity and oxidative stress in gills using a targeted 21-gene PCR array. Gill samples were obtained from experimental and natural Neoparamoeba perurans infections, and sampling points included progressive infection stages and post-freshwater treatment. Up-regulation of genes related to mucin secretion and cell proliferation, and down-regulation of pro-inflammatory and pro-apoptotic genes were associated with AGD severity, while partial restoration of the gill homeostasis was detected post-treatment. Mucins and Th2 cytokines accoun ted for most of the variability observed between groups highlighting their key role in AGD. Two mucins (muc5, muc18) showed differential regulation upon disease. Substantial up-regulation of the secreted muc5 was detected in clinical AGD, and the membrane bound muc18 showed an opposite pattern. Th2 cytokines, il4/13a and il4/13b2, were significantly up-regulated from 2 days post-infection onwards, and changes were lesion-specific. Despite the differences between experimental and natural infections, both yielded comparable results that underline the importance of the studied genes in the respiratory organs of fish, and during AGD progression.

Mucus-Pathogen Interactions in the Gastrointestinal Tract of Farmed Animals

Gastrointestinal infections cause significant challenges and economic losses in animal husbandry. As pathogens becoming resistant to antibiotics are a growing concern worldwide, alternative strategies to treat infections in farmed animals are necessary in order to decrease the risk to human health and increase animal health and productivity. Mucosal surfaces are the most common route used by pathogens to enter the body. The mucosal surface that lines the gastrointestinal tract is covered by a continuously secreted mucus layer that protects the epithelial surface. The mucus layer is the first barrier the pathogen must overcome for successful colonization, and is mainly composed of densely glycosylated proteins called mucins. The vast array of carbohydrate structures present on the mucins provide an important setting for host-pathogen interactions. This review summarizes the current knowledge on gastrointestinal mucins and their role during infections in farmed animals. We examine the interactions between mucins and animal pathogens, with a focus on how pathogenic bacteria can modify the mucin environment in the gut, and how this in turn affects pathogen adhesion and growth. Finally, we discuss analytical challenges and complexities of the mucus-based defense, as well as its potential to control infections in farmed animals.

Low incubation temperature during early development negatively affects survival and related innate immune processes in zebrafish larvae exposed to lipopolysaccharide

In many fish species, the immune system is significantly constrained by water temperature. In spite of its critical importance in protecting the host against pathogens, little is known about the influence of embryonic incubation temperature on the innate immunity of fish larvae. Zebrafish (Danio rerio) embryos were incubated at 24, 28 or 32 °C until first feeding. Larvae originating from each of these three temperature regimes were further distributed into three challenge temperatures and exposed to lipopolysaccharide (LPS) in a full factorial design (3 incubation × 3 challenge temperatures). At 24 h post LPS challenge, mortality of larvae incubated at 24 °C was 1.2 to 2.6-fold higher than those kept at 28 or 32 °C, regardless of the challenge temperature. LPS challenge at 24 °C stimulated similar immune-related processes but at different levels in larvae incubated at 24 or 32 °C, concomitantly with the down-regulation of some chemokine and lysozyme transcripts in the former group. Larvae incubated at 24 °C and LPS-challenged at 32 °C exhibited a limited immune response with up-regulation of hypoxia and oxidative stress processes. Annexin A2a, S100 calcium binding protein A10b and lymphocyte antigen-6, epidermis were identified as promising candidates for LPS recognition and signal transduction.

Structural and functional maturation of skin during metamorphosis in the Atlantic halibut (Hippoglossus hippoglossus)

To establish if the developmental changes in the primary barrier and osmoregulatory capacity of Atlantic halibut skin are modified during metamorphosis, histological, histochemical, gene expression and electrophysiological measurements were made. The morphology of the ocular and abocular skin started to diverge during the metamorphic climax and ocular skin appeared thicker and more stratified. Neutral mucins were the main glycoproteins produced by the goblet cells in skin during metamorphosis. Moreover, the number of goblet cells producing neutral mucins increased during metamorphosis and asymmetry in their abundance was observed between ocular and abocular skin. The increase in goblet cell number and their asymmetric abundance in skin was concomitant with the period that thyroid hormones (THs) increase and suggests that they may be under the control of these hormones. Several mucin transcripts were identified in metamorphosing halibut transcriptomes and Muc18 and Muc5AC were characteristic of the body skin. Na⁺, K⁺-ATPase positive (NKA) cells were observed in skin of all metamorphic stages but their number significantly decreased with the onset of metamorphosis. No asymmetry was observed between ocular and abocular skin in NKA cells. The morphological changes observed were linked to modified skin barrier function as revealed by modifications in its electrophysiological properties. However, the maturation of the skin functional characteristics preceded structural maturation and occurred at stage 8 prior to the metamorphic climax. Treatment of Atlantic halibut with the THs disrupter methimazole (MMI) affected the number of goblet cells producing neutral mucins and the NKA cells. The present study reveals that the asymmetric development of the skin in Atlantic halibut is TH sensitive and is associated with metamorphosis and that this barrier's functional properties mature earlier and are independent of metamorphosis.

Genome-wide analysis of Atlantic salmon (Salmo salar) mucin genes and their role as biomarkers

The aim of this study was to identify potential mucin genes in the Atlantic salmon genome and evaluate tissue-specific distribution and transcriptional regulation in response to aquaculture-relevant stress conditions in post-smolts. Seven secreted gel-forming mucin genes were identified based on several layers of evidence; annotation, transcription, phylogeny and domain structure. Two genes were annotated as muc2 and five genes as muc5. The muc2 genes were predominantly transcribed in the intestinal region while the different genes in the muc5 family were mainly transcribed in either skin, gill or pyloric caeca. In order to investigate transcriptional regulation of mucins during stress conditions, two controlled experiments were conducted. In the first experiment, handling stress induced mucin transcription in the gill, while transcription decreased in the skin and intestine. In the second experiment, long term intensive rearing conditions (fish biomass ~125 kg/m3) interrupted by additional confinement led to increased transcription of mucin genes in the skin at one, seven and fourteen days post-confinement.

Toxic Effects and Molecular Mechanism of Different Types of Silver Nanoparticles to the Aquatic Crustacean Daphnia magna

Silver nanoparticles (AgNPs) have been assessed to have a high exposure risk for humans and aquatic organisms. Toxicity varies considerably between different types of AgNPs. This study aimed to investigate the toxic effects of AgNPs with different particle sizes (40 and 110 nm) and different surface coatings (sodium citrate and polyvinylpyrrolidone, PVP) on Daphnia magna and their mechanisms of action. The results revealed that the citrate-coated AgNPs were more toxic than PVP-coated AgNPs and that the 40 nm AgNPs were more toxic than the 110 nm AgNPs. Transcriptome analysis further revealed that the toxic effects of AgNPs on D. magna were related to the mechanisms of ion binding and several metabolic pathways, such as the "RNA polymerase" pathway and the "protein digestion and absorption" pathway. Moreover, the principal component analysis (PAC) results found that surface coating was the major factor that determines the toxicities compared to particle size. These results could help us better understand the possible mechanism of AgNP toxicity in aquatic invertebrates at the transcriptome level and establish an important foundation for revealing the broad impacts of nanoparticles on aquatic environments.

Omics in fish mucosal immunity

The mucosal immune system of fish is a complex network of immune cells and molecules that are constantly surveilling the environment and protecting the host from infection. A number of "omics" tools are now available and utilized to understand the complexity of mucosal immune systems in non-traditional animal models. This review summarizes recent advances in the implementation of "omics" tools pertaining to the four mucosa-associated lymphoid tissues in teleosts. Genomics, transcriptomics, proteomics, and "omics" in microbiome research require interdisciplinary collaboration and careful experimental design. The data-rich datasets generated are proving really useful at discovering new innate immune players in fish mucosal secretions, identifying novel markers of specific mucosal immune responses, unraveling the diversity of the B and T cell repertoires and characterizing the diversity of the microbial communities present in teleost mucosal surfaces. Bioinformatics, data analysis and storage platforms should be developed to facilitate rapid processing of large datasets, especially when mammalian tools such as bioinformatics analysis software are not available in fishes.

Complexities of gene expression patterns in natural populations of an extremophile fish (Poecilia mexicana, Poeciliidae)

Variation in gene expression can provide insights into organismal responses to environmental stress and physiological mechanisms mediating adaptation to habitats with contrasting environmental conditions. We performed an RNA-sequencing experiment to quantify gene expression patterns in fish adapted to habitats with different combinations of environmental stressors, including the presence of toxic hydrogen sulphide (H2 S) and the absence of light in caves. We specifically asked how gene expression varies among populations living in different habitats, whether population differences were consistent among organs, and whether there is evidence for shared expression responses in populations exposed to the same stressors. We analysed organ-specific transcriptome-wide data from four ecotypes of Poecilia mexicana (nonsulphidic surface, sulphidic surface, nonsulphidic cave and sulphidic cave). The majority of variation in gene expression was correlated with organ type, and the presence of specific environmental stressors elicited unique expression differences among organs. Shared patterns of gene expression between populations exposed to the same environmental stressors increased with levels of organismal organization (from transcript to gene to physiological pathway). In addition, shared patterns of gene expression were more common between populations from sulphidic than populations from cave habitats, potentially indicating that physiochemical stressors with clear biochemical consequences can constrain the diversity of adaptive solutions that mitigate their adverse effects. Overall, our analyses provided insights into transcriptional variation in a unique system, in which adaptation to H2 S and darkness coincide. Functional annotations of differentially expressed genes provide a springboard for investigating physiological mechanisms putatively underlying adaptation to extreme environments.

Aeromonas salmonicida Growth in Response to Atlantic Salmon Mucins Differs between Epithelial Sites, Is Governed by Sialylated and N-Acetylhexosamine-Containing O-Glycans, and Is Affected by Ca2

Aeromonas salmonicida causes furunculosis in salmonids and is a threat to Atlantic salmon aquaculture. The epithelial surfaces that the pathogen colonizes are covered by a mucus layer predominantly comprised of secreted mucins. By using mass spectrometry to identify mucin glycan structures with and without enzymatic removal of glycan residues, coupled to measurements of bacterial growth, we show here that the complex Atlantic salmon intestinal mucin glycans enhance A. salmonicida growth, whereas the more simple skin mucin glycans do not. Of the glycan residues present terminally on the salmon mucins, only N-acetylglucosamine (GlcNAc) enhances growth. Sialic acids, which have an abundance of 75% among terminal glycans from skin and of <50% among intestinal glycans, cannot be removed or used by A. salmonicida for growth-enhancing purposes, and they shield internal GlcNAc from utilization. A Ca²⁺ concentration above 0.1 mM is needed for A. salmonicida to be able to utilize mucins for growth-promoting purposes, and 10 mM further enhances both A. salmonicida growth in response to mucins and binding of the bacterium to mucins. In conclusion, GlcNAc and sialic acids are important determinants of the A. salmonicida interaction with its host at the mucosal surface. Furthermore, since the mucin glycan repertoire affects pathogen growth, the glycan repertoire may be a factor to take into account during breeding and selection of strains for aquaculture.

Transcriptomic variation of locally-infected skin of Epinephelus coioides reveals the mucosal immune mechanism against Cryptocaryon irritans

Fish skin is the largest immunologically active mucosal organ, providing first-line defense against external pathogens. However, the skin-associated immune mechanisms of fish are still unclear. Cryptocaryon irritans is an obligate ectoparasitic ciliated protozoan that infects almost all marine fish, and is believed to be an excellent pathogen model to study fish mucosal immunity. In this study, a de novo transcriptome assembly of Epinephelus coioides skin post C. irritans tail-infection was performed for the first time using the Illumina HiSeq™ 2500 system. Comparative analyses of infected skin (group Isk) and uninfected skin (group Nsk) from the same challenged fish and control skin (group C) from uninfected control fish were conducted. As a result, a total of 91,082 unigenes with an average length of 2880 base pairs were obtained and among them, 38,704 and 48,617 unigenes were annotated based on homology with matches in the non-redundant and zebrafish database, respectively. Pairwise comparison resulted in 10,115 differentially-expressed genes (DEGs) in the Isk/C group comparison (4,983 up-regulated and 5,132 down-regulated), 2,275 DEGs in the Isk/Nsk group comparison (1,319 up-regulated and 956 down-regulated) and 4,566 DEGs in the Nsk/C group comparison (1,534 up-regulated and 3,032 down-regulated). Seven immune-related categories including 91 differentially-expressed immune genes (86 up-regulated and 5 down-regulated) were scrutinized. Both DEGs and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis and immune-related gene expression analysis were used, and both analyses showed that the genes were more significantly altered in the locally-infected skin than in the uninfected skin of the same challenged fish. This suggests the skin's local immune response is important for host defense against this ectoparasite infection. Innate immune molecules, including hepcidin, C-type lectin, transferrin, transferrin receptor protein, serum amyloid A, cathepsin and complement components were significantly up-regulated (fold-change ranged from 3.3 to 12,944) in infected skin compared with control skin. The up-regulation of chemokines and chemokine receptors and activation of the leukocyte transendothelial migration pathway suggested that leucocytes intensively migrated to the local infected sites to mount a local immune defense. Toll-like receptors (TLRs) 1, 2, 5 and 5S were most significantly up-regulated in the infected skin, suggesting that these TLRs may be involved in parasite pathogen-associated molecular pattern (PAMPs) recognition. Up-regulation of the dendritic cell markers CD209 and CD83 and other antigen presentation pathway molecules provided evidence for skin local antigen presentation. Up-regulation of the T cell markers CD4 and CD48, B cell markers CD22 and CD81 and B cell receptor signaling kinase Lyn, showed the presence and population expansion of T/B cells at locally-infected sites, which suggested possible activation of a local specific immune response in the skin. Our results will facilitate in-depth understanding of local immune defense mechanisms in fish skin against ectoparasite infection.

De novo assembly of the sea trout (Salmo trutta m. trutta) skin transcriptome to identify putative genes involved in the immune response and epidermal mucus secretion

In fish, the skin is a multifunctional organ and the first barrier against pathogens. Salmonids differ in their susceptibility to microorganisms due to varied skin morphology and gene expression patterns. The brown trout is a salmonid species with important commercial and ecological value in Europe. However, there is a lack of knowledge regarding the genes involved in the immune response and mucus secretion in the skin of this fish. Thus, we characterized the skin transcriptome of anadromous brown trout using next-generation sequencing (NGS). A total of 1,348,306 filtered reads were obtained and assembled into 75,970 contigs. Of these contigs 48.57% were identified using BLAST tool searches against four public databases. KEGG pathway and Gene Ontology analyses revealed that 13.40% and 34.57% of the annotated transcripts, respectively, represent a variety of biological processes and functions. Among the identified KEGG Orthology categories, the best represented were signal transduction (23.28%) and immune system (8.82%), with a variety of genes involved in immune pathways, implying the differentiation of immune responses in the trout skin. We also identified and transcriptionally characterized 8 types of mucin proteins–the main structural components of the mucosal layer. Moreover, 140 genes involved in mucin synthesis were identified, and 1,119 potential simple sequence repeats (SSRs) were detected in 3,134 transcripts.

Identification and characterization of the expression profile of the microRNAs in the Amazon species Colossoma macropomum by next generation sequencing

Colossoma macropomum is a resistant species native of Amazonas and Orinoco river basins. It is regarded as the second largest finfish of Solimões and Amazon rivers, representing a major fishery resource in Amazonas and an important species in tropical aquaculture. MicroRNAs are non-coding endogenous riboregulators of nearly 22 nucleotides that play a key role in post-transcriptional gene regulation of several organisms. We analyzed samples of liver and skin from specimens of C. macropomum using next generation sequencing. The dataset was evaluated using computational programs to check the quality of sequences, identification of miRNAs, as well as to evaluate the expression levels of these microRNAs and interaction of target genes. We identified 279 conserved miRNAs, being 257 from liver and 272 from skin, with several miRNAs shared between tissues, with divergence in the number of reads. The strands miR-5p and miR-3p were observed in 72 miRNAs, some of them presenting a higher number of 3p reads. The functional annotation of the most expressed miRNAs resulted in 27 pathways for the liver and skin mainly related to the "biological processes" domain. Based on the identified pathways, we visualized a large gene network, suggesting the regulation of selected miRNA over this interactive dataset. We were able to identify and characterize the expression levels of miRNAs in two tissues of great activity in C. macropomum, which stands out as the beginning of several studies that can be carried out to elucidate the influence of miRNAs in this species and their applicability as biotechnological tools.

Transcriptomic responses in the fish intestine

The intestine, being a multifunctional organ central to both nutrient uptake, pathogen recognition and regulating the intestinal microbiome, has been subjected to intense research. This review will focus on the recent studies carried out using high-throughput gene expression approaches, such as microarray and RNA sequencing (RNA-seq). These techniques have advanced greatly in recent years, mainly as a result of the massive changes in sequencing methodologies. At the time of writing, there is a transition between relatively well characterised microarray platforms and the developing RNA-seq, with the prediction that within a few years as costs decrease and computation power increase, RNA-seq related approaches will supersede the microarrays. Comparisons between the approaches are made and specific examples of how the techniques have been used to examine intestinal responses to pathogens, dietary manipulations and osmoregulatory challenges are given.

Sequencing, De Novo Assembly, and Annotation of the Transcriptome of the Endangered Freshwater Pearl Bivalve, Cristaria plicata, Provides Novel Insights into Functional Genes and Marker Discovery

BACKGROUND

The freshwater mussel Cristaria plicata (Bivalvia: Eulamellibranchia: Unionidae), is an economically important species in molluscan aquaculture due to its use in pearl farming. The species have been listed as endangered in South Korea due to the loss of natural habitats caused by anthropogenic activities. The decreasing population and a lack of genomic information on the species is concerning for environmentalists and conservationists. In this study, we conducted a de novo transcriptome sequencing and annotation analysis of C. plicata using Illumina HiSeq 2500 next-generation sequencing (NGS) technology, the Trinity assembler, and bioinformatics databases to prepare a sustainable resource for the identification of candidate genes involved in immunity, defense, and reproduction.

RESULTS

The C. plicata transcriptome analysis included a total of 286,152,584 raw reads and 281,322,837 clean reads. The de novo assembly identified a total of 453,931 contigs and 374,794 non-redundant unigenes with average lengths of 731.2 and 737.1 bp, respectively. Furthermore, 100% coverage of C. plicata mitochondrial genes within two unigenes supported the quality of the assembler. In total, 84,274 unigenes showed homology to entries in at least one database, and 23,246 unigenes were allocated to one or more Gene Ontology (GO) terms. The most prominent GO biological process, cellular component, and molecular function categories (level 2) were cellular process, membrane, and binding, respectively. A total of 4,776 unigenes were mapped to 123 biological pathways in the KEGG database. Based on the GO terms and KEGG annotation, the unigenes were suggested to be involved in immunity, stress responses, sex-determination, and reproduction. A total of 17,251 cDNA simple sequence repeats (cSSRs) were identified from 61,141 unigenes (size of >1 kb) with the most abundant being dinucleotide repeats.

CONCLUSIONS

This dataset represents the first transcriptome analysis of the endangered mollusc, C. plicata. The transcriptome provides a comprehensive sequence resource for the conservation of genetic information in this species and enrichment of the genetic database. The development of molecular markers will assist in the genetic improvement of C. plicata.

De Novo Assembly and Characterization of Early Embryonic Transcriptome of the Horseshoe Crab Tachypleus tridentatus

The horseshoe crab Tachypleus tridentatus is a unique marine species and a potential model for marine invertebrate. Limited genomic and transcriptional data are currently available to understand the molecular mechanisms underlying the embryonic development of T. tridentatus. Here, we reported for the first time the de novo transcriptome assembly for T. tridentatus at embryonic developmental stage using Illumina RNA-seq platform. Approximate 38 million reads were obtained and further assembled into 133,212 unigenes. Sequence homology analysis against public databases revealed that 33,796 unigenes could be annotated with gene descriptions. Of the annotated unigenes, we identified a number of key components of several conserved metazoan signaling pathways (Hedgehog, Wnt, TGF-beta and Notch pathways) and other important regulatory genes involved in embryonic development. Targeted searching of Pax family genes which play critical roles in the formation of tissue and organ during embryonic development identified a complete set of Pax family genes. Moreover, the full length T. tridentatus Pax1/9a (TtPax1/9a) and Pax1/9b (TtPax1/9b) cDNA sequences were determined based on the transcriptome, demonstrating the immediate application of our database. Using quantitative real time PCR, we analyzed the expression patterns of TtPax1/9a and TtPax1/9b in different tissues of horseshoe crab. Taking advantage of Drosophila model, we further found that TtPax1/9b, but not TtPax1/9a, can partly rescue the Drosophila homolog Poxm dysfunction-caused lethality at the larval stage. Our study provides the embryonic transcriptome of T. tridentatus which could be immediately used for gene discovery and characterization, functional genomics studies in T. tridentatus. This transcriptome database will also facilitate the investigations of molecular mechanisms underlying embryonic development of T. tridentatus and other marine arthropods as well.

Effect of temperature and diet on wound healing in Atlantic salmon (Salmo salar L.)

Compromised skin integrity of farmed Atlantic salmon, commonly occurring under low temperature and stressful conditions, has major impacts on animal welfare and economic productivity. Even fish with minimal scale loss and minor wounds can suffer from secondary infections, causing downgrading and mortalities. Wound healing is a complex process, where water temperature and nutrition play key roles. In this study, Atlantic salmon (260 g) were held at different water temperatures (4 or 12 °C) and fed three different diets for 10 weeks, before artificial wounds were inflicted and the wound healing process monitored for 2 weeks. The fish were fed either a control diet, a diet supplemented with zinc (Zn) or a diet containing a combination of functional ingredients in addition to Zn. The effect of diet was assessed through subjective and quantitative skin histology and the transcription of skin-associated chemokines. Histology confirmed that wound healing was faster at 12 °C. The epidermis was more organised, and image analyses of digitised skin slides showed that fish fed diets with added Zn had a significantly larger area of the epidermis covered by mucous cells in the deeper layers after 2 weeks, representing more advanced healing progression. Constitutive levels of the newly described chemokines, herein named CK 11A, B and C, confirmed their preferential expression in skin compared to other tissues. Contrasting modulation profiles at 4 and 12 °C were seen for all three chemokines during the wound healing time course, while the Zn-supplemented diets significantly increased the expression of CK 11A and B during the first 24 h of the healing phase.

Investigating the underlying mechanisms of temperature-related skin diseases in Atlantic salmon, Salmo salar L., as measured by quantitative histology, skin transcriptomics and composition

Skin integrity is recognized as of vital consideration for both animal welfare and final product quality of farmed fish. This study examines the effects of three different rearing temperatures (4, 10 and 16 °C) on the skin of healthy Atlantic salmon post-smolts. Changes in skin condition were assessed by the means of skin composition analyses, quantitative histology assessments and transcriptome analysis. Level of protein, vitamin C and vitamin E was significantly higher at 16 °C compared with 4 °C. Quantitative histology measurements showed that the epidermal thickness decreased from low to high temperature, whereas the epidermal area comprising mucous cells increased. The difference was only significant between 4 and 16 °C. Both high and low temperature exhibited significant changes in the skin transcriptome. A number of immune-related transcripts responded at both temperatures. Contrary to well-described immunosuppressive effects of low water temperature on systemic immunity, a subtle increase in skin-mediated immunity was observed, suggesting a pre-activation of the mucosal system at 4 °C. Upregulation of a number of heat-shock proteins correlating with a decrease in epidermal thickness suggested a stress response in the skin at high temperature. The results demonstrate distinctive temperature-related effects on the skin of Atlantic salmon.

Transcriptomic Analysis of the Adaptation of Listeria monocytogenes to Growth on Vacuum-Packed Cold Smoked Salmon

The foodborne pathogen Listeria monocytogenes is able to survive and grow in ready-to-eat foods, in which it is likely to experience a number of environmental stresses due to refrigerated storage and the physicochemical properties of the food. Little is known about the specific molecular mechanisms underlying survival and growth of L. monocytogenes under different complex conditions on/in specific food matrices. Transcriptome sequencing (RNA-seq) was used to understand the transcriptional landscape of L. monocytogenes strain H7858 grown on cold smoked salmon (CSS; water phase salt, 4.65%; pH 6.1) relative to that in modified brain heart infusion broth (MBHIB; water phase salt, 4.65%; pH 6.1) at 7°C. Significant differential transcription of 149 genes was observed (false-discovery rate [FDR], <0.05; fold change, ≥2.5), and 88 and 61 genes were up- and downregulated, respectively, in H7858 grown on CSS relative to the genes in H7858 grown in MBHIB. In spite of these differences in transcriptomes under these two conditions, growth parameters for L. monocytogenes were not significantly different between CSS and MBHIB, indicating that the transcriptomic differences reflect how L. monocytogenes is able to facilitate growth under these different conditions. Differential expression analysis and Gene Ontology enrichment analysis indicated that genes encoding proteins involved in cobalamin biosynthesis as well as ethanolamine and 1,2-propanediol utilization have significantly higher transcript levels in H7858 grown on CSS than in that grown in MBHIB. Our data identify specific transcriptional profiles of L. monocytogenes growing on vacuum-packaged CSS, which may provide targets for the development of novel and improved strategies to control L. monocytogenes growth on this ready-to-eat food.

Atlantic Salmon Carries a Range of Novel O-Glycan Structures Differentially Localized on Skin and Intestinal Mucins

Aquaculture is a growing industry, increasing the need for understanding host-pathogen interactions in fish. The skin and mucosal surfaces, covered by a mucus layer composed of mucins, is the first point of contact between fish and pathogens. Highly O-glycosylated mucins have been shown to be an important part of the defense against pathogens, and pathogens bind to host surfaces using lectin-like adhesins. However, knowledge of piscine O-glycosylation is very limited. We characterized mucin O-glycosylation of five freshwater acclimated Atlantic salmon, using mass spectrometry. Of the 109 O-glycans found, most were sialylated and differed in distribution among skin, pyloric ceca, and proximal and distal intestine. Skin O-glycans were shorter (2-6 residues) and less diverse (33 structures) than intestinal O-glycans (2-13 residues, 93 structures). Skin mucins carried O-glycan cores 1, 2, 3, and 5 and three types of sialic acids (Neu5Ac, Neu5Gc, and Kdn) and had sialyl-Tn as the predominant structure. Intestinal mucins carried only cores 1, 2, and 5, Neu5Ac was the only sialic acid present, and sialylated core 5 was the most dominant structure. This structural characterization can be used for identifying structures of putative importance in host-pathogen interactions for further testing in biological assays and disease intervention therapies.

The developmental transcriptome of contrasting Arctic charr (Salvelinus alpinus) morphs

Species and populations with parallel evolution of specific traits can help illuminate how predictable adaptations and divergence are at the molecular and developmental level. Following the last glacial period, dwarfism and specialized bottom feeding morphology evolved rapidly in several landlocked Arctic charrSalvelinus alpinuspopulations in Iceland.  To study the genetic divergence between small benthic morphs and limnetic morphs, we conducted RNA-sequencing charr embryos at four stages in early development. We studied two stocks with contrasting morphologies: the small benthic (SB) charr from Lake Thingvallavatn and Holar aquaculture (AC) charr.The data reveal significant differences in expression of several biological pathways during charr development. There was also an expression difference between SB- and AC-charr in genes involved in energy metabolism and blood coagulation genes. We confirmed differing expression of five genes in whole embryos with qPCR, includinglysozymeandnatterin-likewhich was previously identified as a fish-toxin of a lectin family that may be a putative immunopeptide. We also verified differential expression of 7 genes in the developing head that associated consistently with benthic v.s.limnetic morphology (studied in 4 morphs). Comparison of single nucleotide polymorphism (SNP) frequencies reveals extensive genetic differentiation between the SB and AC-charr (~1300 with more than 50% frequency difference). Curiously, three derived alleles in the otherwise conserved 12s and 16s mitochondrial ribosomal RNA genes are found in benthic charr.The data implicate multiple genes and molecular pathways in divergence of small benthic charr and/or the response of aquaculture charr to domestication. Functional, genetic and population genetic studies on more freshwater and anadromous populations are needed to confirm the specific loci and mutations relating to specific ecological traits in Arctic charr.

Transcriptomic analysis and biomarkers (Rag1 and Igμ) for probing the immune system development in Pacific cod, Gadus macrocephalus

Mortality (>90%) is a big concern in larval rearing facilities of Pacific cod, Gadus macrocephalus, limiting its culture presently still in the experimental stages. Understanding the immune system development of G. macrocephalus is crucial to optimize the aquaculture of this species, to improve the use of economic resources and to avoid abuse of antibiotics. For the transcriptome analysis, using an Illumina sequencing platform, 61,775,698 raw reads were acquired. After a de novo assembly, 77,561 unigenes were obtained. We have classified functionally these transcripts by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). 27 genes mainly related to hematopoietic or lymphoid organ development and somatic diversification of immune receptors have been reported for the first time in Pacific cod, and 14 Ig heavy chain (μ chain) locuses were assembled using Trinity. Based on our previous achievement, we have chosen Rag1 and Igμ as immune system development biomarkers. Full length cDNA of Rag1 and Igμ as biomarkers were obtained respectively using RACE PCR. Concerning Rag1, the deduced amino acid of Rag1 and protein immunodetection revealed a Rag1 isoform of 69 kDa, significantly different from other fish orthologs, such as Oncorhynchus mykiss (121 kDa). Phylogenetic analysis reveals a unique immune system for the Gadus genre, not exclusive for Atlantic cod, among vertebrates. Meanwhile, full length cDNA of Igμ included an ORF of 1710 bp and the deduced amino acid was composed of a leader peptide, a variable domain, CH1, CH2, Hinge, CH3, CH4 and C-terminus, which was in accordance with most teleost. Absolute quantification PCR revealed that significant expression of Rag1 appeared earlier than Igμ, 61 and 95 dph compared to 95 dph, respectively. Here we report the first transcriptomic analysis of G. macrocephalus as the starting point for genetic research on immune system development towards improving the Pacific cod aquaculture.

The Listeria monocytogenes strain 10403S BioCyc database

Listeria monocytogenes is a food-borne pathogen of humans and other animals. The striking ability to survive several stresses usually used for food preservation makes L. monocytogenes one of the biggest concerns to the food industry, while the high mortality of listeriosis in specific groups of humans makes it a great concern for public health. Previous studies have shown that a regulatory network involving alternative sigma (σ) factors and transcription factors is pivotal to stress survival. However, few studies have evaluated at the metabolic networks controlled by these regulatory mechanisms. The L. monocytogenes BioCyc database uses the strain 10403S as a model. Computer-generated initial annotation for all genes also allowed for identification, annotation and display of predicted reactions and pathways carried out by a single cell. Further ongoing manual curation based on published data as well as database mining for selected genes allowed the more refined annotation of functions, which, in turn, allowed for annotation of new pathways and fine-tuning of previously defined pathways to more L. monocytogenes-specific pathways. Using RNA-Seq data, several transcription start sites and promoter regions were mapped to the 10403S genome and annotated within the database. Additionally, the identification of promoter regions and a comprehensive review of available literature allowed the annotation of several regulatory interactions involving σ factors and transcription factors. The L. monocytogenes 10403S BioCyc database is a new resource for researchers studying Listeria and related organisms. It allows users to (i) have a comprehensive view of all reactions and pathways predicted to take place within the cell in the cellular overview, as well as to (ii) upload their own data, such as differential expression data, to visualize the data in the scope of predicted pathways and regulatory networks and to carry on enrichment analyses using several different annotations available within the database.

Transcriptome analysis revealed changes of multiple genes involved in immunity in Cynoglossus semilaevis during Vibrio anguillarum infection

Half-smooth tongue sole (Cynoglossus semilaevis) is one of the most valuable marine aquatic species in Northern China. Given to the rapid development of aquaculture industry, the C. semilaevis was subjected to disease-causing bacteria Vibrio anguillarum. It therefore is indispensable and urgent to understand the mechanism of C. semilaevis host defense against V. anguillarum infection. In the present study, the extensively analysis at the transcriptome level for V. Anguillarum disease in tongue sole was carried out. In total, 94,716 high quality contigs were generated from 75,884,572 clean reads in three libraries (HOSG, NOSG, and CG). 22,746 unigenes were identified when compared with SwissProt, an NR protein database and NT nucleotide database. 954 genes exhibiting the differentially expression at least one pair of comparison in all three libraries were identified. GO enrichment for these genes revealed gene response to biotic stimulus, immune system regulation, and immune response and cytokine production. Further, the pathways such as complement and coagulation cascades and Vibrio cholerae infection pathways were enriched in defensing of pathogen. Besides, 13,428 SSRs and 118,239 SNPs were detected in tongue sole, providing further support for genetic variation and marker-assisted selection in future. In summary, this study identifies several putative immune pathways and candidate genes deserving further investigation in the context of development of therapeutic regimens and lays the foundation for selecting resistant lines of C. semilaevis against V. anguillarum.

Sexual maturation and administration of 17β-estradiol and testosterone induce complex gene expression changes in skin and increase resistance of Atlantic salmon to ectoparasite salmon louse

The crustacean ectoparasitic salmon louse (Lepeophtheirus salmonis) is a major problem of Atlantic salmon aquaculture in the Northern hemisphere. Host-pathogen interactions in this system are highly complex. Resistance to the parasite involves variations in genetic background, nutrition, properties of skin, and status of the endocrine and immune systems. This study addressed the relationship between sex hormones and lice infection. Field observation revealed a sharp reduction of lice prevalence during sexual maturation with no difference between male and female fish. To determine if higher resistance against lice was related to sex hormones, post-smolt salmon were administered control feed and feeds containing 17β-estradiol (20 mg/kg) and testosterone (25 mg/kg) during a 3-week pre-challenge period. After challenge with lice, counts were reduced 2-fold and 1.5-fold in fish that received 17β-estradiol and testosterone, respectively. Gene expression analyses were performed from skin of salmon collected in the field trial and from the controlled lab experiment at three time points (end of feeding-before challenge, 3 days post challenge (dpc) and 16 dpc) using oligonucleotide microarray and qPCR. Differential expression was observed in genes associated with diverse biological processes. Both studies revealed similar changes of several antibacterial acute phase proteins; of note was induction of cathelicidin and down-regulation of a defensin gene. Treatment with hormones revealed their ability to modulate T helper cell (Th)-mediated immunity in skin. Enhanced protection achieved by 17β-estradiol administration might in part be due to the skewing of Th responses away from the prototypic anti-parasitic Th2 immunity and towards the more effective Th1 responses. Multiple genes involved in wound healing, differentiation and remodelling of skin tissue were stimulated during maturation but suppressed with sex hormones. Such opposite regulation suggested that these processes were not associated with resistance to the parasite under the studied conditions. Both studies revealed regulation of a suite of genes encoding putative large mucosal proteins found exclusively in fish. Marked decrease of erythrocyte markers indicated reduced circulation while down-regulation of multiple zymogen granule membrane proteins and transporters of cholesterol and other compounds suggested limited availability of nutrients for the parasites.

Aeromonas salmonicida binds differentially to mucins isolated from skin and intestinal regions of Atlantic salmon in an N-acetylneuraminic acid-dependent manner

Aeromonas salmonicida subsp. salmonicida infection, also known as furunculosis disease, is associated with high morbidity and mortality in salmonid aquaculture. The first line of defense the pathogen encounters is the mucus layer, which is predominantly comprised of secreted mucins. Here we isolated and characterized mucins from the skin and intestinal tract of healthy Atlantic salmon and studied how A. salmonicida bound to them. The mucins from the skin, pyloric ceca, and proximal and distal intestine mainly consisted of mucins soluble in chaotropic agents. The mucin density and mucin glycan chain length from the skin were lower than were seen with mucin from the intestinal tract. A. salmonicida bound to the mucins isolated from the intestinal tract to a greater extent than to the skin mucins. The mucins from the intestinal regions had higher levels of sialylation than the skin mucins. Desialylating intestinal mucins decreased A. salmonicida binding, whereas desialylation of skin mucins resulted in complete loss of binding. In line with this, A. salmonicida also bound better to mammalian mucins with high levels of sialylation, and N-acetylneuraminic acid appeared to be the sialic acid whose presence was imperative for binding. Thus, sialylated structures are important for A. salmonicida binding, suggesting a pivotal role for sialylation in mucosal defense. The marked differences in sialylation as well as A. salmonicida binding between the skin and intestinal tract suggest interorgan differences in the host-pathogen interaction and in the mucin defense against A. salmonicida.

De novo transcriptome sequencing of the snail Echinolittorina malaccana: identification of genes responsive to thermal stress and development of genetic markers for population studies

Echinolittorina snails inhabit the upper intertidal rocky shore and face strong selection pressures from thermal extremes and fluctuations. Revealing the molecular processes of adaptive significance is greatly obstructed by the scarcity of genomic resource for these taxa. Here, we reported the first comprehensive transcriptome dataset for the genus Echinolittorina. Using Illumina HiSeq 2000 platform, about 52 M and 54 M paired-end clean reads were, respectively, generated for the control and heat-stressed libraries. Totally, 115,211 unique transcript fragments (unigenes) were assembled, with an average length of 453 bp and a N50 size of 492 bp. Approximately one third of the unigenes could be annotated according to their homology matches against the Nr, Swiss-Prot, COG, or KEGG databases, and they were found to represent 23,098 non-redundant genes. Gene expression comparison revealed that 1,267 and 6,663 annotated genes were, respectively, up- and downregulated with at least twofold changes upon heat stress. Gene Ontology and KEGG pathway analyses indicated that there were overrepresented amount of genes enriched in a broad spectrum of biological processes and pathways, including those associated with cytoskeleton organization, developmental regulation, signaling transduction, infection, and cardiac function. In addition, a transcriptome-wide search for polymorphic loci yielded a total of 11,228 simple sequence repeats (SSRs) from 9,938 unigenes and 138,631 single nucleotide polymorphism (SNP) and insertion/deletion (INDEL) sites among 22,770 unigenes. The large number of transcript sequences acquired, the biological pathways identified, and the candidate microsatellite and SNP/INDEL loci discovered in the study will serve as valuable resources for further investigations of genetic differentiation and thermal adaptation among populations.

Extensive local gene duplication and functional divergence among paralogs in Atlantic salmon

Many organisms can generate alternative phenotypes from the same genome, enabling individuals to exploit diverse and variable environments. A prevailing hypothesis is that such adaptation has been favored by gene duplication events, which generate redundant genomic material that may evolve divergent functions. Vertebrate examples of recent whole-genome duplications are sparse although one example is the salmonids, which have undergone a whole-genome duplication event within the last 100 Myr. The life-cycle of the Atlantic salmon, Salmo salar, depends on the ability to produce alternating phenotypes from the same genome, to facilitate migration and maintain its anadromous life history. Here, we investigate the hypothesis that genome-wide and local gene duplication events have contributed to the salmonid adaptation. We used high-throughput sequencing to characterize the transcriptomes of three key organs involved in regulating migration in S. salar: Brain, pituitary, and olfactory epithelium. We identified over 10,000 undescribed S. salar sequences and designed an analytic workflow to distinguish between paralogs originating from local gene duplication events or from whole-genome duplication events. These data reveal that substantial local gene duplications took place shortly after the whole-genome duplication event. Many of the identified paralog pairs have either diverged in function or become noncoding. Future functional genomics studies will reveal to what extent this rich source of divergence in genetic sequence is likely to have facilitated the evolution of extreme phenotypic plasticity required for an anadromous life-cycle.

Transcriptome analysis reveals a rich gene set related to innate immunity in the Eastern oyster (Crassostrea virginica)

As a benthic filter-feeder of estuaries, the eastern oyster, Crassostrea virginica, faces tremendous exposure to microbial pathogens. How eastern oysters without adaptive immunity survive in pathogen-rich environments is of fundamental interest, but studies on its immune system are hindered by the lack of genomic resources. We sequenced the transcriptome of an adult oyster with short Illumina reads and assembled 66,229 contigs with a N50 length of 1,503 bp. The assembly covered 89.4 % of published ESTs and 97.9 % of mitochondrial genes demonstrating its quality. A set of 39,978 contigs and unigenes (>300 bp) were identified and annotated by searching public databases. Analysis of the gene set yielded a diverse set of 657 genes related to innate immunity, including many pertaining to pattern recognition, effectors, signal transduction, cytokines, and apoptosis. Gene families encoding C1q domain containing proteins, CTLD, IAPs, Ig_I-set, and TRAFs expanded in C. virginica and Crassostrea gigas. Many key genes of the apoptosis system including IAP, BAX, BAC-2, caspase, FADD, and TNFR were identified, suggesting C. virginica possess advanced apoptosis and apoptosis-regulating systems. Our results show that short Illumina reads can produce transcriptomes of highly polymorphic genomes with coverage and integrity comparable to that from longer 454 reads. The expansion and high diversity in gene families related to innate immunity, point to a complex defense system in the lophotrochozoan C. virginica, probably in adaptation to a pathogen-rich environment.

The first transcriptome and genetic linkage map for Asian arowana

Asian arowana or dragonfish (Scleropages formosus) is an important fish species due to its unusual breeding biology and high economic value in the ornamental fish markets. In the present study, we aimed to (i) create the first transcriptome by Roche 454 pyrosequencing of Asian arowana brain and gonad samples; (ii) identify differentially expressed genes between the two sexes and develop microsatellite (SSR) markers; and (iii) construct a first-generation SSR-based genetic linkage map. A total of over 1.3 million reads were obtained from the brain and gonad of adult Asian arowana individuals through pyrosequencing. These reads were assembled into 16,242 contigs that were further grouped into 13,639 isogroups. BLASTX annotation identified a total of 8316 unique proteins from this data set. Many genes with sexually dimorphic expression levels and some putatively involved in sex development were identified. A total of 316 EST-SSRs and 356 new genomic-SSRs were developed by screening through the current transcriptome data set and SSR-enriched genomic libraries. The first genetic linkage map of the species was constructed based on these markers. Linkage analysis allowed for mapping of 308 markers to 28 linkage groups (LGs), ranging in size from 14.9 to 160.6 cM. The potentially sex-associated gene sox9 was mapped to LG4 on the consensus linkage map. Pairwise putative conserved syntenies between the Asian arowana, zebrafish, and three-spined stickleback were also established. These resources will help the conservation of the species through better understanding of its phylogenetics, genomics and biology, and comparative genome analysis within the Osteoglossidae family.

Transcriptome of Atlantic cod (Gadus morhua L.) early embryos from farmed and wild broodstocks

Significant efforts have been made to elucidate factors affecting egg quality in fish. Recently, we have shown that eggs originating from wild broodstock (WB) of Atlantic cod (Gadus morhua L.) are of superior quality to those derived from farmed broodstock (FB), and this is associated with differences in the chemical composition of egg yolk. However, maternal transcripts, accumulated during oogenesis, have not been studied extensively in fish. The aim of the present study was to characterize putative maternal mRNA transcriptome in fertilized eggs of Atlantic cod and to compare transcript pools between WB and FB in order to investigate the relation between egg developmental potential and putative maternal mRNA deposits. We performed high-throughput 454 pyrosequencing. For each WB and FB group, five cDNA libraries were individually tagged and sequenced, resulting in 98,687 (WB) and 119,333 (FB) average reads per library. Sequencing reads were de novo assembled, annotated, and mapped. Out of 13,726 identified isotigs, 238 were differentially expressed between WB and FB, with 155 isotigs significantly upregulated in WB. The sequence reads were mapped to 11,340 different Atlantic cod transcripts and 158 sequences were differentially expressed between the 2 groups. Important transcripts involved in fructose metabolism, fatty acid metabolism, glycerophospholipid metabolism, and oxidative phosphorylation were differentially represented between the two broodstock groups, showing potential as biomarkers of egg quality in teleosts. Our findings contribute to the hypothesis that maternal mRNAs affect egg quality and, consequently, the early development of fish.

The mucosal immune system of fish: the evolution of tolerating commensals while fighting pathogens

The field of mucosal immunology research has grown fast over the past few years, and our understanding on how mucosal surfaces respond to complex antigenic cocktails is expanding tremendously. With the advent of new molecular sequencing techniques, it is easier to understand how the immune system of vertebrates is, to a great extent, orchestrated by the complex microbial communities that live in symbiosis with their hosts. The commensal microbiota is now seen as the "extended self" by many scientists. Similarly, fish immunologist are devoting important research efforts to the field of mucosal immunity and commensals. Recent breakthroughs on our understanding of mucosal immune responses in teleost fish open up the potential of teleosts as animal research models for the study of human mucosal diseases. Additionally, this new knowledge places immunologists in a better position to specifically target the fish mucosal immune system while rationally designing mucosal vaccines and other immunotherapies. In this review, an updated view on how teleost skin, gills and gut immune cells and molecules, function in response to pathogens and commensals is provided. Finally, some of the future avenues that the field of fish mucosal immunity may follow in the next years are highlighted.