An ECG-based machine learning model for predicting new-onset atrial fibrillation is superior to age and clinical features in identifying patients at high stroke risk

BACKGROUND

Several large trials have employed age or clinical features to select patients for atrial fibrillation (AF) screening to reduce strokes. We hypothesized that a machine learning (ML) model trained to predict AF risk from 12‑lead electrocardiogram (ECG) would be more efficient than criteria based on clinical variables in indicating a population for AF screening to potentially prevent AF-related stroke.

METHODS

We retrospectively included all patients with clinical encounters in Geisinger without a prior history of AF. Incidence of AF within 1 year and AF-related strokes within 3 years of the encounter were identified. AF-related stroke was defined as a stroke where AF was diagnosed at the time of stroke or within a year after the stroke. The efficiency of five methods was evaluated for selecting a cohort for AF screening. The methods were selected from four clinical trials (mSToPS, GUARD-AF, SCREEN-AF and STROKESTOP) and the ECG-based ML model. We simulated patient selection for the five methods between the years 2011 and 2014 and evaluated outcomes for 1 year intervals between 2012 and 2015, resulting in a total of twenty 1-year periods. Patients were considered eligible if they met the criteria before the start of the given 1-year period or within that period. The primary outcomes were numbers needed to screen (NNS) for AF and AF-associated stroke.

RESULTS

The clinical trial models indicated large proportions of the population with a prior ECG for AF screening (up to 31%), coinciding with NNS ranging from 14 to 18 for AF and 249-359 for AF-associated stroke. At comparable sensitivity, the ECG ML model indicated a modest number of patients for screening (14%) and had the highest efficiency in NNS for AF (7.3; up to 60% reduction) and AF-associated stroke (223; up to 38% reduction).

CONCLUSIONS

An ECG-based ML risk prediction model is more efficient than contemporary AF-screening criteria based on age alone or age and clinical features at indicating a population for AF screening to potentially prevent AF-related strokes.

Impact of a Population Genomic Screening Program on Health Behaviors Related to Familial Hypercholesterolemia Risk Reduction

Limited information is available regarding clinician and participant behaviors after disclosure of genomic risk variants for familial hypercholesterolemia (FH) from a population genomic screening program.

Predictive Accuracy of a Clinical and Genetic Risk Model for Atrial Fibrillation

BACKGROUND

Atrial fibrillation (AF) risk estimation using clinical factors with or without genetic information may identify AF screening candidates more accurately than the guideline-based age threshold of ≥65 years.

METHODS

We analyzed 4 samples across the United States and Europe (derivation: UK Biobank; validation: FINRISK, Geisinger MyCode Initiative, and Framingham Heart Study). We estimated AF risk using the CHARGE-AF (Cohorts for Heart and Aging Research in Genomic Epidemiology AF) score and a combination of CHARGE-AF and a 1168-variant polygenic score (Predict-AF). We compared the utility of age, CHARGE-AF, and Predict-AF for predicting 5-year AF by quantifying discrimination and calibration.

RESULTS

Among 543 093 individuals, 8940 developed AF within 5 years. In the validation sets, CHARGE-AF (C index range, 0.720-0.824) and Predict-AF (0.749-0.831) had largely comparable discrimination, both favorable to continuous age (0.675-0.801). Calibration was similar using CHARGE-AF (slope range, 0.67-0.87) and Predict-AF (0.65-0.83). Net reclassification improvement using Predict-AF versus CHARGE-AF was modest (net reclassification improvement range, 0.024-0.057) but more favorable among individuals aged <65 years (0.062-0.11). Using Predict-AF among 99 530 individuals aged ≥65 years across each sample, 70 849 had AF risk <5%, of whom 69 067 (97.5%) did not develop AF, whereas 28 681 had AF risk ≥5%, of whom 2264 (7.9%) developed AF. Of 11 379 individuals aged <65 years with AF risk ≥5%, 435 (3.8%) developed AF before age 65 years, with roughly half (46.9%) meeting anticoagulation criteria.

CONCLUSIONS

AF risk estimation using clinical factors may prioritize individuals for AF screening more precisely than the age threshold endorsed in current guidelines. The additional value of genetic predisposition is modest but greatest among younger individuals.

Deep-learning-assisted analysis of echocardiographic videos improves predictions of all-cause mortality

Machine learning promises to assist physicians with predictions of mortality and of other future clinical events by learning complex patterns from historical data, such as longitudinal electronic health records. Here we show that a convolutional neural network trained on raw pixel data in 812,278 echocardiographic videos from 34,362 individuals provides superior predictions of one-year all-cause mortality. The model's predictions outperformed the widely used pooled cohort equations, the Seattle Heart Failure score (measured in an independent dataset of 2,404 patients with heart failure who underwent 3,384 echocardiograms), and a machine learning model involving 58 human-derived variables from echocardiograms and 100 clinical variables derived from electronic health records. We also show that cardiologists assisted by the model substantially improved the sensitivity of their predictions of one-year all-cause mortality by 13% while maintaining prediction specificity. Large unstructured datasets may enable deep learning to improve a wide range of clinical prediction models.

Routinely reported ejection fraction and mortality in clinical practice: where does the nadir of risk lie?

AIMS

We investigated the relationship between clinically assessed left ventricular ejection fraction (LVEF) and survival in a large, heterogeneous clinical cohort.

METHODS AND RESULTS

Physician-reported LVEF on 403 977 echocardiograms from 203 135 patients were linked to all-cause mortality using electronic health records (1998-2018) from US regional healthcare system. Cox proportional hazards regression was used for analyses while adjusting for many patient characteristics including age, sex, and relevant comorbidities. A dataset including 45 531 echocardiograms and 35 976 patients from New Zealand was used to provide independent validation of analyses. During follow-up of the US cohort, 46 258 (23%) patients who had undergone 108 578 (27%) echocardiograms died. Overall, adjusted hazard ratios (HR) for mortality showed a u-shaped relationship for LVEF with a nadir of risk at an LVEF of 60-65%, a HR of 1.71 [95% confidence interval (CI) 1.64-1.77] when ≥70% and a HR of 1.73 (95% CI 1.66-1.80) at LVEF of 35-40%. Similar relationships with a nadir at 60-65% were observed in the validation dataset as well as for each age group and both sexes. The results were similar after further adjustments for conditions associated with an elevated LVEF, including mitral regurgitation, increased wall thickness, and anaemia and when restricted to patients reported to have heart failure at the time of the echocardiogram.

CONCLUSION

Deviation of LVEF from 60% to 65% is associated with poorer survival regardless of age, sex, or other relevant comorbidities such as heart failure. These results may herald the recognition of a new phenotype characterized by supra-normal LVEF.

Deep Neural Networks Can Predict New-Onset Atrial Fibrillation From the 12-Lead ECG and Help Identify Those at Risk of Atrial Fibrillation-Related Stroke

Supplemental Digital Content is available in the text.

Atrial fibrillation (AF) is associated with substantial morbidity, especially when it goes undetected. If new-onset AF could be predicted, targeted screening could be used to find it early. We hypothesized that a deep neural network could predict new-onset AF from the resting 12-lead ECG and that this prediction may help identify those at risk of AF-related stroke.

Electronic health record analysis identifies kidney disease as the leading risk factor for hospitalization in confirmed COVID-19 patients

BACKGROUND

Empirical data on conditions that increase risk of coronavirus disease 2019 (COVID-19) progression are needed to identify high risk individuals. We performed a comprehensive quantitative assessment of pre-existing clinical phenotypes associated with COVID-19-related hospitalization.

METHODS

Phenome-wide association study (PheWAS) of SARS-CoV-2-positive patients from an integrated health system (Geisinger) with system-level outpatient/inpatient COVID-19 testing capacity and retrospective electronic health record (EHR) data to assess pre-COVID-19 pandemic clinical phenotypes associated with hospital admission (hospitalization).

RESULTS

Of 12,971 individuals tested for SARS-CoV-2 with sufficient pre-COVID-19 pandemic EHR data at Geisinger, 1604 were SARS-CoV-2 positive and 354 required hospitalization. We identified 21 clinical phenotypes in 5 disease categories meeting phenome-wide significance (P<1.60x10-4), including: six kidney phenotypes, e.g. end stage renal disease or stage 5 CKD (OR = 11.07, p = 1.96x10-8), six cardiovascular phenotypes, e.g. congestive heart failure (OR = 3.8, p = 3.24x10-5), five respiratory phenotypes, e.g. chronic airway obstruction (OR = 2.54, p = 3.71x10-5), and three metabolic phenotypes, e.g. type 2 diabetes (OR = 1.80, p = 7.51x10-5). Additional analyses defining CKD based on estimated glomerular filtration rate, confirmed high risk of hospitalization associated with pre-existing stage 4 CKD (OR 2.90, 95% CI: 1.47, 5.74), stage 5 CKD/dialysis (OR 8.83, 95% CI: 2.76, 28.27), and kidney transplant (OR 14.98, 95% CI: 2.77, 80.8) but not stage 3 CKD (OR 1.03, 95% CI: 0.71, 1.48).

CONCLUSIONS

This study provides quantitative estimates of the contribution of pre-existing clinical phenotypes to COVID-19 hospitalization and highlights kidney disorders as the strongest factors associated with hospitalization in an integrated US healthcare system.

A Machine Learning Approach to Management of Heart Failure Populations

BACKGROUND

Heart failure is a prevalent, costly disease for which new value-based payment models demand optimized population management strategies.

OBJECTIVES

This study sought to generate a strategy for managing populations of patients with heart failure by leveraging large clinical datasets and machine learning.

METHODS

Geisinger electronic health record data were used to train machine learning models to predict 1-year all-cause mortality in 26,971 patients with heart failure who underwent 276,819 clinical episodes. There were 26 clinical variables (demographics, laboratory test results, medications), 90 diagnostic codes, 41 electrocardiogram measurements and patterns, 44 echocardiographic measurements, and 8 evidence-based "care gaps": flu vaccine, blood pressure of <130/80 mm Hg, A_1c of <8%, cardiac resynchronization therapy, and active medications (active angiotensin-converting enzyme inhibitor/angiotensin II receptor blocker/angiotensin receptor-neprilysin inhibitor, aldosterone receptor antagonist, hydralazine, and evidence-based beta-blocker) were collected. Care gaps represented actionable variables for which associations with all-cause mortality were modeled from retrospective data and then used to predict the benefit of prospective interventions in 13,238 currently living patients.

RESULTS

Machine learning models achieved areas under the receiver-operating characteristic curve (AUCs) of 0.74 to 0.77 in a split-by-year training/test scheme, with the nonlinear XGBoost model (AUC: 0.77) outperforming linear logistic regression (AUC: 0.74). Out of 13,238 currently living patients, 2,844 were predicted to die within a year, and closing all care gaps was predicted to save 231 of these lives. Prioritizing patients for intervention by using the predicted reduction in 1-year mortality risk outperformed all other priority rankings (e.g., random selection or Seattle Heart Failure risk score).

CONCLUSIONS

Machine learning can be used to priority-rank patients most likely to benefit from interventions to optimize evidence-based therapies. This approach may prove useful for optimizing heart failure population health management teams within value-based payment models.

Healthcare Utilization and Costs after Receiving a Positive BRCA1/2 Result from a Genomic Screening Program

Population genomic screening has been demonstrated to detect at-risk individuals who would not be clinically identified otherwise. However, there are concerns about the increased utilization of unnecessary services and the associated increase in costs. The objectives of this study are twofold: (1) determine whether there is a difference in healthcare utilization and costs following disclosure of a pathogenic/likely pathogenic (P/LP) BRCA1/2 variant via a genomic screening program, and (2) measure the post-disclosure uptake of National Comprehensive Cancer Network (NCCN) guideline-recommended risk management. We retrospectively reviewed electronic health record (EHR) and billing data from a female population of BRCA1/2 P/LP variant carriers without a personal history of breast or ovarian cancer enrolled in Geisinger’s MyCode genomic screening program with at least a one-year post-disclosure observation period. We identified 59 women for the study cohort out of 50,726 MyCode participants. We found no statistically significant differences in inpatient and outpatient utilization and average total costs between one-year pre- and one-year post-disclosure periods ($18,821 vs. $19,359, p = 0.76). During the first year post-disclosure, 49.2% of women had a genetic counseling visit, 45.8% had a mammography and 32.2% had an MRI. The uptake of mastectomy and oophorectomy was 3.5% and 11.8%, respectively, and 5% of patients received chemoprevention.

Prevalence and Electronic Health Record-Based Phenotype of Loss-of-Function Genetic Variants in Arrhythmogenic Right Ventricular Cardiomyopathy-Associated Genes

BACKGROUND

Arrhythmogenic right ventricular cardiomyopathy (ARVC) is associated with variants in desmosome genes. Secondary findings of pathogenic/likely pathogenic variants, primarily loss-of-function (LOF) variants, are recommended for clinical reporting; however, their prevalence and associated phenotype in a general clinical population are not fully characterized.

METHODS

From whole-exome sequencing of 61 019 individuals in the DiscovEHR cohort, we screened for putative loss-of-function variants in PKP2, DSC2, DSG2, and DSP. We evaluated measures from prior clinical ECG and echocardiograms, manually over-read to evaluate ARVC diagnostic criteria, and performed a PheWAS (phenome-wide association study). Finally, we estimated expected penetrance using Bayesian inference.

RESULTS

One hundred forty individuals (0.23%; 59±18 years old at last encounter; 33% male) had an ARVC variant (G⁺). None had an existing diagnosis of ARVC in the electronic health record, nor significant differences in prior ECG or echocardiogram findings compared with matched controls without variants. Several G⁺ individuals satisfied major repolarization (n=4) and ventricular function (n=5) criteria, but this prevalence matched controls. PheWAS showed no significant associations of other heart disease diagnoses. Combining our best genetic and disease prevalence estimates yields an estimated penetrance of 6.0%.

CONCLUSIONS

The prevalence of ARVC loss-of-function variants is ≈1:435 in a general clinical population of predominantly European descent, but with limited electronic health record-based evidence of phenotypic association in our population, consistent with a low penetrance estimate. Prospective deep phenotyping and longitudinal follow-up of a large sequenced cohort is needed to determine the true clinical relevance of an incidentally identified ARVC loss-of-function variant.

GSTM1 Copy Number Is Not Associated With Risk of Kidney Failure in a Large Cohort

Deletion of glutathione S-transferase µ1 (GSTM1) is common in populations and has been asserted to associate with chronic kidney disease progression in some research studies. The association needs to be validated. We estimated GSTM1 copy number using whole exome sequencing data in the DiscovEHR cohort. Kidney failure was defined as requiring dialysis or receiving kidney transplant using data from the electronic health record and linkage to the United States Renal Data System, or the most recent eGFR < 15 ml/min/1.73 m². In a cohort of 46,983 unrelated participants, 28.8% of blacks and 52.1% of whites had 0 copies of GSTM1. Over a mean of 9.2 years follow-up, 645 kidney failure events were observed in 46,187 white participants, and 28 in 796 black participants. No significant association was observed between GSTM1 copy number and kidney failure in Cox regression adjusting for age, sex, BMI, smoking status, genetic principal components, or comorbid conditions (hypertension, diabetes, heart failure, coronary artery disease, and stroke), whether using a genotypic, dominant, or recessive model. In sensitivity analyses, GSTM1 copy number was not associated with kidney failure in participants that were 45 years or older at baseline, had baseline eGFR < 60 ml/min/1.73 m², or with baseline year between 1996 and 2002. In conclusion, we found no association between GSTM1 copy number and kidney failure in a large cohort study.

ALG9 Mutation Carriers Develop Kidney and Liver Cysts

BACKGROUND

Mutations in PKD1 or PKD2 cause typical autosomal dominant polycystic kidney disease (ADPKD), the most common monogenic kidney disease. Dominantly inherited polycystic kidney and liver diseases on the ADPKD spectrum are also caused by mutations in at least six other genes required for protein biogenesis in the endoplasmic reticulum, the loss of which results in defective production of the PKD1 gene product, the membrane protein polycystin-1 (PC1).

METHODS

We used whole-exome sequencing in a cohort of 122 patients with genetically unresolved clinical diagnosis of ADPKD or polycystic liver disease to identify a candidate gene, ALG9, and in vitro cell-based assays of PC1 protein maturation to functionally validate it. For further validation, we identified carriers of ALG9 loss-of-function mutations and noncarrier matched controls in a large exome-sequenced population-based cohort and evaluated the occurrence of polycystic phenotypes in both groups.

RESULTS

Two patients in the clinically defined cohort had rare loss-of-function variants in ALG9, which encodes a protein required for addition of specific mannose molecules to the assembling N-glycan precursors in the endoplasmic reticulum lumen. In vitro assays showed that inactivation of Alg9 results in impaired maturation and defective glycosylation of PC1. Seven of the eight (88%) cases selected from the population-based cohort based on ALG9 mutation carrier state who had abdominal imaging after age 50; seven (88%) had at least four kidney cysts, compared with none in matched controls without ALG9 mutations.

CONCLUSIONS

ALG9 is a novel disease gene in the genetically heterogeneous ADPKD spectrum. This study supports the utility of phenotype characterization in genetically-defined cohorts to validate novel disease genes, and provide much-needed genotype-phenotype correlations.

A genome-first approach to aggregating rare genetic variants in LMNA for association with electronic health record phenotypes

PURPOSE

"Genome-first" approaches, in which genetic sequencing is agnostically linked to associated phenotypes, can enhance our understanding of rare variants' contributions to disease. Loss-of-function variants in LMNA cause a range of rare diseases, including cardiomyopathy.

METHODS

We leveraged exome sequencing from 11,451 unselected individuals in the Penn Medicine Biobank to associate rare variants in LMNA with diverse electronic health record (EHR)-derived phenotypes. We used Rare Exome Variant Ensemble Learner (REVEL) to annotate rare missense variants, clustered predicted deleterious and loss-of-function variants into a "gene burden" (N = 72 individuals), and performed a phenome-wide association study (PheWAS). Major findings were replicated in DiscovEHR.

RESULTS

The LMNA gene burden was significantly associated with primary cardiomyopathy (p = 1.78E-11) and cardiac conduction disorders (p = 5.27E-07). Most patients had not been clinically diagnosed with LMNA cardiomyopathy. We also noted an association with chronic kidney disease (p = 1.13E-06). Regression analyses on echocardiography and serum labs revealed that LMNA variant carriers had dilated cardiomyopathy and primary renal disease.

CONCLUSION

Pathogenic LMNA variants are an underdiagnosed cause of cardiomyopathy. We also find that LMNA loss of function may be a primary cause of renal disease. Finally, we show the value of aggregating rare, annotated variants into a gene burden and using PheWAS to identify novel ontologies for pleiotropic human genes.

Genomics-First Evaluation of Heart Disease Associated With Titin-Truncating Variants

BACKGROUND

Truncating variants in the Titin gene (TTNtvs) are common in individuals with idiopathic dilated cardiomyopathy (DCM). However, a comprehensive genomics-first evaluation of the impact of TTNtvs in different clinical contexts, and the evaluation of modifiers such as genetic ancestry, has not been performed.

METHODS

We reviewed whole exome sequence data for >71 000 individuals (61 040 from the Geisinger MyCode Community Health Initiative (2007 to present) and 10 273 from the PennMedicine BioBank (2013 to present) to identify anyone with TTNtvs. We further selected individuals with TTNtvs in exons highly expressed in the heart (proportion spliced in [PSI] >0.9). Using linked electronic health records, we evaluated associations of TTNtvs with diagnoses and quantitative echocardiographic measures, including subanalyses for individuals with and without DCM diagnoses. We also reviewed data from the Jackson Heart Study to validate specific analyses for individuals of African ancestry.

RESULTS

Identified with a TTNtv in a highly expressed exon (hiPSI) were 1.2% individuals in PennMedicine BioBank and 0.6% at Geisinger. The presence of a hiPSI TTNtv was associated with increased odds of DCM in individuals of European ancestry (odds ratio [95% CI]: 18.7 [9.1-39.4] {PennMedicine BioBank} and 10.8 [7.0-16.0] {Geisinger}). hiPSI TTNtvs were not associated with DCM in individuals of African ancestry, despite a high DCM prevalence (odds ratio, 1.8 [0.2-13.7]; P=0.57). Among 244 individuals of European ancestry with DCM in PennMedicine BioBank, hiPSI TTNtv carriers had lower left ventricular ejection fraction (β=-12%, P=3×10^-7), and increased left ventricular diameter (β=0.65 cm, P=9×10^-3). In the Geisinger cohort, hiPSI TTNtv carriers without a cardiomyopathy diagnosis had more atrial fibrillation (odds ratio, 2.4 [1.6-3.6]) and heart failure (odds ratio, 3.8 [2.4-6.0]), and lower left ventricular ejection fraction (β=-3.4%, P=1×10^-7).

CONCLUSIONS

Individuals of European ancestry with hiPSI TTNtv have an abnormal cardiac phenotype characterized by lower left ventricular ejection fraction, irrespective of the clinical manifestation of cardiomyopathy. Associations with arrhythmias, including atrial fibrillation, were observed even when controlling for cardiomyopathy diagnosis. In contrast, no association between hiPSI TTNtvs and DCM was discerned among individuals of African ancestry. Given these findings, clinical identification of hiPSI TTNtv carriers may alter clinical management strategies.

Exome Sequencing-Based Screening for BRCA1/2 Expected Pathogenic Variants Among Adult Biobank Participants

IMPORTANCE

Detection of disease-associated variants in the BRCA1 and BRCA2 (BRCA1/2) genes allows for cancer prevention and early diagnosis in high-risk individuals.

OBJECTIVES

To identify pathogenic and likely pathogenic (P/LP) BRCA1/2 variants in an unselected research cohort, and to characterize the features associated with P/LP variants.

DESIGN, SETTING, AND PARTICIPANTS

This is a cross-sectional study of adult volunteers (n = 50 726) who underwent exome sequencing at a single health care system (Geisinger Health System, Danville, Pennsylvania) from January 1, 2014, to March 1, 2016. Participants are part of the DiscovEHR cohort and were identified through the Geisinger MyCode Community Health Initiative. They consented to a research protocol that included sequencing and return of actionable test results. Clinical data from electronic health records and clinical visits were correlated with variants. Comparisons were made between those with (cases) and those without (controls) P/LP variants in BRCA1/2.

MAIN OUTCOMES

Prevalence of P/LP BRCA1/2 variants in cohort, proportion of variant carriers not previously ascertained through clinical testing, and personal and family history of relevant cancers among BRCA1/2 variant carriers and noncarriers.

RESULTS

Of the 50 726 health system patients who underwent exome sequencing, 50 459 (99.5%) had no expected pathogenic BRCA1/2 variants and 267 (0.5%) were BRCA1/2 carriers. Of the 267 cases (148 [55.4%] were women and 119 [44.6%] were men with a mean [range] age of 58.9 [23-90] years), 183 (68.5%) received clinically confirmed results in their electronic health record. Among the 267 participants with P/LP BRCA1/2 variants, 219 (82.0%) had no prior clinical testing, 95 (35.6%) had BRCA1 variants, and 172 (64.4%) had BRCA2 variants. Syndromic cancer diagnoses were present in 11 (47.8%) of the 23 deceased BRCA1/2 carriers and in 56 (20.9%) of all 267 BRCA1/2 carriers. Among women, 31 (20.9%) of 148 variant carriers had a personal history of breast cancer, compared with 1554 (5.2%) of 29 880 noncarriers (odds ratio [OR], 5.95; 95% CI, 3.88-9.13; P < .001). Ovarian cancer history was present in 15 (10.1%) of 148 variant carriers and in 195 (0.6%) of 29 880 variant noncarriers (OR, 18.30; 95% CI, 10.48-31.4; P < .001). Among 89 BRCA1/2 carriers without prior testing but with comprehensive personal and family history data, 44 (49.4%) did not meet published guidelines for clinical testing.

CONCLUSIONS AND RELEVANCE

This study found that compared with previous clinical care, exome sequencing-based screening identified 5 times as many individuals with P/LP BRCA1/2 variants. These findings suggest that genomic screening may identify BRCA1/2-associated cancer risk that might otherwise remain undetected within health care systems and may provide opportunities to reduce morbidity and mortality in patients.

Functional Invalidation of Putative Sudden Infant Death Syndrome-Associated Variants in the KCNH2-Encoded Kv11.1 Channel

BACKGROUND

Heterologous functional validation studies of putative long-QT syndrome subtype 2-associated variants clarify their pathological potential and identify disease mechanism(s) for most variants studied. The purpose of this study is to clarify the pathological potential for rare nonsynonymous KCNH2 variants seemingly associated with sudden infant death syndrome.

METHODS

Genetic testing of 292 sudden infant death syndrome cases identified 9 KCNH2 variants: E90K, R181Q, A190T, G294V, R791W, P967L, R1005W, R1047L, and Q1068R. Previous studies show R181Q-, P967L-, and R1047L-Kv11.1 channels function similar to wild-type Kv11.1 channels, whereas Q1068R-Kv11.1 channels accelerate inactivation gating. We studied the biochemical and biophysical properties for E90K-, G294V-, R791W-, and R1005W-Kv11.1 channels expressed in human embryonic kidney 293 cells; examined the electronic health records of patients who were genotype positive for the sudden infant death syndrome-linked KCNH2 variants; and simulated their functional impact using computational models of the human ventricular action potential.

RESULTS

Western blot and voltage-clamping analyses of cells expressing E90K-, G294V-, R791W-, and R1005W-Kv11.1 channels demonstrated these variants express and generate peak Kv11.1 current levels similar to cells expressing wild-type-Kv11.1 channels, but R791W- and R1005W-Kv11.1 channels accelerated deactivation and activation gating, respectively. Electronic health records of patients with the sudden infant death syndrome-linked KCNH2 variants showed that the patients had median heart rate-corrected QT intervals <480 ms and none had been diagnosed with long-QT syndrome or experienced cardiac arrest. Simulating the impact of dysfunctional gating variants predicted that they have little impact on ventricular action potential duration.

CONCLUSIONS

We conclude that these rare Kv11.1 missense variants are not long-QT syndrome subtype 2-causative variants and therefore do not represent the pathogenic substrate for sudden infant death syndrome in the variant-positive infants.

PheWAS and Beyond: The Landscape of Associations with Medical Diagnoses and Clinical Measures across 38,662 Individuals from Geisinger

Most phenome-wide association studies (PheWASs) to date have used a small to moderate number of SNPs for association with phenotypic data. We performed a large-scale single-cohort PheWAS, using electronic health record (EHR)-derived case-control status for 541 diagnoses using International Classification of Disease version 9 (ICD-9) codes and 25 median clinical laboratory measures. We calculated associations between these diagnoses and traits with ∼630,000 common frequency SNPs with minor allele frequency > 0.01 for 38,662 individuals. In this landscape PheWAS, we explored results within diseases and traits, comparing results to those previously reported in genome-wide association studies (GWASs), as well as previously published PheWASs. We further leveraged the context of functional impact from protein-coding to regulatory regions, providing a deeper interpretation of these associations. The comprehensive nature of this PheWAS allows for novel hypothesis generation, the identification of phenotypes for further study for future phenotypic algorithm development, and identification of cross-phenotype associations.

Rare variants in drug target genes contributing to complex diseases, phenome-wide

The DrugBank database consists of ~800 genes that are well characterized drug targets. This list of genes is a useful resource for association testing. For example, loss of function (LOF) genetic variation has the potential to mimic the effect of drugs, and high impact variation in these genes can impact downstream traits. Identifying novel associations between genetic variation in these genes and a range of diseases can also uncover new uses for the drugs that target these genes. Phenome Wide Association Studies (PheWAS) have been successful in identifying genetic associations across hundreds of thousands of diseases. We have conducted a novel gene based PheWAS to test the effect of rare variants in DrugBank genes, evaluating associations between these genes and more than 500 quantitative and dichotomous phenotypes. We used whole exome sequencing data from 38,568 samples in Geisinger MyCode Community Health Initiative. We evaluated the results of this study when binning rare variants using various filters based on potential functional impact. We identified multiple novel associations, and the majority of the significant associations were driven by functionally annotated variation. Overall, this study provides a sweeping exploration of rare variant associations within functionally relevant genes across a wide range of diagnoses.

Distribution and clinical impact of functional variants in 50,726 whole-exome sequences from the DiscovEHR study

The DiscovEHR collaboration between the Regeneron Genetics Center and Geisinger Health System couples high-throughput sequencing to an integrated health care system using longitudinal electronic health records (EHRs). We sequenced the exomes of 50,726 adult participants in the DiscovEHR study to identify ~4.2 million rare single-nucleotide variants and insertion/deletion events, of which ~176,000 are predicted to result in a loss of gene function. Linking these data to EHR-derived clinical phenotypes, we find clinical associations supporting therapeutic targets, including genes encoding drug targets for lipid lowering, and identify previously unidentified rare alleles associated with lipid levels and other blood level traits. About 3.5% of individuals harbor deleterious variants in 76 clinically actionable genes. The DiscovEHR data set provides a blueprint for large-scale precision medicine initiatives and genomics-guided therapeutic discovery.

Electronic health record phenotype in subjects with genetic variants associated with arrhythmogenic right ventricular cardiomyopathy: a study of 30,716 subjects with exome sequencing

Purpose

Arrhythmogenic right ventricular cardiomyopathy (ARVC) is an inherited heart disease. Clinical follow-up of incidental findings in ARVC-associated genes is recommended. We aimed to determine the prevalence of disease thus ascertained.

Methods

30,716 individuals underwent exome sequencing. Variants in PKP2, DSG2, DSC2, DSP, JUP, TMEM43, or TGFβ3 that were database-listed as pathogenic or likely pathogenic were identified and evidence-reviewed. For subjects with putative loss-of-function (pLOF) variants or variants of uncertain significance (VUS), electronic health records (EHR) were reviewed for ARVC diagnosis, diagnostic criteria, and International Classification of Diseases (ICD-9) codes.

Results

18 subjects had pLOF variants; none had an EHR diagnosis of ARVC. Of 14 patients with an electrocardiogram (ECG), one had a minor diagnostic criterion, 13 were normal. 184 subjects had VUSs; none had an ARVC diagnosis. In subjects with VUSs, there was no difference in the proportion with major (4%) or minor (13%) ECG diagnostic criteria compared to variant-negative controls. ICD-9 codes showed no difference in defibrillator utilization, electrophysiologic abnormalities or non-ischemic cardiomyopathies in patients with pLOF or VUSs compared to controls.

Conclusion

pLOF variants in an unselected cohort were not associated with ARVC phenotypes based on EHR review. The negative predictive value of EHR review remains uncertain.

Association of Rare and Common Variation in the Lipoprotein Lipase Gene With Coronary Artery Disease

IMPORTANCE

The activity of lipoprotein lipase (LPL) is the rate-determining step in clearing triglyceride-rich lipoproteins from the circulation. Mutations that damage the LPL gene (LPL) lead to lifelong deficiency in enzymatic activity and can provide insight into the relationship of LPL to human disease.

OBJECTIVE

To determine whether rare and/or common variants in LPL are associated with early-onset coronary artery disease (CAD).

DESIGN, SETTING, AND PARTICIPANTS

In a cross-sectional study, LPL was sequenced in 10 CAD case-control cohorts of the multinational Myocardial Infarction Genetics Consortium and a nested CAD case-control cohort of the Geisinger Health System DiscovEHR cohort between 2010 and 2015. Common variants were genotyped in up to 305 699 individuals of the Global Lipids Genetics Consortium and up to 120 600 individuals of the CARDIoGRAM Exome Consortium between 2012 and 2014. Study-specific estimates were pooled via meta-analysis.

EXPOSURES

Rare damaging mutations in LPL included loss-of-function variants and missense variants annotated as pathogenic in a human genetics database or predicted to be damaging by computer prediction algorithms trained to identify mutations that impair protein function. Common variants in the LPL gene region included those independently associated with circulating triglyceride levels.

MAIN OUTCOMES AND MEASURES

Circulating lipid levels and CAD.

RESULTS

Among 46 891 individuals with LPL gene sequencing data available, the mean (SD) age was 50 (12.6) years and 51% were female. A total of 188 participants (0.40%; 95% CI, 0.35%-0.46%) carried a damaging mutation in LPL, including 105 of 32 646 control participants (0.32%) and 83 of 14 245 participants with early-onset CAD (0.58%). Compared with 46 703 noncarriers, the 188 heterozygous carriers of an LPL damaging mutation displayed higher plasma triglyceride levels (19.6 mg/dL; 95% CI, 4.6-34.6 mg/dL) and higher odds of CAD (odds ratio = 1.84; 95% CI, 1.35-2.51; P < .001). An analysis of 6 common LPL variants resulted in an odds ratio for CAD of 1.51 (95% CI, 1.39-1.64; P = 1.1 × 10-22) per 1-SD increase in triglycerides.

CONCLUSIONS AND RELEVANCE

The presence of rare damaging mutations in LPL was significantly associated with higher triglyceride levels and presence of coronary artery disease. However, further research is needed to assess whether there are causal mechanisms by which heterozygous lipoprotein lipase deficiency could lead to coronary artery disease.

Contrasting Association Results between Existing PheWAS Phenotype Definition Methods and Five Validated Electronic Phenotypes

Phenome-Wide Association Studies (PheWAS) comprehensively investigate the association between genetic variation and a wide array of outcome traits. Electronic health record (EHR) based PheWAS uses various abstractions of International Classification of Diseases, Ninth Revision (ICD-9) codes to identify case/control status for diagnoses that are used as the phenotypic variables. However, there have not been comparisons within a PheWAS between results from high quality derived phenotypes and high-throughput but potentially inaccurate use of ICD-9 codes for case/control definition. For this study we first developed a group of high quality algorithms for five phenotypes. Next we evaluated the association of these "gold standard" phenotypes and 4,636,178 genetic variants with minor allele frequency > 0.01 and compared the results from high-throughput associations at the 3 digit, 5 digit, and PheWAS codes for defining case/control status. We found that certain diseases contained similar patient populations across phenotyping methods but had differences in PheWAS.