DNA Methylation Profiling in Genetically Selected Clarias magur (Hamilton, 1822) Provides Insights into the Epigenetic Regulation of Growth and Development.

DNA methylation is an epigenetic alteration that impacts gene expression without changing the DNA sequence affecting an organism's phenotype. This study utilized a reduced representation bisulfite sequencing (RRBS) approach to investigate the patterns of DNA methylation in genetically selected Clarias magur stocks. RRBS generated 249.22 million reads, with an average of 490,120 methylation sites detected in various parts of genes, including exons, introns, and intergenic regions. A total of 896 differentially methylated regions (DMRs) were identified; 356 and 540 were detected as hyper-methylated and hypo-methylated regions, respectively. The DMRs and their association with overlapping genes were explored using whole genome data of magur, which revealed 205 genes in exonic, 210 in intronic, and 480 in intergenic regions. The analysis identified the maximum number of genes enriched in biological processes such as RNA biosynthetic process, response to growth factors, nervous system development, neurogenesis, and anatomical structure morphogenesis. Differentially methylated genes (DMGs) such as myrip, mylk3, mafb, egr3, ndnf, meis2a, foxn3, bmp1a, plxna3, fgf6, sipa1l1, mcu, cnot8, trim55b, and myof were associated with growth and development. The selected DMGs were analyzed using real-time PCR, which showed altered mRNA expression levels. This work offers insights into the epigenetic mechanisms governing growth performance regulation in magur stocks. This work provides a valuable resource of epigenetic data that could be integrated into breeding programs to select high-performing individuals.

Muscle Transcriptome Sequencing Revealed Thermal Stress-Responsive Regulatory Genes in Farmed Rohu, Labeo rohita (Hamilton, 1822).

Rohu, Labeo rohita, is one of the most important aquaculture species in the Indian subcontinent. Understanding the molecular-level physiological responses to thermal stress or climate change is essential. In the present work, transcriptome sequencing was carried out in the muscle tissue of the rohu in response to heat stress (35 °C) in comparison with the control (28 °C). A total of 125 Gb of sequence data was generated, and the raw-reads were filtered and trimmed, which resulted in 484 million quality reads. Reference-based assembly of reads was performed using L. rohita genome, and a total of 90.17% of reads were successfully mapped. A total of 37,462 contigs were assembled with an N50 value of 1854. The differential expression analysis revealed a total of 107 differentially expressed genes (DEGs) (15 up-, 37 down-, and 55 neutrally regulated) as compared to the control group (Log2FC > 2, P < 0.05). Gene enrichment analysis of DEGs indicates that transcripts were associated with molecular, biological, and cellular activities. The randomly selected differentially expressed transcripts were validated by RT-qPCR and found consistent expression patterns in line with the RNA-seq data. Several transcripts such as SERPINE1(HSP47), HSP70, HSP90alpha, Rano class II histocompatibility A beta, PGC-1 and ERR-induced regulator, proto-oncogene c-Fos, myozenin2, alpha-crystallin B chain-like protein, angiopoietin-like protein 8, and acetyl-CoA carboxylases have been identified in muscle tissue of rohu that are associated with stress/immunity. This study identified the key biomarker SERPINE1 (HSP47), which showed significant upregulation (~ 2- to threefold) in muscle tissue of rohu exposed to high temperature. This study can pave a path for the identification of stress-responsive biomarkers linked with thermal adaptations in the farmed carps.

Characterization, Docking and Molecular Dynamics Simulation of Gonadotropin-Inhibitory Hormone Receptor (GnIHR2) in Labeo Catla.

BACKGROUND/AIMS: GnIH receptors (GnIHRs) belong to the family of G-protein coupled receptors (GPCRs) and play a key role in the regulation of reproduction from fishes to mammals, either by inhibiting or stimulating the expression of gonadotropins. The aim of this study was to characterize GnIH receptor (GnIHR2) from Indian Major Carp, Labeo catla and its docking and simulation with GnIH antagonist RF313. METHODS: The full length sequence of GnIHR2 was obtained with RACE PCR. The docking analysis of RF313 with GnIHR2 receptor was performed with AutoDock v. 4.2.6 and molecular dynamics (MD) simulation with GROMACS 5.0. RESULTS: In the present study, we cloned full-length cDNA (1733 bp) of GnIHR2 from the brain of L. catla. The phylogenetic analysis showed clustering of catla GnIHR2 with goldfish and zebrafish in the GPR147 group. L. catla GnIHR2 receptor comprised seven transmembrane domains and the 3D-structure was predicted by I-TASSER tool. The docking analysis revealed high binding affinity (-11.6 kcal/mol) of GnIH antagonist, RF313 towards GnIHR2 receptor. The primary bonds involved were alkyl and hydrogen bonds while the amino acids participated were proline 43, 210, 339, cysteine 214, leucine 211, serine 213 and phenylalanine 338. The MD simulation analysis of docked complex for 100 nano-seconds (ns) in the lipid membrane environment showed the stability of the complex with time. CONCLUSION: Our study showed that GnIH antagonist, RF313 interact tightly with the GnIH receptor, GnIHR2 of L. catla. To the best of our knowledge, this is the first report on computational modelling and MD simulation of GnIH receptor in fishes. This will help in functional characterization studies of GnIH/GnIHR system in vertebrates.

Molecular characterization of gonadotropin-inhibitory hormone (GnIH) gene and effect of intramuscular injection of GnIH peptide on the reproductive axis in Catla catla.

Gonadotropin-inhibitory hormone (GnIH) plays an important role in reproduction by inhibiting the expression of gonadotropins in birds and mammals, but in fishes, it is ambiguous. In this study, we cloned 606 bp long cDNA of GnIH from Catla catla brain (cGnIH). The encoded preproGnIH peptide generated three putative peptides (cGnIH-I, -II, -III) of different size. Phylogenetic analysis of GnIH showed clustering of different peptide sequence with its orthologs in separate clades. The real-time PCR analysis showed the expression of cGnIH in brain, gonads, intestine, stomach, heart, gill and liver with the highest expression in the brain and gonads of both sexes. The basal GnIH mRNA expression was higher in spawning and spent phase of the male brain and spawning phase of the female brain. In testis, the expression was highest in spent phase, while in ovary the expression did not change significantly during reproductive phases. The in vivo experiment of cGnIH-III peptide exhibited the higher expression of HPG axis genes, lhb, fshb, cgnrh, kiss2 and kiss1r and serum hormone level of LH and FSH as soon as 3 h after the intramuscular delivery. These results suggest that the GnIH is positively involved in regulation of reproduction in HPG axis of C. catla.

DNA barcoding of marine ornamental fishes from India.

India has rich marine ornamental fish diversity with 400 fish species distributed in Gulf of Munnar/Palk Bay, Gulf of Kutch, and in reefs around Andaman & Nicobar and Lakshadweep Islands. Marine ornamental fish identification at the field level is very difficult because of their high diversity and profound changes in appearance during their developmental stages and camouflage. To facilitate ornamental fish trading with ease and in compliance with the biodiversity act, DNA barcoding technique could be used to accurately identify species. In this study, DNA barcodes were generated for 31 species of commercially important marine ornamental fishes from India. The average genetic distance (K2P model) within species, genus, and family was 0.446, 13.08, and 20.09%, respectively. Intraspecific variation has increased several folds (15-20 times) after including conspecific sequences from different geographical locations. The presence of allopatric lineages/cryptic species was observed in the Indo-pacific region. The NJ tree constructed based on K2P values showed distinct clusters shared by congeneric species specific to populations.