Citations

BACKGROUND: Anthracyclines are a widely used and effective class of chemotherapy. However, a major limitation for their use is cardiotoxicity, manifesting as systolic dysfunction, congestive heart failure, and arrhythmias. The objective of this study was to identify genes and genetic risk variants for doxorubicin-induced arrhythmia. METHODS: We performed genome-wide association studies and meta-analysis across 2 independent nonoverlapping, genetically homogenous data sets from the Molecular Epidemiology Resource (cases N=77, controls N=1184) and Tapestry (cases N=13, controls N=172). Doxorubicin-related arrhythmia was the primary cardiac outcome, followed by investigation of a broader phenotype of cardiac events. Effects of doxorubicin on top-associated genes were assessed by quantitative polymerase chain reaction on RNA extracted from human cardiomyocytes following treatment with doxorubicin and by echocardiography in a zebrafish model of doxorubicin-induced cardiomyopathy. RESULTS: In the meta-analysis, rs12151014 at the KCTD15 locus was associated with doxorubicin-related arrhythmia at the genome-wide significance (minor allele frequency in Molecular Epidemiology Resource cases 0.195 versus 0.069 in controls; Tapestry cases 0.192 versus 0.079 in controls) P =9.01×10 −9 ; odds ratio, 3.61. MMP 16 rs1094199 was associated with the broader phenotype of arrhythmia and heart failure (minor allele frequency in Molecular Epidemiology Resource cases 0.079 versus 0.019 in controls; Tapestry cases 0.038 versus 0.015) P =1.98×10 − 8 ; odds ratio=6.15. In human cardiomyocytes, MMP16 and KCTD15 gene expression was significantly downregulated by doxorubicin, and in a zebrafish MMP16 knockout was protective of doxorubicin-induced decline in left ventricular ejection fraction. CONCLUSIONS: Genetic variants at KCTD15 and MMP16 are associated with doxorubicin-related arrhythmia and cardiac events. MMP16 is a potential cardioprotective therapeutic target.

Abstract Rationale Idiopathic Pulmonary Fibrosis (IPF) is characterized by chronic progressive pulmonary fibrosis and high mortality. Genetic markers, summarized into a polygenic risk score (PRS), associate with IPF in well-phenotyped research cohorts. Objectives To evaluate the performance of the PRS using real-world data from routinely captured electronic healthcare records. Methods We conducted an observational study evaluating the association of a PRS for IPF with electronic healthcare record diagnosis of IPF as well as lung transplant-free survival in four independent cohorts; the Mass General Brigham Biobank (MGBB), Mayo Clinic Biobank (MCBB), Mayo Clinic Tapestry Cohort (Tapestry), and U.K. Biobank (UKBB). We used multivariable logistic regression and multivariable Cox proportional hazards models adjusting for age, gender, and principal components of ancestry. The cohorts then underwent fixed and random effects meta-analysis. Measurements and Main Results Of 37,709; 44,195; 43,202; and 447,422 participants from MGBB, MCBB, Tapestry and UKBB respectively, 1,015 (2.7%) 2,879 (6.5%), 1,310 (3.0%), and 2,742 (0.6%) participants had an IPF diagnosis. Meta-analysis demonstrated a high-risk PRS associated with IPF diagnosis, OR 2.88 (95%CI 2.41-3.44) compared to all other individuals. A high-risk PRS also associated with the composite endpoint of mortality or lung transplant among those with an IPF diagnosis, HR 1.23(95%CI 1.11-1.35) compared to all other individuals. Conclusions A PRS can identify those at risk for an IPF diagnosis and mortality in biobank-scale data, which may have implications for clinical decisions. Further work is necessary to evaluate the utility of adding genetics in clinical settings.

Abstract Rationale Asthma contributes substantially to healthcare burden and morbidity through frequent exacerbations. Emerging evidence suggests that genetic factors play a crucial role in modulating susceptibility to asthma exacerbations. To better understand the genetic architecture of exacerbations, we conducted a genome-wide association study (GWAS) in the Tapestry cohort. Methods The Tapestry Study, a large-scale, decentralized and clinical grade genomic sequencing study, conducted Exome+® sequencing at Helix, Inc. on 98,075 Mayo Clinic participants. After excluding patients with ≥10 pack-years of smoking and/or emphysema, we performed a GWAS of asthma exacerbations in 3,960 adults with asthma based on ICD codes. Cases (n = 1,059) were defined as individuals who experienced at least two exacerbations between 2020 and 2025; the remaining 2,901 individuals served as controls. Asthma exacerbations were defined by oral or systemic corticosteroid prescriptions with durations ranging from 3 days up to 28 days. Logistic regression models adjusted for age, sex, BMI, and the first six principal components of ancestry. Replication was performed in 1,273 asthma patients (439 African Americans and 834 non-Hispanic Whites) from the NHLBI-sponsored Severe Asthma Research Program (SARP 1-3), with adjustments for age, sex, BMI, and the first 10 principal components of ancestry. Results Based on a minor allele frequency (MAF) threshold of > 0.05, we analyzed 5,256,135 genome-wide variants. Genome-wide significant associations for asthma exacerbations (p < 5.0 × 10⁻8) included loci at a pseudogene in SMPD4P1 (rs680548, MAF = 0.39, OR = 1.37, p = 9.10 × 10⁻9). Additional variants showing suggestive associations included loci within ENDOD1 (rs638776, OR = 0.69, p = 1.20 × 10⁻6), PCSK2 (rs2206454, OR = 1.32, p = 3.53 × 10⁻6), RBFOX1 (rs17667084, OR = 0.60, p = 5.11 × 10⁻6), and CDC42EP3 (rs721660, OR = 1.29, p = 5.68 × 10⁻6). Subthreshold associations in Tapestry including loci at CDC42EP3 (rs4670213, OR = 1.27, p = 1.41x10-5, MAF = 0.36), and SDR39U1 (rs12587663,OR = 0.83, p = 5.86x10-4, MAF = 0.46) were replicated in the SARP longitudinal visit 5 with p-values of 8.523x10-3 (rs4670213, OR = 1.78), and 2.88x10-4 (rs12587663,OR=0.43), respectively. Conclusion Our GWAS for asthma exacerbations in the Tapestry cohort identified variant loci relevant to mitophagy mechanisms (CDC42EP3, PMID: 40783597), immune signaling and RNA splicing regulation (RBFOX1, PMID: 31301374), and neuromodulatory pathways (PCSK2, PMID: 32794589), all previously implicated in asthma pathophysiology. Collectively, these findings highlight potential novel genetic contributors to asthma exacerbations and provide a foundation for improving risk stratification and therapeutic targeting in asthma management. This abstract is funded by: This abstract is funded by: Mayo Clinic Funding: X01 from the NHLBI for TOPMed whole-genome sequencing and NIH grants R01 HL142992, HL69116, HL69167, HL69170, HL69174, U10 HL109164, U10 HL109257, U10 HL109146, U10 HL109172, U10 HL109250, U10 HL109250, U10 HL109250, U10 HL109168, U10 HL109152, U10 HL109086, U10 HL109046.

Abstract Introduction: A subset of patients with chronic lymphocytic leukemia (CLL) experiences disease progression from CLL to Richter transformation (RT), characterized by the occurrence of a high-grade lymphoma. The outcomes of patients with RT are typically dismal and identifying better therapeutic strategies remains critically important for these patients. To discern the biologic underpinnings driving this transformation, we studied immune changes and genetic evolution in samples from CLL patients who developed RT. Methods: We identified 18 patients with CLL who developed RT; and performed single cell (sc) RNA-sequencing(seq) on sorted B and T cells from peripheral blood mononuclear cells (PBMCs) at the time of RT (n=18) and at the earlier time of CLL diagnosis (n=11). The median time between the CLL sample and the one closest to RT was 4.3 years (range, 0.2-7.3); and patients received a median of 1 CLL directed treatment (range, 0 – 7). A total of 749,906 cells and 36,601 genes were captured during scRNAseq and 720,580 cells and 21,892 genes were retained for downstream analysis, respectively. Analysis was performed using the immunopipe (https://github.com/pwwang/immunopipe) that wraps around Seurat (v5.1.0) and identified enriched pathways using Enrichr service. The sc-T cell receptor (TCR)-repertoire(R)-seq data was analyzed and visualized using scplotter (). The whole exome sequencing (WES) analysis with tumor B-cells from same patients of CLL & RT stages were done using the Agilent sureselect Human all exon v5 + UTR capture kit. Reads from paired end sequencing were processed with the sarek pipeline. Results: Analysis of 720,580 cells from 11 CLL and 18 RT samples showed 16 distinct cell types- 3 B and 13 T cell subtypes. We observed an increasing trend of CD8 TEM [mean 0.548 vs. 0.384, p = 0.12 (Wilcoxon test)] but a decreasing trend of CD4 TCM (mean 0.236 vs. 0.337, p = 0.08) and B-memory (mean 0.0085 vs. 0.0195, p = 0.51) subsets in RT compared to their CLL stage. When comparing transcriptional changes, we found that among the top 18 highly variable genes, the expression of FP236383.3 (lncRNA) and TCL1A was significantly higher (p<0.0001) in the B cell subsets of patients with RT vs. CLL stage. Additionally, in T cell subsets—primarily dnT, γδT, and CD8/CD4 naïve cells—the expression of FP236383.3, NKG7, GNLY, and CCL5 was significantly higher (p<0.0001), while the expression of LDHB, GZMK, and IL32 was significantly lower (p<0.0001) in RT vs. CLL stage samples. Comparison of specific gene set enrichment analysis showed significantly increased expression of GZMB in gdT (p<0.001) and decreased expression of GZMA, GZMK, TIGIT and TCF7 in dnT and CD8 Naïve cells (p<0.001) for RT vs. CLL stage. The top 10 differentially expressed pathways in RT included oxidative phosphorylation, allograft rejection, IFN-α/-γ responses, TNF-α signaling via NF-κβ, complement, Myc targets V1, mTORC1, hypoxia, and p53 signaling across various B and T cell subsets compared to the CLL stage. Initial analysis of TCR-R data showed that patients with RT stage have lower number of unique clones (p<0.5) and reduced clonal diversity (p<0.5), and a higher number of large and hyperexpanded clones compared to their CLL stage. Additionally, WES data indicated that the most frequently mutated gene is TP53 in tumor B cells of both CLL and RT stages. However, a higher percentage of RT-stage samples carry mutations in MGA, ATM, and IRF2BP2 genes compared to their CLL-stage. WES data also showed a trend toward differing frequencies of copy number amplification at 19q12 (p=0.10) and 4q13.1 (p=0.52) between RT and CLL stages. Conclusions: The main findings of our study are as follows:The transcriptional landscape showed significant changes in multiple genes and pathway enrichment related to cell growth, differentiation and immune response in RT vs. CLL stage.A decreased fraction of CD4 TCM cells and reduced TCR clonal diversity were observed in samples at the RT stage and may indicate a poor prognosis for the patients at this stage. WES data of tumor B cells also indicated possible genetic alterations in the RT vs. CLL stage. Data from in-depth analyses on understanding of the dynamics of T cell response and additional experiments on these patients are expected to provide valuable insights into tumor genetic heterogeneity, transcriptional and immune regulation in blood-based immune subsets.

Introduction: The lack of understanding about the genetic and molecular basis of young stroke (<55 years of age) has hindered development of diagnostic, management, and therapeutic approaches. Whole exome sequencing (WES) can be leveraged as a tool for studying etiological processes underlying stroke. At the time of study inception, this is the single largest WES study in this unique young stroke population to uncover elusive stroke etiologies. Methods: The Mayo Clinic TAPESTRY database is a research program conducting whole exome sequencing of 100,000 participants. Of these, 406 participants met the criteria for young stroke (18-54 years of age) based on ICD-9 and ICD-10 diagnostic codes. Controls were TAPESTRY participants without a stroke-related ICD diagnosis. Cases were matched for age, gender, and race for a 1:4 case-control ratio. Ischemic stroke etiology was classified using TOAST classification. To uncover genetic mechanisms, a gene centric approach was used to identify variants in our cohort. A list of 66 genes known to be associated with young stroke was curated based on existing genetic databases. Various filtration strategies were used to identify potentially damaging variants (graphic 1). Results: There were 406 cases (363 ischemic and 43 hemorrhagic; 291 female; age at stroke onset 39.8±10.4; 368 were White) and 3,857 controls. TOAST classification yielded 20 patients with large artery atherosclerosis, 58 cardioembolic, 10 small vessel disease, 148 other determined (58 patent foramen ovale, 50 dissection, 48 hypercoagulable state, 11 Moyamoya, 11 vasculitis, 10 fibromuscular dysplasia), and 165 undetermined etiologies. Several individuals had multiple competing etiologies. Of cerebral hemorrhages, there were 6 aneurysmal, 6 arteriovenous malformation, 5 hypertensive, 2 substance abuse, 3 trauma, and 21 other etiologies. Across the entire cohort, 35,054 variants were found within the 66 genes identified from the literature. Filtration identified 48 potentially damaging variants across 30 genes in 44 patients (20 cryptogenic and 21 other: 12 patent foramen ovale, 5 vasculitis, 4 dissection, 4 cardioembolic) that were absent in controls. Conclusion: WES is a powerful tool that can be used to uncover potential genetic mechanisms in the young stroke population. An improved understanding of etiologies would improve diagnostics and pave the way for new preventive and therapeutic approaches, thus improving the overall care of young stroke patients.

770 Background: Pancreatic ductal adenocarcinoma (PDAC) is associated with poor prognosis and limited treatment options. A 584 bp deletion in CTRB2 , which impairs chymotrypsin B2 function, has been linked to increased PDAC risk. This study investigates the impact of this deletion on progression and survival outcomes in patients with PDAC. Methods: There were 633 patients in whom CTRB2ex6 deletion was genotyped; 263 patients received chemotherapy and had time to progression (TTP) information. TTP was calculated from diagnosis date until earliest of 1) progression date (event), 2) surgery, death, or next chemotherapy date (competing risk), or 3) date last known to be alive (censor). Overall survival (OS) was calculated for the cohort from diagnosis until date of death (event) or date last known to be alive (censor). TTP was analyzed using cumulative incidence methods accounting for competing risks; OS was analyzed using Kaplan-Meier methods. Cox proportional hazard models adjusted for age and sex were used to test associations and compute hazard ratios (HR) and 95% confidence intervals (CI). Results: Median TTP in those with CTRB2ex6 deletion negative carrier status (n=209) was 1.3 years vs 5.9 years for those with positive carrier status (n=54). Estimated 3-year progression rates were 54.2% for patients with no CTRB2ex6 deletions, 43.7% for patients with one deletion (n=48), and 33.3% for patients with two CTRB2ex6 deletions (n=6). TTP was not associated with CTRB2 deletion carrier status (HR: 0.8, 95% CI: 0.6-1.3). Median OS for CTRB2ex6 deletion negative carriers (n=514) was 1.0 year vs 0.9 years for positive carriers (n=119); within deletion carriers, heterozygous carriers (n=110) had a median OS of 0.9 years, and homozygous carriers (n=9) had a median OS of 2.6 years. OS was not associated with CTRB2 deletion carrier status (HR: 1.0, 95% CI: 0.9-1.3). Conclusions: Although we did not detect statistically significant differences in TTP and OS by CTRB2ex6 deletion carrier status, there was an interesting trend toward longer TTP and OS in the small number of patients who carried homozygous CTRB2ex6 deletions. Larger studies are warranted to validate these findings.

Abstract Aberrant activation of GLI transcription factors has been implicated in the pathogenesis of different tumor types including pancreatic ductal adenocarcinoma. However, the mechanistic link with established drivers of this disease remains in part elusive. In this study, using a new genetically engineered mouse model overexpressing constitutively active mouse form of GLI2 and a combination of genome-wide assays, we provide evidence of a novel mechanism underlying the interplay between KRAS, a major driver of pancreatic ductal adenocarcinoma development, and GLI2 to control oncogenic gene expression. These mice, also expressing KrasG12D, show significantly reduced median survival rate and accelerated tumorigenesis compared with the KrasG12D only expressing mice. Analysis of the mechanism using RNA sequencing demonstrate higher levels of GLI2 targets, particularly tumor growth–promoting genes, including Ccnd1, N-Myc, and Bcl2, in KrasG12D mutant cells. Furthermore, chromatin immunoprecipitation sequencing studies showed that in these cells KrasG12D increases the levels of trimethylation of lysine 4 of the histone 3 (H3K4me3) at the promoter of GLI2 targets without affecting significantly the levels of other major active chromatin marks. Importantly, Gli2 knockdown reduces H3K4me3 enrichment and gene expression induced by mutant Kras. In summary, we demonstrate that Gli2 plays a significant role in pancreatic carcinogenesis by acting as a downstream effector of KrasG12D to control gene expression.

Abstract We identified resistance mechanisms to abiraterone acetate/prednisone (AA/P) in patients with metastatic castration-resistant prostate cancer (mCRPC) in the Prostate Cancer Medically Optimized Genome-Enhanced Therapy (PROMOTE) study. We analyzed whole-exome sequencing (WES) and RNA-sequencing data from 83 patients with metastatic biopsies before (V1) and after 12 weeks of AA/P treatment (V2). Resistance was determined by time to treatment change (TTTC). At V2, 18 and 11 of 58 patients had either short-term (median 3.6 months; range 1.4–4.5) or long-term (median 29 months; range 23.5–41.7) responses, respectively. Nonresponders had low expression of TGFBR3 and increased activation of the Wnt pathway, cell cycle, upregulation of AR variants, both pre- and posttreatment, with further deletion of AR inhibitor CDK11B posttreatment. Deletion of androgen processing genes, HSD17B11, CYP19A1 were observed in nonresponders posttreatment. Genes involved in cell cycle, DNA repair, Wnt-signaling, and Aurora kinase pathways were differentially expressed between the responder and non-responder at V2. Activation of Wnt signaling in nonresponder and deactivation of MYC or its target genes in responders was detected via SCN loss, somatic mutations, and transcriptomics. Upregulation of genes in the AURKA pathway are consistent with the activation of MYC regulated genes in nonresponders. Several genes in the AKT1 axis had increased mutation rate in nonresponders. We also found evidence of resistance via PDCD1 overexpression in responders. Implications: Finally, we identified candidates drugs to reverse AA/P resistance: topoisomerase inhibitors and drugs targeting the cell cycle via the MYC/AURKA/AURKB/TOP2A and/or PI3K_AKT_MTOR pathways.

P-selectin is a cellular adhesion molecule that enhances procoagulant activity and activation of leukocyte integrins. P-selectin has been shown to be crucial in the development of atherosclerosis. Using a stratified random sample of 2406 participants in the Multi-Ethnic Study of Atherosclerosis (MESA) at Exam 2, we investigated the association of soluble P-selectin and cardiovascular disease. The Tobit model was used to investigate the association of P-selectin levels and coronary artery calcium (CAC) and Cox proportional hazards regression was used to assess the association between protein level and incident cardiovascular disease (CVD) events (myocardial infarction, stroke, angina, or other atherosclerotic death). Mean levels of soluble P-selectin differed significantly by ethnicity; 46±18 ng/mL African Americans, 58±20 ng/mL Caucasians, 43±15 ng/mL Chinese, and 61±20 ng/mL Hispanics (p<0.0001). CAC was present in 46% of African Americans, 61% of Caucasians, 52% of Chinese, and 53% of Hispanics. P-selectin levels were significantly associated with CAC in Caucasians (p < 0.001) and African Americans (p = 0.025), see Table. This association was robust to adjustment for traditional risk factors and a significant interaction of P-selectin and ethnicity was observed (p= 0.02). Over an average of 6.5 years of follow-up there were 111 incident CVD events after MESA exam 2, 26 in African Americans, 33 in Caucasians, 16 in Chinese, and 36 in Hispanics. The risk of CVD events increased 34% per standard deviation increase in P-selectin levels only in Caucasians (p = 0.007) and this finding was consistent after adjustment for CVD risk factors (HR = 1.30, p = 0.017). We have shown in this large cohort that mean circulating P-selectin levels differed by ethnicity. Despite similar percentages of subjects with CAC in each ethnic group, P-selectin levels were significantly associated with CAC only in African Americans and Caucasians. Furthermore, in Caucasians only, P-selectin levels were associated with incident CVD events. Association of P-selectin, per SD (20 ng/mL) for the population and cardiovascular disease Outcome African Americans n = 579 Caucasians n = 620 Chinese n = 602 Hispanics n = 605 CAC, Agatston Score Beta (S.E.) P-value Beta (S.E.) P-value Beta (S.E.) P-value Beta (S.E.) P-value Model 1 85.9 (38.4) 0.025 99.4 (28.7) <0.001 8.2 (33.3) 0.81 32.1 (27.3) 0.24 Model 2 79.3 (39.8) 0.046 94.7 (29.1) 0.001 −19.5 (33.6) 0.56 27.8 (27.6) 0.31 CVD Events HR (CI) HR (CI) HR (CI) HR (CI) Model 1 1.08 (0.75, 1.57) 0.67 1.34 (1.08, 1.65) 0.007 1.09 (0.57, 2.08) 0.80 0.99 (0.71, 1.39) 0.97 Model 2 1.10 (0.76, 1.58) 0.61 1.30 (1.05, 1.62) 0.017 0.82 (0.38, 1.78) 0.62 1.00 (0.71, 1.43) 0.98 Model 1 = age and sex ; Model 2 = age, sex, BMI, smoking status, LDL, HDL, triglycerides, hypertension, diabetes