Bmp7 and Lef1 are the downstream effectors of androgen signaling in androgen-induced sex characteristics development in medaka

Inconspicuous breeding coloration to conceal eggs during mouthbrooding in male cardinalfish

Animals exhibit colorations optimal for their niche, which hides their existence from other organisms. In Apogoninae fishes, the father broods their egg inside their mouth. Since the color of eggs is different from parental fish, it can disrupt the optimal camouflage coloration of parental fish if the lower jaw is transparent. Here, we identified male- and breeding season-specific whitish coloration consisting of iridophores in the lower jaw. Artificial implantation of eggs inside the mouth of females and males showed that iridophores in the lower jaws concealed the conspicuous coloration of eggs only in males. In addition, it was revealed that iridophore development in the lower jaw is induced by androgen through the Alkal-Ltk pathway. These results suggest that androgen-dependent breeding colorations in males, which have been considered to attract females, may serve the opposite function, "inconspicuous breeding coloration" in these species.

Toxic effects of chronic exposure to BPAF and perturbation of gut microbiota homeostasis in marine medaka (Oryzias melastigma)

Bisphenol AF (BPAF), a substitute for bisphenol A (BPA), exhibits potent endocrine-disrupting properties that pose a serious health hazard to organisms. This study employed marine medaka as a model, subjecting them to different concentrations of BPAF (0.61, 6.65, and 91.88 μg/L) from the embryonic stage for a period of 160 days. Findings showed that 91.88 μg/L BPAF reduced survival rates and altered sex ratios. Furthermore, exposure to BPAF at all concentrations led to a significant increase in body length and weight. Behavioral analysis revealed that BPAF exposure impaired the swimming ability of the medaka. Histological changes included disrupted ovarian development, reduced sperm count, liver inflammation, and intestinal damage. Gene expression analysis revealed impacts on nervous system (e.g., gap43, itr, elavl3), HPG axis (e.g., gthα, erα, 3βhsd), and liver genes (e.g., chgl, vtg2). Additionally, BPAF altered the diversity and richness of gut microbes in marine medaka, leading to significant changes in specific bacterial species and intestinal functions. In conclusion, long-term BPAF exposure induced neurotoxicity, reproductive toxicity, and impaired digestive and immune systems in marine medaka, with sex-specific effects. These results provide further evidence of the potential hazards of BPAF as an environmental pollutant.

AOP Report: An Upstream Network for Reduced Androgen Signaling Leading to Altered Gene Expression of Androgen Receptor-Responsive Genes in Target Tissues

Adverse outcome pathways (AOPs) can aid with chemical risk assessment by providing plausible links between chemical activity at the molecular level and effect outcomes in intact organisms. Because AOPs can be used to infer causality between upstream and downstream events in toxicological pathways, the AOP framework can also facilitate increased uptake of alternative methods and new approach methodologies to help inform hazard identification. However, a prevailing challenge is the limited number of fully developed and endorsed AOPs, primarily due to the substantial amount of work required by AOP developers and reviewers. Consequently, a more pragmatic approach to AOP development has been proposed where smaller units of knowledge are developed and reviewed independent of full AOPs. In this context, we have developed an upstream network comprising key events (KEs) and KE relationships related to decreased androgen signaling, converging at a nodal KE that can branch out to numerous adverse outcomes (AOs) relevant to androgen-sensitive toxicological pathways. Androgen signaling represents an extensively studied pathway for endocrine disruption. It is linked to numerous disease outcomes and can be affected by many different endocrine-disrupting chemicals. Still, pathways related to disrupted androgen signaling remain underrepresented in the AOP-wiki, and endorsed AOPs are lacking. Given the pivotal role of androgen signaling in development and function across vertebrate taxa and life stages of both sexes, this upstream AOP network serves as a foundational element for developing numerous AOPs. By connecting the upstream network with various downstream AOs, encompassing different species, it can also facilitate cross-species extrapolations for hazard and risk assessment of chemicals. Environ Toxicol Chem 2024;43:2329-2337. © 2024 The Author(s). Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

De novo assembly and characterization of a highly degenerated ZW sex chromosome in the fish Megaleporinus macrocephalus

BACKGROUND

Megaleporinus macrocephalus (piauçu) is a Neotropical fish within Characoidei that presents a well-established heteromorphic ZZ/ZW sex determination system and thus constitutes a good model for studying W and Z chromosomes in fishes. We used PacBio reads and Hi-C to assemble a chromosome-level reference genome for M. macrocephalus. We generated family segregation information to construct a genetic map, pool sequencing of males and females to characterize its sex system, and RNA sequencing to highlight candidate genes of M. macrocephalus sex determination.

RESULTS

The reference genome of M. macrocephalus is 1,282,030,339 bp in length and has a contig and scaffold N50 of 5.0 Mb and 45.03 Mb, respectively. In the sex chromosome, based on patterns of recombination suppression, coverage, FST, and sex-specific SNPs, we distinguished a putative W-specific region that is highly differentiated, a region where Z and W still share some similarities and is undergoing degeneration, and the PAR. The sex chromosome gene repertoire includes genes from the TGF-β family (amhr2, bmp7) and the Wnt/β-catenin pathway (wnt4, wnt7a), some of which are differentially expressed.

CONCLUSIONS

The chromosome-level genome of piauçu exhibits high quality, establishing a valuable resource for advancing research within the group. Our discoveries offer insights into the evolutionary dynamics of Z and W sex chromosomes in fish, emphasizing ongoing degenerative processes and indicating complex interactions between Z and W sequences in specific genomic regions. Notably, amhr2 and bmp7 are potential candidate genes for sex determination in M. macrocephalus.

Genome-wide analysis of BMP/GDF family and DAP-seq of YY1 suggest their roles in Cynoglossus semilaevis sexual size dimorphism

Sexual size dimorphism (SSD) characterized by different body size between females and males have been reported in various animals. Gonadectomy experiments have implied important regulatory roles of the gonad in SSD. Among multiple factors from the gonad, TGF-β superfamily (especially BMP/GDF family) attracted our interest due to its pleiotropy in growth and reproduction regulations. Thus, whether BMP/GDF family members serve as crucial regulators for SSD was studied in a typically female-biased SSD flatfish named Chinese tongue sole (Cynoglossus semilaevis). Firstly, a total of 26 BMP/GDF family members were identified. The PPI network analysis showed that they may interact with ACVR2a, ACVR2b, ACVR1, BMPR2, SMAD3, BMPR1a, and other proteins. Subsequently, DAP-seq was employed to reveal the binding sites for yin yang 1 (yy1), a transcription factor involved in gonad function and cell growth partly by regulating TGF-β superfamily. The results revealed that two yy1 homologues yy1a and yy1b in C. semilaevis could regulate Hippo signaling pathway, mTOR signaling pathway, and AMPK signaling pathway. Moreover, BMP/GDF family genes including bmp2, bmp4, bmp5, gdf6a, and gdf6b were important components of Hippo pathway. In future, the crosstalk among yy1a, yy1b, and TGF-β family would provide more insight into sexual size dimorphism in C. semilaevis.

A systematic review of the evaluation of endocrine-disrupting chemicals in the Japanese medaka (Oryzias latipes) fish

Japanese medaka (Oryzias latipes) is an acceptable small laboratory fish model for the evaluation and assessment of endocrine-disrupting chemicals (EDCs) found in the environment. In this research, we used this fish as a potential tool for the identification of EDCs that have a significant impact on human health. We conducted an electronic search in PubMed (http://www.ncbi.nlm.nih.gov/pubmed) and Google Scholar (https://scholar.google.com/) using the search terms, Japanese medaka, Oryzias latipes, and endocrine disruptions, and sorted 205 articles consisting of 128 chemicals that showed potential effects on estrogen-androgen-thyroid-steroidogenesis (EATS) pathways of Japanese medaka. From these chemicals, 14 compounds, namely, 17β-estradiol (E2), ethinylestradiol (EE2), tamoxifen (TAM), 11-ketotestosterone (11-KT), 17β-trenbolone (TRB), flutamide (FLU), vinclozolin (VIN), triiodothyronine (T3), perfluorooctanoic acid (PFOA), tetrabromobisphenol A (TBBPA), terephthalic acid (TPA), trifloxystrobin (TRF), ketoconazole (KTC), and prochloraz (PCZ), were selected as references and used for the identification of apical endpoints within the EATS modalities. Among these endpoints, during classification, priorities are given to sex reversal (masculinization of females and feminization of males), gonad histology (testis-ova or ovotestis), secondary sex characteristics (anal fin papillae of males), plasma and liver vitellogenin (VTG) contents in males, swim bladder inflation during larval development, hepatic vitellogenin (vtg) and choriogenin (chg) genes in the liver of males, and several genes, including estrogen-androgen-thyroid receptors in the hypothalamus-pituitary-gonad/thyroid axis (HPG/T). After reviewing 205 articles, we identified 108 (52.68%), 46 (22.43%), 19 (9.26%), 22 (17.18%), and 26 (12.68%) papers that represented studies on estrogen endocrine disruptors (EEDs), androgen endocrine disruptors (AEDs), thyroid endocrine disruptors (TEDs), and/or steroidogenesis modulators (MOS), respectively. Most importantly, among 128 EDCs, 32 (25%), 22 (17.18%), 15 (11.8%), and 14 (10.93%) chemicals were classified as EEDs, AEDs, TEDs, and MOS, respectively. We also identified 43 (33.59%) chemicals as high-priority candidates for tier 2 tests, and 13 chemicals (10.15%) show enough potential to be considered EDCs without any further tier-based studies. Although our literature search was unable to identify the EATS targets of 45 chemicals (35%) studied in 60 (29.26%) of the 205 articles, our approach has sufficient potential to further move the laboratory-based research data on Japanese medaka for applications in regulatory risk assessments in humans.

17β-Trenbolone activates androgen receptor, upregulates transforming growth factor beta/bone morphogenetic protein and Wnt signaling pathways, and induces masculinization of caudal and anal fins in female guppies (Poecilia reticulata)

Sexually mature female guppies (Poecilia reticulata) were exposed to environmentally relevant concentrations (20, 200, and 2000 ng/L) of 17β-trenbolone for four weeks. As evidenced by the increased caudal fin index and anal fins developing into gonopodium-like structures, exposed females displayed masculinized secondary sexual characteristics. Differential gene expression and subsequent pathway analysis of mRNA sequencing data revealed that the transcription of transforming growth factor beta/bone morphogenetic protein signaling pathway and Wnt signaling pathway were upregulated following 17β-trenbolone exposure. Enzyme-linked immunosorbent assays showed that the bone morphogenetic protein 7 protein content was elevated after 17β-trenbolone exposure. Finally, real-time PCR revealed that 17β-trenbolone treatment significantly increased androgen receptor mRNA levels, and molecular docking showed potent interaction between 17β-trenbolone and guppy androgen receptor. Furthermore, 17β-trenbolone-induced masculinization of caudal and anal fins in female guppies, concomitant to the upregulated expression of differentially expressed genes involved in the above-mentioned two signaling pathways, was significantly inhibited by flutamide (androgen receptor antagonist). These findings demonstrated that 17β-trenbolone masculinized fins of female guppies by activating the androgen receptor. This study revealed that 17β-trenbolone could upregulate signaling pathways related to fin growth and differentiation, and eventually cause caudal and anal fin masculinization in female guppies.

Evolutionary differentiation of androgen receptor is responsible for sexual characteristic development in a teleost fish

Teleost fishes exhibit complex sexual characteristics in response to androgens, such as fin enlargement and courtship display. However, the molecular mechanisms underlying their evolutionary acquisition remain largely unknown. To address this question, we analyse medaka (Oryzias latipes) mutants deficient in teleost-specific androgen receptor ohnologs (ara and arb). We discovered that neither ar ohnolog was required for spermatogenesis, whilst they appear to be functionally redundant for the courtship display in males. However, both were required for reproductive success: ara for tooth enlargement and the reproductive behaviour eliciting female receptivity, arb for male-specific fin morphogenesis and sexual motivation. We further showed that differences between the two ar ohnologs in their transcription, cellular localisation of their encoded proteins, and their downstream genetic programmes could be responsible for the phenotypic diversity between the ara and arb mutants. These findings suggest that the ar ohnologs have diverged in two ways: first, through the loss of their roles in spermatogenesis and second, through gene duplication followed by functional differentiation that has likely resolved the pleiotropic roles derived from their ancestral gene. Thus, our results provide insights into how genome duplication impacts the massive diversification of sexual characteristics in the teleost lineage.

Genome-wide identification, evolution and expression analysis of bone morphogenetic protein (BMP) gene family in chinese soft-shell turtle (Pelodiscus sinensis)

Introduction: Bone morphogenetic proteins (BMPs) play a crucial role in bone formation and differentiation. Recent RNA-Seq results suggest that BMPs may be involved in the sex differentiation of P. sinensis, yet more relevant studies about BMPs in P. sinensis are lacking. Methods: Herein, we identified BMP gene family members, analyzed the phylogeny, collinear relationship, scaffold localization, gene structures, protein structures, transcription factors and dimorphic expression by using bioinformatic methods based on genomic and transcriptomic data of P. sinensis. Meanwhile, qRT-PCR was used to verify the RNA-Seq results and initially explore the function of the BMPs in the sex differentiation of P. sinensis. Results: A total of 11 BMP genes were identified, 10 of which were localized to their respective genomic scaffolds. Phylogenetic analysis revealed that BMP genes were divided into eight subfamilies and shared similar motifs ("WII", "FPL", "TNHA", "CCVP", and "CGC") and domain (TGF-β superfamily). The results of the sexually dimorphic expression profile and qRT-PCR showed that Bmp2, Bmp3, Bmp15l, Bmp5, Bmp6 and Bmp8a were significantly upregulated in ovaries, while Bmp2lb, Bmp7, Bmp2bl and Bmp10 were remarkable upregulated in testes, suggesting that these genes may play a role in sex differentiation of P. sinensis. Discussion: Collectively, our comprehensive results enrich the basic date for studying the evolution and functions of BMP genes in P. sinensis.

Bisphenol A alters sexual dimorphism and gene expression in marine medaka Oryzias melastigma

Bisphenol A (BPA) is an endocrine disruptor that is present in freshwater and marine environments. However, conclusive evidence for the toxicity of chronic BPA exposure to marine fishes remains lacking. Therefore, we investigated the influence of BPA on male marine medaka (Oryzias melastigma). BPA exposure induced formation of testis-ova at 2610 µg/L, and male-type anal fins became more female type in a concentration-dependent manner. Some males with female-type anal fins had normal testes, indicating that anal fin shape is more sensitive to BPA. Gonadal soma-derived factor (gsdf) expression decreased after BPA exposure in the 746 and 2610 µg/L exposure groups, although the changes were not statistically significant. Additionally, liver vitellogenin (vtg) expression increased in a dose-dependent manner and was significantly higher in all exposure groups. vtg and gsdf are likely to be useful biomarkers for the impact of estrogenic endocrine disrupters in O. melastigma.

Summary of reference chemicals evaluated by the fish short-term reproduction assay, OECD TG229, using Japanese Medaka, Oryzias latipes

Under the Organisation for Economic Co-operation and Development (OECD), the Ministry of the Environment of Japan (MOE) added Japanese medaka (Oryzias latipes) to the test guideline fish short-term reproduction assay (FSTRA) developed by the United States Environmental Protection Agency (US EPA) using fathead minnow (Pimephales promelas). The FSTRA was designed to detect endocrine disrupting effects of chemicals interacting with the hypothalamic-pituitary-gonadal axis (HPG axis) such as agonists or antagonists on the estrogen receptor (Esr) and/or the androgen receptor (AR) and steroidogenesis inhibitors. We conducted the FSTRA with Japanese medaka, in accordance with OECD test guideline number 229 (TG229), for 16 chemicals including four Esr agonists, two Esr antagonists, three AR agonists, two AR antagonists, two steroidogenesis inhibitors, two progesterone receptor agonists, and a negative substance, and evaluated the usability and the validity of the FSTRA (TG229) protocol. In addition, in vitro reporter gene assays (RGAs) using Esr1 and ARβ of Japanese medaka were performed for the 16 chemicals, to support the interpretation of the in vivo effects observed in the FSTRA. In the present study, all the test chemicals, except an antiandrogenic chemical and a weak Esr agonist, significantly reduced the reproductive status of the test fish, that is, fecundity or fertility, at concentrations where no overt toxicity was observed. Moreover, vitellogenin (VTG) induction in males and formation of secondary sex characteristics (SSC), papillary processes on the anal fin, in females was sensitive endpoints to Esr and AR agonistic effects, respectively, and might be indicators of the effect concentrations in long-term exposure. Overall, it is suggested that the in vivo FSTRA supported by in vitro RGA data can adequately detect effects on the test fish, O. latipes, and probably identify the mode of action (MOA) of the chemicals tested.

Genome editing reveals fitness effects of a gene for sexual dichromatism in Sulawesian fishes

Sexual selection drives rapid phenotypic diversification of mating traits. However, we know little about the causative genes underlying divergence in sexually selected traits. Here, we investigate the genetic basis of male mating trait diversification in the medaka fishes (genus Oryzias) from Sulawesi, Indonesia. Using linkage mapping, transcriptome analysis, and genome editing, we identify csf1 as a causative gene for red pectoral fins that are unique to male Oryzias woworae. A cis-regulatory mutation enables androgen-induced expression of csf1 in male fins. csf1-knockout males have reduced red coloration and require longer for mating, suggesting that coloration can contribute to male reproductive success. Contrary to expectations, non-red males are more attractive to a predatory fish than are red males. Our results demonstrate that integrating genomics with genome editing enables us to identify causative genes underlying sexually selected traits and provides a new avenue for testing theories of sexual selection.

Exposure to 4-nonylphenol induces a shift in the gene expression of gsdf and testis-ova formation and sex reversal in Japanese medaka (Oryzias latipes)

The branched isomer mixture 4-nonylphenol (4-NP) has been used worldwide as a surfactant, and can have endocrine-disrupting effects on aquatic organisms. For instance, 4-NP induces the formation of testis-ova (i.e., testicular and ovarian tissue in the same gonad) or male to female sex reversal of various teleost fishes. Recently, our group revealed that altered gsdf gene expression is associated with disruption of gonadal differentiation in Japanese medaka (Oryzias latipes) embryos exposed to methyltestosterone or bisphenol A, suggesting that gsdf might be useful as a biomarker for predicting the impact of endocrine-disrupting chemicals (EDCs) on gonadal differentiation. Here, we used 4-NP to examine further whether gsdf expression at the embryo stage is useful for predicting EDC impact on gonadal sex differentiation. When fertilized medaka eggs were exposed to 32 or 100 μg/L 4-NP, testis-ova in genetic males and sex reversal from genetic male to phenotypic female were observed. At stage 38 (just before hatching), 4-NP exposure at 1-100 μg/L did not affect gsdf expression in XX embryos compared with the nontreated control; however, in XY embryos, the gsdf expression in the 100 μg/L-exposed group was significantly lower than that in the controls. The 4-NP concentration at which gsdf expression was suppressed was equal to that at which testis-ova and sex reversal were induced. These results indicate that expression of the gsdf gene at the embryonic stage in medaka is a useful biomarker for predicting the impact of EDCs on sexual differentiation.

Functional distinctions associated with the diversity of sex steroid hormone receptors ESR and AR

Sex steroid hormones including estrogens and androgens play fundamental roles in regulating reproductive activities and they act through estrogen and androgen receptors (ESR and AR). These steroid receptors have evolved from a common ancestor in association with several gene duplications. In most vertebrates, this has resulted in two ESR subtypes (ESR1 and ESR2) and one AR, whereas in teleost fish there are at least three ESRs (ESR1, ESR2a and ESR2b) and two ARs (ARα and ARβ) due to a lineage-specific whole genome duplication. Functional distinctions have been suggested among these receptors, but to date their roles have only been characterized in a limited number of species. Sexual differentiation and the development of reproductive organs are indispensable for all animal species and in vertebrates these events depend on the action of sex steroid hormones. Here we review the recent progress in understanding of the functions of the ESRs and ARs in the development and expression of sexually dimorphic characteristics associated with steroid hormone signaling in vertebrates, with representative fish, amphibians, reptiles, birds and mammals.

Summary of the development the US Environmental Protection Agency's Medaka Extended One Generation Reproduction Test (MEOGRT) using data from 9 multigenerational medaka tests

In response to various legislative mandates, the US Environmental Protection Agency (USEPA) formed its Endocrine Disruptor Screening Program (EDSP), which in turn, formed the basis of a tiered testing strategy to determine the potential of pesticides, commercial chemicals, and environmental contaminants to disrupt the endocrine system. The first tier of tests is intended to detect the potential for endocrine disruption mediated through estrogen, androgen, or thyroid pathways, whereas the second tier is intended to further characterize the effects on these pathways and to establish a dose-response relationship for adverse effects. One of these tier 2 tests, the Medaka Extended One Generation Reproduction Test (MEOGRT), was developed by the USEPA for the EDSP and, in collaboration with the Japanese Ministry of the Environment, for the Guidelines for the Testing of Chemicals of the Organisation for Economic Co-operation and Development (OECD). The MEOGRT protocol was iteratively modified based on knowledge gained after the successful completion of 9 tests with variations in test protocols. The present study describes both the final MEOGRT protocol that has been published by the USEPA and the OECD, and the iterations that provided valuable insights into nuances of the protocol. The various tests include exposure to 17β-estradiol, 4-t-octylphenol, o,p'- dichlorodiphenyltrichloroethane, 4-chloro-3-methylphenol, tamoxifen, 17β-trenbolone, vinclozolin, and prochloraz. Environ Toxicol Chem 2017;36:3387-3403. Published 2017 Wiley Periodicals Inc. on behalf of SETAC. This article is a US government work and, as such, is in the public domain in the United States of America.

Cloning, expression profiling and promoter functional analysis of Bone morphogenetic protein 6 and 7 in tongue sole (Cynoglossus semilaevis)

Bone morphogenetic proteins (BMPs) play crucial roles in vertebrate developmental process and are associated with the mechanisms which drive early skeletal development. As a first approach to elucidating the role of BMPs in regulating fish bone formation and growth, we describe the cloning, expression profiling and promoter functional analysis of bmp6 and bmp7 in tongue sole (Cynoglossus semilaevis). The full length of bmp6 and bmp7 cDNA sequences is 1939 and 1836 bp, which encodes a protein of 428 and 427 amino acids, respectively. Tissue expression distribution of bmp6 and bmp7 was examined in 14 tissues of mature individuals by quantitative real-time PCR (qRT-PCR). The results revealed that bmp6 was predominantly expressed in the gonad, and bmp7 exhibited the highest expression level in the dorsal fin. Further comparison of bmp6 expression levels between female and male gonads showed that the expression in the ovary was significantly higher than in the testis. Moreover, bmp6 and bmp7 expression levels were detected at 15 sampling time points of early developmental stages (egg, larva, juvenile and fingerling stages). The highest expression level of bmp6 was observed in the egg stage (multi-cell and gastrula stage); while bmp7 exhibited the highest expression in the larva stage (1-4 days old). The high expression levels of BMP6 in the ovary as well as at early embryonic stages indicated that the maternally stored transcripts of bmp6 might play a role in early embryonic development. Whole-mount in situ hybridization showed that bmp6 and bmp7 exhibited similar spatial expression patterns. Both bmp6 and bmp7 signals were first detected in the head and anterior regions in newly hatched larvae, and then, the mRNAs appeared in the crown-like larval fin, jaw, operculum and fins (pectoral, dorsal, pelvic and anal) along with early development. Subsequently, we characterized the 5'-flanking regions of bmp6 and bmp7 by testing the promoter activity by luciferase reporter assays. Positive regulatory regions were, respectively, detected at the location of -272 to +28 and -740 to -396 in bmp6 and bmp7 gene. The predicted transcription factor binding sites (CREB, AP1 and methyl-CpG-binding protein) in the regions might participate in the transcriptional regulation of these two genes.

Diagnostic Application of Targeted Next-Generation Sequencing of 80 Genes Associated with Disorders of Sexual Development

Disorders of sexual development (DSD) are estimated to occur in 1 of 4500 births. Since the genetic etiology of DSD is highly heterogeneous, obtaining a definitive molecular diagnosis by single gene test is challenging. Utilizing a high-throughput sequencing upfront is proposed as an efficient approach to aid in the diagnosis. This study aimed to examine the diagnostic yield of next-generation sequencing in DSD. 32 DSD patients that previously received clinical examinations and single gene tests were selected, with or without a diagnosis. Prior single gene tests were masked, and then samples went through targeted next-generation sequencing of 80 genes from which the diagnostic yield was assessed. A likely diagnosis, with pathogenic or likely pathogenic variants identified, was obtained from nine of the 32 patients (i.e., 28.1%, versus 10% by single gene tests). In another five patients (15.6%), variants of uncertain significance were found. Among 18 variants identified (i.e., 17 single nucleotide variants and one small deletion), eight had not been previously reported. This study supports the notion that next-generation sequencing can be an efficient tool in the clinical diagnosis and variant discovery in DSD.

Serotonergic neuronal death and concomitant serotonin deficiency curb copulation ability of Drosophila platonic mutants

Drosophila platonic (plt) males court females, but fail to copulate. Here we show that plt is an allele of scribbler (sbb), a BMP signalling component. sbb knockdown in larvae leads to the loss of approximately eight serotonergic neurons, which express the sex-determinant protein Doublesex (Dsx). Genetic deprivation of serotonin (5-HT) from dsx-expressing neurons results in copulation defects. Thus, sbb⁺ is developmentally required for the survival of a specific subset of dsx-expressing neurons, which support the normal execution of copulation in adults by providing 5-HT. Our study highlights the conserved involvement of serotonergic neurons in the control of copulatory mechanisms and the key role of BMP signalling in the formation of a sex-specific circuitry.

Drosophila platonic (plt) mutant males court with females but fail to copulate. Here, the authors find plt is an allele of scribbler and may disrupt courtship behaviour via developmental disruption of a subgroup of serotonergic Doublesex+ neurons in the abdominal ganglion.

Polygenic sex determination in the cichlid fish Astatotilapia burtoni

BACKGROUND

The East African riverine cichlid species Astatotilapia burtoni serves as an important laboratory model for sexually dimorphic physiology and behavior, and also serves as an outgroup species for the explosive adaptive radiations of cichlid species in Lake Malawi and Lake Victoria. An astounding diversity of genetic sex determination systems have been revealed within the adaptive radiation of East African cichlids thus far, including polygenic sex determination systems involving the epistatic interaction of multiple, independently segregating sex determination alleles. However, sex determination has remained unmapped in A. burtoni. Here we present mapping results supporting the presence of multiple, novel sex determination alleles, and thus the presence of polygenic sex determination in A. burtoni.

RESULTS

Using mapping in small families in conjunction with restriction-site associated DNA sequencing strategies, we identify associations with sex at loci on linkage group 13 and linkage group 5-14. Inheritance patterns support an XY sex determination system on linkage group 5-14 (a chromosome fusion relative to other cichlids studied), and an XYW system on linkage group 13, and these associations are replicated in multiple families. Additionally, combining our genetic data with comparative genomic analysis identifies another fusion that is unassociated with sex, with linkage group 8-24 and linkage group 16-21 fused in A. burtoni relative to other East African cichlid species.

CONCLUSIONS

We identify genetic signals supporting the presence of three previously unidentified sex determination alleles at two loci in the species A. burtoni, strongly supporting the presence of polygenic sex determination system in the species. These results provide a foundation for future mapping of multiple sex determination genes and their interactions. A better understanding of sex determination in A. burtoni provides important context for their use in behavioral studies, as well as studies of the evolution of genetic sex determination and sexual conflicts in East African cichlids.

Interpretation of fish biomarker data for identification, classification, risk assessment and testing of endocrine disrupting chemicals

Chemical induced changes in fish biomarkers vitellogenin (VTG), secondary sex characteristics (SSC), and sex ratio indicate modes/mechanisms of action (MOAs) of EAS (estrogen, androgen and steroidogenesis) pathways. These biomarkers could be used for defining MOAs and the causal link between MOAs and adverse effects in fish for the identification of endocrine disrupting chemicals (EDCs). This paper compiled data sets of 150 chemicals for VTG, 57 chemicals for SSC and 38 chemicals for sex ratio in fathead minnow, medaka and zebrafish. It showed 1) changes in fish biomarkers can indicate the MOAs as anticipated; 2) in addition to EAS pathways, chemicals with non-EAS pathways induced changes in fish biomarkers; 3) responses of fish biomarkers did not always follow the anticipated patterns of EAS pathways. These responses may result from the interaction of chemical-induced multiple MOAs and confounding factors like fish diet, infection, culture conditions, general toxicity and stress response. The complex response of fish biomarkers to a chemical of interest requires EDC testing at multiple biological levels. Interpretation of fish biomarker data should be combined with relevant information at different biological levels, which is critical for defining chemical specific MOAs. The utility of fish biomarker data for identification, classification, PBT assessment, risk assessment, and testing of EDCs in the regulatory context was discussed. This paper emphasizes the importance of fish biomarker data in the regulatory context, a weight of evidence approach for the interpretation of fish biomarker data and the need for defining levels of evidence for the identification of EDCs.

Rspo1-activated signalling molecules are sufficient to induce ovarian differentiation in XY medaka (Oryzias latipes)

In contrast to our understanding of testicular differentiation, ovarian differentiation is less well understood in vertebrates. In mammals, R-spondin1 (Rspo1), an activator of Wnt/β-catenin signaling pathway, is located upstream of the female sex determination pathway. However, the functions of Rspo1 in ovarian differentiation remain unclear in non-mammalian species. In order to elucidate the detailed functions of Rspo/Wnt signaling pathway in fish sex determination/differentiation, the ectopic expression of the Rspo1 gene was performed in XY medaka (Oryzias latipes). The results obtained demonstrated that the gain of Rspo1 function induced femininity in XY fish. The overexpression of Rspo1 enhanced Wnt4b and β-catenin transcription, and completely suppressed the expression of male-biased genes (Dmy, Gsdf, Sox9a2 and Dmrt1) as well as testicular differentiation. Gonadal reprograming of Rspo1-over-expressed-XY (Rspo1-OV-XY) fish, induced the production of female-biased genes (Cyp19a1a and Foxl2), estradiol-17β production and further female type secondary sexuality. Moreover, Rspo1-OV-XY females were fertile and produced successive generations. Promoter analyses showed that Rspo1 transcription was directly regulated by DM domain genes (Dmy, the sex-determining gene, and Dmrt1) and remained unresponsive to Foxl2. Taken together, our results strongly suggest that Rspo1 is sufficient to activate ovarian development and plays a decisive role in the ovarian differentiation in medaka.

Variation in the ontogeny of sex steroid levels between latitudinal populations of the medaka

INTRODUCTION

Sex steroids mediate the expression of sexual dimorphism during ontogeny, and populations that differ in the magnitudes of sexual dimorphism may accordingly differ in the ontogenetic patterns of their sex steroid levels. The medaka, Oryzias latipes species complex, shows geographic variation in the magnitude of sexual dimorphism with respect to the lengths of their anal and dorsal fins; dimorphism is greater in low-latitude populations than in high-latitude populations. However, sexual differences in the ontogenetic dynamics of sex steroids, and its interpopulation variation, have not been examined.

RESULTS

We measured testosterone (T), estradiol-17β (E2), and 11-ketotestosterone (11-KT) concentrations throughout ontogeny of laboratory-reared fish from two latitudinal populations: Aomori (northern) and Okinawa (southern). In both populations, the levels of all three steroids were high during early ontogenetic stages and decreased with growth. After reaching about 15 mm in standard length, when sexual dimorphisms in fin lengths became apparent, steroid levels increased and tended to plateau. Sexual differences in the steroid levels were observed only in the later ontogenetic stages; T and 11-KT levels were higher in males, while E2 levels were higher in females. Accordingly, interpopulation differences also became clearer; the southern fish tended to show higher T levels and lower E2 levels than the northern fish.

CONCLUSIONS

The ontogenetic patterns of sex steroid levels paralleled the ontogeny of anal and dorsal fins in the two latitudinal populations, suggesting that interpopulation variation in the degree of sexual dimorphisms in fin lengths is mediated by sex steroid-dependent regulation of fin elongation.

Genetic Architecture of the Variation in Male-Specific Ossified Processes on the Anal Fins of Japanese Medaka

Traits involved in reproduction evolve rapidly and show great diversity among closely related species. However, the genetic mechanisms that underlie the diversification of courtship traits are mostly unknown. Japanese medaka fishes (Oryzias latipes) use anal fins to attract females and to grasp females during courtship; the males have longer anal fins with male-specific ossified papillary processes on the fin rays. However, anal fin morphology varies between populations: the southern populations tend to have longer anal fins and more processes than the northern populations. In the present study, we conducted quantitative trait locus (QTL) mapping to investigate the genetic architecture underlying the variation in the number of papillary processes of Japanese medaka fish and compared the QTL with previously identified QTL controlling anal fin length. First, we found that only a few QTL were shared between anal fin length and papillary process number. Second, we found that the numbers of papillary processes on different fin rays often were controlled by different QTL. Finally, we produced another independent cross and found that some QTL were repeatable between the two crosses, whereas others were specific to only one cross. These results suggest that variation in the number of papillary processes is polygenic and controlled by QTL that are distinct from those controlling anal fin length. Thus, different courtship traits in Japanese medaka share a small number of QTL and have the potential for independent evolution.

Neofunctionalization of Androgen Receptor by Gain-of-Function Mutations in Teleost Fish Lineage

Steroid hormone receptor family provides an example of evolution of diverse transcription factors through whole-genome duplication (WGD). However, little is known about how their functions have been evolved after the duplication. Teleosts present a good model to investigate an accurate evolutionary history of protein function after WGD, because a teleost-specific WGD (TSGD) resulted in a variety of duplicated genes in modern fishes. This study focused on the evolution of androgen receptor (AR) gene, as two distinct paralogs, ARα and ARβ, have evolved in teleost lineage after TSGD. ARα showed a unique intracellular localization with a higher transactivation response than that of ARβ. Using site-directed mutagenesis and computational prediction of protein-ligand interactions, we identified two key substitutions generating a new functionality of euteleost ARα. The substitution in the hinge region contributes to the unique intracellular localization of ARα. The substitution on helices 10/11 in the ligand-binding domain possibly modulates hydrogen bonds that stabilize the receptor-ligand complex leading to the higher transactivation response of ARα. These substitutions were conserved in Acanthomorpha (spiny-rayed fish) ARαs, but not in an earlier branching lineage among teleosts, Japanese eel. Insertion of these substitutions into ARs from Japanese eel recapitulates the evolutionary novelty of euteleost ARα. These findings together indicate that the substitutions generating a new functionality of teleost ARα were fixed in teleost genome after the divergence of the Elopomorpha lineage. Our findings provide a molecular explanation for an adaptation process leading to generation of the hyperactive AR subtype after TSGD.

Natural allelic variations of xenobiotic-metabolizing enzymes affect sexual dimorphism in Oryzias latipes

Sexual dimorphisms, which are phenotypic differences between males and females, are driven by sexual selection. Interestingly, sexually selected traits show geographical variations within species despite strong directional selective pressures. This paradox has eluded many evolutionary biologists for some time, and several models have been proposed (e.g. 'indicator model' and 'trade-off model'). However, disentangling which of these theories explains empirical patterns remains difficult, because genetic polymorphisms that cause variation in sexual differences are still unknown. In this study, we show that polymorphisms in cytochrome P450 (CYP) 1B1, which encodes a xenobiotic-metabolizing enzyme, are associated with geographical differences in sexual dimorphism in the anal fin morphology of medaka fish (Oryzias latipes). Biochemical assays and genetic cross experiments show that high- and low-activity CYP1B1 alleles enhanced and declined sex differences in anal fin shapes, respectively. Behavioural and phylogenetic analyses suggest maintenance of the high-activity allele by sexual selection, whereas the low-activity allele possibly has experienced positive selection due to by-product effects of CYP1B1 in inferred ancestral populations. The present data can elucidate evolutionary mechanisms behind genetic variations in sexual dimorphism and indicate trade-off interactions between two distinct mechanisms acting on the two alleles with pleiotropic effects of xenobiotic-metabolizing enzymes.

Evaluation of the 5α-reductase inhibitor finasteride on reproduction and gonadal development in medaka, Oryzias latipes

5-α reductase (5αR) inhibitors have an anti-androgenic effect in mammals because they inhibit the conversion of testosterone to the potent androgen, dihydrotestosterone. Finasteride is a type-2 5αR inhibitor that is used as a human pharmaceutical for the treatment of prostate cancer, benign prostate hyperplasia and male pattern baldness. This study evaluated the impacts of finasteride (50, 500 and 5000μg/L) on the development and reproduction of medaka (Oryzias latipes) exposed continuously over multiple generations (F0, F1 and F2). The exposure was initiated with reproductively mature fish (F0 generation) and continued until the hatching of the F2 generation. There were no significant effects on survival, fecundity or fertility in the F0 (50, 500, 5000μg/L) and F1 (50, 500μg/L) generations. The F1 generation exposed to 5000μg/L exhibited significant mortality. Histopathology of the gonads demonstrated that medaka and pre-clinical species respond similarly to finasteride exposure. Intersex condition and maldeveloped gonads were observed in F0 generation males exposed to 5000μg/L and F1 generation males exposed to 500μg/L. F1 generation males exposed to 500μg/L displayed reduced gonadosomatic index with an increased incidence of testicular degeneration. Males in both generations exhibited an increased incidence of Leydig cell hyperplasia at concentrations ⩾500μg/L. F0 generation females exposed to 5000μg/L exhibited increased gonadosomatic index. An increased prevalence of accelerated post-ovulatory follicle involution was observed in females at concentrations ⩾500μg/L in both generations. The gonadal changes induced by finasteride support the idea that 5-α reductase inhibition impacts androgen signaling in fish. Results from this study are discussed in the context of differential expression of the androgen receptor between species of fish.

Differential ligand selectivity of androgen receptors α and β from Murray-Darling rainbowfish (Melanotaenia fluviatilis)

Androgen receptors (ARs) mediate the physiological effects of androgens in vertebrates. In fishes, AR-mediated pathways can be modulated by aquatic contaminants, resulting in the masculinisation of female fish or diminished secondary sex characteristics in males. The Murray-Darling rainbowfish (Melanotaenia fluviatilis) is a small-bodied freshwater teleost used in Australia as a test species for environmental toxicology research. We determined concentration-response profiles for selected agonists and antagonists of rainbowfish ARα and ARβ using transient transactivation assays. For both ARα and ARβ, the order of potency of natural agonists was 11-ketotestosterone (11-KT)>5α-dihydrotestosterone>testosterone>androstenedione. Methyltestosterone was a highly potent agonist of both receptors relative to 11-KT. The relative potency of the veterinary growth-promoting androgen, 17β-trenbolone, varied by more than a factor of 5 between ARα and ARβ. The non-steroidal anti-androgen bicalutamide exhibited high inhibitory potency relative to the structurally related model anti-androgen, flutamide. The inhibitory potency of the agricultural fungicide, vinclozolin, was approximately 1.7-fold relative to flutamide for ARα, but over 20-fold in the case of ARβ. Fluorescent protein tagging of ARs showed that the rainbowfish ARα subtype is constitutively localised to the nucleus, while ARβ is cytoplasmic in the absence of ligand, an observation which agrees with the reported subcellular localisation of AR subtypes from other teleost species. Collectively, these data suggest that M. fluviatilis ARα and ARβ respond differently to environmental AR modulators and that in vivo sensitivity to contaminants may depend on the tissue distribution of the AR subtypes at the time of exposure.

Ontogenetic stage-specific quantitative trait loci contribute to divergence in developmental trajectories of sexually dimorphic fins between medaka populations

Sexual dimorphism can evolve when males and females differ in phenotypic optima. Genetic constraints can, however, limit the evolution of sexual dimorphism. One possible constraint is derived from alleles expressed in both sexes. Because males and females share most of their genome, shared alleles with different fitness effects between sexes are faced with intralocus sexual conflict. Another potential constraint is derived from genetic correlations between developmental stages. Sexually dimorphic traits are often favoured at adult stages, but selected against as juvenile, so developmental decoupling of traits between ontogenetic stages may be necessary for the evolution of sexual dimorphism in adults. Resolving intralocus conflicts between sexes and ages is therefore a key to the evolution of age-specific expression of sexual dimorphism. We investigated the genetic architecture of divergence in the ontogeny of sexual dimorphism between two populations of the Japanese medaka (Oryzias latipes) that differ in the magnitude of dimorphism in anal and dorsal fin length. Quantitative trait loci (QTL) mapping revealed that few QTL had consistent effects throughout ontogenetic stages and the majority of QTL change the sizes and directions of effects on fin growth rates during ontogeny. We also found that most QTL were sex-specific, suggesting that intralocus sexual conflict is almost resolved. Our results indicate that sex- and age-specific QTL enable the populations to achieve optimal developmental trajectories of sexually dimorphic traits in response to complex natural and sexual selection.