Th2 single-cell heterogeneity and clonal distribution at distant sites in helminth-infected mice

Version of Record

Accepted for publication after peer review and revision.

Download
Cite
Share
CommentOpen annotations (there are currently 0 annotations on this page).

Version of Record published: August 19, 2022 (This version)
Accepted Manuscript published: August 11, 2022 (Go to version)
Accepted: August 10, 2022
Preprint posted: October 3, 2021 (Go to version)
Received: September 24, 2021

1. Of interest
The CD4 transmembrane GGXXG and juxtamembrane (C/F)CV+C motifs mediate pMHCII-specific signaling independently of CD4-LCK interactions

Mark S Lee, Peter J Tuohy ... Michael S Kuhns

Research Advance Apr 19, 2024
Further reading

Abstract
Editor's evaluation
Introduction
Results
Discussion
Materials and methods
Data availability
References
Article and author information
Metrics

Abstract

Th2 cells provide effector functions in type 2 immune responses to helminths and allergens. Despite knowledge about molecular mechanisms of Th2 cell differentiation, there is little information on Th2 cell heterogeneity and clonal distribution between organs. To address this, we performed combined single-cell transcriptome and T-cell receptor (TCR) clonotype analysis on murine Th2 cells in mesenteric lymph nodes (MLNs) and lung after infection with Nippostrongylus brasiliensis (Nb) as a human hookworm infection model. We find organ-specific expression profiles, but also populations with conserved migration or effector/resident memory signatures that unexpectedly cluster with potentially regulatory Il10^posFoxp3^neg cells. A substantial MLN subpopulation with an interferon response signature suggests a role for interferon signaling in Th2 differentiation or diversification. Further RNA-inferred developmental directions indicate proliferation as a hub for differentiation decisions. Although the TCR repertoire is highly heterogeneous, we identified expanded clones and CDR3 motifs. Clonal relatedness between distant organs confirmed effective exchange of Th2 effector cells, although locally expanded clones dominated the response. We further cloned an Nb-specific TCR from an expanded clone in the lung effector cluster and describe surface markers that distinguish transcriptionally defined clusters. These results provide insights in Th2 cell subset diversity and clonal relatedness in distant organs.

Editor's evaluation

CD4 Th2 effector cells contribute to immune responses to helminths and allergens. There exists little information on Th2 cell heterogeneity and clonal distribution between organs. In this manuscript, Radtke et al. investigate the transcriptional signatures of CD4 Th2 cells in the mesenteric lymph nodes and lungs during helminth infection. By using single cell RNA-sequencing including TCR clonotype analysis, the authors define distinct and overlapping transcriptional signatures and clonal relatedness between CD4 Th2 cells in two different tissues at the peak of a type 2 immune response in vivo.

https://doi.org/10.7554/eLife.74183.sa0

Introduction

Th2 cells are part of the adaptive immune response against helminths and in allergic diseases. They are recruited and differentiate from a pool of naive CD4 T cells with a wide variety of T-cell receptors (TCRs) that are formed during T-cell development and provide clonotypic specificity to antigens. Differentiated Th2 cells produce the key type 2 cytokines IL-4, IL-5, and IL-13 that elevate type 2 immune responses and thereby promote allergic inflammation but also mediate protection against helminths (Walker and McKenzie, 2018). In recent years, several IL-4 producing Th2 subpopulations have been described and point to substantial heterogeneity within the Th2 population. Only a minor fraction of human IL-4⁺ T cells produces IL-5 which defines them as highly differentiated cells (Upadhyaya et al., 2011). In the mouse, Th2 cells in the lung generally appear more activated and coexpress IL-4 and IL-13 as compared to Th2 cells isolated from lymph nodes of helminth-infected mice (Liang et al., 2011). Th2 cells can further differentiate to follicular T helper cells that express IL-4, IL-21, and BCL6 and drive humoral type 2 immune responses in the germinal center (GC) (Glatman Zaretsky et al., 2009; King and Mohrs, 2009; Reinhardt et al., 2009). On the other hand, IL-4 producing T follicular helper (Tfh) cells may develop into strong effector Th2 cells in asthma models (Ballesteros-Tato et al., 2016; Tibbitt et al., 2019). In addition to IL-4 producing Tfh cells, IL-4 secretion by T cells located outside of GCs can be sufficient for GC formation and class switch recombination to IgE (Turqueti-Neves et al., 2014). Tfh13 cells may also develop from Th2 cells in settings or allergic inflammation. These cells coexpress IL-4, IL-5, and IL-13 and promote the generation of high-affinity anaphylactic IgE in response to allergens (Gowthaman et al., 2019). In addition to these subsets with distinct functions, there are different activation and developmental stages present in the Th2 population. A rough model describes Th2 cell development as a process that involves activation of naive T cells, followed by proliferation and subsequent production of effector cytokines associated with a high proliferation rate (Proserpio et al., 2016). Furthermore, fate mapping and adoptive transfer experiments revealed functional plasticity between T helper cell subpopulations which can lead to Th2 cells with remaining or upcoming signatures of other CD4 T-cell subsets like Th1, Th9, or Th17 cells (Panzer et al., 2012; Peine et al., 2013; Tortola et al., 2020; Veldhoen et al., 2008).

Infection of mice with the helminth Nippostrongylus brasiliensis (Nb) is a widely used model for human hookworm infections with a strong induction of Th2 responses in lung and small intestine (Urban et al., 1992). L3 stage larvae are injected subcutaneously and then first migrate into the lung before they are coughed up, swallowed and reposition to the small intestine where they mature to adult worms. (Urban et al., 1992). Using this model, we have previously shown that Nb infection induces a Th2 response with a broad TCR repertoire required for effective worm expulsion (Seidl et al., 2011). Development of single-cell sequencing technology now allowed us to gain a deeper understanding of Th2 cell subsets, TCR clonality, and tissue distribution.

Here, we performed single-cell sequencing of TCR genes combined with transcriptome profiling of Th2 cells isolated from mesenteric lymph nodes (MLNs) and lung of IL-4eGFP reporter mice (4get mice) at day 10 after Nb infection. In this manuscript, we use the term ‘Th2’ for the CD4⁺, IL-4eGFP⁺ cells. This likely includes cells licensed for IL-4 expression that will not develop into terminally differentiated Th2 cells and retain functional plasticity, and NKT cells. By our approach, we revealed heterogeneity and differentiation paths within the Th2 compartment, compared Th2 population similarity at distant sites, analyzed cell exchange between organs by clonal relatedness and characterized expanded clones and their TCR sequences.

Results

Th2 cells show an organ-specific gene expression profile consistent with acquired effector functions

We performed combined transcriptome and TCR clonotype analysis using the chromium 10× Genomics and Illumina single-cell RNA sequencing (scRNAseq) platform on IL-4-expressing Th2 cells isolated from perfused lung (also containing a remaining fraction of intravascular cells) and MLN of two IL-4eGFP reporter (4get) mice (Mohrs et al., 2001) that had been infected 10 days before with Nb (Figure 1A). IL-4-expressing Th2 cells (CD4⁺IL-4eGFP⁺) were sorted from single-cell suspensions of both organs and were directly subjected to scRNA library preparation.

Figure 1 with 3 supplements see all

Download asset Open asset

Th2 cells of mesenteric lymph node (MLN) and lung adopt tissue-specific RNA signatures.

(A) General experimental outline. MLN and lung cells of two individual *Nippostrongylus brasiliensis* (Nb)-infected IL-4eGFP reporter mice (4get B6) were sorted for IL-4eGFP⁺CD4⁺ cells 10 days post infection. Then combined transcriptome and T-cell receptor (TCR) repertoire sequencing were performed. Flow cytrometry plots show the frequency of Th2 cells (IL-4eGFP⁺CD4⁺ cells) in MLN and lung. (B) Uniform Manifold Approximation and Projection (UMAP) representation of MLN and lung cells 10 days post Nb infection. (**C, D**) medLN and MLN cells sequenced in a separate run were plotted on the existing UMAP defined in the initial MLN + lung sequencing run by integration based on shared anchor genes to demonstrate similarity of medLN and MLN cells. (E) De novo unsupervised clustering approach with manually added cell type description. Clusters are indicated on UMAP. (F) Expression of selected CD4 T-cell subset defining genes, (G) genes that are differently expressed between MLN and lung, or (H) gene signature module scores for single cells plotted on top of UMAP representation. Each of the independent single-cell sequencing experiments is based on two mice.

4get mice were chosen as they allow isolation of Th2 cells ex vivo without prior restimulation. In contrast to other IL-4 reporter mice such as the KN2 strain, 4get mice even report the early stages of Th2 differentiation (Mohrs et al., 2005). Sampling of the lung was performed as Th2 cells accumulate in this organ a few days after Nb infection. Complimentary MLNs were included as a distant secondary lymphoid organ associated to the intestine where the worms reside approximately from days 4 to 10 after infection. By choosing a secondary lymphoid organ and a peripheral organ we cover location-specific differences as well as a broad spectrum of developmental stages. In both organs worms of larval stage 4 are present in the course of the infection (Camberis et al., 2003) and provide a common antigen basis in addition to antigens shared between larval stages or systemic dissemination of antigens from larval secretions or dead larvae. This setup enabled us to compare Th2 cell subsets and clonotypes in both organs at single-cell resolution.

In order to restrict the analysis after sequencing to high-quality Th2 cells, we included in total 4710 cells with detected, functional TCR α- and β-chains that also passed our additional QC filters (see Materials and methods) (Figure 1—figure supplement 1). We used an unbiased high-dimensional clustering approach followed by dimensional reduction for simple representation of complex data (Stuart et al., 2019).

Our approach revealed that Th2 cells of lung and MLN have a distinct organ-specific expression profile represented by clear separation of cells from both organs upon dimensional reduction (Figure 1B) and highlighted by differential expression analysis between lung and MLN cells (Figure 1—figure supplement 2A). We mainly focused on the analysis of MLN and lung in order to compare the T-cell repertoire between distant organs later on. Importantly, Th2 cells from lung draining mediastinal LN (medLN) and MLN analyzed in a separate scRNAseq experiment (Figure 1—figure supplement 3) showed a similar expression profile when layered on the lung–MLN UMAP of our initial experiment (Figure 1C,D). This indicates that most of the differences are due to the comparison of a secondary lymphoid organ with a peripheral organ, and only minor differences are related to the comparison of distant sites (medLN vs. MLN).

We continued with unbiased clustering of MLN and lung to identify two lung, seven MLN, and a mixed proliferative cluster that are numbered by size. These are summarized here and further described in the following two sections (Figure 1E): L1 (basic activated/effector), L2 (migrating), L + MLN (proliferating), MLN1 (basic activated), MLN2 (contains Tfh cells), MLN3 (IFN response signature), MLN4 (effector/resident memory like), MLN5 (migrating), MLN6 (innate-like/NKT), and MLN7 (myeloid RNA containing Th2). For every cluster, the 10 marker genes that are upregulated in a cluster compared to all other cells in the dataset were determined as a reference (Figure 1—figure supplement 2B).

We first concentrated on markers for known CD4 subpopulations. As expected, the cells from MLN and lung express the Th2 hallmark genes Il4, Gata3, and Stat6 (Figure 1F). However, almost exclusively the L1 cluster expresses IL-5 which promotes eosinophil development, recruitment, and survival. Similarly, IL-13 is expressed a bit broader in L1 but additionally in lymph node cluster MLN4. IL-13 elicits a broad spectrum of effector type 2 immune functions including eosinophilic inflammation, mucus secretion, and airway hyperresponsiveness (Rothenberg and Hogan, 2006; Takatsu et al., 2009). According to the proinflammatory IL-5 and IL-13 production, double producers are thought to be a strong or pathogenic effector subset of Th2 cells that includes highly differentiated CD27^low, PD-1(Pdcd1)^high memory cells (Upadhyaya et al., 2011), which is also reflected on gene expression level in our data. Enhanced Rgs16 expression of the IL-5⁺/IL-13⁺ cells is also associated with higher cytokine production (Lippert et al., 2003) and further supports effective effector molecule production (Figure 1F).

As expected, very few IL-4-positive cells express the Th1 hallmark genes Ifng and Tbx21 (encodes T-bet). However, cells that express IL-4 and both markers are present in cluster MLN6. An overlapping fraction of cells also expresses Zbtb16 which encodes the NKT cell-associated transcription factor PLZF (Savage et al., 2008) and Klrb1c (encodes NK1.1) or in addition to the TCRα and TCRβ chains Tcrg-1c (encodes the constant region of TCRγ) as an indication that these cells show signs of unsuccessful or not yet successful development into γδ T cells. Hence, MLN6 seems to contain predominantly innate-like/NKT cells (Figure 1F).

The classical marker for Treg cells FOXP3 was hardly found on gene expression level in our dataset. Nevertheless, small fractions of cells in L1 and MLN4 express Il10, which suggests regulatory capacity independent of Foxp3 expression (Figure 1F) that has also been reported as relevant to maintain the Th2 response and to repress Th1 (Balic et al., 2006). Interestingly, the potential regulatory IL-10 producers cluster together with the likely proinflammatory IL-5 and IL-13 producers, suggesting a similar expression profile of cells with regulatory and effector function.

Tfh cells that express IL-4 were detected in MLN2. They express both Tfh markers IL-21 and BCL6 on gene level. Of note, Bcl6 expression seems less restricted to a specific cluster. Where Il21 and Bcl6 expression overlaps cells show expression of Rgs16 associated with enhanced Th2 cytokine production and trafficking (Lippert et al., 2003). In line with a recent publication, we did not observe Tfh13 cells (IL-13^hiIL-4^hiIL-5^hiIL-21^lo) which were reported to be associated with production of high-affinity anaphylactic IgE in Th2 responses to allergens but not helminth infections (Gowthaman et al., 2019).

Increased expression of TGFβ-associated genes further suggested Th2 cell plasticity toward Th9 cells in our dataset (Veldhoen et al., 2008; Figure 1H) and IL-9 is described as a relevant factor for hookworm expulsion (Licona-Limón et al., 2013). However, we only find very few Il9 or Spi1 (encodes the Th9-associated transcription factor PU.1) expressing cells in the lung. Il9 expression was also barely detectable in the MLN, while Spi1 was expressed in the MLN4 population potentially describing an early stage of Th9 development (Figure 1F).

Next, we screened for genes that distinguish Th2 cells from MLN and lung and that might reflect organ-specific differences in gene expression or Th2 development. While most cells from the MLN (MLN1–3) express the gene for the chemokine receptor CXCR5 associated with homing to B cell follicles and recruitment of Tfh cells to GCs (Breitfeld et al., 2000; Schaerli et al., 2000), the majority of Th2 cells from the lung and MLN cells that cluster in proximity to lung cells (MLN4–6) hardly express it (Figure 1G). In contrast, the lung and MLN4–6 cluster cells express high levels of the gene for TAGLN2 which stabilizes the immunological synapse and is relevant for proper T-cell effector function (Na et al., 2015). A stronger effector phenotype of lung cells is also supported by an increase of inflammation signature genes and hypoxia-associated genes in these cells, which are associated with enhanced glycolysis required for late Th2 effector differentiation (Healey et al., 2021; Stark et al., 2019; Figure 1H).

In line with strong effector function, most lung cells and some nearby projected MLN cells express the gene for the IL-33 receptor ST2 (Il1rl1) (Figure 1G). It recognizes the alarmin IL-33 and induces production of key type 2 cytokines IL-5 and IL-13. This is crucial for the clearance of intestinal helminths as mice that lack IL-33 are not able to effectively cope with the infection, likely due to defects in the T cell and ILC2 compartments (Hung et al., 2013). In line with enhanced Il1rl1 expression, we find evidence that ST2 inducing pathways are active in lung cells. As such, there is higher expression of IL-2–STAT5 axis target genes and the IL-2 receptor alpha-encoding gene itself (Il2ra) (Guo et al., 2009; Meisel et al., 2001) as well as enhanced tumor necrosis factor receptor 2 (TNF-R2) transcript expression and correspondingly an increased TNF induced expression signature, suggesting stimulation of the cells by TNF (Kumar et al., 1997; Figure 1GH). Looking further downstream, we find an elevated gene signature for IL-33-stimulated T cells in the lung, which suggests active signaling via the ST2 receptor, associated with pathogen clearance. However, IL-33/ST2 not only invokes effector mechanisms against the worm but also promotes production of Areg, encoding amphiregulin, that we find upregulated and which is involved in tissue repair and resolution of inflammation after Nb infection (Minutti et al., 2019; Figure 1G). However, amphiregulin can also promote reprogramming of eosinophils to develop fibrosis in other setups (Morimoto et al., 2018) with potential implications on longer-term pathology after Nb infection (Marsland et al., 2008).

In our Nb infection model, lung cells and a fraction of proximal clustering MLN cells also express genes associated with asthma or involved in pathways targeted by drugs for asthma treatment like Cysltr1, Plac8, or Adam8 (Naus et al., 2010; Tibbitt et al., 2019; Trinh et al., 2019). Lung cells but only few MLN cells in our dataset also show an increased expression of TGFβ target genes, consistent with the described Th2 cell plasticity toward a Th9 phenotype (Veldhoen et al., 2008) which potentially further broadens the T effector functions (Figure 1F–H).

Conserved expression profiles for migratory and effector/resident memory Th2 cell populations in lung and MLNs

As a next step after concentrating on known CD4 subpopulations and comparison of MLN and lung on organ level, we focused on similarities of MLN and lung and clusters with unexpected signatures.

Cluster L2 and MLN5 visually form a ‘bridge’ between the MLN and lung compartments. Indeed cells in these clusters both express genes coding for CD62L (Sell), CCR7, and S1PR1 involved in cell adhesion and T-cell trafficking suggesting that these clusters contain recent immi-/emmigrants (Figure 1E and Figure 2A). Of note, staining of cells in the blood circulation revealed that a fraction of cells in the L2 cluster likely represents cells located in blood vessels with a gene expression profile similar to tissue-resident L2 cluster cells which further indicates that the L2 cells might come from or migrate into blood vessels (Figure 2—figure supplement 1). The ‘bridge clusters’ also express Tcf7 (associated with self-renewal capacity) and Cd27 (encoding a central memory T-cell marker) (Figure 2A). In line, the whole ‘bridge’ shows a circulating memory signature (Rahimi et al., 2020; Figure 2B). L2 part of the ‘bridge’ is the only lung fraction that expresses reasonable amounts of Cxcr5 and Tox2, which are expressed in a majority of cells in the MLN (Figures 1G and 2A). This suggests that the profile of these lung cells in part reflects the profile found in secondary lymphoid organs and strengthens their identification as migrating cells. There are also differences between the lung- and MLN-associated ‘bridge’ clusters. Cells in L2 expressed more CD44, which is suggestive of cells in later central memory or effector cell state and exhaustion marker genes encoding Tox2 and Pdcd1 compared to MLN5. In contrast, to L2 which shows signs of MLN-associated gene expression, MLN5 does not show clear signs of a lung signature (Figure 1F, H and Figure 2A), potentially reflecting that the visual ‘bridge’ is not a real connection and contains immi-/emmigrants to/from other secondary lymphatic organs like the lung draining lymph nodes or other peripheral organs.

Figure 2 with 3 supplements see all

Download asset Open asset

Conserved expression profiles for migratory and effector Th2 populations between organs and their inferred developmental paths.

(A) Expression of selected marker genes associated with biological processes or (B) gene signature scores for single cells on top of UMAPs. (C) RNA velocity analysis. Arrows present inferred developmental paths. (D) RNA velocity defined root cells and developmental endpoints as well as inferred differentiation speed (velocity length) and velocity confidence for cells are visualized on UMAPs.

Proliferating cells of both organs have a similar expression profile and fall into the same cluster (L + MLN) as proliferation induces a strong gene signature highlighted by a proliferation-associated E2F signature gene set (Figure 2B).

We find three effector clusters in our dataset, L1 lung cells express effector-associated genes (e.g., stronger expression of Il1rl1, Cysltr1 [receptor for cysteinyl-leukotrienes C₄, D₄, and E₄], Il2ra, Il5, and Il13) but MLN4 has a similar signature. Both populations are high for a published signature of lung-resident memory T cells (Rahimi et al., 2020; Figure 2B) and for both the effector-associated genes Plac8 and Adam8 drive the signature (Figure 1—figure supplement 2B and not shown). In contrast to L1, most MLN4-resident memory signature cells express the gene-encoding CCR9 relevant for lamina propria homing (Campbell and Butcher, 2002; Stenstad et al., 2006), while lung cells express the gene for CCR8 hardly found in the MLN (not shown). In addition, to MLN4 there is an innate-like/NKT cluster (MLN6) present in the MLN that also shows expression of effector genes. Both clusters express a resident memory signature and Ccr9, probably reflecting local effector/resident memory populations that participate in the intestinal immune response. Of note, there is also a Ccr9-expressing fraction of cells in L1 with little overlap to the Il5/Il13 secreting cells (Figure 2A). These cells, positioned in the lung might come from or might migrate to the lamina propria. Interestingly, the innate/NKT MLN6 effector cluster expresses the gene for the usually Th1-associated chemokine receptor CXCR3 found before (Kim et al., 2001) in a fraction of Th2 cells but they were not linked to innate/NKT phenotypes. Therefore, there is an expected effector cluster in the lung but cells with an effector signature are also found in the MLN which is mainly recognized to contain earlier stages in T helper cell development. Despite the shared effector characteristic of effector-associated clusters, there is site-specific gene expression likely reflecting local effector function. Similarly, there is also evidence for site-specific effector signatures in our second sequencing run comparing medLN and MLN. medLN contains a number of cells that rather cluster with lung effector cells (L1) and a reduced number of cells that show a lamina propria associated effector signature (like MLN4) compared to MLN (Figure 1C,D).

To highlight the similarities of ‘brigde’ (MLN5, L2) and effector clusters (L1, MLN4) between MLN and lung, we integrated the lung and the MLN data computationally based on shared expression anchors before performing dimensional reduction (Stuart et al., 2019). The approach reveals that bridge clusters as well as effector clusters (L1 and MLN4) show the most similar expression profiles between MLN and lung and hence cluster together in low-dimensional space projection (Figure 2—figure supplement 2).

While Th1 and Th2 cells are often seen as counterparts that can antagonize differentiation of each other, there is also evidence that the Th1 hallmark cytokine IFN-γ promotes proper Th2 heterogeneity and differentiation either directly as suggested by in vitro differentiation studies (Wensky et al., 2001) or indirectly by inducing activation of DCs with Th2-priming capacity (Connor et al., 2017; Webb et al., 2017). In accordance MLN3 expresses an IFN response signature after Nb infection (Figure 2B) likely associated with Th2 priming and differentiation.

We further identified an unusual subset of cells (MLN7) which contains genes for the MHCII-associated invariant chain (CD74), complement component C1q, lysozyme, and CD209b. Some of these genes are rather associated with myeloid cells like DCs or macrophages. This population might therefore represent cells emulsified with exosomes or RNA containing vesicles during library generation, either externally attached or taken up by the cell.

In summary, our dataset with high resolution for the Th2 compartment after Nb infection reveals distinct subsets of Th2 cells in lung and MLN including a ‘bridging’ cluster of cells with closely related gene expression profiles in both organs. Defined clusters of Th2 subsets were later used to link TCR clonotypes to specific clusters.

RNA-inferred developmental directionality of Th2 cells supports proliferation as a hub for differentiation decisions

To infer developmental relatedness of clusters defined above we performed an RNA velocity analysis in which the ratio of spliced to unspliced RNA transcripts is used to calculate and visualize likely developmental directions (Bergen et al., 2020; La Manno et al., 2018; Figure 2C). We used the scVelo algorithm (Bergen et al., 2020) which identifies the most undifferentiated cells as root cells and highly differentiated cells as developmental endpoints, connected by arrows that show likely paths from root to endpoints (Figure 2C,D). The proliferation cluster (L + MLN) reflects the majority of root cells in our data and highlights proliferation, in accordance with the literature (Gett and Hodgkin, 1998; Gulati et al., 2020; Radtke and Bannard, 2018), as a critical branching point at which differentiation decisions are taken. The developmental direction going from proliferating cells toward the other populations is also supported when RNA velocity is performed on PCA level (Figure 2—figure supplement 3A), in paga projection (Figure 2—figure supplement 3B), in the static model of RNA velocity (Figure 2—figure supplement 3C), in developmental analysis performed with the monocle package (Trapnell et al., 2014; Figure 2—figure supplement 3D) or when cell cycle was regressed out before RNA velocity analysis (Figure 2—figure supplement 3E). Another algorithm that is described to robustly detect developmental potential also for proliferating cells additionally indicates that the proliferating cluster has the highest developmental potential (Gulati et al., 2020; Figure 2—figure supplement 3F).

The MLN part of the ‘bridge’ clusters (MLN5), the IFN signature cluster (MLN3), and the main MLN (Basic activated Th2; MLN1) cluster are marked as relatively diffuse endpoints in the MLN, associated with a low differentiation speed and confidence reported by scVelo (Figure 2D). This suggests that wide parts of the MLN Th2 cells are relatively heterogeneous. The effector/resident memory like cluster (MLN4) is in itself heterogeneous and contains cells with a strong root signature which hardly overlap with the also contained strong resident memory signature cells. A relatively high differentiation speed and confidence compared to other MLN clusters suggest that it contains a fast developing effector/resident memory like population. Based on the MLN5 ‘bridge’ cluster definition as an endpoint, MLN5 might rather reflect cells that leave the MLN. The lung cluster of the ‘bridge’ (L2) instead contains cells that differentiate with high confidence and inferred speed toward the main lung cluster of effector cells (L1), which suggests that these cells enter the lung and further differentiate locally. The IL-5/IL-13 double producers previously defined as highly differentiated effector cells (Upadhyaya et al., 2011) reflect the endpoint in the lung (Figure 2C and D).

In conclusion, RNA velocity supports proliferation as a hub for differentiation in the Th2 compartment and supports that migratory Th2 cells rather leave secondary lymphatic organs and enter peripheral organs while the reverse migration path is a rare event.

Longitudinal study of Th2 development on protein level based on a transcriptionally defined markers

In an attempt to relate described transcriptional profiles and cluster definitions to distinguishable surface receptor signatures, we analyzed expression of identified surface marker genes on the protein level by flow cytometry. We included MLN and lung but also medLN to reveal differences in Th2 development at different sites over time. First, we selected differentially expressed surface markers that allow to distinguish defined clusters and analyzed their expression on Th2 cells isolated from Nb-infected 4get mice longitudinally (d0, d6, d8, and d10) (Figure 3). CD4 cells of all samples were computationally gated, extracted, and concatenated before downsampling and dimensional reduction (UMAP) were applied based on Cxcr6, CD62L, Cxcr5, CD279 (PD-1), Ly6a/e, Itgb7, CD127, CD74, CD44, CD4, and CD3e expression, but not on IL-4eGFP expression to resolve similarities between IL-4eGFP⁺ and IL-4eGFP⁻ cells. Staining pattern for single markers (Figure 3—figure supplement 1A) and definition of gates aimed to resemble transcriptionally defined populations on protein level (Figure 3—figure supplement 1B) are given. Data of time points d6, d8, and d10 splitted by organ are shown either for all CD4⁺ cells or only for IL-4eGFP⁺CD4⁺ cells (Figure 3B). IL-4eGFP⁺CD4⁺ cells of the lung show a different distribution over the UMAP than LN cells, with the majority of cells falling into the L1 gate while the L2 population is less well defined by chosen markers. There is also additional heterogeneity with a number of cells also present in gates defined for MLN cells. As on transcriptional level, the distribution of MLN cells seems again comparable between medLN and MLN. Next, concatenated CD4⁺ T cells (Figure 3—figure supplement 1C) or only IL-4eGFP⁺CD4⁺ T cells (Figure 3C) were splitted based on organ and time after infection. In the lung, the L1 effector cluster increases with infection time in accordance with a developing effector response against Nb. Cells in the L2 gate were rather variable and there is no clear trend. The medLN response seems a bit faster than the MLN response but overall similar. Effector cells (in L1, MLN4, and MLN6 gate) moderately increase over time while population MLN2 is increased and MLN3 is strongly increased at d6 and d8 after Nb infection. This pattern is maintained at d10 for MLN but medLN population MLN1, MLN2, and MLN3 are all expressed to the same extend. It might therefore be that MLN3 cells are precursors of MLN2 and MLN1 as we see the proportional increase later on. That the increase in MLN1 and MLN2 signature cells at d10 is not seen in MLN might reflect a delayed development related to the delayed occurrence of worms in the intestine compared to the lung and medLN.

Figure 3 with 1 supplement see all

Download asset Open asset

Th2 subset distribution in *Nippostrongylus brasiliensis* (Nb) infection time course.

CD4-expressing cells from lung, medLN, and mesenteric lymph node (MLN) of IL-4eGFP reporter (4get) mice were stained for markers defined by single-cell sequencing at indicated days after Nb infection. (A) Selection of transcriptionally defined marker genes encoding for surface receptors. (B) Surface staining for selected markers on protein level was performed followed by flow cytometric analysis. CD4⁺ cells were computationally selected, down sampled and a UMAP of the data was generated. Shown are CD4 T cells of all organs after infection (days 6, 8, and 10 post Nb infection) concatenated into one UMAP and splitted by organ (upper row). The same is given for IL-4eGFP⁺CD4⁺ cells (lower row). Gates are aimed to resemble transcriptionally defined cluster on protein level. (C) Quantified distribution of IL-4eGFP⁺CD4⁺ cells in UMAP gates over the course of Nb infection. For (B, C), n = 3 mice per time point; error bars represent standard error (SE).

In summary, we established a surface staining strategy that allows to track Th2 subpopulations during Nb infection by flow cytometry.

Clonal relatedness of Th2 cells in distant organs confirms effective exchange of effector cells

The single-cell immune profiling approach allows for combined RNA expression profiling and TCR repertoire analysis, which made exploration of clonal relatedness between clusters and efficient distributed of clones across organs possible. This analysis also aimed to define potent Nb-specific TCRs based on the hypothesis that these might belong to successfully expanded clones associated with above defined effector clusters.

In line with previous results (Seidl et al., 2011), we find a broad TCR repertoire after Nb infection as the majority of distinct TCRs was only found in one cell. However, 28% of cells expressed a TCR found in at least two different cells (same CDR3 nucleotide sequence, the same variable and joining segments). The most abundant clone has 15 sequenced members in the samples analyzed which translates to about 8600 estimated members in total lung and MLN tissue. As we only analyze a small sample of the whole organs the calculation clearly underestimates the fraction of expanded T-cell clones in the population. A still substantial part of clones is found in both organs, which suggests an effective distribution between MLN and lung. In contrast, only two small clones were identical between the two analyzed mice implicating very few public clones (Figure 4A). The innate-like/NKT innate cluster (MLN6) and the MLN7 cluster hardly contained expanded clones suggesting limited TCR specificity-driven proliferation in these clusters (Figure 4B). Cells of the ‘bridge’ clusters (MLN5 and L2) contain substantially more expanded clones but less than the effector/resident memory like populations (MLN4 and L1), which in turn contain less expanded clones than the more homogeneous majority of MLN clusters (MLN1, MLN2, and MLN3). It might reflect that the ‘bridge’ clusters contain immi-/emmigrated cells from distant sites with less clonal overlap to the local population.

Figure 4

Download asset Open asset

Clonal relatedness of Th2 cells in mesenteric lymph node (MLN) and lung.

(A) UMAP of MLN and lung cells split by cells with unique T-cell receptors (TCRs) (nonexpanded), cells with the same TCRs found in more than one cell (expanded), cells with the same TCRs found in both organs and cells with the same TCRs found in both mice. Cells are colored by cluster. Schematic drawing roughly highlights how MLN and lung cells are separated on UMAP. (B) Stacked bar plot on the presence of expanded clones per cluster. Expansion level relates to overall presence in the dataset. (C) Fraction of cells in MLN, lung, and both organs in relation to clone expansion. Numbers above indicate proportion of expanded cells in total population. (D) Clonal relatedness between clusters. The stacked bar to the left gives the fraction of each cluster in the dataset as a reference (proliferating cluster was excluded). Stacked bar graphs to the right visualize for every expanded clone of a cluster where other members of a clone are found (cluster distribution). Numbers above bars represent the number of cells that each bar represents. Bars for clusters with only few expanded clones are not shown. Black horizontal lines separate the MLN and lung clusters in the bar graphs. Cluster of cells with a migratory signature is highlighted as ‘bridge cluster’. (E) Top 5 expanded clones by occurrence in MLN (left), lung (center), or in total dataset (right).

Next, we find that strongly expanded clones are effectively spread over organs (Figure 4C). The typical caveats of current single-cell technologies (sampling noise and limited sample size) do not allow to draw a similar conclusion for lowly expanded clones (<3 cells per organ). The determination of the clonal relatedness of clusters compared to the overall frequency of a cluster in the dataset again highlights effective distribution of effector Th2 cells between distant organs (Figure 4D). The clones of clusters that are most distant to cells of the other organ (L1, MLN1, MLN2, and MLN3) tend to expand more locally, represented by the higher percentage of related cells found in the same organ compared to the overall distribution of clusters. Directly compared to those clusters, the clones in ‘bridge’ clusters (L2 and MLN5) have a higher frequency of members in the other organ, especially apparent in the other part of the ‘bridge’ in each case. The effector/resident memory like cluster of the MLN (MLN4) also shows increased relatedness to the lung that contains a large number of effector cells in cluster L1. The finding that clusters visualized near the other organ also show enhanced TCR repertoire relatedness to that organ confirms significant migration between organs and that the UMAP efficiently displays real relatedness of clusters.

As a next step, we visualized the five most strongly expanded clones determined for each organ separately or we combined counts to get the most strongly expanded clones in the total dataset (Figure 4E). Members of such clones in the MLN tend to be preferentially found in the MLN1 cluster and to a lower extend in the neighboring Tfh-associated cluster (MLN2) and the IFN signature cluster (MLN3). They are hardly found in the lung proximal clusters (MLN4–MLN7). However, few members of four of the top expanded MLN clones are found in the lung and are therefore successfully spread across organs. Top expanded lung clones do not overlap with the top expanded MLN clones and preferentially show up in the big lung cluster in which effector cells are found (L1). In contrast to the MLN only one of the top expanded lung clones has members in the MLN, indicating that these clones successfully expanded locally but have limited capacity to spread to the MLN. The five most highly expanded clones in the whole dataset strongly overlap with the ones determined for the separate organs. This indicates that despite remarkable exchange between the distant organs, strong local expanders dominate the response and while more evenly distributed clones are present, they do not outnumber locally expanded ones in a combined analysis of MLN and lung cells.

In summary, there is substantial overlap of expanded clones between the MLN and lung during Nb infection, but rather locally restricted clones successfully expanded in an otherwise diverse pool of Th2 cells.

No general preference for specific TCR chain compositions

After analysis of single clones in the last part, where we found expansion but no obvious dominant clones, we determined if there are rather preferentially used TCR segments or segment combinations that could suggest potent effector function of TCRs against Nb. First, we included only one representative member per clone and compared if the same combination of TCRα and TCRβ chain segments is shared between the top 50 most frequently used segment combinations in both organs, both analyzed mice and in nonexpanded versus expanded clones (Figure 5A). For the nonexpanded clones there was hardly any overlap seen between mice or organs, only two of the fifty combinations were found in three of the four analyzed organs (MLN and lung of two mice). For the expanded clones, there was limited overlap in combined segment usage between organs of one mouse but not the other. Combinations of TCRα- or TCRβ variable with joining segments and TCRα with TCRβ variable segments also revealed limited overlap in the top used combinations. Trbv1 was a recurrently used TCRβ variable segment present in frequently used combinations (Figure 5—figure supplement 1A). Similarly, for single segments there was no obvious preferential usage in expanded clones compared to nonexpanded ones. Again, Trbv1 was one of few constituents that was moderately increased in expanded versus nonexpanded clones (Figure 5—figure supplement 1B,E). In addition, there was no striking difference observed in total CDR3 amino acid length/length distribution that could be indicative for changes in specificity (Davis et al., 1998; Rock et al., 1994) between expanded and nonexpanded clones (Figure 5—figure supplement 2A). In a finer grained analysis of single TCRα and TCRβ family members, there was also no change in CDR3 length or length distribution (Figure 5—figure supplement 2B,C). The general TCRα or TCRβ CDR3 length in MLN and lung of the Nb-infected mice is also not altered compared to naive T cells of the peripheral blood (Figure 5B).

Figure 5 with 2 supplements see all

Download asset Open asset

Expanded CDR3 motifs in Th2 cells of *Nippostrongylus brasiliensis* (Nb)-infected mice.

T-cell receptor (TCR) repertoire analysis of mesenteric lymph node (MLN) and lung Th2 cells at day 10 post Nb infection. (A) TCR segment combination analysis for overlap of the hundred most commonly used TCR segment combinations (V, J, and C region for TCRα; V, D, J, and C region for TCRβ) between different mice and organs. Analysis is performed separately for nonexpanded and expanded clones. (B) Amino acid sequence length of TCRα and TCRβ CDR3 regions. We compare CDR3 regions from peripheral blood T cells of naive wild-type C57BL/6 mice (naive) with CDR3 regions of Th2 cells from MLN and lung of Nb-infected mice. (C) Most abundant CDR3 amino acid sequences in cells of dataset presented as percent of each sample. (D) Combinations of TCRα-related CDR3 motifs (x-axis) with indicated CDR3 motifs of TCRβ or a second TCRα chain (y-axis) that occurred more than once were counted. Then TCRα chains were ranked by abundance in these combinations (left). We also present highly expanded TCRα_CDR3 sequences found in combination with various TCRβ or TCRα chains (right). At the bottom, cells that express the corresponding CDR3 sequences from the x-axis are highlighted on top of UMAP representation of the dataset. We also indicate if a CDR3 sequence is associated with the top expanded clones (Figure 4E) below UMAP plots.

Figure 5—source data 1 R script for analyzes of expanded CDR3 motifs from TCR single cell sequencing data.: https://cdn.elifesciences.org/articles/74183/elife-74183-fig5-data1-v2.txt
Download elife-74183-fig5-data1-v2.txt

In conclusion, expanded clones in the Th2 effector population show no evidence for preferential usage of particular TCRα or TCRβ chains or chain combinations.

Definition of abundant CDR3 motifs in Th2 cell of Nb-infected mice

As we only found moderately expanded clones and minor signs for preferentially used TCR chains or chain combinations in our data, we performed a finer grained analysis. T-cell antigen specificity and affinity are mainly confined by variable regions of the TCR (mainly CDR3) (Rock et al., 1994). Therefore, we investigated if specific CDR3 motifs that might indicate potent Nb-specific effector function are repeatedly found in MLN and lung after Nb infection.

We chose the 10 most abundant CDR3 sequences on amino acid level that potentially include motifs relevant for anti-helminth immunity and highlight their abundance in different organs and mice (Figure 5C). TCR data of naive peripheral blood T cells served as a reference to identify germline-associated CDR3 regions. The most abundant sequence CVVGDRGSALGRLHF was found in all samples (Figure 5C), including the peripheral blood and represents the CDR3 motif found in the invariant TCRα chain (Vα14-Jα18) of NKT cells (Lantz and Bendelac, 1994). This chain is coexpressed with a variety of different TCRβ chains (not shown) and cells that express such TCRs are primarily found in the innate-like/NKT cluster (MLN6) (Figure 1E). As this cluster also contains some cells with expression of TCRγ chain segments in addition to functional TCRα and TCRβ chains we were able to compare if expression of the invariant TCRα is correlated to TCRγ expression. Indeed, on the one hand, 38% of cells that expressed any TCRγ constant or variable chain segment also expressed the invariant TCRα chain and on the other hand, 62% of cells with the invariant TCRα chain expressed any TCRγ constant or variable chain segment (detection of TCRγ and the invariant TCRα chain in the same cell: correlation 0.48; p < 10⁻⁶). This might suggest a close relatedness of IL-4-expressing γδ T cells with IL-4-expressing NKT cells in a way that cells unsuccessful or not yet successful to generate a functional γδ TCR preferentially develop into αβNKT cells. Alternatively, NKT cells could induce low level of TCRγ gene expression for other, unknown reasons.

Of the remaining nine most abundant CDR3 motifs of TCRα or TCRβ chains, seven are not found in the naive peripheral blood sample (Figure 5C), which implies an increased probability for them to represent specificity for Nb antigens. Furthermore, only one of these motifs (CAIDPSGSWQLIF) is expressed in both analyzed organs and both mice, which implies that it could be a preferentially selected motif during Nb infection.

We next determined CDR3 motifs that are part of abundant motif combinations (Figure 5D, left panel). As expected, these overlap with the most expanded clones (Figure 4E) as cells of an expanded clone always use the same chain combination. Only one of the five TCRα CDR3 motifs (CAAEAGTGGYKVVF) was associated with expansion in more than one clone (two clones with same TCRα CDR3 motif but different TCRβ CDR3 motifs). In addition, all five depicted TCRα CDR3 motifs present in abundant pairings are also present in unique pairings with other TCRβ CDR3 motifs. This implies that these motifs are not restricted to an exact TCRα/TCRβ combination or a single clone to be recruited to the Th2 compartment.

As others described (He et al., 2002; Padovan et al., 1993; Padovan et al., 1995) we find T-cell clones with expression of two TCRα/TCRβ chain-encoding genes. At least in highly abundant combinations it is unlikely that these are technical artifacts due to contamination with RNA from another cell during library preparation. The clone with the most frequently found combination of CDR3 motifs (clone MLN1) expresses one TCRβ and two TCRα chains, both on average with similar umi counts. Whether both TCRα chains are successfully translated is not known. Of the five depicted TCRα CDR3 motifs, often present in successful CDR3 combinations, four were found in some cells that expressed more than two TCR chains but this frequency is in the expected range for T cells (Alam and Gascoigne, 1998; Balomenos et al., 1995; Davodeau et al., 1995; Dupic et al., 2019).

In addition to CDR3 motifs found frequently in combinations (expanded clones), we also find abundant CDR3 motifs combined with various other unique TCR chains (present in several nonexpanded clones) (Figure 5D, right panel). These could include CDR3 motifs that provide anti-Nb specificity but failed to induce substantial expansion or accumulation of Th2 cells expressing such TCRs.

In line with a slightly preferential usage of the Trbv1 gene in expanded compared to nonexpanded Th2 cells we find that three of the expanded TCRα CDR3 motifs (CAIDPSGSWQLIF, CAIDSSGSWQLIF, and CAASDTNTGKLTF) are preferentially coexpressed with Trbv1, which suggest that Trbv1 might be relevant for the immune response against Nb.

Functional in vivo characterization of effector-associated TCRs

As the TCR response to Nb is highly heterogeneous and as we did not find obvious dominant clones, TCR chain combinations or CDR3 motifs shared between clones, we selected TCRs of three moderately expanded clones with members in the IL-5/IL-13 expressing region of effector-associated cluster L1 (Figure 6A, Figure 6—figure supplement 1A). We hypothesized that TCRs from effector cluster-associated clones might contribute to the development of potent effector cells. The TCRs were retrovirally transduced in hematopoietic stem cells isolated from Rag1-ko mice (not able to rearrange B- and T-cell receptors) that were then transferred into irradiated recipients to generate TCR retrogenic mice (Holst et al., 2006). Eight weeks after transfer, we verified the presence of T cells that contain the introduced TCR α- and β-chain constructs in MLN and spleen of the retrogenic RAG1ko recipient mice (Figure 6—figure supplement 1B). T cells with the transgenic TCR chains (α and β) were used for transfer experiments into mice subsequently infected with Nb (Figure 6B). At day 9 post infection expansion of T cells with the transgenic TCRs was measured by flow cytometry as a readout for Nb reactivity. As controls, mice were either left uninfected or they were infected with another parasitic worm (Heligmosomoides polygyrus) to control for an unspecific expansion in response to worms. We observed a moderate but Nb-specific in vivo expansion for one of the three clones (NB-T2) highlighting that the chosen approach is suitable for the selection and reexpression of Nb-specific TCRs (Figure 6C,D). To our knowledge, this is the first identification of an in vivo validated Nb-specific TCR which might be used in the future to characterize the T-cell response to Nb with a given TCR specificity.

Figure 6 with 1 supplement see all

Download asset Open asset

Expansion of Nb-T2 T-cell receptor (TCR) cells in mesenteric lymph node (MLN) after *Nippostrongylus brasiliensis* (Nb) infection.

Hematopoietic stem cells were transduced with retroviral vectors for expression of potential Nb-specific clones (Nb-T1, Nb-T2, and Nb-T3) and used to generate retrogenic mice. (A) UMAP overlay shows distribution of TCRs selected for retrogenic expression. V gene segments, CDR3 sequence, and distribution between organs for selected lung-expanded TCR clones (Nb-T1, Nb-T2, and Nb-T3). (B) Eight weeks after reconstitution of irradiated Rag1KO mice with TCR chain and fluorescent marker-encoding retrovirus containing retrogenic Rag1 stem cells, T cells were harvested, transferred to Ly5.1 B6 wild-type mice which were then infected with Nb or *Heligmosomoides polygyrus* (Hp) as control. (C) Plots show the percentage of TCR transgenic cells (GFP⁺, Thy1.1⁺) among live CD4⁺ cells from constructs Nb-T1, Nb-T2, and Nb-T3 in the MLN of recipient mice at day 9 post infection. (D) Percentage and total number of Nb-T2 cells in MLN. Quantification is based on six independent experiments. Statistical significance determined by Mann–Whitney U test; **p < 0.001.

Discussion

Th2 heterogeneity, organ crosstalk, and tissue-specific immunity are increasingly appreciated (Schoettler et al., 2019; Szabo et al., 2019a; Szabo et al., 2019b). Here, we applied combined transcriptome and TCR clonotype analysis on Th2 cells across organs upon Nb infection. We identified lung- and MLN-specific gene signatures as well as subpopulations with shared migration and effector/resident memory profiles between organs. We find that expression of tissue damage-associated cytokine coding genes Il13 and Il5 is restricted to the effector/resident memory populations in lung and MLN. Interestingly, these clusters also contain transcriptionally similar cells that express Il10 but widely lack expression of the Treg marker-encoding gene Foxp3. Similar cells have been described in the skin at the infection site after repeated Schistosoma mansoni cercaria infection where these cells have immunosuppressive functions (Sanin et al., 2015). Furthermore, effector/resident memory like cells in the MLN are not homogeneous and are found in two clusters of which one is an innate-like/NKT cluster. Interestingly, the NKT population in this cluster not only coexpresses the invariant NKT cell-associated invariant TCRα chain (Vα14-Jα18) together with a highly diverse repertoire of TCRβ chains but also transcripts for TCRγ chains, which implies shared developmental pathways of NKT and γδ T cells that both tend to express restricted receptor repertoires. The cluster also contains cells that express Cxcr3, encoding a typical Th1 marker. CXCR3 has been noted in a small fraction of Th2 cells (Kim et al., 2001) but was not associated with IL-4 producing NKT or γδ T cells before. These findings reveal a heterogeneous effector/resident memory pool in the Th2 population.

When we searched for activation signatures, we found a population of cells with an IFN response signature (MLN3) present in the MLN. Murine in vitro studies imply that IFN signaling is needed for proper Th2 differentiation (Wensky et al., 2001). Therefore, Th2 cells with an IFN response signature probably reflect cells that undergo priming or differentiation. RNA velocity analysis further defines them as a differentiation endpoint. Of note, the differentiation state not necessarily correlates with time in the Th2 compartment and cells could still be recently selected or undergo plasticity. Nevertheless, RNA velocity also indicates a fraction of proliferating cells that shows low expression of Th2 signature genes and that is likely in a stage where differentiation can be determined. This supports that proliferation and differentiation during Th2 development are linked as suggested before (Proserpio et al., 2016) rather than being independent and orthogonal processes as alternatively hypothesized (Duffy and Hodgkin, 2012). In addition to potential developmental paths, RNA velocity analysis suggests that development is faster and has a stricter directionality in the lung compared to the MLN, consistent with the view that the majority of Th2 cells in the MLN are less terminally differentiated.

Cells in the migratory clusters of both organs show a weaker organ-specific separation after dimensional reduction but rather form a ‘bridge’ on the UMAP that suggests effective exchange between organs. In line, TCR analysis of our Th2 cells revealed effective exchange of expanded clones between organs. However, the most expanded clones in the one organ were not the most expanded in the other organ. This might relate to different immunological preferences in different compartments or to the different larval stages in which Nb is present in lung and MLN. Of note, comparison of human bulk TCR repertoires between the lung and its draining lymph node also showed a higher intraorgan TCR repertoire overlap than between organs. This was interpreted to mean that the T cells originate from different precursor pools and recognize distinct antigens (Schoettler et al., 2019). Our data clearly refine that there is effective spread of Th2 effector T cells from the same pool of cells even across distant organs.

Another factor that likely protects from an overshooting Th2 response and could be subject to modulation by the worm is the expression of TNFR2 as impaired TNFR2 signaling leads to augmented Th2 responses (Li et al., 2017). We observed transcripts for TNFR2 preferentially in the lung compared to the MLN after Nb infection, which has not yet been described to our knowledge.

In summary, combined transcriptome and TCR clonotype analysis at single-cell resolution provides information about Th2 heterogeneity across organs and reveals relatedness of IL-4 producing NKT cells to γδ T cells. RNA velocity combined with knowledge from published data appears compatible with a model, in which poorly differentiated proliferating Th2 cells are at a decision point in their development, with IFN signaling being involved in diversification and differentiation of the Th2 compartment. Despite efficient exchange of expanded Th2 clones between distant organs, the most abundant clones seem to expand locally. As a proof of principle we functionally validated an Nb-specific TCR in vivo. The setup allows further functional characterization of expanded TCR clonotypes and will help to investigate Th2 cell differentiation, plasticity, and memory formation in response to natural helminth pathogens.

Materials and methods

Key resources table

Reagent type (species) or resource	Designation	Source or reference	Identifiers	Additional information
Mus musculus, C57BL/6 background, female	4get mice	Mohrs et al., 2001		IL-4eGFP reporter mouse
Antibody	anti-CD16/32, unconjugated, clone 2.4G2, Rat Monoclonal Antibody	BioXCell		(1:200)
Antibody	anti-CD4, PerCP/Cyanine5.5, clone RM4-5, Rat Monoclonal Antibody	Invitrogen	45-0042-82	(1:200)
Antibody	anti-CD4, BUV737, clone RM4-5, Rat Monoclonal Antibody	BD	Cat#: 612844	(1:100)
Antibody	anti-CXCR6 PerCP-eFluor 710, clone DANID2, Rat Monoclonal Antibody	Invitrogen	Cat#: 46-9186-82	(1:100)
Antibody	anti-CD62L, eFluor 450, clone MEL-14, Rat Monoclonal Antibody	Invitrogen	Cat#: 48-0621-82	(1:200)
Antibody	anti-CD127, BV711, clone SB/199, Rat Monoclonal Antibody	BD	Cat#: 565490	(1:100)
Antibody	anti-Ly-6A/E (Sca-1), BV785, clone D7, Rat Monoclonal Antibody	BioLegend	Cat#: 108139	(1:100)
Antibody	anti-CD3e, BV605, clone 145-2 C11, Armenian Hamster Monoclonal Antibody	BD	Cat#: 563004	(1:100)
Antibody	anti-CD279 (PD-1), Super Bright 780, clone J43, Armenian Hamster Monoclonal Antibody	Invitrogen	Cat#: 78-9985-82	(1:400)
Antibody	anti-CD74, Alexa Fluor 647, clone In1/CD74, Rat Monoclonal Antibody	BioLegend	Cat#: 151004	(1:100)
Antibody	anti-CD44, BUV395, clone IM7, Rat Monoclonal Antibody	BD	Cat#: 740215	(1:100)
Antibody	anti-CD185 (CXCR5), PE/Cyanine7, clone L138D7, Rat Monoclonal Antibody	BioLegend	Cat#: 145516	(1:100)
Antibody	anti-integrin β7, PE, clone FIB27, Rat Monoclonal Antibody	BioLegend	Cat#: 121006	(1:100)
Antibody	TotalSeq-C0301 is a mix of anti-mouse MHC class I and anti-mouse Ly-5 antibodies, clone M1/42 and 30-F11, Rat Monoclonal Antibodies; Hashtag 1 Antibody; Barcode Sequence ACCCACCAGTAAGAC; mouse 1 med LN	BioLegend	Cat#: 155861	(3 μl for 6 × 10⁶ cells)
Antibody	TotalSeq-C0302 is a mix of anti-mouse MHC class I and anti-mouse Ly-5 antibodies, clone M1/42 and 30-F11, Rat Monoclonal Antibodies; Hashtag 2 Antibody; Barcode Sequence GGTCGAGAGCATTCA; mouse 2 med LN	BioLegend	Cat#: 155863	(3 μl for 6 × 10⁶ cells)
Antibody	TotalSeq-C0303 is a mix of anti-mouse MHC class I and anti-mouse Ly-5 antibodies, clone M1/42 and 30-F11, Rat Monoclonal Antibodies; Hashtag 3 Antibody; Barcode Sequence CTTGCCGCATGTCAT; mouse 1 MLN	BioLegend	Cat#: 155865	(3 μl for 6 × 10⁶ cells)
Antibody	TotalSeq-C0304 is a mix of anti-mouse MHC class I and anti-mouse Ly-5 antibodies, clone M1/42 and 30-F11, Rat Monoclonal Antibodies; Hashtag 4 Antibody; Barcode Sequence AAAGCATTCTTCACG; mouse 2 MLN	BioLegend	Cat#: 155867	(3 μl for 6 × 10⁶ cells)
Other	LIVE/DEAD Fixable Lime (506) Viability Kit	Invitrogen	Cat#: L34990	(1:1000)
Software, algorithm	R	R Development Core Team, 2021	Versions 3.5.1, 4.0.0
Software, algorithm	Cell Ranger	Zheng et al., 2017	Versions 2.1.1, 3.0.1, 6.1.2
Software, algorithm	Seurat	Stuart et al., 2019	Versions 3.1.1, 4.1.0
Software, algorithm	kallisto/bustools	Bray et al., 2016	Versions 0.46.2/0.40.0
Software, algorithm	scVelo	Bergen et al., 2020	Version 0.2.3
Software, algorithm	monocle	Trapnell et al., 2014	Version 2.16.0
Software, algorithm	CytoTRACE	Trapnell et al., 2014	Version 0.3.3
Software, algorithm	FlowSOM	Van Gassen et al., 2015	Version 1.20.0
Software, algorithm	CytoTree	Dai et al., 2021	Version 1.3.0

Mice and N. brasiliensis infection

Request a detailed protocol

IL-4eGFP reporter (4get) mice (Mohrs et al., 2001) aged 8–16 weeks were infected with third-stage larvae (L3). Larvae were washed in sterile 0.9% saline (37°C) and 500 organisms were injected subcutaneously (s.c.) into mice 10 days before analysis (if not indicated otherwise). To avoid bacterial infections mice received antibiotics-containing drinking water (2 g/l neomycin sulfate, 100 mg/l polymyxin B sulfate; Sigma-Aldrich, St. Louis, MO) for the first 5 days after infection. Mice were kept under specific pathogen-free conditions and were maintained in the Franz-Penzoldt Center in Erlangen. All experiments were performed in accordance with German animal protection law and European Union guidelines 86/809 and were approved by the Federal Government of Lower Franconia.

Single-cell RNA and TCR sequencing

Request a detailed protocol

For each of the two performed single-cell sequencing experiments two IL-4eGFP reporter (4get) mice (Mohrs et al., 2001) at the age of 13 (MLN, lung comparison) or 10 weeks (medLN, MLN comparison) were infected with Nb. At day 10 after Nb infection lungs and lymph nodes (medLN/MLN) of IL-4eGFP reporter (4get) mice were harvested. Lungs were perfused with phosphate buffered saline, cut into small pieces and digested with the commercial ‘Lung Dissociation kit’ (Miltenyi, Bergisch Gladbach, GER) according to the manufacturer’s instructions. Digested lungs and complete LNs were gently mashed through a 100 μm cell strainer. Next, a 40% percoll purification was applied on lung cells and erythrocytes were lysed with ACK-buffer (0.15 M NH₄Cl, 1 mM KHO₃, 0.1 mM Na₂EDTA). Then samples of both organs were incubated with Fc-receptor blocking antibody (anti-CD16/32, clone 2.4G2, BioXCell, West Lebanon, NH) and stained with anti-CD4-Percp-Cy5.5 antibody (clone: RM4-5). For the medLN/MLN comparison cells of each organ and each mouse were tagged with distinguishable hashtag abs (0.5 µg ab per 2 × 10⁶ cells) added to the fluorescence antibody staining. IL-4eGFP⁺CD4⁺ cells were sorted and for each sample, 5000 cells (in case of MLN and lung comparison) or 6250 cells (in case of medLN and MLN comparison) were sorted. For MLN and lung comparison cells were subjected to 10× Chromium Single Cell 5′ Solution v2 library preparation using the TCR-specific VDJ library kit according to the manufacturer’s instructions (10× Genomics, Pleasanton, CA). Gene expression libraries were sequenced on an Illumina HiSeq 2500 sequencer using the recommended read lengths for 10× Chromium 5’ v2 chemistry to a depth of at least 30000 reads per cell. VDJ libraries were sequenced as paired 150 bp reads to a depth of at least 30,000 reads per cell. For the medLN and MLN comparison, we followed the chromiumNextGEMSingleCell5_v2_Cell SurfaceProtein guide (RevA). IL-4eGFP⁺CD4⁺ cells of different mice and organs were labeled with hashtag antibodies and multiplexed for sequencing. The same multiplexing was performed for IL-4eGFP⁻CD4⁺ cells in parallel.

Computational analysis

Request a detailed protocol

For the MLN and lung sequencing run we used 10× Genomics Cell Ranger to demultiplex sequencing reads, convert them to FASTQ format with mkfastq (Cell Ranger 2.1.1), align them to the murine genome (mm10 v3.0.0) and obtain TCR VDJ clonotypes, consensus sequences, contigs, and CDR3 regions (Cell Ranger 3.0.1). TCR-associated genes (VDJ and constant region genes for α, β, γ, and δ chains) were excluded but kept as metadata to avoid clustering by TCR genes. To be included, cells needed to be defined as such by Cell Ranger and to have >500 UMIs, >500 genes detected per cell, <7% mitochondrial reads, and a novelty >0.8 (log10 of gene number divided by log10 of UMIs). Data normalization, differential expression analysis, clustering (based on top 2000 highly variable genes), and dimensional reduction (UMAP based on top 15 principal components) were performed in Seurat (version 3.1.1) (Stuart et al., 2019) under R (version 3.5.1). For the medLN and MLN sequencing run Cell Ranger 6.1.2 was used for initial aggregation and alignment. Demultiplexing of hashtagged samples and integration analysis was performed in Seurat 4.1.0 under R 4.0.0. Other analysis were performed in the R and Seurat versions used for MLN and lung sequencing if not indicated differently. QC was also performed as for MLN and lung but TCR expression was no criteria. Gene set scores for each cell were calculated in Seurat as published before (Tirosh et al., 2016). Gene sets were taken or generated from published data: resident memory and circulating memory (Rahimi et al., 2020), IL-33 signature (Morimoto et al., 2018), proliferation signature genes for G2M- and S-phase to regress out cell cycle in scVelo (Macosko et al., 2015), other sets were from the ‘Molecular Signatures Database’ (Subramanian et al., 2005). TCR info was added as metadata to Seurat for combined clonotype and RNA-profile analysis. For RNA velocity, sequencing reads were aligned with kallisto/bustools (version 0.46.2/0.40.0) (Bray et al., 2016; Melsted et al., 2021) to a genome reference with unspliced and spliced RNA variants included (version GRCm39). Obtained information was used as input for scVelo (version 0.2.3) (Bergen et al., 2020) under python (version 3.8.5). UMAP information from Seurat was transferred to scVelo for consistency and dynamic and static velocity models with additional paga analysis were calculated with the scVelo package. Additionally a trajectory analysis with monocle (version 2.16.0) (Trapnell et al., 2014) was performed and developmental potential was additionally predicted with CytoTRACE (version 0.3.3) (Gulati et al., 2020). Usage of TCR chains and TCR chain combinations was calculated under R with custom scripts (Figure 5—source data 1). For TCR/CDR3 analysis, we used the Cell Ranger output and followed a recently developed workflow (according to the "CellaRepertorium" R package) with minor modifications. Contigs that missed a sanity check were excluded (needed to be productive, full length, high confidence, supported by >1 UMI, CDR3 length >5 amino acids). Similar CDR3 sequences were not combined (not assuming similar specificity for similar sequences) to maintain higher accuracy. We kept all TCR chains of T-cell clones with two TCRα/TCRβ chain-encoding genes expressed for the same reason. Data are available via GEO (GSE181342) and the 10× Genomics TCR reference dataset via the 10xGenomics website: PBMCs from C57BL/6 mice (v1, 150 × 150), Single Cell Immune Profiling Dataset by Cell Ranger 3.0.0, 10× Genomics (November 19, 2018). Flow cytometric data were pre-gated for living, CD4⁺ singlets and used as input for clustering with FlowSOM (1.20.0) (Van Gassen et al., 2015) using the CytoTree (1.3.0) Toolkit package (Dai et al., 2021) with transformMethod "autoLgcl" and scaling of data to give every channel the same weight. Data were downsampled to about 2.5 × 10⁵ cells preserving FlowSOM defined cluster proportion for UMAP calculation (using fluorophore channel excluding eGFP IL-4 reporter and live/dead discriminator). Then data including UMAP information were exported to FCS format for further analysis.

Generation of TCR retrogenic (Rg) mice

Request a detailed protocol

For analysis of Nb-specific TCR activation, TCR retrogenic mice were generated by transferring retrovirally transduced bone marrow stem cells into irradiated recipient mice as previously described (Holst et al., 2006). In brief, we generated TCR retrogenic mice for expression of the following TCRs: Nb-T1, Nb-T2, and Nb-T3 (Figure 6—figure supplement 1A). For generation of retroviral particles, using the ecotropic Phoenix E cell line and a standard Calcium Phosphate transfection protocol, the target genes were cloned into pMXpie or MSCV-Thy1.1 for expression of the TCR α- or β-chains, respectively. Forty-eight hours before harvesting the bone marrow, donor mice (RAG1KO) were injected with sterile 5-fluorouracil solution i.p. One day before infection, harvested bone marrow cells were seeded into 24-well plates at a concentration of 2–3 × 10⁶ cells in 1 ml medium supplemented with IL-3 (20 ng/ml), IL-6 (50 ng/ml), and SCF (100 ng/ml). Retroviral infection of dividing cells was performed using a spinoculation protocol with RetroNectin and two infection steps. Recipient mice (Rag1KO) were sublethally irradiated (500 rad), followed by reconstitution with 0.5–2 × 10⁶ bone marrow stem cells i.v. Mice were provided with antibiotics-containing drinking water for 8 weeks. After reconstitution (8 weeks) the development of TCR transgenic cells was evaluated by flow cytometry of different organs and cells were further analyzed in a transfer experiment.

Data availability

Single-cell RNA sequencing data are available via GEO with the ID GSE181342.

The following data sets were generated

(2021) NCBI Gene Expression Omnibus
ID GSE181342. Th2 single-cell heterogeneity and clonal interorgan distribution in helminth-infected mice.

https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE181342

References

1. Alam SM
2. Gascoigne NR
(1998)
Posttranslational regulation of TCR valpha allelic exclusion during T cell differentiation

Journal of Immunology 160:3883–3890.
- Google Scholar
(2006) IL-4R signaling is required to induce IL-10 for the establishment of T(h)2 dominance
International Immunology 18:1421–1431.
https://doi.org/10.1093/intimm/dxl075
- PubMed
- Google Scholar
(2016) T follicular helper cell plasticity shapes pathogenic T helper 2 cell-mediated immunity to inhaled house dust mite
Immunity 44:259–273.
https://doi.org/10.1016/j.immuni.2015.11.017
- PubMed
- Google Scholar
(1995)
Incomplete T cell receptor V beta allelic exclusion and dual V beta-expressing cells

Journal of Immunology 155:3308–3312.
- PubMed
- Google Scholar
1. Bergen V
2. Lange M
3. Peidli S
4. Wolf FA
5. Theis FJ
(2020) Generalizing RNA velocity to transient cell states through dynamical modeling
Nature Biotechnology 38:1408–1414.
https://doi.org/10.1038/s41587-020-0591-3
- PubMed
- Google Scholar
(2016) Near-optimal probabilistic RNA-seq quantification
Nature Biotechnology 34:525–527.
https://doi.org/10.1038/nbt.3519
- PubMed
- Google Scholar
1. Breitfeld D
2. Ohl L
3. Kremmer E
4. Ellwart J
5. Sallusto F
6. Lipp M
7. Förster R
(2000) Follicular B helper T cells express CXC chemokine receptor 5, localize to B cell follicles, and support immunoglobulin production
The Journal of Experimental Medicine 192:1545–1552.
https://doi.org/10.1084/jem.192.11.1545
- PubMed
- Google Scholar
(2003) Animal model of nippostrongylus brasiliensis and heligmosomoides polygyrus
Current Protocols in Immunology Chapter 19:12.
https://doi.org/10.1002/0471142735.im1912s55
- PubMed
- Google Scholar
1. Campbell DJ
2. Butcher EC
(2002) Rapid acquisition of tissue-specific homing phenotypes by CD4(+) T cells activated in cutaneous or mucosal lymphoid tissues
The Journal of Experimental Medicine 195:135–141.
https://doi.org/10.1084/jem.20011502
- PubMed
- Google Scholar
1. Connor LM
2. Tang S-C
3. Cognard E
4. Ochiai S
5. Hilligan KL
6. Old SI
7. Pellefigues C
8. White RF
9. Patel D
10. Smith AAT
11. Eccles DA
12. Lamiable O
13. McConnell MJ
14. Ronchese F
(2017) Th2 responses are primed by skin dendritic cells with distinct transcriptional profiles
The Journal of Experimental Medicine 214:125–142.
https://doi.org/10.1084/jem.20160470
- PubMed
- Google Scholar
1. Dai Y
2. Xu A
3. Li J
4. Wu L
5. Yu S
6. Chen J
7. Zhao W
8. Sun X-J
9. Huang J
(2021) CytoTree: an R/bioconductor package for analysis and visualization of flow and mass cytometry data
BMC Bioinformatics 22:138.
https://doi.org/10.1186/s12859-021-04054-2
- PubMed
- Google Scholar
1. Davis MM
2. Boniface JJ
3. Reich Z
4. Lyons D
5. Hampl J
6. Arden B
7. Chien Y
(1998) Ligand recognition by alpha beta T cell receptors
Annual Review of Immunology 16:523–544.
https://doi.org/10.1146/annurev.immunol.16.1.523
- PubMed
- Google Scholar
1. Davodeau F
2. Peyrat MA
3. Romagné F
4. Necker A
5. Hallet MM
6. Vié H
7. Bonneville M
(1995) Dual T cell receptor beta chain expression on human T lymphocytes
The Journal of Experimental Medicine 181:1391–1398.
https://doi.org/10.1084/jem.181.4.1391
- PubMed
- Google Scholar
1. Duffy KR
2. Hodgkin PD
(2012) Intracellular competition for fates in the immune system
Trends in Cell Biology 22:457–464.
https://doi.org/10.1016/j.tcb.2012.05.004
- PubMed
- Google Scholar
1. Dupic T
2. Marcou Q
3. Walczak AM
4. Mora T
(2019) Genesis of the alphabeta T-cell receptor
PLOS Computational Biology 15:e1006874.
https://doi.org/10.1371/journal.pcbi.1006874
- Google Scholar
1. Gett AV
2. Hodgkin PD
(1998) Cell division regulates the T cell cytokine repertoire, revealing a mechanism underlying immune class regulation
PNAS 95:9488–9493.
https://doi.org/10.1073/pnas.95.16.9488
- PubMed
- Google Scholar
(2009) T follicular helper cells differentiate from th2 cells in response to helminth antigens
The Journal of Experimental Medicine 206:991–999.
https://doi.org/10.1084/jem.20090303
- PubMed
- Google Scholar
1. Gowthaman U
2. Chen JS
3. Zhang B
4. Flynn WF
5. Lu Y
6. Song W
7. Joseph J
8. Gertie JA
9. Xu L
10. Collet MA
11. Grassmann JDS
12. Simoneau T
13. Chiang D
14. Berin MC
15. Craft JE
16. Weinstein JS
17. Williams A
18. Eisenbarth SC
(2019) Identification of a T follicular helper cell subset that drives anaphylactic ige
Science 365:eaaw6433.
https://doi.org/10.1126/science.aaw6433
- PubMed
- Google Scholar
1. Gulati GS
2. Sikandar SS
3. Wesche DJ
4. Manjunath A
5. Bharadwaj A
6. Berger MJ
7. Ilagan F
8. Kuo AH
9. Hsieh RW
10. Cai S
11. Zabala M
12. Scheeren FA
13. Lobo NA
14. Qian D
15. Yu FB
16. Dirbas FM
17. Clarke MF
18. Newman AM
(2020) Single-cell transcriptional diversity is a hallmark of developmental potential
Science 367:405–411.
https://doi.org/10.1126/science.aax0249
- PubMed
- Google Scholar
1. Guo L
2. Wei G
3. Zhu J
4. Liao W
5. Leonard WJ
6. Zhao K
7. Paul W
(2009) IL-1 family members and STAT activators induce cytokine production by th2, th17, and th1 cells
PNAS 106:13463–13468.
https://doi.org/10.1073/pnas.0906988106
- PubMed
- Google Scholar
1. He X
2. Janeway CA
3. Levine M
4. Robinson E
5. Preston-Hurlburt P
6. Viret C
7. Bottomly K
(2002) Dual receptor T cells extend the immune repertoire for foreign antigens
Nature Immunology 3:127–134.
https://doi.org/10.1038/ni751
- PubMed
- Google Scholar
1. Healey DCC
2. Cephus JY
3. Barone SM
4. Chowdhury NU
5. Dahunsi DO
6. Madden MZ
7. Ye X
8. Yu X
9. Olszewski K
10. Young K
11. Gerriets VA
12. Siska PJ
13. Dworski R
14. Hemler J
15. Locasale JW
16. Poyurovsky MV
17. Peebles RS
18. Irish JM
19. Newcomb DC
20. Rathmell JC
(2021) Targeting in vivo metabolic vulnerabilities of th2 and th17 cells reduces airway inflammation
Journal of Immunology 206:1127–1139.
https://doi.org/10.4049/jimmunol.2001029
- PubMed
- Google Scholar
(2006) Generation of T-cell receptor retrogenic mice
Nature Protocols 1:406–417.
https://doi.org/10.1038/nprot.2006.61
- PubMed
- Google Scholar
1. Hung LY
2. Lewkowich IP
3. Dawson LA
4. Downey J
5. Yang Y
6. Smith DE
7. Herbert DR
(2013) IL-33 drives biphasic IL-13 production for noncanonical type 2 immunity against hookworms
PNAS 110:282–287.
https://doi.org/10.1073/pnas.1206587110
- PubMed
- Google Scholar
1. Kim CH
2. Rott L
3. Kunkel EJ
4. Genovese MC
5. Andrew DP
6. Wu L
7. Butcher EC
(2001) Rules of chemokine receptor association with T cell polarization in vivo
The Journal of Clinical Investigation 108:1331–1339.
https://doi.org/10.1172/JCI13543
- PubMed
- Google Scholar
1. King IL
2. Mohrs M
(2009) IL-4-producing CD4+ T cells in reactive lymph nodes during helminth infection are T follicular helper cells
The Journal of Experimental Medicine 206:1001–1007.
https://doi.org/10.1084/jem.20090313
- PubMed
- Google Scholar
(1997) Expression of ST2, an interleukin-1 receptor homologue, is induced by proinflammatory stimuli
Biochemical and Biophysical Research Communications 235:474–478.
https://doi.org/10.1006/bbrc.1997.6810
- PubMed
- Google Scholar
1. La Manno G
2. Soldatov R
3. Zeisel A
4. Braun E
5. Hochgerner H
6. Petukhov V
7. Lidschreiber K
8. Kastriti ME
9. Lönnerberg P
10. Furlan A
11. Fan J
12. Borm LE
13. Liu Z
14. van Bruggen D
15. Guo J
16. He X
17. Barker R
18. Sundström E
19. Castelo-Branco G
20. Cramer P
21. Adameyko I
22. Linnarsson S
23. Kharchenko PV
(2018) RNA velocity of single cells
Nature 560:494–498.
https://doi.org/10.1038/s41586-018-0414-6
- PubMed
- Google Scholar
1. Lantz O
2. Bendelac A
(1994) An invariant T cell receptor alpha chain is used by a unique subset of major histocompatibility complex class I-specific CD4+ and CD4-8- T cells in mice and humans
The Journal of Experimental Medicine 180:1097–1106.
https://doi.org/10.1084/jem.180.3.1097
- PubMed
- Google Scholar
1. Li XM
2. Chen X
3. Gu W
4. Guo YJ
5. Cheng Y
6. Peng J
7. Guo XJ
(2017) Impaired TNF/TNFR2 signaling enhances th2 and th17 polarization and aggravates allergic airway inflammation
American Journal of Physiology. Lung Cellular and Molecular Physiology 313:L592–L601.
https://doi.org/10.1152/ajplung.00409.2016
- PubMed
- Google Scholar
1. Liang HE
2. Reinhardt RL
3. Bando JK
4. Sullivan BM
5. Ho IC
6. Locksley RM
(2011) Divergent expression patterns of IL-4 and IL-13 define unique functions in allergic immunity
Nature Immunology 13:58–66.
https://doi.org/10.1038/ni.2182
- PubMed
- Google Scholar
(2013) Th9 cells drive host immunity against gastrointestinal worm infection
Immunity 39:744–757.
https://doi.org/10.1016/j.immuni.2013.07.020
- PubMed
- Google Scholar
1. Lippert E
2. Yowe DL
3. Gonzalo JA
4. Justice JP
5. Webster JM
6. Fedyk ER
7. Hodge M
8. Miller C
9. Gutierrez-Ramos JC
10. Borrego F
11. Keane-Myers A
12. Druey KM
(2003) Role of regulator of G protein signaling 16 in inflammation-induced T lymphocyte migration and activation
Journal of Immunology 171:1542–1555.
https://doi.org/10.4049/jimmunol.171.3.1542
- PubMed
- Google Scholar
1. Macosko EZ
2. Basu A
3. Satija R
4. Nemesh J
5. Shekhar K
6. Goldman M
7. Tirosh I
8. Bialas AR
9. Kamitaki N
10. Martersteck EM
11. Trombetta JJ
12. Weitz DA
13. Sanes JR
14. Shalek AK
15. Regev A
16. McCarroll SA
(2015) Highly parallel genome-wide expression profiling of individual cells using nanoliter droplets
Cell 161:1202–1214.
https://doi.org/10.1016/j.cell.2015.05.002
- PubMed
- Google Scholar
1. Marsland BJ
2. Kurrer M
3. Reissmann R
4. Harris NL
5. Kopf M
(2008) Nippostrongylus brasiliensis infection leads to the development of emphysema associated with the induction of alternatively activated macrophages
European Journal of Immunology 38:479–488.
https://doi.org/10.1002/eji.200737827
- PubMed
- Google Scholar
(2001) Regulation and function of T1/ST2 expression on CD4+ T cells: induction of type 2 cytokine production by T1/ST2 cross-linking
Journal of Immunology 166:3143–3150.
https://doi.org/10.4049/jimmunol.166.5.3143
- PubMed
- Google Scholar
1. Melsted P
2. Booeshaghi AS
3. Liu L
4. Gao F
5. Lu L
6. Min KHJ
7. da Veiga Beltrame E
8. Hjörleifsson KE
9. Gehring J
10. Pachter L
(2021) Modular, efficient and constant-memory single-cell RNA-seq preprocessing
Nature Biotechnology 39:813–818.
https://doi.org/10.1038/s41587-021-00870-2
- PubMed
- Google Scholar
1. Minutti CM
2. Modak RV
3. Macdonald F
4. Li F
5. Smyth DJ
6. Dorward DA
7. Blair N
8. Husovsky C
9. Muir A
10. Giampazolias E
11. Dobie R
12. Maizels RM
13. Kendall TJ
14. Griggs DW
15. Kopf M
16. Henderson NC
17. Zaiss DM
(2019) A macrophage-pericyte axis directs tissue restoration via amphiregulin-induced transforming growth factor beta activation
Immunity 50:645–654.
https://doi.org/10.1016/j.immuni.2019.01.008
- PubMed
- Google Scholar
1. Mohrs M
2. Shinkai K
3. Mohrs K
4. Locksley RM
(2001) Analysis of type 2 immunity in vivo with a bicistronic IL-4 reporter
Immunity 15:303–311.
https://doi.org/10.1016/s1074-7613(01)00186-8
- PubMed
- Google Scholar
1. Mohrs K
2. Wakil AE
3. Killeen N
4. Locksley RM
5. Mohrs M
(2005) A two-step process for cytokine production revealed by IL-4 dual-reporter mice
Immunity 23:419–429.
https://doi.org/10.1016/j.immuni.2005.09.006
- PubMed
- Google Scholar
1. Morimoto Y
2. Hirahara K
3. Kiuchi M
4. Wada T
5. Ichikawa T
6. Kanno T
7. Okano M
8. Kokubo K
9. Onodera A
10. Sakurai D
11. Okamoto Y
12. Nakayama T
(2018) Amphiregulin-producing pathogenic memory T helper 2 cells instruct eosinophils to secrete osteopontin and facilitate airway fibrosis
Immunity 49:134–150.
https://doi.org/10.1016/j.immuni.2018.04.023
- PubMed
- Google Scholar
1. Na B-R
2. Kim H-R
3. Piragyte I
4. Oh H-M
5. Kwon M-S
6. Akber U
7. Lee H-S
8. Park D-S
9. Song WK
10. Park Z-Y
11. Im S-H
12. Rho M-C
13. Hyun Y-M
14. Kim M
15. Jun C-D
(2015) TAGLN2 regulates T cell activation by stabilizing the actin cytoskeleton at the immunological synapse
The Journal of Cell Biology 209:143–162.
https://doi.org/10.1083/jcb.201407130
- PubMed
- Google Scholar
1. Naus S
2. Blanchet MR
3. Gossens K
4. Zaph C
5. Bartsch JW
6. McNagny KM
7. Ziltener HJ
(2010) The metalloprotease-disintegrin ADAM8 is essential for the development of experimental asthma
American Journal of Respiratory and Critical Care Medicine 181:1318–1328.
https://doi.org/10.1164/rccm.200909-1396OC
- PubMed
- Google Scholar
(1993) Expression of two T cell receptor alpha chains: dual receptor T cells
Science 262:422–424.
https://doi.org/10.1126/science.8211163
- PubMed
- Google Scholar
(1995) Normal T lymphocytes can express two different T cell receptor beta chains: implications for the mechanism of allelic exclusion
The Journal of Experimental Medicine 181:1587–1591.
https://doi.org/10.1084/jem.181.4.1587
- PubMed
- Google Scholar
1. Panzer M
2. Sitte S
3. Wirth S
4. Drexler I
5. Sparwasser T
6. Voehringer D
(2012) Rapid in vivo conversion of effector T cells into th2 cells during helminth infection
Journal of Immunology 188:615–623.
https://doi.org/10.4049/jimmunol.1101164
- PubMed
- Google Scholar
1. Peine M
2. Rausch S
3. Helmstetter C
4. Fröhlich A
5. Hegazy AN
6. Kühl AA
7. Grevelding CG
8. Höfer T
9. Hartmann S
10. Löhning M
(2013) Stable T-bet(+)GATA-3(+) th1/th2 hybrid cells arise in vivo, can develop directly from naive precursors, and limit immunopathologic inflammation
PLOS Biology 11:e1001633.
https://doi.org/10.1371/journal.pbio.1001633
- PubMed
- Google Scholar
1. Proserpio V
2. Piccolo A
3. Haim-Vilmovsky L
4. Kar G
5. Lönnberg T
6. Svensson V
7. Pramanik J
8. Natarajan KN
9. Zhai W
10. Zhang X
11. Donati G
12. Kayikci M
13. Kotar J
14. McKenzie ANJ
15. Montandon R
16. Billker O
17. Woodhouse S
18. Cicuta P
19. Nicodemi M
20. Teichmann SA
(2016) Single-cell analysis of CD4+ T-cell differentiation reveals three major cell states and progressive acceleration of proliferation
Genome Biology 17:103.
https://doi.org/10.1186/s13059-016-0957-5
- PubMed
- Google Scholar
Software
1. R Development Core Team
(2021) R: A language and environment for statistical computing
R Foundation for Statistical Computing, Vienna, Austria.

https://www.r-project.org/index.html
1. Radtke D
2. Bannard O
(2018) Expression of the plasma cell transcriptional regulator blimp-1 by dark zone germinal center B cells during periods of proliferation
Frontiers in Immunology 9:3106.
https://doi.org/10.3389/fimmu.2018.03106
- PubMed
- Google Scholar
(2020) Distinct functions of tissue-resident and circulating memory th2 cells in allergic airway disease
The Journal of Experimental Medicine 217:e20190865.
https://doi.org/10.1084/jem.20190865
- PubMed
- Google Scholar
(2009) Cytokine-secreting follicular T cells shape the antibody repertoire
Nature Immunology 10:385–393.
https://doi.org/10.1038/ni.1715
- PubMed
- Google Scholar
1. Rock EP
2. Sibbald PR
3. Davis MM
4. Chien YH
(1994) CDR3 length in antigen-specific immune receptors
The Journal of Experimental Medicine 179:323–328.
https://doi.org/10.1084/jem.179.1.323
- PubMed
- Google Scholar
1. Rothenberg ME
2. Hogan SP
(2006) The eosinophil
Annual Review of Immunology 24:147–174.
https://doi.org/10.1146/annurev.immunol.24.021605.090720
- PubMed
- Google Scholar
(2015) Helminth infection and commensal microbiota drive early IL-10 production in the skin by CD4+ T cells that are functionally suppressive
PLOS Pathogens 11:e1004841.
https://doi.org/10.1371/journal.ppat.1004841
- PubMed
- Google Scholar
1. Savage AK
2. Constantinides MG
3. Han J
4. Picard D
5. Martin E
6. Li B
7. Lantz O
8. Bendelac A
(2008) The transcription factor PLZF directs the effector program of the NKT cell lineage
Immunity 29:391–403.
https://doi.org/10.1016/j.immuni.2008.07.011
- PubMed
- Google Scholar
1. Schaerli P
2. Willimann K
3. Lang AB
4. Lipp M
5. Loetscher P
6. Moser B
(2000) CXC chemokine receptor 5 expression defines follicular homing T cells with B cell helper function
The Journal of Experimental Medicine 192:1553–1562.
https://doi.org/10.1084/jem.192.11.1553
- PubMed
- Google Scholar
(2019) Transcriptional programming and T cell receptor repertoires distinguish human lung and lymph node memory T cells
Communications Biology 2:411.
https://doi.org/10.1038/s42003-019-0657-2
- PubMed
- Google Scholar
(2011) Protective immunity against the gastrointestinal nematode nippostrongylus brasiliensis requires a broad T-cell receptor repertoire
Immunology 134:214–223.
https://doi.org/10.1111/j.1365-2567.2011.03480.x
- PubMed
- Google Scholar
(2019) The metabolic requirements of th2 cell differentiation
Frontiers in Immunology 10:2318.
https://doi.org/10.3389/fimmu.2019.02318
- PubMed
- Google Scholar
1. Stenstad H
2. Ericsson A
3. Johansson-Lindbom B
4. Svensson M
5. Marsal J
6. Mack M
7. Picarella D
8. Soler D
9. Marquez G
10. Briskin M
11. Agace WW
(2006) Gut-associated lymphoid tissue-primed CD4+ T cells display CCR9-dependent and -independent homing to the small intestine
Blood 107:3447–3454.
https://doi.org/10.1182/blood-2005-07-2860
- PubMed
- Google Scholar
1. Stuart T
2. Butler A
3. Hoffman P
4. Hafemeister C
5. Papalexi E
6. Mauck WM
7. Hao Y
8. Stoeckius M
9. Smibert P
10. Satija R
(2019) Comprehensive integration of single-cell data
Cell 177:1888–1902.
https://doi.org/10.1016/j.cell.2019.05.031
- PubMed
- Google Scholar
1. Subramanian A
2. Tamayo P
3. Mootha VK
4. Mukherjee S
5. Ebert BL
6. Gillette MA
7. Paulovich A
8. Pomeroy SL
9. Golub TR
10. Lander ES
11. Mesirov JP
(2005) Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles
PNAS 102:15545–15550.
https://doi.org/10.1073/pnas.0506580102
- PubMed
- Google Scholar
1. Szabo PA
2. Levitin HM
3. Miron M
4. Snyder ME
5. Senda T
6. Yuan J
7. Cheng YL
8. Bush EC
9. Dogra P
10. Thapa P
11. Farber DL
12. Sims PA
(2019a) Single-cell transcriptomics of human T cells reveals tissue and activation signatures in health and disease
Nature Communications 10:4706.
https://doi.org/10.1038/s41467-019-12464-3
- PubMed
- Google Scholar
(2019b) Location, location, location: tissue resident memory T cells in mice and humans
Science Immunology 4:34.
https://doi.org/10.1126/sciimmunol.aas9673
- PubMed
- Google Scholar
(2009) Interleukin 5 in the link between the innate and acquired immune response
Advances in Immunology 101:191–236.
https://doi.org/10.1016/S0065-2776(08)01006-7
- PubMed
- Google Scholar
1. Tibbitt CA
2. Stark JM
3. Martens L
4. Ma J
5. Mold JE
6. Deswarte K
7. Oliynyk G
8. Feng X
9. Lambrecht BN
10. De Bleser P
11. Nylén S
12. Hammad H
13. Arsenian Henriksson M
14. Saeys Y
15. Coquet JM
(2019) Single-cell RNA sequencing of the T helper cell response to house dust mites defines a distinct gene expression signature in airway th2 cells
Immunity 51:169–184.
https://doi.org/10.1016/j.immuni.2019.05.014
- PubMed
- Google Scholar
1. Tirosh I
2. Izar B
3. Prakadan SM
4. Wadsworth MH
5. Treacy D
6. Trombetta JJ
7. Rotem A
8. Rodman C
9. Lian C
10. Murphy G
11. Fallahi-Sichani M
12. Dutton-Regester K
13. Lin J-R
14. Cohen O
15. Shah P
16. Lu D
17. Genshaft AS
18. Hughes TK
19. Ziegler CGK
20. Kazer SW
21. Gaillard A
22. Kolb KE
23. Villani A-C
24. Johannessen CM
25. Andreev AY
26. Van Allen EM
27. Bertagnolli M
28. Sorger PK
29. Sullivan RJ
30. Flaherty KT
31. Frederick DT
32. Jané-Valbuena J
33. Yoon CH
34. Rozenblatt-Rosen O
35. Shalek AK
36. Regev A
37. Garraway LA
(2016) Dissecting the multicellular ecosystem of metastatic melanoma by single-cell RNA-seq
Science 352:189–196.
https://doi.org/10.1126/science.aad0501
- PubMed
- Google Scholar
(2020) High-dimensional T helper cell profiling reveals a broad diversity of stably committed effector states and uncovers interlineage relationships
Immunity 53:597–613.
https://doi.org/10.1016/j.immuni.2020.07.001
- PubMed
- Google Scholar
1. Trapnell C
2. Cacchiarelli D
3. Grimsby J
4. Pokharel P
5. Li S
6. Morse M
7. Lennon NJ
8. Livak KJ
9. Mikkelsen TS
10. Rinn JL
(2014) The dynamics and regulators of cell fate decisions are revealed by pseudotemporal ordering of single cells
Nature Biotechnology 32:381–386.
https://doi.org/10.1038/nbt.2859
- PubMed
- Google Scholar
1. Trinh HKT
2. Nguyen TVT
3. Choi Y
4. Park HS
5. Shin YS
(2019) The synergistic effects of clopidogrel with montelukast may be beneficial for asthma treatment
Journal of Cellular and Molecular Medicine 23:3441–3450.
https://doi.org/10.1111/jcmm.14239
- PubMed
- Google Scholar
(2014) B-cell-intrinsic STAT6 signaling controls germinal center formation
European Journal of Immunology 44:2130–2138.
https://doi.org/10.1002/eji.201344203
- PubMed
- Google Scholar
1. Upadhyaya B
2. Yin Y
3. Hill BJ
4. Douek DC
5. Prussin C
(2011) Hierarchical IL-5 expression defines a subpopulation of highly differentiated human th2 cells
Journal of Immunology 187:3111–3120.
https://doi.org/10.4049/jimmunol.1101283
- PubMed
- Google Scholar
1. Urban JF
2. Madden KB
3. Svetić A
4. Cheever A
5. Trotta PP
6. Gause WC
7. Katona IM
8. Finkelman FD
(1992) The importance of th2 cytokines in protective immunity to nematodes
Immunological Reviews 127:205–220.
https://doi.org/10.1111/j.1600-065x.1992.tb01415.x
- PubMed
- Google Scholar
(2015) FlowSOM: using self-organizing maps for visualization and interpretation of cytometry data
Cytometry. Part A 87:636–645.
https://doi.org/10.1002/cyto.a.22625
- PubMed
- Google Scholar
1. Veldhoen M
2. Uyttenhove C
3. van Snick J
4. Helmby H
5. Westendorf A
6. Buer J
7. Martin B
8. Wilhelm C
9. Stockinger B
(2008) Transforming growth factor-beta “reprograms” the differentiation of T helper 2 cells and promotes an interleukin 9-producing subset
Nature Immunology 9:1341–1346.
https://doi.org/10.1038/ni.1659
- PubMed
- Google Scholar
1. Walker JA
2. McKenzie ANJ
(2018) TH2 cell development and function
Nature Reviews Immunology 18:121–133.
https://doi.org/10.1038/nri.2017.118
- PubMed
- Google Scholar
1. Webb LM
2. Lundie RJ
3. Borger JG
4. Brown SL
5. Connor LM
6. Cartwright AN
7. Dougall AM
8. Wilbers RH
9. Cook PC
10. Jackson-Jones LH
11. Phythian-Adams AT
12. Johansson C
13. Davis DM
14. Dewals BG
15. Ronchese F
16. MacDonald AS
(2017) Type I interferon is required for T helper (th) 2 induction by dendritic cells
The EMBO Journal 36:2404–2418.
https://doi.org/10.15252/embj.201695345
- PubMed
- Google Scholar
(2001) The role of IFN-gamma in the production of th2 subpopulations: implications for variable th2-mediated pathologies in autoimmunity
Journal of Immunology 167:3074–3081.
https://doi.org/10.4049/jimmunol.167.6.3074
- PubMed
- Google Scholar
1. Zheng GXY
2. Terry JM
3. Belgrader P
4. Ryvkin P
5. Bent ZW
6. Wilson R
7. Ziraldo SB
8. Wheeler TD
9. McDermott GP
10. Zhu J
11. Gregory MT
12. Shuga J
13. Montesclaros L
14. Underwood JG
15. Masquelier DA
16. Nishimura SY
17. Schnall-Levin M
18. Wyatt PW
19. Hindson CM
20. Bharadwaj R
21. Wong A
22. Ness KD
23. Beppu LW
24. Deeg HJ
25. McFarland C
26. Loeb KR
27. Valente WJ
28. Ericson NG
29. Stevens EA
30. Radich JP
31. Mikkelsen TS
32. Hindson BJ
33. Bielas JH
(2017) Massively parallel digital transcriptional profiling of single cells
Nature Communications 8:14049.
https://doi.org/10.1038/ncomms14049
- PubMed
- Google Scholar

Article and author information

Author details

Daniel Radtke

Department of Infection Biology, Universitätsklinikum Erlangen and Friedrich-Alexander, University Erlangen-Nürnberg, Erlangen, Germany

Contribution
Conceptualization, Data curation, Software, Investigation, Visualization, Methodology, Writing - original draft, Writing - review and editing

For correspondence
daniel.radtke@uk-erlangen.de

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0003-3241-4542
Natalie Thuma

Department of Infection Biology, Universitätsklinikum Erlangen and Friedrich-Alexander, University Erlangen-Nürnberg, Erlangen, Germany

Contribution
Conceptualization, Investigation, Visualization, Methodology, Writing - original draft, Writing - review and editing

Competing interests
No competing interests declared
Christine Schülein

Institute of Clinical Microbiology, Immunology and Hygiene, University Hospital Erlangen and Friedrich-Alexander-University Erlangen-Nürnberg, Erlangen, Germany

Contribution
Investigation

Competing interests
No competing interests declared
Philipp Kirchner

Institute of Human Genetics, Universitätsklinikum Erlangen and Friedrich-Alexander University Erlangen-Nürnberg, Erlangen, Germany

Contribution
Methodology

Competing interests
No competing interests declared
Arif B Ekici

Institute of Human Genetics, Universitätsklinikum Erlangen and Friedrich-Alexander University Erlangen-Nürnberg, Erlangen, Germany

Contribution
Methodology

Competing interests
No competing interests declared
Kilian Schober

Institute of Clinical Microbiology, Immunology and Hygiene, University Hospital Erlangen and Friedrich-Alexander-University Erlangen-Nürnberg, Erlangen, Germany

Contribution
Conceptualization

Competing interests
No competing interests declared
David Voehringer

Department of Infection Biology, Universitätsklinikum Erlangen and Friedrich-Alexander, University Erlangen-Nürnberg, Erlangen, Germany

Contribution
Conceptualization, Supervision, Funding acquisition, Writing - original draft, Writing - review and editing

For correspondence
David.Voehringer@uk-erlangen.de

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0001-6650-0639

Funding

Deutsche Forschungsgemeinschaft (RTG1660)

David Voehringer

Deutsche Forschungsgemeinschaft (FOR2599_TP4)

David Voehringer

Deutsche Forschungsgemeinschaft (TRR130_TP20)

David Voehringer

The funders had no role in study design, data collection, and interpretation, or the decision to submit the work for publication.

Acknowledgements

The authors thank the Flow cytometry core facility of the FAU for cell sorting, Daniela Döhler and Kirstin Castiglione for technical assistance and members of the Voehringer lab for helpful discussions.

Ethics

All experiments were performed in accordance with German animal protection law and European Union guidelines 86/809 and were approved by the Federal Government of Lower Franconia.

Version history

Received: September 24, 2021
Preprint posted: October 3, 2021 (view preprint)
Accepted: August 10, 2022
Accepted Manuscript published: August 11, 2022 (version 1)
Version of Record published: August 19, 2022 (version 2)

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.