The Penn Medicine BioBank: Towards a Genomics-Enabled Learning Healthcare System to Accelerate Precision Medicine in a Diverse Population

Big data in sarcoidosis

PURPOSE OF REVIEW

This review provides an overview of recent advancements in sarcoidosis research, focusing on collaborative networks, phenotype characterization, and molecular studies. It highlights the importance of collaborative efforts, phenotype characterization, and the integration of multilevel molecular data for advancing sarcoidosis research and paving the way toward personalized medicine.

RECENT FINDINGS

Sarcoidosis exhibits heterogeneous clinical manifestations influenced by various factors. Efforts to define sarcoidosis endophenotypes show promise, while technological advancements enable extensive molecular data generation. Collaborative networks and biobanks facilitate large-scale studies, enhancing biomarker discovery and therapeutic protocols.

SUMMARY

Sarcoidosis presents a complex challenge due to its unknown cause and heterogeneous clinical manifestations. Collaborative networks, comprehensive phenotype delineation, and the utilization of cutting-edge technologies are essential for advancing our understanding of sarcoidosis biology and developing personalized medicine approaches. Leveraging large-scale epidemiological resources and biobanks and integrating multilevel molecular data offer promising avenues for unraveling the disease's heterogeneity and improving patient outcomes.

The common p.Ile291Val variant of ERLIN1 enhances TM6SF2 function and is associated with protection against MASLD

BACKGROUND

The ERLIN1 p.Ile291Val single-nucleotide polymorphism (rs2862954) is associated with protection from steatotic liver disease (SLD), but effects of this variant on metabolic phenotypes remain uncertain.

METHODS

Metabolic phenotypes and outcomes associated with ERLIN1 p.Ile291Val were analyzed by using a genome-first approach in the UK Biobank (UKB), Penn Medicine BioBank (PMBB), and All of Us cohort.

FINDINGS

ERLIN1 p.Ile291Val carriers exhibited significantly lower serum levels of alanine aminotransferase and aspartate aminotransferase as well as higher levels of triglycerides, low-density lipoprotein cholesterol, Apolipoprotein B, high-density lipoprotein cholesterol, and Apolipoprotein A1 in UKB, and these values were affected by ERLIN1 p.Ile291Val in an allele-dose-dependent manner. Homozygous ERLIN1 p.Ile291Val carriers had a significantly reduced risk of developing metabolic dysfunction-associated SLD (MASLD, adjusted odds ratio [aOR] = 0.92, 95% confidence interval [CI], 0.88-0.96). The protective effect of this variant was enhanced in patients with alcoholic liver disease. Our results were replicated in PMBB and the All of Us cohort. Strikingly, the protective effects of ERLIN1 p.Ile291Val were not apparent in individuals carrying the TM6SF2 p.Glu167Lys variant associated with increased risk of SLD. We analyzed the effects of predicted loss-of-function ERLIN1 variants and found that they had opposite effects, namely reduced plasma lipids, suggesting that ERLIN1 p.Ile291Val may be a gain-of-function variant.

CONCLUSION

Our study contributes to a better understanding of ERLIN1 by investigating a coding variant that has emerged as a potential gain-of-function mutation with protective effects against MASLD development.

Interpreting SNP heritability in admixed populations

SNP heritability is defined as the proportion of phenotypic variance explained by genotyped SNPs and is believed to be a lower bound of heritability ( h 2 ), being equal to it if all causal variants are known. Despite the simple intuition behind , its interpretation and equivalence to h 2 is unclear, particularly in the presence of population structure and assortative mating. It is well known that population structure can lead to inflation in estimates because of confounding due to linkage disequilibrium (LD) or shared environment. Here we use analytical theory and simulations to demonstrate that estimates can be biased in admixed populations, even in the absence of confounding and even if all causal variants are known. This is because admixture generates LD, which contributes to the genetic variance, and therefore to heritability. Genome-wide restricted maximum likelihood (GREML) does not capture this contribution leading to under-or over-estimates of relative to h 2 , depending on the genetic architecture. In contrast, Haseman-Elston (HE) regression exaggerates the LD contribution leading to biases in the opposite direction. For the same reason, GREML and HE estimates of local ancestry heritability are also biased. We describe this bias in and as a function of admixture history and the genetic architecture of the trait and show that it can be recovered under some conditions. We clarify the interpretation of in admixed populations and discuss its implication for genome-wide association studies and polygenic prediction.

Assessing the Utility of a Genotype-Guided Tacrolimus Equation in African American Kidney Transplant Recipients: A Single Institution Retrospective Study

Tacrolimus metabolism is heavily influenced by the CYP3A5 genotype, which varies widely among African Americans (AA). We aimed to assess the performance of a published genotype-informed tacrolimus dosing model in an independent set of adult AA kidney transplant (KTx) recipients. CYP3A5 genotypes were obtained for all AA KTx recipients (n = 232) from 2010 to 2019 who met inclusion criteria at a single transplant center in Philadelphia, Pennsylvania, USA. Medical record data were used to calculate predicted tacrolimus clearance using the published AA KTx dosing equation and two modified iterations. Observed and model-predicted trough levels were compared at 3 days, 3 months, and 6 months post-transplant. The mean prediction error at day 3 post-transplant was 3.05 ng/mL, indicating that the model tended to overpredict the tacrolimus trough. This bias improved over time to 1.36 and 0.78 ng/mL at 3 and 6 months post-transplant, respectively. Mean absolute prediction error-a marker of model precision-improved with time to 2.33 ng/mL at 6 months. Limiting genotype data in the model decreased bias and improved precision. The bias and precision of the published model improved over time and were comparable to studies in previous cohorts. The overprediction observed by the published model may represent overfitting to the initial cohort, possibly limiting generalizability.

Genetic predictors of blood pressure traits are associated with preeclampsia

Preeclampsia, a pregnancy complication characterized by hypertension after 20 gestational weeks, is a major cause of maternal and neonatal morbidity and mortality. Mechanisms leading to preeclampsia are unclear; however, there is evidence of high heritability. We evaluated the association of polygenic scores (PGS) for blood pressure traits and preeclampsia to assess whether there is shared genetic architecture. Non-Hispanic Black and White reproductive age females with pregnancy indications and genotypes were obtained from Vanderbilt University's BioVU, Electronic Medical Records and Genomics network, and Penn Medicine Biobank. Preeclampsia was defined by ICD codes. Summary statistics for diastolic blood pressure (DBP), systolic blood pressure (SBP), and pulse pressure (PP) PGS were acquired from Giri et al. Associations between preeclampsia and each PGS were evaluated separately by race and data source before subsequent meta-analysis. Ten-fold cross validation was used for prediction modeling. In 3504 Black and 5009 White included individuals, the rate of preeclampsia was 15.49%. In cross-ancestry meta-analysis, all PGSs were associated with preeclampsia (ORDBP = 1.10, 95% CI 1.02-1.17, p = 7.68 × 10-3; ORSBP = 1.16, 95% CI 1.09-1.23, p = 2.23 × 10-6; ORPP = 1.14, 95% CI 1.07-1.27, p = 9.86 × 10-5). Addition of PGSs to clinical prediction models did not improve predictive performance. Genetic factors contributing to blood pressure regulation in the general population also predispose to preeclampsia.

Population Performance and Individual Agreement of Coronary Artery Disease Polygenic Risk Scores

Importance

Polygenic risk scores (PRSs) for coronary artery disease (CAD) are a growing clinical and commercial reality. Whether existing scores provide similar individual-level assessments of disease liability is a critical consideration for clinical implementation that remains uncharacterized.

Objective

Characterize the reliability of CAD PRSs that perform equivalently at the population level at predicting individual-level risk.

Design

Cross-sectional Study.

Setting

All of Us Research Program (AOU), Penn Medicine Biobank (PMBB), and UCLA ATLAS Precision Health Biobank.

Participants

Volunteers of diverse genetic backgrounds enrolled in AOU, PMBB, and UCLA with available electronic health record and genotyping data.

Exposures

Polygenic risk for CAD from previously published PRSs and new PRSs developed separately from the testing cohorts.

Main Outcomes and Measures

Sets of CAD PRSs that perform population prediction equivalently were identified by comparing calibration and discrimination (Brier score and AUROC) of generalized linear models of prevalent CAD using Bayesian analysis of variance. Among equivalently performing scores, individual-level agreement between risk estimates was tested with intraclass correlation (ICC) and Light's Kappa, measures of inter-rater reliability.

Results

50 PRSs were calculated for 171,095 AOU participants. When included in a model of prevalent CAD, 48 scores had practically equivalent Brier scores and AUROCs (region of practical equivalence = 0.02). Across these scores, 84% of participants had at least one score in both the top and bottom risk quintile. Continuous agreement of individual risk predictions from the 48 scores was poor, with an ICC of 0.351 (95% CI; 0.349, 0.352). Agreement between two statistically equivalent scores was moderate, with an ICC of 0.649 (95% CI; 0.646, 0.652). Light's Kappa, used to evaluate consistency of assignment to high-risk thresholds, did not exceed 0.56 (interpreted as 'fair') across statistically and practically equivalent scores. Repeating the analysis among 41,193 PMBB and 50,748 UCLA participants yielded different sets of statistically and practically equivalent scores which also lacked strong individual agreement.

Conclusions and Relevance

Across three diverse biobanks, CAD PRSs that performed equivalently at the population level produced unreliable individual risk estimates. Approaches to clinical implementation of CAD PRSs must consider the potential for discordant individual risk estimates from otherwise indistinguishable scores.

CT imaging-derived phenotypes for abdominal muscle and their association with age and sex in a medical biobank

The study of muscle mass as an imaging-derived phenotype (IDP) may yield new insights into determining the normal and pathologic variations in muscle mass in the population. This can be done by determining 3D abdominal muscle mass from 12 distinct abdominal muscle regions and groups using computed tomography (CT) in a racially diverse medical biobank. To develop a fully automatic technique for assessment of CT abdominal muscle IDPs and preliminarily determine abdominal muscle IDP variations with age and sex in a clinically and racially diverse medical biobank. This retrospective study was conducted using the Penn Medicine BioBank (PMBB), a research protocol that recruits adult participants during outpatient visits at hospitals in the Penn Medicine network. We developed a deep residual U-Net (ResUNet) to segment 12 abdominal muscle groups including the left and right psoas, quadratus lumborum, erector spinae, gluteus medius, rectus abdominis, and lateral abdominals. 110 CT studies were randomly selected for training, validation, and testing. 44 of the 110 CT studies were selected to enrich the dataset with representative cases of intra-abdominal and abdominal wall pathology. The studies were divided into non-overlapping training, validation and testing sets. Model performance was evaluated using the Sørensen-Dice coefficient. Volumes of individual muscle groups were plotted to distribution curves. To investigate associations between muscle IDPs, age, and sex, deep learning model segmentations were performed on a larger abdominal CT dataset from PMBB consisting of 295 studies. Multivariable models were used to determine relationships between muscle mass, age and sex. The model's performance (Dice scores) on the test data was the following: psoas: 0.85 ± 0.12, quadratus lumborum: 0.72 ± 0.14, erector spinae: 0.92 ± 0.07, gluteus medius: 0.90 ± 0.08, rectus abdominis: 0.85 ± 0.08, lateral abdominals: 0.85 ± 0.09. The average Dice score across all muscle groups was 0.86 ± 0.11. Average total muscle mass for females was 2041 ± 560.7 g with a high of 2256 ± 560.1 g (41-50 year old cohort) and a change of - 0.96 g/year, declining to an average mass of 1579 ± 408.8 g (81-100 year old cohort). Average total muscle mass for males was 3086 ± 769.1 g with a high of 3385 ± 819.3 g (51-60 year old cohort) and a change of - 1.73 g/year, declining to an average mass of 2629 ± 536.7 g (81-100 year old cohort). Quadratus lumborum was most highly correlated with age for both sexes (correlation coefficient of - 0.5). Gluteus medius mass in females was positively correlated with age with a coefficient of 0.22. These preliminary findings show that our CNN can automate detailed abdominal muscle volume measurement. Unlike prior efforts, this technique provides 3D muscle segmentations of individual muscles. This technique will dramatically impact sarcopenia diagnosis and research, elucidating its clinical and public health implications. Our results suggest a peak age range for muscle mass and an expected rate of decline, both of which vary between genders. Future goals are to investigate genetic variants for sarcopenia and malnutrition, while describing genotype-phenotype associations of muscle mass in healthy humans using imaging-derived phenotypes. It is feasible to obtain 3D abdominal muscle IDPs with high accuracy from patients in a medical biobank using fully automated machine learning methods. Abdominal muscle IDPs showed significant variations in lean mass by age and sex. In the future, this tool can be leveraged to perform a genome-wide association study across the medical biobank and determine genetic variants associated with early or accelerated muscle wasting.

Optimising the use of electronic medical records for large scale research in psychiatry

The explosion and abundance of digital data could facilitate large-scale research for psychiatry and mental health. Research using so-called "real world data"-such as electronic medical/health records-can be resource-efficient, facilitate rapid hypothesis generation and testing, complement existing evidence (e.g. from trials and evidence-synthesis) and may enable a route to translate evidence into clinically effective, outcomes-driven care for patient populations that may be under-represented. However, the interpretation and processing of real-world data sources is complex because the clinically important 'signal' is often contained in both structured and unstructured (narrative or "free-text") data. Techniques for extracting meaningful information (signal) from unstructured text exist and have advanced the re-use of routinely collected clinical data, but these techniques require cautious evaluation. In this paper, we survey the opportunities, risks and progress made in the use of electronic medical record (real-world) data for psychiatric research.

Multi-ancestry meta-analysis of tobacco use disorder identifies 461 potential risk genes and reveals associations with multiple health outcomes

Tobacco use disorder (TUD) is the most prevalent substance use disorder in the world. Genetic factors influence smoking behaviours and although strides have been made using genome-wide association studies to identify risk variants, most variants identified have been for nicotine consumption, rather than TUD. Here we leveraged four US biobanks to perform a multi-ancestral meta-analysis of TUD (derived via electronic health records) in 653,790 individuals (495,005 European, 114,420 African American and 44,365 Latin American) and data from UK Biobank (ncombined = 898,680). We identified 88 independent risk loci; integration with functional genomic tools uncovered 461 potential risk genes, primarily expressed in the brain. TUD was genetically correlated with smoking and psychiatric traits from traditionally ascertained cohorts, externalizing behaviours in children and hundreds of medical outcomes, including HIV infection, heart disease and pain. This work furthers our biological understanding of TUD and establishes electronic health records as a source of phenotypic information for studying the genetics of TUD.

Accelerating Genome- and Phenome-Wide Association Studies using GPUs - A case study using data from the Million Veteran Program

The expansion of biobanks has significantly propelled genomic discoveries yet the sheer scale of data within these repositories poses formidable computational hurdles, particularly in handling extensive matrix operations required by prevailing statistical frameworks. In this work, we introduce computational optimizations to the SAIGE (Scalable and Accurate Implementation of Generalized Mixed Model) algorithm, notably employing a GPU-based distributed computing approach to tackle these challenges. We applied these optimizations to conduct a large-scale genome-wide association study (GWAS) across 2,068 phenotypes derived from electronic health records of 635,969 diverse participants from the Veterans Affairs (VA) Million Veteran Program (MVP). Our strategies enabled scaling up the analysis to over 6,000 nodes on the Department of Energy (DOE) Oak Ridge Leadership Computing Facility (OLCF) Summit High-Performance Computer (HPC), resulting in a 20-fold acceleration compared to the baseline model. We also provide a Docker container with our optimizations that was successfully used on multiple cloud infrastructures on UK Biobank and All of Us datasets where we showed significant time and cost benefits over the baseline SAIGE model.

Three-dimensional aortic geometry: clinical correlates, prognostic value and genetic architecture

Aortic structure and function impact cardiovascular health through multiple mechanisms. Aortic structural degeneration increases left ventricular afterload, pulse pressure and promotes target organ damage. Despite the impact of aortic structure on cardiovascular health, aortic 3D-geometry has yet to be comprehensively assessed. Using a convolutional neural network (U-Net) combined with morphological operations, we quantified aortic 3D-geometric phenotypes (AGPs) from 53,612 participants in the UK Biobank and 8,066 participants in the Penn Medicine Biobank. AGPs reflective of structural aortic degeneration, characterized by arch unfolding, descending aortic lengthening and luminal dilation exhibited cross-sectional associations with hypertension and cardiac diseases, and were predictive for new-onset hypertension, heart failure, cardiomyopathy, and atrial fibrillation. We identified 237 novel genetic loci associated with 3D-AGPs. Fibrillin-2 gene polymorphisms were identified as key determinants of aortic arch-3D structure. Mendelian randomization identified putative causal effects of aortic geometry on the risk of chronic kidney disease and stroke.

Combining Transdiagnostic and Disorder-Level GWAS Enhances Precision of Psychiatric Genetic Risk Profiles in a Multi-Ancestry Sample

The etiology of substance use disorders (SUDs) and psychiatric disorders reflects a combination of both transdiagnostic (i.e., common) and disorder-level (i.e., independent) genetic risk factors. We applied genomic structural equation modeling to examine these genetic factors across SUDs, psychotic, mood, and anxiety disorders using genome-wide association studies (GWAS) of European- (EUR) and African-ancestry (AFR) individuals. In EUR individuals, transdiagnostic genetic factors represented SUDs (143 lead single nucleotide polymorphisms [SNPs]), psychotic (162 lead SNPs), and mood/anxiety disorders (112 lead SNPs). We identified two novel SNPs for mood/anxiety disorders that have probable regulatory roles on FOXP1, NECTIN3, and BTLA genes. In AFR individuals, genetic factors represented SUDs (1 lead SNP) and psychiatric disorders (no significant SNPs). The SUD factor lead SNP, although previously significant in EUR- and cross-ancestry GWAS, is a novel finding in AFR individuals. Shared genetic variance accounted for overlap between SUDs and their psychiatric comorbidities, with second-order GWAS identifying up to 12 SNPs not significantly associated with either first-order factor in EUR individuals. Finally, common and independent genetic effects showed different associations with psychiatric, sociodemographic, and medical phenotypes. For example, the independent components of schizophrenia and bipolar disorder had distinct associations with affective and risk-taking behaviors, and phenome-wide association studies identified medical conditions associated with tobacco use disorder independent of the broader SUDs factor. Thus, combining transdiagnostic and disorder-level genetic approaches can improve our understanding of co-occurring conditions and increase the specificity of genetic discovery, which is critical for psychiatric disorders that demonstrate considerable symptom and etiological overlap.

The 14th International Podocyte Conference 2023: from podocyte biology to glomerular medicine

The 14th International Podocyte Conference took place in Philadelphia, Pennsylvania, USA from May 23 to 26, 2023. It commenced with an early-career researchers' meeting on May 23, providing young scientists with a platform to present and discuss their research findings. Throughout the main conference, 29 speakers across 9 sessions shared their insights on podocyte biology, glomerular medicine, novel technologic advancements, and translational approaches. Additionally, the event featured 3 keynote lectures addressing engineered chimeric antigen receptor T cell- and mRNA-based therapies and the use of biobanks for enhanced disease comprehension. Furthermore, 4 brief oral abstract sessions allowed scientists to present their findings to a broad audience. The program also included a panel discussion addressing the challenges of conducting human research within the American Black community. Remarkably, after a 5-year hiatus from in-person conferences, the 14th International Podocyte Conference successfully convened scientists from around the globe, fostering the presentation and discussion of crucial research findings, as summarized in this review. Furthermore, to ensure continuous and sustainable education, research, translation, and trial medicine related to podocyte and glomerular diseases for the benefit of patients, the International Society of Glomerular Disease was officially launched during the conference.

Enhancing Genetic Association Power in Endometriosis through Unsupervised Clustering of Clinical Subtypes Identified from Electronic Health Records

Background

Endometriosis affects 10% of reproductive-age women, and yet, it goes undiagnosed for 3.6 years on average after symptoms onset. Despite large GWAS meta-analyses (N > 750,000), only a few dozen causal loci have been identified. We hypothesized that the challenges in identifying causal genes for endometriosis stem from heterogeneity across clinical and biological factors underlying endometriosis diagnosis.

Methods

We extracted known endometriosis risk factors, symptoms, and concomitant conditions from the Penn Medicine Biobank (PMBB) and performed unsupervised spectral clustering on 4,078 women with endometriosis. The 5 clusters were characterized by utilizing additional electronic health record (EHR) variables, such as endometriosis-related comorbidities and confirmed surgical phenotypes. From four EHR-linked genetic datasets, PMBB, eMERGE, AOU, and UKBB, we extracted lead variants and tag variants 39 known endometriosis loci for association testing. We meta-analyzed ancestry-stratified case/control tests for each locus and cluster in addition to a positive control (Total N endometriosis cases = 10,108).

Results

We have designated the five subtype clusters as pain comorbidities, uterine disorders, pregnancy complications, cardiometabolic comorbidities, and EHR-asymptomatic based on enriched features from each group. One locus, RNLS , surpassed the genome-wide significant threshold in the positive control. Thirteen more loci reached a Bonferroni threshold of 1.3 x 10 -3 (0.05 / 39) in the positive control. The cluster-stratified tests yielded more significant associations than the positive control for anywhere from 5 to 15 loci depending on the cluster. Bonferroni significant loci were identified for four out of five clusters, including WNT4 and GREB1 for the uterine disorders cluster, RNLS for the cardiometabolic cluster, FSHB for the pregnancy complications cluster, and SYNE1 and CDKN2B-AS1 for the EHR-asymptomatic cluster. This study enhances our understanding of the clinical presentation patterns of endometriosis subtypes, showcasing the innovative approach employed to investigate this complex disease.

Distinct genomic and immunologic tumor evolution in germline TP53-driven breast cancers

Pathogenic germline TP53 alterations cause Li-Fraumeni Syndrome (LFS), and breast cancer is the most common cancer in LFS females. We performed first of its kind multimodal analysis of LFS breast cancer (LFS-BC) compared to sporadic premenopausal BC. Nearly all LFS-BC underwent biallelic loss of TP53 with no recurrent oncogenic variants except ERBB2 (HER2) amplification. Compared to sporadic BC, in situ and invasive LFS-BC exhibited a high burden of short amplified aneuploid segments (SAAS). Pro-apoptotic p53 target genes BAX and TP53I3 failed to be up-regulated in LFS-BC as was seen in sporadic BC compared to normal breast tissue. LFS-BC had lower CD8+ T-cell infiltration compared to sporadic BC yet higher levels of proliferating cytotoxic T-cells. Within LFS-BC, progression from in situ to invasive BC was marked by an increase in chromosomal instability with a decrease in proliferating cytotoxic T-cells. Our study uncovers critical events in mutant p53-driven tumorigenesis in breast tissue.

A multi-ancestry genetic study of pain intensity in 598,339 veterans

Chronic pain is a common problem, with more than one-fifth of adult Americans reporting pain daily or on most days. It adversely affects the quality of life and imposes substantial personal and economic costs. Efforts to treat chronic pain using opioids had a central role in precipitating the opioid crisis. Despite an estimated heritability of 25-50%, the genetic architecture of chronic pain is not well-characterized, in part because studies have largely been limited to samples of European ancestry. To help address this knowledge gap, we conducted a cross-ancestry meta-analysis of pain intensity in 598,339 participants in the Million Veteran Program, which identified 126 independent genetic loci, 69 of which are new. Pain intensity was genetically correlated with other pain phenotypes, level of substance use and substance use disorders, other psychiatric traits, education level and cognitive traits. Integration of the genome-wide association studies findings with functional genomics data shows enrichment for putatively causal genes (n = 142) and proteins (n = 14) expressed in brain tissues, specifically in GABAergic neurons. Drug repurposing analysis identified anticonvulsants, β-blockers and calcium-channel blockers, among other drug groups, as having potential analgesic effects. Our results provide insights into key molecular contributors to the experience of pain and highlight attractive drug targets.

Universal Exome Sequencing in Critically Ill Adults: A Diagnostic Yield of 25% and Race-Based Disparities in Access to Genetic Testing

Numerous studies have underscored the diagnostic and therapeutic potential of exome or genome sequencing in critically ill pediatric populations. However, an equivalent investigation in critically ill adults remains conspicuously absent. We retrospectively analyzed whole exome sequencing (WES) data available through the PennMedicine Biobank (PMBB) from all 365 young adult patients, aged 18-40 years, with intensive care unit (ICU) admissions at the University of Pennsylvania Health System who met inclusion criteria for our study. For each participant, two Medical Genetics and Internal Medicine-trained clinicians reviewed WES reports and patient charts for variant classification, result interpretation, and identification of genetic diagnoses related to their critical illness. Of the 365 individuals in our study, 90 (24.7%) were found to have clearly diagnostic results on WES; an additional 40 (11.0%) had a suspicious variant of uncertain significance (VUS) identified; and an additional 16 (4.4%) had a medically actionable incidental finding. The diagnostic rate of exome sequencing did not decrease with increasing patient age. Affected genes were primarily involved in cardiac function (18.8%), vascular health (16.7%), cancer (16.7%), and pulmonary disease (11.5%). Only half of all diagnostic findings were known and documented in the patient chart at the time of ICU admission. Significant disparities emerged in subgroup analysis by EHR-reported race, with genetic diagnoses known/documented for 63.5% of White patients at the time of ICU admission but only for 28.6% of Black or Hispanic patients. There was a trend towards patients with undocumented genetic diagnoses having a 66% increased mortality rate, making these race-based disparities in genetic diagnosis even more concerning. Altogether, universal exome sequencing in ICU-admitted adult patients was found to yield a new definitive diagnosis in 11.2% of patients. Of these diagnoses, 76.6% conferred specific care-altering medical management recommendations. Our study suggests that the diagnostic utility of exome sequencing in critically ill young adults is similar to that observed in neonatal and pediatric populations and is age-independent. The high diagnostic rate and striking race-based disparities we find in genetic diagnoses argue for broad and universal approaches to genetic testing for critically ill adults. The widespread implementation of comprehensive genetic sequencing in the adult population promises to enhance medical care for all individuals and holds the potential to rectify disparities in genetic testing referrals, ultimately promoting more equitable healthcare delivery.

Association between genetic risk and adherence to healthy lifestyle for developing age-related hearing loss

BACKGROUND

Previous studies have shown that lifestyle/environmental factors could accelerate the development of age-related hearing loss (ARHL). However, there has not yet been a study investigating the joint association among genetics, lifestyle/environmental factors, and adherence to healthy lifestyle for risk of ARHL. We aimed to assess the association between ARHL genetic variants, lifestyle/environmental factors, and adherence to healthy lifestyle as pertains to risk of ARHL.

METHODS

This case-control study included 376,464 European individuals aged 40 to 69 years, enrolled between 2006 and 2010 in the UK Biobank (UKBB). As a replication set, we also included a total of 26,523 individuals considered of European ancestry and 9834 individuals considered of African-American ancestry through the Penn Medicine Biobank (PMBB). The polygenic risk score (PRS) for ARHL was derived from a sensorineural hearing loss genome-wide association study from the FinnGen Consortium and categorized as low, intermediate, high, and very high. We selected lifestyle/environmental factors that have been previously studied in association with hearing loss. A composite healthy lifestyle score was determined using seven selected lifestyle behaviors and one environmental factor.

RESULTS

Of the 376,464 participants, 87,066 (23.1%) cases belonged to the ARHL group, and 289,398 (76.9%) individuals comprised the control group in the UKBB. A very high PRS for ARHL had a 49% higher risk of ARHL than those with low PRS (adjusted OR, 1.49; 95% CI, 1.36-1.62; P < .001), which was replicated in the PMBB cohort. A very poor lifestyle was also associated with risk of ARHL (adjusted OR, 3.03; 95% CI, 2.75-3.35; P < .001). These risk factors showed joint effects with the risk of ARHL. Conversely, adherence to healthy lifestyle in relation to hearing mostly attenuated the risk of ARHL even in individuals with very high PRS (adjusted OR, 0.21; 95% CI, 0.09-0.52; P < .001).

CONCLUSIONS

Our findings of this study demonstrated a significant joint association between genetic and lifestyle factors regarding ARHL. In addition, our analysis suggested that lifestyle adherence in individuals with high genetic risk could reduce the risk of ARHL.

Polygenic risk score-based phenome-wide association study of head and neck cancer across two large biobanks

BACKGROUND

Numerous observational studies have highlighted associations of genetic predisposition of head and neck squamous cell carcinoma (HNSCC) with diverse risk factors, but these findings are constrained by design limitations of observational studies. In this study, we utilized a phenome-wide association study (PheWAS) approach, incorporating a polygenic risk score (PRS) derived from a wide array of genomic variants, to systematically investigate phenotypes associated with genetic predisposition to HNSCC. Furthermore, we validated our findings across heterogeneous cohorts, enhancing the robustness and generalizability of our results.

METHODS

We derived PRSs for HNSCC and its subgroups, oropharyngeal cancer and oral cancer, using large-scale genome-wide association study summary statistics from the Genetic Associations and Mechanisms in Oncology Network. We conducted a comprehensive investigation, leveraging genotyping data and electronic health records from 308,492 individuals in the UK Biobank and 38,401 individuals in the Penn Medicine Biobank (PMBB), and subsequently performed PheWAS to elucidate the associations between PRS and a wide spectrum of phenotypes.

RESULTS

We revealed the HNSCC PRS showed significant association with phenotypes related to tobacco use disorder (OR, 1.06; 95% CI, 1.05-1.08; P = 3.50 × 10-15), alcoholism (OR, 1.06; 95% CI, 1.04-1.09; P = 6.14 × 10-9), alcohol-related disorders (OR, 1.08; 95% CI, 1.05-1.11; P = 1.09 × 10-8), emphysema (OR, 1.11; 95% CI, 1.06-1.16; P = 5.48 × 10-6), chronic airway obstruction (OR, 1.05; 95% CI, 1.03-1.07; P = 2.64 × 10-5), and cancer of bronchus (OR, 1.08; 95% CI, 1.04-1.13; P = 4.68 × 10-5). These findings were replicated in the PMBB cohort, and sensitivity analyses, including the exclusion of HNSCC cases and the major histocompatibility complex locus, confirmed the robustness of these associations. Additionally, we identified significant associations between HNSCC PRS and lifestyle factors related to smoking and alcohol consumption.

CONCLUSIONS

The study demonstrated the potential of PRS-based PheWAS in revealing associations between genetic risk factors for HNSCC and various phenotypic traits. The findings emphasized the importance of considering genetic susceptibility in understanding HNSCC and highlighted shared genetic bases between HNSCC and other health conditions and lifestyles.

Genetic Associations of Circulating Cardiovascular Proteins With Gestational Hypertension and Preeclampsia

Importance

Hypertensive disorders of pregnancy (HDPs), including gestational hypertension and preeclampsia, are important contributors to maternal morbidity and mortality worldwide. In addition, women with HDPs face an elevated long-term risk of cardiovascular disease.

Objective

To identify proteins in the circulation associated with HDPs.

Design, Setting, and Participants

Two-sample mendelian randomization (MR) tested the associations of genetic instruments for cardiovascular disease-related proteins with gestational hypertension and preeclampsia. In downstream analyses, a systematic review of observational data was conducted to evaluate the identified proteins' dynamics across gestation in hypertensive vs normotensive pregnancies, and phenome-wide MR analyses were performed to identify potential non-HDP-related effects associated with the prioritized proteins. Genetic association data for cardiovascular disease-related proteins were obtained from the Systematic and Combined Analysis of Olink Proteins (SCALLOP) consortium. Genetic association data for the HDPs were obtained from recent European-ancestry genome-wide association study meta-analyses for gestational hypertension and preeclampsia. Study data were analyzed October 2022 to October 2023.

Exposures

Genetic instruments for 90 candidate proteins implicated in cardiovascular diseases, constructed using cis-protein quantitative trait loci (cis-pQTLs).

Main Outcomes and Measures

Gestational hypertension and preeclampsia.

Results

Genetic association data for cardiovascular disease-related proteins were obtained from 21 758 participants from the SCALLOP consortium. Genetic association data for the HDPs were obtained from 393 238 female individuals (8636 cases and 384 602 controls) for gestational hypertension and 606 903 female individuals (16 032 cases and 590 871 controls) for preeclampsia. Seventy-five of 90 proteins (83.3%) had at least 1 valid cis-pQTL. Of those, 10 proteins (13.3%) were significantly associated with HDPs. Four were robust to sensitivity analyses for gestational hypertension (cluster of differentiation 40, eosinophil cationic protein [ECP], galectin 3, N-terminal pro-brain natriuretic peptide [NT-proBNP]), and 2 were robust for preeclampsia (cystatin B, heat shock protein 27 [HSP27]). Consistent with the MR findings, observational data revealed that lower NT-proBNP (0.76- to 0.88-fold difference vs no HDPs) and higher HSP27 (2.40-fold difference vs no HDPs) levels during the first trimester of pregnancy were associated with increased risk of HDPs, as were higher levels of ECP (1.60-fold difference vs no HDPs). Phenome-wide MR analyses identified 37 unique non-HDP-related protein-disease associations, suggesting potential on-target effects associated with interventions lowering HDP risk through the identified proteins.

Conclusions and Relevance

Study findings suggest genetic associations of 4 cardiovascular disease-related proteins with gestational hypertension and 2 associated with preeclampsia. Future studies are required to test the efficacy of targeting the corresponding pathways to reduce HDP risk.

Novel ancestry-specific primary open-angle glaucoma loci and shared biology with vascular mechanisms and cell proliferation

Primary open-angle glaucoma (POAG), a leading cause of irreversible blindness globally, shows disparity in prevalence and manifestations across ancestries. We perform meta-analysis across 15 biobanks (of the Global Biobank Meta-analysis Initiative) (n = 1,487,441: cases = 26,848) and merge with previous multi-ancestry studies, with the combined dataset representing the largest and most diverse POAG study to date (n = 1,478,037: cases = 46,325) and identify 17 novel significant loci, 5 of which were ancestry specific. Gene-enrichment and transcriptome-wide association analyses implicate vascular and cancer genes, a fifth of which are primary ciliary related. We perform an extensive statistical analysis of SIX6 and CDKN2B-AS1 loci in human GTEx data and across large electronic health records showing interaction between SIX6 gene and causal variants in the chr9p21.3 locus, with expression effect on CDKN2A/B. Our results suggest that some POAG risk variants may be ancestry specific, sex specific, or both, and support the contribution of genes involved in programmed cell death in POAG pathogenesis.

A multi-cohort genome-wide association study in African ancestry individuals reveals risk loci for primary open-angle glaucoma

Primary open-angle glaucoma (POAG), the leading cause of irreversible blindness worldwide, disproportionately affects individuals of African ancestry. We conducted a genome-wide association study (GWAS) for POAG in 11,275 individuals of African ancestry (6,003 cases; 5,272 controls). We detected 46 risk loci associated with POAG at genome-wide significance. Replication and post-GWAS analyses, including functionally informed fine-mapping, multiple trait co-localization, and in silico validation, implicated two previously undescribed variants (rs1666698 mapping to DBF4P2; rs34957764 mapping to ROCK1P1) and one previously associated variant (rs11824032 mapping to ARHGEF12) as likely causal. For individuals of African ancestry, a polygenic risk score (PRS) for POAG from our mega-analysis (African ancestry individuals) outperformed a PRS from summary statistics of a much larger GWAS derived from European ancestry individuals. This study quantifies the genetic architecture similarities and differences between African and non-African ancestry populations for this blinding disease.

Titin-Truncating variants Predispose to Dilated Cardiomyopathy in Diverse Populations

Importance

The effect of high percentage spliced in (hiPSI) TTN truncating variants (TTNtvs) on risk of dilated cardiomyopathy (DCM) has historically been studied among population subgroups defined by genetic similarity to European reference populations. This has raised questions about the effect of TTNtvs in diverse populations, especially among individuals genetically similar to African reference populations.

Objective

To determine the effect of TTNtvs on risk of DCM in diverse population as measured by genetic distance (GD) in principal component (PC) space.

Design

Cohort study.

Setting

Penn Medicine Biobank (PMBB) is a large, diverse biobank.

Participants

Participants were recruited from across the Penn Medicine healthcare system and volunteered to have their electronic health records linked to biospecimen data including DNA which has undergone whole exome sequencing.

Main Outcomes and Measures

Risk of DCM among individuals carrying a hiPSI TTNtv.

Results

Carrying a hiPSI TTNtv was associated with DCM among PMBB participants across a range of GD deciles from the 1000G European centroid; the effect estimates ranged from odds ratio (OR) = 3.29 (95% confidence interval [CI] 1.26 to 8.56) to OR = 9.39 (95% CI 3.82 to 23.13). When individuals were assigned to population subgroups based on genetic similarity to the 1000G reference populations, hiPSI TTNtvs conferred significant risk of DCM among those genetically similar to the 1000G European reference population (OR = 7.55, 95% CI 4.99 to 11.42, P<0.001) and individuals genetically similar to the 1000G African reference population (OR 3.50, 95% CI 1.48 to 8.24, P=0.004).

Conclusions and Relevance

TTNtvs are associated with increased risk of DCM among a diverse cohort. There is no significant difference in effect of TTNtvs on DCM risk across deciles of GD from the 1000G European centroid, suggesting genetic background should not be considered when screening individuals for titin-related DCM.

Penn Healthy Diet survey: pilot validation and scoring

Though diet quality is widely recognised as linked to risk of chronic disease, health systems have been challenged to find a user-friendly, efficient way to obtain information about diet. The Penn Healthy Diet (PHD) survey was designed to fill this void. The purposes of this pilot project were to assess the patient experience with the PHD, to validate the accuracy of the PHD against related items in a diet recall and to explore scoring algorithms with relationship to the Healthy Eating Index (HEI)-2015 computed from the recall data. A convenience sample of participants in the Penn Health BioBank was surveyed with the PHD, the Automated Self-Administered 24-hour recall (ASA24) and experience questions. Kappa scores and Spearman correlations were used to compare related questions in the PHD to the ASA24. Numerical scoring, regression tree and weighted regressions were computed for scoring. Participants assessed the PHD as easy to use and were willing to repeat the survey at least annually. The three scoring algorithms were strongly associated with HEI-2015 scores using National Health and Nutrition Examination Survey 2017-2018 data from which the PHD was developed and moderately associated with the pilot replication data. The PHD is acceptable to participants and at least moderately correlated with the HEI-2015. Further validation in a larger sample will enable the selection of the strongest scoring approach.

Genome-first approach of the prevalence and cancer phenotypes of pathogenic or likely pathogenic germline TP53 variants

Pathogenic or likely pathogenic (P/LP) germline TP53 variants are the primary cause of Li-Fraumeni syndrome (LFS), a hereditary cancer predisposition disorder characterized by early-onset cancers. The population prevalence of P/LP germline TP53 variants is estimated to be approximately one in every 3,500 to 20,000 individuals. However, these estimates are likely impacted by ascertainment biases and lack of clinical and genetic data to account for potential confounding factors, such as clonal hematopoiesis. Genome-first approaches of cohorts linked to phenotype data can further refine these estimates by identifying individuals with variants of interest and then assessing their phenotypes. This study evaluated P/LP germline (variant allele fraction ≥30%) TP53 variants in three cohorts: UK Biobank (UKB, n = 200,590), Geisinger (n = 170,503), and Penn Medicine Biobank (PMBB, n = 43,731). A total of 109 individuals were identified with P/LP germline TP53 variants across the three databases. The TP53 p.R181H variant was the most frequently identified (9 of 109 individuals, 8%). A total of 110 cancers, including 47 hematologic cancers (47 of 110, 43%), were reported in 71 individuals. The prevalence of P/LP germline TP53 variants was conservatively estimated as 1:10,439 in UKB, 1:3,790 in Geisinger, and 1:2,983 in PMBB. These estimates were calculated after excluding related individuals and accounting for the potential impact of clonal hematopoiesis by excluding heterozygotes who ever developed a hematologic cancer. These varying estimates likely reflect intrinsic selection biases of each database, such as healthcare or population-based contexts. Prospective studies of diverse, young cohorts are required to better understand the population prevalence of germline TP53 variants and their associated cancer penetrance.

Polygenic risk scores for cardiometabolic traits demonstrate importance of ancestry for predictive precision medicine

Polygenic risk scores (PRS) have predominantly been derived from genome-wide association studies (GWAS) conducted in European ancestry (EUR) individuals. In this study, we present an in-depth evaluation of PRS based on multi-ancestry GWAS for five cardiometabolic phenotypes in the Penn Medicine BioBank (PMBB) followed by a phenome-wide association study (PheWAS). We examine the PRS performance across all individuals and separately in African ancestry (AFR) and EUR ancestry groups. For AFR individuals, PRS derived using the multi-ancestry LD panel showed a higher effect size for four out of five PRSs (DBP, SBP, T2D, and BMI) than those derived from the AFR LD panel. In contrast, for EUR individuals, the multi-ancestry LD panel PRS demonstrated a higher effect size for two out of five PRSs (SBP and T2D) compared to the EUR LD panel. These findings underscore the potential benefits of utilizing a multi-ancestry LD panel for PRS derivation in diverse genetic backgrounds and demonstrate overall robustness in all individuals. Our results also revealed significant associations between PRS and various phenotypic categories. For instance, CAD PRS was linked with 18 phenotypes in AFR and 82 in EUR, while T2D PRS correlated with 84 phenotypes in AFR and 78 in EUR. Notably, associations like hyperlipidemia, renal failure, atrial fibrillation, coronary atherosclerosis, obesity, and hypertension were observed across different PRSs in both AFR and EUR groups, with varying effect sizes and significance levels. However, in AFR individuals, the strength and number of PRS associations with other phenotypes were generally reduced compared to EUR individuals. Our study underscores the need for future research to prioritize 1) conducting GWAS in diverse ancestry groups and 2) creating a cosmopolitan PRS methodology that is universally applicable across all genetic backgrounds. Such advances will foster a more equitable and personalized approach to precision medicine.

Lymphocyte Count Derived Polygenic Score and Interindividual Variability in CD4 T-cell Recovery in Response to Antiretroviral Therapy

Access to safe and effective antiretroviral therapy (ART) is a cornerstone in the global response to the HIV pandemic. Among people living with HIV, there is considerable interindividual variability in absolute CD4 T-cell recovery following initiation of virally suppressive ART. The contribution of host genetics to this variability is not well understood. We explored the contribution of a polygenic score which was derived from large, publicly available summary statistics for absolute lymphocyte count from individuals in the general population (PGSlymph) due to a lack of publicly available summary statistics for CD4 T-cell count. We explored associations with baseline CD4 T-cell count prior to ART initiation (n=4959) and change from baseline to week 48 on ART (n=3274) among treatment-naïve participants in prospective, randomized ART studies of the AIDS Clinical Trials Group. We separately examined an African-ancestry-derived and a European-ancestry-derived PGSlymph, and evaluated their performance across all participants, and also in the African and European ancestral groups separately. Multivariate models that included PGSlymph, baseline plasma HIV-1 RNA, age, sex, and 15 principal components (PCs) of genetic similarity explained ∼26-27% of variability in baseline CD4 T-cell count, but PGSlymph accounted for <1% of this variability. Models that also included baseline CD4 T-cell count explained ∼7-9% of variability in CD4 T-cell count increase on ART, but PGSlymph accounted for <1% of this variability. In univariate analyses, PGSlymph was not significantly associated with baseline or change in CD4 T-cell count. Among individuals of African ancestry, the African PGSlymph term in the multivariate model was significantly associated with change in CD4 T-cell count while not significant in the univariate model. When applied to lymphocyte count in a general medical biobank population (Penn Medicine BioBank), PGSlymph explained ∼6-10% of variability in multivariate models (including age, sex, and PCs) but only ∼1% in univariate models. In summary, a lymphocyte count PGS derived from the general population was not consistently associated with CD4 T-cell recovery on ART. Nonetheless, adjusting for clinical covariates is quite important when estimating such polygenic effects.

Genomic ascertainment of primary central nervous system cancers in adolescents and young adults

Genomic ascertainment is the inversion of the traditional phenotype-first approach; with a "genome-first" approach, a cohort linked to electronic health records (EHR) undergoes germline sequencing (array, panel, exome, and genome) and deleterious variation of interest in a gene (or set of genes) are identified. Phenotype is then queried from the linked EHR and from call-back investigation and estimates of variant prevalence, disease penetrance, and phenotype can be determined. This should permit a better estimate of the full phenotypic spectrum, severity, and penetrance linked to a deleterious variant. For now, given the modest size, limited EHR, and age of participants in sequenced cohorts, genomic ascertainment approaches to investigate cancer in children and young adults will likely be restricted to descriptive studies and complement traditional phenotype-first work. Another issue is the ascertainment of the cohort itself: Participants need to survive long enough to enroll. Not accounting for this may lead to bias and incorrect estimates of variant prevalence. Adult-focused cohorts with EHR extending back into childhood, linked to cancer registries, and/or studies that permit recontact with participants may facilitate genomic ascertainment in pediatric cancer research. In summary, genomic ascertainment in pediatric primary brain cancer research remains largely untapped and merits further investigation.

Multi-ancestry study of the genetics of problematic alcohol use in over 1 million individuals

Problematic alcohol use (PAU), a trait that combines alcohol use disorder and alcohol-related problems assessed with a questionnaire, is a leading cause of death and morbidity worldwide. Here we conducted a large cross-ancestry meta-analysis of PAU in 1,079,947 individuals (European, N = 903,147; African, N = 122,571; Latin American, N = 38,962; East Asian, N = 13,551; and South Asian, N = 1,716 ancestries). We observed a high degree of cross-ancestral similarity in the genetic architecture of PAU and identified 110 independent risk variants in within- and cross-ancestry analyses. Cross-ancestry fine mapping improved the identification of likely causal variants. Prioritizing genes through gene expression and chromatin interaction in brain tissues identified multiple genes associated with PAU. We identified existing medications for potential pharmacological studies by a computational drug repurposing analysis. Cross-ancestry polygenic risk scores showed better performance of association in independent samples than single-ancestry polygenic risk scores. Genetic correlations between PAU and other traits were observed in multiple ancestries, with other substance use traits having the highest correlations. This study advances our knowledge of the genetic etiology of PAU, and these findings may bring possible clinical applicability of genetics insights-together with neuroscience, biology and data science-closer.

Protocol to evaluate sequential electronic health record-based strategies to increase genetic testing for breast and ovarian cancer risk across diverse patient populations in gynecology practices

BACKGROUND

Germline genetic testing is recommended by the National Comprehensive Cancer Network (NCCN) for individuals including, but not limited to, those with a personal history of ovarian cancer, young-onset (< 50 years) breast cancer, and a family history of ovarian cancer or male breast cancer. Genetic testing is underused overall, and rates are consistently lower among Black and Hispanic populations. Behavioral economics-informed implementation strategies, or nudges, directed towards patients and clinicians may increase the use of this evidence-based clinical practice.

METHODS

Patients meeting eligibility for germline genetic testing for breast and ovarian cancer will be identified using electronic phenotyping algorithms. A pragmatic cohort study will test three sequential strategies to promote genetic testing, two directed at patients and one directed at clinicians, deployed in the electronic health record (EHR) for patients in OB-GYN clinics across a diverse academic medical center. We will use rapid cycle approaches informed by relevant clinician and patient experiences, health equity, and behavioral economics to optimize and de-risk our strategies and methods before trial initiation. Step 1 will send patients messages through the health system patient portal. For non-responders, step 2 will reach out to patients via text message. For non-responders, Step 3 will contact patients' clinicians using a novel "pend and send" tool in the EHR. The primary implementation outcome is engagement with germline genetic testing for breast and ovarian cancer predisposition, defined as a scheduled genetic counseling appointment. Patient data collected through the EHR (e.g., race/ethnicity, geocoded address) will be examined as moderators of the impact of the strategies.

DISCUSSION

This study will be one of the first to sequentially examine the effects of patient- and clinician-directed strategies informed by behavioral economics on engagement with breast and ovarian cancer genetic testing. The pragmatic and sequential design will facilitate a large and diverse patient sample, allow for the assessment of incremental gains from different implementation strategies, and permit the assessment of moderators of strategy effectiveness. The findings may help determine the impact of low-cost, highly transportable implementation strategies that can be integrated into healthcare systems to improve the use of genomic medicine.

TRIAL REGISTRATION

ClinicalTrials.gov. NCT05721326. Registered February 10, 2023. https://www.

CLINICALTRIALS

gov/study/NCT05721326.

Decoding Genetics, Ancestry, and Geospatial Context for Precision Health

Mass General Brigham, an integrated healthcare system based in the Greater Boston area of Massachusetts, annually serves 1.5 million patients. We established the Mass General Brigham Biobank (MGBB), encompassing 142,238 participants, to unravel the intricate relationships among genomic profiles, environmental context, and disease manifestations within clinical practice. In this study, we highlight the impact of ancestral diversity in the MGBB by employing population genetics, geospatial assessment, and association analyses of rare and common genetic variants. The population structures captured by the genetics mirror the sequential immigration to the Greater Boston area throughout American history, highlighting communities tied to shared genetic and environmental factors. Our investigation underscores the potency of unbiased, large-scale analyses in a healthcare-affiliated biobank, elucidating the dynamic interplay across genetics, immigration, structural geospatial factors, and health outcomes in one of the earliest American sites of European colonization.

Phenotypes of undiagnosed adults with actionable OTC and GLA variants

Inherited metabolic disorders (IMDs) are variably expressive, complicating identification of affected individuals. A genotype-first approach can identify individuals at risk for morbidity and mortality from undiagnosed IMDs and can lead to protocols that improve clinical detection, counseling, and management. Using data from 57,340 participants in two hospital biobanks, we assessed the frequency and phenotypes of individuals with pathogenic/likely pathogenic variants (PLPVs) in two IMD genes: GLA, associated with Fabry disease, and OTC, associated with ornithine transcarbamylase deficiency. Approximately 1 in 19,100 participants harbored an undiagnosed PLPV in GLA or OTC. We identified three individuals (2 male, 1 female) with PLPVs in GLA, all of whom were undiagnosed, and three individuals (3 female) with PLPVs in OTC, two of whom were undiagnosed. All three individuals with PLPVs in GLA (100%) had symptoms suggestive of mild Fabry disease, and one individual (14.2%) had an ischemic stroke at age 33, likely indicating the presence of classic disease. No individuals with PLPVs in OTC had documented hyperammonemia despite exposure to catabolic states, but all (100%) had chronic symptoms suggestive of attenuated disease, including mood disorders and migraines. Our findings suggest that GLA and OTC variants identified via a genotype-first approach are of high penetrance and that population screening of these genes can be used to facilitate stepwise phenotyping and appropriate care.

Common- and rare-variant genetic architecture of heart failure across the allele frequency spectrum

Heart failure (HF) is a complex trait, influenced by environmental and genetic factors, that affects over 30 million individuals worldwide. Historically, the genetics of HF have been studied in Mendelian forms of disease, where rare genetic variants have been linked to familial cardiomyopathies. More recently, genome-wide association studies (GWAS) have successfully identified common genetic variants associated with risk of HF. However, the relative importance of genetic variants across the allele-frequency spectrum remains incompletely characterized. Here, we report the results of common- and rare-variant association studies of all-cause heart failure, applying recently developed methods to quantify the heritability of HF attributable to different classes of genetic variation. We combine GWAS data across multiple populations including 207,346 individuals with HF and 2,151,210 without, identifying 176 risk loci at genome-wide significance (p < 5×10-8). Signals at newly identified common-variant loci include coding variants in Mendelian cardiomyopathy genes (MYBPC3, BAG3), as well as regulators of lipoprotein (LPL) and glucose metabolism (GIPR, GLP1R), and are enriched in cardiac, muscle, nerve, and vascular tissues, as well as myocyte and adipocyte cell types. Gene burden studies across three biobanks (PMBB, UKB, AOU) including 27,208 individuals with HF and 349,126 without uncover exome-wide significant (p < 3.15×10-6) associations for HF and rare predicted loss-of-function (pLoF) variants in TTN, MYBPC3, FLNC, and BAG3. Total burden heritability of rare coding variants (2.2%, 95% CI 0.99-3.5%) is highly concentrated in a small set of Mendelian cardiomyopathy genes, and is lower than heritability attributable to common variants (4.3%, 95% CI 3.9-4.7%) which is more diffusely spread throughout the genome. Finally, we demonstrate that common-variant background, in the form of a polygenic risk score (PRS), significantly modifies the risk of HF among carriers of pathogenic truncating variants in the Mendelian cardiomyopathy gene TTN. These findings suggest a significant polygenic component to HF exists that is not captured by current clinical genetic testing.

Identifying genetic loci and phenomic associations of substance use traits: A multi-trait analysis of GWAS (MTAG) study

BACKGROUND AND AIMS

Genome-wide association studies (GWAS) of opioid use disorder (OUD) and cannabis use disorder (CUD) have lagged behind those of alcohol use disorder (AUD) and smoking, where many more loci have been identified. We sought to identify novel loci for substance use traits (SUTs) in both African- (AFR) and European- (EUR) ancestry individuals to enhance our understanding of the traits' genetic architecture.

DESIGN

We used multi-trait analysis of GWAS (MTAG) to analyze four SUTs in EUR subjects (OUD, CUD, AUD and smoking initiation [SMKinitiation]), and three SUTs in AFR subjects (OUD, AUD and smoking trajectory [SMKtrajectory]). We conducted gene-set and protein-protein interaction analyses and calculated polygenic risk scores (PRS) in two independent samples.

SETTING

This study was conducted in the United States.

PARTICIPANTS

A total of 5692 EUR and 4918 AFR individuals in the Yale-Penn sample and 29 054 EUR and 10 265 AFR individuals in the Penn Medicine BioBank sample.

FINDINGS

MTAG identified genome-wide significant (GWS) single nucleotide polymorphisms (SNPs) for all four traits in EUR: 41 SNPs in 36 loci for OUD; 74 SNPs in 60 loci for CUD; 63 SNPs in 52 loci for AUD; and 183 SNPs in 144 loci for SMKinitiation. MTAG also identified GWS SNPs in AFR: 2 SNPs in 2 loci for OUD; 3 SNPs in 3 loci for AUD; and 1 SNP in 1 locus for SMKtrajectory. In the Yale-Penn sample, the MTAG-derived PRS consistently yielded more significant associations with both the corresponding substance use disorder diagnosis and multiple related phenotypes than the GWAS-derived PRS.

CONCLUSIONS

Multi-trait analysis of genome-wide association studies boosted the number of loci found for substance use traits, identifying genes not previously linked to any substance, and increased the power of polygenic risk scores. Multi-trait analysis of genome-wide association studies can be used to identify novel associations for substance use, especially those for which the samples are smaller than those for historically legal substances.

Multi-ancestry meta-analysis of tobacco use disorder prioritizes novel candidate risk genes and reveals associations with numerous health outcomes

Tobacco use disorder (TUD) is the most prevalent substance use disorder in the world. Genetic factors influence smoking behaviors, and although strides have been made using genome-wide association studies (GWAS) to identify risk variants, the majority of variants identified have been for nicotine consumption, rather than TUD. We leveraged five biobanks to perform a multi-ancestral meta-analysis of TUD (derived via electronic health records, EHR) in 898,680 individuals (739,895 European, 114,420 African American, 44,365 Latin American). We identified 88 independent risk loci; integration with functional genomic tools uncovered 461 potential risk genes, primarily expressed in the brain. TUD was genetically correlated with smoking and psychiatric traits from traditionally ascertained cohorts, externalizing behaviors in children, and hundreds of medical outcomes, including HIV infection, heart disease, and pain. This work furthers our biological understanding of TUD and establishes EHR as a source of phenotypic information for studying the genetics of TUD.

Using a polygenic score to account for genomic risk factors in a model to detect individuals with dilated ascending thoracic aortas

Background

Ascending thoracic aortic dilation is a complex trait that involves modifiable and non-modifiable risk factors and can lead to thoracic aortic aneurysm and dissection. Clinical risk factors have been shown to predict ascending thoracic aortic diameter. Polygenic scores (PGS) are increasingly used to assess clinical risk for multifactorial diseases. The degree to which a PGS can improve aortic diameter prediction is not known. In this study we tested the extent to which the addition of a PGS to clinical prediction algorithms improves the prediction of aortic diameter.

Methods

The patient cohort comprised 6,790 Penn Medicine Biobank (PMBB) participants with available echocardiography and clinical data linked to genome-wide genotype data. Linear regression models were used to integrate PGS weights derived from a large genome wide association study of thoracic aortic diameter in the UK biobank and were compared to the performance of the standard and a reweighted variation of the recently published AORTA Score.

Results

Cohort participants were 56% male, had a median age of 61 years (IQR 52-70) with a mean ascending aortic diameter of 3.4 cm (SD 0.5). Compared to the AORTA Score which explained 28.4% (95% CI 28.1% to 29.2%) of the variance in aortic diameter, AORTA Score + PGS explained 28.8%, (95% CI 28.1% to 29.6%), the reweighted AORTA score explained 30.4% (95% CI 29.6% to 31.2%), and the reweighted AORTA Score + PGS explained 31.0% (95% CI 30.2% to 31.8%). The addition of a PGS to either the AORTA Score or the reweighted AORTA Score improved model sensitivity for the identifying individuals with a thoracic aortic diameter ≥ 4 cm. The respective areas under the receiver operator characteristic curve for the AORTA Score + PGS (0.771, 95% CI 0.756 to 0.787) and reweighted AORTA Score + PGS (0.785, 95% CI 0.770 to 0.800) were greater than the standard AORTA Score (0.767, 95% CI 0.751 to 0.783) and reweighted AORTA Score (0.780 95% CI 0.765 to 0.795).

Conclusions

We demonstrated that inclusion of a PGS to the AORTA Score results in a small but clinically meaningful performance enhancement. Further investigation is necessary to determine if combining genetic and clinical risk prediction improves outcomes for thoracic aortic disease.

A genotype-first approach identifies high incidence of NF1 pathogenic variants with distinct disease associations

Loss of function variants in the NF1 gene cause neurofibromatosis type 1 (NF1), a genetic disorder characterized by complete penetrance, prevalence of 1 in 3,000, characteristic physical exam findings, and a substantially increased risk for malignancy. However, our understanding of the disorder is entirely based on patients ascertained through phenotype-first approaches. Leveraging a genotype-first approach in two large patient cohorts, we demonstrate unexpectedly high prevalence (1 in 450-750) of NF1 pathogenic variants. Half were identified in individuals lacking clinical features of NF1, with many appearing to have post-zygotic mosaicism for the identified variant. Incidentally discovered variants were not associated with classic NF1 features but were associated with an increased incidence of malignancy compared to a control population. Our findings suggest that NF1 pathogenic variants are substantially more common than previously thought, often characterized by somatic mosaicism and reduced penetrance, and are important contributors to cancer risk in the general population.

Importance of Diversity in Precision Medicine: Generalizability of Genetic Associations Across Ancestry Groups Toward Better Identification of Disease Susceptibility Variants

Genome-wide association studies (GWAS) revolutionized our understanding of common genetic variation and its impact on common human disease and traits. Developed and adopted in the mid-2000s, GWAS led to searchable genotype-phenotype catalogs and genome-wide datasets available for further data mining and analysis for the eventual development of translational applications. The GWAS revolution was swift and specific, including almost exclusively populations of European descent, to the neglect of the majority of the world's genetic diversity. In this narrative review, we recount the GWAS landscape of the early years that established a genotype-phenotype catalog that is now universally understood to be inadequate for a complete understanding of complex human genetics. We then describe approaches taken to augment the genotype-phenotype catalog, including the study populations, collaborative consortia, and study design approaches aimed to generalize and then ultimately discover genome-wide associations in non-European descent populations. The collaborations and data resources established in the efforts to diversify genomic findings undoubtedly provide the foundations of the next chapters of genetic association studies with the advent of budget-friendly whole-genome sequencing.

Biallelic BRCA Loss and Homologous Recombination Deficiency in Nonbreast/Ovarian Tumors in Germline BRCA1/2 Carriers

PURPOSE

Breast and ovarian tumors in germline BRCA1/2 carriers undergo allele-specific loss of heterozygosity, resulting in homologous recombination deficiency (HRD) and sensitivity to poly-ADP-ribose polymerase (PARP) inhibitors. This study investigated whether biallelic loss and HRD also occur in primary nonbreast/ovarian tumors that arise in germline BRCA1/2 carriers.

METHODS

A clinically ascertained cohort of BRCA1/2 carriers with a primary nonbreast/ovarian cancer was identified, including canonical (prostate and pancreatic cancers) and noncanonical (all other) tumor types. Whole-exome sequencing or clinical sequencing results (n = 45) were analyzed. A pan-cancer analysis of nonbreast/ovarian primary tumors from germline BRCA1/2 carriers from The Cancer Genome Atlas (TCGA, n = 73) was used as a validation cohort.

RESULTS

Ages of nonbreast/ovarian cancer diagnosis in germline BRCA1/2 carriers were similar to controls for the majority of cancer types. Nine of 45 (20%) primary nonbreast/ovarian tumors from germline BRCA1/2 carriers had biallelic loss of BRCA1/2 in the clinical cohort, and 23 of 73 (32%) in the TCGA cohort. In the combined cohort, 35% and 27% of primary canonical and noncanonical BRCA tumor types, respectively, had biallelic loss. High HRD scores (HRDex > 42) were detected in 81% of tumors with biallelic BRCA loss compared with 22% (P < .001) of tumors without biallelic BRCA loss. No differences in genomic profile, including mutational signatures, mutation spectrum, tumor mutational burden, or microsatellite instability, were found in primary nonbreast/ovarian tumors with or without biallelic BRCA1/2 loss.

CONCLUSION

A proportion of noncanonical primary tumors have biallelic loss and evidence of HRD. Our data suggest that assessment of biallelic loss and HRD could supplement identification of germline BRCA1/2 mutations in selection of patients for platinum or PARP inhibitor therapy.

Diversity and Scale: Genetic Architecture of 2,068 Traits in the VA Million Veteran Program

Genome-wide association studies (GWAS) have underrepresented individuals from non-European populations, impeding progress in characterizing the genetic architecture and consequences of health and disease traits. To address this, we present a population-stratified phenome-wide GWAS followed by a multi-population meta-analysis for 2,068 traits derived from electronic health records of 635,969 participants in the Million Veteran Program (MVP), a longitudinal cohort study of diverse U.S. Veterans genetically similar to the respective African (121,177), Admixed American (59,048), East Asian (6,702), and European (449,042) superpopulations defined by the 1000 Genomes Project. We identified 38,270 independent variants associating with one or more traits at experiment-wide P < 4.6 × 10 - 11 significance; fine-mapping 6,318 signals identified from 613 traits to single-variant resolution. Among these, a third (2,069) of the associations were found only among participants genetically similar to non-European reference populations, demonstrating the importance of expanding diversity in genetic studies. Our work provides a comprehensive atlas of phenome-wide genetic associations for future studies dissecting the architecture of complex traits in diverse populations.

The genetic architecture of pain intensity in a sample of 598,339 U.S. veterans

Chronic pain is a common problem, with more than one-fifth of adult Americans reporting pain daily or on most days. It adversely affects quality of life and imposes substantial personal and economic costs. Efforts to treat chronic pain using opioids played a central role in precipitating the opioid crisis. Despite an estimated heritability of 25-50%, the genetic architecture of chronic pain is not well characterized, in part because studies have largely been limited to samples of European ancestry. To help address this knowledge gap, we conducted a cross-ancestry meta-analysis of pain intensity in 598,339 participants in the Million Veteran Program, which identified 125 independent genetic loci, 82 of which are novel. Pain intensity was genetically correlated with other pain phenotypes, level of substance use and substance use disorders, other psychiatric traits, education level, and cognitive traits. Integration of the GWAS findings with functional genomics data shows enrichment for putatively causal genes (n = 142) and proteins (n = 14) expressed in brain tissues, specifically in GABAergic neurons. Drug repurposing analysis identified anticonvulsants, beta-blockers, and calcium-channel blockers, among other drug groups, as having potential analgesic effects. Our results provide insights into key molecular contributors to the experience of pain and highlight attractive drug targets.

Genetic Predictors of Blood Pressure Traits are Associated with Preeclampsia

Background

Preeclampsia, a pregnancy complication characterized by hypertension after 20 gestational weeks, is a major cause of maternal and neonatal morbidity and mortality. The mechanisms leading to preeclampsia are unclear; however, there is evidence that preeclampsia is highly heritable. We evaluated the association of polygenic risk scores (PRS) for blood pressure traits and preeclampsia to assess whether there is shared genetic architecture.

Methods

Participants were obtained from Vanderbilt University's BioVU, the Electronic Medical Records and Genomics network, and the Penn Medicine Biobank. Non-Hispanic Black and White females of reproductive age with indications of pregnancy and genotype information were included. Preeclampsia was defined by ICD codes. Summary statistics for diastolic blood pressure (DBP), systolic blood pressure (SBP), and pulse pressure (PP) PRS were obtained from Giri et al 2019. Associations between preeclampsia and each PRS were evaluated separately by race and study population before evidence was meta-analyzed. Prediction models were developed and evaluated using 10-fold cross validation.

Results

In the 3,504 Black and 5,009 White individuals included, the rate of preeclampsia was 15.49%. The DBP and SBP PRSs were associated with preeclampsia in Whites but not Blacks. The PP PRS was significantly associated with preeclampsia in Blacks and Whites. In trans-ancestry meta-analysis, all PRSs were associated with preeclampsia (OR DBP =1.10, 95% CI=1.02-1.17, p =7.68×10 -3 ; OR SBP =1.16, 95% CI=1.09-1.23, p =2.23×10 -6 ; OR PP =1.14, 95% CI=1.07-1.27, p =9.86×10 -5 ). However, addition of PRSs to clinical prediction models did not improve predictive performance.

Conclusions

Genetic factors contributing to blood pressure regulation in the general population also predispose to preeclampsia.

Multi-ancestry study of the genetics of problematic alcohol use in >1 million individuals

Problematic alcohol use (PAU) is a leading cause of death and disability worldwide. To improve our understanding of the genetics of PAU, we conducted a large cross-ancestry meta-analysis of PAU in 1,079,947 individuals. We observed a high degree of cross-ancestral similarity in the genetic architecture of PAU and identified 110 independent risk variants in within- and cross-ancestry analyses. Cross-ancestry fine-mapping improved the identification of likely causal variants. Prioritizing genes through gene expression and/or chromatin interaction in brain tissues identified multiple genes associated with PAU. We identified existing medications for potential pharmacological studies by drug repurposing analysis. Cross-ancestry polygenic risk scores (PRS) showed better performance in independent sample than single-ancestry PRS. Genetic correlations between PAU and other traits were observed in multiple ancestries, with other substance use traits having the highest correlations. The analysis of diverse ancestries contributed significantly to the findings, and fills an important gap in the literature.