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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.