Search

Neuroblastoma is a childhood cancer that resembles developmental stages of the neural crest. It is not established what developmental processes neuroblastoma cancer cells represent. Here, we sought to reveal the phenotype of neuroblastoma cancer cells by comparing cancer ( n = 19,723) with normal fetal adrenal single-cell transcriptomes ( n = 57,972). Our principal finding was that the neuroblastoma cancer cell resembled fetal sympathoblasts, but no other fetal adrenal cell type. The sympathoblastic state was a universal feature of neuroblastoma cells, transcending cell cluster diversity, individual patients, and clinical phenotypes. We substantiated our findings in 650 neuroblastoma bulk transcriptomes and by integrating canonical features of the neuroblastoma genome with transcriptional signals. Overall, our observations indicate that a pan-neuroblastoma cancer cell state exists, which may be attractive for novel immunotherapeutic and targeted avenues.

Periostin (POSTN) may serve as a biomarker for Type-2 mediated eosinophilic airway inflammation in asthma. We hypothesised that a Type-2 cytokine, interleukin (IL)-13, induces airway epithelial expression of POSTN, which in turn contributes to epithelial changes observed in asthma.We studied the effect of IL-13 on POSTN expression in BEAS-2B and air–liquid interface differentiated primary bronchial epithelial cells (PBECs). Additionally, the effects of recombinant human POSTN on epithelial-to-mesenchymal transition (EMT) markers and mucin genes were assessed. POSTN single cell gene expression and protein levels were analysed in bronchial biopsies and induced sputum from asthma patients and healthy controls.IL-13 increased POSTN expression in both cell types and this was accompanied by EMT-related features in BEAS-2B. In air–liquid interface differentiated PBECs, IL-13 increased POSTN basolateral and apical release. Apical administration of POSTN increased the expression of MMP-9, MUC5B and MUC5AC. In bronchial biopsies, POSTN expression was mainly confined to basal epithelial cells, ionocytes, endothelial cells and fibroblasts, showing higher expression in basal epithelial cells from asthma patients versus those from controls. A higher level of POSTN protein expression in epithelial and subepithelial layers was confirmed in bronchial biopsies from asthma patients when compared to healthy controls. Although sputum POSTN levels were not higher in asthma, levels correlated with eosinophil numbers and with the coughing-up of mucus.POSTN expression is increased by IL-13 in bronchial epithelial cells and is higher in bronchial biopsies from asthma patients. This may have important consequences, as administration of POSTN increases epithelial expression of mucin genes, supporting the relationship of POSTN with Type-2 mediated asthma and mucus secretion.

Cellular beauty is skin deep Human skin works as barrier, preventing the entry of pathogens, among other functions. Reynolds et al. used single-cell sequencing to generate an atlas of the human skin from both developing and adult sources, identifying differences and similarities across heterogeneous populations of skin cells. In this atlas, gene expression in the two disease states studied—atopic dermatitis and psoriasis—varied from that in a healthy adult, suggesting that a fetal skin signature is expressed in adult inflamed skin. Furthermore, differences in immune cell composition between healthy fetal and adult skin and that of individuals suffering from disease were observed. Science , this issue p. eaba6500

The COVID-19 pandemic, caused by SARS coronavirus 2 (SARS-CoV-2), has resulted in excess morbidity and mortality as well as economic decline. To characterise the systemic host immune response to SARS-CoV-2, we performed single-cell RNA-sequencing coupled with analysis of cell surface proteins, providing molecular profiling of over 800,000 peripheral blood mononuclear cells from a cohort of 130 patients with COVID-19. Our cohort, from three UK centres, spans the spectrum of clinical presentations and disease severities ranging from asymptomatic to critical. Three control groups were included: healthy volunteers, patients suffering from a non-COVID-19 severe respiratory illness and healthy individuals administered with intravenous lipopolysaccharide to model an acute inflammatory response. Full single cell transcriptomes coupled with quantification of 188 cell surface proteins, and T and B lymphocyte antigen receptor repertoires have provided several insights into COVID-19: 1. a new non-classical monocyte state that sequesters platelets and replenishes the alveolar macrophage pool; 2. platelet activation accompanied by early priming towards megakaryopoiesis in immature haematopoietic stem/progenitor cells and expansion of megakaryocyte-primed progenitors; 3. increased clonally expanded CD8+ effector:effector memory T cells, and proliferating CD4+ and CD8+ T cells in patients with more severe disease; and 4. relative increase of IgA plasmablasts in asymptomatic stages that switches to expansion of IgG plasmablasts and plasma cells, accompanied with higher incidence of BCR sharing, as disease severity increases. All data and analysis results are available for interrogation and data mining through an intuitive web portal. Together, these data detail the cellular processes present in peripheral blood during an acute immune response to COVID-19, and serve as a template for multi-omic single cell data integration across multiple centers to rapidly build powerful resources to help combat diseases such as COVID-19.

[Abstract]In late 2019 and through 2020, the COVID-19 pandemic swept the world, presenting both scientific and medical challenges associated with understanding and treating a previously unknown disease. To help address the need for great understanding of COVID-19, the scientific community mobilized and banded together rapidly to characterize SARS-CoV-2 infection, pathogenesis and its distinct disease trajectories. The urgency of COVID-19 provided a pressing use-case for leveraging relatively new tools, technologies, and nascent collaborative networks. Single-cell biology is one such example that has emerged over the last decade as a powerful approach that provides unprecedented resolution to the cellular and molecular underpinnings of biological processes. Early foundational work within the single-cell community, including the Human Cell Atlas, utilized published and unpublished data to characterize the putative target cells of SARS-CoV-2 sampled from diverse organs based on expression of the viral receptor ACE2 and associated entry factors TMPRSS2 and CTSL (Muus et al., 2020; Sungnak et al., 2020; Ziegler et al., 2020). This initial characterization of reference data provided an important foundation for framing infection and pathology in the airway as well as other organs. However, initial community analysis was limited to samples derived from uninfected donors and other previously-sampled disease indications. This report provides an overview of a single-cell data resource derived from samples from COVID-19 patients along with initial observations and guidance on data reuse and exploration.

Tumor cells may share some patterns of gene expression with their cell of origin, providing clues into the differentiation state and origin of cancer. Here, we study the differentiation state and cellular origin of 1300 childhood and adult kidney tumors. Using single cell mRNA reference maps of normal tissues, we quantify reference “cellular signals” in each tumor. Quantifying global differentiation, we find that childhood tumors exhibit fetal cellular signals, replacing the presumption of “fetalness” with a quantitative measure of immaturity. By contrast, in adult cancers our assessment refutes the suggestion of dedifferentiation towards a fetal state in most cases. We find an intimate connection between developmental mesenchymal populations and childhood renal tumors. We demonstrate the diagnostic potential of our approach with a case study of a cryptic renal tumor. Our findings provide a cellular definition of human renal tumors through an approach that is broadly applicable to human cancer., Transcriptomic analysis may provide information about the differentiation state and cell of origin of a cancer. Here, the authors assess mRNA signals in 1300 childhood and adult renal tumors and report a fetal origin of childhood tumors and no dedifferentiation of adult tumors.

The human skin confers biophysical and immunological protection through a complex cellular network that is established early in development. We profiled ~500,000 single cells using RNA-sequencing from healthy adult and developing skin, and skin from patients with atopic dermatitis and psoriasis. Our findings reveal a predominance of innate lymphoid cells and macrophages in developing skin in contrast to T cells and migratory dendritic cells in adult skin. We demonstrate dual keratinocyte differentiation trajectories and activated cellular circuits comprising vascular endothelial cells mediating immune cell trafficking, disease-specific clonally expanded IL13/IL22 and IL17A/F-expressing lymphocytes, epidermal IL23-expressing dendritic cells and inflammatory keratinocytes in disease. Our findings provide key insights into the dynamic cellular landscape of human skin in health and disease.One Sentence SummarySingle cell atlas of human skin reveals cell circuits which are quantitatively and qualitatively reconfigured in inflammatory skin disease.

BackgroundDespite its importance in health and disease, the cellular diversity of the lung endothelium has not been systematically characterized in humans. Here we provide a reference atlas of human lung endothelial cells (ECs), to facilitate a better understanding of the phenotypic diversity and composition of cells comprising the lung endothelium, both in health and disease.MethodsWe reprocessed control single cell RNA sequencing (scRNAseq) data from five datasets of whole lungs that were used for the analysis of pan-endothelial markers, we later included a sixth dataset of sorted control EC for the vascular subpopulation analysis. EC populations were characterized through iterative clustering with subsequent differential expression analysis. Marker genes were validated by immunohistochemistry andin situhybridization. Signaling network between different lung cell types was studied using connectomic analysis. For cross species analysis we applied the same methods to scRNAseq data obtained from mouse lungs.ResultsThe six lung scRNAseq datasets were reanalyzed and annotated to identify over 15,000 vascular EC cells from 73 individuals. Differential expression analysis of EC revealed signatures corresponding to endothelial lineage, including pan-endothelial, pan-vascular and subpopulation-specific marker gene sets. Beyond the broad cellular categories of lymphatic, capillary, arterial and venous ECs we found previously indistinguishable subpopulations; among venous EC we identified two previously indistinguishable populations, pulmonary-venous ECs (COL15A1neg) localized to the lung parenchyma and systemic-venous ECs (COL15A1pos) localized to the airways and the visceral pleura; among capillary EC we confirmed their subclassification into recently discovered aerocytes characterized by EDNRB, SOSTDC1 and TBX2 and general capillary EC. We confirmed that all six endothelial cell types, including the systemic-venous EC and aerocytes are present in mice and identified endothelial marker genes conserved in humans and mice. Ligand-Receptor connectome analysis revealed important homeostatic crosstalk of EC with other lung resident cell types. Our manuscript is accompanied by an online data mining tool (www.LungEndothelialCellAtlas.com).ConclusionOur integrated analysis provides the comprehensive and well-crafted reference atlas of lung endothelial cells in the normal lung and confirms and describes in detail previously unrecognized endothelial populations across a large number of humans and mice.

The Human Developmental Cell Atlas (HDCA), as part of the Human Cell Atlas, aims to generate a comprehensive reference map of cells during development. This detailed study of development will be critical for understanding normal organogenesis, the impact of mutations, environmental factors and infectious agents on congenital and childhood disorders, and the molecular cellular basis of ageing, cancer and regenerative medicine. In this perspective, we outline the challenges of mapping and modelling human development using state of the art technologies to create a reference atlas across gestation for scientific and clinical benefit. We discuss the potential value of HDCA to enhance human pluripotent stem cell-derived organoid model systems and, in turn, the use of organoids and animal models to inform HDCA. Finally, we provide a roadmap towards a complete atlas of human development.

Neuroblastoma is an embryonal childhood cancer that arises from aberrant development of the neural crest, mostly within the fetal adrenal medulla. It is not established what developmental processes neuroblastoma cancer cells represent. Here, we sought to reveal the phenotype of neuroblastoma cancer cells by comparing cancer (n=16,591) with fetal adrenal single cell transcriptomes (n=57,972). Our principal finding was that the neuroblastoma cancer cell resembled fetal sympathoblasts, but no other fetal adrenal cell type. The sympathoblastic state was a universal feature of neuroblastoma cells, transcending cell cluster diversity, individual patients and clinical phenotypes. We substantiated our findings in 652 neuroblastoma bulk transcriptomes and by integrating canonical features of the neuroblastoma genome with transcriptional signals. Overall, our observations indicate that there exists a pan-neuroblastoma cancer cell state which may be an attractive target for novel therapeutic avenues.

Identifying the origin of cancer Many cancers are classified on the basis of the organ or tissue from which they originated. However, identifying the specific cells and conditions that precede tumorigenesis can help us understand and better treat the resulting disease. Nowicki-Osuch et al . used a single-cell approach to investigate the cell of origin for Barrett’s esophagus (BE) and the mechanisms leading to the development of esophageal adenocarcinoma (EAC) (see the Perspective by Geboes and Hoorens). Analyses of healthy human esophageal tissues, mutational lineage tracing, and organoid models revealed that BE originates from the gastric cardia and that EAC arises from undifferentiated BE cells. This analysis provides a map of the transcriptional landscape of the healthy esophagus that can be compared with mouse models of disease. —LMZ

Single-cell RNA-Seq (scRNA-Seq) data analysis requires expertise in command-line tools, programming languages and scaling on compute infrastructure. As scRNA-Seq becomes widespread, computational pipelines need to be more accessible, simpler and scalable. We introduce an interactive analysis environment for scRNA-Seq, based on Galaxy, with ~70 functions from major single-cell analysis tools, which can be run on compute clusters, cloud providers or single machines, to bring compute to the data in scRNA-Seq.

Thymus development, cell by cell The human thymus is the organ responsible for the maturation of many types of T cells, which are immune cells that protect us from infection. However, it is not well known how these cells develop with a full immune complement that contains the necessary variation to protect us from a variety of pathogens. By performing single-cell RNA sequencing on more than 250,000 cells, Park et al. examined the changes that occur in the thymus over the course of a human life. They found that development occurs in a coordinated manner among immune cells and with their developmental microenvironment. These data allowed for the creation of models of how T cells with different specific immune functions develop in humans. Science , this issue p. eaay3224

Gastrointestinal microbiota and immune cells interact closely and display regional specificity, but little is known about how these communities differ with location. Here, we simultaneously assess microbiota and single immune cells across the healthy, adult human colon, with paired characterisation of immune cells in the mesenteric lymph nodes, to delineate colonic immune niches at steady-state. We describe distinct T helper cell activation and migration profiles along the colon and characterise the transcriptional adaptation trajectory of T regulatory cells between lymphoid tissue and colon. Finally, we show increasing B cell accumulation, clonal expansion and mutational frequency from caecum to sigmoid colon, and link this to the increasing number of reactive bacterial species.

With the recent transformative developments in single-cell genomics and, in particular, single-cell gene expression analysis, it is now possible to study tissues at the single-cell level, rather than having to rely on data from bulk measurements. Here we review the rapid developments in single-cell RNA sequencing (scRNA-seq) protocols that have the potential for unbiased identification and profiling of all cell types within a tissue or organism. In addition, novel approaches for spatial profiling of gene expression allow us to map individual cells and cell types back into the three-dimensional context of organs. The combination of in-depth single-cell and spatial gene expression data will reveal tissue architecture in unprecedented detail, generating a wealth of biological knowledge and a better understanding of many diseases.

The thymus provides a nurturing environment for the differentiation and selection of T cells, a process orchestrated by their interaction with multiple thymic cell types. We utilised single-cell RNA-sequencing (scRNA-seq) to create a cell census of the human thymus and to reconstruct T-cell differentiation trajectories and T-cell receptor (TCR) recombination kinetics. Using this approach, we identified and located in situ novel CD8αα+ T-cell populations, thymic fibroblast subtypes and activated dendritic cell (aDC) states. In addition, we reveal a bias in TCR recombination and selection, which is attributed to genomic position and suggests later commitment of the CD8+ T-cell lineage. Taken together, our data provide a comprehensive atlas of the human thymus across the lifespan with new insights into human T-cell development.

International audience; Definitive haematopoiesis in the fetal liver supports self-renewal anddifferentiation of haematopoietic stem cells and multipotent progenitors(HSC/MPPs) but remains poorly defined in humans. Here, using single-celltranscriptome profiling of approximately 140,000 liver and 74,000 skin, kidneyand yolk sac cells, we identify the repertoire of human blood and immune cellsduring development. We infer differentiation trajectories from HSC/MPPs andevaluate the influence of the tissue microenvironment on blood and immune celldevelopment. We reveal physiological erythropoiesis in fetal skin and thepresence of mast cells, natural killer and innate lymphoid cell precursors in theyolk sac. We demonstrate a shift in the haemopoietic composition of fetal liverduring gestation away from being predominantly erythroid, accompanied by aparallel change in differentiation potential of HSC/MPPs, which we functionallyvalidate. Our integrated map of fetal liver haematopoiesis provides a blueprintfor the study of paediatric blood and immune disorders, and a reference forharnessing the therapeutic potential of HSC/MPPs.

Background Improving lung cancer risk assessment is required because current early-detection screening criteria miss most cases. We therefore examined the utility for lung cancer risk assessment of a DNA Repair score obtained from OGG1, MPG, and APE1 blood tests. In addition, we examined the relationship between the level of DNA repair and global gene expression. Methods We conducted a blinded case-control study with 150 non–small cell lung cancer case patients and 143 control individuals. DNA Repair activity was measured in peripheral blood mononuclear cells, and the transcriptome of nasal and bronchial cells was determined by RNA sequencing. A combined DNA Repair score was formed using logistic regression, and its correlation with disease was assessed using cross-validation; correlation of expression to DNA Repair was analyzed using Gene Ontology enrichment. Results DNA Repair score was lower in case patients than in control individuals, regardless of the case’s disease stage. Individuals at the lowest tertile of DNA Repair score had an increased risk of lung cancer compared to individuals at the highest tertile, with an odds ratio (OR) of 7.2 (95% confidence interval [CI] = 3.0 to 17.5; P Conclusions The DNA Repair score, previously shown to be a lung cancer risk factor in the Israeli population, was validated in this independent study as a mechanism-based cancer risk biomarker and can substantially improve current lung cancer risk prediction, assisting prevention and early detection by computed tomography scanning.

BackgroundThe Human Cell Atlas is a large international collaborative effort to map all cell types of the human body. Single cell RNA sequencing can generate high quality data for the delivery of such an atlas. However, delays between fresh sample collection and processing may lead to poor data and difficulties in experimental design. Despite this, there has not yet been a systematic assessment of the effect of cold storage time on the quality of scRNAseqResultsThis study assessed the effect of cold storage on fresh healthy spleen, oesophagus and lung from ≥5 donors over 72 hours. We collected 240,000 high quality single cell transcriptomes with detailed cell type annotations and whole genome sequences of donors, enabling future eQTL studies. Our data provide a valuable resource for the study of these three organs and will allow cross-organ comparison of cell types.We see little effect of cold ischaemic time on cell viability, yield, total number of reads per cell and other quality control metrics in any of the tissues within the first 24 hours. However, we observed higher percentage of mitochondrial reads, indicative of cellular stress, and increased contamination by background “ambient RNA” reads in the 72h samples in spleen, which is cell type specific.ConclusionsIn conclusion, we present robust protocols for tissue preservation for up to 24 hours prior to scRNAseq analysis. This greatly facilitates the logistics of sample collection for Human Cell Atlas or clinical studies since it increases the time frames for sample processing.

Adult stem-cells may serve as the cell-of-origin for cancer, yet their unbiased identification in single cell RNA sequencing data is challenging due to the high dropout rate. In the case of breast, the existence of a bipotent stem-like state is also controversial. Here we apply a marker-free algorithm to scRNA-Seq data from the human mammary epithelium, revealing a high-potency cell-state enriched for an independent mammary stem-cell expression module. We validate this stem-like state in independent scRNA-Seq data. Our algorithm further predicts that the stem-like state is bipotent, a prediction we are able to validate using FACS sorted bulk expression data. The bipotent stem-like state correlates with clinical outcome in basal breast cancer and is characterized by overexpression of YBX1 and ENO1, two modulators of basal breast cancer risk. This study illustrates the power of a marker-free computational framework to identify a novel bipotent stem-like state in the mammary epithelium., Chen et al build on a marker-free entropy-based algorithm to assign specific potency states to individual cells using single-cell RNA sequencing data. They identify a putative bipotent stem-cell like state, which associates with basal breast cancer risk and poor clinical outcome.

Human lungs enable efficient gas exchange and form an interface with the environment, which depends on mucosal immunity for protection against infectious agents. Tightly controlled interactions between structural and immune cells are required to maintain lung homeostasis. Here, we use single-cell transcriptomics to chart the cellular landscape of upper and lower airways and lung parenchyma in healthy lungs, and lower airways in asthmatic lungs. We report location-dependent airway epithelial cell states and a novel subset of tissue-resident memory T cells. In the lower airways of patients with asthma, mucous cell hyperplasia is shown to stem from a novel mucous ciliated cell state, as well as goblet cell hyperplasia. We report the presence of pathogenic effector type 2 helper T cells (TH2) in asthmatic lungs and find evidence for type 2 cytokines in maintaining the altered epithelial cell states. Unbiased analysis of cell-cell interactions identifies a shift from airway structural cell communication in healthy lungs to a TH2-dominated interactome in asthmatic lungs.

SummaryDefinitive haematopoiesis in the fetal liver supports self-renewal and differentiation of haematopoietic stem cells/multipotent progenitors (HSC/MPPs), yet remains poorly defined in humans. Using single cell transcriptome profiling of ~133,000 fetal liver and ~65,000 fetal skin and kidney cells, we identify the repertoire of blood and immune cells in first and early second trimesters of development. From this data, we infer differentiation trajectories from HSC/MPPs, and evaluate the impact of tissue microenvironment on blood and immune cell development. We predict coupling of mast cell differentiation with erythro-megakaryopoiesis and identify physiological erythropoiesis in fetal skin. We demonstrate a shift in fetal liver haematopoietic composition during gestation away from being erythroid-predominant, accompanied by a parallel change in HSC/MPP differentiation potential, which we functionally validate. Our integrated map of fetal liver haematopoiesis provides a blueprint for the study of paediatric blood and immune disorders, and a valuable reference for understanding and harnessing the therapeutic potential of HSC/MPPs.

Temporal resolution of cellular features associated with a severe COVID-19 disease trajectory is required for understanding skewed immune responses and finding outcome predictors. Here, we performed a longitudinal multi-omics study using a two-centre cohort of 14 patients. We analysed the bulk transcriptome, bulk DNA methylome, and single-cell transcriptome (>358,000 cells, including BCR profiles) of peripheral blood samples harvested from up to 5 time points. Validation was performed in two independent cohorts of COVID-19 patients. Severe COVID-19 was characterized by an increase of proliferating, metabolically hyperactive plasmablasts. Coinciding with critical illness, we also identified an expansion of IFN-activated circulating megakaryocytes and increased erythropoiesis with features of hypoxic signalling. Megakaryocyte- and erythroid cell-derived co-expression modules were predictive of fatal disease outcome. The study demonstrates broad cellular effects of SARS-CoV-2 infection beyond classical immune cells and may serve as an entry point to develop biomarkers and targeted treatments of patients with COVID-19., Graphical Abstract, Highlights • SARS-CoV2 infection elicits dynamic changes of circulating cells in the blood • Severe COVID-19 is characterized by increased metabolically active plasmablasts • Elevation of IFN-activated megakaryocytes and erythroid cells in severe COVID-19 • Cell-type specific expression signatures are associated with a fatal COVID-19 outcome, Bernardes et al. explore COVID-19 disease trajectories by performing longitudinal multi-omics analyses in peripheral blood samples from hospitalized patients. The analyses identify increased numbers of plasmablasts, interferon-activated megakaryocytes and erythroid cells as hallmarks of severe disease, and define molecular signatures linked to a fatal COVID-19 disease outcome.

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 identification of progenitor and stem like cells in epithelial tissues, as well as those that may serve as the cell of origin for epithelial cancers, is an outstanding challenge. Here we present a novel algorithm, called LandSCENT, which constructs a 3-dimensional integrated landscape of cell-states, encompassing cell-potency and expression subtypes, to facilitate the identification of progenitor and stem-like cells. Application to thousands of single-cell RNA-Seq profiles from the normal mammary epithelium reveals a rare 5% subpopulation of highly potent single-cells. The integrated landscape naturally predicts that these cells define a bi-potent-like state, a result not obtainable via standard methods or without invoking prior assumptions. The bi-potent-like cells are overrepresented within the basal compartment but also overlap with an immature luminal phenotype. We characterize the transcriptome of these cells and show that is enriched for a mammary stem-cell module. We further identifyYBX1, a regulator of breast cancer risk identified from GWAS, as the key transcription factor defining this candidate bi-potent cellular phenotype. We validate the putative bi-potency ofYBX1-marked cells using independent FACS-sorted bulk expression data. In addition,YBX1is overexpressed in basal breast cancer and correlates with clinical outcome. In summary, we here provide a novel computational framework which may serve to identify and prioritize candidate normal or cancer progenitor/stem-like single-cell phenotypes, for subsequent functional studies.

The recent development of single cell gene expression technologies, and especially single cell transcriptomics, have revolutionized the way biologists and clinicians investigate organs and organisms, allowing an unprecedented level of resolution to the description of cell demographics in both healthy and diseased states. Single cell transcriptomics provide information on prevalence, heterogeneity, and gene co-expression at the individual cell level. This enables a cell-centric outlook to define intracellular gene regulatory networks and to bridge toward the definition of intercellular pathways otherwise masked in bulk analysis. The technologies have developed at a fast pace producing a multitude of different approaches, with several alternatives to choose from at any step, including single cell isolation and capturing, lysis, RNA reverse transcription and cDNA amplification, library preparation, sequencing, and computational analyses. Here, we provide guidelines for the experimental design of single cell RNA sequencing experiments, exploring the current options for the crucial steps. Furthermore, we provide a complete overview of the typical data analysis workflow, from handling the raw sequencing data to making biological inferences. Significantly, advancements in single cell transcriptomics have already contributed to outstanding exploratory and functional studies of cardiac development and disease models, as summarized in this review. In conclusion, we discuss achievable outcomes of single cell transcriptomics' applications in addressing unanswered questions and influencing future cardiac clinical applications.

During early human pregnancy the uterine mucosa transforms into the decidua, into which the fetal placenta implants and where placental trophoblast cells intermingle and communicate with maternal cells. Trophoblast-decidual interactions underlie common diseases of pregnancy, including pre-eclampsia and stillbirth. Here we profile the transcriptomes of about 70,000 single cells from first-trimester placentas with matched maternal blood and decidual cells. The cellular composition of human decidua reveals subsets of perivascular and stromal cells that are located in distinct decidual layers. There are three major subsets of decidual natural killer cells that have distinctive immunomodulatory and chemokine profiles. We develop a repository of ligand-receptor complexes and a statistical tool to predict the cell-type specificity of cell-cell communication via these molecular interactions. Our data identify many regulatory interactions that prevent harmful innate or adaptive immune responses in this environment. Our single-cell atlas of the maternal-fetal interface reveals the cellular organization of the decidua and placenta, and the interactions that are critical for placentation and reproductive success.

Background The IRE1a-XBP1 pathway is a conserved adaptive mediator of the unfolded protein response. The pathway is indispensable for the development of secretory cells by facilitating protein folding and enhancing secretory capacity. In the immune system, it is known to function in dendritic cells, plasma cells, and eosinophil development and differentiation, while its role in T helper cell is unexplored. Here, we investigated the role of the IRE1a-XBP1 pathway in regulating activation and differentiation of type-2 T helper cell (Th2), a major T helper cell type involved in allergy, asthma, helminth infection, pregnancy, and tumor immunosuppression. Methods We perturbed the IRE1a-XBP1 pathway and interrogated its role in Th2 cell differentiation. We performed genome-wide transcriptomic analysis of differential gene expression to reveal IRE1a-XBP1 pathway-regulated genes and predict their biological role. To identify direct target genes of XBP1 and define XBP1’s regulatory network, we performed XBP1 ChIPmentation (ChIP-seq). We validated our predictions by flow cytometry, ELISA, and qPCR. We also used a fluorescent ubiquitin cell cycle indicator mouse to demonstrate the role of XBP1 in the cell cycle. Results We show that Th2 lymphocytes induce the IRE1a-XBP1 pathway during in vitro and in vivo activation. Genome-wide transcriptomic analysis of differential gene expression by perturbing the IRE1a-XBP1 pathway reveals XBP1-controlled genes and biological pathways. Performing XBP1 ChIPmentation (ChIP-seq) and integrating with transcriptomic data, we identify XBP1-controlled direct target genes and its transcriptional regulatory network. We observed that the IRE1a-XBP1 pathway controls cytokine secretion and the expression of two Th2 signature cytokines, IL13 and IL5. We also discovered that the IRE1a-XBP1 pathway facilitates activation-dependent Th2 cell proliferation by facilitating cell cycle progression through S and G2/M phase. Conclusions We confirm and detail the critical role of the IRE1a-XBP1 pathway during Th2 lymphocyte activation in regulating cytokine expression, secretion, and cell proliferation. Our high-quality genome-wide XBP1 ChIP and gene expression data provide a rich resource for investigating XBP1-regulated genes. We provide a browsable online database available at http://data.teichlab.org. Electronic supplementary material The online version of this article (10.1186/s13073-018-0589-3) contains supplementary material, which is available to authorized users.

SummaryDuring the early weeks of human pregnancy, the fetal placenta implants into the uterine mucosa (decidua) where placental trophoblast cells intermingle and communicate with maternal cells. Here, we profile transcriptomes of ∼50,000 single cells from this unique microenvironment, sampling matched first trimester maternal blood and decidua, and fetal cells from the placenta itself. We define the cellular composition of human decidua, revealing five distinct subsets of decidual fibroblasts with differing growth factors and hormone production profiles, and show that fibroblast states define two distinct decidual layers. Among decidual NK cells, we resolve three subsets, each with a different immunomodulatory and chemokine profile. We develop a repository of ligand-receptor pairs (www.CellPhoneDB.org) and a statistical tool to predict the probability of cell-cell interactions via these pairs, highlighting specific interactions between decidual NK cells and invading fetal extravillous trophoblast cells, maternal immune and stromal cells. Our single cell atlas of the maternal-fetal interface reveals the cellular organization and interactions critical for placentation and reproductive success.

Motivation Transcription factors (TFs) are key regulators of gene expression, and can activate or repress multiple target genes, forming regulatory units, or regulons. Understanding downstream effects of these regulators includes evaluating how TFs cooperate or compete within regulatory networks. Here we present RTNduals, an R/Bioconductor package that implements a general method for analyzing pairs of regulons. Results RTNduals identifies a dual regulon when the number of targets shared between a pair of regulators is statistically significant. The package extends the RTN (Reconstruction of Transcriptional Networks) package, and uses RTN transcriptional networks to identify significant co-regulatory associations between regulons. The Supplementary Information reports two case studies for TFs using the METABRIC and TCGA breast cancer cohorts. Availability and implementation RTNduals is written in the R language, and is available from the Bioconductor project at http://bioconductor.org/packages/RTNduals/. Supplementary information Supplementary data are available at Bioinformatics online.

Increasing numbers of large scale single cell RNA-Seq projects are leading to a data explosion, which can only be fully exploited through data integration. Therefore, efficient computational tools for combining diverse datasets are crucial for biology in the single cell genomics era. A number of methods have been developed to assist data integration by removing technical batch effects, but most are computationally intensive. To overcome the challenge of enormous datasets, we have developed BBKNN, an extremely fast graph-based data integration method. We illustrate the power of BBKNN for dimensionalityreduced visualisation and clustering in multiple biological scenarios, including a massive integrative study over several murine atlases. BBKNN successfully connects cell populations across experimentally heterogeneous mouse scRNA-Seq datasets, which reveals global markers of cell type and organspecificity and provides the foundation for inferring the underlying transcription factor network. BBKNN is available athttps://github.com/Teichlab/bbknn.

Two opposing clusters of transcription factors (TF) have been associated with the differential risks of estrogen receptor positive or negative breast cancers, but the mechanisms underlying the opposing functions of the two clusters are undefined. In this study, we identified NFIB and YBX1 as novel interactors of the estrogen receptor (ESR1). NFIB and YBX1 are both risk TF associated with progression of ESR1-negative disease. Notably, they both interacted with the ESR1-FOXA1 complex and inhibited the transactivational potential of ESR1. Moreover, signaling through FGFR2, a known risk factor in breast cancer development, augmented these interactions and further repressed ESR1 target gene expression. We therefore show that members of two opposing clusters of risk TFs associated with ESR1-positive and -negative breast cancer can physically interact. We postulate that this interaction forms a toggle between two developmental pathways affected by FGFR2 signaling, possibly offering a junction to exploit therapeutically. Significance: Binding of the transcription factors NFIB and YBX1 to the estrogen receptor can promote an estrogen-independent phenotype that can be reverted by inhibiting FGFR2 signaling. Cancer Res; 78(2); 410–21. ©2017 AACR.

We have recently identified transcription factors (TFs) that are key drivers of breast cancer risk. To better understand the pathways or sub-networks in which these TFs mediate their function we sought to identify upstream modulators of their activity. We applied the MINDy (Modulator Inference by Network Dynamics) algorithm to four TFs (ESR1, FOXA1, GATA3 and SPDEF) that are key drivers of estrogen receptor-positive (ER+) breast cancer risk, as well as cancer progression. Our computational analysis identified over 500 potential modulators. We assayed 189 of these and identified 55 genes with functional characteristics that were consistent with a role as TF modulators. In the future, the identified modulators may be tested as potential therapeutic targets, able to alter the activity of TFs that are critical in the development of breast cancer.

Allele-specific measurements of transcription factor binding from ChIP-seq data are key to dissecting the allelic effects of non-coding variants and their contribution to phenotypic diversity. However, most methods of detecting an allelic imbalance assume diploid genomes. This assumption severely limits their applicability to cancer samples with frequent DNA copy-number changes. Here we present a Bayesian statistical approach called BaalChIP to correct for the effect of background allele frequency on the observed ChIP-seq read counts. BaalChIP allows the joint analysis of multiple ChIP-seq samples across a single variant and outperforms competing approaches in simulations. Using 548 ENCODE ChIP-seq and six targeted FAIRE-seq samples, we show that BaalChIP effectively corrects allele-specific analysis for copy-number variation and increases the power to detect putative cis-acting regulatory variants in cancer genomes. Electronic supplementary material The online version of this article (doi:10.1186/s13059-017-1165-7) contains supplementary material, which is available to authorized users.

T helper type 2 (Th2) cells are important regulators of mammalian adaptive immunity and have relevance for infection, auto-immunity and tumour immunology. Using a newly developed, genome-wide retroviral CRISPR knock-out (KO) library, combined with RNA-seq, ATAC-seq and ChIP-seq, we have dissected the regulatory circuitry governing activation (including proliferation) and differentiation of these cells. Our experiments distinguish cell activationversusdifferentiation in a quantitative framework. We demonstrate that these two processes are tightly coupled and are jointly controlled by many transcription factors, metabolic genes and cytokine/receptor pairs. There is only a small number of genes regulating differentiation without any role in activation. Our study provides an atlas for the T helper cell regulatory network, pinpointing key players of Th2 differentiation and demonstrating remarkable plasticity between the diverse T helper cell fates. We provide an online resource for interactive data querying at:http://data.teichlab.org.

The 10q26 locus in the second intron of FGFR2 is the locus most strongly associated with estrogen-receptor-positive breast cancer in genome-wide association studies. We conducted fine-scale mapping in case-control studies genotyped with a custom chip (iCOGS), comprising 41 studies (n = 89,050) of European ancestry, 9 Asian ancestry studies (n = 13,983), and 2 African ancestry studies (n = 2,028) from the Breast Cancer Association Consortium. We identified three statistically independent risk signals within the locus. Within risk signals 1 and 3, genetic analysis identified five and two variants, respectively, highly correlated with the most strongly associated SNPs. By using a combination of genetic fine mapping, data on DNase hypersensitivity, and electrophoretic mobility shift assays to study protein-DNA binding, we identified rs35054928, rs2981578, and rs45631563 as putative functional SNPs. Chromatin immunoprecipitation showed that FOXA1 preferentially bound to the risk-associated allele (C) of rs2981578 and was able to recruit ERα to this site in an allele-specific manner, whereas E2F1 preferentially bound the risk variant of rs35054928. The risk alleles were preferentially found in open chromatin and bound by Ser5 phosphorylated RNA polymerase II, suggesting that the risk alleles are associated with changes in transcription. Chromatin conformation capture demonstrated that the risk region was able to interact with the promoter of FGFR2, the likely target gene of this risk region. A role for FOXA1 in mediating breast cancer susceptibility at this locus is consistent with the finding that the FGFR2 risk locus primarily predisposes to estrogen-receptor-positive disease. © 2013 The American Society of Human Genetics.

Summary The fibroblast growth factor receptor 2 (FGFR2) locus is the ‘top hit’ in genome-wide association studies for breast cancer. Here, we examine the effect of FGFR2 signalling on transcriptional networks in breast cancer and propose a mechanism for FGFR2 risk single-nucleotide polymorphism function., The fibroblast growth factor receptor 2 (FGFR2) locus is consistently the top hit in genome-wide association studies for oestrogen receptor-positive (ER+) breast cancer. Yet, its mode of action continues to be controversial. Here, we employ a systems biology approach to demonstrate that signalling via FGFR2 counteracts cell activation by oestrogen. In the presence of oestrogen, the oestrogen receptor (ESR1) regulon (set of ESR1 target genes) is in an active state. However, signalling by FGFR2 is able to reverse the activity of the ESR1 regulon. This effect is seen in multiple distinct FGFR2 signalling model systems, across multiple cells lines and is dependent on the presence of FGFR2. Increased oestrogen exposure has long been associated with an increased risk of breast cancer. We therefore hypothesized that risk variants should reduce FGFR2 expression and subsequent signalling. Indeed, transient transfection experiments assaying the three independent variants of the FGFR2 risk locus (rs2981578, rs35054928 and rs45631563) in their normal chromosomal context show that these single-nucleotide polymorphisms (SNPs) map to transcriptional silencer elements and that, compared with wild type, the risk alleles augment silencer activity. The presence of risk variants results in lower FGFR2 expression and increased oestrogen responsiveness. We thus propose a molecular mechanism by which FGFR2 can confer increased breast cancer risk that is consistent with oestrogen exposure as a major driver of breast cancer risk. Our findings may have implications for the clinical use of FGFR2 inhibitors.

Breast cancer exhibits familial aggregation, consistent with variation in genetic susceptibility to the disease. Known susceptibility genes account for less than 25% of the familial risk of breast cancer, and the residual genetic variance is likely to be due to variants conferring more moderate risks. To identify further susceptibility alleles, we conducted a two-stage genome-wide association study in 4,398 breast cancer cases and 4,316 controls, followed by a third stage in which 30 single nucleotide polymorphisms (SNPs) were tested for confirmation in 21,860 cases and 22,578 controls from 22 studies. We used 227,876 SNPs that were estimated to correlate with 77% of known common SNPs in Europeans at r2 > 0.5. SNPs in five novel independent loci exhibited strong and consistent evidence of association with breast cancer (P

The transcription factor E2A plays a crucial role in B cell development, the control of immunoglobulin gene rearrangement and expression. Here we report that in primary mouse B cells lipopolysaccharide (LPS) is able to induce the level of E2A protein by over 50-fold in 3 days of culture. In contrast, CD40 signaling is insufficient to cause an E2A increase and can in fact prevent the LPS-mediated induction of E2A. These results suggest that E2A induction requires both proliferation and differentiation. We find that E2A protein induction is regulated post-transcriptionally since E2A mRNA is not induced by LPS. We have thus identified an important additional layer of regulation affecting the activity of E2A transcription factors.

Upon treatment with lipopolysaccharide (LPS), primary B cells proliferate and differentiate into plasma cells with concomitant up-regulation of immunoglobulin (Ig) gene expression. Here we examine the role of the Igkappa 3' enhancer in this process using a kappa3'-enhancer-driven beta-globin reporter gene in transgenic mice. We find that LPS treatment up-regulates kappa3' enhancer activity as a function of differentiation rather than proliferation, since proliferation only (induced by cross-linking of CD40) is insufficient to activate the element, whilst differentiation with only limited proliferation (LPS + transforming growth factor-beta) does allow activation to occur. The Igkappa 3' enhancer is also induced by cross-linking of surface Ig and this signal can synergize with LPS activation, suggesting that distinct activation pathways are used. Nevertheless, both of these pathways can be inhibited by co-cross-linking of CD40. Thus Ig enhancers in the heavy and light chain loci are differentially regulated in response to CD40.

The Ig kappa 3' enhancer is required for high levels of Ig kappa gene expression. We now show that kappa 3' enhancer function increases five- to eightfold after stimulation of primary murine B cells with phorbol 12-myristate 13-acetate (PMA) and the calcium ionophore ionomycin. In the presence of cyclosporin A this induction is almost halved, suggesting that transcription factors of the NFAT family contribute to kappa 3' enhancer induction. Indeed, we identify a novel NFAT binding site which is required for full enhancer function. We find that this site is transcriptionally active in stimulated B cells, T cells and fibroblasts and that both PMA and ionomycin are required for maximal induction. Time course analysis of the components of the protein-DNA complex in primary lymphocytes reveals that both NFATp and NFATc are present in the complex after 15 min, while only NFATc is detectable after 4 h. This suggests that NFATc plays the dominant role in controlling long-term responses of this transcription factor family. Furthermore, JunB, JunD, FosB and cFos form part of the DNA-protein complex in Bal-17 B cells. Complex formation as well as transcriptional activity can also be induced by crosslinking of surface Ig. We have, thus, identified a unique NFAT complex in B cells that contributes to Ig kappa gene expression.

Background:Breast cancer is one of the most common malignancies in women. Genome-wide association studies have identified FGFR2 as a breast cancer susceptibility gene. Common variation in other fibroblast growth factor (FGF) receptors might also modify risk. We tested this hypothesis by studying genotyped single-nucleotide polymorphisms (SNPs) and imputed SNPs in FGFR1, FGFR3, FGFR4 and FGFRL1 in the Breast Cancer Association Consortium.\ud \ud Methods:Data were combined from 49 studies, including 53 835 cases and 50 156 controls, of which 89 050 (46 450 cases and 42 600 controls) were of European ancestry, 12 893 (6269 cases and 6624 controls) of Asian and 2048 (1116 cases and 932 controls) of African ancestry. Associations with risk of breast cancer, overall and by disease sub-type, were assessed using unconditional logistic regression.\ud \ud Results:Little evidence of association with breast cancer risk was observed for SNPs in the FGF receptor genes. The strongest evidence in European women was for rs743682 in FGFR3; the estimated per-allele odds ratio was 1.05 (95 confidence interval=1.02-1.09, P=0.0020), which is substantially lower than that observed for SNPs in FGFR2.\ud \ud Conclusion:Our results suggest that common variants in the other FGF receptors are not associated with risk of breast cancer to the degree observed for FGFR2. © 2014 Cancer Research UK.

The fibroblast growth factor receptor 2 (FGFR2) locus has been consistently identified as a breast cancer risk locus in independent genome-wide association studies. However, the molecular mechanisms underlying FGFR2-mediated risk are still unknown. Using model systems we show that FGFR2-regulated genes are preferentially linked to breast cancer risk loci in expression quantitative trait loci analysis, supporting the concept that risk genes cluster in pathways. Using a network derived from 2,000 transcriptional profiles we identify SPDEF, ERα, FOXA1, GATA3 and PTTG1 as master regulators of fibroblast growth factor receptor 2 signalling, and show that ERα occupancy responds to fibroblast growth factor receptor 2 signalling. Our results indicate that ERα, FOXA1 and GATA3 contribute to the regulation of breast cancer susceptibility genes, which is consistent with the effects of anti-oestrogen treatment in breast cancer prevention, and suggest that fibroblast growth factor receptor 2 signalling has an important role in mediating breast cancer risk., FGFR2 gene variation is associated with breast cancer risk but the molecular mechanism is unknown. Fletcher et al. provide a link between FGFR2 signalling and breast cancer susceptibility by demonstrating that FGFR2 signalling activates the ERa transcriptional network, which drives transcription of risk genes.

Analysis of 4,405 variants in 89,050 European subjects from 41 case-control studies identified three independent association signals for estrogen-receptor-positive tumors at 11q13. The strongest signal maps to a transcriptional enhancer element in which the G allele of the best candidate causative variant rs554219 increases risk of breast cancer, reduces both binding of ELK4 transcription factor and luciferase activity in reporter assays, and may be associated with low cyclin D1 protein levels in tumors. Another candidate variant, rs78540526, lies in the same enhancer element. Risk association signal 2, rs75915166, creates a GATA3 binding site within a silencer element. Chromatin conformation studies demonstrate that these enhancer and silencer elements interact with each other and with their likely target gene, CCND1. © 2013 The American Society of Human Genetics.

The activity of the immunoglobulin 3' enhancer is restricted to the late stages of B lymphoid development. Here we further examine the molecular basis for the temporally restricted activity of the B-lymphoid IgH 3' enhancer. We demonstrate that a binding site (E5 site) for the E47 and/or E12 proteins is functionally important for enhancer activity. The multimerized E5 site acts as a B cell-specific enhancer and, when assayed in COS cells, can be transactivated by E47/E12 proteins. This transactivation in COS cells, as well as the activity of the full length 3' enhancer in plasma cells, can be repressed by overexpression of the dominant negative nuclear regulator Id3. When examining the tissue distribution of Id3 in murine cell lines, we find that Id3 is expressed throughout the pre-B and B cell stages, but is down-regulated at the plasma cell stage. Thus, Id3 may contribute to the temporal regulation of the IgH 3' enhancer.

The activity of the IgH (E mu) enhancer in the T lymphocyte lineage has been investigated using both transgenic mice and transfection studies. Thymocyte fractionation experiments indicate that a transgene consisting of the bacterial chloramphenicol acetyl transferase (CAT) gene, linked to E mu and the SV40 early promoter (E mu-CAT), is expressed only in thymocytes with a mature medullary phenotype and not in immature cells. Transfection of this same construct into two thymoma cell lines representing different stages of thymocyte development mimics the pattern of activity observed in vivo. Further transfection experiments suggest that this pattern of expression might be attributed to the differential activity of the E2E3 and octanucleotide motifs of E mu during development. In contrast, an Ig lambda transgene (linked to E mu and an Ig V lambda promoter) is expressed in the majority of thymocytes. We envisage that the different patterns of expression of the two transgenes reflect interactions between their respective promoters and the factors which are bound to E mu at different stages of thymocyte development. Although differing in their pattern of expression within the thymus, the two transgenes share the property of extinction in peripheral T lymphocytes. These results indicate that the expression of E mu-linked transgenes in the thymus cannot simply be explained by activation of the enhancer in a lymphoid progenitor cell prior to B/T lineage divergence. Rather, the enhancer (or components of it) must be independently activated (and inactivated) during T lymphocyte development. Furthermore, this activity is consistent with the developmental timing of Ig DH-JH rearrangements in these cells.

Recent genome-wide association studies have identified a breast cancer susceptibility locus on 16q12 with an unknown biological basis. We used a set of single nucleotide polymorphism (SNP) markers to generate a fine-scale map and narrowed the region of association to a 133 kb DNA segment containing the largely uncharacterized hypothetical gene LOC643714, a short intergenic region and the 5' end of TOX3. Re-sequencing this segment in European subjects identified 293 common polymorphisms, including a set of 26 highly correlated candidate causal variants. By evaluation of these SNPs in five breast cancer case-control studies involving more than 23 000 subjects from populations of European and Southeast Asian ancestry, all but 14 variants could be excluded at odds of