A detailed transcriptome study uncovers the epigenetic characteristics associated with Aromatase inhibitor-induced masculinization in Takifugu rubripes larvae gonads.

BACKGROUND: Takifugu rubripes is an economically valuable fish species in Asia. The implementation of all-male culture for T. rubripes is highly anticipated in aquaculture. Aromatase inhibitor (AI, letrozole) treatment was found to be an efficient method to induced masculinization in T. rubripes, as reported in our previous study. Here, to further explore the underlying regulation mechanism of AI-induced masculinization, a whole-transcriptome analysis comparing was conducted between AI-induced masculinized XX (AI-XX) gonads and control (Con) gonads in T. rubripes. RESULTS: In Con-XX/Con-XY comparison, 1,172 differential expression (DE) mRNAs, 129 DEmiRNAs, 210 DElncRNAs, and 4 DEcircRNAs were identified. In the Con-XX/AI-XX comparison, 1,329 DEmRNAs, 174 DEmiRNAs, 6 DEcircRNAs and 280 DElncRNAs were found. Con-XX/Con-XY and Con-XX/AI-XX comparisons shared 690 DEmRNAs, 50 DEmiRNAs, 3 DEcircRNAs, and 105 DElncRNAs. The analyses of protein-protein interaction (PPI) and competitive endogenous RNA (ceRNA) network identified interactions among these shared DERNAs. Kcnh2b, trim27, cnnm2b, reln, cckbra, pkd1l2, steap4, gsg1l, hamp, and foxg1c were predicted as the top ten of hub genes. miRNAs included miRNA-27 family and miRNA-489 family showed targeting relationship with hub genes. GO and KEGG functional enrichment analysis showed that the targeted genes were mainly enriched in GO:0065008 regulation of biological quality and TGF-beta signaling pathway. qPCR validation confirmed the differential expression of selected mRNAs, and ncRNAs. CONCLUSIONS: This research comprehensively reveals the potential regulatory effects of ncRNAs on cellular motility, fate regulation, and hormonal regulation during gonadal masculinization in T. rubripes. It may provide significant insights into the regulation mechanisms underlying sex reversal in fish.