Cellular Transcription Profile Research Articles

Single-nucleus RNA sequencing reveals midgut cellular heterogeneity and transcriptional profiles in Bombyx mori cytoplasmic polyhedrosis virus infection.

The gut is not only used by insects as an organ for the digestion of food and absorption of nutrients but also as an important barrier against the invasion and proliferation of pathogenic microorganisms. Bombyx mori cytoplasmic polyhedrosis virus (BmCPV), an insect-specific virus, predominantly colonizes the midgut epithelial cells of the silkworm, thereby jeopardizing its normal growth. However, there is limited knowledge of the cellular immune responses to viral infection and whether the infection is promoted or inhibited by different types of cells in the silkworm midgut. In this study, we used single-nucleus RNA sequencing to identify representative enteroendocrine cells, enterocytes, and muscle cell types in the silkworm midgut. In addition, by analyzing the transcriptional profiles of various subpopulations in the infected and uninfected groups, we found that BmCPV infection suppresses the response of the antiviral pathways and induces the expression of BmHSP70, which plays a role in promoting BmCPV replication. However, certain immune genes in the midgut of the silkworm, such as BmLebocin3, were induced upon viral infection, and downregulation of BmLEB3 using RNA interference promoted BmCPV replication in the midgut of B. mori. These results suggest that viral immune evasion and active host resistance coexist in BmCPV-infected silkworms. We reveal the richness of cellular diversity in the midgut of B. mori larvae by single-nucleus RNA sequencing analysis and provide new insights into the complex interactions between the host and the virus at the single-cell level.

Single cell RNA sequencing provides novel cellular transcriptional profiles and underlying pathogenesis of presbycusis

Age-related hearing loss (ARHL) or presbycusis is associated with irreversible progressive damage in the inner ear, where the sound is transduced into electrical signal; but the detailed mechanism remains unclear. Here, we sought to determine the potential molecular mechanism involved in the pathogeneses of ARHL with bioinformatics methods. A single-cell transcriptome sequencing study was performed on the cochlear samples from young and aged mice. Detection of identified cell type marker allowed us to screen 18 transcriptional clusters, including myeloid cells, epithelial cells, B cells, endothelial cells, fibroblasts, T cells, inner pillar cells, neurons, inner phalangeal cells, and red blood cells. Cell-cell communications were analyzed between young and aged cochlear tissue samples by using the latest integration algorithms Cellchat. A total of 56 differentially expressed genes were screened between the two groups. Functional enrichment analysis showed these genes were mainly involved in immune, oxidative stress, apoptosis, and metabolic processes. The expression levels of crucial genes in cochlear tissues were further verified by immunohistochemistry. Overall, this study provides new theoretical support for the development of clinical therapeutic drugs.

Mapping herpesvirus-driven impacts on the cellular milieu and transcriptional profile of Kaposi sarcoma in patient-derived mouse models.

Kaposi sarcoma (KS) is defined by aberrant angiogenesis driven by Kaposi sarcoma herpesvirus (KSHV)-infected spindle cells with endothelial characteristics. KS research is hindered by rapid loss of KSHV infection upon explant culture of tumor cells. Here, we establish patient-derived KS xenografts (PDXs) upon orthotopic implantation of cutaneous KS biopsies in immunodeficient mice. KS tumors were maintained in 27/28 PDX until experimental endpoint, up to 272 days in the first passage of recipient mice. KSHV latency associated nuclear antigen (LANA)+ endothelial cell density increased by a mean 4.3-fold in 14/15 PDX analyzed by IHC at passage 1 compared to respective input biopsies, regardless of implantation variables and clinical features of patients. The Ki-67 proliferation marker colocalized with LANA more frequently in PDXs. Spatial transcriptome analysis revealed increased expression of viral transcripts from latent and lytic gene classes in the PDX. The expanded KSHV+ regions of the PDX maintained signature gene expression of KS tumors, with enrichment in pathways associated with angiogenesis and endothelium development. Cells with characteristics of tumor-associated fibroblasts derived from PDX were propagated for 15 passages. These fibroblast-like cells were permissive for de novo KSHV infection, and one lineage produced CXCL12, a cancer-promoting chemokine. Spatial analysis revealed that fibroblasts are a likely source of CXCL12 signaling to CXCR4 that was upregulated in KS regions. The reproducible expansion of KSHV-infected endothelial cells in PDX from multiple donors and recapitulation of a KS tumor gene signature supports the application of patient-derived KS mouse models for studies of pathogenesis and novel therapies.

Meningeal lymphatic supporting cells govern the formation and maintenance of zebrafish mural lymphatic endothelial cells

The meninges are critical for the brain functions, but the diversity of meningeal cell types and intercellular interactions have yet to be thoroughly examined. Here we identify a population of meningeal lymphatic supporting cells (mLSCs) in the zebrafish leptomeninges, which are specifically labeled by ependymin. Morphologically, mLSCs form membranous structures that enwrap the majority of leptomeningeal blood vessels and all the mural lymphatic endothelial cells (muLECs). Based on its unique cellular morphologies and transcriptional profile, mLSC is characterized as a unique cell type different from all the currently known meningeal cell types. Because of the formation of supportive structures and production of pro-lymphangiogenic factors, mLSCs not only promote muLEC development and maintain the dispersed distributions of muLECs in the leptomeninges, but also are required for muLEC regeneration after ablation. This study characterizes a newly identified cell type in leptomeninges, mLSC, which is required for muLEC development, maintenance, and regeneration.

Single-cell atlas of human infrapatellar fat pad and synovium implicates APOE signaling in osteoarthritis pathology.

The infrapatellar fat pad (IPFP) and synovium play essential roles in maintaining knee joint homeostasis and in the progression of osteoarthritis (OA). The cellular and transcriptional mechanisms regulating the function of these specialized tissues under healthy and diseased conditions are largely unknown. Here, single-cell and single-nuclei RNA sequencing of human IPFP and synovial tissues were performed to elucidate the cellular composition and transcriptional profile. Computational trajectory analysis revealed that dipeptidyl peptidase 4+ mesenchymal cells function as a common progenitor for IPFP adipocytes and synovial lining layer fibroblasts, suggesting that IPFP and synovium represent an integrated tissue unit. OA induced a profibrotic and inflammatory phenotype in mesenchymal lineage cells with biglycan+ intermediate fibroblasts as a major contributor to OA fibrosis. Apolipoprotein E (APOE) signaling from intermediate fibroblasts and macrophages was identified as a critical regulatory factor. Ex vivo incubation of human cartilage with soluble APOE accelerated proteoglycan degeneration. Inhibition of APOE signaling by intra-articular injection of an anti-APOE neutralizing antibody attenuated the progression of collagenase-induced OA in mice, demonstrating a detrimental effect of APOE on cartilage. Our studies provide a framework for designing further therapeutic strategies for OA by describing the cellular and transcriptional landscape of human IPFP and synovium in healthy versus OA joints.

SAT355 Single-Cell RNA Sequencing Pipeline for the Human Ovary Using IVF Tissue

Abstract Disclosure: A.J. Hanson: None. S. Shah: None. C. Wai: None. R. Bauer: None. C. Boots: None. M. Urbanek: None. Objective: Polycystic Ovary Syndrome (PCOS) is the most common endocrinopathy of reproductive aged women, affecting 6-15% of women in this age group. Defined by hyperandrogenism and irregular menses, PCOS is a complex trait with a strong genetic component; its underlying cause remains poorly understood. To understand the molecular pathogenesis of PCOS, an appropriate experimental model of the human ovary is needed. The ovary is the primary source of androgen production and an essential regulator of ovulation and menstruation, the key processes disrupted in PCOS. Discard tissue from oocyte retrievals as part of assisted reproduction treatment (ART) is a plentiful source of human ovarian tissue, requires no additional risk to the patient, and is enriched for subjects with PCOS. Therefore, we developed a pipeline using discard ovarian aspirate tissue from reproductively healthy women and women with PCOS undergoing oocyte retrievals for ART to generate a single-cell (sc) transcriptomic atlas of human ovaries. We will use those maps to define the differences in gene expression between the two groups and determine the impact of genetic variation on human ovary function in women with PCOS. Materials and Methods: Our pipeline consists of i. harvesting oocyte retrieval discard tissue including follicular fluid, pieces of ovary, and blood from reproductively healthy women and women with PCOS ii. generating a sc suspension of ovarian cells, iii. constructing and sequencing a scRNA sequencing library, iv. building a transcriptomic atlas of the ovary, and v. identifying differences in cell abundance and expression between reproductively healthy women and women with PCOS. Results: We have processed samples from seven reproductively healthy women and eight women with PCOS using this scRNA sequencing pipeline. The sc isolation protocol resulted in a sc suspension with high viability (95%) and yield (19.6 x106 cells/mL), demonstrating that we can generate a large quantity of high-quality sc suspensions. The sequenced single cells were clustered into distinct groups based on their transcriptomes including cell types predicted to be affected by PCOS including theca, granulosa, and immune cells. Conclusions: We successfully created a scRNA sequencing pipeline from discard oocyte retrieval tissue and generated transcriptomic maps of the human ovary of reproductively healthy women and women with PCOS, demonstrating the feasibility of our approach. Impact Statement: By integrating our transcriptomic data with whole exome sequencing (WES) data to assess the impact of genomic variation on ovarian cellular composition and transcriptional profiles, we will, for the first time, evaluate the impact of an individual woman’s genetic variation on their ovarian functional profile, a critical first step towards understanding the underlying pathology of PCOS and personalized medicine. Presentation Date: Saturday, June 17, 2023

Comparison of cell type distribution between single-cell and single-nucleus RNA sequencing: enrichment of adherent cell types in single-nucleus RNA sequencing

Single-cell ribonucleic acid (RNA) sequencing (scRNA-seq) is an effective technique for estimating the cellular composition and transcriptional profiles of individual cells from fresh tissue. Single-nucleus RNA sequencing (snRNA-seq) is necessary to perform this type of analysis in frozen or difficult-to-dissociate tissues, which cannot be subjected to scRNA-seq. This difference in the state of tissues leads to variation in cell-type distributions among each platform. To identify the characteristics of these methods and their differences, scRNA-seq and snRNA-seq were performed in parallel for colon and liver tissues. The two platforms revealed similar diversity but different proportions of cell types in matched tissues. The proportions of epithelial cells in the colon and hepatocytes in the liver were relatively high in snRNA-seq and that of immune cells was relatively high in scRNA-seq. This difference could be explained by variations in the expression scores of adhesion genes due to the disruption of the cytoplasmic contents during scRNA-seq. The enrichment of epithelial cells in the colon resulted in a discrepancy in the differentiation of epithelial cells. This enrichment was also well matched with the images of hematoxylin and eosin staining and the estimated distribution of cell types in bulk RNA sequencing. These results showed that snRNA-seq could be used to analyze tissues that cannot be subjected to scRNA-seq and provides more information in specific cell type analysis.

Niche deconvolution of the glioblastoma proteome reveals a distinct infiltrative phenotype within the proneural transcriptomic subgroup

Glioblastoma is often subdivided into three transcriptional subtypes (classical, proneural, mesenchymal) based on bulk RNA signatures that correlate with distinct genetic and clinical features. Potential cellular-level differences of these subgroups, such as the relative proportions of glioblastoma’s hallmark histopathologic features (e.g. brain infiltration, microvascular proliferation), may provide insight into their distinct phenotypes but are, however, not well understood. Here we leverage machine learning and reference proteomic profiles derived from micro-dissected samples of these major histomorphologic glioblastoma features to deconvolute and estimate niche proportions in an independent proteogenomically-characterized cohort. This approach revealed a strong association of the proneural transcriptional subtype with a diffusely infiltrating phenotype. Similarly, enrichment of a microvascular proliferation proteomic signature was seen within the mesenchymal subtype. This study is the first to link differences in the cellular pathology signatures and transcriptional profiles of glioblastoma, providing potential new insights into the genetic drivers and poor treatment response of specific subsets of glioblastomas.

Inhibition of BCL9 Modulates the Cellular Landscape of Tumor-Associated Macrophages in the Tumor Immune Microenvironment of Colorectal Cancer

Tumor-associated macrophages (TAMs) are an indispensable part of the tumor microenvironment (TME), and they likely play a negative rather than positive role in cancer treatment. However, the cellular landscape and transcriptional profile regulation of TAMs in the case of tumor gene inactivation or chemical interference remains unclear. The B-cell lymphoma 9/B-cell lymphoma 9-like (BCL9/BCL9L) is a critical transcription co-factor of β-catenin. Suppression of Bcl9 inhibits tumor growth in mouse models of colorectal cancer (CRC). Here, we studied the TAMs of CRC by single-cell sequencing. Bcl9 depletion caused macrophage polarization inhibition from M0 to M2 and changed the CRC TME, which further interferes with the inflammation of M0 and M1. The transcription factor regulating these processes may be related to the Wnt signaling pathway from multiple levels. Furthermore, we also found that the cells delineated from monocyte to NK-like non-functioning cells were significantly different in the BCL9-deprived population. Combining these data, we proposed a TAM-to-NK score to evaluate the dynamic balance in TME of monocyte/TAM cells and NK-like non-functioning cells in The Cancer Genome Atlas (TCGA) clinical samples to verify the clinical significance. We demonstrated that the cell type balance and transcription differences of TAMs regulated by BCL9-driven Wnt signaling affected immune surveillance and inflammation of cancer, ultimately affecting patients’ prognosis. We thereby highlighted the potential of targeting Wnt signaling pathway through cancer immunotherapy.

Novel interferon-sensitive genes unveiled by correlation-driven gene selection and systems biology

Interferons (IFNs) are key cytokines involved in alerting the immune system to viral infection. After IFN stimulation, cellular transcriptional profile critically changes, leading to the expression of several IFN stimulated genes (ISGs) that exert a wide variety of antiviral activities. Despite many ISGs have been already identified, a comprehensive network of coding and non-coding genes with a central role in IFN-response still needs to be elucidated. We performed a global RNA-Seq transcriptome profile of the HCV permissive human hepatoma cell line Huh7.5 and its parental cell line Huh7, upon IFN treatment, to define a network of genes whose coordinated modulation plays a central role in IFN-response. Our study adds molecular actors, coding and non-coding genes, to the complex molecular network underlying IFN-response and shows how systems biology approaches, such as correlation networks, network’s topology and gene ontology analyses can be leveraged to this aim.

Transcription Control of Liver Development.

During liver organogenesis, cellular transcriptional profiles are constantly reshaped by the action of hepatic transcriptional regulators, including FoxA1-3, GATA4/6, HNF1α/β, HNF4α, HNF6, OC-2, C/EBPα/β, Hex, and Prox1. These factors are crucial for the activation of hepatic genes that, in the context of compact chromatin, cannot access their targets. The initial opening of highly condensed chromatin is executed by a special class of transcription factors known as pioneer factors. They bind and destabilize highly condensed chromatin and facilitate access to other “non-pioneer” factors. The association of target genes with pioneer and non-pioneer transcription factors takes place long before gene activation. In this way, the underlying gene regulatory regions are marked for future activation. The process is called “bookmarking”, which confers transcriptional competence on target genes. Developmental bookmarking is accompanied by a dynamic maturation process, which prepares the genomic loci for stable and efficient transcription. Stable hepatic expression profiles are maintained during development and adulthood by the constant availability of the main regulators. This is achieved by a self-sustaining regulatory network that is established by complex cross-regulatory interactions between the major regulators. This network gradually grows during liver development and provides an epigenetic memory mechanism for safeguarding the optimal expression of the regulators.

Immune cell residency in the nasal mucosa may partially explain respiratory disease severity across the age range

Previous studies focusing on the age disparity in COVID-19 severity have suggested that younger individuals mount a more robust innate immune response in the nasal mucosa after infection with SARS-CoV-2. However, it is unclear if this reflects increased immune activation or increased immune residence in the nasal mucosa. We hypothesized that immune residency in the nasal mucosa of healthy individuals may differ across the age range. We applied single-cell RNA-sequencing and measured the cellular composition and transcriptional profile of the nasal mucosa in 35 SARS-CoV-2 negative children and adults, ranging in age from 4 months to 65 years. We analyzed in total of ~ 30,000 immune and epithelial cells and found that age and immune cell proportion in the nasal mucosa are inversely correlated, with little evidence for structural changes in the transcriptional state of a given cell type across the age range. Orthogonal validation by epigenome sequencing indicate that it is especially cells of the innate immune system that underlie the age-association. Additionally, we characterize the predominate immune cell type in the nasal mucosa: a resident T cell like population with potent antiviral properties. These results demonstrate fundamental changes in the immune cell makeup of the uninfected nasal mucosa over the lifespan. The resource we generate here is an asset for future studies focusing on respiratory infection and immunization strategies.

BRLF1-dependent viral and cellular transcriptomes and transcriptional regulation during EBV primary infection in B lymphoma cells

The immediate-early protein BRLF1 plays important roles in lytic infection of Epstein-Barr virus (EBV), in which it activates lytic viral transcription and replication. However, knowledge of the influence of BRLF1 on cellular gene expression and transcriptional reprogramming during the early lytic cycle remains limited. In the present study, deep RNA-sequencing analysis identified all differentially expressed genes (DEGs) and alternative splicing in B lymphoma cells subjected to wild-type and BRLF1-deficient EBV primary infection. The BRLF1-dependent cellular DEGs were annotated, and major differentially enriched pathways were related to DNA replication and transcription, immune and inflammatory responses, cytokine-receptor interactions and chemokine signaling and metabolic processes. Furthermore, analysis of BRLF1-binding proteins by mass spectrometry shows that BRLF1 binds to and cooperates with several transcription factors and components of the spliceosome and then influences both RNA polymerase II-dependent transcription and pre-mRNA splicing. The RTA-binding RRE motifs or specific motifs of unique cooperative transcription factors in viral and cellular DEG promoter regions indicate that BRLF1 employs different strategies for regulating viral and cellular transcription. Thus, our study characterized BRLF1-dependent cellular and viral transcriptional profile during primary infection and then revealed the comprehensive virus-cell interaction and alterations of transcription during EBV primary infection and lytic replication.

Determination of the dynamic cellular transcriptional profiles during kidney development from birth to maturity in rats by single-cell RNA sequencing

Recent single-cell RNA sequencing (scRNA-seq) analyses have offered much insight into the gene expression profiles in early-stage kidney development. However, comprehensive gene expression profiles from mid- and late-stage kidney development are lacking. In the present study, by using the scRNA-seq technique, we analyzed 54,704 rat kidney cells from just after birth to adulthood (six time points: postnatal days 0, 2, 5, 10, 20, and 56) including the mid and late stages of kidney development. Twenty-five original clusters and 13 different cell types were identified during these stages. Gene expression in these 13 cell types was mapped, and single cell atlas of the rat kidney from birth to maturity (http://youngbearlab.com) was built to enable users to search for a gene of interest and to evaluate its expression in different cells. The variation trend of six major types of kidney cells—intercalated cells of the collecting duct (CD-ICs), principal cells of the collecting duct (CD-PCs), cells of the distal convoluted tubules (DCTs), cells of the loop of Henle (LOH), podocytes (PDs), and cells of the proximal tubules (PTs)—during six postnatal time points was demonstrated. The trajectory of rat kidney development and the order of induction of the six major types of kidney cells from just after birth to maturity were determined. In addition, features of the dynamically changing genes as well as transcription factors during postnatal rat kidney development were identified. The present study provides a resource for achieving a deep understanding of the molecular basis of and regulatory events in the mid and late stages of kidney development.

Human splice factors contribute to latent HIV infection in primary cell models and blood CD4+ T cells from ART-treated individuals.

It is unclear what mechanisms govern latent HIV infection in vivo or in primary cell models. To investigate these questions, we compared the HIV and cellular transcription profile in three primary cell models and peripheral CD4+ T cells from HIV-infected ART-suppressed individuals using RT-ddPCR and RNA-seq. All primary cell models recapitulated the block to HIV multiple splicing seen in cells from ART-suppressed individuals, suggesting that this may be a key feature of HIV latency in primary CD4+ T cells. Blocks to HIV transcriptional initiation and elongation were observed more variably among models. A common set of 234 cellular genes, including members of the minor spliceosome pathway, was differentially expressed between unstimulated and activated cells from primary cell models and ART-suppressed individuals, suggesting these genes may play a role in the blocks to HIV transcription and splicing underlying latent infection. These genes may represent new targets for therapies designed to reactivate or silence latently-infected cells.

Up-regulation of hypoxia-inducible factor antisense as a novel approach to treat ovarian cancer.

Ovarian cancer (OC) is estimated to kill ~14,000 women in the United States in 2019. Current chemotherapies to treat OC initially show therapeutic efficacy but frequently drug resistance develops, at which point therapies with alternative targets are needed. Herein, we are describing a novel approach to sensitize these tumors to standard chemotherapies by increasing the transcription of hypoxia-inducible factor antisense.Methods: Genome-wide Bru-seq analysis was performed to fully capture the nascent transcriptional signature of OC cells treated with the gp130 inhibitor, SC144. In vitro and in vivo analysis, including characterization of hypoxia and select protein expression, combination with standard of care chemotherapy and antitumor efficacy were performed to assess the biological activity of SC144 on induction of hypoxia in OC cells.Results: Bru-seq analysis of OVCAR8 cells treated with SC144 shows upregulation of hypoxia related genes. In addition, transcription of hypoxia-inducible factor antisense (HIF1A-AS2) was induced that in turn reduced expression of HIF-1α and simultaneously increased expression of NDRG1. Furthermore, we observed decreased protein levels of EGFR, Met, c-Myc, cyclin D1, MMP-2, MMP-9 and TF, and phosphorylation of Src and P130-cas. SC144-induced alterations of HIF-1α and NDRG1 were also confirmed in prostate cancer cells. Ciclopirox olamine (CPX) induces a cellular transcriptional profile comparable to SC144, suggesting a similar cellular mechanism of action between these two compounds. In addition, SC144 sensitized OC cells to olaparib, carboplatin and cisplatin, and shows better in vivo efficacy than CPX.Conclusion: Induction of hypoxic stress responses through inhibition of gp130 represents a novel approach to design effective anticancer treatments in combination with standard-of-care chemotherapy in OC and the efficacy reported here strongly supports their clinical development.

Blinatumomab-Induced T Cell Activation at Single Cell Transcriptome Resolution

Introduction Blinatumomab, a bi-specific T cell engager (BiTE®) combining the VH and VL domains of two antibodies against human CD19 and CD3, has been approved by U.S. Food and Drug Administration (FDA) for the treatment of patients with relapsed or refractory B precursor ALL (r/r B-ALL) for its clinical benefit demonstrated in different clinical trials. Clinical trials have also shown that there are still patients refractory to blinatumomab. It is thus important to understand the resistance mechanisms. Blinatumomab connects patients' autologous T cells and target cells to form immunologic synapse which potently triggers the activation signaling cascades in T cells and guides T cells to recognize and induce perforin/granzyme-mediated lysis of CD19+ B-ALL cells. Previous studies showed blinatumomab-mediated cytotoxicity involves different T cell subpopulations. But response of each T cell subpopulation to blinatumomab treatment remained largely unknown. Methods and Results In this study, we used 10X Genomics based single cell RNA sequencing to analyze the transcriptome of single T cells before and after blinatumomab treatment. First, ex vivo blinatumomab cytotoxicity model was established, in which healthy PBMCs were used as effector cells and cocultured with target cells (RS4;11 cells or SUP-B15 cells) at an effector-to-target cell ratio of 10:1 with 0.1 ng/mL blinatumomab. Then, PBMCs and BMMCs from 2 B-ALL patients were cultured with 10 ng/mL blinatumomab. Cells from both ex vivo model and patient samples were sequenced using 10X Genomics platform. In total, transcriptome of 17920 single T cells from the ex vivo model and 2271 single T cells from patient sample were analyzed. Based on T cell trajectory analysis, we identified four distinct populations of blinatumomab-activated T cells, which were derived from CD8+ effector memory T (TEM) cells, CD4+ central memory (TCM) cells, naïve T cells and Tregs, respectively. The differentially expressed genes in activated clusters were analyzed to reflect T cell activation mechanisms. The result showed blinatumomab induced the upregulation of aerobic glycolysis pathway (PKM, PGAM1, ENO1, GAPDH and LDHA), cytoskeleton dynamics pathway (ACTD1, ACTB, NME1 and TUBA1B), IFN-responsive pathway (GBP1, PSME2, WARS, CXCL10 and STAT1), and the upregulation of well-known immune-related genes (TNFRSF4, TNFRSF18, LAG3, CD69, IL2RA, MIR155HG, BATF, SH2D2A, LTA, NFKBIA and NDFIP2). We found blinatumomab-activated CD8+ TEM cells showed stronger cytotoxic capability than other activated populations with specific production of cytotoxic factors (PRF1, IFNG and FASLG) and cytokines (CCL2, CCL3, CCL3L1, CCL4, CCL4L2, CCL8, XCL1, XCL2, TNFSF9 and TNFSF14). Last, differential gene expression analysis revealed that co-stimulatory (TNFRSF4,TNFRSF9 and TNFRSF18) and co-inhibitory receptors (LAG3 and TIGIT) were similarly up-regulated in clusters activated from memory and naïve T cells, indicating ligand dependent T cell functional outcomes induced by blinatumomab. Conclusion In summary, we used single cell sequencing to map the blinatumomab-mediated T cell activation state transition and reveal the molecular changes in different T cell subpopulations. Memory T cells, naïve T cells and Tregs were identified functional populations after blinatumomab treatment. CD8+ TEM accounted for the majority of blinatumomab-induced cytotoxicity. Furthermore, T cell co-regulatory receptors were identified as potential targets accountable for blinatumomab sensitivity or resistance mechanisms. The study demonstrated that the collected cellular transcriptional profiles can serve as resource to explore novel strategies to enhance the efficacy of blinatumomab. Disclosures Yin: Amgen: Employment. Huo:Amgen: Employment. Sheng:Amgen: Employment. Li:Amgen: Employment. Ellwanger:Amgen: Employment. Lu:Amgen: Employment. Homann:Amgen: Employment. Wang:Amgen: Employment. Ren:Ruijin hospital: Employment.

Heterogeneity in HIV and cellular transcription profiles in cell line models of latent and productive infection: implications for HIV latency

BackgroundHIV-infected cell lines are widely used to study latent HIV infection, which is considered the main barrier to HIV cure. We hypothesized that these cell lines differ from each other and from cells from HIV-infected individuals in the mechanisms underlying latency.ResultsTo quantify the degree to which HIV expression is inhibited by blocks at different stages of HIV transcription, we employed a recently-described panel of RT-ddPCR assays to measure levels of 7 HIV transcripts (“read-through,” initiated, 5′ elongated, mid-transcribed/unspliced [Pol], distal-transcribed [Nef], polyadenylated, and multiply-sliced [Tat-Rev]) in bulk populations of latently-infected (U1, ACH-2, J-Lat) and productively-infected (8E5, activated J-Lat) cell lines. To assess single-cell variation and investigate cellular genes associated with HIV transcriptional blocks, we developed a novel multiplex qPCR panel and quantified single cell levels of 7 HIV targets and 89 cellular transcripts in latently- and productively-infected cell lines. The bulk cell HIV transcription profile differed dramatically between cell lines and cells from ART-suppressed individuals. Compared to cells from ART-suppressed individuals, latent cell lines showed lower levels of HIV transcriptional initiation and higher levels of polyadenylation and splicing. ACH-2 and J-Lat cells showed different forms of transcriptional interference, while U1 cells showed a block to elongation. Single-cell studies revealed marked variation between/within cell lines in expression of HIV transcripts, T cell phenotypic markers, antiviral factors, and genes implicated in latency. Expression of multiply-spliced HIV Tat-Rev was associated with expression of cellular genes involved in activation, tissue retention, T cell transcription, and apoptosis/survival.ConclusionsHIV-infected cell lines differ from each other and from cells from ART-treated individuals in the mechanisms governing latent HIV infection. These differences in viral and cellular gene expression must be considered when gauging the suitability of a given cell line for future research on HIV. At the same time, some features were shared across cell lines, such as low expression of antiviral defense genes and a relationship between productive infection and genes involved in survival. These features may contribute to HIV latency or persistence in vivo, and deserve further study using novel single cell assays such as those described in this manuscript.

Epstein–Barr virus reprograms human B lymphocytes immediately in the prelatent phase of infection

Epstein-Barr virus (EBV) is a human tumor virus and a model of herpesviral latency. The virus efficiently infects resting human B lymphocytes and induces their continuous proliferation in vitro, which mimics certain aspects of EBV's oncogenic potential in vivo. How lymphoblastoid cell lines (LCLs) evolve from the infected lymphocytes is uncertain. We conducted a systematic time-resolved longitudinal study of cellular functions and transcriptional profiles of newly infected naïve primary B lymphocytes. EBV reprograms the cells comprehensively and globally. Rapid and extensive transcriptional changes occur within 24 h and precede any metabolic and phenotypic changes. Within 72 h, the virus activates the cells, changes their phenotypes with respect to cell size, RNA, and protein content, and induces metabolic pathways to cope with the increased demand for energy, supporting an efficient cell cycle entry on day 3 postinfection. The transcriptional program that EBV initiates consists of 3 waves of clearly discernable clusters of cellular genes that peak on day 2, 3, or 4 and regulate RNA synthesis, metabolic pathways, and cell division, respectively. Upon onset of cell doublings on day 4, the cellular transcriptome appears to be completely reprogrammed to support the proliferating cells, but 3 additional clusters of EBV-regulated genes fine-tune cell signaling, migration, and immune response pathways, eventually. Our study reveals that more than 11,000 genes are regulated upon EBV infection as naïve B cells exit quiescence to enter a germinal center-like differentiation program, which culminates in immortalized, proliferating cells that partially resemble plasmablasts and early plasma cells.

Evolution of metazoan oxygen-sensing involved a conserved divergence of VHL affinity for HIF1\u03b1 and HIF2\u03b1

Duplication of ancestral hypoxia-inducible factor (HIF)α coincided with the evolution of vertebrate species. Paralogs HIF1α and HIF2α are the most well-known factors for modulating the cellular transcriptional profile following hypoxia. However, how the processes of natural selection acted upon the coding region of these two genes to optimize the cellular response to hypoxia during evolution remains unclear. A key negative regulator of HIFα is von Hippel-Lindau (VHL) tumour suppressor protein. Here we show that evolutionarily-relevant substitutions can modulate a secondary contact between HIF1α Met561 and VHL Phe91. Notably, HIF1α binds more tightly than HIF2α to VHL due to a conserved Met to Thr substitution observed in the vertebrate lineage. Similarly, substitution of VHL Phe91 with Tyr, as seen in invertebrate species, decreases VHL affinity for both HIF1α and HIF2α. We propose that vertebrate evolution involved a more complex hypoxia response with fine-tuned divergence of VHL affinity for HIF1α and HIF2α.