Single-cell RNA-seq reveals transcriptomic heterogeneity mediated by host-pathogen dynamics in lymphoblastoid cell lines
- Duke University School of Medicine, United States
- University of Massachusetts Medical School, United States
- Brown University, United States
Abstract
Lymphoblastoid Cell Lines (LCLs) are generated by transforming primary B cells with Epstein-Barr Virus (EBV) and are used extensively as model systems in viral oncology, immunology, and human genetics research. In this study, we characterized single-cell transcriptomic profiles of five LCLs and present a simple discrete-time simulation to explore the influence of stochasticity on LCL clonal evolution. Single-cell RNA sequencing (scRNA-seq) revealed substantial phenotypic heterogeneity within and across LCLs with respect to immunoglobulin isotype; virus-modulated host pathways involved in survival, activation, and differentiation; viral replication state; and oxidative stress. This heterogeneity is likely attributable to intrinsic variance in primary B cells and host-pathogen dynamics. Stochastic simulations demonstrate that initial primary cell heterogeneity, random sampling, time in culture, and even mild differences in phenotype-specific fitness can contribute substantially to dynamic diversity in populations of nominally clonal cells.
Data availability
Raw sequencing data for the three previously unpublished samples (LCL_777_B958, LCL_777_M81, and LCL_461_B958) are deposited in the NCBI Sequence Read Archive (SRA) and can be accessed along with processed data from the NCBI Gene Expression Omnibus (GEO, Series Accession: GSE158275).
-
Single-cell characterization of transcriptomic heterogeneity in lymphoblastoid cell linesNCBI Gene Expression Omnibus, GSE158275.
-
Single cell RNA sequencing of lymphoblastoid cell lines of European and African ancestriesNCBI Gene Expression Omnibus, GSE126321.
Article and author information
Author details
Funding
National Institute of Dental and Craniofacial Research (R01-DE025994)
- Micah A Luftig
National Cancer Institute (T32-CA009111)
- Elliott D SoRelle
- Joanne Dai
- Micah A Luftig
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Copyright
© 2021, SoRelle et al.
This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 7,128
- views
-
- 619
- downloads
-
- 32
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
A two-part list of links to download the article, or parts of the article, in various formats.
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Single-cell RNA-seq reveals transcriptomic heterogeneity mediated by host-pathogen dynamics in lymphoblastoid cell lines
eLife 10:e62586.
https://doi.org/10.7554/eLife.62586
Further reading
-
- Genetics and Genomics
- Immunology and Inflammation
PIK3R1 encodes three regulatory subunits of class IA phosphoinositide 3-kinase (PI3K), each associating with any of three catalytic subunits, namely p110α, p110β, or p110δ. Constitutional PIK3R1 mutations cause diseases with a genotype-phenotype relationship not yet fully explained: heterozygous loss-of-function mutations cause SHORT syndrome, featuring insulin resistance and short stature attributed to reduced p110α function, while heterozygous activating mutations cause immunodeficiency, attributed to p110δ activation and known as APDS2. Surprisingly, APDS2 patients do not show features of p110α hyperactivation, but do commonly have SHORT syndrome-like features, suggesting p110α hypofunction. We sought to investigate this. In dermal fibroblasts from an APDS2 patient, we found no increased PI3K signalling, with p110δ expression markedly reduced. In preadipocytes, the APDS2 variant was potently dominant negative, associating with Irs1 and Irs2 but failing to heterodimerise with p110α. This attenuation of p110α signalling by a p110δ-activating PIK3R1 variant potentially explains co-incidence of gain-of-function and loss-of-function PIK3R1 phenotypes.
