1. Immunology and Inflammation
  2. Microbiology and Infectious Disease

Single-cell RNA-seq reveals transcriptomic heterogeneity mediated by host-pathogen dynamics in lymphoblastoid cell lines

  1. Elliott D SoRelle
  2. Joanne Dai
  3. Emmanuela N Bonglack
  4. Emma M Heckenberg
  5. Jeffrey Y Zhou
  6. Stephanie N Giamberardino
  7. Jeffrey A Bailey
  8. Simon G Gregory
  9. Cliburn Chan
  10. Micah A Luftig  Is a corresponding author
  1. Duke University School of Medicine, United States
  2. University of Massachusetts Medical School, United States
  3. Brown University, United States
https://doi.org/10.7554/eLife.62586
Share this article
Cite this article
  1. Elliott D SoRelle
  2. Joanne Dai
  3. Emmanuela N Bonglack
  4. Emma M Heckenberg
  5. Jeffrey Y Zhou
  6. Stephanie N Giamberardino
  7. Jeffrey A Bailey
  8. Simon G Gregory
  9. Cliburn Chan
  10. Micah A Luftig
(2021)
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
  2. Download BibTeX
  3. Download .RIS

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).

The following data sets were generated
The following previously published data sets were used

Article and author information

Author details

  1. Elliott D SoRelle

    Molecular Genetics and Microbiology, Duke University School of Medicine, Durham, United States
    Competing interests
    No competing interests declared.

  2. Joanne Dai

    Department of Molecular Genetics and Microbiology, Duke University School of Medicine, Durham, United States
    Competing interests
    Joanne Dai, Joanne Dai is affiliated with Amgen Inc. The author has no financial interests to declare..
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-9879-4704

  3. Emmanuela N Bonglack

    Molecular Genetics and Microbiology, Pharmacology and Cancer Biology, Duke University School of Medicine, Durham, United States
    Competing interests
    No competing interests declared.

  4. Emma M Heckenberg

    Molecular Genetics and Microbiology, Duke University School of Medicine, Durham, United States
    Competing interests
    No competing interests declared.

  5. Jeffrey Y Zhou

    Medicine, University of Massachusetts Medical School, Worcester, United States
    Competing interests
    No competing interests declared.

  6. Stephanie N Giamberardino

    Department of Neurology, Duke University School of Medicine, Durham, United States
    Competing interests
    No competing interests declared.

  7. Jeffrey A Bailey

    Department of Pathology and Laboratory Medicine, Brown University, Providence, United States
    Competing interests
    No competing interests declared.

  8. Simon G Gregory

    Department of Neurology, Duke University School of Medicine, Durham, United States
    Competing interests
    No competing interests declared.

  9. Cliburn Chan

    Department of Biostatistics and Bioinformatics, Duke University School of Medicine, Durham, United States
    Competing interests
    No competing interests declared.

  10. Micah A Luftig

    Department of Molecular Genetics and Microbiology, Duke University School of Medicine, Durham, United States
    For correspondence
    micah.luftig@duke.edu
    Competing interests
    No competing interests declared.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-2964-1907

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.

Reviewing Editor

  1. Melanie M Brinkmann, Technische Universität Braunschweig, Germany

Version history

  1. Received: August 29, 2020
  2. Accepted: January 26, 2021
  3. Accepted Manuscript published: January 27, 2021 (version 1)
  4. Version of Record published: February 6, 2021 (version 2)
  5. Version of Record updated: November 11, 2021 (version 3)

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

  • 6,611
    views
  • 574
    downloads
  • 18
    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)

  1. Elliott D SoRelle
  2. Joanne Dai
  3. Emmanuela N Bonglack
  4. Emma M Heckenberg
  5. Jeffrey Y Zhou
  6. Stephanie N Giamberardino
  7. Jeffrey A Bailey
  8. Simon G Gregory
  9. Cliburn Chan
  10. Micah A Luftig
(2021)
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

Share this article

https://doi.org/10.7554/eLife.62586

Categories and tags

Research organisms

Further reading

    1. Evolutionary Biology
    2. Immunology and Inflammation

    The CD4 transmembrane GGXXG and juxtamembrane (C/F)CV+C motifs mediate pMHCII-specific signaling independently of CD4-LCK interactions

    Mark S Lee, Peter J Tuohy ... Michael S Kuhns
    Research Advance

    CD4+ T cell activation is driven by five-module receptor complexes. The T cell receptor (TCR) is the receptor module that binds composite surfaces of peptide antigens embedded within MHCII molecules (pMHCII). It associates with three signaling modules (CD3γε, CD3δε, and CD3ζζ) to form TCR-CD3 complexes. CD4 is the coreceptor module. It reciprocally associates with TCR-CD3-pMHCII assemblies on the outside of a CD4+ T cells and with the Src kinase, LCK, on the inside. Previously, we reported that the CD4 transmembrane GGXXG and cytoplasmic juxtamembrane (C/F)CV+C motifs found in eutherian (placental mammal) CD4 have constituent residues that evolved under purifying selection (Lee et al., 2022). Expressing mutants of these motifs together in T cell hybridomas increased CD4-LCK association but reduced CD3ζ, ZAP70, and PLCγ1 phosphorylation levels, as well as IL-2 production, in response to agonist pMHCII. Because these mutants preferentially localized CD4-LCK pairs to non-raft membrane fractions, one explanation for our results was that they impaired proximal signaling by sequestering LCK away from TCR-CD3. An alternative hypothesis is that the mutations directly impacted signaling because the motifs normally play an LCK-independent role in signaling. The goal of this study was to discriminate between these possibilities. Using T cell hybridomas, our results indicate that: intracellular CD4-LCK interactions are not necessary for pMHCII-specific signal initiation; the GGXXG and (C/F)CV+C motifs are key determinants of CD4-mediated pMHCII-specific signal amplification; the GGXXG and (C/F)CV+C motifs exert their functions independently of direct CD4-LCK association. These data provide a mechanistic explanation for why residues within these motifs are under purifying selection in jawed vertebrates. The results are also important to consider for biomimetic engineering of synthetic receptors.

    1. Genetics and Genomics
    2. Immunology and Inflammation

    Transposable elements regulate thymus development and function

    Jean-David Larouche, Céline M Laumont ... Claude Perreault
    Research Article

    Transposable elements (TEs) are repetitive sequences representing ~45% of the human and mouse genomes and are highly expressed by medullary thymic epithelial cells (mTECs). In this study, we investigated the role of TEs on T-cell development in the thymus. We performed multiomic analyses of TEs in human and mouse thymic cells to elucidate their role in T-cell development. We report that TE expression in the human thymus is high and shows extensive age- and cell lineage-related variations. TE expression correlates with multiple transcription factors in all cell types of the human thymus. Two cell types express particularly broad TE repertoires: mTECs and plasmacytoid dendritic cells (pDCs). In mTECs, transcriptomic data suggest that TEs interact with transcription factors essential for mTEC development and function (e.g., PAX1 and REL), and immunopeptidomic data showed that TEs generate MHC-I-associated peptides implicated in thymocyte education. Notably, AIRE, FEZF2, and CHD4 regulate small yet non-redundant sets of TEs in murine mTECs. Human thymic pDCs homogenously express large numbers of TEs that likely form dsRNA, which can activate innate immune receptors, potentially explaining why thymic pDCs constitutively secrete IFN ɑ/β. This study highlights the diversity of interactions between TEs and the adaptive immune system. TEs are genetic parasites, and the two thymic cell types most affected by TEs (mTEcs and pDCs) are essential to establishing central T-cell tolerance. Therefore, we propose that orchestrating TE expression in thymic cells is critical to prevent autoimmunity in vertebrates.