EB and GEN were the only two compounds that significantly affected both the timing of vaginal opening and the ability to maintain a regular estrous cycle. GnRH activation was strongly impaired in these animals, as well as the PPT animals, suggesting that the mechanism of disruption is ERα dependent. Because GnRH neurons do not express ERα, this observation supports the hypothesis that endocrine disruption occurs within ERα expressing neurons that send efferent projections to GnRH neurons, such as KISS neurons. KISS fiber density was significantly lower in both the AVPV and ARC of the EB treated animals but only in the AVPV of the GEN treated animals. The AVPV population of KISS neurons is thought to regulate steroid positive feedback in mammals, including humans (
Navarro et al., 2004;
Pompolo et al., 2006;
Roa et al., 2006;
Shibata et al., 2007;
Smith et al., 2006b;
Tena-Sempere, 2006). Therefore, collectively the observations in the present study suggest that impaired organization of the steroid positive feedback circuitry on GnRH neurons might be the primary mechanism underling the observed changes in reproductive physiology following neonatal exposure to EDCs.
Prior studies have repeatedly demonstrated that developmental exposure to estrogenic compounds, including phytoestrogens and other EDCs, can alter the timing of puberty, disrupt ovarian cycle function, and masculinize the female hypothalamus (
Colborn et al., 1993;
Jefferson et al., 2006b;
Lephart et al., 2005;
McLachlan et al., 2006;
Patisaul, 2005;
Patisaul and Polston, 2007). To our knowledge, this is the first time that similar effects for the isoflavone metabolite EQ have been reported. Advanced puberty following exposure to GEN, or soy protein containing GEN and other isoflavones, has been observed previously (
Badger et al., 2001). For example, Leaf et al, found that exposure to four days of 10 mg/kg bw per day of GEN, beginning on gestational day 15, (via administration by sc injection to the pregnant dam) also resulted in accelerated vaginal opening and altered ovarian cycle frequency (
Levy et al., 1995). A similar outcome was reported for rats exposed to 40 mg/kg bw GEN per day from birth to weaning (
Lewis et al., 2003) but not 4 mg/kg bw GEN per day. This may indicate that the dose of 10 mg/kg bw used in the present study is approaching the minimum level for which an effect of GEN on these endpoints is observed. In the present study, we found no significant effect of EQ on the timing of vaginal opening, but it remains to be seen whether or not there is a dose effect of EQ. One study failed to find advanced pubertal onset in the rat following oral administration of GEN across the first five days of life in doses ranging from 12.5 mg/kg bw to 100 mg/kg bw (
Nagao et al., 2001), which may indicate that the route of exposure may also affect outcome. In the present study we opted to administer our compounds by sc injection because a thorough literature search indicated that this was the exposure route most frequently used in other studies using PPT and DPN. Therefore we administered all of the test compounds by sc injection to be consistent and ensure effective delivery of PPT and DPN. A newly published study has also now demonstrated that subcutaneous injection of GEN is a suitable model for oral exposure indicating that route of exposure may not be a serious concern (
Jefferson et al., 2007).
Neonatal administration of GEN or EQ also significantly abrogated GnRH neuronal activation following hormone priming in adulthood. This observation is consistent with a prior study which showed that GEN has a defeminizing effect on GnRH secretion (
Faber and Hughes, 1993). The specific mechanism by which this occurs is unclear although activity on one or both ER subtype has long been hypothesized to be crucial. ERα and ERβ are differentially distributed throughout the hypothalamus (
Li et al., 1997;
Shughrue et al., 1996;
Shughrue et al., 1997) suggesting that they play distinct roles in GnRH activation and signaling. Therefore, selective disruption through one ER subtype may produce disparate outcomes on GnRH function. Phytoestrogens generally have a higher relative binding affinity for ERβ than ERα
in vitro (
Kuiper et al., 1998;
Mueller et al., 2004). However,
in vivo they have been shown to alter both ERα and ERβ dependent gene expression in the brain (
Jacob et al., 2001;
Patisaul et al., 2002;
Patisaul et al., 1999) indicating that they at least have the potential to act through either subtype. The degree to which GnRH activation was suppressed by the phytoestrogens was more similar to the results seen with DPN than with PPT suggesting that the phytoestrogens may be operating through an ERβ mediated mechanism. Alternatively, this observation may also indicate that the phytoestrogens are activating as weak agonists on ERα. This distinction is significant because GnRH neurons express ERβ but not ERα throughout development (
Herbison and Pape, 2001;
Hrabovszky et al., 2000;
Hrabovszky et al., 2001). Therefore, activity through ERβ would suggest that phytoestrogens can act on GnRH neurons directly. Further work is necessary to resolve this critical issue.
In the present study, the ERβ agonist DPN produced a less severe outcome than PPT on the ability to maintain a regular estrous cycle and the capacity to activate GnRH neurons in response to hormone priming. The observation that both PPT and DPN suppressed ovarian cycle function is consistent with effects previously reported following post-natal exposure to two other ER specific agonists (
Patchev et al., 2004). In this prior study, agonism of ERβ was found to be more potent on estrous cyclicity, an effect opposite to what was observed here. In the present study, although impairment of GnRH activation by ERβ agonism did not reach statistical significance, it is possible that, within the DPN treated group, GnRH activation was lowest in the individual animals that entered persistent diestrus and higher in the animals that did not. This would then support the hypothesis that EDCs can directly interact with GnRH through ERβ. Unfortunately the sample sizes in the present study were too small to investigate this possibility but future studies will address this issue. The present data suggest that ERα plays a more significant mechanistic role because while the DPN effect did not quite reach statistical significance, ovarian cyclicity and GnRH activation were completely absent in the PPT group. Further studies using a range of DPN and PPT doses are needed to determine if this simply reflects a dose effect and that higher doses of DPN could also eliminate ovarian cyclicity and GnRH activation.
Since GnRH neurons do not express ERα it has been hypothesized that estrogen signaling is conveyed to GnRH neurons by ERα-expressing neurons elsewhere in the hypothalamus, particularly the AVPV (
Simerly, 1998;
Simerly, 2002). In the female AVPV, estrogen administration markedly increases KiSS-1 mRNA expression and the number of KISS-ir cell bodies (
Irwig et al., 2004;
Navarro et al., 2004;
Smith et al., 2005) but this effect can be eliminated by neonatal androgen treatment (
Kauffman et al., 2007b). Here we have shown that neonatal exposure to EB or GEN had a similar effect on AVPV KISS fiber density. AVPV KISS neurons express ERα (
Smith et al., 2005;
Smith et al., 2006b) so presumably ERα plays a mechanistic role in this effect but unfortunately sample sizes were not sufficient to ascertain if DPN or PPT had an effect on KISS fiber density in the AVPV. In the adult female, it has recently been demonstrated that neuronal ERα is required for the generation of the LH surge (
Roa et al., 2008;
Wintermantel et al., 2006) but the relative roles ERα and ERβ play in the organization of the neural circuitry required to generate a gonadotropin surge in response to hormone priming are less well understood. The differential distribution of ERα and ERβ in the neonatal rat AVPV have not yet been definitively described making it difficult to conclude that ERα is playing the most substantial role in mediating the organizing effects of estrogens. For example, a pair of studies conducted by the same research group (
Orikasa et al., 2002;
Orikasa and Sakuma, 2003) have clearly identified both ERβ-ir and ERβ mRNA in the AVPV of male and female rats on PND 7 and beyond (day of birth defined as PND 1 in these studies) but did not look earlier. In contrast a study conducted by a different research group observed no ERβ-ir within the AVPV in the first 2 weeks of life (
Perez et al., 2003). Therefore the neonatal distribution of ERα and ERβ in the rat AVPV remains undefined.
We also found that neonatal exposure to EB reduced KISS fiber density in the ARC following hormone priming in adulthood. ERα appears to play a mechanistic role in this process as KISS fiber density was also reduced in the ARC of the PPT but not the DPN treated animals, an effect which is not unexpected given that ERα is abundant in the ARC while ERβ is absent (
Shughrue and Merchenthaler, 2001;
Simerly et al., 1990). The significant effect of EB was somewhat surprising given that the number of KiSS-1 neurons nor the quantity of KiSS-1 expression in the ARC has not been found to be sexually dimorphic in rats (
Kauffman et al., 2007b). Prior studies have shown that KiSS-1 expression and the number of cell bodies immunoreactive for kisspeptin in the ARC are significantly decreased by steroid hormone administration in adults (
Kauffman et al., 2007b). Very few KISS-ir cell bodies were present in any of our ARC sections, regardless of treatment, an observation that is consistent with these prior studies. Without an untreated control group, it is difficult to ascertain whether the absolute number of KISS neurons in the ARC, or the AVPV, is reduced by EDC treatment or if their response to hormone priming in adulthood is modified. Experiments designed to distinguish between these two alternative hypotheses are underway. It will also be critical to determine if the KISS fibers quantified in the present study make direct contact with OVLT GnRH neurons.
Finally, it is important to note that it is possible that GEN may produce its effects via mechanisms that are not dependent on ERs. GEN is a potent tyrosine kinase inhibitor so some of the observed effects may be attributable, at least in part, to this mechanism of action. GEN has also been found to alter aromatase production and steroid hormone secretion
in vivo (
Akingbemi et al., 2004;
Brooks and Thompson, 2005;
Kishida et al., 2001;
Pelissero et al., 1996;
Whitehead et al., 2002). Although our findings suggest that ERα plays a mechanistic role, particularly in the suppression of the estrous cycle and GnRH activation, further biochemical and pharmacological experiments are required to conclusively determine the specific mechanisms by which GEN, EQ and other EDCs alter HPG organization. It is also essential to acknowledge that premature vaginal opening and estrous cycle loss could result from disruption anywhere within the HPG axis including the ovary. Disrupted ovarian morphology, including absence of corpora lutea, multi-oocyte follicles and inhibited oocyte nest breakdown (
Jefferson et al., 2006a;
Jefferson et al., 2002;
Jefferson et al., 2006b;
Kouki et al., 2003;
Nagao et al., 2001) have been reported in both rats and mice following developmental exposure to GEN. By ovariectomizing and hormone replacing our animals prior to sacrifice we aimed to eliminate this potentially confounding effect of ovarian dysfunction.
Unlike rodents, the human hypothalamus is differentiated during gestation, therefore the critical period akin to the one explored in the present study occurs before birth in humans. The GnRH pulse generator is functional by the end of the first trimester (
Grumbach, 2002), but it is currently unknown when or how KISS neurons differentiate across gestation and early development in humans. GEN, EQ and most other phytoestrogens readily cross the placenta (
Todaka et al., 2005) as do many other EDCs suggesting that human hypothalamic development may be vulnerable to disruption by these compounds however human exposure levels may be too low to produce appreciable effects. Identifying the most sensitive and physiologically relevant critical periods as well as the most salient dose ranges of exposure is essential for environmental risk assessment. Our data support the hypothesis that the time during which the hypothalamus undergoes sexual differentiation is likely one of these critical periods. These findings emphasize the need for further studies aimed at uncovering the mechanisms by which exposure to EDCs, during pre- and postnatal development ultimately impact the health of exposed offspring.