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

A high-quality chromosome-level genome assembly of rohu carp, Labeo rohita, and its utilization in SNP-based exploration of gene flow and sex determination

Labeo rohita (rohu) is a carp important to aquaculture in South Asia, with a production volume close to Atlantic salmon. While genetic improvements to rohu are ongoing, the genomic methods commonly used in other aquaculture improvement programs have historically been precluded in rohu, partially due to the lack of a high-quality reference genome. Here we present a high-quality de novo genome produced using a combination of next-generation sequencing technologies, resulting in a 946 Mb genome consisting of 25 chromosomes and 2,844 unplaced scaffolds. Notably, while approximately half the size of the existing genome sequence, our genome represents 97.9% of the genome size newly estimated here using flow cytometry. Sequencing from 120 individuals was used in conjunction with this genome to predict the population structure, diversity, and divergence in three major rivers (Jamuna, Padma, and Halda), in addition to infer a likely sex determination mechism in rohu. These results demonstrate the utility of the new rohu genome in modernizing some aspects of rohu genetic improvement programs.

De novo assembly, transcriptome characterization and marker discovery in Indian major carp, Labeo rohita through pyrosequencing

Labeo rohita, one of the Indian major carps, is the most popular culture species in Indian subcontinent due to its consumer preference and delicacy. A selective breeding program for harvest body weight has resulted in an average genetic gain of 17% per generation. Transcriptome resource for this species is scanty. Here, we have characterized the liver and muscle transcriptomes of rohu using Roche 454 GS-FLX next generation sequencing platform. In total, 1.2 million reads were generated, de novo assembly and clustering resulted in 4171 transcripts. Out of these, 4171 had significant blast hit against NCBI nr database, and 2130 transcripts were successfully annotated. In total, 289 SSRs were identified with an identification rate of 5.8%, and dinucleotide repeat motifs were observed to be the most abundant SSRs. Further, 2231 putative SNPs were identified with high confidence. Validation of eight putative SNPs using Sanger sequencing resulted in 100% true SNPs. Significant allelic imbalance of M1, M4 and M5 loci between growth selected and control individual were observed. Furthermore, 13 transcription factors were identified in the present study belonging to six different transcription factor families. The present study demonstrated the utility of RNAseq to develop genomics resources in non-model fish species, and the marker resources developed would support the genetic improvement program of this species.

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.

Dynamic Gene Expression and Alternative Splicing Events Demonstrate Co-Regulation of Testicular Differentiation and Maturation by the Brain and Gonad in Common Carp

The common carp (Cyprinus carpio) accounts for approximately 10% of the annual freshwater aquaculture production and is an ideal model to study cyprinidae reproduction. Female common carp grow faster than the males; therefore, related research presents an opportunity with high application value. Although we have a detailed understanding of common carp's early gonadal differentiation process, information about genome-wide gene expression, regulation, and underlying molecular mechanisms during this process remain limited. Here, time-course data comprising six key stages during testicular differentiation and maturation were investigated to further understand the molecular mechanisms underlying the testicular development in cyprinid species. After integrating these time-series data sets, common carp genome, including 98,345 novel transcripts and 3,071 novel genes were re-annotated and precisely updated. Gene co-expression network analysis revealed that the ubiquitin-mediated proteolysis pathway was essential for metabolism during testicular differentiation in the endocrine system of C. carpio. Functional enrichment analyses indicated that genes mainly related to amino acid metabolism and steroid hormone synthesis were relatively highly expressed at the testicular undifferentiation stages, whereas genes associated with cell cycle and meiosis were expressed from the beginning of testicular differentiation until maturation. The dynamics of alternative splicing events demonstrated that exon skipping accounted for majority of the alternative splicing events in the testis and the brain during gonad development. Notably, several potential male-specific genes (fanci and sox30) and brain-specific genes (oxt, gad2, and tac1, etc.) were identified. Importantly, we traversed beyond the level of transcription to test for stage- and gonad-specific alternative splicing patterns between the brain and testis. This study is the first to describe a comprehensive landscape of alternative splicing events and gene expression patterns during gonadogenesis in common carp. This work is extremely valuable to elucidate the mechanisms underlying gonadal differentiation in Cyprinidae as well as other fish species.

Molecular characterization, computational analysis and expression profiling of Dmrt1 gene in Indian major carp, Labeo rohita (Hamilton 1822)

Sexual dimorphism of fish morphology, physiology and behavior is diverse and complex in nature. Doublesex and mab-3 related transcription factor (Dmrt) is a large protein family whose function is sexual development and differentiation in vertebrates. Here, we report a full-length cDNA sequence of Labeo rohita (rohu) Dmrt1 of 907 bp length having 798 bp of open reading frame encoding 265 amino acids. The molecular weight of rohu DMRT1 protein was found to be 28.74 KDa and isoelectric point was 7.53. DMRT1 protein contains 23 positively and 24 negatively charged amino acids with a GRAVY score of -0.618. A characteristic DM domain was found in DMRT1 protein, which is a novel DNA-binding domain. Phylogenetic analysis showed maximum similarity with Cyprinus carpio when compared with DMRT1 of other vertebrates. Molecular docking study identified active sites to be targeted for drug designing. Rohu DMRT1 was observed to interact with other proteins such as FOXL2, CYP19a1a, AMH and SOX9a. Differential expression study revealed higher expression in testis tissue implying its role in male sex differentiation and testicular development. The information generated in the present work could facilitate further research to resolve the issues related to gonadal maturation and reproduction of commercially important aquaculture species.

De novo sequencing and comparative analysis of testicular transcriptome from different reproductive phases in freshwater spotted snakehead Channa punctatus

The spotted snakehead Channa punctatus is a seasonally breeding teleost widely distributed in the Indian subcontinent and economically important due to high nutritional value. The declining population of C. punctatus prompted us to focus on genetic regulation of its reproduction. The present study carried out de novo testicular transcriptome sequencing during the four reproductive phases and correlated differential expression of transcripts with various testicular events in C. punctatus. The Illumina paired-end sequencing of testicular transcriptome from resting, preparatory, spawning and postspawning phases generated 41.94, 47.51, 61.81 and 44.45 million reads, and 105526, 105169, 122964 and 106544 transcripts, respectively. Transcripts annotated using Rattus norvegicus reference protein sequences and classified under various subcategories of biological process, molecular function and cellular component showed that the majority of the subcategories had highest number of transcripts during spawning phase. In addition, analysis of transcripts exhibiting differential expression during the four phases revealed an appreciable increase in upregulated transcripts of biological processes such as cell proliferation and differentiation, cytoskeleton organization, response to vitamin A, transcription and translation, regulation of angiogenesis and response to hypoxia during spermatogenically active phases. The study also identified significant differential expression of transcripts relevant to spermatogenesis (mgat3, nqo1, hes2, rgs4, cxcl2, alcam, agmat), steroidogenesis (star, tkt, gipc3), cell proliferation (eef1a2, btg3, pif1, myo16, grik3, trim39, plbd1), cytoskeletal organization (espn, wipf3, cd276), sperm development (klhl10, mast1, hspa1a, slc6a1, ros1, foxj1, hipk1), and sperm transport and motility (hint1, muc13). Analysis of functional annotation and differential expression of testicular transcripts depending on reproductive phases of C. punctatus helped in developing a comprehensive understanding on genetic regulation of spermatogenic and steroidogenic events in seasonally breeding teleosts. Our findings provide the basis for future investigation on the precise role of testicular genes in regulation of seasonal reproduction in male teleosts.

India's Computational Biology Growth and Challenges

India's computational science is growing swiftly due to the outburst of internet and information technology services. The bioinformatics sector of India has been transforming rapidly by creating a competitive position in global bioinformatics market. Bioinformatics is widely used across India to address a wide range of biological issues. Recently, computational researchers and biologists are collaborating in projects such as database development, sequence analysis, genomic prospects and algorithm generations. In this paper, we have presented the Indian computational biology scenario highlighting bioinformatics-related educational activities, manpower development, internet boom, service industry, research activities, conferences and trainings undertaken by the corporate and government sectors. Nonetheless, this new field of science faces lots of challenges.

Preliminary genetic linkage map of Indian major carp, Labeo rohita (Hamilton 1822) based on microsatellite markers

Linkage map with wide marker coverage is an essential resource for genetic improvement study for any species. Sex-averaged genetic linkage map of Labeo rohita, popularly known as 'rohu', widely cultured in the Indian subcontinent, was developed by placing 68 microsatellite markers generated by a simplified method. The parents and their F1 progeny (92 individuals) were used as segregating populations. The genetic linkage map spans a sex-averaged total length of 1462.2 cM, in 25 linkage groups. The genome length of rohu was estimated to be 3087.9 cM. This genetic linkage map may facilitate systematic searches of the genome to identify genes associated with commercially important characters and marker-assisted selection programmes of this species.

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.