Citations

Abstract Cerebrospinal fluid (CSF) is the biofluid most proximal to the central nervous system (CNS), yet its proteomic organization across neurological disease remains unresolved. We assembled a meta-atlas of 35 CSF proteomic datasets (2016–2026) comprising 10,560 quality-controlled samples spanning neurodegenerative, neuroinflammatory, neuro-oncological, infectious and structural-injury contexts. Meta-analysis of cross-cohort detection defined a 747-protein core mass-spectrometry-accessible CSF proteome that generalizes to held-out cohorts. We derived thirteen consensus meta-programs that transfer across independent cohorts and span the major neural, barrier, immune and plasma compartments. Unexpectedly, these programs are dominated by a single axis of variation opposing complement- and plasma-associated inflammation against neuronal-synaptic homeostasis. This injury–homeostasis axis is diagnosis-spanning and tracks clinical severity and disease activity across diverse neurological conditions. We show that a two-protein C5–NRCAM model can approximate the injury–homeostasis axis, bridging deep proteomic profiling toward bedside assessment.

Abstract The liver is organized into tightly regulated zones with distinct metabolic functions but zonation erodes with age. Cellular senescence contributes to aging and liver diseases, however, its impact on aging biology is ill-defined. As part of The Cellular Senescence Network Consortium, we used multiple spatial transcriptomics approaches (GeoMx, Visium, CosMx) with snRNA-seq to profile senescence signatures, zonation markers, and metabolic pathways in livers from wild-type (WT) mice of multiple ages. We observed a loss of canonical zone signatures in aged mouse livers characterized by “expansion” of midlobular (zone 2) marker gene expression, accompanied by diminished expression of zone 3 marker genes by middle-age (18 months), indicative of loss of cell identity. Multiple analytic approaches identified distinct age-, zone- and sex-specific senescence signatures, which were significantly associated with zonation markers changes. This was recapitulated in Ercc1 mutant models of accelerated senescence, supporting a causal role of senescent cells in liver aging. A “no-zone” hepatocyte-like cluster expanded with age and with the strongest Senescence-Associated Secretory Phenotype (SASP) profile. Gene expression profiles from senescent hepatocytes implicate decreased WNT signaling and increased BMP as contributing to age-related loss of zonation. Together, these data elucidate the role of senescent cells in driving aging biology in non-diseased liver through disruption of cell:cell signaling and the loss of metabolic and cell identity gene expression necessary for hepatocyte function.

G protein-coupled receptors (GPCRs) are among the most common drug targets, with numerous agonists and antagonists approved for clinical use. More recently, it has been appreciated that drugs can also modulate the action of natural agonists of these receptors, thus providing unique clinical advantages. Here, we describe the identification and characterization of a small molecule positive allosteric modulator (PAM) of secretin action at the class B G protein-coupled secretin receptor. This phenyl(trifluoromethyl)-pyrimidine can occupy the secretin receptor without stimulating its internalization, yet priming it to enhance both the potency and efficacy of the action of natural secretin. This is also shown to exhibit its effects on cells expressing low numbers of these receptors, without enhancing the effects of other structurally related hormones acting at other class B GPCRs. The mechanism responsible for this PAM effect is the slowing of the off-rate of receptor-bound secretin. This compound can serve as a lead to the development of other drugs that enhance the action of natural endogenous secretin and can be utilized to explore the potential therapeutic utility of such compounds.

Abstract Spatial transcriptomics has been widely applied to study the spatial distribution of cell types, cell states, and specific gene expression in tissue samples. However, we show that there is a prevalent transcript leakage problem in spatial transcriptomics data, where transcripts expressed by a cell diffuse to its neighborhood and are recurrently detected in the nearby cells. By analyzing published data sets, we show that this problem is general across data produced from different tissues and different species using different imaging-based and sequencing-based spatial transcriptomics platforms. It affects both upstream tasks such as expression quantification as well as downstream tasks such as cell-type annotation and detection of spatially-dependent gene expression. To tackle the transcript leakage problem, we propose a reference-free Bayesian model-based method, DeLeakage, which cleans up the data much more effectively than existing denoising methods. DeLeakage also improves cell-type annotation and avoids false detection of spatially dependent expression.

ABSTRACT Biological aging greatly impacts the body’s ability to handle glucose, and represents a major risk factor the development and progression of type 2 diabetes (T2D). Nonetheless, despite advances in cellular senescence research and the development of new senolytic therapies, the heterogeneity of cellular senescence in the human endocrine pancreas, as well as its roles in normal aging, remains to be elucidated at the single-cell level. Here, we performed single-cell–resolved spatial proteomics and transcriptomics on intact pancreas from 26 donors (ages 20–80) and multiplexed single-cell RNA sequencing and functional assays on dispersed islets from 14 donors (ages 34–69). We identify two discrete SnC subpopulations distinguished by relative expression of CDKN1A and CDKN2A. CDKN1A ⁺ senescent cells (SnCs) exhibit loss of β-cell identity, impaired insulin secretion, and a proinflammatory SASP associated with increased islet immune infiltration. In contrast, CDKN2A ⁺ SnCs retain transcriptional identity and functional competence, with lower inflammatory signaling. Together, these findings identify heterogeneous and functionally divergent senotypes in the human pancreas, distinguishing an adaptive ( CDKN2A ⁺) from a maladaptive ( CDKN1A ⁺) senescence program, thus providing a mechanism-guided framework for senescence-targeted therapies in T2D.

Summary Cancer genome sequencing is essential for understanding tumor evolution and advancing precision medicine. 1 However, reference gaps and germline variants obscure detection of small and large somatic variants and methylation in repetitive regions. 1–3 It is common for tumor cells to gain or lose chromosome arms due to somatic structural changes that occur inside highly repetitive satellite DNA sequences in the centromeres. 4 To identify the full spectrum of somatic variants, including complex rearrangements, we construct and curate near-complete, haplotype-resolved assemblies of the most recent common ancestor of an early-passage broadly-consented hypodiploid pancreatic cancer cell line and matched normal tissues. The tumor assembly completely recapitulates all 35 tumor chromosomes observed with karyotyping, with multiple translocation-induced hybrid chromosomes. The hybrid chromosomes contain putative functional dicentric and fused centromeres, nested foldback inversions causing 14 breakpoints with a haplotype switch in a single event, and centromeric satellite tandem duplications up to 136 kbp. Direct comparison of tumor and normal assembly haplotypes uncovers >7,000 variants altering >1 Mbp of sequence in repetitive regions that have been hidden by reference gaps and germline variants. 44 % of somatic small variants change representation because they alter germline variants on GRCh38, impacting mutational signatures and kataegis/omikli clusters. Most somatic LINE insertions originate from two hypomethylated non-reference germline LINE insertions, highlighting their impact on insertion mutation burden. These assemblies demonstrate that centromeric, acrocentric, and telomeric regions conventionally excluded from analysis harbor extensive somatic and epigenetic changes. Resolving complete tumor genomes enables a deeper understanding of cancer structural plasticity and the endpoints of breakage-fusion-bridge cycles. These assembled, curated paired normal-tumor benchmarks will serve as a critical foundation for developing future algorithms to characterize the most intractable regions of cancer genomes.

Abstract Inherited susceptibility plays a critical role in prostate cancer (PCa) risk. Using large biobank and case-control datasets, we evaluated the contribution of both common variants and rare germline pathogenic variants (PVs) to overall PCa risk. The Prostate Cancer Exome Sequencing Consortium currently includes 427,388 male participants (51,452 PCa cases and 375,936 controls) with whole-exome sequencing data from ten biobanks and studies: UK Biobank (14,669 cases/195,600 controls), All of Us Research Program (7,577/75,226), African Ancestry Prostate Cancer Consortium (7,176/4,675), Mayo Clinic Biobank (6,031/15,084), Mass General Brigham Biobank (3,393/14,095), Geisinger’s MyCode Community Health Initiative (3,026/15,130), UCLA ATLAS Precision Medicine Biobank (2,850/16,904), Penn Medicine Biobank (2,598/16,255), Colorado Center for Personalized Medicine (2,269/13,399), and Malmo Diet and Cancer (1,863/9,568). Based on self-reported race/ethnicity and estimated genetic ancestry, the cases comprise approximately 79% European, 18% African, and 3% other ancestry populations. Single-variant association analyses tested all variants on chromosomes 1-22 and X with a minor allele count ≥ 5. In gene-based analyses, PVs were defined as rare variants (minor allele frequency [MAF] < 1% in controls) that had either a Variant Effect Predictor (VEP) impact score of “high” or a pathogenic or likely pathogenic ClinVar classification. Associations were estimated using Firth logistic regression, adjusting for age and the top ten genetic principal components. Results from individual studies were combined using fixed-effect meta-analysis. In single-variant association analyses, 496 variants reached genome-wide significance (p<5×10-8; MAF>0.02%). Among these, 458 (92%) variants mapped to previously known risk regions, including three rare PVs in HOXB13 (rs138213197), CHEK2 (rs555607708), and FAM111A (rs533676902). Characterization of the remaining 38 variants is ongoing. Gene-based analyses identified significant associations (p<2.4×10-6) for eight genes: HOXB13 (OR=3.7, 95% CI=3.3-4.2), BRCA2 (OR=2.0, 95% CI=1.7-2.3), CHEK2 (OR=1.6, 95% CI=1.5-1.8), ATM (OR=1.6, 95% CI=1.4-1.9), FAM111A (OR=1.4, 95% CI=1.3-1.5), BIK (OR=1.4, 95% CI=1.2-1.6), SAMHD1 (OR=2.1, 95% CI=1.6-2.7), and SMOC2 (OR=3.2, 95% CI=2.0-5.1). All genes except SMOC2 have been previously implicated in PCa susceptibility. Among cancer predisposition and DNA repair genes, nominal associations were also observed for XRCC2 (OR=1.6, 95% CIs=1.2-2.3) and BRCA1 (OR=1.2, 95% CI=1.0-1.4), whereas the association was not significant for PALB2 (OR=1.2, 95% CI=0.9-1.5). These findings reinforce the role of rare germline PVs, particularly in cancer predisposition and DNA repair genes, in PCa susceptibility. As additional studies are incorporated into the Consortium, we expect this work to provide a more comprehensive characterization of the genetic architecture of PCa. Citation Format: Yifan Zhang, Shuyan Cheng, Nicholas Boddicker, Matthew Lebo, Alexander S. Berry, Roni Haas, Ryan Hausler, Tokhir Dadaev, Heena Desai, Alex A. Rodriguez, Ravi K. Madduri, Andrew Hill, Xin Sheng, Susan M. Gundell, Mine Cicek, Penn Medicine Biobank, Olle Melander, Chris R. Gignoux, Isla P. Garraway, Bogdan Pasaniuc, Paul C. Boutros, Matt Oetjens, Adam S. Kibel, Robert J. Klein, Zsofia Kote-Jarai, Fergus J. Couch, Kara N. Maxwell, Burcu F. Darst, David V. Conti, Christopher A. Haiman, Fei Chen. Genetic risk of prostate cancer: Insights from the Prostate Cancer Sequencing Consortium [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2026; Part 2 (Late-Breaking, Clinical Trial, and Invited Abstracts); 2026 Apr 17-22; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2026;86(8_Suppl):Abstract nr LB390.

Abstract Background: Cholangiocarcinoma (CCA) is an aggressive hepatobiliary malignancy with limited treatment options. CCA exhibits a distinct desmoplastic and tolerogenic tumor immune microenvironment (TIME) characterized by elevated bile acid concentrations resulting from obstructive cholestasis and dysregulated bile acid metabolism. TGR5 is an essential component of bile acid signaling. We hypothesize that bile acids promote CCA progression by establishing an immunosuppressive TIME through TGR5 signaling, which can be therapeutically targeted using SBI-364, a novel potent small molecule TGR5 inhibitor. Methods: Single-cell RNA sequencing (scRNA-seq) of immune cells from human intrahepatic (iCCA) and perihilar (pCCA) CCA tumors were conducted. Therapeutic studies were performed using syngeneic murine models of iCCA (PMID: 38458319), a novel syngeneic model of pCCA, and patient-derived xenograft (PDX) models. In collaboration with Sanford Burnham Prebys, SBI-364, a first in-class competitive TGR5 antagonist, was assessed in vivo. Additionally, we performed co-culture experiments with bile acid +/- SBI-364-treated murine bone marrow-derived monocytic myeloid-derived suppressor cells (M-MDSCs) and CD8+ T cells to assess functional immunosuppressive effects. Results: Human scRNA-seq revealed significantly increased TGR5 (GPBAR1) expression predominantly in myeloid cell populations within the CCA TIME. SBI-364 treatment significantly reduced tumor burden in both syngeneic orthotopic iCCA models and a novel pCCA model that recapitulates the anatomic characteristics of perihilar CCA. In contrast, NOD/SCID mice bearing PDX535 (predicted sensitive) and PDX136 (predicted resistant), SBI-364 treatment did not achieve tumor reduction, suggesting that TGR5 inhibition predominantly attenuates tumor progression through modulation of bile acid-immune cell interactions. SBI-364 treated murine tumors had a significant reduction in immunosuppressive macrophages expressing interleukin-10 and M-MDSCs compared to vehicle treated tumors. Furthermore, bile acid-treated M-MDSCs exhibited enhanced suppression of CD8+ T cell proliferation and function, an effect that was abrogated by SBI-364 co-treatment. Conclusion: We demonstrate that TGR5 inhibition represents a promising therapeutic strategy across multiple immunocompetent preclinical CCA models. We demonstrate that the bile-acid TGR5 axis promotes CCA progression via immunosuppressive myeloid cells. Specifically, TGR5 signaling potentiates the immunosuppressive capacity of M-MDSCs, contributing to an immunosuppressive niche that facilitates tumor progression. These results establish TGR5 as a rational therapeutic target for overcoming immune evasion in cholangiocarcinoma. Citation Format: Enis Hikmet Ozmert, Binbin Li, Jingchun Yang, Hannah Elizabeth Stumpf, Jessica R. Willhite, Danielle Carlson, Jennifer L. Tomlinson, Jon Covel, Raghuveer Ramachandra, Steven Olson, Rory L. Smoot, Sumera Ilyas. Therapeutic targeting of TGR5 abrogates bile acid driven myeloid cell immunosuppression in cholangiocarcinoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2026; Part 1 (Regular Abstracts); 2026 Apr 17-22; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2026;86(7 Suppl):Abstract nr 2896.

Abstract Mitochondrial creatine kinase (MtCK) is a key enzyme in energy buffering and homeostasis in cells. It catalyzes transfer of phosphoryl group from ATP to creatine. Overexpression of MtCK occurs in many cancer cells to meet elevated energy demands, which is associated with poor prognosis. This suggests that MtCK may be a promising target for cancer therapeutics. We sought to discover first-in-class selective inhibitors of MtCK with diverse mechanisms of action using high-throughput screening, biochemical characterization and cryo-EM studies. Through these studies, we identified diverse types of compounds that modulate activity of MtCK in vitro , including fast-equilibrium and time-dependent orthosteric and allosteric inhibitors. Select hits were subjected to in vitro enzymatic and binding assays to assess MtCK inhibition and binding. A subset of inhibitors was advanced into structural studies using cryo-EM resulting in the molecular structure of MtCK with and without bound substrates in complex with an allosteric uncompetitive inhibitor that we discovered. These studies identified compound’s unique binding pocket on MtCK and established the molecular steps manifesting in the apparent uncompetitive mode of inhibition. We demonstrate that the compound stabilizes an active site loop in a closed conformation restricting access of the creatine substrate to and the release of phosphocreatine product from its binding pocket. These findings establish chemical tools suitable for validation of MtCK as a promising breast cancer target with high therapeutic potential and build a foundation for future structure-guided optimization of the hit compounds of MtCK we identified and de novo rational design of novel MtCK inhibitors. Abstract Figure

Abstract Background Sex, education and race/ethnicity are all associated with risk of Alzheimer's disease dementia. Here, we assess the effects of self‐reported sex, educational attainment and race/ethnicity on amyloid‐positivity, and tau‐PET‐positivity in 12,048 (7,394 cognitively unimpaired [CU], 2,177 MCI, and 2,477 dementia) individuals from 42 cohorts worldwide. Method Logistic generalized estimating equations were used to estimate frequency of amyloid‐positivity (using cohort‐specific thresholds for amyloid‐PET [84%] or CSF) and tau‐PET‐positivity (cohort‐specific thresholds of 2SD above mean temporal uptake in amyloid‐negative controls). We assessed: i) sex and APOEε4 ( N = 10,098) associations, to complement earlier findings of a higher frequency of tau‐positivity in females, ii) effects of lower/higher education ( N = 10,970; cohort‐specific median‐split), and iii) effects of race/ethnicity (non‐Hispanic White [hereafter: White], N = 4880; Asian, N = 116; Black or African‐American [hereafter: Black], N = 353; Hispanic, N = 356, only from Northern‐American cohorts). Outcomes were frequency of amyloid‐positivity in CU individuals only, and tau‐PET‐positivity in both amyloid‐positive (AB+) CU and cognitively impaired (CI, i.e. MCI and dementia) individuals. Interaction effects on the relationship between age and amyloid/tau‐positivity were assessed and only retained in the models when significant. Result Female sex was associated with an APOEε4 ‐independent increased frequency of amyloid‐positivity (β=0.51[0.22], p = 0.02) in CU and increase of tau‐positivity in both AB+CU (β=0.27[0.08]) and AB+CI (β=0.37[0.08], both p <0.01). Remarkably, tau‐positivity frequencies of female APOEε4 non‐carriers were equivalent to male APOEε4 carriers in AB+CI (Figure 1). No significant sex* APOE interactions were observed. In CU, higher education was associated with lower amyloid‐positivity frequency (β=‐0.12[0.05], p = 0.02). In contrast, among AB+CU, there was an age*education interaction effect that indicated more pronounced age effects on tau‐positivity in individuals with higher education (age*education:β interaction =0.03[0.01], p <0.01). There were no education effects in AB+CI (Figure 2). In CU, an age*race/ethnicity interaction effect was observed across all non‐White groups compared to White (Hispanic:β interaction =‐0.05[0.01], p <0.01; Black:‐0.04[0.01], p <0.01; Asian:‐0.02[0.01], p = 0.04). This suggests that the impact of age on amyloid‐positivity was less pronounced in non‐White groups. Furthermore, in AB+CI, Hispanic ethnicity was related to higher tau‐positivity frequency than White (β=0.51[0.22], p = 0.02; Figure 3). Conclusion In this multi‐center initiative comprised of clinical and community‐based cohorts, we observed that self‐reported sex, educational attainment and race/ethnicity were related to positivity‐frequencies of Alzheimer's disease pathology.