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.

Spatial resolution of the head and neck cancer tumor microenvironment to identify tumor and stromal features associated with therapy response.

Head and neck cancer (HNC) is the seventh most common cancer globally, resulting in 440 000 deaths per year. While there have been advancements in chemoradiotherapy and surgery, relapse occurs in more than half of HNCs, and these patients have a median survival of 10 months and a 2-year survival of < 20%. Only a subset of patients displays durable benefits from immunotherapies in metastatic and recurrent HNC, making it critical to understand the tumor microenvironment (TME) underpinning therapy responses in HNC. To recognize biological differences within the TME that may be predictive of immunotherapy response, we applied cutting-edge geospatial whole-transcriptome profiling (NanoString GeoMx Digital Spatial Profiler) and spatial proteomics profiling (Akoya PhenoCycler-Fusion) on a tumor microarray consisting of 25 cores from 12 patients that included 4 immunotherapy-unresponsive (8 cores) and 2 immunotherapy-responsive patients (5 cores), as well as 6 immunotherapy naïve patients (12 cores). Through high-plex, regional-based transcriptomic mapping of the tumor and TME, pathways involved with the complement system and hypoxia were identified to be differentially expressed in patients who went on to experience a poor immunotherapy response. Single-cell, targeted proteomic analysis found that immune cell infiltration of the cancer cell mass and interactions of CD8 T cells with tumor and other immune cells were associated with positive immunotherapy response. The relative abundance of specific tumor phenotypes and their interactions with various immune cells was identified to be different between response groups. This study demonstrates how spatial transcriptomics and proteomics can resolve novel alterations in the TME of HNC that may contribute to therapy sensitivity and resistance.

Spatial dynamics of tertiary lymphoid aggregates in head and neck cancer: insights into immunotherapy response.

BACKGROUND: Recurrent/metastatic head and neck squamous cell carcinoma (R/M HNSCC) generally has a poor prognosis for patients with limited treatment options. While incorporating immune checkpoint inhibitors (ICIs) has now become the standard of care, the efficacy is variable, with only a subset of patients responding. The complexity of the tumor microenvironment (TME) and the role of tertiary lymphoid structures (TLS) have emerged as critical determinants for immunotherapeutic response. METHODS: In this study, we analyzed two independently collected R/M HNSCC patient tissue cohorts to better understand the role of TLS in response to ICIs. Utilizing a multi-omics approach, we first performed targeted proteomic profiling using the Nanostring GeoMx Digital Spatial Profiler to quantify immune-related protein expression with spatial resolution. This was further characterized by spatially resolved whole transcriptome profiling of TLSs and germinal centers (GCs). Deeper single-cell resolved proteomic profiling of the TLSs was performed using the Akoya Biosciences Phenocycler Fusion platform. RESULTS: Our proteomic analysis revealed the presence of T lymphocyte markers, including CD3, CD45, and CD8, expressing cells and upregulation of immune checkpoint marker PD-L1 within tumor compartments of patients responsive to ICIs, indicative of 'hot tumor' phenotypes. We also observed the presence of antigen-presenting cells marked by expression of CD40, CD68, CD11c, and CD163 with upregulation of antigen-presentation marker HLA-DR, in patients responding to ICIs. Transcriptome analysis of TLS and GCs uncovered a marked elevation in the expression of genes related to immune modulation, diverse immune cell recruitment, and a potent interferon response within the TLS structure. Notably, the distribution of TLS-tumor distance was found to be significantly different across response groups (H = 9.28, p = 0.026). The proximity of TLSs to tumor cells was found to be a critical indicator of ICI response, implying that patients with TLSs located further from tumor cells have worse outcomes. CONCLUSION: The study underscores the multifaceted role of TLSs in modulating the immunogenic landscape of the TME in R/M HNSCC, likely influencing the efficacy of ICIs. Spatially resolved multi-omics approaches offer valuable insights into potential biomarkers for ICI response and highlight the importance of profiling the TME complexity when developing therapeutic strategies and patient stratification.

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.

Transient inhibition of type I interferon enhances CD8 + T cell stemness and vaccine protection.

Developing vaccines that promote CD8 + T cell memory is a challenge for infectious disease and cancer immunotherapy. TCF-1 + stem cell-like memory T (T SCM ) cells are important determinants of long-lived memory. Yet, the developmental requirements for T SCM formation are unclear. Here, we identify the temporal window for type I interferon (IFN-I) receptor (IFNAR) blockade to drive T SCM cell generation. T SCM cells were transcriptionally distinct and emerged from a transitional precursor of exhausted (T PEX ) cellular state concomitant with viral clearance. T SCM differentiation correlated with T cell retention within the lymph node paracortex, due to increased CXCR3 chemokine abundance which disrupted gradient formation. These affects were due a counterintuitive increase in IFNψ, which controlled cell location. Combining IFNAR inhibition with mRNA-LNP vaccination promoted specific T SCM differentiation and enhanced protection against chronic infection. These finding propose a new approach to vaccine design whereby modulation of inflammation promotes memory formation and function. HIGHLIGHTS: Early, transient inhibition of the type I interferon (IFN) receptor (IFNAR) during acute viral infection promotes stem cell-like memory T (T SCM ) cell differentiation without establishing chronic infection. T SCM and precursor of exhausted (T PEX ) cellular states are distinguished transcriptionally and by cell surface markers. Developmentally, T SCM cell differentiation occurs via a transition from a T PEX state coinciding with viral clearance. Transient IFNAR blockade increases IFNψ production to modulate the ligands of CXCR3 and couple T SCM differentiation to cell retention within the T cell paracortex of the lymph node. Specific promotion of T SCM cell differentiation with nucleoside-modified mRNA-LNP vaccination elicits enhanced protection against chronic viral challenge.

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.

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 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.

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.