Background Admixture mapping can be leveraged to identify genomic loci associated with disease traits in admixed individuals such as Hispanic Americans. Methods and Results We characterized fine‐scale population structure of participants (n=105 108) in the multiethnic eMERGE‐III (electronic Medical Records and Genomics) cohort, using reference panels from the 1000 Genomes Project and a Native American cohort. We defined “genetically admixed Hispanic Americans” as those having ≥10% Native American genetic ancestry (n=5025), in addition to European (specifically, Iberian) and African genetic ancestry. This population descriptor was subsequently applied to 3 other cohorts: the MESA (Multi‐Ethnic Study of Atherosclerosis; n=1176), WHI (Women’s Health Initiative; n=2873), and the HCHS/SOL (Hispanic Community Health Study/Study of Latinos; n=7804). We inferred local ancestry along the genome in genetically admixed Hispanic Americans in each cohort using RFMix2 with phased haplotypes from SHAPEIT2. Local ancestry proportions were used for admixture mapping of lipid traits (low‐density lipoprotein cholesterol, high‐density lipoprotein cholesterol, total cholesterol, and triglycerides) using extended linear mixed models that incorporated multiple or single ancestries. After accounting for fixed‐effect covariates and random effects, ancestry proportions at the following loci were associated with lipid traits in the HCHS/SOL cohort: 11q12.2 ( P =6.8×10 −7 ) with low‐density lipoprotein cholesterol, 2p23.3 ( P =5.2×10 −7 ) with total cholesterol, and 11q23.3 ( P =7.3×10 −12 ) with triglycerides. The LDL cholesterol and triglyceride loci replicated in at least one independent cohort, whereas the total cholesterol locus showed suggestive evidence for replication in the eMERGE cohort ( P =0.073). Native American (NAMR) ancestry appeared to drive the LDL‐C association, whereas both NAMR and Iberian population in Spain (IBS) ancestries contributed to the triglyceride signal. Conclusions Admixture mapping in genetically admixed Hispanic Americans identified genomic loci influencing lipid traits and the different ancestral origins driving these associations.
Implementing integrated genomic risk assessments for breast cancer: lessons learned from the Electronic Medical Records and Genomics study
Abstract Objectives To implementation an automated multi-institutional pipeline that delivers breast-cancer risk integrated with polygenic risk scores, monogenic variants, family history, and clinical factors, emphasizing operational challenges and their solutions. Materials and Methods A five-stage process was executed at ten sites. Data streams from REDCap surveys, PRS and monogenic reports, and MeTree pedigrees were normalized and forwarded through a REDCap plug-in to the CanRisk API. Results Integrated risk was returned to >10 000 women; 3.6% were ≥25 % lifetime risk and 0.9% carried pathogenic variants. Pipeline generated score aligns well with manual generated ones. Major barriers such as heterogeneous pedigree formats, missing data, edge-case handling, and evolving model versions were identified and resolved through mapping rules, imputations, and iterative testing. Discussion Cross-platform data harmonization and stakeholder alignment were decisive for success. Borderline-risk communication and model-version drift remain open issues. Conclusion Large-scale PRS-integrated breast-cancer risk reporting is feasible but requires robust interoperability standards and iterative governance.
Author response: Risk factors affecting polygenic score performance across diverse cohorts
Apart from ancestry, personal or environmental covariates may contribute to differences in polygenic score (PGS) performance. We analyzed the effects of covariate stratification and interaction on body mass index (BMI) PGS (PGS BMI ) across four cohorts of European (N = 491,111) and African (N = 21,612) ancestry. Stratifying on binary covariates and quintiles for continuous covariates, 18/62 covariates had significant and replicable R 2 differences among strata. Covariates with the largest differences included age, sex, blood lipids, physical activity, and alcohol consumption, with R 2 being nearly double between best- and worst-performing quintiles for certain covariates. Twenty-eight covariates had significant PGS BMI –covariate interaction effects, modifying PGS BMI effects by nearly 20% per standard deviation change. We observed overlap between covariates that had significant R 2 differences among strata and interaction effects – across all covariates, their main effects on BMI were correlated with their maximum R 2 differences and interaction effects (0.56 and 0.58, respectively), suggesting high-PGS BMI individuals have highest R 2 and increase in PGS effect. Using quantile regression, we show the effect of PGS BMI increases as BMI itself increases, and that these differences in effects are directly related to differences in R 2 when stratifying by different covariates. Given significant and replicable evidence for context-specific PGS BMI performance and effects, we investigated ways to increase model performance taking into account nonlinear effects. Machine learning models (neural networks) increased relative model R 2 (mean 23%) across datasets. Finally, creating PGS BMI directly from GxAge genome-wide association studies effects increased relative R 2 by 7.8%. These results demonstrate that certain covariates, especially those most associated with BMI, significantly affect both PGS BMI performance and effects across diverse cohorts and ancestries, and we provide avenues to improve model performance that consider these effects.
Risk factors affecting polygenic score performance across diverse cohorts
Apart from ancestry, personal or environmental covariates may contribute to differences in polygenic score (PGS) performance. We analyzed the effects of covariate stratification and interaction on body mass index (BMI) PGS (PGS BMI ) across four cohorts of European (N = 491,111) and African (N = 21,612) ancestry. Stratifying on binary covariates and quintiles for continuous covariates, 18/62 covariates had significant and replicable R 2 differences among strata. Covariates with the largest differences included age, sex, blood lipids, physical activity, and alcohol consumption, with R 2 being nearly double between best- and worst-performing quintiles for certain covariates. Twenty-eight covariates had significant PGS BMI –covariate interaction effects, modifying PGS BMI effects by nearly 20% per standard deviation change. We observed overlap between covariates that had significant R 2 differences among strata and interaction effects – across all covariates, their main effects on BMI were correlated with their maximum R 2 differences and interaction effects (0.56 and 0.58, respectively), suggesting high-PGS BMI individuals have highest R 2 and increase in PGS effect. Using quantile regression, we show the effect of PGS BMI increases as BMI itself increases, and that these differences in effects are directly related to differences in R 2 when stratifying by different covariates. Given significant and replicable evidence for context-specific PGS BMI performance and effects, we investigated ways to increase model performance taking into account nonlinear effects. Machine learning models (neural networks) increased relative model R 2 (mean 23%) across datasets. Finally, creating PGS BMI directly from GxAge genome-wide association studies effects increased relative R 2 by 7.8%. These results demonstrate that certain covariates, especially those most associated with BMI, significantly affect both PGS BMI performance and effects across diverse cohorts and ancestries, and we provide avenues to improve model performance that consider these effects.
Polygenic risk alters the penetrance of monogenic kidney disease
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.
Characterizing variability of electronic health record-driven phenotype definitions
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.
Quantifying the phenome‐wide disease burden of obesity using electronic health records and genomics
AbstractObjectiveHigh BMI is associated with many comorbidities and mortality. This study aimed to elucidate the overall clinical risk of obesity using a genome‐ and phenome‐wide approach.MethodsThis study performed a phenome‐wide association study of BMI using a clinical cohort of 736,726 adults. This was followed by genetic association studies using two separate cohorts: one consisting of 65,174 adults in the Electronic Medical Records and Genomics (eMERGE) Network and another with 405,432 participants in the UK Biobank.ResultsClass 3 obesity was associated with 433 phenotypes, representing 59.3% of all billing codes in individuals with severe obesity. A genome‐wide polygenic risk score for BMI, accounting for 7.5% of variance in BMI, was associated with 296 clinical diseases, including strong associations with type 2 diabetes, sleep apnea, hypertension, and chronic liver disease. In all three cohorts, 199 phenotypes were associated with class 3 obesity and polygenic risk for obesity, including novel associations such as increased risk of renal failure, venous insufficiency, and gastroesophageal reflux.ConclusionsThis combined genomic and phenomic systematic approach demonstrated that obesity has a strong genetic predisposition and is associated with a considerable burden of disease across all disease classes.
A research agenda to support the development and implementation of genomics-based clinical informatics tools and resources
Abstract Objective The Genomic Medicine Working Group of the National Advisory Council for Human Genome Research virtually hosted its 13th genomic medicine meeting titled “Developing a Clinical Genomic Informatics Research Agenda”. The meeting’s goal was to articulate a research strategy to develop Genomics-based Clinical Informatics Tools and Resources (GCIT) to improve the detection, treatment, and reporting of genetic disorders in clinical settings. Materials and Methods Experts from government agencies, the private sector, and academia in genomic medicine and clinical informatics were invited to address the meeting's goals. Invitees were also asked to complete a survey to assess important considerations needed to develop a genomic-based clinical informatics research strategy. Results Outcomes from the meeting included identifying short-term research needs, such as designing and implementing standards-based interfaces between laboratory information systems and electronic health records, as well as long-term projects, such as identifying and addressing barriers related to the establishment and implementation of genomic data exchange systems that, in turn, the research community could help address. Discussion Discussions centered on identifying gaps and barriers that impede the use of GCIT in genomic medicine. Emergent themes from the meeting included developing an implementation science framework, defining a value proposition for all stakeholders, fostering engagement with patients and partners to develop applications under patient control, promoting the use of relevant clinical workflows in research, and lowering related barriers to regulatory processes. Another key theme was recognizing pervasive biases in data and information systems, algorithms, access, value, and knowledge repositories and identifying ways to resolve them.
Assessing socioeconomic bias in machine learning algorithms in health care: a case study of the HOUSES index
Abstract Objective Artificial intelligence (AI) models may propagate harmful biases in performance and hence negatively affect the underserved. We aimed to assess the degree to which data quality of electronic health records (EHRs) affected by inequities related to low socioeconomic status (SES), results in differential performance of AI models across SES. Materials and Methods This study utilized existing machine learning models for predicting asthma exacerbation in children with asthma. We compared balanced error rate (BER) against different SES levels measured by HOUsing-based SocioEconomic Status measure (HOUSES) index. As a possible mechanism for differential performance, we also compared incompleteness of EHR information relevant to asthma care by SES. Results Asthmatic children with lower SES had larger BER than those with higher SES (eg, ratio = 1.35 for HOUSES Q1 vs Q2–Q4) and had a higher proportion of missing information relevant to asthma care (eg, 41% vs 24% for missing asthma severity and 12% vs 9.8% for undiagnosed asthma despite meeting asthma criteria). Discussion Our study suggests that lower SES is associated with worse predictive model performance. It also highlights the potential role of incomplete EHR data in this differential performance and suggests a way to mitigate this bias. Conclusion The HOUSES index allows AI researchers to assess bias in predictive model performance by SES. Although our case study was based on a small sample size and a single-site study, the study results highlight a potential strategy for identifying bias by using an innovative SES measure.