Defining the clonality of peripheral T cell lymphomas using RNA-seq

Motivation: In T‐cell lymphoma, malignant T cells arising from a founding clone share an identical T cell receptor (TCR) and can be identified by the over‐representation of this TCR relative to TCRs from the patient's repertoire of normal T cells. Here, we demonstrate that TCR information extracted from RNA‐seq data can provide a higher resolution view of peripheral T cell lymphomas (PTCLs) than that provided by conventional methods. Results: For 60 subjects with PTCL, flow cytometry/FACS was used to identify and sort aberrant T cell populations from diagnostic lymph node cell suspensions. For samples that did not appear to contain aberrant T cell populations, T helper (TH), T follicular helper (TFH) and cytotoxic T lymphocyte (CTL) subsets were sorted. RNA‐seq was performed on sorted T cell populations, and TCR alpha and beta chain sequences were extracted and quantified directly from the RNA‐seq data. 96% of the immunophenotypically aberrant samples had a dominant T cell clone readily identifiable by RNA‐seq. Of the samples where no aberrant population was found by flow cytometry, 80% had a dominant clone by RNA‐seq. This demonstrates the increased sensitivity and diagnostic ability of RNA‐seq over flow cytometry and shows that the presence of a normal immunophenotype does not exclude clonality. Availability and Implementation: R scripts used in the processing of the data are available online at https://www.github.com/scottdbrown/RNAseq‐TcellClonality Contacts: rholt@bcgsc.ca or ksavage@bccancer.bc.ca Supplementary information: Supplementary data are available at Bioinformatics online.

[1]  Shulin Li,et al.  Nonsense surveillance in lymphocytes? , 1998, Immunity.

[2]  Sang Joon Kim,et al.  A Mathematical Theory of Communication , 2006 .

[3]  T. Szczepański,et al.  Heteroduplex PCR analysis of rearranged T cell receptor genes for clonality assessment in suspect T cell proliferations , 1997, Leukemia.

[4]  Steven J. M. Jones,et al.  Evolution of an adenocarcinoma in response to selection by targeted kinase inhibitors , 2010, Genome Biology.

[5]  Rei Watanabe,et al.  TCR sequencing facilitates diagnosis and identifies mature T cells as the cell of origin in CTCL , 2015, Science Translational Medicine.

[6]  M Hummel,et al.  Design and standardization of PCR primers and protocols for detection of clonal immunoglobulin and T-cell receptor gene recombinations in suspect lymphoproliferations: Report of the BIOMED-2 Concerted Action BMH4-CT98-3936 , 2003, Leukemia.

[7]  Steven L Salzberg,et al.  Fast gapped-read alignment with Bowtie 2 , 2012, Nature Methods.

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

[9]  David Wu,et al.  High-Throughput Sequencing Detects Minimal Residual Disease in Acute T Lymphoblastic Leukemia , 2012, Science Translational Medicine.

[10]  R. Gascoyne,et al.  Characterization of peripheral T-cell lymphomas in a single North American institution by the WHO classification. , 2004, Annals of oncology : official journal of the European Society for Medical Oncology.

[11]  C. E. SHANNON,et al.  A mathematical theory of communication , 1948, MOCO.

[12]  D. Weisenburger,et al.  International peripheral T-cell and natural killer/T-cell lymphoma study: pathology findings and clinical outcomes. , 2008, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[13]  John R. Krause,et al.  WHO Classification of Tumours of Haematopoietic and Lymphoid Tissues: An Overview , 2009 .

[14]  Steven J. M. Jones,et al.  Lessons learned from the application of whole-genome analysis to the treatment of patients with advanced cancers , 2015, Cold Spring Harbor molecular case studies.

[15]  Colin N. Dewey,et al.  RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome , 2011, BMC Bioinformatics.

[16]  Ken I. Mills,et al.  Design and standardization of PCR primers and protocols for detection of clonal immunoglobulin and T-cell receptor gene recombinations in suspect lymphoproliferations. Section 9. t(14;18) with BCL2-IGH rearrangement. , 2003 .

[17]  Nitin J. Karandikar,et al.  Flow cytometric features of angioimmunoblastic T‐cell lymphoma , 2006, Cytometry. Part B, Clinical cytometry.

[18]  Robert A Holt,et al.  Sequence analysis of T-cell repertoires in health and disease , 2013, Genome Medicine.

[19]  R. Holt,et al.  Profiling the T-cell receptor beta-chain repertoire by massively parallel sequencing. , 2009, Genome research.

[20]  Mikhail Shugay,et al.  MiTCR: software for T-cell receptor sequencing data analysis , 2013, Nature Methods.

[21]  B. Sleckman,et al.  Allelic exclusion at the TCRβ locus , 2002 .

[22]  Scott D. Brown,et al.  Profiling tissue-resident T cell repertoires by RNA sequencing , 2015, Genome Medicine.