In some fish species, high or low temperature can switch the sex determination mechanisms and induce fish sex reversal when the gonads are undifferentiated. During this high or low temperature-induced sex reversal, the expressions of many genes are altered. However, genome-wide DNA methylation changes in fish gonads after high or low temperature treatment are unclear. Herein, we compared the global DNA methylation changes in the gonads from control females (CF), control males (CM), high temperature-treated females (TF), and high temperature-induced males (IM) from the F8 family of Nile tilapia (Oreochromis niloticus) using methylated DNA immunoprecipitation sequencing. The DNA methylation level in CF was higher than that in CM for various chromosomes. Both females and males showed an increase in methylation levels on various chromosomes after high-temperature induction. We identified 64,438 (CF/CM), 63,437 (TF/IM), 98,675 (TF/CF), 235,270 (IM/CM) and 119,958 (IM/CF) differentially methylated regions (DMRs) in Nile tilapia gonads, representing approximately 0.70% (CF/CM), 0.69% (TF/IM), 1.07% (TF/CF), 2.56% (IM/CM), and 1.30% (IM/CF)of the length of the genome. A total of 89 and 65 genes that exhibited DMRs in their gene bodies and promoters were mapped to the Nile tilapia genome. Furthermore, more than half of the genes with DMRs in the gene body in CF/CM were also included in the IM/CM, TF/CF, TF/IM, and IM/CF groups. Additionally, many important pathways, including neuroactive ligand-receptor interaction, extracellular matrix-receptor interaction, and biosynthesis of unsaturated fatty acids were identified. This study provided an important foundation to investigate the molecular mechanism of high temperature-induced sex reversal in fish species.

Background: Nowadays, the molecular mechanisms governing TSD (temperature-dependent sex determination) or GSD + TE (genotypic sex determination + temperature effects) remain a mystery in fish. Methods: We developed three all-female families of Nile tilapia (Oreochromis niloticus), and the family with the highest male ratio after high-temperature treatment was used for transcriptome analysis. Results: First, gonadal histology analysis indicated that the histological morphology of control females (CF) was not significantly different from that of high-temperature-treated females (TF) at various development stages. However, the high-temperature treatment caused a lag of spermatogenesis in high-temperature-induced neomales (IM). Next, we sequenced the transcriptome of CF, TF, and IM Nile tilapia. 79, 11,117, and 11,000 differentially expressed genes (DEGs) were detected in the CF–TF, CF–IM, and TF–IM comparisons, respectively, and 44 DEGs showed identical expression changes in the CF–TF and CF–IM comparisons. Principal component analysis (PCA) indicated that three individuals in CF and three individuals in TF formed a cluster, and three individuals in IM formed a distinct cluster, which confirmed that the gonad transcriptome profile of TF was similar to that of CF and different from that of IM. Finally, six sex-related genes were validated by qRT-PCR. Conclusions: This study identifies a number of genes that may be involved in GSD + TE, which will be useful for investigating the molecular mechanisms of TSD or GSD + TE in fish.

Northern snakehead (Channa argus) is an indigenous fish species and is one of popularly cultured snakeheads in China and other Asian countries. Unfortunately, Nocardia seriolae infections have caused considerable losses in the snakehead aquaculture industry. However, the infectivity and the immune response induced by N. seriolae in snakehead are unclear. In order to better understand the immune response of Northern snakehead in a series of time points after N. seriolae challenge, we conducted the transcriptomic comparison in snakehead spleen at 48, 96, and 144 h after the challenge of N. seriola against their control counterparts. Gene annotation and pathway analysis of differentially expressed genes (DEGs) were carried out to understand the functions of the DEGs. Additionally, protein-protein interaction networks were conducted to obtain the interaction relationships of immune-related DEGs. These results revealed the expression changes of multiple DEGs and signaling pathways involved in immunity during N. seriolae infection, which will facilitate our comprehensive understanding of the mechanisms involved in the immune response to bacterial infection in the northern snakehead.