Correlates of suicidal behaviors and genetic risk among United States veterans with schizophrenia or bipolar I disorder.

Persons diagnosed with schizophrenia (SCZ) or bipolar I disorder (BPI) are at high risk for self-injurious behavior, suicidal ideation, and suicidal behaviors (SB). Characterizing associations between diagnosed health problems, prior pharmacological treatments, and polygenic scores (PGS) has potential to inform risk stratification. We examined self-reported SB and ideation using the Columbia Suicide Severity Rating Scale (C-SSRS) among 3,942 SCZ and 5,414 BPI patients receiving care within the Veterans Health Administration (VHA). These cross-sectional data were integrated with electronic health records (EHRs), and compared across lifetime diagnoses, treatment histories, follow-up screenings, and mortality data. PGS were constructed using available genomic data for related traits. Genome-wide association studies were performed to identify and prioritize specific loci. Only 20% of the veterans who reported SB had a corroborating ICD-9/10 EHR code. Among those without prior SB, more than 20% reported new-onset SB at follow-up. SB were associated with a range of additional clinical diagnoses, and with treatment with specific classes of psychotropic medications (e.g., antidepressants, antipsychotics, etc.). PGS for externalizing behaviors, smoking initiation, suicide attempt, and major depressive disorder were associated with SB. The GWAS for SB yielded no significant loci. Among individuals with a diagnosed mental illness, self-reported SB were strongly associated with clinical variables across several EHR domains. Analyses point to sequelae of substance-related and psychiatric comorbidities as strong correlates of prior and subsequent SB. Nonetheless, past SB was frequently not documented in health records, underscoring the value of regular screening with direct, in-person assessments, especially among high-risk individuals.

Metformin prescription for U.S. veterans with prediabetes, 2010-2019.

Affecting an estimated 88 million Americans, prediabetes increases the risk for developing type 2 diabetes mellitus (T2DM), and independently, cardiovascular disease, retinopathy, nephropathy, and neuropathy. Nevertheless, little is known about the use of metformin for diabetes prevention among patients in the Veterans Health Administration, the largest integrated healthcare system in the U.S. This is a retrospective observational cohort study of the proportion of Veterans with incident prediabetes who were prescribed metformin at the Veterans Health Administration from October 2010 to September 2019. Among 1,059,605 Veterans with incident prediabetes, 12,009 (1.1%) were prescribed metformin during an average 3.4 years of observation after diagnosis. Metformin prescribing was marginally higher (1.6%) among those with body mass index (BMI) ≥35 kg/m2, age <60 years, HbA1c≥6.0%, or those with a history of gestational diabetes, all subgroups at a higher risk for progression to T2DM. In a multivariable model, metformin was more likely to be prescribed for those with BMI ≥35 kg/m2 incidence rate ratio [IRR] 2.6 [95% confidence intervals (CI): 2.1-3.3], female sex IRR, 2.4 [95% CI: 1.8-3.3], HbA1c≥6% IRR, 1.93 [95% CI: 1.5-2.4], age <60 years IRR, 1.7 [95% CI: 1.3-2.3], hypertriglyceridemia IRR, 1.5 [95% CI: 1.2-1.9], hypertension IRR, 1.5 [95% CI: 1.1-2.1], Major Depressive Disorder IRR, 1.5 [95% CI: 1.1-2.0], or schizophrenia IRR, 2.1 [95% CI: 1.2-3.8]. Over 20% of Veterans with prediabetes attended a comprehensive structured lifestyle modification clinic or program. Among Veterans with prediabetes, metformin was prescribed to 1.1% overall, a proportion that marginally increased to 1.6% in the subset of individuals at highest risk for progression to T2DM.

Effectiveness of Nirmatrelvir-Ritonavir Against the Development of Post-COVID-19 Conditions Among U.S. Veterans : A Target Trial Emulation.

BACKGROUND: COVID-19 has been linked to the development of many post-COVID-19 conditions (PCCs) after acute infection. Limited information is available on the effectiveness of oral antivirals used to treat acute COVID-19 in preventing the development of PCCs. OBJECTIVE: To measure the effectiveness of outpatient treatment of COVID-19 with nirmatrelvir-ritonavir in preventing PCCs. DESIGN: Retrospective target trial emulation study comparing matched cohorts receiving nirmatrelvir-ritonavir versus no treatment. SETTING: Veterans Health Administration (VHA). PARTICIPANTS: Nonhospitalized veterans in VHA care who were at risk for severe COVID-19 and tested positive for SARS-CoV-2 during January through July 2022. INTERVENTION: Nirmatrelvir-ritonavir treatment for acute COVID-19. MEASUREMENTS: Cumulative incidence of 31 potential PCCs at 31 to 180 days after treatment or a matched index date, including cardiac, pulmonary, renal, thromboembolic, gastrointestinal, neurologic, mental health, musculoskeletal, endocrine, and general conditions and symptoms. RESULTS: Eighty-six percent of the participants were male, with a median age of 66 years, and 17.5% were unvaccinated. Baseline characteristics were well balanced between participants treated with nirmatrelvir-ritonavir and matched untreated comparators. No differences were observed between participants treated with nirmatrelvir-ritonavir (n = 9593) and their matched untreated comparators in the incidence of most PCCs examined individually or grouped by organ system, except for lower combined risk for venous thromboembolism and pulmonary embolism (subhazard ratio, 0.65 [95% CI, 0.44 to 0.97]; cumulative incidence difference, -0.29 percentage points [CI, -0.52 to -0.05 percentage points]). LIMITATIONS: Ascertainment of PCCs using International Classification of Diseases, 10th Revision, codes may be inaccurate. Evaluation of many outcomes could have resulted in spurious associations with combined thromboembolic events by chance. CONCLUSION: Out of 31 potential PCCs, only combined thromboembolic events seemed to be reduced by nirmatrelvir-ritonavir. PRIMARY FUNDING SOURCE: U.S. Department of Veterans Affairs.

Anti-SARS-CoV-2 Pharmacotherapies Among Nonhospitalized US Veterans, January 2022 to January 2023.

Importance: Several pharmacotherapies have been authorized to treat nonhospitalized persons with symptomatic COVID-19. Longitudinal information on the use of these therapies is needed. Objective: To analyze trends and factors associated with prescription of outpatient COVID-19 pharmacotherapies within the Veterans Health Administration (VHA). Design, Setting, and Participants: This cohort study evaluated nonhospitalized veterans in VHA care who tested positive for SARS-CoV-2 from January 2022 through January 2023 using VHA and linked Community Care and Medicare databases. Exposures: Demographic characteristics, underlying medical conditions, COVID-19 vaccination, and regional and local systems of care, including Veterans Integrated Services Networks (VISNs). Main Outcomes and Measures: Monthly receipt of any COVID-19 pharmacotherapy (nirmatrelvir-ritonavir, molnupiravir, sotrovimab, or bebtelovimab) was described. Multivariable logistic regression was used to identify factors independently associated with receipt of any vs no COVID-19 pharmacotherapy. Results: Among 285 710 veterans (median [IQR] age, 63.1 [49.9-73.7] years; 247 358 males [86.6%]; 28 444 Hispanic [10.0%]; 61 269 Black [21.4%] and 198 863 White [69.6%]) who tested positive for SARS-CoV-2 between January 2022 and January 2023, the proportion receiving any pharmacotherapy increased from 3285 of 102 343 veterans (3.2%) in January 2022 to 5180 of 21 688 veterans (23.9%) in August 2022. The proportion declined to 2194 of 10 551 veterans (20.8%) by January 2023. Across VISNs, the range in proportion of patients who tested positive who received nirmatrelvir-ritonavir or molnupiravir during January 2023 was 41 of 692 veterans (5.9%) to 106 of 494 veterans (21.4%) and 2.1% to 120 of 1074 veterans (11.1%), respectively. Veterans receiving any treatment were more likely to be older (adjusted odds ratio [aOR] for ages 65-74 vs 50-64 years, 1.18; 95% CI, 1.14-1.22; aOR for ages ≥75 vs 50-64 years, 1.19; 95% CI, 1.15-1.23) and have a higher Charlson Comorbidity Index score (aOR for CCI ≥6 vs 0, 1.52; 95% CI, 1.44-1.59). Compared with White veterans, Black veterans (aOR, 1.06; 95% CI, 1.02-1.09) were more likely to receive treatment, and compared with non-Hispanic veterans, Hispanic veterans (aOR 1.06; 95% CI, 1.01-1.11) were more likely to receive treatment. Conclusions And Relevance: This study found that prescription of outpatient COVID-19 pharmacotherapies in the VHA peaked in August 2022 and declined thereafter. There were large regional differences in patterns of nirmatrelvir-ritonavir and molnupiravir use.

Effectiveness of COVID-19 Treatment With Nirmatrelvir-Ritonavir or Molnupiravir Among U.S. Veterans: Target Trial Emulation Studies With One-Month and Six-Month Outcomes.

BACKGROUND: Information about the effectiveness of oral antivirals in preventing short- and long-term COVID-19-related outcomes in the setting of Omicron variant transmission and COVID-19 vaccination is limited. OBJECTIVE: To measure the effectiveness of nirmatrelvir-ritonavir and molnupiravir for outpatient treatment of COVID-19. DESIGN: Three retrospective target trial emulation studies comparing matched cohorts of nirmatrelvir-ritonavir versus no treatment, molnupiravir versus no treatment, and nirmatrelvir-ritonavir versus molnupiravir. SETTING: Veterans Health Administration (VHA). PARTICIPANTS: Nonhospitalized veterans in VHA care who were at risk for severe COVID-19 and tested positive for SARS-CoV-2 during January through July 2022. INTERVENTION: Nirmatrelvir-ritonavir or molnupiravir pharmacotherapy. MEASUREMENTS: Incidence of any hospitalization or all-cause mortality at 30 days and from 31 to 180 days. RESULTS: Eighty-seven percent of participants were male; the median age was 66 years, and 18% were unvaccinated. Compared with matched untreated control participants, those treated with nirmatrelvir-ritonavir (n = 9607) had lower 30-day risk for hospitalization (22.07 vs. 30.32 per 1000 participants; risk difference [RD], -8.25 [95% CI, -12.27 to -4.23] per 1000 participants) and death (1.25 vs. 5.47 per 1000 participants; RD, -4.22 [CI, -5.45 to -3.00] per 1000 participants). Among persons alive at day 31, reductions were seen in 31- to 180-day incidence of death (hazard ratio, 0.66 [CI, 0.49 to 0.89]) but not hospitalization (subhazard ratio, 0.90 [CI, 0.79 to 1.02]). Molnupiravir-treated participants (n = 3504) had lower 30-day and 31- to 180-day risks for death (3.14 vs. 13.56 per 1000 participants at 30 days; RD, -10.42 [CI, -13.49 to -7.35] per 1000 participants; hazard ratio at 31 to 180 days, 0.67 [CI, 0.48 to 0.95]) but not hospitalization. A difference in 30-day or 31- to 180-day risk for hospitalization or death was not observed between matched nirmatrelvir- or molnupiravir-treated participants. LIMITATION: The date of COVID-19 symptom onset for most veterans was unknown. CONCLUSION: Nirmatrelvir-ritonavir was effective in reducing 30-day hospitalization and death. Molnupiravir was associated with a benefit for 30-day mortality but not hospitalization. Further reductions in mortality from 31 to 180 days were observed with both antivirals. PRIMARY FUNDING SOURCE: U.S. Department of Veterans Affairs.

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

BACKGROUND: The Million Veteran Program (MVP) participants represent 100 years of US history, including significant social and demographic changes over time. Our study assessed two aspects of the MVP: (i) longitudinal changes in population diversity and (ii) how these changes can be accounted for in genome-wide association studies (GWAS). To investigate these aspects, we divided MVP participants into five birth cohorts (N-range = 123,888 [born from 1943 to 1947] to 136,699 [born from 1948 to 1953]). RESULTS: Ancestry groups were defined by (i) HARE (harmonized ancestry and race/ethnicity) and (ii) a random-forest clustering approach using the 1000 Genomes Project and the Human Genome Diversity Project (1kGP + HGDP) reference panels (77 world populations representing six continental groups). In these groups, we performed GWASs of height, a trait potentially affected by population stratification. Birth cohorts demonstrate important trends in ancestry diversity over time. More recent HARE-assigned Europeans, Africans, and Hispanics had lower European ancestry proportions than older birth cohorts (0.010 < Cohen's d < 0.259, p < 7.80 × 10-4). Conversely, HARE-assigned East Asians showed an increase in European ancestry proportion over time. In GWAS of height using HARE assignments, genomic inflation due to population stratification was prevalent across all birth cohorts (linkage disequilibrium score regression intercept = 1.08 ± 0.042). The 1kGP + HGDP-based ancestry assignment significantly reduced the population stratification (mean intercept reduction = 0.045 ± 0.007, p < 0.05) confounding in the GWAS statistics. CONCLUSIONS: This study provides a characterization of ancestry diversity of the MVP cohort over time and compares two strategies to infer genetically defined ancestry groups by assessing differences in controlling population stratification in genome-wide association studies.

Correlates of suicidal behaviors and genetic risk among United States veterans with schizophrenia or bipolar I disorder.

Objective: Persons diagnosed with schizophrenia (SCZ) or bipolar I disorder (BPI) are at high risk for self-injurious behavior, suicidal ideation, and suicidal behaviors (SB). Characterizing associations between diagnosed mental and physical health problems, prior pharmacological treatments, and aggregate genetic factors has potential to inform risk stratification and mitigation strategies. Methods: In this study of 3,942 SCZ and 5,414 BPI patients receiving VA care, self-reported SB and ideation were assessed using the Columbia Suicide Severity Rating Scale (C-SSRS). These cross-sectional data were integrated with electronic health records (EHR), and compared by lifetime diagnoses, treatment histories, follow-up screenings, and mortality data. Polygenic scores (PGS) for traits related to psychiatric disorders, substance use, and cognition were constructed using available genomic data, and exploratory genome-wide association studies were performed to identify and prioritize specific loci. Results: Only 20% of veterans who self-reported SB had a corroborating ICD-9/10 code in their EHR; and among those who denied prior behaviors, more than 20% reported new-onset SB at follow-up. SB were associated with a range of psychiatric and non-psychiatric diagnoses, and with treatment with specific classes of psychotropic medications (e.g., antidepressants, antipsychotics, etc.). PGS for externalizing behaviors, smoking, suicide attempt, and major depressive disorder were also associated with attempt and ideation. Conclusions: Among individuals with a diagnosed mental illness, a GWAS for SB did not yield any significant loci. Self-reported SB were strongly associated with clinical variables across several EHR domains. Overall, clinical and polygenic analyses point to sequelae of substance-use related behaviors and other psychiatric comorbidities as strong correlates of prior and subsequent SB. Nonetheless, past SB was frequently not documented in clinical settings, underscoring the value of regular screening based on direct, in-person assessments, especially among high-risk individuals.

Dome1-JAK-STAT signaling between parasite and host integrates vector immunity and development.

Ancestral signaling pathways serve critical roles in metazoan development, physiology, and immunity. We report an evolutionary interspecies communication pathway involving a central Ixodes scapularis tick receptor termed Dome1, which acquired a mammalian cytokine receptor motif exhibiting high affinity for interferon-gamma (IFN-γ). Host-derived IFN-γ facilitates Dome1-mediated activation of the Ixodes JAK-STAT pathway. This accelerates tick blood meal acquisition and development while upregulating antimicrobial components. The Dome1-JAK-STAT pathway, which exists in most Ixodid tick genomes, regulates the regeneration and proliferation of gut cells-including stem cells-and dictates metamorphosis through the Hedgehog and Notch-Delta networks, ultimately affecting Ixodes vectorial competence. We highlight the evolutionary dependence of I. scapularis on mammalian hosts through cross-species signaling mechanisms that dually influence arthropod immunity and development.

Effectiveness of COVID-19 treatment with nirmatrelvir-ritonavir or molnupiravir among U.S. Veterans: target trial emulation studies with one-month and six-month outcomes.

Background: Information about the effectiveness of oral antivirals in preventing short- and long-term COVID-19-related outcomes during the Omicron surge is limited. We sought to determine the effectiveness of nirmatrelvir-ritonavir and molnupiravir for the outpatient treatment of COVID-19. Methods: We conducted three retrospective target trial emulation studies comparing matched patient cohorts who received nirmatrelvir-ritonavir versus no treatment, molnupiravir versus no treatment, and nirmatrelvir-ritonavir versus molnupiravir in the Veterans Health Administration (VHA). Participants were Veterans in VHA care at risk for severe COVID-19 who tested positive for SARS-CoV-2 in the outpatient setting during January and February 2022. Primary outcomes included all-cause 30-day hospitalization or death and 31-180-day incidence of acute or long-term care admission, death, or post-COVID-19 conditions. For 30-day outcomes, we calculated unadjusted risk rates, risk differences, and risk ratios. For 31-180-day outcomes, we used unadjusted time-to-event analyses. Results: Participants were 90% male with median age 67 years and 26% unvaccinated. Compared to matched untreated controls, nirmatrelvir-ritonavir-treated participants (N=1,587) had a lower 30-day risk of hospitalization (27.10/1000 versus 41.06/1000, risk difference [RD] - 13.97, 95% CI -23.85 to -4.09) and death (3.15/1000 versus 14.86/1000, RD -11.71, 95% CI - 16.07 to -7.35). Among persons who were alive at day 31, further significant reductions in 31-180-day incidence of hospitalization (sub-hazard ratio 1.07, 95% CI 0.83 to 1.37) or death (hazard ratio 0.61, 95% CI 0.35 to 1.08) were not observed. Molnupiravir-treated participants aged ≥65 years (n=543) had a lower combined 30-day risk of hospitalization or death (55.25/1000 versus 82.35/1000, RD -27.10, 95% CI -50.63 to -3.58). A statistically significant difference in 30-day or 31-180-day risk of hospitalization or death was not observed between matched nirmatrelvir- or molnupiravir-treated participants. Incidence of most post-COVID conditions was similar across comparison groups. Conclusions: Nirmatrelvir-ritonavir was highly effective in preventing 30-day hospitalization and death. Short-term benefit from molnupiravir was observed in older groups. Significant reductions in adverse outcomes from 31-180 days were not observed with either antiviral.

Tick transmission of Borrelia burgdorferi to the murine host is not influenced by environmentally acquired midgut microbiota.

BACKGROUND: Ixodes scapularis is the predominant tick vector of Borrelia burgdorferi, the agent of Lyme disease, in the USA. Molecular interactions between the tick and B. burgdorferi orchestrate the migration of spirochetes from the midgut to the salivary glands-critical steps that precede transmission to the vertebrate host. Over the last decade, research efforts have invoked a potential role for the tick microbiome in modulating tick-pathogen interactions. RESULTS: Using multiple strategies to perturb the microbiome composition of B. burgdorferi-infected nymphal ticks, we observe that changes in the microbiome composition do not significantly influence B. burgdorferi migration from the midgut, invasion of salivary glands, or transmission to the murine host. We also show that within 24 and 48 h of the onset of tick feeding, B. burgdorferi spirochetes are within the peritrophic matrix and epithelial cells of the midgut in preparation for exit from the midgut. CONCLUSIONS: This study highlights two aspects of tick-spirochete interactions: (1) environmental bacteria associated with the tick do not influence spirochete transmission to the mammalian host and (2) the spirochete may utilize an intracellular exit route during migration from the midgut to the salivary glands, a strategy that may allow the spirochete to distance itself from microbiota in the midgut lumen effectively. This may explain in part, the inability of environment-acquired midgut microbiota to significantly influence spirochete transmission. Unraveling a molecular understanding of this exit strategy will be critical to gain new insights into the biology of the spirochete and the tick. Video Abstract.

Penetrance and Pleiotropy of Polygenic Risk Scores for Schizophrenia, Bipolar Disorder, and Depression Among Adults in the US Veterans Affairs Health Care System.

Importance: Serious mental illnesses, including schizophrenia, bipolar disorder, and depression, are heritable, highly multifactorial disorders and major causes of disability worldwide. Objective: To benchmark the penetrance of current neuropsychiatric polygenic risk scores (PRSs) in the Veterans Health Administration health care system and to explore associations between PRS and broad categories of human disease via phenome-wide association studies. Design, Setting, and Participants: Extensive Veterans Health Administration's electronic health records were assessed from October 1999 to January 2021, and an embedded cohort of 9378 individuals with confirmed diagnoses of schizophrenia or bipolar 1 disorder were found. The performance of schizophrenia, bipolar disorder, and major depression PRSs were compared in participants of African or European ancestry in the Million Veteran Program (approximately 400 000 individuals), and associations between PRSs and 1650 disease categories based on ICD-9/10 billing codes were explored. Last, genomic structural equation modeling was applied to derive novel PRSs indexing common and disorder-specific genetic factors. Analysis took place from January 2021 to January 2022. Main Outcomes and Measures: Diagnoses based on in-person structured clinical interviews were compared with ICD-9/10 billing codes. PRSs were constructed using summary statistics from genome-wide association studies of schizophrenia, bipolar disorder, and major depression. Results: Of 707 299 enrolled study participants, 459 667 were genotyped at the time of writing; 84 806 were of broadly African ancestry (mean [SD] age, 58 [12.1] years) and 314 909 were of broadly European ancestry (mean [SD] age, 66.4 [13.5] years). Among 9378 individuals with confirmed diagnoses of schizophrenia or bipolar 1 disorder, 8962 (95.6%) were correctly identified using ICD-9/10 codes (2 or more). Among those of European ancestry, PRSs were robustly associated with having received a diagnosis of schizophrenia (odds ratio [OR], 1.81 [95% CI, 1.76-1.87]; P < 10-257) or bipolar disorder (OR, 1.42 [95% CI, 1.39-1.44]; P < 10-295). Corresponding effect sizes in participants of African ancestry were considerably smaller for schizophrenia (OR, 1.35 [95% CI, 1.29-1.42]; P < 10-38) and bipolar disorder (OR, 1.16 [95% CI, 1.11-1.12]; P < 10-10). Neuropsychiatric PRSs were associated with increased risk for a range of psychiatric and physical health problems. Conclusions and Relevance: Using diagnoses confirmed by in-person structured clinical interviews and current neuropsychiatric PRSs, the validity of an electronic health records-based phenotyping approach in US veterans was demonstrated, highlighting the potential of PRSs for disentangling biological and mediated pleiotropy.

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

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

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

Importance: Sickle cell trait (SCT), defined as the presence of 1 hemoglobin beta sickle allele (rs334-T) and 1 normal beta allele, is prevalent in millions of people in the US, particularly in individuals of African and Hispanic ancestry. However, the association of SCT with COVID-19 is unclear. Objective: To assess the association of SCT with the prepandemic health conditions in participants of the Million Veteran Program (MVP) and to assess the severity and sequelae of COVID-19. Design, Setting, and Participants: COVID-19 clinical data include 2729 persons with SCT, of whom 353 had COVID-19, and 129 848 SCT-negative individuals, of whom 13 488 had COVID-19. Associations between SCT and COVID-19 outcomes were examined using firth regression. Analyses were performed by ancestry and adjusted for sex, age, age squared, and ancestral principal components to account for population stratification. Data for the study were collected between March 2020 and February 2021. Exposures: The hemoglobin beta S (HbS) allele (rs334-T). Main Outcomes and Measures: This study evaluated 4 COVID-19 outcomes derived from the World Health Organization severity scale and phenotypes derived from International Classification of Diseases codes in the electronic health records. Results: Of the 132 577 MVP participants with COVID-19 data, mean (SD) age at the index date was 64.8 (13.1) years. Sickle cell trait was present in 7.8% of individuals of African ancestry and associated with a history of chronic kidney disease, diabetic kidney disease, hypertensive kidney disease, pulmonary embolism, and cerebrovascular disease. Among the 4 clinical outcomes of COVID-19, SCT was associated with an increased COVID-19 mortality in individuals of African ancestry (n = 3749; odds ratio, 1.77; 95% CI, 1.13 to 2.77; P = .01). In the 60 days following COVID-19, SCT was associated with an increased incidence of acute kidney failure. A counterfactual mediation framework estimated that on average, 20.7% (95% CI, -3.8% to 56.0%) of the total effect of SCT on COVID-19 fatalities was due to acute kidney failure. Conclusions and Relevance: In this genetic association study, SCT was associated with preexisting kidney comorbidities, increased COVID-19 mortality, and kidney morbidity.

Polygenic Breast Cancer Risk for Women Veterans in the Million Veteran Program.

Accurate breast cancer (BC) risk assessment allows personalized screening and prevention. Prospective validation of prediction models is required before clinical application. Here, we evaluate clinical- and genetic-based BC prediction models in a prospective cohort of women from the Million Veteran Program. MATERIALS AND METHODS: Clinical BC risk prediction models were validated in combination with a genetic polygenic risk score of 313 (PRS313) single-nucleotide polymorphisms in genetic females without prior BC diagnosis (n = 35,130, mean age 49 years) with 30% non-Hispanic African ancestry (AA). Clinical risk models tested were Breast and Prostate Cancer Cohort Consortium, literature review, and Breast Cancer Risk Assessment Tool, and implemented with or without PRS313. Prediction accuracy and association with incident breast cancer was evaluated with area under the receiver operating characteristic curve (AUC), hazard ratios, and proportion with high absolute lifetime risk. RESULTS: Three hundred thirty-eight participants developed incident breast cancers with a median follow-up of 3.9 years (2.5 cases/1,000 person-years), with 196 incident cases in women of European ancestry and 112 incident cases in AA women. Individualized Coherent Absolute Risk Estimator-literature review in combination with PRS313 had an AUC of 0.708 (95% CI, 0.659 to 0.758) in women with European or non-African ancestries and 0.625 (0.539 to 0.711) in AA women. Breast Cancer Risk Assessment Tool with PRS313 had an AUC of 0.695 (0.62 to 0.729) in European or non-AA and 0.675 (0.626 to 0.723) in AA women. Incorporation of PRS313 with clinical models improved prediction in European but not in AA women. Models estimated up to 9% of European and 18% of AA women with absolute lifetime risk > 20%. CONCLUSION: Clinical and genetic BC risk models predict incident BC in a large prospective multiracial cohort; however, more work is needed to improve genetic risk estimation in AA women.

Cooperative Studies Program (CSP) #572: A Study of Serious Mental Illness in Veterans as a Pathway to personalized medicine in Schizophrenia and Bipolar Illness.

BACKGROUND: Personalization of psychiatric treatment includes treatment of symptoms, cognition and functional deficits, suicide, and medical co-morbidities. VA Collaborative Study 572 examined a large sample of male and female veterans with schizophrenia (n=3,942) and with bipolar disorder (n=5,414) with phenotyping and genomic analyses. We present the results to date and future directions. METHODS: All veterans received a structured diagnostic interview and assessments of suicidal ideation and behavior, PTSD, and health. Veterans with schizophrenia were assessed for negative symptoms and lifetime depression. All were assessed with a cognitive and functional capacity assessment. Data for genome wide association studies were collected. Controls came from the VA Million Veteran Program. RESULTS: Suicidal ideation or behavior was present in 66%. Cognitive and functional deficits were consistent with previous studies. 40% of the veterans with schizophrenia had a lifetime major depressive episode and PTSD was present in over 30%. Polygenic risk score (PRS) analyses indicated that cognitive and functional deficits overlapped with PRS for cognition, education, and intelligence in the general population and PRS for suicidal ideation and behavior correlated with previous PRS for depression and suicidal ideation and behavior, as did the PRS for PTSD. DISCUSSION: Results to date provide directions for personalization of treatment in SMI, veterans with SMI, and veterans in general. The results of the genomic analyses suggest that cognitive deficits in SMI may be associated with general population features. Upcoming genomic analyses will reexamine the issues above, as well as genomic factors associated with smoking, substance abuse, negative symptoms, and treatment response.

Genome-Wide Association Studies of Schizophrenia and Bipolar Disorder in a Diverse Cohort of US Veterans.

BACKGROUND: Schizophrenia (SCZ) and bipolar disorder (BIP) are debilitating neuropsychiatric disorders, collectively affecting 2% of the world's population. Recognizing the major impact of these psychiatric disorders on the psychosocial function of more than 200 000 US Veterans, the Department of Veterans Affairs (VA) recently completed genotyping of more than 8000 veterans with SCZ and BIP in the Cooperative Studies Program (CSP) #572. METHODS: We performed genome-wide association studies (GWAS) in CSP #572 and benchmarked the predictive value of polygenic risk scores (PRS) constructed from published findings. We combined our results with available summary statistics from several recent GWAS, realizing the largest and most diverse studies of these disorders to date. RESULTS: Our primary GWAS uncovered new associations between CHD7 variants and SCZ, and novel BIP associations with variants in Sortilin Related VPS10 Domain Containing Receptor 3 (SORCS3) and downstream of PCDH11X. Combining our results with published summary statistics for SCZ yielded 39 novel susceptibility loci including CRHR1, and we identified 10 additional findings for BIP (28 326 cases and 90 570 controls). PRS trained on published GWAS were significantly associated with case-control status among European American (P < 10-30) and African American (P < .0005) participants in CSP #572. CONCLUSIONS: We have demonstrated that published findings for SCZ and BIP are robustly generalizable to a diverse cohort of US veterans. Leveraging available summary statistics from GWAS of global populations, we report 52 new susceptibility loci and improved fine-mapping resolution for dozens of previously reported associations.

Suicidal ideation, behavior, and mortality in male and female US veterans with severe mental illness.

BACKGROUND: We compared male and female American veterans with schizophrenia or bipolar disorder regarding clinical characteristics associated with lifetime suicidal ideation and behavior. Subsequent mortality, including death by suicide, was also assessed. METHODS: Data from questionnaires and face-to-face evaluations were collected during 2011-2014 from 8,049 male and 1,290 female veterans with schizophrenia or bipolar disorder. In addition to comparing male-female characteristics, Cox regression models-adjusted for demographic information, medical-psychiatric comorbidities, and self-reported suicidal ideation and behavior-were used to examine gender differences in associations of putative risk factors with suicide-specific and all-cause mortality during up to six years of follow-up. RESULTS: Women overall were younger, more likely to report a history of suicidal behavior, less likely to be substance abusers, and had lower overall mortality during follow-up. Among women only, psychiatric comorbidity was paradoxically associated with lower all-cause mortality (hazard ratio [HR]=0.53, 95% CI, 0.29-0.96, p = 0.037 for 1 disorder vs. none; HR=0.44, 95% CI, 0.25-0.77, p = 0.004 for ≥2 disorders vs. none). Suicide-specific mortality involved relatively few events, but crude rates were an order of magnitude higher than in the U.S. general and overall veteran populations. LIMITATIONS: Incomplete cause-of-death information and low statistical power for male-female comparisons regarding mortality. CONCLUSIONS: Female veterans with SMI differed from females in the general population by having a higher risk of suicide attempts. They also had more lifetime suicide attempts than male veterans with same diagnoses. These differences should inform public policy and clinical planning.

Genome-wide association study of cognitive performance in U.S. veterans with schizophrenia or bipolar disorder.

Cognitive impairment is a frequent and serious problem in patients with various forms of severe mental illnesses (SMI), including schizophrenia (SZ) and bipolar disorder (BP). Recent research suggests genetic links to several cognitive phenotypes in both SMI and in the general population. Our goal in this study was to identify potential genomic signatures of cognitive functioning in veterans with severe mental illness and compare them to previous findings for cognition across different populations. Veterans Affairs (VA) Cooperative Studies Program (CSP) Study #572 evaluated cognitive and functional capacity measures among SZ and BP patients. In conjunction with the VA Million Veteran Program, 3,959 European American (1,095 SZ, 2,864 BP) and 2,601 African American (1,095 SZ, 2,864 BP) patients were genotyped using a custom Affymetrix Axiom Biobank array. We performed a genome-wide association study of global cognitive functioning, constructed polygenic scores for SZ and cognition in the general population, and examined genetic correlations with 2,626 UK Biobank traits. Although no single locus attained genome-wide significance, observed allelic effects were strongly consistent with previous studies. We observed robust associations between global cognitive functioning and polygenic scores for cognitive performance, intelligence, and SZ risk. We also identified significant genetic correlations with several cognition-related traits in UK Biobank. In a diverse cohort of U.S. veterans with SZ or BP, we demonstrate broad overlap of common genetic effects on cognition in the general population, and find that greater polygenic loading for SZ risk is associated with poorer cognitive performance.

Genome-wide association study of post-traumatic stress disorder reexperiencing symptoms in >165,000 US veterans.

Post-traumatic stress disorder (PTSD) is a major problem among military veterans and civilians alike, yet its pathophysiology remains poorly understood. We performed a genome-wide association study and bioinformatic analyses, which included 146,660 European Americans and 19,983 African Americans in the US Million Veteran Program, to identify genetic risk factors relevant to intrusive reexperiencing of trauma, which is the most characteristic symptom cluster of PTSD. In European Americans, eight distinct significant regions were identified. Three regions had values of P < 5 × 10-10: CAMKV; chromosome 17 closest to KANSL1, but within a large high linkage disequilibrium region that also includes CRHR1; and TCF4. Associations were enriched with respect to the transcriptomic profiles of striatal medium spiny neurons. No significant associations were observed in the African American cohort of the sample. Results in European Americans were replicated in the UK Biobank data. These results provide new insights into the biology of PTSD in a well-powered genome-wide association study.

Genome-wide Association Study of Maximum Habitual Alcohol Intake in >140,000 U.S. European and African American Veterans Yields Novel Risk Loci.

BACKGROUND: Habitual alcohol use can be an indicator of alcohol dependence, which is associated with a wide range of serious health problems. METHODS: We completed a genome-wide association study in 126,936 European American and 17,029 African American subjects in the Veterans Affairs Million Veteran Program for a quantitative phenotype based on maximum habitual alcohol consumption. RESULTS: ADH1B, on chromosome 4, was the lead locus for both populations: for the European American sample, rs1229984 (p = 4.9 × 10-47); for African American, rs2066702 (p = 2.3 × 10-12). In the European American sample, we identified three additional genome-wide-significant maximum habitual alcohol consumption loci: on chromosome 17, rs77804065 (p = 1.5 × 10-12), at CRHR1 (corticotropin-releasing hormone receptor 1); the protein product of this gene is involved in stress and immune responses; and on chromosomes 8 and 10. European American and African American samples were then meta-analyzed; the associated region at CRHR1 increased in significance to 1.02 × 10-13, and we identified two additional genome-wide significant loci, FGF14 (p = 9.86 × 10-9) (chromosome 13) and a locus on chromosome 11. Besides ADH1B, none of the five loci have prior genome-wide significant support. Post-genome-wide association study analysis identified genetic correlation to other alcohol-related traits, smoking-related traits, and many others. Replications were observed in UK Biobank data. Genetic correlation between maximum habitual alcohol consumption and alcohol dependence was 0.87 (p = 4.78 × 10-9). Enrichment for cell types included dopaminergic and gamma-aminobutyric acidergic neurons in midbrain, and pancreatic delta cells. CONCLUSIONS: The present study supports five novel alcohol-use risk loci, with particularly strong statistical support for CRHR1. Additionally, we provide novel insight regarding the biology of harmful alcohol use.