Back

Journal of Clinical Investigation

American Society for Clinical Investigation

Preprints posted in the last 7 days, ranked by how well they match Journal of Clinical Investigation's content profile, based on 164 papers previously published here. The average preprint has a 0.20% match score for this journal, so anything above that is already an above-average fit.

1
OCA-B/Pou2af1 Expression in T Cells Promotes PD-1 Blockade-Induced Autoimmunity but is Dispensable for Anti-Tumor Immunity

Du, J.; Manna, A. K.; Medina-Serpas, M. A.; Hughes, E. P.; Bisoma, P.; Evason, K. J.; Young, A.; Wilson, W. D.; Brusko, T.; Farahat, A. A.; Tantin, D.

2026-04-16 immunology 10.1101/2025.10.22.683978 medRxiv
Top 0.2%
8.3%
Show abstract

The transcription coregulator OCA-B promotes CD4+ T cell memory recall responses and autoimmunity. OCA-B T cell deletion prevents spontaneous type-1 diabetes (T1D) onset in non-obese diabetic (NOD) mice and blunts T1D in a subset of more aggressive models. However, the role of OCA-B in diabetes induced by treatment with immune checkpoint inhibitors (ICIs), and the role of OCA-B in the control of tumors with and without ICI treatment, has not been studied. Here we show that islet and pancreatic lymph node T cells from T1D individuals express measurable POU2AF1 mRNA. Deletion of OCA-B in T cells fully insulates 8-week-old non-obese diabetic (NOD) mice against ICI-induced diabetes and partially protects 12-week-old mice. Salivary and lacrimal gland infiltration and inflammation were also reduced. Protection was associated with a block in the differentiation of progenitor exhausted CD8+ T cells (TPEX) into terminally exhausted CD8+ T cells (TEX). We show that OCA-B T cell loss preserves anti-tumor immune responses following PD-1 blockade in different tumors and mouse strains. These findings point to a potential therapeutic window in which pharmaceuticals targeting OCA-B could be used to block the emergence of both spontaneous and ICI-induced autoimmunity while sparing anti-tumor immunity. We develop first-in-class small molecule inhibitors of Oct1/OCA-B transcription complexes and show that administration into NOD mice also blocks diabetes emergence following PD-1 blockade. These results identify OCA-B as a promising therapeutic target for the prevention of autoimmunity and immune-related adverse events (irAEs).

2
Vaccine-induced antibody and T cell responses in children with acute lymphoblastic leukemia

Shapiro, J. R.; Dorogy, A.; Science, M.; Gupta, S.; Alexander, S.; Bolotin, S.; Watts, T. H.

2026-04-12 oncology 10.64898/2026.04.10.26350531 medRxiv
Top 0.4%
6.1%
Show abstract

Children with acute lymphoblastic leukemia (ALL) are treated with multiagent chemotherapy that causes profound changes to the immune system. There are limited data on how disease and therapy impact antigen-specific immune memory, leading to inconsistent guidelines on best practices for revaccination of this population. Here, to inform vaccine guidance, we investigated whether immunity derived from routine childhood measles and varicella zoster virus (VZV) vaccines is maintained during and after therapy for childhood ALL. We report that antibodies against measles and VZV were significantly reduced in children with ALL (n=45) compared to healthy controls (n=13), particularly in older children in whom a longer time had passed since their most recent vaccine dose. However, the avidity of the measles and VZV-specific antibodies was indistinguishable between groups. Despite changes to the composition of the T cell compartment, both overall and antigen-specific T cell function were preserved in children with ALL. These data provide compelling evidence for revaccination of children following ALL treatment. Intact T cell responses suggest that post-treatment revaccination would be effective.

3
Molecular signature of pediatric B-ALL determines outcomes post CD19 CAR-T cell therapy

Oszer, A.; Pastorczak, A.; Urbanska, Z.; Miarka, K.; Marschollek, P.; Richert-Przygonska, M.; Mielcarek-Siedziuk, M.; Baggott, C.; Schultz, L.; Moon, J.; Aftandilian, C.; Styczynski, J.; Kalwak, K.; Mlynarski, W.; Davis, K. L.

2026-04-13 oncology 10.64898/2026.04.11.26350681 medRxiv
Top 0.5%
4.8%
Show abstract

Chimeric antigen receptor T-cell (CAR-T) therapy targeting CD19 has transformed outcomes for children with relapsed or refractory (R/R) B-cell acute lymphoblastic leukemia (B-ALL), yet the influence of molecular subtype on outcomes remains unclear. We evaluated the impact of cytogenetic and molecular signatures on complete response (CR), overall survival (OS), and leukemia-free survival (LFS) after CD19 CAR-T therapy in eighty-six pediatric patients with R/R B-ALL treated with tisagenlecleucel. CR was assessed 30 days after infusion. Cytogenetic data were available for 84 patients and molecular profiling for 62. Survival analyses included 72 patients who received CD19 CAR-T as the sole cellular therapy. Seventy-seven patients achieved CR (89.5%). Pre-infusion bone marrow blasts of [&ge;]20% were associated with lower CR rates (53.8% vs 95.9%, p<0.0001) and significantly reduced OS and LFS (both p<0.0001). Among molecular markers, RAS mutations correlated with inferior OS (p=0.0222) and LFS (0.0402). In multivariate analysis, bone marrow blasts >20% and RAS mutations independently predicted inferior OS. Post CAR-T, CD19 negative relapses showed almost twice higher prevalence of RAS mutations (66% vs 37.5%). These findings highlight RAS mutations as a key molecular predictor of outcome after CD19 CAR-T therapy and suggest emergence of unique risk stratification for patients receiving CD19-targeting therapy.

4
Gonadotropin-releasing hormone antagonism reduces paedophilic interest through increased cerebellar activity.

Mannfolk, C.; Ertl, N.; Jayasena, C. N.; Liberg, B.; Wall, M. B.; Comninos, A. N.; Rahm, C.

2026-04-13 sexual and reproductive health 10.64898/2026.04.12.26350231 medRxiv
Top 0.6%
4.3%
Show abstract

Mechanistic understanding and biomarkers of gonadotropin-releasing hormone antagonist treatment effect in paedophilic disorder are absent but may enhance outcomes and reduce sexual-offending risk. 52 help-seeking self-referred Swedish men with paedophilic disorder enrolled in a double-blinded, placebo-controlled, randomized clinical trial. Participants underwent task-based fMRI before, and two weeks after, subcutaneous injection of 120mg of degarelix or equal volume of placebo. fMRI blood-oxygen-level-dependent activation was compared between child and adult (child>adult) stimuli in task-derived regions of interest. Primary outcome was within region-of-interest child>adult activation change, whereas secondary outcomes correlated region-of-interest child>adult activation change to change in clinical measurements of risk, paedophilic interest, sexual preoccupation, hyper- and hyposexuality. 19 degarelix and 22 placebo participants had sufficient fMRI data quality. Reductions in paedophilic interest were strongly correlated with increased child>adult cerebellar (vermis) region-of-interest activation following degarelix (r=-0.740, p<0.001) but not placebo (r=0.183, p=0.41; between-group correlation coefficient z=3.347, p<0.001). Treatment did not significantly change child>adult region-of-interest activity. Post hoc analysis indicated that baseline autism symptoms correlated with degarelix-induced changes in paedophilic interest (r=0.717, p<0.001; between-group correlation coefficient z=2.958, p=0.003) and cerebellar activation (r=-0.581, p=0.01; between-group correlation coefficient z=-1.930, p=0.05). Increased child>adult cerebellar activation was associated with degarelix-induced reductions of paedophilic interest, suggesting cerebellar activity as mechanistically important to, and a prospective biomarker of, degarelix treatment effect. Additionally, autism symptoms may inform treatment prediction. Together, these findings have mechanistic and clinical implications for degarelix treatment of paedophilic disorder. EU clinical trials register identifier: 2014-000647-32 https://www.clinicaltrialsregister.eu/ctr-search/trial/2014-000647-32/SE, registered on 05/06/2014.

5
Genome-Wide Association Analysis of Tic Disorders Reveals 6 Independent Risk Loci and Highlights Tic-Associated Cell Types and Brain Circuitry

Yu, D.; Strom, N. I.; Gerring, Z. F.; Topaloudi, A.; Halvorsen, M. W.; Shekhar, S.; Miller-Fleming, T. W.; Tang, M.; Porras, L. M.; Ivankovic, F.; Mahjani, B.; Palviainen, T.; Corfield, E. C.; Androutsos, C.; Apter, A.; Ask, H.; Baglioni, V.; Ball, J.; Barr, C. L.; Barta, C.; Basha, E.; Batterson, J. R.; Benaroya-Milshtein, N.; Benarroch, F.; Boomsma, D. I.; Borglum, A. D.; Budman, C. L.; Buitelaar, J. K.; Buse, J.; Bybjerg-Grauholm, J.; Cardona, F.; Cath, D. C.; Cavallari, L. H.; Cheon, K.-A.; Coffey, B. J.; Dahl, N.; Depienne, C.; Dietrich, A.; Domenech, L.; Drineas, P.; Einarsson, G.; Elste

2026-04-13 genetic and genomic medicine 10.64898/2026.04.09.26350245 medRxiv
Top 0.8%
4.0%
Show abstract

Tourette Syndrome and other tic disorders (TD) are common, highly heritable neurodevelopmental conditions with complex genetic architectures. We conducted a genome-wide association study of 13,247 TD cases and 536,217 European ancestry controls and identified six independent genome-wide significant loci, including a pleiotropic signal at 3p21 shared with attention-deficit/hyperactivity disorder, among other traits. Gene prioritization highlighted 20 genes, including PCDH9, HCN1, NCKIPSD, WDR6, DALRD3, and CELSR3. Integrative analyses provide genetic support for the role of cortico-striato-thalamo-cortical circuits in TD pathophysiology and further localize TD genetic risk to specific cell types, including dopamine D1- and D2-receptor-positive medium spiny neurons, cortical pyramidal neurons, and oligodendrocyte-lineage cells. We further demonstrate extensive genetic correlations with neurodevelopmental and psychiatric traits, but not with neurological disorders. These findings advance our understanding of the genetic basis of TD, pinpointing specific genes and cell types that drive pathophysiology and providing a foundation for future mechanistic studies.

6
Multi-Ancestry Epigenome-Wide Meta-Analysis Identifies Novel Bulk and Cell-Type-Specific Epigenetic Markers of Asthma with Severe Exacerbations

Perez-Garcia, J.; Martin-Gonzalez, E.; Chen, Z. J.; Martin-Almeida, M.; Witonsky, J.; Gorla, A.; Eng, C.; Lorenzo-Diaz, F.; Bozack, A. K.; Elhawary, J.; Hu, D.; Huntsman, S.; Gonzalez-Perez, R.; Hernandez-Perez, J. M.; Poza-Guedes, P.; Mederos-Luis, E.; Sanchez-Machin, I.; Rodriguez-Santana, J.; Villar, J.; Rifas-Shiman, S. L.; Hivert, M.-F.; Oken, E.; Gold, D. R.; Ziv, E.; Rahmani, E.; Gonzalez Burchard, E.; Cardenas, A.; Pino-Yanes, M.

2026-04-18 allergy and immunology 10.64898/2026.04.16.26350345 medRxiv
Top 0.9%
3.9%
Show abstract

Background: Extreme-phenotype comparisons allowed the discovery of novel asthma genetic risk loci. However, this approach remains unexplored in epigenome-wide association studies (EWAS). We aimed to identify bulk and cell-specific methylation markers of asthma with severe exacerbations across diverse ancestry groups. Methods: We conducted a meta-EWAS of 739,543 CpGs in whole blood among 1,192 African American and Latino pediatric populations, comparing non-asthmatics and asthma exacerbators. Genome-wide CpGs were followed up for replication in a meta-analysis across 1,516 ethnically diverse participants and in a cross-tissue evaluation of 393 nasal samples. We conducted differentially methylated region (DMRs), cell-type-deconvoluted, and quantitative trait loci analyses (whole-genome sequencing n=1,668; RNA-seq n=1,209). We examined enrichment in traits, pathways, and druggable genes, and analyzed DNAm predictors of plasma proteins and aging. Results: DNAm at 505 CpGs and 119 DMRs in whole blood were associated with asthma exacerbations (p<9x10-8, {lambda}=1.05). We replicated 25 CpGs in blood cells, cross-validated 7 in nasal samples, and detected 42 cell-specific DNAm markers mainly driven by T cells. DNAm at 134 CpGs was associated with gene expression in whole blood, including 118 associations with T-cell receptor genes, and 446 CpGs were regulated by [&ge;]1 genetic variant. We found enrichment for previous associations with environmental exposures, immune disorders, immune and inflammatory pathways, and druggable genes by developmental drugs. 21 methylation-predicted plasma proteins, involved in host defense, and one lung aging clock were associated with asthma exacerbations. Conclusions: The first meta-EWAS of extreme asthma phenotypes identified hundreds of novel DNAm markers, suggesting novel methylation biomarkers and candidate drugs for asthma and supporting the role of T cells.

7
Comprehensive Immunophenotyping of Monocytes and Dendritic Cells Suggests Distinct Pathophysiology in Chronic Fatigue Syndrome and Long COVID

Petrov, S. I.; Bozhkova, M.; Ivanovska, M.; Kalfova, T.; Dudova, D.; Todorova, Y.; Dimitrova, R.; Murdjeva, M.; Taskov, H.; Nikolova, M.; Maes, M.

2026-04-12 allergy and immunology 10.64898/2026.04.10.26350613 medRxiv
Top 1%
3.6%
Show abstract

Myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS) and long COVID are complex chronic conditions that often follow infectious triggers with overlapping clinical features but poorly defined pathophysiological relationships. This study aimed to identify disease-specific immune signatures through multiparameter immunophenotyping of monocytes, dendritic cells, and T-cell subsets. A total of 207 participants were included (ME/CFS: n = 103; long COVID: n = 63; healthy controls: n = 41). Peripheral blood mononuclear cells were analyzed using multiparameter flow cytometry. Statistical analyses included non-parametric testing, age-adjusted ANCOVA, correlation network analysis, and principal component analysis (PCA). Long COVID was characterized by increased M2-like monocyte polarization, elevated CD80 expression across monocyte subsets, expansion of dendritic cells, and reduced expression of activation markers, indicating persistent immune activation with features of immune exhaustion. In contrast, ME/CFS exhibited reduced costimulatory molecule expression, impaired CCR7-mediated immune cell trafficking, and less coordinated activation patterns, consistent with a state of immune suppression. Correlation network analysis revealed more extensive and integrated immune interactions in long COVID, while PCA identified distinct immunophenotypic components and enabled moderate discrimination between the two conditions. These findings demonstrate that ME/CFS and long COVID are characterized by distinct immune profiles, supporting the concept of divergent immunopathological mechanisms. The identified signatures may contribute to biomarker development and guide targeted therapeutic approaches.

8
NMF Deconvolution of a High-ROS Transcriptional Program Uncovers mTOR-Dependent Therapeutic Sensitivity in Stomach Adenocarcinoma

Roy, R.; Patnaik, J.; Chakraborty, A.; Patnaik, S.; Parija, T.

2026-04-16 oncology 10.64898/2026.04.12.26350699 medRxiv
Top 1%
3.0%
Show abstract

Background: Stomach adenocarcinoma is driven by heterogeneity, limiting therapeutic success. Although ROS acts as a continuous redox rheostat for tumor evolution, it is categorized based on binary models that are masked by tumor-microenvironment (TME) confounders. Here, we have defined a continuous, TME-independent ROS axis to help identify intrinsic vulnerabilities and improve patient stratification. Methods: Non-negative matrix factorization (NMF) defined a ROS-Axis in TCGA-STAD which was validated in ACRG Cohort. Multivariate regression model isolated intrinsic signatures via residual ROS scores by adjusting for TME confounders. Survival was assessed using Cox hazard models. Drug sensitivities were mapped using GDSC2/ElasticNet modeling with cross-cohort replication. Results: Our results define a reproducible ROS gradient, driven by effectors like NQO1 and SOD1, characterizing ROS-high tumors as proliferative, epithelial and immune -cold. High residual ROS score was associated with an improved prognosis, regardless of TNM stage and age. Pharmacogenomic mapping revealed an overlapping sensitivity to mTOR inhibitors in ROS-high gastric cancer tumors which persisted after TME confounder adjustment. Conclusion: The continuous ROS axis provides a functional readout of metabolic dependency that refines traditional anatomical staging. By identifying mTOR dependent cold tumors, our framework offers a precision strategy for immunotherapy-resistant patients like those affected by microsatellite-stable gastric cancer.

9
Prospective Population-Scale Validation of an Electronic Health Record Based Model for Pancreatic Cancer Risk

Lahtinen, E.; Schigiltchoff, N.; Jia, K.; Kundrot, S.; Palchuk, M. B.; Warnick, J.; Chan, L.; Shigiltchoff, N.; Sawhney, M. S.; Rinard, M.; Appelbaum, L.

2026-04-13 oncology 10.64898/2026.04.11.26350318 medRxiv
Top 1%
2.9%
Show abstract

Background and aims: Pancreatic ductal adenocarcinoma (PDAC) surveillance is limited to individuals with familial or genetic risk although most future cases arise outside these groups. In a retrospective study, PRISM, an electronic health record (EHR)-based PDAC risk model, identified individuals in the general population at elevated near-term risk of PDAC. We aimed to prospectively evaluate whether PRISM can identify high-risk individuals beyond current surveillance groups across U.S. health systems. Methods: We performed a prospective multicenter cohort study after deployment of PRISM in April 2023 across 44 U.S. health care organizations. Eligible adults aged [&ge;]40 years without prior PDAC received a single baseline risk score and were assigned to prespecified risk tiers. Patients were followed for incident PDAC for 30 months. We estimated tier-specific 30-month cumulative incidence (positive predictive value, PPV), number needed to screen (NNS), standardized incidence ratios (SIRs), and time from deployment and first high-risk flag to diagnosis. Results: Among 6,282,123 adults assigned a PRISM score, 5,058,067 had follow-up; 3,609 developed PDAC. The highest-risk tier had 30-fold higher PDAC incidence than the study population. At the SIR 5 threshold, 30-month cumulative incidence was 0.35% (NNS, 284.2); at SIR 16, 1.14% (NNS, 87.4); and at SIR 30, 2.19% (NNS, 45.7). Median time from deployment to PDAC diagnosis was 9.5 months, and median time from first high-risk flag to diagnosis at SIR 5 was 3.5 years. Shapley additive explanations (SHAP) analyses supported patient- and tier-level interpretability. Conclusions: Prospective deployment of PRISM across multiple U.S. health care organizations identified individuals at elevated near-term risk for PDAC, with substantial risk enrichment and lead time before diagnosis. These findings support the real-world scalability and generalizability of EHRbased risk stratification for risk-adapted early detection. ClinicalTrials.gov identifier NCT05973331

10
Recombinant zoster vaccination in patients with dementia is associated with improved survival and better cognitive preservation

Soltys, K.; Sara-Buchbut, R.; Ish Shalom, N.; Stokar, J.; Klein, B. Y.; Calderon-Margalit, R.; Greenblatt, C. L.; Ben-Haim, M. S.

2026-04-13 epidemiology 10.64898/2026.04.09.26350509 medRxiv
Top 2%
2.7%
Show abstract

Dementia affects tens of millions of people worldwide, yet disease-modifying treatments remain strikingly limited. Although the recombinant zoster vaccine Shingrix has been associated with reduced dementia incidence, its potential influence on individuals already living with dementia is unknown. Here, we followed a propensity-score matched cohort of 68,960 US dementia patients using a nationwide electronic health record network, comparing Shingrix recipients within two years of diagnosis to recipients of any other vaccine. Shingrix was associated with substantially reduced all-cause mortality across the first three years of follow-up (hazard ratios 0.74, 0.88, and 0.89; P[&le;]0.006), robust across multiple sensitivity analyses. Furthermore, within-individual subgroup analyses of repeated Mini-Mental State Examinations conducted 3-6 years apart revealed significantly divergent cognitive decline rates across groups (time-by-group interaction P=0.002). Interval vaccination was associated with more stable cognition, contrasting with steeper declines in unvaccinated individuals. These findings support prospective evaluation of recombinant zoster vaccination as a potential strategy to improve outcomes in patients with established dementia.

11
Genetic analysis of female genital tract polyps implicates genome stability, estrogen signalling and shared susceptibility with proliferative gynaecological disorders

Ingold, N.; Frankcombe, S.; Bouttle, K.; Moro, E.; Canson, D.; Zoellner, S.; Patil, S.; Dzigurski, J.; Glubb, D. M.; Laisk, T.; O'Mara, T. A.

2026-04-16 genetic and genomic medicine 10.64898/2026.04.13.26350740 medRxiv
Top 2%
2.1%
Show abstract

Female genital tract (FGT) polyps are common benign growths affecting up to half of all women. However, they carry malignant potential, and their genetic architecture remains poorly defined. We conducted a genome-wide association study (GWAS) meta-analysis across four biobanks (48,400 cases, 477,134 controls), identifying 26 risk loci for FGT polyps, 12 of which were previously unreported. Integrative gene prioritisation highlighted 193 candidate genes, revealing a potential convergent biological mechanism: where germline variation in DNA replication and maintenance (e.g., PRIM1, TERT and HMGA1) compromises genomic stability in the context of hormone-driven proliferation (e.g., ESR1 and GREB1). This susceptibility is further modulated by metabolic drivers of estrogen biosynthesis, underscored by specific adiposity-related loci (e.g. RSPO3 and PLCE1) and the aromatase gene CYP19A1. Mendelian randomisation demonstrated bidirectional causal relationships with endometriosis and fibroids, and endometrial cancer. Leveraging the shared genetic architecture of FGT polyps and other gynaecological disorders via multi-trait analysis revealed an additional 26 loci, validating sub-threshold regions encompassing HMGA1 and GREB1. In total, 52 risk loci were identified (36 novel), 39 of which replicated in an independent cohort. These findings reframe polyps not merely as local gynaecological overgrowths but as manifestations of a systemic proliferative syndrome characterised by dysregulated genome stability and estrogen signalling, which may also impact malignant transformation.

12
Independent Genetic Effects of Glucagon-like Peptide-1 Receptor Locus on Body Mass Index and Type 2 Diabetes

Liu, C.; Hui, Q.; Linchangco, G. V.; Dabbs-Brown, A.; Zhou, J. J.; Joseph, J.; Reaven, P. D.; Rhee, M. K.; Djousse, L.; Cho, K.; Gaziano, J. M.; Wilson, P. W.; Phillips, L. S.; The VA Million Veteran Program, ; Sun, Y. V.

2026-04-13 genetic and genomic medicine 10.64898/2026.04.10.26350615 medRxiv
Top 2%
1.8%
Show abstract

Background: The glucagon-like peptide-1 receptor (GLP1R) is a key regulator of glucose metabolism and appetite and a major therapeutic target for type 2 diabetes (T2D) and obesity. Genetic studies have implicated the GLP1R locus in both body mass index (BMI) and T2D, but it remains unclear whether their underlying genetic associations are the same. Methods: We analyzed 431,107 participants of genetically inferred European ancestry from the Million Veteran Program. Within 500 kb of GLP1R, we performed locus-wide linear regression models for BMI and logistic regression models for T2D, adjusted for age, sex, and 10 principal components. We identified primary and secondary BMI sentinel variants using conditional analyses and evaluated their associations with T2D. Bayesian fine-mapping was used to construct credible sets of GLP1R locus for BMI and T2D. Results: Conditioning on the primary sentinel variant rs12213929 (upstream of GLP1R, {beta} = 0.11; 95% CI 0.09-0.14; p = 1.94E-17), we identified a secondary variant (rs13216992, intron of GLP1R) independently associated with BMI ({beta} = 0.10; 95% CI 0.07-0.13; p = 7.88E-14). The two sentinel variants showed low linkage disequilibrium (r2 = 0.03). A two-variant allelic burden score (0-4; sum of the rs12213929 G-allele count and rs13216992 C-allele count) showed that participants with 4 risk alleles had 0.47 kg/m2 higher BMI than those with 0 risk alleles (95% CI 0.39-0.55; p < 2E-16). Both variants were associated with higher T2D risk, but with distinct patterns after BMI adjustment: the rs12213929-T2D association persisted after adjustment for BMI (OR = 1.02; 95% CI 1.01-1.03; p = 0.0004), whereas the rs13216992-T2D association was fully attenuated (OR = 1.00; 95% CI 0.99-1.01; p = 0.68). Fine-mapping identified a compact 95% BMI credible set of 17 variants and a broader 95% T2D credible set of 42 variants, with all BMI credible variants contained within the T2D set. Conclusions: The GLP1R locus harbors at least two independent BMI-associated variants that exhibit heterogeneous relationships with T2D: rs12213929 influences T2D risk partly through BMI-independent pathways, whereas rs13216992 appears to act predominantly via adiposity. These findings refine the genetic architecture at this key therapeutic target gene and provide a foundation for functional and pharmacogenomic studies to determine whether GLP1R variation can inform precision prevention and treatment of obesity and T2D.

13
Functional PD-1/PD-L1 engagement defines a spatial biomarker of immunotherapy response

Ullman, T.; Krantz, D.; Avenel, C.; Lung, M.; Svedman, F. C.; Holmsten, K.; Ostling, P.; Ullen, A.; Stadler, C.

2026-04-17 oncology 10.64898/2026.04.15.26350929 medRxiv
Top 3%
1.7%
Show abstract

Effective predictive biomarkers for immune checkpoint inhibitor (ICI) therapy remain an unmet need across solid tumors. Here, we present an integrated spatial proteomics workflow that combines in situ proximity ligation assay with multiplexed immunofluorescence to directly resolve PD1/PDL1 signaling events at the level of defined cellular phenotypes and their spatial organization within intact tumor tissue. Applied as a proof of concept to tumor samples from patients with metastatic urothelial carcinoma treated with pembrolizumab, this approach reveals that PD1/PDL1 interactions specifically involving cytotoxic CD8CD3 T cells are significantly enriched in complete responders, while such interactions are rare in patients with progressive disease. This interaction defined T cell subset achieves superior discrimination of clinical response compared to single marker PDL1 expression or immune cell abundance alone. By integrating direct detection of protein protein interactions with high dimensional single cell phenotyping, our workflow provides a mechanistically informed, spatially resolved biomarker of functional immune engagement. Beyond urothelial carcinoma, this platform establishes a generalizable framework for translating spatial signaling biology into predictive tools for immunotherapy response across tumor types.

14
Mutation timing, accumulation and selection in the male germline shape inheritance risk for developmental disorders

Neville, M. D. C.; Neuser, S.; Sanghvi, R.; Christopher, J.; Roberts, K.; Smith, K.; ONeill, L.; Hayes, J.; Cagan, A.; Hurles, M. E.; Goriely, A.; Abou Jamra, R.; Rahbari, R.

2026-04-13 genetic and genomic medicine 10.64898/2026.04.09.26350474 medRxiv
Top 3%
1.7%
Show abstract

De novo mutations (DNMs) arising in the parental germline are a major cause of severe developmental disorders. While most DNMs originate in the paternal germline, it remains unclear whether fathers of affected children carry a systematically altered burden of transmissible germline risk, or whether disease largely reflects stochastic outcomes of shared population-wide mutational processes. Here, we combined whole-genome sequencing of 168 parent-child trios with ultra-accurate duplex sequencing of paternal sperm to directly relate transmitted DNMs to the broader mutational and selective landscape of the male germline. In 127 fathers, sperm mutation burden and mutational spectra were indistinguishable from population reference cohorts. Positive selection metrics were likewise concordant, with a global dN/dS of 1.56 (95% CI 1.45-1.67) compared to 1.44 (95% CI 1.17-1.77) in controls and 28 of 32 significantly selected genes overlapping with prior findings. Six fathers harboured a pathogenic early mosaic variant detectable in sperm at allele fractions that ranged from 0.7% to 14.8%. Although these variants generated substantial individual-level risk outliers, they accounted for only [~]11% of the aggregated exome pathogenic burden across the cohort. The remaining burden was distributed across low-VAF mutations, including positively selected driver variants and other rare mutations accumulating with paternal age. Together, these results show that transmissible de novo disease risk is governed primarily by universal germline mutational and selective processes, while early developmental mosaicism produces uncommon but clinically meaningful deviations. This integrated view clarifies how mutation timing, age-associated accumulation and germline selection jointly shape inheritance risk.

15
Inactivating PLEKHA6 Mutations Cause Idiopathic Hypogonadotropic Hypogonadism Through Impaired Kisspeptin Secretion

Topaloglu, A. K.; Plummer, L.; Su, C.-W.; Kotan, L. D.; Celmeli, G.; Simsek, E.; Zhao, Y.; Stamou, M.; Anik, A.; Döger, E.; Altıncık, S. A.; Mengen, E.; Koc, A. F.; Akkus, G.; Balasubramanian, R.; Turan, I.; Seminara, S. B.; Yuksel, B.

2026-04-13 pediatrics 10.64898/2026.04.10.26349358 medRxiv
Top 3%
1.7%
Show abstract

PurposeIdiopathic hypogonadotropic hypogonadism (IHH) is characterized by impaired reproductive maturation, and approximately half of all cases lack an identified genetic cause. We investigated the genetic basis of IHH in two large cohorts to identify novel disease-causing genes. MethodsWe analyzed exome and genome sequencing data from 1,822 patients with IHH from two independent cohorts. Rare variants were filtered using pedigree-informed inheritance models. PLEKHA6 expression in the postmortem human hypothalamus were tested at the mRNA and protein level. Functional studies assessed kisspeptin secretion in cell-based assays. ResultsWe identified 18 distinct PLEKHA6 variants in 24 patients from 20 unrelated families (1.3% of cohort). Variants segregated with disease under autosomal recessive and autosomal dominant (with variable penetrance) inheritance patterns. PLEKHA6 was robustly expressed in the hypothalamus and showed clear colocalization with neurokinin B, which served as the marker for the GnRH pulse generator. Functional studies demonstrated that patient variants significantly impaired kisspeptin secretion. ConclusionPLEKHA6 is a novel IHH gene and the first reported regulator of kisspeptin secretion from the kisspeptin-neurokinin B-dynorphin (KNDy) neurons, which have recently been established as the GnRH pulse generator. These findings establish impaired kisspeptin release as a new disease mechanism in IHH and highlight the critical role of neuropeptide trafficking in reproductive function.

16
Genetic Architecture and Myocardial Fibrotic Remodeling in Mitral Valve Prolapse

Small, A. M.; Yu, M.; Berrandou, T. E.; Georges, A.; Huff, M.; Morningstar, J. E.; Rand, S. A.; Koyama, S.; Lee, J.; Vy, H. M.; Farber-Eger, E.; Jin, S.; Dieterlen, M.-T.; Kontorovich, A. R.; Yang, T.-Y.; Do, R.; Dressen, M.; Krane, M.; Feirer, N.; Doppler, S. A.; Schunkert, H.; Trenkwalder, T.; Wells, Q. S.; Berger, K.; Ostrowski, S. R.; Sorensen, E.; Pedersen, O. B.; Bundgaard, J. S.; Ghouse, J.; Bundgaard, H.; Ganna, A.; Erikstrup, C.; Mikkelsen, C.; Bruun, M. T.; Aagaard, B.; Ullum, H.; Abner, E.; Slaugenhaupt, S. A.; Nadauld, L.; Knowlton, K.; Helgadottir, A.; Sveinbjornsson, G.; Gudbjart

2026-04-13 cardiovascular medicine 10.64898/2026.04.09.26350328 medRxiv
Top 3%
1.6%
Show abstract

Mitral valve prolapse (MVP) is the most common cause of primary mitral regurgitation and is associated with the development of malignant arrhythmias, often in the context of myocardial fibrosis. The genetic architecture of MVP, and whether there are genetic factors explaining why only some individuals with MVP have adverse outcomes, remains poorly understood. We performed a meta-analysis of genome-wide association studies (GWAS) for MVP encompassing 21,517 cases among a total sample size of over 2.2 million individuals. We discovered 89 genomic risk loci for MVP, of which 72 were novel findings. Prioritization of causal genes and pathways using epigenetic and transcriptomic data from mitral valve and extra-valvular tissues replicated known gene associations to MVP including those involved in TGF-{beta} signaling and extracellular matrix biology, but additionally emphasized a role in MVP for biological pathways relevant to cardiomyocyte biology. Accordingly, we identified several MVP risk loci with pleiotropy to cardiomyopathies, especially hypertrophic cardiomyopathy, and demonstrated a significant genetic correlation between MVP and hypertrophic cardiomyopathy. Finally, we interrogated snRNA-seq data in human papillary muscle tissue from two individuals with severe MVP, characterizing genes associated with both risk of papillary muscle fibrosis and MVP.

17
Sex-Stratified Multi-Omics Identifies Sexually Dimorphic Molecular Targets in Parkinsons Disease

Lee, J.-Y.; Lee, J.; Lee, S.; Yoon, J. H.; Park, D. G.; Sung, J.

2026-04-13 genetic and genomic medicine 10.64898/2026.04.10.26350571 medRxiv
Top 3%
1.5%
Show abstract

Parkinsons disease (PD) exhibits well-established sex differences in prevalence and clinical phenotypes, yet the underlying molecular mechanisms remain largely elusive. Here, we conducted a comprehensive sex-stratified multi-omic integration to identify sex-specific causal proteins and biological pathways in PD. We performed gene-based association analysis, transcriptome-wide association studies (TWAS), and proteome-wide Mendelian randomization (PWMR) with colocalization analysis using GWAS summary statistics from the International PD Genetics Consortium (IPDGC; 12,054 male cases/11,999 controls; 7,384 female cases/12,389 controls) for sex-stratified analyses and Global Parkinsons Genetics Program (GP2; 34,933 cases/31,009 controls) for sex-combined analyses. Prioritized candidates were further evaluated through MR with brain expression quantitative trait loci (eQTLs) from MetaBrain and differential protein abundance analysis using the Global Neurodegeneration Proteomics Consortium (GNPC; 704 PD cases/5,629 controls in plasma; 78 cases/1,411 controls in cerebrospinal fluid). Additionally, pathway enrichment analysis was performed for prioritized molecules. Integration across three analytical layers prioritized 102 molecular candidates across 31 unique loci, significant from multiple analyses. Of these, eleven genes reached significance across all three layers, including SNCA, MAPT, and CTSB significant in both sexes; CD160, GPNMB, and LRRC37A2 as male-predominant; STX4 and PRSS53 as female-predominant; and BST1, SCARB2, and LGALS3 significant only in sex-combined analysis. In males, CD160 emerged as a novel candidate with convergent evidence across all three analyses and colocalization, while L3MBTL2 was identified as a novel risk gene from gene-based association and TWAS analyses. In females, STX4 and PRSS53 at the 16p11.2 locus showed female-predominant associations. Pathway enrichment analysis revealed innate immune and SUMOylation pathways in males, with CD160 and L3MBTL2 as key contributors respectively, contrasting with WDR5-mediated chromatin remodeling in females. Brain eQTL-based MR confirmed significant associations for 69 of 86 testable candidates (80.2%) in at least one tissue. Protein abundance analysis confirmed sex-specific patterns, and several candidates showed discordant directions between genetically predicted causal effects and observed protein abundance -- including male-specific plasma elevation of CD160 and female-specific patterns for STX4 -- underscoring the distinction between causal risk mechanisms and disease-state molecular changes. These findings demonstrate that PD is a molecularly heterogeneous disorder with sexually dimorphic pathogenic drivers. While shared axes such as lysosomal dysfunction and vesicle trafficking disruption exist, the divergence into male-specific immune dysregulation and female-specific chromatin remodeling suggests that the primary triggers of neurodegeneration differ by sex. Our results underscore the necessity of sex-stratified approaches in biomarker discovery and the development of precision therapeutic strategies for PD.

18
Characterizing Trainee Workload Using Paired Daily Surveys and EHR Use Data: A Mixed-Methods Pilot Study

Rai, K.; Bianchina, N.; Fischer, C.; Clawson, J.; McBeth, L.; Gottenborg, E.; Keniston, A.; Burden, M.

2026-04-18 health systems and quality improvement 10.64898/2026.04.16.26351062 medRxiv
Top 4%
1.3%
Show abstract

Purpose: High clinical workload is associated with worse patient and hospital outcomes and is a well-established driver of clinician burnout. Trainees may be particularly exposed, shouldering both clinical and educational responsibilities. Evidence-based work design offers a data-driven approach to healthcare work but relies on robust workload measurements. Trainee workload remains poorly characterized, as commonly used metrics (e.g., duty hours, patient census) overlook cognitive and contextual dimensions. This pilot evaluated the feasibility of combining survey-based and electronic health record (EHR) data to characterize internal medicine (IM) trainee workload. Methods: A pilot study was conducted including IM and Medicine-Pediatrics residents (postgraduate years 1-4) between March 31 and June 22, 2025. Participants completed daily surveys during a seven-day inpatient schedule assessing workload and work experience domains, including environment, professional fulfillment, psychological safety, autonomy, and rounding experience, using validated instruments where available. Concurrently, EHR data captured chart review, documentation, orders, and secure messaging activity. Associations between survey and EHR data were assessed. Results: Among 37 eligible residents, 28 (76%) participated in the pilot capturing 166 shifts. Trainees spent 4.4 +/- 1.6 (mean +/- SD) minutes completing daily surveys and 8.6 +/- 2.3 minutes completing the final survey. Trainees reported working 11.6 +/- 1.0 hours/day and a median census of 9.0 (IQR 6.0-11.0). NASA-TLX score was 50.8 +/- 12.6. Positive shift ratings were associated with lower NASA-TLX scores and perceived rounding length. First-to-last EHR login duration was 15 +/- 2 hours/day, and EHR data showed 204 +/- 46 active minutes/day. Login duration correlated with self-reported hours (r=0.43, p<0.0001), and notes signed correlated with self-reported team (r=0.19, p=0.013) and personal census (r=0.34, p<0.0001). Conclusions: Integrating survey-based and EHR-derived workload measures provides multidimensional insight into trainee work. This novel approach supports scalable measurement and evidence-based work design interventions to improve trainee well-being, education, and clinical efficiency.

19
Combined Effects of Severe Immunocompromise and Prolonged Virus Shedding on Within-Host SARS-CoV-2 Evolution in COVID-19

Hirata, Y.; Takahashi, K.; Iwamoto, N.; Dam Jeong, Y.; Miyamoto, S.; Kawasaki, J.; Mine, S.; Iida, S.; Saito, S.; Ainai, A.; Kanno, T.; Katano, H.; Sasaki, N.; Horiba, K.; Ishikane, M.; Kamegai, K.; Harrison, M. T.; Itoh, N.; Akazawa, N.; Okumura, N.; Haraguchi, M.; Sakoh, T.; Morishima, M.; Araoka, H.; Uchida, N.; Hase, R.; Marumo, Y.; Adachi, T.; Matsue, K.; Saito, T.; Ohmagari, N.; Iwami, S.; Suzuki, T.

2026-04-17 infectious diseases 10.64898/2026.04.14.26350918 medRxiv
Top 4%
1.3%
Show abstract

Background: Prolonged SARS-CoV-2 infection in immunocompromised individuals may accelerate virus evolution within the host, raising concerns about the virus evading immunity, developing resistance, and forming novel variants of concern. However, the determinants and public health implications of within-host viral evolution in this population remain incompletely understood. Methods: We performed longitudinal analyses of SARS-CoV-2 genomes from 91 patients with COVID-19 who were classified as being severely or moderately immunocompromised. Using serial measurements of viral RNA loads and infectious titers, we modeled the shedding dynamics of the virus and stratified the infected cases by upper respiratory virus shedding duration to assess associations with within-host evolutionary dynamics. Results: Shedding modeling identified two profiles of shedding duration: intermediate and long. The long shedding profile (shedding lasting >21 days) was found in 14.8% of moderately immunocompromised cases and 72.1% of severely immunocompromised cases. Frequent single-nucleotide variants accumulated specifically in severely immunocompromised individuals with the long shedding phenotype, correlating positively with shedding duration. By contrast, mutations remained limited in moderately immunocompromised individuals with the long shedding phenotype and in severely immunocompromised individuals with the intermediate shedding phenotype. We identified mutations in the spike receptor-binding domain associated with monoclonal antibody resistance; however, we found no fitness-enhancing mutations for inter-host transmission, and antiviral drug resistance mutations were rare. Instead, mutations were introduced frequently and randomly across the entire viral genome. Conclusions: Prolonged upper respiratory virus shedding exceeding 21 days combined with severe immunocompromise is a risk factor of the accumulation of within-host SARS-CoV-2 mutations. Although no variants of concern emerged, the introduction of genome-wide random mutations suggests that the risk for novel variant generation cannot be excluded. These findings highlight the need for intensive antiviral strategies to limit shedding duration to less than 21 days in severely immunocompromised patients, and for immunological investigations to elucidate the host factors underlying residual shedding control in those who achieve clearance within this threshold.

20
A Conversational Artificial Intelligence Framework for Comparative Pathway-Level Profiling of Sezary Syndrome and Primary Cutaneous CD8+ Aggressive Epidermotropic Cytotoxic T-Cell Lymphoma (PCAECTCL)

Diaz, F. C.; Waldrup, B.; Carranza, F. G.; Manjarrez, S.; Velazquez-Villarreal, E.

2026-04-17 oncology 10.64898/2026.04.15.26350992 medRxiv
Top 4%
1.3%
Show abstract

Background: Sezary syndrome (SS) is an aggressive leukemic variant of cutaneous T-cell lymphoma (CTCL) with distinct clinical and biological features compared to rarer entities such as primary cutaneous CD8+ aggressive epidermotropic cytotoxic T-cell lymphoma (PCAECTCL). Although recurrent genomic alterations in CTCL have been described, comparative analyses at the pathway level across biologically divergent subtypes remain limited. Here, we leveraged a conversational artificial intelligence (AI) platform for precision oncology to enable rapid, integrative, and hypothesis-driven interrogation of publicly available genomic datasets. Methods: We conducted a secondary analysis of somatic mutation and clinical data from the Columbia University CTCL cohort accessed via cBioPortal. Cases were stratified into SS (n=26) and PCAECTCL (n=13). High-confidence coding variants were curated and mapped to biologically relevant signaling pathways and functional gene categories implicated in CTCL pathogenesis. Pathway-level mutation frequencies were compared using Chi-square or Fisher's exact tests, with effect sizes quantified as odds ratios. Tumor mutational burden (TMB) was compared using the Wilcoxon rank-sum test. Subtype-specific co-mutation patterns were evaluated using pairwise association analyses and visualized through oncoplots and network heatmaps. Conversational AI agents, AI-HOPE, were used to iteratively refine cohort definitions, prioritize pathway-level signals, and contextualize findings. Results: TMB was comparable between SS and PCAECTCL (p = 0.96), indicating no significant difference in global mutational load. In contrast, pathway-centric analyses revealed marked qualitative differences. SS demonstrated enrichment of alterations in epigenetic regulators, tumor suppressor and cell-cycle control pathways, NFAT signaling, and DNA damage response mechanisms, consistent with transcriptional dysregulation and immune modulation. PCAECTCL exhibited relatively higher frequencies of alterations involving epigenetic regulators and MAPK pathway signaling, suggesting distinct oncogenic dependencies. Co-mutation analysis revealed a more constrained and focused interaction landscape in SS, whereas PCAECTCL displayed broader and more heterogeneous co-mutation networks, indicative of divergent evolutionary trajectories. Notably, ERBB2 mutations were significantly enriched between subtypes (p = 0.031), highlighting a potential subtype-specific therapeutic vulnerability. Conclusions: This study demonstrates that SS is distinguished from PCAECTCL not by increased mutational burden but by distinct pathway-level architectures, particularly involving epigenetic regulation, immune signaling, and transcriptional control. These findings generate biologically grounded, testable hypotheses for subtype-specific therapeutic targeting and underscore the value of conversational AI as a scalable framework for accelerating discovery in translational cancer genomics.