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

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

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

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

Differential Serum Levels of CACNA1C, Circadian Rhythm and Stress Response Molecules in Subjects with Bipolar Disorder: Associations with Genetic and Clinical Factors

Background

Many patients with bipolar disorder (BD) do not respond to or have difficulties tolerating lithium and/or other mood stabilizing agents. There is a need for personalized treatments based on biomarkers in guiding treatment options. The calcium voltage-gated channel CACNA1C is a promising candidate for developing personalized treatments. CACNA1C is implicated in BD by genome-wide association studies and several lines of evidence suggest that targeting L-type calcium channels could be an effective treatment strategy. However, before such individualized treatments can be pursued, biomarkers predicting treatment response need to be developed.

Methods

As a first step in testing the hypothesis that CACNA1C genotype is associated with serum levels of CACNA1C, we conducted ELISA measures on serum samples from 100 subjects with BD and 100 control subjects.

Results

We observed significantly higher CACNA1C (p<0.01) protein levels in subjects with BD. The risk SNP (rs11062170) showed functional significance as subjects homozygous for the risk allele (CC) had significantly greater CACNA1C protein levels compared to subjects with one (p=0.013) or no copies (p=0.009). We observed higher somatostatin (SST) (p<0.003) protein levels and lower levels of the clock protein ARTNL (p<0.03) and stress signaling factor corticotrophin releasing hormone (CRH) (p<0.001) in BD. SST and PER2 protein levels were associated with both alcohol dependence and lithium response.

Conclusions

Our findings represent the first evidence for increased serum levels of CACNA1C in BD. Along with altered levels of SST, ARNTL, and CRH our findings suggest CACNA1C is associated with circadian rhythm and stress response disturbances in BD.

Molecular mechanisms involved in alcohol craving, IRF3, and endoplasmic reticulum stress: a multi-omics study

Alcohol use disorder (AUD) is the most prevalent substance use disorder worldwide. Acamprosate and naltrexone are anti-craving drugs used in AUD pharmacotherapy. However, molecular mechanisms underlying their anti-craving effect remain unclear. This study utilized a patient-derived induced pluripotent stem cell (iPSC)-based model system and anti-craving drugs that are used to treat AUD as "molecular probes" to identify possible mechanisms associated with alcohol craving. We examined the pathophysiology of craving and anti-craving drugs by performing functional genomics studies using iPSC-derived astrocytes and next-generation sequencing. Specifically, RNA sequencing performed using peripheral blood mononuclear cells from AUD patients with extreme values for alcohol craving intensity prior to treatment showed that inflammation-related pathways were highly associated with alcohol cravings. We then performed a genome-wide assessment of chromatin accessibility and gene expression profiles of induced iPSC-derived astrocytes in response to ethanol or anti-craving drugs. Those experiments identified drug-dependent epigenomic signatures, with IRF3 as the most significantly enriched motif in chromatin accessible regions. Furthermore, the activation of IRF3 was associated with ethanol-induced endoplasmic reticulum (ER) stress which could be attenuated by anti-craving drugs, suggesting that ER stress attenuation might be a target for anti-craving agents. In conclusion, we found that craving intensity was associated with alcohol consumption and treatment outcomes. Our functional genomic studies suggest possible relationships among craving, ER stress, IRF3 and the actions of anti-craving drugs.

Polygenic Prediction of Cellular Immune Responses to Mumps Vaccine

In this report, we provide a follow-up analysis of a previously published genome-wide association study of host genetic variants associated with inter-individual variations in cellular immune responses to mumps vaccine. Here we report the results of a polygenic score (PGS) analysis showing how common variants can predict mumps vaccine response. We found higher PGS for IFNγ, IL-2, and TNFα were predictive of higher post-vaccine IFNγ (p-value = 2e-6), IL-2 (p = 2e-7), and TNFα (p = 0.004) levels, respectively. Control of immune responses after vaccination is complex and polygenic in nature. Our results suggest that the PGS-based approach enables better capture of the combined genetic effects that contribute to mumps vaccine-induced immunity, potentially offering a more comprehensive understanding than traditional single-variant GWAS. This approach will likely have broad utility in studying genetic control of immune responses to other vaccines and to infectious diseases.

Pharmacogenomic overlap between antidepressant treatment response in major depression & antidepressant associated treatment emergent mania in bipolar disorder

There is increasing interest in individualizing treatment selection for more than 25 regulatory approved treatments for major depressive disorder (MDD). Despite an inconclusive efficacy evidence base, antidepressants (ADs) are prescribed for the depressive phase of bipolar disorder (BD) with oftentimes, an inadequate treatment response and or clinical concern for mood destabilization. This study explored the relationship between antidepressant response in MDD and antidepressant-associated treatment emergent mania (TEM) in BD. We conducted a genome-wide association study (GWAS) and polygenic score analysis of TEM and tested its association in a subset of BD-type I patients treated with SSRIs or SNRIs. Our results did not identify any genome-wide significant variants although, we found that a higher polygenic score (PGS) for antidepressant response in MDD was associated with higher odds of TEM in BD. Future studies with larger transdiagnostic depressed cohorts treated with antidepressants are encouraged to identify a neurobiological mechanism associated with a spectrum of depression improvement from response to emergent mania.

Patterns of pharmacotherapy for bipolar disorder: A GBC survey

OBJECTIVES

To understand treatment practices for bipolar disorders (BD), this study leveraged the Global Bipolar Cohort collaborative network to investigate pharmacotherapeutic treatment patterns in multiple cohorts of well-characterized individuals with BD in North America, Europe, and Australia.

METHODS

Data on pharmacotherapy, demographics, diagnostic subtypes, and comorbidities were provided from each participating cohort. Individual site and regional pooled proportional meta-analyses with generalized linear mixed methods were conducted to identify prescription patterns.

RESULTS

This study included 10,351 individuals from North America (n = 3985), Europe (n = 3822), and Australia (n = 2544). Overall, participants were predominantly female (60%) with BD-I (60%; vs. BD-II = 33%). Cross-sectionally, mood-stabilizing anticonvulsants (44%), second-generation antipsychotics (42%), and antidepressants (38%) were the most prescribed medications. Lithium was prescribed in 29% of patients, primarily in the Australian (31%) and European (36%) cohorts. First-generation antipsychotics were prescribed in 24% of the European versus 1% in the North American cohort. Antidepressant prescription rates were higher in BD-II (47%) compared to BD-I (35%). Major limitations were significant differences among cohorts based on inclusion/exclusion criteria, data source, and time/year of enrollment into cohort.

CONCLUSIONS

Mood-stabilizing anticonvulsants, second-generation antipsychotics, and antidepressants were the most prescribed medications suggesting prescription patterns that are not necessarily guideline concordant. Significant differences exist in the prescription practices across different geographic regions, especially the underutilization of lithium in the North American cohorts and the higher utilization of first-generation antipsychotics in the European cohorts. There is a need to conduct future longitudinal studies to further explore these differences and their impact on outcomes, and to inform and implement evidence-based guidelines to help improve treatment practices in BD.

Prevalence and associations of multiple hypnotic prescriptions in a clinical sample

STUDY OBJECTIVES

We examined the prevalence of multiple hypnotic prescriptions and its association with clinical and demographic characteristics from the electronic health record (EHR) in the Mayo Clinic Biobank (MCB).

METHODS

Adult subjects enrolled in the MCB with an EHR length of ≥1 year were included (n=52,940). Clinical and demographic characteristics were compared between subjects who were and were not prescribed any hypnotic approved for insomnia by the United States Food and Drug Administration and/or trazodone and in those prescribed a single versus multiple (≥2) hypnotics. A phenotype-based phenome-wide association study (PheWAS) examining associations between hypnotic prescriptions and diagnoses across the EHR was performed adjusting for demographic and other confounders.

RESULTS

17,662 (33%) subjects were prescribed at least one hypnotic and 5,331 (10%) received (≥2) hypnotics. Subjects who were prescribed a hypnotic were more likely to be older, female, White, with a longer EHR, and a greater number of diagnostic codes (all p<0.001). Those with multiple hypnotic prescriptions were more likely to be younger, female, with a longer EHR, and greater number of diagnostic codes (all p<0.001) compared to those prescribed a single hypnotic. The PheWAS revealed that subjects with multiple hypnotic prescriptions had higher rates of mood disorders, anxiety disorders, suicidal ideation, restless legs syndrome, and chronic pain (all p<1e-10).

CONCLUSIONS

Receiving multiple hypnotic prescriptions is common and associated with greater prevalence of psychiatric, chronic pain, and sleep-related movement disorders. Future studies should examine potential genetic associations with multiple hypnotic prescriptions to personalize treatments for chronic insomnia.

A genome-wide association study of antidepressant-induced mania

Antidepressant-induced mania (AIM) is a side effect of antidepressant treatment that is characterized by mania or hypomania after the start of medication. It is likely polygenic, but its genetic component remains largely unexplored. We aim to conduct the first genome-wide association study of AIM in 814 bipolar disorder patients of European ancestry. We report no significant findings from our single-marker or gene-based analyses. Our polygenic risk score analyses also did not yield significant results with bipolar disorder, antidepressant response, or lithium response. Our suggestive findings on the hypothalamic-pituitary-adrenal axis and the opioid system in AIM require independent replications.

Factors associated with suicide attempts in the antecedent illness trajectory of bipolar disorder and schizophrenia

BACKGROUND

Factors associated with suicide attempts during the antecedent illness trajectory of bipolar disorder (BD) and schizophrenia (SZ) are poorly understood.

METHODS

Utilizing the Rochester Epidemiology Project, individuals born after 1985 in Olmsted County, MN, presented with first episode mania (FEM) or psychosis (FEP), subsequently diagnosed with BD or SZ were identified. Patient demographics, suicidal ideation with plan, self-harm, suicide attempts, psychiatric hospitalizations, substance use, and childhood adversities were quantified using the electronic health record. Analyses pooled BD and SZ groups with a transdiagnostic approach given the two diseases were not yet differentiated. Factors associated with suicide attempts were examined using bivariate methods and multivariable logistic regression modeling.

RESULTS

A total of 205 individuals with FEM or FEP (BD = 74, SZ = 131) were included. Suicide attempts were identified in 39 (19%) patients. Those with suicide attempts during antecedent illness trajectory were more likely to be female, victims of domestic violence or bullying behavior, and have higher rates of psychiatric hospitalizations, suicidal ideation with plan and/or self-harm, as well as alcohol, drug, and nicotine use before FEM/FEP onset. Based on multivariable logistic regression, three factors remained independently associated with suicidal attempts: psychiatric hospitalization (OR = 5.84, 95% CI 2.09-16.33, p < 0.001), self-harm (OR = 3.46, 95% CI 1.29-9.30, p = 0.014), and nicotine use (OR = 3.02, 95% CI 1.17-7.76, p = 0.022).

CONCLUSION

Suicidal attempts were prevalent during the antecedents of BD and SZ and were associated with several risk factors before FEM/FEP. Their clinical recognition could contribute to improve early prediction and prevention of suicide during the antecedent illness trajectory of BD and SZ.

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.

Polygenic prediction of bipolar disorder in a Latin American sample

To date, bipolar disorder (BD) genetic studies and polygenic risk scores (PRSs) for BD are based primarily on populations of European descent (EUR) and lack representation from other ancestries including Latin American (LAT). Here, we describe a new LAT cohort from the Mayo Clinic Bipolar Biobank (MCBB), a multisite collaboration with recruitment sites in the United States (EUR; 1,443 cases and 777 controls) and Mexico and Chile (LAT; 211 cases and 161 controls) and use the sample to explore the performance of a BD-PRS in a LAT population. Using results from the largest genome-wide association study of BD in EUR individuals, PRSice2 and LDpred2 were used to compute BD-PRSs in the LAT and EUR samples from the MCBB. PRSs explained up to 1.4% (PRSice) and 4% (LDpred2) of the phenotypic variance on the liability scale in the LAT sample compared to 3.8% (PRSice2) and 3.4% (LDpred2) in the EUR samples. Future larger studies should further explore the differential performance of different PRS approaches across ancestries. International multisite studies, such as this one, have the potential to address diversity-related limitations of prior genomic studies and ultimately contribute to the reduction of health disparities.

Multi-ancestry meta-analysis of tobacco use disorder prioritizes novel candidate risk genes and reveals associations with numerous health outcomes

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

Racial differences in pathways to care preceding first episode mania or psychosis: a historical cohort prodromal study

Background

There is evidence suggesting racial disparities in diagnosis and treatment in bipolar disorder (BD) and schizophrenia (SZ). The purpose of this study is to compare psychiatric diagnoses and psychotropic use preceding a first episode of mania (FEM) or psychosis (FEP) in racially diverse patients.

Methods

Using a comprehensive medical records linkage system (Rochester Epidemiology Project, REP), we retrospectively identified individuals diagnosed with BD or SZ and a documented first episode of mania or psychosis. Illness trajectory before FEP/FEM were characterized as the time from first visit for a mental health complaint to incident case. Pathways to care and clinical events preceding FEP/FEM were compared based on subsequent incident case diagnosis (BD or SZ) and self-reported race (White vs. non-White).

Results

A total of 205 (FEM = 74; FEP = 131) incident cases were identified in the REP. Duration of psychiatric antecedents was significantly shorter in non-White patients, compared to White patients (2.2 ± 4.3 vs. 7.4 ± 6.6 years; p < 0.001) with an older age at time of first visit for a mental health complaint (15.7 ± 6.3 vs. 11.1 ± 6.0 years; p = 0.005). There were no significant differences by race in FEM pathway to care or age of first seeking mental health. Overall non-White patients had lower rates of psychotropic use.

Conclusion

These data are unable to ascertain reasons for shorter duration of psychiatric antecedents and later age of seeking care, and more broadly first age of initial symptom presentation. If symptoms are confirmed to be earlier than first time seeking care in both groups, it would be important to identify barriers that racial minorities face to access timely psychiatric care and optimize early intervention strategies.

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

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

Antidepressant-Associated Treatment Emergent Mania: A Meta-Analysis to Guide Risk Modeling Pharmacogenomic Targets of Potential Clinical Value

BACKGROUND

The purpose of this study was to review the association between the SLC6A4 5-HTTLPR polymorphism and antidepressant (AD)-associated treatment emergent mania (TEM) in bipolar disorder alongside starting a discussion on the merits of developing risk stratification models to guide when not to provide AD treatment for bipolar depression.

METHODS

Studies that examined the association between clinical and genetic risk factors, specifically monoaminergic transporter genetic variation, and TEM were identified. A meta-analysis was performed using the odds ratio to estimate the effect size under the Der-Simonian and Laird model.

RESULTS

Seven studies, referencing the SLC6A4 5-HTTLPR polymorphism and TEM (total N = 1578; TEM+ =594, TEM- = 984), of 142 identified articles were included. The time duration between the start of the AD to emergence of TEM ranged from 4 to 12 weeks. There was a nominally significant association between the s allele of the 5-HTTLPR polymorphism and TEM (odds ratio, 1.434; 95% confidence interval, 1.001-2.055; P = 0.0493; I2 = 52%). No studies have investigated norepinephrine or dopamine transporters.

CONCLUSION

Although the serotonin transporter genetic variation is commercially available in pharmacogenomic decision support tools, greater efforts, more broadly, should focus on complete genome-wide approaches to determine genetic variants that may contribute to TEM. Moreover, these data are exemplary to the merits of developing risk stratification models, which include both clinical and biological risk factors, to guide when not to use ADs in bipolar disorder. Future studies will need to validate new risk models that best inform the development of personalized medicine best practices treating bipolar depression.

Whole Genome Analysis of Venous Thromboembolism: the Trans-Omics for Precision Medicine Program

BACKGROUND

Risk for venous thromboembolism has a strong genetic component. Whole genome sequencing from the TOPMed program (Trans-Omics for Precision Medicine) allowed us to look for new associations, particularly rare variants missed by standard genome-wide association studies.

METHODS

The 3793 cases and 7834 controls (11.6% of cases were individuals of African, Hispanic/Latino, or Asian ancestry) were analyzed using a single variant approach and an aggregate gene-based approach using our primary filter (included only loss-of-function and missense variants predicted to be deleterious) and our secondary filter (included all missense variants).

RESULTS

Single variant analyses identified associations at 5 known loci. Aggregate gene-based analyses identified only PROC (odds ratio, 6.2 for carriers of rare variants; P=7.4×10-14) when using our primary filter. Employing our secondary variant filter led to a smaller effect size at PROC (odds ratio, 3.8; P=1.6×10-14), while excluding variants found only in rare isoforms led to a larger one (odds ratio, 7.5). Different filtering strategies improved the signal for 2 other known genes: PROS1 became significant (minimum P=1.8×10-6 with the secondary filter), while SERPINC1 did not (minimum P=4.4×10-5 with minor allele frequency <0.0005). Results were largely the same when restricting the analyses to include only unprovoked cases; however, one novel gene, MS4A1, became significant (P=4.4×10-7 using all missense variants with minor allele frequency <0.0005).

CONCLUSIONS

Here, we have demonstrated the importance of using multiple variant filtering strategies, as we detected additional genes when filtering variants based on their predicted deleteriousness, frequency, and presence on the most expressed isoforms. Our primary analyses did not identify new candidate loci; thus larger follow-up studies are needed to replicate the novel MS4A1 locus and to identify additional rare variation associated with venous thromboembolism.

The importance of social activity to risk of major depression in older adults

BACKGROUND

Several social determinants of health (SDoH) have been associated with the onset of major depressive disorder (MDD). However, prior studies largely focused on individual SDoH and thus less is known about the relative importance (RI) of SDoH variables, especially in older adults. Given that risk factors for MDD may differ across the lifespan, we aimed to identify the SDoH that was most strongly related to newly diagnosed MDD in a cohort of older adults.

METHODS

We used self-reported health-related survey data from 41 174 older adults (50-89 years, median age = 67 years) who participated in the Mayo Clinic Biobank, and linked ICD codes for MDD in the participants' electronic health records. Participants with a history of clinically documented or self-reported MDD prior to survey completion were excluded from analysis (N = 10 938, 27%). We used Cox proportional hazards models with a gradient boosting machine approach to quantify the RI of 30 pre-selected SDoH variables on the risk of future MDD diagnosis.

RESULTS

Following biobank enrollment, 2073 older participants were diagnosed with MDD during the follow-up period (median duration = 6.7 years). The most influential SDoH was perceived level of social activity (RI = 0.17). Lower level of social activity was associated with a higher risk of MDD [hazard ratio = 2.27 (95% CI 2.00-2.50) for highest v. lowest level].

CONCLUSION

Across a range of SDoH variables, perceived level of social activity is most strongly related to MDD in older adults. Monitoring changes in the level of social activity may help identify older adults at an increased risk of MDD.

Direct Measurement of Hexacontatetrapole, E6 γ Decay from ^{53m}Fe

The only proposed observation of a discrete, hexacontatetrapole (E6) transition in nature occurs from the T_{1/2}=2.54(2)-min decay of ^{53m}Fe. However, there are conflicting claims concerning its γ-decay branching ratio, and a rigorous interrogation of γ-ray sum contributions is lacking. Experiments performed at the Australian Heavy Ion Accelerator Facility were used to study the decay of ^{53m}Fe. For the first time, sum-coincidence contributions to the weak E6 and M5 decay branches have been firmly quantified using complementary experimental and computational methods. Agreement across the different approaches confirms the existence of the real E6 transition; the M5 branching ratio and transition rate have also been revised. Shell model calculations performed in the full fp model space suggest that the effective proton charge for high-multipole, E4 and E6, transitions is quenched to approximately two-thirds of the collective E2 value. Correlations between nucleons may offer an explanation of this unexpected phenomenon, which is in stark contrast to the collective nature of lower-multipole, electric transitions observed in atomic nuclei.

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.

Antidepressants that increase mitochondrial energetics may elevate risk of treatment-emergent mania

Preclinical evidence suggests that antidepressants (ADs) may differentially influence mitochondrial energetics. This study was conducted to investigate the relationship between mitochondrial function and illness vulnerability in bipolar disorder (BD), specifically risk of treatment-emergent mania (TEM). Participants with BD already clinically phenotyped as TEM+ (n = 176) or TEM- (n = 516) were further classified whether the TEM associated AD, based on preclinical studies, increased (Mito+, n = 600) or decreased (Mito-, n = 289) mitochondrial electron transport chain (ETC) activity. Comparison of TEM+ rates between Mito+ and Mito- ADs was performed using generalized estimating equations to account for participants exposed to multiple ADs while adjusting for sex, age at time of enrollment into the biobank and BD type (BD-I/schizoaffective vs. BD-II). A total of 692 subjects (62.7% female, 91.4% White, mean age 43.0 ± 14.0 years) including 176 cases (25.3%) of TEM+ and 516 cases (74.7%) of TEM- with previous exposure to Mito+ and/or Mito- antidepressants were identified. Adjusting for age, sex and BD subtype, TEM+ was more frequent with antidepressants that increased (24.7%), versus decreased (13.5%) mitochondrial energetics (OR = 2.21; p = 0.000009). Our preliminary retrospective data suggests there may be merit in reconceptualizing AD classification, not solely based on monoaminergic conventional drug mechanism of action, but additionally based on mitochondrial energetics. Future prospective clinical studies on specific antidepressants and mitochondrial activity are encouraged. Recognizing pharmacogenomic investigation of drug response may extend or overlap to genomics of disease risk, future studies should investigate potential interactions between mitochondrial mechanisms of disease risk and drug response.

Comorbidity and healthcare utilization in patients with treatment resistant depression: A large-scale retrospective cohort analysis using electronic health records

BACKGROUND

Medical comorbidity and healthcare utilization in patients with treatment resistant depression (TRD) is usually reported in convenience samples, making estimates unreliable. There is only limited large-scale clinical research on comorbidities and healthcare utilization in TRD patients.

METHODS

Electronic Health Record data from over 3.3 million patients from the INSIGHT Clinical Research Network in New York City was used to define TRD as initiation of a third antidepressant regimen in a 12-month period among patients diagnosed with major depressive disorder (MDD). Age and sex matched TRD and non-TRD MDD patients were compared for anxiety disorder, 27 comorbid medical conditions, and healthcare utilization.

RESULTS

Out of 30,218 individuals diagnosed with MDD, 15.2 % of patients met the criteria for TRD (n = 4605). Compared to MDD patients without TRD, the TRD patients had higher rates of anxiety disorder and physical comorbidities. They also had higher odds of ischemic heart disease (OR = 1.38), stroke/transient ischemic attack (OR = 1.57), chronic kidney diseases (OR = 1.53), arthritis (OR = 1.52), hip/pelvic fractures (OR = 2.14), and cancers (OR = 1.41). As compared to non-TRD MDD, TRD patients had higher rates of emergency room visits, and inpatient stays. In relation to patients without MDD, both TRD and non-TRD MDD patients had significantly higher levels of anxiety disorder and physical comorbidities.

LIMITATIONS

The INSIGHT-CRN data lack information on depression severity and medication adherence.

CONCLUSIONS

TRD patients compared to non-TRD MDD patients have a substantially higher prevalence of various psychiatric and medical comorbidities and higher health care utilization. These findings highlight the challenges of developing interventions and care coordination strategies to meet the complex clinical needs of TRD patients.

Genetic predisposition to major depressive disorder differentially impacts alcohol consumption and high-risk drinking situations in men and women with alcohol use disorder

Lifetime history of major depressive disorder (MDD) has a sex-specific association with pretreatment alcohol consumption in patients with alcohol dependence. Here, we investigated the association of genetic load for MDD estimated using a polygenic risk score (PRS) with pretreatment alcohol consumption assessed with Timeline Follow Back in a sample of 287 men and 156 women meeting DSM-IV-TR criteria for alcohol dependence. Preferred drinking situations were assessed using the Inventory of Drug Taking Situations (IDTS). Linear models were used to test for association of normalized alcohol consumption measures with the MDD-PRS, adjusting for ancestry, age, sex, and number of days sober at baseline. We fit models both with and without adjustment for MDD history and alcohol-use-related PRSs as covariates. Higher MDD-PRS was associated with lower 90-day total alcohol consumption in men (β = -0.16, p = 0.0012) but not in women (β = 0.11, p = 0.18). The association of MDD-PRS with IDTS measures was also sex-specific: higher MDD-PRS was associated with higher propensity to drink in temptation-related situations in women, while the opposite (negative association)was found in men. MDD-PRS was not associated with lifetime MDD history in our sample, and adjustment for lifetime MDD and alcohol-related PRSs did not impact the results. Our results suggest that genetic load for MDD impacts pretreatment alcohol consumption in a sex-specific manner, which is similar to, but independent from, the effect of history of MDD. The clinical implications of these findings and contributing biological and psychological factors should be investigated in future studies.

Systematic review: Wearable remote monitoring to detect nonalcohol/nonnicotine-related substance use disorder symptoms

BACKGROUND AND OBJECTIVES

Substance use disorders (SUDs) are chronic relapsing diseases characterized by significant morbidity and mortality. Phenomenologically, patients with SUDs present with a repeating cycle of intoxication, withdrawal, and craving, significantly impacting their diagnosis and treatment. There is a need for better identification and monitoring of these disease states. Remote monitoring chronic illness with wearable devices offers a passive, unobtrusive, constant physiological data assessment. We evaluate the current evidence base for remote monitoring of nonalcohol, nonnicotine SUDs.

METHODS

We performed a systematic, comprehensive literature review and screened 1942 papers.

RESULTS

We found 15 studies that focused mainly on the intoxication stage of SUD. These studies used wearable sensors measuring several physiological parameters (ECG, HR, O2 , Accelerometer, EDA, temperature) and implemented study-specific algorithms to evaluate the data.

DISCUSSION AND CONCLUSIONS

Studies were extracted, organized, and analyzed based on the three SUD disease states. The sample sizes were relatively small, focused primarily on the intoxication stage, had low monitoring compliance, and required significant computational power preventing "real-time" results. Cardiovascular data was the most consistently valuable data in the predictive algorithms. This review demonstrates that there is currently insufficient evidence to support remote monitoring of SUDs through wearable devices.

SCIENTIFIC SIGNIFICANCE

This is the first systematic review to show the available data on wearable remote monitoring of SUD symptoms in each stage of the disease cycle. This clinically relevant approach demonstrates what we know and do not know about the remote monitoring of SUDs within disease states.

Single cell transcriptomics reveals distinct transcriptional responses to oxycodone and buprenorphine by iPSC-derived brain organoids from patients with opioid use disorder

The opioid epidemic represents a national crisis. Oxycodone is one of the most prescribed opioid medications in the United States, whereas buprenorphine is currently the most prescribed medication for opioid use disorder (OUD) pharmacotherapy. Given the extensive use of prescription opioids and the global opioid epidemic, it is essential to understand how opioids modulate brain cell type function at the single-cell level. We performed single nucleus RNA-seq (snRNA-seq) using iPSC-derived forebrain organoids from three male OUD subjects in response to oxycodone, buprenorphine, or vehicle for seven days. We utilized the snRNA-seq data to identify differentially expressed genes following drug treatment using the Seurat integrative analysis pipeline. We utilized iPSC-derived forebrain organoids and single-cell sequencing technology as an unbiased tool to study cell-type-specific and drug-specific transcriptional responses. After quality control filtering, we analyzed 25787 cells and identified sixteen clusters using unsupervised clustering analysis. Our results reveal distinct transcriptional responses to oxycodone and buprenorphine by iPSC-derived brain organoids from patients with OUD. Specifically, buprenorphine displayed a significant influence on transcription regulation in glial cells. However, oxycodone induced type I interferon signaling in many cell types, including neural cells in brain organoids. Finally, we demonstrate that oxycodone, but not buprenorphine activated STAT1 and induced the type I interferon signaling in patients with OUD. These data suggest that elevation of STAT1 expression associated with OUD might play a role in transcriptional regulation in response to oxycodone. In summary, our results provide novel mechanistic insight into drug action at single-cell resolution.

Recent advances and challenges of rare variant association analysis in the biobank sequencing era

Causal variants for rare genetic diseases are often rare in the general population. Rare variants may also contribute to common complex traits and can have much larger per-allele effect sizes than common variants, although power to detect these associations can be limited. Sequencing costs have steadily declined with technological advancements, making it feasible to adopt whole-exome and whole-genome profiling for large biobank-scale sample sizes. These large amounts of sequencing data provide both opportunities and challenges for rare-variant association analysis. Herein, we review the basic concepts of rare-variant analysis methods, the current state-of-the-art methods in utilizing variant annotations or external controls to improve the statistical power, and particular challenges facing rare variant analysis such as accounting for population structure, extremely unbalanced case-control design. We also review recent advances and challenges in rare variant analysis for familial sequencing data and for more complex phenotypes such as survival data. Finally, we discuss other potential directions for further methodology investigation.

Identifying High-Risk Comorbidities Associated with Opioid Use Patterns Using Electronic Health Record Prescription Data

Introduction

Opioid use disorders (OUDs) constitute a major public health issue, and we urgently need alternative methods for characterizing risk for OUD. Electronic health records (EHRs) are useful tools for understanding complex medical phenotypes but have been underutilized for OUD because of challenges related to underdiagnosis, binary diagnostic frameworks, and minimally characterized reference groups. As a first step in addressing these challenges, a new paradigm is warranted that characterizes risk for opioid prescription misuse on a continuous scale of severity, i.e., as a continuum.

Methods

Across sites within the PsycheMERGE network, we extracted prescription opioid data and diagnoses that co-occur with OUD (including psychiatric and substance use disorders, pain-related diagnoses, HIV, and hepatitis C) for over 2.6 million patients across three health registries (Vanderbilt University Medical Center, Mass General Brigham, Geisinger) between 2005 and 2018. We defined three groups based on levels of opioid exposure: no prescriptions, minimal exposure, and chronic exposure and then compared the comorbidity profiles of these groups to the full registries and to those with OUD diagnostic codes.

Results

Our results confirm that EHR data reflects known higher prevalence of substance use disorders, psychiatric disorders, medical, and pain diagnoses in patients with OUD diagnoses and chronic opioid use. Comorbidity profiles that distinguish opioid exposure are strikingly consistent across large health systems, indicating the phenotypes described in this new quantitative framework are robust to health systems differences.

Conclusion

This work indicates that EHR prescription opioid data can serve as a platform to characterize complex risk markers for OUD using existing data.

Plasma TNFSF10 levels associated with acamprosate treatment response in patients with alcohol use disorder

Acamprosate is an anti-craving drug used in alcohol use disorder (AUD) pharmacotherapy. However, only a subset of patients achieves optimal treatment outcomes. The identification of predictive biomarkers of acamprosate treatment response in patients with AUD would be a substantial advance in addiction medicine. We designed this study to use proteomics data as a quantitative biological trait as a step toward identifying inflammatory modulators that might be associated with acamprosate treatment outcomes. The NIAAA-funded Mayo Clinic Center for the Individualized Treatment of Alcoholism study had previously recruited 442 AUD patients who received 3 months of acamprosate treatment. However, only 267 subjects returned for the 3-month follow-up visit and, as a result, had treatment outcome information available. Baseline alcohol craving intensity was the most significant predictor of acamprosate treatment outcomes. We performed plasma proteomics using the Olink target 96 inflammation panel and identified that baseline plasma TNF superfamily member 10 (TNFSF10) concentration was associated with alcohol craving intensity and variation in acamprosate treatment outcomes among AUD patients. We also performed RNA sequencing using baseline peripheral blood mononuclear cells from AUD patients with known acamprosate treatment outcomes which revealed that inflammation-related pathways were highly associated with relapse to alcohol use during the 3 months of acamprosate treatment. These observations represent an important step toward advancing our understanding of the pathophysiology of AUD and molecular mechanisms associated with acamprosate treatment response. In conclusion, applying omics-based approaches may be a practical approach for identifying biologic markers that could potentially predict alcohol craving intensity and acamprosate treatment response.

Electric Monopole Transition from the Superdeformed Band in ^{40}Ca

The electric monopole (E0) transition strength ρ^{2} for the transition connecting the third 0^{+} level, a "superdeformed" band head, to the "spherical" 0^{+} ground state in doubly magic ^{40}Ca is determined via e^{+}e^{-} pair-conversion spectroscopy. The measured value ρ^{2}(E0;0_{3}^{+}→0_{1}^{+})=2.3(5)×10^{-3} is the smallest ρ^{2}(E0;0^{+}→0^{+}) found in A<50 nuclei. In contrast, the E0 transition strength to the ground state observed from the second 0^{+} state, a band head of "normal" deformation, is an order of magnitude larger ρ^{2}(E0;0_{2}^{+}→0_{1}^{+})=25.9(16)×10^{-3}, which shows significant mixing between these two states. Large-scale shell-model (LSSM) calculations are performed to understand the microscopic structure of the excited states and the configuration mixing between them; experimental ρ^{2} values in ^{40}Ca and neighboring isotopes are well reproduced by the LSSM calculations. The unusually small ρ^{2}(E0;0_{3}^{+}→0_{1}^{+}) value is due to destructive interference in the mixing of shape-coexisting structures, which are based on several different multiparticle-multihole excitations. This observation goes beyond the usual treatment of E0 strengths, where two-state shape mixing cannot result in destructive interference.

Improving lithium dose prediction using population pharmacokinetics and pharmacogenomics: a cohort genome-wide association study in Sweden

BACKGROUND

Lithium is the most effective treatment for bipolar disorder, resulting in strong suicide prevention effects. The therapeutic range of lithium, however, is narrow and treatment initiation requires individual titration to address inter-individual variability. We aimed to improve lithium dose prediction using clinical and genomic data.

METHODS

We performed a population pharmacokinetic study followed by a genome-wide association study (GWAS), including two clinical Swedish cohorts. Participants in cohort 1 were from specialised outpatient clinics at Huddinge Hospital, in Stockholm, Sweden, and participants in cohort 2 were identified using the Swedish National Quality Registry for Bipolar disorder (BipoläR). Patients who received a lithium dose corresponding to at least one tablet of lithium sulphate (6 mmol) per day and had clinically relevant plasma concentrations of lithium were included in the study. Data on age, sex, bodyweight, height, creatinine concentration, estimated glomerular filtration rate (eGFR), lithium preparation, number of tablets of lithium per day, serum lithium concentration, and medications affecting kidney function (C09 antihypertensives, C03 [except C03D] sodium-retaining diuretics, and non-steroidal anti-inflammatory drugs) were obtained retrospectively for several timepoints when possible from electronic health records, BipoläR, and the Swedish prescription registry. The median time between timepoints was 1·07 years for cohort 1 and 1·09 years for cohort 2. The primary outcome of interest was the natural logarithm of total body clearance for lithium (CL_Li) associated with the clinical variables. The residual effects after accounting for age and sex, representing the individual-level effects (CL_Li,age/sex), were used as the dependent variable in a GWAS.

FINDINGS

2357 patients who were administered lithium (1423 women [60·4%] and 934 men [39·6%]; mean age 53·6 years [range 17-89], mainly of European descent) were included and 5627 data points were obtained. Age (variance explained [R²]: R²_cohort1=0·41 and R²_cohort2=0·31; both p<0·0001), sex (R²_cohort1=0·0063 [p=0·045] and R²_cohort2=0·026 [p<0·0001]), eGFR (R²_cohort1=0·38 and R²_cohort2=0·20; both p<0·0001), comedication with diuretics (R²_cohort1=0·0058 [p=0·014] and R²_cohort2=0·0026 [p<0·0001]), and agents acting on the renin-aldosterone-angiotensin system (R²_cohort1=0·028 and R²_cohort2=0·015; both p<0·0001) were clinical predictors of CL_Li. Notably, an association between CL_Li and serum lithium was observed, with a lower CL_Li being associated with higher serum lithium (R²_cohort1=0·13 and R²_cohort2=0·15; both p<0·0001). In a GWAS of CL_Li,age/sex, one locus was associated with a change in CL_Li (rs583503; β=-0·053 [95% CI -0·071 to -0·034]; p<0·00000005). We also found enrichment of the associations with genes expressed in the medulla (p=0·0014, corrected FDR=0·04) and cortex of the kidney (p=0·0015, corrected FDR=0·04), as well as associations with polygenic risk scores for eGFR (p value threshold: 0·05, p=0·01), body-mass index (p value threshold: 0·05, p=0·00025), and blood urea nitrogen (p value threshold: 0·001, p=0·00043). The model based on six clinical predictors explained 61·4% of the variance in CL_Li in cohort 1 and 49·8% in cohort 2. Adding genetic markers did not lead to major improvement of the models: within the subsample of genotyped individuals, the variance explained only increased from 59·32% to 59·36% in cohort 1 and from 49·21% to 50·03% in cohort 2 when including rs583503 and the four first principal components.

INTERPRETATION

Our model predictors could be used clinically to better guide lithium dosage, shortening the time to reach therapeutic concentrations, thus improving care. Identification of the first genomic locus and PRS to be associated with CL_Li introduces the opportunity of individualised medicine in lithium treatment.

FUNDING

Stanley Medical Research Institute, Swedish Research Council, Swedish Foundation for Strategic Research, Swedish Brain Foundation, Swedish Research Council, Söderström-Königska Foundation, Bror Gadelius Minnesfond, Swedish Mental Health Fund, Karolinska Institutet and Hospital.

Clinical Phenotype of Tardive Dyskinesia in Bipolar Disorder

PURPOSE

Recognizing the negative impact that antipsychotic-induced movement disorders have on the quality of life and treatment outcomes in bipolar disorder (BD), this study aimed to assess clinical correlates and antipsychotic use patterns of tardive dyskinesia (TD+) in BD.

MATERIALS AND METHODS

Participants with and without TD were included. Clinical variables were compared using t-test and χ2 test. Antipsychotic use patterns in TD+, including number of trials, mean doses, and estimated cumulative exposure, were assessed in a case-only analysis.

RESULTS

The prevalence rate of TD was 5.1%. In comparison to the TD- group (n = 1074), TD+ participants (n = 58) were older, more likely to be female and have type I bipolar illness. There were 60.3% of the TD+ group that continued using antipsychotics at study entry and had a mean cumulative exposure to antipsychotics of 18.2 ± 15.6 years. Average dose, in haloperidol equivalents, was 5.9 ± 3.5 mg and 77.7% of the trials were second-generation antipsychotics.

CONCLUSIONS

This study confirms previously identified TD risk factors, such as age, sex, and bipolar subtype in a large BD cohort. Limitations included a cross-sectional design and the lack of tardive illness severity assessment. As atypical antipsychotics continue to be primary mood stabilization treatment, attempting to harmonize large data sets to identify additional biomarkers of tardive risk will optimize individualized care for patients with BD.

Clinical and Genetic Correlates of Bipolar Disorder With Childhood-Onset Attention Deficit Disorder

BACKGROUND

Bipolar disorder (BD) with co-occurring attention deficit-hyperactivity disorder (ADHD) is associated with an unfavorable course of illness. We aimed to identify potential clinical and genetic correlates of BD with and without ADHD.

METHODS

Among patients with BD (N = 2,198) enrolled in the Mayo Clinic Bipolar Biobank we identified those with ADHD diagnosed in childhood (BD+cADHD; N = 350), those with adult-onset attention deficit symptoms (BD+aAD; N = 254), and those without ADHD (N = 1,594). We compared the groups using linear or logistic regression adjusting for age, sex, and recruitment site. For genotyped patients (N = 1,443), logistic regression was used to compare ADHD and BD polygenic risk scores (PRSs) between the BD groups, as well as to non-BD controls (N = 777).

RESULTS

Compared to the non-ADHD BD group, BD+cADHD patients were younger, more often men and had a greater number of co-occurring anxiety and substance use disorders (all p < 0.001). Additionally, BD+cADHD patients had poorer responses to lithium and lamotrigine (p = 0.005 and p = 0.007, respectively). In PRS analyses, all BD patient subsets had greater genetic risk for BD and ADHD when compared to non-BD controls (p < 0.001 in all comparisons). BD+cADHD patients had a higher ADHD-PRS than non-ADHD BD patients (p = 0.012). However, BD+aAD patients showed no evidence of higher ADHD-PRS than non-ADHD BD patients (p = 0.38).

CONCLUSIONS

BD+cADHD was associated with a greater number of comorbidities and reduced response to mood stabilizing treatments. The higher ADHD PRS for the BD+cADHD group may reflect a greater influence of genetic factors on early presentation of ADHD symptoms.

Relative Age Effect on Problematic Alcohol Use in Adolescents

Objective: Relative age refers to a difference in age between peers in the same academic year. Although younger age of alcohol use is associated with a higher risk of lifetime problematic alcohol use, the potential effects of relative age are poorly understood. We hypothesized that a younger relative age would be associated with a younger chronological age of testing positive for alcohol in a medical setting. Methods: Problematic alcohol use was operationalized and identified as a positive alcohol test (PAT) in a medical setting. This was a retrospective population study of all 12 to 18-year-old residents (n = 4610) of Olmsted County, Minnesota (USA), who were tested for alcohol in a medical care setting from 1998 through 2016. Cox regression models examined the relationship between relative age and the age at testing positive for alcohol. Results: Relative age was not associated with age at first PAT. Results remained nonsignificant after stratifying by gender, and after adjusting for race, number of nonalcohol-related psychiatric comorbidities, and type of alcohol testing. Conclusions: The results did not support a relative age effect as a risk factor for alcohol use in adolescents in Olmsted County, Minnesota. These results contrast with findings from previous studies on this topic, which suggested older relative age increases risk of alcohol use in adolescence.

Prevalence of insomnia symptoms and associated risk factors in UK Biobank participants with hazardous alcohol use and major depression

INTRODUCTION

We aimed to examine the prevalence of insomnia symptoms (IS), sleep duration, and associated risk factors in participants with hazardous/harmful alcohol use (HAU), major depressive disorders (MDD), and HAU+MDD.

METHODS

Data from the UK Biobank (UKB) (n = 55,000) were utilized to categorize participants into those with MDD (n = 5612), HAU (n = 15,893), MDD+HAU (n = 3738), and controls (n = 29,511). We examined whether rates of IS and sleep duration differed among the groups and determined the clinical predictors of IS. Rates of IS and sleep duration were compared using regression analyses accounting for demographic (age, sex, ethnicity, Townsend deprivation index) and clinical (body mass index, neuroticism score, alcohol consumption) factors.

RESULTS

The unadjusted prevalence of IS was 26.5%, 27%, 39.5%, and 43% in control, HAU, MDD, and MDD+HAU categories respectively. Rates of IS in controls versus HAU and MDD versus MDD+HAU did not differ in unadjusted models (p = 0.45 and 0.075, respectively). Prevalence of IS differed in the four groups (p < 0.0001 for all pairwise comparisons) after adjusting for demographic confounders. After further adjustment for clinical factors, effect sizes were reduced, but pairwise comparisons remained significant. After adjusting for demographic and clinical factors, sleep duration did not differ among the groups. After accounting for diagnostic category and demographic/clinical factors, older age (OR=1.33 per 10 year increase; p < 0.0001), female sex (OR=1.39; p < 0.0001), obesity (OR=1.17 compared to normal; p < 0.0001), higher neuroticism score (OR=1.13; p < 0.0001), and alcohol consumption (OR=1.01 per serving increase; p < 0.0001) were associated with IS.

CONCLUSION

Sleep-related morbidity is the greatest in the MDD+HAU group, followed by the MDD group. Demographic and clinical characteristics explain some, but not all of the differences in the prevalence of IS in MDD±HAU. Genetic and other factors capable of influencing IS in those with MDD, HAU, and MDD+HAU merit future investigation.

Association of Attention-Deficit/Hyperactivity Disorder and Depression Polygenic Scores with Lithium Response: A Consortium for Lithium Genetics Study

Response to lithium varies widely between individuals with bipolar disorder (BD). Polygenic risk scores (PRSs) can uncover pharmacogenomics effects and may help predict drug response. Patients (N = 2,510) with BD were assessed for long-term lithium response in the Consortium on Lithium Genetics using the Retrospective Criteria of Long-Term Treatment Response in Research Subjects with Bipolar Disorder score. PRSs for attention-deficit/hyperactivity disorder (ADHD), major depressive disorder (MDD), and schizophrenia (SCZ) were computed using lassosum and in a model including all three PRSs and other covariates, and the PRS of ADHD (β = -0.14; 95% confidence interval [CI]: -0.24 to -0.03; p value = 0.010) and MDD (β = -0.16; 95% CI: -0.27 to -0.04; p value = 0.005) predicted worse quantitative lithium response. A higher SCZ PRS was associated with higher rates of medication nonadherence (OR = 1.61; 95% CI: 1.34-1.93; p value = 2e-7). This study indicates that genetic risk for ADHD and depression may influence lithium treatment response. Interestingly, a higher SCZ PRS was associated with poor adherence, which can negatively impact treatment response. Incorporating genetic risk of ADHD, depression, and SCZ in combination with clinical risk may lead to better clinical care for patients with BD.

Genetic contributions to alcohol use disorder treatment outcomes: a genome-wide pharmacogenomics study

Naltrexone can aid in reducing alcohol consumption, while acamprosate supports abstinence; however, not all patients with alcohol use disorder (AUD) benefit from these treatments. Here we present the first genome-wide association study of AUD treatment outcomes based on data from the COMBINE and PREDICT studies of acamprosate and naltrexone, and the Mayo Clinic CITA study of acamprosate. Primary analyses focused on treatment outcomes regardless of pharmacological intervention and were followed by drug-stratified analyses to identify treatment-specific pharmacogenomic predictors of acamprosate and naltrexone response. Treatment outcomes were defined as: (1) time until relapse to any drinking (TR) and (2) time until relapse to heavy drinking (THR; ≥ 5 drinks for men, ≥4 drinks for women in a day), during the first 3 months of treatment. Analyses were performed within each dataset, followed by meta-analysis across the studies (N = 1083 European ancestry participants). Single nucleotide polymorphisms (SNPs) in the BRE gene were associated with THR (min p = 1.6E-8) in the entire sample, while two intergenic SNPs were associated with medication-specific outcomes (naltrexone THR: rs12749274, p = 3.9E-8; acamprosate TR: rs77583603, p = 3.1E-9). The top association signal for TR (p = 7.7E-8) and second strongest signal in the THR (p = 6.1E-8) analysis of naltrexone-treated patients maps to PTPRD, a gene previously implicated in addiction phenotypes in human and animal studies. Leave-one-out polygenic risk score analyses showed significant associations with TR (p = 3.7E-4) and THR (p = 2.6E-4). This study provides the first evidence of a polygenic effect on AUD treatment response, and identifies genetic variants associated with potentially medication-specific effects on AUD treatment response.

Genetic contributions to bipolar disorder: current status and future directions

Bipolar disorder (BD) is a highly heritable mental disorder and is estimated to affect about 50 million people worldwide. Our understanding of the genetic etiology of BD has greatly increased in recent years with advances in technology and methodology as well as the adoption of international consortiums and large population-based biobanks. It is clear that BD is also highly heterogeneous and polygenic and shows substantial genetic overlap with other psychiatric disorders. Genetic studies of BD suggest that the number of associated loci is expected to substantially increase in larger future studies and with it, improved genetic prediction of the disorder. Still, a number of challenges remain to fully characterize the genetic architecture of BD. First among these is the need to incorporate ancestrally-diverse samples to move research away from a Eurocentric bias that has the potential to exacerbate health disparities already seen in BD. Furthermore, incorporation of population biobanks, registry data, and electronic health records will be required to increase the sample size necessary for continued genetic discovery, while increased deep phenotyping is necessary to elucidate subtypes within BD. Lastly, the role of rare variation in BD remains to be determined. Meeting these challenges will enable improved identification of causal variants for the disorder and also allow for equitable future clinical applications of both genetic risk prediction and therapeutic interventions.

Pain Catastrophizing and Pain Self-Efficacy Mediate Interdisciplinary Pain Rehabilitation Program Outcomes at Posttreatment and Follow-Up

BACKGROUND

Decreasing pain catastrophizing and improving self-efficacy to self-manage chronic pain symptoms are important treatment targets in the context of interdisciplinary cognitive behavioral therapy for chronic pain. Greater pain catastrophizing has been shown to be associated with greater impact of pain symptoms on functioning, while conversely, greater pain self-efficacy has been associated with lower pain intensity and lower levels of disability.

OBJECTIVE

To prospectively evaluate interdisciplinary cognitive behavioral therapy for pain outcomes, as well as to the mediating effects of both pain catastrophizing and pain self-efficacy on outcome.

METHODS

Participants were 315 patients with chronic pain between April 2017 and April 2018 who completed a three-week interdisciplinary pain rehabilitation program. Pain severity, pain interference, pain catastrophizing, pain self-efficacy, quality of life, and depressive symptom questionnaires, and measures of physical performance were assessed at pre- and posttreatment. Follow-up questionnaires were returned by 163 participants. Effect size and reliable change analyses were conducted from pre- to posttreatment and pretreatment to 6-month follow-up. Mediation analyses were conducted to determine the mediating effect of pain catastrophizing and pain self-efficacy on pain outcome.

RESULTS

Significant improvements from pre- to posttreatment in pain outcomes were observed, and over 80% evidenced reliable change in at least one pain-relevant measure. Pain catastrophizing and pain self-efficacy mediated the relationship between changes in pain outcomes.

CONCLUSIONS

Interdisciplinary pain rehabilitation is an effective treatment and decreasing pain catastrophizing and increasing pain self-efficacy can influence maintenance of treatment gains.

Perinatal Mood and Anxiety Disorder Management in Multicenter Community Practices: Clinicians' Training, Current Practices and Perceived Strategies to Improve Future Implementation

Background:

This study aimed to explore clinicians’ perspectives on the current practice of perinatal mood and anxiety disorder (PMAD) management and strategies to improve future implementation.

Methods:

This study had a cross-sectional, descriptive design. A 35-item electronic survey was sent to clinicians (N = 118) who treated perinatal women and practiced at several community clinics at an academic medical center in the United States.

Results:

Among clinicians who provided care for perinatal women, 34.7% reported never receiving PMAD management training and 66.3% had less than 10 years of experience. Out of 10 patients who reported psychiatric symptoms, 47.8% of clinicians on average reported providing PMAD management to 1 to 3 patients and 40.7% noted that they conducted screening only when patient expresses PMAD symptoms. Suggested future improvements were providing training, developing a referral list, and establishing integrated behavioral health services.

Conclusions:

Results from this study indicated that while PMAD screening and management was implemented, improvements are warranted to meet established guidelines. Additionally, clinicians endorsed providing PMAD management to a small percentage of perinatal patients. Suggested strategies to increase adoption and implementation of PMAD management should be explored to improve access to behavioral health services for perinatal women.

BinomiRare: A robust test for association of a rare genetic variant with a binary outcome for mixed models and any case-control proportion

Whole-genome sequencing (WGS) and whole-exome sequencing studies have become increasingly available and are being used to identify rare genetic variants associated with health and disease outcomes. Investigators routinely use mixed models to account for genetic relatedness or other clustering variables (e.g., family or household) when testing genetic associations. However, no existing tests of the association of a rare variant with a binary outcome in the presence of correlated data control the type 1 error where there are (1) few individuals harboring the rare allele, (2) a small proportion of cases relative to controls, and (3) covariates to adjust for. Here, we address all three issues in developing a framework for testing rare variant association with a binary trait in individuals harboring at least one risk allele. In this framework, we estimate outcome probabilities under the null hypothesis and then use them, within the individuals with at least one risk allele, to test variant associations. We extend the BinomiRare test, which was previously proposed for independent observations, and develop the Conway-Maxwell-Poisson (CMP) test and study their properties in simulations. We show that the BinomiRare test always controls the type 1 error, while the CMP test sometimes does not. We then use the BinomiRare test to test the association of rare genetic variants in target genes with small-vessel disease (SVD) stroke, short sleep, and venous thromboembolism (VTE), in whole-genome sequence data from the Trans-Omics for Precision Medicine (TOPMed) program.

Feasibility and Initial Outcomes of a Group-Based Teletherapy Psychiatric Day Program for Adults With Serious Mental Illness: Open, Nonrandomized Trial in the Context of COVID-19

BACKGROUND

In the context of the COVID-19 pandemic, many behavioral health services have transitioned to teletherapy to continue delivering care for patients with mental illness. Studies that evaluate the outcome of this rapid teletherapy adoption and implementation are pertinent.

OBJECTIVE

This single-arm, nonrandomized pilot study aimed to assess the feasibility and initial patient-level outcomes of a psychiatric transitional day program that switched from an in-person group to a video teletherapy group during the COVID-19 pandemic.

METHODS

Patients with transdiagnostic conditions who were at risk of psychiatric hospitalization were referred to the Adult Transitions Program (ATP) at a large academic medical center in the United States. ATP was a 3-week intensive outpatient program that implemented group teletherapy guided by cognitive and behavioral principles delivered daily for 3 hours per day. Feasibility was assessed via retention, attendance rate, and rate of securing aftercare appointments prior to ATP discharge. Patients completed standardized patient-reported outcome measures at admission and discharge to assess the effectiveness of the program for improving quality of mental health, depression, anxiety, and suicide risk.

RESULTS

Patients (N=76) started the program between March and August of 2020. Feasibility was established, with 70 of the 76 patients (92%) completing the program and a mean attendance of 14.43 days (SD 1.22); also, 71 patients (95%) scheduled at least one behavioral health aftercare service prior to ATP discharge. All patient-level reported outcomes demonstrated significant improvements in depression (95% CI -3.6 to -6.2; Cohen d=0.77; P<.001), anxiety (95% CI -3.0 to -4.9; Cohen d=0.74; P<.001), overall suicide risk (95% CI -0.5 to -0.1; Cohen d=0.41; P=.02), wish to live (95% CI 0.3 to 1.0; Cohen d=0.39; P<.001), wish to die (95% CI -0.2 to -1.4; Cohen d=0.52; P=.01), and overall mental health (95% CI 1.5 to 4.5; Cohen d=0.39; P<.001) from admission to discharge.

CONCLUSIONS

Rapid adoption and implementation of a group-based teletherapy day program for adults at risk of psychiatric hospitalization appeared to be feasible and effective. Patients demonstrated high completion and attendance rates and reported significant improvements in psychosocial outcomes. Larger trials should be conducted to further evaluate the efficacy and effectiveness of the program through randomized controlled trials.

Increased Patient Safety-Related Incidents Following the Transition into Daylight Savings Time

BACKGROUND

"Spring forward," the start of daylight savings time (DST), reduces sleep opportunity by an hour. Insufficient sleep in healthcare workers resulting from the spring forward time change could potentially result in an increase in medical errors.

OBJECTIVE

We examined the change in reported patient safety-related incidents (SRIs), in the week following the transition into and out of DST over a period of 8 years.

DESIGN

Observational study SETTING: A US-based large healthcare organization with sites across multiple states MEASUREMENTS: Voluntarily reported SRIs that occurred 7 days prior to and following the spring and fall time changes for years 2010-2017 were ascertained. SRIs likely resulting from human error were identified separately. The changes in the number of SRIs (either all SRIs or SRIs restricted to those likely resulting from human error) from the week before and after the time change (either spring or fall) were modeled using a negative binomial mixed model with a random effect to correct for non-independent observations in consecutive weeks.

RESULTS

Over the 8-year period, we observed 4.2% (95% CI: - 1.1 to 9.7%; p = 0.12) and 8.8% (95% CI: - 2.5 to 21.5%; p = 0.13) increases in overall SRIs in the 7 days following DST when compared with 7 days prior for spring and fall, respectively. By restricting to SRIs likely resulting from human errors, we observed 18.7% (95% CI: 5.6 to 33.6%; p = 0.004) and 4.9% (95% CI: - 1.3 to 11.5%; p = 0.12) increases for spring and fall, respectively.

CONCLUSION

Policy makers and healthcare organizations should evaluate delayed start of shifts or other contingency measures to mitigate the increased risk of SRIs during transition to DST in spring.

TCF7L2 lncRNA: a link between bipolar disorder and body mass index through glucocorticoid signaling

Bipolar disorder (BD) and obesity are highly comorbid. We previously performed a genome-wide association study (GWAS) for BD risk accounting for the effect of body mass index (BMI), which identified a genome-wide significant single-nucleotide polymorphism (SNP) in the gene encoding the transcription factor 7 like 2 (TCF7L2). However, the molecular function of TCF7L2 in the central nervous system (CNS) and its possible role in the BD and BMI interaction remained unclear. In the present study, we demonstrated by studying human induced pluripotent stem cell (hiPSC)-derived astrocytes, cells that highly express TCF7L2 in the CNS, that the BD-BMI GWAS risk SNP is associated with glucocorticoid-dependent repression of the expression of a previously uncharacterized TCF7L2 transcript variant. That transcript is a long non-coding RNA (lncRNA-TCF7L2) that is highly expressed in the CNS but not in peripheral tissues such as the liver and pancreas that are involved in metabolism. In astrocytes, knockdown of the lncRNA-TCF7L2 resulted in decreased expression of the parent gene, TCF7L2, as well as alterations in the expression of a series of genes involved in insulin signaling and diabetes. We also studied the function of TCF7L2 in hiPSC-derived astrocytes by integrating RNA sequencing data after TCF7L2 knockdown with TCF7L2 chromatin-immunoprecipitation sequencing (ChIP-seq) data. Those studies showed that TCF7L2 directly regulated a series of BD risk genes. In summary, these results support the existence of a CNS-based mechanism underlying BD-BMI genetic risk, a mechanism based on a glucocorticoid-dependent expression quantitative trait locus that regulates the expression of a novel TCF7L2 non-coding transcript.

Mood-Stabilizing Antiepileptic Treatment Response in Bipolar Disorder: A Genome-Wide Association Study

Several antiepileptic drugs (AEDs) have US Food and Drug Administration (FDA) approval for use as mood stabilizers in bipolar disorder (BD), but not all BD patients respond to these AED mood stabilizers (AED‐MSs). To identify genetic polymorphisms that contribute to the variability in AED‐MS response, we performed a discovery genome‐wide association study (GWAS) of 199 BD patients from the Mayo Clinic Bipolar Disorder Biobank. Most of these patients had been treated with the AED‐MS valproate/divalproex and/or lamotrigine. AED‐MS response was assessed using the Alda scale, which quantifies clinical improvement while accounting for potential confounding factors. We identified two genome‐wide significant single‐nucleotide polymorphism (SNP) signals that mapped to the THSD7A (rs78835388, P = 7.1E‐09) and SLC35F3 (rs114872993, P = 3.2E‐08) genes. We also identified two genes with statistically significant gene‐level associations: ABCC1 (P = 6.7E‐07; top SNP rs875740, P = 2.0E‐6), and DISP1 (P = 8.9E‐07; top SNP rs34701716, P = 8.9E‐07). THSD7A SNPs were previously found to be associated with risk for several psychiatric disorders, including BD. Both THSD7A and SLC35F3 are expressed in excitatory/glutamatergic and inhibitory/γ‐aminobutyric acidergic (GABAergic) neurons, which are targets of AED‐MSs. ABCC1 is involved in the transport of valproate and lamotrigine metabolites, and the SNPs in ABCC1 and DISP1 with the strongest evidence of association in our GWAS are strong splicing quantitative trait loci in the human gut, suggesting a possible influence on drug absorption. In conclusion, our pharmacogenomic study identified novel genetic loci that appear to contribute to AED‐MS treatment response, and may facilitate precision medicine in BD.

A principal component approach to improve association testing with polygenic risk scores

Polygenic risk scores (PRSs) have become an increasingly popular approach for demonstrating polygenic influences on complex traits and for establishing common polygenic signals between different traits. PRSs are typically constructed using pruning and thresholding (P+T), but the best choice of parameters is uncertain; thus multiple settings are used and the best is chosen. Optimization can lead to inflated Type I error. Permutation procedures can correct this, but they can be computationally intensive. Alternatively, a single parameter setting can be chosen a priori for the PRS, but choosing suboptimal settings results in loss of power. We propose computing PRSs under a range of parameter settings, performing principal component analysis (PCA) on the resulting set of PRSs, and using the first PRS-PC in association tests. The first PC reweights the variants included in the PRS to achieve maximum variation over all PRS settings used. Using simulations and a real data application to study PRS association with bipolar disorder and psychosis in bipolar disorder, we compare the performance of the proposed PRS-PCA approach with a permutation test and an a priori selected p-value threshold. The PRS-PCA approach is simple to implement, outperforms the other strategies in most scenarios, and provides an unbiased estimate of prediction performance.

Potential pharmacogenomic targets in bipolar disorder: considerations for current testing and the development of decision support tools to individualize treatment selection

Background

Treatment in bipolar disorder (BD) is commonly applied as a multimodal therapy based on decision algorithms that lack an integrative understanding of molecular mechanisms or a biomarker associated clinical outcome measure. Pharmacogenetics/genomics study the individual genetic variation associated with drug response. This selective review of pharmacogenomics and pharmacogenomic testing (PGT) in BD will focus on candidate genes and genome wide association studies of pharmacokinetic drug metabolism and pharmacodynamic drug response/adverse event, and the potential role of decision support tools that incorporate multiple genotype/phenotype drug recommendations.

Main body

We searched PubMed from January 2013 to May 2019, to identify studies reporting on BD and pharmacogenetics, pharmacogenomics and PGT. Studies were selected considering their contribution to the field. We summarize our findings in: targeted candidate genes of pharmacokinetic and pharmacodynamic pathways, genome-wide association studies and, PGT platforms, related to BD treatment. This field has grown from studies of metabolizing enzymes (i.e., pharmacokinetics) and drug transporters (i.e., pharmacodynamics), to untargeted investigations across the entire genome with the potential to merge genomic data with additional biological information.

Conclusions

The complexity of BD genetics and, the heterogeneity in BD drug-related phenotypes, are important considerations for the design and interpretation of BD PGT. The clinical applicability of PGT in psychiatry is in its infancy and is far from reaching the robust impact it has in other medical disciplines. Nonetheless, promising findings are discovered with increasing frequency with remarkable relevance in neuroscience, pharmacology and biology.