Weekly reads 25/05/26

Mapping biology across space, time, and scale with AI-powered discovery

This week’s reads focus on how increasingly large-scale datasets and sophisticated computational models are transforming our ability to study biology across space, time, and entire organisms. MouseMapper integrates tissue clearing, light-sheet microscopy, and foundation-model based deep learning to map cellular changes across whole mouse bodies and uncovers surprising obesity-linked degeneration of facial sensory nerves. Alternatively, novel spatial and multimodal approaches such as ALARMIST and MultiTME go beyond static snapshots by reconstructing multicellular communication programs and by combining unpaired spatial and single-cell datasets into cohesive representations of tissue biology. Other studies address the challenge of faithfully modelling complex biological systems. Different immune donors humanised patient-derived xenografts retain surprisingly tumor-specific immune ecosystems, opening new avenues to study immune evasion and therapeutic response. Meanwhile, a time-resolved multiomic atlas of breast cancer reveals that tumours and immune cells run on misaligned circadian schedules, creating temporal windows for potential immune escape. Finally, two papers challenge conventional wisdom on biological robustness and sensing. The cancer cells appear to employ layered “recursive” robustness mechanisms, using alternative splicing and paralog compensation to buffer against deleterious mutations, and homing pigeons may use iron-loaded liver macrophages rather than specialised sensory organs to navigate via Earth’s magnetic field under overcast skies.

Preprints/articles that I managed to read this week

A deep-learning framework reveals whole-body perturbations at cell level

Kaltenecker et al. Nature (2026). 10.1038/s41586-026-10535-2

The paper in one sentence

A deep-learning pipeline (MouseMapper) segments nerves, immune cells, and 31 organs across whole cleared mouse bodies, revealing obesity-induced structural degeneration of facial sensory nerves and systemic immune cell redistribution, with molecular signatures conserved in humans.

Summary

Studying systemic diseases such as obesity requires mapping perturbations across multiple organs simultaneously, but tools for whole-body cellular analysis have been lacking. The authors developed MouseMapper, an ensemble of foundation-model-based deep-learning algorithms that segment peripheral nerves, detect immune cells (Cd68+ macrophages), and map 31 organs and tissues across entire cleared mouse bodies imaged by light-sheet fluorescence microscopy. Using transgenic reporter mice (Uchll-eGFP for nerves, Cd68-eGFP for macrophages) fed a high-fat diet for 16–18 weeks, they quantified obesity-induced changes. Unexpectedly, they identified structural degeneration of the infraorbital branch of the trigeminal nerve, reductions in nerve endings, edges, and vertices by ~58–61%, which correlated with functional sensory deficits in whisker stimulation tests. Spatial proteomics of the trigeminal ganglion revealed downregulation of axon-guidance and cytoskeletal pathways (e.g., SERPINA family proteins) and upregulation of complement and inflammation pathways. These molecular changes were recapitulated in post-mortem trigeminal ganglia from obese humans (BMI >30). Additionally, MouseMapper generated whole-body inflammation maps, showing tissue-specific shifts from small to large Cd68+ macrophage clusters in visceral adipose tissue, liver, and other organs, indicating heightened inflammatory states. The framework generalizes across imaging resolutions and antibody-labeling strategies without retraining.

Personal highlights

1. Whole-body segmentation at cellular resolution: MouseMapper segments peripheral nerves over centimetres, detects individual immune cells, and automatically delineates 31 organs/tissues, far exceeding previous methods (e.g., AIMOS segmented only 6 organs), enabling unbiased system-wide screening.

2. Obesity-induced degeneration of facial sensory nerve: Quantitative nerve graph analysis reveals that high-fat diet reduces infraorbital nerve endings by 60.7%, edges by 57.8%, and vertices by 57.6% without changing main nerve trunk thickness, pointing to impaired axonal branching rather than general degeneration.

3. Functional correlate of structural nerve changes: Obese mice exhibit significantly reduced whisker stimulation responses (mean score ~2.5 vs. ~5.5 in lean mice), linking structural trigeminal nerve pathology to sensory dysfunction.

4. Conserved proteomic signatures in mice and humans: Spatial proteomics of trigeminal ganglia identifies downregulation of SERPINA1/3 (anti-inflammatory proteases) and dysregulation of axon guidance, actin cytoskeleton, and complement pathways, changes mirrored in post-mortem human trigeminal ganglia from obese individuals.

5. Tissue-specific inflammation mapping: The framework quantifies shifts in macrophage cluster size distributions across 12 tissues; visceral adipose tissue and liver show a transition from small to medium/large clusters, while subcutaneous fat shows increased large clusters, revealing spatially resolved inflammatory remodeling.

Why should we care?

This study demonstrates the power of combining tissue clearing, light-sheet microscopy, and deep learning to map systemic disease at single-cell resolution across an entire mammalian body. The discovery that obesity damages a specific facial nerve (infraorbital branch of the trigeminal nerve) and impairs whisker sensation in mice, with analogous proteomic changes in obese humans, suggests that obesity may have previously unrecognized effects on sensory function. However, several caveats warrant attention. First, the Uchll-eGFP reporter line labels only a subset of peripheral nerves; the observed structural changes may not represent all nerve types. Second, the functional whisker test measures gross sensory response but does not isolate mechanosensory versus pain pathways. Third, while the proteomic overlap between mouse and human trigeminal ganglia is intriguing, the human samples were from elderly donors (mean age ~85 years), and the observed changes could reflect age-related neuropathy confounded with obesity. Fourth, the computational pipeline requires massive data storage (up to 50 TB per mouse at 4× resolution) and substantial GPU resources, limiting accessibility for many labs.

Humanized patient-derived xenografts preserve tumour-specific immune microenvironments

Stueckmann et al. bioRxiv (2026). 10.64898/2026.05.15.724697

The paper in one sentence

Patient-derived xenografts grown in mice with a reconstituted human immune system (huNOG-EXL) retain key features of the original tumour’s immune cell composition in a tumour-intrinsic, donor-independent manner, though myeloid representation remains incomplete in some models.

Summary

Preclinical cancer models that faithfully reproduce the human tumour immune microenvironment are essential for studying immune evasion and testing immunotherapies. Humanized immune system (HIS) mice, immunodeficient mice engrafted with human CD34+ hematopoietic stem cells, can support patient-derived xenografts (PDXs), but whether these models recapitulate the immune composition of the parental patient tumour has not been systematically evaluated. The authors generated 82 huNOG-EXL mice (expressing human IL-3 and GM-CSF to support myeloid differentiation) using eight different cord blood HSC donors and implanted them with 15 primary tumours (ovarian, head and neck, renal) previously established as PDXs in NSG mice. Using high-dimensional CyTOF and imaging mass cytometry, they profiled immune populations in PDX tumours, spleen, bone marrow, and matched primary tumour samples. Key findings: (1) Immune composition within PDX tumours is primarily driven by the engrafted tumour, not the HSC donor; (2) PDX tumours derived from the same primary sample cluster together across different donor backgrounds; (3) CD8+ T cell and macrophage phenotypic states are tumour-specific and reproducible; (4) Some distinctive features of primary tumours, such as γδ T cell enrichment in an ovarian cancer model and NK cell infiltration in a renal cancer model, are preserved in matched PDX tumours; (5) However, overall immune infiltration tends to be lower in PDX tumours than in primary tumours, and monocyte/macrophage representation is incomplete in certain models (e.g., HNSCC-5, OV-6). The authors also show that a patient-derived cell line xenograft (CDX) from a renal tumour generated an immune composition similar to its matched PDX. The study supports huNOG-EXL PDX models as a platform for studying tumour-intrinsic determinants of immune infiltration, while highlighting remaining limitations in myeloid cell fidelity.

Personal highlights

1. Tumour-intrinsic immune composition across HSC donors: PDX tumours derived from the same primary sample show highly similar immune profiles (median Pearson correlation 0.909) regardless of which human stem cell donor was used, whereas tumours from different primaries are much less similar (0.537). This establishes that the engrafted tumour, not donor variation, is the dominant determinant.

2. Preservation of rare immune features: An ovarian primary tumour (OV-2) with >20% γδ T cells gave rise to six PDX tumours across two HSC donors that all retained this enrichment. Similarly, a renal primary (RCC-3) with abundant NK cells produced eight PDX tumours, six of which showed NK cell infiltrates validated by imaging mass cytometry.

3. Reproducible T cell and macrophage states: CD8+ T cell exhaustion marker expression (PD-1, TIM-3) and macrophage subtype markers (FOLR2, TREM2, PD-L1) were consistent across PDX replicates of the same tumour, indicating that not just cell abundance but also functional states are tumour-driven.

4. Cell line xenografts resemble matched PDXs: A patient-derived cell line from the RCC-2 tumour (RCC-2-CL) generated immune infiltrates (monocytes/macrophages, CD4+ T cells) highly similar to those in five RCC-2 PDX tumours, suggesting that established cell lines may be usable for certain immune profiling studies.

5. Systemic immune composition reflects both HSC donor and tumour: In bone marrow and spleen, samples from mice with the same HSC donor were more similar than those from different donors. However, blood chimerism did not correlate with subsequent tumour immune infiltration, reinforcing that intratumoural immune composition is tumour-determined.

Why should we care?

The development of immunotherapies, drugs that help the immune system kill cancer, has been hampered by a lack of animal models that accurately reflect how human tumours interact with a human immune system. Traditional mouse models either lack a human immune system (so they cannot test human-specific therapies) or use human immune cells but lose the complex, patient-specific features of real tumours. This study systematically tests whether a popular humanized mouse model (huNOG-EXL with PDXs) reproduces the immune cell types found in the original patient tumours. The encouraging finding is that several tumour-specific immune patterns—including unusual enrichments of γδ T cells or NK cells—were preserved across multiple mice, even when the human stem cell donor varied. This suggests that the model can capture biology driven by the tumour itself rather than by individual immune system variation.

Decoding multicellular communication motifs from Spatial Transcriptomics with ALARMIST

Fan et al. bioRxiv (2026). 10.64898/2026.05.21.726986

The paper in one sentence

ALARMIST uses Bayesian Poisson tensor factorization to discover recurring multicellular communication motifs, coordinated sets of cell types and ligand-receptor interactions, from spatial transcriptomics data, revealing higher-order signaling programs missed by pairwise LRI analysis.

Summary

Cell-cell communication in tissues involves multiple cell types sending and receiving multiple signals simultaneously, yet existing computational methods analyze ligand-receptor interactions (LRIs) in isolation, losing the coordinated, higher-order structure of tissue microenvironments. The authors introduce ALARMIST, a probabilistic framework that decomposes a patch-by-LRI count matrix (constructed from spatial neighborhoods) into latent communication motifs. Each motif captures a recurrent pattern: which cell types signal to which others, which LRIs mediate those signals, and how these patterns are spatially organized. The framework uses Bayesian Poisson tensor factorization (BPTF) to handle extreme sparsity, projects motif activities to single-cell resolution, and links motif activation to downstream gene expression changes via Poisson GLMs (excluding ligand/receptor genes to avoid circularity). Benchmarking on semi-synthetic data (generated from scRNA-seq reference with ground-truth motifs) shows ALARMIST recovers cell-type interaction networks with ~0.8 cosine similarity, substantially outperforming COMPOTES, Tensor-cell2cell, and NICHES (<0.6). It also achieves higher F1 scores for motif-associated gene recovery and better spatial reconstruction (Lee’s L) than baseline methods. Cross-platform validation on matched Xenium 5K and CosMx 6K sections (COAD, OV, HCC) shows motif LRI compositions and spatial niches are concordant, with discordance attributable to platform differences in cell-type-specific LRI detection rather than fundamental method instability. Applied to lung adenocarcinoma (LUAD) progression (AIS to invasive), ALARMIST identifies a "healthy vasculature" motif (fibroblast-epithelial support signals, HGF→MET, FGF7→FGFR2) and a "tumor vasculature" motif (VEGFA→KDR, DLL4→NOTCH1). pDCs in healthy vasculature regions upregulate IRF7, TLR7/9, and type I interferon response genes, while tumor vasculature pDCs show downregulation, implicating pDCs as drivers of early inflammatory niches. In glioblastoma (GBM) low-grade to high-grade transformation, ALARMIST identifies an mGAM (malignancy-associated glioma macrophage)-centered hub-and-spoke motif, with GRN→SORT1 signaling from mGAMs to MES-like tumor cells and reciprocal ANXA1→FPR1 signaling.

Personal highlights

1. Higher-order communication motifs, not pairwise LRIs: ALARMIST jointly models all LRI co-occurrence within spatial patches using Bayesian Poisson factorization, capturing recurring multicellular signaling programs (e.g., “healthy vasculature” vs “tumor vasculature”) rather than ranking individual interactions, addressing a fundamental limitation of existing methods.

2. Superior performance on synthetic benchmarks: With ground-truth motifs embedded in scRNA-seq reference data (10 motifs, 16 cell types), ALARMIST recovers cell-type interaction networks with median cosine similarity ~0.8 versus <0.6 for NICHES, <0.5 for Tensor-cell2cell, and <0.4 for COMPOTES. It also achieves higher F1 scores for motif-associated gene recovery and better spatial reconstruction (Lee’s L) across varying cell densities.

3. Cross-platform concordance across Xenium and CosMx: On matched consecutive sections (COAD, OV, HCC), ALARMIST motifs show cosine similarities ~0.7 between platforms, and spatial niche proportions correlate strongly (Pearson r up to ~0.85). Discordance arises from cell-type-level LRI detection differences, not algorithmic instability, providing practical guidance for multi-platform studies.

4. pDCs as drivers of early lung cancer inflammation: In AIS-to-LUAD progression, ALARMIST reveals that “healthy vasculature” motif-active pDCs upregulate IRF7, TLR7/9, and type I interferon response genes, while T cells and macrophages show elevated IFN-γ response and glycolysis. Tumor vasculature motif-active pDCs show downregulated IRF7 and METTL3, with suppressed effector functions. This suggests chronic interferon stimulation may enable immune escape even before tumor invasion, a hypothesis the authors link to ICB resistance mechanisms.

5. mGAM-centered motif and GRN-SORT1 axis in glioma transformation: In low-grade glioma with high-grade regions, ALARMIST identifies a macrophage-centered hub-and-spoke motif where mGAMs (malignancy-associated glioma macrophages) engage MES-like tumor cells via GRN→SORT1 (macrophage→tumor) and receive ANXA1→FPR1 signals (tumor→macrophage). A 20-gene signature from motif-active mGAMs significantly stratifies TCGA LGG patient survival (log-rank p < 0.05), linking a spatial communication program to clinical outcomes.

Why should we care?

ALARMIST addresses a genuine gap: tumors are not just collections of pairwise cell interactions but complex networks where multiple cell types coordinate through multiple signals. By extracting "motifs", repeating patterns of multicellular signaling, the method distills spatial transcriptomics data into interpretable programs that can be mapped across disease stages. This is conceptually powerful. However, several limitations warrant caution. First, ALARMIST assumes mRNA expression correlates with functional protein signaling, ignoring post-transcriptional regulation, secretion dynamics, and receptor internalization. Second, the patch size (default 50 µm) is arbitrary and user-tuned; too large mixes distinct microenvironments, too small yields sparse counts. The authors acknowledge this but provide no data-driven default selection. Third, the binary classification of cells into "motif-active" via GMM thresholding may oversimplify continuous signaling gradients; the authors retain continuous loadings for impact analysis but threshold for spatial maps. Fourth, the cross-platform validation shows non-trivial discordance (e.g., 5–8 CosMx motifs poorly aligned per cancer type), attributed to cell-type LRI detection differences, but this also implies that motif discovery is sensitive to platform-specific biases in gene detection and cell typing, users switching platforms may obtain different biological conclusions. Fifth, the biological findings: pDCs driving early inflammation and GRN-SORT1 in glioma, are intriguing but entirely hypothesis-generating; no experimental validation (e.g., pDC depletion, GRN neutralization) is provided, and the survival signature, while statistically significant, was derived from the same data used to identify the motif (though tested in independent TCGA data).

Cycle-consistent deep generative modeling unifies cellular states across unpaired spatial and single-cell modalities

Zhang et al. bioRxiv (2026). 10.64898/2026.05.25.727736

The paper in one sentence

MultiTME uses cycle-consistent deep generative modeling to integrate unpaired spatial and single-cell datasets (e.g., Xenium, CODEX, scRNA-seq) into a shared latent space, enabling cross-modal cell typing, transcriptome panel completion, and correction of platform-specific technical biases without requiring paired measurements.

Summary

The authors present MultiTME, a variational autoencoder framework that learns a shared latent representation across unpaired modalities. Key innovations include: (1) modality-specific projection layers that map heterogeneous feature spaces (e.g., 50 proteins vs 5,000 genes) to a common intermediate dimension; (2) a shared encoder with no modality identifiers, forcing the latent space to capture biological state rather than platform artifacts; (3) cycle consistency losses that enforce bidirectional translation between modalities at both latent and observation levels, aligning distributions without paired data; (4) a spatial regularizer that encourages cells of the same type in local tissue neighborhoods to have similar latent representations; and (5) optional semi-supervision using marker genes or expert annotations. Benchmarking on a human tonsil dataset (scRNA-seq + CODEX) shows MultiTME achieves 94.7% cell typing accuracy, outperforming MaxFuse (89.1%), Celesta (70.9%), and Astir (73.7%). On colorectal cancer data (scRNA-seq + Xenium), MultiTME imputes held-out genes with median Pearson correlation ~0.5, significantly better than ENVI, Harmony, and StabMap (<0.3). The imputed transcriptome reveals a spatially organized proliferative–invasive tumor axis not visible from Xenium alone. On serial sections of Visium HD (spatially resolved bulk transcriptomics) and CODEX (single-cell proteomics), MultiTME assigns whole-transcriptome profiles to individual CODEX cells, achieving per-gene correlations ~0.7 with ground-truth Visium HD, compared to ~0.1 for iStar. Finally, across five cancer types, MultiTME translates CosMx measurements to match Xenium, correcting platform-specific background biases and improving cross-platform concordance (R² increases from ~0.42 to near 1.0 after translation).

Personal highlights

1. Cycle consistency for unpaired multimodal integration: Enforces bidirectional translation between modalities without requiring paired cells or shared features. This allows MultiTME to align dissociated scRNA-seq with spatially resolved CODEX or Xenium data, where no cell-level correspondence exists.

2. Spatially regularized latent space: Uses k-nearest neighbors in physical tissue coordinates weighted by cell-type probabilities to encourage that neighboring cells of the same type have similar latent representations. This preserves spatial organization in the integrated embedding, a feature absent from most single-cell integration methods.

3. State-of-the-art cross-modal cell typing: On tonsil scRNA-seq + CODEX, MultiTME achieves 94.7% accuracy in transferring scRNA-seq annotations to spatial proteomic cells, substantially outperforming MaxFuse (89.1%) and proteomic-only classifiers (70–74%). The confusion matrix shows reduced misclassification among closely related lymphocyte subsets (CD4 vs CD8 T cells, germinal-center vs proliferating B cells).

4. Whole-transcriptome panel completion and spatial super-resolution: MultiTME imputes missing genes in Xenium panels (median per-gene Pearson r ~0.5) and assigns full transcriptomes to CODEX cells from adjacent Visium HD sections, recovering spatial expression patterns of EPCAM and MUC2 at single-cell resolution, sharper than Visium HD bins and far better than the iStar baseline.

5. Platform bias correction across CosMx and Xenium: On five cancer types, MultiTME translation from CosMx to Xenium dramatically improves cross-platform gene expression concordance (slope from 0.42 to near 1.0, R² increases substantially). The model generalizes to held-out fields of view and, to a lesser extent, to unseen cancer types, suggesting it learns transferable technical correction.

Why should we care?

MultiTME addresses a conceptually big problem: biological measurements are always incomplete and biased by the technology used. A cell’s transcriptome measured by Xenium differs systematically from the same cell measured by CosMx; a protein panel captures only a fraction of functional state; scRNA-seq loses all spatial context. MultiTME’s core assumption, that there exists a shared underlying biological state that generates all these modality-specific observations, is reasonable but unprovable. The cycle consistency mechanism is elegant: if translating from modality A to B and back returns the original, then the translation must preserve biologically relevant information. However, this does not guarantee that the latent space captures true biology rather than a technically convenient but biologically arbitrary alignment. Several limitations temper enthusiasm. First, the model requires either expert annotations or high-confidence marker-based pseudo-labels to anchor cell types. In datasets where marker specificity is poor or cell states are not well captured by known markers, the semi-supervised regularization may propagate errors. Second, the spatial regularizer assumes that consecutive tissue sections preserve local cell-type composition, which holds approximately but ignores tissue deformation, cell migration, and differences in sectioning plane. Third, the panel completion results, while statistically significant, show median correlations around 0.5, meaningful for pathway-level analyses but insufficient for confident single-gene conclusions, especially for low-abundance transcripts. The authors demonstrate that aggregated pathway scores are better preserved than individual genes, which is the appropriate use case. Fourth, the platform bias correction is impressive but the leave-one-disease-out generalization drops noticeably, indicating that some platform-specific biases are tissue-dependent; a universal “CosMx-to-Xenium” transformer does not yet exist

Circadian misalignment underlies immune escape in breast cancer

Liang et al. bioRxiv (2026). 10.64898/2026.05.26.726543

The paper in one sentence

Time-resolved single-nucleus multiomic mapping of the breast cancer tumor microenvironment reveals that cancer epithelial cells and immune cells are rhythmically but desynchronized across the circadian cycle, creating temporal gaps in antigen presentation, T cell activation, and checkpoint signaling that promote immune evasion.

Summary

Circadian rhythms regulate cellular processes, but how they are coordinated across cell types within the tumor microenvironment (TME) has remained largely unexplored. The authors performed snRNA-seq and snATAC-seq on 4T1 mouse triple-negative breast cancer (TNBC) tumors collected at four circadian time points (CT4, CT10, CT16, CT22), generating a time-resolved multiomic atlas of 101,321 nuclei spanning 14 cell types. They identified 6,224 circadian genes across the TME, organized into 20 functional modules covering cell cycle, immune activation, metabolism, and extracellular matrix remodeling. Cancer epithelial cells exhibited a striking temporal switch: a proliferative state (high cell cycle, DNA repair) peaked around the early daytime (CT0–3), while an inflammatory state (high antigen presentation, interferon response) peaked at night. In contrast, CD4⁺ and CD8⁺ T cells and macrophages showed peak activation and effector programs during the nighttime (murine active phase). This created three axes of circadian misalignment: (1) tumor proliferation decoupled from immune activation; (2) antigen presentation (peaking at night) out of phase with T cell recognition (peaking during the day); and (3) PD-1 expression in T cells (peak at night) asynchronous with PD-L1 expression in cancer cells and macrophages (peak during the day). Within T cells, activation and exhaustion programs overlapped temporally, potentially promoting dysfunction. Using CYCLOPS 2.0 to infer circadian phase from human TNBC scRNA-seq data, the authors found conserved phase relationships and similar misalignment patterns, with immune activation peaking during the inferred daytime (active phase for diurnal humans). The study proposes that circadian desynchronization across TME compartments is a previously underappreciated mechanism of tumor immune evasion and suggests that timing of immunotherapy may influence efficacy.

Personal highlights

1. Comprehensive circadian atlas of the TME: Time-resolved single-nucleus multiomics across 14 cell types reveals 6,224 rhythmic genes organized into 20 functional modules, including cell-type-specific and shared circadian programs in cancer cells, immune populations, and stromal cells.

2. Cancer epithelial cells oscillate between proliferative and inflammatory states: A time-of-day–dependent shift from a proliferative state (peaking at early daytime, CT0–3) to an inflammatory state (peaking at night, CT20–24), characterized by differential expression of cell cycle vs. antigen presentation and interferon-response genes. Both state composition and intrinsic subtype-specific rhythms contribute to this oscillation.

3. Three axes of circadian misalignment: (i) Tumor proliferation peaks opposite to T cell/macrophage activation; (ii) Antigen presentation (cancer and myeloid cells) peaks at night, while T cell recognition peaks during the day; (iii) PD-1 expression in T cells peaks at night, whereas PD-L1 in cancer/myeloid cells peaks during the day, creating prolonged checkpoint signaling across the cycle.

4. Temporal overlap of T cell activation and exhaustion: In CD4⁺ T cells, activation/costimulation genes peak alongside exhaustion markers (Ctla4, Tox) during the daytime, suggesting that circadian-driven activation may simultaneously engage inhibitory circuits, constraining effective antitumor responses.

5. Conserved circadian architecture in human TNBC: Using CYCLOPS 2.0 to infer phase from single-cell data, the authors demonstrate preserved phase relationships among immune populations and similar misalignment between antigen presentation and T cell recognition, supporting translational relevance for chronotherapy.

Why should we care?

The main takeaway of this study is that the time of day matters for tumor–immune interactions, and this has potential implications for when cancer immunotherapies are administered. But the field is still early: robust prospective clinical trials testing timed immunotherapy are needed before any practice change. The paper is a strong hypothesis-generating resource, not a therapeutic guideline. It also highlights a broader principle: biological systems are not static snapshots; time is a dimension that must be integrated into our understanding of disease mechanisms and treatment design.

Recursive mutational robustness in cancer through intra- and inter-genic compensation

Dandage et al. bioRxiv (2026). 10.64898/2026.05.26.727768

The paper in one sentence

Cancer cells tolerate deleterious mutations by upregulating alternatively spliced isoforms that skip the mutated region, a form of intra-genic compensation that often works together with paralog-mediated inter-genic compensation in a recursive manner, driven by nonsense-mediated decay and transcriptional adaptation.

Summary

Cancer cells harbor hundreds to thousands of somatic mutations, yet most are tolerated without catastrophic fitness loss. Known mechanisms include paralog buffering (inter-genic compensation), but many essential genes lack paralogs. This study explores whether alternative splicing provides intra-genic functional redundancy, where one isoform can compensate for another that carries a deleterious mutation. Using pan-cancer genomics and transcriptomics data (30 tumor types from TCGA, cancer cell lines from CCLE), the authors systematically identified cases where mutations occur in exons that can be skipped in alternative isoforms. They found that mutation-skipping isoforms are frequently expressed (median ~33% of isoforms per gene) and often upregulated in samples carrying the mutation compared to matched controls. They developed a "mutational robustness score" combining upregulation magnitude and compensation extent. Strong intra-genic robustness was associated with higher perturbation tolerance (higher deleterious mutation frequency) and was more context-specific than inter-genic robustness. Notably, genes exhibiting intra-genic robustness also showed stronger inter-genic compensation through paralogs, a "recursive" architecture. Mechanistically, they implicate nonsense-mediated mRNA decay (NMD) of mutation-bearing isoforms triggering self-transcriptional adaptation (self-TA), leading to upregulation of the gene's pre-mRNA and consequently relative increase of mutation-skipping isoforms. The study also shows that mutation-skipping isoform upregulation correlates with differential isoform usage in protein interaction partners, suggesting restorative rewiring of protein complexes. Tumor suppressor genes (TSGs) showed significantly lower robustness than non-TSGs, and high robustness in non-TSGs was associated with worse patient survival across most cancer types, indicating that cancer cells exploit this buffering to preserve pro-tumorigenic functions.

Personal highlights

1. Widespread expression of mutation-skipping isoforms: Across 30 cancer types, ~33% of a gene’s isoforms that skip deleterious mutations are expressed, comparable to the overall fraction of expressed isoforms (~39%). Exon skipping is the predominant splicing event enabling mutation avoidance.

2. Compensatory upregulation of skipping isoforms: In tumors carrying mutations that perturb specific isoforms, the mutation-skipping isoforms are frequently upregulated (162 significant events at FDR<0.1, far outnumbering downregulation). This trend holds in CRISPR perturbation data, where perturbation-skipping isoforms increase expression after sgRNA targeting.

3. Recursive robustness between isoform and paralog levels: Genes with strong intra-genic robustness (score >0.5) show significantly higher inter-genic robustness scores (p=3e-5), and vice versa. This suggests layered compensatory architectures where alternative splicing and gene duplication provide nested buffering against mutational insults.

4. Self-transcriptional adaptation as mechanism: Mutations that generate premature termination codons (frameshift, stop-gain) are most associated with intra-genic robustness. NMD of these isoforms correlates with increased pre-mRNA abundance (rs=0.29, p<1e-10), supporting a model where degradation intermediates trigger transcriptional upregulation of the same gene, selectively enriching mutation-skipping isoforms.

5. Tumor suppressors evade robustness, non-TSGs exploit it: TSGs have significantly lower intra-genic and inter-genic robustness scores compared to non-TSGs, suggesting their inactivation requires not just mutation but also circumvention of compensatory mechanisms. Conversely, high robustness in non-TSGs correlates with worse survival in most cancer types (hazard ratio >1 in 16/20 cancer types), indicating cancer cells depend on protecting non-essential but pro-growth functions.

Why should we care?

The main takeaway is that cancer's mutational tolerance is not merely passive, it involves active transcriptional adaptation through alternative splicing and paralog upregulation. This opens potential therapeutic angles: if cancer cells depend on specific mutation-skipping isoforms for survival, those isoforms could be targeted (e.g., with antisense oligonucleotides) to create synthetic lethality. However, the translational gap is large; the study is hypothesis-generating and would require extensive functional validation before any clinical application. As a conceptual advance, it reframes how we think about genetic robustness—not just as duplication-driven but as a layered, recursive property of gene architecture. But the mechanistic details remain speculative, and the clinical relevance awaits prospective testing.

Homing pigeon navigation relies on superparamagnetic macrophages under overcast conditions

Lisowski et al. Science (2026). 10.1126/science.ady2486

The paper in one sentence

Superparamagnetic iron-accumulating macrophages in the pigeon liver are required for magnetic orientation when solar cues are unavailable, suggesting a novel organ-level mechanism for magnetoreception.

Summary

For decades, the mechanisms by which birds sense Earth’s magnetic field have remained controversial, with competing hypotheses implicating cryptochromes in the retina, magnetite particles in the beak, or ion-channel perturbations in the vestibular system. This study reports a fourth, unexpected mechanism: macrophages in the pigeon liver accumulate ferric iron (Fe³⁺) within ferritin protein nanocages, rendering them superparamagnetic—a property previously described in mammalian splenic red pulp macrophages. Using vibrating sample magnetometry, the authors show that pigeon liver (and to a lesser extent spleen) exhibits a magnetic blocking temperature (T_B) of 12 K and hysteresis loops at low temperatures, consistent with superparamagnetic ferritin nanoparticles, whereas muscle, beak, and eye lack such signals. Histological staining (Prussian blue) reveals iron-positive cells exclusively in liver and spleen, colocalizing with MHC II⁺ macrophages. Magnetic column separation and single-cell RNA sequencing confirm these cells express macrophage signature genes (including Spi-c, involved in erythrocyte clearance) and phagocytose dextran. Clodronate liposome treatment eliminates magnetic liver cells, reduces Spi-c expression, and depletes MHC II⁺ macrophages, without affecting heterophils (avian neutrophils). Electron microscopy shows that iron-laden macrophages reside within 2 μm of unmyelinated nerve fibers in the hepatic portal triad, suggesting potential neuro-immune signaling. Crucially, the authors performed a behavioral experiment: 34 homing pigeons trained over a 19 km route were randomly assigned to clodronate or control liposomes. Under completely overcast conditions (no sun or polarized light cues), all 16 control pigeons homed within 70 minutes, whereas none of the 18 clodronate-treated birds returned that day, showing random spatial orientation. When cloud cover cleared and the sun became visible, clodronate-treated pigeons homed normally, indicating intact flight ability and reliance on solar cues when available. The authors propose that superparamagnetic hepatic macrophages collectively sense the geomagnetic field and transmit directional information to the brain via afferent vagal or sympathetic innervation.

Personal highlights

1. Superparamagnetic macrophages in the liver: Vibrating sample magnetometry identifies a magnetic blocking temperature (T_B = 12 K) and hysteresis loops in pigeon liver and spleen, but not in muscle, beak, or eye. Prussian blue staining confirms ferric iron accumulation exclusively in these organs, localized to cells with macrophage morphology and MHC II expression.

2. Macrophage depletion abolishes magnetic orientation: Under overcast skies (no solar cues), 0/18 clodronate-treated pigeons homed vs 16/16 controls. The same depleted birds homed normally when the sun emerged, proving the effect is specific to magnetic orientation, not general flight impairment.

3. Macrophage–nerve proximity: Electron microscopy reveals iron-positive macrophages within 2 μm of unmyelinated nerve bundles in the liver portal triad, providing an anatomical substrate for signal transmission to the brain via the autonomic nervous system.

4. Mechanistic link to iron metabolism: Pigeon liver macrophages express Spi-c and other genes involved in erythrocyte clearance and ferritin storage, similar to red pulp macrophages in mammals. The superparamagnetic property arises naturally from the ferritin nanocages that sequester iron from hemoglobin degradation.

5. A fourth magnetoreception mechanism: Distinct from cryptochrome-based (light-dependent), beak magnetite, and vestibular ion-channel hypotheses, this liver macrophage-based mechanism operates under overcast conditions and may generalize to other animals (e.g., bats, sharks) that navigate without visual cues.

Why should we care?

This paper challenges a long-standing orthodoxy in sensory biology: that magnetic field detection resides in specialized cells within the head (eye, beak, inner ear). Instead, it implicates a peripheral organ, the liver, and a cell type better known for immune defense than for sensing physical forces. The finding that iron-accumulating macrophages, which exist in many vertebrates, can become superparamagnetic and potentially influence navigation is conceptually striking. The behavioral experiment is clean: clodronate treatment removes the cells, and pigeons lose their ability to orient under overcast conditions but not under sunny skies, ruling out non-specific toxicity.

Other papers that peeked my interest and were added to the purgatory of my “to read” pile

• Scalable subclonal reconstruction of cancer cells in DNA sequencing data using a penalized likelihood model

• A targetable opioid/cancer associated fibroblast axis drives extracellular matrix remodeling and tumor aggressiveness in pancreatic cancer

• Mapping Tumor-Microenvironment dependencies with TMEformer: A spatial foundation framework enabling in silico perturbation

• Tumor-derived branched-chain α-keto acids activate Notch signaling in tumor-associated macrophages to limit immunity

• Decoding sequence determinants of gene expression in diverse cellular and disease states

• Multiplexed perturbation enables scalable pooled screens

• Spatially resolved, multimodal in vivo Perturb-seq using antibody-based cell hashing

• MetaboRamics: Highly multiplexed metabolic imaging by stimulated Raman for spatial metabolomics in live cells

• Many Needles in a Haystack: Active Hit Discovery for Perturbation Experiments

• What Makes a Representation Good for Single-Cell Perturbation Prediction?

• SCALLOPS: a scalable, integrated computational framework for Optical Pooled Screens

• Scoring gene importance by interpreting single-cell foundation models

• Spatiotemporal transcriptome atlas of human embryos after gastrulation

• Geometry-First Generative Spatial Single-Cell Reconstruction

• Massively multiplex multimodal chemical screens at single-cell resolution

• gRely: Relyability for genome trained sequence-to-expression models

• Unbiased niche labeling maps immune-excluded niche in bone metastasis

• BioArc: Discovering Optimal Neural Architectures for Biological Foundation Models

• Pan-cancer spatial atlas of tertiary lymphoid structures

• Language Modeling Materializes a World Model of Protein Biology

• Cophenetic Spatial Topology Embedding reveals multiscale tissue architecture in spatial omics

• IL-1-activated cancer-associated fibroblasts promote STAT1-driven transcriptional reprogramming of pancreatic tumour cells

• Fibroblast-derived thrombospondin-1 shapes macrophage polarization in advanced human co-culture models

Thanks for reading.

Cheers,

Seb.