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.

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.

Single-cell RNA-seq reveals dynamic paracrine control of cellular variation.

High-throughput single-cell transcriptomics offers an unbiased approach for understanding the extent, basis and function of gene expression variation between seemingly identical cells. Here we sequence single-cell RNA-seq libraries prepared from over 1,700 primary mouse bone-marrow-derived dendritic cells spanning several experimental conditions. We find substantial variation between identically stimulated dendritic cells, in both the fraction of cells detectably expressing a given messenger RNA and the transcript's level within expressing cells. Distinct gene modules are characterized by different temporal heterogeneity profiles. In particular, a 'core' module of antiviral genes is expressed very early by a few 'precocious' cells in response to uniform stimulation with a pathogenic component, but is later activated in all cells. By stimulating cells individually in sealed microfluidic chambers, analysing dendritic cells from knockout mice, and modulating secretion and extracellular signalling, we show that this response is coordinated by interferon-mediated paracrine signalling from these precocious cells. Notably, preventing cell-to-cell communication also substantially reduces variability between cells in the expression of an early-induced 'peaked' inflammatory module, suggesting that paracrine signalling additionally represses part of the inflammatory program. Our study highlights the importance of cell-to-cell communication in controlling cellular heterogeneity and reveals general strategies that multicellular populations can use to establish complex dynamic responses.

Single molecule quantitation and sequencing of rare translocations using microfluidic nested digital PCR.

Cancers are heterogeneous and genetically unstable. New methods are needed that provide the sensitivity and specificity to query single cells at the genetic loci that drive cancer progression, thereby enabling researchers to study the progression of individual tumors. Here, we report the development and application of a bead-based hemi-nested microfluidic droplet digital PCR (dPCR) technology to achieve 'quantitative' measurement and single-molecule sequencing of somatically acquired carcinogenic translocations at extremely low levels (<10(-6)) in healthy subjects. We use this technique in our healthy study population to determine the overall concentration of the t(14;18) translocation, which is strongly associated with follicular lymphoma. The nested dPCR approach improves the detection limit to 1×10(-7) or lower while maintaining the analysis efficiency and specificity. Further, the bead-based dPCR enabled us to isolate and quantify the relative amounts of the various clonal forms of t(14;18) translocation in these subjects, and the single-molecule sensitivity and resolution of dPCR led to the discovery of new clonal forms of t(14;18) that were otherwise masked by the conventional quantitative PCR measurements. In this manner, we created a quantitative map for this carcinogenic mutation in this healthy population and identified the positions on chromosomes 14 and 18 where the vast majority of these t(14;18) events occur.

Modulation of Ras signaling alters the toxicity of hydroquinone, a benzene metabolite and component of cigarette smoke.

BACKGROUND: Benzene is an established human leukemogen, with a ubiquitous environmental presence leading to significant population exposure. In a genome-wide functional screen in the yeast Saccharomyces cerevisiae, inactivation of IRA2, a yeast ortholog of the human tumor suppressor gene NF1 (Neurofibromin), enhanced sensitivity to hydroquinone, an important benzene metabolite. Increased Ras signaling is implicated as a causal factor in the increased pre-disposition to leukemia of individuals with mutations in NF1. METHODS: Growth inhibition of yeast by hydroquinone was assessed in mutant strains exhibiting varying levels of Ras activity. Subsequently, effects of hydroquinone on both genotoxicity (measured by micronucleus formation) and proliferation of WT and Nf1 null murine hematopoietic precursors were assessed. RESULTS: Here we show that the Ras status of both yeast and mammalian cells modulates hydroquinone toxicity, indicating potential synergy between Ras signaling and benzene toxicity. Specifically, enhanced Ras signaling increases both hydroquinone-mediated growth inhibition in yeast and genotoxicity in mammalian hematopoetic precursors as measured by an in vitro erythroid micronucleus assay. Hydroquinone also increases proliferation of CFU-GM progenitor cells in mice with Nf1 null bone marrow relative to WT, the same cell type associated with benzene-associated leukemia. CONCLUSIONS: Together our findings show that hydroquinone toxicity is modulated by Ras signaling. Individuals with abnormal Ras signaling could be more vulnerable to developing myeloid diseases after exposure to benzene. We note that hydroquinone is used cosmetically as a skin-bleaching agent, including by individuals with cafe-au-lait spots (which may be present in individuals with neurofibromatosis who have a mutation in NF1), which could be unadvisable given our findings.

High-performance single cell genetic analysis using microfluidic emulsion generator arrays.

High-throughput genetic and phenotypic analysis at the single cell level is critical to advance our understanding of the molecular mechanisms underlying cellular function and dysfunction. Here we describe a high-performance single cell genetic analysis (SCGA) technique that combines high-throughput microfluidic emulsion generation with single cell multiplex polymerase chain reaction (PCR). Microfabricated emulsion generator array (MEGA) devices containing 4, 32, and 96 channels are developed to confer a flexible capability of generating up to 3.4 x 10(6) nanoliter-volume droplets per hour. Hybrid glass-polydimethylsiloxane diaphragm micropumps integrated into the MEGA chips afford uniform droplet formation, controlled generation frequency, and effective transportation and encapsulation of primer functionalized microbeads and cells. A multiplex single cell PCR method is developed to detect and quantify both wild type and mutant/pathogenic cells. In this method, microbeads functionalized with multiple forward primers targeting specific genes from different cell types are used for solid-phase PCR in droplets. Following PCR, the droplets are lysed and the beads are pooled and rapidly analyzed by multicolor flow cytometry. Using Escherichia coli bacterial cells as a model, we show that this technique enables digital detection of pathogenic E. coli O157 cells in a high background of normal K12 cells, with a detection limit on the order of 1/10(5). This result demonstrates that multiplex SCGA is a promising tool for high-throughput quantitative digital analysis of genetic variation in complex populations.

Analysis of individual cells identifies cell-to-cell variability following induction of cellular senescence.

Senescent cells play important roles in both physiological and pathological processes, including cancer and aging. In all cases, however, senescent cells comprise only a small fraction of tissues. Senescent phenotypes have been studied largely in relatively homogeneous populations of cultured cells. In vivo, senescent cells are generally identified by a small number of markers, but whether and how these markers vary among individual cells is unknown. We therefore utilized a combination of single-cell isolation and a nanofluidic PCR platform to determine the contributions of individual cells to the overall gene expression profile of senescent human fibroblast populations. Individual senescent cells were surprisingly heterogeneous in their gene expression signatures. This cell-to-cell variability resulted in a loss of correlation among the expression of several senescence-associated genes. Many genes encoding senescence-associated secretory phenotype (SASP) factors, a major contributor to the effects of senescent cells in vivo, showed marked variability with a subset of highly induced genes accounting for the increases observed at the population level. Inflammatory genes in clustered genomic loci showed a greater correlation with senescence compared to nonclustered loci, suggesting that these genes are coregulated by genomic location. Together, these data offer new insights into how genes are regulated in senescent cells and suggest that single markers are inadequate to identify senescent cells in vivo.

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.

Corrigendum: Fluidic Logic Used in a Systems Approach to Enable Integrated Single-Cell Functional Analysis.

[This corrects the article on p. 70 in vol. 4, PMID: 27709111.].

Massively parallel digital transcriptional profiling of single cells.

Characterizing the transcriptome of individual cells is fundamental to understanding complex biological systems. We describe a droplet-based system that enables 3' mRNA counting of tens of thousands of single cells per sample. Cell encapsulation, of up to 8 samples at a time, takes place in ∼6 min, with ∼50% cell capture efficiency. To demonstrate the system's technical performance, we collected transcriptome data from ∼250k single cells across 29 samples. We validated the sensitivity of the system and its ability to detect rare populations using cell lines and synthetic RNAs. We profiled 68k peripheral blood mononuclear cells to demonstrate the system's ability to characterize large immune populations. Finally, we used sequence variation in the transcriptome data to determine host and donor chimerism at single-cell resolution from bone marrow mononuclear cells isolated from transplant patients.

Fluidic Logic Used in a Systems Approach to Enable Integrated Single-Cell Functional Analysis.

The study of single cells has evolved over the past several years to include expression and genomic analysis of an increasing number of single cells. Several studies have demonstrated wide spread variation and heterogeneity within cell populations of similar phenotype. While the characterization of these populations will likely set the foundation for our understanding of genomic- and expression-based diversity, it will not be able to link the functional differences of a single cell to its underlying genomic structure and activity. Currently, it is difficult to perturb single cells in a controlled environment, monitor and measure the response due to perturbation, and link these response measurements to downstream genomic and transcriptomic analysis. In order to address this challenge, we developed a platform to integrate and miniaturize many of the experimental steps required to study single-cell function. The heart of this platform is an elastomer-based integrated fluidic circuit that uses fluidic logic to select and sequester specific single cells based on a phenotypic trait for downstream experimentation. Experiments with sequestered cells that have been performed include on-chip culture, exposure to various stimulants, and post-exposure image-based response analysis, followed by preparation of the mRNA transcriptome for massively parallel sequencing analysis. The flexible system embodies experimental design and execution that enable routine functional studies of single cells.

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. Although cellular barcoding strategies allow parallel sequencing of single cells at ultra-low depths, the limitations of shallow sequencing have not been investigated directly. 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 the 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.

Bone marrow genotoxicity of 2,5-dimethylfuran, a green biofuel candidate.

2,5-Dimethylfuran (DMF) is being considered as a potential green transportation biofuel, but there is limited information about its toxicity and safety. We examined DMF toxicity in the bone marrow using a murine in vitro erythropoietic micronucleus assay and found that exposure to DMF (0.1 mM, 1 hr) induced an increase in micronuclei frequency compared with controls. These data suggest that DMF may be genotoxic to hematopoietic cells and that more thorough toxicological studies on DMF should be conducted to ensure public and occupational safety before it is considered a viable biofuel and produced in mass quantities. As well as specific data on DMF, our study further validates an in vitro cell culture system that captures the essential features of the in vivo mammalian micronucleus genotoxicity assay, enabling increased throughput and controlled studies on hematopoietic DNA damage response, while reducing animal sacrifice. In vitro assays, such as the in vitro micronucleus assay, will be essential as international chemical policy is increasingly utilizing green chemistry principles that require more toxicological testing.

Integrated molecular analysis indicates undetectable change in DNA damage in mice after continuous irradiation at ~ 400-fold natural background radiation.

BACKGROUND: In the event of a nuclear accident, people are exposed to elevated levels of continuous low dose-rate radiation. Nevertheless, most of the literature describes the biological effects of acute radiation. OBJECTIVES: DNA damage and mutations are well established for their carcinogenic effects. We assessed several key markers of DNA damage and DNA damage responses in mice exposed to low dose-rate radiation to reveal potential genotoxic effects associated with low dose-rate radiation. METHODS: We studied low dose-rate radiation using a variable low dose-rate irradiator consisting of flood phantoms filled with 125Iodine-containing buffer. Mice were exposed to 0.0002 cGy/min (~ 400-fold background radiation) continuously over 5 weeks. We assessed base lesions, micronuclei, homologous recombination (HR; using fluorescent yellow direct repeat mice), and transcript levels for several radiation-sensitive genes. RESULTS: We did not observe any changes in the levels of the DNA nucleobase damage products hypoxanthine, 8-oxo-7,8-dihydroguanine, 1,N6-ethenoadenine, or 3,N4-ethenocytosine above background levels under low dose-rate conditions. The micronucleus assay revealed no evidence that low dose-rate radiation induced DNA fragmentation, and there was no evidence of double strand break-induced HR. Furthermore, low dose-rate radiation did not induce Cdkn1a, Gadd45a, Mdm2, Atm, or Dbd2. Importantly, the same total dose, when delivered acutely, induced micronuclei and transcriptional responses. CONCLUSIONS: These results demonstrate in an in vivo animal model that lowering the dose-rate suppresses the potentially deleterious impact of radiation and calls attention to the need for a deeper understanding of the biological impact of low dose-rate radiation.

Assessment of genetic damage in healthy and diseased tissue.

DNA, along with other cellular components, is under constant attack by chemical, physical, and infectious agents present in the human environment, as well as by reactive metabolites generated by physiological processes. Mutations occur as the consequence of this damage, but may also be caused by improper DNA repair of alterations occurring during normal DNA replication and transcription. Genetic damage can occur at the level of the gene (e.g. point mutations, insertions, and deletions) or at the level of the chromosome (e.g. aneuploidy, translocations). Further, mutations can also take place in mitochondrial DNA. Another form of DNA modification is epigenetic methylation of CpG islands, which affects the dynamics of chromatin as well as the expression of a large panel of genes. Recent technical advances have improved the capacity to detect and quantify genetic and epigenetic changes. This chapter summarizes current knowledge on mechanisms of DNA damage and mutagenesis, laying out the concepts for interpreting mutations as biomarkers in investigating the causes and consequences of cancer. It also outlines both established and novel methods for detecting genetic and epigenetic changes in normal and diseased tissues, and then discusses their application in the realm of molecular epidemiology.

Selected technologies for measuring acquired genetic damage in humans.

Technical advances have improved the capacity to detect and quantify genetic variants, providing novel methods for the detection of rare mutations and for better understanding the underlying environmental factors and biological mechanisms contributing to mutagenesis. The polymerase chain reaction (PCR) has revolutionized genetic testing and remains central to many of these new techniques for mutation detection. Millions of genetic variations have been discovered across the genome. These variations include germline mutations and polymorphisms, which are inherited in a Mendelian manner and present in all cells, as well as acquired, somatic mutations that differ widely by type and size [from single-base mutations to whole chromosome rearrangements, and including submicroscopic copy number variations (CNVs)]. This review focuses on current methods for assessing acquired somatic mutations in the genome, and it examines their application in molecular epidemiology and sensitive detection and analysis of disease. Although older technologies have been exploited for detecting acquired mutations in cancer and other disease, the high-throughput and high-sensitivity offered by next-generation sequencing (NGS) systems are transforming the discovery of disease-associated acquired mutations by enabling comparative whole-genome sequencing of diseased and healthy tissues from the same individual. Emerging microfluidic technologies are beginning to facilitate single-cell genetic analysis of target variable regions for investigating cell heterogeneity within tumors as well as preclinical detection of disease. The technologies discussed in this review will significantly expand our knowledge of acquired genetic mutations and causative mechanisms.