Stricturing Crohn's Disease Single-Cell RNA Sequencing Reveals Fibroblast Heterogeneity and Intercellular Interactions.

BACKGROUND & AIMS: Fibroblasts play a key role in stricture formation in Crohn's disease (CD) but understanding its pathogenesis requires a systems-level investigation to uncover new treatment targets. We studied full-thickness CD tissues to characterize fibroblast heterogeneity and function by generating the first single-cell RNA sequencing (scRNAseq) atlas of strictured bowel and providing proof of principle for therapeutic target validation. METHODS: We performed scRNAseq of 13 fresh full-thickness CD resections containing noninvolved, inflamed nonstrictured, and strictured segments as well as 7 normal non-CD bowel segments. Each segment was separated into mucosa/submucosa or muscularis propria and analyzed separately for a total of 99 tissue samples and 409,001 cells. We validated cadherin-11 (CDH11) as a potential therapeutic target by using whole tissues, isolated intestinal cells, NanoString nCounter, next-generation sequencing, proteomics, and animal models. RESULTS: Our integrated dataset revealed fibroblast heterogeneity in strictured CD with the majority of stricture-selective changes detected in the mucosa/submucosa, but not the muscle layer. Cell-cell interaction modeling revealed CXCL14+ as well as MMP/WNT5A+ fibroblasts displaying a central signaling role in CD strictures. CDH11, a fibroblast cell-cell adhesion molecule, was broadly expressed and up-regulated, and its profibrotic function was validated using NanoString nCounter, RNA sequencing, tissue target expression, in vitro gain- and loss-of-function experiments, proteomics, and knock-out and antibody-mediated CDH11 blockade in experimental colitis. CONCLUSIONS: A full-thickness bowel scRNAseq atlas revealed previously unrecognized fibroblast heterogeneity and interactions in CD strictures and CDH11 was validated as a potential therapeutic target. These results provide a new resource for a better understanding of CD stricture formation and open potential therapeutic developments. This work has been posted as a preprint on Biorxiv under doi: 10.1101/2023.04.03.534781.

Stricturing Crohn's disease single-cell RNA sequencing reveals fibroblast heterogeneity and intercellular interactions.

Background: Fibroblasts play a key role in stricture formation in Crohn's disease (CD) but understanding it's pathogenesis requires a systems-level investigation to uncover new treatment targets. We studied full thickness CD tissues to characterize fibroblast heterogeneity and function by generating the first single cell RNA sequencing (scRNAseq) atlas of strictured bowel and providing proof of principle for therapeutic target validation. Methods: We performed scRNAseq of 13 fresh full thickness CD resections containing non-involved, inflamed non-strictured, and strictured segments as well as 7 normal non-CD bowel segments. Each segment was separated into mucosa/submucosa or muscularis propria and analyzed separately for a total of 99 tissue samples and 409,001 cells. We validated cadherin-11 (CDH11) as a potential therapeutic target by using whole tissues, isolated intestinal cells, NanoString nCounter, next generation sequencing, proteomics and animal models. Results: Our integrated dataset revealed fibroblast heterogeneity in strictured CD with the majority of stricture-selective changes detected in the mucosa/submucosa, but not the muscle layer. Cell-cell interaction modeling revealed CXCL14+ as well as MMP/WNT5A+ fibroblasts displaying a central signaling role in CD strictures. CDH11, a fibroblast cell-cell adhesion molecule, was broadly expressed and upregulated, and its pro-fibrotic function was validated by NanoString nCounter, RNA sequencing, tissue target expression, in vitro gain- and loss-of-function experiments, proteomics, and two animal models of experimental colitis. Conclusion: A full-thickness bowel scRNAseq atlas revealed previously unrecognized fibroblast heterogeneity and interactions in CD strictures and CDH11 was validated as a potential therapeutic target. These results provide a new resource for a better understanding of CD stricture formation and opens potential therapeutic developments.

Epigenome erosion and SOX10 drive neural crest phenotypic mimicry in triple-negative breast cancer.

Intratumoral heterogeneity is caused by genomic instability and phenotypic plasticity, but how these features co-evolve remains unclear. SOX10 is a neural crest stem cell (NCSC) specifier and candidate mediator of phenotypic plasticity in cancer. We investigated its relevance in breast cancer by immunophenotyping 21 normal breast and 1860 tumour samples. Nuclear SOX10 was detected in normal mammary luminal progenitor cells, the histogenic origin of most TNBCs. In tumours, nuclear SOX10 was almost exclusive to TNBC, and predicted poorer outcome amongst cross-sectional (p = 0.0015, hazard ratio 2.02, n = 224) and metaplastic (p = 0.04, n = 66) cases. To understand SOX10's influence over the transcriptome during the transition from normal to malignant states, we performed a systems-level analysis of co-expression data, de-noising the networks with an eigen-decomposition method. This identified a core module in SOX10's normal mammary epithelial network that becomes rewired to NCSC genes in TNBC. Crucially, this reprogramming was proportional to genome-wide promoter methylation loss, particularly at lineage-specifying CpG-island shores. We propose that the progressive, genome-wide methylation loss in TNBC simulates more primitive epigenome architecture, making cells vulnerable to SOX10-driven reprogramming. This study demonstrates potential utility for SOX10 as a prognostic biomarker in TNBC and provides new insights about developmental phenotypic mimicry-a major contributor to intratumoral heterogeneity.

Chemotherapy induces dynamic immune responses in breast cancers that impact treatment outcome.

To elucidate the effects of neoadjuvant chemotherapy (NAC), we conduct whole transcriptome profiling coupled with histopathology analyses of a longitudinal breast cancer cohort of 146 patients including 110 pairs of serial tumor biopsies collected before treatment, after the first cycle of treatment and at the time of surgery. Here, we show that cytotoxic chemotherapies induce dynamic changes in the tumor immune microenvironment that vary by subtype and pathologic response. Just one cycle of treatment induces an immune stimulatory microenvironment harboring more tumor infiltrating lymphocytes (TILs) and up-regulation of inflammatory signatures predictive of response to anti-PD1 therapies while residual tumors are immune suppressed at end-of-treatment compared to the baseline. Increases in TILs and CD8+ T cell proportions in response to NAC are independently associated with pathologic complete response. Further, on-treatment immune response is more predictive of treatment outcome than immune features in paired baseline samples although these are strongly correlated.

LobSig is a multigene predictor of outcome in invasive lobular carcinoma.

Invasive lobular carcinoma (ILC) is the most common special type of breast cancer, and is characterized by functional loss of E-cadherin, resulting in cellular adhesion defects. ILC typically present as estrogen receptor positive, grade 2 breast cancers, with a good short-term prognosis. Several large-scale molecular profiling studies have now dissected the unique genomics of ILC. We have undertaken an integrative analysis of gene expression and DNA copy number to identify novel drivers and prognostic biomarkers, using in-house (n = 25), METABRIC (n = 125) and TCGA (n = 146) samples. Using in silico integrative analyses, a 194-gene set was derived that is highly prognostic in ILC (P = 1.20 × 10-5)-we named this metagene 'LobSig'. Assessing a 10-year follow-up period, LobSig outperformed the Nottingham Prognostic Index, PAM50 risk-of-recurrence (Prosigna), OncotypeDx, and Genomic Grade Index (MapQuantDx) in a stepwise, multivariate Cox proportional hazards model, particularly in grade 2 ILC cases (χ 2, P = 9.0 × 10-6), which are difficult to prognosticate clinically. Importantly, LobSig status predicted outcome with 94.6% accuracy amongst cases classified as 'moderate-risk' according to Nottingham Prognostic Index in the METABRIC cohort. Network analysis identified few candidate pathways, though genesets related to proliferation were identified, and a LobSig-high phenotype was associated with the TCGA proliferative subtype (χ 2, P < 8.86 × 10-4). ILC with a poor outcome as predicted by LobSig were enriched with mutations in ERBB2, ERBB3, TP53, AKT1 and ROS1. LobSig has the potential to be a clinically relevant prognostic signature and warrants further development.

Breast cancer metastasis to gynaecological organs: a clinico-pathological and molecular profiling study.

Breast cancer metastasis to gynaecological organs is an understudied pattern of tumour spread. We explored clinico-pathological and molecular features of these metastases to better understand whether this pattern of dissemination is organotropic or a consequence of wider metastatic dissemination. Primary and metastatic tumours from 54 breast cancer patients with gynaecological metastases were analysed using immunohistochemistry, DNA copy-number profiling, and targeted sequencing of 386 cancer-related genes. The median age of primary tumour diagnosis amongst patients with gynaecological metastases was significantly younger compared to a general breast cancer population (46.5 versus 60 years; p < 0.0001). Median age at metastatic diagnosis was 54.4, time to progression was 4.8 years (range 0-20 years), and survival following a diagnosis of metastasis was 1.95 years (range 0-18 years). Patients had an average of five involved sites (most frequently ovary, fallopian tube, omentum/peritoneum), with fewer instances of spread to the lungs, liver, or brain. Invasive lobular histology and luminal A-like phenotype were over-represented in this group (42.8 and 87.5%, respectively) and most patients had involved axillary lymph nodes (p < 0.001). Primary tumours frequently co-expressed oestrogen receptor cofactors (GATA3, FOXA1) and harboured amplifications at 8p12, 8q24, and 11q13. In terms of phenotype conversion, oestrogen receptor status was generally maintained in metastases, FOXA1 increased, and expression of progesterone receptor, androgen receptor, and GATA3 decreased. ESR1 and novel AR mutations were identified. Metastasis to gynaecological organs is a complication frequently affecting young women with invasive lobular carcinoma and luminal A-like breast cancer, and hence may be driven by sustained hormonal signalling. Molecular analyses reveal a spectrum of factors that could contribute to de novo or acquired resistance to therapy and disease progression.

Multi-omics profiling of younger Asian breast cancers reveals distinctive molecular signatures.

Breast cancer (BC) in the Asia Pacific regions is enriched in younger patients and rapidly rising in incidence yet its molecular bases remain poorly characterized. Here we analyze the whole exomes and transcriptomes of 187 primary tumors from a Korean BC cohort (SMC) enriched in pre-menopausal patients and perform systematic comparison with a primarily Caucasian and post-menopausal BC cohort (TCGA). SMC harbors higher proportions of HER2+ and Luminal B subtypes, lower proportion of Luminal A with decreased ESR1 expression compared to TCGA. We also observe increased mutation prevalence affecting BRCA1, BRCA2, and TP53 in SMC with an enrichment of a mutation signature linked to homologous recombination repair deficiency in TNBC. Finally, virtual microdissection and multivariate analyses reveal that Korean BC status is independently associated with increased TIL and decreased TGF-ß signaling expression signatures, suggesting that younger Asian BCs harbor more immune-active microenvironment than western BCs.

Molecular signatures in breast cancer.

The use of molecular signatures to add value to standard clinical and pathological parameters has impacted clinical practice in many cancer types, but perhaps most notably in the breast cancer field. This is, in part, due to the considerable complexity of the disease at the clinical, morphological and molecular levels. The adoption of molecular profiling of DNA, RNA and protein continues to reveal important differences in the intrinsic biology between molecular subtypes and has begun to impact the way patients are managed. Several bioinformatic tools have been developed using DNA or RNA-based signatures to stratify the disease into biologically and/or clinically meaningful subgroups. Here, we review the approaches that have been used to develop gene expression signatures into currently available diagnostic assays (e.g., OncotypeDX® and Mammaprint®), plus we describe the latest work on genome sequencing, the methodologies used in the discovery process of mutational signatures, and the potential of these signatures to impact the clinic.

Understanding the functional impact of copy number alterations in breast cancer using a network modeling approach.

UNLABELLED: Copy number alterations (CNAs) are thought to account for 85% of the variation in gene expression observed among breast tumours. The expression of cis-associated genes is impacted by CNAs occurring at proximal loci of these genes, whereas the expression of trans-associated genes is impacted by CNAs occurring at distal loci. While a majority of these CNA-driven genes responsible for breast tumourigenesis are cis-associated, trans-associated genes are thought to further abet the development of cancer and influence disease outcomes in patients. Here we present a network-based approach that integrates copy-number and expression profiles to identify putative cis- and trans-associated genes in breast cancer pathogenesis. We validate these cis- and trans-associated genes by employing them to subtype a large cohort of breast tumours obtained from the METABRIC consortium, and demonstrate that these genes accurately reconstruct the ten subtypes of breast cancer. We observe that individual breast cancer subtypes are driven by distinct sets of cis- and trans-associated genes. Among the cis-associated genes, we recover several known drivers of breast cancer (e.g. CCND1, ERRB2, MDM2 and ZNF703) and some novel putative drivers (e.g. BRF2 and SF3B3). siRNA-mediated knockdown of BRF2 across a panel of breast cancer cell lines showed significant reduction in cell viability for ER-/HER2+ (MDA-MB-453) cells, but not in normal (MCF10A) cells thereby indicating that BRF2 could be a viable therapeutic target for estrogen receptor-negative/HER2-enriched (ER-/HER2+) cancers. Among the trans-associated genes, we identify modules of immune response (CD2, CD19, CD38 and CD79B), mitotic/cell-cycle kinases (e.g. AURKB, MELK, PLK1 and TTK), and DNA-damage response genes (e.g. RFC4 and FEN1). siRNA-mediated knockdown of RFC4 significantly reduced cell proliferation in ER-negative normal breast and cancer lines, thereby indicating that RFC4 is essential for both normal and cancer cell survival but could be a useful biomarker for aggressive (ER-negative) breast tumours. AVAILABILITY: under NetStrat.

Personalised pathway analysis reveals association between DNA repair pathway dysregulation and chromosomal instability in sporadic breast cancer.

The Homologous Recombination (HR) pathway is crucial for the repair of DNA double-strand breaks (DSBs) generated during DNA replication. Defects in HR repair have been linked to the initiation and development of a wide variety of human malignancies, and exploited in chemical, radiological and targeted therapies. In this study, we performed a personalised pathway analysis independently for four large sporadic breast cancer cohorts to investigate the status of HR pathway dysregulation in individual sporadic breast tumours, its association with HR repair deficiency and its impact on tumour characteristics. Specifically, we first manually curated a list of HR genes according to our recent review on this pathway (Liu et al., 2014), and then applied a personalised pathway analysis method named Pathifier (Drier et al., 2013) on the expression levels of the curated genes to obtain an HR score quantifying HR pathway dysregulation in individual tumours. Based on the score, we observed a great diversity in HR dysregulation between and within gene expression-based breast cancer subtypes, and by using two published HR-defect signatures, we found HR pathway dysregulation reflects HR repair deficiency. Furthermore, we identified a novel association between HR pathway dysregulation and chromosomal instability (CIN) in sporadic breast cancer. Although CIN has long been considered as a hallmark of most solid tumours, with recent extensive studies highlighting its importance in tumour evolution and drug resistance, the molecular basis of CIN in sporadic cancers remains poorly understood. Our results imply that HR pathway dysregulation might contribute to CIN in sporadic breast cancer.