Modeling gene expression cascades during cell state transitions.

During cellular processes such as differentiation or response to external stimuli, cells exhibit dynamic changes in their gene expression profiles. Single-cell RNA sequencing (scRNA-seq) can be used to investigate these dynamic changes. To this end, cells are typically ordered along a pseudotemporal trajectory which recapitulates the progression of cells as they transition from one cell state to another. We infer transcriptional dynamics by modeling the gene expression profiles in pseudotemporally ordered cells using a Bayesian inference approach. This enables ordering genes along transcriptional cascades, estimating differences in the timing of gene expression dynamics, and deducing regulatory gene interactions. Here, we apply this approach to scRNA-seq datasets derived from mouse embryonic forebrain and pancreas samples. This analysis demonstrates the utility of the method to derive the ordering of gene dynamics and regulatory relationships critical for proper cellular differentiation and maturation across a variety of developmental contexts.

Inferring early genetic progression in cancers with unobtainable premalignant disease.

Analysis of premalignant tissue has identified the typical order of somatic events leading to invasive tumors in several cancer types. For other cancers, premalignant tissue is unobtainable, leaving genetic progression unknown. Here, we demonstrate how to infer progression from exome sequencing of primary tumors. Our computational method, PhylogicNDT, recapitulated the previous experimentally determined genetic progression of human papillomavirus-negative (HPV-) head and neck squamous cell carcinoma (HNSCC). We then evaluated HPV+ HNSCC, which lacks premalignant tissue, and uncovered its previously unknown progression, identifying early drivers. We converted relative timing estimates of driver mutations and HPV integration to years before diagnosis based on a clock-like mutational signature. We associated the timing of transitions to aneuploidy with increased intratumor genetic heterogeneity and shorter overall survival. Our approach can establish previously unknown early genetic progression of cancers with unobtainable premalignant tissue, supporting development of experimental models and methods for early detection, interception and prognostication.

Evolutionary history of transformation from chronic lymphocytic leukemia to Richter syndrome.

Richter syndrome (RS) arising from chronic lymphocytic leukemia (CLL) exemplifies an aggressive malignancy that develops from an indolent neoplasm. To decipher the genetics underlying this transformation, we computationally deconvoluted admixtures of CLL and RS cells from 52 patients with RS, evaluating paired CLL-RS whole-exome sequencing data. We discovered RS-specific somatic driver mutations (including IRF2BP2, SRSF1, B2M, DNMT3A and CCND3), recurrent copy-number alterations beyond del(9p21)(CDKN2A/B), whole-genome duplication and chromothripsis, which were confirmed in 45 independent RS cases and in an external set of RS whole genomes. Through unsupervised clustering, clonally related RS was largely distinct from diffuse large B cell lymphoma. We distinguished pathways that were dysregulated in RS versus CLL, and detected clonal evolution of transformation at single-cell resolution, identifying intermediate cell states. Our study defines distinct molecular subtypes of RS and highlights cell-free DNA analysis as a potential tool for early diagnosis and monitoring.

The immune microenvironment shapes transcriptional and genetic heterogeneity in chronic lymphocytic leukemia.

In chronic lymphocytic leukemia (CLL), B-cell receptor signaling, tumor-microenvironment interactions, and somatic mutations drive disease progression. To better understand the intersection between the microenvironment and molecular events in CLL pathogenesis, we integrated bulk transcriptome profiling of paired peripheral blood (PB) and lymph node (LN) samples from 34 patients. Oncogenic processes were upregulated in LN compared with PB and in immunoglobulin heavy-chain variable (IGHV) region unmutated compared with mutated cases. Single-cell RNA sequencing (scRNA-seq) distinguished 3 major cell states: quiescent, activated, and proliferating. The activated subpopulation comprised only 2.2% to 4.3% of the total tumor bulk in LN samples. RNA velocity analysis found that CLL cell fate in LN is unidirectional, starts in the proliferating state, transitions to the activated state, and ends in the quiescent state. A 10-gene signature derived from activated tumor cells was associated with inferior treatment-free survival (TFS) and positively correlated with the proportion of activated CD4+ memory T cells and M2 macrophages in LN. Whole exome sequencing (WES) of paired PB and LN samples showed subclonal expansion in LN in approximately half of the patients. Since mouse models have implicated activation-induced cytidine deaminase in mutagenesis, we compared AICDA expression between cases with and without clonal evolution but did not find a difference. In contrast, the presence of a T-cell inflamed microenvironment in LN was associated with clonal stability. In summary, a distinct minor tumor subpopulation underlies CLL pathogenesis and drives the clinical outcome. Clonal trajectories are shaped by the LN milieu, where T-cell immunity may contribute to suppressing clonal outgrowth. The clinical study is registered at clinicaltrials.gov as NCT00923507.

Enhanced cortical neural stem cell identity through short SMAD and WNT inhibition in human cerebral organoids facilitates emergence of outer radial glial cells.

Cerebral organoids exhibit broad regional heterogeneity accompanied by limited cortical cellular diversity despite the tremendous upsurge in derivation methods, suggesting inadequate patterning of early neural stem cells (NSCs). Here we show that a short and early Dual SMAD and WNT inhibition course is necessary and sufficient to establish robust and lasting cortical organoid NSC identity, efficiently suppressing non-cortical NSC fates, while other widely used methods are inconsistent in their cortical NSC-specification capacity. Accordingly, this method selectively enriches for outer radial glia NSCs, which cyto-architecturally demarcate well-defined outer sub-ventricular-like regions propagating from superiorly radially organized, apical cortical rosette NSCs. Finally, this method culminates in the emergence of molecularly distinct deep and upper cortical layer neurons, and reliably uncovers cortex-specific microcephaly defects. Thus, a short SMAD and WNT inhibition is critical for establishing a rich cortical cell repertoire that enables mirroring of fundamental molecular and cyto-architectural features of cortical development and meaningful disease modelling.

Combined tumor and immune signals from genomes or transcriptomes predict outcomes of checkpoint inhibition in melanoma.

Immune checkpoint blockade (CPB) improves melanoma outcomes, but many patients still do not respond. Tumor mutational burden (TMB) and tumor-infiltrating T cells are associated with response, and integrative models improve survival prediction. However, integrating immune/tumor-intrinsic features using data from a single assay (DNA/RNA) remains underexplored. Here, we analyze whole-exome and bulk RNA sequencing of tumors from new and published cohorts of 189 and 178 patients with melanoma receiving CPB, respectively. Using DNA, we calculate T cell and B cell burdens (TCB/BCB) from rearranged TCR/Ig sequences and find that patients with TMBhigh and TCBhigh or BCBhigh have improved outcomes compared to other patients. By combining pairs of immune- and tumor-expressed genes, we identify three gene pairs associated with response and survival, which validate in independent cohorts. The top model includes lymphocyte-expressed MAP4K1 and tumor-expressed TBX3. Overall, RNA or DNA-based models combining immune and tumor measures improve predictions of melanoma CPB outcomes.

Characterizing genetic intra-tumor heterogeneity across 2,658 human cancer genomes.

Intra-tumor heterogeneity (ITH) is a mechanism of therapeutic resistance and therefore an important clinical challenge. However, the extent, origin, and drivers of ITH across cancer types are poorly understood. To address this, we extensively characterize ITH across whole-genome sequences of 2,658 cancer samples spanning 38 cancer types. Nearly all informative samples (95.1%) contain evidence of distinct subclonal expansions with frequent branching relationships between subclones. We observe positive selection of subclonal driver mutations across most cancer types and identify cancer type-specific subclonal patterns of driver gene mutations, fusions, structural variants, and copy number alterations as well as dynamic changes in mutational processes between subclonal expansions. Our results underline the importance of ITH and its drivers in tumor evolution and provide a pan-cancer resource of comprehensively annotated subclonal events from whole-genome sequencing data.

Massively parallel single-cell mitochondrial DNA genotyping and chromatin profiling.

Natural mitochondrial DNA (mtDNA) mutations enable the inference of clonal relationships among cells. mtDNA can be profiled along with measures of cell state, but has not yet been combined with the massively parallel approaches needed to tackle the complexity of human tissue. Here, we introduce a high-throughput, droplet-based mitochondrial single-cell assay for transposase-accessible chromatin with sequencing (scATAC-seq), a method that combines high-confidence mtDNA mutation calling in thousands of single cells with their concomitant high-quality accessible chromatin profile. This enables the inference of mtDNA heteroplasmy, clonal relationships, cell state and accessible chromatin variation in individual cells. We reveal single-cell variation in heteroplasmy of a pathologic mtDNA variant, which we associate with intra-individual chromatin variability and clonal evolution. We clonally trace thousands of cells from cancers, linking epigenomic variability to subclonal evolution, and infer cellular dynamics of differentiating hematopoietic cells in vitro and in vivo. Taken together, our approach enables the study of cellular population dynamics and clonal properties in vivo.

The evolutionary history of 2,658 cancers.

Cancer develops through a process of somatic evolution1,2. Sequencing data from a single biopsy represent a snapshot of this process that can reveal the timing of specific genomic aberrations and the changing influence of mutational processes3. Here, by whole-genome sequencing analysis of 2,658 cancers as part of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA)4, we reconstruct the life history and evolution of mutational processes and driver mutation sequences of 38 types of cancer. Early oncogenesis is characterized by mutations in a constrained set of driver genes, and specific copy number gains, such as trisomy 7 in glioblastoma and isochromosome 17q in medulloblastoma. The mutational spectrum changes significantly throughout tumour evolution in 40% of samples. A nearly fourfold diversification of driver genes and increased genomic instability are features of later stages. Copy number alterations often occur in mitotic crises, and lead to simultaneous gains of chromosomal segments. Timing analyses suggest that driver mutations often precede diagnosis by many years, if not decades. Together, these results determine the evolutionary trajectories of cancer, and highlight opportunities for early cancer detection.

Mitochondrial Reprogramming Underlies Resistance to BCL-2 Inhibition in Lymphoid Malignancies.

Mitochondrial apoptosis can be effectively targeted in lymphoid malignancies with the FDA-approved B cell lymphoma 2 (BCL-2) inhibitor venetoclax, but resistance to this agent is emerging. We show that venetoclax resistance in chronic lymphocytic leukemia is associated with complex clonal shifts. To identify determinants of resistance, we conducted parallel genome-scale screens of the BCL-2-driven OCI-Ly1 lymphoma cell line after venetoclax exposure along with integrated expression profiling and functional characterization of drug-resistant and engineered cell lines. We identified regulators of lymphoid transcription and cellular energy metabolism as drivers of venetoclax resistance in addition to the known involvement by BCL-2 family members, which were confirmed in patient samples. Our data support the implementation of combinatorial therapy with metabolic modulators to address venetoclax resistance.

RNA sequence analysis reveals macroscopic somatic clonal expansion across normal tissues.

How somatic mutations accumulate in normal cells is poorly understood. A comprehensive analysis of RNA sequencing data from ~6700 samples across 29 normal tissues revealed multiple somatic variants, demonstrating that macroscopic clones can be found in many normal tissues. We found that sun-exposed skin, esophagus, and lung have a higher mutation burden than other tested tissues, which suggests that environmental factors can promote somatic mosaicism. Mutation burden was associated with both age and tissue-specific cell proliferation rate, highlighting that mutations accumulate over both time and number of cell divisions. Finally, normal tissues were found to harbor mutations in known cancer genes and hotspots. This study provides a broad view of macroscopic clonal expansion in human tissues, thus serving as a foundation for associating clonal expansion with environmental factors, aging, and risk of disease.

Growth dynamics in naturally progressing chronic lymphocytic leukaemia.

How the genomic features of a patient's cancer relate to individual disease kinetics remains poorly understood. Here we used the indolent growth dynamics of chronic lymphocytic leukaemia (CLL) to analyse the growth rates and corresponding genomic patterns of leukaemia cells from 107 patients with CLL, spanning decades-long disease courses. We found that CLL commonly demonstrates not only exponential expansion but also logistic growth, which is sigmoidal and reaches a certain steady-state level. Each growth pattern was associated with marked differences in genetic composition, the pace of disease progression and the extent of clonal evolution. In a subset of patients, whose serial samples underwent next-generation sequencing, we found that dynamic changes in the disease course of CLL were shaped by the genetic events that were already present in the early slow-growing stages. Finally, by analysing the growth rates of subclones compared with their parental clones, we quantified the growth advantage conferred by putative CLL drivers in vivo.

Mechanisms of acquired resistance to rapalogs in metastatic renal cell carcinoma.

The mechanistic target of rapamycin (mTOR) is an established therapeutic target in renal cell carcinoma (RCC). Mechanisms of secondary resistance to rapalog therapy in RCC have not been studied previously. We identified six patients with metastatic RCC who initially responded to mTOR inhibitor therapy and then progressed, and had pre-treatment and post-treatment tumor samples available for analysis. We performed deep whole exome sequencing on the paired tumor samples and a blood sample. Sequence data was analyzed using Mutect, CapSeg, Absolute, and Phylogic to identify mutations, copy number changes, and their changes over time. We also performed in vitro functional assays on PBRM1 in RCC cell lines. Five patients had clear cell and one had chromophobe RCC. 434 somatic mutations in 416 genes were identified in the 12 tumor samples. 201 (46%) of mutations were clonal in both samples while 129 (30%) were acquired in the post-treatment samples. Tumor heterogeneity or sampling issues are likely to account for some mutations that were acquired in the post-treatment samples. Three samples had mutations in TSC1; one in PTEN; and none in MTOR. PBRM1 was the only gene in which mutations were acquired in more than one post-treatment sample. We examined the effect of PBRM1 loss in multiple RCC cell lines, and could not identify any effect on rapalog sensitivity in in vitro culture assays. We conclude that mTOR pathway gene mutations did not contribute to rapalog resistance development in these six patients with advanced RCC. Furthermore, mechanisms of resistance to rapalogs in RCC remain unclear and our results suggest that PBRM1 loss may contribute to sensitivity through complex transcriptional effects.

Widespread Chromosomal Losses and Mitochondrial DNA Alterations as Genetic Drivers in Hürthle Cell Carcinoma.

Hürthle cell carcinoma of the thyroid (HCC) is a form of thyroid cancer recalcitrant to radioiodine therapy that exhibits an accumulation of mitochondria. We performed whole-exome sequencing on a cohort of primary, recurrent, and metastatic tumors, and identified recurrent mutations in DAXX, TP53, NRAS, NF1, CDKN1A, ARHGAP35, and the TERT promoter. Parallel analysis of mtDNA revealed recurrent homoplasmic mutations in subunits of complex I of the electron transport chain. Analysis of DNA copy-number alterations uncovered widespread loss of chromosomes culminating in near-haploid chromosomal content in a large fraction of HCC, which was maintained during metastatic spread. This work uncovers a distinct molecular origin of HCC compared with other thyroid malignancies.

Author Correction: Molecular subtypes of diffuse large B cell lymphoma are associated with distinct pathogenic mechanisms and outcomes.

In the version of this article originally published, an asterisk was omitted from Fig. 1a. The asterisk has been added to the figure. Additionally, a "NOTCH2" label was erroneously included in Fig. 4a. The label has been removed. The errors have been corrected in the PDF and HTML versions of this article.

Publisher Correction: Molecular subtypes of diffuse large B cell lymphoma are associated with distinct pathogenic mechanisms and outcomes.

In the version of this article originally published, some text above the "Tri-nucleotide sequence motifs" label in Fig. 2a appeared incorrectly. The text was garbled and should have appeared as nucleotide codes.Additionally, the labels on the bars in Fig. 2c were not italicized in the original publication. These are gene symbols, and they should have been italicized.The colored labels above the graphs in Fig. 4b were also erroneously not italicized. These labels represent gene names and loci, and they should have been italicized.

Molecular subtypes of diffuse large B cell lymphoma are associated with distinct pathogenic mechanisms and outcomes.

Diffuse large B cell lymphoma (DLBCL), the most common lymphoid malignancy in adults, is a clinically and genetically heterogeneous disease that is further classified into transcriptionally defined activated B cell (ABC) and germinal center B cell (GCB) subtypes. We carried out a comprehensive genetic analysis of 304 primary DLBCLs and identified low-frequency alterations, captured recurrent mutations, somatic copy number alterations, and structural variants, and defined coordinate signatures in patients with available outcome data. We integrated these genetic drivers using consensus clustering and identified five robust DLBCL subsets, including a previously unrecognized group of low-risk ABC-DLBCLs of extrafollicular/marginal zone origin; two distinct subsets of GCB-DLBCLs with different outcomes and targetable alterations; and an ABC/GCB-independent group with biallelic inactivation of TP53, CDKN2A loss, and associated genomic instability. The genetic features of the newly characterized subsets, their mutational signatures, and the temporal ordering of identified alterations provide new insights into DLBCL pathogenesis. The coordinate genetic signatures also predict outcome independent of the clinical International Prognostic Index and suggest new combination treatment strategies. More broadly, our results provide a roadmap for an actionable DLBCL classification.

On Rhythms in Neuronal Networks with Recurrent Excitation.

We investigate rhythms in networks of neurons with recurrent excitation, that is, with excitatory cells exciting each other. Recurrent excitation can sustain activity even when the cells in the network are driven below threshold, too weak to fire on their own. This sort of "reverberating" activity is often thought to be the basis of working memory. Recurrent excitation can also lead to "runaway" transitions, sudden transitions to high-frequency firing; this may be related to epileptic seizures. Not all fundamental questions about these phenomena have been answered with clarity in the literature. We focus on three questions here: (1) How much recurrent excitation is needed to sustain reverberating activity? How does the answer depend on parameters? (2) Is there a positive minimum frequency of reverberating activity, a positive "onset frequency"? How does it depend on parameters? (3) When do runaway transitions occur? For reduced models, we give mathematical answers to these questions. We also examine computationally to which extent our findings are reflected in the behavior of biophysically more realistic model networks. Our main results can be summarized as follows. (1) Reverberating activity can be fueled by extremely weak slow recurrent excitation, but only by sufficiently strong fast recurrent excitation. (2) The onset of reverberating activity, as recurrent excitation is strengthened or external drive is raised, occurs at a positive frequency. It is faster when the external drive is weaker (and the recurrent excitation stronger). It is slower when the recurrent excitation has a longer decay time constant. (3) Runaway transitions occur only with fast, not with slow, recurrent excitation. We also demonstrate that the relation between reverberating activity fueled by recurrent excitation and runaway transitions can be visualized in an instructive way by a (generalized) cusp catastrophe surface.