ABSTRACT Immunotherapy is now an essential tool for cancer treatment, and the unique features of an individual’s T cell receptor repertoire are known to play a key role in its effectiveness.

The repertoire, famously vast due to a cascade of cellular mechanisms, can be quantified using repertoire sequencing. In this study, we sampled the repertoire over several time points

following treatment with anti-CTLA-4, in a syngeniec mouse model for colorectal cancer, generating a longitudinal dataset of T cell clones and their abundance. The dynamics of the repertoire

can be observed in response to treatment over a period of four weeks, as clonal expansion of specific clones ascends and descends. The data made available here can be used to determine

October 2022 DYNAMICS AND SPECIFICITIES OF T CELLS IN CANCER IMMUNOTHERAPY Article 12 April 2023 BACKGROUND & SUMMARY The immune repertoire is a collection of trans-membrane proteins

that assist with the recognition of different antigens and initiate an immune response. These proteins, located on the surface of T cells and known as T cell receptors (TCRs), are made of

several chains (alpha through delta) and while sharing a general molecular structure, each TCR is slightly different. This repertoire has the potential to generate more than 1015 1 different

T cell receptors in humans. Once activated, the immune system sees major changes in its T cell composition; some T cells undergo clonal expansion2 and proliferation3, which later subsides.

These dynamics govern the outcome of the immune response, and are key to understanding both the disease itself and its effective treatment4. Fig. 1 is a schematic view of the immune response

in the face of infection, which also highlights the temporal nature of the response, rather than a binary state. The immune system must face a varied list of potential threats; be it

viruses, bacteria, or the body’s own malfunctioning cells. Recruitment and activation of immune cells is needed. For this to happen, the immune system must be diverse enough to interact with

 different antigens. This is achieved with three cellular and molecular mechanisms: 1) the recombination of the T cell alpha genes and beta genes by the RAG1–2 protein complex during T cell

development, 2) the stochastic nature of the formation of the third complementarity-determining region (CDR3) of the T cell, and 3) the combination of alpha and beta chains. This entire

process is termed V(D)J recombination5,6 (Fig. 2) and is key to the development of the adaptive immune system, and the correct formation of the rearranged receptors is critical to their 

association with  antigens. Immune checkpoint inhibitors (“ICI”)7,8,9 are used as treatment for several types of cancer. This kind of treatment will suppress natural checkpoint and

regulators of the immune system, resulting in an amplified immune response, one that is more likely to deal with cancerous cells. This type of treatment has recently been employed in

colorectal cancer (“CRC”), with several positive results reported, specifically in metastatic CRC that is mismatch-repair-deficient and microsatellite instability-high (dMMR–MSI-H)10. The

fundamental involvement of the immune response in this type of cancer makes ICI and its effect on the T cell repertoire especially important. CTLA4 is a negative regulator of the immune

response11,12, aimed at preventing hyperactivation and adjusting the intensity of the immune response by causing clonal anergy of lymphocytes. CTLA4 operates on several cellular mechanisms,

such as CD28 antagonism in binding to B7 ligands and other co-stimulatory ligands, recruiting T cell inhibitory effectors, and preventing immune stimulating antibody conjugate3. Anti-CTLA4

will block all of the above, and mediate a stronger immune attack13. Unlike the genetic code, the TCR repertoire is dynamic; its state, made readable with RNA sequencing, indicates the

current status of the immune system. As such, several methods have been developed to quantify the TCR repertoire and its diverse components14,15,16, and many computational and analytical

tools are constantly being developed and used to analyze these repertoires17,18,19. Many works have highlighted the importance of a diverse TCR repertoire for the immune system to perform

optimally1,20, and while these works track and document changes to the repertoire, the supporting data is limited to pre- and post- treatment time points at best. Additionally, since a

phenotypical change, such as reduction in tumor size, can occur some time after the peak of the immunological response, it can be difficult to create a snapshot of the immune state at the

most relevant time. Moreover, the high variability in the favorable reaction of cancers to immunotherapy, makes it even challenging for researcher to quantify the true dynamics of the immune

response21. In light of the above, we designed and carried out an experiment aimed at creating a temporal profile of the dynamics of the immune response as it occurs. These data, made

available here, would allow researchers from all backgrounds to observe the changes over time as they develop and eventually subside in response to treatment with anti-CTLA-4. Analyzing the

behavior of the immune response, presented here in several forms, from raw reads to clonotype listings, would benefit efforts to better understand the immune response in cancer, and

hopefully lead to better treatments. METHODS EXPERIMENTAL PROCEDURE In this experiment we used the mouse tumor cell line MC38, a grade III mouse colon adenocarcinoma, induce in a syngeneic

mouse model. An estimated 5·105 MC38 cells were subcutaneously implanted to immunocompetent C57BL/6 N mice of the C57Bl/6NAnNCrl strain. The untreated group (5 mice) was dosed with 10 

ml/kg/day of vehicle on day 1, 3 and 6 of the experiment. The treated group (15 mice) was dosed with the immune checkpoint inhibitor antibody anti-mCTLA4 (5/2.5 mg/kg/day at 1/3, 6).

Following treatment, blood samples for the isolation of PBMCs were taken at days 0, 7, 14, and 21. The experimental timeline is shown in Fig. 3. ANIMAL HEALTH AND ETHICS All experiments and

protocols were approved by the animal welfare body at CR Discovery Research Services Germany and the local authorities, and were conducted according to all applicable international, national

and local laws and guidelines. Only animals with unobjectionable health were selected to enter testing procedures. The animals were delivered at the age of four to six weeks and were used

for experiments after at least one week of acclimatization. Animals were arbitrarily numbered during tumor implantation using radio frequency identification transponders. Each cage was

labeled with a record card indicating all relevant experimental details. ANIMAL HOUSING CONDITIONS Animals were housed in individually ventilated type II and III cages, based on group size,

and were kept under a 14 L:10D artificial light cycle. The temperature inside the cages was maintained at 25 °C with a relative humidity of 40–70% and 60–65 air changes/hour in the cage.

Dust-free bedding consisting of aspen wood chips with approximate dimensions of 5 × 5 × 1 mm and additional nesting material was used. The cages, bedding, and the nesting material were

changed weekly. The animals were fed autoclaved Teklad Global 19% Protein Extruded Diet (T.2019S.12) from Envigo RMS SARL and had access to sterile filtered and acidified (pH 2.5) tap water

that was changed per week. Feed and water were provided ad libitum. All materials were autoclaved prior to use. IMMUNOTHERAPY TREATMENT A dosing solution with a concentration of 0.5 mg/ml

for dosing at 5 mg/kg/day was freshly prepared on dosing day 1 by diluting 203 _μ_L of anti-mCTLA4 stock solution with 2797 _μ_L of PBS. A dosing solution with a concentration of 0.25 mg/ml

for dosing at 2.5 mg/kg/day was freshly prepared on dosing day 3 and 6 by diluting 101 _μ_L of anti-mCTLA4 with 2899 _μ_L of PBS. All dosing solutions were administered in a dose volume of

10 ml/kg. PBMCS PURIFICATION FROM BLOOD SAMPLES 100_μ_L of blood was collected by retro-orbital sinus puncture under isoflurane anesthesia. PBMCs were prepared from 100 _μ_L blood by

spinning down cells and lysing pellets with ACK lysis buffer (150 mM ammonium chloride, 10 mM potassium bicarbonate, 0.1 mM EDTA, pH 7.2–7.4) to remove erythrocytes. The resulting PBMC were

frozen in 10% DMSO/90% FCS at 80 °C until shipment. SOLID SAMPLE COLLECTION Tumors (where still visible) were collected immediately after euthanasia, and directly transferred to liquid

nitrogen (snap-frozen samples). Spleens were collected immediately after euthanasia, and directly transferred to liquid nitrogen (snap-frozen samples). RNA PURIFICATION FROM BLOOD, TUMOR AND

SPLEEN SAMPLES PBMCs samples were thawed, and their total RNA content extracted using the RNeasy mini kit (Qiagen). Tumor and spleen samples were thawed, weighted and applied to Trizol

(Invitrogen) purification procedure according to the manufacturer protocol. The obtained RNA product was then evaluated for its integrity value (RIN) to better standardize RNA quality. Next,

RNA concentrations were determined by qubit, which provides an accurate method for the quantitation of low abundance RNA samples. LIBRARY PREPARATION A fixed total RNA concentration of 250 

ng from each sample was subjected to the SMARTer TCR profiling kit (Takara); a single cDNA strand is synthesized at the first round of reverse transcription, using poly T primer for the 3’

end and a smarter oligo (Switching Mechanism at 5’ End of RNA Template) for the 5’ ends, catalyzed by MMLV-derived SMARTScibe Reverse Transcriptase. By employing 5’ RACE-like approach, the

SMARTer TCR profiling kit allows capturing complete V(D)J variable regions of TCR alpha and beta chains and in the following step, each fragment was exponentially increased by PCR. On each

dsDNA fragment, dedicated adapters containing flow cell-binding sites, P5 and P7, were inserted to allow the enriched library fragments to attach to the cell flow surface. The TCR sequencing

library was size-selected and purified using AMPure XP beads. The generated libraries were measured for their DNA concentration by qubit and assessed for their size by Tapestation or

Bioanalyzer. SEQUENCING AND READS PROCESSING TCR sequencing was performed on an Illumina Miseq sequencer using the 600-cycle Miseq reagent kit v3 (Illumina) with pair-end, 2 × 300 base pair

reads. Library pool molarity was 4 nM. Transforming raw reads to alpha and beta CDR3 sequences was done by using MiXCR, a universal framework that processes big immunome data from raw

sequences and output quantitated clonotypes (Bolotin _et al_.18). For the alpha chain we obtained 13·106 reads with an average of 171628 copies and 25159 unique amino-acid clones, and for

the beta chain we obtained 18.6·106 with an average of 245867 copies and 49725 unique amino-acid clones. In order to avoid bias originating from different sample sizes, we sub-sampled 97048

sequences from each sample, as the number of clones in the smallest sample. RANDOMIZATION AND BLINDING The initial division of treatment (15 mice, labeled F to T) and control groups (5 mice,

labeled A to E) was done randomly. At the end of the experiment, each mouse was assigned to a response group according to the changes to its tumor size in response to treatment: Early

responders (mice L, N, P, Q, R) – mice that showed early response to treatment and reduction in tumor volume; Mid-responders – mice that showed medium response to treatment and reduction in

tumor volume (mice J, M, O, S, T), and Late responders - mice that showed later response to treatment and reduction in tumor volume (mice F, G, H, I, K). CLONOTYPING To generate the alpha

and beta chains clonotyping lists, the tool MiXCR was used18. MiXCR is a universal framework that processes big immunome data from raw sequences to quantitated clonotypes. MiXCR was used

with the following arguments: −−_speciesmmu_−−_libraryde f ault_ −−_threads_12−−_reportreport_ − _prna_−_seq_−_OvParameters.geneFeatureToAlign_ = 

_VTranscriptWithP_−_OvParameters_._parameters_. _f loatingLe f tBound_ = true−_OjParameters.parameters. f loatingRightBound_ = _f alse_−_OcParameters.parameters_. _f loatingRightBound_ = 

_true_. Sample code for similar use is provided in the code availability section of this work. DATA RECORDS The TCR-Seq data were deposited in the NCBI Sequence Read Archive (SRA) under the

accession number PRJNA65834822. This data appears as zipped FASTQ files, as is standard for next generation sequencing results. This data represents the TCR sequencing of all data points:

blood samples taken on days 0, 7, 14, and 21 from treated and control mice, as well as TCR sequencing of tumors and spleens taken at the end of the experiment. The Alpha and Beta chain TCRs

were deposited on FigShare under the following https://doi.org/10.6084/m9.figshare.1928750623 (alpha) and https://doi.org/10.6084/m9.figshare.1928751224 (beta). These datasets appear as

tabular text files, and include information on each clonotype in the data. This includes a unique clone ID, clone count, the clone’s fraction out of the total number of clones, the actual

clonal sequence in base-pairs, and more. For complete details of the fields populating this data, we suggest reading the documentation at https://mixcr.readthedocs.io/en/master/export.html.

Note that this table is sorted according to the clone count by default, with the most prevalent clonotypes appearing on top. The metadata for each clonotype dataset were deposited on

FigShare under the following https://doi.org/10.6084/m9.figshare.1931226225 (alpha) and https://doi.org/10.6084/m9.figshare.1931226526 (beta). These data sets appear as tabular text files,

and include meta-information on each data point. This includes the sample name (matching the RNA-Seq data files), the sample ID, Time point, a mouse identifier, the group (either treated or

control), the tissue from which the sample was taken (Blood, Tumor or Spleen), and how quickly, if at all, a given mouse responded to treatment (as determined by tumor size measurements

throughout the experiment) The entire data collection process, as well as the data made available in this work, is presented as a work flow in Fig. 4. TECHNICAL VALIDATION TCR-SEQ VALIDATION

To validate the raw sequencing data, we used FASTQC27, A quality control tool for high throughput sequence data. Results showed that all sequences had high per sequence quality scores”.

Additionally, all samples had valid percentage of base calls at each position for which an N was called, and all samples were of the same length of 301 basepairs, as expected. It should be

noted that as these sequences are of T cell receptors only, and do not align a standard mouse genome, as is evident in the low alignment score of all samples, as well as the presence of

promoters in the data. CLONOTYPING VALIDATION To demonstrate that the clonotype data is valid, we used Immunarch28 to measure the number and length of the quantified clonotypes for each

chain. As can be seen in Fig. 5, there is no significant difference in number of clonotypes among between the treatment groups, nor is there a significant difference between the number of

clonotypes between the response groups. Similarly, as can be seen in Fig. 6, there is no significant difference in the length of the CDR3 section of the clonotypes between the treatment

groups, nor is there a significant difference in the length of the CDR3 section of the clonotypes between the response groups. The above emphasizes the validity of the data, and allows for

unbiased comparison between different treatment or response groups. USAGE NOTES Reuse of the data provided herein is possible either by re-construction of the clonotypes from the raw

sequencing files, or using our built clonotyping list and continue the analysis from there. Users should note that when constructing their own clonotyping data, different arguments can be

used to obtain fields of interest in the data. The raw repertoire sequencing data could be used in several ways; first, clonotype quantification can be performed in several ways and under

different conditions29. We recommend the use of MiXCR (see code examples section for our application), but other tools such as IgBLAST30 and ImRep (https://github.com/mandricigor/imrep) have

useful features. Alternatively, the data can be used like any other short-read RNA sequencing data, although generated from T cells only; such a work-flow would include alignment, assembly

and quantification of gene expression. While tools for each step exist, we recommend using tools able to perform all of these steps, as exemplified in a recent review by Stark _et al_.31.

Once clonotyping is performed, these data, both of the alpha and beta chains, could be used with many existing, specialized pipelines. These include Immunarch28, Immunoseq analyzer32,

Immcantation Framework33 and VDJtools17. A recent useful review by Arunkumar _et al_. compares these tools and others19. Analysis such as this can be done to estimate the number of copies a

specific clonotype has, compare clonotypes to those found in external databases, or attempt to tie the presence of clonotypes with response to treatment, or lack thereof. Finally, the data

presented here could be used with machine learning approaches to predict response to treatment. This could be done either with dedicated platforms such as ImmuneML34, or with any other tool

or approach that could uncover a connection between response to treatment and the composition of the repertoire. Any such methods would benefit from the temporal nature of the data, and

could support approaches dedicated to precise treatment of individual patients. We would recommend the use of time-series analysis to the data presented here. CODE AVAILABILITY Several

scripts are made available to allow users of all backgrounds and experience levels to access this data. After downloading the raw sequencing data, users can use the following code35 to

perform clonotyping and obtain a list of alpha and beta chains from each sample, as well as a report file with details on the success of the process. Note that the user should have a recent

version of MiXCR installed, as detailed in the18 MiXCR documentation. When using the clonotype data as input, either after following the step above or by downloading the clonotyping data

made available here, the user can use an R script36 to gather some statistics and perform an initial analysis of the clonotypes. Note that the user should have a recent version of the coding

language R installed, and the appropriate packages up to date. The code also demonstrates performing sampling of the data, to prevent bias caused by samples with a larger number of

clonotypes. REFERENCES * Nikolich-Žugich, J., Slifka, M. K. & Messaoudi, I. The many important facets of T-cell repertoire diversity. _Nat. Rev. Immunol_. 4, 123–132. URL

https://www.nature.com/articles/nri1292. https://doi.org/10.1038/nri1292. Number: 2 Publisher: Nature Publishing Group (2004). * Adams, N. M., Grassmann, S. & Sun, J. C. Clonal expansion

of innate and adaptive lymphocytes. _Nat. Rev. Immunol_. 20, 694–707. Publisher: Nature Publishing Group (2020). * Heinzel, S., Marchingo, J. M., Horton, M. B. & Hodgkin, P. D. The

regulation of lymphocyte activation and proliferation. _Curr. Opin. Immunol_. 51, 32–38 URL https://www.sciencedirect.com/science/article/pii/S0952791517300948.

https://doi.org/10.1016/j.coi.2018.01.002 (2018). * Dash, P. _et al_. Quantifiable predictive features define epitope-specific T cell receptor repertoires. _Nat_. 547, 89–93. URL

https://www.nature.com/articles/nature22383. https://doi.org/10.1038/nature22383. Number: 7661 Publisher: Nature Publishing Group (2017). * Chi, X., Li, Y. & Qiu, X. V(D)J recombination,

somatic hypermutation and class switch recombination of immunoglobulins: mechanism and regulation. _Immunol_. 160, 233–247 URL https://onlinelibrary.wiley.com/doi/abs/10.1111/imm.13176.

https://doi.org/10.1111/imm.13176. _eprint: (2020). * Alt, F. W. _et al_. VDJ recombination. _Immunol. Today_ 13, 306–314. URL

https://www.sciencedirect.com/science/article/pii/0167569992900437. https://doi.org/10.1016/0167-5699(92)90043-7 (1992). * Darvin, P., Toor, S. M., Sasidharan Nair, V. & Elkord, E.

Immune checkpoint inhibitors: recent progress and potential biomarkers. _Exp. & Mol. Medicine_ 50, 1–11. URL https://www.nature.com/articles/s12276-018-0191-1.

https://doi.org/10.1038/s12276-018-0191-1. Number: 12 Publisher: Nature Publishing Group (2018). * Esfahani, K. _et al_. A Review of Cancer Immunotherapy: From the Past, to the Present, to

the Future. _Curr. Oncol_. 27, 87–97. URL https://www.mdpi.com/1718-7729/27/12/5223. https://doi.org/10.3747/co.27.5223. Number: s2 Publisher: Multidisciplinary Digital Publishing Institute

(2020). * Kidman, J. _et al_. Characteristics of TCR Repertoire Associated With Successful Immune Checkpoint Therapy Responses. _Front. Immunol_. 11. URL

https://www.frontiersin.org/article/10.3389/fimmu.2020.587014 (2020). * Ganesh, K. _et al_. Immunotherapy in colorectal cancer: rationale, challenges and potential. _Nat. Rev. Gastroenterol.

& Hepatol_. 16, 361–375. URL https://www.nature.com/articles/s41575-019-0126-x. https://doi.org/10.1038/s41575-019-0126-x. Number: 6 Publisher: Nature Publishing Group (2019). * Chae,

Y. K. _et al_. Current landscape and future of dual anti-CTLA4 and PD-1/PD-L1 blockade immunotherapy in cancer; lessons learned from clinical trials with melanoma and non-small cell lung

cancer (NSCLC). _J. for ImmunoTherapy Cancer_ 6, 39. URL https://doi.org/10.1186/s40425-018-0349-3 (2018). * Leach, D. R., Krummel, M. F. & Allison, J. P. Enhancement of Antitumor

Immunity by CTLA-4 Blockade. _Sci_. 271, 1734–1736. URL https://www.science.org/doi/abs/10.1126/science.271.5256.1734. https://doi.org/10.1126/science.271.5256.1734. Publisher: American

Association for the Advancement of Science (1996). * Waldman, A. D., Fritz, J. M. & Lenardo, M. J. A guide to cancer immunotherapy: from T cell basic science to clinical practice. _Nat.

Rev. Immunol_. 20, 651–668. URL https://www.nature.com/articles/s41577-020-0306-5. https://doi.org/10.1038/s41577-020-0306-5. Number: 11 Publisher: Nature Publishing Group (2020). *

Benichou, J., Ben-Hamo, R., Louzoun, Y. & Efroni, S. Rep-Seq: uncovering the immunological repertoire through next-generation sequencing. _Immunol._ 135, 183–191,

https://doi.org/10.1111/j.1365-2567.2011.03527.x (2012). Article  CAS  Google Scholar  * Warren, R. L. _et al_. Exhaustive T-cell repertoire sequencing of human peripheral blood samples

reveals signatures of antigen selection and a directly measured repertoire size of at least 1 million clonotypes. _Genome Res_. 21, 790–797. URL https://genome.cshlp.org/content/21/5/790.

https://doi.org/10.1101/gr.115428.110. Company: Cold Spring Harbor Laboratory Press Distributor: Cold Spring Harbor Laboratory Press Institution: Cold Spring Harbor Laboratory Press Label:

Cold Spring Harbor Laboratory Press Publisher: Cold Spring Harbor Lab (2011). * Georgiou, G. _et al_. The promise and challenge of high-throughput sequencing of the antibody repertoire.

_Nat. Biotechnol_. 32, 158–168. URL https://www.nature.com/articles/nbt.2782. https://doi.org/10.1038/nbt.2782. Number: 2 Publisher: Nature Publishing Group (2014). * Shugay, M. _et al_.

VDJtools: Unifying Post-analysis of T Cell Receptor Repertoires. _PLOS Comput. Biol_. 11, e1004503. URL https://journals.plos.org/ploscompbiol/article?id=10.1371/journal.pcbi.1004503.

https://doi.org/10.1371/journal.pcbi.1004503. Publisher: Public Library of Science (2015). * Bolotin, D. A. _et al_. MiXCR: software for comprehensive adaptive immunity profiling. _Nat.

Methods_ 12, 380–381, https://doi.org/10.1038/nmeth.3364 (2015). Article  CAS  PubMed  Google Scholar  * Arunkumar, M. & Zielinski, C. E. T-Cell Receptor Repertoire Analysis with

Computational Tools—An Immunologist’s Perspective. _Cells_ 10, 3582. URL https://www.mdpi.com/2073-4409/10/12/3582. https://doi.org/10.3390/cells10123582. Number: 12 Publisher:

Multidisciplinary Digital Publishing Institute (2021). * Robins, H. Immunosequencing: applications of immune repertoire deep sequencing. _Curr. Opin. Immunol_. 25, 646–652. URL

https://www.sciencedirect.com/science/article/pii/S0952791513001520. https://doi.org/10.1016/j.coi.2013.09.017 (2013). * Philip, H. _et al_. A T cell repertoire timestamp is at the core of

responsiveness to CTLA-4 blockade. _iScience_ 24, 102100. URL https://www.sciencedirect.com/science/article/pii/S2589004221000687. https://doi.org/10.1016/j.isci.2021.102100 (2021). * A

collection of SRA experiments of TCR sequencing of Mus Musculus - collectively referenced as SRP278268. _NCBI SRA_https://identifiers.org/insdc.sra:SRP278268 (2023). * Alpha Chain

Clonotypes. _figshare_ https://doi.org/10.6084/m9.figshare.19287506.v1 (2022). * Beta Chain Clonotypes. _figshare_ https://doi.org/10.6084/m9.figshare.19287512.v1 (2022). * Metadata for

Alpha chain clonotypes. _figshare_ https://doi.org/10.6084/m9.figshare.19312262.v2 (2022). * Metadata for Beta chain clonotypes. _figshare_ https://doi.org/10.6084/m9.figshare.19312265.v1

(2022). * Babraham Bioinformatics - FastQC A Quality Control tool for High Throughput Sequence Data. URL https://www.bioinformatics.babraham.ac.uk/projects/fastqc/. * Popov, A. _et al_.

immunomind/immunarch: Immunarch 0.6.9. _Zenodo_ https://doi.org/10.5281/zenodo.6599744 (2022). * Miho, E. _et al_. Computational Strategies for Dissecting the High-Dimensional Complexity of

Adaptive Immune Repertoires. _Front. Immunol_. 9. URL https://www.frontiersin.org/article/10.3389/fimmu.2018.00224 (2018). * Ye, J., Ma, N., Madden, T. L. & Ostell, J.M. IgBLAST: an

immunoglobulin variable domain sequence analysis tool. _Nucleic Acids Res_. 41, W34–W40. URL https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3692102/. https://doi.org/10.1093/nar/gkt382 (2013).

* Stark, R., Grzelak, M. & Hadfield, J. RNA sequencing: the teenage years. _Nat. Rev. Genet_. 20, 631–656. URL https://www.nature.com/articles/s41576-019-0150-2.

https://doi.org/10.1038/s41576-019-0150-2. Number: 11 Publisher: Nature Publishing Group (2019). * Morin, A. _et al_. Immunoseq: the identification of functionally relevant variants through

targeted capture and sequencing of active regulatory regions in human immune cells. _BMC Med. Genomics_ 9, 59. URL https://doi.org/10.1186/s12920-016-0220-7 (2016). * Vander Heiden, J. A.

_et al_. pRESTO: a toolkit for processing high-throughput sequencing raw reads of lymphocyte receptor repertoires. _Bioinforma_. 30, 1930–1932. URL

https://doi.org/10.1093/bioinformatics/btu138 (2014). * Pavlović, M. _et al_. The immuneML ecosystem for machine learning analysis of adaptive immune receptor repertoires. _Nat. Mach.

Intell_. 3, 936–944. URL https://www.nature.com/articles/s42256-021-00413-z. https://doi.org/10.1038/s42256-021-00413-z. Number: 11 Publisher: Nature Publishing Group (2021). * Clonotyping

of RNAseq. _figshare_ https://doi.org/10.6084/m9.figshare.19308539.v1 (2022). * R Markdown file to load and analyze clonotype data. _figshare_ https://doi.org/10.6084/m9.figshare.19342052.v3

(2022). Download references AUTHOR INFORMATION AUTHORS AND AFFILIATIONS * The Mina & Everard Goodman Faculty of Life Sciences, Bar-Ilan University, Ramat-Gan, Israel Tom Snir, Hagit

Philip, Miri Gordin, Alona Zilberberg & Sol Efroni Authors * Tom Snir View author publications You can also search for this author inPubMed Google Scholar * Hagit Philip View author

publications You can also search for this author inPubMed Google Scholar * Miri Gordin View author publications You can also search for this author inPubMed Google Scholar * Alona Zilberberg

View author publications You can also search for this author inPubMed Google Scholar * Sol Efroni View author publications You can also search for this author inPubMed Google Scholar

CONTRIBUTIONS A.Z. and S.E. designed the experiments, A.Z. and M.G. conducted the biological experiment, T.S. and H.P. performed the computational analysis. A.Z., S.E. and T.S. wrote the

manuscript. All authors reviewed the manuscript. CORRESPONDING AUTHOR Correspondence to Sol Efroni. ETHICS DECLARATIONS COMPETING INTERESTS The authors declare no competing interests.

ADDITIONAL INFORMATION PUBLISHER’S NOTE Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. RIGHTS AND PERMISSIONS OPEN

ACCESS This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format,

as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third

party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the

article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright

holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/. Reprints and permissions ABOUT THIS ARTICLE CITE THIS ARTICLE Snir, T., Philip, H., Gordin, M. _et

al._ The temporal behavior of the murine T cell receptor repertoire following Immunotherapy. _Sci Data_ 10, 108 (2023). https://doi.org/10.1038/s41597-023-01982-x Download citation *

Received: 08 August 2022 * Accepted: 24 January 2023 * Published: 23 February 2023 * DOI: https://doi.org/10.1038/s41597-023-01982-x SHARE THIS ARTICLE Anyone you share the following link

with will be able to read this content: Get shareable link Sorry, a shareable link is not currently available for this article. Copy to clipboard Provided by the Springer Nature SharedIt

content-sharing initiative