Abstract
Autism spectrum disorder (ASD) is a complex disorder that is underpinned by numerous dysregulated biological pathways, including pathways that affect mitochondrial function. Epigenetic mechanisms contribute to this dysregulation and DNA methylation is an important factor in the etiology of ASD. We measured DNA methylation of peroxisome proliferator-activated receptor-gamma coactivator-1 alpha (PGC-1α), as well as five genes involved in regulating mitochondrial homeostasis to examine mitochondrial dysfunction in an ASD cohort of South African children. Using targeted Next Generation bisulfite sequencing, we found differential methylation (p < 0.05) at six key genes converging on mitochondrial biogenesis, fission and fusion in ASD, namely PGC-1α, STOML2, MFN2, FIS1, OPA1, and GABPA. PGC-1α, the transcriptional regulator of biogenesis, was significantly hypermethylated at eight CpG sites in the gene promoter, one of which contained a putative binding site for CAMP response binding element 1 (CREB1) (p = 1 × 10–6). Mitochondrial DNA (mtDNA) copy number, a marker of mitochondrial function, was elevated (p = 0.002) in ASD compared to controls and correlated significantly with DNA methylation at the PGC-1α promoter and there was a positive correlation between methylation at PGC-1α CpG#1 and mtDNA copy number (Spearman’s r = 0.2, n = 49, p = 0.04) in ASD. Furthermore, DNA methylation at PGC-1α CpG#1 and mtDNA copy number correlated significantly (p < 0.05) with levels of urinary organic acids associated with mitochondrial dysfunction, oxidative stress, and neuroendocrinology. Our data show differential methylation in ASD at six key genes converging on PGC-1α-dependent regulation of mitochondrial biogenesis and function. We demonstrate that methylation at the PGC-1α promoter is associated with elevated mtDNA copy number and metabolomic evidence of mitochondrial dysfunction in ASD. This highlights an unexplored role for DNA methylation in regulating specific pathways involved in mitochondrial biogenesis, fission and fusion contributing to mitochondrial dysfunction in ASD.
Introduction
Autism spectrum disorder (ASD) is defined by the presence of behavioral traits (Lord et al., 2020) despite being a highly heritable neurodevelopmental disorder (Sandin et al., 2017). ASD is characterized by deficits in social communication and restrictive, repetitive behaviors (American Psychiatric Association [APA], 2013). ASD is a complex disorder that affects the central nervous system as well as other organ systems, implying the dysregulation of pleiotropic biological and developmental pathways. ASD is underpinned by a heterogeneous genetic architecture that includes rare de novo genetic variations and low risk, common single nucleotide mutations (Gaugler et al., 2014; Feliciano et al., 2019). There is increasing evidence for the role of DNA methylation in modulating ASD phenotypes. This is evident from discordant identical ASD twin studies (Wong et al., 2014), studies using brain tissue from individuals with ASD (Ladd-Acosta et al., 2014; Nardone et al., 2014), with recent reviews collating numerous reports on ASD epigenetics (Rylaarsdam and Guemez-Gamboa, 2019; Tremblay and Jiang, 2019; Wiśniowiecka-Kowalnik and Nowakowska, 2019; Yoon et al., 2020).
Given the varied phenotypes and co-morbidities observed in individuals with ASD, numerous and diverse biological pathways have been implicated in ASD etiology. These include gene regulatory-, signaling-, synaptic-, and mitochondrial-pathways (Frye and Rossignol, 2011; Voineagu, 2012; Ayhan and Konopka, 2019; Kumar et al., 2019). Mitochondrial dysfunction is emerging as a key contributor to ASD etiology; deficiencies in oxidative phosphorylation (OXPHOS) can decrease the production of ATP, which is an essential requirement for brain function and neurodevelopment (Valenti et al., 2014). The observation that congenital errors of mitochondrial metabolism contribute to >5% of ASD cases (Manzi et al., 2008) first implied a role for mitochondrial dysfunction in ASD. This has since been supported by clinical (Oliveira et al., 2005), biochemical (Melnyk et al., 2012; Tang et al., 2013), molecular (Dhillon et al., 2011) and more recently, epigenetic data (Stathopoulos et al., 2020). Citrigno et al. (2020) comprehensively reviewed the recent experimental data that support a role for mitochondrial dysfunction in ASD.
Mitochondrial homeostasis is dynamic and is regulated by interdependent pathways that govern mitochondrial biogenesis, mitophagy, mitochondrial fission and fusion (Ventura-Clapier et al., 2008; Ploumi et al., 2017). These dynamic mechanisms enable cellular adaptation to change energy demands, nutrient availability, and oxidative stress by regulating mitochondrial DNA (mtDNA) copy number (Gaziev et al., 2014; Santos et al., 2014). Abnormal or fluctuating levels of mtDNA copy number is a marker of mitochondrial dysfunction (Malik and Czajka, 2013; Sun et al., 2019). Oxidative stress and mitochondrial dysfunction can lead to increased mtDNA copy number as a compensatory mechanism to maintain cellular homeostasis (Lee et al., 2000; Michel et al., 2012; Petersen et al., 2014; Kim et al., 2019). An essential transcriptional regulator of mitochondrial homeostasis is peroxisome proliferator-activated receptor-gamma coactivator-1 alpha (PGC-1α) which regulates fatty acid β-oxidation, OXPHOS, gluconeogenesis and antioxidant defense responses (Barone et al., 2019). PGC-1α catalyzes mitochondrial biogenesis by upregulating nuclear respiratory factors 1 (Nrf-1) and 2 (Nrf-2), also known as GA binding protein transcription factor subunit alpha (GABPA), which promotes the transcription of mitochondrial transcription factors A (TFAM) and B2 (TFB2M) (Ventura-Clapier et al., 2008; Ploumi et al., 2017). Importantly, numerous studies report a correlation between DNA methylation of the PGC-1α promoter and PGC-1α transcription, mtDNA copy number and metabolic disease (Barrès et al., 2009; Sookoian et al., 2010; Heinonen et al., 2015; Kresovich et al., 2018), suggesting that DNA methylation regulates PGC-1α-driven mitochondrial biogenesis.
Similarly, mitochondrial fission and fusion maintain mitochondrial homeostasis by limiting intracellular reactive oxygen species (ROS; Karbowski and Youle, 2003). Mitochondrial fusion occurs in response to mild oxidative stress and is mediated by mitofusins 1 and 2 (MFN1 and MFN2) and optic atrophy 1 (OPA1) (Chen et al., 2003; Duvezin-Caubet et al., 2006; Ishihara et al., 2006). These proteins work in conjunction with accessory proteins, such as stomatin-like protein 2 (STOML2), which maintains the long isoforms of OPA1 needed to fuse the inner mitochondrial membranes (Tondera et al., 2009). Fission occurs during severe oxidative stress and separates damaged mitochondrial components from the healthy mitochondrial network. Fission is mediated by dynamin-related protein 1 (DRP1) which works with mitochondrial fission protein 1 (FIS1) to divide the outer mitochondrial membrane (Smirnova et al., 2001; Stojanovski et al., 2004).
Recently, genes involved in mitochondrial biogenesis, fission and fusion have been implicated in neuropathology. PGC-1α is reported to play an important role in excitatory neurotransmitter signaling, neuroprotection, neuroinflammation and neurogenesis and has been implicated in bipolar disorder, Parkinson’s disease, Huntington’s disease, schizophrenia, and Alzheimer’s disease (Bu et al., 2017). In addition, recent studies have implicated PGC-1α (Barone et al., 2019), fusion and fission genes (Carrasco et al., 2019; Pecorelli et al., 2020) in ASD etiology (Barone et al., 2019), while STOML2 was the most significantly differentially methylated (DM) gene in a South African ASD population (Stathopoulos et al., 2020).
Although there has been a massive increase in ASD genetic data, it has been primarily from ASD samples from Northern hemisphere countries. Our recent study identified DM genes in a South African cohort that were enriched for mitochondrial pathways, supporting a role for epigenetic dysregulation of energy metabolism in ASD. In view of our previous data (Stathopoulos et al., 2020) and the evidence from literature implicating both differential DNA methylation and mitochondrial dysfunction in ASD, we examined the relationship between these two processes in South African children with ASD. First, we determined whether PGC-1α was DM between ASD and matched controls from a South African population. Secondly, we measured DNA methylation of five additional genes that are central to mitochondrial biogenesis (GABPA), fission (FIS1) and fusion (MFN2, STOML2, and OPA1). Subsequently, we addressed whether DNA methylation modified mitochondrial function, which was measured using mtDNA copy number and urinary metabolomics. Overall, our aim was to test the hypothesis that DNA methylation of mitochondrial biogenesis genes changes mitochondrial function, thereby contributing to the etiology of ASD.
Materials and Methods
Participants and Sample Collection
Given that South Africa has no centralized ASD network with an associated repository of DNA, RNA and/or ASD phenotypes, we recruited children with ASD and age- and gender-matched typically developing controls (Stathopoulos et al., 2020). We recruited both male and female participants (6–17 years), but only pre-pubertal boys (6–12 years old) were used given that we recruited fewer than five females. We screened 145 participants (93 ASD; 52 controls) from three demographic groups: African-, European-, and Mixed-ancestry. An Autism Diagnostic Observation Schedule, Second Edition (ADOS-2) assessment was used to phenotype the ASD group, as well as to ensure the absence of ASD traits in the control group. The study protocol was approved by the University of Cape Town, as well as the Western Cape Government approval to recruit participants at schools. Buccal cells were collected from participants for DNA extraction; this is a minimally invasive collection method suited for DNA methylation studies (Lowe et al., 2013; Berko et al., 2014; Smith et al., 2015). Participants who did not consent to providing buccal samples or those samples which yielded poor DNA quality were excluded from the study. Urine samples were collected for organic acids extraction for metabolomic analysis using gas chromatography–mass spectrometry (GC–MS).
Targeted Next-Generation Bisulfite Sequencing
DNA was extracted from buccal cells as previously described (Stathopoulos et al., 2020). DNA methylation was quantified using targeted Next-Generation Bisulfite Sequencing (tNGBS). After excluding participants who would not provide a DNA sample or whose buccal cells yielded DNA of poor quality, we proceeded to DNA methylation quantification using tNGBS. First, we quantified methylation for PGC-1α and STOML2, both of which were DM in our previous discovery cohort, in a larger validation cohort of ASD (n = 55) and controls (n = 44). We also quantified methylation for additional genes required for mitochondrial fusion and fission FIS1, MFN2, OPA1, and GABPA in a subset of samples (ASD, n = 22; controls n = 22). The tNGBS was completed by EpigenDx, Inc. (MA, United States) who designed a total of 32 tNGBS assays to analyze 171 CpG sites across six genes. tNGBS assays were designed by (i) obtaining and annotating each gene sequence using the Ensembl genome browser, (ii) re-evaluating target EPIC array probe sequences against the UCSC genome browser to identify LINE, SINE, LTR elements and other DNA repeat sequences; (iii) excluding sequences containing repetitive elements, low sequences complexity, high thymidine content and overall CpG density, and (iv) identifying assays that passed PCR optimization. This process was used to designed seven assays that analyzed 26 CpG sites for PGC-1α, eight assays that analyzed 45 CpG sites for STOML2, four assays covering 30 CpG sites for FIS1, five assays covering 26 CpG sites for MFN2, five assays covering 25 CpG sites for OPA1 and three assays that analyzed 19 CpG sites for GABPA. Methylation levels were calculated by dividing the number of methylated reads by the number of total reads. Unpaired two-tailed t-tests with unequal variance were used to determine the significantly DM CpG sites between ASD and control (p < 0.05).
Mitochondrial DNA Copy Number and Deletion
Mitochondrial DNA copy number and mitochondrial deletions were measured using multiplex real-time quantitative polymerase chain reaction (RT-qPCR) in 108 participants (ASD n = 68; controls n = 40). Mitochondrial genes mitochondrially encoded NADH: ubiquinone oxidoreductase core subunit 1 (MT-ND1) and mitochondrially encoded NADH: ubiquinone oxidoreductase core subunit 4 (MT-ND4), were amplified by RT-qPCR and normalized to the nuclear gene, beta-2-microglobulin (B2M), in the same PCR reaction. The probes were coupled to non-fluorescent quenchers (BHQ®, LGC BioSearch) and both the primers and probes used were previously reported by Grady et al. (2014). Each DNA sample (20 ng/μl) was amplified in triplicate, with standard curves set for each gene using equimolar pooled DNA from ASD and controls in a 10-fold dilution series. A two-step thermal profile was used with denaturation at 95°C for 10 min, followed by 40 cycles of 10 s at 95°C, 30 s at 60°C on the Rotor-Gene Q 6-plex (QIAGEN). The raw data was checked for the presence of outliers that were the result of poor qPCR detection. Samples with a triplicate threshold cycle Ct standard deviation greater than 1 were removed before analysis. A total of 99 samples (ASD n = 59; controls n = 40) were used in subsequent analyzes. Mitochondrial copy number was calculated using the equation 2−ΔΔCt, where ΔCt = Ct (MT-ND1) - Ct (B2M). Mitochondrial deletion was calculated using the equation 2−ΔΔCt, where ΔCt = Ct (MT-ND4) - Ct (MT-ND1). Significance was determined using a two-tailed unpaired t-test p < 0.05 was considered statistically significant.
Metabolomic Correlation Analysis
We examined the correlation between PGC-1α methylation, mtDNA copy number and levels of urinary organic acids associated with mitochondrial dysfunction. Urinary organic acids for 35 participants (ASD n = 21 and controls n = 13) had been previously extracted and quantified by GC–MS (Stathopoulos et al., 2020) to screen for a panel of metabolites that were indicative of mitochondrial respiratory disease in South African children (Reinecke et al., 2012). The GC–MS data had been analyzed and deconvoluted using a standard metabolomics-based data processing workflow (Reinecke et al., 2012) and was log 2 transformed before statistical analysis. The Shapiro-Wilks test was used to test for normality, after which the Spearman Rank Correlation analysis was performed in IBM SPSS Statistics (v26).
Statistical Analysis
Data analysis was conducted in Excel (2019) and IBM SPSS Statistics (v26). Scatter plots were generated in Spyder (Python 3.7) (Raybaut, 2009) using the matplotlib package (Hunter, 2007), REL: v3.4.2 (Caswell et al., 2021). All tests with p value less than 0.05 were considered statistically significant. For detailed statistical analyzes for each dataset, see the relevant sections in “Materials and Methods” and “Results.”
Results
Autism Spectrum Disorder Cohort Phenotype and Demography
The individuals in our study with ASD spanned the full range of developmental phenotypes observed in ASD. This is reflected in both the number of different Autism Diagnostic Observation Schedule, Second Edition (ADOS-2) Modules used for assessments and the autism severity scores (Supplementary Table 1). Each ADOS-2 Module is tailored to match different developmental levels, ranging from pre-verbal individuals with ASD to those with fluent speech. Although our cohort comprised of different demographic groups, demography did not correlate with any ASD trait or any molecular marker (data not shown).
PGC-1α Promoter Is Hypermethylated in ASD
Our previous study implicated mitochondrial dysfunction in our South African ASD cohort, therefore we investigated whether DNA methylation contributed specifically to the regulation of mitochondrial biogenesis by measuring the methylation of PGC-1α, a central transcriptional regulator of mitochondrial biogenesis, in a larger sample of participants. We defined highly variable CpGs as those sites where the methylation range exceeded 5% across all samples to identify functionally significant DM genes. This threshold is consistent with in vitro (Ziller et al., 2013) and in vivo (D’Aquila et al., 2020) methylation studies. There were 12 highly variable CpG sites in PGC-1α that were significantly DM (p < 0.05) between ASD and controls (ASD n = 55, controls n = 43) (Supplementary Table 2). Of these, eight CpG sites were hypermethylated in ASD and clustered around the transcriptional start site (TSS), between the 5′ untranslated region (UTR) and intron 1 (Figures 1A,B), while four sites, located at intron 2, intron 12, and the 3′UTR were hypomethylated in ASD (Supplementary Figure 1).
FIGURE 1
The significant DM sites of PGC-1α in our study are consistent with a role for DNA methylation in regulating mitochondrial biogenesis. To examine this hypothesis, we quantified DNA methylation of GABPA, the transcriptional regulator of mitochondrial biogenesis that acts directly downstream of PGC-1α, as well as four genes involved in mitochondrial fission and fusion (STOML2, MFN2, OPA1, and FIS1) in a subset of our cohort (n = 22 ASD, n = 22 controls). STOML2 contained two DM CpG sites, located in intron 2 and exon 5 downstream of the TSS; these sites were hypermethylated in ASD (Figure 2 and Supplementary Table 2). Significant DM sites between ASD and controls were also identified at one CpG site in GABPA and FIS1, and at two sites in MFN2 and OPA1 (Figure 2 and Supplementary Table 2). Therefore, we observed multiple DM mitochondrial biogenesis, fission and fusion genes converging on the regulation of mitochondrial homeostasis in ASD (Figure 3), congruent with a role for DNA methylation in the dysregulation of mitochondrial function in our cohort.
FIGURE 2
FIGURE 3
PGC-1α Hypermethylation Is Associated With Altered Mitochondrial Function
We examined whether methylation affected mitochondrial biogenesis and function because DNA methylation was altered at several key regulators of mitochondrial biogenesis in our ASD cohort. In silico transcription factor binding site analysis of DM CpG sites in the PGC-1α promoter predicted a putative binding site for the transcription factor CAMP response binding element 1 (CREB1) at the CpG#1site (p = 1 × 10–6). This suggests that the DM site in ASD may have functional significance, which we examined by assessing mitochondrial function in our cohort. This was done by quantifying mtDNA copy number and deletions, as well as the levels of urinary metabolites typically associated with mitochondrial disease in ASD relative to controls.
Mitochondrial DNA Copy Number Is Increased in ASD
Increased mtDNA copy number is a compensatory response to mild oxidative stress (Gaziev et al., 2014; Santos et al., 2014) and is an established biomarker of mitochondrial function (Castellani et al., 2020). mtDNA copy number of MT-ND1 relative to B2M was significantly elevated in the ASD compared to the controls (p = 0.002) (Figure 4A). Mitochondrial deletions are typical of mitochondrial disease, therefore we examined whether mitochondrial deletions differed in ASD compared to controls by quantifying the copy number of the mitochondrial gene MT-ND1 relative to MT-ND4; the latter resides in the major mitochondrial deletion arc. While mtDNA deletions were not significantly elevated in ASD relative to controls (p = 0.162), we observed markedly elevated mtDNA deletions in some ASD individuals (Figure 4B). Notably, mtDNA copy number correlated significantly (Spearman’s r = 0.9, n = 49, p = 8.814 × 10–10) with mtDNA deletions (Figure 4C), which suggests that elevated mtDNA copy number is associated with mitochondrial dysfunction in our cohort.
FIGURE 4
PGC-1α Promoter Methylation Correlates Significantly With mtDNA Copy Number and This Relationship Is Altered in ASD
We examined whether PGC-1α promoter methylation is associated with mtDNA copy number and/or deletions, and thus mitochondrial function in ASD. We observed a significant correlation between DNA methylation at the PGC-1α promoter and mtDNA copy number, and this relationship differed between ASD and control groups. In the control group, PGC-1α methylation at CpG#4 correlated negatively with both mtDNA copy number (Spearman’s r = −0.4, n = 42, p = 0.045) and mtDNA deletions (Spearman’s r = −0.4, n = 42, p = 0.032; Figure 4D). However, in the ASD group, there was a significant positive correlation between PGC-1α methylation at CpG#1 and mtDNA copy number (Spearman’s r = 0.9, n = 49, p = 0.04; Figure 4D), with no correlation between PGC-1α methylation and mtDNA deletions. This suggests that PGC-1α methylation is associated with mitochondrial biogenesis and function and that this relationship is disrupted in our ASD group.
DNA Methylation and mtDNA Copy Number Are Associated With Metabolomic Markers of Mitochondrial Dysfunction
We investigated whether the differential methylation of PGC-1α and elevated mtDNA copy number observed in ASD were associated with metabolomic evidence of mitochondrial dysfunction. We examined the correlation between PGC-1α methylation, mtDNA copy number and levels of urinary organic acids which had previously been measured (ASD n = 20, controls n = 13) using GC-MS (Stathopoulos et al., 2020). DNA methylation at PGC-1α CpG#1, which is DM in our samples, correlated significantly with three of the 55 urinary organic acids tested (Supplementary Table 3B). MtDNA copy number correlated with 22 urinary metabolites associated with mitochondrial dysfunction (Supplementary Table 3A). Notably, both PGC-1α methylation and mtDNA copy number are associated with metabolites derived from Branched-chain amino acid (BCAA) metabolism, including 3-hydroxy-3-methylglutaric acid (3-H-3-MGA), which correlated significantly with both PGC-1α CpG#1 (p = 0.043) and with mtDNA copy number (p = 0.005) (Figures 5A,B). In addition, mtDNA copy number is associated most significantly (r < −0.3 or >0.3; p < 0.01) with metabolites derived from four metabolic pathways: fatty acid oxidation, phenylalanine-dopamine synthesis pathway, glycine-glutamine metabolism, and BCAAs (Supplementary Table 3A). These pathways converged on mitochondrial OXPHOS, one-carbon metabolism and neurotransmitter synthesis (Figure 6). Our data are consistent with an established metabolomic model for altered mitochondrial metabolism and neuroendocrinology (Yu et al., 2016) and supports an association between DNA methylation, elevated mtDNA copy number and mitochondrial dysfunction in ASD.
FIGURE 5
FIGURE 6
Discussion
Autism spectrum disorder is a heritable, complex phenotype with numerous molecular pathways contributing to its etiology (Sandin et al., 2017). Despite the high heritability of ASD, there is no single or simple genetic mutation that accounts for ASD, and the disorder is characterized by phenotypic and clinical heterogeneity. This implies that epigenetic mechanisms may be important in ASD, and DNA methylation is known to contribute to ASD etiology (Wiśniowiecka-Kowalnik and Nowakowska, 2019). DNA methylation dysregulates many biological pathways in ASD, including, but not limited to, immune function (Nardone et al., 2017), chromatin remodeling (Andrews et al., 2018), synaptic signaling and neuronal regulation (Wong et al., 2019). Increasingly, mitochondrial biological pathways have also been implicated in ASD (Oliveira et al., 2005; Dhillon et al., 2011; Tang et al., 2013). The link between mitochondrial dysfunction and ASD is unsurprising given that efficient ATP production is essential for brain development and function. However, the relationship between DNA methylation and mitochondrial function is not fully understood.
Our data examined this relationship by testing the hypothesis that mitochondrial biogenesis and fusion genes are DM between ASD and controls and that this differential methylation affects mitochondrial function in ASD. Previously, we performed a whole-epigenome DNA methylation screen which identified 898 DM genes in a cohort of South African children (n = 48) which converged on nine canonical pathways involved in mitochondrial metabolism (Stathopoulos et al., 2020). Targeted next-generation bisulfite sequencing of a subset of genes and DNA pyrosequencing of two DM genes confirmed the DM observed in the whole-epigenome screen. Urinary metabolomic data in a subset of this cohort supported our hypothesis of mitochondrial dysfunction in ASD. Given this association between DM and general mitochondrial metabolism (Stathopoulos et al., 2020), this paper focused on the role of DNA methylation in regulating mitochondrial biogenesis in a larger South African ASD cohort by examining the link between PGC-1α methylation, mtDNA copy number and metabolomic evidence of mitochondrial dysfunction.
The transcriptional regulator of biogenesis, PGC-1α, was significantly DM between ASD and controls in our cohort, with the promoter region being hypermethylated in ASD. The promoter region included a CpG site (CpG#1) containing a putative transcription binding site for CREB1, which is a potent activator of PGC-1α transcription (Wu et al., 2006). Although we were not able to directly examine PGC-1α transcription in our cohort, this CREB1 site is reported to be DM in metabolic disease (Sookoian et al., 2010) which suggests that DM sites at the PGC-1α promoter region and TSS in our ASD cohort could affect gene transcription and subsequently, mitochondrial biogenesis and function.
Mitochondria are adaptive to changing cellular metabolic demands, thus they are dynamic organelles that are remodeled by biogenesis, fission, and fusion (Wai and Langer, 2016). Therefore, we measured the methylation of additional genes involved in mitochondrial biogenesis, fission and fusion and found that GABPA, which facilitates mtDNA replication downstream of PGC-1α in ASD, was DM (Ploumi et al., 2017). We also found that genes involved in mitochondrial fission (STOML2, MFN2, and OPA1) and fusion (FIS1) were DM in ASD. Mitochondrial fusion is induced in response to ROS stress to maintain optimal mitochondrial function and homeostasis (Karbowski and Youle, 2003). Collectively, these DM genes converge on the pathways regulating mitochondrial homeostasis in response to metabolic and oxidative stress (Figure 3). Of note, STOML2 was hypermethylated at two CpG sites downstream of the TSS in ASD and it was the most significantly DM gene in our ASD cohort in our earlier study (Stathopoulos et al., 2020). The STOML2 protein plays an important role in mitochondrial fusion by stabilizing the OPA1 protein, which facilitates fusion of the inner mitochondrial membranes (Duvezin-Caubet et al., 2006; Ishihara et al., 2006; Tondera et al., 2009). STOML2-deficient cells fail to undergo mitochondrial fusion during stress, leading to mitochondrial fragmentation (Tondera et al., 2009). STOML2 is also well-established as an anti-apoptotic gene in cancer cells, highlighting the importance of fusion to re-establish mitochondrial homeostasis and prevent mitophagy (autophagy of mitochondria) under stress. Our data is consistent with previous work showing that both fission and fusion genes are differentially expressed in ASD (Tang et al., 2013; Carrasco et al., 2019; Pecorelli et al., 2020). This supports the link between ASD and mitochondrial fusion and fission which highlights the differential methylation of mitochondrial genes on an integrated pathway level.
Our data show genes that regulate mitochondrial biogenesis, fission and fusion are DM in our cohort. The relationship between DM and mRNA expression is complex; hyper- or hypo-methylation can either increase or decrease gene expression depending on the gene as well as where in the gene region the methylation occurs (Arechederra et al., 2018; Tremblay and Jiang, 2019; Rauluseviciute et al., 2020). The low integrity of RNA extracted from the buccal cells used in our study meant we were not able to determine how the DM we observed impacted gene expression in our cohort. However, previous studies have shown that DM alters the expression of these genes (Tang et al., 2013; Carrasco et al., 2019; Pecorelli et al., 2020; Yang et al., 2020). Additionally, we investigated whether mtDNA copy number, a marker of mitochondrial function (Castellani et al., 2020), was altered in our cohort. Changes in mtDNA copy number have been reported in ASD, with both increases (Gu et al., 2013; Chen et al., 2015; Yoo et al., 2017) and decreases (Valiente-Pallejà et al., 2018; Singh et al., 2020) observed in ASD. These discrepancies can be attributed to several factors that differed across studies, including the age of the participants studied, the presence of co-morbidities and the degree of mitochondrial dysfunction in participants. We observed a significant increase in mtDNA copy number in our ASD group compared to controls, which represents an established compensatory mechanism in response to mitochondrial dysfunction (Ventura-Clapier et al., 2008). We also observed a significant positive correlation of mtDNA copy number with mtDNA deletions in ASD, suggesting that elevated mtDNA copy number is indicative of mitochondrial dysfunction in our cohort. Increased mtDNA copy number is observed as a response to oxidative stress in animal models and in vitro studies (Al-Kafaji and Golbahar, 2013; Masser et al., 2016; Zeng et al., 2019), as well as in human clinical studies using buccal samples (Kolbasina et al., 2020). This compensatory mechanism has also been reported in mitochondrial diseases (Bai et al., 2004; Grady et al., 2014; Thompson et al., 2016), neuropsychiatric and neurodevelopmental disorders (Montier et al., 2009; Giulivi et al., 2010; Chung et al., 2019; Chen et al., 2020; Ridout et al., 2020).
A correlation between PGC-1α methylation and mtDNA copy number is established with numerous studies reporting that PGC-1α promoter methylation is consistently associated with reduced mRNA expression and lower mtDNA copy numbers in a range of different tissues (Lillycrop et al., 2005; Ling et al., 2008; Barrès et al., 2009; Sookoian et al., 2010; Chen et al., 2012; Gillberg et al., 2014; Heinonen et al., 2015; Kelstrup et al., 2016; Kresovich et al., 2018; Liu et al., 2019; Petrie et al., 2020; Yang et al., 2020). These include studies that focused on both neurological (Yang et al., 2020) and metabolic disorders (Sookoian et al., 2010; Zhao et al., 2017; Kresovich et al., 2018). Consistent with this, we report a negative correlation between PGC-1α promoter hypermethylation and both mtDNA copy number and deletion in our control group. Thus, hypomethylation of the PGC-1α promoter is associated with increased mtDNA copy number and deletion. This is consistent with hypomethylation leading to an upregulation of PGC-1α transcription, and subsequently, PGC-1α-dependent mitochondrial biogenesis under conditions of mild mitochondrial dysfunction.
However, this relationship between PGC-1α promoter methylation and mtDNA copy number was perturbed in ASD. We observed a significant positive correlation between PGC-1α promoter methylation at CpG#1 and mtDNA copy number. This is consistent with observed PGC-1α promoter methylation being associated with the mtDNA copy number alterations seen in our ASD cohort. However, given that PGC-1α methylation is associated with reduced mRNA expression, this does not necessarily mean that PGC-1α methylation directly alters mtDNA copy number in ASD by dysregulating mitogenesis. Instead, it shows that hypermethylation of PGC-1α is associated with a dysregulation of mitochondrial homeostasis in our ASD cohort. This is consistent with the systemic metabolic dysfunction associated with ASD, and the central role played by PGC-1α in regulating mitochondrial metabolism and homeostasis. Hypermethylation at the PGC-1α promoter could inhibit the PGC-1α-dependent activation of antioxidant genes (St-Pierre et al., 2006) which is congruent with the evidence for elevated oxidative stress in ASD (Bjørklund et al., 2020). Oxidative stress can induce mitochondrial biogenesis via nuclear factor erythroid 2-related factor 2 (Nfe2l2) (Erlank et al., 2011) which upregulates GABPA independently of PGC-1α in a redox-sensitive manner (Kolbasina et al., 2020). We found no relationship between PGC-1α methylation and mtDNA deletion in ASD, suggesting an absence of adaptive hypomethylation to compensate for metabolic and putative oxidative stress.
To further explore the link between differential methylation, mtDNA copy number and mitochondrial function, we used metabolomics analysis, which directly reflects the biochemical activity, including mitochondrial activity, of a biological sample (Glinton and Elsea, 2019). Our previous study examined which metabolites were significantly altered in ASD compared to controls, using a panel of urinary organic acids that are indicative of mitochondrial disease in South African children (Reinecke et al., 2012). Here, we used this metabolomic data in conjunction with the new data generated in this study to investigate whether PGC-1α methylation or mtDNA copy number correlated with metabolomic evidence of mitochondrial dysfunction in our cohort (Stathopoulos et al., 2020).
We found that mtDNA copy number was associated with a metabolomic profile that was consistent with a link between DNA methylation, mitochondrial dysfunction, and neuropathology. Metabolites derived from tyrosine, tryptophan, and glycine significantly (p < 0.01) correlated with mtDNA copy number. These are precursors of dopamine, serotonin, melatonin, and glutamate synthesis, which are all important neurotransmitters implicated in ASD etiology (Marotta et al., 2020). The notable enrichment of metabolites derived from the dopamine pathway, which is closely tied to the cellular oxidation state (Ben-Shachar, 2020), is consistent with altered redox homeostasis in ASD. Moreover, glycine serves as the precursor to one-carbon metabolism (cysteine, methionine, and glutathione pathways). These are essential regulators of an oxidative state which have been implicated in the metabolomic profile associated with ASD severity in two independent cohorts (Smith et al., 2020; Needham et al., 2021) and have been identified as a link between DNA methylation and mitochondrial dysfunction (Maddocks et al., 2016; Shen et al., 2020). In addition, both PGC-1α methylation and mtDNA copy number correlate significantly with metabolites derived from BCAA catabolism. BCAAs provide nitrogen to the glutamate-glutamine cycle and have been implicated as important regulators of glutamatergic neurotransmission (Yudkoff, 2017), which contribute to ASD etiology (Essa et al., 2013). Two metabolites [3-H-3-MGA and 3-methylglutaconic acid (3-MGA)] were significantly elevated in our ASD cohort (Stathopoulos et al., 2020) and are characterized as urinary biomarkers of mitochondrial respiratory chain deficiencies (Reinecke et al., 2012; Smuts et al., 2013). Altered BCAA metabolism has also been linked to oxidative stress resulting from perturbed NAD+/NADH redox ratios (Esterhuizen et al., 2017). Of note, mtDNA copy number correlated significantly with the direct precursor of de novo NAD+ synthesis, QA (p = 0.008), which is a neurotoxin (Notarangelo and Pocivavsek, 2017) that is also implicated in ASD etiology (Savino et al., 2020). NAD+ has also been identified as a metabolic link between changes in mtDNA copy number, methionine metabolism and DNA methylation (Lozoya et al., 2018). Therefore, both PGC-1α methylation and mtDNA copy number are associated with metabolomic evidence of mitochondrial dysfunction, oxidative stress, and neuropathology. Together, our metabolic data are consistent with a dysregulation of the link between methionine metabolism, mitochondrial dysfunction, and neurotransmitter synthesis in ASD.
Our study is one of the few molecular studies from Sub-Saharan Africa to examine ASD in an understudied African population. We present the first report of an association between the methylation of genes involved in mitochondrial biogenesis and remodeling, and mitochondrial function in a South African ASD cohort. While our results are correlative and cannot establish causality, they contribute to and are supported by, a growing body of evidence that point to aberrant DNA methylation and mitochondrial dysfunction in ASD etiology. Our results highlight the value of epigenetic research in under-studied populations to highlight novel associations. The central role of DNA methylation in modulating mitochondrial function highlights the potential to explore existing mitochondrial medications as putative therapeutic interventions for ASD symptomatology.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Ethics statement
The studies involving human participants were reviewed and approved by University of Cape Town, FSREC076-2014. Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin.
Author contributions
CO’R conceptualized the overall study design, was responsible for the phenotype data, supervised the laboratory work and data analysis, and was the major contributor in writing the manuscript. SB assisted with the design, data acquisition, and analysis of the mitochondrial DNA copy number data, and contributed to writing the manuscript. EB assisted with the design and analysis of the methylation data for the mitochondrial fission and fusion genes and contributed to writing the manuscript. CM assisted with the design and analysis of the methylation data for the mitochondrial biogenesis genes, analyzed the urinary metabolic data, and contributed to writing the manuscript. All authors read and approved the final manuscript.
Funding
This research was supported by the National Research Foundation, South Africa (Grant number 118524). The content hereof is the sole responsibility of the authors and does not necessarily represent the official views of the funder.
Acknowledgments
We thank the National Research Foundation, South Africa for funding this research. We are grateful to the families and the staff of all the schools who participated in our research.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2021.696428/full#supplementary-material
Supplementary Figure 1PGC-1α is differentially methylated in ASD. Targeted Next-Generation Bisulfite Sequencing (n = 55 ASD, n = 44 controls) identified differential methylation 12 CpG sites on PGC-1α. Differential methylation was identified using a two-tailed unpaired t-test with unequal variance (p < 0.05). Data represent the percent methylation in each individual at each site; box plots represent the mean percent methylation across either ASD or controls at each site; error bars represent standard deviations.
Supplementary Table 1Summary of demographic and phenotypic characteristics of participants used in this study. ASD, autism spectrum disorder.
Supplementary Table 2Differentially methylated CpG sites (p < 0.05, methylation range > 5%) of PGC-1α and STOML2 (n = 55 ASD, 44 controls) and GABPA, MFN2, FIS1, and OPA1 (n = 22 ASD, n = 22 controls).
Supplementary Table 3ASpearman correlation coefficients for urinary organic acids that correlate significantly with mtDNA copy number (n = 20 ASD, n = 13 controls). The most significant (p < 0.01) correlations with r < −0.3 or >0.3 are shown in bold, non-significant correlations are not shown. BCAA, branched chain amino acids; TCA, tricarboxylic acid cycle.
Supplementary Table 3BSpearman correlation coefficients for urinary organic acids (n = 20 ASD, n = 13 controls) that correlate significantly with methylation at promoter of PGC-1α at the CpG#1 site.
Abbreviations
- ADOS-2
Autism Diagnostic Observation Schedule Second Edition
- ASD
autism spectrum disorder
- B2M
Beta-2-microglobulin
- BCAAs
Branched-chain amino acids
- CREB1
CAMP response binding element 1
- DM
Differentially methylated
- FIS1
Mitochondrial fission protein 1
- GABPA
GA binding protein transcription factor subunit alpha
- mtDNA
Mitochondrial DNA
- MFN2
Mitofusin 2
- Nrf-2A
Nuclear respiratory factor 2 alpha subunit
- ND1
Mitochondrially encoded NADH ubiquinone oxidoreductase core subunit 1
- ND4
Mitochondrially encoded NADH ubiquinone oxidoreductase core subunit 4
- OPA1
Optic atrophy 1
- PGC-1 α
Peroxisome proliferator-activated receptor-gamma coactivator-1 alpha
- STOML2
Stomatin-like protein 2
- TSS
Transcriptional start site.
References
1
Al-KafajiG.GolbaharJ. (2013). High glucose-induced oxidative stress increases the copy number of mitochondrial DNA in human mesangial cells.Biomed. Res. Int.2013:754946. 10.1155/2013/754946
2
American Psychiatric Association [APA] (2013). Diagnostic and Statistical Manual of Mental Disorders, 5th Edn. Arlington, VA: American Psychiatric Association.
3
AndrewsS. V.SheppardB.WindhamG. C.SchieveL. A.SchendelD. E.CroenL. A.et al (2018). Case-control meta-analysis of blood DNA methylation and autism spectrum disorder.Mol. Autism9:40. 10.1186/s13229-018-0224-6
4
ArechederraM.DaianF.YimA.BazaiS. K.RichelmeS.DonoR.et al (2018). Hypermethylation of gene body CpG islands predicts high dosage of functional oncogenes in liver cancer.Nat. Commun.9:3164. 10.1038/s41467-018-05550-5
5
AyhanF.KonopkaG. (2019). Regulatory genes and pathways disrupted in autism spectrum disorders.Prog. Neuropsychopharmacol. Biol. Psychiatry8957–64. 10.1016/j.pnpbp.2018.08.017
6
BaiR. K.PerngC. L.HsuC. H.WongL. J. (2004). Quantitative PCR analysis of mitochondrial DNA content in patients with mitochondrial disease.Ann. N.Y. Acad. Sci.1011304–309. 10.1007/978-3-662-41088-2_29
7
BaroneR.RizzoR.TabbíG.MalaguarneraM.FryeR. E.BastinJ. (2019). Nuclear Peroxisome Proliferator-Activated Receptors (PPARs) as therapeutic targets of resveratrol for autism spectrum disorder.Int. J. Mol. Sci.20:1878. 10.3390/ijms20081878
8
BarrèsR.OslerM. E.YanJ.RuneA.FritzT.CaidahlK.et al (2009). Non-CpG methylation of the PGC-1α promoter through DNMT3B controls mitochondrial density.Cell Metab.10189–198. 10.1016/j.cmet.2009.07.011
9
Ben-ShacharD. (2020). The bimodal mechanism of interaction between dopamine and mitochondria as reflected in Parkinson’s disease and in schizophrenia.J. Neural Transm.127159–168. 10.1007/s00702-019-02120-x
10
BerkoE. R.SuzukiM.BerenF.LemetreC.AlaimoC. M.CalderR. B.et al (2014). Mosaic epigenetic dysregulation of ectodermal cells in autism spectrum disorder.PLoS Genet.10:e1004402. 10.1371/journal.pgen.1004402
11
BjørklundG.MeguidN. A.El-BanaM. A.TinkovA. A.SaadK.DadarM.et al (2020). oxidative stress in autism spectrum disorder.Mol. Neurobiol.572314–2332. 10.1007/s12035-019-01742-2
12
BuX.WuD.LuX.YangL.XuX.WangJ.et al (2017). Role of SIRT1/PGC-1α in mitochondrial oxidative stress in autistic spectrum disorder.Neuropsychiatr. Dis. Treat.131633–1645. 10.2147/NDT.S129081
13
CarrascoM.SalazarC.TiznadoW.RuizL. M. (2019). Alterations of mitochondrial biology in the oral mucosa of chilean children with Autism Spectrum Disorder (ASD).Cells8:367. 10.3390/cells8040367
14
CastellaniC. A.LongchampsR. J.SunJ.GuallarE.ArkingD. E. (2020). Thinking outside the nucleus: mitochondrial DNA copy number in health and disease.Mitochondrion53214–223. 10.1016/j.mito.2020.06.004
15
CaswellT. A.DroettboomM.LeeA.Sales de AndradeE.HunterJ.HoffmannT.et al (2021). Matplotlib/Matplotlib: REL: v3.4.2 (Version v3.4.2).Geneva: Zenodo.
16
ChenC.VincentA. E.BlainA. P.SmithA. L.TurnbullD. M.ReeveA. K. (2020). Investigation of mitochondrial biogenesis defects in single substantia nigra neurons using post-mortem human tissues.Neurobiol. Dis.134:104631. 10.1016/j.nbd.2019.104631
17
ChenH.DetmerS. A.EwaldA. J.GriffinE. E.FraserS. E.ChanD. C. (2003). Mitofusins Mfn1 and Mfn2 coordinately regulate mitochondrial fusion and are essential for embryonic development.J. Cell Biol.160189–200. 10.1083/jcb.200211046
18
ChenM.MacphersonA.OwensJ.WittertG.HeilbronnL. K. (2012). Obesity alone or with type 2 diabetes is associated with tissue specific alterations in DNA methylation and gene expression of PPARGC1A and IGF2.J. Diab. Res. Clin. Metab.1:16. 10.7243/2050-0866-1-16
19
ChenS.LiZ.HeY.ZhangF.LiH.LiaoY.et al (2015). Elevated mitochondrial DNA copy number in peripheral blood cells is associated with childhood autism.BMC Psychiatry15:50. 10.1186/s12888-015-0432-y
20
ChungJ. K.LeeS. Y.ParkM.JooE. J.KimS. A. (2019). Investigation of mitochondrial DNA copy number in patients with major depressive disorder.Psychiatry Res.282:112616. 10.1016/j.psychres.2019.112616
21
CitrignoL.MugliaM.QualtieriA.SpadaforaP.CavalcantiF.PioggiaG.et al (2020). The mitochondrial dysfunction hypothesis in autism spectrum disorders: current status and future perspectives.Int. J. Mol. Sci.2:5785. 10.3390/ijms21165785
22
D’AquilaP.De RangoF.GuarasciF.MandalàM.CorsonelloA.BellizziD.et al (2020). Erratum: multi-tissue DNA methylation remodeling at mitochondrial quality control genes according to diet in rat ageing models.Nutrients12:1494. 10.3390/nu12051494
23
DhillonS.HellingsJ. A.ButlerM. G. (2011). Genetics and mitochondrial abnormalities in autism spectrum disorders: a review.Curr. Genomics12322–332. 10.2174/138920211796429745
24
Duvezin-CaubetS.JagasiaR.WagenerJ.HofmannS.TrifunovicA.HanssonA.et al (2006). Proteolytic processing of OPA1 links mitochondrial dysfunction to alterations in mitochondrial morphology.J. Biol. Chem.28137972–37979. 10.1074/jbc.M606059200
25
ErlankH.ElmannA.KohenR.KannerJ. (2011). Polyphenols activate Nrf2 in astrocytes via H2O2, semiquinones, and quinones.Free Radic. Biol. Med.512319–2327. 10.1016/j.freeradbiomed.2011.09.033
26
EssaM. M.BraidyN.VijayanK. R.SubashS.GuilleminG. J. (2013). Excitotoxicity in the pathogenesis of autism.Neurotox. Res.23393–400. 10.1007/s12640-012-9354-3
27
EsterhuizenK.van der WesthuizenF. H.LouwR. (2017). Metabolomics of mitochondrial disease.Mitochondrion3597–110. 10.1016/j.mito.2017.05.012
28
FelicianoP.ZhouX.AstrovskayaI.TurnerT. N.WangT.BrueggemanL.et al (2019). Exome sequencing of 457 autism families recruited online provides evidence for autism risk genes.NPJ Genom. Med.4:19. 10.1038/s41525-019-0093-8
29
FryeR. E.RossignolD. A. (2011). Mitochondrial dysfunction can connect the diverse medical symptoms associated with autism spectrum disorders.Pediatr. Res.69(5 Pt. 2), 41–47. 10.1203/PDR.0b013e318212f16b
30
GauglerT.KleiL.SandersS. J.BodeaC. A.GoldbergA. P.LeeA. B.et al (2014). Most genetic risk for autism resides with common variation.Nat. Genet.46881–885. 10.1038/ng.3039.Most
31
GazievA. I.AbdullaevS.PodlutskyA. (2014). Mitochondrial function and mitochondrial DNA maintenance with advancing age.Biogerontology15417–438. 10.1007/s10522-014-9515-2
32
GillbergL.JacobsenS. C.RönnT.BrønsC.VaagA. (2014). PPARGC1A DNA methylation in subcutaneous adipose tissue in low birth weight subjects - impact of 5 days of high-fat overfeeding.Metabolism63263–271. 10.1016/j.metabol.2013.10.003
33
GiuliviC.ZhangY.Omanska-KlusekA.Ross-IntaC.WongS.Hertz-PicciottoI.et al (2010). Mitochondrial Dysfunction in Autism.JAMA3042389–2396. 10.1001/jama.2010.1706
34
GlintonK. E.ElseaS. H. (2019). Untargeted metabolomics for autism spectrum disorders: current status and future directions.Front. Psychiatry10:647. 10.3389/fpsyt.2019.00647
35
GradyJ. P.MurphyJ. L.BlakelyE. L.HallerR. G.TaylorR. W.TurnbullD. M.et al (2014). Accurate measurement of mitochondrial DNA deletion level and copy number differences in human skeletal muscle.PLoS One9:e114462. 10.1371/journal.pone.0114462
36
GuF.ChauhanV.KaurK.BrownW. T.LaFauciG.WegielJ.et al (2013). Alterations in mitochondrial DNA copy number and the activities of electron transport chain complexes and pyruvate dehydrogenase in the frontal cortex from subjects with autism.Transl. Psychiatry3:e299. 10.1038/tp.2013.68
37
HeinonenS.BuzkovaJ.MuniandyM.KaksonenR.OllikainenM.IsmailK.et al (2015). Impaired mitochondrial biogenesis in adipose tissue in acquired obesity.Diabetes643135–3145. 10.2337/db14-1937
38
HunterJ. D. (2007). Matplotlib: a 2D graphics environment.Comput. Sci. Eng.990–95. 10.1109/MCSE.2007.55
39
IshiharaN.FujitaY.OkaT.MiharaK. (2006). Regulation of mitochondrial morphology through proteolytic cleavage of OPA1.EMBO J.252966–2977. 10.1038/sj.emboj.7601184
40
KarbowskiM.YouleR. J. (2003). Dynamics of mitochondrial morphology in healthy cells and during apoptosis.Cell Death Differ.10870–880. 10.1038/sj.cdd.4401260
41
KelstrupL.HjortL.Houshmand-OeregaardA.TineC. D.HansenN. S.BroholmC.et al (2016). Gene expression and DNA methylation of PPARGC1A in muscle and adipose tissue from adult offspring of women with diabetes in pregnancy.Diabetes652900–2910. 10.2337/db16-0227
42
KimJ. I.LeeS. Y.ParkM.KimS. Y.KimJ. W.KimS. A.et al (2019). Peripheral mitochondrial DNA copy number is increased in Korean attention-deficit hyperactivity disorder patients.Front. Psychiatry.10:506. 10.3389/fpsyt.2019.00506
43
KolbasinaN. A.GureevA. P.SerzhantovaO. V.MikhailovA. A.MoshurovI. P.StarkovA. A.et al (2020). Lung cancer increases H2O2 concentration in the exhaled breath condensate, extent of mtDNA damage, and mtDNA copy number in buccal mucosa.Heliyon6:e04303. 10.1016/j.heliyon.2020.e04303
44
KresovichJ. K.JoyceB. T.GaoT.ZhengY.ZhangZ.AchenbachC. J.et al (2018). Promoter methylation of PGC1A and PGC1B predicts cancer incidence in a veteran cohort.Epigenomics10733–743. 10.2217/epi-2017-0141
45
KumarS.ReynoldsK.JiY.GuR.RaiS.ZhouC. J. (2019). Impaired neurodevelopmental pathways in autism spectrum disorder: a review of signaling mechanisms and crosstalk.J. Neurodev. Disord.11:10. 10.1186/s11689-019-9268-y
46
Ladd-AcostaC.HansenK. D.BriemE.FallinM. D.KaufmannW. E.FienbergA. P. (2014). Common DNA methylation alterations in multiple brain regions in autism.Mol. Psychiatry19862–871. 10.1038/mp.2013.114
47
LeeH.YinP.LuC.ChiC.WeiY. (2000). Increase of mitochondria and mitochondrial DNA in response to oxidative stress in human cells.Biochem. J.348425–432. 10.1042/bj3480425
48
LillycropK. A.PhillipsE. S.JacksonA. A.HansonM. A.BurdgeG. C. (2005). Dietary protein restriction of pregnant rats induces and folic acid supplementation prevents epigenetic modification of hepatic gene expression in the offspring.J. Nutr.1351382–1386. 10.1093/jn/135.6.1382
49
LingC.Del GuerraS.LupiR.RönnT.GranhallC.LuthmanH.et al (2008). Epigenetic regulation of PPARGC1A in human type 2 diabetic islets and effect on insulin secretion.Diabetologia51615–622. 10.1007/s00125-007-0916-5
50
LiuC.MaJ.ZhangJ.ZhaoH.ZhuY.QiJ.et al (2019). Testosterone deficiency caused by castration modulates mitochondrial biogenesis through the AR/PGC1α/TFAM pathway.Front. Genet.10:505. 10.3389/fgene.2019.00505
51
LordC.BrughaT. S.CharmanT.CusackJ.DumasG.FrazierT.et al (2020). Autism spectrum disorder.Nat. Rev. Dis. Primers6:5. 10.1038/s41572-019-0138-4
52
LoweR.GemmaC.BeyanH.HawaM. I.BazeosA.David LeslieR.et al (2013). Buccals are likely to be a more informative surrogate tissue than blood for epigenome-wide association studies.Epigenetics8445–454. 10.4161/epi.24362
53
LozoyaO. A.Martinez-ReyesI.WangT.GrenetD.BushelP.LiJ.et al (2018). Mitochondrial nicotinamide adenine dinucleotide reduced (n.d.) oxidation links the tricarboxylic acid (TCA) cycle with methionine metabolism and nuclear DNA methylation.PLoS Biol.16:e2005707. 10.1371/journal.pbio.2005707
54
MaddocksO. D. K.LabuschagneC. F.AdamsP. D.VousdenK. H. (2016). Serine metabolism supports the methionine cycle and DNA/RNA methylation through de novo ATP synthesis in cancer cells.Mol. Cell61210–221. 10.1016/j.molcel.2015.12.014
55
MalikA. N.CzajkaA. (2013). Is mitochondrial DNA content a potential biomarker of mitochondrial dysfunction?Mitochondrion13481–492. 10.1016/j.mito.2012.10.011
56
ManziB.LoizzoA. L.GianaG.CuratoloP. (2008). Autism and metabolic diseases.J. Child Neurol.23307–314. 10.1177/0883073807308698
57
MarottaR.RisoleoM. C.MessinaG.ParisiL.CarotenutoM.VetriL.et al (2020). The neurochemistry of Autism.Brain Sci.10:163. 10.3390/brainsci10030163
58
MasserD. R.ClarkN. W.Van RemmenH.FreemanW. M. (2016). Loss of the antioxidant enzyme CuZnSOD (Sod1) mimics an age-related increase in absolute mitochondrial DNA copy number in the skeletal muscle.Age Dordr.38323–333. 10.1007/s11357-016-9930-1
59
MelnykS.FuchsG. J.SchulzE.LopezM.KahlerS. G.FussellJ. J.et al (2012). Metabolic imbalance associated with methylation dysregulation and oxidative damage in children with autism.J. Autism Dev. Disord.42367–377. 10.1007/s10803-011-1260-7
60
MichelS.WanetA.De PauwA.RommelaereG.ArnouldT.RenardP. (2012). Crosstalk between mitochondrial (dys)function and mitochondrial abundance.J. Cell Physiol.2272297–2310. 10.1002/jcp.23021
61
MontierL. L. C.DengJ.BaiY. (2009). Numbers matters: control of mammalian mitochondrial DNA copy number.J. Genet. Genomics36125–131. 10.1016/S1673-8527(08)60099-5
62
NardoneS.SamsD. S.ReuveniE.GetselterD.OronO.KarpujM.et al (2014). DNA methylation analysis of the autistic brain reveals multiple dysregulated biological pathways.Transl. Psychiatry4:e433. 10.1038/tp.2014.70
63
NardoneS.SamsD. S.ZitoA.ReuveniE.ElliottE. (2017). Dysregulation of cortical neuron DNA methylation profile in autism spectrum disorder.Cereb. Cortex275739–5754. 10.1093/cercor/bhx250
64
NeedhamB. D.AdameM. D.SerenaG.RoseD. R.PrestonG. M.ConradM. C.et al (2021). Plasma and fecal metabolite profiles in autism spectrum disorder.Biol. Psychiatry89451–462. 10.1016/j.biopsych.2020.09.025
65
NotarangeloF. M.PocivavsekA. (2017). Elevated kynurenine pathway metabolism during neurodevelopment: implications for brain and behavior.Neuropharmacology112(Pt. B), 275–285. 10.1016/j.neuropharm.2016.03.00
66
OliveiraG.DiogoL.GrazinaM.GarciaP.AtaídeA.MarquesC.et al (2005). Mitochondrial dysfunction in autism spectrum disorders: a population-based study.Dev. Med. Child. Neurol.47185–189. 10.1017/S0012162205000332
67
PecorelliA.FerraraF.MessanoN.CordoneV.SchiavoneM. L.CervellatiF.et al (2020). Alterations of mitochondrial bioenergetics, dynamics, and morphology support the theory of oxidative damage involvement in autism spectrum disorder.FASEB J.346521–6538. 10.1096/fj.201902677r
68
PetersenM. H.Budtz-JørgensenE.SørensenS. A.NielsenJ. E.HjermindL. E.Vinther-JensenT.et al (2014). Reduction in mitochondrial DNA copy number in peripheral leukocytes after onset of Huntington’s disease.Mitochondrion1714–21. 10.1016/j.mito.2014.05.001
69
PetrieM. A.SharmaA.TaylorE. B.SunejaM.ShieldsR. K. (2020). Impact of short-and long-term electrically induced muscle exercise on gene signaling pathways, gene expression, and PGC1a methylation in men with spinal cord injury.Physiol. Genomics5271–80. 10.1152/physiolgenomics.00064.2019
70
PloumiC.DaskalakiI.TavernarakisN. (2017). Mitochondrial biogenesis and clearance: a balancing act.FEBS J.284183–195. 10.1111/febs.13820
71
RauluseviciuteI.DrabløsF.RyeM. B. (2020). DNA hypermethylation associated with upregulated gene expression in prostate cancer demonstrates the diversity of epigenetic regulation.BMC Med. Genomics13:6. 10.1186/s12920-020-0657-6
72
RaybautP. (2009). Spyder-Documentation. Available online at: pythonhosted.org
73
ReineckeC. J.KoekemoerG.van der WesthuizenF. H.LouwR.LindequeJ. Z.MienieL. J.et al (2012). Metabolomics of urinary organic acids in respiratory chain deficiencies in children.Metabolomics8264–283. 10.1007/s11306-011-0309-0
74
RidoutK. K.CoeJ. L.ParadeS. H.MarsitC. J.KaoH. T.PortonB.et al (2020). Molecular markers of neuroendocrine function and mitochondrial biogenesis associated with early life stress.Psychoneuroendocrinology116:104632. 10.1016/j.psyneuen.2020.104632
75
RylaarsdamL.Guemez-GamboaA. (2019). Genetic causes and modifiers of autism spectrum disorder.Front. Cell Neurosci.13:385. 10.3389/fncel.2019.00385
76
SandinS.LichtensteinP.Kuja-HalkolaR.HultmanC.LarssonH.ReichenbergA. (2017). The heritability of autism spectrum disorder.JAMA3181182–1184. 10.1001/jama.2017.12141
77
SantosR. X.CorreiaS. C.AlvesM. G.OliveiraP. F.CardosoS.CarvalhoC.et al (2014). Insulin therapy modulates mitochondrial dynamics and biogenesis, autophagy and tau protein phosphorylation in the brain of type 1 diabetic rats.Biochim. Biophys. Acta18421154–1166. 10.1016/j.bbadis.2014.04.011
78
SavinoR.CarotenutoM.PolitoA. N.Di NoiaS.AlbenzioM.ScarinciA.et al (2020). Analyzing the potential biological determinants of autism spectrum disorder: from neuroinflammation to the kynurenine pathway.Brain Sci.10:631. 10.3390/brainsci10090631
79
ShenW.GaoC.CuetoR.LiuL.FuH.ShaoY.et al (2020). Homocysteine-methionine cycle is a metabolic sensor system controlling methylation-regulated pathological signaling.Redox Biol.28:101322. 10.1016/j.redox.2019.101322
80
SinghK.SinghI. N.DigginsE.ConnorsS. L.KarimM. A.LeeD.et al (2020). Developmental regression and mitochondrial function in children with autism.Ann. Clin. Transl. Neurol.7683–694. 10.1002/acn3.51034
81
SmirnovaE.GriparicL.ShurlandD. L.Van Der BliekA. M. (2001). Dynamin-related protein Drp1 is required for mitochondrial division in mammalian cells.Mol. Biol. Cell122245–2256. 10.1091/mbc.12.8.2245
82
SmithA. K.KilaruV.KlengelT.MercerK. B.BradleyB.ConneelyK. N.et al (2015). DNA extracted from saliva for methylation studies of psychiatric traits: evidence tissue specificity and relatedness to brain.Am. J. Med. Genet. B Neuropsychiatr. Genet.168B36–44. 10.1002/ajmg.b.32278
83
SmithA. M.NatowiczM. R.BraasD.LudwigM. A.NeyD. M.DonleyE. L. R.et al (2020). A metabolomics approach to screening for autism risk in the children’s autism metabolome project.Autism Res.131270–1285. 10.1002/aur.2330
84
SmutsI.van der WesthuizenF. H.LouwR.MienieL. J.EngelkeU. F. H.WeversR. A.et al (2013). Disclosure of a putative biosignature for respiratory chain disorders through a metabolomics approach.Metabolomics9379–391. 10.1007/s11306-012-0455-z
85
SookoianS.RosselliM. S.GemmaC.BurgueñoA. L.GianottiT. F.CastañoG. O.et al (2010). Epigenetic regulation of insulin resistance in nonalcoholic fatty liver disease: impact of liver methylation of the peroxisome proliferator-activated receptor γ coactivator 1α promoter.Hepatology521992–2000. 10.1002/hep.23927
86
StathopoulosS.GaujouxR.LindequeZ.MahonyC.Van Der ColffR.Van Der WesthuizenF.et al (2020). DNA methylation associated with mitochondrial dysfunction in a South African autism spectrum disorder cohort.Autism Res.131079–1093. 10.1002/aur.2310
87
StojanovskiD.KoutsopoulosO. S.OkamotoK.RyanM. T. (2004). Levels of human Fis1 at the mitochondrial outer membrane regulate mitochondrial morphology.J. Cell Sci.1171201–1210. 10.1242/jcs.01058
88
St-PierreJ.DroriS.UldryM.SilvaggiJ. M.RheeJ.JägerS.et al (2006). Suppression of reactive oxygen species and neurodegeneration by the PGC-1 transcriptional coactivators.Cell127397–408. 10.1016/j.cell.2006.09.024
89
SunJ.LongchampsR. J.PiggottD. A.CastellaniC. A.SumpterJ. A.BrownT. T.et al (2019). Association between HIV infection and mitochondrial DNA copy number in peripheral blood: a population-based, prospective cohort study.J. Infect. Dis.2191285–1293. 10.1093/infdis/jiy658
90
TangG.RiosP. G.KuoS. H.AkmanH. O.RosoklijaG.TanjiK.et al (2013). Mitochondrial abnormalities in temporal lobe of autistic brain.Neurobiol. Dis.54349–361. 10.1016/j.nbd.2013.01.006
91
ThompsonK.MajdH.DallabonaC.ReinsonK.KingM. S.AlstonC. L.et al (2016). Recurrent de novo dominant mutations in SLC25A4 cause severe early-onset mitochondrial disease and loss of mitochondrial DNA copy number.Am. J. Hum. Genet.99860–876. 10.1016/j.ajhg.2016.08.014
92
TonderaD.GrandemangeS.JourdainA.KarbowskiM.MattenbergerY.HerzigS.et al (2009). SLP-2 is required for stress-induced mitochondrial hyperfusion.EMBO J.281589–1600. 10.1038/emboj.2009.89
93
TremblayM. W.JiangY. H. (2019). DNA methylation and susceptibility to autism spectrum disorder.Annu. Rev. Med.70151–166. 10.1146/annurev-med-120417-091431
94
ValentiD.de BariL.De FilippisB.Henrion-CaudeA.VaccaR. A. (2014). Mitochondrial dysfunction as a central actor in intellectual disability-related diseases: an overview of down syndrome, autism, fragile X and rett syndrome.Neurosci. Biobehav. Rev.46202–217. 10.1111/j.1469-8749.2005.tb01113.x
95
Valiente-PallejàA.TorrellH.MuntanéG.CortésM. J.Martínez-LealR.AbasoloN.et al (2018). Genetic and clinical evidence of mitochondrial dysfunction in autism spectrum disorder and intellectual disability.Hum. Mol. Genet.27891–900. 10.1093/hmg/ddy009
96
Ventura-ClapierR.GarnierA.VekslerV. (2008). Transcriptional control of mitochondrial biogenesis: the central role of PGC-1α.Cardiovasc. Res.79208–217. 10.1093/cvr/cvn098
97
VoineaguI. (2012). Gene expression studies in autism: moving from the genome to the transcriptome and beyond.Neurobiol. Dis.4569–75. 10.1016/j.nbd.2011.07.017
98
WaiT.LangerT. (2016). Mitochondrial dynamics and metabolic regulation.Trends Endocrinol. Metab.27105–117. 10.1016/j.tem.2015.12.001
99
Wiśniowiecka-KowalnikB.NowakowskaB. A. (2019). Genetics and epigenetics of autism spectrum disorder - current evidence in the field.J. Appl. Genet.6037–47. 10.1007/s13353-018-00480-w
100
WongC. C. Y.MeaburnE. L.RonaldA.PriceT. S.JeffriesA. R.SchalkwykL. C.et al (2014). Methylomic analysis of monozygotic twins discordant for autism spectrum disorder and related behavioural traits.Mol. Psychiatry19495–503. 10.1038/mp.2013.41
101
WongC. C. Y.SmithR. G.HannonE.RamaswamiG.ParikshakN. N.AssaryE.et al (2019). Genome-wide DNA methylation profiling identifies convergent molecular signatures associated with idiopathic and syndromic autism in post-mortem human brain tissue.Hum. Mol. Genet.282201–2211. 10.1093/hmg/ddz052
102
WuZ.HuangX.FengY.HandschinC.FengY.GullicksenP. S.et al (2006). Transducer of regulated CREB-binding proteins (TORCs) induce PGC-1α transcription and mitochondrial biogenesis in muscle cells.Proc. Natl. Acad. Sci. U.S.A.10314379–14384. 10.1073/pnas.0606714103
103
YangX.XuS.QianY.HeX.ChenS.XiaoQ. (2020). Hypermethylation of the gene coding for PGC-1α in peripheral blood leukocytes of patients with Parkinson’s disease.Front. Neurosci.14:97. 10.3389/fnins.2020.00097
104
YooH. J.ParkM.KimS. A. (2017). Difference in mitochondrial DNA copy number in peripheral blood cells between probands with autism spectrum disorders and their unaffected siblings.World J. Biol. Psychiatry18151–156. 10.1080/15622975.2016.1234069
105
YoonS. H.ChoiJ.LeeW. J.DoJ. T. (2020). Genetic and epigenetic etiology underlying autism spectrum disorder.J. Clin. Med.9:966. 10.3390/jcm9040966
106
YuX.QiaoS.WangD.DaiJ.WangJ.ZhangR.et al (2016). A metabolomics-based approach for ranking the depressive level in a chronic unpredictable mild stress rat model.RSC Adv.625751–25765. 10.1039/c6ra00665e
107
YudkoffM. (2017). Interactions in the metabolism of glutamate and the branched-chain amino acids and ketoacids in the CNS.Neurochem. Res.4210–18. 10.1007/s11064-016-2057-z.INTERACTIONS
108
ZengY.PanQ.WangX.LiD.LinY.ManF.et al (2019). Impaired Mitochondrial fusion and oxidative phosphorylation triggered by high glucose is mediated by Tom22 in endothelial cells.Oxid. Med. Cell Longev.2019:4508762. 10.1155/2019/4508762
109
ZhaoH.ZhaoY.RenY.LiM.LiT.LiR.et al (2017). Epigenetic regulation of an adverse metabolic phenotype in polycystic ovary syndrome: the impact of the leukocyte methylation of PPARGC1A promoter.Fertil. Steril.107467–474.e5. 10.1016/j.fertnstert.2016.10.039
110
ZillerM. J.GuH.MüllerF.DonagheyJ.TsaiL. T. Y.KohlbacherO.et al (2013). Charting a dynamic DNA methylation landscape of the human genome.Nature500477–481. 10.1038/nature12433
Summary
Keywords
Autism Spectrum Disorder, methylation, PGC-1α, mtDNA copy number, metabolomics, mitochondrial dysfunction
Citation
Bam S, Buchanan E, Mahony C and O’Ryan C (2021) DNA Methylation of PGC-1α Is Associated With Elevated mtDNA Copy Number and Altered Urinary Metabolites in Autism Spectrum Disorder. Front. Cell Dev. Biol. 9:696428. doi: 10.3389/fcell.2021.696428
Received
16 April 2021
Accepted
05 July 2021
Published
26 July 2021
Volume
9 - 2021
Edited by
Alaguraj Veluchamy, St. Jude Children’s Research Hospital, United States
Reviewed by
Adam R. Wende, University of Alabama at Birmingham, United States; Hashim Islam, University of British Columbia Okanagan, Canada; Rick Schnellmann, University of Arizona, United States
Updates
Copyright
© 2021 Bam, Buchanan, Mahony and O’Ryan.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Colleen O’Ryan, colleen.oryan@uct.ac.za
†These authors have contributed equally to this work and share first authorship
This article was submitted to Epigenomics and Epigenetics, a section of the journal Frontiers in Cell and Developmental Biology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.