Ebinur Lake virus (EBIV) is a recently identified orthobunyavirus with broad host range and zoonotic potential, posing a public health risk. However, the mechanisms underlying EBIV pathogenesis and host innate immune responses remain unclear. Here, we investigated the pattern recognition receptors (PRRs) responsible for sensing EBIV infection and subsequent pathogenesis. EBIV infects diverse cell types and exhibits broad tissue tropism in vivo. In vitro, RIG-I was essential for type I interferon (IFN-I) and inflammatory responses in HEK293 and A549 cells. In contrast, both RIG-I and MDA5 contributed to IFN-I induction in Huh-7 and HCT116 cells, correlating with the specific accumulation of viral dsRNA intermediates in these cell types. Both RIG-I and MDA5 preferentially recognize RNA derived from the viral S segment; however, they have different abilities in sensing incoming viral genomic RNA bearing a 5'-phosphate motif and the replication intermediates. In vivo, RIG-I deficiency severely impairs host defense, while MDA5 deficiency has a more restricted effect in the spleen and liver. In addition to RIG-I and MDA5, TLR7, which is predominantly expressed in dendritic cells, also plays a crucial role for host defense by mediating systemic inflammatory cytokine production without significantly impacting IFN-I response. Our findings suggest that multiple innate sensing receptors, including RIG-I, MDA5, and TLR7, are differentially involved in host defense against EBIV by mediating IFN-I and inflammatory responses, respectively, in a cell-specific manner.IMPORTANCEThis study elucidates the complex mechanisms by which host RIG-I, MDA5, and TLR7 sense the emerging EBIV and trigger cell-specific immune responses. These findings not only clarify crucial aspects of EBIV-host interactions, particularly the differential sensing of viral RNA by distinct PRRs, but also underscore how this differential sensing dictates cell-specific innate immune activation (IFN-I vs. inflammatory responses) and viral pathogenesis, providing critical insights for understanding and combating EBIV and related emerging bunyaviruses.

Crimean-Congo hemorrhagic fever virus (CCHFV) is the most widespread tick-born zoonotic bunyavirus that causes severe hemorrhagic fever and death in humans. CCHFV enters the cell via clathrin-mediated endocytosis which is dependent on its surface glycoproteins. However, the cellular receptors that are required for CCHFV entry are unknown. Here we show that the low density lipoprotein receptor (LDLR) is an entry receptor for CCHFV. Genetic knockout of LDLR impairs viral infection in various CCHFV-susceptible human, monkey and mouse cells, which is restored upon reconstitution with ectopically-expressed LDLR. Mutagenesis studies indicate that the ligand binding domain (LBD) of LDLR is necessary for CCHFV infection. LDLR binds directly to CCHFV glycoprotein Gc with high affinity, which supports virus attachment and internalization into host cells. Consistently, a soluble sLDLR–Fc fusion protein or anti-LDLR blocking antibodies impair CCHFV infection into various susceptible cells. Furthermore, genetic knockout of LDLR or administration of an LDLR blocking antibody significantly reduces viral loads, pathological effects and death following CCHFV infection in mice. Our findings suggest that LDLR is an entry receptor for CCHFV and pharmacological targeting of LDLR may provide a strategy to prevent and treat Crimean-Congo hemorrhagic fever.

Interferons (IFNs) activate JAK-STAT pathways to induce downstream effector genes for host defense against invaded pathogens and tumors. Here both type I (β) and II (γ) IFNs are shown that can activate the transcription factor IRF3 in parallel with STAT1. IRF3-deficiency impairs transcription of a subset of downstream effector genes induced by IFN-β and IFN-γ. Mechanistically, IFN-induced activation of IRF3 is dependent on the cGAS-STING-TBK1 axis. Both IFN-β and IFN-γ cause mitochondrial DNA release into the cytosol. In addition, IFNs induce JAK1-mediated tyrosine phosphorylation of cGAS at Y214/Y215, which is essential for its DNA binding activity and signaling. Furthermore, deficiency of cGAS, STING, or IRF3 impairs IFN-β- or IFN-γ-mediated antiviral and antitumor activities. The findings reveal a novel IRF3 activation pathway parallel with the canonical STAT1/2 activation pathways triggered by IFNs and provide an explanation for the pleiotropic roles of the cGAS-STING-IRF3 axis in host defense.

Growth is one of the most important traits of commercial fishes. Bighead carp are intensively cultured in China and some Asian countries; however, the genetic basis of growth is less understood. In this study, we calculated the heritability of growth-related traits in a mixed population of 840 bighead carp and estimated growth-related QTL in a family of bighead carp using bulked segregation analysis of resequencing data. High phenotypic (0.70-0.95) and genetic (0.77-0.97) correlations were found between all growth-related traits. Estimated heri-tability of 0.20 +/- 0.11, 0.28 +/- 0.14, 0.32 +/- 0.16, 0.14 +/- 0.11 and 0.12 +/- 0.09 were obtained for body weight, body height, head length, total length, and body length in bighead carp with seven month old. In the BSA analysis for progeny with extreme segregated body weight of a single family, three regions located in LG01, LG11 and LG22, respectively with a total length of 17.84 Mb which were obtained according to correlation thresholds. Annotation of genes within this region yielded 120 and 61 nonsynonymous and code-shifting mutations, respectively, by comparing SNPs and InDel between the two growth difference groups. Finally, 14 candidate genes, including T-box transcription factor 3, transcription factor E2F2, paired box gene 9, and fibroblast growth factor 6b were identified. Association analysis and mRNA levels revealed that paired box gene 9 and fibroblast growth factor 6b were significantly associated with growth-related traits in extreme growth samples. These results provide insights into the genetic mechanisms underlying the regulation of growth and development in the bighead carp.

Growth is an important economically trait for aquatic animals. The popularity of farmed channel catfish (Ictalurus punctatus) in China has recently surged, prompting a need for research into the genetic mechanisms that drive growth and development to expedite the selection of fast-growing variants. In this study, the brain, liver and muscle transcriptomes of channel catfish between fast-growing and slow-growing groups were analyzed using RNA-Seq. Totally, 63, 110 and 86 differentially expressed genes (DEGs) were from brain, liver and muscle tissues. DEGs are primarily involved in growth, development, metabolism and immunity, which are related to the growth regulation of channel catfish, such as growth hormone receptor b (ghrb), fibroblast growth factor receptor 4 (fgfr4), bone morphogenetic protein 1a (bmp1a), insulin-like growth factor 2a (igf2a), collagen, type I, alpha 1a (col1a1a), acyl-CoA synthetase long chain family member 2 (acsl2) and caveolin 1 (cav1). This study advances our knowledge of the genetic mechanisms accounting for differences in growth rate and offers crucial gene resources for future growth-related molecular breeding programs in channel catfish.

Growth is an economically important trait in bighead carp and other aquaculture species that affects production efficiency. Interestingly, the head of the bighead carp has a high market value in China; therefore, it is important to study the genetic bases of both growth and body shape traits. A genome-wide association study was performed based on 2b-RAD sequencing of 776 individuals to identify SNPs associated with growth and body shape traits, including body weight, body length, body height, and deheaded body length. In total, 26 significant and 19 suggestive SNPs were identified, and more than half of these significant SNPs were clustered in LG16. Two LGs (LG16 and LG21) contained QTLs associated with body weight. Fourteen SNPs of LG16 and two LG21 SNPs were found to be associated with body length. For body height, 12 significantly associated SNPs were identified in LG16. Additionally, 12 SNPs of LG16 and 3 SNPs of LG21 were found to be associated with deheaded body length. Forty-three genes were significantly or suggestively associated with body shape/growth traits based on GWAS results, 18 of which were candidate genes for all BW, BL, BH, and DBL traits. One of these genes, fndc5b , was selected for further analyses. Association analysis revealed that one SNP (g.245 C > T) in the introns of fndc5b was significantly associated with growth-related traits in growth-extreme samples. The mRNA levels of fndc5b in the brains of the lightweight group were significantly higher than those of the heavy-weight group. This study helps to reveal the genetic structure of growth and body development in fish and provides candidate genes for future molecular marker-assisted selection for fast growth and better body conformation in bighead carp.

Sensing of cytosolic DNA of microbial or cellular/mitochondrial origin by cGAS initiates innate immune responses via the adaptor protein STING. It remains unresolved how the activity of STING is balanced between a productive innate immune response and induction of autoimmunity. Here we show that interferon regulatory factor 8 (IRF8) is essential for efficient activation of STING-mediated innate immune responses in monocytes. This function of IRF8 is independent of its transcriptional role in monocyte differentiation. In uninfected cells, IRF8 remains inactive via sequestration of its IRF-associated domain by its N- and C-terminal tails, which reduces its association with STING. Upon triggering the DNA sensing pathway, IRF8 is phosphorylated at Serine 151 to allow its association with STING via the IRF-associated domain. This is essential for STING polymerization and TBK1-mediated STING and IRF3 phosphorylation. Consistently, IRF8-deficiency impairs host defense against the DNA virus HSV-1, and blocks DNA damage-induced cellular senescence. Bone marrow-derived mononuclear cells which have an autoimmune phenotype due to deficiency of Trex1, respond to IRF-8 deletion with reduced pro-inflammatory cytokine production. Peripheral blood mononuclear cells from systemic lupus erythematosus patients are characterized by elevated phosphorylation of IRF8 at the same Serine residue we find to be important in STING activation, and in these cells STING is hyper-active. Taken together, the transcription-independent function of IRF8 we describe here appears to mediate STING activation and represents an important regulatory step in the cGAS/STING innate immune pathway in monocytes.

Gynogenesis is an effective technique for establishing pure line and validating possible chromosomal mode of sex determination in fish. As one of the major aquaculture fish historically cultured in China, bighead carp (Hypophthalmichthys nobilis) is famous for its market value of the head, but slow in its genetic improvement due to a long period of sexual maturation. In this study, females of bighead carp (F0 dams) with fast growth and bigger head were chose to induce first and second generation of meiotic gynogenesis (GF1 and GF2 for short) by using UV-irradiated heterologous sperm of common carp (Cyprinus carpio). On average, diploid gynogens were obtained by cold shock with a survival rate of 7.36% for GF1 and 6.24% for GF2, respectively. Results of morphology comparison and parentage assignment showed that all survival gynogens were genetically derived from the dams. Compared with their F0 parents, heterozygosities of GF1 and GF2 gynogens decreased by 42% and 60% when analyzed using 50 SSRs, and by 35% and 53% when using 1505 SNP markers, respectively. Sexing of the GF1 and GF2 gynogens showed all-female genotypes and phenotypes when examined by male-specific DNA markers, histology section and morphology of gonads. Overall, large number of meiotic gynogens have been successfully induced in two consecutive generations, providing valuable resources for fast domestication and genetic improvement in bighead carp. Our evidence of sex ratios in both GF1 and GF2 gynogens also strongly supports a proposed chromosomal mechanism of sex determination (XX/XY) for this large-size cyprinid fish.

MITA (also known as STING) is an ER-located adaptor protein, which mediates DNA-triggered innate immune response and is critically involved in autoimmune diseases and tumorigenesis. MITA is regulated by post-translational modifications, but how post-transcriptional mechanisms are involved in the regulation of MITA is still largely unknown. Here, we identified the RNA-binding protein LUC7L2 as a negative regulator of DNA virus-triggered innate immune response. LUC7L2-deficient mice exhibited resistance to lethal herpes simplex virus 1 (HSV-1) infection and reduced HSV-1 loads in the brain. Mechanistically, LUC7L2 directly bound to intron 3 of MITA precursor messenger RNA, inhibited its splicing and promoted its nonsense-mediated decay, leading to its downregulation at protein level. LUC7L2-deficient cells had markedly increased MITA level, leading to heightened innate antiviral response. Finally, LUC7L2 was induced following HSV-1 infection. Our findings reveal a feedback negative post-transcriptional regulatory mechanism for regulation of MITA-mediated innate immune response to viral and aberrant cellular DNA.

Head-related traits are economically important in fish because they will essentially affect fillet yields. Bighead carp (Hypophthalmichthys nobilis) is an important aquaculture fish, and interestingly its head has higher culinary and market values than its fillet in China. In this study, using 776 bighead carp individuals and high-throughput SNP genotyping technology, we performed genome-wide association study (GWAS) to identify genomic regions and candidate genes potentially associated with head-size and head-shape traits such as head length (HL), head width (HW), head height (HH), head length/head height ratio (HL/HH), head length/head width ratio (HL/HW) and head length/body length ratio (HL/BL). The morphometrics of most of these traits were significantly correlated (P < 0.01). GWAS detected eleven significant SNPs and twelve suggestive SNPs associated with head size of bighead carp, which were mainly located on linkage group (LG) 16 and four of these SNPs were commonly shared among three head-size traits. For HL, four significant and three suggestive SNPs were on LG16. For HW, five significant and four suggestive SNPs were identified on LG3, LG16 and LG21, respectively. For HH, two significant and five suggestive SNPs were detected on LG11 and LG16. Genes involved in cell proliferation and bone development, such as cc2b, prdm1b, fndc5b, lbr, bmp8a and ATP6v1cb, were found surrounding SNPs in the QTL regions. For head shape, seventeen and six SNPs were significantly related to HL/HH and HL/HW, respectively, and underlying candidate genes were identified from genomic regions of bighead carp surrounding these SNPs, including ptch1, col9a1a, tgfbr2, hecw2a, zbtb42 and sema7a. This study shed lights on understanding genetic architecture of head size/shape traits in fish and facilitated to identify candidate genes for markerassisted selection towards breeding bigger headed carp in future.

Growth, one of the most important traits monitored in domestic animals, is essentially associated with bone development. To date, no large-scale transcriptome studies investigating bone development in bighead carp have been reported. In this study, we applied Isoform-sequencing technology to uncover the entire transcriptomic landscape of the bighead carp (Hypophthalmichthys nobilis) in early growth stage, and obtained 63,873 non-redundant transcripts, 20,907 long non-coding RNAs, and 1,579 transcription factors. A total of 381 alternative splicing events were seen in the frontal and parietal bones with another 784 events simultaneously observed in the vertebral bones. Coupling this to RNA sequencing (RNA-seq) data, we identified 27 differentially expressed unigenes (DEGs) in the frontal and parietal bones and 45 DEGs in the vertebral bones in the fast-growing group of fish, when compared to the slow-growing group of fish. Finally, 15 key pathways and 20 key DEGs were identified and found to be involved in regulation of early growth such as energy metabolism, immune function, and cytoskeleton function and important cellular pathways such as the arginine and proline metabolic pathway (p4ha1), FoxO signaling pathway (sgk1), cell adhesion molecules (b2m, ptprc, and mhcII), and peroxisome proliferator-activated receptor signaling pathway (scd). We established a novel full-length transcriptome resource and combined it with RNA-seq to elucidate the mechanism of genetic regulation of differential growth in bighead carp. The key DEGs identified in this study could fuel further studies investigating associations between growth and bone development and serve as a source of potential candidate genes for marker-assisted breeding programs.

Cyclic GMP-AMP synthase (cGAS) senses double-stranded DNA and synthesizes the second messenger cyclic GMP-AMP (cGAMP), which binds to mediator of IRF3 activation (MITA) and initiates MITA-mediated signaling, leading to induction of type I interferons (IFNs) and other antiviral effectors. Human cytomegalovirus (HCMV), a widespread and opportunistic pathogen, antagonizes the host antiviral immune response to establish latent infection. Here, we identified HCMV tegument protein UL94 as an inhibitor of the cGAS-MITA-mediated antiviral response. Ectopic expression of UL94 impaired cytosolic double-stranded DNA (dsDNA)- and DNA virustriggered induction of type I IFNs and enhanced viral replication. Conversely, UL94 deficiency potentiated HCMV-induced transcription of type I IFNs and downstream antiviral effectors and impaired viral replication. UL94 interacted with MITA, disrupted the dimerization and translocation of MITA, and impaired the recruitment of TBK1 to the MITA signalsome. These results suggest that UL94 plays an important role in the immune evasion of HCMV. IMPORTANCE Human cytomegalovirus (HCMV), a large double-stranded DNA (dsDNA) virus, encodes more than 200 viral proteins. HCMV infection causes irreversible abnormalities of the central nervous system in newborns and severe syndromes in organ transplantation patients or AIDS patients. It has been demonstrated that HCMV has evolved multiple immune evasion strategies to establish latent infection. Previous studies pay more attention to the mechanism by which HCMV evades immune response in the early phase of infection. In this study, we identified UL94 as a negative regulator of the innate immune response, which functions in the late phase of HCMV infection.

As aquatic animals, fishes often encounter various situations of low oxygen, and they have evolved the ability to respond to hypoxia stress. Studies of physiological and molecular responses to hypoxia stress are essential to clarify genetic mechanisms underlying hypoxia tolerance in fish. In this study, we performed acute hypoxia treatment in juvenile bighead carp (Hypophthalmicthys nobilis) by decreasing water O2 from 6.5 mg/L to 0.5 mg/L in three hours. This hypoxia stress resulted in a significant increase in blood lactate and serum glucose. Comparisons of heart transcriptome among hypoxia tolerant (HT), hypoxia sensitive (HS), and normoxia control (NC) groups showed that 820, 273, and 301 differentially expressed genes (DEGs) were identified in HS vs. HT, NC vs. HS, and NC vs. HT (false discovery rate (FDR) < 0.01, Fold Change> 2), respectively. KEGG pathway enrichment showed that DEGs between HS and HT groups were mainly involved in mitogen-activated protein kinase (MAPK) signaling, insulin signaling, apoptosis, tight junction and adrenergic signaling in cardiomyocytes pathways, and DEGs in MAPK signaling pathway played a key role in cardiac tolerance to hypoxia. Combined with the results of our previous cDNA-amplified fragment length polymorphism (cDNA-AFLP) analysis of hypoxia stress in this species, such genes as stbp2, ttn, mapk, kcnh, and tnfrsf were identified in both studies, representing the significance of these DEGs in hypoxia tolerance in bighead carp. These results provide insights into the understanding of genetic modulations for fish heart coping with hypoxia stress and generate basic resources for future breeding studies of hypoxia resistance in bighead carp.

Mediator of IRF3 activation (MITA, also known as stimulator of interferon genes, STING) senses the second messenger cyclic GMP-AMP (cGAMP) which is synthesized upon DNA virus infection and activates innate antiviral immune response. It has been demonstrated that the activity of MITA is delicately regulated by various post-translational modifications including polyubiquitination. In this study, we identified the deubiquitinating enzyme USP44 as a positive regulator of MITA. USP44 is recruited to MITA following DNA virus infection and removes K48-linked polyubiquitin moieties from MITA at K236, therefore prevents MITA from proteasome mediated degradation. USP44-deficiency results in acceleration of HSV-1-induced degradation of MITA and reduced induction of type I interferons (IFNs) and proinflammatory cytokines. Consistently, Usp44-/- mice are more susceptible to HSV-1 infection as indicated by higher tissue viral titers, greater tissue damage and lower survival rate. These findings suggest that USP44 plays a specific and critical role in the regulation of innate immune response against DNA viruses.

Body deformity occurs both in wild and farmed fishes, which is one of the most challenging problems for aquaculture industry. In most cases, such body deformities are linked to skeletal deformities. Currently, very limited information is available on skeletal deformities of farmed fish species which may be caused by genetic factor. In this study, we performed muscle and vertebra transcriptome analyses in body deformity and normality of bighead carp Hypophthalmichthys nobilis (from one meiotic gynogenesis family) using RNA-Seq. A total of 43,923 and 44,416 unigenes were predicted in muscles and vertebrae, respectively. Based on these data, we further explored the gene expression profiles in gynogenetic normal and abnormal bighead carp. No differentially expressed gene (DEG) was found in transcriptome data of muscles. Totally, 20 key DEGs were identified in transcriptome data of vertebrae, such as low density lipoprotein-related protein 2 (lrp2), bone morphogenetic protein 2B (bmp2b) and collagen alpha-1(IV) (col4a1). 12 potential pathways were also identified in vertebra transcriptome data, which were mainly involved in development, growth, cytoskeleton and energy metabolism, such as MAPK signaling pathway, regulation of actin cytoskeleton and TGF-beta signaling pathway. Results of this study will be informative for the understanding of genetic mechanisms for body shape formation and also provide potential candidate genes for selection program involved in body shape and skeletal development in H. nobilis.

Trafficking of toll-like receptor 3 (TLR3) from the endoplasmic reticulum (ER) to endolysosomes and its subsequent proteolytic cleavage are required for it to sense viral double-stranded RNA (dsRNA) and trigger antiviral response, yet the underlying mechanisms remain enigmatic. We show that the E3 ubiquitin ligase TRIM3 is mainly located in the Golgi apparatus and transported to the early endosomes upon stimulation with the dsRNA analog poly(I:C). TRIM3 mediates K63-linked polyubiquitination of TLR3 at K831, which is enhanced following poly(I:C) stimulation. The polyubiquitinated TLR3 is recognized and sorted by the ESCRT (endosomal sorting complex required for transport) complexes to endolysosomes. Deficiency of TRIM3 impairs TLR3 trafficking from the Golgi apparatus to endosomes and its subsequent activation. Trim3(-/-) cells and mice express lower levels of antiviral genes and show lower levels of inflammatory response following poly(I:C) but not lipopolysaccharide (LPS) stimulation. These findings suggest that TRIM3-mediated polyubiquitination of TLR3 represents a feedback-positive regulatory mechanism for TLR3-mediated innate immune and inflammatory responses.

Recognition of viral RNA by the retinoic acid-inducible gene-I (RIG-I)-like receptors (RLRs), including RIG-I and MDA5, initiates innate antiviral responses. Although regulation of RLR-mediated signal transduction has been extensively investigated, how the recognition of viral RNA by RLRs is regulated remains enigmatic. In this study, we identified heterogeneous nuclear ribonucleoprotein M (hnRNPM) as a negative regulator of RLR-mediated signaling. Overexpression of hnRNPM markedly inhibited RNA virus-triggered innate immune responses. Conversely, hnRNPM-deficiency increased viral RNA-triggered innate immune responses and inhibited replication of RNA viruses. Viral infection caused translocation of hnRNPM from the nucleus to the cytoplasm. hnRNPM interacted with RIG-I and MDA5, and impaired the binding of the RLRs to viral RNA, leading to inhibition of innate antiviral response. Our findings suggest that hnRNPM acts as an important decoy for excessive innate antiviral immune response.

STING plays central roles in the innate immune response to pathogens that contain DNA. Sensing cytoplasmic DNA by cyclic GMP-AMP synthase produces cyclic GMP-AMP, which binds to and activates STING and induces STING translocation from the endoplasmic reticulum to the perinuclear microsome. However, this trafficking process has not been fully elucidated yet. In this study, we identified YIPF5 as a positive regulator of STING trafficking. YIPF5 is essential for DNA virus- or intracellular DNA-triggered production of type I IFNs. Consistently, knockdown of YIPF5 impairs cellular antiviral responses to DNA virus. Mechanistically, YIPF5 interacts with both STING and components of COPII, facilitating STING recruitment to COPII in the presence of cytoplasmic dsDNA. Furthermore, knockdown of components of COPII inhibits DNA virus-triggered production of type I IFNs, suggesting that COPII is involved in innate immune responses to DNA viruses. Collectively, our findings demonstrate that YIPF5 positively regulates STING-mediated innate immune responses by recruiting STING to COPII-coated vesicles and facilitating STING trafficking from the endoplasmic reticulum to Golgi, providing important insights into the molecular mechanisms of intracellular DNA-stimulated STING trafficking and activation.

MITA (also called STING) is a central adaptor protein in innate immune response to cytosolic DNA. Cellular trafficking of MITA from the ER to perinuclear microsomes after DNA virus infection is critical for MITA activation and onset of innate antiviral response. Here we found that SNX8 is a component of DNA-triggered induction of downstream effector genes and innate immune response. Snx8(-/-) mice infected with the DNA virus HSV-1 exhibited lower serum cytokine levels and higher viral titers in the brains, resulting in higher lethality. Mechanistically, SNX8 recruited the class III phosphatylinositol 3-kinase VPS34 to MITA, which is required for trafficking of MITA from the ER to perinuclear microsomes. Our findings suggest that SNX8 is a critical component in innate immune response to cytosolic DNA and DNA virus.

Feed efficiency is an economically crucial trait for cultured animals, however, progress has been scarcely made in the genetic analyses of feed conversion efficiency (FCE) in fish because of the difficulties in measurement of trait phenotypes. In the present investigation, we present the first application of RNA sequencing (RNA-Seq) combined with differentially expressed genes (DEGs) analysis for identification of functional determinants related to FCE at the gene level in an aquaculture fish, crucian carp (Carassius auratus). Brain tissues of six crucian carp with extreme FCE performances were subjected to transcriptome analysis. A total of 544,612 unigenes with a mean size of 644.38 bp were obtained from Low- and High-FCE groups, and 246 DEGs that may be involved in FCE traits were identified in these two groups. qPCR confirmed that genes previously identified as up- or down-regulated by RNA-Seq were effectively up- or down-regulated under the studied conditions. Thirteen key genes, whose functions are associated with metabolism (Dgkk, Mgst3 and Guk1b), signal transduction (Vdnccsa1b, Tgfα, Nr4a1 and Tacr2) and growth (Endog, Crebrtc2, Myh7, Myh1, Myh14 and Igfbp7) were identified according to GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) annotations. Our novel findings provide useful pathway information and candidate genes for future studies of genetic mechanisms underlying FCE in crucian carp.