scholarly journals Quartz-Seq2: a high-throughput single-cell RNA-sequencing method that effectively uses limited sequence reads

2017 ◽  
Author(s):  
Yohei Sasagawa ◽  
Hiroki Danno ◽  
Hitomi Takada ◽  
Masashi Ebisawa ◽  
Kaori Tanaka ◽  
...  

AbstractHigh-throughput single-cell RNA-seq methods assign limited unique molecular identifier (UMI) counts as gene expression values to single cells from shallow sequence reads and detect limited gene counts. We thus developed a high-throughput single-cell RNA-seq method, Quartz-Seq2, to overcome these issues. Our improvements in the reaction steps make it possible to effectively convert initial reads to UMI counts (at a rate of 30%–50%) and detect more genes. To demonstrate the power of Quartz-Seq2, we analyzed approximately 10,000 transcriptomes in total from in vitro embryonic stem cells and an in vivo stromal vascular fraction with a limited number of reads.

Blood ◽  
2018 ◽  
Vol 132 (Supplement 1) ◽  
pp. 3887-3887
Author(s):  
Moosa Qureshi ◽  
Fernando Calero-Nieto ◽  
Iwo Kucinski ◽  
Sarah Kinston ◽  
George Giotopoulos ◽  
...  

Abstract The C/EBPα transcription factor plays a pivotal role in myeloid differentiation and E2F-mediated cell cycle regulation. Although CEBPA mutations are common in acute myeloid leukaemia (AML), little is known regarding pre-leukemic alterations caused by mutated CEBPA. Here, we investigated early events involved in pre-leukemic transformation driven by CEBPA N321D in the LMPP-like cell line Hoxb8-FL (Redecke et al., Nat Methods 2013), which can be maintained in vitro as a self-renewing LMPP population using Flt3L and estradiol, as well as differentiated both in vitro and in vivo into myeloid and lymphoid cell types. Hoxb8-FL cells were retrovirally transduced with Empty Vector (EV), wild-type CEBPA (CEBPA WT) or its N321D mutant form (CEBPA N321D). CEBPA WT-transduced cells showed increased expression of cd11b and SIRPα and downregulation of c-kit, suggesting that wild-type CEBPA was sufficient to promote differentiation even under LMPP growth conditions. Interestingly, we did not observe the same phenotype in CEBPA N321D-transduced cells. Upon withdrawal of estradiol, both EV and CEBPA WT-transduced cells differentiated rapidly into a conventional dendritic cell (cDC) phenotype by day 7 and died within 12 days. By contrast, CEBPA N321D-transduced cells continued to grow for in excess of 56 days, with an initial cDC phenotype but by day 30 demonstrating a plasmacytoid dendritic cell precursor phenotype. CEBPA N321D-transduced cells were morphologically distinct from EV-transduced cells. To test leukemogenic potential in vivo, we performed transplantation experiments in lethally irradiated mice. Serial monitoring of peripheral blood demonstrated that Hoxb8-FL derived cells had disappeared by 4 weeks, and did not reappear. However, at 6 months CEBPA N321D-transduced cells could still be detected in bone marrow in contrast to EV-transduced cells but without any leukemic phenotype. To identify early events involved in pre-leukemic transformation, the differentiation profiles of EV, CEBPA WT and CEBPA N321D-transduced cells were examined with single cell RNA-seq (scRNA-seq). 576 single cells were taken from 3 biological replicates at days 0 and 5 post-differentiation, and analysed using the Automated Single-Cell Analysis Pipeline (Gardeux et al., Bioinformatics 2017). Visualisation by t-SNE (Fig 1) demonstrated: (i) CEBPA WT-transduced cells formed a distinct cluster at day 0 before withdrawal of estradiol; (ii) CEBPA N321D-transduced cells separated from EV and CEBPA WT-transduced cells after 5 days of differentiation, (iii) two subpopulations could be identified within the CEBPA N321D-transduced cells at day 5, with a cluster of five CEBPA N321D-transduced single cells distributed amongst or very close to the day 0 non-differentiated cells. Differential expression analysis identified 224 genes upregulated and 633 genes downregulated specifically in the CEBPA N321D-transduced cells when compared to EV cells after 5 days of differentiation. This gene expression signature revealed that CEBPA N321D-transduced cells switched on a HSC/MEP/CMP transcriptional program and switched off a myeloid dendritic cell program. Finally, in order to further dissect the effect of the N321D mutation, the binding profile of endogenous and CEBPA N321D was compared by ChIP-seq before and after 5 days of differentiation. Integration with scRNA-seq data identified 160 genes specifically downregulated in CEBPA N321D-transduced cells which were associated with the binding of the mutant protein. This list of genes included genes previously implicated in dendritic cell differentiation (such as NOTCH2, JAK2), as well as a number of genes not previously implicated in the evolution of AML, representing potentially novel therapeutic targets. Disclosures No relevant conflicts of interest to declare.


2021 ◽  
Author(s):  
Xiaozhong Shen ◽  
Gangcai Xie

AbstractN(6)-methyladenosine (m(6)a) is the most common internal modification of messenger RNA (mRNA) in higher eukaryotes. According to previous literature reports, alkbh5, as another demethylase in mammals, can reverse the expression of m(6)a gene in vivo and in vitro. In order to reveal the effect of Alkbh5 deletion on the level of single cells in the testis during spermatogenesis in mice, the data were compared using single-cell sequencing. In this article, we discussed the transcription profile and cell type identification of mouse testis, the expression of mitochondrial and ribosomal genes in mice, the analysis of differential gene expression, and the effects of Alkbh5 deletion, and try to explain the role and influence of Alkbh5 on reproduction at the level of single-cell sequencing.


BMC Genomics ◽  
2021 ◽  
Vol 22 (1) ◽  
Author(s):  
Tracy M. Yamawaki ◽  
Daniel R. Lu ◽  
Daniel C. Ellwanger ◽  
Dev Bhatt ◽  
Paolo Manzanillo ◽  
...  

Abstract Background Elucidation of immune populations with single-cell RNA-seq has greatly benefited the field of immunology by deepening the characterization of immune heterogeneity and leading to the discovery of new subtypes. However, single-cell methods inherently suffer from limitations in the recovery of complete transcriptomes due to the prevalence of cellular and transcriptional dropout events. This issue is often compounded by limited sample availability and limited prior knowledge of heterogeneity, which can confound data interpretation. Results Here, we systematically benchmarked seven high-throughput single-cell RNA-seq methods. We prepared 21 libraries under identical conditions of a defined mixture of two human and two murine lymphocyte cell lines, simulating heterogeneity across immune-cell types and cell sizes. We evaluated methods by their cell recovery rate, library efficiency, sensitivity, and ability to recover expression signatures for each cell type. We observed higher mRNA detection sensitivity with the 10x Genomics 5′ v1 and 3′ v3 methods. We demonstrate that these methods have fewer dropout events, which facilitates the identification of differentially-expressed genes and improves the concordance of single-cell profiles to immune bulk RNA-seq signatures. Conclusion Overall, our characterization of immune cell mixtures provides useful metrics, which can guide selection of a high-throughput single-cell RNA-seq method for profiling more complex immune-cell heterogeneity usually found in vivo.


2019 ◽  
Author(s):  
Ugur M. Ayturk ◽  
Joseph P. Scollan ◽  
Alexander Vesprey ◽  
Christina M. Jacobsen ◽  
Paola Divieti Pajevic ◽  
...  

ABSTRACTSingle cell RNA-seq (scRNA-seq) is emerging as a powerful technology to examine transcriptomes of individual cells. We determined whether scRNA-seq could be used to detect the effect of environmental and pharmacologic perturbations on osteoblasts. We began with a commonly used in vitro system in which freshly isolated neonatal mouse calvarial cells are expanded and induced to produce a mineralized matrix. We used scRNA-seq to compare the relative cell type abundances and the transcriptomes of freshly isolated cells to those that had been cultured for 12 days in vitro. We observed that the percentage of macrophage-like cells increased from 6% in freshly isolated calvarial cells to 34% in cultured cells. We also found that Bglap transcripts were abundant in freshly isolated osteoblasts but nearly undetectable in the cultured calvarial cells. Thus, scRNA-seq revealed significant differences between heterogeneity of cells in vivo and in vitro. We next performed scRNA-seq on freshly recovered long bone endocortical cells from mice that received either vehicle or Sclerostin-neutralizing antibody for 1 week. Bone anabolism-associated transcripts were also not significantly increased in immature and mature osteoblasts recovered from Sclerostin-neutralizing antibody treated mice; this is likely a consequence of being underpowered to detect modest changes in gene expression, since only 7% of the sequenced endocortical cells were osteoblasts, and a limited portion of their transcriptomes were sampled. We conclude that scRNA-seq can detect changes in cell abundance, identity, and gene expression in skeletally derived cells. In order to detect modest changes in osteoblast gene expression at the single cell level in the appendicular skeleton, larger numbers of osteoblasts from endocortical bone are required.


Blood ◽  
2021 ◽  
Vol 138 (Supplement 1) ◽  
pp. 1541-1541
Author(s):  
Mary T Scott ◽  
Wei Liu ◽  
Rebecca Mitchell ◽  
Cassie Clarke ◽  
Hassan Almasoudi ◽  
...  

Abstract Although it has been recognized for many years that cancer stem cells and embryonic stem cells (ESC) share molecular features, identifying ways to exploit this therapeutically has proved challenging. To date, these shared features have not been examined in the leukemic stem cells (LSC) found in patients with chronic myeloid leukemia (CML). By integrating known ES regulatory circuitry with transcriptomics datasets, including deep single-cell RNA-seq profiling of 15,670 LSC from five patients with CML, we identified a core ESC regulome in the LSC containing 1243 genes. The significant majority of this regulome (1102 genes) was up-regulated in cycling LSC, whilst quiescent LSC showed up-regulation of a characteristic set of 101 genes, unique to cells with high ESC identity and with regulatory circuitry enriched for c-Myc and Nanog modules. Membership of the ESC regulome included the TP53 gene which was transcriptionally repressed and detected at a lower frequency in quiescent LSC compared to cycling ones (11.8% vs 43.6%). We also demonstrated that tyrosine kinase inhibitors (TKI) repress the ESC regulome and TP53 expression in LSC, suggesting that the regulome safeguards against high levels of TP53 expression, thus promoting survival of quiescent LSC in the presence of TKI. We hypothesized that overcoming the influence of the regulome on TP53 expression would provide an opportunity to eradicate quiescent LSC. To this end, we used an MDM2 inhibitor (MDM2i), RG7388 (idasanutlin) or RG7112, to stabilize the p53 protein, examining its potential in combination with nilotinib (NIL) to eradicate CML LSC in vitro and in vivo, with RG7388 being the most optimized and furthest in development. The combination of NIL plus MDM2i in vitro was more effective at targeted LSC from primary patient samples than NIL treatment alone, as evidenced by reduced CFC and LTC-IC outputs (p<0.05, 0.01 respectively). Intriguingly, the combination of NIL plus MDM2i did not result in significant reductions in the number of LSC compared to NIL only, when we quantified them at the end of drug treatments in pre-clinical mouse models. Instead, we observed a functional decline of the LSC as evidenced by diminished engraftment potential in 2 o recipient mice (p<0.05; SCLtTA x BCR-ABL1 transgenic model) or diminished colony-plating potential (p<0.05). This was followed by near complete depletion of the LSC population (p<0.05) 28 days after cessation of combination drug treatment (patient-derived xenografts/PDX in immunocompromised mice). In order to understand the molecular events underpinning these drug effects on LSC, we performed RNA-seq analysis of drug-treated CD34 + cells in vitro (bulk cells), or of human CD34 + cells obtained from PDX (single cell RNA-seq). CD34 + cells treated with NIL plus MDM2i in vitro showed evidence of increased p53 stabilization and activation of p53 target genes, and this was accompanied by repression of the ESC regulome beyond that normally observed with NIL only. Similarly, in PDX we observed increased repression of the ESC regulome in human CD34 + cells exposed to the combination of NIL plus MDM2i that included repression of HIF1alpha and a signature of genes required for cellular adaptations to hypoxia, and growth factor-mediated resistance to TKI therapy. Further, single cell analysis of differentiated human CD45 + cells from our PDX model, provided compelling evidence that acquisition of this repressive signature in the LSC, through combined NIL plus MDM2i treatment, re-wires them towards a basophilic fate, consistent with functional exhaustion of the LSC compartment. In conclusion, we have identified an ESC regulome in CML LSC and demonstrate that a combination of a TKI plus an MDM2i leads to p53 upregulation which antagonizes this regulome, providing a highly effective strategy to target near complete loss of functional LSC in pre-clinical models. Our study has revealed a new therapeutic paradigm to examine in other cancer stem cell populations that utilize ESC regulatory programs. Disclosures Higgins: Roche/Genentech: Current Employment, Current equity holder in publicly-traded company. Copland: Astellas: Honoraria, Speakers Bureau; Novartis: Honoraria, Speakers Bureau; Pfizer: Honoraria, Speakers Bureau; Incyte: Honoraria, Research Funding, Speakers Bureau; Cyclacel Ltd: Research Funding; Jazz: Honoraria, Speakers Bureau.


2021 ◽  
Vol 3 (2) ◽  
Author(s):  
Xin Wang ◽  
Jane Frederick ◽  
Hongbin Wang ◽  
Sheng Hui ◽  
Vadim Backman ◽  
...  

Abstract The transcriptional plasticity of cancer cells promotes intercellular heterogeneity in response to anticancer drugs and facilitates the generation of subpopulation surviving cells. Characterizing single-cell transcriptional heterogeneity after drug treatments can provide mechanistic insights into drug efficacy. Here, we used single-cell RNA-seq to examine transcriptomic profiles of cancer cells treated with paclitaxel, celecoxib and the combination of the two drugs. By normalizing the expression of endogenous genes to spike-in molecules, we found that cellular mRNA abundance shows dynamic regulation after drug treatment. Using a random forest model, we identified gene signatures classifying single cells into three states: transcriptional repression, amplification and control-like. Treatment with paclitaxel or celecoxib alone generally repressed gene transcription across single cells. Interestingly, the drug combination resulted in transcriptional amplification and hyperactivation of mitochondrial oxidative phosphorylation pathway linking to enhanced cell killing efficiency. Finally, we identified a regulatory module enriched with metabolism and inflammation-related genes activated in a subpopulation of paclitaxel-treated cells, the expression of which predicted paclitaxel efficacy across cancer cell lines and in vivo patient samples. Our study highlights the dynamic global transcriptional activity driving single-cell heterogeneity during drug response and emphasizes the importance of adding spike-in molecules to study gene expression regulation using single-cell RNA-seq.


2020 ◽  
Author(s):  
Alina Isakova ◽  
Norma Neff ◽  
Stephen R. Quake

ABSTRACTThe ability to interrogate total RNA content of single cells would enable better mapping of the transcriptional logic behind emerging cell types and states. However, current RNA-seq methods are unable to simultaneously monitor both short and long, poly(A)+ and poly(A)-transcripts at the single-cell level, and thus deliver only a partial snapshot of the cellular RNAome. Here, we describe Smart-seq-total, a method capable of assaying a broad spectrum of coding and non-coding RNA from a single cell. Built upon the template-switch mechanism, Smart-seq-total bears the key feature of its predecessor, Smart-seq2, namely, the ability to capture full-length transcripts with high yield and quality. It also outperforms current poly(A)–independent total RNA-seq protocols by capturing transcripts of a broad size range, thus, allowing us to simultaneously analyze protein-coding, long non-coding, microRNA and other non-coding RNA transcripts from single cells. We used Smart-seq-total to analyze the total RNAome of human primary fibroblasts, HEK293T and MCF7 cells as well as that of induced murine embryonic stem cells differentiated into embryoid bodies. We show that simultaneous measurement of non-coding RNA and mRNA from the same cell enables elucidation of new roles of non-coding RNA throughout essential processes such as cell cycle or lineage commitment. Moreover, we show that cell types can be distinguished based on the abundance of non-coding transcripts alone.


2021 ◽  
Vol 118 (51) ◽  
pp. e2113568118
Author(s):  
Alina Isakova ◽  
Norma Neff ◽  
Stephen R. Quake

The ability to interrogate total RNA content of single cells would enable better mapping of the transcriptional logic behind emerging cell types and states. However, current single-cell RNA-sequencing (RNA-seq) methods are unable to simultaneously monitor all forms of RNA transcripts at the single-cell level, and thus deliver only a partial snapshot of the cellular RNAome. Here we describe Smart-seq-total, a method capable of assaying a broad spectrum of coding and noncoding RNA from a single cell. Smart-seq-total does not require splitting the RNA content of a cell and allows the incorporation of unique molecular identifiers into short and long RNA molecules for absolute quantification. It outperforms current poly(A)-independent total RNA-seq protocols by capturing transcripts of a broad size range, thus enabling simultaneous analysis of protein-coding, long-noncoding, microRNA, and other noncoding RNA transcripts from single cells. We used Smart-seq-total to analyze the total RNAome of human primary fibroblasts, HEK293T, and MCF7 cells, as well as that of induced murine embryonic stem cells differentiated into embryoid bodies. By analyzing the coexpression patterns of both noncoding RNA and mRNA from the same cell, we were able to discover new roles of noncoding RNA throughout essential processes, such as cell cycle and lineage commitment during embryonic development. Moreover, we show that independent classes of short-noncoding RNA can be used to determine cell-type identity.


Author(s):  
Eszter Posfai ◽  
John Paul Schell ◽  
Adrian Janiszewski ◽  
Isidora Rovic ◽  
Alexander Murray ◽  
...  

AbstractTotipotency is the ability of a single cell to give rise to all the differentiated cells that build the conceptus, yet how to capture this property in vitro remains incompletely understood. Defining totipotency relies upon a variety of assays of variable stringency. Here we describe criteria to define totipotency. We illustrate how distinct criteria of increasing stringency can be used to judge totipotency by evaluating candidate totipotent cell types in the mouse, including early blastomeres and expanded or extended pluripotent stem cells. Our data challenge the notion that expanded or extended pluripotent states harbor increased totipotent potential relative to conventional embryonic stem cells under in vivo conditions.


2020 ◽  
Author(s):  
Chi-Ming Kevin Li ◽  
Tracy M Yamawaki ◽  
Daniel R Lu ◽  
Daniel C Ellwanger ◽  
Dev Bhatt ◽  
...  

Abstract Background: Elucidation of immune populations with single-cell RNA-seq has greatly benefited the fieldof immunology by deepening the characterization of immune heterogeneity and leading to thediscovery of new subtypes. However, single-cell methods inherently suffer from limitations in therecovery of complete transcriptomes due to the prevalence of cellular and transcriptional dropoutevents. This issue is often compounded by limited sample availability and limited prior knowledge ofheterogeneity, which can confound data interpretation.Results: Here, we systematically benchmarked seven high-throughput single-cell RNA-seq methods. Weprepared 21 libraries under identical conditions of a defined mixture of two human and two murinelymphocyte cell lines, simulating heterogeneity across immune-cell types and cell sizes. We evaluatemethods by their cell recovery rate, library efficiency, sensitivity, and ability to recover expressionsignatures for each cell type. We observed higher mRNA detection sensitivity with the 10x Genomics 5’v1 and 3’ v3 methods. We demonstrate that these methods have fewer drop-out events whichfacilitates the identification of differentially-expressed genes and improves the concordance of singlecellprofiles to immune bulk RNA-seq signatures.Conclusion: Overall, our characterization of immune cell mixtures provides useful metrics, which canguide selection of a high-throughput single-cell RNA-seq method for profiling more complex immunecellheterogeneity usually found in vivo.


Sign in / Sign up

Export Citation Format

Share Document