Transcriptome profiling and digital gene expression analysis of Nile tilapia (Oreochromis niloticus) infected by Streptococcus agalactiae

Mechanisms of immune responses in Acanthopagrus latus to Streptococcus agalactiae infection revealed by transcriptomic analysis

Acanthopagrus latus (yellowfin seabream) is an economically important fish in the southeast coastal sea of China. Its slower growth rate makes it more prone to diseases in the cultivation period, leading to substantial economic losses. Epidemiological investigations have indicated that Streptococcus agalactiae is one of the most common Gram-positive pathogens, which have garnered increasing attention due to its high contagion and lethality rates in A. latus. In this work, an infection model of yellowfin seabream was established with an intraperitoneal injection of S. agalactiae. Clinical sign observations and various analyses, including histological examination, serum biochemical index assessment, immune-related enzyme level measurement, and transcriptome analysis of tissues (liver and intestine) with obvious clinical signs, were conducted for revealing the effects of S. agalactiae infection and immune response mechanisms in yellowfin seabream. The results indicate that evident clinical signs and multi-tissue damages with the notable changes in indices and significant increase in immune-related enzyme levels in the serum occurred in infected fish. RNA sequencing analysis identified 1130 differentially expressed genes (DEGs) in the liver and 1218 DEGs in the intestine, which were involved in multiple immune- and metabolism-related pathways via KEGG enrichment analysis. The transcriptomic results were further corroborated by quantitative real-time RT-PCR (qRT-PCR) tests of some specific immune-related genes. These findings provide new insights into the molecular immune mechanisms in yellowfin seabream following S. agalactiae infection and offer valuable reference data for disease prevention and molecular breeding (i.e., selective breeding through developing molecular markers of key genes).

Comparative transcriptome analysis reveals a serotype-specific immune response in Nile tilapia (Oreochromis niloticus) infected with Streptococcus agalactiae

Streptococcus agalactiae is a major causative agent of streptococcosis in Nile tilapia (Oreochromis niloticus) and understanding its etiology is important to ensure the sustainable development of global tilapia farming. Our research group recently observed contrasting disease patterns in animals infected with two different S. agalactiae serotypes (Ib and III). To better understand the basis for these divergent responses, we analyzed the brain transcriptome of Nile tilapia following bacterial exposure. Our findings revealed significant variation in the expression of genes involved in immune (e.g., CD209 antigen, granulin, C-X-C motif chemokine 10, prostacyclin synthase, and interleukins) and neuroendocrine (e.g., mmp13a, mmp9, brain aromatase, and pmch) pathways. The serotype Ib strain seems promptly recognized by the host, triggering a potent inflammatory response, whereas the serotype III strain elicited a less immediate response, resulting in more pronounced central nervous system (CNS) symptoms and behavioral effects. To the best of our knowledge, this is the first study to show serotype-specific immune responses to S. agalactiae in Nile tilapia. These findings are important for advancing disease management and control strategies in aquaculture. Identifying different immune reactions triggered by serotypes Ib and III may assist the development of more specific approaches for preventive measures, early detection, and effective treatment against streptococcosis.

Natural infection of hybrid sturgeon (Acipenser baerii♀×Acipenser schrenckii♂) with Nocardia seriolae and white sturgeon iridovirus: pathological and transcriptomic analyses

Introduction

In August 2023, hybrid sturgeons (Acipenser baerii♀×Acipenser schrenckii♂) cultured in Sichuan, China, showed infectious disease symptoms, including ulcers, liver and spleen nodules, and high mortality rates.

Methods

Pathogenic bacteria were isolated from the liver of diseased sturgeons and analyzed for their phenotypic and molecular traits. Furthermore, iridovirus-specific TaqMan real-time PCR (RT-PCR) analyses were conducted. The histopathological characteristics were analyzed using paraffin sectioning and transmission electron microscopy (TEM). Transcriptome sequencing was performed to elucidate the impact of pathogen exposure and immune response profiles in infected sturgeon.

Results

Pathogenic bacteria isolation and phylogenetic analyses of the 16S rRNA gene confirmed that the isolated bacteria clustered within the Nocardia seriolae group. The TaqMan RT-PCR assay was performed to detect the presence of white sturgeon iridovirus (WSIV), indicating a weakly positive signal. Histopathological examination revealed severe damage to various tissues, and a notable presence of bacteria was observed through acid-fast staining. Transmission electron microscopy analysis showed the presence of abundant bacteria and virus particles, indicating cellular invasion and subsequent damage. In summary, the disease in hybrid sturgeons was diagnosed as infection of N. seriolae and WSIV. To investigate the immune response of hybrid sturgeons to this infection, spleen transcriptomes were analyzed. Numerous immune-related genes and pathways, including the "Toll-like receptor", "B-cell receptor", and "T-cell receptor" signaling pathways, were altered in response to pathogenic threats. Significantly downregulated of key components of TCR and BCR signaling pathways, such as ZAP70, BTK, and CD79A, suggested a temporary inhibition of these pathways critical for cellular immunity post-infection. Gene set enrichment analysis revealed significant suppression of the apoptosis signaling pathway and activation of autophagy and mitophagy signaling pathways following infection. Specifically, in the death receptor-mediated apoptosis signaling pathway, downregulation of TNFα, TRAIL, CASP6, and CASP8 was observed, while several genes in the autophagy and mitophagy pathways showed upregulated expression post-infection.

Discussion

We report the initial occurrence of N. seriolae infection in cultured sturgeons. These findings could provide a theoretical basis for diagnosing and preventing this disease, as well as enhance the understanding of host-pathogen interactions in fish.

The role of HMGB2 in the immune response of Nile tilapia (Oreochromis niloticus) to streptococcal infection

High mobility group protein B2 (HMGB2) is an abundant chromatin-associated protein with pivotal roles in transcription, cell proliferation, differentiation, inflammation, and tumorigenesis. However, its immune function in Nile tilapia (Oreochromis niloticus) remains unclear. In this study, we identified a homologue of HMGB2 from Nile tilapia (On-HMGB2) and investigated its functions in the immune response against streptococcus infection. The open reading frame (ORF) of On-HMGB2 spans 642 bp, encoding 213 amino acids, and contains two conserved HMG domains. On-HMGB2 shares over 80 % homology with other fish species and 74%-76 % homology with mammals. On-HMGB2 was widely distributed in various tissues, with its highest transcript levels in the liver and the lowest in the intestine. Knockdown of On-HMGB2 promoted the inflammatory response in Nile tilapia, increased the bacterial load in the tissues, and led to elevated mortality in Nile tilapia following Streptococcus agalactiae infection. Taken together, On-HMGB2 significantly influences the immune system of Nile tilapia in response to streptococcus infection.

Comparative transcriptomes reveal different tolerance mechanisms to Streptococcus agalactiae in hybrid tilapia, nile tilapia, and blue tilapia

Tilapia is one of the most economically important freshwater fish farmed in China. Streptococcosis outbreaks have been extensively documented in farmed tilapia species. Hybrid tilapia (Oreochromis niloticus ♀ × O. aureus ♂) exhibit greater disease resistance than Nile tilapia (O. niloticus) and blue tilapia (O. aureus). However, the molecular mechanism underlying the enhanced tolerance of hybrid tilapia is still poorly understood. In this study, comparative transcriptome analysis was performed to reveal the different tolerance mechanisms to Streptococcus agalactiae in the three tilapia lines. In total, 1982, 2355, and 2076 differentially expressed genes were identified at 48 h post-infection in hybrid tilapia, Nile tilapia, and blue tilapia, respectively. Functional enrichment analysis indicated that numerous metabolic and immune-related pathways were activated in all three tilapia lines. The differential expression of specific genes associated with phagosome, focal adhesion, cytokine-cytokine receptor interaction, and toll-like receptor signaling pathways contributed to the resistance of hybrid tilapia. Notably, immune response genes in hybrid tilapia, such as P38, TLR5, CXCR3, CXCL12, PSTPIP1, and TFR, were generally suppressed under normal conditions but selectively induced following pathogen challenge. These results expand our knowledge of the molecular mechanisms underlying S. agalactiae tolerance in hybrid tilapia and provide valuable insights for tilapia breeding programs.

Transcriptome analysis reveals the neuroactive receptor genes response to Streptococcus agalactiae infection in tilapia, Oreochromis niloticus

The detailed crosstalk between the neuroendocrine and immune systems in Oreochromis niloticus, an economically important fish, in response to pathogenic infections, remains unclear. This study revealed the head kidney transcriptional profiles of O. niloticus upon infections with Streptococcus agalactiae, a prevalent pathogen known to cause severe meningitis. Twelve cDNA libraries of O. niloticus head kidney, representing four treatment time points (0, 6, 24, and 48 h), were constructed and a total of 2,528 differentially expressed genes were identified based on pairwise comparisons. KEGG pathway analysis revealed a significant enrichment of the 'neuroactive ligand-receptor interaction' pathway (ko04080), with 13 genes exhibiting differential expression during S. agalactiae infection. Among these, six neuroactive receptor genes (lepr, nr3c1, ptger4, thrb, tspo, and β2-ar) were selected, cloned, and characterized. Although these genes are ubiquitously expressed, and in head kidney leukocytes, their expression was mainly observed in T cells, Mo/Mφ, and NCCs, which are characterized by antimicrobial responses. Furthermore, we examined the response patterns of these six neuroactive receptor genes to gram-positive (S. agalactiae) and gram-negative (Aeromonas hydrophila) bacteria in four different tissues. Notably, lepr, ptger4, tspo, and β2-ar were upregulated in all selected tissues in response to S. agalactiae and A. hydrophila infections. However, nr3c1 and thrb were downregulated in response to S. agalactiae infection in the head kidney and spleen, whereas nr3c1 was upregulated, and thrb was unresponsive to A. hydrophila infection. Our findings provide a theoretical foundation for understanding new links between the neuroendocrine and immune systems during bacterial infection in teleost fish.

Genome-Wide Identification of Trachinotus ovatus Antimicrobial Peptides and Their Immune Response against Two Pathogen Challenges

Golden pompano, Trachinotus ovatus, as a highly nutritious commercially valuable marine fish, has become one of the preferred species for many fish farmers due to its rapid growth, wide adaptability, and ease of feeding and management. However, with the expansion of aquaculture scale, bacterial and parasitic diseases have also become major threats to the golden pompano industry. This study, based on comparative genomics, shows the possibility of preferential evolution of freshwater fish over marine fish by analyzing the phylogenetic relationships and divergence times of 14 marine fish and freshwater fish. Furthermore, we identified antimicrobial peptide genes from 14 species at the genomic level and found that the number of putative antimicrobial peptides may be related to species evolution. Subsequently, we classified the 341 identified AMPs from golden pompano into 38 categories based on the classification provided by the APD3. Among them, TCP represented the highest proportion, accounting for 23.2% of the total, followed by scolopendin, lectin, chemokine, BPTI, and histone-derived peptides. At the same time, the distribution of AMPs in chromosomes varied with type, and covariance analysis showed the frequency of its repeat events. Enrichment analysis and PPI indicated that AMP was mainly concentrated in pathways associated with disease immunity. In addition, our transcriptomic data measured the expression of putative AMPs of golden pompano in 12 normal tissues, as well as in the liver, spleen, and kidney infected with Streptococcus agalactiae and skin infected with Cryptocaryon irritans. As the infection with S. agalactiae and C. irritans progressed, we observed tissue specificity in the number and types of responsive AMPs. Positive selection of AMP genes may participate in the immune response through the MAPK signaling pathway. The genome-wide identification of antimicrobial peptides in the golden pompano provided a complete database of potential AMPs that can contribute to further understanding the immune mechanisms in pathogens. AMPs were expected to replace traditional antibiotics and be developed into targeted drugs against specific bacterial and parasitic pathogens for more precise and effective treatment to improve aquaculture production.

Applying genetic technologies to combat infectious diseases in aquaculture

Disease and parasitism cause major welfare, environmental and economic concerns for global aquaculture. In this review, we examine the status and potential of technologies that exploit genetic variation in host resistance to tackle this problem. We argue that there is an urgent need to improve understanding of the genetic mechanisms involved, leading to the development of tools that can be applied to boost host resistance and reduce the disease burden. We draw on two pressing global disease problems as case studies-sea lice infestations in salmonids and white spot syndrome in shrimp. We review how the latest genetic technologies can be capitalised upon to determine the mechanisms underlying inter- and intra-species variation in pathogen/parasite resistance, and how the derived knowledge could be applied to boost disease resistance using selective breeding, gene editing and/or with targeted feed treatments and vaccines. Gene editing brings novel opportunities, but also implementation and dissemination challenges, and necessitates new protocols to integrate the technology into aquaculture breeding programmes. There is also an ongoing need to minimise risks of disease agents evolving to overcome genetic improvements to host resistance, and insights from epidemiological and evolutionary models of pathogen infestation in wild and cultured host populations are explored. Ethical issues around the different approaches for achieving genetic resistance are discussed. Application of genetic technologies and approaches has potential to improve fundamental knowledge of mechanisms affecting genetic resistance and provide effective pathways for implementation that could lead to more resistant aquaculture stocks, transforming global aquaculture.

Comparative Transcriptomic Immune Responses of Mullet (Mugil cephalus) Infected by Planktonic and Biofilm Lactococcus Garvieae

Lactococcus garvieae is an important pathogen of fish, associated with high rates of mortality and infection recurrence in summer or stressful conditions. Chronic infection and disease recurrence have also been reported to be associated with biofilms. However, the impact of biofilm and planktonic bacterial infection on fish immune responses remains unclear. In this study, de novo sequencing was used to compare differences of the spleen transcriptome in planktonic- and biofilm-infected mullets. Among the 181,024 unigenes obtained, 3,392 unigenes were associated with immune response genes. Comparative analysis of the gene expression between infection with the L. garvieae planktonic type and biofilm type identified a total of 3,120 and 3,489 differentially expressed genes in response to planktonic and biofilm infection, respectively, of which 1,366 and 1,458 genes were upregulated, and 1,754 and 1,458 genes were downregulated, respectively. Gene ontology enrichment analysis of immune genes identified genes involved in the complement system, toll-like receptor signaling, and antigen processing, which were further verified by qPCR. Additionally, genes encoding TLR2, IL-1β, TNF-α, C7, and MHC class II peptides were downregulated in response to biofilm infection. Importantly, the results show that biofilm infection induces a different immune pathway response compared with planktonic bacterial infection and, furthermore, illustrates that the prevention of biofilm formation may be a necessary and new strategy for controlling bacterial infection in aquaculture.

Molecular Identification of Nocardia seriolae and Comparative Analysis of Spleen Transcriptomes of Hybrid Snakehead (Channa maculata Female × Channa argus Male) With Nocardiosis Disease

Hybrid snakehead (Channa maculata female × Channa argus male) is a new freshwater aquaculture fish species in southern China. During intensive aquaculture, hybrid snakeheads are often infected by Nocardia seriolae. In this study, hybrid snakehead infected suspiciously by N. seriolae in an artificial breeding pond were examined. Diseased hybrid snakeheads swam slowly without food intake, and the clinical symptoms included skin wound, anal swelling and ascites, and white granulomatous in liver, spleen, and kidney of fish. Through bacterial isolation, 16S rDNA sequencing, fluorescence in situ hybridization (FISH) and artificial infection experiment, the pathogen was identified as N. seriolae. Furthermore, the spleen samples from diseased and healthy male hybrid snakeheads in the same pond were used for RNA-Seq analysis. A total of 3,512 unique transcripts (unigenes) were identified as differentially expressed genes (DEGs), and 1,886 of them were up-regulated in diseased fish. The expression patterns of 20 DEGs were verified by quantitative polymerase chain reaction (qPCR). Several immune-related pathways and many immune-related genes were identified. qPCR results showed that the expression patterns of immune-related genes in the liver and kidney of diseased fish were comparable to that in the spleen. This study provides deep-sequencing data of hybrid snakehead spleen and will help understand the immune response of hybrid snakehead to N. seriolae. It is also helpful for the biomarker screening of fish-borne Nocardia spp. and the breeding of nocardiosis-resistant fish species.

Comparative Transcriptomic Analysis Reveals the Regulated Expression Profiles in Oreochromis niloticus in Response to Coinfection of Streptococcus agalactiae and Streptococcus iniae

Tilapia (Oreochromis sp.) is one of the important economical fishes in the world. Streptococcosis is commonly found in tilapia, causing severe and devastating effects in tilapia cultures. Streptococcus agalactiae and Streptococcus iniae are the predominant pathogens causing tilapia streptococcosis. To understand the molecular mechanisms underlying differential streptococcal infection patterns, Nile tilapias (Oreochromis niloticus) were infected by 1 × 107 CFU/mL S. agalactiae, 1 × 107 CFU/mL S. iniae, and 1 × 107 CFU/mL S. agalactiae and S. iniae (1:1), respectively, and transcriptome analysis was conducted to the intestine samples of Nile tilapia (Oreochromis niloticus) at 6, 12, 24 h, and 7 days post-infection. A total of 6,185 genes that differentially expressed among groups were identified. Eight differentially expressed genes (DEGs) including E3 ubiquitin-protein ligase TRIM39-like, C-X-C motif chemokine 10-like(CXCL 10), C-C motif chemokine 19-like, interleukin-1 beta-like, IgM heavy chain VH region, partial, IgG Fc-binding protein, proteasome subunit beta type-8 (PSMB8), and ATP synthase F(0) complex subunit B1, mitochondrial that involved in the immune system were selected, and their expression levels in the coinfection group were significantly higher than those in either of the single infection groups. These genes were associated with four different KEGG pathways. Additionally, the differential expression of eight DEGs was validated by using the RT-qPCR approach, and their immunological importance was discussed. The results provided insights into the responses of tilapia against S. agalactiae and S. iniae at the transcriptome level, promoting our better understanding of immune responses for aquatic animal against Streptococcus.

Gene expression profiling of DNA methyltransferase genes in Siniperca chuatsi based on transcriptome sequencing

The mandarin fish (Siniperca chuatsi) DNA methyltransferase gene 1 (dnmt1) was highly expressed in the mesonephros, head kidney and gonad, whereas dnmt2 was expressed in most tissues. dnmt3a was highly expressed in the brain and spleen, but dnmt3b was mainly expressed in the brain and head kidney. The genes dnmt1 and dnmt2 were highly expressed in the early stages of embryonic development, and dnmt3a and dnmt3b were expressed later. These genes also showed certain changes after artificial diet acclimation, salinity adaptation and immune stress.

Granulomatous bacterial diseases in fish: An overview of the host's immune response

Bacterial diseases represent the main impediment to the development of fish aquaculture. Granulomatous diseases caused by bacteria lead to fish culture losses by high mortality rates and slow growth. Bacteria belonging to genera Streptococcus spp., Mycobacterium sp., Nocardia sp., Francisella sp., and Staphylococcus sp. have been implicated in the development of granulomatous processes. The granuloma formation and the fish's immune response continue to be the subject of scientific research. In fish, the first defense line is constituted by non-specific humoral factors through growth-inhibiting substances such as transferrin and antiproteases, or lytic effectors as lysozyme and antimicrobial peptides, and linking with non-specific phagocyte responses. If the first line is breached, fish produce antibody constituents for a specific humoral defense inhibiting bacterial adherence, as well as the mobilization of non-phagocytic host cells and counteracting toxins from bacteria. However, bacteria causing granulomatous diseases can be persistent microorganisms, difficult to eliminate that can cause chronic diseases, even using some immune system components to survive. Understanding the infectious process leading to granulomatosis and how the host's immune system responds against granulomatous diseases is crucial to know more about fish immunology and develop strategies to overcome granulomatous diseases.

LECT2 Protects Nile Tilapia (Oreochromis niloticus) Against Streptococcus agalatiae Infection

Leukocyte cell-derived chemotaxin 2 (LECT2) is a multifunctional cytokine that especially plays an important role in innate immune. However, the roles of LECT2 in the immune response of the economically important fish Nile tilapia (Oreochromis niloticus) against bacterial infection remains unclear. In this study, a lect2 gene from Nile tilapia (On-lect2) was identified, and its roles in the fish’s immune response against bacterial infection were determined and characterised. On-lect2 contains an open reading frame of 456 bp that encodes a peptide of 151 amino acids, as well as the conservative peptidase M23 domain. On-LECT2 is 62%–84% identical to other fish species and about 50% identical to mammals. The highest transcriptional level of On-lect2 was detected in the liver, whereas the lowest levels were detected in the other tissues. Moreover, the On-LECT2 protein is located mainly in the brain and head kidney. The transcriptional levels of On-lect2 substantially increased in the head kidney, brain, liver and spleen after Streptococcus agalactiae infection. Knockdown On-lect2 led to higher mortality due to liver necrosis or haemorrhage and splenomegaly. In vitro analysis indicated that the recombinant protein of On-LECT2 improved phagocytic activity of head kidney-derived macrophages. In vivo challenge experiments revealed several functions of On-LECT2 in the immune response of Nile tilapia against bacterial infection, including promotion of inflammation, reduction of tissue damages and improvement of survival rate.

Single molecule real-time sequencing and RNA-seq unravel the role of long non-coding and circular RNA in the regulatory network during Nile tilapia (Oreochromis niloticus) infection with Streptococcus agalactiae

BACKGROUND

The tilapia aquaculture industry is facing heavy economic losses due to Streptococcus agalactiae (S. agalactiae) infections. While progress has been made in past years, the lack of a high-quality tilapia genome and transcript annotations makes systematic and comprehensive exploration for a non-coding RNA regulatory network associated with the infection process unfeasible, and it stunts further research focused on disease defense and treatment. Herein, single molecular real time sequencing (SMRT-Seq) and RNA-seq data were utilized to generate a high-quality transcript annotation. In addition, Changes in mRNA and non-coding RNA expression were also analyzed during a S. agalactiae infection in tilapia.

FINDINGS

In total, 16.79 Gb of clean data were obtained by sequencing on six SMRT cells, with 712,294 inserts (326,645 full-length non-chimeric reads and 354,188 non-full-length reads). A total of 197,952 consensus transcripts were obtained. Additionally, 55,857 transcript sequences were acquired, with 12,297 previously annotated and 43,560 newly identified transcripts. To further examine the immune response in Oreochromis niloticus following a S. agalactiae infection, a total of 470.62 Gb of clean data was generated by sequencing a library containing 18 S. agalactiae infected tilapia samples. Of the identified genes, 9911 were newly exploited, of which 7102 were functional annotated. Furthermore, 7874 mRNAs, 1281 long non-coding RNAs (out of 21,860 long non-coding RNAs), and 61 circular RNAs (out of 1026 circular RNAs) were found to be differentially expressed during infection, with the 1026 circRNAs not previously identified in tilapia. Moreover, k-means clustering and WGCNA analyses revealed that the immune response of tilapia to a S. agalactiae infection can be divided into three stages: cytokines driven rapid immune response, energy metabolism promotion, and the production of lysosomes and phagosomes. During this response, the head kidney and spleen have synergistic effects, while maintaining independent characteristics. Finally, lncRNA-mRNA (trans and cis), lncRNA-miRNA-mRNA, and circRNA-miRNA-mRNA regulatory networks were constructed and revealed that non-coding RNA is involved in the regulation of immune-related genes.

CONCLUSIONS

This study generated a greatly-improved transcript annotation for tilapia using long-read PacBio sequencing technology, and revealed the presence of a regulatory network comprised of non-coding RNAs in Nile tilapia infected with S. agalactiae.

Genomics to accelerate genetic improvement in tilapia

Selective breeding of tilapia populations started in the early 1990s and over the past three decades tilapia has become one of the most important farmed freshwater species, being produced in more than 125 countries around the globe. Although genome assemblies have been available since 2011, most of the tilapia industry still depends on classical selection techniques using mass spawning or pedigree information to select for growth traits with reported genetic gains of up to 20% per generation. The involvement of international breeding companies and research institutions has resulted in the rapid development and application of genomic resources in the last few years. GWAS and genomic selection are expected to contribute to uncovering the genetic variants involved in economically relevant traits and increasing the genetic gain in selective breeding programs, respectively. Developments over the next few years will probably focus on achieving a deep understanding of genetic architecture of complex traits, as well as accelerating genetic progress in the selection for growth-, quality- and robustness-related traits. Novel phenotyping technologies (i.e. phenomics), lower-cost whole-genome sequencing approaches, functional genomics and gene editing tools will be crucial in future developments for the improvement of tilapia aquaculture.

An effective live attenuated vaccine against Streptococcus agalactiae infection in farmed Nile tilapia (Oreochromis niloticus)

Streptococcus agalactiae is an important pathogen associated with various aquatic animals, especially tilapia. Streptococcosis has greatly limited the healthy development of tilapia aquaculture in recent times. The development of novel effective vaccines is important for the prevention and control of streptococcosis in fish. We previously constructed a non-encapsulated S. agalactiae strain △cps by the in-frame deletion method. Here, we evaluated whether this mutant △cps is safe for tilapia and suitable for protection against streptococcosis. We observed that the △cps strain was non-pathogenic to tilapia, and there was no reversion of virulence when it was passaged in tilapia. Moreover, the △cps strain survived for at least 11 d in the main immune organs of tilapia. The tilapia vaccinated via intraperitoneal (IP) injection with △cps strain induced a high antibody titer, and the IgM antibody levels were significantly higher in the vaccinated group than in the control group. The vaccination protected tilapia against the S. agalactiae challenge with a relative percent survival of 90.47%. In addition, tilapia immunized with the △cps strain showed significantly higher expression level of IFN-γ, IL-1β, MyD88, IgM, and MHC-Iα in the head kidney than those in the control during the entire observation period. The expression of MHC-IIβ was inhibited during 1-7 d of immunization. These results revealed that the △cps strain is able to induce humoral and cell-mediated immune response in tilapia. Therefore, the strain △cps has a broad application prospect as a target for attenuation in vaccine development.

Inorganic nitrite increases the susceptibility of tilapia (Oreochromis niloticus) leucocytes to Streptococcus agalactiae

Deteriorating water quality, especially from high concentrations of nitrite, is currently largely blamed for disease outbreaks in farmed tilapia (Oreochromis niloticus). In this study, the underlying mechanism of nitrite on the susceptibility of tilapia leucocytes to Streptococcus agalactiae (S. agalactiae) was studied. We found that a high dose of heat-killed S. agalactiae decreased tilapia leucocytes cell viability, whereas nitrite decreased the cell viability of leucocytes exposed to a low dose of bacteria. Bacterial challenge increased the production of nitric oxide (NO), whereas nitrite and bacteria coexposure caused higher NO production than nitrite or bacterial exposure alone. Cell viability increased after elimination of NO, and negative correlations existed between cell viability and the NO content, suggesting that nitrite increased the susceptibility of the leucocytes against S. agalactiae was NO-dependent. For a more comprehensive understanding of the mechanism of nitrite affecting disease resistance in tilapia leucocytes, an RNA-Seq-based transcriptome was generated. The results showed that 6173 transcripts were differently expressed, and the differentially expressed transcripts (DETs) of the bacterial group, nitrite group and bacteria-nitrite co-treatment group compared to the control group were selected for GO and KEGG analyses. The DETs in the bacterial group and bacteria-nitrite cotreatment group were highly involved with the membrane component, signal transduction, and immune responses. KEGG analysis showed that the protein processing in the endoplasmic reticulum and the AMPK signaling pathway, which are related to autophagy, were significantly enriched in the cotreatment group but not in bacterial group. In addition, the mRNA expression of ten DETs and several autophagy and apoptosis related genes validated by q-PCR showed the high reliability of the RNA-seq. Taken together, the results of this study suggest that nitrite may increase the susceptibility of tilapia leucocytes to S. agalactiae by generating excess NO to affect the autophagy and apoptosis process.

Digital gene expression analysis in the liver of ScpB-vaccinated and Streptococcus agalactiae-challenged Nile tilapia

In recent years, streptococcal diseases have severely threatened the development of tilapia aquaculture, but effective prevention and control methods have not yet been established. To understand the immune responses of vaccinated Nile tilapia (Oreochromis niloticus), digital gene expression (DGE) technology was applied in this study to detect the gene expression profile of the Nile tilapia (O. niloticus) liver in response to ScpB (Streptococcal C5a peptidase from group B Streptococcus, ScpB) vaccination and a Streptococcus agalactiae-challenge. The control and the ScpB-vaccinated Nile tilapia yielded a total of 25,788,734 and 27,088,598 clean reads, respectively. A total of 1234 significant differentially expressed unigenes were detected (P < 0.05), of which 236 were significantly up-regulated, and 269 were significantly down-regulated (P < 0.05, |fold|>2, FDR<0.05). Of the differentially expressed gene, the identified genes which were enriched using databases of GO and KEGG could be categorized into a total of 67 functional groups and were mapped to 153 signaling pathways including 15 immune-related pathways. The differentially expressed genes (TLR1, TLR2, TLR3, TLR5, TLR9, MyD88, C3, IL-1β, IL-10) were detected in the expression profiles, and this was subsequently verified via quantitative real-time PCR (qPCR). The results of this study can serve as a basis for future research not only on the molecular mechanism of S. agalactiae invasion, but also on the anti-S. agalactiae mechanism in targeted tissues of Nile tilapia.

Transcriptome analysis of juvenile genetically improved farmed tilapia (Oreochromis niloticus) livers by dietary resveratrol supplementation

Here we used RNA-Seq to explore the transcriptomic response and specific involvement of hepatic mRNA of juvenile Oreochromis niloticus (GIFT) as a result of dietary resveratrol supplementation (0.05 g/kg RES). More than 24,513,018 clean reads were reference genome guided assembly into 23,417 unigenes. 12,596 unigenes (29.64%) were annotated to GO database. There were 5, 179 and 1526 genes significantly differentially expressed genes at 15, 30 and 45 d respectively, and 8 KEGG pathways were enriched associated with this immune response. Hyperemia and compressed hepatic sinusoid, fibrosis of liver cell and abnormal hepatic epidermal cell revealed by H&E and SEM analysis respectively. Genes related with cytokine production (il12rb2, scfr), immune system (ig8l, hlfl, cd226, prf1l), autophagy regulation (atg4b), foxo signaling (ccnb2), steroid hormone biosynthesis (cyp3a40), fatty acid metabolism (scd1), metabolism (cacna1b) have been significantly decreased, while genes associated with such pathways above (leap-2, prdx4, mb, homer1, mif, sat1, cytbc1_8) and the pathway of protein processing in endoplasmic reticulum (cne1, tram1) have been significantly increased. These findings suggested RES activated some immune and biological process-related genes to enhance GIFT's innate immunity. It also suggested high concentration addition or long-time administration may bring negative effect in tilapia liver.

The important role of phagocytosis and interleukins for nile tilapia (Oreochromis niloticus) to defense infection of Aeromonas hydrophila based on transcriptome analysis

Tilapia is an important economic fish worldwide. It is vital to understand the mechanism of immune response for the prevention and treatment the infection of Aeromonas hydrophila. Based on high-throughput sequencing of Illumina HiSeq™, we found differentially expressed genes in the immune-related pathway were classified into phagosome, cytokine-cytokine receptor interaction and toll-like receptor signaling pathway. Gene Ontology terms were divided into three categories of transporting function, DNA replication activity and energy supply activity. The first one was related to phagocytosis and the process or transporting of antigen driven by tubulins; the second one was to differentiation and proliferation of lymphocyte activated by cytokines; and the former two both needed energy provided by the third one. According to colchicine assay, cross-immune assay, ELISA of interleukins and classical phagocytosis assay, phagocytosis and interleukins were verified to be most important to defense the infection of A. hydrophila.

Molecular cloning, characterization and expression analysis of MHCI and chemokines CXCR3 and CXCR4 gene from freshwater carp, Catla catla

The immune system with large number of molecules protects the host against a plethora of continuously evolving microbes. Major histocompatibility complex (MHC) molecules serve as cardinal elements of the adaptive immune system responsible for the activation of the adaptive immunity in the host. The present study reports MHCI molecule in freshwater carp, Catla catla, and its differential expression in immunologically relevant tissues post-infection with Gram-negative and Gram-positive bacteria. The MHCI sequence of C. catla had 502 bp nucleotides encoding putative 146 amino acids. The phylogenetic analysis exhibited its evolutionary conservation within the Cyprinidae family and formed a different clade with the higher vertebrates. Simultaneously, CXCR3 and CXCR4 chemokines were cloned and characterized for their expression in infected tissues. Analysis of immunologically relevant tissues of the infected fish exhibited an increase of MHCI gene expression and the down-regulation of CXCR3 and CXCR4 chemokines, indicating a tricky interaction between the innate and adaptive immune system. It was found that intestine, skin and spleen played a crucial role in the contribution of the defense activity which instigated the self-immunity. These immune activities can provide useful information to understand the interaction of self and non-self- immune system in freshwater fish, Catla catla.

Identification and expression profiling analysis of microRNAs in Nile tilapia (Oreochromis niloticus) in response to Streptococcus agalactiae infection

MicroRNAs (miRNAs) play vital regulatory roles in various biological processes, including in immune responses. Nile tilapia (Oreochromis niloticus) is an important commercial fish species in China. To identify immune-related miRNAs of O. niloticus, 4 libraries from liver during S. agalactiae infection (0 h, 5 h, 50 h, and 7 d) were sequenced by high-throughput sequencing technology in tilapia. We obtained 10,703,531, 11,507,163, 11,180,179 and 13,408,414 clean reads per library, respectively. In our results, a total of 482 miRNAs were identified through bioinformatic analysis, including 220 conserved miRNAs and 262 putative novel miRNAs. Moreover, 21 (4.36%), 50 (10.37%), and 46 (9.54%) miRNAs were significantly differentially expressed at 5 h, 50 h and 7 d, respectively. In addition, 6939 target genes regulated by these differentially expressed miRNAs were predicted, and their functional annotations were predicted by Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, which revealed that a majority of differentially expressed miRNAs were involved in apoptotic process, metabolic process, and immune responses. Finally, Real-time quantitative PCR experiments were performed for 7 miRNAs by stem-loop RT-PCR, and a general agreement was confirmed between the sequencing and RT-qPCR data. To our understanding, this is the first report of comprehensive identification of O. niloticus miRNAs being differentially regulated in liver related to S. agalactiae infection. This work provides an opportunity for further understanding of the molecular mechanisms of miRNA regulation in O. niloticus host-pathogen interactions, and genetic resources for molecular assistant selection for disease resistant breeding program.

Transcriptome profiles in the spleen of African catfish (Clarias gariepinus) challenged with Aeromonas veronii

The African catfish, Clarias gariepinus, an important cultured freshwater species in many countries, possess the characteristic of high disease resistance. However, little genomic information for this character of the fish is available up to now. To address the shortfall and to better understand C. gariepinus immune response to pathogen infection at molecular level, C. gariepinus were challenged with potent A. veronii and the high-throughput RNA sequencing (RNA-seq) technology were employed to produce transcriptomes from spleen. In total, an average of 46,073,372 clean reads obtained were de novo assembled into 156,955 unigenes with an average length of 1082 bp. All of unigenes were annotated to seven public databases. Three comparisons were separately conducted between the infected groups at 3 h, 24 h, 48 h post-challenge and control group. A total of 2482 differentially expressed unigenes (DEGs) were identified. Among these, 114 immune-related DEGs were captured, including 88, 42, and 31 genes at 3 h, 24 h and 48 h after infection respectively, for analysis of expression pattern and enrichment. The 114 DEGs displayed four expression patterns by cluster analysis and they were significantly enriched in 38 pathways (q < 0.01) related to the immune or disease, five of which were NF-kappa B, TNF, NLR, TLR and RLR pathways. Finally, the expression levels of twelve selected immune-related DEGs involved in above five pathways were scrutinized. Seven of which were up-regulated at 3 h after infection, afterward, their expression dropped to control level. In summary, this study provides valuable transcriptome resource for understanding the defense mechanisms of C. gariepinus in resistance to pathogens from the gene expression viewpoint, which also open up the possibility to study the immune complexity and to better comprehend the interrelationships between some immune pathways in C. gariepinus.

Transcriptome analysis of immune response in fat greenling (Hexagrammos otakii) against Vibrio harveyi infection

Fat greenling (Hexagrammos otakii) is an important aquaculture fish species in northern China. Unfortunately, Vibrio infections have caused considerable losses to the fat greenling aquaculture industry. However, the study on immune response of fat greenling against Vibrio species has not been reported yet. In this paper, the immune response of fat greenling against V. harveyi at gene expression level was studied by transcriptome analysis. A total of 189753 high-quality unigenes with a N50 length of 672bp were obtained by transcriptome profiling, which provided abundant data for the future study of fat greenling. Comparative analysis showed that 5425 differentially expressed genes (DEGs) were identified on day 3 post-infection (3dpi), containing 1837 up-regulated and 3588 down-regulated genes. Further annotation and analysis revealed that the DEGs were enriched in complement and coagulation cascades, ribosome, oxidative phosphorylation, glycine, serine and threonine metabolism and peroxisome proliferator-activated receptor (PPAR) signaling pathway. These pathways were mainly associated with phagocytosis and pathogen clearance, rarely involved in bacteria adhesion and pathogen identification, which suggested that the host might begin to clear and kill the invading bacteria on 3dpi. The research might provide a valuable resource to further study immune response and suggest strategies against V. harveyi infection in fat greenling.

Fipronil-induced toxic effects in zebrafish (Danio rerio) larvae by using digital gene expression profiling

Fipronil residue has caused widespread concern around the world, especially after the recent "toxic eggs" event in seven European countries. To evaluate the effects of fipronil on vertebrates, zebrafish larvae were used as an animal model to examine the lethal effect, developmental phenotypes at high doses, and possible mechanisms of toxicity by employing digital gene expression (DGE) profiling at environmentally relevant doses. The results of acute toxicity test indicated that treatment with fipronil from 75 h post-fertilization (hpf) led to the death of larvae with a 96-h LC50 value of 459 μg/L, as well as abnormal development including bent spine and shortened body length. Besides, we obtained high-quality-sequencing DGE profilings at fipronil concentrations of 0.5, 5, and 50 μg/L, respectively. The results revealed that 44 differentially expressed genes, 10 GO terms, and 3 KEGG pathways were overlapped among the three concentrations. MIDN, one of the 44 differentially expressed genes, showed dose-dependent responses at the transcriptional level, indicating that it was possibly a potential biomarker to reflect fipronil toxicity in zebrafish. Furthermore, we presumed that the changing transcriptional level of AP-1 family was possibly a reason for bent spine and shortened body length in larvae exposed to fipronil. Concurrently, altered abundance of transcripts of the ELOVL family in a key step of fatty acid elongation could possibly lead to the accumulation of long-chain fatty acids. Collectively, our results suggested that exposure to fipronil caused lethal and developmental toxicity in zebrafish larvae, and demonstrated the need for a comprehensive understanding of the potential mechanisms of fipronil toxicity due to fipronil's frequent presence in the environment and its potential threat to human health.

Interaction of Group B Streptococcus sialylated capsular polysaccharides with host Siglec-like molecules dampens the inflammatory response in tilapia

Group B Streptococcus (GBS, S. agalactiae) infection in tilapia (Oreochromis niloticus) causes widespread death of this species and is a significant issue for the aquaculture industry. The major virulence factor for GBS is its sialylated capsular polysaccharides (CPs). These CPs interact with sialic acid-binding immunoglobulin-like lectins (Siglecs) on the host immune cells to regulate the downstream inflammatory response and evade detection. Previously, we cloned multiple Siglec-like molecules from an O. niloticus cDNA library, all of which were shown to interact with the sialylated CPs of GBS. In the present study, we investigated the effects of GBS infection on the expression of pro- and anti-inflammatory cytokines in O. niloticus as well as OnSiglec-like-transfected macrophage cells. Eukaryotic expression vectors containing full-length OnSiglec-1-like, -4b-like, -14-like were constructed and used to transfect RAW264 macrophages in vitro as well as live tilapia in vivo prior to GBS infection. The expression of the anti-inflammatory cytokine interleukin (IL)-10 and the pro-inflammatory cytokines tumor necrosis factor (TNF)-α, IL-6, and interferon (INF)-β were then analyzed by qPCR. Our results indicate that as infection progressed, IL-10 expression was significantly upregulated, while that of TNF-α and IL-6 were significantly downregulated in the OnSiglec-like-transfected cells. INF-β expression was also downregulated in cells transfected with OnSiglec-1-like and -4b-like, but was not significantly effected in OnSiglec-14-like-transfected cells. Notably, the magnitude of these cytokine expression changes was greatly decreased when a ΔneuA GBS mutant was used to infect the OnSiglec-1-like-transfected cells. In GBS-infected tilapia, IL-10 expression was significantly upregulated in all tissues, whereas INF-β expression in the spleen, kidney, and gills was significantly downregulated at 12 hpi. While the expression of TNF-α was slightly upregulated, this change was not significant. In GBS ΔneuA mutant-infected O. niloticus, IL-10 expression in all of the tissues was significantly lower than that observed for the wild-type GBS group, while TNF-α expression was higher in the mutant infected group. There was no significant difference in INF-β expression between the two groups. Taken together, sialylated CPs on GBS appear to interact with host OnSiglec-like molecules to transmit negative regulatory signals via enhanced anti-inflammatory cytokine IL-10 production and reduced pro-inflammatory cytokine production, ultimately leading to dampening of the host immune response. The results of this study further elucidate the molecular mechanism underlying GBS infection in tilapia and also provide candidate drug target molecules.

Immune-Related Functional Differential Gene Expression in Koi Carp (Cyprinus carpio) after Challenge with Aeromonas sobria

In order to understand the molecular basis underlying the host immune response of koi carp (Cyprinus carpio), Illumina HiSeq^TM 2000 is used to analyze the muscle and spleen transcriptome of koi carp infected with Aeromonas sobria (A. sobria). De novo assembly of paired-end reads yielded 69,480 unigenes, of which the total length, average length, N50, and GC content are 70,120,028 bp, 1037 bp, 1793 bp, and 45.77%, respectively. Annotation is performed by comparison against various databases, yielding 42,229 (non-redundant protein sequence (NR): 60.78%), 59,255 (non-redundant nucleotide (NT): 85.28%), 35,900 (Swiss-Prot: 51.67%), 11,772 (clusters of orthologous groups (COG): 16.94%), 33,057 (Kyoto Encyclopedia of Genes and Genomes (KEGG): 47.58%), 18,764 (Gene Ontology (GO): 27.01%), and 32,085 (Interpro: 46.18%) unigenes. Comparative analysis of the expression profiles between bacterial challenge fish and control fish identifies 7749 differentially expressed genes (DEGs) from the muscle and 7846 DEGs from the spleen. These DEGs are further categorized with KEGG. Enrichment analysis of the DEGs and unigenes reveals major immune-related functions, including up-regulation of genes related with Toll-like receptor signaling, complement and coagulation cascades, and antigen processing and presentation. The results from RNA-Seq data are also validated and confirmed the consistency of the expression levels of seven immune-related genes after 24 h post infection with qPCR. Microsatellites (11,534), including di-to hexa nucleotide repeat motifs, are also identified. Altogether, this work provides valuable insights into the underlying immune mechanisms elicited during bacterial infection in koi carp that may aid in the future development of disease control measures in protection against A. sobria.

Applications of transcriptomics and proteomics in understanding fish immunity

With the development of intensive aquaculture, economic losses increasingly result from fish mortality due to pathogen infection. In recent years, a growing number of researchers have used transcriptomic and proteomic analyses to study fish immune responses to exogenous pathogen infection. Integrating transcriptomic and proteomic analyses provides a better understanding of the fish immune system including gene expression, regulation, and the intricate biological processes underlying immune responses against infection. This review focuses on the recent advances in the fields of transcriptomics and proteomics, which have contributed to our understanding of fish immunity to exogenous pathogens.

Molecular cloning and functional characterization of cathepsin B from Nile tilapia (Oreochromis niloticus)

Cathepsin B (CatB) has been widely known for its hydrolytic ability and involvement in the innate immunity. However, the mechanism of CatB from teleosts participating in immunoregulation remains poorly understood; and the sequence of CatB from Nile tilapia (NtCatB) has not been cloned and characterized. In this study, the coding sequence of NtCatB was cloned, and then characterized by bioinformatic analysis and heterologous expression. The deduced amino acid sequence (330-aa) of NtCatB contains the representative features of CatB. Quantitative real-time PCR revealed the extensive mRNA expression of NtCatB in six tissues of healthy Nile tilapia, and its transcription level was significantly up-regulated after Streptococcus agalactiae challenge. NtCatB may interact with some immunological function proteins and take part in the regulatory pathway. These results suggest that NtCatB is likely to be involved in the immune reaction. The mature region (residues 79-328, mNtCatB) of NtCatB was cloned and transferred to pET-28a for expressing the recombinant protein. The purified recombinant mNtCatB was verified with the activity of 992.34 U mg^-1 min^-1 under the optimal condition using a substrate hydrolyzing assay. The recombinant cystatin-A1-like can effectively inhibit the activity of the recombinant mNtCatB, and their binding form was predicted by molecular docking. Our results contribute to elucidating the immunological functions of NtCatB.

Basal polarization of the immune responses to Streptococcus agalactiae susceptible and resistant tilapia (Oreochromis niloticus)

One of the highest priority areas for improvement is the development of effective strategies for decreasing disease mortality levels in aquaculture production, a better understanding of the components of the fish immune system and their functions in the context of pathogen invasion is needed. Tilapia is the most common fish in South China, and Streptococcus agalactiae has become the most serious disease problem for tilapia industry in China. Here, we profiled gene expression differences between tilapia differing in their susceptibility to S. agalactiae both basally (before infection) and at three early timepoints post-infection (5 h, 50 h, and 7 d). Between group comparisons revealed 5756 unique genes differentially expressed greater than 2-fold at one or more timepoints. And the resistant fish showed much more strong ability in pathogen recognition, antigen presentation, immune activation, while the susceptible fish showed fast activation of apoptosis. Taken together, the immune profiles expand our knowledge for molecular mechanisms for disease resistance, as well as provide solid molecular resources for further identification of the candidate markers for disease-resistant selection and evaluation of disease prevention and treatment options for tilapia industry.

Identifying a Major QTL Associated with Salinity Tolerance in Nile Tilapia Using QTL-Seq

Selection of new lines with high salinity tolerance allows for economically feasible production of tilapias in brackish water areas. Mapping QTLs and identifying the markers linked to salinity-tolerant traits are the first steps in the improvement of the tolerance in tilapia through marker-assisted selection techniques. By using QTL-seq strategy and linkage-based analysis, two significant QTL intervals (chrLG4 and chrLG18) on salinity-tolerant traits were firstly identified in the Nile tilapia. Fine mapping with microsatellite and SNP markers suggested a major QTL region that located at 23.0 Mb of chrLG18 and explained 79% of phenotypic variation with a LOD value of 95. Expression analysis indicated that at least 10 genes (e.g., LACTB2, KINH, NCOA2, DIP2C, LARP4B, PEX5R, and KCNJ9) near or within the QTL interval were significantly differentially expressed in intestines, brains, or gills under 10, 15, or 20 ppt challenges. Our findings suggest that QTL-seq can be effectively utilized in QTL mapping of salinity-tolerant traits in fish. The identified major QTL is a promising locus to improve our knowledge on the genetic mechanism of salinity tolerance in tilapia.

Hepatic transcriptome analysis of juvenile GIFT tilapia (Oreochromis niloticus), fed diets supplemented with different concentrations of resveratrol

The GIFT (Genetically Improved Farmed Tilapia) tilapia, Oreochromis niloticus, is cultured widely for the production of freshwater fish in China. Streptococcosis, which is related to pathogenic infections, occurs frequently in juvenile and adult female GIFT individuals. Resveratrol (RES) has been used in feed to control these infections in freshwater tilapia. To address the effects of RES on tilapia, we used high-throughput RNA sequencing technology (RNA-Seq, HiSeq. 2500) to explore the global transcriptomic response and specific involvement of hepatic mRNA of juvenile O. niloticus fed with diets containing different concentrations of (0, 0.025, 0.05, and 0.1g/kg) RES. A total of > 24,513,018 clean reads were generated and then assembled into 23,244 unigenes. The unigenes were annotated by comparing them against non-redundant protein sequence (Nr), non-redundant nucleotide (Nt), Swiss-Prot, Pfam, Gene Ontology database (GO), Clusters of Orthologous Groups (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases, and 12,578 unigenes were annotated to the GO database. A total of 1444 (0.025g/kg RES), 1526 (0.05g/kg RES), and 3135 (0.1g/kg RES) genes were detected as significant differentially expressed genes (DEGs), when compared with the controls. A total of 6 (0.025 vs 0.05g/kg RES), 19 (0.025 vs 0.1g/kg RES), and 124 (0.05 vs 0.1g/kg RES) genes were detected as significant DEGs. Six genes, including dnah7x1, sox4, fam46a, hsp90a, ddit4, and nmrk2, were associated with an immune response. These findings provide information on the innate immunity of GIFT and might contribute to the development of strategies for the effective management of diseases and long-term sustainability of O. niloticus culture.

l-rhamnose-binding lectins (RBLs) in Nile tilapia, Oreochromis niloticus: Characterization and expression profiling in mucosal tissues

Rhamnose-binding lectins (RBLs) are crucial elements associated with innate immune responses to infections and have been characterized from a variety of teleost fishes. Given the importance of RBL in teleost fishes, we sought to study the diversity and expression profiles of RBLs in an important cultured fish, Nile tilapia (Oreochromis niloticus) following experimental infection with Streptococcus agalactiae, a major cause of streptococcosis in farmed tilapia. In this study, four predicted RBL genes were identified from Nile tilapia and were designated as OnRBL3a, OnRBL3b, OnRBL3c, and OnRBL3d. These OnRBLs were composed of two tandem-repeated type five carbohydrate recognition domains (CRDs), classified as type IIIc, and all clustered together phylogenetically. OnRBL-CRDs shared conserved topology of eight cysteine residues, characteristic peptide motifs of -YGR- and -DPC- (or -FGR- and -DTC-), and similar exon/intron organization. OnRBLs had the highest expression in immune-related tissues, gills, intestine or liver. However, the changes of OnRBL expression in the gills and intestine at 2 h, 4 h and 24 h post S. agalactiae challenge were modest, suggesting that tilapia may not mediate the entry or confront the infection of S. agalactiae through induction of RBL genes. The observed expression pattern may be related to the RBL type and CRD composition, S. agalactiae pathogenesis, the accessibility of ligands on the bacterial surface, and/or the species of fish. OnRBLs characterized in this study were the first RBL members identified in Nile tilapia and their characterization will expand our knowledge of RBLs in immunity.

Transcriptome analysis of hybrid tilapia (Oreochromis spp.) with Streptococcus agalactiae infection identifies Toll-like receptor pathway-mediated induction of NADPH oxidase complex and piscidins as primary immune-related responses

Streptococcus agalactiae infection is one of the most significant bacterial diseases in tilapia aquaculture. Identification of immune-related genes associated with Streptococcus agalactiae infection may provide a basis for breeding selection or therapeutics to augment disease resistance. Therefore, we utilized transcriptome profiling to study the host response in tilapia following Streptococcus agalactiae infection. Based on GO and KEGG enrichment analyses, we found that differentially expressed genes are widely involved in immune-related pathways, including the induction of antimicrobial peptides. Moreover, the main components of two immune-related pathways (Toll-like receptor signaling and leukocyte transendothelial migration) and four environmental information processing pathways (TNF, PI3K-Akt, Jak-STAT and MAPK) were identified. Finally, a time-course expression profile for several of the identified transcripts including tilapia piscidin 3 (TP3), tilapia piscidin 4 (TP4), TLR2, TLR5, MyD88, TRAF6, p38, and interleukin components was performed by qRT-PCR. Collectively, these results provide a starting point to study molecular mechanisms of tilapia immune response to Streptococcus agalactiae infection and may be applied as a basis for developing disease resistant strains by breeding selection.

Transcriptome analysis of immune response against Vibrio harveyi infection in orange-spotted grouper (Epinephelus coioides)

Vibrio harveyi is a gram-negative bacterium reported as found in many aquaculture species. To increase knowledge of the immune response against V. harveyi, in this study we performed transcriptome analysis of head kidney and spleen in orange-spotted grouper (Epinephelus coioides) at 1 and 2 days post-infection (dpi), using the Illumina sequencing platform. After de novo assembly, a total of 79,128 unigenes was detected with an N50 of 2511 bp. After alignments with sequences recorded in the major databases (NT, NR, Swiss-Prot COG, KEGG, Interpro and GO), based on sequence similarity, 61,208 (77.4%) of the unigene total could be annotated using at least one database. Comparison of gene expression levels between V. harveyi and a control group at each time point revealed differentially expressed genes (DEGs) (P < 0.05): a total of 7918 (5536 upregulated and 2282 downregulated genes) from head kidney at 1 day post infection (dpi), 4260 (1444 upregulated and 2816 downregulated genes) from head kidney at 2 dpi, 7887 (4892 upregulated and 2995 downregulated genes) from spleen at 1 dpi, and 8952 (7388 upregulated and 1564 downregulated genes) from spleen at 2 dpi. The DEGs were mainly annotated into signal transduction and immune system categories, based on the KEGG database. The DEGs were enriched in immune-related pathway functions, NOD-like receptor signaling pathways, Toll-like receptor signaling pathways, NF-κB signaling pathways, and Jak-STAT signaling pathways. Additionally, we selected several DEGs and validated their expression level by RT-qPCR. The data generated in this study may provide a valuable resource for further immune response research and offer improved strategies against V. harveyi infection in teleost fishes.

Genome-Wide QTL Analysis Identified Significant Associations Between Hypoxia Tolerance and Mutations in the GPR132 and ABCG4 Genes in Nile Tilapia

Exposure to hypoxia induces both acute and chronic stress responses, which plays an important role in health of cultured organisms including growth, reproduction, immunity, and other energy demanding activities. Application of advanced genomic technologies allows rapid identification of hypoxia trait-associated genes and precise selection of superior brood stocks with high tolerance in tilapia. By applying QTL-seq and double-digest restriction-site associated DNA sequencing (ddRAD-seq) techniques, we identified four genome-wide significant quantitative trait loci (QTLs) for hypoxia tolerance and many suggestive QTLs in Nile tilapia. These QTLs explained 6.6-14.7% of the phenotypic variance. Further analysis revealed that single nucleotide polymorphisms (SNPs) in exons of both GPR132 and ABCG4 genes located in genome-wide QTL intervals were significantly associated with hypoxia-tolerant traits. Expression analysis of both genes suggested that they were strong candidate genes involved into hypoxia tolerance in tilapia. Our findings suggest that both QTL-seq and ddRAD-seq techniques can be effectively utilized in QTL mapping of hypoxia traits in fish. Our data supply a basis for further marker-assisted selection of super lines with a high level of tolerance against low oxygen stress in the tilapia.

Omics and cytokine discovery in fish: Presenting the Yellowtail kingfish (Seriola lalandi) as a case study

A continued programme of research is essential to overcome production bottlenecks in any aquacultured fish species. Since the introduction of genetic and molecular techniques, the quality of immune research undertaken in fish has greatly improved. Thousands of species specific cytokine genes have been discovered, which can be used to conduct more sensitive studies to understand how fish physiology is affected by aquaculture environments or disease. Newly available transcriptomic technologies, make it increasingly easier to study the immunogenetics of farmed species for which little data exists. This paper reviews how the application of transcriptomic procedures such as RNA Sequencing (RNA-Seq) can advance fish research. As a case study, we present some preliminary findings using RNA-Seq to identify cytokine related genes in Seriola lalandi. These will allow in-depth investigations to understand the immune responses of these fish in response to environmental change or disease and help in the development of therapeutic approaches.

Transcriptomic analysis of the differentiating ovary of the protogynous ricefield eel Monopterus albus

BACKGROUND

The ricefield eel is a protogynous hermaphroditic Synbranchiform species that changes sex naturally from female to male, which offers an interesting model for studying gonadal (particularly ovarian) differentiation in vertebrates. In the present study, transcriptome sequencing of the gonad of ricefield eel larvae was performed to explore the molecular mechanisms underlying the ovarian differentiation and development.

RESULTS

A total of 301,267,988 clean reads were generated from cDNA libraries of gonadal tissues of ricefield eel larvae at 6, 9, 12, and 20 days post hatching (dph), which contained undifferentiated gonads, differentiating ovaries, ovaries with oogonia, and ovaries with meiotic oocytes, respectively. De-novo assembly of all the clean reads generated a total of 265,896 unigenes with a mean size of 720 bp and a N50 of 1107 bp. RT-qPCR analysis of the developmental expression of 13 gonadal development-related functional genes indicated that RNA-seq data are reliable. Transcriptome data suggest that high expression of female development-related genes and low expression of male development-related genes in the early gonads of ricefield eel larvae participate in the cascade of sex differentiation leading to the final female phenotype. The contrasting expression patterns of genes involved in retinoid acid (RA) synthesis and degradation might result in peak production of RA at 12 dph in the gonad of ricefield eel larvae, and induce molecular events responsible for the initiation of meiosis before the meiotic signs could be observed at 20 dph. In addition, only stra6 but not stra8 could be identified in gonadal transcriptome data of ricefield eel larvae, and the expression pattern of stra6 paralleled those of genes involved in RA synthesis, suggesting that stra6 may be a downstream target of RA and play a role in RA metabolism and/or meiotic initiation in the gonad of ricefield eel larvae.

CONCLUSIONS

The present study depicted the first large-scale RNA sequencing of the gonad of ricefield eel larvae, and identified many important functional genes, GO terms and KEGG pathways involved in gonadal development and germ cell meiosis. Results of the present study will facilitate future study on the ovarian differentiation of ricefield eels and other teleosts as well.

Transcriptomic profiling analysis of tilapia (Oreochromis niloticus) following Streptococcus agalactiae challenge

Innate immune system is the primary defense mechanism against pathogen infection in teleost, which are living in pathogen-rich aquatic environment. It has been long hypothesized that the disease resistance in teleost are strongly correlated to the activities of innate immune genes. Tilapia is an important economical fish around the world, especially in China, where the production accounts for nearly half of the global production. Recently, S. agalactiae has become one of the most serious bacterial diseases in southern China, resulted in high cumulative mortality and economic loss to tilapia industry. Therefore, we sought here to characterize the expression profiles of tilapia against S. agalactiae infection at whole transcriptome level by RNA-seq technology. A total of 2822 genes were revealed significantly expressed in tilapia spleen with a general trend of induction. Notably, most of the genes were rapidly the most induced at the early timepoint. The significantly changed genes highlighted the function of pathogen attachment and recognition, antioxidant/apoptosis, cytoskeletal rearrangement, and immune activation. Collectively, the induced expression patterns suggested the strong ability of tilapia to rapidly recognize the invasive bacteria, and activation of downstream immune signaling pathways to clear the bacteria and prevent the tissue damage and bacteria triggered cell apoptosis. Our results heighted important roles of novel candidate genes which were often missed in previous tilapia studies. Further studies are needed to characterize the molecular relationships between key immune genes and disease resistance, and to identify the candidate genes for molecular-assistant selection of disease-resistant broodstock and evaluation of disease prevention and treatment measures.

Transcriptome analysis of grey mullet (Mugil cephalus) after challenge with Lactococcus garvieae

Grey mullet (Mugil cephalus) is an economically important fish species in Taiwan mariculture industry. Moreover, grey mullet are common hosts of a bacterial infection by Lactococcus garvieae. However, until now the information related to the immune system of grey mullet is unclear. Therefore, to understand the molecular basis underlying the host immune response to L. garvieae infection, Illumina HiSeq™ 2000 was used to analyse the head kidney and spleen transcriptome of infected grey mullet. De novo assembly of paired-end reads yielded 55,203 unigenes. Comparative analysis of the expression profiles between bacterial challenge fish and control fish identified a total of 7192 from head kidney and 7280 in spleen differentially expressed genes (P < 0.05), including 4211 upregulated genes and 2981 downregulated genes in head kidney, while in spleen 3598 genes were upregulated and 3682 downregulated. A significant enrichment analysis of these differentially expressed genes (DEG) in spleen and head kidney revealed major immune-related pathways, including complement and coagulation cascades, Toll-like receptor signalling, and antigen processing and presentation. Moreover, selected DEGs were validated using qPCR. Altogether, the results obtained on immune-related genes may allow for a better understanding of immunity in grey mullet to Lactococcus garvieae, carrying out detailed functional analysis of these genes and developing strategies for efficient immune protection against infections in grey mullet.

Analysis of digital gene expression profiling in hemocytes of white shrimp Litopenaeus vannamei under nitrite stress

Accumulation of nitrite in water is highly toxic to aquatic animals. To understand immune responses in shrimp under such environmental stress, a digital gene expression (DGE) technology was applied to detect the gene expression profile of the Litopenaeus vannamei hemocytes in response to nitrite for 48 h. A total of 1922 differently expressed unigenes were generated. Of these transcripts, 1269 and 653 genes were up- or down-regulated respectively. Functional categorization and pathways of the differentially expressed genes revealed that immune defense, xenobiotics biodegradation and metabolism, amino acid and nucleobase metabolic process, apoptosis were the differentially regulated processes occurring during nitrite stress. We selected 19 differential expression transcripts (DETs) to validate the sequencing results by real time quantitative PCR (qPCR). The Pearson's correlation coefficient (R) of the 19 DETs was 0.843, which confirmed the consistency and accuracy between these two approaches. Subsequently, we screened 10 genes to examine the changes in the time course of gene expression in more detail. The results indicated that expressions of ATP-binding cassette transporter (ABC transporter), caspase10, QM protein, C type lectin 4 (CTL4), protein disulfide isomerase (PDI), serine protease inhibitor 8 (SPI8), transglutaminase (TGase), chitinase1, inhibitors of apoptosis proteins (IAP) and cytochrome P450 enzyme (CYP450) were induced to participate in the anti-stress defense against nitrite. These results will provide a reference for follow-up study of molecular toxicology and valuable gene information for better understanding of immune response in L. vannamei under environmental stress.

De Novo Transcriptome Analysis of Differential Functional Gene Expression in Largemouth Bass (Micropterus salmoides) after Challenge with Nocardia seriolae

Largemouth bass (Micropterus salmoides) are common hosts of an epizootic bacterial infection by Nocardia seriolae. We conducted transcriptome profiling of M. salmoides to understand the host immune response to N. seriolae infection, using the Illumina sequencing platform. De novo assembly of paired-end reads yielded 47,881 unigenes, the total length, average length, N50, and GC content of which were 49,734,288, 1038, 1983 bp, and 45.94%, respectively. Annotation was performed by comparison against non-redundant protein sequence (NR), non-redundant nucleotide (NT), Swiss-Prot, Clusters of Orthologous Groups (COG), Kyoto Encyclopaedia of Genes and Genomes (KEGG), Gene Ontology (GO), and Interpro databases, yielding 28,964 (NR: 60.49%), 36,686 (NT: 76.62%), 24,830 (Swissprot: 51.86%), 8913 (COG: 18.61%), 20,329 (KEGG: 42.46%), 835 (GO: 1.74%), and 22,194 (Interpro: 46.35%) unigenes. Additionally, 8913 unigenes were classified into 25 Clusters of Orthologous Groups (KOGs) categories, and 20,329 unigenes were assigned to 244 specific signalling pathways. RNA-Seq by Expectation Maximization (RSEM) and PossionDis were used to determine significantly differentially expressed genes (False Discovery Rate (FDR) < 0.05) and we found that 1384 were upregulated genes and 1542 were downregulated genes, and further confirmed their regulations using reverse transcription quantitative PCR (RT-qPCR). Altogether, these results provide information on immune mechanisms induced during bacterial infection in largemouth bass, which may facilitate the prevention of nocardiosis.

Comprehensive identification and profiling of Nile tilapia (Oreochromis niloticus) microRNAs response to Streptococcus agalactiae infection through high-throughput sequencing

MicroRNAs are a kind of small non-coding RNAs that participate in various biological processes. Deregulated microRNA expression is associated with several types of diseases. Tilapia (Oreochromis niloticus) is an important commercial fish species in China. To identify miRNAs and investigate immune-related miRNAs of O. niloticus, we applied high-throughput sequencing technology to identify and analyze miRNAs from tilapia infected with Streptococcus agalactiae at a timescale of 72 h divided into six different time points. The results showed that a total of 3009 tilapia miRNAs were identified, including in 1121 miRNAs which have homologues in the currently available databases and 1878 novel miRNAs. The expression levels of 218 tilapia miRNAs were significantly altered at 6 h-72 h post-bacterial infection (pi), and these miRNAs were therefore classified as differentially expressed tilapia miRNAs. For the 1121 differentially expressed tilapia miRNAs target 41961 genes. GO and KEGG enrichment analysis revealed that some target genes of tilapia miRNAs were grouped mainly into the categories of apoptotic process, signal pathway, and immune response. This is the first report of comprehensive identification of O. niloticus miRNAs being differentially regulated in spleen in normal conditions relating to S. agalactiae infection. This work provides an opportunity for further understanding of the molecular mechanisms of miRNA regulation in O. niloticus host-pathogen interactions.

De novo characterization of the liver transcriptome of javelin goby Synechogobius hasta and analysis of its transcriptomic profile following waterborne copper exposure

Previous studies have investigated the physiological responses to chronic copper (Cu) exposure in the liver of Synechogobius hasta; however, little information is available on the underlying molecular mechanisms. In an effort to better understand the mechanisms of Cu toxicity and to illuminate global gene expression patterns modulated by Cu exposure, we obtained the liver transcriptome information of S. hasta by RNA sequencing (RNA-seq) technology and also investigated the differential expression of genes following waterborne Cu exposure. Using the Illumina sequencing platform, as many as 60,217 unigenes were generated, with 815 bp of average length and 1298 bp of unigene N50 after filtering and assembly. For functional annotation analysis, 34,860, 31,526, 31,576, 25,808, 11,542, and 21,721 unigenes were annotated to the NR, NT, Swiss-Prot, KEGG, COG, and GO databases, respectively, and total annotation unigenes were 37,764. After 30 days of exposure to 55 μg Cu/l, a total of 292 and 1076 genes were significantly up- and down-regulated, respectively. By KEGG analysis, 660 had a specific pathway annotation. Subsequent bioinformatics analysis revealed that the differentially expressed genes were mainly related to lipid metabolism, immune system, apoptosis, and signal transduction, suggesting that these signaling pathways may be regulated by Cu exposure. The present study provides comprehensive sequence information for subsequent gene expression studies regarding S. hasta, and the transcriptome profiling after Cu exposure is also expected to improve our understanding of the molecular toxicology of Cu.

De Novo Transcriptome Analysis of Two Seahorse Species (Hippocampus erectus and H. mohnikei) and the Development of Molecular Markers for Population Genetics

Seahorse conservation has been performed utilizing various strategies for many decades, and the deeper understanding of genomic information is necessary to more efficiently protect the germplasm resources of seahorse species. However, little genetic information about seahorses currently exists in the public databases. In this study, high-throughput RNA sequencing for two seahorse species, Hippocampus erectus and H. mohnikei, was carried out, and de novo assembly generated 37,506 unigenes for H. erectus and 36,113 unigenes for H. mohnikei. Among them, 17,338 (46.23%) unigenes for H. erectus and 17,900 (49.57%) for H. mohnikei were successfully annotated based on the information available from the public databases. Through comparing the unigenes of two seahorse species, 7,802 candidate orthologous genes were identified and 5,268 genes among them could be annotated. In addition, gene ontology analysis of two species was similarly performed on biological processes, cellular components, and molecular functions. Twenty-four and twenty-one unigenes in H. erectus and H. mohnikei were annotated in the biosynthesis of unsaturated fatty acids pathways, and both seahorses lacked the Δ12 and Δ15 desaturases. Total of 8,992 and 9,116 SSR loci were obtained from H. erectus and H. mohnikei unigenes, respectively. Dozens of SSR were developed and then applied to assess the population genetic diversity, as well as cross-amplified in a related species, H. trimaculatus. The HO and HE values of the tested populations for H. erectus, H. mohnikei, and H. trimaculatus were medium. These resources would facilitate the conservation of the species through a better understanding of the genomics and comparative genome analysis within the Hippocampus genus.

Trancriptomic profiling revealed the signatures of acute immune response in tilapia (Oreochromis niloticus) following Streptococcus iniae challenge

Streptococcus iniae is the most significant bacterial disease of tilapia throughout the world, and commonly leads to tremendous economic losses. In contrast to other important fish species, our knowledge about the molecular mechanisms of tilapia in response to bacterial infection is still limited. Here, therefore, we utilized RNA-seq to first profiling of host responses in tilapia spleen following S. iniae infection at transcriptome level. A total of 223 million reads were obtained and assembled into 192,884 contigs with average length 844 bp. Gene expression analysis between control and infected samples at 5 h, 50 h, and 7 d revealed 1475 differentially expressed genes. In particular, the differentially expressed gene set was dramatically induced as early as 5 h, and rapidly declined to basal levels at 50 h. Enrichment and pathway analysis of the differentially expressed genes revealed the centrality of the pathogen attachment and recognition, cytoskeletal rearrangement and immune activation/inflammation in the pathogen entry and host inflammatory responses. Understanding of these responses can highlight mechanisms of tilapia host defense, and expand our knowledge of teleost immunology. Our findings will set a foundation of valuable biomarkers for future individual, strain, and family-level studies to evaluate immune effect of vaccine and individual response in host defense mechanisms to S. iniae infection, to select disease resistant families and strains.

Transcriptome analysis of the pearl oyster (Pinctada fucata) hemocytes in response to Vibrio alginolyticus infection

The pearl oyster Pinctada fucata is cultured widely for production of marine pearls in China, while mass mortalities, likely related to pathogenic infections, have occurred frequently in juvenile, mother and operated oysters. To address this issue, understanding host defense mechanisms of P. fucata against pathogenic challenge is extremely important. In the present study, a comparative analysis of hemocyte transcriptomes of P. fucata before and after Vibrio alginolyticus infection was conducted using the Illumina/Hiseq-2000 RNA-Seq technology. A total of 56,345,139 clean reads were generated and then assembled into 74,007 unigenes with an average length of 680 bp and an N50 of 1197 bp. Unigenes were annotated by comparing against non-redundant protein sequence (nr), non-redundant nucleotide (nt), Swiss-Prot, Pfam, Gene Ontology database (GO), Clusters of Orthologous Groups (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases, and 29,615 unigenes (40.01%) were annotated in at least one database. There were 636 genes (518 up-regulated and 118 down-regulated) that were significantly differentially expressed after bacterial challenge, and among which 369 were associated with 122 pathways, including classical immune-related pathways, such as 'MAPK signaling pathway', 'Chemokine signaling pathway', 'Apoptosis' and 'Wnt signaling pathway'. These findings provide information on the pearl oyster innate immunity and may contribute to developing strategies for management of diseases and long-term sustainability of P. fucata culture.