Search

Angiopoietin-like 4 (ANGPTL4) is an endogenous inhibitor of lipoprotein lipase that modulates lipid levels, coronary atherosclerosis risk, and nutrient partitioning. We hypothesize that loss of ANGPTL4 function might improve glucose homeostasis and decrease risk of type 2 diabetes (T2D). We investigate protein-altering variants in ANGPTL4 among 58,124 participants in the DiscovEHR human genetics study, with follow-up studies in 82,766 T2D cases and 498,761 controls. Carriers of p.E40K, a variant that abolishes ANGPTL4 ability to inhibit lipoprotein lipase, have lower odds of T2D (odds ratio 0.89, 95% confidence interval 0.85-0.92, p = 6.3 × 10 -10 ), lower fasting glucose, and greater insulin sensitivity. Predicted loss-of-function variants are associated with lower odds of T2D among 32,015 cases and 84,006 controls (odds ratio 0.71, 95% confidence interval 0.49-0.99, p = 0.041). Functional studies in Angptl4-deficient mice confirm improved insulin sensitivity and glucose homeostasis. In conclusion, genetic inactivation of ANGPTL4 is associated with improved glucose homeostasis and reduced risk of T2D.

Large-scale human genetics studies are ascertaining increasing proportions of populations as they continue growing in both number and scale. As a result, the amount of cryptic relatedness within these study cohorts is growing rapidly and has significant implications on downstream analyses. We demonstrate this growth empirically among the first 92,455 exomes from the DiscovEHR cohort and, via a custom simulation framework we developed called SimProgeny, show that these measures are in line with expectations given the underlying population and ascertainment approach. For example, within DiscovEHR we identified ∼66,000 close (first- and second-degree) relationships, involving 55.6% of study participants. Our simulation results project that >70% of the cohort will be involved in these close relationships, given that DiscovEHR scales to 250,000 recruited individuals. We reconstructed 12,574 pedigrees by using these relationships (including 2,192 nuclear families) and leveraged them for multiple applications. The pedigrees substantially improved the phasing accuracy of 20,947 rare, deleterious compound heterozygous mutations. Reconstructed nuclear families were critical for identifying 3,415 de novo mutations in ∼1,783 genes. Finally, we demonstrate the segregation of known and suspected disease-causing mutations, including a tandem duplication that occurs in LDLR and causes familial hypercholesterolemia, through reconstructed pedigrees. In summary, this work highlights the prevalence of cryptic relatedness expected among large healthcare population-genomic studies and demonstrates several analyses that are uniquely enabled by large amounts of cryptic relatedness., (Copyright © 2018 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

PurposeThe clinical utility of screening unselected individuals for pathogenic BRCA1/2 variants has not been established. Data on cancer risk management behaviors and diagnoses of BRCA1/2-associated cancers can help inform assessments of clinical utility.MethodsWhole-exome sequences of participants in the MyCode Community Health Initiative were reviewed for pathogenic/likely pathogenic BRCA1/2 variants. Clinically confirmed variants were disclosed to patient-participants and their clinicians. We queried patient-participants' electronic health records for BRCA1/2-associated cancer diagnoses and risk management that occurred within 12 months after results disclosure, and calculated the percentage of patient-participants of eligible age who had begun risk management.ResultsThirty-seven MyCode patient-participants were unaware of their pathogenic/likely pathogenic BRCA1/2 variant, had not had a BRCA1/2-associated cancer, and had 12 months of follow-up. Of the 33 who were of an age to begin BRCA1/2-associated risk management, 26 (79%) had performed at least one such procedure. Three were diagnosed with an early-stage, BRCA1/2-associated cancer-including a stage 1C fallopian tube cancer-via these procedures.ConclusionScreening for pathogenic BRCA1/2 variants among unselected individuals can lead to occult cancer detection shortly after disclosure. Comprehensive outcomes data generated within our learning healthcare system will aid in determining whether population-wide BRCA1/2 genomic screening programs offer clinical utility.

Background: Elucidation of the genetic factors underlying chronic liver disease may reveal new therapeutic targets., Methods: We used exome sequence data and electronic health records from 46,544 participants in the DiscovEHR human genetics study to identify genetic variants associated with serum levels of alanine aminotransferase (ALT) and aspartate aminotransferase (AST). Variants that were replicated in three additional cohorts (12,527 persons) were evaluated for association with clinical diagnoses of chronic liver disease in DiscovEHR study participants and two independent cohorts (total of 37,173 persons) and with histopathological severity of liver disease in 2391 human liver samples., Results: A splice variant (rs72613567:TA) in HSD17B13, encoding the hepatic lipid droplet protein hydroxysteroid 17-beta dehydrogenase 13, was associated with reduced levels of ALT (P=4.2×10 -12 ) and AST (P=6.2×10 -10 ). Among DiscovEHR study participants, this variant was associated with a reduced risk of alcoholic liver disease (by 42% [95% confidence interval {CI}, 20 to 58] among heterozygotes and by 53% [95% CI, 3 to 77] among homozygotes), nonalcoholic liver disease (by 17% [95% CI, 8 to 25] among heterozygotes and by 30% [95% CI, 13 to 43] among homozygotes), alcoholic cirrhosis (by 42% [95% CI, 14 to 61] among heterozygotes and by 73% [95% CI, 15 to 91] among homozygotes), and nonalcoholic cirrhosis (by 26% [95% CI, 7 to 40] among heterozygotes and by 49% [95% CI, 15 to 69] among homozygotes). Associations were confirmed in two independent cohorts. The rs72613567:TA variant was associated with a reduced risk of nonalcoholic steatohepatitis, but not steatosis, in human liver samples. The rs72613567:TA variant mitigated liver injury associated with the risk-increasing PNPLA3 p.I148M allele and resulted in an unstable and truncated protein with reduced enzymatic activity., Conclusions: A loss-of-function variant in HSD17B13 was associated with a reduced risk of chronic liver disease and of progression from steatosis to steatohepatitis. (Funded by Regeneron Pharmaceuticals and others.).

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.

Atrial fibrillation (AF) is a common cardiac arrhythmia and a major risk factor for stroke, heart failure, and premature death. The pathogenesis of AF remains poorly understood, which contributes to the current lack of highly effective treatments. To understand the genetic variation and biology underlying AF, we undertook a genome-wide association study (GWAS) of 6,337 AF individuals and 61,607 AF-free individuals from Norway, including replication in an additional 30,679 AF individuals and 278,895 AF-free individuals. Through genotyping and dense imputation mapping from whole-genome sequencing, we tested almost nine million genetic variants across the genome and identified seven risk loci, including two novel loci. One novel locus (lead single-nucleotide variant [SNV] rs12614435; p = 6.76 × 10 -18 ) comprised intronic and several highly correlated missense variants situated in the I-, A-, and M-bands of titin, which is the largest protein in humans and responsible for the passive elasticity of heart and skeletal muscle. The other novel locus (lead SNV rs56202902; p = 1.54 × 10 -11 ) covered a large, gene-dense chromosome 1 region that has previously been linked to cardiac conduction. Pathway and functional enrichment analyses suggested that many AF-associated genetic variants act through a mechanism of impaired muscle cell differentiation and tissue formation during fetal heart development., (Copyright © 2017 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

As the impact of genetics, genomics, and bioinformatics on drug discovery has been increasingly recognized, this session of the 2018 Pacific Symposium on Biocomputing (PSB) aims to facilitate scientific discussions between academia and pharmaceutical industry on how to best apply genetics, genomics and bioinformatics to enable drug discovery. The selected papers focus on developing and applying computational approaches to understand drug mechanisms of action and develop drug combination strategies, to enable in silico drug screening, and to further delineate disease pathways for target identification and validation.

PurposeWe integrated whole-exome sequencing (WES) and chromosomal microarray analysis (CMA) into a clinical workflow to serve an endogamous, uninsured, agrarian community.MethodsSeventy-nine probands (newborn to 49.8 years) who presented between 1998 and 2015 remained undiagnosed after biochemical and molecular investigations. We generated WES data for probands and family members and vetted variants through rephenotyping, segregation analyses, and population studies.ResultsThe most common presentation was neurological disease (64%). Seven (9%) probands were diagnosed by CMA. Family WES data were informative for 37 (51%) of the 72 remaining individuals, yielding a specific genetic diagnosis (n = 32) or revealing a novel molecular etiology (n = 5). For five (7%) additional subjects, negative WES decreased the likelihood of genetic disease. Compared to trio analysis, "family" WES (average seven exomes per proband) reduced filtered candidate variants from 22 ± 6 to 5 ± 3 per proband. Nineteen (51%) alleles were de novo and 17 (46%) inherited; the latter added to a population-based diagnostic panel. We found actionable secondary variants in 21 (4.2%) of 502 subjects, all of whom opted to be informed.ConclusionCMA and family-based WES streamline and economize diagnosis of rare genetic disorders, accelerate novel gene discovery, and create new opportunities for community-based screening and prevention in underserved populations.

PurposeArrhythmogenic 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.MethodsIndividuals (n = 30,716) 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.ResultsEighteen subjects had pLOF variants; none of these had an EHR diagnosis of ARVC. Of 14 patients with an electrocardiogram, one had a minor diagnostic criterion; the rest were normal. A total of 184 subjects had VUS, none of whom had an ARVC diagnosis. The proportion of subjects with VUS with major (4%) or minor (13%) electrocardiogram diagnostic criteria did not differ from that of variant-negative controls. ICD-9 codes showed no difference in defibrillator use, electrophysiologic abnormalities or nonischemic cardiomyopathies in patients with pLOF or VUSs compared with controls.ConclusionpLOF variants in an unselected cohort were not associated with ARVC phenotypes based on EHR review. The negative predictive value of EHR review remains uncertain.

Background: Loss-of-function variants in the angiopoietin-like 3 gene (ANGPTL3) have been associated with decreased plasma levels of triglycerides, low-density lipoprotein (LDL) cholesterol, and high-density lipoprotein (HDL) cholesterol. It is not known whether such variants or therapeutic antagonism of ANGPTL3 are associated with a reduced risk of atherosclerotic cardiovascular disease., Methods: We sequenced the exons of ANGPTL3 in 58,335 participants in the DiscovEHR human genetics study. We performed tests of association for loss-of-function variants in ANGPTL3 with lipid levels and with coronary artery disease in 13,102 case patients and 40,430 controls from the DiscovEHR study, with follow-up studies involving 23,317 case patients and 107,166 controls from four population studies. We also tested the effects of a human monoclonal antibody, evinacumab, against Angptl3 in dyslipidemic mice and against ANGPTL3 in healthy human volunteers with elevated levels of triglycerides or LDL cholesterol., Results: In the DiscovEHR study, participants with heterozygous loss-of-function variants in ANGPTL3 had significantly lower serum levels of triglycerides, HDL cholesterol, and LDL cholesterol than participants without these variants. Loss-of-function variants were found in 0.33% of case patients with coronary artery disease and in 0.45% of controls (adjusted odds ratio, 0.59; 95% confidence interval, 0.41 to 0.85; P=0.004). These results were confirmed in the follow-up studies. In dyslipidemic mice, inhibition of Angptl3 with evinacumab resulted in a greater decrease in atherosclerotic lesion area and necrotic content than a control antibody. In humans, evinacumab caused a dose-dependent placebo-adjusted reduction in fasting triglyceride levels of up to 76% and LDL cholesterol levels of up to 23%., Conclusions: Genetic and therapeutic antagonism of ANGPTL3 in humans and of Angptl3 in mice was associated with decreased levels of all three major lipid fractions and decreased odds of atherosclerotic cardiovascular disease. (Funded by Regeneron Pharmaceuticals and others; ClinicalTrials.gov number, NCT01749878 .).

Rationale: Therapies that inhibit CETP (cholesteryl ester transfer protein) have failed to demonstrate a reduction in risk for coronary heart disease (CHD). Human DNA sequence variants that truncate the CETP gene may provide insight into the efficacy of CETP inhibition., Objective: To test whether protein-truncating variants (PTVs) at the CETP gene were associated with plasma lipid levels and CHD., Methods and Results: We sequenced the exons of the CETP gene in 58 469 participants from 12 case-control studies (18 817 CHD cases, 39 652 CHD-free controls). We defined PTV as those that lead to a premature stop, disrupt canonical splice sites, or lead to insertions/deletions that shift frame. We also genotyped 1 Japanese-specific PTV in 27561 participants from 3 case-control studies (14 286 CHD cases, 13 275 CHD-free controls). We tested association of CETP PTV carrier status with both plasma lipids and CHD. Among 58 469 participants with CETP gene-sequencing data available, average age was 51.5 years and 43% were women; 1 in 975 participants carried a PTV at the CETP gene. Compared with noncarriers, carriers of PTV at CETP had higher high-density lipoprotein cholesterol (effect size, 22.6 mg/dL; 95% confidence interval, 18-27; P <1.0×10 -4 ), lower low-density lipoprotein cholesterol (-12.2 mg/dL; 95% confidence interval, -23 to -0.98; P =0.033), and lower triglycerides (-6.3%; 95% confidence interval, -12 to -0.22; P =0.043). CETP PTV carrier status was associated with reduced risk for CHD (summary odds ratio, 0.70; 95% confidence interval, 0.54-0.90; P =5.1×10 -3 )., Conclusions: Compared with noncarriers, carriers of PTV at CETP displayed higher high-density lipoprotein cholesterol, lower low-density lipoprotein cholesterol, lower triglycerides, and lower risk for CHD., (© 2017 American Heart Association, Inc.)

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.

Chronic obstructive pulmonary disease (COPD) is characterized by reduced lung function and is the third leading cause of death globally. Through genome-wide association discovery in 48,943 individuals, selected from extremes of the lung function distribution in UK Biobank, and follow-up in 95,375 individuals, we increased the yield of independent signals for lung function from 54 to 97. A genetic risk score was associated with COPD susceptibility (odds ratio per 1 s.d. of the risk score (∼6 alleles) (95% confidence interval) = 1.24 (1.20-1.27), P = 5.05 × 10 -49 ), and we observed a 3.7-fold difference in COPD risk between individuals in the highest and lowest genetic risk score deciles in UK Biobank. The 97 signals show enrichment in genes for development, elastic fibers and epigenetic regulation pathways. We highlight targets for drugs and compounds in development for COPD and asthma (genes in the inositol phosphate metabolism pathway and CHRM3) and describe targets for potential drug repositioning from other clinical indications.

A wide range of patient health data is recorded in Electronic Health Records (EHR). This data includes diagnosis, surgical procedures, clinical laboratory measurements, and medication information. Together this information reflects the patient's medical history. Many studies have efficiently used this data from the EHR to find associations that are clinically relevant, either by utilizing International Classification of Diseases, version 9 (ICD-9) codes or laboratory measurements, or by designing phenotype algorithms to extract case and control status with accuracy from the EHR. Here we developed a strategy to utilize longitudinal quantitative trait data from the EHR at Geisinger Health System focusing on outpatient metabolic and complete blood panel data as a starting point. Comprehensive Metabolic Panel (CMP) as well as Complete Blood Counts (CBC) are parts of routine care and provide a comprehensive picture from high level screening of patients' overall health and disease. We randomly split our data into two datasets to allow for discovery and replication. We first conducted a genome-wide association study (GWAS) with median values of 25 different clinical laboratory measurements to identify variants from Human Omni Express Exome beadchip data that are associated with these measurements. We identified 687 variants that associated and replicated with the tested clinical measurements at p<5×10-08. Since longitudinal data from the EHR provides a record of a patient's medical history, we utilized this information to further investigate the ICD-9 codes that might be associated with differences in variability of the measurements in the longitudinal dataset. We identified low and high variance patients by looking at changes within their individual longitudinal EHR laboratory results for each of the 25 clinical lab values (thus creating 50 groups - a high variance and a low variance for each lab variable). We then performed a PheWAS analysis with ICD-9 diagnosis codes, separately in the high variance group and the low variance group for each lab variable. We found 717 PheWAS associations that replicated at a p-value less than 0.001. Next, we evaluated the results of this study by comparing the association results between the high and low variance groups. For example, we found 39 SNPs (in multiple genes) associated with ICD-9 250.01 (Type-I diabetes) in patients with high variance of plasma glucose levels, but not in patients with low variance in plasma glucose levels. Another example is the association of 4 SNPs in UMOD with chronic kidney disease in patients with high variance for aspartate aminotransferase (discovery p-value: 8.71×10-09 and replication p-value: 2.03×10-06). In general, we see a pattern of many more statistically significant associations from patients with high variance in the quantitative lab variables, in comparison with the low variance group across all of the 25 laboratory measurements. This study is one of the first of its kind to utilize quantitative trait variance from longitudinal laboratory data to find associations among genetic variants and clinical phenotypes obtained from an EHR, integrating laboratory values and diagnosis codes to understand the genetic complexities of common diseases.

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., (Copyright © 2016, American Association for the Advancement of Science.)

Familial hypercholesterolemia (FH) remains underdiagnosed despite widespread cholesterol screening. Exome sequencing and electronic health record (EHR) data of 50,726 individuals were used to assess the prevalence and clinical impact of FH-associated genomic variants in the Geisinger Health System. The estimated FH prevalence was 1:256 in unselected participants and 1:118 in participants ascertained via the cardiac catheterization laboratory. FH variant carriers had significantly increased risk of coronary artery disease. Only 24% of carriers met EHR-based presequencing criteria for probable or definite FH diagnosis. Active statin use was identified in 58% of carriers; 46% of statin-treated carriers had a low-density lipoprotein cholesterol level below 100 mg/dl. Thus, we find that genomic screening can prompt the diagnosis of FH patients, most of whom are receiving inadequate lipid-lowering therapy., (Copyright © 2016, American Association for the Advancement of Science.)

Somatic mosaicism, the occurrence and propagation of genetic variation in cell lineages after fertilization, is increasingly recognized to play a causal role in a variety of human diseases. We investigated the case of life-threatening arrhythmia in a 10-day-old infant with long QT syndrome (LQTS). Rapid genome sequencing suggested a variant in the sodium channel Na V 1.5 encoded by SCN5A, NM&#95;000335:c.5284G > T predicting p.(V1762L), but read depth was insufficient to be diagnostic. Exome sequencing of the trio confirmed read ratios inconsistent with Mendelian inheritance only in the proband. Genotyping of single circulating leukocytes demonstrated the mutation in the genomes of 8% of patient cells, and RNA sequencing of cardiac tissue from the infant confirmed the expression of the mutant allele at mosaic ratios. Heterologous expression of the mutant channel revealed significantly delayed sodium current with a dominant negative effect. To investigate the mechanism by which mosaicism might cause arrhythmia, we built a finite element simulation model incorporating Purkinje fiber activation. This model confirmed the pathogenic consequences of cardiac cellular mosaicism and, under the presenting conditions of this case, recapitulated 2:1 AV block and arrhythmia. To investigate the extent to which mosaicism might explain undiagnosed arrhythmia, we studied 7,500 affected probands undergoing commercial gene-panel testing. Four individuals with pathogenic variants arising from early somatic mutation events were found. Here we establish cardiac mosaicism as a causal mechanism for LQTS and present methods by which the general phenomenon, likely to be relevant for all genetic diseases, can be detected through single-cell analysis and next-generation sequencing., Competing Interests: At the time of this work K.M.K., S.R., and L.B. were employed by Gilead Sciences. M.J.C., S.T.K.G., S.B., J.W., and R.C. were employed by Personalis, Inc. J.W. and E.A.A. are founders of Personalis, Inc. which offers clinical genetic testing but does not offer Clinical Laboratory Improvement Amendments (CLIA)-certified rapid-turnaround whole-genome sequencing. N.A.T. is a co-founder of CardioSolv LLC.

Myosin motors are the fundamental force-generating elements of muscle contraction. Variation in the human β-cardiac myosin heavy chain gene (MYH7) can lead to hypertrophic cardiomyopathy (HCM), a heritable disease characterized by cardiac hypertrophy, heart failure, and sudden cardiac death. How specific myosin variants alter motor function or clinical expression of disease remains incompletely understood. Here, we combine structural models of myosin from multiple stages of its chemomechanical cycle, exome sequencing data from two population cohorts of 60,706 and 42,930 individuals, and genetic and phenotypic data from 2,913 patients with HCM to identify regions of disease enrichment within β-cardiac myosin. We first developed computational models of the human β-cardiac myosin protein before and after the myosin power stroke. Then, using a spatial scan statistic modified to analyze genetic variation in protein 3D space, we found significant enrichment of disease-associated variants in the converter, a kinetic domain that transduces force from the catalytic domain to the lever arm to accomplish the power stroke. Focusing our analysis on surface-exposed residues, we identified a larger region significantly enriched for disease-associated variants that contains both the converter domain and residues on a single flat surface on the myosin head described as the myosin mesa. Notably, patients with HCM with variants in the enriched regions have earlier disease onset than patients who have HCM with variants elsewhere. Our study provides a model for integrating protein structure, large-scale genetic sequencing, and detailed phenotypic data to reveal insight into time-shifted protein structures and genetic disease.

Background: Higher-than-normal levels of circulating triglycerides are a risk factor for ischemic cardiovascular disease. Activation of lipoprotein lipase, an enzyme that is inhibited by angiopoietin-like 4 (ANGPTL4), has been shown to reduce levels of circulating triglycerides., Methods: We sequenced the exons of ANGPTL4 in samples obtain from 42,930 participants of predominantly European ancestry in the DiscovEHR human genetics study. We performed tests of association between lipid levels and the missense E40K variant (which has been associated with reduced plasma triglyceride levels) and other inactivating mutations. We then tested for associations between coronary artery disease and the E40K variant and other inactivating mutations in 10,552 participants with coronary artery disease and 29,223 controls. We also tested the effect of a human monoclonal antibody against ANGPTL4 on lipid levels in mice and monkeys., Results: We identified 1661 heterozygotes and 17 homozygotes for the E40K variant and 75 participants who had 13 other monoallelic inactivating mutations in ANGPTL4. The levels of triglycerides were 13% lower and the levels of high-density lipoprotein (HDL) cholesterol were 7% higher among carriers of the E40K variant than among noncarriers. Carriers of the E40K variant were also significantly less likely than noncarriers to have coronary artery disease (odds ratio, 0.81; 95% confidence interval, 0.70 to 0.92; P=0.002). K40 homozygotes had markedly lower levels of triglycerides and higher levels of HDL cholesterol than did heterozygotes. Carriers of other inactivating mutations also had lower triglyceride levels and higher HDL cholesterol levels and were less likely to have coronary artery disease than were noncarriers. Monoclonal antibody inhibition of Angptl4 in mice and monkeys reduced triglyceride levels., Conclusions: Carriers of E40K and other inactivating mutations in ANGPTL4 had lower levels of triglycerides and a lower risk of coronary artery disease than did noncarriers. The inhibition of Angptl4 in mice and monkeys also resulted in corresponding reductions in these values. (Funded by Regeneron Pharmaceuticals.).

Electronic health records (EHR) provide a comprehensive resource for discovery, allowing unprecedented exploration of the impact of genetic architecture on health and disease. The data of EHRs also allow for exploration of the complex interactions between health measures across health and disease. The discoveries arising from EHR based research provide important information for the identification of genetic variation for clinical decision-making. Due to the breadth of information collected within the EHR, a challenge for discovery using EHR based data is the development of high-throughput tools that expose important areas of further research, from genetic variants to phenotypes. Phenome-Wide Association studies (PheWAS) provide a way to explore the association between genetic variants and comprehensive phenotypic measurements, generating new hypotheses and also exposing the complex relationships between genetic architecture and outcomes, including pleiotropy. EHR based PheWAS have mainly evaluated associations with case/control status from International Classification of Disease, Ninth Edition (ICD-9) codes. While these studies have highlighted discovery through PheWAS, the rich resource of clinical lab measures collected within the EHR can be better utilized for high-throughput PheWAS analyses and discovery. To better use these resources and enrich PheWAS association results we have developed a sound methodology for extracting a wide range of clinical lab measures from EHR data. We have extracted a first set of 21 clinical lab measures from the de-identified EHR of participants of the Geisinger MyCodeTM biorepository, and calculated the median of these lab measures for 12,039 subjects. Next we evaluated the association between these 21 clinical lab median values and 635,525 genetic variants, performing a genome-wide association study (GWAS) for each of 21 clinical lab measures. We then calculated the association between SNPs from these GWAS passing our Bonferroni defined p-value cutoff and 165 ICD-9 codes. Through the GWAS we found a series of results replicating known associations, and also some potentially novel associations with less studied clinical lab measures. We found the majority of the PheWAS ICD-9 diagnoses highly related to the clinical lab measures associated with same SNPs. Moving forward, we will be evaluating further phenotypes and expanding the methodology for successful extraction of clinical lab measurements for research and PheWAS use. These developments are important for expanding the PheWAS approach for improved EHR based discovery.

Rationale: Idiopathic pulmonary arterial hypertension (IPAH) is a life-threatening disorder characterized by progressive loss of pulmonary microvessels. Although mutations in the bone morphogenetic receptor 2 (BMPR2) are found in 80% of heritable and ∼15% of patients with IPAH, their low penetrance (∼20%) suggests that other unidentified genetic modifiers are required for manifestation of the disease phenotype. Use of whole-exome sequencing (WES) has recently led to the discovery of novel susceptibility genes in heritable PAH, but whether WES can also accelerate gene discovery in IPAH remains unknown., Objectives: To determine whether WES can help identify novel gene modifiers in patients with IPAH., Methods: Exome capture and sequencing was performed on genomic DNA isolated from 12 unrelated patients with IPAH lacking BMPR2 mutations. Observed genetic variants were prioritized according to their pathogenic potential using ANNOVAR., Measurements and Main Results: A total of nine genes were identified as high-priority candidates. Our top hit was topoisomerase DNA binding II binding protein 1 (TopBP1), a gene involved in the response to DNA damage and replication stress. We found that TopBP1 expression was reduced in vascular lesions and pulmonary endothelial cells isolated from patients with IPAH. Although TopBP1 deficiency made endothelial cells susceptible to DNA damage and apoptosis in response to hydroxyurea, its restoration resulted in less DNA damage and improved cell survival., Conclusions: WES led to the discovery of TopBP1, a gene whose deficiency may increase susceptibility to small vessel loss in IPAH. We predict that use of WES will help identify gene modifiers that influence an individual's risk of developing IPAH.

Importance: Whole-genome sequencing (WGS) is increasingly applied in clinical medicine and is expected to uncover clinically significant findings regardless of sequencing indication., Objectives: To examine coverage and concordance of clinically relevant genetic variation provided by WGS technologies; to quantitate inherited disease risk and pharmacogenomic findings in WGS data and resources required for their discovery and interpretation; and to evaluate clinical action prompted by WGS findings., Design, Setting, and Participants: An exploratory study of 12 adult participants recruited at Stanford University Medical Center who underwent WGS between November 2011 and March 2012. A multidisciplinary team reviewed all potentially reportable genetic findings. Five physicians proposed initial clinical follow-up based on the genetic findings., Main Outcomes and Measures: Genome coverage and sequencing platform concordance in different categories of genetic disease risk, person-hours spent curating candidate disease-risk variants, interpretation agreement between trained curators and disease genetics databases, burden of inherited disease risk and pharmacogenomic findings, and burden and interrater agreement of proposed clinical follow-up., Results: Depending on sequencing platform, 10% to 19% of inherited disease genes were not covered to accepted standards for single nucleotide variant discovery. Genotype concordance was high for previously described single nucleotide genetic variants (99%-100%) but low for small insertion/deletion variants (53%-59%). Curation of 90 to 127 genetic variants in each participant required a median of 54 minutes (range, 5-223 minutes) per genetic variant, resulted in moderate classification agreement between professionals (Gross κ, 0.52; 95% CI, 0.40-0.64), and reclassified 69% of genetic variants cataloged as disease causing in mutation databases to variants of uncertain or lesser significance. Two to 6 personal disease-risk findings were discovered in each participant, including 1 frameshift deletion in the BRCA1 gene implicated in hereditary breast and ovarian cancer. Physician review of sequencing findings prompted consideration of a median of 1 to 3 initial diagnostic tests and referrals per participant, with fair interrater agreement about the suitability of WGS findings for clinical follow-up (Fleiss κ, 0.24; P < 001)., Conclusions and Relevance: In this exploratory study of 12 volunteer adults, the use of WGS was associated with incomplete coverage of inherited disease genes, low reproducibility of detection of genetic variation with the highest potential clinical effects, and uncertainty about clinically reportable findings. In certain cases, WGS will identify clinically actionable genetic variants warranting early medical intervention. These issues should be considered when determining the role of WGS in clinical medicine.

Genetic variation at the chromosome 9p21 risk locus promotes cardiovascular disease; however, it is unclear how or which proteins encoded at this locus contribute to disease. We have previously demonstrated that loss of one candidate gene at this locus, cyclin-dependent kinase inhibitor 2B (Cdkn2b), in mice promotes vascular SMC apoptosis and aneurysm progression. Here, we investigated the role of Cdnk2b in atherogenesis and found that in a mouse model of atherosclerosis, deletion of Cdnk2b promoted advanced development of atherosclerotic plaques composed of large necrotic cores. Furthermore, human carriers of the 9p21 risk allele had reduced expression of CDKN2B in atherosclerotic plaques, which was associated with impaired expression of calreticulin, a ligand required for activation of engulfment receptors on phagocytic cells. As a result of decreased calreticulin, CDKN2B-deficient apoptotic bodies were resistant to efferocytosis and not efficiently cleared by neighboring macrophages. These uncleared SMCs elicited a series of proatherogenic juxtacrine responses associated with increased foam cell formation and inflammatory cytokine elaboration. The addition of exogenous calreticulin reversed defects associated with loss of Cdkn2b and normalized engulfment of Cdkn2b-deficient cells. Together, these data suggest that loss of CDKN2B promotes atherosclerosis by increasing the size and complexity of the lipid-laden necrotic core through impaired efferocytosis.

Personalized medicine is expected to benefit from combining genomic information with regular monitoring of physiological states by multiple high-throughput methods. Here, we present an integrative personal omics profile (iPOP), an analysis that combines genomic, transcriptomic, proteomic, metabolomic, and autoantibody profiles from a single individual over a 14 month period. Our iPOP analysis revealed various medical risks, including type 2 diabetes. It also uncovered extensive, dynamic changes in diverse molecular components and biological pathways across healthy and diseased conditions. Extremely high-coverage genomic and transcriptomic data, which provide the basis of our iPOP, revealed extensive heteroallelic changes during healthy and diseased states and an unexpected RNA editing mechanism. This study demonstrates that longitudinal iPOP can be used to interpret healthy and diseased states by connecting genomic information with additional dynamic omics activity., (Copyright © 2012 Elsevier Inc. All rights reserved.)

Whole-genome sequencing harbors unprecedented potential for characterization of individual and family genetic variation. Here, we develop a novel synthetic human reference sequence that is ethnically concordant and use it for the analysis of genomes from a nuclear family with history of familial thrombophilia. We demonstrate that the use of the major allele reference sequence results in improved genotype accuracy for disease-associated variant loci. We infer recombination sites to the lowest median resolution demonstrated to date (< 1,000 base pairs). We use family inheritance state analysis to control sequencing error and inform family-wide haplotype phasing, allowing quantification of genome-wide compound heterozygosity. We develop a sequence-based methodology for Human Leukocyte Antigen typing that contributes to disease risk prediction. Finally, we advance methods for analysis of disease and pharmacogenomic risk across the coding and non-coding genome that incorporate phased variant data. We show these methods are capable of identifying multigenic risk for inherited thrombophilia and informing the appropriate pharmacological therapy. These ethnicity-specific, family-based approaches to interpretation of genetic variation are emblematic of the next generation of genetic risk assessment using whole-genome sequencing., Competing Interests: JVT and AWZ are founders, consultants, and equity holders in Clinical Future; GMC has advisory roles in and research sponsorships from several companies involved in genome sequencing technology and personal genomics (see http://arep.med.harvard.edu/gmc/tech.html); MS is on the scientific advisory board of DNA Nexus and holds stock in Personalis; RBA has received consultancy fees from Novartis and 23andMe and holds stock in Personalis; AJB is a scientific advisory board member and founder for NuMedii and Genstruct, a scientific advisory board member for Johnson and Johnson, has received consultancy fees from Lilly, NuMedii, Johnson and Johnson, Genstruct, Tercica, and Prevendia and honoraria from Lilly and Siemens, and holds stock in NuMedii, Genstruct, and Personalis. EAA holds stock in Personalis.

Background: The cost of genomic information has fallen steeply, but the clinical translation of genetic risk estimates remains unclear. We aimed to undertake an integrated analysis of a complete human genome in a clinical context., Methods: We assessed a patient with a family history of vascular disease and early sudden death. Clinical assessment included analysis of this patient's full genome sequence, risk prediction for coronary artery disease, screening for causes of sudden cardiac death, and genetic counselling. Genetic analysis included the development of novel methods for the integration of whole genome and clinical risk. Disease and risk analysis focused on prediction of genetic risk of variants associated with mendelian disease, recognised drug responses, and pathogenicity for novel variants. We queried disease-specific mutation databases and pharmacogenomics databases to identify genes and mutations with known associations with disease and drug response. We estimated post-test probabilities of disease by applying likelihood ratios derived from integration of multiple common variants to age-appropriate and sex-appropriate pre-test probabilities. We also accounted for gene-environment interactions and conditionally dependent risks., Findings: Analysis of 2.6 million single nucleotide polymorphisms and 752 copy number variations showed increased genetic risk for myocardial infarction, type 2 diabetes, and some cancers. We discovered rare variants in three genes that are clinically associated with sudden cardiac death-TMEM43, DSP, and MYBPC3. A variant in LPA was consistent with a family history of coronary artery disease. The patient had a heterozygous null mutation in CYP2C19 suggesting probable clopidogrel resistance, several variants associated with a positive response to lipid-lowering therapy, and variants in CYP4F2 and VKORC1 that suggest he might have a low initial dosing requirement for warfarin. Many variants of uncertain importance were reported., Interpretation: Although challenges remain, our results suggest that whole-genome sequencing can yield useful and clinically relevant information for individual patients., Funding: National Institute of General Medical Sciences; National Heart, Lung And Blood Institute; National Human Genome Research Institute; Howard Hughes Medical Institute; National Library of Medicine, Lucile Packard Foundation for Children's Health; Hewlett Packard Foundation; Breetwor Family Foundation., (Copyright 2010 Elsevier Ltd. All rights reserved.)

Background: Atrial fibrillation (AF) is the most prevalent arrhythmia in the United States and accounts for more than 750,000 strokes per year. Noninvasive predictors of AF may help identify patients at risk of developing AF. Our objective was to identify the electrocardiographic characteristics associated with onset of AF., Methods: This was a retrospective cohort analysis of 42,751 patients with electrocardiograms (ECGs) ordered by physician's discretion and analyzed using a computerized system. The population was followed for detection of AF on subsequent ECGs. Cox proportional hazard regression analysis was performed to test the association between these ECG characteristics and development of AF., Results: For a mean follow-up of 5.3 years, 1,050 (2.4%) patients were found to have AF on subsequent ECG recordings. Several ECG characteristics, such as P-wave dispersion (the difference between the widest and narrowest P waves), premature atrial contractions, and an abnormal P axis, were predictive of AF with hazard ratio of approximately 2 after correcting for age and sex. P-wave index, the SD of P-wave duration across all leads, was one of the strongest predictors of AF with a concordance index of 0.62 and a hazard ratio of 2.7 (95% CI 2.1-3.3) for a P-wave index >35. These were among the several independently predictive markers identified on multivariate analysis., Conclusions: Several ECG markers are independently predictive of future onset of AF. The P index, a measurement of disorganized atrial depolarization, is one of the strongest predictors of AF. The ECG contains valuable prognostic information that can identify patients at risk of AF.

Objective: The aim of this study is to evaluate the predictive accuracy of a cardiopulmonary exercise test (CPX) score., Background: Cardiopulmonary exercise test responses, including peak VO(2), markers of ventilatory inefficiency (eg, the VE/VCO(2) slope and oxygen uptake efficiency slope [OUES]), and hemodynamic responses, such as heart rate recovery (HRR) and chronotropic incompetence (CRI) are strong predictors of outcomes in patients with heart failure (HF). However, there is a need for simplified approaches that integrate the additive prognostic information from CPX., Methods: At 4 institutions, 710 patients with HF (568 male/142 female, mean age 56 +/- 13 years, resting left ventricular ejection fraction 33 +/- 14%) underwent CPX and were followed for cardiac-related mortality and separately for major cardiac events (death, hospitalization for HF, transplantation, left ventricular assist device implantation) for a mean of 29 +/- 25 months. The age-adjusted prognostic power of peak VO(2), VE/VCO(2) slope, OUES (VO(2) = a log(10)VE + b), resting end-tidal carbon dioxide pressure (PetCO(2)), HRR, and CRI were determined using Cox proportional hazards analysis, optimal cutpoints were determined, the variables were weighted, and a multivariate score was derived., Results: There were 175 composite outcomes. The VE/VCO(2) slope (> or =34) was the strongest predictor of risk and was attributed a relative weight of 7, with weighted scores for abnormal HRR (< or =6 beats at 1 minute), OUES (>1.4), PetCO(2) (<33 mm Hg), and peak VO(2) (< or =14 mL kg(-1) min(-1)) having scores of 5, 3, 3, and 2, respectively. Chronotropic incompetence was not a significant predictor and was excluded from the score. A summed score >15 was associated with an annual mortality rate of 27% and a relative risk of 7.6, whereas a score <5 was associated with a mortality rate of 0.4%. The composite score was the most accurate predictor of cardiovascular events among all CPX responses considered (concordance indexes 0.77 for mortality and 0.75 for composite outcome composed of mortality, transplantation, left ventricular assist device implantation, and HF-related hospitalization). The summed score remained significantly associated with increased risk after adjusting for age, gender, body mass index, ejection fraction, and cardiomyopathy type., Conclusion: A multivariable score based on readily available CPX responses provides a simple and integrated method that powerfully predicts outcomes in patients with HF.

The Duke Treadmill Score (DTS) is an established clinical tool for risk stratification. Our aim was to determine if other variables could improve the prognostic power of the DTS and if so, to modify the DTS nomogram. From a total of 1,959 patients referred for exercise testing at the Palo Alto VA Medical Center from 1997 to 2006 (a mean follow-up of 5.4 years), we studied 1,759 male veterans (age 57 +/- 12 years) free of heart failure. Double product (DP) was calculated by multiplying systolic blood pressure and heart rate; variables and their products were subtracted to obtain the differences between at rest and maximal exercise (reserve) and recovery. Of all the hemodynamic measurements, DP reserve was the strongest predictor of cardiovascular death (CVD) (Wald Z-score -3.84, p <0.001) after adjustment for potential confounders. When the components of DTS were entered in the Cox hazard model with DP reserve and age, only DP reserve and age were chosen (p <0.00001). Using the Cox coefficients, a score calculated by [age - DTS - 3 x (DP reserve/1,000)] yielded an area under the curve of 0.84 compared with 0.76 for the DTS. Using this equation, a nomogram was constructed by adding age and DP reserve to the original DTS nomogram improving estimation of annual CVD. In conclusion, we propose an age and DP reserve-adjusted DTS nomogram that improves the prognostic estimates of average annual CVD over the DTS alone.

Background: Although the prognostic power of exercise capacity has been demonstrated, the relative prognostic potential of other hemodynamic responses has not been thoroughly investigated. We aimed to assess the prognostic power of double product (DP) parameters in patients referred for standard exercise testing., Design: A retrospective cohort study., Methods: Analyses were performed on 1959 patients referred for exercise testing at the Palo Alto Veterans Affairs Medical Center from 1997 to 2006. After removal of female and heart failure patients, 1759 male veterans (mean age 57+/-12 years) remained. DP was calculated by multiplying systolic blood pressure (SBP) and heart rate (HR); variables and their products were subtracted to obtain the differences between rest and maximal exercise (reserve), and recovery., Results: Multivariable Cox survival analysis was performed for 157 all-cause and 53 cardiovascular deaths during a mean follow-up of 5.4+/-2.1 years. Although most of the hemodynamic variables were individually significant in Cox survival analysis, when age, DP reserve, exercise capacity (in metabolic equivalents), and HR recovery were entered together, only age and DP reserve were chosen. Of all hemodynamic measurements considered, DP reserve was the strongest predictor of cardiovascular prognosis after adjustment for age (Wald Z-score,-5.12; P<0.0001) and when considering other potential confounders such as age, beta-blocker use, and the Duke treadmill Score (Wald Z-score,-3.84; P<0.0001)., Conclusion: In this study population, DP reserve had greater prognostic power than metabolic equivalents, maximal HR or systolic blood pressure, or HR recovery.

Background: The Duke Treadmill Score (DTS) is an established clinical tool for risk stratification of patients referred for exercise testing, but it does not consider age. We aimed to determine if age could improve the prognostic power of the DTS and if so, to modify the DTS nomogram to include age., Methods: Of 1,959 patients referred for exercise testing from 1997 to 2006, 1,759 male veterans (age range 23-86 years) remained after exclusion of female and patients with heart failure. Cardiovascular mortality was the main outcome considered., Results: Cox survival analysis was performed entering age and the DTS; both were significant (P

Extensive evidence is available that cardiovascular structure and function, along with other biological properties that span the range of organism size and speciation, scale with body size. Although appreciation of such factors is commonplace in pediatrics, cardiovascular measurements in the adult population, with similarly wide variation in body size, are rarely corrected for body size. In this review, we describe the critical role of body size measurements in cardiovascular medicine. Using examples, we illustrate the confounding effects of body size. Current cardiovascular scaling practices are reviewed, as are limitations and alternative relationships between body and cardiovascular dimensions. The experimental evidence, theoretical basis, and clinical application of scaling of various functional parameters are presented. Appropriately scaled parameters aid diagnostic and therapeutic decision making in specific disease states such as hypertrophic cardiomyopathy and congestive heart failure. Large-scale studies in clinical populations are needed to define normative relationships for this purpose. Lack of appropriate consideration of body size in the evaluation of cardiovascular structure and function may adversely affect recognition and treatment of cardiovascular disease states in the adult patient.