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.

The relationship between high density lipoprotein cholesterol and sepsis: A clinical and genetic approach.

Sepsis accounts for one in three hospital deaths. Higher concentrations of high-density lipoprotein cholesterol (HDL-C) are associated with apparent protection from sepsis, suggesting a potential therapeutic role for HDL-C or drugs, such as cholesteryl ester transport protein (CETP) inhibitors that increase HDL-C. However, these beneficial clinical associations might be due to confounding; genetic approaches can address this possibility. We identified 73,406 White adults admitted to Vanderbilt University Medical Center with infection; 11,612 had HDL-C levels, and 12,377 had genotype information from which we constructed polygenic risk scores (PRS) for HDL-C and the effect of CETP on HDL-C. We tested the associations between predictors (measured HDL-C, HDL-C PRS, CETP PRS, and rs1800777) and outcomes: sepsis, septic shock, respiratory failure, and in-hospital death. In unadjusted analyses, lower measured HDL-C concentrations were significantly associated with increased risk of sepsis (p = 2.4 × 10-23 ), septic shock (p = 4.1 × 10-12 ), respiratory failure (p = 2.8 × 10-8 ), and in-hospital death (p = 1.0 × 10-8 ). After adjustment (age, sex, electronic health record length, comorbidity score, LDL-C, triglycerides, and body mass index), these associations were markedly attenuated: sepsis (p = 2.6 × 10-3 ), septic shock (p = 8.1 × 10-3 ), respiratory failure (p = 0.11), and in-hospital death (p = 4.5 × 10-3 ). HDL-C PRS, CETP PRS, and rs1800777 significantly predicted HDL-C (p < 2 × 10-16 ), but none were associated with sepsis outcomes. Concordant findings were observed in 13,254 Black patients hospitalized with infections. Lower measured HDL-C levels were significantly associated with increased risk of sepsis and related outcomes in patients with infection, but a causal relationship is unlikely because no association was found between the HDL-C PRS or the CETP PRS and the risk of adverse sepsis outcomes.

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.

Genome-wide association analyses of common infections in a large practice-based biobank.

INTRODUCTION: Infectious diseases are common causes of morbidity and mortality worldwide. Susceptibility to infection is highly heritable; however, little has been done to identify the genetic determinants underlying common infectious diseases. One GWAS was performed using 23andMe information about self-reported infections; we set out to confirm previous loci and identify new ones using medically diagnosed infections. METHODS: We used the electronic health record (EHR)-based biobank at Vanderbilt and diagnosis codes to identify cases of 12 infectious diseases in white patients: urinary tract infection, pneumonia, chronic sinus infections, otitis media, candidiasis, streptococcal pharyngitis, herpes zoster, herpes labialis, hepatitis B, infectious mononucleosis, tuberculosis (TB) or a positive TB test, and hepatitis C. We selected controls from patients with no diagnosis code for the candidate disease and matched by year of birth, sex, and calendar year at first and last EHR visits. We conducted GWAS using SAIGE and transcriptome-wide analysis (TWAS) using S-PrediXcan. We also conducted phenome-wide association study to understand associations between identified genetic variants and clinical phenotypes. RESULTS: We replicated three 23andMe loci (p ≤ 0.05): herpes zoster and rs7047299-A (p = 2.6 × 10-3) and rs2808290-C (p = 9.6 × 10-3;); otitis media and rs114947103-C (p = 0.04). We also identified 2 novel regions (p ≤ 5 × 10-8): rs113235453-G for otitis media (p = 3.04 × 10-8), and rs10422015-T for candidiasis (p = 3.11 × 10-8). In TWAS, four gene-disease associations were significant: SLC30A9 for otitis media (p = 8.06 × 10-7); LRP3 and WDR88 for candidiasis (p = 3.91 × 10-7 and p = 1.95 × 10-6); and AAMDC for hepatitis B (p = 1.51 × 10-6). CONCLUSION: We conducted GWAS and TWAS for 12 infectious diseases and identified novel genetic contributors to the susceptibility of infectious diseases.

Genetic Determinants of Body Mass Index and Fasting Glucose Are Mediators of Grade 1 Diastolic Dysfunction.

Background Early (grade 1) cardiac left ventricular diastolic dysfunction (G1DD) increases the risk for heart failure with preserved ejection fraction and may improve with aggressive risk factor modification. Type 2 diabetes, obesity, hypertension, and coronary heart disease are associated with increased incidence of diastolic dysfunction. The genetic drivers of G1DD are not defined. Methods and Results We curated genotyped European ancestry G1DD cases (n=668) and controls with normal diastolic function (n=1772) from Vanderbilt's biobank. G1DD status was explored through (1) an additive model genome-wide association study, (2) shared polygenic risk through logistic regression, and (3) instrumental variable analysis using 2-sample Mendelian randomization (the inverse-variance weighted method, Mendelian randomization-Egger, and median) to determine potential modifiable risk factors. There were no common single nucleotide polymorphisms significantly associated with G1DD status. A polygenic risk score for BMI was significantly associated with increased G1DD risk (odds ratio [OR], 1.20 for 1-SD increase in BMI [95% CI, 1.08-1.32]; P=0.0003). The association was confirmed by the inverse-variance weighted method (OR, 1.89 [95% CI, 1.37-2.61]). Among the candidate mediators for BMI, only fasting glucose was significantly associated with G1DD status by the inverse-variance weighted method (OR, 4.14 for 1-SD increase in fasting glucose [95% CI, 1.55-11.02]; P=0.005). Multivariable Mendelian randomization showed a modest attenuation of the BMI association (OR, 1.84 [95% CI, 1.35-2.52]) when adjusting for fasting glucose. Conclusions These data suggest that a genetic predisposition to elevated BMI increases the risk for G1DD. Part of this effect may be mediated through altered glucose homeostasis.

Genome-wide association study of platelet factor 4/heparin antibodies in heparin-induced thrombocytopenia.

Heparin, a widely used anticoagulant, carries the risk of an antibody-mediated adverse drug reaction, heparin-induced thrombocytopenia (HIT). A subset of heparin-treated patients produces detectable levels of antibodies against complexes of heparin bound to circulating platelet factor 4 (PF4). Using a genome-wide association study (GWAS) approach, we aimed to identify genetic variants associated with anti-PF4/heparin antibodies that account for the variable antibody response seen in HIT. We performed a GWAS on anti-PF4/heparin antibody levels determined via polyclonal enzyme-linked immunosorbent assays. Our discovery cohort (n = 4237) and replication cohort (n = 807) constituted patients with European ancestry and clinical suspicion of HIT, with cases confirmed via functional assay. Genome-wide significance was considered at α = 5 × 10-8. No variants were significantly associated with anti-PF4/heparin antibody levels in the discovery cohort at a genome-wide significant level. Secondary GWAS analyses included the identification of variants with suggestive associations in the discovery cohort (α = 1 × 10-4). The top variant in both cohorts was rs1555175145 (discovery ß = -0.112 [0.018], P = 2.50 × 10-5; replication ß = -0.104 [0.051], P = .041). In gene set enrichment analysis, 3 gene sets reached false discovery rate-adjusted significance (q < 0.05) in both discovery and replication cohorts: "Leukocyte Transendothelial Migration," "Innate Immune Response," and "Lyase Activity." Our results indicate that genomic variation is not significantly associated with anti-PF4/heparin antibody levels. Given our power to identify variants with moderate frequencies and effect sizes, this evidence suggests genetic variation is not a primary driver of variable antibody response in heparin-treated patients with European ancestry.

Composite CYP3A phenotypes influence tacrolimus dose-adjusted concentration in lung transplant recipients.

OBJECTIVES: Interpatient variability in tacrolimus pharmacokinetics is attributed to metabolism by cytochrome P-450 3A4/5 isoenzymes (encoded by CYP3A4 and CYP3A5). Guidelines for adjusting tacrolimus based on CYP3A5 test results are published; however, CYP3A4 variants also contribute to the variability in tacrolimus pharmacokinetics. The effects of composite phenotypes incorporating CYP3A5 and CYP3A4 increased (*1G, *1B) and decreased (*22) function variants have not been evaluated. The objective of this study is to investigate the impact of both increased and decreased function CYP3A variants on weight and dose-adjusted tacrolimus concentration (C0/D). METHODS: We performed a single-center retrospective cohort study of lung transplant recipients to evaluate the median tacrolimus C0/D by composite CYP3A phenotype groups during the index transplant hospitalization. CYP3A4 and CYP3A5 alleles were used to classify patients into four CYP3A groups from least to most CYP3A activity. Exploratory analyses of ABCB1 and additional candidate genes were also assessed. RESULTS: Of the 92 included individuals, most (58) were CYP3A Group 2. The median tacrolimus C0/D differed significantly between CYP3A groups (P = 0.0001). CYP3A Group 2 median tacrolimus C0/D was 190.5 (interquartile range: 147.6-267.5) (ng/ml)/(mg/kg/d) and significantly higher than Group 4 [107.9 (90.4-116.1), P = 0.0001)]. Group 2 median tacrolimus C0/D did not significantly differ from Group 1 and Group 3 [373.5 (149.2-490.3) and 81.4 (62.6-184.1), respectively]. No significant differences in tacrolimus C0/D were found for the ABCB1 diplotypes. CONCLUSION: These data indicate that a composite CYP3A phenotype incorporating both increase and decrease variant information from CYP3A4 in addition to CYP3A5 may significantly influence tacrolimus C0/D during the early postoperative period.

ABO O blood group as a risk factor for platelet reactivity in heparin-induced thrombocytopenia.

Heparin-induced thrombocytopenia (HIT) is an unpredictable, potentially catastrophic adverse effect resulting from an immune response to platelet factor 4 (PF4)/heparin complexes. We performed a genome-wide association study (GWAS) with positive functional assay as the outcome in a large discovery cohort of patients divided into 3 groups: (1) functional assay-positive cases (n = 1269), (2) antibody-positive (functional assay-negative) controls (n = 1131), and (3) antibody-negative controls (n = 1766). Significant associations (α = 5 × 10-8) were investigated in a replication cohort (α = 0.05) of functional assay-confirmed HIT cases (n = 177), antibody-positive (function assay-negative) controls (n = 258), and antibody-negative controls (n = 351). We observed a strong association for positive functional assay with increasing PF4/heparin immunoglobulin-G (IgG) level (odds ratio [OR], 16.53; 95% confidence interval [CI], 13.83-19.74; P = 1.51 × 10-209) and female sex (OR, 1.15; 95% CI, 1.01-1.32; P = .034). The rs8176719 C insertion variant in ABO was significantly associated with positive functional assay status in the discovery cohort (frequency = 0.41; OR, 0.751; 95% CI, 0.682-0.828; P = 7.80 × 10-9) and in the replication cohort (OR, 0.467; 95% CI, 0.228-0.954; P = .0367). The rs8176719 C insertion, which encodes all non-O blood group alleles, had a protective effect, indicating that the rs8176719 C deletion and the O blood group were risk factors for HIT (O blood group OR, 1.42; 95% CI, 1.26-1.61; P = 3.09 × 10-8). Meta-analyses indicated that the ABO association was independent of PF4/heparin IgG levels and was stronger when functional assay-positive cases were compared with antibody-positive (functional assay-negative) controls than with antibody-negative controls. Sequencing and fine-mapping of ABO demonstrated that rs8176719 was the causal single nucleotide polymorphism (SNP). Our results clarify the biology underlying HIT pathogenesis with ramifications for prediction and may have important implications for related conditions, such as vaccine-induced thrombotic thrombocytopenia.

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.

Early-Onset Atrial Fibrillation and the Prevalence of Rare Variants in Cardiomyopathy and Arrhythmia Genes.

Importance: Early-onset atrial fibrillation (AF) can be the initial manifestation of a more serious underlying inherited cardiomyopathy or arrhythmia syndrome. Objective: To examine the results of genetic testing for early-onset AF. Design, Setting, and Participants: This prospective, observational cohort study enrolled participants from an academic medical center who had AF diagnosed before 66 years of age and underwent whole genome sequencing through the National Heart, Lung, and Blood Institute's Trans-Omics for Precision Medicine program. Participants were enrolled from November 23, 1999, to June 2, 2015. Data analysis was performed from October 24, 2020, to March 11, 2021. Exposures: Rare variants identified in a panel of 145 genes that are included on cardiomyopathy and arrhythmia panels used by commercial clinical genetic testing laboratories. Main Outcomes and Measures: Sequencing data were analyzed using an automated process followed by manual review by a panel of independent, blinded reviewers. The primary outcome was classification of rare variants using American College of Medical Genetics and Genomics criteria: benign, likely benign, variant of undetermined significance, likely pathogenic, or pathogenic. Disease-associated variants were defined as pathogenic/likely pathogenic variants in genes associated with autosomal dominant or X-linked dominant disorders. Results: Among 1293 participants (934 [72.2%] male; median [interquartile range] age at enrollment, 56 [48-61] years; median [interquartile range] age at AF diagnosis, 50 [41-56] years), genetic testing identified 131 participants (10.1%) with a disease-associated variant, 812 (62.8%) with a variant of undetermined significance, 92 (7.1%) as heterozygous carriers for an autosomal recessive disorder, and 258 (20.0%) with no suspicious variant. The likelihood of a disease-associated variant was highest in participants with AF diagnosed before the age of 30 years (20 of 119 [16.8%; 95% CI, 10.0%-23.6%]) and lowest after the age of 60 years (8 of 112 [7.1%; 95% CI, 2.4%-11.9%]). Disease-associated variants were more often associated with inherited cardiomyopathy syndromes compared with inherited arrhythmias. The most common genes were TTN (n = 38), MYH7 (n = 18), MYH6 (n = 10), LMNA (n = 9), and KCNQ1 (n = 8). Conclusions and Relevance: In this cohort study, genetic testing identified a disease-associated variant in 10% of patients with early-onset AF (the percentage was higher if diagnosed before the age of 30 years and lower if diagnosed after the age of 60 years). Most pathogenic/likely pathogenic variants are in genes associated with cardiomyopathy. These results support the use of genetic testing in early-onset AF.

The genetic architecture of plasma kynurenine includes cardiometabolic disease mechanisms associated with the SH2B3 gene.

Inflammation increases the risk of cardiometabolic disease. Delineating specific inflammatory pathways and biomarkers of their activity could identify the mechanistic underpinnings of the increased risk. Plasma levels of kynurenine, a metabolite involved in inflammation, associates with cardiometabolic disease risk. We used genetic approaches to identify inflammatory mechanisms associated with kynurenine variability and their relationship to cardiometabolic disease. We identified single-nucleotide polymorphisms (SNPs) previously associated with plasma kynurenine, including a missense-variant (rs3184504) in the inflammatory gene SH2B3/LNK. We examined the association between rs3184504 and plasma kynurenine in independent human samples, and measured kynurenine levels in SH2B3-knock-out mice and during human LPS-evoked endotoxemia. We conducted phenome scanning to identify clinical phenotypes associated with each kynurenine-related SNP and with a kynurenine polygenic score using the UK-Biobank (n = 456,422), BioVU (n = 62,303), and Electronic Medical Records and Genetics (n = 32,324) databases. The SH2B3 missense variant associated with plasma kynurenine levels and SH2B3-/- mice had significant tissue-specific differences in kynurenine levels.LPS, an acute inflammatory stimulus, increased plasma kynurenine in humans. Mendelian randomization showed increased waist-circumference, a marker of central obesity, associated with increased kynurenine, and increased kynurenine associated with C-reactive protein (CRP). We found 30 diagnoses associated (FDR q < 0.05) with the SH2B3 variant, but not with SNPs mapping to genes known to regulate tryptophan-kynurenine metabolism. Plasma kynurenine may be a biomarker of acute and chronic inflammation involving the SH2B3 pathways. Its regulation lies upstream of CRP, suggesting that kynurenine may be a biomarker of one inflammatory mechanism contributing to increased cardiometabolic disease risk.

Polygenic Risk Score to Identify Subclinical Coronary Heart Disease Risk in Young Adults.

BACKGROUND: Polygenic risk scores (PRS) may enhance risk stratification for coronary heart disease among young adults. Whether a coronary heart disease PRS improves prediction beyond modifiable risk factors in this population is not known. METHODS: Genotyped adults aged 18 to 35 years were selected from the CARDIA study (Coronary Artery Risk Development in Young Adults; n=1132) and FOS (Framingham Offspring Study; n=663). Systolic blood pressure, total and HDL (high-density lipoprotein) cholesterol, triglycerides, smoking, and waist circumference or body mass index were measured at the visit 1 exam of each study, and coronary artery calcium, a measure of coronary atherosclerosis, was assessed at year 15 (CARDIA) or year 30 (FOS). A previously validated PRS for coronary heart disease was computed for each subject. The C statistic and integrated discrimination improvement were used to compare improvements in prediction of elevated coronary artery calcium between models containing the PRS, risk factors, or both. RESULTS: There were 62 (5%) and 93 (14%) participants with a coronary artery calcium score >20 (CARDIA) and >300 (FOS), respectively. At these thresholds, the C statistic changes of adding the PRS to a risk factor-based model were 0.015 (0.004-0.028) and 0.020 (0.001-0.039) in CARDIA and FOS, respectively. When adding risk factors to a PRS-based model, the respective changes were 0.070 (0.033-0.109) and 0.051 (0.017-0.079). The integrated discrimination improvement, when adding the PRS to a risk factor model, was 0.027 (-0.006 to 0.054) in CARDIA and 0.039 (0.0005-0.072) in FOS. CONCLUSIONS: Among young adults, a PRS improved model discrimination for coronary atherosclerosis, but improvements were smaller than those associated with modifiable risk factors.

High-throughput framework for genetic analyses of adverse drug reactions using electronic health records.

Understanding the contribution of genetic variation to drug response can improve the delivery of precision medicine. However, genome-wide association studies (GWAS) for drug response are uncommon and are often hindered by small sample sizes. We present a high-throughput framework to efficiently identify eligible patients for genetic studies of adverse drug reactions (ADRs) using "drug allergy" labels from electronic health records (EHRs). As a proof-of-concept, we conducted GWAS for ADRs to 14 common drug/drug groups with 81,739 individuals from Vanderbilt University Medical Center's BioVU DNA Biobank. We identified 7 genetic loci associated with ADRs at P < 5 × 10-8, including known genetic associations such as CYP2D6 and OPRM1 for CYP2D6-metabolized opioid ADR. Additional expression quantitative trait loci and phenome-wide association analyses added evidence to the observed associations. Our high-throughput framework is both scalable and portable, enabling impactful pharmacogenomic research to improve precision medicine.

Association Between a Common, Benign Genotype and Unnecessary Bone Marrow Biopsies Among African American Patients.

Importance: Up to two-thirds of African American individuals carry the benign rs2814778-CC genotype that lowers total white blood cell (WBC) count. Objective: To examine whether the rs2814778-CC genotype is associated with an increased likelihood of receiving a bone marrow biopsy (BMB) for an isolated low WBC count. Design, Setting, and Participants: This retrospective genetic association study assessed African American patients younger than 90 years who underwent a BMB at Vanderbilt University Medical Center, Mount Sinai Health System, or Children's Hospital of Philadelphia from January 1, 1998, to December 31, 2020. Exposure: The rs2814778-CC genotype. Main Outcomes and Measures: The proportion of individuals with the CC genotype who underwent BMB for an isolated low WBC count and had a normal biopsy result compared with the proportion of individuals with the CC genotype who underwent BMB for other indications and had a normal biopsy result. Results: Among 399 individuals who underwent a BMB (mean [SD] age, 41.8 [22.5] years, 234 [59%] female), 277 (69%) had the CC genotype. A total of 35 patients (9%) had clinical histories of isolated low WBC counts, and 364 (91%) had other histories. Of those with a clinical history of isolated low WBC count, 34 of 35 (97%) had the CC genotype vs 243 of 364 (67%) of those without a low WBC count history. Among those with the CC genotype, 33 of 34 (97%) had normal results for biopsies performed for isolated low WBC counts compared with 134 of 243 individuals (55%) with biopsies performed for other histories (P < .001). Conclusions and Relevance: In this genetic association study, among patients of African American race who had a BMB with a clinical history of isolated low WBC counts, the rs2814778-CC genotype was highly prevalent, and 97% of these BMBs identified no hematologic abnormality. Accounting for the rs2814778-CC genotype in clinical decision-making could avoid unnecessary BMB procedures.

Genetic Thyrotropin Regulation of Atrial Fibrillation Risk Is Mediated Through an Effect on Height.

CONTEXT: A genetic predisposition to lower thyrotropin (TSH) levels is associated with increased atrial fibrillation (AF) risk through undefined mechanisms. OBJECTIVE: Defining the genetic mediating mechanisms could lead to improved targeted therapies to mitigate AF risk. METHODS: We used 2-sample mendelian randomization (MR) to test associations between TSH-associated single-nucleotide variations and 16 candidate mediators. We then performed multivariable mendelian randomization (MVMR) to test for a significant attenuation of the genetic association between TSH and AF, after adjusting for each mediator significantly associated with TSH. RESULTS: Four candidate mediators (free thyroxine, systolic blood pressure, heart rate, and height) were significantly inversely associated with genetically predicted TSH after adjusting for multiple testing. In MVMR analyses, adjusting for height significantly decreased the magnitude of the association between TSH and AF from -0.12 (SE 0.02) occurrences of AF per SD change in height to -0.06 (0.02) (P = .005). Adjusting for the other candidate mediators did not significantly attenuate the association. CONCLUSION: The genetic association between TSH and increased AF risk is mediated, in part, by taller stature. Thus, some genetic mechanisms underlying TSH variability may contribute to AF risk through mechanisms determining height occurring early in life that differ from those driven by thyroid hormone-level elevations in later life.

Variants in ADRB1 and CYP2C9: Association with Response to Atenolol and Losartan in Marfan Syndrome.

OBJECTIVE: To test whether variants in ADRB1 and CYP2C9 genes identify subgroups of individuals with differential response to treatment for Marfan syndrome through analysis of data from a large, randomized trial. STUDY DESIGN: In a subset of 250 white, non-Hispanic participants with Marfan syndrome in a prior randomized trial of atenolol vs losartan, the common variants rs1801252 and rs1801253 in ADRB1 and rs1799853 and rs1057910 in CYP2C9 were analyzed. The primary outcome was baseline-adjusted annual rate of change in the maximum aortic root diameter z-score over 3 years, assessed using mixed effects models. RESULTS: Among 122 atenolol-assigned participants, the 70 with rs1801253 CC genotype had greater rate of improvement in aortic root z-score compared with 52 participants with CG or GG genotypes (Time × Genotype interaction P = .005, mean annual z-score change ± SE -0.20 ± 0.03 vs -0.09 ± 0.03). Among participants with the CC genotype in both treatment arms, those assigned to atenolol had greater rate of improvement compared with the 71 of the 121 assigned to losartan (interaction P = .002; -0.20 ± 0.02 vs -0.07 ± 0.02; P < .001). There were no differences in atenolol response by rs1801252 genotype or in losartan response by CYP2C9 metabolizer status. CONCLUSIONS: In this exploratory study, ADRB1-rs1801253 was associated with atenolol response in children and young adults with Marfan syndrome. If these findings are confirmed in future studies, ADRB1 genotyping has the potential to guide therapy by identifying those who are likely to have greater therapeutic response to atenolol than losartan.

The polygenic architecture of left ventricular mass mirrors the clinical epidemiology.

Left ventricular (LV) mass is a prognostic biomarker for incident heart disease and all-cause mortality. Large-scale genome-wide association studies have identified few SNPs associated with LV mass. We hypothesized that a polygenic discovery approach using LV mass measurements made in a clinical population would identify risk factors and diseases associated with adverse LV remodeling. We developed a polygenic single nucleotide polymorphism-based predictor of LV mass in 7,601 individuals with LV mass measurements made during routine clinical care. We tested for associations between this predictor and 894 clinical diagnoses measured in 58,838 unrelated genotyped individuals. There were 29 clinical phenotypes associated with the LV mass genetic predictor at FDR q < 0.05. Genetically predicted higher LV mass was associated with modifiable cardiac risk factors, diagnoses related to organ dysfunction and conditions associated with abnormal cardiac structure including heart failure and atrial fibrillation. Secondary analyses using polygenic predictors confirmed a significant association between higher LV mass and body mass index and, in men, associations with coronary atherosclerosis and systolic blood pressure. In summary, these analyses show that LV mass-associated genetic variability associates with diagnoses of cardiac diseases and with modifiable risk factors which contribute to these diseases.

Phenome-wide association analysis suggests the APOL1 linked disease spectrum primarily drives kidney-specific pathways.

The relationship between commonly occurring genetic variants (G1 and G2) in the APOL1 gene in African Americans and different disease traits, such as kidney disease, cardiovascular disease, and pre-eclampsia, remains the subject of controversy. Here we took a genotype-first approach, a phenome-wide association study, to define the spectrum of phenotypes associated with APOL1 high-risk variants in 1,837 African American participants of Penn Medicine Biobank and 4,742 African American participants of Vanderbilt BioVU. In the Penn Medicine Biobank, outpatient creatinine measurement-based estimated glomerular filtration rate and multivariable regression models were used to evaluate the association between high-risk APOL1 status and renal outcomes. In meta-analysis of both cohorts, the strongest phenome-wide association study associations were for the high-risk APOL1 variants and diagnoses codes were highly significant for "kidney dialysis" (odds ratio 3.75) and "end stage kidney disease" (odds ratio 3.42). A number of phenotypes were associated with APOL1 high-risk genotypes in an analysis adjusted only for demographic variables. However, no associations were detected with non-renal phenotypes after controlling for chronic/end stage kidney disease status. Using calculated estimated glomerular filtration rate -based phenotype analysis in the Penn Medicine Biobank, APOL1 high-risk status was associated with prevalent chronic/end stage kidney disease /kidney transplant (odds ratio 2.27, 95% confidence interval 1.67-3.08). In high-risk participants, the estimated glomerular filtration rate was 15.4 mL/min/1.73m2; significantly lower than in low-risk participants. Thus, although APOL1 high-risk variants are associated with a range of phenotypes, the risks for other associated phenotypes appear much lower and in our dataset are driven by a primary effect on renal disease.

Association of FADS1/2 Locus Variants and Polyunsaturated Fatty Acids With Aortic Stenosis.

Importance: Aortic stenosis (AS) has no approved medical treatment. Identifying etiological pathways for AS could identify pharmacological targets. Objective: To identify novel genetic loci and pathways associated with AS. Design, Setting, and Participants: This genome-wide association study used a case-control design to evaluate 44 703 participants (3469 cases of AS) of self-reported European ancestry from the Genetic Epidemiology Research on Adult Health and Aging (GERA) cohort (from January 1, 1996, to December 31, 2015). Replication was performed in 7 other cohorts totaling 256 926 participants (5926 cases of AS), with additional analyses performed in 6942 participants from the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium. Follow-up biomarker analyses with aortic valve calcium (AVC) were also performed. Data were analyzed from May 1, 2017, to December 5, 2019. Exposures: Genetic variants (615 643 variants) and polyunsaturated fatty acids (ω-6 and ω-3) measured in blood samples. Main Outcomes and Measures: Aortic stenosis and aortic valve replacement defined by electronic health records, surgical records, or echocardiography and the presence of AVC measured by computed tomography. Results: The mean (SD) age of the 44 703 GERA participants was 69.7 (8.4) years, and 22 019 (49.3%) were men. The rs174547 variant at the FADS1/2 locus was associated with AS (odds ratio [OR] per C allele, 0.88; 95% CI, 0.83-0.93; P = 3.0 × 10-6), with genome-wide significance after meta-analysis with 7 replication cohorts totaling 312 118 individuals (9395 cases of AS) (OR, 0.91; 95% CI, 0.88-0.94; P = 2.5 × 10-8). A consistent association with AVC was also observed (OR, 0.91; 95% CI, 0.83-0.99; P = .03). A higher ratio of arachidonic acid to linoleic acid was associated with AVC (OR per SD of the natural logarithm, 1.19; 95% CI, 1.09-1.30; P = 6.6 × 10-5). In mendelian randomization, increased FADS1 liver expression and arachidonic acid were associated with AS (OR per unit of normalized expression, 1.31 [95% CI, 1.17-1.48; P = 7.4 × 10-6]; OR per 5-percentage point increase in arachidonic acid for AVC, 1.23 [95% CI, 1.01-1.49; P = .04]; OR per 5-percentage point increase in arachidonic acid for AS, 1.08 [95% CI, 1.04-1.13; P = 4.1 × 10-4]). Conclusions and Relevance: Variation at the FADS1/2 locus was associated with AS and AVC. Findings from biomarker measurements and mendelian randomization appear to link ω-6 fatty acid biosynthesis to AS, which may represent a therapeutic target.