Skip to main content

Possible sexually dimorphic role of miRNA and other sncRNA in ASD brain

Abstract

Background

Autism spectrum disorder (ASD) is sexually dimorphic in brain structure, genetics, and behaviors. In studies of brain tissue, the age of the population is clearly a factor in interpreting study outcome, yet sex is rarely considered. To begin to address this issue, we extend our previously published microarray analyses to examine expression of small noncoding RNAs (sncRNAs), including microRNAs (miRNAs), in ASD and in the control temporal cortex in males and females. Predicted miRNA targets were identified as well as the pathways they overpopulate.

Findings

After considering age, sexual dimorphism in ASD sncRNA expression persists in the temporal cortex and in the patterning that distinguishes regions. Among the sexually dimorphic miRNAs are miR-219 and miR-338, which promote oligodendrocyte differentiation, miR-125, implicated in neuronal differentiation, and miR-488, implicated in anxiety. Putative miRNA targets are significantly over-represented in immune and nervous system pathways in both sexes, consistent with previous mRNA studies. Even for common pathways, the specific target mRNAs are often sexually dimorphic. For example, both male and female target genes significantly populate the Axonal Guidance Signaling pathway, yet less than a third of the targets are common to both sexes.

Conclusions

Our findings of sexual dimorphism in sncRNA levels underscore the importance of considering sex, in addition to age, when interpreting molecular findings on ASD brain.

Background

Autism spectrum disorder (ASD) is one of a number of neurodevelopmental disorders that display sexual dimorphism, occurring more frequently in males, which affect brain structure, gene expression, pathways, function, and ultimately behaviors that will require individualized treatments [13]. Extensive evidence demonstrates sex differences in ASD brain [47]. Sex chromosomes may play a role, with the Y chromosome being a possible risk factor for ASD and X chromosome perhaps having a protective effect [3]. Females appear to have a higher threshold for being affected by genetic factors than males, thus requiring a greater genetic burden, and may have greater brain plasticity [8, 9]. Females carry a higher proportion of de novo CNVs (copy number variants) than males, the CNVs in females disrupt a larger number of genes than in males, and females carry a greater number of de novo single nucleotide variants (SNVs) than males [10, 11]. Environmental and hormone factors that differ between the sexes, like testosterone, may impact the time course and severity of symptoms [12].

Few molecular studies of ASD brain tissue to date have considered potential sexual dimorphism often because of limited tissue availability of female cases. In typical brain development, male-biased gene expression changes are enriched for extracellular matrix, immune response, chromatin, and cell cytoskeleton pathways that have been implicated in ASD [13]. Sex differences in microRNA (miRNA) expression in the frontal cortex have also been described in typical neurodevelopment [14]. Specific genetic mechanisms, such as the expression of retinoic acid-related orphan receptor alpha (RORA) in the frontal cortex, which regulates CYFIP1, may be related to elevated testosterone levels and a potential contributor to the sex bias [15]. There is no clear evidence to date for systematic sex-differential expression of ASD risk genes in human brain; however, genes expressed at higher levels in males are significantly enriched for genes upregulated in postmortem autistic brain, including astrocyte and microglia markers [16].

Here, we extend our analyses of previously published microarray data [17, 18] to examine sexual dimorphism of microRNA and other small noncoding RNA (sncRNA) in male and female ASD and control brain tissue. We focused on two temporal cortical regions: the superior temporal sulcus (STS), a region implicated in social impairments in ASD [1922], and the primary auditory cortex (PAC). Predicted miRNA targets were identified as well as pathways in which they over-populated. We then evaluated our findings in relation to the SFARI database and two gene expression studies in ASD brain: Ziats and Rennert, 2013 [13] and Werling et al., 2015 [16].

Methods

Brain samples

The methods used in the current study are similar to our recent publications [17, 18] (Table 1). Briefly, a total of 34 samples were obtained from 10 ASD and 8 control subjects. The PAC sample was taken from the crown of Heschl’s gyrus and included Brodmann’s areas, 41 and 42 (Fig. 1). The STS samples included Brodmann’s area 22 and were taken from the upper wall of the STS opposite Heschl’s gyrus (Fig. 1). Postmortem tissue integrity, processing, and microarrays are as described in our prior publications on this data set [17, 18].

Table 1 Subject characteristics for ASD (autism spectrum disorders) and controls
Fig. 1
figure 1

Dysregulated sncRNA in ASD females vs control females and ASD males vs control males for the STS (superior temporal sulcus) and the PAC (primary auditory cortex)

Statistics

We used a mixed regression model including diagnosis, sex, age, region, and subject, and group*sex and group*sex*region interactions. To obtain estimates of variance components for our mixed model, we used restricted maximum likelihood estimation (REML), which is suitable for unbalanced designs (W.A. Thompson, 1962). REML optimizes the parameter estimates for the effects in the model. We used repeated measures to account for multiple regions measured within the same subjects (STS and PAC), and we made the SubjectID a random effect. Considering that our sample size was unbalanced between males and females, some identified differences may be due to this limitation. Statistical significance was defined as P ≤ 0.005 and absolute fold change ≥1.2. To account for normal sexual dimorphism, we compared ASD female to control female and ASD male to control male and overlapped the results to identify sexually dimorphic sncRNA expression specifically in ASD.

Pathway analysis and overlap with other studies

Ingenuity pathway analysis (IPA) was used to identify predicted miRNA targets of mature miRNA of major isoforms with sexual dimorphism and to identify over-represented pathways and biological functions. Targets of stem-loop precursor miRNAs (pre-miRNA) were not considered in this pathway analysis, since only mature miRNA can affect gene expression. A Fisher exact test with a Benjamini-Hochberg corrected P < 0.05 was used for determining significance of pathway enrichment. The sexually dimorphic predicted mRNA target genes and pathways were compared to three other data sets: (1) autism risk genes from SFARI using the hypergeometric probability function phyper in R with population size set to all protein-coding human genes of 20,687 [23], (2) Ziats and Rennert [13] based upon a total pool of 656 pathways and 90 bio-functions and disease categories in the IPA Knowledge Database, and (3) Werling et al. [16], an adult brain (102 mRNA/14,246), a replication adult brain data set (50 mRNA/14,869), and a prenatal brain (303 mRNA/9,840) (P < 0.005, FC > |1.2|).

Results

A total of 17 STS and 17 PAC samples were analyzed from 10 ASD and 8 control subjects (Table 1). There were no significant differences in average ages or postmortem interval (PMI) between groups (Table 1). All data and analyses are available in supplementary materials or upon request.

Sexual dimorphism of sncRNA in ASD brain

There are 20 sncRNAs differentially expressed (DE) in STS of ASD females compared to control females, which is significantly more (P = 0.04) than the 8 sncRNAs regulated in STS of ASD males compared to control males (Fig. 1, Table 2, Additional file 1: Table S1). There are eight sncRNAs DE in PAC of ASD females compared to control females, which was not significantly more than the three sncRNAs (P = 0.23) regulated in PAC of ASD male compared to control males (Fig. 1, Table 2, Additional file 1: Table S2). There are significantly more dysregulated sncRNAs in STS compared to PAC in ASD females, but not in ASD males (Fig. 1, Table 2). Regional analyses for STS vs PAC are shown in Fig. 2. There are 55 combined dysregulated sncRNAs in ASD females and 34 in ASD males (P = 0.03) (Fig. 2, Table 2, Additional file 1: Tables S3–S6). One overlapping sncRNA (mir-455) is DE in an opposite direction in ASD females and ASD males compared to control females and control males, respectively.

Table 2 Dysregulated sncRNA in the temporal cortex of ASD male and female brains
Fig. 2
figure 2

Regional sncRNA dysregulation in the superior temporal sulcus (STS) vs the primary auditory cortex (PAC) for female ASD and female control subjects and for male ASD and male control subjects

Predicted mRNA targets and pathways

Pathway analyses are summarized in Additional file 1: Table S7 and presented in Additional file 1: Tables S8–S17. There are significantly more over-represented pathways in ASD females than those in ASD males (P < 0.0001) in regional analysis (STS vs PAC), and the numbers of common pathways are also significant (P < 0.05) (Fig. 3). Of the total over-represented pathways in regional analyses, ASD females have significantly more immune, neurotransmitter, and other nervous system signaling pathways than those in ASD males (P < 0.0001) (Fig. 3). It is notable that even for common pathways, the specific target mRNAs are often sexually dimorphic. For example, both male and female target genes significantly populate the axonal guidance signaling pathway, yet less than a third of the targets are common to both sexes (Fig. 4).

Fig. 3
figure 3

Predicted mRNA target genes and their enriched pathways for the sncRNA dysregulated in the regional STS to PAC comparison for female ASD and female control subjects combined (pink) compared to male ASD and male controls subjects combined (blue). The predicted targets and enriched pathways that are common to males and females are in green

Fig. 4
figure 4

Sexual dimorphism in the axonal signaling pathway. Predicted targets of sexually dimorphic mature miRNAs that were dysregulated in regional analysis comparison of STS to PAC. Pink indicates targets predicted to be regionally dysregulated in ASD female only, blue are dysregulated in ASD males only, and green regulated in both males and females. Note that the majority of predicted targets are sexually dimorphic

Overlap of the sexually dimorphic pathways with others implicated in ASD

We found a significant overlap between the predicted targets and pathways with the findings of Ziats and Rennert [13, 14] and SFARI risk genes (Additional file 2: Figure S1 and Additional file Tables 1: S18-S22). Seventy-four genes overlap between SFARI genes and our female regional analysis (STS vs PAC) and are significantly enriched in nervous system signaling, including glutamate receptor signaling, reelin signaling in neurons, and CREB signaling in neurons, as well as in p53 signaling. There is no significant overlap between our results and the three groups of Werling et al. [16]. However, we performed IPA pathway analysis on their 433 combined sexually dimorphic mRNA, and they are over-represented in 14 pathways. Of these, three overlap with our PAC pathways and nine overlap with our STS vs PAC pathways (trend towards significance, P < 0.1). In addition, miR-125 is sexually dimorphically dysregulated in the prefrontal cortex in Ziats and Rennert [14] and in the STS in our study (Table 2).

Discussion

This short report underscores the importance of considering biological sex as a factor when interpreting gene expression studies of postmortem brains from individuals with ASD. Our findings from this small sample suggest that the expression of miRNA and other sncRNA is sexually dimorphic in the temporal cortex of ASD individuals. There are generally more dysregulated sncRNA, miRNA target genes, and pathways in ASD females compared to ASD males. There is a significant overlap between the male miRNA putatively dysregulated pathways in our study and that which occurs during normal male neural development [13]. Though there is an overlap between the female ASD regional STS-PAC miRNA targets and the SFARI genes, most of the sncRNA-regulated sexually dimorphic target genes are not enriched for autism risk genes, similar to the findings of Werling et al. [16]. The sexually dimorphic over-enrichment of miRNA target genes in the immune and nervous system pathways is consistent with prior gene expression studies of ASD brain [24, 25] that postulated that the immune pathways were related to environment and the neuronal pathways were related to genetics [24].

The greater number of dysregulated sncRNA, miRNA target genes, and pathways in females compared to male ASD subjects supports a body of evidence, suggesting that there is a greater genetic load in ASD females [9] [3, 26], and is consistent with a possible female protective effect [3, 27, 28]. The greater sncRNA dysregulation in females in our study might also support a recent proposal that female brains are less vulnerable to ASD because they are more plastic [8]. Our data could be interpreted to mean that the female ASD brain mounts a greater protective molecular response compared to males, factors that may contribute to the so-called “female camouflage effect” [29].

There are a number of sexually dimorphic miRNAs from our study that have been associated with ASD-relevant diseases and processes (Table 3). For example, one study of miRNA in serum of children with ASD [30] found two miRNAs, miR-151 and miR-181 that were also differentially expressed in the current study. miR-181 promotes synaptogenesis and decreases in axon growth [31, 32]. It is expressed in the brain and previously associated with autistic phenotypes [33] and schizophrenia [34]. The microRNA is also involved in an inflammatory response [35], influences apoptosis and mitochondrial function [36] in astrocytes, and targets GABA receptors [37].

Table 3 Relevant literature on example miRNAs with sexual dimorphism in ASD relative

In STS of ASD female compared to control female, miR-219 and miR-338 showed the highest level of downregulation. Both miRNAs are involved in regulating oligodendrocyte differentiation and likely myelin production [38]. Though abnormalities of white matter have been observed in ASD brains for some time, recent studies point to important sex differences. For example, a recent anatomical study of autism brains showed large white matter regions showing significant sex × diagnosis interactions [5]. This was supported by sex differences found in the corpus callosum of young children with ASD [6]. While no alterations were observed in ASD males compared to control males, mean diffusivity, axial diffusivity, and radial diffusivity were all increased in ASD females compared to control females [6].

miR-219 and miR-338 both promote oligodendrocyte differentiation [38]. Inhibition of both miRNAs inhibits oligodendrocyte maturation and function in part by directly repressing negative regulators of oligodendrocyte differentiation, including transcription factors Sox6 and Hes5 [38]. miR-338 also attenuates cortical neuronal outgrowth by modulating expression of axon guidance genes and axonal mitochondrial genes [3941]. However, it is the role of miR-219 and miR-338 on differentiation of oligodendrocyte precursors during development (with downregulation over fourfold in STS of ASD females, but not of ASD males) that may significantly contribute to the sexually dimorphic changes of white matter tracts including the corpus callosum seen in ASD brain [5, 6]. It is notable that decreased miR-219 and miR-338 were detected in female STS of ASD brains and not PAC, an important finding since STS is an association cortex implicated in social behavior, which is a core symptom of ASD, whereas PAC is not generally associated with ASD [42]. It will be important to examine miR-219 and miR-338 in amygdala and other brain regions implicated in ASD core symptoms, particularly since the number of oligodendrocytes in amygdala is decreased in ASD brain [43]. It will also be important to quantify oligodendrocyte numbers in STS of ASD brain given the current miRNA results.

It is important to note that sex differences in the miRNAome are prominent in health and disease [4447]. Several miRNAs we identified in this study are also sexually dimorphic in other tissues and disease conditions. For example, miR-100, miR-196a, and miR-31 are sexually dimorphic in human amniotic mesenchymal stem cells (hA-MSCs) from obese versus normal weight women who gave birth to females, but not males [48]. This suggests that these miRNAs may be involved in metabolic changes. Additionally, miR-31 was sexually dimorphic in an animal model of systemic lupus erythematosus (SLE) [49]. Several mature miRNAs map to cytoband 14q32, where there is an imprinted miRNA cluster. In our study, miRNA-151b was downregulated in ASD-STS female vs typical STS female, while miR-342-3p, miR-432, and miR-485-5p were upregulated in ASD male in the inter-regional analysis compared to typical male. Loss of imprinting in this region has been associated with multiple diseases [50, 51]. Our results suggest involvement of miRNAs within or near this cluster in ASD as well.

As we have reported previously [17], there are changes of expression of sncRNA in ASD-PAC, though there are more changes in ASD-STS. In fact, some of the miRNAs we previously reported are included in the present study, indicating they may have been driven by sex differences as well. Interestingly, even though there were attenuated regional differences (STS vs PAC) in both male and female ASD compared to control (Fig. 2), the greatest number of miRNA targets by far was identified in the regional comparisons of STS to PAC in both females and males, though there were significantly more targets in females. These large differences of expression between regions, particularly between primary sensory cortices like PAC and association cortices like STS, have been noted previously in neurotypical adult and developing human brain [52, 53] and are consistent with the present data both for ASD and for controls. Integrative functional analyses of ASD risk genes implicate cortical layers II/III and V/VI pyramidal neurons [54], and co-expression network analyses of ASD risk genes also implicate layers V/VI cortical projection neurons [55]. The current findings of dysregulated miRNA that control oligodendrocyte differentiation in female STS, but not in male STS, adds complexity to the picture particularly since miR-219 and miR-338 are not dysregulated in either females or males in PAC. Thus, if projection neurons are involved in ASD, they are likely selectively affected in specific cortical regions and it is possible that there are regional differences of miRNA-mediated oligodendrocyte-myelination of the projection neuron axons that are also sexually dimorphic.

In general, there is little overlap of the specific miRNA regulated in ASD and control brains in this study and the miRNA described by Ziats and Rennert for typically developing human brain [14]. However, there is tremendous overlap of the predicted pathways that are sexually dimorphic. There is significant overlap for the pathways for males in our regional analyses and male pathways in neurotypical brains [13] and a trend for overlap between the sexually dimorphic pathways reported by Werling et al. [16] and our study.

There are important limitations to consider when interpreting this data that are common to studies of postmortem human brains, including small sample size and variation in cause of death, postmortem interval, age at death, agonal state, postmortem RNA integrity, and tissue preparation at different brain banks. In addition, the changes in miRNA and predicted targets and pathways studied here could lead to some aspects of ASD, but they could just as well be a consequence of the condition and even represent compensatory mechanisms. Future studies will need to assess both miRNA and mRNA in the same samples, so that one can determine if specific miRNA target mRNA are present and are regulated in an inverse direction as occurs for most miRNA-mRNA interactions. Small sample size is also limitation of this report and findings need to be confirmed in larger cohorts and with alternative gene expression profiling techniques such as qRT-PCR. However, this study clearly indicates that sex needs to be considered when interpreting data on postmortem human brains in ASD.

Conclusions

Our own molecular studies, and many others, have demonstrated clear age-related changes in the brains of both typically developing and autistic individuals [17, 18, 24, 56]. Beyond age considerations, the current data also support a very important role for biological sex and suggest that pharmacological treatments would likely have to be evaluated separately in both sexes [3, 13, 15, 16]. A greater number of dysregulated sncRNAs and their gene targets and pathways in female ASD brain are consistent with a greater genetic load in females, a female protective effect, and possibly greater plasticity of female ASD brain. Although speculative, there are specific miRNAs dysregulated in STS of female ASD brain associated with oligodendrocyte differentiation (miR-219 and miR-338) that could relate to sexual dimorphism of white matter tracts, miR-488 that could relate to more anxiety in females, and miR-125 and miR-181 implicated in neuronal development that may be sexually dimorphic.

Abbreviations

ASD:

Autism spectrum disorders

BA:

Brodmann’s Area

CNV:

Copy number variant

F:

Female

M:

Male

miRNA or miR-:

microRNA

PAC:

Primary auditory cortex of temporal lobe

PMI:

Postmortem interval

RNAseq:

RNA sequencing

sncRNA:

Small noncoding RNA

SNV:

Single nucleotide variant

STS:

Superior temporal sulcus of temporal lobe

References

  1. Lai MC, Baron-Cohen S, Buxbaum JD. Understanding autism in the light of sex/gender. Mol Autism. 2015;6:24.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Tylee DS, Espinoza AJ, Hess JL, Tahir MA, McCoy SY, Rim JK, Dhimal T, Cohen OS, Glatt SJ. RNA sequencing of transformed lymphoblastoid cells from siblings discordant for autism spectrum disorders reveals transcriptomic and functional alterations: evidence for sex-specific effects. Autism Res. 2016.

  3. Werling DM, Geschwind DH. Sex differences in autism spectrum disorders. Curr Opin Neurol. 2013;26(2):146–53.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Campbell DJ, Chang J, Chawarska K. Early generalized overgrowth in autism spectrum disorder: prevalence rates, gender effects, and clinical outcomes. J Am Acad Child Adolesc Psychiatry. 2014;53(10):1063–73. e1065.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Lai MC, Lombardo MV, Suckling J, Ruigrok AN, Chakrabarti B, Ecker C, Deoni SC, Craig MC, Murphy DG, Bullmore ET, et al. Biological sex affects the neurobiology of autism. Brain. 2013;136(Pt 9):2799–815.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Nordahl CW, Iosif AM, Young GS, Perry LM, Dougherty R, Lee A, Li D, Buonocore MH, Simon T, Rogers S, et al. Sex differences in the corpus callosum in preschool-aged children with autism spectrum disorder. Mol Autism. 2015;6:26.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Schumann CM, Bloss CS, Barnes CC, Wideman GM, Carper RA, Akshoomoff N, Pierce K, Hagler D, Schork N, Lord C, et al. Longitudinal magnetic resonance imaging study of cortical development through early childhood in autism. J Neurosci. 2010;30(12):4419–27.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Mottron L, Duret P, Mueller S, Moore RD, Forgeot d’Arc B, Jacquemont S, Xiong L. Sex differences in brain plasticity: a new hypothesis for sex ratio bias in autism. Mol Autism. 2015;6:33.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Robinson EB, Lichtenstein P, Anckarsater H, Happe F, Ronald A. Examining and interpreting the female protective effect against autistic behavior. Proc Natl Acad Sci U S A. 2013;110(13):5258–62.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Iossifov I, O’Roak BJ, Sanders SJ, Ronemus M, Krumm N, Levy D, Stessman HA, Witherspoon KT, Vives L, Patterson KE, et al. The contribution of de novo coding mutations to autism spectrum disorder. Nature. 2014;515(7526):216–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Neale BM, Kou Y, Liu L, Ma’ayan A, Samocha KE, Sabo A, Lin CF, Stevens C, Wang LS, Makarov V, et al. Patterns and rates of exonic de novo mutations in autism spectrum disorders. Nature. 2012;485(7397):242–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Lombardo MV, Chakrabarti B, Lai MC, Consortium MA, Baron-Cohen S. Self-referential and social cognition in a case of autism and agenesis of the corpus callosum. Mol Autism. 2012;3(1):14.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Ziats MN, Rennert OM. Sex-biased gene expression in the developing brain: implications for autism spectrum disorders. Mol Autism. 2013;4(1):10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Ziats MN, Rennert OM. Identification of differentially expressed microRNAs across the developing human brain. Mol Psychiatry. 2014;19(7):848–52.

    Article  CAS  PubMed  Google Scholar 

  15. Hu VW, Sarachana T, Sherrard RM, Kocher KM. 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. Mol Autism. 2015;6:7.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Werling DM, Parikshak NN, Geschwind DH. Gene expression in human brain implicates sexually dimorphic pathways in autism spectrum disorders. Nat Commun. 2016;7:10717.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Ander BP, Barger N, Stamova B, Sharp FR, Schumann CM. Atypical miRNA expression in temporal cortex associated with dysregulation of immune, cell cycle, and other pathways in autism spectrum disorders. Mol Autism. 2015;6:37.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Stamova B, Ander BP, Barger N, Sharp FR, Schumann CM. Specific regional and age-related small noncoding RNA expression patterns within superior temporal gyrus of typical human brains are less distinct in autism brains. J Child Neurol. 2015;30(14):1930–46.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Alaerts K, Nayar K, Kelly C, Raithel J, Milham MP, Di Martino A. Age-related changes in intrinsic function of the superior temporal sulcus in autism spectrum disorders. Soc Cogn Affect Neurosci. 2015;10(10):1413–23.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Bigler ED, Mortensen S, Neeley ES, Ozonoff S, Krasny L, Johnson M, Lu J, Provencal SL, McMahon W, Lainhart JE. Superior temporal gyrus, language function, and autism. Dev Neuropsychol. 2007;31(2):217–38.

    Article  PubMed  Google Scholar 

  21. Pantelis PC, Byrge L, Tyszka JM, Adolphs R, Kennedy DP. A specific hypoactivation of right temporo-parietal junction/posterior superior temporal sulcus in response to socially awkward situations in autism. Soc Cogn Affect Neurosci. 2015;10(10):1348–56.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Redcay E. The superior temporal sulcus performs a common function for social and speech perception: implications for the emergence of autism. Neurosci Biobehav Rev. 2008;32(1):123–42.

    Article  PubMed  Google Scholar 

  23. ENCODE Project Consortium. An integrated encyclopedia of DNA elements in the human genome. Nature. 2012;489(7414):57–74.

  24. Voineagu I, Wang X, Johnston P, Lowe JK, Tian Y, Horvath S, Mill J, Cantor RM, Blencowe BJ, Geschwind DH. Transcriptomic analysis of autistic brain reveals convergent molecular pathology. Nature. 2011;474(7351):380–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Garbett K, Ebert PJ, Mitchell A, Lintas C, Manzi B, Mirnics K, Persico AM. Immune transcriptome alterations in the temporal cortex of subjects with autism. Neurobiol Dis. 2008;30(3):303–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Jacquemont S, Coe BP, Hersch M, Duyzend MH, Krumm N, Bergmann S, Beckmann JS, Rosenfeld JA, Eichler EE. A higher mutational burden in females supports a “female protective model” in neurodevelopmental disorders. Am J Hum Genet. 2014;94(3):415–25.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Gockley J, Willsey AJ, Dong S, Dougherty JD, Constantino JN, Sanders SJ. The female protective effect in autism spectrum disorder is not mediated by a single genetic locus. Mol Autism. 2015;6:25.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Constantino JN. Data from the Baby Siblings Research Consortium confirm and specify the nature of the female protective effect in autism: a commentary on Messinger et al. Mol Autism. 2016;7:32.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Rynkiewicz A, Lucka I, Grabowski K. Letter to editor. report “An investigation of the ‘female camouflage effect’ in autism using a new computerized test showing sex/gender differences during ADOS-2”. IMFAR 2016, Baltimore, USA. Psychiatr Pol. 2016;50(3):663–6.

    Article  PubMed  Google Scholar 

  30. Mundalil Vasu M, Anitha A, Thanseem I, Suzuki K, Yamada K, Takahashi T, Wakuda T, Iwata K, Tsujii M, Sugiyama T, et al. Serum microRNA profiles in children with autism. Mol Autism. 2014;5:40.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Kos A, Olde Loohuis N, Meinhardt J, Van Bokhoven H, Kaplan BB, Martens GJ, Aschrafi A. MicroRNA-181 promotes synaptogenesis and attenuates axonal outgrowth in cortical neurons. Cell Mol Life Sci. 2016.

  32. Carrella S, Barbato S, D’Agostino Y, Salierno FG, Manfredi A, Banfi S, Conte I. TGF-beta controls miR-181/ERK regulatory network during retinal axon specification and growth. PLoS ONE. 2015;10(12), e0144129.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Stepniak B, Kastner A, Poggi G, Mitjans M, Begemann M, Hartmann A, Van der Auwera S, Sananbenesi F, Krueger-Burg D, Matuszko G, et al. Accumulated common variants in the broader fragile X gene family modulate autistic phenotypes. EMBO Mol Med. 2015;7(12):1565–79.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Rizos E, Siafakas N, Katsantoni E, Skourti E, Salpeas V, Rizos I, Tsoporis JN, Kastania A, Filippopoulou A, Xiros N, et al. Correction: Let-7, Mir-98 and Mir-181 as biomarkers for cancer and schizophrenia. PLoS ONE. 2015;10(8), e0135863.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Hutchison ER, Kawamoto EM, Taub DD, Lal A, Abdelmohsen K, Zhang Y, Wood 3rd WH, Lehrmann E, Camandola S, Becker KG, et al. Evidence for miR-181 involvement in neuroinflammatory responses of astrocytes. Glia. 2013;61(7):1018–28.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Ouyang YB, Lu Y, Yue S, Giffard RG. miR-181 targets multiple Bcl-2 family members and influences apoptosis and mitochondrial function in astrocytes. Mitochondrion. 2012;12(2):213–9.

    Article  CAS  PubMed  Google Scholar 

  37. Zhao C, Huang C, Weng T, Xiao X, Ma H, Liu L. Computational prediction of MicroRNAs targeting GABA receptors and experimental verification of miR-181, miR-216 and miR-203 targets in GABA-A receptor. BMC Res Notes. 2012;5:91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Zhao X, He X, Han X, Yu Y, Ye F, Chen Y, Hoang T, Xu X, Mi QS, Xin M, et al. MicroRNA-mediated control of oligodendrocyte differentiation. Neuron. 2010;65(5):612–26.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Kos A, Klein-Gunnewiek T, Meinhardt J, Loohuis NF, Van Bokhoven H, Kaplan BB, Martens GJ, Kolk SM, Aschrafi A. MicroRNA-338 attenuates cortical neuronal outgrowth by modulating the expression of axon guidance genes. Mol Neurobiol. 2016.

  40. Vargas JN, Kar AN, Kowalak JA, Gale JR, Aschrafi A, Chen CY, Gioio AE, Kaplan BB. Axonal localization and mitochondrial association of precursor microRNA 338. Cell Mol Life Sci. 2016.

  41. Aschrafi A, Kar AN, Natera-Naranjo O, MacGibeny MA, Gioio AE, Kaplan BB. MicroRNA-338 regulates the axonal expression of multiple nuclear-encoded mitochondrial mRNAs encoding subunits of the oxidative phosphorylation machinery. Cell Mol Life Sci. 2012;69(23):4017–27.

    Article  CAS  PubMed  Google Scholar 

  42. Saitovitch A, Bargiacchi A, Chabane N, Brunelle F, Samson Y, Boddaert N, Zilbovicius M. Social cognition and the superior temporal sulcus: implications in autism. Rev Neurol (Paris). 2012;168(10):762–70.

    Article  CAS  Google Scholar 

  43. Morgan JT, Barger N, Amaral DG, Schumann CM. Stereological study of amygdala glial populations in adolescents and adults with autism spectrum disorder. PLoS ONE. 2014;9(10), e110356.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Flanagan KL. Sexual dimorphism in biomedical research: a call to analyse by sex. Trans R Soc Trop Med Hyg. 2014;108(7):385–7.

    Article  PubMed  Google Scholar 

  45. Guebel DV, Torres NV. Sexual dimorphism and aging in the human hyppocampus: identification, validation, and impact of differentially expressed genes by factorial microarray and network analysis. Front Aging Neurosci. 2016;8:229.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Pak TR, Rao YS, Prins SA, Mott NN. An emerging role for microRNAs in sexually dimorphic neurobiological systems. Pflugers Arch. 2013;465(5):655–67.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Sharma A, Hoeffer CA, Takayasu Y, Miyawaki T, McBride SM, Klann E, Zukin RS. Dysregulation of mTOR signaling in fragile X syndrome. J Neurosci. 2010;30(2):694–702.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Nardelli C, Granata I, Iaffaldano L, D’Argenio V, Del Monaco V, Maruotti GM, Del Vecchio L, Martinelli P, Salvatore F, Guarracino MR, et al. Sex-comparative analysis of the miRNome of human amniotic stem cells during obesity. Stem Cells Dev. 2016.

  49. Dai R, McReynolds S, Leroith T, Heid B, Liang Z, Ahmed SA. Sex differences in theexpression of lupus-associated miRNAs in splenocytes from lupus-prone NZB/WF1 mice. Biol Sex Differ. 2013;4(1):19.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Benetatos L, Hatzimichael E, Londin E, Vartholomatos G, Loher P, Rigoutsos I, Briasoulis E. The microRNAs within the DLK1-DIO3 genomic region: involvement in disease pathogenesis. Cell Mol Life Sci. 2013;70(5):795–814.

    Article  CAS  PubMed  Google Scholar 

  51. Glazov EA, McWilliam S, Barris WC, Dalrymple BP. Origin, evolution, and biological role of miRNA cluster in DLK-DIO3 genomic region in placental mammals. Mol Biol Evol. 2008;25(5):939–48.

    Article  CAS  PubMed  Google Scholar 

  52. Hawrylycz MJ, Lein ES, Guillozet-Bongaarts AL, Shen EH, Ng L, Miller JA, van de Lagemaat LN, Smith KA, Ebbert A, Riley ZL, et al. An anatomically comprehensive atlas of the adult human brain transcriptome. Nature. 2012;489(7416):391–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Miller JA, Ding SL, Sunkin SM, Smith KA, Ng L, Szafer A, Ebbert A, Riley ZL, Royall JJ, Aiona K, et al. Transcriptional landscape of the prenatal human brain. Nature. 2014;508(7495):199–206.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Parikshak NN, Luo R, Zhang A, Won H, Lowe JK, Chandran V, Horvath S, Geschwind DH. Integrative functional genomic analyses implicate specific molecular pathways and circuits in autism. Cell. 2013;155(5):1008–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Willsey AJ, Sanders SJ, Li M, Dong S, Tebbenkamp AT, Muhle RA, Reilly SK, Lin L, Fertuzinhos S, Miller JA, et al. Coexpression networks implicate human midfetal deep cortical projection neurons in the pathogenesis of autism. Cell. 2013;155(5):997–1007.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Chow ML, Pramparo T, Winn ME, Barnes CC, Li HR, Weiss L, Fan JB, Murray S, April C, Belinson H, et al. Age-dependent brain gene expression and copy number anomalies in autism suggest distinct pathological processes at young versus mature ages. PLoS Genet. 2012;8(3), e1002592.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Hudish LI, Galati DF, Ravanelli AM, Pearson CG, Huang P, Appel B. miR-219 regulates neural progenitors by dampening apical Par protein-dependent Hedgehog signaling. Development. 2016;143(13):2292–304.

    Article  CAS  PubMed  Google Scholar 

  58. Zhang Y, Fan M, Wang Q, He G, Fu Y, Li H, Yu S. Polymorphisms in microRNA genes and genes involving in NMDAR signaling and schizophrenia: a case-control study in Chinese Han population. Sci Rep. 2015;5:12984.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Pusic AD, Kraig RP. Youth and environmental enrichment generate serum exosomes containing miR-219 that promote CNS myelination. Glia. 2014;62(2):284–99.

    Article  PubMed  Google Scholar 

  60. Ebrahimi-Barough S, Kouchesfehani HM, Ai J, Mahmoodinia M, Tavakol S, Massumi M. Programming of human endometrial-derived stromal cells (EnSCs) into pre-oligodendrocyte cells by overexpression of miR-219. Neurosci Lett. 2013;537:65–70.

    Article  CAS  PubMed  Google Scholar 

  61. Ebrahimi-Barough S, Massumi M, Kouchesfahani HM, Ai J. Derivation of pre-oligodendrocytes from human endometrial stromal cells by using overexpression of microRNA 338. J Mol Neurosci. 2013;51(2):337–43.

    Article  CAS  PubMed  Google Scholar 

  62. Muinos-Gimeno M, Espinosa-Parrilla Y, Guidi M, Kagerbauer B, Sipila T, Maron E, Pettai K, Kananen L, Navines R, Martin-Santos R, et al. Human microRNAs miR-22, miR-138-2, miR-148a, and miR-488 are associated with panic disorder and regulate several anxiety candidate genes and related pathways. Biol Psychiatry. 2011;69(6):526–33.

    Article  CAS  PubMed  Google Scholar 

  63. Stappert L, Borghese L, Roese-Koerner B, Weinhold S, Koch P, Terstegge S, Uhrberg M, Wernet P, Brustle O. MicroRNA-based promotion of human neuronal differentiation and subtype specification. PLoS ONE. 2013;8(3), e59011.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors acknowledge the technical assistance of Alicja Omanska and Dr. Nicole Barger of the UC Davis MIND Institute, as well as Ryan Davis and Stephanie Liu of the Genomics Shared Resource at UC Davis. The authors also express their gratitude to the brain donors and their families who made this research possible.

Funding

R21 MH108909-01 (CMS, BS), the UC Davis Children’s Miracle Network (FRS, BS), R01 MH 097236 (CMS), and NICHD U54 HD079125 (MIND Institute).

Availability of data and materials

Data analyzed during this study are included in this published article (and its additional files) and/or available from the corresponding author on reasonable request.

Authors’ contributions

BPA and CMS performed the experiments. CMS, BPA, FRS, and BS conceived and designed the studies. BS performed the data analyses. CMS, FRS, and BS wrote the manuscript. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Consent was not required to participate. This study was deemed exempt by the Internal Review Board at the UC Davis School of Medicine.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Boryana Stamova.

Additional files

Additional file 1:

Supplementary Data Tables. Table S1. sncRNA differentially expressed in ASD-STS vs CTRL-STS. Table S2. sncRNA differentially expressed in ASD-PAC vs CTRL-PAC. Table S3. Female (F) ASD Regional Analysis. F_ASD-STS vs F_ASD-PAC. Table S4. Female (F) controls (CTRL) regional analysis. F_CTRL-STS vs F_CTRL-PAC. Table S5. Male (M) ASD regional analysis. M_ASD-STS vs M_ASD-PAC. Table S6. Male (M) controls (CTRL) regional analysis. M_CTRL-STS vs M_CTRL-PAC. Table S7. Summary of putative target analyses. Table S8. Pathway analysis of putative targets of the mature miRNA from the PAC analysis in females. Table S9. Pathway analysis of putative targets of the mature miRNA from the PAC analysis in males. Table S10. Pathway analysis of putative targets of the mature miRNA from the STS analysis in females. Table S11. 227 common targets between female and male regional analyses. Table S12. Pathway analysis of putative targets of the mature miRNA from the female regional analysis. Table S13. Pathway analysis of putative targets of the mature miRNA from the male regional analysis. Table S14. Pathway analysis of putative targets of the mature miRNA from the female regional analysis: immune pathways. Table S15. Pathway analysis of putative targets of the mature miRNA from the male regional analysis: immune pathways. Table S16. Pathway analysis of putative targets of the mature miRNA from the female regional analysis: Neurotransmitter and other nervous system signaling pathways. Table S17. Pathway analysis of putative targets of the mature miRNA from the male regional analysis: neurotransmitter and other nervous system signaling pathways. Table S18. Pathway overlapping between our male regional analysis and the male-enriched pathways from Ziats and Rennert (Ref 15). Table S19. Seventy-four ASD-implicated genes amongst the 1382 putative targets of the miRNAs from the female regional analysis. Table S20. Canonical pathway enriched in the 74 ASD-implicated genes from the predicted targets in the female regional analysis Table S21. Targets overlapping between our analysis and the sexually dimorphic mRNA from Werling et al, Nature Communications, 2016; Supplementary Material: ncomms10717-s2. Table S22. Pathway overlapping between our male and female regional analysis and pathways were derived based on combined sexually dimorphic mRNAs from Werling et al, Nature Communications, 2016; Supplementary Material: ncomms10717-s2. (XLSX 130 kb)

Additional file 2:

Supplementary Figure. Overlap with relevant studies. A Overlap between our male regional analysis (118 pathways) and Ziats and Rennert, Mol. Autism, 2013, male-enriched pathways (10 pathways); P of overlap = 0.0034. B Overlap between all our sexually dimorphic pathways (290 pathways) and Ziats and Rennert, Mol. Psychiatry, 2014, sexually dimorphic pathways (204 pathways); P of overlap < 0.00001. C Overlap between all our sexually dimorphic pathways (287 pathways) and Werling et al, Nat. Communications, 2016, sexually dimorphic pathways (14 pathways). P of overlap between Werling et al, 2016 and our regional analysis = 0.098, P of overlap between Werling et al, 2016, and our PAC analysis = 0.092. D Overlap between the predicted targets from our female regional analysis (STS vs PAC) (1382 targets) and SFARI ASD-implicated genes (768). P of overlap < 0.001. (TIF 48 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Schumann, C.M., Sharp, F.R., Ander, B.P. et al. Possible sexually dimorphic role of miRNA and other sncRNA in ASD brain. Molecular Autism 8, 4 (2017). https://doi.org/10.1186/s13229-017-0117-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13229-017-0117-0

Keywords