Large-Scale Exome Sequencing Study Implicates Both Developmental and Functional Changes in the Neurobiology of Autism

We present the largest exome sequencing study of autism spectrum disorder (ASD) to date (n = 35,584 total samples, 11,986 with ASD). Using an enhanced analytical framework to integrate de novo and case-control rare variation, we identify 102 risk genes at a false discovery rate of 0.1 or less. Of these genes, 49 show higher frequencies of disruptive de novo variants in individuals ascertained to have severe neurodevelopmental delay, whereas 53 show higher frequencies in individuals ascertained to have ASD; comparing ASD cases with mutations in these groups reveals phenotypic differences. Expressed early in brain development, most risk genes have roles in regulation of gene expression or neuronal communication (i.e., mutations effect neurodevelopmental and neurophysiological changes), and 13 fall within loci recurrently hit by copy number variants. In cells from the human cortex, expression of risk genes is enriched in excitatory and inhibitory neuronal lineages, consistent with multiple paths to an excitatory-inhibitory imbalance underlying ASD.

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.

Single-cell genomics identifies cell type-specific molecular changes in autism

Despite the clinical and genetic heterogeneity of autism, bulk gene expression studies show that changes in the neocortex of autism patients converge on common genes and pathways. However, direct assessment of specific cell types in the brain affected by autism has not been feasible until recently. We used single-nucleus RNA sequencing of cortical tissue from patients with autism to identify autism-associated transcriptomic changes in specific cell types. We found that synaptic signaling of upper-layer excitatory neurons and the molecular state of microglia are preferentially affected in autism. Moreover, our results show that dysregulation of specific groups of genes in cortico-cortical projection neurons correlates with clinical severity of autism. These findings suggest that molecular changes in upper-layer cortical circuits are linked to behavioral manifestations of autism.

Multimodal Analysis of Composition and Spatial Architecture in Human Squamous Cell Carcinoma

To define the cellular composition and architecture of cutaneous squamous cell carcinoma (cSCC), we combined single-cell RNA sequencing with spatial transcriptomics and multiplexed ion beam imaging from a series of human cSCCs and matched normal skin. cSCC exhibited four tumor subpopulations, three recapitulating normal epidermal states, and a tumor-specific keratinocyte (TSK) population unique to cancer, which localized to a fibrovascular niche. Integration of single-cell and spatial data mapped ligand-receptor networks to specific cell types, revealing TSK cells as a hub for intercellular communication. Multiple features of potential immunosuppression were observed, including T regulatory cell (Treg) co-localization with CD8 T cells in compartmentalized tumor stroma. Finally, single-cell characterization of human tumor xenografts and in vivo CRISPR screens identified essential roles for specific tumor subpopulation-enriched gene networks in tumorigenesis. These data define cSCC tumor and stromal cell subpopulations, the spatial niches where they interact, and the communicating gene networks that they engage in cancer.

Single-cell profiling of human gliomas reveals macrophage ontogeny as a basis for regional differences in macrophage activation in the tumor microenvironment

Background

Tumor-associated macrophages (TAMs) are abundant in gliomas and immunosuppressive TAMs are a barrier to emerging immunotherapies. It is unknown to what extent macrophages derived from peripheral blood adopt the phenotype of brain-resident microglia in pre-treatment gliomas. The relative proportions of blood-derived macrophages and microglia have been poorly quantified in clinical samples due to a paucity of markers that distinguish these cell types in malignant tissue.

Results

We perform single-cell RNA-sequencing of human gliomas and identify phenotypic differences in TAMs of distinct lineages. We isolate TAMs from patient biopsies and compare them with macrophages from non-malignant human tissue, glioma atlases, and murine glioma models. We present a novel signature that distinguishes TAMs by ontogeny in human gliomas. Blood-derived TAMs upregulate immunosuppressive cytokines and show an altered metabolism compared to microglial TAMs. They are also enriched in perivascular and necrotic regions. The gene signature of blood-derived TAMs, but not microglial TAMs, correlates with significantly inferior survival in low-grade glioma. Surprisingly, TAMs frequently co-express canonical pro-inflammatory (M1) and alternatively activated (M2) genes in individual cells.

Conclusions

We conclude that blood-derived TAMs significantly infiltrate pre-treatment gliomas, to a degree that varies by glioma subtype and tumor compartment. Blood-derived TAMs do not universally conform to the phenotype of microglia, but preferentially express immunosuppressive cytokines and show an altered metabolism. Our results argue against status quo therapeutic strategies that target TAMs indiscriminately and in favor of strategies that specifically target immunosuppressive blood-derived TAMs.

Electronic supplementary material

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

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.

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.

Single-Cell Analyses Identify Brain Mural Cells Expressing CD19 as Potential Off-Tumor Targets for CAR-T Immunotherapies

CD19-directed immunotherapies are clinically effective for treating B cell malignancies but also cause a high incidence of neurotoxicity. A subset of patients treated with chimeric antigen receptor (CAR) T cells or bispecific T cell engager (BiTE) antibodies display severe neurotoxicity, including fatal cerebral edema associated with T cell infiltration into the brain. Here, we report that mural cells, which surround the endothelium and are critical for blood-brain-barrier integrity, express CD19. We identify CD19 expression in brain mural cells using single-cell RNA sequencing data and confirm perivascular staining at the protein level. CD19 expression in the brain begins early in development alongside the emergence of mural cell lineages and persists throughout adulthood across brain regions. Mouse mural cells demonstrate lower levels of Cd19 expression, suggesting limitations in preclinical animal models of neurotoxicity. These data suggest an on-target mechanism for neurotoxicity in CD19-directed therapies and highlight the utility of human single-cell atlases for designing immunotherapies.

Development and Arealization of the Cerebral Cortex

Adult cortical areas consist of specialized cell types and circuits that support unique higher-order cognitive functions. How this regional diversity develops from an initially uniform neuroepithelium has been the subject of decades of seminal research, and emerging technologies, including single-cell transcriptomics, provide a new perspective on area-specific molecular diversity. Here, we review the early developmental processes that underlie cortical arealization, including both cortex intrinsic and extrinsic mechanisms as embodied by the protomap and protocortex hypotheses, respectively. We propose an integrated model of serial homology whereby intrinsic genetic programs and local factors establish early transcriptomic differences between excitatory neurons destined to give rise to broad "proto-regions," and activity-dependent mechanisms lead to progressive refinement and formation of sharp boundaries between functional areas. Finally, we explore the potential of these basic developmental processes to inform our understanding of the emergence of functional neural networks and circuit abnormalities in neurodevelopmental disorders.

Single-cell atlas of early human brain development highlights heterogeneity of human neuroepithelial cells and early radial glia

The human cortex comprises diverse cell types that emerge from an initially uniform neuroepithelium that gives rise to radial glia, the neural stem cells of the cortex. To characterize the earliest stages of human brain development, we performed single-cell RNA-sequencing across regions of the developing human brain, including the telencephalon, diencephalon, midbrain, hindbrain and cerebellum. We identify nine progenitor populations physically proximal to the telencephalon, suggesting more heterogeneity than previously described, including a highly prevalent mesenchymal-like population that disappears once neurogenesis begins. Comparison of human and mouse progenitor populations at corresponding stages identifies two progenitor clusters that are enriched in the early stages of human cortical development. We also find that organoid systems display low fidelity to neuroepithelial and early radial glia cell types, but improve as neurogenesis progresses. Overall, we provide a comprehensive molecular and spatial atlas of early stages of human brain and cortical development.

Genomic analysis of mycosis fungoides and Sézary syndrome identifies recurrent alterations in TNFR2

Mycosis fungoides and Sézary syndrome comprise the majority of cutaneous T cell lymphomas (CTCLs), disorders notable for their clinical heterogeneity that can present in skin or peripheral blood. Effective treatment options for CTCL are limited, and the genetic basis of these T cell lymphomas remains incompletely characterized. Here we report recurrent point mutations and genomic gains of TNFRSF1B, encoding the tumor necrosis factor receptor TNFR2, in 18% of patients with mycosis fungoides and Sézary syndrome. Expression of the recurrent TNFR2 Thr377Ile mutant in T cells leads to enhanced non-canonical NF-κB signaling that is sensitive to the proteasome inhibitor bortezomib. Using an integrative genomic approach, we additionally discovered a recurrent CTLA4-CD28 fusion, as well as mutations in downstream signaling mediators of these receptors.

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.

Origins and Proliferative States of Human Oligodendrocyte Precursor Cells

Human cerebral cortex size and complexity has increased greatly during evolution. While increased progenitor diversity and enhanced proliferative potential play important roles in human neurogenesis and gray matter expansion, the mechanisms of human oligodendrogenesis and white matter expansion remain largely unknown. Here, we identify EGFR-expressing "Pre-OPCs" that originate from outer radial glial cells (oRGs) and undergo mitotic somal translocation (MST) during division. oRG-derived Pre-OPCs provide an additional source of human cortical oligodendrocyte precursor cells (OPCs) and define a lineage trajectory. We further show that human OPCs undergo consecutive symmetric divisions to exponentially increase the progenitor pool size. Additionally, we find that the OPC-enriched gene, PCDH15, mediates daughter cell repulsion and facilitates proliferation. These findings indicate properties of OPC derivation, proliferation, and dispersion important for human white matter expansion and myelination.

An atlas of cortical arealization identifies dynamic molecular signatures

The human brain is subdivided into distinct anatomical structures, including the neocortex, which in turn encompasses dozens of distinct specialized cortical areas. Early morphogenetic gradients are known to establish early brain regions and cortical areas, but how early patterns result in finer and more discrete spatial differences remains poorly understood¹. Here we use single-cell RNA sequencing to profile ten major brain structures and six neocortical areas during peak neurogenesis and early gliogenesis. Within the neocortex, we find that early in the second trimester, a large number of genes are differentially expressed across distinct cortical areas in all cell types, including radial glia, the neural progenitors of the cortex. However, the abundance of areal transcriptomic signatures increases as radial glia differentiate into intermediate progenitor cells and ultimately give rise to excitatory neurons. Using an automated, multiplexed single-molecule fluorescent in situ hybridization approach, we find that laminar gene-expression patterns are highly dynamic across cortical regions. Together, our data suggest that early cortical areal patterning is defined by strong, mutually exclusive frontal and occipital gene-expression signatures, with resulting gradients giving rise to the specification of areas between these two poles throughout successive developmental timepoints.

Quantitative analysis of mammalian translation initiation sites by FACS-seq

An approach combining fluorescence-activated cell sorting and high-throughput DNA sequencing (FACS-seq) was employed to determine the efficiency of start codon recognition for all possible translation initiation sites (TIS) utilizing AUG start codons. Using FACS-seq, we measured translation from a genetic reporter library representing all 65,536 possible TIS sequences spanning the −6 to +5 positions. We found that the motif RYMRMVAUGGC enhanced start codon recognition and translation efficiency. However, dinucleotide interactions, which cannot be conveyed by a single motif, were also important for modeling TIS efficiency. Our dataset combined with modeling allowed us to predict genome-wide translation initiation efficiency for all mRNA transcripts. Additionally, we screened somatic TIS mutations associated with tumorigenesis to identify candidate driver mutations consistent with known tumor expression patterns. Finally, we implemented a quantitative leaky scanning model to predict alternative initiation sites that produce truncated protein isoforms and compared predictions with ribosome footprint profiling data. The comprehensive analysis of the TIS sequence space enables quantitative predictions of translation initiation based on genome sequence.

The Functional Proximal Proteome of Oncogenic Ras Includes mTORC2

Proximity-dependent biotin labeling (BioID) may identify new targets for cancers driven by difficult-to-drug oncogenes such as Ras. Therefore, BioID was used with wild-type (WT) and oncogenic mutant (MT) H-, K-, and N-Ras, identifying known interactors, including Raf and PI3K, as well as a common set of 130 novel proteins proximal to all Ras isoforms. A CRISPR screen of these proteins for Ras dependence identified mTOR, which was also found proximal to MT Ras in human tumors. Oncogenic Ras directly bound two mTOR complex 2 (mTORC2) components, mTOR and MAPKAP1, to promote mTORC2 kinase activity at the plasma membrane. mTORC2 enabled the Ras pro-proliferative cell cycle transcriptional program, and perturbing the Ras-mTORC2 interaction impaired Ras-dependent neoplasia in vivo. Combining proximity-dependent proteomics with CRISPR screening identified a new set of functional Ras-associated proteins, defined mTORC2 as a new direct Ras effector, and offers a strategy for finding new proteins that cooperate with dominant oncogenes.

Transient epileptic amnesia: regional brain atrophy and its relationship to memory deficits

Transient epileptic amnesia (TEA) is a recently recognised form of epilepsy of which the principle manifestation is recurrent, transient episodes of isolated memory loss. In addition to the amnesic episodes, many patients describe significant interictal memory difficulties. Performance on standard neuropsychological tests is often normal. However, two unusual forms of memory deficit have recently been demonstrated in TEA: (i) accelerated long-term forgetting (ALF): the excessively rapid loss of newly acquired memories over a period of days or weeks and (ii) remote autobiographical memory loss: a loss of memories for salient, personally experienced events of the past few decades. The neuroanatomical bases of TEA and its associated memory deficits are unknown. In this study, we first assessed the relationship between subjective and objective memory performance in 41 patients with TEA. We then analysed MRI data from these patients and 20 matched healthy controls, using manual volumetry and voxel-based morphometry (VBM) to correlate regional brain volumes with clinical and neuropsychological data. Subjective memory estimates were unrelated to performance on standard neuropsychological tests but were partially predicted by mood, ALF and remote autobiographical memory. Manual volumetry identified subtle hippocampal volume loss in the patient group. Both manual volumetry and VBM revealed correlations between medial temporal lobe atrophy and standard anterograde memory scores, but no relation between atrophy and ALF or remote autobiographical memory. These results add weight to the hypothesis that TEA is a syndrome of mesial temporal lobe epilepsy. Furthermore, they suggest that although standard anterograde memory test performance is related to the degree of mesial temporal lobe damage, this is not true for ALF and autobiographical amnesia. It is possible that these unusual memory deficits have a more diffuse physiological basis rather than being a consequence of discrete structural damage.

Recurrent point mutations in the kinetochore gene KNSTRN in cutaneous squamous cell carcinoma

Here we report the discovery of recurrent mutations concentrated at an ultraviolet signature hotspot in KNSTRN, which encodes a kinetochore protein, in 19% of cutaneous squamous cell carcinomas (SCCs). Cancer-associated KNSTRN mutations, most notably those encoding p.Ser24Phe, disrupt chromatid cohesion in normal cells, occur in SCC precursors, correlate with increased aneuploidy in primary tumors and enhance tumorigenesis in vivo. These findings suggest a role for KNSTRN mutagenesis in SCC development.

Rapid identification of non-human sequences in high-throughput sequencing datasets

Rapid identification of non-human sequences (RINS) is an intersection-based pathogen detection workflow that utilizes a user-provided custom reference genome set for identification of non-human sequences in deep sequencing datasets. In <2 h, RINS correctly identified the known virus in the dataset SRR73726 and is compatible with any computer capable of running the prerequisite alignment and assembly programs. RINS accurately identifies sequencing reads from intact or mutated non-human genomes in a dataset and robustly generates contigs with these non-human sequences (Supplementary Material).

AVAILABILITY

RINS is available for free download at http://khavarilab.stanford.edu/resources.html.

Mechanism of action of amphotericin B on Leishmania donovani promastigotes

The growth of Leishmania donovani promastigotes in a liquid medium was completely inhibited by amphotericin B at a concentration of 0.3 microgram ml-1 (0.3 microM). Continuous release of small molecules that absorb at 260 nm and 280 nm was observed after contact with the drug. Uptake of [U-14C]glucose was inhibited in cells treated with the drug. An immediate release of isotopic glucose and its metabolites from preloaded cells could be demonstrated after incubation with amphotericin B (0.4 microM). Inhibition of respiration by the drug was a comparatively slower process. All the above effects could be effectively prevented in the presence of either cholesterol or ergosterol. The primary site of action of amphotericin B on L. donovani promastigote cells appears to be membrane sterols that result in a loss of the permeability barrier to small metabolites. An interesting biochemical similarity, thus, emerges between flagellated protozoa and fungi.

Shared and derived features of cellular diversity in the human cerebral cortex

The cerebral cortex is the hallmark of the mammalian nervous system, and its large size and cellular diversity in humans support our most sophisticated cognitive abilities. Although the basic cellular organization of the cortex is conserved across mammals, cells have diversified during evolution. An increasingly integrated taxonomy of cell types, especially with the advent of single-cell transcriptomic data, has revealed an unprecedented variety of human cortical cell subtypes. Here, we broadly review the cellular composition and diversity of the mammalian brain, and how progenitor pools generate cell subtypes during development. We then discuss human cortical cells that are distinct from rodent cells, as well as the challenges and advantages of using model systems to study human cell types in health and disease.

Ensembles of endothelial and mural cells promote angiogenesis in prenatal human brain

Interactions between angiogenesis and neurogenesis regulate embryonic brain development. However, a comprehensive understanding of the stages of vascular cell maturation is lacking, especially in the prenatal human brain. Using fluorescence-activated cell sorting, single-cell transcriptomics, and histological and ultrastructural analyses, we show that an ensemble of endothelial and mural cell subtypes tile the brain vasculature during the second trimester. These vascular cells follow distinct developmental trajectories and utilize diverse signaling mechanisms, including collagen, laminin, and midkine, to facilitate cell-cell communication and maturation. Interestingly, our results reveal that tip cells, a subtype of endothelial cells, are highly enriched near the ventricular zone, the site of active neurogenesis. Consistent with these observations, prenatal vascular cells transplanted into cortical organoids exhibit restricted lineage potential that favors tip cells, promotes neurogenesis, and reduces cellular stress. Together, our results uncover important mechanisms into vascular maturation during this critical period of human brain development.

Black tea is a powerful chemopreventor of reactive oxygen and nitrogen species: comparison with its individual catechin constituents and green tea

Production of black tea [BT] results in biotransformation of catechins of green tea [GT] to theaflavins and thearubigins. BT was found to be more efficient than GT and its individual catechin constituents in proportionate amounts in abrogating production of NO and O2(-) in activated murine peritoneal macrophages. In a reconstitution system of BT that is free of all catechins, stepwise addition of catechins showed that though all the constituents contributed to the overall effect of BT, theaflavin was the most powerful in abrogating NO production. RT-PCR analysis also showed theaflavin to be the most important constituent in down-regulating synthesis of iNOS. Clearly, BT containing theaflavin is an excellent chemopreventor against reactive oxygen and nitrogen species.

Next-generation sequencing of idiopathic multicentric and unicentric Castleman disease and follicular dendritic cell sarcomas

Castleman disease (CD) is a rare lymphoproliferative disorder subclassified as unicentric CD (UCD) or multicentric CD (MCD) based on clinical features and the distribution of enlarged lymph nodes with characteristic histopathology. MCD can be further subtyped based on human herpes virus 8 (HHV8) infection into HHV8-associated MCD, HHV8-/idiopathic MCD (iMCD), and polyneuropathy, organomegaly, endocrinopathy, monoclonal gammopathy, and skin change (POEMS)-associated MCD. In a subset of cases of UCD, an associated follicular dendritic cell sarcoma (FDCS) may be seen. Although numerous reports of the clinical and histologic features of UCD, MCD, and FDCS exist, an understanding of the genetic and epigenetic landscape of these rare diseases is lacking. Given this paucity of knowledge, we analyzed 15 cases of UCD and 3 cases of iMCD by targeted next-generation sequencing (NGS; 405 genes) and 3 cases of FDCS associated with UCD hyaline vascular variant (UCD-HVV) by whole-exome sequencing. Common amplifications of ETS1, PTPN6, and TGFBR2 were seen in 1 iMCD and 1 UCD case; the iMCD case also had a somatic DNMT3A L295Q mutation. This iMCD patient also showed clinicopathologic features consistent with a specific subtype known as Castleman-Kojima disease (thrombocytopenia, anasarca, fever, reticulin fibrosis, and organomegaly [TAFRO] clinical subtype). Additionally, 1 case of UCD-HVV showed amplification of the cluster of histone genes on chromosome 6p. FDCS associated with UCD-HVV showed mutations and copy number changes in known oncogenes, tumor suppressors, and chromatin structural-remodeling proteins.