Identification of cell-type-specific mutations in nodal T-cell lymphomas

<jats:title>Abstract</jats:title><jats:p>Recent genetic analysis has identified frequent mutations in <jats:italic>ten-eleven translocation 2</jats:italic> (<jats:italic>TET2</jats:italic>), <jats:italic>DNA methyltransferase 3A</jats:italic> (<jats:italic>DNMT3A</jats:italic>), <jats:italic>isocitrate dehydrogenase 2</jats:italic> (<jats:italic>IDH2</jats:italic>) and <jats:italic>ras homolog family member A</jats:italic> (<jats:italic>RHOA</jats:italic>) in nodal T-cell lymphomas, including angioimmunoblastic T-cell lymphoma and peripheral T-cell lymphoma, not otherwise specified. We examined the distribution of mutations in these subtypes of mature T-/natural killer cell neoplasms to determine their clonal architecture. Targeted sequencing was performed for 71 genes in tumor-derived DNA of 87 cases. The mutations were then analyzed in a programmed death-1 (PD1)-positive population enriched with tumor cells and CD20-positive B cells purified by laser microdissection from 19 cases. <jats:italic>TET2</jats:italic> and <jats:italic>DNMT3A</jats:italic> mutations were identified in both the PD1+ cells and the CD20+ cells in 15/16 and 4/7 cases, respectively. All the <jats:italic>RHOA</jats:italic> and <jats:italic>IDH2</jats:italic> mutations were confined to the PD1+ cells, indicating that some, including <jats:italic>RHOA</jats:italic> and <jats:italic>IDH2</jats:italic> mutations, being specific events in tumor cells. Notably, we found that all <jats:italic>NOTCH1</jats:italic> mutations were detected only in the CD20+ cells. In conclusion, we identified both B- as well as T-cell-specific mutations, and mutations common to both T and B cells. These findings indicate the expansion of a clone after multistep and multilineal acquisition of gene mutations.</jats:p>

section_content: Nodal T-cell lymphomas are subtypes of mature T-/natural killer-cell neoplasms, including angioimmunoblastic T-cell lymphoma (AITL); nodal peripheral T-cell lymphoma (PTCL) with T follicular helper (TFH) phenotype; peripheral T-cell lymphoma, not otherwise specified (PTCL-NOS), and follicular T-cell lymphoma. 2] [3] AITL is characterized by specific clinical features, including generalized lymphadenopathy, high fever, skin rash and autoimmune-like manifestations. AITL tumor cells share characteristics with TFH cells, expressing B-cell lymphoma protein 6, a transcription factor; C-C motif chemokine receptor 5, a chemokine receptor; C-X-C motif ligand 13, a chemokine; and programmed death-1 (PD1), a member of the CD28 costimulatory membrane receptor family. 4, 5 AITL tissues display prominent infiltration of inflammatory cells, follicular dendritic cell meshwork formation and branching vascular structures. Some nodal T-cell lymphomas exhibit several features reminiscent of AITL, although they do not show the typical morphology of AITL (nodal PTCL with TFH phenotype). 6, 7 The massive infiltration of inflammatory cells in AITL has been explained by cytokines and chemokines being released from TFH-like tumor cells. 4 ] [10] [11] [12] [13] Mutations in TET2 encoding a methylcytosine dioxygenase and those in DNMT3A encoding a DNA methyltransferase presumably result in epigenetic abnormalities in nodal T-cell lymphomas. IDH2 mutations also affect epigenetic modifications by inhibiting TET and histone demethylation enzymes through production of 2hydroxyglutarate. 14 Mutations in RHOA encoding a small GTPase are almost always located at the hotspot site, resulting in conversion from glycine to valine at the seventeenth position of the RHOA protein (G17V RHOA mutation). The G17V RHOA mutants could not be converted to an active GTP-bound form, although the downstream signaling of the G17V RHOA mutants in nodal T-cell lymphomas development has yet to be clarified. 8, 9, 13 ET2 and DNMT3A mutations are proposed to arise in hematopoietic stem/progenitors upstream of T-lineage commitment. 6] [17] In contrast, the origins of the G17V RHOA mutation remain to be elucidated: it may be a tumor-specific event, considering that the allele frequencies of G17V RHOA mutations were lower than those of TET2 mutations and that G17V RHOA mutations were found in only CD4+T lymphocytes in 1 AITL and 1 PTCL-NOS case. 8 ere we describe the clonal architecture of nodal T-cell lymphomas by determining the distribution of mutations in enriched tumor cells and infiltrated B cells.
section_name: INTRODUCTION

section_name: MATERIALS AND METHODS

section_content: Samples, obtained from 87 patients (Supplementary Table S1 ) with AITL (n = 48), nodal PTCL with TFH phenotype (n = 5) and either of PTCL-NOS or of nodal PTCL with TFH phenotype (PTCL-NOS/nodal PTCL with TFH phenotype, n = 34), 5 were used after approval was obtained from the local ethics committees of all the participating institutes. 

Genomic DNA was extracted from 56 fresh frozen samples using the Puregene DNA Blood Kit (Qiagen, Hilden, Germany) and 31 periodate-lysineparaformaldehyde (PLP)-fixed frozen samples using the QIAamp DNA FFPE Tissue Kit (Qiagen).
section_name: Patients and samples

section_content: Targeted sequencing was performed for 71 genes, which are listed in Supplementary Table S2. Sixty-one of the genes were previously screened by whole-exome sequencing, 8 while 6 were the family genes of those whose mutations were identified by the whole-exome sequencing. 9] [20] [21] [22] [23] All the exons of the selected genes were captured by use of a SureSelect Target Enrichment Kit (Agilent, Santa Clara, CA, USA) and then massively sequenced using HiSeq2000 (Illumina, Santa Clara, CA, USA). For each sample, all the sequencing reads were aligned to hg19 using BWA version 0. 5. 8 with default parameters. After all the duplicated reads and the lowquality reads and bases were removed, the allele frequencies of singlenucleotide variants and indels at each genomic position were calculated by enumerating the relevant reads using SAMtools (http://www. htslib. org). Initially, all the variants showing allele frequencies 40. 02 were extracted and annotated using ANNOVAR 24, 25 for further consideration, if they were found in 46 reads of 410 total reads and appeared in both the positiveand the negative-strand reads. All synonymous variants and known singlenucleotide polymorphisms in public and private databases, including dbSNP131, the 1000 genomes project as of 2012/05/21 and our in-house database, were removed. To exclude germline variants, nonsynonymous variants were excluded when the allele frequencies were from 0. 45 to 0. 55. Candidate mutations were validated by amplicon-based deep sequencing using Ion PGM (Life Technologies, Carlsbad, CA, USA) and/or Sanger sequencing (see below). 

In the cohort of 87 cases, 79 were analyzed for RHOA, TET2, DNMT3A and IDH2 mutations, and the results of this analysis were described in the previous paper. 7 Now, eight were new cases. We re-analyzed all the 87 samples for targeted sequencing of 71 genes.
section_name: Targeted sequencing

section_content: The libraries were prepared using the Ion Plus Fragment Library Kit according to the protocol for preparing short amplicon libraries (Life Technologies). Briefly, PCR amplicons were ligated to the barcode adapters and P1 adapters and then amplified. The amplified libraries were quantitated by quantitative PCR with the Ion Library Quantitation Kit according to the manufacturer's instructions (Life Technologies). The libraries were then subjected to deep sequencing on the Ion Torrent PGM platform according to the standard protocol for 300 base-pair single-end reads (Life Technologies). The data were analyzed using Variant Caller 3. 4 (Life Technologies). 

Immunohistochemistry PLP-fixed frozen samples were cut in a cryostat at -22 °C into 5-μm sections and mounted on PEN-Membrane slides (Leica, Wetzlar, Germany). The tissue sections were stained with mouse anti-human PD1 (NAT105 ab52587, Abcam, Cambridge, UK) and anti-human CD20cy (clone L26, Dako, Michigan, MI, USA) antibodies, diluted 1:2000 and 1:1000, respectively, and detected by use of the Envision + Dual Link System-HRP (Dako). The tissue sections were then counterstained with hematoxylin (Mayer's hematoxylin, Muto Pure Chemical, Tokyo, Japan) for 20 s at room temperature. After staining, tissue sections were dehydrated with ethanol and dried at room temperature before laser microdissection (LMD).
section_name: Amplicon-based sequencing

section_content: Nineteen of the 87 cases (13 AITL, 1 nodal PTCL with TFH phenotype and PTCL-NOS/nodal PTCL with TFH phenotype) were analyzed by LMD, which was performed using LMD7000 (Leica). The cells being positive for either PD1 or CD20 were dissected and collected into 0. 2-ml PCR tubes (Takara, Shiga, Japan) containing 20 μl of distilled water. Stained cells at approximately 100 000 μm 2 were dissected and collected for each sample. Genomic DNA was extracted using the QIAamp DNA PFFE Tissue Kit (Qiagen) following the manufacturer's protocol. Then 1 μl of DNA was used for PCR under the following conditions: 95 °C for 15 min, 60 °C for 4 min, 72 °C for 4 min, 35 to 40 cycles at 95 °C for 1 min, 60 °C for 1 min, 72 °C for 1 min, and 72 °C for 10 min using the AmpliTaq Gold 360 Kit (Applied Biosystems, Foster City, CA, USA) with each primer set (Supplementary Table S6 ). PCR amplicons were used for amplicon-based sequencing and Sanger sequencing.
section_name: LMD, DNA extraction and PCR

section_content: Multiplex PCR assays were used to detect the clonality of B cells according to the European BIOMED-2 collaborative study. 26 PCR products migrating at the expected size were extracted and sequenced using the Sanger method. Subcloning was performed if the Sanger sequencing indicated a polyclonal background by use of the pGEM-T Easy Vector System I (Promega, Madison, WI, USA). At least 12 colonies were picked up and sequenced to confirm the clonal expansion. The sequence results were analyzed using the IMGT tools 27 and aligned to the closest match with the germline IGHV segment. 

Sequencing results with a germline identity of o98% were regarded as mutated and vice versa according to previous study. 28
section_name: IgH gene rearrangement analysis and subcloning of the PCR product

section_content: Novel recurrent mutations in nodal T-cell lymphomas Targeted sequencing for 71 genes was performed in 87 samples (Supplementary Table S1 ), including AITL (n = 48), nodal PTCL with TFH phenotype (n = 5) and PTCL-NOS/nodal PTCL with TFH phenotype (n = 34). TET2, DNMT3A, RHOA and IDH2 mutations were identified in 60 (68. 7%), 23 (26. 4%), 41 (47. 1%) and 13 (14. 9%) of 87 cases, respectively (Figure 1, Table 1, Supplementary Table S3 ). The mutational profiles of these 4 genes in 79 of the 87 samples are described elsewhere. 7 hirty-four novel recurrent mutations were identified in 13 of the 71 (18. 3%) genes and in 24 of the 87 (26. 4%) cases (Figure 1, Table 1 and Supplementary Table S4 ). Mutations in genes associated with lymphoid malignancies, for example, Notch homolog 1, translocation-associated (NOTCH1), β2 microglobulin (B2M) and mixed-lineage leukemia 2 (MLL2) were identified in 3, 2 and 2 cases, respectively. Mutations in FAT atypical cadherin 2 (FAT2), a gene associating with several cancers, and those in TET3, a member of the TET gene family, were identified in 3 and 2 cases, respectively. 

Nineteen of the 87 samples were analyzed by the use of LMD (Figure 1 ). The frequencies of TET2, DNMT3A, RHOA and IDH2 mutations in the laser-microdissected samples were similar to those found in the entire cohort (TET2 mutations, 16 S5 ). NOTCH1 and COL19A1 mutations were identified in 3 and 2 cases, respectively. Other gene mutations, including NAV2, ODZ1, FAT2, MTERFD3, B2M, HMCN1, MLL2, TET3 and LYN, were identified in a single case each.
section_name: RESULTS

section_content: Previously, we reported that G17V RHOA mutations were detected by flow cytometry only in CD4-positive (CD4+) cells but not in other cell lineages purified from the skin tumor of a PTCL-NOS patient and the pleural effusion cells of an AITL patient. 8 These preliminary results suggested that the G17V RHOA mutation may specifically exist in mature CD4+ T cells in PTCL-NOS and AITL. To gain further insight into the origin of the G17V RHOA mutation, we examined the mutation in laser-microdissected PD1+ and CD20+ B cells, which were assumed to be enriched and depleted in tumor cells, respectively, in 10 nodal T-cell lymphomas (1 nodal PTCL with TFH phenotype and 9 AITL cases). The G17V RHOA mutation was detected only in the PD1+ cells but not in the CD20+ cells in all 10 cases (Figure 2 ). The allele frequencies of the G17V RHOA mutations in the dissected PD1+ cells were substantially higher than those in the matched whole tumor samples in 7 of the 10 cases. The efficiency of mutation allele enrichment was not substantial in three cases (PTCL63, PTCL78 and PTCL127). In these cases, PD1+ cell selection was not sucessful enough to purify the tumor cells because of the presumed abundance of PD1+ nontumor cells or the very high tumor cell content before the selection. Additionally, using flow cytometry, we found an AITL case showing that the G17V RHOA mutation existed in PD1+CD4+ cells sorted from bone marrow mononuclear cells (Supplementary Figure S2 ). This finding strongly supports our hypothesis that the acquisition of the G17V RHOA mutation is a specific event in TFH cells. 6] [17] We examined the distribution of TET2, IDH2 and DNMT3A mutations in PD1+ and CD20+ cells. Twenty of the 26 TET2 mutations were identified in both the PD1+ and the CD20+ cells (Supplementary Table S7 ), and 15 of the 16 TET2mutated samples had at least one mutation in both the PD1+ and the CD20+ cells (Figure 3 ). Concomitantly, DNMT3A mutations were identified in both the PD1+ and CD20+ cells in four of the seven DNMT3A-mutated samples (Figure 3, Supplementary Table S8 ). In myeloid malignancies, TET2 and IDH2 mutations are known to be mutually exclusive. 14, 29 However, we and others reported that IDH2 mutations often coexist with TET2 mutations in PTCL. 8, 10, 30 DH2 mutations were identified in PD1+ cells but not in CD20+ cells in all 4 TET2-and IDH2-comutated samples (PTCL8, PTCL61, PTCL63 and PTCL70) (Figure 3 ). Each of these samples had at least one TET2 mutation in both the PD1+ and CD20+ cells and the G17V RHOA mutation only in the PD1+ cells. That is, TET2, IDH2 and G17V RHOA mutations coexisted in the PD1+ cells in these cases. In addition, we also found the coexistence of IDH2, TET2 and G17V RHOA mutations in PD1+CD4+ cells sorted from the bone marrow mononuclear cells of an AITL patient (Supplementary Figure S2 ).
section_name: Specific existence of G17V RHOA mutations in tumor cell-enriched cells of nodal T-cell lymphomas

section_content: To clarify the cellular origin of newly identified gene mutations, we also checked the distribution of these mutations in PD1+ and CD20+ cells (Table 2 ). We identified B2M, COL19A1, HMCN1, MTERFD3 and TET3 mutations only in PD1+ cells but not in CD20+ cells. COL19A1, LYN, NAV2 and NOTCH2NL mutations were identified in both the PD1+ and CD20+ cells (Figure 4 ). 

Interestingly, three NOTCH1 and one FAT2, MLL2 and ODZ1 mutations each were found only in the CD20+ but not in the PD1+ cells in four samples (PTCL 63, PTCL70, PTCL78 and PTCL128) (Figure 5 ). Especially, all three NOTCH1 mutations identified by targeted sequencing were identified only in the CD20+ cells with high allele frequencies. The NOTCH1 gene encodes a transmembrane protein. One of the NOTCH1 mutations was a frameshift mutation residing in the PEST domain of the Notch1 protein. This would be an activating mutation, because deletion of the PEST domain enhances Notch signaling after ligand binding. 19 The other two mutations were located in one of the epidermal growth factor-like and in the ankyrin repeat domains (Supplementary Figure S1 ). One of the NOTCH1-mutated samples simultaneously had two TET2 mutations and G17V RHOA mutation (PTCL 63, Supplementary Table S9 ). In this case, both TET2 mutations were detected in both the PD1+ and CD20+ cells, while the G17V RHOA mutation was confined to the PD1+ cells. We used the multiplex PCR method 26 to also check the clonality of immunoglobulin genes in the samples with B-cell-specific mutations. Interestingly, only one sample showed monoclonal rearrangement while the others showed oligoclonal rearrangement (Table 3 ).
section_name: B-cell-specific mutations in nodal T-cell lymphomas

section_content: By determining the distribution of the mutations, we elucidated the clonal architecture of nodal T-cell lymphomas. RHOA mutations were identified only in PD1+ cells in 100% cases, while TET2 and DNMT3A mutations were identified in both the PD1+ cells and CD20+, tumor-cell-depleted cells in the majority of cases. In addition, IDH2 mutations were actually found only in the PD1+ cells and coexisted with TET2 mutations. These data suggest that, in nodal T-cell lymphoma development, multistep tumorigenesis may progress in association with the differentiation of blood cells/ lymphocytes. Surprisingly, some of the mutations resided in a B-cell-specific manner. 

Recent genetic studies have revealed that, in several hematological cancers, several gene mutations existed in preleukemic hematopoietic cells as well as in tumor cells; 31 examples are TET2 and/or DNMT3A mutations in acute myeloid leukemia [32] [33] [34] and NOTCH1 and SF3B1 mutations in chronic lymphocytic leukemia. 35, 36 8] [39] [40] Similarly, our data indicated that in nodal T-cell lymphomas, premalignant cells having TET2 and/or DNMT3A mutations may differentiate not only into T-lineage tumor cells but also into B cells. In contrast, the G17V RHOA mutations specifically existed in the T cells of nodal T-cell lymphomas in all 13 cases (11 cases have been described in this paper, while 2 were previously described elsewhere 8 ), indicating that the G17V RHOA mutation is the event after the B-and T-cell specification. This could happen right after the T/B specification, after differentiation into TFH cells or even after malignant transformation establishing a subclone. One possibility is that the G17V RHOA mutation occurs in TET2-mutated premalignant cells and facilitates the selective differentiation of TET2-mutated premalignant cells into tumor cells with the TFH phenotype. This needs to be proven in the future. 

IDH2 mutations were also specifically identified in the tumorcell-enriched cells, suggesting that IDH2 mutations are also tumorcell-specific events in AITL, although the number of samples was not large enough to allow a definite conclusion. We have previously showed that the IHD2-mutated cases were almost a subcohort of G17V RHOA-mutated cases. 8 This result could be interpreted that the acquisition of IDH2 mutations may be the event occurring after the acquisition of RHOA mutation and thus the IDH2-mutated cells may, at least in some cases such as PTCL70, constitute a subclone in the RHOA-mutated clone. TET2-and IDH2comutated AITL samples were reported to have more extensive histone modification profiles than those with TET2 mutations without an IDH2 mutation, while the difference in genome-wide cytosine methylation profiles between these samples was only moderate. 30 ur data showed that B cells that have infiltrated AITL tissues also have gene mutations: the multilineal mutations represented by those in TET2 and DNMT3A, and B-cell-specific mutations represented by those in NOTCH1 and other genes. ] [46] [47] Some lymphoma cells are infected by Epstein-Barr virus. In such cases, Epstein-Barr virus is proposed to contribute to the transformation of B cells. 45, 46 This hypothesis, however, needs to be re-evaluated because Epstein-Barr virus was not detected in a substantial proportion of B-cell malignancies accompanying AITL. 48 TET2 mutations are found in diffuse large B-cell lymphomas. 49 Tet2-deficient mice show the expansion of both B-and T-cell lineages in addition to prominent myeloproliferation. 15 ombinational loss of Tet1 and Tet2 provokes B-cell malignancies in mice. 50 Activating NOTCH1 mutations were reported in diffuse large B-cell lymphomas, 51 chronic lymphocytic leukemia, 52 mantle cell lymphoma 22 and follicular lymphoma. 53 In our cohort, all three NOTCH1 mutations were defined only in B cells with very high allele frequencies (Figure 5, Supplementary Table S9 ) and two of the three samples showed oligoclonality of B cells (Table 3 ). This implies that the origin of NOTCH1 mutation is earlier than the acquisition of hypermutation of the CDR3 region in the immunoglobulin gene. Anyway, acquisition of these mutations in B-cell lineage may account for the frequent occurrence of B-cell lymphomas in AITL. Moreover, our data alert us to the need for careful interpretation of the mutational profiles of PTCLs because some of the mutations may not exist in tumor cells. 

In conclusion, our findings illustrate the concept of multistep and multilineal tumorigenesis in nodal T-cell lymphomas (Supplementary Figure S3 ). Understanding the pathogenesis will lead us to better management of nodal T-cell lymphomas in future.
section_name: DISCUSSION

section_name: ACKNOWLEDGEMENTS

section_content: We thank Tamaki Takahashi and Yukari Sakashita for their technical assistances. We also thank Ms Flaminia Miyamasu, associate professor of English for Medical Purposes, Medical English Communications Center, University of Tsukuba, for her editorial assistance. This work was supported by Grants-in-Aid for Scientific Research (KAKENHI: 25461407 to MS-Y; and 25112703, 15H01504 and 16H02660 to SC) from the Ministry of Education, Culture, Sports, and Science of Japan. This work was also supported by The Project for Cancer Research and Therapeutic Evolution (P-CREATE) from The Japan Agency for Medical Research and Development, AMED and by research grants from the Daiichi Sankyo Foundation of Life Science, the Naito Foundation, the Cell Science Research Foundation, KANAE Foundation for the Promotion of Medical Science and Shiseido Female Researcher Science Grant to MS-Y; and The Uehara Memorial Foundation, Leukemia Research Fund, Takeda Science Foundation, and Kobayashi Foundation for Cancer Research to SC.
section_name: ACKNOWLEDGEMENTS

section_name: CONFLICT OF INTEREST

section_content: The authors declare no conflict of interest.
section_name: CONFLICT OF INTEREST

section_content: TBN, MS-Y and YA performed experiments and data analysis; KY, YS, KC, HT, SM and SO performed the targeted sequencing; NN, KT, HM, MN and KO were responsible for pathological diagnosis; DM, JK, TM, NN, KT, HM and MN established the LMD procedure; KI and KO collected specimens and were involved in planning the project; TBN, MS-Y and SC generated figures and tables and wrote the manuscript; SC lad the entire project; and all authors participated in discussions and interpretation of the data and results. 

Supplementary Information accompanies this paper on Blood Cancer Journal website (http://www. nature. com/bcj)
section_name: AUTHOR CONTRIBUTIONS