Orton F , Säfholm M , Jansson E , Carlsson Y , Eriksson A , Fick J , Uren Webster T , McMillan T , Leishman M , Verbruggen B , Economou T , Tyler CR , Berg C .

	Figure 1. Experimental design. Fertilised eggs were exposed to flutamide or linuron until life stage 66 when they were transferred to clean water. Gene expression analysis on a set of target genes was conducted to identify their utility for distinguishing gonadal sex and to investigate the effects of anti-androgens on their expression in developing tadpoles (during sensitive windows of sexual development). Effects of anti-androgens on gonadal development and sex ratio were assessed in metamorphs via histopathology. Frogs were then maintained to adulthood in clean water when they were assessed for reproductive competence and morphology. Breeding studies were conducted to assess the effects of early life exposure to anti-androgens on the gonadal sexual development, secondary sex characters in males (nuptial pads and forearm) and breeding success (fecundity, fertility and time in amplexus). a=sex determination, b=sex differentiation, c=completed sex differentiation. SSC=secondary sexual characteristic, exp.=experimental.
	Figure 2. Sex ratio (a) in stage 51–53 tadpoles (control/linuron low n = 24, linuron high n = 23, flutamide n = 16), stage 55–58 tadpoles (control/linuron high n = 24, linuron low/flutamide n = 16), metamorphs (control n = 37, linuron low n = 24, linuron high n = 33, flutamide n = 11) and adults (control n = 82, linuron low n = 84, linuron high n = 85, flutamide n = 47). Solid bars = males, dotted bars = females. Gonadal differentiation in female (b) control n = 12, linuron low n = 10, linuron high n = 9, flutamide n = 4) and male (c) control n = 23, linuron low n = 10, linuron high n = 22, flutamide n = 5) metamorphs, and female (D: control/linuron low/linuron high n = 9, flutamide n = 7) and male (E: control/linuron low/linuron high/ flutamide: n = 10) adults. *indicates significant difference compared with control group.
	Figure 3. Relationship between cyp19 and amh for distinguishing gonadal sex in tadpole gonado-mesonephros complex (GMC: a,b), gonadal ontogeny of gene expression (dmrt1, cyp17, amh, ar, foxl2, cyp19) in control tadpoles (c,d) and the interrelationships between these genes in control tadpole GMC (e,f). Unsupervised cluster analysis was used to delineate ‘males’ (black) and ‘females’ (red) during sex determination (a) and sex differentiation (b) for all treatment groups collectively (data not shown for dmrt1 v. foxl2/cyp19 or amh v. foxl2 as outcomes of the analysis did not differ). Each letter represents one individual and the asterices indicate the centre of each cluster (C=control, L=linuron low, H=linuron high, F=flutamide). Ontogeny of gene expression in males (normal, filled symbols) and females (italic, open symbols) tadpole GMC (c,d), bold indicates significant change during ontogeny using linear regression analysis *expression is relative to gonad pool standards (see methods for details). Interrelationships between gene expression in male (e) and female (f) GMC analysed using Pearsons’ (or Spearmans’ for non-parametric data), *p < 0.05, **p < 0.01, ***p < 0.001.
	Figure 4. Effects of anti-androgens on expression of target genes during gonad ontogeny (NF stage 51–58). Linear regression of dmrt1 (a), cyp17 (b), amh (c), ar (d), foxl2 (e) and cyp19 (f). Black lines = control males (dmrt1 r2 = 0.59, cyp17 r2 = 0.62, amh r2 = 0.64, foxl2/cyp19 r2 = ns) green lines = linuron low males (dmrt1 r2 = 0.57, cyp17 r2 = 0.67, amh r2 = 0.53, foxl2/cyp19 r2 = ns), blue lines = linuron high males (dmrt1 r2 = 0.79, cyp17 r2 = 0.51, amh r2 = 0.65, foxl2 r2 = 0.58, cyp19 r2 = ns), red lines = flutamide males (dmrt1 r2 = 0.81, cyp17 r2 = 0.83, amh r2 = 0.9, foxl2 r2 = 0.33, cyp19 r2 = ns), grey lines = control females (dmrt1 r2 = 0.41, cyp17 r2 = ns, amh r2 = 0.33, foxl2 r2 = 0.52, cyp19 r2 = ns), light green = linuron low females (dmrt1 r2 = 0.37, cyp17 r2 = ns, amh r2 = 0.79, foxl2 r2 = 0.81, cyp19 r2 = ns), light blue = linuron high females (dmrt1 r2 = 0.42, cyp17 r2 = ns, amh r2 = ns, foxl2 r2 = 0.31, cyp19 r2 = 0.35), pink = flutamide females (dmrt1 r2 = ns, cyp17 r2 = ns, amh r2 = 0.32, foxl2 r2 = 0.50, cyp19 r2 = ns); with matching filled (male) or open (female) symbols. For MANCOVA analysis during sex determination (stage 51–53) or sex differentiation (stage 55–58) effects of flutamide were observed on foxl2 levels in ‘females’ (see Supplementary Table S5 for numerical values/more details). Dashed lines represent significant differences in the linear regression slope compared with controls (ANCOVA, stage as a covariate). *relative gene expression compared to gonad pool standards.
	Figure 5. Effects of antiandrogens on expression of target genes in the brain in both male (circles) and female (diamonds) metamorphs (NF stage 66). Dmrt1 (a), cyp17 (b), amh (c), ar (d), foxl2 (e) and cyp19 (f) expression in controls (filled symbols), linuron low (open symbols), linuron high (dotted symbols) and flutamide (crossed symbols). Control male, n = 22; control female, n = 10; linuron low male, n = 9; linuron low female, n = 7; linuron high male, n = 22; linuron high female, n = 10; flutamide male, n = 5; flutamide female, n = 4. *relative gene expression compared to brain pool standards (see methods for details).

Al-Asaad, Müllerian inhibiting substance in the caudate amphibian Pleurodeles waltl. 2013, Pubmed

Al-Asaad, Müllerian inhibiting substance in the caudate amphibian Pleurodeles waltl. 2013, Pubmed
Aoyama, Expression of Dmrt1 protein in developing and in sex-reversed gonads of amphibians. 2003, Pubmed
Baron, Expression profiling of candidate genes during ovary-to-testis trans-differentiation in rainbow trout masculinized by androgens. 2008, Pubmed
Berg, Xenopus tropicalis as a test system for developmental and reproductive toxicity. 2009, Pubmed , Xenbase
Bhatia, Effects of the commercial antiandrogen flutamide on the biomarkers of reproduction in male Murray rainbowfish (Melanotaenia fluviatilis). 2014, Pubmed
Chakrabarty, Endosulfan and flutamide, alone and in combination, target ovarian growth in juvenile catfish, Clarias batrachus. 2012, Pubmed
Duarte-Guterman, Transcript profiles and triiodothyronine regulation of sex steroid- and thyroid hormone-related genes in the gonad-mesonephros complex of Silurana tropicalis. 2011, Pubmed , Xenbase
Egea-Serrano, Understanding of the impact of chemicals on amphibians: a meta-analytic review. 2012, Pubmed
El Jamil, Early aspects of gonadal sex differentiation in Xenopus tropicalis with reference to an antero-posterior gradient. 2008, Pubmed , Xenbase
Fujii, Involvement of androgen receptor in sex determination in an amphibian species. 2014, Pubmed
Gyllenhammar, Clotrimazole exposure modulates aromatase activity in gonads and brain during gonadal differentiation in Xenopus tropicalis frogs. 2009, Pubmed , Xenbase
Gyllenhammar, Reproductive toxicity in Xenopus tropicalis after developmental exposure to environmental concentrations of ethynylestradiol. 2009, Pubmed , Xenbase
Hayes, Atrazine induces complete feminization and chemical castration in male African clawed frogs (Xenopus laevis). 2010, Pubmed , Xenbase
Hayes, Characterization of atrazine-induced gonadal malformations in African clawed frogs (Xenopus laevis) and comparisons with effects of an androgen antagonist (cyproterone acetate) and exogenous estrogen (17beta-estradiol): Support for the demasculinization/feminization hypothesis. 2006, Pubmed , Xenbase
Iwade, Elevated expression of P450c17 (CYP17) during testicular formation in the frog. 2008, Pubmed
Jansson, Sex-dependent expression of anti-Müllerian hormone (amh) and amh receptor 2 during sex organ differentiation and characterization of the Müllerian duct development in Xenopus tropicalis. 2016, Pubmed , Xenbase
Jobling, Statistical modeling suggests that antiandrogens in effluents from wastewater treatment works contribute to widespread sexual disruption in fish living in English rivers. 2009, Pubmed
Kalt, Morphology and kinetics of spermatogenesis in Xenopus laevis. 1976, Pubmed , Xenbase
Knechtges, A novel amphibian tier 2 testing protocol: a 30-week exposure of Xenopus tropicalis to the antiandrogen flutamide. 2007, Pubmed , Xenbase
Kvarnryd, Early life progestin exposure causes arrested oocyte development, oviductal agenesis and sterility in adult Xenopus tropicalis frogs. 2011, Pubmed , Xenbase
Lange, Environmental chemicals active as human antiandrogens do not activate a stickleback androgen receptor but enhance a feminising effect of oestrogen in roach. 2015, Pubmed
Lee, Subchronic exposure to low concentrations of di-n-butyl phthalate disrupts spermatogenesis in Xenopus laevis frogs. 2005, Pubmed , Xenbase
Levy, Bisphenol A induces feminization in Xenopus laevis tadpoles. 2004, Pubmed , Xenbase
Main, Genital anomalies in boys and the environment. 2010, Pubmed
Mann, Amphibians and agricultural chemicals: review of the risks in a complex environment. 2009, Pubmed
Monosson, Peripubertal exposure to the antiandrogenic fungicide, vinclozolin, delays puberty, inhibits the development of androgen-dependent tissues, and alters androgen receptor function in the male rat. 1999, Pubmed
Muller, Effect of castration and testosterone treatment on fiber width of the flexor carpi radialis muscle in the male frog (Rana temporaria L.). 1969, Pubmed
Navarro-Martín, Sexing frogs by real-time PCR: using aromatase (cyp19) as an early ovarian differentiation marker. 2012, Pubmed
Ohtani, Role of aromatase and androgen receptor expression in gonadal sex differentiation of ZW/ZZ-type frogs, Rana rugosa. 2003, Pubmed
Olmstead, Sex reversal of the amphibian, Xenopus tropicalis, following larval exposure to an aromatase inhibitor. 2009, Pubmed , Xenbase
Orton, Do hormone-modulating chemicals impact on reproduction and development of wild amphibians? 2015, Pubmed , Xenbase
Orton, Body size, nuptial pad size and hormone levels: potential non-destructive biomarkers of reproductive health in wild toads (Bufo bufo). 2014, Pubmed
Orton, Agricultural intensity in ovo affects growth, metamorphic development and sexual differentiation in the common toad (Bufo bufo). 2011, Pubmed
Orton, Endocrine disrupting effects of herbicides and pentachlorophenol: in vitro and in vivo evidence. 2009, Pubmed , Xenbase
Orton, Widely used pesticides with previously unknown endocrine activity revealed as in vitro antiandrogens. 2011, Pubmed
Oshima, Molecular cloning and gene expression of Foxl2 in the frog Rana rugosa. 2008, Pubmed , Xenbase
Pettersson, Persistent sex-reversal and oviducal agenesis in adult Xenopus (Silurana) tropicalis frogs following larval exposure to the environmental pollutant ethynylestradiol. 2006, Pubmed , Xenbase
Pettersson, Environmentally relevant concentrations of ethynylestradiol cause female-biased sex ratios in Xenopus tropicalis and Rana temporaria. 2007, Pubmed , Xenbase
Piprek, Sexual dimorphism of AMH, DMRT1 and RSPO1 localization in the developing gonads of six anuran species. 2013, Pubmed , Xenbase
Rhen, Molecular mechanisms of sex determination in reptiles. 2010, Pubmed
Stuart, Status and trends of amphibian declines and extinctions worldwide. 2004, Pubmed
Säfholm, Molecular and histological endpoints for developmental reproductive toxicity in Xenopus tropicalis: Levonorgestrel perturbs anti-Müllerian hormone and progesterone receptor expression. 2016, Pubmed , Xenbase
Säfholm, Risks of hormonally active pharmaceuticals to amphibians: a growing concern regarding progestagens. 2014, Pubmed , Xenbase
Säfholm, Disrupted oogenesis in the frog Xenopus tropicalis after exposure to environmental progestin concentrations. 2012, Pubmed , Xenbase
Vandenberg, Hormones and endocrine-disrupting chemicals: low-dose effects and nonmonotonic dose responses. 2012, Pubmed
Villalpando, Determination of the sensitive stages for gonadal sex-reversal in Xenopus laevis tadpoles. 1990, Pubmed , Xenbase
Wake, Ecology. Facing extinction in real time. 2012, Pubmed
Wilson, The herbicide linuron reduces testosterone production from the fetal rat testis during both in utero and in vitro exposures. 2009, Pubmed
Woudneh, Pesticide multiresidues in waters of the Lower Fraser Valley, British Columbia, Canada. Part I. Surface water. 2009, Pubmed
Yoshimoto, A W-linked DM-domain gene, DM-W, participates in primary ovary development in Xenopus laevis. 2008, Pubmed , Xenbase
van OORDT, The effect of testosterone on the spermatogenesis of the common frog, Rana temporaria. 1960, Pubmed
van Wyk, The effects of anti-androgenic and estrogenic disrupting contaminants on breeding gland (nuptial pad) morphology, plasma testosterone levels, and plasma vitellogenin levels in male Xenopus laevis (African clawed frog). 2003, Pubmed , Xenbase