Two regulatory T cell populations in the visceral adipose tissue shape systemic metabolism.

Visceral adipose tissue (VAT) is an energy store and endocrine organ critical for metabolic homeostasis. Regulatory T (Treg) cells restrain inflammation to preserve VAT homeostasis and glucose tolerance. Here, we show that the VAT harbors two distinct Treg cell populations: prototypical serum stimulation 2-positive (ST2+) Treg cells that are enriched in males and a previously uncharacterized population of C-X-C motif chemokine receptor 3-positive (CXCR3+) Treg cells that are enriched in females. We show that the transcription factors GATA-binding protein 3 and peroxisome proliferator-activated receptor-γ, together with the cytokine interleukin-33, promote the differentiation of ST2+ VAT Treg cells but repress CXCR3+ Treg cells. Conversely, the differentiation of CXCR3+ Treg cells is mediated by the cytokine interferon-γ and the transcription factor T-bet, which also antagonize ST2+ Treg cells. Finally, we demonstrate that ST2+ Treg cells preserve glucose homeostasis, whereas CXCR3+ Treg cells restrain inflammation in lean VAT and prevent glucose intolerance under high-fat diet conditions. Overall, this study defines two molecularly and developmentally distinct VAT Treg cell types with unique context- and sex-specific functions.

A single-cell RNA expression atlas of normal, preneoplastic and tumorigenic states in the human breast.

To examine global changes in breast heterogeneity across different states, we determined the single-cell transcriptomes of > 340,000 cells encompassing normal breast, preneoplastic BRCA1+/- tissue, the major breast cancer subtypes, and pairs of tumors and involved lymph nodes. Elucidation of the normal breast microenvironment revealed striking changes in the stroma of post-menopausal women. Single-cell profiling of 34 treatment-naive primary tumors, including estrogen receptor (ER)+ , HER2+ , and triple-negative breast cancers, revealed comparable diversity among cancer cells and a discrete subset of cycling cells. The transcriptomes of preneoplastic BRCA1+/- tissue versus tumors highlighted global changes in the immune microenvironment. Within the tumor immune landscape, proliferative CD8+ T cells characterized triple-negative and HER2+ cancers but not ER+ tumors, while all subtypes comprised cycling tumor-associated macrophages, thus invoking potentially different immunotherapy targets. Copy number analysis of paired ER+ tumors and lymph nodes indicated seeding by genetically distinct clones or mass migration of primary tumor cells into axillary lymph nodes. This large-scale integration of patient samples provides a high-resolution map of cell diversity in normal and cancerous human breast.

cellCounts: an R function for quantifying 10x Chromium single-cell RNA sequencing data.

SUMMARY: The 10x Genomics Chromium single-cell RNA sequencing technology is a powerful gene expression profiling platform, which is capable of profiling expression of thousands of genes in tens of thousands of cells simultaneously. This platform can produce hundreds of million reads in a single experiment, making it a very challenging task to quantify expression of genes in individual cells due to the massive data volume. Here, we present cellCounts, a new tool for efficient and accurate quantification of Chromium data. cellCounts employs the seed-and-vote strategy to align reads to a reference genome, collapses reads to Unique Molecular Identifiers (UMIs) and then assigns UMIs to genes based on the featureCounts program. Using both simulation and real datasets for evaluation, cellCounts was found to compare favourably to cellRanger and STARsolo. cellCounts is implemented in R, making it easily integrated with other R programs for analysing Chromium data. AVAILABILITY AND IMPLEMENTATION: cellCounts was implemented as a function in R package Rsubread that can be downloaded from http://bioconductor.org/packages/release/bioc/html/Rsubread.html. Data and analysis code used in this study can be freely accessed via La Trobe University's Institutional Repository at https://doi.org/10.26181/21588276.

Single cell transcriptome atlas of mouse mammary epithelial cells across development.

BACKGROUND: Heterogeneity within the mouse mammary epithelium and potential lineage relationships have been recently explored by single-cell RNA profiling. To further understand how cellular diversity changes during mammary ontogeny, we profiled single cells from nine different developmental stages spanning late embryogenesis, early postnatal, prepuberty, adult, mid-pregnancy, late-pregnancy, and post-involution, as well as the transcriptomes of micro-dissected terminal end buds (TEBs) and subtending ducts during puberty. METHODS: The single cell transcriptomes of 132,599 mammary epithelial cells from 9 different developmental stages were determined on the 10x Genomics Chromium platform, and integrative analyses were performed to compare specific time points. RESULTS: The mammary rudiment at E18.5 closely aligned with the basal lineage, while prepubertal epithelial cells exhibited lineage segregation but to a less differentiated state than their adult counterparts. Comparison of micro-dissected TEBs versus ducts showed that luminal cells within TEBs harbored intermediate expression profiles. Ductal basal cells exhibited increased chromatin accessibility of luminal genes compared to their TEB counterparts suggesting that lineage-specific chromatin is established within the subtending ducts during puberty. An integrative analysis of five stages spanning the pregnancy cycle revealed distinct stage-specific profiles and the presence of cycling basal, mixed-lineage, and 'late' alveolar intermediates in pregnancy. Moreover, a number of intermediates were uncovered along the basal-luminal progenitor cell axis, suggesting a continuum of alveolar-restricted progenitor states. CONCLUSIONS: This extended single cell transcriptome atlas of mouse mammary epithelial cells provides the most complete coverage for mammary epithelial cells during morphogenesis to date. Together with chromatin accessibility analysis of TEB structures, it represents a valuable framework for understanding developmental decisions within the mouse mammary gland.

Transcriptome sequencing and multi-plex imaging of prostate cancer microenvironment reveals a dominant role for monocytic cells in progression.

BACKGROUND: Prostate cancer is caused by genomic aberrations in normal epithelial cells, however clinical translation of findings from analyses of cancer cells alone has been very limited. A deeper understanding of the tumour microenvironment is needed to identify the key drivers of disease progression and reveal novel therapeutic opportunities. RESULTS: In this study, the experimental enrichment of selected cell-types, the development of a Bayesian inference model for continuous differential transcript abundance, and multiplex immunohistochemistry permitted us to define the transcriptional landscape of the prostate cancer microenvironment along the disease progression axis. An important role of monocytes and macrophages in prostate cancer progression and disease recurrence was uncovered, supported by both transcriptional landscape findings and by differential tissue composition analyses. These findings were corroborated and validated by spatial analyses at the single-cell level using multiplex immunohistochemistry. CONCLUSIONS: This study advances our knowledge concerning the role of monocyte-derived recruitment in primary prostate cancer, and supports their key role in disease progression, patient survival and prostate microenvironment immune modulation.

Canonical PRC2 function is essential for mammary gland development and affects chromatin compaction in mammary organoids.

Distinct transcriptional states are maintained through organization of chromatin, resulting from the sum of numerous repressive and active histone modifications, into tightly packaged heterochromatin versus more accessible euchromatin. Polycomb repressive complex 2 (PRC2) is the main mammalian complex responsible for histone 3 lysine 27 trimethylation (H3K27me3) and is integral to chromatin organization. Using in vitro and in vivo studies, we show that deletion of Suz12, a core component of all PRC2 complexes, results in loss of H3K27me3 and H3K27 dimethylation (H3K27me2), completely blocks normal mammary gland development, and profoundly curtails progenitor activity in 3D organoid cultures. Through the application of mammary organoids to bypass the severe phenotype associated with Suz12 loss in vivo, we have explored gene expression and chromatin structure in wild-type and Suz12-deleted basal-derived organoids. Analysis of organoids led to the identification of lineage-specific changes in gene expression and chromatin structure, inferring cell type-specific PRC2-mediated gene silencing of the chromatin state. These expression changes were accompanied by cell cycle arrest but not lineage infidelity. Together, these data indicate that canonical PRC2 function is essential for development of the mammary gland through the repression of alternate transcription programs and maintenance of chromatin states.

Control of mammary stem cell function by steroid hormone signalling.

The ovarian hormones oestrogen and progesterone profoundly influence breast cancer risk, underpinning the benefit of endocrine therapies in the treatment of breast cancer. Modulation of their effects through ovarian ablation or chemoprevention strategies also significantly decreases breast cancer incidence. Conversely, there is an increased risk of breast cancer associated with pregnancy in the short term. The cellular mechanisms underlying these observations, however, are poorly defined. Here we demonstrate that mouse mammary stem cells (MaSCs) are highly responsive to steroid hormone signalling, despite lacking the oestrogen and progesterone receptors. Ovariectomy markedly diminished MaSC number and outgrowth potential in vivo, whereas MaSC activity increased in mice treated with oestrogen plus progesterone. Notably, even three weeks of treatment with the aromatase inhibitor letrozole was sufficient to reduce the MaSC pool. In contrast, pregnancy led to a transient 11-fold increase in MaSC numbers, probably mediated through paracrine signalling from RANK ligand. The augmented MaSC pool indicates a cellular basis for the short-term increase in breast cancer incidence that accompanies pregnancy. These findings further indicate that breast cancer chemoprevention may be achieved, in part, through suppression of MaSC function.

Integration of microRNA signatures of distinct mammary epithelial cell types with their gene expression and epigenetic portraits.

INTRODUCTION: MicroRNAs (miRNAs) have been implicated in governing lineage specification and differentiation in multiple organs; however, little is known about their specific roles in mammopoiesis. We have determined the global miRNA expression profiles of functionally distinct epithelial subpopulations in mouse and human mammary tissue, and compared these to their cognate transcriptomes and epigenomes. Finally, the human miRNA signatures were used to interrogate the different subtypes of breast cancer, with a view to determining miRNA networks deregulated during oncogenesis. METHODS: RNA from sorted mouse and human mammary cell subpopulations was subjected to miRNA expression analysis using the TaqMan MicroRNA Array. Differentially expressed (DE) miRNAs were correlated with gene expression and histone methylation profiles. Analysis of miRNA signatures of the intrinsic subtypes of breast cancer in The Cancer Genome Atlas (TCGA) database versus those of normal human epithelial subpopulations was performed. RESULTS: Unique miRNA signatures characterized each subset (mammary stem cell (MaSC)/basal, luminal progenitor, mature luminal, stromal), with a high degree of conservation across species. Comparison of miRNA and transcriptome profiles for the epithelial subtypes revealed an inverse relationship and pinpointed key developmental genes. Interestingly, expression of the primate-specific miRNA cluster (19q13.4) was found to be restricted to the MaSC/basal subset. Comparative analysis of miRNA signatures with H3 lysine modification maps of the different epithelial subsets revealed a tight correlation between active or repressive marks for the top DE miRNAs, including derepression of miRNAs in Ezh2-deficient cellular subsets. Interrogation of TCGA-identified miRNA profiles with the miRNA signatures of different human subsets revealed specific relationships. CONCLUSIONS: The derivation of global miRNA expression profiles for the different mammary subpopulations provides a comprehensive resource for understanding the interplay between miRNA networks and target gene expression. These data have highlighted lineage-specific miRNAs and potential miRNA-mRNA networks, some of which are disrupted in neoplasia. Furthermore, our findings suggest that key developmental miRNAs are regulated by global changes in histone modification, thus linking the mammary epigenome with genome-wide changes in the expression of genes and miRNAs. Comparative miRNA signature analyses between normal breast epithelial cells and breast tumors confirmed an important linkage between luminal progenitor cells and basal-like tumors.

PRMT1 acts as a suppressor of MHC-I and anti-tumor immunity.

Cancer immunotherapies have demonstrated remarkable success; however, the majority of patients do not respond or develop resistance. Here, we conduct epigenetic gene-targeted CRISPR-Cas9 screens to identify epigenomic factors that limit CD8+ T cell-mediated anti-tumor immunity. We identify that PRMT1 suppresses interferon gamma (Ifnγ)-induced MHC-I expression, thus dampening CD8+ T cell-mediated killing. Indeed, PRMT1 knockout or pharmacological targeting of type I PRMT with the clinical inhibitor GSK3368715 enhances Ifnγ-induced MHC-I expression through elevated STAT1 expression and activation, while re-introduction of PRMT1 in PRMT1-deficient cells reverses this effect. Importantly, loss of PRMT1 enhances the efficacy of anti-PD-1 immunotherapy, and The Cancer Genome Atlas analysis reveals that PRMT1 expression in human melanoma is inversely correlated with expression of human leukocyte antigen molecules, infiltration of CD8+ T cells, and overall survival. Taken together, we identify PRMT1 as a negative regulator of anti-tumor immunity, unveiling clinical type I PRMT inhibitors as immunotherapeutic agents or as adjuncts to existing immunotherapies.

Mechanisms of cellular crosstalk in the gastric tumor microenvironment are mediated by YAP1 and STAT3.

Deregulation of the Hippo pathway is a driver for cancer progression and treatment resistance. In the context of gastric cancer, YAP1 is a biomarker for poor patient prognosis. Although genomic tumor profiling provides information of Hippo pathway activation, the present study demonstrates that inhibition of Yap1 activity has anti-tumor effects in gastric tumors driven by oncogenic mutations and inflammatory cytokines. We show that Yap1 is a key regulator of cell metabolism, proliferation, and immune responses in normal and neoplastic gastric epithelium. We propose that the Hippo pathway is targetable across gastric cancer subtypes and its therapeutic benefits are likely to be mediated by both cancer cell-intrinsic and -extrinsic mechanisms.

Ehf controls mammary alveolar lineage differentiation and is a putative suppressor of breast tumorigenesis.

The transcription factor EHF is highly expressed in the lactating mammary gland, but its role in mammary development and tumorigenesis is not fully understood. Utilizing a mouse model of Ehf deletion, herein, we demonstrate that loss of Ehf impairs mammary lobuloalveolar differentiation at late pregnancy, indicated by significantly reduced levels of milk genes and milk lipids, fewer differentiated alveolar cells, and an accumulation of alveolar progenitor cells. Further, deletion of Ehf increased proliferative capacity and attenuated prolactin-induced alveolar differentiation in mammary organoids. Ehf deletion also increased tumor incidence in the MMTV-PyMT mammary tumor model and increased the proliferative capacity of mammary tumor organoids, while low EHF expression was associated with higher tumor grade and poorer outcome in luminal A and basal human breast cancers. Collectively, these findings establish EHF as a non-redundant regulator of mammary alveolar differentiation and a putative suppressor of mammary tumorigenesis.

FOXC1 is enriched in the mammary luminal progenitor population, but is not necessary for mouse mammary ductal morphogenesis.

Expression of FOXC1, a forkhead box transcription factor, correlates with the human basal-like breast cancer (BLBC) subtype, and functional analyses have revealed its importance for in vitro invasiveness of BLBC cells. Women diagnosed with this breast tumor subtype have a poorer outcome because of the lack of targeted therapies; thus, continued investigation of factors driving these tumors is critical to uncover novel therapeutic targets. Several processes that dictate normal mammary morphogenesis parallel cancer progression, and enforced expression of FOXC1 can induce a progenitor state in more-differentiated mammary epithelial cells. Consequently, evaluating how FOXC1 functions in the normal gland is critical to further understand BLBC biology. Although FOXC1 is well known to control normal development of a number of tissues, its role in the mammary gland has not yet been investigated. Herein, we describe FOXC1 expression patterning in the normal breast, where it is localized to the basal/myoepithelium, suggesting that FOXC1 would be required for normal development. However, mammary glands lacking Foxc1 have no overt defect in ductal outgrowth, alveologenesis, or lineage specification. Of significant interest, we found that expression of FOXC1 is enriched in the normal luminal progenitor population, which is the postulated cell of origin of BLBC. These results indicate that FOXC1 is unnecessary for mammary morphogenesis and that its role in BLBC likely involves processes that are unrelated to cell lineage specification.

Single-cell RNA sequencing captures patient-level heterogeneity and associated molecular phenotypes in breast cancer pleural effusions.

BACKGROUND: Malignant pleural effusions (MPEs) are a common complication of advanced cancers, particularly those adjacent to the pleura, such as lung and breast cancer. The pathophysiology of MPE formation remains poorly understood, and although MPEs are routinely used for the diagnosis of breast cancer patients, their composition and biology are poorly understood. It is difficult to distinguish invading malignant cells from resident mesothelial cells and to identify the directionality of interactions between these populations in the pleura. There is a need to characterize the phenotypic diversity of breast cancer cell populations in the pleural microenvironment, and investigate how this varies across patients. METHODS: Here, we used single-cell RNA-sequencing to study the heterogeneity of 10 MPEs from seven metastatic breast cancer patients, including three Miltenyi-enriched samples using a negative selection approach. This dataset of almost 65 000 cells was analysed using integrative approaches to compare heterogeneous cell populations and phenotypes. RESULTS: We identified substantial inter-patient heterogeneity in the composition of cell types (including malignant, mesothelial and immune cell populations), in expression of subtype-specific gene signatures and in copy number aberration patterns, that captured variability across breast cancer cell populations. Within individual MPEs, we distinguished mesothelial cell populations from malignant cells using key markers, the presence of breast cancer subtype expression patterns and copy number aberration patterns. We also identified pleural mesothelial cells expressing a cancer-associated fibroblast-like transcriptomic program that may support cancer growth. CONCLUSIONS: Our dataset presents the first unbiased assessment of breast cancer-associated MPEs at a single cell resolution, providing the community with a valuable resource for the study of MPEs. Our work highlights the molecular and cellular diversity captured in MPEs and motivates the potential use of these clinically relevant biopsies in the development of targeted therapeutics for patients with advanced breast cancer.

Induction of endoplasmic reticulum stress is associated with the anti-tumor activity of monepantel across cancer types.

BACKGROUND: Monepantel is an anti-helminthic drug that also has anti-cancer properties. Despite several studies over the years, the molecular target of monepantel in mammalian cells is still unknown, and its mechanism-of-action is not fully understood, though effects on cell cycle, mTOR signalling and autophagy have been implicated. METHODS: Viability assays were performed on >20 solid cancer cell cells, and apoptosis assays were performed on a subset of these, including 3D cultures. Genetic deletion of BAX/BAK and ATG were used to establish roles of apoptosis and autophagy in killing activity. RNA-sequencing was performed on four cell lines after monepantel treatment, and differentially regulated genes were confirmed by Western blotting. RESULTS: We showed that monepantel has anti-proliferative activity on a broad range of cancer cell lines. In some, this was associated with induction of apoptosis which was confirmed using a BAX/BAK-deficient cell line. However, proliferation is still inhibited in these cells following monepantel treatment, indicating cell-cycle disruption as the major anti-cancer effect. Previous studies have also indicated autophagic cell death occurs following monepantel treatment. We showed autophagy induction in multiple cell lines; however, deletion of a key autophagy regulator ATG7 had minimal impact on monepantel's anti-proliferative activity, suggesting autophagy is associated with, but not required for its anti-tumour effects. Transcriptomic analysis of four cell lines treated with monepantel revealed downregulation of many genes involved in the cell cycle, and upregulation of genes linked to ATF4-mediated ER stress responses, especially those involved in amino-acid metabolism and protein synthesis. CONCLUSIONS: As these outcomes are all associated with mTOR signalling, cell cycle and autophagy, we now provide a likely triggering mechanism for the anti-cancer activity of monepantel.

A tuft cell - ILC2 signaling circuit provides therapeutic targets to inhibit gastric metaplasia and tumor development.

Although gastric cancer is a leading cause of cancer-related deaths, systemic treatment strategies remain scarce. Here, we report the pro-tumorigenic properties of the crosstalk between intestinal tuft cells and type 2 innate lymphoid cells (ILC2) that is evolutionarily optimized for epithelial remodeling in response to helminth infection. We demonstrate that tuft cell-derived interleukin 25 (IL25) drives ILC2 activation, inducing the release of IL13 and promoting epithelial tuft cell hyperplasia. While the resulting tuft cell - ILC2 feed-forward circuit promotes gastric metaplasia and tumor formation, genetic depletion of tuft cells or ILC2s, or therapeutic targeting of IL13 or IL25 alleviates these pathologies in mice. In gastric cancer patients, tuft cell and ILC2 gene signatures predict worsening survival in intestinal-type gastric cancer where ~40% of the corresponding cancers show enriched co-existence of tuft cells and ILC2s. Our findings suggest a role for ILC2 and tuft cells, along with their associated cytokine IL13 and IL25 as gatekeepers and enablers of metaplastic transformation and gastric tumorigenesis, thereby providing an opportunity to therapeutically inhibit early-stage gastric cancer through repurposing antibody-mediated therapies.

TCF-1 limits intraepithelial lymphocyte antitumor immunity in colorectal carcinoma.

Intraepithelial lymphocytes (IELs), including αß and γδ T cells (T-IELs), constantly survey and play a critical role in maintaining the gastrointestinal epithelium. We show that cytotoxic molecules important for defense against cancer were highly expressed by T-IELs in the small intestine. In contrast, abundance of colonic T-IELs was dependent on the microbiome and displayed higher expression of TCF-1/TCF7 and a reduced effector and cytotoxic profile, including low expression of granzymes. Targeted deletion of TCF-1 in γδ T-IELs induced a distinct effector profile and reduced colon tumor formation in mice. In addition, TCF-1 expression was significantly reduced in γδ T-IELs present in human colorectal cancers (CRCs) compared with normal healthy colon, which strongly correlated with an enhanced γδ T-IEL effector phenotype and improved patient survival. Our work identifies TCF-1 as a colon-specific T-IEL transcriptional regulator that could inform new immunotherapy strategies to treat CRC.

Silicone finger prosthesis. A clinical report.

The loss of all or part of a finger following traumatic amputation may have a negative impact on physical and psychological well being. An esthetic prosthesis can offer psychological, functional, and rehabilitative advantages. The success of a prosthetic restoration primarily depends on its retention. This clinical report describes an alternative method of retention by scoring the master cast of a partially amputated finger, thus enhancing the vacuum effect for the retention of the prosthesis. The methodology of treatment is also explained. Silicone material was used to provide function and esthetics.

R code and downstream analysis objects for the scRNA-seq atlas of normal and tumorigenic human breast tissue.

Breast cancer is a common and highly heterogeneous disease. Understanding cellular diversity in the mammary gland and its surrounding micro-environment across different states can provide insight into cancer development in the human breast. Recently, we published a large-scale single-cell RNA expression atlas of the human breast spanning normal, preneoplastic and tumorigenic states. Single-cell expression profiles of nearly 430,000 cells were obtained from 69 distinct surgical tissue specimens from 55 patients. This article extends the study by providing quality filtering thresholds, downstream processed R data objects, complete cell annotation and R code to reproduce all the analyses. Data quality assessment measures are presented and details are provided for all the bioinformatic analyses that produced results described in the study.

Computational Screening of Anti-Cancer Drugs Identifies a New BRCA Independent Gene Expression Signature to Predict Breast Cancer Sensitivity to Cisplatin.

The development of therapies that target specific disease subtypes has dramatically improved outcomes for patients with breast cancer. However, survival gains have not been uniform across patients, even within a given molecular subtype. Large collections of publicly available drug screening data matched with transcriptomic measurements have facilitated the development of computational models that predict response to therapy. Here, we generated a series of predictive gene signatures to estimate the sensitivity of breast cancer samples to 90 drugs, comprising FDA-approved drugs or compounds in early development. To achieve this, we used a cell line-based drug screen with matched transcriptomic data to derive in silico models that we validated in large independent datasets obtained from cell lines and patient-derived xenograft (PDX) models. Robust computational signatures were obtained for 28 drugs and used to predict drug efficacy in a set of PDX models. We found that our signature for cisplatin can be used to identify tumors that are likely to respond to this drug, even in absence of the BRCA-1 mutation routinely used to select patients for platinum-based therapies. This clinically relevant observation was confirmed in multiple PDXs. Our study foreshadows an effective delivery approach for precision medicine.

A diverse fibroblastic stromal cell landscape in the spleen directs tissue homeostasis and immunity.

The spleen is a compartmentalized organ that serves as a blood filter and safeguard of systemic immune surveillance. Labyrinthine networks of fibroblastic stromal cells construct complex niches within the white pulp and red pulp that are important for tissue homeostasis and immune activation. However, the identity and roles of the global splenic fibroblastic stromal cells in homeostasis and immune responses are poorly defined. Here, we performed a cellular and molecular dissection of the splenic reticular stromal cell landscape. We found that white pulp fibroblastic reticular cells (FRCs) responded robustly during acute viral infection, but this program of gene regulation was suppressed during persistent viral infection. Single-cell transcriptomic analyses in mice revealed diverse fibroblast cell niches and unexpected heterogeneity among podoplanin-expressing cells that include glial, mesothelial, and adventitial cells in addition to FRCs. We found analogous fibroblastic stromal cell diversity in the human spleen. In addition, we identify the transcription factor SpiB as a critical regulator required to support white pulp FRC differentiation, homeostatic chemokine expression, and antiviral T cell responses. Together, our study provides a comprehensive map of fibroblastic stromal cell types in the spleen and defines roles for red and white pulp fibroblasts for splenic function and orchestration of immune responses.