Transcriptome and TCR Repertoire Measurements of CXCR3+ T Follicular Helper Cells Within HIV-Infected Human Lymph Nodes.

Follicular-helper T cells (TFH) are an essential arm of the adaptive immune system. Although TFH were first discovered through their ability to contribute to antibody affinity maturation through co-stimulatory interactions with B cells, new light has been shed on their ability to remain a complex and functionally plastic cell type. Due to a lack sample availability, however, many studies have been limited to characterizing TFH in mice or non-canonical tissue types, such as peripheral blood. Such constraints have resulted in a limited, and sometimes contradictory, understanding of this fundamental cell type. One subset of TFH receiving attention in chronic infection are CXCR3-expressing TFH cells (CXCR3+TFH) due to their abnormal accumulation in secondary lymphoid tissues. Their function and clonal relationship with other TFH subsets in lymphoid tissues during infection, however, remains largely unclear. We thus systematically investigated this and other subsets of TFH within untreated HIV-infected human lymph nodes using Mass CyTOF and a combination of RNA and TCR repertoire sequencing. We show an inflation of the CXCR3+TFH compartment during HIV infection that correlates with a lower HIV burden. Deeper analysis into this population revealed a functional shift of CXCR3+TFH away from germinal center TFH (GC-TFH), including the altered expression of several important transcription factors and cytokines. CXCR3+TFH also upregulated cell migration transcriptional programs and were clonally related to peripheral TFH populations. In combination, these data suggest that CXCR3+TFH have a greater tendency to enter circulation than their CXCR3- counterparts, potentially functioning through distinct modalities that may lead to enhanced defense.

Introduction

Follicular-helper T cells (TFH) provide B cells with the necessary costimulatory signals for their affinity maturation of broadly neutralizing antibodies (BNAbs) (1), a critical component of host defense. Despite high levels of circulating IgG in HIV-infected patients (2, 3), BNAbs and natural viral control are uncommon (4, 5). Given the importance of TFH contribution in a successful antibody response, it is possible their inadequacy during HIV infection leads to a lack of effective and timely B cell help.

Surprisingly, sheer numbers of TFH do not appear to be the problem. In fact, TFH exist at high frequencies during untreated HIV infection (6, 7), suggesting a qualitative, rather than quantitative, impotence. In fact, we and others have shown that TFH in HIV+ patients are functionally and phenotypically distinct compared to healthy donors (8, 9), even during effective anti-retroviral therapy (ART). In addition to their altered cell state, HIV-specific TFH appear to be slow-acting, requiring prolonged exposure to antigen likely during a period of uncontrolled viral replication (10). Further complicating matters, TFH are readily infected by HIV. Several studies have implicated their role as an important reservoir for latent infection (11–14) and subsequent viral persistence. This complex interplay between host and virus have made it a great challenge to understand the pitfalls of the immune system during HIV infection.

A high degree of cellular plasticity in the TFH compartment (15–17) has also led to our incomplete understanding of its role in HIV infection. TFH are generally identified by their expression of CXCR5, a chemokine receptor necessary for proper follicular localization in the lymph node (LN) (18, 19). However, several less understood subsets have been characterized based on other cell surface markers. For example, in experimental autoimmune encephalomyelitis (EAE), a mouse disease model for multiple sclerosis, IL17a-expressing TFH have been found to infiltrate the central nervous system and exacerbate inflammation and disease (20). In contrast, antigen-specific CCR6-expressing TFH were observed at high levels in SARS-COV-2 infection, correlating with less severe disease (21). Interestingly, CCR6+TFH in this context, normally described as TH17-like, seldom expressed IL17a. In an opposite functional context, regulatory TFH have also been described as FOXP3+CXCR5+TFH that appear to suppress the germinal center reaction (22). Recently, we described an activated CXCR5- TFH-like population in LNs that accumulates in HIV infection. These cells specifically co-express high levels of ICOS and PD-1, are clonally related to TFH and can promote the production of HIV-specific antibodies from B cells in vitro (9). We speculate that these cells, having a heightened migratory gene program, may infiltrate inflamed LNs and modify the B cell response; however, our understanding of their role in infection and relationship with other CXCR5+ TFH populations remains insufficient.

Another TFH subset that has received a fair amount of attention expresses the bona fide TH1 marker, CXCR3. A recent study depicting the clonal relationships between TFH subsets in the tonsils and blood of healthy humans found a strong clonal connection between CXCR3-expressing TFH (CXCR3+TFH) and circulating TFH (cTFH) (23). In agreement with previous publications (24, 25), the authors posit that the origin of germinal center TFH (GC-TFH) may be CXCR3+TFH that initially encountered antigen in the periphery. Interestingly, other studies have shown CXCR3+TFH to be enriched during chronic SIV infection in non-human primates (26, 27), in the blood during HIV infection (28, 29), and after influenza vaccination (30). Given their inflation in non-steady state conditions, it is possible that studying CXCR3+TFH in these contexts may reveal a clearer picture of their role leading up to and during their B cell interactions. Furthermore, with an altered cell state appearing to be a hallmark of TFH in HIV disease progression (8), it is possible that understanding these alternative TFH subsets will bring light to novel mechanisms of disease and/or host defense. A systematic study of CXCR3+TFH during perturbation, across multiple tissues in humans is thus of critical need to fully understand their function within the adaptive immune system.

Although blood studies in humans and work done with model organisms have laid the foundation for our understanding of TFH, their results have been contradictory at times. For example, Velu et al. observed similar helper function between LN CXCR3+TFH and CXCR3-TFH to induce SIV-specific IgG in vitro (27), while other studies suggested that CXCR3+TFH from human blood provide poor B cell help (31, 32). It is likely that context, such as tissue type and disease state, has a great impact on the role and function of TFH. As such, our grasp remains limited without clear delineation of the relationship between cTFH and their native, lymphoid-resident counterparts in both health and disease. In this study, we aim at narrowing this deficit in knowledge using a combination of Mass CyTOF and TCR repertoire and RNA sequencing on several different T cell populations from paired LNs and peripheral blood samples in untreated HIV+ patients. Importantly, we find that CXCR3+TFH are both inflated, phenotypically distinct, and correlate with lower HIV burden. We thus emphasize our analysis on this population and its relationship to several better characterized T cell populations within the human LN. Our multi-parameter approach revealed enhanced proliferative potential, an upregulation of cell migration pathways, and strong cross-tissue and cross-phenotype clonal relationships within CXCR3+TFH. These data suggest the potential of CXCR3+TFH to expand within the LN, enter circulation and possibly contribute to host defense through alternative processes than their canonical GC-TFH counterparts. Further investigation into the impact of CXCR3+TFH on viral load and spread over time with a larger, longitudinally tracked cohort will be important to unravel specific mechanisms that these may operate through.

Materials and Methods

Study Subjects

Subjects included in this study are from a subset of patients recruited for our previous study (8). Briefly, LN samples from HIV+ donors were excised from palpable cervical LNs for clinical diagnostic workup in Mexico. HD samples were de-identified mesenteric or inguinal lymph nodes from the Cooperative Human Tissue Network (CHTN). Twenty-two HIV+ donor LNs and nine HD LNs were used for Mass CyTOF experiments. Considering sample availability and size, LN and PBMC samples from an independent, but semi-overlapping cohort () were collected from LN and PBMC samples were collected from seven non-ART and one ART-treated HIV+ patients. Sample sizes were not pre-specified but dictated by the availability of the sample. All samples from HIV+ patients were de-identified and were obtained in accordance with the Declaration of Helsinki after obtaining written informed consent of participants and as part of protocol B03-16, which was reviewed and approved by the Research Committee and the Ethics in Research Committee of the National Institute of Respiratory Diseases “Ismael Cosío Villegas”, Mexico City.

Mass CyTOF and Analysis

Mass CyTOF data were collected as described in a previous publication (8), except that the cells were not stimulated prior to staining. Data were collected and normalized as described there. In analysis specific to this paper, cells were subset by drilling down on the gates defined in with FlowJo™ v10.8.0 Software (BD Life Sciences) and exported as individual FCS files. FCS files were then read into R using the “flowCore” package and analyzed using custom scripts. Dimensional reduction and clustering were done using the R package “Seurat”. Cluster annotations and rational are described in . Boxplots and UMAP plots were generating using the R package “ggplot2”. The heatmap in was made with the R package “pheatmap” (33). Statistical analyses were done using the R package “rstatix”. Final plots were compiled in Biorender.

Figure 1

Mass CyTOF reveals an enrichment of non-canonical CXCR3-expressing TFH unique to HIV+ lymph nodes. (A) TFH populations were in silico sorted based on the gating illustrated in . The frequencies of HIV+ and HD LN cells within each population were compared with a t-test. (B) CXCR3+TFH, CXCR3-TFH, CXCR5-TFH and T cells were in silico sorted, visualized with UMAP, and clustered (left). Clusters were annotated based on the heatmap of the markers shared across each experiment (right). (C) Log pVL and the ratio of CXCR3+TFH to CXCR3-TFH frequencies were compared for each donor using a Spearman’s rank correlation coefficient. (D) The frequencies of HIV+ patient cells versus HD cells for each donor within each cluster from B were compared using a Wilcox test. (E) Differential expression HIV+ versus HD derived cells in cluster 5 was calculated for each marker based on a likelihood-ratio test. (F) Pairwise frequencies of in silico sorted TFH within each cluster (HIV+ donors only) was compared using a t-test. (G) The median UMAP coordinates were calculated for each donor in silico sorted populations and projected onto the same coordinates as in (B) Significance: ****p < 0.0001; ***p < 0.001; **p < 0.01; *p < 0.05; ns p ≥ 0.05.

Cell Staining and Sorting by Flow Cytometry

Cell staining and sorting were performed as previously described (9). Briefly, cryopreserved LN cell suspensions were freshly thawed. A BD FACSARIA was used to sort live CD4+ T cells as CD3+CD4+CD8-CD14-CD19-. Then, from the live CD4+ T cell population, naïve (CD45RO-CXCR5-CD57-CCR7+), CXCR3+TFH (CD45RO+CXCR5+CCR7-PD1+CXCR3+), CXCR3-TFH (CD45RO+CXCR5+CCR7-PD1+CXCR3-), GC-TFH (CD45RO+CXCR5+CD57+) and CXCR5-TFH (CD45RO+CXCR5-PD1+ICOS+) populations were sorted directly into lysis buffer and stored at -80°C until further use.

RNA Purification

Total RNA was purified using AllPrep DNA/RNA Mini and RNeasy MiniElute kits (Qiagen) according to the manufacturer’s instructions. Purified RNA quality was evaluated using Agilent RNA 6000 Pico Kit. Samples were then either supplemented with RNase inhibitor (RNase OUT, Thermofisher Scientific) and stored at -80°C or taken directly to reverse transcription.

Bulk TCRβ Sequencing and Analysis

TCRβ repertoire sequencing libraries were prepared and sequenced as described previously (8). Briefly, 30% of purified RNA was used for reverse transcription. Molecular identifiers (MIDs) were added to cDNA templates during reverse transcription using a TRBC-binding primer with 12 random nucleotides and a partial Illumina adapter (). PCR1 was carried out using multiplexed TRBV-binding primers. A second round of PCR was then used to append final Illumina sequencing adapters to the TCRβ junction-containing inserts. Final libraries were sequenced on an Illumina paired end 150x150 configuration at a minimum depth of 10 reads/cell. After sequencing, reads were clustered based on their MID and TCRβ sequence similarity. Consensus sequences were then built to correct PCR and sequencing errors as described (34). TCRβ sequencing information is summarized in . The CDR3 blast module MIGEC (35) was used for TRBV/J and CDR3 annotation. Bhattacharyya coefficient (36) was used to measure TCR repertoire similarity. Circos plots were generated using circlize R package (37).

RNA Sequencing and Analysis

RNA-Seq libraries were prepared using a protocol modified from SMART-seq2 (SSII) (38). Briefly, 2ul of purified total RNA was reverse transcribed using a poly-T primer fused to the ISPCR handle (; RT_dT30VN). Second-strand synthesis was then done using the SSII template-switching oligo (; SSII_TSO). The following program was used: 42°C RT for 90min, 10 cycles of 50°C for 2min, 42°C for 2min and then 70°C for 15min. cDNA amplification was then done using KAPA HIFI and the ISPCR handle (; RT_TSPCR) with the following program: 98°C initial denature for 3min, 16 cycles of 98°C denature for 20sec, 64°C annealing for 30sec, 72°C extension for 6min and 72°C final extension for 5min. The PCR product was then purified using AmpureXP beads (Agencourt), according to the manufacturer’s protocol. Purified PCR product was then diluted to 0.1 – 0.3 ng/ul and tagmented using the Nextera XT kit (Illumina) with a reduced reaction volume. Briefly, 1.25 ul of diluted sample was used in a 5ul total reaction volume and fragmented for 10min at 55°C, and then held at 10°C. The reaction was then neutralized by adding 1.25 ul of NT buffer. Final libraries were then generated from the tagmentation product using Nextera adapters () and the following program: 72°C for 3 min, 95°C denature for 3sec, 12 cycles of 95°C denature for 10sec, 55°C annealing for 30sec, 72°C extension for 1min and 72°C final extension for 5min. Indexed libraries were then purified using Ampure XP beads according to manufacturer’s instructions and quantified with an Agilent High Sensitivity DNA kit. The final libraries were pooled and sequenced on Illumina HiSeq (150x150) with a minimum depth of about 200 reads per cell. Each library was split into two technique replicates and sequenced independently to ensure the reliability of RNA sequencing results.

Sequencing reads were aligned to human reference genome GRCh38.p5 using RSEM (39). Differentially expressed genes were analyzed using DESeq2 (40). GSVA, a non-parametric unsupervised method that quantifies the relative enrichment of selected pathways, was performed using the R package, GSVA (41). Gene Set Enrichment Analysis (GSEA) on selected pathways or gene sets were performed with the R package fgsea (42). Tonsil Tfh signatures (GSE50391, CXCR5highCD45RO+ versus CXCR5- tonsil samples) were created using the GEO2R online tool (https://www.ncbi.nlm.nih.gov/geo/geo2r/). Metascape (43) and NetworkAnalyst (44) was used to determine GO pathway enrichment.

Quantification of HIV Transcripts From Bulk RNA-Seq

HIV transcripts were quantified from bulk RNA-Seq following the stepwise procedures () similar to the method outlined in Chang et al. (45). Briefly, reads that could not be aligned to the human genome were mapped to the HIV genome (NC_001802.1) with the gapped aligner software HISAT2 (46), accounting for splice junction events. Quantification of HIV transcript expression was calculated by normalizing the HIV genome mappable reads to human genome mappable reads and the input cell number. For a more reliable quantification of HIV copies in each library, the correlation between two replicates was first evaluated (), then the HIV copy numbers of both technical replicates were averaged.

Statistical Analysis

A paired sample two-tailed Student’s t-test was used for pairwise comparisons, with a P value less than 0.05 as a cutoff to determine statistical significance. These comparisons include CXCR3+/CXCR3-TFH percentages in , TCR repertoire similarities in and HIV copies in . For gene set enrichment analysis, CAMERA (47) (Correlation Adjusted MEan RAnk gene set test) was used to calculate enrichment significance. Pearson correlation was used to determine the degree of association.

Figure 4

TCR repertoire analysis demonstrates the clonal relationship of CXCR3+TFH with the CXCR5-TFH outside GC and cTFH in the blood. (A) Average BHA similarity among different T cell subsets. (B) Dot plot shows BHA similarity of CXCR3+/CXCR3-TFH with the CXCR5-TFH population (n = 7). (C) Dot plot shows BHA similarity comparing CXCR3+/CXCR3-TFH with GC-TFH (n = 7); n.s.: not significant. (D) Circos plots show the overlap of cTFH TCR with LN T cell populations. Two donors were analyzed. Each small slice of the arc represents one TCR clonotype, sorted by its size (darker green for larger clones, inner circle). The gray curves link overlapping TCR nucleotide clonotypes in naïve (black, outer circle), CXCR5-TFH (dark blue, outer circle), CXCR3+TFH (blue, outer circle), CXCR3-TFH (red, outer circle) and cTFH (dark green, outer circle). The most expanded 30% clonotypes were highlighted with darker gray curves. (E) Histogram plot shows percentage of TCR clonotype overlapping comparing CXCR3+/CXCR3-TFH with cTFH. (F) Histogram plot shows BHA similarity comparing CXCR3+/CXCR3-TFH with cTFH.

Figure 5

HIV transcript mapping shows CXCR3+TFH is one major reservoir of HIV viruses. (A) Dot plot shows the distribution of normalized HIV copy number across donors and T cell populations (n = 5); n.s.: not significant; (B) Scatter plot shows the correlation between normalized HIV copy number versus Normalized Shannon Entropy (NSE). Blue line represents the regression line, while the dashed area shows the 95% confident interval.

Results

Mass CyTOF Reveals an Enrichment of Non-Canonical CXCR3-Expressing TFH Unique to HIV+ Lymph Nodes

We first set out to determine if, in concordance with SIV infection, CXCR3+TFH were inflated in human LN during HIV infection. We used a 29 marker Mass CyTOF panel to survey T cells from 22 HIV+ and 9 healthy donor LNs (Characteristics of the donors are listed in ). We compared the frequencies of four TFH populations using the gating strategies described in . Two of these populations, denoted as MEM TFH and CXCR3-TFH, are defined as CD45RO+CXCR5+ and CD45RO+CXCR5+PD1+CXCR3- CD4+ T cells, respectively. These two gating schemes are frequently used when characterizing TFH (8, 32, 48). The remaining two populations, CXCR3+TFH and CXCR5-TFH, are defined as CD45RO+CXCR5+PD1+CXCR3+ and CD45RO+CXCR5-PD1+ICOS+ CD4+ T cells, respectively. CXCR3+TFH have been characterized in rhesus macaques (RMs) and human blood (27, 31, 32), while CXCR5-TFH were recently described as functional TFH in human LN and blood during HIV infection (9). Amongst the four populations analyzed, CXCR3+TFH and CXCR5-TFH were significantly enriched in HIV+ patient LNs versus healthy donor LNs (). It should be noted, however, that all LNs from HIV+ donors were cervical, whereas HD LNs were derived from a mix of distinct anatomical locations (). Regardless, given their high frequency over steady state, we suspected these two TFH subsets may play a role in HIV infection. To directly address this question, we compared the frequencies of TFH populations with each patient’s plasma viral load (pVL) (). Interestingly, the only metric correlated to pVL is the ratio of CXCR3+TFH to CXCR3-TFH (), suggesting that donors with greater CXCR3+TFH populations may be better at controlling virus.

To gain insight into the role and interplay of these TFH populations in HIV infection, we sought to delve deeper into the Mass CyTOF data. To do so, we subset the data by drilling down on four T cell populations (CXCR3+TFH, CXCR3-TFH, CXCR5-TFH and Naïve) described in , effectively sorting these populations in silico. We then projected the in silico sorted populations onto two dimensions using Uniform Manifold and Projection (UMAP). At a first glance, cells derived from HIV+ donors occupied several distinct locations on the UMAP (). Additionally, although modest heterogeneity was observed between the three in silico sorted TFH populations, each tended to occupy different regions at different densities (). We thus decided to use unsupervised clustering to gain unbiased insight into each of the populations and their relationships with one another. Twelve clusters emerged, several representing canonical CD4+ T cell and TFH phenotypes (). Select signature markers (), as well as the heatmap in were used to annotate clusters 0-8. Clusters 9-11 were not considered in downstream analyses due to their insignificant sizes. To evaluate the impact of ART on phenotype distribution, we next calculated the frequencies of ART+ and ART- patients in each cluster (). Although some clusters show some trending, we found no significant phenotypic skewing and thus proceeded to group both ART+ and ART- patients as HIV+ in subsequent analyses of this dataset.

To determine unique features of HIV infection, we compared the distribution of T cells in HIV+ and healthy LN-derived cells in each cluster (). In agreement the elevated expansion of the TFH compartment seen in the literature (6, 7), we found that cluster 5, denoted as proliferating TFH based on high expression of Ki67 and BLyS, was significantly enriched in HIV+ donor LN cells (p < 0.001). Surprisingly, even though cluster 2 had the highest and most homogeneous expression of CXCR3 (), we found no difference between HIV+ and healthy cells. This might reflect that many of the CXCR3+TFH in HIV infection take on diverse phenotypic programs. Specifically, given the proximity of cluster 2 with the NAÏVE cluster, it is possible that the CXCR3+TFH within cluster 2 represent early TFH entering the LN as previously described (20), whereas CXCR3+TFH taking on alternative phenotypic programs may represent other stages of the TFH lifecycle. In line with this argument, we reasoned that cluster 5, in addition to being enriched in HIV+ donors, might also have signatures unique to HIV infection. To test this hypothesis, we subsampled cluster 5 to have equal numbers of HIV+ and healthy donor cells and ran a likelihood-ratio test on each marker (49). We found that CXCR3 expression in cluster 5 was in fact significantly higher in HIV+ than healthy donors (). Given the HIV-intrinsic upregulation of CXCR3 in cluster 5, we reasoned that the in silico sorted populations might have differential distributions within cluster 5. As expected, CXCR3+TFH were significantly enriched in cluster 5 over all other populations in HIV+ patients, but not healthy donors (). Concordantly, CXCR3+TFH from HIV+ donors occupied cluster 5 at a significantly higher frequency than from healthy donors, whose predominant population was contrastingly CXCR3-TFH (). In addition to CXCR3, Granzyme A, a marker indicative of cytotoxic function in T-helper cells (50), was also enriched in cells from HIV+ donors within cluster 5 (). It is possible that, although CXCR3+TFH generally reside in a more quiescent, immature state in the steady state, HIV infection triggers them to activate, proliferate in the LN, and take on phenotypic signatures unique to HIV infection.

To better understand the relationship between the phenotypic signatures that the in silico sorted TFH populations () take on specifically in HIV infection, we did pairwise comparisons of the distributions of CXCR3+, CXCR3-, and CXCR5- TFH populations within HIV+ donors across each cluster (). Reflecting the cellular plasticity of TFH, gating down on only a handful of markers unsurprisingly revealed modest heterogeneity within each population. As expected, however, cluster 2, defined by the highest expression of CXCR3, bore a significantly large proportion of CXCR3+TFH compared to all other populations. In line with previous studies (8, 9), CXCR5-TFH were significantly skewed toward activated states in clusters 5 and 6 (Proliferating TFH and EOMES+TFH, respectively), while CXCR3+ and CXCR3-TFH both tended to exist in cluster 7 (Activated GC-TFH). Although statistical significance was only reached between CXCR3+TFH and CXCR3-TFH in cluster 2, we noticed that CXCR3+TFH often existed at a frequency between CXCR3-TFH and CXCR5-TFH, suggesting this population might exist as a transitional state between the canonical TFH and the recently characterized CXCR5-TFH. Given this possibility, we next evaluated the median position of each HIV+ donor’s populations on the same UMAP coordinates as in (). At first glance, the Naïve population is clearly distinct from the TFH populations, which group together on the right side of the UMAP. We also noted that, as suggested in , CXCR3+TFH appear at an intermediate position between CXCR3-TFH and CXCR5-TFH. Statistical analysis of each population also revealed that all populations, except CXCR3+TFH and CXCR5-TFH, were distinct over UMAP_1 (). Additionally, hierarchical clustering of each population from each donor revealed a grouping of CXCR3+TFH with CXCR5-TFH more than with CXCR3-TFH (). In concordance with our hypothesis, these data also posit CXCR3+TFH as an interim population between CXCR3-TFH and CXCR5-TFH.

Transcriptomic Analysis Reveals CXCR3+ and CXCR3-TFH Are Transcriptionally Distinct but Similar on Canonical TFH Marker Genes

Our findings in the Mass CyTOF dataset prompted us to gain deeper insight into the similarities and differences of these TFH populations in HIV infection. Although Mass CyTOF can be useful for delineating cell states due to its ability to analyze large numbers of single cells rapidly and robustly for a selected panel of surface markers, it failed to provide an unbiased analysis of the transcriptome differences between samples. Additionally, due to the destructive nature of Mass CyTOF, TCR sequences, and thus clonal relationships, cannot be obtained. To circumvent these limitations, we sorted five populations of T cells for bulk RNA sequencing (RNA-Seq) from 7 HIV+, non-ART-treated donors using the gating strategy illustrated in . We excluded the one ART-treated donor in this section as we were unsure how treatment would affect a highly sensitive assay such as RNA-seq. We focused specifically on five populations. The first two, CXCR3+TFH (CD45RO+CXCR5+CCR7-PD1+CXCR3+), and CXCR3-TFH (CD45RO+CXCR5+CCR7-PD1+CXCR3-), were sorted to delineate the specific effects of CXCR3 expression on the state of LN TFH. GC-TFH (CD45RO+CXCR5+PD1+CD57+) and CXCR5-TFH (CD45RO+CXCR5-ICOS+PD1+) were sorted as two distinct reference points between two subsets that have been demonstrated to provide B cell help (8, 9). CD57-expressing TFH were sorted as a proxy for GC-TFH as they have been shown to be a major subset of active GC cells (51–53), which serves as a positive control for TFH function. An initial analysis on CyTOF and Flow Cytometry data shows comparable expression of CD57 between CXCR3+TFH and CXCR3-TFH populations (). Naïve CD4+ T cells CD45RO-CXCR5-CCR7+) were sorted as a non-TFH control.

Given a clear enrichment of CXCR3+TFH within HIV+ patient LNs in the Mass CyTOF dataset, we wanted to understand their transcriptional characteristics in HIV infection. Of particular concern, existing studies with animal LN (27) and human blood (31) show contradicting results on their potential for B cell help. We thus set out to explore the transcriptional landscape of CXCR3+TFH in relation to the other, better characterized, TFH populations described above using RNA sequencing (RNA-Seq). We projected the RNA-Seq data from each population onto 2-dimensional space using Principal Component Analysis (PCA). As expected, the naïve population was distinct from all the other populations based on PC1. The TFH populations, although similar on PC1, were separated based on PC2 (). In concordance with our analysis on the Mass CyTOF data, CXCR5-TFH were distinguishable from GC-TFH and CXCR3-TFH on PC2, while CXCR3+TFH situated in the middle. These data again suggest CXCR3 expression on TFH may be indicative of a transitional state between the canonical GC-TFH and the recently described motile CXCR5-TFH (9). Further advocating for CXCR3+TFH as a transitional state, we also noticed CXCR3 expression in both flow cytometry and Mass CyTOF was higher in CXCR5-TFH than CXCR5+TFH (). To better understand the transcriptome differences among various TFH populations, we next used Metascape (43) to evaluate the pathways enriched in genes contributing to PC2 (). Interestingly, genes that mostly negatively correlated with PC2 (upregulated in CXCR5-TFH while downregulated in GC-TFH) were enriched in cell division and migration. However, genes that most positively correlated with PC2 (upregulated in GC-TFH while downregulated in CXCR5-TFH) were enriched on pathways inhibiting cell proliferation and migration (). For example, CTLA4 transmits inhibitory signal to T cells proliferation (54), SMAD3 mediates inhibition of CD4 T-cell proliferation (55), MCC blocks cell cycle progression from G0/G1 to S phase (56). In support of our CyTOF dataset, CXCR3+TFH and CXCR5-TFH appear to be in a more activated, motile state than their canonical GC-TFH counterpart.

Figure 2

CXCR3+TFH population and CXCR3-TFH are functionally distinct. (A) PCA visualization of the transcriptional data of different T cell subsets. Square-plus dots represent the well-established T cell populations while circle dots represent the TFH subsets. Numbers on the axis represents the proportion of variance. (B) Volcano plot of DEGs between CXCR3+TFH and CXCR3-TFH. Highlighted genes are upregulated genes of interest within CXCR3+TFH population. DEGs were identified as log2(Fold Change) larger than 1 and adjusted P value (Benjamini-Hochberg adjustment) < 0.05. (C) Heatmap showing the expressing of DEGs.

To evaluate the differences between CXCR3+ and CXCR3-TFH more directly, we next analyzed differentially expressed genes (DEGs, ) between the two populations. In total, 86 out of 12,515 genes were identified as significant DEGs (). Hierarchical clustering based on those DEGs further demonstrated the similarity of CXCR3-TFH with GC-TFH and CXCR3+TFH with CXCR5-TFH (). A significant upregulation of canonical TH1 genes (CXCR3, IFNG, CXCR6 and GBP1) in CXCR3+TFH also indicated their unique TH1-like program. Additionally, upregulation of NELL2 and MKI67 in CXCR3+TFH suggested an increased proliferative capacity compared to their CXCR3- counterpart. Interestingly, CXCR3+TFH also upregulated ITGB7, encoding a subunit of the integrin α4β7, that plays a role in leukocyte adhesion and can serve as a homing receptor bound by HIV (27). Thus, it is possible that CXCR3+TFH could be similarly, or even more, susceptible to HIV infection than other TFH subsets. Surprisingly, IL10, encoding an unconventional cytokine (IL10) in TH1 cells, was also highly expressed by CXCR3+TFH. Since previous studies support IL10 as a key player in the establishment and perpetuation of HIV persistence (57), the upregulation of IL10 in CXCR3+TFH may result in enhanced persistence of HIV after infection.

To infer the functional potential of human LN-derived CXCR3+TFH, we next compared them to CXCR3-TFH specifically on TFH-related genes. We created a TFH signature gene set using previously published RNA-Seq data from human tonsil samples (31). Both CXCR3+TFH, CXCR3-TFH and GC-TFH appeared similar when hierarchically clustered on tonsil TFH signature genes (). GSEA and GSVA based on the same gene set also revealed no significant differences between these two cell populations (). Detailed analysis on several manually curated TFH-related genes also suggested CXCR3+ and CXCR3-TFH bear similar TFH marker gene expression patterns (). For example: FOXO1, a negative regulator of BCL6 ( 58), was downregulated in both populations. Additionally, BCL6, MAF and CD84, key transcriptional regulators of TFH differentiation, were similar in both CXCR3+ and CXCR3-TFH. B3GAT1 (an enzyme necessary for the production of CD57), which is specifically expressed by active GC-TFH (8, 51–53), was also comparably expressed in both populations. As expected, CXCR5 and its upstream regulator ASCL were both similarly expressed in CXCR3+ and CXCR3-TFH. Taken together, these observations reveal that CXCR3+TFH have a similar TFH transcriptional program to CXCR3-TFH, suggesting a functional overlap between the two populations.

In summary, global transcriptome analysis depicts CXCR3+TFH unique from GC-TFH, biasing toward a TH1-like program. It also revealed CXCR3+TFH are similar to the recently described CXCR5-TFH. Focusing specifically on TFH marker genes, significant differences between CXCR3+ and CXCR3- TFH were not observed, supporting the paradigm that human LN-derived CXCR3-expressing TFH may still bear TFH function.

CXCR3+TFH Upregulate Cell Migratory Pathways and Uniquely Express a Signature of CXCR5-TFH

To investigate deeper into the functional program of CXCR3+TFH, we used NetworkAnalyst (44) to identify GO biological pathways enriched within this population. Only two GO pathways were significantly enriched within CXCR3+TFH when compared with CXCR3-TFH. Interestingly, both pathways were related to T cell migration (). We visualized the differentially expression genes of these two pathways in detail, as expected, most of the genes were upregulated in CXCR3+TFH (; ), which suggested a possibility of this cell subset to be more motile.

Figure 3

CXCR3+TFH uniquely express cell migratory genes and transcriptomic-ally similar with CXCR5-TFH. DEGs between CXCR3+TFH and CXCR3-TFH were enriched on T cell migration pathways (A) with migratory genes upregulated in CXCR3+TFH (Leukocyte chemotaxis, GO:0030595) (B) GSEA analysis indicates genes elevated in CXCR5-TFH (Gene set 1) were upregulated in CXCR3+TFH (C) while genes upregulated in GC-TFH (Gene set 2) were also elevated in CXCR3-TFH (D) GSEA were performed with R package ‘fgsea’, specifically, genes were pre-ranked by their p-values, which indicates whether one gene is highly expressed in CXCR3+TFH or CXCR3-TFH. The bar on the X-axis indicates one gene from a selected gene set (i.e., either Gene set 1 or Gene set 2). The curve is a running sum of the bars on the X-axis.

Our previous study (9) suggested that TFH can downregulate CXCR5 expression and accumulate as CXCR5-TFH in LNs during HIV infection. The same study also found that these CXCR5-TFH provide B cell help and have a propensity to migrate into the periphery. The combination of CXCR3+TFH also bearing a cell migration signature, as well as existing in an intermediate cell state between CXCR5-TFH and GC-TFH based on both Mass CyTOF and RNA-Seq data, led us to hypothesize that CXCR3+TFH may be the immediate relative of CXCR5-TFH. To test this hypothesis more directly, we compared the transcriptomes of CXCR3+ and CXCR3-TFH using GC-TFH and CXCR5-TFH as reference gene sets. Specifically, we first identified 938 DEGs from 12,559 total shared genes between GC-TFH and CXCR5-TFH. Among these DEGs, 514 were upregulated in CXCR5-TFH (Gene set 1) and 424 were upregulated in GC-TFH (Gene set 2). We then performed GSEA analysis to compare CXCR3+ and CXCR3-TFH on Gene Set 1 and 2 (). As hypothesized, CXCR3+TFH followed the CXCR5-TFH program, while CXCR3-TFH more closely followed the GC-TFH program.

Together with our observations in Mass CyTOF, global transcriptome, differential gene expression, gene set enrichment, and pathway analyses, we conclude that CXCR3+TFH likely exist in a phenotypically intermediate state between canonical GC-TFH and the more motile and peripheral CXCR5-TFH.

T Cell Receptor Repertoire Demonstrates That CXCR3+TFH Are Clonally Related to CXCR5-TFH and Peripheral TFH

After observing that CXCR3+TFH upregulated cell migratory pathways and skew phenotypically toward CXCR5-TFH, we suspected that changes in CXCR3 expression on TFH might facilitate transitions to and from GC-TFH and CXCR5-TFH cell states. To test this hypothesis, we used bulk T cell receptor (TCR) repertoire sequencing to measure their in vivo clonal relationship with various T cell populations. Because of the immense diversity generated by V(D)J recombination, TCR sequences can be thought as unique ‘ID-Cards’ (59), where cells sharing a given TCR sequence are likely to be the progenies of the same ancestor. By comparing the overlap of TCR sequences among samples, a clonal lineage across tissues and cell states can be inferred.

We thus compared the TCR repertoire similarity across the five LN T cell populations using Bhattacharyya Coefficient (36) as the similarity index. Bhattacharyya coefficient (BHA) measures the overlap of clonotypes between two T cell populations, while taking size of clones into consideration. In contrast to previous findings in healthy human tonsils (23), CXCR3+TFH, CXCR3-TFH and GC-TFH were strikingly similar in their repertoires across all donors (; ). This suggests that these three populations are likely capable of transitioning in and out of each of these states with ease, which may be a phenomenon intrinsic to HIV infection. In concordance with transcriptomic data, however, CXCR3+TFH were significantly more related to CXCR5-TFH than CXCR3-TFH across all donors (). For two out of the seven donors, we were able to compare the CXCR3+ and CXCR3-TFH repertoires with the cTFH repertoires in the blood (60) (; ). Circos plots representing the overlap of TCRs between cTFH and LN populations indicate that CXCR3+TFH are also more clonally related to cTFH than their CXCR3-TFH counterpart (), similar to a recent finding in healthy tonsils (23). A similar trend was observed considering clone size weighted BHA index (). In combination, high clonal overlap amongst all CXCR5-expressing TFH in the LN, and the elevated relationship of CXCR3+TFH with CXCR5-TFH and cTFH, suggests that this population sits as a bridge between secondary lymphoid tissue and the periphery.

CXCR3+TFH in the Human LN Are a Reservoir for HIV

TFH are major reservoirs for HIV (61), complicating their role in controlling virus through the coordination of GC reactions. Furthermore, CXCR3+TFH have been speculated to maintain a dynamic HIV reservoir due to high expression of HIV co-receptors (i.e. CCR5 and α4β7) (62). High viral load in this highly motile cell type could either exacerbate systemic spread of virus or expose viral antigens to the immune system. Given our result that CXCR3+TFH frequency negatively correlates with pVL (), we speculated a high level of HIV copies in CXCR3+TFH might be a mechanism for improving clearance of virus.

To investigate this idea, we mapped the bulk RNA-Seq reads from each T cell population to the HIV genome to quantify their relative HIV infection intensities (see methods for details). We detected comparable levels of HIV transcripts in both CXCR3+ and CXCR3-TFH, although both were significantly higher than Naïve cells (; ). Similar to experiments in RMs (27), HIV abundance also negatively correlated with repertoire diversity (), reflecting the propensity of proliferating and antigen-experienced cells to bear a higher HIV burden. Given that the two TFH subsets discussed here share a similarly high HIV burden, the consequence of the migratory potential of CXCR3+TFH is likely a complex balance between limiting systemic spread of virus and exposing it to novel immune compartments.

Discussion

TFH are paramount to the elicitation of BNAbs. Whether induced by vaccination or natural infection, BNAbs targeting invading pathogens have been shown in a variety of contexts to mitigate, eliminate, and even prevent disease (63, 64). Some viruses, however, can escape from or disrupt this response, often leading to uncontrolled viral replication and chronic infection. In HIV infection, natural BNAbs are seldom produced and almost never lead to viral control (4, 5). Given that TFH are a major site of HIV infection (11–13), it is possible that alternative phenotypic and functional programs within the TFH compartment are taken on to overcome this deficit in humoral immunity. Unfortunately, sampling TFH from humans in their native lymphoid tissue in untreated HIV infection poses an immense challenge in the field. Here, we have overcome this obstacle, combining several high-throughput assays and analyses on TFH from HIV+ patient LNs to gain deeper insight into their role in disease.

Mounting evidence points to a high degree of cellular plasticity within the TFH compartment (15–17). We show here the inflation of CXCR3+TFH in HIV+ patient derived LNs compared to healthy donors that correlates with a lower viral load. Further investigation into this population revealed its abnormal cell state specifically within HIV+ donors that positions it in an intermediate phenotypic state between GC-TFH (8, 31) and the recently described CXCR5-TFH (9). Deeper analysis into the clonal relationship of CXCR3+TFH corroborated these results, revealing its heightened repertoire similarity with both CXCR5-TFH and cTFH. Additionally, transcriptome analysis revealed a propensity for CXCR3+TFH to upregulate cell migration pathways. Taken together, we posit CXCR3+TFH as a bridge between lymphoid tissue and the periphery. Given that these CXCR3+TFH bear a high HIV burden, may be primed for cellular migration, and are affiliated with a lower viral load, it is possible that their movement in and out of lymphoid tissues leads to better viral control by exposing antigens to novel immune compartments. If so, it is possible that the upregulation of CXCR3 in TFH may be an effort of the immune system to drive viral reservoirs out of hiding. Given the propensity of low levels of HIV to remain through ART (12) and the increased likelihood that BNAbs are produced during high antigen exposure (10), a mass exodus of TFH may be beneficial to the host during early infection or ART. Furthermore, targeting these migration pathways has the potential to be a viable therapeutic. Of course, it is also possible CXCR3 upregulation in TFH is not HIV-intrinsic, but rather a product of chronic infection. In fact, previous research has implicated CXCR3+TFH in HCV infection (60), Zika virus infection (65), and acute febrile malaria (66) and vaccination (67) in children. Interestingly, however, while CXCR3+TFH are affiliated with neutralizing antibodies in HCV and Zika virus infection, they appear to lead to poor prognosis in malaria. Understanding the differences and similarities between CXCR3+TFH in each of these disease contexts, specifically in lymphoid tissues as addressed in this study, will shed light on different immune-intrinsic properties of chronic infection.

An important limitation of this study was our inability to directly measure TFH function on the populations of interest. Although we were able to point to specific genes, pathways, and signatures that suggest their ability to provide B cell help, definitive knowledge of the functional capacity of CXCR3+TFH in human LNs will be important to understanding their role in HIV infection. Future studies should aim at accomplishing this. Additionally, the spectral limitations of FACS also prevented us from being able to sort more specific TFH subsets that may be important to resolve the unique phenotypic programs of CXCR3+TFH in HIV infection. For example, our analyses on the Mass CyTOF data suggested that CXCR3+TFH take on a unique proliferative program specific to HIV+ patients that was not robustly seen in healthy donors, a phenomenon that may explain some of the repertoire discrepancies between this study and others (20). However, since we were unable to accommodate more markers when sorting for TCR and RNA sequencing, especially with the limited number of cells within each LN sample, these unique HIV-intrinsic factors may have been diluted. Recent innovations in antibody barcoding and single cell RNA-seq coupled with creative approaches to isolate these unique TFH populations may lead to a clearer understanding of their phenotypic and clonal relationships.

In summary, we have evaluated several LN-blood-matched TFH from a small, albeit rare cohort of untreated HIV+ patients using Mass CyTOF and RNA and bulk TCR repertoire sequencing. Our analyses revealed a phenotypic shift of CXCR3+TFH from a GC-TFH cell state toward an unconventional CXCR5-TFH state. CXCR3+TFH also upregulated migratory transcriptional programs and were clonally related to peripheral TFH cell populations. Altogether, these data suggest that CXCR3+TFH may be transitional state between their CXCR3- lymphoid and peripheral counterparts. Future work aimed at delineating the temporal relationships of these TFH populations, drilling deeper into the function of more specific phenotypic niches, will be pertinent to fully understanding their implication in HIV infection.

Data Availability Statement

Publicly available datasets along with new sequencing data were analyzed and can be found here: https://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_id=phs001548.v2.p1.

Ethics Statement

Subjects included in this study are from a subset of patients recruited for our previous study (8). All samples were de-identified and obtained with IRB regulatory approval from the University of Pennsylvania. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

CH and MJM analyzed data, performed research, and wrote the manuscript; MJM, BSW and K-YM performed sequencing experiments; DDA performed cell sorting experiments; DBW and PLDJ helped with data interpretation; PMDR-E, YA-T, and GR-T established the infrastructure to recruit HIV+ patients, provided HIV-infected samples, and associated clinical information. NJ and LFS designed the study; NJ directed the study. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by NIH grants S10OD020072 (NJ), R56AG064801 (NJ), R01AI134879 (LFS), NIH IPCAVD grant U19 Al109646-04(DBW), Chan Zuckerberg Initiative Neurodegeneration Challenge Network Ben Barres Early Career Acceleration Awards 191856 (NJ), VA Merit Award IMMA-020-15F (LFS), NIH/NIAID Collaborative Influenza Vaccine Innovation Centers (CIVICs) contract 75N93019C00051(DBW). DBW is the W.W. Smith Charitable Trust Professor at the Wistar Institute.

Conflict of Interest

NJ is a Scientific Advisor and holds equity interest in ImmuDX, LLC, and Immune Arch, Inc. In the interest of full disclosure, DBW reports that he serves on Advisories for AstraZeneca, Geneos, Advaccinepharma, he participates in BOD service for Inovio. Remuneration received by DBW for these services includes SRA funding, direct payments, stock or stock options disclosed here.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.