scholarly journals Integrated single cell and bulk multi-omics reveals heterogeneity and early changes in pathways associated with cetuximab resistance in HNSCC sensitive cell lines

2019 ◽  
Author(s):  
Luciane Tsukamoto Kagohara ◽  
Fernando Zamuner ◽  
Michael Considine ◽  
Jawara Allen ◽  
Srinivasan Yegnasubramanian ◽  
...  

Identifying potential mechanisms of resistance while tumor cells still respond to therapy is critical to delay acquired resistance. We generated the first comprehensive multi-omics, bulk and single cell data in sensitive head and neck squamous cell carcinoma (HNSCC) cells to identify immediate responses to cetuximab. Two pathways potentially associated with resistance were focus of the study: regulation of receptor tyrosine kinases through the transcription factor TFAP2A, and epithelial-to-mesenchymal transition (EMT) process. Single cell RNA-seq demonstrates heterogeneity, with cell specific TFAP2A and VIM expression profiles in response to treatment. RNA-seq and ATAC-seq reveal global changes within five days of cetuximab therapy, suggesting early onset of mechanisms of resistance; and corroborates cell line heterogeneity, with different TFAP2A targets or EMT markers affected by therapy. Lack of TFAP2A reduces HNSCC growth and is enhanced by cetuximab and JQ1. Regarding the EMT process, short term cetuximab therapy has the strongest effect on inhibiting migration. TFAP2A silencing does not affect cell migration, supporting an independent role for both mechanisms in resistance. Overall, we show that immediate adaptive transcriptional and epigenetic changes induced by cetuximab are heterogeneous and cell type dependent; and independent mechanisms of resistance arise while tumor cells are still sensitive to therapy.

2020 ◽  
Vol 123 (1) ◽  
pp. 101-113 ◽  
Author(s):  
Luciane T. Kagohara ◽  
Fernando Zamuner ◽  
Emily F. Davis-Marcisak ◽  
Gaurav Sharma ◽  
Michael Considine ◽  
...  

Abstract Background Identifying potential resistance mechanisms while tumour cells still respond to therapy is critical to delay acquired resistance. Methods We generated the first comprehensive multi-omics, bulk and single-cell data in sensitive head and neck squamous cell carcinoma (HNSCC) cells to identify immediate responses to cetuximab. Two pathways potentially associated with resistance were focus of the study: regulation of receptor tyrosine kinases by TFAP2A transcription factor, and epithelial-to-mesenchymal transition (EMT). Results Single-cell RNA-seq demonstrates heterogeneity, with cell-specific TFAP2A and VIM expression profiles in response to treatment and also with global changes to various signalling pathways. RNA-seq and ATAC-seq reveal global changes within 5 days of therapy, suggesting early onset of mechanisms of resistance; and corroborates cell line heterogeneity, with different TFAP2A targets or EMT markers affected by therapy. Lack of TFAP2A expression is associated with HNSCC decreased growth, with cetuximab and JQ1 increasing the inhibitory effect. Regarding the EMT process, short-term cetuximab therapy has the strongest effect on inhibiting migration. TFAP2A silencing does not affect cell migration, supporting an independent role for both mechanisms in resistance. Conclusion Overall, we show that immediate adaptive transcriptional and epigenetic changes induced by cetuximab are heterogeneous and cell type dependent; and independent mechanisms of resistance arise while tumour cells are still sensitive to therapy.


2021 ◽  
Vol 118 (19) ◽  
pp. e2102050118
Author(s):  
Abhijeet P. Deshmukh ◽  
Suhas V. Vasaikar ◽  
Katarzyna Tomczak ◽  
Shubham Tripathi ◽  
Petra den Hollander ◽  
...  

The epithelial-to-mesenchymal transition (EMT) plays a critical role during normal development and in cancer progression. EMT is induced by various signaling pathways, including TGF-β, BMP, Wnt–β-catenin, NOTCH, Shh, and receptor tyrosine kinases. In this study, we performed single-cell RNA sequencing on MCF10A cells undergoing EMT by TGF-β1 stimulation. Our comprehensive analysis revealed that cells progress through EMT at different paces. Using pseudotime clustering reconstruction of gene-expression profiles during EMT, we found sequential and parallel activation of EMT signaling pathways. We also observed various transitional cellular states during EMT. We identified regulatory signaling nodes that drive EMT with the expression of important microRNAs and transcription factors. Using a random circuit perturbation methodology, we demonstrate that the NOTCH signaling pathway acts as a key driver of TGF-β–induced EMT. Furthermore, we demonstrate that the gene signatures of pseudotime clusters corresponding to the intermediate hybrid EMT state are associated with poor patient outcome. Overall, this study provides insight into context-specific drivers of cancer progression and highlights the complexities of the EMT process.


2020 ◽  
Vol 8 (Suppl 3) ◽  
pp. A799-A799
Author(s):  
Dhiraj Kumar ◽  
Sreeharsha Gurrapu ◽  
Hyunho Han ◽  
Yan Wang ◽  
Seongyeon Bae ◽  
...  

BackgroundLong non-coding RNAs (lncRNAs) are involved in various biological processes and diseases. Malat1 (metastasis-associated lung adenocarcinoma transcript 1), also known as Neat2, is one of the most abundant and highly conserved nuclear lncRNAs. Several studies have shown that the expression of lncRNA Malat1 is associated with metastasis and serving as a predictive marker for various tumor progression. Metastatic relapse often develops years after primary tumor removal as a result of disseminated tumor cells undergoing a period of latency in the target organ.1–4 However, the correlation of tumor intrinsic lncRNA in regulation of tumor dormancy and immune evasion is largely unknown.MethodsUsing an in vivo screening platform for the isolation of genetic entities involved in either dormancy or reactivation of breast cancer tumor cells, we have identified Malat1 as a positive mediator of metastatic reactivation. To functionally uncover the role of Malat1 in metastatic reactivation, we have developed a knock out (KO) model by using paired gRNA CRISPR-Cas9 deletion approach in metastatic breast and other cancer types, including lung, colon and melanoma. As proof of concept we also used inducible knockdown system under in vivo models. To delineate the immune micro-environment, we have used 10X genomics single cell RNA-seq, ChIRP-seq, multi-color flowcytometry, RNA-FISH and immunofluorescence.ResultsOur results reveal that the deletion of Malat1 abrogates the tumorigenic and metastatic potential of these tumors and supports long-term survival without affecting their ploidy, proliferation, and nuclear speckles formation. In contrast, overexpression of Malat1 leads to metastatic reactivation of dormant breast cancer cells. Moreover, the loss of Malat1 in metastatic cells induces dormancy features and inhibits cancer stemness. Our RNA-seq and ChIRP-seq data indicate that Malat1 KO downregulates several immune evasion and stemness associated genes. Strikingly, Malat1 KO cells exhibit metastatic outgrowth when injected in T cells defective mice. Our single-cell RNA-seq cluster analysis and multi-color flow cytometry data show a greater proportion of T cells and reduce Neutrophils infiltration in KO mice which indicate that the immune microenvironment playing an important role in Malat1-dependent immune evasion. Mechanistically, loss of Malat1 is associated with reduced expression of Serpinb6b, which protects the tumor cells from cytotoxic killing by the T cells. Indeed, overexpression of Serpinb6b rescued the metastatic potential of Malat1 KO cells by protecting against cytotoxic T cells.ConclusionsCollectively, our data indicate that targeting this novel cancer-cell-initiated domino effect within the immune system represents a new strategy to inhibit tumor metastatic reactivation.Trial RegistrationN/AEthics ApprovalFor all the animal studies in the present study, the study protocols were approved by the Institutional Animal Care and Use Committee(IACUC) of UT MD Anderson Cancer Center.ConsentN/AReferencesArun G, Diermeier S, Akerman M, et al., Differentiation of mammary tumors and reduction in metastasis upon Malat1 lncRNA loss. Genes Dev 2016 Jan 1;30(1):34–51.Filippo G. Giancotti, mechanisms governing metastatic dormancy and reactivation. Cell 2013 Nov 7;155(4):750–764.Gao H, Chakraborty G, Lee-Lim AP, et al., The BMP inhibitor Coco reactivates breast cancer cells at lung metastatic sites. Cell 2012b;150:764–779.Gao H, Chakraborty G, Lee-Lim AP, et al., Forward genetic screens in mice uncover mediators and suppressors of metastatic reactivation. Proc Natl Acad Sci U S A 2014 Nov 18; 111(46): 16532–16537.


2020 ◽  
Author(s):  
Haoyu Ruan ◽  
Yihang Zhou ◽  
Jie Shen ◽  
Yue Zhai ◽  
Ying Xu ◽  
...  

AbstractMetastatic lung cancer accounts for about half of the brain metastases (BM). Development of leptomeningeal metastases (LM) are becoming increasingly common, and its prognosis is still poor despite the advances in systemic and local approaches. Cytology analysis in the cerebrospinal fluid (CSF) remains the diagnostic gold standard. Although several previous studies performed in CSF have offered great promise for the diagnostics and therapeutics of LM, a comprehensive characterization of circulating tumor cells (CTCs) in CSF is still lacking. To fill this critical gap of lung adenocarcinoma LM (LUAD-LM), we analyzed the transcriptomes of 1,375 cells from 5 LUAD-LM patient and 3 control samples using single-cell RNA sequencing technology. We defined CSF-CTCs based on abundant expression of epithelial markers and genes with lung origin, as well as the enrichment of metabolic pathway and cell adhesion molecules, which are crucial for the survival and metastases of tumor cells. Elevated expression of CEACAM6 and SCGB3A2 was discovered in CSF-CTCs, which could serve as candidate biomarkers of LUAD-LM. We identified substantial heterogeneity in CSF-CTCs among LUAD-LM patients and within patient among individual cells. Cell-cycle gene expression profiles and the proportion of CTCs displaying mesenchymal and cancer stem cell properties also vary among patients. In addition, CSF-CTC transcriptome profiling identified one LM case as cancer of unknown primary site (CUP). Our results will shed light on the mechanism of LUAD-LM and provide a new direction of diagnostic test of LUAD-LM and CUP cases from CSF samples.


Author(s):  
Jinfen Wei ◽  
Zixi Chen ◽  
Meiling Hu ◽  
Ziqing He ◽  
Dawei Jiang ◽  
...  

Hypoxia is a characteristic of tumor microenvironment (TME) and is a major contributor to tumor progression. Yet, subtype identification of tumor-associated non-malignant cells at single-cell resolution and how they influence cancer progression under hypoxia TME remain largely unexplored. Here, we used RNA-seq data of 424,194 single cells from 108 patients to identify the subtypes of cancer cells, stromal cells, and immune cells; to evaluate their hypoxia score; and also to uncover potential interaction signals between these cells in vivo across six cancer types. We identified SPP1+ tumor-associated macrophage (TAM) subpopulation potentially enhanced epithelial–mesenchymal transition (EMT) by interaction with cancer cells through paracrine pattern. We prioritized SPP1 as a TAM-secreted factor to act on cancer cells and found a significant enhanced migration phenotype and invasion ability in A549 lung cancer cells induced by recombinant protein SPP1. Besides, prognostic analysis indicated that a higher expression of SPP1 was found to be related to worse clinical outcome in six cancer types. SPP1 expression was higher in hypoxia-high macrophages based on single-cell data, which was further validated by an in vitro experiment that SPP1 was upregulated in macrophages under hypoxia-cultured compared with normoxic conditions. Additionally, a differential analysis demonstrated that hypoxia potentially influences extracellular matrix remodeling, glycolysis, and interleukin-10 signal activation in various cancer types. Our work illuminates the clearer underlying mechanism in the intricate interaction between different cell subtypes within hypoxia TME and proposes the guidelines for the development of therapeutic targets specifically for patients with high proportion of SPP1+ TAMs in hypoxic lesions.


Author(s):  
Noriyuki Satoh ◽  
Hitoshi Tominaga ◽  
Masato Kiyomoto ◽  
Kanako Hisata ◽  
Jun Inoue ◽  
...  

Among chordate taxa, the cephalochordates diverged earlier than urochordates and vertebrates; thus, they retain unique, primitive developmental features. In particular, the amphioxus notochord has muscle-like properties, a feature not seen in urochordates or vertebrates. Amphioxus contains two Brachyury genes, Bra1 and Bra2. Bra2 is reportedly expressed in the blastopore, notochord, somites, and tail bud, in contrast to a low level of Bra1 expression only in notochord. To distinguish the expression profiles of the two Brachyury genes at the single-cell level, we carried out single-cell RNA-seq (scRNA-seq) analysis using the amphioxus, Branchiostoma japonicum. This scRNA-seq analysis classified B. japonicum embryonic cells into 15 clusters at developmental stages from midgastrula to early swimming larva. Brachyury was expressed in cells of clusters 4, 5, 8, and 9. We first confirmed that cluster 8 comprises cells that form somites since this cluster specifically expresses four myogenic factor genes. Cluster 9 contains a larger number of cells with high levels of Bra2 expression and a smaller number of cells with Bra1 expression. Simultaneous expression in cluster 9 of tool-kit genes, including FoxA, Goosecoid, and hedgehog, showed that this cluster comprises cells that form the notochord. Expression of Bra2, but not Bra1, in cells of clusters 4 and 5 at the gastrula stage together with expression of Wnt1 and Caudal indicates that clusters 4 and 5 comprise cells of the blastopore, which contiguously form the tail bud. In addition, Hox1, Hox3, and Hox4 were highly expressed in Bra2-expressing clusters 4, 5, 8, and 9 in a temporally coordinated manner, suggesting roles of anterior Hox genes in specification of mesodermal organs, including somites, notochord, and tail bud. This scRNA-seq analysis therefore highlights differences between the two Brachyury genes in relation to embryonic regions in which they are expressed and their levels of expression. Bra2 is the ancestral Brachyury in amphioxus, since expression in the blastopore is shared with other deuterostomes. On the other hand, Bra1 is a duplicate copy and likely evolved a supplementary function in notochord and somite formation in the Branchiostoma lineage.


2021 ◽  
Vol 5 (1) ◽  
Author(s):  
Yingmei Li ◽  
Dina Polyak ◽  
Layton Lamsam ◽  
Ian David Connolly ◽  
Eli Johnson ◽  
...  

AbstractNon-small cell lung cancer (NSCLC) metastatic to the brain leptomeninges is rapidly fatal, cannot be biopsied, and cancer cells in the cerebrospinal fluid (CSF) are few; therefore, available tissue samples to develop effective treatments are severely limited. This study aimed to converge single-cell RNA-seq and cell-free RNA (cfRNA) analyses to both diagnose NSCLC leptomeningeal metastases (LM), and to use gene expression profiles to understand progression mechanisms of NSCLC in the brain leptomeninges. NSCLC patients with suspected LM underwent withdrawal of CSF via lumbar puncture. Four cytology-positive CSF samples underwent single-cell capture (n = 197 cells) by microfluidic chip. Using robust principal component analyses, NSCLC LM cell gene expression was compared to immune cells. Massively parallel qPCR (9216 simultaneous reactions) on human CSF cfRNA samples compared the relative gene expression of patients with NSCLC LM (n = 14) to non-tumor controls (n = 7). The NSCLC-associated gene, CEACAM6, underwent in vitro validation in NSCLC cell lines for involvement in pathologic behaviors characteristic of LM. NSCLC LM gene expression revealed by single-cell RNA-seq was also reflected in CSF cfRNA of cytology-positive patients. Tumor-associated cfRNA (e.g., CEACAM6, MUC1) was present in NSCLC LM patients’ CSF, but not in controls (CEACAM6 detection sensitivity 88.24% and specificity 100%). Cell migration in NSCLC cell lines was directly proportional to CEACAM6 expression, suggesting a role in disease progression. NSCLC-associated cfRNA is detectable in the CSF of patients with LM, and corresponds to the gene expression profile of NSCLC LM cells. CEACAM6 contributes significantly to NSCLC migration, a hallmark of LM pathophysiology.


2020 ◽  
Author(s):  
Jinxue Zhang ◽  
Yuan Zhang ◽  
Yongming Liu ◽  
Xin Yi ◽  
Shiyang Cheng ◽  
...  

Abstract Background: Leukocyte-associated immunoglobulin-like receptor-1 (LAIR-1) is a collagen receptor belonging to the immunoglobulin superfamily. Although prior studies have evaluated the biological role of LAIR in solid tumors, the precise mechanisms underlying LAIR-1 functions as a regulator of tumor biological functions remains unclear. Methods: LAIR-1 expression was evaluated using an osteosarcoma (OS) tissue microarray by immunohistochemical analysis. Wound healing and Transwell assays were performed to evaluate tumor cell migration. Quantitative PCR and western blotting were conducted to detect the expression of epithelial-mesenchymal transition (EMT)-related molecules. RNA-sequencing (RNA-seq) was conducted to evaluate the mRNA expression profiles after overexpressing LAIR-1 in OS cells. Glucose uptake and glucose transporter (Glut) 1 expression in OS cells in vitro were evaluated by flow cytometry and western blotting. Results: LAIR-1 expression significantly differed between the T1 and T2 stages of OS tumors, and LAIR-1 overexpression inhibited OS cell migration. LAIR-1 expression was inversely correlated with the expression of EMT-associated transcription factors via the Forkhead box O1/Twist1 signal transduction pathway. Furthermore, RNA-seq and quantitative PCR demonstrated that EMT energy metabolism-related molecules were significantly reduced after LAIR-1 overexpression. Conclusions: Notably, overexpression of LAIR-1 in OS cells decreased Glut1 expression. These findings provide insight into the molecular mechanism underlying OS progression.


2019 ◽  
Author(s):  
Gabriela S. Kinker ◽  
Alissa C. Greenwald ◽  
Rotem Tal ◽  
Zhanna Orlova ◽  
Michael S. Cuoco ◽  
...  

AbstractCultured cell lines are the workhorse of cancer research, but it is unclear to what extent they recapitulate the cellular heterogeneity observed among malignant cells in tumors, given the absence of a native tumor microenvironment. Here, we used multiplexed single cell RNA-seq to profile ~200 cancer cell lines. We uncovered expression programs that are recurrently heterogeneous within many cancer cell lines and are largely independent of observed genetic diversity. These programs of heterogeneity are associated with diverse biological processes, including cell cycle, senescence, stress and interferon responses, epithelial-to-mesenchymal transition, and protein maturation and degradation. Notably, some of these recurrent programs recapitulate those seen in human tumors, suggesting a prominent role of intrinsic plasticity in generating intra-tumoral heterogeneity. Moreover, the data allowed us to prioritize specific cell lines as model systems of cellular plasticity. We used two such models to demonstrate the dynamics, regulation and drug sensitivities associated with a cancer senescence program also observed in human tumors. Our work describes the landscape of cellular heterogeneity in diverse cancer cell lines, and identifies recurrent patterns of expression heterogeneity that are shared between tumors and specific cell lines and can thus be further explored in follow up studies.


Sign in / Sign up

Export Citation Format

Share Document