Recruitment of CXCR4+ type 1 innate lymphoid cells distinguishes sarcoidosis from other skin granulomatous diseases.

Sarcoidosis is a multiorgan granulomatous disease that lacks diagnostic biomarkers and targeted treatments. Using blood and skin from patients with sarcoid and non-sarcoid skin granulomas, we discovered that skin granulomas from different diseases exhibit unique immune cell recruitment and molecular signatures. Sarcoid skin granulomas were specifically enriched for type 1 innate lymphoid cells (ILC1s) and B cells and exhibited molecular programs associated with formation of mature tertiary lymphoid structures (TLSs), including increased CXCL12/CXCR4 signaling. Lung sarcoidosis granulomas also displayed similar immune cell recruitment. Thus, granuloma formation was not a generic molecular response. In addition to tissue-specific effects, patients with sarcoidosis exhibited an 8-fold increase in circulating ILC1s, which correlated with treatment status. Multiple immune cell types induced CXCL12/CXCR4 signaling in sarcoidosis, including Th1 T cells, macrophages, and ILCs. Mechanistically, CXCR4 inhibition reduced sarcoidosis-activated immune cell migration, and targeting CXCR4 or total ILCs attenuated granuloma formation in a noninfectious mouse model. Taken together, our results show that ILC1s are a tissue and circulating biomarker that distinguishes sarcoidosis from other skin granulomatous diseases. Repurposing existing CXCR4 inhibitors may offer a new targeted treatment for this devastating disease.

UBE2N Is Essential for Maintenance of Skin Homeostasis and Suppression of Inflammation.

UBE2N, a Lys63 ubiquitin-conjugating enzyme, plays critical roles in embryogenesis and immune system development and function. However, its roles in adult epithelial tissue homeostasis and pathogenesis are unclear. We generated conditional mouse models that deleted Ube2n in skin cells in a temporally and spatially controlled manner. We found that Ube2n knockout in the adult skin keratinocytes induced a range of inflammatory skin defects characteristic of psoriatic and actinic keratosis. These included inflammation, epidermal and dermal thickening, parakeratosis, and increased immune cell infiltration as well as signs of edema and blistering. Single-cell transcriptomic analyses and RT-qPCR showed that Ube2n-knockout keratinocytes expressed elevated myeloid cell chemoattractants such as Cxcl1 and Cxcl2 and decreased the homeostatic T lymphocyte chemoattractant Ccl27a. Consistently, the infiltrating immune cells were predominantly myeloid-derived cells, including neutrophils and M1-like macrophages, which expressed high levels of inflammatory cytokines such as Il1ß and Il24. Pharmacological blockade of the IL-1 receptor associated kinases (IRAK1/4) alleviated inflammation, epidermal and dermal thickening, and immune infiltration of the Ube2n-mutant skin. Together, these findings highlight a key role of keratinocyte UBE2N in maintenance of epidermal homeostasis and skin immunity and identify IRAK1/4 as potential therapeutic target for inflammatory skin disorders.

A neural tract-inspired conduit for facile, on-demand biopsy of glioblastoma.

Background: A major hurdle to effectively treating glioblastoma (GBM) patients is the lack of longitudinal information about tumor progression, evolution, and treatment response. Methods: In this study, we report the use of a neural tract-inspired conduit containing aligned polymeric nanofibers (i.e., an aligned nanofiber device) to enable on-demand access to GBM tumors in 2 rodent models. Depending on the experiment, a humanized U87MG xenograft and/or F98-GFP+ syngeneic rat tumor model was chosen to test the safety and functionality of the device in providing continuous sampling access to the tumor and its microenvironment. Results: The aligned nanofiber device was safe and provided a high quantity of quality genomic materials suitable for omics analyses and yielded a sufficient number of live cells for in vitro expansion and screening. Transcriptomic and genomic analyses demonstrated continuity between material extracted from the device and that of the primary, intracortical tumor (in the in vivo model). Conclusions: The results establish the potential of this neural tract-inspired, aligned nanofiber device as an on-demand, safe, and minimally invasive access point, thus enabling rapid, high-throughput, longitudinal assessment of tumor and its microenvironment, ultimately leading to more informed clinical treatment strategies.

Spatial transcriptomics reveals segregation of tumor cell states in glioblastoma and marked immunosuppression within the perinecrotic niche.

Glioblastoma (GBM) remains an untreatable malignant tumor with poor patient outcomes, characterized by palisading necrosis and microvascular proliferation. While single-cell technology made it possible to characterize different lineage of glioma cells into neural progenitor-like (NPC-like), oligodendrocyte-progenitor-like (OPC-like), astrocyte-like (AC-like) and mesenchymal like (MES-like) states, it does not capture the spatial localization of these tumor cell states. Spatial transcriptomics empowers the study of the spatial organization of different cell types and tumor cell states and allows for the selection of regions of interest to investigate region-specific and cell-type-specific pathways. Here, we obtained paired 10x Chromium single-nuclei RNA-sequencing (snRNA-seq) and 10x Visium spatial transcriptomics data from three GBM patients to interrogate the GBM microenvironment. Integration of the snRNA-seq and spatial transcriptomics data reveals patterns of segregation of tumor cell states. For instance, OPC-like tumor and NPC-like tumor significantly segregate in two of the three samples. Our differentially expressed gene and pathway analyses uncovered significant pathways in functionally relevant niches. Specifically, perinecrotic regions were more immunosuppressive than the endogenous GBM microenvironment, and perivascular regions were more pro-inflammatory. Our gradient analysis suggests that OPC-like tumor cells tend to reside in areas closer to the tumor vasculature compared to tumor necrosis, which may reflect increased oxygen requirements for OPC-like cells. In summary, we characterized the localization of cell types and tumor cell states, the gene expression patterns, and pathways in different niches within the GBM microenvironment. Our results provide further evidence of the segregation of tumor cell states and highlight the immunosuppressive nature of the necrotic and perinecrotic niches in GBM.

UBE2N is essential for maintenance of skin homeostasis and suppression of inflammation.

UBE2N, a Lys63-ubiquitin conjugating enzyme, plays critical roles in embryogenesis and immune system development and function. However, its roles in adult epithelial tissue homeostasis and pathogenesis are unclear. We generated conditional mouse models that deleted Ube2n in skin cells in a temporally and spatially controlled manner. We found that Ube2n-knockout (KO) in the adult skin keratinocytes induced a range of inflammatory skin defects characteristic of psoriatic and actinic keratosis. These included eczematous inflammation, epidermal and dermal thickening, parakeratosis, and increased immune cell infiltration, as well as signs of edema and blistering. Single cell transcriptomic analyses and RT-qPCR showed that Ube2n KO keratinocytes expressed elevated myeloid cell chemo-attractants such as Cxcl1 and Cxcl2 and decreased the homeostatic T lymphocyte chemo-attractant, Ccl27a. Consistently, the infiltrating immune cells of Ube2n-KO skin were predominantly myeloid-derived cells including neutrophils and M1-like macrophages that were highly inflammatory, as indicated by expression of Il1ß and Il24. Pharmacological blockade of the IL-1 receptor associated kinases (IRAK1/4) alleviated eczema, epidermal and dermal thickening, and immune infiltration of the Ube2n mutant skin. Together, these findings highlight a key role of keratinocyte-UBE2N in maintenance of epidermal homeostasis and skin immunity and identify IRAK1/4 as potential therapeutic target for inflammatory skin disorders.

Ataxia-telangiectasia mutated ( Atm ) disruption sensitizes spatially-directed H3.3K27M/TP53 diffuse midline gliomas to radiation therapy.

Diffuse midline gliomas (DMGs) are lethal brain tumors characterized by p53-inactivating mutations and oncohistone H3.3K27M mutations that rewire the cellular response to genotoxic stress, which presents therapeutic opportunities. We used RCAS/tv-a retroviruses and Cre recombinase to inactivate p53 and induce K27M in the native H3f3a allele in a lineage- and spatially-directed manner, yielding primary mouse DMGs. Genetic or pharmacologic disruption of the DNA damage response kinase Ataxia-telangiectasia mutated (ATM) enhanced the efficacy of focal brain irradiation, extending mouse survival. This finding suggests that targeting ATM will enhance the efficacy of radiation therapy for p53-mutant DMG but not p53-wildtype DMG. We used spatial in situ transcriptomics and an allelic series of primary murine DMG models with different p53 mutations to identify transactivation-independent p53 activity as a key mediator of such radiosensitivity. These studies deeply profile a genetically faithful and versatile model of a lethal brain tumor to identify resistance mechanisms for a therapeutic strategy currently in clinical trials.

The RNA helicase DDX39B activates FOXP3 RNA splicing to control T regulatory cell fate.

Genes associated with increased susceptibility to multiple sclerosis (MS) have been identified, but their functions are incompletely understood. One of these genes codes for the RNA helicase DExD/H-Box Polypeptide 39B (DDX39B), which shows genetic and functional epistasis with interleukin-7 receptor-α gene (IL7R) in MS-risk. Based on evolutionary and functional arguments, we postulated that DDX39B enhances immune tolerance thereby decreasing MS risk. Consistent with such a role we show that DDX39B controls the expression of many MS susceptibility genes and important immune-related genes. Among these we identified Forkhead Box P3 (FOXP3), which codes for the master transcriptional factor in CD4+/CD25+ T regulatory cells. DDX39B knockdown led to loss of immune-regulatory and gain of immune-effector expression signatures. Splicing of FOXP3 introns, which belong to a previously unrecognized type of introns with C-rich polypyrimidine tracts, was exquisitely sensitive to DDX39B levels. Given the importance of FOXP3 in autoimmunity, this work cements DDX39B as an important guardian of immune tolerance.

Ganglioglioma deep transcriptomics reveals primitive neuroectoderm neural precursor-like population.

Gangliogliomas are brain tumors composed of neuron-like and macroglia-like components that occur in children and young adults. Gangliogliomas are often characterized by a rare population of immature astrocyte-appearing cells expressing CD34, a marker expressed in the neuroectoderm (neural precursor cells) during embryogenesis. New insights are needed to refine tumor classification and to identify therapeutic approaches. We evaluated five gangliogliomas with single nucleus RNA-seq, cellular indexing of transcriptomes and epitopes by sequencing, and/or spatially-resolved RNA-seq. We uncovered a population of CD34+ neoplastic cells with mixed neuroectodermal, immature astrocyte, and neuronal markers. Gene regulatory network interrogation in these neuroectoderm-like cells revealed control of transcriptional programming by TCF7L2/MEIS1-PAX6 and SOX2, similar to that found during neuroectodermal/neural development. Developmental trajectory analyses place neuroectoderm-like tumor cells as precursor cells that give rise to neuron-like and macroglia-like neoplastic cells. Spatially-resolved transcriptomics revealed a neuroectoderm-like tumor cell niche with relative lack of vascular and immune cells. We used these high resolution results to deconvolute clinically-annotated transcriptomic data, confirming that CD34+ cell-associated gene programs associate with gangliogliomas compared to other glial brain tumors. Together, these deep transcriptomic approaches characterized a ganglioglioma cellular hierarchy-confirming CD34+ neuroectoderm-like tumor precursor cells, controlling transcription programs, cell signaling, and associated immune cell states. These findings may guide tumor classification, diagnosis, prognostication, and therapeutic investigations.

Aging and obesity prime the methylome and transcriptome of adipose stem cells for disease and dysfunction.

The epigenome of stem cells occupies a critical interface between genes and environment, serving to regulate expression through modification by intrinsic and extrinsic factors. We hypothesized that aging and obesity, which represent major risk factors for a variety of diseases, synergistically modify the epigenome of adult adipose stem cells (ASCs). Using integrated RNA- and targeted bisulfite-sequencing in murine ASCs from lean and obese mice at 5- and 12-months of age, we identified global DNA hypomethylation with either aging or obesity, and a synergistic effect of aging combined with obesity. The transcriptome of ASCs in lean mice was relatively stable to the effects of age, but this was not true in obese mice. Functional pathway analyses identified a subset of genes with critical roles in progenitors and in diseases of obesity and aging. Specifically, Mapt, Nr3c2, App, and Ctnnb1 emerged as potential hypomethylated upstream regulators in both aging and obesity (AL vs. YL and AO vs. YO), and App, Ctnnb1, Hipk2, Id2, and Tp53 exhibited additional effects of aging in obese animals. Furthermore, Foxo3 and Ccnd1 were potential hypermethylated upstream regulators of healthy aging (AL vs. YL), and of the effects of obesity in young animals (YO vs. YL), suggesting that these factors could play a role in accelerated aging with obesity. Finally, we identified candidate driver genes that appeared recurrently in all analyses and comparisons undertaken. Further mechanistic studies are needed to validate the roles of these genes capable of priming ASCs for dysfunction in aging- and obesity-associated pathologies.

Adolescent peer struggles predict accelerated epigenetic aging in midlife.

This study examined struggles to establish autonomy and relatedness with peers in adolescence and early adulthood as predictors of advanced epigenetic aging assessed at age 30. Participants (N = 154; 67 male and 87 female) were observed repeatedly, along with close friends and romantic partners, from ages 13 through 29. Observed difficulty establishing close friendships characterized by mutual autonomy and relatedness from ages 13 to 18, an interview-assessed attachment state of mind lacking autonomy and valuing of attachment at 24, and self-reported difficulties in social integration across adolescence and adulthood were all linked to greater epigenetic age at 30, after accounting for chronological age, gender, race, and income. Analyses assessing the unique and combined effects of these factors, along with lifetime history of cigarette smoking, indicated that each of these factors, except for adult social integration, contributed uniquely to explaining epigenetic age acceleration. Results are interpreted as evidence that the adolescent preoccupation with peer relationships may be highly functional given the relevance of such relationships to long-term physical outcomes.

Resistance of mitochondrial DNA to cadmium and Aflatoxin B1 damage-induced germline mutation accumulation in C. elegans.

Mitochondrial DNA (mtDNA) is prone to mutation in aging and over evolutionary time, yet the processes that regulate the accumulation of de novo mtDNA mutations and modulate mtDNA heteroplasmy are not fully elucidated. Mitochondria lack certain DNA repair processes, which could contribute to polymerase error-induced mutations and increase susceptibility to chemical-induced mtDNA mutagenesis. We conducted error-corrected, ultra-sensitive Duplex Sequencing to investigate the effects of two known nuclear genome mutagens, cadmium and Aflatoxin B1, on germline mtDNA mutagenesis in Caenorhabditis elegans. Detection of thousands of mtDNA mutations revealed pervasive heteroplasmy in C. elegans and that mtDNA mutagenesis is dominated by C:G → A:T mutations generally attributed to oxidative damage. However, there was no effect of either exposure on mtDNA mutation frequency, spectrum, or trinucleotide context signature despite a significant increase in nuclear mutation rate after aflatoxin B1 exposure. Mitophagy-deficient mutants pink-1 and dct-1 accumulated significantly higher levels of mtDNA damage compared to wild-type C. elegans after exposures. However, there were only small differences in mtDNA mutation frequency, spectrum, or trinucleotide context signature compared to wild-type after 3050 generations, across all treatments. These findings suggest mitochondria harbor additional previously uncharacterized mechanisms that regulate mtDNA mutational processes across generations.

Lifetime marijuana use and epigenetic age acceleration: A 17-year prospective examination.

AIMS: This study was designed to assess links between lifetime levels of marijuana use and accelerated epigenetic aging. DESIGN: Prospective longitudinal study, following participants annually from age 13 to age 30. SETTING AND PARTICIPANTS: A community sample of 154 participants recruited from a small city in the Southeastern United States. MEASUREMENTS: Participants completed annual assessments of marijuana use from age 13 to age 29 and provided blood samples that yielded two indices of epigenetic aging (DNAmGrimAge and DunedinPoAm) at age 30. Additional covariates examined included history of cigarette smoking, anxiety and depressive symptoms, childhood illness, gender, adolescent-era family income, and racial/ethnic minority status. FINDINGS: Lifetime marijuana use predicted accelerated epigenetic aging, with effects remaining even after covarying cell counts, demographic factors and chronological age (ß's = 0.32 & 0.27, p's < 0.001, 95% CI's = 0.21-0.43 & 0.16-0.39 for DNAmGrimAge and DunedinPoAm, respectively). Predictions remained after accounting for cigarette smoking (ß's = 0.25 & 0.21, respectively, p's < 0.001, 95% CI's = 0.14-0.37 & 0.09-0.32 for DNAmGrimAge and DunedinPoAm, respectively). A dose-response effect was observed and there was also evidence that effects were dependent upon recency of use. Effects of marijuana use appeared to be fully mediated by hypomethylation of a site linked to effects of hydrocarbon inhalation (cg05575921). CONCLUSIONS: Marijuana use predicted epigenetic changes linked to accelerated aging, with evidence suggesting that effects may be primarily due to hydrocarbon inhalation among marijuana smokers. Further research is warranted to explore mechanisms underlying this linkage.

Characterization of a castrate-resistant prostate cancer xenograft derived from a patient of West African ancestry.

BACKGROUND: Prostate cancer is a clinically and molecularly heterogeneous disease, with highest incidence and mortality among men of African ancestry. To date, prostate cancer patient-derived xenograft (PCPDX) models to study this disease have been difficult to establish because of limited specimen availability and poor uptake rates in immunodeficient mice. Ancestrally diverse PCPDXs are even more rare, and only six PCPDXs from self-identified African American patients from one institution were recently made available. METHODS: In the present study, we established a PCPDX from prostate cancer tissue from a patient of estimated 90% West African ancestry with metastatic castration resistant disease, and characterized this model's pathology, karyotype, hotspot mutations, copy number, gene fusions, gene expression, growth rate in normal and castrated mice, therapeutic response, and experimental metastasis. RESULTS: This PCPDX has a mutation in TP53 and loss of PTEN and RB1. We have documented a 100% take rate in mice after thawing the PCPDX tumor from frozen stock. The PCPDX is castrate- and docetaxel-resistant and cisplatin-sensitive, and has gene expression patterns associated with such drug responses. After tail vein injection, the PCPDX tumor cells can colonize the lungs of mice. CONCLUSION: This PCPDX, along with others that are established and characterized, will be useful pre-clinically for studying the heterogeneity of prostate cancer biology and testing new therapeutics in models expected to be reflective of the clinical setting.

Circulating Tumor Cell Chromosomal Instability and Neuroendocrine Phenotype by Immunomorphology and Poor Outcomes in Men with mCRPC Treated with Abiraterone or Enzalutamide.

PURPOSE: While the detection of AR-V7 in circulating tumor cells (CTC) is associated with resistance to abiraterone or enzalutamide in men with metastatic castration-resistant prostate cancer (mCRPC), it only accounts for a minority of this resistance. Neuroendocrine (NE) differentiation or chromosomal instability (CIN) may be additional mechanisms that mediate resistance. EXPERIMENTAL DESIGN: PROPHECY was a multicenter prospective study of men with high-risk mCRPC starting abiraterone or enzalutamide. A secondary objective was to assess Epic CTC CIN and NE phenotypes before abiraterone or enzalutamide and at progression. The proportional hazards (PH) model was used to investigate the prognostic importance of CIN and NE in predicting progression-free survival and overall survival (OS) adjusting for CTC number (CellSearch), AR-V7, prior therapy, and clinical risk score. The PH model was utilized to validate this association of NE with OS in an external dataset of patients treated similarly at Memorial Sloan Kettering Cancer Center (MSKCC; New York, NY). RESULTS: We enrolled 118 men with mCRPC starting on abiraterone or enzalutamide; 107 were evaluable on the Epic platform. Of these, 36.4% and 8.4% were CIN positive and NE positive, respectively. CIN and NE were independently associated with worse OS [HR, 2.2; 95% confidence interval (CI), 1.2-4.0 and HR 3.8; 95% CI, 1.2-12.3, respectively] when treated with abiraterone/enzalutamide. The prognostic significance of NE positivity for worse OS was confirmed in the MSKCC dataset (n = 173; HR, 5.7; 95% CI, 2.6-12.7). CONCLUSIONS: A high CIN and NE CTC phenotype is independently associated with worse survival in men with mCRPC treated with abiraterone/enzalutamide, warranting further prospective controlled predictive studies to inform treatment decisions.

Circulating Tumor Cell Genomic Evolution and Hormone Therapy Outcomes in Men with Metastatic Castration-Resistant Prostate Cancer.

Men with circulating tumor cell (CTC) AR-V7-positive metastatic castration-resistant prostate cancer (mCRPC) have worse outcomes when treated with enzalutamide/abiraterone. However, most men lack CTC AR-V7 detection, and additional predictive biomarkers are needed. We conducted a retrospective secondary analysis of the prospective PROPHECY trial (NCT02269982) of men with mCRPC undergoing treatment with enzalutamide/abiraterone, analyzing pooled CTC and germline DNA for whole-genome copy-number alterations (CNA) in 73 samples from 48 men over time along with pooled CTC and germline whole-exome sequencing on 22 paired samples before and following progression on androgen receptor (AR) inhibitor therapy to identify somatic genomic alterations associated with acquired resistance. We observed broad interpatient and longitudinal CTC genomic heterogeneity from AR-V7-negative men with mCRPC, including common gains of KDM6A, MYCN, and AR, and loss of ZFHX3, BRCA1, and PTEN. Men who had progression-free survival of ≤3 months despite enzalutamide/abiraterone treatment were more likely to have baseline CTC genomic loss of CHD1, PTEN, PHLPP1, and ZFHX3 and gains of BRCA2, KDM5D, MYCN, and SPARC. After progression on abiraterone/enzalutamide, we observed clonal evolution of CTCs harboring TP53 mutations and gain of ATM, KDM6A, and MYC, and loss of NCOR1, PTEN, RB1, and RUNX2. CTC genomic findings were independently confirmed in a separate cohort of mCRPC men who progressed despite prior treatment with abiraterone/enzalutamide (NCT02204943). IMPLICATIONS: We identified common and reproducible genomic alterations in CTCs from AR-V7-negative mCRPC men associated with poor outcomes during enzalutamide/abiraterone treatment, including CNAs in genes linked to lineage plasticity and epigenetic signaling, DNA repair, AR, TP53/RB1, PTEN, and WNT pathways.

A non-canonical type 2 immune response coordinates tuberculous granuloma formation and epithelialization.

The central pathogen-immune interface in tuberculosis is the granuloma, a complex host immune structure that dictates infection trajectory and physiology. Granuloma macrophages undergo a dramatic transition in which entire epithelial modules are induced and define granuloma architecture. In tuberculosis, relatively little is known about the host signals that trigger this transition. Using the zebrafish-Mycobacterium marinum model, we identify the basis of granuloma macrophage transformation. Single-cell RNA-sequencing analysis of zebrafish granulomas and analysis of Mycobacterium tuberculosis-infected macaques reveal that, even in the presence of robust type 1 immune responses, countervailing type 2 signals associate with macrophage epithelialization. We find that type 2 immune signaling, mediated via stat6, is absolutely required for epithelialization and granuloma formation. In mixed chimeras, stat6 acts cell autonomously within macrophages, where it is required for epithelioid transformation and incorporation into necrotic granulomas. These findings establish the signaling pathway that produces the hallmark structure of mycobacterial infection.

U2AF2 binds IL7R exon 6 ectopically and represses its inclusion.

Interleukin 7 receptor α-chain is crucial for the development and maintenance of T cells and is genetically associated with autoimmune disorders including multiple sclerosis (MS), a demyelinating disease of the CNS. Exon 6 of IL7R encodes for the transmembrane domain of the receptor and is regulated by alternative splicing: inclusion or skipping of IL7R exon 6 results in membrane-bound or soluble IL7R isoforms, respectively. We previously identified a SNP (rs6897932) in IL7R exon 6, strongly associated with MS risk and showed that the risk allele (C) increases skipping of the exon, resulting in elevated levels of sIL7R. This has important pathological consequences as elevated levels of sIL7R has been shown to exacerbate the disease in the experimental autoimmune encephalomyelitis mouse model of MS. Understanding the regulation of exon 6 splicing provides important mechanistic insights into the pathogenesis of MS. Here we report two mechanisms by which IL7R exon 6 is controlled. First, a competition between PTBP1 and U2AF2 at the polypyrimidine tract (PPT) of intron 5, and second, an unexpected U2AF2-mediated assembly of spicing factors in the exon. We noted the presence of a branchpoint sequence (BPS) (TACTAAT or TACTAAC) within exon 6, which is stronger with the C allele. We also noted that the BPS is followed by a PPT and conjectured that silencing could be mediated by the binding of U2AF2 to that tract. In support of this model, we show that evolutionary conservation of the exonic PPT correlates well with the degree of alternative splicing of exon 6 in two non-human primate species and that U2AF2 binding to this PPT recruits U2 snRNP components to the exon. These observations provide the first explanation for the stronger silencing of IL7R exon 6 with the disease associated C allele at rs6897932.

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

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.