Revelation of candidate genes and molecular mechanism of reproductive seasonality in female rohu (Labeo rohita Ham.) by RNA sequencing.

BACKGROUND: Carp fish, rohu (Labeo rohita Ham.) is important freshwater aquaculture species of South-East Asia having seasonal reproductive rhythm. There is no holistic study at transcriptome level revealing key candidate genes involved in such circannual rhythm regulated by biological clock genes (BCGs). Seasonality manifestation has two contrasting phases of reproduction, i.e., post-spawning resting and initiation of gonadal activity appropriate for revealing the associated candidate genes. It can be deciphered by RNA sequencing of tissues involved in BPGL (Brain-Pituitary-Gonad-Liver) axis controlling seasonality. How far such BCGs of this fish are evolutionarily conserved across different phyla is unknown. Such study can be of further use to enhance fish productivity as seasonality restricts seed production beyond monsoon season. RESULT: A total of ~ 150 Gb of transcriptomic data of four tissues viz., BPGL were generated using Illumina TruSeq. De-novo assembled BPGL tissues revealed 75,554 differentially expressed transcripts, 115,534 SSRs, 65,584 SNPs, 514 pathways, 5379 transcription factors, 187 mature miRNA which regulates candidate genes represented by 1576 differentially expressed transcripts are available in the form of web-genomic resources. Findings were validated by qPCR. This is the first report in carp fish having 32 BCGs, found widely conserved in fish, amphibian, reptile, birds, prototheria, marsupials and placental mammals. This is due to universal mechanism of rhythmicity in response to environment and earth rotation having adaptive and reproductive significance. CONCLUSION: This study elucidates evolutionary conserved mechanism of photo-periodism sensing, neuroendocrine secretion, metabolism and yolk synthesis in liver, gonadal maturation, muscular growth with sensory and auditory perception in this fish. Study reveals fish as a good model for research on biological clock besides its relevance in reproductive efficiency enhancement.

Revealing liver specific microRNAs linked with carbohydrate metabolism of farmed carp, Labeo rohita (Hamilton, 1822).

The role of microRNA in gene regulation during developmental biology has been well depicted in several organisms. The present study was performed to investigate miRNAs role in the liver tissues during carbohydrate metabolism and their targets in the farmed carp rohu, Labeo rohita, which is economically important species in aquaculture. Using Illumina-HiSeq technology, a total of 22,612,316; 44,316,046 and 13,338,434 clean reads were obtained from three small-RNA libraries. We have identified 138 conserved and 161 novel miRNAs and studies revealed that miR-22, miR-122, miR-365, miR-200, and miR-146 are involved in carbohydrate metabolism. Further analysis depicted mature miRNA and their predicted target sites in genes that were involved in developmental biology, cellular activities, transportation, etc. This is the first report of the presence of miRNAs in liver tissue of rohu and their comparative profile linked with metabolism serves as a vital resource as a biomarker.

Revealing Alteration in the Hepatic Glucose Metabolism of Genetically Improved Carp, Jayanti Rohu Labeo rohita Fed a High Carbohydrate Diet Using Transcriptome Sequencing.

Although feed cost is the greatest concern in aquaculture, the inclusion of carbohydrates in the fish diet, and their assimilation, are still not well understood in aquaculture species. We identified molecular events that occur due to the inclusion of high carbohydrate levels in the diets of genetically improved 'Jayanti rohu' Labeo rohita. To reveal transcriptional changes in the liver of rohu, a feeding experiment was conducted with three doses of gelatinized starch (20% (control), 40%, and 60%). Transcriptome sequencing revealed totals of 15,232 (4464 up- and 4343 down-regulated) and 15,360 (4478 up- and 4171 down-regulated) differentially expressed genes. Up-regulated transcripts associated with glucose metabolisms, such as hexokinase, PHK, glycogen synthase and PGK, were found in fish fed diets with high starch levels. Interestingly, a de novo lipogenesis mechanism was found to be enriched in the livers of treated fish due to up-regulated transcripts such as FAS, ACCα, and PPARγ. The insulin signaling pathways with enriched PPAR and mTOR were identified by Kyoto Encyclopedia of Genes and Genome (KEGG) as a result of high carbohydrates. This work revealed for the first time the atypical regulation transcripts associated with glucose metabolism and lipogenesis in the livers of Jayanti rohu due to the inclusion of high carbohydrate levels in the diet. This study also encourages the exploration of early nutritional programming for enhancing glucose efficiency in carp species, for sustainable and cost-effective aquaculture production.

RNA-seq analysis of mucosal immune responses reveals signatures of intestinal barrier disruption and pathogen entry following Edwardsiella ictaluri infection in channel catfish, Ictalurus punctatus.

The mucosal surfaces of fish (gill, skin, gastrointestinal tract) are important sites of bacterial exposure and host defense mechanisms. In mammalian systems, the intestinal epithelium is well characterized as both a selectively permeable barrier regulated by junctional proteins and as a primary site of infection for a number of enteric pathogens including viruses, bacteria, and parasites. The causative bacterium of enteric septicemia of catfish, Edwardsiella ictaluri, is believed to gain entry through the intestinal epithelium, with previous research using a rat intestinal epithelial cell line (IEC-6) indicating actin polymerization and receptor-mediated endocytosis as potential mechanisms of uptake. Here, we utilized high-throughput RNA-seq to characterize the role of the intestinal epithelial barrier following E. ictaluri challenge. A total of 197.6 million reads were obtained and assembled into 176,481 contigs with an average length of 893.7 bp and N50 of 1676 bp. The assembled contigs contained 14,457 known unigenes, including 2719 genes not previously identified in other catfish transcriptome studies. Comparison of digital gene expression between challenged and control samples revealed 1633 differentially expressed genes at 3 h, 24 h, and 3 day following exposure. Gene pathway analysis of the differentially expressed gene set indicated the centrality of actin cytoskeletal polymerization/remodelling and junctional regulation in pathogen entry and subsequent inflammatory responses. The expression patterns of fifteen differentially expressed genes related to intestinal epithelial barrier dysfunction were validated by quantitative real-time RT-PCR (average correlation coeff. 0.92, p < 0.001). Our results set a foundation for future studies comparing mechanisms of pathogen entry and mucosal immunity across several important catfish pathogens including E. ictaluri, Edwardsiellatarda, Flavobacterium columnare, and virulent atypical Aeromonas hydrophila. Understanding of molecular mechanisms of pathogen entry during infection will provide insight into strategies for selection of resistant catfish brood stocks against various diseases.

Microarray analysis of gene expression in the blue catfish liver reveals early activation of the MHC class I pathway after infection with Edwardsiella ictaluri.

The acute nature of disease outbreaks in aquaculture settings has served to emphasize the importance of the innate immune response of fish for survival and led to the recent identification and characterization of many of its components. Catfish, the predominant aquaculture species in the United States, is an important model for the study of the teleost immune system. However, transcriptomic-level studies of disease-related gene expression in catfish have only recently been initiated, and understanding of immune responses to pathogen infections is limited. Here, we have developed and utilized a 28K in situ oligonucleotide microarray composed of blue catfish (Ictalurus furcatus) and channel catfish (Ictalurus punctatus) transcripts. While channel catfish accounts for the majority of commercial production, the closely related blue catfish possesses several economically important phenotypic traits. Microarray analysis of gene expression changes in blue catfish liver after infection with Gram-negative bacterium Edwardsiella ictaluri indicated the strong upregulation of several pathways involved in the inflammatory immune response and potentially in innate disease resistance. A multifaceted response to infection could be observed, encompassing the complement cascade, iron regulation, inflammatory cell signaling, and antigen processing and presentation. The induction of several components of the MHC class I-related pathway following infection with an intracellular bacterium is reported here for the first time in fish. A comparison with previously published expression profiles in the channel catfish liver was also made and the microarray results extended by use of quantitative RT-PCR. Our results add to the understanding of the teleost immune responses and provide a solid foundation for future functional characterization, genetic mapping, and QTL analysis of immunity-related genes from catfish.

Liver-Specific microRNA Identification in Farmed Carp, Labeo bata (Hamilton, 1822), Fed with Starch Diet Using High-Throughput Sequencing.

The liver is an important central organ, which controls carbohydrate metabolism through maintaining glucose homeostasis by a tightly regulated system of genes or enzymes. The microRNAs are small non-coding RNAs playing an important role in the regulation of genes associated with developmental biology, physiology, metabolism, etc. Thus, in this study, we have intended to detect liver-specific microRNAs in farmed carp, Labeo bata, upon being fed a diet with different levels of carbohydrates. Here, we have conducted the experiment for 45 days using fingerlings of farmed carp fed with 20% (control), 40%, and 60% gelatinized starch levels. The liver tissues were collected from each treatment and processed for RNA isolation, small RNA library preparation, and high-throughput sequencing using Illumina NexSeq500. Through sequencing, 15,779,417 reads in 20% CHO, 13,959,039 in 40% CHO, and 13,661,950 in 60% CHO reads were generated for control and treated fishes using three small RNA libraries. We have investigated 445 novel and 231 conserved microRNAs in 20%, 40%, and 60% carbohydrate (CHO), respectively, through computational analysis. The differential expression analysis of miRNAs was carried out between different treatments compared with control and this study depicted 117 known and 114 novel miRNA genes involved in carbohydrate metabolic pathways. Further, target prediction and gene ontology analysis revealed that miRNAs were involved in several pathways such as signaling pathway, G protein pathway, complement receptor-mediated pathway, dopamine receptor signaling pathway, epidermal growth factor pathway, and notch signaling pathway. The predicted miRNA sites in targeted genes were associated with cellular activities, developmental biology, DNA binding, Golgi apparatus, extracellular region, catalytic activity, MAPK cascade, etc. Overall, we have generated a vital resource of liver-specific miRNAs involved in metabolic gene regulation. These studies further will help develop miRNA inhibitors to study their role during carbohydrate metabolism in farmed carp.

Towards the ictalurid catfish transcriptome: generation and analysis of 31,215 catfish ESTs.

BACKGROUND: EST sequencing is one of the most efficient means for gene discovery and molecular marker development, and can be additionally utilized in both comparative genome analysis and evaluation of gene duplications. While much progress has been made in catfish genomics, large-scale EST resources have been lacking. The objectives of this project were to construct primary cDNA libraries, to conduct initial EST sequencing to generate catfish EST resources, and to obtain baseline information about highly expressed genes in various catfish organs to provide a guide for the production of normalized and subtracted cDNA libraries for large-scale transcriptome analysis in catfish. RESULTS: A total of 17 cDNA libraries were constructed including 12 from channel catfish (Ictalurus punctatus) and 5 from blue catfish (I. furcatus). A total of 31,215 ESTs, with average length of 778 bp, were generated including 20,451 from the channel catfish and 10,764 from blue catfish. Cluster analysis indicated that 73% of channel catfish and 67% of blue catfish ESTs were unique within the project. Over 53% and 50% of the channel catfish and blue catfish ESTs, respectively, had significant similarities to known genes. All ESTs have been deposited in GenBank. Evaluation of the catfish EST resources demonstrated their potential for molecular marker development, comparative genome analysis, and evaluation of ancient and recent gene duplications. Subtraction of abundantly expressed genes in a variety of catfish tissues, identified here, will allow the production of low-redundancy libraries for in-depth sequencing. CONCLUSION: The sequencing of 31,215 ESTs from channel catfish and blue catfish has significantly increased the EST resources in catfish. The EST resources should provide the potential for microarray development, polymorphic marker identification, mapping, and comparative genome analysis.

Expression analysis of the acute phase response in channel catfish (Ictalurus punctatus) after infection with a Gram-negative bacterium.

The acute phase response (APR) is a set of metabolic and physiological reactions occurring in the host in response to tissue infection or injury and is a crucial component of the larger innate immune response. The APR is best characterized by dramatic changes in the concentration of a group of plasma proteins known as acute phase proteins (APPs) which are synthesized in the liver and function in a wide range of immunity-related activities. Utilizing a new high-density in situ oligonucleotide microarray, we have evaluated the APR in channel catfish liver following infection with Edwardsiella ictaluri, a bacterial pathogen that causes enteric septicemia of catfish. Our catfish microarray design (28K) builds upon a previous 19K channel catfish array by adding recently sequenced immune transcripts from channel catfish along with 7159 unique sequences from closely related blue catfish. The analysis of microarray results using a traditional 2-fold change in gene expression cutoff and a 10% false-discovery rate revealed a well-developed APR in catfish, with particularly high upregulation (>50-fold) of genes involved in iron homeostasis (i.e. intelectin, hemopexin, haptoglobin, ferritin, and transferrin). Other classical APP genes upregulated greater than 2-fold included coagulation factors, proteinase inhibitors, transport proteins, and complement components. Upregulation of the majority of the complement cascade was observed including the membrane attack complex components and complement inhibitors. A number of pathogen recognition receptors (PRRs) and chemokines were also differentially expressed in the liver following infection. Independent testing of a selection of differentially expressed genes with real-time RT-PCR confirmed microarray results.

Construction, De-Novo Assembly and Analysis of Transcriptome for Identification of Reproduction-Related Genes and Pathways from Rohu, Labeo rohita (Hamilton).

Rohu is a leading candidate species for freshwater aquaculture in South-East Asia. Unlike common carp the monsoon breeding habit of rohu restricts its seed production beyond season indicating strong genetic control over spawning. Genetic information is limited in this regard. The problem is exacerbated by the lack of genomic-resources. We identified 182 reproduction-related genes previously by Sanger-sequencing which were less to address the issue of seasonal spawning behaviour of this important carp. Therefore, the present work was taken up to generate transcriptome profile by mRNAseq. 16 GB, 72 bp paired end (PE) data was generated from the pooled-RNA of twelve-tissues from pre-spawning rohu using IlluminaGA-II-platform. There were 64.97 million high-quality reads producing 62,283 contigs and 88,612 numbers of transcripts using velvet and oases programs, respectively. Gene ontology annotation identified 940 reproduction-related genes consisting of 184 mainly associated with reproduction, 223 related to hormone-activity and receptor-binding, 178 receptor-activity and 355 embryonic-development related-proteins. The important reproduction-relevant pathways found in KEGG analysis were GnRH-signaling, oocyte-meiosis, steroid-biosynthesis, steroid-hormone biosynthesis, progesterone-mediated oocyte-maturation, retinol-metabolism, neuroactive-ligand-receptor interaction, neurotrophin-signaling and photo-transduction. Twenty nine simple sequence repeat containing sequences were also found out of which 12 repeat loci were polymorphic with mean expected-&-observed heterozygosity of 0.471 and 0.983 respectively. Quantitative RT-PCR analyses of 13-known and 6-unknown transcripts revealed differences in expression level between preparatory and post-spawning phase. These transcriptomic sequences have significantly increased the genetic-&-genomic resources for reproduction-research in Labeo rohita.

Identification of reproduction-related genes and SSR-markers through expressed sequence tags analysis of a monsoon breeding carp rohu, Labeo rohita (Hamilton).

Labeo rohita (Ham.) also called rohu is the most important freshwater aquaculture species on the Indian sub continent. Monsoon dependent breeding restricts its seed production beyond season indicating a strong genetic control about which very limited information is available. Additionally, few genomic resources are publicly available for this species. Here we sought to identify reproduction-relevant genes from normalized cDNA libraries of the brain-pituitary-gonad-liver (BPGL-axis) tissues of adult L. rohita collected during post preparatory phase. 6161 random clones sequenced (Sanger-based) from these libraries produced 4642 (75.34%) high-quality sequences. They were assembled into 3631 (78.22%) unique sequences composed of 709 contigs and 2922 singletons. A total of 182 unique sequences were found to be associated with reproduction-related genes, mainly under the GO term categories of reproduction, neuro-peptide hormone activity, hormone and receptor binding, receptor activity, signal transduction, embryonic development, cell-cell signaling, cell death and anti-apoptosis process. Several important reproduction-related genes reported here for the first time in L. rohita are zona pellucida sperm-binding protein 3, aquaporin-12, spermine oxidase, sperm associated antigen 7, testis expressed 261, progesterone receptor membrane component, Neuropeptide Y and Pro-opiomelanocortin. Quantitative RT-PCR-based analyses of 8 known and 8 unknown transcripts during preparatory and post-spawning phase showed increased expression level of most of the transcripts during preparatory phase (except Neuropeptide Y) in comparison to post-spawning phase indicating possible roles in initiation of gonad maturation. Expression of unknown transcripts was also found in prolific breeder common carp and tilapia, but levels of expression were much higher in seasonal breeder rohu. 3631 unique sequences contained 236 (6.49%) putative microsatellites with the AG (28.16%) repeat as the most frequent motif. Twenty loci showed polymorphism in 36 unrelated individuals with allele frequency ranging from 2 to 7 per locus. The observed heterozygosity ranged from 0.096 to 0.774 whereas the expected heterozygosity ranged from 0.109 to 0.801. Identification of 182 important reproduction-related genes and expression pattern of 16 transcripts in preparatory and post-spawning phase along with 20 polymorphic EST-SSRs should be highly useful for the future reproductive molecular studies and selection program in Labeo rohita.

Assembly of 500,000 inter-specific catfish expressed sequence tags and large scale gene-associated marker development for whole genome association studies.

BACKGROUND: Through the Community Sequencing Program, a catfish EST sequencing project was carried out through a collaboration between the catfish research community and the Department of Energy's Joint Genome Institute. Prior to this project, only a limited EST resource from catfish was available for the purpose of SNP identification. RESULTS: A total of 438,321 quality ESTs were generated from 8 channel catfish (Ictalurus punctatus) and 4 blue catfish (Ictalurus furcatus) libraries, bringing the number of catfish ESTs to nearly 500,000. Assembly of all catfish ESTs resulted in 45,306 contigs and 66,272 singletons. Over 35% of the unique sequences had significant similarities to known genes, allowing the identification of 14,776 unique genes in catfish. Over 300,000 putative SNPs have been identified, of which approximately 48,000 are high-quality SNPs identified from contigs with at least four sequences and the minor allele presence of at least two sequences in the contig. The EST resource should be valuable for identification of microsatellites, genome annotation, large-scale expression analysis, and comparative genome analysis. CONCLUSIONS: This project generated a large EST resource for catfish that captured the majority of the catfish transcriptome. The parallel analysis of ESTs from two closely related Ictalurid catfishes should also provide powerful means for the evaluation of ancient and recent gene duplications, and for the development of high-density microarrays in catfish. The inter- and intra-specific SNPs identified from all catfish EST dataset assembly will greatly benefit the catfish introgression breeding program and whole genome association studies.

Repeat structure of the catfish genome: a genomic and transcriptomic assessment of Tc1-like transposon elements in channel catfish (Ictalurus punctatus).

We have assessed the distribution and diversity of members of the Tc1/mariner superfamily of transposable elements in the channel catfish (Ictalurus punctatus) genome as well as evaluating the extent of transcription of Tc1 transposases in the species. Through use of PCR amplification and sequencing, assessment of random BAC end sequences (BES) equivalent to 1.2% genome coverage, and screening of over 45,000 catfish ESTs, a significant proportion of Tc1-like elements and their associated transcripts were captured. Up to 4.2% of the catfish genome in base pairs appears to be composed of Tc1-like transposon-related sequences and a significant fraction of the catfish cellular mRNA, approximately 0.6%, was transcribed from transposon-related sequences in both sense and antisense orientations. Based on results of repeat-masking, as much as 10% of BAC end sequences from catfish, which is a random survey of the genome, contain some remnant of Tc1 elements, suggesting that these elements are present in the catfish genome as numerous, small remnants of the transposons. Phylogenetic analysis allowed comparison of catfish Tc1 transposase types with those found in other vertebrate and invertebrate species. In spite of the existence of many types of Tc1-like sequences that are not yet able to be placed in clades with strong statistical support, it is clear that multiple families of Tc1-like elements exist in channel catfish.