Messenger RNA transport on lysosomal vesicles maintains axonal mitochondrial homeostasis and prevents axonal degeneration.

In neurons, RNA granules are transported along the axon for local translation away from the soma. Recent studies indicate that some of this transport involves hitchhiking of RNA granules on lysosome-related vesicles. In the present study, we leveraged the ability to prevent transport of these vesicles into the axon by knockout of the lysosome-kinesin adaptor BLOC-one-related complex (BORC) to identify a subset of axonal mRNAs that depend on lysosome-related vesicles for transport. We found that BORC knockout causes depletion of a large group of axonal mRNAs mainly encoding ribosomal and mitochondrial/oxidative phosphorylation proteins. This depletion results in mitochondrial defects and eventually leads to axonal degeneration in human induced pluripotent stem cell (iPSC)-derived and mouse neurons. Pathway analyses of the depleted mRNAs revealed a mechanistic connection of BORC deficiency with common neurodegenerative disorders. These results demonstrate that mRNA transport on lysosome-related vesicles is critical for the maintenance of axonal homeostasis and that its failure causes axonal degeneration.

Loss of Ezh2 in the medial ganglionic eminence alters interneuron fate, cell morphology and gene expression profiles.

Introduction: Enhancer of zeste homolog 2 (Ezh2) is responsible for trimethylation of histone 3 at lysine 27 (H3K27me3), resulting in repression of gene expression. Here, we explore the role of Ezh2 in forebrain GABAergic interneuron development. Methods: We removed Ezh2 in the MGE by generating Nkx2-1Cre;Ezh2 conditional knockout mice. We then characterized changes in MGE-derived interneuron fate and electrophysiological properties in juvenile mice, as well as alterations in gene expression, chromatin accessibility and histone modifications in the MGE. Results: Loss of Ezh2 increases somatostatin-expressing (SST+) and decreases parvalbumin-expressing (PV+) interneurons in the forebrain. We observe fewer MGE-derived interneurons in the first postnatal week, indicating reduced interneuron production. Intrinsic electrophysiological properties in SST+ and PV+ interneurons are normal, but PV+ interneurons display increased axonal complexity in Ezh2 mutant mice. Single nuclei multiome analysis revealed differential gene expression patterns in the embryonic MGE that are predictive of these cell fate changes. Lastly, CUT&Tag analysis revealed that some genomic loci are particularly resistant or susceptible to shifts in H3K27me3 levels in the absence of Ezh2, indicating differential selectivity to epigenetic perturbation. Discussion: Thus, loss of Ezh2 in the MGE alters interneuron fate, morphology, and gene expression and regulation. These findings have important implications for both normal development and potentially in disease etiologies.

Multimodal spatiotemporal transcriptomic resolution of embryonic palate osteogenesis.

The terminal differentiation of osteoblasts and subsequent formation of bone marks an important phase in palate development that leads to the separation of the oral and nasal cavities. While the morphogenetic events preceding palatal osteogenesis are well explored, major gaps remain in our understanding of the molecular mechanisms driving the formation of this bony union of the fusing palate. Through bulk, single-nucleus, and spatially resolved RNA-sequencing analyses of the developing secondary palate, we identify a shift in transcriptional programming between embryonic days 14.5 and 15.5 pinpointing the onset of osteogenesis. We define spatially restricted expression patterns of key osteogenic marker genes that are differentially expressed between these developmental timepoints. Finally, we identify genes in the palate highly expressed by palate nasal epithelial cells, also enriched within palatal osteogenic mesenchymal cells. This investigation provides a relevant framework to advance palate-specific diagnostic and therapeutic biomarker discovery.

Extensive diversity in RNA termination and regulation revealed by transcriptome mapping for the Lyme pathogen Borrelia burgdorferi.

Transcription termination is an essential and dynamic process that can tune gene expression in response to diverse molecular signals. Yet, the genomic positions, molecular mechanisms, and regulatory consequences of termination have only been studied thoroughly in model bacteria. Here, we use several RNA-seq approaches to map RNA ends for the transcriptome of the spirochete Borrelia burgdorferi - the etiological agent of Lyme disease. We identify complex gene arrangements and operons, untranslated regions and small RNAs. We predict intrinsic terminators and experimentally test examples of Rho-dependent transcription termination. Remarkably, 63% of RNA 3' ends map upstream of or internal to open reading frames (ORFs), including genes involved in the unique infectious cycle of B. burgdorferi. We suggest these RNAs result from premature termination, processing and regulatory events such as cis-acting regulation. Furthermore, the polyamine spermidine globally influences the generation of truncated mRNAs. Collectively, our findings provide insights into transcription termination and uncover an abundance of potential RNA regulators in B. burgdorferi.

Germline loss-of-function PAM variants are enriched in subjects with pituitary hypersecretion.

Introduction: Pituitary adenomas (PAs) are common, usually benign tumors of the anterior pituitary gland which, for the most part, have no known genetic cause. PAs are associated with major clinical effects due to hormonal dysregulation and tumoral impingement on vital brain structures. PAM encodes a multifunctional protein responsible for the essential C-terminal amidation of secreted peptides. Methods: Following the identification of a loss-of-function variant (p.Arg703Gln) in the peptidylglycine a-amidating monooxygenase (PAM) gene in a family with pituitary gigantism, we investigated 299 individuals with sporadic PAs and 17 familial isolated PA kindreds for PAM variants. Genetic screening was performed by germline and tumor sequencing and germline copy number variation (CNV) analysis. Results: In germline DNA, we detected seven heterozygous, likely pathogenic missense, truncating, and regulatory SNVs. These SNVs were found in sporadic subjects with growth hormone excess (p.Gly552Arg and p.Phe759Ser), pediatric Cushing disease (c.-133T>C and p.His778fs), or different types of PAs (c.-361G>A, p.Ser539Trp, and p.Asp563Gly). The SNVs were functionally tested in vitro for protein expression and trafficking by Western blotting, splicing by minigene assays, and amidation activity in cell lysates and serum samples. These analyses confirmed a deleterious effect on protein expression and/or function. By interrogating 200,000 exomes from the UK Biobank, we confirmed a significant association of the PAM gene and rare PAM SNVs with diagnoses linked to pituitary gland hyperfunction. Conclusion: The identification of PAM as a candidate gene associated with pituitary hypersecretion opens the possibility of developing novel therapeutics based on altering PAM function.

Integrated spatiotemporal transcriptomic resolution of embryonic palate osteogenesis.

The differentiation of osteoblasts and the subsequent formation of bone marks an important terminal phase in palate formation that leads to the separation of the oral and nasal cavities. While the developmental events that precede palatal osteogenesis are well explored, major gaps remain in our understanding of the molecular mechanisms that lead to the bony union of fusing palatal shelves. Herein, the timeline of osteogenic transcriptional programming is unveiled in the embryonic palate by way of integrated bulk, single-cell, and spatially resolved RNA-seq analyses. We define spatially restricted expression patterns of key marker genes, both regulatory and structural, that are differentially expressed during palatal fusion, including the identification of several novel genes ( Deup1, Dynlrb2, Lrrc23 ) spatially restricted in expression to the palate, providing a relevant framework for future studies that identify new candidate genes for cleft palate anomalies in humans as well as the timing of mammalian embryonic palatal osteogenesis.

Cerebrospinal Fluid Protein Biomarker Discovery in CLN3.

Syndromic CLN3-Batten is a fatal, pediatric, neurodegenerative disease caused by variants in CLN3, which encodes the endolysosomal transmembrane CLN3 protein. No approved treatment for CLN3 is currently available. The protracted and asynchronous disease presentation complicates the evaluation of potential therapies using clinical disease progression parameters. Biomarkers as surrogates to measure the progression and effect of potential therapeutics are needed. We performed proteomic discovery studies using cerebrospinal fluid (CSF) samples from 28 CLN3-affected and 32 age-similar non-CLN3 individuals. Proximal extension assay (PEA) of 1467 proteins and untargeted data-dependent mass spectrometry [MS; MassIVE FTP server (ftp://MSV000090147@massive.ucsd.edu)] were used to generate orthogonal lists of protein marker candidates. At an adjusted p-value of <0.1 and threshold CLN3/non-CLN3 fold-change ratio of 1.5, PEA identified 54 and MS identified 233 candidate biomarkers. Some of these (NEFL, CHIT1) have been previously linked with other neurologic conditions. Others (CLPS, FAM217B, QRICH2, KRT16, ZNF333) appear to be novel. Both methods identified 25 candidate biomarkers, including CHIT1, NELL1, and ISLR2 which had absolute fold-change ratios >2. NELL1 and ISLR2 regulate axonal development in neurons and are intriguing new candidates for further investigation in CLN3. In addition to identifying candidate proteins for CLN3 research, this study provides a comparison of two large-scale proteomic discovery methods in CSF.

Absence of TRIC-B from type XIV Osteogenesis Imperfecta osteoblasts alters cell adhesion and mitochondrial function - A multi-omics study.

Osteogenesis Imperfecta (OI) is a heritable collagen-related bone dysplasia characterized by bone fractures, growth deficiency and skeletal deformity. Type XIV OI is a recessive OI form caused by null mutations in TMEM38B, which encodes the ER membrane intracellular cation channel TRIC-B. Previously, we showed that absence of TMEM38B alters calcium flux in the ER of OI patient osteoblasts and fibroblasts, which further disrupts collagen synthesis and secretion. How the absence of TMEM38B affects osteoblast function is still poorly understood. Here we further investigated the role of TMEM38B in human osteoblast differentiation and mineralization. TMEM38B-null osteoblasts showed altered expression of osteoblast marker genes and decreased mineralization. RNA-Seq analysis revealed that cell-cell adhesion was one of the most downregulated pathways in TMEM38B-null osteoblasts, with further validation by real-time PCR and Western blot. Gap and tight junction proteins were also decreased by TRIC-B absence, both in patient osteoblasts and in calvarial osteoblasts of Tmem38b-null mice. Disrupted cell adhesion decreased mutant cell proliferation and cell cycle progression. An important novel finding was that TMEM38B-null osteoblasts had elongated mitochondria with altered fusion and fission markers, MFN2 and DRP1. In addition, TMEM38B-null osteoblasts exhibited a significant increase in superoxide production in mitochondria, further supporting mitochondrial dysfunction. Together these results emphasize the novel role of TMEM38B/TRIC-B in osteoblast differentiation, affecting cell-cell adhesion processes, gap and tight junction, proliferation, cell cycle, and mitochondrial function.

Identification of cerebral spinal fluid protein biomarkers in Niemann-Pick disease, type C1.

BACKGROUND: Niemann-Pick disease, type C1 (NPC1) is an ultrarare, recessive, lethal, lysosomal disease characterized by progressive cerebellar ataxia and cognitive impairment. Although the NPC1 phenotype is heterogeneous with variable age of onset, classical NPC1 is a pediatric disorder. Currently there are no therapies approved by the FDA and therapeutics trials for NPC1 are complicated by disease rarity, heterogeneity, and the relatively slow rate of neurological decline. Thus, identification of disease relevant biomarkers is necessary to provide tools that can support drug development efforts for this devastating neurological disease. METHODS: Proximal extension assays (O-link® Explore 1536) were used to compare cerebrospinal fluid (CSF) samples from individuals with NPC1 enrolled in a natural history study and non-NPC1 comparison samples. Relative expression levels of 1467 proteins were determined, and candidate protein biomarkers were identified by evaluating fold-change and adjusted Kruskal-Wallis test p-values. Selected proteins were orthogonally confirmed using ELISA. To gain insight into disease progression and severity we evaluated the altered protein expression with respect to clinically relevant phenotypic aspects: NPC Neurological Severity Score (NPC1 NSS), Annual Severity Increment Score (ASIS) and age of neurological onset. RESULTS: This study identified multiple proteins with altered levels in CSF from individuals with NPC1 compared to non-NPC1 samples. These included proteins previously shown to be elevated in NPC1 (NEFL, MAPT, CHIT1, CALB1) and additional proteins confirmed by orthogonal assays (PARK7, CALB2/calretinin, CHI3L1/YKL-40, MIF, CCL18 and ENO2). Correlations with clinically relevant phenotypic parameters demonstrated moderate negative (p = 0.0210, r = -0.41) and possible moderate positive (p = 0.0631, r = 0.33) correlation of CSF CALB2 levels with age of neurological onset and ASIS, respectively. CSF CHI3L1 levels showed a moderate positive (p = 0.0183, r = 0.40) correlation with the concurrent NPC1 NSS. A strong negative correlation (p = 0.0016, r = -0.648) was observed between CSF CCL18 and age of neurological onset for childhood/adolescent cases. CSF CCL18 levels also showed a strong positive correlation (p = 0.0017, r = 0.61) with ASIS. CONCLUSION: Our study identified and validated multiple proteins in CSF from individuals with NPC1 that are candidates for further investigation in a larger cohort. These analytes may prove to be useful as supportive data in therapeutic trials. TRIAL REGISTRATIONS: NCT00344331, NCT00001721, NCT02931682.

Germline loss-of-function PAM variants are enriched in subjects with pituitary hypersecretion.

Pituitary adenomas (PAs) are common, usually benign tumors of the anterior pituitary gland which, for the most part, have no known genetic cause. PAs are associated with major clinical effects due to hormonal dysregulation and tumoral impingement on vital brain structures. Following the identification of a loss-of-function variant (p.Arg703Gln) in the PAM gene in a family with pituitary gigantism, we investigated 299 individuals with sporadic PAs and 17 familial isolated pituitary adenomas kindreds for PAM variants. PAM encodes a multifunctional protein responsible for the essential C-terminal amidation of secreted peptides. Genetic screening was performed by germline and tumor sequencing and germline copy number variation (CNV) analysis. No germline CNVs or somatic single nucleotide variants (SNVs) were identified. We detected seven likely pathogenic heterozygous missense, truncating, and regulatory SNVs. These SNVs were found in sporadic subjects with GH excess (p.Gly552Arg and p.Phe759Ser), pediatric Cushing disease (c.-133T>C and p.His778fs), or with different types of PAs (c.-361G>A, p.Ser539Trp, and p.Asp563Gly). The SNVs were functionally tested in vitro for protein expression and trafficking by Western blotting, for splicing by minigene assays, and for amidation activity in cell lysates and serum samples. These analyses confirmed a deleterious effect on protein expression and/or function. By interrogating 200,000 exomes from the UK Biobank, we confirmed a significant association of the PAM gene and rare PAM SNVs to diagnoses linked to pituitary gland hyperfunction. Identification of PAM as a candidate gene associated with pituitary hypersecretion opens the possibility of developing novel therapeutics based on altering PAM function.

Extensive diversity in RNA termination and regulation revealed by transcriptome mapping for the Lyme pathogen B. burgdorferi.

Transcription termination is an essential and dynamic process that can tune gene expression in response to diverse molecular signals. Yet, the genomic positions, molecular mechanisms, and regulatory consequences of termination have only been studied thoroughly in model bacteria. We employed complementary RNA-seq approaches to map RNA ends for the transcriptome of the spirochete Borrelia burgdorferi - the etiological agent of Lyme disease. By systematically mapping B. burgdorferi RNA ends at single nucleotide resolution, we delineated complex gene arrangements and operons and mapped untranslated regions (UTRs) and small RNAs (sRNAs). We experimentally tested modes of B. burgdorferi transcription termination and compared our findings to observations in E. coli , P. aeruginosa , and B. subtilis . We discovered 63% of B. burgdorferi RNA 3' ends map upstream or internal to open reading frames (ORFs), suggesting novel mechanisms of regulation. Northern analysis confirmed the presence of stable 5' derived RNAs from mRNAs encoding gene products involved in the unique infectious cycle of B. burgdorferi . We suggest these RNAs resulted from premature termination and regulatory events, including forms of cis- acting regulation. For example, we documented that the polyamine spermidine globally influences the generation of truncated mRNAs. In one case, we showed that high spermidine concentrations increased levels of RNA fragments derived from an mRNA encoding a spermidine import system, with a concomitant decrease in levels of the full- length mRNA. Collectively, our findings revealed new insight into transcription termination and uncovered an abundance of potential RNA regulators.

Enhancer-promoter interactions can bypass CTCF-mediated boundaries and contribute to phenotypic robustness.

How enhancers activate their distal target promoters remains incompletely understood. Here we dissect how CTCF-mediated loops facilitate and restrict such regulatory interactions. Using an allelic series of mouse mutants, we show that CTCF is neither required for the interaction of the Sox2 gene with distal enhancers, nor for its expression. Insertion of various combinations of CTCF motifs, between Sox2 and its distal enhancers, generated boundaries with varying degrees of insulation that directly correlated with reduced transcriptional output. However, in both epiblast and neural tissues, enhancer contacts and transcriptional induction could not be fully abolished, and insertions failed to disrupt implantation and neurogenesis. In contrast, Sox2 expression was undetectable in the anterior foregut of mutants carrying the strongest boundaries, and these animals fully phenocopied loss of SOX2 in this tissue. We propose that enhancer clusters with a high density of regulatory activity can better overcome physical barriers to maintain faithful gene expression and phenotypic robustness.

Profiles of Wnt pathway gene expression during tooth morphogenesis.

Mouse and human genetic studies indicate key roles of the Wnt10a ligand in odontogenesis. Previous studies have identified effectors and regulators of the Wnt signaling pathway actively expressed during key stages of tooth morphogenesis. However, limitations in multiplexing and spatial resolution hindered a more comprehensive analysis of these signaling molecules. Here, profiling of transcriptomes using fluorescent multiplex in situ hybridization and single-cell RNA-sequencing (scRNA-seq) provide robust insight into the synchronized expression patterns of Wnt10a, Dkk1, and Sost simultaneously during tooth development. First, we identified Wnt10a transcripts restricted to the epithelium at the stage of tooth bud morphogenesis, contrasting that of Sost and Dkk1 localization to the dental mesenchyme. By embryonic day 15.5 (E15.5), a marked shift of Wnt10a expression from dental epithelium to mesenchyme was noted, while Sost and Dkk1 expression remained enriched in the mesenchyme. By postnatal day 0 (P0), co-localization patterns of Wnt10a, Dkk1, and Sost were observed in both terminally differentiating and secreting odontoblasts of molars and incisors. Interestingly, Wnt10a exhibited robust expression in fully differentiated ameloblasts at the developing cusp tip of both molars and incisors, an observation not previously noted in prior studies. At P7 and 14, after the mineralization of dentin and enamel, Wnt10a expression was limited to odontoblasts. Meanwhile, Wnt modulators showed reduced or absent signals in molars. In contrast, strong signals persisted in ameloblasts (for Wnt10a) and odontoblasts (for Wnt10a, Sost, and Dkk1) towards the proximal end of incisors, near the cervical loop. Our scRNA-seq analysis used CellChat to further contextualize Wnt pathway-mediated communication between cells by examining ligand-receptor interactions among different clusters. The co-localization pattern of Wnt10a, Dkk1, and Sost in both terminally differentiating and secreting odontoblasts of molars and incisors potentially signifies the crucial ligand-modulator interaction along the gradient of cytodifferentiation starting from each cusp tip towards the apical region. These data provide cell type-specific insight into the role of Wnt ligands and mediators during epithelial-mesenchymal interactions in odontogenesis.

SMAD3 mutation in LDS3 causes bone fragility by impairing the TGF-ß pathway and enhancing osteoclastogenesis.

Loss-of-function mutations in SMAD3 cause Loeys-Dietz syndrome type 3 (LDS3), a rare autosomal-dominant connective tissue disorder characterized by vascular pathology and skeletal abnormalities. Dysregulation of TGF-ß/SMAD signaling is associated with abnormal skeletal features and bone fragility. To date, histomorphometric and ultrastructural characteristics of bone with SMAD3 mutations have not been reported in humans and the exact mechanism by which SMAD3 mutations cause the LDS3 phenotype is poorly understood. Here, we investigated bone histomorphometry and matrix mineralization in human bone with a SMAD3 mutation and explored the associated cellular defect in the TGF-ß/SMAD pathway in vitro. The index patient had recurrent fractures, mild facial dysmorphism, arachnodactyly, pectus excavatum, chest asymmetry and kyphoscoliosis. Bone histomorphometry revealed markedly reduced cortical thickness (-68 %), trabecular thickness (-32 %), bone formation rate (-50 %) and delayed mineralization. Quantitative backscattered electron imaging demonstrated undermineralized bone matrix with increased heterogeneity in mineralization. The patient's SMAD3 mutation (c.200 T > G; p.I67S), when expressed from plasmid vectors in HEK293 cells, showed reduced phosphorylation and transcription factor activity compared to normal control and SMAD3 (p.S264Y), a gain-of-function mutation, somatic mosaicism of which causes melorheostosis. Transfection study of the patients' SMAD3 (p.I67S) mutation displayed lower luciferase reporter activity than normal SMAD3 and reduced expression of TGF-ß signaling target genes. Patient fibroblasts also demonstrated impaired SMAD3 protein stability. Osteoclastogenic differentiation significantly increased and osteoclast-associated genes, including ACP5 (encoding TRAP), ATP6V0D2, and DCSTAMP, were up-regulated in CD14 (+) peripheral blood mononuclear cells (PBMCs) with the SMAD3 (p.I67S) mutation. Upregulation of osteoclastogenic genes was associated with decreased expression of TGF-ß signaling target genes. We conclude that bone with the SMAD3 (p.I67S) mutation features reduced bone formation, and our functional studies revealed decreased SMAD3 activation and protein stability as well as increased osteoclastogenesis. These findings enhance our understanding of the pathophysiology of LDS3 caused by SMAD3 mutations. Emerging therapies targeting in the TGF-ß/SMAD pathway also raise hope for treatment of LDS3.

An epigenome atlas of neural progenitors within the embryonic mouse forebrain.

A comprehensive characterization of epigenomic organization in the embryonic mouse forebrain will enhance our understanding of neurodevelopment and provide insight into mechanisms of neurological disease. Here we collected single-cell chromatin accessibility profiles from four distinct neurogenic regions of the embryonic mouse forebrain using single nuclei ATAC-Seq (snATAC-Seq). We identified thousands of differentially accessible peaks, many restricted to distinct progenitor cell types or brain regions. We integrated snATAC-Seq and single cell transcriptome data to characterize changes of chromatin accessibility at enhancers and promoters with associated transcript abundance. Multi-modal integration of histone modifications (CUT&Tag and CUT&RUN), promoter-enhancer interactions (Capture-C) and high-order chromatin structure (Hi-C) extended these initial observations. This dataset reveals a diverse chromatin landscape with region-specific regulatory mechanisms and genomic interactions in distinct neurogenic regions of the embryonic mouse brain and represents an extensive public resource of a 'ground truth' epigenomic landscape at this critical stage of neurogenesis.

Extensive co-binding and rapid redistribution of NANOG and GATA6 during emergence of divergent lineages.

Fate-determining transcription factors (TFs) can promote lineage-restricted transcriptional programs from common progenitor states. The inner cell mass (ICM) of mouse blastocysts co-expresses the TFs NANOG and GATA6, which drive the bifurcation of the ICM into either the epiblast (Epi) or the primitive endoderm (PrE), respectively. Here, we induce GATA6 in embryonic stem cells-that also express NANOG-to characterize how a state of co-expression of opposing TFs resolves into divergent lineages. Surprisingly, we find that GATA6 and NANOG co-bind at the vast majority of Epi and PrE enhancers, a phenomenon we also observe in blastocysts. The co-bound state is followed by eviction and repression of Epi TFs, and quick remodeling of chromatin and enhancer-promoter contacts thus establishing the PrE lineage while repressing the Epi fate. We propose that co-binding of GATA6 and NANOG at shared enhancers maintains ICM plasticity and promotes the rapid establishment of Epi- and PrE-specific transcriptional programs.

Correlation of age of onset and clinical severity in Niemann-Pick disease type C1 with lysosomal abnormalities and gene expression.

Niemann-Pick disease type C1 (NPC1) is a rare, prematurely fatal lysosomal storage disorder which exhibits highly variable severity and disease progression as well as a wide-ranging age of onset, from perinatal stages to adulthood. This heterogeneity has made it difficult to obtain prompt diagnosis and to predict disease course. In addition, small NPC1 patient sample sizes have been a limiting factor in acquiring genome-wide transcriptome data. In this study, primary fibroblasts from an extensive cohort of 41 NPC1 patients were used to validate our previous findings that the lysosomal quantitative probe LysoTracker can be used as a predictor for age of onset and disease severity. We also examined the correlation between these clinical parameters and RNA expression data from primary fibroblasts and identified a set of genes that were significantly associated with lysosomal defects or age of onset, in particular neurological symptom onset. Hierarchical clustering showed that these genes exhibited distinct expression patterns among patient subgroups. This study is the first to collect transcriptomic data on such a large scale in correlation with clinical and cellular phenotypes, providing a rich genomic resource to address NPC1 clinical heterogeneity and discover potential biomarkers, disease modifiers, or therapeutic targets.

Transcriptional heterogeneity of ventricular zone cells in the ganglionic eminences of the mouse forebrain.

The ventricular zone (VZ) of the nervous system contains radial glia cells that were originally considered relatively homogenous in their gene expression, but a detailed characterization of transcriptional diversity in these VZ cells has not been reported. Here, we performed single-cell RNA sequencing to characterize transcriptional heterogeneity of neural progenitors within the VZ and subventricular zone (SVZ) of the ganglionic eminences (GEs), the source of all forebrain GABAergic neurons. By using a transgenic mouse line to enrich for VZ cells, we characterize significant transcriptional heterogeneity, both between GEs and within spatial subdomains of specific GEs. Additionally, we observe differential gene expression between E12.5 and E14.5 VZ cells, which could provide insights into temporal changes in cell fate. Together, our results reveal a previously unknown spatial and temporal genetic diversity of VZ cells in the ventral forebrain that will aid our understanding of initial fate decisions in the forebrain.

Transcriptome of HPßCD-treated Niemann-Pick disease type C1 cells highlights GPNMB as a biomarker for therapeutics.

The rare, fatal neurodegenerative disorder Niemann-Pick disease type C1 (NPC1) arises from lysosomal accumulation of unesterified cholesterol and glycosphingolipids. These subcellular pathologies lead to phenotypes of hepatosplenomegaly, neurological degeneration and premature death. The timing and severity of NPC1 clinical presentation is extremely heterogeneous. This study analyzed RNA-Seq data from 42 NPC1 patient-derived, primary fibroblast cell lines to determine transcriptional changes induced by treatment with 2-hydroxypropyl-ß-cyclodextrin (HPßCD), a compound currently under investigation in clinical trials. A total of 485 HPßCD-responsive genes were identified. Pathway enrichment analysis of these genes showed significant involvement in cholesterol and lipid biosynthesis. Furthermore, immunohistochemistry of the cerebellum as well as measurements of plasma from Npc1m1N null mice treated with HPßCD and adeno-associated virus gene therapy suggests that one of the identified genes, GPNMB, may serve as a useful biomarker of treatment response in NPC1 disease. Overall, this large NPC1 patient-derived dataset provides a comprehensive foundation for understanding the genomic response to HPßCD treatment.