Gene expression profiles during human CD4+ T cell differentiation

Slc43a2+ T cell metastasis from spleen to brain in RGNNV infected teleost

The origin of T cells in the teleost's brain is unclear. While viewing the central nervous system (CNS) as immune privileged has been widely accepted, previous studies suggest that T cells residing in the thymus but not in the spleen of the teleost play an essential role in communicating with the peripheral organs. Here, we identified nine T cell subpopulations in the thymus and spleen of orange-spotted grouper (Epinephelus coioices) through single-cell RNA-sequencing analysis. After viral CNS infection with red-spotted grouper nervous necrosis virus (RGNNV), the number of slc43a2+ T cells synchronously increased in the spleen and brain. During the infection tests in asplenic zebrafish (tlx1▲ zebrafish model), no increase in the number of slc43a2+ T cells was observed in the brain. Single-cell transcriptomic analysis indicated that slc43a2+ T cells mature and functionally differentiate within the spleen and then migrate into the brain to trigger an immune response. This study suggests a novel route for T cell migration from the spleen to the brain during viral infection in fish.

The gene expression profile and cell of origin of canine peripheral T-cell lymphoma

BACKGROUND

Peripheral T-cell lymphoma (PTCL) refers to a heterogenous group of T-cell neoplasms with poor treatment responses and survival times. Canine PTCL clinically and immunophenotypically resembles the most common human subtype, PTCL-not otherwise specified (PTCL-NOS), leading to interest in this canine disease as a naturally occurring model for human PTCL. Gene expression profiling in human PTCL-NOS has helped characterize this ambiguous diagnosis into distinct subtypes, but similar gene expression profiling in canine PTCL is lacking.

METHODS

Bulk RNA-sequencing was performed on tumor samples from 33 dogs with either CD4+ (26/33), CD8+ (4/33), or CD4-CD8- (3/33) PTCL as diagnosed by flow cytometry, and sorted CD4+ and CD8+ lymphocytes from healthy control dogs. Following normalization of RNA-seq data, we performed differential gene expression and unsupervised clustering methods. Gene set enrichment analysis was performed to determine the enrichment of canine CD4+ PTCL for human PTCL-NOS, oncogenic pathways, and various stages of T-cell development gene signatures. We utilized gene set variation analysis to evaluate individual canine CD4+ PTCLs for various human and murine T-cell and thymocyte gene signatures. Cultured canine PTCL cells were treated with a pan-PI3K inhibitor, and cell survival and proliferation were compared to DMSO-treated controls. Expression of GATA3 and phosphorylated AKT was validated by immunohistochemistry.

RESULTS

While the canine CD4+ PTCL phenotype exhibited a consistent gene expression profile, the expression profiles of CD8+ and CD4-CD8- canine PTCLs were more heterogeneous. Canine CD4+ PTCL had increased expression of GATA3, upregulation of its target genes, enrichment for PI3K/AKT/mTOR signaling, and downregulation of PTEN, features consistent with the more aggressive GATA3-PTCL subtype of human PTCL-NOS. In vitro assays validated the reliance of canine CD4+ PTCL cells on PI3K/AKT/mTOR signaling for survival and proliferation. Canine CD4+ PTCL was enriched for thymic precursor gene signatures, exhibited increased expression of markers of immaturity (CD34, KIT, DNTT, and CCR9), and downregulated genes associated with the T-cell receptor, MHC class II associated genes (DLA-DQA1, DLA-DRA, HLA-DQB1, and HLA-DQB2), and CD25.

CONCLUSIONS

Canine CD4+ PTCL most closely resembled the GATA3-PTCL subtype of PTCL-NOS and may originate from an earlier stage of T-cell development than the more conventionally posited mature T-helper cell origin.

DISCERN: deep single-cell expression reconstruction for improved cell clustering and cell subtype and state detection

BACKGROUND

Single-cell sequencing provides detailed insights into biological processes including cell differentiation and identity. While providing deep cell-specific information, the method suffers from technical constraints, most notably a limited number of expressed genes per cell, which leads to suboptimal clustering and cell type identification.

RESULTS

Here, we present DISCERN, a novel deep generative network that precisely reconstructs missing single-cell gene expression using a reference dataset. DISCERN outperforms competing algorithms in expression inference resulting in greatly improved cell clustering, cell type and activity detection, and insights into the cellular regulation of disease. We show that DISCERN is robust against differences between batches and is able to keep biological differences between batches, which is a common problem for imputation and batch correction algorithms. We use DISCERN to detect two unseen COVID-19-associated T cell types, cytotoxic CD4+ and CD8+ Tc2 T helper cells, with a potential role in adverse disease outcome. We utilize T cell fraction information of patient blood to classify mild or severe COVID-19 with an AUROC of 80% that can serve as a biomarker of disease stage. DISCERN can be easily integrated into existing single-cell sequencing workflow.

CONCLUSIONS

Thus, DISCERN is a flexible tool for reconstructing missing single-cell gene expression using a reference dataset and can easily be applied to a variety of data sets yielding novel insights, e.g., into disease mechanisms.

Transcriptional signatures associated with persisting CD19 CAR-T cells in children with leukemia

In the context of relapsed and refractory childhood pre-B cell acute lymphoblastic leukemia (R/R B-ALL), CD19-targeting chimeric antigen receptor (CAR)-T cells often induce durable remissions, which requires the persistence of CAR-T cells. In this study, we systematically analyzed CD19 CAR-T cells of 10 children with R/R B-ALL enrolled in the CARPALL trial via high-throughput single-cell gene expression and T cell receptor sequencing of infusion products and serial blood and bone marrow samples up to 5 years after infusion. We show that long-lived CAR-T cells developed a CD4/CD8 double-negative phenotype with an exhausted-like memory state and distinct transcriptional signature. This persistence signature was dominant among circulating CAR-T cells in all children with a long-lived treatment response for which sequencing data were sufficient (4/4, 100%). The signature was also present across T cell subsets and clonotypes, indicating that persisting CAR-T cells converge transcriptionally. This persistence signature was also detected in two adult patients with chronic lymphocytic leukemia with decade-long remissions who received a different CD19 CAR-T cell product. Examination of single T cell transcriptomes from a wide range of healthy and diseased tissues across children and adults indicated that the persistence signature may be specific to long-lived CAR-T cells. These findings raise the possibility that a universal transcriptional signature of clinically effective, persistent CD19 CAR-T cells exists.

Multimodal human thymic profiling reveals trajectories and cellular milieu for T agonist selection

To prevent autoimmunity, thymocytes expressing self-reactive T cell receptors (TCRs) are negatively selected, however, divergence into tolerogenic, agonist selected lineages represent an alternative fate. As thymocyte development, selection, and lineage choices are dependent on spatial context and cell-to-cell interactions, we have performed Cellular Indexing of Transcriptomes and Epitopes by sequencing (CITE-seq) and spatial transcriptomics on paediatric human thymu​​s. Thymocytes expressing markers of strong TCR signalling diverged from the conventional developmental trajectory prior to CD4⁺ or CD8⁺ lineage commitment, while markers of different agonist selected T cell populations (CD8αα(I), CD8αα(II), T_(agonist), T_reg(diff), and T_reg) exhibited variable timing of induction. Expression profiles of chemokines and co-stimulatory molecules, together with spatial localisation, supported that dendritic cells, B cells, and stromal cells contribute to agonist selection, with different subsets influencing thymocytes at specific developmental stages within distinct spatial niches. Understanding factors influencing agonist T cells is needed to benefit from their immunoregulatory effects in clinical use.

CCL2: a Chemokine Potentially Promoting Early Seeding of the Latent HIV Reservoir

HIV infects long-lived CD4 memory T cells, establishing a latent viral reservoir that necessitates lifelong antiretroviral therapy (ART). How this reservoir is formed so quickly after infection remains unclear. We now show the innate inflammatory response to HIV infection results in CCL2 chemokine release, leading to recruitment of cells expressing the CCR2 receptor, including a subset of central memory CD4 T cells. Supporting a role for the CCL2/CCR2 axis in rapid reservoir formation, we find (i) treatment of humanized mice with anti-CCL2 antibodies during early HIV infection decreases reservoir seeding and preserves CCR2/5+ cells and (ii) CCR2/5+ cells from the blood of HIV-infected individuals on long-term ART contain significantly more integrated provirus than CCR2/5-negative memory or naive cells. Together, these studies support a model where the host's innate inflammatory response to HIV infection, including CCL2 production, leads to the recruitment of CCR2/5+ central memory CD4 T cells to zones of virus-associated inflammation, likely contributing to rapid formation of the latent HIV reservoir. IMPORTANCE There are currently over 35 million people living with HIV worldwide, and we still have no vaccine or scalable cure. One of the difficulties with HIV is its ability to rapidly establish a viral reservoir in lymphoid tissues that allows it to elude antivirals and the immune system. Thus, it is important to understand how HIV accomplishes this so we can develop preventive strategies. Our current results show that an early inflammatory response to HIV infection includes production of the chemokine CCL2, which recruits a unique subset of CCR2/5+ CD4+ T cells that become infected and form a significant reservoir for latent infection. Furthermore, we show that blockade of CCL2 in humanized mice significantly reduces persistent HIV infection. This information is relevant to the development of therapeutics to prevent and/or treat chronic HIV infections.

Sequential Single-Cell Transcriptional and Protein Marker Profiling Reveals TIGIT as a Marker of CD19 CAR-T Cell Dysfunction in Patients with Non-Hodgkin Lymphoma

Chimeric antigen receptor T-cell (CAR-T cell) therapy directed at CD19 produces durable remissions in the treatment of relapsed/refractory non-Hodgkin lymphoma (NHL). Nonetheless, many patients receiving CD19 CAR-T cells fail to respond for unknown reasons. To reveal changes in 4-1BB-based CD19 CAR-T cells and identify biomarkers of response, we used single-cell RNA sequencing and protein surface marker profiling of patient CAR-T cells pre- and postinfusion into patients with NHL. At the transcriptional and protein levels, we note the evolution of CAR-T cells toward a nonproliferative, highly differentiated, and exhausted state, with an enriched exhaustion profile in CAR-T cells of patients with poor response marked by TIGIT expression. Utilizing in vitro and in vivo studies, we demonstrate that TIGIT blockade alone improves the antitumor function of CAR-T cells. Altogether, we provide evidence of CAR-T cell dysfunction marked by TIGIT expression driving a poor response in patients with NHL.

SIGNIFICANCE

This is the first study investigating the mechanisms linked to CAR-T patient responses based on the sequential analysis of manufactured and infused CAR-T cells using single-cell RNA and protein expression data. Furthermore, our findings are the first to demonstrate an improvement of CAR-T cell efficacy with TIGIT inhibition alone. This article is highlighted in the In This Issue feature, p. 1825.

Assessing Consistency Across Functional Screening Datasets in Cancer Cells

Motivation

Many high-throughput screening studies have been carried out in cancer cell lines to identify therapeutic agents and targets. Existing consistency assessment studies only examined two datasets at a time, with conclusions based on a subset of carefully selected features rather than considering global consistency of all the data. However, poor concordance can still be observed for a large part of the data even when selected features are highly consistent.

Results

In this study, we assembled nine compound screening datasets and three functional genomics datasets. We derived direct measures of consistency as well as indirect measures of consistency based on association between functional data and copy number-adjusted gene expression data. These results have been integrated into a web application—the Functional Data Consistency Explorer (FDCE), to allow users to make queries and generate interactive visualizations so that functional data consistency can be assessed for individual features of interest.

Availability and implementation

The FDCE web tool and we have developed and the functional data consistency measures we have generated are available at https://lccl.shinyapps.io/FDCE/.

Supplementary information

Supplementary data are available at Bioinformatics online.

HCV eradication with IFN-based therapy does not completely restore gene expression in PBMCs from HIV/HCV-coinfected patients

Objective

To evaluate the impact of hepatitis C virus (HCV) elimination via interferon (IFN)-based therapy on gene expression profiles related to the immune system in HIV/HCV-coinfected patients.

Methods

We conducted a prospective study in 28 HIV/HCV-coinfected patients receiving IFN-based therapy at baseline (HIV/HCV-b) and week 24 after sustained virological response (HIV/HCV-f). Twenty-seven HIV-monoinfected patients (HIV-mono) were included as a control. RNA-seq analysis was performed on peripheral blood mononuclear cells (PBMCs). Genes with a fold-change (FC) ≥ 1.5 (in either direction) and false discovery rate (FDR) ≤ 0.05 were identified as significantly differentially expressed (SDE).

Results

HIV/HCV-b showed six SDE genes compared to HIV-mono group, but no significantly enriched pathways were observed. For HIV/HCV-f vs. HIV/HCV-b, we found 58 SDE genes, 34 upregulated and 24 downregulated in the HIV/HCV-f group. Of these, the most overexpressed were CXCL2, PDCD6IP, ATP5B, IGSF9, RAB26, and CSRNP1, and the most downregulated were IFI44 and IFI44L. These 58 SDE genes revealed two significantly enriched pathways (FDR < 0.05), one linked to Epstein-Barr virus infection and another related to p53 signaling. For HIV/HCV-f vs. HIV-mono group, we found 44 SDE genes that revealed 31 enriched pathways (FDR < 0.05) related to inflammation, cancer/cell cycle alteration, viral and bacterial infection, and comorbidities associated with HIV/HCV-coinfection. Five genes were overrepresented in most pathways (JUN, NFKBIA, PIK3R2, CDC42, and STAT3).

Conclusion

HIV/HCV-coinfected patients who eradicated hepatitis C with IFN-based therapy showed profound gene expression changes after achieving sustained virological response. The altered pathways were related to inflammation and liver-related complications, such as non-alcoholic fatty liver disease and hepatocellular carcinoma, underscoring the need for active surveillance for these patients.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12929-021-00718-6.

ShinyButchR: Interactive NMF-based decomposition workflow of genome-scale datasets

Non-negative matrix factorization (NMF) has been widely used for the analysis of genomic data to perform feature extraction and signature identification due to the interpretability of the decomposed signatures. However, running a basic NMF analysis requires the installation of multiple tools and dependencies, along with a steep learning curve and computing time. To mitigate such obstacles, we developed ShinyButchR, a novel R/Shiny application that provides a complete NMF-based analysis workflow, allowing the user to perform matrix decomposition using NMF, feature extraction, interactive visualization, relevant signature identification, and association to biological and clinical variables. ShinyButchR builds upon the also novel R package ButchR, which provides new TensorFlow solvers for algorithms of the NMF family, functions for downstream analysis, a rational method to determine the optimal factorization rank and a novel feature selection strategy.

International Union of Basic and Clinical Pharmacology. CX. Classification of Receptors for 5-hydroxytryptamine; Pharmacology and Function

5-HT receptors expressed throughout the human body are targets for established therapeutics and various drugs in development. Their diversity of structure and function reflects the important role 5-HT receptors play in physiologic and pathophysiological processes. The present review offers a framework for the official receptor nomenclature and a detailed understanding of each of the 14 5-HT receptor subtypes, their roles in the systems of the body, and, where appropriate, the (potential) utility of therapeutics targeting these receptors. SIGNIFICANCE STATEMENT: This review provides a comprehensive account of the classification and function of 5-hydroxytryptamine receptors, including how they are targeted for therapeutic benefit.

ZipSeq: barcoding for real-time mapping of single cell transcriptomes

Spatial transcriptomics seeks to integrate single cell transcriptomic data within the three-dimensional space of multicellular biology. Current methods to correlate a cell's position with its transcriptome in living tissues have various limitations. We developed an approach, called 'ZipSeq', that uses patterned illumination and photocaged oligonucleotides to serially print barcodes ('zipcodes') onto live cells in intact tissues, in real time and with an on-the-fly selection of patterns. Using ZipSeq, we mapped gene expression in three settings: in vitro wound healing, live lymph node sections and a live tumor microenvironment. In all cases, we discovered new gene expression patterns associated with histological structures. In the tumor microenvironment, this demonstrated a trajectory of myeloid and T cell differentiation from the periphery inward. A combinatorial variation of ZipSeq efficiently scales in the number of regions defined, providing a pathway for complete mapping of live tissues, subsequent to real-time imaging or perturbation.

The Brain-Skin Connection and the Pathogenesis of Psoriasis: A Review with a Focus on the Serotonergic System

Psoriasis is a common non-communicable chronic immune-mediated skin disease, affecting approximately 125 million people in the world. Its pathogenesis results from a combination of genetic and environmental factors. The pathogenesis of psoriasis seems to be driven by the interaction between innate immune cells, adaptive immune cells and keratinocytes, in a process mediated by cytokines (including interleukins (IL)-6, IL-17 and IL-22, interferon and tumor necrosis factor) and other signaling molecules. This leads to an inflammatory process with increased proliferation of epidermal cells, neo-angiogenesis and infiltration of dendritic cells in the skin. Dysfunctional de novo glucocorticoid synthesis in psoriatic keratinocytes and the skin microbiome have also been suggested as mediators in the pathogenesis of this disease. To understand psoriasis, it is essential to comprehend the processes underlying the skin immunity and neuroendocrinology. This review paper focuses on the skin as a neuroendocrine organ and summarizes what is known about the skin immune system, the brain-skin connection and the role played by the serotonergic system in skin. Subsequently, the alterations of neuroimmune processes and of the serotonergic system in psoriatic skin are discussed, as well as, briefly, the genetic basis of psoriasis.

Glucocorticoids paradoxically facilitate steroid resistance in T cell acute lymphoblastic leukemias and thymocytes

Glucocorticoids (GCs) are a central component of therapy for patients with T cell acute lymphoblastic leukemia (T-ALL), and although resistance to GCs is a strong negative prognostic indicator in T-ALL, the mechanisms of GC resistance remain poorly understood. Using diagnostic samples from patients enrolled in the frontline Children's Oncology Group (COG) T-ALL clinical trial AALL1231, we demonstrated that one-third of primary T-ALLs were resistant to GCs when cells were cultured in the presence of IL-7, a cytokine that is critical for normal T cell function and that plays a well-established role in leukemogenesis. We demonstrated that in these T-ALLs and in distinct populations of normal developing thymocytes, GCs paradoxically induced their own resistance by promoting upregulation of IL-7 receptor (IL-7R) expression. In the presence of IL-7, this augmented downstream signal transduction, resulting in increased STAT5 transcriptional output and upregulation of the prosurvival protein BCL-2. Taken together, we showed that IL-7 mediates an intrinsic and physiologic mechanism of GC resistance in normal thymocyte development that is retained during leukemogenesis in a subset of T-ALLs and is reversible with targeted inhibition of the IL-7R/JAK/STAT5/BCL-2 axis.

An ERG Enhancer-Based Reporter Identifies Leukemia Cells with Elevated Leukemogenic Potential Driven by ERG-USP9X Feed-Forward Regulation

Acute leukemia is a rapidly progressing blood cancer with low survival rates. Unfavorable prognosis is attributed to insufficiently characterized subpopulations of leukemia stem cells (LSC) that drive chemoresistance and leukemia relapse. Here we utilized a genetic reporter that assesses stemness to enrich and functionally characterize LSCs. We observed heterogeneous activity of the ERG+85 enhancer-based fluorescent reporter in human leukemias. Cells with high reporter activity (tagBFP^High) exhibited elevated expression of stemness and chemoresistance genes and demonstrated increased clonogenicity and resistance to chemo- and radiotherapy as compared with their tagBFP^Neg counterparts. The tagBFP^High fraction was capable of regenerating the original cellular heterogeneity and demonstrated increased invasive ability. Moreover, the tagBFP^High fraction was enriched for leukemia-initiating cells in a xenograft assay. We identified the ubiquitin hydrolase USP9X as a novel ERG transcriptional target that sustains ERG+85-positive cells by controlling ERG ubiquitination. Therapeutic targeting of USP9X led to preferential inhibition of the ERG-dependent leukemias. Collectively, these results characterize human leukemia cell functional heterogeneity and suggest that targeting ERG via USP9X inhibition may be a potential treatment strategy in patients with leukemia. SIGNIFICANCE: This study couples a novel experimental tool with state-of-the-art approaches to delineate molecular mechanisms underlying stem cell-related characteristics in leukemia cells.

Activation of Th1 Immunity within the Tumor Microenvironment Is Associated with Clinical Response to Lenalidomide in Chronic Lymphocytic Leukemia

Immune stimulation contributes to lenalidomide's antitumor activity. Chronic lymphocytic leukemia (CLL) is characterized by the accumulation of mature, autoreactive B cells in secondary lymphoid tissues, blood, and bone marrow and progressive immune dysfunction. Previous studies in CLL indicated that lenalidomide can repair defective T cell function in vitro. Whether T cell activation is required for clinical response to lenalidomide remains unclear. In this study, we report changes in the immune microenvironment in patients with CLL treated with single-agent lenalidomide and associate the immunologic effects of lenalidomide with antitumor response. Within days of starting lenalidomide, T cells increased in the tumor microenvironment and showed Th1-type polarization. Gene expression profiling of pretreatment and on-treatment lymph node biopsy specimens revealed upregulation of IFN-γ and many of its target genes in response to lenalidomide. The IFN-γ-mediated Th1 response was limited to patients achieving a clinical response defined by a reduction in lymphadenopathy. Deep sequencing of TCR genes revealed decreasing diversity of the T cell repertoire and an expansion of select clonotypes in responders. To validate our observations, we stimulated T cells and CLL cells with lenalidomide in culture and detected lenalidomide-dependent increases in T cell proliferation. Taken together, our data demonstrate that lenalidomide induced Th1 immunity in the lymph node that is associated with clinical response.

Classification of gene signatures for their information value and functional redundancy

Gene signatures are more and more used to interpret results of omics data analyses but suffer from compositional (large overlap) and functional (correlated read-outs) redundancy. Moreover, many gene signatures rarely come out as significant in statistical tests. Based on pan-cancer data analysis, we construct a restricted set of 962 signatures defined as informative and demonstrate that they have a higher probability to appear enriched in comparative cancer studies. We show that the majority of informative signatures conserve their weights for the genes composing the signature (eigengenes) from one cancer type to another. We finally construct InfoSigMap, an interactive online map of these signatures and their cross-correlations. This map highlights the structure of compositional and functional redundancies between informative signatures, and it charts the territories of biological functions. InfoSigMap can be used to visualize the results of omics data analyses and suggests a rearrangement of existing gene sets.

An informative collection of gene signatures for transcriptomic data analysis is constructed. The number of transcriptomic signatures grows fast and their collections are highly redundant that hampers omics data analyses interpretation. A computational biology team from Institut Curie led by Andrei Zinovyev selected a collection of 962 gene signatures shown to be informative for cancer studies and reflecting mechanisms of cancer progression. The signatures were filtered from a large compendium without requiring any manual curation by experts through a large-scale unbiased analysis of pancancer data. They have much higher chance to obtain significant enrichment scores in a comparative trancriptomic study. The authors integrated the 962 signatures into InfoSigMap, a new data visualization resource for the interpretation of the results of omics data analyses, which facilitates getting an insight into the mechanisms driving cancer.

Clinical Evidence Supports a Protective Role for CXCL5 in Coronary Artery Disease

Our goal was to measure the association of CXCL5 and molecular phenotypes associated with coronary atherosclerosis severity in patients at least 65 years old. CXCL5 is classically defined as a proinflammatory chemokine, but its role in chronic inflammatory diseases, such as coronary atherosclerosis, is not well defined. We enrolled individuals who were at least 65 years old and undergoing diagnostic cardiac catheterization. Coronary artery disease (CAD) severity was quantified in each subject via coronary angiography by calculating a CAD score. Circulating CXCL5 levels were measured from plasma, and both DNA genotyping and mRNA expression levels in peripheral blood mononuclear cells were quantified via microarray gene chips. We observed a negative association of CXCL5 levels with CAD at an odds ratio (OR) of 0.46 (95% CI, 0.27-0.75). Controlling for covariates, including sex, statin use, hypertension, hyperlipidemia, obesity, self-reported race, smoking, and diabetes, the OR was not significantly affected [OR, 0.54 (95% CI, 0.31-0.96)], consistent with a protective role for CXCL5 in coronary atherosclerosis. We also identified 18 genomic regions with expression quantitative trait loci of genes correlated with both CAD severity and circulating CXCL5 levels. Our clinical findings are consistent with the emerging link between chemokines and atherosclerosis and suggest new therapeutic targets for CAD.

Immune activation is associated with decreased thymic function in asymptomatic, untreated HIV-infected individuals

Background

Reduced thymic function causes poor immunological reconstitution in human immunodeficiency virus (HIV)-positive patients on combined antiretroviral therapy (cART). The association between immune activation and thymic function in asymptomatic HIV-positive treatment-naive individuals has thus far not been investigated.

Aims and objectives

To optimise a five-colour flow cytometric assay for measurement of thymic function by measuring recent thymic emigrants (RTEs) in treatment-naive HIV-infected patients and healthy controls and correlate results with levels of immune activation, CD4 counts and viral load.

Methods

Blood obtained from 53 consenting HIV-positive individuals and 32 controls recruited from HIV prevention and testing clinic in Cape Town, South Africa. RTEs were measured (CD3⁺/CD4⁺/CD45RA⁺/CD31⁺/CD62L⁺) and levels were correlated with CD4 counts of HIV-infected individuals, log viral load and levels of immune activation (CD8⁺/CD38⁺ T-cells).

Results

HIV-infected individuals had reduced frequencies of RTEs when compared to controls (p = 0.0035). Levels of immune activation were inversely correlated with thymic function (p = 0.0403), and the thymic function in HIV-infected individuals showed no significant correlation with CD4 counts (p = 0.31559) and viral load (p = 0.20628).

Conclusions

There was impaired thymic function in HIV-infected individuals, which was associated with increased levels of immune activation. The thymic dysfunction was not associated with CD4 counts and viral load. Immune activation may result in inflammatory damage to the thymus and subsequent thymic dysfunction, and CD4 counts and viral load may not necessarily reflect thymic dysfunction in HIV.

A curated transcriptome dataset collection to investigate the functional programming of human hematopoietic cells in early life

Compendia of large-scale datasets made available in public repositories provide an opportunity to identify and fill gaps in biomedical knowledge. But first, these data need to be made readily accessible to research investigators for interpretation. Here we make available a collection of transcriptome datasets to investigate the functional programming of human hematopoietic cells in early life. Thirty two datasets were retrieved from the NCBI Gene Expression Omnibus (GEO) and loaded in a custom web application called the Gene Expression Browser (GXB), which was designed for interactive query and visualization of integrated large-scale data. Quality control checks were performed. Multiple sample groupings and gene rank lists were created allowing users to reveal age-related differences in transcriptome profiles, changes in the gene expression of neonatal hematopoietic cells to a variety of immune stimulators and modulators, as well as during cell differentiation. Available demographic, clinical, and cell phenotypic information can be overlaid with the gene expression data and used to sort samples. Web links to customized graphical views can be generated and subsequently inserted in manuscripts to report novel findings. GXB also enables browsing of a single gene across projects, thereby providing new perspectives on age- and developmental stage-specific expression of a given gene across the human hematopoietic system. This dataset collection is available at: http://developmentalimmunology.gxbsidra.org/dm3/geneBrowser/list.

Complexity and Specificity of the Neutrophil Transcriptomes in Juvenile Idiopathic Arthritis

NIH projects such as ENCODE and Roadmap Epigenomics have revealed surprising complexity in the transcriptomes of mammalian cells. In this study, we explored transcriptional complexity in human neutrophils, cells generally regarded as nonspecific in their functions and responses. We studied distinct human disease phenotypes and found that, at the gene, gene isoform, and miRNA level, neutrophils exhibit considerable specificity in their transcriptomes. Thus, even cells whose responses are considered non-specific show tailoring of their transcriptional repertoire toward specific physiologic or pathologic contexts. We also found that miRNAs had a global impact on neutrophil transcriptome and are associated with innate immunity in juvenile idiopathic arthritis (JIA). These findings have important implications for our understanding of the link between genes, non-coding transcripts and disease phenotypes.

Long Noncoding RNA Expression during Human B-Cell Development

Long noncoding RNAs (lncRNAs) have emerged as important regulators of diverse cellular processes, but their roles in the developing immune system are poorly understood. In this study, we analysed lncRNA expression during human B-cell development by array-based expression profiling of eleven distinct flow-sorted B-cell subsets, comprising pre-B1, pre-B2, immature, naive, memory, and plasma cells from bone marrow biopsies (n = 7), and naive, centroblast, centrocyte, memory, and plasmablast cells from tonsil tissue samples (n = 6), respectively. A remapping strategy was used to assign the array probes to 37630 gene-level probe sets, reflecting recent updates in genomic and transcriptomic databases, which enabled expression profiling of 19579 long noncoding RNAs, comprising 3947 antisense RNAs, 5277 lincRNAs, 7625 pseudogenes, and 2730 additional lncRNAs. As a first step towards inferring the functions of the identified lncRNAs in developing B-cells, we analysed their co-expression with well-characterized protein-coding genes, a method known as “guilt by association”. By using weighted gene co-expression network analysis, we identified 272 lincRNAs, 471 antisense RNAs, 376 pseudogene RNAs, and 64 lncRNAs within seven sub-networks associated with distinct stages of B-cell development, such as early B-cell development, B-cell proliferation, affinity maturation of antibody, and terminal differentiation. These data provide an important resource for future studies on the functions of lncRNAs in development of the adaptive immune response, and the pathogenesis of B-cell malignancies that originate from distinct B-cell subpopulations.

Neonatal regulatory T cells have reduced capacity to suppress dendritic cell function

Regulatory T cells (Treg cells) limit contact between dendritic cells (DCs) and conventional T cells (Tcons), decreasing the formation of aggregates as well as down-modulating the expression of co-stimulatory molecules by DCs, thus decreasing DC immunogenicity and abrogating T-cell activation. Despite the importance of this Treg-cell function, the capacity of Treg cells from term and preterm neonates to suppress DCs, and the suppressive mechanisms they use, are still undefined. We found that, relative to adult Treg cells, activated Treg cells from human neonates expressed lower FOXP3 and CTLA-4, but contained higher levels of cAMP. We developed an in vitro model in which Treg function was measured at a physiological ratio of 1 Treg for 10 Tcon and 1 monocyte-derived DC, as Treg target. Term and preterm Treg cells failed to suppress the formation of DC-Tcon aggregates, in contrast to naïve and memory Treg cells from adults. However, neonatal Treg cells diminished DC and Tcon activation as well as actin polymerization at the immunological synapses. In addition, CTLA-4 and cAMP were the main suppressive molecules used by neonatal Treg. Altogether, both preterm and term neonatal Treg cells appear less functional than adult Treg cells, and this defect is consistent with the general impairment of CD4 cell function in newborns.

A common gene signature across multiple studies relate biomarkers and functional regulation in tolerance to renal allograft

Patients tolerant to a kidney graft display a specific blood cell transcriptional pattern but results from five different studies were inconsistent, raising the question of relevance for future clinical application. To resolve this, we sought to identify a common gene signature, specific functional and cellular components, and discriminating biomarkers for tolerance following kidney transplantation. A meta-analysis of studies identified a robust gene signature involving proliferation of B and CD4 T cells, and inhibition of CD14 monocyte related functions among 96 tolerant samples. This signature was further supported through a cross-validation approach, yielding 92.5% accuracy independent of the study of origin. Experimental validation, performed on new tolerant samples and using a selection of the top-20 biomarkers, returned 91.7% of good classification. Beyond the confirmation of B-cell involvement, our data also indicated participation of other cell subsets in tolerance. Thus, the use of the top 20 biomarkers, mostly centered on B cells, may provide a common and standardized tool towards personalized medicine for the monitoring of tolerant or low-risk patients among kidney allotransplant recipients. These data point to a global preservation of genes favoring the maintenance of a homeostatic and ‘healthy' environment in tolerant patients and may contribute to a better understanding of tolerance maintenance mechanisms.

The H3K27me3 demethylase UTX is a gender-specific tumor suppressor in T-cell acute lymphoblastic leukemia

T-cell acute lymphoblastic leukemia (T-ALL) is an aggressive form of leukemia that is mainly diagnosed in children and shows a skewed gender distribution toward males. In this study, we report somatic loss-of-function mutations in the X-linked histone H3K27me3 demethylase ubiquitously transcribed X (UTX) chromosome, in human T-ALL. Interestingly, UTX mutations were exclusively present in male T-ALL patients and allelic expression analysis revealed that UTX escapes X-inactivation in female T-ALL lymphoblasts and normal T cells. Notably, we demonstrate in vitro and in vivo that the H3K27me3 demethylase UTX functions as a bona fide tumor suppressor in T-ALL. Moreover, T-ALL driven by UTX inactivation exhibits collateral sensitivity to pharmacologic H3K27me3 inhibition. All together, our results show how a gender-specific and therapeutically relevant defect in balancing H3K27 methylation contributes to T-cell leukemogenesis.

NUP98-PHF23 is a chromatin-modifying oncoprotein that causes a wide array of leukemias sensitive to inhibition of PHD histone reader function

In this report, we show that expression of a NUP98-PHF23 (NP23) fusion, associated with acute myeloid leukemia (AML) in humans, leads to myeloid, erythroid, T-cell, and B-cell leukemia in mice. The leukemic and preleukemic tissues display a stem cell-like expression signature, including Hoxa, Hoxb, and Meis1 genes. The PHF23 plant homeodomain (PHD) motif is known to bind to H3K4me3 residues, and chromatin immunoprecipitation experiments demonstrated that the NP23 protein binds to chromatin at a specific subset of H3K4me3 sites, including at Hoxa, Hoxb, and Meis1. Treatment of NP23 cells with disulfiram, which inhibits the binding of PHD motifs to H3K4me3, rapidly and selectively killed NP23-expressing myeloblasts; cell death was preceded by decreased expression of Hoxa, Hoxb, and Meis1. Furthermore, AML driven by a related fusion gene, NUP98-JARID1A (NJL), was also sensitive to disulfiram. Thus, the NP23 mouse provides a platform to evaluate compounds that disrupt binding of oncogenic PHD proteins to H3K4me3.

The ubiquitin ligase FBXW7 modulates leukemia-initiating cell activity by regulating MYC stability

Sequencing efforts led to the identification of somatic mutations that could affect the self-renewal and differentiation of cancer-initiating cells. One such recurrent mutation targets the binding pocket of the ubiquitin ligase Fbxw7. Missense FBXW7 mutations are prevalent in various tumors, including T cell acute lymphoblastic leukemia (T-ALL). To study the effects of such lesions, we generated animals carrying regulatable Fbxw7 mutant alleles. Here, we show that these mutations specifically bolster cancer-initiating cell activity in collaboration with Notch1 oncogenes but spare normal hematopoietic stem cell function. We were also able to show that FBXW7 mutations specifically affect the ubiquitylation and half-life of c-Myc protein, a key T-ALL oncogene. Using animals carrying c-Myc fusion alleles, we connected Fbxw7 function to c-Myc abundance and correlated c-Myc expression to leukemia-initiating activity. Finally, we demonstrated that small-molecule-mediated suppression of MYC activity leads to T-ALL remission, suggesting an effective therapeutic strategy.

CD4⁺ follicular helper T cell infiltration predicts breast cancer survival

CD4⁺ T cells are critical regulators of immune responses, but their functional role in human breast cancer is relatively unknown. The goal of this study was to produce an image of CD4⁺ T cells infiltrating breast tumors using limited ex vivo manipulation to better understand the in vivo differences associated with patient prognosis. We performed comprehensive molecular profiling of infiltrating CD4⁺ T cells isolated from untreated invasive primary tumors and found that the infiltrating T cell subpopulations included follicular helper T (Tfh) cells, which have not previously been found in solid tumors, as well as Th1, Th2, and Th17 effector memory cells and Tregs. T cell signaling pathway alterations included a mixture of activation and suppression characterized by restricted cytokine/chemokine production, which inversely paralleled lymphoid infiltration levels and could be reproduced in activated donor CD4⁺ T cells treated with primary tumor supernatant. A comparison of extensively versus minimally infiltrated tumors showed that CXCL13-producing CD4⁺ Tfh cells distinguish extensive immune infiltrates, principally located in tertiary lymphoid structure germinal centers. An 8-gene Tfh signature, signifying organized antitumor immunity, robustly predicted survival or preoperative response to chemotherapy. Our identification of CD4⁺ Tfh cells in breast cancer suggests that they are an important immune element whose presence in the tumor is a prognostic factor.

Genome-wide profiling of in vivo LPS-responsive genes in splenic myeloid cells

Lipopolysaccharide (LPS), the major causative agent of bacterial sepsis, has been used by many laboratories in genome-wide expression profiling of the LPS response. However, these studies have predominantly used in vitro cultured macrophages (Macs), which may not accurately reflect the LPS response of these innate immune cells in vivo. To overcome this limitation and to identify inflammatory genes in vivo, we have profiled genome-wide expression patterns in non-lymphoid, splenic myeloid cells extracted directly from LPS-treated mice. Genes encoding factors known to be involved in mediating or regulating inflammatory processes, such as cytokines and chemokines, as well as many genes whose immunological functions are not well known, were strongly induced by LPS after 3 h or 8 h of treatment. Most of the highly LPS-responsive genes that we randomly selected from the microarray data were independently confirmed by quantitative RT-PCR, implying that our microarray data are quite reliable. When our in vivo data were compared to previously reported microarray data for in vitro LPS-treated Macs, a significant proportion (∼20%) of the in vivo LPS-responsive genes defined in this study were specific to cells exposed to LPS in vivo, but a larger proportion of them (∼60%) were influenced by LPS in both in vitro and in vivo settings. This result indicates that our in vivo LPS-responsive gene set includes not only previously identified in vitro LPS-responsive genes but also novel LPS-responsive genes. Both types of genes would be a valuable resource in the future for understanding inflammatory responses in vivo.

Biomolecular events in cancer revealed by attractor metagenes

Mining gene expression profiles has proven valuable for identifying signatures serving as surrogates of cancer phenotypes. However, the similarities of such signatures across different cancer types have not been strong enough to conclude that they represent a universal biological mechanism shared among multiple cancer types. Here we present a computational method for generating signatures using an iterative process that converges to one of several precise attractors defining signatures representing biomolecular events, such as cell transdifferentiation or the presence of an amplicon. By analyzing rich gene expression datasets from different cancer types, we identified several such biomolecular events, some of which are universally present in all tested cancer types in nearly identical form. Although the method is unsupervised, we show that it often leads to attractors with strong phenotypic associations. We present several such multi-cancer attractors, focusing on three that are prominent and sharply defined in all cases: a mesenchymal transition attractor strongly associated with tumor stage, a mitotic chromosomal instability attractor strongly associated with tumor grade, and a lymphocyte-specific attractor.

Cancer is known to be characterized by several unifying biological capabilities or “hallmarks.” However, attempts to computationally identify patterns, such as gene expression signatures, shared across many different cancer types have been largely unsuccessful. A typical approach has been to classify samples into mutually exclusive subtypes, each of which is characterized by a particular gene signature. Although occasional similarities of such signatures in different cancer types exist, these similarities have not been sufficiently strong to conclude that they reflect the same biological event. By contrast, we have developed a computational methodology that has identified some signatures of co-expressed genes exhibiting remarkable similarity across many different cancer types. These signatures appear as stable “attractors” of an iterative computational procedure that tends to collect mutually associated genes, so that its convergence can point to the core (“heart”) of the underlying biological co-expression mechanism. One of these “pan-cancer” attractors corresponds to a transdifferentiation of cancer cells empowering them with invasiveness and motility. Another represents a mitotic chromosomal instability of cancer cells. A third attractor is lymphocyte-specific.

CRNDE: A Long Non-Coding RNA Involved in CanceR, Neurobiology, and DEvelopment

CRNDE is the gene symbol for Colorectal Neoplasia Differentially Expressed (non-protein-coding), a long non-coding RNA (lncRNA) gene that expresses multiple splice variants and displays a very tissue-specific pattern of expression. CRNDE was initially identified as a lncRNA whose expression is highly elevated in colorectal cancer, but it is also upregulated in many other solid tumors and in leukemias. Indeed, CRNDE is the most upregulated lncRNA in gliomas and here, as in other cancers, it is associated with a "stemness" signature. CRNDE is expressed in specific regions within the human and mouse brain; the mouse ortholog is high in induced pluripotent stem cells and increases further during neuronal differentiation. We suggest that CRNDE is a multifunctional lncRNA whose different splice forms provide specific functional scaffolds for regulatory complexes, such as the polycomb repressive complex 2 (PRC2) and CoREST chromatin-modifying complexes, which CRNDE helps pilot to target genes.

New molecular insights into peripheral T cell lymphomas

Peripheral T cell lymphomas (PTCLs) are heterogeneous neoplasms and represent about 12% of all lymphoid malignancies. They are often regarded as "orphan diseases," a designation that does not reflect their real incidence but rather signifies the difficulties encountered in their classification, diagnosis, and treatment. Here we revise the current understanding of the pathobiological characteristics of the most common nodal PTCLs by focusing on the contribution given by high-throughput technologies and the identification of potential therapeutic targets proposed by translational studies.

Dissecting the HIV-specific immune response: a systems biology approach

PURPOSE OF REVIEW

Several unique HIV-infected or HIV-resistant cohorts have been studied over the years to try and delineate the correlates of protection. Although several mechanisms have been put forward, studies aiming to integrate the different mechanisms into a comprehensive model are still lacking. Current systems biology approaches emphasize the importance of unifying independent datasets, provide tools that facilitate hypothesis formulation and testing, and direct us toward uncovering novel therapeutic targets by defining molecular networks perturbed during disease. This review will focus on the current findings that utilized systems biology techniques in order to identify correlates of protection from HIV disease progression and resistance to infection in unique cohorts of individuals as well as in nonhuman primate models of SIV infection.

RECENT FINDINGS

Using systems biology technologies and data analysis tools, the studies described herein have found that pathways implicated in survival, cell cycling, inflammation, and oxidative stress work in unison to limit pathology caused by chronic immune activation. This situation favors the survival of effector lymphocytes and limits the dissemination of viral particles in HIV elite controllers, exposed-uninfected individuals, and natural hosts of SIV infection.

SUMMARY

Systems and computational biology tools have clearly expanded our understanding of HIV pathogenesis by unifying independent observations and by giving us novel molecular targets to pursue. These molecular signatures have the potential to uncover correlates of protection in HIV disease and, in the era of personalized medicine, to determine predictive signatures of treatment efficacy and/or failure.

Integrated Analysis Reveals hsa-miR-142 as a Representative of a Lymphocyte-Specific Gene Expression and Methylation Signature

Gene expression profiling has provided insights into different cancer types and revealed tissue-specific expression signatures. Alterations in microRNA expression contribute to the pathogenesis of many types of human diseases. Few studies have integrated all levels of gene expression, miRNA and methylation to uncover correlations between these data types. We performed an integrated profiling to discover instances of miRNAs associated with a gene expression and DNA methylation signature across multiple cancer types. Using data from The Cancer Genome Atlas (TCGA), we revealed a concordant gene expression and methylation signature associated with the microRNA hsa-miR-142 across the same samples. In all cancer types examined, we found a signature of co-expression of a gene set R and methylated sites M, which correlate positively (M+) or negatively (M−) with the expression of hsa-miR-142. The set R consistently contains many genes, such as TRAF3IP3, NCKAP1L, CD53, LAPTM5, PTPRC, EVI2B, DOCK2, LCP2, CYBB and FYB. The signature is preserved across glioblastoma, ovarian, breast, colon, kidney, lung, uterine and rectum cancer. There is 28% overlap of methylation sites in M between glioblastoma (GBM) and ovarian cancer. There is 60% overlap of genes in R between GBM and ovarian (P = 1.3e⁻¹¹). Most of the genes in R are known to be expressed in lymphocytes and haematopoietic stem cells, while M reflects membrane proteins involved in cell-cell adhesion functions. We speculate that the hsa-miR-142 associated signature may signal haematopoietic-specific processes and an accumulation of methylation events triggering a progressive loss of cell-cell adhesion. We also observed that GBM samples belonging to the proneural subtype tend to have underexpressed hsa-miR-142 and R genes, hypomethylated M+ and hypermethylated M−, while the mesenchymal samples have the opposite profile.

HMGA1 drives stem cell, inflammatory pathway, and cell cycle progression genes during lymphoid tumorigenesis

BACKGROUND

Although the high mobility group A1 (HMGA1) gene is widely overexpressed in diverse cancers and portends a poor prognosis in some tumors, the molecular mechanisms that mediate its role in transformation have remained elusive. HMGA1 functions as a potent oncogene in cultured cells and induces aggressive lymphoid tumors in transgenic mice. Because HMGA1 chromatin remodeling proteins regulate transcription, HMGA1 is thought to drive malignant transformation by modulating expression of specific genes. Genome-wide studies to define HMGA1 transcriptional networks during tumorigenesis, however, are lacking. To define the HMGA1 transcriptome, we analyzed gene expression profiles in lymphoid cells from HMGA1a transgenic mice at different stages in tumorigenesis.

RESULTS

RNA from lymphoid samples at 2 months (before tumors develop) and 12 months (after tumors are well-established) was screened for differential expression of > 20,000 unique genes by microarray analysis (Affymetrix) using a parametric and nonparametric approach. Differential expression was confirmed by quantitative RT-PCR in a subset of genes. Differentially expressed genes were analyzed for cellular pathways and functions using Ingenuity Pathway Analysis. Early in tumorigenesis, HMGA1 induced inflammatory pathways with NFkappaB identified as a major node. In established tumors, HMGA1 induced pathways involved in cell cycle progression, cell-mediated immune response, and cancer. At both stages in tumorigenesis, HMGA1 induced pathways involved in cellular development, hematopoiesis, and hematologic development. Gene set enrichment analysis showed that stem cell and immature T cell genes are enriched in the established tumors. To determine if these results are relevant to human tumors, we knocked-down HMGA1 in human T-cell leukemia cells and identified a subset of genes dysregulated in both the transgenic and human lymphoid tumors.

CONCLUSIONS

We found that HMGA1 induces inflammatory pathways early in lymphoid tumorigenesis and pathways involved in stem cells, cell cycle progression, and cancer in established tumors. HMGA1 also dyregulates genes and pathways involved in stem cells, cellular development and hematopoiesis at both early and late stages of tumorigenesis. These results provide insight into HMGA1 function during tumor development and point to cellular pathways that could serve as therapeutic targets in lymphoid and other human cancers with aberrant HMGA1 expression.

Cell surface markers in HTLV-1 pathogenesis

The phenotype of HTLV-1-transformed CD4⁺ T lymphocytes largely depends on defined viral effector molecules such as the viral oncoprotein Tax. In this review, we exemplify the expression pattern of characteristic lineage markers, costimulatory receptors and ligands of the tumor necrosis factor superfamily, cytokine receptors, and adhesion molecules on HTLV-1-transformed cells. These molecules may provide survival signals for the transformed cells. Expression of characteristic surface markers might therefore contribute to persistence of HTLV-1-transformed lymphocytes and to the development of HTLV-1-associated disease.

Peripheral T-cell lymphoma classification: the matter of cellular derivation

Peripheral T-cell lymphomas (PTCLs) represent approximately 12% of all non-Hodgkin's lymphomas in Western countries. They are quite heterogeneous as far as morphology and phenotype are concerned. Furthermore, until now, PTCLs could not be referred to specific normal counterparts, in contrast to B-cell-derived non-Hodgkin's lymphomas. In particular, in the last edition of the WHO classification of Tumors of the Hematopoietic and Lymphoid Tissues, for the majority of nodal PTCLs (including the not otherwise specified type and anaplastic large-cell lymphoma), the postulated cell of origin remained undefined. However, in the last few years, high-throughput genomic techniques, especially gene-expression profiling, have allowed us to better define the relationship between some entities and the different T-cell subpopulations. Consequently, it has become possible to clearly define, for example, the association between angioimmunoblastic T-cell lymphoma and T-follicular helper cells. In addition, within PTCLs/not otherwise specified, different subgroups were identified based on their similarity to different cellular counterparts, including T-helper, T-cytotoxic and T-follicular helper cells. In this article, based on their own experience as well as up-to-date literature, the authors revise the concept of PTCL classification by specially focusing on their cellular counterparts and discuss the possible clinical relevance of such an approach.

An integrated ChIP-seq analysis platform with customizable workflows

BACKGROUND

Chromatin immunoprecipitation followed by next generation sequencing (ChIP-seq), enables unbiased and genome-wide mapping of protein-DNA interactions and epigenetic marks. The first step in ChIP-seq data analysis involves the identification of peaks (i.e., genomic locations with high density of mapped sequence reads). The next step consists of interpreting the biological meaning of the peaks through their association with known genes, pathways, regulatory elements, and integration with other experiments. Although several programs have been published for the analysis of ChIP-seq data, they often focus on the peak detection step and are usually not well suited for thorough, integrative analysis of the detected peaks.

RESULTS

To address the peak interpretation challenge, we have developed ChIPseeqer, an integrative, comprehensive, fast and user-friendly computational framework for in-depth analysis of ChIP-seq datasets. The novelty of our approach is the capability to combine several computational tools in order to create easily customized workflows that can be adapted to the user's needs and objectives. In this paper, we describe the main components of the ChIPseeqer framework, and also demonstrate the utility and diversity of the analyses offered, by analyzing a published ChIP-seq dataset.

CONCLUSIONS

ChIPseeqer facilitates ChIP-seq data analysis by offering a flexible and powerful set of computational tools that can be used in combination with one another. The framework is freely available as a user-friendly GUI application, but all programs are also executable from the command line, thus providing flexibility and automatability for advanced users.

Densely interconnected transcriptional circuits control cell states in human hematopoiesis

Though many individual transcription factors are known to regulate hematopoietic differentiation, major aspects of the global architecture of hematopoiesis remain unknown. Here, we profiled gene expression in 38 distinct purified populations of human hematopoietic cells and used probabilistic models of gene expression and analysis of cis-elements in gene promoters to decipher the general organization of their regulatory circuitry. We identified modules of highly coexpressed genes, some of which are restricted to a single lineage but most of which are expressed at variable levels across multiple lineages. We found densely interconnected cis-regulatory circuits and a large number of transcription factors that are differentially expressed across hematopoietic states. These findings suggest a more complex regulatory system for hematopoiesis than previously assumed.

Differential genomic targeting of the transcription factor TAL1 in alternate haematopoietic lineages

Expression of the basic helix-loop-helix transcription factor TAL1/SCL is required for erythrocyte differentiation; aberrant expression in lymphoid cells leads to oncogenic transformation. Here, global analysis of TAL1 binding in erythroid and malignant T cells identifies cell type specific functional interaction with the transcription factors RUNX and ETS1.

TAL1/SCL is a master regulator of haematopoiesis whose expression promotes opposite outcomes depending on the cell type: differentiation in the erythroid lineage or oncogenesis in the T-cell lineage. Here, we used a combination of ChIP sequencing and gene expression profiling to compare the function of TAL1 in normal erythroid and leukaemic T cells. Analysis of the genome-wide binding properties of TAL1 in these two haematopoietic lineages revealed new insight into the mechanism by which transcription factors select their binding sites in alternate lineages. Our study shows limited overlap in the TAL1-binding profile between the two cell types with an unexpected preference for ETS and RUNX motifs adjacent to E-boxes in the T-cell lineage. Furthermore, we show that TAL1 interacts with RUNX1 and ETS1, and that these transcription factors are critically required for TAL1 binding to genes that modulate T-cell differentiation. Thus, our findings highlight a critical role of the cellular environment in modulating transcription factor binding, and provide insight into the mechanism by which TAL1 inhibits differentiation leading to oncogenesis in the T-cell lineage.

TLX1 and NOTCH coregulate transcription in T cell acute lymphoblastic leukemia cells

BACKGROUND

The homeobox gene TLX1 (for T-cell leukemia homeobox 1, previously known as HOX11) is inappropriately expressed in a major subgroup of T cell acute lymphoblastic leukemia (T-ALL) where it is strongly associated with activating NOTCH1 mutations. Despite the recognition that these genetic lesions cooperate in leukemogenesis, there have been no mechanistic studies addressing how TLX1 and NOTCH1 functionally interact to promote the leukemic phenotype.

RESULTS

Global gene expression profiling after downregulation of TLX1 and inhibition of the NOTCH pathway in ALL-SIL cells revealed that TLX1 synergistically regulated more than 60% of the NOTCH-responsive genes. Structure-function analysis demonstrated that TLX1 binding to Groucho-related TLE corepressors was necessary for maximal transcriptional regulation of the NOTCH-responsive genes tested, implicating TLX1 modulation of the NOTCH-TLE regulatory network. Comparison of the dataset to publicly available biological databases indicated that the TLX1/NOTCH-coregulated genes are frequently targeted by MYC. Gain- and loss-of-function experiments confirmed that MYC was an essential mediator of TLX1/NOTCH transcriptional output and growth promotion in ALL-SIL cells, with TLX1 contributing to the NOTCH-MYC regulatory axis by posttranscriptional enhancement of MYC protein levels. Functional classification of the TLX1/NOTCH-coregulated targets also showed enrichment for genes associated with other human cancers as well as those involved in developmental processes. In particular, we found that TLX1, NOTCH and MYC coregulate CD1B and RAG1, characteristic markers of early cortical thymocytes, and that concerted downregulation of the TLX1 and NOTCH pathways resulted in their irreversible repression.

CONCLUSIONS

We found that TLX1 and NOTCH synergistically regulate transcription in T-ALL, at least in part via the sharing of a TLE corepressor and by augmenting expression of MYC. We conclude that the TLX1/NOTCH/MYC network is a central determinant promoting the growth and survival of TLX1+ T-ALL cells. In addition, the TLX1/NOTCH/MYC transcriptional network coregulates genes involved in T cell development, such as CD1 and RAG family members, and therefore may prescribe the early cortical stage of differentiation arrest characteristic of the TLX1 subgroup of T-ALL.

Integrating genomic signatures for immunologic discovery

Understanding heterogeneity in adaptive immune responses is essential to dissect pathways of memory B and T cell differentiation and to define correlates of protective immunity. Traditionally, immunologists have deconvoluted this heterogeneity with flow cytometry--with combinations of markers to define signatures that represent specific lineages, differentiation states, and functions. Genome-scale technologies have become widely available and provide the ability to define expression signatures--sets of genes--that represent discrete biological properties of cell populations. Because genomic signatures can serve as surrogates of a phenotype, function, or cell state, they can integrate phenotypic information between experiments, cell types, and species. Here, we discuss how integration of well-defined expression signatures across experimental conditions together with functional analysis of their component genes could provide new opportunities to dissect the complexity of the adaptive immune response and map the immune response to vaccines and pathogens.

Gene coexpression analysis in single cells indicates lymphomyeloid copriming in short-term hematopoietic stem cells and multipotent progenitors

Progressive restriction to a differentiation pathway results from both activation and silencing of particular gene expression programs. To identify the coexpression and the expression levels of regulatory genes during hematopoietic stem cell (HSC) differentiation toward the T cell branch, we applied a new single-cell RT-PCR technique to analyze the simultaneous expression of 13 genes in 9 functionally purified populations from the bone marrow and the thymus. We report in this paper that Lin(-)Sca1(+)ckit(+) HSCs display, at the single-cell level, a homogeneous and high transcriptional activity as do early thymic progenitors. Moreover, the coexpression of lymphoid and myeloid genes is an early event detected in approximately 30% of short-term HSC and most multipotent progenitors, suggesting novel sources for the generation of early thymic progenitors, common lymphoid progenitors (CLPs), and common myeloid progenitors. Loss of multipotency in Lin(-)Sca1(+)ckit(+) cells directed to the lymphoid branch is characterized by Lmo2 and Gata2 gene expression downregulation. Indeed, highest levels of Gata2 expression are detected only in long-term and short-term HSC populations. Complete shutdown of Pu1 gene expression in all triple-negative (TN)3 stage thymic pre-T cells is indicative of total T cell commitment. Interestingly, this is also observed in 30% of TN2 cells and 25% of CLP in the bone marrow, suggesting a possible initiation of T cell engagement in TN2 and CLP. Also, our strategy highlights similar gene patterns among HSCs and intrathymic progenitors, proposing, therefore, that identical activation signals are maintained until further maturation and generation of CD4 and CD8 coreceptors bearing thymocytes.

A whole genome transcriptional analysis of the early immune response induced by live attenuated and inactivated influenza vaccines in young children

The protective mechanisms of influenza vaccines in young children are not completely understood. A phase 2 clinical study was conducted in 85 children 12-35 months of age to describe and compare the immune responses to live attenuated influenza vaccine (LAIV) with trivalent inactivated influenza vaccine (TIV). To better understand the biology of vaccine effects, oligonucleotide microarrays were employed to measure the genome-wide changes in transcript profiles in whole blood at approximately 7 days after 1 dose of LAIV or TIV. Of the total 265 differentially expressed genes identified in this study, 6 clusters of genes were identified to be tightly coexpressed, many of which are likely modulated by cytokines including type 1 interferons (IFNs) and granulocyte-macrophage colony-stimulating factor. Additional functional analyses revealed that the type 1 IFN pathway and cell cycle regulation-related genes are enriched in the 6 coexpressed gene sets. Promoter characterization of these coexpressed genes also supported this conclusion. Moreover, it is suggested that the IFN-stimulated response element is likely to be a potential bidirectional promoter, and the CCAAT/enhancer-binding protein might cooperate with the E2F transcription factor family in the regulation of the cell cycle in the early immune response induced by the influenza vaccine. Overall, our study clearly indicates that the expression profile changes induced by LAIV are significantly different from those induced by TIV. These results suggest that the pattern of overexpression of type 1 IFN-stimulated genes can potentially be used as a biomarker to identify the early vaccination response of LAIV and may also explain, to a certain extent, previous clinical study observations of LAIV-induced protection against influenza-like illness in the first 2 weeks after administration.

Transcriptional response in the peripheral blood of patients infected with Salmonella enterica serovar Typhi

We used microarrays and transcriptional profiling of peripheral blood to investigate the host response of 29 individuals who contracted typhoid fever in the Mekong Delta region of Vietnam. Samples were taken over a nine month period encompassing acute disease, convalescence, and recovery. We found that typhoid fever induced a distinct and highly reproducible signature in the peripheral blood that changed during treatment and convalescence, returning in the majority of cases to the "normal" profile as measured in healthy uninfected controls. Unexpectedly, there was a strong, distinct signature of convalescence present at day 9 after infection that remained virtually unchanged one month after acute infection and in some cases persisted as long as nine months despite a complete clinical recovery in all patients. Patients who retain the convalescent signature may be genetically or temporarily incapable of developing an effective immune response and may be more susceptible to reinfection, relapse, or the establishment of a carrier state.

Human late memory CD8+ T cells have a distinct cytokine signature characterized by CC chemokine production without IL-2 production

Late memory T cell skewing is observed in the setting of immune recovery after cord blood transplantation, and may be associated with inferior control of viral reactivation and cancers. Therefore, we sought to understand how late memory cells differ functionally from earlier stage memory T cells, and whether surface phenotypes associated with differentiation stages were predictably associated with functional signatures. Higher order cytokine flow cytometry allows characterization of human T cells based on complex phenotypic markers and their differential capacity to simultaneously secrete effector proteins, including cytokines and chemokines. We used 8-color, 10-parameter cytokine flow cytometry to characterize the functional activation of human late memory CD8(+) T cells defined by CD45RA and CD27 expression (CD27(-)CD45RA(+)). We assessed the 15 possible functional signatures of cells defined by production of IL-2, IFN-gamma, TNF-alpha, and MIP-1beta alone or in combination, following activation with Ags stimulating bypassing surface proteins (PMA:ionomycin) or through the TCR (e.g., viral Ags). Late memory CD8(+) T cells produced abundant amounts of CC chemokines (MIP-1beta, MIP-1alpha, and RANTES) but not IL-2. IL-2/IFN-gamma coproduction, characteristic of protective immune responses to viral infections, was absent in late memory CD8(+) T cells. These data demonstrate that functional cytokine signatures are predictably associated with CD8(+) maturation stages, and that the polarization of late memory CD8(+) T cells toward CC chemokine production and away from IL-2 production suggests a unique functional role for this subset.

Meta-analysis of gene expression in human pancreatic islets after in vitro expansion

Pancreatic islet transplantation as a potential cure for type 1 diabetes (T1D) cannot be scaled up due to a scarcity of human pancreas donors. In vitro expansion of beta-cells from mature human pancreatic islets provides an alternative source of insulin-producing cells. The exact nature of the expanded cells produced by diverse expansion protocols and their potential for differentiation into functional beta-cells remain elusive. We performed a large-scale meta-analysis of gene expression in human pancreatic islet cells, which were processed using three different previously described protocols for expansion and for which redifferentiation was attempted. All three expansion protocols induced dramatic changes in the expression profiles of pancreatic islets; many of these changes are shared among the three protocols. Attempts at redifferentiation of expanded cells induce a limited number of gene expression changes. Nevertheless, these fail to restore a pancreatic islet-like gene expression pattern. Comparison with a collection of public microarray datasets confirmed that expanded cells are highly comparable to mesenchymal stem cells. Genes induced in expanded cells are also enriched for targets of transcription factors important for pluripotency induction. The present data increase our understanding of the active pathways in expanded and redifferentiated islets. Knowledge of the mesenchymal stem cell potential may help development of drug therapeutics to restore beta-cell mass in T1D patients.