Body mass index stratified meta-analysis of genome-wide association studies of polycystic ovary syndrome in women of European ancestry.

BACKGROUND: Polycystic ovary syndrome (PCOS) is a complex multifactorial disorder with a substantial genetic component. However, the clinical manifestations of PCOS are heterogeneous with notable differences between lean and obese women, implying a different pathophysiology manifesting in differential body mass index (BMI). We performed a meta-analysis of genome-wide association study (GWAS) data from six well-characterised cohorts, using a case-control study design stratified by BMI, aiming to identify genetic variants associated with lean and overweight/obese PCOS subtypes. RESULTS: The study comprised 254,588 women (5,937 cases and 248,651 controls) from individual studies performed in Australia, Estonia, Finland, the Netherlands and United States of America, and separated according to three BMI stratifications (lean, overweight and obese). Genome-wide association analyses were performed for each stratification within each cohort, with the data for each BMI group meta-analysed using METAL software. Almost half of the total study population (47%, n = 119,584) were of lean BMI (≤ 25 kg/m2). Two genome-wide significant loci were identified for lean PCOS, led by rs12000707 within DENND1A (P = 1.55 × 10-12) and rs2228260 within XBP1 (P = 3.68 × 10-8). One additional locus, LINC02905, was highlighted as significantly associated with lean PCOS through gene-based analyses (P = 1.76 × 10-6). There were no significant loci observed for the overweight or obese sub-strata when analysed separately, however, when these strata were combined, an association signal led by rs569675099 within DENND1A reached genome-wide significance (P = 3.22 × 10-9) and a gene-based association was identified with ERBB4 (P = 1.59 × 10-6). Nineteen of 28 signals identified in previous GWAS, were replicated with consistent allelic effect in the lean stratum. There were less replicated signals in the overweight and obese groups, and only 4 SNPs were replicated in each of the three BMI strata. CONCLUSIONS: Genetic variation at the XBP1, LINC02905 and ERBB4 loci were associated with PCOS within unique BMI strata, while DENND1A demonstrated associations across multiple strata, providing evidence of both distinct and shared genetic features between lean and overweight/obese PCOS-affected women. This study demonstrated that PCOS-affected women with contrasting body weight are not only phenotypically distinct but also show variation in genetic architecture; lean PCOS women typically display elevated gonadotrophin ratios, lower insulin resistance, higher androgen levels, including adrenal androgens, and more favourable lipid profiles. Overall, these findings add to the growing body of evidence supporting a genetic basis for PCOS as well as differences in genetic patterns relevant to PCOS BMI-subtype.

Phenomewide Association Study of Health Outcomes Associated With the Genetic Correlates of 25 Hydroxyvitamin D Concentration and Vitamin D Binding Protein Concentration.

While it is known that vitamin D deficiency is associated with adverse bone outcomes, it remains unclear whether low vitamin D status may increase the risk of a wider range of health outcomes. We had the opportunity to explore the association between common genetic variants associated with both 25 hydroxyvitamin D (25OHD) and the vitamin D binding protein (DBP, encoded by the GC gene) with a comprehensive range of health disorders and laboratory tests in a large academic medical center. We used summary statistics for 25OHD and DBP to generate polygenic scores (PGS) for 66,482 participants with primarily European ancestry and 13,285 participants with primarily African ancestry from the Vanderbilt University Medical Center Biobank (BioVU). We examined the predictive properties of PGS25OHD, and two scores related to DBP concentration with respect to 1322 health-related phenotypes and 315 laboratory-measured phenotypes from electronic health records. In those with European ancestry: (a) the PGS25OHD and PGSDBP scores, and individual SNPs rs4588 and rs7041 were associated with both 25OHD concentration and 1,25 dihydroxyvitamin D concentrations; (b) higher PGS25OHD was associated with decreased concentrations of triglycerides and cholesterol, and reduced risks of vitamin D deficiency, disorders of lipid metabolism, and diabetes. In general, the findings for the African ancestry group were consistent with findings from the European ancestry analyses. Our study confirms the utility of PGS and two key variants within the GC gene (rs4588 and rs7041) to predict the risk of vitamin D deficiency in clinical settings and highlights the shared biology between vitamin D-related genetic pathways a range of health outcomes.

Medical and genetic correlates of long-term buprenorphine treatment in the electronic health records.

Despite the benefits associated with longer buprenorphine treatment duration (i.e., >180 days) (BTD) for opioid use disorder (OUD), retention remains poor. Research on the impact of co-occurring psychiatric issues on BTD has yielded mixed results. It is also unknown whether the genetic risk in the form of polygenic scores (PGS) for OUD and other comorbid conditions, including problematic alcohol use (PAU) are associated with BTD. We tested the association between somatic and psychiatric comorbidities and long BTD and determined whether PGS for OUD-related conditions was associated with BTD. The study included 6686 individuals with a buprenorphine prescription that lasted for less than 6 months (short-BTD) and 1282 individuals with a buprenorphine prescription that lasted for at least 6 months (long-BTD). Recorded diagnosis of substance addiction and disorders (Odds Ratio (95% CI) = 22.14 (21.88-22.41), P = 2.8 × 10-116), tobacco use disorder (OR (95% CI) = 23.4 (23.13-23.68), P = 4.5 × 10-111), and bipolar disorder (OR(95% CI) = 9.70 (9.48-9.92), P = 1.3 × 10-91), among others, were associated with longer BTD. The PGS of OUD and several OUD co-morbid conditions were associated with any buprenorphine prescription. A higher PGS for OUD (OR per SD increase in PGS (95%CI) = 1.43(1.16-1.77), P = 0.0009), loneliness (OR(95% CI) = 1.39(1.13-1.72), P = 0.002), problematic alcohol use (OR(95%CI) = 1.47(1.19-1.83), P = 0.0004), and externalizing disorders (OR(95%CI) = 1.52(1.23 to 1.89), P = 0.0001) was significantly associated with long-BTD. Associations between BTD and the PGS of depression, chronic pain, nicotine dependence, cannabis use disorder, and bipolar disorder did not survive correction for multiple testing. Longer BTD is associated with diagnoses of psychiatric and somatic conditions in the EHR, as is the genetic score for OUD, loneliness, problematic alcohol use, and externalizing disorders.

Developing a phenotype risk score for tic disorders in a large, clinical biobank.

Tics are a common feature of early-onset neurodevelopmental disorders, characterized by involuntary and repetitive movements or sounds. Despite affecting up to 2% of children and having a genetic contribution, the underlying causes remain poorly understood. In this study, we leverage dense phenotype information to identify features (i.e., symptoms and comorbid diagnoses) of tic disorders within the context of a clinical biobank. Using de-identified electronic health records (EHRs), we identified individuals with tic disorder diagnosis codes. We performed a phenome-wide association study (PheWAS) to identify the EHR features enriched in tic cases versus controls (n = 1406 and 7030; respectively) and found highly comorbid neuropsychiatric phenotypes, including: obsessive-compulsive disorder, attention-deficit/hyperactivity disorder, autism spectrum disorder, and anxiety (p < 7.396 × 10-5). These features (among others) were then used to generate a phenotype risk score (PheRS) for tic disorder, which was applied across an independent set of 90,051 individuals. A gold standard set of tic disorder cases identified by an EHR algorithm and confirmed by clinician chart review was then used to validate the tic disorder PheRS; the tic disorder PheRS was significantly higher among clinician-validated tic cases versus non-cases (p = 4.787 × 10-151; ß = 1.68; SE = 0.06). Our findings provide support for the use of large-scale medical databases to better understand phenotypically complex and underdiagnosed conditions, such as tic disorders.

Cross-ancestry genetic investigation of schizophrenia, cannabis use disorder, and tobacco smoking.

Individuals with schizophrenia frequently experience co-occurring substance use, including tobacco smoking and heavy cannabis use, and substance use disorders. There is interest in understanding the extent to which these relationships are causal, and to what extent shared genetic factors play a role. We explored the relationships between schizophrenia (Scz), cannabis use disorder (CanUD), and ever-regular tobacco smoking (Smk) using the largest available genome-wide studies of these phenotypes in individuals of African and European ancestries. All three phenotypes were positively genetically correlated (rgs = 0.17 - 0.62). Causal inference analyses suggested the presence of horizontal pleiotropy, but evidence for bidirectional causal relationships was also found between all three phenotypes even after correcting for horizontal pleiotropy. We identified 439 pleiotropic loci in the European ancestry data, 150 of which were novel (i.e., not genome-wide significant in the original studies). Of these pleiotropic loci, 202 had lead variants which showed convergent effects (i.e., same direction of effect) on Scz, CanUD, and Smk. Genetic variants convergent across all three phenotypes showed strong genetic correlations with risk-taking, executive function, and several mental health conditions. Our results suggest that both horizontal pleiotropy and causal mechanisms may play a role in the relationship between CanUD, Smk, and Scz, but longitudinal, prospective studies are needed to confirm a causal relationship.

Genome-wide association studies of lifetime and frequency cannabis use in 131,895 individuals.

Cannabis is one of the most widely used drugs globally. Decriminalization of cannabis is further increasing cannabis consumption. We performed genome-wide association studies (GWASs) of lifetime (N=131,895) and frequency (N=73,374) of cannabis use. Lifetime cannabis use GWAS identified two loci, one near CADM2 (rs11922956, p=2.40E-11) and another near GRM3 (rs12673181, p=6.90E-09). Frequency of use GWAS identified one locus near CADM2 (rs4856591, p=8.10E-09; r2 =0.76 with rs11922956). Both traits were heritable and genetically correlated with previous GWASs of lifetime use and cannabis use disorder (CUD), as well as other substance use and cognitive traits. Polygenic scores (PGSs) for lifetime and frequency of cannabis use associated cannabis use phenotypes in AllofUs participants. Phenome-wide association study of lifetime cannabis use PGS in a hospital cohort replicated associations with substance use and mood disorders, and uncovered associations with celiac and infectious diseases. This work demonstrates the value of GWASs of CUD transition risk factors.

Development and multi-site external validation of a generalizable risk prediction model for bipolar disorder.

Bipolar disorder is a leading contributor to disability, premature mortality, and suicide. Early identification of risk for bipolar disorder using generalizable predictive models trained on diverse cohorts around the United States could improve targeted assessment of high risk individuals, reduce misdiagnosis, and improve the allocation of limited mental health resources. This observational case-control study intended to develop and validate generalizable predictive models of bipolar disorder as part of the multisite, multinational PsycheMERGE Network across diverse and large biobanks with linked electronic health records (EHRs) from three academic medical centers: in the Northeast (Massachusetts General Brigham), the Mid-Atlantic (Geisinger) and the Mid-South (Vanderbilt University Medical Center). Predictive models were developed and valid with multiple algorithms at each study site: random forests, gradient boosting machines, penalized regression, including stacked ensemble learning algorithms combining them. Predictors were limited to widely available EHR-based features agnostic to a common data model including demographics, diagnostic codes, and medications. The main study outcome was bipolar disorder diagnosis as defined by the International Cohort Collection for Bipolar Disorder, 2015. In total, the study included records for 3,529,569 patients including 12,533 cases (0.3%) of bipolar disorder. After internal and external validation, algorithms demonstrated optimal performance in their respective development sites. The stacked ensemble achieved the best combination of overall discrimination (AUC = 0.82-0.87) and calibration performance with positive predictive values above 5% in the highest risk quantiles at all three study sites. In conclusion, generalizable predictive models of risk for bipolar disorder can be feasibly developed across diverse sites to enable precision medicine. Comparison of a range of machine learning methods indicated that an ensemble approach provides the best performance overall but required local retraining. These models will be disseminated via the PsycheMERGE Network website.

Comorbidity Profiles of Posttraumatic Stress Disorder Across the Medical Phenome.

Background: Previous epidemiological research has linked posttraumatic stress disorder (PTSD) with specific physical health problems, but the comprehensive landscape of medical conditions associated with PTSD remains uncharacterized. Electronic health records provide an opportunity to overcome clinical knowledge gaps and uncover associations with biological relevance that potentially vary by sex. Methods: PTSD was defined among biobank participants (N = 145,959) in 3 major healthcare systems using 2 ICD code-based definitions: broad (≥1 PTSD or acute stress codes vs. 0; n cases = 16,706) and narrow (≥2 PTSD codes vs. 0; n cases = 3325). Using a phenome-wide association study design, we tested associations between each PTSD definition and all prevalent disease umbrella categories, i.e., phecodes. We also conducted sex-stratified phenome-wide association study analyses including a sex × diagnosis interaction term in each logistic regression. Results: A substantial number of phecodes were significantly associated with PTSDNarrow (61%) and PTSDBroad (83%). While the strongest associations were shared between the 2 definitions, PTSDBroad captured 334 additional phecodes not significantly associated with PTSDNarrow and exhibited a wider range of significantly associated phecodes across various categories, including respiratory, genitourinary, and circulatory conditions. Sex differences were observed in that PTSDBroad was more strongly associated with osteoporosis, respiratory failure, hemorrhage, and pulmonary heart disease among male patients and with urinary tract infection, acute pharyngitis, respiratory infections, and overweight among female patients. Conclusions: This study provides valuable insights into a diverse range of comorbidities associated with PTSD, including both known and novel associations, while highlighting the influence of sex differences and the impact of defining PTSD using electronic health records.

Posttraumatic stress disorder (PTSD) is a debilitating psychiatric disorder that some people develop following a traumatic event. In addition to mental symptoms, PTSD can impact human health in ways other ways; for example, there are many known conditions that co-occur with PTSD, such as cardiovascular conditions. In this study, we set out to understand the breadth and degree to which PTSD co-occurs with medical outcomes in a sample of over 146,000 patients across 3 large medical systems. We found that both narrowly and broadly defined PTSD diagnosis co-occurred with hundreds of medical conditions, and the strongest associations were with other psychiatric disorders, respiratory conditions (asthma, GERD), sleep-related conditions, and pain. These results provide insights into future genetic studies of PTSD in large-scale biobanks and deepen our understanding of the complex needs of patients with PTSD.

Cross-ancestry genetic investigation of schizophrenia, cannabis use disorder, and tobacco smoking.

Individuals with schizophrenia frequently experience co-occurring substance use, including tobacco smoking and heavy cannabis use, and substance use disorders. There is interest in understanding the extent to which these relationships are causal, and to what extent shared genetic factors play a role. We explored the relationships between schizophrenia (Scz; European ancestry N = 161,405; African ancestry N = 15,846), cannabis use disorder (CanUD; European ancestry N = 886,025; African ancestry N = 120,208), and ever-regular tobacco smoking (Smk; European ancestry N = 805,431; African ancestry N = 24,278) using the largest available genome-wide studies of these phenotypes in individuals of African and European ancestries. All three phenotypes were positively genetically correlated (rgs = 0.17-0.62). Genetic instrumental variable analyses suggested the presence of shared heritable factors, but evidence for bidirectional causal relationships was also found between all three phenotypes even after correcting for these shared genetic factors. We identified 327 pleiotropic loci with 439 lead SNPs in the European ancestry data, 150 of which were novel (i.e., not genome-wide significant in the original studies). Of these pleiotropic loci, 202 had lead variants which showed convergent effects (i.e., same direction of effect) on Scz, CanUD, and Smk. Genetic variants convergent across all three phenotypes showed strong genetic correlations with risk-taking, executive function, and several mental health conditions. Our results suggest that both shared genetic factors and causal mechanisms may play a role in the relationship between CanUD, Smk, and Scz, but longitudinal, prospective studies are needed to confirm a causal relationship.

Genome-wide association studies of coffee intake in UK/US participants of European ancestry uncover cohort-specific genetic associations.

Coffee is one of the most widely consumed beverages. We performed a genome-wide association study (GWAS) of coffee intake in US-based 23andMe participants (N = 130,153) and identified 7 significant loci, with many replicating in three multi-ancestral cohorts. We examined genetic correlations and performed a phenome-wide association study across hundreds of biomarkers, health, and lifestyle traits, then compared our results to the largest available GWAS of coffee intake from the UK Biobank (UKB; N = 334,659). We observed consistent positive genetic correlations with substance use and obesity in both cohorts. Other genetic correlations were discrepant, including positive genetic correlations between coffee intake and psychiatric illnesses, pain, and gastrointestinal traits in 23andMe that were absent or negative in the UKB, and genetic correlations with cognition that were negative in 23andMe but positive in the UKB. Phenome-wide association study using polygenic scores of coffee intake derived from 23andMe or UKB summary statistics also revealed consistent associations with increased odds of obesity- and red blood cell-related traits, but all other associations were cohort-specific. Our study shows that the genetics of coffee intake associate with substance use and obesity across cohorts, but also that GWAS performed in different populations could capture cultural differences in the relationship between behavior and genetics.

Multi-ancestry meta-analysis of tobacco use disorder identifies 461 potential risk genes and reveals associations with multiple health outcomes.

Tobacco use disorder (TUD) is the most prevalent substance use disorder in the world. Genetic factors influence smoking behaviours and although strides have been made using genome-wide association studies to identify risk variants, most variants identified have been for nicotine consumption, rather than TUD. Here we leveraged four US biobanks to perform a multi-ancestral meta-analysis of TUD (derived via electronic health records) in 653,790 individuals (495,005 European, 114,420 African American and 44,365 Latin American) and data from UK Biobank (ncombined = 898,680). We identified 88 independent risk loci; integration with functional genomic tools uncovered 461 potential risk genes, primarily expressed in the brain. TUD was genetically correlated with smoking and psychiatric traits from traditionally ascertained cohorts, externalizing behaviours in children and hundreds of medical outcomes, including HIV infection, heart disease and pain. This work furthers our biological understanding of TUD and establishes electronic health records as a source of phenotypic information for studying the genetics of TUD.

Multi-ancestry genome-wide association study of major depression aids locus discovery, fine mapping, gene prioritization and causal inference.

Most genome-wide association studies (GWAS) of major depression (MD) have been conducted in samples of European ancestry. Here we report a multi-ancestry GWAS of MD, adding data from 21 cohorts with 88,316 MD cases and 902,757 controls to previously reported data. This analysis used a range of measures to define MD and included samples of African (36% of effective sample size), East Asian (26%) and South Asian (6%) ancestry and Hispanic/Latin American participants (32%). The multi-ancestry GWAS identified 53 significantly associated novel loci. For loci from GWAS in European ancestry samples, fewer than expected were transferable to other ancestry groups. Fine mapping benefited from additional sample diversity. A transcriptome-wide association study identified 205 significantly associated novel genes. These findings suggest that, for MD, increasing ancestral and global diversity in genetic studies may be particularly important to ensure discovery of core genes and inform about transferability of findings.