A Highly Scalable Method for Joint Whole-Genome Sequencing and Gene-Expression Profiling of Single Cells.

To address how genetic variation alters gene expression in complex cell mixtures, we developed direct nuclear tagmentation and RNA sequencing (DNTR-seq), which enables whole-genome and mRNA sequencing jointly in single cells. DNTR-seq readily identified minor subclones within leukemia patients. In a large-scale DNA damage screen, DNTR-seq was used to detect regions under purifying selection and identified genes where mRNA abundance was resistant to copy-number alteration, suggesting strong genetic compensation. mRNA sequencing (mRNA-seq) quality equals RNA-only methods, and the low positional bias of genomic libraries allowed detection of sub-megabase aberrations at ultra-low coverage. Each cell library is individually addressable and can be re-sequenced at increased depth, allowing multi-tiered study designs. Additionally, the direct tagmentation protocol enables coverage-independent estimation of ploidy, which can be used to identify cell singlets. Thus, DNTR-seq directly links each cell's state to its corresponding genome at scale, enabling routine analysis of heterogeneous tumors and other complex tissues.

Toward the functional interpretation of somatic structural variations: bulk- and single-cell approaches.

Structural variants (SVs) are genomic rearrangements that can take many different forms such as copy number alterations, inversions and translocations. During cell development and aging, somatic SVs accumulate in the genome with potentially neutral, deleterious or pathological effects. Generation of somatic SVs is a key mutational process in cancer development and progression. Despite their importance, the detection of somatic SVs is challenging, making them less studied than somatic single-nucleotide variants. In this review, we summarize recent advances in whole-genome sequencing (WGS)-based approaches for detecting somatic SVs at the tissue and single-cell levels and discuss their advantages and limitations. First, we describe the state-of-the-art computational algorithms for somatic SV calling using bulk WGS data and compare the performance of somatic SV detectors in the presence or absence of a matched-normal control. We then discuss the unique features of cutting-edge single-cell-based techniques for analyzing somatic SVs. The advantages and disadvantages of bulk and single-cell approaches are highlighted, along with a discussion of their sensitivity to copy-neutral SVs, usefulness for functional inferences and experimental and computational costs. Finally, computational approaches for linking somatic SVs to their functional readouts, such as those obtained from single-cell transcriptome and epigenome analyses, are illustrated, with a discussion of the promise of these approaches in health and diseases.

Single-cell whole-genome sequencing reveals the functional landscape of somatic mutations in B lymphocytes across the human lifespan.

Accumulation of mutations in somatic cells has been implicated as a cause of aging since the 1950s. However, attempts to establish a causal relationship between somatic mutations and aging have been constrained by the lack of methods to directly identify mutational events in primary human tissues. Here we provide genome-wide mutation frequencies and spectra of human B lymphocytes from healthy individuals across the entire human lifespan using a highly accurate single-cell whole-genome sequencing method. We found that the number of somatic mutations increases from <500 per cell in newborns to >3,000 per cell in centenarians. We discovered mutational hotspot regions, some of which, as expected, were located at Ig genes associated with somatic hypermutation (SHM). B cell-specific mutation signatures associated with development, aging, or SHM were found. The SHM signature strongly correlated with the signature found in human B cell tumors, indicating that potential cancer-causing events are already present even in B cells of healthy individuals. We also identified multiple mutations in sequence features relevant to cellular function (i.e., transcribed genes and gene regulatory regions). Such mutations increased significantly during aging, but only at approximately one-half the rate of the genome average, indicating selection against mutations that impact B cell function. This full characterization of the landscape of somatic mutations in human B lymphocytes indicates that spontaneous somatic mutations accumulating with age can be deleterious and may contribute to both the increased risk for leukemia and the functional decline of B lymphocytes in the elderly.

Multi-region and single-cell sequencing reveal variable genomic heterogeneity in rectal cancer.

BACKGROUND: Colorectal cancer is a heterogeneous group of malignancies with complex molecular subtypes. While colon cancer has been widely investigated, studies on rectal cancer are very limited. Here, we performed multi-region whole-exome sequencing and single-cell whole-genome sequencing to examine the genomic intratumor heterogeneity (ITH) of rectal tumors. METHODS: We sequenced nine tumor regions and 88 single cells from two rectal cancer patients with tumors of the same molecular classification and characterized their mutation profiles and somatic copy number alterations (SCNAs) at the multi-region and the single-cell levels. RESULTS: A variable extent of genomic heterogeneity was observed between the two patients, and the degree of ITH increased when analyzed on the single-cell level. We found that major SCNAs were early events in cancer development and inherited steadily. Single-cell sequencing revealed mutations and SCNAs which were hidden in bulk sequencing. In summary, we studied the ITH of rectal cancer at regional and single-cell resolution and demonstrated that variable heterogeneity existed in two patients. The mutational scenarios and SCNA profiles of two patients with treatment naïve from the same molecular subtype are quite different. CONCLUSIONS: Our results suggest each tumor possesses its own architecture, which may result in different diagnosis, prognosis, and drug responses. Remarkable ITH exists in the two patients we have studied, providing a preliminary impression of ITH in rectal cancer.

Somatic mutations in aging and disease.

Time always leaves its mark, and our genome is no exception. Mutations in the genome of somatic cells were first hypothesized to be the cause of aging in the 1950s, shortly after the molecular structure of DNA had been described. Somatic mutation theories of aging are based on the fact that mutations in DNA as the ultimate template for all cellular functions are irreversible. However, it took until the 1990s to develop the methods to test if DNA mutations accumulate with age in different organs and tissues and estimate the severity of the problem. By now, numerous studies have documented the accumulation of somatic mutations with age in normal cells and tissues of mice, humans, and other animals, showing clock-like mutational signatures that provide information on the underlying causes of the mutations. In this review, we will first briefly discuss the recent advances in next-generation sequencing that now allow quantitative analysis of somatic mutations. Second, we will provide evidence that the mutation rate differs between cell types, with a focus on differences between germline and somatic mutation rate. Third, we will discuss somatic mutational signatures as measures of aging, environmental exposure, and activities of DNA repair processes. Fourth, we will explain the concept of clonally amplified somatic mutations, with a focus on clonal hematopoiesis. Fifth, we will briefly discuss somatic mutations in the transcriptome and in our other genome, i.e., the genome of mitochondria. We will end with a brief discussion of a possible causal contribution of somatic mutations to the aging process.

Negative Selection Allows DNA Mismatch Repair-Deficient Mouse Fibroblasts In Vitro to Tolerate High Levels of Somatic Mutations.

Substantial numbers of somatic mutations have been found to accumulate with age in different human tissues. Clonal cellular amplification of some of these mutations can cause cancer and other diseases. However, it is as yet unclear if and to what extent an increased burden of random mutations can affect cellular function without clonal amplification. We tested this in cell culture, which avoids the limitation that an increased mutation burden in vivo typically leads to cancer. We performed single-cell whole-genome sequencing of primary fibroblasts from DNA mismatch repair (MMR) deficient Msh2-/- mice and littermate control animals after long-term passaging. Apart from analyzing somatic mutation burden we analyzed clonality, mutational signatures, and hotspots in the genome, characterizing the complete landscape of somatic mutagenesis in normal and MMR-deficient mouse primary fibroblasts during passaging. While growth rate of Msh2-/- fibroblasts was not significantly different from the controls, the number of de novo single-nucleotide variants (SNVs) increased linearly up until at least 30,000 SNVs per cell, with the frequency of small insertions and deletions (INDELs) plateauing in the Msh2-/- fibroblasts to about 10,000 INDELS per cell. We provide evidence for negative selection and large-scale mutation-driven population changes, including significant clonal expansion of preexisting mutations and widespread cell-strain-specific hotspots. Overall, our results provide evidence that increased somatic mutation burden drives significant cell evolutionary changes in a dynamic cell culture system without significant effects on growth. Since similar selection processes against mutations preventing organ and tissue dysfunction during aging are difficult to envision, these results suggest that increased somatic mutation burden can play a causal role in aging and diseases other than cancer.

Single cell analyses of cancer cells identified two regulatorily and functionally distinct categories in differentially expressed genes among tumor subclones.

To explore the feature of cancer cells and tumor subclones, we analyzed 101,065 single-cell transcriptomes from 12 colorectal cancer (CRC) patients and 92 single cell genomes from one of these patients. We found cancer cells, endothelial cells and stromal cells in tumor tissue expressed much more genes and had stronger cell-cell interactions than their counterparts in normal tissue. We identified copy number variations (CNVs) in each cancer cell and found correlation between gene copy number and expression level in cancer cells at single cell resolution. Analysis of tumor subclones inferred by CNVs showed accumulation of mutations in each tumor subclone along lineage trajectories. We found differentially expressed genes (DEGs) between tumor subclones had two populations: DEGCNV and DEGreg. DEGCNV, showing high CNV-expression correlation and whose expression differences depend on the differences of CNV level, enriched in housekeeping genes and cell adhesion associated genes. DEGreg, showing low CNV-expression correlation and mainly in low CNV variation regions and regions without CNVs, enriched in cytokine signaling genes. Furthermore, cell-cell communication analyses showed that DEGCNV tends to involve in cell-cell contact while DEGreg tends to involve in secreted signaling, which further support that DEGCNV and DEGreg are two regulatorily and functionally distinct categories.

Single-cell DNA replication dynamics in genomically unstable cancers.

Dysregulated DNA replication is both a cause and a consequence of aneuploidy, yet the dynamics of DNA replication in aneuploid cell populations remains understudied. We developed a new method, PERT, for inferring cell-specific DNA replication states from single-cell whole genome sequencing, and investigated clone-specific DNA replication dynamics in >50,000 cells obtained from a collection of aneuploid and clonally heterogeneous cell lines, xenografts and primary cancer tissues. Clone replication timing (RT) profiles correlated with future copy number changes in serially passaged cell lines. Cell type was the strongest determinant of RT heterogeneity, while whole genome doubling and mutational process were associated with accumulation of late S-phase cells and weaker RT associations. Copy number changes affecting chromosome X had striking impact on RT, with loss of the inactive X allele shifting replication earlier, and loss of inactive Xq resulting in reactivation of Xp. Finally, analysis of time series xenografts illustrate how cell cycle distributions approximate clone proliferation, recapitulating expected relationships between proliferation and fitness in treatment-naive and chemotherapeutic contexts.

Integration of multiple lineage measurements from the same cell reconstructs parallel tumor evolution.

Organoid evolution models complemented with integrated single-cell sequencing technology provide a powerful platform to characterize intra-tumor heterogeneity (ITH) and tumor evolution. Here, we conduct a parallel evolution experiment to mimic the tumor evolution process by evolving a colon cancer organoid model over 100 generations, spanning 6 months in time. We use single-cell whole-genome sequencing (WGS) in combination with viral lineage tracing at 12 time points to simultaneously monitor clone size, CNV states, SNV states, and viral lineage barcodes for 1,641 single cells. We integrate these measurements to construct clonal evolution trees with high resolution. We characterize the order of events in which chromosomal aberrations occur and identify aberrations that recur multiple times within the same tumor sub-population. We observe recurrent sequential loss of chromosome 4 after loss of chromosome 18 in four unique tumor clones. SNVs and CNVs identified in our organoid experiments are also frequently reported in colorectal carcinoma samples, and out of 334 patients with chromosome 18 loss in a Memorial Sloan Kettering colorectal cancer cohort, 99 (29.6%) also harbor chromosome 4 loss. Our study reconstructs tumor evolution in a colon cancer organoid model at high resolution, demonstrating an approach to identify potentially clinically relevant genomic aberrations in tumor evolution.

Whole-Chromosome Karyotyping of Fetal Nucleated Red Blood Cells Using the Ion Proton Sequencing Platform.

The current gold standard for the definitive diagnosis of fetal aneuploidy uses either chorionic villus sampling (CVS) or amniocentesis, both of which are which are invasive procedures carrying a procedure-related risk of miscarriage of up to 0.1-0.2%. Non-invasive prenatal diagnosis using fetal nucleated red blood cells (FNRBCs) isolated from maternal peripheral venous blood would remove this risk of miscarriage since these cells can be isolated from the mother's blood. We aimed to detect whole-chromosome aneuploidies from single nucleated fetal red blood cells using whole-genome amplification followed by massively parallel sequencing performed on a semiconductor sequencing platform. Twenty-six single cells were picked from the placental villi of twelve patients thought to have a normal fetal genotype and who were undergoing elective first-trimester surgical termination of pregnancy. Following karyotyping, it was subsequently found that two of these cases were also abnormal (one trisomy 15 and one mosaic genotype). One single cell from chorionic villus samples for two patients carrying a fetus with trisomy 21 and two single cells from women carrying fetuses with T18 were also picked. Pooled libraries were sequenced on the Ion Proton and data were analysed using Ion Reporter software. We correctly classified fetal genotype in all 24 normal cells, as well as the 2 T21 cells, the 2 T18 cells, and the two T15 cells. The two cells picked from the fetus with a mosaic result by CVS were classified as unaffected, suggesting that this was a case of confined placental mosaicism. Fetal sex was correctly assigned in all cases. We demonstrated that semiconductor sequencing using commercially available software for data analysis can be achieved for the non-invasive prenatal diagnosis of whole-chromosome aneuploidy with 100% accuracy.

Somatic copy number variant load in neurons of healthy controls and Alzheimer's disease patients.

The possible role of somatic copy number variations (CNVs) in Alzheimer's disease (AD) aetiology has been controversial. Although cytogenetic studies suggested increased CNV loads in AD brains, a recent single-cell whole-genome sequencing (scWGS) experiment, studying frontal cortex brain samples, found no such evidence. Here we readdressed this issue using low-coverage scWGS on pyramidal neurons dissected via both laser capture microdissection (LCM) and fluorescence activated cell sorting (FACS) across five brain regions: entorhinal cortex, temporal cortex, hippocampal CA1, hippocampal CA3, and the cerebellum. Among reliably detected somatic CNVs identified in 1301 cells obtained from the brains of 13 AD patients and 7 healthy controls, deletions were more frequent compared to duplications. Interestingly, we observed slightly higher frequencies of CNV events in cells from AD compared to similar numbers of cells from controls (4.1% vs. 1.4%, or 0.9% vs. 0.7%, using different filtering approaches), although the differences were not statistically significant. On the technical aspects, we observed that LCM-isolated cells show higher within-cell read depth variation compared to cells isolated with FACS. To reduce within-cell read depth variation, we proposed a principal component analysis-based denoising approach that significantly improves signal-to-noise ratios. Lastly, we showed that LCM-isolated neurons in AD harbour slightly more read depth variability than neurons of controls, which might be related to the reported hyperploid profiles of some AD-affected neurons.

Predicting the Number of Bases to Attain Sufficient Coverage in High-Throughput Sequencing Experiments.

For many types of high-throughput sequencing experiments, success in downstream analysis depends on attaining sufficient coverage for individual positions in the genome. For example, when identifying single-nucleotide variants de novo, the number of reads supporting a particular variant call determines our confidence in that variant call. If sequenced reads are distributed uniformly along the genome, the coverage of a nucleotide position is easily approximated by a Poisson distribution, with rate equal to average sequencing depth. Unfortunately, as has become well known, high-throughput sequencing data are never uniform. The numerous factors contributing to variation in coverage have resisted attempts at direct modeling and change along with minor adjustments in the underlying technology. We propose a new nonparametric method to predict the portion of a genome that will attain some specified minimum coverage, as a function of sequencing effort, using information from a shallow sequencing experiment from the same library. Simulations show our approach performs well under an array of distributional assumptions that deviate from uniformity. We applied this approach to estimate coverage at varying depths in single-cell whole-genome sequencing data from multiple protocols. These resulted in highly accurate predictions, demonstrating the effectiveness of our approach in analyzing complexity of sequencing libraries and optimizing design of sequencing experiments.

SCCNV: A Software Tool for Identifying Copy Number Variation From Single-Cell Whole-Genome Sequencing.

Identification of de novo copy number variations (CNVs) across the genome in single cells requires single-cell whole-genome amplification (WGA) and sequencing. Although many experimental protocols of amplification methods have been developed, all suffer from uneven distribution of read depth across the genome after sequencing of DNA amplicons, which constrains the usage of conventional CNV calling methodologies. Here, we present SCCNV, a software tool for detecting CNVs from whole genome-amplified single cells. SCCNV is a read-depth based approach with adjustment for the WGA bias. We demonstrate its performance by analyzing data obtained with most of the single-cell amplification methods that have been employed for CNV analysis, including DOP-PCR, MDA, MALBAC, and LIANTI. SCCNV is freely available at https://github.com/biosinodx/SCCNV.

Quantification of Aneuploidy in Mammalian Systems.

High-throughput next generation sequencing karyotyping has emerged as a powerful tool for the detection of genomic heterogeneity in normal tissues and cancers. Here we describe a single-cell whole genome sequencing (scWGS) platform to assess whole-chromosome aneuploidy, structural aneuploidies involving only chromosome fragments and more local small copy number alterations in individual cells. We provide a detailed protocol for the isolation, library preparation, low coverage sequencing and data analysis of single cells. Since our approach does not involve a whole-genome preamplification step, our method allows for acquisition of reliable high-resolution single-cell copy number profiles. Moreover, the protocol allows multiplexing of 384 single-cell libraries in one sequencing run, thereby significantly reducing sequencing costs and can be completed in 3-4 days starting from single cell isolation to analysis of sequencing data.