Genome-Wide DNA Methylation and RNA Analysis Reveal Potential Mechanism of Resistance to Streptococcus agalactiae in GIFT Strain of Nile Tilapia (Oreochromis niloticus )

DNA methylation regulates growth traits by influencing metabolic pathways in Pacific white shrimp (Litopenaeus vannamei)

BACKGROUND

DNA methylation is a critical epigenetic modification that dynamically regulates gene expression associated with economic traits. Pacific white shrimp (Litopenaeus vannamei) is one of the most important aquatic species for culturing, and growth trait is one of the most important economic traits for its production. However, research on DNA methylation regulation of growth traits is still at an early stage. This study explored DNA methylome dynamics and their associations with the regulatory mechanism behind growth traits using full-subfamily individuals with discrepant growth performance.

RESULTS

The DNA methylation-related genes in L. vannamei were identified, and the expression of DNA methylation genes showed significantly higher levels in the slow growth (SG) group compared to the fast-growing (FG) individuals. The Whole Genome Bisulfite Sequencing (WGBS) analysis revealed that the methylation levels in the muscles of shrimp were notably decreased in SG individuals compared to FG individuals. A total of 532 differentially methylated promoters and 2,067 differentially methylated regions were identified. Through integrative analysis of DNA methylation and transcriptomic data from SG and FG group shrimp, a total of 47 genes were screened out with differential methylation levels (DMGs) and expression levels (DEGs). Functional enrichment analysis revealed that the overlapping DEGs/DMGs were enriched mainly in metabolic pathways, starch and sucrose metabolism, linoleic acid metabolism, ascorbate and aldarate metabolism, pentose and glucuronate interconversions.

CONCLUSIONS

DNA methylation plays a role in the regulation of growth traits in L. vannamei. The level of DNA methylation was found to be negatively correlated with growth traits. Through comprehensive analysis, it was discovered that DNA methylation predominantly affects growth performance by up-regulating the expression of genes involved in metabolic pathways, such as glucose metabolism and amino acid metabolism in L. vannamei. This suggests a higher metabolism activity in SG individuals derived DNA methylation to cope with some unknown internal stress or environmental stress rather than being allocated for growth.

Insights into the DNA methylation of Portunus trituberculatus in response to Vibrio parahaemolyticus infection

Vibrio parahaemolyticus is the main pathogen causing acute hepatopancreatic necrotic disease in crustaceans. To elucidate the epigenetic regulatory mechanism of crustacean resistance to V. parahaemolyticus infection, we conducted artificial infection studies on Portunus trituberculatus. The results showed that the mortality rate reached the highest at 12 h of artificial infection, which was 23.69 %. At 72 h after V parahaemolyticus infection, the expression level of DNA demethylase (ten-eleven-translocation protein) Tet was significantly decreased, the expression of DNA methyltransferase Dnmt3B fluctuated significantly. Based on the differential expression levels of Tet and Dnmt3B. We depict for DNA methylation profiles of the whole genome of P. trituberculatus at single-base resolution by using whole-genome bisulfite sequencing (WGBS) on hemolymph tissues. The overall DNA methylation level was low at 2.16 % in P. trituberculatus hemolymph. A total of 2590 differentially methylated regions (DMRs) were identified, of which 1329 were hypermethylated and 1261 were hypomethylated, and 1389 genes were annotated in these DMRs. Differently methylated genes (DMGs) were significantly enriched in ribosomes (KO03010), protein kinases (KO01001), cell cycle (HSA04110), endocrine resistance (HSA01522) and FoxO signaling pathway (KO04068). Finally, we selected six differentially methylated genes for quantitative analysis. The results showed that DNA methylation not only has a negative regulatory effect on gene expression, but also has a positive regulatory effect. These results indicated that DNA methylation in the regulation of genes involved in immune responses contributes to the resistance of P. trituberculatus to V. parahaemolyticus, which is valuable for understanding how crustaceans regulate the innate immune system to defend against bacterial infections.

Comparative analysis of T-cell immunity between Streptococcus agalactiae susceptible and resistant tilapia (Oreochromis niloticus)

Nile tilapia (Oreochromis niloticus) is one of the important economic fish species cultured worldwide. However, Streptococcus agalactiae has emerged as a significant bacterial threat, severely impacting the economy of tilapia industry. The immune response underlying the resistance of tilapia to S. agalactiae are not well understood, hindering the reasonable evaluation of breeding and the formulation of effective strategies. In this study, we investigated the differences in T-cell immunity between S. agalactiae-resistant and -susceptible tilapia. Compared with susceptible tilapia, resistant tilapia exhibited a higher percentage of T cells and BrdU+ T cells during infection, indicating a superior proliferative capacity. Whether infected or not, T cells from resistant fish demonstrated a greater ability to resist apoptosis. Additionally, T cell effector genes, including interleukin (IL)-2, interferon (IFN)-γ, perforin A, and granzyme B were expressed at higher levels in resistant tilapia after infection. Along with these T-cell immune responses, resistant fish showed more effective clearance of infection. Our study elucidates the T-cell immune responses in resistant tilapia, which may contribute to the high resistance of tilapia to S. agalactiae, and provide valuable theoretical references for the selection and evaluation of disease-resistant fish strains in the future.

Whole-genome DNA methylation profiling revealed epigenetic regulation of NF-κB signaling pathway involved in response to Vibrio alginolyticus infection in the Pacific oyster, Crassostrea gigas

DNA methylation, an essential epigenetic alteration, is tightly linked to a variety of biological processes, such as immune response. To identify the epigenetic regulatory mechanism in Pacific oyster (Crassostrea gigas), whole-genome bisulfite sequencing (WGBS) was conducted on C. gigas at 0 h, 6 h, and 48 h after infection with Vibrio alginolyticus. At 6 h and 48 h, a total of 11,502 and 14,196 differentially methylated regions (DMRs) were identified (p<0.05, FDR<0.001) compared to 0 h, respectively. Gene ontology (GO) analysis showed that differentially methylated genes (DMGs) were significantly enriched in various biological pathways including immunity, cytoskeleton, epigenetic modification, and metabolic processes. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that transcription machinery (ko03021) is one of the most important pathways. Integrated transcriptome and methylome analyses allowed the identification of 167 and 379 DMG-related DEGs at 6 h and 48 h, respectively. These genes were significantly enriched in immune-related pathways, including nuclear factor kappa B (NF-κB) signaling pathway (ko04064) and tumor necrosis factor (TNF) signaling pathway (ko04668). Interestingly, it's observed that the NF-κB pathway could be activated jointly by TNF Receptor Associated Factor 2 (TRAF2) and Baculoviral IAP Repeat Containing 3 (BIRC3, the homolog of human BIRC2) which were regulated by DNA methylation in response to the challenge posed by V. alginolyticus infection. Through this study, we provided insightful information about the epigenetic regulation of immunity-related genes in the C. gigas, which will be valuable for the understanding of the innate immune system modulation and defense mechanism against bacterial infection in invertebrates.

Hypogonadotropic hypogonadism in male tilapia lacking a functional rln3b gene

Relaxin 3 is a neuropeptide that plays a crucial role in reproductive functions of mammals. Previous studies have confirmed that rln3a plays an important role in the male reproduction of tilapia. To further understand the significance of its paralogous gene rln3b in male fertility, we generated a homozygous mutant line of rln3b in Nile tilapia. Our findings indicated that rln3b mutation delayed spermatogenesis and led to abnormal testes structure. Knocking out rln3b gene resulted in a decrease in sperm count, sperm motility and male fish fertility. TUNEL detection revealed a small amount of apoptosis in the testes of rln3b-/- male fish at 390 days after hatching (dah). RT-qPCR analysis demonstrated that mutation of rln3b gene caused a significant downregulation of steroid synthesis-related genes such as cyp17a1, cyp11b2, germ cell marker gene, Vasa, and gonadal somatic cell marker genes of amh and amhr2. Furthermore, we found a significant down-regulation of hypothalamic-pituitary-gonadal (HPG) axis-related genes, while a significantly up-regulation of the dopamine synthetase gene in the rln3b-/- male fish. Taken together, our data strongly suggested that Rln3b played a crucial role in the fertility of XY tilapia by regulating HPG axis genes.

DNA Methylation Analysis Reveals Potential Mechanism in Takifugu rubripes Against Cryptocaryon irritans Infection

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.

DNA Methylation Profiling of Ovarian Tissue of Climbing Perch (Anabas testudienus) in Response to Monocrotophos Exposure

Epigenetic modifications like DNA methylation can alter an organism's phenotype without changing its DNA sequence. Exposure to environmental toxicants has the potential to change the resilience of aquatic species. However, little information is available on the dynamics of DNA methylation in fish gonadal tissues in response to organophosphates. In the present work, reduced-representation bisulfite sequencing was performed to identify DNA methylation patterns in the ovarian tissues of Anabas testudienus exposed to organophosphates, specifically monocrotophos (MCP). Through sequencing, an average of 41,087 methylated cytosine sites were identified and distributed in different parts of genes, i.e., in transcription start sites (TSS), promoters, exons, etc. A total of 1058 and 1329 differentially methylated regions (DMRs) were detected as hyper-methylated and hypo-methylated in ovarian tissues, respectively. Utilizing whole-genome data of the climbing perch, the DMRs, and their associated overlapping genes revealed a total of 22 genes within exons, 45 genes at transcription start sites (TSS), and 218 genes in intergenic regions. Through gene ontology analysis, a total of 16 GO terms particularly involved in ovarian follicular development, response to oxidative stress, oocyte maturation, and multicellular organismal response to stress associated with reproductive biology were identified. After functional enrichment analysis, relevant DMGs such as steroid hormone biosynthesis (Cyp19a, 11-beta-HSD, 17-beta-HSD), hormone receptors (ar, esrrga), steroid metabolism (StAR), progesterone-mediated oocyte maturation (igf1ar, pgr), associated with ovarian development in climbing perch showed significant differential methylation patterns. The differentially methylated genes (DMGs) were subjected to analysis using real-time PCR, which demonstrated altered gene expression levels. This study revealed a molecular-level alteration in genes associated with ovarian development in response to chemical exposure. This work provides evidence for understanding the relationship between DNA methylation and gene regulation in response to chemicals that affect the reproductive fitness of aquatic animals.

Genome-wide DNA methylation mediates the resistance to vibriosis in Cynoglossus semilaevis

Chinese tongue sole (Cynoglossus semilaevis) is an economically important marine fish in China. However, vibriosis has caused huge mortality and economic losses in its culturing industry. To reveal the effect of DNA methylation on the resistance to vibriosis in tongue sole, we conducted RNA sequencing and whole genome bisulfite sequencing (WGBS), and compared the gene expressions and DNA methylation patterns between the resistant and susceptible families. We identified a total of 741 significantly differentially expressed genes (DEGs) in kidney and 17460 differentially methylated genes (DMGs), which were both enriched in immune-related pathways, such as "cAMP signaling pathway" and "NOD-like receptor signaling pathway". Through the correlation analysis of DEGs and DMGs, we identified two important immune pathways, including "complement and coagulation cascades", and "cytokine-cytokine receptor interaction", which played important roles in regulating the inflammation level and immune homeostasis. For example, the expression of proinflammatory cytokine il17c was down-regulated under the regulation of DNA methylation; in addition, the expression of protease-activated receptor 3 (par3) was up-regulated, which could induce the up-expressionof il8. These results demonstrated that the regulation of DNA methylation on the genes involved in immune responses might contribute to the resistance to vibriosis in tongue sole, and provided a basis for the control of diseases in fish aquaculture.

Applying genetic technologies to combat infectious diseases in aquaculture

Disease and parasitism cause major welfare, environmental and economic concerns for global aquaculture. In this review, we examine the status and potential of technologies that exploit genetic variation in host resistance to tackle this problem. We argue that there is an urgent need to improve understanding of the genetic mechanisms involved, leading to the development of tools that can be applied to boost host resistance and reduce the disease burden. We draw on two pressing global disease problems as case studies-sea lice infestations in salmonids and white spot syndrome in shrimp. We review how the latest genetic technologies can be capitalised upon to determine the mechanisms underlying inter- and intra-species variation in pathogen/parasite resistance, and how the derived knowledge could be applied to boost disease resistance using selective breeding, gene editing and/or with targeted feed treatments and vaccines. Gene editing brings novel opportunities, but also implementation and dissemination challenges, and necessitates new protocols to integrate the technology into aquaculture breeding programmes. There is also an ongoing need to minimise risks of disease agents evolving to overcome genetic improvements to host resistance, and insights from epidemiological and evolutionary models of pathogen infestation in wild and cultured host populations are explored. Ethical issues around the different approaches for achieving genetic resistance are discussed. Application of genetic technologies and approaches has potential to improve fundamental knowledge of mechanisms affecting genetic resistance and provide effective pathways for implementation that could lead to more resistant aquaculture stocks, transforming global aquaculture.

Immunity of turbot Induced by inactivated vaccine of Aeromonas salmonicida from the perspective of DNA methylation

Introduction

DNA methylation was one of the most important modification in epigenetics and played an important role in immune response. Since the introduction of Scophthalmus maximus, the scale of breeding has continued to expand, during which diseases caused by various bacteria, viruses and parasites have become increasingly serious. Therefore, the inactivated vaccines have been widely researched and used in the field of aquatic products with its unique advantages. However, the immune mechanism that occurred in turbot after immunization with inactivated vaccine of Aeromonas salmonicida was not clear.

Methods

In this study, differentially methylated regions (DMRs) were screened by Whole Genome Bisulfite Sequencing (WGBS) and significantly differentially expressed genes (DEGs) were screened by Transcriptome sequencing. Double luciferase report assay and DNA pull-down assay were further verified the DNA methylation state of the gene promoter region affected genes transcriptional activity after immunization with inactivated vaccine of Aeromonas salmonicida.

Results

A total of 8149 differentially methylated regions (DMRs) were screened, in which there were many immune-related genes with altered DNA methylation status. Meanwhile, 386 significantly differentially expressed genes (DEGs) were identified, many of which were significantly enriched in Toll-like receptor signaling pathway, NOD-like receptor signaling pathway and C-type lectin receptor signaling pathway. Combined analysis of WGBS results and RNA-seq results, a total of 9 DMRs of negatively regulated genes are located in the promoter region, including 2 hypermethylated genes with lower expression and 7 hypomethylated genes with higher expression. Then, two immune-related genes C5a anaphylatoxin chemotactic receptor 1-like (C5ar1-Like) and Eosinophil peroxidase-like (EPX-Like), were screened to explore the regulation mechanism of DNA methylation modification on their expression level. Moreover, the DNA methylation state of the gene promoter region affected genes transcriptional activity by inhibiting the binding of transcription factors, which lead to changes in the expression level of the gene.

Discussion

We jointly analyzed WGBS and RNA-seq results and revealed the immune mechanism that occurred in turbot after immunized with inactivated vaccine of A. salmonicida from the perspective of DNA methylation.

Integrated chromatin accessibility and DNA methylation analysis to reveal the critical epigenetic modification and regulatory mechanism in gonadal differentiation of the sequentially hermaphroditic fish, Monopterus albus

BACKGROUND

Monopterus albus is a hermaphroditic and economically farmed fish that undergoes sex reversal from ovary to testis via ovotestis during gonadal development. The epigenetic changes that are associated with gonadal development in this species remain unclear.

METHODS

We produced DNA methylome, transcriptome, and chromatin accessibility maps of the key stages of gonad development: ovary, ovotestis, and testis. The expression of the key candidate genes was detected using qRT-PCR and in situ hybridization and the methylation levels were analysed using bisulphite sequencing PCR. Promoter activity and regulation were assessed using dual-luciferase reporter assays.

RESULTS

Gonadal development exhibits highly dynamic transcriptomic, DNA methylation, and chromatin accessibility changes. We found that DNA methylation status, especially of the transcription start site, was significantly negatively correlated with gene expression while chromatin accessibility exhibited no correlation with gene expression during gonadal development. The epigenetic signatures revealed many novel regulatory elements and genes involved in sex reversal, which were validated. DNA methylation detection and site mutation of plastin-2 promoter, as a candidate gene, revealed that DNA methylation could impact the binding of transcription factor dmrt1 and foxl2 through methylation and demethylation to regulate plastin-2 expression during gonadal development.

CONCLUSIONS

These data provide novel insights into epigenetic modification and help elucidate the potential molecular mechanism by which dynamic modification of DNA methylation plays a crucial role in gonadal development.

Genome-wide DNA methylation reveals potential epigenetic mechanism of age-dependent viral susceptibility in grass carp

Background

Grass carp are an important farmed fish in China that are infected by many pathogens, especially grass carp reovirus (GCRV). Notably, grass carp showed age-dependent susceptibility to GCRV; that is, grass carp not older than one year were sensitive to GCRV, while those over three years old were resistant to this virus. However, the underlying mechanism remains unclear. Herein, whole genome-wide DNA methylation and gene expression variations between susceptible five-month-old (FMO) and resistant three-year-old (TYO) grass carp were investigated aiming to uncover potential epigenetic mechanisms.

Results

Colorimetric quantification revealed that the global methylation level in TYO fish was higher than that in FMO fish. Whole-genome bisulfite sequencing (WGBS) of the two groups revealed 6214 differentially methylated regions (DMRs) and 4052 differentially methylated genes (DMGs), with most DMRs and DMGs showing hypermethylation patterns in TYO fish. Correlation analysis revealed that DNA hypomethylation in promoter regions and DNA hypermethylation in gene body regions were associated with gene expression. Enrichment analysis revealed that promoter hypo-DMGs in TYO fish were significantly enriched in typical immune response pathways, whereas gene body hyper-DMGs in TYO fish were significantly enriched in terms related to RNA transcription, biosynthesis, and energy production. RNA-seq analysis of the corresponding samples indicated that most of the genes in the above terms were upregulated in TYO fish. Moreover, gene function analysis revealed that the two genes involved in energy metabolism displayed antiviral effects.

Conclusions

Collectively, these results revealed genome-wide variations in DNA methylation between grass carp of different ages. DNA methylation and gene expression variations in genes involved in immune response, biosynthesis, and energy production may contribute to age-dependent susceptibility to GCRV in grass carp. Our results provide important information for disease-resistant breeding programs for grass carp and may also benefit research on age-dependent diseases in humans.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12979-022-00285-w.

Sequence, Expression, and Anti-GCRV Function of the Ferritin from the Grass Carp, Ctenopharyngodon idellus

Ferritin possesses an immune function to defend against pathogen infection. To elucidate the immunity-protecting roles of ferritin from Ctenopharyngodon idellus (Ciferritin) against virus infection, the cDNA and promoter sequences of Ciferritin were determined, and the correlations between Ciferrtin expressions and promoter methylation levels were analyzed. In addition, the functional role of Ciferrtin on GCRV (grass carp reovirus) infection was assessed. The full-length cDNA of Ciferritin is 1053 bp, consists of a 531 bp open-reading frame, and encodes 176 amino acids. Ciferritin showed the highest sequence identity with the ferritin middle subunit of Mylopharyngodon piceus (93.56%), followed by the subunits of Megalobrama amblycephala and Sinocyclocheilus rhinocerous. Ciferritin contains a conserved ferritin domain (interval: 10−94 aa), and the caspase recruitment domain (CARD) and Rubrerythrin domain were also predicted. In the spleen and kidney, significantly higher Ciferritin expressions were observed at 6, 12, 24, or 168 h post GCRV infection than those in the PBS injection group (p < 0.05). The Ciferrtin expression level in the progeny of maternal-immunized grass carp was significantly higher than that in the progeny of common grass carp (p < 0.05). Ciferritin promoter methylation level in the progeny from common grass carp was 1.27 ± 0.15, and in the progeny of the maternal-immunized group was 1.00 ± 0.14. In addition, methylation levels of “CpG9” and “CpG10” loci were significantly lower in the progeny of maternal-immunized fish than those in the common group. Except for the “CpG5”, methylation levels of all other detected “CpG” loci negatively correlated with Ciferritin expression levels. Furthermore, the total methylation level of “CpG1−10” negatively correlated with the Ciferritin expressions. The Ciferritin expression level was significantly up-regulated, and the VP7 protein levels were significantly reduced, at 24 h post GCRV infection in the Ciferritin over-expression cells (p < 0.05). The results from the present study provide sequence, epigenetic modification and expression, and anti-GCRV functional information of Ciferritin, which provide a basis for achieving resistance to GCRV in grass carp breeding.

Coiled-Coil Domain-Containing (CCDC) Proteins: Functional Roles in General and Male Reproductive Physiology

The coiled-coil domain-containing (CCDC) proteins have been implicated in a variety of physiological and pathological processes. Their functional roles vary from their interaction with molecular components of signaling pathways to determining the physiological functions at the cellular and organ level. Thus, they govern important functions like gametogenesis, embryonic development, hematopoiesis, angiogenesis, and ciliary development. Further, they are implicated in the pathogenesis of a large number of cancers. Polymorphisms in CCDC genes are associated with the risk of lifetime diseases. Because of their role in many biological processes, they have been extensively studied. This review concisely presents the functional role of CCDC proteins that have been studied in the last decade. Studies on CCDC proteins continue to be an active area of investigation because of their indispensable functions. However, there is ample opportunity to further understand the involvement of CCDC proteins in many more functions. It is anticipated that basing on the available literature, the functional role of CCDC proteins will be explored much further.

Multi-omics integration strategies for animal epigenetic studies - A review

Genome-wide studies provide considerable insights into the genetic background of animals; however, the inheritance of several heritable factors cannot be elucidated. Epigenetics explains these heritabilities, including those of genes influenced by environmental factors. Knowledge of the mechanisms underlying epigenetics enables understanding the processes of gene regulation through interactions with the environment. Recently developed next-generation sequencing (NGS) technologies help understand the interactional changes in epigenetic mechanisms. There are large sets of NGS data available; however, the integrative data analysis approaches still have limitations with regard to reliably interpreting the epigenetic changes. This review focuses on the epigenetic mechanisms and profiling methods and multi-omics integration methods that can provide comprehensive biological insights in animal genetic studies.