ABSTRACT
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.
Subject(s)
Aortic Valve Stenosis , Dyslipidemias , Humans , Genome-Wide Association Study/methods , Adiposity/genetics , Genetic Predisposition to Disease , Aortic Valve Stenosis/genetics , Obesity , Risk Factors , Inflammation , Dyslipidemias/complications , Dyslipidemias/genetics , Apolipoproteins/genetics , Mendelian Randomization Analysis , Polymorphism, Single Nucleotide/geneticsABSTRACT
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.
Subject(s)
Diabetes Mellitus, Type 2 , Fasting , Body Mass Index , Diabetes Mellitus, Type 2/epidemiology , Genome-Wide Association Study , Glucose , Humans , Mendelian Randomization Analysis , Polymorphism, Single Nucleotide , Risk FactorsABSTRACT
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.
Subject(s)
Platelet Factor 4 , Thrombocytopenia , Antibodies , Genome-Wide Association Study , Heparin/adverse effects , Humans , Immunologic Factors/adverse effects , Platelet Factor 4/genetics , Thrombocytopenia/chemically induced , Thrombocytopenia/geneticsABSTRACT
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.
Subject(s)
Kidney Transplantation , Tacrolimus , Cytochrome P-450 CYP3A/genetics , Cytochrome P-450 CYP3A/metabolism , Genotype , Humans , Immunosuppressive Agents/pharmacokinetics , Lung , Phenotype , Polymorphism, Single Nucleotide , Retrospective Studies , Tacrolimus/pharmacokinetics , Transplant RecipientsABSTRACT
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.
Subject(s)
Genome-Wide Association Study , Thrombocytopenia , ABO Blood-Group System/genetics , Anticoagulants/adverse effects , Female , Heparin/adverse effects , Humans , Immunoglobulin G , Male , Platelet Factor 4/genetics , Risk Factors , Thrombocytopenia/chemically induced , Thrombocytopenia/geneticsABSTRACT
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.
Subject(s)
Arrhythmias, Cardiac/chemically induced , Ion Channels/metabolism , HumansABSTRACT
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.
Subject(s)
Arrhythmias, Cardiac , Genetic Testing , Arrhythmias, Cardiac/diagnosis , Arrhythmias, Cardiac/genetics , Genetic Predisposition to Disease , Genetic Testing/methods , Genomics , HEK293 Cells , Humans , Phenotype , Prospective StudiesABSTRACT
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.
Subject(s)
Atrial Fibrillation/genetics , Cardiomyopathies/genetics , Registries , Whole Genome Sequencing/methods , Adult , Aged , Atrial Fibrillation/epidemiology , Cardiomyopathies/epidemiology , Female , Follow-Up Studies , Humans , Male , Middle Aged , Precision Medicine , Prevalence , Prospective Studies , United States/epidemiologyABSTRACT
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.
Subject(s)
Calcium/blood , Cholesterol, HDL/blood , Coronary Disease/blood , Triglycerides/blood , Adolescent , Adult , Body Mass Index , Cohort Studies , Female , Humans , Male , Risk , Risk Assessment , Smoking/adverse effects , Smoking/blood , Young AdultABSTRACT
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.
Subject(s)
Kynurenine , Polymorphism, Single Nucleotide , Animals , Biomarkers , C-Reactive Protein , Inflammation , Male , Mice , Tryptophan/metabolismABSTRACT
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.
Subject(s)
Biopsy , Black or African American/genetics , Bone Marrow Examination , Duffy Blood-Group System/genetics , Neutropenia , Receptors, Cell Surface/genetics , Adult , Biopsy/methods , Biopsy/statistics & numerical data , Bone Marrow Examination/methods , Bone Marrow Examination/statistics & numerical data , Female , Gene Expression Profiling/statistics & numerical data , Genetic Profile , Genome-Wide Association Study , Humans , Leukocyte Count , Male , Neutropenia/diagnosis , Neutropenia/ethnology , Neutropenia/genetics , Polymorphism, Single Nucleotide , United States/epidemiology , Unnecessary Procedures/methods , Unnecessary Procedures/statistics & numerical dataABSTRACT
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.
Subject(s)
Drug-Related Side Effects and Adverse Reactions/genetics , Electronic Health Records , High-Throughput Nucleotide Sequencing/methods , Genome-Wide Association Study , Humans , Pharmacogenetics , Precision MedicineABSTRACT
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.
Subject(s)
Atrial Fibrillation/genetics , Body Height/genetics , Genetic Predisposition to Disease/genetics , Thyrotropin/blood , Blood Pressure/genetics , Genome-Wide Association Study , Heart Disease Risk Factors , Humans , Mendelian Randomization Analysis , Polymorphism, Single Nucleotide , Thyroid Function Tests , Thyroid Gland/metabolismABSTRACT
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.
Subject(s)
Atenolol/therapeutic use , Cytochrome P-450 CYP2C9/genetics , Gene Expression Regulation , Losartan/therapeutic use , Marfan Syndrome/drug therapy , Receptors, Adrenergic, beta-1/genetics , Adolescent , Adrenergic beta-1 Receptor Antagonists/therapeutic use , Adult , Angiotensin II Type 1 Receptor Blockers/therapeutic use , Child , Child, Preschool , Cytochrome P-450 CYP2C9/biosynthesis , DNA/genetics , Female , Follow-Up Studies , Genotype , Humans , Infant , Male , Marfan Syndrome/genetics , Marfan Syndrome/metabolism , Receptors, Adrenergic, beta-1/biosynthesis , Retrospective Studies , Young AdultABSTRACT
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.
Subject(s)
Genetic Predisposition to Disease , Heart Diseases/epidemiology , Heart Diseases/etiology , Heart Ventricles/pathology , Multifactorial Inheritance , Ventricular Remodeling , Echocardiography , Female , Genome-Wide Association Study , Genotype , Heart Diseases/diagnosis , Heart Diseases/mortality , Heart Ventricles/anatomy & histology , Heart Ventricles/diagnostic imaging , Humans , Male , Odds Ratio , Organ Size , Phenotype , Polymorphism, Single Nucleotide , Risk Assessment , Risk Factors , Ventricular Remodeling/geneticsABSTRACT
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.
Subject(s)
Apolipoprotein L1 , Kidney , Apolipoprotein L1/genetics , Creatinine , Genetic Predisposition to Disease , Genotype , Glomerular Filtration Rate , Humans , Risk FactorsABSTRACT
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.
Subject(s)
Aortic Valve Stenosis/genetics , DNA/genetics , Fatty Acid Desaturases/genetics , Fatty Acids, Unsaturated/genetics , Genetic Predisposition to Disease , Polymorphism, Single Nucleotide , Aged , Alleles , Aortic Valve Stenosis/metabolism , Case-Control Studies , Delta-5 Fatty Acid Desaturase , Fatty Acid Desaturases/metabolism , Fatty Acids, Unsaturated/metabolism , Female , Genome-Wide Association Study , Humans , MaleABSTRACT
Importance: Polygenic risk scores comprising millions of single-nucleotide polymorphisms (SNPs) could be useful for population-wide coronary heart disease (CHD) screening. Objective: To determine whether a polygenic risk score improves prediction of CHD compared with a guideline-recommended clinical risk equation. Design, Setting, and Participants: A retrospective cohort study of the predictive accuracy of a previously validated polygenic risk score was assessed among 4847 adults of white European ancestry, aged 45 through 79 years, participating in the Atherosclerosis Risk in Communities (ARIC) study and 2390 participating in the Multi-Ethnic Study of Atherosclerosis (MESA) from 1996 through December 31, 2015, the final day of follow-up. The performance of the polygenic risk score was compared with that of the 2013 American College of Cardiology and American Heart Association pooled cohort equations. Exposures: Genetic risk was computed for each participant by summing the product of the weights and allele dosage across 6â¯630â¯149 SNPs. Weights were based on an international genome-wide association study. Main Outcomes and Measures: Prediction of 10-year first CHD events (including myocardial infarctions, fatal coronary events, silent infarctions, revascularization procedures, or resuscitated cardiac arrest) assessed using measures of model discrimination, calibration, and net reclassification improvement (NRI). Results: The study population included 4847 adults from the ARIC study (mean [SD] age, 62.9 [5.6] years; 56.4% women) and 2390 adults from the MESA cohort (mean [SD] age, 61.8 [9.6] years; 52.2% women). Incident CHD events occurred in 696 participants (14.4%) and 227 participants (9.5%), respectively, over median follow-up of 15.5 years (interquartile range [IQR], 6.3 years) and 14.2 (IQR, 2.5 years) years. The polygenic risk score was significantly associated with 10-year CHD incidence in ARIC with hazard ratios per SD increment of 1.24 (95% CI, 1.15 to 1.34) and in MESA, 1.38 (95% CI, 1.21 to 1.58). Addition of the polygenic risk score to the pooled cohort equations did not significantly increase the C statistic in either cohort (ARIC, change in C statistic, -0.001; 95% CI, -0.009 to 0.006; MESA, 0.021; 95% CI, -0.0004 to 0.043). At the 10-year risk threshold of 7.5%, the addition of the polygenic risk score to the pooled cohort equations did not provide significant improvement in reclassification in either ARIC (NRI, 0.018, 95% CI, -0.012 to 0.036) or MESA (NRI, 0.001, 95% CI, -0.038 to 0.076). The polygenic risk score did not significantly improve calibration in either cohort. Conclusions and Relevance: In this analysis of 2 cohorts of US adults, the polygenic risk score was associated with incident coronary heart disease events but did not significantly improve discrimination, calibration, or risk reclassification compared with conventional predictors. These findings suggest that a polygenic risk score may not enhance risk prediction in a general, white middle-aged population.
Subject(s)
Coronary Disease/genetics , Genetic Predisposition to Disease , Multifactorial Inheritance , Risk Assessment/methods , Aged , Cohort Studies , Coronary Disease/epidemiology , Female , Genome-Wide Association Study , Genotype , Humans , Incidence , Male , Middle Aged , Myocardial Infarction/epidemiology , Odds Ratio , Phenotype , Polymorphism, Single Nucleotide , Predictive Value of Tests , Proportional Hazards Models , Retrospective Studies , RiskABSTRACT
BACKGROUND: Abiraterone and enzalutamide are recently-approved androgen deprivation therapies (ADTs) for metastatic prostate cancer, with unknown cardiac safety profiles. Abiraterone has a propensity to hypermineralocorticism on top of androgen deprivation, so might carry an additional risk for atrial tachyarrhythmia (AT) and heart failure (HF) compared with other ADTs. AIM: To determine if abiraterone was associated with an increased proportion of AT and HF reports among all suspected adverse drug reactions (ADRs) reported in several pharmacovigilance databases compared with enzalutamide, other ADTs and all other drugs. METHODS: In this observational retrospective pharmacovigilance study, we performed a disproportionality analysis of reports of suspected ADRs in men in the French pharmacovigilance database, the European pharmacovigilance database and the international pharmacovigilance database VigiBase, to evaluate the reporting odds ratios (RORs) of AT and HF for abiraterone compared with enzalutamide, other ADTs and all other drugs. RESULTS: In the 5,759,781 ADR reports in men in VigiBase, 55,070 pertained to ADTs. The RORs for AT for abiraterone versus enzalutamide, other ADTs and all other drugs were 4.1 (95% confidence interval 3.1-5.3), 3.7 (3-4.5) and 3.2 (2.7-3.7), respectively (P<0.0001 for all). The corresponding RORs for HF were 2.5 (2-3), 1.5 (1.3-1.7) and 2 (1.7-2.3), respectively (P<0.0001 for all). These results were concordant with the French and European pharmacovigilance databases. Mean times to AT and HF onset were shorter with abiraterone (5.2±0.8 and 4.5±0.6 months, respectively) versus other ADTs (13.3±3.2 and 9.2±1.1 months, respectively) (both P<0.05). Cases on abiraterone versus other ADTs were more frequently associated with at least two ADR terms, including AT, HF, hypokalaemia, hypertension and oedema (13.6% vs 6%; P<0.0001). For abiraterone, age, but not dose, was associated with reporting of AT and HF versus any other ADR. CONCLUSIONS: Compared with other ADTs, abiraterone was associated with higher reporting of AT and HF, associated with hypokalaemia, hypertension and oedema. These findings are consistent with the hypermineralocorticism induced by abiraterone, but not by other ADTs.
Subject(s)
Androgen Antagonists/adverse effects , Androstenes/adverse effects , Antineoplastic Agents, Hormonal/adverse effects , Heart Failure/chemically induced , Pharmacovigilance , Phenylthiohydantoin/analogs & derivatives , Prostatic Neoplasms/drug therapy , Tachycardia, Supraventricular/chemically induced , Adverse Drug Reaction Reporting Systems , Aged , Aged, 80 and over , Benzamides , Cardiotoxicity , Databases, Factual , Heart Failure/diagnosis , Heart Failure/epidemiology , Humans , Male , Middle Aged , Nitriles , Phenylthiohydantoin/adverse effects , Prostatic Neoplasms/diagnosis , Prostatic Neoplasms/epidemiology , Retrospective Studies , Risk Assessment , Risk Factors , Tachycardia, Supraventricular/diagnosis , Tachycardia, Supraventricular/epidemiology , Time FactorsABSTRACT
Importance: Whether the PCSK9 gene is associated with the progress from infection to sepsis is unknown to date. Objective: To test the associations between PCSK9 genetic variants, a PCSK9 genetic risk score (GRS), or genetically estimated PCSK9 expression levels and the risk of sepsis among patients admitted to a hospital with infection. Design, Setting, and Participants: This retrospective cohort study used deidentified electronic health records to identify patients admitted to Vanderbilt University Medical Center, Nashville, Tennessee, with infection. Patients were white adults, had a code indicating infection from the International Classification of Diseases, Ninth Revision, Clinical Modification, or the International Statistical Classification of Diseases, Tenth Revision, Clinical Modification, and received an antibiotic within 1 day of hospital admission (N = 61 502). Data were collected from January 1, 1993, through December 31, 2017, and analyzed from April 1, 2018, to March 16, 2019. Exposures: Four known PCSK9 functional variants, a GRS for PCSK9, and genetically estimated PCSK9 expression. Main Outcomes and Measures: The primary outcome was sepsis; secondary outcomes included cardiovascular failure and in-hospital death. Results: Of patients with infection, genotype information was available in 10 922 white patients for PCSK9 functional variants (5628 men [51.5%]; mean [SD] age, 60.1 [15.7] years), including 7624 patients with PCSK9 GRS and 6033 patients with estimated PCSK9 expression. Of these, 3391 developed sepsis, 835 developed cardiovascular failure, and 366 died during hospitalization. None of the 4 functional PCSK9 variants were significantly associated with sepsis, cardiovascular failure, or in-hospital death, with or without adjustment for (1) age and sex or (2) age, sex, and Charlson-Deyo comorbidities (in model adjusted for age, sex, and comorbidities, odds ratios for any loss-of function variant were 0.96 [95% CI, 0.88-1.04] for sepsis, 1.05 [95% CI, 0.90-1.22] for cardiovascular failure, and 0.89 [95% CI, 0.72-1.11] for death). Similarly, neither the PCSK9 GRS nor genetically estimated PCSK9 expression were significantly associated with sepsis, cardiovascular failure, or in-hospital death in any of the analysis models. For GRS, in the full model adjusted for age, sex, and comorbidities, the odds ratios were 1.01 for sepsis (95% CI, 0.96-1.06; P = .70), 1.03 for cardiovascular failure (95% CI, 0.95-1.12; P = .48), and 1.05 for in-hospital death (95% CI, 0.92-1.19; P = .50). For genetically estimated PCSK9 expression, in the full model adjusted for age, sex, and comorbidities, the odds ratios were 1.01 for sepsis (95% CI, 0.95-1.06; P = .86), 0.96 for cardiovascular failure (95% CI, 0.88-1.05; P = .41), and 0.99 for in-hospital death (95% CI, 0.87-1.14; P = .94). Conclusions and Relevance: In this study, PCSK9 genetic variants were not significantly associated with risk of sepsis or the outcomes of sepsis in patients hospitalized with infection.