Next-generation phenotyping: introducing phecodeX for enhanced discovery research in medical phenomics

MOTIVATION

Phecodes are widely used and easily adapted phenotypes based on International Classification of Diseases codes. The current version of phecodes (v1.2) was designed primarily to study common/complex diseases diagnosed in adults; however, there are numerous limitations in the codes and their structure.

RESULTS

Here, we present phecodeX, an expanded version of phecodes with a revised structure and 1,761 new codes. PhecodeX adds granularity to phenotypes in key disease domains that are under-represented in the current phecode structure-including infectious disease, pregnancy, congenital anomalies, and neonatology-and is a more robust representation of the medical phenome for global use in discovery research.

AVAILABILITY AND IMPLEMENTATION

phecodeX is available at https://github.com/PheWAS/phecodeX.

Correction: Hypertension genetic risk score is associated with burden of coronary heart disease among patients referred for coronary angiography

[This corrects the article DOI: 10.1371/journal.pone.0208645.].

Determinants of mosaic chromosomal alteration fitness

Clonal hematopoiesis (CH) is characterized by the acquisition of a somatic mutation in a hematopoietic stem cell that results in a clonal expansion. These driver mutations can be single nucleotide variants in cancer driver genes or larger structural rearrangements called mosaic chromosomal alterations (mCAs). The factors that influence the variations in mCA fitness and ultimately result in different clonal expansion rates are not well-understood. We used the Passenger-Approximated Clonal Expansion Rate (PACER) method to estimate clonal expansion rate for 6,381 individuals in the NHLBI TOPMed cohort with gain, loss, and copy-neutral loss of heterozygosity mCAs. Our estimates of mCA fitness were correlated (R 2 = 0.49) with an alternative approach that estimated fitness of mCAs in the UK Biobank using a theoretical probability distribution. Individuals with lymphoid-associated mCAs had a significantly higher white blood cell count and faster clonal expansion rate. In a cross-sectional analysis, genome-wide association study of estimates of mCA expansion rate identified TCL1A , NRIP1 , and TERT locus variants as modulators of mCA clonal expansion rate.

Proteomic and genetic analyses of influenza A viruses identify pan-viral host targets

Influenza A Virus (IAV) is a recurring respiratory virus with limited availability of antiviral therapies. Understanding host proteins essential for IAV infection can identify targets for alternative host-directed therapies (HDTs). Using affinity purification-mass spectrometry and global phosphoproteomic and protein abundance analyses using three IAV strains (pH1N1, H3N2, H5N1) in three human cell types (A549, NHBE, THP-1), we map 332 IAV-human protein-protein interactions and identify 13 IAV-modulated kinases. Whole exome sequencing of patients who experienced severe influenza reveals several genes, including scaffold protein AHNAK, with predicted loss-of-function variants that are also identified in our proteomic analyses. Of our identified host factors, 54 significantly alter IAV infection upon siRNA knockdown, and two factors, AHNAK and coatomer subunit COPB1, are also essential for productive infection by SARS-CoV-2. Finally, 16 compounds targeting our identified host factors suppress IAV replication, with two targeting CDK2 and FLT3 showing pan-antiviral activity across influenza and coronavirus families. This study provides a comprehensive network model of IAV infection in human cells, identifying functional host targets for pan-viral HDT.

Detection of distant relatedness in biobanks for identification of undiagnosed carriers of a Mendelian disease variant: application to Long QT Syndrome

Rare genetic diseases are typically studied in referral populations, resulting in underdiagnosis and biased assessment of penetrance and phenotype. To address this, we developed a generalizable method of genotype inference based on distant relatedness and deployed this to identify undiagnosed Type 5 Long QT Syndrome (LQT5) rare variant carriers in a non-referral population. We identified 9 LQT5 families referred to a single specialty clinic, each carrying p.Asp76Asn, the most common LQT5 variant. We uncovered recent common ancestry and a single shared haplotype among probands. Application to a non-referral population of 69,819 BioVU biobank subjects identified 22 additional subjects sharing this haplotype, subsequently confirmed to carry p.Asp76Asn. Referral and non-referral carriers had prolonged QTc compared to controls, and, among carriers, QTc polygenic score additively associated with QTc prolongation. Thus, our novel analysis of shared chromosomal segments identified undiagnosed cases of genetic disease and refined the understanding of LQT5 penetrance and phenotype.

ParSE-seq: A Calibrated Multiplexed Assay to Facilitate the Clinical Classification of Putative Splice-altering Variants

Background

Interpreting the clinical significance of putative splice-altering variants outside 2-base pair canonical splice sites remains difficult without functional studies.

Methods

We developed Parallel Splice Effect Sequencing (ParSE-seq), a multiplexed minigene-based assay, to test variant effects on RNA splicing quantified by high-throughput sequencing. We studied variants in SCN5A, an arrhythmia-associated gene which encodes the major cardiac voltage-gated sodium channel. We used the computational tool SpliceAI to prioritize exonic and intronic candidate splice variants, and ClinVar to select benign and pathogenic control variants. We generated a pool of 284 barcoded minigene plasmids, transfected them into Human Embryonic Kidney (HEK293) cells and induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs), sequenced the resulting pools of splicing products, and calibrated the assay to the American College of Medical Genetics and Genomics scheme. Variants were interpreted using the calibrated functional data, and experimental data were compared to SpliceAI predictions. We further studied some splice-altering missense variants by cDNA-based automated patch clamping (APC) in HEK cells and assessed splicing and sodium channel function in CRISPR-edited iPSC-CMs.

Results

ParSE-seq revealed the splicing effect of 224 SCN5A variants in iPSC-CMs and 244 variants in HEK293 cells. The scores between the cell types were highly correlated (R2=0.84). In iPSCs, the assay had concordant scores for 21/22 benign/likely benign and 24/25 pathogenic/likely pathogenic control variants from ClinVar. 43/112 exonic variants and 35/70 intronic variants with determinate scores disrupted splicing. 11 of 42 variants of uncertain significance were reclassified, and 29 of 34 variants with conflicting interpretations were reclassified using the functional data. SpliceAI computational predictions correlated well with experimental data (AUC = 0.96). We identified 20 unique SCN5A missense variants that disrupted splicing, and 2 clinically observed splice-altering missense variants of uncertain significance had normal function when tested with the cDNA-based APC assay. A splice-altering intronic variant detected by ParSE-seq, c.1891-5C>G, also disrupted splicing and sodium current when introduced into iPSC-CMs at the endogenous locus by CRISPR editing.

Conclusions

ParSE-seq is a calibrated, multiplexed, high-throughput assay to facilitate the classification of candidate splice-altering variants.

Demonstrating paths for unlocking the value of cloud genomics through cross cohort analysis

Recently, large scale genomic projects such as All of Us and the UK Biobank have introduced a new research paradigm where data are stored centrally in cloud-based Trusted Research Environments (TREs). To characterize the advantages and drawbacks of different TRE attributes in facilitating cross-cohort analysis, we conduct a Genome-Wide Association Study of standard lipid measures using two approaches: meta-analysis and pooled analysis. Comparison of full summary data from both approaches with an external study shows strong correlation of known loci with lipid levels (R2 ~ 83-97%). Importantly, 90 variants meet the significance threshold only in the meta-analysis and 64 variants are significant only in pooled analysis, with approximately 20% of variants in each of those groups being most prevalent in non-European, non-Asian ancestry individuals. These findings have important implications, as technical and policy choices lead to cross-cohort analyses generating similar, but not identical results, particularly for non-European ancestral populations.

Laboratory and demographic predictors of functional assay positive status in suspected heparin-induced thrombocytopenia: A multicenter retrospective cohort study

Heparin-induced thrombocytopenia (HIT) is an antibody-mediated immune response against platelet factor 4 (PF4) bound to heparin anticoagulants. A priori identification of patients at-risk for HIT remains elusive and a number of risk factors have been identified, but these associations and their effect sizes have limited validation in large cohorts of suspected HIT patients. The aim of this study was to investigate existing anti-PF4/heparin antibody thresholds and model the relationship of demographic variables and anti-PF4/heparin antibody levels with functional assay positivity across multiple institutions in the absence of detailed clinical data. In a large collection of suspected HIT patients (n = 8904), we tested for associations between laboratory and demographic variables and functional assay positive status as well as anti-PF4/heparin antibody levels. We also tested for correlation between IgG-specific and polyspecific (IgG/IgA/IgM) anti-PF4/heparin antibody values and their ability to predict functional assay positive status using area under the receiver operating characteristic (AUROC). Logistic regression identified increasing anti-PF4/heparin antibody OD levels (OR = 51.84 [37.27-74.34], p < 2.0 × 10-16) and female sex (OR = 1.47 [1.19-1.82], p = 3.5 × 10-4) as risk factors for positive functional assay in the largest cohort with consistent effect sizes in two other cohorts. In a subset of 1175 patients, polyspecific and IgG-specific anti-PF4/heparin antibody values were heterogeneous (mean coefficient of variation = 31.9 %), but strongly correlated (rho = 0.878; p < 2 × 10-16) with similar prediction of functional assay positivity (polyspecific AUROC = 0.976 and IgG-specific AUROC = 0.980). Thus, we recapitulate previously identified risk factors of functional assay positivity, providing precise effect sizes in a large observational population of suspected HIT patients. Our data reinforce the necessity of functional assay confirmation and suggest that, despite heterogeneity, polyspecific and IgG-specific anti-PF4/heparin antibody assays predict functional assay positive status similarly, even in the absence of 4Ts scores and detailed clinical data.

WHOLE GENOME SEQUENCING ANALYSIS OF BODY MASS INDEX IDENTIFIES NOVEL AFRICAN ANCESTRY-SPECIFIC RISK ALLELE

Obesity is a major public health crisis associated with high mortality rates. Previous genome-wide association studies (GWAS) investigating body mass index (BMI) have largely relied on imputed data from European individuals. This study leveraged whole-genome sequencing (WGS) data from 88,873 participants from the Trans-Omics for Precision Medicine (TOPMed) Program, of which 51% were of non-European population groups. We discovered 18 BMI-associated signals (P < 5 × 10-9). Notably, we identified and replicated a novel low frequency single nucleotide polymorphism (SNP) in MTMR3 that was common in individuals of African descent. Using a diverse study population, we further identified two novel secondary signals in known BMI loci and pinpointed two likely causal variants in the POC5 and DMD loci. Our work demonstrates the benefits of combining WGS and diverse cohorts in expanding current catalog of variants and genes confer risk for obesity, bringing us one step closer to personalized medicine.

Clinical consequences of a polygenic predisposition to benign lower white blood cell counts: Consequences of benign WBC count genetics

Polygenic variation unrelated to disease contributes to interindividual variation in baseline white blood cell (WBC) counts, but its clinical significance is undefined. We investigated the clinical consequences of a genetic predisposition toward lower WBC counts among 89,559 biobank participants from tertiary care centers using a polygenic score for WBC count (PGSWBC) comprising single nucleotide polymorphisms not associated with disease. A predisposition to lower WBC counts was associated with a decreased risk of identifying pathology on a bone marrow biopsy performed for a low WBC count (odds-ratio=0.55 per standard deviation increase in PGSWBC [95%CI, 0.30 - 0.94], p=0.04), an increased risk of leukopenia (a low WBC count) when treated with a chemotherapeutic (n=1,724, hazard ratio [HR]=0.78 [0.69 - 0.88], p=4.0×10-5) or immunosuppressant (n=354, HR=0.61 [0.38 - 0.99], p=0.04). A predisposition to benign lower WBC counts was associated with an increased risk of discontinuing azathioprine treatment (n=1,466, HR=0.62 [0.44 - 0.87], p=0.006). Collectively, these findings suggest that a WBC count polygenic score identifies individuals who are susceptible to escalations or alterations in clinical care that may be harmful or of little benefit.

The All of Us Data and Research Center: Creating a Secure, Scalable, and Sustainable Ecosystem for Biomedical Research

The All of Us Research Program's Data and Research Center (DRC) was established to help acquire, curate, and provide access to one of the world's largest and most diverse datasets for precision medicine research. Already, over 500,000 participants are enrolled in All of Us, 80% of whom are underrepresented in biomedical research, and data are being analyzed by a community of over 2,300 researchers. The DRC created this thriving data ecosystem by collaborating with engaged participants, innovative program partners, and empowered researchers. In this review, we first describe how the DRC is organized to meet the needs of this broad group of stakeholders. We then outline guiding principles, common challenges, and innovative approaches used to build the All of Us data ecosystem. Finally, we share lessons learned to help others navigate important decisions and trade-offs in building a modern biomedical data platform.

Multicenter clinical and functional evidence reclassifies a recurrent noncanonical filamin C splice-altering variant

BACKGROUND

Truncating variants in filamin C (FLNC) can cause arrhythmogenic cardiomyopathy (ACM) through haploinsufficiency. Noncanonical splice-altering variants may contribute to this phenotype.

OBJECTIVE

The purpose of this study was to investigate the clinical and functional consequences of a recurrent FLNC intronic variant of uncertain significance (VUS), c.970-4A>G.

METHODS

Clinical data in 9 variant heterozygotes from 4 kindreds were obtained from 5 tertiary health care centers. We used in silico predictors and functional studies with peripheral blood and patient-specific induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs). Isolated RNA was studied by reverse transcription polymerase chain reaction. iPSC-CMs were further characterized at baseline and after nonsense-mediated decay (NMD) inhibition, using quantitative polymerase chain reaction (qPCR), RNA-sequencing, and cellular electrophysiology. American College of Medical Genetics and Genomics (ACMG) criteria were used to adjudicate variant pathogenicity.

RESULTS

Variant heterozygotes displayed a spectrum of disease phenotypes, spanning from mild ventricular dysfunction with palpitations to severe ventricular arrhythmias requiring device shocks or progressive cardiomyopathy requiring heart transplantation. Consistent with in silico predictors, the c.970-4A>G FLNC variant activated a cryptic splice acceptor site, introducing a 3-bp insertion containing a premature termination codon. NMD inhibition upregulated aberrantly spliced transcripts by qPCR and RNA-sequencing. Patch clamp studies revealed irregular spontaneous action potentials, increased action potential duration, and increased sodium late current in proband-derived iPSC-CMs. These findings fulfilled multiple ACMG criteria for pathogenicity.

CONCLUSION

Clinical, in silico, and functional evidence support the prediction that the intronic c.970-4A>G VUS disrupts splicing and drives ACM, enabling reclassification from VUS to pathogenic.

PheMIME: An Interactive Web App and Knowledge Base for Phenome-Wide, Multi-Institutional Multimorbidity Analysis

Motivation

Multimorbidity, characterized by the simultaneous occurrence of multiple diseases in an individual, is an increasing global health concern, posing substantial challenges to healthcare systems. Comprehensive understanding of disease-disease interactions and intrinsic mechanisms behind multimorbidity can offer opportunities for innovative prevention strategies, targeted interventions, and personalized treatments. Yet, there exist limited tools and datasets that characterize multimorbidity patterns across different populations. To bridge this gap, we used large-scale electronic health record (EHR) systems to develop the Phenome-wide Multi-Institutional Multimorbidity Explorer (PheMIME), which facilitates research in exploring and comparing multimorbidity patterns among multiple institutions, potentially leading to the discovery of novel and robust disease associations and patterns that are interoperable across different systems and organizations.

Results

PheMIME integrates summary statistics from phenome-wide analyses of disease multimorbidities. These are currently derived from three major institutions: Vanderbilt University Medical Center, Mass General Brigham, and the UK Biobank. PheMIME offers interactive exploration of multimorbidity through multi-faceted visualization. Incorporating an enhanced version of associationSubgraphs, PheMIME enables dynamic analysis and inference of disease clusters, promoting the discovery of multimorbidity patterns. Once a disease of interest is selected, the tool generates interactive visualizations and tables that users can delve into multimorbidities or multimorbidity networks within a single system or compare across multiple systems. The utility of PheMIME is demonstrated through a case study on schizophrenia.

Availability and implementation

The PheMIME knowledge base and web application are accessible at https://prod.tbilab.org/PheMIME/. A comprehensive tutorial, including a use-case example, is available at https://prod.tbilab.org/PheMIME_supplementary_materials/. Furthermore, the source code for PheMIME can be freely downloaded from https://github.com/tbilab/PheMIME.

Data availability statement

The data underlying this article are available in the article and in its online web application or supplementary material.

Leveraging Generative AI to Prioritize Drug Repurposing Candidates: Validating Identified Candidates for Alzheimer's Disease in Real-World Clinical Datasets

Drug repurposing represents an attractive alternative to the costly and time-consuming process of new drug development, particularly for serious, widespread conditions with limited effective treatments, such as Alzheimer's disease (AD). Emerging generative artificial intelligence (GAI) technologies like ChatGPT offer the promise of expediting the review and summary of scientific knowledge. To examine the feasibility of using GAI for identifying drug repurposing candidates, we iteratively tasked ChatGPT with proposing the twenty most promising drugs for repurposing in AD, and tested the top ten for risk of incident AD in exposed and unexposed individuals over age 65 in two large clinical datasets: 1) Vanderbilt University Medical Center and 2) the All of Us Research Program. Among the candidates suggested by ChatGPT, metformin, simvastatin, and losartan were associated with lower AD risk in meta-analysis. These findings suggest GAI technologies can assimilate scientific insights from an extensive Internet-based search space, helping to prioritize drug repurposing candidates and facilitate the treatment of diseases.

Pulmonary Vein Myocardial Sleeve Length and its Association With Sex and 4q25/PITX2 Genotype

BACKGROUND

Experimental evidence suggests genetic variation in 4q25/PITX2 modulates pulmonary vein (PV) myocardial sleeve length. Although PV sleeves are the main target of atrial fibrillation (AF) ablation, little is known about the association between different PV sleeve characteristics with ablation outcomes.

OBJECTIVES

This study sought to evaluate the association between clinical and genetic (4q25) risk factors with PV sleeve length in humans, and to evaluate the association between PV sleeve length and recurrence after AF ablation.

METHODS

In a prospective, observational study of patients undergoing de novo AF ablation, PV sleeve length was measured using electroanatomic voltage mapping before ablation. The sentinel 4q25 AF susceptibility single nucleotide polymorphism, rs2200733, was genotyped. The primary analysis tested the association between clinical and genetic (4q25) risk factors with PV sleeve length using a multivariable linear regression model. Covariates included age, sex, body mass index, height, and persistent AF. The association between PV sleeve length and atrial arrhythmia recurrence (>30 seconds) was tested using a multivariable Cox proportional hazards model.

RESULTS

Between 2014 and 2019, 197 participants were enrolled (median age 63 years [IQR: 55 to 70 years], 133 male [67.5%]). In multivariable modeling, men were found to have PV sleeves 2.94 mm longer than women (95% CI: 0.99-4.90 mm; P < 0.001). Sixty participants (30.5%) had one 4q25 risk allele and 6 (3.1%) had 2 alleles. There was no association between 4q25 genotype and PV sleeve length. Forty-six participants (23.4%) experienced arrhythmia recurrence within 3 to 12 months, but there was no association between recurrence and PV sleeve length.

CONCLUSIONS

Common genetic variation at 4q25 was not associated with PV sleeve length and PV sleeve length was not associated with ablation outcomes. Men did have longer PV sleeves than women, but more research is needed to define the potential clinical significance of this observation.

Cardiac Sarcoidosis and a Likely Pathogenic TTN Variant in a Patient Presenting With Ventricular Tachycardia

Rare variants in TTN are the most common monogenic cause of early-onset atrial fibrillation and dilated cardiomyopathy. Whereas cardiac sarcoidosis is very underdiagnosed, a common presentation can be ventricular arrhythmias. This report presents a patient with a likely pathogenic TTN variant and cardiac sarcoidosis. (Level of Difficulty: Intermediate.).

High-throughput functional mapping of variants in an arrhythmia gene, KCNE1, reveals novel biology

Background

KCNE1 encodes a 129-residue cardiac potassium channel (IKs) subunit. KCNE1 variants are associated with long QT syndrome and atrial fibrillation. However, most variants have insufficient evidence of clinical consequences and thus limited clinical utility.

Results

Here, we demonstrate the power of variant effect mapping, which couples saturation mutagenesis with high-throughput sequencing, to ascertain the function of thousands of protein coding KCNE1 variants. We comprehensively assayed KCNE1 variant cell surface expression (2,554/2,709 possible single amino acid variants) and function (2,539 variants). We identified 470 loss-of-surface expression and 588 loss-of-function variants. Out of the 588 loss-of-function variants, only 155 had low cell surface expression. The latter half of the protein is dispensable for protein trafficking but essential for channel function. 22 of the 30 KCNE1 residues (73%) highly intolerant of variation were in predicted close contact with binding partners KCNQ1 or calmodulin. Our data were highly concordant with gold standard electrophysiological data (ρ = -0.65), population and patient cohorts (32/38 concordant variants), and computational metrics (ρ = -0.55). Our data provide moderate-strength evidence for the ACMG/AMP functional criteria for benign and pathogenic variants.

Conclusions

Comprehensive variant effect maps of KCNE1 can both provide insight into IKs channel biology and help reclassify variants of uncertain significance.

Detection of distant familial relatedness in biobanks for identification of undiagnosed carriers of a Mendelian disease variant: application to Long QT syndrome

Importance

The diagnosis and study of rare genetic disease is often limited to referral populations, leading to underdiagnosis and a biased assessment of penetrance and phenotype.

Objective

To develop a generalizable method of genotype inference based on distant relatedness and to deploy this to identify undiagnosed Type 5 Long QT Syndrome (LQT5) rare variant carriers in a non-referral population.

Participants

We identified 9 LQT5 probands and 3 first-degree relatives referred to a single Genetic Arrhythmia clinic, each carrying D76N (p.Asp76Asn), the most common variant implicated in LQT5. The non-referral population consisted of 69,879 ancestry-matched subjects in BioVU, a large biobank that links electronic health records to dense array data. Participants were enrolled from 2007-2022. Data analysis was performed in 2022.

Exposures

We developed and applied a novel approach to genotype inference (Distant Relatedness for Identification and Variant Evaluation, or DRIVE) to identify shared, identical-by-descent (IBD) large chromosomal segments in array data.

Main Outcomes and Measures

We sought to establish genetic relatedness among the probands and to use genomic segments underlying D76N to identify other potential carriers in BioVU. We then further studied the role of D76N in LQT5 pathogenesis.

Results

Genetic reconstruction of pedigrees and distant relatedness detection among clinic probands using DRIVE revealed shared recent common ancestry and identified a single long shared haplotype. Interrogation of the non-referral population in BioVU identified a further 23 subjects sharing this haplotype, and sequencing confirmed D76N carrier status in 22, all previously undiagnosed with LQT5. The QTc was prolonged in D76N carriers compared to BioVU controls, with 40% penetrance of QTc ≥ 480 msec. Among D76N carriers, a QTc polygenic score was additively associated with QTc prolongation.

Conclusions and Relevance

Detection of IBD shared chromosomal segments around D76N enabled identification of distantly related and previously undiagnosed rare-variant carriers, demonstrated the contribution of polygenic risk to monogenic disease penetrance, and further established LQT5 as a primary arrhythmia disorder. Analysis of shared chromosomal regions spanning disease-causing mutations can identify undiagnosed cases of genetic diseases.

Association of Longitudinal Activity Measures and Diabetes Risk: An Analysis From the National Institutes of Health All of Us Research Program

CONTEXT

Prior studies of the relationship between physical activity and incident type 2 diabetes mellitus (T2DM) relied primarily on questionnaires at a single time point.

OBJECTIVE

We sought to investigate the relationship between physical activity and incident T2DM with an innovative approach using data from commercial wearable devices linked to electronic health records in a real-world population.

METHODS

Using All of Us participants' accelerometer data from their personal Fitbit devices, we used a time-varying Cox proportional hazards models with repeated measures of physical activity for the outcome of incident T2DM. We evaluated for effect modification with age, sex, body mass index (BMI), and sedentary time using multiplicative interaction terms.

RESULTS

From 5677 participants in the All of Us Research Program (median age 51 years; 74% female; 89% White), there were 97 (2%) cases of incident T2DM over a median follow-up period of 3.8 years between 2010 to 2021. In models adjusted for age, sex, and race, the hazard of incident diabetes was reduced by 44% (95% CI, 15%-63%; P = 0.01) when comparing those with an average daily step count of 10 700 to those with 6000. Similar benefits were seen comparing groups based on average duration of various intensities of activity (eg, lightly active, fairly active, very active). There was no evidence for effect modification by age, sex, BMI, or sedentary time.

CONCLUSION

Greater time in any type of physical activity intensity was associated with lower risk of T2DM irrespective of age, sex, BMI, or sedentary time.

Returning integrated genomic risk and clinical recommendations: The eMERGE study

PURPOSE

Assessing the risk of common, complex diseases requires consideration of clinical risk factors as well as monogenic and polygenic risks, which in turn may be reflected in family history. Returning risks to individuals and providers may influence preventive care or use of prophylactic therapies for those individuals at high genetic risk.

METHODS

To enable integrated genetic risk assessment, the eMERGE (electronic MEdical Records and GEnomics) network is enrolling 25,000 diverse individuals in a prospective cohort study across 10 sites. The network developed methods to return cross-ancestry polygenic risk scores, monogenic risks, family history, and clinical risk assessments via a genome-informed risk assessment (GIRA) report and will assess uptake of care recommendations after return of results.

RESULTS

GIRAs include summary care recommendations for 11 conditions, education pages, and clinical laboratory reports. The return of high-risk GIRA to individuals and providers includes guidelines for care and lifestyle recommendations. Assembling the GIRA required infrastructure and workflows for ingesting and presenting content from multiple sources. Recruitment began in February 2022.

CONCLUSION

Return of a novel report for communicating monogenic, polygenic, and family history-based risk factors will inform the benefits of integrated genetic risk assessment for routine health care.

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.

Aberrant activation of TCL1A promotes stem cell expansion in clonal haematopoiesis

Mutations in a diverse set of driver genes increase the fitness of haematopoietic stem cells (HSCs), leading to clonal haematopoiesis1. These lesions are precursors for blood cancers2-6, but the basis of their fitness advantage remains largely unknown, partly owing to a paucity of large cohorts in which the clonal expansion rate has been assessed by longitudinal sampling. Here, to circumvent this limitation, we developed a method to infer the expansion rate from data from a single time point. We applied this method to 5,071 people with clonal haematopoiesis. A genome-wide association study revealed that a common inherited polymorphism in the TCL1A promoter was associated with a slower expansion rate in clonal haematopoiesis overall, but the effect varied by driver gene. Those carrying this protective allele exhibited markedly reduced growth rates or prevalence of clones with driver mutations in TET2, ASXL1, SF3B1 and SRSF2, but this effect was not seen in clones with driver mutations in DNMT3A. TCL1A was not expressed in normal or DNMT3A-mutated HSCs, but the introduction of mutations in TET2 or ASXL1 led to the expression of TCL1A protein and the expansion of HSCs in vitro. The protective allele restricted TCL1A expression and expansion of mutant HSCs, as did experimental knockdown of TCL1A expression. Forced expression of TCL1A promoted the expansion of human HSCs in vitro and mouse HSCs in vivo. Our results indicate that the fitness advantage of several commonly mutated driver genes in clonal haematopoiesis may be mediated by TCL1A activation.

The impact of damaging epilepsy and cardiac genetic variant burden in sudden death in the young

Background

Sudden unexpected death in children is a tragic event. Understanding the genetics of sudden death in the young (SDY) enables family counseling and cascade screening. The objective of this study was to characterize genetic variation in an SDY cohort using whole genome sequencing.

Methods

The SDY Case Registry is a National Institutes of Health/Centers for Disease Control surveillance effort to discern the prevalence, causes, and risk factors for SDY. The SDY Case Registry prospectively collected clinical data and DNA biospecimens from SDY cases <20 years of age. SDY cases were collected from medical examiner and coroner offices spanning 13 US jurisdictions from 2015-2019. The cohort included 211 children (mean age 1 year; range 0-20 years), determined to have died suddenly and unexpectedly and in whom DNA biospecimens and next-of-kin consent were ascertained. A control cohort consisted of 211 randomly sampled, sex-and ancestry-matched individuals from the 1000 Genomes Project. Genetic variation was evaluated in epilepsy, cardiomyopathy and arrhythmia genes in the SDY and control cohorts. American College of Medical Genetics/Genomics guidelines were used to classify variants as pathogenic or likely pathogenic. Additionally, genetic variation predicted to be damaging was identified using a Bayesian-based artificial intelligence (AI) tool.

Results

The SDY cohort was 42% European, 30% African, 17% Hispanic, and 11% with mixed ancestries, and 39% female. Six percent of the cohort was found to harbor a pathogenic or likely pathogenic genetic variant in an epilepsy, cardiomyopathy or arrhythmia gene. The genomes of SDY cases, but not controls, were enriched for rare, damaging variants in epilepsy, cardiomyopathy and arrhythmia-related genes. A greater number of rare epilepsy genetic variants correlated with younger age at death.

Conclusions

While damaging cardiomyopathy and arrhythmia genes are recognized contributors to SDY, we also observed an enrichment in epilepsy-related genes in the SDY cohort, and a correlation between rare epilepsy variation and younger age at death. These findings emphasize the importance of considering epilepsy genes when evaluating SDY.