Self-patterning of human stem cells into post-implantation lineages.

Investigating human development is a substantial scientific challenge due to the technical and ethical limitations of working with embryonic samples. In the face of these difficulties, stem cells have provided an alternative to experimentally model inaccessible stages of human development in vitro1-13. Here we show that human pluripotent stem cells can be triggered to self-organize into three-dimensional structures that recapitulate some key spatiotemporal events of early human post-implantation embryonic development. Our system reproducibly captures spontaneous differentiation and co-development of embryonic epiblast-like and extra-embryonic hypoblast-like lineages, establishes key signalling hubs with secreted modulators and undergoes symmetry breaking-like events. Single-cell transcriptomics confirms differentiation into diverse cell states of the perigastrulating human embryo14,15 without establishing placental cell types, including signatures of post-implantation epiblast, amniotic ectoderm, primitive streak, mesoderm, early extra-embryonic endoderm, as well as initial yolk sac induction. Collectively, our system captures key features of human embryonic development spanning from Carnegie stage16 4-7, offering a reproducible, tractable and scalable experimental platform to understand the basic cellular and molecular mechanisms that underlie human development, including new opportunities to dissect congenital pathologies with high throughput.

T-bet suppresses proliferation of malignant B cells in chronic lymphocytic leukemia.

The T-box transcription factor T-bet is known as a master regulator of T-cell response but its role in malignant B cells is not sufficiently explored. Here, we conducted single-cell resolved multi-omics analyses of malignant B cells from patients with chronic lymphocytic leukemia (CLL) and studied a CLL mouse model with genetic knockout of TBX21. We found that T-bet acts as a tumor suppressor in malignant B cells by decreasing their proliferation rate. NF-κB activity induced by inflammatory signals provided by the microenvironment, triggered T-bet expression which impacted on promoter proximal and distal chromatin co-accessibility and controlled a specific gene signature by mainly suppressing transcription. Gene set enrichment analysis identified a positive regulation of interferon signaling, and a negative control of proliferation by T-bet. In line, we showed that T-bet represses cell cycling and is associated with longer overall survival of CLL patients. Our study uncovers a novel tumor suppressive role of T-bet in malignant B cells via its regulation of inflammatory processes and cell cycling which has implications for stratification and therapy of CLL patients. Linking T-bet activity to inflammation explains the good prognostic role of genetic alterations in inflammatory signaling pathways in CLL.

Epigenetic regulator function through mouse gastrulation.

During ontogeny, proliferating cells become restricted in their fate through the combined action of cell-type-specific transcription factors and ubiquitous epigenetic machinery, which recognizes universally available histone residues or nucleotides in a context-dependent manner1,2. The molecular functions of these regulators are generally well understood, but assigning direct developmental roles to them is hampered by complex mutant phenotypes that often emerge after gastrulation3,4. Single-cell RNA sequencing and analytical approaches have explored this highly conserved, dynamic period across numerous model organisms5-8, including mouse9-18. Here we advance these strategies using a combined zygotic perturbation and single-cell RNA-sequencing platform in which many mutant mouse embryos can be assayed simultaneously, recovering robust  morphological and transcriptional information across a panel of ten essential regulators. Deeper analysis of central Polycomb repressive complex (PRC) 1 and 2 components indicates substantial cooperativity, but distinguishes a dominant role for PRC2 in restricting the germline. Moreover, PRC mutant phenotypes emerge after gross epigenetic and transcriptional changes within the initial conceptus prior to gastrulation. Our experimental framework may eventually lead to a fully quantitative view of how cellular diversity emerges using an identical genetic template and from a single totipotent cell.

Lambda3: homology search for protein, nucleotide, and bisulfite-converted sequences.

MOTIVATION: Local alignments of query sequences in large databases represent a core part of metagenomic studies and facilitate homology search. Following the development of NCBI Blast, many applications aimed to provide faster and equally sensitive local alignment frameworks. Most applications focus on protein alignments, while only few also facilitate DNA-based searches. None of the established programs allow searching DNA sequences from bisulfite sequencing experiments commonly used for DNA methylation profiling, for which specific alignment strategies need to be implemented. RESULTS: Here, we introduce Lambda3, a new version of the local alignment application Lambda. Lambda3 is the first solution that enables the search of protein, nucleotide as well as bisulfite-converted nucleotide query sequences. Its protein mode achieves comparable performance to that of the highly optimized protein alignment application Diamond, while the nucleotide mode consistently outperforms established local nucleotide aligners. Combined, Lambda3 presents a universal local alignment framework that enables fast and sensitive homology searches for a wide range of use-cases. AVAILABILITY AND IMPLEMENTATION: Lambda3 is free and open-source software publicly available at https://github.com/seqan/lambda/.

Molecular recording of mammalian embryogenesis.

Ontogeny describes the emergence of complex multicellular organisms from single totipotent cells. This field is particularly challenging in mammals, owing to the indeterminate relationship between self-renewal and differentiation, variation in progenitor field sizes, and internal gestation in these animals. Here we present a flexible, high-information, multi-channel molecular recorder with a single-cell readout and apply it as an evolving lineage tracer to assemble mouse cell-fate maps from fertilization through gastrulation. By combining lineage information with single-cell RNA sequencing profiles, we recapitulate canonical developmental relationships between different tissue types and reveal the nearly complete transcriptional convergence of endodermal cells of extra-embryonic and embryonic origins. Finally, we apply our cell-fate maps to estimate the number of embryonic progenitor cells and their degree of asymmetric partitioning during specification. Our approach enables massively parallel, high-resolution recording of lineage and other information in mammalian systems, which will facilitate the construction of a quantitative framework for understanding developmental processes.

Focal structural variants revealed by whole genome sequencing disrupt the histone demethylase KDM4C in B-cell lymphomas.

Histone methylation-modifiers, such as EZH2 and KMT2D, are recurrently altered in B-cell lymphomas. To comprehensively describe the landscape of alterations affecting genes encoding histone methylation-modifiers in lymphomagenesis we investigated whole genome and transcriptome data of 186 mature B-cell lymphomas sequenced in the ICGC MMML-Seq project. Besides confirming common alterations of KMT2D (47% of cases), EZH2 (17%), SETD1B (5%), PRDM9 (4%), KMT2C (4%), and SETD2 (4%), also identified by prior exome or RNA-sequencing studies, we here found recurrent alterations to KDM4C in chromosome 9p24, encoding a histone demethylase. Focal structural variation was the main mechanism of KDM4C alterations, and was independent from 9p24 amplification. We also identified KDM4C alterations in lymphoma cell lines including a focal homozygous deletion in a classical Hodgkin lymphoma cell line. By integrating RNA-sequencing and genome sequencing data we predict that KDM4C structural variants result in loss-offunction. By functional reconstitution studies in cell lines, we provide evidence that KDM4C can act as a tumor suppressor. Thus, we show that identification of structural variants in whole genome sequencing data adds to the comprehensive description of the mutational landscape of lymphomas and, moreover, establish KDM4C as a putative tumor suppressive gene recurrently altered in subsets of B-cell derived lymphomas.

T-cell prolymphocytic leukemia is associated with deregulation of oncogenic microRNAs on transcriptional and epigenetic level.

Deregulation of micro(mi)-RNAs is a common mechanism in tumorigenesis. We investigated the expression of 2083 miRNAs in T-cell prolymphocytic leukemia (T-PLL). Compared to physiologic CD4+ and CD8+ T-cell subsets, 111 miRNAs were differentially expressed in T-PLL. Of these, 33 belonged to miRNA gene clusters linked to cancer. Genomic variants affecting miRNAs were infrequent with the notable exception of copy number aberrations. Remarkably, we found strong upregulation of the miR-200c/-141 cluster in T-PLL to be associated with DNA hypomethylation and active promoter marks. Our findings suggest that copy number aberrations and epigenetic changes could contribute to miRNA deregulation in T-PLL.

RLM: fast and simplified extraction of read-level methylation metrics from bisulfite sequencing data.

SUMMARY: Bisulfite sequencing data provide value beyond the straightforward methylation assessment by analyzing single-read patterns. Over the past years, various metrics have been established to explore this layer of information. However, limited compatibility with alignment tools, reference genomes or the measurements they provide present a bottleneck for most groups to routinely perform read-level analysis. To address this, we developed RLM, a fast and scalable tool for the computation of several frequently used read-level methylation statistics. RLM supports standard alignment tools, works independently of the reference genome and handles most sequencing experiment designs. RLM can process large input files with a billion reads in just a few hours on common workstations. AVAILABILITY AND IMPLEMENTATION: https://github.com/sarahet/RLM. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Nanopype: a modular and scalable nanopore data processing pipeline.

SUMMARY: Long-read third-generation nanopore sequencing enables researchers to now address a range of questions that are difficult to tackle with short read approaches. The rapidly expanding user base and continuously increasing throughput have sparked the development of a growing number of specialized analysis tools. However, streamlined processing of nanopore datasets using reproducible and transparent workflows is still lacking. Here we present Nanopype, a nanopore data processing pipeline that integrates a diverse set of established bioinformatics software while maintaining consistent and standardized output formats. Seamless integration into compute cluster environments makes the framework suitable for high-throughput applications. As a result, Nanopype facilitates comparability of nanopore data analysis workflows and thereby should enhance the reproducibility of biological insights. AVAILABILITY AND IMPLEMENTATION: https://github.com/giesselmann/nanopype, https://nanopype.readthedocs.io. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

metilene: fast and sensitive calling of differentially methylated regions from bisulfite sequencing data.

The detection of differentially methylated regions (DMRs) is a necessary prerequisite for characterizing different epigenetic states. We present a novel program, metilene, to identify DMRs within whole-genome and targeted data with unrivaled specificity and sensitivity. A binary segmentation algorithm combined with a two-dimensional statistical test allows the detection of DMRs in large methylation experiments with multiple groups of samples in minutes rather than days using off-the-shelf hardware. metilene outperforms other state-of-the-art tools for low coverage data and can estimate missing data. Hence, metilene is a versatile tool to study the effect of epigenetic modifications in differentiation/development, tumorigenesis, and systems biology on a global, genome-wide level. Whether in the framework of international consortia with dozens of samples per group, or even without biological replicates, it produces highly significant and reliable results.

MINCR is a MYC-induced lncRNA able to modulate MYC's transcriptional network in Burkitt lymphoma cells.

Despite the established role of the transcription factor MYC in cancer, little is known about the impact of a new class of transcriptional regulators, the long noncoding RNAs (lncRNAs), on MYC ability to influence the cellular transcriptome. Here, we have intersected RNA-sequencing data from two MYC-inducible cell lines and a cohort of 91 B-cell lymphomas with or without genetic variants resulting in MYC overexpression. We identified 13 lncRNAs differentially expressed in IG-MYC-positive Burkitt lymphoma and regulated in the same direction by MYC in the model cell lines. Among them, we focused on a lncRNA that we named MYC-induced long noncoding RNA (MINCR), showing a strong correlation with MYC expression in MYC-positive lymphomas. To understand its cellular role, we performed RNAi and found that MINCR knockdown is associated with an impairment in cell cycle progression. Differential gene expression analysis after RNAi showed a significant enrichment of cell cycle genes among the genes down-regulated after MINCR knockdown. Interestingly, these genes are enriched in MYC binding sites in their promoters, suggesting that MINCR acts as a modulator of the MYC transcriptional program. Accordingly, MINCR knockdown was associated with a reduction in MYC binding to the promoters of selected cell cycle genes. Finally, we show that down-regulation of Aurora kinases A and B and chromatin licensing and DNA replication factor 1 may explain the reduction in cellular proliferation observed on MINCR knockdown. We, therefore, suggest that MINCR is a newly identified player in the MYC transcriptional network able to control the expression of cell cycle genes.

Alterations of microRNA and microRNA-regulated messenger RNA expression in germinal center B-cell lymphomas determined by integrative sequencing analysis.

MicroRNA are well-established players in post-transcriptional gene regulation. However, information on the effects of microRNA deregulation mainly relies on bioinformatic prediction of potential targets, whereas proof of the direct physical microRNA/target messenger RNA interaction is mostly lacking. Within the International Cancer Genome Consortium Project "Determining Molecular Mechanisms in Malignant Lymphoma by Sequencing", we performed miRnome sequencing from 16 Burkitt lymphomas, 19 diffuse large B-cell lymphomas, and 21 follicular lymphomas. Twenty-two miRNA separated Burkitt lymphomas from diffuse large B-cell lymphomas/follicular lymphomas, of which 13 have shown regulation by MYC. Moreover, we found expression of three hitherto unreported microRNA. Additionally, we detected recurrent mutations of hsa-miR-142 in diffuse large B-cell lymphomas and follicular lymphomas, and editing of the hsa-miR-376 cluster, providing evidence for microRNA editing in lymphomagenesis. To interrogate the direct physical interactions of microRNA with messenger RNA, we performed Argonaute-2 photoactivatable ribonucleoside-enhanced cross-linking and immunoprecipitation experiments. MicroRNA directly targeted 208 messsenger RNA in the Burkitt lymphomas and 328 messenger RNA in the non-Burkitt lymphoma models. This integrative analysis discovered several regulatory pathways of relevance in lymphomagenesis including Ras, PI3K-Akt and MAPK signaling pathways, also recurrently deregulated in lymphomas by mutations. Our dataset reveals that messenger RNA deregulation through microRNA is a highly relevant mechanism in lymphomagenesis.

Reconstructing axial progenitor field dynamics in mouse stem cell-derived embryoids.

Embryogenesis requires substantial coordination to translate genetic programs to the collective behavior of differentiating cells, but understanding how cellular decisions control tissue morphology remains conceptually and technically challenging. Here, we combine continuous Cas9-based molecular recording with a mouse embryonic stem cell-based model of the embryonic trunk to build single-cell phylogenies that describe the behavior of transient, multipotent neuro-mesodermal progenitors (NMPs) as they commit into neural and somitic cell types. We find that NMPs show subtle transcriptional signatures related to their recent differentiation and contribute to downstream lineages through a surprisingly broad distribution of individual fate outcomes. Although decision-making can be heavily influenced by environmental cues to induce morphological phenotypes, axial progenitors intrinsically mature over developmental time to favor the neural lineage. Using these data, we present an experimental and analytical framework for exploring the non-homeostatic dynamics of transient progenitor populations as they shape complex tissues during critical developmental windows.

TET activity safeguards pluripotency throughout embryonic dormancy.

Dormancy is an essential biological process for the propagation of many life forms through generations and stressful conditions. Early embryos of many mammals are preservable for weeks to months within the uterus in a dormant state called diapause, which can be induced in vitro through mTOR inhibition. Cellular strategies that safeguard original cell identity within the silent genomic landscape of dormancy are not known. Here we show that the protection of cis-regulatory elements from silencing is key to maintaining pluripotency in the dormant state. We reveal a TET-transcription factor axis, in which TET-mediated DNA demethylation and recruitment of methylation-sensitive transcription factor TFE3 drive transcriptionally inert chromatin adaptations during dormancy transition. Perturbation of TET activity compromises pluripotency and survival of mouse embryos under dormancy, whereas its enhancement improves survival rates. Our results reveal an essential mechanism for propagating the cellular identity of dormant cells, with implications for regeneration and disease.

Basal delamination during mouse gastrulation primes pluripotent cells for differentiation.

The blueprint of the mammalian body plan is laid out during gastrulation, when a trilaminar embryo is formed. This process entails a burst of proliferation, the ingression of embryonic epiblast cells at the primitive streak, and their priming toward primitive streak fates. How these different events are coordinated remains unknown. Here, we developed and characterized a 3D culture of self-renewing mouse embryonic cells that captures the main transcriptional and architectural features of the early gastrulating mouse epiblast. Using this system in combination with microfabrication and in vivo experiments, we found that proliferation-induced crowding triggers delamination of cells that express high levels of the apical polarity protein aPKC. Upon delamination, cells become more sensitive to Wnt signaling and upregulate the expression of primitive streak markers such as Brachyury. This mechanistic coupling between ingression and differentiation ensures that the right cell types become specified at the right place during embryonic development.

Targeting oncogenic TERT promoter variants by allele-specific epigenome editing.

BACKGROUND: Activation of dominant oncogenes by small or structural genomic alterations is a common driver mechanism in many cancers. Silencing of such dominantly activated oncogenic alleles, thus, is a promising strategy to treat cancer. Recently, allele-specific epigenome editing (ASEE) has been described as a means to reduce transcription of genes in an allele-specific manner. In cancer, specificity to an oncogenic allele can be reached by either targeting directly a pathogenic single-nucleotide variant or a polymorphic single-nucleotide variant linked to the oncogenic allele. To investigate the potential of ASEE in cancer, we here explored this approach by targeting variants at the TERT promoter region. The TERT promoter region has been described as one of the most frequently mutated non-coding cancer drivers. RESULTS: Sequencing of the TERT promoter in cancer cell lines showed 53% (41/77) to contain at least one heterozygous sequence variant allowing allele distinction. We chose the hepatoblastoma cell line Hep-G2 and the lung cancer cell line A-549 for this proof-of-principle study, as they contained two different kinds of variants, namely the activating mutation C228T in the TERT core promoter and the common SNP rs2853669 in the THOR region, respectively. These variants were targeted in an allele-specific manner using sgRNA-guided dCas9-DNMT3A-3L complexes. In both cell lines, we successfully introduced DNA methylation specifically to the on-target allele of the TERT promoter with limited background methylation on the off-target allele or an off-target locus (VEGFA), respectively. We observed a maximum CpG methylation gain of 39% and 76% on the target allele when targeting the activating mutation and the common SNP, respectively. The epigenome editing translated into reduced TERT RNA expression in Hep-G2. CONCLUSIONS: We applied an ASEE-mediated approach to silence TERT allele specifically. Our results show that the concept of dominant oncogene inactivation by allele-specific epigenome editing can be successfully translated into cancer models. This new strategy may have important advantages in comparison with existing therapeutic approaches, e.g., targeting telomerase, especially with regard to reducing adverse side effects.

Spatiotemporal transcriptomic maps of whole mouse embryos at the onset of organogenesis.

Spatiotemporal orchestration of gene expression is required for proper embryonic development. The use of single-cell technologies has begun to provide improved resolution of early regulatory dynamics, including detailed molecular definitions of most cell states during mouse embryogenesis. Here we used Slide-seq to build spatial transcriptomic maps of complete embryonic day (E) 8.5 and E9.0, and partial E9.5 embryos. To support their utility, we developed sc3D, a tool for reconstructing and exploring three-dimensional 'virtual embryos', which enables the quantitative investigation of regionalized gene expression patterns. Our measurements along the main embryonic axes of the developing neural tube revealed several previously unannotated genes with distinct spatial patterns. We also characterized the conflicting transcriptional identity of 'ectopic' neural tubes that emerge in Tbx6 mutant embryos. Taken together, we present an experimental and computational framework for the spatiotemporal investigation of whole embryonic structures and mutant phenotypes.

Dynamic antagonism between key repressive pathways maintains the placental epigenome.

DNA and Histone 3 Lysine 27 methylation typically function as repressive modifications and operate within distinct genomic compartments. In mammals, the majority of the genome is kept in a DNA methylated state, whereas the Polycomb repressive complexes regulate the unmethylated CpG-rich promoters of developmental genes. In contrast to this general framework, the extra-embryonic lineages display non-canonical, globally intermediate DNA methylation levels, including disruption of local Polycomb domains. Here, to better understand this unusual landscape's molecular properties, we genetically and chemically perturbed major epigenetic pathways in mouse trophoblast stem cells. We find that the extra-embryonic epigenome reflects ongoing and dynamic de novo methyltransferase recruitment, which is continuously antagonized by Polycomb to maintain intermediate, locally disordered methylation. Despite its disorganized molecular appearance, our data point to a highly controlled equilibrium between counteracting repressors within extra-embryonic cells, one that can seemingly persist indefinitely without bistable features typically seen for embryonic forms of epigenetic regulation.

T-REX17 is a transiently expressed non-coding RNA essential for human endoderm formation.

Long non-coding RNAs (lncRNAs) have emerged as fundamental regulators in various biological processes, including embryonic development and cellular differentiation. Despite much progress over the past decade, the genome-wide annotation of lncRNAs remains incomplete and many known non-coding loci are still poorly characterized. Here, we report the discovery of a previously unannotated lncRNA that is transcribed 230 kb upstream of the SOX17 gene and located within the same topologically associating domain. We termed it T-REX17 (Transcript Regulating Endoderm and activated by soX17) and show that it is induced following SOX17 activation but its expression is more tightly restricted to early definitive endoderm. Loss of T-REX17 affects crucial functions independent of SOX17 and leads to an aberrant endodermal transcriptome, signaling pathway deregulation and epithelial to mesenchymal transition defects. Consequently, cells lacking the lncRNA cannot further differentiate into more mature endodermal cell types. Taken together, our study identified and characterized T-REX17 as a transiently expressed and essential non-coding regulator in early human endoderm differentiation.