Library size confounds biology in spatial transcriptomics data.

Spatial molecular data has transformed the study of disease microenvironments, though, larger datasets pose an analytics challenge prompting the direct adoption of single-cell RNA-sequencing tools including normalization methods. Here, we demonstrate that library size is associated with tissue structure and that normalizing these effects out using commonly applied scRNA-seq normalization methods will negatively affect spatial domain identification. Spatial data should not be specifically corrected for library size prior to analysis, and algorithms designed for scRNA-seq data should be adopted with caution.

Spatial omics techniques and data analysis for cancer immunotherapy applications.

In-depth profiling of cancer cells/tissues is expanding our understanding of the genomic, epigenomic, transcriptomic, and proteomic landscape of cancer. However, the complexity of the cancer microenvironment, particularly its immune regulation, has made it difficult to exploit the potential of cancer immunotherapy. High-throughput spatial omics technologies and analysis pipelines have emerged as powerful tools for tackling this challenge. As a result, a potential revolution in cancer diagnosis, prognosis, and treatment is on the horizon. In this review, we discuss the technological advances in spatial profiling of cancer around and beyond the central dogma to harness the full benefits of immunotherapy. We also discuss the promise and challenges of spatial data analysis and interpretation and provide an outlook for the future.

In situ single-cell profiling sheds light on IFI27 localisation during SARS-CoV-2 infection.

The utilization of single-cell resolved spatial transcriptomics to delineate immune responses during SARS-CoV-2 infection was able to identify M1 macrophages to have elevated expression of IFI27 in areas of infection.

vissE: a versatile tool to identify and visualise higher-order molecular phenotypes from functional enrichment analysis.

Functional analysis of high throughput experiments using pathway analysis is now ubiquitous. Though powerful, these methods often produce thousands of redundant results owing to knowledgebase redundancies upstream. This scale of results hinders extensive exploration by biologists and can lead to investigator biases due to previous knowledge and expectations. To address this issue, we present vissE, a flexible network-based analysis and visualisation tool that organises information into semantic categories and provides various visualisation modules to characterise them with respect to the underlying data, thus providing a comprehensive view of the biological system. We demonstrate vissE's versatility by applying it to three different technologies: bulk, single-cell and spatial transcriptomics. Applying vissE to a factor analysis of a breast cancer spatial transcriptomic data, we identified stromal phenotypes that support tumour dissemination. Its adaptability allows vissE to enhance all existing gene-set enrichment and pathway analysis workflows, empowering biologists during molecular discovery.

standR: spatial transcriptomic analysis for GeoMx DSP data.

To gain a better understanding of the complexity of gene expression in normal and diseased tissues it is important to account for the spatial context and identity of cells in situ. State-of-the-art spatial profiling technologies, such as the Nanostring GeoMx Digital Spatial Profiler (DSP), now allow quantitative spatially resolved measurement of the transcriptome in tissues. However, the bioinformatics pipelines currently used to analyse GeoMx data often fail to successfully account for the technical variability within the data and the complexity of experimental designs, thus limiting the accuracy and reliability of the subsequent analysis. Carefully designed quality control workflows, that include in-depth experiment-specific investigations into technical variation and appropriate adjustment for such variation can address this issue. Here, we present standR, an R/Bioconductor package that enables an end-to-end analysis of GeoMx DSP data. With four case studies from previously published experiments, we demonstrate how the standR workflow can enhance the statistical power of GeoMx DSP data analysis and how the application of standR enables scientists to develop in-depth insights into the biology of interest.

Spatial proteomics identifies a spectrum of immune dysregulation in acquired bone marrow failure syndromes.

Poor graft function (PGF), manifested by multilineage cytopenias and complete donor chimerism post-allogeneic stem cell transplantation (alloSCT), and acquired aplastic anaemia (AA) are immune-mediated acquired bone marrow (BM) failure syndromes with a similar clinical presentation. In this study, we used spatial proteomics to compare the immunobiology of the BM microenvironment and identify common mechanisms of immune dysregulation under these conditions. Archival BM trephines from patients exhibited downregulation of the immunoregulatory protein VISTA and the M2 macrophage marker and suppressor of T-cell activation ARG1 with increased expression of the immune checkpoint B7-H3 compared to normal controls. Increased CD163 and CD14 expression suggested monocyte/macrophage skewing, which, combined with dysregulation of STING and VISTA, is indicative of an environment of reduced immunoregulation resulting in the profound suppression of hematopoiesis in these two conditions. There were no changes in the immune microenvironment between paired diagnostic AA and secondary MDS/AML samples suggesting that leukaemic clones develop in the impaired immune microenvironment of AA without the need for further alterations. Of the eight proteins with dysregulated expression shared by diagnostic AA and PGF, the diagnostic AA samples had a greater fold change in expression than PGF, suggesting that these diseases represent a spectrum of immune dysregulation. Unexpectedly, analysis of samples from patients with good graft function post-alloSCT demonstrated significant changes in the immune microenvironment compared to normal controls, with downregulation of CD44, STING, VISTA, and ARG1, suggesting that recovery of multilineage haematopoiesis post-alloSCT does not reflect recovery of immune function and may prime patients for the development of PGF upon further inflammatory insult. The demonstrable similarities in the immunopathology of AA and PGF will allow the design of clinical interventions that include both patient cohorts to accelerate therapeutic discovery and translation.

Discovery of a highly potent, selective, orally bioavailable inhibitor of KAT6A/B histone acetyltransferases with efficacy against KAT6A-high ER+ breast cancer.

KAT6A, and its paralog KAT6B, are histone lysine acetyltransferases (HAT) that acetylate histone H3K23 and exert an oncogenic role in several tumor types including breast cancer where KAT6A is frequently amplified/overexpressed. However, pharmacologic targeting of KAT6A to achieve therapeutic benefit has been a challenge. Here we describe identification of a highly potent, selective, and orally bioavailable KAT6A/KAT6B inhibitor CTx-648 (PF-9363), derived from a benzisoxazole series, which demonstrates anti-tumor activity in correlation with H3K23Ac inhibition in KAT6A over-expressing breast cancer. Transcriptional and epigenetic profiling studies show reduced RNA Pol II binding and downregulation of genes involved in estrogen signaling, cell cycle, Myc and stem cell pathways associated with CTx-648 anti-tumor activity in ER-positive (ER+) breast cancer. CTx-648 treatment leads to potent tumor growth inhibition in ER+ breast cancer in vivo models, including models refractory to endocrine therapy, highlighting the potential for targeting KAT6A in ER+ breast cancer.

Compartmentalized spatial profiling of the tumor microenvironment in head and neck squamous cell carcinoma identifies immune checkpoint molecules and tumor necrosis factor receptor superfamily members as biomarkers of response to immunotherapy.

Mucosal head and neck squamous cell carcinoma (HNSCC) are the seventh most common cancer, with approximately 50% of patients living beyond 5 years. Immune checkpoint inhibitors (ICIs) have shown promising results in patients with recurrent or metastatic (R/M) disease, however, only a subset of patients benefit from immunotherapy. Studies have implicated the tumor microenvironment (TME) of HNSCC as a major factor in therapy response, highlighting the need to better understand the TME, particularly by spatially resolved means to determine cellular and molecular components. Here, we employed targeted spatial profiling of proteins on a cohort of pre-treatment tissues from patients with R/M disease to identify novel biomarkers of response within the tumor and stromal margins. By grouping patient outcome categories into response or non-response, based on Response Evaluation Criteria in Solid Tumors (RECIST) we show that immune checkpoint molecules, including PD-L1, B7-H3, and VISTA, were differentially expressed. Patient responders possessed significantly higher tumor expression of PD-L1 and B7-H3, but lower expression of VISTA. Analysis of response subgroups indicated that tumor necrosis factor receptor (TNFR) superfamily members including OX40L, CD27, 4-1BB, CD40, and CD95/Fas, were associated with immunotherapy outcome. CD40 expression was higher in patient-responders than non responders, while CD95/Fas expression was lower in patients with partial response (PR) relative to those with stable disease (SD) and progressive disease (PD). Furthermore, we found that high 4-1BB expression in the tumor compartment, but not in the stroma, was associated with better overall survival (OS) (HR= 0.28, p-adjusted= 0.040). Moreover, high CD40 expression in tumor regions (HR= 0.27, p-adjusted= 0.035), and high CD27 expression in the stroma (HR= 0.2, p-adjusted=0.032) were associated with better survival outcomes. Taken together, this study supports the role of immune checkpoint molecules and implicates the TNFR superfamily as key players in immunotherapy response in our cohort of HNSCC. Validation of these findings in a prospective study is required to determine the robustness of these tissue signatures.

vissE.cloud: a webserver to visualise higher order molecular phenotypes from enrichment analysis.

Gene-set analysis (GSA) dominates the functional interpretation of omics data and downstream hypothesis generation. Despite its ability to summarise thousands of measurements into semantically interpretable components, GSA often results in hundreds of significantly enriched gene-sets. However, summarisation and effective visualisation of GSA results to facilitate hypothesis generation is still lacking. While some webservers provide gene-set visualization tools, there is still a need for tools that can effectively summarize and guide exploration of GSA results. To enable versatility, webservers accept gene lists as input, however, none provide end-to-end solutions for emerging data types such as single-cell and spatial omics. Here, we present vissE.Cloud, a webserver for end-to-end gene-set analysis, offering gene-set summarisation and highly interactive visualisation. vissE.Cloud uses algorithms from our earlier R package vissE to summarise GSA results by identifying biological themes. We maintain versatility by allowing analysis of gene lists, as well as, analysis of raw single-cell and spatial omics data, including CosMx and Xenium data, making vissE.Cloud the first webserver to provide end-to-end gene-set analysis of sub-cellular localised spatial data. Structuring the results hierarchically allows swift interactive investigations of results at the gene, gene-set, and clusters level. vissE.Cloud is freely available at https://www.vissE.Cloud.

Transcriptomic profiling of cardiac tissues from SARS-CoV-2 patients identifies DNA damage.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is known to present with pulmonary and extra-pulmonary organ complications. In comparison with the 2009 pandemic (pH1N1), SARS-CoV-2 infection is likely to lead to more severe disease, with multi-organ effects, including cardiovascular disease. SARS-CoV-2 has been associated with acute and long-term cardiovascular disease, but the molecular changes that govern this remain unknown. In this study, we investigated the host transcriptome landscape of cardiac tissues collected at rapid autopsy from seven SARS-CoV-2, two pH1N1, and six control patients using targeted spatial transcriptomics approaches. Although SARS-CoV-2 was not detected in cardiac tissue, host transcriptomics showed upregulation of genes associated with DNA damage and repair, heat shock, and M1-like macrophage infiltration in the cardiac tissues of COVID-19 patients. The DNA damage present in the SARS-CoV-2 patient samples, were further confirmed by γ-H2Ax immunohistochemistry. In comparison, pH1N1 showed upregulation of interferon-stimulated genes, in particular interferon and complement pathways, when compared with COVID-19 patients. These data demonstrate the emergence of distinct transcriptomic profiles in cardiac tissues of SARS-CoV-2 and pH1N1 influenza infection supporting the need for a greater understanding of the effects on extra-pulmonary organs, including the cardiovascular system of COVID-19 patients, to delineate the immunopathobiology of SARS-CoV-2 infection, and long term impact on health.

Profiling of lung SARS-CoV-2 and influenza virus infection dissects virus-specific host responses and gene signatures.

BACKGROUND: The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) which emerged in late 2019 has spread globally, causing a pandemic of respiratory illness designated coronavirus disease 2019 (COVID-19). A better definition of the pulmonary host response to SARS-CoV-2 infection is required to understand viral pathogenesis and to validate putative COVID-19 biomarkers that have been proposed in clinical studies. METHODS: Here, we use targeted transcriptomics of formalin-fixed paraffin-embedded tissue using the NanoString GeoMX platform to generate an in-depth picture of the pulmonary transcriptional landscape of COVID-19, pandemic H1N1 influenza and uninfected control patients. RESULTS: Host transcriptomics showed a significant upregulation of genes associated with inflammation, type I interferon production, coagulation and angiogenesis in the lungs of COVID-19 patients compared to non-infected controls. SARS-CoV-2 was non-uniformly distributed in lungs (emphasising the advantages of spatial transcriptomics) with the areas of high viral load associated with an increased type I interferon response. Once the dominant cell type present in the sample, within patient correlations and patient-patient variation, had been controlled for, only a very limited number of genes were differentially expressed between the lungs of fatal influenza and COVID-19 patients. Strikingly, the interferon-associated gene IFI27, previously identified as a useful blood biomarker to differentiate bacterial and viral lung infections, was significantly upregulated in the lungs of COVID-19 patients compared to patients with influenza. CONCLUSION: Collectively, these data demonstrate that spatial transcriptomics is a powerful tool to identify novel gene signatures within tissues, offering new insights into the pathogenesis of SARS-COV-2 to aid in patient triage and treatment.

IFI27 transcription is an early predictor for COVID-19 outcomes, a multi-cohort observational study.

Purpose: Robust biomarkers that predict disease outcomes amongst COVID-19 patients are necessary for both patient triage and resource prioritisation. Numerous candidate biomarkers have been proposed for COVID-19. However, at present, there is no consensus on the best diagnostic approach to predict outcomes in infected patients. Moreover, it is not clear whether such tools would apply to other potentially pandemic pathogens and therefore of use as stockpile for future pandemic preparedness. Methods: We conducted a multi-cohort observational study to investigate the biology and the prognostic role of interferon alpha-inducible protein 27 (IFI27) in COVID-19 patients. Results: We show that IFI27 is expressed in the respiratory tract of COVID-19 patients and elevated IFI27 expression in the lower respiratory tract is associated with the presence of a high viral load. We further demonstrate that the systemic host response, as measured by blood IFI27 expression, is associated with COVID-19 infection. For clinical outcome prediction (e.g., respiratory failure), IFI27 expression displays a high sensitivity (0.95) and specificity (0.83), outperforming other known predictors of COVID-19 outcomes. Furthermore, IFI27 is upregulated in the blood of infected patients in response to other respiratory viruses. For example, in the pandemic H1N1/09 influenza virus infection, IFI27-like genes were highly upregulated in the blood samples of severely infected patients. Conclusion: These data suggest that prognostic biomarkers targeting the family of IFI27 genes could potentially supplement conventional diagnostic tools in future virus pandemics, independent of whether such pandemics are caused by a coronavirus, an influenza virus or another as yet-to-be discovered respiratory virus.

Low-viscosity matrix suspension culture enables scalable analysis of patient-derived organoids and tumoroids from the large intestine.

Cell embedment into a solid support matrix is considered essential for the culture of intestinal epithelial organoids and tumoroids, but this technique presents challenges that impede scalable culture expansion, experimental manipulation, high-throughput screening and diagnostic applications. We have developed a low-viscosity matrix (LVM) suspension culture method that enables efficient establishment and propagation of organoids and tumoroids from the human large intestine. Organoids and tumoroids cultured in LVM suspension recapitulate the morphological development observed in solid matrices, with tumoroids reflecting the histological features and genetic heterogeneity of primary colorectal cancers. We demonstrate the utility of LVM suspension culture for organoid and tumoroid bioreactor applications and biobanking, as well as tumoroid high-throughput drug sensitivity testing. These methods provide opportunities for the study and use of patient-derived organoids and tumoroids from the large intestine.

Dissection of the bone marrow microenvironment in hairy cell leukaemia identifies prognostic tumour and immune related biomarkers.

Hairy cell leukaemia (HCL) is a rare CD20+ B cell malignancy characterised by rare "hairy" B cells and extensive bone marrow (BM) infiltration. Frontline treatment with the purine analogue cladribine (CDA) results in a highly variable response duration. We hypothesised that analysis of the BM tumour microenvironment would identify prognostic biomarkers of response to CDA. HCL BM immunology pre and post CDA treatment and healthy controls were analysed using Digital Spatial Profiling to assess the expression of 57 proteins using an immunology panel. A bioinformatics pipeline was developed to accommodate the more complex experimental design of a spatially resolved study. Treatment with CDA was associated with the reduction in expression of HCL tumour markers (CD20, CD11c) and increased expression of myeloid markers (CD14, CD68, CD66b, ARG1). Expression of HLA-DR, STING, CTLA4, VISTA, OX40L were dysregulated pre- and post-CDA. Duration of response to treatment was associated with greater reduction in tumour burden and infiltration by CD8 T cells into the BM post-CDA. This is the first study to provide a high multiplex analysis of HCL BM microenvironment demonstrating significant immune dysregulation and identify biomarkers of response to CDA. With validation in future studies, prospective application of these biomarkers could allow early identification and increased monitoring in patients at increased relapse risk post CDA.

SOX9 Defines Distinct Populations of Cells in SHH Medulloblastoma but Is Not Required for Math1-Driven Tumor Formation.

Medulloblastoma is the most common malignant pediatric brain tumor and there is an urgent need for molecularly targeted and subgroup-specific therapies. The stem cell factor SOX9, has been proposed as a potential therapeutic target for the treatment of Sonic Hedgehog medulloblastoma (SHH-MB) subgroup tumors, given its role as a downstream target of Hedgehog signaling and in functionally promoting SHH-MB metastasis and treatment resistance. However, the functional requirement for SOX9 in the genesis of medulloblastoma remains to be determined. Here we report a previously undocumented level of SOX9 expression exclusively in proliferating granule cell precursors (GCP) of the postnatal mouse cerebellum, which function as the medulloblastoma-initiating cells of SHH-MBs. Wild-type GCPs express comparatively lower levels of SOX9 than neural stem cells and mature astroglia and SOX9low GCP-like tumor cells constitute the bulk of both infant (Math1Cre:Ptch1lox/lox ) and adult (Ptch1LacZ/+ ) SHH-MB mouse models. Human medulloblastoma single-cell RNA data analyses reveal three distinct SOX9 populations present in SHH-MB and noticeably absent in other medulloblastoma subgroups: SOX9 + MATH1 + (GCP), SOX9 + GFAP + (astrocytes) and SOX9 + MATH1 + GFAP + (potential tumor-derived astrocytes). To functionally address whether SOX9 is required as a downstream effector of Hedgehog signaling in medulloblastoma tumor cells, we ablated Sox9 using a Math1Cre model system. Surprisingly, targeted ablation of Sox9 in GCPs (Math1Cre:Sox9lox/lox ) revealed no overt phenotype and loss of Sox9 in SHH-MB (Math1Cre:Ptch1lox/lox;Sox9lox/lox ) does not affect tumor formation. IMPLICATIONS: Despite preclinical data indicating SOX9 plays a key role in SHH-MB biology, our data argue against SOX9 as a viable therapeutic target.

Preclinical small molecule WEHI-7326 overcomes drug resistance and elicits response in patient-derived xenograft models of human treatment-refractory tumors.

Targeting cell division by chemotherapy is a highly effective strategy to treat a wide range of cancers. However, there are limitations of many standard-of-care chemotherapies: undesirable drug toxicity, side-effects, resistance and high cost. New small molecules which kill a wide range of cancer subtypes, with good therapeutic window in vivo, have the potential to complement the current arsenal of anti-cancer agents and deliver improved safety profiles for cancer patients. We describe results with a new anti-cancer small molecule, WEHI-7326, which causes cell cycle arrest in G2/M, cell death in vitro, and displays efficacious anti-tumor activity in vivo. WEHI-7326 induces cell death in a broad range of cancer cell lines, including taxane-resistant cells, and inhibits growth of human colon, brain, lung, prostate and breast tumors in mice xenografts. Importantly, the compound elicits tumor responses as a single agent in patient-derived xenografts of clinically aggressive, treatment-refractory neuroblastoma, breast, lung and ovarian cancer. In combination with standard-of-care, WEHI-7326 induces a remarkable complete response in a mouse model of high-risk neuroblastoma. WEHI-7326 is mechanistically distinct from known microtubule-targeting agents and blocks cells early in mitosis to inhibit cell division, ultimately leading to apoptotic cell death. The compound is simple to produce and possesses favorable pharmacokinetic and toxicity profiles in rodents. It represents a novel class of anti-cancer therapeutics with excellent potential for further development due to the ease of synthesis, simple formulation, moderate side effects and potent in vivo activity. WEHI-7326 has the potential to complement current frontline anti-cancer drugs and to overcome drug resistance in a wide range of cancers.

Modelling the effect of subcellular mutations on the migration of cells in the colorectal crypt.

BACKGROUND: Many cancers arise from mutations in cells within epithelial tissues. Mutations manifesting at the subcellular level influence the structure and function of the tissue resulting in cancer. Previous work has proposed how cell level properties can lead to mutant cell invasion, but has not incorporated detailed subcellular modelling RESULTS: We present a framework that allows the straightforward integration and simulation of SBML representations of subcellular dynamics within multiscale models of epithelial tissues. This allows us to investigate the effect of mutations in subcellular pathways on the migration of cells within the colorectal crypt. Using multiple models we find that mutations in APC, a key component in the Wnt signalling pathway, can bias neutral drift and can also cause downward invasion of mutant cells in the crypt. CONCLUSIONS: Our framework allows us to investigate how subcellular mutations, i.e. knockouts and knockdowns, affect cell-level properties and the resultant migration of cells within epithelial tissues. In the context of the colorectal crypt, we see that mutations in APC can lead directly to mutant cell invasion.

Imbalance of the renin-angiotensin system may contribute to inflammation and fibrosis in IBD: a novel therapeutic target?

OBJECTIVE: We evaluated the influence of the renin-angiotensin system (RAS) on intestinal inflammation and fibrosis. DESIGN: Cultured human colonic myofibroblast proliferation and collagen secretion were assessed following treatment with angiotensin (Ang) II and Ang (1-7), their receptor antagonists candesartan and A779, and the ACE inhibitor captopril. Circulating and intestinal RAS components were evaluated in patients with and without IBD. Disease outcomes in patients with IBD treated with ACE inhibitors and angiotensin receptor blockers (ARBs) were assessed in retrospective studies. RESULTS: Human colonic myofibroblast proliferation was reduced by Ang (1-7) in a dose-dependent manner (p<0.05). Ang II marginally but not significantly increased proliferation, an effect reversed by candesartan (p<0.001). Colonic myofibroblast collagen secretion was reduced by Ang (1-7) (p<0.05) and captopril (p<0.001), and was increased by Ang II (p<0.001). Patients with IBD had higher circulating renin (mean 25.4 vs 18.6 mIU/L, p=0.026) and ACE2:ACE ratio (mean 0.92 vs 0.69, p=0.015) than controls without IBD. RAS gene transcripts and peptides were identified in healthy and diseased bowels. Colonic mucosal Masson's trichrome staining correlated with Ang II (r=0.346, p=0.010) and inversely with ACE2 activity (r=-0.373, p=0.006). Patients with IBD who required surgery (1/37 vs 12/75, p=0.034) and hospitalisation (0/34 vs 8/68, p=0.049) over 2 years were less often treated with ACE inhibitors and ARBs than patients not requiring surgery or hospitalisation. CONCLUSIONS: The RAS mediates fibrosis in human cell cultures, is expressed in the intestine and perturbed in intestinal inflammation, and agents targeting this system are associated with improved disease outcomes.

Loss of Bcl-G, a Bcl-2 family member, augments the development of inflammation-associated colorectal cancer.

Gastrointestinal epithelial cells provide a selective barrier that segregates the host immune system from luminal microorganisms, thereby contributing directly to the regulation of homeostasis. We have shown that from early embryonic development Bcl-G, a Bcl-2 protein family member with unknown function, was highly expressed in gastrointestinal epithelial cells. While Bcl-G was dispensable for normal growth and development in mice, the loss of Bcl-G resulted in accelerated progression of colitis-associated cancer. A label-free quantitative proteomics approach revealed that Bcl-G may contribute to the stability of a mucin network, which when disrupted, is linked to colon tumorigenesis. Consistent with this, we observed a significant reduction in Bcl-G expression in human colorectal tumors. Our study identifies an unappreciated role for Bcl-G in colon cancer.