Scholarly article on topic 'An embryonic atrazine exposure results in reproductive dysfunction in adult zebrafish and morphological alterations in their offspring'

An embryonic atrazine exposure results in reproductive dysfunction in adult zebrafish and morphological alterations in their offspring Academic research paper on "Biological sciences"

0
0
Share paper
Academic journal
Sci. Rep.
OECD Field of science
Keywords
{""}

Academic research paper on topic "An embryonic atrazine exposure results in reproductive dysfunction in adult zebrafish and morphological alterations in their offspring"

SCIENTIFIC REPpRTS

Received: 01 December 2015 Accepted: 21 January 2016 Published: 19 February 2016

An embryonic atrazine exposure results in reproductive dysfunction in adult zebrafish and morphological alterations in their offspring

Sara E. Wirbisky1, Gregory J. Weber1, Maria S. Sepulveda1,2, Tsang-Long Lin3, Amber S. Jannasch4 & Jennifer L. Freeman1

The herbicide atrazine, a suspected endocrine disrupting chemical (EDC), frequently contaminates potable water supplies. Studies suggest alterations in the neuroendocrine system along the hypothalamus-pituitary-gonadal axis; however, most studies address either developmental, pubertal, or adulthood exposures, with few investigations regarding a developmental origins hypothesis. In this study, zebrafish were exposed to 0, 0.3, 3, or 30 parts per billion (ppb) atrazine through embryogenesis and then allowed to mature with no additional chemical exposure. Reproductive function, histopathology, hormone levels, offspring morphology, and the ovarian transcriptome were assessed. Embryonic atrazine exposure resulted in a significant increase in progesterone levels in the 3 and 30 ppb groups. A significant decrease in spawning and a significant increase in follicular atresia in the 30 ppb group were observed. In offspring, a decrease in the head length to body ratio in the 30 ppb group, along with a significant increase in head width to body ratio in the 0.3 and 3 ppb groups occurred. Transcriptomic alterations involved genes associated with endocrine system development and function, tissue development, and behavior. This study provides evidence to support atrazine as an EDC causing reproductive dysfunction and molecular alterations in adults exposed only during embryogenesis and morphological alterations in their offspring.

Studies investigating the effects of early life exposure to environmental stressors or stimuli have increased dramatically over the past decade. These studies seek to investigate the developmental origin of health and adult disease (DOHaD) hypothesis which states that exposure to stressors during sensitive times during an organism's life, specifically during developmental stages, can cause changes to the genome and epigenome thereby resulting in an increased susceptibility to the development of health issues or diseases later on in life1'2. A key element complicating the establishment of a link between exposure and a disease state is the time that elapses between exposure and outward response or development of a disease1'3. Thus, it may take years for an individual to present a disease state and in addition may pass on these adverse health effects to future generations4.

Endocrine disrupting chemicals (EDCs) are exogenous agents that alter endocrine system functions and are associated with a myriad of diseases. In recent years, public concern about the effects of EDCs on human health has increased substantially and heightened the need for further research into the underlying molecular mechanisms of toxicity of these compounds5,6. EDCs are diverse in structure and are present in many products such as pharmaceuticals, plasticizers, and pesticides, making human exposure to these potentially harmful compounds a likely event. Evidence suggests that EDCs do not adhere to classic dose-response toxicological principles; rather they are part of the 'low dose hypothesis' due to their ability to disrupt hormonal homeostasis at low concentrations7. Studies show that EDCs can cause irreversible changes in tissue formation, decreased reproductive

1School of Health Sciences, West Lafayette, IN, 47907, USA. 2Department of Forestry and Natural Resources, West

Lafayette, IN, 47907, USA. 3Department of Comparative Pathobiology, West Lafayette, IN, 47907, USA. 4Bindley

Bioscience Center, Purdue University, West Lafayette, IN, 47907, USA. Correspondence and requests for materials

should be addressed to J.L.F. (email: jfreema@purdue.edu)

potential, obesity, and cancer8-12. Moreover, evidence suggests that exposure to EDCs can cause adverse effects not only in organisms that come into contact with them, but also to future progeny of exposed individuals13.

Atrazine (2-chloro-4-ethylamino-6-isopropylamino-1,3,5-triazine) is a pre-emergent herbicide predominately used in the Midwestern United States to control broadleaf and grassy weeds on a variety of field crops14. Exposure to atrazine can occur through different routes including ingestion of contaminated drinking water and in occupational settings via inhalation15-17. Furthermore, atrazine is often reported to exceed the Maximum Contaminant Level (MCL) of 3 parts per billion (ppb; |i g/L) set by the U.S. Environmental Protection Agency (EPA) in potable water supplies18'19. As such the European Union banned the use of atrazine in 200320'21. Epidemiological studies show several potential adverse health effects associated with maternal atrazine exposure including an increased risk of babies born small for their gestational age (SGA), intrauterine growth retardation (IUGR), and birth defects22-24.

Reproductive dysfunction caused by atrazine exposure through the hypothalamus-pituitary-gonadal (HPG) axis has been investigated in female rodent models. Developmental and peripubertal studies report a delay in sexual maturation and mammary gland development25. Adult studies report an inhibition of gonadotropin releasing hormone (GnRH) and a reduction in the pre-ovulatory surge of luteinizing hormone (LH), follicle stimulating hormone (FSH), and prolactin (PRL)8,26,27. Furthermore, atrazine has been reported to increase progesterone (P4) levels and is hypothesized to contribute to ovarian degeneration and decreased levels of LH and FSH; potentially leading to early reproductive senescence and dysfunction28,29.

When investigating developmental toxicant exposure and the developmental origins paradigm, the zebrafish provides a strong complementary vertebrate model. There are multiple strengths associated with utilizing the zebrafish including ex utero fertilization and embryonic development, rapid embryogenesis, and a relatively short life span. Paired with these biological strengths are the structural and functional homology of the zebrafish central nervous system (CNS) to humans and the conserved genetic, molecular, and endocrine pathways making the zebrafish a powerful model to assess the later-in-life alterations caused by an embryonic atrazine exposure30,31.

We previously reported that an embryonic atrazine exposure of 0.3, 3, or 30 ppb in zebrafish larvae resulted in immediate alterations to the transcriptome with gene ontology analysis showing enrichment for genes associated with reproductive system function and development, cell cycle regulation, and cancer32. In addition, our previous study examining the DOHaD paradigm showed that an embryonic atrazine exposure alters serotonin turnover and its metabolite 5-hydroxyindoleacetic acid (5-HIAA) in adult female zebrafish with brain transcriptomic profiles indicating enrichment of genes associated with nervous system development and function, behavior, and tissue development33.

In the current study, we aimed to address the later-in-life consequences of an embryonic atrazine exposure by assessing effects and function of the reproductive system, transcriptomic analysis of adult female ovarian tissue, and morphological alterations in the exposed generation's offspring (Supplementary Material, Figure S1). Gonad tissue transcriptomic analysis was compared to the previously completed analysis of brain tissue33 for further assessment of alterations within the HPG axis.

Results

Assessment of an embryonic atrazine exposure on adult zebrafish reproductive function and offspring viability and morphology. We did not observe a skew in sex ratios in any of the treatment groups (p = 0.64; Supplementary Material, Figure S2). The average number of breeding pairs that spawned was significantly lower in the 30 ppb treatment group as compared to other treatment groups (p = 0.008; Fig. 1A), but the average number of embryos per pair (p = 0.21; Fig. 1B) and total number of live embryos in each treatment were not statistically different among treatments (p = 0.08; Fig. 1C). In addition, there were no statistically significant differences in mortality at 24, 48, or 72 hpf (at 24 hpf: p = 0.62; Fig. 1D; 48 and 72 hpf data not shown as there were no additional deaths) or hatching rates among the treatment groups (48 hpf: p = 0.64; Fig. 1E; 72 hpf: p = 0.43; Fig. 1F). Morphological characteristics of the offspring were also measured. While no significant alterations occurred in total body length (p = 0.43; Fig. 2A), a significant decrease in the ratio of head length to total body length in the 30 ppb breeding group (p = 0.0061; Fig. 2B) and a significant increase in the ratio of head width to total body length in the 0.3 and 3 ppb breeding groups were observed (p = 0.0011; Fig. 2C).

Effects of an embryonic atrazine exposure on adult female zebrafish. Approximately 5% of the females from the 30 ppb treatment groups displayed an increase in abdominal swelling (Fig. 3A,B). Two of these individuals had severe swelling to the point of rupture. Acid-fast Ziehl-Neelsen staining revealed absence of mycobacterial organisms in these females indicating abdominal swelling was likely not due to infection (data not shown). Pathological assessment indicated swelling was due to the inability to release eggs. Several endpoints were then assessed to further investigate this observation. No significant differences were observed in the total weight of females in the 30 ppb treatment groups compared to the control treatment group (p = 0.09; Fig. 3C), but there was a significant increase in ovarian weight (p = 0.03; Fig. 3D). There was also no significant difference in GSI (p = 0.11; Supplementary Material, Figure S3). Ten individual females were then analyzed for differences in follicular staging in each of the different treatment groups in each of the four replicates (40 total female fish assessed). The percent follicles in different stages (perinuclear, cortical alveoli, early and late vitellogenin and post-ovulatory) and the percent of atretic follicles did not differ across treatments (Supplementary Material, Figure S4A), but when specifically evaluating females exhibiting swollen abdomens in comparison to those that were not, a significant increase in the number of atretic follicles was observed (p = 0.0002; Fig. 3E,F; Supplementary Material, Figure S4B).

Estradiol and progesterone levels in adult female ovarian tissue. No significant alterations were observed in ovarian tissue concentrations of estradiol in any of the atrazine treatments (p = 0.5337; Fig. 3G).

Figure 1. Assessment of embryonic atrazine exposure on adult zebrafish reproductive function and offspring viability. Adults were individually paired in mating experiments to assess mating success (16 pairs from each of the 4 biological replicates). Average number of pairs that bred was decreased in the group exposed to 30 ppb atrazine during embryogenesis (a). There were no significant differences observed for the number of embryos per pair or the total number of embryos per treatment (b,c, respectively). In addition, no significant changes were observed in mortality of the offspring (d) or in hatching rates at 48 and 72 hpf (e,f, respectively). Error bars are expressed as ± SD. (*p < 0.05).

However, a significant increase in progesterone concentrations were observed in the 3 and 30 ppb atrazine treatment groups (p = 0.0043; Fig. 3H).

Transcriptome analysis of adult female gonad tissue. Transcriptomic profiling of ovaries was assessed to investigate the genetic mechanisms underlying alterations in reproductive success. Results showed expression alterations in 2,024 mapped genes in the 0.3 ppb treatment group, 843 genes in the 3 ppb treatment group, and 696 genes in the 30 ppb treatment group (GSE73740). Of these differentially expressed genes, 383 were common among all three atrazine treatments (Supplementary Material, Table Si). Gene ontology analysis revealed embryonic development, behavior, and organismal survival were enriched for the 0.3 ppb treatment group (Supplementary Material, Table S2). The 3 ppb treatment group showed a similar response as pathways enriched showed changes to genes involved in organ morphology, behavior, and organismal survival (See Supplementary Material, Table S3). Pathways enriched in the 30 ppb treatment group included tissue development, behavior, and

Figure 2. Morphological assessment of offspring. Twenty embryos (considered subsamples; n = 4) of the adult zebrafish population exposed to atrazine during embryogenesis were assessed at 72 hpf to obtain whole larvae total length (measured from snout to end of tail), head length, and head width. No significant alterations were found in total body length (a); however, a significant decrease in the head length to body length ratio in the 30 ppb atrazine treatment group was observed (b). In addition, a significant increase in the head width to body length ratio in the 0.3 and 3 ppb atrazine treatments was observed (c). Error bars are expressed as ± SD. (*p < 0.05).

organismal survival (See Supplementary Material, Table S4). Gene ontology analysis of the 383 genes altered in all three treatment groups revealed effects in endocrine system development and function, tissue development, and behavior (Table 1).

qPCR confirmation of microarray data. A subset of genes detected to be significantly altered by microar-rays in all atrazine treatment groups was independently confirmed by qPCR with all seven target genes (STAR, ACE, CYP19A1, CYP1B1, CRHBP, HPGD, and VIP) significantly changed in the qPCR assessment (Table 2).

Discussion

In order to assess how an embryonic atrazine exposure affects adult reproductive system function we performed paired breeding experiments to assess mating success, spawning, and survival of the progeny. We observed a significant decrease in the number of females that spawned in the 30 ppb treatment group. Several females in

Figure 3. Effects of embryonic atrazine exposure on adult female zebrafish. A representative image of an adult female from the control group (a) compared to an adult female with a swollen abdomen in the 30 ppb developmentally exposed group (b). There was a ~5% incidence rate observed in this treatment group (~24-37 females assessed in each treatment of the four replicates). No significant difference in overall body weight was seen between the control and 30 ppb treatment group (c), but a significant increase in ovary weight was observed (d) (4 replicates with 10 female fish assessed per treatment replicate). Moreover, in comparison to the normal adult female zebrafish (e) an increase in atretic ovarian follicles (black arrow) was observed in those fish that presented with abdominal swelling (f) determined to be from an inability to release eggs (4 replicates with 10 female fish assessed per treatment replicate). Furthermore, while no significant differences were observed in estradiol (g), a significant increase in progesterone (P4) was present in adult females exposed to 3 or 30 ppb atrazine during embryogenesis (h) [4 to 7 adult females were pooled for each biological replicate (n = 4)]. Error bars are expressed as ± SD. (*p < 0.05).

the 30 ppb treatment group presented with swollen abdomens due to an inability to spawn supporting the lack of successful mating with an increase in atretic ovarian follicles present in these individuals. Previous studies have shown that atrazine exposure causes a decrease in spawning events34,35. Tillitt et al.34 exposed adult fathead

Physiological System Development and Function p-valuea Number of Genesb

NERVOUS SYSTEM DEVELOPMENT AND FUNCTION 7.96E-06 - 1.85E-02 106

Growth of neurites 7.96E-06 27

Morphology of CNS 3.35E-05 28

Synaptic Transmission 1.16E-03 17

Differentiation of neurons 3.36E-04 20

ENDOCRINE SYSTEM DEVELOPMENT AND FUNCTION 1.52E-06 - 1.61E-02 22

Hormone Metabolism 6.75E-04 12

Steroidogenesis 6.40E-04 9

Accumulation of progesterone 6.76E-03 2

Concentration of corticosterone 1.23E-05 11

Synthesis of estrogen 6.26E-03 4

BEHAVIOR 2.85E-07 - 1.85E - 02 59

Behavior 2.85E-07 48

Cognition 2.91E-05 23

Learning 7.43E-04 19

Locomotion 7.91E-03 14

TISSUE DEVELOPMENT 7.96E-06 - 1.85E-02 78

Growth of neurites 7.96E-06 27

Gonadogenesis 4.99E-03 20

Development of genital organ 5.51E-03 21

Development of ovary 6.94E-03 8

Table 1. Genes altered in all three atrazine treatments. aDerived from the likelihood of observing the degree of enrichment in a gene set of a given size by chance alone. bClassified as being differentially expressed that relate to the specified function category; a gene may be present in more than one category.

Gene SEQJD ANOVA p-value (p < 0.05) Up/Down

STAR NM_131663.1 0.0134 Up

ACE XM_005169447.1 0.019 Up

CYP19A1 NM_131154.2 0.0312 Up

CYP1B1 NM_001045256.1 0.0015 Up

CRHBP XM_689239.5 0.0023 Up

HPGD NM_001045256.1 0.0027 Up

VIP NM_001114553.2 0.009 Down

Table 2. qPCR array confirmation of adult female gonad tissue.

minnows to 0.5, 5, or 50 ppb atrazine for 14 or 30 days resulting in a decrease in spawning events and cumulative mean egg production with increasing atrazine concentrations34. A second study conducted by the same group exposed adult Japanese medaka to similar atrazine concentrations and confirmed the decrease in mean egg production; however, a decrease in spawning events was not observed35. Interestingly, we observed a similar effect in our current study in adult zebrafish exposed to atrazine only during embryogenesis. Reductions in spawning could be attributed to either a male or female specific effect on reproductive behavior. Although direct breeding behavior was not assessed in our study, previous studies have noted alterations in reproductive behavior in various fish species with atrazine exposure36'37.

Offspring from the successful matings were collected and grown to the end of the embryonic phase to observe mortality, hatching rates, and morphology. While no significant differences were observed for mortality or hatching rate, a decrease in head length-to-body ratio in offspring from the 30 ppb treatment group and an increase in head width-to-body ratio in offspring from the 0.3 and 3 ppb treatment groups was observed. A previous study from our laboratory reported an increase in head length and head-to-body ratio in zebrafish larvae exposed to 0.3, 3, or 30 ppb atrazine through embryogenesis32. In addition, epidemiological studies report a link between in utero atrazine exposure and an increased risk for babies born small for their gestational age (SGA) as well as impairments to fetal growth and birth defects22,24,38. While these studies assess the immediate effects of atrazine exposure during development, the results of the current study indicate that an embryonic atrazine exposure can also stimulate and repress the development of offspring of the exposed generation.

Due to our observed abdominal swelling and decrease in spawning events in the 30 ppb treatment group, gonadal tissue was collected from adult females in each of the treatment groups and subsequent histological

analysis of ovaries was performed to observe frequency of different oocyte stages. Although we found no significant differences in oocyte stage across the treatment groups, an increase in atretic ovarian follicles was observed in females with abdominal swelling. Moreover, the fish with this phenotype tested negative for mycobacteria, which is also reported to elicit swollen abdomens in zebrafish. It should be noted that we did observe females with larger abdomens in the other atrazine treatment groups, albeit not to the extent that was observed in the 30 ppb treatment group. In addition, there were no females that presented with this phenotype in the control groups and we have never observed a similar effect in our breeding groups. A previous study exposed female Sprague-Dawley rats to 3, 30, or 300 mg/kg atrazine for two or four weeks reported a prolonged diestrus period, large-sized atretic follicles, and uterine atrophy in the 300 mg/kg treatment group39. Furthermore, a study conducted by Gojmerac et al.28 reported that an atrazine exposure of 2 mg/kg in feed for 19 days in adult female pigs caused cystic ovarian degeneration and an increase in the persistence of the corpus luteum28. Reproductive hormones play a key role in maintaining reproductive homeostasis. Of primary importance is LH which plays a key role in follicular development and normal reproductive function26,40. Studies show that atrazine reduces the release of LH from the pituitary gland26'41. Therefore, the results observed here and in previously mentioned studies suggest an anovulatory effect of atrazine which could be due to the previously observed reduction in LH, providing a mechanism behind the observed reduction in spawning, atretic follicles, and cystic ovarian degeneration39.

In addition, this study is innovative in that we are the first to show an increase in ovarian progesterone (P4) in adult female zebrafish exposed to 3 or 30 ppb atrazine only during embryogenesis. A study conducted by Gojmerac et al.28 reported that an atrazine exposure of 2 mg/kg in feed for 19 days in adult female pigs caused a significant increase in P428. As previously discussed, results from this study also included an increase in cystic ovarian degeneration. It was shown that higher levels of P4 inhibit follicle development and can potentiate follicular atresia29. Additional in vivo studies report that atrazine exposure in Long-Evans Hooded and Sprague Dawley rats also caused significant increases in P442-44. As previously noted, atrazine exposure decreased LH and FSH leading to reproductive dysfunction8,26,45. Furthermore, increases in P4 may contribute to the inhibition of the pre-ovulatory surge of LH from the pituitary enhancing reproductive dysfunction. Furthermore, in vitro studies utilizing swine and rat granulosa cells have reported increases in P4 following atrazine exposure46,47. An in vitro study conducted by Pogrmic-Majkic et al.47 provided mechanistic data regarding the observed increase in P4 elicited by atrazine exposure. Results suggest that atrazine causes an increase in P4 and overexpression of luteal markers (STAR and CYPllAl) through the stimulation of the cAMP, AKT, and CEBPB-signaling pathways47. The rise in P4 (a critical step in early luteinization) was also demonstrated in H295R adrenal cortical carcinoma cells48. The findings that atrazine elevates P4 combined with the observed follicular atresia indicates the ability of atrazine to alter ovarian function and reduce fertility following only a developmental exposure. Results from our study also revealed no significant alterations in ovarian levels of 17ß-estradiol. Current literature surrounding the effects of atrazine on estrogen is primarily conducted in male animal models; while levels of LH, FSH, PRL, P4, and gonadotropin releasing hormone (GnRH) were the forefront in female models8,27,44,45.

Transcriptomic analysis was also completed on the ovarian tissue of adult female zebrafish exposed to atrazine during embryogenesis in order to identify genetic alterations behind the observed decrease in spawning and increase in P4. Analysis revealed a subset of genes altered that are involved in endocrine system development and function. The first key gene of interest is steroidogenic acute regulatory protein (STAR). This protein is responsible for the transport of cholesterol from the outer to inner mitochondrial membrane where it is then metabolized to pregnenolone49. A study conducted by Pogrmic et al.49 showed a decrease in STAR and other steroidogenic genes (scavenger receptor B1 (SR-B1), steroidogenic factor 1 (SF-1), cytochrome P450 17A1 (CYP17A1), and 17ß-hydroxysteroid dehydrogenase (17ß-HSD)) in Leydig cells after exposure from PND 23-51 to 50 or 200 mg/ kg/day atrazine49. A second study conducted by Pogrmic-Majkic et al.50 only treated Leydig cells for 24 hours with 20 |iM atrazine which caused an increase in expression of STAR, SF-1, and 17ß-HSDso. The results from these studies led to the hypothesis of a potential transient stimulatory effect of atrazine. This up regulation of STAR was also observed in FSH stimulated granulosa cells after atrazine exposure47. Results from these previous studies indicate that atrazine can elicit alterations in STAR. Our results provide new insight into the long-term effects of atrazine on female reproductive function due to the increase in STAR expression.

One of the most controversial genes associated with atrazine exposure is aromatase (CYP19A1). Aromatase is responsible for catalyzing the conversion of testosterone to estradiol as well as androstenedione to estrone in the steroidogenic pathway51. Transcriptomic results from our study reveal an increase in CYP19A1 expression in all three atrazine treatments approximately six months after exposure completion. Although an increase in CYP19A1 was observed in ovarian tissue, no alterations were observed in levels of gonadal estradiol as is hypothesized52,53. However, further investigation is warranted to determine if the increase in CYP19A1 elicits a response in plasma levels of estradiol and a corresponding decrease in testosterone. In 2002 and 2010, Hayes et al. reported that atrazine exposure resulted in hermaphroditic, demasculinized, and chemically castrated X. laevis at concentrations as low as 0.1 and 2.5 ppb (|ig/L)52,53. Hayes et al.53 hypothesized that the observed morphological alterations were due to an increase in CYP19A1 expression which would therefore lead to an increase in estrogen and a decrease in testosterone53. This hypothesis has been under investigation as other studies have not demonstrated that atrazine causes increases in CYP19A1 expression in X. laevis and other aquatic species54-58. In vitro studies report contrasting results indicating that atrazine exposure alters aromatase mRNA and activity in human cancer cell lines48,59-61 and luteinized granulosa cells51.

An additional gene altered in ovarian tissue was the neuropeptide adenylate cyclase-activating peptide 1 (ADCYAP1) which regulates gonadotropin gene expression. ADCYAP1 has widespread distribution and function which includes expression in the central and peripheral nervous system, adrenal glands, placenta, ovaries, and testes62. Estradiol and P4 stimulate the expression of ADCYAP1 in the hypothalamus and ovary63-66. The observed increase in ADCYAP1 expression in ovary tissue corresponds to the observed increase in P4 levels. Additionally, our previous study examining the effects of a developmental atrazine exposure on female brain

Atrazine exposure during embryogenesis (1-72 hpf)

•I©

■Sol

Figure 4. Schematic diagram of the effects of an embryonic atrazine exposure on the HPG axis in female zebrafish. Diagram representing the effects of an embryonic atrazine exposure at 72 hpf as well as in brain and ovarian tissue of adult females exposed during embryogenesis compiling the data from this study, Weber et al.32, and Wirbisky et al.33. The genes in green are representative of changes observed at 72 hpf; genes in purple are altered in adult female brain tissue; genes listed in pink are those altered in adult female gonad tissue; and red signifies the morphological alterations observed.

tissue also revealed an alteration in ADCY upon transcriptomic analysis33 as it works synergistically with GnRH, therefore playing a role in the release of LH and FSH62. All three are needed for proper reproductive function and each is shown to be decreased by atrazine exposure8,27,40. Furthermore, as previously stated, atrazine increases P4 through the stimulation of the cAMP, AKT, and CEBPB-signaling pathways47. CREB signaling proteins are also identified as a stimulator of ADCYAP167.

An increase was also observed in CYP1B1 in all three atrazine treatments. CYP1B1 is expressed constitutively in steroidogenic tissues such as the adrenal, ovary, and testes68. Therefore, its regulation is primarily by hormones that elevate cAMP. Steroidogenic cytochrome P450's undergo transcriptional regulation by numerous transcription factors such as SF-1. SF-1 is a transcriptional regulator for CYP19A1, STAR, and CYP11A1. Although this gene is heavily integrated within the steroidogenic pathway, its up regulation may provide further insight into the mechanisms behind the adverse effects of atrazine exposure, especially those effects lasting into adulthood.

An additional gene altered by all three atrazine treatments was vasoactive intestinal peptide (VIP). VIP mRNA is found in the ovaries of rodents, specifically, in the granulosa layer of the pre-ovulatory follicles. VIP acts on specific targets which activates adenylyl cyclase pathway, therefore stimulating the cAMP production. This gene stimulates progesterone, androgens, and estradiol secretion. It is also associated with the initiation of follicular growth and the inhibition of follicular atresia69. VIP was down regulated in our transcriptomic analysis; this down regulation could play a role in the observed follicular atresia observed. The observed upregulation of ADCYAP1 may be a compensatory mechanism to the down regulation of VIP to continue the release of steroid hormones.

Additionally, it is known that EDCs can affect many aspects of the hormonal pathway including their synthesis, transport, excretion, and metabolism7. Our microarray and qPCR analysis show that embryonic atrazine exposure elicits an up-regulation of the CRHBP gene. This gene is a potent stimulator of the synthesis and secretion of preopiomelanocortin-derived peptides such as a-Melanocyte-stimulating hormone (a-MSH) and adren-ocorticotropic hormone (ACTH). Although these additional hormones are highly associated with the body's primary stress pathway (hypothalamus-pituitary-adrenal (HPA) axis); this axis works in tandem with the HPG axis. Gene ontology revealed CRHBP to aid in the cellular response to GnRH. GnRH is a key regulator of the HPG axis and has been the focus of studies aiming to elucidate the mechanism of atrazine toxicity27,40.

Our data suggests a developmental origin of reproductive dysfunction in adult female zebrafish caused by an embryonic atrazine exposure. Reproductive traits, morphological and hormone data, and gene expression analysis provide a conceptual working model of atrazine toxicity when coupled to our previous studies32,33 (Fig. 4). Indeed, many adult diseases and disorders are believed to originate from exposure to environmental stressors during development by altering and reprogramming normal cellular processes that lead to an altered physiological state. The disruption of the reproductive system through the developing neuroendocrine system was identified in our initial gene expression assessment32. From our transcriptomic analysis, we reported a disruption in the expression of CYP17A1, LH, and ADCY1. These genes are critical for steroidogenesis and proper reproductive function. Additionally, alterations in PER1 and PER3 were observed which play a role in regulating circadian rhythm; a process vital for maintaining homeostatic reproductive function70. We also reported alterations in THRA; a gene responsible for proper thyroid function at 72 hpf and in adult female brain tissue following the developmental atrazine exposure32,33. Although the primary effect of atrazine exposure is reported throughout the HPG axis, results from our developmental and adult studies show that the hypothalamus-pituitary-thyroid (HPT) axis may also be disrupted, although further studies are necessary. In addition, we observed alterations in PDE10A following embryonic atrazine exposure as well as PDE1A in adult brain tissue. Phosphodiesterase genes are reported to be affected by atrazine exposure71,72. Alterations in these upstream cellular regulators can have an

adverse effect on steroidogenesis through cAMP. Transcriptomic profiling of the brain tissue of females exposed to atrazine only during embryogenesis also revealed enrichment in genes associated with neurological system development and function, organ and tissue development, body size, and behavior as well as a decrease in the serotonin metabolite 5-HIAA and serotonin turnover33. Although neurotransmitter systems were not addressed in the results shown here, neurotransmitters and neuropeptides that are involved in GnRH release are reported to be altered by atrazine exposure73-75. Furthermore, transcriptomic profiling of ovarian tissue of females exposed to atrazine during embryogenesis revealed enrichment of genes associated with endocrine system development and function and steroidogenesis (CYP19A1, CYP1B1, ADCYAP2, STAR, VIP) providing a genetic link to the observed reproductive dysfunction, follicular atresia, and increased levels of progesterone. The results from these and previous studies provide support for atrazine exposure causing adverse effects to the neuroendocrine system immediately following atrazine exposure as well as later in life.

Methods

Zebrafish husbandry and experimental design. Zebrafish (wild-type AB strain) were housed in a Z-Mod System (Aquatic Habitats) on a 14:10 hour light:dark cycle and maintained at 28 °C with a pH of 6.9-7.2 and conductivity range of 470-520 |iS. Adult zebrafish were bred in cages and embryos were collected, staged, and rinsed with system water as described previously for experimental use76. Embryos were dosed with 0, 0.3, 3, or 30 ppb atrazine (CAS #1912-24-9; Chem Service, 98% purity) from 1-72 hours post fertilization (hpf) as previously described32,33. Atrazine sample concentrations were verified using an U.S. EPA approved immunoassay kit for atrazine (Abraxis Atrazine ELISA kit, Warminster, PA) as previously described77. After exposure, larvae were rinsed with clean fish system water, housed in 4-liter tanks in the system with 20-30 fish per tank, and allowed to mature under normal growing conditions (Supplementary Material, Figure S1). Biological replicates are defined by groups started on different days. All animal protocols were approved and performed in accordance with Purdue University's Institutional Animal Care and Use Committee guidelines.

Assessment of reproductive potential and morphology of offspring. At approximately 5 months post fertilization (mpf) zebrafish from each treatment group (0, 0.3, 3, or 30 ppb atrazine) were sexed and bred within respective treatment groups in 4-liter tanks once a week for three continuous weeks in order to initiate breeding cycles (Supplementary Material, Figure S1). To assess individual breeding potential, one male and one female from each treatment group were randomly selected and paired for breeding. A total of 64 individual pairs were tested for each treatment group among the four biological replicates (16 pairs per replicate). The number of pairs that spawned was noted and embryos from each pair were collected and observed through embryogenesis (72 hpf) noting: the number of fertilized embryos and unfertilized eggs per pair; mortality at 24, 48 and 72 hpf; and number of embryos hatched at 48 and 72 hpf. In addition, twenty larvae from each treatment dish (considered subsamples) of each biological replicate (n = 4) were analyzed with light microscopy using a Nikon SMZ1500 dissecting scope with NIS Elements imaging software to attain total body length (measured from snout to end of tail), head length, and head width.

Fish gonadal tissue collection and histology. Adult female zebrafish (approximately 5-8 mpf) were euthanized in MS-222 (Ethyl 3-aminobenzoate methanesulfonate, Sigma, St. Louis, MO) (4 mg/mL) and subsequently weighed to obtain body weight (4 replicates with 10 fish per replicate). Ovaries were collected and weighed for calculation of gonadal somatic index (GSI) determined as (gonad weight/total tissue weight)*100. A sample of ovarian tissue was collected from both ovaries and fixed in Davidson's Fixative78 overnight at room temperature and then transferred to histology grade 70% EtOH (Sigma). Hematoxylin and eosin (H & E) sections were prepared following standard procedures and sections examined under a light microscope (10-40X) for staging and evaluation of any abnormalities as outlined in Johnson et al.78. All follicles were counted and classified as perinuclear, cortical alveoli, early vitellogenic, late vitellogenic, post-ovulatory, or atretic. Frequencies of occurrence of each stage were then calculated from each section and averaged.

Estradiol and progesterone analysis in adult female ovary tissue. Ovarian tissue from adult females (~6 mpf) was collected from each treatment group. Four to seven adult females were pooled for each biological replicate (n = 4). Tissue was stored at — 80 °C prior to extraction and analysis. Samples were weighed and 100 mg of tissue placed in a 1.7 mL tube. 1000 ng of d5-estradiol or 100 ng of d9-progesterone was added as an internal standard to each tube. Glass homogenizer beads and water were added, tubes briefly vortexed, and then 1 mL of hexane/ethyl acetate (60/40 v/v) added to each. The samples were vortexed for 15 minutes on maximum speed and then centrifuged at 15,000 rpm for 5 minutes. The top organic layer was transferred to a new tube, dried in a rotary evaporation device, and stored at —20 °C until derivatization.

For estradiol analysis the extracts were derivatized with dansyl chloride just prior to LC/MS/MS analysis79. Acetonitrile, 10mM sodium carbonate, and freshly prepared dansyl chloride was added to each sample for a reaction time of 10 minutes at 60 °C. For progesterone analysis the extracts were derivatized with the AB Sciex Keto derivatization kit (AB Sciex, Framingham, MA) just prior to LC/MS/MS analysis. 50 |iL of reagent was added to each sample with a reaction time of 60 minutes at room temperature. The samples were transferred to autosampler vials and immediately analyzed on the instrument.

An Agilent 1200 Rapid Resolution liquid chromatography (LC) system coupled to an Agilent 6460 series QQQ mass spectrometer (MS) was used to analyze estradiol in each sample. A Waters Xbridge C18 2.1 mm x 100 mm, 3 |im column was used for LC separation. The buffers were (A) water + 0.1% formic acid and (B) acetonitrile + 0.1% formic acid. The linear LC gradient was as follows: time 0 minutes, 10% B; time 1 minute, 10% B; time 5 minutes, 100% B; time 15 minutes, 100% B; time 15.5 minutes, 10% B; time 18 minutes, 10% B. The flow rate was 0.3 mL/min. Multiple reaction monitoring was used for MS analysis. The data were acquired

in positive electrospray ionization (ESI) mode by monitoring the following transitions: Estradiol-dansyl chloride 506.1^155.8 (40V), 171 (30 V); d5-Estradiol-dansyl chloride 511.1^155.8 (40V), 171 (30 V); Progesterone 429.1^370 (20V), 126 (30V); d9-Progesterone 438.1^379 (20V), 132 (30V). The jet stream ESI interface had a gas temperature of 325 °C, gas flow rate of 8 L/minute, nebulizer pressure of 40 psi, sheath gas temperature of 250 °C, sheath gas flow rate of 7 L/minute, capillary voltage of 4000 V, and nozzle voltage of 1000 V.

Transcriptome microarray analysis of adult female gonad tissue. Ovarian tissue from six adult females (approximately 5-8 mpf) was collected among the groups exposed to an atrazine treatment (0, 0.3, 3, or 30 ppb) during embryogenesis, homogenized in Trizol (Life Technologies), and flash frozen in liquid nitrogen (n = 6). Total RNA was isolated by the RNeasy Mini Kit (Qiagen). Transcriptomic microarray analysis was conducted using the one-color hybridization strategy to compare gene expression profiles among the atrazine treatments with a custom zebrafish 4 x 180K expression platform (Agilent Technologies). This microarray is a multiplex format of 4 arrays each consisting of 180,000 probes interrogating 36,000 known and predicted targets with approximately 3-5 probes per target and is based on the Ensembl and UCSC Genome Databases. Following hybridization, arrays were washed and then scanned on an Agilent Technologies SureScan Microarray Scanner (Agilent Technologies). Array image data was extracted using Agilent Feature Extraction Software 11.5 (Agilent Technologies,). Data was uploaded to GeneSpring 12.5 (Agilent Technologies) for statistical analysis. Microarray analysis was performed following MIAME guidelines80. Each gene list was imported into Ingenuity Pathway Analysis (IPA) for gene ontology and molecular pathway analysis. Genes referred to in the results and discussion sections are reported as the human homologs of the genes identified to be altered by microarrays.

Quantitative polymerase chain reaction (qPCR) confirmation of microarray. qPCR was performed on a subset of selected genes altered in the gene expression array analysis: STAR, ACE, CYP19A1, CYP1B1, CRHBP, HPGD, and VIP using the BioRad SSOAdvance SYBR Green Supermix kit according to the manufacturer's recommendations. Probes specific to target genes were designed using the Primer3 website (Supplementary Material, Table S5). qPCR was performed following similar methods as previously described32,33 following the minimum information for publication of quantitative real-time PCR experiments (MIQE) guidelines81. Similar to as performed in previous studies in our laboratory several genes were assessed to determine the best reference gene to be used for this data set (data not shown)32,33. ¡3-ACTIN was found to be the most consistent and least variable for this analysis. qPCR was performed on the same samples as used in the microarray analysis (n = 6). Experimental samples were run in triplicate (technical replicates) and gene expression was normalized to (3-ACTIN. Efficiency and specificity were checked with melting and dilution curve analysis and no-template controls.

Statistical analysis. The average number of embryos per treatment group, number of embryos per pair, survival, hatching success, spawning rates, sex ratios, morphology measurements of offspring, and gene expression array confirmation by qPCR were analyzed using an Analysis of Variance (ANOVA) with a least significant difference (LSD) post-hoc test when a significant ANOVA was observed (a = 0.05) using SAS statistical software. Total body weight and ovary weights between the control and 30 ppb treatment groups was analyzed with a T-test (p < 0.05). Counts of ovarian follicles and corresponding stages of development were assessed using a combination of Chi-square and a generalized linear mixed model (p < 0.05) in SAS statistical software. Hormone analysis was performed using an ANOVA and Tukey's post-hoc test when a significant ANOVA was observed (a = 0.05). Microarray analysis was completed with an ANOVA and a Tukey's post-hoc test when a significant ANOVA was observed (a = 0.05). In addition, a mean absolute log2 expression ratio of at least 0.585 (50% increase or decrease in expression) must be satisfied.

References

1. Barouki, R., Gluckman, P. D., Grandjean, P., Hanson, M. & Heindel, J. J. Developmental origins of non-communicable disease: implications for research and public health. Environ Health. 11, 1-9 (2012).

2. Zhang, X. & Ho, S. M. Epigenetics meets endocrinology. J. Mol Endocrinol. 46, 11-32 (2011).

3. Gluckman, P. D., Hanson, M. A., Cooper, C. & Thornburg, K. L. Effect of in utero and early-life conditions on adult health and disease. N Engl J Med. 359, 61-73 (2008).

4. Jirtle, R. L. & Skinner, M. K. Environmental epigenomics and disease susceptibility. Nat Rev Genet. 8, 253-262 (2007).

5. Birnbaum, L. S. & Fenton, S. E. Cancer and developmental exposure to endocrine disruptors. Environ Health Perspect. 111, 389-394 (2003).

6. Ma, Y., Zhang, P., Wang, F., Yang, J., Yang, Z. & Qin, H. The relationship between early embryo development and tumourigenesis. J. Cell Mol Med. 14, 2697-2701 (2010).

7. Vandenberg, L. N. Low-dose effects of hormones and endocrine disruptors. Vitam Horm. 94, 129-165 (2014).

8. Cooper, R. L., Stoker, T. E., Tyrey, L., Goldman, J. M. & McElroy, W. K. Atrazine disrupts the hypothalamic control of pituitary-ovarian function. Toxicol Sci. 53, 297-307 (2000).

9. Hatch, E. E., Troisi, R., Wise, L. A., Titus-Ernstoff, L., Hyer, M., Palmer, J. R., Strohsnitter, W. C., Robboy, S. J., Anderson, D., Kaufman, R., Adam, E. & Hoover, R. N. Preterm birth, fetal growth, and age at menarche among women exposed prenatally to diethylstilbestrol (DES). Reprod Toxicol. 31, 151-157 (2011).

10. Roy, J. R., Chakraborty, S. & Chakraborty, T. R. Estrogen-like endocrine disrupting chemicals affecting puberty in humans-a review. Med SciMonit. 15, 137-145 (2009).

11. Swedenborg, E., Ruegg, J., Makela, S. & Pongratz, I. Endocrine disruptive chemicals: mechanisms of action and involvement in metabolic disorders. J. Mol Endocrinol. 43, 1-10 (2009).

12. Wetzel, L. T., Luempert, L. G. 3rd., Breckenridge, C. B., Tisdel, M. O., Stevens, J. T., Thakur, A. K., Extrom, P. J. & Eldridge, J. C. Chronic effects of atrazine on estrus and mammary tumor formation in female Sprague-Dawley and Fischer 344 rats. J. Toxicol Environ Health. 43, 169-182 (1994).

13. Guerrero-Bosagna, C., Settles, M., Lucker, B. & Skinner, M. K. Epigenetic transgenerational actions of vinclozolin on promoter regions of the sperm epigenome. PLoS One 5, Pii: e13100 (2010).

14. Solomon, K. R., Carr, J. A., Du Preez, L. H., Giesy, J. P., Kendall, R. J., Smith, E. E. & Van Der Kraak, G. J. Effects of atrazine on fish, amphibians, and aquatic reptiles: a critical review. Crit Rev Toxicol. 38, 721-772 (2008).

15. Freeman, L. E., Rusiecki, J. A., Hoppin, J. A., Lubin, J. H., Koutros, S., Andreotti, G., Zahm, S. H., Hines, C. J., Coble, J. B., Barone-Adesi, F., Sloan, J., Sandler, D. P., Blair, A. & Alavanja, M. C. Atrazine and cancer incidence among pesticide applicators in the Agricultural Health Study (1994-2007). Environ Health Perspect. 119, 1253-1259 (2011).

16. Gammon, D. W., Aldous, C. N., Carr, W. C. Jr., Sanborn, J. R. & Pfeifer, K. F. A risk assessment of atrazine use in California: human health and ecological aspects. PestManag Sci. 61, 331-355 (2005).

17. Rinsky, J. L., Hopenhayn, C., Golla, V., Browning, S. & Bush, H. M. Atrazine exposure in public drinking water and preterm birth. Public Health Rep. 127, 72-80 (2012).

18. Rohr, J. R. & McCoy, K. A. A qualitative meta-analysis reveals consistent effects of atrazine on freshwater fish and amphibians. Environ Health Perspect. 118, 20-32 (2010).

19. U.S. EPA. 2002. List of Contaminants and Their MCLs. EPA 816-F-02-013. Washington, DC: U.S. Environmental Protection Agency.

20. Commission E. Review report for the active substance atrazine; Finalized in the Standing Committee on the Food Chain and Animal Health at its meeting on 3 October 2003 in support of a decision concerning the non-inclusion of atrazine in Annex I of Directive 91/414/EEC and the withdrawal of authorization for plant protection products containing this active substance:European Commission Health and Consumer Protection Direc-torate-General; 2003. SANC0/10496/2003-final.

21. Sass, J. B. & Colangelo, A. European Union bans atrazine, while the United States negotiates its continued use. Int. J. Occup. Environ. Health. 12, 260-267 (2006).

22. Munger, R., Isacson, P., Hu, S., Burns, T., Hanson, J., Lynch, C. F., Cherryholms, K., Van Drope, P. & Hausler, W. J. Intrauterine growth retardation in Iowa communities with herbicide-contaminated drinking water supplies. Environ Health Perspect. 105, 308-314 (1997).

23. Ochoa-Acuna, H., Frankenberger, J., Hahn, L. & Carbajo, C. Drinking-water herbicide exposure in Indiana and Prevalence of small-for-gestational-age and preterm delivery. Environ. Health. Perspect. 117, 1619-1624 (2009).

24. Winchester, P. D., Huskins, J. & Ying, J. Agrichemicals in surface water and birth defects in the United States. Acta Paediatrica. 98, 664-669 (2009).

25. Rayner, J. L. Enoch, R. R. & Fenton, S. E. Adverse effects of prenatal exposure to atrazine during a critical period of mammary gland growth. Toxicol. Sci. 87, 255-266 (2005).

26. Foradori, C. D., Hinds, L. R., Hanneman, W. H., Legare, M. E., Clay, C. M. & Handa, R. J. Atrazine inhibits pulsatile luteinizing hormone release without altering pituitary sensitivity to a gonadotropin-releasing hormone receptor agonist in female Wistar rats. Biol Reprod. 81, 40-45 (2009).

27. Foradori, C. D., Zimmerman, A. D., Hinds, L. R., Zuloaga, K. L., Breckenridge, C. B. & Handa, R. J. Atrazine inhibits pulsatile gonadotropin-releasing hormone (GnRH) release without altering GnRH messenger RNA or protein levels in the female rat. Biol Reprod. 88, 1-7 (2013).

28. Gojmerac, T., Kartal, B., Curie, S., Zurie, M., Kusevic, S. & Cventi, Z. Serum biochemical changes associated with cystic ovarian degeneration in pigs after atrazine treatment. Toxicol Lett. 85, 9-15 (1996).

29. Telfer, E., Gosden, R. G. & Faddy, M. J. Impact of exogenous progesterone on ovarian follicular dynamics and function in mice. J. Reprod Fertil. 93, 263-269 (1991).

30. de Esch, C., Slieker, R., Wolterbeek, A., Woutersen, R. & De Groot, D. Zebrafish as potential model for developmental neurotoxicity testing: A mini review. Neurotoxicol Teratol. 34, 545-553 (2012).

31. Howe, K. et al. The zebrafish reference genome sequence and its relationship to the human genome. Nature. 496, 498-503 (2013).

32. Weber, G. J., Sepulveda, M. S., Peterson, S. M., Lewis, S. L. & Freeman, J. L. Transcriptome alterations following developmental atrazine exposure in zebrafish are associated with disruption of neuroendocrine and reproductive system function, cell cycle, and carcinogenesis. Toxicol Sci. 132, 458-466 (2013).

33. Wirbisky, S. E., Weber, G. J., Sepulveda, M. S., Xiao, C., Cannon, J. R. & Freeman, J. L. Developmental origins of neurotransmitter and transcriptome alterations in adult female zebrafish exposed to atrazine during embryogenesis. Toxicology. 333, 156-167 (2015).

34. Tillitt, D. E., Papoulias, D. M., Whyte, J. J. & Richter, C. A. Atrazine reduces reproduction in fathead minnow (Pimephales promelas). Aquat Toxicol. 99, 149-159 (2010).

35. Papoulais, D. M., Tillittm, D. E., Talykinam, M. G., Whyte, J. J. & Richter, C. A. Atrazine reduces reproduction in Japanese medaka. Aquatic Toxicol. 154, 230-239 (2014).

36. Moore, A. & Lower, N. The impact of two pesticides on olfactory-mediated endocrine function in mature male Atlantic salmon (Salmo salar L.) parr. Comp Biochem Physiol B Biochem Mol Bio. 129, 269-276 (2001).

37. Shenoy, K. Prenatal exposure to low doses of atrazine affects mating behaviors in male guppies. Horm Behav. 66, 439-448 (2014).

38. Chevrier, C., Limon, G., Monfort, C., Rouget, F., Garlantezec, R., Petit, C., Durand, G. & Cordier, S. Urinary biomarkers of prenatal atrazine exposure and adverse birth outcomes in the PELAGIE birth cohort. Environ Health Perspect. 119, 1034-1041 (2011).

39. Shibayama, H., Kotera, T., Shinoda, Y., Hanada, T., Kajihara, T., Ueda, M., Tamura, H., Ishibashi, S., Yamashita, Y. & Ochi, S. Collaborative work on evaluation of ovarian toxicity. Two- or four-week repeated-dose studies and fertility study of atrazine in female rats. Toxicol Sci. 34, SP147-155 (2009).

40. Foradori, C. D., Hinds, L. R., Hanneman, W. H. & Handa, R. J. Effects of atrazine and its withdrawal on gonadotropin-releasing hormone neuroendocrine function in the adult female Wistar rat. Biol Reprod. 81, 1099-1105 (2009).

41. Foradori, C. D., Coder, P. S., Tisdel, M., Yi, K. D., Simpkins, J. W., Handa, R. J. & Breckenridge, C. B. The effect of atrazine administered by gavage or in diet on the LH surge and reproductive performance in intact female Sprague-Dawley and Long-Evans rats. Birth Defects Res B Dev Reprod Toxicol. 101, 262-275 (2014).

42. Fraites, M. J., Cooper, R. L., Buckalew, A., Jayaraman, S., Mills, L. & Laws, S. C. Characterization of the hypothalamic-pituitary-adrenal axis response to atrazine and metabolites in the female rat. Toxicol Sci. 112, 88-99 (2009).

43. Goldman, J. M., Davis, L. K., Murr, A. S. & Cooper, R. L. Atrazine-induced elevation or attenuation of the LH surge in the ovariectomized, estrogen-primed female rat: role of adrenal progesterone. Reproduction. 146, 305-314 (2013).

44. Taketa, Y., Yoshida, M., Inoue, K., Takahashi, M., Sakamoto, Y., Watanabe, G., Taya, K., Yamate, J. & Nishikawa, A. Differential stimulation pathways of progesterone secretion from newly formed corpora lutea in rats treated with ethylene glycol monomethyl ether, sulpiride, or atrazine. Toxicol Sci. 121, 267-287 (2011).

45. Foradori, C. D., Hinds, L. R., Quihuis, A. M., Lacagnina, A. F., Breckenridge, C. B. & Handa, R. J. The differential effect of atrazine on luteinizing hormone release in adrenalectomized adult female Wistar rats. Biol Reprod. 85, 684-689 (2011).

46. Basini, G., Bianchi, F., Bussolati, S., Baioni, L., Ramoni, R., Grolli, S., Conti, V., Bianchi, F. & Grasselli, F. Atrazine disrupts steroidogenesis, VEGF and NO production in swine granulosa cells. Ecotoxicol Environ Saf. 85, 59-63 (2012).

47. Pogrmic-Majkic, K., Samardzija, D., Fa, S., Hrubik, J., Glisic, B., Kaisarevic, S. & Andric, N. Atrazine enhances progesterone production through activation of multiple signaling pathways in FSH-stimulated rat granulosa cells: evidence for premature luteinization. Biol Reprod. 91, 1-10 (2014).

48. Tinfo, N. S., Hotchkiss, M. G., Buckalew, A. R., Zorrilla, L. M., Cooper, R. L. & Laws, S. C. Understanding the effects of atrazine on steroidogenesis in rat granulosa and H295R adrenal cortical carcinoma cells. Reprod Toxicol. 31, 184-193 (2011).

49. Pogrmic, K., Fa, S., Dakic, V., Kaisarevic, S. & Kovacevic, R. Atrazine oral exposure of peripubertal male rats down regulates steroidogenesis gene expression in Leydig cells. Toxicol Sci. 111, 189-197 (2009).

50. Pogrmic-Majkic, K., Fa, S., Dakic, V., Kaisarevic, S. & Kovacevic, R. Upregulation of peripubertal rat Leydig cell steroidogenesis following 24 h in vitro and in vivo exposure to atrazine. Toxicol. Sci. 118, 52-60 (2010).

51. Holloway, A. C., Anger, D. A., Crankshaw, D. J., Wu, M. & Foster, W. G. Atrazine-induced changes in aromatase activity in estrogen sensitive target tissues. J. Appl Toxicol. 28, 260-270 (2008).

52. Hayes, T. B., Collins, A., Lee, M., Mendoza, M., Noriega, N., Stuart, A. A. & Vonk, A. Hermaphroditic, demasculinized frogs after exposure to the herbicide atrazine at low ecologically relevant doses. Proc Natl Acad Sci USA 99, 5476-5480 (2002).

53. Hayes, T. B., Khoury, V., Narayan, A., Nazir, M., Park, A., Brown, T., Adame, L., Chan, E., Buchholz, D., Stueve, T. & Gallipeau, S. Atrazine induces complete feminization and chemical castration in male African clawed frogs (Xenopus laevis). Proc Natl Acad Sci USA 107, 4612-4617 (2010).

54. Coady, K. K., Murphy, M. B., Villeneuve, D. L., Hecker, M., Jones, P. D., Carr, J. A., Solomon, K. R., Smith, E. E., Van Der Kraak, G., Kendall, R. J. & Giesy, J. P. Effects of atrazine on metamorphosis, growth, laryngeal and gonadal development, aromatase activity, and sex steroid concentrations in Xenopus laevis. Ecotoxicol Environ Saf. 62, 160-173 (2005).

55. Kazeto, Y., Place, A. R. & Trant, J. M. Effects of endocrine disrupting chemicals on the expression of CYP19 genes in zebrafish (Danio rerio) juveniles. Aquat Toxicol. 69, 25-34 (2004).

56. Hecker, M., Park, J. W., Murphy, M. B., Jones, P. D., Solomon, K. R., Van Der Kraak, G., Carr, J. A., Smith, E. E., du Preez, L., Kendall, R. J. & Giesy, J. P. Effects of atrazine on CYP19 gene expression and aromatase activity in testes and on plasma sex steroid concentrations of male African clawed frogs (Xenopus laevis). Toxicol Sci. 86, 273-280 (2005).

57. Murphy, M. B., Hecker, M., Coady, K. K., Tompsett, A. R., Higley, E. B., Jones, P. D., Du Preez, L. H., Solomon, K. R., Carr, J. A., Smith, E. E., Kendall, R. J., Van Der Kraak, G. & Giesy, J. P. Plasma steroid hormone concentrations, aromatase activities and GSI in ranid frogs collected from agricultural and non-agricultural sites in Michigan (USA). Aquat Toxicol. 77, 153-166 (2006).

58. Thibaut, R. & Porte, C. Effects of endocrine disrupters on sex steroid synthesis and metabolism pathways in fish. J. Steriod Biochem Mol Biol. 92, 485-494 (2004).

59. Sanderson, J. T., Seinen, W., Giesy, J. P. & Van dan Berg, M. 2-Chloro-s-triazine herbicides induce aromatase (CYP19) activity in H295R human adrenocortical carcinoma cells: a novel mechanism for estrogenicity? Toxicol Sci. 1, 121-127 (2000).

60. Sanderson, J. T., Letcher, R. J., Heneweer, M., Giesy, J. P. & Van den Berg, M. Effects of chloro-s-triazine herbicides and metabolites on aromatase activity in various human cell lines and on vitellogenin production in male carp hepatocytes. Environ Health Perspect.

109, 1027-1031 (2001).

61. Heneweer, M., Van den Berg, M. & Sanderson, J. T. A comparison of human H295R and rat R2C cell line as in vitro screening tools for effects on aromatase. Toxicol Lett. 142, 183-194 (2004).

62. Thomas, R. L., Crawford, N. M., Grafer, C. M., Zheng, W. & Halvorson, L. M. GATA augments GnRH-mediated increases in Adcyapl gene expression in pituitary gonadotrope cells. J. Mol Endocrinol. 51, 313-324 (2013).

63. Apostolakis, E. M., Lanz, R. & O'Malley, B. W. Pituitary adenylate cyclase-activating peptide: a pivotal modulator of steroid-induced reproductive behavior in female rodents. Mol Endocrinol. 18, 173-183 (2004).

64. Ha, C. M., Kang, J. H., Choi, E. J., Kim, M. S., Park, J. W., Kim, Y., Choi, W. S., Chun, S. Y., Kwon, H. B. & Lee, B. J. Progesterone increases mRNA levels of pituitary adenylate cyclase-activating polypeptide (PACAP) and type I PACAP receptor (PAC(1)) in the rat hypothalamus. Brain Res Mol Brain Res. 78, 59-68 (2000).

65. Moore, J. P. Jr., Burger, L. L., Dalkin, A. C. & Winters, S. J. Pituitary adenylate cyclase activating polypeptide messenger RNA in the paraventricular nucleus and anterior pituitary during the rat estrous cycle. Biol Reprod. 73, 491-499 (2005).

66. Park, J. I., Kim, W. J., Wang, L., Park, H. J., Lee, J., Park, J. H., Kwon, H. B., Tsafriri, A. & Chun, S. Y. Involvement of progesterone in gonadotrophin-induced pituitary adenylate cyclase-activating polypeptide gene expression in pre-ovulatory follicles of rat ovary. Mol Hum Reprod. 6, 238-245 (2000).

67. Grafer, C. M., Thomas, R., Lambrakos, L., Montoya, I., White, S. & Halvorson, L. M. GnRH stimulates expression of PACAP in the pituitary gonadotropes via both PKA and PKC signaling systems. Mol Endocrinol. 23, 1022-1032 (2009).

68. Zheng, W. & Jefcoate, C. R. Steroidogenic factor-1 interacts with cAMP response element binding protein to mediate cAMP stimulation of CYP1B1 via an upstream enhancer. Mol Pharmacol. 67, 499-512 (2005).

69. Rosas, G., Ramírez, M. I., Linares, R., Trujillo, A., Domínguez, R. & Morales-Ledesma, L. Asymmetric steroidogenic response by the ovaries to the vasoactive intestinal peptide. Endocrine. 48, 968-977 (2015).

70. Boden, M. J. & Kennaway, D. J. Circadian rhythms and reproduction. Reproduction 132, 379-392 (2006).

71. Kucka, M., Pogrmic-Majkic, K., Fa, S., Stojilkovic, S. S. & Kovacevic, R. Atrazine acts as an endocrine disrupter by inhibiting cAMP-specific phosphodiesterase-4. Toxicol Appl Pharmacol. 265, 19-26 (2012).

72. Roberge, M., Hakk, H. & Larsen, G. Atrazine is a competitive inhibitor of phosphodiesterase but does not affect the estrogen receptor. Toxicol Lett. 154, 61-68 (2004).

73. Das, P. C., McElroy, W. K. & Cooper, R. L. Differential modulation of catecholamines by chlorotriazine herbicides in pheochromocytoma (PC12) cells in vitro. Toxicol Sci. 56, 324-331 (2000).

74. Das, P. C., McElroy, W. K. & Cooper, R. L. Potential mechanisms responsible for chlorotriazine-induced alteratiosn in catecholamines in pheochromocytoma (PC12) cells. Life Sci. 73, 3123-3138 (2003).

75. Shafer, T. J., Ward, T. R., Meacham, C. A. & Cooper, R. L. Effects of the chlorotriazine herbicide, cyanazine, on GABA(A) receptors in cortical tissue from rat brain. Toxicology. 142, 57-68 (1999).

76. Peterson, S. M., Zhang, J., Weber, G. J. & Freeman, J. L. Global gene expression analysis reveals dynamic and developmental stage-dependent enrichment of lead-induced neurological gene alterations. Environ Health Perspect. 119, 615-621 (2011).

77. Freeman, J. L., Beccue, N. & Rayburn, A. L. Differential metamorphosis alters the endocrine response in anuran larvae exposed to T3 and atrazine. Aquat. Toxicol. 75, 263-276 (2005).

78. Johnson, R. J., Wolf, J. & Braunbeck, T. OECD guidance document for the diagnosis of endocrine-related histopathology of fish gonads. 96 pp (2009).

79. Zhang, F., Bartles, M. J., Geter, D. R., Carr, M. S., McClaymount, L. E., Marino, T. A. & Klecka, G. M. Simultaneous quantitation of testosterone, estradiol, ethinyl estradiol, and 11-ketotestosterone in fathead minnow fish plasma by liquid chromatography/positive atmospheric pressure photoionization tandem mass spectrometry. Rapid CommunMass Spectrom. 23, 3637-3646 (2009).

80. Brazma, A., Hingamp, P., Quackenbush, J., Sherlock, G., Spellman, P., Stoeckert, C., Aach, J., Ansorge, W., Ball, C. A., Causton, H. C., Gaasterland, T., Glenisson, P., Holstege, F. C., Kim, I. F., Markowitz, V., Matese, J. C., Parkinson, H., Robinson, A., Sarkans, U., Schulze-Kremer, S., Stewart, J., Taylor, R., Vilo, J. & Vingron, M. Minimum information about microarray experiment (MIAME)-toward standards for microarray data. Nat. Genet. 29, 365-371 (2001).

81. Bustin, S. A., Benes, V., Garson, J. A., Hellemans, J., Huggett, J., Kubista, M., Mueller, R., Nolan, T., Pfaffl, M. W., Shipley, G. L., Vandesompele, J. & Wittwer, C. T. The MIQE Guidelines: Minimum information for publication of quantitative real-time PCR experiments. Clin Chem. 55, 611-622.

Acknowledgements

This work was supported by the National Institutes of Health, National Institute of Environmental Health Sciences [R15 ES019137 to J.L.F. and M.S.S.], a Purdue Center for Cancer Research Innovative Research Pilot Project [J.L.F. and M.S.S.], a Purdue Research Foundation Grant [J.L.F. and G.J.W.], and a Bilsland Dissertation Fellowship [S.E.W. and J.L.F].

Author Contributions

S.E.W. wrote the manuscript, designed and performed microarray and qPCR experiments, and collected samples for hormone analysis; G.J.W. designed and performed breeding experiments and reviewed the manuscript; M.S.S. performed histological assessment and reviewed manuscript; T.L.L. performed pathology assessment of female ovarian tissue and reviewed the manuscript; A.S.J. performed the LC/MS for hormone analysis and reviewed the manuscript; J.L.F. was involved in the conception, experimental design, interpretation of the data, and manuscript preparation.

Additional Information

Accession Codes: Transcriptomic data discussed in this publication has been deposited in NCBI's Gene Expression Omnibus and is accessible through GEO Series accession number GSE73740. (http://www.ncbi.nlm. nih.gov/geo/query/acc.cgi?acc=GSE73740).

Supplementary information accompanies this paper at http://www.nature.com/srep

Competing financial interests: The authors declare no competing financial interests.

How to cite this article: Wirbisky, S. E. et al. An embryonic atrazine exposure results in reproductive dysfunction in adult zebrafish and morphological alterations in their offspring. Sci. Rep. 6, 21337; doi: 10.1038/ srep21337 (2016).

l/JJv ® I This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/