Recalibrating the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO) Environmental Risk Score for Use in US Veterans.

BACKGROUND: Risk for colorectal cancer (CRC) may accumulate through multiple environmental factors. Understanding their effects, along with genetics, age and family history, could allow improvements in clinical decisions for screening protocols. We aimed to extend previous work by recalibrating an environmental risk score (e-Score) for CRC among a sample of US Veteran participants of the Million Veteran Program (MVP). METHODS: Demographic, lifestyle, and CRC data from 2011-2022 were abstracted from survey responses and health records of 227,504 male MVP participants. Weighting for each environmental factor's effect size was recalculated using VA training data to create a recalibrated e-Score. This recalibrated score was compared with the original weighted e-Score in a validation sample of 113,752 (n cases=590). Nested multiple logistic regression models tested associations between quintiles for recalibrated and original e-Scores. Likelihood Ratio Tests were used to compare model performance. RESULTS: Age (p<0.0001), education (p<0.0001), diabetes (p<0.0001), physical activity (p<0.0001), smoking (p<0.0001), NSAID use (p<0.0001), calcium (p=0.015), folate (p=0.020), and fruit consumption (p=0.019) were significantly different between CRC case and control groups. In the validation sample, the recalibrated e-Score model significantly improved the base model performance (p<0.001), but the original e-Score model did not (p=0.07). The recalibrated e-Score model quintile 5 was associated with significantly higher odds for CRC compared with quintile 1 (Q5 vs Q1: 1.79, 95% CI: 1.38-2.33). CONCLUSIONS: Multiple environmental factors, and the recalibrated e-Score quintiles were significantly associated with CRC cases. IMPACT: A recalibrated, Veteran-specific e-Score could be used to help personalize CRC screening and prevention strategies.

Protein Identification for Stroke Progression via Mendelian Randomization in Million Veteran Program and UK Biobank.

BACKGROUND: Individuals who have experienced a stroke, or transient ischemic attack, face a heightened risk of future cardiovascular events. Identification of genetic and molecular risk factors for subsequent cardiovascular outcomes may identify effective therapeutic targets to improve prognosis after an incident stroke. METHODS: We performed genome-wide association studies for subsequent major adverse cardiovascular events (MACE; ncases=51 929; ncontrols=39 980) and subsequent arterial ischemic stroke (AIS; ncases=45 120; ncontrols=46 789) after the first incident stroke within the Million Veteran Program and UK Biobank. We then used genetic variants associated with proteins (protein quantitative trait loci) to determine the effect of 1463 plasma protein abundances on subsequent MACE using Mendelian randomization. RESULTS: Two variants were significantly associated with subsequent cardiovascular events: rs76472767 near gene RNF220 (odds ratio, 0.75 [95% CI, 0.64-0.85]; P=3.69×10-8) with subsequent AIS and rs13294166 near gene LINC01492 (odds ratio, 1.52 [95% CI, 1.37-1.67]; P=3.77×10-8) with subsequent MACE. Using Mendelian randomization, we identified 2 proteins with an effect on subsequent MACE after a stroke: CCL27 ([C-C motif chemokine 27], effect odds ratio, 0.77 [95% CI, 0.66-0.88]; adjusted P=0.05) and TNFRSF14 ([tumor necrosis factor receptor superfamily member 14], effect odds ratio, 1.42 [95% CI, 1.24-1.60]; adjusted P=0.006). These proteins are not associated with incident AIS and are implicated to have a role in inflammation. CONCLUSIONS: We found evidence that 2 proteins with little effect on incident stroke appear to influence subsequent MACE after incident AIS. These associations suggest that inflammation is a contributing factor to subsequent MACE outcomes after incident AIS and highlights potential novel targets.

Estimating the Effect of Bariatric Surgery on Cardiovascular Events Using Observational Data?

BACKGROUND: Observational studies have reported strongly protective effects of bariatric surgery on cardiovascular disease, but with oversimplified definitions of the intervention, eligibility criteria, and follow-up, which deviate from those in a randomized trial. We describe an attempt to estimate the effect of bariatric surgery on cardiovascular disease without introducing these sources of bias, which may not be entirely possible with existing observational data. METHODS: We propose two target trials among persons with diabetes: (1) bariatric operation (vs. no operation) among individuals who have undergone preoperative preparation (lifestyle modifications and screening) and (2) preoperative preparation and a bariatric operation (vs. neither preoperative nor operative component). We emulated both target trials using observational data of US veterans. RESULTS: Comparing bariatric surgery with no surgery (target trial #1; 8,087 individuals), the 7-year cardiovascular risk was 18.0% (95% CI = 6.9, 32.7) in the surgery group and 18.9% (95% CI = 17.7, 20.1) in the no-surgery group (risk difference -0.9, 95% CI = -12.0, 14.0). Comparing preoperative components plus surgery vs. neither (target trial #2; 10,065 individuals), the 7-year cardiovascular risk was 17.4% (95% CI = 13.6, 22.0) in the surgery group and 18.8% (95% CI = 17.8, 19.9) in the no-surgery group (risk difference -1.4, 95% CI = -5.1, 3.2). Body mass index and hemoglobin A1c were reduced with bariatric interventions in both emulations. CONCLUSIONS: Within limitations of available observational data, our estimates do not provide evidence that bariatric surgery reduces cardiovascular disease and support equipoise for a randomized trial of bariatric surgery for cardiovascular disease prevention.

Long-Term Body Mass Index Variability and Adverse Cardiovascular Outcomes.

Importance: Body mass index (BMI; calculated as weight in kilograms divided by height in meters squared) is a commonly used estimate of obesity, which is a complex trait affected by genetic and lifestyle factors. Marked weight gain and loss could be associated with adverse biological processes. Objective: To evaluate the association between BMI variability and incident cardiovascular disease (CVD) events in 2 distinct cohorts. Design, Setting, and Participants: This cohort study used data from the Million Veteran Program (MVP) between 2011 and 2018 and participants in the UK Biobank (UKB) enrolled between 2006 and 2010. Participants were followed up for a median of 3.8 (5th-95th percentile, 3.5) years. Participants with baseline CVD or cancer were excluded. Data were analyzed from September 2022 and September 2023. Exposure: BMI variability was calculated by the retrospective SD and coefficient of variation (CV) using multiple clinical BMI measurements up to the baseline. Main Outcomes and Measures: The main outcome was incident composite CVD events (incident nonfatal myocardial infarction, acute ischemic stroke, and cardiovascular death), assessed using Cox proportional hazards modeling after adjustment for CVD risk factors, including age, sex, mean BMI, systolic blood pressure, total cholesterol, high-density lipoprotein cholesterol, smoking status, diabetes status, and statin use. Secondary analysis assessed whether associations were dependent on the polygenic score of BMI. Results: Among 92 363 US veterans in the MVP cohort (81 675 [88%] male; mean [SD] age, 56.7 [14.1] years), there were 9695 Hispanic participants, 22 488 non-Hispanic Black participants, and 60 180 non-Hispanic White participants. A total of 4811 composite CVD events were observed from 2011 to 2018. The CV of BMI was associated with 16% higher risk for composite CVD across all groups (hazard ratio [HR], 1.16; 95% CI, 1.13-1.19). These associations were unchanged among subgroups and after adjustment for the polygenic score of BMI. The UKB cohort included 65 047 individuals (mean [SD] age, 57.30 (7.77) years; 38 065 [59%] female) and had 6934 composite CVD events. Each 1-SD increase in BMI variability in the UKB cohort was associated with 8% increased risk of cardiovascular death (HR, 1.08; 95% CI, 1.04-1.11). Conclusions and Relevance: This cohort study found that among US veterans, higher BMI variability was a significant risk marker associated with adverse cardiovascular events independent of mean BMI across major racial and ethnic groups. Results were consistent in the UKB for the cardiovascular death end point. Further studies should investigate the phenotype of high BMI variability.

Red Meat Intake and the Risk of Cardiovascular Diseases: A Prospective Cohort Study in the Million Veteran Program.

BACKGROUND: Red meat consumption was associated with an increased risk of cardiovascular disease (CVD) in prospective cohort studies and a profile of biomarkers favoring high CVD risk in short-term controlled trials. However, several recent systematic reviews and meta-analyses concluded with no or weak evidence for limiting red meat intake. OBJECTIVES: To prospectively examine the associations between red meat intake and incident CVD in an ongoing cohort study with diverse socioeconomic and racial or ethnic backgrounds. METHODS: Our study included 148,506 participants [17,804 female (12.0%)] who were free of cancer, diabetes, and CVD at baseline from the Million Veteran Program. A food frequency questionnaire measured red meat intakes at baseline. Nonfatal myocardial infarction and acute ischemic stroke were identified through a high-throughput phenotyping algorithm, and fatal CVD events were identified by searching the National Death Index. RESULTS: Comparing the extreme categories of intake, the multivariate-adjusted relative risks of CVD was 1.18 (95% CI: 1.01, 1.38; P-trend < 0.0001) for total red meat, 1.14 (95% CI: 0.96, 1.36; P-trend = 0.01) for unprocessed red meat, and 1.29 (95% CI: 1.04, 1.60; P-trend = 0.003) for processed red meat. We observed a more pronounced positive association between red meat intake and CVD in African American participants than in White participants (P-interaction = 0.01). Replacing 0.5 servings/d of red meat with 0.5 servings/d of nuts, whole grains, and skimmed milk was associated with 14% (RR: 0.86; 95% CI: 0.83, 0.90), 7% (RR: 0.93; 95% CI: 0.89, 0.96), and 4% (RR: 0.96; 95% CI: 0.94, 0.99) lower risks of CVD, respectively. CONCLUSIONS: Red meat consumption is associated with an increased risk of CVD. Our findings support lowering red meat intake and replacing red meat with plant-based protein sources or low-fat dairy foods as a key dietary recommendation for the prevention of CVD.

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.

Impact of 8 lifestyle factors on mortality and life expectancy among United States veterans: The Million Veteran Program.

BACKGROUND: Lifestyle medicine has been proposed as a way to address the root causes of chronic disease and their associated health care costs. OBJECTIVE: This study aimed to estimate mortality risk and longevity associated with individual lifestyle factors and comprehensive lifestyle therapy. METHODS: Age- and sex-specific mortality rates were calculated on the basis of 719,147 veterans aged 40-99 y enrolled in the Veteran Affairs Million Veteran Program (2011-2019). Hazard ratios and estimated increase in life expectancy were examined among a subgroup of 276,132 veterans with complete data on 8 lifestyle factors at baseline. The 8 lifestyle factors included never smoking, physical activity, no excessive alcohol consumption, restorative sleep, nutrition, stress management, social connections, and no opioid use disorder. RESULTS: On the basis of 1.12 million person-years of follow-up, 34,247 deaths were recorded. Among veterans who adopted 1, 2, 3, 4, 5, 6, 7, and 8 lifestyle factors, the adjusted hazard ratios for mortality were 0.74 (0.60-0.90), 0.60 (95% CI: 0.49, 0.73), 0.50 (95% CI: 0.41, 0.61), 0.43 (95% CI: 0.35, 0.52), 0.35 (95% CI: 0.29, 0.43), 0.27 (95% CI: 0.22, 0.33), 0.21 (95% CI: 0.17, 0.26), and 0.13 (95% CI: 0.10, 0.16), respectively, as compared with veterans with no adopted lifestyle factors. The estimated life expectancy at age 40 y was 23.0, 26.5, 28.8, 30.8, 32.7, 35.1, 38.3, 41.3, and 47.0 y among males and 27.0, 28.8, 33.1, 38.0, 39.2, 41.4, 43.8, 46.3, and 47.5 y for females who adopted 0, 1, 2, 3, 4, 5, 6, 7, and 8 lifestyle factors, respectively. The difference in life expectancy at age 40 y was 24.0 y for male veterans and 20.5 y for female veterans when comparing adoption of 8-9 lifestyle factors. CONCLUSIONS: A combination of 8 lifestyle factors is associated with a significantly lower risk of premature mortality and an estimated prolonged life expectancy.

Moderate alcohol consumption on the risk of stroke in the Million Veteran Program.

BACKGROUND: There is inconsistent evidence on the association of moderate alcohol consumption and stroke risk in the general population and is not well studied among U.S. Veterans. Furthermore, it is unclear whether primarily drinking beer, wine, or liquor is associated with a difference in stroke risk. METHODS: The study included 185,323 Million Veteran Program participants who self-reported alcohol consumption on the Lifestyle Survey. Moderate consumption was defined as 1-2 drinks/day and beverage preference of beer, wine or liquor was defined if ≥ 50% of total drinks consumed were from a single type of beverage. Strokes were defined using ICD-9 and ICD-10 codes from the participants' electronic health record. RESULTS: The mean (sd) age of the sample was 64 (13) years and 11% were women. We observed 4,339 (94% ischemic; 6% hemorrhagic) strokes over a median follow-up of 5.2 years. In Cox models adjusted for age, sex, race, education, income, body mass index, smoking, exercise, diet, cholesterol, prevalent diabetes, prevalent hypertension, lipid-lowering medication, antihypertensive medication, and diabetes medication, moderate alcohol consumption (1-2 drinks/day) was associated with a 22% lower risk of total stroke compared with never drinking [Hazards ratio (HR) 95% confidence interval (CI): 0.78 (0.67, 0.92)]. When stratifying by stroke type, we observed a similar protective association with moderate consumption and ischemic stroke [HR (95% CI): 0.76 (0.65, 0.90)], but a non-statistically significant higher risk of hemorrhagic stroke [HR (95% CI): 1.29 (0.64, 2.61)]. We did not observe a difference in ischemic or hemorrhagic stroke risk among those who preferred beer, liquor or wine vs. no beverage preference. When stratifying by prior number of hospital visits (≤ 15, 16-33, 34-64, ≥ 65) as a proxy for health status, we observed attenuation of the protective association with greater number of visits [HR (95% CI): 0.87 (0.63, 1.19) for ≥ 65 visits vs. 0.80 (0.59, 1.08) for ≤ 15 visits]. CONCLUSIONS: We observed a lower risk of ischemic stroke, but not hemorrhagic stroke with moderate alcohol consumption and did not observe substantial differences in risk by beverage preference among a sample of U.S. Veterans. Healthy user bias of moderate alcohol consumption may be driving some of the observed protective association.

Characterizing prostate cancer risk through multi-ancestry genome-wide discovery of 187 novel risk variants.

The transferability and clinical value of genetic risk scores (GRSs) across populations remain limited due to an imbalance in genetic studies across ancestrally diverse populations. Here we conducted a multi-ancestry genome-wide association study of 156,319 prostate cancer cases and 788,443 controls of European, African, Asian and Hispanic men, reflecting a 57% increase in the number of non-European cases over previous prostate cancer genome-wide association studies. We identified 187 novel risk variants for prostate cancer, increasing the total number of risk variants to 451. An externally replicated multi-ancestry GRS was associated with risk that ranged from 1.8 (per standard deviation) in African ancestry men to 2.2 in European ancestry men. The GRS was associated with a greater risk of aggressive versus non-aggressive disease in men of African ancestry (P = 0.03). Our study presents novel prostate cancer susceptibility loci and a GRS with effective risk stratification across ancestry groups.

Semi-supervised calibration of noisy event risk (SCANER) with electronic health records.

OBJECTIVE: Electronic health records (EHR), containing detailed longitudinal clinical information on a large number of patients and covering broad patient populations, open opportunities for comprehensive predictive modeling of disease progression and treatment response. However, since EHRs were originally constructed for administrative purposes not for research, in the EHR-linked studies, it is often not feasible to capture reliable information for analytical variables, especially in the survival setting, when both accurate event status and event times are needed for model building. For example, progression-free survival (PFS), a commonly used survival outcome for cancer patients, often involves complex information embedded in free-text clinical notes and cannot be extracted reliably. Proxies of PFS time such as time to the first mention of progression in the notes are at best good approximations to the true event time. This leads to difficulty in efficiently estimating event rates for an EHR patient cohort. Estimating survival rates based on error-prone outcome definitions can lead to biased results and hamper the power in the downstream analysis. On the other hand, extracting accurate event time information via manual annotation is time and resource intensive. The objective of this study is to develop a calibrated survival rate estimator using noisy outcomes from EHR data. MATERIALS AND METHODS: In this paper, we propose a two-stage semi-supervised calibration of noisy event rate (SCANER) estimator that can effectively overcome censoring induced dependency and attains more robust performance (i.e., not sensitive to misspecification of the imputation model) by fully utilizing both a small-labeled set of gold-standard survival outcomes annotated via manual chart review and a set of proxy features automatically captured via EHR in the unlabeled set. We validate the SCANER estimator by estimating the PFS rates for a virtual cohort of lung cancer patients from one large tertiary care center and the ICU-free survival rates for COVID patients from two large tertiary care centers. RESULTS: In terms of survival rate estimates, the SCANER had very similar point estimates compared to the complete-case Kaplan Meier estimator. On the other hand, other benchmark methods for comparison, which fail to account for the induced dependency between event time and the censoring time conditioning on surrogate outcomes, produced biased results across all three case studies. In terms of standard errors, the SCANER estimator was more efficient than the KM estimator, with up to 50% efficiency gain. CONCLUSION: The SCANER estimator achieves more efficient, robust, and accurate survival rate estimates compared to existing approaches. This promising new approach can also improve the resolution (i.e., granularity of event time) by using labels conditioning on multiple surrogates, particularly among less common or poorly coded conditions.

Cardiovascular Disease Risk Assessment Using Traditional Risk Factors and Polygenic Risk Scores in the Million Veteran Program.

Importance: Primary prevention of atherosclerotic cardiovascular disease (ASCVD) relies on risk stratification. Genome-wide polygenic risk scores (PRSs) are proposed to improve ASCVD risk estimation. Objective: To determine whether genome-wide PRSs for coronary artery disease (CAD) and acute ischemic stroke improve ASCVD risk estimation with traditional clinical risk factors in an ancestrally diverse midlife population. Design, Setting, and Participants: This was a prognostic analysis of incident events in a retrospectively defined longitudinal cohort conducted from January 1, 2011, to December 31, 2018. Included in the study were adults free of ASCVD and statin naive at baseline from the Million Veteran Program (MVP), a mega biobank with genetic, survey, and electronic health record data from a large US health care system. Data were analyzed from March 15, 2021, to January 5, 2023. Exposures: PRSs for CAD and ischemic stroke derived from cohorts of largely European descent and risk factors, including age, sex, systolic blood pressure, total cholesterol, high-density lipoprotein (HDL) cholesterol, smoking, and diabetes status. Main Outcomes and Measures: Incident nonfatal myocardial infarction (MI), ischemic stroke, ASCVD death, and composite ASCVD events. Results: A total of 79 151 participants (mean [SD] age, 57.8 [13.7] years; 68 503 male [86.5%]) were included in the study. The cohort included participants from the following harmonized genetic ancestry and race and ethnicity categories: 18 505 non-Hispanic Black (23.4%), 6785 Hispanic (8.6%), and 53 861 non-Hispanic White (68.0%) with a median (5th-95th percentile) follow-up of 4.3 (0.7-6.9) years. From 2011 to 2018, 3186 MIs (4.0%), 1933 ischemic strokes (2.4%), 867 ASCVD deaths (1.1%), and 5485 composite ASCVD events (6.9%) were observed. CAD PRS was associated with incident MI in non-Hispanic Black (hazard ratio [HR], 1.10; 95% CI, 1.02-1.19), Hispanic (HR, 1.26; 95% CI, 1.09-1.46), and non-Hispanic White (HR, 1.23; 95% CI, 1.18-1.29) participants. Stroke PRS was associated with incident stroke in non-Hispanic White participants (HR, 1.15; 95% CI, 1.08-1.21). A combined CAD plus stroke PRS was associated with ASCVD deaths among non-Hispanic Black (HR, 1.19; 95% CI, 1.03-1.17) and non-Hispanic (HR, 1.11; 95% CI, 1.03-1.21) participants. The combined PRS was also associated with composite ASCVD across all ancestry groups but greater among non-Hispanic White (HR, 1.20; 95% CI, 1.16-1.24) than non-Hispanic Black (HR, 1.11; 95% CI, 1.05-1.17) and Hispanic (HR, 1.12; 95% CI, 1.00-1.25) participants. Net reclassification improvement from adding PRS to a traditional risk model was modest for the intermediate risk group for composite CVD among men (5-year risk >3.75%, 0.38%; 95% CI, 0.07%-0.68%), among women, (6.79%; 95% CI, 3.01%-10.58%), for age older than 55 years (0.25%; 95% CI, 0.03%-0.47%), and for ages 40 to 55 years (1.61%; 95% CI, -0.07% to 3.30%). Conclusions and Relevance: Study results suggest that PRSs derived predominantly in European samples were statistically significantly associated with ASCVD in the multiancestry midlife and older-age MVP cohort. Overall, modest improvement in discrimination metrics were observed with addition of PRSs to traditional risk factors with greater magnitude in women and younger age groups.

Evidence of Novel Susceptibility Variants for Prostate Cancer and a Multiancestry Polygenic Risk Score Associated with Aggressive Disease in Men of African Ancestry.

BACKGROUND: Genetic factors play an important role in prostate cancer (PCa) susceptibility. OBJECTIVE: To discover common genetic variants contributing to the risk of PCa in men of African ancestry. DESIGN, SETTING, AND PARTICIPANTS: We conducted a meta-analysis of ten genome-wide association studies consisting of 19378 cases and 61620 controls of African ancestry. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: Common genotyped and imputed variants were tested for their association with PCa risk. Novel susceptibility loci were identified and incorporated into a multiancestry polygenic risk score (PRS). The PRS was evaluated for associations with PCa risk and disease aggressiveness. RESULTS AND LIMITATIONS: Nine novel susceptibility loci for PCa were identified, of which seven were only found or substantially more common in men of African ancestry, including an African-specific stop-gain variant in the prostate-specific gene anoctamin 7 (ANO7). A multiancestry PRS of 278 risk variants conferred strong associations with PCa risk in African ancestry studies (odds ratios [ORs] >3 and >5 for men in the top PRS decile and percentile, respectively). More importantly, compared with men in the 40-60% PRS category, men in the top PRS decile had a significantly higher risk of aggressive PCa (OR = 1.23, 95% confidence interval = 1.10-1.38, p = 4.4 × 10-4). CONCLUSIONS: This study demonstrates the importance of large-scale genetic studies in men of African ancestry for a better understanding of PCa susceptibility in this high-risk population and suggests a potential clinical utility of PRS in differentiating between the risks of developing aggressive and nonaggressive disease in men of African ancestry. PATIENT SUMMARY: In this large genetic study in men of African ancestry, we discovered nine novel prostate cancer (PCa) risk variants. We also showed that a multiancestry polygenic risk score was effective in stratifying PCa risk, and was able to differentiate risk of aggressive and nonaggressive disease.

Plant-based diets and the incidence of cardiovascular disease: the Million Veteran Program.

Background: A healthful plant-based diet was associated with lower risks of coronary heart disease and type 2 diabetes, and a favourable profile of adiposity-associated biomarkers, while an unhealthful plant-based diet was associated with elevated risk of cardiometabolic disease in health professional populations. However, little is known about the associations between plant-based dietary patterns and risk of cardiovascular disease (CVD) in US veterans. Methods: The study population consisted of 148 506 participants who were free of diabetes, CVD and cancer at baseline in the Veterans Affairs (VA) Million Veteran Program. Diet was assessed using a Food Frequency Questionnaire at baseline. We calculated an overall Plant-Based Diet Index (PDI), a healthful PDI (hPDI) and an unhealthful PDI (uPDI). The CVD endpoints included non-fatal myocardial infarction (MI) and acute ischaemic stroke (AIS) identified through high-throughput phenotyping algorithms approach and fatal CVD events identified by searching the National Death Index. Results: With up to 8 years of follow-up, we documented 5025 CVD cases. After adjustment for confounding factors, a higher PDI was significantly associated with a lower risk of CVD (HR comparing extreme quintiles=0.75, 95% CI 0.68 to 0.82, P trend<0.0001). We observed an inverse association between hPDI and the risk of CVD (HR comparing extreme quintiles=0.71, 95% CI 0.64 to 0.78, P trend<0.001), whereas uPDI was positively associated with the risk of CVD (HR comparing extreme quintiles=1.12, 95% CI 1.02 to 1.24, P trend<0.001). We found similar associations of hPDI with subtypes of CVD; a 10-unit increment in hPDI was associated with HRs (95% CI) of 0.81 (0.75 to 0.87) for fatal CVD, 0.86 (0.79 to 0.94) for non-fatal MI and 0.86 (0.78 to 0.95) for non-fatal AIS. Conclusions: Plant-based dietary pattern enriched with healthier plant foods was associated with a substantially lower CVD risk in US veterans.

Mild-to-Moderate Kidney Dysfunction and Cardiovascular Disease: Observational and Mendelian Randomization Analyses.

BACKGROUND: End-stage renal disease is associated with a high risk of cardiovascular events. It is unknown, however, whether mild-to-moderate kidney dysfunction is causally related to coronary heart disease (CHD) and stroke. METHODS: Observational analyses were conducted using individual-level data from 4 population data sources (Emerging Risk Factors Collaboration, EPIC-CVD [European Prospective Investigation into Cancer and Nutrition-Cardiovascular Disease Study], Million Veteran Program, and UK Biobank), comprising 648 135 participants with no history of cardiovascular disease or diabetes at baseline, yielding 42 858 and 15 693 incident CHD and stroke events, respectively, during 6.8 million person-years of follow-up. Using a genetic risk score of 218 variants for estimated glomerular filtration rate (eGFR), we conducted Mendelian randomization analyses involving 413 718 participants (25 917 CHD and 8622 strokes) in EPIC-CVD, Million Veteran Program, and UK Biobank. RESULTS: There were U-shaped observational associations of creatinine-based eGFR with CHD and stroke, with higher risk in participants with eGFR values <60 or >105 mL·min-1·1.73 m-2, compared with those with eGFR between 60 and 105 mL·min-1·1.73 m-2. Mendelian randomization analyses for CHD showed an association among participants with eGFR <60 mL·min-1·1.73 m-2, with a 14% (95% CI, 3%-27%) higher CHD risk per 5 mL·min-1·1.73 m-2 lower genetically predicted eGFR, but not for those with eGFR >105 mL·min-1·1.73 m-2. Results were not materially different after adjustment for factors associated with the eGFR genetic risk score, such as lipoprotein(a), triglycerides, hemoglobin A1c, and blood pressure. Mendelian randomization results for stroke were nonsignificant but broadly similar to those for CHD. CONCLUSIONS: In people without manifest cardiovascular disease or diabetes, mild-to-moderate kidney dysfunction is causally related to risk of CHD, highlighting the potential value of preventive approaches that preserve and modulate kidney function.

Egg consumption and risk of acute stroke in the Million Veteran Program.

BACKGROUND: While observational data suggest a higher risk of coronary artery disease with frequent egg consumption, only limited and inconsistent data are available on the relation of egg consumption with stroke. OBJECTIVE: The primary objective was to assess whether egg consumption is associated with a higher risk of ischemic and hemorrhagic stroke among US veterans. METHODS: In a prospective cohort study of US veterans from the Million Veteran Program (MVP), egg intake was collected through a self-reported food frequency questionnaire at baseline. Incidence of stroke was ascertained via ICD9/ICD10 codes from the electronic health records. Multivariable Cox proportional hazard models were used to estimate the multivariable adjusted hazard ratios. RESULTS: A total of 233,792 veterans (91.6% men) were studied with a mean age of 65.6 ± 11.7 years. During a mean follow-up of 3.3 years, a total of 5740 cases of fatal and non-fatal ischemic stroke and 423 cases of fatal and non-fatal hemorrhagic stroke occurred. Median egg consumption was 2-4 eggs/week. Crude incidence rates for acute ischemic stroke were 6.5, 7.2, 7.1, 7.4, 8.0, 8.1, and 8.6 cases per 1000 person-years for egg consumption of <1/month, 1-3/month, 1/week, 2-4/week, 5-6/week, 1/day, and ≥2/day, respectively. Corresponding multivariable adjusted hazard ratios (95% CI) were 1.00 (ref), 1.10 (0.96-1.25), 1.09 (0.96-1.23), 1.10 (0.98-1.25), 1.16 (1.01-1.33), 1.20 (1.03-1.40), and 1.22 (1.03-1.45) controlling for age, sex, ethnicity, body mass index, diabetes, smoking, alcohol intake, modified DASH score, and education (p linear trend 0.0085). For hemorrhagic stroke, multivariable adjusted hazard ratios (95% CI) after controlling for age, sex, ethnicity, body mass index, modified DASH score, and level of education were 1.00 (ref), 1.28 (0.96-1.72), 1.22 (0.93-1.61), and 1.19 (0.88-1.61) for egg consumption of <1/week, 1/week, 2-4/week, 5+/week, respectively. CONCLUSION: Our data are consistent with a positive association of egg consumption with acute ischemic stroke but not hemorrhagic stroke among veterans.

Validation of a multi-ancestry polygenic risk score and age-specific risks of prostate cancer: A meta-analysis within diverse populations.

Background: We recently developed a multi-ancestry polygenic risk score (PRS) that effectively stratifies prostate cancer risk across populations. In this study, we validated the performance of the PRS in the multi-ancestry Million Veteran Program and additional independent studies. Methods: Within each ancestry population, the association of PRS with prostate cancer risk was evaluated separately in each case-control study and then combined in a fixed-effects inverse-variance-weighted meta-analysis. We further assessed the effect modification by age and estimated the age-specific absolute risk of prostate cancer for each ancestry population. Results: The PRS was evaluated in 31,925 cases and 490,507 controls, including men from European (22,049 cases, 414,249 controls), African (8794 cases, 55,657 controls), and Hispanic (1082 cases, 20,601 controls) populations. Comparing men in the top decile (90-100% of the PRS) to the average 40-60% PRS category, the prostate cancer odds ratio (OR) was 3.8-fold in European ancestry men (95% CI = 3.62-3.96), 2.8-fold in African ancestry men (95% CI = 2.59-3.03), and 3.2-fold in Hispanic men (95% CI = 2.64-3.92). The PRS did not discriminate risk of aggressive versus nonaggressive prostate cancer. However, the OR diminished with advancing age (European ancestry men in the top decile: ≤55 years, OR = 7.11; 55-60 years, OR = 4.26; >70 years, OR = 2.79). Men in the top PRS decile reached 5% absolute prostate cancer risk ~10 years younger than men in the 40-60% PRS category. Conclusions: Our findings validate the multi-ancestry PRS as an effective prostate cancer risk stratification tool across populations. A clinical study of PRS is warranted to determine whether the PRS could be used for risk-stratified screening and early detection. Funding: This work was supported by the National Cancer Institute at the National Institutes of Health (grant numbers U19 CA214253 to C.A.H., U01 CA257328 to C.A.H., U19 CA148537 to C.A.H., R01 CA165862 to C.A.H., K99 CA246063 to B.F.D, and T32CA229110 to F.C), the Prostate Cancer Foundation (grants 21YOUN11 to B.F.D. and 20CHAS03 to C.A.H.), the Achievement Rewards for College Scientists Foundation Los Angeles Founder Chapter to B.F.D, and the Million Veteran Program-MVP017. This research has been conducted using the UK Biobank Resource under application number 42195. This research is based on data from the Million Veteran Program, Office of Research and Development, and the Veterans Health Administration. This publication does not represent the views of the Department of Veteran Affairs or the United States Government.

Semi-supervised approach to event time annotation using longitudinal electronic health records.

Large clinical datasets derived from insurance claims and electronic health record (EHR) systems are valuable sources for precision medicine research. These datasets can be used to develop models for personalized prediction of risk or treatment response. Efficiently deriving prediction models using real world data, however, faces practical and methodological challenges. Precise information on important clinical outcomes such as time to cancer progression are not readily available in these databases. The true clinical event times typically cannot be approximated well based on simple extracts of billing or procedure codes. Whereas, annotating event times manually is time and resource prohibitive. In this paper, we propose a two-step semi-supervised multi-modal automated time annotation (MATA) method leveraging multi-dimensional longitudinal EHR encounter records. In step I, we employ a functional principal component analysis approach to estimate the underlying intensity functions based on observed point processes from the unlabeled patients. In step II, we fit a penalized proportional odds model to the event time outcomes with features derived in step I in the labeled data where the non-parametric baseline function is approximated using B-splines. Under regularity conditions, the resulting estimator of the feature effect vector is shown as root-n consistent. We demonstrate the superiority of our approach relative to existing approaches through simulations and a real data example on annotating lung cancer recurrence in an EHR cohort of lung cancer patients from Veteran Health Administration.

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.

Genome-wide and phenome-wide analysis of ideal cardiovascular health in the VA Million Veteran Program.

BACKGROUND: Genetic studies may help identify causal pathways; therefore, we sought to identify genetic determinants of ideal CVH and their association with CVD outcomes in the multi-population Veteran Administration Million Veteran Program. METHODS: An ideal health score (IHS) was calculated from 3 clinical factors (blood pressure, total cholesterol, and blood glucose levels) and 3 behavioral factors (smoking status, physical activity, and BMI), ascertained at baseline. Multi-population genome-wide association study (GWAS) was performed on IHS and binary ideal health using linear and logistic regression, respectively. Using the genome-wide significant SNPs from the IHS GWAS, we created a weighted IHS polygenic risk score (PRSIHS) which was used (i) to conduct a phenome-wide association study (PheWAS) of associations between PRSIHS and ICD-9 phenotypes and (ii) to further test for associations with mortality and selected CVD outcomes using logistic and Cox regression and, as an instrumental variable, in Mendelian Randomization. RESULTS: The discovery and replication cohorts consisted of 142,404 (119,129 European American (EUR); 16,495 African American (AFR)), and 45,766 (37,646 EUR; 5,366 AFR) participants, respectively. The mean age was 65.8 years (SD = 11.2) and 92.7% were male. Overall, 4.2% exhibited ideal CVH based on the clinical and behavioral factors. In the multi-population meta-analysis, variants at 17 loci were associated with IHS and each had known GWAS associations with multiple components of the IHS. PheWAS analysis in 456,026 participants showed that increased PRSIHS was associated with a lower odds ratio for many CVD outcomes and risk factors. Both IHS and PRSIHS measures of ideal CVH were associated with significantly less CVD outcomes and CVD mortality. CONCLUSION: A set of high interest genetic variants contribute to the presence of ideal CVH in a multi-ethnic cohort of US Veterans. Genetically influenced ideal CVH is associated with lower odds of CVD outcomes and mortality.

Degree of Adherence to Based Diet and Total and Cause-Specific Mortality: Prospective Cohort Study in the Million Veteran Program.

OBJECTIVE: To examine the associations between adherence to plant-based diets and mortality. DESIGN: prospective study. We calculated a plant-based diet index (PDI) by assigning positive scores to plant foods and reverse scores to animal foods. We also created a healthful PDI (hPDI) and an unhealthful PDI (uPDI) by further separate the healthy plant foods from less-healthy plant foods. SETTING: the VA Million Veteran Program. PARTICIPANTS: 315,919 men and women aged 19 to 104 years who completed a food frequency questionnaire at the baseline. RESULTS: We documented 31,136 deaths during the follow-up. A higher PDI was significantly associated with lower total mortality [hazard ratio (HR) comparing extreme deciles =0.75, 95% confidence interval (CI): 0.71 to 0.79, Ptrend <0.001]. We observed an inverse association between hPDI and total mortality (HR comparing extreme deciles =0.64, 95% CI: 0.61 to 0.68, Ptrend <0.001), whereas uPDI was positively associated with total mortality (HR comparing extreme deciles =1.41, 95% CI: 1.33 to 1.49, Ptrend <0.001). Similar significant associations of PDI, hPDI, and uPDI were also observed for CVD and cancer mortality. The associations between the plant-based diet indices and total mortality were consistent among African and European American participants, and participants free from CVD and cancer and those who were diagnosed with major chronic disease at baseline. CONCLUSIONS: A greater adherence to a plant-based diet was associated with substantially lower total mortality in this large population of veterans. These findings support recommending plant-rich dietary patterns for the prevention of major chronic diseases.