17α-Methyltestosterone (MT) is a widely used androgen for all-male fish production in aquaculture. However, the molecular mechanism underlying MT-induced masculinization remains unclear. In this study, we aim to identify the key gene responsible for MT-induced masculinization using the Nile tilapia (Oreochromis niloticus) amhy, dmrt1, and gsdf mutants, which exhibit male-to-female sex reversal. Nile tilapia fry from these three mutant lines were treated with 50 μg/g MT from 5 to 30 days after hatching (dah). The results showed that amhy and gsdf mutants, but not dmrt1 mutants, were masculinized by the MT treatment. Gonadal transcriptome analysis revealed that genes involved in steroidogenesis and germ cell development in MT-treated dmrt1 mutants exhibited a similar expression pattern to that of the wild type (WT) XX. In addition, the dmrt1 mutants cannot be masculinized by co-treatment with MT and the aromatase inhibitor fadrozole. The MT treatment completely blocked early steroidogenic enzyme (Star2, Cyp17a2, and Cyp19a1a) expression independent of amhy, gsdf, and dmrt1. A luciferase analysis showed that MT directly suppressed basal and Sf-1-activated cyp19a1a promoter activity through ara and arb in cultured HEK293 cells. Furthermore, MT treatment inhibited germ cell proliferation in amhy and gsdf mutants but not in dmrt1 mutants. Consistently, dmrt1 expression was induced in MT-treated WT XX, -amhy, and -gsdf mutants. Taken together, these results suggest that dmrt1 is indispensable for MT-induced masculinization in Nile tilapia and that MT functions by inhibiting early steroid synthesis and activating dmrt1 to promote testis development.

Neural crest cells give rise to the craniofacial skeleton and multiple pigment cell lineages, yet how duplicated developmental regulators partition their ancestral functions after teleost-specific whole-genome duplication remains unclear. Here, we employed CRISPR/Cas9 to generate alx4a and alx4b single and double mutants in Nile tilapia (Oreochromis niloticus). By integrating phenotype, skeleton, transcriptome, quantitative PCR, and AlphaFold-based structural modeling analyses, we revealed their functional divergence. Loss of alx4a caused a regionally restricted reduction in iridophore-derived reflectance and abnormal cranial morphology, whereas alx4b single mutants showed no obvious phenotype under the conditions examined. By contrast, double mutants exhibited an almost complete loss of iridophore-derived structural coloration and substantially more severe cranial defects, accompanied by reduced calcein labeling in the opercular region, consistent with altered cranial mineralization. Skin transcriptomic and quantitative PCR analyses revealed marked downregulation of pnp4a and tfec, which are associated with iridophore differentiation and coloration, whereas no significant expression differences were detected for the iridophore survival-related genes ltk and mpv17. AlphaFold2-assisted HDOCK protein–DNA modeling yielded more favorable docking metrics for Alx4a than for Alx4b with the pnp4a promoter, supporting a potential Alx4a–pnp4a promoter interaction that requires experimental validation. In contrast, no significant genotype-dependent differences were detected in the measured abundance of melanophores, xanthophores, or erythrophores, and no obvious difference in gross dorsal-fin spine formation was observed under the conditions examined. Together, these findings reveal unequal functional contributions of alx4a and alx4b, with alx4a acting as the dominant paralog in iridophore-associated structural coloration and both paralogs contributing unequally to cranial development, and support pnp4a as a candidate downstream gene associated with Alx4a activity.