A single-cell time-lapse of mouse prenatal development from gastrula to birth.

The house mouse (Mus musculus) is an exceptional model system, combining genetic tractability with close evolutionary affinity to humans1,2. Mouse gestation lasts only 3 weeks, during which the genome orchestrates the astonishing transformation of a single-cell zygote into a free-living pup composed of more than 500 million cells. Here, to establish a global framework for exploring mammalian development, we applied optimized single-cell combinatorial indexing3 to profile the transcriptional states of 12.4 million nuclei from 83 embryos, precisely staged at 2- to 6-hour intervals spanning late gastrulation (embryonic day 8) to birth (postnatal day 0). From these data, we annotate hundreds of cell types and explore the ontogenesis of the posterior embryo during somitogenesis and of kidney, mesenchyme, retina and early neurons. We leverage the temporal resolution and sampling depth of these whole-embryo snapshots, together with published data4-8 from earlier timepoints, to construct a rooted tree of cell-type relationships that spans the entirety of prenatal development, from zygote to birth. Throughout this tree, we systematically nominate genes encoding transcription factors and other proteins as candidate drivers of the in vivo differentiation of hundreds of cell types. Remarkably, the most marked temporal shifts in cell states are observed within one hour of birth and presumably underlie the massive physiological adaptations that must accompany the successful transition of a mammalian fetus to life outside the womb.

A single-cell multi-omic atlas spanning the adult rhesus macaque brain.

Cataloging the diverse cellular architecture of the primate brain is crucial for understanding cognition, behavior, and disease in humans. Here, we generated a brain-wide single-cell multimodal molecular atlas of the rhesus macaque brain. Together, we profiled 2.58 M transcriptomes and 1.59 M epigenomes from single nuclei sampled from 30 regions across the adult brain. Cell composition differed extensively across the brain, revealing cellular signatures of region-specific functions. We also identified 1.19 M candidate regulatory elements, many previously unidentified, allowing us to explore the landscape of cis-regulatory grammar and neurological disease risk in a cell type-specific manner. Altogether, this multi-omic atlas provides an open resource for investigating the evolution of the human brain and identifying novel targets for disease interventions.

Cross-species imputation and comparison of single-cell transcriptomic profiles.

Cross-species comparison and prediction of gene expression profiles are important to understand regulatory changes during evolution and to transfer knowledge learned from model organisms to humans. Single-cell RNA-seq (scRNA-seq) profiles enable us to capture gene expression profiles with respect to variations among individual cells; however, cross-species comparison of scRNA-seq profiles is challenging because of data sparsity, batch effects, and the lack of one-to-one cell matching across species. Moreover, single-cell profiles are challenging to obtain in certain biological contexts, limiting the scope of hypothesis generation. Here we developed Icebear, a neural network framework that decomposes single-cell measurements into factors representing cell identity, species, and batch factors. Icebear enables accurate prediction of single-cell gene expression profiles across species, thereby providing high-resolution cell type and disease profiles in under-characterized contexts. Icebear also facilitates direct cross-species comparison of single-cell expression profiles for conserved genes that are located on the X chromosome in eutherian mammals but on autosomes in chicken. This comparison, for the first time, revealed evolutionary and diverse adaptations of X-chromosome upregulation in mammals.

A single-cell transcriptional timelapse of mouse embryonic development, from gastrula to pup.

The house mouse, Mus musculus, is an exceptional model system, combining genetic tractability with close homology to human biology. Gestation in mouse development lasts just under three weeks, a period during which its genome orchestrates the astonishing transformation of a single cell zygote into a free-living pup composed of >500 million cells. Towards a global framework for exploring mammalian development, we applied single cell combinatorial indexing (sci-*) to profile the transcriptional states of 12.4 million nuclei from 83 precisely staged embryos spanning late gastrulation (embryonic day 8 or E8) to birth (postnatal day 0 or P0), with 2-hr temporal resolution during somitogenesis, 6-hr resolution through to birth, and 20-min resolution during the immediate postpartum period. From these data (E8 to P0), we annotate dozens of trajectories and hundreds of cell types and perform deeper analyses of the unfolding of the posterior embryo during somitogenesis as well as the ontogenesis of the kidney, mesenchyme, retina, and early neurons. Finally, we leverage the depth and temporal resolution of these whole embryo snapshots, together with other published data, to construct and curate a rooted tree of cell type relationships that spans mouse development from zygote to pup. Throughout this tree, we systematically nominate sets of transcription factors (TFs) and other genes as candidate drivers of the in vivo differentiation of hundreds of mammalian cell types. Remarkably, the most dramatic shifts in transcriptional state are observed in a restricted set of cell types in the hours immediately following birth, and presumably underlie the massive changes in physiology that must accompany the successful transition of a placental mammal to extrauterine life.

A human cell atlas of fetal chromatin accessibility.

The chromatin landscape underlying the specification of human cell types is of fundamental interest. We generated human cell atlases of chromatin accessibility and gene expression in fetal tissues. For chromatin accessibility, we devised a three-level combinatorial indexing assay and applied it to 53 samples representing 15 organs, profiling ~800,000 single cells. We leveraged cell types defined by gene expression to annotate these data and cataloged hundreds of thousands of candidate regulatory elements that exhibit cell type-specific chromatin accessibility. We investigated the properties of lineage-specific transcription factors (such as POU2F1 in neurons), organ-specific specializations of broadly distributed cell types (such as blood and endothelial), and cell type-specific enrichments of complex trait heritability. These data represent a rich resource for the exploration of in vivo human gene regulation in diverse tissues and cell types.

A human cell atlas of fetal gene expression.

The gene expression program underlying the specification of human cell types is of fundamental interest. We generated human cell atlases of gene expression and chromatin accessibility in fetal tissues. For gene expression, we applied three-level combinatorial indexing to >110 samples representing 15 organs, ultimately profiling ~4 million single cells. We leveraged the literature and other atlases to identify and annotate hundreds of cell types and subtypes, both within and across tissues. Our analyses focused on organ-specific specializations of broadly distributed cell types (such as blood, endothelial, and epithelial), sites of fetal erythropoiesis (which notably included the adrenal gland), and integration with mouse developmental atlases (such as conserved specification of blood cells). These data represent a rich resource for the exploration of in vivo human gene expression in diverse tissues and cell types.

Whole-Genome and RNA Sequencing Reveal Variation and Transcriptomic Coordination in the Developing Human Prefrontal Cortex.

Gene expression levels vary across developmental stage, cell type, and region in the brain. Genomic variants also contribute to the variation in expression, and some neuropsychiatric disorder loci may exert their effects through this mechanism. To investigate these relationships, we present BrainVar, a unique resource of paired whole-genome and bulk tissue RNA sequencing from the dorsolateral prefrontal cortex of 176 individuals across prenatal and postnatal development. Here we identify common variants that alter gene expression (expression quantitative trait loci [eQTLs]) constantly across development or predominantly during prenatal or postnatal stages. Both "constant" and "temporal-predominant" eQTLs are enriched for loci associated with neuropsychiatric traits and disorders and colocalize with specific variants. Expression levels of more than 12,000 genes rise or fall in a concerted late-fetal transition, with the transitional genes enriched for cell-type-specific genes and neuropsychiatric risk loci, underscoring the importance of cataloging developmental trajectories in understanding cortical physiology and pathology.

Marijuana use differentially affects cannabinoid receptor expression in early gestational human endometrium and placenta.

Marijuana is one of the most abused drugs among pregnant women leading to maternal and fetal abnormalities. Cannabinoids are the active ingredients of marijuana, which interact with cannabinoid receptors such as CNR1 and CNR2 to activate cellular signaling pathways. Human endometrium and placenta are known to express CNR1 and CNR2 and can respond to cannabinoid signaling. In this study, we show that marijuana use significantly increases mRNA or protein expression of CNR1 and CNR2 in human endometrium from the first and early second trimester pregnancies, with minor effects on placental expression of CNRs.

Keratinocytes produce IL-17c to protect peripheral nervous systems during human HSV-2 reactivation.

Despite frequent herpes simplex virus (HSV) reactivation, peripheral nerve destruction and sensory anesthesia are rare. We discovered that skin biopsies obtained during asymptomatic human HSV-2 reactivation exhibit a higher density of nerve fibers relative to biopsies during virological and clinical quiescence. We evaluated the effects of HSV infection on keratinocytes, the initial target of HSV replication, to better understand this observation. Keratinocytes produced IL-17c during HSV-2 reactivation, and IL-17RE, an IL-17c-specific receptor, was expressed on nerve fibers in human skin and sensory neurons in dorsal root ganglia. In ex vivo experiments, exogenous human IL-17c provided directional guidance and promoted neurite growth and branching in microfluidic devices. Exogenous murine IL-17c pretreatment reduced apoptosis in HSV-2-infected primary neurons. These results suggest that IL-17c is a neurotrophic cytokine that protects peripheral nerve systems during HSV reactivation. This mechanism could explain the lack of nerve damage from recurrent HSV infection and may provide insight to understanding and treating sensory peripheral neuropathies.

Abnormal glycosylation in Joubert syndrome type 10.

BACKGROUND: The discovery of disease pathogenesis requires systematic agnostic screening of multiple homeostatic processes that may become deregulated. We illustrate this principle in the evaluation and diagnosis of a 5-year-old boy with Joubert syndrome type 10 (JBTS10). He carried the OFD1 mutation p.Gln886Lysfs*2 (NM_003611.2: c.2656del) and manifested features of Joubert syndrome. METHODS: We integrated exome sequencing, MALDI-TOF mass spectrometry analyses of plasma and cultured dermal fibroblasts glycomes, and full clinical evaluation of the proband. Analyses of cilia formation and lectin staining were performed by immunofluorescence. Measurement of cellular nucleotide sugar levels was performed with high-performance anion-exchange chromatography with pulsed amperometric detection. Statistical analyses utilized the Student's and Fisher's exact t tests. RESULTS: Glycome analyses of plasma and cultured dermal fibroblasts identified abnormal N- and O-linked glycosylation profiles. These findings replicated in two unrelated males with OFD1 mutations. Cultured fibroblasts from affected individuals had a defect in ciliogenesis. The proband's fibroblasts also had an abnormally elevated nuclear sialylation signature and increased total cellular levels of CMP-sialic acid. Ciliogenesis and each glycosylation anomaly were rescued by expression of wild-type OFD1. CONCLUSIONS: The rescue of ciliogenesis and glycosylation upon reintroduction of WT OFD1 suggests that both contribute to the pathogenesis of JBTS10.

Mutations in LAMA1 cause cerebellar dysplasia and cysts with and without retinal dystrophy.

Cerebellar dysplasia with cysts (CDC) is an imaging finding typically seen in combination with cobblestone cortex and congenital muscular dystrophy in individuals with dystroglycanopathies. More recently, CDC was reported in seven children without neuromuscular involvement (Poretti-Boltshauser syndrome). Using a combination of homozygosity mapping and whole-exome sequencing, we identified biallelic mutations in LAMA1 as the cause of CDC in seven affected individuals (from five families) independent from those included in the phenotypic description of Poretti-Boltshauser syndrome. Most of these individuals also have high myopia, and some have retinal dystrophy and patchy increased T2-weighted fluid-attenuated inversion recovery (T2/FLAIR) signal in cortical white matter. In one additional family, we identified two siblings who have truncating LAMA1 mutations in combination with retinal dystrophy and mild cerebellar dysplasia without cysts, indicating that cysts are not an obligate feature associated with loss of LAMA1 function. This work expands the phenotypic spectrum associated with the lamininopathy disorders and highlights the tissue-specific roles played by different laminin-encoding genes.

Hypomorphism for RPGRIP1L, a ciliary gene vicinal to the FTO locus, causes increased adiposity in mice.

Common polymorphisms in the first intron of FTO are associated with increased body weight in adults. Previous studies have suggested that a CUX1-regulatory element within the implicated FTO region controls expression of FTO and the nearby ciliary gene, RPGRIP1L. Given the role of ciliary genes in energy homeostasis, we hypothesized that mice hypomorphic for Rpgrip1l would display increased adiposity. We find that Rpgrip1l⁺/⁻ mice are hyperphagic and fatter, and display diminished suppression of food intake in response to leptin administration. In the hypothalamus of Rpgrip1l⁺/⁻ mice, and in human fibroblasts with hypomorphic mutations in RPGRIP1L, the number of AcIII-positive cilia is diminished, accompanied by impaired convening of the leptin receptor to the vicinity of the cilium, and diminished pStat3 in response to leptin. These findings suggest that RPGRIP1L may be partly or exclusively responsible for the obesity susceptibility signal at the FTO locus.

Mutations in CSPP1 cause primary cilia abnormalities and Joubert syndrome with or without Jeune asphyxiating thoracic dystrophy.

Joubert syndrome (JBTS) is a recessive ciliopathy in which a subset of affected individuals also have the skeletal dysplasia Jeune asphyxiating thoracic dystrophy (JATD). Here, we have identified biallelic truncating CSPP1 (centrosome and spindle pole associated protein 1) mutations in 19 JBTS-affected individuals, four of whom also have features of JATD. CSPP1 mutations explain ∼5% of JBTS in our cohort, and despite truncating mutations in all affected individuals, the range of phenotypic severity is broad. Morpholino knockdown of cspp1 in zebrafish caused phenotypes reported in other zebrafish models of JBTS (curved body shape, pronephric cysts, and cerebellar abnormalities) and reduced ciliary localization of Arl13b, further supporting loss of CSPP1 function as a cause of JBTS. Fibroblasts from affected individuals with CSPP1 mutations showed reduced numbers of primary cilia and/or short primary cilia, as well as reduced axonemal localization of ciliary proteins ARL13B and adenylyl cyclase III. In summary, CSPP1 mutations are a major cause of the Joubert-Jeune phenotype in humans; however, the mechanism by which these mutations lead to both JBTS and JATD remains unknown.

Multiplex targeted sequencing identifies recurrently mutated genes in autism spectrum disorders.

Exome sequencing studies of autism spectrum disorders (ASDs) have identified many de novo mutations but few recurrently disrupted genes. We therefore developed a modified molecular inversion probe method enabling ultra-low-cost candidate gene resequencing in very large cohorts. To demonstrate the power of this approach, we captured and sequenced 44 candidate genes in 2446 ASD probands. We discovered 27 de novo events in 16 genes, 59% of which are predicted to truncate proteins or disrupt splicing. We estimate that recurrent disruptive mutations in six genes-CHD8, DYRK1A, GRIN2B, TBR1, PTEN, and TBL1XR1-may contribute to 1% of sporadic ASDs. Our data support associations between specific genes and reciprocal subphenotypes (CHD8-macrocephaly and DYRK1A-microcephaly) and replicate the importance of a ß-catenin-chromatin-remodeling network to ASD etiology.