scholarly journals Single-Cell Gene Network Analysis and Transcriptional Landscape of MYCN-Amplified Neuroblastoma Cell Lines

Biomolecules ◽  
2021 ◽  
Vol 11 (2) ◽  
pp. 177
Author(s):  
Daniele Mercatelli ◽  
Nicola Balboni ◽  
Alessandro Palma ◽  
Emanuela Aleo ◽  
Pietro Paolo Sanna ◽  
...  

Neuroblastoma (NBL) is a pediatric cancer responsible for more than 15% of cancer deaths in children, with 800 new cases each year in the United States alone. Genomic amplification of the MYC oncogene family member MYCN characterizes a subset of high-risk pediatric neuroblastomas. Several cellular models have been implemented to study this disease over the years. Two of these, SK-N-BE-2-C (BE2C) and Kelly, are amongst the most used worldwide as models of MYCN-Amplified human NBL. Here, we provide a transcriptome-wide quantitative measurement of gene expression and transcriptional network activity in BE2C and Kelly cell lines at an unprecedented single-cell resolution. We obtained 1105 Kelly and 962 BE2C unsynchronized cells, with an average number of mapped reads/cell of roughly 38,000. The single-cell data recapitulate gene expression signatures previously generated from bulk RNA-Seq. We highlight low variance for commonly used housekeeping genes between different cells (ACTB, B2M and GAPDH), while showing higher than expected variance for metallothionein transcripts in Kelly cells. The high number of samples, despite the relatively low read coverage of single cells, allowed for robust pathway enrichment analysis and master regulator analysis (MRA), both of which highlight the more mesenchymal nature of BE2C cells as compared to Kelly cells, and the upregulation of TWIST1 and DNAJC1 transcriptional networks. We further defined master regulators at the single cell level and showed that MYCN is not constantly active or expressed within Kelly and BE2C cells, independently of cell cycle phase. The dataset, alongside a detailed and commented programming protocol to analyze it, is fully shared and reusable.

2019 ◽  
Author(s):  
JM Robinson

AbstractThis brief report details results from a comparative analysis of Nanostring expression data between cell lines HEPG2, Caco-2, HT-29, and colon fibroblasts. Raw and normalized data are available publicly in the NCBI GEO/Bioproject databases. Results identify cell-line specific variations in gene expression relevant to intestinal epithelial function.


2013 ◽  
Vol 15 (4) ◽  
pp. 363-372 ◽  
Author(s):  
Victoria Moignard ◽  
Iain C. Macaulay ◽  
Gemma Swiers ◽  
Florian Buettner ◽  
Judith Schütte ◽  
...  

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

2019 ◽  
Author(s):  
Chiaowen Joyce Hsiao ◽  
PoYuan Tung ◽  
John D. Blischak ◽  
Jonathan E. Burnett ◽  
Kenneth A. Barr ◽  
...  

AbstractCellular heterogeneity in gene expression is driven by cellular processes such as cell cycle and cell-type identity, and cellular environment such as spatial location. The cell cycle, in particular, is thought to be a key driver of cell-to-cell heterogeneity in gene expression, even in otherwise homogeneous cell populations. Recent advances in single-cell RNA-sequencing (scRNA-seq) facilitate detailed characterization of gene expression heterogeneity, and can thus shed new light on the processes driving heterogeneity. Here, we combined fluorescence imaging with scRNA-seq to measure cell cycle phase and gene expression levels in human induced pluripotent stem cells (iPSCs). Using these data, we developed a novel approach to characterize cell cycle progression. While standard methods assign cells to discrete cell cycle stages, our method goes beyond this, and quantifies cell cycle progression on a continuum. We found that, on average, scRNA-seq data from only five genes predicted a cell’s position on the cell cycle continuum to within 14% of the entire cycle, and that using more genes did not improve this accuracy. Our data and predictor of cell cycle phase can directly help future studies to account for cell-cycle-related heterogeneity in iPSCs. Our results and methods also provide a foundation for future work to characterize the effects of the cell cycle on expression heterogeneity in other cell types.


Blood ◽  
2015 ◽  
Vol 126 (23) ◽  
pp. 4249-4249
Author(s):  
Amit Kumar Mitra ◽  
Ujjal Mukherjee ◽  
Taylor Harding ◽  
Holly Stessman ◽  
Ying Li ◽  
...  

Abstract Multiple myeloma (MM) is characterized by significant genetic diversity at subclonal levels that likely plays a defining role in the heterogeneity of tumor progression, clinical aggressiveness and drug sensitivity. Such heterogeneity is a driving factor in the evolution of MM, from founder clones through outgrowth of subclonal fractions. DNA Sequencing studies on MM samples have indeed demonstrated such heterogeneity in subclonal architecture at diagnosis based on recurrent mutations in pathologically relevant genes that may ultimately to lead to relapse. However, no study so far has reported a predictive gene expression signature that can identify, distinguish and quantify drug sensitive and drug-resistant subpopulations within a bulk population of myeloma cells. In recent years, our laboratory has successfully developed a gene expression profile (GEP)-based signature that could not only distinguish drug response of MM cell lines, but also was effective in stratifying patient outcomes when applied to GEP profiles from MM clinical trials using proteasome inhibitors (PI) as chemotherapeutic agents. Further, we noted myeloma cell lines that responded to the drug often contained residual sub-population of cells that did not respond, and likely were selectively propagated during drug treatment in vitro, and in patients. In this study, we performed targeted qRT-PCR analysis of single cells using a gene panel that included PI sensitivity genes and gene signatures that could discriminate between low and high-risk myeloma followed by intensive bioinformatics and statistical analysis for the classification and prediction of PI response in individual cells within bulk multiple myeloma tumors. Fluidigm's C1 Single-Cell Auto Prep System was used to perform automated single-cell capture, processing and cDNA synthesis on 576 pre-treatment cells from 12 cell lines representing a wide range of PI-sensitivity and 370 cells from 7 patient samples undergoing PI treatment followed by targeted gene expression profiling of single cells using automated, high-throughput on-chip qRT-PCR analysis using 96.96 Dynamic Array IFCs on the BioMark HD System. Probability of resistance for each individual cell was predicted using a pipeline that employed the machine learning methods Random Forest, Support Vector Machine (radial and sigmoidal), LASSO and kNN (k Nearest Neighbor) for making single-cell GEP data-driven predictions/ decisions. The weighted probabilities from each of the algorithms were used to quantify resistance of each individual cell and plotted using Ensemble forecasting algorithm. Using our drug response GEP signature at the single cell level, we could successfully identify distinct subpopulations of tumor cells that were predicted to be sensitive or resistant to PIs. Subsequently, we developed a R Statistical analysis package (http://cran.r-project.org), SCATTome (Single Cell Analysis of Targeted Transcriptome), that can restructure data obtained from Fluidigm qPCR analysis run, filter missing data, perform scaling of filtered data, build classification models and successfully predict drug response of individual cells and classify each cell's probability of response based on the targeted transcriptome. We will present the program output as graphical displays of single cell response probabilities. This package provides a novel classification method that has the potential to predict subclonal response to a variety of therapeutic agents. Disclosures Kumar: Skyline: Consultancy, Honoraria; BMS: Consultancy; Onyx: Consultancy, Research Funding; Sanofi: Consultancy, Research Funding; Janssen: Consultancy, Research Funding; Novartis: Research Funding; Takeda: Consultancy, Research Funding; Celgene: Consultancy, Research Funding.


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.


2021 ◽  
Vol 23 (Supplement_2) ◽  
pp. ii10-ii10
Author(s):  
K Joseph ◽  
L Vollmer ◽  
V M Ravi ◽  
J Beck ◽  
U G Hofmann ◽  
...  

Abstract BACKGROUND Owing to recent advances in understanding of the active functional states exhibited within glioblastoma (GBM), intra-tumoral cellular signaling has moved into focus of neuro-oncological research. In our study, we aim to explore the diversity of transcellular signaling and investigate correlations to transcriptional dynamics and cellular behavior. MATERIAL AND METHODS Electrophysiological mapping of primary GBM cultures was performed by planar microelectrodes, in conjunction with calcium imaging in a human neocortical section based GBM model. Exposure to conditions that are physiologically present within the tumor was carried out to identify specific signaling cells of interest and signaling diversity presented as response to specific environmental conditions. Transcriptional dynamics and plasticity were examined by means of scRNA-sequencing with CRISPR based perturbation, spatial transcriptomics and deep long-read RNA-sequencing. RESULTS Electrophysiological profiles of primary GBM cell lines revealed highly variable network activity. Despite these different characteristics, all profiled primary cell-lines exhibited characteristics of scale-free networks, confirmed in a human neocortical GBM model. When the GBM was allowed to grow in “in-vivo” like environment, basal activity was significantly increased, owing to interactions with elements within the neural environment. Cellular signaling was directly correlated to changes in the environment, like hypoxia or glutamatergic activation, and total inhibition of electrical signaling was achieved only with a combination of both gap junction and synaptic inhibitors. Using single-cell sequencing and proteomics, we identified several genes related to synaptogenesis that plays a crucial role in network formation and consequently transcellular signaling. CRISPR based perturbation of these genes resulted in alterations in cellular morphology and decreased cellular connectivity, with electrical signaling being significantly attenuated. Single-cell sequencing of perturbed tumor cells in the GBM model revealed a loss of developmental lineages and significant reduction of cellular stress response state. CONCLUSION Our findings highlight the role of electrical signaling in glioblastoma. Cellular stressors induce intercellular signaling, leading to transcriptional adaptation suggesting that there exists a highly complex and powerful mechanism for dynamic transcriptional state adaptation.


Blood ◽  
2014 ◽  
Vol 124 (21) ◽  
pp. 3385-3385 ◽  
Author(s):  
Amit Kumar Mitra ◽  
Holly Stessman ◽  
Michael A. Linden ◽  
Brian Van Ness

Abstract Multiple myeloma (MM) is a plasma cell neoplasm with significant complexity and heterogeneity. Proteasome inhibitors (PI) including bortezomib (Velcade/Bz), carfilzomib (Kyprolis/Cz) and Ixazomib are effective chemotherapeutic agents in the treatment of MM, used alone or in combination with other anti-cancer agents. However, in spite of the recent improvements in treatment strategies, MM still remains a difficult disease to cure with median survival rate of around 7 years. In a recently published study, we have shown that the heterogeneity in response to proteasome inhibitor (PI)-based treatment in MM is governed by underlying molecular characteristics of the subclones within tumor population (Stessman et al. 2013). We confirmed the presence of residual resistant sub-population comprising up to 15% of the bulk Bz-sensitive cell population in drug-naïve MM tumors. We hypothesize that this pre-existing resistant sub-population may give rise to emerging resistance in course of treatment with PIs. In the current study, we used single cell transcriptomics analysis to identify tumor subclones within Human Myeloma Cell Lines (HMCLs) based on a 48-gene model of predictive genetic signature for baseline PI response. Automated single-cell capture and cDNA synthesis from cellular RNA were performed using Fluidigm’s C1TM Single-Cell Auto Prep System. The cDNA was then harvested and transferred to BioMark HD System for single-cell targeted high-throughput qPCR-based gene expression analysis of a 48 gene-panel using Fluidigm DELTAgene assays. Our 48-gene model combines our previously published 23 gene expression profiling (GEP) signature that could discriminate between sensitive and resistant responsiveness to Bz, and the Shaughnessy et al prognostic 17-gene GEP model along with control genes, including cell cycle genes, anti-apoptotic genes, proteasome subunit genes, house-keeping genes and internal negative controls. Based on the differential expression of these 48 genes used in the modeling, distinct subclonal populations were then identified using a combination of Fluidigm’s analysis software and the R Statistical analysis package. Further, a principal component analysis (PCA) score plot was generated as a two-dimensional grid to visualize the separate populations associated with resistant profiles. Finally, hierarchical clustering (HC) analysis was used to generate heat maps that group expression patterns associated with response. Our results demonstrated the presence of pre-existing subclones of cells within untreated myeloma cells with a characteristic genetic signature profile distinct from the pre-treatment overall (bulk) profile of myeloma cells. As an additional validation of subclonal architecture, we demonstrated the presence of subclones within HMCLs using multi-color flow cytometry. The results presented will help identify the presence and extent of intra-tumor heterogeneity in MM by single cell transcriptomics and may define residual pre-existing subclones resistant to PI therapies. Disclosures No relevant conflicts of interest to declare.


2021 ◽  
Author(s):  
◽  
Daniel R. Kick

Neural networks produce critical rhythmic behaviors throughout an animal's lifespan, despite growth, differing environments, and changes in physiological state. This requires networks which balance stability in their properties with the plasticity necessary to respond to altered demands or perturbations. Studying the mechanisms which confer these properties requires a well characterized system with a known network topology and identifiable neurons that are amenable to both electrophysiological and molecular characterization and manipulation. Here, we use two networks from Cancer borealis to explore activity dependent regulation of cell connectivity, changes in cell properties with prolonged perturbation, and reliability of gene expression as a means for cell identification. For the first two topics we use the cardiac ganglion alone. The cardiac ganglion consists of a kernel of four interneurons that drive five motor neurons (termed large cells, LCs) which innervate the heart musculature. LCs burst synchronously due to simultaneous stimulation and electrical coupling through gap junctions. Depolarizing pharmacological perturbations have been shown to result in hyperexcitability (Ransdell et al., 2012a) and disrupt synchrony between LCs (Lane et al., 2016) eliciting rapid plasticity in ionic currents and electrical coupling which restores synchrony and excitability (Ransdell et al., 2012a; Lane et al., 2016). The salient electrophysiological signal which elicits coupling plasticity has not been identified. Using voltage clamp we directly control LC depolarizations to vary amplitude and timing of activity between LCs. We find that timing between cells, rather than depolarization elicits plasticity with the direction, i.e., potentiation or depression, being determined by the degree of desynchronization. With dynamic clamp we artificially couple networks from two animals and show that strong coupling with sufficient desynchronization can compromise a cell's output. These results suggest that coupling strength is tuned promoting synchrony or baseline cellular activity in a degree dependent manner. While rapid compensatory plasticity to hyperexcitability has been shown, it is unknown whether the changes are solely post-transcriptional and whether the short-term changes persist over longer time scales. We perturb networks for one or twenty-four hours and compare LCs' excitability, membrane properties, and abundances of ion channel and gap junction transcripts. We find evidence of rapid transcriptional changes at one hour, which may be maintained or regress at twenty-four hours. Additionally, we find that membrane properties and excitability are not maintained from one to twenty-four hours, suggesting a failure to maintain homeostasis or that additional compensatory changes are occurring at the network level. To address our third topic, we use LCs in addition to neurons collected form the stomatogastric ganglion which coordinates mastication and filtering in the digestive track. Both systems allow for unambiguous identification of cells based on anatomy or neuronal projections. We use this to evaluate the efficacy of cluster estimation procedures, clustering methods, and classification algorithms to determine the number of cell types present, group like cells together, and identify cells based on gene expression alone. We use single cell RNA-seq and single cell qRT-PCR to measure all contigs or a select set of ion channel, receptor, and gap junction mRNAs. We find these methods do not reproduce the known number of cell types present. Furthermore, although clustering and classification both outperform chance, we are unable to recapitulate cell type with complete accuracy from these data. These results indicate that, while promising, determining cell type by molecular profiling should not be relied on as the sole metric of cell type determination.


Blood ◽  
2008 ◽  
Vol 111 (2) ◽  
pp. 806-815 ◽  
Author(s):  
Unn-Merete Fagerli ◽  
Randi U. Holt ◽  
Toril Holien ◽  
Thea K. Vaatsveen ◽  
Fenghuang Zhan ◽  
...  

Multiple myeloma (MM) is characterized by accumulation and dissemination of malignant plasma cells (PCs) in the bone marrow (BM). Gene expression profiling of 2 MM cell lines (OH-2 and IH-1) indicated that expression of PRL-3, a metastasis-associated tyrosine phosphatase, was induced by several mitogenic cytokines. Cytokine-driven PRL-3 expression could be shown in several myeloma cell lines at both the mRNA and protein levels. There was significantly higher expression of the PRL-3 gene in PCs from patients with monoclonal gammopathy of undetermined significance (MGUS), smoldering myeloma (SMM), and myeloma than in PCs from healthy persons. Among 7 MM subgroups identified by unsupervised hierarchical cluster analysis, PRL-3 gene expression was significantly higher in the 3 groups denoted as “proliferation,” “low bone disease,” and “MMSET/FGFR3.” PRL-3 protein was detected in 18 of 20 BM biopsies from patients with MM. Silencing of the PRL-3 gene by siRNA reduced cell migration in the MM cell line INA-6, but had no detectable effect on proliferation and cell-cycle phase distribution of the cells. In conclusion, PRL-3 is a gene product specifically expressed in malignant plasma cells and may have a role in migration of these cells.


Sign in / Sign up

Export Citation Format

Share Document