Search

SUMMARYBrain endothelial cells (BECs) are key elements of the blood-brain barrier (BBB), protecting the brain from pathogens and restricting access to circulatory factors. Recent studies have demonstrated that the circulatory environment can modulate brain aging, yet, the underlying processes remain largely unknown. Given the BBB’s intermediary position, we hypothesized that BECs sense, adapt to, and relay signals between the aging blood and brain. We sequenced single endothelial cells from the hippocampus—a brain region key to learning, memory, and neurogenesis— of healthy young and aged mice as well as post-exposure to inflammatory and age-related circulatory factors. We discovered that aged capillary BECs, compared with arterial and venous cells, exhibit the greatest transcriptional changes, upregulating innate immunity, antigen presentation, TGF-β signaling and oxidative stress response pathways. Remarkably, short-term infusions of aged plasma into young mice recapitulated key aspects of this aging transcriptome, while infusions of young plasma into aged mice reversed select aging signatures, essentially rejuvenating the BBB endothelium transcriptome. We identify candidate pathways mediating blood-borne brain rejuvenation by comparing age-upregulated genes with those modulated by plasma exposure. Together, these findings suggest that the transcriptional age of BECs is exquisitely sensitive to age-related circulatory cues and pinpoint the BBB itself as a promising therapeutic target to treat brain disease.HighlightsSingle-cell RNA sequencing of brain endothelial cells (BECs) reveals transcriptional segmentation into distinct arterial, capillary, and venous identities with age and experimental interventionsChanges with age are heterogenous across vessel segments, with aged capillaries enriched in signatures of innate immunity, TGF-β and VEGF signaling, hypoxia and oxidative stressBECs sense and respond transcriptionally to diverse circulatory cues: inflammatory, proaging, or rejuvenatingAged plasma exposure recapitulates—and young plasma reverses—key transcriptomic signatures of normal BEC agingBEC response to aged and young plasma reveals cell non-autonomous mechanisms of blood-brain-barrier aging

The regulatory mechanisms by which neurons coordinate their physiology and connectivity are not well understood. The Drosophila olfactory receptor neurons (ORNs) provide an excellent system to investigate this question. Each ORN type expresses a unique olfactory receptor or a combination thereof, and sends their axons to a stereotyped glomerulus. Using single-cell RNA-sequencing, we identified 33 transcriptomic clusters for ORNs and mapped 20 to their glomerular types, demonstrating that transcriptomic clusters correspond well with anatomically and physiologically defined ORN types. Each ORN type expresses hundreds of transcription factors. Transcriptome-instructed genetic analyses revealed that 1) one broadly expressed transcription factor (Acj6) only regulates olfactory receptor expression in one ORN type and only wiring specificity in another type; 2) one type-restricted transcription factor (Forkhead) only regulates receptor expression; and 3) another type-restricted transcription factor (Unplugged) regulates both events. Thus, ORNs utilize diverse strategies and complex regulatory networks to coordinate their physiology and connectivity.

Molecular interactions at the cellular interface mediate organized assembly of single cells into tissues, and thus govern the development and physiology of multicellular organisms. Here, we developed a cell-type-specific, spatiotemporally-resolved approach to profile cell-surface proteomes in intact tissues. Quantitative profiling of cell-surface proteomes of Drosophila olfactory projection neurons (PNs) in pupae and adults revealed a global down-regulation of wiring molecules and an up-regulation of synaptic molecules in the transition from developing to mature PNs. A proteome-instructed in vivo screen identified 20 new cell-surface molecules regulating neural circuit assembly, many of which belong to evolutionarily conserved protein families not previously linked to neural development. Genetic analysis further revealed that the lipoprotein receptor LRP1 cell-autonomously controls PN dendrite targeting, contributing to the formation of a precise olfactory map. These findings highlight the power of temporally-resolved in situ cell-surface proteomic profiling in discovering new regulators of brain wiring.

This work presents the workflow for generating single cell transcriptomes derived from primary human uterine and cervical tissue obtained during planned cesarean hysterectomies. In total, a catalogue of 310 single cell transcriptomes are obtained, cell types present in these biopsies are inferred, and specific genes defining each of the cellular types present in the tissue are identified. Further validation of the inferred cell identity is also demonstrated via meta-analysis of independent repositories in literature generated by bulk sequenced data of fluorescence-activated cell sorting sorted cells.

Significance A fraction of the 400 million people infected with dengue annually progresses to severe dengue (SD). Yet, there are currently no biomarkers to predict disease progression. We profiled the landscape of host transcripts and viral RNA in thousands of single blood cells from dengue patients prior to progressing to SD. We discovered cell type-specific immune activation and candidate predictive biomarkers. We also determined preferential virus association with specific cell populations, particularly naive B cells and monocytes. We explored immune activation of bystander cells, clonality and somatic evolution of adaptive immune repertoires, as well as viral genomics. This multifaceted approach could advance understanding of pathogenesis of any viral infection, map an atlas of infected cells, and promote the development of prognostics., Dengue virus (DENV) infection can result in severe complications. However, the understanding of the molecular correlates of severity is limited, partly due to difficulties in defining the peripheral blood mononuclear cells (PBMCs) that contain DENV RNA in vivo. Accordingly, there are currently no biomarkers predictive of progression to severe dengue (SD). Bulk transcriptomics data are difficult to interpret because blood consists of multiple cell types that may react differently to infection. Here, we applied virus-inclusive single-cell RNA-seq approach (viscRNA-Seq) to profile transcriptomes of thousands of single PBMCs derived early in the course of disease from six dengue patients and four healthy controls and to characterize distinct leukocyte subtypes that harbor viral RNA (vRNA). Multiple IFN response genes, particularly MX2 in naive B cells and CD163 in CD14+ CD16+ monocytes, were up-regulated in a cell-specific manner before progression to SD. The majority of vRNA-containing cells in the blood of two patients who progressed to SD were naive IgM B cells expressing the CD69 and CXCR4 receptors and various antiviral genes, followed by monocytes. Bystander, non-vRNA–containing B cells also demonstrated immune activation, and IgG1 plasmablasts from two patients exhibited clonal expansions. Lastly, assembly of the DENV genome sequence revealed diversity at unexpected sites. This study presents a multifaceted molecular elucidation of natural dengue infection in humans with implications for any tissue and viral infection and proposes candidate biomarkers for prediction of SD.

Venezuelan Equine Encephalitis Virus (VEEV) is a major biothreat agent that naturally causes outbreaks in humans and horses particularly in tropical areas of the western hemisphere, for which no antiviral therapy is currently available. The host response to VEEV and the cellular factors this alphavirus hijacks to support its effective replication or evade cellular immune responses are largely uncharacterized. We have previously demonstrated tremendous cell-to-cell heterogeneity in viral RNA (vRNA) and cellular transcript levels during flaviviral infection using a novel virus-inclusive single-cell RNA-Seq approach. Here, we used this unbiased, genome-wide approach to simultaneously profile the host transcriptome and vRNA in thousands of single cells during infection of human astrocytes with the live-attenuated vaccine strain of VEEV (TC-83). Host transcription was profoundly suppressed, yet “superproducer cells” with extremely high vRNA abundance emerged during the first viral life cycle and demonstrated an altered transcriptome relative to both uninfected cells and cells with high vRNA abundance harvested at later time points. Additionally, cells with increased structural-to-nonstructural transcript ratio exhibited upregulation of intracellular membrane trafficking genes at later time points. Loss- and gain-of-function experiments confirmed pro- and antiviral activities in both vaccine and virulent VEEV infections among the products of transcripts that positively or negatively correlated with vRNA abundance, respectively. Lastly, comparison with single cell transcriptomic data from other viruses highlighted common and unique pathways perturbed by infection across evolutionary scales. This study provides a high-resolution characterization of the VEEV (TC-83)-host interplay, identifies candidate targets for antivirals, and establishes a comparative single-cell approach to study the evolution of virus-host interactions., Author summary Little is known about the host response to Venezuelan Equine Encephalitis Virus (VEEV) and the cellular factors this alphavirus hijacks to support effective replication or evade cellular immune responses. Monitoring dynamics of host and viral RNA (vRNA) during viral infection at a single-cell level can provide insight into the virus-host interplay at a high resolution. Here, a single-cell RNA sequencing technology that detects host and viral RNA was used to investigate the interactions between TC-83, the vaccine strain of VEEV, and the human host during the course of infection of U-87 MG cells (human astrocytoma). Virus abundance and host transcriptome were heterogeneous across cells from the same culture. Subsets of differentially expressed genes, positively or negatively correlating with vRNA abundance, were identified and subsequently in vitro validated as candidate proviral and antiviral factors, respectively, in TC-83 and/or virulent VEEV infections. In the first replication cycle, “superproducer” cells exhibited rapid increase in vRNA abundance and unique gene expression patterns. At later time points, cells with increased structural-to-nonstructural transcript ratio demonstrated upregulation of intracellular membrane trafficking genes. Lastly, comparing the VEEV dataset with published datasets on other RNA viruses revealed unique and overlapping responses across viral clades. Overall, this study improves the understanding of VEEV-host interactions, reveals candidate targets for antiviral approaches, and establishes a comparative single-cell approach to study the evolution of virus-host interactions.

Somatic mosaicism, the occurrence and propagation of genetic variation in cell lineages after fertilization, is increasingly recognized to play a causal role in a variety of human diseases. We investigated the case of life-threatening arrhythmia in a 10-day-old infant with long QT syndrome (LQTS). Rapid genome sequencing suggested a variant in the sodium channel NaV1.5 encoded by SCN5A, NM_000335:c.5284G > T predicting p.(V1762L), but read depth was insufficient to be diagnostic. Exome sequencing of the trio confirmed read ratios inconsistent with Mendelian inheritance only in the proband. Genotyping of single circulating leukocytes demonstrated the mutation in the genomes of 8% of patient cells, and RNA sequencing of cardiac tissue from the infant confirmed the expression of the mutant allele at mosaic ratios. Heterologous expression of the mutant channel revealed significantly delayed sodium current with a dominant negative effect. To investigate the mechanism by which mosaicism might cause arrhythmia, we built a finite element simulation model incorporating Purkinje fiber activation. This model confirmed the pathogenic consequences of cardiac cellular mosaicism and, under the presenting conditions of this case, recapitulated 2:1 AV block and arrhythmia. To investigate the extent to which mosaicism might explain undiagnosed arrhythmia, we studied 7,500 affected probands undergoing commercial gene-panel testing. Four individuals with pathogenic variants arising from early somatic mutation events were found. Here we establish cardiac mosaicism as a causal mechanism for LQTS and present methods by which the general phenomenon, likely to be relevant for all genetic diseases, can be detected through single-cell analysis and next-generation sequencing.

SummaryThe striatum contributes to many cognitive processes and disorders, but its cell types are incompletely characterized. We show that microfluidic and FACS-based single-cell RNA sequencing of mouse striatum provides a well-resolved classification of striatal cell type diversity. Transcriptome analysis revealed ten differentiated, distinct cell types, including neurons, astrocytes, oligodendrocytes, ependymal, immune, and vascular cells, and enabled the discovery of numerous marker genes. Furthermore, we identified two discrete subtypes of medium spiny neurons (MSNs) that have specific markers and that overexpress genes linked to cognitive disorders and addiction. We also describe continuous cellular identities, which increase heterogeneity within discrete cell types. Finally, we identified cell type-specific transcription and splicing factors that shape cellular identities by regulating splicing and expression patterns. Our findings suggest that functional diversity within a complex tissue arises from a small number of discrete cell types, which can exist in a continuous spectrum of functional states.

Direct lineage reprogramming represents a remarkable conversion of cellular and transcriptome states. However, the intermediate stages through which individual cells progress during reprogramming are largely undefined. Here we use single-cell RNA sequencing at multiple time points to dissect direct reprogramming from mouse embryonic fibroblasts to induced neuronal cells. By deconstructing heterogeneity at each time point and ordering cells by transcriptome similarity, we find that the molecular reprogramming path is remarkably continuous. Overexpression of the proneural pioneer factor Ascl1 results in a well-defined initialization, causing cells to exit the cell cycle and re-focus gene expression through distinct neural transcription factors. The initial transcriptional response is relatively homogeneous among fibroblasts, suggesting that the early steps are not limiting for productive reprogramming. Instead, the later emergence of a competing myogenic program and variable transgene dynamics over time appear to be the major efficiency limits of direct reprogramming. Moreover, a transcriptional state, distinct from donor and target cell programs, is transiently induced in cells undergoing productive reprogramming. Our data provide a high-resolution approach for understanding transcriptome states during lineage differentiation.

In a primitive chordate model of natural chimerism, one chimeric partner is often eliminated in a process of allogeneic resorption. Here, we identify the cellular framework underlying loss of tolerance to one partner within a natural Botryllus schlosseri chimera. We show that the principal cell type mediating chimeric partner elimination is a cytotoxic morula cell (MC). Proinflammatory, developmental cell death programs render MCs cytotoxic and, in collaboration with activated phagocytes, eliminate chimeric partners during the "takeover" phase of blastogenic development. Among these genes, the proinflammatory cytokine IL-17 enhances cytotoxicity in allorecognition assays. Cellular transfer of FACS-purified MCs from allogeneic donors into recipients shows that the resorption response can be adoptively acquired. Transfer of 1 × 10(5) allogeneic MCs eliminated 33 of 78 (42%) recipient primary buds and 20 of 76 (20.5%) adult parental adult organisms (zooids) by 14 d whereas transfer of allogeneic cell populations lacking MCs had only minimal effects on recipient colonies. Furthermore, reactivity of transferred cells coincided with the onset of developmental-regulated cell death programs and disproportionately affected developing tissues within a chimera. Among chimeric partner "losers," severe developmental defects were observed in asexually propagating tissues, reflecting a pathologic switch in gene expression in developmental programs. These studies provide evidence that elimination of one partner in a chimera is an immune cell-based rejection that operates within histocompatible pairs and that maximal allogeneic responses involve the coordination of both phagocytic programs and the "arming" of cytotoxic cells.

Cerebral spinal fluid (CSF) from brain tumor patients contains tumor cellular and cell-free DNA (cfDNA), which provides a less-invasive and routinely accessible method to obtain tumor genomic information. In this report, we used droplet digital PCR to test mutant tumor DNA in CSF of a patient to monitor the treatment response of metastatic melanoma leptomeningeal disease (LMD). The primary melanoma was known to have a BRAF (V600E) mutation, and the patient was treated with whole brain radiotherapy and BRAF inhibitors. We collected 9 CSF samples over 6 months. The mutant cfDNA fraction gradually decreased from 53 % (time of diagnosis) to 0 (time of symptom alleviation) over the first 6 time points. Three months after clinical improvement, the patient returned with severe symptoms and the mutant cfDNA was again detected in CSF at high levels. The mutant DNA fraction corresponded well with the patient's clinical response. We used whole exome sequencing to examine the mutation profiles of the LMD tumor DNA in CSF before therapeutic response and after disease relapse, and discovered a canonical cancer mutation PTEN (R130*) at both time points. The cellular and cfDNA revealed similar mutation profiles, suggesting cfDNA is representative of LMD cells. This study demonstrates the potential of using cellular or cfDNA in CSF to monitor treatment response for LMD.

Nature Biotechnology asks a selection of researchers about the most exciting frontier in their field and the most needed technologies for advancing knowledge and applications.

The vast majority of microbial life remains uncatalogued due to the inability to cultivate these organisms in the laboratory. This “microbial dark matter” represents a substantial portion of the tree of life and of the populations that contribute to chemical cycling in many ecosystems. In this work, we leveraged an existing single-cell genomic data set representing the candidate bacterial phylum “ Calescamantes ” (EM19) to calibrate machine learning algorithms and define metagenomic bins directly from pyrosequencing reads derived from Great Boiling Spring in the U.S. Great Basin. Compared to other assembly-based methods, taxonomic binning with a read-based machine learning approach yielded final assemblies with the highest predicted genome completeness of any method tested. Read-first binning subsequently was used to extract Calescamantes bins from all metagenomes with abundant Calescamantes populations, including metagenomes from Octopus Spring and Bison Pool in Yellowstone National Park and Gongxiaoshe Spring in Yunnan Province, China. Metabolic reconstruction suggests that Calescamantes are heterotrophic, facultative anaerobes, which can utilize oxidized nitrogen sources as terminal electron acceptors for respiration in the absence of oxygen and use proteins as their primary carbon source. Despite their phylogenetic divergence, the geographically separate Calescamantes populations were highly similar in their predicted metabolic capabilities and core gene content, respiring O 2 , or oxidized nitrogen species for energy conservation in distant but chemically similar hot springs.

The field of single-cell genomics is advancing rapidly and is generating many new insights into complex biological systems, ranging from the diversity of microbial ecosystems to the genomics of human cancer. In this Review, we provide an overview of the current state of the field of single-cell genome sequencing. First, we focus on the technical challenges of making measurements that start from a single molecule of DNA, and then explore how some of these recent methodological advancements have enabled the discovery of unexpected new biology. Areas highlighted include the application of single-cell genomics to interrogate microbial dark matter and to evaluate the pathogenic roles of genetic mosaicism in multicellular organisms, with a focus on cancer. We then attempt to predict advances we expect to see in the next few years.

Current digital biology platforms lack the ability to perform continuous measurements for transient analysis. A fundamental challenge is to perform complex fluidic manipulation processes such as washing and mixing in individual reaction volumes. Here, we present a reusable digital biology platform where the reaction compartmentalization and commencement are controlled by micromechanical valves fabricated in high density through microfluidic very large scale integration (mVLSI) technology. Background noise correction enabled by the platform improves signal-to-noise ratio and thus eliminates the need for sophisticated imaging technologies. We have used the detection platform for probing single molecules of the β-galactosidase enzyme. The measurements were repeated hundreds of times at concentrations as low as 0.8 fM (resulting in a theoretical detection limit of 3 aM). We have also demonstrated multiple TNF-α measurements with a magnetic bead based digital ELISA assay. The assay showed that the average number of enzymes per magnetic bead is 0.55 for 10 pM TNF-α (compared to 0.08 for negative control). This automated and reusable digital platform allows on-chip assay preparation and continuous measurements; as a result, it will enable single cell/enzyme studies and clinical diagnostic tests (i.e. digital ELISA) to be performed in shorter time scales and with lower detection limits.

In-depth phenotyping of human intestinal antibody secreting cells (ASCs) and their precursors is important for developing improved mucosal vaccines. We used single-cell mass cytometry to simultaneously analyze 34 differentiation and trafficking markers on intestinal and circulating B cells. In addition, we labeled rotavirus (RV) double-layered particles with a metal isotope and characterized B cells specific to the RV VP6 major structural protein. We describe the heterogeneity of the intestinal B-cell compartment, dominated by ASCs with some phenotypic and transcriptional characteristics of long-lived plasma cells. Using principal component analysis, we visualized the phenotypic relationships between major B-cell subsets in the intestine and blood, and revealed that IgM(+) memory B cells (MBCs) and naive B cells were phenotypically related as were CD27(-) MBCs and switched MBCs. ASCs in the intestine and blood were highly clonally related, but associated with distinct trajectories of phenotypic development. VP6-specific B cells were present among diverse B-cell subsets in immune donors, including naive B cells, with phenotypes representative of the overall B-cell pool. These data provide a high dimensional view of intestinal B cells and the determinants regulating humoral memory to a ubiquitous, mucosal pathogen at steady-state.

Objective: ACE2 enzyme serves as SARS-CoV-2 human receptor through binding of the viral S protein and subsequent trimming of S protein between S1 and S2 units by host serine proteases as TMPRSS2, CTSB or CTSL. Here, we aim to investigate the expression of the different cell entry proteins involved in SARS-CoV-2 infection in the different cell types of the human endometrium throughout the menstrual cycle using single-cell RNAseq (scRNAseq). Design: Gene expression patterns for SARS-CoV-2 entry molecules were analyzed by scRNAseq in a total of 73,181 endometrial single cells obtained from endometrial biopsies from 19 reproductive-age women across the full menstrual cycle (Fluidigm C1: 2,149 cells) and 10 women from the same cohort (10x: 71,032 cells). For two women, both C1 and 10x data were collected as anchors for comparison. Materials and Methods: After tissue dissociation, single cell capture was performed on Fluidigm C1 system (n= 2,149 cells) or Chromium 10x system (Chromium Next GEM Chip G, 10x Genomics) (n= 71,032 cells) followed by reverse-transcription, cDNA generation and library construction. Barcoded libraries were sequenced in pair-end reads on Nextseq (Illumina) for the C1 dataset or Novaseq (Illumina) for the 10x dataset. Data pre-processing, quality filtering, and statistical analyses were performed using custom Python, R, and Java scripts. Results: Expression analysis across the menstrual cycle showed no significant expression of ACE2 in stromal or unciliated epithelial cells in any cycle phase. TMPRSS2 was expressed more highly in glandular epithelial cells during the early proliferative phase and towards the end of the cycle. Interestingly, expression of CTSB and CTSL was observed in both stromal and epithelial cells across all phases of the menstrual cycle, with CTSB the more abundant of the two. All four genes were simultaneously expressed in less than 0.7% of glandular epithelial cells. Expression analysis during the secretory phase did not detect significant expression of ACE2 (less than 2% of epithelial or stromal cells). TMPRSS2 showed mild expression in about 12% of unciliated epithelial cells. In contrast, CTSB and CTSLwere highly expressed in ∼80% and ∼40% of cells during the mid-late secretory phase, aligning with what we detected via Fluidigm. In addition, while CTSB was highly expressed in both epithelial and stromal cells, CTSL was more highly expressed in stromal cells across the menstrual cycle. Conclusions: Percentages of endometrial cells expressing ACE2, TMPRSS2, CTSB, or CTSL were

Background: Our feasibility study employed a novel genomic detection methodology that enriches 5-hydroxymethylcytosine (5hmC) loci in cell free DNA (cfDNA) from the plasma of cancer patients using click chemistry coupled with sequencing and machine learning based classification methods. These classification methods were developed to detect the presence of disease in the plasma of cancer and control subjects. Cancer and control patient cfDNA cohorts were accrued from multiple sites consisting of 48 breast, 55 lung, 32 prostate and 2 pancreatic datasets consisting of 41 and 53 cancer subjects (Set 1 and 2). In addition, a control cohort of 260 subjects (non-cancer) was employed to match cancer patient demographics (age, sex and smoking status) in a case-control study design. Methods: Machine learning methods, applied to each cancer case cohort individually, with a balancing non-cancer cohort, were able to classify cancer and control samples. Measures of predictive performance using 5-fold cross validation coupled with out-of-fold Area Under the Receiver Operating Characteristic Curve (AUROC) measures were employed. Gene sets selected as part of biomarker discovery were further analyzed for disease relevance using pathway analysis tools (GSEA, mSigDB). Results: 260 controls and 229 cancers from four disease types (breast, lung, pancreas and prostate) were analyzed; more than 60% of cancer patients had early stage disease (I or II). Predictive performance, employing AUROC measures, was established for breast (0.89), lung (0.84), pancreas (set 1 - 0.95 and 2 - 0.93) and prostate (0.83). The genes defining each of these predictive models were enriched for pathways relevant to disease specific etiology, notably in the control of gene regulation in these same pathways. The breast cancer cohort consisted primarily of stage I and II patients including tumors < 2 cm and these samples exhibited a higher prediction probability score. The prostate cancer cohort consisted of both indolent and aggressive disease sample and prediction performance was equally high for both (AUROC for indolent vs aggressive was 0.81 and 0.77, respectively). Conclusions: These findings suggest that 5hmC changes in cfDNA enable non-invasive detection of early stage breast, pancreatic, prostate, and lung cancers. Furthermore, 5hmC profiling in cfDNA may enable the prediction of clinically relevant features such as tumor size in breast adenocarcinoma or indolent disease in prostate cancer. Finally, this study identified a suite of 5hmC biomarkers that may be further validated in larger, and more diverse, patient cohorts. Citation Format: Anna Bergamaschi, Jeremy Ku, Yuhong Ning, Francois Collin, Chris Ellison, Tierney Phillips, Erin McCarthy, Wendy Wang, Michael Antoine, David Haan, Aaron Scott, Paul Lloyd, Gulfem Guler, Alan Ashworth, Stephen Quake, Samuel Levy. Epigenomic detection of multiple cancers in plasma derived cell free DNA [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 783.

SMAD pathways govern epithelial proliferation, and transforming growth factor β (TGF-β and BMP signaling through SMAD members has distinct effects on mammary development and homeostasis. Here, we show that LEFTY1, a secreted inhibitor of NODAL/SMAD2 signaling, is produced by mammary progenitor cells and, concomitantly, suppresses SMAD2 and SMAD5 signaling to promote long-term proliferation of normal and malignant mammary epithelial cells. In contrast, BMP7, a NODAL antagonist with context-dependent functions, is produced by basal cells and restrains progenitor cell proliferation. In normal mouse epithelium, LEFTY1 expression in a subset of luminal cells and rare basal cells opposes BMP7 to promote ductal branching. LEFTY1 binds BMPR2 to suppress BMP7-induced activation of SMAD5, and this LEFTY1-BMPR2 interaction is specific to tumor-initiating cells in triple-negative breast cancer xenografts that rely on LEFTY1 for growth. These results suggest that LEFTY1 is an endogenous dual-SMAD inhibitor and that suppressing its function may represent a therapeutic vulnerability in breast cancer.

Summary There is pressing urgency to understand the pathogenesis of the severe acute respiratory syndrome coronavirus clade 2 (SARS-CoV-2), which causes the disease COVID-19. SARS-CoV-2 spike (S) protein binds angiotensin-converting enzyme 2 (ACE2), and in concert with host proteases, principally transmembrane serine protease 2 (TMPRSS2), promotes cellular entry. The cell subsets targeted by SARS-CoV-2 in host tissues and the factors that regulate ACE2 expression remain unknown. Here, we leverage human, non-human primate, and mouse single-cell RNA-sequencing (scRNA-seq) datasets across health and disease to uncover putative targets of SARS-CoV-2 among tissue-resident cell subsets. We identify ACE2 and TMPRSS2 co-expressing cells within lung type II pneumocytes, ileal absorptive enterocytes, and nasal goblet secretory cells. Strikingly, we discovered that ACE2 is a human interferon-stimulated gene (ISG) in vitro using airway epithelial cells and extend our findings to in vivo viral infections. Our data suggest that SARS-CoV-2 could exploit species-specific interferon-driven upregulation of ACE2, a tissue-protective mediator during lung injury, to enhance infection., Graphical Abstract, Highlights • Meta-analysis of human, non-human primate, and mouse single-cell RNA-seq datasets for putative SARS-CoV-2 targets • Type II pneumocytes, nasal secretory cells, and absorptive enterocytes are ACE2 + TMPRSS2 + • Interferon and influenza increase ACE2 in human nasal epithelia and lung tissue • Mouse Ace2 is not upregulated by interferon, raising implications for disease modeling, Analysis of single-cell RNA-seq datasets from human, non-human primate, and mouse barrier tissues identifies putative cellular targets of SARS-CoV-2 on the basis of ACE2 and TMPRSS2 expression. ACE2 represents a previously unappreciated interferon-stimulated gene in human, but not mouse, epithelial tissues, identifying anti-viral induction of a host tissue-protective mechanism, but also a potential means for viral exploitation of the host response.

The mammalian striatum is involved in many complex behaviors and yet is composed largely of a single neuron class: the spiny projection neuron (SPN). It is unclear to what extent the functional specialization of the striatum is due to the molecular specialization of SPN subtypes. We sought to define the molecular and anatomical diversity of adult SPNs using single-cell RNA sequencing (scRNA-seq) and quantitative RNA in situ hybridization (ISH). We computationally distinguished discrete versus continuous heterogeneity in scRNA-seq data and found that SPNs in the striatum can be classified into four major discrete types with no implied spatial relationship between them. Within these discrete types, we find continuous heterogeneity encoding spatial gradients of gene expression and defining anatomical location in a combinatorial mechanism. Our results suggest that neuronal circuitry has a substructure at far higher resolution than is typically interrogated, which is defined by the precise identity and location of a neuron.

SUMMARY Antibody memory protects humans from many diseases. Protective antibody memory responses require activation of transcriptional programs, cell proliferation, and production of antigen-specific antibodies, but how these aspects of the response are coordinated is poorly understood. We profile the molecular and cellular features of the antibody response to influenza vaccination by integrating single-cell transcriptomics, longitudinal antibody repertoire sequencing, and antibody binding measurements. Single-cell transcriptional profiling reveals a program of memory B cell activation characterized by CD11c and T-bet expression associated with clonal expansion and differentiation toward effector function. Vaccination elicits an antibody clone, which rapidly acquired broad high-affinity hemagglutinin binding during affinity maturation. Unexpectedly, many antibody clones elicited by vaccination do not bind vaccine, demonstrating non-specific activation of bystander antibodies by influenza vaccination. These results offer insight into how molecular recognition, transcriptional programs, and clonal proliferation are coordinated in the human B cell repertoire during memory recall., Graphical Abstract, In Brief Antibody memory requires coordination of molecular recognition, gene expression programs, and clonal dynamics. Horns et al. study the human antibody memory response using single-cell and repertoire sequencing, revealing a transcriptional program of memory B cell activation, broadly binding anti-influenza antibodies, and widespread bystander activation of non-vaccine-binding antibodies after influenza vaccination.

Summary Molecular interactions at the cellular interface mediate organized assembly of single cells into tissues and, thus, govern the development and physiology of multicellular organisms. Here, we developed a cell-type-specific, spatiotemporally resolved approach to profile cell-surface proteomes in intact tissues. Quantitative profiling of cell-surface proteomes of Drosophila olfactory projection neurons (PNs) in pupae and adults revealed global downregulation of wiring molecules and upregulation of synaptic molecules in the transition from developing to mature PNs. A proteome-instructed in vivo screen identified 20 cell-surface molecules regulating neural circuit assembly, many of which belong to evolutionarily conserved protein families not previously linked to neural development. Genetic analysis further revealed that the lipoprotein receptor LRP1 cell-autonomously controls PN dendrite targeting, contributing to the formation of a precise olfactory map. These findings highlight the power of temporally resolved in situ cell-surface proteomic profiling in discovering regulators of brain wiring.

Background Deer antlers are bony structures that re-grow at very high rates, making them an attractive model for studying rapid bone regeneration. Methods To identify the genes that are involved in this fast pace of bone growth, an in vitro RNA-seq model that paralleled the sharp differences in bone growth between deer antlers and humans was established. Subsequently, RNA-seq (> 60 million reads per library) was used to compare transcriptomic profiles. Uniquely expressed deer antler proliferation as well as mineralization genes were identified via a combination of differential gene expression and subtraction analysis. Thereafter, the physiological relevance as well as contributions of these identified genes were determined by immunofluorescence, gene overexpression, and gene knockdown studies. Results Cell characterization studies showed that in vitro-cultured deer antler-derived reserve mesenchyme (RM) cells exhibited high osteogenic capabilities and cell surface markers similar to in vivo counterparts. Under identical culture conditions, deer antler RM cells proliferated faster (8.6–11.7-fold increase in cell numbers) and exhibited increased osteogenic differentiation (17.4-fold increase in calcium mineralization) compared to human mesenchymal stem cells (hMSCs), paralleling in vivo conditions. Comparative RNA-seq identified 40 and 91 previously unknown and uniquely expressed fallow deer (FD) proliferation and mineralization genes, respectively, including uhrf1 and s100a10. Immunofluorescence studies showed that uhrf1 and s100a10 were expressed in regenerating deer antlers while gene overexpression and gene knockdown studies demonstrated the proliferation contributions of uhrf1 and mineralization capabilities of s100a10. Conclusion Using a simple, in vitro comparative RNA-seq approach, novel genes pertinent to fast bony antler regeneration were identified and their proliferative/osteogenic function was verified via gene overexpression, knockdown, and immunostaining. This combinatorial approach may be applicable to discover unique gene contributions between any two organisms for a given phenomenon-of-interest. Electronic supplementary material The online version of this article (10.1186/s13287-018-1027-6) contains supplementary material, which is available to authorized users.

SUMMARYSmall non-coding RNAs (ncRNAs) play a vital role in a broad range of biological processes both in health and disease. A comprehensive quantitative reference of small ncRNA expression would significantly advance our understanding of ncRNA roles in shaping tissue functions. Here, we systematically profiled the levels of five ncRNA classes (miRNA, snoRNA, snRNA, scaRNA and tRNA fragments) across eleven mouse tissues by deep sequencing. Using fourteen biological replicates spanning both sexes, we identified that ~ 30% of small ncRNAs are distributed across the body in a tissue-specific manner with some are also being sexually dimorphic. We found that miRNAs are subject to “arm switching” between healthy tissues and that tRNA fragments are retained within tissues in both a gene- and a tissue-specific manner. Out of eleven profiled tissues we confirmed that brain contains the largest number of unique small ncRNA transcripts, some of which were previously annotated while others are identified for the first time in this study. Furthermore, by combining these findings with single-cell ATAC-seq data, we were able to connect identified brain-specific ncRNA with their cell types of origin. These results yield the most comprehensive characterization of specific and ubiquitous small RNAs in individual murine tissues to date, and we expect that this data will be a resource for the further identification of ncRNAs involved in tissue-function in health and dysfunction in disease.HIGHLIGHTS-An atlas of tissue levels of multiple small ncRNA classes generated from 14 biological replicates of both sexes across 11 tissues-Distinct distribution patterns of miRNA arms and tRNA fragments across tissues suggest the existence of tissue-specific mechanisms of ncRNA cleavage and retention-miRNA expression is sex specific in healthy tissues-Small RNA-seq and scATAC-seq data integration produce a detailed map of cell-type specific ncRNA profiles in the mouse brain

Pancreatic cancers are typically diagnosed at late stage where disease prognosis is poor as exemplified by a 5-year survival rate of 8.2%. Earlier diagnosis would be beneficial by enabling surgical resection or earlier application of therapeutic regimens. We investigated the detection of pancreatic ductal adenocarcinoma (PDAC) in a non-invasive manner by interrogating changes in 5-hydroxymethylation cytosine status (5hmC) of circulating cell free DNA in the plasma of a PDAC cohort (n=51) in comparison with a non-cancer cohort (n=41). We found that 5hmC sites are enriched in a disease and stage specific manner in exons, 3’UTRs and transcription termination sites. Our data show that 5hmC density is reduced in promoters and histone H3K4me3-associated sites with progressive disease suggesting increased transcriptional activity. 5hmC density is differentially represented in thousands of genes, and a stringently filtered set of the most significant genes points to biology related to pancreas (GATA4, GATA6, PROX1, ONECUT1) and/or cancer development (YAP1, TEAD1, PROX1, ONECUT1, ONECUT2, IGF1 and IGF2). Regularized regression models were built using 5hmC densities in statistically filtered genes or a comprehensive set of highly variable 5hmC counts in genes and performed with an AUC = 0.94-0.96 on training data. We were able to test the ability to classify PDAC and non-cancer samples with the Elastic net and Lasso models on two external pancreatic cancer 5hmC data sets and found validation performance to be AUC = 0.74-0.97. The findings suggest that 5hmC changes enable classification of PDAC patients with high fidelity and are worthy of further investigation on larger cohorts of patient samples.

SummaryMicroglia are increasingly recognized for their major contributions during brain development and neurodegenerative disease. It is currently unknown if these functions are carried out by subsets of microglia during different stages of development and adulthood or within specific brain regions. Here, we performed deep single-cell RNA sequencing (scRNA-seq) of microglia and related myeloid cells sorted from various regions of embryonic, postnatal, and adult mouse brains. We found that the majority of adult microglia with homeostatic signatures are remarkably similar in transcriptomes, regardless of brain region. By contrast, postnatal microglia represent a more heterogeneous population. We discovered that postnatal white matter-associated microglia (WAM) are strikingly different from microglia in other regions and express genes enriched in degenerative disease-associated microglia. These postnatal WAM have distinct amoeboid morphology, are metabolically active, and phagocytose newly formed oligodendrocytes. This scRNA-seq atlas will be a valuable resource for dissecting innate immune functions in health and disease.HighlightsMyeloid scRNA-seq atlas across brain regions and developmental stagesLimited transcriptomic heterogeneity of homeostatic microglia in the adult brainPhase-specific gene sets of proliferating microglia along cell cycle pseudotimePhagocytic postnatal white matter-associated microglia sharing DAM gene signatures

Recently high-throughput image-based transcriptomic methods were developed and enabled researchers to spatially resolve gene expression variation at the molecular level for the first time. In this work, we develop a general analysis tool to quantitatively study the spatial correlations of gene expression in fixed tissue sections. As an illustration, we analyze the spatial distribution of single mRNA molecules measured by in situ sequencing on human fetal pancreas at three developmental time points–80, 87 and 117 days post-fertilization. We develop a density profile-based method to capture the spatial relationship between gene expression and other morphological features of the tissue sample such as position of nuclei and endocrine cells of the pancreas. In addition, we build a statistical model to characterize correlations in the spatial distribution of the expression level among different genes. This model enables us to infer the inhibitory and clustering effects throughout different time points. Our analysis framework is applicable to a wide variety of spatially-resolved transcriptomic data to derive biological insights.

Dengue virus (DENV) infection can result in severe complications. Yet, the understanding of the molecular correlates of severity is limited, partly due to difficulties in defining the peripheral blood mononuclear cells (PBMCs) that are associated with DENV in vivo. Additionally, there are currently no biomarkers predictive of progression to severe dengue (SD). Bulk transcriptomics data are difficult to interpret because blood consists of multiple cell types that may react differently to infection. Here we applied virus-inclusive single cell RNA-seq approach (viscRNA-Seq) to profile transcriptomes of thousands of single PBMCs derived early in the course of disease from six dengue patients and four healthy controls, and to characterize distinct DENV-associated leukocytes. Multiple genes, particularly interferon response genes, were upregulated in a cell-specific manner prior to progression to SD. Expression of MX2 in naive B cells and CD163 in CD14+ CD16+ monocytes was predictive of SD. The majority of DENV-associated cells in the blood of two patients who progressed to SD were naive IgM B cells expressing the CD69 and CXCR4 receptors and antiviral genes, followed by monocytes. Bystander uninfected B cells also demonstrated immune activation, and plasmablasts from two patients exhibited antibody lineages with convergently hypermutated heavy chain sequences. Lastly, assembly of the DENV genome revealed diversity at unexpected genomic sites. This study presents a multi-faceted molecular elucidation of natural dengue infection in humans and proposes biomarkers for prediction of SD, with implications for profiling any tissue and viral infection, and for the development of a dengue prognostic assay.SignificanceA fraction of the 400 million people infected with dengue annually progresses to severe dengue (SD). Yet, there are currently no biomarkers to effectively predict disease progression. We profiled the landscape of host transcripts and viral RNA in thousands of single blood cells from dengue patients prior to progressing to SD. We discovered cell-type specific immune activation and candidate predictive biomarkers. We also revealed preferential virus association with specific cell populations, particularly naive B cells and monocytes. We then explored immune activation of bystander cells, clonality and somatic evolution of adaptive immune repertoires, and viral genomics. This multi-faceted approach could advance understanding of pathogenesis of any viral infection, map an atlas of infected cells and promote the development of prognostics.

SummaryIn a human menstrual cycle, the endometrium undergoes remodeling, shedding, and regeneration, all of which are driven by substantial gene expression changes in the underlying cellular hierarchy. Despite its importance in human fertility and regenerative biology, mechanistic understanding of this unique type of tissue homeostasis remains rudimentary. We characterized the transcriptomic transformation of human endometrium at single cell resolution, dissecting the multidimensional cellular heterogeneity of this tissue across the entire natural menstrual cycle. We profiled the behavior of 6 endometrial cell types, including a previously uncharacterized ciliated epithelial cell type, during four major phases of endometrial transformation, and found characteristic signatures for each cell type and phase. We discovered that human window of implantation opens with an abrupt and discontinuous transcriptomic activation in the epithelia, accompanied with widespread decidualized feature in the stromal fibroblasts. These data reveal signatures in the luminal and glandular epithelia during epithelial gland reconstruction, and suggest a mechanism for adult gland formation.

IgE B cells unmasked Immunoglobulin E (IgE) antibodies play a central role in immune responses against helminth and protozoan parasites; however, they also contribute to allergies. IgE antibodies (and the B cells generating them) are rare and thus poorly characterized. Croote et al. performed single-cell RNA sequencing of peripheral blood B cells from patients with peanut allergies and delineated each cell's gene expression, splice variants, and antibody sequences (see the Perspective by Gould and Ramadani). Unlike other isotypes, circulating IgE B cells were mostly immature plasmablasts. Surprisingly, certain IgE antibodies manifested identical gene rearrangements in unrelated individuals. These IgE antibodies showed high affinity and unexpected cross-reactivity to peanut allergens. Science , this issue p. 1306 ; see also p. 1247

IgE antibodies provide defense against helminth infections, but can also cause life-threatening allergic reactions. Despite their importance to human health, these antibodies and the cells that produce them remain enigmatic due to their scarcity in humans; much of our knowledge of their properties is derived from model organisms. Here we describe the isolation of IgE producing B cells from the blood of individuals with food allergies, followed by a detailed study of their properties by single cell RNA sequencing (scRNA-seq). We discovered that IgE B cells are deficient in membrane immunoglobulin expression and that the IgE plasmablast state is more immature than that of other antibody producing cells. Through recombinant expression of monoclonal antibodies derived from single cells, we identified IgE antibodies which had unexpected cross-reactive specificity for major peanut allergens Ara h 2 and Ara h 3; not only are these among the highest affinity native human antibodies discovered to date, they represent a surprising example of convergent evolution in unrelated individuals who independently evolved nearly identical antibodies. Finally, we discovered that splicing within B cells of all isotypes reveals polarized germline transcription of the IgE, but not IgG4, isotype as well as several examples of biallelic expression of germline transcripts. Our results offer insights into IgE B cell transcriptomics, clonality and regulation, provide a striking example of adaptive immune convergence, and offer an approach for accelerating mechanistic disease understanding by characterizing a rare B cell population underlying IgE-mediated disease at single cell resolution.

The liver is largely composed of hepatocytes and bile duct epithelial cells (BECs). Controversy exists as to whether a liver stem/progenitor cell capable of renewing both hepatocytes and BECs exists. Single cell RNA sequencing of freshly isolated human foetal and healthy adult liver identified hepatocyte, hepatoblast and liver progenitor cell (hLPC) populations. hLPCs, found at the interface between hepatocytes and bile ducts in both foetal and adult tissue, were distinguishable from BECs by their negative expression of TROP-2. Prospective isolation followed by in vitro culture demonstrated their potential for expansion and bi-lineage differentiation. The hLPC expression signature was also conserved within expanded cell populations specific to certain cases of liver injury and cancer. These data support the idea of a true progenitor existing within healthy adult liver that can be activated upon injury. Further work to define the mechanisms regulating hLPC behaviour could advance understanding of human development and disease.

The recent advent of methods for high-throughput single-cell molecular profiling has catalyzed a growing sense in the scientific community that the time is ripe to complete the 150-year-old effort to identify all cell types in the human body. The Human Cell Atlas Project is an international collaborative effort that aims to define all human cell types in terms of distinctive molecular profiles (such as gene expression profiles) and to connect this information with classical cellular descriptions (such as location and morphology). An open comprehensive reference map of the molecular state of cells in healthy human tissues would propel the systematic study of physiological states, developmental trajectories, regulatory circuitry and interactions of cells, and also provide a framework for understanding cellular dysregulation in human disease. Here we describe the idea, its potential utility, early proofs-of-concept, and some design considerations for the Human Cell Atlas, including a commitment to open data, code, and community.

Dengue and Zika viral infections affect millions of people annually and can be complicated by hemorrhage or neurological manifestations, respectively. However, a thorough understanding of the host response to these viruses is lacking, partly because conventional approaches ignore heterogeneity in virus abundance across cells. We present viscRNA-Seq (virus-inclusive single cell RNA-Seq), an approach to probe the host transcriptome together with intracellular viral RNA at the single cell level. We applied viscRNA-Seq to monitor dengue and Zika virus infection in cultured cells and discovered extreme heterogeneity in virus abundance. We exploited this variation to identify host factors that show complex dynamics and a high degree of specificity for either virus, including proteins involved in the endoplasmic reticulum translocon, signal peptide processing, and membrane trafficking. We validated the viscRNA-Seq hits and discovered novel proviral and antiviral factors. viscRNA-Seq is a powerful approach to assess the genome-wide virus-host dynamics at single cell level.

Neuroendocrine tumors (NETs) arising from tailgut cysts are a rare but increasingly reported entity with gene expression profiles that may be indicative of the gastrointestinal cell of origin. We present a case report describing the unique pathological and genomic characteristics of a tailgut cyst NET that metastasized to liver. The histologic and immunohistochemical findings were consistent with a well-differentiated NET. Genomic testing indicates a germline frameshift in BRCA1 and a few somatic mutations of unknown significance. Transcriptomic analysis suggests an enteroendocrine L cell in the tailgut as a putative cell of origin. Genomic profiling of a rare NET and metastasis provides insight into its origin, development, and potential therapeutic options.

How a neuronal cell type is defined and how this relates to its transcriptome are still open questions. The Drosophila olfactory projection neurons (PNs) are among the best-characterized neuronal types: Different PN classes target dendrites to distinct olfactory glomeruli and PNs of the same class exhibit indistinguishable anatomical and physiological properties. Using single-cell RNA-sequencing, we comprehensively characterized the transcriptomes of 40 PN classes and unequivocally identified transcriptomes for 6 classes. We found a new lineage-specific transcription factor that instructs PN dendrite targeting. Transcriptomes of closely related PN classes exhibit the largest difference during circuit assembly, but become indistinguishable in adults, suggesting that neuronal subtype diversity peaks during development. Genes encoding transcription factors and cell-surface molecules are the most differentially expressed, indicating their central roles in specifying neuronal identity. Finally, we show that PNs use highly redundant combinatorial molecular codes to distinguish subtypes, enabling robust specification of cell identity and circuit assembly.

SummaryHematopoiesis is an essential process that evolved in multicellular animals. At the heart of this process are hematopoietic stem cells (HSCs), which are multipotent, self-renewing and generate the entire repertoire of blood and immune cells throughout life. Here we studied the hematopoietic system of Botryllus schlosseri, a colonial tunicate that has vasculature, circulating blood cells, and interesting characteristics of stem cell biology and immunity. Self-recognition between genetically compatible B. schlosseri colonies leads to the formation of natural parabionts with shared circulation, whereas incompatible colonies reject each other. Using flow-cytometry, whole-transcriptome sequencing of defined cell populations, and diverse functional assays, we identified HSCs, progenitors, immune-effector cells, the HSC niche, and demonstrated that self-recognition inhibits cytotoxic reaction. Our study implies that the HSC and myeloid lineages emerged in a common ancestor of tunicates and vertebrates and suggests that hematopoietic bone marrow and the B. schlosseri endostyle niche evolved from the same origin.

The Tabula Muris ConsortiumWe have created a compendium of single cell transcriptome data from the model organism Mus musculus comprising more than 100,000 cells from 20 organs and tissues. These data represent a new resource for cell biology, revealing gene expression in poorly characterized cell populations and allowing for direct and controlled comparison of gene expression in cell types shared between tissues, such as T-lymphocytes and endothelial cells from distinct anatomical locations. Two distinct technical approaches were used for most tissues: one approach, microfluidic droplet-based 3’-end counting, enabled the survey of thousands of cells at relatively low coverage, while the other, FACS-based full length transcript analysis, enabled characterization of cell types with high sensitivity and coverage. The cumulative data provide the foundation for an atlas of transcriptomic cell biology.

The frequent use of precautionary food allergen labeling statements such as “may contain” poses challenges to allergic individuals who rely on such labeling to determine whether a food is safe to consume. We developed a multiplexed liquid chromatography-mass spectrometry assay targeting 14 common allergens in order to survey how frequently these precautionary statements indicate allergen contamination and to assess whether variations in precautionary phrasing affect the likelihood of allergen contamination. A survey of 84 foods revealed how scheduled multiple reaction monitoring (MRM) transition interference derived from complex and heterogeneous sample matrices hinders sensitive analyte detection. As a solution, we developed MAtrix-Dependent Interference Correction (MADIC), an approach to sensitively detect trace peptide quantities through interference identification and stringent peptide quality control criteria. Applying this method, we find frequent contamination of soy in breads and corn flour, and observe additional instances of food contamination with tree nuts, wheat, milk, and egg. In some of these cases, the food had no precautionary labeling for the offending allergen. We also find that only 10% of warning labels are indicative of contamination, and that products with “same facility” precautionary labeling are not necessarily less likely to contain trace amounts of allergens than products labeled “may contain.”

From bacteria to humans, individual cells within isogenic populations can show significant variation in stress tolerance, but the nature of this heterogeneity is not clear. To investigate this, we used single-cell RNA sequencing to quantify transcript heterogeneity in single Saccharomyces cerevisiae cells treated with and without salt stress to explore population variation and identify cellular covariates that influence the stress-responsive transcriptome. Leveraging the extensive knowledge of yeast transcriptional regulation, we uncovered significant regulatory variation in individual yeast cells, both before and after stress. We also discovered that a subset of cells appears to decouple expression of ribosomal protein genes from the environmental stress response in a manner partly correlated with the cell cycle but unrelated to the yeast ultradian metabolic cycle. Live-cell imaging of cells expressing pairs of fluorescent regulators, including the transcription factor Msn2 with Dot6, Sfp1, or MAP kinase Hog1, revealed both coordinated and decoupled nucleocytoplasmic shuttling. Together with transcriptomic analysis, our results suggest that cells maintain a cellular filter against decoupled bursts of transcription factor activation but mount a stress response upon coordinated regulation, even in a subset of unstressed cells., Author summary Genetically identical cells growing in the same environment can vary in their cellular state and behavior. Such heterogeneity may explain why some cells in an isogenic population can survive sudden severe environmental stress whereas other cells succumb. Cell-to-cell variation in gene expression has been linked to variable stress survival, but how and why transcript levels vary across the transcriptome in single cells is only beginning to emerge. Here, we used single-cell RNA sequencing (scRNA-seq) to measure cell-to-cell heterogeneity in the transcriptome of budding yeast (Saccharomyces cerevisiae). We find surprising patterns of variation across known sets of transcription factor targets, indicating that cells vary in their transcriptome profile both before and after stress exposure. scRNA-seq analysis combined with live-cell imaging of transcription factor activation dynamics revealed some cells in which the stress response was coordinately activated and other cells in which the traditional response was decoupled, suggesting unrecognized regulatory nuances that expand our understanding of stress response and survival.

Human tissue in vitro models on-chip are highly desirable to dissect the complexity of a physio-pathological in vivo response because of their advantages compared to traditional static culture systems in terms of high control of microenvironmental conditions, including accurate perturbations and high temporal resolution analyses of medium outflow. Human adipose tissue (hAT) is a key player in metabolic disorders, such as Type 2 Diabetes Mellitus (T2DM). It is involved in the overall energy homeostasis not only as passive energy storage but also as an important metabolic regulator. Here, we aim at developing a large scale microfluidic platform for generating high temporal resolution of glucose uptake profiles, and consequently insulin sensitivity, under physio-pathological stimulations in ex vivo adipose tissues from nondiabetic and T2DM individuals. A multiscale mathematical model that integrates fluid dynamics and an intracellular insulin signaling pathway description was used for assisting microfluidic design in order to maximize measurement accuracy of tissue metabolic activity in response to perturbations. An automated microfluidic injection system was included on-chip for performing precise dynamic biochemical stimulations. The temporal evolution of culture conditions could be monitored for days, before and after perturbation, measuring glucose concentration in the outflow with high temporal resolution. As a proof of concept for detection of insulin resistance, we measured insulin-dependent glucose uptake by hAT from nondiabetic and T2DM subjects, mimicking the postprandial response. The system presented thus represents an important tool in dissecting the role of single tissues, such as hAT, in the complex interwoven picture of metabolic diseases.

The human brain is a tissue of vast complexity in terms of the cell types it comprises. Conventional approaches to classifying cell types in the human brain at single cell resolution have been limited to exploring relatively few markers and therefore have provided a limited molecular characterization of any given cell type. We used single cell RNA sequencing on 466 cells to capture the cellular complexity of the adult and fetal human brain at a whole transcriptome level. Healthy adult temporal lobe tissue was obtained during surgical procedures where otherwise normal tissue was removed to gain access to deeper hippocampal pathology in patients with medical refractory seizures. We were able to classify individual cells into all of the major neuronal, glial, and vascular cell types in the brain. We were able to divide neurons into individual communities and show that these communities preserve the categorization of interneuron subtypes that is typically observed with the use of classic interneuron markers. We then used single cell RNA sequencing on fetal human cortical neurons to identify genes that are differentially expressed between fetal and adult neurons and those genes that display an expression gradient that reflects the transition between replicating and quiescent fetal neuronal populations. Finally, we observed the expression of major histocompatibility complex type I genes in a subset of adult neurons, but not fetal neurons. The work presented here demonstrates the applicability of single cell RNA sequencing on the study of the adult human brain and constitutes a first step toward a comprehensive cellular atlas of the human brain.

Previous studies have proposed that epithelial to mesenchymal transition (EMT) in breast cancer cells regulates metastasis, stem cell properties and chemo-resistance; most studies were based on in vitro culture of cell lines and mouse transgenic cancer models. However, the identity and function of cells expressing EMT-associated genes in normal murine mammary gland homeostasis and human breast cancer still remains under debate. Using in vivo lineage tracing and triple negative breast cancer (TNBC) patient derived xenografts we demonstrate that the repopulating capacity in normal mammary epithelial cells and tumorigenic capacity in TNBC is independent of expression of EMT-associated genes. In breast cancer, while a subset of cells with epithelial and mesenchymal phenotypes have stem cell activity, in many cells that have lost epithelial characteristics with increased expression of mesenchymal genes, have decreased tumor-initiating capacity and plasticity. These findings have implications for the development of effective therapeutic agents targeting tumor-initiating cells., The contribution of EMT in mammary gland homeostasis and human breast cancer is still unclear. Here, using in vivo lineage tracing and breast cancer PDXs the authors demonstrate that the repopulating capacity in normal mammary epithelial cells and tumorigenic capacity in breast cancer is independent of expression of EMT-associated genes.

Significance The recent advances in cancer immunotherapy motivated us to investigate the clonal structure of the T cell receptor repertoire in breast tumors, normal breast, and blood in the same individuals. We found quantitatively distinct clonal structures in all three tissues, which enabled us to predict whether tissue is normal or tumor solely by comparing the repertoire of the tissue with blood. T cell receptor sequences shared between patients’ tumors are rare and, in general, do not appear to be specific to the cancer., Tumor-infiltrating T cells play an important role in many cancers, and can improve prognosis and yield therapeutic targets. We characterized T cells infiltrating both breast cancer tumors and the surrounding normal breast tissue to identify T cells specific to each, as well as their abundance in peripheral blood. Using immune profiling of the T cell beta-chain repertoire in 16 patients with early-stage breast cancer, we show that the clonal structure of the tumor is significantly different from adjacent breast tissue, with the tumor containing ∼2.5-fold greater density of T cells and higher clonality compared with normal breast. The clonal structure of T cells in blood and normal breast is more similar than between blood and tumor, and could be used to distinguish tumor from normal breast tissue in 14 of 16 patients. Many T cell sequences overlap between tissue and blood from the same patient, including ∼50% of T cells between tumor and normal breast. Both tumor and normal breast contain high-abundance “enriched” sequences that are absent or of low abundance in the other tissue. Many of these T cells are either not detected or detected with very low frequency in the blood, suggesting the existence of separate compartments of T cells in both tumor and normal breast. Enriched T cell sequences are typically unique to each patient, but a subset is shared between many different patients. We show that many of these are commonly generated sequences, and thus unlikely to play an important role in the tumor microenvironment.

We performed a high time-resolution, longitudinal study of normal pregnancy development by measuring cell-free RNA (cfRNA) in blood from women during each week of pregnancy. Analysis of tissue-specific transcripts in these samples enabled us to follow fetal and placental development with high resolution and sensitivity, and also to detect gene-specific responses of the maternal immune system to pregnancy. We established a “clock” for normal pregnancy development and enabled a direct molecular approach to determine expected delivery dates with comparable accuracy to ultrasound, creating the basis for a portable, inexpensive fetal dating method. We also identified a related gene set that accurately discriminated women at risk for spontaneous preterm delivery up to two months in advance of labor, forming the basis of a potential screening test for risk of preterm delivery.

SummaryHow transcription factors (TFs) reprogram one cell lineage to another remains unclear. Here, we define chromatin accessibility changes induced by the proneural TF Ascl1 throughout conversion of fibroblasts into induced neuronal (iN) cells. Thousands of genomic loci are affected as early as 12 hr after Ascl1 induction. Surprisingly, over 80% of the accessibility changes occur between days 2 and 5 of the 3-week reprogramming process. This chromatin switch coincides with robust activation of endogenous neuronal TFs and nucleosome phasing of neuronal promoters and enhancers. Subsequent morphological and functional maturation of iN cells is accomplished with relatively little chromatin reconfiguration. By integrating chromatin accessibility and transcriptome changes, we built a network model of dynamic TF regulation during iN cell reprogramming and identified Zfp238, Sox8, and Dlx3 as key TFs downstream of Ascl1. These results reveal a singular, coordinated epigenomic switch during direct reprogramming, in contrast to stepwise cell fate transitions in development.

Systemic sclerosis with pulmonary arterial hypertension (SSc-PAH) is a debilitating and frequently lethal disease of unknown cause lacking effective treatment options. Lymphocyte anomalies and autoantibodies observed in systemic sclerosis have suggested an autoimmune character. We study the clonal structure of the B cell repertoire in SSc-PAH using immunoglobulin heavy chain (IGH) sequencing before and after B cell depletion. We found SSc-PAH to be associated with anomalies in B cell development, namely, altered VDJ rearrangement frequencies (reduced IGHV2-5 segment usage) and an increased somatic mutation–fixation probability in expanded B cell lineages. SSc-PAH was also characterized by anomalies in B cell homeostasis, namely, an expanded immunoglobulin D–positive (IgD+) proportion with reduced mutation loads and an expanded proportion of highly antibody-secreting cells. Disease signatures pertaining to IGHV2-5 segment usage, IgD proportions, and mutation loads were temporarily reversed after B cell depletion. Analyzing the time course of B cell depletion, we find that the kinetics of naïve replenishment are predictable from baseline measurements alone, that release of plasma cells into the periphery can precede naïve replenishment, and that modes of B cell receptor diversity are highly elastic. Our findings reveal humoral immune signatures of SSc-PAH and uncover determinism in the effects of B cell depletion on the antibody repertoire.