Transcriptome assembly and functional annotation
In the present study, kidney transcriptome generated from two different temperature conditions represented a complete complement of expressed genes, including the major immune genes as well as the important pathways. A total of 84.47 Gb reads were generated and after trimming and removal of contaminants, high-quality paired-end reads and clustering resulted in 1.275 million transcripts, which confirmed 98.60% completeness of the transcriptome.
DEGs and analysis for functional enrichment
Differential gene expression analysis of common carp kidney tissues from 20oC experimental temperature versus 30oC water temperature (control) conditions revealed 588 DEGs and 253 and 335 genes were found upregulated and downregulated, respectively. Under biological processes, most significant GO terms (FDR<0.05) were DNA integration (GO:0015074), followed by transmembrane transport (GO:0055085) and intracellular signal transduction (GO:0035556), while under molecular functions, GO:0003676-nucleic acid binding followed by GO:0005509-calcium ion binding and GO:0005085-guanyl-nucleotide exchange factor activity and under Cellular component, GO:0016020-membranes and GO:0005886 plasma membrane.
Molecular pathway analysis of these DEGs identified 118 genes in Organismal Systems, with a maximum number of 41 genes under immune system; 76 in environmental information processing with 57 genes under signal transduction, 29 in genetic information processing with 9 genes under folding, sorting and degradation and 59 genes under metabolism with 16 genes under amino acid metabolism. A total of 56 genes were classified under cellular processes, with 26 genes under transport and catabolism.
Under pathway analysis, maximum DEGs were observed in signal transduction (57) and signalling molecules and interaction (28) under environmental information processing, followed by endocrine system (47) and immune system (41) under organismal systems. Under signal transduction and signaling molecules and interaction, maximum number of DEGs was observed in Calcium signaling pathway [PATH: ko04020], MAPK signaling pathway [PATH: ko04010], Apelin signaling pathway [PATH: ko04371], cAMP signaling pathway [PATH: ko04024], PI3K-Akt signaling pathway [PATH:ko04151] with DEGs 13, 13, 8, 8 and 8, respectively. Under metabolism, pathways of Arginine and proline metabolism [PATH: ko00330] and Purine metabolism [PATH: ko00230] contained 5 and 4 DEGs, respectively; and under cell growth and death category, apoptosis [path: ko04210] 6 DEGs. Thus, the present study revealed an increase in purine, arginine and proline metabolisms, which may play a vital role in countering cold stress. These pathways were reported to be involved in energy production and mitigating the effects of cold stress
(Melis et al., 2017, Xu et al., 2018) and the results are in line with the earlier studies that fish tolerate cold temperature stress by regulating its metabolism in response to increased energy demand
(Schleger et al., 2021). Since the current study found that the cold temperature had no effect on the glycolysis or glycogenesis pathways, the genes involved in the purine, arginine and proline metabolism pathways appear to be important for energy production and to combat cold stress in common carp.
Molecular pathways identified under immune system category of organismal systems identified 40 DEGs, out of which 7 genes were under complement and coagulation cascades pathway, Chemokine signalling pathway and Platelet activation followed by 5 genes each under B cell receptor signaling, C-type lectin receptor signaling and T cell receptor signaling pathways and 5 genes under IL-17 signaling pathway (Table 2).
DEGs classification under abiotic stress
With reference to
Danio rerio database, 5 GO terms with 22 DEGs were identified, with 10 DEGs under response to oxidative stress (GO:0001666, GO:0055093, GO:0071456), followed by 4 under response to heat (GO:0009408, GO:0034605), 5 under response to mechanical stimulus (GO:0009612, GO:0071260) and 2 under response to cold (GO:0009409, GO:0070417) (Fig 1). The most significant pathways identified were MAPK signalling, Focal adhesion and necroptosis for these genes. The s
tress-regulated MAPK signaling pathway mainly involves kinases, which are reported to be mainly the mitogen-activated protein kinase (MAPK) family members
(Cross et al., 2000), that is involved in numerous significant cellular functions, including proliferation, differentiation and migration of cells (https://www.kegg.jp/pathway/hsa04010).. Interesting points of concern are significant upregulation of genes: Mitogen-activated protein kinase 4 (MAP4K4), cytosolic phospholipase A2 (cPLA2a) and Thrombospondin-1 (tsp-1). Up-regulation of MAP4K4 gene is indicative of positive stress response
(Singh et al., 2023) andcPLA2a activation in renal epithelium has been suggested to promote renal repair by initiating partial de-differentiation, proliferation and migration
(Montford et al., 2016). and Third up-regulated gene under abiotic stress under the present study, that is
tsp-1 under colder conditions, which is a component of the PI3K-Akt signaling pathway, has been reported with a role in the immune system regulation
(Bao et al., 2022) and in adaptation to varying temperatures
(Zhang et al., 2021).
Differentially expressed immune and related genes
Molecular pathways identified under immune system category of organismal systems identified 40 DEGs, out of which 7 genes were under complement and coagulation cascades pathway, Chemokine signaling pathway and Platelet activation followed by 5 genes each under B cell receptor signaling, C-type lectin receptor signaling and T cell receptor signaling pathways and 5 genes under IL-17 signaling pathway (Table 2). Out of 40 DEGs, 15 were immune genes, majority (12) under innate immune categorywhile only 3 adaptive immune DEG were identified (Table 3). These genes fall under 9 GO terms, of which under biological process, immune response (GO:0006955); under cellular component under membrane (GO:0016020); under molecular function signaling receptor activity (GO:0038023), carbohydrate binding (GO:0030246) and GTP binding (GO:0005525) (Fig 2). Three DEGs were found under under cytokine-cytokine receptor interaction (CCRI) pathway and 2 each under Chemokine signaling pathway and Viral-protein interaction with CCR (Fig 3). In the current study, these results involving DE immune genes point to highly significant groups, with the largest number of genes are under Response to stimulus and CCRI pathways (Fig 4), which are interlinked and share most of the associated DEGs. CCRI, one of the vital immunity pathways, has a crucial role in the inflammatory and defense process of the host
(Cheng et al., 2017). In an earlier study in transgenic zebrafish under cold stress
(Wang et al., 2014), where CCRI genes were mostly down-regulated. However, in the current study, up-regulated CCRI genes (
ccr1, 2) at 20oC might be providing the innate immune response in common carp.
Verma et al., (2021) also reported the chemokine signaling pathway to be a critical pathway, with the innate immune gene
ccr2 up-regulation in the early and middle stage infection of oomycete,
A. invadans in common carp. It is fascinating to find in the present study, a very high expression of another innate immune gene, complement component 7b, identified under C-C chemokine binding, which performs function in host defence against pathogens and promotion of inflammation through complement and coagulation pathway
(Bossi et al., 2009). Thus, the relative higher expression of the immune genes indicates enhanced pathogen defense
(Philominal et al., 2025).
Validation of expression pattern of immune genes
Out of the 11 genes
, rps20, hprt, eif5a, rlp7, b-actin, 18s-1, myb-2, anxa1a, errif1a, kat2b-1 and
ep300a-1 tested for their stability in expression,
eif5a was found to show most stable expression, so it was utilized as the house keeping gene, for normalization of Ct values of target genes. The target genes,
ccr2, si:dkey-241l7.3, cc7 and
ifi44/loc795887 were tested and all these genes showed similar trend in the expression patterns for transcriptome and qRT-PCR experiments, except
ccr7 at 20
oC.
Thus, the present transcriptomic study points out that cold resistant common carp can modulate the cold stress positively and also its immune responses, which may facilitate the resistant responses to the disease, to which cold susceptible species succumbed to. However, deep studies are required for the combinations of these factors and processes, to reveal overall modulatory effects.