Associations among NMR-measured inflammatory and metabolic biomarkers and accelerated aging in cardiac catheterization patients

Research into aging has grown substantially with the creation of molecular biomarkers of biological age that can be used to determine age acceleration. Concurrently, nuclear magnetic resonance (NMR) assessment of biomarkers of inflammation and metabolism provides researchers with new ways to examine intermediate risk factors for chronic disease. We used data from a cardiac catheterization cohort to examine associations between biomarkers of cardiometabolic health and accelerated aging assessed using both gene expression (Transcriptomic Age) and DNA methylation (Hannum Age, GrimAge, Horvath Age, and Phenotypic Age). Linear regression models were used to associate accelerated aging with each outcome (cardiometabolic health biomarkers) while adjusting for chronological age, sex, race, and neighborhood socioeconomic status. Our study shows a robust association between GlycA and GrimAge (5.71, 95% CI = 4.36, 7.05, P = 7.94 × 10-16), Hannum Age (1.81, 95% CI = 0.65, 2.98, P = 2.30 × 10-3), and Phenotypic Age (2.88, 95% CI = 1.91, 3.87, P = 1.21 × 10-8). We also saw inverse associations between apolipoprotein A-1 and aging biomarkers. These associations provide insight into the relationship between aging and cardiometabolic health that may be informative for vulnerable populations.

Genome-Wide Genetic Analysis of Dropout in a Controlled Exercise Intervention in Sedentary Adults With Overweight or Obesity and Cardiometabolic Disease

BACKGROUND

Despite the benefits of exercise, many individuals are unable or unwilling to adopt an exercise intervention.

PURPOSE

The purpose of this analysis was to identify putative genetic variants associated with dropout from exercise training interventions among individuals in the STRRIDE trials.

METHODS

We used a genome-wide association study approach to identify genetic variants in 603 participants initiating a supervised exercise intervention. Exercise intervention dropout occurred when a subject withdrew from further participation in the study or was otherwise lost to follow-up.

RESULTS

Exercise intervention dropout was associated with a cluster of single-nucleotide polymorphisms with the top candidate being rs722069 (T/C, risk allele = C) (unadjusted p = 2.2 × 10-7, odds ratio = 2.23) contained within a linkage disequilibrium block on chromosome 16. In Genotype-Tissue Expression, rs722069 is an expression quantitative trait locus of the EARS2, COG7, and DCTN5 genes in skeletal muscle tissue. In subsets of the STRRIDE genetic cohort with available muscle gene expression (n = 37) and metabolic data (n = 82), at baseline the C allele was associated with lesser muscle expression of EARS2 (p < .002) and COG7 (p = .074) as well as lesser muscle concentrations of C2- and C3-acylcarnitines (p = .026).

CONCLUSIONS

Our observations imply that exercise intervention dropout is genetically moderated through alterations in gene expression and metabolic pathways in skeletal muscle. Individual genetic traits may allow the development of a biomarker-based approach for identifying individuals who may benefit from more intensive counseling and other interventions to optimize exercise intervention adoption.

CLINICAL TRIAL INFORMATION

STRRIDE I = NCT00200993; STRRIDE AT/RT = NCT00275145; STRRIDE-PD = NCT00962962.

Longitudinal Patterns of Multimorbidity in Gulf War Era Veterans With and Without Gulf War Illness

Objectives: To examine whether severe Gulf War illness (SGWI) case status was associated with longitudinal multimorbidity patterns. Methods: Participants were users of the Veteran Health Administration Health Care System drawn from the Gulf War Era Cohort and Biorepository (n = 840). Longitudinal measures of multimorbidity were constructed using (1) electronic health records (Charlson Comorbidity Index; Elixhauser; and Veterans Affairs Frailty Index) from 10/1/1999 to 6/30/2023 and (2) self-reported medical conditions (Deficit Accumulation Index) since the war until the survey date. Accelerated failure time models examined SGWI case status as a predictor of time until threshold level of multimorbidity was reached, adjusted for age and sociodemographic and military characteristics. Results: Models, adjusted for covariates, revealed that (1) relative to the SWGI- group, the SGWI+ group was associated with an accelerated time for reaching each threshold and (2) the relationship between SGWI and each threshold was not moderated by age. Discussion: Findings suggest that veterans with SGWI experienced accelerated aging.

Military exposures and Gulf War illness in veterans with and without posttraumatic stress disorder

Gulf War illness (GWI) is a chronic multisymptom disorder of unknown etiology that is believed to be caused by neurotoxicant exposure experienced during deployment to the Gulf War. Posttraumatic stress disorder (PTSD) covaries with GWI and is believed to play a role in GWI symptoms. The present study examined the association between self-reported military exposures and GWI, stratified by PTSD status, in veterans from the Gulf War Era Cohort and Biorepository who were deployed to the Persian Gulf during the war. Participants self-reported current GWI and PTSD symptoms as well as military exposures (e.g., pyridostigmine [PB] pills, pesticides/insecticides, combat, chemical attacks, and oil well fires) experienced during the Gulf War. Deployed veterans' (N = 921) GWI status was ascertained using the Centers for Disease Control and Prevention definition. Individuals who met the GWI criteria were stratified by PTSD status, yielding three groups: GWI-, GWI+/PTSD-, and GWI+/PTSD+. Multivariable logistic regression, adjusted for covariates, was used to examine associations between GWI/PTSD groups and military exposures. Apart from insect bait use, the GWI+/PTSD+ group had higher odds of reporting military exposures than the GWI+/PTSD- group, adjusted odds ratio (aOR) = 2.15, 95% CI [1.30, 3.56]-aOR = 6.91, 95% CI [3.39, 14.08]. Except for PB pills, the GWI+/PTSD- group had a higher likelihood of reporting military exposures than the GWI- group, aOR = 2.03, 95% CI [1.26, 3.26]-aOR = 4.01, 95% CI [1.57, 10.25]. These findings are consistent with roles for both PTSD and military exposures in the etiology of GWI.

The Million Veteran Program 1990-1991 Gulf War Era Survey: An Evaluation of Veteran Response, Characteristics, and Representativeness of the Gulf War Era Veteran Population

To address gaps in understanding the pathophysiology of Gulf War Illness (GWI), the VA Million Veteran Program (MVP) developed and implemented a survey to MVP enrollees who served in the U.S. military during the 1990-1991 Persian Gulf War (GW). Eligible Veterans were invited via mail to complete a survey assessing health conditions as well as GW-specific deployment characteristics and exposures. We evaluated the representativeness of this GW-era cohort relative to the broader population by comparing demographic, military, and health characteristics between respondents and non-respondents, as well as with all GW-era Veterans who have used Veterans Health Administration (VHA) services and the full population of U.S. GW-deployed Veterans. A total of 109,976 MVP GW-era Veterans were invited to participate and 45,270 (41%) returned a completed survey. Respondents were 84% male, 72% White, 8% Hispanic, with a mean age of 61.6 years (SD = 8.5). Respondents were more likely to be older, White, married, better educated, slightly healthier, and have higher socioeconomic status than non-respondents, but reported similar medical conditions and comparable health status. Although generally similar to all GW-era Veterans using VHA services and the full population of U.S. GW Veterans, respondents included higher proportions of women and military officers, and were slightly older. In conclusion, sample characteristics of the MVP GW-era cohort can be considered generally representative of the broader GW-era Veteran population. The sample represents the largest research cohort of GW-era Veterans established to date and provides a uniquely valuable resource for conducting in-depth studies to evaluate health conditions affecting 1990-1991 GW-era Veterans.

Association of Atherosclerotic Cardiovascular Disease, Hypertension, Diabetes, and Hyperlipidemia With Gulf War Illness Among Gulf War Veterans

BACKGROUND

Approximately 30% of the 700 000 Gulf War veterans report a chronic symptom-based illness of varying severity referred to as Gulf War illness (GWI). Toxic deployment-related exposures have been implicated in the cause of GWI, some of which contribute to metabolic dysregulation and lipid abnormalities. As this cohort ages, the relationship between GWI and atherosclerotic cardiovascular disease (ASCVD) is a growing concern. We evaluated associations between GWI and ASCVD, diabetes, hyperlipidemia, and hypertension in veterans of the Gulf War (1990-1991).

METHODS AND RESULTS

Analysis of survey data collected in 2014 to 2016 from a national sample of deployed Gulf War veterans (n=942) and Veterans Health Administration electronic health record data (n=669). Multivariable logistic regression models tested for associations of GWI with self-reported ASCVD, diabetes, hyperlipidemia, and hypertension, controlling for confounding factors. Separate models tested for GWI associations with ASCVD and risk factors documented in the electronic health record. GWI was associated with self-reported hypertension (adjusted odds ratio [aOR], 1.67 [95% CI, 1.18-2.36]), hyperlipidemia (aOR, 1.46 [95% CI, 1.03-2.05]), and ASCVD (aOR, 2.65 [95% CI, 1.56-4.51]). In the subset of veterans with electronic health record data, GWI was associated with documented diabetes (aOR, 2.34 [95% CI, 1.43-3.82]) and hypertension (aOR, 2.84 [95% CI, 1.92-4.20]). Hyperlipidemia and hypertension served as partial mediators of the association between GWI and self-reported ASCVD.

CONCLUSIONS

Gulf War veterans with GWI had higher odds of hyperlipidemia, hypertension, diabetes, and ASCVD compared with Gulf War veterans without GWI. Further examination of the mechanisms underlying this association, including a possible shared exposure-related mechanism, is necessary.

GWAS Meta-Analysis of Suicide Attempt: Identification of 12 Genome-Wide Significant Loci and Implication of Genetic Risks for Specific Health Factors

OBJECTIVE

Suicidal behavior is heritable and is a major cause of death worldwide. Two large-scale genome-wide association studies (GWASs) recently discovered and cross-validated genome-wide significant (GWS) loci for suicide attempt (SA). The present study leveraged the genetic cohorts from both studies to conduct the largest GWAS meta-analysis of SA to date. Multi-ancestry and admixture-specific meta-analyses were conducted within groups of significant African, East Asian, and European ancestry admixtures.

METHODS

This study comprised 22 cohorts, including 43,871 SA cases and 915,025 ancestry-matched controls. Analytical methods across multi-ancestry and individual ancestry admixtures included inverse variance-weighted fixed-effects meta-analyses, followed by gene, gene-set, tissue-set, and drug-target enrichment, as well as summary-data-based Mendelian randomization with brain expression quantitative trait loci data, phenome-wide genetic correlation, and genetic causal proportion analyses.

RESULTS

Multi-ancestry and European ancestry admixture GWAS meta-analyses identified 12 risk loci at p values <5×10-8. These loci were mostly intergenic and implicated DRD2, SLC6A9, FURIN, NLGN1, SOX5, PDE4B, and CACNG2. The multi-ancestry SNP-based heritability estimate of SA was 5.7% on the liability scale (SE=0.003, p=5.7×10-80). Significant brain tissue gene expression and drug set enrichment were observed. There was shared genetic variation of SA with attention deficit hyperactivity disorder, smoking, and risk tolerance after conditioning SA on both major depressive disorder and posttraumatic stress disorder. Genetic causal proportion analyses implicated shared genetic risk for specific health factors.

CONCLUSIONS

This multi-ancestry analysis of suicide attempt identified several loci contributing to risk and establishes significant shared genetic covariation with clinical phenotypes. These findings provide insight into genetic factors associated with suicide attempt across ancestry admixture populations, in veteran and civilian populations, and in attempt versus death.

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.

Declining autozygosity over time: An exploration in over 1 million individuals from three diverse cohorts

Previous studies have hypothesized that autozygosity is decreasing over generational time. However, these studies were limited to relatively small samples (n < 11,000) lacking in diversity, which may limit the generalizability of their findings. We present data that partially support this hypothesis from three large cohorts of diverse ancestries, two from the US (All of Us, n = 82,474; the Million Veteran Program, n = 622,497) and one from the UK (UK Biobank, n = 380,899). Our results from a mixed-effect meta-analysis demonstrate an overall trend of decreasing autozygosity over generational time (meta-analyzed slope = -0.029, SE = 0.009, p = 6.03e-4). On the basis of our estimates, we would predict FROH to decline 0.29% for every 20-year increase in birth year. We determined that a model including an ancestry-by-country interaction term fit the data best, indicating that ancestry differences in this trend differ by country. We found further evidence to suggest a difference between the US and UK cohorts by meta-analyzing within country, observing a significant negative estimate in the US cohorts (meta-analyzed slope = -0.058, SE = 0.015, p = 1.50e-4) but a non-significant estimate in the UK (meta-analyzed slope = -0.001, SE = 0.008, p = 0.945). The association between autozygosity and birth year was substantially attenuated when accounting for educational attainment and income (meta-analyzed slope = -0.011, SE = 0.008, p = 0.167), suggesting they may partially account for decreasing autozygosity over time. Overall, we demonstrate decreasing autozygosity over time in a large, modern sample and speculate that this trend can be attributed to increases in urbanization and panmixia and differences in sociodemographic processes lead to country-specific differences in the rate of decline.

Self-reported gastrointestinal disorders among veterans with gulf war illness with and without posttraumatic stress disorder

BACKGROUND

Gulf War Illness (GWI) is a chronic, multi-symptom disorder affecting 25%-32% of Gulf War veterans. Veterans with GWI disproportionately suffer from gastrointestinal (GI) disorders. Given the increasing evidence supporting a gut-brain axis, we explore the relationship between post-traumatic stress disorder (PTSD), GWI, and self-reported GI disorders among GW veterans.

METHODS

Veterans from the Gulf War Era Cohort and Biorepository responded to a mail-based survey (N = 1058). They were stratified by GWI (Centers for Disease Control definition) and PTSD status. This yielded three groups: GWI-, GWI+/PTSD-, and GWI+/PTSD+. Multivariable logistic regression adjusting for demographic and military characteristics examined associations between GWI/PTSD groups and GI disorders. Results were expressed as adjusted odds ratios (aOR) with 95% confidence intervals (95% CI).

KEY RESULTS

The most frequently reported GI disorders were irritable bowel syndrome (IBS), gastroesophageal reflux disease (GERD), and colon polyps (CP). The GWI+/PTSD+ group had a higher odds of these disorders than the GWI+/PTSD- group (aOR_IBS  = 3.12, 95% CI: 1.93-5.05; aOR_GERD  = 2.04, 95% CI: 1.44-2.90; aOR_CP  = 1.85, 95% CI: 1.23-2.80), which had a higher odds of these disorders than the GWI- group (aOR_IBS  = 4.38, 95% CI: 1.55-12.36; aOR_GERD  = 2.51 95% CI: 1.63-3.87; aOR_CP  = 2.57, 95% CI: 1.53-4.32).

CONCLUSIONS & INFERENCES

GW veterans with GWI and PTSD have significantly higher odds of specific self-reported GI disorders than the other groups. Given the known bidirectional influences of the gut and brain, these veterans may benefit from a holistic healthcare approach that considers biopsychosocial contributors to the assessment and management of disease.

Genetic and epigenetic signatures associated with plasma oxytocin levels in children and adolescents with autism spectrum disorder

Oxytocin (OT), the brain's most abundant neuropeptide, plays an important role in social salience and motivation. Clinical trials of the efficacy of OT in autism spectrum disorder (ASD) have reported mixed results due in part to ASD's complex etiology. We investigated whether genetic and epigenetic variation contribute to variable endogenous OT levels that modulate sensitivity to OT therapy. To carry out this analysis, we integrated genome-wide profiles of DNA-methylation, transcriptional activity, and genetic variation with plasma OT levels in 290 participants with ASD enrolled in a randomized controlled trial of OT. Our analysis identified genetic variants with novel association with plasma OT, several of which reside in known ASD risk genes. We also show subtle but statistically significant association of plasma OT levels with peripheral transcriptional activity and DNA-methylation profiles across several annotated gene sets. These findings broaden our understanding of the effects of the peripheral oxytocin system and provide novel genetic candidates for future studies to decode the complex etiology of ASD and its interaction with OT signaling and OT-based interventions. LAY SUMMARY: Oxytocin (OT) is an abundant chemical produced by neurons that plays an important role in social interaction and motivation. We investigated whether genetic and epigenetic factors contribute to variable OT levels in the blood. To this, we integrated genetic, gene expression, and non-DNA regulated (epigenetic) signatures with blood OT levels in 290 participants with autism enrolled in an OT clinical trial. We identified genetic association with plasma OT, several of which reside in known autism risk genes. We also show statistically significant association of plasma OT levels with gene expression and epigenetic across several gene pathways. These findings broaden our understanding of the factors that influence OT levels in the blood for future studies to decode the complex presentation of autism and its interaction with OT and OT-based treatment.

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

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

Effect of behavioral weight-loss program on biomarkers of cardiometabolic disease risk: Heart Health Study randomized trial

OBJECTIVE

This study aimed to determine whether novel biomarkers of cardiometabolic health improve in response to a 12-month behavioral weight-loss intervention and to compare benefits of diet alone with diet plus physical activity for these biomarkers.

METHODS

Participants (N = 374) were randomized to either diet alone (DIET), diet plus 150 min/wk of prescribed moderate-intensity physical activity (DIET + PA150), or diet plus 250 min/wk of prescribed moderate-intensity physical activity (DIET + PA250). Biomarker concentrations were determined using nuclear magnetic resonance spectroscopy. Mixed models assessed for a time effect, group effect, or group by time interaction.

RESULTS

All groups significantly improved body weight (time: p < 0.0001), Lipoprotein Insulin Resistance Index score (time: p < 0.0001), Diabetes Risk Index score (time: p < 0.0001), branched-chain amino acid concentration (time: p < 0.0001), and GlycA concentration (time: p < 0.0001), with no group effect or group by time interactions.

CONCLUSIONS

All intervention groups prompted a notable beneficial change among biomarkers of insulin resistance and cardiometabolic health. However, the addition of at least moderate-intensity physical activity to a diet-only intervention did not provide any additional benefit. These findings highlight that an average weight loss of approximately 10% profoundly impacts biomarkers of insulin resistance and cardiometabolic disease in adults with overweight or obesity.

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

Importance

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

Objective

To identify novel, replicable genomic risk loci for SITB.

Design, Setting, and Participants

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

Main Outcome and Measures

SITB.

Results

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

Conclusions and Relevance

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

Association of Gulf War Illness with Characteristics in Deployed vs. Non-Deployed Gulf War Era Veterans in the Cooperative Studies Program 2006/Million Veteran Program 029 Cohort: A Cross-Sectional Analysis

Gulf War Illness (GWI), a chronic multisymptom illness with a complex and uncertain etiology and pathophysiology, is highly prevalent among veterans deployed to the 1990-1991 GW. We examined how GWI phenotypes varied by demographic and military characteristics among GW-era veterans. Data were from the VA's Cooperative Studies Program 2006/Million Veteran Program (MVP) 029 cohort, Genomics of GWI. From June 2018 to March 2019, 109,976 MVP enrollees (out of a total of over 676,000) were contacted to participate in the 1990-1991 GW-era Survey. Of 109,976 eligible participants, 45,169 (41.1%) responded to the 2018-2019 survey, 35,902 respondents met study inclusion criteria, 13,107 deployed to the GW theater. GWI phenotypes were derived from Kansas (KS) and Centers for Disease Control and Prevention (CDC) GWI definitions: (a) KS Symptoms (KS Sym+), (b) KS GWI (met symptom criteria and without exclusionary health conditions) [KS GWI: Sym+/Dx-], (c) CDC GWI and (d) CDC GWI Severe. The prevalence of each phenotype was 67.1% KS Sym+, 21.5% KS Sym+/Dx-, 81.1% CDC GWI, and 18.6% CDC GWI severe. These findings affirm the persistent presence of GWI among GW veterans providing a foundation for further exploration of biological and environmental underpinnings of this condition.

Establishment of multi-stage intravenous self-administration paradigms in mice

Genetically tractable animal models provide needed strategies to resolve the biological basis of drug addiction. Intravenous self-administration (IVSA) is the gold standard for modeling psychostimulant and opioid addiction in animals, but technical limitations have precluded the widespread use of IVSA in mice. Here, we describe IVSA paradigms for mice that capture the multi-stage nature of the disorder and permit predictive modeling. In these paradigms, C57BL/6J mice with long-standing indwelling jugular catheters engaged in cocaine- or remifentanil-associated lever responding that was fixed ratio-dependent, dose-dependent, extinguished by withholding the drug, and reinstated by the presentation of drug-paired cues. The application of multivariate analysis suggested that drug taking in both paradigms was a function of two latent variables we termed incentive motivation and discriminative control. Machine learning revealed that vulnerability to drug seeking and relapse were predicted by a mouse's a priori response to novelty, sensitivity to drug-induced locomotion, and drug-taking behavior. The application of these behavioral and statistical-analysis approaches to genetically-engineered mice will facilitate the identification of neural circuits driving addiction susceptibility and relapse and focused therapeutic development.

The association of accelerated epigenetic age with all-cause mortality in cardiac catheterization patients as mediated by vascular and cardiometabolic outcomes

BACKGROUND

Epigenetic age is a DNA methylation-based biomarker of aging that is accurate across the lifespan and a range of cell types. The difference between epigenetic age and chronological age, termed age acceleration (AA), is a strong predictor of lifespan and healthspan. The predictive capabilities of AA for all-cause mortality have been evaluated in the general population; however, its utility is less well evaluated in those with chronic conditions. Additionally, the pathophysiologic pathways whereby AA predicts mortality are unclear. We hypothesized that AA predicts mortality in individuals with underlying cardiovascular disease; and the association between AA and mortality is mediated, in part, by vascular and cardiometabolic measures.

METHODS

We evaluated 562 participants in an urban, three-county area of central North Carolina from the CATHGEN cohort, all of whom received a cardiac catheterization procedure. We analyzed three AA biomarkers, Horvath epigenetic age acceleration (HAA), phenotypic age acceleration (PhenoAA), and Grim age acceleration (GrimAA), by Cox regression models, to assess whether AAs were associated with all-cause mortality. We also evaluated if these associations were mediated by vascular and cardiometabolic outcomes, including left ventricular ejection fraction (LVEF), blood cholesterol concentrations, angiopoietin-2 (ANG2) protein concentration, peripheral artery disease, coronary artery disease, diabetes, and hypertension. The total effect, direct effect, indirect effect, and percentage mediated were estimated using pathway mediation tests with a regression adjustment approach.

RESULTS

PhenoAA (HR = 1.05, P < 0.0001), GrimAA (HR = 1.10, P < 0.0001) and HAA (HR = 1.03, P = 0.01) were all associated with all-cause mortality. The association of mortality and PhenoAA was partially mediated by ANG2, a marker of vascular function (19.8%, P = 0.016), and by diabetes (8.2%, P = 0.043). The GrimAA-mortality association was mediated by ANG2 (12.3%, P = 0.014), and showed weaker evidence for mediation by LVEF (5.3%, P = 0.065).

CONCLUSIONS

Epigenetic age acceleration remains strongly predictive of mortality even in individuals already burdened with cardiovascular disease. Mortality associations were mediated by ANG2, which regulates endothelial permeability and angiogenic functions, suggesting that specific vascular pathophysiology may link accelerated epigenetic aging with increased mortality risks.

Exposures to low-levels of fine particulate matter are associated with acute changes in heart rate variability, cardiac repolarization, and circulating blood lipids in coronary artery disease patients

Exposure to air pollution is a major risk factor for cardiovascular disease, disease risk factors, and mortality. Specifically, particulate matter (PM), and to some extent ozone, are contributors to these effects. In addition, exposures to these pollutants may be especially dangerous for susceptible populations. In this repeated-visit panel study, cardiovascular markers were collected from thirteen male participants with stable coronary artery disease. For 0-4 days prior to the health measurement collections, daily concentrations of fine PM (PM_2.5) and ozone were obtained from local central monitoring stations located near the participant's homes. Then, single (PM_2.5) and two-pollutant (PM_2.5 and ozone) models were used to assess whether there were short-term changes in cardiovascular health markers. Per interquartile range increase in PM_2.5, there were decrements in several heart rate variability metrics, including the standard deviation of the normal-to-normal intervals (lag 3, -5.8%, 95% confidence interval (CI) = -11.5, 0.3) and root-mean squared of successive differences (five day moving average, -8.1%, 95% CI = -15.0, -0.7). In addition, increases in PM_2.5 were also associated with changes in P complexity (lag 1, 4.4%, 95% CI = 0.5, 8.5), QRS complexity (lag 1, 4.9%, 95% CI = 1.4, 8.5), total cholesterol (five day moving average, -2.1%, 95% CI = -4.1, -0.1), and high-density lipoprotein cholesterol (lag 2, -1.6%, 95% CI = -3.1, -0.1). Comparisons to our previously published work on ozone were conducted. We found that ozone affected inflammation and endothelial function, whereas PM_2.5 influenced heart rate variability, repolarization, and lipids. All the health changes from these two studies were found at concentrations below the United States Environmental Protection Agency's National Ambient Air Quality Standards. Our results imply clear differences in the cardiovascular outcomes observed with exposure to the two ubiquitous air pollutants PM_2.5 and ozone; this observation suggests different mechanisms of toxicity for these exposures.

Longitudinal assessment of colonoscopy adverse events in the prospective Cooperative Studies Program no. 380 colorectal cancer screening and surveillance cohort

BACKGROUND AND AIMS

Data are limited regarding colonoscopy risk during long-term, programmatic colorectal cancer screening and follow-up. We aimed to describe adverse events during follow-up in a colonoscopy screening program after the baseline examination and examine factors associated with increased risk.

METHODS

Cooperative Studies Program no. 380 includes 3121 asymptomatic veterans aged 50 to 75 years who underwent screening colonoscopy between 1994 and 1997. Periprocedure adverse events requiring significant intervention were defined as major events (other events were minor) and were tracked during follow-up for at least 10 years. Multivariable odds ratios (ORs) were calculated for factors associated with risk of follow-up adverse events.

RESULTS

Of 3727 follow-up examinations in 1983 participants, adverse events occurred in 105 examinations (2.8%) in 93 individuals, including 22 major and 87 minor events (examinations may have had >1 event). Incidence of major events (per 1000 examinations) remained relatively stable over time, with 6.1 events at examination 2, 4.8 at examination 3, and 7.2 at examination 4. Examinations with major events included 1 perforation, 3 GI bleeds requiring intervention, and 17 cardiopulmonary events. History of prior colonoscopic adverse events was associated with increased risk of events (major or minor) during follow-up (OR, 2.7; 95% confidence interval, 1.6-4.6).

CONCLUSIONS

Long-term programmatic screening and surveillance was safe, as major events were rare during follow-up. However, serious cardiopulmonary events were the most common major events. These results highlight the need for detailed assessments of comorbid conditions during routine clinical practice, which could help inform individual decisions regarding the utility of ongoing colonoscopy follow-up.

Large-scale genome-wide association study of coronary artery disease in genetically diverse populations

We report a genome-wide association study (GWAS) of coronary artery disease (CAD) incorporating nearly a quarter of a million cases, in which existing studies are integrated with data from cohorts of white, Black and Hispanic individuals from the Million Veteran Program. We document near equivalent heritability of CAD across multiple ancestral groups, identify 95 novel loci, including nine on the X chromosome, detect eight loci of genome-wide significance in Black and Hispanic individuals, and demonstrate that two common haplotypes at the 9p21 locus are responsible for risk stratification in all populations except those of African origin, in which these haplotypes are virtually absent. Moreover, in the largest GWAS for angiographically derived coronary atherosclerosis performed to date, we find 15 loci of genome-wide significance that robustly overlap with established loci for clinical CAD. Phenome-wide association analyses of novel loci and polygenic risk scores (PRSs) augment signals related to insulin resistance, extend pleiotropic associations of these loci to include smoking and family history, and precisely document the markedly reduced transferability of existing PRSs to Black individuals. Downstream integrative analyses reinforce the critical roles of vascular endothelial, fibroblast, and smooth muscle cells in CAD susceptibility, but also point to a shared biology between atherosclerosis and oncogenesis. This study highlights the value of diverse populations in further characterizing the genetic architecture of CAD.

Ten or More Cumulative Lifetime Adenomas Are Associated with Increased Risk for Advanced Neoplasia and Colorectal Cancer

BACKGROUND

Colorectal cancer (CRC) screening guidelines recommend frequent colonoscopies and consideration of genetic testing in individuals with ≥10 cumulative adenomas. However, it is unclear how these guidelines apply to routine practice.

AIMS

We estimated the proportion of participants found to have ≥10 cumulative adenomas in a screening population and described their outcomes of advanced neoplasia (AN), CRC, and extra-colonic malignancy.

METHODS

We performed a secondary analysis of VA CSP#380, which includes 3121 veterans aged 50-75 who were followed up to 10 years after screening colonoscopy. We calculated the cumulative risk of ≥10 cumulative adenomas by Kaplan-Meier method. We compared baseline risk factors in those with and without ≥10 cumulative adenomas as well as the risk for AN (adenoma ≥1 cm, villous adenoma or high-grade dysplasia, or CRC) and extra-colonic malignancy by multivariate logistic regression.

RESULTS

The cumulative risk of ≥10 cumulative adenomas over 10.5 years was 6.51% (95% CI 4.38%-9.62%). Age 60-69 or 70-75 at baseline colonoscopy was the only factors associated with the finding of ≥10 cumulative adenomas. Compared to those with 0-9 cumulative adenomas, participants with ≥10 cumulative adenomas were more likely to have had AN (OR 17.03; 95% CI 9.41-30.84), including CRC (OR 7.00; 95% CI 2.84-17.28), but not extra-colonic malignancies.

CONCLUSIONS

Approximately 6.5% of participants in this screening population were found to have ≥10 cumulative adenomas over 10.5 years, which was uncommon before age 60. These participants were found to have AN and CRC significantly more often compared to those with lower cumulative adenomas.

Rheumatoid arthritis T cell and muscle oxidative metabolism associate with exercise-induced changes in cardiorespiratory fitness

Rheumatoid arthritis (RA) T cells drive autoimmune features via metabolic reprogramming that reduces oxidative metabolism. Exercise training improves cardiorespiratory fitness (i.e., systemic oxidative metabolism) and thus may impact RA T cell oxidative metabolic function. In this pilot study of RA participants, we took advantage of heterogeneous responses to a high-intensity interval training (HIIT) exercise program to identify relationships between improvements in cardiorespiratory fitness with changes in peripheral T cell and skeletal muscle oxidative metabolism. In 12 previously sedentary persons with seropositive RA, maximal cardiopulmonary exercise tests, fasting blood, and vastus lateralis biopsies were obtained before and after 10 weeks of HIIT. Following HIIT, improvements in RA cardiorespiratory fitness were associated with changes in RA CD4 + T cell basal and maximal respiration and skeletal muscle carnitine acetyltransferase (CrAT) enzyme activity. Further, changes in CD4 + T cell respiration were associated with changes in naïve CD4 + CCR7 + CD45RA + T cells, muscle CrAT, and muscle medium-chain acylcarnitines and fat oxidation gene expression profiles. In summary, modulation of cardiorespiratory fitness and molecular markers of skeletal muscle oxidative metabolism during exercise training paralleled changes in T cell metabolism. Exercise training that improves RA cardiorespiratory fitness may therefore be valuable in managing pathologically related immune and muscle dysfunction.Trial registration: ClinicalTrials.gov, NCT02528344. Registered on 19 August 2015.

Sex-dimorphic gene effects on survival outcomes in people with coronary artery disease

Background:

Ischemic coronary heart disease (IHD) is the leading cause of death worldwide. Genetic variation is presumed to be a major factor underlying sex differences for IHD events, including mortality. The purpose of this study was to identify sex-specific candidate genes associated with all-cause mortality among people diagnosed with coronary artery disease (CAD).

Methods:

We performed a sex-stratified, exploratory genome-wide association (GWAS) screen using existing data from CAD-diagnosed males (n = 510) and females (n = 174) who reported European ancestry from the Duke Catheterization Genetics biorepository. Extant genotype data for 785,945 autosomal SNPs generated with the Human Omni1-Quad BeadChip (Illumina, CA, USA) were analyzed using an additive inheritance model. We estimated instantaneous risk of all-cause mortality by genotype groups across the 11-year follow-up using Cox multivariate regression, covarying for age and genomic ancestry.

Results:

The top GWAS hits associated with all-cause mortality among people with CAD included 8 SNPs among males and 15 among females (p = 1 × 10⁻⁶ or 10⁻⁷), adjusted for covariates. Cross-sex comparisons revealed distinct candidate genes. Biologically relevant candidates included rs9932462 (EMP2/TEKT5) and rs2835913 (KCNJ6) among males and rs7217169 (RAP1GAP2), rs8021816 (PRKD1), rs8133010 (PDE9A), and rs12145981 (LPGAT1) among females.

Conclusions:

We report 20 sex-specific candidate genes having suggestive association with all-cause mortality among CAD-diagnosed subjects. Findings demonstrate proof of principle for identifying sex-associated genetic factors that may help explain differential mortality risk in people with CAD. Replication and meta-analyses in larger studies with more diverse samples will strengthen future work in this area.

Screening Colonoscopy Findings Are Associated With Noncolorectal Cancer Mortality

INTRODUCTION

Controversy exists regarding the impact of various risk factors on noncolorectal cancer (CRC) mortality in healthy screening populations. We examined the impact of known CRC risk factors, including baseline colonoscopy findings, on non-CRC mortality in a screening population.

METHODS

Cooperative Studies Program (CSP) #380 is comprised of 3,121 veterans aged 50-75 years who underwent screening colonoscopy from 1994 to 97 and were then followed for at least 10 years or until death. Hazard ratios (HRs) for risk factors on non-CRC mortality were estimated by multivariate Cox proportional hazards.

RESULTS

Current smoking (HR 2.12, 95% confidence interval [CI] 1.78-2.52, compared with nonsmokers) and physical activity (HR 0.89, 95% CI 0.84-0.93) were the modifiable factors most associated with non-CRC mortality in CSP#380. In addition, compared with no neoplasia at baseline colonoscopy, non-CRC mortality was higher in participants with ≥3 small adenomas (HR 1.43, 95% CI 1.06-1.94), advanced adenomas (HR 1.32, 95% CI 0.99-1.75), and CRC (HR 2.95, 95% CI 0.98-8.85). Those with 1-2 small adenomas were not at increased risk for non-CRC mortality (HR 1.15, 95% CI 0.94-1.4).

DISCUSSION

In a CRC screening population, known modifiable risk factors were significantly associated with 10-year non-CRC mortality. Furthermore, those who died from non-CRC causes within 10 years were more likely to have had high-risk findings at baseline colonoscopy. These results suggest that advanced colonoscopy findings may be a risk marker of poor health outcomes. Integrated efforts are needed to motivate healthy lifestyle changes during CRC screening, particularly in those with high-risk colonoscopy findings and unaddressed risk factors.

A genome-wide association study of suicide attempts in the million veterans program identifies evidence of pan-ancestry and ancestry-specific risk loci

To identify pan-ancestry and ancestry-specific loci associated with attempting suicide among veterans, we conducted a genome-wide association study (GWAS) of suicide attempts within a large, multi-ancestry cohort of U.S. veterans enrolled in the Million Veterans Program (MVP). Cases were defined as veterans with a documented history of suicide attempts in the electronic health record (EHR; N = 14,089) and controls were defined as veterans with no documented history of suicidal thoughts or behaviors in the EHR (N = 395,064). GWAS was performed separately in each ancestry group, controlling for sex, age and genetic substructure. Pan-ancestry risk loci were identified through meta-analysis and included two genome-wide significant loci on chromosomes 20 (p = 3.64 × 10-9) and 1 (p = 3.69 × 10-8). A strong pan-ancestry signal at the Dopamine Receptor D2 locus (p = 1.77 × 10-7) was also identified and subsequently replicated in a large, independent international civilian cohort (p = 7.97 × 10-4). Additionally, ancestry-specific genome-wide significant loci were also detected in African-Americans, European-Americans, Asian-Americans, and Hispanic-Americans. Pathway analyses suggested over-representation of many biological pathways with high clinical significance, including oxytocin signaling, glutamatergic synapse, cortisol synthesis and secretion, dopaminergic synapse, and circadian rhythm. These findings confirm that the genetic architecture underlying suicide attempt risk is complex and includes both pan-ancestry and ancestry-specific risk loci. Moreover, pathway analyses suggested many commonly impacted biological pathways that could inform development of improved therapeutics for suicide prevention.

Dissecting the Shared Genetic Architecture of Suicide Attempt, Psychiatric Disorders, and Known Risk Factors

BACKGROUND

Suicide is a leading cause of death worldwide, and nonfatal suicide attempts, which occur far more frequently, are a major source of disability and social and economic burden. Both have substantial genetic etiology, which is partially shared and partially distinct from that of related psychiatric disorders.

METHODS

We conducted a genome-wide association study (GWAS) of 29,782 suicide attempt (SA) cases and 519,961 controls in the International Suicide Genetics Consortium (ISGC). The GWAS of SA was conditioned on psychiatric disorders using GWAS summary statistics via multitrait-based conditional and joint analysis, to remove genetic effects on SA mediated by psychiatric disorders. We investigated the shared and divergent genetic architectures of SA, psychiatric disorders, and other known risk factors.

RESULTS

Two loci reached genome-wide significance for SA: the major histocompatibility complex and an intergenic locus on chromosome 7, the latter of which remained associated with SA after conditioning on psychiatric disorders and replicated in an independent cohort from the Million Veteran Program. This locus has been implicated in risk-taking behavior, smoking, and insomnia. SA showed strong genetic correlation with psychiatric disorders, particularly major depression, and also with smoking, pain, risk-taking behavior, sleep disturbances, lower educational attainment, reproductive traits, lower socioeconomic status, and poorer general health. After conditioning on psychiatric disorders, the genetic correlations between SA and psychiatric disorders decreased, whereas those with nonpsychiatric traits remained largely unchanged.

CONCLUSIONS

Our results identify a risk locus that contributes more strongly to SA than other phenotypes and suggest a shared underlying biology between SA and known risk factors that is not mediated by psychiatric disorders.

Determinants of Dropout from and Variation in Adherence to an Exercise Intervention: The STRRIDE Randomized Trials

Purpose

This study aimed to characterize the timing and self-reported determinants of exercise dropout among sedentary adults with overweight or obesity. We also sought to explore variations in adherence among individuals who completed a 6- to 8-month structured exercise intervention.

Methods

A total of 947 adults with dyslipidemia [STRRIDE I, STRRIDE AT/RT] or prediabetes [STRRIDE-PD] were enrolled to either control or to one of 10 exercise interventions, ranging from doses of 8-23 kcal/kg/week; intensities of 50%-75% V̇O2 peak; and durations of 6-8 months. Two groups included resistance training and one included dietary intervention (7% weight loss goal). Dropout was defined as an individual who withdrew from the study due a variety of determinants. Timing of intervention dropout was defined as the last session attended and categorized into phases. Exercise training adherence was calculated by dividing weekly minutes or total sets of exercise completed by weekly minutes or total sets of exercise prescribed. General linear models were used to characterize the associations between timing of dropout and determinant category.

Results

Compared to exercise intervention completers (n=652), participants who dropped out (n=295) were on average non-white (98% vs. 80%, p<0.01), had higher body mass index (31.0 kg/m2 vs. 30.2 kg/m2; p<0.01), and were less fit at baseline (25.0 mg/kg/min vs. 26.7 ml/kg/min, p<0.01). Of those who dropped out, 67% did so prior to the start of or while ramping up to the prescribed exercise volume and intensity. The most commonly reported reason for dropout was lack of time (40%). Notably, among individuals who completed the ramp training period, subsequent exercise intervention adherence did not waiver over the ensuing 6-8 months of training.

Conclusion

These findings are some of the first to delineate associations between the timing of dropout and dropout determinants, providing guidance to future exercise interventions to better support individuals at-risk for dropout.

Associations between neighborhood socioeconomic cluster and hypertension, diabetes, myocardial infarction, and coronary artery disease within a cohort of cardiac catheterization patients

Background

Neighborhood-level socioeconomic status (SES) is associated with health outcomes, including cardiovascular disease and diabetes, but these associations are rarely studied across large, diverse populations.

Methods

We used Ward’s Hierarchical clustering to define eight neighborhood clusters across North Carolina using 11 census-based indicators of SES, race, housing, and urbanicity and assigned 6992 cardiac catheterization patients at Duke University Hospital from 2001 to 2010 to clusters. We examined associations between clusters and coronary artery disease index > 23 (CAD), history of myocardial infarction, hypertension, and diabetes using logistic regression adjusted for age, race, sex, body mass index, region of North Carolina, distance to Duke University Hospital, and smoking status.

Results

Four clusters were urban, three rural, and one suburban higher-middle-SES (referent). We observed greater odds of myocardial infarction in all six clusters with lower or middle-SES. Odds of CAD were elevated in the rural cluster that was low-SES and plurality Black (OR 1.16, 95% CI 0.94-1.43) and in the rural cluster that was majority American Indian (OR 1.31, 95% CI 0.91-1.90). Odds of diabetes and hypertension were elevated in two urban and one rural low- and lower-middle SES clusters with large Black populations.

Conclusions

We observed higher prevalence of cardiovascular disease and diabetes in neighborhoods that were predominantly rural, low-SES, and non-White, highlighting the importance of public health and healthcare system outreach into these communities to promote cardiometabolic health and prevent and manage hypertension, diabetes and coronary artery disease.

Gene-Toxicant Interactions in Gulf War Illness: Differential Effects of the PON1 Genotype

About 25-35% of United States veterans who fought in the 1990-1991 Gulf War report several moderate or severe chronic systemic symptoms, defined as Gulf War illness (GWI). Thirty years later, there is little consensus on the causes or biological underpinnings of GWI. The Gulf War Era Cohort and Biorepository (GWECB) was designed to investigate genetic and environmental associations with GWI and consists of 1343 veterans. We investigate candidate gene-toxicant interactions that may be associated with GWI based on prior associations found in human and animal model studies, focusing on SNPs in or near ACHE, BCHE, and PON1 genes to replicate results from prior studies. SOD1 was also considered as a candidate gene. CDC Severe GWI, the primary outcome, was observed in 26% of the 810 deployed veterans included in this study. The interaction between the candidate SNP rs662 and pyridostigmine bromide (PB) pills was found to be associated with CDC Severe GWI. Interactions between PB pill exposure and rs3917545, rs3917550, and rs2299255, all in high linkage disequilibrium in PON1, were also associated with respiratory symptoms. These SNPs could point toward biological pathways through which GWI may develop, which could lead to biomarkers to detect GWI or to better treatment options for veterans with GWI.

Genomics of Gulf War Illness in U.S. Veterans Who Served during the 1990-1991 Persian Gulf War: Methods and Rationale for Veterans Affairs Cooperative Study #2006

Background: Approximately 697,000 members of the U.S. Armed Forces were deployed to the Persian Gulf in support of the 1990–1991 Persian Gulf War (GW). Subsequently, many deployed and some non-deployed veterans developed a chronic multi-symptom illness, now named Gulf War Illness (GWI). This manuscript outlines the methods and rationale for studying the genomics of GWI within the Million Veteran Program (MVP), a VA-based national research program that has linked medical records, surveys, and genomic data, enabling genome-wide association studies (GWASs). Methods: MVP participants who served in the military during the GW era were contacted by mail and invited to participate in the GWI study. A structured health questionnaire, based on a previously tested instrument, was also included in the mailing. Data on deployment locations and exposures, symptoms associated with GWI, clinical diagnoses, personal habits, and health care utilization were collected. Self-reported data will be augmented with chart reviews and structured international classification of disease codes, to classify participants by GWI case status. We will develop a phenotyping algorithm, based on two commonly used case definitions, to determine GWI status, and then conduct a nested case-control GWAS. Genetic variants associated with GWI will be investigated, and gene–gene and gene–environment interactions studied. The genetic overlap of GWI with, and causative mechanisms linking this illness to, other health conditions and the effects of genomic regulatory mechanisms on GWI risk will also be explored. Conclusions: The proposed initial GWAS described in this report will investigate the genomic underpinnings of GWI with a large sample size and state-of-the-art genomic analyses and phenotyping. The data generated will provide a rich and expansive foundation on which to build additional analyses.

Genome-Wide Variants Associated With Longitudinal Survival Outcomes Among Individuals With Coronary Artery Disease

Objective

Coronary artery disease (CAD) is an age-associated condition that greatly increases the risk of mortality. The purpose of this study was to identify gene variants associated with all-cause mortality among individuals with clinically phenotyped CAD using a genome-wide screening approach.

Approach and Results

We performed discovery (n = 684), replication (n = 1,088), and meta-analyses (N = 1,503) for association of genomic variants with survival outcome using secondary data from White participants with CAD from two GWAS sub-studies of the Duke Catheterization Genetics Biorepository. We modeled time from catheterization to death or last follow-up (median 7.1 years, max 12 years) using Cox multivariable regression analysis. Target statistical screening thresholds were p × 10^–8 for the discovery phase and Bonferroni-calculated p-values for the replication (p < 5.3 × 10^–4) and meta-analysis (p < 1.6 × 10^–3) phases. Genome-wide analysis of 785,945 autosomal SNPs revealed two SNPs (rs13007553 and rs587936) that had the same direction of effect across all three phases of the analysis, with suggestive p-value association in discovery and replication and significant meta-analysis association in models adjusted for clinical covariates. The rs13007553 SNP variant, LINC01250, which resides between MYTIL and EIPR1, conferred increased risk for all-cause mortality even after controlling for clinical covariates [HR 1.47, 95% CI 1.17–1.86, p(adj) = 1.07 × 10^–3 (discovery), p(adj) = 0.03 (replication), p(adj) = 9.53 × 10^–5 (meta-analysis)]. MYT1L is involved in neuronal differentiation. TSSC1 is involved in endosomal recycling and is implicated in breast cancer. The rs587936 variant annotated to DAB2IP was associated with increased survival time [HR 0.65, 95% CI 0.51–0.83, p(adj) = 4.79 × 10^–4 (discovery), p(adj) = 0.02 (replication), p(adj) = 2.25 × 10^–5 (meta-analysis)]. DAB2IP is a ras/GAP tumor suppressor gene which is highly expressed in vascular tissue. DAB2IP has multiple lines of evidence for protection against atherosclerosis.

Conclusion

Replicated findings identified two candidate genes for further study regarding association with survival in high-risk CAD patients: novel loci LINC01250 (rs13007553) and biologically relevant candidate DAB2IP (rs587936). These candidates did not overlap with validated longevity candidate genes. Future research could further define the role of common variants in survival outcomes for people with CAD and, ultimately, improve longitudinal outcomes for these patients.

Genome-wide analysis identifies novel susceptibility loci for myocardial infarction

AIMS

While most patients with myocardial infarction (MI) have underlying coronary atherosclerosis, not all patients with coronary artery disease (CAD) develop MI. We sought to address the hypothesis that some of the genetic factors which establish atherosclerosis may be distinct from those that predispose to vulnerable plaques and thrombus formation.

METHODS AND RESULTS

We carried out a genome-wide association study for MI in the UK Biobank (n∼472 000), followed by a meta-analysis with summary statistics from the CARDIoGRAMplusC4D Consortium (n∼167 000). Multiple independent replication analyses and functional approaches were used to prioritize loci and evaluate positional candidate genes. Eight novel regions were identified for MI at the genome wide significance level, of which effect sizes at six loci were more robust for MI than for CAD without the presence of MI. Confirmatory evidence for association of a locus on chromosome 1p21.3 harbouring choline-like transporter 3 (SLC44A3) with MI in the context of CAD, but not with coronary atherosclerosis itself, was obtained in Biobank Japan (n∼165 000) and 16 independent angiography-based cohorts (n∼27 000). Follow-up analyses did not reveal association of the SLC44A3 locus with CAD risk factors, biomarkers of coagulation, other thrombotic diseases, or plasma levels of a broad array of metabolites, including choline, trimethylamine N-oxide, and betaine. However, aortic expression of SLC44A3 was increased in carriers of the MI risk allele at chromosome 1p21.3, increased in ischaemic (vs. non-diseased) coronary arteries, up-regulated in human aortic endothelial cells treated with interleukin-1β (vs. vehicle), and associated with smooth muscle cell migration in vitro.

CONCLUSIONS

A large-scale analysis comprising ∼831 000 subjects revealed novel genetic determinants of MI and implicated SLC44A3 in the pathophysiology of vulnerable plaques.

An atlas connecting shared genetic architecture of human diseases and molecular phenotypes provides insight into COVID-19 susceptibility

BACKGROUND

While genome-wide associations studies (GWAS) have successfully elucidated the genetic architecture of complex human traits and diseases, understanding mechanisms that lead from genetic variation to pathophysiology remains an important challenge. Methods are needed to systematically bridge this crucial gap to facilitate experimental testing of hypotheses and translation to clinical utility.

RESULTS

Here, we leveraged cross-phenotype associations to identify traits with shared genetic architecture, using linkage disequilibrium (LD) information to accurately capture shared SNPs by proxy, and calculate significance of enrichment. This shared genetic architecture was examined across differing biological scales through incorporating data from catalogs of clinical, cellular, and molecular GWAS. We have created an interactive web database (interactive Cross-Phenotype Analysis of GWAS database (iCPAGdb)) to facilitate exploration and allow rapid analysis of user-uploaded GWAS summary statistics. This database revealed well-known relationships among phenotypes, as well as the generation of novel hypotheses to explain the pathophysiology of common diseases. Application of iCPAGdb to a recent GWAS of severe COVID-19 demonstrated unexpected overlap of GWAS signals between COVID-19 and human diseases, including with idiopathic pulmonary fibrosis driven by the DPP9 locus. Transcriptomics from peripheral blood of COVID-19 patients demonstrated that DPP9 was induced in SARS-CoV-2 compared to healthy controls or those with bacterial infection. Further investigation of cross-phenotype SNPs associated with both severe COVID-19 and other human traits demonstrated colocalization of the GWAS signal at the ABO locus with plasma protein levels of a reported receptor of SARS-CoV-2, CD209 (DC-SIGN). This finding points to a possible mechanism whereby glycosylation of CD209 by ABO may regulate COVID-19 disease severity.

CONCLUSIONS

Thus, connecting genetically related traits across phenotypic scales links human diseases to molecular and cellular measurements that can reveal mechanisms and lead to novel biomarkers and therapeutic approaches. The iCPAGdb web portal is accessible at http://cpag.oit.duke.edu and the software code at https://github.com/tbalmat/iCPAGdb .

Epigenome-wide association study of kidney function identifies trans-ethnic and ethnic-specific loci

BACKGROUND

DNA methylation (DNAm) is associated with gene regulation and estimated glomerular filtration rate (eGFR), a measure of kidney function. Decreased eGFR is more common among US Hispanics and African Americans. The causes for this are poorly understood. We aimed to identify trans-ethnic and ethnic-specific differentially methylated positions (DMPs) associated with eGFR using an agnostic, genome-wide approach.

METHODS

The study included up to 5428 participants from multi-ethnic studies for discovery and 8109 participants for replication. We tested the associations between whole blood DNAm and eGFR using beta values from Illumina 450K or EPIC arrays. Ethnicity-stratified analyses were performed using linear mixed models adjusting for age, sex, smoking, and study-specific and technical variables. Summary results were meta-analyzed within and across ethnicities. Findings were assessed using integrative epigenomics methods and pathway analyses.

RESULTS

We identified 93 DMPs associated with eGFR at an FDR of 0.05 and replicated 13 and 1 DMPs across independent samples in trans-ethnic and African American meta-analyses, respectively. The study also validated 6 previously published DMPs. Identified DMPs showed significant overlap enrichment with DNase I hypersensitive sites in kidney tissue, sites associated with the expression of proximal genes, and transcription factor motifs and pathways associated with kidney tissue and kidney development.

CONCLUSIONS

We uncovered trans-ethnic and ethnic-specific DMPs associated with eGFR, including DMPs enriched in regulatory elements in kidney tissue and pathways related to kidney development. These findings shed light on epigenetic mechanisms associated with kidney function, bridging the gap between population-specific eGFR-associated DNAm and tissue-specific regulatory context.

Drebrin attenuates atherosclerosis by limiting smooth muscle cell transdifferentiation

AIMS

The F-actin-binding protein Drebrin inhibits smooth muscle cell (SMC) migration, proliferation and pro-inflammatory signaling. Therefore, we tested the hypothesis that Drebrin constrains atherosclerosis.

METHODS AND RESULTS

SM22-Cre+/Dbnflox/flox/Ldlr-/- (SMC-Dbn-/-/Ldlr-/-) and control mice (SM22-Cre+/Ldlr-/-, Dbnflox/flox/Ldlr-/-, and Ldlr-/-) were fed a Western diet for 14-20 weeks. Brachiocephalic arteries of SMC-Dbn-/-/Ldlr-/- mice exhibited 1.5- or 1.8-fold greater cross-sectional lesion area than control mice at 14 or 20 wk, respectively. Aortic atherosclerotic lesion surface area was 1.2-fold greater in SMC-Dbn-/-/Ldlr-/- mice. SMC-Dbn-/-/Ldlr-/- lesions comprised necrotic cores that were two-fold greater in size than those of control mice. Consistent with their bigger necrotic core size, lesions in SMC-Dbn-/- arteries also showed more transdifferentiation of SMCs to macrophage-like cells: 1.5- to 2.5-fold greater, assessed with BODIPY or with CD68, respectively. In vitro data were concordant: Dbn-/- SMCs had 1.7-fold higher levels of KLF4 and transdifferentiated to macrophage-like cells more readily than Dbnflox/flox SMCs upon cholesterol loading, as evidenced by greater up-regulation of CD68 and galectin-3. Adenovirally mediated Drebrin rescue produced equivalent levels of macrophage-like transdifferentiation in Dbn-/- and Dbnflox/flox SMCs. During early atherogenesis, SMC-Dbn-/-/Ldlr-/- aortas demonstrated 1.6-fold higher levels of reactive oxygen species than control mouse aortas. The 1.8-fold higher levels of Nox1 in Dbn-/- SMCs was reduced to WT levels with KLF4 silencing. Inhibition of Nox1 chemically or with siRNA produced equivalent levels of macrophage-like transdifferentiation in Dbn-/- and Dbnflox/flox SMCs.

CONCLUSIONS

We conclude that SMC Drebrin limits atherosclerosis by constraining SMC Nox1 activity and SMC transdifferentiation to macrophage-like cells.

TRANSLATIONAL PERSPECTIVE

Drebrin is abundantly expressed in vascular smooth muscle cells (SMCs) and is up-regulated in human atherosclerosis. A hallmark of atherosclerosis is the accumulation of foam cells that secrete pro-inflammatory cytokines and contribute to plaque instability. A large proportion of these foam cells in humans derive from SMCs. We found that SMC Drebrin limits atherosclerosis by reducing SMC transdifferentiation to macrophage-like foam cells in a manner dependent on Nox1 and KLF4. For this reason, strategies aimed at augmenting SMC Drebrin expression in atherosclerotic plaques may limit atherosclerosis progression and enhance plaque stability by bridling SMC-to-foam-cell transdifferentiation.

Th17 Immunity in the Colon Is Controlled by Two Novel Subsets of Colon-Specific Mononuclear Phagocytes

Intestinal immunity is coordinated by specialized mononuclear phagocyte populations, constituted by a diversity of cell subsets. Although the cell subsets constituting the mononuclear phagocyte network are thought to be similar in both small and large intestine, these organs have distinct anatomy, microbial composition, and immunological demands. Whether these distinctions demand organ-specific mononuclear phagocyte populations with dedicated organ-specific roles in immunity are unknown. Here we implement a new strategy to subset murine intestinal mononuclear phagocytes and identify two novel subsets which are colon-specific: a macrophage subset and a Th17-inducing dendritic cell (DC) subset. Colon-specific DCs and macrophages co-expressed CD24 and CD14, and surprisingly, both were dependent on the transcription factor IRF4. Novel IRF4-dependent CD14+CD24+ macrophages were markedly distinct from conventional macrophages and failed to express classical markers including CX3CR1, CD64 and CD88, and surprisingly expressed little IL-10, which was otherwise robustly expressed by all other intestinal macrophages. We further found that colon-specific CD14+CD24+ mononuclear phagocytes were essential for Th17 immunity in the colon, and provide definitive evidence that colon and small intestine have distinct antigen presenting cell requirements for Th17 immunity. Our findings reveal unappreciated organ-specific diversity of intestine-resident mononuclear phagocytes and organ-specific requirements for Th17 immunity.

Association between short-term exposure to ambient fine particulate matter and myocardial injury in the CATHGEN cohort

Exposure to fine particulate matter (PM_2.5) has been associated with a higher risk for coronary events. Elevated circulating cardiac troponins (cTn) are suggestive of myocardial injury in both ischemic and non-ischemic conditions. However, little is known about the association between PM_2.5 and cTn. In this study, we investigated short-term PM_2.5 effects on cardiac troponin T (cTnT), as well as N-terminal-pro brain natriuretic peptide (NT-pro BNP) and inflammatory biomarkers among cardiac catheterized participants. We analyzed 7444 plasma cTnT measurements in 2732 participants who presented to Duke University Hospital with myocardial infarction symptoms between 2001 and 2012, partly along with measurements of NT-pro BNP and inflammatory biomarkers. Daily PM_2.5 concentrations were predicted by a neural network-based hybrid model and were assigned to participants' residential addresses. We applied generalized estimating equations to assess associations of PM_2.5 with biomarker levels and the risk of a positive cTnT test (cTnT > 0.1 ng/mL). The median plasma cTnT concentration at presentation was 0.05 ng/mL and the prevalence of a positive cTnT test was 35.4%. For an interquartile range (7.6 μg/m³) increase in PM_2.5 on the previous day, cTnT concentrations increased by 7.7% (95% CI: 3.4-12.3) and the odds ratio of a positive cTnT test was 1.08 (1.01-1.16). Participants under 60 years (effect estimate: 15.2%; 95% CI: 7.4-23.5) or living in rural areas (12.3%; 95% CI: 4.8-20.3) were more susceptible. There was evidence for increases in fibrinogen and NT-pro BNP associated with elevated PM_2.5 on the concurrent and previous two days. Our study suggests that acute PM_2.5 exposure may elevate indicators of myocardial tissue damage. This finding substantiates the association of air pollution exposure with adverse cardiovascular events.

Characterizing chronological accumulation of comorbidities in healthy veterans: a computational approach

Understanding patient accumulation of comorbidities can facilitate healthcare strategy and personalized preventative care. We applied a directed network graph to electronic health record (EHR) data and characterized comorbidities in a cohort of healthy veterans undergoing screening colonoscopy. The Veterans Affairs Cooperative Studies Program #380 was a prospective longitudinal study of screening and surveillance colonoscopy. We identified initial instances of three-digit ICD-9 diagnoses for participants with at least 5 years of linked EHR history (October 1999 to December 2015). For diagnoses affecting at least 10% of patients, we calculated pairwise chronological relative risk (RR). iGraph was used to produce directed graphs of comorbidities with RR > 1, as well as summary statistics, key diseases, and communities. A directed graph based on 2210 patients visualized longitudinal development of comorbidities. Top hub (preceding) diseases included ischemic heart disease, inflammatory and toxic neuropathy, and diabetes. Top authority (subsequent) diagnoses were acute kidney failure and hypertensive chronic kidney failure. Four communities of correlated comorbidities were identified. Close analysis of top hub and authority diagnoses demonstrated known relationships, correlated sequelae, and novel hypotheses. Directed network graphs portray chronologic comorbidity relationships. We identified relationships between comorbid diagnoses in this aging veteran cohort. This may direct healthcare prioritization and personalized care.

An atlas connecting shared genetic architecture of human diseases and molecular phenotypes provides insight into COVID-19 susceptibility

While genome-wide associations studies (GWAS) have successfully elucidated the genetic architecture of complex human traits and diseases, understanding mechanisms that lead from genetic variation to pathophysiology remains an important challenge. Methods are needed to systematically bridge this crucial gap to facilitate experimental testing of hypotheses and translation to clinical utility. Here, we leveraged cross-phenotype associations to identify traits with shared genetic architecture, using linkage disequilibrium (LD) information to accurately capture shared SNPs by proxy, and calculate significance of enrichment. This shared genetic architecture was examined across differing biological scales through incorporating data from catalogs of clinical, cellular, and molecular GWAS. We have created an interactive web database (interactive Cross-Phenotype Analysis of GWAS database (iCPAGdb); http://cpag.oit.duke.edu) to facilitate exploration and allow rapid analysis of user-uploaded GWAS summary statistics. This database revealed well-known relationships among phenotypes, as well as the generation of novel hypotheses to explain the pathophysiology of common diseases. Application of iCPAGdb to a recent GWAS of severe COVID-19 demonstrated unexpected overlap of GWAS signals between COVID-19 and human diseases, including with idiopathic pulmonary fibrosis driven by the DPP9 locus. Transcriptomics from peripheral blood of COVID-19 patients demonstrated that DPP9 was induced in SARS-CoV-2 compared to healthy controls or those with bacterial infection. Further investigation of cross-phenotype SNPs with severe COVID-19 demonstrated colocalization of the GWAS signal of the ABO locus with plasma protein levels of a reported receptor of SARS-CoV-2, CD209 (DC-SIGN), pointing to a possible mechanism whereby glycosylation of CD209 by ABO may regulate COVID-19 disease severity. Thus, connecting genetically related traits across phenotypic scales links human diseases to molecular and cellular measurements that can reveal mechanisms and lead to novel biomarkers and therapeutic approaches.

Dysregulated lipid and fatty acid metabolism link perfluoroalkyl substances exposure and impaired glucose metabolism in young adults

BACKGROUND

Per- and polyfluoroalkyl substances (PFASs) exposure is ubiquitous among the US population and has been linked to adverse health outcomes including cardiometabolic diseases, immune dysregulation and endocrine disruption. However, the metabolic mechanism underlying the adverse health effect of PFASs exposure is unknown.

OBJECTIVE

The aim of this project is to investigate the association between PFASs exposure and altered metabolic pathways linked to increased cardiometabolic risk in young adults.

METHODS

A total of 102 young adults with 82% overweight or obese participants were enrolled from Southern California between 2014 and 2017. Cardiometabolic outcomes were assessed including oral glucose tolerance test (OGTT) measures, body fat and lipid profiles. High-resolution metabolomics was used to quantify plasma exposure levels of three PFAS congeners and intensity profiles of the untargeted metabolome. Fasting concentrations of 45 targeted metabolites involved in fatty acid and lipid metabolism were used to verify untargeted metabolomics findings. Bayesian Kernel Machine Regression (BKMR) was used to examine the associations between PFAS exposure mixture and cardiometabolic outcomes adjusting for covariates. Mummichog pathway enrichment analysis was used to explore PFAS-associated metabolic pathways. Moreover, the effect of PFAS exposure on the metabolic network, including metabolomic profiles and cardiometabolic outcomes, was investigated.

RESULTS

Higher exposure to perfluorooctanoic acid (PFOA) was associated with higher 30-minute glucose levels and glucose area under the curve (AUC) during the OGTT (p < 0.001). PFAS exposure was also associated with altered lipid pathways, which contributed to the metabolic network connecting PFOA and higher glucose levels following the OGTT. Targeted metabolomics analysis indicated that higher PFOA exposure was associated with higher levels of glycerol (p = 0.006), which itself was associated with higher 30-minute glucose (p = 0.006).

CONCLUSIONS

Increased lipolysis and fatty acid oxidation could contribute to the biological mechanisms linking PFAS exposure and impaired glucose metabolism among young adults. Findings of this study warrants future experimental studies and epidemiological studies with larger sample size to replicate.

Evaluating the precision of EBF1 SNP x stress interaction association: sex, race, and age differences in a big harmonized data set of 28,026 participants

In prior work, we identified a novel gene-by-stress association of EBF1's common variation (SNP rs4704963) with obesity (i.e., hip, waist) in Whites, which was further strengthened through multiple replications using our synthetic stress measure. We now extend this prior work in a precision medicine framework to find the risk group using harmonized data from 28,026 participants by evaluating the following: (a) EBF1 SNPxSTRESS interaction in Blacks; (b) 3-way interaction of EBF1 SNPxSTRESS with sex, race, and age; and (c) a race and sex-specific path linking EBF1 and stress to obesity to fasting glucose to the development of cardiometabolic disease risk. Our findings provided additional confirmation that genetic variation in EBF1 may contribute to stress-induced human obesity, including in Blacks (P = 0.022) that mainly resulted from race-specific stress due to "racism/discrimination" (P = 0.036) and "not meeting basic needs" (P = 0.053). The EBF1 gene-by-stress interaction differed significantly (P = 1.01e-03) depending on the sex of participants in Whites. Race and age also showed tentative associations (Ps = 0.103, 0.093, respectively) with this interaction. There was a significant and substantially larger path linking EBF1 and stress to obesity to fasting glucose to type 2 diabetes for the EBF1 minor allele group (coefficient = 0.28, P = 0.009, 95% CI = 0.07-0.49) compared with the same path for the EBF1 major allele homozygotes in White females and also a similar pattern of the path in Black females. Underscoring the race-specific key life-stress indicators (e.g., racism/discrimination) and also the utility of our synthetic stress, we identified the potential risk group of EBF1 and stress-induced human obesity and cardiometabolic disease.

Genetic Colorectal Cancer and Adenoma Risk Variants Are Associated with Increasing Cumulative Adenoma Counts

BACKGROUND

The genetic basis for most individuals with high cumulative lifetime colonic adenomas is unknown. We investigated associations between known colorectal cancer-risk single-nucleotide polymorphisms (SNP) and increasing cumulative adenoma counts.

METHODS

The Cooperative Studies Program #380 screening colonoscopy cohort includes 612 selected participants age 50 to 75 with genotyped blood samples and 10 years of clinical follow-up. We evaluated 41 published "colorectal cancer-risk SNPs" for associations with individual cumulative adenoma counts or having ≥10 cumulative adenomas. SNPs were analyzed singly or combined in a polygenic risk score (PRS). The PRS was constructed from eight published SNPs associated with multiple adenomas, termed "adenoma-risk SNPs."

RESULTS

Four colorectal cancer-risk SNPs were associated with increasing cumulative adenoma counts (P < 0.05): rs12241008 (gene: VTI1A), rs2423279 (BMP2/HAO1), rs3184504 (SH2B3), and rs961253 (FERMT1/BMP2), with risk allele risk ratios of 1.31, 1.29, 1.24, and 1.23, respectively. Three colorectal cancer-risk SNPs were associated with ≥10 cumulative adenomas (P < 0.05), with risk allele odds ratios of 2.09 (rs3184504), 2.30 (rs961253), and 1.94 (rs3217901). A weighted PRS comprised of adenoma-risk SNPs was associated with higher cumulative adenomas (weighted rate ratio = 1.57; P = 0.03).

CONCLUSIONS

In this mostly male veteran colorectal cancer screening cohort, several known colorectal cancer-risk SNPs were associated with increasing cumulative adenoma counts and the finding of ≥10 cumulative adenomas. In addition, an increasing burden of adenoma-risk SNPs, measured by a weighted PRS, was associated with higher cumulative adenomas.

IMPACT

Future work will seek to validate these findings in different populations and then augment current colorectal cancer risk prediction tools with precancerous, adenoma genetic data.

Baseline Colonoscopy Findings Associated With 10-Year Outcomes in a Screening Cohort Undergoing Colonoscopy Surveillance

BACKGROUND & AIMS

Few studies have evaluated long-term outcomes of ongoing colonoscopic screening and surveillance in a screening population. We aimed to determine the 10-year risk for advanced neoplasia (defined as adenomas ≥10mm, adenomas with villous histology or high-grade dysplasia, or colorectal cancer [CRC]) and assessed whether baseline colonoscopy findings were associated with long-term outcomes.

METHODS

We collected data from the Department of Veterans Affairs Cooperative Studies Program Study on 3121 asymptomatic veterans (50-75 years old) who underwent a screening colonoscopy from 1994 through 1997 at 13 medical centers and were then followed for 10 years or until death. We included 1915 subjects with at least 1 surveillance colonoscopy and estimated cumulative incidence of advanced neoplasia by Kaplan-Meier curves. We then fit a longitudinal joint model to estimate risk of advanced neoplasia at each subsequent examination after baseline, adjusting for multiple colonoscopies within individuals.

RESULTS

Through 10 years of follow-up, there were 146 individuals among all baseline colonoscopy groups found to have at least 1 incident advanced neoplasia. The cumulative 10-year incidence of advanced neoplasia was highest among those with baseline CRC (43.7%; 95% CI 13.0%-74.4%), followed by those with baseline advanced adenoma (AA) (21.9%; 95% CI 15.7-28.1). The cumulative 10-year incidence of advanced neoplasia was 6.3% (95% CI 4.1%-8.5%) and 4.1% (95% CI 2.7%-5.4%) for baseline 1 to 2 small adenomas (<1cm, and without villous histology or high-grade dysplasia) and no neoplasia, respectively (log-rank P = .10). After adjusting for prior surveillance, the risk of advanced neoplasia at each subsequent examination was not significantly increased in veterans with 1 or 2 small adenomas at baseline (odds ratio 0.96; 95% CI 0.67-1.41) compared with veterans with no baseline neoplasia.

CONCLUSIONS

Baseline screening colonoscopy findings associate with advanced neoplasia within 10 years. Individuals with only 1 or 2 small adenomas at baseline have a low risk of advanced neoplasia over 10 years. Alternative surveillance strategies, could be considered for these individuals.

Advanced neoplasia in Veterans at screening colonoscopy using the National Cancer Institute Risk Assessment Tool

BACKGROUND

Adapting screening strategy to colorectal cancer (CRC) risk may improve efficiency for all stakeholders however limited tools for such risk stratification exist. Colorectal cancers usually evolve from advanced neoplasms that are present for years. We applied the National Cancer Institute (NCI) CRC Risk Assessment Tool, which calculates future risk of CRC, to determine whether it could be used to predict current advanced neoplasia (AN) in a veteran cohort undergoing a baseline screening colonoscopy.

METHODS

This was a prospective assessment of the relationship between future CRC risk predicted by the NCI tool, and the presence of AN at screening colonoscopy. Family, medical, dietary and physical activity histories were collected at the time of screening colonoscopy and used to calculate absolute CRC risk at 5, 10 and 20 years. Discriminatory accuracy was assessed.

RESULTS

Of 3121 veterans undergoing screening colonoscopy, 94% had complete data available to calculate risk (N = 2934, median age 63 years, 100% men, and 15% minorities). Prevalence of AN at baseline screening colonoscopy was 11 % (N = 313). For tertiles of estimated absolute CRC risk at 5 years, AN prevalences were 6.54% (95% CI, 4.99, 8.09), 11.26% (95% CI, 9.28-13.24), and 14.21% (95% CI, 12.02-16.40). For tertiles of estimated risk at 10 years, the prevalences were 6.34% (95% CI, 4.81-7.87), 11.25% (95% CI, 9.27-13.23), and 14.42% (95% CI, 12.22-16.62). For tertiles of estimated absolute CRC risk at 20 years, current AN prevalences were 7.54% (95% CI, 5.75-9.33), 10.53% (95% CI, 8.45-12.61), and 12.44% (95% CI, 10.2-14.68). The area under the curve for predicting current AN was 0.60 (95% CI; 0.57-0.63, p < 0.0001) at 5 years, 0.60 (95% CI, 0.57-0.63, p < 0.0001) at 10 years and 0.58 (95% CI, 0.54-0.61, p < 0.0001) at 20 years.

CONCLUSION

The NCI tool had modest discriminatory function for estimating the presence of current advanced neoplasia in veterans undergoing a first screening colonoscopy. These findings are comparable to other clinically utilized cancer risk prediction models and may be used to inform the benefit-risk assessment of screening, particularly for patients with competing comorbidities and lower risk, for whom a non-invasive screening approach is preferred.

Systolic Blood Pressure and Socioeconomic Status in a large multi-study population

The present study used harmonized data from eight studies (N = 28,891) to examine the association between socioeconomic status (SES) and resting systolic blood pressure (SBP). The study replicates and extends our prior work on this topic by examining potential moderation of this association by race and gender. We also examined the extent to which body mass index (BMI), waist circumference (WC), and smoking might explain the association between SES and SBP. Data were available from six race/gender groups: 9200 Black women; 2337 Black men; 7248 White women; 6519 White men; 2950 Hispanic women; and 637 Hispanic men. Multivariable regression models showed that greater annual household income was associated with lower SBP in all groups except Hispanic men. The magnitude and form of this negative association differed across groups, with White women showing the strongest linear negative association. Among Black men and Hispanic women, the association was curvilinear: relatively flat among lower income levels, but then negative among higher income ranges. Education also was independently, negatively related to SBP, though evidence was weaker for race and gender differences in the strength of the association. Higher BMI and WC were associated with higher SBP, and current smoking with lower SBP. Inclusion of these risk factors resulted in only a modest change in the magnitude of the SBP and SES relation, accounting on average about 0.4 mmHg of the effect of income and 0.2 mmHg of the effect of education—effects unlikely to be clinically significant. Further understanding of mechanisms underlying the association between SBP and SES may improve risk stratification in clinical settings and potentially inform interventions aimed at reductions in social disparities in health.

•

Harmonized data (n = 28,291) examined association between blood pressure and socioeconomic status.

•

Replication and extension of prior work by examining potential moderation of this association by race and gender.

•

Further understanding of this association inform interventions aimed at reductions in social disparities in health.

An electronic family health history tool to identify and manage patients at increased risk for colorectal cancer: protocol for a randomized controlled trial

Background

Colorectal cancer is the fourth most commonly diagnosed cancer in the United States. Approximately 3–10% of the population has an increased risk for colorectal cancer due to family history and warrants more frequent or intensive screening. Yet, < 50% of that high-risk population receives guideline-concordant care. Systematic collection of family health history and decision support may improve guideline-concordant screening for patients at increased risk of colorectal cancer. We seek to test the effectiveness of a web-based, systematic family health history collection tool and decision support platform (MeTree) to improve risk assessment and appropriate management of colorectal cancer risk among patients in the Department of Veterans Affairs primary care practices.

Methods

In this ongoing randomized controlled trial, primary care providers at the Durham Veterans Affairs Health Care System and the Madison VA Medical Center are randomized to immediate intervention or wait-list control. Veterans are eligible if assigned to enrolled providers, have an upcoming primary care appointment, and have no conditions that would place them at increased risk for colorectal cancer (such as personal history, adenomatous polyps, or inflammatory bowel disease). Those with a recent lower endoscopy (e.g. colonoscopy, sigmoidoscopy) are excluded. Immediate intervention patients put their family health history information into a web-based platform, MeTree, which provides both patient- and provider-facing decision support reports. Wait-list control patients access MeTree 12 months post-consent. The primary outcome is the risk-concordant colorectal cancer screening referral rate obtained via chart review. Secondary outcomes include patient completion of risk management recommendations (e.g. colonoscopy) and referral for genetic consultation. We will also conduct an economic analysis and an assessment of providers’ experience with MeTree clinical decision support recommendations to inform future implementation efforts if the intervention is found to be effective.

Discussion

This trial will assess the feasibility and effectiveness of patient-collected family health history linked to decision support to promote risk-appropriate screening in a large healthcare system such as the Department of Veterans Affairs.

Trial registration

ClinicalTrials.gov, NCT02247336. Registered on 25 September 2014.

Electronic supplementary material

The online version of this article (10.1186/s13063-019-3659-y) contains supplementary material, which is available to authorized users.

Near-roadway air pollution exposure and altered fatty acid oxidation among adolescents and young adults - The interplay with obesity

BACKGROUND

Air pollution exposure has been shown to increase the risk of obesity and metabolic dysfunction in animal models and human studies. However, the metabolic pathways altered by air pollution exposure are unclear, especially in adolescents and young adults who are at a critical period in the development of cardio-metabolic diseases.

OBJECTIVES

The aim of this study was to examine the associations between air pollution exposure and indices of fatty acid and amino acid metabolism.

METHODS

A total of 173 young adults (18-23 years) from eight Children's Health Study (CHS) Southern California communities were examined from 2014 to 2018. Near-roadway air pollution (NRAP) exposure (freeway and non-freeway) and regional air pollution exposure (nitrogen dioxide, ozone and particulate matter) during one year before the study visit were estimated based on participants' residential addresses. Serum concentrations of 64 targeted metabolites including amino acids, acylcarnitines, non-esterified fatty acid (NEFA) and glycerol were measured in fasting serum samples. Principal component analysis of metabolites was performed to identify metabolite clusters that represent key metabolic pathways. Mixed effects models were used to analyze the associations of air pollution exposure with metabolomic principal component (PC) scores and individual metabolite concentrations adjusting for potential confounders.

RESULTS

Higher lagged one-year averaged non-freeway NRAP exposure was associated with higher concentrations of NEFA oxidation byproducts and higher NEFA-related PC score (all p's ≤ 0.038). The effect sizes were larger among obese individuals (interaction p = 0.047). Among females, higher freeway NRAP exposure was also associated with a higher NEFA-related PC score (p = 0.042). Among all participants, higher freeway NRAP exposure was associated with a lower PC score for lower concentrations of short- and median-chain acylcarnitines (p = 0.044).

CONCLUSIONS

Results of this study indicate that NRAP exposure is associated with altered fatty acid metabolism, which could contribute to the metabolic perturbation in obese youth.

Characterization of temporal relationships of comorbidities developed following cancer diagnoses in veterans

e18049

Background: Understanding patient trajectories and common sequences of comorbidity accrual among those newly diagnosed with cancer is critical for precision approaches to care and prevention. The Veterans Affairs (VA) Cooperative Studies Program (CSP) #380 cohort includes 3,121 healthy asymptomatic veterans who underwent screening colonoscopy and were followed for at least ten years. The current analysis leverages computational approaches to characterize the temporal relationships of diagnoses in CSP #380 participants following diagnosis of colorectal or other cancers. Methods: Patients enrolled in CSP #380 with at least 5 years of linked electronic health record data from the VA Corporate Data Warehouse (October 1999-December 2015) were included. Cancer diagnoses and their most common subsequent new diagnoses were identified per patient by the first instance of each three-digit ICD-9 diagnosis affecting at least 50 patients. Pairwise chronological relative risks (RR) between subsequent diagnoses were represented as a directed network graph, which maps the probability of developing a diagnosis following a prior diagnosis. Results: A total of 2,210 patients were included. The most common cancer diagnoses were prostate (436), thoracic (169), bladder (120), colon (72), and kidney cancers (65). Most first diagnoses following a cancer diagnosis were related to progressive cancer or acute/subacute treatment toxicity. For prostate cancer, comorbidities with greatest RR were carcinoma in situ (RR 6.85), unspecified (NOS) metastases (2.75), and urethral stricture (2.53). For lung cancer, they were metastases of respiratory and digestive sites (12.24), lymph nodes (6.47), and NOS (5.68), pneumothorax and air leak (4.16), and convalescence and palliative care (3.07). In bladder cancer, they were carcinoma in situ (9.00), cystitis (6.78), kidney or other urinary cancer (6.19), attention to artificial openings (3.40), and urethral stricture (2.78). These and other results were visualized with network graphs. Conclusions: Computational techniques can identify and visualize future health concerns following cancer diagnoses. In this cohort of initially healthy and asymptomatic veterans on a prospective screening colonoscopy study, most subsequent diagnoses were related to cancer or toxicities of therapy, as might be expected in an aging cohort. Future work may focus on streamlining in-clinic identification of potential high likelihood comorbidities.

Evaluating DNA methylation age on the Illumina MethylationEPIC Bead Chip

DNA methylation age (DNAm age) has become a widely utilized epigenetic biomarker for the aging process. The Horvath method for determining DNAm age is perhaps the most widely utilized and validated DNA methylation age assessment measure. Horvath DNAm age is calculated based on methylation measurements at 353 loci, present on Illumina's 450k and 27k DNA methylation microarrays. With increasing use of the more recently developed Illumina MethylationEPIC (850k) microarray, it is worth revisiting this aging measure to evaluate estimation differences due to array design. Of the requisite 353 loci, 17 are missing from the 850k microarray. Similarly, an alternate, 71 loci DNA methylation age assessment measure created by Hannum et al. is missing 6 requisite loci. Using 17 datasets with 27k, 450k, and/or 850k methylation data, we compared each sample's epigenetic age estimated from all 353 loci required by the Horvath DNAm age calculator, and using only the 336 loci available on the 850k array. In 450k/27k data, removing loci not on the 850k array resulted in underestimation of Horvath's DNAm age. Underestimation of Horvath DNAm age increased from ages 0 to ~20, remaining stable thereafter (mean deviation = -3.46 y, SD = 1.13 for individuals ≥20 years). Underestimation of Horvath's DNAm age by the reduced 450k/27k data was similar to the underestimation observed in the 850k data indicating it is driven by missing probes. In analogous examination of Hannum's DNAm age, the magnitude and direction of epigenetic age misestimation varied with chronological age. In conclusion, inter-array deviations in DNAm age estimations may be largely driven by missing probes between arrays, despite default probe imputation procedures. Though correlations and associations based on Horvath's DNAm age may be unaffected, researchers should exercise caution when interpreting results based on absolute differences in DNAm age or when mixing samples assayed on different arrays.