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.

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.

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.

Enrichment of self-domestication and neural crest function loci in the heritability of neurodevelopmental disorders.

Self-domestication could contribute to shaping the biology of human brain and consequently the predisposition to neurodevelopmental disorders. Leveraging genome-wide data from the Psychiatric Genomics Consortium, we tested the enrichment of self-domestication and neural crest function loci with respect to the heritability of autism spectrum disorder, schizophrenia (SCZ in East Asian and European ancestries, EAS and EUR, respectively), attention-deficit/hyperactivity disorder, obsessive-compulsive disorder, and Tourette's syndrome (TS). Considering only self-domestication and neural-crest-function annotations in the linkage disequilibrium score regression (LDSC) model, our partitioned heritability analysis revealed statistically significant enrichments across all disorders investigated. The estimates of the heritability enrichments for self-domestication loci were similar across neurodevelopmental disorders, ranging from 0.902 (EAS SCZ, p = 4.55 × 10-20) to 1.577 (TS, p = 5.85 × 10-5). Conversely, a wider spectrum of heritability enrichment estimates was present for neural crest function with the highest enrichment observed for TS (enrichment = 3.453, p = 2.88 × 10-3) and the lowest for EAS SCZ (enrichment = 1.971, p = 3.81 × 10-3). Although these estimates appear to be strong, the enrichments for self-domestication and neural crest function were null once we included additional annotations related to different genomic features. This indicates that the effect of self-domestication on the polygenic architecture of neurodevelopmental disorders is not independent of other functions of human genome.

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 (gcp) linking ALS genetic liability to seven traits. While the genetic component of "self-reported other gastritis including duodenitis" showed a causal effect on ALS (gcp = 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" (gcp = 0.79, p = 5.011 × 10-15) and inverse effects on "other leisure/social group activities" (gcp = 0.66, p = 1 × 10-4) and prospective memory result (gcp = 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.

Genetic liability to suicidal thoughts and behaviors and risk of suicide attempt in US military veterans: moderating effects of cumulative trauma burden.

BACKGROUND: Little is known about environmental factors that may influence associations between genetic liability to suicidality and suicidal behavior. METHODS: This study examined whether a suicidality polygenic risk score (PRS) derived from a large genome-wide association study (N = 122,935) was associated with suicide attempts in a population-based sample of European-American US military veterans (N = 1664; 92.5% male), and whether cumulative lifetime trauma exposure moderated this association. RESULTS: Eighty-five veterans (weighted 6.3%) reported a history of suicide attempt. After adjusting for sociodemographic and psychiatric characteristics, suicidality PRS was associated with lifetime suicide attempt (odds ratio 2.65; 95% CI 1.37-5.11). A significant suicidality PRS-by-trauma exposure interaction emerged, such that veterans with higher levels of suicidality PRS and greater trauma burden had the highest probability of lifetime suicide attempt (16.6%), whereas the probability of attempts was substantially lower among those with high suicidality PRS and low trauma exposure (1.4%). The PRS-by-trauma interaction effect was enriched for genes implicated in cellular and developmental processes, and nervous system development, with variants annotated to the DAB2 and SPNS2 genes, which are implicated in inflammatory processes. Drug repurposing analyses revealed upregulation of suicide gene-sets in the context of medrysone, a drug targeting chronic inflammation, and clofibrate, a triacylglyceride level lowering agent. CONCLUSION: Results suggest that genetic liability to suicidality is associated with increased risk of suicide attempt among veterans, particularly in the presence of high levels of cumulative trauma exposure. Additional research is warranted to investigate whether incorporation of genomic information may improve suicide prediction models.

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"; NTraining = 50%; NTest = 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.

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.

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.

The second-generation antipsychotic drug aripiprazole modulates the serotonergic system in pancreatic islets and induces beta cell dysfunction in female mice.

AIMS/HYPOTHESIS: Second-generation antipsychotic (SGA) drugs have been associated with the development of type 2 diabetes and the metabolic syndrome in patients with schizophrenia. In this study, we aimed to investigate the effects of two different SGA drugs, olanzapine and aripiprazole, on metabolic state and islet function and plasticity. METHODS: We analysed the functional adaptation of beta cells in 12-week-old B6;129 female mice fed an olanzapine- or aripiprazole-supplemented diet (5.5-6.0 mg kg-1 day-1) for 6 months. Glucose and insulin tolerance tests, in vivo glucose-stimulated insulin secretion and indirect calorimetry were performed at the end of the study. The effects of SGAs on beta cell plasticity and islet serotonin levels were assessed by transcriptomic analysis and immunofluorescence. Insulin secretion was assessed by static incubations and Ca2+ fluxes by imaging techniques. RESULTS: Treatment of female mice with olanzapine or aripiprazole for 6 months induced weight gain (p<0.01 and p<0.05, respectively), glucose intolerance (p<0.01) and impaired insulin secretion (p<0.05) vs mice fed a control chow diet. Aripiprazole, but not olanzapine, induced serotonin production in beta cells vs controls, likely by increasing tryptophan hydroxylase 1 (TPH1) expression, and inhibited Ca2+ flux. Of note, aripiprazole increased beta cell size (p<0.05) and mass (p<0.01) vs mice fed a control chow diet, along with activation of mechanistic target of rapamycin complex 1 (mTORC1)/S6 signalling, without preventing beta cell dysfunction. CONCLUSIONS/INTERPRETATION: Both SGAs induced weight gain and beta cell dysfunction, leading to glucose intolerance; however, aripiprazole had a more potent effect in terms of metabolic alterations, which was likely a result of its ability to modulate the serotonergic system. The deleterious metabolic effects of SGAs on islet function should be considered while treating patients as these drugs may increase the risk for development of the metabolic syndrome and diabetes.

HCP5 rs2395029 is a rapid and inexpensive alternative to HLA-B*57:01 genotyping to predict abacavir hypersensitivity reaction in Spain.

Abacavir (ABC) is an HIV nucleotide-analogue reverse transcriptase inhibitor that can produce a severe hypersensitivity reaction (ABC-HSR) in about 5% of the patients. The HLA-B*57:01 allele is associated with the development of ABC-HSR. Therefore, HLA-B*57:01 genotyping is required prior to the prescription of ABC. The technique routinely used in our laboratory is the sequence-specific oligonucleotide probes (SSOP) reverse hybridization method followed by Sanger sequencing. This technique is time-consuming and expensive. The single-nucleotide polymorphism (SNP) HCP5 rs2395029 was described to be in complete linkage disequilibrium with HLA-B*57:01. In this study, we aimed to assess the linkage disequilibrium between HCP5 rs2395029 and HLA-B*57:01 in patients receiving medical assistance at our hospital. We selected 226 HIV-infected patients from our hospital who had been routinely genotyped since 2009 with the SSOP and Sanger sequencing method: 49 HLA-B*57:01 positives and 177 negatives. We genotyped them for HCP5 rs2395019 by real time PCR (qPCR). We exploratory performed two copy number variation assays flanking HCP5 rs2395019 to explore possible deletions that could break the linkage disequilibrium with HLA-B*57:01. The concordance between HLA-B*57:01 and the HCP5 rs2395029 G allele was absolute, with a specificity and sensitivity of 100% (95% confidence interval: 93.0-100.0% and 98.0-100.0%, respectively) and estimated positive and negative predictive values of 84.4% (48.1-93.9%) and 99.9% (99.4-100.0%), respectively. No deletions were found at HCP5 flanking regions. The duration and cost of the SSOP-based method was considerably higher than the SNP-based method. Therefore, the HCP5 rs2395029 genotyping method may be alternatively used in the clinical practice.

Safety and cardiovascular effects of multiple-dose administration of aripiprazole and olanzapine in a randomised clinical trial.

OBJECTIVE: To assess adverse events (AEs) and safety of aripiprazole (ARI) and olanzapine (OLA) treatment. METHODS: Twenty-four healthy volunteers receiving five daily oral doses of 10 mg ARI and 5 mg OLA in a crossover clinical trial were genotyped for 46 polymorphisms in 14 genes by qPCR. Drug plasma concentrations were measured by high-performance liquid chromatography tandem mass spectrometry. Blood pressure (BP) and 12-lead electrocardiogram were measured in supine position. AEs were also recorded. RESULTS: ARI decreased diastolic BP on the first day and decreased QTc on the third and fifth day. OLA had a systolic and diastolic BP, heart rate and QTc lowering effect on the first day. Polymorphisms in ADRA2A, COMT, DRD3 and HTR2A genes were significantly associated to these changes. The most frequent adverse drug reactions (ADRs) to ARI were somnolence, headache, insomnia, dizziness, restlessness, palpitations, akathisia and nausea while were somnolence, dizziness, asthenia, constipation, dry mouth, headache and nausea to OLA. Additionally, HTR2A, HTR2C, DRD2, DRD3, OPRM1, UGT1A1 and CYP1A2 polymorphisms had a role in the development of ADRs. CONCLUSIONS: OLA induced more cardiovascular changes; however, more ADRs were registered to ARI. In addition, some polymorphisms may explain the difference in the incidence of these effects among subjects.

The molecular effect of a polymorphic microRNA binding site of Wolfram syndrome 1 gene in dogs.

BACKGROUND: Although the molecular function of wolframin remains unclear, the lack of this protein is known to cause stress in the endoplasmic reticulum. Some variants in the Wolfram Syndrome 1 gene (WFS1) were associated with various neuropsychiatric disorders in humans, such as aggressiveness, impulsivity and anxiety. RESULTS: Here we present an in silico study predicting a single nucleotide polymorphism (rs852850348) in the canine WFS1 gene which was verified by direct sequencing and was genotyped by a PCR-based technique. We found that the rs852850348 polymorphism is located in a putative microRNA (cfa-miR-8834a and cfa-miR-1838) binding site. Therefore, the molecular effect of allelic variants was studied in a luciferase reporter system that allowed assessing gene expression. We demonstrated that the variant reduced the activity of the reporter protein expression in an allele-specific manner. Additionally, we performed a behavioral experiment and investigated the association with this locus to different performance in this test. Association was found between food possessivity and the studied WFS1 gene polymorphism in the Border collie breed. CONCLUSIONS: Based on our findings, the rs852850348 locus might contribute to the genetic risk of possessivity behavior of dogs in at least one breed and might influence the regulation of wolframin expression.

The effects of aripiprazole and olanzapine on pupillary light reflex and its relationship with pharmacogenetics in a randomized multiple-dose trial.

AIMS: Pupillography is a noninvasive and cost-effective method to determine autonomic nerve activity. Genetic variants in cytochrome P450 (CYP), dopamine receptor (DRD2, DRD3), serotonin receptor (HTR2A, HTR2C) and ATP-binding cassette subfamily B (ABCB1) genes, among others, were previously associated with the pharmacokinetics and pharmacodynamics of antipsychotic drugs. Our aim was to evaluate the effects of aripiprazole and olanzapine on pupillary light reflex related to pharmacogenetics. METHODS: Twenty-four healthy volunteers receiving 5 oral doses of 10 mg aripiprazole and 5 mg olanzapine tablets were genotyped for 46 polymorphisms by quantitative polymerase chain reaction. Pupil examination was performed by automated pupillometry. Aripiprazole, dehydro-aripiprazole and olanzapine plasma concentrations were measured by high-performance liquid chromatography-tandem mass spectrometry. RESULTS: Aripiprazole affected pupil contraction: it caused dilatation after the administration of the first dose, then caused constriction after each dosing. It induced changes in all pupillometric parameters (P < .05). Olanzapine only altered minimum pupil size (P = .046). Polymorphisms in CYP3A, HTR2A, UGT1A1, DRD2 and ABCB1 affected pupil size, the time of onset of constriction, pupil recovery and constriction velocity. Aripiprazole, dehydro-aripiprazole and olanzapine pharmacokinetics were significantly affected by polymorphisms in CYP2D6, CYP3A, CYP1A2, ABCB1 and UGT1A1 genes. CONCLUSIONS: In conclusion, aripiprazole and its main metabolite, dehydro-aripiprazole altered pupil contraction, but olanzapine did not have such an effect. Many polymorphisms may influence pupillometric parameters and several polymorphisms had an effect on aripiprazole, dehydro-aripiprazole and olanzapine pharmacokinetics. Pupillography could be a useful tool for the determination of autonomic nerve activity during antipsychotic treatment.

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.

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.

Distinguishing vulnerability and resilience to posttraumatic stress disorder evaluating traumatic experiences, genetic risk and electronic health records.

What distinguishes vulnerability and resilience to posttraumatic stress disorder (PTSD) remains unclear. Levering traumatic experiences reporting, genetic data, and electronic health records (EHR), we investigated and predicted the clinical comorbidities (co-phenome) of PTSD vulnerability and resilience in the UK Biobank (UKB) and All of Us Research Program (AoU), respectively. In 60,354 trauma-exposed UKB participants, we defined PTSD vulnerability and resilience considering PTSD symptoms, trauma burden, and polygenic risk scores. EHR-based phenome-wide association studies (PheWAS) were conducted to dissect the co-phenomes of PTSD vulnerability and resilience. Significant diagnostic endpoints were applied as weights, yielding a phenotypic risk score (PheRS) to conduct PheWAS of PTSD vulnerability and resilience PheRS in up to 95,761 AoU participants. EHR-based PheWAS revealed three significant phenotypes positively associated with PTSD vulnerability (top association "Sleep disorders") and five outcomes inversely associated with PTSD resilience (top association "Irritable Bowel Syndrome"). In the AoU cohort, PheRS analysis showed a partial inverse relationship between vulnerability and resilience with distinct comorbid associations. While PheRSvulnerability associations were linked to multiple phenotypes, PheRSresilience showed inverse relationships with eye conditions. Our study unveils phenotypic differences in PTSD vulnerability and resilience, highlighting that these concepts are not simply the absence and presence of PTSD.

Psychosocial moderators of polygenic risk scores of inflammatory biomarkers in relation to GrimAge.

GrimAge acceleration has previously predicted age-related morbidities and mortality. In the current study, we sought to examine how GrimAge is associated with genetic predisposition for systemic inflammation and whether psychosocial factors moderate this association. Military veterans from the National Health and Resilience in Veterans study, which surveyed a nationally representative sample of European American male veterans, provided saliva samples for genotyping (N = 1135). We derived polygenic risk scores (PRS) from the UK Biobank as markers of genetic predisposition to inflammation. Results revealed that PRS for three inflammatory PRS markers-HDL (lower), apolipoprotein B (lower), and gamma-glutamyl transferase (higher)-were associated with accelerated GrimAge. Additionally, these PRS interacted with a range of potentially modifiable psychosocial variables, such as exercise and gratitude, previously identified as associated with accelerated GrimAge. Using gene enrichment, we identified anti-inflammatory and antihistamine drugs that perturbate pathways of genes highly represented in the inflammatory PRS, laying the groundwork for future work to evaluate the potential of these drugs in mitigating epigenetic aging.

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.

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.