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.

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.

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.

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.

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.

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.

TGF-ß and IL-6 family signalling crosstalk: an integrated model.

Mathematical models for TGF-ß and IL-6 signalling have been linked, providing a platform for analyzing the crosstalk between the systems. An integrated IL-6:TGF-ß model was developed via a reduced set of reaction equations which incorporate both feedback loops and appropriate time-delays for transcription and translation processes. The model simulates stable, robust and realistic responses to both ligands. Pulsatile (multiple pulses) inputs for both TGF-ß and IL-6 have been simulated to investigate the effects of each ligand on the sensitivity, equilibrium and dynamic responses of the integrated signalling system. In our simulations the crosstalk between constant IL-6 and TGF-ß signalling via SMAD7 does not appear to be sufficient to render the cells resistant to TGF-ß inhibition. However, the simulations predict that pulsatile IL-6 stimulation would increase SMAD7 levels substantially and consequentially, lead to resistance to TGF-ß. The model also allows the prediction of the integrated signalling pathway responses to the mutation of key components, e.g. Gp130 F/F.

The nuage mediates retrotransposon silencing in mouse primordial ovarian follicles.

Mobilization of endogenous retrotransposons can destabilize the genome, an imminent danger during epigenetic reprogramming of cells in the germline. The P-element-induced wimpy testis (PIWI)-interacting RNA (piRNA) pathway is known to silence retrotransposons in the mouse testes. Several piRNA pathway components localize to the unique, germline structure known as the nuage. In this study, we surveyed mouse ovaries and found, for the first time, transient appearance of nuage-like structures in oocytes of primordial follicles. Mouse vasa homolog (MVH), Piwi-like 2 (PIWIL2/MILI) and tudor domain-containing 9 (TDRD9) are present in these structures, whereas aggregates of germ cell protein with ankyrin repeats, sterile alpha motif and leucine zipper (GASZ) localize separately in the cytoplasm. Retrotransposons are silenced in primordial ovarian follicles, and de-repressed upon reduction of piRNA expression in Mvh, Mili or Gasz mutants. However, these null-mutant females, unlike their male counterparts, are fertile, uncoupling retrotransposon activation from sterility.

Discrete Element Framework for Modelling Extracellular Matrix, Deformable Cells and Subcellular Components.

This paper presents a framework for modelling biological tissues based on discrete particles. Cell components (e.g. cell membranes, cell cytoskeleton, cell nucleus) and extracellular matrix (e.g. collagen) are represented using collections of particles. Simple particle to particle interaction laws are used to simulate and control complex physical interaction types (e.g. cell-cell adhesion via cadherins, integrin basement membrane attachment, cytoskeletal mechanical properties). Particles may be given the capacity to change their properties and behaviours in response to changes in the cellular microenvironment (e.g., in response to cell-cell signalling or mechanical loadings). Each particle is in effect an 'agent', meaning that the agent can sense local environmental information and respond according to pre-determined or stochastic events. The behaviour of the proposed framework is exemplified through several biological problems of ongoing interest. These examples illustrate how the modelling framework allows enormous flexibility for representing the mechanical behaviour of different tissues, and we argue this is a more intuitive approach than perhaps offered by traditional continuum methods. Because of this flexibility, we believe the discrete modelling framework provides an avenue for biologists and bioengineers to explore the behaviour of tissue systems in a computational laboratory.

Colonic myofibroblast cell line stimulates colonoid formation.

A stable and efficient system for the culture of murine colon epithelial cells or crypts is required to facilitate studies of the dynamics and factors affecting colon stem cell niche and crypt formation. Survival of colonic epithelial cells or crypts in vitro was not established until recently, when it was found that exogenous Wnt3A and R-spondin could promote cell survival and formation of spheroids (colonospheres) or some advanced organoids with well-developed crypts (colonoids). However, after 6-8 days in these culture conditions, only small numbers of colonospheres form organoids with crypt-like structures (colonoids). This study describes the use of a myofibroblast cell line and a coculture system that increases the efficiency of colonoid formation from isolated crypts. The enhanced coculture system has significantly improved colonoid-forming efficiency compared with results from previous systems. Crypt formation can be detected as early as day 2. The coculture system will facilitate the characterization of the colon stem cell niche and the changes that occur as a result of perturbations or mutations in colon stem or epithelial cells, such as those that favor precancerous adenoma or cancer.

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.

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.

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.

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.

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.

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.