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.
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.
Familial Hypercholesterolemia in the Electronic Medical Records and Genomics Network: Prevalence, Penetrance, Cardiovascular Risk, and Outcomes After Return of Results
Background: The implications of secondary findings detected in large-scale sequencing projects remain uncertain. We assessed prevalence and penetrance of pathogenic familial hypercholesterolemia (FH) variants, their association with coronary heart disease (CHD), and 1-year outcomes following return of results in phase III of the electronic medical records and genomics network. Methods: Adult participants (n=18 544) at 7 sites were enrolled in a prospective cohort study to assess the clinical impact of returning results from targeted sequencing of 68 actionable genes, including LDLR , APOB , and PCSK9 . FH variant prevalence and penetrance (defined as low-density lipoprotein cholesterol >155 mg/dL) were estimated after excluding participants enrolled on the basis of hypercholesterolemia. Multivariable logistic regression was used to estimate the odds of CHD compared to age- and sex-matched controls without FH-associated variants. Process (eg, referral to a specialist or ordering new tests), intermediate (eg, new diagnosis of FH), and clinical (eg, treatment modification) outcomes within 1 year after return of results were ascertained by electronic health record review. Results: The prevalence of FH-associated pathogenic variants was 1 in 188 (69 of 13,019 unselected participants). Penetrance was 87.5%. The presence of an FH variant was associated with CHD (odds ratio, 3.02 [2.00–4.53]) and premature CHD (odds ratio, 3.68 [2.34–5.78]). At least 1 outcome occurred in 92% of participants; 44% received a new diagnosis of FH and 26% had treatment modified following return of results. Conclusions: In a multisite cohort of electronic health record–linked biobanks, monogenic FH was prevalent, penetrant, and associated with presence of CHD. Nearly half of participants with an FH-associated variant received a new diagnosis of FH and a quarter had treatment modified after return of results. These results highlight the potential utility of sequencing electronic health record–linked biobanks to detect FH.
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.
Clinical, biochemical and genetic characteristics of MOGS-CDG: a rare congenital disorder of glycosylation
Purpose To summarise the clinical, molecular and biochemical phenotype of mannosyl-oligosaccharide glucosidase-related congenital disorders of glycosylation (MOGS-CDG), which presents with variable clinical manifestations, and to analyse which clinical biochemical assay consistently supports diagnosis in individuals with bi-allelic variants in MOGS. Methods Phenotypic characterisation was performed through an international and multicentre collaboration. Genetic testing was done by exome sequencing and targeted arrays. Biochemical assays on serum and urine were performed to delineate the biochemical signature of MOGS-CDG. Results Clinical phenotyping revealed heterogeneity in MOGS-CDG, including neurological, immunological and skeletal phenotypes. Bi-allelic variants in MOGS were identified in 12 individuals from 11 families. The severity in each organ system was variable, without definite genotype correlation. Urine oligosaccharide analysis was consistently abnormal for all affected probands, whereas other biochemical analyses such as serum transferrin analysis was not consistently abnormal. Conclusion The clinical phenotype of MOGS-CDG includes multisystemic involvement with variable severity. Molecular analysis, combined with biochemical testing, is important for diagnosis. In MOGS-CDG, urine oligosaccharide analysis via matrix-assisted laser desorption/ionisation time-of-flight mass spectrometry can be used as a reliable biochemical test for screening and confirmation of disease.
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.
Arrhythmia Variant Associations and Reclassifications in the eMERGE-III Sequencing Study
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.
Generalizability of Polygenic Risk Scores for Breast Cancer Among Women With European, African, and Latinx Ancestry
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.