Skip to main content

Advertisement

Springer Nature is making SARS-CoV-2 and COVID-19 research free. View research | View latest news | Sign up for updates

Integrated genome-wide Alu methylation and transcriptome profiling analyses reveal novel epigenetic regulatory networks associated with autism spectrum disorder

Molecular Autism volume 9, Article number: 27 (2018) Cite this article

Abstract

Background

Alu elements are a group of repetitive elements that can influence gene expression through CpG residues and transcription factor binding. Altered gene expression and methylation profiles have been reported in various tissues and cell lines from individuals with autism spectrum disorder (ASD). However, the role of Alu elements in ASD remains unclear. We thus investigated whether Alu elements are associated with altered gene expression profiles in ASD.

Methods

We obtained five blood-based gene expression profiles from the Gene Expression Omnibus database and human Alu-inserted gene lists from the TranspoGene database. Differentially expressed genes (DEGs) in ASD were identified from each study and overlapped with the human Alu-inserted genes. The biological functions and networks of Alu-inserted DEGs were then predicted by Ingenuity Pathway Analysis (IPA). A combined bisulfite restriction analysis of lymphoblastoid cell lines (LCLs) derived from 36 ASD and 20 sex- and age-matched unaffected individuals was performed to assess the global DNA methylation levels within Alu elements, and the Alu expression levels were determined by quantitative RT-PCR.

Results

In ASD blood or blood-derived cells, 320 Alu-inserted genes were reproducibly differentially expressed. Biological function and pathway analysis showed that these genes were significantly associated with neurodevelopmental disorders and neurological functions involved in ASD etiology. Interestingly, estrogen receptor and androgen signaling pathways implicated in the sex bias of ASD, as well as IL-6 signaling and neuroinflammation signaling pathways, were also highlighted. Alu methylation was not significantly different between the ASD and sex- and age-matched control groups. However, significantly altered Alu methylation patterns were observed in ASD cases sub-grouped based on Autism Diagnostic Interview-Revised scores compared with matched controls. Quantitative RT-PCR analysis of Alu expression also showed significant differences between ASD subgroups. Interestingly, Alu expression was correlated with methylation status in one phenotypic ASD subgroup.

Conclusion

Alu methylation and expression were altered in LCLs from ASD subgroups. Our findings highlight the association of Alu elements with gene dysregulation in ASD blood samples and warrant further investigation. Moreover, the classification of ASD individuals into subgroups based on phenotypes may be beneficial and could provide insights into the still unknown etiology and the underlying mechanisms of ASD.

Background

Autism spectrum disorder (ASD) refers to a group of complex neurodevelopmental disorders that are characterized according to the Diagnostic and Statistical Manual of Mental Disorders, Fifth Edition (DSM-5) criteria by two domains: (i) behavioral impairment, including significant deficits in social interactions and communication, and (ii) restricted interests and repetitive behaviors [1]. Recent data released from the Autism and Developmental Disabilities Monitoring Network have shown a 78% increase in ASD prevalence over the past decade, and approximately 1 in 68 children in the USA have ASD [2]. This increase in ASD prevalence leads to a large economic burden, including costs for healthcare, ASD-related therapy, family-coordinated services, and special education systems [3].

A number of studies have supported the hypothesis that genetic factors are strongly associated with the etiology and susceptibility of ASD. However, abnormalities in genomic DNA are found in only 10–20% of ASD cases accumulatively, partly due to the etiological heterogeneity of ASD, which has a wide variety of different risk factors in addition to genetic factors [4]. A broad variability in clinical phenotypes of ASD individuals are thought to result from complicated interactions between genetic and environmental factors that increase ASD risk [5, 6]. Although the concordance rate among monozygotic (MZ) twins was found to range from 60% to as high as 90%, notable discordance in the ASD diagnosis within monozygotic twin pairs and significant differences in the ASD severity within ASD-concordant monozygotic twin pairs have also been observed [7, 8]. This evidence strongly suggests that environmental factors may play an important role in the etiology and/or the susceptibility of ASD.

DNA methylation is a major epigenetic regulator of gene expression and associated phenotypes. Methylation patterns in genomic DNA are generated during embryogenesis and early fetal development and are altered throughout life in response to endogenous or exogenous environmental signals. Recently, differentially methylated variants (DMVs) at specific CpG sites or differentially methylated regions (DMRs) have been investigated in ASD individuals using various types of tissues, including lymphoblastoid cell lines (LCLs) [9], whole blood [10], and brain tissue samples [11,12,13]. Nguyen and colleagues performed a large-scale methylation profiling analysis of LCLs derived from discordantly diagnosed (i.e., ASD and non-ASD) monozygotic twins and sibling pairs [9]. The results revealed differentially methylated genes associated with several biological functions, including gene transcription, nervous system development, and other biological mechanisms implicated in ASD. Moreover, the mRNA expression level of one ASD candidate gene (e.g., retinoic acid-related orphan receptor-alpha gene, RORA), which exhibits a differential methylation pattern in ASD, was also decreased in LCLs. Reduced RORA protein levels were also observed in the brain tissues of ASD individuals. This finding suggests that molecular changes in ASD peripheral cells may reflect at least some pathobiological conditions in the brain.

Several studies involving the gene expression profiling of blood or blood-derived cell lines from ASD and non-ASD subjects reported different (but somewhat overlapping) sets of differentially expressed genes (DEGs), most likely due to the heterogeneity within the ASD population. To reduce phenotypic variability and increase statistical power, Hu and Steinberg defined phenotypes within ASD based on cluster analyses of Autism Diagnostic Interview-Revised (ADI-R) scores [14]. Interestingly, gene expression profiling among the LCLs of ASD subgroups also showed differential expression [15]. This result suggested that the sub-classification of ASD patients could help identify subphenotype-specific risk factors in heterogeneous ASD populations. However, most global methylation and gene expression profiling studies have focused on protein-coding regions rather than noncoding regions that include repetitive sequences.

Alu elements are a group of repetitive sequences or mobile genetic elements with copy numbers in excess of one million in the human genome, thus contributing to almost 11% of the human genome [16]. Alu elements belong to a class of retrotransposons termed SINEs (short interspersed elements). Several reports have demonstrated that Alu elements can influence gene expression via insertion into the gene structure and can attract transcription factor binding to regulate gene expression [17]. Alu elements have many CpG residues in their sequences, and these CpG residues are common methylation sites. Methylated CpGs represent approximately 23% of all methylated residues in the human genome [18]. DNA methylation can be altered by exposure to environmental factors that can reduce Alu element methylation in human tissues [19], and there is evidence that Alu methylation plays an important role in cell proliferation and resistance to DNA damage [20]. However, the role of Alu elements and their methylation in ASD remains unclear.

In this study, we therefore aimed to investigate the association between Alu elements and altered gene expression profiles in ASD. First, publicly available gene expression data from previously published ASD transcriptome profiling studies were downloaded from the GEO DataSets database. DEGs from each ASD transcriptome study were identified, and the association between the DEGs from each study and human Alu-inserted genes were determined by Fisher’s exact test. Moreover, the DEGs with Alu insertion were identified and subsequently subjected to Ingenuity Pathway Analysis (IPA) to predict the biological functions, canonical pathways, and gene regulatory networks associated with ASD. A combined bisulfite restriction analysis (COBRA) of Alu was then performed to assess DNA methylation within Alu elements and associated CpGs, which might regulate the expression of Alu elements and Alu-inserted genes. The schematic diagram showing the experimental workflow of this study is illustrated in Fig. 1.

Fig. 1
figure1

Schematic diagram of experimental workflow. Our workflow initiated with the acquisition of blood-based gene expression profiles from GEO DataSets and human Alu-inserted gene lists. Fisher’s exact test was then used to identify differentially expressed genes (DEGs) with Alu insertions. A total of 320 overlapping genes among the selected study results were used to predict biological functions, diseases, and gene regulatory networks. Fifty-six LCLs were used as a model to investigate the association between the Alu methylation status and Alu expression profiles in LCLs

Full size image

Methods

Data collection

Gene expression profiles of ASD and control individuals were obtained from the NCBI Gene Expression Omnibus database (GEO DataSets: http://www.ncbi.nlm.nih.gov/gds) using the following criteria: the included studies must be ASD studies; the studies must include blood-based gene expression profiles from microarray experiments; and the sample sizes must be greater than or equal to 40 samples. All supplementary information, including series matrix files and related platforms, was freely available in GEO DataSets [21, 22]. The details of each study are provided in Table 1.

Table 1 Details of the gene expression profiles obtained from GEO DataSets
Full size table

Alu subfamily-inserted gene lists of Human Genome 18 (UCSC hg18, NCBI build 36.1) were downloaded from the TranspoGene database (http://transpogene.tau.ac.il), which is a publicly available database of transposed elements (TEs) located within protein-coding genes [23]. For the human Alu-inserted genes, a single gene can be inserted by multiple Alu elements with different insertion types. In this study, we obtained all genes with at least one instance of one of the four types of Alu insertions, namely, the exonic, exonized, intronic, and promoter inserts (Table 2). All human Alu subfamily consensus sequences were obtained from Repbase, which is a database of repetitive elements in eukaryotic genomes [24].

Table 2 Total number of the Alu-inserted gene lists from TranspoGene database
Full size table

Identification of DEGs and association with Alu-inserted genes

The transcriptome profile from each study was analyzed separately by Multiple Experiment Viewer (MEV) [25]. All transcriptome data were filtered using a 70% cutoff, which removes transcripts for which intensity values are missing in > 30% of the samples. The available transcripts were used for identifying DEGs in ASD with two-tailed t tests with adjusted Bonferroni correction. All studies obtained from GEO DataSets employed ASD samples vs. sex- and/or age-matched controls in the analyses. The DEGs and non-DEGs were intersected with the human Alu-inserted gene lists, and the number of intersected genes was classified based on a crosstab 2 × 2 table into four categories, namely, DEGs with Alu insertion, DEGs without Alu insertion, non-DEGs with Alu insertion, and non-DEGs without Alu insertion. Fisher’s exact test was then used to determine whether the DEG distributions were dependent on the human Alu-inserted gene lists. Moreover, the DEGs were classified as downregulated or upregulated and used for comparison with the human Alu-inserted gene lists. These processes were repeated with each type of human Alu insertion, including intronic, exonized, exonic, and promoter type. A Fisher’s exact test P value with Benjamini-Hochberg correction (FDR = 0.05) of less than 0.05 was considered significant.

To identify the reproducibility of the human Alu-inserted genes that were differentially expressed in peripheral blood and cell lines derived from ASD individuals, the significant gene lists from the individual studies were used to create Venn diagrams (http://bioinformatics.psb.ugent.be/webtools/Venn). The DEGs found in at least two studies were selected to identify biological functions and gene regulatory networks that were enriched by Alu elements in ASD.

Identification of biological functions, canonical pathways, and gene regulatory networks associated with Alu elements in ASD

IPA (QIAGEN Inc., https://www.qiagenbioinformatics.com/products/ingenuity-pathway-analysis/) is a powerful bioinformatics tool that is helpful for understanding complex “-omics” data, such as data from microarray experiments. In this study, IPA was used to identify the biological functions, diseases, canonical pathways, and gene regulatory networks of genes that were identified in at least two studies using Fisher’s exact tests with Benjamini-Hochberg correction (FDR = 0.05, P value < 0.05).

Experimental models and cell culture

LCLs derived from the peripheral lymphocytes of male individuals were obtained from the Autism Genetic Resource Exchange Repository (AGRE, Los Angeles, CA, USA). Our subjects included LCLs from 36 ASD individuals and 20 sex- and age-matched unaffected controls. These individuals were previously used in large-scale clustering analysis for identification of phenotypic subgroups within ASD based on the Autism Diagnostic Interview-Revised (ADI-R) scores, as previously described in detail [14]. As a result, a total of 1954 individual ASD probands were subdivided into four phenotypic groups based on the scores of the ADI-R questionnaire. Then, we selected representative samples from the groups after excluding ASD individuals with cognitive impairment (Raven’s scores < 70), known genetic or chromosomal abnormalities (e.g., Fragile X, Rett syndrome, tuberous sclerosis, chromosome 15q11–q13 duplication), or diagnosed comorbid psychiatric disorders (e.g., bipolar disorder, obsessive compulsive disorder, severe anxiety). Impairment in spoken language was also confirmed based on low standard scores (< 80) on the Peabody Picture Vocabulary Test. Individuals born prematurely (< 35 weeks gestation) were also excluded from this study. These exclusion criteria are expected to reduce the heterogeneity of subjects to study idiopathic autism. Our LCLs represented individuals from three phenotypic groups: severe language impairment (subgroup L), milder symptoms (subgroup M), and savant skills (subgroup S). The demographic information of the LCLs used in this study is shown in Additional file 1. The LCLs were cultured according to the protocol of the Rutgers University Cell and DNA Repository, which maintains the AGRE collection of biological materials from autistic individuals and their relatives. Briefly, cells were cultured in RPMI 1640 medium supplemented with 15% fetal bovine serum and 1% penicillin/streptomycin. The cultures were split 1:2 every 3–4 days and harvested for DNA and RNA isolation 3 days after splitting, when the cultures were in the logarithmic growth phase.

Quantitative reverse transcription-PCR analysis

Total RNA was isolated from LCLs of ASD individuals and sex- and age-matched unaffected controls using the GENEzol Reagent (Geneaid, Taiwan) according to the manufacturer’s recommended protocol. The RNA concentration was determined using a NanoDrop 1000 spectrophotometer (Thermo Scientific, USA). Quantitative reverse-transcription-PCR (RT-PCR) analysis was used to determine AluS subfamily expression in the LCLs of ASD individuals and controls. First, 5 μg of extracted RNA was treated with DNase enzyme in a 10-μl reaction (RQ1 RNase-Free DNase, Promega), and then, 2 μl of the DNase-treated RNA was reverse-transcribed to complementary DNA (cDNA) using AccuPower® RT PreMix (Bioneer, Korea) and oligo dT18 primer in a volume of 20 μl, according to the manufacturer’s instructions. The quantitative PCR assay was performed in triplicate using 1 μl of the cDNA in master mix reactions according to the manufacturer’s instructions (AccuPower® 2X GreenStar™ qPCR MasterMix, Bioneer, Korea). The amplification cycles consisted of an initial denaturing cycle at 95 °C for 15 min followed by 40 cycles of 45 s at 95 °C for denaturing and 45 s at 60 °C for annealing/extension. Product formation was confirmed by melting curve analysis (55 to 94 °C). The AluS transcript-specific primers were as follows: forward 5′-GTGGCTCACGCCTGTAATC-3′ and reverse 5′-GTAGAGACGGGGTTTCACCA-3′. The number of AluS transcripts was normalized to the housekeeping gene GAPDH whose expression was measured using the following primer sequences (forward 5′-ATGTTCGTCATGGGTGTGAA-3′ and reverse 5′-ACAGTCTTCTGGGTGGCAGT-3′), and the AluS expression level was calculated using the 2−ΔΔCt method.

Determination of the Alu methylation levels and patterns in LCLs

Genomic DNA was isolated from LCLs of ASD individuals and sex-/age-matched unaffected controls using the GENEzol Reagent (Geneaid, Taiwan) according to the manufacturer’s recommended protocol. The DNA concentration was determined using a NanoDrop 1000 spectrophotometer (Thermo Scientific, USA). COBRA is designed to determine methylation levels and patterns of two CpG loci within AluS subfamilies with the highest number of copies and CpG loci in the human genome (Fig. 2a). Briefly, 1 μg of genomic DNA from each sample was treated with sodium bisulfite using EZ DNA Methylation-Gold™ Kit (Zymo, Irving, CA, USA). The bisulfite-treated DNA was subjected to 45 cycles of PCR (Hot Start Taq DNA polymerase, QIAGEN, USA) with two specific primers for AluS subfamilies, AluS-F (5′-GGRGRGGTGGTTTARGTTTGTAA-3′) and AluS-R (5′-CTAACTTTTTATATTTTTAATAAAAACRAAATTTCACCA-3′), at 53 °C for annealing. Then, the AluS amplicons were digested with TaqI restriction enzyme (Thermo Scientific, USA) in TaqI buffer and incubated at 65 °C overnight. Finally, the digested products were electrophoresed on an 8% non-denaturing polyacrylamide gel, and band intensities were determined for assessment of the AluS methylation levels and patterns as described previously [26].

Fig. 2
figure2

Alu element structure and illustration of COBRA for determining AluS methylation levels and patterns. a Alu elements are approximately 300 bp in length and have a dimeric structure that is separated by an A-rich region (A5TACA6) and ends with a poly-A tail. The left half of the Alu contains the A and B boxes, which are internal promoters for RNA polymerase III. b Illustration of the COBRA method designed to assess methylation of two CpGs at the internal promoter of AluS subfamilies. The four different methylation patterns of AluS were calculated from the percentages of differently digested products of 133, 75, 58, 43, and 32 bp. c Representative gel image from the COBRA for AluS subfamilies

Full size image

Upon gel electrophoresis, the digested Alu-amplicons were resolved into six fragments of 133, 90, 75, 58, 43, and 32 bp, which represented different methylation states (Fig. 2b, c), including the uCuC methylation state (represented by the 133-bp fragment). The mCuC methylation state was represented by the 90-bp fragment. The uCmC methylation state was represented by the 75- and 58-bp fragments. The mCmC methylation state was represented by the 43- and 32-bp fragments. The calculation for percent AluS methylation was performed as follows. First, the percentage of band intensity was divided by the length (bp) of each DNA fragment: %133/133 = A, %58/58 = B, %75/75 = C, %90/90 = D, %43/43 = E, and %32/32 = F. The percentages of Alu methylation levels and patterns were then calculated using the following formulas: percentage methylated loci (%mC) = 100 × (E + B)/(2A + E + B + C + D), percentage of hypermethylated pattern (%mCmC) = 100 × F/(A + C + D + F), percentage of partially methylated pattern (%uCmC) = 100 × C/(A + C + D + F), percentage of partially methylated pattern (%mCuC) = 100 × D/(A + C + D + F), and percentage of partially hypomethylated pattern (%uCuC) = 100 × A/(A + C + D + F).

Statistical analyses

DEGs were determined using two-tailed t tests with adjusted Bonferroni correction, and adjusted P values less than 0.05 were considered significant. Fisher’s exact test with Benjamini-Hochberg (FDR = 0.05) correction was used to investigate the association between the DEG lists and the human Alu-inserted gene lists. P values less than 0.05 were considered significant. Pathway and function analyses were performed with IPA using Fisher’s exact test with Benjamini-Hochberg correction for multiple testing (FDR = 0.05); P values less than 0.05 were considered significant. Two-tailed t tests with Benjamini-Hochberg correction (FDR = 0.05) were used to analyze the differences in AluS methylation and expression level; the adjusted P value threshold was 0.05.

Results

Genes containing Alu elements are differentially expressed in ASD blood based on the integration of data from multiple studies

We hypothesized that Alu elements are associated with altered gene expression in the peripheral blood and blood-derived cell lines of ASD individuals and provided CpG sites for DNA methylation within Alu-inserted genes. To test this hypothesis, we obtained five transcriptome profiles from ASD studies available in GEO DataSets and identified a list of DEGs from each study. We then overlapped these DEG lists with the human Alu-inserted genes using Fisher’s exact test; the results are shown in Table 3. These results showed that the DEGs in the peripheral blood of ASD individuals were significantly associated with the Alu-inserted gene lists. For the group of all Alu insertion types, the lists of DEGs (up- and downregulated genes) were significantly associated with Alu elements. These lists consist of 388, 869, 1001, and 1492 DEGs from the GSE6575, GSE25507, GSE42133, and GSE18123 studies, respectively. We subsequently identified the up- and downregulated DEGs and assessed their overlap with the “all Alu insertion” gene list type; the results showed a strong association for downregulated genes (adjusted P value < 0.0005). However, the “all Alu insertion” type was only weakly associated with upregulated genes (adjusted P value = 0.015) in the GSE42133 study and showed a weaker association in the other studies. DEGs were also compared with the other types of Alu insertion gene lists, including the exonic, exonized, intronic, and promoter types. The results are shown in Table 3. These results showed that the DEGs were more strongly associated with the Alu intronic lists than with the other types of Alu insertions. Intronic Alu elements were also strongly associated with the downregulated gene lists. Moreover, because the GSE15402 study has also reduced the heterogeneity of ASD by subgrouping ASD individuals based on their clinical phenotypes using supervised and unsupervised clustering analyses of ADI-R scores, we further investigated whether Alu insertions were associated with DEGs in each ASD subgroup. Interestingly, we found that genes with Alu insertion were significantly associated with DEGs in ASD subgroups but not when all ASD cases were combined (Additional file 2).

Table 3 Association analyses between the differentially expressed genes (DEGs) in ASD and the human Alu-inserted gene lists
Full size table

We subsequently selected the lists of downregulated genes that were significantly associated with all Alu insertion lists from each study. These lists were used to identify overlapping genes among the selected studies using Venn diagram analysis (Fig. 3). Interestingly, the diagram revealed the reproducibility of Alu-inserted genes that were differentially expressed in ASD whole blood/blood cells. Significant genes that were identified in at least two studies (320 overlapping genes; see Additional file 3) were selected to identify gene regulatory networks, diseases, and biological functions that were associated with Alu elements in ASD.

Fig. 3
figure3

Venn diagram of genes containing Alu that are differentially expressed in ASD. The significant DEGs with Alu insertions from each study based on Fisher’s exact test overlapped. The diagram shows the reproducibility of Alu-inserted genes that were differentially expressed in peripheral blood and blood-derived cell lines from ASD individuals. A total of 320 genes were selected to identify biological functions and gene regulatory networks through an Ingenuity Pathway Analysis (IPA)

Full size image

Biological functions, canonical pathways, and gene regulatory networks of dysregulated genes containing Alu elements are associated with neurological functions and neurodevelopmental disorders

IPA was used to predict the diseases, biological functions, canonical pathways, and gene regulatory networks that were associated with Alu elements in ASD. The results showed that the overlapping genes (320 genes) were significantly associated with neurological diseases and nervous system development and function (adjusted Fisher’s exact test, P value < 0.05) (Table 4). Interestingly, 21 overlapping genes were associated with autism or intellectual disability. Eighteen and 20 of the overlapping genes were also associated with mental retardation and cognitive impairment, respectively, which were classified as neurodevelopmental disorders comorbid with ASD. The IPA results also revealed significant canonical pathways associated with the overlapping genes (Table 5), including neurotrophin/TRK signaling, ERK/MAPK signaling, axonal guidance signaling, CREB signaling in neurons, estrogen receptor signaling, androgen signaling, IL-6 signaling, and neuroinflammation signaling, all of which have been associated with ASD.

Table 4 Diseases and biological functions associated with reproducible DEGs with Alu insertion predicted by the Ingenuity Pathway Analysis (IPA)
Full size table
Table 5 Canonical pathways associated with reproducible DEGs with Alu insertion predicted by the Ingenuity Pathway Analysis (IPA)
Full size table

Gene regulatory networks are sets of genes that interact to control a specific cell function, including cell differentiation, metabolism, cell cycle, and response to environmental cues. A representative gene network associated with neurological disease (Fig. 4) revealed the interaction of overlapping genes with additional molecules from the IPA database involved in estrogen receptor and androgen signaling, axonal guidance signaling, CREB signaling in neurons, and neurotrophin/TRK signaling. This network suggested that some ASD-related biological pathways are associated with Alu-inserted DEGs in the blood of ASD individuals.

Fig. 4
figure4

Predicted gene regulatory network of the overlapping genes associated with neurological disease. This network revealed interactions or relationships among the overlapping molecules (gray background) and with other molecules from the IPA database (white background) that play a role in several mechanisms associated with neurological disease and estrogen receptor and androgen signaling, which is known to be associated with sex bias in ASD (labeled pink)

Full size image

AluS methylation levels and patterns in LCLs

To determine the levels of DNA methylation within Alu elements, which contributes to a large proportion of the CpG sites in the human genome, we used LCLs as a model and conducted a COBRA. Specifically, a COBRA of AluS subfamilies was performed to determine the percentages of AluS methylation levels (mC) and patterns (mCmC, mCuC, uCmC, uCuC) in LCLs, and the results showed that comparisons of the AluS methylation levels and patterns between the ASD (n = 36) and sex- and age-matched unaffected control (n = 20) groups were not significantly different (Table 6). The levels and patterns of Alu methylation in individual LCLs are shown in Additional file 4. Due to heterogeneity within ASD individuals, our LCLs were categorized into three subgroups, namely, savant (S), mild (M), and language-impaired (L), based on the ADI-R interview scores according to Hu and Steinberg [14]. In addition, Hu and colleagues reported that the gene expression profiles of LCLs were significantly different among these ASD subgroups [15].

Table 6 COBRA-derived percentages of AluS methylation and patterns in LCLs from ASD individuals and sex- and age-matched controls
Full size table

We then compared the AluS methylation levels and patterns between each ASD subgroup and the sex- and age-matched groups. Interestingly, the results showed significant AluS methylation patterns associated with specific ASD subgroups. The percentage of the partially methylated pattern uCmC (20.06% ± 0.92%, adjusted P value = 0.043) was significantly increased in ASD subgroup M (Table 6, Fig. 5). In addition, the percentage of the partially methylated pattern mCuC (20.43% ± 1.17%, adjusted P value = 0.010) was significantly decreased in ASD subgroup S compared with sex- and age-matched controls. There were no significant patterns of AluS methylation in ASD subgroup L. These findings suggested that methylation of Alu elements might play a role in AluS expression and/or transcriptional profile of some ASD subgroups but not all ASD individuals.

Fig. 5
figure5

Box plot of the Alu methylation patterns in the LCLs of ASD subgroup M. In ASD subgroup M, the percentage of the partially methylated pattern uCmC (20.06% ± 0.92%) was significantly increased. In ASD subgroup S, the partially methylated pattern mCuC was significantly decreased. *adjusted P value < 0.05

Full size image

Quantitative reverse transcription-PCR analysis of the AluS expression levels in 56 LCLs

To understand the influence of DNA methylation within AluS elements on AluS expression levels, a quantitative RT-PCR analysis was performed using the LCLs of ASD and sex- and age-matched unaffected control groups. The results showed no significant differences in AluS expression between the ASD and the matched control groups (fold change (FC) = 1.75, adjusted P value 0.316, Table 7), which was similar to the results obtained from the AluS methylation level and pattern analyses. However, comparisons of the AluS expression levels between the ASD phenotypic subgroups and sex- and age-matched control groups revealed that these were significantly different. Specifically, the AluS expression levels were significantly decreased in ASD subgroup L (FC = 0.29, adjusted P value = 0.032) and significantly increased in ASD subgroup S (FC = 3.68, adjusted P value = 0.038). However, in ASD subgroup M, the AluS expression level was not significantly different.

Table 7 Quantitative RT-PCR analyses of AluS expression levels in the LCLs of ASD and control groups
Full size table

Correlation analysis between AluS methylation levels and AluS expression levels

To demonstrate that AluS methylation regulates the expression of AluS elements, we correlated the DNA methylation and expression levels of AluS subfamilies in the LCLs from ASD individuals and sex- and age-matched controls. The correlation analyses revealed non-linear relationships between DNA methylation and the expression of AluS subfamilies in all groups (ASD + control, Fig. 6). However, when analyzed within the ASD subgroups, we found that the partially methylated pattern uCmC showed a moderate positive relationship with AluS expression in ASD subgroup M (coefficient r = 0.5149, Fig. 7), whereas other AluS methylation patterns did not correlate with ASD subgroup M. Similarly, DNA methylation and expression of AluS subfamilies were not correlated in ASD subgroups L and S (Additional files 5 and 6). These findings suggested that the partially methylated pattern uCmC might regulate the expression of AluS elements in ASD subgroup M.

Fig. 6
figure6

Correlation analysis between AluS methylation and expression level for all LCL samples. The AluS expression for each LCL was normalized with the average GAPDH dCt of the control group. The Alu expression levels were then calculated using the 2−ΔΔCt method

Full size image
Fig. 7
figure7

Correlation analysis between AluS methylation and expression levels in ASD subgroup M and sex- and age-matched controls. The AluS expression of each LCL was normalized to the average GAPDH dCt of the control group. The Alu expression levels were then calculated using the 2−ΔΔCt method

Full size image

Discussion

The mechanisms through which human Alu elements are involved in ASD are unclear. In a previous study, Mbarek et al. found a polymorphism of an Alu element in the blood of ASD individuals with more severe clinical symptoms [27]. This allele is located in intron 27b of the neurofibromatosis type 1 (NF1) gene, which is reported to be an ASD risk gene. However, substantially more information was needed to conclude that Alu elements play important roles in ASD. Our findings provide preliminary information regarding the association between DEGs and Alu elements in the peripheral blood samples of ASD individuals. It is important to note that Alu elements have the potential to influence gene expression through insertion into the gene structure and through the contribution of gene regulatory elements, such as transcription binding sites and CpG sites [17]. Moreover, Alu elements can act as enhancers, alternate promoters, transcription start sites, and inhibitors of transcription via heterochromatin formation [28, 29].

Several studies have reported altered gene expression profiles in the peripheral blood or blood-derived cell lines of ASD individuals [15, 30,31,32,33]. However, no studies have identified an association between DEGs in ASD and Alu elements that could influence gene expression. We first performed multiple gene expression profile comparisons with human Alu-inserted genes in ASD samples using Fisher’s exact test. Interestingly, we found that the Alu-inserted genes were associated with DEGs in ASD individuals in four studies (Table 3). These results indicated a strong association of Alu insertion with downregulated genes in ASD, suggesting that Alu elements could affect gene downregulation. Another interesting question is whether the specific Alu-inserted positions are associated with gene regulation. For this analysis, the Alu-inserted gene list was categorized into four types of Alu insertions: intronic, exonized, exonic, and promoter inserts. Interestingly, we found that intronic Alu insertion was significantly associated with DEGs in ASD invididuals. According to a previous study by Tsirigos et al., a genome-wide computational analysis showed that Alu elements selectively retained in the intronic region of inserted genes were associated with specific functions, including the regulation of transcription, RNA processing and splicing, and translation [34]. However, it remains unclear how intronic Alu insertions can influence the expression of genes associated with ASD-related biological functions, and this topic should be investigated further in future research.

We then explored the biological functions associated with 320 DEGs with Alu insertions that overlapped between the selected studies. These overlapping genes were significantly associated with autism or intellectual disability and ASD co-morbid disorders, including mental retardation and cognitive impairment. It is noteworthy that sex hormone signaling pathways, including estrogen receptor signaling and androgen signaling pathways, were also associated with DEGs containing Alu elements. ASD is biased towards males with a male-to-female ratio of at least 4:1, and there is accumulating evidence showing that sex hormones and related pathways may play an important role in the sex bias of ASD [35,36,37]. Recent studies have found that high testosterone exposure during pregnancy significantly correlates with the development of ASD, social development, and language development [38,39,40]. Moreover, the expression levels of estrogen receptor mRNA and protein were found to be reduced in the brain tissues and serum of ASD individuals [41, 42]. Interestingly, one of the DEGs with Alu insertion related to estrogen receptor and androgen signaling is NCOA1, which is a gene encoding the nuclear receptor coactivator 1 protein. This coregulatory protein has been reported to interact with estrogen receptor and androgen receptor, which oppositely regulate the transcription of the RORA gene encoding the retinoic acid-related (RAR) orphan receptor-alpha (RORA) protein [36]. RORA is a hormone-dependent transcription factor that regulates many genes, including CYP19A1 (aromatase), an enzyme that converts testosterone to estrogen [37, 43,44,45]. Our previous study also found that RORA binds NCOA1 and regulates the transcription of the CYP19A1 gene [36]. Thus, it is possible that Alu elements are associated with the sex bias of ASD by disrupting the androgen receptor/estrogen receptor-mediated regulation of RORA and CYP19A1. Other interesting pathways that are known to be associated with ASD include neurotrophin/TRK signaling [46], ERK/MAPK signaling [47], axonal guidance signaling [48], CREB signaling in neurons [49], IL-6 signaling [50], and neuroinflammation signaling [51]. Although these biological functions and pathways associated with DEGs in ASD with Alu element insertion have been strongly implicated in ASD, we cannot exclude the possibility that all Alu-inserted genes might be enriched for these functions or pathways, regardless of whether they exhibit differential expression in ASD. The full list of Alu-inserted genes could have been used for IPA analysis as a control to address this question. However, IPA does not allow the use of a gene list larger than 8000 genes. This issue should be investigated in future research studies using other pathway analysis programs.

To date, several DMVs at specific CpG sites or DMRs have been identified in many ASD tissue types, including LCLs [9], whole blood [10], and brain [11]. However, these studies did not cover noncoding regions, including repetitive sequences and retrotransposons. DNA methylation is the major epigenetic mechanism that represses retrotransposons in the human genome, particularly Alu elements due to their relatively high CpG density [18, 52, 53]. Kochanek and colleagues found that DNA methylation at CpG sites within an internal promoter (B box) of Alu elements could inhibit their transcriptional activity [54]. In this study, we also assessed the DNA methylation and expression level of AluS subfamilies that contributed the most Alu element copies and CpG sites in the human genome [18]. COBRA of AluS was designed to measure CpG methylation at the internal promoter, including the B box, in LCLs from ASD patients and sex- and age-matched controls. The LCLs were representative samples from three phenotypic groups based on previous multivariate cluster analyses of ADI-R scores of 1954 individuals with ASD. Our results revealed that the AluS methylation levels and patterns in the LCLs of the combined group of ASD individuals compared with sex- and age-control groups were not significantly different. However, when LCLs were divided into phenotypic subgroups, the AluS methylation patterns in two of the ASD subgroups were significantly different compared with those in the sex- and age-matched controls. In the ASD subgroup with mild symptoms (M), the AluS methylation pattern uCmC was increased compared with that in the sex- and age-matched controls. Furthermore, the AluS methylation pattern mCuC was reduced in the ASD subgroup with savant skill (S) compared with the sex- and age-matched controls. Although the abovementioned methylation patterns are statistically significant within a specific ASD subtype, it is also worth noting that other AluS methylation patterns were not much altered in ASD in comparison with control groups. One possible explanation is that AluS elements may be dysregulated only at certain sites in the genome, such as in exonic regions. However, the COBRA analysis used in this study does not measure methylation of AluS at specific genomic locations. Thus, it is likely that signals from significant methylation sites were dampened by noise from a large number of Alu elements located in other non-significant sites in the genome.

We then measured the expression of AluS in the LCLs from the subgroups of ASD and found that AluS was over-expressed in ASD subgroup S and downregulated in the ASD subgroup with severe language impairment (L) relative to their respective sex- and age-matched controls. To investigate whether DNA methylation patterns regulated the expression of Alu elements, we analyzed the correlations between AluS methylation and expression levels. A positive correlation between AluS methylation and expression was only observed in ASD subgroup M and in sex- and age-matched controls in which the methylation pattern uCmC showed a moderate correlation with AluS expression. It is interesting to note that the methylation pattern uCmC was also significantly increased in subgroup M in the COBRA analysis. These findings suggest that AluS expression may be regulated by the AluS methylation pattern uCmC. However, such relationships were not found when all ASD individuals were combined or in other ASD subgroups. In ASD subgroup S, the AluS methylation patterns and expression levels were not correlated, although the AluS methylation pattern mCuC was significantly reduced, and AluS expression was over-expressed. Moreover, ASD subgroup L showed a significant reduction in AluS expression, but not DNA methylation, indicating that the decrease in AluS expression in ASD subgroup L might result from other epigenetic or gene regulatory mechanisms, such as the disrupted transcription of genes containing an Alu element by altered transcription factor binding at the promoter of the genes but not at the Alu promoter or the other gene regulatory mechanism that have been implicated in LCLs of ASD individuals [55]. These findings suggest that the AluS regulatory mechanism related to methylation might be unique for a specific ASD subpopulation, and there might be other regulatory mechanisms involved in the regulation of AluS methylation and expression, which should be investigated in future research.

The changes of Alu RNAs in LCLs from some ASD subgroups might reflect molecular function because Alu RNAs are involved in transcriptome diversity by contributing recognition sites for RNA editing and alternative splicing within Alu sequences [56]. Altered DNA methylation within Alu elements may have a negative impact on gene regulatory networks, which in turn, may affect the biological functions associated with ASD (Fig. 8). The exact mechanisms remain unclear, but it is possible that altered methylation within the Alu promoter, together with other gene regulatory mechanisms, may lead to changes in Alu transcript expression and retrotransposition. Because Alu elements contain transcription factor binding sites within their sequences, these elements could disrupt gene structure and functions by serving as an alternative transcription start site, enhancer, or promoter for cis-/trans-regulation of the target genes once inserted into a new genomic location by the retrotransposition [28, 29, 57,58,59].

Fig. 8
figure8

Schematic diagram illustrating a possible mechanism of Alu elements in ASD. Our model suggests that exposure to environmental factors or dysregulation of other DNA methylation regulatory mechanisms lead to changes in CpG methylation patterns in Alu elements. Such changes alter transcription factor binding and, possibly in combination with other Alu regulatory mechanisms, cause the dysregulation of the expression and retrotransposition of Alu elements. Disrupted Alu retrotransposition results in changes in target genes via cis-/trans-regulatory mechanisms, which, in turn, dysregulate gene expression and gene regulatory networks known to be negatively impacted in ASD

Full size image

We obtained transcriptomic data from five different studies that used different sample types (e.g., whole blood, leukocytes, and LCLs) as well as different cohorts of subjects with ASD. It remains unclear whether ASD gene expression or methylation signatures will be similar across tissue sources. There might be variations in gene expression and methylation patterns among different sample types and among different cells in the same tissue type. This issue is a potential limitation of this study and should be investigated further in future research. In addition, some of the five selected transcriptome studies did not subgroup ASD individuals before performing DEG analysis and there were no ADI-R scores available for meta-analysis, whereas others used different criteria and strategies to divide ASD into subgroups. It is therefore difficult to relate our methylation analysis findings to previous analyses of Alu-inserted DEGs in which all ASD individuals were combined or studies using different subgrouping strategies. However, we performed Fisher’s exact tests using the microarray data from a previous study (GSE15402) that used the same subgrouping criteria and the same cell model (LCLs) used in this study. Interestingly, we found that genes with Alu element insertion were associated with specific subgroups of ASD rather than with all individuals with ASD. The results of Alu methylation analysis using COBRA and Alu expression analysis using qRT-PCR also showed that the dysregulation of Alu methylation and Alu expression were observed only when ASD cases were divided into subgroups. This finding suggests that the subgrouping of ASD individuals will help reduce heterogeneity and may lead to the discovery of novel mechanisms associated with Alu element in ASD subpopulations. A better understanding of the molecular mechanisms specific for each ASD subgroup might allow the identification of biomarkers and treatment strategies personalized to each subgroup in the future. However, the sample size of this study is relatively small, especially when individuals with ASD were further divided into subgroups. The role of Alu in the context of ASD deserves further studies using a larger cohort. Most importantly, more experiments are needed to verify that Alu elements play important roles in ASD. Brain tissue samples may be used to confirm these findings, and single-cell RNA sequencing might help reduce cell heterogeneity.

Conclusion

Our findings show that the DEGs in males with ASD from several studies with different sample groups are associated with human Alu elements. Differentially expressed genes with Alu insertions were associated with neurodevelopmental disorders and neurological functions involved in the etiology of ASD. In particular, genes involved in estrogen receptor and androgen signaling pathways, which have been reported to be related with sex bias in ASD, were identified as DEGs with Alu insertions. In addition, the global methylation of AluS subfamilies in LCLs was investigated and revealed different AluS methylation patterns within specific ASD phenotypic subgroups. Our findings suggest that the classification of ASD patients into subgroups based on the clinical or behavioral phenotypes of individual patients might help improve our understanding of ASD etiology by reducing the inherent heterogeneity within the ASD population. In addition, this study provides suggestive evidence for an association between Alu elements and ASD.

Abbreviations

ADI-R:

Autism Diagnostic Interview-Revised

ASD:

Autism spectrum disorder

COBRA:

Combined bisulfite restriction analysis

DEGs:

Differentially expressed genes

DMRs:

Differentially methylated regions

DMVs:

Differentially methylated variants

DSM-5:

The Diagnostic and Statistical Manual of Mental Disorders, Fifth Edition

ERK:

Extracellular signal-regulated kinases

GEO:

Gene Expression Omnibus

IPA:

Ingenuity Pathway Analysis

LCLs:

Lymphoblastoid cell lines

MAPK:

Mitogen-activated protein kinases

RORA:

Retinoic acid receptor-related orphan receptor alpha

RPMI:

Roswell Park Memorial Institute

RT-PCR:

Reverse transcription-polymerase chain reaction

SINEs:

Short interspersed elements

TEs:

Transposed elements

TRK:

Tyrosine kinases

References

  1. 1.

    American Psychiatric Association: Diagnostic and Statistical Manual of Mental Disorders. 5th edition. Arlington, VA: 2013.

  2. 2.

    Christensen DL, Baio J, Braun KV, Bilder D, Charles J, Constantino JN, Daniels J, Durkin MS, Fitzgerald RT, Kurzius-Spencer M, et al. Prevalence and characteristics of autism spectrum disorder among children aged 8 years—autism and developmental disabilities monitoring network, 11 sites, United States, 2012. Morbidity Mortal Weekl Repo Surveillance Summ (Washington, DC : 2002). 2016;65(3):1–23.

  3. 3.

    Lavelle TA, Weinstein MC, Newhouse JP, Munir K, Kuhlthau KA, Prosser LA. Economic burden of childhood autism spectrum disorders. Pediatrics. 2014;133(3):e520–9.

  4. 4.

    Schaefer GB, Mendelsohn NJ. Genetics evaluation for the etiologic diagnosis of autism spectrum disorders. Genet Med. 2008;10(1):4–12.

  5. 5.

    Moosa A, Shu H, Sarachana T, Hu VW. Are endocrine disrupting compounds environmental risk factors for autism spectrum disorder? Horm Behav. 2017. https://doi.org/10.1016/j.yhbeh.2017.10.003.

  6. 6.

    Tordjman S, Somogyi E, Coulon N, Kermarrec S, Cohen D, Bronsard G, Bonnot O, Weismann-Arcache C, Botbol M, Lauth B, et al. Gene x environment interactions in autism spectrum disorders: role of epigenetic mechanisms. Front Psychiatry. 2014;5:53.

  7. 7.

    Bailey A, Le Couteur A, Gottesman I, Bolton P, Simonoff E, Yuzda E, Rutter M. Autism as a strongly genetic disorder: evidence from a British twin study. Psychol Med. 1995;25(1):63–77.

  8. 8.

    Hallmayer J, Cleveland S, Torres A, Phillips J, Cohen B, Torigoe T, Miller J, Fedele A, Collins J, Smith K, et al. Genetic heritability and shared environmental factors among twin pairs with autism. Arch Gen Psychiatry. 2011;68(11):1095–102.

  9. 9.

    Nguyen A, Rauch TA, Pfeifer GP, Hu VW. Global methylation profiling of lymphoblastoid cell lines reveals epigenetic contributions to autism spectrum disorders and a novel autism candidate gene, RORA, whose protein product is reduced in autistic brain. FASEB J. 2010;24(8):3036–51.

  10. 10.

    Wong CC, Meaburn EL, Ronald A, Price TS, Jeffries AR, Schalkwyk LC, Plomin R, Mill J. Methylomic analysis of monozygotic twins discordant for autism spectrum disorder and related behavioural traits. Mol Psychiatry. 2014;19(4):495–503.

  11. 11.

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

  12. 12.

    Ladd-Acosta C, Hansen KD, Briem E, Fallin MD, Kaufmann WE, Feinberg AP. Common DNA methylation alterations in multiple brain regions in autism. Mol Psychiatry. 2014;19(8):862–71.

  13. 13.

    Nardone S, Sams DS, Reuveni E, Getselter D, Oron O, Karpuj M, Elliott E. DNA methylation analysis of the autistic brain reveals multiple dysregulated biological pathways. Transl Psychiatry. 2014;4:e433.

  14. 14.

    Hu VW, Steinberg ME. Novel clustering of items from the Autism Diagnostic Interview-Revised to define phenotypes within autism spectrum disorders. Autism Res. 2009;2(2):67–77.

  15. 15.

    Hu VW, Sarachana T, Kim KS, Nguyen A, Kulkarni S, Steinberg ME, Luu T, Lai Y, Lee NH. Gene expression profiling differentiates autism case-controls and phenotypic variants of autism spectrum disorders: evidence for circadian rhythm dysfunction in severe autism. Autism Res. 2009;2(2):78–97.

  16. 16.

    Batzer MA, Deininger PL. Alu repeats and human genomic diversity. Nat Rev Genet. 2002;3(5):370–9.

  17. 17.

    Polak P, Domany E. Alu elements contain many binding sites for transcription factors and may play a role in regulation of developmental processes. BMC Genomics. 2006;7:133.

  18. 18.

    Luo Y, Lu X, Xie H. Dynamic Alu methylation during normal development, aging, and tumorigenesis. Biomed Res Int. 2014;2014:784706.

  19. 19.

    Baccarelli A, Bollati V. Epigenetics and environmental chemicals. Curr Opin Pediatr. 2009;21(2):243–51.

  20. 20.

    Patchsung M, Settayanon S, Pongpanich M, Mutirangura D, Jintarith P, Mutirangura A. Alu siRNA to increase Alu element methylation and prevent DNA damage. Epigenomics. 2018;10(2):175–85.

  21. 21.

    Barrett T, Troup DB, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, Marshall KA, Phillippy KH, Sherman PM, et al. NCBI GEO: archive for functional genomics data sets—10 years on. Nucleic Acids Res. 2011;39(Database issue):D1005–10.

  22. 22.

    Edgar R, Domrachev M, Lash AE. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002;30(1):207–10.

  23. 23.

    Levy A, Sela N, Ast G. TranspoGene and microTranspoGene: transposed elements influence on the transcriptome of seven vertebrates and invertebrates. Nucleic Acids Res. 2008;36(Database issue):D47–52.

  24. 24.

    Bao W, Kojima KK, Kohany O. Repbase update, a database of repetitive elements in eukaryotic genomes. Mob DNA. 2015;6:11.

  25. 25.

    Saeed AI, Sharov V, White J, Li J, Liang W, Bhagabati N, Braisted J, Klapa M, Currier T, Thiagarajan M, et al. TM4: a free, open-source system for microarray data management and analysis. BioTechniques. 2003;34(2):374–8.

  26. 26.

    Tiwawech D, Srisuttee R, Rattanatanyong P, Puttipanyalears C, Kitkumthorn N, Mutirangura A. Alu methylation in serum from patients with nasopharyngeal carcinoma. Asian Pac J Cancer Prev. 2014;15(22):9797–800.

  27. 27.

    Mbarek O, Marouillat S, Martineau J, Barthelemy C, Muh JP, Andres C: Association study of the NF1 gene and autistic disorder. Am J Med Genet 1999, 88(6):729-732.

  28. 28.

    Elbarbary RA, Lucas BA, Maquat LE. Retrotransposons as regulators of gene expression. Science (New York, NY). 2016;351(6274):aac7247.

  29. 29.

    Ichiyanagi K. Epigenetic regulation of transcription and possible functions of mammalian short interspersed elements, SINEs. Genes Genet Syst. 2013;88(1):19–29.

  30. 30.

    Alter MD, Kharkar R, Ramsey KE, Craig DW, Melmed RD, Grebe TA, Bay RC, Ober-Reynolds S, Kirwan J, Jones JJ, et al. Autism and increased paternal age related changes in global levels of gene expression regulation. PLoS One. 2011;6(2):e16715.

  31. 31.

    Gregg JP, Lit L, Baron CA, Hertz-Picciotto I, Walker W, Davis RA, Croen LA, Ozonoff S, Hansen R, Pessah IN, et al. Gene expression changes in children with autism. Genomics. 2008;91(1):22–9.

  32. 32.

    Kong SW, Collins CD, Shimizu-Motohashi Y, Holm IA, Campbell MG, Lee IH, Brewster SJ, Hanson E, Harris HK, Lowe KR, et al. Characteristics and predictive value of blood transcriptome signature in males with autism spectrum disorders. PLoS One. 2012;7(12):e49475.

  33. 33.

    Pramparo T, Lombardo MV, Campbell K, Barnes CC, Marinero S, Solso S, Young J, Mayo M, Dale A, Ahrens-Barbeau C, et al. Cell cycle networks link gene expression dysregulation, mutation, and brain maldevelopment in autistic toddlers. Mol Syst Biol. 2015;11(12):841.

  34. 34.

    Tsirigos A, Rigoutsos I. Alu and b1 repeats have been selectively retained in the upstream and intronic regions of genes of specific functional classes. PLoS Comput Biol. 2009;5(12):e1000610.

  35. 35.

    Baron-Cohen S, Auyeung B, Norgaard-Pedersen B, Hougaard DM, Abdallah MW, Melgaard L, Cohen AS, Chakrabarti B, Ruta L, Lombardo MV. Elevated fetal steroidogenic activity in autism. Mol Psychiatry. 2015;20(3):369–76.

  36. 36.

    Sarachana T, Hu VW. Differential recruitment of coregulators to the RORA promoter adds another layer of complexity to gene (dys) regulation by sex hormones in autism. Mol Autism. 2013;4(1):39.

  37. 37.

    Sarachana T, Xu M, Wu RC, Hu VW. Sex hormones in autism: androgens and estrogens differentially and reciprocally regulate RORA, a novel candidate gene for autism. PLoS One. 2011;6(2):e17116.

  38. 38.

    Knickmeyer R, Baron-Cohen S, Raggatt P, Taylor K, Hackett G. Fetal testosterone and empathy. Horm Behav. 2006;49(3):282–92.

  39. 39.

    Tordjman S, Ferrari P, Sulmont V, Duyme M, Roubertoux P. Androgenic activity in autism. Am J Psychiatry. 1997;154(11):1626–7.

  40. 40.

    Whitehouse AJ, Maybery MT, Hart R, Mattes E, Newnham JP, Sloboda DM, Stanley FJ, Hickey M. Fetal androgen exposure and pragmatic language ability of girls in middle childhood: implications for the extreme male-brain theory of autism. Psychoneuroendocrinology. 2010;35(8):1259–64.

  41. 41.

    Crider A, Thakkar R, Ahmed AO, Pillai A. Dysregulation of estrogen receptor beta (ERbeta), aromatase (CYP19A1), and ER co-activators in the middle frontal gyrus of autism spectrum disorder subjects. Mol Autism. 2014;5(1):46.

  42. 42.

    Altun H, Kurutas EB, Sahin N, Sinir H, Findikli E. Decreased levels of G protein-coupled estrogen receptor in children with autism spectrum disorders. Psychiatry Res. 2017;257:67–71.

  43. 43.

    Hu VW, Nguyen A, Kim KS, Steinberg ME, Sarachana T, Scully MA, Soldin SJ, Luu T, Lee NH. Gene expression profiling of lymphoblasts from autistic and nonaffected sib pairs: altered pathways in neuronal development and steroid biosynthesis. PLoS One. 2009;4(6):e5775.

  44. 44.

    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.

  45. 45.

    Sarachana T, Hu VW. Genome-wide identification of transcriptional targets of RORA reveals direct regulation of multiple genes associated with autism spectrum disorder. Mol Autism. 2013;4(1):14.

  46. 46.

    Correia CT, Coutinho AM, Sequeira AF, Sousa IG, Lourenco Venda L, Almeida JP, Abreu RL, Lobo C, Miguel TS, Conroy J, et al. Increased BDNF levels and NTRK2 gene association suggest a disruption of BDNF/TrkB signaling in autism. Genes Brain Behav. 2010;9(7):841–8.

  47. 47.

    Mukaetova-Ladinska EB, Arnold H, Jaros E, Perry R, Perry E. Depletion of MAP2 expression and laminar cytoarchitectonic changes in dorsolateral prefrontal cortex in adult autistic individuals. Neuropathol Appl Neurobiol. 2004;30(6):615–23.

  48. 48.

    Suda S, Iwata K, Shimmura C, Kameno Y, Anitha A, Thanseem I, Nakamura K, Matsuzaki H, Tsuchiya KJ, Sugihara G, et al. Decreased expression of axon-guidance receptors in the anterior cingulate cortex in autism. Mol Autism. 2011;2(1):14.

  49. 49.

    Ebert DH, Greenberg ME. Activity-dependent neuronal signalling and autism spectrum disorder. Nature. 2013;493(7432):327–37.

  50. 50.

    Wei H, Zou H, Sheikh AM, Malik M, Dobkin C, Brown WT, Li X. IL-6 is increased in the cerebellum of autistic brain and alters neural cell adhesion, migration and synaptic formation. J Neuroinflammation. 2011;8:52.

  51. 51.

    El-Ansary A, Al-Ayadhi L. Neuroinflammation in autism spectrum disorders. J Neuroinflammation. 2012;9:265.

  52. 52.

    Liu WM, Maraia RJ, Rubin CM, Schmid CW. Alu transcripts: cytoplasmic localisation and regulation by DNA methylation. Nucleic Acids Res. 1994;22(6):1087–95.

  53. 53.

    Xing J, Hedges DJ, Han K, Wang H, Cordaux R, Batzer MA. Alu element mutation spectra: molecular clocks and the effect of DNA methylation. J Mol Biol. 2004;344(3):675–82.

  54. 54.

    Kochanek S, Renz D, Doerfler W. Transcriptional silencing of human Alu sequences and inhibition of protein binding in the box B regulatory elements by 5′-CG-3′ methylation. FEBS Lett. 1995;360(2):115–20.

  55. 55.

    Sarachana T, Zhou R, Chen G, Manji HK, Hu VW. Investigation of post-transcriptional gene regulatory networks associated with autism spectrum disorders by microRNA expression profiling of lymphoblastoid cell lines. Genome Med. 2010;2(4):23.

  56. 56.

    Sorek R, Ast G, Graur D. Alu-containing exons are alternatively spliced. Genome Res. 2002;12(7):1060–7.

  57. 57.

    Deininger P. Alu elements: know the SINEs. Genome Biol. 2011;12(12):236.

  58. 58.

    Deininger PL, Batzer MA. Alu repeats and human disease. Mol Genet Metab. 1999;67(3):183–93.

  59. 59.

    Hasler J, Strub K. Alu elements as regulators of gene expression. Nucleic Acids Res. 2006;34(19):5491–7.

Download references

Acknowledgements

TSae, CT, and ST are graduate students in the M.Sc. Program in Clinical Biochemistry and Molecular Medicine, Faculty of Allied Health Sciences, Chulalongkorn University. This work is part of TSae’s thesis research to be presented in partial fulfillment of the requirements for the M.Sc. degree. This research was supported by the Ratchadaphiseksomphot Endowment Fund part of the “Research Grant for New Scholar CU Researcher’s Project” (GDNS 57-034-37-002) and the Faculty of Allied Health Sciences Research Fund (AHS-CU 58003) to TS. TSae, CT, and ST were financially supported by The 90th Anniversary Chulalongkorn University Fund (Ratchadaphiseksomphot Endowment Fund: TSae-GCUGR1125601055M 53-1; CT-GCUGR1125601058M 53-4; ST-GCUGR1125601056M 53-2). TSae and ST received additional financial support from “The Scholarship from the Graduate School, Chulalongkorn University to commemorate the 72nd anniversary of His Majesty King Bhumibala Aduladeja,” while CT received the Research Assistant Scholarship, Chulalongkorn University. VWH is supported by NIEHS grant R21 ES023061. AM is supported by the Thailand Research Fund (DPG5980005). KS is supported by the Thailand Research Fund (BRG5980001). We wish to thank Dr. Charoenchai Puttipanyalears for teaching TSae and CT to conduct the COBRA. We also gratefully acknowledge the resources provided by the Autism Genetic Resource Exchange (AGRE) Consortium* and the participating AGRE families. The Autism Genetic Resource Exchange is a program of Autism Speaks and is supported, in part, by grant 1U24MH081810 from the National Institute of Mental Health to Clara M. Lajonchere (PI). *The AGRE Consortium: Dan Geschwind, M.D., Ph.D., UCLA, Los Angeles, CA; Maja Bucan, Ph.D., University of Pennsylvania, Philadelphia, PA; W.Ted Brown, M.D., Ph.D., F.A.C.M.G., N.Y.S. Institute for Basic Research in Developmental Disabilities, Long Island, NY; Rita M. Cantor, Ph.D., UCLA School of Medicine, Los Angeles, CA; John N. Constantino, M.D., Washington University School of Medicine, St. Louis, MO; T.Conrad Gilliam, Ph.D., University of Chicago, Chicago, IL; Martha Herbert, M.D., Ph.D., Harvard Medical School, Boston, MA; Clara Lajonchere, Ph.D, Cure Autism Now, Los Angeles, CA; David H. Ledbetter, Ph.D., Emory University, Atlanta, GA; Christa Lese-Martin, Ph.D., Emory University, Atlanta, GA; Janet Miller, J.D., Ph.D., Cure Autism Now, Los Angeles, CA; Stanley F. Nelson, M.D., UCLA School of Medicine, Los Angeles, CA; Gerard D. Schellenberg, Ph.D., University of Washington, Seattle, WA; Carol A. Samango-Sprouse, Ed.D., George Washington University, Washington, D.C.; Sarah Spence, M.D., Ph.D., UCLA, Los Angeles, CA; Matthew State, M.D., Ph.D., Yale University, New Haven, CT. Rudolph E. Tanzi, Ph.D., Massachusetts General Hospital, Boston, MA.

Funding

This research was supported by the Ratchadaphiseksomphot Endowment Fund part of the “Research Grant for New Scholar CU Researcher’s Project” (GDNS 57-034-37-002) and the Faculty of Allied Health Sciences Research Fund (AHS-CU 58003) to TS. TSae, CT, and ST were financially supported by The 90th Anniversary Chulalongkorn University Fund (Ratchadaphiseksomphot Endowment Fund: TSae-GCUGR1125601055M 53-1; CT-GCUGR1125601058M 53-4; ST-GCUGR1125601056M 53-2). TSae and ST received additional financial support from “The Scholarship from the Graduate School, Chulalongkorn University to commemorate the 72nd anniversary of His Majesty King Bhumibala Aduladeja,” while CT received the Research Assistant Scholarship, Chulalongkorn University. VWH is supported by NIEHS grant R21 ES023061. AM is supported by the Thailand Research Fund (DPG5980005). KS is supported by the Thailand Research Fund (BRG5980001).

Availability of data and materials

The datasets generated and/or analyzed during the current study are available from the corresponding author on reasonable request.

Author information

Affiliations

  1. M.Sc. Program in Clinical Biochemistry and Molecular Medicine, Department of Clinical Chemistry, Faculty of Allied Health Sciences, Chulalongkorn University, Bangkok, Thailand
    • Thanit Saeliw
    • , Chayanin Tangsuwansri
    •  & Surangrat Thongkorn
  2. Maximizing Thai Children’s Developmental Potential Research Unit, Department of Pediatrics, Faculty of Medicine, Chulalongkorn University and King Chulalongkorn Memorial Hospital, the Thai Red Cross Society, Bangkok, Thailand
    • Weerasak Chonchaiya
  3. Center of Excellence for Medical Genetics, Department of Pediatrics, Faculty of Medicine, Chulalongkorn University, Bangkok, Thailand
    • Kanya Suphapeetiporn
  4. Excellence Center for Medical Genetics, King Chulalongkorn Memorial Hospital, Thai Red Cross Society, Bangkok, Thailand
    • Kanya Suphapeetiporn
  5. Center of Excellence in Molecular Genetics of Cancer and Human Diseases, Department of Anatomy, Faculty of Medicine, Chulalongkorn University, Bangkok, Thailand
    • Apiwat Mutirangura
  6. Age-related Inflammation and Degeneration Research Unit, Department of Clinical Chemistry, Faculty of Allied Health Sciences, Chulalongkorn University, 154 Soi Chula 12, Rama 1 Road, Wangmai, Pathumwan, Bangkok, 10330, Thailand
    • Tewin Tencomnao
    •  & Tewarit Sarachana
  7. Department of Biochemistry and Molecular Medicine, The George Washington University School of Medicine and Health Sciences, The George Washington University, Washington, DC, USA
    • Valerie W. Hu
Authors
  1. Search for Thanit Saeliw in:
  2. Search for Chayanin Tangsuwansri in:
  3. Search for Surangrat Thongkorn in:
  4. Search for Weerasak Chonchaiya in:
  5. Search for Kanya Suphapeetiporn in:
  6. Search for Apiwat Mutirangura in:
  7. Search for Tewin Tencomnao in:
  8. Search for Valerie W. Hu in:
  9. Search for Tewarit Sarachana in:

Contributions

TSae performed all the experiments, analyzed the data, and drafted the manuscript under the supervision of TS, AM, TT, WC, and KS. CT assisted TSae in performing COBRA and real-time RT-PCR experiments. ST conducted LCL culture and genomic DNA and RNA extraction. AM and TT provided the enzymes and reagents required for the COBRA analysis. VWH provided all LCLs and the Ingenuity Pathway Analysis program. TS conceived of the study, designed the experiments, and participated in the writing and editing of this manuscript. All the authors read and approved the final manuscript.

Corresponding author

Correspondence to Tewarit Sarachana.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Demographic information of LCLs used in this study. (DOC 73 kb)

Additional file 2:

Association analysis between the DEGs in ASD from GSE15402 and the human Alu-inserted genes lists when ASD individuals were sub-grouped based on ADI-R scores. The comparisons were performed for five types of Alu insertions, including exonized, exonic, intronic, promoter and all types. The Fisher’s exact test P value and number of DEGs are shown. (DOC 46 kb)

Additional file 3:

List of the overlapping genes. (DOC 118 kb)

Additional file 4:

Levels and patterns of Alu methylation in each individual. (DOC 98 kb)

Additional file 5:

Correlation analysis between AluS methylation and expression level of ASD subgroup L. (PDF 329 kb)

Additional file 6:

Correlation analysis between AluS methylation and expression level of ASD subgroup S. (PDF 333 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

Verify currency and authenticity via CrossMark

Cite this article

Saeliw, T., Tangsuwansri, C., Thongkorn, S. et al. Integrated genome-wide Alu methylation and transcriptome profiling analyses reveal novel epigenetic regulatory networks associated with autism spectrum disorder. Molecular Autism 9, 27 (2018). https://doi.org/10.1186/s13229-018-0213-9

Download citation

Keywords

  • Autism spectrum disorder
  • Alu elements
  • Retrotransposon
  • DNA methylation
  • Epigenetic regulation
  • Gene expression profiles
  • Subgrouping
  • Lymphoblastoid cell lines
  • Sex bias
  • Neuroinflammation

Molecular Autism

ISSN: 2040-2392