Accelerated epigenetic aging and mitochondrial DNA copy number in bipolar disorder

1Accelerated aging is commonly observed in several chronic illnesses, increasing the rate of aging-related medical conditions and shortening the life span in patients–. In particular, bipolar disorder (BD), an often severe and highly disabling psychiatric condition that affects around 1% of the population, has been characterized by many features of aging that can complicate disease outcomes and converge to premature aging (commonly described as “accelerated aging” in the field),.

2As a multidimensional construct, aging includes physical, psychosocial, and biological changes. The latter are commonly referred to as “senescence”, a concept traditionally used for individual cells to explain their limited capacity to proliferate, but that can also be applied to whole organisms. The study of biological aging (i.e., the molecular mechanisms involved in the age-related decline in physiological functions and components of an organism) in clinical populations generally relies on the measure of relevant “biological clocks”, among which telomere length has been the most commonly investigated. Specifically, studies of telomere length in BD have yielded contrasting results. The majority of findings showed telomere shortening–, while few others revealed either no significant differences compared to controls– or a longer telomere length in BD patients possibly due to the effects of medications. Interestingly, first-degree relatives of BD patients have been recently shown to present shorter telomeres compared to patients and healthy controls,, although multiple aging markers have yet to be measured in this high-risk population. Accordingly, data regarding telomeres in patients with severe mental illness are not robust enough to determine that they are a reliable marker of premature aging. Although several studies reported telomere shortening,, some actually found the opposite direction,, even after controlling for medication exposure. This suggests that novel, more specific aging markers are needed, with a potential approach being an integrative analysis involving multiple markers.

3In this sense, DNA methylation, an epigenetic marker associated with BD pathophysiology and treatment, has been shown to accurately predict chronological age in multiple tissues. This so-called “epigenetic aging” has been recently shown to be accelerated with cumulative lifetime stress, and a faster-running epigenetic clock has been significantly linked to higher mortality risk. Another potential biological clock is the number of copies of mitochondrial DNA (mtDNA), which is tightly associated with mitochondrial function. This marker is of relevance for biological aging as aging has been consistently associated with a progressive dysfunction in respiratory chain activity and cumulative mitochondrial dysfunction. Accordingly, some studies have found significant correlations between mtDNA copy number and chronological age, although not in all populations,. Importantly, while mtDNA copy number and telomere length have been significantly correlated–, recent studies suggest that DNA methylation age and telomere length are uncorrelated and independently predict chronological age,.

4To date, no study has ever investigated the cross talk between these three biological clocks in the same clinical sample. Thus, this study aimed to assess these aging markers in patients with BD and siblings of BD patients, with the ultimate goal of identifying mechanisms involved in the accelerated aging observed in BD and its risk.

5Twenty-two patients with BD I, 16 siblings of BD patients, and 20 healthy controls were enrolled for this study, all matched by age, sex, ethnicity, and race. Subjects were recruited at the University of Texas Health Science Center at Houston, where the local institutional review board approved the protocol. Informed consent was obtained from all participants upon enrollment and prior to any procedure. Participants had no current medical disorder including neurological disorders and traumatic brain injury, schizophrenia, developmental disorders, eating disorders, and intellectual disability. Siblings were enrolled provided that they had at least one relative who met the criteria for BD as determined via a detailed family history assessment, and 14 of them (87.5%) were related to the patients recruited in this study. Healthy controls were excluded if they had a history of any Axis I disorder in first-degree relatives or if they had taken a prescribed psychotropic medication at any point in their lives. Female participants of reproductive age underwent a urine pregnancy test. All participants underwent a urine drug screen to exclude illegal drug use.

6Psychiatric diagnosis of individuals with BD and their siblings was based on the Structured Clinical Interview for the DSM-IV Axis I Disorders (SCID I) interview. All interviews were administered to participants by trained evaluators and later reviewed by a board-certified psychiatrist. Individuals’ current affective state was assessed with the Young Mania Rating Scale (YMRS) and the Montgomery–Åsberg Depression Rating Scale (MADRS). Global functioning was assessed by the Global Assessment of Functioning Scale (GAF).

7Peripheral blood was drawn from fasting subjects in EDTA-containing tubes, followed by DNA isolation from buffy coat with the DNeasy Blood & Tissue Mini Kit (Qiagen, Hilden, Germany), according to the manufacturer’s instructions. Isolated DNA samples were quantified on NanoDrop (Thermo, Waltham, MA, USA) and bisulfite-converted (500 ng) using the EZ DNA Methylation™ Kit (Zymo Research, Irvine, CA, USA), followed by hybridization on the Infinium HumanMethylation450 BeadChip Kit (Illumina, San Diego, CA, USA). Beadchips were scanned in an iScan microarray reader (Illumina), and analyses were performed using the RnBeads R package. Details of the DNA methylation analyses can be found in the Supplementary Methods. Genome-wide DNA methylation levels measured in postmortem cerebellum tissue from BD patients (n = 47) and controls (n = 47) were also extracted from a publicly available data set (GSE38873) at the Gene Expression Omnibus (GEO) platform (ncbi.nlm.nih.gov/geo/) for further analyses, as described below.

8The DNA methylation age (DNAm age) was calculated in blood samples and in the data from GSE38873 using the Horvath age estimation algorithm with a freely available online tool (http://labs.genetics.ucla.edu/horvath/htdocs/dnamage/). The algorithm predicts DNAm age based on the methylation levels of 353 CpGs from the Illumina 450 K Beadchip and is applicable to all types of tissues and cells. The online tool uses an elastic net-penalized regression model (implemented in the R package) and provides an estimate of accelerated epigenetic aging (in years) by regressing DNAm age on chronological age, forming residuals that can be compared between groups.

9Real-time quantitative polymerase chain reactions (PCRs) were performed to measure the amount of mitochondrial DNA with a modified protocol from Tyrka and collaborators, and telomere length as previously described. Please see Supplementary Methods for more details.

10Categorical variables were compared between groups with chi-square tests. Parametric distribution of the continuous variables was tested by Shapiro–Wilk test and histogram visualization. Age, years of education, body mass index, GAF, and cell count estimates (monocytes and CD4 + T-lymphocytes) were analyzed with one-way analysis of variance (ANOVA) followed by post hoc test of Tukey, when appropriate, with a p-value < 0.05 deemed as significant. MADRS, YMRS, and cell count estimates (B-lymphocytes, granulocytes, natural killer cells, and CD8 + T-lymphocytes) showed nonparametric distributions and were therefore analyzed by Kruskal–Wallis tests followed by Mann–Whitney tests with a Bonferroni-corrected p-value < 0.016 deemed as significant. Epigenetic aging acceleration (measured in blood), mtDNA copy number, and telomere length were compared between patients, siblings, and controls by linear regression models using chronological age, sex, and cell count estimates (monocytes, CD4 + T-lymphocytes, B-lymphocytes, granulocytes, natural killer cells, and CD8 + T-lymphocytes) as covariates. Epigenetic aging acceleration in cerebellum samples was compared between BD patients and nonpsychiatric controls by a linear regression model controlling for postmortem interval (PMI). Correlations between variables were tested with either Pearson or Spearman tests, depending on their distribution.

11Controls, BD patients, and siblings showed no differences for age, sex, ethnicity, race, years of education, smoking status, or body mass index (Table 1). Chronological age ranged between 22 (youngest) and 51 (oldest) years old among controls, 20 and 51 among patients, and 19 and 58 among siblings (p = 0.325). Patients showed significantly higher scores on the MADRS and YMRS compared to controls (although of mild severity, attesting to their euthymic status), as well as lower functioning scores, as assessed by the GAF compared to siblings and controls. In total 10 out of the 16 siblings enrolled were asymptomatic (unaffected) and presented with no current or lifetime history of mental illness. The other six were diagnosed with alcohol use disorder (n = 3), generalized anxiety disorder (n = 1), specific phobia (n = 1), and posttraumatic stress disorder (n = 1). Seventy-two percent of the patients were on medications (Table 1), while none of the siblings or controls were medicated at the time of enrollment.

12We initially performed a preliminary methylome analysis looking for differences between controls, BD patients, and siblings in individual CpG probes and regions. No locus or regions withstood false-discovery rate correction for multiple testing (adjusted p-value > 0.05 for all comparisons), possibly due to a lack of statistical power. Considering the cell-specific patterns of methylation, we also used the methylation data to predict the cell-type composition in blood. No differences were found in the predicted frequency of B-lymphocytes, granulocytes, monocytes, natural killer cells, CD4+ T-lymphocytes, or CD8+ T-lymphocytes between groups (Supplementary Table S1).

13DNAm age and a measure of accelerated epigenetic aging were calculated for each individual based on the 450K data. As expected, there was a strong positive correlation between individuals’ DNAm age and chronological age (r = 0.944, p < 0.001; Fig. 1). Accelerated epigenetic aging did not differ between men and women, either when analyzing the entire sample or when checking for differences within each diagnosis group (p < 0.05). When analyzing the entire sample, accelerated epigenetic aging showed no difference between groups after controlling for chronological age, sex, and blood cell-type composition (F(2, 46) = 2.043, p = 0.141; η2 p = 0.082) (Fig. 2a). Nonetheless, based on previous findings showing a cumulative effect of lifetime stress on epigenetic aging, we sought to investigate whether a difference would be detected when limiting our analysis specifically to older subjects. We therefore split our sample into two smaller groups of younger and older subjects based on the median of the chronological age (33 years old) and performed our analyses a second time. As hypothesized, accelerated epigenetic aging was significantly different between groups in the older subsample (F(2, 17) = 3.924, p = 0.04, and η2 p = 0.316, Fig. 2b), with BD patients presenting a higher accelerated aging compared to controls (pairwise comparison’s Bonferroni-corrected p = 0.047). Siblings did not significantly differ from controls (p = 0.159). No difference between groups was detected among the younger subjects (F(2,18) = 0.095, p = 0.910, and η2 p = 0.010).

14We used publicly available DNA methylation data from postmortem cerebellum samples from BD patients and controls (GSE38873) to check whether accelerated epigenetic aging would be detected in brain tissue, as well. Specifically, several cerebellar changes have been reported in patients with BD, including decreased cerebellar volume, cerebellar atrophy, aberrant cerebellar connectivity, and microstructure abnormalities. Age, DNAm age, and postmortem interval (PMI) of the samples analyzed are shown in Supplementary Table S2. As expected, DNAm age and chronological age were significantly correlated (r = 0.830, p < 0.001, Supplementary Fig. S1A). We found no statistically significant difference in accelerated epigenetic aging between patients and controls when analyzing the entire sample and controlling for PMI (F(1, 91) = 1.854, p = 0.177, and η2 p = 0.020; Supplementary Fig. S1B). We further split the sample into younger and older groups based on the median age (45 years old) and reran the analysis between groups, as similarly performed with the blood samples, but no significant differences were found within the younger (F(1, 34) = 0.063, p = 0.803, and η2 p = 0.002, Supplementary Fig. S1C) or older individuals (F(1, 54) = 2.894, p = 0.095, and η2 p = 0.051, Supplementary Fig. S1D).

15No correlation was observed between mtDNA copy number and chronological age (r s = 0.054, p = 0.693), and a significantly higher mtDNA copy number index was observed in women compared to men (U = 172, p = 0.003). When analyzing the entire sample controlling for chronological age, sex, and cell count estimates, mtDNA copy number was significantly different between groups (F(2, 43) = 3.963, p = 0.026, and η2 p = 0.156, Fig. 2c). Specifically, BD patients presented a higher mtDNA copy number index compared to controls (p = 0.04), while controls and siblings did not show any difference (p > 0.05). When we compared mtDNA copy number within younger and older subjects (median 33 years old, as performed for the epigenetic clock analysis), we found no significant differences between groups (younger subjects—F(2, 18) = 0.845, p = 0.446, and η2 p = 0.086; older subjects— F(2, 14) = 3.570, p = 0.056, and η2 p = 0.338).

16Although mtDNA copy number and the accelerated epigenetic aging score were not significantly correlated when considering the entire sample (r s = 0.038, p = 0.780, Fig. 2e), they were positively correlated within the older group (r s = 0.697, p < 0.001, Fig. 2f). As an exploratory analysis to further explore this association and based on previous evidence of DNA methylation control of nuclear-encoded mitochondrial genes (i.e., defined as genes coding for mitochondrial proteins that are encoded in the nuclear DNA), we checked whether some of the genes that compose the epigenetic clock included mitochondrial genes. Our analyses showed that out of the 353 clock CpGs, 17 of them are located within nuclear-encoded mitochondrial genes (according to the latest update on the Human MitoCarta2.0 data set) (Table S3). Further details of these genes, including their exact genomic location and annotation, can be found in Supplementary Table S3.

17Contrary to our hypothesis, we did not find a significant correlation between telomere length and chronological age in our sample (r s = –0.117, p = 0.394). We also found no statistical differences in telomere length between men and women (U = 269, p = 0.203). The groups did not show any significant differences with regard to the telomere length when controlling for chronological age, sex, and cell count estimates (F(2, 43) = 0.486, p = 0.619, and η2 p = 0.022, Fig. 3). Telomere length did not significantly correlate with epigenetic accelerated aging (r s = 0.005, p = 0.971) or mtDNA copy number (r s = 0.197, p = 0.148) either. These correlations remained to be not significant even when we split our sample into younger and older subjects (younger subjects—F(2, 18) = 0.692, p = 0.513, and η2 p = 0.071; older subjects—F(2, 14) = 1.222, p = 0.324, and η2 p = 0.149).

18Patients with BD presented significantly lower global functioning scores as assessed by the GAF scale compared to both controls and siblings (F(2, 47) = 71.382, p < 0.001, and η2 p = 0.752, Table 1). Based on this, we sought to investigate whether global functioning would be associated with the aging markers previously measured in our sample. GAF scores showed a positive and moderately strong correlation with mtDNA copy number (r s = –0.364, p = 0.013), and a higher mtDNA copy number was associated with a greater number of episodes within the BD group (r s = 0.353, p = 0.02). By contrast, GAF scores did not correlate with chronological age (r = 0.032, p = 802), epigenetic age acceleration (r s = –0.024, p = 0.872), or telomere length (r s = –0.017, p = 0.910). Telomere length did not show any significant correlations with clinical variables, including illness duration (r = –0.206, p = 0.383) or total number of episodes (r = 0.068, p = 0.777). Further correlations between mtDNA copy number and epigenetic age with other clinical variables (mood symptoms scores, age of illness onset, and length of illness) did not reach statistical significance (p > 0.05).

19Given the dearth of findings on the role of accelerated biological aging in the pathophysiology and progression of complex psychiatric disorders, this study aimed to compare multiple biological clocks in BD patients, siblings, and healthy controls. The most compelling finding of this study is that patients displayed an accelerated epigenetic aging that seems to be at least partly associated with an increase in mtDNA copy number. Furthermore, both of these aging markers were independent of telomere shortening. To our knowledge, this is the first study to examine epigenetic aging and its association with other biological aging markers in BD.

20The association between DNA methylation and aging has been suggested by several studies, and so has the use of epigenetic clocks in predicting chronological age,. In line with previous results, our current findings suggest that DNAm age might represent a better aging marker than the commonly used telomere length, particularly in older individuals. By being dynamically modulated, epigenetic aging has been shown to be influenced by a plethora of lifestyle factors and to accelerate with the cumulative effects of life stressors. This might be one of the reasons why we only detected significant differences in epigenetic aging in older subjects, especially considering the key role proposed for chronic stress and the high rates of unhealthy lifestyles seen in BD patients,.

21Our findings suggest that accelerated epigenetic aging in BD may be at least partially associated with an increase in mtDNA copy number. Studies specifically addressing mtDNA levels in BD show conflicting results. While some authors report no differences in patients–, others report both decreased and increased levels compared to controls. Although mtDNA copy number has been shown to decline with age in specific populations, its content is thought to slightly increase in the fifth decade of life and later decline in older subjects. Moreover, Tyrka and collaborators have recently shown that early life stress is associated with an increase in mtDNA copy number, which has been hypothesized as a compensatory mechanism to the aging-related increases in energy demand or reduced mitochondrial function. A similar interpretation could be adopted for our findings based on the mitochondrial dysfunction consistently reported for BD patients, and the strong association between BD and early life stress.

22To the best of our knowledge, this is the first study to suggest a link between mtDNA copy number and the epigenetic clock. Accordingly, previous studies have reported an important mitochondrial–nuclear cross talk and an association between DNA methylation and mtDNA,. Nuclear DNA methylation significantly interferes with mitochondrial function since most mitochondrial proteins are encoded by the nuclear genome. For instance, polymerase γ catalytic subunit, which is a mitochondrial protein encoded in the nucleus, is known to modulate mtDNA replication, and its methylation status has been negatively correlated with mtDNA copy number. In particular, our exploratory analysis showed that 17 CpGs included in the clock CpGs used to predict DNAm age are actually located within nuclear-encoded mitochondrial genes. Even though this finding does not confirm mitochondrial function as a top mechanism related to the epigenetic clock, it does suggest an interesting loop of modulation between nuclear DNA methylation and mitochondrial function that might underlie the accelerated aging process observed in our sample.

23Contrary to our hypothesis, we were unable to find a significant correlation between telomere length and chronological age in our sample. This suggests a nonlinear association between aging and telomere shortening, at least in specific populations. Accordingly, a recent study in a large cohort showed that telomere length declines over midlife but remains stable (and actually subtly increases) in later life. In addition, some medications used for BD treatment have been shown to increase telomere length, including lithium,. Along this line, our sample was medicated and 23.8% took lithium. Finally, our results showing no differences in telomere length are consistent with a recent meta-analysis that suggests that telomere shortening is not seen in all BD samples.

24Our findings suggest independent mechanisms regulating the association between telomere length with the epigenetic clock and mtDNA copy number. In fact, telomere length and epigenetic clock have been previously shown to be uncorrelated and to independently predict chronological age,. Furthermore, some authors speculate that epigenetic aging and cellular senescence are independent mechanisms that do not necessarily happen together. Indeed, DNAm age does not measure mitotic age or cellular senescence and has been hypothesized to contain information that is complementary to that of the telomere clock. Of note, the lack of association between telomere length and mtDNA copy number in our sample contradicts recent evidence of a cross talk between both markers,. This finding may indicate the presence of mechanisms specific to our population of medicated BD patients.

25There are several possible explanations for the accelerated aging seen in our sample of BD patients. So far, it is unclear whether such an accelerated epigenetic aging is specifically related to BD’s genetic makeup or to the direct effects of environmental parameters, such as general medical comorbidities, use of tobacco and alcohol, obesity, exposure to different medications, toxins, and lifetime-stressful events, among others. All of these could be at least partly linked to an increase in oxidative stress (as consistently reported in BD patients,), which might drive an aging acceleration through damage accumulation. The lack of differences in siblings (individuals at familial risk for BD) suggests that genetic predisposition alone may not be enough to determine epigenetic accelerated aging, at least in our particular sample. The most likely hypothesis is that both genetic background and environmental exposure interact to determine such acceleration. Also, BD patients have been reported to present a dysfunctional HPA axis and therefore be less resilient to stress, which could potentially feed a vicious cycle of stress-induced epigenetic aging acceleration.

26Finally, contrary to our hypothesis, we did not find a statistically significant epigenetic aging acceleration in the cerebellum of patients with BD. In accordance with these results, a recent study found no evidence of accelerated brain aging in BD based on the analysis of aging-related imaging structural trajectories. Nonetheless, a previous study assessing the epigenetic clock has shown that the cerebellum ages more slowly than other tissues possibly due to mechanisms involving RNA helicases. Based on this, we cannot rule out the possibility that epigenetic acceleration may be detectable in other brain regions from patients with BD that are more tightly related to mood symptoms than the cerebellum, such as the prefrontal cortex or the hippocampus.

27One of the limitations of this study is our small sample size, which precluded us from analyzing the effects of medications and other important demographic variables in our results, including the effects of comorbidities and the different diagnoses presented by some of the siblings included in the study. Also, the cross-sectional design of this study did not allow us to directly assess aging ‘acceleration’ per se, which will only be confirmed by longitudinal studies. Moreover, information on participants’ history of lifetime trauma exposure and early life stress was not available, which limited the analysis of stress influences on age acceleration in our sample. In this sense, childhood trauma and stress may play a role in the mechanisms regulating the association between telomere length and mtDNA copy number, and have been associated with a shorter telomere length and disease risk. Also, medication use and its heterogeneity among patients may represent potential important confounders in our results.

28In summary, the findings of this study provide novel evidence of accelerated epigenetic aging in BD patients. We showed that accelerated epigenetic aging measured in peripheral blood is partly associated with a higher mtDNA copy number, but is independent of telomere shortening. Future longitudinal studies focused on the mechanisms underlying the link between accelerated aging and vulnerability to BD are warranted, as well as the assessment of aging markers in other tissues of relevance to the pathophysiology of BD. The complex interplay between different aging markers in BD highlighted in this study provides evidence that the epigenetic clock may be a relevant novel target for future treatment and prophylactic interventions for BD.