Germinal centers (GCs) are transient organized microstructures within secondary lymphoid tissues, which support the development of high-affinity antibodies by B cells. Even though GCs are mainly B cell driven, their formation and maintenance depends on follicular T helper cells (Tfh), a specialized subset of effector CD4 T cells. The differentiation into Tfh initiates in T cell zones of lymphoid organs before GC formation starts. After their interaction with B cells at the T-B border, Tfh locate within GC (GC-Tfh) and regulate the survival of proliferating GC-B cells, which compete for antigen (Ag) and for signals from GC-Tfh to undergo further somatic hypermutation and to mature into high-affinity-antibody-producing plasma cells and memory B cells (Crotty, 2019; Qi, 2016). Thus, GC-Tfh are central players in the regulation of humoral immune responses. Both (i) their differentiation into GC-Tfh and (ii) their survival and clonal expansion within GCs are constantly driven by competition for Ag contacts (Baumjohann et al., 2013; Fazilleau et al., 2009; Hwang et al., 2015; Knowlden and Sant, 2016; Tubo et al., 2013; Merkenschlager et al., 2021).

Here, we asked whether this clonal competition of Tfh for the limited space in GCs involves GC reactions of one or more activated lymphoid organs. This is of special interest for autoimmune conditions when multiple inflamed tissue sites and autoreactive GC reactions exist in distinct lymphoid organs. Tfh can shuttle between GCs of one lymph node or spleen indicating a local competition within one lymphoid organ (Merkenschlager et al., 2021; Shulman et al., 2013). The finding that CXCR5+ Tfh circulate in blood (Brenna et al., 2020; He et al., 2013) opens the possibility for a systemic exchange, in which each Tfh clone could participate in the clonal competition within distinct lymphoid organs in one individual. However, the diversity of individual T cell receptor (TCR) repertoires is extremely high, and it has been shown that individual T cell responses towards Ag are unique and polyclonal (Textor et al., 2018; Thomas et al., 2014; Zarnitsyna et al., 2013). This might be especially valid for auto-Ag because TCR bind with low affinities (Dolton et al., 2018; Rius et al., 2018). In line, the number of autoreactive T cell clones within one individual is particularly low due to the thymic negative selection process, which could support the involvement of multiple cross-reactive T cell clones specific for the same auto-Ag within one individual as it has been suggested previously (Ritvo et al., 2018).

To find out how polyclonal autoreactive Tfh responses are, we studied the distribution of GC-Tfh clones between two activated lymph nodes and challenged the endogenous TCR repertoire by injecting one auto-Ag, directed against a structural protein of murine skin, into the two hind footpads of the same mouse and analyzed the TCRβ sequences of GC-Tfh in both draining popliteal lymph nodes (pln). In contrast to our expectations, we found that the GC-Tfh clonotypes were highly shared between the two pln. This high overlap between both pln was restricted to the dominant GC-Tfh clonotypes and disappeared after establishment of skin pathology even though GC reaction continued. These temporal variations in the distribution of Tfh clones should be considered when utilizing Tfh repertoires to monitor and diagnose autoimmune diseases.

CD4 T cells are pivotal in adaptive immune responses during infection, autoimmune diseases, cancer, and vaccinations. Each CD4 T cell expresses one specific αβ TCR interacting with specific peptide-Ag presented in MHCII complexes (Itano and Jenkins, 2003). Upon interaction of the TCR with peptide:MHCII complexes, naive CD4 T cells undergo proliferation and functional differentiation into distinct T helper subsets such as Th1, Th2, Th17 effector cells, Tfh, and immunosuppressive Treg cells (Zhu et al., 2010). The main factor that regulates CD4 T cell differentiation is the TCR signaling strength, which is controlled by TCR/peptide:MHCII interactions, co-stimulation, and/or optimal dwell times. The differentiation into Tfh depends on strong TCR signals. It is suggested that they are strictly Ag-specific (Hwang et al., 2015; Knowlden and Sant, 2016; Tubo et al., 2013; Merkenschlager et al., 2021). Consistently, it has been shown that the maintenance and expansion of Tfh requires sustained Ag stimulation (Baumjohann et al., 2013; Merkenschlager et al., 2021).

In this study, we asked how this particular CD4 T cell subset is distributed within one individual at the clonal level. We used an autoantibody-mediated disease model of the skin, epidermolysis bullosa acquisita, and induced GC-Tfh in two separate draining pln by injecting auto-Ag in adjuvant (Ag1 or Ag2 group) or adjuvant only (PBS group).

Our study shows two major findings. First, we found a high overlap of dominant GC-Tfh clonotypes in both pln of one mouse (Figure 2a, b). This data demonstrates that even though the response to the auto-Ag was polyclonal in each mouse, a group of almost identical GC-Tfh clonotypes became dominant in each draining pln. Thus, exposures to the auto-Ag (Ag1 or Ag2) decreased the diversity of the endogenous Tfh repertoire compared to the PBS group. To establish this finding, we used multiplex iRepertiore PCRs and subsequent MiTCR analysis for preparing libraries and annotating TCRβ sequences (Bolotin et al., 2013), which ensures to yield the highest possible diversity of TCRβ clonotypes (Afzal et al., 2019). To confirm this decrease in diversity of the endogenous Tfh repertoire, we reanalyzed our data for Ag1 with a more precise analysis tool MiXCR that is advantageous for error corrections and adjusts for a more accurate clonal composition (Bolotin et al., 2013; Bolotin et al., 2015; Team I, 2019). The high overlap of the dominant Tfh clonotypes in the Ag1 group was clearly confirmed (Figure 2—figure supplement 1). To further compare both analysis tools, we calculated the MHI from samples of the Ag1 and PBS groups analyzed either with MiTCR or MiXCR. Data show a high similarity in both the Ag1 group (0.70 ± 0.04) and the PBS group (0.69 ± 0.039). Obviously, our finding that Tfh clonotypes highly overlap in contralateral pln is very robust independent of the analysis tools. On the other side, the suitability of the multiplex PCRs can be seen by the correct detections of the deletions in the TRBV genes of SJL mice (Figure 3c; Behlke et al., 1986). In summary, this data that the endogenous Tfh repertoire decreased after Ag exposure indicates that the number of T cell clones specific for the injected auto-Ag is restricted per mouse and supports previous reports demonstrating that the precursor frequency specific for auto-Ag within the T cell population of an individual is rather low (Jenkins and Moon, 2012). However, the precursor frequency of naive T cells specific for different peptide:MHC complexes varies in size. Therefore, it will be interesting whether this high clonal overlap would be found after injection of xenogeneic or pathogenic Ag. Moreover, because both Ag used in this study are polypeptides of a relatively large size (200–400 aa) and clearly differ from the size of peptides (13–25 aa) that are usually presented in MHCII (Natarajan et al., 2018), it is likely that they contain more than one antigenic epitope. It is reasonable to assume that injecting smaller peptides will further narrow the diversity of the individual GC-Tfh repertoire. In an initial experiment as an example for a recombinant xenogeneic protein antigen, we injected the Schistosoma mansoni-derived glutathione-s-transferase (GST)-tag of Ag1 (Rao et al., 2003), emulsified in TM into both footpads, and analyzed the distribution of Tfh clonotypes in both draining pln 4 weeks p.i. (Figure 2—figure supplement 2, Supplementary file 6). We found that the distribution of the GST-specific Tfh clonotypes between both pln was less controlled compared to the auto-Ag1. Dot plot diagrams showed a higher variability between the mice. In line, the MHI was not significantly different from the PBS group (Figure 2—figure supplement 2) and the number of overlapping Tfh clonotypes is lower in the GST group compared to the Ag1 group (Figure 3—figure supplement 1). Analysis of the V/J usage revealed significant differences in expression of TRBV3. The J2-4-segment was significantly less abundant after GST/TM immunization compared to the PBS/-group and the J2-5 segment was jointly higher abundant after GST/TM and Ag1 immunization compared to the control group (Figure 3—figure supplement 1). These data confirm the strict Ag specificity of Tfh cells in GCs, especially during the first 4 weeks p.i. The higher inconsistency in the number of overlapping Tfh clonotypes in the GST group could be explained by an earlier consumption of the Ag or the presence of more high-affinity T cell clones. These findings might be strongly influenced by the time point and the Ag concentration, which should be further investigated in future studies.

Another arising question is when this overlap of T cell clonotypes between both pln of one mouse starts. For example, T cells could travel between pln and the clonal selection could be a continuous organism-wide systemic process. The presence of circulating Tfh cells has been described in humans and mice (Brenna et al., 2020; He et al., 2013); however, it is not clear whether these circulating Tfh cells can enter ongoing GC reactions. Alternatively, it could take place during the priming phase. In this case, either the naive TCR repertoire is sufficiently broad that the Ag-driven selection for the GC reactions in both pln would robustly yield groups of similar clones independently or recently activated CD4 T cell clones would exchange between both pln after priming before GC reaction starts. To address this, we performed an initial experiment and identified TCRβ sequences of entire left and right pln of the same mouse (naive mice and mice that had been exposed to Ag1 for 1 or 3 days; please find detailed information in Supplementary file 5) and compared the relative percentage of overlapping T cell clonotypes between left and right pln by using the MHI. We found that the MHI did not differ between naive mice and Ag-exposed mice at 1 day p.i. (0.2035 + 0.01962 in naive mice, 0.2113 + 0.009 at 1 day p.i.), but a significant increased MHI at 3 days p.i. (0.3716 + 0.031; p<0.01, Kruskal–Wallis test, mean + SD for n = 3 mice, Figure 4—figure supplement 1). In parallel, the number of proliferating T cells and the mRNA expression for the proliferation marker Ki67 increased significantly (Figure 4—figure supplement 1). One may speculate that a high number of progenies emerges from the high-affinity T cell clones 3 days p.i. (Cho et al., 2017), which then migrate into other activated lymphoid tissues before differentiating into GC-Tfh. In this case, the initial progeny of activated T cell clones could originate from identical precursors, which would enter the circulation and migrate into the other pln, in which they outcompete low-affinity competitors and develop into GC-Tfh. However, on the other side it is assumed that Tfh stay resident in lymph nodes to provide B cell help (Crotty, 2019). One possibility to observe potential exchanges of Tfh cells between lymph nodes would be treatment with the T cell trafficking blocker FTY720. However, previous studies revealed no effect in this autoimmune mouse model (Niebuhr et al., 2017; Thieme et al., 2019). The second major finding of our study is that the PBS group, which received only adjuvant, develops GCs to a similar extent as the Ag group. Here, in contrast to the Ag group, the GC-Tfh repertoire remains diverse. No overlap of dominant GC-Tfh clonotypes was found between both pln. TM is a water-in-oil adjuvant that contains squalene among other components but no peptides that could be presented in MHCII (TiterMax Inc). It is proposed that adjuvants stress or kill local cells, which leads to the release of endogenous peptide-epitopes that are subsequently presented to T cells (Riteau et al., 2016) and induce the formation of GC reactions as observed in this study (Figure 1c). In this scenario, numerous endogenous peptide-epitopes at low concentrations might be released by the local inflammation, which leads to the diverse GC-Tfh repertoire in both pln (Figure 2a, b). However, to break the tolerance in our model, both auto-Ag had to be injected at very high concentrations (60 µg/Ag1 or 120 µg/Ag2, respectively) (Sitaru et al., 2006). This data leads to the hypothesis that the concentration of the Ag is critically involved in the emergence of overlapping GC-Tfh repertoires in separate pln (summarized in Figure 6). It will be interesting to find out how low Ag concentrations at distinct tissue sites affect the GC-Tfh repertoires. Moreover, it will be important to study the adjuvant-stimulated peptide-epitopes in future studies, especially because squalene compounds are frequently used in human vaccines and might evoke unwanted autoimmune responses (Pellegrini et al., 2009).

Figure 6

Download asset Open asset

Another conclusion of our study is that the total number of GC-Tfh within one individual is regulated by the number of lymphoid tissues that drain Ag-exposed tissue sites. This data suggests that Ag injection at two or more tissue sites could improve the development of high-affinity antibodies. We believe that this information will be highly useful for the development of new and more efficient vaccination strategies.

Finally, we wish to highlight that this in vivo approach could help to identify new peptide-specific TCR sequences and confirm already established datasets (Teraguchi et al., 2020). Even though TCR datasets are growing, the number of TCRs with known Ag specificity and function is extremely low (Shugay et al., 2018; Zvyagin et al., 2020). Especially the detection of self-reactive CD4 TCR is problematic due to the low binding affinities between TCR and autoantigenic-peptide:MHCII complexes (Dolton et al., 2018; Rius et al., 2018). Thus, even though our approach is restricted to mouse models, the use of HLA-transgenic mice might enable to link TCR sequences also to human antigenic epitopes.

T cell receptor-RNA sequencing data generated in this study are deposited in the sequence read archive hosted at https://www.ncbi.nlm.nih.gov/sra under the primary accession code PRJNA731654. Some of the data (skin effector T cells) have been deposited under the accession code PRJNA586880.

1. 1. Kalies K
   2. Niebuhr M

   (2021) NCBI BioProject

   ID PRJNA731654. TRBV sequences of murine follicular T helper cells in experimental epidermolysis bullosa acquisita.

   https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA731654
2. 1. Kalies K
   2. Niebuhr M

   (2021) NCBI BioProject

   ID PRJNA586880. TCRb-seq data from dermal T cells in SJL mice upon immunization.

   https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA586880

1. 1. Afzal S
   2. Gil-Farina I
   3. Gabriel R
   4. Ahmad S
   5. von Kalle C
   6. Schmidt M
   7. Fronza R

   (2019) Systematic comparative study of computational methods for T-cell receptor sequencing data analysis

   Briefings in Bioinformatics 20:222–234.

   https://doi.org/10.1093/bib/bbx111

   • PubMed
   • Google Scholar
2. 1. Baumjohann D
   2. Preite S
   3. Reboldi A
   4. Ronchi F
   5. Ansel KM
   6. Lanzavecchia A
   7. Sallusto F

   (2013) Persistent antigen and germinal center B cells sustain T follicular helper cell responses and phenotype

   Immunity 38:596–605.

   https://doi.org/10.1016/j.immuni.2012.11.020

   • PubMed
   • Google Scholar
3. 1. Behlke MA
   2. Chou HS
   3. Huppi K
   4. Loh DY

   (1986) Murine T-cell receptor mutants with deletions of beta-chain variable region genes

   PNAS 83:767–771.

   https://doi.org/10.1073/pnas.83.3.767

   • PubMed
   • Google Scholar
4. 1. Bolotin DA
   2. Shugay M
   3. Mamedov IZ
   4. Putintseva EV
   5. Turchaninova MA
   6. Zvyagin IV
   7. Britanova OV
   8. Chudakov DM

   (2013) MiTCR: software for T-cell receptor sequencing data analysis

   Nature Methods 10:813–814.

   https://doi.org/10.1038/nmeth.2555

   • PubMed
   • Google Scholar
5. 1. Bolotin DA
   2. Poslavsky S
   3. Mitrophanov I
   4. Shugay M
   5. Mamedov IZ
   6. Putintseva EV
   7. Chudakov DM

   (2015) MiXCR: software for comprehensive adaptive immunity profiling

   Nature Methods 12:380–381.

   https://doi.org/10.1038/nmeth.3364

   • PubMed
   • Google Scholar
6. 1. Brenna E
   2. Davydov AN
   3. Ladell K
   4. McLaren JE
   5. Bonaiuti P
   6. Metsger M
   7. Ramsden JD
   8. Gilbert SC
   9. Lambe T
   10. Price DA
   11. Campion SL
   12. Chudakov DM
   13. Borrow P
   14. McMichael AJ

   (2020) CD4⁺ T follicular helper cells in human tonsils and blood are clonally convergent but divergent from Non-Tfh CD4⁺Cells

   Cell Reports 30:137–152.

   https://doi.org/10.1016/j.celrep.2019.12.016

   • PubMed
   • Google Scholar
7. 1. Cho YL
   2. Flossdorf M
   3. Kretschmer L
   4. Höfer T
   5. Busch DH
   6. Buchholz VR

   (2017) TCR signal quality modulates fate decisions of single CD4⁺ T Cells in a Probabilistic Manner

   Cell Reports 20:806–818.

   https://doi.org/10.1016/j.celrep.2017.07.005

   • PubMed
   • Google Scholar
8. 1. Crotty S

   (2019) T follicular helper cell biology: a decade of discovery and diseases

   Immunity 50:1132–1148.

   https://doi.org/10.1016/j.immuni.2019.04.011

   • PubMed
   • Google Scholar
9. 1. Dolton G
   2. Zervoudi E
   3. Rius C
   4. Wall A
   5. Thomas HL
   6. Fuller A
   7. Yeo L
   8. Legut M
   9. Wheeler S
   10. Attaf M
   11. Chudakov DM
   12. Choy E
   13. Peakman M
   14. Sewell AK

   (2018) Optimized Peptide-MHC multimer protocols for detection and isolation of autoimmune T-Cells

   Frontiers in Immunology 9:1378–1391.

   https://doi.org/10.3389/fimmu.2018.01378

   • PubMed
   • Google Scholar
10. 1. Ellebrecht CT
    2. Srinivas G
    3. Bieber K
    4. Banczyk D
    5. Kalies K
    6. Künzel S
    7. Hammers CM
    8. Baines JF
    9. Zillikens D
    10. Ludwig RJ
    11. Westermann J

    (2016) Skin microbiota-associated inflammation precedes autoantibody induced tissue damage in experimental epidermolysis bullosa acquisita

    Journal of Autoimmunity 68:14–22.

    https://doi.org/10.1016/j.jaut.2015.08.007

    • PubMed
    • Google Scholar
11. 1. Fähnrich A
    2. Krebbel M
    3. Decker N
    4. Leucker M
    5. Lange FD
    6. Kalies K
    7. Möller S

    (2017) ClonoCalc and ClonoPlot: immune repertoire analysis from raw files to publication figures with graphical user interface

    BMC Bioinformatics 18:164–170.

    https://doi.org/10.1186/s12859-017-1575-2

    • PubMed
    • Google Scholar
12. 1. Fähnrich A
    2. Klein S
    3. Sergé A
    4. Nyhoegen C
    5. Kombrink S
    6. Möller S
    7. Keller K
    8. Westermann J
    9. Kalies K

    (2018) CD154 costimulation shifts the local T-Cell receptor repertoire not only during thymic selection but also during peripheral T-Dependent humoral immune responses

    Frontiers in Immunology 9:1019–1030.

    https://doi.org/10.3389/fimmu.2018.01019

    • PubMed
    • Google Scholar
13. 1. Fazilleau N
    2. McHeyzer-Williams LJ
    3. Rosen H
    4. McHeyzer-Williams MG

    (2009) The function of follicular helper T cells is regulated by the strength of T cell antigen receptor binding

    Nature Immunology 10:375–384.

    https://doi.org/10.1038/ni.1704

    • PubMed
    • Google Scholar
14. 1. Gaide O
    2. Emerson RO
    3. Jiang X
    4. Gulati N
    5. Nizza S
    6. Desmarais C
    7. Robins H
    8. Krueger JG
    9. Clark RA
    10. Kupper TS

    (2015) Common clonal origin of central and resident memory T cells following skin immunization

    Nature Medicine 21:647–653.

    https://doi.org/10.1038/nm.3860

    • PubMed
    • Google Scholar
15. 1. Hammers CM
    2. Bieber K
    3. Kalies K
    4. Banczyk D
    5. Ellebrecht CT
    6. Ibrahim SM
    7. Zillikens D
    8. Ludwig RJ
    9. Westermann J

    (2011) Complement-fixing anti-type VII collagen antibodies are induced in Th1-polarized lymph nodes of epidermolysis bullosa acquisita-susceptible mice

    The Journal of Immunology 187:5043–5050.

    https://doi.org/10.4049/jimmunol.1100796

    • PubMed
    • Google Scholar
16. 1. He J
    2. Tsai LM
    3. Leong YA
    4. Hu X
    5. Ma CS
    6. Chevalier N
    7. Sun X
    8. Vandenberg K
    9. Rockman S
    10. Ding Y
    11. Zhu L
    12. Wei W
    13. Wang C
    14. Karnowski A
    15. Belz GT
    16. Ghali JR
    17. Cook MC
    18. Riminton DS
    19. Veillette A
    20. Schwartzberg PL
    21. Mackay F
    22. Brink R
    23. Tangye SG
    24. Vinuesa CG
    25. Mackay CR
    26. Li Z
    27. Yu D

    (2013) Circulating precursor CCR7(lo)PD-1(hi) CXCR5⁺ CD4⁺ T cells indicate tfh cell activity and promote antibody responses upon antigen reexposure

    Immunity 39:770–781.

    https://doi.org/10.1016/j.immuni.2013.09.007

    • PubMed
    • Google Scholar
17. 1. Horn HS

    (1966) Measurement of "Overlap" in Comparative Ecological Studies

    The American Naturalist 100:419–424.

    https://doi.org/10.1086/282436

    • Google Scholar
18. 1. Hwang S
    2. Palin AC
    3. Li L
    4. Song KD
    5. Lee J
    6. Herz J
    7. Tubo N
    8. Chu H
    9. Pepper M
    10. Lesourne R
    11. Zvezdova E
    12. Pinkhasov J
    13. Jenkins MK
    14. McGavern D
    15. Love PE

    (2015) TCR ITAM multiplicity is required for the generation of follicular helper T-cells

    Nature Communications 6:6982–6995.

    https://doi.org/10.1038/ncomms7982

    • PubMed
    • Google Scholar
19. 1. Irla M
    2. Guenot J
    3. Sealy G
    4. Reith W
    5. Imhof BA
    6. Sergé A

    (2013) Three-dimensional visualization of the mouse Thymus organization in health and immunodeficiency

    The Journal of Immunology 190:586–596.

    https://doi.org/10.4049/jimmunol.1200119

    • PubMed
    • Google Scholar
20. 1. Itano AA
    2. Jenkins MK

    (2003) Antigen presentation to naive CD4 T cells in the lymph node

    Nature Immunology 4:733–739.

    https://doi.org/10.1038/ni957

    • PubMed
    • Google Scholar
21. 1. Iwata H
    2. Bieber K
    3. Tiburzy B
    4. Chrobok N
    5. Kalies K
    6. Shimizu A
    7. Leineweber S
    8. Ishiko A
    9. Vorobyev A
    10. Zillikens D
    11. Köhl J
    12. Westermann J
    13. Seeger K
    14. Manz R
    15. Ludwig RJ

    (2013) B cells, dendritic cells, and macrophages are required to induce an autoreactive CD4 helper T cell response in experimental epidermolysis bullosa acquisita

    The Journal of Immunology 191:2978–2988.

    https://doi.org/10.4049/jimmunol.1300310

    • Google Scholar
22. 1. Jenkins MK
    2. Moon JJ

    (2012) The role of naive T cell precursor frequency and recruitment in dictating immune response magnitude

    The Journal of Immunology 188:4135–4140.

    https://doi.org/10.4049/jimmunol.1102661

    • Google Scholar
23. 1. Knowlden ZA
    2. Sant AJ

    (2016) CD4 T cell epitope specificity determines follicular versus non-follicular helper differentiation in the polyclonal response to influenza infection or vaccination

    Scientific Reports 6:28287–28300.

    https://doi.org/10.1038/srep28287

    • PubMed
    • Google Scholar
24. 1. Lefranc MP
    2. Giudicelli V
    3. Busin C
    4. Malik A
    5. Mougenot I
    6. Déhais P
    7. Chaume D

    (1995) LIGM-DB/IMGT: an integrated database of ig and TcR, part of the immunogenetics database

    Annals of the New York Academy of Sciences 764:47–49.

    https://doi.org/10.1111/j.1749-6632.1995.tb55805.x

    • PubMed
    • Google Scholar
25. 1. Li KP
    2. Fähnrich A
    3. Roy E
    4. Cuda CM
    5. Grimes HL
    6. Perlman HR
    7. Kalies K
    8. Hildeman DA

    (2017) Temporal expression of bim limits the development of Agonist-Selected thymocytes and skews their TCRbeta repertoire

    Journal of Immunology 198:257–269.

    https://doi.org/10.4049/jimmunol.1601200

    • PubMed
    • Google Scholar
26. 1. Madi A
    2. Shifrut E
    3. Reich-Zeliger S
    4. Gal H
    5. Best K
    6. Ndifon W
    7. Chain B
    8. Cohen IR
    9. Friedman N

    (2014) T-cell receptor repertoires share a restricted set of public and abundant CDR3 sequences that are associated with self-related immunity

    Genome Research 24:1603–1612.

    https://doi.org/10.1101/gr.170753.113

    • PubMed
    • Google Scholar
27. 1. Madi A
    2. Poran A
    3. Shifrut E
    4. Reich-Zeliger S
    5. Greenstein E
    6. Zaretsky I
    7. Arnon T
    8. Laethem FV
    9. Singer A
    10. Lu J
    11. Sun PD
    12. Cohen IR
    13. Friedman N

    (2017) T cell receptor repertoires of mice and humans are clustered in similarity networks around conserved public CDR3 sequences

    eLife 6:e22057.

    https://doi.org/10.7554/eLife.22057

    • PubMed
    • Google Scholar
28. 1. Meli AP
    2. King IL

    (2015) Identification of mouse T follicular helper cells by flow cytometry

    Methods in Molecular Biology 1291:3–11.

    https://doi.org/10.1007/978-1-4939-2498-1_1

    • PubMed
    • Google Scholar
29. 1. Merkenschlager J
    2. Finkin S
    3. Ramos V
    4. Kraft J
    5. Cipolla M
    6. Nowosad CR
    7. Hartweger H
    8. Zhang W
    9. Olinares PDB
    10. Gazumyan A
    11. Oliveira TY
    12. Chait BT
    13. Nussenzweig MC

    (2021) Dynamic regulation of T_FH selection during the germinal centre reaction

    Nature 591:458–463.

    https://doi.org/10.1038/s41586-021-03187-x

    • PubMed
    • Google Scholar
30. 1. Morisita M

    (1962) I σ-Index, a measure of dispersion of individuals

    Researches on Population Ecology 4:1–7.

    https://doi.org/10.1007/BF02533903

    • Google Scholar
31. 1. Natarajan K
    2. Jiang J
    3. May NA
    4. Mage MG
    5. Boyd LF
    6. McShan AC
    7. Sgourakis NG
    8. Bax A
    9. Margulies DH

    (2018) The role of molecular flexibility in antigen presentation and T cell Receptor-Mediated signaling

    Frontiers in Immunology 9:1657.

    https://doi.org/10.3389/fimmu.2018.01657

    • PubMed
    • Google Scholar
32. 1. Nazarov VI
    2. Pogorelyy MV
    3. Komech EA
    4. Zvyagin IV
    5. Bolotin DA
    6. Shugay M
    7. Chudakov DM
    8. Lebedev YB
    9. Mamedov IZ

    (2015) tcR: an R package for T cell receptor repertoire advanced data analysis

    BMC Bioinformatics 16:175–180.

    https://doi.org/10.1186/s12859-015-0613-1

    • PubMed
    • Google Scholar
33. 1. Niebuhr M
    2. Kasperkiewicz M
    3. Maass S
    4. Hauenschild E
    5. Bieber K
    6. Ludwig RJ
    7. Westermann J
    8. Kalies K

    (2017) Evidence for a contributory role of a xenogeneic immune response in experimental epidermolysis bullosa acquisita

    Experimental Dermatology 26:1207–1213.

    https://doi.org/10.1111/exd.13439

    • PubMed
    • Google Scholar
34. 1. Niebuhr M
    2. Bieber K
    3. Banczyk D
    4. Maass S
    5. Klein S
    6. Becker M
    7. Ludwig R
    8. Zillikens D
    9. Westermann J
    10. Kalies K

    (2020) Epidermal damage induces Th1 polarization and defines the site of inflammation in murine epidermolysis bullosa acquisita

    Journal of Investigative Dermatology 140:1713–1722.

    https://doi.org/10.1016/j.jid.2020.01.022

    • Google Scholar
35. 1. Pellegrini M
    2. Nicolay U
    3. Lindert K
    4. Groth N
    5. Della Cioppa G

    (2009) MF59-adjuvanted versus non-adjuvanted influenza vaccines: integrated analysis from a large safety database

    Vaccine 27:6959–6965.

    https://doi.org/10.1016/j.vaccine.2009.08.101

    • PubMed
    • Google Scholar
36. 1. Qi H

    (2016) T follicular helper cells in space-time

    Nature Reviews Immunology 16:612–625.

    https://doi.org/10.1038/nri.2016.94

    • PubMed
    • Google Scholar
37. 1. Rao KV
    2. He YX
    3. Kalyanasundaram R

    (2003) Expression of a 28-kilodalton glutathione S-transferase antigen of Schistosoma mansoni on the surface of filamentous phages and evaluation of its vaccine potential

    Clinical and Vaccine Immunology 10:536–541.

    https://doi.org/10.1128/CDLI.10.4.536-541.2003

    • PubMed
    • Google Scholar
38. 1. Riteau N
    2. Radtke AJ
    3. Shenderov K
    4. Mittereder L
    5. Oland SD
    6. Hieny S
    7. Jankovic D
    8. Sher A

    (2016) Water-in-Oil-Only adjuvants selectively promote T follicular helper cell polarization through a type I IFN and IL-6-Dependent pathway

    The Journal of Immunology 197:3884–3893.

    https://doi.org/10.4049/jimmunol.1600883

    • PubMed
    • Google Scholar
39. 1. Ritvo PG
    2. Saadawi A
    3. Barennes P
    4. Quiniou V
    5. Chaara W
    6. El Soufi K
    7. Bonnet B
    8. Six A
    9. Shugay M
    10. Mariotti-Ferrandiz E
    11. Klatzmann D

    (2018) High-resolution repertoire analysis reveals a major bystander activation of tfh and tfr cells

    PNAS 115:9604–9609.

    https://doi.org/10.1073/pnas.1808594115

    • PubMed
    • Google Scholar
40. 1. Rius C
    2. Attaf M
    3. Tungatt K
    4. Bianchi V
    5. Legut M
    6. Bovay A
    7. Donia M
    8. Thor Straten P
    9. Peakman M
    10. Svane IM
    11. Ott S
    12. Connor T
    13. Szomolay B
    14. Dolton G
    15. Sewell AK

    (2018) Peptide-MHC class I tetramers can fail to detect relevant functional T cell clonotypes and underestimate Antigen-Reactive T cell populations

    The Journal of Immunology 200:2263–2279.

    https://doi.org/10.4049/jimmunol.1700242

    • PubMed
    • Google Scholar
41. 1. Sergé A
    2. Bailly AL
    3. Aurrand-Lions M
    4. Imhof BA
    5. Irla M

    (2015) For3D: full organ reconstruction in 3D, an automatized tool for deciphering the complexity of lymphoid organs

    Journal of Immunological Methods 424:32–42.

    https://doi.org/10.1016/j.jim.2015.04.019

    • PubMed
    • Google Scholar
42. 1. Shugay M
    2. Bagaev DV
    3. Zvyagin IV
    4. Vroomans RM
    5. Crawford JC
    6. Dolton G
    7. Komech EA
    8. Sycheva AL
    9. Koneva AE
    10. Egorov ES
    11. Eliseev AV
    12. Van Dyk E
    13. Dash P
    14. Attaf M
    15. Rius C
    16. Ladell K
    17. McLaren JE
    18. Matthews KK
    19. Clemens EB
    20. Douek DC
    21. Luciani F
    22. van Baarle D
    23. Kedzierska K
    24. Kesmir C
    25. Thomas PG
    26. Price DA
    27. Sewell AK
    28. Chudakov DM

    (2018) VDJdb: a curated database of T-cell receptor sequences with known antigen specificity

    Nucleic Acids Research 46:D419–D427.

    https://doi.org/10.1093/nar/gkx760

    • PubMed
    • Google Scholar
43. 1. Shulman Z
    2. Gitlin AD
    3. Targ S
    4. Jankovic M
    5. Pasqual G
    6. Nussenzweig MC
    7. Victora GD

    (2013) T follicular helper cell dynamics in germinal centers

    Science 341:673–677.

    https://doi.org/10.1126/science.1241680

    • PubMed
    • Google Scholar
44. 1. Sitaru C
    2. Chiriac MT
    3. Mihai S
    4. Büning J
    5. Gebert A
    6. Ishiko A
    7. Zillikens D

    (2006) Induction of complement-fixing autoantibodies against type VII collagen results in subepidermal blistering in mice

    The Journal of Immunology 177:3461–3468.

    https://doi.org/10.4049/jimmunol.177.5.3461

    • PubMed
    • Google Scholar
45. 1. Stamm C
    2. Barthelmann J
    3. Kunz N
    4. Toellner KM
    5. Westermann J
    6. Kalies K

    (2013) Dose-dependent induction of murine Th1/Th2 responses to sheep red blood cells occurs in two steps: antigen presentation during second encounter is decisive

    PLOS ONE 8:e67746.

    https://doi.org/10.1371/journal.pone.0067746

    • PubMed
    • Google Scholar
46. Software

    1. Team I

    (2019) immunarch: An R Package for Painless Bioinformatics Analysis of T-Cell and B-Cell Immune Repertoires, version 0.4.1

    Zenodo.

    https://doi.org/10.5281/zenodo.3383240
47. 1. Teraguchi S
    2. Saputri DS
    3. Llamas-Covarrubias MA
    4. Davila A
    5. Diez D
    6. Nazlica SA
    7. Rozewicki J
    8. Ismanto HS
    9. Wilamowski J
    10. Xie J
    11. Xu Z
    12. Loza-Lopez MJ
    13. van Eerden FJ
    14. Li S
    15. Standley DM

    (2020) Methods for sequence and structural analysis of B and T cell receptor repertoires

    Computational and Structural Biotechnology Journal 18:2000–2011.

    https://doi.org/10.1016/j.csbj.2020.07.008

    • PubMed
    • Google Scholar
48. 1. Textor J
    2. Fähnrich A
    3. Meinhardt M
    4. Tune C
    5. Klein S
    6. Pagel R
    7. König P
    8. Kalies K
    9. Westermann J

    (2018) Deep sequencing reveals transient segregation of T cell repertoires in splenic T cell zones during an immune response

    The Journal of Immunology 201:350–358.

    https://doi.org/10.4049/jimmunol.1800091

    • PubMed
    • Google Scholar
49. 1. Thieme M
    2. Bieber K
    3. Sezin T
    4. Wannick M
    5. Gupta Y
    6. Kalies K
    7. Ludwig RJ
    8. Mousavi S
    9. Zillikens D
    10. Sadik CD

    (2019) The Sphingosine-1-Phosphate receptor modulator fingolimod aggravates murine epidermolysis bullosa Acquisita

    Journal of Investigative Dermatology 139:2381–2384.

    https://doi.org/10.1016/j.jid.2019.03.1159

    • Google Scholar
50. 1. Thomas N
    2. Best K
    3. Cinelli M
    4. Reich-Zeliger S
    5. Gal H
    6. Shifrut E
    7. Madi A
    8. Friedman N
    9. Shawe-Taylor J
    10. Chain B

    (2014) Tracking global changes induced in the CD4 T-cell receptor repertoire by immunization with a complex antigen using short stretches of CDR3 protein sequence

    Bioinformatics 30:3181–3188.

    https://doi.org/10.1093/bioinformatics/btu523

    • PubMed
    • Google Scholar
51. 1. Tubo NJ
    2. Pagán AJ
    3. Taylor JJ
    4. Nelson RW
    5. Linehan JL
    6. Ertelt JM
    7. Huseby ES
    8. Way SS
    9. Jenkins MK

    (2013) Single naive CD4+ T cells from a diverse repertoire produce different effector cell types during infection

    Cell 153:785–796.

    https://doi.org/10.1016/j.cell.2013.04.007

    • PubMed
    • Google Scholar
52. 1. Wittenbrink N
    2. Weber TS
    3. Klein A
    4. Weiser AA
    5. Zuschratter W
    6. Sibila M
    7. Schuchhardt J
    8. Or-Guil M

    (2010) Broad volume distributions indicate nonsynchronized growth and suggest sudden collapses of germinal center B cell populations

    The Journal of Immunology 184:1339–1347.

    https://doi.org/10.4049/jimmunol.0901040

    • Google Scholar
53. 1. Zarnitsyna VI
    2. Evavold BD
    3. Schoettle LN
    4. Blattman JN
    5. Antia R

    (2013) Estimating the diversity, completeness, and cross-reactivity of the T cell repertoire

    Frontiers in Immunology 4:485.

    https://doi.org/10.3389/fimmu.2013.00485

    • PubMed
    • Google Scholar
54. 1. Zhu J
    2. Yamane H
    3. Paul WE

    (2010) Differentiation of effector CD4 T cell populations (*)

    Annual Review of Immunology 28:445–489.

    https://doi.org/10.1146/annurev-immunol-030409-101212

    • PubMed
    • Google Scholar
55. 1. Zvyagin IV
    2. Tsvetkov VO
    3. Chudakov DM
    4. Shugay M

    (2020) An overview of immunoinformatics approaches and databases linking T cell receptor repertoires to their antigen specificity

    Immunogenetics 72:77–84.

    https://doi.org/10.1007/s00251-019-01139-4

    • PubMed
    • Google Scholar

[Editors' note: this paper was reviewed by Review Commons.]

Acceptance summary:

This paper makes use of an elegant and technically complex approach to study the T cell receptor (TCR) clonotype dynamics of follicular helper T (Tfh) cell during an auto-antigen challenge, which is done by laser capturing of germinal centres inside popliteal lymph nodes, and combined with TCR sequencing, clearly isolating differentiated Tfh cells.

Reviewer #1 (Evidence, reproducibility and clarity (Required)):

Niebuhr et al. use a mouse model of the autoimmune disease epidermolysis bullosa acquisita to investigate the T cell receptor repertoire of follicular helper T cells (TFH) isolated from the light zone of germinal centers by laser capture as determined using high throughput sequencing. The key question the authors address is whether such repertoires differ when the same mouse is immunized in two different locations in parallel. Two weeks after immunization dominant T cell clones were comparable in germinal centers in the two lymph nodes that drain the two different injection sites. The size of such clones was contracting at 7 and 10 weeks after injection and similarities of repertoires between the two lymph nodes were lost. These are well executed complex experiments. Welcome are the repertoire characterization of the PBS control, the comparison with effector T cells, and the use of a second mouse model for corroboration. The data support the idea that the immune response to the autoantigens is comparable across the entire mouse rather than distinct in each draining lymph node.

We thank the reviewer for this valuable note.

The compact manuscript lacks causal experiments. The authors could for example prevent lymphocyte trafficking between lymph nodes to distinguish whether the similar TCR repertoires in the two draining lymph nodes two weeks after injection arise from independent parallel priming or from lymphocyte trafficking between the lymph nodes. Such experiments would substantially strengthen the manuscript.

Thank you for raising this point. We agree with the reviewer that is important to examine lymphocyte trafficking. We surgically removed one 2wk lymph node and compared Tfh clonotypes to the 10wk contralateral lymph node of the same mouse, which reflects an inhibition of T cell trafficking between both lymph nodes for a period of 8 wks. The data show that especially the size of the dominant 2wk Tfh clonotypes decreased in the 10 wk lymph nodes and other Tfh clonotypes became dominant but most of the 2 wk Tfh clonotypes were maintained within the GC over this period of 8 wk. We agree with the reviewer that the changes observed indicate an exchange of dominant clones within one lymph node, but it cannot be concluded that Tfh clonotypes exchange between lymph nodes. To clarify this, we rewrote the manuscript substantially (title, abstract, second paragraph on page 7, line 166 and 181).

In addition, we performed experiments and compared T cell receptor sequences in left and right lymph nodes of the same mouse 1d and 3d after priming. The number of shared dominant T cell clonotypes increased 3d p.i.. (Morisita Horn index). We included these data in the discussion (page 11, line 267).

Looking at the week 7 and 10 TFH TCR repertoires, Niebuhr et al. interpret the lack of significant correlation between the repertoires in the two draining lymph nodes as evidence for clonotype replacement (line 157). There may be a simpler explanation that should be considered. Clone sizes contract over the course of the immune response as seems evident in Figure 4A. This should be quantified. As the repertoire of non-antigen-specific follicular helper T cells is different in each lymph node, a contraction of the antigen-specific clones to the size of the non-antigen-specific ones may simply make it technically impossible to follow the antigen-specific clones on the variable background of the non-antigen-specific ones. In this context, the linear relation analysis used by the authors is largely driven by the relatively small number of large clones.

Thank you for raising this basic concern. We agree that our data focus especially on the high frequent clones that are shared between left and right lymph nodes. In our opinion these are the most relevant Tfh clonotypes to look at due to the limited size in GC and the high clonal competition between Tfh for space within GC (Merkenschlager et al. Nature 2021, doi: 10.1038/s41586-021-03187-x). Our data show that the size of the overlapping dominant 2 wk Tfh clones declines over time and the size of other potential bystander activated Tfh clones becomes superior (Figure 4). This does not mean that the initial large clones would be completely replaced. Instead, they are still there but at lower frequencies. We rewrote the manuscript accordingly (title, abstract, second paragraph on page 7, line 166 and 181). In addition, to study the clone distribution we performed box plot analysis from the 2wks and 10wks GC-Tfh clonotypes (Figure 4c and d) that shows that the overall clone size of GC-Tfh does not generally contract over time.

As minor comments, statistics for Figure 1C should be given. Arrows in Figure 1D may be shifted.

We appreciate this comment and added statistics and changed the arrows

accordingly.

Reviewer #1 (Significance (Required)):

The manuscript leaves important questions unresolved.

A similar TCR repertoire in the distinct lymph nodes draining two injection sites could in principle be caused by two different mechanisms. The naïve T cell repertoire in each lymph node could be sufficiently broad that the antigen-driven selection in the germinal centers robustly yields independent groups of similar clones. Alternatively, lymphocytes could travel between lymph nodes thus setting up organism-wide clonal selection. These scenarios can be experimentally distinguished as discussed above.

Thank you for raising this point. The surgical removal of one lymph node resembles an inhibition of T cell trafficking within one mouse. In addition, we performed new experiments that show the synchronization between both lymph nodes starts during initial T cell proliferation 3d p.i.. However, we cannot conclude whether these shared T cell clonotypes would differentiate into GC-Tfh cells. Even though the idea to inhibit lymphocyte trafficking is appealing we feel that this kind of kinetic analysis that would be the required is beyond the scope of this paper. We rewrote the manuscript substantially (title, abstract, second paragraph on page 7, line 166 and 181).

As is, the manuscript describes that TFH TCR repertoires are similar in different draining lymph nodes. Are there any practical consequences to that? There should not be, one injection should be sufficient. Are there any practical consequences to the mechanism ensuring this similarity, as outlined in the first question? If there are, this should be discussed.

Thank you for asking these clarifying questions. To better emphasis the practical consequences of our manuscript we changed the focus towards Tfh clonotypes during the course of an autoimmune disease and the importance of timing. Usually in autoimmune diseases, multiple tissues sites are affected, and it has been suggested that autoreactive clonotypes accumulate in the most adjacent lymphoid organ, which was simple judged by V/J gene segment expression (Oftedahl et al. 2017, doi: 10.1016/j.jaut.2017.03.002). In contrast to this, another report showed that high numbers of bystander-activated Tfh clonotypes accumulate also in autoimmune models (Ritvo et al., 2018, doi: 10.1073/pnas.1808594115). Our data fills the gap by showing that indeed identical Tfh clonotypes accumulate in GC of lymph nodes that drain autoantigen-exposed skin sites but that this high overlap of prominent Tfh clonotypes is only transiently. In addition, our approach has the potential to identify antigen-specific T cell clones in vivo and could help to improve vaccination strategies. To emphasize these practical consequences we rewrote the manuscript profoundly (title, abstract, main text and discussion).

The key expertise of this reviewer is in T cell signal transduction. This review thus presents a more general immunological view from outside of the core question of TFH TCR repertoires in the germinal centers.

Referees cross-commenting

I share the conceptual concerns of reviewer 2. I regard laser capture of individual germinal centers and whole lymph node cytometry as complimentary, each having its strengths and weaknesses. I can't comment of the choice of sequencing approach as discussed by reviewer 3.

Reviewer #2 (Evidence, reproducibility and clarity (Required)):

The manuscript by Niebuhr et al. used laser dissection and TCR sequencing to address the role of follicular T helper (Tfh) cell clonality after immunization in mice. While the topic of the manuscript is of interest, the experimental data and several technical concerns do not fully support the conclusions drawn from the study.

Major points:

1. Mice were immunized s.c. in both hind foot pads with either a model antigen and Titermax adjuvant or with PBS and Titermax adjuvant, i.e. one mouse received either the Ag or the PBS control. They found that 90% of the clones overlapped in the LNs and that this was transient. The authors claim that TCRs failed to remain enriched and conclude from that the selection of germinal center (GC)-Tfh is controlled by systemic clonal competition throughout the response as stated in the abstract.

The connection between the enrichment of clones (which is induced by the expansion of activated T cells, most likely clonal, which is not surprising) and its control by "systemic clonal competition" remains questionable though.

We thank the reviewer for drawing our attention to the lack of clarity. We agree with the reviewer that an accumulation of highly expanded T cell clones after immunizations is not surprising. However, we showed previously that high expanded T cell clones even within adjacent T cell zones of murine spleen do not substantially overlap after antigen exposure (Textor et al. doi: 10.4049/jimmunol.1800091). Therefore, the huge difference in the number of shared Tfh clones in both lymph nodes between the Ag and PBS group is surprising. We used the term “systemic versus local selection” to explain this difference between the PBS and Ag group. The adjuvant in the PBS group induces a local inflammation in both footpads leading to the presentation of multiple local epitopes. Each lymph node develops its local Tfh repertoire that does not overlap with the other lymph node. In contrast, the Tfh repertoire of both lymph nodes overlaps in the Ag group. We called this “systemic” because both lymph nodes are involved. We understand that this wording might be misleading. To clarify this, we rewrote the abstract, and changed it throughout the text.

Conceptually, one big question remains: Why did the authors not use both Ag challenge and control treatment in the same mouse? The foot pad/popl. LNs are ideally suited to test this in one mouse. The presented data further do not allow the conclusion that "the GC response is shared between separate LNs and underlies a constant systemic competition" (line 270/271 and similarly in the abstract). Another point is the timing of 2,4, and 7 weeks, which appears quite late given that GCs (admittingly dependent on the Ag/adjuvant used) often reach their peak magnitude already around 7-14 days post immunization.

Thank you for suggesting this experiment. This kind of experiment should be clearly done. It would be also interesting to inject different antigens and repeat the injections. However, even though it would substantially add interesting information, we feel that this would be beyond the scope of this manuscript and would extend it too much. Additionally, we rewrote substantial parts of the manuscript and avoided terms like “shared GC responses and systemic competition”. Regarding the timing, the time points were chosen depending on the onset of skin pathologies. It takes at least 4 wks for the development of initial skin lesions. This long period is probably required for the loss of tolerance and the recruitment of autoreactive T cell clones. To address earlier time points, we compared the T cell repertoire of complete lymph nodes 1d p.i. and 3d p.i..(discussion page 11, line 267, supplemental Figure 2).

2. The approach of dissecting GCs with laser capture is appealing, however, it would be much easier to use sorting by flow cytometry, as this would not only allow single cell sequencing analysis (which was not performed here, only bulk in Figure 2.) but flow cytometry would also provide a much more powerful analysis of the Tfh cells in this setting. While the authors argue that the observed skin phenotype is depending on high-affinity Abs derived from GCs, GC B cells were not assessed here at all, even though a B cell read-out would further help in dissecting the Tfh differences between Ag/TM vs. PBS/TM conditions.

We appreciate this comment. In this approach we aimed to analyze Tfh in GC without disturbing the organized structure of the lymph nodes. This approach has the advantage that Tfh cell can be analyzed in vivo, which avoids any bias such as unwanted T cell activations or deaths during cell isolation pro. Flow cytometry was used as control. Interestingly, the similarity was higher in FACS-sorted GC-Tfh clonotypes compared to lasercaptured GC-Tfh clonotypes (MHI, 0.46 ± 0.07 versus 0.71 ± 0.11, mean ± SD, Figure 2f). This data indicates that each GC within one lymph node shares the majority of Tfh clonotypes. We included this result in the text (page 6, line 137). We agree with the reviewer that FACS should be the method of choice for further studies. We have not expected before that lasercaptured GC and FACS sorting would give such similar results. Regarding B cell data, we feel that this kind of new experiments would be beyond the scope of this manuscript.

3. Principally, sorting works (Figure 2e-h). Nevertheless, there are some concerns here: A) The gates in Figure 2e appear (hopefully) not correctly positioned. Did it shift upward during the preparation of the figure? B) From the plot and the legend, it seems that all CXCR5+PD-1+ cells were sorted. These cannot be regarded as GC-Tfh cells, since GC-Tfh cells are only those cells with very high expression of both markers.

Thank you very much pointing towards this important issue. We apologize for the

bad quality of this FACS plot. We now added the complete gating strategy (Figure 2e).

4. The authors use very low numbers of replicates for their experiments. While complex assays such as laser capture and subsequent sequencing require a lot of effort, it needs to be made sure that the results are robust and reproducible. Claiming in Figure 1b that there is no difference between left and right pLN if there are only two mice shown and in mouse 2 there are ~35 vs. ~55 GCs in left vs. right LN, underscores this (with n=2 being not enough for statistical use). Same for 1c.

Thank you! We understand the concern of the reviewer, changed Figure 1b and

1c and rewrote this sentences (page 5, lines 97 and 113). It is more important that GC reactions emerge in both pln instead that the number of GC are identical.

5. Lines 117-119: “This data demonstrate that distinct T cell clones contribute to the responses in GC and skins, which is unexpected considering the fact that all T cells were activated with the same Ag.” I would argue that this is actually expected, since those cells that emigrate are most likely not the same cells that end up in GCs, i.e. they are potentially primed at different sites within the T zone or at the T-B border, with different signal strengths, and additionally at different time points (Marc Jenkins’ work). Many variables here. Furthermore, mice were immunized with complex antigens most likely containing several different epitopes.

We appreciate this comment and apologize for this inconclusive conclusion. We

agree with the reviewer and rewrote this paragraph (page 6, lines 133).

6. It appears that during cell sorting by flow cytometry dead cells were not excluded with a viability dye. That should always be done to reduce background/false-positive staining, and is regarded as good scientific practice. This is particularly important for sensitive down-stream applications such as sequencing.

To clarify this issue we included the complete gating strategy (Figure 2e).

7. The authors injected up to 120µg of Ag in 60µl total volume in the hind foot pads. That is in both regards a lot. Is this due to the autoantigen-driven immunization model that would otherwise not precipitate a phenotype? Is the observed data dependent on such high antigen-load? Would the results be different if other exogenous (model)antigens would be used at lower concentrations?

Thank you for raising this interesting questions. This high dosage of antigen is required to induce pathologies in this skin blistering autoimmune model. It will be important in further studies to use different kind of antigens and antigenic peptides.

Minor points

8. Some words should be revised, e.g. "skins" line 118, "maintenances" line 144, "Therefore" line 162, "avoid" line 262, "scarified" line 312.

Thank you. We revised all points accordingly.

9. Line 335, provide a number or rough estimate of the serial sections used for the reconstruction.

Thank you. We included this information on page 15, line 367.

10. What is the time point of the PBS/TM data in Figure 5? Also week 2?

Thank you. The time point for PBS/TM is 4 wks p.i.. To clarify this, we changed Figure 5 and rewrote the respective paragraph (page 8, line 204).

Reviewer #2 (Significance (Required)):

The question of clonality of endogenous Tfh cell response is relevant and open questions remain. The antigen-specificity in this system is not entirely clear, since complex antigens were used. The clonality of the Tfh cells response has been adressed in humans before, e.g. in influenza infection: Brenna et al. Cell Rep 2020 (PMID: 31914381) found that TCR sequences showed similarities between tonsils and the periphery. This study shouldbe be discussed.

Thank you for these comments. We agree with the reviewer that peptides should be used for further analysis. Regarding the antigen-specificity recent data demonstrated that GC-Tfh cells are selected by affinity and TCR receptor signaling strength and underly a high clonal competition for space within GC (Merkenschlager et al. Nature 2021, doi: 10.1038/s41586-021-03187-x). Therefore, especially the dominant Tfh clones should be antigen-specific. In addition, the study of Brenna et al. shows that also circulating Tfh cells are antigen-specific and overlap with GC-Tfh clones in tonsils. Both studies were included into the introduction and discussion (page 3, line 59; page 10, line 231 and page 11, line 262)

Reviewer #3 (Evidence, reproducibility and clarity (Required)):

The article "Receptor repertoires of murine follicular T helper cells reveal shared responses in separate lymph nodes" summarizes the results of a beautiful and technically complex study devoted to the Tfh clonotype dynamics in antigenic response. The general design of the experiment is fascinating. Authors managed to perform both short-term and long-term experiments on littermate mice, performed beautiful work of laser capturing of germinal centers inside popliteal lymph nodes, and combined the technology with TCR sequencing, clearly isolating differentiated Tfh cells. The reviewer would access this study from the point of view of TCR sequencing expert.

Thank you very much for this comment.

Major comments

1. TCR sequencing data analysis is out-of date and not sufficiently justified. Instead of

MiTCR software more recent and customizable MiXCR software should be used. MiTCR was well-known to produce overestimation of diversity in samples, which is detrimental in conjunction with multiplex PCR approach of iRepertoire libraries. At the same time, MiXCR could be optimized for murine samples and multiplex PCR-based DNA libraries and provide adequate error correction. Data engineering is extremely important to check for overamplification in iRepertoire data, filter the cross-contamination between adjacent samples and avoid batch effects such as extensive clonal overlaps in top clonotypes in samples which were sequenced together on MiSeq. Authors should consider re-analysis of raw data with additional attention to batch effects as just considering clonotypes with at least 2 reads as real could be not enough to get a clean dataset. Usually with iRepertoire data you need to look at the suspicious expanded clones in sorted naïve T cells, which represent contamination from adjacent memory samples. With technical and biological replicates you can estimate the level of contamination and discard all clones with this threshold.

We appreciate these critical comments about our TCR analysis approach. However, many new TCR analysis approaches have been developed and computational comparative studies identified MiTCR as solid choice (Azfal et al. 2019, doi:

10.1093/bib/bbx111). Additionally, the higher diversity obtained in multiplex PCR enrichmentbased libraries makes sequencing easier compared to usage of 5’RACE based libraries. In this approach it is wanted and it is not problematic to yield high numbers of T cell clonotypes because:

(1) Due to the number of Tfh cells within GC, the number of Tfh clonotypes is considerable low when compared to entire lymph nodes (see table S1 and table S4). Each GC reaction emerged just in response to Ag or to adjuvant. This excludes the presence of potential false expanded memory T cells among naïve T cells.

(2) We focused only on dominant Tfh clones. All those TCR sequences below the median of the respective sample were removed. Thus, all clonotypes with 2 reads only do not play a role at all.

(3) Our data does not focus on the diversity within one sample instead we compared dominant Tfh clonotypes at two tissue sites. In our opinion, using MiTCR and iRepertoire rather strengthen our data due to the high yields of TCR sequences and our finding that the overall diversity decreases in the Ag exposed groups. In addition, a comparative study of computational methods for T-cell receptor sequencing data revealed that MiTCR and MiXCR do not profoundly differ in clonality measures (Azfal et al. 2019, doi: 10.1093/bib/bbx111). To address this point, we included this reference and wrote a paragraph in the discussion (page 10, line 235).

2. The manuscript contains inconsistent terminology, for example, TRBV segment naming. Axes on plots either lack labels or are incomprehensible in many cases. Therefore, the manuscript needs some editing.

Thank you for this comment! We edited the manuscript accordingly.

3. Authors describe that the time points in this study were selected accordingly to the development of experimental autoimmune process, the dynamics of its phenotype and its pathology – e.g. autoantibody deposition at the epidermis-dermis junction. However, it would be interesting to explore earlier processes of Tfh differentiation, competition for the antigen and establishment of the new GC. Do authors have the data on Tfh TCR dynamics in time points 4-14 days post injection of the antigen? The clonal competition could be more apparent and show if lateral lymph nodes develop the response independently. It may happen that in the first days the competition is restricted to lymph node and therefore private, and later the competition is systemic and both contralateral LN exchange large numbers of T cell clones. However this requires a whole set of experiments and may be out of scope of this particular short manuscript.

Thank you very much for these appealing suggestions. We performed new experiments and included the TCR analysis of contralateral draining lymph nodes 1d and 3d after priming. We found that the number of dominant T cell clonotypes that are shared between both lymph nodes increased 3d p.i.. (Morisita-Horn-index, Supplemental Figure S2). However, before this at 1d p.i. the number of shared clones between both pln decreased significantly (Jaccard Index, naïve: 0.082 ± 0.02, 1d: 0.053 ± 0.015 (p<0.05) and 3d 0.12 ± 0.01 (p<0.01), n = 3, 2 pln each, Mann-Whitney-U test). It might be that this decrease reflects the recruitment of new naïve T cells into the Ag-exposed lymph nodes before the Agactivated clones start to expand (Figure S2a). Subsequently, those clones that bound with highest affinity and therefore produced the highest number of progenies could distribute between both pln before they differentiate into Tfh cells. We did not include this data into the manuscript because further experiments will be required to answer this question.

Minor comments

1. Authors need to choose 1 strategy of naming and keep it consistent. Mice groups are named Ag and control group, Ag and PBS group and, finally, + and – groups in Figure 3C, while meaning the same groups.

Thank you! We changed our manuscript accordingly.

2. TRBV segment names are used in 3 different ways. Authors need to choose 1 strategy of naming and keep it consistent, for example, keep IMGT recommended gene segment names.

Thank you! We changed our manuscript accordingly.

3. Figure 2A should have clonotype frequency annotated on axis label, r2 coefficient labeled on plot. Figure 3A, 4A, 4C – same comments.

Thank you! We changed the legends accordingly.

4. Figure 2C should have axes labeled. Effector T cells – why are they derived from ear? Text lacks the explanation; it only could be found deep in methods section. Labeling of axes should be more consistent.

Thank you! We changed the legends accordingly and included an explanation into the manuscript on page 6, line 124.

5. Figure 2A,2C – it's not clear from figure legend, what exactly does this percentage mean. Weighted by frequency of clonotype or not. Percentage of top N clones, of total repertoire volume?

Data are weighted by the frequency of the clonotypes that were used for analysis and includes all clonotypes with frequencies above the median (see table S1-S5).

We changed legends accordingly.

6. Figure 2b. If n=2, authors should show both mice, not one representative for both animals.

We changed Figure 1b accordingly and rewrote this paragraph on page 5, lines

97-100.

7. Page 6. "TCRbeta seq directly from GC" – what does it mean? Is there a separate step of nucleic acid isolation and library preparation?

Directly from GC means that we did not disrupted the lymph nodes into a cell suspension before RNA isolation and library preparation. To clarify this, we rewrote this sentence on page 6, line 110.

8. Page 6. 1 million of CDR3-containing reads: this number is collected from 1 GC or from the whole lymph node?

The 1 mio TCR sequences were obtained from 4-6 GCs collected by laser microdissection from histological sections to yield sufficient amount of RNA. This information is described on page 5 line 107.

9. Page 6. "This data shows that GC reactions occur simultaneously at two separate lymph nodes independent whether Ag is present or not." Authors should consider moving this phrase forward to the Discussion section or removing it as it is not justified by the data in this paragraph.

To clarify this data, we rewrote this phrase but kept it at the end of this paragraph on page 5 line 113. The histological data clearly show that GC reactions take place in both lymph nodes and in both, the Ag1 and the PBS group.

10. In the literature analysis, authors have missed recent important research in the field of Tfh TCR repertoires, notably https://pubmed.ncbi.nlm.nih.gov/31914381/.

We added this reference in the introduction on page 3, line 59 and in the discussion on page 11, line 262.

11. Figure 1C The phrase "% area of GC of B cell follicles area" is unclear at least for nonprofessionals in imaging: first the area of B cell zone is calculated, and then the percentage of GC?

The cumulative area of the complete B cell zone of one cryosection was determined. Likewise, the cumulative area of all germinal centers of this cryosection was determined. The relative proportion of the germinal center area within the B cell zone area was calculated as quotient. For each sample three individual cryosection were evaluated and the mean proportion was used for evaluation. To clarify this, we included this explanation into the Material and Method section (page 16, line 346) and in the legend of Figure 1.

12. “This data demonstrate that distinct T cell clones contribute to the responses in GC and skins, which is unexpected considering the fact that all T cells were activated with the same Ag.” Why authors imply the unexpected result if the consensus is that Teff and Tfh frequently recognize different set of antigens? Again, this part of discussion may benefit from more recent works on Tfh TCR repertoires.

We apologize for this lack of clarity. We rewrote this paragraph on page 6, line 133

13. Page 7. The similarity of clonotypes is not defined clearly. May be replaced with "sequence overlaps in Tfh TCR repertoires" or "shared identical sequences", or defined through edit distance. It would be useful to see clonotypes tables of all samples to compute different metrics of repertoire pairwise distance. Could authors provide the tables or the raw sequencing data in open access, e.g., through GEO dataset?

Thank you for raising this point. We replaced the phrase “similarity” accordingly. Of course, the data sets will be published in the Sequence Read Archive as done for previous work (bioproject PRJNA586880).

14. Page 9. The usage of sharing metrics in repertoire comparison should be explained in more details. Why a certain sharing index was used, why another edit distance measure is not applicable? Repertoire sharing strongly depends on the size of repertoire, therefore the data cleaning and downsampling needs to be described in details. 18 shared clonotypes and 65 shared clonotypes: are these nucleotide sequences or aminoacid CDR3 sequences? Do they have the same pattern of recombination – is there a match in both TRBV and TRBJ segments used? It's 18 and 65 of 500 clonotypes? Of n clonotypes if sizes of compared repertoires differ?

We apologize for the lack of clarity. We used the term “sharing index” not as a special metrics but just as a simple enumeration of shared Tfh clonotypes between the mice either of the Ag group and or the PBS group as described by Madi et al. 2014, doi:

10.1101/gr.170753.113. First, all the Tfh clonotypes that were detected once per mouse were assessed by combining all Tfh clonotypes of both contralateral lymph nodes and removing all duplicates. These number of unique Tfh clonotypes were combined from all three mice in each group. Thus, 4 wks after immunization, 7449 and 6445 Tfh clonotypes were unique in all three mice of the Ag1 group or PBS group, respectively. From those, 424 (Ag1) or 244 (PBS) Tfh clonotypes were present in 2 mice and 65 versus 18 in 3 mice of the Ag1 group or PBS group. This data demonstrates that Ag-exposure increases the number of the shared Tfh clonotypes between all mice. For a better clarification we rewrote this paragraph on page 8, line 204 and removed the 2 wks and 7 wks groups on Figure 5. For all other samples, the Morisita-Horn-index were used for quantifying shared clonotypes between the groups. This index considers not only the quantity of clonotypes within one sample but also the frequency with which each clonotype exist. We included an explanation on page 6 line128.

15. Figure 5. There are no significant changes on the plot, why the name of the figure strongly claims changes over time? Figure 5B. It's not convincing without statistics that TRBV3 is the most prominent gene family which changes frequency after Ag challenge. Moreover, it is certainly not enough without TCR transfection experiments to claim that this TRBV is specific to the Ag1-derived epitopes. These claims should be removed, or additional TCR cloning and specificity confirmation experiments are required.

Changes in the V/J usage in autoimmune models have been described from other groups before without transfection experiments. In this Figure, we demonstrate that among the Tfh clonotypes that are shared between all three mice especially the percentage of those that bears the TRBV3 gene segment accumulated in the Ag1 group compared to the PBS group (from 16.6% to 33.84%). To make this point clearer, we removed the 2wks and 7wks time points as done for Figure 5A. Due to the fact that only the shared Tfh clonotypes are analyzed statistical analysis is not applicable. Accordingly, we changed our wording and rewrote this paragraph on page 8, line 214.

16. Page 16. "the repertoire of TCR β clonotypes was identified as described above" Should be below instead of above?

Thank you for pointing that out! We changed the wording accordingly.

Reviewer #3 (Significance (Required)):

The article "Receptor repertoires of murine follicular T helper cells reveal shared responses in separate lymph nodes" summarizes the results of a beautiful and technically complex study. However, reviewer argues that both the experimental data, the methods, and result presentation need substantial revision and improvement to be considered in a peer-reviewed journal.

Reviewer's expertise covers T cell biology, TCRseq data analysis in both human and rodents' immune repertoires datasets, as well as features of different library preparation protocols and sequencing platforms. Unfortunately, TCRseq data analysis in this particular study doesn't stand up to the modern standards and recommendations of AIRR community. Histology and related analysis sections are out of the scope of the reviewer's expertise.

To find out whether the high overlap of Tfh clonotypes in separate pln is biased by the TCR analysis tool MiTCR, which might overestimate TCR repertoire diversities, we reanalyzed the raw data of some key samples with MiXCR and Immunarch according to the reviewer suggestions. Data reveal an almost identical distribution of Tfh clonotypes between both pln regardless of whether MiTCR and MiXCR were used for analysis. Results are included as Figure 2-figure supplement 1 in the discussion section on page 11 and 12, lines 245-259. In addition, the sharing index, which displays the number of identical (public) Tfh clonotypes between the mice of one group (Figure 5) is now complemented with venn diagrams in supplementary Figure 2-figure supplement 2. This information is included in the result section page 10, lines 211-214 and in the Material and Method section page 19 lines 445-446.

Author details

1. Markus Niebuhr

   Institute for Anatomy, University of Lübeck, Lübeck, Germany

   Contribution

   Data curation, Software, Formal analysis, Validation, Investigation, Methodology

   Competing interests

   No competing interests declared

2. Julia Belde

   Institute for Anatomy, University of Lübeck, Lübeck, Germany

   Contribution

   Data curation, Investigation

   Competing interests

   No competing interests declared

3. Anke Fähnrich

   1. Institute for Anatomy, University of Lübeck, Lübeck, Germany
   2. Lübeck Institute of Experimental Dermatology, University of Lübeck, Lübeck, Germany

   Contribution

   Data curation, Formal analysis, Validation, Methodology

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0003-1904-9544

4. Arnauld Serge

   Laboratoire Adhésion et Inflammation, Inserm U1067 CNRS, Aix-Marseille Université, Marseille, France

   Contribution

   Data curation, Methodology

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0002-1280-6277

5. Magali Irla

   Centre d'Immunologie de Marseille Luminy (CIML), INSERM U1104, Aix-Marseille Université UM2, Marseille, France

   Contribution

   Data curation, Investigation, Methodology

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0001-8803-9708

6. Christoph T Ellebrecht

   1. Institute for Anatomy, University of Lübeck, Lübeck, Germany
   2. Department of Dermatology, University of Pennsylvania, Philadelphia, United States

   Contribution

   Data curation, Methodology

   Competing interests

   No competing interests declared

7. Christoph M Hammers

   1. Institute for Anatomy, University of Lübeck, Lübeck, Germany
   2. Department of Dermatology, University of Lübeck, Lübeck, Germany

   Contribution

   Data curation, Investigation

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0001-5631-2415

8. Katja Bieber

   Lübeck Institute of Experimental Dermatology, University of Lübeck, Lübeck, Germany

   Contribution

   Data curation, Investigation

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0002-3855-6683

9. Jürgen Westermann

   Institute for Anatomy, University of Lübeck, Lübeck, Germany

   Contribution

   Resources, Funding acquisition

   Competing interests

   No competing interests declared

10. Kathrin Kalies

    Institute for Anatomy, University of Lübeck, Lübeck, Germany

    Contribution

    Conceptualization, Supervision, Funding acquisition, Investigation, Methodology, Writing - original draft, Project administration, Writing - review and editing

    For correspondence

    kathrin.kalies@uni-luebeck.de

    Competing interests

    No competing interests declared

    "This ORCID iD identifies the author of this article:" 0000-0002-8929-4249

• 797

  Page views

• 81

  Downloads

• 3

  Citations

Article citation count generated by polling the highest count across the following sources: PubMed Central, Crossref, Scopus.