Daniel Seung Kim

Actionable exomic incidental findings in 6503 participants: challenges of variant classification

Recommendations for laboratories to report incidental findings from genomic tests have stimulated interest in such results. In order to investigate the criteria and processes for assigning the pathogenicity of specific variants and to estimate the frequency of such incidental findings in patients of European and African ancestry, we classified potentially actionable pathogenic single-nucleotide variants (SNVs) in all 4300 European- and 2203 African-ancestry participants sequenced by the NHLBI Exome Sequencing Project (ESP). We considered 112 gene-disease pairs selected by an expert panel as associated with medically actionable genetic disorders that may be undiagnosed in adults. The resulting classifications were compared to classifications from other clinical and research genetic testing laboratories, as well as with in silico pathogenicity scores. Among European-ancestry participants, 30 of 4300 (0.7%) had a pathogenic SNV and six (0.1%) had a disruptive variant that was expected to be pathogenic, whereas 52 (1.2%) had likely pathogenic SNVs. For African-ancestry participants, six of 2203 (0.3%) had a pathogenic SNV and six (0.3%) had an expected pathogenic disruptive variant, whereas 13 (0.6%) had likely pathogenic SNVs. Genomic Evolutionary Rate Profiling mammalian conservation score and the Combined Annotation Dependent Depletion summary score of conservation, substitution, regulation, and other evidence were compared across pathogenicity assignments and appear to have utility in variant classification. This work provides a refined estimate of the burden of adult onset, medically actionable incidental findings expected from exome sequencing, highlights challenges in variant classification, and demonstrates the need for a better curated variant interpretation knowledge base.

Design and anticipated outcomes of the eMERGE-PGx project: a multicenter pilot for preemptive pharmacogenomics in electronic health record systems.

We describe here the design and initial implementation of the eMERGE‐PGx project. eMERGE‐PGx, a partnership of the Electronic Medical Records and Genomics Network and the Pharmacogenomics Research Network, has three objectives: (i) to deploy PGRNseq, a next‐generation sequencing platform assessing sequence variation in 84 proposed pharmacogenes, in nearly 9,000 patients likely to be prescribed drugs of interest in a 1‐ to 3‐year time frame across several clinical sites; (ii) to integrate well‐established clinically validated pharmacogenetic genotypes into the electronic health record with associated clinical decision support and to assess process and clinical outcomes of implementation; and (iii) to develop a repository of pharmacogenetic variants of unknown significance linked to a repository of electronic health record–based clinical phenotype data for ongoing pharmacogenomics discovery. We describe site‐specific project implementation and anticipated products, including genetic variant and phenotype data repositories, novel variant association studies, clinical decision support modules, clinical and process outcomes, approaches to managing incidental findings, and patient and clinician education methods.

Quantification of HDL Particle Concentration by Calibrated Ion Mobility Analysis

BACKGROUND It is critical to develop new metrics to determine whether HDL is cardioprotective in humans. One promising approach is HDL particle concentration (HDL-P), the size and concentration of HDL in plasma. However, the 2 methods currently used to determine HDL-P yield concentrations that differ >5-fold. We therefore developed and validated an improved approach to quantify HDL-P, termed calibrated ion mobility analysis (calibrated IMA). METHODS HDL was isolated from plasma by ultracentrifugation, introduced into the gas phase with electrospray ionization, separated by size, and quantified by particle counting. We used a calibration curve constructed with purified proteins to correct for the ionization efficiency of HDL particles. RESULTS The concentrations of gold nanoparticles and reconstituted HDLs measured by calibrated IMA were indistinguishable from concentrations determined by orthogonal methods. In plasma of control (n = 40) and cerebrovascular disease (n = 40) participants, 3 subspecies of HDL were reproducibility measured, with an estimated total HDL-P of 13.4 (2.4) μmol/L. HDL-C accounted for 48% of the variance in HDL-P. HDL-P was significantly lower in participants with cerebrovascular disease (P = 0.002), and this difference remained significant after adjustment for HDL cholesterol concentrations (P = 0.02). CONCLUSIONS Calibrated IMA accurately determined the concentration of gold nanoparticles and synthetic HDL, strongly suggesting that the method could accurately quantify HDL particle concentration. The estimated stoichiometry of apolipoprotein A-I determined by calibrated IMA was 3-4 per HDL particle, in agreement with current structural models. Furthermore, HDL-P was associated with cardiovascular disease status in a clinical population independently of HDL cholesterol.

HDL-3 is a Superior Predictor of Carotid Artery Disease in a Case-Control Cohort of 1725 Participants

Background Recent data suggest that high‐density lipoprotein cholesterol (HDL‐C) levels are likely not in the causative pathway of atheroprotection, shifting focus from HDL‐C to its subfractions and associated proteins. This studys goal was to determine which HDL phenotype was the better predictor of carotid artery disease (CAAD). Methods and Results HDL‐2 and HDL‐3 were measured in 1725 participants of European ancestry in a prevalent case‐control cohort study of CAAD. Stratified analyses were conducted for men (n=1201) and women (n=524). Stepwise linear regression was used to determine whether HDL‐C, HDL‐2, HDL‐3, or apolipoprotein A1 was the best predictor of CAAD, while adjusting for the confounders of censored age, diabetes, and current smoking status. In both men and women, HDL‐3 was negatively associated with CAAD (P=0.0011 and 0.033 for men and women, respectively); once HDL‐3 was included in the model, no other HDL phenotype was significantly associated with CAAD. Addition of paraoxonase 1 activity to the aforementioned regression model showed a significant and independent (of HDL‐3) association with CAAD in men (P=0.001) but not in the smaller female subgroup. Conclusions This study is the first to contrast the associations of HDL‐2 and HDL‐3 with CAAD. We found that HDL‐3 levels were more predictive of CAAD status than HDL‐2, HDL‐C, or apolipoprotein A1. In addition, for men, paraoxonase 1 activity improved the overall model prediction for CAAD independently and additively with HDL‐3 levels. Further investigation into the molecular mechanisms through which HDL‐3 is associated with protection from CAAD is warranted.

Pharmacogenetics of paraoxonase activity: elucidating the role of high-density lipoprotein in disease

PON1 is a key component of high-density lipoproteins (HDLs) and is at least partially responsible for HDLs antioxidant/atheroprotective properties. PON1 is also associated with numerous human diseases, including cardiovascular disease, Parkinsons disease and cancer. In addition, PON1 metabolizes a broad variety of substrates, including toxic organophosphorous compounds, statin adducts, glucocorticoids, the likely atherogenic L-homocysteine thiolactone and the quorum-sensing factor of Pseudomonas aeruginosa. Numerous cardiovascular and antidiabetic pharmacologic agents, dietary macronutrients, lifestyle factors and antioxidant supplements affect PON1 expression and enzyme activity levels. Owing to the importance of PON1 to HDL function and its individual association with diverse human diseases, pharmacogenomic interactions between PON1 and the various factors that alter its expression and activity may represent an important therapeutic target for future investigation.

Genetic variation in the HLA region is associated with susceptibility to herpes zoster

Herpes zoster, commonly referred to as shingles, is caused by the varicella zoster virus (VZV). VZV initially manifests as chicken pox, most commonly in childhood, can remain asymptomatically latent in nerve tissues for many years and often re-emerges as shingles. Although reactivation may be related to immune suppression, aging and female sex, most inter-individual variability in re-emergence risk has not been explained to date. We performed a genome-wide association analyses in 22 981 participants (2280 shingles cases) from the electronic Medical Records and Genomics Network. Using Cox survival and logistic regression, we identified a genomic region in the combined and European ancestry groups that has an age of onset effect reaching genome-wide significance (P>1.0 × 10−8). This region tags the non-coding gene HCP5 (HLA Complex P5) in the major histocompatibility complex. This gene is an endogenous retrovirus and likely influences viral activity through regulatory functions. Variants in this genetic region are known to be associated with delay in development of AIDS in people infected by HIV. Our study provides further suggestion that this region may have a critical role in viral suppression and could potentially harbor a clinically actionable variant for the shingles vaccine.

Validation of association of the apolipoprotein E ε2 allele with neurodevelopmental dysfunction after cardiac surgery in neonates and infants.

OBJECTIVE Apolipoprotein E (APOE) genotype is a determinant of neurologic recovery after brain ischemia and traumatic brain injury. The APOE ε2 allele has been associated with worse neurodevelopmental (ND) outcome after repair of congenital heart defects (CHD) in infancy. Replication of this finding in an independent cohort is essential to validate the observed genotype-phenotype association. METHODS The association of APOE genotype with ND outcomes was assessed in a combined cohort of patients with single-ventricle CHD enrolled in the Single Ventricle Reconstruction and Infant Single Ventricle trials. ND outcome was assessed at 14 months using the Psychomotor Development Index (PDI) and Mental Development Index (MDI) of the Bayley Scales of Infant Development-II. Stepwise multivariable regression was performed to develop predictive models for PDI and MDI scores. RESULTS Complete data were available for 298 of 435 patients. After adjustment for preoperative and postoperative covariates, the APOE ε2 allele was associated with a lower PDI score (P = .038). Patients with the ε2 allele had a PDI score approximately 6 points lower than those without the risk allele, explaining 1.04% of overall PDI variance, because the ε2 allele was present in only 11% of the patients. There was a marginal effect of the ε2 allele on MDI scores (P = .058). CONCLUSIONS These data validate the association of the APOE ε2 allele with adverse early ND outcomes after cardiac surgery in infants, independent of patient and operative factors. Genetic variants that decrease neuroresilience and impair neuronal repair after brain injury are important risk factors for ND dysfunction after surgery for CHD.

Dietary cholesterol increases paraoxonase 1 enzyme activity

HDL-associated paraoxonase 1 (PON1) activity has been consistently associated with cardiovascular and other diseases. Vitamins C and E intake have previously been positively associated with PON1 in a subset of the Carotid Lesion Epidemiology and Risk (CLEAR) cohort. The goal of this study was to replicate these findings and determine whether other nutrient intake affected PON1 activity. To predict nutrient and mineral intake values, 1,402 subjects completed a standardized food frequency survey of their dietary habits over the past year. Stepwise regression was used to evaluate dietary and covariate effects on PON1 arylesterase activity. Five dietary components, cholesterol (P < 2.0 × 10−16), alcohol (P = 8.51 × 10−8), vitamin C (P = 7.97 × 10−5), iron (P = 0.0026), and folic acid (0.037) were independently predictive of PON1 activity. Dietary cholesterol was positively associated and predicted 5.5% of PON1 activity, second in variance explained. This study presents a novel finding of dietary cholesterol, iron, and folic acid predicting PON1 activity in humans and confirms prior reported associations, including that with vitamin C. Identifying and understanding environmental factors that affect PON1 activity is necessary to understand its role and that of HDL in human disease.

Rare coding variation in paraoxonase-1 is associated with ischemic stroke in the NHLBI Exome Sequencing Project

HDL-associated paraoxonase-1 (PON1) is an enzyme whose activity is associated with cerebrovascular disease. Common PON1 genetic variants have not been consistently associated with cerebrovascular disease. Rare coding variation that likely alters PON1 enzyme function may be more strongly associated with stroke. The National Heart, Lung, and Blood Institute Exome Sequencing Project sequenced the coding regions (exomes) of the genome for heart, lung, and blood-related phenotypes (including ischemic stroke). In this sample of 4,204 unrelated participants, 496 had verified, noncardioembolic ischemic stroke. After filtering, 28 nonsynonymous PON1 variants were identified. Analysis with the sequence kernel association test, adjusted for covariates, identified significant associations between PON1 variants and ischemic stroke (P = 3.01 × 10−3). Stratified analyses demonstrated a stronger association of PON1 variants with ischemic stroke in African ancestry (AA) participants (P = 5.03 × 10−3). Ethnic differences in the association between PON1 variants with stroke could be due to the effects of PON1Val109Ile (overall P = 7.88 × 10−3; AA P = 6.52 × 10−4), found at higher frequency in AA participants (1.16% vs. 0.02%) and whose protein is less stable than the common allele. In summary, rare genetic variation in PON1 was associated with ischemic stroke, with stronger associations identified in those of AA. Increased focus on PON1 enzyme function and its role in cerebrovascular disease is warranted.

Burden of potentially pathologic copy number variants is higher in children with isolated congenital heart disease and significantly impairs covariate-adjusted transplant-free survival

OBJECTIVES Copy number variants (CNVs) are duplications or deletions of genomic regions. Large CNVs are potentially pathogenic and are overrepresented in children with congenital heart disease (CHD). We sought to determine the frequency of large CNVs in children with isolated CHD, and to evaluate the relationship of these potentially pathogenic CNVs with transplant-free survival. METHODS These cases are derived from a prospective cohort of patients with nonsyndromic CHD (n = 422) identified before first surgery. Healthy pediatric controls (n = 500) were obtained from the electronic Medical Records and Genetic Epidemiology Network, and CNV frequency was contrasted for CHD cases and controls. CNVs were determined algorithmically; subsequently screened for >95% overlap between 2 methods, size (>300 kb), quality score, overlap with a gene, and novelty (absent from databases of known, benign CNVs); and separately validated by quantitative polymerase chain reaction. Survival likelihoods for cases were calculated using Cox proportional hazards modeling to evaluate the joint effect of CNV burden and known confounders on transplant-free survival. RESULTS Children with nonsyndromic CHD had a higher burden of potentially pathogenic CNVs compared with pediatric controls (12.1% vs 5.0%; P = .00016). Presence of a CNV was associated with significantly decreased transplant-free survival after surgery (hazard ratio, 3.42; 95% confidence interval, 1.66-7.09; P = .00090) with confounder adjustment. CONCLUSIONS We confirm that children with isolated CHD have a greater burden of rare/large CNVs. We report a novel finding that these CNVs are associated with an adjusted 2.55-fold increased risk of death or transplant. These data suggest that CNV burden is an important modifier of survival after surgery for CHD.