Chenghang Zong

Genome-wide detection of single-nucleotide and copy-number variations of a single human cell.

Single-Cell Sequencing With the rapid progress in sequencing technologies, single-cell sequencing is now possible, promising insight into how cell-to-cell heterogeneity affects biological behavior. Achieving adequate genome coverage remains a challenge because single-cell sequencing relies on genome amplification that is prone to sequence bias. Zong et al. (p. 1622) report a new amplification method: multiple annealing and looping-based amplification cycles that allowed 93% genome coverage for a human cell. This coverage facilitated accurate detection of point mutations and copy number variations. Lu et al. (p. 1627) used the method to sequence 99 sperm cells from a single individual. Mapping the meiotic crossovers revealed a nonrandom distribution with a reduced recombination rate near transcription start sites. A whole-genome amplification method with reduced bias compares a single cell with its descendants. Kindred cells can have different genomes because of dynamic changes in DNA. Single-cell sequencing is needed to characterize these genomic differences but has been hindered by whole-genome amplification bias, resulting in low genome coverage. Here, we report on a new amplification method—multiple annealing and looping-based amplification cycles (MALBAC)—that offers high uniformity across the genome. Sequencing MALBAC-amplified DNA achieves 93% genome coverage ≥1x for a single human cell at 25x mean sequencing depth. We detected digitized copy-number variations (CNVs) of a single cancer cell. By sequencing three kindred cells, we were able to identify individual single-nucleotide variations (SNVs), with no false positives detected. We directly measured the genome-wide mutation rate of a cancer cell line and found that purine-pyrimidine exchanges occurred unusually frequently among the newly acquired SNVs.

Reproducible copy number variation patterns among single circulating tumor cells of lung cancer patients

Significance In a few milliliters of blood from a cancer patient, one can isolate a few circulating tumor cells (CTCs). Originating from the primary tumor, CTCs seed metastases, which account for the majority of cancer-related deaths. We demonstrate the analyses of the whole genome of single CTCs, which are highly needed for personalized treatment. We discovered that copy number variations (CNVs), one of the major genomic variations, are specific to cancer types, reproducible from cell to cell, and even from patient to patient. We hypothesize that CNVs at certain genomic loci are selected for and lead to metastasis. Our work shows the prospect of noninvasive CTC-based cancer diagnostics. Circulating tumor cells (CTCs) enter peripheral blood from primary tumors and seed metastases. The genome sequencing of CTCs could offer noninvasive prognosis or even diagnosis, but has been hampered by low single-cell genome coverage of scarce CTCs. Here, we report the use of the recently developed multiple annealing and looping-based amplification cycles for whole-genome amplification of single CTCs from lung cancer patients. We observed characteristic cancer-associated single-nucleotide variations and insertions/deletions in exomes of CTCs. These mutations provided information needed for individualized therapy, such as drug resistance and phenotypic transition, but were heterogeneous from cell to cell. In contrast, every CTC from an individual patient, regardless of the cancer subtypes, exhibited reproducible copy number variation (CNV) patterns, similar to those of the metastatic tumor of the same patient. Interestingly, different patients with the same lung cancer adenocarcinoma (ADC) shared similar CNV patterns in their CTCs. Even more interestingly, patients of small-cell lung cancer have CNV patterns distinctly different from those of ADC patients. Our finding suggests that CNVs at certain genomic loci are selected for the metastasis of cancer. The reproducibility of cancer-specific CNVs offers potential for CTC-based cancer diagnostics.

General properties of transcriptional time series in Escherichia coli

Gene activity is described by the time series of discrete, stochastic mRNA production events. This transcriptional time series shows intermittent, bursty behavior. One consequence of this temporal intricacy is that gene expression can be tuned by varying different features of the time series. Here we quantify copy-number statistics of mRNA from 20 Escherichia coli promoters using single-molecule fluorescence in situ hybridization in order to characterize the general properties of these transcriptional time series. We find that the degree of burstiness is correlated with gene expression level but is largely independent of other parameters of gene regulation. The observed behavior can be explained by the underlying variation in the duration of bursting events. Using Shannons mutual information function, we estimate the mutual information transmitted between an outside stimulus, such as the extracellular concentration of inducer molecules, and intracellular levels of mRNA. This suggests that the outside stimulus transmits information reflected in the properties of transcriptional time series.

Probing meiotic recombination and aneuploidy of single sperm cells by whole-genome sequencing.

Single-Cell Sequencing With the rapid progress in sequencing technologies, single-cell sequencing is now possible, promising insight into how cell-to-cell heterogeneity affects biological behavior. Achieving adequate genome coverage remains a challenge because single-cell sequencing relies on genome amplification that is prone to sequence bias. Zong et al. (p. 1622) report a new amplification method: multiple annealing and looping-based amplification cycles that allowed 93% genome coverage for a human cell. This coverage facilitated accurate detection of point mutations and copy number variations. Lu et al. (p. 1627) used the method to sequence 99 sperm cells from a single individual. Mapping the meiotic crossovers revealed a nonrandom distribution with a reduced recombination rate near transcription start sites. A whole-genome amplification method with reduced bias yields a personal meiotic recombination map. Meiotic recombination creates genetic diversity and ensures segregation of homologous chromosomes. Previous population analyses yielded results averaged among individuals and affected by evolutionary pressures. We sequenced 99 sperm from an Asian male by using the newly developed amplification method—multiple annealing and looping-based amplification cycles—to phase the personal genome and map recombination events at high resolution, which are nonuniformly distributed across the genome in the absence of selection pressure. The paucity of recombination near transcription start sites observed in individual sperm indicates that such a phenomenon is intrinsic to the molecular mechanism of meiosis. Interestingly, a decreased crossover frequency combined with an increase of autosomal aneuploidy is observable on a global per-sperm basis.

Mutual regulation of tumour vessel normalization and immunostimulatory reprogramming

Blockade of angiogenesis can retard tumour growth, but may also paradoxically increase metastasis. This paradox may be resolved by vessel normalization, which involves increased pericyte coverage, improved tumour vessel perfusion, reduced vascular permeability, and consequently mitigated hypoxia. Although these processes alter tumour progression, their regulation is poorly understood. Here we show that type 1 T helper (TH1) cells play a crucial role in vessel normalization. Bioinformatic analyses revealed that gene expression features related to vessel normalization correlate with immunostimulatory pathways, especially T lymphocyte infiltration or activity. To delineate the causal relationship, we used various mouse models with vessel normalization or T lymphocyte deficiencies. Although disruption of vessel normalization reduced T lymphocyte infiltration as expected, reciprocal depletion or inactivation of CD4+ T lymphocytes decreased vessel normalization, indicating a mutually regulatory loop. In addition, activation of CD4+ T lymphocytes by immune checkpoint blockade increased vessel normalization. TH1 cells that secrete interferon-γ are a major population of cells associated with vessel normalization. Patient-derived xenograft tumours growing in immunodeficient mice exhibited enhanced hypoxia compared to the original tumours in immunocompetent humans, and hypoxia was reduced by adoptive TH1 transfer. Our findings elucidate an unexpected role of TH1 cells in vasculature and immune reprogramming. TH1 cells may be a marker and a determinant of both immune checkpoint blockade and anti-angiogenesis efficacy.

Lysogen stability is determined by the frequency of activity bursts from the fate‐determining gene

The ability of living cells to maintain an inheritable memory of their gene‐expression state is key to cellular differentiation. Bacterial lysogeny serves as a simple paradigm for long‐term cellular memory. In this study, we address the following question: in the absence of external perturbation, how long will a cell stay in the lysogenic state before spontaneously switching away from that state? We show by direct measurement that lysogen stability exhibits a simple exponential dependence on the frequency of activity bursts from the fate‐determining gene, cI. We quantify these gene‐activity bursts using single‐molecule‐resolution mRNA measurements in individual cells, analyzed using a stochastic mathematical model of the gene‐network kinetics. The quantitative relation between stability and gene activity is independent of the fine details of gene regulation, suggesting that a quantitative prediction of cell‐state stability may also be possible in more complex systems.

Effective detection of variation in single-cell transcriptomes using MATQ-seq.

The quantification of transcriptional variation in single cells, particularly within the same cell population, is currently limited by the low sensitivity and high technical noise of single-cell RNA-seq assays. We report multiple annealing and dC-tailing-based quantitative single-cell RNA-seq (MATQ-seq), a highly sensitive and quantitative method for single-cell sequencing of total RNA. By systematically determining technical noise, we show that MATQ-seq captures genuine biological variation between whole transcriptomes of single cells.

Rare event of histone demethylation can initiate singular gene expression of olfactory receptors

Significance In mammals, the sense of odors relies on the peculiar expression pattern of olfactory receptors (ORs). Each single neuron chooses one, and only one, from all ∼1,400 OR genes that are present in a mouse genome. In neurobiology, a long-standing mystery is how such singularity can be achieved. We show theoretically that a simple kinetic scheme of OR activation followed by feedback can be solely responsible for the observed singularity, as long as the two timescales—slow activation by epigenetic modification and fast feedback by transcriptional regulation—are well separated. Our work provides the theoretical underpinning behind the choice of ORs, and demonstrates how the nervous system utilizes the kinetics of epigenetic changes to direct neurogenesis. Mammals sense odors through the gene family of olfactory receptors (ORs). Despite the enormous number of OR genes (∼1,400 in mouse), each olfactory sensory neuron expresses one, and only one, of them. In neurobiology, it remains a long-standing mystery how this singularity can be achieved despite intrinsic stochasticity of gene expression. Recent experiments showed an epigenetic mechanism for maintaining singular OR expression: Once any ORs are activated, their expression inhibits further OR activation by down-regulating a histone demethylase Lsd1 (also known as Aof2 or Kdm1a), an enzyme required for the removal of the repressive histone marker H3K9me3 on OR genes. However, it remains unclear at a quantitative level how singularity can be initiated in the first place. In particular, does a simple activation/feedback scheme suffice to generate singularity? Here we show theoretically that rare events of histone demethylation can indeed produce robust singularity by separating two timescales: slow OR activation by stepwise H3K9me3 demethylation, and fast feedback to turn off Lsd1. Given a typical 1-h response of transcriptional feedback, to achieve the observed extent of singularity (only 2% of neurons express more than one ORs), we predict that OR activation must be as slow as 5–10 d—a timescale compatible with experiments. Our model further suggests H3K9me3-to-H3K9me2 demethylation as an additional rate-limiting step responsible for OR singularity. Our conclusions may be generally applicable to other systems where monoallelic expression is desired, and provide guidelines for the design of a synthetic system of singular expression.

Bone-in-culture array as a platform to model early-stage bone metastases and discover anti-metastasis therapies

The majority of breast cancer models for drug discovery are based on orthotopic or subcutaneous tumours. Therapeutic responses of metastases, especially microscopic metastases, are likely to differ from these tumours due to distinct cancer-microenvironment crosstalk in distant organs. Here, to recapitulate such differences, we established an ex vivo bone metastasis model, termed bone-in-culture array or BICA, by fragmenting mouse bones preloaded with breast cancer cells via intra-iliac artery injection. Cancer cells in BICA maintain features of in vivo bone micrometastases regarding the microenvironmental niche, gene expression profile, metastatic growth kinetics and therapeutic responses. Through a proof-of-principle drug screening using BICA, we found that danusertib, an inhibitor of the Aurora kinase family, preferentially inhibits bone micrometastases. In contrast, certain histone methyltransferase inhibitors stimulate metastatic outgrowth of indolent cancer cells, specifically in the bone. Thus, BICA can be used to investigate mechanisms involved in bone colonization and to rapidly test drug efficacies on bone micrometastases.

Endocrine lineage biases arise in temporally distinct endocrine progenitors during pancreatic morphogenesis

Decoding the molecular composition of individual Ngn3 + endocrine progenitors (EPs) during pancreatic morphogenesis could provide insight into the mechanisms regulating hormonal cell fate. Here, we identify population markers and extensive cellular diversity including four EP subtypes reflecting EP maturation using high-resolution single-cell RNA-sequencing of the e14.5 and e16.5 mouse pancreas. While e14.5 and e16.5 EPs are constantly born and share select genes, these EPs are overall transcriptionally distinct concomitant with changes in the underlying epithelium. As a consequence, e16.5 EPs are not the same as e14.5 EPs: e16.5 EPs have a higher propensity to form beta cells. Analysis of e14.5 and e16.5 EP chromatin states reveals temporal shifts, with enrichment of beta cell motifs in accessible regions at later stages. Finally, we provide transcriptional maps outlining the route progenitors take as they make cell fate decisions, which can be applied to advance the in vitro generation of beta cells.Endocrine progenitors form early in pancreatic development but the diversity of this cell population is unclear. Here, the authors use single cell RNA sequencing of the mouse pancreas at e14.5 and e16.5 to show that endocrine progenitors are temporally distinct and those formed later are more likely to become beta cells