<figure><img src="image_1.png" /></figure>

HAL
open science

# Immune fingerprinting through repertoire similarity

Thomas Dupic, Meriem Bensouda Koraichi, Anastasia Minervina, Mikhail Pogorelyy, Thierry Mora, Aleksandra Walczak

<figure><img src="image_2.png" /></figure>

► To cite this version:

Thomas Dupic, Meriem Bensouda Koraichi, Anastasia Minervina, Mikhail Pogorelyy, Thierry Mora, et al.. Immune fingerprinting through repertoire similarity. PLoS Genetics, 2021, 17 (1), pp.e1009301. 10.1371/journal.pgen.1009301. hal-03367815

HAL Id: hal-03367815

https://hal.science/hal-03367815

Submitted on 7 Oct 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.
RESEARCH ARTICLE

# Immune fingerprinting through repertoire similarity

Thomas Dupic<sup>1,2</sup>, Meriem Bensouda Koraichi<sup>2</sup>, Anastasia A. Minervina<sup>3</sup>, Mikhail V. Pogorelyy<sup>3,4</sup>, Thierry Mora<sup>2,*</sup>, Aleksandra M. Walczak<sup>2,*</sup>

¹ Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, Massachusetts, USA,

² Laboratoire de physique de l'École Normale Supérieure, CNRS, Sorbonne Université, Université de Paris, and École normale supérieure (PSL), Paris, France, ³ Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry, Moscow, Russia, ⁴ Pirogov Russian National Research Medical University, Moscow, Russia

* These authors contributed equally to this work.

* tmora@phys.ens.fr (TM); awalczak@phys.ens.fr (AMW)

<figure><img src="image_3.png" /></figure>

## Abstract

Immune repertoires provide a unique fingerprint reflecting the immune history of individuals, with potential applications in precision medicine. However, the question of how personal that information is and how it can be used to identify individuals has not been explored. Here, we show that individuals can be uniquely identified from repertoires of just a few thousands lymphocytes. We present “Immprint,” a classifier using an information-theoretic measure of repertoire similarity to distinguish pairs of repertoire samples coming from the same versus different individuals. Using published T-cell receptor repertoires and statistical modeling, we tested its ability to identify individuals with great accuracy, including identical twins, by computing false positive and false negative rates < 10⁻⁶ from samples composed of 10,000 T-cells. We verified through longitudinal datasets that the method is robust to acute infections and that the immune fingerprint is stable for at least three years. These results emphasize the private and personal nature of repertoire data.

<figure><img src="image_4.png" /></figure>

OPEN ACCESS

**Citation:** Dupic T, Bensouda Koraichi M, Minervina AA, Pogorelyy MV, Mora T. Walczak AM (2021) Immune fingerprinting through repertoire similarity. PLoS Genet 17(1): e1009301. https://doi.org/10.1371/journal.pgen.1009301

**Editor:** Derry C. Roopenian, The Jackson Laboratory, UNITED STATES

**Received:** July 13, 2020

**Accepted:** December 7, 2020

**Published:** January 4, 2021

**Copyright:** © 2021 Dupic et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting information files except for the sequence data: the Emerson et al. dataset is available from ImmuneAccess (https://clients.adaptivebiotech.com) with DOI 10.21417/B7001Z. The Chu et al. dataset is also available from ImmuneAccess, with DOI 10.21417/B7J01X; the Pogorely et al. dataset can be found on the Sequence Read Archive (SRA: https://www.ncbi.nlm.nih.gov/sra), under the accession number PRJNA493983; the Briney et al. dataset can also be found on the SRA, under the accession number PRJNA406949; and lastly the Britanova et al.

## Author summary

Immune repertoires are a trove of personal information: unique to each individual, they are also windows into their past and future health. Thanks to their potential for personalized medicine and progress of sequencing technologies, these repertoires are now routinely sequenced. As a consequence they raise the question of identifiability of samples. In this paper, we estimate the quantity of immune cells needed to associate two samples from the same individual: as little as a finger prick worth of blood can serve as an immune fingerprint that can distinguish even identical twins, without giving away genetic information about non-consenting relatives. We show that this fingerprint is stable through time, and is not erased during infections or vaccinations.
dataset can be extracted from the SRA under
accession number PRJNA316572.

**Funding:** The study was supported by the European Research Council COG 724208 (AMW). A.A.M and M.V.P. are supported by RSF-20-15-00351, Russian Science Foundation. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing interests:** The authors have declared that no competing interests exist.

# 1 Introduction

Personalized medicine is a frequent promise of next-generation sequencing. These high-throughput and low-cost sequencing technologies hold the potential of tailored treatment for each individual. However, progress comes with privacy concerns. Genome sequences cannot be anonymized: a genetic fingerprint is in itself enough to fully identify an individual, with the rare exception of monozygotic twins. The privacy risks brought by these pseudonymized genomes have been highlighted by multiple studies [1–3], and the approach is now routinely used by law enforcement. Sequencing experiments that focus on a limited number of expressed genes should be less prone to these concerns. However, as we will show, B- and T-cell receptor (BCR and TCR) genes are an exception to this rule.

BCR and TCR are randomly generated through somatic recombination [4], and the fate of each B- or T-cell clone depends on the environment and immune history. The immune T-cell repertoire, defined as the set of TCR expressed in an individual, has been hailed a faithful, personalized medical record, and repertoire sequencing (RepSeq) as a potential tool of choice in personalized medicine [5–9]. In this report, we describe how, from small quantities of blood (blood spot or heel prick), one can extract enough information to uniquely identify an individual, providing an immune fingerprint. The “Immprint” classifier analyzes this immune fingerprint to decide whether two samples were sampled from the same individual.

# 2 Results

Given two samples of peripheral blood containing respectively $M_1$ and $M_2$ T-cells, we want to distinguish between two hypothetical scenarios: either the two samples come from the same individual (“autologous” scenario), or they were obtained from two different individuals (“heterologous” scenario), see Fig 1A.

TCR are formed by two protein chains $\alpha$ and $\beta$. They each present a region of high somatic variability, labeled CDR3$\alpha$ and CDR3$\beta$, randomly generated during the recombination process. These regions are coded by short sequences (around 50 nucleotides), which are captured by RepSeq experiments. The two chains are usually not sequenced together so that the pairing information between $\alpha$ and $\beta$ is lost. Most experiments focus on the $\beta$ chain, and when not otherwise specified, the term “receptor sequence” in this paper will refer only to the nucleotide sequence of the TRB gene coding for this $\beta$ chain (which include CDR3$\beta$). Similarly, as most cells expressing the same beta chain are clonally related, we will be using the terms “clone” and “clonotype” to refer to set of cells with the same nucleotide TRB sequence, even if they were produced in separate generation events and are not real biological clones (since we have no means of distinguishing the two cases). CDR3$\beta$ sequences are very diverse, with more than $10^{40}$ possible sequences [10]. For comparison, the TCR$\beta$ repertoire of a given individual is composed of $10^8$ to $10^{10}$ unique clonotypes [11, 12]. As a result, most of the sequences found in a repertoire are “private”.

## 2.1 Immprint scores

To discriminate between the autologous and heterologous scenarios, one can simply count the number of unique nucleotide receptor sequences shared between the two samples, which we call $S$. Samples coming from the same individual should have more receptors in common because T-cells are organized in clones of cells carrying the same TCR. By contrast, $S$ should be low in pairs of samples from different individuals, in which sharing is due to rare convergent recombinations. Appropriately setting a threshold to jointly minimize the rates of false positives and false negatives (Fig 1B), we can use $S$ as a classifier to distinguish autologous
<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /><figcaption>**Fig 1.** A) The two samples A and B can either originate from the same individual (autologous) or two different individuals (heterologous). In both scenarios, sequences can be shared between the two samples, but their quantity and quality vary. B) Schematic representation of the distribution of the *S* or *I* scores for multiple pairs of samples extracted from the same individual (in the autologous scenario) or the same pair of individuals (heterologous). The dashed vertical line represents the threshold value. C) Expected value of *S* and *I* for different pairs of samples, sampled from the same individual (in blue) or different ones (yellow). Green dots represent samples extracted from pairs of identical twins. The dashed lines represents the theoretical upper bound for heterologous repertoires (see Methods) for both *S* and *I* (γ = 12). D) Sampling distributions of *S* for 6 different patients (autologous case, each blue curve is one patient) or 6 different pairs of patients (heterologous, each yellow or green curve is a pair of patients) for *M* = 5000. The y-axis scale on the left is adapted to the heterologous distributions while the scale on the right corresponds to the (much wider) autologous ones. The 3 sampling distributions in green correspond to a pair of samples extracted from identical twins. E) Detection Error Trade-off (DET) graph for both summary statistics and different sample sizes *M*. *I* (γ = 12) outperforms *S* in all scenarios. F) AUROC (Area Under Receiver Operating Characteristic), as a function of *M*. The AUROC is a traditional measure of the quality of a binary classifier (a score closer to one indicates a better classifier). The results are shown for *S* and *I* both in the default case (only the *β* chain considered) or for the full (*α*—*β*) receptor.</figcaption></figure>

https://doi.org/10.1371/journal.pgen.1009301.g001
from heterologous samples. S is not normalized for sequencing depth and values of S should not be compared between samples of different size.

The S score can be improved upon by exploiting the fact that some receptors are much more likely than others to be generated during V(D)J-recombination, with variations in generation probability ($P_{gen}$, [13–15]) spanning 15 orders of magnitude. Public sequences (with high $P_{gen}$) are likely to be found in multiple individuals [16], while rare sequences (low $P_{gen}$) are unlikely to be shared by different individuals, and thus provide strong evidence for the autologous scenario when found in both samples. To account for this information, we define the score:

$$ \mathcal{I} = \sum_{\text{shared } s} [\ln(1/P_{\text{gen}}(s)) - \gamma], \quad (1) $$

which accounts for Shannon’s “surprise” $\ln(1/P_{gen})$—a measure of unexpectedness—associated with each shared sequence s, so that rare shared sequences count more than public ones. The constant $\gamma$ depends on the repertoire’s clonal structure and is set to 12 in the following (see Methods for an information-theoretic derivation). $P_{gen}$ is computed using models previously trained on data from multiple individuals [14]. Small differences reported between the $P_{gen}$ of distinct individuals justify the use of a universal model [15].

## 2.2 Measuring error rates

We tested the classifiers based on the S and I scores on TCRβ RepSeq datasets from 656 individuals [17]—labeled according to their cytomegalovirus (CMV) serological status. Sequences were downsampled to mimic experiments where $M_1 = M_2 = M$ cells were sampled and their receptors sequenced. The frequency of a particular clonotype in the sample (the proportion of cells expressing a particular beta chain) was estimated using the read counts of unique TCRβ sequences, and the mean values of S and I computed with a procedure designed to correct for the limited diversity of the sampled repertoire relative to the full repertoire, see Methods 4.3. Similar results may be obtained when $M_1$ and $M_2$ are different (see Methods). The clones most often shared between two autologous samples are also the most clonally expanded—and hence are probably antigen experienced. We verified that the sequence statistics of those expanded clonotypes did not differ from generic ones (S1 Fig).

In Fig 1C, we plot the mean value of S (over many draws of M = 5000 cells) for each individual (autologous scenario, in blue) and between pairs of different individuals (heterologous scenario, in yellow). The two scenarios are clearly discernable under both scores. This result holds for pairs of monozygotic twins obtained from a distinct dataset [18] (green dots), consistent with previous reports that twins differ almost as much in their repertoires as unrelated individuals [18–20]. Heterologous scores (yellow dots) vary little, and may be bounded from above by a theoretical prediction (dashed line) based on a model of recombination [21] (see Methods). On the other hand, autologous scores (blue dots) show several orders of magnitude of variability across individuals. These variations stem from the clonal structure of the repertoire, and correlate with measures of diversity (S2 Fig), which is known to vary a lot between individuals and correlates with age [22], serological status, and infectious disease history [23, 24]. To explore the worst case scenario of discriminability, hereafter we will focus on the individual with the lowest autologous S found in the dataset.

The sampling process introduces an additional source of variability within each individual. Two samples of blood from the same individual do not contain the exact same receptors, and the values of S and I is expected to vary between replicates. Examples of these variations are shown in Fig 1D. The blue (respectively yellow) curves correspond to the sample distributions
in the autologous (heterologous) scenario for different individuals (pairs of individuals). The distribution of *S* is well-approximated by a Poisson distribution, while *I* follows approximately a compound distribution of a normal and Poisson distributions (see Methods for details). Armed with these statistical models of variations, we can predict upper bounds for the false negative and false positive rates. As seen from the detection error trade-off (DET) graph Fig 1E, the Immprint classifier performs very well for a few thousand receptors with an advantage for *I*.

With 10,000 cells, corresponding to ~10 µL of blood, Immprint may simultaneously achieve a false positive rate of < 10<sup>-16</sup> and false negative rate of < 10<sup>-6</sup>, allowing for the near-certain identification of an individual based on the *I* score in pairwise comparisons against the world population ~10<sup>10</sup>. When a large reference repertoire has been collected (M<sub>1</sub> = 1, 000, 000, corresponding to ~1mL of blood), an individual can be identified with just 100 cells (S3 Fig).

The AUROC estimator (Area Under the Curve of the Receiver Operating Characteristic), a typical measure of a binary classifier performance, can be used to score the quality of the classifier with a number between 0.5 (chance) and 1 (perfect classification). The *I* score outperforms the *S* score (Fig 1F), particularly above moderate sample sizes (M ≈ 5000). Both scores can be readily generalized to the case of paired receptors, if the pairing of the two chains is available through single-cell sequencing [25–27] or computational pairing [28], we can consider the sequence of the full receptor TCRαβ, instead of just TCRβ. The generation probability of the full TCR is then given by P<sub>gen</sub>(α, β) = P<sub>gen</sub>(α) × P<sub>gen</sub>(β)[29]. Because coincidental sharing of both chains is substantially rarer than with the β chain alone, using the paired chain information greatly improves the classifier.

While this paper focuses on T-cells and TCR sequences, the structure of the B-cell receptors (BCR) repertoire is very similar to the TCR repertoire and we expect to find qualitatively similar results. As an example we use the dataset obtained in Ref. [30] to measure *S* and *I* for IGH chains (forming half of the BCR receptor), S4 Fig. We see that 5000 IgG+ memory B-cells are enough to reliably identify the individuals in the study. However, B-cells are 5 times less common in peripheral blood than T-cells, and somatic hypermutation tends to distort the statistics of receptors, reducing the reliability of our classifier. Hence, for practical applications, T-cells are better means of identification.

## 2.3 Evolution with time

The previous results used samples obtained at the same time. However, immune repertoires are not static: interaction with pathogens and natural aging modify their composition. The evolution of clonal frequencies will decrease Immprint's reliability with time, especially if the individual has experienced immune challenges in the meantime.

To study the effect of short-term infections, we analyzed an experiment where 6 individuals were vaccinated with the yellow fever vaccine, which is regarded as a good model of acute infection, and their immune system was monitored regularly through blood draws [18]. We observe an only moderate drop in *I* caused by vaccination (Fig 2A).

This is consistent with the fact that infections lead to the strong expansion of only a limited number of clones, while the rest of the immune system stays stable [31–34]. While other types of infections, auto-immune diseases, and cancers may affect Immprint in more substantial ways, our result suggests that it is relatively robust to changes in condition.

We then asked how stable Immprint is over long times. Using longitudinal datasets [18, 22, 35], we show in Fig 2 that the Immprint score *I* is only slightly reduced for samples collected at intervals of up to three years. For longer timescale, given the lack of longitudinal datasets,
<figure><img src="image_7.png" /><figcaption>**Fig 2.** A) Evolution of *I* (M = 5000) during vaccination, between a sample taken at day 0 (vaccination date) and at a later timepoint. Each color represents a different individual. Each pair timepoint/individual has two biological replicates. The dashed line represents the threshold value. B) Evolution of *I* between a sample taken at year 0 and a later timepoint. The red histogram corresponds to one of the individuals sampled in [18] and the blue curves show theoretical estimates, fitted to match (τ = 0.66). C) Evolution of the (normalized) mean of *I* (M = 5000) as a function of time for different values of the turnover rate τ. The dashed line represents the threshold value divided by the smallest value of *I*<sub>t</sub> = 0 (M = 5000) in the data. The data points were obtained from the datasets [35] (yellow), [18] (green) and [22] (orange). Different markers indicate different individuals.</figcaption></figure>

https://doi.org/10.1371/journal.pgen.1009301.g002

we turn to mathematical models [12, 36–39] to describe the dynamics of the repertoire. Following the model of fluctuating growth rate described in Ref. [38], we define two typical evolutionary timescales for the immune system: τ, the typical turnover rate of T-cell clones, and θ, which represents the typical time for a clonotype to grow or shrink by a factor two as its growth rate fluctuates. The model predicts a power-law distribution for the clone-size distribution, with exponent −1 − τ/2θ. This exponent has been experimentally measured to be ≈ −2 [10], which leaves us with a single parameter τ, and θ = τ/2. An example of simulated evolution of Immprint with time is shown in Fig 2B The highlighted histogram represents a data point at two years obtained from [18]. While a fit is possible for this specific individual (τ = 0.66 in Fig 2B), the τ parameter is not universal, and we expect it to vary between individuals, especially as a function of age. In Fig 2C we explore the effect on the stability of Immprint for a range of reasonable values for the clone turn-over rate τ, from 6 months to 10 years, encompassing both previous estimates of the parameter [38] and measured turnover rates for different types of T-cells [40]. While Fig 2 focuses on *I*, the behaviour is similar for *S* (see S5 Fig). We observe that under this model, for most individuals and bar exceptional events, Immprint should conserve its accuracy for years or even decades.
## 3 Discussion

In summary, the T-cells present in small blood samples provide a somatic and long-lived bar-code of human individuality, which is robust to immune challenges and stable over time. While the uniqueness of the repertoire was a well known fact, we demonstrated that the most common T-cells clones are still diverse enough to uniquely define an individual and frequent enough to be reliably sampled multiple times. Unlike genome sequencing, repertoire sequencing can discriminate monozygotic twins with the same accuracy as unrelated individuals. Additionally, a person's unique immune fingerprint can be completely wiped out by a hematopoietic stem cell transplant [41].

A potential complication in applying Immprint is the convergent evolution of repertoires: individuals who encounter similar pathogenic environments could share many receptors. While this phenomenon occurs [17, 42], its influence on the immune repertoire is low. For example, in the context of cytomegalovirus infection, shared TCR clones are only slightly more common in co-infected individuals [17], and the result of Immprint does not seem to be affected (S6 Fig). The possibility to discriminate between twins—who shared a common environment for part of their lives—also hints that in most cases the effects of environment-driven convergence is small. Nonetheless we cannot reject the possibility that this effect is stronger for some specific pathogens, or long and strongly overlapping infection histories with pathogens that severely modify immune repertoires. A limit case study to quantitatively investigate this effect would require looking at data from mice living in otherwise sterile environments that are exposed in a controlled way to the same pathogens at the same time, bearing in mind that the diversity of mice repertoires is smaller than that of humans.

The different datasets used cover a range of different sequencing methods (see 4.1), but different approaches may lead to slightly different threshold choices. In particular, in practical implementations, sequencing depth is an important concern. One needs enough coverage to sequence TCRβ genes from as many as possible of the T-cells present in the sample, in order to measure a more precise immune fingerprint. In addition, the specific calculations presented here only apply to peripheral blood cells. Specific cell types or cells extracted from tissue samples may have different clonal distributions and potentially different receptor statistics. For example the value of *S* in the autologous case varies between CD4+ and CD8+ T-cells (S7 Fig), although different individuals remain distinguishable using each subset.

Immprint is implemented in a python package and webapp (see Methods) allowing the user to determine the autologous or heterologous origin of a pair of repertoires. Beyond identifying individuals, the tool could be used to check for contamination or labelling errors between samples containing TCR information. The repertoire information used by Immprint can be garnered not only from RepSeq experiments, but also from RNA-Seq experiments, which contain thousands of immune receptor transcripts [43, 44]. Relatively small samples of immune repertoires are enough to uniquely identify an individual even among twins, with potential forensics applications. At the same time, unlike genetic data from genomic or mRNA sequencing, Immprint provides no information about kin relationships, very much like classical fingerprints, and avoids privacy concerns about disclosing genetic information shared with non consenting relatives.

## 4 Methods

### 4.1 Datasets & pre-processing

We use five independent RepSeq datasets in this study: (i) genomic DNA from Peripheral blood mononuclear cells (PBMCs) from 656 healthy donors [17]; (ii) cDNA of PBMCs
sampled from three pairs of twins, before and after a yellow-fever vaccination [18]; (iii), (iv) two longitudinal studies of healthy adults [22, 35];(v) cDNA dataset of IGH genes (B-cells) from 9 individuals, with multiple replicates [30].

CDR3 nucleotide sequences were extracted with MIGEC [45] (for the second dataset) coupled with MiXCR [46]. We also extract the frequency of reads from the three datasets. The non-productive sequences were discarded (out-of-frame, non-functional V gene, or presence of a stop codon). The generation probability ($P_{gen}$) was computed using OLGA [47], with the default TCRβ model. The frequency of each clone was estimated by summing the frequencies of all reads that shared the same nucleotide CDR3 sequence and identical V,J genes.

The preprocessing code is distributed on the Git repository associated with the paper. We also developed a command-line tool (https://github.com/statbiophys/immprint) that discriminates between sample origins, and a companion webapp (https://immprint.herokuapp.com).

## 4.2 Discrimination scores

To discriminate between the autologous and heterologous scenarios, we introduce a log-likelihood ratio test between the two possibilities:

$$ \mathcal{I} = \sum_s \ln \frac{P(y_1(s), y_2(s)|\text{autologous})}{P(y_1(s), y_2(s)|\text{heterologous})}, \quad (2) $$

where $y_1(s) = 1$ if the sequence $s$ is found in sample 1, and 0 otherwise; likewise $y_2(s) = 1$ if $s$ is in sample 2. The sum runs over all potential sequences $s$, including unseen ones. To be present in a sample, a sequence $s$ first has to be present in the repertoire. This occurs with probability $1 - (1 - p(s))^{N_c}$, where $N_c$ is the total number of clonotypes in the repertoire, and $p(s)$ is the probability of occurrence of sequence $s$ (resulting from generation and selection, see below). Second, it must be picked in a sample of size $M$, with probability $1 - (1 - f)^M \approx Mf$ (assuming $Mf \ll 1$) depending on its frequency $f$, which is distributed according to the clone size distribution $\rho(f)$. We checked that $f(s)$ and $P_{gen}(s)$ were not correlated (S1 Fig) and that the effects of a shared infection between different individuals were limited to a handful of clones (S8 Fig). Then one can write

$$ P(y_1(s) = 1, y_2(s) = 1 | \text{autologous}) \approx (1 - e^{-N_c p(s)}) M_1 M_2 \int df \rho(f) f^2, \quad (3) $$

$$ P(y_1(s) = 1, y_2(s) = 0 | \text{autologous}) \approx (1 - e^{-N_c p(s)}) \frac{M_1}{N_c} \quad \text{and} \quad 1 \leftrightarrow 2, \quad (4) $$

$$ P(y_1(s) = 0, y_2(s) = 0 | \text{autologous}) \approx 1 - (1 - e^{-N_c p(s)}) \frac{M_1 + M_2}{N_c}, \quad (5) $$

where we've used $\int df \rho(f) f = 1/N_c$. For the heterologous case the probability factorizes as:

$$ P(y_1(s), y_2(s) | \text{heterologous}) = P_1(y_1(s)) P_2(y_2(s)), \quad (6) $$

with

$$ P_a(y_a(s) = 1) \approx (1 - e^{-N_c p(s)}) \frac{M_a}{N_c}, \quad a = 1, 2. \quad (7) $$
Since only the term $\gamma_1(s) = \gamma_2(s) = 1$ (shared sequences) is different between the autologous and heterologous cases, we obtain:

$$ \mathcal{I} = \sum_{\text{shared } s} [\ln(N_c^2 \langle f^2 \rangle) - \ln(1 - e^{-N_c p(s)})]. \qquad (8) $$

Further assuming $N_c p(s) \ll 1$, and $p(s) = P_{\text{gen}}(s)q^{-1}$ (where $q$ accounts for selection [21] and $P_{\text{gen}}(s)$ is the probability of sequence generation [14]), the score simplifies to Eq 1, with $\gamma = -\ln(qN_c\langle f^2 \rangle) = \ln(q^{-1}\langle f \rangle/\langle f^2 \rangle)$. The factor $\gamma$ depends on unknown parameters of the model, but can be estimated assuming a power-law for the clone size distribution [48], $\rho(f) \propto f^{-2}$ extending from $f = 10^{-11}$ to $f = 0.01$, and $q = 0.01$ [21], yielding $\gamma \approx 12.24$. Alternatively we optimized $\gamma$ to minimize the AUROC, yielding $\gamma \approx 15$ (S9 Fig). Since performance degrades quickly for larger values, we conservatively set $\gamma = 12$.

### 4.3 Estimating mean scores from RepSeq datasets

The sampling of $M$ cells from blood is simulated using large repertoire datasets. In a bulk repertoire sequencing dataset, the absolute number of cells for each clonotype (cells with a specific receptor) is unknown, but the fraction of each clonotype can be estimated using the proportion of reads that are associated with this specific receptor. To estimate the autologous $S$ and $\mathcal{I}$ of two samples of size $M_1$ and $M_2$ in the absence of true replicates, we computed their expected values from a single dataset containing $N$ reads, from which two random subsamples of sizes $M_1$ and $M_2$ were taken. The mean value of $S$ is equal to $\langle S \rangle = \sum_s (1 - (1 - f(s))^{M_1})(1 - (1 - f(s))^{M_2})$, where $f(s)$ is the true (and unknown) frequency of sequence $s$. A naive estimate of $\langle mS \rangle$ may be obtained by repeatedly resampling subsets of sizes $M_1$ and $M_2$ from the observed repertoire, calculate $S$ for each draw, and average. One gets the same result by replacing $f(s)$ by $\hat{f}_s = n(s)/N$ in the previous formula, where $n(s)$ is the number of $s$ reads in the full dataset, and $N = \sum_s n(s)$. However, this naive estimate leads to a systematic overestimate of the sharing (visible when compared with biological replicates, see S10 Fig), simply because this procedure overestimates the probability of resampling rare sequences, in particular singletons whose true frequency may be much lower than $1/N$. A similar bias occurs when computing $\mathcal{I}$. To correct for this bias, we look for a function $h(n)$ that satisfies for all $f$:

$$ \langle h(n) \rangle = \sum_n \binom{N}{n} f^n (1-f)^{N-n} h(n) = (1 - (1-f)^{M_1}) (1 - (1-f)^{M_2}), \qquad (9) $$

so that $\langle S \rangle$ and $\langle \mathcal{I} \rangle$ can be well approximated by:

$$ \langle S \rangle \approx \sum_s h(n(s)), \qquad (10) $$

$$ \langle \mathcal{I} \rangle \approx -\sum_s h(n(s)) [\ln(1/P_{\text{gen}}(s)) - \gamma]. \qquad (11) $$

Expanding the right-hand side of Eq 9 into 4 terms, we find that $h(n) = 1 - g_{M_1}(n) - g_{M_2}(n) + g_{M_1+M_2}(n)$ satisfies Eq 9 provided that:

$$ \sum_n \binom{N}{n} f^n (1-f)^{N-n} g_M(n) = (1-f)^M. \qquad (12) $$
Under the change of variable $x = f(1-f)$, the expression becomes:

$$ \sum_{n} \binom{N}{n} x^n g_M(n) = (1+x)^{N-M} = \sum_{n} \binom{N-M}{n} x^n. \quad (13) $$

Identifying the polynomial coefficients in $x^n$ on both sides yields:

$$ g_M(n) = \binom{N-M}{n} / \binom{N}{n}. \quad (14) $$

These corrected estimates agree with the direct estimates using biological replicates (S10 Fig).

Similarly, $\langle S \rangle$ and $\langle I \rangle$ in heterologous samples can be estimated using:

$$ \langle S \rangle \approx \sum_s [1 - g_{M_1}(n(s))] [1 - g_{M_2}(n'(s))], \quad (15) $$

$$ \langle I \rangle \approx \sum_s [1 - g_{M_1}(n(s))] [1 - g_{M_2}(n'(s))] [\ln(1/P_{\text{gen}}(s)) - \gamma]. \quad (16) $$

where $n(s)$ and $n'(s)$ are the empirical counts of sequence $s$ in the two samples.

## 4.4 Theoretical upper bound on heterologous scores

When the two samples were extracted from two different individuals (heterologous scenario), we can use the universality of the recombination process to give upper bounds on the values of $S$ and $I$. These bounds are represented by the dashed lines in Fig 1C). If two samples of respectively $M_1$ and $M_2$ cells, containing $T_1 \le M_1$ and $T_2 \le M_2$ unique sequences are extracted from two different individuals, the number of shared sequences between them is given by [21]:

$$ \langle S \rangle_{\text{heterologous}} \le \sum_s (1 - (1 - p(s))^{T_1}) (1 - (1 - p(s))^{T_2}) \quad (17) $$

$$ \lesssim T_1 T_2 \sum_s p(s)^2 = T_1 T_2 \langle p(s) \rangle \le M_1 M_2 \langle p(s) \rangle. \quad (18) $$

$p(s)$ is the probability of finding a sequence $s$ in the blood. Following [21], we make the approximation $p(s) = P_{\text{gen}}(s)q^{-1}$, where the $q = 0.01$ factor is the probability that a generated sequence passes selection. Then $\langle p(s) \rangle$ can be estimated from the mean over generated sequences. Similarly, we can estimate $I$ as

$$ \langle I \rangle_{\text{heterologous}} \lesssim T_1 T_2 \sum_s p(s)^2 [\ln(1/P_{\text{gen}}(s)) - \gamma] \quad (19) $$

$$ = -T_1 T_2 \langle p(s)[\gamma + \ln(qp(s))] \rangle \le -M_1 M_2 \langle p(s)[\gamma + \ln(qp(s))] \rangle, \quad (20) $$

which is also estimated from the mean over generated sequences.

## 4.5 Error rate estimates

To make the quantitative predictions shown in Fig 1, we need to constrain the tail behavior of the distributions of $S$ and $I$, for the two scenarios.

The $S$ statistic can be rewritten as a sum of Bernoulli variables over all possible sequences, each with a parameter corresponding to its probability of being present in both samples, either
in the autologous or the heterologous case. Therefore S is a Poisson binomial distribution, a sum of independent Bernoulli variables with potentially different parameters. The variance and tails of that distribution are bounded by those of the Poisson distribution with the same mean, denoted by $m_a$ for the autologous case, and $m_h$ for the heterologous case (S11 Fig).

Thanks to that inequality, the rates of false negatives and false positives for a given threshold r are bounded by:

$$P(S < r|\text{autologous}) \le Q(r + 1, m_a), \quad P(S > r|\text{heterologous}) \le 1 - Q(r + 1, m_h), \quad (21)$$

where Q is the regularized gamma function, which appears in the cumulative distribution function of the Poisson distribution. The mean autologous score $m_a$ is estimated from experimental data: we use the smallest value of $\langle S \rangle$ in the Emerson dataset and Eq 10. To compute $m_h$, we use the semi-theoretical prediction made using the universality of the recombination process, Eq 17.

Similarly, $\mathcal{I}$ can be viewed as a sum of S independent random variables, all following the distribution of $\ln(1/P_{\text{gen}})-\gamma$. However, this distribution differs in the two scenarios. Sequences shared between more than one donor have an higher probability of being generated, their $\ln(P_{\text{gen}})$ distribution has higher mean and smaller variance (S12 Fig).

The sum is composed of a relatively large number of variables in most realistic scenarios. Hence, we rely on the central limit theorem to approximate it by a normal distribution, of mean and variance proportional to S. Explicitly:

$$P(\mathcal{I} < r|\text{autologous}) = \frac{1}{2} \sum_{S=0}^{\infty} \frac{(m_a)^S e^{-m_a}}{S!} \left( 1 + \text{erf} \left( \frac{r - S \langle \ln[1/P_{\text{gen}}) - \gamma \rangle}{\sqrt{2S\text{Var}(\langle \ln(1/P_{\text{gen}}) - \gamma \rangle)}} \right) \right) \quad (22)$$

$$P(\mathcal{I} > r|\text{heterologous}) = \frac{1}{2} \sum_{S=0}^{\infty} \frac{(m_h)^S e^{-m_h}}{S!} \left( 1 - \text{erf} \left( \frac{r - S \langle \ln(1/P_{\text{gen}}) - \gamma \rangle_{\text{shared}}}{\sqrt{2S\text{Var}(\langle \ln(1/P_{\text{gen}}) - \gamma \rangle_{\text{shared}})}} \right) \right). \quad (23)$$

The AUROC are computed based on these estimates, by numerically integrating the true positive rate $P(S, \mathcal{I} < r|\text{heterologous})$ with respect to the false negative rate $P(S, \mathcal{I} < r|\text{autologous})$ as the threshold r is varied.

## 4.6 Modeling the evolution of autologous scores

We use the model of Ref. [38] to describe the dynamics of individual T-cell clone frequencies f under a fluctuating growth rate reflecting the changing state of the environment and the random nature of immune stimuli:

$$\frac{df}{dt} = \left[ -\frac{1}{\tau} + \frac{1}{2\theta} + \frac{1}{\sqrt{\theta}} \eta(t) \right] f(t), \quad (24)$$

where $\eta(t)$ is a Gaussian white noise with $\langle \eta(t) \rangle = 0$ and $\langle \eta(t) \eta(t') \rangle = \delta(t-t')$.

With the change of variable $x = \ln(f)$, these dynamics simplify to a simple Brownian motion in log-frequency: $\partial_t x = -\tau^{-1} + \theta^{-1/2} \eta(t)$. In that equation, $\tau$ appears as the decay rate of the frequency, while $\theta$ is the timescale of the noise, interpreted as the typical time it takes for the frequency to rise or fall by a logarithmic unit owing to fluctuations. Considering a large population of clones, each with their independent frequency evolving according to Eq 24, and a source term at small f corresponding to thymic exports, one can show that the steady-state probability density function of f follows a power-law [38], $\rho(f) \propto f^{-\alpha}$, with exponent
$α = 1 + 2θ/τ$. $α$ was empirically found to be ≈2 in a wide variety of immune repertoires [10, 48–50], implying $2θ ≈ τ$. The turn-over time $τ$ is unknown, and was varied from 0.5 years to 10 years in the simulations.

We simulated the evolution of human TRB repertoires by starting with the empirical values of the frequencies of each observed clones, $f(s, 0) = \hat{f}(s, 0) = n(s, 0)/N$ from the analysed datasets. A sample of size $M$ was randomly selected with respect to these frequencies, and the frequencies of the clones captured in that sample were then evolved with a time-step of 2 days using Euler-Maruyama's method, which is exact in the case of Brownian motion. Clones with frequencies falling below $10^{-11}$ (corresponding to a single cell in the organism) were removed. At each time $t > 0$, we measured the mean value of $S$ with the formula $\sum_s (1 - (1 - f(s, t))^M)$ where the sum runs over the sequences captured in the initial sample.

## Supporting information

S1 Fig. Mean value of log $P_{gen}$ as a function of the rank of the clone type (from most abundant to least abundant) from [17]. The black line represents the mean of log $P_{gen}$ for naive clones. The statistic for the top-clones (low rank) is similar to the one for the naive clones.
(TIF)

S2 Fig. Comparison between the mean of $S$ (autologous case), and three common diversity measures. The number of unique sequences found in the dataset (top left), the Shannon index, $-\sum \hat{f}_s \ln \hat{f}_s$ (top right), the Simpson index (bottom left), and the total number of reads in each datasets (bottom right). All the diversity measures show a strong correlation with $S$, but the correlation with the sequencing depth is low.
(TIF)

S3 Fig. Detection Error Trade-off (DET) graph for both summary statistics, between a large sample (full dataset, $M_1 = 10^6$) and a smaller one, of size $M_2 = M$.
(TIF)

S4 Fig. Mean value of $S$ (left) and $\mathcal{I}$ (right) for the IGH chain (part of the B-cell receptor) for sample sizes $M = 5000$. The IGH sequences used are restricted to IgG+ B-cells (selected according to their CH gene). The sequences were obtained from 8 different individuals (6 biological replicates each) in the dataset from [30]. Autologous (blue) and heterologous (yellow) are well separated.
(TIF)

S5 Fig. Longitudinal evolution of $S$. A) Evolution of $S$ ($M = 5000$) during vaccination, between a sample taken at day 0 (vaccination date) and at a later timepoint. Each color represents a different individual. Each pair timepoint/individual has two biological replicates. The dashed line represents the threshold value. B) Evolution of $S$ between a sample taken at year 0 and a later timepoint. The red histogram corresponds to one of the individuals sampled in [18] and the blue curves show theoretical estimates, fitted to match ($τ = 0.66$). C) Evolution of the (normalized) mean of $S$ ($M = 5000$) as a function of time for different values of the turnover rate $τ$. The dashed line represents the threshold value divided by the smallest value of $S_{t=0}$ ($M = 5000$) in the data. The data points were obtained from the datasets [35] (yellow), [18] (green) and [22] (orange). Different markers indicate different individuals.
(TIF)

S6 Fig. For each pair of different individuals from the dataset [17], the values of $S$ (blue) and $\mathcal{I}$ (yellow) are shown as a function of the serological status of the two individuals—CMV+/CMV+ (resp. CMV-/CMV-) if both are positive (negative) to cytomegalovirus and
CMV+/CMV- if one is positive while the other one is not. Co-infection with cytomegalovirus does not hamper the capacity of Immprint to distinguish between different people.
(TIF)

S7 Fig. Distribution of S (left) and I (right) for CD4+ and CD8+ T-cells in the dataset from [18] (6 individuals). In the autologous case (blue), there is a significant difference between the two cell types, caused by the differences between their clonal distributions.
(TIF)

S8 Fig. Left: Mean value of $P_{gen}$ as a function of the rank of the clonal type, for generic sequences (blue) and sequences shared between more than two donors (orange). The mean stays flat indicating that the probability of being generated does not generally depend on the clonal type size. There is an exception (black rectangle), shown as a close-up on the right panel. The top twenty clones, when shared between donors, have a smaller probability of being generated than expected by chance. This difference is likely to be driven by convergent selection against common pathogens, since CMV positive donors show a more pronounced effect than CMV negative ones.
(TIF)

S9 Fig. Left panel: AUROC (Area Under Receiver Operating Characteristic) of $\mathcal{I}$, as a function of $\gamma$ ($M = M_1 = M_2 = 5000$). We observe an optimum near $\gamma = 15$. Right panel: AUROC as a function of $M$, for $S$, $\mathcal{I}(\gamma = 0)$, and $\mathcal{I}(\gamma = 15)$.
(TIF)

S10 Fig. Naive and corrected estimates of the autologous $S$ from single datasets, versus its values computed using true biological replicates from Ref. [18].
(TIF)

S11 Fig. Comparison between the distribution of $S$ obtained by computationally and repeatedly downsampling a single repertoire from Ref. [17] with $M = 5,000$ (histogram), and a Poisson distribution of the same mean (full line).
(TIF)

S12 Fig. Distribution of $\ln(P_{gen})$ for generic sequences, and for sequences shared between heterologous samples.
(TIF)

**Acknowledgments**

The authors are grateful for the help of Natanael Spisak with the analysis of BCR repertoire datasets.

**Author Contributions**

**Conceptualization:** Thierry Mora, Aleksandra M. Walczak.

**Data curation:** Thomas Dupic, Meriem Bensouda Koraichi.

**Formal analysis:** Thomas Dupic, Meriem Bensouda Koraichi, Thierry Mora, Aleksandra M. Walczak.

**Funding acquisition:** Anastasia A. Minervina, Mikhail V. Pogorelyy, Thierry Mora, Aleksandra M. Walczak.
**Investigation:** Thomas Dupic, Meriem Bensouda Koraichi, Anastasia A. Minervina, Mikhail V. Pogorelyy, Thierry Mora, Aleksandra M. Walczak.

**Methodology:** Thomas Dupic, Meriem Bensouda Koraichi, Thierry Mora, Aleksandra M. Walczak.

**Project administration:** Thierry Mora, Aleksandra M. Walczak.

**Software:** Thomas Dupic, Meriem Bensouda Koraichi.

**Supervision:** Thierry Mora, Aleksandra M. Walczak.

**Visualization:** Thomas Dupic, Meriem Bensouda Koraichi, Anastasia A. Minervina, Mikhail V. Pogorelyy, Thierry Mora, Aleksandra M. Walczak.

**Writing – original draft:** Thomas Dupic, Thierry Mora, Aleksandra M. Walczak.

**Writing – review & editing:** Thomas Dupic, Meriem Bensouda Koraichi, Anastasia A. Minervina, Mikhail V. Pogorelyy, Thierry Mora, Aleksandra M. Walczak.

## References

1. Homer N, Szelinger S, Redman M, Duggan D, Tembe W, Muehling J, et al. Resolving Individuals Contributing Trace Amounts of DNA to Highly Complex Mixtures Using High-Density SNP Genotyping Microarrays. PLOS Genetics. 2008; 4(8):e1000167. https://doi.org/10.1371/journal.pgen.1000167 PMID: 18769715

2. Naveed M, Ayday E, Clayton EW, Fellay J, Gunter CA, Hubaux JP, et al. Privacy in the Genomic Era. ACM computing surveys. 2015; 48(1). https://doi.org/10.1145/2767007 PMID: 26640318

3. Sweeney L, Abu A, Winn J. Identifying Participants in the Personal Genome Project by Name. Rochester, NY: Social Science Research Network; 2013. ID 2257732.

4. Hozumi N, Tonegawa S. Evidence for somatic rearrangement of immunoglobulin genes coding for variable and constant regions. Proc Natl Acad Sci USA. 1976; 73(10):3628-3632. https://doi.org/10.1073/pnas.73.10.3628

5. Robins H. Immunosequencing: applications of immune repertoire deep sequencing. Current Opinion in Immunology. 2013; 25(5):646-52. https://doi.org/10.1016/j.coi.2013.09.017

6. Attaf M, Huseby E, Sewell AK. αβ T cell receptors as predictors of health and disease. Cellular & molecular immunology. 2015; 12(4):391-399. https://doi.org/10.1038/cmi.2014.134

7. Woodsworth DJ, Castellarin M, Holt Ra. Sequence analysis of T-cell repertoires in health and disease. Genome medicine. 2013; 5(10):98. https://doi.org/10.1186/gm502

8. Bradley P, Thomas PG. Using T Cell Receptor Repertoires to Understand the Principles of Adaptive Immune Recognition. Annual Review of Immunology. 2019; 37(1):547-570. https://doi.org/10.1146/annurev-immunol-042718-041757

9. Davis MM, Boyd SD. Recent progress in the analysis of αβ T cell and B cell receptor repertoires. Current Opinion in Immunology. 2019; 59:109-114. https://doi.org/10.1016/j.coi.2019.05.012

10. Mora T, Walczak A. Quantifying Lymphocyte Receptor Diversity. é. 2016;

11. Qi Q, Liu Y, Cheng Y, Glanville J, Zhang D, Lee JY, et al. Diversity and clonal selection in the human T-cell repertoire. Proceedings of the National Academy of Sciences of the United States of America. 2014; 111(36):13139-44. https://doi.org/10.1073/pnas.1409155111 PMID: 25157137

12. Lythe G, Callard RE, Hoare RL, Molina-París C. How many TCR clonotypes does a body maintain? Journal of Theoretical Biology. 2016; 389:214-224.

13. Murugan A, Mora T, Walczak AM, Callan CG. Statistical inference of the generation probability of T-cell receptors from sequence repertoires. Proceedings of the National Academy of Sciences. 2012; 109 (40):16161-16166. https://doi.org/10.1073/pnas.1212755109

14. Marcou Q, Mora T, Walczak AM. High-Throughput Immune Repertoire Analysis with IGoR. Nature Communications. 2018; 9(1):561. https://doi.org/10.1038/s41467-018-02832-w PMID: 29422654

15. Sethna Z, Isacchini G, Dupic T, Mora T, Walczak AM, Elhanati Y. Population Variability in the Generation and Thymic Selection of T-Cell Repertoires. bioRxiv. 2020; p. 2020.01.08.899682.

16. Venturi V, Price DA, Douek DC, Davenport MP. The molecular basis for public T-cell responses? Nature Reviews Immunology. 2008; 8(3):231-238.
17. Emerson RO, DeWitt WS, Vignali M, Gravley J, Hu JK, Osborne EJ, et al. Immunosequencing identifies signatures of cytomegalovirus exposure history and HLA-mediated effects on the T cell repertoire. Nature Genetics. 2017; 49(5):659–665 https://doi.org/10.1038/ng.3822 PMID: 28369038

18. Pogorelyy MV, Minervina AA, Touzel MP, Sycheva AL, Komech EA, Kovalenko EI, et al. Precise Tracking of Vaccine-Responding T Cell Clones Reveals Convergent and Personalized Response in Identical Twins. Proceedings of the National Academy of Sciences. 2018; p. 201809642. https://doi.org/10.1073/pnas.1809642115 PMID: 30459272

19. Zvyagin IV, Pogorelyy MV, Ivanova ME, Komech EA, Shugay M, Bolotin DA, et al. Distinctive Properties of Identical Twins' TCR Repertoires Revealed by High-Throughput Sequencing. Proceedings of the National Academy of Sciences of the United States of America. 2014; 111(16):5980–5985. https://doi.org/10.1073/pnas.1319389111 PMID: 24711416

20. Tanno H, Gould TM, McDaniel JR, Cao W, Tanno Y, Durrett RE, et al. Determinants Governing T Cell Receptor α/β-Chain Pairing in Repertoire Formation of Identical Twins. Proceedings of the National Academy of Sciences. 2020; 117(1):532–540. https://doi.org/10.1073/pnas.1915008117 PMID: 31879353

21. Elhanati Y, Sethna Z, Callan CG, Mora T, Walczak AM. Predicting the spectrum of TCR repertoire sharing with a data-driven model of recombination. Immunological Reviews. 2018; 284(1):167–179. https://doi.org/10.1111/imr.12665

22. Britanova OV, Shugay M, Merzlyak EM, Staroverov DB, Putintseva EV, Turchaninova MA, et al. Dynamics of Individual T Cell Repertoires: From Cord Blood to Centenarians. The Journal of Immunology. 2016; 196(12):5005–5013. https://doi.org/10.4049/jimmunol.1600005 PMID: 27183615

23. Sylwester AW, Mitchell BL, Edgar JB, Taormina C, Pelte C, Ruchti F, et al. Broadly targeted human cytomegalovirus-specific CD4+ and CD8+ T cells dominate the memory compartments of exposed subjects. Journal of Experimental Medicine. 2005; 202(5):673–685. https://doi.org/10.1084/jem.20050882 PMID: 16147978

24. Khan N, Shariff N, Cobbold M, Bruton R, Ainsworth JA, Sinclair AJ, et al. Cytomegalovirus Seropositivity Drives the CD8 T Cell Repertoire Toward Greater Clonality in Healthy Elderly Individuals. The Journal of Immunology. 2002; 169(4):1984–1992. https://doi.org/10.4049/jimmunol.169.4.1984 PMID: 12165524

25. Dash P, McClaren JL, Oguin TH, Rothwell W, Todd B, Morris MY, et al. Paired analysis of TCRα and TCRβ chains at the single-cell level in mice. Journal of Clinical Investigation. 2011; 121(1):288–295. https://doi.org/10.1172/JCI44752 PMID: 21135507

26. Redmond D, Poran A, Elemento O. Single-Cell TCRseq: Paired Recovery of Entire T-Cell Alpha and Beta Chain Transcripts in T-Cell Receptors from Single-Cell RNAseq. Genome Medicine. 2016; 8(1):80. https://doi.org/10.1186/s13073-016-0335-7

27. Grigaityte K, Carter JA, Goldfless SJ, Jeffery EW, Ronald J, Jiang Y, et al. Single-cell sequencing reveals αβ chain pairing shapes the T cell repertoire. bioRxiv:213462. 2017;

28. Howie B, Sherwood AM, Berkebile AD, Berka J, Emerson RO, Williamson DW, et al. High-throughput pairing of T cell receptor a and b sequences. Science translational medicine. 2015; 7(301):301ra131. https://doi.org/10.1126/scitranslmed.aac5624 PMID: 26290413

29. Dupic T, Marcou Q, Walczak AM, Mora T. Genesis of the Aβ T-Cell Receptor. PLOS Computational Biology. 2019; 15(3):e1006874. https://doi.org/10.1371/journal.pcbi.1006874

30. Briney B, Inderbitzin A, Joyce C, Burton DR. Commonality despite Exceptional Diversity in the Baseline Human Antibody Repertoire. Nature. 2019; 566(7744):393–397. https://doi.org/10.1038/s41586-019-0879-y PMID: 30664748

31. DeWitt WS, Emerson RO, Lindau P, Vignali M, Snyder TM, Desmarais C, et al. Dynamics of the Cytotoxic T Cell Response to a Model of Acute Viral Infection. Journal of Virology. 2015; 249(February): JVI.03474–14. https://doi.org/10.1128/JVI.03474-14 PMID: 25653453

32. Wolf K, Hether T, Gilchuk P, Kumar A, Rajeh A, Schiebout C, et al. Identifying and Tracking Low-Frequency Virus-Specific TCR Clonotypes Using High-Throughput Sequencing. Cell Reports. 2018; 25(9):2369–2378.e4. https://doi.org/10.1016/j.celrep.2018.11.009 PMID: 30485806

33. Qi Q, Cavanagh MM, Le Saux S, NamKoong H, Kim C, Turgano E, et al. Diversification of the antigen-specific T cell receptor repertoire after varicella zoster vaccination. Science Translational Medicine. 2016; 8(332):332ra46–332ra46. https://doi.org/10.1126/scitranslmed.aaf1725 PMID: 27030598

34. Sycheva AL, Pogorelyy MV, Komech EA, Minervina AA, Zvyagin IV, Staroverov DB, et al. Quantitative profiling reveals minor changes of T cell receptor repertoire in response to subunit inactivated influenza vaccine. Vaccine. 2018; 36(12):1599–1605. https://doi.org/10.1016/j.vaccine.2018.02.027 PMID: 29454515
35. Chu ND, Bi HS, Emerson RO, Sherwood AM, Birnbaum ME, Robins HS, et al. Longitudinal Immunose-quencing in Healthy People Reveals Persistent T Cell Receptors Rich in Highly Public Receptors. BMC Immunology. 2019; 20(1):19. https://doi.org/10.1186/s12865-019-0300-5 PMID: 31226930

36. Borghans JAM, De Boer RJ. Quantification of T-cell dynamics: From telomeres to DNA labeling. Immunological Reviews. 2007; 216(1):35-47. https://doi.org/10.1111/j.1600-065X.2007.00497.x PMID: 17367333

37. Thomas-Vaslin V, Altes HK, de Boer RJ, Klatzmann D. Comprehensive Assessment and Mathematical Modeling of T Cell Population Dynamics and Homeostasis. The Journal of Immunology. 2008; 180(4):2240-2250. https://doi.org/10.4049/jimmunol.180.4.2240

38. Desponds J, Mora T, Walczak AM. Fluctuating Fitness Shapes the Clone-Size Distribution of Immune Repertoires. Proceedings of the National Academy of Sciences. 2016; 113(2):274-279. https://doi.org/10.1073/pnas.1512977112

39. de Greef PC, Oakes T, Gerritsen B, Ismail M, Heather JM, Hermsen R, et al. The naive t-cell receptor repertoire has an extremely broad distribution of clone sizes. eLife. 2020; 9:1-24. https://doi.org/10.7554/eLife.49900 PMID: 32187010

40. De Boer RJ, Perelson AS. Quantifying T Lymphocyte Turnover. Journal of theoretical biology. 2013; 327:45-87. https://doi.org/10.1016/j.jtbi.2012.12.025

41. Buhler S, Bettens F, Dantin C, Ferrari-Lacraz S, Ansari M, Mamez AC, et al. Genetic T-cell receptor diversity at 1 year following allogeneic hematopoietic stem cell transplantation. Leukemia. 2020; 34(5):1422-1432. https://doi.org/10.1038/s41375-019-0654-y PMID: 31772297

42. Pogorelyy MV, et al. (2018) Precise tracking of vaccine-responding T cell clones reveals convergent and personalized response in identical twins. Proceedings of the National Academy of Sciences 115:12704-12709. https://doi.org/10.1073/pnas.1809642115 PMID: 30459272

43. Li B, Li T, Wang B, Dou R, Zhang J, Liu JS, et al. Ultrasensitive detection of TCR hypervariable-region sequences in solid-tissue RNA-seq data. Nature Genetics. 2017; 49(4):482-483. https://doi.org/10.1038/ng.3820 PMID: 28358132

44. Bolotin DA, Poslavsky S, Davydov AN, Frenkel FE, Fanchi L, Zolotareva OI, et al. Antigen receptor repertoire profiling from RNA-seq data. Nature Biotechnology. 2017; 35(10):908-911. https://doi.org/10.1038/nbt.3979 PMID: 29020005

45. Shugay M, Britanova OV, Merzlyak EM, Turchaninova MA, Mamedov IZ, Tuganbaev TR, et al. Towards Error-Free Profiling of Immune Repertoires. Nature Methods. 2014; 11(6):653-655. https://doi.org/10.1038/nmeth.2960 PMID: 24793455

46. Bolotin DA, Poslavsky S, Mitrophanov I, Shugay M, Mamedov IZ, Putintseva EV, et al. MiXCR: Software for Comprehensive Adaptive Immunity Profiling. Nature Methods. 2015; 12(5):380-381. https://doi.org/10.1038/nmeth.3364 PMID: 25924071

47. Sethna Z, Elhanati Y, Callan C, Walczak AM, Mora T. OLGA: Fast Computation of Generation Probabilities of B- and T-Cell Receptor Amino Acid Sequences and Motifs. Bioinformatics. 2019; 35(17):2974-2981. https://doi.org/10.1093/bioinformatics/btz035

48. Touzel MP, Walczak AM, Mora T. Inferring the immune response from repertoire sequencing. PLoS Computational Biology. 2020; 16(4):1-21.

49. Weinstein JA, Jiang N, White RA, Fisher DS, Quake SR. High-throughput sequencing of the zebrafish antibody repertoire. Science. 2009; 324(5928):807-810. https://doi.org/10.1126/science.1170020

50. Oakes T, Heather JM, Best K, Byng-Maddick R, Husovsky C, Ismail M, et al. Quantitative characterization of the T cell receptor repertoire of naïve and memory subsets using an integrated experimental and computational pipeline which is robust, economical, and versatile. Frontiers in Immunology. 2017; 8(OCT):1-17.