Shang, Ning

Abstract Importance Chronic kidney disease (CKD) is a genetically complex disease determined by an interplay of monogenic, polygenic, and environmental risks. The most common forms of monogenic kidney disorders include autosomal dominant polycystic kidney disease (ADPKD), caused by mutations in the PKD1 or PKD2 genes, and COL4A-associated nephropathy (COL4A-AN), caused by mutations in COL4A3, COL4A4, or COL4A5 genes. These disorders have incomplete penetrance and variable expressivity. It is presently unknown if some of the variability in penetrance can be attributed to polygenic factors. Objective Our objectives were to investigate the interplay of polygenic and monogenic risk of kidney disease and to test the utility of polygenic risk scores for risk stratification in ADPKD and COL4-AN. Design We combined exome/genome sequencing, SNP micro-array, and electronic health record data from the UK Biobank (N = 469,835 participants) and the All of Us (N = 98,622 participants) datasets to test the effects of the genome-wide polygenic score (GPS) for CKD in ADPKD and COL4-AN variant carriers defined by strict variant classification criteria. We used the eMERGE-III electronic CKD phenotype to define cases (estimated glomerular filtration rate (eGFR) < 60 mL/min/1.73m2 or kidney failure) and controls (eGFR > 90 mL/min/1.73m2 in the absence of kidney disease diagnoses). The GPS was tested as a predictor of CKD in qualifying variant carriers and non-carriers using logistic regression after adjustment for age, sex, diabetes, and genetic ancestry. Results The GPS was predictive of CKD in ADPKD variant carriers (ORmeta=2.28 per SD, 95%CI: 1.55–3.37, P = 2.6E-05) as well as non-carriers (ORmeta=1.72 per SD, 95% CI = 1.69–1.76, P < E-300) independent of age, sex, diabetes, and genetic ancestry. Compared to the middle tertile of the GPS distribution for non-carriers, ADPKD variant carriers in the top tertile had a 54-fold increased risk of CKD, while ADPKD variant carriers in the bottom tertile had only a 3-fold increased risk of CKD. Similarly, the GPS was predictive of CKD in both COL4-AN carriers (ORmeta=1.78, 95% CI = 1.22–2.58, P = 2.38E-03) and non-carriers (ORmeta=1.70, 95%CI: 1.68–1.73 P < E-300). The carriers in the top GPS tertile had a 2.5-fold higher risk of CKD while the risk for carriers in the bottom GPS tertile was not different from the average population risk. Conclusions and Relevance Variable penetrance of kidney disease in ADPKD and COL4-AN is partially explained by the differences in polygenic risk profiles. Accounting for polygenic factors improves risk stratification in monogenic kidney disease and may have clinical implications for genetic counseling.

AbstractObjectiveThe aim of this study was to analyze a publicly available sample of rule-based phenotype definitions to characterize and evaluate the variability of logical constructs used.Materials and MethodsA sample of 33 preexisting phenotype definitions used in research that are represented using Fast Healthcare Interoperability Resources and Clinical Quality Language (CQL) was analyzed using automated analysis of the computable representation of the CQL libraries.ResultsMost of the phenotype definitions include narrative descriptions and flowcharts, while few provide pseudocode or executable artifacts. Most use 4 or fewer medical terminologies. The number of codes used ranges from 5 to 6865, and value sets from 1 to 19. We found that the most common expressions used were literal, data, and logical expressions. Aggregate and arithmetic expressions are the least common. Expression depth ranges from 4 to 27.DiscussionDespite the range of conditions, we found that all of the phenotype definitions consisted of logical criteria, representing both clinical and operational logic, and tabular data, consisting of codes from standard terminologies and keywords for natural language processing. The total number and variety of expressions are low, which may be to simplify implementation, or authors may limit complexity due to data availability constraints.ConclusionsThe phenotype definitions analyzed show significant variation in specific logical, arithmetic, and other operators but are all composed of the same high-level components, namely tabular data and logical expressions. A standard representation for phenotype definitions should support these formats and be modular to support localization and shared logic.

Background: Sequencing Mendelian arrhythmia genes in individuals without an indication for arrhythmia genetic testing can identify carriers of pathogenic or likely pathogenic (P/LP) variants. However, the extent to which these variants are associated with clinically meaningful phenotypes before or after return of variant results is unclear. In addition, the majority of discovered variants are currently classified as variants of uncertain significance, limiting clinical actionability. Methods: The eMERGE-III study (Electronic Medical Records and Genomics Phase III) is a multicenter prospective cohort that included 21 846 participants without previous indication for cardiac genetic testing. Participants were sequenced for 109 Mendelian disease genes, including 10 linked to arrhythmia syndromes. Variant carriers were assessed with electronic health record–derived phenotypes and follow-up clinical examination. Selected variants of uncertain significance (n=50) were characterized in vitro with automated electrophysiology experiments in HEK293 cells. Results: As previously reported, 3.0% of participants had P/LP variants in the 109 genes. Herein, we report 120 participants (0.6%) with P/LP arrhythmia variants. Compared with noncarriers, arrhythmia P/LP carriers had a significantly higher burden of arrhythmia phenotypes in their electronic health records. Fifty-four participants had variant results returned. Nineteen of these 54 participants had inherited arrhythmia syndrome diagnoses (primarily long-QT syndrome), and 12 of these 19 diagnoses were made only after variant results were returned (0.05%). After in vitro functional evaluation of 50 variants of uncertain significance, we reclassified 11 variants: 3 to likely benign and 8 to P/LP. Conclusions: Genome sequencing in a large population without indication for arrhythmia genetic testing identified phenotype-positive carriers of variants in congenital arrhythmia syndrome disease genes. As the genomes of large numbers of people are sequenced, the disease risk from rare variants in arrhythmia genes can be assessed by integrating genomic screening, electronic health record phenotypes, and in vitro functional studies. Registration: URL: https://www.clinicaltrials.gov ; Unique identifier; NCT03394859.

Abstract IntroductionCurrently, one of the commonly used methods for disseminating electronic health record (EHR)-based phenotype algorithms is providing a narrative description of the algorithm logic, often accompanied by flowcharts. A challenge with this mode of dissemination is the potential for under-specification in the algorithm definition, which leads to ambiguity and vagueness. MethodsThis study examines incidents of under-specification that occurred during the implementation of 34 narrative phenotyping algorithms in the electronic Medical Record and Genomics (eMERGE) network. We reviewed the online communication history between algorithm developers and implementers within the Phenotype Knowledge Base (PheKB) platform, where questions could be raised and answered regarding the intended implementation of a phenotype algorithm. ResultsWe developed a taxonomy of under-specification categories via an iterative review process between two groups of annotators. Under-specifications that lead to ambiguity and vagueness were consistently found across narrative phenotype algorithms developed by all involved eMERGE sites. Discussion & ConclusionOur findings highlight that under-specification is an impediment to the accuracy and efficiency of the implementation of current narrative phenotyping algorithms, and we propose approaches for mitigating these issues and improved methods for disseminating EHR phenotyping algorithms.

Background: Lp(a) (lipoprotein [a]) levels are higher in individuals of African ancestry (AA) than in individuals of European ancestry (EA). We examined associations of genetically predicted Lp(a) levels with (1) atherosclerotic cardiovascular disease subtypes: coronary heart disease, cerebrovascular disease, peripheral artery disease, and abdominal aortic aneurysm and (2) nonatherosclerotic cardiovascular disease phenotypes, stratified by ancestry. Methods: We performed (1) Mendelian randomization analyses for previously reported cardiovascular associations and (2) Mendelian randomization–phenome-wide association analyses for novel associations. Analyses were stratified by ancestry in electronic Medical Records and Genomics, United Kingdom Biobank, and Million Veteran Program cohorts separately and in a combined cohort of 804 507 EA and 103 580 AA participants. Results: In Mendelian randomization analyses using the combined cohort, a 1-SD genetic increase in Lp(a) level was associated with atherosclerotic cardiovascular disease subtypes in EA—odds ratio and 95% CI for coronary heart disease 1.28 (1.16–1.41); cerebrovascular disease 1.14 (1.07–1.21); peripheral artery disease 1.22 (1.11–1.34); abdominal aortic aneurysm 1.28 (1.17–1.40); in AA, the effect estimate was lower than in EA and nonsignificant for coronary heart disease 1.11 (0.99–1.24) and cerebrovascular disease 1.06 (0.99–1.14) but similar for peripheral artery disease 1.16 (1.01–1.33) and abdominal aortic aneurysm 1.34 (1.11–1.62). In EA, a 1-SD genetic increase in Lp(a) level was associated with aortic valve disorders 1.34 (1.10–1.62), mitral valve disorders 1.18 (1.09–1.27), congestive heart failure 1.12 (1.05–1.19), and chronic kidney disease 1.07 (1.01–1.14). In AA, no significant associations were noted for aortic valve disorders 1.08 (0.94–1.25), mitral valve disorders 1.02 (0.89–1.16), congestive heart failure 1.02 (0.95–1.10), or chronic kidney disease 1.05 (0.99–1.12). Mendelian randomization–phenome-wide association analyses identified novel associations in EA with arterial thromboembolic disease, nonaortic aneurysmal disease, atrial fibrillation, cardiac conduction disorders, and hypertension. Conclusions: Many cardiovascular associations of genetically increased Lp(a) that were significant in EA were not significant in AA. Lp(a) was associated with atherosclerotic cardiovascular disease in four major arterial beds in EA but only with peripheral artery disease and abdominal aortic aneurysm in AA. Additionally, novel cardiovascular associations were detected in EA.

Background: The pathophysiology of lipoprotein (a) [Lp(a)] remains obscure. We conducted a phenome wide analysis study (PheWAS) to identify pleiotropic effects of LPA variants. Methods: Among 51,843 European-ancestry participants (age 58±16 years, 54% female) of the electronic MEdical Records and Genomics (eMERGE) network we assessed whether LPA variants exhibited pleiotropic affects through an electronic health record-based PheWAS. We adjusted significant associations by presence of atherosclerotic cardiovascular disease (ASCVD; defined as a composite of coronary heart disease, cerebrovascular disease, and peripheral artery disease) to determine whether those associations were independent of ASCVD. Results: In PheWAS analysis, we tested 864 phecodes with 36 independent LPA variants. We identified 21 significant associations with non-ASCVD phenotypes including heart failure, aortic valve stenosis, and chronic kidney disease of which 14 were replicated across independent cohorts. The strength of associations between the LPA variant rs10455872 and both heart failure and aortic valve stenosis related phecodes were significantly attenuated after adjustment for ASCVD. However, the association with chronic kidney disease phecode remained independent following adjustment for ASCVD without any attenuation in the estimated odds ratio (Table). Conclusions: LPA genetic variants were associated with aortic valve stenosis, heart failure and chronic kidney disease. The observed association of LPA variant rs10455872 with aortic stenosis and heart failure appear to be partially mediated by ASCVD, while the association with chronic kidney disease appears to be independent of ASCVD. Additional studies are needed to elucidate potential mechanisms by which Lp(a) may contribute to the pathogenesis of non-ASCVD disease phenotypes.

Abstract Electronic health record (EHR) algorithms for defining patient cohorts are commonly shared as free-text descriptions that require human intervention both to interpret and implement. We developed the Phenotype Execution and Modeling Architecture (PhEMA, http://projectphema.org) to author and execute standardized computable phenotype algorithms. With PhEMA, we converted an algorithm for benign prostatic hyperplasia, developed for the electronic Medical Records and Genomics network (eMERGE), into a standards-based computable format. Eight sites (7 within eMERGE) received the computable algorithm, and 6 successfully executed it against local data warehouses and/or i2b2 instances. Blinded random chart review of cases selected by the computable algorithm shows PPV ≥90%, and 3 out of 5 sites had &gt;90% overlap of selected cases when comparing the computable algorithm to their original eMERGE implementation. This case study demonstrates potential use of PhEMA computable representations to automate phenotyping across different EHR systems, but also highlights some ongoing challenges.

Background: Coronary heart disease (CHD) is a leading cause of death globally. Although therapy with statins decreases circulating levels of low-density lipoprotein cholesterol and the incidence of CHD, additional events occur despite statin therapy in some individuals. The genetic determinants of this residual cardiovascular risk remain unknown. Methods: We performed a 2-stage genome-wide association study of CHD events during statin therapy. We first identified 3099 cases who experienced CHD events (defined as acute myocardial infarction or the need for coronary revascularization) during statin therapy and 7681 controls without CHD events during comparable intensity and duration of statin therapy from 4 sites in the Electronic Medical Records and Genomics Network. We then sought replication of candidate variants in another 160 cases and 1112 controls from a fifth Electronic Medical Records and Genomics site, which joined the network after the initial genome-wide association study. Finally, we performed a phenome-wide association study for other traits linked to the most significant locus. Results: The meta-analysis identified 7 single nucleotide polymorphisms at a genome-wide level of significance within the LPA/PLG locus associated with CHD events on statin treatment. The most significant association was for an intronic single nucleotide polymorphism within LPA/PLG (rs10455872; minor allele frequency, 0.069; odds ratio, 1.58; 95% confidence interval, 1.35–1.86; P =2.6×10 − 10 ). In the replication cohort, rs10455872 was also associated with CHD events (odds ratio, 1.71; 95% confidence interval, 1.14–2.57; P =0.009). The association of this single nucleotide polymorphism with CHD events was independent of statin-induced change in low-density lipoprotein cholesterol (odds ratio, 1.62; 95% confidence interval, 1.17–2.24; P =0.004) and persisted in individuals with low-density lipoprotein cholesterol ≤70 mg/dL (odds ratio, 2.43; 95% confidence interval, 1.18–4.75; P =0.015). A phenome-wide association study supported the effect of this region on coronary heart disease and did not identify noncardiovascular phenotypes. Conclusions: Genetic variations at the LPA locus are associated with CHD events during statin therapy independently of the extent of low-density lipoprotein cholesterol lowering. This finding provides support for exploring strategies targeting circulating concentrations of lipoprotein(a) to reduce CHD events in patients receiving statins.