A Metabolic Basis for Endothelial-to-Mesenchymal Transition.

Endothelial-to-mesenchymal transition (EndoMT) is a cellular process often initiated by the transforming growth factor ß (TGF-ß) family of ligands. Although required for normal heart valve development, deregulated EndoMT is linked to a wide range of pathological conditions. Here, we demonstrate that endothelial fatty acid oxidation (FAO) is a critical in vitro and in vivo regulator of EndoMT. We further show that this FAO-dependent metabolic regulation of EndoMT occurs through alterations in intracellular acetyl-CoA levels. Disruption of FAO via conditional deletion of endothelial carnitine palmitoyltransferase II (Cpt2E-KO) augments the magnitude of embryonic EndoMT, resulting in thickening of cardiac valves. Consistent with the known pathological effects of EndoMT, adult Cpt2E-KO mice demonstrate increased permeability in multiple vascular beds. Taken together, these results demonstrate that endothelial FAO is required to maintain endothelial cell fate and that therapeutic manipulation of endothelial metabolism could provide the basis for treating a growing number of EndoMT-linked pathological conditions.

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.

Multiancestry Genome-Wide Association Study of Aortic Stenosis Identifies Multiple Novel Loci in the Million Veteran Program.

BACKGROUND: Calcific aortic stenosis (CAS) is the most common valvular heart disease in older adults and has no effective preventive therapies. Genome-wide association studies (GWAS) can identify genes influencing disease and may help prioritize therapeutic targets for CAS. METHODS: We performed a GWAS and gene association study of 14 451 patients with CAS and 398 544 controls in the Million Veteran Program. Replication was performed in the Million Veteran Program, Penn Medicine Biobank, Mass General Brigham Biobank, BioVU, and BioMe, totaling 12 889 cases and 348 094 controls. Causal genes were prioritized from genome-wide significant variants using polygenic priority score gene localization, expression quantitative trait locus colocalization, and nearest gene methods. CAS genetic architecture was compared with that of atherosclerotic cardiovascular disease. Causal inference for cardiometabolic biomarkers in CAS was performed using Mendelian randomization and genome-wide significant loci were characterized further through phenome-wide association study. RESULTS: We identified 23 genome-wide significant lead variants in our GWAS representing 17 unique genomic regions. Of the 23 lead variants, 14 were significant in replication, representing 11 unique genomic regions. Five replicated genomic regions were previously known risk loci for CAS (PALMD, TEX41, IL6, LPA, FADS) and 6 were novel (CEP85L, FTO, SLMAP, CELSR2, MECOM, CDAN1). Two novel lead variants were associated in non-White individuals (P<0.05): rs12740374 (CELSR2) in Black and Hispanic individuals and rs1522387 (SLMAP) in Black individuals. Of the 14 replicated lead variants, only 2 (rs10455872 [LPA], rs12740374 [CELSR2]) were also significant in atherosclerotic cardiovascular disease GWAS. In Mendelian randomization, lipoprotein(a) and low-density lipoprotein cholesterol were both associated with CAS, but the association between low-density lipoprotein cholesterol and CAS was attenuated when adjusting for lipoprotein(a). Phenome-wide association study highlighted varying degrees of pleiotropy, including between CAS and obesity at the FTO locus. However, the FTO locus remained associated with CAS after adjusting for body mass index and maintained a significant independent effect on CAS in mediation analysis. CONCLUSIONS: We performed a multiancestry GWAS in CAS and identified 6 novel genomic regions in the disease. Secondary analyses highlighted the roles of lipid metabolism, inflammation, cellular senescence, and adiposity in the pathobiology of CAS and clarified the shared and differential genetic architectures of CAS with atherosclerotic cardiovascular diseases.

Transcriptional Dysregulation Underlies Both Monogenic Arrhythmia Syndrome and Common Modifiers of Cardiac Repolarization.

BACKGROUND: Brugada syndrome (BrS) is an inherited arrhythmia syndrome caused by loss-of-function variants in the cardiac sodium channel gene SCN5A (sodium voltage-gated channel alpha subunit 5) in ≈20% of subjects. We identified a family with 4 individuals diagnosed with BrS harboring the rare G145R missense variant in the cardiac transcription factor TBX5 (T-box transcription factor 5) and no SCN5A variant. METHODS: We generated induced pluripotent stem cells (iPSCs) from 2 members of a family carrying TBX5-G145R and diagnosed with Brugada syndrome. After differentiation to iPSC-derived cardiomyocytes (iPSC-CMs), electrophysiologic characteristics were assessed by voltage- and current-clamp experiments (n=9 to 21 cells per group) and transcriptional differences by RNA sequencing (n=3 samples per group), and compared with iPSC-CMs in which G145R was corrected by CRISPR/Cas9 approaches. The role of platelet-derived growth factor (PDGF)/phosphoinositide 3-kinase (PI3K) pathway was elucidated by small molecule perturbation. The rate-corrected QT (QTc) interval association with serum PDGF was tested in the Framingham Heart Study cohort (n=1893 individuals). RESULTS: TBX5-G145R reduced transcriptional activity and caused multiple electrophysiologic abnormalities, including decreased peak and enhanced "late" cardiac sodium current (INa), which were entirely corrected by editing G145R to wild-type. Transcriptional profiling and functional assays in genome-unedited and -edited iPSC-CMs showed direct SCN5A down-regulation caused decreased peak INa, and that reduced PDGF receptor (PDGFRA [platelet-derived growth factor receptor α]) expression and blunted signal transduction to PI3K was implicated in enhanced late INa. Tbx5 regulation of the PDGF axis increased arrhythmia risk due to disruption of PDGF signaling and was conserved in murine model systems. PDGF receptor blockade markedly prolonged normal iPSC-CM action potentials and plasma levels of PDGF in the Framingham Heart Study were inversely correlated with the QTc interval (P<0.001). CONCLUSIONS: These results not only establish decreased SCN5A transcription by the TBX5 variant as a cause of BrS, but also reveal a new general transcriptional mechanism of arrhythmogenesis of enhanced late sodium current caused by reduced PDGF receptor-mediated PI3K signaling.

Interactive network-based clustering and investigation of multimorbidity association matrices with associationSubgraphs.

MOTIVATION: Making sense of networked multivariate association patterns is vitally important to many areas of high-dimensional analysis. Unfortunately, as the data-space dimensions grow, the number of association pairs increases in O(n2); this means that traditional visualizations such as heatmaps quickly become too complicated to parse effectively. RESULTS: Here, we present associationSubgraphs: a new interactive visualization method to quickly and intuitively explore high-dimensional association datasets using network percolation and clustering. The goal is to provide an efficient investigation of association subgraphs, each containing a subset of variables with stronger and more frequent associations among themselves than the remaining variables outside the subset, by showing the entire clustering dynamics and providing subgraphs under all possible cutoff values at once. Particularly, we apply associationSubgraphs to a phenome-wide multimorbidity association matrix generated from an electronic health record and provide an online, interactive demonstration for exploring multimorbidity subgraphs. AVAILABILITY AND IMPLEMENTATION: An R package implementing both the algorithm and visualization components of associationSubgraphs is available at https://github.com/tbilab/associationsubgraphs. Online documentation is available at https://prod.tbilab.org/associationsubgraphs_info/. A demo using a multimorbidity association matrix is available at https://prod.tbilab.org/associationsubgraphs-example/.

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.

Impact of metabolic and weight components on incident asthma using a real-world cohort.

BACKGROUND: Obesity and metabolic dysregulation (MetD) have increasing prevalence and adversely affect asthma morbidity and therapeutic response. OBJECTIVE: To determine the role of weight and MetD on incident asthma in adulthood. METHODS: In a retrospective, longitudinal cohort of patients, we performed a time-to-asthma diagnosis analysis after a 3-year landmark period (t0-t3) during which weight and MetD components were evaluated. We assessed incident asthma risk with MetD components and weight. RESULTS: In total, 90,081 patients met the inclusion criteria, with 836 cases (0.93%) of incident asthma in our primary cohort. Diabetes present at t0, but no other MetD components, was associated with increased risk of asthma (adjusted hazard ratio = 1.85, 95% CI: 1.27-2.71, P = .0002). The effect of weight on asthma risk, independent of other MetD components, identified individuals with overweight or obesity as having a 10-year attributable risk of 15.4%. Metformin was prescribed more frequently, and hemoglobin A1c levels were lower in patients with diabetes in whom asthma did not develop (P < .0001). CONCLUSION: Weight and diabetes prevention and management represent modifiable risk factors for adult asthma development.

Dyslipidemia, inflammation, calcification, and adiposity in aortic stenosis: a genome-wide study.

AIMS: Although highly heritable, the genetic etiology of calcific aortic stenosis (AS) remains incompletely understood. The aim of this study was to discover novel genetic contributors to AS and to integrate functional, expression, and cross-phenotype data to identify mechanisms of AS. METHODS AND RESULTS: A genome-wide meta-analysis of 11.6 million variants in 10 cohorts involving 653 867 European ancestry participants (13 765 cases) was performed. Seventeen loci were associated with AS at P ≤ 5 × 10-8, of which 15 replicated in an independent cohort of 90 828 participants (7111 cases), including CELSR2-SORT1, NLRP6, and SMC2. A genetic risk score comprised of the index variants was associated with AS [odds ratio (OR) per standard deviation, 1.31; 95% confidence interval (CI), 1.26-1.35; P = 2.7 × 10-51] and aortic valve calcium (OR per standard deviation, 1.22; 95% CI, 1.08-1.37; P = 1.4 × 10-3), after adjustment for known risk factors. A phenome-wide association study indicated multiple associations with coronary artery disease, apolipoprotein B, and triglycerides. Mendelian randomization supported a causal role for apolipoprotein B-containing lipoprotein particles in AS (OR per g/L of apolipoprotein B, 3.85; 95% CI, 2.90-5.12; P = 2.1 × 10-20) and replicated previous findings of causality for lipoprotein(a) (OR per natural logarithm, 1.20; 95% CI, 1.17-1.23; P = 4.8 × 10-73) and body mass index (OR per kg/m2, 1.07; 95% CI, 1.05-1.9; P = 1.9 × 10-12). Colocalization analyses using the GTEx database identified a role for differential expression of the genes LPA, SORT1, ACTR2, NOTCH4, IL6R, and FADS. CONCLUSION: Dyslipidemia, inflammation, calcification, and adiposity play important roles in the etiology of AS, implicating novel treatments and prevention strategies.

Common Ancestry-Specific Ion Channel Variants Predispose to Drug-Induced Arrhythmias.

BACKGROUND: Multiple reports associate the cardiac sodium channel gene (SCN5A) variants S1103Y and R1193Q with type 3 congenital long QT syndrome and drug-induced long QT syndrome. These variants are too common in ancestral populations to be highly arrhythmogenic at baseline, however: S1103Y allele frequency is 8.1% in African Americans and R1193Q 6.1% in East Asians. R1193Q is known to increase late sodium current (INa-L) in cardiomyocytes derived from induced pluripotent stem cells but the role of these variants in modulating repolarization remains poorly understood. METHODS: We determined the effect of S1103Y on QT intervals among African-American participants in a large electronic health record. Using cardiomyocytes derived from induced pluripotent stem cells carrying naturally occurring or genome-edited variants, we studied action potential durations (APDs) at baseline and after challenge with the repolarizing potassium current (IKr) blocker dofetilide and INa-L and IKr at baseline. RESULTS: In 1479 African-American participants with no confounding medications or diagnoses of heart disease, QT intervals in S1103Y carriers was no different from that in noncarriers. Baseline APD was no different in cells expressing the Y allele (SY, YY cells) compared with isogenic cells with the reference allele (SS cells). However, INa-L was increased in SY and YY cells and the INa-L blocker GS967 shortened APD in SY/YY but not SS cells (P<0.001). IKr was increased almost 2-fold in SY/YY cells compared with SS cells (tail current: 0.66±0.1 versus 1.2±0.1 pA/pF; P<0.001). Dofetilide challenge prolonged APD at much lower concentrations in SY (4.1 nmol/L [interquartile range, 1.5-9.3]; n=11) and YY (4.2 nmol/L [1.7-5.0]; n=5) than in SS cells (249 nmol/L [22.3-2905]; n=14; P<0.001 and P<0.01, respectively) and elicited afterdepolarizations in 8/16 SY/YY cells but only in 1/14 SS cells. R1193Q cells similarly displayed no difference in baseline APD but increased IKr and increased dofetilide sensitivity. CONCLUSIONS: These common ancestry-specific variants do not affect baseline repolarization, despite generating increased INa-L. We propose that increased IKr serves to maintain normal repolarization but increases the risk of manifest QT prolongation with IKr block in variant carriers. Our findings emphasize the need for inclusion of diverse populations in the study of adverse drug reactions.

Arrhythmia Variant Associations and Reclassifications in the eMERGE-III Sequencing Study.

BACKGROUND: Sequencing Mendelian arrhythmia genes in individuals without an indication for arrhythmia genetic testing can identify carriers of pathogenic or likely pathogenic (P/LP) variants. However, the extent to which these variants are associated with clinically meaningful phenotypes before or after return of variant results is unclear. In addition, the majority of discovered variants are currently classified as variants of uncertain significance, limiting clinical actionability. METHODS: The eMERGE-III study (Electronic Medical Records and Genomics Phase III) is a multicenter prospective cohort that included 21 846 participants without previous indication for cardiac genetic testing. Participants were sequenced for 109 Mendelian disease genes, including 10 linked to arrhythmia syndromes. Variant carriers were assessed with electronic health record-derived phenotypes and follow-up clinical examination. Selected variants of uncertain significance (n=50) were characterized in vitro with automated electrophysiology experiments in HEK293 cells. RESULTS: As previously reported, 3.0% of participants had P/LP variants in the 109 genes. Herein, we report 120 participants (0.6%) with P/LP arrhythmia variants. Compared with noncarriers, arrhythmia P/LP carriers had a significantly higher burden of arrhythmia phenotypes in their electronic health records. Fifty-four participants had variant results returned. Nineteen of these 54 participants had inherited arrhythmia syndrome diagnoses (primarily long-QT syndrome), and 12 of these 19 diagnoses were made only after variant results were returned (0.05%). After in vitro functional evaluation of 50 variants of uncertain significance, we reclassified 11 variants: 3 to likely benign and 8 to P/LP. CONCLUSIONS: Genome sequencing in a large population without indication for arrhythmia genetic testing identified phenotype-positive carriers of variants in congenital arrhythmia syndrome disease genes. As the genomes of large numbers of people are sequenced, the disease risk from rare variants in arrhythmia genes can be assessed by integrating genomic screening, electronic health record phenotypes, and in vitro functional studies. REGISTRATION: URL: https://www. CLINICALTRIALS: gov; Unique identifier; NCT03394859.

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.

Expression of CD70 Modulates Nitric Oxide and Redox Status in Endothelial Cells.

BACKGROUND: Endothelial dysfunction is a critical component in the pathogenesis of cardiovascular diseases and is closely associated with nitric oxide (NO) levels and oxidative stress. Here, we report on novel findings linking endothelial expression of CD70 (also known as CD27 ligand) with alterations in NO and reactive oxygen species. METHODS: CD70 expression was genetically manipulated in human aortic and pulmonary artery endothelial cells. Intracellular NO and hydrogen peroxide (H2O2) were measured using genetically encoded biosensors, and cellular phenotypes were assessed. RESULTS: An unbiased phenome-wide association study demonstrated that polymorphisms in CD70 associate with vascular phenotypes. Endothelial cells treated with CD70-directed short-interfering RNA demonstrated impaired wound closure, decreased agonist-stimulated NO levels, and reduced eNOS (endothelial nitric oxide synthase) protein. These changes were accompanied by reduced NO bioactivity, increased 3-nitrotyrosine levels, and a decrease in the eNOS binding partner heat shock protein 90. Following treatment with the thioredoxin inhibitor auranofin or with agonist histamine, intracellular H2O2 levels increased up to 80% in the cytosol, plasmalemmal caveolae, and mitochondria. There was increased expression of NADPH oxidase 1 complex and gp91phox; expression of copper/zinc and manganese superoxide dismutases was also elevated. CD70 knockdown reduced levels of the H2O2 scavenger catalase; by contrast, glutathione peroxidase 1 expression and activity were increased. CD70 overexpression enhanced endothelial wound closure, increased NO levels, and attenuated the reduction in eNOS mRNA induced by TNFα. CONCLUSIONS: Taken together, these data establish CD70 as a novel regulatory protein in endothelial NO and reactive oxygen species homeostasis, with implications for human vascular disease.

Systems Approach to Integrating Preclinical Apolipoprotein E-Knockout Investigations Reveals Novel Etiologic Pathways and Master Atherosclerosis Network in Humans.

OBJECTIVE: Animal models of atherosclerosis are used extensively to interrogate molecular mechanisms in serial fashion. We tested whether a novel systems biology approach to integration of preclinical data identifies novel pathways and regulators in human disease. Approach and Results: Of 716 articles published in ATVB from 1995 to 2019 using the apolipoprotein E knockout mouse to study atherosclerosis, data were extracted from 360 unique studies in which a gene was experimentally perturbed to impact plaque size or composition and analyzed using Ingenuity Pathway Analysis software. TREM1 (triggering receptor expressed on myeloid cells) signaling and LXR/RXR (liver X receptor/retinoid X receptor) activation were identified as the top atherosclerosis-associated pathways in mice (both P<1.93×10-4, TREM1 implicated early and LXR/RXR in late atherogenesis). The top upstream regulatory network in mice (sc-58125, a COX2 inhibitor) linked 64.0% of the genes into a single network. The pathways and networks identified in mice were interrogated by testing for associations between the genetically predicted gene expression of each mouse pathway-identified human homolog with clinical atherosclerosis in a cohort of 88 660 human subjects. Homologous human pathways and networks were significantly enriched for gene-atherosclerosis associations (empirical P<0.01 for TREM1 and LXR/RXR pathways and COX2 network). This included 12(60.0%) TREM1 pathway genes, 15(53.6%) LXR/RXR pathway genes, and 67(49.3%) COX2 network genes. Mouse analyses predicted, and human study validated, the strong association of COX2 expression (PTGS2) with increased likelihood of atherosclerosis (odds ratio, 1.68 per SD of genetically predicted gene expression; P=1.07×10-6). CONCLUSIONS: PRESCIANT (Preclinical Science Integration and Translation) leverages published preclinical investigations to identify high-confidence pathways, networks, and regulators of human disease.

Can Treatment Support Mitigate Nicotine Metabolism-Based Disparities in Smoking Abstinence? Secondary Analysis of the Helping HAND 4 Trial.

INTRODUCTION: The nicotine metabolite ratio (NMR), a biomarker of CYP2A6-mediated nicotine metabolism, predicts the efficacy of nicotine replacement therapy (NRT), with fast metabolizers benefiting less than slow metabolizers. Whether treatment support to optimize NRT use (henceforth "treatment support") modifies this pharmacogenetic relationship is unknown. METHODS: Hospitalized adult daily smokers were assigned to one of two post-discharge smoking cessation interventions offering NRT and counseling: (1) Transitional Tobacco Care Management, which delivered enhanced treatment support via free combination NRT at discharge and automated counseling, and (2) a quitline-based approach representing usual care (UC). The primary outcome was biochemically verified 7-day point prevalence abstinence 6 months after discharge. Secondary outcomes were the use of NRT and counseling during the 3-month intervention period. Logistic regression models tested for interactions between NMR and intervention, controlling for sex, race, alcohol use, and BMI. RESULTS: Participants (N = 321) were classified as slow (n = 80) or fast (n = 241) metabolizers relative to the first quartile of NMR (0.012-0.219 vs. 0.221-3.455, respectively). Under UC, fast (vs. slow) metabolizers had lower odds of abstinence at 6 months (aOR 0.35, 95% CI 0.13-0.95) and similar odds of NRT and counseling use. Compared to UC, enhanced treatment support increased abstinence (aOR 2.13, 95% CI 0.98-4.64) and use of combination NRT (aOR 4.62, 95% CI 2.57-8.31) in fast metabolizers, while reducing abstinence in slow metabolizers (aOR 0.21, 95% CI 0.05-0.87; NMR-by-intervention interaction p = .004). CONCLUSIONS: Treatment support increased abstinence and optimal use of NRT among fast nicotine metabolizers, thereby mitigating the gap in abstinence between fast and slow metabolizers. IMPLICATIONS: In this secondary analysis of two smoking cessation interventions for recently hospitalized smokers, fast nicotine metabolizers quit at lower rates than slow metabolizers, but providing fast metabolizers with enhanced treatment support doubled the odds of quitting in this group and mitigated the disparity in abstinence between fast and slow metabolizers. If validated, these findings could lead to personalized approaches to smoking cessation treatment that improve outcomes by targeting treatment support to those who need it most.

A MUC5B Gene Polymorphism, rs35705950-T, Confers Protective Effects Against COVID-19 Hospitalization but Not Severe Disease or Mortality.

Rationale: A common MUC5B gene polymorphism, rs35705950-T, is associated with idiopathic pulmonary fibrosis (IPF), but its role in severe acute respiratory syndrome coronavirus 2 infection and disease severity is unclear. Objectives: To assess whether rs35705950-T confers differential risk for clinical outcomes associated with coronavirus disease (COVID-19) infection among participants in the Million Veteran Program (MVP). Methods: The MUC5B rs35705950-T allele was directly genotyped among MVP participants; clinical events and comorbidities were extracted from the electronic health records. Associations between the incidence or severity of COVID-19 and rs35705950-T were analyzed within each ancestry group in the MVP followed by transancestry meta-analysis. Replication and joint meta-analysis were conducted using summary statistics from the COVID-19 Host Genetics Initiative (HGI). Sensitivity analyses with adjustment for additional covariates (body mass index, Charlson comorbidity index, smoking, asbestosis, rheumatoid arthritis with interstitial lung disease, and IPF) and associations with post-COVID-19 pneumonia were performed in MVP subjects. Measurements and Main Results: The rs35705950-T allele was associated with fewer COVID-19 hospitalizations in transancestry meta-analyses within the MVP (Ncases = 4,325; Ncontrols = 507,640; OR = 0.89 [0.82-0.97]; P = 6.86 × 10-3) and joint meta-analyses with the HGI (Ncases = 13,320; Ncontrols = 1,508,841; OR, 0.90 [0.86-0.95]; P = 8.99 × 10-5). The rs35705950-T allele was not associated with reduced COVID-19 positivity in transancestry meta-analysis within the MVP (Ncases = 19,168/Ncontrols = 492,854; OR, 0.98 [0.95-1.01]; P = 0.06) but was nominally significant (P < 0.05) in the joint meta-analysis with the HGI (Ncases = 44,820; Ncontrols = 1,775,827; OR, 0.97 [0.95-1.00]; P = 0.03). Associations were not observed with severe outcomes or mortality. Among individuals of European ancestry in the MVP, rs35705950-T was associated with fewer post-COVID-19 pneumonia events (OR, 0.82 [0.72-0.93]; P = 0.001). Conclusions: The MUC5B variant rs35705950-T may confer protection in COVID-19 hospitalizations.

Genome-wide association study reveals novel genetic loci: a new polygenic risk score for mitral valve prolapse.

AIMS: Mitral valve prolapse (MVP) is a common valvular heart disease with a prevalence of >2% in the general adult population. Despite this high incidence, there is a limited understanding of the molecular mechanism of this disease, and no medical therapy is available for this disease. We aimed to elucidate the genetic basis of MVP in order to better understand this complex disorder. METHODS AND RESULTS: We performed a meta-analysis of six genome-wide association studies that included 4884 cases and 434 649 controls. We identified 14 loci associated with MVP in our primary analysis and 2 additional loci associated with a subset of the samples that additionally underwent mitral valve surgery. Integration of epigenetic, transcriptional, and proteomic data identified candidate MVP genes including LMCD1, SPTBN1, LTBP2, TGFB2, NMB, and ALPK3. We created a polygenic risk score (PRS) for MVP and showed an improved MVP risk prediction beyond age, sex, and clinical risk factors. CONCLUSION: We identified 14 genetic loci that are associated with MVP. Multiple analyses identified candidate genes including two transforming growth factor-ß signalling molecules and spectrin ß. We present the first PRS for MVP that could eventually aid risk stratification of patients for MVP screening in a clinical setting. These findings advance our understanding of this common valvular heart disease and may reveal novel therapeutic targets for intervention.

Genetic risk for major depressive disorder and loneliness in sex-specific associations with coronary artery disease.

Major depressive disorder (MDD) and loneliness are phenotypically and genetically correlated with coronary artery disease (CAD), but whether these associations are explained by pleiotropic genetic variants or shared comorbidities is unclear. To tease apart these scenarios, we first assessed the medical morbidity pattern associated with genetic risk factors for MDD and loneliness by conducting a phenome-wide association study in 18,385 European-ancestry individuals in the Vanderbilt University Medical Center biobank, BioVU. Polygenic scores for MDD and loneliness were developed for each person using previously published meta-GWAS summary statistics, and were tested for association with 882 clinical diagnoses ascertained via billing codes in electronic health records. We discovered strong associations with heart disease diagnoses, and next embarked on targeted analyses of CAD in 3893 cases and 4197 controls. We found odds ratios of 1.11 (95% CI, 1.04-1.18; P 8.43 × 10-4) and 1.13 (95% CI, 1.07-1.20; P 4.51 × 10-6) per 1-SD increase in the polygenic scores for MDD and loneliness, respectively. Results were similar in patients without psychiatric symptoms, and the increased risk persisted in females even after adjusting for multiple conventional risk factors and a polygenic score for CAD. In a final sensitivity analysis, we statistically adjusted for the genetic correlation between MDD and loneliness and re-computed polygenic scores. The polygenic score unique to loneliness remained associated with CAD (OR 1.09, 95% CI 1.03-1.15; P 0.002), while the polygenic score unique to MDD did not (OR 1.00, 95% CI 0.95-1.06; P 0.97). Our replication sample was the Atherosclerosis Risk in Communities (ARIC) cohort of 7197 European-ancestry participants (1598 incident CAD cases). In ARIC, polygenic scores for MDD and loneliness were associated with hazard ratios of 1.07 (95% CI, 0.99-1.14; P = 0.07) and 1.07 (1.01-1.15; P = 0.03), respectively, and we replicated findings from the BioVU sensitivity analyses. We conclude that genetic risk factors for MDD and loneliness act pleiotropically to increase CAD risk in females.

Deep learning analysis of electrocardiogram for risk prediction of drug-induced arrhythmias and diagnosis of long QT syndrome.

AIMS: Congenital long-QT syndromes (cLQTS) or drug-induced long-QT syndromes (diLQTS) can cause torsade de pointes (TdP), a life-threatening ventricular arrhythmia. The current strategy for the identification of drugs at the high risk of TdP relies on measuring the QT interval corrected for heart rate (QTc) on the electrocardiogram (ECG). However, QTc has a low positive predictive value. METHODS AND RESULTS: We used convolutional neural network (CNN) models to quantify ECG alterations induced by sotalol, an IKr blocker associated with TdP, aiming to provide new tools (CNN models) to enhance the prediction of drug-induced TdP (diTdP) and diagnosis of cLQTS. Tested CNN models used single or multiple 10-s recordings/patient using 8 leads or single leads in various cohorts: 1029 healthy subjects before and after sotalol intake (n = 14 135 ECGs); 487 cLQTS patients (n = 1083 ECGs: 560 type 1, 456 type 2, 67 type 3); and 48 patients with diTdP (n = 1105 ECGs, with 147 obtained within 48 h of a diTdP episode). CNN models outperformed models using QTc to identify exposure to sotalol [area under the receiver operating characteristic curve (ROC-AUC) = 0.98 vs. 0.72, P ≤ 0.001]. CNN models had higher ROC-AUC using multiple vs. single 10-s ECG (P ≤ 0.001). Performances were comparable for 8-lead vs. single-lead models. CNN models predicting sotalol exposure also accurately detected the presence and type of cLQTS vs. healthy controls, particularly for cLQT2 (AUC-ROC = 0.9) and were greatest shortly after a diTdP event and declining over time (P ≤ 0.001), after controlling for QTc and intake of culprit drugs. ECG segment analysis identified the J-Tpeak interval as the best discriminator of sotalol intake. CONCLUSION: CNN models applied to ECGs outperform QTc measurements to identify exposure to drugs altering the QT interval, congenital LQTS, and are greatest shortly after a diTdP episode.

A retrospective approach to evaluating potential adverse outcomes associated with delay of procedures for cardiovascular and cancer-related diagnoses in the context of COVID-19.

OBJECTIVE: During the COVID-19 pandemic, health systems postponed non-essential medical procedures to accommodate surge of critically-ill patients. The long-term consequences of delaying procedures in response to COVID-19 remains unknown. We developed a high-throughput approach to understand the impact of delaying procedures on patient health outcomes using electronic health record (EHR) data. MATERIALS AND METHODS: We used EHR data from Vanderbilt University Medical Center's (VUMC) Research and Synthetic Derivatives. Elective procedures and non-urgent visits were suspended at VUMC between March 18, 2020 and April 24, 2020. Surgical procedure data from this period were compared to a similar timeframe in 2019. Potential adverse impact of delay in cardiovascular and cancer-related procedures was evaluated using EHR data collected from January 1, 1993 to March 17, 2020. For surgical procedure delay, outcomes included length of hospitalization (days), mortality during hospitalization, and readmission within six months. For screening procedure delay, outcomes included 5-year survival and cancer stage at diagnosis. RESULTS: We identified 416 surgical procedures that were negatively impacted during the COVID-19 pandemic compared to the same timeframe in 2019. Using retrospective data, we found 27 significant associations between procedure delay and adverse patient outcomes. Clinician review indicated that 88.9% of the significant associations were plausible and potentially clinically significant. Analytic pipelines for this study are available online. CONCLUSION: Our approach enables health systems to identify medical procedures affected by the COVID-19 pandemic and evaluate the effect of delay, enabling them to communicate effectively with patients and prioritize rescheduling to minimize adverse patient outcomes.

Prognostic implications of pre-existing medical comorbidity in Takotsubo cardiomyopathy.

Takotsubo cardiomyopathy (TC) is associated with significant short-term morbidity and mortality. Several risk factors for poor outcomes have been identified; however, the prognostic implications of pre-existing comorbidity in TC are poorly delineated. We sought to assess the association of aggregate pre-existing comorbidity with short-term outcomes in TC. We performed a retrospective observational study of adult subjects diagnosed with TC at two academic tertiary care hospitals between 2005 and 2018. Overall burden of medical comorbidity was estimated using the Charlson comorbidity index (CCI). Multivariable logistic regression was used to test for independent association of CCI with 30-day mortality and severe shock at index presentation. Multivariable poisson regression was performed to assess the association of CCI with duration of hospitalization. Five-hundred and thirty-eight subjects were diagnosed with TC during the study period. The median CCI score of all subjects was 2 (IQR 1-4). Among subjects with physical triggers of TC, the median CCI score was 2 (IQR 1-4) compared to a median CCI score of 1 (IQR 0-1) in subjects with non-physical triggers of TC (P < 0.001). Seventy-six (14%) subjects died within 30 days of index diagnosis and 185 (34%) subjects experienced severe shock. The median duration of hospitalization was 7 days (IQR 3-14 days). In multivariable logistic regression, CCI was not associated with 30-day mortality or severe shock. In multivariable Poisson regression, CCI (IRR 1.17, 95% CI 1.16-1.18, P < 0.001) was associated with duration of hospitalization. Increased burden of pre-existing medical comorbidity was not independently associated with 30-day mortality or severe shock at index presentation, but was associated with increased duration of hospitalization after diagnosis of TC.