Unraveling COVID-19 relationship with anxiety disorders and symptoms using genome-wide data

BACKGROUND

There is still a limited understanding of the dynamics contributing to the comorbidity of COVID-19 and anxiety outcomes.

METHODS

To dissect the pleiotropic mechanisms contributing to COVID-19/anxiety comorbidity, we used genome-wide data from UK Biobank (up to 420,531 participants), FinnGen Project (up to 329,077 participants), Million Veteran Program (175,163 participants), and COVID-19 Host Genetics Initiative (up to 122,616 cases and 2,475,240 controls). Specifically, we assessed global and local genetic correlation and genetically inferred effects linking COVID-19 outcomes (infection, hospitalization, and severe respiratory symptoms) to anxiety disorders and symptoms.

RESULTS

We observed a strong genetic correlation of anxiety disorder with COVID-19 positive status (rg = 0.35, p = 2×10-4) and COVID-19 hospitalization (rg = 0.31, p = 7.2×10-4). Among anxiety symptoms, "Tense, sore, or aching muscles during worst period of anxiety" was genetically correlated with COVID-19 positive status (rg = 0.33, p = 0.001), while "Frequent trouble falling or staying asleep during worst period of anxiety" was genetically correlated with COVID-19 hospitalization (rg = 0.24, p = 0.004). Through a latent causal variable analysis, we observed that COVID-19 outcomes have statistically significant genetic causality proportion (gcp) on anxiety symptoms (e.g., COVID-19 positive status→"Recent easy annoyance or irritability" │gcp│ = 0.18, p = 6.72×10-17). Conversely, anxiety disorders appear to have a possible causal effect on COVID-19 (│gcp│ = 0.38, p = 3.17×10-9). Additionally, we also identified multiple loci with evidence of local genetic correlation between anxiety and COVID-19. These appear to be related to genetic effects shared with lung function, brain morphology, alcohol and tobacco use, and hematologic parameters.

CONCLUSIONS

This study provided insights into the pleiotropic mechanisms linking COVID-19 and anxiety outcomes, suggesting differences between dynamics related to anxiety disorders and those related to anxiety symptoms.

Genome-wide association analyses identify 95 risk loci and provide insights into the neurobiology of post-traumatic stress disorder

Post-traumatic stress disorder (PTSD) genetics are characterized by lower discoverability than most other psychiatric disorders. The contribution to biological understanding from previous genetic studies has thus been limited. We performed a multi-ancestry meta-analysis of genome-wide association studies across 1,222,882 individuals of European ancestry (137,136 cases) and 58,051 admixed individuals with African and Native American ancestry (13,624 cases). We identified 95 genome-wide significant loci (80 new). Convergent multi-omic approaches identified 43 potential causal genes, broadly classified as neurotransmitter and ion channel synaptic modulators (for example, GRIA1, GRM8 and CACNA1E), developmental, axon guidance and transcription factors (for example, FOXP2, EFNA5 and DCC), synaptic structure and function genes (for example, PCLO, NCAM1 and PDE4B) and endocrine or immune regulators (for example, ESR1, TRAF3 and TANK). Additional top genes influence stress, immune, fear and threat-related processes, previously hypothesized to underlie PTSD neurobiology. These findings strengthen our understanding of neurobiological systems relevant to PTSD pathophysiology, while also opening new areas for investigation.

Longitudinal Patterns of Multimorbidity in Gulf War Era Veterans With and Without Gulf War Illness

Objectives: To examine whether severe Gulf War illness (SGWI) case status was associated with longitudinal multimorbidity patterns. Methods: Participants were users of the Veteran Health Administration Health Care System drawn from the Gulf War Era Cohort and Biorepository (n = 840). Longitudinal measures of multimorbidity were constructed using (1) electronic health records (Charlson Comorbidity Index; Elixhauser; and Veterans Affairs Frailty Index) from 10/1/1999 to 6/30/2023 and (2) self-reported medical conditions (Deficit Accumulation Index) since the war until the survey date. Accelerated failure time models examined SGWI case status as a predictor of time until threshold level of multimorbidity was reached, adjusted for age and sociodemographic and military characteristics. Results: Models, adjusted for covariates, revealed that (1) relative to the SWGI- group, the SGWI+ group was associated with an accelerated time for reaching each threshold and (2) the relationship between SGWI and each threshold was not moderated by age. Discussion: Findings suggest that veterans with SGWI experienced accelerated aging.

Clinical spectrum of Transthyretin amyloidogenic mutations among diverse population origins

PURPOSE

Coding mutations in the Transthyretin (TTR) gene cause a hereditary form of amyloidosis characterized by a complex genotype-phenotype correlation with limited information regarding differences among worldwide populations.

METHODS

We compared 676 diverse individuals carrying TTR amyloidogenic mutations (rs138065384, Phe44Leu; rs730881165, Ala81Thr; rs121918074, His90Asn; rs76992529, Val122Ile) to 12,430 non-carriers matched by age, sex, and genetically-inferred ancestry to assess their clinical presentations across 1,693 outcomes derived from electronic health records in UK biobank.

RESULTS

In individuals of African descent (AFR), Val122Ile mutation was linked to multiple outcomes related to the circulatory system (fold-enrichment = 2.96, p = 0.002) with the strongest associations being cardiac congenital anomalies (phecode 747.1, p = 0.003), endocarditis (phecode 420.3, p = 0.006), and cardiomyopathy (phecode 425, p = 0.007). In individuals of Central-South Asian descent (CSA), His90Asn mutation was associated with dermatologic outcomes (fold-enrichment = 28, p = 0.001). The same TTR mutation was linked to neoplasms in European-descent individuals (EUR, fold-enrichment = 3.09, p = 0.003). In EUR, Ala81Thr showed multiple associations with respiratory outcomes related (fold-enrichment = 3.61, p = 0.002), but the strongest association was with atrioventricular block (phecode 426.2, p = 2.81 × 10- 4). Additionally, the same mutation in East Asians (EAS) showed associations with endocrine-metabolic traits (fold-enrichment = 4.47, p = 0.003). In the cross-ancestry meta-analysis, Val122Ile mutation was associated with peripheral nerve disorders (phecode 351, p = 0.004) in addition to cardiac congenital anomalies (fold-enrichment = 6.94, p = 0.003).

CONCLUSIONS

Overall, these findings highlight that TTR amyloidogenic mutations present ancestry-specific and ancestry-convergent associations related to a range of health domains. This supports the need to increase awareness regarding the range of outcomes associated with TTR mutations across worldwide populations to reduce misdiagnosis and delayed diagnosis of TTR-related amyloidosis.

A genetic exploration of the relationship between Posttraumatic Stress Disorder and cardiovascular diseases

Background and Aims

Experiencing a traumatic event may lead to Posttraumatic Stress Disorder (PTSD), including symptoms such as flashbacks and hyperarousal. Individuals suffering from PTSD are at increased risk of cardiovascular disease (CVD), but it is unclear why. This study assesses shared genetic liability and potential causal pathways between PTSD and CVD.

Methods

We leveraged summary-level data of genome-wide association studies (PTSD: N= 1,222,882; atrial fibrillation (AF): N=482,409; coronary artery disease (CAD): N=1,165,690; hypertension: N=458,554; heart failure (HF): N=977,323). First, we estimated genetic correlations and utilized genomic structural equation modeling to identify a common genetic factor for PTSD and CVD. Next, we assessed biological, behavioural, and psychosocial factors as potential mediators. Finally, we employed multivariable Mendelian randomization to examine causal pathways between PTSD and CVD, incorporating the same potential mediators.

Results

Significant genetic correlations were found between PTSD and CAD, HT, and HF (rg =0.21-0.32, p≤ 3.08 · 10-16), but not between PTSD and AF. Insomnia, smoking, alcohol dependence, waist-to-hip ratio, and inflammation (IL6, C-reactive protein) partly mediated these associations. Mendelian randomization indicated that PTSD causally increases CAD (IVW OR=1.53, 95% CIs=1.19-1.96, p=0.001), HF (OR=1.44, CIs=1.08-1.92, p=0.012), and to a lesser degree hypertension (OR=1.25, CIs=1.05-1.49, p=0.012). While insomnia, smoking, alcohol, and inflammation were important mediators, independent causal effects also remained.

Conclusions

In addition to shared genetic liability between PTSD and CVD, we present strong evidence for causal effects of PTSD on CVD. Crucially, we implicate specific lifestyle and biological mediators (insomnia, substance use, inflammation) which has important implications for interventions to prevent CVD in PTSD patients.

Sex differences in the pleiotropy of hearing difficulty with imaging-derived phenotypes: a brain-wide investigation

Hearing difficulty (HD) is one of the major health burdens in older adults. While aging-related changes in the peripheral auditory system play an important role, genetic variation associated with brain structure and function could also be involved in HD predisposition. We analyzed a large-scale HD genome-wide association study (GWAS; Ntotal = 501,825, 56% females) and GWAS data related to 3,935 brain imaging-derived phenotypes (IDPs) assessed in up to 33,224 individuals (52% females) using multiple magnetic resonance imaging modalities. To investigate HD pleiotropy with brain structure and function, we conducted genetic correlation, latent causal variable, Mendelian randomization, and multivariable generalized linear regression analyses. Additionally, we performed local genetic correlation and multi-trait colocalization analyses to identify genomic regions and loci implicated in the pleiotropic mechanisms shared between HD and brain IDPs. We observed a widespread genetic correlation of HD with 120 IDPs in females, 89 IDPs in males, and 171 IDPs in the sex-combined analysis. The latent causal variable analysis showed that some of these genetic correlations could be due to cause-effect relationships. For seven correlations, the causal effects were also confirmed by the Mendelian randomization approach: vessel volume→HD in the sex-combined analysis; hippocampus volume→HD, cerebellum grey matter volume→HD, primary visual cortex volume→HD, and HD→fluctuation amplitudes of node 46 in resting-state functional MRI dimensionality 100 in females; global mean thickness→HD and HD→mean orientation dispersion index in superior corona radiata in males. The local genetic correlation analysis identified 13 pleiotropic regions between HD and these seven IDPs. We also observed a colocalization signal for the rs13026575 variant between HD, primary visual cortex volume, and SPTBN1 transcriptomic regulation in females. Brain structure and function may have a role in the sex differences in HD predisposition via possible cause-effect relationships and shared regulatory mechanisms.

Genetic contribution to the comorbidity between attention-deficit/hyperactivity disorder and substance use disorders

We characterized the genetic architecture of the attention-deficit hyperactivity disorder-substance use disorder (ADHD-SUD) relationship by investigating genetic correlation, causality, pleiotropy, and common polygenic risk. Summary statistics from genome-wide association studies (GWAS) were used to investigate ADHD (Neff = 51,568), cannabis use disorder (CanUD, Neff = 161,053), opioid use disorder (OUD, Neff = 57,120), problematic alcohol use (PAU, Neff = 502,272), and problematic tobacco use (PTU, Neff = 97,836). ADHD, CanUD, and OUD GWAS meta-analyses included cohorts with case definitions based on different diagnostic criteria. PAU GWAS combined information related to alcohol use disorder, alcohol dependence, and the items related to alcohol problematic consequences assessed by the alcohol use disorders identification test. PTU GWAS was generated a multi-trait analysis including information regarding Fagerström Test for Nicotine Dependence and cigarettes per day. Linkage disequilibrium score regression analyses indicated positive genetic correlation with CanUD, OUD, PAU, and PTU. Genomic structural equation modeling showed that these genetic correlations were related to two latent factors: one including ADHD, CanUD, and PTU and the other with OUD and PAU. The evidence of a causal effect of PAU and PTU on ADHD was stronger than the reverse in the two-sample Mendelian randomization analysis. Conversely, similar strength of evidence was found between ADHD and CanUD. CADM2 rs62250713 was a pleiotropic SNP between ADHD and all SUDs. We found seven, one, and twenty-eight pleiotropic variants between ADHD and CanUD, PAU, and PTU, respectively. Finally, OUD, CanUD, and PAU PRS were associated with increased odds of ADHD. Our findings demonstrated the contribution of multiple pleiotropic mechanisms to the comorbidity between ADHD and SUDs.

The impact of assortative mating, participation bias and socioeconomic status on the polygenic risk of behavioural and psychiatric traits

To investigate assortative mating (AM), participation bias and socioeconomic status (SES) with respect to the genetics of behavioural and psychiatric traits, we estimated AM signatures using gametic phase disequilibrium and within-spouses and within-siblings polygenic risk score correlation analyses, also performing a SES conditional analysis. The cross-method meta-analysis identified AM genetic signatures for multiple alcohol-related phenotypes, bipolar disorder, major depressive disorder, schizophrenia and Tourette syndrome. Here, after SES conditioning, we observed changes in the AM genetic signatures for maximum habitual alcohol intake, frequency of drinking alcohol and Tourette syndrome. We also observed significant gametic phase disequilibrium differences between UK Biobank mental health questionnaire responders versus non-responders for major depressive disorder and alcohol use disorder. These results highlight the impact of AM, participation bias and SES on the polygenic risk of behavioural and psychiatric traits, particularly in alcohol-related traits.

Gene Discovery and Biological Insights into Anxiety Disorders from a Multi-Ancestry Genome-wide Association Study of >1.2 Million Participants

We leveraged information from more than 1.2 million participants to investigate the genetics of anxiety disorders across five continental ancestral groups. Ancestry-specific and cross-ancestry genome-wide association studies identified 51 anxiety-associated loci, 39 of which are novel. Additionally, polygenic risk scores derived from individuals of European descent were associated with anxiety in African, Admixed-American, and East Asian groups. The heritability of anxiety was enriched for genes expressed in the limbic system, the cerebral cortex, the cerebellum, the metencephalon, the entorhinal cortex, and the brain stem. Transcriptome- and proteome-wide analyses highlighted 115 genes associated with anxiety through brain-specific and cross-tissue regulation. We also observed global and local genetic correlations with depression, schizophrenia, and bipolar disorder and putative causal relationships with several physical health conditions. Overall, this study expands the knowledge regarding the genetic risk and pathogenesis of anxiety disorders, highlighting the importance of investigating diverse populations and integrating multi-omics information.

Pleiotropy and genetically inferred causality linking multisite chronic pain to substance use disorders

Individuals suffering from chronic pain develop substance use disorders (SUDs) more often than others. Understanding the shared genetic influences underlying the comorbidity between chronic pain and SUDs will lead to a greater understanding of their biology. Genome-wide association statistics were obtained from the UK Biobank for multisite chronic pain (MCP, Neffective = 387,649) and from the Million Veteran Program and the Psychiatric Genomics Consortium meta-analyses for alcohol use disorder (AUD, Neffective = 296,974), cannabis use disorder (CanUD, Neffective = 161,053), opioid use disorder (OUD, Neffective = 57,120), and problematic tobacco use (PTU, Neffective = 270,120). SNP-based heritability was estimated for each of the traits and genetic correlation (rg) analyses were performed to assess MCP-SUD pleiotropy. Bidirectional Mendelian Randomization analyses evaluated possible causal relationships. Finally, to identify and characterize individual loci, we performed a genome-wide pleiotropy analysis and a brain-wide analysis using imaging phenotypes available from the UK Biobank. MCP was positively genetically correlated with AUD (rg = 0.26, p = 7.55 × 10-18), CanUD (rg = 0.37, p = 8.21 × 10-37), OUD (rg = 0.20, p = 1.50 × 10-3), and PTU (rg = 0.29, p = 8.53 × 10-12). Although the MR analyses supported bi-directional relationships, MCP had larger effects on AUD (pain-exposure: beta = 0.18, p = 8.21 × 10-4; pain-outcome: beta = 0.07, p = 0.018), CanUD (pain-exposure: beta = 0.58, p = 2.70 × 10-6; pain-outcome: beta = 0.05, p = 0.014) and PTU (pain-exposure: beta = 0.43, p = 4.16 × 10-8; pain-outcome: beta = 0.09, p = 3.05 × 10-6) than the reverse. The genome-wide analysis identified two SNPs pleiotropic between MCP and all SUD investigated: IHO1 rs7652746 (ppleiotropy = 2.69 × 10-8), and CADM2 rs1248857 (ppleiotropy = 1.98 × 10-5). In the brain-wide analysis, rs7652746 was associated with multiple cerebellum and amygdala imaging phenotypes. When analyzing MCP pleiotropy with each SUD separately, we found 25, 22, and 4 pleiotropic variants for AUD, CanUD, and OUD, respectively. To our knowledge, this is the first large-scale study to provide evidence of potential causal relationships and shared genetic mechanisms underlying MCP-SUD comorbidity.

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.

Genetic Liability to Posttraumatic Stress Disorder Symptoms and Its Association With Cardiometabolic and Respiratory Outcomes

Importance

Posttraumatic stress disorder (PTSD) has been reported to be a risk factor for several physical and somatic symptoms. However, the genetics of PTSD and its potential association with medical outcomes remain unclear.

Objective

To examine disease categories and laboratory tests from electronic health records (EHRs) that are associated with PTSD polygenic scores.

Design, Setting, and Participants

This genetic association study was conducted from July 15, 2021, to January 24, 2023, using EHR data from participants across 4 biobanks. The polygenic scores of PTSD symptom severity (PGS-PTSD) were tested with all available phecodes in Vanderbilt University Medical Center's biobank (BioVU), Mass General Brigham (MGB), Michigan Genomics Initiative (MGI), and UK Biobank (UKBB). The significant medical outcomes were tested for overrepresented disease categories and subsequently tested for genetic correlation and 2-sample mendelian randomization (MR) to determine genetically informed associations. Multivariable MR was conducted to assess whether PTSD associations with health outcomes were independent of the genetic effect of body mass index and tobacco smoking.

Exposures

Polygenic score of PTSD symptom severity.

Main Outcomes and Measures

A total of 1680 phecodes (ie, International Classification of Diseases, Ninth Revision- and Tenth Revision-based phenotypic definitions of health outcomes) across 4 biobanks and 490 laboratory tests across 2 biobanks (BioVU and MGB).

Results

In this study including a total of 496 317 individuals (mean [SD] age, 56.8 [8.0] years; 263 048 female [53%]) across the 4 EHR sites, meta-analyzing associations of PGS-PTSD with 1680 phecodes from 496 317 individuals showed significant associations to be overrepresented from mental health disorders (fold enrichment = 3.15; P = 5.81 × 10-6), circulatory system (fold enrichment = 3.32; P = 6.39 × 10-12), digestive (fold enrichment = 2.42; P = 2.16 × 10-7), and respiratory outcomes (fold enrichment = 2.51; P = 8.28 × 10-5). The laboratory measures scan with PGS-PTSD in BioVU and MGB biobanks revealed top associations in metabolic and immune domains. MR identified genetic liability to PTSD symptom severity as an associated risk factor for 12 health outcomes, including alcoholism (β = 0.023; P = 1.49 × 10-4), tachycardia (β = 0.045; P = 8.30 × 10-5), cardiac dysrhythmias (β = 0.016, P = 3.09 × 10-5), and acute pancreatitis (β = 0.049, P = 4.48 × 10-4). Several of these associations were robust to genetic effects of body mass index and smoking. We observed a bidirectional association between PTSD symptoms and nonspecific chest pain and C-reactive protein.

Conclusions and Relevance

Results of this study suggest the broad health repercussions associated with the genetic liability to PTSD across 4 biobanks. The circulatory and respiratory systems association was observed to be overrepresented in all 4 biobanks.

Brainwide Mendelian Randomization Study of Anxiety Disorders and Symptoms

BACKGROUND

To gain insights into the role of brain structure and function on anxiety (ANX), we conducted a genetically informed investigation leveraging information from ANX genome-wide association studies available from the UK Biobank (n = 380,379), the FinnGen Program (n = 290,361), and the Million Veteran Program (n = 175,163) together with UK Biobank genome-wide data (n = 33,224) related to 3935 brain imaging-derived phenotypes (IDPs).

METHODS

A genetic correlation analysis between ANX and brain IDPs was performed using linkage disequilibrium score regression. To investigate ANX-brain associations, a 2-sample Mendelian randomization was performed considering multiple methods and sensitivity analyses. A subsequent multivariable Mendelian randomization was conducted to distinguish between direct and indirect effects. Finally, a generalized linear model was used to explore the associations of brain IDPs with ANX symptoms.

RESULTS

After false discovery rate correction (q < .05), we identified 41 brain IDPs genetically correlated with ANX without heterogeneity among the datasets investigated (i.e., UK Biobank, FinnGen, and Million Veteran Program). Six of these IDPs showed genetically inferred causal effects on ANX. In the subsequent multivariable Mendelian randomization analysis, reduced area of the right posterior middle cingulate gyrus (β = -0.09, p = 8.01 × 10-4) and reduced gray matter volume of the right anterior superior temporal gyrus (β = -0.09, p = 1.55 × 10-3) had direct effects on ANX. In the ANX symptom-level analysis, the right posterior middle cingulate gyrus was negatively associated with "tense, sore, or aching muscles during the worst period of anxiety" (β = -0.13, p = 8.26 × 10-6).

CONCLUSIONS

This study identified genetically inferred effects that are generalizable across large cohorts, thereby contributing to our understanding of how changes in brain structure and function can lead to ANX.

Household income does not affect the pleiotropy of schizophrenia genetic liability with mental and physical health outcomes

Background and Hypothesis

Individuals with schizophrenia (SCZ) suffer from comorbidities that substantially reduce their life expectancy. Socioeconomic inequalities could contribute to many of the negative health outcomes associated with SCZ.

Study Design

We investigated genome-wide datasets related to SCZ (52,017 cases and 75,889 controls) from the Psychiatric Genomics Consortium, household income (HI; N=361,687) from UK Biobank, and 2,202 medical endpoints assessed in up to 342,499 FinnGen participants. A phenome-wide genetic correlation analysis of SCZ and HI was performed, also assessing whether SCZ genetic correlations were influenced by HI effect on SCZ. Additionally, SCZ and HI direct effects on medical endpoints were estimated using multivariable Mendelian randomization (MR).

Study Results

SCZ and HI showed overlapping genetic correlations with 70 traits (p<2.89×10 -5 ), including mental health, substance use, gastrointestinal illnesses, reproductive outcomes, liver diseases, respiratory problems, and musculoskeletal phenotypes. SCZ genetic correlations with these traits were not affected by HI effect on SCZ. Considering Bonferroni multiple testing correction (p<7.14×10 -4 ), MR analysis indicated that SCZ and HI may affect medical abortion (SCZ odds ratio, OR=1.07; HI OR=0.78), panic disorder (SCZ OR=1.20; HI OR=0.60), personality disorders (SCZ OR=1.31; HI OR=0.67), substance use (SCZ OR=1.2; HI OR=0.68), and adjustment disorders (SCZ OR=1.18; HI OR=0.78). Multivariable MR analysis confirmed that SCZ effects on these outcomes were independent of HI.

Conclusions

The effect of SCZ genetic liability on mental and physical health may not be strongly affected by socioeconomic differences. This suggests that SCZ-specific strategies are needed to reduce negative health outcomes affecting patients and high-risk individuals.

Brain-Wide Mendelian Randomization Study of Anxiety Disorders and Symptoms

Background

To gain insights into the role of brain structure and function on anxiety (ANX), we conducted a genetically informed investigation leveraging information from ANX genome-wide association studies available from UK Biobank (UKB; N=380,379), FinnGen Program (N=290,361), and Million Veteran Program (MVP; N=199,611) together with UKB genome-wide data (N=33,224) related to 3,935 brain imaging-derived phenotypes (IDP).

Methods

A genetic correlation analysis between ANX and brain IDPs was performed using linkage disequilibrium score regression. To investigate ANX-brain associations, a two-sample Mendelian randomization (MR) was performed considering multiple methods and sensitivity analyses. A subsequent multivariable MR (MVMR) was executed to distinguish between direct and indirect effects. Finally, a generalized linear model was used to explore the associations of brain IDPs with ANX symptoms.

Results

After false discovery rate correction (FDR q<0.05), we identified 41 brain IDPs genetically correlated with ANX without heterogeneity among the datasets investigated (i.e., UKB, FinnGen, and MVP). Six of these IDPs showed genetically inferred causal effects on ANX. In the subsequent MVMR analysis, reduced area of the right posterior middle-cingulate gyrus (rpMCG; beta=-0.09, P= 8.01×10 -4 ) and reduced gray-matter volume of the right anterior superior temporal gyrus (raSTG; beta=-0.09, P=1.55×10 -3 ) had direct effects on ANX. In the ANX symptom-level analysis, rpMCG was negatively associated with "tense sore oraching muscles during the worst period of anxiety" (beta=-0.13, P=8.26×10 -6 ).

Conclusions

This study identified genetically inferred effects generalizable across large cohorts, contributing to understand how changes in brain structure and function can lead to ANX.

Discovery of 95 PTSD loci provides insight into genetic architecture and neurobiology of trauma and stress-related disorders

Posttraumatic stress disorder (PTSD) genetics are characterized by lower discoverability than most other psychiatric disorders. The contribution to biological understanding from previous genetic studies has thus been limited. We performed a multi-ancestry meta-analysis of genome-wide association studies across 1,222,882 individuals of European ancestry (137,136 cases) and 58,051 admixed individuals with African and Native American ancestry (13,624 cases). We identified 95 genome-wide significant loci (80 novel). Convergent multi-omic approaches identified 43 potential causal genes, broadly classified as neurotransmitter and ion channel synaptic modulators (e.g., GRIA1, GRM8, CACNA1E ), developmental, axon guidance, and transcription factors (e.g., FOXP2, EFNA5, DCC ), synaptic structure and function genes (e.g., PCLO, NCAM1, PDE4B ), and endocrine or immune regulators (e.g., ESR1, TRAF3, TANK ). Additional top genes influence stress, immune, fear, and threat-related processes, previously hypothesized to underlie PTSD neurobiology. These findings strengthen our understanding of neurobiological systems relevant to PTSD pathophysiology, while also opening new areas for investigation.

Sex differences in the pleiotropy of hearing difficulty with imaging-derived phenotypes: a brain-wide investigation

Background

Hearing difficulty (HD) is one of the major health burdens in older adults. While aging-related changes in the peripheral auditory system play an important role, genetic variation associated with brain structure and function could also be involved in HD predisposition.

Methods

We analyzed a large-scale HD genome-wide association study (GWAS; N total = 501,825, 56% females) and GWAS data related to 3,935 brain imaging-derived phenotypes (IDPs) assessed in up to 33,224 individuals (52% females) using multiple magnetic resonance imaging (MRI) modalities. To investigate HD pleiotropy with brain structure and function, we conducted genetic correlation, latent causal variable (LCV), Mendelian randomization (MR), and multivariable generalized linear regression analyses. Additionally, we performed local genetic correlation and multi-trait colocalization analyses to identify genomic regions and loci implicated in the pleiotropic mechanisms shared between HD and brain IDPs.

Results

We observed a widespread genetic correlation of HD with 120 IDPs in females, 89 IDPs in males, and 171 IDPs in the sex-combined analysis. The LCV analyses showed that some of these genetic correlations could be due to cause-effect relationships. For seven correlations, the causal effects were also confirmed by the MR approach: vessel volume→HD in the sex-combined analysis; hippocampus volume→HD, cerebellum grey matter volume→HD, primary visual cortex volume→HD, and HD→rfMRI-ICA100 node 46 in females; global mean thickness→HD and HD→mean orientation dispersion index in superior corona radiata in males. The local genetic correlation analyses identified 13 pleiotropic regions between HD and these seven IDPs. We also observed a colocalization signal for the rs13026575 variant between HD, primary visual cortex volume, and SPTBN1 transcriptomic regulation in females.

Conclusion

Brain structure and function may have a role in the sex differences in HD predisposition via possible cause-effect relationships and shared regulatory mechanisms.

Integrating Genome-wide information and Wearable Device Data to Explore the Link of Anxiety and Antidepressants with Heart Rate Variability

Background

Anxiety disorders are associated with decreased heart rate variability (HRV), but the underlying mechanisms remain elusive.

Methods

We selected individuals with whole-genome sequencing, Fitbit, and electronic health record data (N=920; 61,333 data points) from the All of Us Research Program. Anxiety PRS were derived with PRS-CS after meta-analyzing anxiety genome-wide association studies from three major cohorts-UK Biobank, FinnGen, and the Million Veterans Program (N Total =364,550). The standard deviation of average RR intervals (SDANN) was calculated using five-minute average RR intervals over full 24-hour heart rate measurements. Antidepressant exposure was defined as an active antidepressant prescription at the time of the HRV measurement in the EHR. The associations of daily SDANN measurements with the anxiety PRS, antidepressant classes, and antidepressant substances were tested. Participants with lifetime diagnoses of cardiovascular disorders, diabetes mellitus, and major depression were excluded in sensitivity analyses. One-sample Mendelian randomization (MR) was employed to assess potential causal effect of anxiety on SDANN.

Results

Anxiety PRS was independently associated with reduced SDANN (beta=-0.08; p=0.003). Of the eight antidepressant medications and four classes tested, venlafaxine (beta=-0.12, p=0.002) and bupropion (beta=-0.071, p=0.01), tricyclic antidepressants (beta=-0.177, p=0.0008), selective serotonin reuptake inhibitors (beta=-0.069; p=0.0008) and serotonin and norepinephrine reuptake inhibitors (beta=-0.16; p=2×10 -6 ) were associated with decreased SDANN. One-sample MR indicated an inverse effect of anxiety on SDANN (beta=-2.22, p=0.03).

Conclusions

Anxiety and antidepressants are independently associated with decreased HRV, and anxiety appears to exert a causal effect on HRV. Our observational findings provide novel insights into the impact of anxiety on HRV.

Phenome-wide genetic-correlation analysis and genetically informed causal inference of amyotrophic lateral sclerosis

Leveraging genome-wide association statistics generated from a large study of amyotrophic lateral sclerosis (ALS; 29,612 cases and 122,656 controls) and UK Biobank (UKB; 4,024 phenotypes, up to 361,194 participants), we conducted a phenome-wide analysis of ALS genetic liability and identified 46 genetically correlated traits, such as fluid intelligence score (rg = - 0.21, p = 1.74 × 10-6), "spending time in pub or social club" (rg = 0.24, p = 2.77 × 10-6), non-work related walking (rg = - 0.25, p = 1.95 × 10-6), college education (rg = - 0.15, p = 7.08 × 10-5), "ever diagnosed with panic attacks (rg = 0.39, p = 4.24 × 10-5), and "self-reported other gastritis including duodenitis" (rg = 0.28, p = 1.4 × 10-3). To assess the putative directionality of these genetic correlations, we conducted a latent causal variable analysis, identifying significant genetic causality proportions (gĉp) linking ALS genetic liability to seven traits. While the genetic component of "self-reported other gastritis including duodenitis" showed a causal effect on ALS (gĉp = 0.50, p = 1.26 × 10-29), the genetic liability to ALS is potentially causal for multiple traits, also including an effect on "ever being diagnosed with panic attacks" (gĉp = 0.79, p = 5.011 × 10-15) and inverse effects on "other leisure/social group activities" (gĉp = 0.66, p = 1 × 10-4) and prospective memory result (gĉp = 0.35, p = 0.005). Our subsequent Mendelian randomization analysis indicated that some of these associations may be due to bidirectional effects. In conclusion, this phenome-wide investigation of ALS polygenic architecture highlights the widespread pleiotropy linking this disorder with several health domains.

Unraveling COVID-19 Relationship with Anxiety Disorders and Symptoms

Background

While COVID-19 outcomes are associated with increased anxiety, individuals affected by anxiety disorders are more likely to develop severe COVID-19 outcomes.

Methods

We used genome-wide data from UK Biobank (up to 420,531 participants), FinnGen Project (up to 329,077 participants), Million Veteran Program (175,163 participants), and COVID-19 Host Genetics Initiative (up to 122,616 cases and 2,475,240 controls) to investigate possible causal effects and shared genetic mechanisms linking COVID-19 outcomes to anxiety disorders and symptoms.

Results

We observed a strong genetic correlation of anxiety disorder with COVID-19 positive status (rg=0.35, p=2 × 10 -4 ) and COVID-19 hospitalization (rg=0.31, p=7.2 × 10 -4 ). Among anxiety symptoms, "Tense, sore, or aching muscles during worst period of anxiety" was genetically correlated with COVID-19 positive status (rg=0.33, p=0.001), while "Frequent trouble falling or staying asleep during worst period of anxiety" was genetically correlated with COVID-19 hospitalization (rg=0.24, p=0.004). Through a latent causal variable analysis, we observed that COVID-19 outcomes have statistically significant genetic causality proportion (gcp) on anxiety symptoms (e.g., COVID-19 positive status→"Recent easy annoyance or irritability" │gcp│=0.18, p=6.72 × 10 -17 ). Conversely, anxiety disorders appear to have a possible causal effect on COVID-19 (│gcp│=0.38, p=3.17 × 10 -9 ). Additionally, we also identified multiple loci with evidence of local genetic correlation between anxiety and COVID-19. These appear to be related to genetic effects shared with lung function, brain morphology, alcohol and tobacco use, and hematologic parameters.

Conclusions

This study provided important insights into the relationship between COVID-19 and mental health, differentiating the dynamics linking anxiety disorders to COVID-19 from the effect of COVID-19 on anxiety symptoms.

Characterizing the polygenic architecture of complex traits in populations of East Asian and European descent

To investigate the polygenicity of complex traits in populations of East Asian (EAS) and European (EUR) descents, we leveraged genome-wide data from Biobank Japan, UK Biobank, and FinnGen cohorts. Specifically, we analyzed up to 215 outcomes related to 18 health domains, assessing their polygenic architecture via descriptive statistics, such as the proportion of susceptibility SNPs per trait (πc). While we did not observe EAS-EUR differences in the overall distribution of polygenicity parameters across the phenotypes investigated, there were ancestry-specific patterns in the polygenicity differences between health domains. In EAS, pairwise comparisons across health domains showed enrichment for πc differences related to hematological and metabolic traits (hematological fold-enrichment = 4.45, p = 2.15 × 10-7; metabolic fold-enrichment = 4.05, p = 4.01 × 10-6). For both categories, the proportion of susceptibility SNPs was lower than that observed for several other health domains (EAS-hematological median πc = 0.15%, EAS-metabolic median πc = 0.18%) with the strongest πc difference with respect to respiratory traits (EAS-respiratory median πc = 0.50%; hematological-p = 2.26 × 10-3; metabolic-p = 3.48 × 10-3). In EUR, pairwise comparisons showed multiple πc differences related to the endocrine category (fold-enrichment = 5.83, p = 4.76 × 10-6), where these traits showed a low proportion of susceptibility SNPs (EUR-endocrine median πc = 0.01%) with the strongest difference with respect to psychiatric phenotypes (EUR-psychiatric median πc = 0.50%; p = 1.19 × 10-4). Simulating sample sizes of 1,000,000 and 5,000,000 individuals, we also showed that ancestry-specific polygenicity patterns translate into differences across health domains in the genetic variance explained by susceptibility SNPs projected to be genome-wide significant (e.g., EAS hematological-neoplasm p = 2.18 × 10-4; EUR endocrine-gastrointestinal p = 6.80 × 10-4). These findings highlight that traits related to the same health domains may present ancestry-specific variability in their polygenicity.

Characterizing the polygenic architecture of complex traits in populations of East Asian and European descent

To investigate the polygenicity of complex traits in populations of East Asian (EAS) and European (EUR) descents, we leveraged genome-wide data from Biobank Japan, UK Biobank, and FinnGen cohorts. Specifically, we analyzed up to 215 outcomes related to 18 health domains, assessing their polygenic architecture via descriptive statistics, such as the proportion of susceptibility SNPs per trait (π _c ). While we did not observe EAS-EUR differences in the overall distribution of polygenicity parameters across the phenotypes investigated, there were ancestry-specific patterns in the polygenicity differences between health domains. In EAS, pairwise comparisons across health domains showed enrichment for π _c differences related to hematological and metabolic traits (hematological fold-enrichment=4.45, p=2.15×10 ^-7 ; metabolic fold-enrichment=4.05, p=4.01×10 ^-6 ). For both categories, the proportion of susceptibility SNPs was lower than that observed for several other health domains (EAS-hematological median π _c =0.15%, EAS-metabolic median π _c =0.18%) with the strongest π _c difference with respect to respiratory traits (EAS-respiratory median π _c =0.50%; Hematological-p=2.26×10 ^-3 ; Metabolic-p=3.48×10 ^-3 ). In EUR, pairwise comparisons showed multiple π _c differences related to the endocrine category (fold-enrichment=5.83, p=4.76×10 ^-6 ), where these traits showed a low proportion of susceptibility SNPs (EUR-endocrine median π _c =0.01%) with the strongest difference with respect to psychiatric phenotypes (EUR-psychiatric median π _c =0.50%; p=1.19×10 ^-4 ). Simulating sample sizes of 1,000,000 and 5,000,000 individuals, we also showed that ancestry-specific polygenicity patterns translate into differences across health domains in the genetic variance explained by susceptibility SNPs projected to be genome-wide significant (e.g., EAS hematological-neoplasm p=2.18×10 ^-4 ; EUR endocrine-gastrointestinal p=6.80×10 ^-4 ). These findings highlight that traits related to the same health domains may present ancestry-specific variability in their polygenicity.

Modeling the longitudinal changes of ancestry diversity in the Million Veteran Program

BACKGROUND

The Million Veteran Program (MVP) participants represent 100 years of US history, including significant social and demographic changes over time. Our study assessed two aspects of the MVP: (i) longitudinal changes in population diversity and (ii) how these changes can be accounted for in genome-wide association studies (GWAS). To investigate these aspects, we divided MVP participants into five birth cohorts (N-range = 123,888 [born from 1943 to 1947] to 136,699 [born from 1948 to 1953]).

RESULTS

Ancestry groups were defined by (i) HARE (harmonized ancestry and race/ethnicity) and (ii) a random-forest clustering approach using the 1000 Genomes Project and the Human Genome Diversity Project (1kGP + HGDP) reference panels (77 world populations representing six continental groups). In these groups, we performed GWASs of height, a trait potentially affected by population stratification. Birth cohorts demonstrate important trends in ancestry diversity over time. More recent HARE-assigned Europeans, Africans, and Hispanics had lower European ancestry proportions than older birth cohorts (0.010 < Cohen's d < 0.259, p < 7.80 × 10-4). Conversely, HARE-assigned East Asians showed an increase in European ancestry proportion over time. In GWAS of height using HARE assignments, genomic inflation due to population stratification was prevalent across all birth cohorts (linkage disequilibrium score regression intercept = 1.08 ± 0.042). The 1kGP + HGDP-based ancestry assignment significantly reduced the population stratification (mean intercept reduction = 0.045 ± 0.007, p < 0.05) confounding in the GWAS statistics.

CONCLUSIONS

This study provides a characterization of ancestry diversity of the MVP cohort over time and compares two strategies to infer genetically defined ancestry groups by assessing differences in controlling population stratification in genome-wide association studies.

Declining autozygosity over time: An exploration in over 1 million individuals from three diverse cohorts

Previous studies have hypothesized that autozygosity is decreasing over generational time. However, these studies were limited to relatively small samples (n < 11,000) lacking in diversity, which may limit the generalizability of their findings. We present data that partially support this hypothesis from three large cohorts of diverse ancestries, two from the US (All of Us, n = 82,474; the Million Veteran Program, n = 622,497) and one from the UK (UK Biobank, n = 380,899). Our results from a mixed-effect meta-analysis demonstrate an overall trend of decreasing autozygosity over generational time (meta-analyzed slope = -0.029, SE = 0.009, p = 6.03e-4). On the basis of our estimates, we would predict FROH to decline 0.29% for every 20-year increase in birth year. We determined that a model including an ancestry-by-country interaction term fit the data best, indicating that ancestry differences in this trend differ by country. We found further evidence to suggest a difference between the US and UK cohorts by meta-analyzing within country, observing a significant negative estimate in the US cohorts (meta-analyzed slope = -0.058, SE = 0.015, p = 1.50e-4) but a non-significant estimate in the UK (meta-analyzed slope = -0.001, SE = 0.008, p = 0.945). The association between autozygosity and birth year was substantially attenuated when accounting for educational attainment and income (meta-analyzed slope = -0.011, SE = 0.008, p = 0.167), suggesting they may partially account for decreasing autozygosity over time. Overall, we demonstrate decreasing autozygosity over time in a large, modern sample and speculate that this trend can be attributed to increases in urbanization and panmixia and differences in sociodemographic processes lead to country-specific differences in the rate of decline.

Sex differences in the polygenic architecture of hearing problems in adults

BACKGROUND

Hearing problems (HP) in adults are common and are associated with several comorbid conditions. Its prevalence increases with age, reflecting the cumulative effect of environmental factors and genetic predisposition. Although several risk loci have been already identified, HP biology and epidemiology are still insufficiently investigated by large-scale genetic studies.

METHODS

Leveraging the UK Biobank, the Nurses' Health Studies (I and II), the Health Professionals Follow-up Study, and the Million Veteran Program, we conducted a comprehensive genome-wide investigation of HP in 748,668 adult participants (discovery N = 501,825; replication N = 226,043; cross-ancestry replication N = 20,800). We leveraged the GWAS findings to characterize HP polygenic architecture, exploring sex differences, polygenic risk across ancestries, tissue-specific transcriptomic regulation, cause-effect relationships with genetically correlated traits, and gene interactions with HP environmental risk factors.

RESULTS

We identified 54 risk loci and demonstrated that HP polygenic risk is shared across ancestry groups. Our transcriptomic regulation analysis highlighted the potential role of the central nervous system in HP pathogenesis. The sex-stratified analyses showed several additional associations related to peripheral hormonally regulated tissues reflecting a potential role of estrogen in hearing function. This evidence was supported by the multivariate interaction analysis that showed how genes involved in brain development interact with sex, noise pollution, and tobacco smoking in relation to their HP associations. Additionally, the genetically informed causal inference analysis showed that HP is linked to many physical and mental health outcomes.

CONCLUSIONS

The results provide many novel insights into the biology and epidemiology of HP in adults. Our sex-specific analyses and transcriptomic associations highlighted molecular pathways that may be targeted for drug development or repurposing. Additionally, the potential causal relationships identified may support novel preventive screening programs to identify individuals at risk.

Sex-Specific Genetic and Transcriptomic Liability to Neuroticism

BACKGROUND

The presentation, etiology, and relative risk of psychiatric disorders are strongly influenced by biological sex. Neuroticism is a transdiagnostic feature of psychiatric disorders displaying prominent sex differences. We performed genome-wide association studies of neuroticism separately in males and females to identify sex-specific genetic and transcriptomic profiles.

METHODS

Neuroticism scores were derived from the Eysenck Personality Inventory Neuroticism scale. Genome-wide association studies were performed in 145,669 females and 129,229 males from the UK Biobank considering autosomal and X chromosomal variation. Two-sided z tests were used to test for sex-specific effects of discovered loci, genetic correlates (n = 673 traits), tissue and gene transcriptomic profiles, and polygenic associations across health outcomes in the Vanderbilt University Biobank (39,692 females and 31,268 males).

RESULTS

The single nucleotide polymorphism heritability of neuroticism was not statistically different between males (h2 = 10.6%) and females (h2 = 11.85%). Four female-specific (rs10736549-CNTN5, rs6507056-ASXL3, rs2087182-MMS22L, and rs72995548-HSPB2) and 2 male-specific (rs10507274-MED13L and rs7984597) neuroticism risk loci reached genome-wide significance. Male- and female-specific neuroticism polygenic scores were most significantly associated with mood disorders (males: odds ratio = 1.11, p = 1.40 × 10-9; females: odds ratio = 1.14, p = 6.05 × 10-22). They also associated with sex-specific laboratory measurements related to erythrocyte count, distribution, and hemoglobin concentration. Gene expression variation in the pituitary was enriched for neuroticism loci in males (male: b = 0.026, p = .002), and genetically regulated transcriptomic changes highlighted the effect of SHISHA9, TEX26, and NCOA6.

CONCLUSIONS

Through a comprehensive assessment of genetic risk for neuroticism and the associated biological processes, this study identified several molecular pathways that can partially explain the known sex differences in neurotic symptoms and their psychiatric comorbidities.

Epidemiologic and Genetic Associations of Endometriosis With Depression, Anxiety, and Eating Disorders

IMPORTANCE

Endometriosis is a common chronic gynecologic pathology with a large negative impact on women's health. Beyond severe physical symptoms, endometriosis is also associated with several psychiatric comorbidities, including depression and anxiety.

OBJECTIVE

To investigate whether pleiotropy contributes to the association of endometriosis with depression, anxiety, and eating disorders.

DESIGN, SETTING, AND PARTICIPANTS

This genetic association study was performed between September 13, 2021, and June 24, 2022, in 202 276 unrelated female participants. Genotypic and phenotypic information from the UK Biobank was combined with genome-wide association statistics available from the Psychiatric Genomics Consortium (11 countries), the Million Veteran Program (US), the FinnGen study (Finland), and the CHARGE (Cohorts for Heart and Aging Research in Genomic Epidemiology) consortium (5 countries).

MAIN OUTCOMES AND MEASURES

The main outcomes were the phenotypic and genetic associations of endometriosis with anxiety, depression, and eating disorders.

RESULTS

A total of 8276 women with endometriosis (mean [SD] age, 53.1 [7.9] years) and 194 000 female controls (mean [SD] age, 56.7 [7.9] years) were included in the study. In a multivariate regression analysis accounting for age, body mass index, socioeconomic status, chronic pain-related phenotypes, irritable bowel syndrome, and psychiatric comorbidities, endometriosis was associated with increased odds of depression (odds ratio [OR], 3.61; 95% CI, 3.32-3.92), eating disorders (OR, 2.94; 95% CI, 1.96-4.41), and anxiety (OR, 2.61; 95% CI, 2.30-2.97). These associations were supported by consistent genetic correlations (rg) (depression rg, 0.36, P = 1.5 × 10-9; anxiety rg, 0.33, P = 1.17 × 10-5; and eating disorders rg, 0.61, P = .02). With the application of a 1-sample mendelian randomization, the genetic liabilities to depression and anxiety were associated with increased odds of endometriosis (depression: OR, 1.09; 95% CI, 1.08-1.11; anxiety: OR, 1.39; 95% CI, 1.13-1.65). A genome-wide analysis of pleiotropic associations shared between endometriosis and psychiatric disorders identified 1 locus, DGKB rs12666606, with evidence of pleiotropy between endometriosis and depression after multiple testing correction (z = -9.46 for endometriosis, z = 8.10 for depression, P = 5.56 × 10-8; false discovery rate q = 4.95 × 10-4).

CONCLUSIONS AND RELEVANCE

These findings highlight that endometriosis is associated with women's mental health through pleiotropic mechanisms. To our knowledge, this is the first large-scale study to provide genetic and phenotypic evidence of the processes underlying the psychiatric comorbidities of endometriosis.

Phenome-wide association study of loci harboring de novo tandem repeat mutations in UK Biobank exomes

When present in coding regions, tandem repeats (TRs) may have large effects on protein structure and function contributing to health and disease. We use a family-based design to identify de novo TRs and assess their impact at the population level in 148,607 European ancestry participants from the UK Biobank. The 427 loci with de novo TR mutations are enriched for targets of microRNA-184 (21.1-fold, P = 4.30 × 10^-5, FDR = 9.50 × 10^-3). There are 123 TR-phenotype associations with posterior probabilities > 0.95. These relate to body structure, cognition, and cardiovascular, metabolic, psychiatric, and respiratory outcomes. We report several loci with large likely causal effects on tissue microstructure, including the FAN1-[TG]_N and carotid intima-media thickness (mean thickness: beta = 5.22, P = 1.22 × 10^-6, FDR = 0.004; maximum thickness: beta = 6.44, P = 1.12 × 10^-6, FDR = 0.004). Two exonic repeats FNBP4-[GGT]_N and BTN2A1-[CCT]_N alter protein structure. In this work, we contribute clear and testable hypotheses of dose-dependent TR implications linking genetic variation and protein structure with health and disease outcomes.

A MUC5B Gene Polymorphism, rs35705950-T, Confers Protective Effects Against COVID-19 Hospitalization but Not Severe Disease or Mortality

Rationale: A common MUC5B gene polymorphism, rs35705950-T, is associated with idiopathic pulmonary fibrosis (IPF), but its role in severe acute respiratory syndrome coronavirus 2 infection and disease severity is unclear. Objectives: To assess whether rs35705950-T confers differential risk for clinical outcomes associated with coronavirus disease (COVID-19) infection among participants in the Million Veteran Program (MVP). Methods: The MUC5B rs35705950-T allele was directly genotyped among MVP participants; clinical events and comorbidities were extracted from the electronic health records. Associations between the incidence or severity of COVID-19 and rs35705950-T were analyzed within each ancestry group in the MVP followed by transancestry meta-analysis. Replication and joint meta-analysis were conducted using summary statistics from the COVID-19 Host Genetics Initiative (HGI). Sensitivity analyses with adjustment for additional covariates (body mass index, Charlson comorbidity index, smoking, asbestosis, rheumatoid arthritis with interstitial lung disease, and IPF) and associations with post-COVID-19 pneumonia were performed in MVP subjects. Measurements and Main Results: The rs35705950-T allele was associated with fewer COVID-19 hospitalizations in transancestry meta-analyses within the MVP (N_cases = 4,325; N_controls = 507,640; OR = 0.89 [0.82-0.97]; P = 6.86 × 10^-3) and joint meta-analyses with the HGI (N_cases = 13,320; N_controls = 1,508,841; OR, 0.90 [0.86-0.95]; P = 8.99 × 10^-5). The rs35705950-T allele was not associated with reduced COVID-19 positivity in transancestry meta-analysis within the MVP (N_cases = 19,168/N_controls = 492,854; OR, 0.98 [0.95-1.01]; P = 0.06) but was nominally significant (P < 0.05) in the joint meta-analysis with the HGI (N_cases = 44,820; N_controls = 1,775,827; OR, 0.97 [0.95-1.00]; P = 0.03). Associations were not observed with severe outcomes or mortality. Among individuals of European ancestry in the MVP, rs35705950-T was associated with fewer post-COVID-19 pneumonia events (OR, 0.82 [0.72-0.93]; P = 0.001). Conclusions: The MUC5B variant rs35705950-T may confer protection in COVID-19 hospitalizations.

Denisovan and Neanderthal archaic introgression differentially impacted the genetics of complex traits in modern populations

BACKGROUND

Introgression from extinct Neanderthal and Denisovan human species has been shown to contribute to the genetic pool of modern human populations and their phenotypic spectrum. Evidence of how Neanderthal introgression shaped the genetics of human traits and diseases has been extensively studied in populations of European descent, with signatures of admixture reported for instance in genes associated with pigmentation, immunity, and metabolic traits. However, limited information is currently available about the impact of archaic introgression on other ancestry groups. Additionally, to date, no study has been conducted with respect to the impact of Denisovan introgression on the health and disease of modern populations. Here, we compare the way evolutionary pressures shaped the genetics of complex traits in East Asian and European populations, and provide evidence of the impact of Denisovan introgression on the health of East Asian and Central/South Asian populations.

RESULTS

Leveraging genome-wide association statistics from the Biobank Japan and UK Biobank, we assessed whether Denisovan and Neanderthal introgression together with other evolutionary genomic signatures were enriched for the heritability of physiological and pathological conditions in populations of East Asian and European descent. In EAS, Denisovan-introgressed loci were enriched for coronary artery disease heritability (1.69-fold enrichment, p=0.003). No enrichment for archaic introgression was observed in EUR. We also performed a phenome-wide association study of Denisovan and Neanderthal alleles in six ancestry groups available in the UK Biobank. In EAS, the Denisovan-introgressed SNP rs62391664 in the major histocompatibility complex region was associated with albumin/globulin ratio (beta=-0.17, p=3.57×10^-7). Neanderthal-introgressed alleles were associated with psychiatric and cognitive traits in EAS (e.g., "No Bipolar or Depression"-rs79043717 beta=-1.5, p=1.1×10^-7), and with blood biomarkers (e.g., alkaline phosphatase-rs11244089 beta=0.1, p=3.69×10^-116) and red hair color (rs60733936 beta=-0.86, p=4.49×10^-165) in EUR. In the other ancestry groups, Neanderthal alleles were associated with several traits, also including the use of certain medications (e.g., Central/South East Asia: indapamide - rs732632 beta=-2.38, p=5.22×10^-7).

CONCLUSIONS

Our study provides novel evidence regarding the impact of archaic introgression on the genetics of complex traits in worldwide populations, highlighting the specific contribution of Denisovan introgression in EAS populations.

Genetic variants associated with psychiatric disorders are enriched at epigenetically active sites in lymphoid cells

Multiple psychiatric disorders have been associated with abnormalities in both the innate and adaptive immune systems. The role of these abnormalities in pathogenesis, and whether they are driven by psychiatric risk variants, remains unclear. We test for enrichment of GWAS variants associated with multiple psychiatric disorders (cross-disorder or trans-diagnostic risk), or 5 specific disorders (cis-diagnostic risk), in regulatory elements in immune cells. We use three independent epigenetic datasets representing multiple organ systems and immune cell subsets. Trans-diagnostic and cis-diagnostic risk variants (for schizophrenia and depression) are enriched at epigenetically active sites in brain tissues and in lymphoid cells, especially stimulated CD4⁺ T cells. There is no evidence for enrichment of either trans-risk or cis-risk variants for schizophrenia or depression in myeloid cells. This suggests a possible model where environmental stimuli activate T cells to unmask the effects of psychiatric risk variants, contributing to the pathogenesis of mental health disorders.

The association of obesity-related traits on COVID-19 severity and hospitalization is affected by socio-economic status: a multivariable Mendelian randomization study

BACKGROUND

Due to its large impact on human health, socio-economic status (SES) could at least partially influence the established association between obesity and coronavirus disease 2019 (COVID-19) severity. To estimate the independent effect of body size and SES on the clinical manifestations of COVID-19, we conducted a Mendelian randomization (MR) study.

METHODS

Applying two-sample MR approaches, we evaluated the effects of body mass index (BMI, n = 322 154), waist circumference (WC, n = 234 069), hip circumference (n = 213 019) and waist-hip ratio (n = 210 088) with respect to three COVID-19 outcomes: severe respiratory COVID-19 (cases = 8779, controls = 1 000 875), hospitalized COVID-19 (cases = 17 992, controls = 1 810 493) and COVID-19 infection (cases = 87 870, controls = 2 210 804). Applying a multivariable MR (MVMR) approach, we estimated the effect of these anthropometric traits on COVID-19 outcomes accounting for the effect of SES assessed as household income (n = 286 301).

RESULTS

BMI and WC were associated with severe respiratory COVID-19 [BMI: odds ratio (OR) = 1.51, CI = 1.24-1.84, P = 3.01e-05; WC: OR = 1.48, 95% CI = 1.15-1.91, P = 0.0019] and hospitalized COVID-19 (BMI: OR = 1.50, 95% CI = 1.32-1.72, P = 8.83e-10; WC: OR = 1.41, 95% CI = 1.20-1.67, P = 3.72e-05). Conversely, income was associated with lower odds of severe respiratory (OR = 0.70, 95% CI = 0.53-0.93, P = 0.015) and hospitalized COVID-19 (OR = 0.78, 95% CI = 0.66-0.92, P = 0.003). MVMR analyses showed that the effect of these obesity-related traits on increasing the odds of COVID-19 negative outcomes becomes null when accounting for income. Conversely, the association of income with lower odds of COVID-19 negative outcomes is not affected when including the anthropometric traits in the multivariable model.

CONCLUSION

Our findings indicate that SES contributes to the effect of obesity-related traits on COVID-19 severity and hospitalization.

Understanding the comorbidity between posttraumatic stress severity and coronary artery disease using genome-wide information and electronic health records

The association between coronary artery disease (CAD) and posttraumatic stress disorder (PTSD) contributes to the high morbidity and mortality observed for these conditions. To understand the dynamics underlying PTSD-CAD comorbidity, we investigated large-scale genome-wide association (GWA) statistics from the Million Veteran Program (MVP), the UK Biobank (UKB), the Psychiatric Genomics Consortium, and the CARDIoGRAMplusC4D Consortium. We observed a genetic correlation of CAD with PTSD case-control and quantitative outcomes, ranging from 0.18 to 0.32. To investigate possible cause-effect relationships underlying these genetic correlations, we performed a two-sample Mendelian randomization (MR) analysis, observing a significant bidirectional relationship between CAD and PTSD symptom severity. Genetically-determined PCL-17 (PTSD 17-item Checklist) total score was associated with increased CAD risk (odds ratio = 1.04; 95% confidence interval, 95% CI = 1.01-1.06). Conversely, CAD genetic liability was associated with reduced PCL-17 total score (beta = -0.42; 95% CI = -0.04 to -0.81). Because of these opposite-direction associations, we conducted a pleiotropic meta-analysis to investigate loci with concordant vs. discordant effects on PCL-17 and CAD, observing that concordant-effect loci were enriched for molecular pathways related to platelet amyloid precursor protein (beta = 1.53, p = 2.97 × 10-7) and astrocyte activation regulation (beta = 1.51, p = 2.48 × 10-6) while discordant-effect loci were enriched for biological processes related to lipid metabolism (e.g., triglyceride-rich lipoprotein particle clearance, beta = 2.32, p = 1.61 × 10-10). To follow up these results, we leveraged MVP and UKB electronic health records (EHR) to assess longitudinal changes in the association between CAD and posttraumatic stress severity. This EHR-based analysis highlighted that earlier CAD diagnosis is associated with increased PCL-total score later in life, while lower PCL total score was associated with increased risk of a later CAD diagnosis (Mann-Kendall trend test: MVP tau = 0.932, p < 2 × 10-16; UKB tau = 0.376, p = 0.005). In conclusion, both our genetically-informed analyses and our EHR-based follow-up investigation highlighted a bidirectional relationship between PTSD and CAD where multiple pleiotropic mechanisms are likely to be involved.

Association of Kidney Comorbidities and Acute Kidney Failure With Unfavorable Outcomes After COVID-19 in Individuals With the Sickle Cell Trait

Importance

Sickle cell trait (SCT), defined as the presence of 1 hemoglobin beta sickle allele (rs334-T) and 1 normal beta allele, is prevalent in millions of people in the US, particularly in individuals of African and Hispanic ancestry. However, the association of SCT with COVID-19 is unclear.

Objective

To assess the association of SCT with the prepandemic health conditions in participants of the Million Veteran Program (MVP) and to assess the severity and sequelae of COVID-19.

Design, Setting, and Participants

COVID-19 clinical data include 2729 persons with SCT, of whom 353 had COVID-19, and 129 848 SCT-negative individuals, of whom 13 488 had COVID-19. Associations between SCT and COVID-19 outcomes were examined using firth regression. Analyses were performed by ancestry and adjusted for sex, age, age squared, and ancestral principal components to account for population stratification. Data for the study were collected between March 2020 and February 2021.

Exposures

The hemoglobin beta S (HbS) allele (rs334-T).

Main Outcomes and Measures

This study evaluated 4 COVID-19 outcomes derived from the World Health Organization severity scale and phenotypes derived from International Classification of Diseases codes in the electronic health records.

Results

Of the 132 577 MVP participants with COVID-19 data, mean (SD) age at the index date was 64.8 (13.1) years. Sickle cell trait was present in 7.8% of individuals of African ancestry and associated with a history of chronic kidney disease, diabetic kidney disease, hypertensive kidney disease, pulmonary embolism, and cerebrovascular disease. Among the 4 clinical outcomes of COVID-19, SCT was associated with an increased COVID-19 mortality in individuals of African ancestry (n = 3749; odds ratio, 1.77; 95% CI, 1.13 to 2.77; P = .01). In the 60 days following COVID-19, SCT was associated with an increased incidence of acute kidney failure. A counterfactual mediation framework estimated that on average, 20.7% (95% CI, -3.8% to 56.0%) of the total effect of SCT on COVID-19 fatalities was due to acute kidney failure.

Conclusions and Relevance

In this genetic association study, SCT was associated with preexisting kidney comorbidities, increased COVID-19 mortality, and kidney morbidity.

An Atlas of Genetic Correlations and Genetically Informed Associations Linking Psychiatric and Immune-Related Phenotypes

Importance

Certain psychiatric and immune-related disorders are reciprocal risk factors. However, the nature of these associations is unclear.

Objective

To characterize the pleiotropy between psychiatric and immune-related traits, as well as risk factors of hypothesized relevance.

Design, Setting, and Participants

This genetic association study was conducted from July 10, 2020, to January 15, 2022. Analyses used genome-wide association (GWA) statistics related to 14 psychiatric traits; 13 immune-related phenotypes, ie, allergic, autoimmune, and inflammatory disorders; and 15 risk factors related to health-related behaviors, social determinants of health, and stress response. Genetically correlated psychiatric-immune pairs were assessed using 2-sample mendelian randomization (MR) with sensitivity analyses and multivariable adjustment for genetic associations of third variables. False discovery rate correction (Q value < .05) was applied for each analysis.

Exposures

Genetic associations.

Main Outcomes and Measures

Genetic correlations and MR association estimates with SEs and P values. A data-driven approach was used that did not test a priori planned hypotheses.

Results

A total of 44 genetically correlated psychiatric-immune pairs were identified, including 31 positive correlations (most consistently involving asthma, Crohn disease, hypothyroidism, and ulcerative colitis) and 13 negative correlations (most consistently involving allergic rhinitis and type 1 diabetes). Correlations with third variables were especially strong for psychiatric phenotypes. MR identified 7 associations of psychiatric phenotypes on immune-related phenotypes that were robust to multivariable adjustment, including the positive association of (1) the psychiatric cross-disorder phenotype with asthma (odds ratio [OR], 1.04; 95% CI, 1.02-1.06), Crohn disease (OR, 1.09; 95% CI, 1.05-1.14), and ulcerative colitis (OR, 1.09; 95% CI, 1.05-1.14); (2) major depression with asthma (OR, 1.25; 95% CI, 1.13-1.37); (3) schizophrenia with Crohn disease (OR, 1.12; 95% CI, 1.05-1.18) and ulcerative colitis (OR, 1.14; 95% CI, 1.07-1.21); and a negative association of risk tolerance with allergic rhinitis (OR, 0.77; 95% CI, 0.67-0.92).

Conclusions and Relevance

Results of this genetic association study suggest that genetic liability for psychiatric disorders was associated with liability for several immune disorders, suggesting that vertical pleiotropy related to behavioral traits (or correlated third variables) contributes to clinical associations observed in population-scale data.

The integration of genetically-regulated transcriptomics and electronic health records highlights a pattern of medical outcomes related to increased hepatic transthyretin expression

Transthyretin (TTR) is the precursor of the fibrils that compromise organ function in hereditary and sporadic systemic amyloidoses (ATTR). RNA-interference and anti-sense therapeutics targeting TTR hepatic transcription have been shown to reduce TTR amyloid formation. In the present study, we leveraged genetic and phenotypic information from the UK Biobank and transcriptomic profiles from the Genotype-Tissue Expression project to test the association of genetically regulated TTR gene expression with 7149 traits assessed in 420,531 individuals. We conducted a multi-tissue analysis of TTR transcription and identified an association with a operational procedure related to bone fracture (p = 5.46×10^-6). Using tissue-specific TTR expression information, we demonstrated that the association is driven by the genetic regulation of TTR hepatic expression (odds ratio [OR] = 3.46, p = 9.51×10^-5). Using the UK Biobank electronic health records (EHRs), we investigated the comorbidities affecting individuals undergoing this surgical procedure. Excluding bone fracture EHRs, we identified a pattern of health outcomes previously associated with ATTR manifestations. These included osteoarthritis (OR = 3.18, p = 9.18×10^-8), carpal tunnel syndrome (OR = 2.15, p = .002), and a history of gastrointestinal diseases (OR = 2.01, p = 8.07×10^-4). In conclusion, our study supports that TTR hepatic expression can affect health outcomes linked to physiological and pathological processes presumably related to the encoded protein.

A Phenome-Wide Association Study of genes associated with COVID-19 severity reveals shared genetics with complex diseases in the Million Veteran Program

The study aims to determine the shared genetic architecture between COVID-19 severity with existing medical conditions using electronic health record (EHR) data. We conducted a Phenome-Wide Association Study (PheWAS) of genetic variants associated with critical illness (n = 35) or hospitalization (n = 42) due to severe COVID-19 using genome-wide association summary data from the Host Genetics Initiative. PheWAS analysis was performed using genotype-phenotype data from the Veterans Affairs Million Veteran Program (MVP). Phenotypes were defined by International Classification of Diseases (ICD) codes mapped to clinically relevant groups using published PheWAS methods. Among 658,582 Veterans, variants associated with severe COVID-19 were tested for association across 1,559 phenotypes. Variants at the ABO locus (rs495828, rs505922) associated with the largest number of phenotypes (nrs495828 = 53 and nrs505922 = 59); strongest association with venous embolism, odds ratio (ORrs495828 1.33 (p = 1.32 x 10-199), and thrombosis ORrs505922 1.33, p = 2.2 x10-265. Among 67 respiratory conditions tested, 11 had significant associations including MUC5B locus (rs35705950) with increased risk of idiopathic fibrosing alveolitis OR 2.83, p = 4.12 × 10-191; CRHR1 (rs61667602) associated with reduced risk of pulmonary fibrosis, OR 0.84, p = 2.26× 10-12. The TYK2 locus (rs11085727) associated with reduced risk for autoimmune conditions, e.g., psoriasis OR 0.88, p = 6.48 x10-23, lupus OR 0.84, p = 3.97 x 10-06. PheWAS stratified by ancestry demonstrated differences in genotype-phenotype associations. LMNA (rs581342) associated with neutropenia OR 1.29 p = 4.1 x 10-13 among Veterans of African and Hispanic ancestry but not European. Overall, we observed a shared genetic architecture between COVID-19 severity and conditions related to underlying risk factors for severe and poor COVID-19 outcomes. Differing associations between genotype-phenotype across ancestries may inform heterogenous outcomes observed with COVID-19. Divergent associations between risk for severe COVID-19 with autoimmune inflammatory conditions both respiratory and non-respiratory highlights the shared pathways and fine balance of immune host response and autoimmunity and caution required when considering treatment targets.

Using phenotype risk scores to enhance gene discovery for generalized anxiety disorder and posttraumatic stress disorder

UK Biobank (UKB) is a key contributor in mental health genome-wide association studies (GWAS) but only ~31% of participants completed the Mental Health Questionnaire ("MHQ responders"). We predicted generalized anxiety disorder (GAD), posttraumatic stress disorder (PTSD), and major depression symptoms using elastic net regression in the ~69% of UKB participants lacking MHQ data ("MHQ non-responders"; N_Training = 50%; N_Test = 50%), maximizing the informative sample for these traits. MHQ responders were more likely to be female, from higher socioeconomic positions, and less anxious than non-responders. Genetic correlation of GAD and PTSD between MHQ responders and non-responders ranged from 0.636 to 1.08; both were predicted by polygenic scores generated from independent cohorts. In meta-analyses of GAD (N = 489,579) and PTSD (N = 497,803), we discovered many novel genomic risk loci (13 for GAD and 40 for PTSD). Transcriptomic analyses converged on altered regulation of prenatal dorsolateral prefrontal cortex in these disorders. Our results provide one roadmap by which sample size and statistical power may be improved for gene discovery of incompletely ascertained traits in the UKB and other biobanks with limited mental health assessment.

Attachment Style Moderates Polygenic Risk for Incident Posttraumatic Stress in U.S. Military Veterans: A 7-Year, Nationally Representative, Prospective Cohort Study

BACKGROUND

Posttraumatic stress disorder (PTSD) develops consequent to complex gene-by-environment interactions beyond the precipitating trauma. To date, however, no known study has used a prospective design to examine how polygenic risk scores (PRSs) interact with social-environmental factors such as attachment style to predict PTSD development.

METHODS

PRSs were derived from a genome-wide association study of PTSD symptoms (N = 186,689; Million Veteran Program cohort). We evaluated combined effects of PRS and attachment style in predicting incident PTSD in a 7-year, nationally representative cohort of trauma-exposed, European-American U.S. military veterans without PTSD (N = 1083). We also conducted multivariate gene-by-environment interaction and drug repositioning analyses to identify loci that interact with multiple environmental factors and potential pharmacotherapies that may be repurposed for this disorder.

RESULTS

Veterans with higher PTSD PRS were more likely to have an incident-positive screen for PTSD over 7 years. A gene-by-environment interaction was also observed, such that higher PRS only predicted incident PTSD in veterans with an insecure attachment style and not those with a secure attachment style. At an individual locus level, the strongest gene-by-environment interaction was observed for the rs4702 variant of the FURIN gene with cumulative lifetime trauma burden. Drug repositioning revealed that genes implicated in PRS are perturbated by the drug doxylamine.

CONCLUSIONS

Attachment style moderates polygenic risk for the development of PTSD in European-American veterans. These findings may inform PTSD prevention and treatment for veterans with high polygenic risk for PTSD and suggest a potential pharmacotherapeutic target for risk genes moderated by social-environmental factors.

Genetically regulated multi-omics study for symptom clusters of posttraumatic stress disorder highlights pleiotropy with hematologic and cardio-metabolic traits

Posttraumatic stress disorder (PTSD) is a psychiatric disorder that may arise in response to severe traumatic event and is diagnosed based on three main symptom clusters (reexperiencing, avoidance, and hyperarousal) per the Diagnostic Manual of Mental Disorders (version DSM-IV-TR). In this study, we characterized the biological heterogeneity of PTSD symptom clusters by performing a multi-omics investigation integrating genetically regulated gene, splicing, and protein expression in dorsolateral prefrontal cortex tissue within a sample of US veterans enrolled in the Million Veteran Program (N total = 186,689). We identified 30 genes in 19 regions across the three PTSD symptom clusters. We found nine genes to have cell-type specific expression, and over-representation of miRNA-families - miR-148, 30, and 8. Gene-drug target prioritization approach highlighted cyclooxygenase and acetylcholine compounds. Next, we tested molecular-profile based phenome-wide impact of identified genes with respect to 1678 phenotypes derived from the Electronic Health Records of the Vanderbilt University biorepository (N = 70,439). Lastly, we tested for local genetic correlation across PTSD symptom clusters which highlighted metabolic (e.g., obesity, diabetes, vascular health) and laboratory traits (e.g., neutrophil, eosinophil, tau protein, creatinine kinase). Overall, this study finds comprehensive genomic evidence including clinical and regulatory profiles between PTSD, hematologic and cardiometabolic traits, that support comorbidities observed in epidemiologic studies of PTSD.

Polygenic scores for empathy associate with posttraumatic stress severity in response to certain traumatic events

Background

Posttraumatic stress disorder (PTSD) is triggered by environmental stressors. Empathy may predispose an individual to respond to life events differently if high empathizers are emotionally more sensitive to trauma. For the first time, we test this hypothesis using genetic information.

Methods

We applied polygenic scoring (PGS) to investigate the shared genetics linking empathy (measured using the Empathy Quotient (EQ), a self-report measure of empathy; N = 46,861) and PTSD symptom severity (measured using the 6-item PTSD Checklist 6-item (PCL-6)) in the UK Biobank (N = 126,219). Follow-up analyses were performed in the context of (1) experiencing any of 16 potential traumas, (2) the total number of traumas endorsed, and (3) the context of trauma. Autism, depression, generalized anxiety, and PCL-17 PGS were included as covariates to verify the specificity of the effect.

Results

EQ_PGS associated with PCL-6 (R² = 0.012%, P = 9.35 × 10⁻⁵). This effect remained significant after accounting for autism, depression, PTSD, and anxiety PGS but was observed only in those who endorsed experiencing at least one traumatic event. EQ_PGS showed the strongest effect on PCL-6 (β = 2.32, s.e. = 0.762, P = 0.002) among those who endorsed childhood neglect/abuse (felt hated as a child). With respect to case status, the highest probability of PTSD was 17.93% and 10.04% for those who endorsed “feeling hated as a child” and those who did not, respectively (P_diff = 0.011; Cohen's d = 1.951, 95%CI 1.70–2.20).

Conclusions

A genetic predisposition to higher empathy, which may index greater emotional sensitivity, predisposes an individual to more severe PTSD symptoms, especially after early-life adversity.

Host Genetic Liability for Severe COVID-19 Associates with Alcohol Drinking Behavior and Diabetic Outcomes in Participants of European Descent

Risk factors and long-term consequences of COVID-19 infection are unclear but can be investigated with large-scale genomic data. To distinguish correlation from causation, we performed in-silico analyses of three COVID-19 outcomes (N > 1,000,000). We show genetic correlation and putative causality with depressive symptoms, metformin use (genetic causality proportion (gĉp) with severe respiratory COVID-19 = 0.576, p = 1.07 × 10⁻⁵ and hospitalized COVID-19 = 0.713, p = 0.003), and alcohol drinking status (gĉp with severe respiratory COVID-19 = 0.633, p = 7.04 × 10⁻⁵ and hospitalized COVID-19 = 0.848, p = 4.13 × 10⁻¹³). COVID-19 risk loci associated with several hematologic biomarkers. Comprehensive findings inform genetic contributions to COVID-19 epidemiology, molecular mechanisms, and risk factors and potential long-term health effects of severe response to infection.

Drinking and smoking polygenic risk is associated with childhood and early-adulthood psychiatric and behavioral traits independently of substance use and psychiatric genetic risk

Alcohol drinking and tobacco smoking are hazardous behaviors associated with a wide range of adverse health outcomes. In this study, we explored the association of polygenic risk scores (PRS) related to drinks per week, age of smoking initiation, smoking initiation, cigarettes per day, and smoking cessation with 433 psychiatric and behavioral traits in 4498 children and young adults (aged 8-21) of European ancestry from the Philadelphia neurodevelopmental cohort. After applying a false discovery rate multiple testing correction accounting for the number of PRS and traits tested, we identified 36 associations related to psychotic symptoms, emotion and age recognition social competencies, verbal reasoning, anxiety-related traits, parents' education, and substance use. These associations were independent of the genetic correlations among the alcohol-drinking and tobacco-smoking traits and those with cognitive performance, educational attainment, risk-taking behaviors, and psychopathology. The removal of participants endorsing substance use did not affect the associations of each PRS with psychiatric and behavioral traits identified as significant in the discovery analyses. Gene-ontology enrichment analyses identified several neurobiological processes underlying mechanisms of the PRS associations we report. In conclusion, we provide novel insights into the genetic overlap of smoking and drinking behaviors in children and young adults, highlighting their independence from psychopathology and substance use.

Disentangling sex differences in the shared genetic architecture of posttraumatic stress disorder, traumatic experiences, and social support with body size and composition

There is a well-known association of traumatic experiences and posttraumatic stress disorder (PTSD) with body size and composition, including consistent differences between sexes. However, the biology underlying these associations is unclear. To understand the genetic underpinnings of this complex relationship, we investigated genome-wide datasets informative of African and European ancestries from the Psychiatric Genomic Consortium, the UK Biobank, the GIANT Consortium, and the Million Veteran Program. We used genome-wide association statistics to estimate sex-specific genetic correlations (r g ) of traumatic experiences, social support, and PTSD with multiple anthropometric traits. After multiple testing corrections (false discovery rate, FDR q < 0.05), we observed 58 significant r g relationships in females (e.g., childhood physical abuse and body mass index, BMI r g  = 0.245, p = 3.88 × 10-10) and 21 significant r g relationships in males (e.g., been involved in combat or exposed to warzone and leg fat percentage; r g  = 0.405, p = 4.42 × 10-10). We performed causal inference analyses of these genetic overlaps using Mendelian randomization and latent causal variable approaches. Multiple female-specific putative causal relationships were observed linking body composition/size with PTSD (e.g., leg fat percentage→PTSD; beta = 0.319, p = 3.13 × 10-9), traumatic experiences (e.g., childhood physical abuse→waist circumference; beta = 0.055, p = 5.07 × 10-4), and childhood neglect (e.g., "someone to take you to doctor when needed as a child"→BMI; beta = -0.594, p = 1.09 × 10-5). In males, we observed putative causal effects linking anthropometric-trait genetic liabilities to traumatic experiences (e.g., BMI→childhood physical abuse; beta = 0.028, p = 8.19 × 10-3). Some of these findings were replicated in individuals of African descent although the limited sample size available did not permit us to conduct a sex-stratified analysis in this ancestry group. In conclusion, our findings provide insights regarding sex-specific causal networks linking anthropometric traits to PTSD, traumatic experiences, and social support.

ACE2 Netlas: In silico Functional Characterization and Drug-Gene Interactions of ACE2 Gene Network to Understand Its Potential Involvement in COVID-19 Susceptibility

Angiotensin-converting enzyme-2 (ACE2) receptor has been identified as the key adhesion molecule for the transmission of the SARS-CoV-2. However, there is no evidence that human genetic variation in ACE2 is singularly responsible for COVID-19 susceptibility. Therefore, we performed an integrative multi-level characterization of genes that interact with ACE2 (ACE2-gene network) for their statistically enriched biological properties in the context of COVID-19. The phenome-wide association of 51 genes including ACE2 with 4,756 traits categorized into 26 phenotype categories, showed enrichment of immunological, respiratory, environmental, skeletal, dermatological, and metabolic domains (p < 4e-4). Transcriptomic regulation of ACE2-gene network was enriched for tissue-specificity in kidney, small intestine, and colon (p < 4.7e-4). Leveraging the drug-gene interaction database we identified 47 drugs, including dexamethasone and spironolactone, among others. Considering genetic variants within ± 10 kb of ACE2-network genes we identified miRNAs whose binding sites may be altered as a consequence of genetic variation. The identified miRNAs revealed statistical over-representation of inflammation, aging, diabetes, and heart conditions. The genetic variant associations in RORA, SLC12A6, and SLC6A19 genes were observed in genome-wide association study (GWAS) of COVID-19 susceptibility. We also report the GWAS-identified variant in 3p21.31 locus, serves as trans-QTL for RORA and RORC genes. Overall, functional characterization of ACE2-gene network highlights several potential mechanisms in COVID-19 susceptibility. The data can also be accessed at https://gpwhiz.github.io/ACE2Netlas/.

Integrative genomic analyses identify susceptibility genes underlying COVID-19 hospitalization

Despite rapid progress in characterizing the role of host genetics in SARS-Cov-2 infection, there is limited understanding of genes and pathways that contribute to COVID-19. Here, we integrate a genome-wide association study of COVID-19 hospitalization (7,885 cases and 961,804 controls from COVID-19 Host Genetics Initiative) with mRNA expression, splicing, and protein levels (n = 18,502). We identify 27 genes related to inflammation and coagulation pathways whose genetically predicted expression was associated with COVID-19 hospitalization. We functionally characterize the 27 genes using phenome- and laboratory-wide association scans in Vanderbilt Biobank (n = 85,460) and identified coagulation-related clinical symptoms, immunologic, and blood-cell-related biomarkers. We replicate these findings across trans-ethnic studies and observed consistent effects in individuals of diverse ancestral backgrounds in Vanderbilt Biobank, pan-UK Biobank, and Biobank Japan. Our study highlights and reconfirms putative causal genes impacting COVID-19 severity and symptomology through the host inflammatory response.

Posttraumatic Stress Disorder and Cardiovascular Disease: State of the Science, Knowledge Gaps, and Research Opportunities

Posttraumatic stress disorder (PTSD) is characterized by a persistent maladaptive reaction after exposure to severe psychological trauma. Traumatic events that may precipitate PTSD include violent personal assaults, natural and human-made disasters, and exposure to military combat or warfare. There is a growing body of evidence for associations of PTSD with major risk factors for cardiovascular disease (CVD), such as hypertension and diabetes, as well as with major CVD outcomes, such as myocardial infarction and heart failure. However, it is unclear whether these associations are causal or confounded. Furthermore, the biological and behavioral mechanisms underlying these associations are poorly understood. Here, the available evidence on the association of PTSD with CVD from population, basic, and genomic research as well as from clinical and translational research are reviewed, seeking to identify major research gaps, barriers, and opportunities in knowledge acquisition and technology as well as research tools to support and accelerate critical research for near-term and longer-term translational research directions. Large-scale, well-designed prospective studies, capturing diverse and high-risk populations, are warranted that include uniform phenotyping of PTSD as well as broad assessment of biological and behavioral risk factors and CVD outcomes. Available evidence from functional brain imaging studies demonstrates that PTSD pathophysiology includes changes in specific anatomical brain regions and circuits, and studies of immune system function in individuals with PTSD suggest its association with enhanced immune inflammatory activity. However, establishment of animal models and human tissue biobanks is also warranted to elucidate the potential causal connection of PTSD-induced brain changes and/or inflammation with CVD pathophysiology. Emerging large-scale genome-wide association studies of PTSD will provide an opportunity to conduct mendelian randomization studies that test hypotheses regarding the presence, magnitude, and direction of causal associations between PTSD and CVD outcomes. By identifying research gaps in epidemiology and genomics, animal, and human translational research, opportunities to better justify and design future interventional trials are highlighted that may test whether treatment of PTSD or underlying neurobiological or immune dysregulation may improve or prevent CVD risk or outcomes.

Cross-ancestry genome-wide association studies identified heterogeneous loci associated with differences of allele frequency and regulome tagging between participants of European descent and other ancestry groups from the UK Biobank

To investigate cross-ancestry genetics of complex traits, we conducted a phenome-wide analysis of loci with heterogeneous effects across African, Admixed-American, Central/South Asian, East Asian, European and Middle Eastern participants of the UK Biobank (N = 441 331). Testing 843 phenotypes, we identified 82 independent genomic regions mapping variants showing genome-wide significant (GWS) associations (P < 5 × 10-8) in the trans-ancestry meta-analysis and GWS heterogeneity among the ancestry-specific effects. These included (i) loci with GWS association in one ancestry and concordant but heterogeneous effects among the other ancestries and (ii) loci with a GWS association in one ancestry group and an experiment-wide significant discordant effect (P < 6.1 × 10-4) in at least another ancestry. Since the trans-ancestry GWS associations were mostly driven by the European ancestry sample size, we investigated the differences of the allele frequency (ΔAF) and linkage disequilibrium regulome tagging (ΔLD) between European populations and the other ancestries. Within loci with concordant effects, the degree of heterogeneity was associated with European-Middle Eastern ΔAF (P = 9.04 × 10-6) and ΔLD of European populations with respect to African, Admixed-American and Central/South Asian groups (P = 8.21 × 10-4, P = 7.17 × 10-4 and P = 2.16 × 10-3, respectively). Within loci with discordant effects, ΔAF and ΔLD of European populations with respect to African and Central/South Asian ancestries were associated with the degree of heterogeneity (ΔAF: P = 7.69 × 10-3 and P = 5.31 × 10-3, ΔLD: P = 0.016 and P = 2.65 × 10-4, respectively). Considering the traits associated with cross-ancestry heterogeneous loci, we observed enrichments for blood biomarkers (P = 5.7 × 10-35) and physical appearance (P = 1.38 × 10-4). This suggests that these specific phenotypic classes may present considerable cross-ancestry heterogeneity owing to large allele frequency and LD variation among worldwide populations.

Role of microbes in the pathogenesis of neuropsychiatric disorders

Microbes inhabit different anatomical sites of the human body including oral cavity, gut, and skin. A growing literature highlights how microbiome variation is associated with human health and disease. There is strong evidence of bidirectional communication between gut and brain mediated by neurotransmitters and microbial metabolites. Here, we review the potential involvement of microbes residing in the gut and in other body sites in the pathogenesis of eight neuropsychiatric disorders, discussing findings from animal and human studies. The data reported provide a comprehensive overview of the current state of the microbiome research in neuropsychiatry, including hypotheses about the mechanisms underlying the associations reported and the translational potential of probiotics and prebiotics.