Spatiotemporal analysis of human intestinal development at single-cell resolution.

Development of the human intestine is not well understood. Here, we link single-cell RNA sequencing and spatial transcriptomics to characterize intestinal morphogenesis through time. We identify 101 cell states including epithelial and mesenchymal progenitor populations and programs linked to key morphogenetic milestones. We describe principles of crypt-villus axis formation; neural, vascular, mesenchymal morphogenesis, and immune population of the developing gut. We identify the differentiation hierarchies of developing fibroblast and myofibroblast subtypes and describe diverse functions for these including as vascular niche cells. We pinpoint the origins of Peyer's patches and gut-associated lymphoid tissue (GALT) and describe location-specific immune programs. We use our resource to present an unbiased analysis of morphogen gradients that direct sequential waves of cellular differentiation and define cells and locations linked to rare developmental intestinal disorders. We compile a publicly available online resource, spatio-temporal analysis resource of fetal intestinal development (STAR-FINDer), to facilitate further work.

Apc-mutant cells act as supercompetitors in intestinal tumour initiation.

A delicate equilibrium of WNT agonists and antagonists in the intestinal stem cell (ISC) niche is critical to maintaining the ISC compartment, as it accommodates the rapid renewal of the gut lining. Disruption of this balance by mutations in the tumour suppressor gene APC, which are found in approximately 80% of all human colon cancers, leads to unrestrained activation of the WNT pathway1,2. It has previously been established that Apc-mutant cells have a competitive advantage over wild-type ISCs3. Consequently, Apc-mutant ISCs frequently outcompete all wild-type stem cells within a crypt, thereby reaching clonal fixation in the tissue and initiating cancer formation. However, whether the increased relative fitness of Apc-mutant ISCs involves only cell-intrinsic features or whether Apc mutants are actively involved in the elimination of their wild-type neighbours remains unresolved. Here we show that Apc-mutant ISCs function as bona fide supercompetitors by secreting WNT antagonists, thereby inducing differentiation of neighbouring wild-type ISCs. Lithium chloride prevented the expansion of Apc-mutant clones and the formation of adenomas by rendering wild-type ISCs insensitive to WNT antagonists through downstream activation of WNT by inhibition of GSK3ß. Our work suggests that boosting the fitness of healthy cells to limit the expansion of pre-malignant clones may be a powerful strategy to limit the formation of cancers in high-risk individuals.

A KMT2A-AFF1 gene regulatory network highlights the role of core transcription factors and reveals the regulatory logic of key downstream target genes.

Regulatory interactions mediated by transcription factors (TFs) make up complex networks that control cellular behavior. Fully understanding these gene regulatory networks (GRNs) offers greater insight into the consequences of disease-causing perturbations than can be achieved by studying single TF binding events in isolation. Chromosomal translocations of the lysine methyltransferase 2A (KMT2A) gene produce KMT2A fusion proteins such as KMT2A-AFF1 (previously MLL-AF4), causing poor prognosis acute lymphoblastic leukemias (ALLs) that sometimes relapse as acute myeloid leukemias (AMLs). KMT2A-AFF1 drives leukemogenesis through direct binding and inducing the aberrant overexpression of key genes, such as the anti-apoptotic factor BCL2 and the proto-oncogene MYC However, studying direct binding alone does not incorporate possible network-generated regulatory outputs, including the indirect induction of gene repression. To better understand the KMT2A-AFF1-driven regulatory landscape, we integrated ChIP-seq, patient RNA-seq, and CRISPR essentiality screens to generate a model GRN. This GRN identified several key transcription factors such as RUNX1 that regulate target genes downstream of KMT2A-AFF1 using feed-forward loop (FFL) and cascade motifs. A core set of nodes are present in both ALL and AML, and CRISPR screening revealed several factors that help mediate response to the drug venetoclax. Using our GRN, we then identified a KMT2A-AFF1:RUNX1 cascade that represses CASP9, as well as KMT2A-AFF1-driven FFLs that regulate BCL2 and MYC through combinatorial TF activity. This illustrates how our GRN can be used to better connect KMT2A-AFF1 behavior to downstream pathways that contribute to leukemogenesis, and potentially predict shifts in gene expression that mediate drug response.

Low-cost and clinically applicable copy number profiling using repeat DNA.

BACKGROUND: Somatic copy number alterations (SCNAs) are an important class of genomic alteration in cancer. They are frequently observed in cancer samples, with studies showing that, on average, SCNAs affect 34% of a cancer cell's genome. Furthermore, SCNAs have been shown to be major drivers of tumour development and have been associated with response to therapy and prognosis. Large-scale cancer genome studies suggest that tumours are driven by somatic copy number alterations (SCNAs) or single-nucleotide variants (SNVs). Despite the frequency of SCNAs and their clinical relevance, the use of genomics assays in the clinic is biased towards targeted gene panels, which identify SNVs but provide limited scope to detect SCNAs throughout the genome. There is a need for a comparably low-cost and simple method for high-resolution SCNA profiling. RESULTS: We present conliga, a fully probabilistic method that infers SCNA profiles from a low-cost, simple, and clinically-relevant assay (FAST-SeqS). When applied to 11 high-purity oesophageal adenocarcinoma samples, we obtain good agreement (Spearman's rank correlation coefficient, rs=0.94) between conliga's inferred SCNA profiles using FAST-SeqS data (approximately £14 per sample) and those inferred by ASCAT using high-coverage WGS (gold-standard). We find that conliga outperforms CNVkit (rs=0.89), also applied to FAST-SeqS data, and is comparable to QDNAseq (rs=0.96) applied to low-coverage WGS, which is approximately four-fold more expensive, more laborious and less clinically-relevant. By performing an in silico dilution series experiment, we find that conliga is particularly suited to detecting SCNAs in low tumour purity samples. At two million reads per sample, conliga is able to detect SCNAs in all nine samples at 3% tumour purity and as low as 0.5% purity in one sample. Crucially, we show that conliga's hidden state information can be used to decide when a sample is abnormal or normal, whereas CNVkit and QDNAseq cannot provide this critical information. CONCLUSIONS: We show that conliga provides high-resolution SCNA profiles using a convenient, low-cost assay. We believe conliga makes FAST-SeqS a more clinically valuable assay as well as a useful research tool, enabling inexpensive and fast copy number profiling of pre-malignant and cancer samples.

A Diffusion-like Process Accommodates New Crypts During Clonal Expansion in Human Colonic Epithelium.

BACKGROUND & AIMS: Colorectal cancer (CRC) is thought to arise when the cumulative mutational burden within colonic crypts exceeds a certain threshold that leads to clonal expansion and ultimately neoplastic transformation. Therefore, quantification of the fixation and subsequent expansion of somatic mutations in normal epithelium is key to understanding colorectal cancer initiation. The aim of the present study was to determine how advantaged expansions can be accommodated in the human colon. METHODS: Immunohistochemistry was used to visualize loss of the cancer driver KDM6A in formalin-fixed paraffin-embedded (FFPE) normal human colonic epithelium. Combining microscopy with neural network-based image analysis, we determined the frequencies of KDM6A-mutant crypts and fission/fusion intermediates as well as the spatial distribution of clones. Mathematical modeling then defined the dynamics of their fixation and expansion. RESULTS: Interpretation of the age-related behavior of KDM6A-negative clones revealed significant competitive advantage in intracrypt dynamics as well as a 5-fold increase in crypt fission rate. This was not accompanied by an increase in crypt fusion. Mathematical modeling of crypt spacing identifies evidence for a crypt diffusion process. We define the threshold fission rate at which diffusion fails to accommodate new crypts, which can be exceeded by KRAS activating mutations. CONCLUSIONS: Advantaged gene mutations in KDM6A expand dramatically by crypt fission but not fusion. The crypt diffusion process enables accommodation of the additional crypts up to a threshold value, beyond which polyp growth may occur. The fission rate associated with KRAS mutations offers a potential explanation for KRAS-initiated polyps.

Spatiotemporal regulation of clonogenicity in colorectal cancer xenografts.

Cancer evolution is predominantly studied by focusing on differences in the genetic characteristics of malignant cells within tumors. However, the spatiotemporal dynamics of clonal outgrowth that underlie evolutionary trajectories remain largely unresolved. Here, we sought to unravel the clonal dynamics of colorectal cancer (CRC) expansion in space and time by using a color-based clonal tracing method. This method involves lentiviral red-green-blue (RGB) marking of cell populations, which enabled us to track individual cells and their clonal outgrowth during tumor initiation and growth in a xenograft model. We found that clonal expansion largely depends on the location of a clone, as small clones reside in the center and large clones mostly drive tumor growth at the border. These dynamics are recapitulated in a computational model, which confirms that the clone position within a tumor rather than cell-intrinsic features, is crucial for clonal outgrowth. We also found that no significant clonal loss occurs during tumor growth and clonal dispersal is limited in most models. Our results imply that, in addition to molecular features of clones such as (epi-)genetic differences between cells, clone location and the geometry of tumor growth are crucial for clonal expansion. Our findings suggest that either microenvironmental signals on the tumor border or differences in physical properties within the tumor, are major contributors to explain heterogeneous clonal expansion. Thus, this study provides further insights into the dynamics of solid tumor growth and progression, as well as the origins of tumor cell heterogeneity in a relevant model system.

Oncogenic Gata1 causes stage-specific megakaryocyte differentiation delay.

The megakaryocyte/erythroid Transient Myeloproliferative Disorder (TMD) in newborns with Down Syndrome (DS) occurs when N-terminal truncating mutations of the hemopoietic transcription factor GATA1, that produce GATA1short protein (GATA1s), are acquired early in development. Prior work has shown that murine GATA1s, by itself, causes a transient yolk sac myeloproliferative disorder. However, it is unclear where in the hemopoietic cellular hierarchy GATA1s exerts its effects to produce this myeloproliferative state. Here, through a detailed examination of hemopoiesis from murine GATA1s ES cells and GATA1s embryos we define defects in erythroid and megakaryocytic differentiation that occur relatively late in hemopoiesis. GATA1s causes an arrest late in erythroid differentiation in vivo, and even more profoundly in ES-cell derived cultures, with a marked reduction of Ter-119 cells and reduced erythroid gene expression. In megakaryopoiesis, GATA1s causes a differentiation delay at a specific stage, with accumulation of immature, kit-expressing CD41hi megakaryocytic cells. In this specific megakaryocytic compartment, there are increased numbers of GATA1s cells in S-phase of cell cycle and reduced number of apoptotic cells compared to GATA1 cells in the same cell compartment. There is also a delay in maturation of these immature GATA1s megakaryocytic lineage cells compared to GATA1 cells at the same stage of differentiation. Finally, even when GATA1s megakaryocytic cells mature, they mature aberrantly with altered megakaryocyte-specific gene expression and activity of the mature megakaryocyte enzyme, acetylcholinesterase. These studies pinpoint the hemopoietic compartment where GATA1s megakaryocyte myeloproliferation occurs, defining where molecular studies should now be focussed to understand the oncogenic action of GATA1s.

Single-cell fate decisions of bipotential hematopoietic progenitors.

PURPOSE OF REVIEW: In hematopoiesis, rapid cell fate decisions are necessary for timely responses to environmental stimuli resulting in the production of diverse types of blood cells. Early studies have led to a hierarchical, tree-like view of hematopoiesis with hematopoietic stem cells residing at the apex and serially branching out to give rise to bipotential progenitors with increasingly restricted lineage potential. Recent single-cell studies have challenged some aspects of the classical model of hematopoiesis. Here, we review the latest articles on cell fate decision in hematopoietic progenitors, highlighting single-cell studies that have questioned previously established concepts and those that have reaffirmed them. RECENT FINDINGS: The hierarchical organization of hematopoiesis and the importance of transcription factors have been largely validated at the single-cell level. In contrast, single-cell studies have shown that lineage commitment is progressive rather than switch-like as originally proposed. Furthermore, the reconstruction of cell fate paths suggested the existence of a gradient of hematopoietic progenitors that are in a continuum of changing fate probabilities rather than in a static bipotential state, leading us to reconsider the notion of bipotential progenitors. SUMMARY: Single-cell transcriptomic and proteomic studies have transformed our view of lineage commitment and offer a drastically different perspective on hematopoiesis.

Metabolic switches and adaptations deduced from the proteomes of Streptomyces coelicolor wild type and phoP mutant grown in batch culture.

Bacteria in the genus Streptomyces are soil-dwelling oligotrophs and important producers of secondary metabolites. Previously, we showed that global messenger RNA expression was subject to a series of metabolic and regulatory switches during the lifetime of a fermentor batch culture of Streptomyces coelicolor M145. Here we analyze the proteome from eight time points from the same fermentor culture and, because phosphate availability is an important regulator of secondary metabolite production, compare this to the proteome of a similar time course from an S. coelicolor mutant, INB201 (ΔphoP), defective in the control of phosphate utilization. The proteomes provide a detailed view of enzymes involved in central carbon and nitrogen metabolism. Trends in protein expression over the time courses were deduced from a protein abundance index, which also revealed the importance of stress pathway proteins in both cultures. As expected, the ΔphoP mutant was deficient in expression of PhoP-dependent genes, and several putatively compensatory metabolic and regulatory pathways for phosphate scavenging were detected. Notably there is a succession of switches that coordinately induce the production of enzymes for five different secondary metabolite biosynthesis pathways over the course of the batch cultures.

Inferring the time-invariant topology of a nonlinear sparse gene regulatory network using fully Bayesian spline autoregression.

We propose a semiparametric Bayesian model, based on penalized splines, for the recovery of the time-invariant topology of a causal interaction network from longitudinal data. Our motivation is inference of gene regulatory networks from low-resolution microarray time series, where existence of nonlinear interactions is well known. Parenthood relations are mapped by augmenting the model with kinship indicators and providing these with either an overall or gene-wise hierarchical structure. Appropriate specification of the prior is crucial to control the flexibility of the splines, especially under circumstances of scarce data; thus, we provide an informative, proper prior. Substantive improvement in network inference over a linear model is demonstrated using synthetic data drawn from ordinary differential equation models and gene expression from an experimental data set of the Arabidopsis thaliana circadian rhythm.

Pomelo II: finding differentially expressed genes.

Pomelo II (http://pomelo2.bioinfo.cnio.es) is an open-source, web-based, freely available tool for the analysis of gene (and protein) expression and tissue array data. Pomelo II implements: permutation-based tests for class comparisons (t-test, ANOVA) and regression; survival analysis using Cox model; contingency table analysis with Fisher's exact test; linear models (of which t-test and ANOVA are especial cases) that allow additional covariates for complex experimental designs and use empirical Bayes moderated statistics. Permutation-based and Cox model analysis use parallel computing, which permits taking advantage of multicore CPUs and computing clusters. Access to, and further analysis of, additional biological information and annotations (PubMed references, Gene Ontology terms, KEGG and Reactome pathways) are available either for individual genes (from clickable links in tables and figures) or sets of genes. The source code is available, allowing for extending and reusing the software. A comprehensive test suite is also available, and covers both the user interface and the numerical results. The possibility of including additional covariates, parallelization of computation, open-source availability of the code and comprehensive testing suite make Pomelo II a unique tool.

Heterogeneity in clone dynamics within and adjacent to intestinal tumours identified by Dre-mediated lineage tracing.

Somatic models of tissue pathology commonly use induction of gene-specific mutations in mice mediated by spatiotemporal regulation of Cre recombinase. Subsequent investigation of the onset and development of disease can be limited by the inability to track changing cellular behaviours over time. Here, a lineage-tracing approach based on ligand-dependent activation of Dre recombinase that can be employed independently of Cre is described. The clonal biology of the intestinal epithelium following Cre-mediated stabilisation of ß-catenin reveals that, within tumours, many new clones rapidly become extinct. Surviving clones show accelerated population of tumour glands compared to normal intestinal crypts but in a non-uniform manner, indicating that intra-tumour glands follow heterogeneous dynamics. In tumour-adjacent epithelia, clone sizes are smaller than in the background epithelia, as a whole. This suggests a zone of ∼seven crypt diameters within which clone expansion is inhibited by tumours and that may facilitate their growth.

Continuous clonal labeling reveals uniform progenitor potential in the adult exocrine pancreas.

The tissue dynamics that govern maintenance and regeneration of the pancreas remain largely unknown. In particular, the presence and nature of a cellular hierarchy remains a topic of debate. Previous lineage tracing strategies in the pancreas relied on specific marker genes for clonal labeling, which left other populations untested and failed to account for potential widespread phenotypical plasticity. Here we employed a tracing system that depends on replication-induced clonal marks. We found that, in homeostasis, steady acinar replacement events characterize tissue dynamics, to which all acinar cells have an equal ability to contribute. Similarly, regeneration following pancreatitis was best characterized by an acinar self-replication model because no evidence of a cellular hierarchy was detected. In particular, rapid regeneration in the pancreas was found to be driven by an accelerated rate of acinar fission-like events. These results provide a comprehensive and quantitative model of cell dynamics in the exocrine pancreas.

The onset of circulation triggers a metabolic switch required for endothelial to hematopoietic transition.

Hematopoietic stem cells (HSCs) emerge during development from the vascular wall of the main embryonic arteries. The onset of circulation triggers several processes that provide critical external factors for HSC generation. Nevertheless, it is not fully understood how and when the onset of circulation affects HSC emergence. Here we show that in Ncx1-/- mouse embryos devoid of circulation the HSC lineage develops until the phenotypic pro-HSC stage. However, these cells reside in an abnormal microenvironment, fail to activate the hematopoietic program downstream of Runx1, and are functionally impaired. Single-cell transcriptomics shows that during the endothelial-to-hematopoietic transition, Ncx1-/- cells fail to undergo a glycolysis to oxidative phosphorylation metabolic switch present in wild-type cells. Interestingly, experimental activation of glycolysis results in decreased intraembryonic hematopoiesis. Our results suggest that the onset of circulation triggers metabolic changes that allow HSC generation to proceed.

The dynamic architecture of the metabolic switch in Streptomyces coelicolor.

BACKGROUND: During the lifetime of a fermenter culture, the soil bacterium S. coelicolor undergoes a major metabolic switch from exponential growth to antibiotic production. We have studied gene expression patterns during this switch, using a specifically designed Affymetrix genechip and a high-resolution time-series of fermenter-grown samples. RESULTS: Surprisingly, we find that the metabolic switch actually consists of multiple finely orchestrated switching events. Strongly coherent clusters of genes show drastic changes in gene expression already many hours before the classically defined transition phase where the switch from primary to secondary metabolism was expected. The main switch in gene expression takes only 2 hours, and changes in antibiotic biosynthesis genes are delayed relative to the metabolic rearrangements. Furthermore, global variation in morphogenesis genes indicates an involvement of cell differentiation pathways in the decision phase leading up to the commitment to antibiotic biosynthesis. CONCLUSIONS: Our study provides the first detailed insights into the complex sequence of early regulatory events during and preceding the major metabolic switch in S. coelicolor, which will form the starting point for future attempts at engineering antibiotic production in a biotechnological setting.

Asterias: integrated analysis of expression and aCGH data using an open-source, web-based, parallelized software suite.

Asterias (http://www.asterias.info) is an open-source, web-based, suite for the analysis of gene expression and aCGH data. Asterias implements validated statistical methods, and most of the applications use parallel computing, which permits taking advantage of multicore CPUs and computing clusters. Access to, and further analysis of, additional biological information and annotations (PubMed references, Gene Ontology terms, KEGG and Reactome pathways) are available either for individual genes (from clickable links in tables and figures) or sets of genes. These applications cover from array normalization to imputation and preprocessing, differential gene expression analysis, class and survival prediction and aCGH analysis. The source code is available, allowing for extention and reuse of the software. The links and analysis of additional functional information, parallelization of computation and open-source availability of the code make Asterias a unique suite that can exploit features specific to web-based environments.

Single-Cell Proteomics Reveal that Quantitative Changes in Co-expressed Lineage-Specific Transcription Factors Determine Cell Fate.

Hematopoiesis provides an accessible system for studying the principles underlying cell-fate decisions in stem cells. Proposed models of hematopoiesis suggest that quantitative changes in lineage-specific transcription factors (LS-TFs) underlie cell-fate decisions. However, evidence for such models is lacking as TF levels are typically measured via RNA expression rather than by analyzing temporal changes in protein abundance. Here, we used single-cell mass cytometry and absolute quantification by mass spectrometry to capture the temporal dynamics of TF protein expression in individual cells during human erythropoiesis. We found that LS-TFs from alternate lineages are co-expressed, as proteins, in individual early progenitor cells and quantitative changes of LS-TFs occur gradually rather than abruptly to direct cell-fate decisions. Importantly, upregulation of a megakaryocytic TF in early progenitors is sufficient to deviate cells from an erythroid to a megakaryocyte trajectory, showing that quantitative changes in protein abundance of LS-TFs in progenitors can determine alternate cell fates.

Quantifying Microsatellite Mutation Rates from Intestinal Stem Cell Dynamics in Msh2-Deficient Murine Epithelium.

Microsatellite sequences have an enhanced susceptibility to mutation, and can act as sentinels indicating elevated mutation rates and increased risk of cancer. The probability of mutant fixation within the intestinal epithelium is dictated by a combination of stem cell dynamics and mutation rate. Here, we exploit this relationship to infer microsatellite mutation rates. First a sensitive, multiplexed, and quantitative method for detecting somatic changes in microsatellite length was developed that allowed the parallel detection of mutant [CA]n sequences from hundreds of low-input tissue samples at up to 14 loci. The method was applied to colonic crypts in Mus musculus, and enabled detection of mutant subclones down to 20% of the cellularity of the crypt (∼50 of 250 cells). By quantifying age-related increases in clone frequencies for multiple loci, microsatellite mutation rates in wild-type and Msh2-deficient epithelium were established. An average 388-fold increase in mutation per mitosis rate was observed in Msh2-deficient epithelium (2.4 × 10-2) compared to wild-type epithelium (6.2 × 10-5).

Phospho-regulation of ATOH1 Is Required for Plasticity of Secretory Progenitors and Tissue Regeneration.

The intestinal epithelium is largely maintained by self-renewing stem cells but with apparently committed progenitors also contributing, particularly following tissue damage. However, the mechanism of, and requirement for, progenitor plasticity in mediating pathological response remain unknown. Here we show that phosphorylation of the transcription factor Atoh1 is required for both the contribution of secretory progenitors to the stem cell pool and for a robust regenerative response. As confirmed by lineage tracing, Atoh1+ cells (Atoh1(WT)CreERT2 mice) give rise to multilineage intestinal clones both in the steady state and after tissue damage. In a phosphomutant Atoh1(9S/T-A)CreERT2 line, preventing phosphorylation of ATOH1 protein acts to promote secretory differentiation and inhibit the contribution of progenitors to self-renewal. Following chemical colitis, Atoh1+ cells of Atoh1(9S/T-A)CreERT2 mice have reduced clonogenicity that affects overall regeneration. Progenitor plasticity maintains robust self-renewal in the intestinal epithelium, and the balance between stem and progenitor fate is directly coordinated by ATOH1 multisite phosphorylation.