Reproducible Genetic Risk Loci for Anxiety: Results From ∼200,000 Participants in the Million Veteran Program

A genome-wide investigation into the underlying genetic architecture of personality traits and overlap with psychopathology

Personality is influenced by both genetic and environmental factors and is associated with other psychiatric traits such as anxiety and depression. The "Big Five" personality traits, which include neuroticism, extraversion, agreeableness, conscientiousness, and openness, are a widely accepted and influential framework for understanding and describing human personality. Of the big five personality traits, neuroticism has most often been the focus of genetic studies and is linked to various mental illnesses including depression, anxiety, and schizophrenia. Our knowledge of the genetic architecture of the other four personality traits is more limited. Utilizing the Million Veteran Program (MVP) cohort we conducted a genome-wide association study (GWAS) in individuals of European and African ancestry. Adding other published data, we performed GWAS meta-analysis for each of the five personality traits with sample sizes ranging from 237,390 to 682,688. We identified 158, 14, 3, 2, and 7 independent genome-wide significant (GWS) loci associated with neuroticism, extraversion, agreeableness, conscientiousness, and openness, respectively. These findings represent 55 novel loci for neuroticism, as well as the first GWS loci discovered for extraversion and agreeableness. Gene-based association testing revealed 254 genes showing significant association with at least one of the five personality traits. Transcriptome-wide and proteome-wide analysis identified altered expression of genes and proteins such as CRHR1, SLC12A5, MAPT, and STX4. Pathway enrichment and drug perturbation analyses identified complex biology underlying human personality traits. We also studied the inter-relationship of personality traits with 1,437 other traits in a phenome-wide genetic correlation analysis, identifying new associations. Mendelian randomization showed positive bidirectional effects between neuroticism and depression and anxiety while a negative bidirectional effect was observed for agreeableness and these psychiatric traits. This study improves our comprehensive understanding of the genetic architecture underlying personality traits and their relationship to other complex human traits.

Apically localized PANX1 impacts neuroepithelial expansion in human cerebral organoids

Dysfunctional paracrine signaling through Pannexin 1 (PANX1) channels is linked to several adult neurological pathologies and emerging evidence suggests that PANX1 plays an important role in human brain development. It remains unclear how early PANX1 influences brain development, or how loss of PANX1 alters the developing human brain. Using a cerebral organoid model of early human brain development, we find that PANX1 is expressed at all stages of organoid development from neural induction through to neuroepithelial expansion and maturation. Interestingly, PANX1 cellular distribution and subcellular localization changes dramatically throughout cerebral organoid development. During neural induction, PANX1 becomes concentrated at the apical membrane domain of neural rosettes where it co-localizes with several apical membrane adhesion molecules. During neuroepithelial expansion, PANX1-/- organoids are significantly smaller than control and exhibit significant gene expression changes related to cell adhesion, WNT signaling and non-coding RNAs. As cerebral organoids mature, PANX1 expression is significantly upregulated and is primarily localized to neuronal populations outside of the ventricular-like zones. Ultimately, PANX1 protein can be detected in all layers of a 21-22 post conception week human fetal cerebral cortex. Together, these results show that PANX1 is dynamically expressed by numerous cell types throughout embryonic and early fetal stages of human corticogenesis and loss of PANX1 compromises neuroepithelial expansion due to dysregulation of cell-cell and cell-matrix adhesion, perturbed intracellular signaling, and changes to gene regulation.

Effects of physical activity and sedentary time on depression, anxiety and well-being: a bidirectional Mendelian randomisation study

BACKGROUND

Mental health conditions represent one of the major groups of non-transmissible diseases. Physical activity (PA) and sedentary time (ST) have been shown to affect mental health outcomes in opposite directions. In this study, we use accelerometery-derived measures of PA and ST from the UK Biobank (UKB) and depression, anxiety and well-being data from the UKB mental health questionnaire as well as published summary statistics to explore the causal associations between these phenotypes.

METHODS

We used MRlap to test if objectively measured PA and ST associate with mental health outcomes using UKB data and summary statistics from published genome-wide association studies. We also tested for bidirectional associations. We performed sex stratified as well as sensitivity analyses.

RESULTS

Genetically instrumented higher PA was associated with lower odds of depression (OR = 0.92; 95% CI: 0.88, 0.97) and depression severity (beta =  - 0.11; 95% CI: - 0.18, - 0.04), Genetically instrumented higher ST was associated higher odds of anxiety (OR = 2.59; 95% CI: 1.10, 4.60). PA was associated with higher well-being (beta = 0.11, 95% CI: 0.04; 0.18) and ST with lower well-being (beta =  - 0.18; 95% CI: - 0.32, - 0.03). Similar findings were observed when stratifying by sex. There was evidence for a bidirectional relationship, with higher genetic liability to depression associated with lower PA (beta =  - 0.25, 95% CI: - 0.42; - 0.08) and higher well-being associated with higher PA (beta = 0.15; 95% CI: 0.05, 0.25).

CONCLUSIONS

We have demonstrated the bidirectional effects of both PA and ST on a range of mental health outcomes using objectively measured predictors and MR methods for causal inference. Our findings support a causal role for PA and ST in the development of mental health problems and in affecting well-being.

Assessing the validity of a self-reported clinical diagnosis of schizophrenia

Background

Diagnoses in psychiatric research can be derived from various sources. This study assesses the validity of a self-reported clinical diagnosis of schizophrenia.

Methods

The study included 3,029 clinically ascertained participants with schizophrenia or psychotic disorders diagnosed by self-report and/or research interview and 1,453 UK Biobank participants with self-report and/or medical record diagnosis of schizophrenia or schizoaffective disorder depressed-type (SA-D). We assessed positive predictive values (PPV) of self-reported clinical diagnoses against research interview and medical record diagnoses. We compared polygenic risk scores (PRS) and phenotypes across diagnostic groups, and compared the variance explained by schizophrenia PRS to samples in the Psychiatric Genomics Consortium (PGC).

Results

In the clinically ascertained sample, the PPV of self-reported schizophrenia to a research diagnosis of schizophrenia was 0.70, which increased to 0.81 when benchmarked against schizophrenia or SA-D. In UK Biobank, the PPV of self-reported schizophrenia to a medical record diagnosis was 0.74. Compared to self-report participants, those with a research diagnosis were younger and more likely to have a high school qualification (clinically ascertained sample) and those with a medical record diagnosis were less likely to be employed or have a high school qualification (UK Biobank). Schizophrenia PRS did not differ between participants that had a diagnosis from self-report, research diagnosis or medical record diagnosis. Polygenic liability r2, for all diagnosis definitions, fell within the distribution of PGC schizophrenia cohorts.

Conclusions

Self-report measures of schizophrenia are justified in research to maximise sample size and representativeness, although within sample validation of diagnoses is recommended.

Multi-ancestry genome-wide association study of cannabis use disorder yields insight into disease biology and public health implications

As recreational use of cannabis is being decriminalized in many places and medical use widely sanctioned, there are growing concerns about increases in cannabis use disorder (CanUD), which is associated with numerous medical comorbidities. Here we performed a genome-wide association study of CanUD in the Million Veteran Program (MVP), followed by meta-analysis in 1,054,365 individuals (ncases = 64,314) from four broad ancestries designated by the reference panel used for assignment (European n = 886,025, African n = 123,208, admixed American n = 38,289 and East Asian n = 6,843). Population-specific methods were applied to calculate single nucleotide polymorphism-based heritability within each ancestry. Statistically significant single nucleotide polymorphism-based heritability for CanUD was observed in all but the smallest population (East Asian). We discovered genome-wide significant loci unique to each ancestry: 22 in European, 2 each in African and East Asian, and 1 in admixed American ancestries. A genetically informed causal relationship analysis indicated a possible effect of genetic liability for CanUD on lung cancer risk, suggesting potential unanticipated future medical and psychiatric public health consequences that require further study to disentangle from other known risk factors such as cigarette smoking.

Dimensional and transdiagnostic phenotypes in psychiatric genome-wide association studies

Genome-wide association studies (GWAS) provide biological insights into disease onset and progression and have potential to produce clinically useful biomarkers. A growing body of GWAS focuses on quantitative and transdiagnostic phenotypic targets, such as symptom severity or biological markers, to enhance gene discovery and the translational utility of genetic findings. The current review discusses such phenotypic approaches in GWAS across major psychiatric disorders. We identify themes and recommendations that emerge from the literature to date, including issues of sample size, reliability, convergent validity, sources of phenotypic information, phenotypes based on biological and behavioral markers such as neuroimaging and chronotype, and longitudinal phenotypes. We also discuss insights from multi-trait methods such as genomic structural equation modelling. These provide insight into how hierarchical 'splitting' and 'lumping' approaches can be applied to both diagnostic and dimensional phenotypes to model clinical heterogeneity and comorbidity. Overall, dimensional and transdiagnostic phenotypes have enhanced gene discovery in many psychiatric conditions and promises to yield fruitful GWAS targets in the years to come.

Anxiety and post-traumatic stress disorders in patients with chronic respiratory diseases

PURPOSE OF REVIEW

This brief critical review evaluates recent research on the impact of anxiety disorders and post-traumatic stress disorders (PTSD) on patients with chronic obstructive pulmonary disease (COPD) and asthma. In these patients, untreated anxiety and PTSD exacerbate acute symptoms, increase disability and impair quality of life. Therefore, effective interventions are also reviewed.

RECENT FINDINGS

Anxiety disorders and PTSD are prevalent in COPD and asthma, worsen symptoms in acute exacerbations, and increase morbidity and healthcare utilization. Pulmonary rehabilitation (PR), cognitive behavioural therapy and pharmacological therapy are each effective in COPD patients with anxiety or PTSD. However, very little work examined therapy for combined anxiety and PTSD in patients with either COPD or asthma.

SUMMARY

Despite the high prevalence of anxiety disorder or PTSD in patients with chronic respiratory diseases, a paucity of literature remains demonstrating the efficacy of pharmacological therapies for these conditions. This review highlights the promising benefits of PR on anxiety, but prospective trials are needed to demonstrate the efficacy of interventions with PTSD alone, or with concomitant anxiety.

Testing associations between human anxiety and genes previously implicated by mouse anxiety models

Anxiety disorders are common and can be debilitating, with effective treatments remaining hampered by an incomplete understanding of the underlying genetic etiology. Improvements have been made in understanding the genetic influences on mouse behavioral models of anxiety, yet it is unclear the extent to which genes identified in these experimental systems contribute to genetic variation in human anxiety phenotypes. Leveraging new and existing large-scale human genome-wide association studies, we tested whether sets of genes previously identified in mouse anxiety-like behavior studies contribute to a range of human anxiety disorders. When tested as individual genes, 13 mouse-identified genes were associated with human anxiety phenotypes, suggesting an overlap of individual genes contributing to both mouse models of anxiety-like behaviors and human anxiety traits. When genes were tested as sets, we did identify 14 significant associations between mouse gene sets and human anxiety, but the majority of gene sets showed no significant association with human anxiety phenotypes. These few significant associations indicate a need to identify and develop more translatable mouse models by identifying sets of genes that "match" between model systems and specific human phenotypes of interest. We suggest that continuing to develop improved behavioral paradigms and finer-scale experimental data, for instance from individual neuronal subtypes or cell-type-specific expression data, is likely to improve our understanding of the genetic etiology and underlying functional changes in anxiety disorders.

Investigating the role of the relaxin-3/RXFP3 system in neuropsychiatric disorders and metabolic phenotypes: A candidate gene approach

The relaxin-3/RXFP3 system has been implicated in the modulation of depressive- and anxiety-like behaviour in the animal literature; however, there is a lack of human studies investigating this signalling system. We seek to bridge this gap by leveraging the large UK Biobank study to retrospectively assess genetic risk variants linked with this neuropeptidergic system. Specifically, we conducted a candidate gene study in the UK Biobank to test for potential associations between a set of functional, candidate single nucleotide polymorphisms (SNPs) pertinent to relaxin-3 signalling, determined using in silico tools, and several outcomes, including depression, atypical depression, anxiety and metabolic syndrome. For each outcome, we used several rigorously defined phenotypes, culminating in subsample sizes ranging from 85,881 to 386,769 participants. Across all outcomes, there were no associations between any candidate SNP and any outcome phenotype, following corrections for multiple testing burden. Regression models comprising several SNPs per relevant candidate gene as exploratory variables further exhibited no prediction of outcome. Our findings corroborate conclusions from previous literature about the limitations of candidate gene approaches, even when based on firm biological hypotheses, in the domain of genetic research for neuropsychiatric disorders.

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.

Assessment of the Genetic Characteristics of a Generation Born during a Long-Term Socioeconomic Crisis

BACKGROUND

A socioeconomic crisis in Russia lasted from 1991 to 1998 and was accompanied by a sharp drop in the birth rate. The main factor that influenced the refusal to have children during this period is thought to be prolonged social stress.

METHODS

comparing frequencies of common gene variants associated with stress-induced diseases among generations born before, after, and during this crisis may show which genes may be preferred under the pressure of natural selection during periods of increased social stress in urban populations.

RESULTS

In the "crisis" group, a statistically significant difference from the other two groups was found in rs6557168 frequency (p = 0.001); rs4522666 was not in the Hardy-Weinberg equilibrium in this group, although its frequency did not show a significant difference from the other groups (p = 0.118). Frequencies of VNTRs in SLC6A3 and MAOA as well as common variants rs17689918 in CRHR1, rs1360780 in FKBP5, rs53576 in OXTR, rs12720071 and rs806377 in CNR1, rs4311 in ACE, rs1800497 in ANKK1, and rs7412 and rs429358 in APOE did not differ among the groups.

CONCLUSIONS

a generation born during a period of prolonged destructive events may differ from the rest of the gene pool of the population in some variants associated with personality traits or stress-related disorders.

Sex differences in avoidance behavior and cued threat memory dynamics in mice: Interactions between estrous cycle and genetic background

Anxiety disorders are the most prevalent mental illnesses worldwide, exhibit high heritability, and affect twice as many women as men. To evaluate potential interactions between genetic background and cycling ovarian hormones on sex differences in susceptibility to negative valence behaviors relevant to anxiety disorders, we assayed avoidance behavior and cued threat memory dynamics in gonadally-intact adult male and female mice across four common inbred mouse strains: C57Bl/6J, 129S1/SVlmJ, DBA/2J, and BALB/cJ. Independent of sex, C57Bl/6J mice exhibited low avoidance but high threat memory, 129S1/SvlmJ mice high avoidance and high threat memory, DBA/2J mice low avoidance and low threat memory, and BALB/cJ mice high avoidance but low threat memory. Within-strain comparisons revealed reduced avoidance behavior in the high hormone phase of the estrous cycle (proestrus) compared to all other estrous phases in all strains except DBA/2J, which did not exhibit cycle-dependent behavioral fluctuations. Robust and opposing sex differences in threat conditioning and extinction training were found in the C57Bl/6J and 129S1/SvlmJ lines, whereas no sex differences were observed in the DBA/2J or BALB/cJ lines. C57Bl/6J males exhibited enhanced acute threat memory, whereas 129S1/SvlmJ females exhibited enhanced sustained threat memory, compared to their sex-matched littermates. These effects were not mediated by estrous cycle stage or sex differences in active versus passive defensive behavioral responses. Our data demonstrate that core features of behavioral endophenotypes relevant to anxiety disorders, such as avoidance and threat memory, are genetically driven yet dissociable and can be influenced further by cycling ovarian hormones.

Examining the relationship between genetic risk for depression and youth episodic stress exposure

BACKGROUND

Offspring of depressed mothers have elevated risk of developing depression because they are exposed to greater stress. While generally assumed that youth's increased exposure to stress is due to the environmental effects of living with a depressed parent, youth's genes may influence stress exposure through gene-environment correlations (rGEs). To understand the relationship between risk for depression and stress, we examined the effects of polygenic risk for depression on youth stress exposure.

METHODS

We examined the relations of a polygenic risk score (PRS) for depression (DEP-PRS), as well as PRSs for 5 other disorders, with youth stress exposure. Data were from a longitudinal study of a community sample of youth and their parents (n = 377) focusing on data collected at youth's aged 12 and 15 assessments.

RESULTS

Elevated youth DEP-PRS was robustly associated with increased dependent stress, particularly interpersonal events. Exploratory analyses indicated that findings were driven by major stress and were not moderated by maternal nor paternal history of depression, and of the 5 additional PRSs tested, only elevated genetic liability for bipolar I was associated with increased dependent stress-particularly non-interpersonal events.

LIMITATIONS

Like other PRS studies, we focused on those of European ancestry thus, generalizability of findings is limited.

CONCLUSION

Polygenic risk contributes to youth experiencing stressful life events which are dependent on their behavior. This rGE appears to be specific to genetic risk for mood disorders.

The genetic contribution to the comorbidity of depression and anxiety: a multi-site electronic health records study of almost 178 000 people

BACKGROUND

Depression and anxiety are common and highly comorbid, and their comorbidity is associated with poorer outcomes posing clinical and public health concerns. We evaluated the polygenic contribution to comorbid depression and anxiety, and to each in isolation.

METHODS

Diagnostic codes were extracted from electronic health records for four biobanks [N = 177 865 including 138 632 European (77.9%), 25 612 African (14.4%), and 13 621 Hispanic (7.7%) ancestry participants]. The outcome was a four-level variable representing the depression/anxiety diagnosis group: neither, depression-only, anxiety-only, and comorbid. Multinomial regression was used to test for association of depression and anxiety polygenic risk scores (PRSs) with the outcome while adjusting for principal components of ancestry.

RESULTS

In total, 132 960 patients had neither diagnosis (74.8%), 16 092 depression-only (9.0%), 13 098 anxiety-only (7.4%), and 16 584 comorbid (9.3%). In the European meta-analysis across biobanks, both PRSs were higher in each diagnosis group compared to controls. Notably, depression-PRS (OR 1.20 per s.d. increase in PRS; 95% CI 1.18-1.23) and anxiety-PRS (OR 1.07; 95% CI 1.05-1.09) had the largest effect when the comorbid group was compared with controls. Furthermore, the depression-PRS was significantly higher in the comorbid group than the depression-only group (OR 1.09; 95% CI 1.06-1.12) and the anxiety-only group (OR 1.15; 95% CI 1.11-1.19) and was significantly higher in the depression-only group than the anxiety-only group (OR 1.06; 95% CI 1.02-1.09), showing a genetic risk gradient across the conditions and the comorbidity.

CONCLUSIONS

This study suggests that depression and anxiety have partially independent genetic liabilities and the genetic vulnerabilities to depression and anxiety make distinct contributions to comorbid depression and anxiety.

Adverse Childhood Events, Mood and Anxiety Disorders, and Substance Dependence: Gene x Environment Effects and Moderated Mediation

Background

Adverse childhood events (ACEs) contribute to the development of mood and anxiety disorders and substance dependence. However, the extent to which these effects are direct or indirect and whether genetic risk moderates them is unclear.

Methods

We examined associations among ACEs, mood/anxiety disorders, and substance dependence in 12,668 individuals (44.9% female, 42.5% African American/Black, 42.1% European American/White). We generated latent variables for each phenotype and modeled direct and indirect effects of ACEs on substance dependence, mediated by mood/anxiety disorders (forward or "self-medication" model) and of ACEs on mood/anxiety disorders, mediated by substance dependence (reverse or "substance-induced" model). In a sub-sample, we also generated polygenic scores for substance dependence and mood/anxiety disorder factors, which we tested as moderators in the mediation models.

Results

Although there were significant indirect effects in both directions, mediation by mood/anxiety disorders (forward model) was greater than by substance dependence (reverse model). Greater genetic risk for substance dependence was associated with a weaker direct effect of ACEs on substance dependence in both the African- and European-ancestry groups (i.e., gene-environment interaction) and a weaker indirect effect in European-ancestry individuals (i.e., moderated mediation).

Conclusion

We found greater evidence that substance dependence results from self-medication of mood/anxiety disorders than that mood/anxiety disorders are substance induced. Among individuals at higher genetic risk for substance dependence who are more likely to develop a dependence diagnosis, ACEs exert less of an effect in promoting that outcome. Following exposure to ACEs, multiple pathways lead to mood/anxiety disorders and substance dependence. Specification of these pathways could inform individually targeted prevention and treatment approaches.

Adverse Childhood Events, Mood and Anxiety Disorders, and Substance Dependence: Gene X Environment Effects and Moderated Mediation

Background

Adverse childhood events (ACEs) contribute to the development of mood and anxiety disorders and substance dependence. However, the extent to which these effects are direct or indirect and whether genetic risk moderates them is unclear.

Methods

We examined associations among ACEs, mood/anxiety disorders, and substance dependence in 12,668 individuals (44.9% female, 42.5% African American/Black, 42.1% European American/White). We generated latent variables for each phenotype and modeled direct and indirect effects of ACEs on substance dependence, mediated by mood/anxiety disorders (forward or "self-medication" model) and of ACEs on mood/anxiety disorders, mediated by substance dependence (reverse or "substance-induced" model). In a sub-sample, we also generated polygenic scores for substance dependence and mood/anxiety disorder factors, which we tested as moderators in the mediation models.

Results

Although there were significant indirect effects in both directions, mediation by mood/anxiety disorders (forward model) was greater than by substance dependence (reverse model). Greater genetic risk for substance dependence was associated with a weaker direct effect of ACEs on substance dependence in both the African- and European-ancestry groups (i.e., gene-environment interaction) and a weaker indirect effect in European-ancestry individuals (i.e., moderated mediation).

Conclusion

We found greater evidence that substance dependence results from self-medication of mood/anxiety disorders than that mood/anxiety disorders are substance induced. Among individuals at higher genetic risk for substance dependence who are more likely to develop a dependence diagnosis, ACEs exert less of an effect in promoting that outcome. Following exposure to ACEs, multiple pathways lead to mood/anxiety disorders and substance dependence. Specification of these pathways could inform individually targeted prevention and treatment approaches.

Impact of Copy Number Variants and Polygenic Risk Scores on Psychopathology in the UK Biobank

BACKGROUND

Our understanding of the impact of copy number variants (CNVs) on psychopathology and their joint influence with polygenic risk scores (PRSs) remains limited.

METHODS

The UK Biobank recruited 502,534 individuals ages 37 to 73 years living in the United Kingdom between 2006 and 2010. After quality control, genotype data from 459,855 individuals were available for CNV calling. A total of 61 commonly studied recurrent neuropsychiatric CNVs were selected for analyses and examined individually and in aggregate (any CNV, deletion, or duplication). CNV risk scores were used to quantify intolerance of CNVs to haploinsufficiency. Major depressive disorder and generalized anxiety disorder PRSs were generated for White British individuals (N = 408,870). Mood/anxiety factor scores were generated using item-level questionnaire data (N = 501,289).

RESULTS

CNV carriers showed higher mood/anxiety scores than noncarriers, with the largest effects seen for intolerant deletions. A total of 11 individual deletions and 8 duplications were associated with higher mood/anxiety. Carriers of the 9p24.3 (DMRT1) duplication showed lower mood/anxiety. Associations remained significant for most CNVs when excluding individuals with psychiatric diagnoses. Nominally significant CNV × PRS interactions provided preliminary evidence that associations between select individual CNVs, but not CNVs in aggregate, and mood/anxiety may be modulated by PRSs.

CONCLUSIONS

CNVs associated with risk for psychiatric disorders showed small to large effects on dimensional mood/anxiety scores in a general population cohort, even when excluding individuals with psychiatric diagnoses. CNV × PRS interactions showed that associations between select CNVs and mood/anxiety may be modulated by PRSs.

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.

Genome-wide association studies and cross-population meta-analyses investigating short and long sleep duration

Sleep duration has been linked to a wide range of negative health outcomes and to reduced life expectancy. We present genome-wide association studies of short ( ≤ 5 h) and long ( ≥ 10 h) sleep duration in adults of European (N = 445,966), African (N = 27,785), East Asian (N = 3141), and admixed-American (N = 16,250) ancestry from UK Biobank and the Million Veteran Programme. In a cross-population meta-analysis, we identify 84 independent loci for short sleep and 1 for long sleep. We estimate SNP-based heritability for both sleep traits in each ancestry based on population derived linkage disequilibrium (LD) scores using cov-LDSC. We identify positive genetic correlation between short and long sleep traits (rg = 0.16 ± 0.04; p = 0.0002), as well as similar patterns of genetic correlation with other psychiatric and cardiometabolic phenotypes. Mendelian randomisation reveals a directional causal relationship between short sleep and depression, and a bidirectional causal relationship between long sleep and depression.

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 and cell-specific gene expression in corticolimbic brain regions associated with psychiatric disorders revealed by bulk and single-nuclei RNA sequencing

BACKGROUND

There are sex-specific differences in the prevalence, symptomology and course of psychiatric disorders. However, preclinical models have primarily used males, such that the molecular mechanisms underlying sex-specific differences in psychiatric disorders are not well established.

METHODS

In this study, we compared transcriptome-wide gene expression profiles in male and female rats within the corticolimbic system, including the cingulate cortex, nucleus accumbens medial shell (NAcS), ventral dentate gyrus and the basolateral amygdala (n = 22-24 per group/region).

FINDINGS

We found over 3000 differentially expressed genes (DEGs) in the NAcS between males and females. Of these DEGs in the NAcS, 303 showed sex-dependent conservation DEGs in humans and were significantly enriched for gene ontology terms related to blood vessel morphogenesis and regulation of cell migration. Single nuclei RNA sequencing in the NAcS of male and female rats identified widespread sex-dependent expression, with genes upregulated in females showing a notable enrichment for synaptic function. Female upregulated genes in astrocytes, Drd3+MSNs and oligodendrocyte were also enriched in several psychiatric genome-wide association studies (GWAS).

INTERPRETATION

Our data provide comprehensive evidence of sex- and cell-specific molecular profiles in the NAcS. Importantly these differences associate with anxiety, bipolar disorder, schizophrenia, and cross-disorder, suggesting an intrinsic molecular basis for sex-based differences in psychiatric disorders that strongly implicates the NAcS.

FUNDING

This work was supported by funding from the Hope for Depression Research Foundation (MJM).

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.

Identification of drug repurposing candidates for the treatment of anxiety: A genetic approach

Anxiety disorders are a group of prevalent and heritable neuropsychiatric diseases. We previously conducted a genome-wide association study (GWAS) which identified genomic loci associated with anxiety; however, the biological consequences underlying the genetic associations are largely unknown. Integrating GWAS and functional genomic data may improve our understanding of the genetic effects on intermediate molecular phenotypes such as gene expression. This can provide an opportunity for the discovery of drug targets for anxiety via drug repurposing. We used the GWAS summary statistics to determine putative causal genes for anxiety using MAGMA and colocalization analyses. A transcriptome-wide association study was conducted to identify genes with differential genetically regulated levels of gene expression in human brain tissue. The genes were integrated with a large drug-gene expression database (Connectivity Map), discovering compounds that are predicted to "normalise" anxiety-associated expression changes. The study identified 64 putative causal genes associated with anxiety (35 genes upregulated; 29 genes downregulated). Drug mechanisms adrenergic receptor agonists, sigma receptor agonists, and glutamate receptor agonists gene targets were enriched in anxiety-associated genetic signal and exhibited an opposing effect on the anxiety-associated gene expression signature. The significance of the project demonstrated genetic links for novel drug candidates to potentially advance anxiety therapeutics.

Methylome-wide association study of anxiety disorders

Anxiety Disorders (ANX) such as panic disorder, generalized anxiety disorder, and phobias, are highly prevalent conditions that are moderately heritable. Evidence suggests that DNA methylation may play a role, as it is involved in critical adaptations to changing environments. Applying an enrichment-based sequencing approach covering nearly 28 million autosomal CpG sites, we conducted a methylome-wide association study (MWAS) of lifetime ANX in 1132 participants (618 cases/514 controls) from the Netherlands Study of Depression and Anxiety. Using epigenomic deconvolution, we performed MWAS for the main cell types in blood: granulocytes, T-cells, B-cells and monocytes. Cell-type specific analyses identified 280 and 82 methylome-wide significant associations (q-value < 0.1) in monocytes and granulocytes, respectively. Our top finding in monocytes was located in ZNF823 on chromosome 19 (p = 1.38 × 10-10) previously associated with schizophrenia. We observed significant overlap (p < 1 × 10-06) with the same direction of effect in monocytes (210 sites), T-cells (135 sites), and B-cells (727 sites) between this Discovery MWAS signal and a comparable replication dataset from the Great Smoky Mountains Study (N = 433). Overlapping Discovery-Replication MWAS signal was enriched for findings from published GWAS of ANX, major depression, and post-traumatic stress disorder. In monocytes, two specific sites in the FZR1 gene showed significant replication after Bonferroni correction with an additional 15 nominally replicated sites in monocytes and 4 in T-cells. FZR1 regulates neurogenesis in the hippocampus, and its knockout leads to impairments in associative fear memory and long-term potentiation in mice. In the largest and most extensive methylome-wide study of ANX, we identified replicable methylation sites located in genes of potential relevance for brain mechanisms of psychiatric conditions.

The nature and neurobiology of fear and anxiety: State of the science and opportunities for accelerating discovery

Fear and anxiety play a central role in mammalian life, and there is considerable interest in clarifying their nature, identifying their biological underpinnings, and determining their consequences for health and disease. Here we provide a roundtable discussion on the nature and biological bases of fear- and anxiety-related states, traits, and disorders. The discussants include scientists familiar with a wide variety of populations and a broad spectrum of techniques. The goal of the roundtable was to take stock of the state of the science and provide a roadmap to the next generation of fear and anxiety research. Much of the discussion centered on the key challenges facing the field, the most fruitful avenues for future research, and emerging opportunities for accelerating discovery, with implications for scientists, funders, and other stakeholders. Understanding fear and anxiety is a matter of practical importance. Anxiety disorders are a leading burden on public health and existing treatments are far from curative, underscoring the urgency of developing a deeper understanding of the factors governing threat-related emotions.

Shared genetic architecture between irritable bowel syndrome and psychiatric disorders reveals molecular pathways of the gut-brain axis

BACKGROUND

Irritable bowel syndrome (IBS) often co-occurs with psychiatric and gastrointestinal disorders. A recent genome-wide association study (GWAS) identified several genetic risk variants for IBS. However, most of the heritability remains unidentified, and the genetic overlap with psychiatric and somatic disorders is not quantified beyond genome-wide genetic correlations. Here, we characterize the genetic architecture of IBS, further, investigate its genetic overlap with psychiatric and gastrointestinal phenotypes, and identify novel genomic risk loci.

METHODS

Using GWAS summary statistics of IBS (53,400 cases and 433,201 controls), and psychiatric and gastrointestinal phenotypes, we performed bivariate casual mixture model analysis to characterize the genetic architecture and genetic overlap between these phenotypes. We leveraged identified genetic overlap to boost the discovery of genomic loci associated with IBS, and to identify specific shared loci associated with both IBS and psychiatric and gastrointestinal phenotypes, using the conditional/conjunctional false discovery rate (condFDR/conjFDR) framework. We used functional mapping and gene annotation (FUMA) for functional analyses.

RESULTS

IBS was highly polygenic with 12k trait-influencing variants. We found extensive polygenic overlap between IBS and psychiatric disorders and to a lesser extent with gastrointestinal diseases. We identified 132 independent IBS-associated loci (condFDR < 0.05) by conditioning on psychiatric disorders (n = 127) and gastrointestinal diseases (n = 24). Using conjFDR, 70 unique loci were shared between IBS and psychiatric disorders. Functional analyses of shared loci revealed enrichment for biological pathways of the nervous and immune systems. Genetic correlations and shared loci between psychiatric disorders and IBS subtypes were different.

CONCLUSIONS

We found extensive polygenic overlap of IBS and psychiatric and gastrointestinal phenotypes beyond what was revealed with genetic correlations. Leveraging the overlap, we discovered genetic loci associated with IBS which implicate a wide range of biological pathways beyond the gut-brain axis. Genetic differences may underlie the clinical subtype of IBS. These results increase our understanding of the pathophysiology of IBS which may form the basis for the development of individualized interventions.

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.

Depression pathophysiology, risk prediction of recurrence and comorbid psychiatric disorders using genome-wide analyses

Depression is a common psychiatric disorder and a leading cause of disability worldwide. Here we conducted a genome-wide association study meta-analysis of six datasets, including >1.3 million individuals (371,184 with depression) and identified 243 risk loci. Overall, 64 loci were new, including genes encoding glutamate and GABA receptors, which are targets for antidepressant drugs. Intersection with functional genomics data prioritized likely causal genes and revealed new enrichment of prenatal GABAergic neurons, astrocytes and oligodendrocyte lineages. We found depression to be highly polygenic, with ~11,700 variants explaining 90% of the single-nucleotide polymorphism heritability, estimating that >95% of risk variants for other psychiatric disorders (anxiety, schizophrenia, bipolar disorder and attention deficit hyperactivity disorder) were influencing depression risk when both concordant and discordant variants were considered, and nearly all depression risk variants influenced educational attainment. Additionally, depression genetic risk was associated with impaired complex cognition domains. We dissected the genetic and clinical heterogeneity, revealing distinct polygenic architectures across subgroups of depression and demonstrating significantly increased absolute risks for recurrence and psychiatric comorbidity among cases of depression with the highest polygenic burden, with considerable sex differences. The risks were up to 5- and 32-fold higher than cases with the lowest polygenic burden and the background population, respectively. These results deepen the understanding of the biology underlying depression, its disease progression and inform precision medicine approaches to treatment.

Towards precision medicine for anxiety disorders: objective assessment, risk prediction, pharmacogenomics, and repurposed drugs

Anxiety disorders are increasingly prevalent, affect people's ability to do things, and decrease quality of life. Due to lack of objective tests, they are underdiagnosed and sub-optimally treated, resulting in adverse life events and/or addictions. We endeavored to discover blood biomarkers for anxiety, using a four-step approach. First, we used a longitudinal within-subject design in individuals with psychiatric disorders to discover blood gene expression changes between self-reported low anxiety and high anxiety states. Second, we prioritized the list of candidate biomarkers with a Convergent Functional Genomics approach using other evidence in the field. Third, we validated our top biomarkers from discovery and prioritization in an independent cohort of psychiatric subjects with clinically severe anxiety. Fourth, we tested these candidate biomarkers for clinical utility, i.e. ability to predict anxiety severity state, and future clinical worsening (hospitalizations with anxiety as a contributory cause), in another independent cohort of psychiatric subjects. We showed increased accuracy of individual biomarkers with a personalized approach, by gender and diagnosis, particularly in women. The biomarkers with the best overall evidence were GAD1, NTRK3, ADRA2A, FZD10, GRK4, and SLC6A4. Finally, we identified which of our biomarkers are targets of existing drugs (such as a valproate, omega-3 fatty acids, fluoxetine, lithium, sertraline, benzodiazepines, and ketamine), and thus can be used to match patients to medications and measure response to treatment. We also used our biomarker gene expression signature to identify drugs that could be repurposed for treating anxiety, such as estradiol, pirenperone, loperamide, and disopyramide. Given the detrimental impact of untreated anxiety, the current lack of objective measures to guide treatment, and the addiction potential of existing benzodiazepines-based anxiety medications, there is a urgent need for more precise and personalized approaches like the one we developed.

Candidate biomarkers in psychiatric disorders: state of the field

The field of psychiatry is hampered by a lack of robust, reliable and valid biomarkers that can aid in objectively diagnosing patients and providing individualized treatment recommendations. Here we review and critically evaluate the evidence for the most promising biomarkers in the psychiatric neuroscience literature for autism spectrum disorder, schizophrenia, anxiety disorders and post-traumatic stress disorder, major depression and bipolar disorder, and substance use disorders. Candidate biomarkers reviewed include various neuroimaging, genetic, molecular and peripheral assays, for the purposes of determining susceptibility or presence of illness, and predicting treatment response or safety. This review highlights a critical gap in the biomarker validation process. An enormous societal investment over the past 50 years has identified numerous candidate biomarkers. However, to date, the overwhelming majority of these measures have not been proven sufficiently reliable, valid and useful to be adopted clinically. It is time to consider whether strategic investments might break this impasse, focusing on a limited number of promising candidates to advance through a process of definitive testing for a specific indication. Some promising candidates for definitive testing include the N170 signal, an event-related brain potential measured using electroencephalography, for subgroup identification within autism spectrum disorder; striatal resting-state functional magnetic resonance imaging (fMRI) measures, such as the striatal connectivity index (SCI) and the functional striatal abnormalities (FSA) index, for prediction of treatment response in schizophrenia; error-related negativity (ERN), an electrophysiological index, for prediction of first onset of generalized anxiety disorder, and resting-state and structural brain connectomic measures for prediction of treatment response in social anxiety disorder. Alternate forms of classification may be useful for conceptualizing and testing potential biomarkers. Collaborative efforts allowing the inclusion of biosystems beyond genetics and neuroimaging are needed, and online remote acquisition of selected measures in a naturalistic setting using mobile health tools may significantly advance the field. Setting specific benchmarks for well-defined target application, along with development of appropriate funding and partnership mechanisms, would also be crucial. Finally, it should never be forgotten that, for a biomarker to be actionable, it will need to be clinically predictive at the individual level and viable in clinical settings.

Examining the shared etiology of psychopathology with genome-wide association studies

Genome-wide association studies (GWASs) have ushered in a new era of reproducible discovery in psychiatric genetics. The field has now identified hundreds of common genetic variants that are associated with mental disorders, and many of them influence more than one disorder. By advancing the understanding of causal biology underlying psychopathology, GWAS results are poised to inform the development of novel therapeutics, stratification of at-risk patients, and perhaps even the revision of top-down classification systems in psychiatry. Here, we provide a concise review of GWAS findings with an emphasis on findings that have elucidated the shared genetic etiology of psychopathology, summarizing insights at three levels of analysis: 1) genome-wide architecture; 2) networks, pathways, and gene sets; and 3) individual variants/genes. Three themes emerge from these efforts. First, all psychiatric phenotypes are heritable, highly polygenic, and influenced by many pleiotropic variants with incomplete penetrance. Second, GWAS results highlight the broad etiological roles of neuronal biology, system-wide effects over localized effects, and early neurodevelopment as a critical period. Third, many loci that are robustly associated with multiple forms of psychopathology harbor genes that are involved in synaptic structure and function. Finally, we conclude our review by discussing the implications that GWAS results hold for the field of psychiatry, as well as expected challenges and future directions in the next stage of psychiatric genetics.

Genetic insights into the neurobiology of anxiety

Anxiety and fear are evolutionarily conserved emotions that increase the likelihood of an organism surviving threatening situations. Anxiety and vigilance states are regulated by neural networks involving multiple brain regions. In anxiety disorders, this intricate regulatory system is disturbed, leading to excessive or prolonged anxiety or fear. Anxiety disorders have both genetic and environmental risk factors. Genetic research has the potential to identify specific genetic variants causally associated with specific phenotypes. In recent decades, genome-wide association studies (GWASs) have revealed variants predisposing to neuropsychiatric disorders, suggesting novel neurobiological pathways in the etiology of these disorders. Here, we review recent human GWASs of anxiety disorders, and genetic studies of anxiety-like behavior in rodent models. These studies are paving the way for a better understanding of the neurobiological mechanisms underlying anxiety disorders.

Pharmacotherapy exposure as a marker of disease complexity in bipolar disorder: Associations with clinical & genetic risk factors

Individuals with bipolar disorder (BD) require chronic pharmacotherapy, typically including medication switches or polypharmacy due to persisting symptoms or intolerable side effects. Here, we quantified pharmacotherapy exposure (PE) of Mayo Clinic BD Biobank participants using the number of cross-sectional (at enrollment) and lifetime BD-specific medications and medication classes, to understand the relationship between PE and markers of disease severity or treatment failure, psychiatric comorbidities, and polygenic risk scores (PRS) for six major psychiatric disorders. Being female (p < 0.05), older (p < 0.01), having history of suicide attempts (p < 0.0001), and comorbid attention-deficit/hyperactivity disorder (p < 0.05) or generalized anxiety disorder (p < 0.05) were uniformly associated with higher PE. Lifetime exposure to unique medication classes among participants with BD-I was significantly lower than for those with schizoaffective disorder (estimate = -2.1, p < 0.0001) while significantly higher than for those with BD-II (estimate = 0.5, p < 0.01). Further, higher PRS for schizophrenia (SCZ) and anxiety resulted in greater lifetime medication counts (p < 0.01), both driven by antipsychotic (p < 0.001) and anxiolytic use (p < 0.05). Our results provide initial evidence of the utility of PE as a measure of disease complexity or treatment resistance, and that PE may be predicted by higher genetic risk for SCZ and anxiety.

Loss of Tyro3 causes anxiety-relevant behavioural changes in female mice

White-matter brain abnormalities have been found across a variety of psychiatric disorders. The extent of white matter pathology is proposed to be predictive of the severity of anxiety disorders. However, it is still unknown whether disruptions of white matter integrity precede, and are sufficient to give rise to, the behavioural symptoms. Interestingly, mood disturbances feature prominently in central demyelinating diseases such as multiple sclerosis. It is unclear whether the greater frequency of neuropsychiatric symptoms is linked to underlying neuropathology. In this study, we characterised male and female Tyro3 knockout (KO) mice using a variety of behavioural paradigms. Anxiety-related behaviours were assessed with the elevated-plus maze and light-dark box. Fear memory processing was assessed using fear conditioning and extinction paradigms. Finally, we assessed immobility time in the Porsolt swim test as a measure of depression-related behavioural despair. Surprisingly, loss of Tyro3 did not lead to manifestation of major shifts in baseline behaviour. We noted significant differences in habituation to novel environments and post-conditioning freezing levels of female Tyro3 KO mice, which are consistent with the female bias in anxiety disorders and could be indicative of maladaptive stress-responses. This study has demonstrated that white matter pathology related to a loss of Tyro3 is associated with pro-anxiety behavioural responses of female mice. Future studies could probe their contribution to increased risk for neuropsychiatric disorders when combined with stressful triggering events.

Functional and molecular characterization of suicidality factors using phenotypic and genome-wide data

Genome-wide association studies (GWAS) of suicidal thoughts and behaviors support the existence of genetic contributions. Continuous measures of psychiatric disorder symptom severity can sometimes model polygenic risk better than binarized definitions. We compared two severity measures of suicidal thoughts and behaviors at the molecular and functional levels using genome-wide data. We used summary association data from GWAS of four traits analyzed in 122,935 individuals of European ancestry: thought life was not worth living (TLNWL), thoughts of self-harm, actual self-harm, and attempted suicide. A new trait for suicidal thoughts and behaviors was constructed first, phenotypically, by aggregating the previous four traits (termed "suicidality") and second, genetically, by using genomic structural equation modeling (gSEM; termed S-factor). Suicidality and S-factor were compared using SNP-heritability (h2) estimates, genetic correlation (rg), partitioned h2, effect size distribution, transcriptomic correlations (ρGE) in the brain, and cross-population polygenic scoring (PGS). The S-factor had good model fit (χ2 = 0.21, AIC = 16.21, CFI = 1.00, SRMR = 0.024). Suicidality (h2 = 7.6%) had higher h2 than the S-factor (h2 = 2.54, Pdiff = 4.78 × 10-13). Although the S-factor had a larger number of non-null susceptibility loci (πc = 0.010), these loci had small effect sizes compared to those influencing suicidality (πc = 0.005, Pdiff = 0.045). The h2 of both traits was enriched for conserved biological pathways. The rg and ρGE support highly overlapping genetic and transcriptomic features between suicidality and the S-factor. PGS using European-ancestry SNP effect sizes strongly associated with TLNWL in Admixed Americans: Nagelkerke's R2 = 8.56%, P = 0.009 (PGSsuicidality) and Nagelkerke's R2 = 7.48%, P = 0.045 (PGSS-factor). An aggregate suicidality phenotype was statistically more heritable than the S-factor across all analyses and may be more informative for future genetic study designs interested in common genetic factors among different suicide related phenotypes.

DNA methylation and general psychopathology in childhood: an epigenome-wide meta-analysis from the PACE consortium

The general psychopathology factor (GPF) has been proposed as a way to capture variance shared between psychiatric symptoms. Despite a growing body of evidence showing both genetic and environmental influences on GPF, the biological mechanisms underlying these influences remain unclear. In the current study, we conducted epigenome-wide meta-analyses to identify both probe- and region-level associations of DNA methylation (DNAm) with school-age general psychopathology in six cohorts from the Pregnancy And Childhood Epigenetics (PACE) Consortium. DNAm was examined both at birth (cord blood; prospective analysis) and during school-age (peripheral whole blood; cross-sectional analysis) in total samples of N = 2178 and N = 2190, respectively. At school-age, we identified one probe (cg11945228) located in the Bromodomain-containing protein 2 gene (BRD2) that negatively associated with GPF (p = 8.58 × 10-8). We also identified a significant differentially methylated region (DMR) at school-age (p = 1.63 × 10-8), implicating the SHC Adaptor Protein 4 (SHC4) gene and the EP300-interacting inhibitor of differentiation 1 (EID1) gene that have been previously implicated in multiple types of psychiatric disorders in adulthood, including obsessive compulsive disorder, schizophrenia, and major depressive disorder. In contrast, no prospective associations were identified with DNAm at birth. Taken together, results of this study revealed some evidence of an association between DNAm at school-age and GPF. Future research with larger samples is needed to further assess DNAm variation associated with GPF.

Genome-wide association study identifies four pan-ancestry loci for suicidal ideation in the Million Veteran Program

Suicidal ideation (SI) often precedes and predicts suicide attempt and death, is the most common suicidal phenotype and is over-represented in veterans. The genetic architecture of SI in the absence of suicide attempt (SA) is unknown, yet believed to have distinct and overlapping risk with other suicidal behaviors. We performed the first GWAS of SI without SA in the Million Veteran Program (MVP), identifying 99,814 SI cases from electronic health records without a history of SA or suicide death (SD) and 512,567 controls without SI, SA or SD. GWAS was performed separately in the four largest ancestry groups, controlling for sex, age and genetic substructure. Ancestry-specific results were combined via meta-analysis to identify pan-ancestry loci. Four genome-wide significant (GWS) loci were identified in the pan-ancestry meta-analysis with loci on chromosomes 6 and 9 associated with suicide attempt in an independent sample. Pan-ancestry gene-based analysis identified GWS associations with DRD2, DCC, FBXL19, BCL7C, CTF1, ANNK1, and EXD3. Gene-set analysis implicated synaptic and startle response pathways (q's<0.05). European ancestry (EA) analysis identified GWS loci on chromosomes 6 and 9, as well as GWS gene associations in EXD3, DRD2, and DCC. No other ancestry-specific GWS results were identified, underscoring the need to increase representation of diverse individuals. The genetic correlation of SI and SA within MVP was high (rG = 0.87; p = 1.09e-50), as well as with post-traumatic stress disorder (PTSD; rG = 0.78; p = 1.98e-95) and major depressive disorder (MDD; rG = 0.78; p = 8.33e-83). Conditional analysis on PTSD and MDD attenuated most pan-ancestry and EA GWS signals for SI without SA to nominal significance, with the exception of EXD3 which remained GWS. Our novel findings support a polygenic and complex architecture for SI without SA which is largely shared with SA and overlaps with psychiatric conditions frequently comorbid with suicidal behaviors.

Cross-continental environmental and genome-wide association study on children and adolescent anxiety and depression

Anxiety and depression in children and adolescents warrant special attention as a public health issue given their devastating and long-term effects on development and mental health. Multiple factors, ranging from genetic vulnerabilities to environmental stressors, influence the risk for the disorders. This study aimed to understand how environmental factors and genomics affect children and adolescents anxiety and depression across three cohorts: Adolescent Brain and Cognitive Development Study (US, age of 9-10), Consortium on Vulnerability to Externalizing Disorders and Addictions (INDIA, age of 6-17) and IMAGEN (EUROPE, age of 14). We performed data harmonization and identified the environmental impact on anxiety/depression using a linear mixed-effect model, recursive feature elimination regression, and the LASSO regression model. Subsequently, genome-wide association analyses with consideration of significant environmental factors were performed for all three cohorts by mega-analysis and meta-analysis, followed by functional annotations. The results showed that multiple environmental factors contributed to the risk of anxiety and depression during development, where early life stress and school risk had the most significant and consistent impact across all three cohorts. Both meta and mega-analysis identified a novel SNP rs79878474 in chr11p15 to be the most promising SNP associated with anxiety and depression. Gene set analysis on the common genes mapped from top promising SNPs of both meta and mega analyses found significant enrichment in regions of chr11p15 and chr3q26, in the function of potassium channels and insulin secretion, in particular Kv3, Kir-6.2, SUR potassium channels encoded by the KCNC1, KCNJ11, and ABCCC8 genes respectively, in chr11p15. Tissue enrichment analysis showed significant enrichment in the small intestine and a trend of enrichment in the cerebellum. Our findings provide evidence of consistent environmental impact from early life stress and school risks on anxiety and depression during development and also highlight the genetic association between mutations in potassium channels along with the potential role of the cerebellum region, which are worthy of further investigation.

New insights from the last decade of research in psychiatric genetics: discoveries, challenges and clinical implications

Psychiatric genetics has made substantial progress in the last decade, providing new insights into the genetic etiology of psychiatric disorders, and paving the way for precision psychiatry, in which individual genetic profiles may be used to personalize risk assessment and inform clinical decision-making. Long recognized to be heritable, recent evidence shows that psychiatric disorders are influenced by thousands of genetic variants acting together. Most of these variants are commonly occurring, meaning that every individual has a genetic risk to each psychiatric disorder, from low to high. A series of large-scale genetic studies have discovered an increasing number of common and rare genetic variants robustly associated with major psychiatric disorders. The most convincing biological interpretation of the genetic findings implicates altered synaptic function in autism spectrum disorder and schizophrenia. However, the mechanistic understanding is still incomplete. In line with their extensive clinical and epidemiological overlap, psychiatric disorders appear to exist on genetic continua and share a large degree of genetic risk with one another. This provides further support to the notion that current psychiatric diagnoses do not represent distinct pathogenic entities, which may inform ongoing attempts to reconceptualize psychiatric nosology. Psychiatric disorders also share genetic influences with a range of behavioral and somatic traits and diseases, including brain structures, cognitive function, immunological phenotypes and cardiovascular disease, suggesting shared genetic etiology of potential clinical importance. Current polygenic risk score tools, which predict individual genetic susceptibility to illness, do not yet provide clinically actionable information. However, their precision is likely to improve in the coming years, and they may eventually become part of clinical practice, stressing the need to educate clinicians and patients about their potential use and misuse. This review discusses key recent insights from psychiatric genetics and their possible clinical applications, and suggests future directions.

Identification of Novel, Replicable Genetic Risk Loci for Suicidal Thoughts and Behaviors Among US Military Veterans

Importance

Suicide is a leading cause of death; however, the molecular genetic basis of suicidal thoughts and behaviors (SITB) remains unknown.

Objective

To identify novel, replicable genomic risk loci for SITB.

Design, Setting, and Participants

This genome-wide association study included 633 778 US military veterans with and without SITB, as identified through electronic health records. GWAS was performed separately by ancestry, controlling for sex, age, and genetic substructure. Cross-ancestry risk loci were identified through meta-analysis. Study enrollment began in 2011 and is ongoing. Data were analyzed from November 2021 to August 2022.

Main Outcome and Measures

SITB.

Results

A total of 633 778 US military veterans were included in the analysis (57 152 [9%] female; 121 118 [19.1%] African ancestry, 8285 [1.3%] Asian ancestry, 452 767 [71.4%] European ancestry, and 51 608 [8.1%] Hispanic ancestry), including 121 211 individuals with SITB (19.1%). Meta-analysis identified more than 200 GWS (P < 5 × 10-8) cross-ancestry risk single-nucleotide variants for SITB concentrated in 7 regions on chromosomes 2, 6, 9, 11, 14, 16, and 18. Top single-nucleotide variants were largely intronic in nature; 5 were independently replicated in ISGC, including rs6557168 in ESR1, rs12808482 in DRD2, rs77641763 in EXD3, rs10671545 in DCC, and rs36006172 in TRAF3. Associations for FBXL19 and AC018880.2 were not replicated. Gene-based analyses implicated 24 additional GWS cross-ancestry risk genes, including FURIN, TSNARE1, and the NCAM1-TTC12-ANKK1-DRD2 gene cluster. Cross-ancestry enrichment analyses revealed significant enrichment for expression in brain and pituitary tissue, synapse and ubiquitination processes, amphetamine addiction, parathyroid hormone synthesis, axon guidance, and dopaminergic pathways. Seven other unique European ancestry-specific GWS loci were identified, 2 of which (POM121L2 and METTL15/LINC02758) were replicated. Two additional GWS ancestry-specific loci were identified within the African ancestry (PET112/GATB) and Hispanic ancestry (intergenic locus on chromosome 4) subsets, both of which were replicated. No GWS loci were identified within the Asian ancestry subset; however, significant enrichment was observed for axon guidance, cyclic adenosine monophosphate signaling, focal adhesion, glutamatergic synapse, and oxytocin signaling pathways across all ancestries. Within the European ancestry subset, genetic correlations (r > 0.75) were observed between the SITB phenotype and a suicide attempt-only phenotype, depression, and posttraumatic stress disorder. Additionally, polygenic risk score analyses revealed that the Million Veteran Program polygenic risk score had nominally significant main effects in 2 independent samples of veterans of European and African ancestry.

Conclusions and Relevance

The findings of this analysis may advance understanding of the molecular genetic basis of SITB and provide evidence for ESR1, DRD2, TRAF3, and DCC as cross-ancestry candidate risk genes. More work is needed to replicate these findings and to determine if and how these genes might impact clinical care.

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.

Higher polygenic risk scores for anxiety disorders are associated with reduced area in the anterior cingulate gyrus

Anxiety disorders are heterogeneous, show a moderate genetic contribution and are associated with inconsistent cortical structure alterations. Here, we investigated whether genetic factors for anxiety disorders contribute to cortical alterations by conducting polygenic risk score (PRS) analyses. We calculated PRSs for anxiety disorders at several P value thresholds (from P_T ≤ 5.0 × 10^-8 to P_T ≤ 1.0) based on the latest large-scale genome-wide association study of anxiety disorders from the UK biobank (25,453 cases; 58,113 controls) in an independent sample of psychiatrically and physically healthy subjects (n = 174). Using regression after adjusting for confounding factors, we tested whether these PRSs were associated with the surface area and cortical thickness in 34 bilateral brain regions extracted using FreeSurfer. A higher PRS for anxiety disorders at P_T ≤ 1.0 was significantly associated with a reduced right caudal anterior cingulate area (beta = -0.25, p_uncorrected = 9.51 × 10^-4, p_corrected = 0.032). PRSs based on more common SNPs, especially from P_T ≤ 0.01 to P_T ≤ 1.0, were associated with the right caudal anterior cingulate area (a maximum at P_T ≤ 0.5: R² = 0.066, beta = -0.27, p_uncorr = 3.81 × 10^-4, p_corr = 0.013). Furthermore, individuals in the highest quartile for anxiety disorder PRS had lower surface area and volume in the right anterior cingulate gyrus than those in the lowest quartile. We suggest a shared genetic etiology between anxiety disorders and structural features of the anterior cingulate gyrus, possibly contributing to the pathogenesis of anxiety disorders via emotional dysregulations. Our findings suggest the potential usefulness of PRS to reduce pathological heterogeneity among anxiety disorders.

Adolescent anxiety and pain problems: A joint, genome-wide investigation and pathway-based analysis

Both common pain and anxiety problems are widespread, debilitating and often begin in childhood-adolescence. Twin studies indicate that this co-occurrence is likely due to shared elements of risk, rather than reciprocal causation. A joint genome-wide investigation and pathway/network-based analysis of adolescent anxiety and pain problems can identify genetic pathways that subserve shared etiopathogenetic mechanisms. Pathway-based analyses were performed in the independent samples of: The Quebec Newborn Twin Study (QNTS; 246 twin pairs and 321 parents), the Longitudinal Study of Child Development in Quebec (QLSCD; n = 754), and in the combined QNTS and QLSCD sample. Multiple suggestive associations (p<1×10-5), and several enriched pathways were found after FDR correction for both phenotypes in the QNTS; many nominally-significant enriched pathways overlapped between pain problems and anxiety symptoms (uncorrected p<0.05) and yielded results consistent with previous studies of pain or anxiety. The QLSCD and the combined QNTS and QLSCD sample yielded similar findings. We replicated an association between the pathway involved in the regulation of myotube differentiation (GO:0010830) and both pain and anxiety problems in the QLSDC and the combined QNTS and QLSCD sample. Although limited by sample size and thus power, these data provide an initial support to conjoint molecular investigations of adolescent pain and anxiety problems. Understanding the etiology underlying pain and anxiety co-occurrence in this age range is relevant to address the nature of comorbidity and its developmental pathways, and shape intervention. The replication across samples implies that these effects are reliable and possess external validity.

Childhood inhibition predicts adolescent social anxiety: Findings from a longitudinal twin study

An enduring issue in the study of mental health is identifying developmental processes that explain how childhood characteristics progress to maladaptive forms. We examine the role that behavioral inhibition (BI) has on social anxiety (SA) during adolescence in 868 families of twins assessed at ages 8, 13, and 15 years. Multimodal assessments of BI and SA were completed at each phase, with additional measures (e.g., parenting stress) for parents and twins. Analyses were conducted in several steps: first, we used a cross-lagged panel model to demonstrate bidirectional paths between BI and SA; second a biometric Cholesky decomposition showed that both genetic and environmental influences on childhood BI also affect adolescent SA; next, multilevel phenotypic models tested moderation effects between BI and SA. We tested seven potential moderators of the BI to SA prediction in individual models and included only those that emerged as significant in a final conditional model examining predictors of SA. Though several main effects emerged as significant, only parenting stress had a significant interaction with BI to predict SA, highlighting the importance of environmental moderators in models examining temperamental effects on later psychological symptoms. This comprehensive assessment continues to build the prototype for such developmental psychopathology models.

Exploring the genetic overlap between twelve psychiatric disorders

The widespread comorbidity among psychiatric disorders demonstrated in epidemiological studies1-5 is mirrored by non-zero, positive genetic correlations from large-scale genetic studies6-10. To identify shared biological processes underpinning this observed phenotypic and genetic covariance and enhance molecular characterization of general psychiatric disorder liability11-13, we used several strategies aimed at uncovering pleiotropic, that is, cross-trait-associated, single-nucleotide polymorphisms (SNPs), genes and biological pathways. We conducted cross-trait meta-analysis on 12 psychiatric disorders to identify pleiotropic SNPs. The meta-analytic signal was driven by schizophrenia, hampering interpretation and joint biological characterization of the cross-trait meta-analytic signal. Subsequent pairwise comparisons of psychiatric disorders identified substantial pleiotropic overlap, but mainly among pairs of psychiatric disorders, and mainly at less stringent P-value thresholds. Only annotations related to evolutionarily conserved genomic regions were significant for multiple (9 out of 12) psychiatric disorders. Overall, identification of shared biological mechanisms remains challenging due to variation in power and genetic architecture between psychiatric disorders.

Collective genomic segments with differential pleiotropic patterns between cognitive dimensions and psychopathology

Cognitive deficits are known to be related to most forms of psychopathology. Here, we perform local genetic correlation analysis as a means of identifying independent segments of the genome that show biologically interpretable pleiotropic associations between cognitive dimensions and psychopathology. We identify collective segments of the genome, which we call "meta-loci", showing differential pleiotropic patterns for psychopathology relative to either cognitive task performance (CTP) or performance on a non-cognitive factor (NCF) derived from educational attainment. We observe that neurodevelopmental gene sets expressed during the prenatal-early childhood period predominate in CTP-relevant meta-loci, while post-natal gene sets are more involved in NCF-relevant meta-loci. Further, we demonstrate that neurodevelopmental gene sets are dissociable across CTP meta-loci with respect to their spatial distribution across the brain. Additionally, we find that GABA-ergic, cholinergic, and glutamatergic genes drive pleiotropic relationships within dissociable meta-loci.

Quality Control Procedures for Genome-Wide Association Studies

Genome-wide association studies (GWAS) are being conducted at an unprecedented rate in population-based cohorts and have increased our understanding of the pathophysiology of many complex diseases. Regardless of the context, the practical utility of this information ultimately depends upon the quality of the data used for statistical analyses. Quality control (QC) procedures for GWAS are constantly evolving. Here, we enumerate some of the challenges in QC of genotyped GWAS data and describe the approaches involving genotype imputation of a sample dataset along with post-imputation quality assurance, thereby minimizing potential bias and error in GWAS results. We discuss common issues associated with QC of the GWAS data (genotyped and imputed), including data file formats, software packages for data manipulation and analysis, sex chromosome anomalies, sample identity, sample relatedness, population substructure, batch effects, and marker quality. We provide detailed guidelines along with a sample dataset to suggest current best practices and discuss areas of ongoing and future research. © 2022 Wiley Periodicals LLC.

Genome-wide Scan of Dental Fear and Anxiety Nominates Novel Genes

Dental care-related fear and anxiety (DFA) is prevalent, affects oral health care utilization, and is related to poor oral health and decreased quality of life. In addition to learned and cultural factors, genetics is hypothesized to contribute to DFA. Therefore, we performed a genome-wide association study to identify genetic variants contributing to DFA. Adult and adolescent participants were from 4 cohorts (3 from the US-based Center for Oral Health Research in Appalachia, n = 1,144, 1,164, and 535, and the UK-based Avon Longitudinal Study of Parents and Children [ALSPAC], n = 2,078). Two self-report instruments were used to assess DFA: the Dental Fear Survey (US cohorts) and Corah's Dental Anxiety Scale (ALSPAC). Genome-wide scans were performed for the DFA total scores and subscale scores (avoidance, physiological arousal, fear of dental treatment-specific stimuli), adjusting for age, sex, educational attainment, recruitment site, and genetic ancestry. Results across cohorts were combined using meta-analysis. Heritability estimates for DFA total and subscale scores were similar across cohorts and ranged from 23% to 59%. The meta-analysis revealed 3 significant (P < 5E-8) associations between genetic loci and 2 DFA subscales: physiological arousal and avoidance. Nearby genes included NTSR1 (P = 3.05E-8), DMRTA1 (P = 4.40E-8), and FAM84A (P = 7.72E-9). Of these, NTSR1, which was associated with the avoidance subscale, mediates neurotensin function, and its deficiency may lead to altered fear memory in mice. Gene enrichment analyses indicated that loci associated with the DFA total score and physiological arousal subscale score were enriched for genes associated with severe and persistent mental health (e.g., schizophrenia) and neurocognitive (e.g., autism) disorders. Heritability analysis indicated that DFA is partly explained by genetic factors, and our association results suggested shared genetic underpinnings with other psychological conditions.

Causality of genetically determined metabolites on anxiety disorders: a two-sample Mendelian randomization study

BACKGROUND

Although anxiety disorders are one of the most prevalent mental disorders, their underlying biological mechanisms have not yet been fully elucidated. In recent years, genetically determined metabolites (GDMs) have been used to reveal the biological mechanisms of mental disorders. However, this strategy has not been applied to anxiety disorders. Herein, we explored the causality of GDMs on anxiety disorders through Mendelian randomization study, with the overarching goal of unraveling the biological mechanisms.

METHODS

A two-sample Mendelian randomization (MR) analysis was implemented to assess the causality of GDMs on anxiety disorders. A genome-wide association study (GWAS) of 486 metabolites was used as the exposure, whereas four different GWAS datasets of anxiety disorders were the outcomes. Notably, all datasets were acquired from publicly available databases. A genetic instrumental variable (IV) was used to explore the causality between the metabolite and anxiety disorders for each metabolite. The MR Steiger filtering method was implemented to examine the causality between metabolites and anxiety disorders. The standard inverse variance weighted (IVW) method was first used for the causality analysis, followed by three additional MR methods (the MR-Egger, weighted median, and MR-PRESSO (pleiotropy residual sum and outlier) methods) for sensitivity analyses in MR analysis. MR-Egger intercept, and Cochran's Q statistical analysis were used to evaluate possible heterogeneity and pleiotropy. Bonferroni correction was used to determine the causative association features (P < 1.03 × 10^-4). Furthermore, metabolic pathways analysis was performed using the web-based MetaboAnalyst 5.0 software. All statistical analysis were performed in R software. The STROBE-MR checklist for the reporting of MR studies was used in this study.

RESULTS

In MR analysis, 85 significant causative relationship GDMs were identified. Among them, 11 metabolites were overlapped in the four different datasets of anxiety disorders. Bonferroni correction showing1-linoleoylglycerophosphoethanolamine (OR_{fixed-effect IVW} = 1.04; 95% CI 1.021-1.06; P_{fixed-effect IVW} = 4.3 × 10^-5) was the most reliable causal metabolite. Our results were robust even without a single SNP because of a "leave-one-out" analysis. The MR-Egger intercept test indicated that genetic pleiotropy had no effect on the results (intercept = - 0.0013, SE = 0.0006, P = 0.06). No heterogeneity was detected by Cochran's Q test (MR-Egger. Q = 7.68, P = 0.742; IVW. Q = 12.12, P = 0.436). A directionality test conducted by MR Steiger confirmed our estimation of potential causal direction (P < 0.001). In addition, two significant pathways, the "primary bile acid biosynthesis" pathway (P = 0.008) and the "valine, leucine, and isoleucine biosynthesis" pathway (P = 0.03), were identified through metabolic pathway analysis.

CONCLUSION

This study provides new insights into the causal effects of GDMs on anxiety disorders by integrating genomics and metabolomics. The metabolites that drive anxiety disorders may be suited to serve as biomarkers and also will help to unravel the biological mechanisms of anxiety disorders.

Genome-Wide Investigation of Maximum Habitual Alcohol Intake in US Veterans in Relation to Alcohol Consumption Traits and Alcohol Use Disorder

Importance

Alcohol genome-wide association studies (GWASs) have generally focused on alcohol consumption and alcohol use disorder (AUD); few have examined habitual drinking behaviors like maximum habitual alcohol intake (MaxAlc).

Objectives

To identify genetic loci associated with MaxAlc and to elucidate the genetic architecture across alcohol traits.

Design, Setting, and Participants

This MaxAlc genetic association study was performed among Million Veteran Program participants enrolled from January 10, 2011, to September 30, 2020. Ancestry-specific GWASs were conducted in participants with European (n = 218 623) and African (n = 29 132) ancestry, then meta-analyzed (N = 247 755). Linkage-disequilibrium score regression was used to estimate single nucleotide variant (SNV)-heritability and genetic correlations (rg) with other alcohol and psychiatric traits. Genomic structural equation modeling (gSEM) was used to evaluate genetic associations between MaxAlc and other alcohol traits. Mendelian randomization was used to examine potential causal relationships between MaxAlc and liver enzyme levels. MTAG (multitrait analysis of GWAS) was used to analyze MaxAlc and problematic alcohol use (PAU) jointly.

Exposures

Genetic associations.

Main Outcomes and Measures

MaxAlc was defined from the following survey item: "in a typical month, what is/was the largest number of drinks of alcohol you may have had in one day?" with ordinal responses from 0 to 15 or more drinks.

Results

GWASs were conducted on sample sizes of as many as 247 455 US veterans. Participants were 92.68% male and had mean (SD) age of 65.92 (11.70) years. The MaxAlc GWAS resulted in 15 genome-wide significant loci. Top associations in European-ancestry and African-ancestry participants were with known functional variants in the ADH1B gene, namely rs1229984 (P = 3.12 × 10-101) and rs2066702 (P = 6.30 × 10-17), respectively. Novel associations were also found. SNV-heritability was 6.65% (SE, 0.41) in European-ancestry participants and 3.42% (SE, 1.46) in African-ancestry participants. MaxAlc was positively correlated with PAU (rg = 0.79; P = 3.95 × 10-149) and AUD (rg = 0.76; P = 1.26 × 10-127) and had negative rg with the UK Biobank "alcohol usually taken with meals" (rg = -0.53; P = 1.40 × 10-50). For psychiatric traits, MaxAlc had the strongest genetic correlation with suicide attempt (rg = 0.40; P = 3.02 × 10-21). gSEM supported a 2-factor model with MaxAlc loading on a factor with PAU and AUD and other alcohol consumption measures loading on a separate factor. Mendelian randomization supported an association between MaxAlc and the liver enzyme gamma-glutamyltransferase (β = 0.012; P = 2.66 × 10-10). MaxAlc MTAG resulted in 31 genome-wide significant loci.

Conclusions and Relevance

The findings suggest that MaxAlc closely aligns genetically with PAU traits. This study improves understanding of the mechanisms associated with normative alcohol consumption vs problematic habitual use and AUD as well as how MaxAlc relates to psychiatric and medical conditions genetically and biologically.