- Open Access
Investigation of sex differences in the expression of RORA and its transcriptional targets in the brain as a potential contributor to the sex bias in autism
Molecular Autism volume 6, Article number: 7 (2015)
Autism spectrum disorder (ASD) is a neurodevelopmental condition characterized by significant impairment in reciprocal social interactions and communication coupled with stereotyped, repetitive behaviors and restricted interests. Although genomic and functional studies are beginning to reveal some of the genetic complexity and underlying pathobiology of ASD, the consistently reported male bias of ASD remains an enigma. We have recently proposed that retinoic acid-related orphan receptor alpha (RORA), which is reduced in the brain and lymphoblastoid cell lines of multiple cohorts of individuals with ASD and oppositely regulated by male and female hormones, might contribute to the sex bias in autism by differentially regulating target genes, including CYP19A1 (aromatase), in a sex-dependent manner that can also lead to elevated testosterone levels, a proposed risk factor for autism.
In this study, we examine sex differences in RORA and aromatase protein levels in cortical tissues of unaffected and affected males and females by re-analyzing pre-existing confocal immunofluorescence data from our laboratory. We further investigated the expression of RORA and its correlation with several of its validated transcriptional targets in the orbital frontal cortex and cerebellum as a function of development using RNAseq data from the BrainSpan Atlas of the Developing Human Brain. In a pilot study, we also analyzed the expression of Rora and the same transcriptional targets in the cortex and cerebellum of adult wild-type male and female C57BL/6 mice.
Our findings suggest that Rora/RORA and several of its transcriptional targets may exhibit sexually dimorphic expression in certain regions of the brain of both mice and humans. Interestingly, the correlation coefficients between Rora expression and that of its targets are much higher in the cortex of male mice relative to that of female mice. A strong positive correlation between the levels of RORA and aromatase proteins is also seen in the cortex of control human males and females as well as ASD males, but not ASD females.
Based on these studies, we suggest that disruption of Rora/RORA expression may have a greater impact on males, since sex differences in the correlation of RORA and target gene expression indicate that RORA-deficient males may experience greater dysregulation of genes relevant to ASD in certain brain regions during development.
Autism spectrum disorder (ASD) refers to a group of neurodevelopmental disorders that are diagnosed on the basis of impaired social interactions and communication, especially social reciprocity, and the presence of aberrant, repetitive, and stereotyped behaviors . Because of the strong heritability of ASD based on monozygotic twin and sibling studies [2, 3], there have been extensive searches for genetic mutations and variations that may cause ASD [4–13]. To date, there are hundreds to thousands of genes that are considered autism candidate or susceptibility genes in the autism gene databases [14, 15], with the number of genes continually growing. However, the combined genetic variations account for only approximately 20% of the cases, with no single gene or genetic variation associated with more than 1% of the cases. Thus, the etiology of the majority of ASD cases is still unknown.
Interestingly, ASD is consistently reported as having a higher incidence (approximately 4:1) in males than in females . Yet, the mechanism for the sex bias is unknown. Several hypotheses for the sex bias in ASD include: (1) genetic mechanisms which might involve (a) genes on the X or Y chromosomes; (b) skewed X-inactivation, (c) sex-specific imprinting defects on either the X chromosome or autosomes; (2) the extreme male brain hypothesis which posits that elevated fetal testosterone is a risk factor for ASD; and (3) gene-environment interactions that predispose an individual to ASD. All of these hypotheses have been comprehensively described in several recent reviews [17–19]. To date, none of these hypotheses, discussed briefly in the next paragraph, have been either proven or disproven. Given the clinical and genetic heterogeneity of ASD, it is possible that each of these mechanisms for sex bias may apply to specific cohorts of individuals with ASD. What is clear, however, is that a validated mechanism for sex bias will reveal a fundamental process inherent to the core biology of ASD.
With regard to genetic mechanisms for sex bias, there are a few candidate genes for ASD on the sex chromosomes, but the reported genetic variations in them cannot account for the majority of cases. On the other hand, the extreme male brain hypothesis proposed by Baron-Cohen and colleagues focuses on elevated testosterone levels as a risk factor for ASD. This hypothesis has been investigated primarily by correlating autistic behaviors and traits with fetal testosterone levels in amniotic fluid during gestation of the individuals [17, 20–23]. More recently, this group has investigated differences in brain morphology linked to circulating testosterone levels in adults. Their studies demonstrated both morphological and volumetric changes in the brain that are both sex-dependent  as well as dependent on diagnosis of ASD [25, 26]. Direct investigation of the third hypothesis regarding gene by environment (GxE) interactions that may predispose an individual to ASD has been hampered by lack of knowledge regarding the genes that may be involved and the environmental factors that are relevant to ASD. We suggest that the reciprocal regulatory interaction between retinoic acid-related orphan receptor alpha (RORA) and sex hormones which includes RORA-mediated transcriptional regulation of CYP19A1, described below, may be a candidate for GxE interactions that modulate risk for ASD.
Our recent integrative genomic analyses of ASD involving gene expression and methylation profiling of lymphoblastoid cell lines (LCL) from monozygotic twins and pairs of siblings discordant for autism diagnosis revealed the dysregulated expression of many genes in ASD [27, 28], some of which have been associated with aberrant methylation of the gene promoters . One of the downregulated and hypermethylated genes was the nuclear receptor RORA. Further analysis of this gene (which exhibits steroid hormone receptor binding sites in its promoter region) revealed that its expression was upregulated by estradiol but downregulated by the androgen dihydrotestosterone (DHT) . Moreover, RORA was shown to transcriptionally regulate CYP19A1, a gene coding for aromatase, an enzyme that converts testosterone to estradiol. This intriguing relationship between RORA and sex hormones led us to propose a model in which RORA deficiency may lead to an elevation of testosterone and depletion of estradiol through suppression of CYP19A1 expression. This model provides a plausible biochemical explanation for the elevated testosterone levels seen in studies by the Baron-Cohen group. It also suggests a more direct mechanism for the sex bias in ASD in which normal females, with higher estrogen levels, might exhibit higher expression of RORA, thus buffering against agents that induce RORA deficiency. This model further predicts that during development, RORA expression may be sexually dimorphic at least in some regions of the brain.
This pilot study was undertaken to examine the possibility of sex differences in the expression of RORA in several brain regions at different stages of development of the normal human brain using RNAseq data from the BrainSpan Atlas. In addition, sex differences in RORA and aromatase protein levels were investigated by reanalysis of our published confocal immunofluorescence data from the cortex of both male and female controls and age-matched male and female ASD donors. To reduce heterogeneity in the expression patterns due to the genetics as well as age of the brain donors, we also conducted expression analyses of the cortex and cerebellum of a strain of mice at 3 months of age. Together, these studies suggest sexual dimorphism in the expression of RORA/Rora in some brain regions during certain stages of development as well as a high correlation with the expression of its target genes, especially CYP19A1.
Re-analysis of pre-existing confocal immunofluorescence data from tissue arrays
Data from a previous study  employing confocal immunofluorescence to quantitate protein levels of RORA and aromatase on a tissue array containing specimens of the frontal cortex of ASD cases and age-matched controls were reanalyzed on the basis of gender. The tissue array was obtained through the Autism Tissue Program (San Diego, CA, USA) and kindly provided by Dr. Janine LaSalle (UC Davis). Each array contained 600 μm in diameter × 5 μm thick sections in triplicate from the BA9 region of the frontal cortex of autistic individuals and age- and sex-matched controls as well as samples from individuals with a variety of other neurodevelopmental disorders as previously described . Only specimens from neurotypical controls (without any developmental disorder) were used as ‘Controls’ for this study. The amounts of protein are determined by calculating the ‘mean fluorescence’ for RORA and aromatase across 40 to 50 neurons per sample. Furthermore, comparisons of protein levels between females and males and between cases and controls are made for samples that are matched for donor age within ± 3 years. Additional file 1 contains the ages and mean fluorescence levels of antibody staining for RORA and aromatase in all of the ASD donors and controls from the tissue array.
Analysis of RNAseq data from the BrainSpan atlas of the developing human brain
RNAseq data (RPKM) for gene expression of RORA and that of selected transcriptional targets of RORA in specific brain regions were downloaded from the BrainSpan Atlas of the Developing Human Brain . The downloaded data for the orbital frontal cortex and the cerebellum are presented in Additional file 2. The samples, divided by gender, were grouped into three developmental periods: (1) before birth (BB), which was designated in BrainSpan as ‘pcw’ for post-conception weeks; (2) birth to 18 years of age; and (3) older than 19 years of age. Within each group, the samples were matched for age ± 3 years (or ± 3 pcw for the prenatal samples) for comparisons between males and females or between cases and controls. The validated transcriptional targets of RORA selected for RNAseq data analyses were A2BP1, ITPR1, and NLGN1, which are among the autism susceptibility genes represented in the SFARI gene and AutismKB databases [14, 15]. CYP19A1 was not included in these analyses because the RNAseq values were either zero or too low for confidence. The functions of these genes in the context of ASD will be discussed later.
Extraction of frontal cortex and cerebellum from wild-type C57BL/6 mice
Wild type C57BL/6 mice were obtained from our colony at the IFR 83 Biologie Integrative, UPMC (Paris, France). Animal housing and all handling procedures for this study followed ethical guidelines established by Le Comité National d’Ethique pour les Sciences de la Vie et de la Santé (animal ethics committee for France), in accordance with the European Communities Council Directive 2010/63/EU. Male and female mice (three per group) aged 3 months were euthanized with an overdose of sodium pentobarbital (300 mg/Kg i.p.) and the brain rapidly dissected into RNase-free normal saline (0.9% NaCl). Female mice were not synchronized with respect to estrous cycle before the brains were harvested. The frontal cortex and cerebellum were isolated and placed into 100 μL RNALater (Qiagen, France) and maintained at 4°C for 24 h prior to freezing at -80°C.
RNA isolation and quantitative RT-PCR analysis
Mouse brain tissues were homogenized in a Bullet Blender Homogenizer (Next Advance, Averill Park, NY, USA), after which total RNA was isolated using an RNeasy Mini Kit (Qiagen, Gaithersburg, MD, USA). A total of 1 μg purified total RNA was used for cDNA synthesis using the iScript cDNA Synthesis Kit (BioRad, Hercules, CA, USA) according to the manufacturer’s protocols. The reaction (20 μL) was incubated at 25°C for 5 min, followed by 42°C for 30 min, and ending with 85°C for 5 min. After reverse transcription, the cDNA reaction mixture was diluted to a volume of 50 μL with nuclease-free water and used as a template for qPCR analyses. Real-time PCR analyses were conducted using the Applied Biosystems 7300 Real-Time PCR System (Applied Biosystems, Foster City, CA, USA). Each sample was run in triplicate and the average deviation of the CT values was calculated to assess the consistency of the assays. An average CT deviation of 0.25 or less was considered acceptable for replicates. Primers for RT-qPCR analyses were designed using Primer3 software for mouse Rora as well as A2bp1, Cyp19a1, Itpr1, and Nlgn1. The primer sequences for these mouse genes are listed in Additional file 3. Cyp19a1 expression was not determined for the cerebellum because the expression levels were either zero or too low for confidence. Primers for the human genes have been previously published . The relative quantity of transcripts in each sample was calculated using the standard curve method with 18S RNA expression as a reference.
The Student’s t-test within the StatPac Statistical Program (StatPac, Inc., Bloomington, MN, USA) was used to determine t-values and P values for differences between the means of gene expression in age-matched males and females and between the means of age-matched cases and controls. T-values derived from the t-tests and degrees of freedom (sum of the number of samples in two groups minus 2) were also used to determine Cohen’s d and effect size-r to provide an additional measure of the magnitude of the effect of sex (or autism diagnosis) on protein or gene expression level. Cohen’s d and effect size-r were calculated using the open-access online Effect Size Calculator http://www.uccs.edu/lbecker/index.html made available by Dr. Lee A. Becker (University of Colorado, Colorado Springs, CO, USA). StatPac was also used to determine correlation coefficients and P values for comparisons of gene expression levels for RORA/Rora and each of the transcriptional targets. Two-tailed P values are reported for all statistical analyses.
Power and sample size analyses for the experiments involving comparisons of independent group means (for example, case vs. control or female vs. male) were performed using the open-access online Power/Sample Size Calculator http://www.stat.ubc.ca/~rollin/stats/ssize/n2.html which is kindly provided by Dr. Rollin Brant (University of British Columbia, Vancouver, Canada). For power calculations involving a given number of samples in each group (n), a two-sided test was selected with α set to 0.05. The parameters used for sample size calculation were two-sided test, with α = 0.05 and power = 0.80.
All boxplots were generated using the open-source software BoxPlotR (http://boxplot.tyerslab.com) kindly made available online by the Tyers (IRIC-Universite de Montreal, Montreal, Quebec, Canada) and Rappsilber (University of Edinburgh, UK) Laboratories.
Confocal immunofluorescence data on the human frontal cortex suggests lower RORAand aromatase protein levels in males
We have previously demonstrated reduced expression of RORA and aromatase protein in ASD cases relative to controls using confocal immunofluorescence analyses of tissue arrays containing postmortem frontal cortex specimens . Here, we reanalyze the confocal immunofluorescence data from age-matched male (n = 9) and female (n = 8) controls and show that there is a nominally higher level (1.22-fold) of RORA protein in females relative to males (Figure 1). Similarly, the female-to-male (F/M) ratio for aromatase is 1.21. However, these modest differences are not statistically significant (two-tailed P >0.47 for RORA). There is no apparent sex difference in RORA protein in the cortex of male and female individuals with ASD (F/M = 1.01, P = 0.97, effect size = 0.01), while the F/M ratio for aromatase is 1.61 (P = 0.28). Table 1 summarizes the case-control comparisons of both RORA and aromatase protein levels for the combined (males + females) set of ASD and control specimens as well as for age-matched females and males. While the lower protein levels for the combined cases vs. combined controls are statistically significant as previously reported , the nominally lower levels of RORA and aromatase in tissues from cases are not statistically significant from that of controls when the samples are divided into age-matched females and age-matched males. Because the correlation between RORA and aromatase protein levels was previously found to be highly correlated for the combined samples (R2 = 0.91), we examined the correlation between these protein levels in both control samples and ASD samples as a function of sex. Figure 2 shows that the correlation between RORA and aromatase is very high for both control and ASD male samples (R2 ≥ 0.96; P <0.01) as well as for female controls (R2 = 0.96; P <0.0001), while R2 is only 0.62 (P = 0.11) for samples from ASD females.
Correlation of RORAand target gene expression in human postmortem brain tissues
Because of the strong correlation between RORA and aromatase protein levels on tissue arrays of cortical specimens from males and females, we examined the correlation between the mRNA expression of RORA and that of several validated transcriptional targets of RORA in postmortem tissues from the prefrontal cortex of male controls and ASD donors using RT-qPCR data from an earlier study which investigated genome-wide transcriptional targets of RORA . Table 2 summarizes the results of the correlation analyses for the expression of CYP19A1, A2BP1, ITPR1, and NGLN1 vs. RORA expression. As shown, the correlation coefficients are in the range of 0.72 to 0.99 for the control samples, indicating a relatively strong positive correlation between expression of each of these four genes and RORA expression. Interestingly, the correlation between CYP19A1 and RORA expression is identically high and significant for both male controls and age-matched ASD males (not shown) as was seen for the protein correlation in these groups on tissue arrays (Figure 2), while the correlation between the other three target genes and RORA are much lower for cases (R2 ≤ 0.36).
Expression of RORAand selected transcriptional targets of RORA across brain development based on RNAseq data from BrainSpan
As ASD is a developmental disorder, we examined gender-related developmental changes in RORA through meta-analyses of RORA expression using RNAseq data from the BrainSpan Atlas of Human Development . Here, we analyzed the expression of RORA and its correlation with that of three of its transcriptional targets (A2BP1, ITPR1, and NLGN1) in two brain regions which are involved in ASD: the orbital prefrontal cortex and the cerebellum.
Orbital prefrontal cortex
Although there is suggestive evidence for sexually dimorphic expression of RORA in the frontal cortex (presented above), at present, there is no known biological or anatomical correlate for this difference between females and males. However, because a recent magnetic resonance imaging study on brain morphometric differences between adult females and males reported increased volume in the orbital frontal cortex (OFC) of unaffected females in comparison to age-matched unaffected males , we examined RORA expression in this region across several developmental periods. Figure 3 shows that there are developmental differences but no significant sex differences in the average level of RORA expression in this brain region. There is, however, a suggestive difference in expression of RORA (F/M = 1.40; P = 0.21, effect size = 0.6) in the adult group (≥19 years of age), which is the age group represented in the imaging study by Lai et al.. Table 3 and Additional file 4 show the correlation coefficients and R2 plots, respectively, for the expression of RORA and each of its three target genes in the OFC in the three defined developmental periods. For both females and males, there is very little correlation between RORA expression and that of the three target genes before birth, but relatively high correlation between RORA and all three genes in females in the 0 to 18 years age group. While the high correlation with NLGN1 is maintained in the adult females, the correlation for two genes, A2BP1 and ITPR1, decreases in this group. Males, on the other hand, have a much more variable pattern of correlations for RORA and these target genes in the OFC, in which each of the genes exhibits a high correlation with RORA expression, but in different postnatal periods and, in one case (A2BP1), in a negative direction.
In order to further probe the sexually dimorphic expression of RORA, we also studied during several developmental periods the expression of RORA in the cerebellar cortex where it is known to play key roles in development and function [34–36]. Moreover, cerebellar pathology has been consistently reported in ASD . Figure 4 shows that there were no significant sex differences in RORA expression in the cerebellum of age-matched males and females at any of these developmental periods. Table 4 and Additional file 5 further show the correlation data and graphs, respectively, for the expression of RORA and three of its target genes both before and after birth. Interestingly, the correlation between the prenatal expression levels of RORA and all three genes is significantly high (R2 ≥ 0.96, P ≤0.04) for females, but only significantly high for RORA and A2BP1 in males. The postnatal female pattern of correlations between RORA and the three genes approximately mirrors the trend seen for correlations in the OFC in females, with high and significant correlations in the 0 to 18 years age group, and high correlations of lower significance in the adult group. The pattern of correlations in the male cerebellum after birth is more variable, as it was in the OFC.
Analysis of Roraand transcriptional target gene expression in the frontal cortex and cerebellum of wild-type male and female C57BL/6 mice
Because the genetic heterogeneity of human subjects most likely influences gene expression , thus rendering direct comparisons of expression data from a limited number of postmortem samples (n = 3 to 6) insignificant or at best marginally significant, we conducted pilot studies with male and female wild-type C57BL/6 mice to address the question of sexually dimorphic expression of Rora and that of several of its validated transcriptional targets . Figure 5 shows that, even with only three mice of each sex, there is a significant difference between the expression of Rora, Cyp19a1, and Nlgn1 in the cortex of male and female mice (P = 0.008 - 0.026), with expression of all three genes being higher in the females. The strength of this association between gene expression and sex is also reflected in the large effect sizes (0.85 - 0.93) for these three genes. The expression of A2bp1 and Itpr1 is also nominally higher in females (F/M ratio approximately 1.2), but the P values for sex differences fall below the standard level of significance.
Table 5 and Figure 6 show the results of correlation analyses between the expression of Rora and that of three transcriptional targets in the mouse cortex. It is interesting to note that when separated by sex, the correlation of expression between Rora and each of the four targets represented here is much higher in males than in females, with correlation coefficients in the range of 0.92 to 0.99 (P values: 0.06 to 0.25), and R2 values (an indicator of effect size) in the range of 0.85 to 0.99. Also of interest is the strongly negative correlation coefficient (-0.97) between Cyp19a1 and Rora in females.
With respect to the cerebellum, there were no significant sex differences in the expression of Rora and any of its target genes. This is comparable to the lack of significant sex differences in RORA expression in the human cerebellum. Table 6 and Figure 7 show that the correlation of Rora expression with that of the target genes is high for both males and females (R2 ≥ 0.94), with the exception of A2bp1 in females, for which R2 is 0.64. It should be noted that the high correlation between the expression of Rora and its targets in both brain regions of adult male mice is distinctly different from the more variable patterns seen in the OFC and cerebellum of adult human males. On the other hand, the correlation of Rora expression with all three gene targets in the adult female mouse cerebellum is remarkably similar to that observed for the analogous correlations in the cerebellum of adult human females.
This exploratory study on sex differences in the ASD brain examines the sexually dimorphic expression of RORA, a functionally relevant candidate gene for autism, in the postmortem brain tissues of humans and mice, focusing on the frontal cortex and cerebellum. In addition, we investigated the correlation of RORA (or Rora) expression with several of its validated transcriptional targets in the same brain regions. For this study, we evaluated data for potential sex differences in RORA/Rora expression and/or its correlation with the expression of the transcriptional target genes from four sources: (1) our published confocal immunofluorescence analyses of tissue arrays containing cortical specimens from both male and female donors with ASD and age-matched controls ; (2) our prior gene expression analyses on frozen postmortem brain tissues from male controls and ASD donors ; (3) the collection of RNAseq data from the BrainSpan Atlas of the Developing Human Brain ; and (4) de novo expression analyses of Rora and the orthologous target genes in mouse brain tissues.
Although there was suggestive evidence for sexually dimorphic expression of RORA (F/M ratio approximately 1.2 - 1.4) in the human frontal cortex, these sex differences were not statistically significant. The lack of significance for sex differences in RORA expression may not be surprising given the genetic heterogeneity of human donors and the limited number of tissues available for analyses. In this regard, the effect sizes for sex and ASD influences on RORA expression and protein levels in the frontal cortex may be more informative because they are independent of sample size. Notably, the effect size for RORA protein level in male ASD cases vs. controls is lower (0.22) than the effect size for the female case–control comparison (0.36) (see Table 1). This difference suggests a smaller effect of ASD on RORA levels in males than in females, which may reflect a lower basal level of RORA in unaffected males relative to that in unaffected females. Interestingly, the effect size is virtually zero (0.01) for sex differences in RORA protein in ASD cases (see Figure 1), suggesting that both males and females affected by ASD have a comparably low level of RORA protein in the frontal cortex. In addition, the effect size for sex differences in RORA expression in the OFC of control males and females is at least 0.6 (see Figure 3). Thus, several determinations of effect size suggest that there may be a modest difference in the expression of RORA in the frontal cortex of control males and females which vanishes in males and females with ASD.
Perhaps of more relevance to the sex bias in ASD are sex differences in the correlation of RORA expression with that of some of its target genes in both the cortex and cerebellum. Studies with wild-type adult male and female C57BL/6 mice lend support to the existence of sex differences in RORA expression in the cortex and its correlation with target gene expression in these two brain regions. However, although the correlations between cortical expression of Rora and that of its transcriptional targets are higher in male mice in comparison to female mice, the analogous correlations using RNAseq data for human tissues are more complex and dependent on target gene and developmental period. It is also possible that the variability in correlations for human tissues is in part due to the genetic heterogeneity of the brain donors as well as the fact that the majority of age-matched samples used for the RORA-target expression correlations are not from the same donors. This donor disparity in the BrainSpan Atlas is in contrast to the samples involved in the confocal immunofluorescence and gene expression analyses performed by our laboratory, in which correlations of protein or gene expression levels (Figure 2 and Table 2, respectively) were determined for RORA and targets in the same samples, resulting in consistently higher and positive correlations. The following sections first highlight the principal findings of this study, and then discuss the implications of these findings with respect to the potential contribution of RORA to the sex bias and pathobiology of ASD.
Comparison of RORA and aromatase protein levels in the frontal cortex
Although there were no statistically significant sex differences in RORA and aromatase protein levels in the cortical specimens on the tissue array (Figure 1), the correlation analyses show that the protein levels for both RORA and aromatase in ASD cases (age-matched females and males) fall at the lower end of the range of protein levels for control females and males (Figure 2). However, the correlation between these protein levels is higher for ASD males than for ASD females, which may relate to sex differences in the response to RORA deficiency. Notably, there is a greater percent reduction in aromatase protein in ASD males in comparison to ASD females, suggesting that females may be better able to offset aromatase deficiency caused by RORA deficiency through compensatory mechanisms.
Correlation of gene expression levels in the frontal cortex
Table 2 reveals a significantly high correlation between the expression of RORA and that of CYP19A1 (the gene coding for aromatase) in frozen postmortem cortical tissues from both control and ASD males. This high correlation at the level of gene expression thus reflects the high correlation seen at the protein level for male cases and controls. The strong and positive correlation between CYP19A1 and RORA protein levels suggests that a deficiency in aromatase is likely to result in higher testosterone (and lower estradiol) levels in the presence of RORA deficiency. Physiological disturbances in either of these hormones during critical periods of development may result in impaired neurological development and function, as discussed later.
Analysis of RNAseq data for RORA and several of its transcriptional targets in the orbital frontal cortex (OFC) and cerebellum
RNAseq data from the BrainSpan Atlas was used as an experimental tool to explore sexually dimorphic expression of RORA and that of its target genes in different human brain regions across different developmental periods. Although there were no statistically significant sex differences in RORA expression in the OFC and cerebellum (Figures 3 and 4), there was a suggestively higher expression (and a correspondingly larger effect size) in the OFC of adult females in comparison to age-matched males, with an F/M ratio comparable to that observed for RORA protein in the frontal cortex (Figure 1). Of potentially greater interest are the sex differences in the correlation between RORA expression and that of its transcriptional targets, specifically A2BP1 (RBFOX1), ITPR1, and NLGN1, summarized in Tables 3 and 4. However, the sex-dependent correlations are complex, being dependent on both brain region and period of development. In the OFC before birth where the expression of RORA is very low (approximately 10% of that in the cerebellum), there is virtually no correlation between RORA expression and that of any of the three target genes in either females or males. This suggests that RORA may not play a major role in the prenatal development of this brain region, in contrast to the critical role of RORA in development of the cerebellum , where the level of RORA expression is consistently high in both males and females across all developmental stages. The correlation patterns in the cerebellum are also different from the OFC, with females exhibiting significantly high correlations between RORA expression and that of the target genes through age 18, while the correlation pattern is more variable for males (Table 4), with each target gene showing a high correlation in at least one developmental period. Because some of this variability may be due to the mismatch of tissue donors (though age-matched) for the RNAseq studies which may confound correlation analyses, we performed gene expression analyses of Rora and the orthologous target genes in the cortex and cerebellum of a strain of adult male and female mice.
Sexually dimorphic expression of Rora and its transcriptional targets in the mouse brain
The strongest evidence for sexually dimorphic expression of RORA in the brain is provided by the pilot studies using wild-type adult male and female C57BL/6 mice. Here, we show sex-dependent differences in the cortical expression of Rora and two of its transcriptional targets, Cyp19a1 and Nlgn1, with expression of all three being significantly higher in females, which is also consistent with the large effect sizes for these genes. Power analyses, based on the expression data for these three genes, indicate sufficient power (≥0.98) to detect significant differences, even with only three mice per group. In part, these results using genetically homogeneous female and male mice support the sex differences in RORA protein that were suggested by our confocal immunofluorescence analyses of postmortem tissue specimens from the frontal cortex of males and females as well as that suggested by RNAseq data from the OFC of donors ≥19 years of age. However, there were no significant differences in Rora and target gene expression in the mouse cerebellum, mirroring the results of our analyses of the RNAseq data from the human cerebellum.
Correlation of Rora and target gene expression in the mouse brain
Unlike the variable correlations of RORA and target gene expression in the brain of genetically heterogeneous human males and females, the correlations derived from the adult male and female mouse brain may be more informative and relevant to sex differences that potentially impact brain function. Interestingly, the correlation coefficients are high between Rora and all targets (including Cyp19a1) in the male cortex but variable and generally lower in the female cortex (Table 5). In contrast, the Rora-target gene correlations in the cerebellum are comparably high for both males and females with respect to all targets except for A2bp1 in females (Table 6). Based on the more than 2,500 putative transcriptional targets of RORA that were identified by our previous chromatin immunoprecipitation-promoter chip hybridization (ChIP-on-chip) analyses , the higher correlation between the expression levels of Rora and each of its transcriptional targets suggests that Rora deficiency may have a greater impact on neurological development and function in males than in females, especially in the cortex. This is particularly relevant since over 400 of the identified transcriptional targets of RORA are listed in SFARI Gene and AutismKB databases as autism candidate/susceptibility genes [14, 15]. Moreover, RORA and estrogen receptor (ER) are both nuclear receptors that share the same DNA consensus binding sites and regulate the transcription of some of the same genes . Thus, the higher concentration of estrogen, which has been demonstrated even in prepubescent females versus age-matched males , may offset the loss of transcriptional control due to RORA deficiency by activating ER-mediated transcription, resulting in a lesser impact on neurological development in females. Also intriguing is the strong negative correlation coefficient between Rora and Cyp19a1 expression levels in female mice. In contrast to the expected increase in Cyp19a1 (aromatase) expression with increasing Rora expression which is seen in males, it is possible that aromatase expression in vivo is instead tightly regulated to maintain estrogen homeostasis (that is, to prevent overproduction of estrogen) in females under normal conditions. Such homeostatic regulation of estrogen in females is particularly important inasmuch as RORA and ER share overlapping transcriptional targets.
Complex regulation of RORAexpression in the brain
Regardless of sex differences, the underlying mechanisms driving regional and developmental differences in RORA expression, such as those seen in the OFC, are unclear. Recently, we described differential recruitment of co-regulator proteins that are involved in the sex hormone-dependent regulation of RORA in the SH-SY5Y neuronal cell model . We demonstrated that estrogen-mediated upregulation of RORA expression required both the estrogen receptor (ER-alpha) and the coactivator NCOA5, while androgen (DHT)-mediated downregulation of RORA required the androgen receptor (AR) and the corepressor function of SUMO1. Thus, the regional and temporal patterns of RORA expression in males and females during development may reflect both the tissue-dependent and developmental patterns of co-regulator and/or sex hormone receptor expression in the respective brain regions. Likewise, transcriptional regulation of target genes by RORA, a nuclear receptor whose function requires interaction with co-regulator proteins, is also dependent on the regional and temporal patterns of co-regulator expression.
A model connecting RORAto the sex bias in autism
Based on our previous studies showing the reciprocal regulatory relationship between RORA and the sex hormones, we proposed a mechanistic model in which RORA deficiency may biochemically explain the increase in testosterone in ASD, which has been reported by Baron-Cohen and colleagues [21, 23]. The strong correlation between RORA and aromatase protein and gene expression levels in the human brain from our earlier study  and this study lends support to this model explaining, in part, the higher testosterone levels associated with some cases of ASD. However, the question regarding a direct role for RORA in the sex bias of ASD remains unresolved since there were no significant sex differences in RORA expression, despite a marginally higher level of RORA expression in the frontal cortex of normal adult females. Indeed, whether a 20% to 40% higher level of RORA expression in females can protect against genetically, physiologically, or environmentally-induced RORA deficiency during critical periods of development remains to be determined, possibly by using gene-knockdown animal models. It is noted that the levels of RORA protein in the cortex of male and female individuals with ASD are quite similar (Figure 1), perhaps reinforcing the idea of a ‘higher threshold for ASD in females’ which may include greater tolerance to RORA deficiency down to a level at which some of the pathobiological processes associated with ASD are triggered. At present, it is unclear what the functional consequences might be for sex differences in RORA expression, particularly in the adult OFC. However, given the known neuroprotective function of RORA [42–46], it is worth noting that the reported volume of the OFC is greater in female than in male adults . With respect to sex bias in ASD, perhaps of more importance than RORA expression per se is the correlation between RORA expression and that of its transcriptional targets, the functions of which are described in the next section.
RORA, a molecular link between sex hormones, neurodevelopment, and autism pathobiology
The involvement of sex hormones in brain development has been well-documented [47–50], yet the precise molecular mechanisms through which the hormones affect neuronal processes, such as neurogenesis, neuron migration, synaptogenesis, synaptic plasticity, axon guidance, and dendrite formation, are not well understood. Equally unclear are the mechanisms through which the sex hormones (in particular estradiol) mediate the observed sexually dimorphic effects on dendritic spine formation and synaptogenesis in certain brain regions during development [51, 52]. A reasonable assumption is that these hormones act through their respective hormone receptors. However, although there is little evidence for genetic or functional alterations of the primary estrogen and androgen receptors (ERα and AR) in a significant number of ASD cases, homeostatic imbalance of their circulating hormones (for example, through altered aromatase expression) will inevitably alter receptor activation.
We suggest that RORA may be a molecular link between the sex hormones and neurodevelopment as well as a mediator of at least some of the pathobiological processes associated with autism. While RORA is a known regulator of circadian rhythm  which in turn has been linked to synaptic regulation , it is also a nuclear hormone receptor that is involved in the transcriptional regulation of many genes in different tissues, thus exhibiting tissue-dependent pleiotropic effects. In a neuronal cell model, we have demonstrated that RORA can bind to the promoter regions of over 2,500 genes, 438 of which are included in autism gene databases [14, 15]. Gene ontology analyses of the putative gene targets of RORA revealed significant overrepresentation of genes involved in neuronal differentiation, neuron projection morphogenesis, axonogenesis, and axon guidance. We validated six transcriptional targets of RORA (A2BP1, CYP19A1, HSD17B10, ITPR1, NLGN1, and NTRK2) and, in this study, investigated the correlation between the expression of RORA and that of four of these targets in both the human and mouse brain. With regard to functional relevance, A2BP1, also known as RBFOX1, codes for a neuron-specific splicing factor associated with synaptic transmission, neurodevelopment, and developmental delay [55, 56]. Furthermore, A2BP1/RBFOX1 has been recently highlighted by a whole exome sequencing study as an evolutionarily constrained gene that regulates the splicing of a large number of other genes identified with high-confidence as ASD risk genes . Interestingly, this study, which utilized a novel statistical model that integrated transmission and de novo association (TADA) analyses, also identified NLGN1 as one of the synaptic proteins encoded by the TADA genes, which included a number of other still unvalidated transcriptional targets of RORA (specifically, ANK2, APH1A, CACNA1D, HOMER, MYO9B, NR3C2, and TRIO). These strong genetic associations of putative RORA targets with ASD risk further suggest that RORA deficiency may have a large impact on neuronal functions disrupted in autism. Other neurologically relevant validated targets of RORA include ITPR1, a calcium signaling molecule involved in synaptogenesis, plasticity, dendritic contact, and long-term depression [58, 59], and NTRK2, a neurotrophin kinase also involved in axon guidance, synaptogenesis, plasticity, mood disorder, and learning [60, 61]. On the other hand, CYP19A1 and HSD17B10 both code for metabolic enzymes involved in the conversion of testosterone to estradiol, suggesting that RORA deficiency may inhibit both biochemical pathways, thus exacerbating the expected increase in testosterone or depletion of estradiol. At present, it is not known what homeostatic mechanisms might restore the hormonal balance when both of these pathways are inhibited. Moreover, although our model predicts that testosterone levels would increase when RORA (and aromatase) are decreased, there is no information to date on the sex hormone levels in brain tissues that are RORA-deficient. It should be mentioned that, although the sex hormone status and possible hormonal imbalance in the heterozygous staggerer (Rora+/sg) mouse has been discussed at length by Doulazmi et al. , our study on the regulation of RORA by sex hormones as well as the regulation of CYP19A1 by RORA only recently demonstrated these associations in a neuronal cell model . With respect to neurological functions, CYP19A1 (aromatase) is associated with neurogenesis, neuronal differentiation, synaptic plasticity, and social cognition [62–65], while HSD17B10 is associated with mitochondrial integrity, mental retardation, and language impairment [66–70], therefore extending the impact of their dysregulation beyond endocrine metabolism. Thus, the neuronal activities and high level neurological functions associated with these validated transcriptional targets of RORA suggest a mechanism for induction of autism brain pathology driven by sex hormones under conditions of RORA deficiency. Finally, we suggest that RORA deficiency may be the direct result of genetic  and epigenetic modifications of RORA, and/or gene-environment interactions. With respect to GxE interactions, RORA expression may be dysregulated by intrinsic metabolic or physiological conditions (for example, sex hormone imbalance due to alterations in the steroid hormone biosynthetic pathways [20, 28, 72]) or by extrinsic environmental factors, such as endocrine disrupting compounds, which are known to interfere with normal hormonal signaling [73–76].
Limitations and future directions
While this exploratory study provides suggestive evidence for sexually dimorphic expression of RORA in certain brain regions during development, the major limitation is a lack of sufficient postmortem brain samples to reveal statistically significant expression differences between males and females. Power analyses, based on the human data reported in this study, indicate that at least 121 samples per group would be necessary to significantly detect a 20% difference in the level of RORA expression between males and females with a power of 0.8. The lack of sufficient samples is further compounded by the need to divide the available samples into subgroups to evaluate sex-dependent differences in gene expression across development.
Another confounding factor is the genetic and phenotypic heterogeneity of the brain donors as there is ample evidence that genetics influences gene expression which, in turn, controls phenotype. Indeed, aside from not having ASD or any other diagnosed neurodevelopmental disorder, the behavioral phenotypes of the control brain donors are unknown. This pilot study using BrainSpan data will help to direct attention to certain brain regions and developmental stages in which sex differences in the expression of RORA and/or regulation of its target genes are suggested. This study also suggests that any study of sex-dependent differences in gene expression in the brain must take into account both regional and developmental changes in gene expression.
Regarding the issue of genetic heterogeneity, mouse strains are much more homogeneous with respect to genotype, but individual expression differences are still apparent even within a strain. Some of these differences, especially in post-pubertal females, may be due to hormonal cycling, so future studies on sex differences in gene expression should utilize more hormonally-synchronized female mice. At present, it is not possible to predict how estrous synchronization of female mice would affect sex hormone levels in the brain and the impact on neuronal expression of Rora. It is further noted that the effects of sex steroids in the brain are also mediated by the local synthesis of neurosteroids , and may not be directly correlated to circulating hormonal levels. Interestingly, there is greater variance in RORA expression among the brain tissue samples from human females included in this study, especially in the later stages of development. In analogy to studies with human brain tissues, studies using animal models to study sex differences in gene expression should also include different brain regions and developmental periods.
Finally, despite the obvious advantages of studying sex-dependent differences in gene expression in an animal model where genotype and other physiological conditions can be better controlled to reduce heterogeneity, it should be kept in mind that regulation of gene expression in the human brain is likely to be much more complex than in mouse brain. This complexity is in part suggested by the differences in correlation coefficients for the expression of RORA/Rora and its respective target genes in male humans and male mice, where the correlation coefficients are generally higher in the mouse. Some of these differences may be due to species-dependent differential expression and recruitment of co-regulators or hormone receptors, and/or species differences in epigenetic regulation and alternative splicing.
Through analyses of pre-existing confocal immunofluorescence and gene expression data from our laboratory and publicly available RNAseq data, we present suggestive evidence that RORA may exhibit sex-dependent differences in gene expression in the human brain that are dependent on both brain region as well as stage of development. Significant sex differences in Rora and target gene expression are more readily observed in the cortex of a genetically homogeneous mouse model. With respect to gene expression in the mouse cortex, the stronger correlation between Rora and target gene expression in male mice in comparison to female mice suggests that Rora deficiency may have a greater impact on downstream events affecting neurological development and function in males in comparison to females. While sex differences in the correlation of RORA expression with that of its gene targets in the human brain are more complex, the consistently strong correlation between RORA and CYP19A1 protein and gene expression levels in our small number of postmortem human samples supports our proposed model for the involvement of RORA deficiency in the higher testosterone levels associated with increased risk for ASD.
Ataxin 2 binding protein 1 (Also known as RBFOX1) *(mouse homologue of genes in parenthesis)
Autism spectrum disorder
Chromatin immunoprecipitation followed by microarray analysis
Cytochrome P450, family 19, subfamily A, polypeptide 1
Hydroxysteroid (17-beta) dehydrogenase 10
Inositol 1,4,5-trisphosphate receptor, type 1
Lymphoblastoid cell line
Neurotrophic tyrosine kinase, receptor, type 2
Quantitative polymerase chain reaction
Retinoic acid receptor-related orphan receptor alpha
Human neuroblastoma cell line.
American Psychiatric Association: Diagnostic and Statistical Manual of Mental Disorders. 5th edition. Arlington, VA: American Psychiatric Association; 2013.
Bailey A, Le Couteur A, Gottesman I, Bolton P, Simonoff E, Yuzda E, et al.: Autism as a strongly genetic disorder: evidence from a British twin study.Psychol Med. 1995, 25:63–77. 10.1017/S0033291700028099
Folstein S, Rutter M: Infantile autism: a genetic study of 21 twin pairs.J Child Psychol Psychiatry. 1977, 18:297–321. 10.1111/j.1469-7610.1977.tb00443.x
State MW, Levitt P: The conundrums of understanding genetic risks for autism spectrum disorders.Nat Neurosci. 2011, 14:1499–506. 10.1038/nn.2924
Weiss LA, Arking DE, Daly MJ, Chakravarti A, Brune CW, West K, et al.: A genome-wide linkage and association scan reveals novel loci for autism.Nature. 2009, 461:802–8. 10.1038/nature08490
Ma D, Salyakina D, Jaworski JM, Konidari I, Whitehead PL, Andersen AN, et al.: A genome-wide association study of autism reveals a common novel risk locus at 5p14.1.Ann Hum Genet. 2009, 73:263–73. 10.1111/j.1469-1809.2009.00523.x
Anney R, Klei L, Pinto D, Regan R, Conroy J, Magalhaes TR, et al.: A genome-wide scan for common alleles affecting risk for autism.Hum Mol Genet. 2010, 19:4072–82. 10.1093/hmg/ddq307
Wang K, Zhang H, Ma D, Bucan M, Glessner JT, Abrahams BS, et al.: Common genetic variants on 5p14.1 associate with autism spectrum disorders.Nature. 2009, 459:528–33. 10.1038/nature07999
Sanders SJ, Murtha MT, Gupta AR, Murdoch JD, Raubeson MJ, Willsey AJ, et al.: De novo mutations revealed by whole-exome sequencing are strongly associated with autism.Nature. 2012, 484:237–41. 10.1038/nature10923
O’Roak BJ, Vives L, Girirajan S, Karakoc E, Krumm N, Coe BP, et al.: Sporadic autism exomes reveal a highly interconnected protein network of de novo mutations.Nature. 2012, 484:246–50. 10.1038/nature10897
Neale BM, Kou Y, Liu L, Ma’Ayan A, Samocha KE, Sabo A, et al.: Patterns and rates of exonic de novo mutations in autism spectrum disorders.Nature. 2012, 484:242–5. 10.1038/nature10928
Chahrour MH, Yu TW, Lim ET, Ataman B, Coulter ME, Hill RS, et al.: Whole-exome sequencing and homozygosity analysis implicate depolarization-regulated neuronal genes in autism.PLoS Genet. 2012, 8:e1002635. 10.1371/journal.pgen.1002635
Jiang Y-H, Yuen RKC, Jin X, Wang M, Chen N, Wu X, et al.: Detection of clinically relevant genetic variants in autism spectrum disorder by whole-genome sequencing.Am J Hum Genet. 2013, 93:249–63. 10.1016/j.ajhg.2013.06.012
Basu SN, Kollu R, Banerjee-Basu S: AutDB: A gene reference resource for autism research.Nucleic Acids Res. 2009,37(Suppl 1):D832–6.
Xu L-M, Li J-R, Huang Y, Zhao M, Tang X, Wei L: AutismKB: An evidence-based knowledgebase of autism genetics.Nucleic Acids Res. 2012, 40:D1016–22. 10.1093/nar/gkr1145
Fombonne E: The epidemiology of pervasive developmental disorders. In The neurobiology of autism. Edited by: Bauman ML, Kemper T. Baltimore, MD: The Johns Hopkins University Press; 2005:3–22.
Baron-Cohen S, Lombardo MV, Auyeung B, Ashwin E, Chakrabarti B, Knickmeyer R: Why are autism spectrum conditions more prevalent in males?PLoS Biol. 2011, 9:e1001081. 10.1371/journal.pbio.1001081
Schaafsma SM, Pfaff DW: Etiologies underlying sex differences in autism spectrum disorders.Front Neuroendocrinol. 2014, 35:255–71. 10.1016/j.yfrne.2014.03.006
Werling DM, Geschwind DH: Sex differences in autism spectrum disorders.Curr Opin Neurol. 2013, 26:146–53. 10.1097/WCO.0b013e32835ee548
Baron-Cohen S, Auyeung B, Nørgaard-Pedersen B, Hougaard DM, Abdallah MW, Melgaard L, et al.: Elevated fetal steroidogenic activity in autism.Mol Psychiatry 2014. doi:10.1038/mp.2014.48
Knickmeyer RC, Baron-Cohen S: Fetal testosterone and sex differences in typical social development and in autism.J Child Neurol. 2006, 21:825–45. 10.1177/08830738060210101601
Auyeung B, Baron-Cohen S, Chapman E, Knickmeyer R, Taylor K, Hackett G: Foetal testosterone and the child systemizing quotient.Eur J Endocrinol. 2006, 155:S123–30. 10.1530/eje.1.02260
Auyeung B, Baron-Cohen S, Ashwin E, Knickmeyer R, Taylor K, Hackett G: Fetal testosterone and autistic traits.Br J Psychol. 2009, 100:1–22. 10.1348/000712608X311731
Lombardo MV, Ashwin E, Auyeung B, Chakrabarti B, Taylor K, Hackett G, et al.: Fetal testosterone influences sexually dimorphic gray matter in the human brain.J Neurosci. 2012, 32:674–80. 10.1523/JNEUROSCI.4389-11.2012
Chura LR, Lombardo MV, Ashwin E, Auyeung B, Chakrabarti B, Bullmore ET, et al.: Organizational effects of fetal testosterone on human corpus callosum size and asymmetry.Psychoneuroendocrinology. 2010, 35:122–32. 10.1016/j.psyneuen.2009.09.009
Lai M-C, Lombardo MV, Suckling J, Ruigrok ANV, Chakrabarti B, Ecker C, et al.: Biological sex affects the neurobiology of autism.Brain. 2013, 136:2799–815. 10.1093/brain/awt216
Hu VW, Frank BC, Heine S, Lee NH, Quackenbush J: Gene expression profiling of lymphoblastoid cell lines from monozygotic twins discordant in severity of autism reveals differential regulation of neurologically relevant genes.BMC Genomics. 2006, 7:118. 10.1186/1471-2164-7-118
Hu VW, Nguyen A, Kim KS, Steinberg ME, Sarachana T, Scully MA, et al.: Gene expression profiling of lymphoblasts from autistic and nonaffected sib pairs: altered pathways in neuronal development and steroid biosynthesis.PLoS One. 2009, 4:e5775. 10.1371/journal.pone.0005775
Nguyen A, Rauch TA, Pfeifer GP, Hu VW: Global methylation profiling of lymphoblastoid cell lines reveals epigenetic contributions to autism spectrum disorders and a novel autism candidate gene, RORA, whose protein product is reduced in autistic brain.FASEB J. 2010, 24:3036–51. 10.1096/fj.10-154484
Sarachana T, Xu M, Wu R-C, Hu VW: Sex hormones in autism: androgens and estrogens differentially and reciprocally regulate RORA, a novel candidate gene for autism.PLoS One. 2011, 6:e17116. 10.1371/journal.pone.0017116
Nagarajan RP, Hogart AR, Gwye Y, Martin MR, LaSalle JM: Reduced MeCP2 expression is frequent in autism frontal cortex and correlates with aberrant MECP2 promoter methylation.Epigenetics. 2006, 1:e1–11. 10.4161/epi.1.1.2701
Miller JA, Ding SL, Sunkin SM, Smith KA, Ng L, Szafer A, et al.: Transcriptional landscape of the prenatal human brain.Nature. 2014, 508:199–206. 10.1038/nature13185
Sarachana T, Hu VW: Genome-wide identification of transcriptional targets of RORA reveals direct regulation of multiple genes associated with autism spectrum disorder.Mol Autism. 2013, 4:14. 10.1186/2040-2392-4-14
Gold DA, Gent PM, Hamilton BA: ROR alpha in genetic control of cerebellum development: 50 staggering years.Brain Res. 2007, 1140:19–25.
Steinmayr M, Andre E, Conquet F, Rondi-Reig L, Delhaye-Bouchaud N, Auclair N, et al.: Staggerer phenotype in retinoid-related orphan receptor alpha-deficient mice.Proc Natl Acad Sci U S A. 1998, 95:3960–5. 10.1073/pnas.95.7.3960
Mariani J, Changeux JP: Multiple innervation of Purkinje cells by climbing fibers in the cerebellum of the adult staggerer mutant mouse.J Neurobiol. 1980, 11:41–50. 10.1002/neu.480110106
Fatemi SH, Aldinger KA, Ashwood P, Bauman ML, Blaha CD, Blatt GJ, et al.: Consensus paper: pathological role of the cerebellum in autism.Cerebellum. 2012, 11:777–807. 10.1007/s12311-012-0355-9
Cheung VG, Spielman RS: The genetics of variation in gene expression.Nat Genet 2002,32(Suppl 5):522–5.
Jin VX, Leu Y-W, Liyanarachchi S, Sun H, Fan M, Nephew KP, et al.: Identifying estrogen receptor α target genes using integrated computational genomics and chromatin immunoprecipitation microarray.Nucleic Acids Res. 2004, 32:6627–35. 10.1093/nar/gkh1005
Klein KO, Baron J, Colli MJ, McDonnell DP, Cutler GB Jr: Estrogen levels in childhood determined by an ultrasensitive recombinant cell bioassay.J Clin Invest. 1994, 94:2475–80. 10.1172/JCI117616
Sarachana T, Hu VW: Differential recruitment of coregulators to the RORA promoter adds another layer of complexity to gene (dys) regulation by sex hormones in autism.Mol Autism. 2013, 4:39. 10.1186/2040-2392-4-39
Boukhtouche F, Vodjdani G, Jarvis CI, Bakouche J, Staels B, Mallet J, et al.: Human retinoic acid receptor-related orphan receptor α1 overexpression protects neurones against oxidative stress-induced apoptosis.J Neurochem. 2006, 96:1778–89. 10.1111/j.1471-4159.2006.03708.x
Janmaat S, Akwa Y, Doulazmi M, Bakouche J, Gautheron V, Liere P, et al.: Age-related Purkinje cell death is steroid dependent: RORα haplo-insufficiency impairs plasma and cerebellar steroids and Purkinje cell survival.Age. 2011, 33:565–78. 10.1007/s11357-010-9203-3
Boukhtouche F, Doulazmi M, Frederic F, Dusart I, Brugg B, Mariani J: RORα, a pivotal nuclear receptor for Purkinje neuron survival and differentiation: from development to ageing.Cerebellum. 2006, 5:97–104. 10.1080/14734220600750184
Doulazmi M, Frederic F, Lemaigre-Dubreuil Y, Hadj-Sahraoui N, Delhaye-Bouchaud N, Mariani J: Cerebellar Purkinje cell loss during life span of the heterozygous staggerer mouse (Rora(+)/Rora(sg)) is gender-related.J Comp Neurol. 1999, 411:267–73. 10.1002/(SICI)1096-9861(19990823)411:2<267::AID-CNE7>3.0.CO;2-4
Jolly S, Journiac N, Naudet F, Gautheron V, Mariani J, Vernet-der GB: Cell-autonomous and non-cell-autonomous neuroprotective functions of RORa in neurons and astrocytes during hypoxia.J Neurosci. 2011, 31:14314–23. 10.1523/JNEUROSCI.1443-11.2011
McCarthy MM: Estradiol and the developing brain.Physiol Rev. 2008, 88:91–124. 10.1152/physrev.00010.2007
Wright CL, Schwarz JS, Dean SL, McCarthy MM: Cellular mechanisms of estradiol-mediated sexual differentiation of the brain.Trends Endocrinol Metab. 2010, 21:553–61. 10.1016/j.tem.2010.05.004
Haraguchi S, Sasahara K, Shikimi H, Honda S-I, Harada N, Tsutsui K: Estradiol promotes purkinje dendritic growth, spinogenesis, and synaptogenesis during neonatal life by inducing the expression of BDNF.Cerebellum. 2012, 11:416–7. 10.1007/s12311-011-0342-6
Sasahara K, Shikimi H, Haraguchi S, Sakamoto H, Honda S, Harada N, et al.: Mode of action and functional significance of estrogen-inducing dendritic growth, spinogenesis, and synaptogenesis in the developing Purkinje cell.J Neurosci. 2007, 27:7408–17. 10.1523/JNEUROSCI.0710-07.2007
Brandt N, Vierk R, Rune GM: Sexual dimorphism in estrogen-induced synaptogenesis in the adult hippocampus.Int J Dev Biol. 2013, 57:351–6. 10.1387/ijdb.120217gr
Matsumoto A: Synaptic plasticity in the sex steroid-sensitive neuroendocrine brain.Kaibogaku Zasshi. 1991, 66:81–98.
Sato TK, Panda S, Miraglia LJ, Reyes TM, Rudic RD, McNamara P, et al.: A functional genomics strategy reveals Rora as a component of the mammalian circadian clock.Neuron. 2004, 43:527–37. 10.1016/j.neuron.2004.07.018
Bourgeron T: The possible interplay of synaptic and clock genes in autism spectrum disorders.Cold Spring Harb Symp Quant Biol. 2007, 72:645–54. 10.1101/sqb.2007.72.020
Gehman LT, Stoilov P, Maguire J, Damianov A, Lin C, Shiue L, et al.: The splicing regulator Rbfox1 (A2BP1) controls neuronal excitation in the mammalian brain.Nat Genet. 2011, 43:706–11. 10.1038/ng.841
Fogel BL, Wexler E, Wahnich A, Friedrich T, Vijayendran C, Gao F, et al.: RBFOX1 regulates both splicing and transcriptional networks in human neuronal development.Hum Mol Genet. 2012, 21:4171–86. 10.1093/hmg/dds240
De Rubeis S, He X, Goldberg AP, Poultney CS, Samocha K, Cicek AE, et al.: Synaptic, transcriptional and chromatin genes disrupted in autism.Nature. 2014, 515:209–15. 10.1038/nature13772
Banerjee S, Hasan G: The InsP3 receptor: its role in neuronal physiology and neurodegeneration.Bioessays. 2005, 27:1035–47. 10.1002/bies.20298
Hernjak N, Slepchenko BM, Fernald K, Fink CC, Fortin D, Moraru II, et al.: Modeling and analysis of calcium signaling events leading to long-term depression in cerebellar Purkinje cells.Biophys J. 2005, 89:3790–806. 10.1529/biophysj.105.065771
Huang EJ, Reichardt LF: Trk receptors: roles in neuronal signal transduction.Annu Rev Biochem. 2003, 72:609–42. 10.1146/annurev.biochem.72.121801.161629
Correia CT, Coutinho AM, Sequeira AF, Sousa IG, Lourenco Venda L, Almeida JP, et al.: Increased BDNF levels and NTRK2 gene association suggest a disruption of BDNF/TrkB signaling in autism.Genes Brain Behav. 2010, 9:841–8. 10.1111/j.1601-183X.2010.00627.x
Garcia-Segura LM: Aromatase in the brain: not just for reproduction anymore.J Neuroendocrinol. 2008, 20:705–12. 10.1111/j.1365-2826.2008.01713.x
Garcia-Segura LM, Balthazart J: Steroids and neuroprotection: new advances.Front Neuroendocrinol. 2009, 30:v-ix. 10.1016/j.yfrne.2009.04.006
Rune GM, Frotscher M: Neurosteroid synthesis in the hippocampus: role in synaptic plasticity.Neuroscience. 2005, 136:833–42. 10.1016/j.neuroscience.2005.03.056
El-Emam Dief A, Caldwell JD, Jirikowski GF: Colocalization of P450 aromatase and oxytocin immunostaining in the rat hypothalamus.Horm Metab Res. 2013, 45:273–6.
Yang S-Y, He X-Y, Miller D: Hydroxysteroid (17β) dehydrogenase X in human health and disease.Mol Cell Endocrinol. 2011, 343:1–6. 10.1016/j.mce.2011.06.011
Lenski C, Kooy RF, Reyniers E, Loessner D, Wanders RJ, Winnepenninckx B, et al.: The reduced expression of the HADH2 protein causes X-linked mental retardation, choreoathetosis, and abnormal behavior.Am J Hum Genet. 2007, 80:372–7. 10.1086/511527
Ofman R, Ruiter JP, Feenstra M, Duran M, Poll-The BT, Zschocke J, et al.: 2-Methyl-3-hydroxybutyryl-CoA dehydrogenase deficiency is caused by mutations in the HADH2 gene.Am J Hum Genet. 2003, 72:1300–7. 10.1086/375116
Reyniers E, Van Bogaert P, Peeters N, Vits L, Pauly F, Fransen E, et al.: A new neurological syndrome with mental retardation, choreoathetosis, and abnormal behavior maps to chromosome Xp11.Am J Hum Genet. 1999, 65:1406–12. 10.1086/302638
Edens AC, Lyons MJ, Duron RM, Dupont BR, Holden KR: Autism in two females with duplications involving Xp11.22-p11.23.Dev Med Child Neurol. 2011, 53:463–6. 10.1111/j.1469-8749.2010.03909.x
Tsang KM, Croen LA, Torres AR, Kharrazi M, Delorenze GN, Windham GC, et al.: A genome-wide survey of transgenerational genetic effects in autism.PLoS One. 2013, 8:e76978. 10.1371/journal.pone.0076978
Chakrabarti B, Dudbridge F, Kent L, Wheelwright S, Hill-Cawthorne G, Allison C, et al.: Genes related to sex steroids, neural growth, and social-emotional behavior are associated with autistic traits, empathy, and Asperger syndrome.Autism Res. 2009, 2:157–77. 10.1002/aur.80
Colborn T: A case for revisiting the safety of pesticides: a closer look at neurodevelopment.Environ Health Perspect. 2006, 114:10–7. 10.1289/ehp.7940
Guerrero-Bosagna C, Skinner MK: Environmentally induced epigenetic transgenerational inheritance of phenotype and disease.Mol Cell Endocrinol. 2012, 354:3–8. 10.1016/j.mce.2011.10.004
Vandenberg LN, Colborn T, Hayes TB, Heindel JJ, Jacobs DR, Lee D, et al.: Hormones and endocrine-disrupting chemicals: low-dose effects and nonmonotonic dose responses.Endocr Rev. 2012, 33:378–455. 10.1210/er.2011-1050
Hu VW: Is retinoic acid-related orphan receptor-alpha (RORA) a target for gene-environment interactions contributing to autism?Neurotoxicology. 2012, 33:1434–5. 10.1016/j.neuro.2012.07.009
We wish to thank the families of all of the brain donors for their invaluable gift to autism research and the Autism Tissue Program (San Diego, CA, USA) for making frozen human brain tissues available for this work through the Harvard Brain Tissue Resource Center (Harvard University, Boston, MA, USA). We also thank the Allen Institute for Brain Science (Seattle, WA, USA) and all the members of the consortium working on the brain mapping project for making the RNAseq expression data of the human brain publicly available.
This study was supported in part by generous gifts from the Ireland Family Foundation (Chapel Hill, NC, USA) and the LIFE Foundation (Aspen, CO, USA), and by a grant from the NIEHS (R21 ES023061-01). None of the funding sources played any role in the study design, collection, analysis, and interpretation of data, writing of the manuscript, or decision to submit this study for publication. TS was a predoctoral student in the Institute for Biomedical Sciences at the George Washington University, who was supported by the Higher Educational Strategic Scholarship for Frontier Research Network (SFR scholarship) from the Office of the Commission on Higher Education of the Royal Thai Government, Thailand, through the Faculty of Allied Health Sciences, Chulalongkorn University.
The authors declare that they have no competing interests.
VWH conducted the re-analyses of the confocal immunofluorescence data, the meta-analyses of RNAseq data from BrainSpan, and prepared the manuscript, with input and feedback from TS, RMS, and KMK, who all approved the final manuscript. TS performed the original confocal analyses and RT-qPCR analyses of postmortem human brain tissues as well as the RT-qPCR analyses of Rora and target genes in the mouse brain tissues. RMS contributed the mouse tissues and KMK assisted with the downloading and organization of the BrainSpan data. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Demographic and mean fluorescence data for tissue samples analyzed by confocal immunofluorescence. (PDF 59 KB)
Additional file 2: RNAseq data from the BrainSpan Atlas of the Developing Human Brain which was used for the meta-analyses of gene expression of RORA and transcriptional targets in the OFC and cerebellum. (XLSX 18 KB)
Additional file 3: Primer sequences for qPCR analyses of Rora and its transcriptional targets in mouse tissues. (PDF 56 KB)
Additional file 4:RORA-target gene expression in the OFC.(PDF 539 KB)
Additional file 5:(PDF 559 KB)
About this article
Cite this article
Hu, V.W., Sarachana, T., Sherrard, R.M. et al. Investigation of sex differences in the expression of RORA and its transcriptional targets in the brain as a potential contributor to the sex bias in autism. Molecular Autism 6, 7 (2015). https://doi.org/10.1186/2040-2392-6-7
- RORA expression
- Transcriptional targets
- Sex differences
- Postmortem brain tissues
- Mouse brain