Skip to main content

Atypical miRNA expression in temporal cortex associated with dysregulation of immune, cell cycle, and other pathways in autism spectrum disorders

Abstract

Background

Autism spectrum disorders (ASDs) likely involve dysregulation of multiple genes related to brain function and development. Abnormalities in individual regulatory small non-coding RNA (sncRNA), including microRNA (miRNA), could have profound effects upon multiple functional pathways. We assessed whether a brain region associated with core social impairments in ASD, the superior temporal sulcus (STS), would evidence greater transcriptional dysregulation of sncRNA than adjacent, yet functionally distinct, primary auditory cortex (PAC).

Methods

We measured sncRNA expression levels in 34 samples of postmortem brain from STS and PAC to find differentially expressed sncRNA in ASD compared with control cases. For differentially expressed miRNA, we further analyzed their predicted mRNA targets and carried out functional over-representation analysis of KEGG pathways to examine their functional significance and to compare our findings to reported alterations in ASD gene expression.

Results

Two mature miRNAs (miR-4753-5p and miR-1) were differentially expressed in ASD relative to control in STS and four (miR-664-3p, miR-4709-3p, miR-4742-3p, and miR-297) in PAC. In both regions, miRNA were functionally related to various nervous system, cell cycle, and canonical signaling pathways, including PI3K-Akt signaling, previously implicated in ASD. Immune pathways were only disrupted in STS. snoRNA and pre-miRNA were also differentially expressed in ASD brain.

Conclusions

Alterations in sncRNA may underlie dysregulation of molecular pathways implicated in autism. sncRNA transcriptional abnormalities in ASD were apparent in STS and in PAC, a brain region not directly associated with core behavioral impairments. Disruption of miRNA in immune pathways, frequently implicated in ASD, was unique to STS.

Background

The autism spectrum disorders (ASDs) are heritable neurodevelopmental conditions characterized by impairments in social interaction and communication and restricted, repetitive, and stereotyped patterns of behavior. Several hundred different genetic loci, including single gene mutations, chromosome abnormalities, and copy number variations, have now been associated with a growing number of diagnosed cases [1]. The genetic underpinnings of the remainder of the cases are unclear, but they are postulated as being related to the interactions of multiple genes with environmental factors.

One of the questions central to the study of ASD is how so many different genetic conditions can produce a similar clinical phenotype. Part of the answer to this problem appears to be that a variety of etiologically heterogeneous disorders converge on and disrupt key stages of neurodevelopment that affect specific brain regions and molecular pathways. Evidence for this has come in part from pathway analyses of ASD candidate genes and changes of gene expression in brain tissue [2, 3]. Several genomic and gene expression studies of brain tissue have shown distinct and reproducible changes in transcriptome organization [4], genetic pathways related to neuron function and development [57], and immune pathways [5] in ASD.

Given that changes in gene expression must ultimately underlie the clinical phenotype at least to some degree, what might be driving these changes that appear to be region and perhaps even cell specific? Among the many factors that can regulate gene expression—including DNA methylation, transcription factors, and others—we opted to examine small non-coding RNA (sncRNA) in this study, including microRNA (miRNA) and small nucleolar RNA (snoRNA). miRNA are of particular interest because, although they do not code for protein, they bind sequence-specific sites in target transcripts to regulate expression levels of mRNA and/or modulate protein translation [8]. There is growing evidence suggesting miRNA may play a role in ASD supported through mining of SNP and CNV data [9, 10], epidemiological studies of genetic variants [11], and serum biomarker studies [12, 13]. It is estimated that over 60 % of all genes are regulated by miRNA because each miRNA has dozens to hundreds of targets. For this reason, we postulated that miRNA could play an important role in regulating the transcriptional changes that have been reported in ASD brain. Moreover, because miRNA (and other small and non-coding RNA) are often cell or tissue specific and are important in timed molecular events, they can play critical roles in very specific aspects of neurodevelopment [8]. Disturbances in a single miRNA could affect hundreds of genes, directly, and many hundreds of genes, indirectly, and would be one way where a single disturbance could dramatically affect complex cellular organization and function.

Recent whole genome expression studies show that gene expression in brain tissue is region specific and can reflect the distributions of major cell types including neurons and glia [14]. In cases with no known developmental disorders, neocortical transcriptional patterns are relatively homogenous and physically closer regions typically exhibit more similar transcriptomes. However, the primary sensory cortices appear to show the most distinct transcriptional profiles [14], suggesting that they are transcriptionally as well as functionally distinct from other cortical subtypes. To begin to address the underlying genetics of neuropathology, examination of transcriptionally distinct brain regions may be crucial in developing a precise understanding of developmental or pathological disturbances in regulatory mRNA modulators like miRNA.

To examine the expression of the small non-coding RNAs in ASD versus typically developing brains, we focused on two functionally distinct but physically close areas in the temporal lobe, the superior temporal sulcus (STS) and the primary auditory cortex (PAC). The STS is instrumental in functions critical for reciprocal social interaction [1517] and language [1820] that are intimately related to core features of ASD [15]. Abnormalities in the STS have been detected in ASD including decreased gray matter, hypo-perfusion at resting state, and abnormal activation during social tasks [21, 22]. In contrast, the PAC is a primary sensory cortex that processes rudimentary components of sounds, e.g., auditory tone, and may function similarly in ASD and control subjects [23]. As such, we predicted that miRNA in PAC would be less affected in ASD than STS. Through careful dissection of these two regions, isolation and examination of expressed sncRNA, identification of regulated gene targets, and functional interpretation of RNA, we found, contrary to our expectation, miRNA in both regions were altered, while no single miRNA was commonly affected in both regions.

Methods

Brain tissue

This study was exempt from human subjects review by the Internal Review Board at the UC Davis School of Medicine. Donor brain tissue was obtained from the Autism Tissue Program collection previously housed at the Harvard Brain Tissue Resource Center (http://www.mcleanhospital.org/research-programs/harvard-brain-tissue-resource-center). Postmortem Confirmation of Consent was provided by next-of-kin and held along with identifiable personal health information by the Autism Tissue Program. A total of 36 samples were obtained from ten ASD and eight control subjects. Individual and summarized subject demographics are presented in Table 1. One STS and one PAC sample were taken from the same frozen coronal brain section of each donor brain. The primary auditory cortex (PAC) sample included the primary sensory areas of Brodmann’s areas 41 and 42 and was taken from the crown of Heschl’s gyrus (Fig. 1), where it is reliably found as indicated by functional and cytoarchitectonic maps [24]. Samples from the superior temporal sulcus (STS) included the polymodal association cortex of Brodmann’s area 22 and were taken from the upper wall of the superior temporal sulcus immediately opposite Heschl’s gyrus [25, 26] (Fig. 1).

Table 1 Individual demographic information for the brain donors of samples used in this study and postmortem interval of samples
Fig. 1
figure 1

Sampling of brain regions. Schematic of the human brain indicating regions of the superior temporal gyrus sampled—superior temporal sulcus (STS - orange) and primary auditory cortex (PAC - teal). A section from the caudal portion of the temporal lobe was utilized for each case (the approximate position is noted by the dashed lines). In the same coronal section, tissue blocks for the PAC (dark teal) were sampled from Heschl’s gyrus (light teal) and blocks for the STS (dark orange) were sampled from the upper wall of the STS (light orange)

RNA isolation and arrays

Total RNA was isolated from ~10 to 75 mg of dissected brain tissue using the RecoverAll Total Nucleic Acid Isolation Kit (Ambion) with omission of the heat de-paraffinization step. Yields of RNA were an average of 200 ng/mg brain tissue. Quality of RNA was assessed using the Agilent 2100 Bioanalyzer with the RNA 6000 Nano kit, and RNA quantity was measured using the NanoDrop ND-1000. Samples had A260/A280 absorbance ratios ≥1.6 [27, 28], and RNA integrity numbers ranged between 2.2 and 5.3.

Two hundred nanograms of total RNA were processed for Affymetrix miRNA 3.0 microarrays according to the manufacturer’s protocol and our previous publication [28]. These arrays contain probes for 5607 small non-coding RNA including 1733 mature human miRNA, 1658 human pre-miRNA, and 2216 human snoRNA, CDBox RNA, H/ACA Box RNA, and scaRNA. Briefly, 200 ng of RNA is ligated to a biotinylated marker using the FlashTag Biotin HSR RNA labeling kit (Affymetrix). No amplification of the RNA was required during this procedure. The samples are hybridized to the microarrays, and biotin is detected with streptavidin-PE. The microarrays are scanned and fluorescent signals measured and transferred into a CEL file for subsequent analyses. Standard quality control measures are applied, such as excluding arrays that are significant outliers or showing skewing in the overall intensity images. Two samples were excluded after these tests (see Table 1). All within chip quality control hybridization metrics on remaining samples were within normal ranges indicating successful sample processing.

Statistics

Analyses were performed using Partek Genomics Suite 6.6 (Partek Inc., St. Louis). Cel files were imported and normalized using Robust Multiarray Averaging (RMA) [29]. The primary comparisons were expression differences of each sncRNA for ASD versus control cases in the PAC and expression differences of each sncRNA for ASD versus control in the STS. Age and postmortem interval (PMI) were included in the following statistical model: Y ijk = μ + Primary Diagnosis i + age + PMI + Region j + Primary Diagnosis*Region ij + ε ijk , where Y ijk represents the k th observation on the i th Primary Diagnosis j th Region, μ is the common effect for the whole experiment, and ε ijk represents the random error component.

Statistical significance was defined as P ≤ 0.005 combined with an absolute fold change ≥1.2. Inclusion of a fold change threshold helps identify changes of biological significance [27, 30, 31].

Unsupervised hierarchical clustering

Unsupervised hierarchical clustering was performed in the Partek Genomics Suite in order to visualize differences of sncRNA expression between ASD and controls. Hierarchical clustering grouped similar sncRNA expression levels and samples into “clusters” using Euclidean dissimilarity and average linkage method and drawing a tree (dendrogram) that shows the hierarchy of the clusters.

Pathway analysis

In order to examine the functional significance of the altered miRNA, computationally derived targets of the miRNA were determined. These predicted mRNA targets were then included in an over-representation pathway analysis to determine biological functions from the Kyoto Encyclopedia of Genes and Genomes (KEGG). This was done using the web-based software tool DIANA miRPath v2.0 [32]. Target genes of miRNA lists were computationally predicted using the DIANA microT-CDS with threshold set to 0.9. These mRNA targets were then used to determine if there were more represented genes in a given biological pathway than would be expected by chance using a Fisher exact test corrected for multiple comparisons (FDR, P < 0.05) [32].

In order to contrast the biological and functional significance of our study with studies previously reporting gene or gene expression lists associated with autism spectrum disorder, we performed functional over-representation analysis on KEGG pathways using reported data from two prior analyses of temporal cortex in ASD, Voineagu et al. (135 probes) [4] and Garbett et al. (152 probe sets) [33], and two datasets broadly focusing on ASD-associated genes, Pinto et al. (139 genes) [34] and the SFARI database of positively associated ASD genes (573 genes) [35]. Datasets were input into the online WebGestalt tool [36] to derive the KEGG pathway data for Fig. 5. Predicted gene targets from our study were taken from microT-CDS after combining all brain samples (STS and PAC) and then deriving miRNA that were different between ASD and CTRL groups (P < 0.005, fold change > |1.2|). Combining temporal regions allowed for a more appropriate comparison to the previous studies that did not target specific functional territories. Euler diagrams of area proportional overlap were generated using eulerAPE [37].

Results

Small non-coding RNA significantly differed between ASD and controls

There were no significant differences in the sex or age of ASD compared to control subjects (P > 0.05, Table 1). Table 2 lists the sncRNA that were differentially expressed in ASD compared to controls (P < 0.005 and fold change > |1.2|) in the temporal lobe for the superior temporal sulcus (STS) and for the primary auditory cortex (PAC) (detailed fold change and P values provided in Additional file 1: Tables S1 and S2). When subjected to hierarchical clustering, these sncRNA completely separated ASD from control subjects for both the PAC and STS brain regions (Fig. 2), lending further support that these sncRNA were significantly different in ASD compared to controls in both brain regions.

Table 2 Differentially expressed sncRNA between ASD and CTRL
Fig. 2
figure 2

Relative expression of sncRNA in STS and PAC of ASD and CTRL subjects. Heat maps showing hierarchical clustering of dysregulated small non-coding RNAs on the X axis for ASD and control (CTRL) subjects on the Y axis for the superior temporal sulcus (STS) (upper panel) and the primary auditory cortex (PAC) (lower panel). Note that the ASD and CTRL subjects cluster together in both the PAC and STS. Bright red color indicates a relative twofold increase in expression for the particular sncRNA, and a bright green color indicates a relative twofold decrease in expression across ASD and control subjects (see bar at bottom of lower panel)

For the STS, there were 14 sncRNA differentially expressed between ASD and control brains (Table 2). Of the two mature miRNA, miR-4753-5p increased expression and miR-1 decreased expression (Fig. 3, Table 2). For the PAC, there were 16 sncRNA differentially expressed between ASD and control brains (Table 2). Of the four mature miRNA, miR-664-3p and miR-4709-3p increased expression and miR-297 and miR-4742-3p decreased expression (Fig. 3, Table 2).

Fig. 3
figure 3

Expression levels of significant miRNA is STS and PAC. Expression levels of microRNA (miRNA) that had significantly increased or decreased expression in ASD compared to control subjects in the superior temporal sulcus and in the primary auditory cortex (P < 0.005 and fold change > |1.2|). ASD values are in red, and control (CTRL) values are in blue. Note that individual subjects are plotted along with box and whisker plots, and the size of the diamonds for each individual corresponds to age of the subject

Small nucleolar RNA (snoRNA), another class of small non-coding RNA, were also differentially expressed between ASD and control brains. ACA39, SNORA11C, and SNORA27 increased in STS while the U13 paralog on chromosome 1 decreased expression in STS (Table 2). In PAC, no snoRNA had increased expression, whereas SNORA22, SCARNA6, and three U13 paralogs on chromosomes 2, 3, and 11 showed decreased expression in ASD compared to controls (Table 2). Eight pre-miRNA changed expression in the STS, and seven pre-miRNA changed expression in the PAC in ASD compared to controls (Table 2).

Functional roles of miRNA in the STS and PAC

Based on the microT-CDS algorithm (which includes coding and 3′ UTR sequences, threshold = 0.9), we found miRNA differentially expressed in ASD had 1214 predicted gene targets in PAC and 445 in STS (Additional file 1: Table S3). Over represented KEGG pathways that contained more gene targets than expected by chance (Fisher exact test, FDR corrected P < 0.05) are summarized in Fig. 4. KEGG pathways regulated by these miRNA showed little overlap between regions (Fig. 4). Indeed, only two pathways were shared between PAC and STS: “pathways in cancer”, which are cell cycle, cell differentiation, mitochondrion, and gene regulation focused, and the PI3K-Akt signaling pathway (Fig. 4). Another major finding is that immune pathways under control of the differentially expressed miRNA appear to be dysregulated only in STS whereas neuronal, cell cycle, signaling and cell processes were dysregulated in both PAC and STS in ASD compared to control subjects (Fig. 4).

Fig. 4
figure 4

Functional categorization of miRNA target genes. KEGG pathways of the predicted gene targets of the microRNA that showed altered expression in ASD compared to control subjects in the PAC and in the STS. These pathways had more predicted regulated target genes in the pathway than would have been expected by chance (P < 0.05). Pathways are grouped according to categories in vertical text. Only miRNA differentially expressed in STS had target genes over-represented in immune related pathways and functions

Figure 5 shows that many of the predicted brain miRNA-regulated pathways in the STS and PAC-combined analysis in this study were also predicted to be dysregulated based on genes identified in the expression studies of Voineagu et al. [4] and Garbett et al. [33] as well as pathways implicated in ASD gene sets from Pinto et al. [34] and the SFARI gene database [35]. The significance of the overlap of our data to each of these data sets was P < 0.05 using the hypergeometric probability function phyper in R.

Fig. 5
figure 5

Functional similarities to other studies looking at ASD-related genes. Overlap of KEGG pathways identified through analysis of mRNA targets of brain miRNA in ASD and control brain from the STS and PAC combined analysis with previously published and compiled gene lists associated with ASD in (a) studies with transcriptional measurements of ASD brain (Garbett et al. [33] and Voineagu et al. [4]) and (b) sets of genes associated with ASD (Pinto et al. [34] and Abrahams et al. (SFARI) [35]). Numbers indicate KEGG pathways significantly over-represented by mRNA in the respective data set. The 44 KEGG pathways overlapping between this study, SFARI, and Pinto et al. appear in Additional file 1: Table S4

Discussion

This study provides the first evidence of alterations in expression of specific small regulatory non-coding RNA in the temporal cortex of autism spectrum disorder (ASD) brains. The major findings include the following: 1) There are differences in miRNA, snoRNA, and precursor miRNA in both the association cortex of the superior temporal sulcus (STS) and the primary auditory cortex (PAC). While the STS is involved in functions associated with the core features of ASD, the many changes in sncRNAs in primary auditory cortex were not expected as it is not functionally associated with core deficits of ASD. 2) Functional pathways for the gene targets of dysregulated miRNA across both regions were associated with nervous system, cell cycle, cell signaling, and other processes. 3) miRNAs that regulate immune functions were only altered in STS, supporting a role for immune dysfunction in STS associated with the core social deficits of autism. Even a few dysregulated miRNA and other sncRNAs have the potential to affect transcription and/or translation of hundreds of target genes and affect brain signaling, connectivity, and behavior.

Differentially expressed sncRNA in both primary and association cortex in ASD

As hypothesized, based on functional and structural differences known to be present in ASD STS, we found that the STS contained differentially expressed regulatory sncRNA. One mature miRNA, miR-4753-5p, exhibited increased expression levels and another, miR-1, had decreased expression in ASD relative to control tissue. Five immature, stem-loop precursor miRNA showed increased expression, and three had decreased expression in the STS. Given their regulatory role, miRNA can impact hundreds of gene targets and can greatly influence expression levels of mRNA and protein translation [8], potentially making a substantial contribution to differences in overall gene expression levels in ASD.

There were also differences in a number of small nucleolar RNAs (snoRNAs) in STS in ASD. snoRNAs are small RNA molecules that act to guide chemical modifications of other RNAs including ribosomal RNAs, transfer RNAs, and small nuclear RNAs. Three snoRNAs, including ACA39, SNORA11C, SNORA27, exhibited increased expression, and one, U13 (Chr1), had decreased levels. Functionally, sequences within snoRNAs target them to specific mRNA targets, with the C/D box snoRNAs methylating target bases and H/ACA box snoRNAs being associated with pseudouridylation of target bases. These modifications likely affect translation as well as RNA silencing, telomerase maintenance, and alternative splicing [38]. For example, U13 is involved in proper formation of the 3′ end of the 18S ribosomal subunit [39], potentially impacting protein translation and alternative splicing. Aberrant expression of U13 could have widespread cellular effects in the brain. We find that non-coding snoRNAs along with non-coding miRNA are dysregulated in ASD brain, while other studies have also implicated long non-coding RNAs in ASD neuropathology [40, 41].

To our knowledge, there have not been any previous studies targeting molecular abnormalities specifically in primary sensory cortex of ASD brains, making this the first report suggesting regulatory abnormalities in a primary sensory region. In PAC, two mature miRNAs exhibit increased expression, miR-664-3p and miR-4709-3p, and two decreased expression, miR-4742-3p and miR-297. Five snoRNAs showed decreased expression, including SNORA22, SCARNA6, and three paralogs of U13. Four stem-loop precursors were additionally up-regulated, while three had decreased expression.

It is notable that three of the identified precursor miRNA have been previously associated with ASD. miR-4436 has been described as a candidate susceptibility miRNA gene in autism [42]. miR-548 is dysregulated in ASD lymphoblastoid cell lines [43] and targets genes associated with several pathways implicated in ASD including ubiquitin, Wnt, axon guidance, LTP, and natural killer cells [30]. Additionally, it can regulate host antiviral responses via direct targeting of IFN-λ1 [44], which may have relevance for maternal immune activation models of ASD and other neurodevelopmental disorders [45, 46]. miR-146 is reported to have altered expression in ASD cerebellar cortex [47] and in ASD lymphoblastoid cell lines. In the nervous system, miR-146 (a,b) modulates AMPA receptor endocytosis [48] and glial proliferation via Notch signaling [49]. miR-146 also influences the activation state of brain microglia [50]. Microglia are reported to increase in number, show evidence of activation, and show increased gene expression in ASD post-mortem brain tissue [5155]. Moreover, deficiencies in neuron-microglia signaling result in impaired functional brain connectivity and can influence social behavior [56]. Although the functional importance of immature miRNA is not clear, the mature forms of the three stem-loop precursors found in the PAC have previously been associated with ASD.

A number of the predicted gene targets of the miRNA we found different between ASD and typically developing STS have previously been identified in studies looking at gene expression in the temporal cortex [4, 33]. These target genes (INHBA, WDR1, TAGLN2, FNDC3B, ZFP36L1, CNN3, CPLX2, ZNF365, and GLS2) have important roles including those related to “cell communication and motility” and “cell fate and differentiation.” As we observed decreased miR-1 expression and increased miR-4753-5p expression in our study, we would expect to find higher and lower expression of the respective gene transcripts they target in STS. The regulation of these genes and others by altered patterns of miRNA may be at least in part responsible for regional and cellular differences that shape the social and behavioral phenotypes associated with ASD.

Dysregulation of sncRNAs in the PAC of ASD brain is a novel molecular finding we did not predict. Although sensory dysfunction is not considered a core diagnostic feature, neurophysiological studies have reported auditory abnormalities in ASD patients. Abnormalities in cortical auditory-evoked potentials have been reported when sound is presented to the left ear [57]. In contrast, speech stimulation does not produce fMRI abnormalities in PAC of individuals with ASD [23]. Atypical vertical sound localization and sound-onset sensitivity evidenced in ASD cases likely involves primary auditory areas [58]. In addition, a magnetoencephalography study examining differences in auditory cortex function in young ASD and TD children found that there is atypical brain function in the auditory cortex in young children with ASD, regardless of language development [59]. Coupled with neurophysiological findings, the dysregulated sncRNA in PAC in this study provide additional evidence for PAC dysfunction in ASD. Increased attention to primary cortices may be warranted in future investigations of ASD, as such dysfunction could conceivably contribute to cognitive impairments in complex social and linguistic tasks.

miRNA targets incorporated into pathways associated with ASD pathogenesis

We postulated that miRNA play an important role in regulating the transcriptional changes that have been reported in ASD brain tissue. Several genomic and gene expression studies of ASD brain have shown distinct and reproducible differences from typically developing cases that involve neurons/synapses and the immune system. In one study, Voineagu et al. showed consistent differences in transcriptome organization between ASD and typically developing brain by gene co-expression analysis [4]. Moreover, two main modules of co-expressed genes in ASD brain were identified including a neuronal/synapse module that was associated with genes identified by GWAS in ASD and an immune module that was not associated with ASD GWAS genes [4, 60]. These data supported prior studies of neuron, glia, and immune abnormalities in ASD brain [33, 51, 6167]. Another recent whole genome study showed age-related differences of gene expression in ASD brain that also affected neuron and immune pathways including neurogenesis, cell cycle, cell differentiation, cell death, and immune response in young ASD brain [5]. A recent study used weighted gene co-expression network analysis (WGCNA) combined with ASD candidate genes to suggest that cortical projection neurons in layers V and VI of prefrontal/frontal cortex during 10–24 weeks post-conception were key cells in ASD pathogenesis [6]. In a similar WGCNA analysis using different ASD candidate genes, the superficial (pyramidal) glutamatergic cortical projection neurons were implicated in ASD pathogenesis [7].

We found many pathways associated with differentially expressed miRNA in ASD temporal cortex were similarly associated with nervous system and immune function, as well as basic cellular processes previously reported to be dysregulated in ASD. Of the 34 predicted pathways and processes that were represented in the target genes of miRNA in STS and PAC (Fig. 4), a number have been implicated in previous ASD studies including glutamatergic, dopaminergic, and cholinergic synapse pathways, axon guidance, synaptic vesicle cycling, neurotrophins, immune pathways, hedgehog and Wnt signaling, spliceosome, and RNA degradation [24, 27, 60]. Comparing our broader temporal cortex findings with previous analyses, it is notable that 44 of the 75 pathways regulated by differentially expressed miRNA in our study are shared with the 57 pathways reported to be affected by copy number variants in ASD [34] (Fig. 5, Additional file 1: Table S4). Moreover, we found 7 and 15 pathways, respectively, are shared with the 16 and 24 pathways reported to demonstrate alterations in the two prior analyses of mRNA expression ASD temporal cortex [4, 33] (Fig. 5). Common functions associated with pathways across all studies of ASD temporal cortex include cell cycle, metabolic, and neuron-related pathways and processes. Given the substantial number of convergent pathways between studies of temporal cortex mRNA and miRNA expression, it is highly plausible that miRNA are exerting a regulatory influence over mRNA associated with atypical cortical function in ASD.

Within the temporal lobe, two pathways were shared between PAC and STS, the “PI3K-Akt pathway” and “pathways in cancer.” The PI3K-Akt pathway is of interest as it has been implicated in a number of previous ASD studies. The PI3K-Akt/mTOR pathway is dysregulated in monogenic disorders associated with ASD including Fragile X, tuberous sclerosis, neurofibromatosis, and PTEN [68]. Rare and even some common genetic variants associated with ASD affect PI3K-Akt signaling [2, 69]. Studies of gene expression in blood of ASD subjects have found abnormal PI3K-Akt/mTOR signaling [70]. Our studies of alternative splicing of mRNA expressed in blood leukocytes detected abnormalities of PI3K-Akt signaling in about half the 2–4 year old ASD children with idiopathic ASD [27]. Thus, the current study suggests abnormal PI3K-Akt signaling in temporal cortex of ASD subjects in both primary auditory and STS association cortex, and previous studies suggest PI3K-Akt abnormalities in blood and differential alternative splicing of PI3K-Akt pathways in blood of many ASD subjects. These results would support this as a potential therapeutic target in at least a subgroup of ASD subjects as has been suggested by others [71].

One of the major findings of this study is that few predicted miRNA-regulated pathways were shared between the PAC and STS, despite the fact that they are adjacent cortical structures. Because many pathways are specific for PAC or STS, there is evidence that differential regulation of miRNA, their mRNA targets, and associated pathways can be region specific in ASD.

Differentially expressed miRNA having a significant number of gene targets in immune and most of the cell cycle KEGG pathways (except “pathways in cancer”) were restricted only to STS. Cell cycle genes might be related to neuronal cell death and/or to glial proliferation and perhaps linked to increased microglial and neuron numbers in ASD brain [5153, 63, 72]. Observed immune pathways in this study, related to both viruses and bacteria, are consistent with many reports of immune pathway activation in ASD brain and blood [4, 5, 33, 60, 7375]. Since the immune regulated genes in ASD brain were not among the GWAS-associated ASD genes [4], these immune genes might reflect the environmental factor in ASD pathogenesis. ASD is considered to be a complex genetic disorder in which several genes interact with an environmental precipitant. Activation of immune and cell cycle pathways in STS might suggest that this is a region of primary pathology, and lack of immune and cell cycle activation in PAC might suggest gene changes in PAC are secondary to those in STS association cortex because of interconnections between the two regions. At the very least, these findings suggest molecular and immune disruption pathologies present in the temporal association cortex that are not shared with primary sensory cortex. This further emphasizes the importance of assessing molecular abnormalities in distinct cortical regions.

Limitations of the study

Future studies will need to confirm the miRNA changes discovered here with microarrays using qRT-PCR, though in other studies we have shown most miRNA regulated on arrays can be confirmed using qRT-PCR [28]. In addition, future studies will need to confirm the changes in a separate cohort, to examine mRNA targets of the miRNA and their associated pathways, and to examine other cortical regions which may also be associated with core ASD symptoms. Future studies of primarily young brains are needed to assess changes around the time of diagnosis, though postmortem studies like this by their nature tend to have older subjects. Though age, gender, and other factors were included in the statistical model, some covariates like medications, cause of death, and others could not be included because of the limited number of samples.

Conclusions

This study provides the first demonstration of dysregulation of sncRNA in superior temporal cortex of ASD brains that could regulate hundreds of target transcripts and many predicted pathways and processes that have been implicated in previous ASD studies. In addition to changes in sncRNA expression in the STS, transcriptional dysregulation in the PAC could play a role in social and language dysfunction in ASD. This has implications for understanding the pathogenesis of ASD and ascertaining whether the “primary ASD pathology” is only in association cortex and whether abnormalities of both primary auditory cortex and STS temporal association cortex contribute to social and language behavioral dysfunction observed in ASD.

Abbreviations

ASD:

autism spectrum disorders

CNV:

copy number variant

GWAS:

genome-wide association study

KEGG:

Kyoto Encyclopedia of Genes and Genomes

miRNA:

microRNA

PAC:

primary auditory cortex of temporal lobe

PMI:

postmortem interval

sncRNA:

small non-coding RNA

snoRNA:

small nucleolar RNA

SNP:

single nucleotide polymorphism

STS:

superior temporal sulcus of temporal lobe

WGCNA:

weighted gene co-expression network analysis

References

  1. Holt R, Monaco AP. Links between genetics and pathophysiology in the autism spectrum disorders. EMBO Mol Med. 2011;3:438–50.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  2. Levitt P, Campbell DB. The genetic and neurobiologic compass points toward common signaling dysfunctions in autism spectrum disorders. J Clin Invest. 2009;119:747–54.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  3. Geschwind DH. Genetics of autism spectrum disorders. Trends Cogn Sci. 2011;15:409–16.

    Article  PubMed Central  PubMed  Google Scholar 

  4. Voineagu I, Wang X, Johnston P, Lowe JK, Tian Y, Horvath S, et al. Transcriptomic analysis of autistic brain reveals convergent molecular pathology. Nature. 2011;474:380–4.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  5. Chow ML, Pramparo T, Winn ME, Barnes CC, Li HR, Weiss L, 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:e1002592.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  6. Willsey AJ, Sanders SJ, Li M, Dong S, Tebbenkamp AT, Muhle RA, et al. Coexpression networks implicate human midfetal deep cortical projection neurons in the pathogenesis of autism. Cell. 2013;155:997–1007.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  7. Parikshak NN, Luo R, Zhang A, Won H, Lowe JK, Chandran V, et al. Integrative functional genomic analyses implicate specific molecular pathways and circuits in autism. Cell. 2013;155:1008–21.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  8. Follert P, Cremer H, Beclin C. MicroRNAs in brain development and function: a matter of flexibility and stability. Front Mol Neurosci. 2014;7:5.

    Article  PubMed Central  PubMed  Google Scholar 

  9. Vaishnavi V, Manikandan M, Munirajan AK. Mining the 3′UTR of autism-implicated genes for SNPs perturbing microRNA regulation. Genomics Proteomics Bioinformatics. 2014;12:92–104.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  10. Vaishnavi V, Manikandan M, Tiwary BK, Munirajan AK. Insights on the functional impact of microRNAs present in autism-associated copy number variants. PLoS One. 2013;8:e56781.

    Article  PubMed Central  PubMed  Google Scholar 

  11. Salem AM, Ismail S, Zarouk WA, Abdul Baky O, Sayed AA, Abd El-Hamid S, et al. Genetic variants of neurotransmitter-related genes and miRNAs in Egyptian autistic patients. The Scientific World JOURNAL. 2013;2013:670621.

    Article  PubMed Central  Google Scholar 

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

    Article  PubMed Central  PubMed  Google Scholar 

  13. Mellios N, Sur M. The emerging role of microRNAs in schizophrenia and autism spectrum disorders. Frontiers Psychiatry. 2012;3:39.

    Article  Google Scholar 

  14. Hawrylycz MJ, Lein ES, Guillozet-Bongaarts AL, Shen EH, Ng L, Miller JA, et al. An anatomically comprehensive atlas of the adult human brain transcriptome. Nature. 2012;489:391–9.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  15. 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:123–42.

    Article  PubMed  Google Scholar 

  16. Nummenmaa L, Calder AJ. Neural mechanisms of social attention. Trends Cogn Sci. 2009;13:135–43.

    Article  PubMed  Google Scholar 

  17. Rushworth MF, Mars RB, Sallet J. Are there specialized circuits for social cognition and are they unique to humans? Curr Opin Neurobiol. 2013;23:436–42.

    Article  CAS  PubMed  Google Scholar 

  18. Hein G, Knight RT. Superior temporal sulcus—It’s my area: or is it? J Cogn Neurosci. 2008;20:2125–36.

    Article  PubMed  Google Scholar 

  19. Emmorey K, McCullough S. The bimodal bilingual brain: effects of sign language experience. Brain Lang. 2009;109:124–32.

    Article  PubMed Central  PubMed  Google Scholar 

  20. Stefanatos GA, Baron IS. The ontogenesis of language impairment in autism: a neuropsychological perspective. Neuropsychol Rev. 2011;21:252–70.

    Article  PubMed  Google Scholar 

  21. Sugranyes G, Kyriakopoulos M, Corrigall R, Taylor E, Frangou S. Autism spectrum disorders and schizophrenia: meta-analysis of the neural correlates of social cognition. PLoS One. 2011;6:e25322.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  22. Saitovitch A, Bargiacchi A, Chabane N, Brunelle F, Samson Y, Boddaert N, et al. Social cognition and the superior temporal sulcus: implications in autism. Rev Neurol. 2012;168:762–70.

    Article  CAS  PubMed  Google Scholar 

  23. Lai G, Schneider HD, Schwarzenberger JC, Hirsch J. Speech stimulation during functional MR imaging as a potential indicator of autism. Radiology. 2011;260:521–30.

    Article  PubMed  Google Scholar 

  24. Da Costa S, van der Zwaag W, Marques JP, Frackowiak RS, Clarke S, Saenz M. Human primary auditory cortex follows the shape of Heschl’s gyrus. J Neurosci. 2011;31:14067–75.

    Article  PubMed  Google Scholar 

  25. Fullerton BC, Pandya DN. Architectonic analysis of the auditory-related areas of the superior temporal region in human brain. J Comp Neurol. 2007;504:470–98.

    Article  PubMed  Google Scholar 

  26. Barger N, Sheley MF, Schumann CM. Stereological study of pyramidal neurons in the human superior temporal gyrus from childhood to adulthood. J Comp Neurol. 2014:523(7). doi:10.1002/cne.23707

  27. Stamova BS, Tian Y, Nordahl CW, Shen MD, Rogers S, Amaral DG, et al. Evidence for differential alternative splicing in blood of young boys with autism spectrum disorders. Mol Autism. 2013;4:30.

    Article  PubMed Central  PubMed  Google Scholar 

  28. Jickling GC, Ander BP, Zhan X, Noblett D, Stamova B, Liu D. microRNA expression in peripheral blood cells following acute ischemic stroke and their predicted gene targets. PLoS One. 2014;9:e99283.

    Article  PubMed Central  PubMed  Google Scholar 

  29. Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, et al. Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003;4:249–64.

    Article  PubMed  Google Scholar 

  30. Enstrom AM, Lit L, Onore CE, Gregg JP, Hansen RL, Pessah IN, et al. Altered gene expression and function of peripheral blood natural killer cells in children with autism. Brain Behav Immun. 2009;23:124–33.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  31. Stamova B, Green PG, Tian Y, Hertz-Picciotto I, Pessah IN, Hansen R, et al. Correlations between gene expression and mercury levels in blood of boys with and without autism. Neurotox Res. 2011;19:31–48.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  32. Vlachos IS, Kostoulas N, Vergoulis T, Georgakilas G, Reczko M, Maragkakis M, et al. DIANA miRPath v. 2.0: investigating the combinatorial effect of microRNAs in pathways. Nucleic Acids Res. 2012;40:W498–504.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  33. Garbett K, Ebert PJ, Mitchell A, Lintas C, Manzi B, Mirnics K, et al. Immune transcriptome alterations in the temporal cortex of subjects with autism. Neurobiol Dis. 2008;30:303–11.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  34. Pinto D, Pagnamenta AT, Klei L, Anney R, Merico D, Regan R, et al. Functional impact of global rare copy number variation in autism spectrum disorders. Nature. 2010;466:368–72.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  35. Abrahams BS, Arking DE, Campbell DB, Mefford HC, Morrow EM, Weiss LA, et al. SFARI Gene 2.0: a community-driven knowledgebase for the autism spectrum disorders (ASDs). Mol. Autism. 2013;4:36.

    Article  Google Scholar 

  36. Wang J, Duncan D, Shi Z, Zhang B. WEB-based GEne SeT AnaLysis Toolkit (WebGestalt): update 2013. Nucleic Acids Res. 2013;41:W77–83.

    Article  PubMed Central  PubMed  Google Scholar 

  37. Micallef L, Rodgers P. eulerAPE: drawing area-proportional 3-Venn diagrams using ellipses. PLoS One. 2014;9:e101717.

    Article  PubMed Central  PubMed  Google Scholar 

  38. Lui L, Lowe T. Small nucleolar RNAs and RNA-guided post-transcriptional modification. Essays Biochem. 2013;54:53–77.

    Article  CAS  PubMed  Google Scholar 

  39. Cavaille J, Hadjiolov AA, Bachellerie JP. Processing of mammalian rRNA precursors at the 3′ end of 18S rRNA. Identification of cis-acting signals suggests the involvement of U13 small nucleolar RNA. Eur J Biochem. 1996;242:206–13.

    Article  CAS  PubMed  Google Scholar 

  40. Wilkinson B, Campbell DB. Contribution of long noncoding RNAs to autism spectrum disorder risk. Int Rev Neurobiol. 2013;113:35–59.

    Article  CAS  PubMed  Google Scholar 

  41. Ziats MN, Rennert OM. Aberrant expression of long noncoding RNAs in autistic brain. J Mol Neurosci. 2013;49:589–93.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  42. Marrale M, Albanese NN, Cali F, Romano V. Assessing the impact of copy number variants on miRNA genes in autism by Monte Carlo simulation. PLoS One. 2014;9:e90947.

    Article  PubMed Central  PubMed  Google Scholar 

  43. Ghahramani Seno MM, Hu P, Gwadry FG, Pinto D, Marshall CR, Casallo G, et al. Gene and miRNA expression profiles in autism spectrum disorders. Brain Res. 2011;1380:85–97.

    Article  CAS  PubMed  Google Scholar 

  44. Li Y, Xie J, Xu X, Wang J, Ao F, Wan Y, et al. MicroRNA-548 down-regulates host antiviral response via direct targeting of IFN-lambda1. Protein Cell. 2013;4:130–41.

    Article  CAS  PubMed  Google Scholar 

  45. Patterson PH. Maternal infection and immune involvement in autism. Trends Mol Med. 2011;17:389–94.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  46. Bauman MD, Iosif AM, Smith SE, Bregere C, Amaral DG, Patterson PH. Activation of the maternal immune system during pregnancy alters behavioral development of rhesus monkey offspring. Biol Psychiatry. 2014;75:332–41.

    Article  CAS  PubMed  Google Scholar 

  47. Abu-Elneel K, Liu T, Gazzaniga FS, Nishimura Y, Wall DP, Geschwind DH, et al. Heterogeneous dysregulation of microRNAs across the autism spectrum. Neurogenetics. 2008;9:153–61.

    Article  CAS  PubMed  Google Scholar 

  48. Chen YL, Shen CK. Modulation of mGluR-dependent MAP1B translation and AMPA receptor endocytosis by microRNA miR-146a-5p. J Neurosci. 2013;33:9013–20.

    Article  CAS  PubMed  Google Scholar 

  49. Mei J, Bachoo R, Zhang CL. MicroRNA-146a inhibits glioma development by targeting Notch1. Mol Cell Biol. 2011;31:3584–92.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  50. Saba R, Gushue S, Huzarewich RL, Manguiat K, Medina S, Robertson C, et al. MicroRNA 146a (miR-146a) is over-expressed during prion disease and modulates the innate immune response and the microglial activation state. PLoS One. 2012;7:e30832.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  51. Vargas DL, Nascimbene C, Krishnan C, Zimmerman AW, Pardo CA. Neuroglial activation and neuroinflammation in the brain of patients with autism. Ann Neurol. 2005;57:67–81.

    Article  CAS  PubMed  Google Scholar 

  52. Morgan JT, Chana G, Pardo CA, Achim C, Semendeferi K, Buckwalter J, et al. Microglial activation and increased microglial density observed in the dorsolateral prefrontal cortex in autism. Biol Psychiatry. 2010;68:368–76.

    Article  PubMed  Google Scholar 

  53. Kern JK, Geier DA, Sykes LK, Geier MR. Evidence of neurodegeneration in autism spectrum disorder. Translational Neurodegeneration. 2013;2:17.

    Article  PubMed Central  PubMed  Google Scholar 

  54. Suzuki K, Sugihara G, Ouchi Y, Nakamura K, Futatsubashi M, Takebayashi K, et al. Microglial activation in young adults with autism spectrum disorder. JAMA Psychiatry. 2013;70:49–58.

    Article  PubMed  Google Scholar 

  55. Edmonson C, Ziats MN, Rennert OM. Altered glial marker expression in autistic post-mortem prefrontal cortex and cerebellum. Mol Autism. 2014;5:3.

    Article  PubMed Central  PubMed  Google Scholar 

  56. Zhan Y, Paolicelli RC, Sforazzini F, Weinhard L, Bolasco G, Pagani F, et al. Deficient neuron-microglia signaling results in impaired functional brain connectivity and social behavior. Nat Neurosci. 2014;17:400–6.

    Article  CAS  PubMed  Google Scholar 

  57. Stroganova TA, Kozunov VV, Posikera IN, Galuta IA, Gratchev VV, Orekhova EV. Abnormal pre-attentive arousal in young children with autism spectrum disorder contributes to their atypical auditory behavior: an ERP study. PLoS One. 2013;8:e69100.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  58. Visser E, Zwiers MP, Kan CC, Hoekstra L, van Opstal AJ, Buitelaar JK. Atypical vertical sound localization and sound-onset sensitivity in people with autism spectrum disorders. J Psychiatry Neurosci. 2013;38:398–406.

    Article  PubMed Central  PubMed  Google Scholar 

  59. Yoshimura Y, Kikuchi M, Shitamichi K, Ueno S, Munesue T, Ono Y, et al. Atypical brain lateralisation in the auditory cortex and language performance in 3- to 7-year-old children with high-functioning autism spectrum disorder: a child-customised magnetoencephalography (MEG) study. Mol Autism. 2013;4:38.

    Article  PubMed Central  PubMed  Google Scholar 

  60. Voineagu I, Eapen V. Converging pathways in autism spectrum disorders: interplay between synaptic dysfunction and immune responses. Front Hum Neurosci. 2013;7:738.

    Article  PubMed Central  PubMed  Google Scholar 

  61. Purcell AE, Jeon OH, Zimmerman AW, Blue ME, Pevsner J. Postmortem brain abnormalities of the glutamate neurotransmitter system in autism. Neurology. 2001;57:1618–28.

    Article  CAS  PubMed  Google Scholar 

  62. Schumann CM, Amaral DG. Stereological analysis of amygdala neuron number in autism. J Neurosci. 2006;26:7674–9.

    Article  CAS  PubMed  Google Scholar 

  63. 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:e110356.

    Article  PubMed Central  PubMed  Google Scholar 

  64. van Kooten IA, Palmen SJ, von Cappeln P, Steinbusch HW, Korr H, Heinsen H, et al. Neurons in the fusiform gyrus are fewer and smaller in autism. Brain. 2008;131:987–99.

    Article  PubMed  Google Scholar 

  65. Uppal N, Wicinski B, Buxbaum JD, Heinsen H, Schmitz C, Hof PR. Neuropathology of the anterior midcingulate cortex in young children with autism. J Neuropathol Exp Neurol. 2014;73:891–902.

    Article  PubMed  Google Scholar 

  66. Wegiel J, Flory M, Kuchna I, Nowicki K, Ma S, Imaki H, et al. Stereological study of the neuronal number and volume of 38 brain subdivisions of subjects diagnosed with autism reveals significant alterations restricted to the striatum, amygdala and cerebellum. Acta Neuropathologica Communications. 2014;2:141.

    Article  PubMed Central  PubMed  Google Scholar 

  67. Wegiel J, Flory M, Kuchna I, Nowicki K, Ma SY, Imaki H, et al. Brain-region-specific alterations of the trajectories of neuronal volume growth throughout the lifespan in autism. Acta Neuropathologica Communications. 2014;2:28.

    Article  PubMed Central  PubMed  Google Scholar 

  68. Chen J, Alberts I, Li X. Dysregulation of the IGF-I/PI3K/AKT/mTOR signaling pathway in autism spectrum disorders. Int J Dev Neurosci. 2014;35:35–41.

    Article  PubMed  Google Scholar 

  69. Lee TL, Raygada MJ, Rennert OM. Integrative gene network analysis provides novel regulatory relationships, genetic contributions and susceptible targets in autism spectrum disorders. Gene. 2012;496:88–96.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  70. Kong SW, Collins CD, Shimizu-Motohashi Y, Holm IA, Campbell MG, Lee IH, et al. Characteristics and predictive value of blood transcriptome signature in males with autism spectrum disorders. PLoS One. 2012;7:e49475.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  71. Kleijer KT, Schmeisser MJ, Krueger DD, Boeckers TM, Scheiffele P, Bourgeron T, et al. Neurobiology of autism gene products: towards pathogenesis and drug targets. Psychopharmacology. 2014;231:1037–62.

    Article  CAS  PubMed  Google Scholar 

  72. Schumann CM, Amaral DG. Stereological estimation of the number of neurons in the human amygdaloid complex. J Comp Neurol. 2005;491:320–9.

    Article  PubMed Central  PubMed  Google Scholar 

  73. Ziats MN, Rennert OM. Expression profiling of autism candidate genes during human brain development implicates central immune signaling pathways. PLoS One. 2011;6:e24691.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  74. Ashwood P, Corbett BA, Kantor A, Schulman H, Van de Water J, Amaral DG. In search of cellular immunophenotypes in the blood of children with autism. PLoS One. 2011;6:e19299.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  75. Ginsberg MR, Rubin RA, Falcone T, Ting AH, Natowicz MR. Brain transcriptional and epigenetic associations with autism. PLoS One. 2012;7:e44736.

    Article  PubMed Central  PubMed  Google Scholar 

Download references

Acknowledgements

The authors appreciate the technical assistance of Alicja Omanska, Kenneth Lam, and Carolyn Ho of the UC Davis MIND Institute, as well Ryan Davis and Stephenie Liu of the Genomics Shared Resource at the UC Davis Comprehensive Cancer Center. Financial support came from the MIND Institute (CMS, FRS), IDDRC U54 HD079125, UC Davis Children’s Miracle Network (FRS, BS), and R01 MH 097236 (CMS). The authors would also like to express their gratitude to the brain donors and their families who made this research possible. Biological specimens and clinical data used in this research was obtained from Autism BrainNet, previously Autism Tissue Program, sponsored by the Simons Foundation and Autism Speaks. Postmortem Confirmation of Consent was obtained from next-of-kin for use of donor brain tissue and publication of individual details presented in this manuscript. Consent forms are held by the Harvard Brain Tissue Resource Center and are available for review by the Editor-in-Chief.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Cynthia M. Schumann.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

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

FRS and CMS are Co-senior author.

Additional file

Additional file 1: Table S1-S4.

Table 1: P values and fold changes for the 14 sncRNA different between ASD and CTRL STS. Table S2: P values and fold changes for the 16 sncRNA different between ASD and CTRL PAC. Table S3: The 445 predicted gene targets of STS and 1214 predicted gene targets of miRNA in PAC based on the microT-CDS algorithm. Table S4: Pathways common to our study, SFARI ASD implicated genes and Pinto et al. 2010 (the 44 common KEGG pathways of Fig. 5B).

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.

The Creative Commons Public Domain Dedication waiver (https://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ander, B.P., Barger, N., Stamova, B. et al. Atypical miRNA expression in temporal cortex associated with dysregulation of immune, cell cycle, and other pathways in autism spectrum disorders. Molecular Autism 6, 37 (2015). https://doi.org/10.1186/s13229-015-0029-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13229-015-0029-9

Keywords