Studies have linked sex-biased neurodevelopmental disorders, including autism and schizophrenia, with fetal antecedents such as prenatal stress. Further, these outcomes can persist into subsequent generations, raising the possibility that aspects of heritability in these diseases involve epigenetic mechanisms. Utilizing a mouse model in which we previously identified a period in early gestation when stress results in dysmasculinized and stress-sensitive male offspring, we have examined programming effects in second-generation offspring of prenatally stressed (F2-S) or control (F2-C) sires. Examination of gene expression patterns during the perinatal sensitive period, when organizational gonadal hormones establish the sexually dimorphic brain, confirmed dysmasculinization in F2-S males, where genes important in neurodevelopment showed a female-like pattern. Analyses of the epigenomic miRNA environment detected significant reductions in miR-322, miR-574, and miR-873 in the F2-S male brain, levels that were again more similar to those of control females. Increased expression of a common gene target for these three miRNAs, β-glycan, was confirmed in these males. These developmental effects were associated with the transmission of a stress-sensitive phenotype and shortened anogenital distance in adult F2-S males. As confirmation that the miRNA environment is responsive to organizational testosterone, neonatal males administered the aromatase inhibitor formestane exhibited dramatic changes in brain miRNA patterns, suggesting that miRNAs may serve a previously unappreciated role in organizing the sexually dimorphic brain. Overall, these data support the existence of a sensitive period of early gestation when epigenetic programming of the male germline can occur, permitting transmission of specific phenotypes into subsequent generations.
Epidemiological studies have linked prenatal stress to increases in the incidence of neurodevelopmental disorders, including schizophrenia and autism spectrum disorders, associations that are often sex dependent (Huttunen and Niskanen, 1978; van Os and Selten, 1998; Khashan et al., 2008; Kinney et al., 2008). These diseases often display sex differences in prevalence, presentation, or therapeutic outcomes (Bale et al., 2010). While many factors likely contribute to these differences, sex-specific responses to fetal antecedents are likely involved (Weinstock, 2007).
We have previously identified early gestation as a specific window of sensitivity during which male mice were susceptible to the programming effects of maternal stress. These males exhibited physiological and behavioral stress sensitivity and cognitive deficits, endophenotypes associated with human neuropsychiatric disease. In addition, these changes reduced or disrupted established sex differences by dysmasculinizing male offspring measures of stress responsivity (Mueller and Bale, 2007, 2008). Similar disruptions of sex differences in behavior, morphology, and gene expression profiles have previously been reported in studies using prenatal stress paradigms across multiple species (Ward, 1972; Meisel et al., 1979; Reznikov et al., 1999; Kapoor and Matthews, 2005; Biala et al., 2010). The organizational/activational hypothesis of brain development suggests that a surge of gonadal hormones organize the brain in a sexually dimorphic manner during the perinatal sensitive period. Then in adulthood, gonadal hormones can activate this organized neurocircuitry to express appropriate sex-specific behavioral phenotypes, including stress axis responsivity (Phoenix et al., 1959; Arnold and Gorski, 1984; Seale et al., 2005; Bingham and Viau, 2008). The disruption of sex differences identified in our model suggests that early prenatal stress alters the trajectory of neurodevelopment during the perinatal period.
Fetal antecedents likely contribute to adult disease through programming changes in the epigenome. Examples of this phenomenon are emerging in human studies. For example, infants with prenatal exposure to maternal depression or anxious mood exhibited increased glucocorticoid methylation, which was associated with a heightened cortisol response to a mild stressor (Oberlander et al., 2008). Such programming effects may transmit to subsequent generations, predisposing offspring to disease. Animal models have clearly established a role for epigenetics in transgenerational phenotypic inheritance following exposure to environmental factors such as maternal stress, diet, and endocrine disruptors (Dunn and Bale, 2009; Franklin et al., 2010; Guerrero-Bosagna et al., 2010; Skinner, 2011). Importantly, in rodent models, inheritance of a phenotype through the paternal lineage excludes confounding effects of the maternal intrauterine environment and postnatal rearing behavior, and indicates gametic epigenetic transmission (Youngson and Whitelaw, 2008). These phenomena raise the interesting possibility that the heritability observed in neurodevelopmental diseases could be the product of both classic genetic and non-Mendelian, or epigenetic, mechanisms (Maher, 2008; Slatkin, 2009). Therefore, we used our mouse model to determine epigenetic programming effects and the transmission of a dysmasculinized phenotype to second-generation (F2) male offspring via the paternal lineage, eventually focusing mechanistically on analysis of the miRNA environment as a form of transcriptional regulation with broad potential to impact developmental processes.
Materials and Methods
All dams bred for gestational stress studies were virgin, experimentally naive C57BL/6:129 F1 hybrid 5-week-old mice purchased from The Jackson Laboratory. Offspring from these breedings were used to generate the second-generation (F2) litters for transgenerational studies. Justification for using a hybrid background strain in these studies is related to stress responsivity phenotypes and physiology. C57BL/6 are extremely low-stress responders and display low levels of maternal care, making them poor choices for studies focusing on neurodevelopment. While 129 mice are great stress responders and show high levels of quality maternal care, they frequently lack a fully formed corpus callosum and are poor performers in behavioral tests, especially learning and memory tasks. However, the combination of these two strains produces a hybrid vigor that has served our research well with predictable stress responses, behavioral outcomes, and sex differences in stress physiology and behavioral tests (Mueller and Bale, 2006, 2007, 2008). Pregnancy was established by confirmation of a copulation plug (checked for each morning between 7:00 and 8:00 A.M.). Presence of a copulation plug denoted experimental day 1 for early prenatal stress exposure. The pregnant female was individually housed, given a cotton nestlet, and randomly assigned to a stress treatment or control group. Food (Purina Rodent Chow; 28.1% protein, 59.8% carbohydrate, 12.1% fat) and water was provided ad libitum throughout the study. All studies were performed according to experimental protocols approved by the University of Pennsylvania Institutional Animal Care and Use Committee, and all procedures were conducted in accordance with institutional guidelines.
Early prenatal stress
Administration of chronic variable stress was performed as described previously (Mueller and Bale, 2006). Briefly, pregnant dams were randomly assigned to either an experimental treatment group to receive chronic variable stress during gestation days 1–7, or to a control nonstressed treatment group. Pregnant dams assigned to the stress group experienced a different stressor on each of the seven days. Stressors included: 36 h constant light, 15 min of fox odor exposure (1:10,000 2,4,5-trimethylthiazole, Acros Organics), novel objects in cage overnight (eight marbles of similar size and color), 5 min restraint stress in a 50 ml conical tube, novel white noise overnight (Sleep machine, Brookstone), multiple cage changes, and water saturated bedding overnight. These mild stressors were selected to be nonhabituating and to not induce pain. We previously demonstrated that this chronic variable stress paradigm does not affect maternal food or water intake, weight gain, gestation length, litter size, or postpartum maternal behaviors (Mueller and Bale, 2006).
To determine whether the dysmasculinized stress-sensitive phenotype identified in first-generation (F1) prenatally stressed males could be transmitted through the paternal lineage to second-generation (F2) offspring, F1 control or prenatally stressed males were bred with F1 control females to generate F2 control litters (F2-C n = 9 litters, n = 69 total animals) or F2 stress litters (F2-S n = 9 litters, n = 59 total animals). Molecular, physiological, morphological, and behavioral correlates of the dysmasculinized stress-sensitive phenotype were then tested in these F2 offspring. All experimental n values described throughout the manuscript refer to litter n values.
Postnatal day 1 brain expression analysis
To explore mechanisms through which the F1 dysmasculinized stress-sensitive phenotype may be programmed in the F2 generation, we examined whole-brain gene expression in male and female neonates at birth, a period critical for the organization of the sexually dimorphic brain by gonadal hormones. One male and one female neonate per litter were killed on the day of parturition. Whole brains were dissected, frozen in liquid nitrogen, and stored at −80°C before assay. Neonate sex was confirmed by SRY genotyping of tail genomic DNA as previously described (Mueller and Bale, 2008). Whole brains were sonicated in TRIzol reagent (Invitrogen) and total RNA was isolated according to manufacturer's protocol.
Custom Taqman qRT-PCR Array.
Two hundred fifty nanograms of total RNA from postnatal day 1 (PN1) brains was reverse transcribed to cDNA using the High-Capacity cDNA reverse transcriptase kit (Applied Biosystems). Expression levels of 93 genes important in neurodevelopment were determined using a custom Taqman array (F2-C♂ n = 4, F2-S♂ n = 3, and F2-C♀ n = 4) (Applied Biosystems). Analysis was performed using the comparative Ct method. Eighteen svedberg rRNA was used as an endogenous loading control. Expression levels of each sample were normalized to the average F2-C♂ expression level.
miRNA Taqman qRT-PCR Array.
Five hundred nanograms of total RNA from PN1 brain was reverse transcribed to cDNA using Megaplex RT pool A primers and Multiscribe reverse transcriptase (Applied Biosystems). Expression levels of 239 miRNAs were determined using the Taqman Array MicroRNA card A Array (F2-C♂ n = 5, F2-S♂ n = 3, and F2-C♀ n = 5) (Applied Biosystems). Analysis was performed using the comparative Ct method. For each sample, the average of the Ct values of sno135 and sno202 was used as an endogenous loading control. Expression levels of each sample were normalized to the average F2-C♂ expression level.
qRT-PCR of predicted miRNA target transcripts.
Predicted miRNA targets were identified using the web-based algorithm miRDB (Wang, 2008; Wang and El Naqa, 2008). Expression of mRNAs that were shared predicted targets of miRNAs identified as significantly changed by early prenatal stress, and with potential relevance to disease mechanisms, were assayed by Taqman qRT-PCR (Applied Biosystems). GAPDH was used as an endogenous loading control. Samples were run in triplicate, and expression was determined using the comparative Ct method. Expression levels of each sample were normalized to the average F2-C♂ expression level.
F2 physiological and behavioral measures
For each test, one female and one male per litter were examined to avoid litter effects. An investigator blind to animal treatment group and sex performed scoring of behavioral tests.
Tail suspension test.
The tail suspension test (TST) was performed as previously described to examine active versus passive stress coping strategies (Steru et al., 1985). Mice were secured to a rod by adhesive tape placed ∼1 cm from the tip of the tail, and suspended 50 cm from the bench-top in a visually isolated area (F2-C♂ n = 7, F2-S♂ n = 7, F2-C♀ n = 8, and F2-S♀ n = 6). Immobility time, defined as the absence of all movement except for whisker movement or respiration, was scored for the 6 min test session.
As a test with predictable sex differences in outcome measures of learning and memory, we conducted the Barnes maze (F2-C♂ n = 7, F2-S♂ n = 7, F2-C♀ n = 9, and F2-C♀ n = 9) as previously described (Mueller and Bale, 2007). Briefly, the maze consists of a black circular disk (90 cm in diameter) with 24 holes evenly spaced around its perimeter. An escape box is located under one of the holes. The location of the escape box remains constant throughout training. The disk is elevated 70 cm above the floor and situated in a room with white walls containing three distinct visual cues around the perimeter of the maze. All mice were trained on the maze 2 trials/d for 3 d. Trials within each day were separated by 4 h. To begin each trial, the mouse was placed under a glass beaker in the center of the maze for 15 s before trial start. Latency to identify the target escape box was determined. Each trial was terminated when the mouse located and entered the target escape box or after 4 min elapsed. If the mouse did not successfully locate the target box, the investigator guided the mouse to the target and a latency to target of 240 s was assigned.
HPA response to restraint stress.
The HPA axis response to acute stress was measured by exposing adult F2 mice to a 15 min restraint in a 50 ml conical tube (F2-C♂ n = 7, F2-S♂ n = 7, F2-C♀ n = 7, and F2-C♀ n = 9). Testing occurred 2–5 h following lights on. Blood samples were collected from a tail nick at four time points: (1) time 0, immediately upon removal from the cage, (2) time 15, immediately after the restraint stress, (3) time 30, following 15 min recovery in the home cage, and (4) time 90, following 75 min recovery in the home cage. Samples were collected into EDTA-treated tubes, centrifuged, and plasma was stored at −80°C until corticosterone levels were measured by radioimmunoassay (MP Biomedicals). The minimum detection limit of the assay was 7.7 ng/ml, and the intra-assay coefficient of variation was 7.1%.
One experimentally naive male and female per litter were killed as adults (F2-C♂ n = 7, F2-S♂ n = 7, and F2-C♀ n = 9). Anogenital distances and left testis weights were obtained in males.
Effects of an aromatase inhibitor on the neonatal brain miRNA environment
To examine the influence of organizational estradiol on the neonatal brain miRNA environment, a separate cohort of C57BL/6:129 control pups were treated with an aromatase inhibitor or vehicle on the morning following parturition. Male pups were randomly assigned to receive 20 μg of formestane (Sigma-Aldrich) in 20 μl of sesame oil with 10% ethanol (n = 7 from 6 litters) or vehicle injections (n = 7 from 6 litters). This dose, after adapting for differences in rat versus mouse neonate weight, was previously shown to reduce male hypothalamic estrogen to female levels (Amateau et al., 2004). All female pups received vehicle injections (n = 8 from 5 litters). Injections were administered subcutaneously between the shoulders, and the injection site was treated with New Skin liquid bandage to prevent leakage.
miRNA expression analysis.
Pups were killed 24 h after treatment (PN2). Whole brains were dissected, frozen in liquid nitrogen, and stored at −80°C before assay. Neonate sex was confirmed by SRY genotyping of tail genomic DNA as previously described (Mueller and Bale, 2008). Whole brains were sonicated in TRIzol reagent (Invitrogen) and total RNA was isolated according to the manufacturer's protocol.
Five hundred nanograms of total RNA from PN2 brain was reverse transcribed to cDNA using Megaplex RT pool A primers and Multiscribe reverse transcriptase (Applied Biosystems). Expression levels of 239 miRNAs were determined using the Taqman Array MicroRNA card A Array (Applied Biosystems). Analysis was performed using the comparative Ct method. For each sample, the average of the Ct values of sno135 and sno202 was used as an endogenous loading control. Expression levels of each sample were normalized to the average expression level of vehicle-treated males.
An investigator blind to animal treatment group and sex conducted all studies and analyses. To control for litter effects in prenatal manipulations, only one male or female from a litter was used for a test or assay. We tested the null hypothesis that early prenatal stress exposed males could not transmit their dysmasculinized stress-sensitive phenotype to their second-generation offspring. Morphological measures and immobility in the TST were analyzed by one-tailed Student's t test, testing for effects identified in early prenatal stress exposed (F1) sires. Barnes maze performance and corticosterone response to restraint was analyzed by one-way ANOVA for paternal (F1) prenatal stress exposure, and trial or time as repeated measure, respectively. miRNA and gene target expression in F2-C♂, F2-S♂, and F2-C♀ PN1 whole brains was analyzed by one-way ANOVA, using Student's t tests for post hoc analysis. miRNA expression data in control or formestane-treated brains were analyzed by hierarchical clustering of samples, using Pearson correlation as a metric. Treatment effects on individual miRNAs were identified by two-tailed Student's t tests. Heat maps and hierarchical clustering was performed using MultiExperiment Viewer (TM4.org). All other statistics were performed using JMP8 (SAS).
F2 gene expression analysis during the perinatal sensitive period
F2 PN1 brain gene expression analysis
To explore mechanisms through which the F1 dysmasculinized stress-sensitive phenotype may be programmed in the F2 generation, we examined brain gene expression in male and female neonates at birth, a period critical for the organization of the sexually dimorphic brain by gonadal hormones. Expression data for individual genes are presented in Table 1. These data are also displayed in a heat map to highlight the broad similarities in gene expression that exist between F2-S♂ and F2-C♀ relative to F2-C♂ (Fig. 1A). Mean F2-S♂ expression of the 17 genes that display significant basal sex differences in PN1 whole brain are plotted on a continuum between average F2-C♂ and F2-C♀ expression (Fig. 1B). In F2-S♂, 13 of these 17 genes displayed expression levels closer to F2-C♀ than to F2-C♂ levels.
To identify potential mediators of the program of dysmasculinized gene expression, we assayed the F2 PN1 brain expression of aromatase, ERα, and ERβ, known effectors of masculinization during the perinatal sensitive period (Fig. 1C). While there was no significant effect of group on aromatase (F(2,10) = 0.69; p = 0.70) or ERα levels (F(2,10) = 1.51; p = 0.27), ERβ expression was elevated in F2-S♂ compared to F2-C♂ (F(2,10) = 4.23; p = 0.05).
F2 postnatal day 1 (PN1) brain expression of miRNA and predicted targets
To identify potential alternative mediators of a program of dysmasculinized gene expression in F2-S males, the F2 PN1 brain miRNA environment was assayed. miRNAs with expression levels that displayed a statistically significant effect of group are displayed in Table 2. The expression of three of these miRNAs (miR-322, miR-574–3p, and miR-873) appeared dysmasculinized in F2-S♂ (Fig. 1D).
To determine whether miRNA changes in F2 PN1 brains were associated with altered expression of target transcripts, expression of genes that were shared predicted targets of miRNAs identified as significantly changed by early prenatal stress, and with potential relevance to disease mechanisms, were assayed by Taqman qRT-PCR. β-Glycan was the only predicted target of all three dysmasculinized miRNAs (miR-322, miR-574–3p, and miR-873). There was significantly greater expression of β-glycan in F2-S♂ than in F2-C♂ (F(2,10) = 4.99; p = 0.03) (Fig. 1D). There was also significantly greater expression of Reep3, the shared predicted target of miR-302b and miR-28, in F2-S♂ (1.23 ± 0.11) compared to F2-C♀ (0.91 ± 0.06) (F(2,10) = 4.64; p = 0.04). There were no significant differences in expression of the additional predicted targets of miR-322 and miR-873, Plxna2 (F(2,10) = 0.72; p = 0.51) and Prkar2a (F(2,10) = 0.23; p = 0.23), or of the predicted targets of miR-28 and miR-302b, Unk (F(2,10) = 0.37; p = 0.70) and Hif1an (F(2,10) = 0.68; p = 0.53).
Analysis of adult F2 behavior and physiology
F2 adult male morphology
To assess the degree to which paternal (F1) prenatal stress exposure affects morphological measures of masculinization, we examined male anogenital distances and testis weights (Fig. 2A,B). F2-S♂ had reduced anogenital distances (one-tail t = 2.7; p = 0.01) and reduced testis weights (one-tail t = 1.97; p = 0.04).
F2 adult behavior
To determine whether prenatal stressed males (F1) could transmit their dysmasculinized stress-sensitive phenotype to F2 offspring, we examined F2 adult performance in the TST and Barnes maze. F2-S♂ spent significantly more time immobile than F2-C♂ (one-tailed t(12) = 1.85, p = 0.04) (Fig. 2C). There was no corresponding increase in F2-S♀ immobility relative to F2-C♀ (one-tailed t(12) = 0.17, p = 0.44) (Fig. 2D). Analyzing Barnes maze performance, there was no statistically significant between-subjects effect of paternal (F1) prenatal stress exposure in males (F(1,9) = 0.09; p = 0.40). As expected, there was a significant within subjects effect of time (F(5,5) = 6.91; p = 0.03) (Fig. 2E). In females, there was also no significant between-subjects effect of paternal (F1) prenatal stress exposure (F(1,11) = 0.03; p = 0.57), and again there was a significant within-subjects effect of time (F(5,7) = 4.34; p = 0.02) (Fig. 2F).
F2 adult HPA stress axis
To examine the impact of paternal (F1) prenatal stress exposure on F2 offspring HPA axis sensitivity, we examined corticosterone levels in response to a 15 min restraint stress. In males, there was no significant between-subjects effect of paternal prenatal stress exposure on corticosterone levels (F(1,12) = 0.96; p = 0.35), but there was a significant within-subjects effect of time (F(3,10) = 51.0; p < 0.0001) (Fig. 2G). In females, as in males, there was no significant between-subjects effect of paternal prenatal stress exposure on corticosterone response (F(1,14) = 0.89; p = 0.36), though there was a significant within-subjects effect of time (F(3,12) = 33.5; p < 0.0001) (Fig. 2H).
Effects of formestane treatment on the neonatal miRNA environment
To determine the role of organizational gonadal hormones in the regulation of the neonatal brain miRNA environment, we assayed changes in brain miRNA complement 24 h following a PN1 injection of the aromatase inhibitor, formestane. These data are displayed in a heat map (Fig. 3). Hierarchical clustering analysis using Pearson correlation as a metric successfully segregated male vehicle samples from formestane-treated male and female vehicle samples, while it was unable to distinguish between formestane-treated male and vehicle-treated female samples (Fig. 3). Formestane significantly increased expression of miR-143 (t(12) = 3.0, p = 0.02), miR-152 (t(12) = 3.69, p = 0.005), miR-18a (t(12) = 2.35, p = 0.04), miR-298 (t(12) = 2.24, p = 0.05), miR-301b (t(12) = 2.29, p = 0.04), miR-34a (t(12) = 2.65, p = 0.03), miR-362–3p (t(12) = 2.79, p = 0.02), miR-365 (t(12) = 2.75, p = 0.02), miR-384–3p (t(12) = 3.34, p = 0.007), miR-448 (t(12) = 2.47, p = 0.03), miR-451 (t(12) = 2.57, p = 0.03), and miR-674 (t(12) = 2.38, p = 0.04). Formestane significantly reduced expression of miR-133b (t(12) = 4.16, p = 0.002), miR-15a (t(12) = 2.23, p = 0.05), miR-467c (t(12) = 2.47, p = 0.03), and miR-671–3p (t(12) = 2.54, p = 0.04).
Epidemiological studies have linked sex-biased neurodevelopmental disorders, including autism and schizophrenia, with prenatal stress (Huttunen and Niskanen, 1978; van Os and Selten, 1998; Khashan et al., 2008; Kinney et al., 2008). Animal models of prenatal and postnatal stress have provided insight into sensitive periods and sex-specific vulnerabilities related to neurodevelopmental disorder etiology (Champagne and Meaney, 2007; Mueller and Bale, 2007; Kapoor and Matthews, 2008; Mueller and Bale, 2008; Cottrell and Seckl, 2009; Biala et al., 2010; Eiland and McEwen, 2010; Ivy et al., 2010; Korosi et al., 2010). We previously identified early gestation as a period sensitive to the sex-specific programming effects of prenatal stress in which male offspring showed a dysmasculinized phenotype in behavioral and physiological stress measures as adults (Mueller and Bale, 2007, 2008). As certain disease outcomes persist into subsequent generations, we examined the paternal transmission and programming of the prenatal stress induced dysmasculinized phenotype in second-generation (F2) offspring.
F2 brain gene expression was examined during the perinatal sensitive period to identify mechanisms of a disruption in masculinization in the F2 male brain. This period is critical for the organization of the sexually dimorphic brain by gonadal hormones. Using a custom Taqman qRT-PCR Array for genes involved in neurodevelopment, we observed a broad shift in expression from a male-typical to a more female-typical pattern in the F2 male offspring of prenatally stressed sires (F2-S). In F2-S male PN1 brains, 13 of 17 genes with statistical sex differences displayed expression levels closer to F2-C females than to F2 control (F2-C) male levels. These data correlate with diminished sex differences previously reported in adult hippocampal gene expression of prenatally stressed rats, supporting the hypothesis that disrupted masculinization during the perinatal sensitive period may be a mechanism through which paternal (F1) prenatal stress exposure impacts F2 offspring development (Biala et al., 2010).
Sex differences in gene expression result from combinations of chromosomal and hormonal effects. The male brain is organized in a sex-specific manner by a surge of testosterone during the perinatal sensitive period (Phoenix et al., 1959; McCarthy et al., 2009a). Testosterone is converted to estradiol by a neuronal-specific aromatase where it alters gene expression to masculinize and defeminize neurocircuitry through the estrogen receptors ERα and ERβ. We examined the expression of these primary effectors to determine whether their dysregulation was associated with the broad shift in gene expression observed in F2-S males. While aromatase expression was unchanged, both ERα and ERβ appeared upregulated, an effect suggestive of reduced ligand availability supporting a hypothesis for decreased perinatal testosterone in F2-S males. To identify potential alternative mediators of the dysmasculinized gene expression in F2-S males, we examined the PN1 brain miRNA environment. miRNAs are small non-protein-coding RNAs involved in the posttranscriptional regulation of genes (Bartel, 2009). Interestingly, a single miRNA may interact with up to a hundred target transcripts, potentially regulating critical gene families involved in early neurodevelopment. We identified 3 miRNAs whose expression appeared dysmasculinized in F2-S males, and 2 miRNAs that showed a significant effect of paternal (F1) prenatal stress. Several of these miRNAs have known functions in peripheral tissues (Caruso et al., 2010; Ghosh et al., 2010; Qin et al., 2010). To determine whether these changes had functional consequences on gene expression, we identified predicted gene targets using the web-based algorithm miRDB (Wang, 2008; Wang and El Naqa, 2008). Only one gene, β-glycan (TGFβr3), was a shared predicted target of all three dysmasculinized miRNAs. As would be predicted based on the reduced expression of miR-322, miR-574-3p, and miR-873, β-glycan expression was significantly increased in the F2-S male PN1 brain. Beta-glycan is a member of the TGFβ superfamily expressed in adult brain, pituitary, and gonadal tissues where it acts as an accessory protein, binding other TGFβ isoforms, such as inhibin A, and increasing their receptor affinity (Lewis et al., 2000; MacConell et al., 2002). Interestingly, in pituitary gonadotrophs and gonadal Leydig or theca cells, β-glycan is involved in regulating the release of gonadal hormones (MacConell et al., 2002; Chapman and Woodruff, 2003; Wiater et al., 2009). As a role for β-glycan in neurodevelopment has not been identified, our data suggest that it may serve an unappreciated role in the organization of the sexually dimorphic brain.
As an additional physiological marker programmed by perinatal testosterone, adult male anogenital distances were measured (Scott et al., 2008). As predicted, F2-S males showed a significantly reduced anogenital distance and adult testis weights, supportive of decreased testosterone exposure during the perinatal sensitive period. Interestingly, studies examining prenatal stress during late pregnancy have also reported decreased perinatal testosterone, adult anogenital distance, and testis weight in rats (Dahlöf et al., 1978; Ward and Weisz, 1980). It is important to note that we are examining these measures in F2 animals that were not themselves exposed to any prenatal manipulation. Thus, these data suggest that reduced exposure to organizational hormones during the perinatal sensitive period is a common marker of gestational stress experience, and that this effect can be transmitted along the paternal lineage to F2 male offspring. Of translational importance, male schizophrenics have been reported to display reduced circulating testosterone and disruptions in brain masculinization (Gur et al., 2004; Goldstein et al., 2007). Further, boys with prepsychotic prodromal symptoms had significantly lower testosterone levels during adolescents, a period of increased psychotic disorder onset (van Rijn et al., 2011).
To examine F2-S male adult dysmasculinized and stress-sensitive phenotypes, we measured their stress responsivity in the tail suspension test, Barnes maze, and HPA stress axis. These tests were selected as they measure predictable sex differences in stress-provoking environments, and performance in these tests was previously found to be significantly dysmasculinized in F1 prenatally stressed males (Mueller and Bale, 2007, 2008). In the tail suspension test, F2-S males spent significantly more time immobile than F2-C males. No effect of paternal prenatal stress was detected in females. These results are similar to those reported in the first generation (Mueller and Bale, 2008). While there was not a statistically significant effect of F2-S on overall performance in the Barnes maze or HPA axis sensitivity, males did show a general trend for a pattern of stress responsivity similar to that identified in F1 prenatally stressed males (Mueller and Bale, 2008). Thus, it appears that aspects of the adult dysmasculinized stress-sensitive phenotype were transmitted from F1 prenatally stressed sires to their F2 male offspring. It is also possible that increased numbers of litters may have provided sufficient statistical power to identify significant effects in additional stress tests.
As our data point to a likely reduction in testosterone-mediated developmental organization in F2-S male brains, we hypothesized that miRNAs in the brain are responsive to organizational gonadal hormones. Therefore, in a subsequent study examining the influences of estradiol on the neonatal brain miRNA environment, we administered the aromatase inhibitor, formestane, to PN1 male neonates. miRNA expression was then assayed using a miRNA Taqman qRT-PCR Array. Aromatase inhibition dramatically dysmasculinized the brain miRNA environment where statistical hierarchical clustering was unable to distinguish between formestane-treated males and control females based on miRNA expression patterns, while completely segregating control males from these groups. Thus, these data confirm the dynamic response of the miRNA environment during this critical window. Gonadal hormones have previously been shown to regulate miRNAs in peripheral target tissues (Klinge, 2009; Delić et al., 2010; Narayanan et al., 2010). However, our data appear to demonstrate a novel impact of organizational hormones on brain miRNA expression during the perinatal sensitive period. Epigenetic mechanisms have been attributed to gonadal hormone status and shown to influence brain sexual differentiation and may intersect with miRNAs to program the sexually dimorphic brain (McCarthy et al., 2009b; Auger and Auger, 2011; Auger et al., 2011).
Our studies provide intriguing evidence for the paternal transmission of prenatal stress effects on neurodevelopmental processes including programming of the miRNA environment and adult stress responsivity. Transmission through the paternal lineage excludes confounds associated with maternal transmission, such as the intrauterine environment or maternal behaviors, and implicates gametic epigenetic mechanisms (Youngson and Whitelaw, 2008). However, we cannot completely discount paternal experience effects that, while unlikely, could occur during the brief time the males are in the cage with females where stress-sensitive F1 males may impart some aspect of their behavior upon the pregnant dam. A recent report examining a postnatal stress model has also demonstrated the ability of early-life maternal separation to alter adult behavior and methylation patterns of several genes in the germ line of male mice, with effects persisting into second-generation offspring (Franklin et al., 2010). In addition, the dysmasculinization we observed in F2-S male offspring importantly points to a developmental window of susceptibility during which the programming effects of early prenatal stress exposure may manifest. As such, identifying developmental processes affected during this window, such as the dynamic changes in miRNAs detected, may lead to critical therapeutic targets or biomarkers predictive for neurodevelopmental diseases, particularly in at-risk pregnancies. Overall, these data support an early gestational period vulnerable to prenatal stress epigenetic programming of the male germline, permitting paternal transmission into subsequent generations.
This work was supported by NIH Grants MH087597 and MH091258. We thank A. Scott, D. Fisher, and C. Taylor for technical assistance, A. Morgan for graphics support, and G. Dunn for scientific discussion.
- Correspondence should be addressed to Dr. Tracy L. Bale, University of Pennsylvania, School of Veterinary Medicine, 3800 Spruce Street, Philadelphia, PA 19104-6046.