Cortisol induces masculinization of XX medaka through gonadal soma-derived growth factor (GSDF) and anti-Müllerian hormone receptor type 2 (AMHR2)

The medaka Oryzias latipes is a teleost fish with an XX/XY sex determination system similar to that of mammals. However, under high-temperature conditions, XX medaka are masculinized by an elevation of the key teleost glucocorticoid, cortisol. Cortisol inhibits female-type proliferation of germ cells and induces masculinization of XX medaka during gonadal sex differentiation. To identify masculinization mechanisms downstream of cortisol, we analysed the functions of gonadal soma-derived growth factor (gsdf) and anti-Müllerian hormone receptor type 2 (amhr2); these genes are known to play important roles in the inhibition of germ cell proliferation and male differentiation. We investigated the impact of gsdf and amhr2 on the proliferation of germ cells using gsdf knockout (KO) and amhr2 KO medaka. At hatching stage, loss of gsdf or amhr2 function recovered female-type proliferation in germ cells under cortisol treatment. Moreover, cortisol treatment of gsdf KO or amhr2 KO medaka did not induce masculinization of XX medaka. These results suggest that cortisol inhibits female-type proliferation of germ cells and induces masculinization of XX medaka through GSDF and AMHR2. This study thereby provides the first evidence that GSDF and AMHR2 are involved in cortisol-induced masculinization.


Introduction
Many vertebrates reproduce sexually and have separate male and female sexes determined by genotype; however, sex determination in poikilothermic vertebrates-such as reptiles, amphibians, and fish-is also highly influenced by ambient temperature (Baroiller et al. 1999). Furthermore, in some teleost fishes, sex can be affected by environmental variables such as pH (Rubin 1985), density (Francis 1984), and social factors (Francis and Barlow 1933). The molecular mechanisms of environmental sex determination in such species are not well understood.
The medaka Oryzias latipes is a small gonochoristic teleost fish that offers many advantages for researchers, including the availability of a range of useful strains (Ishikawa 2000). Applications of transgenic techniques and gene knockout (KO) systems, using transcription activator-like effector nucleases (TALENs) and clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9), have also been established in medaka (Ozato et al. 1986;Ansai et al. 2013;Sawamura et al. 2017). Medaka has an XX/XY sex determination system, and the sex-determining gene dmy (also known as dmrt1bY) has been identified on the medaka Y chromosome (Matsuda et al. 2002). Medaka is therefore a suitable model vertebrate for the genetic analysis of sex determination and differentiation. The first appearance of morphological sex differentiation in medaka is a difference in the number of germ cells between sexes prior to hatching. Germ cells in genetic females (XX) proliferate rapidly and subsequently initiate oogenesis, while germ cells in genetic males (XY) remain quiescent (Satoh and Egami 1972;Kobayashi 1 3 et al. 2004). Notably, exposure to high-temperature (HT) treatment during sex differentiation inhibits germ cell proliferation and induces masculinization in XX medaka (Sato et al. 2005;Hattori et al. 2007;Selim et al. 2009). HT also induces masculinization by elevating the levels of cortisol, which is the major glucocorticoid in teleost fish (Hayashi et al. 2010). In fish, cortisol is secreted from interrenal tissues by environmental stimuli such as acidic water and rapid temperature changes (Wendelaar Bonga 1997). The frequencies of XX males caused by cortisol treatment varies depending on experimental conditions including treatment periods and medaka strains (21.9-95.7%) (Hayashi et al. 2010;Kitano et al. 2012;Adolfi et al. 2019;Hara et al. 2020). Moreover, exposure to cortisol or HT in XX medaka is also known to elevate the expression of gonadal soma-derived growth factor (gsdf) by 5 days posthatching (dph) (Kitano et al. 2012); however, the exact role of HT-elevated cortisol in the masculinization of XX medaka remains unknown.
GSDF plays an essential role in the development of both the testis and the ovary, as demonstrated by gain-and lossof-function experiments in medaka (Zhang et al. 2016). The expression of gsdf co-localizes with dmy, and gsdf is regulated by dmy during the initiation of morphological testicular differentiation (Shibata et al. 2010;Zhang et al. 2016). Importantly, gsdf depletion leads to an excessive proliferation of germ cells during sexual differentiation regardless of sex, and some XY gonads in gsdf mutants develop into ovaries at the adult stage (Imai et al. 2015;Zhang et al. 2016). These phenotypes resemble those of anti-Müllerian hormone receptor 2 (amhr2) KO medaka (Morinaga et al. 2007). Anti-Müllerian hormone (AMH) signalling regulates proliferation at a specific stage of germ cell development, and this regulation is crucial for the proper manifestation of gonadal sex directed by sex determination genes (Nakamura et al. 2012;Pfennig et al. 2015). Thus, the importance of gsdf and amhr2 in sexual development has been established; however, our understanding of how these genes function in cortisol-induced masculinization is limited.
In this study, we investigated whether GSDF and AMHR2 are involved in masculinization by cortisol. First, we generated gsdf KO medaka using the CRISPR/Cas9 system. Then, we investigated the number of germ cells at the hatching stage in the gsdf KO medaka, and in previously generated amhr2 KO medaka (Furukawa et al. 2019), following treatment with cortisol. We also examined the sex ratios of gsdf KO and amhr2 KO adult fish by histological observation of the gonads after cortisol treatment.

Ethics statement
The study was performed using protocols approved by the Animal Care and Use Committee of Kumamoto University (Approval number: 30-022). All experiments were performed in accordance with the relevant guidelines and regulations.

Animals
Female leucophore-free (FLFII) medaka stock was used (Wakamatsu et al. 2003); this stock allows the identification of genotypic sex before the onset of sex differentiation, via the appearance of leucophores at 2 days post-fertilization (dpf). Fish embryos and larvae were maintained in embryo rearing medium (ERM) (17 mM NaCl, 0.4 mM KCl, 0.27 mM CaCl 2 2H 2 O, 0.66 mM MgSO 4 , pH 7), at a water temperature of 26 °C, with a 14 h light and 10 h dark cycle.
For genotyping, genomic DNA was extracted from adult caudal fins, as described previously (Hayashi et al. 2010). Genetic sex was then determined by genomic polymerase chain reaction (PCR) using primers specific for dmy, as also described previously (Hayashi et al. 2010). The gsdf genotypes were determined by PCR using primers specific for gsdf (5′-CAT GGG TGT CAC GGA CAC TAC-3′ and 5′-ATC AGC CTC ACA GGT GTT TG-3′) and AmpliTaq Gold (Applied Biosystems). The PCR conditions were as follows: preheating at 95 °C for 10 min, 40 cycles of PCR at 94 °C for 30 s, 59 °C for 30 s, 72 °C for 1 min, and a final extension at 72 °C for 5 min. PCR products were subcloned into the pT7Blue vector (Novagen) and sequenced with the Genom-eLab GeXP™ Genetic Analysis System (Beckman Coulter, Brea, CA, USA).

Histological analysis
Fish at 0 dph and 2 mph were fixed in Bouin's solution overnight at 4 °C, then dehydrated, embedded in paraffin, and serially sectioned at 5 μm, as described previously (Hayashi et al. 2010). Sections were stained with haematoxylin and eosin as also described previously (Hayashi et al. 2010). Germ cells were counted using an MZ FLIII microscope (Leica Microsystems).

Statistical analysis
Experimental results were tested using Levene's test for homogeneity of variance. Data were analysed by Student's t-test or by one-way ANOVA followed by Tukey's multiple comparison test using SPSS statistics 20 (IBM Corp., Armonk, NY).

Generation of gsdf KO medaka
We knocked out medaka gsdf using the CRISPR/Cas9 system (Sawamura et al. 2017). As shown in Fig. 1, deletion of the transforming growth factor-β (TGF-β) domain in gsdf was achieved by simultaneous injection of two crRNAs, a tracrRNA, and the Cas9 protein into 1-cell stage embryos. After mating injected medaka with wild-type fish, gsdf genotypes in the F1 generation were determined by sequence analysis of PCR products amplified from genomic DNA extracted from adult caudal fins. The gsdf KO fish were obtained by crossing F1 heterozygous mutants.

Number of germ cells at the hatching stage in cortisol-treated gsdf KO and amhr2 KO medaka
We confirmed the effects of cortisol on gsdf KO and amhr2 KO medaka (Furukawa et al. 2019). First, we counted the number of germ cells at the hatching stage in medaka treated with cortisol from 0 dpf. In wild-type fish, XX individuals treated with cortisol had around 50 germ cells, which corresponds to the number of germ cells in males. In contrast, there were high numbers of germ cells in gsdf KO medaka and amhr2 KO medaka, in both XX and XY individuals,  with and without cortisol treatment (Fig. 2); gsdf KO XX medaka with and without cortisol treatment had significantly more germ cells than wild-type XX fish at the hatching stage.

Sex ratios in cortisol-treated gsdf KO and amhr2 KO medaka
We examined the sex ratios of adult fish by histological observation of the gonads at 2 mph after treatment with cortisol during the sexual differentiation period. In wild-type fish, cortisol treatment induced masculinization of XX individuals in some cases (Fig. 3a-e, Table 1). Conversely, gsdf KO medaka and amhr2 KO medaka were not masculinized by cortisol. In addition, under non-cortisol conditions, ovaries were found in several XY individuals of the KO medaka ( Fig. 3f-o, Table 1). Histological examination revealed testes or ovaries only; there were no cases of intersex gonads.

Discussion
To investigate whether gsdf is essential for cortisolinduced masculinization, we knocked out the TGF-β domain of medaka gsdf using the CRISPR/Cas9 system (Sawamura et al. 2017). The gsdf KO medaka underwent female-type proliferation of germ cells, and 7.1% of XY gsdf KO fish differentiated into females. In addition, gsdf KO females had larger ovaries than wild-type females and were fertile (data not shown). As these data are similar to previously reported phenotypes of gsdf KO medaka (Imai et al. 2015;Zhang et al. 2016;Guan et al. 2017) in which the location of the gsdf mutation differed from that in our study, we inferred that the gsdf KO line generated here also lacks gsdf function.
In this study, we showed that cortisol induces masculinization of wild-type XX medaka, but not of gsdf KO and amhr2 KO XX fish, suggesting important roles of GSDF and AMHR2 in cortisol-induced masculinization. The frequency of wild-type XX males induced by cortisol in this study was lower than that in previous reports (Hayashi et al. 2010;Kitano et al. 2012;Adolfi et al. 2019;Hara et al. 2020). This may be due to the death of the masculinized fish during breeding, because the final survival rate (48.9%) of the cortisol-treated individuals in this study was lower than that (59.0%) in a recent study (Hara et al. 2020). On the other hand, cortisol inhibited female-type proliferation of germ cells at hatching stage in wild-type XX medaka, but not in gsdf KO and amhr2 KO XX fish. This suggests that cortisol inhibits germ cell proliferation via gsdf and amhr2. These results indicate that GSDF and AMHR2 are indispensable for cortisol-induced masculinization and male-type proliferation of germ cells.
All gsdf KO and amhr2 KO XY medaka showed a femaletype proliferation of germ cells at hatching stage. However, as also seen previously, some gsdf KO and amhr2 KO XY medaka developed into males by adulthood, under non-cortisol treatment (Imai et al. 2015;Zhang et al. 2016;Morinaga et al. 2007). It is thought that doublesex-and mab-3-related transcription factor 1 (dmrt1), which acts downstream of gsdf, affects the differentiation of gsdf KO XY medaka into males (Imai et al. 2015). It was recently reported that the expression of dmrt1 was dramatically down-regulated in gsdf KO XY gonads compared with wild-type XY gonads, indicating that GSDF may directly or indirectly activate or maintain its expression (Zhang et al. 2016). Moreover, exposure to cortisol or HT in XX medaka induced dmrt1 expression during sex differentiation (Adolfi et al. 2019). Therefore, cortisol appears to activate and maintain masculinization of XX medaka by inducing the expression of dmrt1 via GSDF.
Our previous study provided evidence that peroxisome proliferator-activated receptor α (PPARα) is involved in masculinization by cortisol (Hara et al. 2020). Interestingly, PPARα agonist treatment induced masculinization of XX medaka in some cases. In addition, pparaa KO XX medaka were not masculinized by cortisol, as observed in gsdf KO and amhr2 KO medaka in this study, suggesting that activation of PPARα may directly or indirectly regulate gsdf and amhr2 expression (Fig. 4). Currently, little is known about the relationship between PPARα and male-promoting genes.
To elucidate the masculinization mechanism, it will be necessary to investigate whether PPARα directly induces the expression of gsdf and amhr2, for example, by the analysis of DNA-protein complexes via chromatin immunoprecipitation assays.
In summary, this study provides the first evidence that GSDF and AMHR2 are involved in cortisol-induced masculinization in medaka. Elucidation of the molecular mechanisms by which this stress hormone affects sex is expected to greatly contribute to artificial sex control. In medaka, gsdf and amhr2 are located on autosomal chromosomes, but in some fish species these genes are sex-determining genes: for example, gsdf in Oryzias luzonensis (Myosho et al. 2012), and amhr2 in Takifugu rubripes (Kamiya et al. 2012). This suggests that genetic female fish of such species have lost the masculinization functions of gsdf or amhr2. Furthermore, if gsdf or amhr2 are essential for masculinization by cortisol in fish, then cortisol may not be able to induce masculinization of genetic females in these species. Future studies will therefore focus on the effects of environmental stress or cortisol on fish species in which gsdf and amhr2 are sexdetermining genes.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.