The Development of an Infrastructure to Facilitate the Use of Whole Genome Sequencing for Population Health

The clinical use of genomic analysis has expanded rapidly resulting in an increased availability and utility of genomic information in clinical care. We have developed an infrastructure utilizing informatics tools and clinical processes to facilitate the use of whole genome sequencing data for population health management across the healthcare system. Our resulting framework scaled well to multiple clinical domains in both pediatric and adult care, although there were domain specific challenges that arose. Our infrastructure was complementary to existing clinical processes and well-received by care providers and patients. Informatics solutions were critical to the successful deployment and scaling of this program. Implementation of genomics at the scale of population health utilizes complicated technologies and processes that for many health systems are not supported by current information systems or in existing clinical workflows. To scale such a system requires a substantial clinical framework backed by informatics tools to facilitate the flow and management of data. Our work represents an early model that has been successful in scaling to 29 different genes with associated genetic conditions in four clinical domains. Work is ongoing to optimize informatics tools; and to identify best practices for translation to smaller healthcare systems.

Sex-specific association of RAGE and HMGB1 genotype variations with susceptibility to ischemic stroke in Caucasians

BACKGROUND

The role of genotype variants of HMGB1 and RAGE in susceptibility to acute ischemic stroke remains inconclusive.

METHODS

Caucasian acute ischemic stroke patients admitted to three hospitals within a large healthcare system in the U.S. between 2009 and 2017 were reviewed. For each stroke case, three age and sex-matched non-stroke patients were identified as controls. Associations of phased-genotype data for RAGE (rs1035798, rs2070600, rs1800624, rs1800625) and HMGB1 (rs1360485, rs1045411, rs3742305, rs2249825, rs1412125) single-nucleotide-polymorphisms (SNPs) and haplotypes with stroke susceptibility were analyzed. The Benjamini-Hochberg procedure was performed.

RESULTS

Collectively, 4,264 patients, 1,066 acute ischemic stroke and 3,198 controls were identified. Genotype distributions were in Hardy-Weinberg equilibrium. None of the SNPs alternate allele frequencies differed from the NCBI SNP database. No differences were found in the genotype distributions when analyzing each SNP and the two most common haplotypes in a covariate adjusted model. In a sex-specific stratification, males harboring the RAGE SNP rs1800625 AG or GG genotype had an independently increased risk for ischemic strokes compared to controls (adjusted OR = 1.27,95%CI 1.03-1.57, p_a = 0.0276). After the Benjamini-Hochberg procedure, a trend towards this association remained (p_BH = 0.1104).

CONCLUSION

No association of RAGE and HMGB1 genotypes variations with risk for overall ischemic stroke or specific stroke subtypes could be observed. Congruent with the literature, a sex-specific role of RAGE SNPs might associate with stroke susceptibility. The functional role of the HMGB1-RAGE axis in this context warrants further exploration.

User testing of a diagnostic decision support system with machine-assisted chart review to facilitate clinical genomic diagnosis

Objectives

There is a need in clinical genomics for systems that assist in clinical diagnosis, analysis of genomic information and periodic reanalysis of results, and can use information from the electronic health record to do so. Such systems should be built using the concepts of human-centred design, fit within clinical workflows and provide solutions to priority problems.

Methods

We adapted a commercially available diagnostic decision support system (DDSS) to use extracted findings from a patient record and combine them with genomic variant information in the DDSS interface. Three representative patient cases were created in a simulated clinical environment for user testing. A semistructured interview guide was created to illuminate factors relevant to human factors in CDS design and organisational implementation.

Results

Six individuals completed the user testing process. Tester responses were positive and noted good fit with real-world clinical genetics workflow. Technical issues related to interface, interaction and design were minor and fixable. Testers suggested solving issues related to terminology and usability through training and infobuttons. Time savings was estimated at 30%–50% and additional uses such as in-house clinical variant analysis were suggested for increase fit with workflow and to further address priority problems.

Conclusion

This study provides preliminary evidence for usability, workflow fit, acceptability and implementation potential of a modified DDSS that includes machine-assisted chart review. Continued development and testing using principles from human-centred design and implementation science are necessary to improve technical functionality and acceptability for multiple stakeholders and organisational implementation potential to improve the genomic diagnosis process.

Cysteine-Altering NOTCH3 Variants Are a Risk Factor for Stroke in the Elderly Population

Supplemental Digital Content is available in the text.

Cysteine altering NOTCH3 variants, which have previously been exclusively associated with the rare hereditary small vessel disease cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy, have a population frequency of 1:300 worldwide. Using a large population database, and taking genotype as a starting point, we aimed to determine whether individuals harboring a NOTCH3 cysteine altering variant have a higher load of small vessel disease markers on brain magnetic resonance imaging than controls, as well as a higher risk of stroke and cognitive impairment.

Matrix Gla protein polymorphism rs1800801 associates with recurrence of ischemic stroke

The MGP single nucleotide polymorphism (SNP) rs1800801 has previously been associated with recurrent ischemic stroke in a Spanish cohort. Here, we tested for association of this SNP with ischemic stroke recurrence in a North American Caucasian cohort. Acute ischemic stroke patients admitted between 10/2009 and 12/2016 at three hospitals within a large healthcare system in the northeastern United States that were enrolled in a healthcare system-wide exome sequencing program were retrospectively reviewed. Patients with recurrent stroke within 1 year after index event were compared to those without recurrence. Of 9,348 suspected acute ischemic strokes admitted between 10/2009 and 12/2016, 1,727 (18.5%) enrolled in the exome-sequencing program. Among those, 1,068 patients had exome sequencing completed and were eligible for inclusion. Recurrent stroke within the first year of stroke was observed in 79 patients (7.4%). In multivariable analysis, stroke prior to the index stroke (OR 9.694, 95% CI 5.793-16.224, p ≤ 0.001), pro-coagulant status (OR = 3.563, 95% CI 1.504-8.443, p = 0.004) and the AA genotype of SNP rs1800801 (OR = 2.408, 95% CI 1.079-4.389, p = 0.004) were independently associated with recurrent stroke within the first year. The AA genotype of the MGP SNP rs1800801 is associated with recurrence within the first year after ischemic stroke in North American Caucasians. Study of stroke subtypes and additional populations will be required to determine if incorporation of allelic status at this SNP into current risk scores improves prediction of recurrent ischemic stroke.

Risk Factors for Acute Ischemic Stroke Caused by Anterior Large Vessel Occlusion

Background and Purpose- Accurate prediction of acute ischemic stroke (AIS) caused by anterior large vessel occlusion (LVO) that is amendable to mechanical thrombectomy remains a challenge. We developed and validated a prediction model for anterior circulation LVO stroke using past medical history elements present on admission and neurological examination. Methods- We retrospectively reviewed AIS patients admitted between 2009 and 2017 to 3 hospitals within a large healthcare system in the United States. Patients with occlusions of the internal carotid artery or M1 or M2 segments of the middle cerebral artery were randomly split into 2/3 derivation and 1/3 validation cohorts for development of an anterior circulation LVO prediction model and score that was further curtailed for potential use in the prehospital setting. Results- A total of 1654 AIS were reviewed, including 248 (15%) with proximal anterior circulation LVO AIS. In the derivation cohort, National Institutes of Health Stroke Scale score at the time of cerebrovascular imaging, current smoking status, type 2 diabetes mellitus, extracranial carotid, and intracranial atherosclerotic stenosis was significantly associated with anterior circulation LVO stroke. The prehospital score was curtailed to National Institutes of Health Stroke Scale score, current smoking status, and type 2 diabetes mellitus. The areas under the curve for the prediction model, prehospital score, and National Institutes of Health Stroke Scale score alone were 0.796, 0.757, and 0.725 for the derivation cohort and 0.770, 0.689, and 0.665 for the validation cohort, respectively. The Youden index J was 0.46 for a score of >6 with 84.7% sensitivity and 62.0% specificity for the prediction model. Conclusions- Previously reported LVO stroke prediction scores focus solely on elements of the neurological examination. In addition to stroke severity, smoking, diabetes mellitus, extracranial carotid, and intracranial atherosclerotic stenosis were associated with anterior circulation LVO AIS. Although atherosclerotic stenosis may not be known until imaging is obtained, smoking and diabetes mellitus history can be readily obtained in the field and represent important elements of the prehospital score supplementing National Institutes of Health Stroke Scale score.

Harmonizing Clinical Sequencing and Interpretation for the eMERGE III Network

The advancement of precision medicine requires new methods to coordinate and deliver genetic data from heterogeneous sources to physicians and patients. The eMERGE III Network enrolled >25,000 participants from biobank and prospective cohorts of predominantly healthy individuals for clinical genetic testing to determine clinically actionable findings. The network developed protocols linking together the 11 participant collection sites and 2 clinical genetic testing laboratories. DNA capture panels targeting 109 genes were used for testing of DNA and sample collection, data generation, interpretation, reporting, delivery, and storage were each harmonized. A compliant and secure network enabled ongoing review and reconciliation of clinical interpretations, while maintaining communication and data sharing between clinicians and investigators. A total of 202 individuals had positive diagnostic findings relevant to the indication for testing and 1,294 had additional/secondary findings of medical significance deemed to be returnable, establishing data return rates for other testing endeavors. This study accomplished integration of structured genomic results into multiple electronic health record (EHR) systems, setting the stage for clinical decision support to enable genomic medicine. Further, the established processes enable different sequencing sites to harmonize technical and interpretive aspects of sequencing tests, a critical achievement toward global standardization of genomic testing. The eMERGE protocols and tools are available for widespread dissemination.

Publisher Correction: Protein-altering variants associated with body mass index implicate pathways that control energy intake and expenditure in obesity

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

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.

Systematic characterization of germline variants from the DiscovEHR study endometrial carcinoma population

Background

Endometrial cancer (EMCA) is the fifth most common cancer among women in the world. Identification of potentially pathogenic germline variants from individuals with EMCA will help characterize genetic features that underlie the disease and potentially predispose individuals to its pathogenesis.

Methods

The Geisinger Health System’s (GHS) DiscovEHR cohort includes exome sequencing on over 50,000 consenting patients, 297 of whom have evidence of an EMCA diagnosis in their electronic health record. Here, rare variants were annotated as potentially pathogenic.

Results

Eight genes were identified as having increased burden in the EMCA cohort relative to the non-cancer control cohort. None of the eight genes had an increased burden in the other hormone related cancer cohort from GHS, suggesting they can help characterize the underlying genetic variation that gives rise to EMCA. Comparing GHS to the cancer genome atlas (TCGA) EMCA germline data illustrated 34 genes with potentially pathogenic variation and eight unique potentially pathogenic variants that were present in both studies. Thus, similar germline variation among genes can be observed in unique EMCA cohorts and could help prioritize genes to investigate for future work.

Conclusion

In summary, this systematic characterization of potentially pathogenic germline variants describes the genetic underpinnings of EMCA through the use of data from a single hospital system.

Electronic supplementary material

The online version of this article (10.1186/s12920-019-0504-9) contains supplementary material, which is available to authorized users.

Publisher Correction: Protein-altering variants associated with body mass index implicate pathways that control energy intake and expenditure in obesity

In the version of this article originally published, one of the two authors with the name Wei Zhao was omitted from the author list and the affiliations for both authors were assigned to the single Wei Zhao in the author list. In addition, the ORCID for Wei Zhao (Department of Biostatistics and Epidemiology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA) was incorrectly assigned to author Wei Zhou. The errors have been corrected in the HTML and PDF versions of the article.

Publisher Correction: Protein-altering variants associated with body mass index implicate pathways that control energy intake and expenditure in obesity

In the published version of this paper, the name of author Emanuele Di Angelantonio was misspelled. This error has now been corrected in the HTML and PDF versions of the article.

Protein-altering variants associated with body mass index implicate pathways that control energy intake and expenditure in obesity

Genome-wide association studies (GWAS) have identified >250 loci for body mass index (BMI), implicating pathways related to neuronal biology. Most GWAS loci represent clusters of common, noncoding variants from which pinpointing causal genes remains challenging. Here we combined data from 718,734 individuals to discover rare and low-frequency (minor allele frequency (MAF) < 5%) coding variants associated with BMI. We identified 14 coding variants in 13 genes, of which 8 variants were in genes (ZBTB7B, ACHE, RAPGEF3, RAB21, ZFHX3, ENTPD6, ZFR2 and ZNF169) newly implicated in human obesity, 2 variants were in genes (MC4R and KSR2) previously observed to be mutated in extreme obesity and 2 variants were in GIPR. The effect sizes of rare variants are ~10 times larger than those of common variants, with the largest effect observed in carriers of an MC4R mutation introducing a stop codon (p.Tyr35Ter, MAF = 0.01%), who weighed ~7 kg more than non-carriers. Pathway analyses based on the variants associated with BMI confirm enrichment of neuronal genes and provide new evidence for adipocyte and energy expenditure biology, widening the potential of genetically supported therapeutic targets in obesity.

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.

Rare and low-frequency coding variants alter human adult height

Height is a highly heritable, classic polygenic trait with approximately 700 common associated variants identified through genome-wide association studies so far. Here, we report 83 height-associated coding variants with lower minor-allele frequencies (in the range of 0.1-4.8%) and effects of up to 2 centimetres per allele (such as those in IHH, STC2, AR and CRISPLD2), greater than ten times the average effect of common variants. In functional follow-up studies, rare height-increasing alleles of STC2 (giving an increase of 1-2 centimetres per allele) compromised proteolytic inhibition of PAPP-A and increased cleavage of IGFBP-4 in vitro, resulting in higher bioavailability of insulin-like growth factors. These 83 height-associated variants overlap genes that are mutated in monogenic growth disorders and highlight new biological candidates (such as ADAMTS3, IL11RA and NOX4) and pathways (such as proteoglycan and glycosaminoglycan synthesis) involved in growth. Our results demonstrate that sufficiently large sample sizes can uncover rare and low-frequency variants of moderate-to-large effect associated with polygenic human phenotypes, and that these variants implicate relevant genes and pathways.

SeqHBase: a big data toolset for family based sequencing data analysis

Background

Whole-genome sequencing (WGS) and whole-exome sequencing (WES) technologies are increasingly used to identify disease-contributing mutations in human genomic studies. It can be a significant challenge to process such data, especially when a large family or cohort is sequenced. Our objective was to develop a big data toolset to efficiently manipulate genome-wide variants, functional annotations and coverage, together with conducting family based sequencing data analysis.

Methods

Hadoop is a framework for reliable, scalable, distributed processing of large data sets using MapReduce programming models. Based on Hadoop and HBase, we developed SeqHBase, a big data-based toolset for analysing family based sequencing data to detect de novo, inherited homozygous, or compound heterozygous mutations that may contribute to disease manifestations. SeqHBase takes as input BAM files (for coverage at every site), variant call format (VCF) files (for variant calls) and functional annotations (for variant prioritisation).

Results

We applied SeqHBase to a 5-member nuclear family and a 10-member 3-generation family with WGS data, as well as a 4-member nuclear family with WES data. Analysis times were almost linearly scalable with number of data nodes. With 20 data nodes, SeqHBase took about 5 secs to analyse WES familial data and approximately 1 min to analyse WGS familial data.

Conclusions

These results demonstrate SeqHBase's high efficiency and scalability, which is necessary as WGS and WES are rapidly becoming standard methods to study the genetics of familial disorders.