OCA-B promotes autoimmune demyelination through control of stem-like CD4 + T cells

Stem-like T cell populations can selectively contribute to autoimmunity, but the activities that promote and sustain these populations are incompletely understood. Here, we show that T cell-intrinsic loss of the transcription cofactor OCA-B protects mice from experimental autoimmune encephalomyelitis (EAE) while preserving responses to CNS infection. In adoptive transfer EAE models driven by multiple antigen encounters, OCA-B deletion nearly eliminates CNS infiltration, proinflammatory cytokine production and clinical disease. OCA-B-expressing CD4 + T cells within the CNS of mice with EAE comprise a minority of the population but display a memory phenotype and preferentially confer disease. In a relapsing-remitting EAE model, OCA-B T cell deficiency specifically protects mice from relapse. During remission, OCA-B promotes the expression of Tcf7 , Slamf6 , and Sell in proliferating T cell populations. At relapse, OCA-B loss results in both the accumulation of an immunomodulatory CD4 + T cell population expressing Ccr9 and Bach2 , and the loss of pro-inflammatory gene expression from Th17 cells. These results identify OCA-B as a driver of pathogenic stem-like T cells.

Mcam stabilizes luminal progenitor breast cancer phenotypes via Ck2 control and Src/Akt/Stat3 attenuation

Breast cancers are categorized into subtypes with distinctive therapeutic vulnerabilities and prognoses based on their expression of clinically targetable receptors and gene expression patterns mimicking different cell types of the normal gland. Here, we tested the role of Mcam in breast cancer cell state control and tumorigenicity in a luminal progenitor-like murine tumor cell line (Py230) that exhibits lineage and tumor subtype plasticity. Mcam knockdown Py230 cells show augmented Stat3 and Pi3K/Akt activation associated with a lineage state switch away from a hormone-sensing/luminal progenitor state toward alveolar and basal cell related phenotypes that were refractory to growth inhibition by the anti-estrogen therapeutic, tamoxifen. Inhibition of Stat3, or the upstream activator Ck2, reversed these cell state changes. Mcam binds Ck2 and acts as a regulator of Ck2 substrate utilization across multiple mammary tumor cell lines. In Py230 cells this activity manifests as increased mesenchymal morphology, migration, and Src/Fak/Mapk/Paxillin adhesion complex signaling in vitro, in contrast to Mcam's reported roles in promoting mesenchymal phenotypes. In vivo, Mcam knockdown reduced tumor growth and take rate and inhibited cell state transition to Sox10+/neural crest like cells previously been associated with tumor aggressiveness. This contrasts with human luminal breast cancers where MCAM copy number loss is highly coupled to Cyclin D amplification, increased proliferation, and the more aggressive Luminal B subtype. Together these data indicate a critical role for Mcam and its regulation of Ck2 in control of breast cancer cell state plasticity with implications for progression, evasion of targeted therapies and combination therapy design.

Abstract P3-09-16: Validating alternative splicing events between bulk and scRNA sequencing data: A bioinformatic approach

Background: As bulk and single-cell RNA (scRNA) sequencing studies and data continue to accrue to publicly accessible databases, bioinformatic tools continue being developed to analyze these datasets. The resulting pipelines often focus on either bulk or scRNA sequencing but rarely both. scRNA sequencing allows for in-depth characterization of transcription level events in distinct cell populations while bulk data offers a more global view, and for now, is more widely available in databases. Alternative splicing events in particular have been underexplored in a single-cell genomic architecture. The purpose of this study was to develop a bioinformatic pipeline to employ both bulk and single-cell mammary datasets to identify and validate alternative splicing events.Methods: Bulk fastq’s were utilized to identify AS events in luminal progenitor (LP) and mouse mammary stem cell (MSC) lineages in ovariectomized FVB mice that had been randomized into three treatment groups (SHAM [C], estradiol + progesterone [EP], EP + the selective progesterone receptor inhibitor telapristone acetate, TPA [EPT]). scRNA sequencing data from (Bach et al., 2017) was downloaded from the Gene Expression Omnibus (GEO) and utilized for validation of the AS events identified in the bulk data. To separate the data by progesterone level, the scRNA data was parsed into nulliparous, gestational, lactational, and post-involution groups. Further, to determine cell lineage level effects, the data was parsed into luminal and basal compartments using krt18 and krt5 expression. To accomplish the overall comparison between bulk and scRNA sequencing data, we developed a pipeline that would process and parse the scRNA data, identify the alternative splicing in bulk sequencing using rMATS and in the scRNA data using Outrigger (part of the Expedition suite). The resulting data from rMATS [bulk] LP/MSC: CvEP, CvEPT, & EPvEPT were filtered for events meeting statistical significance (p-value < 0.05, FDR < 0.01). After filtering, the events from each treatment comparison were then filtered so only events unique to LP or MSC CvEP, respectively, remain. The remaining significant and unique AS events are then compared to the Krt18-high and Krt5-high Outrigger results from each development stage based on genomic event coordinates (with a buffer +/- 20 base pairs).Results: 12 alternative splicing (AS) events were identified to be shared between the bulk and scRNA sequencing data in the context of increased progesterone exposure (EP and gestational, respectively) as well as cell lineage (LP and luminal respectively). Among the genes identified as having been alternatively spiced with an exact nucleotide match for the splice were Eif4a2, Pik3c2a, Brd4, Cdh1, and Enah. Conclusions: Alternative splicing events can be identified and validated between genomic sequencing methods, specifically bulk and single-cell RNA seq. We hypothesize that as scRNA-seq becomes more developed and quantification errors caused by short-read lengths fixed by full length scRNA seq methods currently in development, comparing data between single-cell and bulk may be made easier. However, we have shown that, when orthogonal methods for validation may not be feasible, events as specific as alternative splices can be validated using publicly available data and in-silico methods.

Citation Format: Gannon Cottone, Benjamin T Spike, Elnaz Mirzaei Mehrabad, Seema A Khan, Susan Clare. Validating alternative splicing events between bulk and scRNA sequencing data: A bioinformatic approach [abstract]. In: Proceedings of the 2021 San Antonio Breast Cancer Symposium; 2021 Dec 7-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2022;82(4 Suppl):Abstract nr P3-09-16.

Progesterone receptor antagonists reverse stem cell expansion and the paracrine effectors of progesterone action in the mouse mammary gland

Background

The ovarian hormones estrogen and progesterone (EP) are implicated in breast cancer causation. A specific consequence of progesterone exposure is the expansion of the mammary stem cell (MSC) and luminal progenitor (LP) compartments. We hypothesized that this effect, and its molecular facilitators, could be abrogated by progesterone receptor (PR) antagonists administered in a mouse model.

Methods

Ovariectomized FVB mice were randomized to 14 days of treatment: sham, EP, EP + telapristone (EP + TPA), EP + mifepristone (EP + MFP). Mice were then sacrificed, mammary glands harvested, and mammary epithelial cell lineages separated by flow cytometry using cell surface markers. RNA from each lineage was sequenced and differential gene expression was analyzed using DESeq. Quantitative PCR was performed to confirm the candidate genes discovered in RNA seq. ANOVA with Tukey post hoc analysis was performed to compare relative expression. Alternative splicing events were examined using the rMATs multivariate analysis tool.

Results

Significant increases in the MSC and luminal mature (LM) cell fractions were observed following EP treatment compared to control (p < 0.01 and p < 0.05, respectively), whereas the LP fraction was significantly reduced (p < 0.05). These hormone-induced effects were reversed upon exposure to TPA and MFP (p < 0.01 for both). Gene Ontology analysis of RNA-sequencing data showed EP-induced enrichment of several pathways, with the largest effect on Wnt signaling in MSC, significantly repressed by PR inhibitors. In LP cells, significant induction of Wnt4 and Rankl, and Wnt pathway intermediates Lrp2 and Axin2 (confirmed by qRTPCR) were reversed by TPA and MFP (p < 0.0001). Downstream signaling intermediates of these pathways (Lrp5, Mmp7) showed similar effects. Expression of markers of epithelial-mesenchymal transition (Cdh1, Cdh3) and the induction of EMT regulators (Zeb1, Zeb2, Gli3, Snai1, and Ptch2) were significantly responsive to progesterone. EP treatment was associated with large-scale alternative splicing events, with an enrichment of motifs associated with Srsf, Esrp, and Rbfox families. Exon skipping was observed in Cdh1, Enah, and Brd4.

Conclusions

PR inhibition reverses known tumorigenic pathways in the mammary gland and suppresses a previously unknown effect of progesterone on RNA splicing events. In total, our results strengthen the case for reconsideration of PR inhibitors for breast cancer prevention.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13058-021-01455-2.

An NKX2-1/ERK/WNT feedback loop modulates gastric identity and response to targeted therapy in lung adenocarcinoma

Cancer cells undergo lineage switching during natural progression and in response to therapy. NKX2-1 loss in human and murine lung adenocarcinoma leads to invasive mucinous adenocarcinoma (IMA), a lung cancer subtype that exhibits gastric differentiation and harbors a distinct spectrum of driver oncogenes. In murine BRAFV600E-driven lung adenocarcinoma, NKX2-1 is required for early tumorigenesis, but dispensable for established tumor growth. NKX2-1-deficient, BRAFV600E-driven tumors resemble human IMA and exhibit a distinct response to BRAF/MEK inhibitors. Whereas BRAF/MEK inhibitors drive NKX2-1-positive tumor cells into quiescence, NKX2-1-negative cells fail to exit the cell cycle after the same therapy. BRAF/MEK inhibitors induce cell identity switching in NKX2-1-negative lung tumors within the gastric lineage, which is driven in part by WNT signaling and FoxA1/2. These data elucidate a complex, reciprocal relationship between lineage specifiers and oncogenic signaling pathways in the regulation of lung adenocarcinoma identity that is likely to impact lineage-specific therapeutic strategies.

CRIPTO antagonist ALK4L75A-Fc inhibits breast cancer cell plasticity and adaptation to stress

BACKGROUND

CRIPTO is a multi-functional signaling protein that promotes stemness and oncogenesis. We previously developed a CRIPTO antagonist, ALK4L75A-Fc, and showed that it causes loss of the stem cell phenotype in normal mammary epithelia suggesting it may similarly inhibit CRIPTO-dependent plasticity in breast cancer cells.

METHODS

We focused on two triple negative breast cancer cell lines (MDA-MB-231 and MDA-MB-468) to measure the effects of ALK4L75A-Fc on cancer cell behavior under nutrient deprivation and endoplasmic reticulum stress. We characterized the proliferation and migration of these cells in vitro using time-lapse microscopy and characterized stress-dependent changes in the levels and distribution of CRIPTO signaling mediators and cancer stem cell markers. We also assessed the effects of ALK4L75A-Fc on proliferation, EMT, and stem cell markers in vivo as well as on tumor growth and metastasis using inducible lentiviral delivery or systemic administration of purified ALK4L75A-Fc, which represents a candidate therapeutic approach.

RESULTS

ALK4L75A-Fc inhibited adaptive responses of breast cancer cells under conditions of nutrient and ER stress and reduced their proliferation, migration, clonogenicity, and expression of EMT and cancer stem cell markers. ALK4L75A-Fc also inhibited proliferation of human breast cancer cells in stressed tumor microenvironments in xenografts and reduced both primary tumor size and metastatic burden.

CONCLUSIONS

Cancer cell adaptation to stresses such as nutrient deprivation, hypoxia, and chemotherapy can critically contribute to dormancy, metastasis, therapy resistance, and recurrence. Identifying mechanisms that govern cellular adaptation, plasticity, and the emergence of stem-like cancer cells may be key to effective anticancer therapies. Results presented here indicate that targeting CRIPTO with ALK4L75A-Fc may have potential as such a therapy since it inhibits breast cancer cell adaptation to microenvironmental challenges and associated stem-like and EMT phenotypes.

Abstract ES10-2: Understanding breast cancer using a developmental perspective

Parallels among embryonic development, stem cells, and cancer have long been recognized. We identified, isolated, and characterized stem cells that first become committed to a mammary fate during embryogenesis; we refer to these cells as fetal mammary stem cells (fMaSCs). Lineage tracing, in vitro sphere formation, and in vivo transplantation studies by our group and many others all confirm that cells in the embryo are the bipotent progenitors of the mammary gland. There is debate, however, on whether such bipotent cells persist into the adult, or whether the luminal and basal lineages are maintained by unipotent progenitors. To gain insight into the relationships between fMaSCs and breast cancer, and to investigate their potential persistence in the adult, we have applied bulk and single cell RNA-sequencing (sc-RNA-seq) and single nucleus ATAC-sequencing (snATAC-seq) throughout mammary development. The results to be discussed demonstrate that fMaSC transcriptomes are heterogeneous, but all share co-expression of genes associated with luminal and basal cell fates. This fits a model in which the bipotent state is created by a balance of lineage specifiers. We also find that the fMaSC transcriptome is highly enriched in basal-like human breast cancers and identify potential embryonic pathways that correlate with poor prognosis. We used a variety of computational tools to infer the gene expression programs that ensue when fMaSCs commit to luminal and basal states. The data from scRNA-seq and snATAC-seq demonstrate that the transitions are gradual, not precipitous, and that luminal and basal cells exhibit significant transcriptomic and epigenetic heterogeneity. This challenges the notion that the mammary gland consists of discrete cell types defined by rigid transcriptomic parameters, and reveals a potential for intrinsic phenotypic plasticity of normal mammary cells. Using the combined databases, we identified Sox10 as a significantly differentially expressed cell state regulator. We show that tumors are heterogeneous with regard to Sox10 expression, and that locally invasive cells tend to express high Sox10 levels. Elevated Sox10 correlates with acquisition of a neural-crest like, EMT-related state. Implications for interception of metastasis by targeting neural crest-like cells will be discussed. Finally, we have generated a web resource that is available to the scientific community to enable the transcription and epigenetic characteristics of any gene of interest to be tracked through mammary development (https://wahl-labsalk.shinyapps.io/Mammary_snATAC/).

Citation Format: GM Wahl, Z Ma, C Chung, C Dravis, BT Spike, RR Giraddi, O Balcioglu, C Fan, B Hagos, R Heinz, Herrera-Valdez J, X Hou, J Hwang, R Lasken, G Luna, NE Lytle, EM Mehrabad, M Novotny, CM Perou, O Poirion, S Preissl, B Ren, T Reya, CL Trejo, KT Varley. Understanding breast cancer using a developmental perspective [abstract]. In: Proceedings of the 2019 San Antonio Breast Cancer Symposium; 2019 Dec 10-14; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2020;80(4 Suppl):Abstract nr ES10-2.

Single-Cell Transcriptomes Distinguish Stem Cell State Changes and Lineage Specification Programs in Early Mammary Gland Development

The mammary gland consists of cells with gene expression patterns reflecting their cellular origins, function, and spatiotemporal context. However, knowledge of developmental kinetics and mechanisms of lineage specification is lacking. We address this significant knowledge gap by generating a single-cell transcriptome atlas encompassing embryonic, postnatal, and adult mouse mammary development. From these data, we map the chronology of transcriptionally and epigenetically distinct cell states and distinguish fetal mammary stem cells (fMaSCs) from their precursors and progeny. fMaSCs show balanced co-expression of factors associated with discrete adult lineages and a metabolic gene signature that subsides during maturation but reemerges in some human breast cancers and metastases. These data provide a useful resource for illuminating mammary cell heterogeneity, the kinetics of differentiation, and developmental correlates of tumorigenesis.

Single-cell RNA sequencing of developing mouse mammary epithelia reveals the timing of lineage specification. Giraddi et al. find that fetal mammary stem cells co-express factors that define distinct lineages in their progeny and bear functionally relevant metabolic program signatures that change with differentiation and are resurrected in human breast cancers and metastases.

PMLPR: A novel method for predicting subcellular localization based on recommender systems

The importance of protein subcellular localization problem is due to the importance of protein's functions in different cell parts. Moreover, prediction of subcellular locations helps to identify the potential molecular targets for drugs and has an important role in genome annotation. Most of the existing prediction methods assign only one location for each protein. But, since some proteins move between different subcellular locations, they can have multiple locations. In recent years, some multiple location predictors have been introduced. However, their performances are not accurate enough and there is much room for improvement. In this paper, we introduced a method, PMLPR, to predict locations for a protein. PMLPR predicts a list of locations for each protein based on recommender systems and it can properly overcome the multiple location prediction problem. For evaluating the performance of PMLPR, we considered six datasets RAT, FLY, HUMAN, Du et al., DBMLoc and Höglund. The performance of this algorithm is compared with six state-of-the-art algorithms, YLoc, WOLF-PSORT, prediction channel, MDLoc, Du et al. and MultiLoc2-HighRes. The results indicate that our proposed method is significantly superior on RAT and Fly proteins, and decent on HUMAN proteins. Moreover, on the datasets introduced by Du et al., DBMLoc and Höglund, PMLPR has comparable results. For the case study, we applied the algorithms on 8 proteins which are important in cancer research. The results of comparison with other methods indicate the efficiency of PMLPR.