Sophia Bam
Erin Buchanan
Caitlyn Mahony
Colleen O’Ryan- Department of Molecular and Cell Biology, University of Cape Town, Cape Town, South Africa
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 PGC-1α promoter is hypermethylated in autism spectrum disorder (ASD). (A) Diagrammatic representation of the PCG-1α gene promoter region (Chr4:23889974) showing the location of the eight hypermethylated CpG sites (black circles) relative to the transcription start site (TSS) and the sequence of the binding site for the transcription factor, CREB1. (B) Percentage methylation of the eight differentially methylated PGC-1α promoter CpG sites measured using Targeted Next-Generation Bisulfite Sequencing (ASD n = 55, controls n = 44) (UTR = untranslated region). 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 in either controls (black) or ASD (red) at each site; error bars represent standard deviations.
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. Key gene regulators of mitochondrial biogenesis, fission and fusion are differentially methylated in ASD. Box plots show the percentage methylation of the differentially methylated CpG sites (p < 0.05), measured using Targeted Next-Generation Bisulfite Sequencing of STOML2 (ASD n = 55, controls 44) and FIS, MFN2, OPA1, and GABPA (ASD n = 22, controls n = 22). 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 controls (black) or ASD (red) at each site; error bars represent standard deviations.
Figure 3. Differentially methylated genes converge on pathways regulating mitochondrial homeostasis in response to metabolic and oxidative stress. Mitochondrial homeostasis is maintained by differentially methylated genes (in red) involved in mitochondrial biogenesis, fission and fusion that converge on the regulation of mitochondrial copy number in response to metabolic and oxidative stress. Metabolic stress, which decreases ATP production, activates cAMP and AMPK signaling leading to the transcription and activation of PGC-1α. PGC-1α upregulates the expression of Nrf1 and GABPA which induce the transcription of the mitochondrial transcription factors to facilitate mitochondrial biogenesis. PGC-1α also upregulates the transcription of antioxidant enzymes to modulate oxidative stress or ROS levels. Oxidative stress activates the redox-sensitive Nfe2l2 pathway, which upregulates both GABPA and the Pink-Parkin pathway that controls mitophagy, fission and fusion (Nfe2l2 = nuclear factor erythroid 2-related factor 2; NRF1 = Nuclear respiratory factor 1; p62 = ubiquitin-binding protein p62; Pink = PTEN-induced kinase 1; ROS = reactive oxygen species).
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. Relationship between mitochondrial DNA (mtDNA) copy number, mtDNA deletion and PGC-1α promoter methylation. (A) mtDNA copy number is significantly elevated in ASD. Relative quantification of mtDNA copy number was performed by multiplex real-time qPCR of MT-ND1 and B2M (ASD n = 59, controls n = 40). (B) A non-significant increase in mtDNA deletions is observed in ASD. Relative quantification of mtDNA deletion was performed by multiplex real-time qPCR of MT-ND1 and MT-ND4 (ASD n = 59, controls n = 40). Significance was established using Student t-tests where (p < 0.05). (C) mtDNA copy number correlates significantly with mtDNA deletion in ASD and controls (n = 99), where Spearman’s rho = 0.882, p < 0.001. (D) The relationship between PGC-1α promoter methylation, mtDNA copy number and mtDNA deletions is altered in ASD. Heatmap shows Spearman rank correlations between mtDNA deletions, mtDNA copy number and percentage methylation at differentially methylated CpG sites in the PGC-1α promoter (ASD n = 49, control n = 42). Spearman’s rho is represented according to the color key provided, ∗ indicates p < 0.05. PGC-1α CpG#1 correlated positively with mtDNA copy number in ASD (Spearman’s r = 0.9, p = 8.814 10-10) while CpG#4 correlated negatively with both mtDNA copy number and (Spearman’s r = –0.4, p = 0.045) and mtDNA deletions (Spearman’s r = –0.4, p = 0.032) in controls.
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. Urinary metabolite, 3-hydroxy-3-methylglutaric acid (3-H-3-MGA), correlates with PGC-1α methylation and mitochondrial DNA (mtDNA) copy number. Normal linear regression analysis shows that (A) DNA methylation at PGC-1α CpG#1 correlates with 3-H-3-MGA levels, RSQ = 0,2163, Spearman Rho = 0.386, p = 0,043 and (B) mtDNA copy number correlates with 3-H-3-MGA, RSQ = 0,2282, Spearman Rho = 0.355, p = 0,004; (ASD n = 20, control n = 13).
Figure 6. Mitochondrial DNA (mtDNA) copy number correlates with urinary metabolites implicating mitochondrial metabolism and neurotransmitter synthesis. Urinary metabolites that most significantly correlate with mtDNA copy number (p < 0.01) converge on pathways involved in the metabolism of fatty acids (ethylmalonic acid); Branched Chain Amino Acids (BCAAS) (hydroxy-isovaleric acid, 3-hydroxy-3-methylglutaric acid); and neurotransmitters including Dopamine (2-hydroxyphenylacetic acid, 4-hydroxyhippuric acid, vanillylmandelic acid, and homovanillic acid), Serotonin (Quinolinic acid), Glutamate (2-ketoglutaric acid, D-L-hydroxyglutaric acid; 2,5-furandicarboxylic acid, hydroxymethyl-2-furoic acid), and Glycine (N-2Furanylcarbonylglycine, oxalic acid). Metabolites that are associated with elevated mtDNA copy number are shown in red; pathways involved in neurotransmitter metabolism are highlighted in blue.
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.
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.
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.
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.
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.
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 1 | PGC-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 1 | Summary of demographic and phenotypic characteristics of participants used in this study. ASD, autism spectrum disorder.
Supplementary Table 2 | Differentially 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 3A | Spearman 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 3B | Spearman 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
Al-Kafaji, G., and Golbahar, J. (2013). High glucose-induced oxidative stress increases the copy number of mitochondrial DNA in human mesangial cells. Biomed. Res. Int. 2013:754946. doi: 10.1155/2013/754946
American Psychiatric Association [APA] (2013). Diagnostic and Statistical Manual of Mental Disorders, 5th Edn. Arlington, VA: American Psychiatric Association.
Andrews, S. V., Sheppard, B., Windham, G. C., Schieve, L. A., Schendel, D. E., Croen, L. A., et al. (2018). Case-control meta-analysis of blood DNA methylation and autism spectrum disorder. Mol. Autism 9:40. doi: 10.1186/s13229-018-0224-6
Arechederra, M., Daian, F., Yim, A., Bazai, S. K., Richelme, S., Dono, R., et al. (2018). Hypermethylation of gene body CpG islands predicts high dosage of functional oncogenes in liver cancer. Nat. Commun. 9:3164. doi: 10.1038/s41467-018-05550-5
Ayhan, F., and Konopka, G. (2019). Regulatory genes and pathways disrupted in autism spectrum disorders. Prog. Neuropsychopharmacol. Biol. Psychiatry 89, 57–64. doi: 10.1016/j.pnpbp.2018.08.017
Bai, R. K., Perng, C. L., Hsu, C. H., and Wong, L. J. (2004). Quantitative PCR analysis of mitochondrial DNA content in patients with mitochondrial disease. Ann. N.Y. Acad. Sci. 1011, 304–309. doi: 10.1007/978-3-662-41088-2_29
Barone, R., Rizzo, R., Tabbí, G., Malaguarnera, M., Frye, R. E., and Bastin, J. (2019). Nuclear Peroxisome Proliferator-Activated Receptors (PPARs) as therapeutic targets of resveratrol for autism spectrum disorder. Int. J. Mol. Sci. 20:1878. doi: 10.3390/ijms20081878
Barrès, R., Osler, M. E., Yan, J., Rune, A., Fritz, T., Caidahl, K., et al. (2009). Non-CpG methylation of the PGC-1α promoter through DNMT3B controls mitochondrial density. Cell Metab. 10, 189–198. doi: 10.1016/j.cmet.2009.07.011
Ben-Shachar, D. (2020). The bimodal mechanism of interaction between dopamine and mitochondria as reflected in Parkinson’s disease and in schizophrenia. J. Neural Transm. 127, 159–168. doi: 10.1007/s00702-019-02120-x
Berko, E. R., Suzuki, M., Beren, F., Lemetre, C., Alaimo, C. M., Calder, R. B., et al. (2014). Mosaic epigenetic dysregulation of ectodermal cells in autism spectrum disorder. PLoS Genet. 10:e1004402. doi: 10.1371/journal.pgen.1004402
Bjørklund, G., Meguid, N. A., El-Bana, M. A., Tinkov, A. A., Saad, K., Dadar, M., et al. (2020). oxidative stress in autism spectrum disorder. Mol. Neurobiol. 57, 2314–2332. doi: 10.1007/s12035-019-01742-2
Bu, X., Wu, D., Lu, X., Yang, L., Xu, X., Wang, J., et al. (2017). Role of SIRT1/PGC-1α in mitochondrial oxidative stress in autistic spectrum disorder. Neuropsychiatr. Dis. Treat. 13, 1633–1645. doi: 10.2147/NDT.S129081
Carrasco, M., Salazar, C., Tiznado, W., and Ruiz, L. M. (2019). Alterations of mitochondrial biology in the oral mucosa of chilean children with Autism Spectrum Disorder (ASD). Cells 8:367. doi: 10.3390/cells8040367
Castellani, C. A., Longchamps, R. J., Sun, J., Guallar, E., and Arking, D. E. (2020). Thinking outside the nucleus: mitochondrial DNA copy number in health and disease. Mitochondrion 53, 214–223. doi: 10.1016/j.mito.2020.06.004
Caswell, T. A., Droettboom, M., Lee, A., Sales de Andrade, E., Hunter, J., Hoffmann, T., et al. (2021). Matplotlib/Matplotlib: REL: v3.4.2 (Version v3.4.2). Geneva: Zenodo.
Chen, C., Vincent, A. E., Blain, A. P., Smith, A. L., Turnbull, D. M., and Reeve, A. K. (2020). Investigation of mitochondrial biogenesis defects in single substantia nigra neurons using post-mortem human tissues. Neurobiol. Dis. 134:104631. doi: 10.1016/j.nbd.2019.104631
Chen, H., Detmer, S. A., Ewald, A. J., Griffin, E. E., Fraser, S. E., and Chan, D. C. (2003). Mitofusins Mfn1 and Mfn2 coordinately regulate mitochondrial fusion and are essential for embryonic development. J. Cell Biol. 160, 189–200. doi: 10.1083/jcb.200211046
Chen, M., Macpherson, A., Owens, J., Wittert, G., and Heilbronn, L. 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. doi: 10.7243/2050-0866-1-16
Chen, S., Li, Z., He, Y., Zhang, F., Li, H., Liao, Y., et al. (2015). Elevated mitochondrial DNA copy number in peripheral blood cells is associated with childhood autism. BMC Psychiatry 15:50. doi: 10.1186/s12888-015-0432-y
Chung, J. K., Lee, S. Y., Park, M., Joo, E. J., and Kim, S. A. (2019). Investigation of mitochondrial DNA copy number in patients with major depressive disorder. Psychiatry Res. 282:112616. doi: 10.1016/j.psychres.2019.112616
Citrigno, L., Muglia, M., Qualtieri, A., Spadafora, P., Cavalcanti, F., Pioggia, G., et al. (2020). The mitochondrial dysfunction hypothesis in autism spectrum disorders: current status and future perspectives. Int. J. Mol. Sci. 2:5785. doi: 10.3390/ijms21165785
D’Aquila, P., De Rango, F., Guarasci, F., Mandalà, M., Corsonello, A., Bellizzi, D., et al. (2020). Erratum: multi-tissue DNA methylation remodeling at mitochondrial quality control genes according to diet in rat ageing models. Nutrients 12:1494. doi: 10.3390/nu12051494
Dhillon, S., Hellings, J. A., and Butler, M. G. (2011). Genetics and mitochondrial abnormalities in autism spectrum disorders: a review. Curr. Genomics 12, 322–332. doi: 10.2174/138920211796429745
Duvezin-Caubet, S., Jagasia, R., Wagener, J., Hofmann, S., Trifunovic, A., Hansson, A., et al. (2006). Proteolytic processing of OPA1 links mitochondrial dysfunction to alterations in mitochondrial morphology. J. Biol. Chem. 281, 37972–37979. doi: 10.1074/jbc.M606059200
Erlank, H., Elmann, A., Kohen, R., and Kanner, J. (2011). Polyphenols activate Nrf2 in astrocytes via H2O2, semiquinones, and quinones. Free Radic. Biol. Med. 51, 2319–2327. doi: 10.1016/j.freeradbiomed.2011.09.033
Essa, M. M., Braidy, N., Vijayan, K. R., Subash, S., and Guillemin, G. J. (2013). Excitotoxicity in the pathogenesis of autism. Neurotox. Res. 23, 393–400. doi: 10.1007/s12640-012-9354-3
Esterhuizen, K., van der Westhuizen, F. H., and Louw, R. (2017). Metabolomics of mitochondrial disease. Mitochondrion 35, 97–110. doi: 10.1016/j.mito.2017.05.012
Feliciano, P., Zhou, X., Astrovskaya, I., Turner, T. N., Wang, T., Brueggeman, L., et al. (2019). Exome sequencing of 457 autism families recruited online provides evidence for autism risk genes. NPJ Genom. Med. 4:19. doi: 10.1038/s41525-019-0093-8
Frye, R. E., and Rossignol, D. A. (2011). Mitochondrial dysfunction can connect the diverse medical symptoms associated with autism spectrum disorders. Pediatr. Res. 69(5 Pt. 2), 41–47. doi: 10.1203/PDR.0b013e318212f16b
Gaugler, T., Klei, L., Sanders, S. J., Bodea, C. A., Goldberg, A. P., Lee, A. B., et al. (2014). Most genetic risk for autism resides with common variation. Nat. Genet. 46, 881–885. doi: 10.1038/ng.3039.Most
Gaziev, A. I., Abdullaev, S., and Podlutsky, A. (2014). Mitochondrial function and mitochondrial DNA maintenance with advancing age. Biogerontology 15, 417–438. doi: 10.1007/s10522-014-9515-2
Gillberg, L., Jacobsen, S. C., Rönn, T., Brøns, C., and Vaag, A. (2014). PPARGC1A DNA methylation in subcutaneous adipose tissue in low birth weight subjects - impact of 5 days of high-fat overfeeding. Metabolism 63, 263–271. doi: 10.1016/j.metabol.2013.10.003
Giulivi, C., Zhang, Y., Omanska-Klusek, A., Ross-Inta, C., Wong, S., Hertz-Picciotto, I., et al. (2010). Mitochondrial Dysfunction in Autism. JAMA 304, 2389–2396. doi: 10.1001/jama.2010.1706
Glinton, K. E., and Elsea, S. H. (2019). Untargeted metabolomics for autism spectrum disorders: current status and future directions. Front. Psychiatry 10:647. doi: 10.3389/fpsyt.2019.00647
Grady, J. P., Murphy, J. L., Blakely, E. L., Haller, R. G., Taylor, R. W., Turnbull, D. M., et al. (2014). Accurate measurement of mitochondrial DNA deletion level and copy number differences in human skeletal muscle. PLoS One 9:e114462. doi: 10.1371/journal.pone.0114462
Gu, F., Chauhan, V., Kaur, K., Brown, W. T., LaFauci, G., Wegiel, J., 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. Psychiatry 3:e299. doi: 10.1038/tp.2013.68
Heinonen, S., Buzkova, J., Muniandy, M., Kaksonen, R., Ollikainen, M., Ismail, K., et al. (2015). Impaired mitochondrial biogenesis in adipose tissue in acquired obesity. Diabetes 64, 3135–3145. doi: 10.2337/db14-1937
Hunter, J. D. (2007). Matplotlib: a 2D graphics environment. Comput. Sci. Eng. 9, 90–95. doi: 10.1109/MCSE.2007.55
Ishihara, N., Fujita, Y., Oka, T., and Mihara, K. (2006). Regulation of mitochondrial morphology through proteolytic cleavage of OPA1. EMBO J. 25, 2966–2977. doi: 10.1038/sj.emboj.7601184
Karbowski, M., and Youle, R. J. (2003). Dynamics of mitochondrial morphology in healthy cells and during apoptosis. Cell Death Differ. 10, 870–880. doi: 10.1038/sj.cdd.4401260
Kelstrup, L., Hjort, L., Houshmand-Oeregaard, A., Tine, C. D., Hansen, N. S., Broholm, C., et al. (2016). Gene expression and DNA methylation of PPARGC1A in muscle and adipose tissue from adult offspring of women with diabetes in pregnancy. Diabetes 65, 2900–2910. doi: 10.2337/db16-0227
Kim, J. I., Lee, S. Y., Park, M., Kim, S. Y., Kim, J. W., Kim, S. A., et al. (2019). Peripheral mitochondrial DNA copy number is increased in Korean attention-deficit hyperactivity disorder patients. Front. Psychiatry. 10:506. doi: 10.3389/fpsyt.2019.00506
Kolbasina, N. A., Gureev, A. P., Serzhantova, O. V., Mikhailov, A. A., Moshurov, I. P., Starkov, A. 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. Heliyon 6:e04303. doi: 10.1016/j.heliyon.2020.e04303
Kresovich, J. K., Joyce, B. T., Gao, T., Zheng, Y., Zhang, Z., Achenbach, C. J., et al. (2018). Promoter methylation of PGC1A and PGC1B predicts cancer incidence in a veteran cohort. Epigenomics 10, 733–743. doi: 10.2217/epi-2017-0141
Kumar, S., Reynolds, K., Ji, Y., Gu, R., Rai, S., and Zhou, C. J. (2019). Impaired neurodevelopmental pathways in autism spectrum disorder: a review of signaling mechanisms and crosstalk. J. Neurodev. Disord. 11:10. doi: 10.1186/s11689-019-9268-y
Ladd-Acosta, C., Hansen, K. D., Briem, E., Fallin, M. D., Kaufmann, W. E., and Fienberg, A. P. (2014). Common DNA methylation alterations in multiple brain regions in autism. Mol. Psychiatry 19, 862–871. doi: 10.1038/mp.2013.114
Lee, H., Yin, P., Lu, C., Chi, C., and Wei, Y. (2000). Increase of mitochondria and mitochondrial DNA in response to oxidative stress in human cells. Biochem. J. 348, 425–432. doi: 10.1042/bj3480425
Lillycrop, K. A., Phillips, E. S., Jackson, A. A., Hanson, M. A., and Burdge, G. 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. 135, 1382–1386. doi: 10.1093/jn/135.6.1382
Ling, C., Del Guerra, S., Lupi, R., Rönn, T., Granhall, C., Luthman, H., et al. (2008). Epigenetic regulation of PPARGC1A in human type 2 diabetic islets and effect on insulin secretion. Diabetologia 51, 615–622. doi: 10.1007/s00125-007-0916-5
Liu, C., Ma, J., Zhang, J., Zhao, H., Zhu, Y., Qi, J., et al. (2019). Testosterone deficiency caused by castration modulates mitochondrial biogenesis through the AR/PGC1α/TFAM pathway. Front. Genet. 10:505. doi: 10.3389/fgene.2019.00505
Lord, C., Brugha, T. S., Charman, T., Cusack, J., Dumas, G., Frazier, T., et al. (2020). Autism spectrum disorder. Nat. Rev. Dis. Primers 6:5. doi: 10.1038/s41572-019-0138-4
Lowe, R., Gemma, C., Beyan, H., Hawa, M. I., Bazeos, A., David Leslie, R., et al. (2013). Buccals are likely to be a more informative surrogate tissue than blood for epigenome-wide association studies. Epigenetics 8, 445–454. doi: 10.4161/epi.24362
Lozoya, O. A., Martinez-Reyes, I., Wang, T., Grenet, D., Bushel, P., Li, J., 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. doi: 10.1371/journal.pbio.2005707
Maddocks, O. D. K., Labuschagne, C. F., Adams, P. D., and Vousden, K. H. (2016). Serine metabolism supports the methionine cycle and DNA/RNA methylation through de novo ATP synthesis in cancer cells. Mol. Cell 61, 210–221. doi: 10.1016/j.molcel.2015.12.014
Malik, A. N., and Czajka, A. (2013). Is mitochondrial DNA content a potential biomarker of mitochondrial dysfunction? Mitochondrion 13, 481–492. doi: 10.1016/j.mito.2012.10.011
Manzi, B., Loizzo, A. L., Giana, G., and Curatolo, P. (2008). Autism and metabolic diseases. J. Child Neurol. 23, 307–314. doi: 10.1177/0883073807308698
Marotta, R., Risoleo, M. C., Messina, G., Parisi, L., Carotenuto, M., Vetri, L., et al. (2020). The neurochemistry of Autism. Brain Sci. 10:163. doi: 10.3390/brainsci10030163
Masser, D. R., Clark, N. W., Van Remmen, H., and Freeman, W. 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. 38, 323–333. doi: 10.1007/s11357-016-9930-1
Melnyk, S., Fuchs, G. J., Schulz, E., Lopez, M., Kahler, S. G., Fussell, J. J., et al. (2012). Metabolic imbalance associated with methylation dysregulation and oxidative damage in children with autism. J. Autism Dev. Disord. 42, 367–377. doi: 10.1007/s10803-011-1260-7
Michel, S., Wanet, A., De Pauw, A., Rommelaere, G., Arnould, T., and Renard, P. (2012). Crosstalk between mitochondrial (dys)function and mitochondrial abundance. J. Cell Physiol. 227, 2297–2310. doi: 10.1002/jcp.23021
Montier, L. L. C., Deng, J., and Bai, Y. (2009). Numbers matters: control of mammalian mitochondrial DNA copy number. J. Genet. Genomics 36, 125–131. doi: 10.1016/S1673-8527(08)60099-5
Nardone, S., Sams, D. S., Reuveni, E., Getselter, D., Oron, O., Karpuj, M., et al. (2014). DNA methylation analysis of the autistic brain reveals multiple dysregulated biological pathways. Transl. Psychiatry 4:e433. doi: 10.1038/tp.2014.70
Nardone, S., Sams, D. S., Zito, A., Reuveni, E., and Elliott, E. (2017). Dysregulation of cortical neuron DNA methylation profile in autism spectrum disorder. Cereb. Cortex 27, 5739–5754. doi: 10.1093/cercor/bhx250
Needham, B. D., Adame, M. D., Serena, G., Rose, D. R., Preston, G. M., Conrad, M. C., et al. (2021). Plasma and fecal metabolite profiles in autism spectrum disorder. Biol. Psychiatry 89, 451–462. doi: 10.1016/j.biopsych.2020.09.025
Notarangelo, F. M., and Pocivavsek, A. (2017). Elevated kynurenine pathway metabolism during neurodevelopment: implications for brain and behavior. Neuropharmacology 112(Pt. B), 275–285. doi: 10.1016/j.neuropharm.2016.03.00
Oliveira, G., Diogo, L., Grazina, M., Garcia, P., Ataíde, A., Marques, C., et al. (2005). Mitochondrial dysfunction in autism spectrum disorders: a population-based study. Dev. Med. Child. Neurol. 47, 185–189. doi: 10.1017/S0012162205000332
Pecorelli, A., Ferrara, F., Messano, N., Cordone, V., Schiavone, M. L., Cervellati, F., et al. (2020). Alterations of mitochondrial bioenergetics, dynamics, and morphology support the theory of oxidative damage involvement in autism spectrum disorder. FASEB J. 34, 6521–6538. doi: 10.1096/fj.201902677r
Petersen, M. H., Budtz-Jørgensen, E., Sørensen, S. A., Nielsen, J. E., Hjermind, L. E., Vinther-Jensen, T., et al. (2014). Reduction in mitochondrial DNA copy number in peripheral leukocytes after onset of Huntington’s disease. Mitochondrion 17, 14–21. doi: 10.1016/j.mito.2014.05.001
Petrie, M. A., Sharma, A., Taylor, E. B., Suneja, M., and Shields, R. 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. Genomics 52, 71–80. doi: 10.1152/physiolgenomics.00064.2019
Ploumi, C., Daskalaki, I., and Tavernarakis, N. (2017). Mitochondrial biogenesis and clearance: a balancing act. FEBS J. 284, 183–195. doi: 10.1111/febs.13820
Rauluseviciute, I., Drabløs, F., and Rye, M. B. (2020). DNA hypermethylation associated with upregulated gene expression in prostate cancer demonstrates the diversity of epigenetic regulation. BMC Med. Genomics 13:6. doi: 10.1186/s12920-020-0657-6
Reinecke, C. J., Koekemoer, G., van der Westhuizen, F. H., Louw, R., Lindeque, J. Z., Mienie, L. J., et al. (2012). Metabolomics of urinary organic acids in respiratory chain deficiencies in children. Metabolomics 8, 264–283. doi: 10.1007/s11306-011-0309-0
Ridout, K. K., Coe, J. L., Parade, S. H., Marsit, C. J., Kao, H. T., Porton, B., et al. (2020). Molecular markers of neuroendocrine function and mitochondrial biogenesis associated with early life stress. Psychoneuroendocrinology 116:104632. doi: 10.1016/j.psyneuen.2020.104632
Rylaarsdam, L., and Guemez-Gamboa, A. (2019). Genetic causes and modifiers of autism spectrum disorder. Front. Cell Neurosci. 13:385. doi: 10.3389/fncel.2019.00385
Sandin, S., Lichtenstein, P., Kuja-Halkola, R., Hultman, C., Larsson, H., and Reichenberg, A. (2017). The heritability of autism spectrum disorder. JAMA 318, 1182–1184. doi: 10.1001/jama.2017.12141
Santos, R. X., Correia, S. C., Alves, M. G., Oliveira, P. F., Cardoso, S., Carvalho, C., 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. Acta 1842, 1154–1166. doi: 10.1016/j.bbadis.2014.04.011
Savino, R., Carotenuto, M., Polito, A. N., Di Noia, S., Albenzio, M., Scarinci, A., et al. (2020). Analyzing the potential biological determinants of autism spectrum disorder: from neuroinflammation to the kynurenine pathway. Brain Sci. 10:631. doi: 10.3390/brainsci10090631
Shen, W., Gao, C., Cueto, R., Liu, L., Fu, H., Shao, Y., et al. (2020). Homocysteine-methionine cycle is a metabolic sensor system controlling methylation-regulated pathological signaling. Redox Biol. 28:101322. doi: 10.1016/j.redox.2019.101322
Singh, K., Singh, I. N., Diggins, E., Connors, S. L., Karim, M. A., Lee, D., et al. (2020). Developmental regression and mitochondrial function in children with autism. Ann. Clin. Transl. Neurol. 7, 683–694. doi: 10.1002/acn3.51034
Smirnova, E., Griparic, L., Shurland, D. L., and Van Der Bliek, A. M. (2001). Dynamin-related protein Drp1 is required for mitochondrial division in mammalian cells. Mol. Biol. Cell 12, 2245–2256. doi: 10.1091/mbc.12.8.2245
Smith, A. K., Kilaru, V., Klengel, T., Mercer, K. B., Bradley, B., Conneely, K. 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. 168B, 36–44. doi: 10.1002/ajmg.b.32278
Smith, A. M., Natowicz, M. R., Braas, D., Ludwig, M. A., Ney, D. M., Donley, E. L. R., et al. (2020). A metabolomics approach to screening for autism risk in the children’s autism metabolome project. Autism Res. 13, 1270–1285. doi: 10.1002/aur.2330
Smuts, I., van der Westhuizen, F. H., Louw, R., Mienie, L. J., Engelke, U. F. H., Wevers, R. A., et al. (2013). Disclosure of a putative biosignature for respiratory chain disorders through a metabolomics approach. Metabolomics 9, 379–391. doi: 10.1007/s11306-012-0455-z
Sookoian, S., Rosselli, M. S., Gemma, C., Burgueño, A. L., Gianotti, T. F., Castaño, G. 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. Hepatology 52, 1992–2000. doi: 10.1002/hep.23927
Stathopoulos, S., Gaujoux, R., Lindeque, Z., Mahony, C., Van Der Colff, R., Van Der Westhuizen, F., et al. (2020). DNA methylation associated with mitochondrial dysfunction in a South African autism spectrum disorder cohort. Autism Res. 13, 1079–1093. doi: 10.1002/aur.2310
Stojanovski, D., Koutsopoulos, O. S., Okamoto, K., and Ryan, M. T. (2004). Levels of human Fis1 at the mitochondrial outer membrane regulate mitochondrial morphology. J. Cell Sci. 117, 1201–1210. doi: 10.1242/jcs.01058
St-Pierre, J., Drori, S., Uldry, M., Silvaggi, J. M., Rhee, J., Jäger, S., et al. (2006). Suppression of reactive oxygen species and neurodegeneration by the PGC-1 transcriptional coactivators. Cell 127, 397–408. doi: 10.1016/j.cell.2006.09.024
Sun, J., Longchamps, R. J., Piggott, D. A., Castellani, C. A., Sumpter, J. A., Brown, T. 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. 219, 1285–1293. doi: 10.1093/infdis/jiy658
Tang, G., Rios, P. G., Kuo, S. H., Akman, H. O., Rosoklija, G., Tanji, K., et al. (2013). Mitochondrial abnormalities in temporal lobe of autistic brain. Neurobiol. Dis. 54, 349–361. doi: 10.1016/j.nbd.2013.01.006
Thompson, K., Majd, H., Dallabona, C., Reinson, K., King, M. S., Alston, C. 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. 99, 860–876. doi: 10.1016/j.ajhg.2016.08.014
Tondera, D., Grandemange, S., Jourdain, A., Karbowski, M., Mattenberger, Y., Herzig, S., et al. (2009). SLP-2 is required for stress-induced mitochondrial hyperfusion. EMBO J. 28, 1589–1600. doi: 10.1038/emboj.2009.89
Tremblay, M. W., and Jiang, Y. H. (2019). DNA methylation and susceptibility to autism spectrum disorder. Annu. Rev. Med. 70, 151–166. doi: 10.1146/annurev-med-120417-091431
Valenti, D., de Bari, L., De Filippis, B., Henrion-Caude, A., and Vacca, R. 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. 46, 202–217. doi: 10.1111/j.1469-8749.2005.tb01113.x
Valiente-Pallejà, A., Torrell, H., Muntané, G., Cortés, M. J., Martínez-Leal, R., Abasolo, N., et al. (2018). Genetic and clinical evidence of mitochondrial dysfunction in autism spectrum disorder and intellectual disability. Hum. Mol. Genet. 27, 891–900. doi: 10.1093/hmg/ddy009
Ventura-Clapier, R., Garnier, A., and Veksler, V. (2008). Transcriptional control of mitochondrial biogenesis: the central role of PGC-1α. Cardiovasc. Res. 79, 208–217. doi: 10.1093/cvr/cvn098
Voineagu, I. (2012). Gene expression studies in autism: moving from the genome to the transcriptome and beyond. Neurobiol. Dis. 45, 69–75. doi: 10.1016/j.nbd.2011.07.017
Wai, T., and Langer, T. (2016). Mitochondrial dynamics and metabolic regulation. Trends Endocrinol. Metab. 27, 105–117. doi: 10.1016/j.tem.2015.12.001
Wiśniowiecka-Kowalnik, B., and Nowakowska, B. A. (2019). Genetics and epigenetics of autism spectrum disorder - current evidence in the field. J. Appl. Genet. 60, 37–47. doi: 10.1007/s13353-018-00480-w
Wong, C. C. Y., Meaburn, E. L., Ronald, A., Price, T. S., Jeffries, A. R., Schalkwyk, L. C., et al. (2014). Methylomic analysis of monozygotic twins discordant for autism spectrum disorder and related behavioural traits. Mol. Psychiatry 19, 495–503. doi: 10.1038/mp.2013.41
Wong, C. C. Y., Smith, R. G., Hannon, E., Ramaswami, G., Parikshak, N. N., Assary, E., 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. 28, 2201–2211. doi: 10.1093/hmg/ddz052
Wu, Z., Huang, X., Feng, Y., Handschin, C., Feng, Y., Gullicksen, P. 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. 103, 14379–14384. doi: 10.1073/pnas.0606714103
Yang, X., Xu, S., Qian, Y., He, X., Chen, S., and Xiao, Q. (2020). Hypermethylation of the gene coding for PGC-1α in peripheral blood leukocytes of patients with Parkinson’s disease. Front. Neurosci. 14:97. doi: 10.3389/fnins.2020.00097
Yoo, H. J., Park, M., and Kim, S. 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. Psychiatry 18, 151–156. doi: 10.1080/15622975.2016.1234069
Yoon, S. H., Choi, J., Lee, W. J., and Do, J. T. (2020). Genetic and epigenetic etiology underlying autism spectrum disorder. J. Clin. Med. 9:966. doi: 10.3390/jcm9040966
Yu, X., Qiao, S., Wang, D., Dai, J., Wang, J., Zhang, R., et al. (2016). A metabolomics-based approach for ranking the depressive level in a chronic unpredictable mild stress rat model. RSC Adv. 6, 25751–25765. doi: 10.1039/c6ra00665e
Yudkoff, M. (2017). Interactions in the metabolism of glutamate and the branched-chain amino acids and ketoacids in the CNS. Neurochem. Res. 42, 10–18. doi: 10.1007/s11064-016-2057-z.INTERACTIONS
Zeng, Y., Pan, Q., Wang, X., Li, D., Lin, Y., Man, F., 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. doi: 10.1155/2019/4508762
Zhao, H., Zhao, Y., Ren, Y., Li, M., Li, T., Li, R., 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. 107, 467–474.e5. doi: 10.1016/j.fertnstert.2016.10.039
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.
Edited by:
Alaguraj Veluchamy, St. Jude Children’s Research Hospital, United StatesReviewed by:
Adam R. Wende, University of Alabama at Birmingham, United StatesHashim Islam, University of British Columbia Okanagan, Canada
Rick Schnellmann, University of Arizona, United States
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, Y29sbGVlbi5vcnlhbkB1Y3QuYWMuemE=
†These authors have contributed equally to this work and share first authorship