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

The T cell receptor (TCR) repertoire can provide a personalized biomarker for infectious and non-infectious diseases. We describe a protocol for amplifying, sequencing, and analyzing TCRs which is robust, sensitive, and versatile. The key experimental step is ligation of a single-stranded oligonucleotide to the 3′ end of the TCR cDNA. This allows amplification of all possible rearrangements using a single set of primers per locus. It also introduces a unique molecular identifier to label each starting cDNA molecule. This molecular identifier is used to correct for sequence errors and for effects of differential PCR amplification efficiency, thus producing more accurate measures of the true TCR frequency within the sample. This integrated experimental and computational pipeline is applied to the analysis of human memory and naive subpopulations, and results in consistent measures of diversity and inequality. After error correction, the distribution of TCR sequence abundance in all subpopulations followed a power law over a wide range of values. The power law exponent differed between naïve and memory populations, but was consistent between individuals. The integrated experimental and analysis pipeline we describe is appropriate to studies of T cell responses in a broad range of physiological and pathological contexts.

[1]  John Shawe-Taylor,et al.  Computational analysis of stochastic heterogeneity in PCR amplification efficiency revealed by single molecule barcoding , 2015, Scientific Reports.

[2]  Christian Stemberger,et al.  A single naive CD8+ T cell precursor can develop into diverse effector and memory subsets. , 2007, Immunity.

[3]  A. Lanzavecchia,et al.  Delineation of several DR-restricted tetanus toxin T cell epitopes. , 1989, Journal of immunology.

[4]  Alfred V. Aho,et al.  Efficient string matching , 1975, Commun. ACM.

[5]  J. Mallet,et al.  Oligodeoxyribonucleotide ligation to single-stranded cDNAs: a new tool for cloning 5' ends of mRNAs and for constructing cDNA libraries by in vitro amplification. , 1991, Nucleic acids research.

[6]  Arne N. Akbar,et al.  Cytomegalovirus-Specific CD4+ T Cells in Healthy Carriers Are Continuously Driven to Replicative Exhaustion1 , 2005, The Journal of Immunology.

[7]  Mikhail Shugay,et al.  MiXCR: software for comprehensive adaptive immunity profiling , 2015, Nature Methods.

[8]  Benjamin Chain,et al.  High-throughput sequencing of the T-cell receptor repertoire: pitfalls and opportunities , 2017, Briefings Bioinform..

[9]  M. Ivanchenko,et al.  Assessing T Cell Clonal Size Distribution: A Non-Parametric Approach , 2014, PloS one.

[10]  C. Vaquero,et al.  Correlation between up-regulation of lymphokine mRNA and down-regulation of TcR, CD4, CD8 and lck mRNA as shown by the effect of CsA on activated T lymphocytes. , 1992, Biochemical and biophysical research communications.

[11]  T. Mak,et al.  Genes of the T-cell antigen receptor in normal and malignant T cells. , 1987, Annual review of immunology.

[12]  G. Lythe,et al.  A new mechanism shapes the naïve CD8+ T cell repertoire: the selection for full diversity , 2017, Molecular Immunology.

[13]  B. Chain,et al.  The sequence of sequencers: The history of sequencing DNA , 2016, Genomics.

[14]  Abigail Wacher,et al.  Comprehensive assessment of T-cell receptor beta-chain diversity in alphabeta T cells. , 2009, Blood.

[15]  Liran Carmel,et al.  System-wide Analysis of the T Cell Response. , 2016, Cell reports.

[16]  Thierry Mora,et al.  The Past, Present, and Future of Immune Repertoire Biology – The Rise of Next-Generation Repertoire Analysis , 2013, Front. Immunol..

[17]  Michael J. Bevan,et al.  The precursors of memory: models and controversies , 2009, Nature Reviews Immunology.

[18]  Thierry Mora,et al.  Fluctuating fitness shapes the clone-size distribution of immune repertoires , 2015, Proceedings of the National Academy of Sciences.

[19]  Olga V. Britanova,et al.  Preparing Unbiased T-Cell Receptor and Antibody cDNA Libraries for the Deep Next Generation Sequencing Profiling , 2013, Front. Immunol..

[20]  Frank Baas,et al.  Human T-cell memory consists mainly of unexpanded clones. , 2010, Immunology letters.

[21]  H. Bauke.  Parameter estimation for power-law distributions by maximum likelihood methods , 2007, 0704.1867.

[22]  J. D. Burgos,et al.  Zipf-scaling behavior in the immune system. , 1996, Bio Systems.

[23]  Olivier Elemento,et al.  Single-cell TCRseq: paired recovery of entire T-cell alpha and beta chain transcripts in T-cell receptors from single-cell RNAseq , 2016, Genome Medicine.

[24]  C. Thompson,et al.  Transcription of T cell antigen receptor genes is induced by protein kinase C activation. , 1988, Journal of immunology.

[25]  Rob J. De Boer,et al.  RTCR: a pipeline for complete and accurate recovery of T cell repertoires from high throughput sequencing data , 2016, Bioinform..

[26]  H. Robins.  Immunosequencing: applications of immune repertoire deep sequencing. , 2013, Current opinion in immunology.

[27]  A. Cumano,et al.  The expansion and selection of T cell receptor αβ intestinal intraepithelial T cell clones , 1996 .

[28]  A. Sewell,et al.  αβ T cell receptors as predictors of health and disease , 2015, Cellular and Molecular Immunology.

[29]  Marion C Lanteri,et al.  Human memory T cells with a naïve phenotype accumulate with aging and respond to persistent viruses , 2016, Nature Immunology.

[30]  Richard A. Olshen,et al.  Diversity and clonal selection in the human T-cell repertoire , 2014, Proceedings of the National Academy of Sciences.

[31]  Wilfred Ndifon,et al.  Chromatin conformation governs T-cell receptor Jβ gene segment usage , 2012, Proceedings of the National Academy of Sciences.

[32]  C. Desmarais,et al.  Fractal organization of the human T cell repertoire in health and after stem cell transplantation. , 2013, Biology of blood and marrow transplantation : journal of the American Society for Blood and Marrow Transplantation.

[33]  Anneliese O. Speak,et al.  T cell fate and clonality inference from single cell transcriptomes , 2016, Nature Methods.

[34]  H. Diao,et al.  Current status and recent advances of next generation sequencing techniques in immunological repertoire , 2016, Genes and Immunity.

[35]  Baback Gharizadeh,et al.  High throughput sequencing reveals a complex pattern of dynamic interrelationships among human T cell subsets , 2010, Proceedings of the National Academy of Sciences.

[36]  B. Chain,et al.  In Vivo Molecular Dissection of the Effects of HIV-1 in Active Tuberculosis , 2016, PLoS pathogens.

[37]  R. White,et al.  High-Throughput Sequencing of the Zebrafish Antibody Repertoire , 2009, Science.

[38]  Grant Lythe,et al.  How many TCR clonotypes does a body maintain? , 2016, Journal of theoretical biology.

[39]  S. J. Griffiths,et al.  Age-Associated Increase of Low-Avidity Cytomegalovirus-Specific CD8+ T Cells That Re-Express CD45RA , 2013, The Journal of Immunology.

[40]  D. Tessier,et al.  Ligation of single-stranded oligodeoxyribonucleotides by T4 RNA ligase. , 1986, Analytical biochemistry.

[41]  Mikhail Shugay,et al.  Towards error-free profiling of immune repertoires , 2014, Nature Methods.

[42]  B. Pulendran,et al.  Ligation-anchored PCR: a simple amplification technique with single-sided specificity. , 1992, Proceedings of the National Academy of Sciences of the United States of America.

[43]  Alvin Shi,et al.  Sequence and Structural Analyses Reveal Distinct and Highly Diverse Human CD8+ TCR Repertoires to Immunodominant Viral Antigens. , 2017, Cell reports.

[44]  A. Khoruts,et al.  Naïve and Memory CD4+ T Cell Survival Controlled by Clonal Abundance , 2006, Science.

[45]  R. Hagedoorn,et al.  Efficiency of T-cell receptor expression in dual-specific T cells is controlled by the intrinsic qualities of the TCR chains within the TCR-CD3 complex. , 2007, Blood.

[46]  B. Chain,et al.  Transcriptional profiling of innate and adaptive human immune responses to mycobacteria in the tuberculin skin test , 2011, European journal of immunology.

[47]  L. Battistini,et al.  Reversible Senescence in Human CD4+CD45RA+CD27− Memory T Cells , 2011, Journal of Immunology.

[48]  John Shawe-Taylor,et al.  Decombinator: a tool for fast, efficient gene assignment in T-cell receptor sequences using a finite state machine , 2013, Bioinform..

[49]  Emmanuel Beaudoing,et al.  Size Estimate of the αβ TCR Repertoire of Naive Mouse Splenocytes1 , 2000, The Journal of Immunology.

[50]  G. Foucras,et al.  Tracking T cell clonotypes in complex T lymphocyte populations by real-time quantitative PCR using fluorogenic complementarity-determining region-3-specific probes. , 2002, Journal of immunological methods.

[51]  Thomas Höfer,et al.  Disparate Individual Fates Compose Robust CD8+ T Cell Immunity , 2013, Science.

[52]  Mark E. J. Newman,et al.  Power-Law Distributions in Empirical Data , 2007, SIAM Rev..

[53]  D. Laydon,et al.  Estimating T-cell repertoire diversity: limitations of classical estimators and a new approach , 2015, Philosophical Transactions of the Royal Society B: Biological Sciences.

[54]  C. Vaquero,et al.  Transcriptional and post‐transcriptional regulation of TcR, CD4 and CD8 gene expression during activation of normal human T lymphocytes. , 1990, The EMBO journal.