Polygenic Liability to Depression Is Associated With Multiple Medical Conditions in the Electronic Health Record: Phenome-wide Association Study of 46,782 Individuals

Application of polygenic scores to a deeply phenotyped sample enriched for substance use disorders reveals extensive pleiotropy with psychiatric and somatic traits

Co-occurring psychiatric, medical, and substance use disorders (SUDs) are common, but the complex pathways leading to such comorbidities are poorly understood. A greater understanding of genetic influences on this phenomenon could inform precision medicine efforts. We used the Yale-Penn dataset, a cross-sectional sample enriched for individuals with SUDs, to examine pleiotropic effects of genetic liability for psychiatric and somatic traits. Participants completed an in-depth interview that provides information on demographics, environment, medical illnesses, and psychiatric and SUDs. Polygenic scores (PGS) for psychiatric disorders and somatic traits were calculated in European-ancestry (EUR; n = 5691) participants and, when discovery datasets were available, for African-ancestry (AFR; n = 4918) participants. Phenome-wide association studies (PheWAS) were then conducted. In AFR participants, the only PGS with significant associations was bipolar disorder (BD), all of which were with substance use phenotypes. In EUR participants, PGS for major depressive disorder (MDD), generalized anxiety disorder (GAD), post-traumatic stress disorder (PTSD), schizophrenia (SCZ), body mass index (BMI), coronary artery disease (CAD), and type 2 diabetes (T2D) all showed significant associations, the majority of which were with phenotypes in the substance use categories. For instance, PGSMDD was associated with over 200 phenotypes, 15 of which were depression-related (e.g., depression criterion count), 55 of which were other psychiatric phenotypes, and 126 of which were substance use phenotypes; and PGSBMI was associated with 138 phenotypes, 105 of which were substance related. Genetic liability for psychiatric and somatic traits is associated with numerous phenotypes across multiple categories, indicative of the broad genetic liability of these traits.

Identifying data-driven subtypes of major depressive disorder with electronic health records

BACKGROUND

Efforts to reduce the heterogeneity of major depressive disorder (MDD) by identifying subtypes have not yet facilitated treatment personalization or investigation of biology, so novel approaches merit consideration.

METHODS

We utilized electronic health records drawn from 2 academic medical centers and affiliated health systems in Massachusetts to identify data-driven subtypes of MDD, characterizing sociodemographic features, comorbid diagnoses, and treatment patterns. We applied Latent Dirichlet Allocation (LDA) to summarize diagnostic codes followed by agglomerative clustering to define patient subgroups.

RESULTS

Among 136,371 patients (95,034 women [70 %]; 41,337 men [30 %]; mean [SD] age, 47.0 [14.0] years), the 15 putative MDD subtypes were characterized by comorbidities and distinct patterns in medication use. There was substantial variation in rates of selective serotonin reuptake inhibitor (SSRI) use (from a low of 62 % to a high of 78 %) and selective norepinephrine reuptake inhibitor (SNRI) use (from 4 % to 21 %).

LIMITATIONS

Electronic health records lack reliable symptom-level data, so we cannot examine the extent to which subtypes might differ in clinical presentation or symptom dimensions.

CONCLUSION

These data-driven subtypes, drawing on representative clinical cohorts, merit further investigation for their utility in identifying more homogeneous patient populations for basic as well as clinical investigation.

Polygenic risk scores for mood and related disorders and environmental factors: Interaction effects on wellbeing in the UK biobank

Mood disorders have a genetic and environmental component and interactions (GxE) on the risk of psychiatric diseases have been investigated. The same GxE interactions may affect wellbeing measures, which go beyond categorical diagnoses and reflect the health-disease continuum. We evaluated GxE effects in the UK Biobank, considering as outcomes subjective wellbeing (feeling good and functioning well) and objective measures (education and income). We estimated the polygenic risk scores (PRSs) of major depressive disorder, bipolar disorder, schizophrenia, and attention deficit hyperactivity disorder. Stressful/traumatic events during adulthood or childhood were considered as E variables, as well as social support. The addition of the PRSxE interaction to PRS and E variables was tested in linear or multinomial regression models, adjusting for confounders. We included 33 k-380 k participants, depending on the variables considered. Most PRSs and E factors showed additive effects on outcomes, with effect sizes generally 3-5 times larger for E variables than PRSs. We found some interaction effects, particularly when considering recent stress, history of a long illness/disability/infirmity, and social support. Higher PRSs increased the negative effects of stress on wellbeing, but they also increased the positive effects of social support, with interaction effects particularly for the outcomes health satisfaction, loneliness, and income (p < Bonferroni corrected threshold of 1.92e-4). PRSxE terms usually added ∼0.01-0.02% variance explained to the corresponding additive model. PRSxE effects on wellbeing involve both positive and negative E factors. Despite small variance explained at the population level, preventive/therapeutic interventions that modify E factors could be beneficial at the individual level.

Clinical and genetic contributions to medical comorbidity in bipolar disorder: a study using electronic health records-linked biobank data

Bipolar disorder is a chronic and complex polygenic disease with high rates of comorbidity. However, the independent contribution of either diagnosis or genetic risk of bipolar disorder to the medical comorbidity profile of individuals with the disease remains unresolved. Here, we conducted a multi-step phenome-wide association study (PheWAS) of bipolar disorder using phenomes derived from the electronic health records of participants enrolled in the Mayo Clinic Biobank and the Mayo Clinic Bipolar Disorder Biobank. First, we explored the conditions associated with a diagnosis of bipolar disorder by conducting a phenotype-based PheWAS followed by LASSO-penalized regression to account for correlations within the phenome. Then, we explored the conditions associated with bipolar disorder polygenic risk score (BD-PRS) using a PRS-based PheWAS with a sequential exclusion approach to account for the possibility that diagnosis, instead of genetic risk, may drive such associations. 53,386 participants (58.7% women) with a mean age at analysis of 67.8 years (SD = 15.6) were included. A bipolar disorder diagnosis (n = 1479) was associated with higher rates of psychiatric conditions, injuries and poisonings, endocrine/metabolic and neurological conditions, viral hepatitis C, and asthma. BD-PRS was associated with psychiatric comorbidities but, in contrast, had no positive associations with general medical conditions. While our findings warrant confirmation with longitudinal-prospective studies, the limited associations between bipolar disorder genetics and medical conditions suggest that shared environmental effects or environmental consequences of diagnosis may have a greater impact on the general medical comorbidity profile of individuals with bipolar disorder than its genetic risk.

Mendelian randomization study of the relationship between blood and urine biomarkers and schizophrenia in the UK Biobank cohort

BACKGROUND

The identification of suitable biomarkers is of crucial clinical importance for the early diagnosis of treatment-resistant schizophrenia (TRS). This study aims to comprehensively analyze the association between TRS and blood and urine biomarkers.

METHODS

Candidate TRS-related single nucleotide polymorphisms (SNPs) were obtained from a recent genome-wide association study. The UK Biobank cohort, comprising 376,807 subjects with blood and urine biomarker testing data, was used to calculate the polygenic risk score (PRS) for TRS. Pearson correlation analyses were performed to evaluate the correlation between TRS PRS and each of the biomarkers, using calculated TRS PRS as the instrumental variables. Bidirectional two-sample Mendelian randomization (MR) was used to assess potential causal associations between candidate biomarkers with TRS.

RESULTS

Here we identify a significant association between TRS PRS and phosphate (r = 0.007, P = 1.96 × 10-4). Sex subgroup analyses identify seven and three candidate biomarkers associated with TRS PRS in male and female participants, respectively. For example, total protein and phosphate for males, creatinine and phosphate for females. Bidirectional two-sample MR analyses indicate that TRS is negatively associated with cholesterol (estimate = -0.363, P  = 0.008). Conversely, TRS is positively associated with total protein (estimate = 0.137, P  = 0.027), mean corpuscular volume (estimate = 0.032, P  = 2.25 × 10-5), and mean corpuscular hemoglobin (estimate = 0.018, P  = 0.007).

CONCLUSIONS

Our findings provide insights into the roles of blood and urine biomarkers in the early detection and treatment of TRS.

Pangenomics: A new era in the field of neurodegenerative diseases

A pangenome is composed of all the genetic variability of a group of individuals, and its application to the study of neurodegenerative diseases may provide valuable insights into the underlying aspects of genetic heterogenetiy for these complex ailments, including gene expression, epigenetics, and translation mechanisms. Furthermore, a reference pangenome allows for the identification of previously undetected structural commonalities and differences among individuals, which may help in the diagnosis of a disease, support the prediction of what will happen over time (prognosis) and aid in developing novel treatments in the perspective of personalized medicine. Therefore, in the present review, the application of the pangenome concept to the study of neurodegenerative diseases will be discussed and analyzed for its potential to enable an improvement in diagnosis and prognosis for these illnesses, leading to the development of tailored treatments for individual patients from the knowledge of the genomic composition of a whole population.

Application of polygenic scores to a deeply phenotyped sample enriched for substance use disorders reveals extensive pleiotropy with psychiatric and medical traits

Co-occurring psychiatric, medical, and substance use disorders (SUDs) are common, but the complex pathways leading to such comorbidities are poorly understood. A greater understanding of genetic influences on this phenomenon could inform precision medicine efforts. We used the Yale-Penn dataset, a cross-sectional sample enriched for individuals with SUDs, to examine pleiotropic effects of genetic liability for psychiatric and medical traits. Participants completed an in-depth interview that provides information on demographics, environment, medical illnesses, and psychiatric and SUDs. Polygenic scores (PGS) for psychiatric disorders and medical traits were calculated in European-ancestry (EUR; n=5,691) participants and, when discovery datasets were available, for African-ancestry (AFR; n=4,918) participants. Phenome-wide association studies (PheWAS) were then conducted. In AFR participants, the only PGS with significant associations was bipolar disorder (BD), all of which were with substance use phenotypes. In EUR participants, PGS for major depressive disorder (MDD), generalized anxiety disorder (GAD), post-traumatic stress disorder (PTSD), schizophrenia (SCZ), body mass index (BMI), coronary artery disease (CAD), and type 2 diabetes (T2D) all showed significant associations, the majority of which were with phenotypes in the substance use categories. For instance, PGS MDD was associated with over 200 phenotypes, 15 of which were depression-related (e.g., depression criterion count), 55 of which were other psychiatric phenotypes, and 126 of which were substance use phenotypes; and PGS BMI was associated with 138 phenotypes, 105 of which were substance related. Genetic liability for psychiatric and medical traits is associated with numerous phenotypes across multiple categories, indicative of the broad genetic liability of these traits.

Enhanced polygenic risk score incorporating gene-environment interaction suggests the association of major depressive disorder with cardiac and lung function

BACKGROUND

Depression has been linked to an increased risk of cardiovascular and respiratory diseases; however, its impact on cardiac and lung function remains unclear, especially when accounting for potential gene-environment interactions.

METHODS

We developed a novel polygenic and gene-environment interaction risk score (PGIRS) integrating the major genetic effect and gene-environment interaction effect of depression-associated loci. The single nucleotide polymorphisms (SNPs) demonstrating major genetic effect or environmental interaction effect were obtained from genome-wide SNP association and SNP-environment interaction analyses of depression. We then calculated the depression PGIRS for non-depressed individuals, using smoking and alcohol consumption as environmental factors. Using linear regression analysis, we assessed the associations of PGIRS and conventional polygenic risk score (PRS) with lung function (N = 42 886) and cardiac function (N = 1791) in the subjects with or without exposing to smoking and alcohol drinking.

RESULTS

We detected significant associations of depression PGIRS with cardiac and lung function, contrary to conventional depression PRS. Among smokers, forced vital capacity exhibited a negative association with PGIRS (β = -0.037, FDR = 1.00 × 10-8), contrasting with no significant association with PRS (β = -0.002, FDR = 0.943). In drinkers, we observed a positive association between cardiac index with PGIRS (β = 0.088, FDR = 0.010), whereas no such association was found with PRS (β = 0.040, FDR = 0.265). Notably, in individuals who both smoked and drank, forced expiratory volume in 1-second demonstrated a negative association with PGIRS (β = -0.042, FDR = 6.30 × 10-9), but not with PRS (β = -0.003, FDR = 0.857).

CONCLUSIONS

Our findings underscore the profound impact of depression on cardiac and lung function, highlighting the enhanced efficacy of considering gene-environment interactions in PRS-based studies.