Meta-analysis of epigenetic aging in schizophrenia reveals multifaceted relationships with age, sex, illness duration, and polygenic risk

BACKGROUND

The study of biological age acceleration may help identify at-risk individuals and reduce the rising global burden of age-related diseases. Using DNA methylation (DNAm) clocks, we investigated biological aging in schizophrenia (SCZ), a mental illness that is associated with an increased prevalence of age-related disabilities and morbidities. In a whole blood DNAm sample of 1090 SCZ cases and 1206 controls across four European cohorts, we performed a meta-analysis of differential aging using three DNAm clocks (i.e., Hannum, Horvath, and Levine). To dissect how DNAm aging contributes to SCZ, we integrated information on duration of illness and SCZ polygenic risk, as well as stratified our analyses by chronological age and biological sex.

RESULTS

We found that blood-based DNAm aging is significantly altered in SCZ independent from duration of the illness since onset. We observed sex-specific and nonlinear age effects that differed between clocks and point to possible distinct age windows of altered aging in SCZ. Most notably, intrinsic cellular age (Horvath clock) is decelerated in SCZ cases in young adulthood, while phenotypic age (Levine clock) is accelerated in later adulthood compared to controls. Accelerated phenotypic aging was most pronounced in women with SCZ carrying a high polygenic burden with an age acceleration of + 3.82 years (CI 2.02-5.61, P = 1.1E-03). Phenotypic aging and SCZ polygenic risk contributed additively to the illness and together explained up to 14.38% of the variance in disease status.

CONCLUSIONS

Our study contributes to the growing body of evidence of altered DNAm aging in SCZ and points to intrinsic age deceleration in younger adulthood and phenotypic age acceleration in later adulthood in SCZ. Since increased phenotypic age is associated with increased risk of all-cause mortality, our findings indicate that specific and identifiable patient groups are at increased mortality risk as measured by the Levine clock. Our study did not find that DNAm aging could be explained by the duration of illness of patients, but we did observe age- and sex-specific effects that warrant further investigation. Finally, our results show that combining genetic and epigenetic predictors can improve predictions of disease outcomes and may help with disease management in schizophrenia.

Longitudinal analyses of depression, anxiety, and suicidal ideation highlight greater prevalence in the northern Dutch population during the COVID-19 lockdowns

BACKGROUND

The pandemic of the coronavirus disease 2019 (COVID-19) has led to an increased burden on mental health.

AIMS

To investigate the development of major depressive disorder (MDD), generalized anxiety disorder (GAD), and suicidal ideation in the Netherlands during the first fifteen months of the pandemic and three nation-wide lockdowns.

METHOD

Participants of the Lifelines Cohort Study -a Dutch population-based sample-reported current symptoms of MDD and GAD, including suicidal ideation, according to DSM-IV criteria. Between March 2020 and June 2021, 36,106 participants (aged 18-96) filled out a total of 629,811 questionnaires across 23 time points. Trajectories over time were estimated using generalized additive models and analyzed in relation to age, sex, and lifetime history of MDD/GAD.

RESULTS

We found non-linear trajectories for MDD and GAD with a higher number of symptoms and prevalence rates during periods of lockdown. The point prevalence of MDD and GAD peaked during the third hard lockdown at 2.88 % (95 % CI: 2.71 %-3.06 %) and 2.92 % (95 % CI: 2.76 %-3.08 %), respectively, in March 2021. Women, younger adults, and participants with a history of MDD/GAD reported significantly more symptoms. For suicidal ideation, we found a significant linear increase over time in younger participants. For example, 20-year-old participants reported 4.14× more suicidal ideation at the end of June 2021 compared to the start of the pandemic (4.64 % (CI: 3.09 %-6.96 %) versus 1.12 % (CI: 0.76 %-1.66 %)).

LIMITATIONS

Our findings should be interpreted in relation to the societal context of the Netherlands and the public health response of the Dutch government during the pandemic, which may be different in other regions in the world.

CONCLUSIONS

Our study showed greater prevalence of MDD and GAD during COVID-19 lockdowns and a continuing increase in suicidal thoughts among young adults suggesting that the pandemic and government enacted restrictions impacted mental health in the population. Our findings provide actionable insights on mental health in the population during the pandemic, which can guide policy makers and clinical care during future lockdowns and epi/pandemics.

An Ethical Framework for Research Using Genetic Ancestry

A wide range of research uses patterns of genetic variation to infer genetic similarity between individuals, typically referred to as genetic ancestry. This research includes inference of human demographic history, understanding the genetic architecture of traits, and predicting disease risk. Researchers are not just structuring an intellectual inquiry when using genetic ancestry, they are also creating analytical frameworks with broader societal ramifications. This essay presents an ethics framework in the spirit of virtue ethics for these researchers: rather than focus on rule following, the framework is designed to build researchers' capacities to react to the ethical dimensions of their work. The authors identify one overarching principle of intellectual freedom and responsibility, noting that freedom in all its guises comes with responsibility, and they identify and define four principles that collectively uphold researchers' intellectual responsibility: truthfulness, justice and fairness, anti-racism, and public beneficence. Researchers should bring their practices into alignment with these principles, and to aid this, the authors name three common ways research practices infringe these principles, suggest a step-by-step process for aligning research choices with the principles, provide rules of thumb for achieving alignment, and give a worked case. The essay concludes by identifying support needed by researchers to act in accord with the proposed framework.

Significant variation in the performance of DNA methylation predictors across data preprocessing and normalization strategies

Background

DNA methylation (DNAm)-based predictors hold great promise to serve as clinical tools for health interventions and disease management. While these algorithms often have high prediction accuracy, the consistency of their performance remains to be determined. We therefore conduct a systematic evaluation across 101 different DNAm data preprocessing and normalization strategies and assess how each analytical strategy affects the consistency of 41 DNAm-based predictors.

Results

Our analyses are conducted in a large EPIC DNAm array dataset from the Jackson Heart Study (N = 2053) that included 146 pairs of technical replicate samples. By estimating the average absolute agreement between replicate pairs, we show that 32 out of 41 predictors (78%) demonstrate excellent consistency when appropriate data processing and normalization steps are implemented. Across all pairs of predictors, we find a moderate correlation in performance across analytical strategies (mean rho = 0.40, SD = 0.27), highlighting significant heterogeneity in performance across algorithms. Successful or unsuccessful removal of technical variation furthermore significantly impacts downstream phenotypic association analysis, such as all-cause mortality risk associations.

Conclusions

We show that DNAm-based algorithms are sensitive to technical variation. The right choice of data processing strategy is important to achieve reproducible estimates and improve prediction accuracy in downstream phenotypic association analyses. For each of the 41 DNAm predictors, we report its degree of consistency and provide the best performing analytical strategy as a guideline for the research community. As DNAm-based predictors become more and more widely used, our work helps improve their performance and standardize their implementation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02793-w.

Using symptom-based case predictions to identify host genetic factors that contribute to COVID-19 susceptibility

Epidemiological and genetic studies on COVID-19 are currently hindered by inconsistent and limited testing policies to confirm SARS-CoV-2 infection. Recently, it was shown that it is possible to predict COVID-19 cases using cross-sectional self-reported disease-related symptoms. Here, we demonstrate that this COVID-19 prediction model has reasonable and consistent performance across multiple independent cohorts and that our attempt to improve upon this model did not result in improved predictions. Using the existing COVID-19 prediction model, we then conducted a GWAS on the predicted phenotype using a total of 1,865 predicted cases and 29,174 controls. While we did not find any common, large-effect variants that reached genome-wide significance, we do observe suggestive genetic associations at two SNPs (rs11844522, p = 1.9x10-7; rs5798227, p = 2.2x10-7). Explorative analyses furthermore suggest that genetic variants associated with other viral infectious diseases do not overlap with COVID-19 susceptibility and that severity of COVID-19 may have a different genetic architecture compared to COVID-19 susceptibility. This study represents a first effort that uses a symptom-based predicted phenotype as a proxy for COVID-19 in our pursuit of understanding the genetic susceptibility of the disease. We conclude that the inclusion of symptom-based predicted cases could be a useful strategy in a scenario of limited testing, either during the current COVID-19 pandemic or any future viral outbreak.

The association of sex, age and FKBP5 genotype with common somatic symptoms: A replication study in the lifelines cohort study

OBJECTIVE

Our aim was to replicate a recent study that reported an association between the rs9470080 CC-genotype and common somatic symptoms in women, but not in men. Additionally, we quantified the genetic contribution to phenotypic variation in common somatic symptom levels.

METHODS

We used data from the Lifelines Cohort Study, including 28,299 participants (60.0% female; 44.2% CC-genotype; mean age 42.9 (14.2) years). Common somatic symptoms were measured with the SCL-90 SOM subscale. To assess the association between the rs9470080 genotype and SCL-90 SOM scores we applied similar analyses as the original study, including independent t-tests, two-way ANOVAs and a mixed ANOVA. To estimate the proportion of phenotypic variance in SCL-90 SOM scores explained by single nucleotide polymorphisms (SNPs), we used a genomic-relatedness-based restricted maximum-likelihood method.

RESULTS

We could not replicate the original study's findings. We found no association between the rs9470080 genotype and common somatic symptom levels in either female or male participants (F(1, 8775) = 1.07, p = 0.30 and F(1,13,903) = 0.01, p = 0.93, respectively). Genome-wide heritability analyses show that 12.1% (p = 2.1e-08) of the phenotypic variance in common somatic symptom levels in Lifelines can be explained by SNPs. The genetic contribution to common somatic symptom levels was higher in male participants (SNP-h² = 20.5%; p = 9.1e-08) than in female participants (SNP-h² = 12.0%, p = 2.8e-05).

CONCLUSION

Our findings of significant SNP-h² and the sex-specific differences herein, does warrant further sex-stratified research of individual genetic variants associated with common somatic symptoms. Preferably, further research should be performed within the analytic framework of a genome-wide association study.

Accelerated aging in the brain, epigenetic aging in blood, and polygenic risk for schizophrenia

Schizophrenia patients show signs of accelerated aging in cognitive and physiological domains. Both schizophrenia and accelerated aging, as measured by MRI brain images and epigenetic clocks, are correlated with increased mortality. However, the association between these aging measures have not yet been studied in schizophrenia patients. In schizophrenia patients and healthy subjects, accelerated aging was assessed in brain tissue using a longitudinal MRI (N = 715 scans; mean scan interval 3.4 year) and in blood using two epigenetic age clocks (N = 172). Differences ('gaps') between estimated ages and chronological ages were calculated, as well as the acceleration rate of brain aging. The correlations between these aging measures as well as with polygenic risk scores for schizophrenia (PRS; N = 394) were investigated. Brain aging and epigenetic aging were not significantly correlated. Polygenic risk for schizophrenia was significantly correlated with brain age gap, brain age acceleration rate, and negatively correlated with DNAmAge gap, but not with PhenoAge gap. However, after controlling for disease status and multiple comparisons correction, these effects were no longer significant. Our results imply that the (accelerated) aging observed in the brain and blood reflect distinct biological processes. Our findings will require replication in a larger cohort.

Lifelines COVID-19 cohort: investigating COVID-19 infection and its health and societal impacts in a Dutch population-based cohort

PURPOSE

The Lifelines COVID-19 cohort was set up to assess the psychological and societal impacts of the COVID-19 pandemic and investigate potential risk factors for COVID-19 within the Lifelines prospective population cohort.

PARTICIPANTS

Participants were recruited from the 140 000 eligible participants of Lifelines and the Lifelines NEXT birth cohort, who are all residents of the three northern provinces of the Netherlands. Participants filled out detailed questionnaires about their physical and mental health and experiences on a weekly basis starting in late March 2020, and the cohort consists of everyone who filled in at least one questionnaire in the first 8 weeks of the project. FINDINGS TO DATE: >71 000 unique participants responded to the questionnaires at least once during the first 8 weeks, with >22 000 participants responding to seven questionnaires. Compiled questionnaire results are continuously updated and shared with the public through the Corona Barometer website. Early results included a clear signal that younger people living alone were experiencing greater levels of loneliness due to lockdown, and subsequent results showed the easing of anxiety as lockdown was eased in June 2020.

FUTURE PLANS

Questionnaires were sent on a (bi)weekly basis starting in March 2020 and on a monthly basis starting July 2020, with plans for new questionnaire rounds to continue through 2020 and early 2021. Questionnaire frequency can be increased again for subsequent waves of infections. Cohort data will be used to address how the COVID-19 pandemic developed in the northern provinces of the Netherlands, which environmental and genetic risk factors predict disease susceptibility and severity and the psychological and societal impacts of the crisis. Cohort data are linked to the extensive health, lifestyle and sociodemographic data held for these participants by Lifelines, a 30-year project that started in 2006, and to data about participants held in national databases.

Genetic and clinical analyses of psychosis spectrum symptoms in a large multiethnic youth cohort reveal significant link with ADHD

Psychotic symptoms are not only an important feature of severe neuropsychiatric disorders, but are also common in the general population, especially in youth. The genetic etiology of psychosis symptoms in youth remains poorly understood. To characterize genetic risk for psychosis spectrum symptoms (PS), we leverage a community-based multiethnic sample of children and adolescents aged 8-22 years, the Philadelphia Neurodevelopmental Cohort (n = 7225, 20% PS). Using an elastic net regression model, we aim to classify PS status using polygenic scores (PGS) based on a range of heritable psychiatric and brain-related traits in a multi-PGS model. We also perform univariate PGS associations and evaluate age-specific effects. The multi-PGS analyses do not improve prediction of PS status over univariate models, but reveal that the attention deficit hyperactivity disorder (ADHD) PGS is robustly and uniquely associated with PS (OR 1.12 (1.05, 1.18) P = 0.0003). This association is driven by subjects of European ancestry (OR = 1.23 (1.14, 1.34), P = 4.15 × 10^-7) but is not observed in African American subjects (P = 0.65). We find a significant interaction of ADHD PGS with age (P = 0.01), with a stronger association in younger children. The association is independent of phenotypic overlap between ADHD and PS, not indirectly driven by substance use or childhood trauma, and appears to be specific to PS rather than reflecting general psychopathology in youth. In an independent sample, we replicate an increased ADHD PGS in 328 youth at clinical high risk for psychosis, compared to 216 unaffected controls (OR 1.06, CI(1.01, 1.11), P = 0.02). Our findings suggest that PS in youth may reflect a different genetic etiology than psychotic symptoms in adulthood, one more akin to ADHD, and shed light on how genetic risk can be investigated across early disease trajectories.

Whole blood transcriptome analysis in bipolar disorder reveals strong lithium effect

BACKGROUND

Bipolar disorder (BD) is a highly heritable mood disorder with complex genetic architecture and poorly understood etiology. Previous transcriptomic BD studies have had inconsistent findings due to issues such as small sample sizes and difficulty in adequately accounting for confounders like medication use.

METHODS

We performed a differential expression analysis in a well-characterized BD case-control sample (Nsubjects = 480) by RNA sequencing of whole blood. We further performed co-expression network analysis, functional enrichment, and cell type decomposition, and integrated differentially expressed genes with genetic risk.

RESULTS

While we observed widespread differential gene expression patterns between affected and unaffected individuals, these effects were largely linked to lithium treatment at the time of blood draw (FDR < 0.05, Ngenes = 976) rather than BD diagnosis itself (FDR < 0.05, Ngenes = 6). These lithium-associated genes were enriched for cell signaling and immune response functional annotations, among others, and were associated with neutrophil cell-type proportions, which were elevated in lithium users. Neither genes with altered expression in cases nor in lithium users were enriched for BD, schizophrenia, and depression genetic risk based on information from genome-wide association studies, nor was gene expression associated with polygenic risk scores for BD.

CONCLUSIONS

These findings suggest that BD is associated with minimal changes in whole blood gene expression independent of medication use but emphasize the importance of accounting for medication use and cell type heterogeneity in psychiatric transcriptomic studies. The results of this study add to mounting evidence of lithium's cell signaling and immune-related mechanisms.

Contribution of common and rare variants to bipolar disorder susceptibility in extended pedigrees from population isolates

Current evidence from case/control studies indicates that genetic risk for psychiatric disorders derives primarily from numerous common variants, each with a small phenotypic impact. The literature describing apparent segregation of bipolar disorder (BP) in numerous multigenerational pedigrees suggests that, in such families, large-effect inherited variants might play a greater role. To identify roles of rare and common variants on BP, we conducted genetic analyses in 26 Colombia and Costa Rica pedigrees ascertained for bipolar disorder 1 (BP1), the most severe and heritable form of BP. In these pedigrees, we performed microarray SNP genotyping of 838 individuals and high-coverage whole-genome sequencing of 449 individuals. We compared polygenic risk scores (PRS), estimated using the latest BP1 genome-wide association study (GWAS) summary statistics, between BP1 individuals and related controls. We also evaluated whether BP1 individuals had a higher burden of rare deleterious single-nucleotide variants (SNVs) and rare copy number variants (CNVs) in a set of genes related to BP1. We found that compared with unaffected relatives, BP1 individuals had higher PRS estimated from BP1 GWAS statistics (P = 0.001 ~ 0.007) and displayed modest increase in burdens of rare deleterious SNVs (P = 0.047) and rare CNVs (P = 0.002 ~ 0.033) in genes related to BP1. We did not observe rare variants segregating in the pedigrees. These results suggest that small-to-moderate effect rare and common variants are more likely to contribute to BP1 risk in these extended pedigrees than a few large-effect rare variants.

Genome-wide association study identifies 30 loci associated with bipolar disorder

Bipolar disorder is a highly heritable psychiatric disorder. We performed a genome-wide association study (GWAS) including 20,352 cases and 31,358 controls of European descent, with follow-up analysis of 822 variants with P < 1 × 10-4 in an additional 9,412 cases and 137,760 controls. Eight of the 19 variants that were genome-wide significant (P < 5 × 10-8) in the discovery GWAS were not genome-wide significant in the combined analysis, consistent with small effect sizes and limited power but also with genetic heterogeneity. In the combined analysis, 30 loci were genome-wide significant, including 20 newly identified loci. The significant loci contain genes encoding ion channels, neurotransmitter transporters and synaptic components. Pathway analysis revealed nine significantly enriched gene sets, including regulation of insulin secretion and endocannabinoid signaling. Bipolar I disorder is strongly genetically correlated with schizophrenia, driven by psychosis, whereas bipolar II disorder is more strongly correlated with major depressive disorder. These findings address key clinical questions and provide potential biological mechanisms for bipolar disorder.

Transcriptome analysis in whole blood reveals increased microbial diversity in schizophrenia

The role of the human microbiome in health and disease is increasingly appreciated. We studied the composition of microbial communities present in blood across 192 individuals, including healthy controls and patients with three disorders affecting the brain: schizophrenia, amyotrophic lateral sclerosis, and bipolar disorder. By using high-quality unmapped RNA sequencing reads as candidate microbial reads, we performed profiling of microbial transcripts detected in whole blood. We were able to detect a wide range of bacterial and archaeal phyla in blood. Interestingly, we observed an increased microbial diversity in schizophrenia patients compared to the three other groups. We replicated this finding in an independent schizophrenia case-control cohort. This increased diversity is inversely correlated with estimated cell abundance of a subpopulation of CD8+ memory T cells in healthy controls, supporting a link between microbial products found in blood, immunity and schizophrenia.

A gene co-expression network in whole blood of schizophrenia patients is independent of antipsychotic-use and enriched for brain-expressed genes

Despite large-scale genome-wide association studies (GWAS), the underlying genes for schizophrenia are largely unknown. Additional approaches are therefore required to identify the genetic background of this disorder. Here we report findings from a large gene expression study in peripheral blood of schizophrenia patients and controls. We applied a systems biology approach to genome-wide expression data from whole blood of 92 medicated and 29 antipsychotic-free schizophrenia patients and 118 healthy controls. We show that gene expression profiling in whole blood can identify twelve large gene co-expression modules associated with schizophrenia. Several of these disease related modules are likely to reflect expression changes due to antipsychotic medication. However, two of the disease modules could be replicated in an independent second data set involving antipsychotic-free patients and controls. One of these robustly defined disease modules is significantly enriched with brain-expressed genes and with genetic variants that were implicated in a GWAS study, which could imply a causal role in schizophrenia etiology. The most highly connected intramodular hub gene in this module (ABCF1), is located in, and regulated by the major histocompatibility (MHC) complex, which is intriguing in light of the fact that common allelic variants from the MHC region have been implicated in schizophrenia. This suggests that the MHC increases schizophrenia susceptibility via altered gene expression of regulatory genes in this network.