Comparative Methods for Demystifying Spatial Transcriptomics

Spatial Transcriptomics (ST), coined as the term for parallel RNA-Seq on cell populations ordered spatially on a histological tissue section, has recently become increasingly popular, especially in experiments where microfluidics-based single-cell sequencing fails, such as assays on neurons. ST platforms, like the 10x Visium technology investigated herein, therefore produce in a single experiment simultaneously thousands of RNA readouts, captured by an array of micrometer scale spots under the histological section. Therefore, a central challenge of analyzing ST experiments consists of analyzing the gene expression morphology of all spots to delineate clusters of similar cell mixtures, which are then compared to each other to identify up- or down-regulated marker genes. Moreover, another level of complexity in ST experiments, compared to traditional RNA-Seq, is imposed by staining the tissue section with protein markers of cells or cell components to identify spots providing relevant information afterward. The corresponding microscopy images need to be analyzed in addition to the RNA-Seq read mappings on the reference genome and transcriptome sequences. Focusing on the software suite provided by the Visium platform manufacturer, we break down the ST analysis pipeline into its four essential steps-the image analysis, the read alignment, the gene quantification, and the spot clustering-and compare results obtained when using reads from different subsets of spots and/or when employing alternative genome or transcriptome references. Our comparative analyses demonstrate the impact of spot selection and the choice of genome/transcriptome references on the analysis results when employing the manufacturer's pipeline.

Designing Oligonucleotide-Based FISH Probe Sets with PaintSHOP

Fluorescent in situ hybridization (FISH) enables the visualization of the position and abundance of nucleic acid molecules in fixed cell and tissue samples. Many FISH-based methods employ sets of synthetic, computationally designed DNA oligonucleotide (oligo) FISH probes, including massively multiplexed imaging spatial transcriptomics and spatial genomics technologies. Oligo probes can either be designed de novo or accessed from an existing database of pre-discovered probe sequences. This chapter describes the use of PaintSHOP, a user-friendly, web-based platform for the design of sets of oligo-based FISH probes. PaintSHOP hosts large collections of pre-discovered probes from many model organisms and also provides collections of functional sequences such as primers and readout domains and interactive tools to add these functional sequences to selected probes. Detailed examples are provided for three common experimental scenarios.

ENGEP: advancing spatial transcriptomics with accurate unmeasured gene expression prediction

Imaging-based spatial transcriptomics techniques provide valuable spatial and gene expression information at single-cell resolution. However, their current capability is restricted to profiling a limited number of genes per sample, resulting in most of the transcriptome remaining unmeasured. To overcome this challenge, we develop ENGEP, an ensemble learning-based tool that predicts unmeasured gene expression in spatial transcriptomics data by using multiple single-cell RNA sequencing datasets as references. ENGEP outperforms current state-of-the-art tools and brings biological insight by accurately predicting unmeasured genes. ENGEP has exceptional efficiency in terms of runtime and memory usage, making it scalable for analyzing large datasets.

Exploring distinct properties of endometrial stem cells through advanced single-cell analysis platforms

The endometrium is a dynamic tissue that undergoes cyclic changes in response to ovarian hormones during the menstrual cycle. These changes are crucial for pregnancy establishment and maintenance. Endometrial stem cells play a pivotal role in endometrial regeneration and repair by differentiating into various cell types within the endometrium. However, their involvement in endometrial disorders such as endometriosis, infertility, and endometrial cancer is still not fully understood yet. Traditional bulk sequencing methods have limitations in capturing heterogeneity and complexity of endometrial stem cell populations. To overcome these limitations, recent single-cell analysis techniques, including single-cell RNA sequencing (scRNA-Seq), single-cell ATAC sequencing (scATAC-Seq), and spatial transcriptomics, have emerged as valuable tools for studying endometrial stem cells. In this review, although there are still many technical limitations that require improvement, we will summarize the current state-of-the-art single-cell analysis techniques for endometrial stem cells and explore their relevance to related diseases. We will discuss studies utilizing various single-cell analysis platforms to identify and characterize distinct endometrial stem cell populations and investigate their dynamic changes in gene expression and epigenetic patterns during menstrual cycle and differentiation processes. These techniques enable the identification of rare cell populations, capture heterogeneity of cell populations within the endometrium, and provide potential targets for more effective therapies.

Unveiling the hidden AP-1: revealing the crucial role of AP-1 in ccRCC at single-cell resolution

Clear cell renal cell carcinoma (ccRCC), as the most common histological subtype of kidney cancer, has been reported to originate primarily from proximal tubule (PT) cells in the kidney. However, the current research on its associated molecular mechanisms remains relatively limited. In our study, we analyzed multiple single-cell multi-omics datasets obtained from various research teams, revealing the significant role of the activator protein 1 (AP-1) in ccRCC tumorigenesis. The motif activity analysis of transcription factors (TFs) showed a predominant activation of AP-1 in ccRCC cancer cells compared to PT cells. Furthermore, our findings at single-cell resolution revealed a notable absence of AP-1 expression in PT cells when compared to ccRCC cancer cells. In bulk-RNA of discovery cohort, no differential expression of AP-1 was detected in normal kidney and ccRCC samples, which may be attributed to confounding effects in bulk-RNA sequencing. Meanwhile, spatial transcriptomics analysis demonstrated a broader expression range of the AP-1 compared to the ccRCC marker CA9. Moreover, we observed chromatin accessibility of the AP-1 in various cell-types, including PT cells, suggesting that the transcriptional expression of AP-1 in PT cells may be influenced by subsequent transcriptional modifications, reflecting the complex regulatory mechanism of AP-1 transcription. These findings provide important insights for a deeper understanding of the function and regulatory mechanisms of AP-1 in ccRCC, thereby establishing a theoretical foundation for future clinical research and the development of treatment strategies.

Acute ischemia induces spatially and transcriptionally distinct microglial subclusters

BACKGROUND

Damage in the ischemic core and penumbra after stroke affects patient prognosis. Microglia immediately respond to ischemic insult and initiate immune inflammation, playing an important role in the cellular injury after stroke. However, the microglial heterogeneity and the mechanisms involved remain unclear.

METHODS

We first performed single-cell RNA-sequencing (scRNA-seq) and spatial transcriptomics (ST) on middle cerebral artery occlusion (MCAO) mice from three time points to determine stroke-associated microglial subclusters and their spatial distributions. Furthermore, the expression of microglial subcluster-specific marker genes and the localization of different microglial subclusters were verified on MCAO mice through RNAscope and immunofluorescence. Gene set variation analysis (GSVA) was performed to reveal functional characteristics of microglia sub-clusters. Additionally, ingenuity pathway analysis (IPA) was used to explore upstream regulators of microglial subclusters, which was confirmed by immunofluorescence, RT-qPCR, shRNA-mediated knockdown, and targeted metabolomics. Finally, the infarct size, neurological deficits, and neuronal apoptosis were evaluated in MCAO mice after manipulation of specific microglial subcluster.

RESULTS

We discovered stroke-associated microglial subclusters in the brains of MCAO mice. We also identified novel marker genes of these microglial subclusters and defined these cells as ischemic core-associated (ICAM) and ischemic penumbra-associated (IPAM) microglia, according to their spatial distribution. ICAM, induced by damage-associated molecular patterns, are probably fueled by glycolysis, and exhibit increased pro-inflammatory cytokines and chemokines production. BACH1 is a key transcription factor driving ICAM generation. In contrast, glucocorticoids, which are enriched in the penumbra, likely trigger IPAM formation, which are presumably powered by the citrate cycle and oxidative phosphorylation and are characterized by moderate pro-inflammatory responses, inflammation-alleviating metabolic features, and myelinotrophic properties.

CONCLUSIONS

ICAM could induce excessive neuroinflammation, aggravating brain injury, whereas IPAM probably exhibit neuroprotective features, which could be essential for the homeostasis and survival of cells in the penumbra. Our findings provide a biological basis for targeting specific microglial subclusters as a potential therapeutic strategy for ischemic stroke.

An interactive resource of molecular signalling in the developing human haematopoietic stem cell niche

The emergence of definitive human haematopoietic stem cells (HSCs) from Carnegie Stage (CS) 14 to CS17 in the aorta-gonad-mesonephros (AGM) region is a tightly regulated process. Previously, we conducted spatial transcriptomic analysis of the human AGM region at the end of this period (CS16/CS17) and identified secreted factors involved in HSC development. Here, we extend our analysis to investigate the progression of dorso-ventral polarised signalling around the dorsal aorta over the entire period of HSC emergence. Our results reveal a dramatic increase in ventral signalling complexity from the CS13-CS14 transition, coinciding with the first appearance of definitive HSCs. We further observe stage-specific changes in signalling up to CS17, which may underpin the step-wise maturation of HSCs described in the mouse model. The data-rich resource is also presented in an online interface enabling in silico analysis of molecular interactions between spatially defined domains of the AGM region. This resource will be of particular interest for researchers studying mechanisms underlying human HSC development as well as those developing in vitro methods for the generation of clinically relevant HSCs from pluripotent stem cells.

Hexagonal image segmentation on spatially resolved transcriptomics

Spatial transcriptomics is a rapidly evolving field that enables researchers to capture comprehensive molecular profiles while preserving information about the physical locations. One major challenge in this research area involves the identification of spatial domains, which are distinct regions characterized by unique gene expression patterns. However, current unsupervised methods have struggled to perform well in this regard due to the presence of high levels of noise and dropout events in spatial transcriptomic profiles. In this paper, we propose a novel hexagonal Convolutional Neural Network (hexCNN) for hexagonal image segmentation on spatially resolved transcriptomics. To address the problem of noise and dropout occurrences within spatial transcriptomics data, we first extend an unsupervised algorithm to a supervised learning method that can identify useful features and reduce noise hindrance. Then, inspired by the classical convolution in convolutional neural networks (CNNs), we designed a regular hexagonal convolution to compensate for the missing gene expression patterns from adjacent spots. We evaluated the performance of hexCNN by applying it to the DLPFC dataset. The results show that hexCNN achieves a classification accuracy of 86.8% and an average Rand index (ARI) of 77.1% (1.4% and 2.5% higher than those of GNNs). The results also demonstrate that hexCNN is capable of removing the noise caused by batch effect while preserving the biological signal differences.

Transformation of non-neuritic into neuritic plaques during AD progression drives cortical spread of tau pathology via regenerative failure

Extracellular amyloid-β (Aβ) plaques and intracellular aggregates of tau protein in form of neurofibrillary tangles (NFT) are pathological hallmarks of Alzheimer's disease (AD). The exact mechanism how these two protein aggregates interact in AD is still a matter of debate. Neuritic plaques (NP), a subset of Aβ plaques containing dystrophic neurites (DN), are suggested to be unique to AD and might play a role in the interaction of Aβ and tau. Quantifying NP and non-NP in postmortem brain specimens from patients with increasing severity of AD neuropathological changes (ADNC), we demonstrate that the total number of Aβ plaques and NP increase, while the number of non-NP stagnates. Furthermore, investigating the correlation between NP and NFT, we identified unexpected brain region-specific differences when comparing cases with increasingly more severe ADNC. In neocortical regions NFT counts increase in parallel with NP counts during the progression of ADNC, while this correlation is not observed in hippocampus. These data support the notion that non-NP are transformed into NP during the progression of ADNC and indicate that NP might drive cortical NFT formation. Next, using spatial transcriptomics, we analyzed the gene expression profile of the microenvironment around non-NP and NP. We identified an upregulation of neuronal systems and Ca-dependent event pathways around NP compared to non-NP. We speculate that the upregulation of these transcripts may hint at a compensatory mechanism underlying NP formation. Our studies suggest that the transformation of non-NP to NP is a key event in ADNC progression and points to regenerative failure as a potential driving force of this process.

The role of LSM1 in breast cancer: Shaping metabolism and tumor-associated macrophage infiltration

LSM1 is part of the cytoplasmic protein complex Lsm1-7-Pat1 and is likely involved in pre-mRNA degradation by aiding U4/U6 snRNP formation. More research is needed to uncover LSM1's potential in breast cancer (BRCA) clinical pathology, the tumor immune microenvironment, and precision oncology. We discovered LSM1 as a diagnostic marker for advanced BRCA with poor survival, using a multi-omics approach. We studied LSM1 expression across BRCA regions and its link to immune cells through various methods, including spatial transcriptomics and single-cell RNA-sequencing. We also examined how silencing LSM1 affects mitochondrial function and energy metabolism in the tumor environment. These findings were confirmed using 54 BRCA patient biopsies and tissue microarrays. Immunofluorescence and bioinformatics assessed LSM1's connection to clinicopathological features and prognosis. This study uncovers gene patterns linked to breast cancer, with LSM1 linked to macrophage energy processes. Silencing LSM1 in breast cancer cells disrupts mitochondria and energy metabolism. Spatial analysis aligns with previous results, showing LSM1's connection to macrophages. Biopsies confirm LSM1 elevation in advanced breast cancer with increased macrophage presence. To summarize, LSM1 changes may drive BRCA progression, making it a potential diagnostic and prognostic marker. It also influences energy metabolism and the tumor's immune environment during metastasis, showing promise for precision medicine and drug screening in BRCA.

Talniflumate abrogates mucin immune suppressive barrier improving efficacy of gemcitabine and nab-paclitaxel treatment in pancreatic cancer

BACKGROUND

Pancreatic ductal adenocarcinoma (PDAC) is a lethal disease. This is due to its aggressive course, late diagnosis and its intrinsic drugs resistance. The complexity of the tumor, in terms of cell components and heterogeneity, has led to the approval of few therapies with limited efficacy. The study of the early stages of carcinogenesis provides the opportunity for the identification of actionable pathways that underpin therapeutic resistance.

METHODS

We analyzed 43 Intraductal papillary mucinous neoplasms (IPMN) (12 Low-grade and 31 High-grade) by Spatial Transcriptomics. Mouse and human pancreatic cancer organoids and T cells interaction platforms were established to test the role of mucins expression on T cells activity. Syngeneic mouse model of PDAC was used to explore the impact of mucins downregulation on standard therapy efficacy.

RESULTS

Spatial transcriptomics showed that mucin O-glycosylation pathway is increased in the progression from low-grade to high-grade IPMN. We identified GCNT3, a master regulator of mucins expression, as an actionable target of this pathway by talniflumate. We showed that talniflumate impaired mucins expression increasing T cell activation and recognition using both mouse and human organoid interaction platforms. In vivo experiments showed that talniflumate was able to increase the efficacy of the chemotherapy by boosting immune infiltration.

CONCLUSIONS

Finally, we demonstrated that combination of talniflumate, an anti-inflammatory drug, with chemotherapy effectively improves anti-tumor effect in PDAC.

A SELECTIVE REVIEW OF RECENT DEVELOPMENTS IN SPATIALLY VARIABLE GENE DETECTION FOR SPATIAL TRANSCRIPTOMICS

With the emergence of advanced spatial transcriptomic technologies, there has been a surge in research papers dedicated to analyzing spatial transcriptomics data, resulting in significant contributions to our understanding of biology. The initial stage of downstream analysis of spatial transcriptomic data has centered on identifying spatially variable genes (SVGs) or genes expressed with specific spatial patterns across the tissue. SVG detection is an important task since many downstream analyses depend on these selected SVGs. Over the past few years, a plethora of new methods have been proposed for the detection of SVGs, accompanied by numerous innovative concepts and discussions. This article provides a selective review of methods and their practical implementations, offering valuable insights into the current literature in this field.

Spatial transcriptomics and neurofilament light chain reveal changes in lesion patterns in murine autoimmune neuroinflammation

OBJECTIVE

Ongoing neuroaxonal damage is a major contributor to disease progression and long-term disability in multiple sclerosis. However, spatio-temporal distribution and pathophysiological mechanisms of neuroaxonal damage during acute relapses and later chronic disease stages remain poorly understood.

METHODS

Here, we applied immunohistochemistry, single-molecule array, spatial transcriptomics, and microglia/axon co-cultures to gain insight into spatio-temporal neuroaxonal damage in experimental autoimmune encephalomyelitis (EAE).

RESULTS

Association of spinal cord white matter lesions and blood-based neurofilament light (sNfL) levels revealed a distinct, stage-dependent anatomical pattern of neuroaxonal damage: in chronic EAE, sNfL levels were predominately associated with anterolateral lumbar lesions, whereas in early EAE sNfL showed no correlation with lesions in any anatomical location. Furthermore, neuroaxonal damage in late EAE was largely confined to white matter lesions but showed a widespread distribution in early EAE. Following this pattern of neuroaxonal damage, spatial transcriptomics revealed a widespread cyto- and chemokine response at early disease stages, whereas late EAE was characterized by a prominent glial cell accumulation in white matter lesions. These findings were corroborated by immunohistochemistry and microglia/axon co-cultures, which further revealed a strong association between CNS myeloid cell activation and neuroaxonal damage both in vivo and in vitro.

INTERPRETATION

Our findings indicate that CNS myeloid cells may play a crucial role in driving neuroaxonal damage in EAE. Moreover, neuroaxonal damage can progress in a stage-dependent centripetal manner, transitioning from normal-appearing white matter to focal white matter lesions. These insights may contribute to a better understanding of neurodegeneration and elevated sNfL levels observed in multiple sclerosis patients at different disease stages.

Unraveling temporal and spatial biomarkers of epithelial-mesenchymal transition in colorectal cancer: insights into the crucial role of immunosuppressive cells

The occurrence and progression of tumors can be established through a complex interplay among tumor cells undergoing epithelial-mesenchymal transition (EMT), invasive factors and immune cells. In this study, we employed single-cell RNA sequencing (scRNA-seq) and spatially resolved transcriptomics (ST) to evaluate the pseudotime trajectory and spatial interactive relationship between EMT-invasive malignant tumors and immune cells in primary colorectal cancer (CRC) tissues at different stages (stage I/II and stage III with tumor deposit). Our research characterized the spatiotemporal relationship among different invasive tumor programs by constructing pseudotime endpoint-EMT-invasion tumor programs (EMTPs) located at the edge of ST, utilizing evolution trajectory analysis integrated with EMT-invasion genes. Strikingly, the invasive and expansive process of tumors undergoes remarkable spatial reprogramming of regulatory and immunosuppressive cells, such as myeloid-derived suppressor cells (MDSCs), tumor-associated macrophages (TAMs), regulatory T cells (Treg), and exhausted T cells (Tex). These EMTP-adjacent cell are linked to EMT-related invasion genes, especially the C-X-C motif ligand 1 (CXCL1) and CXCL8 genes that are important for CRC prognosis. Interestingly, the EMTPs in stage I mainly produce an inflammatory margin invasive niche, while the EMTPs in stage III tissues likely produce a hypoxic pre-invasive niche. Our data demonstrate the crucial role of regulatory and immunosuppressive cells in tumor formation and progression of CRC. This study provides a framework to delineate the spatiotemporal invasive niche in CRC samples.

Architecture of skin inflammation in psoriasis revealed by spatial transcriptomics

Psoriasis is a chronic inflammatory skin disease, thought to be predominantly mediated by TH17 cells. Significance of other inflammatory pathways and the innate immune system is not well understood and the spatial heterogeneity of inflammation in the skin has largely been overlooked. Our aim was to create a comprehensive map of skin inflammation in psoriasis, exploring the tissue patterning of inflammation. In situ whole transcriptome sequencing (spatial sequencing) was performed on lesional psoriatic skin in four patients with moderate-to-severe disease to quantify all expressed genes within a tissue section. Transcriptional analysis revealed three major inflammatory niches in psoriasis skin, each with distinct cytokine circuits and chemokines: the hyperplastic epidermis, upper (papillary) dermis, and reticular dermis. Interestingly, key cytokines such as IL-23, IL-17 s, and TNFα were not notably present in the skin's transcriptomic signature. Unexpectedly, IL-32 showed strong expression in the dermis. Our findings underscore the complexity of psoriatic inflammation, highlighting its architectural heterogeneity and the roles of innate cytokines. Both IL-32 and IL-1 family cytokines appear to play critical roles in the dermal and epidermal inflammation, respectively, and may provide pharmacological targets to improve the control of the inflammatory process.

Microglial contribution to the pathology of neurodevelopmental disorders in humans

Microglia are the brain's resident macrophages, which guide various developmental processes crucial for brain maturation, activity, and plasticity. Microglial progenitors enter the telencephalic wall by the 4th postconceptional week and colonise the fetal brain in a manner that spatiotemporally tracks key neurodevelopmental processes in humans. However, much of what we know about how microglia shape neurodevelopment comes from rodent studies. Multiple differences exist between human and rodent microglia warranting further focus on the human condition, particularly as microglia are emerging as critically involved in the pathological signature of various cognitive and neurodevelopmental disorders. In this article, we review the evidence supporting microglial involvement in basic neurodevelopmental processes by focusing on the human species. We next concur on the neuropathological evidence demonstrating whether and how microglia contribute to the aetiology of two neurodevelopmental disorders: autism spectrum conditions and schizophrenia. Next, we highlight how recent technologies have revolutionised our understanding of microglial biology with a focus on how these tools can help us elucidate at unprecedented resolution the links between microglia and neurodevelopmental disorders. We conclude by reviewing which current treatment approaches have shown most promise towards targeting microglia in neurodevelopmental disorders and suggest novel avenues for future consideration.

Combined near infrared photoacoustic imaging and ultrasound detects vulnerable atherosclerotic plaque

Atherosclerosis is an inflammatory process resulting in the deposition of cholesterol and cellular debris, narrowing of the vessel lumen and clot formation. Characterization of the morphology and vulnerability of the lesion is essential for effective clinical management. Here, near-infrared auto-photoacoustic (NIRAPA) imaging is shown to detect plaque components and, when combined with ultrasound imaging, to differentiate stable and vulnerable plaque. In an ex vivo study of photoacoustic imaging of excised plaque from 25 patients, 88.2% sensitivity and 71.4% specificity were achieved using a clinically-relevant protocol. In order to determine the origin of the NIRAPA signal, immunohistochemistry, spatial transcriptomics and spatial proteomics were co-registered with imaging and applied to adjacent plaque sections. The highest NIRAPA signal was spatially correlated with bilirubin and associated blood-based residue and with the cytoplasmic contents of inflammatory macrophages bearing CD74, HLA-DR, CD14 and CD163 markers. In summary, we establish the potential to apply the NIRAPA-ultrasound imaging combination to detect vulnerable carotid plaque and a methodology for fusing molecular imaging with spatial transcriptomic and proteomic methods.

Spatial enrichment of the type 1 interferon signature in the brain of a neuropsychiatric lupus murine model

Among systemic lupus erythematosus (SLE) patients, neuropsychiatric symptoms are highly prevalent, being observed in up to 80% of adult and 95% of pediatric patients. Type 1 interferons, particularly interferon alpha (IFNα), have been implicated in the pathogenesis of SLE and its associated neuropsychiatric symptoms (NPSLE). However, it remains unclear how type 1 interferon signaling in the central nervous system (CNS) might result in neuropsychiatric sequelae. In this study, we validate an NPSLE mouse model and find an elevated peripheral type 1 interferon signature alongside clinically relevant NPSLE symptoms such as anxiety and fatigue. Unbiased single-nucleus sequencing of the hindbrain and hippocampus revealed that interferon-stimulated genes (ISGs) were among the most highly upregulated genes in both regions and that gene pathways involved in cellular interaction and neuronal development were generally repressed among astrocytes, oligodendrocytes, and neurons. Using image-based spatial transcriptomics, we found that the type 1 interferon signature is enriched as spatially distinct patches within the brain parenchyma of these mice. Our results suggest that type 1 interferon in the CNS may play an important mechanistic role in mediating NPSLE behavioral phenotypes by repressing general cellular communication pathways, and that type 1 interferon signaling modulators are a potential therapeutic option for NPSLE.

GeneSegNet: a deep learning framework for cell segmentation by integrating gene expression and imaging

When analyzing data from in situ RNA detection technologies, cell segmentation is an essential step in identifying cell boundaries, assigning RNA reads to cells, and studying the gene expression and morphological features of cells. We developed a deep-learning-based method, GeneSegNet, that integrates both gene expression and imaging information to perform cell segmentation. GeneSegNet also employs a recursive training strategy to deal with noisy training labels. We show that GeneSegNet significantly improves cell segmentation performances over existing methods that either ignore gene expression information or underutilize imaging information.

Spatial transcriptomics reveal basal sex differences in supraoptic nucleus gene expression of adult rats related to cell signaling and ribosomal pathways

BACKGROUND

The supraoptic nucleus (SON) of the hypothalamus contains magnocellular neurosecretory cells that secrete the hormones vasopressin and oxytocin. Sex differences in SON gene expression have been relatively unexplored. Our study used spatially resolved transcriptomics to visualize gene expression profiles in the SON of adult male (n = 4) and female (n = 4) Sprague-Dawley rats using Visium Spatial Gene Expression (10x Genomics).

METHODS

Briefly, 10-μm coronal sections (~ 4 × 4 mm) containing the SON were collected from each rat and processed using Visium slides and recommended protocols. Data were analyzed using 10x Genomics' Space Ranger and Loupe Browser applications and other bioinformatic tools. Two unique differential expression (DE) analysis methods, Loupe Browser and DESeq2, were used.

RESULTS

Loupe Browser DE analysis of the SON identified 116 significant differentially expressed genes (DEGs) common to both sexes (e.g., Avp and Oxt), 31 significant DEGs unique to the males, and 73 significant DEGs unique to the females. DESeq2 analysis revealed 183 significant DEGs between the two groups. Gene Ontology (GO) enrichment and pathway analyses using significant genes identified via Loupe Browser revealed GO terms and pathways related to (1) neurohypophyseal hormone activity, regulation of peptide hormone secretion, and regulation of ion transport for the significant genes common to both males and females, (2) Gi signaling/G-protein mediated events for the significant genes unique to males, and (3) potassium ion transport/voltage-gated potassium channels for the significant genes unique to females, as some examples. GO/pathway analyses using significant genes identified via DESeq2 comparing female vs. male groups revealed GO terms/pathways related to ribosomal structure/function. Ingenuity Pathway Analysis (IPA) identified additional sex differences in canonical pathways (e.g., 'Mitochondrial Dysfunction', 'Oxidative Phosphorylation') and upstream regulators (e.g., CSF3, NFKB complex, TNF, GRIN3A).

CONCLUSION

There was little overlap in the IPA results for the two different DE methods. These results suggest sex differences in SON gene expression that are associated with cell signaling and ribosomal pathways.

Dual spatially resolved transcriptomics for human host-pathogen colocalization studies in FFPE tissue sections

Technologies to study localized host-pathogen interactions are urgently needed. Here, we present a spatial transcriptomics approach to simultaneously capture host and pathogen transcriptome-wide spatial gene expression information from human formalin-fixed paraffin-embedded (FFPE) tissue sections at a near single-cell resolution. We demonstrate this methodology in lung samples from COVID-19 patients and validate our spatial detection of SARS-CoV-2 against RNAScope and in situ sequencing. Host-pathogen colocalization analysis identified putative modulators of SARS-CoV-2 infection in human lung cells. Our approach provides new insights into host response to pathogen infection through the simultaneous, unbiased detection of two transcriptomes in FFPE samples.

Pan-cancer spatially resolved single-cell analysis reveals the crosstalk between cancer-associated fibroblasts and tumor microenvironment

Cancer-associated fibroblasts (CAFs) are a heterogeneous cell population that plays a crucial role in remodeling the tumor microenvironment (TME). Here, through the integrated analysis of spatial and single-cell transcriptomics data across six common cancer types, we identified four distinct functional subgroups of CAFs and described their spatial distribution characteristics. Additionally, the analysis of single-cell RNA sequencing (scRNA-seq) data from three additional common cancer types and two newly generated scRNA-seq datasets of rare cancer types, namely epithelial-myoepithelial carcinoma (EMC) and mucoepidermoid carcinoma (MEC), expanded our understanding of CAF heterogeneity. Cell-cell interaction analysis conducted within the spatial context highlighted the pivotal roles of matrix CAFs (mCAFs) in tumor angiogenesis and inflammatory CAFs (iCAFs) in shaping the immunosuppressive microenvironment. In patients with breast cancer (BRCA) undergoing anti-PD-1 immunotherapy, iCAFs demonstrated heightened capacity in facilitating cancer cell proliferation, promoting epithelial-mesenchymal transition (EMT), and contributing to the establishment of an immunosuppressive microenvironment. Furthermore, a scoring system based on iCAFs showed a significant correlation with immune therapy response in melanoma patients. Lastly, we provided a web interface ( https://chenxisd.shinyapps.io/pancaf/ ) for the research community to investigate CAFs in the context of pan-cancer.

Spatial transcriptomic analysis of virtual prostate biopsy reveals confounding effect of tissue heterogeneity on genomic signatures

Genetic signatures have added a molecular dimension to prognostics and therapeutic decision-making. However, tumour heterogeneity in prostate cancer and current sampling methods could confound accurate assessment. Based on previously published spatial transcriptomic data from multifocal prostate cancer, we created virtual biopsy models that mimic conventional biopsy placement and core size. We then analysed the gene expression of different prognostic signatures (OncotypeDx®, Decipher®, Prostadiag®) using a step-wise approach with increasing resolution from pseudo-bulk analysis of the whole biopsy, to differentiation by tissue subtype (benign, stroma, tumour), followed by distinct tumour grade and finally clonal resolution. The gene expression profile of virtual tumour biopsies revealed clear differences between grade groups and tumour clones, compared to a benign control, which were not reflected in bulk analyses. This suggests that bulk analyses of whole biopsies or tumour-only areas, as used in clinical practice, may provide an inaccurate assessment of gene profiles. The type of tissue, the grade of the tumour and the clonal composition all influence the gene expression in a biopsy. Clinical decision making based on biopsy genomics should be made with caution while we await more precise targeting and cost-effective spatial analyses.

The promising application of cell-cell interaction analysis in cancer from single-cell and spatial transcriptomics

Cell-cell interactions instruct cell fate and function. These interactions are hijacked to promote cancer development. Single-cell transcriptomics and spatial transcriptomics have become powerful new tools for researchers to profile the transcriptional landscape of cancer at unparalleled genetic depth. In this review, we discuss the rapidly growing array of computational tools to infer cell-cell interactions from non-spatial single-cell RNA-sequencing and the limited but growing number of methods for spatial transcriptomics data. Downstream analyses of these computational tools and applications to cancer studies are highlighted. We finish by suggesting several directions for further extensions that anticipate the increasing availability of multi-omics cancer data.

Emerging Concepts in Precision Medicine in Axial Spondyloarthritis

PURPOSE OF REVIEW

Axial spondyloarthritis (AxSpA) is among the rheumatology's most heritable complex diseases, yet precision medicine at clinics still needs to be explored. We reviewed the emerging concepts and recent developments in polygenic risk scores, Mendelian randomization, pharmacogenomics, single-cell sequencing, and spatial transcriptomics.

RECENT FINDINGS

Polygenic risk score has resulted in encouraging results with potential diagnostic utility as it appears to outperform current diagnostic tools. Its performance and generalizability vary with ethnicity. Mendelian randomization has elucidated multiple causal associations, particularly between inflammatory bowel disease and AxSpA. Single-cell transcriptomics (particularly scRNA-seq and scATAC-seq) has identified numerous cell types, including synovial and blood immunological cells, to understand the contribution of both innate and adaptative immunity in AxSpA. Current molecular tools provide an exciting opportunity to advance precision medicine for AxSpA patients. However, extensive research and implementation strategies are still required before they can have an impact in the clinic.