Transcriptome dynamics of CD4+ T cells during malaria maps gradual transit from effector to memory.

The dynamics of CD4+ T cell memory development remain to be examined at genome scale. In malaria-endemic regions, antimalarial chemoprevention protects long after its cessation and associates with effects on CD4+ T cells. We applied single-cell RNA sequencing and computational modelling to track memory development during Plasmodium infection and treatment. In the absence of central memory precursors, two trajectories developed as T helper 1 (TH1) and follicular helper T (TFH) transcriptomes contracted and partially coalesced over three weeks. Progeny of single clones populated TH1 and TFH trajectories, and fate-mapping suggested that there was minimal lineage plasticity. Relationships between TFH and central memory were revealed, with antimalarials modulating these responses and boosting TH1 recall. Finally, single-cell epigenomics confirmed that heterogeneity among effectors was partially reset in memory. Thus, the effector-to-memory transition in CD4+ T cells is gradual during malaria and is modulated by antiparasitic drugs. Graphical user interfaces are presented for examining gene-expression dynamics and gene-gene correlations ( http://haquelab.mdhs.unimelb.edu.au/cd4_memory/ ).

Single-Cell Survey of Human Lymphatics Unveils Marked Endothelial Cell Heterogeneity and Mechanisms of Homing for Neutrophils.

Lymphatic vessels form a critical component in the regulation of human health and disease. While their functional significance is increasingly being recognized, the comprehensive heterogeneity of lymphatics remains uncharacterized. Here, we report the profiling of 33,000 lymphatic endothelial cells (LECs) in human lymph nodes (LNs) by single-cell RNA sequencing. Unbiased clustering revealed six major types of human LECs. LECs lining the subcapsular sinus (SCS) of LNs abundantly expressed neutrophil chemoattractants, whereas LECs lining the medullary sinus (MS) expressed a C-type lectin CD209. Binding of a carbohydrate Lewis X (CD15) to CD209 mediated neutrophil binding to the MS. The neutrophil-selective homing by MS LECs may retain neutrophils in the LN medulla and allow lymph-borne pathogens to clear, preventing their spread through LNs in humans. Our study provides a comprehensive characterization of LEC heterogeneity and unveils a previously undefined role for medullary LECs in human immunity.

Dissecting the polygenic basis of atherosclerosis via disease-associated cell state signatures.

Coronary artery disease (CAD) is a pandemic disease where up to half of the risk is explained by genetic factors. Advanced insights into the genetic basis of CAD require deeper understanding of the contributions of different cell types, molecular pathways, and genes to disease heritability. Here, we investigate the biological diversity of atherosclerosis-associated cell states and interrogate their contribution to the genetic risk of CAD by using single-cell and bulk RNA sequencing (RNA-seq) of mouse and human lesions. We identified 12 disease-associated cell states that we characterized further by gene set functional profiling, ligand-receptor prediction, and transcription factor inference. Importantly, Vcam1+ smooth muscle cell state genes contributed most to SNP-based heritability of CAD. In line with this, genetic variants near smooth muscle cell state genes and regulatory elements explained the largest fraction of CAD-risk variance between individuals. Using this information for variant prioritization, we derived a hybrid polygenic risk score (PRS) that demonstrated improved performance over a classical PRS. Our results provide insights into the biological mechanisms associated with CAD risk, which could make a promising contribution to precision medicine and tailored therapeutic interventions in the future.

Single-cell characterisation of tissue homing CD4 + and CD8 + T cell clones in immune-mediated refractory arthritis.

BACKGROUND: Immune-mediated arthritis is a group of autoinflammatory diseases, where the patient's own immune system attacks and destroys synovial joints. Sustained remission is not always achieved with available immunosuppressive treatments, warranting more detailed studies of T cell responses that perpetuate synovial inflammation in treatment-refractory patients. METHODS: In this study, we investigated CD4 + and CD8 + T lymphocytes from the synovial tissue and peripheral blood of patients with treatment-resistant immune-mediated arthritis using paired single-cell RNA and TCR-sequencing. To gain insights into the trafficking of clonal families, we compared the phenotypes of clones with the exact same TCRß amino acid sequence between the two tissues. RESULTS: Our results show that both CD4 + and CD8 + T cells display a more activated and inflamed phenotype in the synovial tissue compared to peripheral blood both at the population level and within individual T cell families. Furthermore, we found that both cell subtypes exhibited clonal expansion in the synovial tissue. CONCLUSIONS: Our findings suggest that the local environment in the synovium drives the proliferation of activated cytotoxic T cells, and both CD4 + and CD8 + T cells may contribute to tissue destruction and disease pathogenesis.

Transinteractome analysis reveals distinct niche requirements for isotype-based plasma cell subsets in the bone marrow.

Bone marrow (BM) long-lived plasma cells (PCs) are essential for long-term protection against infection, and their persistence within this organ relies on interactions with Cxcl12-expressing stromal cells that are still not clearly identified. Here, using single cell RNAseq and in silico transinteractome analyses, we identified Leptin receptor positive (LepR+ ) mesenchymal cells as the stromal cell subset most likely to interact with PCs within the BM. Moreover, we demonstrated that depending on the isotype they express, PCs may use different sets of integrins and adhesion molecules to interact with these stromal cells. Altogether, our results constitute an unprecedented characterization of PC subset stromal niches and open new avenues for the specific targeting of BM PCs based on their isotype.

Single-Cell Epigenomics and Functional Fine-Mapping of Atherosclerosis GWAS Loci.

[Figure: see text].

Single-cell characterization of dog allergen-specific T cells reveals TH2 heterogeneity in allergic individuals.

BACKGROUND: Allergen-specific type 2 CD4+ TH2 cells are critically involved in the pathogenesis of IgE-mediated allergic diseases. However, the heterogeneity of the TH2 response has only recently been appreciated. OBJECTIVE: We sought to characterize at the single-cell level the ex vivo phenotype, transcriptomic profile, and T-cell receptor (TCR) repertoire of circulating CD4+ T cells specific to the major dog allergens Can f 1, Can f 4, and Can f 5 in subjects with and without dog allergy. METHODS: Dog allergen-specific memory CD4+ T cells were detected ex vivo by flow cytometry using a CD154-based enrichment assay and single-cell sorted for targeted gene expression analysis and TCR sequencing. RESULTS: Dog allergen-specific T-cell responses in allergic subjects were dominantly of TH2 type. TH2 cells could be phenotypically further divided into 3 subsets, which consisted of TH2-like (CCR6-CXCR3-CRTH2-), TH2 (CCR6-CXCR3-CRTH2+CD161-), and TH2A (CCR6-CXCR3-CRTH2+CD161+CD27-) cells. All these subsets were nonexistent within the allergen-specific T-cell repertoire of healthy subjects. Single-cell transcriptomic profiling confirmed the TH2-biased signature in allergen-specific T cells from allergic subjects and revealed a TH1/TH17 signature in nonallergic subjects. TCR repertoire analyses showed that dog allergen-specific T cells were diverse and allergic subjects demonstrated less clonality compared to nonallergic donors. Finally, TCR and transcriptomic analyses revealed a close relationship between TH2-like, TH2, and TH2A cells, with the last ones representing the most terminally differentiated and highly polarized subtype. CONCLUSIONS: Our study demonstrates heterogeneity within allergen-specific TH2 cells at the single-cell level. The results may be utilized for improving immune monitoring after allergen immunotherapy and for designing targeted immunomodulatory approaches.

Cell type markers indicate distinct contributions of decidual stromal cells and natural killer cells in preeclampsia.

In brief: Preeclampsia is a common serious disorder that can occur during pregnancy. This study uses integrative analysis of preeclampsia transcriptomes and single-cell transcriptomes to predict cell type-specific contributions to preeclampsia. Abstract: Preeclampsia is a devastating pregnancy disorder and a major cause of maternal and perinatal mortality. By combining previous transcriptomic results on preeclampsia with single-cell sequencing data, we here predict distinct and partly unanticipated contributions of decidual stromal cells and uterine natural killer cells in early- and late-onset preeclampsia.

Global chromatin state analysis reveals lineage-specific enhancers during the initiation of human T helper 1 and T helper 2 cell polarization.

Naive CD4⁺ T cells can differentiate into specific helper and regulatory T cell lineages in order to combat infection and disease. The correct response to cytokines and a controlled balance of these populations is critical for the immune system and the avoidance of autoimmune disorders. To investigate how early cell-fate commitment is regulated, we generated the first human genome-wide maps of histone modifications that reveal enhancer elements after 72 hr of in vitro polarization toward T helper 1 (Th1) and T helper 2 (Th2) cell lineages. Our analysis indicated that even at this very early time point, cell-specific gene regulation and enhancers were at work directing lineage commitment. Further examination of lineage-specific enhancers identified transcription factors (TFs) with known and unknown T cell roles as putative drivers of lineage-specific gene expression. Lastly, an integrative analysis of immunopathogenic-associated SNPs suggests a role for distal regulatory elements in disease etiology.

Microanatomy of the Human Atherosclerotic Plaque by Single-Cell Transcriptomics.

RATIONALE: Atherosclerotic lesions are known for their cellular heterogeneity, yet the molecular complexity within the cells of human plaques has not been fully assessed. OBJECTIVE: Using single-cell transcriptomics and chromatin accessibility, we gained a better understanding of the pathophysiology underlying human atherosclerosis. METHODS AND RESULTS: We performed single-cell RNA and single-cell ATAC sequencing on human carotid atherosclerotic plaques to define the cells at play and determine their transcriptomic and epigenomic characteristics. We identified 14 distinct cell populations including endothelial cells, smooth muscle cells, mast cells, B cells, myeloid cells, and T cells and identified multiple cellular activation states and suggested cellular interconversions. Within the endothelial cell population, we defined subsets with angiogenic capacity plus clear signs of endothelial to mesenchymal transition. CD4+ and CD8+ T cells showed activation-based subclasses, each with a gradual decline from a cytotoxic to a more quiescent phenotype. Myeloid cells included 2 populations of proinflammatory macrophages showing IL (interleukin) 1B or TNF (tumor necrosis factor) expression as well as a foam cell-like population expressing TREM2 (triggering receptor expressed on myeloid cells 2) and displaying a fibrosis-promoting phenotype. ATACseq data identified specific transcription factors associated with the myeloid subpopulation and T cell cytokine profiles underlying mutual activation between both cell types. Finally, cardiovascular disease susceptibility genes identified using public genome-wide association studies data were particularly enriched in lesional macrophages, endothelial, and smooth muscle cells. CONCLUSIONS: This study provides a transcriptome-based cellular landscape of human atherosclerotic plaques and highlights cellular plasticity and intercellular communication at the site of disease. This detailed definition of cell communities at play in atherosclerosis will facilitate cell-based mapping of novel interventional targets with direct functional relevance for the treatment of human disease.

T cell fate and clonality inference from single-cell transcriptomes.

We developed TraCeR, a computational method to reconstruct full-length, paired T cell receptor (TCR) sequences from T lymphocyte single-cell RNA sequence data. TraCeR links T cell specificity with functional response by revealing clonal relationships between cells alongside their transcriptional profiles. We found that T cell clonotypes in a mouse Salmonella infection model span early activated CD4(+) T cells as well as mature effector and memory cells.

Single-cell technologies are revolutionizing the approach to rare cells.

In the last lustrum single-cell techniques such as single-cell quantitative PCR, RNA and DNA sequencing, and the state-of-the-art cytometry by time of flight (CyTOF) mass cytometer have allowed a detailed analysis of the sub-composition of different organs from the bone marrow hematopoietic compartment to the brain. These fine-grained analyses have highlighted the great heterogeneity within each cell compartment revealing previously unknown subpopulations of cells. In this review, we analyze how this fast technological evolution has improved our understanding of the biological processes with a particular focus on rare cells of the immune system.

Gene Regulatory Network Analysis of Decidual Stromal Cells and Natural Killer Cells.

Human reproductive success relies on the proper differentiation of the uterine endometrium to facilitate implantation, formation of the placenta, and pregnancy. This process involves two critical types of decidual uterine cells: endometrial/decidual stromal cells (dS) and uterine/decidual natural killer (dNK) cells. To better understand the transcription factors governing the in vivo functions of these cells, we analyzed single-cell transcriptomics data from first-trimester terminations of pregnancy, and for the first time conducted gene regulatory network analysis of dS and dNK cell subpopulations. Our analysis revealed stromal cell populations that corresponded to previously described in vitro decidualized cells and senescent decidual cells. We discovered new decidualization driving transcription factors of stromal cells for early pregnancy, including DDIT3 and BRF2, which regulate oxidative stress protection. For dNK cells, we identified transcription factors involved in the immunotolerant (dNK1) subpopulation, including IRX3 and RELB, which repress the NFKB pathway. In contrast, for the less immunotolerant (dNK3) population we predicted TBX21 (T-bet) and IRF2-mediated upregulation of the interferon pathway. To determine the clinical relevance of our findings, we tested the overrepresentation of the predicted transcription factors target genes among cell type-specific regulated genes from pregnancy disorders, such as recurrent pregnancy loss and preeclampsia. We observed that the predicted decidualized stromal and dNK1-specific transcription factor target genes were enriched with the genes downregulated in pregnancy disorders, whereas the predicted dNK3-specific targets were enriched with genes upregulated in pregnancy disorders. Our findings emphasize the importance of stress tolerance pathways in stromal cell decidualization and immunotolerance promoting regulators in dNK differentiation.

Single-cell transcriptome analysis of the early immune response in the lymph nodes of Borrelia burgdorferi-infected mice.

Lyme borreliosis is a disease caused by Borrelia burgdorferi sensu lato bacteria. Borrelia burgdorferi is known to induce prolonged extrafollicular immune responses and abnormal germinal centre formation. The infection fails to generate a neutralizing type of immunity, eventually establishing a persistent infection. Here, we performed single-cell RNA sequencing to characterize the immune landscape of lymph node lymphocytes during the early Borrelia burgdorferi infection in a murine model. Our results indicate key features of an extrafollicular immune response four days after Borrelia burgdorferi infection, including notable B cell proliferation, immunoglobulin class switching to IgG3 and IgG2b isotypes, plasmablast differentiation, and the presence of extrafollicular B cells identified through immunohistochemistry. Additionally, we found infection-derived upregulation of suppressor of cytokine signalling genes Socs1 and Socs3, along with downregulation of genes associated with MHC II antigen presentation in B cells. Our results support the central role of B cells in the immune response of a Borrelia burgdorferi infection, and provide cues of mechanisms behind the determination between extrafollicular and germinal centre responses during Borrelia burgdorferi infection.

PIM kinases regulate early human Th17 cell differentiation.

The serine/threonine-specific Moloney murine leukemia virus (PIM) kinase family (i.e., PIM1, PIM2, and PIM3) has been extensively studied in tumorigenesis. PIM kinases are downstream of several cytokine signaling pathways that drive immune-mediated diseases. Uncontrolled T helper 17 (Th17) cell activation has been associated with the pathogenesis of autoimmunity. However, the detailed molecular function of PIMs in human Th17 cell regulation has yet to be studied. In the present study, we comprehensively investigated how the three PIMs simultaneously alter transcriptional gene regulation during early human Th17 cell differentiation. By combining PIM triple knockdown with bulk and scRNA-seq approaches, we found that PIM deficiency promotes the early expression of key Th17-related genes while suppressing Th1-lineage genes. Further, PIMs modulate Th cell signaling, potentially via STAT1 and STAT3. Overall, our study highlights the inhibitory role of PIMs in human Th17 cell differentiation, thereby suggesting their association with autoimmune phenotypes.

Single-cell characterization of anti-LAG-3 and anti-PD-1 combination treatment in patients with melanoma.

BackgroundRelatlimab plus nivolumab (anti-lymphocyte-activation gene 3 plus anti-programmed death 1 [anti-LAG-3+anti-PD-1]) has been approved by the FDA as a first-line therapy for stage III/IV melanoma, but its detailed effect on the immune system is unknown.MethodsWe evaluated blood samples from 40 immunotherapy-naive or prior immunotherapy-refractory patients with metastatic melanoma treated with anti-LAG-3+anti-PD-1 in a phase I trial using single-cell RNA and T cell receptor sequencing (scRNA+TCRαß-Seq) combined with other multiomics profiling.ResultsThe highest LAG3 expression was noted in NK cells, Tregs, and CD8+ T cells, and these cell populations underwent the most significant changes during the treatment. Adaptive NK cells were enriched in responders and underwent profound transcriptomic changes during the therapy, resulting in an active phenotype. LAG3+ Tregs expanded, but based on the transcriptome profile, became metabolically silent during the treatment. Last, higher baseline TCR clonality was observed in responding patients, and their expanding CD8+ T cell clones gained a more cytotoxic and NK-like phenotype.ConclusionAnti-LAG-3+anti-PD-1 therapy has profound effects on NK cells and Tregs in addition to CD8+ T cells.Trial registrationClinicalTrials.gov (NCT01968109)FundingCancer Foundation Finland, Sigrid Juselius Foundation, Signe and Ane Gyllenberg Foundation, Relander Foundation, State funding for university-level health research in Finland, a Helsinki Institute of Life Sciences Fellow grant, Academy of Finland (grant numbers 314442, 311081, 335432, and 335436), and an investigator-initiated research grant from BMS.

An integrative computational systems biology approach identifies differentially regulated dynamic transcriptome signatures which drive the initiation of human T helper cell differentiation.

BACKGROUND: A proper balance between different T helper (Th) cell subsets is necessary for normal functioning of the adaptive immune system. Revealing key genes and pathways driving the differentiation to distinct Th cell lineages provides important insight into underlying molecular mechanisms and new opportunities for modulating the immune response. Previous computational methods to quantify and visualize kinetic differential expression data of three or more lineages to identify reciprocally regulated genes have relied on clustering approaches and regression methods which have time as a factor, but have lacked methods which explicitly model temporal behavior. RESULTS: We studied transcriptional dynamics of human umbilical cord blood T helper cells cultured in absence and presence of cytokines promoting Th1 or Th2 differentiation. To identify genes that exhibit distinct lineage commitment dynamics and are specific for initiating differentiation to different Th cell subsets, we developed a novel computational methodology (LIGAP) allowing integrative analysis and visualization of multiple lineages over whole time-course profiles. Applying LIGAP to time-course data from multiple Th cell lineages, we identified and experimentally validated several differentially regulated Th cell subset specific genes as well as reciprocally regulated genes. Combining differentially regulated transcriptional profiles with transcription factor binding site and pathway information, we identified previously known and new putative transcriptional mechanisms involved in Th cell subset differentiation. All differentially regulated genes among the lineages together with an implementation of LIGAP are provided as an open-source resource. CONCLUSIONS: The LIGAP method is widely applicable to quantify differential time-course dynamics of many types of datasets and generalizes to any number of conditions. It summarizes all the time-course measurements together with the associated uncertainty for visualization and manual assessment purposes. Here we identified novel human Th subset specific transcripts as well as regulatory mechanisms important for the initiation of the Th cell subset differentiation.

Quantitative proteomics analysis of the nuclear fraction of human CD4+ cells in the early phases of IL-4-induced Th2 differentiation.

We used stable isotope labeling with 4-plex iTRAQ (isobaric tags for relative and absolute quantification) reagents and LC-MS/MS to investigate proteomic changes in the nucleus of activated human CD4(+) cells during the early stages of Th2 cell differentiation. The effects of IL-4 stimulation upon activated naïve CD4(+) cells were measured in the nuclear fractions from 6 and 24 h in three biological replicates, each using pooled cord blood samples derived from seven or more individuals. In these analyses, in the order of 800 proteins were detected with two or more peptides and quantified in three biological replicates. In addition to consistent differences observed with the nuclear localization/expression of established human Th2 and Th1 markers, there were changes that suggested the involvement of several proteins either only recently reported or otherwise not known in this context. These included SATB1 and among the novel changes detected and validated an IL-4-induced increase in the level of YB1. This unique data set from human cord blood CD4(+) T cells details an extensive list of protein determinations that compares with and complements previous data determined from the Jurkat cell nucleus.

Evolution and modulation of antigen-specific T cell responses in melanoma patients.

Analyzing antigen-specific T cell responses at scale has been challenging. Here, we analyze three types of T cell receptor (TCR) repertoire data (antigen-specific TCRs, TCR-repertoire, and single-cell RNA + TCRαß-sequencing data) from 515 patients with primary or metastatic melanoma and compare it to 783 healthy controls. Although melanoma-associated antigen (MAA) -specific TCRs are restricted to individuals, they share sequence similarities that allow us to build classifiers for predicting anti-MAA T cells. The frequency of anti-MAA T cells distinguishes melanoma patients from healthy and predicts metastatic recurrence from primary melanoma. Anti-MAA T cells have stem-like properties and frequent interactions with regulatory T cells and tumor cells via Galectin9-TIM3 and PVR-TIGIT -axes, respectively. In the responding patients, the number of expanded anti-MAA clones are higher after the anti-PD1(+anti-CTLA4) therapy and the exhaustion phenotype is rescued. Our systems immunology approach paves the way for understanding antigen-specific responses in human disorders.