Cell Types of the Human Retina and Its Organoids at Single-Cell Resolution.

Human organoids recapitulating the cell-type diversity and function of their target organ are valuable for basic and translational research. We developed light-sensitive human retinal organoids with multiple nuclear and synaptic layers and functional synapses. We sequenced the RNA of 285,441 single cells from these organoids at seven developmental time points and from the periphery, fovea, pigment epithelium and choroid of light-responsive adult human retinas, and performed histochemistry. Cell types in organoids matured in vitro to a stable "developed" state at a rate similar to human retina development in vivo. Transcriptomes of organoid cell types converged toward the transcriptomes of adult peripheral retinal cell types. Expression of disease-associated genes was cell-type-specific in adult retina, and cell-type specificity was retained in organoids. We implicate unexpected cell types in diseases such as macular degeneration. This resource identifies cellular targets for studying disease mechanisms in organoids and for targeted repair in human retinas.

Convergence of multiple RNA-silencing pathways on GW182/TNRC6.

The RNA-binding protein TRIM71/LIN-41 is a phylogenetically conserved developmental regulator that functions in mammalian stem cell reprogramming, brain development, and cancer. TRIM71 recognizes target mRNAs through hairpin motifs and silences them through molecular mechanisms that await identification. Here, we uncover that TRIM71 represses its targets through RNA-supported interaction with TNRC6/GW182, a core component of the miRNA-induced silencing complex (miRISC). We demonstrate that AGO2, TRIM71, and UPF1 each recruit TNRC6 to specific sets of transcripts to silence them. As cellular TNRC6 levels are limiting, competition occurs among the silencing pathways, such that the loss of AGO proteins or of AGO binding to TNRC6 enhances the activities of the other pathways. We conclude that a miRNA-like silencing activity is shared among different mRNA silencing pathways and that the use of TNRC6 as a central hub provides a means to integrate their activities.

Phenotypic landscape of intestinal organoid regeneration.

The development of intestinal organoids from single adult intestinal stem cells in vitro recapitulates the regenerative capacity of the intestinal epithelium1,2. Here we unravel the mechanisms that orchestrate both organoid formation and the regeneration of intestinal tissue, using an image-based screen to assay an annotated library of compounds. We generate multivariate feature profiles for hundreds of thousands of organoids to quantitatively describe their phenotypic landscape. We then use these phenotypic fingerprints to infer regulatory genetic interactions, establishing a new approach to the mapping of genetic interactions in an emergent system. This allows us to identify genes that regulate cell-fate transitions and maintain the balance between regeneration and homeostasis, unravelling previously unknown roles for several pathways, among them retinoic acid signalling. We then characterize a crucial role for retinoic acid nuclear receptors in controlling exit from the regenerative state and driving enterocyte differentiation. By combining quantitative imaging with RNA sequencing, we show the role of endogenous retinoic acid metabolism in initiating transcriptional programs that guide the cell-fate transitions of intestinal epithelium, and we identify an inhibitor of the retinoid X receptor that improves intestinal regeneration in vivo.

Cooperation between HDAC3 and DAX1 mediates lineage restriction of embryonic stem cells.

Mouse embryonic stem cells (mESCs) are biased toward producing embryonic rather than extraembryonic endoderm fates. Here, we identify the mechanism of this barrier and report that the histone deacetylase Hdac3 and the transcriptional corepressor Dax1 cooperatively limit the lineage repertoire of mESCs by silencing an enhancer of the extraembryonic endoderm-specifying transcription factor Gata6. This restriction is opposed by the pluripotency transcription factors Nr5a2 and Esrrb, which promote cell type conversion. Perturbation of the barrier extends mESC potency and allows formation of 3D spheroids that mimic the spatial segregation of embryonic epiblast and extraembryonic endoderm in early embryos. Overall, this study shows that transcriptional repressors stabilize pluripotency by biasing the equilibrium between embryonic and extraembryonic lineages that is hardwired into the mESC transcriptional network.

Characterizing light-regulated retinal microRNAs reveals rapid turnover as a common property of neuronal microRNAs.

Adaptation to different levels of illumination is central to the function of the retina. Here, we demonstrate that levels of the miR-183/96/182 cluster, miR-204, and miR-211 are regulated by different light levels in the mouse retina. Concentrations of these microRNAs were downregulated during dark adaptation and upregulated in light-adapted retinas, with rapid decay and increased transcription being responsible for the respective changes. We identified the voltage-dependent glutamate transporter Slc1a1 as one of the miR-183/96/182 targets in photoreceptor cells. We found that microRNAs in retinal neurons decay much faster than microRNAs in nonneuronal cells. The high turnover is also characteristic of microRNAs in hippocampal and cortical neurons, and neurons differentiated from ES cells in vitro. Blocking activity reduced turnover of microRNAs in neuronal cells while stimulation with glutamate accelerated it. Our results demonstrate that microRNA metabolism in neurons is higher than in most other cells types and linked to neuronal activity.

Mammalian ISWI and SWI/SNF selectively mediate binding of distinct transcription factors.

Chromatin remodelling complexes evict, slide, insert or replace nucleosomes, which represent an intrinsic barrier for access to DNA. These remodellers function in most aspects of genome utilization including transcription-factor binding, DNA replication and repair1,2. Although they are frequently mutated in cancer3, it remains largely unclear how the four mammalian remodeller families (SWI/SNF, ISWI, CHD and INO80) orchestrate the global organization of nucleosomes. Here we generated viable embryonic stem cells that lack SNF2H, the ATPase of ISWI complexes, enabling study of SNF2H cellular function, and contrast it to BRG1, the ATPase of SWI/SNF. Loss of SNF2H decreases nucleosomal phasing and increases linker lengths, providing in vivo evidence for an ISWI function in ruling nucleosomal spacing in mammals. Systematic analysis of transcription-factor binding reveals that these remodelling activities have specific effects on binding of different transcription factors. One group critically depends on BRG1 and contains the transcriptional repressor REST, whereas a non-overlapping set of transcription factors, including the insulator protein CTCF, relies on SNF2H. This selectivity readily explains why chromosomal folding and insulation of topologically associated domains requires SNF2H, but not BRG1. Collectively, this study shows that mammalian ISWI is critical for nucleosomal periodicity and nuclear organization and that transcription factors rely on specific remodelling pathways for correct genomic binding.

Self-organization and symmetry breaking in intestinal organoid development.

Intestinal organoids are complex three-dimensional structures that mimic the cell-type composition and tissue organization of the intestine by recapitulating the self-organizing ability of cell populations derived from a single intestinal stem cell. Crucial in this process is a first symmetry-breaking event, in which only a fraction of identical cells in a symmetrical sphere differentiate into Paneth cells, which generate the stem-cell niche and lead to asymmetric structures such as the crypts and villi. Here we combine single-cell quantitative genomic and imaging approaches to characterize the development of intestinal organoids from single cells. We show that their development follows a regeneration process that is driven by transient activation of the transcriptional regulator YAP1. Cell-to-cell variability in YAP1, emerging in symmetrical spheres, initiates Notch and DLL1 activation, and drives the symmetry-breaking event and formation of the first Paneth cell. Our findings reveal how single cells exposed to a uniform growth-promoting environment have the intrinsic ability to generate emergent, self-organized behaviour that results in the formation of complex multicellular asymmetric structures.

The Dynamics of mRNA Turnover Revealed by Single-Molecule Imaging in Single Cells.

RNA degradation plays a fundamental role in regulating gene expression. In order to characterize the spatiotemporal dynamics of RNA turnover in single cells, we developed a fluorescent biosensor based on dual-color, single-molecule RNA imaging that allows intact transcripts to be distinguished from stabilized degradation intermediates. Using this method, we measured mRNA decay in single cells and found that individual degradation events occur independently within the cytosol and are not enriched within processing bodies. We show that slicing of an mRNA targeted for endonucleolytic cleavage by the RNA-induced silencing complex can be observed in real time in living cells. This methodology provides a framework for investigating the entire life history of individual mRNAs from birth to death in single cells.

Zfp281 orchestrates interconversion of pluripotent states by engaging Ehmt1 and Zic2.

Developmental cell fate specification is a unidirectional process that can be reverted in response to injury or experimental reprogramming. Whether differentiation and de-differentiation trajectories intersect mechanistically is unclear. Here, we performed comparative screening in lineage-related mouse naïve embryonic stem cells (ESCs) and primed epiblast stem cells (EpiSCs), and identified the constitutively expressed zinc finger transcription factor (TF) Zfp281 as a bidirectional regulator of cell state interconversion. We showed that subtle chromatin binding changes in differentiated cells translate into activation of the histone H3 lysine 9 (H3K9) methyltransferase Ehmt1 and stabilization of the zinc finger TF Zic2 at enhancers and promoters. Genetic gain-of-function and loss-of-function experiments confirmed a critical role of Ehmt1 and Zic2 downstream of Zfp281 both in driving exit from the ESC state and in restricting reprogramming of EpiSCs. Our study reveals that cell type-invariant chromatin association of Zfp281 provides an interaction platform for remodeling the cis-regulatory network underlying cellular plasticity.

Single-cell Analysis Reveals Inter- and Intratumour Heterogeneity in Metastatic Breast Cancer.

Metastasis is the leading cause of cancer-related deaths of breast cancer patients. Some cancer cells in a tumour go through successive steps, referred to as the metastatic cascade, and give rise to metastases at a distant site. We know that the plasticity and heterogeneity of cancer cells play critical roles in metastasis but the precise underlying molecular mechanisms remain elusive. Here we aimed to identify molecular mechanisms of metastasis during colonization, one of the most important yet poorly understood steps of the cascade. We performed single-cell RNA-Seq (scRNA-Seq) on tumours and matched lung macrometastases of patient-derived xenografts of breast cancer. After correcting for confounding factors such as the cell cycle and the percentage of detected genes (PDG), we identified cells in three states in both tumours and metastases. Gene-set enrichment analysis revealed biological processes specific to proliferation and invasion in two states. Our findings suggest that these states are a balance between epithelial-to-mesenchymal (EMT) and mesenchymal-to-epithelial transitions (MET) traits that results in so-called partial EMT phenotypes. Analysis of the top differentially expressed genes (DEGs) between these cell states revealed a common set of partial EMT transcription factors (TFs) controlling gene expression, including ZNF750, OVOL2, TP63, TFAP2C and HEY2. Our data suggest that the TFs related to EMT delineate different cell states in tumours and metastases. The results highlight the marked interpatient heterogeneity of breast cancer but identify common features of single cells from five models of metastatic breast cancer.

monaLisa: an R/Bioconductor package for identifying regulatory motifs.

SUMMARY: Proteins binding to specific nucleotide sequences, such as transcription factors, play key roles in the regulation of gene expression. Their binding can be indirectly observed via associated changes in transcription, chromatin accessibility, DNA methylation and histone modifications. Identifying candidate factors that are responsible for these observed experimental changes is critical to understand the underlying biological processes. Here, we present monaLisa, an R/Bioconductor package that implements approaches to identify relevant transcription factors from experimental data. The package can be easily integrated with other Bioconductor packages and enables seamless motif analyses without any software dependencies outside of R. AVAILABILITY AND IMPLEMENTATION: monaLisa is implemented in R and available on Bioconductor at https://bioconductor.org/packages/monaLisa with the development version hosted on GitHub at https://github.com/fmicompbio/monaLisa. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

OBF1 and Oct factors control the germinal center transcriptional program.

OBF1 is a specific coactivator of the POU family transcription factors OCT1 and OCT2. OBF1 and OCT2 are B cell-specific and indispensable for germinal center (GC) formation, but their mechanism of action is unclear. Here, we show by chromatin immunoprecipitation-sequencing that OBF1 extensively colocalizes with OCT1 and OCT2. We found that these factors also often colocalize with transcription factors of the ETS family. Furthermore, we showed that OBF1, OCT2, and OCT1 bind widely to the promoters or enhancers of genes involved in GC formation in mouse and human GC B cells. Short hairpin RNA knockdown experiments demonstrated that OCT1, OCT2, and OBF1 regulate each other and are essential for proliferation of GC-derived lymphoma cell lines. OBF1 downregulation disrupts the GC transcriptional program: genes involved in GC maintenance, such as BCL6, are downregulated, whereas genes related to exit from the GC program, such as IRF4, are upregulated. Ectopic expression of BCL6 does not restore the proliferation of GC-derived lymphoma cells depleted of OBF1 unless IRF4 is also depleted, indicating that OBF1 controls an essential regulatory node in GC differentiation.

The Hippo kinases LATS1 and 2 control human breast cell fate via crosstalk with ERα.

Cell fate perturbations underlie many human diseases, including breast cancer. Unfortunately, the mechanisms by which breast cell fate are regulated are largely unknown. The mammary gland epithelium consists of differentiated luminal epithelial and basal myoepithelial cells, as well as undifferentiated stem cells and more restricted progenitors. Breast cancer originates from this epithelium, but the molecular mechanisms that underlie breast epithelial hierarchy remain ill-defined. Here, we use a high-content confocal image-based short hairpin RNA screen to identify tumour suppressors that regulate breast cell fate in primary human breast epithelial cells. We show that ablation of the large tumour suppressor kinases (LATS) 1 and 2 (refs 5, 6), which are part of the Hippo pathway, promotes the luminal phenotype and increases the number of bipotent and luminal progenitors, the proposed cells-of-origin of most human breast cancers. Mechanistically, we have identified a direct interaction between Hippo and oestrogen receptor-α (ERα) signalling. In the presence of LATS, ERα was targeted for ubiquitination and Ddb1-cullin4-associated-factor 1 (DCAF1)-dependent proteasomal degradation. Absence of LATS stabilized ERα and the Hippo effectors YAP and TAZ (hereafter YAP/TAZ), which together control breast cell fate through intrinsic and paracrine mechanisms. Our findings reveal a non-canonical (that is, YAP/TAZ-independent) effect of LATS in the regulation of human breast cell fate.

Preprocessing choices affect RNA velocity results for droplet scRNA-seq data.

Experimental single-cell approaches are becoming widely used for many purposes, including investigation of the dynamic behaviour of developing biological systems. Consequently, a large number of computational methods for extracting dynamic information from such data have been developed. One example is RNA velocity analysis, in which spliced and unspliced RNA abundances are jointly modeled in order to infer a 'direction of change' and thereby a future state for each cell in the gene expression space. Naturally, the accuracy and interpretability of the inferred RNA velocities depend crucially on the correctness of the estimated abundances. Here, we systematically compare five widely used quantification tools, in total yielding thirteen different quantification approaches, in terms of their estimates of spliced and unspliced RNA abundances in five experimental droplet scRNA-seq data sets. We show that there are substantial differences between the quantifications obtained from different tools, and identify typical genes for which such discrepancies are observed. We further show that these abundance differences propagate to the downstream analysis, and can have a large effect on estimated velocities as well as the biological interpretation. Our results highlight that abundance quantification is a crucial aspect of the RNA velocity analysis workflow, and that both the definition of the genomic features of interest and the quantification algorithm itself require careful consideration.

PIK3CA(H1047R) induces multipotency and multi-lineage mammary tumours.

The adult mouse mammary epithelium contains self-sustained cell lineages that form the inner luminal and outer basal cell layers, with stem and progenitor cells contributing to its proliferative and regenerative potential. A key issue in breast cancer biology is the effect of genomic lesions in specific mammary cell lineages on tumour heterogeneity and progression. The impact of transforming events on fate conversion in cancer cells of origin and thus their contribution to tumour heterogeneity remains largely elusive. Using in situ genetic lineage tracing and limiting dilution transplantation, we have unravelled the potential of PIK3CA(H1047R), one of the most frequent mutations occurring in human breast cancer, to induce multipotency during tumorigenesis in the mammary gland. Here we show that expression of PIK3CA(H1047R) in lineage-committed basal Lgr5-positive and luminal keratin-8-positive cells of the adult mouse mammary gland evokes cell dedifferentiation into a multipotent stem-like state, suggesting this to be a mechanism involved in the formation of heterogeneous, multi-lineage mammary tumours. Moreover, we show that the tumour cell of origin influences the frequency of malignant mammary tumours. Our results define a key effect of PIK3CA(H1047R) on mammary cell fate in the pre-neoplastic mammary gland and show that the cell of origin of PIK3CA(H1047R) tumours dictates their malignancy, thus revealing a mechanism underlying tumour heterogeneity and aggressiveness.

The CSR-1 endogenous RNAi pathway ensures accurate transcriptional reprogramming during the oocyte-to-embryo transition in Caenorhabditis elegans.

Endogenous RNAi (endoRNAi) is a conserved mechanism for fine-tuning gene expression. In the nematode Caenorhabditis elegans, several endoRNAi pathways are required for the successful development of reproductive cells. The CSR-1 endoRNAi pathway promotes germ cell development, primarily by facilitating the expression of germline genes. In this study, we report a novel function for the CSR-1 pathway in preventing premature activation of embryonic transcription in the developing oocytes, which is accompanied by a general Pol II activation. This CSR-1 function requires its RNase activity, suggesting that, by controlling the levels of maternal mRNAs, CSR-1-dependent endoRNAi contributes to an orderly reprogramming of transcription during the oocyte-to-embryo transition.

Structure of the DDB1-CRBN E3 ubiquitin ligase in complex with thalidomide.

In the 1950s, the drug thalidomide, administered as a sedative to pregnant women, led to the birth of thousands of children with multiple defects. Despite the teratogenicity of thalidomide and its derivatives lenalidomide and pomalidomide, these immunomodulatory drugs (IMiDs) recently emerged as effective treatments for multiple myeloma and 5q-deletion-associated dysplasia. IMiDs target the E3 ubiquitin ligase CUL4-RBX1-DDB1-CRBN (known as CRL4(CRBN)) and promote the ubiquitination of the IKAROS family transcription factors IKZF1 and IKZF3 by CRL4(CRBN). Here we present crystal structures of the DDB1-CRBN complex bound to thalidomide, lenalidomide and pomalidomide. The structure establishes that CRBN is a substrate receptor within CRL4(CRBN) and enantioselectively binds IMiDs. Using an unbiased screen, we identified the homeobox transcription factor MEIS2 as an endogenous substrate of CRL4(CRBN). Our studies suggest that IMiDs block endogenous substrates (MEIS2) from binding to CRL4(CRBN) while the ligase complex is recruiting IKZF1 or IKZF3 for degradation. This dual activity implies that small molecules can modulate an E3 ubiquitin ligase and thereby upregulate or downregulate the ubiquitination of proteins.

The methyltransferase Setdb1 is essential for meiosis and mitosis in mouse oocytes and early embryos.

Oocytes develop the competence for meiosis and early embryogenesis during their growth. Setdb1 is a histone H3 lysine 9 (H3K9) methyltransferase required for post-implantation development and has been implicated in the transcriptional silencing of genes and endogenous retroviral elements (ERVs). To address its role in oogenesis and pre-implantation development, we conditionally deleted Setdb1 in growing oocytes. Loss of Setdb1 expression greatly impaired meiosis. It delayed meiotic resumption, altered the dynamics of chromatin condensation, and impaired kinetochore-spindle interactions, bipolar spindle organization and chromosome segregation in more mature oocytes. The observed phenotypes related to changes in abundance of specific transcripts in mutant oocytes. Setdb1 maternally deficient embryos arrested during pre-implantation development and showed comparable defects during cell cycle progression and in chromosome segregation. Finally, transcriptional profiling data indicate that Setdb1 downregulates rather than silences expression of ERVK and ERVL-MaLR retrotransposons and associated chimearic transcripts during oogenesis. Our results identify Setdb1 as a newly discovered meiotic and embryonic competence factor safeguarding genome integrity at the onset of life.

XRN2 Autoregulation and Control of Polycistronic Gene Expresssion in Caenorhabditis elegans.

XRN2 is a conserved 5'→3' exoribonuclease that complexes with proteins that contain XRN2-binding domains (XTBDs). In Caenorhabditis elegans (C. elegans), the XTBD-protein PAXT-1 stabilizes XRN2 to retain its activity. XRN2 activity is also promoted by 3'(2'),5'-bisphosphate nucleotidase 1 (BPNT1) through hydrolysis of an endogenous XRN inhibitor 3'-phosphoadenosine-5'-phosphate (PAP). Here, we find through unbiased screening that loss of bpnt-1 function suppresses lethality caused by paxt-1 deletion. This unexpected finding is explained by XRN2 autoregulation, which occurs through repression of a cryptic promoter activity and destabilization of the xrn-2 transcript. De-repression appears to be triggered such that more robust XRN2 perturbation, by elimination of both PAXT-1 and BPNT1, is less detrimental to worm viability than absence of PAXT-1 alone. Indeed, we find that two distinct XRN2 repression mechanisms are alleviated at different thresholds of XRN2 inactivation. Like more than 15% of C. elegans genes, xrn-2 occurs in an operon, and we identify additional operons under its control, consistent with a broader function of XRN2 in polycistronic gene regulation. Regulation occurs through intercistronic regions that link genes in an operon, but a part of the mechanisms may allow XRN2 to operate on monocistronic genes in organisms lacking operons.

Silencing of X-Linked MicroRNAs by Meiotic Sex Chromosome Inactivation.

During the pachytene stage of meiosis in male mammals, the X and Y chromosomes are transcriptionally silenced by Meiotic Sex Chromosome Inactivation (MSCI). MSCI is conserved in therian mammals and is essential for normal male fertility. Transcriptomics approaches have demonstrated that in mice, most or all protein-coding genes on the X chromosome are subject to MSCI. However, it is unclear whether X-linked non-coding RNAs behave in a similar manner. The X chromosome is enriched in microRNA (miRNA) genes, with many exhibiting testis-biased expression. Importantly, high expression levels of X-linked miRNAs (X-miRNAs) have been reported in pachytene spermatocytes, indicating that these genes may escape MSCI, and perhaps play a role in the XY-silencing process. Here we use RNA FISH to examine X-miRNA expression in the male germ line. We find that, like protein-coding X-genes, X-miRNAs are expressed prior to prophase I and are thereafter silenced during pachynema. X-miRNA silencing does not occur in mouse models with defective MSCI. Furthermore, X-miRNAs are expressed at pachynema when present as autosomally integrated transgenes. Thus, we conclude that silencing of X-miRNAs during pachynema in wild type males is MSCI-dependent. Importantly, misexpression of X-miRNAs during pachynema causes spermatogenic defects. We propose that MSCI represents a chromosomal mechanism by which X-miRNAs, and other potential X-encoded repressors, can be silenced, thereby regulating genes with critical late spermatogenic functions.