Samiran Nandi

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.

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

BackgroundThrough the Community Sequencing Program, a catfish EST sequencing project was carried out through a collaboration between the catfish research community and the Department of Energys Joint Genome Institute. Prior to this project, only a limited EST resource from catfish was available for the purpose of SNP identification.ResultsA 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.ConclusionsThis 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.

Towards the Ictalurid Catfish Transcriptome: Generation and Analysis of 31,215 Catfish ESTs.

BackgroundEST 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.ResultsA 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.ConclusionThe 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.

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.

Simple sequence repeats (SSRs) markers in fish genomic research and their acceleration via next-generation sequencing and computational approaches

AbstractSimple sequence repeats (SSRs) are becoming a choice of markers in fish genetic research due to their abundance in the genome, co-dominant nature, high polymorphism and ability to reproduce. Thus, in this review, we have discussed regarding SSRs markers developed in fishes using different techniques. These markers have been used for revealing genetic variability, strain and species identification, genetic linkage map construction and parentage assignment in fish genetic research. Recently, high-throughput sequencing platform has been widely used in non-model fishes for genome/transcriptome sequencing to understand genomic information. The rapid progress in fish genomic research has been made due to sequencing platform along with their low cost for sequencing and use of the advanced computational tools for generated data analysis. We have shown that different next-generation sequencing platforms have been applied in the genomic studies for SSRs markers identification in fishes with evidence. We have depicted the use of various computational tools/algorithms for the SSRs identification from genome/transcriptome data. However, we also highlighted existing challenges in high-throughput sequencing data analysis as well as the current need of computationally deep analysis software/tools/expertise. The purpose of this review is to get envisage on the various possibilities, which can be harnessed via these new technologies and advanced computational tools for SSRs marker development via genome/transcriptome sequencing of aquaculture species.

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.

Survey of the transcriptome of Brevibacillus borstelensis exposed to low temperature shock

Molecular mechanisms underlying the ability of Brevibacillus borstelensis to survive and adapt to various environmentally relevant stresses are poorly understood. To define organisms molecular response to low temperature, gene expression profile of B. borstelensis at 20 °C was carried out by high-throughput sequencing technology. A total of 4579 transcripts with a maximum transcript length of 9919 bp were annotated. Gene expression profiling identified 712 genes that were significantly up- or down-regulated during cold shock. Functional categorization of the differentially expressed genes revealed that response to stress, regulation of transcription, transport, signal transduction and cytoplasm were the differentially regulated processes. The microbial stress responsive genes (hsp90, hslU, grpE, dnaK, dnaJ, hslV) and genes under regulatory adaptive responses (rpoN) were identified. The gene encoding cold shock protein purine nucleoside phosphorylase was found to be remarkably up-regulated. RT-PCR experiments carried out on genes expressed under cold shock independently verified the transcriptome data results. In addition, a large number of genes encoding hypothetical protein were identified. The brief survey of the transcripts obtained in response to cold shock underlines the survival strategy of thermophilic bacteria exposed to low temperature environment, which is further helpful in generating genetic information associated with this bacteria.

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. 16GB, 72bp 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.