Abstract
Background Mental disorders are associated with substantially increased morbidity and reduced life expectancy. Accelerated biological ageing might contribute to excess mortality of individuals with mental disorders. The aim of this study was to characterise telomere length, a biological hallmark of ageing, in individuals with mental disorders, and to examine associations between telomere length, age and genetic risk for mental disorders.
Methods The UK Biobank is a multicentre, community-based observational study that recruited >500,000 middle-aged and older adults across England, Scotland and Wales. Average leukocyte telomere length (T/S ratio) was measured using quantitative polymerase chain reaction. Polygenic risk scores (PRS) were calculated for individuals of European ancestry. We estimated differences in T/S ratio and age-related changes in T/S ratio between individuals with anxiety disorder, depression or bipolar disorder and people without mental disorders. We also estimated associations between T/S ratio and PRS for these three disorders.
Results The analyses included up to 308,725 participants. Individuals with depression had shorter telomeres than people without mental disorders (adjusted β = -0.011, 95% CI -0.019 to -0.004, pBonf. = 0.027). There was only limited evidence of case-control differences in telomere length for anxiety disorders or bipolar disorders. Age-related changes in telomere length did not differ between individuals with and without mental disorders. PRS for depression were associated with shorter telomeres (β = -0.006, 95% CI -0.010 to -0.003, pBonf. = 0.001). There was no evidence that PRS for anxiety disorder or bipolar disorder were associated with telomere length.
Conclusion Although telomere length is a biological hallmark of ageing, we observed limited evidence that leukocyte telomere length is a clinically useful marker to quantify accelerated biological ageing in middle-aged and older adults with a lifetime history of anxiety disorder, depression or bipolar disorder.
Introduction
Telomeres are repetitive nucleoprotein complexes at the chromosome ends that play an important role in maintaining genomic stability. Telomeres shorten with each cell division and therefore represent a biological marker of replicative history and cellular age.1,2 Although telomere length is highly heritable,3 age-related attrition results from biological and environmental factors, including lifestyle and chronic stress.4 Telomere attrition has been associated with an increased risk of age-related diseases. Mendelian randomisation analyses in the UK Biobank, a major biomedical database, suggested that telomere length had widespread influence on biomedical traits, disease risk, multiple body systems and life expectancy.5
Individuals with mental disorders have an increased prevalence of age-related diseases and a lower life expectancy.6 They also show signs of accelerated biological ageing, including advanced brain ageing,7 changes in DNA methylation,8 greater levels of inflammation,9 elevated frailty10 and differences in physiological measures such as grip strength.11-13 Telomere length as a molecular marker of biological senescence could provide insight into the relationship between mental health and accelerated biological ageing. Data from a meta-analysis suggested that individuals with anxiety disorders, depressive disorders and posttraumatic stress disorder had shorter telomeres than people without these disorders.14 Findings regarding bipolar disorder have been inconsistent,15 with some studies observing longer telomeres in patients,14 likely due to lithium treatment.16 Bipolar disorder patients not exposed to lithium had shorter telomeres than patients who had been treated with lithium.17
Most previous studies have had limited sample sizes and few studies have included cross-disorder comparisons within the same database. There has also been little exploration of associations between telomere length and genetic risk for mental disorders. The UK Biobank provides an unprecedented resource to investigate health and ageing trajectories, with the world’s largest database of leukocyte telomere length measurements. The aim of this study was to characterise telomere length in individuals with a history of anxiety disorder, depression or bipolar disorder, and to examine associations between telomere length, age and polygenic risk scores for these disorders.
Methods
Study population
The UK Biobank is a prospective study of >500,000 UK residents aged 37–73 at baseline who were recruited between 2006 and 2010. The study rationale and design have been described elsewhere.18 Briefly, individuals registered with the UK National Health Service (NHS) and living within a 25-mile (∼40 km) radius of one of 22 assessment centres were invited to participate. Participants provided data on their sociodemographic characteristics, health behaviours and medical history, underwent physical examination and had blood and urine samples taken. Linked hospital inpatient records are available for most participants and primary care data are available for half of the participants. A third of participants completed an online follow-up mental health questionnaire (MHQ) between 2016 and 2017.
Leukocyte telomere length
Details of the measurement of leukocyte telomere length (UK Biobank data fields 22191 and 22192), including technical adjustments, have been reported elsewhere.19 Briefly, relative telomere length was measured using a validated quantitative polymerase chain reaction (PCR) assay that expresses telomere length as the ratio of the telomere repeat copy number relative to the reference single-copy hemoglobin subunit beta gene (T/S ratio). This ratio is proportional to the individual’s average telomere length.20 Measurements were adjusted for operational and technical parameters (PCR machine, staff member, enzyme batch, primer batch, temperature, humidity, primer batch × PCR machine, primer batch × staff member, A260/A280 ratio of the DNA sample, and A260/A280 ratio squared), loge transformed and Z-standardised. T/S ratio was also converted to base pairs using the formula: base pairs = 3274+2413×((T/S-0.0545)/1.16).21
Mental disorders
We identified individuals with lifetime anxiety disorders, depression or bipolar disorder using our previously reported criteria.11-13 Data sources included the modified Composite International Diagnostic Interview Short Form (CIDI-SF), self-report questions on (hypo)mania and a question on psychiatric diagnoses (field 20544) which were assessed as part of the MHQ, the nurse-led baseline interview in which participants reported medical diagnoses (field 20002), hospital inpatient records (ICD-10 codes), primary care records (Read v2 or CTV3 codes) and self-report questions on mood disorders from the baseline assessment (field 20126). Individuals with psychosis were excluded from all groups and individuals with bipolar disorder were excluded from the anxiety disorder group due to their increased risk of physical multimorbidity.22,23 The depression and bipolar disorder groups were mutually exclusive. Individuals could be included in both the anxiety disorder and the depression group.
Individuals in the control group had no mental disorders: (i) did not report “schizophrenia”, “depression”, “mania / bipolar disorder / manic depression”, “anxiety / panic attacks”, “obsessive compulsive disorder”, “anorexia / bulimia / other eating disorder”, “post-traumatic stress disorder” at the baseline interview; (ii) no psychiatric diagnoses on the MHQ; (iii) no current psychotropic medication use at baseline (field 20003);24 (iv) no ICD-10 Chapter V code in their hospital inpatient record (F20-F99), except for organic causes or substance use; (v) no primary care record containing diagnostic codes for mental disorders;25 (vi) not classified as individuals with probable mood disorder; (vii) no Patient Health Questionnaire-9 (PHQ-9) or Generalised Anxiety Disorder Assessment (GAD-7) sum score of ≥5; (viii) never felt worried, tense, or anxious for most of a month or longer (field 20421); (ix) not identified as cases based on the CIDI-SF and questions on (hypo)manic symptoms.12,13
Genetic quality control
Genetic quality control was performed as described previously.26 Individuals were excluded where recommended by the UK Biobank for unusual levels of missingness (>5%) or heterozygosity.18 Using the genotyped single nucleotide polymorphisms (SNPs), individuals with call rate of less than 98%, who were genetically related to another individual in the dataset (KING r < 0.044, equivalent to removing third-degree relatives and closer)27 or whose self-reported and genotypic sex did not match (X-chromosome homozygosity (FX) < 0.9 for phenotypic males, FX > 0.5 for phenotypic females) were also excluded. To account for familial correlation, removal of relatives was performed using a “greedy” algorithm, which minimises exclusions (for example, by excluding the child in a mother–father–child trio). All analyses were limited to individuals of European ancestry, as defined by 4-means clustering on the first two genetic principal components (PCs) provided by the UK Biobank.28 Principal components analysis was also performed on the European-only subset of the data using FlashPCA2.29
Polygenic risk scores
Polygenic risk scores (PRS) for anxiety disorder, depression and bipolar disorder were calculated using PRSice v.2.30 This method involves calculating PRS as the sum of risk alleles weighted by SNP effect sizes from independent genome-wide association study (GWAS) summary statistics (Supplementary Table 1). Clumping was performed to remove SNPs in high linkage disequilibrium (defined as r2 ≥ 0.1 within 250 kilobases on each side) as linkage disequilibrium can falsely inflate polygenic scores. PRS were calculated at 11 p-value thresholds (5×10−8, 1×10−5, 1×10−3, 0.01, 0.05, 0.1, 0.2, 0.3, 0.4, 0.5 and 1) and the PRS most predictive threshold was selected for regression analyses. All individual-level PRS were standardized prior to analyses.
Covariates
Covariates were identified from previous research and included age, sex, white blood cell count, Townsend deprivation index, physical activity, smoking status, body mass index, body fat percentage and C-reactive protein. For the analyses of PRS, covariates included the first six ancestry-informative population PCs, batch number and assessment centre.
Statistical analyses
Regression analyses were performed in R (version 3.6.2).
Sample characteristics were summarised using means and standard deviations or counts and percentages. Differences in T/S ratio (log z adjusted) between individuals with anxiety disorder, depression or bipolar disorder and the control group were estimated using ordinary least squares regression (± 95% confidence intervals). For these analyses, we fitted minimally adjusted models that included age and sex and fully adjusted models that included all covariates. Finally, associations between T/S ratio (log z adjusted) and PRS for anxiety, depression and bipolar disorder were estimated using ordinary least squares regression. These models included six PCs, batch number and assessment centre.
We calculated adjusted P-values to correct for multiple testing. Two methods were used: (1) Bonferroni and (2) Benjamini & Hochberg,31 all two-tailed, with α = .05 and a false discovery rate of 5%, respectively. P-values were corrected for six to eight tests (see tables for details). We have opted for this approach because the Bonferroni correction may be too conservative and potentially leads to a high number of false negatives.
Additional analyses
We repeated our main analyses (i) with the bipolar disorder group stratified by current lithium use and (ii) to assess the association between PRS independent of diagnosis and treatment-related confounders, we stratified the PRS analyses by case status. As a sensitivity analysis, we excluded individuals with comorbid depression and anxiety disorder. Finally, we repeated the main and sensitivity analyses by antidepressant medication use.
Results
After quality control exclusions and restricting our sample to individuals of European ancestry, 458,078 participants (out of 502,476) had data on both telomere length and polygenic risk scores. We retained up to 308,725 participants with complete data on all covariates. 41,524 individuals had lifetime anxiety disorders, 84,965 had lifetime depression and 3449 had bipolar disorders (Table 1).
Average telomere length (T/S ratio log z adjusted) was slightly greater in individuals with mental disorders relative to the non-psychiatric control group (Supplement Figure 1). However, this difference could be explained by group differences in age and sex. After adjusting for age and sex (Model 1) and other potential confounders (Model 2), we observed that individuals with mental disorders had slightly shorter telomeres (Figure 1). This difference was only statistically significant for the comparison between individuals with depression and the control group (fully adjusted β = -0.011, 95% CI -0.019 to -0.004, pBonf. = 0.027) (Table 2).
When stratifying individuals with bipolar disorder by current lithium use, we found that after adjusting for age and sex, telomeres were shorter in individuals who did not use lithium (adjusted β = -0.045, 95% CI -0.079 to -0.011, pBH = 0.027), relative to the control group. Individuals with bipolar disorder who used lithium had slightly longer telomeres than individuals in the control group (fully adjusted β = 0.121, 95% CI 0.011 to 0.230, pBH = 0.061), although this difference was not statistically significant after multiple testing correction. Individuals with anxiety disorder or depression who reported antidepressant medication use had shorter telomeres than the control group (Supplement Figure 2 and Supplement Table 2).
As expected, telomere length reduced with age (Figure 2). There was no evidence that age-related changes in T/S ratio differed between individuals with mental disorders and people without mental disorders.
The distribution of polygenic risk scores in individuals with mental disorders and the non-psychiatric control group are shown in Supplement Figure 3, confirming that small to moderate increased PRS in individuals with mental disorders. There was little evidence of an association between the PRS for anxiety disorder, depression or bipolar disorder and telomere length (Supplement Figure 4). In a regression model, the PRS for depression was associated with shorter telomeres (adjusted β = -0.006, 95% CI -0.010 to -0.003, pBonf. = 0.001). However, the magnitude of this association was negligible. There was no evidence that the PRS for anxiety disorder or bipolar disorder were associated with telomere length (Figure 3). When stratifying these analyses by case status, the PRS for depression was only associated with shorter telomeres in individuals without mental disorders (Supplement Table 3).
Sensitivity analysis
53,780 individuals had a history of depression without comorbid anxiety disorders, while 14,829 individuals had anxiety disorders without comorbid depression. Individuals with depression had shorter telomeres (fully adjusted β = -0.015, 95% CI -0.025 to -0.006, pBonf. = 0.004) (Supplement Figure 5), and this difference was slightly greater than in the main analysis. There was no evidence of a difference in telomere length between individuals with anxiety disorders and the non-psychiatric control group (fully adjusted β = 0.004, 95% CI -0.012 to 0.021, pBonf. > 0.999) (Supplement Table 4).
Finally, individuals with depression who reported antidepressant medication use had shorter telomeres than individuals with depression not reporting medication use, relative to the control group (Supplement Figure 6 and Supplement Table 5). There was no evidence of an association between telomere length and anxiety disorder status, irrespective of antidepressant medication use.
Discussion
Individuals with depression had slightly shorter telomeres than people without mental disorders. There was only limited evidence that telomere length differed between individuals with anxiety disorder or bipolar disorder and people without mental disorders. Notably, lithium use was associated with elongated telomeres in individuals with bipolar disorder. Antidepressant medication use was associated with reduced telomere length in individuals with depression. Age-related changes in telomere length did not differ between individuals with and without mental disorders. Polygenic risk scores for depression were associated with shorter telomeres. There was no evidence that polygenic risk scores for anxiety disorder or bipolar disorder were associated with telomere length.
The observation that depression was associated with shorter telomeres in the UK Biobank is consistent with data from meta-analyses.14,32,33 Although meta-analyses also provided evidence of an association between anxiety disorders and shorter telomeres,14,34 we did not observe a statistically significant difference between individuals with anxiety disorders and people without mental disorders. This discrepancy could be due to adjustment for multiple potential confounders in the present study. Although many previous studies examining telomere length in individuals with mental disorders adjusted for age and sex, most did not adjust for other confounders. Data from two meta-analyses found no association between bipolar disorder and telomere length.14,35 However, the most recent meta-analysis suggested that bipolar disorder patients had shorter telomeres compared to participants in the control group.36 Inconsistencies between studies could relate to differences in sample characteristics. For example, a recent study found that patients with bipolar disorder type I, but not bipolar disorder type II, had shorter telomeres than healthy controls.37 Another study did not observe group differences in telomere length between bipolar disorder subtypes but was likely underpowered (n=119 vs n=12, respectively).17 The possibility of subgroup differences should be further investigated.
Our finding that lithium use modified the direction of association between bipolar disorder and telomere length is consistent with previous observations that lithium treatment was associated with increased telomere length38,39 and that telomere length positively correlated with duration of lithium treatment16,38. A recent study found that bipolar disorder patients who had never been treated with lithium had shorter telomeres than healthy controls, while patients treated with lithium had longer telomeres than the never treated patients, although not compared to healthy controls.17 Our finding that antidepressant medication use was associated with reduced telomere length in individuals with depression confirms a preliminary study showing that antidepressant use (n=40) was associated with shorter telomeres, independent of depression diagnosis and current depression severity.40
Although previous research suggested that age-related decline in telomere length was greater in individuals with chronic stress or comorbidities,41 we observed similar association between telomere length and age in individuals with mental disorders and the control group.
Multiple studies have examined polygenic scores for telomere length to predict mental disorders,42 but there has been limited research on polygenic risk scores for mental disorders to predict telomere length. Preliminary studies found that unaffected first-degree relatives of individuals with bipolar disorder had shorter telomeres than healthy controls.39,43 Similarly, a small cross-sectional study found that daughters of mothers with depression had shorter telomeres than daughters of never depressed mothers.44 Although these findings suggest that increased genetic risk for mental disorders may affect telomere length, these studies were limited by modest sample sizes and cannot fully disentangle genetic and environmental risk factors. Depression polygenic risk scores were not associated with telomere length or telomere attrition rate in 2032 adults aged 18 to 65 years.42 A study of 290 adults without depression also found no evidence that polygenic risk scores for depression, bipolar disorder or schizophrenia were associated with telomere length.40 We found that polygenic risk scores for depression, but not for anxiety disorder or bipolar disorder, were associated with shorter telomeres, although the strength of this association was negligible. Our finding that depression polygenic risk scores were associated with telomere length only in individuals without mental disorders could be explained by the lower sample size in the case group and warrants replication.
In contrast to most previous studies, we examined associations between telomere length and mental disorders in the same data using a shared control group, allowing for cross-disorder comparison. Our study also included a considerably larger number of individuals with mental disorders. Indeed, the UK Biobank is by far the largest data resource with measured telomere length, which allowed us to adjust for a range of potential confounders.
Our observational study has limitations. The age range was limited to middle-aged and older adults (most between 40 to 69 years old). A previous Dutch study found no difference in telomere length between individuals with current depression aged 60 or older compared to never depressed individuals, suggesting that our findings might not extrapolate to late-life depression.45 For the cross-sectional analyses, our ability to draw causal conclusions was limited. Longitudinal studies have found that major depressive disorder46 or persistent internalizing disorders in men47 predicted reduced telomere length, although not all studies found evidence of a prospective association between mental disorders and telomere length.33,48 A large study of 67,306 individuals aged 20-100 from the Danish general population found no evidence that telomere length predicted depression prospectively or that genetically shorter telomeres predicted depression. Nevertheless, the authors observed that depression was associated with shorter telomeres cross-sectionally, which could be explained by depression causing shorter telomeres or residual confounding.49 Finally, telomeres were measured from leukocyte DNA and findings might differ when examining other tissues. However, research suggests that leukocyte telomeres correlate well with telomere length measured in other tissues.50
Conclusion
Although telomere length is a biological hallmark of ageing, we observed limited evidence that leukocyte telomere length is a clinically useful marker to quantify accelerated biological ageing in middle-aged and older adults with a lifetime history of anxiety disorder, depression or bipolar disorder in the UK Biobank.
Data Availability
The data used are available to all bona fide researchers for health-related research that is in the public interest, subject to an application process and approval criteria. Study materials are publicly available online at http://www.ukbiobank.ac.uk.
Ethics
Ethical approval for the UK Biobank study has been granted by the National Information Governance Board for Health and Social Care and the NHS North West Multicentre Research Ethics Committee (11/NW/0382). No project-specific ethical approval is needed.
Authorship contributions
JM conceived the idea of the study, acquired the data, carried out the statistical analysis, interpreted the findings, wrote the manuscript and revised the manuscript for final submission. CML acquired the studentship funding, interpreted the findings and critically reviewed the manuscript. Both authors read and approved the final manuscript.
Funding and disclosure
JM receives studentship funding from the Biotechnology and Biological Sciences Research Council (BBSRC) (ref: 2050702) and Eli Lilly and Company Limited. CML is part-funded by the National Institute for Health Research (NIHR) Maudsley Biomedical Research Centre at South London and Maudsley NHS Foundation Trust and King’s College London. The views expressed are those of the authors and not necessarily those of the NHS, the NIHR or the Department of Health and Social Care. CML is a member of the Scientific Advisory Board of Myriad Neuroscience.
Data sharing statement
The data used are available to all bona fide researchers for health-related research that is in the public interest, subject to an application process and approval criteria. Study materials are publicly available online at http://www.ukbiobank.ac.uk.
Supplementary material
Supplementary information is available online.
Acknowledgments
This research has been conducted using data from UK Biobank, a major biomedical database. This project made use of time on Rosalind HPC, funded by Guy’s & St Thomas’ Hospital NHS Trust Biomedical Research Centre (GSTT-BRC), South London & Maudsley NHS Trust Biomedical Research Centre (SLAM-BRC), and Faculty of Natural Mathematics & Science (NMS) at King’s College London. Data access permission has been granted under UK Biobank application 45514.