Proteomic analysis of cardiorespiratory fitness for prediction of mortality and multisystem disease risks.

Despite the wide effects of cardiorespiratory fitness (CRF) on metabolic, cardiovascular, pulmonary and neurological health, challenges in the feasibility and reproducibility of CRF measurements have impeded its use for clinical decision-making. Here we link proteomic profiles to CRF in 14,145 individuals across four international cohorts with diverse CRF ascertainment methods to establish, validate and characterize a proteomic CRF score. In a cohort of around 22,000 individuals in the UK Biobank, a proteomic CRF score was associated with a reduced risk of all-cause mortality (unadjusted hazard ratio 0.50 (95% confidence interval 0.48-0.52) per 1 s.d. increase). The proteomic CRF score was also associated with multisystem disease risk and provided risk reclassification and discrimination beyond clinical risk factors, as well as modulating high polygenic risk of certain diseases. Finally, we observed dynamicity of the proteomic CRF score in individuals who undertook a 20-week exercise training program and an association of the score with the degree of the effect of training on CRF, suggesting potential use of the score for personalization of exercise recommendations. These results indicate that population-based proteomics provides biologically relevant molecular readouts of CRF that are additive to genetic risk, potentially modifiable and clinically translatable.

Overcome the Limitation of Phenome-Wide Association Studies (PheWAS): Extension of PheWAS to Efficient and Robust Large-Scale ICD Codes Analysis.

The Phenome-wide association studies (PheWAS) have become widely used for efficient, high-throughput evaluation of relationship between a genetic factor and a large number of disease phenotypes, typically extracted from a DNA biobank linked with electronic medical records (EMR). Phecodes, billing code-derived disease case-control status, are usually used as outcome variables in PheWAS and logistic regression has been the standard choice of analysis method. Since the clinical diagnoses in EMR are often inaccurate with errors which can lead to biases in the odds ratio estimates, much effort has been put to accurately define the cases and controls to ensure an accurate analysis. Specifically in order to correctly classify controls in the population, an exclusion criteria list for each Phecode was manually compiled to obtain unbiased odds ratios. However, the accuracy of the list cannot be guaranteed without extensive data curation process. The costly curation process limits the efficiency of large-scale analyses that take full advantage of all structured phenotypic information available in EMR. Here, we proposed to estimate relative risks (RR) instead. We first demonstrated the desired nature of RR that overcomes the inaccuracy in the controls via theoretical formula. With simulation and real data application, we further confirmed that RR is unbiased without compiling exclusion criteria lists. With RR as estimates, we are able to efficiently extend PheWAS to a larger-scale, phenome construction agnostic analysis of phenotypes, using ICD 9/10 codes, which preserve much more disease-related clinical information than Phecodes.

Interoperability of phenome-wide multimorbidity patterns: a comparative study of two large-scale EHR systems.

Background: Electronic health records (EHR) are increasingly used for studying multimorbidities. However, concerns about accuracy, completeness, and EHRs being primarily designed for billing and administrative purposes raise questions about the consistency and reproducibility of EHR-based multimorbidity research. Methods: Utilizing phecodes to represent the disease phenome, we analyzed pairwise comorbidity strengths using a dual logistic regression approach and constructed multimorbidity as an undirected weighted graph. We assessed the consistency of the multimorbidity networks within and between two major EHR systems at local (nodes and edges), meso (neighboring patterns), and global (network statistics) scales. We present case studies to identify disease clusters and uncover clinically interpretable disease relationships. We provide an interactive web tool and a knowledge base combining data from multiple sources for online multimorbidity analysis. Findings: Analyzing data from 500,000 patients across Vanderbilt University Medical Center and Mass General Brigham health systems, we observed a strong correlation in disease frequencies (Kendall's τ = 0.643) and comorbidity strengths (Pearson ρ = 0.79). Consistent network statistics across EHRs suggest similar structures of multimorbidity networks at various scales. Comorbidity strengths and similarities of multimorbidity connection patterns align with the disease genetic correlations. Graph-theoretic analyses revealed a consistent core-periphery structure, implying efficient network clustering through threshold graph construction. Using hydronephrosis as a case study, we demonstrated the network's ability to uncover clinically relevant disease relationships and provide novel insights. Interpretation: Our findings demonstrate the robustness of large-scale EHR data for studying phenome-wide multimorbidities. The alignment of multimorbidity patterns with genetic data suggests the potential utility for uncovering shared biology of diseases. The consistent core-periphery structure offers analytical insights to discover complex disease interactions. This work also sets the stage for advanced disease modeling, with implications for precision medicine. Funding: VUMC Biostatistics Development Award, the National Institutes of Health, and the VA CSRD.

Impact of Insurance Status and Region on Angiotensin Receptor-Neprilysin Inhibitor Prescription During Heart Failure Hospitalizations.

BACKGROUND: An angiotensin receptor-neprilysin inhibitor (ARNI) is the preferred renin-angiotensin system (RAS) inhibitor for heart failure with reduced ejection fraction (HFrEF). Among eligible patients, insurance status and prescriber concern regarding out-of-pocket costs may constrain early initiation of ARNI and other new therapies. OBJECTIVES: In this study, the authors sought to evaluate the association of insurance and other social determinants of health with ARNI initiation at discharge from HFrEF hospitalization. METHODS: The authors analyzed ARNI initiation from January 2017 to June 2020 among patients with HFrEF eligible to receive RAS inhibitor at discharge from hospitals in the Get With The Guidelines-Heart Failure registry. The primary outcome was the proportion of ARNI prescription at discharge among those prescribed RAS inhibitor who were not on ARNI on admission. A logistic regression model was used to determine the association of insurance status, U.S. region, and their interaction, as well as self-reported race, with ARNI initiation at discharge. RESULTS: From 42,766 admissions, 24,904 were excluded for absolute or relative contraindications to RAS inhibitors. RAS inhibitors were prescribed for 16,817 (94.2%) of remaining discharges, for which ARNI was prescribed in 1,640 (9.8%). Self-reported Black patients were less likely to be initiated on ARNI compared to self-reported White patients (OR: 0.64; 95% CI: 0.50-0.81). Compared to Medicare beneficiaries, patients with third-party insurance, Medicaid, or no insurance were less likely to be initiated on ARNI (OR: 0.47 [95% CI: 0.31-0.72], OR: 0.41 [95% CI: 0.25-0.67], and OR: 0.20 [95% CI: 0.08-0.47], respectively). ARNI therapy varied by hospital region, with lowest utilization in the Mountain region. An interaction was demonstrated between the impact of insurance disparities and hospital region. CONCLUSIONS: Among patients hospitalized between 2017 and 2020 for HFrEF who were prescribed RAS inhibitor therapy at discharge, insurance status, geographic region, and self-reported race were associated with ARNI initiation.

Development and Implementation of an Integrated Preclinical Atherosclerosis Database.

BACKGROUND: Basic scientists have used preclinical animal models to explore mechanisms driving human diseases for decades, resulting in thousands of publications, each supporting causative inferences. Despite substantial advances in the mechanistic construct of disease, there has been limited translation from individual studies to advances in clinical care. An integrated approach to these individual studies has the potential to improve translational success. METHODS: Using atherosclerosis as a test case, we extracted data from the 2 most common mouse models of atherosclerosis (ApoE [apolipoprotein E]-knockout and LDLR [low-density lipoprotein receptor]-knockout). We restricted analyses to manuscripts published in 2 well-established journals, Arteriosclerosis, Thrombosis, and Vascular Biology and Circulation, as of query in 2021. Predefined variables including experimental conditions, intervention, and outcomes were extracted from each publication to produce a preclinical atherosclerosis database. RESULTS: Extracted data include animal sex, diet, intervention type, and distinct plaque pathologies (size, inflammation, and lipid content). Procedures are provided to standardize data extraction, attribute interventions to specific genes, and transform the database for use with available transcriptomics software. The database integrates hundreds of genes, each directly tested in vivo for causation in a murine atherosclerosis model. The database is provided to allow the research community to perform integrated analyses that reflect the global impact of decades of atherosclerosis investigation. CONCLUSIONS: This database is provided as a resource for future interrogation of sub-data sets associated with distinct plaque pathologies, cell type, or sex. We also provide the methods and software needed to expand this data set and apply this approach to the extensive repository of peer-reviewed data utilizing preclinical models to interrogate mechanisms of diverse human diseases.

Feasibility of precision smoking treatment in a low-income community setting: results of a pilot randomized controlled trial in The Southern Community Cohort Study.

BACKGROUND: The feasibility of precision smoking treatment in socioeconomically disadvantaged communities has not been studied. METHODS: Participants in the Southern Community Cohort Study who smoked daily were invited to join a pilot randomized controlled trial of three smoking cessation interventions: guideline-based care (GBC), GBC plus nicotine metabolism-informed care (MIC), and GBC plus counseling guided by a polygenic risk score (PRS) for lung cancer. Feasibility was assessed by rates of study enrollment, engagement, and retention, targeting > 70% for each. Using logistic regression, we also assessed whether feasibility varied by age, sex, race, income, education, and attitudes toward precision smoking treatment. RESULTS: Of 92 eligible individuals (79.3% Black; 68.2% with household income < $15,000), 67 (72.8%; 95% CI 63.0-80.9%) enrolled and were randomized. Of these, 58 (86.6%; 95% CI 76.4-92.8%) engaged with the intervention, and of these engaged participants, 43 (74.1%; 95% CI 61.6-83.7%) were retained at 6-month follow-up. Conditional on enrollment, older age was associated with lower engagement (OR 0.83, 95% CI 0.73-0.95, p = 0.008). Conditional on engagement, retention was significantly lower in the PRS arm than in the GBC arm (OR 0.18, 95% CI 0.03-1.00, p = 0.050). No other selection effects were observed. CONCLUSIONS: Genetically informed precision smoking cessation interventions are feasible in socioeconomically disadvantaged communities, exhibiting high enrollment, engagement, and retention irrespective of race, sex, income, education, or attitudes toward precision smoking treatment. Future smoking cessation interventions in this population should take steps to engage older people and to sustain participation in interventions that include genetic risk counseling. TRIAL REGISTRATION: ClinicalTrials.gov No. NCT03521141, Registered 27 April 2018, https://www. CLINICALTRIALS: gov/study/NCT03521141.

Recurrence After Atrial Fibrillation Ablation and Investigational Biomarkers of Cardiac Remodeling.

BACKGROUND: Recurrence after atrial fibrillation (AF) ablation remains common. We evaluated the association between recurrence and levels of biomarkers of cardiac remodeling, and their ability to improve recurrence prediction when added to a clinical prediction model. METHODS AND RESULTS: Blood samples collected before de novo catheter ablation were analyzed. Levels of bone morphogenetic protein-10, angiopoietin-2, fibroblast growth factor-23, insulin-like growth factor-binding protein-7, myosin-binding protein C3, growth differentiation factor-15, interleukin-6, N-terminal pro-brain natriuretic peptide, and high-sensitivity troponin T were measured. Recurrence was defined as ≥30 seconds of an atrial arrhythmia 3 to 12 months postablation. Multivariable logistic regression was performed using biomarker levels along with clinical covariates: APPLE score (Age >65 years, Persistent AF, imPaired eGFR [<60 ml/min/1.73m2], LA diameter ≥43 mm, EF <50%; which includes age, left atrial diameter, left ventricular ejection fraction, persistent atrial fibrillation, and estimated glomerular filtration rate), preablation rhythm, sex, height, body mass index, presence of an implanted continuous monitor, year of ablation, and additional linear ablation. A total of 1873 participants were included. A multivariable logistic regression showed an association between recurrence and levels of angiopoietin-2 (odds ratio, 1.08 [95% CI, 1.02-1.15], P=0.007) and interleukin-6 (odds ratio, 1.02 [95% CI, 1.003-1.03]; P=0.02). The area under the receiver operating characteristic curve of a model that only contained clinical predictors was 0.711. The addition of any of the 9 studied biomarkers to the predictive model did not result in a statistically significant improvement in the area under the receiver operating characteristic curve. CONCLUSIONS: Higher angiopoietin-2 and interleukin-6 levels were associated with recurrence after atrial fibrillation ablation in multivariable modeling. However, the addition of biomarkers to a clinical prediction model did not significantly improve recurrence prediction.

Volumetric brain MRI signatures of heart failure with preserved ejection fraction in the setting of dementia.

Heart failure with preserved ejection fraction (HFpEF) is an important, emerging risk factor for dementia, but it is not clear whether HFpEF contributes to a specific pattern of neuroanatomical changes in dementia. A major challenge to studying this is the relative paucity of datasets of patients with dementia, with/without HFpEF, and relevant neuroimaging. We sought to demonstrate the feasibility of using modern data mining tools to create and analyze clinical imaging datasets and identify the neuroanatomical signature of HFpEF-associated dementia. We leveraged the bioinformatics tools at Vanderbilt University Medical Center to identify patients with a diagnosis of dementia with and without comorbid HFpEF using the electronic health record. We identified high resolution, clinically-acquired neuroimaging data on 30 dementia patients with HFpEF (age 76.9 ± 8.12 years, 61% female) as well as 301 age- and sex-matched patients with dementia but without HFpEF to serve as comparators (age 76.2 ± 8.52 years, 60% female). We used automated image processing pipelines to parcellate the brain into 132 structures and quantify their volume. We found six regions with significant atrophy associated with HFpEF: accumbens area, amygdala, posterior insula, anterior orbital gyrus, angular gyrus, and cerebellar white matter. There were no regions with atrophy inversely associated with HFpEF. Patients with dementia and HFpEF have a distinct neuroimaging signature compared to patients with dementia only. Five of the six regions identified in are in the temporo-parietal region of the brain. Future studies should investigate mechanisms of injury associated with cerebrovascular disease leading to subsequent brain atrophy.

LNK/SH2B3 loss of function increases susceptibility to murine and human atrial fibrillation.

AIMS: The lymphocyte adaptor protein (LNK) is a negative regulator of cytokine and growth factor signalling. The rs3184504 variant in SH2B3 reduces LNK function and is linked to cardiovascular, inflammatory, and haematologic disorders, including stroke. In mice, deletion of Lnk causes inflammation and oxidative stress. We hypothesized that Lnk-/- mice are susceptible to atrial fibrillation (AF) and that rs3184504 is associated with AF and AF-related stroke in humans. During inflammation, reactive lipid dicarbonyls are the major components of oxidative injury, and we further hypothesized that these mediators are critical drivers of the AF substrate in Lnk-/- mice. METHODS AND RESULTS: Lnk-/- or wild-type (WT) mice were treated with vehicle or 2-hydroxybenzylamine (2-HOBA), a dicarbonyl scavenger, for 3 months. Compared with WT, Lnk-/- mice displayed increased AF duration that was prevented by 2-HOBA. In the Lnk-/- atria, action potentials were prolonged with reduced transient outward K+ current, increased late Na+ current, and reduced peak Na+ current, pro-arrhythmic effects that were inhibited by 2-HOBA. Mitochondrial dysfunction, especially for Complex I, was evident in Lnk-/- atria, while scavenging lipid dicarbonyls prevented this abnormality. Tumour necrosis factor-α (TNF-α) and interleukin-1 beta (IL-1ß) were elevated in Lnk-/- plasma and atrial tissue, respectively, both of which caused electrical and bioenergetic remodelling in vitro. Inhibition of soluble TNF-α prevented electrical remodelling and AF susceptibility, while IL-1ß inhibition improved mitochondrial respiration but had no effect on AF susceptibility. In a large database of genotyped patients, rs3184504 was associated with AF, as well as AF-related stroke. CONCLUSION: These findings identify a novel role for LNK in the pathophysiology of AF in both experimental mice and humans. Moreover, reactive lipid dicarbonyls are critical to the inflammatory AF substrate in Lnk-/- mice and mediate the pro-arrhythmic effects of pro-inflammatory cytokines, primarily through electrical remodelling.

Unbiased characterization of atrial fibrillation phenotypic architecture provides insight to genetic liability and clinically relevant outcomes.

Background: Atrial Fibrillation (AF) is a common and clinically heterogeneous arrythmia. Machine learning (ML) algorithms can define data-driven disease subtypes in an unbiased fashion, but whether the AF subgroups defined in this way align with underlying mechanisms, such as high polygenic liability to AF or inflammation, and associate with clinical outcomes is unclear. Methods: We identified individuals with AF in a large biobank linked to electronic health records (EHR) and genome-wide genotyping. The phenotypic architecture in the AF cohort was defined using principal component analysis of 35 expertly curated and uncorrelated clinical features. We applied an unsupervised co-clustering machine learning algorithm to the 35 features to identify distinct phenotypic AF clusters. The clinical inflammatory status of the clusters was defined using measured biomarkers (CRP, ESR, WBC, Neutrophil %, Platelet count, RDW) within 6 months of first AF mention in the EHR. Polygenic risk scores (PRS) for AF and cytokine levels were used to assess genetic liability of clusters to AF and inflammation, respectively. Clinical outcomes were collected from EHR up to the last medical contact. Results: The analysis included 23,271 subjects with AF, of which 6,023 had available genome-wide genotyping. The machine learning algorithm identified 3 phenotypic clusters that were distinguished by increasing prevalence of comorbidities, particularly renal dysfunction, and coronary artery disease. Polygenic liability to AF across clusters was highest in the low comorbidity cluster. Clinically measured inflammatory biomarkers were highest in the high comorbid cluster, while there was no difference between groups in genetically predicted levels of inflammatory biomarkers. Subgroup assignment was associated with multiple clinical outcomes including mortality, stroke, bleeding, and use of cardiac implantable electronic devices after AF diagnosis. Conclusion: Patient subgroups identified by unsupervised clustering were distinguished by comorbidity burden and associated with risk of clinically important outcomes. Polygenic liability to AF across clusters was greatest in the low comorbidity subgroup. Clinical inflammation, as reflected by measured biomarkers, was lowest in the subgroup with lowest comorbidities. However, there were no differences in genetically predicted levels of inflammatory biomarkers, suggesting associations between AF and inflammation is driven by acquired comorbidities rather than genetic predisposition.

Identification of Clinical Drivers of Left Atrial Enlargement Through Genomics of Left Atrial Size.

BACKGROUND: Greater left atrial size is associated with a higher incidence of cardiovascular disease and mortality, but the full spectrum of diagnoses associated with left atrial enlargement in sex-stratified clinical populations is not well known. Our study sought to identify genetic risk mechanisms affecting left atrial diameter (LAD) in a clinical cohort. METHODS: Using Vanderbilt deidentified electronic health record, we studied 6163 females and 5993 males of European ancestry who had at least 1 LAD measure and available genotyping. A sex-stratified polygenic score was constructed for LAD variation and tested for association against 1680 International Classification of Diseases code-based phenotypes. Two-sample univariable and multivariable Mendelian randomization approaches were used to assess etiologic relationships between candidate associations and LAD. RESULTS: A phenome-wide association study identified 25 International Classification of Diseases code-based diagnoses in females and 11 in males associated with a polygenic score of LAD (false discovery rate q<0.01), 5 of which were further evaluated by Mendelian randomization (waist circumference [WC], atrial fibrillation, heart failure, systolic blood pressure, and coronary artery disease). Sex-stratified differences in the genetic associations between risk factors and a polygenic score for LAD were observed (WC for females; heart failure, systolic blood pressure, atrial fibrillation, and WC for males). By multivariable Mendelian randomization, higher WC remained significantly associated with larger LAD in females, whereas coronary artery disease, WC, and atrial fibrillation remained significantly associated with larger LAD in males. CONCLUSIONS: In a clinical population, we identified, by genomic approaches, potential etiologic risk factors for larger LAD. Further studies are needed to confirm the extent to which these risk factors may be modified to prevent or reverse adverse left atrial remodeling and the extent to which sex modifies these risk factors.

Next-generation phenotyping: introducing phecodeX for enhanced discovery research in medical phenomics.

MOTIVATION: Phecodes are widely used and easily adapted phenotypes based on International Classification of Diseases codes. The current version of phecodes (v1.2) was designed primarily to study common/complex diseases diagnosed in adults; however, there are numerous limitations in the codes and their structure. RESULTS: Here, we present phecodeX, an expanded version of phecodes with a revised structure and 1,761 new codes. PhecodeX adds granularity to phenotypes in key disease domains that are under-represented in the current phecode structure-including infectious disease, pregnancy, congenital anomalies, and neonatology-and is a more robust representation of the medical phenome for global use in discovery research. AVAILABILITY AND IMPLEMENTATION: phecodeX is available at https://github.com/PheWAS/phecodeX.

Development and implementation of an integrated preclinical atherosclerosis database.

Background: Basic scientists have used preclinical animal models to explore mechanisms driving human diseases for decades, resulting in thousands of publications, each supporting causative inferences. Despite substantial advances in the mechanistic construct of disease, there has been limited translation from individual studies to advances in clinical care. An integrated approach to these individual studies has the potential to improve translational success. Methods: Using atherosclerosis as a test case, we extracted data from the two most common mouse models of atherosclerosis (ApoE and LDLR knockout). We restricted analyses to manuscripts published in two well-established journals, Arteriosclerosis, Thrombosis, and Vascular Biology and Circulation, as of query in 2021. Predefined variables including experimental conditions, intervention and outcomes were extracted from each publication to produce a preclinical atherosclerosis database. Results: Extracted data include animal sex, diet, intervention type and distinct plaque pathologies (size, inflammation, lipid content). Procedures are provided to standardize data extraction, attribute interventions to specific genes and transform the database for use with available transcriptomics software. The database integrates hundreds of genes, each directly tested in vivo for causation in a murine atherosclerosis model. The database is provided to allow the research community to perform integrated analyses that reflect the global impact of decades of atherosclerosis investigation. Conclusions: Future database uses include interrogation of sub-datasets associated with distinct plaque pathologies, cell-type or sex. We provide the methods and software needed to apply this approach to the extensive repository of peer-reviewed data utilizing preclinical models to interrogate mechanisms of diverse human diseases.

Clonal hematopoiesis of indeterminate potential and outcomes after heart transplantation: A multicenter study.

Cardiac allograft vasculopathy (CAV) is a leading cause of late graft failure and mortality after heart transplantation (HT). Sharing some features with atherosclerosis, CAV results in diffuse narrowing of the epicardial coronaries and microvasculature, with consequent graft ischemia. Recently, clonal hematopoiesis of indeterminate potential (CHIP) has emerged as a risk factor for cardiovascular disease and mortality. We aimed to investigate the relationship between CHIP and posttransplant outcomes, including CAV. We analyzed 479 HT recipients with stored DNA samples at 2 high-volume transplant centers, Vanderbilt University Medical Center and Columbia University Irving Medical Center. We explored the association between the presence of CHIP mutations with CAV and mortality after HT. In this case-control analysis, carriers of CHIP mutations were not at increased risk of CAV or mortality after HT. In a large multicenter genomics study of the heart transplant population, the presence of CHIP mutations was not associated with an increased risk of CAV or posttransplant mortality.

Association of HIV Infection and Incident Abdominal Aortic Aneurysm Among 143 001 Veterans.

BACKGROUND: People with HIV (PWH) have an increased risk of cardiovascular disease. Previous cross-sectional data suggest there is a higher prevalence of abdominal aortic aneurysm (AAA) in PWH than in those without HIV. Whether PWH have an increased risk of incident AAA compared with those without HIV is unknown. METHODS: We analyzed data among participants without prevalent AAA from the Veterans Aging Cohort Study, a prospective, observational, longitudinal cohort of veterans with HIV matched 1:2 with veterans without HIV infection. We calculated AAA rates by HIV status and assessed the association between HIV infection and incident AAA using Cox proportional hazards models. We defined AAA using the International Classification of Diseases, 9th or 10th revision, or Current Procedural Terminology codes and adjusted all models for demographic characteristics, cardiovascular disease risk factors, and substance use. Secondary analyses examined the association between time-varying CD4+ T-cell count or HIV viral load and incident AAA. RESULTS: Among 143 001 participants (43 766 with HIV), over a median follow-up of 8.7 years, there were 2431 incident AAA events (26.4% among PWH). Rates of incident AAA per 1000 person-years were similar among PWH (2.0 [95% CI, 1.9-2.2]) and people without HIV (2.2 [95% CI, 2.1-2.3]). There was no evidence that HIV infection increased the risk of incident AAA compared with no HIV infection (adjusted hazard ratio, 1.02 [95% CI, 0.92-1.13]). In adjusted analyses with time-varying CD4+ T-cell counts or HIV viral load, PWH with CD4+ T-cell counts <200 cells/mm3 (adjusted hazard ratio, 1.29 [95% CI, 1.02-1.65]) or HIV viral load ≥500 copies/mL (adjusted hazard ratio, 1.29 [95% CI, 1.09-1.52]) had an increased risk of AAA compared with those without HIV. CONCLUSIONS: HIV infection is associated with an increased risk of AAA among those with low CD4+ T-cell counts or elevated HIV viral load over time.