Structurally divergent and recurrently mutated regions of primate genomes

We sequenced and assembled using multiple long-read sequencing technologies the genomes of chimpanzee, bonobo, gorilla, orangutan, gibbon, macaque, owl monkey, and marmoset. We identified 1,338,997 lineage-specific fixed structural variants (SVs) disrupting 1,561 protein-coding genes and 136,932 regulatory elements, including the most complete set of human-specific fixed differences. We estimate that 819.47 Mbp or ∼27% of the genome has been affected by SVs across primate evolution. We identify 1,607 structurally divergent regions wherein recurrent structural variation contributes to creating SV hotspots where genes are recurrently lost (e.g., CARD, C4, and OLAH gene families) and additional lineage-specific genes are generated (e.g., CKAP2, VPS36, ACBD7, and NEK5 paralogs), becoming targets of rapid chromosomal diversification and positive selection (e.g., RGPD gene family). High-fidelity long-read sequencing has made these dynamic regions of the genome accessible for sequence-level analyses within and between primate species.

Thalamocortical organoids enable in vitro modeling of 22q11.2 microdeletion associated with neuropsychiatric disorders

Thalamic dysfunction has been implicated in multiple psychiatric disorders. We sought to study the mechanisms by which abnormalities emerge in the context of the 22q11.2 microdeletion, which confers significant genetic risk for psychiatric disorders. We investigated early stages of human thalamus development using human pluripotent stem cell-derived organoids and show that the 22q11.2 microdeletion underlies widespread transcriptional dysregulation associated with psychiatric disorders in thalamic neurons and glia, including elevated expression of FOXP2. Using an organoid co-culture model, we demonstrate that the 22q11.2 microdeletion mediates an overgrowth of thalamic axons in a FOXP2-dependent manner. Finally, we identify ROBO2 as a candidate molecular mediator of the effects of FOXP2 overexpression on thalamic axon overgrowth. Together, our study suggests that early steps in thalamic development are dysregulated in a model of genetic risk for schizophrenia and contribute to neural phenotypes in 22q11.2 deletion syndrome.

Human neuronal maturation comes of age: cellular mechanisms and species differences

The delayed and prolonged postmitotic maturation of human neurons, compared with neurons from other species, may contribute to human-specific cognitive abilities and neurological disorders. Here we review the mechanisms of neuronal maturation, applying lessons from model systems to understand the specific features of protracted human cortical maturation and species differences. We cover cell-intrinsic features of neuronal maturation, including transcriptional, epigenetic and metabolic mechanisms, as well as cell-extrinsic features, including the roles of activity and synapses, the actions of glial cells and the contribution of the extracellular matrix. We discuss evidence for species differences in biochemical reaction rates, the proposed existence of an epigenetic maturation clock and the contributions of both general and modular mechanisms to species-specific maturation timing. Finally, we suggest approaches to measure, improve and accelerate the maturation of human neurons in culture, examine crosstalk and interactions among these different aspects of maturation and propose conceptual models to guide future studies.

Mapping cis- and trans-regulatory target genes of human-specific deletions

Deletion of functional sequence is predicted to represent a fundamental mechanism of molecular evolution1,2. Comparative genetic studies of primates2,3 have identified thousands of human-specific deletions (hDels), and the cis-regulatory potential of short (≤31 base pairs) hDels has been assessed using reporter assays4. However, how structural variant-sized (≥50 base pairs) hDels influence molecular and cellular processes in their native genomic contexts remains unexplored. Here, we design genome-scale libraries of single-guide RNAs targeting 7.2 megabases of sequence in 6,358 hDels and present a systematic CRISPR interference (CRISPRi) screening approach to identify hDels that modify cellular proliferation in chimpanzee pluripotent stem cells. By intersecting hDels with chromatin state features and performing single-cell CRISPRi (Perturb-seq) to identify their cis- and trans-regulatory target genes, we discovered 19 hDels controlling gene expression. We highlight two hDels, hDel_2247 and hDel_585, with tissue-specific activity in the liver and brain, respectively. Our findings reveal a molecular and cellular role for sequences lost in the human lineage and establish a framework for functionally interrogating human-specific genetic variants.

A cross-species proteomic map reveals neoteny of human synapse development

The molecular mechanisms and evolutionary changes accompanying synapse development are still poorly understood1,2. Here we generate a cross-species proteomic map of synapse development in the human, macaque and mouse neocortex. By tracking the changes of more than 1,000 postsynaptic density (PSD) proteins from midgestation to young adulthood, we find that PSD maturation in humans separates into three major phases that are dominated by distinct pathways. Cross-species comparisons reveal that human PSDs mature about two to three times slower than those of other species and contain higher levels of Rho guanine nucleotide exchange factors (RhoGEFs) in the perinatal period. Enhancement of RhoGEF signalling in human neurons delays morphological maturation of dendritic spines and functional maturation of synapses, potentially contributing to the neotenic traits of human brain development. In addition, PSD proteins can be divided into four modules that exert stage- and cell-type-specific functions, possibly explaining their differential associations with cognitive functions and diseases. Our proteomic map of synapse development provides a blueprint for studying the molecular basis and evolutionary changes of synapse maturation.

Human-specific genetics: new tools to explore the molecular and cellular basis of human evolution

Our ancestors acquired morphological, cognitive and metabolic modifications that enabled humans to colonize diverse habitats, develop extraordinary technologies and reshape the biosphere. Understanding the genetic, developmental and molecular bases for these changes will provide insights into how we became human. Connecting human-specific genetic changes to species differences has been challenging owing to an abundance of low-effect size genetic changes, limited descriptions of phenotypic differences across development at the level of cell types and lack of experimental models. Emerging approaches for single-cell sequencing, genetic manipulation and stem cell culture now support descriptive and functional studies in defined cell types with a human or ape genetic background. In this Review, we describe how the sequencing of genomes from modern and archaic hominins, great apes and other primates is revealing human-specific genetic changes and how new molecular and cellular approaches - including cell atlases and organoids - are enabling exploration of the candidate causal factors that underlie human-specific traits.

Comparative landscape of genetic dependencies in human and chimpanzee stem cells

Comparative studies of great apes provide a window into our evolutionary past, but the extent and identity of cellular differences that emerged during hominin evolution remain largely unexplored. We established a comparative loss-of-function approach to evaluate whether human cells exhibit distinct genetic dependencies. By performing genome-wide CRISPR interference screens in human and chimpanzee pluripotent stem cells, we identified 75 genes with species-specific effects on cellular proliferation. These genes comprised coherent processes, including cell-cycle progression and lysosomal signaling, which we determined to be human-derived by comparison with orangutan cells. Human-specific robustness to CDK2 and CCNE1 depletion persisted in neural progenitor cells and cerebral organoids, supporting the G1-phase length hypothesis as a potential evolutionary mechanism in human brain expansion. Our findings demonstrate that evolutionary changes in human cells reshaped the landscape of essential genes and establish a platform for systematically uncovering latent cellular and molecular differences between species.

Genetic architecture of human brain evolution

Comparative studies of hominids have long sought to identify mutational events that shaped the evolution of the human nervous system. However, functional genetic differences are outnumbered by millions of nearly neutral mutations, and the developmental mechanisms underlying human nervous system specializations are difficult to model and incompletely understood. Candidate-gene studies have attempted to map select human-specific genetic differences to neurodevelopmental functions, but it remains unclear how to contextualize the relative effects of genes that are investigated independently. Considering these limitations, we discuss scalable approaches for probing the functional contributions of human-specific genetic differences. We propose that a systems-level view will enable a more quantitative and integrative understanding of the genetic, molecular and cellular underpinnings of human nervous system evolution.

Comparative landscape of genetic dependencies in human and chimpanzee stem cells

Comparative studies of great apes provide a window into our evolutionary past, but the extent and identity of cellular differences that emerged during hominin evolution remain largely unexplored. We established a comparative loss-of-function approach to evaluate whether changes in human cells alter requirements for essential genes. By performing genome-wide CRISPR interference screens in human and chimpanzee pluripotent stem cells, we identified 75 genes with species-specific effects on cellular proliferation. These genes comprised coherent processes, including cell cycle progression and lysosomal signaling, which we determined to be human-derived by comparison with orangutan cells. Human-specific robustness to CDK2 and CCNE1 depletion persisted in neural progenitor cells, providing support for the G1-phase length hypothesis as a potential evolutionary mechanism in human brain expansion. Our findings demonstrate that evolutionary changes in human cells can reshape the landscape of essential genes and establish a platform for systematically uncovering latent cellular and molecular differences between species.

Machine learning dissection of human accelerated regions in primate neurodevelopment

Using machine learning (ML), we interrogated the function of all human-chimpanzee variants in 2,645 human accelerated regions (HARs), finding 43% of HARs have variants with large opposing effects on chromatin state and 14% on neurodevelopmental enhancer activity. This pattern, consistent with compensatory evolution, was confirmed using massively parallel reporter assays in chimpanzee and human neural progenitor cells. The species-specific enhancer activity of HARs was accurately predicted from the presence and absence of transcription factor footprints in each species. Despite these striking cis effects, activity of a given HAR sequence was nearly identical in human and chimpanzee cells. This suggests that HARs did not evolve to compensate for changes in the trans environment but instead altered their ability to bind factors present in both species. Thus, ML prioritized variants with functional effects on human neurodevelopment and revealed an unexpected reason why HARs may have evolved so rapidly.

Filling an ARHGAP in our knowledge of human brain evolution

The human cerebral cortex has tripled in size since our divergence from a common ancestor with chimpanzees. This cortical expansion is driven by the increased proliferative capacity of radial glia (RG), a neural progenitor cell (NPC) population that generates cortical neurons. RG along the ventricular zone (VZ) produce neurons and also give rise to basal progenitors (BPs), which migrate to the embryonic subventricular zone (SVZ). Comparative studies suggest that the increased proliferative capacity of human NPCs involves cell-intrinsic mechanisms (Otani et al, 2016), and a number of human-specific genetic changes have recently been linked to NPC proliferation. In particular, overexpression studies in model organisms indicate that the human-specific gene ARHGAP11B is sufficient to increase BP abundance when introduced into the developing brain of non-human model organisms (Florio et al, 2015; Kalebic et al, 2018; Heide et al, 2020). However, studying human-specific mutations in a hominid genetic and developmental context, rather than in more divergent model organisms, could provide further insight into the evolutionary consequences and effect size of human mutations. Fischer et al (2022) now developed a novel organoid electroporation technique to establish the necessity and sufficiency of ARHGAP11B for BP proliferation in cells from humans and our closest living relative, chimpanzees (Fig 1).

IoT cloud laboratory: Internet of Things architecture for cellular biology

The Internet of Things (IoT) provides a simple framework to control online devices easily. IoT is now a commonplace tool used by technology companies but is rarely used in biology experiments. IoT can benefit cloud biology research through alarm notifications, automation, and the real-time monitoring of experiments. We developed an IoT architecture to control biological devices and implemented it in lab experiments. Lab devices for electrophysiology, microscopy, and microfluidics were created from the ground up to be part of a unified IoT architecture. The system allows each device to be monitored and controlled from an online web tool. We present our IoT architecture so other labs can replicate it for their own experiments.

Tropism of SARS-CoV-2 for human cortical astrocytes

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) readily infects a variety of cell types impacting the function of vital organ systems, with particularly severe impact on respiratory function. Neurological symptoms, which range in severity, accompany as many as one-third of COVID-19 cases, indicating a potential vulnerability of neural cell types. To assess whether human cortical cells can be directly infected by SARS-CoV-2, we utilized stem-cell-derived cortical organoids as well as primary human cortical tissue, both from developmental and adult stages. We find significant and predominant infection in cortical astrocytes in both primary tissue and organoid cultures, with minimal infection of other cortical populations. Infected and bystander astrocytes have a corresponding increase in inflammatory gene expression, reactivity characteristics, increased cytokine and growth factor signaling, and cellular stress. Although human cortical cells, particularly astrocytes, have no observable ACE2 expression, we find high levels of coronavirus coreceptors in infected astrocytes, including CD147 and DPP4. Decreasing coreceptor abundance and activity reduces overall infection rate, and increasing expression is sufficient to promote infection. Thus, we find tropism of SARS-CoV-2 for human astrocytes resulting in inflammatory gliosis-type injury that is dependent on coronavirus coreceptors.

The development and evolution of inhibitory neurons in primate cerebrum

Neuroanatomists have long speculated that expanded primate brains contain an increased morphological diversity of inhibitory neurons (INs)¹, and recent studies have identified primate-specific neuronal populations at the molecular level². However, we know little about the developmental mechanisms that specify evolutionarily novel cell types in the brain. Here, we reconstruct gene expression trajectories specifying INs generated throughout the neurogenic period in macaques and mice by analysing the transcriptomes of 250,181 cells. We find that the initial classes of INs generated prenatally are largely conserved among mammals. Nonetheless, we identify two contrasting developmental mechanisms for specifying evolutionarily novel cell types during prenatal development. First, we show that recently identified primate-specific TAC3 striatal INs are specified by a unique transcriptional programme in progenitors followed by induction of a distinct suite of neuropeptides and neurotransmitter receptors in new-born neurons. Second, we find that multiple classes of transcriptionally conserved olfactory bulb (OB)-bound precursors are redirected to expanded primate white matter and striatum. These classes include a novel peristriatal class of striatum laureatum neurons that resemble dopaminergic periglomerular cells of the OB. We propose an evolutionary model in which conserved initial classes of neurons supplying the smaller primate OB are reused in the enlarged striatum and cortex. Together, our results provide a unified developmental taxonomy of initial classes of mammalian INs and reveal multiple developmental mechanisms for neural cell type evolution.

Evolutionary modelling shows that an initial set of inhibitory neurons serving olfactory bulbs may have been repurposed to diversify the taxonomy of interneurons found in the expanded striata and cortices in primates.

Light-weight electrophysiology hardware and software platform for cloud-based neural recording experiments

Objective.Neural activity represents a functional readout of neurons that is increasingly important to monitor in a wide range of experiments. Extracellular recordings have emerged as a powerful technique for measuring neural activity because these methods do not lead to the destruction or degradation of the cells being measured. Current approaches to electrophysiology have a low throughput of experiments due to manual supervision and expensive equipment. This bottleneck limits broader inferences that can be achieved with numerous long-term recorded samples.Approach.We developed Piphys, an inexpensive open source neurophysiological recording platform that consists of both hardware and software. It is easily accessed and controlled via a standard web interface through Internet of Things (IoT) protocols.Main results.We used a Raspberry Pi as the primary processing device along with an Intan bioamplifier. We designed a hardware expansion circuit board and software to enable voltage sampling and user interaction. This standalone system was validated with primary human neurons, showing reliability in collecting neural activity in near real-time.Significance.The hardware modules and cloud software allow for remote control of neural recording experiments as well as horizontal scalability, enabling long-term observations of development, organization, and neural activity at scale.

Picroscope: low-cost system for simultaneous longitudinal biological imaging

Simultaneous longitudinal imaging across multiple conditions and replicates has been crucial for scientific studies aiming to understand biological processes and disease. Yet, imaging systems capable of accomplishing these tasks are economically unattainable for most academic and teaching laboratories around the world. Here, we propose the Picroscope, which is the first low-cost system for simultaneous longitudinal biological imaging made primarily using off-the-shelf and 3D-printed materials. The Picroscope is compatible with standard 24-well cell culture plates and captures 3D z-stack image data. The Picroscope can be controlled remotely, allowing for automatic imaging with minimal intervention from the investigator. Here, we use this system in a range of applications. We gathered longitudinal whole organism image data for frogs, zebrafish, and planaria worms. We also gathered image data inside an incubator to observe 2D monolayers and 3D mammalian tissue culture models. Using this tool, we can measure the behavior of entire organisms or individual cells over long-time periods.

Distinct nuclear compartment-associated genome architecture in the developing mammalian brain

Nuclear compartments are thought to play a role in three-dimensional genome organization and gene expression. In mammalian brain, the architecture and dynamics of nuclear compartment-associated genome organization is not known. In this study, we developed Genome Organization using CUT and RUN Technology (GO-CaRT) to map genomic interactions with two nuclear compartments-the nuclear lamina and nuclear speckles-from different regions of the developing mouse, macaque and human brain. Lamina-associated domain (LAD) architecture in cells in vivo is distinct from that of cultured cells, including major differences in LADs previously considered to be cell type invariant. In the mouse and human forebrain, dorsal and ventral neural precursor cells have differences in LAD architecture that correspond to their regional identity. LADs in the human and mouse cortex contain transcriptionally highly active sub-domains characterized by broad depletion of histone-3-lysine-9 dimethylation. Evolutionarily conserved LADs in human, macaque and mouse brain are enriched for transcriptionally active neural genes associated with synapse function. By integrating GO-CaRT maps with genome-wide association study data, we found speckle-associated domains to be enriched for schizophrenia risk loci, indicating a physical relationship between these disease-associated genetic variants and a specific nuclear structure. Our work provides a framework for understanding the relationship between distinct nuclear compartments and genome function in brain development and disease.

UCSC Cell Browser: Visualize Your Single-Cell Data

Summary

As the use of single-cell technologies has grown, so has the need for tools to explore these large, complicated datasets. The UCSC Cell Browser is a tool that allows scientists to visualize gene expression and metadata annotation distribution throughout a single-cell dataset or multiple datasets.

Availability and implementation

We provide the UCSC Cell Browser as a free website where scientists can explore a growing collection of single-cell datasets and a freely available python package for scientists to create stable, self-contained visualizations for their own single-cell datasets. Learn more at https://cells.ucsc.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

Tropism of SARS-CoV-2 for Developing Human Cortical Astrocytes

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) readily infects a variety of cell types impacting the function of vital organ systems, with particularly severe impact on respiratory function. It proves fatal for one percent of those infected. Neurological symptoms, which range in severity, accompany a significant proportion of COVID-19 cases, indicating a potential vulnerability of neural cell types. To assess whether human cortical cells can be directly infected by SARS-CoV-2, we utilized primary human cortical tissue and stem cell-derived cortical organoids. We find significant and predominant infection in cortical astrocytes in both primary and organoid cultures, with minimal infection of other cortical populations. Infected astrocytes had a corresponding increase in reactivity characteristics, growth factor signaling, and cellular stress. Although human cortical cells, including astrocytes, have minimal ACE2 expression, we find high levels of alternative coronavirus receptors in infected astrocytes, including DPP4 and CD147. Inhibition of DPP4 reduced infection and decreased expression of the cell stress marker, ARCN1. We find tropism of SARS-CoV-2 for human astrocytes mediated by DPP4, resulting in reactive gliosis-type injury.

Direct analysis of brain phenotypes via neural blastocyst complementation

We provide a protocol for generating forebrain structures in vivo from mouse embryonic stem cells (ESCs) via neural blastocyst complementation (NBC). We developed this protocol for studies of development and function of specific forebrain regions, including the cerebral cortex and hippocampus. We describe a complete workflow, from methods for modifying a given genomic locus in ESCs via CRISPR-Cas9-mediated editing to the generation of mouse chimeras with ESC-reconstituted forebrain regions that can be directly analyzed. The procedure begins with genetic editing of mouse ESCs via CRISPR-Cas9, which can be accomplished in ~4-8 weeks. We provide protocols to achieve fluorescent labeling of ESCs in ~2-3 weeks, which allows tracing of the injected, ESC-derived donor cells in chimeras generated via NBC. Once modified ESCs are ready, NBC chimeras are generated in ~3 weeks via injection of ESCs into genetically programmed blastocysts that are subsequently transferred into pseudo-pregnant fosters. Our in vivo brain organogenesis platform is efficient, allowing functional and systematic analysis of genes and other genomic factors in as little as 3 months, in the context of a whole organism.

Reverse engineering human brain evolution using organoid models

Primate brains vary dramatically in size and organization, but the genetic and developmental basis for these differences has been difficult to study due to lack of experimental models. Pluripotent stem cells and brain organoids provide a potential opportunity for comparative and functional studies of evolutionary differences, particularly during the early stages of neurogenesis. However, many challenges remain, including isolating stem cell lines from additional great ape individuals and species to capture the breadth of ape genetic diversity, improving the reproducibility of organoid models to study evolved differences in cell composition and combining multiple brain regions to capture connectivity relationships. Here, we describe strategies for identifying evolved developmental differences between humans and non-human primates and for isolating the underlying cellular and genetic mechanisms using comparative analyses, chimeric organoid culture, and genome engineering. In particular, we focus on how organoid models could ultimately be applied beyond studies of progenitor cell evolution to decode the origin of recent changes in cellular organization, connectivity patterns, myelination, synaptic development, and physiology that have been implicated in human cognition.

Human-specific tandem repeat expansion and differential gene expression during primate evolution

Short tandem repeats (STRs) and variable number tandem repeats (VNTRs) are important sources of natural and disease-causing variation, yet they have been problematic to resolve in reference genomes and genotype with short-read technology. We created a framework to model the evolution and instability of STRs and VNTRs in apes. We phased and assembled 3 ape genomes (chimpanzee, gorilla, and orangutan) using long-read and 10x Genomics linked-read sequence data for 21,442 human tandem repeats discovered in 6 haplotype-resolved assemblies of Yoruban, Chinese, and Puerto Rican origin. We define a set of 1,584 STRs/VNTRs expanded specifically in humans, including large tandem repeats affecting coding and noncoding portions of genes (e.g., MUC3A, CACNA1C). We show that short interspersed nuclear element-VNTR-Alu (SVA) retrotransposition is the main mechanism for distributing GC-rich human-specific tandem repeat expansions throughout the genome but with a bias against genes. In contrast, we observe that VNTRs not originating from retrotransposons have a propensity to cluster near genes, especially in the subtelomere. Using tissue-specific expression from human and chimpanzee brains, we identify genes where transcript isoform usage differs significantly, likely caused by cryptic splicing variation within VNTRs. Using single-cell expression from cerebral organoids, we observe a strong effect for genes associated with transcription profiles analogous to intermediate progenitor cells. Finally, we compare the sequence composition of some of the largest human-specific repeat expansions and identify 52 STRs/VNTRs with at least 40 uninterrupted pure tracts as candidates for genetically unstable regions associated with disease.

Physiological Models of Human Neuronal Development and Disease

Human neural network development occurs at stages inaccessible to longitudinal monitoring. By transplanting human neurons to the adult mouse brain, recent studies explore human neural circuit formation in realistic cellular and physiological environments, establishing new models for human neurodevelopmental disorders.

Postmitotic Fate Refinement in the Subplate

How is the astonishing diversity of cortical neurons specified? In this issue of Cell Stem Cell, Ozair et al. (2018) leverage hPSC neural differentiation to show that projection neurons undergo prolonged sojourns in the subplate before migrating to deep layers, suggesting that pausing in the subplate may enable integration of intrinsic and extrinsic cues during postmitotic fate refinement.

Paired involvement of human-specific Olduvai domains and NOTCH2NL genes in human brain evolution

Sequences encoding Olduvai (DUF1220) protein domains show the largest human-specific increase in copy number of any coding region in the genome and have been linked to human brain evolution. Most human-specific copies of Olduvai (119/165) are encoded by three NBPF genes that are adjacent to three human-specific NOTCH2NL genes that have been shown to promote cortical neurogenesis. Here, employing genomic, phylogenetic, and transcriptomic evidence, we show that these NOTCH2NL/NBPF gene pairs evolved jointly, as two-gene units, very recently in human evolution, and are likely co-regulated. Remarkably, while three NOTCH2NL paralogs were added, adjacent Olduvai sequences hyper-amplified, adding 119 human-specific copies. The data suggest that human-specific Olduvai domains and adjacent NOTCH2NL genes may function in a coordinated, complementary fashion to promote neurogenesis and human brain expansion in a dosage-related manner.

Establishing Cerebral Organoids as Models of Human-Specific Brain Evolution

Direct comparisons of human and non-human primate brains can reveal molecular pathways underlying remarkable specializations of the human brain. However, chimpanzee tissue is inaccessible during neocortical neurogenesis when differences in brain size first appear. To identify human-specific features of cortical development, we leveraged recent innovations that permit generating pluripotent stem cell-derived cerebral organoids from chimpanzee. Despite metabolic differences, organoid models preserve gene regulatory networks related to primary cell types and developmental processes. We further identified 261 differentially expressed genes in human compared to both chimpanzee organoids and macaque cortex, enriched for recent gene duplications, and including multiple regulators of PI3K-AKT-mTOR signaling. We observed increased activation of this pathway in human radial glia, dependent on two receptors upregulated specifically in human: INSR and ITGB8. Our findings establish a platform for systematic analysis of molecular changes contributing to human brain development and evolution.

Regulation of cell-type-specific transcriptomes by microRNA networks during human brain development

MicroRNAs (miRNAs) regulate many cellular events during brain development by interacting with hundreds of mRNA transcripts. However, miRNAs operate nonuniformly upon the transcriptional profile with an as yet unknown logic. Shortcomings in defining miRNA-mRNA networks include limited knowledge of in vivo miRNA targets and their abundance in single cells. By combining multiple complementary approaches, high-throughput sequencing of RNA isolated by cross-linking immunoprecipitation with an antibody to AGO2 (AGO2-HITS-CLIP), single-cell profiling and computational analyses using bipartite and coexpression networks, we show that miRNA-mRNA interactions operate as functional modules that often correspond to cell-type identities and undergo dynamic transitions during brain development. These networks are highly dynamic during development and over the course of evolution. One such interaction is between radial-glia-enriched ORC4 and miR-2115, a great-ape-specific miRNA, which appears to control radial glia proliferation rates during human brain development.

Identification of cell types in a mouse brain single-cell atlas using low sampling coverage

Background

High throughput methods for profiling the transcriptomes of single cells have recently emerged as transformative approaches for large-scale population surveys of cellular diversity in heterogeneous primary tissues. However, the efficient generation of such atlases will depend on sufficient sampling of diverse cell types while remaining cost-effective to enable a comprehensive examination of organs, developmental stages, and individuals.

Results

To examine the relationship between sampled cell numbers and transcriptional heterogeneity in the context of unbiased cell type classification, we explored the population structure of a publicly available 1.3 million cell dataset from E18.5 mouse brain and validated our findings in published data from adult mice. We propose a computational framework for inferring the saturation point of cluster discovery in a single-cell mRNA-seq experiment, centered around cluster preservation in downsampled datasets. In addition, we introduce a “complexity index,” which characterizes the heterogeneity of cells in a given dataset. Using Cajal-Retzius cells as an example of a limited complexity dataset, we explored whether the detected biological distinctions relate to technical clustering. Surprisingly, we found that clustering distinctions carrying biologically interpretable meaning are achieved with far fewer cells than the originally sampled, though technical saturation of rare populations such as Cajal-Retzius cells is not achieved. We additionally validated these findings with a recently published atlas of cell types across mouse organs and again find using subsampling that a much smaller number of cells recapitulates the cluster distinctions of the complete dataset.

Conclusions

Together, these findings suggest that most of the biologically interpretable cell types from the 1.3 million cell database can be recapitulated by analyzing 50,000 randomly selected cells, indicating that instead of profiling few individuals at high “cellular coverage,” cell atlas studies may instead benefit from profiling more individuals, or many time points at lower cellular coverage and then further enriching for populations of interest. This strategy is ideal for scenarios where cost and time are limited, though extremely rare populations of interest (< 1%) may be identifiable only with much higher cell numbers.

Electronic supplementary material

The online version of this article (10.1186/s12915-018-0580-x) contains supplementary material, which is available to authorized users.

Transcriptional fates of human-specific segmental duplications in brain

Despite the importance of duplicate genes for evolutionary adaptation, accurate gene annotation is often incomplete, incorrect, or lacking in regions of segmental duplication. We developed an approach combining long-read sequencing and hybridization capture to yield full-length transcript information and confidently distinguish between nearly identical genes/paralogs. We used biotinylated probes to enrich for full-length cDNA from duplicated regions, which were then amplified, size-fractionated, and sequenced using single-molecule, long-read sequencing technology, permitting us to distinguish between highly identical genes by virtue of multiple paralogous sequence variants. We examined 19 gene families as expressed in developing and adult human brain, selected for their high sequence identity (average >99%) and overlap with human-specific segmental duplications (SDs). We characterized the transcriptional differences between related paralogs to better understand the birth-death process of duplicate genes and particularly how the process leads to gene innovation. In 48% of the cases, we find that the expressed duplicates have changed substantially from their ancestral models due to novel sites of transcription initiation, splicing, and polyadenylation, as well as fusion transcripts that connect duplication-derived exons with neighboring genes. We detect unannotated open reading frames in genes currently annotated as pseudogenes, while relegating other duplicates to nonfunctional status. Our method significantly improves gene annotation, specifically defining full-length transcripts, isoforms, and open reading frames for new genes in highly identical SDs. The approach will be more broadly applicable to genes in structurally complex regions of other genomes where the duplication process creates novel genes important for adaptive traits.

Human-Specific NOTCH2NL Genes Affect Notch Signaling and Cortical Neurogenesis

Genetic changes causing brain size expansion in human evolution have remained elusive. Notch signaling is essential for radial glia stem cell proliferation and is a determinant of neuronal number in the mammalian cortex. We find that three paralogs of human-specific NOTCH2NL are highly expressed in radial glia. Functional analysis reveals that different alleles of NOTCH2NL have varying potencies to enhance Notch signaling by interacting directly with NOTCH receptors. Consistent with a role in Notch signaling, NOTCH2NL ectopic expression delays differentiation of neuronal progenitors, while deletion accelerates differentiation into cortical neurons. Furthermore, NOTCH2NL genes provide the breakpoints in 1q21.1 distal deletion/duplication syndrome, where duplications are associated with macrocephaly and autism and deletions with microcephaly and schizophrenia. Thus, the emergence of human-specific NOTCH2NL genes may have contributed to the rapid evolution of the larger human neocortex, accompanied by loss of genomic stability at the 1q21.1 locus and resulting recurrent neurodevelopmental disorders.

Spatiotemporal gene expression trajectories reveal developmental hierarchies of the human cortex

Systematic analyses of spatiotemporal gene expression trajectories during organogenesis have been challenging because diverse cell types at different stages of maturation and differentiation coexist in the emerging tissues. We identified discrete cell types as well as temporally and spatially restricted trajectories of radial glia maturation and neurogenesis in developing human telencephalon. These lineage-specific trajectories reveal the expression of neurogenic transcription factors in early radial glia and enriched activation of mammalian target of rapamycin signaling in outer radial glia. Across cortical areas, modest transcriptional differences among radial glia cascade into robust typological distinctions among maturing neurons. Together, our results support a mixed model of topographical, typological, and temporal hierarchies governing cell-type diversity in the developing human telencephalon, including distinct excitatory lineages emerging in rostral and caudal cerebral cortex.

Transformation of the Radial Glia Scaffold Demarcates Two Stages of Human Cerebral Cortex Development

The classic view of cortical development, embodied in the radial unit hypothesis, highlights the ventricular radial glia (vRG) scaffold as a key architectonic feature of the developing neocortex. The scaffold includes continuous fibers spanning the thickness of the developing cortex during neurogenesis across mammals. However, we find that in humans, the scaffold transforms into a physically discontinuous structure during the transition from infragranular to supragranular neuron production. As a consequence of this transformation, supragranular layer neurons arrive at their terminal positions in the cortical plate along outer radial glia (oRG) cell fibers. In parallel, the radial glia that contact the ventricle develop distinct gene expression profile and "truncated" morphology. We propose a supragranular layer expansion hypothesis that posits a deterministic role of oRG cells in the radial and tangential expansion of supragranular layers in primates, with implications for patterns of neuronal migration, area patterning, and cortical folding.

Human iPSC-Derived Cerebral Organoids Model Cellular Features of Lissencephaly and Reveal Prolonged Mitosis of Outer Radial Glia

Classical lissencephaly is a genetic neurological disorder associated with mental retardation and intractable epilepsy, and Miller-Dieker syndrome (MDS) is the most severe form of the disease. In this study, to investigate the effects of MDS on human progenitor subtypes that control neuronal output and influence brain topology, we analyzed cerebral organoids derived from control and MDS-induced pluripotent stem cells (iPSCs) using time-lapse imaging, immunostaining, and single-cell RNA sequencing. We saw a cell migration defect that was rescued when we corrected the MDS causative chromosomal deletion and severe apoptosis of the founder neuroepithelial stem cells, accompanied by increased horizontal cell divisions. We also identified a mitotic defect in outer radial glia, a progenitor subtype that is largely absent from lissencephalic rodents but critical for human neocortical expansion. Our study, therefore, deepens our understanding of MDS cellular pathogenesis and highlights the broad utility of cerebral organoids for modeling human neurodevelopmental disorders.

Single-cell sequencing maps gene expression to mutational phylogenies in PDGF- and EGF-driven gliomas

Glioblastoma multiforme (GBM) is the most common and aggressive type of primary brain tumor. Epidermal growth factor (EGF) and platelet-derived growth factor (PDGF) receptors are frequently amplified and/or possess gain-of-function mutations in GBM However, clinical trials of tyrosine-kinase inhibitors have shown disappointing efficacy, in part due to intra-tumor heterogeneity. To assess the effect of clonal heterogeneity on gene expression, we derived an approach to map single-cell expression profiles to sequentially acquired mutations identified from exome sequencing. Using 288 single cells, we constructed high-resolution phylogenies of EGF-driven and PDGF-driven GBMs, modeling transcriptional kinetics during tumor evolution. Descending the phylogenetic tree of a PDGF-driven tumor corresponded to a progressive induction of an oligodendrocyte progenitor-like cell type, expressing pro-angiogenic factors. In contrast, phylogenetic analysis of an EGFR-amplified tumor showed an up-regulation of pro-invasive genes. An in-frame deletion in a specific dimerization domain of PDGF receptor correlates with an up-regulation of growth pathways in a proneural GBM and enhances proliferation when ectopically expressed in glioma cell lines. In-frame deletions in this domain are frequent in public GBM data.

Primate Neurons Flex Their Musclin

Sensory experience evokes long-lasting changes in neural circuits through activity-dependent gene expression. Ataman et al. (2016) report in Nature that primates evolved novel transcriptional responses to neuronal activity, including induction of musclin/osteocrin (OSTN), which may regulate specialized aspects of primate neural circuits.

Expression Analysis Highlights AXL as a Candidate Zika Virus Entry Receptor in Neural Stem Cells

The recent outbreak of Zika virus (ZIKV) in Brazil has been linked to substantial increases in fetal abnormalities and microcephaly. However, information about the underlying molecular and cellular mechanisms connecting viral infection to these defects remains limited. In this study we have examined the expression of receptors implicated in cell entry of several enveloped viruses including ZIKV across diverse cell types in the developing brain. Using single-cell RNA-seq and immunohistochemistry, we found that the candidate viral entry receptor AXL is highly expressed by human radial glial cells, astrocytes, endothelial cells, and microglia in developing human cortex and by progenitor cells in developing retina. We also show that AXL expression in radial glia is conserved in developing mouse and ferret cortex and in human stem cell-derived cerebral organoids, highlighting multiple experimental systems that could be applied to study mechanisms of ZIKV infectivity and effects on brain development.

Single-cell analysis of long non-coding RNAs in the developing human neocortex

Background

Long non-coding RNAs (lncRNAs) comprise a diverse class of transcripts that can regulate molecular and cellular processes in brain development and disease. LncRNAs exhibit cell type- and tissue-specific expression, but little is known about the expression and function of lncRNAs in the developing human brain. Furthermore, it has been unclear whether lncRNAs are highly expressed in subsets of cells within tissues, despite appearing lowly expressed in bulk populations.

Results

We use strand-specific RNA-seq to deeply profile lncRNAs from polyadenylated and total RNA obtained from human neocortex at different stages of development, and we apply this reference to analyze the transcriptomes of single cells. While lncRNAs are generally detected at low levels in bulk tissues, single-cell transcriptomics of hundreds of neocortex cells reveal that many lncRNAs are abundantly expressed in individual cells and are cell type-specific. Notably, LOC646329 is a lncRNA enriched in single radial glia cells but is detected at low abundance in tissues. CRISPRi knockdown of LOC646329 indicates that this lncRNA regulates cell proliferation.

Conclusion

The discrete and abundant expression of lncRNAs among individual cells has important implications for both their biological function and utility for distinguishing neural cell types.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-016-0932-1) contains supplementary material, which is available to authorized users.

Molecular identity of human outer radial glia during cortical development

Radial glia, the neural stem cells of the neocortex, are located in two niches: the ventricular zone and outer subventricular zone. Although outer subventricular zone radial glia may generate the majority of human cortical neurons, their molecular features remain elusive. By analyzing gene expression across single cells, we find that outer radial glia preferentially express genes related to extracellular matrix formation, migration, and stemness, including TNC, PTPRZ1, FAM107A, HOPX, and LIFR. Using dynamic imaging, immunostaining, and clonal analysis, we relate these molecular features to distinctive behaviors of outer radial glia, demonstrate the necessity of STAT3 signaling for their cell cycle progression, and establish their extensive proliferative potential. These results suggest that outer radial glia directly support the subventricular niche through local production of growth factors, potentiation of growth factor signals by extracellular matrix proteins, and activation of self-renewal pathways, thereby enabling the developmental and evolutionary expansion of the human neocortex.

Radial glia require PDGFD-PDGFRβ signalling in human but not mouse neocortex

Evolutionary expansion of the human neocortex underlies many of our unique mental abilities. This expansion has been attributed to the increased proliferative potential of radial glia (RG; neural stem cells) and their subventricular dispersion from the periventricular niche during neocortical development. Such adaptations may have evolved through gene expression changes in RG. However, whether or how RG gene expression varies between humans and other species is unknown. Here we show that the transcriptional profiles of human and mouse neocortical RG are broadly conserved during neurogenesis, yet diverge for specific signalling pathways. By analysing differential gene co-expression relationships between the species, we demonstrate that the growth factor PDGFD is specifically expressed by RG in human, but not mouse, corticogenesis. We also show that the expression domain of PDGFRβ, the cognate receptor for PDGFD, is evolutionarily divergent, with high expression in the germinal region of dorsal human neocortex but not in the mouse. Pharmacological inhibition of PDGFD-PDGFRβ signalling in slice culture prevents normal cell cycle progression of neocortical RG in human, but not mouse. Conversely, injection of recombinant PDGFD or ectopic expression of constitutively active PDGFRβ in developing mouse neocortex increases the proportion of RG and their subventricular dispersion. These findings highlight the requirement of PDGFD-PDGFRβ signalling for human neocortical development and suggest that local production of growth factors by RG supports the expanded germinal region and progenitor heterogeneity of species with large brains.

Low-coverage single-cell mRNA sequencing reveals cellular heterogeneity and activated signaling pathways in developing cerebral cortex

Large-scale surveys of single-cell gene expression have the potential to reveal rare cell populations and lineage relationships, but require efficient methods for cell capture and mRNA sequencing^1–4. Although cellular barcoding strategies allow parallel sequencing of single cells at ultra-low depths⁵, the limitations of shallow sequencing have not been directly investigated. By capturing 301 single cells from 11 populations using microfluidics and analyzing single-cell transcriptomes across downsampled sequencing depths, we demonstrate that shallow single-cell mRNA sequencing (~50,000 reads per cell) is sufficient for unbiased cell-type classification and biomarker identification. In developing cortex we identify diverse cell types including multiple progenitor and neuronal subtypes, and we identify EGR1 and FOS as previously unreported candidate targets of Notch signaling in human but not mouse radial glia. Our strategy establishes an efficient method for unbiased analysis and comparison of cell populations from heterogeneous tissue by microfluidic single-cell capture and low-coverage sequencing of many cells.

cis-Regulatory changes in Kit ligand expression and parallel evolution of pigmentation in sticklebacks and humans

Dramatic pigmentation changes have evolved within most vertebrate groups, including fish and humans. Here we use genetic crosses in sticklebacks to investigate the parallel origin of pigmentation changes in natural populations. High-resolution mapping and expression experiments show that light gills and light ventrums map to a divergent regulatory allele of the Kit ligand (Kitlg) gene. The divergent allele reduces expression in gill and skin tissue and is shared by multiple derived freshwater populations with reduced pigmentation. In humans, Europeans and East Asians also share derived alleles at the KITLG locus. Strong signatures of selection map to regulatory regions surrounding the gene, and admixture mapping shows that the KITLG genomic region has a significant effect on human skin color. These experiments suggest that regulatory changes in Kitlg contribute to natural variation in vertebrate pigmentation, and that similar genetic mechanisms may underlie rapid evolutionary change in fish and humans.