Protocol for mapping cell lineage and cell-type identity of clonally-related cells in situ using MADM-CloneSeq.

The lineage relationship of clonally-related cells offers important insights into the ontogeny and cytoarchitecture of the brain in health and disease. Here, we provide a protocol to concurrently assess cell lineage relationship and cell-type identity among clonally-related cells in situ. We first describe the preparation and screening of acute brain slices containing clonally-related cells labeled using mosaic analysis with double markers (MADM). We then outline steps to collect RNA from individual cells for downstream applications and cell-type identification using RNA sequencing. For complete details on the use and execution of this protocol, please refer to Cheung et al.1.

Multipotent progenitors instruct ontogeny of the superior colliculus.

The superior colliculus (SC) in the mammalian midbrain is essential for multisensory integration and is composed of a rich diversity of excitatory and inhibitory neurons and glia. However, the developmental principles directing the generation of SC cell-type diversity are not understood. Here, we pursued systematic cell lineage tracing in silico and in vivo, preserving full spatial information, using genetic mosaic analysis with double markers (MADM)-based clonal analysis with single-cell sequencing (MADM-CloneSeq). The analysis of clonally related cell lineages revealed that radial glial progenitors (RGPs) in SC are exceptionally multipotent. Individual resident RGPs have the capacity to produce all excitatory and inhibitory SC neuron types, even at the stage of terminal division. While individual clonal units show no pre-defined cellular composition, the establishment of appropriate relative proportions of distinct neuronal types occurs in a PTEN-dependent manner. Collectively, our findings provide an inaugural framework at the single-RGP/-cell level of the mammalian SC ontogeny.

Large neutral amino acid levels tune perinatal neuronal excitability and survival.

Little is known about the critical metabolic changes that neural cells have to undergo during development and how temporary shifts in this program can influence brain circuitries and behavior. Inspired by the discovery that mutations in SLC7A5, a transporter of metabolically essential large neutral amino acids (LNAAs), lead to autism, we employed metabolomic profiling to study the metabolic states of the cerebral cortex across different developmental stages. We found that the forebrain undergoes significant metabolic remodeling throughout development, with certain groups of metabolites showing stage-specific changes, but what are the consequences of perturbing this metabolic program? By manipulating Slc7a5 expression in neural cells, we found that the metabolism of LNAAs and lipids are interconnected in the cortex. Deletion of Slc7a5 in neurons affects the postnatal metabolic state, leading to a shift in lipid metabolism. Additionally, it causes stage- and cell-type-specific alterations in neuronal activity patterns, resulting in a long-term circuit dysfunction.

Tissue-wide genetic and cellular landscape shapes the execution of sequential PRC2 functions in neural stem cell lineage progression.

The generation of a correctly sized cerebral cortex with all-embracing neuronal and glial cell-type diversity critically depends on faithful radial glial progenitor (RGP) cell proliferation/differentiation programs. Temporal RGP lineage progression is regulated by Polycomb repressive complex 2 (PRC2), and loss of PRC2 activity results in severe neurogenesis defects and microcephaly. How PRC2-dependent gene expression instructs RGP lineage progression is unknown. Here, we use mosaic analysis with double markers (MADM)-based single-cell technology and demonstrate that PRC2 is not cell-autonomously required in neurogenic RGPs but rather acts at the global tissue-wide level. Conversely, cortical astrocyte production and maturation is cell-autonomously controlled by PRC2-dependent transcriptional regulation. We thus reveal highly distinct and sequential PRC2 functions in RGP lineage progression that are dependent on complex interplays between intrinsic and tissue-wide properties. In a broader context, our results imply a critical role for the genetic and cellular niche environment in neural stem cell behavior.

Simultaneous brain cell type and lineage determined by scRNA-seq reveals stereotyped cortical development.

Mutations are acquired frequently, such that each cell's genome inscribes its history of cell divisions. Common genomic alterations involve loss of heterozygosity (LOH). LOH accumulates throughout the genome, offering large encoding capacity for inferring cell lineage. Using only single-cell RNA sequencing (scRNA-seq) of mouse brain cells, we found that LOH events spanning multiple genes are revealed as tracts of monoallelically expressed, constitutionally heterozygous single-nucleotide variants (SNVs). We simultaneously inferred cell lineage and marked developmental time points based on X chromosome inactivation and the total number of LOH events while identifying cell types from gene expression patterns. Our results are consistent with progenitor cells giving rise to multiple cortical cell types through stereotyped expansion and distinct waves of neurogenesis. This type of retrospective analysis could be incorporated into scRNA-seq pipelines and, compared with experimental approaches for determining lineage in model organisms, is applicable where genetic engineering is prohibited, such as humans.

Tissue-Wide Effects Override Cell-Intrinsic Gene Function in Radial Neuron Migration.

The mammalian neocortex is composed of diverse neuronal and glial cell classes that broadly arrange in six distinct laminae. Cortical layers emerge during development and defects in the developmental programs that orchestrate cortical lamination are associated with neurodevelopmental diseases. The developmental principle of cortical layer formation depends on concerted radial projection neuron migration, from their birthplace to their final target position. Radial migration occurs in defined sequential steps, regulated by a large array of signaling pathways. However, based on genetic loss-of-function experiments, most studies have thus far focused on the role of cell-autonomous gene function. Yet, cortical neuron migration in situ is a complex process and migrating neurons traverse along diverse cellular compartments and environments. The role of tissue-wide properties and genetic state in radial neuron migration is however not clear. Here we utilized mosaic analysis with double markers (MADM) technology to either sparsely or globally delete gene function, followed by quantitative single-cell phenotyping. The MADM-based gene ablation paradigms in combination with computational modeling demonstrated that global tissue-wide effects predominate cell-autonomous gene function albeit in a gene-specific manner. Our results thus suggest that the genetic landscape in a tissue critically affects the overall migration phenotype of individual cortical projection neurons. In a broader context, our findings imply that global tissue-wide effects represent an essential component of the underlying etiology associated with focal malformations of cortical development in particular, and neurological diseases in general.

An inhibitor-mediated beta-cell dedifferentiation model reveals distinct roles for FoxO1 in glucagon repression and insulin maturation.

OBJECTIVE: The loss of forkhead box protein O1 (FoxO1) signaling in response to metabolic stress contributes to the etiology of type II diabetes, causing the dedifferentiation of pancreatic beta cells to a cell type reminiscent of endocrine progenitors. Lack of methods to easily model this process in vitro, however, have hindered progress into the identification of key downstream targets and potential inhibitors. We therefore aimed to establish such an in vitro cellular dedifferentiation model and apply it to identify novel agents involved in the maintenance of beta-cell identity. METHODS: The murine beta-cell line, Min6, was used for primary experiments and high-content screening. Screens encompassed a library of small-molecule drugs representing the chemical and target space of all FDA-approved small molecules with an automated immunofluorescence readout. Validation experiments were performed in a murine alpha-cell line as well as in primary murine and human diabetic islets. Developmental effects were studied in zebrafish and C. elegans models, while diabetic db/db mouse models were used to elucidate global glucose metabolism outcomes. RESULTS: We show that short-term pharmacological FoxO1 inhibition can model beta-cell dedifferentiation by downregulating beta-cell-specific transcription factors, resulting in the aberrant expression of progenitor genes and the alpha-cell marker glucagon. From a high-content screen, we identified loperamide as a small molecule that can prevent FoxO inhibitor-induced glucagon expression and further stimulate insulin protein processing and secretion by altering calcium levels, intracellular pH, and FoxO1 localization. CONCLUSIONS: Our study provides novel models, molecular targets, and drug candidates for studying and preventing beta-cell dedifferentiation.

LINC01133 Inhibits Invasion and Promotes Proliferation in an Endometriosis Epithelial Cell Line.

Endometriosis is a common gynecological disorder characterized by ectopic growth of endometrium outside the uterus and is associated with chronic pain and infertility. We investigated the role of the long intergenic noncoding RNA 01133 (LINC01133) in endometriosis, an lncRNA that has been implicated in several types of cancer. We found that LINC01133 is upregulated in ectopic endometriotic lesions. As expression appeared higher in the epithelial endometrial layer, we performed a siRNA knockdown of LINC01133 in an endometriosis epithelial cell line. Phenotypic assays indicated that LINC01133 may promote proliferation and suppress cellular migration, and affect the cytoskeleton and morphology of the cells. Gene ontology analysis of differentially expressed genes indicated that cell proliferation and migration pathways were affected in line with the observed phenotype. We validated upregulation of p21 and downregulation of Cyclin A at the protein level, which together with the quantification of the DNA content using fluorescence-activated cell sorting (FACS) analysis indicated that the observed effects on cellular proliferation may be due to changes in cell cycle. Further, we found testis-specific protein kinase 1 (TESK1) kinase upregulation corresponding with phosphorylation and inactivation of actin severing protein Cofilin, which could explain changes in the cytoskeleton and cellular migration. These results indicate that endometriosis is associated with LINC01133 upregulation, which may affect pathogenesis via the cellular proliferation and migration pathways.

Inducible uniparental chromosome disomy to probe genomic imprinting at single-cell level in brain and beyond.

Genomic imprinting is an epigenetic mechanism that results in parental allele-specific expression of ~1% of all genes in mouse and human. Imprinted genes are key developmental regulators and play pivotal roles in many biological processes such as nutrient transfer from the mother to offspring and neuronal development. Imprinted genes are also involved in human disease, including neurodevelopmental disorders, and often occur in clusters that are regulated by a common imprint control region (ICR). In extra-embryonic tissues ICRs can act over large distances, with the largest surrounding Igf2r spanning over 10 million base-pairs. Besides classical imprinted expression that shows near exclusive maternal or paternal expression, widespread biased imprinted expression has been identified mainly in brain. In this review we discuss recent developments mapping cell type specific imprinted expression in extra-embryonic tissues and neocortex in the mouse. We highlight the advantages of using an inducible uniparental chromosome disomy (UPD) system to generate cells carrying either two maternal or two paternal copies of a specific chromosome to analyze the functional consequences of genomic imprinting. Mosaic Analysis with Double Markers (MADM) allows fluorescent labeling and concomitant induction of UPD sparsely in specific cell types, and thus to over-express or suppress all imprinted genes on that chromosome. To illustrate the utility of this technique, we explain how MADM-induced UPD revealed new insights about the function of the well-studied Cdkn1c imprinted gene, and how MADM-induced UPDs led to identification of highly cell type specific phenotypes related to perturbed imprinted expression in the mouse neocortex. Finally, we give an outlook on how MADM could be used to probe cell type specific imprinted expression in other tissues in mouse, particularly in extra-embryonic tissues.

Generation and isolation of single cells from mouse brain with mosaic analysis with double markers-induced uniparental chromosome disomy.

Mosaic analysis with double markers (MADM) technology enables concomitant fluorescent cell labeling and induction of uniparental chromosome disomy (UPD) with single-cell resolution. In UPD, imprinted genes are either overexpressed 2-fold or are not expressed. Here, the MADM platform is utilized to probe imprinting phenotypes at the transcriptional level. This protocol highlights major steps for the generation and isolation of projection neurons and astrocytes with MADM-induced UPD from mouse cerebral cortex for downstream single-cell and low-input sample RNA-sequencing experiments. For complete details on the use and execution of this protocol, please refer to Laukoter et al. (2020b).

Correction: The Airn lncRNA does not require any DNA elements within its locus to silence distant imprinted genes.

[This corrects the article DOI: 10.1371/journal.pgen.1008268.].

Cell-Type Specificity of Genomic Imprinting in Cerebral Cortex.

In mammalian genomes, a subset of genes is regulated by genomic imprinting, resulting in silencing of one parental allele. Imprinting is essential for cerebral cortex development, but prevalence and functional impact in individual cells is unclear. Here, we determined allelic expression in cortical cell types and established a quantitative platform to interrogate imprinting in single cells. We created cells with uniparental chromosome disomy (UPD) containing two copies of either the maternal or the paternal chromosome; hence, imprinted genes will be 2-fold overexpressed or not expressed. By genetic labeling of UPD, we determined cellular phenotypes and transcriptional responses to deregulated imprinted gene expression at unprecedented single-cell resolution. We discovered an unexpected degree of cell-type specificity and a novel function of imprinting in the regulation of cortical astrocyte survival. More generally, our results suggest functional relevance of imprinted gene expression in glial astrocyte lineage and thus for generating cortical cell-type diversity.

Imprinted Cdkn1c genomic locus cell-autonomously promotes cell survival in cerebral cortex development.

The cyclin-dependent kinase inhibitor p57KIP2 is encoded by the imprinted Cdkn1c locus, exhibits maternal expression, and is essential for cerebral cortex development. How Cdkn1c regulates corticogenesis is however not clear. To this end we employ Mosaic Analysis with Double Markers (MADM) technology to genetically dissect Cdkn1c gene function in corticogenesis at single cell resolution. We find that the previously described growth-inhibitory Cdkn1c function is a non-cell-autonomous one, acting on the whole organism. In contrast we reveal a growth-promoting cell-autonomous Cdkn1c function which at the mechanistic level mediates radial glial progenitor cell and nascent projection neuron survival. Strikingly, the growth-promoting function of Cdkn1c is highly dosage sensitive but not subject to genomic imprinting. Collectively, our results suggest that the Cdkn1c locus regulates cortical development through distinct cell-autonomous and non-cell-autonomous mechanisms. More generally, our study highlights the importance to probe the relative contributions of cell intrinsic gene function and tissue-wide mechanisms to the overall phenotype.

The Airn lncRNA does not require any DNA elements within its locus to silence distant imprinted genes.

Long non-coding (lnc) RNAs are numerous and found throughout the mammalian genome, and many are thought to be involved in the regulation of gene expression. However, the majority remain relatively uncharacterised and of uncertain function making the use of model systems to uncover their mode of action valuable. Imprinted lncRNAs target and recruit epigenetic silencing factors to a cluster of imprinted genes on the same chromosome, making them one of the best characterized lncRNAs for silencing distant genes in cis. In this study we examined silencing of the distant imprinted gene Slc22a3 by the lncRNA Airn in the Igf2r imprinted cluster in mouse. Previously we proposed that imprinted lncRNAs may silence distant imprinted genes by disrupting promoter-enhancer interactions by being transcribed through the enhancer, which we called the enhancer interference hypothesis. Here we tested this hypothesis by first using allele-specific chromosome conformation capture (3C) to detect interactions between the Slc22a3 promoter and the locus of the Airn lncRNA that silences it on the paternal chromosome. In agreement with the model, we found interactions enriched on the maternal allele across the entire Airn gene consistent with multiple enhancer-promoter interactions. Therefore, to test the enhancer interference hypothesis we devised an approach to delete the entire Airn gene. However, the deletion showed that there are no essential enhancers for Slc22a2, Pde10a and Slc22a3 within the Airn gene, strongly indicating that the Airn RNA rather than its transcription is responsible for silencing distant imprinted genes. Furthermore, we found that silent imprinted genes were covered with large blocks of H3K27me3 on the repressed paternal allele. Therefore we propose an alternative hypothesis whereby the chromosome interactions may initially guide the lncRNA to target imprinted promoters and recruit repressive chromatin, and that these interactions are lost once silencing is established.

Norepinephrine triggers an immediate-early regulatory network response in primary human white adipocytes.

BACKGROUND: Norepinephrine (NE) signaling has a key role in white adipose tissue (WAT) functions, including lipolysis, free fatty acid liberation and, under certain conditions, conversion of white into brite (brown-in-white) adipocytes. However, acute effects of NE stimulation have not been described at the transcriptional network level. RESULTS: We used RNA-seq to uncover a broad transcriptional response. The inference of protein-protein and protein-DNA interaction networks allowed us to identify a set of immediate-early genes (IEGs) with high betweenness, validating our approach and suggesting a hierarchical control of transcriptional regulation. In addition, we identified a transcriptional regulatory network with IEGs as master regulators, including HSF1 and NFIL3 as novel NE-induced IEG candidates. Moreover, a functional enrichment analysis and gene clustering into functional modules suggest a crosstalk between metabolic, signaling, and immune responses. CONCLUSIONS: Altogether, our network biology approach explores for the first time the immediate-early systems level response of human adipocytes to acute sympathetic activation, thereby providing a first network basis of early cell fate programs and crosstalks between metabolic and transcriptional networks required for proper WAT function.

Mapping the mouse Allelome reveals tissue-specific regulation of allelic expression.

To determine the dynamics of allelic-specific expression during mouse development, we analyzed RNA-seq data from 23 F1 tissues from different developmental stages, including 19 female tissues allowing X chromosome inactivation (XCI) escapers to also be detected. We demonstrate that allelic expression arising from genetic or epigenetic differences is highly tissue-specific. We find that tissue-specific strain-biased gene expression may be regulated by tissue-specific enhancers or by post-transcriptional differences in stability between the alleles. We also find that escape from X-inactivation is tissue-specific, with leg muscle showing an unexpectedly high rate of XCI escapers. By surveying a range of tissues during development, and performing extensive validation, we are able to provide a high confidence list of mouse imprinted genes including 18 novel genes. This shows that cluster size varies dynamically during development and can be substantially larger than previously thought, with the Igf2r cluster extending over 10 Mb in placenta.

Mosaic Analysis with Double Markers Reveals Distinct Sequential Functions of Lgl1 in Neural Stem Cells.

The concerted production of neurons and glia by neural stem cells (NSCs) is essential for neural circuit assembly. In the developing cerebral cortex, radial glia progenitors (RGPs) generate nearly all neocortical neurons and certain glia lineages. RGP proliferation behavior shows a high degree of non-stochasticity, thus a deterministic characteristic of neuron and glia production. However, the cellular and molecular mechanisms controlling RGP behavior and proliferation dynamics in neurogenesis and glia generation remain unknown. By using mosaic analysis with double markers (MADM)-based genetic paradigms enabling the sparse and global knockout with unprecedented single-cell resolution, we identified Lgl1 as a critical regulatory component. We uncover Lgl1-dependent tissue-wide community effects required for embryonic cortical neurogenesis and novel cell-autonomous Lgl1 functions controlling RGP-mediated glia genesis and postnatal NSC behavior. These results suggest that NSC-mediated neuron and glia production is tightly regulated through the concerted interplay of sequential Lgl1-dependent global and cell intrinsic mechanisms.

Artemisinins Target GABAA Receptor Signaling and Impair α Cell Identity.

Type 1 diabetes is characterized by the destruction of pancreatic ß cells, and generating new insulin-producing cells from other cell types is a major aim of regenerative medicine. One promising approach is transdifferentiation of developmentally related pancreatic cell types, including glucagon-producing α cells. In a genetic model, loss of the master regulatory transcription factor Arx is sufficient to induce the conversion of α cells to functional ß-like cells. Here, we identify artemisinins as small molecules that functionally repress Arx by causing its translocation to the cytoplasm. We show that the protein gephyrin is the mammalian target of these antimalarial drugs and that the mechanism of action of these molecules depends on the enhancement of GABAA receptor signaling. Our results in zebrafish, rodents, and primary human pancreatic islets identify gephyrin as a druggable target for the regeneration of pancreatic ß cell mass from α cells.

Long non-coding RNAs display higher natural expression variation than protein-coding genes in healthy humans.

BACKGROUND: Long non-coding RNAs (lncRNAs) are increasingly implicated as gene regulators and may ultimately be more numerous than protein-coding genes in the human genome. Despite large numbers of reported lncRNAs, reference annotations are likely incomplete due to their lower and tighter tissue-specific expression compared to mRNAs. An unexplored factor potentially confounding lncRNA identification is inter-individual expression variability. Here, we characterize lncRNA natural expression variability in human primary granulocytes. RESULTS: We annotate granulocyte lncRNAs and mRNAs in RNA-seq data from 10 healthy individuals, identifying multiple lncRNAs absent from reference annotations, and use this to investigate three known features (higher tissue-specificity, lower expression, and reduced splicing efficiency) of lncRNAs relative to mRNAs. Expression variability was examined in seven individuals sampled three times at 1- or more than 1-month intervals. We show that lncRNAs display significantly more inter-individual expression variability compared to mRNAs. We confirm this finding in two independent human datasets by analyzing multiple tissues from the GTEx project and lymphoblastoid cell lines from the GEUVADIS project. Using the latter dataset we also show that including more human donors into the transcriptome annotation pipeline allows identification of an increasing number of lncRNAs, but minimally affects mRNA gene number. CONCLUSIONS: A comprehensive annotation of lncRNAs is known to require an approach that is sensitive to low and tight tissue-specific expression. Here we show that increased inter-individual expression variability is an additional general lncRNA feature to consider when creating a comprehensive annotation of human lncRNAs or proposing their use as prognostic or disease markers.

A human haploid gene trap collection to study lncRNAs with unusual RNA biology.

Many thousand long non-coding (lnc) RNAs are mapped in the human genome. Time consuming studies using reverse genetic approaches by post-transcriptional knock-down or genetic modification of the locus demonstrated diverse biological functions for a few of these transcripts. The Human Gene Trap Mutant Collection in haploid KBM7 cells is a ready-to-use tool for studying protein-coding gene function. As lncRNAs show remarkable differences in RNA biology compared to protein-coding genes, it is unclear if this gene trap collection is useful for functional analysis of lncRNAs. Here we use the uncharacterized LOC100288798 lncRNA as a model to answer this question. Using public RNA-seq data we show that LOC100288798 is ubiquitously expressed, but inefficiently spliced. The minor spliced LOC100288798 isoforms are exported to the cytoplasm, whereas the major unspliced isoform is nuclear localized. This shows that LOC100288798 RNA biology differs markedly from typical mRNAs. De novo assembly from RNA-seq data suggests that LOC100288798 extends 289kb beyond its annotated 3' end and overlaps the downstream SLC38A4 gene. Three cell lines with independent gene trap insertions in LOC100288798 were available from the KBM7 gene trap collection. RT-qPCR and RNA-seq confirmed successful lncRNA truncation and its extended length. Expression analysis from RNA-seq data shows significant deregulation of 41 protein-coding genes upon LOC100288798 truncation. Our data shows that gene trap collections in human haploid cell lines are useful tools to study lncRNAs, and identifies the previously uncharacterized LOC100288798 as a potential gene regulator.