Human NOP2/NSUN1 regulates ribosome biogenesis through non-catalytic complex formation with box C/D snoRNPs.

5-Methylcytosine (m5C) is a base modification broadly found on various RNAs in the human transcriptome. In eukaryotes, m5C is catalyzed by enzymes of the NSUN family composed of seven human members (NSUN1-7). NOP2/NSUN1 has been primarily characterized in budding yeast as an essential ribosome biogenesis factor required for the deposition of m5C on the 25S ribosomal RNA (rRNA). Although human NOP2/NSUN1 has been known to be an oncogene overexpressed in several types of cancer, its functions and substrates remain poorly characterized. Here, we used a miCLIP-seq approach to identify human NOP2/NSUN1 RNA substrates. Our analysis revealed that NOP2/NSUN1 catalyzes the deposition of m5C at position 4447 on the 28S rRNA. We also find that NOP2/NSUN1 binds to the 5'ETS region of the pre-rRNA transcript and regulates pre-rRNA processing through non-catalytic complex formation with box C/D snoRNAs. We provide evidence that NOP2/NSUN1 facilitates the recruitment of U3 and U8 snoRNAs to pre-90S ribosomal particles and their stable assembly into snoRNP complexes. Remarkably, expression of both WT and catalytically inactive NOP2/NSUN1 in knockdown background rescues the rRNA processing defects and the stable assembly of box C/D snoRNP complexes, suggesting that NOP2/NSUN1-mediated deposition of m5C on rRNA is not required for ribosome synthesis.

Knowledge-based classification of fine-grained immune cell types in single-cell RNA-Seq data.

Single-cell RNA sequencing (scRNA-Seq) is an emerging strategy for characterizing immune cell populations. Compared to flow or mass cytometry, scRNA-Seq could potentially identify cell types and activation states that lack precise cell surface markers. However, scRNA-Seq is currently limited due to the need to manually classify each immune cell from its transcriptional profile. While recently developed algorithms accurately annotate coarse cell types (e.g. T cells versus macrophages), making fine distinctions (e.g. CD8+ effector memory T cells) remains a difficult challenge. To address this, we developed a machine learning classifier called ImmClassifier that leverages a hierarchical ontology of cell type. We demonstrate that its predictions are highly concordant with flow-based markers from CITE-seq and outperforms other tools (+15% recall, +14% precision) in distinguishing fine-grained cell types with comparable performance on coarse ones. Thus, ImmClassifier can be used to explore more deeply the heterogeneity of the immune system in scRNA-Seq experiments.

Circulating immune cell phenotype dynamics reflect the strength of tumor-immune cell interactions in patients during immunotherapy.

The extent to which immune cell phenotypes in the peripheral blood reflect within-tumor immune activity prior to and early in cancer therapy is unclear. To address this question, we studied the population dynamics of tumor and immune cells, and immune phenotypic changes, using clinical tumor and immune cell measurements and single-cell genomic analyses. These samples were serially obtained from a cohort of advanced gastrointestinal cancer patients enrolled in a trial with chemotherapy and immunotherapy. Using an ecological population model, fitted to clinical tumor burden and immune cell abundance data from each patient, we find evidence of a strong tumor-circulating immune cell interaction in responder patients but not in those patients that progress on treatment. Upon initiation of therapy, immune cell abundance increased rapidly in responsive patients, and once the peak level is reached tumor burden decreases, similar to models of predator-prey interactions; these dynamic patterns were absent in nonresponder patients. To interrogate phenotype dynamics of circulating immune cells, we performed single-cell RNA sequencing at serial time points during treatment. These data show that peripheral immune cell phenotypes were linked to the increased strength of patients' tumor-immune cell interaction, including increased cytotoxic differentiation and strong activation of interferon signaling in peripheral T cells in responder patients. Joint modeling of clinical and genomic data highlights the interactions between tumor and immune cell populations and reveals how variation in patient responsiveness can be explained by differences in peripheral immune cell signaling and differentiation soon after the initiation of immunotherapy.

Kras mutation rate precisely orchestrates ductal derived pancreatic intraepithelial neoplasia and pancreatic cancer.

Pancreatic ductal adenocarcinoma (PDAC) is the third leading cause of cancer-related death in the United States. Despite the high prevalence of Kras mutations in pancreatic cancer patients, murine models expressing the oncogenic mutant Kras (Krasmut) in mature pancreatic cells develop PDAC at a low frequency. Independent of cell of origin, a second genetic hit (loss of tumor suppressor TP53 or PTEN) is important for development of PDAC in mice. We hypothesized ectopic expression and elevated levels of oncogenic mutant Kras would promote PanIN arising in pancreatic ducts. To test our hypothesis, the significance of elevating levels of K-Ras and Ras activity has been explored by expression of a CAG driven LGSL-KrasG12V allele (cKras) in pancreatic ducts, which promotes ectopic Kras expression. We predicted expression of cKras in pancreatic ducts would generate neoplasia and PDAC. To test our hypothesis, we employed tamoxifen dependent CreERT2 mediated recombination. Hnf1b:CreERT2;KrasG12V (cKrasHnf1b/+) mice received 1 (Low), 5 (Mod) or 10 (High) mg per 20 g body weight to recombine cKras in low (cKrasLow), moderate (cKrasMod), and high (cKrasHigh) percentages of pancreatic ducts. Our histologic analysis revealed poorly differentiated aggressive tumors in cKrasHigh mice. cKrasMod mice had grades of Pancreatic Intraepithelial Neoplasia (PanIN), recapitulating early and advanced PanIN observed in human PDAC. Proteomics analysis revealed significant differences in PTEN/AKT and MAPK pathways between wild type, cKrasLow, cKrasMod, and cKrasHigh mice. In conclusion, in this study, we provide evidence that ectopic expression of oncogenic mutant K-Ras in pancreatic ducts generates early and late PanIN as well as PDAC. This Ras rheostat model provides evidence that AKT signaling is an important early driver of invasive ductal derived PDAC.

EMTome: a resource for pan-cancer analysis of epithelial-mesenchymal transition genes and signatures.

BACKGROUND: The epithelial-mesenchymal transition (EMT) enables dissociation of tumour cells from the primary tumour mass, invasion through the extracellular matrix, intravasation into blood vessels and colonisation of distant organs. Cells that revert to the epithelial state via the mesenchymal-epithelial transition cause metastases, the primary cause of death in cancer patients. EMT also empowers cancer cells with stem-cell properties and induces resistance to chemotherapeutic drugs. Understanding the driving factors of EMT is critical for the development of effective therapeutic interventions. METHODS: This manuscript describes the generation of a database containing EMT gene signatures derived from cell lines, patient-derived xenografts and patient studies across cancer types and multiomics data and the creation of a web-based portal to provide a comprehensive analysis resource. RESULTS: EMTome incorporates (i) EMT gene signatures; (ii) EMT-related genes with multiomics features across different cancer types; (iii) interactomes of EMT-related genes (miRNAs, transcription factors, and proteins); (iv) immune profiles identified from The Cancer Genome Atlas (TCGA) cohorts by exploring transcriptomics, epigenomics, and proteomics, and drug sensitivity and (iv) clinical outcomes of cancer cohorts linked to EMT gene signatures. CONCLUSION: The web-based EMTome portal is a resource for primary and metastatic tumour research publicly available at www.emtome.org .

CD8+ T cells inhibit metastasis and CXCL4 regulates its function.

BACKGROUND: The mechanism by which immune cells regulate metastasis is unclear. Understanding the role of immune cells in metastasis will guide the development of treatments improving patient survival. METHODS: We used syngeneic orthotopic mouse tumour models (wild-type, NOD/scid and Nude), employed knockout (CD8 and CD4) models and administered CXCL4. Tumours and lungs were analysed for cancer cells by bioluminescence, and circulating tumour cells were isolated from blood. Immunohistochemistry on the mouse tumours was performed to confirm cell type, and on a tissue microarray with 180 TNBCs for human relevance. TCGA data from over 10,000 patients were analysed as well. RESULTS: We reveal that intratumoral immune infiltration differs between metastatic and non-metastatic tumours. The non-metastatic tumours harbour high levels of CD8+ T cells and low levels of platelets, which is reverse in metastatic tumours. During tumour progression, platelets and CXCL4 induce differentiation of monocytes into myeloid-derived suppressor cells (MDSCs), which inhibit CD8+ T-cell function. TCGA pan-cancer data confirmed that CD8lowPlatelethigh patients have a significantly lower survival probability compared to CD8highPlateletlow. CONCLUSIONS: CD8+ T cells inhibit metastasis. When the balance between CD8+ T cells and platelets is disrupted, platelets produce CXCL4, which induces MDSCs thereby inhibiting the CD8+ T-cell function.

Lipidomic atlas of mammalian cell membranes reveals hierarchical variation induced by culture conditions, subcellular membranes, and cell lineages.

Lipid membranes are ubiquitous biological organizers, required for structural and functional compartmentalization of the cell and sub-cellular organelles. Membranes in living cells are compositionally complex, comprising hundreds of dynamically regulated, distinct lipid species. Cellular physiology requires tight regulation of these lipidomic profiles to achieve proper membrane functionality. While some general features of tissue- and organelle-specific lipid complements have been identified, less is known about detailed lipidomic variations caused by cell-intrinsic or extrinsic factors. Here, we use shotgun lipidomics to report detailed, comprehensive lipidomes of a variety of cultured and primary mammalian membrane preparations to identify trends and sources of variation. Unbiased principle component analysis (PCA) shows clear separation between cultured and primary cells, with primary erythrocytes, synaptic membranes, and other mammalian tissue lipidomes sharply diverging from all cultured cell lines and also from one other. Most broadly, cultured cell membrane preparations were distinguished by their paucity of polyunsaturated lipids. Cultured mammalian cell lines were comparatively similar to one another, although we detected clear, highly reproducible lipidomic signatures of individual cell lines and plasma membrane (PM) isolations thereof. These measurements begin to establish a comprehensive lipidomic atlas of mammalian cells and tissues, identifying some major sources of variation. These observations will allow investigation of the regulation and functional significance of mammalian lipidomes in various contexts.

Methodological Advancements for Investigating Intra-tumoral Heterogeneity in Breast Cancer at the Bench and Bedside.

There is a major need to overcome therapeutic resistance and metastasis that eventually arises in many breast cancer patients. Therapy resistant and metastatic tumors are increasingly recognized to possess intra-tumoral heterogeneity (ITH), a diversity of cells within an individual tumor. First hypothesized in the 1970s, the possibility that this complex ITH may endow tumors with adaptability and evolvability to metastasize and evade therapies is now supported by multiple lines of evidence. Our understanding of ITH has been driven by recent methodological advances including next-generation sequencing, computational modeling, lineage tracing, single-cell technologies, and multiplexed in situ approaches. These have been applied across a range of specimens, including patient tumor biopsies, liquid biopsies, cultured cell lines, and mouse models. In this review, we discuss these approaches and how they have deepened our understanding of the mechanistic origins of ITH amongst tumor cells, including stem cell-like differentiation hierarchies and Darwinian evolution, and the functional role for ITH in breast cancer progression. While ITH presents a challenge for combating tumor evolution, in-depth analyses of ITH in clinical biopsies and laboratory models hold promise to elucidate therapeutic strategies that should ultimately improve outcomes for breast cancer patients.

GSK3ß regulates epithelial-mesenchymal transition and cancer stem cell properties in triple-negative breast cancer.

BACKGROUND: Triple-negative breast cancers (TNBCs), which lack receptors for estrogen, progesterone, and amplification of epidermal growth factor receptor 2, are highly aggressive. Consequently, patients diagnosed with TNBCs have reduced overall and disease-free survival rates compared to patients with other subtypes of breast cancer. TNBCs are characterized by the presence of cancer cells with mesenchymal properties, indicating that the epithelial to mesenchymal transition (EMT) plays a major role in the progression of this disease. The EMT program has also been implicated in chemoresistance, tumor recurrence, and induction of cancer stem cell (CSC) properties. Currently, there are no targeted therapies for TNBC, and hence, it is critical to identify the novel targets to treat TNBC. METHODS: A library of compounds was screened for their ability to inhibit EMT in cells with mesenchymal phenotype as assessed using the previously described Z-cad reporters. Of the several drugs tested, GSK3ß inhibitors were identified as EMT inhibitors. The effects of GSK3ß inhibitors on the properties of TNBC cells with a mesenchymal phenotype were assessed using qRT-PCR, flow cytometry, western blot, mammosphere, and migration and cell viability assays. Publicly available datasets also were analyzed to examine if the expression of GSK3ß correlates with the overall survival of breast cancer patients. RESULTS: We identified a GSK3ß inhibitor, BIO, in a drug screen as one of the most potent inhibitors of EMT. BIO and two other GSK3ß inhibitors, TWS119 and LiCl, also decreased the expression of mesenchymal markers in several different cell lines with a mesenchymal phenotype. Further, inhibition of GSK3ß reduced EMT-related migratory properties of cells with mesenchymal properties. To determine if GSK3ß inhibitors target mesenchymal-like cells by affecting the CSC population, we employed mammosphere assays and profiled the stem cell-related cell surface marker CD44+/24- in cells after exposure to GSK3ß inhibitors. We found that GSK3ß inhibitors indeed decreased the CSC properties of cell types with mesenchymal properties. We treated cells with epithelial and mesenchymal properties with GSK3ß inhibitors and found that GSK3ß inhibitors selectively kill cells with mesenchymal attributes while sparing cells with epithelial properties. We analyzed patient data to identify genes predictive of poor clinical outcome that could serve as novel therapeutic targets for TNBC. The Wnt signaling pathway is critical to EMT, but among the various factors known to be involved in Wnt signaling, only the higher expression of GSK3ß correlated with poorer overall patient survival. CONCLUSIONS: Taken together, our data demonstrate that GSK3ß is a potential target for TNBCs and suggest that GSK3ß inhibitors could serve as selective inhibitors of EMT and CSC properties for the treatment of a subset of aggressive TNBC. GSK3ß inhibitors should be tested for use in combination with standard-of-care drugs in preclinical TNBC models.

Antagonistic actions of two human Pan3 isoforms on global mRNA turnover.

Deadenylation is a fundamental process that regulates eukaryotic gene expression. Mammalian deadenylation exhibits biphasic kinetics, with the Pan2-Pan3 and Ccr4-Caf1 deadenylase complexes mediating the first and second phase, respectively; however, the significance of the biphasic nature of deadenylation in mRNA turnover remains unclear. In this study, we discovered that two distinct isoforms of human Pan3 display opposing properties necessary for coordinating the two phases of deadenylation. The shorter isoform (Pan3S) interacts more strongly with PABP than the longer isoform (Pan3L) does. Pan2 deadenylase activity is enhanced by Pan3S but suppressed by Pan3L. Knocking down individual Pan3 isoforms has opposing effects on the global poly(A) tail length profile, P-body formation, and different mRNA decay pathways. Transcriptome-wide analysis of Pan3 knockdown effects on mRNA turnover shows that depleting either Pan3 isoform causes profound and extensive changes in mRNA stability globally. These results reveal a new fundamental step governing mammalian mRNA metabolism. We propose that the first phase of deadenylation, coordinated through the interplay among the two Pan3 isoforms, Pan2, and PABP, represents a cytoplasmic mRNA maturation step important for proper mRNA turnover.

Bioinformatics service center projects go beyond service.

To probe the complexity of modern diseases, multidisciplinary approaches are increasingly applied. Typically underpinning such studies are collaborations between wet bench experimentalists and dry lab bioinformaticians. Despite the need, bioinformatics collaborators remain difficult to find. Therefore, we undertook a study to understand the nature of this research, so that we may better understand how to meet the needs of future multidisciplinary projects. To accomplish this, we have performed a retrospective study of data from three years of projects performed by the UTHealth Bioinformatics Service Center. Based on this, we found that the bioinformatics in these collaborative projects are extremely diverse and require a high degree of intellectual engagement, while requiring only a small amount of publishable methods development. Very few of the specific skills, the strength of a service core, could be recycled across projects, which were generally exploratory and open-ended and required cycles of biological hypothesis development and (in silico) testing. We find that biomedical research requires bioinformaticians that are highly trained, having the ability to think biologically, but investigating using computational rather than bench experiments. This is in contrast to the activities that are typically the basis for an independent career in biomedical informatics, namely developing new software and algorithms. These findings suggest that to foster team-based multidisciplinary research, institutions must adopt policies that recognize contributions to research by applied bioinformatics scientists.

PCSK9 deficiency reduces atherosclerosis, apolipoprotein B secretion, and endothelial dysfunction.

Proprotein convertase subtilisin/kexin type 9 (PCSK9) interacts directly with cytoplasmic apoB and prevents its degradation via the autophagosome/lysosome pathway. This process affects VLDL and LDL production and influences atherogenesis. Here, we investigated the molecular machinery by which PCSK9 modulates autophagy and affects atherogenesis. We backcrossed Pcsk9-/- mice with atherosclerosis-prone Ldlr-/-Apobec1-/- (LDb) mice to generate Ldlr-/-Apobec1-/-Pcsk9-/- (LTp) mice. Deletion of PCSK9 resulted in decreased hepatic apoB secretion, increased autophagic flux, and decreased plasma levels of IDL and LDL particles. The LDLs from LTp mice (LTp-LDLs) were less atherogenic and contained less cholesteryl ester and phospholipids than LDb-LDLs. Moreover LTp-LDLs induced lower endothelial expression of the genes encoding TLR2, Lox-1, ICAM-1, CCL2, CCL7, IL-6, IL-1ß, Beclin-1, p62, and TRAF6 Collectively, these effects were associated with substantially less atherosclerosis development (>4-fold) in LTp mice. The absence of PCSK9 in LDb mice results in decreased lipid and apoB levels, fewer atherogenic LDLs, and marked reduction of atherosclerosis. The effect on atherogenesis may be mediated in part by the effects of modified LDLs on endothelial cell receptors and proinflammatory and autophagy molecules. These findings suggest that there may be clinical benefits of PCSK9 inhibition due to mechanisms unrelated to increased LDL receptor activity.

Contributions of the RhoA guanine nucleotide exchange factor Net1 to polyoma middle T antigen-mediated mammary gland tumorigenesis and metastasis.

BACKGROUND: The RhoA activating protein Net1 contributes to breast cancer cell proliferation, motility, and invasion in vitro, yet little is known about its roles in mammary gland tumorigenesis and metastasis. METHODS: Net1 knockout (KO) mice were bred to mice with mammary gland specific expression of the polyoma middle T antigen (PyMT) oncogene. Mammary gland tumorigenesis and lung metastasis were monitored. Individual tumors were assessed for proliferation, apoptosis, angiogenesis, RhoA activation, and activation of PyMT-dependent signaling pathways. Primary tumor cells from wild-type and Net1 KO mice were transplanted into the mammary glands of wild-type, nontumor-bearing mice, and tumor growth and metastasis were assessed. Gene expression in wild-type and Net1 KO tumors was analyzed by gene ontology enrichment and for relative activation of gene expression signatures indicative of signaling pathways important for breast cancer initiation and progression. A gene expression signature indicative of Net1 function was identified. Human breast cancer gene expression profiles were screened for the presence of a Net1 gene expression signature. RESULTS: We show that Net1 makes fundamental contributions to mammary gland tumorigenesis and metastasis. Net1 deletion delays tumorigenesis and strongly suppresses metastasis in PyMT-expressing mice. Moreover, we observe that loss of Net1 reduces cancer cell proliferation, inhibits tumor angiogenesis, and promotes tumor cell apoptosis. Net1 is required for maximal RhoA activation within tumors and for primary tumor cell motility. Furthermore, the ability of PyMT to initiate oncogenic signaling to ERK1/2 and PI3K/Akt1 is inhibited by Net1 deletion. Primary tumor cell transplantation indicates that the reduction in tumor angiogenesis and lung metastasis observed upon Net1 deletion are tumor cell autonomous effects. Using a gene expression signature indicative of Net1 activity, we show that Net1 signaling is activated in 10% of human breast cancers, and that this correlates with elevated proliferation and PI3K pathway activity. We also demonstrate that human breast cancer patients with a high Net1 gene expression signature experience shorter distant metastasis-free survival. CONCLUSIONS: These data indicate that Net1 is required for tumor progression in the PyMT mouse model and suggest that Net1 may contribute to breast cancer progression in humans.

Planning bioinformatics workflows using an expert system.

Motivation: Bioinformatic analyses are becoming formidably more complex due to the increasing number of steps required to process the data, as well as the proliferation of methods that can be used in each step. To alleviate this difficulty, pipelines are commonly employed. However, pipelines are typically implemented to automate a specific analysis, and thus are difficult to use for exploratory analyses requiring systematic changes to the software or parameters used. Results: To automate the development of pipelines, we have investigated expert systems. We created the Bioinformatics ExperT SYstem (BETSY) that includes a knowledge base where the capabilities of bioinformatics software is explicitly and formally encoded. BETSY is a backwards-chaining rule-based expert system comprised of a data model that can capture the richness of biological data, and an inference engine that reasons on the knowledge base to produce workflows. Currently, the knowledge base is populated with rules to analyze microarray and next generation sequencing data. We evaluated BETSY and found that it could generate workflows that reproduce and go beyond previously published bioinformatics results. Finally, a meta-investigation of the workflows generated from the knowledge base produced a quantitative measure of the technical burden imposed by each step of bioinformatics analyses, revealing the large number of steps devoted to the pre-processing of data. In sum, an expert system approach can facilitate exploratory bioinformatic analysis by automating the development of workflows, a task that requires significant domain expertise. Availability and Implementation: https://github.com/jefftc/changlab. Contact: jeffrey.t.chang@uth.tmc.edu.

Lipin-1 regulation of phospholipid synthesis maintains endoplasmic reticulum homeostasis and is critical for triple-negative breast cancer cell survival.

Cancer cells reprogram their metabolism to increase the synthesis of macromolecules for rapid proliferation. Compared to fatty acids, much less is known about the synthesis of phospholipids, which is essential for membrane biogenesis in cancer cells. We found that LPIN1, which encodes lipin-1, a phosphatidic acid phosphatase (PAP) controlling the rate-limiting step in the phospholipid synthesis pathway, is highly up-regulated in basal-like triple-negative breast cancer (TNBC). Moreover, high LPIN1 expression correlates with the poor prognosis of these patients. Knockdown of LPIN1 increases apoptosis in basal-like TNBC cell lines, whereas it has minimal or less effect on normal human mammary gland epithelial cells (HMECs) and estrogen receptor-positive breast cancer cell lines. Fatty acid incorporation and lipidomics analyses showed that LPIN1 knockdown blocks phospholipid synthesis and changes membrane lipid compositions that ultimately induce the activation of 1 of the 3 branches of unfolded protein responses, the inositol-requiring enzyme-1α pathway. We also show for the first time, to our knowledge, that lipin-1 knockdown significantly inhibits tumor growth in vivo using an orthotopic xenograft breast mouse model. Our results suggest that lipin-1 is a potential target for cancer therapy.-He, J., Zhang, F., Tay, L. W. R., Boroda, S., Nian, W., Levental, K. R., Levental, I., Harris, T. E., Chang, J. T., Du, G. Lipin-1 regulation of phospholipid synthesis maintains endoplasmic reticulum homeostasis and is critical for triple-negative breast cancer cell survival.

MiR-26 down-regulates TNF-α/NF-κB signalling and IL-6 expression by silencing HMGA1 and MALT1.

MiR-26 has emerged as a key tumour suppressor in various cancers. Accumulating evidence supports that miR-26 regulates inflammation and tumourigenicity largely through down-regulating IL-6 production, but the underlying mechanism remains obscure. Here, combining a transcriptome-wide approach with manipulation of cellular miR-26 levels, we showed that instead of directly targeting IL-6 mRNA for gene silencing, miR-26 diminishes IL-6 transcription activated by TNF-α through silencing NF-κB signalling related factors HMGA1 and MALT1. We demonstrated that miR-26 extensively dampens the induction of many inflammation-related cytokine, chemokine and tissue-remodelling genes that are activated via NF-κB signalling pathway. Knocking down both HMGA1 and MALT1 by RNAi had a silencing effect on NF-κB-responsive genes similar to that caused by miR-26. Moreover, we discovered that poor patient prognosis in human lung adenocarcinoma is associated with low miR-26 and high HMGA1 or MALT1 levels and not with levels of any of them individually. These new findings not only unravel a novel mechanism by which miR-26 dampens IL-6 production transcriptionally but also demonstrate a direct role of miR-26 in down-regulating NF-κB signalling pathway, thereby revealing a more critical and broader role of miR-26 in inflammation and cancer than previously realized.

A genomic strategy to elucidate modules of oncogenic pathway signaling networks.

Recent studies have emphasized the importance of pathway-specific interpretations for understanding the functional relevance of gene alterations in human cancers. Although signaling activities are often conceptualized as linear events, in reality, they reflect the activity of complex functional networks assembled from modules that each respond to input signals. To acquire a deeper understanding of this network structure, we developed an approach to deconstruct pathways into modules represented by gene expression signatures. Our studies confirm that they represent units of underlying biological activity linked to known biochemical pathway structures. Importantly, we show that these signaling modules provide tools to dissect the complexity of oncogenic states that define disease outcomes as well as response to pathway-specific therapeutics. We propose that this model of pathway structure constitutes a framework to study the processes by which information propogates through cellular networks and to elucidate the relationships of fundamental modules to cellular and clinical phenotypes.

p38 MAPK inhibits breast cancer metastasis through regulation of stromal expansion.

p38 MAPK signaling controls cell growth, proliferation and the cell cycle under stress conditions. However, the function of p38 activation in tumor metastasis is still not well understood. We report that p38 activation in breast cancer cells inhibits tumor metastasis but does not substantially modulate primary tumor growth. Stable p38 knockdown in breast cancer cells suppressed NF-κB p65 activation, inhibiting miR-365 expression and resulting in increased IL-6 secretion. The inhibitory effect of p38 signaling on metastasis was mediated by suppression of mesenchymal stem cell (MSC) migration to the primary tumor and sites of metastasis, where MSCs can differentiate into cancer-associated fibroblasts to promote tumor metastasis. The migration of MSCs to these sites relies on CXCR4-SDF1 signaling in the tumor microenvironment. Analysis of human primary and metastatic breast cancer tumors showed that p38 activation was inversely associated with IL-6 and vimentin expression. This study suggests that combination analysis of p38 MAPK and IL-6 signaling in patients with breast cancer may improve prognosis and treatment of metastatic breast cancer.