{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,3,12]],"date-time":"2026-03-12T00:36:37Z","timestamp":1773275797734,"version":"3.50.1"},"reference-count":24,"publisher":"Oxford University Press (OUP)","issue":"4","content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2009,2,15]]},"abstract":"<jats:title>Abstract<\/jats:title>\n               <jats:p>Motivation: T-cell receptor (TCR) diversity in peripheral blood has not yet been fully profiled with sequence level resolution. Each T-cell clonotype expresses a unique receptor, generated by somatic recombination of TCR genes and the enormous potential for T-cell diversity makes repertoire analysis challenging. We developed a sequencing approach and assembly software (immuno-SSAKE or iSSAKE) for profiling T-cell metagenomes using short reads from the massively parallel sequencing platforms.<\/jats:p>\n               <jats:p>Results: Models of sequence diversity for the TCR \u03b2-chain CDR3 region were built using empirical data and used to simulate, at random, distinct TCR clonotypes at 1\u201320 p.p.m. Using simulated TCR\u03b2 (sTCR\u03b2) sequences, we randomly created 20 million 36 nt reads having 1\u20132% random error, 20 million 42 or 50 nt reads having 1% random error and 20 million 36 nt reads with 1% error modeled on real short read data. Reads aligning to the end of known TCR variable (V) genes and having consecutive unmatched bases in the adjacent CDR3 were used to seed iSSAKE de novo assemblies of CDR3. With assembled 36 nt reads, we detect over 51% and 63% of rare (1 p.p.m.) clonotypes using a random or modeled error distribution, respectively. We detect over 99% of more abundant clonotypes (6 p.p.m. or higher) using either error distribution. Longer reads improve sensitivity, with assembled 42 and 50 nt reads identifying 82.0% and 94.7% of rare 1 p.p.m. clonotypes, respectively. Our approach illustrates the feasibility of complete profiling of the TCR repertoire using new massively parallel short read sequencing technology.<\/jats:p>\n               <jats:p>Availability: \u00a0ftp:\/\/ftp.bcgsc.ca\/supplementary\/iSSAKE<\/jats:p>\n               <jats:p>Contact: \u00a0rwarren@bcgsc.ca<\/jats:p>\n               <jats:p>Supplementary information: \u00a0Supplementary methods and data are available at Bioinformatics online.<\/jats:p>","DOI":"10.1093\/bioinformatics\/btp010","type":"journal-article","created":{"date-parts":[[2009,1,10]],"date-time":"2009-01-10T01:34:33Z","timestamp":1231551273000},"page":"458-464","source":"Crossref","is-referenced-by-count":37,"title":["Profiling model T-cell metagenomes with short reads"],"prefix":"10.1093","volume":"25","author":[{"given":"Ren\u00e9 L.","family":"Warren","sequence":"first","affiliation":[{"name":"1 BC Cancer Agency, Michael Smith Genome Sciences Centre, 675 West 10th Avenue, Vancouver, BC V5Z 1L3 Canada and 2BC Cancer Agency, Deeley Research Centre, 2410 Lee Ave, Victoria, BC V8R 6V5 Canada"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Brad H.","family":"Nelson","sequence":"additional","affiliation":[{"name":"1 BC Cancer Agency, Michael Smith Genome Sciences Centre, 675 West 10th Avenue, Vancouver, BC V5Z 1L3 Canada and 2BC Cancer Agency, Deeley Research Centre, 2410 Lee Ave, Victoria, BC V8R 6V5 Canada"}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Robert A.","family":"Holt","sequence":"additional","affiliation":[{"name":"1 BC Cancer Agency, Michael Smith Genome Sciences Centre, 675 West 10th Avenue, Vancouver, BC V5Z 1L3 Canada and 2BC Cancer Agency, Deeley Research Centre, 2410 Lee Ave, Victoria, BC V8R 6V5 Canada"}],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"286","published-online":{"date-parts":[[2009,1,9]]},"reference":[{"key":"2023013110100545900_B1","doi-asserted-by":"crossref","first-page":"958","DOI":"10.1126\/science.286.5441.958","article-title":"A direct estimate of the human \u03b1\u03b2 T cell receptor diversity","volume":"286","author":"Arstila","year":"1999","journal-title":"Science"},{"key":"2023013110100545900_B2","doi-asserted-by":"crossref","first-page":"433","DOI":"10.1517\/14622416.5.4.433","article-title":"Solexa Ltd","volume":"5","author":"Bennett","year":"2004","journal-title":"Pharmacogenomics"},{"key":"2023013110100545900_B3","doi-asserted-by":"crossref","first-page":"523","DOI":"10.1146\/annurev.immunol.16.1.523","article-title":"Ligand recognition by alpha beta T cell receptors","volume":"16","author":"Davis","year":"1998","journal-title":"Annu. Rev. Immunol."},{"key":"2023013110100545900_B4","doi-asserted-by":"crossref","first-page":"1697","DOI":"10.1101\/gr.6435207","article-title":"SHARCGS, a fast and highly accurate short-read assembly algorithm for de novo genomic sequencing","volume":"17","author":"Dohm","year":"2007","journal-title":"Genome Res."},{"key":"2023013110100545900_B5","doi-asserted-by":"crossref","first-page":"D707","DOI":"10.1093\/nar\/gkm988","article-title":"Ensembl 2008","volume":"36","author":"Flicek","year":"2008","journal-title":"Nucleic Acids Res"},{"key":"2023013110100545900_B6","doi-asserted-by":"crossref","first-page":"387","DOI":"10.1038\/310387a0","article-title":"Genomic organization and sequence of T-cell receptor beta-chain constant- and joining-region genes","volume":"310","author":"Gascoigne","year":"1984","journal-title":"Nature"},{"key":"2023013110100545900_B7","doi-asserted-by":"crossref","first-page":"5109","DOI":"10.4049\/jimmunol.152.10.5109","article-title":"Circulating T cell repertoire complexity in normal individuals and bone marrow recipients analyzed by CDR3 size spectratyping. Correlation with immune status","volume":"152","author":"Gorski","year":"1994","journal-title":"J. Immunol."},{"key":"2023013110100545900_B8","doi-asserted-by":"crossref","first-page":"149","DOI":"10.1038\/308149a0","article-title":"Isolation of cDNA clones encoding T cell-specific membrane-associated proteins","volume":"308","author":"Hedrick","year":"1984","journal-title":"Nature"},{"key":"2023013110100545900_B9","doi-asserted-by":"crossref","first-page":"1851","DOI":"10.1101\/gr.078212.108","article-title":"Mapping short DNA sequencing reads and calling variants using mapping quality scores","volume":"18","author":"Heng","year":"2008","journal-title":"Genome Res."},{"key":"2023013110100545900_B10","doi-asserted-by":"crossref","first-page":"839","DOI":"10.1101\/gr.073262.107","article-title":"The new paradigm of flow cell sequencing","volume":"18","author":"Holt","year":"2008","journal-title":"Genome Res."},{"key":"2023013110100545900_B11","doi-asserted-by":"crossref","first-page":"719","DOI":"10.1002\/j.1460-2075.1988.tb02868.x","article-title":"Variable region genes in the human T-cell rearranging gamma (TRG) locus: V-J junction and homology with the mouse genes","volume":"7","author":"Huck","year":"1988","journal-title":"EMBO J"},{"key":"2023013110100545900_B12","volume-title":"Immunobiology.","author":"Janeway","year":"2001","edition":"6th edn"},{"key":"2023013110100545900_B13","doi-asserted-by":"crossref","first-page":"2942","DOI":"10.1093\/bioinformatics\/btm451","article-title":"Extending assembly of short DNA sequences to handle error","volume":"23","author":"Jeck","year":"2007","journal-title":"Bioinformatics"},{"key":"2023013110100545900_B14","doi-asserted-by":"crossref","first-page":"421","DOI":"10.1038\/310421a0","article-title":"Localization of a T-cell receptor diversity-region element","volume":"310","author":"Kavaler","year":"1984","journal-title":"Nature"},{"key":"2023013110100545900_B15","doi-asserted-by":"crossref","first-page":"5836","DOI":"10.1073\/pnas.81.18.5836","article-title":"Cloning of terminal transferase cDNA by antibody screening","volume":"81","author":"Landau","year":"1984","journal-title":"Proc. Natl Acad. Sci. USA"},{"key":"2023013110100545900_B16","volume-title":"The T cell Receptor Facts-Book.","author":"Lefranc","year":"2001"},{"key":"2023013110100545900_B17","doi-asserted-by":"crossref","first-page":"376","DOI":"10.1038\/nature03959","article-title":"Genome sequencing in microfabricated high-density picolitre reactors","volume":"437","author":"Margulies","year":"2005","journal-title":"Nature"},{"key":"2023013110100545900_B18","doi-asserted-by":"crossref","first-page":"820","DOI":"10.1016\/j.bbrc.2008.01.011","article-title":"Comprehensive analysis of the functional TCR repertoire at the single-cell level","volume":"367","author":"Ozawa","year":"2008","journal-title":"Biochem. Biophys. Res. Commun."},{"key":"2023013110100545900_B19","doi-asserted-by":"crossref","first-page":"4319","DOI":"10.1073\/pnas.90.9.4319","article-title":"The sizes of the CDR3 hypervariable regions of the murine T-cell receptor beta chains vary as a function of the recombined germ-line segments","volume":"90","author":"Pannetier","year":"1993","journal-title":"Proc. Natl Acad. Sci. USA"},{"key":"2023013110100545900_B20","doi-asserted-by":"crossref","first-page":"6611","DOI":"10.4049\/jimmunol.180.10.6611","article-title":"Administration of PLP139-151 primes T cells distinct from those spontaneously responsive in vitro to this antigen","volume":"180","author":"Penitente","year":"2008","journal-title":"J. Immunol"},{"key":"2023013110100545900_B21","doi-asserted-by":"crossref","first-page":"288","DOI":"10.1038\/280288a0","article-title":"Sequences at the somatic recombination sites of immunoglobulin light-chain genes","volume":"280","author":"Sakano","year":"1979","journal-title":"Nature"},{"key":"2023013110100545900_B22","doi-asserted-by":"crossref","first-page":"31","DOI":"10.1186\/1471-2105-6-31","article-title":"Automated generation of heuristics for biological sequence comparison","volume":"6","author":"Slater","year":"2005","journal-title":"BMC Bioinformatics"},{"key":"2023013110100545900_B23","doi-asserted-by":"crossref","first-page":"500","DOI":"10.1093\/bioinformatics\/btl629","article-title":"Assembling millions of short DNA sequences using SSAKE","volume":"23","author":"Warren","year":"2007","journal-title":"Bioinformatics"},{"key":"2023013110100545900_B24","doi-asserted-by":"crossref","first-page":"314","DOI":"10.1038\/labinvest.3700381","article-title":"High throughput analysis of TCR-b rearrangement and gene expression in single T cells","volume":"86","author":"Zhou","year":"2006","journal-title":"Lab. Invest."}],"container-title":["Bioinformatics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/25\/4\/458\/48984260\/bioinformatics_25_4_458.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article-pdf\/25\/4\/458\/48984260\/bioinformatics_25_4_458.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2023,1,31]],"date-time":"2023-01-31T19:42:48Z","timestamp":1675194168000},"score":1,"resource":{"primary":{"URL":"https:\/\/academic.oup.com\/bioinformatics\/article\/25\/4\/458\/251049"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2009,1,9]]},"references-count":24,"journal-issue":{"issue":"4","published-print":{"date-parts":[[2009,2,15]]}},"URL":"https:\/\/doi.org\/10.1093\/bioinformatics\/btp010","relation":{},"ISSN":["1367-4811","1367-4803"],"issn-type":[{"value":"1367-4811","type":"electronic"},{"value":"1367-4803","type":"print"}],"subject":[],"published-other":{"date-parts":[[2009,2,15]]},"published":{"date-parts":[[2009,1,9]]}}}