Estrodiol-17β and aromatase inhibitor treatment induced alternations of genome-wide DNA methylation pattern in Takifugu rubripes gonads.

Estradiol-17β (E2) and aromatase inhibitor (AI) exposure can change the phenotypic sex of fish gonads. To investigated whether alterations in DNA methylation is involved in this process, the level of genome-wide DNA methylation in Takifugu rubripes gonads was quantitatively analyzed during the E2-induced feminization and AI-induced masculinization processes in this study. The methylation levels of the total cytosine (C) in control-XX(C-XX), control-XY (C-XY), E2-treated-XY (E-XY) and AI-treated-XX (AI-XX) were 9.11%, 9.19%, 8.63% and 9.23%, respectively. In the C-XX vs C-XY comparison, 4,196 differentially methylated regions (DMRs) overlapped with the gene body of 2,497 genes and 608 DMRs overlapped with the promoter of 575 genes. In the E-XY vs C-XY comparison, 6,539 DMRs overlapped with the gene body of 3,416 genes and 856 DMRs overlapped with the promoter of 776 genes. In the AI-XX vs C-XX comparison, 2,843 DMRs overlapped with the gene body of 1,831 genes and 461 DMRs overlapped with the promoter of 421 genes. Gonadal genomic methylation mainly occurred at CG sites and the genes that overlapped with DMRs on CG context were most enriched in the signaling pathways related to gonad differentiation, such as the Wnt, TGF-β, MAPK, CAM and GnRH pathways. The DNA methylation levels of steroid synthesis genes and estrogen receptor genes promoter or gene body were negative correlated with their expression. After bisulfite sequencing verification, the DNA methylation level of the amhr2 promoter in XY was increased after E2 treatment, which consistent with the data from the genome-wide DNA methylation sequencing. In C-XY group, the expression of amhr2 was significantly higher than that in E-XY (p < 0.05). Additionally, dnmt1, which is responsible for methylation maintenance, expressed at similar level in four groups (p > 0.05). dnmt3, tet2, and setd1b, which were responsible for methylation modification, expressed at significantly higher levels in E-XY compared to the C-XY (p < 0.05). Dnmt3 and tet2 were expressed at significantly higher levels in AI-XX than that in C-XX (p < 0.05). These results indicated that E2 and AI treatment lead to the aberrant genome-wide DNA methylation level and expression level of dnmt3, tet2, and setd1b in T. rubripes gonad.