Systems toxicology of male reproductive development: profiling 774 chemicals for molecular targets and adverse outcomes.
Exposure to chemicals during prenatal development may increase the risk of adverse outcomes, and biomonitoring studies suggest pregnant women are exposed to multiple environmental chemicals (Woodruff et al. 2011). Adverse trends in male developmental reproductive health have been reported for rates of testicular germ cell tumors (TGCT), low semen quality, and relatively common human developmental defects such as undescended testes (cryptorchidism) and malformations of the genital tubercle (e.g., hypospadias) (Sharpe and Skakkebaek 2008; Virtanen et al. 2005).
The testicular dysgenesis syndrome (TDS) hypothesis posits an interrelationship among these adverse outcomes, as manifestations of altered prenatal testicular development in humans (Aschim et al. 2004; Bay et al. 2006; Biggs et al. 2002; Kaleva and Toppari 2005; Scott et al. 2007; Sharpe and Skakkebaek 2008; Skakkebaek et al. 2001). Epidemiological studies, however, provide scant support for a shared mechanistic origin of the four elements contributing to the TDS hypothesis (decreased sperm counts/infertility/subfertility, and common developmental defects such as cryptorchidism and hypospadias, and increasing incidences of TGCT) (Akre and Richiardi 2009). On the other hand, human studies have reported associations between at least some adverse outcomes of male reproductive tract development: shortened anogenital distance (AGD) in boys with undescended testis (Jain and Singal 2013); reduced AGD/ penile length in patients with hypospadias or cryptorchidism (Thankamony et al. 2014); and TGCT with cryptorchidism/hypospadias/ genital malformations (Trabert et al. 2013). The TDS hypothesis is difficult to test experimentally because of the inaccessibility of the fetal testis during formation and organization (8-12 weeks gestation), the lengthy period between induction and manifestation of some adult outcomes (20-45 years), the lack of definitive etiology for prenatal studies underlying hypospadias/cryptorchidism, and the lack of animal models for TGCT (Sharpe and Skakkebaek 2008).
Evidence for reduced androgenicity in the human TDS hypothesis comes from linkage studies correlating TDS elements to familial mutations in the androgen receptor (AR) (Lottrup et al. 2013). Studies using animal models of prenatal endocrine dysfunction have shown increased incidence rates for cryptorchidism, hypospadias, and low sperm quality following exposure to environmental compounds that may be acting through an anti-androgenic mode of action (Gray et al. 2006; Hu et al. 2014; Mylchreest et al. 1998, 2000, 2002; Sharpe and Skakkebaek 2008; Virtanen et al. 2005; Wilson et al. 2008). For example, exposure of pregnant rats to dibutyl phthalate (DBP) during the critical period of male reproductive development resulted in reduced AGD, increased cryptorchidism and hypospadias, testicular atrophy with germ cell loss, and weight reductions in the epididymis, seminal vesicles, and prostate (Mylchreest et al. 1998, 2000). Exposure of pregnant rats to di(2-ethylhexyl) phthalate (DEHP) resulted in abnormal testes development with large aggregates or clusters of Leydig cells in the interstitial spaces, multinucleated germ cells in the seminiferous cords, and significant reductions in testosterone levels in the fetal rat testis (Parks et al. 2000). Phthalates may increase testicular testosterone production (low-dose effect) in fetal rats, but higher dosages decrease testosterone production in the fetal testis, thereby leading to reduced AGD and increased rates of cryptorchidism (Ge et al. 2007). The effect on fetal testosterone production has been linked to reduced expression of steriodogenic enzymes (CYP17, CYP11a, StAR) and insulin-like 3 (INSL3) (Wilson et al. 2008). INSL3, a Leydig cell hormone that drives development of the gubernaculum, may underlie phthalate-induced cryptorchidism (Johnson et al. 2010). These findings support the connectivity to prenatal disruption of androgen signaling (AR mutations in the human TDS hypothesis) or androgen homeostasis (phthalate syndrome in rats) during the critical period of sexual differentiation.
Further evaluation of the fetal rat testis proteome following prenatal exposure to DEHP revealed altered cellular differentiation and cell migration pathways, as well as decreased testosterone production consistent with an anti-androgenic effect (Veeramachaneni and Klinefelter 2014); however, those authors also found that plasma estradiol was elevated approximately 2-fold in rat fetuses exposed to DEHP (Veeramachaneni and Klinefelter 2014). Alterations in estrogen pathways may also contribute to the prenatal origins of at least some elements of TDS (humans) or phthalate syndrome (rats). Although the synthetic estrogen diethylstilbestrol (DES) can impair testosterone production in fetal rats, the estrogen receptor (ESR1) is not expressed in human fetal Leydig cells (Mitchell et al. 2013). As such, our understanding of male reproductive developmental toxicity following estrogenic and/or anti-androgenic exposures is incomplete (Ono and Harley 2013).
Abnormal testis development could have numerous other primary causes beyond prenatal disruption of AR or ESR1 pathways. A human genome-wide association study (GWAS) on a cohort of 488 patients with symptoms of TDS, comprising 212 patients with TGCT, 138 with cryptorchidism, 31 with hypospadias, and 107 with infertility found an association with genetic markers located in the region of TGF$R3 (transforming growth factor, beta receptor III) and BMP7 (bone morphogenetic protein 7) (Dalgaard et al. 2012). GWAS for developmental defects of the male reproductive system have also identified candidate genes for hormonal (AR, steroidogenesis, ESR1, and ESR2) and nonhormonal (FGF8, FGFR2, BMP4, and BMP7) pathways underlying hypospadias (Carmichael et al. 2013; Geller et al. 2014). The association of hypospadias to fibroblast growth factor (FGF) signaling pathways in humans is consistent with mouse studies in which various degrees of hypospadias have been observed in mutations that affect FGF signaling, either alone or in interaction with pathways directing the spatial patterning of male genital tubercle morphogenesis (e.g., SHH, BMP, WNT signals) (Kojima et al. 2010).
Male reproductive developmental disorders that involve elements in the hypothesized TDS (humans) or a phthalate syndrome (rats) thus reflect a complex disease process in which fundamental and varied genetic and environmental factors interact to perturb morphogenetic signals that spatially pattern the male reproductive tract, or otherwise disrupt endocrine signals that further regulate its growth and maturation during sexual diversification. Given this complexity, a systems analysis may help unravel the molecular and functional changes leading to male reproductive developmental toxicity (MRDT). The recent availability of high-throughput screening (HTS) data, quantitative structure-activity relationship (QSAR) models, and systems modeling tools enables mining of toxicity information across a vast chemical landscape (Collins et al. 2008; Kavlock et al. 2012; NRC 2007; Sturla et al. 2014).
In the present study we used an HTS-driven approach to unravel the potential connectivity of environmental exposures to diverse molecular pathways and cellular processes that are disrupted in male reproductive developmental toxicity. The aims of this study were to a) provide a comprehensive evaluation of structurally diverse chemicals with male reproductive developmental effects in laboratory animal studies; b) identify chemical-bioactivity and chemical-end point relationships that are anchored to altered developmental phenotypes of the male reproductive system; and c) provide a cohesive network of interactions underlying the molecular and functional changes predicted to invoke male reproductive developmental toxicity.
All in vivo data for the current work was extracted from ToxRefDB (Toxicity Reference Database, http://actor.epa.gov/; Figure 1A,B). This collection of publicly available in vivo data built from thousands of guideline animal studies includes guideline chronic/ cancer rat and mouse studies (Martin et al. 2009a), multigenerational reproduction rat studies (Martin et al. 2009b), and prenatal developmental studies (mostly pregnant rats and rabbits) (Knudsen et al. 2009). The ToxRefDB version used for the present study (release date 17 December 2012) had 6,021 animal studies on 1,113 chemicals, including 3,684 studies on pesticides and industrial chemicals, including perfluorinated compounds and phthalates; 151 studies with mammalian preclinical data for failed pharmaceutical agents; 434 studies from the National Toxicology Program (NTP); and 392 studies from open literature (PubMed). Figure 1A summarizes the ToxRefDB chemical coverage by study type and test species. Chemistry, biology, and usage features of these chemicals are accessible in the Aggregated Computational Toxicology Resource (http://actor.epa.gov/).
Because many ToxRefDB studies were not designed to assess male reproductive developmental end points, we searched the biomedical literature in PubMed for additional evidence of developmental toxicity for chemicals in ToxRefDB and ToxCastDB (Figure 2A,B; see also Excel File S1). Raw data were acquired from MEDLINE/PubMed Baseline Repository using MySQL on 15 May 2015 (http://mbr.nlm.nih.gov/) and processed in Microsoft Excel as Pivot Table. The chemicals were matched to the MeSH vocabulary to find the MeSH name and ID, either by name, by synonym, or CAS (Chemical Abstracts Service) number. Annotations of male reproductive end points were retrieved by searching on the MeSH categories "Urogenital Abnormalities," "Spermatogenesis," "Spermatozoa," "Semen Analysis," "Genitalia," "Male," and "Genital Neoplasms, Male." Annotations of male reproductive developmental effects were retrieved by searching on the MESH terms "Cryptorchidism," "Hypospadias," and "Urogenital Abnormalities." Because an end point could be co-annotated with a chemical when a lack of effect was discussed (i.e., studies with null findings or studies that did not assess any given effect), the annotation was manually curated to make sure whether or not the effect was found. To be inclusive, a chemical that showed an effect in one study would be counted as "positive" even if other comparable studies did not report that effect.
Strategy to Define Male Reproductive Developmental Toxicity
The present study defined MRDT as all male reproductive birth defects reported in ToxRefDB and/or PubMed, as well as the sperm effects, histological effects, and tissue weight change reported exclusively in developmental studies in ToxRefDB. In other words, if an effect was reported in a developmental study and also reported in a nondevelopmental study, it would not be counted as MRDT. A total of 1,113 chemicals in ToxRefDB were filtered for in vivo phenotypic analysis based on the workflow described in Figure 1B. For chemicals that had a relevant study in ToxRefDB, the workflow selected those positive for any adverse male reproductive outcome. A chemical was considered positive for an effect when one or more study type in ToxRefDB recorded a relevant end point, regardless of dosage. That list was further parsed into a subset of chemicals where male reproductive developmental toxicity was evident in ToxRefDB for an early life stage exposure scenario. These data were culled from ToxRefDB studies of any type, but specifically for end points that could be linked to a prenatal developmental toxicity study or generational reproductive study. This strategy pointed to developmental defects including malformations, relative weight change, and histological effects for male reproductive organs (penis, testis, prostate, seminal vesicle), and any effects on sperm parameters or testicular cancers that were not strictly invoked in an adult-only chronic/cancer study type.
Statistical Analysis for Male Reproductive End points
Comparisons for ToxRefDB studies across species (rat vs. mouse, rat vs. rabbit) and across generations (P, parental vs. F1, filial) used Pearson's correlation in R-version 3.1.0 (R Core Team 2013). Co-occurrence of male reproductive end points was visualized in Microsoft Excel 2013 by grouping chemicals on a logarithmic scale ([log.sub.10]) by the least common end point. Chemicals passing through the workflow filter (Figure 1B) were subjected to two-directional hierarchical clustering in R version 3.1.0 using Euclidean distance for measure and Ward's method for linkage analysis (Ward 1963). When the end point reflected a change to some parameter of male reproductive development in ToxRefDB, the lowest effect level (LEL) was determined for each category (cLELs for malformations, relative weight change, histological effects, and sperm parameters). If, for a particular chemical, the effects data were missing for a single category, the highest cLEL value for that chemical was multiplied by 10 to impute a value for clustering purposes. This does not imply that an actual effect is predicted, but rather enables the clustering across features in a hypothetical or virtual syndrome, where the missing value might be observed had a high enough dose been tested. Next, an "effect score" was calculated for each chemical and end point category to compute relative potency based on the cLEL values. The effect score was computed from the -[log.sub.10] (cLEL mg/kg/day) dosage and normalizing the data set from 0 to 100 across the entire matrix. On this scale, 100 reflects the most potent chemical-end point category pair and 0 the least potent.
All in vitro data for the current work was extracted from 1,860 chemical-bioactivity effects data in ToxCastDB (release date 15 December 2014) containing assay data and assay information (http://actor.epa.gov/ dashboard/). The large battery of molecular and cellular HTS assays (Kavlock et al. 2012) has been used in earlier versions to develop predictive signatures of diverse chemical-bioactivity relationships in prenatal developmental toxicity (Sipes et al. 2011). The ToxCastDB version for the present study included 328 HTS assays pointing to 795 cellular features and 293 annotated molecular targets (BioAssay Ontology version 1.6). Some molecular targets among the 293 were assayed in multiple orthogonal platforms, whereas others may have been unique to one technology platform. The data processing pipeline (http://epa.gov/ncct/toxcast/data. html) included raw data quality assurance checks, normalization, automated curve fitting, visual plots of the concentration-response relationships, and calculation of concentration causing half maximal response in activity ([AC.sub.50]) or the lowest effect concentration (LEC).
Gene Score Calculation
We used a panel of cytotoxicity assays in ToxCastDB to account for cytotoxicity-related assay interference. This accounted for nonspecific bioactivity that was presumably caused by integrated stress responses (Harding et al. 2003). The assay interference was assessed with a collection of 37 cellular features that indicated cytotoxicity, other forms of cell loss, or decrease in proliferation in several cell lines and primary cell types. For any given pairing of a chemical and a cytotoxicity end point, we defined a "hit" as a situation in which a) the best fit to the concentration-response data was given by either a Hill model or gain-loss model, and b) the maximum response predicted by the best-fitting model was greater than some critical threshold.
For a chemical X with two or more cytotoxicity hits, we analyzed the set of logarithms ([log.sub.10]) of the [AC.sub.50] values corresponding to those hits. We let m(X) denote the median of this set and computed the median absolute deviation (MAD) of the set. Then, we computed [M.sub.global] as the median of the MADs for all chemicals. Finally, for each chemical X and each cytotoxicity end point A, we computed a standardized score:
Z (X,A) = - [log.sub.10][[AC.sub.50] (X,A)] - m (X)/[M.sub.global]. 
For a chemical X with fewer than two cytotoxicity hits, we set m(X) = 3 for all A.
Because we considered multiple cellular features and many molecular targets, we computed "gene score" to summarize the ToxCastDB data. Gene score for chemical X and gene g was given by
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] 
and [n.sub.assay](g) was the number of cellular features corresponding to gene g, and [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] was the set that contained all of cellular features. A hit with a large value of G occurred at concentrations significantly below cytotoxicity and therefore indicative of selective bioactivity. This excluded most of the cytotoxicity-driven hits, but allowed borderline cases to be included. The modified gene score as a micromolar concentration can be computed as
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. 
Bipartite Network Visualization
The workflow for in vitro network analysis focused on an expanded list of ToxCast[TM] chemicals with male reproductive developmental effects identified from ToxRefDB (Figure 1B) and PubMed literature search to capture male developmental toxicants that might be missed or excluded from the chemical list (Figure 1B). Hierarchical clustering was performed on ToxCastDB effects using the same R metrics for the ToxRefDB analysis (Ward 1963; R Core Team 2013). A bipartite network for chemical-assay effects was then constructed using ForceAtlas layout in Gephi version 0.8.2 beta (Bastian et al. 2009). Functional annotation of the nodes (assay target) and edges (chemical target) in the network graph used the NIH/NIAID (National Institutes of Health/National Institute of Allergy and Infectious Diseases) Database for Annotation, Visualization and Integrated Discovery (DAVID) version 6.7 (Huang et al. 2009a, 2009b), selecting terms from significant (p < 0.05) functional annotation clusters generally based on pathways and biological processes.
ToxRefDB Workflow for Male Reproductive Developmental Toxicity
A total of 774 ToxRefDB chemicals were identified in five study types that employed three test species (Figure 1A,B). ToxRefDB study types included prenatal developmental studies (maternal exposure during pregnancy; fetuses evaluated at term), one-generational studies (parental exposure, preconception through lactation; F1 filial generation evaluated postnatally); multigenerational reproductive studies (exposure continued to the F2 filial generation), subchronic studies (90-day adult exposure); and chronic studies (adult exposure up to 2 years). The relevant 3,479 studies came from testing on rats (2,274 studies), mice (672 studies), and rabbits (533 studies).
Twenty-three male reproductive end points, regardless of dosage, were reported for a total of 281 chemicals (36% of the 774) across 1,653 ToxRefDB studies. These 23 end points were classified into five categories for the male reproductive system, whether developmental or not: a) malformations, b) testicular tumors, c) sperm effects, d) histological effects, and e) relative weight change in male reproductive organs (Table 1; see also Excel File S2). Malformations encompassed reduced AGD, hypospadias, cryptorchidism, and abnormal nipple retention, all of which were linked to early life stage (prenatal, early postnatal) exposure scenarios. The testicular tumor types included adenoma, interstitial cell tumor, mesothelioma, and other (i.e., uncategorized in ToxRefDB). Although subchronic and chronic study types identified 31 chemicals that produced testicular tumors, seminomas (TGCT) were not reported in ToxRefDB studies. Therefore, the TGCT component of TDS could not be considered further. Sperm effects included alterations in morphology, motility, epididymal sperm count, testicular spermatid count, as well as semen sperm count. Histological effects were reported for the prostate, seminal vesicles, epididymis, and testes based on general features (e.g., inflammation, necrosis, and pigmentation) and spermatogenesis-specific abnormalities (e.g., aspermia, immature sperm, and retained spermatids). Tumor effects, sperm effects, and histological effects could result from early lifestage as well as adult exposure scenarios in ToxRefDB study designs; therefore, those adverse outcomes may indicate an adult manifestation of developmental toxicity and/or toxicity to the adult reproductive tract.
Table 1 categorizes 248 (88%) of the 281 chemicals that are both positive for male reproductive effects and have dosage information for prenatal and generational (P1, F1) studies. All 23 male reproductive effects registered by the ToxRefDB chemicals were reported in the rat as a test species. Most effects were also reported in mouse; however, no malformation and only one of the four tumor types were reported in the mouse. Only 41 (18%) of the 224 male reproductive toxicants in rats are also male reproductive toxicants in mice. These 41 chemicals accounted for 64% of the 64 male reproductive toxicants in mice. Despite this discordance for the rat, the total number of chemicals associated with sperm effects, histological effects, and relative weight change correlated strongly between these rodent species (Pearson's R = 0.93 across 15 effects). The mean cLEL was generally lower for the rat (50 [+ or -] 1 mg/kg/day) versus the mouse (120 [+ or -] 1 mg/kg/day), indicating that the rat was the more sensitive species. Rabbit prenatal studies reported male reproductive effects for two chemicals only (spiroxamine and sulfluramid) and could not be directly compared. Testicular tumors and malformations were specific to evaluation of either the parental (P1) or filial (F1) generation, respectively, in rat studies (mouse not tested). The number of chemicals associated with sperm effects, histological effects, and relative weight change correlated strongly between the P1 and F1 generations (R = 0.89). Mean cLELs for these end points was also similar (53 [+ or -] 1 mg/kg/day and 54 [+ or -] 1 mg/kg/day, respectively; Table 1).
The similarities between the human TDS hypothesis and the rodent phthalate syndrome motivated an examination of the co-occurrence of male reproductive effects across the 281 chemical list (Figure 3). A broken hierarchy was evident, with relative weight change accounting for the broadest response (227 chemicals, or 81% of 281). This was followed in turn by histological effects (126 chemicals, or 45%), sperm effects (58 chemicals, or 21%), testicular tumors (31 chemicals, or 11%), and malformations (18 chemicals, or 6.4%). Only vinclozolin produced both malformations and testicular tumors; and only 6 chemicals produced both malformations and sperm effects (Figure 3). The hierarchy of male reproductive effects resembled what might be expected in the human TDS hypothesis (i.e., testicular atrophy > sperm effects > tumors > malformations) (Akre and Richiardi 2009), with the caveat that seminomas (TGCT) could not be evaluated.
Figure 4 illustrates the hierarchical data structure of 63 MRDT chemicals by effect score in a heat map annotated by chemical class and adverse outcome. This analysis excluded 218 out of the 281 chemicals where adult male reproductive toxicity was observed in a chronic/cancer or subchronic study. This resulted in a subset of 63 MRDT chemicals that included 8 of 13 phthalates from the original (774 chemical) data set. The arrangement by effect score was similar to the broken hierarchy for the 281 male reproductive toxicants (Figure 3). Butylbenzyl phthalate stood out among the 63 chemicals in general, and the 13 phthalates in particular, as the sole MRDT with a high effect score across all four phenotypic categories (center of heat map in Figure 4). An enriched "phthalate cluster" (right side of the heat map in Figure 4) included diisobutyl phthalate, dicyclohexyl phthalate, dihexyl phthalate, and monobenzyl phthalate driven by a high effect score for malformations (primarily) and relative organ weight change (secondarily). In contrast to butylbenzyl phthalate, DES was the sole MRDT with a low effect score across all four phenotypic categories (left side of heat map in Figure 4). It was distinct from other phenolic compounds with estrogenic activity (4-nonylphenol, 4-octylphenol, bisphenol A, and isoeugenol) (Laws et al. 2000), which showed a high effect score in at least one phenotypic category. The phenolic estrogens 4-nonylphenol and 4-octylphenol co-clustered strongly with SAR115740 (TRPV1 antagonist) and chlorpromazine hydrochloride (dopamine antagonist), respectively.
ToxCastDB Bioactivity Profile of Male Reproductive Toxicants
To visualize the ToxCastDB bioactivity profiles on the 281 male reproductive toxicants, we performed hierarchical clustering of the chemicals based on gene score (Figure 5). ToxCast HTS data were available for 212 (75%) of the 281 toxicants (Figure 2B). Strong bioactivity was detected across a broad set of assays, most of which were co-culture cell-cell signaling platforms (BioSeek, BSK), multiplex transcriptional reporter platforms (Attagene, ATG), and HepG2 cell stress assays at 24 hr or 72 hr (APR and Tox21). The heatmap shown in Figure 5 is a hierarchical cluster of the 212 chemicals by 213 molecular targets based on gene score. The following are notable features in the heat map:
1. NFE2L2 and NR1I2: nuclear receptors in xenobiotic metabolism and response to oxidative stress
2. CREB3, FOS, and JUN: FOS/Jun (AP-1 transcription factor) family members that crosstalk with CREB
3. ESR1 and ESR2: estrogen receptors
4. CYP3A1, CYP2J2, CYP2B6, CYP2C11, and CYP2A1: cytochrome-P450s
5. SERPINE1, ICAM1, F3, MMP3, CCL2, CXCL5, CXCL9, and VCAM1; and CD69, PLAU, PLAT, SELP, SELE, HLA-DRA, IL1 A, PLAUR, SAA1, TIMP2, CCL26, TGFB1, NR4A2, CSF1, and CXCL10: tissue remodeling during inflammation and angiogenesis
6. RXRA and NR1H4: retinoid X receptor and a binding partner; RORB, RORA, and RORC: retinoic acid receptor-related orphan receptors; RARB, RARG, and VDR: retinoic acid and vitamin D receptors; and CEBPB, POU2F1, and PPARA
7. HTR2C, CHRM2, ADRA2B, HTR6, HTR5A, ADRA2C, HTR7, DRD2, and DRD4; SLC6A4, ADRA1A, HTR1A, ADRA1B, HTR2A, and DRD1; OPRL1, CHRM5, CHRM1, ADRB1, ADRA2A, and HRH1; and CHRM4, CHRM3, KCNH2, SLC6A4, CYP2D6, SCN1A: four clusters of neuronal GPCRs (G-proteincoupled receptors: adrenergic, cholinergic, histamine, and serotonin receptors).
In general, the promiscuity of chemical-target hits was greater on the chemical side than the molecular target side. A subset of chemicals in the NFE2L2 and NR1I2 cluster could also be found in the RAR, RXR, and AP-1 clusters. Many neuronal GPCR chemical actives were found in the HTR2C cluster, followed by the SLC6A4, OPRL1, and CHRM4 clusters. Some clusters and molecular targets responded to specific chemicals, such as the SERPINE1 cluster (tributyltin chloride and triphenyltin hydroxide), FOX01 and FOX03 (simazine), PTPN1 and PTPN4 (sulfasalazine), and TFAP2B and TFAP2D (17 alpha-ethinylestradiol). Thiram and ziram produced similar gene scores over six molecular targets. We did not observe an evident clustering pattern for MRDT versus male reproductive toxicants (Figure 5).
Literature Search for Male Reproductive Developmental Toxicants
Of the 212 male reproductive toxicants in ToxRefDB with ToxCastDB data (Figure 2B), we found 88 annotated by MeSH terms describing male reproductive effects in 2,099 articles in PubMed (see Excel File S1). We filtered this list by selecting chemicals with developmental toxicity and found 15 chemicals that were co-annotated with the MeSH terms "Cryptorchidism," "Hypospadias," and "Urogenital Abnormalities." Eight of these chemicals were not previously identified as MRDT chemicals in ToxRefDB, thus increasing the total number of MRDT chemicals from 63 to 71 (of which 54 had ToxCast[TM] data).
Human Homologs of ToxCast[TM] Assay Molecular Targets
Fifty-four MRDT chemicals from ToxRefDB and PubMed were mapped to 156 molecular targets in ToxCastDB (Figure 6). A list of the official gene symbols for these molecular targets was queried in the DAVID website (http://david.abcc.ncifcrf.gov/). Functional annotation by Homo sapiens captured 149 of the 156 gene symbols. Of the 7 IDs not captured or unmapped: a) CXCL8, formerly interleukin-8 (IL-8), was renamed by the Chemokine Nomenclature Subcommittee of the International Union of Immunological Societies, and the approved HUGO gene symbol is CXCL8; it functions as a major mediator of the inflammatory response and potent pro-angiogenic factor; b) CYP2C6V1 is an allelic variant of CYP2C19, which is associated with drug (anticonvulsant) metabolism; c) Cyp2a1 is the rat homolog of human CYP2A13 and functions in the steroid/lipid/cholesterol/retinol metabolic pathways; d) rat Cyp2a2 is a sexually dimorphic CYP (male-specific liver expression starting at puberty) that has steroid hydroxylase activity (e.g., testosterone 15-alphahydroxylase); homologs are Cyp2a4 (mouse) and CYP2A6, CYP2A7, and CYP2A13 (human); e) rat Cyp2c11 is a male-specific steroid (testosterone) alpha-hydroxylase homologous to CYP2C9 (human), which converts testosterone to 2-alpha-hydroxyl and 16-alpha-hydroxyl metabolites; f) rat Cyp2c12 is a female-predominant liver P450 that is active in the 15-beta-hydroxylation of steroid sulfates (prominent in rat preovulatory ovarian follicles), but no information is available on a human homolog; and g) Cyp3a2 is another male-specific steroid (testosterone) alpha-hydroxylase in rats, and in humans the CYP3A cluster consists of four genes (CYP3A43, CYP3A4, CYP3A7, and CYP3A5); these enzymes hydroxylate testosterone to less active metabolites.
Bipartite Network of Male Reproductive Developmental Toxicity
A bipartite network of the 54 MRDT chemicals and the 156 molecular targets in ToxCastDB was constructed in Gephi (Figure 6). In this network visualization of the target heat map (Figure 5), global inferences are extended to local relationships between different groups of molecular targets. Edge thickness of each chemical-target vector is proportional to the corresponding gene score, and the spatial distribution of nodes is based on connectivity. Nodes positioned centrally in the network represented the more commonly hit molecular targets: PPARs (peroxisome proliferator-activated receptors), RARs (retinoic acid receptors), RXRs (retinoid X receptors) in developmental regulation (yellow domain in Figure 6); as well as many nuclear receptors and crytochrome P450s in steroid homeostasis and xenobiotic metabolism (blue domain in Figure 6), such as NR1I2, ERs, AR, CYP19A1, and CYP3A5. Nodes positioned near the network's rim had limited connectivity, implying greater specificity: GPCRs in neural activity (green domain in Figure 6), inflammation and angiogenesis (purple and red domains).
Mining ToxRefDB revealed a phenotypic hierarchy of male reproductive developmental disorders. The pleiotropic in vitro bioactivity profiles suggested that complex arrays of molecular targets contributed to adverse outcomes of the developing male reproductive system in rodents. Estrogenic and anti-androgenic activities, although extensively studied in male reproductive toxicology, were only a fraction of the potential molecular targets found (e.g., nuclear receptors, GPCRs, vascular remodeling proteins, and cytochrome P450s). We conclude from this analysis that phenotypic features of a human TDS-like hierarchy inferred from the pattern of MRDT evident in ToxRefDB animal studies connect with multiple pathways in addition to estrogen and androgen, leading to adverse male developmental outcomes. This highlights a need for additional research, including use of computational systems-based models, to better characterize the mechanistic causes of male reproductive developmental disorders.
Different modes of action have been ascribed to various aspects of MRDT (Albert and Jegou 2014; Gray et al. 2006; Kapp 2010; Kongsbak et al. 2014; Wilson et al. 2008). The present study confirmed the connections in rodents to ESR1 and ESR2, and to a lesser degree AR. The inferred role of ER signaling is consistent with recent findings that estrogen synthesis is a key target of prenatal exposure to phthalates in the fetal rat testis proteome (Klinefelter et al. 2012; Veeramachaneni and Klinefelter 2014). Steroid hormone metabolism involves multiple nuclear receptors and cytochrome P450s, many of which were found to be molecular targets of MRDT chemicals. CYP1A1, CYP19A1, CYP3A4, and COMT are involved in the synthetic pathway of estrogen and testosterone; NR1I2 (PXR, pregnane X receptor) and NR1I3 (CAR, constitutive androstane receptor) regulate the transcription of cytochrome P450s, such as CYP3A4, CYP2Cs, and CYP2B6 (Wang and LeCluyse 2003; Willson and Kliewer 2002). NR1I2, NR1I3, as well as PPARs are all involved in fatty acid and steroid metabolism (Handschin and Meyer 2005). Although NR1I2 and NR1I3 are xenobiotic receptors investigated in cancer, drug interaction, obesity, and metabolic diseases (Chai et al. 2013; Gao and Xie 2012), in the present study we infer a potential connectivity with many of the MRDT chemicals.
Our results support a developmental connection to male reproductive developmental disorders in rodents beyond AR and ER (Cohn 2011; Kilcoyne et al. 2014; Kojima et al. 2010; van den Driesche et al. 2015). Retinoic acid signaling is involved in morphogenesis, growth, differentiation, and tissue homeostasis, which include genital tubercle development and different phases of spermatogenesis (Liu et al. 2012; Mark et al. 2006, 2015). RXRs are binding partners of PPARs, RARs, thyroid hormone receptor, and vitamin D receptor (Abbott 2009), all of which were located at the central region of the bipartite network. AHR (aryl hydrocarbon receptor), PPARG, and CEBPB are all highly expressed in the placenta, which represents a potential MRDT target given its role as an endocrine organ (Scott et al. 2009; Storvik et al. 2014). Crosstalk between the AHR and ESR1/ESR2 signaling pathways (Grans et al. 2010; Swedenborg et al. 2012) requires further investigation. Although AHR, ESR1, and ESR2 are promiscuous targets in the bipartite network, the measured AHR downstream targets (CYP1A1) and Per-Arnt-Sim (PAS) domain-containing proteins (HIF1A and KCNH2) show limited connectivity with one another, consistent with the recent molecular genetic screening studies (Geller et al. 2014; Madak-Erdogan et al. 2013). Further HTS experiments are needed to expand coverage for the TDS-related candidate genes, such as SHH, FGF10, BMP7, HOXA4, DGKK, and TGFBR3 (Carmichael et al. 2013; Dalgaard et al. 2012; Geller et al. 2014).
Vascularization plays an important role in the early stage of testis and placenta development (Brennan et al. 2002; McClelland et al. 2012; Reynolds and Redmer 2001). Although the testis and placenta play different roles in prenatal regulation of steroidogenesis in humans and rodents (Scott et al. 2009), the present study supports the hypothesis that vascular disruption is a key event in MRDT, similar to pregnancy loss and other birth defects (Abe et al. 2013; Kleinstreuer et al. 2011; Morita et al. 2014). Several molecular targets were identified in the IL-6 pathway (CSNK2A1, MAPK3, IL6, JUN, STAT3, and FOS), but previous knockout studies failed to show that the IL-6 pathway was critical to rodent embryo development (Sakurai et al. 2012, 2013; Yoshida et al. 2011). In comparison, the finding of molecular targets in the MAPK pathway (MAPK3, IL1A, JUN, NFKB1, TGFB1, TNF, TP53, and FOS) is consistent with previous toxicity studies in Sertoli cells (Qi et al. 2014) and embryo limb bud cells (Liu and Kapron 2010). Pathways underlying developmental angiogenesis and inflammation have crosstalk with one another as well as other biological processes regulated by growth factors, GPCRs, and nuclear receptors (Knudsen and Kleinstreuer 2011). For example, PPARs play a regulatory role in angiogenesis (BishopBailey 2011), and NR1H4 (FXR, farnesoid X receptor) down-regulates genes involved in inflammation (Hollman et al. 2012). These observations indicate multiple molecular initiating events and intermediate key events leading to MRDT adverse outcomes, supporting the use of the adverse outcome pathways (AOPs) approach (Ankley et al. 2010; Sturla et al. 2014) as a model framework to assess potential MRDT chemicals.
Dopamine, serotonin, and gonadotropinreleasing hormone receptors have been studied in the context of "neuroendocrine disruption" (Waye and Trudeau 2011). GPCR signaling is also involved in estradiol regulation of steroidogenesis (Vaucher et al. 2014). In the present study we identified several neuronal GPCRs targeted by a small number of MRDT chemicals in rodents, clustering at the lower rim of the bipartite network. Many of these chemicals are pharmaceutical agents that display strong neural or inflammatory activity, such as chlorpromazine hydrochloride (dopamine antagonist), SSR150106 (chemokine receptor antagonist), and SSR240612 (bradykinin B1 receptor antagonist). The molecular targets of neural activity were located apart from other predicted modes of action in the network graph, suggesting a unique bioactivity profile for a separate class of MRDT chemicals.
Unraveling the multiple pathways and processes underlying adverse prenatal developmental outcomes predicted from ToxCastDB-ToxRefDB signatures is a multi-scale problem (Sipes et al. 2011). The concept of bipartite networks (Fuxman Bass et al. 2013) has been used in the functional annotation of protein interaction networks (Kelley and Ideker 2005; Sharan et al. 2007) and in statistical modeling of complex social networks (Easley and Kleinberg 2010; Latapy et al. 2008). We applied the network concept to quantitatively visualize one-way interactions between chemical(s) and molecular target(s) in the ToxCastDB HTS data set and revealed novel interactions between different clusters of molecular targets across chemical classes, pathways, target score, and connectivity. This concept can also be applied for interpreting other high-dimensional data sets in computational toxicology.
The pattern of male reproductive effects in ToxRefDB varies by developmental stages and by test species, which is consistent with the differences in male reproductive development between rats, mice, rabbits, and humans. For instance, both mice and rats exposed to DBP exhibit male reproductive tract defects, including enlarged seminiferous cords and multinucleated germ cell formation, but fetal testis steroidogenesis is inhibited only in rats (Johnson et al. 2012). The human fetal testis is also susceptible to the seminiferous cord and germ cell effects of DBP, but it is relatively insensitive to the anti-steroidogenic effect (Spade et al. 2014; van den Driesche et al. 2015). The placenta plays an important role in regulating fetal testis steroidogenesis in humans, but not in rodents (Scott et al. 2009). Although rabbits have been reported as a sensitive model for many male reproductive effects (Foote and Carney 2000; Veeramachaneni 2008), almost all rabbit studies in ToxRefDB belong to the prenatal developmental study type and thus were not designed to capture male reproductive end points, particularly sperm effects. Therefore, it is not surprising that ToxRefDB records only two male reproductive toxicants in rabbits.
According to the TDS hypothesis, TDS pathogenesis begins during fetal development, involving symptoms of both newborn (cryptorchidism, hypospadias) and young adult males (impaired spermatogenesis, TGCT) (Akre and Richiardi 2009; Sharpe and Skakkebaek 2008). Thus, animal studies may miss the TDS-related effects of a chemical if they are not specifically designed to capture an extended scenario of early-life exposure and later-life effects, such as TGCT (a failure of some gonocytes to differentiate during postnatal/pubertal development). In contrast, testicular adenoma, interstitial cell tumor, and mesothelioma can all be detected in ToxRefDB (Table 1). These tumors originate from somatic cell mutations, potentially occurring at later life stages (Clegg et al. 1997; Meeks et al. 2012). On the other hand, the aggregate distribution of malformations, sperm effects, and histological effects were better correlated to the TDS hypothesis. Taken together, an extended one-generation reproductive toxicity study (OECD 2011) may be a suitable tool for validating present HTS findings for MRDT (Mitchard et al. 2012; Spielmann 2009). Sperm effects strongly correlate with malformations and histological effects, both of which are important indicators of early-stage testicular toxicity. Sensitivity of the rat as a species in detecting sperm effects and histological effects was stronger than the mouse based on number of positive chemicals and mean cLELs. With the caveat that poor semen quality assessed in adult men cannot be attributed to early-life perturbation in the absence of sufficient maternal exposure history, this supports the case for using sperm effects as a noninvasive assessment to improve surveillance of human reproductive function at a population level (CDC 2014).
In conclusion, ToxRefDB and ToxCastDB provide a novel resource for exploring phenotypic linkages of MRDT to potential mechanistic pathways. This study reveals a phenotypic hierarchy of male reproductive effects across animal studies that were consistent with the human TDS hypothesis. Although a number of recent studies have demonstrated that both hormonal and nonhormonal signaling is involved in male reproductive developmental defects (Carmichael et al. 2013; Dalgaard et al. 2012; Geller et al. 2014), the present study extends these findings, showing androgen and estrogen receptors to be a subset of the potential landscape of molecular targets in rodents. Steroid homeostasis and xenobiotic metabolism, developmental regulation, vascularization and angiogenesis, inflammation, and neural activity provide a more comprehensive set of predicted modes of action underlying the TDS hypothesis. Further multi-scale modeling of MRDT will be needed to capture the spatiotemporal dynamics of male reproductive tract development and disruption. In combination with further HTS experiments that cover more molecular targets in embryonic development, this approach provides novel mechanistic insights for male reproductive developmental disorders.
Abbott BD. 2009. Review of the expression of peroxisome proliferator-activated receptors alpha (PPAR[alpha]), beta (PPAR[beta]), and gamma (PPAR[gamma]) in rodent and human development. Reprod Toxicol 27:246-257.
Abe N, Nakahara T, Morita A, Wada Y, Mori A, Sakamoto K, et al. 2013. KRN633, an inhibitor of vascular endothelial growth factor receptor tyrosine kinase, induces intrauterine growth restriction in mice. Birth Defects Res B Dev Reprod Toxicol 98:297-303.
Akre O, Richiardi L. 2009. Does a testicular dysgenesis syndrome exist? Hum Reprod 24:2053-2060.
Albert O, Jegou B. 2014. A critical assessment of the endocrine susceptibility of the human testis to phthalates from fetal life to adulthood. Hum Reprod Update 20:231-249.
Ankley GT, Bennett RS, Erickson RJ, Hoff DJ, Hornung MW, Johnson RD, et al. 2010. Adverse outcome pathways: a conceptual framework to support ecotoxicology research and risk assessment. Environ Toxicol Chem 29:730-741.
Aschim EL, Nordenskjold A, Giwercman A, Lundin KB, Ruhayel Y, Haugen TB, et al. 2004. Linkage between cryptorchidism, hypospadias, and GGN repeat length in the androgen receptor gene. J Clin Endocrinol Metab 89:5105-5109.
Bastian M, Heymann S, Jacomy M. 2009. Gephi: an open source software for exploring and manipulating networks. In: Proceedings of the Third International Conference on Weblogs and Social Media, ICWSM, 17-20 May 2009, San Jose, California, USA.
Bay K, Asklund C, Skakkebaek NE, Andersson AM. 2006. Testicular dysgenesis syndrome: possible role of endocrine disrupters. Best Pract Res Clin Endocrinol Metab 20:77-90.
Biggs ML, Baer A, Critchlow CW. 2002. Maternal, delivery, and perinatal characteristics associated with cryptorchidism: a population-based case-control study among births in Washington State. Epidemiology 13:197-204.
Bishop-Bailey D. 2011. PPARs and angiogenesis. Biochem Soc Trans 39:1601-1605.
Brennan J, Karl J, Capel B. 2002. Divergent vascular mechanisms downstream of Sry establish the arterial system in the XY gonad. Dev Biol 244:418-428.
Carmichael SL, Ma C, Choudhry S, Lammer EJ, Witte JS, Shaw GM. 2013. Hypospadias and genes related to genital tubercle and early urethral development. J Urol 190:1884-1892.
CDC (Centers for Disease Control and Prevention). 2014. National Public Health Action Plan for the Detection, Prevention, and Management of Infertility. Atlanta, GA:CDC. Available: http:// www.cdc.gov/reproductivehealth/Infertility/PDF/ DRH_NAP_Final_508.pdf [accessed 24 June 2015].
Chai X, Zeng S, Xie W. 2013. Nuclear receptors PXR and CAR: implications for drug metabolism regulation, pharmacogenomics and beyond. Expert Opin Drug Metab Toxicol 9:253-266.
Clegg ED, Cook JC, Chapin RE, Foster PM, Daston GP. 1997. Leydig cell hyperplasia and adenoma formation: mechanisms and relevance to humans. Reprod Toxicol 11:107-121.
Cohn MJ. 2011. Development of the external genitalia: conserved and divergent mechanisms of appendage patterning. Dev Dyn 240:1108-1115.
Collins FS, Gray GM, Bucher JR. 2008. Toxicology. Transforming environmental health protection. Science 319:906-907.
Dalgaard MD, Weinhold N, Edsgard D, Silver JD, Pers TH, Nielsen JE, et al. 2012. A genome-wide association study of men with symptoms of testicular dysgenesis syndrome and its network biology interpretation. J Med Genet 49:58-65.
Easley DA, Kleinberg JM. 2010. Networks, Crowds, and Markets: Reasoning About a Highly Connected World. Cambridge, UK:Cambridge University Press.
Foote RH, Carney EW. 2000. The rabbit as a model for reproductive and developmental toxicity studies. Reprod Toxicol 14:477-493.
Fuxman Bass JI, Diallo A, Nelson J, Soto JM, Myers CL, Walhout AJ. 2013. Using networks to measure similarity between genes: association index selection. Nat Methods 10:1169-1176.
Gao J, Xie W. 2012. Targeting xenobiotic receptors PXR and CAR for metabolic diseases. Trends Pharmacol Sci 33:552-558.
Ge RS, Chen GR, Dong Q, Akingbemi B, Sottas CM, Santos M, et al. 2007. Biphasic effects of postnatal exposure to diethylhexylphthalate on the timing of puberty in male rats. J Androl 28:513-520.
Geller F, Feenstra B, Carstensen L, Pers TH, van Rooij IA, Korberg IB, et al. 2014. Genome-wide association analyses identify variants in developmental genes associated with hypospadias. Nat Genet 46:957-963.
Grans J, Wassmur B, Celander MC. 2010. One-way inhibiting cross-talk between arylhydrocarbon receptor (AhR) and estrogen receptor (ER) signaling in primary cultures of rainbow trout hepatocytes. Aquat Toxicol 100:263-270.
Gray LE Jr, Wilson VS, Stoker T, Lambright C, Furr J, Noriega N, et al. 2006. Adverse effects of environmental antiandrogens and androgens on reproductive development in mammals. Int J Androl 29:96-104.
Handschin C, Meyer UA. 2005. Regulatory network of lipid-sensing nuclear receptors: roles for CAR, PXR, LXR, and FXR. Arch Biochem Biophys 433:387-396.
Harding HP, Zhang Y, Zeng H, Novoa I, Lu PD, Calfon M, et al. 2003. An integrated stress response regulates amino acid metabolism and resistance to oxidative stress. Mol Cell 11:619-633.
Hollman DA, Milona A, van Erpecum KJ, van Mil SW. 2012. Anti-inflammatory and metabolic actions of FXR: insights into molecular mechanisms. Biochim Biophys Acta 1821:1443-1452.
Hu Y, Wang R, Xiang Z, Qian W, Han X, Li D. 2014. Antagonistic effects of a mixture of low-dose nonylphenol and di-n-butyl phthalate (monobutyl phthalate) on the Sertoli cells and serum reproductive hormones in prepubertal male rats in vitro and in vivo. PLoS One 9:e93425, doi:10.1371/journal. pone.0093425.
Huang DW, Sherman BT, Lempicki RA. 2009a. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res 37:1-13.
Huang DW, Sherman BT, Lempicki RA. 2009b. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc 4:44-57.
Jain VG, Singal AK. 2013. Shorter anogenital distance correlates with undescended testis: a detailed genital anthropometric analysis in human newborns. Hum Reprod 28:2343-2349.
Johnson KJ, Heger NE, Boekelheide K. 2012. Of mice and men (and rats): phthalate-induced fetal testis endocrine disruption is species-dependent. Toxicol Sci 129:235-248.
Johnson KJ, Robbins AK, Wang Y, McCahan SM, Chacko JK, Barthold JS. 2010. Insulin-like 3 exposure of the fetal rat gubernaculum modulates expression of genes involved in neural pathways. Biol Reprod 83:774-782.
Kaleva M, Toppari J. 2005. Cryptorchidism: an indicator of testicular dysgenesis? Cell Tissue Res 322:167-172.
Kapp RW Jr. 2010. Male reproductive toxicity. In: Reproductive Toxicology (Kapp RW Jr, Tyl RW, eds). 3rd ed. New York:Informa Healthcare, 249-269.
Kavlock R, Chandler K, Houck K, Hunter S, Judson R, Kleinstreuer N, et al. 2012. Update on EPA's ToxCast program: providing high throughput decision support tools for chemical risk management. Chem Res Toxicol 25:1287-1302.
Kelley R, Ideker T. 2005. Systematic interpretation of genetic interactions using protein networks. Nat Biotechnol 23:561-566.
Kilcoyne KR, Smith LB, Atanassova N, Macpherson S, McKinnell C, van den Driesche S, et al. 2014. Fetal programming of adult Leydig cell function by androgenic effects on stem/progenitor cells. Proc Natl Acad Sci USA 111:E1924-E1932.
Kleinstreuer NC, Judson RS, Reif DM, Sipes NS, Singh AV, Chandler KJ, et al. 2011. Environmental impact on vascular development predicted by high-throughput screening. Environ Health Perspect 119:1596-1603, doi:10.1289/ehp.1103412.
Klinefelter GR, Laskey JW, Winnik WM, Suarez JD, Roberts NL, Strader LF, et al. 2012. Novel molecular targets associated with testicular dysgenesis induced by gestational exposure to diethylhexyl phthalate in the rat: a role for estradiol. Reproduction 144:747-761.
Knudsen TB, Kleinstreuer NC. 2011. Disruption of embryonic vascular development in predictive toxicology. Birth Defects Res C Embryo Today 93:312-323.
Knudsen TB, Martin MT, Kavlock RJ, Judson RS, Dix DJ, Singh AV. 2009. Profiling the activity of environmental chemicals in prenatal developmental toxicity studies using the U.S. EPA's ToxRefDB. Reprod Toxicol 28:209-219.
Kojima Y, Kohri K, Hayashi Y. 2010. Genetic pathway of external genitalia formation and molecular etiology of hypospadias. J Pediatr Urol 6:346-354.
Kongsbak K, Vinggaard AM, Hadrup N, Audouze K. 2014. A computational approach to mechanistic and predictive toxicology of pesticides. ALTEX 31:11-22.
Latapy M, Magnien C, Del Vecchio N. 2008. Basic notions for the analysis of large two-mode networks. Soc Networks 30:31-48.
Laws SC, Carey SA, Ferrell JM, Bodman GJ, Cooper RL. 2000. Estrogenic activity of octylphenol, nonylphenol, bisphenol A and methoxychlor in rats. Toxicol Sci 54:154-167.
Liu J, Kapron CM. 2010. Differential induction of MAP kinase signalling pathways by cadmium in primary cultures of mouse embryo limb bud cells. Reprod Toxicol 29:286-291.
Liu L, Suzuki K, Nakagata N, Mihara K, Matsumaru D, Ogino Y, et al. 2012. Retinoic acid signaling regulates sonic hedgehog and bone morphogenetic protein signaling during genital tubercle development. Birth Defects Res B Dev Reprod Toxicol 95:79-88.
Lottrup G, Jorgensen A, Nielsen JE, Jorgensen N, Duno M, Vinggaard AM, et al. 2013. Identification of a novel androgen receptor mutation in a family with multiple components compatible with the testicular dysgenesis syndrome. J Clin Endocrinol Metab 98:2223-2229.
Madak-Erdogan Z, Charn TH, Jiang Y, Liu ET, Katzenellenbogen JA, Katzenellenbogen BS. 2013. Integrative genomics of gene and metabolic regulation by estrogen receptors a and 3, and their coregulators. Mol Syst Biol 9:676, doi:10.1038/msb.2013.28.
Mark M, Ghyselinck NB, Chambon P. 2006. Function of retinoid nuclear receptors: lessons from genetic and pharmacological dissections of the retinoic acid signaling pathway during mouse embryogenesis. Annu Rev Pharmacol Toxicol 46:451-480.
Mark M, Teletin M, Vernet N, Ghyselinck NB. 2015. Role of retinoic acid receptor (RAR) signaling in post-natal male germ cell differentiation. Biochim Biophys Acta 1849:84-93.
Martin MT, Judson RS, Reif DM, Kavlock RJ, Dix DJ. 2009a. Profiling chemicals based on chronic toxicity results from the U.S. EPA ToxRef Database. Environ Health Perspect 117:392-399, doi:10.1289/ ehp.0800074.
Martin MT, Mendez E, Corum DG, Judson RS, Kavlock RJ, Rotroff DM, et al. 2009b. Profiling the reproductive toxicity of chemicals from multigeneration studies in the toxicity reference database. Toxicol Sci 110:181-190.
McClelland K, Bowles J, Koopman P. 2012. Male sex determination: insights into molecular mechanisms. Asian J Androl 14:164-171.
Meeks JJ, Sheinfeld J, Eggener SE. 2012. Environmental toxicology of testicular cancer. Urol Oncol 30:212-215.
Mitchard T, Jarvis P, Stewart J. 2012. Assessment of male rodent fertility in general toxicology 6-month studies. Birth Defects Res B Dev Reprod Toxicol 95:410-420.
Mitchell RT, Sharpe RM, Anderson RA, McKinnell C, Macpherson S, Smith LB, et al. 2013. Diethylstilboestrol exposure does not reduce testosterone production in human fetal testis xenografts. PLoS One 8:e61726, doi:10.1371/journal.pone.0061726.
Morita A, Nakahara T, Abe N, Kurauchi Y, Mori A, Sakamoto K, et al. 2014. Treatment of mid-pregnant mice with KRN633, an inhibitor of vascular endothelial growth factor receptor tyrosine kinase, induces abnormal retinal vascular patterning in their newborn pups. Birth Defects Res B Dev Reprod Toxicol 101:293-299.
Mylchreest E, Cattley RC, Foster PM. 1998. Male reproductive tract malformations in rats following gestational and lactational exposure to di(n-butyl) phthalate: an antiandrogenic mechanism? Toxicol Sci 43:47-60.
Mylchreest E, Sar M, Wallace DG, Foster PM. 2002. Fetal testosterone insufficiency and abnormal proliferation of Leydig cells and gonocytes in rats exposed to di(n-butyl) phthalate. Reprod Toxicol 16:19-28.
Mylchreest E, Wallace DG, Cattley RC, Foster PM. 2000. Dose-dependent alterations in androgenregulated male reproductive development in rats exposed to di(n-butyl) phthalate during late gestation. Toxicol Sci 55:143-151.
NRC (National Research Council). 2007. Toxicity Testing in the 21st Century: a Vision and a Strategy. Washington, DC:National Academies Press.
OECD (Organization for Economic Co-operation and Development). 2011. Test No. 443: Extended One-Generation Reproductive Toxicity Study, OECD Guidelines for the Testing of Chemicals, Section 4. Paris, France:OECD Publishing. Available: http:// dx.doi.org/10.1787/9789264122550-en [accessed 24 June 2015].
Ono M, Harley VR. 2013. Disorders of sex development: new genes, new concepts. Nat Rev Endocrinol 9:79-91.
Parks LG, Ostby JS, Lambright CR, Abbott BD, Klinefelter GR, Barlow NJ, et al. 2000. The plasticizer diethylhexyl phthalate induces malformations by decreasing fetal testosterone synthesis during sexual differentiation in the male rat. Toxicol Sci 58:339-349.
Qi S, Fu W, Wang C, Liu C, Quan C, Kourouma A, et al. 2014. BPA-induced apoptosis of rat Sertoli cells through Fas/FasL and JNKs/p38 MAPK pathways. Reprod Toxicol 50:108-116.
R Core Team. 2013. R: A Language and Environment for Statistical Computing. Vienna, Austria:R Foundation for Statistical Computing. Available: http://www.R-project.org/ [accessed 16 November 2015].
Reynolds LP, Redmer DA. 2001. Angiogenesis in the placenta. Biol Reprod 64:1033-1040.
Sakurai T, Takai R, Burgin H, Ishihara K, Sakamoto Y, Amano J, et al. 2012. The effects of interleukin-6 signal blockade on fertility, embryo-fetal development, and immunization in vivo. Birth Defects Res B Dev Reprod Toxicol 95:304-317.
Sakurai T, Takai R, Burgin H, Shioda A, Sakamoto Y, Amano J, et al. 2013. The effects of interleukin-6 signal blockade on immune system, reproductive and skeletal development in juvenile mice. Birth Defects Res B Dev Reprod Toxicol 98:170-182.
Scott HM, Hutchison GR, Mahood IK, Hallmark N, Welsh M, De Gendt K, et al. 2007. Role of androgens in fetal testis development and dysgenesis. Endocrinology 148:2027-2036.
Scott HM, Mason JI, Sharpe RM. 2009. Steroidogenesis in the fetal testis and its susceptibility to disruption by exogenous compounds. Endocr Rev 30:883-925.
Sharan R, Ulitsky I, Shamir R. 2007. Network-based prediction of protein function. Mol Syst Biol 3:88, doi:10.1038/msb4100129.
Sharpe RM, Skakkebaek NE. 2008. Testicular dysgenesis syndrome: mechanistic insights and potential new downstream effects. Fertil Steril 89:e33-e38.
Sipes NS, Martin MT, Reif DM, Kleinstreuer NC, Judson RS, Singh AV, et al. 2011. Predictive models of prenatal developmental toxicity from ToxCast high-throughput screening data. Toxicol Sci 124:109-127.
Skakkebaek NE, Rajpert-De Meyts E, Main KM. 2001. Testicular dysgenesis syndrome: an increasingly common developmental disorder with environmental aspects. Hum Reprod 16:972-978.
Spade DJ, Hall SJ, Saffarini CM, Huse SM, McDonnell EV, Boekelheide K. 2014. Differential response to abiraterone acetate and di-n-butyl phthalate in an androgen-sensitive human fetal testis xenograft bioassay. Toxicol Sci 138:148-160.
Spielmann H. 2009. The way forward in reproductive/ developmental toxicity testing. Altern Lab Anim 37:641-656.
Storvik M, Huuskonen P, Pehkonen P, Pasanen M. 2014. The unique characteristics of the placental transcriptome and the hormonal metabolism enzymes in placenta. Reprod Toxicol 47:9-14.
Sturla SJ, Boobis AR, FitzGerald RE, Hoeng J, Kavlock RJ, Schirmer K, et al. 2014. Systems toxicology: from basic research to risk assessment. Chem Res Toxicol 27:314-329.
Swedenborg E, Kotka M, Seifert M, Kanno J, Pongratz I, Ruegg J. 2012. The aryl hydrocarbon receptor ligands 2,3,7,8-tetrachlorodibenzo-pdioxin and 3-methylcholanthrene regulate distinct genetic networks. Mol Cell Endocrinol 362:39-47.
Thankamony A, Lek N, Carroll D, Williams M, Dunger DB, Acerini CL, et al. 2014. Anogenital distance and penile length in infants with hypospadias or cryptorchidism: comparison with normative data. Environ Health Perspect 122:207-211, doi:10.1289/ehp.1307178.
Trabert B, Zugna D, Richiardi L, McGlynn KA, Akre O. 2013. Congenital malformations and testicular germ cell tumors. Int J Cancer 133:1900-1904.
van den Driesche S, McKinnell C, Calarrao A, Kennedy L, Hutchison GR, Hrabalkova L, et al. 2015. Comparative effects of di(n-butyl) phthalate exposure on fetal germ cell development in the rat and in human fetal testis xenografts. Environ Health Perspect 123:223-230, doi:10.1289/ehp.1408248.
Vaucher L, Funaro MG, Mehta A, Mielnik A, Bolyakov A, Prossnitz ER, et al. 2014. Activation of GPER-1 estradiol receptor downregulates production of testosterone in isolated rat Leydig cells and adult human testis. PLoS One 9:e92425, doi:10.1371/journal.pone.0092425.
Veeramachaneni DN. 2008. Impact of environmental pollutants on the male: effects on germ cell differentiation. Anim Reprod Sci 105:144-157.
Veeramachaneni DN, Klinefelter GR. 2014. Phthalate-induced pathology in the foetal testis involves more than decreased testosterone production. Reproduction 147:435-442.
Virtanen HE, Rajpert-De Meyts E, Main KM, Skakkebaek NE, Toppari J. 2005. Testicular dysgenesis syndrome and the development and occurrence of male reproductive disorders. Toxicol Appl Pharmacol 207(2 suppl):501-505.
Wang H, LeCluyse EL. 2003. Role of orphan nuclear receptors in the regulation of drug-metabolising enzymes. Clin Pharmacokinet 42:1331-1357.
Ward JH. 1963. Hierarchical grouping to optimize an objective function. J Am Stat Assoc 58:236-244.
Waye A, Trudeau VL. 2011. Neuroendocrine disruption: more than hormones are upset. J Toxicol Environ Health B Crit Rev 14:270-291.
Willson TM, Kliewer SA. 2002. PXR, CAR and drug metabolism. Nat Rev Drug Discov 1:259-266.
Wilson VS, Blystone CR, Hotchkiss AK, Rider CV, Gray LE Jr. 2008. Diverse mechanisms of anti-androgen action: impact on male rat reproductive tract development. Int J Androl 31:178-187.
Woodruff TJ, Zota AR, Schwartz JM. 2011. Environmental chemicals in pregnant women in the United States: NHANES 2003-2004. Environ Health Perspect 119:878-885, doi:10.1289/ehp.1002727.
Yoshida H, Hashizume M, Suzuki M, Mihara M. 2011. Induction of high-dose tolerance to the rat anti-mouse IL-6 receptor antibody in NZB/NZW F1 mice. Rheumatol Int 31:1445-1449.
Maxwell C.K. Leung, (1,2) Jimmy Phuong, (2) Nancy C. Baker, (3) Nisha S. Sipes, (12) Gary R. Klinefelter, (4) Matthew T. Martin, (2) Keith W. McLaurin, (1,2) R. Woodrow Setzer, (2) Sally Perreault Darney, (4)* Richard S. Judson, (2) and Thomas B. Knudsen (2)
(1) Oak Ridge Institute for Science and Education, Oak Ridge, Tennessee; (2) National Center for Computational Toxicology, U.S. Environmental Protection Agency (EPA), Research Triangle Park, North Carolina; (3) Lockheed Martin, Research Triangle Park, North Carolina; (4) National Health and Environmental Effects Research Laboratory, U.S. EPA, Research Triangle Park, North Carolina
Address correspondence to M.C.K. Leung, U.S. Environmental Protection Agency, 109 T.W. Alexander Dr., Research Triangle Park, NC 27711 USA. Telephone: (919) 541-2721. E-mail: email@example.com, or T.B. Knudsen, U.S. Environmental Protection Agency, 109 T.W. Alexander Dr., Research Triangle Park, NC 27711 USA. Telephone: (919) 541-9776. E-mail: firstname.lastname@example.org
* Current address for S.P.D.: Environmental Health Perspectives, National Institute of Environmental Health Sciences, National Institutes of Health, Department of Health and Human Services, Research Triangle Park, NC.
Supplemental Material is available online (http:// dx.doi.org/10.1289/ehp.1510385).
We thank D.F. Kapraun, L. Truong, and A.M. Zoss for technical consultation; A.L. Karmaus, A.M. Frame, and P. Volarath for database queries; and A.H. Piersma for technical review.
N.C.B. is employed by Lockheed Martin, Research Triangle Park, North Carolina.
The views expressed in this article are those of the authors and do not necessarily reflect the views or policies of the U.S. EPA. Mention of trade names or commercial products does not constitute endorsement or recommendation for use.
The authors declare they have no actual or potential competing financial interests.
Received: 24 June 2015; Accepted: 24 November 2015; Advance Publication: 11 December 2015; Final Publication: 1 July 2016.
Caption: Figure 1. In vivo phenotypic analysis: chemical coverage and workflow. (A) Experimental protocols in general follow EPA Health Effects Test Guidelines OPPTS (a) 870.3700, (b) 870.3100, 870.4100, 870.4200, 870.4300, (c) 870.3800, and (d) 870.3550. Test strains include (e) Alderly Park, CD, Fischer 344, Long-Evans, Sprague-Dawley, and Wistar for rats; (f) B6C3F1, CD, and Swiss for mice; and (g) Chinchilla, Dutch, Himalayan, and New Zealand White for rabbits. The prenatal developmental study type (OPPTS 870.3700 or equivalent) consisted mostly of rabbit (527) and rat studies (731), whereas both one-generational and multigenerational reproductive study types (OPPTS 870.3800, 870.3550, or equivalent) consisted primarily of rat studies (104 and 455, respectively). The ToxRefDB subchronic and chronic study types (OPPT 870.3100, 870.4100, 870.4200, 870.4300, or equivalent) consisted mostly of mouse (568) and rat studies (984). (B) Workflow for data-mining ToxRefDB.
Caption: Figure 2. In vitro network analysis: workflow and evidence of effects. (A) Workflow for data-mining ToxCastDB. (B) A total of 2,099 articles were co-annotated with 88 of the 212 male reproductive toxicants in ToxRefDB (Figure 5), including 15 chemicals co-annotated with cryptorchidism, hypospadias, and urogenital abnormalities. The citations were listed in Excel File S1.
Caption: Figure 4. Phenotypic hierarchy of 63 male reproductive developmental toxicants in ToxRefDB. An unsupervised two-directional heat map was constructed based on hierarchical clustering of 63 chemicals by their effect scores for malformations, sperm effects, histological effects, and relative weight change in male reproductive tissues, using Euclidean distance for measure and Ward's method for linkage analysis.
Caption: Figure 5. Unsupervised two-directional heatmap from hierarchical clustering of 212 male reproductive toxicants in ToxCastDB. Data for 212 male reproductive toxicants in ToxCastDB were used to construct a heat map to visualize the complex chemical-effects relationships on 213 molecular targets. The heat map was generated with the "gplot" package in R version 3.1.0, using Euclidean distance for measure and Ward's method for linkage analysis. The intensity of blue color indicated the value of gene score.
Caption: Figure 6. Quantitative bipartite network of 54 male reproductive developmental toxicants and 156 molecular targets in ToxCastDB. A bipartite network consisting of 54 male reproductive developmental toxicants (source nodes) and 156 molecular targets (target nodes; bold) was constructed with ForceAtlas layout using Gephi version 0.8.2 beta, with 5 predicted modes of action in different colors. The thickness of each chemical-to- molecular target arrow was proportional to the corresponding gene score.
Table 1. ToxRefDB positive chemicals for male reproductive effects. Number of positive chemicals (mean lowest effect level; cLEL) (mg/kg/day) Male reproductive effect Rabbit Mouse Rat P1 Malformation 1 0 14 NA Hypospadias 0 0 2 (98) NA Nipple retention 0 0 3 (70) NA Decreased anogenital 0 0 9 (57) NA distance Cryptorchidism 1 (80) 0 7(241) NA Testicular tumor 0 1 27 28 Mesothelioma 0 0 5 (6) 5 (6) Adenoma 0 0 8 (25) 8 (25) Interstitial cell tumor 0 1 (113) 15 (35) 16 (38) Other 0 0 1 (11) 1 (11) Sperm effect 1 15 42 34 Morphology 1 (3) 3 (39) 10 (64) 9 (67) Motility 0 8 (102) 21 (44) 21 (56) Testicular spermatid 0 4 (40) 9 (84) 8 (111) count Epididymal sperm count 0 4 (79) 12 (64) 7 (29) Semen sperm count 0 5 (58) 11 (90) 12 (65) Histological effect 0 34 95 103 Prostate 0 2 (398) 17 (30) 15 (44) Seminal vesicle 0 3 (158) 8 (25) 10 (41) Epididymis 0 6 (153) 25 (20) 25 (23) Epididymis 0 5 (62) 22 (67) 22 (68) (spermatogenesis- specific) Testis 0 11 (159) 29 (56) 35 (69) Testis (spermatogenesis- 0 23 (96) 65 (56) 71 (57) specific) Relative weight change in 0 42 176 174 tissue Prostate 0 3 (342) 40 (64) 29 (100) Seminal vesicle 0 4 (274) 36 (41) 29 (36) Epididymis 0 5 (294) 34 (54) 25 (45) Testis 0 36 (88) 153 (56) 154 (56) Total 2 64 224 221 Number of positive chemicals (mean lowest effect level; cLEL) (mg/kg/day) Male reproductive effect F1 P1 or F1 MRDT Malformation 12 12 12 Hypospadias 2 (98) 2 2 Nipple retention 3 (70) 3 3 Decreased anogenital 8 (232) 8 8 distance Cryptorchidism 6 (285) 6 6 Testicular tumor NA 28 NA Mesothelioma NA 5 NA Adenoma NA 8 NA Interstitial cell tumor NA 16 NA Other NA 1 NA Sperm effect 24 51 36 Morphology 10 (42) 14 12 Motility 6 (31) 25 14 Testicular spermatid 5 (58) 12 11 count Epididymal sperm count 10 (159) 15 13 Semen sperm count 4 (36) 15 8 Histological effect 27 111 20 Prostate 6 (31) 19 7 Seminal vesicle 1 (50) 11 2 Epididymis 9 (17) 27 7 Epididymis 4 (50) 26 8 (spermatogenesis- specific) Testis 4 (127) 39 8 Testis (spermatogenesis- 17 (32) 77 16 specific) Relative weight change in 58 196 50 tissue Prostate 16 (42) 42 16 Seminal vesicle 14 (138) 39 20 Epididymis 21 (81) 37 27 Testis 47 (83) 172 36 Total 84 248 63 NA, not applicable. Includes 248 of the 281 chemicals positive for male reproductive effects with dosage Information for prenatal and generational (P1, F1) studies (F2 and F3 data not shown); mean cLELs are listed in parentheses. A chemical was found to produce a male reproductive effect if produced any of the listed search terms in ToxRefDB in Excel File S2, regardless of dosage. MRDT (male reproductive developmental toxicity): sperm effects, histological effects, and tissue weight change were reported exclusively in developmental studies, plus all malformations; tumor location and methods used for sperm count were not specified in ToxRefDB; general histological effects include inflammation, necrosis, and pigmentation, in contrast with spermatogenesis-specific histological effects such as aspermia, immature sperm, and retained spermatids. Figure 3. Co-occurrence of male reproductive effects in ToxRefDB. A chemical was found to produce a male reproductive effect if it was associated with any of the listed search terms in Excel File S2 in ToxRefDB, regardless of dosage. Chemicals may be associated with multiple effects. Chemicals were ranked in the order of association with the least common effect (i.e., malformations), followed by the next (i.e., testicular tumors), and so on. Chemical numbers Malformation 18 (2%) Testicular tumor 31 (4%) Sperm effect 58 (8%) Histological effect 126 (16%) Relative weight change in tissue 227 (29%) Note: Table made from bar graph.
Please note: Illustration(s) are not available due to copyright restrictions.
|Printer friendly Cite/link Email Feedback|
|Author:||Leung, Maxwell C.K.; Phuong, Jimmy; Baker, Nancy C.; Sipes, Nisha S.; Klinefelter, Gary R.; Martin,|
|Publication:||Environmental Health Perspectives|
|Date:||Jul 1, 2016|
|Previous Article:||Ingested nitrate and breast cancer in the Spanish multicase-control study on cancer (MCC-Spain).|
|Next Article:||The genetic basis for variation in sensitivity to lead toxicity in Drosophila melanogaster.|