Abstract

Takifugu rubripes (T. rubripes) is a valuable commercial fish, and Cryptocaryon irritans (C. irritans) has a significant impact on its aquaculture productivity. DNA methylation is one of the earliest discovered ways of gene epigenetic modification and also an important form of modification, as well as an essential type of alteration that regulates gene expression, including immune response. To further explore the anti-infection mechanism of T. rubripes in inhibiting this disease, we determined genome-wide DNA methylation profiles in the gill of T. rubripes using whole-genome bisulfite sequencing (WGBS) and combined with RNA sequence (RNA-seq). A total of 4659 differentially methylated genes (DMGs) in the gene body and 1546 DMGs in the promoter between the infection and control group were identified. And we identified 2501 differentially expressed genes (DEGs), including 1100 upregulated and 1401 downregulated genes. After enrichment analysis, we identified DMGs and DEGs of immune-related pathways including MAPK, Wnt, ErbB, and VEGF signaling pathways, as well as node genes prkcb, myca, tp53, and map2k2a. Based on the RNA-Seq results, we plotted a network graph to demonstrate the relationship between immune pathways and functional related genes, in addition to gene methylation and expression levels. At the same time, we predicted the CpG island and transcription factor of four immune-related key genes prkcb and mapped the gene structure. These unique discoveries could be helpful in the understanding of C. irritans pathogenesis, and the candidate genes screened may serve as optimum methylation-based biomarkers that can be utilized for the correct diagnosis and therapy T. rubripes in the development of the ability to resist C. irritans infection.