DNA methylation-based telomere length is associated with HIV infection, physical frailty, cancer, and all-cause mortality.

Telomere length (TL) is an important indicator of cellular aging. Shorter TL is associated with several age-related diseases including coronary heart disease, heart failure, diabetes, osteoporosis, and cancer. Recently, a DNA methylation-based TL (DNAmTL) estimator has been developed as an alternative method for directly measuring TL. In this study, we examined the association of DNAmTL with cancer prevalence and mortality risk among people with and without HIV in the Veterans Aging Cohort Study Biomarker Cohort (VACS, N = 1917) and Women's Interagency HIV Study Cohort (WIHS, N = 481). We profiled DNAm in whole blood (VACS) or in peripheral blood mononuclear cells (WIHS) using an array-based method. Cancer prevalence was estimated from electronic medical records and cancer registry data. The VACS Index was used as a measure of physiologic frailty. Models were adjusted for self-reported race and ethnicity, batch, smoking status, alcohol consumption, and five cell types (CD4, CD8, NK, B cell, and monocyte). We found that people with HIV had shorter average DNAmTL than those without HIV infection [beta = -0.25, 95% confidence interval (-0.32, -0.18), p = 1.48E-12]. Greater value of VACS Index [beta = -0.002 (-0.003, -0.001), p = 2.82E-05] and higher cancer prevalence [beta = -0.07 (-0.10, -0.03), p = 1.37E-04 without adjusting age] were associated with shortened DNAmTL. In addition, one kilobase decrease in DNAmTL was associated with a 40% increase in mortality risk [hazard ratio: 0.60 (0.44, 0.82), p = 1.42E-03]. In summary, HIV infection, physiologic frailty, and cancer are associated with shortening DNAmTL, contributing to an increased risk of all-cause mortality.

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.

Computationally inferred cell-type specific epigenome-wide DNA methylation analysis unveils distinct methylation patterns among immune cells for HIV infection in three cohorts.

BACKGROUND: Epigenome-wide association studies (EWAS) have identified CpG sites associated with HIV infection in blood cells in bulk, which offer limited knowledge of cell-type specific methylation patterns associated with HIV infection. In this study, we aim to identify differentially methylated CpG sites for HIV infection in immune cell types: CD4+ T-cells, CD8+ T-cells, B cells, Natural Killer (NK) cells, and monocytes. METHODS: Applying a computational deconvolution method, we performed a cell-type based EWAS for HIV infection in three independent cohorts (Ntotal = 1,382). DNA methylation in blood or in peripheral blood mononuclear cells (PBMCs) was profiled by an array-based method and then deconvoluted by Tensor Composition Analysis (TCA). The TCA-computed CpG methylation in each cell type was first benchmarked by bisulfite DNA methylation capture sequencing in a subset of the samples. Cell-type EWAS of HIV infection was performed in each cohort separately and a meta-EWAS was conducted followed by gene set enrichment analysis. RESULTS: The meta-analysis unveiled a total of 2,021 cell-type unique significant CpG sites for five inferred cell types. Among these inferred cell-type unique CpG sites, the concordance rate in the three cohorts ranged from 96% to 100% in each cell type. Cell-type level meta-EWAS unveiled distinct patterns of HIV-associated differential CpG methylation, where 74% of CpG sites were unique to individual cell types (false discovery rate, FDR <0.05). CD4+ T-cells had the largest number of unique HIV-associated CpG sites (N = 1,624) compared to any other cell type. Genes harboring significant CpG sites are involved in immunity and HIV pathogenesis (e.g. CD4+ T-cells: NLRC5, CX3CR1, B cells: IFI44L, NK cells: IL12R, monocytes: IRF7), and in oncogenesis (e.g. CD4+ T-cells: BCL family, PRDM16, monocytes: PRDM16, PDCD1LG2). HIV-associated CpG sites were enriched among genes involved in HIV pathogenesis and oncogenesis that were enriched among interferon-α and -γ, TNF-α, inflammatory response, and apoptotic pathways. CONCLUSION: Our findings uncovered computationally inferred cell-type specific modifications in the host epigenome for people with HIV that contribute to the growing body of evidence regarding HIV pathogenesis.

Postoperative Outcomes Associated with the Timing of Surgery After SARS-CoV-2 Infection.

OBJECTIVE: Examine the association between prior SARS-CoV-2 infection, interval from infection to surgery, and adverse surgical outcomes. SUMMARY BACKGROUND DATA: Earlier series have reported worse outcomes for surgery after COVID-19 illness, and these findings have led to routinely deferring surgery seven weeks after infection. METHODS: We created a retrospective cohort of patients from US Veterans Health Administration facilities nationwide, April 2020-September 2022, undergoing surgical procedures. Primary outcomes were 90-day all-cause mortality and 30-day complications. Within surgical procedure groupings, SARS-CoV-2 infected and uninfected patients were matched in a 1:4 ratio. We categorized patients by two-week intervals from SARS-CoV-2 positive test to surgery. Hierarchical multilevel multivariable logistic regression models were used to estimate the association between infection to surgery interval versus no infection and primary endpoints. RESULTS: We identified 82,815 veterans undergoing eligible operations (33% general, 27% orthopedic, 13% urologic, 9% vascular), of whom 16,563 (20%) had laboratory confirmed SARS-CoV-2 infection prior to surgery. The multivariable models demonstrated an association between prior SARS-CoV-2 infection and increased 90-day mortality (odds ratio (OR) 1.42, 95% CI 1.08, 1.86) and complications (OR 1.32, 95% CI 1.11, 1.57) only for patients having surgery within 14 days of infection. ASA-stratified multivariable models showed that the associations between increased 90-day mortality (OR 1.40, 95% CI 1.12, 1.75) and complications (OR 1.73, 95% CI 1.34, 2.24) for patients having surgery within 14 days of infection were confined to those with ASA 4-5. CONCLUSIONS: In a contemporary surgical cohort, patients with prior SARS-CoV-2 infection only had increased post-operative mortality or complications when they had surgery within 14 days after positive test. These findings support revising timing recommendations between surgery and prior SARS-CoV-2 infection.

Cardiovascular Disease Risk Factor Control in People With and Without HIV.

BACKGROUND: Management of hypertension, dyslipidemia, diabetes and other modifiable factors may mitigate the cardiovascular disease (CVD) risk in people with human immunodeficiency virus (HIV, PWH) compared with people without HIV (PWoH). METHODS: This was a retrospective cohort study of 8285 PWH and 170 517 PWoH from an integrated health system. Risk factor control was measured using a novel disease management index (DMI) accounting for amount/duration above treatment goals (0% to 100% [perfect control]), including 2 DMIs for hypertension (diastolic and systolic blood pressure), 3 for dyslipidemia (low-density lipoprotein, total cholesterol, triglycerides), and 1 for diabetes (HbA1c). CVD risk by HIV status was evaluated overall and in subgroups defined by DMIs, smoking, alcohol use, and overweight/obesity in adjusted Cox proportional hazards models. RESULTS: PWH and PWoH had similar DMIs (80%-100%) except for triglycerides (worse for PWH) and HbA1c (better for PWH). In adjusted models, PWH had an elevated risk of CVD compared with PWoH (hazard ratio [HR], 1.18; 95% confidence interval [CI], 1.07-1.31). This association was attenuated in subgroups with controlled dyslipidemia and diabetes but remained elevated for PWH with controlled hypertension or higher total cholesterol. The strongest HIV status association with CVD was seen in the subgroup with frequent unhealthy alcohol use (HR, 2.13; 95% CI, 1.04-4.34). CONCLUSIONS: Control of dyslipidemia and diabetes, but not hypertension, attenuated the HIV status association with CVD. The strong association of HIV and CVD with frequent unhealthy alcohol use suggests enhanced screening and treatment of alcohol problems in PWH is warranted.

Adaption and National Validation of a Tool for Predicting Mortality from Other Causes Among Men with Nonmetastatic Prostate Cancer.

BACKGROUND: An electronic health record-based tool could improve accuracy and eliminate bias in provider estimation of the risk of death from other causes among men with nonmetastatic cancer. OBJECTIVE: To recalibrate and validate the Veterans Aging Cohort Study Charlson Comorbidity Index (VACS-CCI) to predict non-prostate cancer mortality (non-PCM) and to compare it with a tool predicting prostate cancer mortality (PCM). DESIGN, SETTING, AND PARTICIPANTS: An observational cohort of men with biopsy-confirmed nonmetastatic prostate cancer, enrolled from 2001 to 2018 in the national US Veterans Health Administration (VA), was divided by the year of diagnosis into the development (2001-2006 and 2008-2018) and validation (2007) sets. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: Mortality (all cause, non-PCM, and PCM) was evaluated. Accuracy was assessed using calibration curves and C statistic in the development, validation, and combined sets; overall; and by age (<65 and 65+ yr), race (White and Black), Hispanic ethnicity, and treatment groups. RESULTS AND LIMITATIONS: Among 107 370 individuals, we observed 24 977 deaths (86% non-PCM). The median age was 65 yr, 4947 were Black, and 5010 were Hispanic. Compared with CCI and age alone (C statistic 0.67, 95% confidence interval [CI] 0.67-0.68), VACS-CCI demonstrated improved validated discrimination (C statistic 0.75, 95% CI 0.74-0.75 for non-PCM). The prostate cancer mortality tool also discriminated well in validation (C statistic 0.81, 95% CI 0.78-0.83). Both were well calibrated overall and within subgroups. Owing to missing data, 18 009/125 379 (14%) were excluded, and VACS-CCI should be validated outside the VA prior to outside application. CONCLUSIONS: VACS-CCI is ready for implementation within the VA. Electronic health record-assisted calculation is feasible, improves accuracy over age and CCI alone, and could mitigate inaccuracy and bias in provider estimation. PATIENT SUMMARY: Veterans Aging Cohort Study Charlson Comorbidity Index is ready for application within the Veterans Health Administration. Electronic health record-assisted calculation is feasible, improves accuracy over age and Charlson Comorbidity Index alone, and might help mitigate inaccuracy and bias in provider estimation of the risk of non-prostate cancer mortality.

The forecasted prevalence of comorbidities and multimorbidity in people with HIV in the United States through the year 2030: A modeling study.

BACKGROUND: Estimating the medical complexity of people aging with HIV can inform clinical programs and policy to meet future healthcare needs. The objective of our study was to forecast the prevalence of comorbidities and multimorbidity among people with HIV (PWH) using antiretroviral therapy (ART) in the United States (US) through 2030. METHODS AND FINDINGS: Using the PEARL model-an agent-based simulation of PWH who have initiated ART in the US-the prevalence of anxiety, depression, stage ≥3 chronic kidney disease (CKD), dyslipidemia, diabetes, hypertension, cancer, end-stage liver disease (ESLD), myocardial infarction (MI), and multimorbidity (≥2 mental or physical comorbidities, other than HIV) were forecasted through 2030. Simulations were informed by the US CDC HIV surveillance data of new HIV diagnosis and the longitudinal North American AIDS Cohort Collaboration on Research and Design (NA-ACCORD) data on risk of comorbidities from 2009 to 2017. The simulated population represented 15 subgroups of PWH including Hispanic, non-Hispanic White (White), and non-Hispanic Black/African American (Black/AA) men who have sex with men (MSM), men and women with history of injection drug use and heterosexual men and women. Simulations were replicated for 200 runs and forecasted outcomes are presented as median values (95% uncertainty ranges are presented in the Supporting information). In 2020, PEARL forecasted a median population of 670,000 individuals receiving ART in the US, of whom 9% men and 4% women with history of injection drug use, 60% MSM, 8% heterosexual men, and 19% heterosexual women. Additionally, 44% were Black/AA, 32% White, and 23% Hispanic. Along with a gradual rise in population size of PWH receiving ART-reaching 908,000 individuals by 2030-PEARL forecasted a surge in prevalence of most comorbidities to 2030. Depression and/or anxiety was high and increased from 60% in 2020 to 64% in 2030. Hypertension decreased while dyslipidemia, diabetes, CKD, and MI increased. There was little change in prevalence of cancer and ESLD. The forecasted multimorbidity among PWH receiving ART increased from 63% in 2020 to 70% in 2030. There was heterogeneity in trends across subgroups. Among Black women with history of injection drug use in 2030 (oldest demographic subgroup with median age of 66 year), dyslipidemia, CKD, hypertension, diabetes, anxiety, and depression were most prevalent, with 92% experiencing multimorbidity. Among Black MSM in 2030 (youngest demographic subgroup with median age of 42 year), depression and CKD were highly prevalent, with 57% experiencing multimorbidity. These results are limited by the assumption that trends in new HIV diagnoses, mortality, and comorbidity risk observed in 2009 to 2017 will persist through 2030; influences occurring outside this period are not accounted for in the forecasts. CONCLUSIONS: The PEARL forecasts suggest a continued rise in comorbidity and multimorbidity prevalence to 2030, marked by heterogeneities across race/ethnicity, gender, and HIV acquisition risk subgroups. HIV clinicians must stay current on the ever-changing comorbidities-specific guidelines to provide guideline-recommended care. HIV clinical directors should ensure linkages to subspecialty care within the clinic or by referral. HIV policy decision-makers must allocate resources and support extended clinical capacity to meet the healthcare needs of people aging with HIV.

High dimensional predictions of suicide risk in 4.2 million US Veterans using ensemble transfer learning.

We present an ensemble transfer learning method to predict suicide from Veterans Affairs (VA) electronic medical records (EMR). A diverse set of base models was trained to predict a binary outcome constructed from reported suicide, suicide attempt, and overdose diagnoses with varying choices of study design and prediction methodology. Each model used twenty cross-sectional and 190 longitudinal variables observed in eight time intervals covering 7.5 years prior to the time of prediction. Ensembles of seven base models were created and fine-tuned with ten variables expected to change with study design and outcome definition in order to predict suicide and combined outcome in a prospective cohort. The ensemble models achieved c-statistics of 0.73 on 2-year suicide risk and 0.83 on the combined outcome when predicting on a prospective cohort of [Formula: see text] 4.2 M veterans. The ensembles rely on nonlinear base models trained using a matched retrospective nested case-control (Rcc) study cohort and show good calibration across a diversity of subgroups, including risk strata, age, sex, race, and level of healthcare utilization. In addition, a linear Rcc base model provided a rich set of biological predictors, including indicators of suicide, substance use disorder, mental health diagnoses and treatments, hypoxia and vascular damage, and demographics.

Genome-wide association study of prostate-specific antigen levels in 392,522 men identifies new loci and improves cross-ancestry prediction.

We conducted a multi-ancestry genome-wide association study of prostate-specific antigen (PSA) levels in 296,754 men (211,342 European ancestry; 58,236 African ancestry; 23,546 Hispanic/Latino; 3,630 Asian ancestry; 96.5% of participants were from the Million Veteran Program). We identified 318 independent genome-wide significant (p≤5e-8) variants, 184 of which were novel. Most demonstrated evidence of replication in an independent cohort (n=95,768). Meta-analyzing discovery and replication (n=392,522) identified 447 variants, of which a further 111 were novel. Out-of-sample variance in PSA explained by our new polygenic risk score reached 16.9% (95% CI=16.1%-17.8%) in European ancestry, 9.5% (95% CI=7.0%-12.2%) in African ancestry, 18.6% (95% CI=15.8%-21.4%) in Hispanic/Latino, and 15.3% (95% CI=12.7%-18.1%) in Asian ancestry, and lower for higher age. Our study highlights how including proportionally more participants from underrepresented populations improves genetic prediction of PSA levels, with potential to personalize prostate cancer screening.

Integrase strand-transfer inhibitor use and cardiovascular events in adults with HIV: an emulation of target trials in the HIV-CAUSAL Collaboration and the Antiretroviral Therapy Cohort Collaboration.

BACKGROUND: A recent observational study suggested that the risk of cardiovascular events could be higher among antiretroviral therapy (ART)-naive individuals with HIV who receive integrase strand-transfer inhibitor (INSTI)-based ART than among those who receive other ART regimens. We aimed to emulate target trials separately in ART-naive and ART-experienced individuals with HIV to examine the effect of using INSTI-based regimens versus other ART regimens on the 4-year risk of cardiovascular events. METHODS: We used routinely recorded clinical data from 12 cohorts that collected information on cardiovascular events, BMI, and blood pressure from two international consortia of cohorts of people with HIV from Europe and North America. For the target trial in individuals who had previously never used ART (ie, ART-naive), eligibility criteria were aged 18 years or older, a detectable HIV-RNA measurement while ART-naive (>50 copies per mL), and no history of a cardiovascular event or cancer. Eligibility criteria for the target trial in those with previous use of non-INSTI-based ART (ie, ART-experienced) were the same except that individuals had to have been on at least one non-INSTI-based ART regimen and be virally suppressed (≤50 copies per mL). We assessed eligibility for both trials for each person-month between January, 2013, and January, 2023, and assigned individuals to the treatment strategy that was compatible with their data. We estimated the standardised 4-year risks of cardiovascular events (myocardial infarction, stroke, or invasive cardiovascular procedure) via pooled logistic regression models adjusting for time and baseline covariates. In per-protocol analyses, we censored individuals if they deviated from their assigned treatment strategy for more than 2 months and weighted uncensored individuals by the inverse of their time-varying probability of remaining uncensored. The denominator of the weight was estimated via a pooled logistic model that included baseline and time-varying covariates. FINDINGS: The analysis in ART-naive individuals included 10 767 INSTI initiators and 8292 non-initiators of INSTI. There were 43 cardiovascular events in INSTI initiators (median follow-up of 29 months; IQR 15-45) and 52 in non-initiators (39 months; 18-47): standardised 4-year risks were 0·76% (95% CI 0·51 to 1·04) in INSTI initiators and 0·75% (0·54 to 0·98) in non-INSTI initiators; risk ratio 1·01 (0·57 to 1·57); risk difference 0·0089% (-0·43 to 0·36). The analysis in ART-experienced individuals included 7875 INSTI initiators and 373 965 non-initiators. There were 56 events in INSTI initiators (median follow-up 18 months; IQR 9-29) and 3103 events (808 unique) in non-INSTI initiators (26 months; 15-37) in non-initiators: standardised 4-year risks 1·41% (95% CI 0·88 to 2·03) in INSTI initiators and 1·48% (1·28 to 1·71) in non-initiators; risk ratio 0·95 (0·60 to 1·36); risk difference -0·068% (-0·60 to 0·52). INTERPRETATION: We estimated that INSTI use did not result in a clinically meaningful increase of cardiovascular events in ART-naive and ART-experienced individuals with HIV. FUNDING: National Institute of Allergy and Infectious Diseases and National Institute on Alcohol Abuse and Alcoholism.

Do Polygenic Risk Scores Add to Clinical Data in Predicting Pancreatic Cancer? A Scoping Review.

BACKGROUND: Polygenic risk scores (PRS) summarize an individual's germline genetic risk, but it is unclear whether PRS offer independent information for pancreatic cancer risk prediction beyond routine clinical data. METHODS: We searched 8 databases from database inception to March 10, 2023 to identify studies evaluating the independent performance of pancreatic cancer-specific PRS for pancreatic cancer beyond clinical risk factors. RESULTS: Twenty-one studies examined associations between a pancreatic cancer-specific PRS and pancreatic cancer. Seven studies evaluated risk factors beyond age and sex. Three studies evaluated the change in discrimination associated with the addition of PRS to routine risk factors and reported improvements (AUCs: 0.715 to 0.745; AUC 0.791 to 0.830; AUC from 0.694 to 0.711). Limitations to clinical applicability included using source populations younger/healthier than those at risk for pancreatic cancer (n = 10), exclusively of European ancestry (n = 13), or controls without relevant exposures (n = 1). CONCLUSIONS: While most studies of pancreatic cancer-specific PRS did not evaluate the independent discrimination of PRS for pancreatic cancer beyond routine risk factors, three that did showed improvements in discrimination. IMPACT: For pancreatic cancer PRS to be clinically useful, they must demonstrate substantial improvements in discrimination beyond established risk factors, apply to diverse ancestral populations representative of those at risk for pancreatic cancer, and use appropriate controls.

Relative Burden of Cancer and Noncancer Mortality Among Long-Term Survivors of Breast, Prostate, and Colorectal Cancer in the US.

Importance: Improvements in cancer outcomes have led to a need to better understand long-term oncologic and nononcologic outcomes and quantify cancer-specific vs noncancer-specific mortality risks among long-term survivors. Objective: To assess absolute and relative cancer-specific vs noncancer-specific mortality rates among long-term survivors of cancer, as well as associated risk factors. Design, Setting, and Participants: This cohort study included 627 702 patients in the Surveillance, Epidemiology, and End Results cancer registry with breast, prostate, or colorectal cancer who received a diagnosis between January 1, 2003, and December 31, 2014, who received definitive treatment for localized disease and who were alive 5 years after their initial diagnosis (ie, long-term survivors of cancer). Statistical analysis was conducted from November 2022 to January 2023. Main Outcomes and Measures: Survival time ratios (TRs) were calculated using accelerated failure time models, and the primary outcome of interest examined was death from index cancer vs alternative (nonindex cancer) mortality across breast, prostate, colon, and rectal cancer cohorts. Secondary outcomes included subgroup mortality in cancer-specific risk groups, categorized based on prognostic factors, and proportion of deaths due to cancer-specific vs noncancer-specific causes. Independent variables included age, sex, race and ethnicity, income, residence, stage, grade, estrogen receptor status, progesterone receptor status, prostate-specific antigen level, and Gleason score. Follow-up ended in 2019. Results: The study included 627 702 patients (mean [SD] age, 61.1 [12.3] years; 434 848 women [69.3%]): 364 230 with breast cancer, 118 839 with prostate cancer, and 144 633 with colorectal cancer who survived 5 years or more from an initial diagnosis of early-stage cancer. Factors associated with shorter median cancer-specific survival included stage III disease for breast cancer (TR, 0.54; 95% CI, 0.53-0.55) and colorectal cancer (colon: TR, 0.60; 95% CI, 0.58-0.62; rectal: TR, 0.71; 95% CI, 0.69-0.74), as well as a Gleason score of 8 or higher for prostate cancer (TR, 0.61; 95% CI, 0.58-0.63). For all cancer cohorts, patients at low risk had at least a 3-fold higher noncancer-specific mortality compared with cancer-specific mortality at 10 years of diagnosis. Patients at high risk had a higher cumulative incidence of cancer-specific mortality than noncancer-specific mortality in all cancer cohorts except prostate. Conclusions and Relevance: This study is the first to date to examine competing oncologic and nononcologic risks focusing on long-term adult survivors of cancer. Knowledge of the relative risks facing long-term survivors may help provide pragmatic guidance to patients and clinicians regarding the importance of ongoing primary and oncologic-focused care.

Racial, Ethnic, and Rural Disparities in US Veteran COVID-19 Vaccine Rates.

Background: Race, ethnicity, and rurality-related disparities in coronavirus disease 2019 (COVID-19) vaccine uptake have been documented in the United States (US). Objective: We determined whether these disparities existed among patients at the Department of Veterans Affairs (VA), the largest healthcare system in the US. Design Settings Participants Measurements: Using VA Corporate Data Warehouse data, we included 5,871,438 patients (9.4% women) with at least one primary care visit in 2019 in a retrospective cohort study. Each patient was assigned a single race/ethnicity, which were mutually exclusive, self-reported categories. Rurality was based on 2019 home address at the zip code level. Our primary outcome was time-to-first COVID-19 vaccination between December 15, 2020-June 15, 2021. Additional covariates included age (in years), sex, geographic region (North Atlantic, Midwest, Southeast, Pacific, Continental), smoking status (current, former, never), Charlson Comorbidity Index (based on ≥1 inpatient or two outpatient ICD codes), service connection (any/none, using standardized VA-cutoffs for disability compensation), and influenza vaccination in 2019-2020 (yes/no). Results: Compared with unvaccinated patients, those vaccinated (n=3,238,532; 55.2%) were older (mean age in years vaccinated=66.3, (standard deviation=14.4) vs. unvaccinated=57.7, (18.0), p<.0001)). They were more likely to identify as Black (18.2% vs. 16.1%, p<.0001), Hispanic (7.0% vs. 6.6% p<.0001), or Asian American/Pacific Islander (AA/PI) (2.0% vs. 1.7%, P<.0001). In addition, they were more likely to reside in urban settings (68.0% vs. 62.8, p<.0001). Relative to non-Hispanic White urban Veterans, the reference group for race/ethnicity-urban/rural hazard ratios reported, all urban race/ethnicity groups were associated with increased likelihood for vaccination except American Indian/Alaskan Native (AI/AN) groups. Urban Black groups were 12% more likely (Hazard Ratio (HR)=1.12 [CI 1.12-1.13]) and rural Black groups were 6% more likely to receive a first vaccination (HR=1.06 [1.05-1.06]) relative to white urban groups. Urban Hispanic, AA/PI and Mixed groups were more likely to receive vaccination while rural members of these groups were less likely (Hispanic: Urban HR=1.17 [1.16-1.18], Rural HR=0.98 [0.97-0.99]; AA/PI: Urban HR=1.22 [1.21-1.23], Rural HR=0.86 [0.84-0.88]). Rural White Veterans were 21% less likely to receive an initial vaccine compared with urban White Veterans (HR=0.79 [0.78-0.79]). AI/AN groups were less likely to receive vaccination regardless of rurality: Urban HR=0.93 [0.91-0.95]; AI/AN-Rural HR=0.76 [0.74-0.78]. Conclusions: Urban Black, Hispanic, and AA/PI Veterans were more likely than their urban White counterparts to receive a first vaccination; all rural race/ethnicity groups except Black patients had lower likelihood for vaccination compared with urban White patients. A better understanding of disparities and rural outreach will inform equitable vaccine distribution.

Identifying genetic loci and phenomic associations of substance use traits: A multi-trait analysis of GWAS (MTAG) study.

BACKGROUND AND AIMS: Genome-wide association studies (GWAS) of opioid use disorder (OUD) and cannabis use disorder (CUD) have lagged behind those of alcohol use disorder (AUD) and smoking, where many more loci have been identified. We sought to identify novel loci for substance use traits (SUTs) in both African- (AFR) and European- (EUR) ancestry individuals to enhance our understanding of the traits' genetic architecture. DESIGN: We used multi-trait analysis of GWAS (MTAG) to analyze four SUTs in EUR subjects (OUD, CUD, AUD and smoking initiation [SMKinitiation]), and three SUTs in AFR subjects (OUD, AUD and smoking trajectory [SMKtrajectory]). We conducted gene-set and protein-protein interaction analyses and calculated polygenic risk scores (PRS) in two independent samples. SETTING: This study was conducted in the United States. PARTICIPANTS: A total of 5692 EUR and 4918 AFR individuals in the Yale-Penn sample and 29 054 EUR and 10 265 AFR individuals in the Penn Medicine BioBank sample. FINDINGS: MTAG identified genome-wide significant (GWS) single nucleotide polymorphisms (SNPs) for all four traits in EUR: 41 SNPs in 36 loci for OUD; 74 SNPs in 60 loci for CUD; 63 SNPs in 52 loci for AUD; and 183 SNPs in 144 loci for SMKinitiation. MTAG also identified GWS SNPs in AFR: 2 SNPs in 2 loci for OUD; 3 SNPs in 3 loci for AUD; and 1 SNP in 1 locus for SMKtrajectory. In the Yale-Penn sample, the MTAG-derived PRS consistently yielded more significant associations with both the corresponding substance use disorder diagnosis and multiple related phenotypes than the GWAS-derived PRS. CONCLUSIONS: Multi-trait analysis of genome-wide association studies boosted the number of loci found for substance use traits, identifying genes not previously linked to any substance, and increased the power of polygenic risk scores. Multi-trait analysis of genome-wide association studies can be used to identify novel associations for substance use, especially those for which the samples are smaller than those for historically legal substances.

Multi-trait genome-wide association analyses leveraging alcohol use disorder findings identify novel loci for smoking behaviors in the Million Veteran Program.

Smoking behaviors and alcohol use disorder (AUD), both moderately heritable traits, commonly co-occur in the general population. Single-trait genome-wide association studies (GWAS) have identified multiple loci for smoking and AUD. However, GWASs that have aimed to identify loci contributing to co-occurring smoking and AUD have used small samples and thus have not been highly informative. Applying multi-trait analysis of GWASs (MTAG), we conducted a joint GWAS of smoking and AUD with data from the Million Veteran Program (N = 318,694). By leveraging GWAS summary statistics for AUD, MTAG identified 21 genome-wide significant (GWS) loci associated with smoking initiation and 17 loci associated with smoking cessation compared to 16 and 8 loci, respectively, identified by single-trait GWAS. The novel loci for smoking behaviors identified by MTAG included those previously associated with psychiatric or substance use traits. Colocalization analysis identified 10 loci shared by AUD and smoking status traits, all of which achieved GWS in MTAG, including variants on SIX3, NCAM1, and near DRD2. Functional annotation of the MTAG variants highlighted biologically important regions on ZBTB20, DRD2, PPP6C, and GCKR that contribute to smoking behaviors. In contrast, MTAG of smoking behaviors and alcohol consumption (AC) did not enhance discovery compared with single-trait GWAS for smoking behaviors. We conclude that using MTAG to augment the power of GWAS enables the identification of novel genetic variants for commonly co-occuring phenotypes, providing new insights into their pleiotropic effects on smoking behavior and AUD.

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.

Assessing the contributions of modifiable risk factors to serious falls and fragility fractures among older persons living with HIV.

BACKGROUND: Although 50 years represents middle age among uninfected individuals, studies have shown that persons living with HIV (PWH) begin to demonstrate elevated risk for serious falls and fragility fractures in the sixth decade; the proportions of these outcomes attributable to modifiable factors are unknown. METHODS: We analyzed 21,041 older PWH on antiretroviral therapy (ART) from the Veterans Aging Cohort Study from 01/01/2010 through 09/30/2015. Serious falls were identified by Ecodes and a machine-learning algorithm applied to radiology reports. Fragility fractures (hip, vertebral, and upper arm) were identified using ICD9 codes. Predictors for both models included a serious fall within the past 12 months, body mass index, physiologic frailty (VACS Index 2.0), illicit substance and alcohol use disorders, and measures of multimorbidity and polypharmacy. We separately fit multivariable logistic models to each outcome using generalized estimating equations. From these models, the longitudinal extensions of average attributable fraction (LE-AAF) for modifiable risk factors were estimated. RESULTS: Key risk factors for both outcomes included physiologic frailty (VACS Index 2.0) (serious falls [15%; 95% CI 14%-15%]; fractures [13%; 95% CI 12%-14%]), a serious fall in the past year (serious falls [7%; 95% CI 7%-7%]; fractures [5%; 95% CI 4%-5%]), polypharmacy (serious falls [5%; 95% CI 4%-5%]; fractures [5%; 95% CI 4%-5%]), an opioid prescription in the past month (serious falls [7%; 95% CI 6%-7%]; fractures [9%; 95% CI 8%-9%]), and diagnosis of alcohol use disorder (serious falls [4%; 95% CI 4%-5%]; fractures [8%; 95% CI 7%-8%]). CONCLUSIONS: This study confirms the contributions of risk factors important in the general population to both serious falls and fragility fractures among older PWH. Successful prevention programs for these outcomes should build on existing prevention efforts while including risk factors specific to PWH.

Cell-type specific EWAS identifies genes involved in HIV pathogenesis and oncogenesis among people with HIV infection.

Epigenome-wide association studies (EWAS) of heterogenous blood cells have identified CpG sites associated with chronic HIV infection, which offer limited knowledge of cell-type specific methylation patterns associated with HIV infection. Applying a computational deconvolution method validated by capture bisulfite DNA methylation sequencing, we conducted a cell type-based EWAS and identified differentially methylated CpG sites specific for chronic HIV infection among five immune cell types in blood: CD4+ T-cells, CD8+ T-cells, B cells, Natural Killer (NK) cells, and monocytes in two independent cohorts (N total =1,134). Differentially methylated CpG sites for HIV-infection were highly concordant between the two cohorts. Cell-type level meta-EWAS revealed distinct patterns of HIV-associated differential CpG methylation, where 67% of CpG sites were unique to individual cell types (false discovery rate, FDR <0.05). CD4+ T-cells had the largest number of HIV-associated CpG sites (N=1,472) compared to any other cell type. Genes harboring statistically significant CpG sites are involved in immunity and HIV pathogenesis (e.g. CX3CR1 in CD4+ T-cells, CCR7 in B cells, IL12R in NK cells, LCK in monocytes). More importantly, HIV-associated CpG sites were overrepresented for hallmark genes involved in cancer pathology ( FDR <0.05) (e.g. BCL family, PRDM16, PDCD1LGD, ESR1, DNMT3A, NOTCH2 ). HIV-associated CpG sites were enriched among genes involved in HIV pathogenesis and oncogenesis such as Kras-signaling, interferon-α and -γ, TNF-α, inflammatory, and apoptotic pathways. Our findings are novel, uncovering cell-type specific modifications in the host epigenome for people with HIV that contribute to the growing body of evidence regarding pathogen-induced epigenetic oncogenicity, specifically on HIV and its comorbidity with cancers.

Enhanced Identification of Hispanic Ethnicity Using Clinical Data: A Study in the Largest Integrated United States Health Care System.

BACKGROUND: Collection of accurate Hispanic ethnicity data is critical to evaluate disparities in health and health care. However, this information is often inconsistently recorded in electronic health record (EHR) data. OBJECTIVE: To enhance capture of Hispanic ethnicity in the Veterans Affairs EHR and compare relative disparities in health and health care. METHODS: We first developed an algorithm based on surname and country of birth. We then determined sensitivity and specificity using self-reported ethnicity from the 2012 Veterans Aging Cohort Study survey as the reference standard and compared this to the research triangle institute race variable from the Medicare administrative data. Finally, we compared demographic characteristics and age-adjusted and sex-adjusted prevalence of conditions in Hispanic patients among different identification methods in the Veterans Affairs EHR 2018-2019. RESULTS: Our algorithm yielded higher sensitivity than either EHR-recorded ethnicity or the research triangle institute race variable. In 2018-2019, Hispanic patients identified by the algorithm were more likely to be older, had a race other than White, and foreign born. The prevalence of conditions was similar between EHR and algorithm ethnicity. Hispanic patients had higher prevalence of diabetes, gastric cancer, chronic liver disease, hepatocellular carcinoma, and human immunodeficiency virus than non-Hispanic White patients. Our approach evidenced significant differences in burden of disease among Hispanic subgroups by nativity status and country of birth. CONCLUSIONS: We developed and validated an algorithm to supplement Hispanic ethnicity information using clinical data in the largest integrated US health care system. Our approach enabled clearer understanding of demographic characteristics and burden of disease in the Hispanic Veteran population.

Delirium among people aging with and without HIV: Role of alcohol and Neurocognitively active medications.

BACKGROUND: People aging with and without HIV (PWH and PWoH) want to avoid neurocognitive dysfunction, especially delirium. Continued use of alcohol in conjunction with neurocognitively active medications (NCAMs) may be a largely underappreciated cause, especially for PWH who experience polypharmacy a decade earlier than PWoH. We compare absolute and relative risk of delirium among PWH and PWoH by age, level of alcohol use, and exposure to NCAMs. METHODS: Using the VACS cohort, we compare absolute and relative risk of inpatient delirium among PWH and PWoH by age, level of alcohol use, and exposure to NCAMs between 2007 and 2019. We matched each case based on age, race/ethnicity, sex, HIV, baseline year, and observation time with up to 5 controls. The case/control date was defined as date of admission for cases and the date corresponding to the same length of time on study for controls. Level of alcohol use was defined using Alcohol Use Disorder Identification Test-Consumption (AUDIT-C). Medication exposure was measured from 45 to 3 days prior to index date; medications were classified as anticholinergic NCAM, non-anticholinergic NCAM, or non NCAM and counts generated. We used logistic regression to determine odds ratios (ORs) for delirium associated with medication counts stratified by HIV status and adjusted for demographics, severity of illness, and related diagnoses. RESULTS: PWH experienced a higher incidence of delirium (5.6, [95% CI 5.3-5.9/1000 PY]) than PWoH (5.0, [95% CI 4.8-5.1/1000 PY]). In multivariable analysis, anticholinergic and non-anticholinergic NCAM counts and level of alcohol use demonstrated strong independent dose-response associations with delirium. CONCLUSIONS: Decreasing alcohol use and limiting the use of neurocognitively active medications may help decrease excess rates of delirium, especially among PWH.