Microarray technology is an effective tool for identifying genes related to the aquacultural improvement of Japanese flounder, Paralichthys olivaceus

Genetics Responses to Hypoxia and Reoxygenation Stress in Larimichthys crocea Revealed via Transcriptome Analysis and Weighted Gene Co-Expression Network

Simple Summary

Hypoxia, which occurs frequently in aquaculture, can cause serious harm to all aspects of the growth, reproduction and metabolism of cultured fish. Due to the intolerance of Larimichthys crocea to hypoxia, Larimichthys crocea often floats head or even dies under hypoxic environment. However, the molecular mechanism of hypoxia tolerance in Larimichthys crocea has not been fully described. Therefore, the aim of this study was to explore the hub regulatory genes under hypoxic stress environment by transcriptome analysis of three key tissues (liver, blood and gill) in Larimichthys crocea. We identified a number of important genes that exercise different regulatory functions. Overall, this study will provide important clues to the molecular mechanisms of hypoxia tolerance in Larimichthys crocea.

Abstract

The large yellow croaker (Larimichthys crocea) is an important marine economic fish in China; however, its intolerance to hypoxia causes widespread mortality. To understand the molecular mechanisms underlying hypoxia tolerance in L. crocea, the transcriptome gene expression profiling of three different tissues (blood, gills, and liver) of L. crocea exposed to hypoxia and reoxygenation stress were performed. In parallel, the gene relationships were investigated based on weighted gene co-expression network analysis (WGCNA). Accordingly, the Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis showed that several pathways (e.g., energy metabolism, signal transduction, oxygen transport, and osmotic regulation) may be involved in the response of L. crocea to hypoxia and reoxygenation stress. In addition, also, four key modules (darkorange, magenta, saddlebrown, and darkolivegreen) that were highly relevant to the samples were identified by WGCNA. Furthermore, some hub genes within the association module, including RPS16, EDRF1, KCNK5, SNAT2, PFKL, GSK-3β, and PIK3CD, were found. This is the first study to report the co-expression patterns of a gene network after hypoxia stress in marine fish. The results provide new clues for further research on the molecular mechanisms underlying hypoxia tolerance in L. crocea.

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.

Microarray analysis of gene expression in olive flounder liver infected with viral haemorrhagic septicaemia virus (VHSV)

The most fatal viral pathogen in olive flounder Paralichthys olivaceus, is viral hemorrhagic septicemia virus, which afflicts over 48 species of freshwater and marine fish. Here, we performed gene expression profiling on transcripts isolated from VHSV-infected olive flounder livers using a 13 K cDNA microarray chip. A total of 1832 and 1647 genes were upregulated and down-regulated over two-fold, respectively, after infection. A variety of immune-related genes showing significant changes in gene expression were identified in upregulated genes through gene ontology annotation. These genes were grouped into categories such as antibacterial peptide, antigen-recognition and adhesion molecules, apoptosis, cytokine-related pathway, immune system, stress response, and transcription factor and regulatory factors. To verify the cDNA microarray data, we performed quantitative real-time PCR, and the results were similar to the microarray data. In conclusion, these results may be useful for the identification of specific genes or for the diagnosis of VHSV infection in flounder.

De Novo assembly of the Japanese flounder (Paralichthys olivaceus) spleen transcriptome to identify putative genes involved in immunity

Background

Japanese flounder (Paralichthys olivaceus) is an economically important marine fish in Asia and has suffered from disease outbreaks caused by various pathogens, which requires more information for immune relevant genes on genome background. However, genomic and transcriptomic data for Japanese flounder remain scarce, which limits studies on the immune system of this species. In this study, we characterized the Japanese flounder spleen transcriptome using an Illumina paired-end sequencing platform to identify putative genes involved in immunity.

Methodology/Principal Findings

A cDNA library from the spleen of P. olivaceus was constructed and randomly sequenced using an Illumina technique. The removal of low quality reads generated 12,196,968 trimmed reads, which assembled into 96,627 unigenes. A total of 21,391 unigenes (22.14%) were annotated in the NCBI Nr database, and only 1.1% of the BLASTx top-hits matched P. olivaceus protein sequences. Approximately 12,503 (58.45%) unigenes were categorized into three Gene Ontology groups, 19,547 (91.38%) were classified into 26 Cluster of Orthologous Groups, and 10,649 (49.78%) were assigned to six Kyoto Encyclopedia of Genes and Genomes pathways. Furthermore, 40,928 putative simple sequence repeats and 47, 362 putative single nucleotide polymorphisms were identified. Importantly, we identified 1,563 putative immune-associated unigenes that mapped to 15 immune signaling pathways.

Conclusions/Significance

The P. olivaceus transciptome data provides a rich source to discover and identify new genes, and the immune-relevant sequences identified here will facilitate our understanding of the mechanisms involved in the immune response. Furthermore, the plentiful potential SSRs and SNPs found in this study are important resources with respect to future development of a linkage map or marker assisted breeding programs for the flounder.

Analysis of immune-related ESTs and differential expression analysis of few important genes in lines of rohu (Labeo rohita) selected for resistance and susceptibility to Aeromonas hydrophila infection

A total of 137,629 contigs generated via de novo transcriptome assembly from resistant and susceptible lines of rohu (first generation) raised against aeromoniasis were further analyzed in terms of defence-related genes. Out of 1,939 contigs showing homology to genes involved in immune processes, 1,866 were further categorised into different functional subgroups. Comparative analysis revealed five genes for the first time in any carp species out of which apolipoprotein h, septin 4 isoform 3 and septin isoform cra_c were identified for the first time in fish. Differential expression analysis of ten genes viz., heat shock proteins (Hsps) (Hsp30, Hsp70 and Hsp90), serum lectin isoform 1 (SLI1), linker histone H1M (LHH1M), NAD(P)H quinone 1 (NQO1), zona pellucida 2 (ZP2) and three unknown genes that were highly up-expressed in first generation resistant line fish from mRNA-seq coverage data, was carried out using susceptible and resistant individuals of the second generation selected populations in eight different tissues viz. liver, kidney, intestine, gill, brain, spleen, skin and muscle using qPCR. Significant up-regulation in Hsp90, NQO1, C_116914 and C_22454 in specific tissues of resistant line and variable expression in Hsp30 and LHH1M genes in different tissues of both lines were noticed. The expression of Hsp70 was lower in many tissues of the resistant line than in susceptible line rohu. The expression of ZP2, SLI1 and C_94589 genes was not significantly different in terms of fold difference between the two lines. Differentially expressed genes need further characterisation to explore their role in resistance to Aeromonas hydrophila infection in rohu.

Geneious plugins for the access of PlasmoDB and PiroplasmaDB databases

The bioinformatics software, Geneious, provides a useful platform for researchers to retrieve and analyse genomic and functional genomics information. However, the main databases that the software is able to access are hosted by NCBI (National Center for Biotechnology Information). The databases of EuPathDB (Eukaryotic Pathogen Database Resources), such as PlasmoDB and PiroplasmaDB, collect more specific and detailed information about eukaryotic pathogens than those kept in NCBI databases. Two plugins for Geneious, one for PlasmaDB and one for PiroplasmaDB were developed. When installed, users can use search facilities to find and import gene and protein sequences from the EuPathDB databases. Users can then use the functions of Geneious to process the sequence information. When information unique to PlasmoDB and PiroplasmaDB is required, the user can access results linked with the gene/protein sequence via the default web browser. The plugins are freely available from the Victorian Bioinformatics Consortium website. The plugins can be modified to access any of the databases of EuPathDB.

Pathogenesis of and strategies for preventing Edwardsiella tarda infection in fish

Edwardsiella tarda is one of the serious fish pathogens, infecting both cultured and wild fish species. Research on edwardsiellosis has revealed that E. tarda has a broad host range and geographic distribution, and contains important virulence factors that enhance bacterial survival and pathogenesis in hosts. Although recent progress in edwardsiellosis research has enabled the development of numerous, highly effective vaccine candidates, these efforts have not been translated into a commercialized vaccine. The present review aims to provide an overview of the identification, pathology, diagnosis and virulence factors of E. tarda in fish, and describe recent strategies for developing vaccines against edwardsiellosis. The hope is that this presentation will be useful not only from the standpoint of understanding the pathogenesis of E. tarda, but also from the perspective of facilitating the development of effective vaccines.

Expressed sequences and polymorphisms in rohu carp (Labeo rohita, Hamilton) revealed by mRNA-seq

Expressed genes and polymorphisms were identified in lines of rohu Labeo rohita selected for resistance or susceptibility to Aeromonas hydrophila, an important bacterial pathogen causing aeromoniasis. All animals were grown in a common environment and RNA from ten individuals from each line pooled for Illumina mRNA-seq. De novo transcriptome assembly produced 137,629 contigs with 40× average coverage.Forty-four percent of the assembled sequences were annotated with gene names and ontology terms. Of these, 3,419 were assigned biological process terms related to "stress response" and 1,939 "immune system". Twenty-six contigs containing 38 single nucleotide polymorphisms (SNPs) were found to map to the Cyprinus carpio mitochondrial genome and over 26,000 putative SNPs and 1,700 microsatellite loci were detected. Seventeen percent of the 100 transcripts with coverage data most indicative of higher-fold expression(>5.6 fold) in the resistant line pool showed homology to major histocompatibility (MH), heat shock proteins (HSP)30, 70 and 90, glycoproteins or serum lectin genes with putative functions affecting immune response. Forty-one percent of these 100 transcripts showed no or low homology to known genes. Of the SNPs identified, 96 showing the highest allele frequency differences between susceptible and resistant line fish included transcripts with homology to MH class I and galactoside-binding soluble lectin, also with putative functions affecting innate and acquired immune response. A comprehensive sequence resource for L. rohita, including annotated microsatellites and SNPs from a mixture of A. hydrophila-susceptible and -resistant individuals, was created for subsequent experiments aiming to identify genes associated with A. hydrophila resistance.

Microarray-based identification of differentially expressed genes in families of turbot (Scophthalmus maximus) after infection with viral haemorrhagic septicaemia virus (VHSV)

Viral haemorrhagic septicaemia virus (VHSV) is one of the major threats to the development of the aquaculture industry worldwide. The present study was aimed to identify genes differentially expressed in several turbot (Scophthalmus maximus) families showing different mortality rates after VHSV. The expression analysis was conducted through genome-wide expression profiling with an oligo-microarray in the head kidney. A significant proportion of the variation in the gene expression profiles seemed to be explained by the genetic background, indicating that the mechanisms by which particular species and/or populations can resist a pathogen(s) are complex and multifactorial. Before the experimental infections, fish from resistant families (low mortality rates after VHSV infection) showed high expression of different antimicrobial peptides, suggesting that their pre-immune state may be stronger than fish of susceptible families (high mortality rates after VHSV infection). After infection, fish from both high- and low-mortality families showed an up-modulation of the interferon-induced Mx2 gene, the IL-8 gene and the VHSV-induced protein 5 gene compared with control groups. Low levels of several molecules secreted in the mucus were observed in high-mortality families, but different genes involved in viral entrance into target cells were down-regulated in low-mortality families. Moreover, these families also showed a strong down-modulation of marker genes related to VHSV target organs, including biochemical markers of renal dysfunction and myocardial injury. In general, the expression of different genes involved in the metabolism of sugars, lipids and proteins were decreased in both low- and high-mortality families after infection. The present study serves as an initial screen for genes of interest and provides an extensive overview of the genetic basis underlying the differences between families that are resistant or susceptible to VHSV infection.

Transcriptional analysis of orange-spotted grouper reacting to experimental grouper iridovirus infection

Disease caused by grouper iridovirus (GIV) has resulted in economic losses due to high mortality in grouper culture. Thirty-eight up- and 48 down-regulated known entities have been identified using a GIV-infected grouper kidney cDNA microarray chip. Further quantitative validation was executed in the head-kidney and spleen for 24 candidate genes and 7 immune factors following GIV inoculation. Significant induction with various patterns could be seen in 30 tested genes in the spleen. However, only 23 genes had induction in the head-kidney and meanwhile 5 genes showed reduction. Transcriptional expression profiles of selected genes in response to lipopolysaccharide (LPS) or polyinosinic:polycytidylic acid (PIC) were also established to compare with the GIV-stimulated expression. The results indicated that the responses of most genes facing GIV invasion have more similarities to PIC treatment than LPS. Seven genes are thought to be interferon-related factors: RNA helicase DHX58, ISG15, viperin, HECT E3 ligase (HECT), CD9, urokinase plasminogen activator surface receptor (PLAUR) and Mx-1. Following immunization with inactivated GIV, significant induction could be seen in DHX58, viperin, IL-1β, IL-8, COX-2, HECT, PLAUR, IgM, Mx-1, very large inducible GTPase-1 (VLIG1) and TNF-α in the head-kidney or spleen, and the latter 6 genes also had a gradual increasing pattern by a boosting immunization. These factors might play important roles in adaptive antiviral protection. Thus, we have characterized the temporal response patterns of virus responsive genes and have also identified several potential immune markers to further investigate host antiviral defense mechanisms.

Gene expression profiling in respond to TBT exposure in small abalone Haliotis diversicolor

In this study, we investigated the gene expression profiling of small abalone, Haliotis diversicolor by tributyltin (TBT) exposure using a cDNA microarray containing 2473 unique transcripts. Totally, 107 up-regulated genes and 41 down-regulated genes were found. For further investigation of candidate genes from microarray data and EST analysis, quantitative real-time PCR was performed at 6 h, 24 h, 48 h, 96 h and 192 h TBT exposure. 26 genes were found to be significantly differentially expressed in different time course, 3 of them were unknown. Some gene homologues like cellulose, endo-beta-1,4-glucanase, ferritin subunit 1 and thiolester containing protein II CG7052-PB might be the good biomarker candidate for TBT monitor. The identification of stress response genes and their expression profiles will permit detailed investigation of the defense responses of small abalone genes.