Single-cell Gene Expression Analysis Research Articles

Increased Expression of Interleukin-1 Receptor Characterizes Anti-estrogen-Resistant ALDH+ Breast Cancer Stem Cells.

SummaryEstrogen-receptor-positive breast tumors are treated with anti-estrogen (AE) therapies but frequently develop resistance. Cancer stem cells (CSCs) with high aldehyde dehydrogenase activity (ALDH+ cells) are enriched following AE treatment. Here, we show that the interleukin-1β (IL-1β) signaling pathway is activated in ALDH+ cells, and data from single cells reveals that AE treatment selects for IL-1 receptor (IL1R1)-expressing ALDH+ cells. Importantly, CSC activity is reduced by an IL1R1 inhibitor in AE-resistant models. Moreover, IL1R1 expression is increased in the tumors of patients treated with AE therapy and predicts treatment failure. Single-cell gene expression analysis revealed that at least two subpopulations exist within the ALDH+ population, one proliferative and one quiescent. Following AE therapy the quiescent population is expanded, which suggests CSC dormancy as an adaptive strategy that facilitates treatment resistance. Targeting of ALDH+IL1R1+ cells merits testing as a strategy to combat AE resistance in patients with residual disease.

Bulk and single-cell gene expression analyses reveal aging human choriocapillaris has pro-inflammatory phenotype

The human choroidal vasculature is subject to age-related structural and gene expression changes implicated in age-related macular degeneration (AMD). In this study, we performed both bulk and single-cell RNA sequencing on infant (n = 4 for bulk experiments, n = 2 for single-cell experiments) and adult (n = 13 for bulk experiments, n = 6 for single-cell experiments) human donors to characterize how choroidal gene expression changes with age. Differential expression analysis revealed that aged choroidal samples were enriched in genes encoding pro-inflammatory transcription factors and leukocyte transendothelial cell migration adhesion proteins. Such genes were observed to be differentially expressed specifically within choroidal endothelial cells at the single-cell level. Immunohistochemistry experiments support transcriptional findings that CD34 is elevated in infant choriocapillaris endothelial cells while ICAM-1 is enriched in adults. These results suggest several potential drivers of the pro-inflammatory vascular phenotype observed with advancing age.

Single-cell RNA-seq analysis of human CSF microglia and myeloid cells in neuroinflammation

ObjectiveTo identify and characterize myeloid cell populations within the CSF of patients with MS and anti-myelin oligodendrocyte glycoprotein (MOG) disorder by high-resolution single-cell gene expression analysis.MethodsSingle-cell RNA sequencing (scRNA-seq) was used to profile individual cells of CSF and blood from 2 subjects with relapsing-remitting MS (RRMS) and one with anti-MOG disorder. Publicly available scRNA-seq data from the blood and CSF of 2 subjects with HIV were also analyzed. An informatics pipeline was used to cluster cell populations by transcriptomic profiling. Based on gene expression by CSF myeloid cells, a flow cytometry panel was devised to examine myeloid cell populations from the CSF of 11 additional subjects, including individuals with RRMS, anti-MOG disorder, and control subjects without inflammatory demyelination.ResultsCommon myeloid populations were identified within the CSF of subjects with RRMS, anti-MOG disorder, and HIV. These included monocytes, conventional and plasmacytoid dendritic cells, and cells with a transcriptomic signature matching microglia. Microglia could be discriminated from other myeloid cell populations in the CSF by flow cytometry.ConclusionsHigh-resolution single-cell gene expression analysis clearly distinguishes distinct myeloid cell types present within the CSF of subjects with neuroinflammation. A population of microglia exists within the human CSF, which is detectable by surface protein expression. The function of these cells during immunity and disease requires further investigation.

High-throughput multi-dimensional single antigen-specific T cell profiling

Abstract Recent advances in single cell analysis has made it possible to interrogate the transcriptome of single cells in high-throughput manner. However, multi-dimensional profiling of antigen-specific T cells at single cell level for hundreds of antigens has not been realized, which becomes a bottleneck to the understanding of these important cells in immunological diseases. Here we describe a new technology, TetTCR-SeqHD, that extends our previously developed technology, TetTCR-Seq for Tetramer associated TCR sequencing, to additional two dimensions, single cell gene expression and phenotype. In addition, we expand the throughput to tens of thousands of single T cells. We achieved these by using DNA-barcoded pMHC tetramers and antibodies in combination with single cell gene expression analysis, thus, paralleling all read-outs using next-gen sequencing. To validate the technology, we have used previously generated polyclonal T cells to 6 different antigens with known TCR sequences. This analysis showed that we could achieve less than 5% on false discovery rate. We then used TetTCR-SeqHD on a panel of type 1 diabetes related peptides and showed unexpected features of antigen-specific T cells that separate healthy individuals from type 1 diabetes patients. In addition, TCR transgenics T cells confirmed their antigen specificity. In summary, TetTCR-SeqHD enables us to analyze antigen-specific T cells with an extraordinary depth and breadth, which enables discoveries directly in humans. TetTCR-SeqHD would be a useful technology for studying many different diseases where screening for T cells antigens and deep profiling of antigen-specific T cells are highly sought after.

Expression of RUNX1-ETO Rapidly Alters the Chromatin Landscape and Growth of Early Human Myeloid Precursor Cells

SummaryAcute myeloid leukemia (AML) is a hematopoietic malignancy caused by recurrent mutations in genes encoding transcriptional, chromatin, and/or signaling regulators. The t(8;21) translocation generates the aberrant transcription factor RUNX1-ETO (RUNX1-RUNX1T1), which by itself is insufficient to cause disease. t(8;21) AML patients show extensive chromatin reprogramming and have acquired additional mutations. Therefore, the genomic and developmental effects directly and solely attributable to RUNX1-ETO expression are unclear. To address this, we employ a human embryonic stem cell differentiation system capable of forming definitive myeloid progenitor cells to express RUNX1-ETO in an inducible fashion. Induction of RUNX1-ETO causes extensive chromatin reprogramming by interfering with RUNX1 binding, blocks differentiation, and arrests cellular growth, whereby growth arrest is reversible following RUNX1-ETO removal. Single-cell gene expression analyses show that RUNX1-ETO induction alters the differentiation of early myeloid progenitors, but not of other progenitor types, indicating that oncoprotein-mediated transcriptional reprogramming is highly target cell specific.

Memory of Divisional History Directs the Continuous Process of Primitive Hematopoietic Lineage Commitment.

SummaryHematopoietic stem cells (HSCs) exist in a dormant state and progressively lose regenerative potency as they undergo successive divisions. Why this functional decline occurs and how this information is encoded is unclear. To better understand how this information is stored, we performed RNA sequencing on HSC populations differing only in their divisional history. Comparative analysis revealed that genes upregulated with divisions are enriched for lineage genes and regulated by cell-cycle-associated transcription factors, suggesting that proliferation itself drives lineage priming. Downregulated genes are, however, associated with an HSC signature and targeted by the Polycomb Repressive Complex 2 (PRC2). The PRC2 catalytic subunits Ezh1 and Ezh2 promote and suppress the HSC state, respectively, and successive divisions cause a switch from Ezh1 to Ezh2 dominance. We propose that cell divisions drive lineage priming and Ezh2 accumulation, which represses HSC signature genes to consolidate information on divisional history into memory.

Spatiotemporal single-cell analysis of gene expression in the mouse suprachiasmatic nucleus.

Mammalian circadian behaviors are orchestrated by the suprachiasmatic nucleus (SCN) in the ventral hypothalamus, but the number of SCN cell types and their functional roles remain unclear. We have used single-cell RNA-sequencing to identify the basic cell types in the mouse SCN and to characterize their circadian and light-induced gene expression patterns. We identified eight major cell types, with each type displaying a specific pattern of circadian gene expression. Five SCN neuronal subtypes, each with specific combinations of markers, differ in their spatial distribution, circadian rhythmicity and light responsiveness. Through a complete three-dimensional reconstruction of the mouse SCN at single-cell resolution, we obtained a standardized SCN atlas containing the spatial distribution of these subtypes and gene expression. Furthermore, we observed heterogeneous circadian gene expression between SCN neuron subtypes. Such a spatiotemporal pattern of gene regulation within the SCN may have an important function in the circadian pacemaker.

Single-Cell Gene Expression Analyses Reveal Distinct Self-Renewing and Proliferating Subsets in the Leukemia Stem Cell Compartment in Acute Myeloid Leukemia.

Standard chemotherapy for acute myeloid leukemia (AML) targets proliferative cells and efficiently induces complete remission; however, many patients relapse and die of their disease. Relapse is caused by leukemia stem cells (LSC), the cells with self-renewal capacity. Self-renewal and proliferation are separate functions in normal hematopoietic stem cells (HSC) in steady-state conditions. If these functions are also separate functions in LSCs, then antiproliferative therapies may fail to target self-renewal, allowing for relapse. We investigated whether proliferation and self-renewal are separate functions in LSCs as they often are in HSCs. Distinct transcriptional profiles within LSCs of Mll-AF9/NRASG12V murine AML were identified using single-cell RNA sequencing. Single-cell qPCR revealed that these genes were also differentially expressed in primary human LSCs and normal human HSPCs. A smaller subset of these genes was upregulated in LSCs relative to HSPCs; this subset of genes constitutes "LSC-specific" genes in human AML. To assess the differences between these profiles, we identified cell surface markers, CD69 and CD36, whose genes were differentially expressed between these profiles. In vivo mouse reconstitution assays resealed that only CD69High LSCs were capable of self-renewal and were poorly proliferative. In contrast, CD36High LSCs were unable to transplant leukemia but were highly proliferative. These data demonstrate that the transcriptional foundations of self-renewal and proliferation are distinct in LSCs as they often are in normal stem cells and suggest that therapeutic strategies that target self-renewal, in addition to proliferation, are critical to prevent relapse and improve survival in AML. SIGNIFICANCE: These findings define and functionally validate a self-renewal gene profile of leukemia stem cells at the single-cell level and demonstrate that self-renewal and proliferation are distinct in AML. GRAPHICAL ABSTRACT: http://cancerres.aacrjournals.org/content/canres/80/3/458/F1.large.jpg.

Treg Heterogeneity, Function, and Homeostasis.

T-regulatory cells (Tregs) represent a unique subpopulation of helper T-cells by maintaining immune equilibrium using various mechanisms. The role of T-cell receptors (TCR) in providing homeostasis and activation of conventional T-cells is well-known; however, for Tregs, this area is understudied. In the last two decades, evidence has accumulated to confirm the importance of the TCR in Treg homeostasis and antigen-specific immune response regulation. In this review, we describe the current view of Treg subset heterogeneity, homeostasis and function in the context of TCR involvement. Recent studies of the TCR repertoire of Tregs, combined with single-cell gene expression analysis, revealed the importance of TCR specificity in shaping Treg phenotype diversity, their functions and homeostatic maintenance in various tissues. We propose that Tregs, like conventional T-helper cells, act to a great extent in an antigen-specific manner, which is provided by a specific distribution of Tregs in niches.

Increased mesodermal and mesendodermal populations by BMP4 treatment facilitates human iPSC line differentiation into a cardiac lineage

Human induced pluripotent stem cell-derived cardiomyocytes (hiPS-CMs) have attracted attention as a novel tool for drug safety screening and several differentiation protocols of hiPSC lines into cardiomyocytes have been reported; the standardization of these protocols will expand their applications for safety assessments such as "clinical safety trial-on-dish". Bone morphogenetic protein 4 (BMP4) is an important factor in promoting mesoderm differentiation and BMP4 treatment has been used at the early stage of cardiac differentiation into different hiPSCs. In the present study, we evaluated the effects of BMP4 treatment at the early stage of cardiac differentiation. We performed gene expression profiling of the germ layer during mesoderm differentiation of hiPSCs derived from three different donors. The expression of T (a mesoderm marker) and GATA6 (an endoderm marker) increased and that of PAX6 (a neuroectoderm marker) decreased in pooled embryoid bodies (EBs) after BMP4 treatment. Single-cell gene expression analysis revealed that mesodermal and mesendodermal populations increased in EBs derived from 253G1. Finally, BMP4 treatment increased mesodermal and mesendodermal populations compared with that without BMP4 in two other hiPSC lines, confirming the reproducibility of multiple hiPSC lines. Thus, our results suggest that BMP4 treatment increases mesodermal and mesendodermal populations at the early stage of cardiac differentiation in different hiPSC lines.

Meta-analysis of Caenorhabditis elegans single-cell developmental data reveals multi-frequency oscillation in gene activation.

MotivationThe advent of in vivo automated techniques for single-cell lineaging, sequencing and analysis of gene expression has begun to dramatically increase our understanding of organismal development. We applied novel meta-analysis and visualization techniques to the EPIC single-cell-resolution developmental gene expression dataset for Caenorhabditis elegans from Bao, Murray, Waterston et al. to gain insights into regulatory mechanisms governing the timing of development.ResultsOur meta-analysis of the EPIC dataset revealed that a simple linear combination of the expression levels of the developmental genes is strongly correlated with the developmental age of the organism, irrespective of the cell division rate of different cell lineages. We uncovered a pattern of collective sinusoidal oscillation in gene activation, in multiple dominant frequencies and in multiple orthogonal axes of gene expression, pointing to the existence of a coordinated, multi-frequency global timing mechanism. We developed a novel method based on Fisher’s Discriminant Analysis to identify gene expression weightings that maximally separate traits of interest, and found that remarkably, simple linear gene expression weightings are capable of producing sinusoidal oscillations of any frequency and phase, adding to the growing body of evidence that oscillatory mechanisms likely play an important role in the timing of development. We cross-linked EPIC with gene ontology and anatomy ontology terms, employing Fisher’s Discriminant Analysis methods to identify previously unknown positive and negative genetic contributions to developmental processes and cell phenotypes. This meta-analysis demonstrates new evidence for direct linear and/or sinusoidal mechanisms regulating the timing of development. We uncovered a number of previously unknown positive and negative correlations between developmental genes and developmental processes or cell phenotypes. Our results highlight both the continued relevance of the EPIC technique, and the value of meta-analysis of previously published results. The presented analysis and visualization techniques are broadly applicable across developmental and systems biology.Availability and implementationAnalysis software available upon request.Supplementary information Supplementary data are available at Bioinformatics online.

Characterization of Leukemic Stem Cells Heterogeneity in Chronic Myeloid Leukemia

In chronic myeloid leukemia (CML), a rare subset of leukemic stem cells (LSC) persists in patients responding to conventional tyrosine kinase inhibitor (TKI) therapy. The failure to eradicate these LSCs results in indefinite therapy dependence and a risk of leukemic relapse. However, the conventional LSC compartment (Lin-CD34+CD38-) is highly heterogeneous where only a subpopulation is believed to be functional, TKI-insensitive LSCs. Previously, using single-cell gene expression analysis we characterized the heterogeneity within the LSC population (Lin-CD34+CD38-) in CML patients using a selected panel of 96 primers. Interestingly, by comparing LSC heterogeneity at diagnosis with the heterogeneity following 3 months of TKI therapy we uncovered a therapy-insensitive, quiescent subpopulation, which could be isolated at high-purity using a combination of the surface markers: Lin-CD34+CD38-CD45RA-cKIT-CD26+ (Warfvinge, Geironson, Sommarin et al., 2017). Here, we expand the single-cell analysis of CML LSC populations to include combined immunophenotype-/RNA sequencing analysis (CITE-seq). CITE-seq allows for unbiased, further in-depth transcriptome analysis as wells as immunophenotypic characterization by pre-staining cells with a panel of DNA-barcoded antibodies prior to sequencing. DNA-barcoded antibodies convert the protein expression into readable sequences through unique oligo-conjugates as identifiers. Using CITE-seq with a panel of 44 distinct surface markers designed to immunophenotypically differentiate between stem/progenitors cells and leukemic clones we simultaneously characterize the molecular and immunophenotypic heterogeneity within Lin-CD34+/Lin-CD34+CD38- CML stem/progenitor compartment at diagnosis. Additionally by comparing the LSCs transcriptome from patients with different therapeutic outcome after 12 months of therapy we describe how differences in heterogeneity and the presence of immunophenotypic therapy-insensitive LSCs at diagnosis (Lin-CD34+CD38-CD45RA-cKIT-CD26+) contribute to therapy response. Disclosures Richter: Novartis: Consultancy; Pfizer: Consultancy, Research Funding.

SB-3 Combining live imaging and single-cell whole gene expression analysis by developing an automated cell picking system

SB-3 Combining live imaging and single-cell whole gene expression analysis by developing an automated cell picking system

Using single-cell transcriptomics to understand functional states and interactions in microbial eukaryotes.

Understanding the diversity and evolution of eukaryotic microorganisms remains one of the major challenges of modern biology. In recent years, we have advanced in the discovery and phylogenetic placement of new eukaryotic species and lineages, which in turn completely transformed our view on the eukaryotic tree of life. But we remain ignorant of the life cycles, physiology and cellular states of most of these microbial eukaryotes, as well as of their interactions with other organisms. Here, we discuss how high-throughput genome-wide gene expression analysis of eukaryotic single cells can shed light on protist biology. First, we review different single-cell transcriptomics methodologies with particular focus on microbial eukaryote applications. Then, we discuss single-cell gene expression analysis of protists in culture and what can be learnt from these approaches. Finally, we envision the application of single-cell transcriptomics to protist communities to interrogate not only community components, but also the gene expression signatures of distinct cellular and physiological states, as well as the transcriptional dynamics of interspecific interactions. Overall, we argue that single-cell transcriptomics can significantly contribute to our understanding of the biology of microbial eukaryotes. This article is part of a discussion meeting issue 'Single cell ecology'.

Tensor Decomposition-Based Unsupervised Feature Extraction Applied to Single-Cell Gene Expression Analysis.

Although single-cell RNA sequencing (scRNA-seq) technology is newly invented and a promising one, but because of lack of enough information that labels individual cells, it is hard to interpret the obtained gene expression of each cell. Because of insufficient information available, unsupervised clustering, for example, t-distributed stochastic neighbor embedding and uniform manifold approximation and projection, is usually employed to obtain low-dimensional embedding that can help to understand cell–cell relationship. One possible drawback of this strategy is that the outcome is highly dependent upon genes selected for the usage of clustering. In order to fulfill this requirement, there are many methods that performed unsupervised gene selection. In this study, a tensor decomposition (TD)-based unsupervised feature extraction (FE) was applied to the integration of two scRNA-seq expression profiles that measure human and mouse midbrain development. TD-based unsupervised FE could select not only coincident genes between human and mouse but also biologically reliable genes. Coincidence between two species as well as biological reliability of selected genes is increased compared with that using principal component analysis (PCA)-based FE applied to the same data set in the previous study. Since PCA-based unsupervised FE outperformed the other three popular unsupervised gene selection methods, highly variable genes, bimodal genes, and dpFeature, TD-based unsupervised FE can do so as well. In addition to this, 10 transcription factors (TFs) that might regulate selected genes and might contribute to midbrain development were identified. These 10 TFs, BHLHE40, EGR1, GABPA, IRF3, PPARG, REST, RFX5, STAT3, TCF7L2, and ZBTB33, were previously reported to be related to brain functions and diseases. TD-based unsupervised FE is a promising method to integrate two scRNA-seq profiles effectively.

A Microwell Cell Capture Device Reveals Variable Response to Dobutamine in Isolated Cardiomyocytes

Isolated ventricular cardiomyocytes exhibit substantial cell-to-cell variability, even when obtained from the same small volume of myocardium. In this study, we investigated the possibility that cardiomyocyte responses to β-adrenergic stimulus are also highly heterogeneous. To achieve the throughput and measurement duration desired for these experiments, we designed and validated a novel microwell system that immobilizes and uniformly orients isolated adult cardiomyocytes. In this configuration, detailed drug responses of dozens of cells can be followed for intervals exceeding 1 h. At the conclusion of an experiment, specific cells can also be harvested via a precision aspirator for single-cell gene expression profiling. Using this system, we followed changes in Ca2+ signaling and contractility of individual cells under sustained application of either dobutamine or omecamtiv mecarbil. Both compounds increased average cardiomyocyte contractility over the course of an hour, but responses of individual cells to dobutamine were significantly more variable. Surprisingly, some dobutamine-treated cardiomyocytes augmented Ca2+ release without increasing contractility. Other cells responded with increased contractility despite unchanged Ca2+ release. Single-cell gene expression analysis revealed significant co-expression of β-adrenergic pathway genes PKA regulatory subunit type I, PKA regulatory subunit type II, and Ca2+/calmodulin-dependent protein kinase II across cardiomyocytes. Other data supported a connection between the effects of dobutamine on relaxation rate and the expression of protein phosphatase 2. These findings suggest that variable drug responses among cells are not merely experimental artifacts. By enabling direct comparison of the functional behavior of an individual cell and the genes it expresses, this new system constitutes a unique tool for interrogating cardiomyocyte drug responses and discovering the genes that modulate them.

RUNX1-ETO RAPIDLY REPROGRAMS THE EPIGENETIC LANDSCAPE IN HUMAN EMBRYONIC DEFINITIVE HAEMATOPOIETIC PROGENITORS

Acute Myeloid Leukaemia results from the successive acquisition of multiple genetic aberrations. Therefore, the sole study of patient samples, which represent the end-point of this process, hinders the understanding of the effect of each individual mutation on reprogramming the epigenome. In order to dissect tumour development, we generated human Embryonic Stem Cell lines carrying an inducible version of RUNX1-ETO – the product of the t(8;21) translocation – which we subsequently differentiated into haematopoietic cells. This mutation is observed in utero and in cells of long-term remission patients, where it consists a reservoir for relapse. When induced at physiological levels at early embryonic stages, RUNX1-ETO disrupts vasculogenesis and prevents blood formation. However, induction in already formed progenitors (i) blocks differentiation at an immature stage, (ii) reduces colony-forming capacity, (iii) induces a reversible cell-cycle block caused by the down-regulation of haematopoietic, cell cycle and replication genes, and (iv) induces the up-regulation of multiple signalling genes. Genome wide chromatin immunoprecipitation (ChIP) and ATAC-Seq analyses show that RUNX1-ETO displaces RUNX1 from a subset of binding sites, resulting in reduced chromatin accessibility and gene silencing. RUNX1-ETO induction dramatically alters gene expression with some changes appearing to be irreversible. Single cell gene expression analyses show that induction of RUNX1-ETO in the CD45+CD34+RUNX1C+ population results in the emergence of a new cell population unrelated to normal cells that may represent the pre-leukaemic clone. Our data are consistent with the idea that RUNX1-ETO establishes a precondition for leukaemic transformation by maintaining a reservoir of quiescent pre-leukaemic progenitors with susceptibility to expand upon acquisition of an additional oncogenic event.

CHARACTERIZATION OF THERAPY-INSENSITIVE LSC SUBPOPULATIONS IN CHRONIC MYELOID LEUKEMIA

In chronic myeloid leukemia (CML), a rare subset of leukemic stem cells (LSC) persists in patients responding to conventional tyrosine kinase inhibitor (TKI) therapy. The failure to eradicate these LSCs results in indefinite therapy dependence and a risk of leukemic relapse. However, the conventional LSC compartment (Lin-CD34+CD38-) is highly heterogeneous where only a sub population is believed to be functional, TKI-insensitive LSCs. Previously, using single-cell gene expression analysis we characterized the heterogeneity within the LSC population (Lin-CD34+CD38-) in CML patients. Interestingly, by comparing LSC heterogeneity at diagnosis with the heterogeneity following 3 months of TKI therapy we uncovered a therapy-insensitive, quiescent LSC-subpopulation, which could be isolated at high-purity using a combination of the surface markers: Lin-CD34+CD38-CD45RA-cKIT-CD26+ (Warfvinge, Geironson, Sommarin et al., 2017). Here, we expand the single-cell analysis of CML LSC populations to include combined immunophenotype-/RNA sequencing analysis (CITE-seq). By comparing LSCs from patients with different therapeutic outcome we describe how LSC heterogeneity at diagnosis correlate with long-term therapy response. We furthermore evaluate the potential for humanized ossicle-niche xenotransplantation (Reinisch et al., 2016) as a model system for functional characterization of LSC subpopulations in vivo.

GENE REGULATORY NETWORK IN THE CONTEXT OF GENOME INSTABILITY IN HUMAN EMBRYO DEVELOPMENT

Genomic instability is common in preimplantation embryos of a variety of species. Up to 80% of human cleavage stage embryos following in vitro fertilization acquire genetic mosaicism, but some may develop to normal individuals after uterine transfer. The causative mechanisms of this chromosome instability as well as the impact of acquired genetic anomalies on embryo development are not understood and remain speculative. Using novel methodologies for single-cell genome-plus-transcriptome sequencing (G&T-seq) of all cells of human preimplantation embryos from the fertilized egg to the blastocyst stage, we disclose insight in the gene regulatory network of preimplantation embryo development and how this is impacted by aneuploidy. Single-cell DNA CNV analysis revealed frequent missegregations of whole chromosomes as well as segmental rearrangements in all embryonic cell stages. Using these single-cell genomic profiles in combination with cell cleavage imaging of the developing embryos, we could deduce the origin of the abnormalities and construct cell lineages. Single-cell gene expression analysis of the same cells, classified the cells according to the embryonic developmental stage, the expression activation of the embryonic genome and disclosed the functional impact of acquired numerical and structural chromosome aberrations on development of the human embryo. Furthermore, the data reveals which genetic anomalies contribute to the epiblast cells in the blastocyst that provide the ectodermal, mesodermal and endodermal cell lineages, the building blocks of our organs. Genomic instability is common in preimplantation embryos of a variety of species. Up to 80% of human cleavage stage embryos following in vitro fertilization acquire genetic mosaicism, but some may develop to normal individuals after uterine transfer. The causative mechanisms of this chromosome instability as well as the impact of acquired genetic anomalies on embryo development are not understood and remain speculative. Using novel methodologies for single-cell genome-plus-transcriptome sequencing (G&T-seq) of all cells of human preimplantation embryos from the fertilized egg to the blastocyst stage, we disclose insight in the gene regulatory network of preimplantation embryo development and how this is impacted by aneuploidy. Single-cell DNA CNV analysis revealed frequent missegregations of whole chromosomes as well as segmental rearrangements in all embryonic cell stages. Using these single-cell genomic profiles in combination with cell cleavage imaging of the developing embryos, we could deduce the origin of the abnormalities and construct cell lineages. Single-cell gene expression analysis of the same cells, classified the cells according to the embryonic developmental stage, the expression activation of the embryonic genome and disclosed the functional impact of acquired numerical and structural chromosome aberrations on development of the human embryo. Furthermore, the data reveals which genetic anomalies contribute to the epiblast cells in the blastocyst that provide the ectodermal, mesodermal and endodermal cell lineages, the building blocks of our organs.

Identification of Host Biomarkers of Epstein-Barr Virus Latency IIb and Latency III.

Deciphering the molecular pathogenesis of virally induced cancers is challenging due, in part, to the heterogeneity of both viral gene expression and host gene expression. Epstein-Barr virus (EBV) is a ubiquitous herpesvirus prevalent in B-cell lymphomas of immune-suppressed individuals. EBV infection of primary human B cells leads to their immortalization into lymphoblastoid cell lines (LCLs), serving as a model of these lymphomas. In previous studies, reports from our laboratory have described a temporal model for immortalization with an initial phase characterized by expression of Epstein-Barr nuclear antigens (EBNAs), high levels of c-Myc activity, and hyperproliferation in the absence of the latent membrane proteins (LMPs), called latency IIb. This is followed by the long-term outgrowth of LCLs expressing the EBNAs along with the LMPs, particularly NFκB-activating LMP1, defining latency III. However, LCLs express a broad distribution of LMP1 such that a subset of these cells express LMP1 at levels similar to those seen in latency IIb, making it difficult to distinguish these two latency states. In this study, we performed mRNA sequencing (mRNA-Seq) on early EBV-infected latency IIb cells and latency III LCLs sorted by NFκB activity. We found that latency IIb transcriptomes clustered independently from latency III independently of NFκB. We identified and validated mRNAs defining these latency states. Indeed, we were able to distinguish latency IIb cells from LCLs expressing low levels of LMP1 using multiplex RNA-fluorescence in situ hybridization (RNA-FISH) targeting EBV EBNA2 or LMP1 and human CCR7 or MGST1 This report defines latency IIb as a bona fide latency state independent from latency III and identifies biomarkers for understanding EBV-associated tumor heterogeneity.IMPORTANCE EBV is a ubiquitous pathogen, with >95% of adults harboring a life-long latent infection in memory B cells. In immunocompromised individuals, latent EBV infection can result in lymphoma. The established expression profile of these lymphomas is latency III, which includes expression of all latency genes. However, single-cell analysis of EBV latent gene expression in these lymphomas suggests heterogeneity where most cells express the transcription factor, EBNA2, and only a fraction of the cells express membrane protein LMP1. Our work describes an early phase after infection where the EBNAs are expressed without LMP1, called latency IIb. However, LMP1 levels within latency III vary widely, making these states hard to discriminate. This may have important implications for therapeutic responses. It is crucial to distinguish these states to understand the molecular pathogenesis of these lymphomas. Ultimately, better tools to understand the heterogeneity of these cancers will support more-efficacious therapies in the future.