Efficient assembly and annotation of the transcriptome of catfish by RNA-Seq analysis of a doubled haploid homozygote

Transcript level of telomerase reverse-transcriptase (TERT) gene in the rainbow trout (Oncorhynchus mykiss) eggs with different developmental competence for gynogenesis

Expression of the telomerase reverse-transcriptase (TERT) gene and activity of telomerase have been reported in the somatic tissues and gonads in fish irrespective of their age and size. Nevertheless, little is known about TERT expression in the fish eggs. In the current study, the presence of the TERT transcripts was confirmed in the rainbow trout ovulated eggs before and after activation with nonirradiated and UV-irradiated (gynogenesis) sperm. Eggs originating from eight females had high and comparable quality expressed by similar hatching rates. However, survival of the gynogenetic larvae that hatched from eggs activated with UV-irradiated sperm and further exposed to the high hydrostatic pressure (HHP) shock for duplication of the maternal chromosomes varied between females from 2.1 ± 0.4 to 40.5 ± 2.2%. Increased level of TERT transcripts was observed in eggs originating from two females, and gametes from only one of them showed improved competence for gynogenesis (27.3 ± 1.9%). In turn, eggs from the female that exhibited the highest survival after gynogenetic activation were characterized by the lowest expression of the TERT gene. Telomerase in rainbow trout eggs may compensate erosion of the telomeres during early embryonic development; however, its upregulation does not assure better development after gynogenetic activation.

New reports on iron related proteins: Molecular characterization of two ferroportin genes in common carp (Cyprinus carpio L.) and its expression pattern

Iron uptake, transport, and storage require the involvement of several proteins, including ferroportin (fpn), the sole known iron efflux transporter. Due to its critical function fpn has been studied, particularly in humans. Here, we characterized the ferroportin gene in common carp (Cyprinus carpio L.) and performed RNA-seq analysis to evaluate its constitutive transcription levels across different tissues. Our results indicate that C. carpio possesses two functional fpns with distinct expression patterns, highlighting the potential for functional divergence and expression differentiation among fpns in this species.

Quality of fish eggs and production of androgenetic and gynogenetic doubled haploids (DHs)

Induced development of haploid embryos (H) with only paternal (androgenesis) or maternal (gynogenesis) chromosomes requires irradiation of eggs before fertilization or activation of eggs with irradiated spermatozoa, respectively. To provide doubled haploids (DHs), androgenetic and gynogenetic haploid zygotes need to be subjected to the thermal or high hydrostatic pressure (HHP) shock to suppress the first mitotic cleavage and to double paternal or maternal haploid set of chromosomes. Androgenesis and mitotic gynogenesis (mito-gynogenesis) result in the generation of fully homozygous individuals in a single generation. DHs have been utilized in selective breeding programs, in studies concerning the phenotypic consequences of recessive alleles and to evaluate the impact of sex chromosomes on the early ontogeny. Moreover, the use of DHs for the NGS approach radically improves de novo the assembly of the genomes. However, reduced survival of the doubled haploids limits the wide application of androgenotes and gynogenotes. The high mortality of DHs may be only partly explained by the expression of recessive traits. Observed inter-clutch variation in the survival of DHs developing in eggs originating from different females make it necessary to take a closer look at the quality of the eggs used during induced androgenesis and gynogenesis. Moreover, the developmental competence of eggs that are subjected to irradiation before fertilization in order to deactivate maternal chromosomes when undergoing induced androgenesis and exposed to the physical shock after fertilization that leads to the duplication of the zygotes in both mito-gynogenesis and androgenesis may be also altered as irradiation and sublethal values of temperatures and hydrostatic pressure are considered as harmful for the cell organelles and biomolecules. Here, recently provided results concerning the morphological, biochemical, genomic, and transcriptomic characteristics of fish eggs showing high and low competence for androgenesis and mito-gynogenesis are reviewed.

Molecular identification and functional characterization of IRF4 from common carp (Cyprinus carpio. L) in immune response: a negative regulator in the IFN and NF-κB signalling pathways

Background

The interferon (IFN) regulatory factors (IRFs) were originally identified as transcription factors playing critical roles in the regulation of IFN-related genes in the signal pathway. In mammals, IRF4 plays a vital role in both the innate and adaptive immune system. This study aims to reveal the molecular characterization, phylogenetic analysis, expression profiles and the regulatory role in the IFN and NF-κB signalling pathways of IRF4 in common carp (Cyprinus carpio. L) (abbreviation, ccIRF4).

Results

Here, ccIRF4 was identified and characterized, it contained a DNA binding domain (DBD) which possess five tryptophans and an IRF-associated domain (IAD). The predicted protein sequence of the ccIRF4 showed higher identities with grass carp (Ctenopharyngodon idella) and zebrafish (Danio rerio). Phylogenetic analysis suggested that ccIRF4 has the closest relationship with zebrafish IRF4. Quantitative real-time PCR analysis showed that ccIRF4 was constitutively expressed in all investigated tissues with the highest expression level in the gonad. Polyinosinic:polycytidylic acid (poly I:C) stimulation up-regulated the ccIRF4 expressions in the liver, spleen, head kidney, skin, foregut and hindgut. Upon Aeromonas hydrophila injection, the expression level of ccIRF4 was up-regulated in all tissues with the exception of spleen. In addition, ccIRF4 was induced by lipopolysaccharide (LPS), peptidoglycan (PGN) and Flagellin in head kidney leukocytes (HKLs). Overexpression of the ccIRF4 gene in epithelioma papulosum cyprini cells (EPC) down regulated the expressions of IFN-related genes and proinflammatory factors. Dual-luciferase reporter assay revealed that ccIRF4 decreased the activation of NF-κB through MyD88.

Conclusions

These results indicate that ccIRF4 participates in both antiviral and antibacterial immune response and negatively regulates the IFN and NF-κB response. Overall, our study on ccIRF4 provides more new insights into the innate immune system of common carp as well as a theoretical basis for investigating the pathogenesis and prevention of fish disease.

Toll-Like Receptor Evolution: Does Temperature Matter?

Toll-like receptors (TLRs) recognize conserved pathogen-associated molecular patterns (PAMPs) and are an ancient and well-conserved group of pattern recognition receptors (PRRs). The isolation of the Antarctic continent and its unique teleost fish and microbiota prompted the present investigation into Tlr evolution. Gene homologues of tlr members in teleosts from temperate regions were present in the genome of Antarctic Nototheniidae and the non-Antarctic sister lineage Bovichtidae. Overall, in Nototheniidae apart from D. mawsoni, no major tlr gene family expansion or contraction occurred. Instead, lineage and species-specific changes in the ectodomain and LRR of Tlrs occurred, particularly in the Tlr11 superfamily that is well represented in fish. Positive selective pressure and associated sequence modifications in the TLR ectodomain and within the leucine-rich repeats (LRR), important for pathogen recognition, occurred in Tlr5, Tlr8, Tlr13, Tlr21, Tlr22, and Tlr23 presumably associated with the unique Antarctic microbiota. Exposure to lipopolysaccharide (Escherichia coli O111:B4) Gram negative bacteria did not modify tlr gene expression in N. rossii head–kidney or anterior intestine, although increased water temperature (+4°C) had a significant effect.

A Comprehensive Annotation of the Channel Catfish (Ictalurus punctatus) T Cell Receptor Alpha/Delta, Beta, and Gamma Loci

The complete germline repertoires of the channel catfish, Ictalurus punctatus, T cell receptor (TR) loci, TRAD, TRB, and TRG were obtained by analyzing genomic data from PacBio sequencing. The catfish TRB locus spans 214 kb, and contains 112 TRBV genes, a single TRBD gene, 31 TRBJ genes and two TRBC genes. In contrast, the TRAD locus is very large, at 1,285 kb. It consists of four TRDD genes, one TRDJ gene followed by the exons for TRDC, 125 TRAJ genes and the exons encoding the TRAC. Downstream of the TRAC, are 140 TRADV genes, and all of them are in the opposite transcriptional orientation. The catfish TRGC locus spans 151 kb and consists of four diverse V-J-C cassettes. Altogether, this locus contains 15 TRGV genes and 10 TRGJ genes. To place our data into context, we also analyzed the zebrafish TR germline gene repertoires. Overall, our findings demonstrated that catfish possesses a more restricted repertoire compared to the zebrafish. For example, the 140 TRADV genes in catfish form eight subgroups based on members sharing 75% nucleotide identity. However, the 149 TRAD genes in zebrafish form 53 subgroups. This difference in subgroup numbers between catfish and zebrafish is best explained by expansions of catfish TRADV subgroups, which likely occurred through multiple, relatively recent gene duplications. Similarly, 112 catfish TRBV genes form 30 subgroups, while the 51 zebrafish TRBV genes are placed into 36 subgroups. Notably, several catfish and zebrafish TRB subgroups share ancestor nodes. In addition, the complete catfish TR gene annotation was used to compile a TR gene segment database, which was applied in clonotype analysis of an available gynogenetic channel catfish transcriptome. Combined, the TR annotation and clonotype analysis suggested that the expressed TRA, TRB, and TRD repertoires were generated by different mechanisms. The diversity of the TRB repertoire depends on the number of TRBV subgroups and TRBJ genes, while TRA diversity relies on the many different TRAJ genes, which appear to be only minimally trimmed. In contrast, TRD diversity relies on nucleotide additions and the utilization of up to four TRDD segments.

Nile tilapia Toll-like receptor 7 subfamily: Intracellular TLRs that recruit MyD88 as an adaptor and activate the NF-κB pathway in the immune response

Toll-like receptor 7 (TLR7) subfamily members are important pattern recognition receptors that participate in the recognition of pathogen-associated molecular patterns. In the present study, three TLR family members, OnTLR7, OnTLR8 and OnTLR9, were identified in the Nile tilapia Oreochromis niloticus. TLR7-, TLR8-and TLR9-deduced proteins have typical structural characteristics of TLRs, including Toll/interleukin-1 receptor (TIR), leucine-rich repeat (LRR) and transmembrane region (TM). OnTLR7, OnTLR8 and OnTLR9 were broadly expressed in all of the tissues tested, with the highest expression levels in the brain (TLR7) and spleen (TLR8 and TLR9). Moreover, the expression levels of OnTLR7, OnTLR8 and OnTLR9 were significantly increased in most tested tissues after Streptococcus agalactiae infection in vivo. After LPS stimulation, OnTLR7 and OnTLR9 mRNA expression levels were downregulated in the intestine and upregulated in the liver, spleen and kidney; however, OnTLR8 mRNA expression levels were upregulated in the kidney only after LPS stimulation for 5 d. After Poly I:C stimulation, OnTLR7 and OnTLR9 mRNA expression levels were upregulated in the intestine, liver, spleen and kidney, and the highest expression was found in the liver, while OnTLR8 mRNA expression levels were upregulated in the intestine, liver and kidney and downregulated in the spleen. Subcellular localization of OnTLR7, OnTLR8, and OnTLR9 in 293T cells showed that OnTLR9 was distributed in both the cytoplasm and nucleus while OnTLR8 and OnTLR7 were distributed mainly in the cytoplasm. Overexpression of OnTLR7, OnTLR8 and OnTLR9 in 293T cells had no significant effect on the activity of NF-κB, but they could significantly enhance MyD88-mediated NF-κB activity after cotransfection with MyD88. Pulldown assays showed that OnTLR7, OnTLR8, and OnTLR9 could interact with OnMyD88. Taken together, these results indicate that TLR7 subfamily genes play a role in the immune response to pathogen invasion of Nile tilapia.

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

Atlantic salmon (Salmo salar) is a major species produced in world aquaculture and an important vertebrate model organism for studying the process of rediploidization following whole genome duplication events (Ss4R, 80 mya). The current Salmo salar transcriptome is largely generated from genome sequence based in silico predictions supported by ESTs and short-read sequencing data. However, recent progress in long-read sequencing technologies now allows for full-length transcript sequencing from single RNA-molecules. This study provides a de novo full-length mRNA transcriptome from liver, head-kidney and gill materials. A pipeline was developed based on Iso-seq sequencing of long-reads on the PacBio platform (HQ reads) followed by error-correction of the HQ reads by short-reads from the Illumina platform. The pipeline successfully processed more than 1.5 million long-reads and more than 900 million short-reads into error-corrected HQ reads. A surprisingly high percentage (32%) represented expressed interspersed repeats, while the remaining were processed into 71 461 full-length mRNAs from 23 071 loci. Each transcript was supported by several single-molecule long-read sequences and at least three short-reads, assuring a high sequence accuracy. On average, each gene was represented by three isoforms. Comparisons to the current Atlantic salmon transcripts in the RefSeq database showed that the long-read transcriptome validated 25% of all known transcripts, while the remaining full-length transcripts were novel isoforms, but few were transcripts from novel genes. A comparison to the current genome assembly indicates that the long-read transcriptome may aid in improving transcript annotation as well as provide long-read linkage information useful for improving the genome assembly. More than 80% of transcripts were assigned GO terms and thousands of transcripts were from genes or splice-variants expressed in an organ-specific manner demonstrating that hybrid error-corrected long-read transcriptomes may be applied to study genes and splice-variants expressed in certain organs or conditions (e.g., challenge materials). In conclusion, this is the single largest contribution of full-length mRNAs in Atlantic salmon. The results will be of great value to salmon genomics research, and the pipeline outlined may be applied to generate additional de novo transcriptomes in Atlantic Salmon or applied for similar projects in other species.

Identifying the Related Genes of Muscle Growth and Exploring the Functions by Compensatory Growth in Mandarin Fish (Siniperca chuatsi)

How organisms display many different biochemical, physiological processes through genes expression and regulatory mechanisms affecting muscle growth is a central issue in growth and development. In Siniperca chuatsi, the growth-related genes and underlying relevant mechanisms are poorly understood, especially for difference of body sizes and compensatory growth performance. Muscle from 3-month old individuals of different sizes was used for transcriptome analysis. Results showed that 8,942 different expression genes (DEGs) were identified after calculating the RPKM. The DEGs involved in GH-IGF pathways, protein synthesis, ribosome synthesis and energy metabolisms, which were expressed significantly higher in small individuals (S) than large fish (L). In repletion feeding and compensatory growth experiments, eight more significant DEGs were used for further research (GHR2, IGFR1, 4ebp, Mhc, Mlc, Myf6, MyoD, troponin). When food was plentiful, eight genes participated in and promoted growth and muscle synthesis, respectively. Starvation can be shown to inhibit the expression of Mhc, Mlc and troponin, and high expression of GHR2, IGFR1, and 4ebp inhibited growth. Fasting promoted the metabolic actions of GHR2, IGFR1, and 4ebp rather than the growth-promoting actions. MyoD can sense and regulate the hunger, which also worked with Mhc and Mlc to accelerate the compensatory growth of S. chuatsi. This study is helpful to understand the regulation mechanisms of muscle growth-related genes. The elected genes will contribute to the selective breeding in future as candidate genes.

Comparative analysis of growth performance and liver transcriptome response of juvenile Ancherythroculter nigrocauda fed diets with different protein levels

This study aimed at investigating the effects of dietary protein levels on the growth and liver transcriptome in juvenile Ancherythroculter nigrocauda. Six semi-purified diets were formulated containing 25 (control), 30, 35, 40, 45, and 50% protein. Each diet was fed to three groups of 35 fish (mean initial weight: 5.86 ± 0.10 g) for 56 days. The rate of weight gain and specific growth rate increased with dietary protein levels from 25% to 40%, but remained unchanged when fed with 45 or 50% dietary protein. The feed conversion ratio was significantly influenced by the dietary protein levels, being the lowest in fish fed 40% protein. Illumina RNA-seq analysis was performed to investigate liver gene expression changes under different dietary protein treatments. A total of 367.78 million clean reads were obtained from the six libraries. Compared with 25% protein treatment library, there were 734, 1946, 1755, 2726, and 1523 upregulated genes, and 407, 1882, 1865, 2216 and 1624 downregulated genes in the 30, 35, 40, 45, and 50% protein treatment libraries, respectively. Trend analysis of these differentially expressed genes (DEGs) identified six statistically significant trends. A series of DEGs that related to protein metabolism, growth and development, lipid metabolism and immune and stress response were identified. Moreover, gene ontology enrichment analysis of the DEGs demonstrated that cellular process, single-organism process, metabolic process and biological regulation were the most highly overrepresented biological processes. Kyoto Encyclopedia of Genes and Genomes enrichment analysis revealed that protein processing in endoplasmic reticulum, PPAR signaling pathway, complement and coagulation cascades, and cytochrome P450 (CYP450s) were significantly enriched in the dietary protein treatment groups. Furthermore, qPCR results showed excellent agreement on those of RNA-seq for both up- and down-regulated genes (including fasn, accα, SCD, CPT-I, igf1, ST, AST, trdmt1, hsp70, cyp450, MHC-II, C4, tgfβ, ube4b, apoE and abcb7). Thus, our results provide the baseline information for the feed formulation and nutritional research for A. nigrocauda.

Divergence, evolution and adaptation in ray-finned fish genomes

With the rapid development of next-generation sequencing technologies and bioinformatics, over 50 ray-finned fish genomes by far have been sequenced with high quality. The genomic work provides abundant genetic resources for deep understanding of divergence, evolution and adaptation in the fish genomes. They are also instructive for identification of candidate genes for functional verification, molecular breeding, and development of novel marine drugs. As an example of other omics data, the Fish-T1K project generated a big database of fish transcriptomes to integrate with these published fish genomes for potential applications. In this review, we highlight the above-mentioned recent investigations and core topics on the ray-finned fish genome research, with a main goal to obtain a deeper understanding of fish biology for theoretical and practical applications.

Polyadenylation sites and their characteristics in the genome of channel catfish (Ictalurus punctatus) as revealed by using RNA-Seq data

Polyadenylation plays important roles in gene expression regulation in eukaryotes, which typically involves cleavage and poly(A) tail addition at the polyadenylation site (PAS) of the pre-mature mRNA. Many eukaryotic genes contain more than one PASs, termed as alternative polyadenylation (APA). As a crucial post-transcriptional regulation, polyadenylation affects various aspects of RNA metabolism such as mRNA stability, translocation, and translation. However, polyadenylation has been rarely studied in teleosts. Here we conducted polyadenylation analysis in channel catfish, a commercially important aquaculture species around the world. Using RNA-Seq data, we identified 20,320 PASs which were classified into 14,500 clusters by merging adjacent PASs. Most of the PASs were found in 3' UTRs, followed by intron regions based on the annotation of channel catfish reference genome. No apparent difference in PAS distribution was observed between the sense and antisense strand of the channel catfish genome. The sequence analysis of nucleotide composition and motif around PASs yielded a highly similar profile among various organisms, suggesting the conservation and importance of polyadenylation in evolution. Using APA genes with more than two PASs, gene ontology enrichment revealed genes particularly involved in RNA binding. Reactome pathway analysis showed the enrichment of the innate immune system, especially neutrophil degranulation.

A high-throughput SNP discovery strategy for RNA-seq data

Background

Single nucleotide polymorphisms (SNP) have been applied as important molecular markers in genetics and breeding studies. The rapid advance of next generation sequencing (NGS) provides a high-throughput means of SNP discovery. However, SNP development is limited by the availability of reliable SNP discovery methods. Especially, the optimum assembler and SNP caller for accurate SNP prediction from next generation sequencing data are not known.

Results

Herein we performed SNP prediction based on RNA-seq data of peach and mandarin peel tissue under a comprehensive comparison of two paired-end read lengths (125 bp and 150 bp), five assemblers (Trinity, IDBA, oases, SOAPdenovo, Trans-abyss) and two SNP callers (GATK and GBS). The predicted SNPs were compared with the authentic SNPs identified via PCR amplification followed by gene cloning and sequencing procedures. A total of 40 and 240 authentic SNPs were presented in five anthocyanin biosynthesis related genes in peach and in nine carotenogenic genes in mandarin. Putative SNPs predicted from the same RNA-seq data with different strategies led to quite divergent results. The rate of false positive SNPs was significantly lower when the paired-end read length was 150 bp compared with 125 bp. Trinity was superior to the other four assemblers and GATK was substantially superior to GBS due to a low rate of missing authentic SNPs. The combination of assembler Trinity, SNP caller GATK, and the paired-end read length 150 bp had the best performance in SNP discovery with 100% accuracy both in peach and in mandarin cases. This strategy was applied to the characterization of SNPs in peach and mandarin transcriptomes.

Conclusions

Through comparison of authentic SNPs obtained by PCR cloning strategy and putative SNPs predicted from different combinations of five assemblers, two SNP callers, and two paired-end read lengths, we provided a reliable and efficient strategy, Trinity-GATK with 150 bp paired-end read length, for SNP discovery from RNA-seq data. This strategy discovered SNP at 100% accuracy in peach and mandarin cases and might be applicable to a wide range of plants and other organisms.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5533-4) contains supplementary material, which is available to authorized users.

Liver transcriptome analysis of the Sparus macrocephlus in response to Vibrio parahaemolyticus infection

The black seabream (Sparus macrocephlus) is an economically pivotal aquaculture species cultured in China and Southeast Asian countries. To understand the molecular immune mechanisms underlying the response to Vibrio parahaemolyticus, a comparative gene transcription analysis were performed with utilized fresh livers of V. parahaemolyticus-immunized Sparus macrocephlus with a control group through RNA-Seq technology. A total of 256663 contigs were obtained after excluded the low-quality sequences and assembly. The average length of contigs collected from this research is 1066.93 bp. Furthermore, blast analysis indicates 30747 contigs were annotated based on homology with matches in the NT, NR, gene, and string databases. A gene ontology analysis was employed to classify 21598 genes according to three major functional categories: molecular function, cellular component, and biological process. A total of 14470 genes were discovered in 303 KEGG pathways. RSEM and EdgeR were introduced to estimate 3841 genes significantly different expressed (False Discovery Rate<0.001) which includes 4072 up-regulated genes and 3771 down-regulated genes. A significant enrichment analysis of these differentially expressed genes and isogenes were conducted to reveal the major immune-related pathways which refer to the toll-like receptor, complement, coagulation cascades, and chemokine signaling pathways. In addition, 92175 potential simple sequence repeats (SSRs) and 121912 candidate single nucleotide polymorphisms (SNPs) were detected and identified sequencely in the Sparus macrocephlus liver transcriptome. This research characterized a gene expression pattern for normal and the V. parahaemolyticus -immunized Sparus macrocephlus for the first time and not only sheds new light on the molecular mechanisms underlying the host-V. parahaemolyticus interaction but contribute to facilitate future studies on Sparus macrocephlus gene expression and functional genomics.

Increased Alternative Splicing as a Host Response to Edwardsiella ictaluri Infection in Catfish

Alternative splicing is the process of generating multiple transcripts from a single pre-mRNA used by eukaryotes to regulate gene expression and increase proteomic complexity. Although alternative splicing profiles have been well studied in mammalian species, they have not been well studied in aquatic species, especially after biotic stresses. In the present study, genomic information and RNA-Seq datasets were utilized to characterize alternative splicing profiles and their induced changes after bacterial infection with Edwardsiella ictaluri in channel catfish (Ictalurus punctatus). A total of 27,476 alternative splicing events, derived from 9694 genes, were identified in channel catfish. Exon skipping was the most abundant while mutually exclusive exon was the least abundant type of alternative splicing. Alternative splicing was greatly induced by E. ictaluri infection with 21.9% increase in alternative splicing events. Interestingly, genes involved in RNA binding and RNA splicing themselves were significantly enriched in differentially alternatively spliced genes after infection. Sequence analyses of splice variants of a representative alternatively spliced gene, splicing factor srsf2, revealed that certain spliced transcripts may undergo nonsense-mediated decay (NMD), suggesting functional significance of the induced alternative splicing. Although statistical analysis was not possible with such large datasets, results from quantitative real-time PCR from representative differential alternative splicing events provided general validation of the bacterial infection-induced alternative splicing. This is the first comprehensive study of alternative splicing and its changes in response to bacterial infection in fish species, providing insights into the molecular mechanisms of host responses to biotic stresses.

Molecular cloning, expression and antibacterial activity of goose-type lysozyme gene in Microptenus salmoides

It is well known that lysozymes are key proteins to teleosts in the innate immune system and possess high bactericidal properties. In the present study, a g-type lysozyme gene was cloned from Microptenus salmoides. The g-type sequence consisted of 582 bp, which translated into a 193 amino acid (AA) protein (GenBank accession no: MH087462). The predicted molecular weight and theoretical isoelectric point were 21.36 kDa and 6.91 respectively and no signal peptide was observed. The qRT-PCR analysis showed that the g-type lysozyme gene was differentially expressed in various tissues under normal conditions and the highest g-type lysozyme level was observed in liver, gill and spleen while there seemed to be low expression in the muscle, heart and head-kidney. The expression of g-type lysozyme was differentially upregulated in the spleen, gill and intestine after stimulation with heat stress and Aeromonas hydrophila (A. hydrophila). Under heat stress and A. hydrophila injection, the g-type lysozyme mRNA levels all in spleens, gill and intestine tissues increased significantly (P < 0.05), with the maximum levels attained at 12 h, 24 h (or 12 h) and 24 h. Thereafter, they all decreased significantly (P < 0.01) and the expression in gill returned to nearly the basal value within 72 h. Those results suggested that g-type lysozyme was involved in the immune response to heat stress and bacterial challenge. The cloning and expression analysis of the g-type lysozyme provide theoretical basis to further study the mechanism of anti-adverseness in Microptenus salmoides. The g-type lysozyme gene perhaps also played an important role in the immune responses against bacterial invasion.

Hypothalamus and pituitary transcriptome profiling of male and female Hong Kong grouper (Epinephelus akaara)

Hong Kong grouper (Epinephelus akaara) is an important commercially cultured marine fish in Asia, and a protogynous hermaphrodite with the "diandry" pattern. In order to explore the gene expression patterns of hypothalamus and pituitary between male and female Hong Kong grouper, we used RNA-seq technology to investigate transcriptomes of both tissues in immature and mature male and female adults. This produced 227,227,148 and 215,858,948 high quality reads from hypothalamus and pituitary, which were jointly assembled into 199,203 unigenes. Among them, 30,786 unigenes were mapped to known genes. Differential expression analysis revealed 275 unigenes that were differentially expressed between immature male and female adults and 561 between mature male and female adults. According to annotation and KEGG information, these differentially expressed genes (DEGs) were involved in development, metabolism, and regulation of transcription. One DEG, amino-terminal enhancer of split (AES), a member of the Groucho/transducin-like enhancer of split family of transcriptional regulators that played important roles in neurogenesis, segmentation, and sex determination, was significantly upregulated in male individuals in both immature and mature adult comparisons, indicating it may be involved in male reproductive function during development. Our report, for the first time, uses RNA-seq technology to investigate transcriptomes of both hypothalamus and pituitary in teleost fish, and provides a basis for further studies of molecular mechanism of sex determination and development in Hong Kong grouper.

MicroRNA regulation of Toll-like receptor signaling pathways in teleost fish

The innate immune system is the first line defense mechanism that recognizes, responds to, controls or eliminates invading pathogens. Toll-like receptors (TLRs) are a critical family of pattern recognition receptors (PRRs) tightly regulated by complex mechanisms involving many molecules to ensure a beneficial outcome in response to foreign invaders. MicroRNAs (miRNAs), a transcriptional and posttranscriptional regulator family in a wide range of biological processes, have been identified as new molecules related to the regulation of TLR-signaling pathways in immune responses. To date, at least 22 TLR types have been identified in more than a dozen different fish species. However, the functions and underlying mechanisms of miRNAs in the regulation of inflammatory responses related to the TLR-signaling pathway in fish is lacking. In this review, we summarize the regulation of miRNA expression profiles in the presence of TLR ligands or pathogen infections in teleost fish. We focus on the effects of miRNAs in regulating TLR-signaling pathways by targeting multiple molecules, including TLRs themselves, TLR-associated signaling proteins, and TLR-induced cytokines. An understanding of the relationship between the TLR-signaling pathways and miRNAs may provide new insights for drug intervention to manipulate immune responses in fish.

JAK and STAT members in channel catfish: Identification, phylogenetic analysis and expression profiling after Edwardsiella ictaluri infection

The Janus kinase/signal transducers and activators of transcription (JAK/STAT) signaling pathway is one of the main pleiotropic cascades used to transmit information from extracellular receptors to the nucleus, which results in DNA transcription and expression of genes involved in immunity, proliferation, differentiation, migration, apoptosis, and cell survival. Members of JAK family and STAT family have been extensively studied in different mammalian species because of their important roles in innate and adaptive immune responses. However, they have not been systematically studied among teleost fish species. In this study, five JAK family members and eight STAT family members were identified and characterized from channel catfish. Phylogenetic analysis was conducted to properly annotate these genes. Syntenic analysis was also conducted to establish orthology, and confirm the results from phylogenetic analysis. Compared to mammals, more members of the JAK and STAT family were identified in channel catfish genome. Expression of JAK and STAT family members was detected in healthy catfish tissues, but was induced in gill, liver, and intestine after bacterial challenge. Notably, the significant upregulation of STAT1b gene in catfish liver, gill and intestine after Edwardsiella ictaluri infection supported the notion that high STAT1 expression are involved in defense against pathogens. Collectively, the increased expression of JAK and STAT members in tested tissues suggested their crucial function in defending the host against pathogen invasion.

Cloning, characterisation, docking and expression analysis of 3-beta-hydroxysteroid dehydrogenase during ontogenetic development and annual reproductive cycles in catfish, Clarias batrachus

Fish like higher animals, have a well-defined mechanism to produce sex steroids that play a critical role in gonadal development and maturation. In this study, we aimed to analyse the expression pattern of 3β-HSD in different tissues, during ontogenetic development and gonadal recrudescence of Clarias batrachus. A full-length cDNA of 1617 bp including an open reading frame (ORF) of 1125 bp encoding 374 amino acids was isolated from testes of C. batrachus. The docking analysis between C. batrachus 3β-HSD protein and eurycomanone exhibited high binding affinity toward each other with total energy of -108.292 kcal/mol and van der Waals (VDW) interaction of -84.2838 kcal/mol. The 3β-HSD transcript level during ontogeny was detected in all the stages starting from the fertilized egg. The mature C. batrachus showed more expression of 3β-HSD mRNA in gonads and brain while weak expression was detected in the remaining tissues analysed. The 3β-HSD mRNA expression during annual reproductive phases of gonads was more in preparatory and pre-spawning stages than that of spawning and post-spawning phases. The mRNA expression results together suggest that 3β-HSD plays an important role in gonadal development. Furthermore, the active binding sites on 3β-HSD protein could be targeted in pharmacological drug designing to cope with reproductive dysfunctions in fish.

The CC and CXC chemokine receptors in channel catfish (Ictalurus punctatus) and their involvement in disease and hypoxia responses

Chemokines are vital regulators of cell mobilization for immune surveillance, inflammation, and development. Chemokines signal through binding to their receptors that are a superfamily of seven-transmembrane domain G-coupled receptors. Recently, a complete repertoire of both CC and CXC chemokines have been identified in channel catfish, but nothing is known about their receptors. In this study, a set of 29 CC chemokine receptor (CCR) genes and 8 CXC chemokine receptor (CXCR) genes were identified and annotated from the channel catfish genome. Extensive phylogenetic and comparative genomic analyses were conducted to annotate these genes, revealing fish-specific CC chemokine receptors, and lineage-specific tandem duplications of chemokine receptors in the teleost genomes. With 29 genes, the channel catfish genome harbors the largest numbers of CC chemokine receptors among all the genomes characterized. Analysis of gene expression after bacterial infections indicated that the chemokine receptors were regulated in a gene-specific manner. Most differentially expressed chemokine receptors were up-regulated after Edwardsiella ictaluri and Flavobacterium columnare infection. Among which, CXCR3 and CXCR4 were observed to participate in immune responses to both bacterial infections, indicating their potential roles in catfish immune activities. In addition, CXCR3.2 was significantly up-regulated in ESC-susceptible fish, and CXCR4b was mildly induced in ESC-resistant fish, further supporting the significant roles of CXCR3 and CXCR4 in catfish immune responses. CXCR4b and CCR9a were both up-regulated not only after bacterial infection, but also after hypoxia stress, providing the linkage between bacterial infection and low oxygen stresses. These results should be valuable for comparative immunological studies and provide insights into their roles in disease and stress responses.

Selection shapes the patterns of codon usage in three closely related species of genus Misgurnus

Neutrality plots revealed that selection probably dominates codon bias, whereas mutation plays only a minor role, in shaping the codon bias in three loaches, Misgurnus anguillicaudatus, M. mohoity, and M. bipartitus. These three species also clearly showed similar tendencies in the preferential usage of codons. Nineteen, nine, and 14 preferred codon pairs and 179, 182, and 173 avoided codon pairs were also detected in M. anguillicaudatus, M. bipartitus, and M. mohoity, respectively, and the most frequently avoided type of cP3-cA1 dinucleotide in these species was nnUAnn. The expression-linked patterns of codon usage revealed that higher expression was associated with higher GC₃, lower ENC, and a smaller proportion of amino acids with high size/complexity (S/C) scores in these three species. These results elucidate selectively driven codon bias in Misgurnus species, and reveal the potential importance of expression-mediated selection in shaping the genome evolution of fish.

The NCK and ABI adaptor genes in catfish and their involvement in ESC disease response

Adaptor proteins non-catalytic region of tyrosine kinase (NCK) and Abelson interactor (ABI) are crucial for disease response. NCK1 was identified to be a candidate gene for enteric septicemia of catfish (ESC) disease resistance, and was speculated to play similar roles during ESC and enteropathogenic Escherichia coli (EPEC) pathogenicity. ABI1 was reported as a positional candidate gene for bacterial cold water disease (BCWD) resistance in rainbow trout. In this study, three NCK genes and six ABI genes were identified in the channel catfish (Ictalurus punctatus) genome and blue catfish (I. furcatus) transcriptome, and annotated by domain structures, phylogenetic and syntenic analyses. Their expression patterns were examined in the intestine and liver of catfish after challenge with Edwardsiella ictaluri. In the intestine, NCK1, ABI2a, ABI2b, ABI3a were differentially expressed after E. ictaluri infection. In the liver, NCK2a, NCK2b, ABI1b, ABI2a, ABI2b were significantly upregulated in ESC susceptible fish. In general, the NCK and ABI genes, with exception of ABI3a gene and NCK1 gene, were expressed at higher levels in susceptible fish after infection than in control fish, but were expressed at lower levels in resistant fish than in the control fish. Taken together, these results support the notion that NCK and ABI genes are involved in disease processes facilitating pathogenesis of the E. ictaluri bacteria.

The chemokinome superfamily: II. The 64 CC chemokines in channel catfish and their involvement in disease and hypoxia responses

Chemokines are a superfamily of structurally related chemotactic cytokines exerting significant roles in regulating cell migration and activation. Based on the arrangement of the first four cysteine residues, they are classified into CC, CXC, C and CX3C subfamilies. In this study, a complete set of 64 CC chemokine ligand (CCL) genes was systematically identified, annotated, and characterized from the channel catfish genome. Extensive phylogenetic and comparative genomic analyses supported their annotations, allowing establishment of their orthologies, revealing fish-specific CC chemokines and the expansion of CC chemokines in the teleost genomes through lineage-specific tandem duplications. With 64 genes, the channel catfish genome harbors the largest numbers of CC chemokines among all the genomes characterized to date, however, they fall into 11 distinct CC chemokine groups. Analysis of gene expression after bacterial infections indicated that the CC chemokines were regulated in a gene-specific and time-dependent manner. While only one member of CCL19 (CCL19a.1) was significantly up-regulated after Edwardsiella ictaluri infection, all CCL19 members (CCL19a.1, CCL19a.2 and CCL19b) were significantly induced after Flavobacterium columnare infection. In addition, CCL19a.1, CCL19a.2 and CCL19b were also drastically up-regulated in ESC-susceptible fish, but not in resistant fish, suggesting potential significant roles of CCL19 in catfish immune responses. High expression levels of certain CC appeared to be correlated with susceptibility to diseases and intolerance to hypoxia.

Analysis on the expression and function of a chicken-type and goose-type lysozymes in Chinese giant salamanders Andrias davidianus

Lysozymes as an important immune factor, play vital roles in innate immune response against pathogen infection. In the present study, one c-type and g-type lysozymes were identified from Chinese giant salamander (Andrias davidianus). They shared highly conserved structural features with lysozymes from other species. Spatial expression analysis revealed that AdlysC transcript was most abundant in liver and stomach, and least in muscle and brain. In contrast, the expression level of AdlysG was most abundant in liver and least in muscle and skin. The transcription level of c-type and g-type lysozymes were up-regulated after Aeromonas hydrophila infection in liver and spleen, indicating their participations in the immune response. Moreover, the recombinant AdlysC and AdlysG protein were produced and purified, and were used to investigate the lysozyme activity at different pH and temperatures. The optimal lytic activity was determined at pH 6.0 and at a temperature of 30 °C. Through the minimal inhibitory concentration test, the rAdlysC and rAdlysG exhibited apparent antibacterial activity against both Gram-positive and Gram-negative bacteria with a variable concentration. In conclusion, it is the first report of lysozymes in A. davidianus, and c-type and g-type lysozymes should be involved in the innate immune response of A. davidianus.

Transcriptome of Pterospermum kingtungense provides implications on the mechanism underlying its rapid vegetative growth and limestone adaption

Pterospermum kingtungense C.Y.Wu ex Hsue is a typical tree species living in the relatively adverse limestone habitat. Due to its excellent wood quality and big size, it is an important timber resource which caused its endangered. We firstly provide the data resources by reporting an annotated transcriptome assembly. 203 million unique Illumina RNA-seq reads were produced with totally 50,333 transcripts, among which 48,778 transcripts were annotated. By a global comparison of homology between P. kingtungense and cacao, we identified 9,507 single copy orthologues and 990 P. kingtungense specific genes. GO enrichment analyses indicate that P. kingtungense specific genes are enriched in defense response, implying potential adaptation to limestone environment. As to cell compartment, the genes are enriched in thylakoid component. Consistently, KEGG enrichment indicates that genes are enriched in photosynthesis. In addition, we identified two genes under positive selection in P. kingtungense species. These results suggest that P. kingtungense have strong photosynthetic capacity, which related to vegetation growth. Our work provides the genomic resources of a limestone specific tree with economic importance to local society and suggests possible mechanism on its characteristics on the limestone adaption and excellent wood properties, which will be important for its conservation and sustainable utilization.

Identification of NF-κB related genes in channel catfish and their expression profiles in mucosal tissues after columnaris bacterial infection

Interactions of NF-κB family, IκB family and IKK complex are the key components of NF-κB pathway that is essential for many biological processes including innate and adaptive immunity, inflammation and stress responses. In spite of their importance, systematic analysis of these genes in fish has been lacking. Here we report a systematic study of the NF-κB related genes in channel catfish. Five NF-κB family genes, five IκB family genes and three IKK complex genes were identified in the channel catfish genome. Annotation of these 13 NF-κB related genes was further confirmed by phylogenetic and syntenic analysis. Negative selection was found to play a crucial role in the adaptive evolution of these genes. Expression profiles of NF-κB related genes after Flavobacterium columnare (columnaris) infection were determined by analysis of the existing RNA-Seq dataset. The majority of NF-κB related genes were significantly regulated in mucosal tissues of gill, skin and intestine after columnaris infection, indicating their potential involvement in host defense responses. Distinct expression patterns of NF-κB related genes were observed in susceptible and resistant catfish in response to columnaris infection, suggesting that expression of these genes may contribute to the variations in disease resistance/susceptibility of catfish.

An isoform of Taiman that contains a PRD-repeat motif is indispensable for transducing the vitellogenic juvenile hormone signal in Locusta migratoria

Taiman (Tai) has been recently identified as the dimerizing partner of juvenile hormone (JH) receptor, Methoprene-tolerant (Met). However, the role of Tai isoforms in transducing vitellogenic signal of JH has not been determined. In this study, we show that the migratory locust Locusta migratoria has two Tai isoforms, which differ in an INDEL-1 domain with the PRD-repeat motif rich in histidine and proline at the C-terminus. Tai-A with the INDEL-1 is expressed at levels about 50-fold higher than Tai-B without the INDEL-1 in the fat body of vitellogenic adult females. Knockdown of Tai-A but not Tai-B results in a substantial reduction of vitellogenin expression in the fat body accompanied by the arrest of ovarian development and oocyte maturation, similar to that caused by depletion of both Tai isoforms. Either Tai-A or Tai-B combined with Met can induce target gene transcription in response to JH, but Tai-A appears to mediate a significantly higher transactivation. Our data suggest that the INDEL-1 domain plays a critical role in Tai function during reproduction as Tai-A appears be more active than Tai-B in transducing the vitellogenic JH signal in L. migratoria.

Aquaculture genomics, genetics and breeding in the United States: current status, challenges, and priorities for future research

Advancing the production efficiency and profitability of aquaculture is dependent upon the ability to utilize a diverse array of genetic resources. The ultimate goals of aquaculture genomics, genetics and breeding research are to enhance aquaculture production efficiency, sustainability, product quality, and profitability in support of the commercial sector and for the benefit of consumers. In order to achieve these goals, it is important to understand the genomic structure and organization of aquaculture species, and their genomic and phenomic variations, as well as the genetic basis of traits and their interrelationships. In addition, it is also important to understand the mechanisms of regulation and evolutionary conservation at the levels of genome, transcriptome, proteome, epigenome, and systems biology. With genomic information and information between the genomes and phenomes, technologies for marker/causal mutation-assisted selection, genome selection, and genome editing can be developed for applications in aquaculture. A set of genomic tools and resources must be made available including reference genome sequences and their annotations (including coding and non-coding regulatory elements), genome-wide polymorphic markers, efficient genotyping platforms, high-density and high-resolution linkage maps, and transcriptome resources including non-coding transcripts. Genomic and genetic control of important performance and production traits, such as disease resistance, feed conversion efficiency, growth rate, processing yield, behaviour, reproductive characteristics, and tolerance to environmental stressors like low dissolved oxygen, high or low water temperature and salinity, must be understood. QTL need to be identified, validated across strains, lines and populations, and their mechanisms of control understood. Causal gene(s) need to be identified. Genetic and epigenetic regulation of important aquaculture traits need to be determined, and technologies for marker-assisted selection, causal gene/mutation-assisted selection, genome selection, and genome editing using CRISPR and other technologies must be developed, demonstrated with applicability, and application to aquaculture industries.Major progress has been made in aquaculture genomics for dozens of fish and shellfish species including the development of genetic linkage maps, physical maps, microarrays, single nucleotide polymorphism (SNP) arrays, transcriptome databases and various stages of genome reference sequences. This paper provides a general review of the current status, challenges and future research needs of aquaculture genomics, genetics, and breeding, with a focus on major aquaculture species in the United States: catfish, rainbow trout, Atlantic salmon, tilapia, striped bass, oysters, and shrimp. While the overall research priorities and the practical goals are similar across various aquaculture species, the current status in each species should dictate the next priority areas within the species. This paper is an output of the USDA Workshop for Aquaculture Genomics, Genetics, and Breeding held in late March 2016 in Auburn, Alabama, with participants from all parts of the United States.

The liver transcriptome of suckermouth armoured catfish (Pterygoplichthys anisitsi, Loricariidae): Identification of expansions in defensome gene families

Pterygoplichthys is a genus of related suckermouth armoured catfishes native to South America, which have invaded tropical and subtropical regions worldwide. Physiological features, including an augmented resistance to organic xenobiotics, may have aided their settlement in foreign habitats. The liver transcriptome of Pterygoplichthys anisitsi was sequenced and used to characterize the diversity of mRNAs potentially involved in the responses to natural and anthropogenic chemicals. In total, 66,642 transcripts were assembled. Among the identified defensome genes, cytochromes P450 (CYP) were the most abundant, followed by sulfotransferases (SULT), nuclear receptors (NR) and ATP binding cassette transporters (ABC). A novel expansion in the CYP2Y subfamily was identified, as well as an independent expansion of the CYP2AAs. Two expansions were also observed among SULT1. Thirty-two transcripts were classified into twelve subfamilies of NR, while 21 encoded ABC transporters. The diversity of defensome transcripts sequenced herein could contribute to this species' resistance to organic xenobiotics.

Identification, annotation and expression analysis of 29 Rho GTPase genes from channel catfish (Ictalurus punctatus) after bacterial infections

The Rho family GTPases are a group of small monomeric G proteins, which are molecular switches in signaling pathways. They have been known to regulate a diverse range of cellular processes including actin cytoskeleton rearrangement and microtubule dynamics. In particular, their participations in immune responses are also significant. However, little information of the Rho GTPases is available in teleost including channel catfish, an economically important species and one of the best teleost models forimmunological research. In this study, Rho GTPase genes were identified from channel catfish and well annotated by phylogenetic and syntenic analyses. Their expression profiles were determined in channel catfish healthy tissues and infected tissues. Altogether seven Rho GTPase genes were significantly regulated after bacterial infection, with six genes in the gill after Flavobacterium columnare challenge and two genes in the intestine in response to Edwardsiella ictaluri. All the differentially expressed genes were up-regulated soon after bacterial infection. Different expression patterns between the two experiments were observed, which may be attributed to tissue-specific regulation or pathogen-specific regulation. These results suggested that Rho GTPases play important roles in immune responses to bacterial pathogens, setting a foundation for future investigation on Rho GTPases.

Development of a 690 K SNP array in catfish and its application for genetic mapping and validation of the reference genome sequence

Single nucleotide polymorphisms (SNPs) are capable of providing the highest level of genome coverage for genomic and genetic analysis because of their abundance and relatively even distribution in the genome. Such a capacity, however, cannot be achieved without an efficient genotyping platform such as SNP arrays. In this work, we developed a high-density SNP array with 690,662 unique SNPs (herein 690 K array) that were relatively evenly distributed across the entire genome, and covered 98.6% of the reference genome sequence. Here we also report linkage mapping using the 690 K array, which allowed mapping of over 250,000 SNPs on the linkage map, the highest marker density among all the constructed linkage maps. These markers were mapped to 29 linkage groups (LGs) with 30,591 unique marker positions. This linkage map anchored 1,602 scaffolds of the reference genome sequence to LGs, accounting for over 97% of the total genome assembly. A total of 1,007 previously unmapped scaffolds were placed to LGs, allowing validation and in few instances correction of the reference genome sequence assembly. This linkage map should serve as a valuable resource for various genetic and genomic analyses, especially for GWAS and QTL mapping for genes associated with economically important traits.

Expression of Bcl-2 genes in channel catfish after bacterial infection and hypoxia stress

Bcl-2 proteins are of vital importance in regulation of apoptosis, and are involved in a number of biological processes such as carcinogenesis and immune responses. Bcl-2 genes have been well studied in mammals, while they are not well investigated in teleost fish including channel catfish, the major aquaculture species in the United States. In this study, we identified 34 bcl-2 genes from the channel catfish genome, and verified their identities by conducting phylogenetic and syntenic analyses. The expression profiles of the bcl-2 genes in response to bacterial infections (Edwardsiella ictaluri and Flavobacterium columnare) and hypoxia stress were determined by performing meta-analysis using the existing RNA-Seq datasets. Differential expressions of bcl-2 genes were observed after bacterial infections and hypoxia treatment, including 22 bcl-2 genes after E. ictaluri infection, 22 bcl-2 genes after F. columnare infection, and 19 bcl-2 genes after hypoxia stress. Overall, the expression of the pro-apoptotic bcl-2 genes were repressed after bacterial infection and hypoxia stress, indicating that bcl-2 genes are potentially involved in the stress response by reducing cell apoptosis. Some bcl-2 genes, such as bcl2b, mcl1a, bmf1, and bnip3, showed different expression pattern during the E. ictaluri and F. columnare infection, suggesting the difference in the pathogenicity of diseases. This work presented the first systematic identification and annotation of bcl-2 genes in catfish, providing essential genomic resources for further immune and physiological studies.

De novo assembly of Sockeye salmon kidney transcriptomes reveal a limited early response to piscine reovirus with or without infectious hematopoietic necrosis virus superinfection

Background

Piscine reovirus (PRV) has been associated with the serious disease known as Heart and Skeletal Muscle Inflammation (HSMI) in cultured Atlantic salmon Salmo salar in Norway. PRV is also prevalent in wild and farmed salmon without overt disease manifestations, suggesting multifactorial triggers or PRV variant-specific factors are required to initiate disease. In this study, we explore the head kidney transcriptome of Sockeye salmon Oncorhynchus nerka during early PRV infection to identify host responses in the absence of disease in hopes of elucidating mechanisms by which PRV may directly alter host functions and contribute to the development of a disease state. We further investigate the role of PRV as a coinfecting agent following superinfection with infectious hematopoietic necrosis virus (IHNV) – a highly pathogenic rhabdovirus endemic to the west coast of North America.

Results

Challenge of Sockeye salmon with PRV resulted in high quantities of viral transcripts to become present in the blood and kidney of infected fish without manifestations of disease. De novo transcriptome assembly of over 2.3 billion paired RNA-seq reads from the head kidneys of 36 fish identified more than 320,000 putative unigenes, of which less than 20 were suggested to be differentially expressed in response to PRV at either 2 or 3 weeks post challenge by DESeq2 and edgeR analysis. Of these, only one, Ependymin, was confirmed to be differentially expressed by qPCR in an expanded sample set. In contrast, IHNV induced substantial transcriptional changes (differential expression of > 20,000 unigenes) which included transcripts involved in antiviral and inflammatory response pathways. Prior infection with PRV had no significant effect on host responses to superinfecting IHNV, nor did host responses initiated by IHNV exposure influence increasing PRV loads.

Conclusions

PRV does not substantially alter the head kidney transcriptome of Sockeye salmon during early (2 to 3 week) infection and dissemination in a period of significant increasing viral load, nor does the presence of PRV change the host transcriptional response to an IHNV superinfection. Further, concurrent infections of PRV and IHNV do not appear to significantly influence the infectivity or severity of IHNV associated disease, or conversely, PRV load.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3196-y) contains supplementary material, which is available to authorized users.

A transcriptomic scan for potential candidate genes involved in osmoregulation in an obligate freshwater palaemonid prawn (Macrobrachium australiense)

Background

Understanding the genomic basis of osmoregulation (candidate genes and/or molecular mechanisms controlling the phenotype) addresses one of the fundamental questions in evolutionary ecology. Species distributions and adaptive radiations are thought to be controlled by environmental salinity levels, and efficient osmoregulatory (ionic balance) ability is the main mechanism to overcome the problems related to environmental salinity gradients.

Methods

To better understand how osmoregulatory performance in freshwater (FW) crustaceans allow individuals to acclimate and adapt to raised salinity conditions, here we (i), reviewed the literature on genes that have been identified to be associated with osmoregulation in FW crustaceans, and (ii), performed a transcriptomic analysis using cDNA libraries developed from mRNA isolated from three important osmoregulatory tissues (gill, antennal gland, hepatopancreas) and total mRNA from post larvae taken from the freshwater prawn, Macrobrachium australiense using Illumina deep sequencing technology. This species was targeted because it can complete its life cycle totally in freshwater but, like many Macrobrachium sp., can also tolerate brackish water conditions and hence should have genes associated with tolerance of both FW and saline conditions.

Results

We obtained between 55.4 and 65.2 million Illumina read pairs from four cDNA libraries. Overall, paired end sequences assembled into a total of 125,196 non-redundant contigs (≥200 bp) with an N50 length of 2,282 bp and an average contig length of 968 bp. Transcriptomic analysis of M. australiense identified 32 different gene families that were potentially involved with osmoregulatory capacity. A total of 32,597 transcripts were specified with gene ontology (GO) terms identified on the basis of GO categories. Abundance estimation of expressed genes based on TPM (transcript per million) ≥20 showed 1625 transcripts commonly expressed in all four libraries. Among the top 10 genes expressed in four tissue libraries associated with osmoregulation, arginine kinase and Na+/K+- ATPase showed the highest transcript copy number with 7098 and 660, respectively in gill which is considered to be the most important organ involved in osmoregulation.

Discussion

The current study provides the first broad transcriptome from M. australiense using next generation sequencing and identifies potential candidate genes involved in salinity tolerance and osmoregulation that can provide a foundation for investigating osmoregulatory capacity in a wide variety of freshwater crustaceans.

An 'Omics' Perspective on Cardiomyopathies and Heart Failure

Pathological enlargement of the heart, represented by hypertrophic cardiomyopathy (HCM) and dilated cardiomyopathy (DCM), occurs in response to many genetic and non-genetic factors. The clinical course of cardiac hypertrophy is remarkably variable, ranging from lifelong absence of symptoms to rapidly declining heart function and sudden cardiac death (SCD). Unbiased omics studies have begun to provide a glimpse into the molecular framework underpinning altered mechanotransduction, mitochondrial energetics, oxidative stress, and extracellular matrix in the heart undergoing physiological and pathological hypertrophy. Omics analyses indicate that post-transcriptional regulation of gene expression plays an overriding role in the normal and diseased heart. Studies to date highlight a need for more effective bioinformatics to better integrate patient omics data with their comprehensive clinical histories.

Characterization of common carp (Cyprinus carpio L.) interferon regulatory factor 5 (IRF5) and its expression in response to viral and bacterial challenges

Background

Common carp (Cyprinus carpio L.), one of the most economically valuable commercial farming fish species in China, is often infected by a variety of viruses. As the first line of defence against microbial pathogens, the innate immune system plays a crucial role in teleost fish, which are lower vertebrates. Interferon (IFN) regulatory factor 5 (IRF5) is a key molecule in antiviral immunity that regulating the expression of IFN and other pro-inflammatory cytokines. It is necessary to gain more insight into the common carp IFN system and the function of fish IRF5 in the antiviral and antibacterial response.

Results

In the present study, we characterized the cDNA and genomic sequence of the IRF5 gene in common carp, and analysed tissue distribution and expression profile of this gene in response to polyinosinic:polycytidylic acid (poly I:C) and lipopolysaccharides (LPS) treatment. The common carp IRF5 (ccIRF5) gene is 5790 bp in length and is composed of 9 exons and 8 introns. The open reading frame (ORF) of ccIRF5 is 1554 bp, and encodes 517 amino acid protein. The putative ccIRF5 protein shares identity (65.4–90.0 %) with other fish IRF5s and contains a DNA binding domain (DBD), a middle region (MR), an IRF-associated domain (IAD), a virus activated domain (VAD) and two nuclear localization signals (NLSs) similar to those found in vertebrate IRF5. Phylogenetic analysis clustered ccIRF5 into the IRF5 subfamily with other vertebrate IRF5 and IRF6 genes. Real-time PCR analysis revealed that ccIRF5 mRNA was expressed in all examined tissues of healthy carps, with high levels observed in the gills and the brain. After poly I:C challenge, expression levels of ccIRF5, tumour-necrosis factor α (ccTNFα) and two IFN stimulated genes [ISGs (ccISG5 and ccPKR)] were up-regulated in seven immune-related tissues (liver, spleen, head kidney, foregut, hindgut, skin and gills). Furthermore, all four genes were up-regulated in vitro upon poly I:C and LPS challenges.

Conclusions

Our findings suggest that IRF5 might play an important role in regulating the antiviral and antibacterial response in fish. These results could provide a clue for preventing common carp infection by pathogenic microorganisms present in the aquatic environment.

Electronic supplementary material

The online version of this article (doi:10.1186/s12917-016-0750-4) contains supplementary material, which is available to authorized users.

The channel catfish genome sequence provides insights into the evolution of scale formation in teleosts

Catfish represent 12% of teleost or 6.3% of all vertebrate species, and are of enormous economic value. Here we report a high-quality reference genome sequence of channel catfish (Ictalurus punctatus), the major aquaculture species in the US. The reference genome sequence was validated by genetic mapping of 54,000 SNPs, and annotated with 26,661 predicted protein-coding genes. Through comparative analysis of genomes and transcriptomes of scaled and scaleless fish and scale regeneration experiments, we address the genomic basis for the most striking physical characteristic of catfish, the evolutionary loss of scales and provide evidence that lack of secretory calcium-binding phosphoproteins accounts for the evolutionary loss of scales in catfish. The channel catfish reference genome sequence, along with two additional genome sequences and transcriptomes of scaled catfishes, provide crucial resources for evolutionary and biological studies. This work also demonstrates the power of comparative subtraction of candidate genes for traits of structural significance.

Optimizing Hybrid de Novo Transcriptome Assembly and Extending Genomic Resources for Giant Freshwater Prawns (Macrobrachium rosenbergii): The Identification of Genes and Markers Associated with Reproduction

The giant freshwater prawn, Macrobrachium rosenbergii, a sexually dimorphic decapod crustacean is currently the world’s most economically important cultured freshwater crustacean species. Despite its economic importance, there is currently a lack of genomic resources available for this species, and this has limited exploration of the molecular mechanisms that control the M. rosenbergii sex-differentiation system more widely in freshwater prawns. Here, we present the first hybrid transcriptome from M. rosenbergii applying RNA-Seq technologies directed at identifying genes that have potential functional roles in reproductive-related traits. A total of 13,733,210 combined raw reads (1720 Mbp) were obtained from Ion-Torrent PGM and 454 FLX. Bioinformatic analyses based on three state-of-the-art assemblers, the CLC Genomic Workbench, Trans-ABySS, and Trinity, that use single and multiple k-mer methods respectively, were used to analyse the data. The influence of multiple k-mers on assembly performance was assessed to gain insight into transcriptome assembly from short reads. After optimisation, de novo assembly resulted in 44,407 contigs with a mean length of 437 bp, and the assembled transcripts were further functionally annotated to detect single nucleotide polymorphisms and simple sequence repeat motifs. Gene expression analysis was also used to compare expression patterns from ovary and testis tissue libraries to identify genes with potential roles in reproduction and sex differentiation. The large transcript set assembled here represents the most comprehensive set of transcriptomic resources ever developed for reproduction traits in M. rosenbergii, and the large number of genetic markers predicted should constitute an invaluable resource for future genetic research studies on M. rosenbergii and can be applied more widely on other freshwater prawn species in the genus Macrobrachium.

A Leukocyte Immune-Type Receptor Subset Is a Marker of Antiviral Cytotoxic Cells in Channel Catfish, Ictalurus punctatus

Channel catfish, Ictalurus punctatus, leukocyte immune type receptors (LITRs) represent a multigene family that encodes Ig superfamily proteins that mediate activating or inhibitory signaling. In this study, we demonstrate the use of mAb CC41 to monitor viral cytotoxic responses in catfish and determine that CC41 binds to a subset of LITRs on the surface of catfish clonal CTLs. Homozygous gynogenetic catfish were immunized with channel catfish virus (CCV)-infected MHC-matched clonal T cells (G14D-CCV), and PBL were collected at various times after immunization for flow cytometric analyses. The percentage of CC41(+) cells was significantly increased 5 d after primary immunization with G14D-CCV and at 3 d after a booster immunization as compared with control fish only injected with G14D. Moreover, CC41(+) cells magnetically isolated from the PBL specifically killed CCV-infected targets as measured by (51)Cr release assays and expressed messages for CD3γδ, perforin, and at least one of the CD4-like receptors as analyzed by RNA flow cytometry. When MLC effector cells derived from a G14D-CCV-immunized fish were preincubated with CC41 mAb, killing of G14D-CCV targets was reduced by ∼40%, suggesting that at least some LITRs have a role in target cell recognition and/or cytotoxicity. The availability of a LITR-specific mAb has allowed, to our knowledge for the first time, functional characterization of LITRs in an autologous system. In addition, the identification of an LITR subset as a cytotoxic cell marker will allow for more effective monitoring of catfish immune responses to pathogens.

Phosphoinositide 3-kinase family in channel catfish and their regulated expression after bacterial infection

The phosphoinositide-3-kinase (PI3Ks) family of lipid kinases is widely conserved from yeast to mammals. In this work, we identified a total of 14 members of the PI3Ks from the channel catfish genome and transcriptome and conducted phylogenetic and syntenic analyses of these genes. The expression profiles after infection with Edwardsiella ictaluri and Flavobacterium columnare were examined to determine the involvement of PI3Ks in immune responses after bacterial infection in catfish. The results indicated that PI3Ks genes including all of the catalytic subunit and several regulatory subunits genes were widely regulated after bacterial infection. The expression patterns were quite different when challenged with different bacteria. The PI3Ks were up-regulated rapidly at the early stage after ESC infection, but their induced expression was much slower, at the middle stage after columnaris infection. RNA-Seq datasets indicated that PI3K genes may be expressed at different levels in different catfish differing in their resistance levels against columnaris. Future studies are required to confirm and validate these observations. Taken together, this study indicated that PI3K genes may be involved as a part of the defense responses of catfish after infections, and they could be one of the determinants for disease resistance.

Genome-wide identification of Hsp70 genes in channel catfish and their regulated expression after bacterial infection

Heat shock proteins 70/110 (Hsp70/110) are a family of conserved ubiquitously expressed heat shock proteins which are produced by cells in response to exposure to stressful conditions. Besides the chaperone and housekeeping functions, they are also known to be involved in immune response during infection. In this study, we identified 16 Hsp70/110 geness in channel catfish (Ictalurus punctatus) through in silico analysis using RNA-Seq and genome databases. Among them 12 members of Hsp70 (Hspa) family and 4 members of Hsp110 (Hsph) family were identified. Phylogenetic and syntenic analyses provided strong evidence in supporting the orthologies of these HSPs. In addition, we also determined the expression patterns of Hsp70/110 genes after Flavobacterium columnare and Edwardsiella ictaluri infections by meta-analyses, for the first time in channel catfish. Ten out of sixteen genes were significantly up/down-regulated after bacterial challenges. Specifically, nine genes were found significantly expressed in gill after F. columnare infection. Two genes were found significantly expressed in intestine after E. ictaluri infection. Pathogen-specific pattern and tissue-specific pattern were found in the two infections. The significantly regulated expressions of catfish Hsp70 genes after bacterial infections suggested their involvement in immune response in catfish.

Transcriptome analysis of the Larimichthys crocea liver in response to Cryptocaryon irritans

The large yellow croaker (Larimichthys crocea) is an economically important marine fish cultured in China and East Asian countries and is facing a serious threat from Cryptocaryon irritans, which is a protozoan ectoparasite that infects most reared marine fish species. To understand the molecular immune mechanisms underlying the response to C. irritans, we first performed a comparative gene transcription analysis using livers from C. irritans-immunized L. croceas and from a control group through RNA-Seq technology. After the removal of low-quality sequences and assembly, 51360 contigs were obtained, with an average length of 1066.93 bp. Further, a blast analysis indicates that 30747 contigs can be annotated based on homology with matches in the NT, NR, gene, and string databases. A gene ontology analysis was used to classify 21598 genes according to three major functional categories: molecular function, cellular component, and biological process. Moreover, 14470 genes were found in 303 KEGG pathways. We used RSEM and EdgeR to determine that 3841 genes were significantly differentially expressed (FDR < 0.001), including 2129 up-regulated genes and 1712 down-regulated genes. A significant enrichment analysis of these differentially expressed genes and isogenes revealed major immune-related pathways, including the toll-like receptor, complement and coagulation cascades, and chemokine signaling pathways. In addition, 28748 potential simple sequence repeats (SSRs) were detected from 12776 transcripts, and 62992 candidate single nucleotide polymorphisms (SNPs) were identified in the L. croceas liver transcriptome. This study characterized a gene expression pattern for normal and C. irritans-immunized L. croceas for the first time and not only sheds new light on the molecular mechanisms underlying the host-C. irritans interaction but also facilitates future studies on L. croceas gene expression and functional genomics.

Large-scale transcriptome sequencing reveals novel expression patterns for key sex-related genes in a sex-changing fish

Background

Teleost fishes exhibit remarkably diverse and plastic sexual developmental patterns. One of the most astonishing is the rapid socially controlled female-to-male (protogynous) sex change observed in bluehead wrasses (Thalassoma bifasciatum). Such functional sex change is widespread in marine fishes, including species of commercial importance, yet its underlying molecular basis remains poorly explored.

Methods

RNA sequencing was performed to characterize the transcriptomic profiles and identify genes exhibiting sex-biased expression in the brain (forebrain and midbrain) and gonads of bluehead wrasses. Functional annotation and enrichment analysis were carried out for the sex-biased genes in the gonad to detect global differences in gene products and genetic pathways between males and females.

Results

Here we report the first transcriptomic analysis for a protogynous fish. Expression comparison between males and females reveals a large set of genes with sex-biased expression in the gonad, but relatively few such sex-biased genes in the brain. Functional annotation and enrichment analysis suggested that ovaries are mainly enriched for metabolic processes and testes for signal transduction, particularly receptors of neurotransmitters and steroid hormones. When compared to other species, many genes previously implicated in male sex determination and differentiation pathways showed conservation in their gonadal expression patterns in bluehead wrasses. However, some critical female-pathway genes (e.g., rspo1 and wnt4b) exhibited unanticipated expression patterns. In the brain, gene expression patterns suggest that local neurosteroid production and signaling likely contribute to the sex differences observed.

Conclusions

Expression patterns of key sex-related genes suggest that sex-changing fish predominantly use an evolutionarily conserved genetic toolkit, but that subtle variability in the standard sex-determination regulatory network likely contributes to sexual plasticity in these fish. This study not only provides the first molecular data on a system ideally suited to explore the molecular basis of sexual plasticity and tissue re-engineering, but also sheds some light on the evolution of diverse sex determination and differentiation systems.

Electronic supplementary material

The online version of this article (doi:10.1186/s13293-015-0044-8) contains supplementary material, which is available to authorized users.

Complement regulatory protein genes in channel catfish and their involvement in disease defense response

Complement system is one of the most important defense systems of innate immunity, which plays a crucial role in disease defense responses in channel catfish. However, inappropriate and excessive complement activation could lead to potential damage to the host cells. Therefore the complement system is controlled by a set of complement regulatory proteins to allow normal defensive functions, but prevent hazardous complement activation to host tissues. In this study, we identified nine complement regulatory protein genes from the channel catfish genome. Phylogenetic and syntenic analyses were conducted to determine their orthology relationships, supporting their correct annotation and potential functional inferences. The expression profiles of the complement regulatory protein genes were determined in channel catfish healthy tissues and after infection with the two main bacterial pathogens, Edwardsiella ictaluri and Flavobacterium columnare. The vast majority of complement regulatory protein genes were significantly regulated after bacterial infections, but interestingly were generally up-regulated after E. ictaluri infection while mostly down-regulated after F. columnare infection, suggesting a pathogen-specific pattern of regulation. Collectively, these findings suggested that complement regulatory protein genes may play complex roles in the host immune responses to bacterial pathogens in channel catfish.

Effective de novo assembly of fish genome using haploid larvae

Recent improvements in next-generation sequencing technology have made it possible to do whole genome sequencing, on even non-model eukaryote species with no available reference genomes. However, de novo assembly of diploid genomes is still a big challenge because of allelic variation. The aim of this study was to determine the feasibility of utilizing the genome of haploid fish larvae for de novo assembly of whole-genome sequences. We compared the efficiency of assembly using the haploid genome of yellowtail (Seriola quinqueradiata) with that using the diploid genome obtained from the dam. De novo assembly from the haploid and the diploid sequence reads (100 million reads per each datasets) generated by the Ion Proton sequencer (200 bp) was done under two different assembly algorithms, namely overlap-layout-consensus (OLC) and de Bruijn graph (DBG). This revealed that the assembly of the haploid genome significantly reduced (approximately 22% for OLC, 9% for DBG) the total number of contigs (with longer average and N50 contig lengths) when compared to the diploid genome assembly. The haploid assembly also improved the quality of the scaffolds by reducing the number of regions with unassigned nucleotides (Ns) (total length of Ns; 45,331,916 bp for haploids and 67,724,360 bp for diploids) in OLC-based assemblies. It appears clear that the haploid genome assembly is better because the allelic variation in the diploid genome disrupts the extension of contigs during the assembly process. Our results indicate that utilizing the genome of haploid larvae leads to a significant improvement in the de novo assembly process, thus providing a novel strategy for the construction of reference genomes from non-model diploid organisms such as fish.