Spatial Expression Research Articles

ResSAT: Enhancing Spatial Transcriptomics Prediction from H&E- Stained Histology Images with Interactive Spot Transformer.

Spatial transcriptomics (ST) revolutionizes RNA quantification with high spatial resolution. Hematoxylin and eosin (H&E) images, the gold standard in medical diagnosis, offer insights into tissue structure, correlating with gene expression patterns. Current methods for predicting spatial gene expression from H&E images often overlook spatial relationships. We introduce ResSAT (Residual networks - Self-Attention Transformer), a framework generating spatially resolved gene expression profiles from H&E images by capturing tissue structures and using a self-attention transformer to enhance prediction.Benchmarking on 10× Visium datasets, ResSAT significantly outperformed existing methods, promising reduced ST profiling costs and rapid acquisition of numerous profiles.

Combinatorial regulatory states define cell fate diversity during embryogenesis

Cell fate specification occurs along invariant species-specific trajectories that define the animal body plan. This process is controlled by gene regulatory networks that regulate the expression of the limited set of transcription factors encoded in animal genomes. Here we globally assess the spatial expression of ~90% of expressed transcription factors during sea urchin development from embryo to larva to determine the activity of gene regulatory networks and their regulatory states during cell fate specification. We show that >200 embryonically expressed transcription factors together define >70 cell fates that recapitulate the morphological and functional organization of this organism. Most cell fate-specific regulatory states consist of ~15–40 transcription factors with similarity particularly among functionally related cell types regardless of developmental origin. Temporally, regulatory states change continuously during development, indicating that progressive changes in regulatory circuit activity determine cell fate specification. We conclude that the combinatorial expression of transcription factors provides molecular definitions that suffice for the unique specification of cell states in time and space during embryogenesis.

A spatial expression atlas of the adult human proximal small intestine.

The mouse small intestine shows profound variability in gene expression along the crypt-villus axis1,2. Whether similar spatial heterogeneity exists in the adult human gut remains unclear. Here we use spatial transcriptomics, spatial proteomics and single-molecule fluorescence in situ hybridization to reconstruct a comprehensive spatial expression atlas of the adult human proximal small intestine. We describe zonated expression and cell type representation for epithelial, mesenchymal and immune cell types. We find that migrating enterocytes switch from lipid droplet assembly and iron uptake at the villus bottom to chylomicron biosynthesis and iron release at the tip. Villus tip cells are pro-immunogenic, recruiting γδ T cells and macrophages to the tip, in contrast to their immunosuppressive roles in mouse. We also show that the human small intestine contains abundant serrated and branched villi that are enriched at the tops of circular folds. Our study presents a detailed resource for understanding the biology of the adult human small intestine.

18 Unraveling the Impact of Tumor Microenvironment on Immunotherapy Response/Resistance in Human Sarcomatoid Clear Cell Renal Cell Carcinoma: Insights from a Novel Mouse Model

Abstract Background Renal cell carcinoma (RCC) with histological sarcomatoid de-differentiation (sRCC) has a historically poor prognosis across all RCC subtypes. Recently, immunotherapy with anti-PD1 + anti-CTLA4 combo has significantly improved disease-specific survival, and complete response are seen in upwards of 20% of patients. However, half of sRCC patients still have no response to this immunotherapy, and the mechanism of response remains unclear. Methods Our preliminary bulk RNA sequencing data reveal a distinct tumor microenvironment (TME) in human sRCC tumors, characterized by heightened immune cell infiltration, particularly enriched for regulatory T cells (Tregs). Subsequently, we delineated the spatial gene expression topography within sRCC and clear cell renal cell carcinoma (ccRCC) regions, culminating in the development of the sRCC signature. Concurrently, we established an innovative immunocompetent murine model of sRCC. Leveraging CRISPR technology, we replicated the human sRCC genotype by knocking out Vhl, Bap1, and Cdkn2a/2b genes in a mouse telomerase reverse transcriptase (mTERT) overexpressing renal proximal tubule epithelial cell line. This effort yielded a successfully engineered mouse implantable sRCC cell line bearing VhlMutBap1DelCdkn2aDelCdkn2bDel (HJRCC68N). Results The syngeneic mouse sRCC model we developed faithfully replicates human sRCC histology and tumor microenvironment and recapitulates the sRCC signature identified in human sRCC regions. Importantly, like in humans, mouse sRCC tumors exhibit a heterogeneous response to anti-PD1 + anti-CTLA4 combination, and anti-CTLA4 monotherapy, reflecting the clinical variability observed in human patients. Furthermore, transcriptomic analysis from our mouse model reveals the involvement of type 1 immunity in responder, providing insights into overcoming immunotherapy resistance. Conclusions In conclusion, our study unveiled distinct changes in the TME at the transcriptomic and spatial gene expression levels in human sRCC. Leveraging our novel murine model, which faithfully replicates human sRCC histology, TME characteristics, and mixed response to immunotherapy, we elucidated the involvement of type 1 immunity in treatment response. These insights offer valuable pathways for overcoming immunotherapy resistance and pave the way for the development of tailored treatment strategies aimed at improving outcomes for sRCC patients.

Single-cell spatial transcriptomics reveals distinct patterns of dysregulation in non-neuronal and neuronal cells induced by the Trem2R47H Alzheimer's risk gene mutation.

The R47H missense mutation of the TREM2 gene is a known risk factor for development of Alzheimer's Disease. In this study, we analyze the impact of the Trem2R47H mutation on specific cell types in multiple cortical and subcortical brain regions in the context of wild-type and 5xFAD mouse background. We profile 19 mouse brain sections consisting of wild-type, Trem2R47H, 5xFAD and Trem2R47H; 5xFAD genotypes using MERFISH spatial transcriptomics, a technique that enables subcellular profiling of spatial gene expression. Spatial transcriptomics and neuropathology data are analyzed using our custom pipeline to identify plaque and Trem2R47H-induced transcriptomic dysregulation. We initially analyzecell type-specific transcriptomic alterations induced by plaque proximity. Next, we analyze spatial distributions of disease associated microglia and astrocytes, and how they vary between 5xFAD and Trem2R47H; 5xFAD mouse models. Finally, we analyze the impact of the Trem2R47H mutationon neuronal transcriptomes. The Trem2R47H mutation induces consistent upregulation of Bdnf and Ntrk2 across many cortical excitatory neuron types, independent of amyloid pathology. Spatial investigation of genotype enriched subclusters identified spatially localized neuronal subpopulations reduced in 5xFAD and Trem2R47H; 5xFAD mice. Overall, our MERFISH spatial transcriptomics analysis identifies glial and neuronal transcriptomic alterations induced independently by 5xFAD and Trem2R47H mutations, impacting inflammatory responses in microglia and astrocytes, and activity and BDNF signaling in neurons.

Characterization of c-Jun N-terminal kinase (JNK) gene reveals involvement of immune defense against Vibrio splendidus infection in Apostichopus japonicus

The c-Jun N-terminal kinase (JNK) constitutes an evolutionarily conserved family of serine/threonine protein kinases, pivotal in regulating various physiological processes in vertebrates, encompassing apoptosis and antibacterial immunity. Nevertheless, the involvement of JNK in the innate immune response remains largely unexplored in pathogen-induced echinoderms. We isolated and characterized the JNK gene from Apostichopus japonicus (AjJNK) in our investigation. The full-length cDNA sequences of AjJNK spanned 1806 bp, comprising a 1299 bp open reading frame (ORF) encoding 432 amino acids, a 274 bp 5′-untranslated region (UTR), and a 233 bp 3′-UTR. Structural analysis revealed the presence of a classical S_TKc domain (37–335 amino acids) within AjJNK and contains several putative immune-related transcription factor-binding sites, including Elk-1, NF-κB, AP-1, and STAT5. Spatial expression analysis indicated ubiquitous expression of AjJNK across all examined tissues, with the highest expression noted in coelomocytes. The mRNA, protein, and phosphorylation levels of AjJNK were obviously induced in coelomocytes upon V. splendidus challenge and lipopolysaccharide stimulation. Immunofluorescence analysis demonstrated predominant cytoplasmic localization of AjJNK in coelomocytes with subsequent nuclear translocation following the V. splendidus challenge in vivo. Moreover, siRNA-mediated knockdown of AjJNK led to a significant increase in intracellular bacterial load, as well as elevated levels of Ajcaspase 3 and coelomocyte apoptosis post V. splendidus infection. Furthermore, the phosphorylation levels of AjJNK inhibited by its specific inhibitor SP600125 and also significantly suppressed the expression of Ajcaspase 3 and coelomocyte apoptosis during pathogen infection. Collectively, these data underscored the pivotal role of AjJNK in immune defense, specifically in the regulation of coelomocyte apoptosis in V. splendidus-challenged A. japonicus.

Mutual exclusivity and co-occurrence patterns of immune checkpoints indicate NKG2A relates to anti-PD-1 resistance in gastric cancer

BackgroundAn increasing number of clinical studies have begun to explore combination strategies with immune checkpoint inhibitors, aiming to present new opportunities for overcoming anti-PD-1 treatment resistance in gastric cancer. Unfortunately, the exploration of certain immune checkpoint inhibitor combination strategies has yielded suboptimal results. Therefore, it is necessary to comprehensively analyze the expression patterns of immune checkpoints and identify optimal combination regimens of anti-PD-1 inhibitors with other immune checkpoint inhibitors.MethodsLeveraging single-cell RNA sequencing (scRNA-seq) and multivariate linear regression interaction models, we dissected the immune checkpoint expression characteristics of CD8+ T cells in gastric cancer and the immune checkpoint expression pattern (ICEP) mediating anti-PD-1 treatment resistance. Furthermore, we employed transcription factor analysis and CellOracle to explore the transcriptional regulatory mechanisms governing CD8+ T cell differentiation fates. Finally, we utilized Nichenet and spatial transcriptomic analysis to investigate the spatial expression patterns of immune checkpoints.ResultsInteraction analysis indicated that, among the known immune checkpoints, co-expression of NKG2A and PD-1 might exert a more profound inhibitory effect on the proliferative capacity of CD8+ T cells. The co-expression analysis revealed differential co-expression pattern of PD-1 and NKG2A, defined as ICEP1 (CD8+ T cells co-expressing PD-1, CTLA-4, TIGIT, LAG-3 or CD38) and ICEP2 (CD8+ T cells solely expressing NKG2A or co-expressing with other immune checkpoints), reflecting the co-occurrence pattern of PD-1 and the mutual exclusivity of NKG2A. Further, these two ICEP CD8+ T cell subsets represented distinct CD8+ T cell differentiation fates governed by MSC and RUNX3. Notably, ICEP2 CD8+ T cells were associated with anti-PD-1 therapy resistance in gastric cancer. This phenomenon may be attributed to the recruitment of LGMN+ macrophages mediated by the CXCL16-CXCR6 signaling pathway.ConclusionThis study unveiled two distinct ICEPs and the mutually exclusivity and co-occurrence characteristics of CD8+ T cells in gastric cancer. The ICEP2 CD8+ T cell subset, highly expressed in gastric cancer patients resistant to anti-PD-1 therapy, may be recruited by LGMN+ macrophages through CXCL16-CXCR6 axis. These findings provide evidence for NKG2A as a novel immunotherapeutic target in gastric cancer and offer new insights into combination strategies for immune checkpoint inhibitors in gastric cancer.Graphical

Spatial transcriptomics reveals organized and distinct immune activation in cutaneous granulomatous disorders

BackgroundNon-infectious (inflammatory) cutaneous granulomatous disorders include cutaneous sarcoidosis (CS), granuloma annulare (GA), necrobiosis lipoidica (NL), and necrobiotic xanthogranuloma (NXG). These disorders share macrophage predominant inflammation histologically, but the inflammatory architecture and the pattern of extracellular matrix alteration varies. The underlying molecular explanations for these differences remain unclear. ObjectiveTo understand spatial gene expression characteristics in these disorders. MethodsWe performed spatial transcriptomics in cases of CS, GA, NL, and NXG to compare patterns of immune activation and other molecular features in a spatially resolved fashion. ResultsCS is characterized by a polarized, spatially organized T helper (Th) 1 predominant response with classical macrophage activation. GA is characterized by a mixed, but spatially organized pattern of Th1 and Th2 polarization with both classical and alternative macrophage activation. NL showed concomitant activation of Th1, Th2, and Th17 immunity with a mixed pattern of macrophage activation. Activation of type 1 immunity was shared among, CS, GA, and NL and included upregulation of IL-32. NXG showed upregulation of CXCR4-CXCL12/14 chemokine signaling and exaggerated alternative macrophage polarization. Histologic alteration of extracellular matrix correlated with hypoxia and glycolysis programs and type 2 immune activation. ConclusionsInflammatory cutaneous granulomatous disorders show distinct and spatially organized immune activation that correlate with hallmark histologic changes.

Cellular changes and PR-1 gene expression accompanying infection-related morphogenesis and cell progression during invasive growth Pyricularia oryzae pathotype Triticum in resistant and susceptible wheat cultivars

Pyricularia oryzae pathotype Triticum (PoT) causes Wheat Blast, a devastating crop disease present in many wheat-producing countries. With the aim to contribute with the understanding of the mechanisms underlying pathogenesis of the fungus on wheat (Triticum aestivum), cytological and molecular studies were performed here.Susceptible and resistant interactions of wheat-PoT in cultivars of different levels of resistance were investigated by cellular changes and the molecular responses of wheat leaves according PR-1 expression. There was significant difference between susceptible (Apogee-PY15-PY34, Buck, MSINTA-PY15) and resistant (Buck-PY34, MSINTA-PY34) interactions in the responses tested. The two isolates of PoT analysed showed different ability to infect wheat plants at the cellular level during early tissue penetration analysed and documented by epi-fluorescence microscopy, patterns of PR-1 genes, and wheat plant defense reaction.Microscopic time-course analysis (24, 48, 72 and 96 h after inoculation) of individual interaction sites per leaf revealed different leaf invasion strategies of the more and less aggressive PoT isolates during the first stages of the course of pathogenesis.The results and conclusions of the present study provides the first evidence that different PR-1 spatial and temporal expression patterns occur in the early infective process of wheat-PoT, underlying two different defence mechanisms in plants according susceptible or resistant interactions. Thus, the results here laid a theoretical foundation for the future control of wheat blast using the different pattern of PR-1 gene as markers for disease resistance during the first stages of infection providing a significant contribution to a more efficient selection of wheat genotypes in breeding studies.

Spatial enrichment and genomic analyses reveal the link of NOMO1 with amyotrophic lateral sclerosis.

Amyotrophic lateral sclerosis (ALS) is a severe motor neuron disease with uncertain genetic predisposition in most sporadic cases. The spatial architecture of cell types and gene expression are the basis of cell-cell interactions, biological function and disease pathology, but are not well investigated in the human motor cortex, a key ALS-relevant brain region. Recent studies indicated single nucleus transcriptomic features of motor neuron vulnerability in ALS motor cortex. However, the brain regional vulnerability of ALS-associated genes and the genetic link between region-specific genes and ALS risk remain largely unclear. Here, we developed an entropy-weighted differential gene expression matrix-based tool (SpatialE) to identify the spatial enrichment of gene sets in spatial transcriptomics. We benchmarked SpatialE against another enrichment tool (multimodal intersection analysis) using spatial transcriptomics data from both human and mouse brain tissues. To investigate regional vulnerability, we analysed three human motor cortex and two dorsolateral prefrontal cortex tissues for spatial enrichment of ALS-associated genes. We also used Cell2location to estimate the abundance of cell types in ALS-related cortex layers. To dissect the link of regionally expressed genes and ALS risk, we performed burden analyses of rare loss-of-function variants detected by whole-genome sequencing in ALS patients and controls, then analysed differential gene expression in the TargetALS RNA-sequencing dataset. SpatialE showed more accurate and specific spatial enrichment of regional cell type markers than multimodal intersection analysis in both mouse brain and human dorsolateral prefrontal cortex. Spatial transcriptomic analyses of human motor cortex showed heterogeneous cell types and spatial gene expression profiles. We found that 260 manually curated ALS-associated genes are significantly enriched in layer 5 of the motor cortex, with abundant expression of upper motor neurons and layer 5 excitatory neurons. Burden analyses of rare loss-of-function variants in Layer 5-associated genes nominated NOMO1 as a novel ALS-associated gene in a combined sample set of 6814 ALS patients and 3324 controls (P = 0.029). Gene expression analyses in CNS tissues revealed downregulation of NOMO1 in ALS, which is consistent with a loss-of-function disease mechanism. In conclusion, our integrated spatial transcriptomics and genomic analyses identified regional brain vulnerability in ALS and the association of a layer 5 gene (NOMO1) with ALS risk.

The complex molecular epileptogenesis landscape of glioblastoma

The cortical microenvironment surrounding malignant glioblastoma is a source of depolarizing crosstalk favoring hyperexcitability, tumor expansion, and immune evasion. Neosynaptogenesis, excess glutamate, and altered intrinsic membrane currents contribute to excitability dyshomeostasis, yet only half of the cases develop seizures, suggesting that tumor and host genomics, along with location, rather than mass effect, play a critical role. We analyzed the spatial contours and expression of 358 clinically validated human epilepsy genes in the human glioblastoma transcriptome compared to non-tumor adult and developing cortex datasets. Nearly half, including dosage-sensitive genes whose expression levels are securely linked to monogenic epilepsy, are strikingly enriched and aberrantly regulated at the leading edge, supporting a complex epistatic basis for peritumoral epileptogenesis. Surround hyperexcitability induced by complex patterns of proepileptic gene expression may explain the limited efficacy of narrowly targeted antiseizure medicines and the persistence of epilepsy following tumor resection and clarify why not all brain tumors provoke seizures.

Heterogeneity of hepatocellular carcinoma that responds differently to combination therapy with TACE and Sorafenib as determined by digital spatial gene expression profiling.

The combination of Sorafenib and transcatheter arterial chemoembolization (TACE) exhibits limited efficacy in the treatment of certain advanced hepatocellular carcinomas (HCC), and the molecular mechanisms underlying resistance to this combination remain unclear. This study aims to underscore the distinctive contribution of GeoMx DSP technology in elucidating the molecular intricacies of HCC resistance to the Sorafenib and TACE combination. Patients with advanced HCC during the waiting period before liver transplantation were classified into sensitive and resistant groups based on their response to Sorafenib and TACE combination therapy. Employing GeoMx DSP technology for comprehensive gene expression profiling, we identified pivotal molecular targets linked to resistance against combination therapy. The investigation scrutinized intra-tumoral and inter-individual variances, unveiling a spectrum of crucial molecular targets, such as PLG, PLVAP, immunoglobulin genes, ORM1, and NR4A1, among others. Additionally, we explored signaling pathways associated with treatment responsiveness, including the PPAR signaling pathway. Notably, we emphasized the significance of the immune microenvironment characterized by heightened SPP1 expression in HCC resistance to combination therapy. In the resistant group, SPP1+ tumor-associated macrophage (TAM) infiltration was notably pronounced (p = 0.037), while T-cell depletion showed a mitigated presence (p = 0.013). The study reveals intra- and inter-individual heterogeneity in HCC that is differentially responsive to the combination of Sorafenib and TACE, highlighting multiple key molecular targets associated with treatment resistance. The immune microenvironment is important, and in particular, SPP1+ TAM infiltration may play a key role. Meanwhile, the introduction of immunotherapy in patients resistant to combination therapy may lead to positive results.

Calreticulin from Apostichopus japonicus relieves endoplasmic reticulum stress induced by Vibrio splendidus through autophagy

When organisms are exposed to external stimuli, misfolded proteins accumulate continuously, resulting in endoplasmic reticulum (ER) stress. Autophagy is of great significance for eliminating aggregated proteins and maintaining cellular homeostasis. However, the molecular mechanism of activating autophagy in response to ER stress in sea cucumber is remain unclear. In the current study, we demonstrated that the pathogen Vibrio splendidus can cause ER stress in Apostichopus japonicus coelomocytes and identified a Ca2+ binding partner calreticulin (designated as AjCRT), which increased with the occurrence of ER stress. The nucleotide sequence analysis showed that the open reading frame of AjCRT was 1242 bp and encoded a 413-amino-acid residue polyprotein with calreticulin domains. The spatial expression analysis revealed that AjCRT was ubiquitously expressed in all examined tissues with large magnitude in the coelomocytes and was minimally expressed in muscle. Furthermore, silencing AjCRT in vivo could significantly exacerbate ER stress induced by V. splendidus and resulted in the significant reduction of coelomocyte autophagy. These findings indicate a calreticulin-based mechanism that positively regulates autophagy in response to ER stress induced by pathogen infection. The results will provide a basis for understanding the way of host alleviating ER stress through autophagy, and pharmacological approaches may have potential for managing ER stress induced by pathogen and related cellular disorders.

Interneuron diversity in the human dorsal striatum

Deciphering the striatal interneuron diversity is key to understanding the basal ganglia circuit and to untangling the complex neurological and psychiatric diseases affecting this brain structure. We performed snRNA-seq and spatial transcriptomics of postmortem human caudate nucleus and putamen samples to elucidate the diversity and abundance of interneuron populations and their inherent transcriptional structure in the human dorsal striatum. We propose a comprehensive taxonomy of striatal interneurons with eight main classes and fourteen subclasses, providing their full transcriptomic identity and spatial expression profile as well as additional quantitative FISH validation for specific populations. We have also delineated the correspondence of our taxonomy with previous standardized classifications and shown the main transcriptomic and class abundance differences between caudate nucleus and putamen. Notably, based on key functional genes such as ion channels and synaptic receptors, we found matching known mouse interneuron populations for the most abundant populations, the recently described PTHLH and TAC3 interneurons. Finally, we were able to integrate other published datasets with ours, supporting the generalizability of this harmonized taxonomy.

Nuclear factors NF-YC3 and NF-YBs positively regulate arbuscular mycorrhizal symbiosis in tomato.

The involvement of nuclear factor Y (NF-Y) in transcriptional reprogramming during arbuscular mycorrhizal symbiosis has been demonstrated in several plant species. However, a comprehensive picture is lacking. We showed that the spatial expression of NF-YC3 was observed in cortical cells containing arbuscules via the cis-regulatory element GCC boxes. Moreover, the NF-YC3 promoter was transactivated by the combination of CYCLOPS and autoactive calcium and calmodulin-dependent kinase (CCaMK) via GCC boxes. Knockdown of NF-YC3 significantly reduced the abundance of all intraradical fungal structures and affected arbuscule size. BCP1, SbtM1, and WRI5a, whose expression associated with NF-YC3 levels, might be downstream of NF-YC3. NF-YC3 interacted with NF-YB3a, NF-YB5c, or NF-YB3b, in yeast (Saccharomyces cerevisiae) and in planta, and interacted with NF-YA3a in yeast. Spatial expression of three NF-YBs was observed in all cell layers of roots under both mock and mycorrhizal conditions. Simultaneous knockdown of three NF-YBs, but not individually, reduced the fungal colonization level, suggesting that there might be functional redundancy of NF-YBs to regulate AM symbiosis. Collectively, our data suggest that NF-YC3 and NF-YBs positively regulate AM symbiosis in tomato, and arbuscule-related NF-YC3 may be an important downstream gene of the common symbiosis signaling pathway.

Spatial and temporal gene expression patterns during early human odontogenesis process.

Studies on odontogenesis are of great importance to treat dental abnormalities and tooth loss. However, the odontogenesis process was poorly studied in humans, especially at the early developmental stages. Here, we combined RNA sequencing (RNA-seq) with Laser-capture microdissection (LCM) to establish a spatiotemporal transcriptomic investigation for human deciduous tooth germs at the crucial developmental stage to offer new perspectives to understand tooth development and instruct tooth regeneration. Several hallmark events, including angiogenesis, ossification, axonogenesis, and extracellular matrix (ECM) organization, were identified during odontogenesis in human dental epithelium and mesenchyme from the cap stage to the early bell stage. ECM played an essential role in the shift of tooth-inductive capability. Species comparisons demonstrated these hallmark events both in humans and mice. This study reveals the hallmark events during odontogenesis, enriching the transcriptomic research on human tooth development at the early stage.

Self-supervised learning for characterising histomorphological diversity and spatial RNA expression prediction across 23 human tissue types

As vast histological archives are digitised, there is a pressing need to be able to associate specific tissue substructures and incident pathology to disease outcomes without arduous annotation. Here, we learn self-supervised representations using a Vision Transformer, trained on 1.7 M histology images across 23 healthy tissues in 838 donors from the Genotype Tissue Expression consortium (GTEx). Using these representations, we can automatically segment tissues into their constituent tissue substructures and pathology proportions across thousands of whole slide images, outperforming other self-supervised methods (43% increase in silhouette score). Additionally, we can detect and quantify histological pathologies present, such as arterial calcification (AUROC = 0.93) and identify missing calcification diagnoses. Finally, to link gene expression to tissue morphology, we introduce RNAPath, a set of models trained on 23 tissue types that can predict and spatially localise individual RNA expression levels directly from H&E histology (mean genes significantly regressed = 5156, FDR 1%). We validate RNAPath spatial predictions with matched ground truth immunohistochemistry for several well characterised control genes, recapitulating their known spatial specificity. Together, these results demonstrate how self-supervised machine learning when applied to vast histological archives allows researchers to answer questions about tissue pathology, its spatial organisation and the interplay between morphological tissue variability and gene expression.

Melanoma progression and prognostic models drawn from single-cell, spatial maps of benign and malignant tumors.

Melanoma clinical outcomes emerge from incompletely understood genetic mechanisms operating within the tumor and its microenvironment. Here, we used single-cell RNA-based spatial molecular imaging (RNA-SMI) in patient-derived archival tumors to reveal clinically relevant markers of malignancy progression and prognosis. We examined spatial gene expression of 203,472 cells inside benign and malignant melanocytic neoplasms, including melanocytic nevi and primary invasive and metastatic melanomas. Algorithmic cell clustering paired with intratumoral comparative two-dimensional analyses visualized synergistic, spatial gene signatures linking cellular proliferation, metabolism, and malignancy, validated by protein expression. Metastatic niches included up-regulation of CDK2 and FABP5, which independently predicted poor clinical outcome in 473 patients with melanoma via Cox regression analysis. More generally, our work demonstrates a framework for applying single-cell RNA-SMI technology toward identifying gene regulatory landscapes pertinent to cancer progression and patient survival.

From star field affiliation to sundial theory: spatial expression in Qing dynasty Taiwan gazetteers

Star field affiliation is an important part of ancient Chinese gazetteers. The gazetteers of Taiwan in the Qing Dynasty presented many ‘different opinions’ about the star field affiliation of Taiwan. These different opinions illustrate the difficulty of self-consistency in traditional Chinese star field affiliation theory and reflect the influence of the concept of grand unification on traditional Chinese geographical knowledge. During the Qianlong period, with the introduction of Western surveying and mapping knowledge, the sundial began to replace star field affiliation and became a new geographical positioning theory. Taiwan, whose geographical positioning also completed the transformation from the star field affiliation to sundial theory at the end of the Qing Dynasty, witnessed the replacement of traditional Chinese place knowledge with cosmopolitan knowledge. Looking back at the writings regarding the star field affiliation in Taiwan gazetteers in the Qing Dynasty, there was both debate and consensus.

Decoding the Epigenetic Landscape: Insights into 5mC and 5hmC Patterns in Mouse Cortical Cell Types.

The DNA modifications, 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC), represent powerful epigenetic regulators of temporal and spatial gene expression. Yet, how the cooperation of these genome-wide, epigenetic marks determine unique transcriptional signatures across different brain cell populations is unclear. Here we applied Nanopore sequencing of native DNA to obtain a complete, genome-wide, single-base resolution atlas of 5mC and 5hmC modifications in neurons, astrocytes and microglia in the mouse cortex (99% genome coverage, 40 million CpG sites). In tandem with RNA sequencing, analysis of 5mC and 5hmC patterns across cell types reveals astrocytes drive uniquely high brain 5hmC levels and support two decades of research regarding methylation patterns, gene expression and alternative splicing, benchmarking this resource. As such, we provide the most comprehensive DNA methylation data in mouse brain as an interactive, online tool (NAM-Me, https://olsenlab.shinyapps.io/NAMME/) to serve as a resource dataset for those interested in the methylome landscape.