Accounting for Errors in Data Improves Divergence Time Estimates in Single-cell Cancer Evolution

Single-cell sequencing provides a new way to explore the evolutionary history of cells. Compared to traditional bulk sequencing, where a population of heterogeneous cells is pooled to form a single observation, single-cell sequencing isolates and amplifies genetic material from individual cells, thereby preserving the information about the origin of the sequences. However, single-cell data are more error-prone than bulk sequencing data due to the limited genomic material available per cell. Here, we present error and mutation models for evolutionary inference of single-cell data within a mature and extensible Bayesian framework, BEAST2. Our framework enables integration with biologically informative models such as relaxed molecular clocks and population dynamic models. Our simulations show that modeling errors increase the accuracy of relative divergence times and substitution parameters. We reconstruct the phylogenetic history of a colorectal cancer patient and a healthy patient from single-cell DNA sequencing data. We find that the estimated times of terminal splitting events are shifted forward in time compared to models which ignore errors. We observed that not accounting for errors can overestimate the phylogenetic diversity in single-cell DNA sequencing data. We estimate that 30–50% of the apparent diversity can be attributed to error. Our work enables a full Bayesian approach capable of accounting for errors in the data within the integrative Bayesian software framework BEAST2.     

Hunting down the dominating subclone of cancer stem cells as a potential new therapeutic target in multiple myeloma: An artificial intelligence perspective

The development of single-cell subclones, which can rapidly switch from dormant to dominant subclones, occur in the natural pathophysiology of multiple myeloma(MM) but is often "pressed" by the standard treatment of MM. These emerging subclones present a challenge, providing reservoirs for chemoresistant mutations. Technological advancement is required to track MM subclonal changes, as understanding MM's mechanism of evolution at the cellular level can prompt the development of new targeted ways of treating this disease. Current methods to study the evolution of subclones in MM rely on technologies capable of phenotypically and genotypically characterizing plasma cells, which include immunohistochemistry, flow cytometry, or cytogenetics. Still, all of these technologies may be limited by the sensitivity for picking up rare events. In contrast, more incisive methods such as RNA sequencing, comparative genomic hybridization, or whole-genome sequencing are not yet commonly used in clinical practice. Here we introduce the epidemiological diagnosis and prognosis of MM and review current methods for evaluating MM subclone evolution, such as minimal residual disease/multiparametric flow cytometry/next-generation sequencing, and their respective advantages and disadvantages. In addition, we propose our new single-cell method of evaluation to understand MM's mechanism of evolution at the molecular and cellular level and to prompt the development of new targeted ways of treating this disease, which has a broad prospect.     

Comparison of Multiple Displacement Amplification (MDA) and Multiple Annealing and Looping-Based Amplification Cycles (MALBAC) in Single-Cell Sequencing

Single-cell sequencing promotes our understanding of the heterogeneity of cellular populations, including the haplotypes and genomic variability among different generation of cells. Whole-genome amplification is crucial to generate sufficient DNA fragments for single-cell sequencing projects. Using sequencing data from single sperms, we quantitatively compare two prevailing amplification methods that extensively applied in single-cell sequencing, multiple displacement amplification (MDA) and multiple annealing and looping-based amplification cycles (MALBAC). Our results show that MALBAC, as a combination of modified MDA and tweaked PCR, has a higher level of uniformity, specificity and reproducibility.     

Clonal Architecture of Secondary Acute Myeloid Leukemia Defined by Single-Cell Sequencing

Next-generation sequencing has been used to infer the clonality of heterogeneous tumor samples. These analyses yield specific predictions—the population frequency of individual clones, their genetic composition, and their evolutionary relationships—which we set out to test by sequencing individual cells from three subjects diagnosed with secondary acute myeloid leukemia, each of whom had been previously characterized by whole genome sequencing of unfractionated tumor samples. Single-cell mutation profiling strongly supported the clonal architecture implied by the analysis of bulk material. In addition, it resolved the clonal assignment of single nucleotide variants that had been initially ambiguous and identified areas of previously unappreciated complexity. Accordingly, we find that many of the key assumptions underlying the analysis of tumor clonality by deep sequencing of unfractionated material are valid. Furthermore, we illustrate a single-cell sequencing strategy for interrogating the clonal relationships among known variants that is cost-effective, scalable, and adaptable to the analysis of both hematopoietic and solid tumors, or any heterogeneous population of cells.     

单细胞转录组测序在药物成瘾研究中的应用

药物成瘾是复杂的中枢神经系统疾病，相关基础与临床研究均证实药物成瘾的神经机制及神经环路在成瘾行为形成的不同阶段逐渐发生改变。利用全基因组关联研究、全基因组测序、全外显子测序或高通量转录组测序等技术的组学研究对包括药物成瘾在内的精神疾病遗传的脆弱性进行了深入研究。上述单核苷酸多态性检测技术或测序技术主要预测疾病的遗传风险位点。然而，许多中枢神经系统疾病的发生与环境因素密切相关，而且在疾病发展的不同阶段，相关基因的表达存在脑区特异性的细胞异质性信息。因此，传统研究对发病机制的解释存在一定的局限性。单细胞转录组测序技术是针对单个细胞进行转录水平的测定，规避了传统测序对细胞群体平均转录水平检测的缺点，可以定量描述细胞异质性。近年来，单细胞转录测序技术在神经精神科学研究中的应用逐渐受到关注，本文总结了该技术在神经科学研究中的重要应用，并以药物成瘾为例，重点阐述说明其在中枢神经系统疾病中的应用价值。     

scLM: Automatic Detection of Consensus Gene Clusters Across Multiple Single-cell Datasets

In gene expression profiling studies, including single-cell RNA sequencing(sc RNA-seq)analyses, the identification and characterization of co-expressed genes provides critical information on cell identity and function. Gene co-expression clustering in sc RNA-seq data presents certain challenges. We show that commonly used methods for single-cell data are not capable of identifying co-expressed genes accurately, and produce results that substantially limit biological expectations of co-expressed genes. Herein, we present single-cell Latent-variable Model(sc LM), a gene coclustering algorithm tailored to single-cell data that performs well at detecting gene clusters with significant biologic context. Importantly, sc LM can simultaneously cluster multiple single-cell datasets, i.e., consensus clustering, enabling users to leverage single-cell data from multiple sources for novel comparative analysis. sc LM takes raw count data as input and preserves biological variation without being influenced by batch effects from multiple datasets. Results from both simulation data and experimental data demonstrate that sc LM outperforms the existing methods with considerably improved accuracy. To illustrate the biological insights of sc LM, we apply it to our in-house and public experimental sc RNA-seq datasets. sc LM identifies novel functional gene modules and refines cell states, which facilitates mechanism discovery and understanding of complex biosystems such as cancers. A user-friendly R package with all the key features of the sc LM method is available at https://github.com/QSong-github/sc LM.     

Methods for predicting single-cell miRNA in breast cancer

It has been demonstrated that miRNAs are involved in many biological processes including cell proliferation and differentiation, apoptosis, and stress responses. Although single-cell RNA sequencing technology is prevailing nowadays, it still remains challenging in quantifying miRNA at the single-cell level. Herein, we present the computational methods to infer the single-cell miRNA expression level using its target gene abundances. Firstly, we developed an enrichment-based approach in estimating miRNA expression considering miRNA-mRNA regulation information and miRNA-mRNA correlation signal captured from existing TCGA datasets. Further efforts were made to infer the miRNA expression with machine learning models. The methods were applied to compare the accuracy and robustness with the simulated single-cell data. Finally, we applied the method in single-cell RNA-seq triple negative breast cancer (TNBC) patients to further discover miRNA marker at the single-cell level for the malignant cells. Our tool is available online at: https://github.com/ChengkuiZhao/Single-cell-miRNA-prediction.     

Revealing allele-specific gene expression by single-cell transcriptomics

Single-cell sequencing has emerged as a revolutionary method that reveals biological processes with unprecedented resolution and scale, and has already greatly impacted biology and medicine. To investigate processes such as alternative splicing, novel exon detection and allele-specific expression (ASE), full-length based single-cell RNA-seq methods are required for broad sequence coverage and single nucleotide polymorphism (SNP) identification. In this review, we revisit recent achievements from studies that used single-cell RNA-seq to advance our understanding of ASE in the context of both autosomal and X-chromosome genes. We also recapitulate useful bioinformatic tools developed to identify haplotype phase.     

A statistical framework for multiparameter analysis at the single-cell level

Phenotypic characterization of individual cells provides crucial insights into intercellular heterogeneity and enables access to information that is unavailable from ensemble averaged, bulk cell analyses. Single-cell studies have attracted significant interest in recent years and spurred the development of a variety of commercially available and research-grade technologies. To quantify cell-to-cell variability of cell populations, we have developed an experimental platform for real-time measurements of oxygen consumption (OC) kinetics at the single-cell level. Unique challenges inherent to these single-cell measurements arise, and no existing data analysis methodology is available to address them. Here we present a data processing and analysis method that addresses challenges encountered with this unique type of data in order to extract biologically relevant information. We applied the method to analyze OC profiles obtained with single cells of two different cell lines derived from metaplastic and dysplastic human Barrett's esophageal epithelium. In terms of method development, three main challenges were considered for this heterogeneous dynamic system: (i) high levels of noise, (ii) the lack of a priori knowledge of single-cell dynamics, and (iii) the role of intercellular variability within and across cell types. Several strategies and solutions to address each of these three challenges are presented. The features such as slopes, intercepts, breakpoint or change-point were extracted for every OC profile and compared across individual cells and cell types. The results demonstrated that the extracted features facilitated exposition of subtle differences between individual cells and their responses to cell-cell interactions. With minor modifications, this method can be used to process and analyze data from other acquisition and experimental modalities at the single-cell level, providing a valuable statistical framework for single-cell analysis.     

Diversity of Gene Expression in Hepatocellular Carcinoma Cells

Understanding tumor diversity has been a long-lasting and challenging question for researchers in the field of cancer heterogeneity or tumor evolution. Studies have reported that compared to normal cells, there is a higher genetic diversity in tumor cells, while higher genetic diversity is associated with higher progression risks of tumor. We thus hypothesized that tumor diversity also holds true at the gene expression level. To test this hypothesis, we used t-test to compare the means of Simpson's diversity index for gene expression(SDIG) between tumor and non-tumor samples.We found that the mean SDIG in tumor tissues is significantly higher than that in the non-tumor or normal tissues(P 0.05) for most datasets. We also combined microarrays and next-generation sequencing data for validation. This cross-platform and cross-experimental validation greatly increased the reliability of our results.     

Cancer Cell Analyses at the Single Cell-Level Using Electroactive Microwell Array Device

Circulating tumor cells (CTCs), shed from primary tumors and disseminated into peripheral blood, are playing a major role in metastasis. Even after isolation of CTCs from blood, the target cells are mixed with a population of other cell types. Here, we propose a new method for analyses of cell mixture at the single-cell level using a microfluidic device that contains arrayed electroactive microwells. Dielectrophoretic (DEP) force, induced by the electrodes patterned on the bottom surface of the microwells, allows efficient trapping and stable positioning of single cells for high-throughput biochemical analyses. We demonstrated that various on-chip analyses including immunostaining, viability/apoptosis assay and fluorescent in situ hybridization (FISH) at the single-cell level could be conducted just by applying specific reagents for each assay. Our simple method should greatly help discrimination and analysis of rare cancer cells among a population of blood cells.     

RNA sequencing: new technologies and applications in cancer research

Acute myeloid leukemia (AML) is a fatal hematopoietic malignancy and has a prognosis that varies with its genetic complexity. However, there has been no appropriate integrative analysis on the hierarchy of different AML subtypes. Using Microwell-seq, a high-throughput single-cell mRNA sequencing platform, we analyzed the cellular hierarchy of bone marrow samples from 40 patients and 3 healthy donors. We also used single-cell single-molecule real-time (SMRT) sequencing to investigate the clonal heterogeneity of AML cells. From the integrative analysis of 191727 AML cells, we established a single-cell AML landscape and identified an AML progenitor cell cluster with novel AML markers. Patients with ribosomal protein high progenitor cells had a low remission rate. We deduced two types of AML with diverse clinical outcomes. We traced mitochondrial mutations in the AML landscape by combining Microwell-seq with SMRT sequencing. We propose the existence of a phenotypic “cancer attractor” that might help to define a common phenotype for AML progenitor cells. Finally, we explored the potential drug targets by making comparisons between the AML landscape and the Human Cell Landscape. We identified a key AML progenitor cell cluster. A high ribosomal protein gene level indicates the poor prognosis. We deduced two types of AML and explored the potential drug targets. Our results suggest the existence of a cancer attractor.     

单个肿瘤细胞测序研究进展

肿瘤是机体在各种致癌因素作用下,遗传物质受损导致特定细胞失去对正常生长调控,引起其克隆性异常增生而形成的恶性赘生物。随着单细胞分离技术及单细胞测序技术日益成熟,单个肿瘤细胞全基因组测序已成为肿瘤研究的一个崭新的领域。该文就对单个肿瘤细胞的获取及单个肿瘤细胞测序研究进展进行综述。