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.