Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

TP53 alterations in primary and secondary Sézary syndrome: A diagnostic tool for the assessment of malignancy in patients with erythroderma

  • Audrey Gros,

    Affiliations INSERM U1053, Bordeaux Research in Translational Oncology University Bordeaux, Bordeaux, France, Tumor Bank and Tumor Biology Laboratory, CHU de Bordeaux, Pessac, France

  • Elodie Laharanne ,

    Contributed equally to this work with: Elodie Laharanne, Marie Vergier

    Affiliations INSERM U1053, Bordeaux Research in Translational Oncology University Bordeaux, Bordeaux, France, Tumor Bank and Tumor Biology Laboratory, CHU de Bordeaux, Pessac, France

  • Marie Vergier ,

    Contributed equally to this work with: Elodie Laharanne, Marie Vergier

    Affiliation Tumor Bank and Tumor Biology Laboratory, CHU de Bordeaux, Pessac, France

  • Martina Prochazkova-Carlotti,

    Affiliation INSERM U1053, Bordeaux Research in Translational Oncology University Bordeaux, Bordeaux, France

  • Anne Pham-Ledard,

    Affiliations INSERM U1053, Bordeaux Research in Translational Oncology University Bordeaux, Bordeaux, France, Dermatology Department, CHU de Bordeaux, Bordeaux, France

  • Thomas Bandres,

    Affiliation Tumor Bank and Tumor Biology Laboratory, CHU de Bordeaux, Pessac, France

  • Sandrine Poglio,

    Affiliation INSERM U1053, Bordeaux Research in Translational Oncology University Bordeaux, Bordeaux, France

  • Sabine Berhouet,

    Affiliation Tumor Bank and Tumor Biology Laboratory, CHU de Bordeaux, Pessac, France

  • Béatrice Vergier,

    Affiliations INSERM U1053, Bordeaux Research in Translational Oncology University Bordeaux, Bordeaux, France, Pathology Department, CHU de Bordeaux, Pessac, France

  • Jean-Philippe Vial,

    Affiliation Hematology Laboratory, CHU de Bordeaux, Pessac, France

  • Edith Chevret,

    Affiliation INSERM U1053, Bordeaux Research in Translational Oncology University Bordeaux, Bordeaux, France

  • Marie Beylot-Barry,

    Affiliations INSERM U1053, Bordeaux Research in Translational Oncology University Bordeaux, Bordeaux, France, Dermatology Department, CHU de Bordeaux, Bordeaux, France

  • Jean-Philippe Merlio

    jp.merlio@u-bordeaux.fr

    Affiliations INSERM U1053, Bordeaux Research in Translational Oncology University Bordeaux, Bordeaux, France, Tumor Bank and Tumor Biology Laboratory, CHU de Bordeaux, Pessac, France

Abstract

Recent massive parallel sequencing data have evidenced the genetic diversity and complexity of Sézary syndrome mutational landscape with TP53 alterations being the most prevalent genetic abnormality. We analyzed a cohort of 35 patients with SS and a control group of 8 patients with chronic inflammatory dermatoses. TP53 status was analyzed at different clinical stages especially in 9 patients with a past-history of mycosis fungoides (MF), coined secondary SS. TP53 mutations were only detected in 10 patients with either primary or secondary SS (29%) corresponding to point mutations, small insertions and deletions which were unique in each case. Interestingly, TP53 mutations were both detected in sequential unselected blood mononuclear cells and in skin specimens. Cytogenetic analysis of blood specimens of 32 patients with SS showed a TP53 deletion in 27 cases (84%). Altogether 29 out of 35 cases exhibited TP53 mutation and/or deletion (83%). No difference in prognosis was observed according to TP53 status while patients with secondary SS had a worse prognosis than patients with primary SS. Interestingly, patients with TP53 alterations displayed a younger age and the presence of TP53 alteration at initial diagnosis stage supports a pivotal oncogenic role for TP53 mutation in SS as well as in erythrodermic MF making TP53 assessment an ancillary method for the diagnosis of patients with erythroderma as patients with inflammatory dermatoses did not display TP53 alteration.

Introduction

Epidermotropic cutaneous T-cell lymphomas, i.e. mycosis fungoides (MF) and Sézary syndrome (SS), account for approximately 65% of cutaneous T-cell lymphomas [1, 2]. Early stages MF (Ia-IIa) are characterized by limited and fixed skin patches or plaques and a low proportion (10 to 20%) of neoplastic T-cells in the upper dermis and epidermis. At advanced stages (IIb-IVb), corresponding to tumor or erythroderma, tumor cells become more abundant and epidermotropism may be lost [2, 3].

Conversely, Sézary syndrome occurs de novo at an advanced stage (IV) and involves more than 90% of the skin with erythroderma or T3, generalized lymphadenopathy and neoplastic T-cells, coined Sézary cells, in skin, lymph nodes and peripheral blood [2, 3]. While most patients present with de novo SS or primary SS, patients meeting the criteria for SS but with past history of MF are also called as secondary SS [2].

The International Society for Cutaneous Lymphoma (ISCL) and the European Organization of Research and Treatment of Cancer (EORTC) have provided consensus staging criteria to evaluate the blood burden in CTCL which are especially useful for the diagnosis of patients with erythroderma [3]. B0 stage or absence of blood involvement correspond to 5% or less Sézary cells; B1 or low blood burden is defined as more than 5% of Sézary cells but either less than 1000 Sézary cells /μL or absence of a clonal rearrangement of the T-cell receptor genes or both; B2 or leukemic phase is defined by clonal rearrangement of the TCR in the blood and either more than 1000 Sézary cells/μL or one of two of the following phenotypic criteria (1) expanded CD4+ or CD3+ cells with CD4/CD8 ratio of 10 or more (2) expanded CD4+ cells with abnormal phenotype including loss of CD7 (>40%) or CD26 (>30%) [3]. Recently, these cutoff levels were reevaluated in a multicenter study at 40% for CD7 loss and 80% for CD26 loss achieving a 83% sensitivity and 100% specificity to distinguish patients with SS from patients with erythrodermic inflammatory dermatoses [4]. Other phenotypic biomarkers such as CD158k expression in more than 5% of CD4+ cells [5] had only 33% sensitivity and 95% specificity [4]. Genetic biomarkers such as MYC gain at 8q24.21 (40%) and MNT loss at 17p13 (30%) were found useful along with as the demonstration of upregulated expression of STAT4, TWIST1 and DNM3 or PLS3, especially for cases with erythroderma and lack of abnormal phenotype [4].

The presence of an identical circulating and cutaneous T-cell clone was found to be associated with the diagnosis of CTCL in most patients having erythroderma [6]. The presence of a clonal T-cell population in the peripheral blood but not in the skin has per se limited diagnostic or prognostic value as it can be observed in reactive conditions, especially in the elderly [68]. Differentiating SS from clinical mimickers such as drug induced erythroderma or erythrodermic inflammatory dermatoses is crucial to formulate an effective diagnosis and appropriate therapy (for a commentary see [9])

In patients with MF or SS, we and others observed that the presence of an identical clone in skin and blood specimens is an independent and adverse prognostic factor whatever the histological subtype and disease stage [1012]. However, there are general difficulties of statistical interpretations of biological or genetic parameters in patients with CTCL, especially when MF and SS cases are mixed together as a single entity.

Therefore, the demonstration of additional phenotypic and/or genetic alterations is necessary to assess blood involvement in CTCL and the malignancy of circulating monoclonal T-cells [13, 14].

Cytogenetic analysis using array comparative genomic hybridization (aCGH) have identified a specific cytogenetic profile of SS cells [1517]. Indeed, 17p loss involving the TP53 gene was observed by several groups in patients with SS but appeared rare in patients with MF except in cases with large cell transformation [15, 1720]. These studies detected genetic copy number variation (CNV) with losses of genes such as CDKN2A, E2A, A20 or TP53 being more frequent than gains such as those of MYC [15, 17, 21].

Recent next generation sequencing (NGS) analysis of CD4+ selected blood cells in patients with SS have evidenced a variability of genetic aberrations resulting in a unique combination per patient of chromosomal imbalances, fusion genes and point mutations [2225, 26, 27]. Interestingly, TP53 alterations (mutation and/or copy number loss) were found as the most common genetic abnormality in patients with SS in several studies [2426]. However, in the largest series of 101 SS cases studied so far [28], TP53 mutation or deletion were detected in only 15% and 6% of patients, respectively. Indeed, the prevalence of TP53 mutations was found highly variable ranging from 15% to 43% of cases by most exome sequencing studies (for a review see [22]). A first pivotal focused study identified a heterozygous loss of TP53 was detected in 8 out of 9 patients with SS but TP53 deleterious mutations were only detected in 3 out of 8 patients using Sanger sequencing [29]. Previous studies did not determine if TP53 alterations were present at onset of the disease or detected in patients under treatment or during progression. The aim of the present study was to evaluate the incidence of TP53 mutations in patients with primary and secondary SS at onset and during the course of the disease. Therefore, NGS analysis of TP53 status was performed on unselected peripheral blood mononuclear cells and in skin specimens in the same patients. Then, we analyzed whether TP53 mutations were associated with specific clinical features such as age, past-history of MF or evolutive features such as tumor formation and/or large cell transformation. We also evaluated TP53 status along with several biological parameters such as the amount of Sézary cells (> or < 1000 cells/mm3), the presence of TP53 genetic loss as determined by FISH analysis, and the length of telomeres.

Materials and methods

Patients and material

The cohort contained peripheral blood samples of 35 patients with SS at B2 stage according to the ISCL/EORTC criteria. According to the French Public Health and Bioethical Law, the present study was considered by our research direction lawyer as a non-interventional study without the need for ethics committee approval (Article L1121-1 and Article R1121-3). Moreover, informed consent was given to all patients with written consent for the use of the clinical data and research on their biological material. All data were collected anonymously.

In the 35 patients with SS, peripheral blood mononuclear cells (PBMC) were isolated using Pancoll (PAN-Biotech, Aidenbach, Germany) according manufacturer’s recommendations at SS diagnosis stage prior to treatment and for 28 patients during follow-up with a total of 63 analyzed blood samples. In patients with TP53 mutation within blood samples, we also analyzed the initial skin specimen either at SS diagnosis time in 7 patients with primary SS and at MF stage in the 3 patients with secondary SS. Genomic DNA was extracted using standard phenol/chloroform protocol from peripheral blood mononuclear cells or frozen skin specimens.

A control group of 8 patients with idiopathic chronic inflammatory erythroderma was referred to the Dermatology department. There was neither evidence of drug induced erythroderma nor history of atopic dermatitis or psoriasis. Skin biopsy showed non-specific features and both PBMC count with phenotypic analysis and clonality analysis were performed for differential diagnosis with SS. In five patients, a paired analysis of skin and blood specimens was performed while in three patients only a skin specimen was evaluated for TP53 status.

Flow cytometry analyzes and sorting

In peripheral blood monuclear cells, TCR Vβ was determined using flow cytometry with the IOTest Beta Mark TCR Repertoire Kit (Beckman Coulter, Marseille, France) designed to identify 24 distinct TCR Vβ families according to manufacturer’s recommendations. Data acquisition was performed using a BD FACSCANTO II flow cytometer and analyzed using BD FACSDIVA software (BD Biosciences, le Pont de Claix, France). Tumor and non-tumor populations were sorted based on TCR Vβ22, CD4 and CD8 expression using BD FACSAriaTM II cell sorter (BD Biosciences).

Targeted deep sequencing

Genomic mutation of TP53 was studied by Next Generation Sequencing (NGS) using the Ion Ampliseq TP53 panel (ThermoFisher Scientific, Life Technologies, Les Ullis France). The design covers 100% of the coding sequence of TP53 including mutations at canonical splice sites and generates 24 amplicons on a length of 1 280 bases. Libraries were amplified by emulsion PCR and enriched using automatic system IonChef (ThermoFisher Scientific). Ion sphere particles were then sequenced with the Ion torrent personal genome machine (ThermoFisher Scientific) on 316v2 or 318v2 Chips (ThermoFisher Scientific) with high mean coverage > 3000X. Torrent Suite version 5.0 software (ThermoFisher Scientific) was used to perform data analysis. Reads were mapped to the human hg19 reference genome. Data processing, alignment and mutation calling were performed using the Torrent Suite. The Variant Caller detected point mutations with a variant frequency ≥2% for Single Nucleotide Variation (SNV) and ≥5% for short insertion/deletion (INDEL). VCF files generated by Variant Caller were annotated by ANNOVAR [30].

In the 10 patients with TP53 mutation, BAM of wild-type samples for TP53 were also checked using Alamut Software (Interactive Biosoftware, Rouen, France).

TCRG gene rearrangement analysis

TCRG gene rearrangement was studied by using a GC-clamp multiplex PCR denaturing gradient gel electrophoresis (DGGE) as originally described and used in further reports [6, 12, 31]. The procedure allows a direct clonality assessment and comparison of clonal rearrangement detected in blood and skin samples as migration of clonal bands depends on size and sequence-specific denaturation [6, 12]. For comparison of tumoral and non tumoral flow sorted cells, we performed TCRG gene rearrangement analysis using a standardized Biomed-2 protocol with commercially available reagents (Invivoscribe, San Diego, CA, USA) according to the EuroClonality/BIOMED-2 guidelines [32]. GeneScanning of PCR products was performed on a Applied 3100 genetic analyzer (Life Sciences, Les Ullis, France).

FISH analysis

TP53 FISH assay was performed using Metasystems FISH XL P53 Deletion Probe (D-5017-025-RG, Altussheim, Germany) according to the manufacturer’s protocol. Unselected fresh blood cells isolated by gradient separation or after short-term culture were fixed in Methanol:Acetic Acid (3:1) were spread on SuperFrost Plus positively charged slides, pretreated in 2x SSC at 37°C, and passed through graded ethanol and air-dried, as described [33, 34]. Formalin-fixed sections were also used and pretreated, as described [35]. Slides were co-denaturated with probe at 75°C for 2 minutes and hybridized at 37°C overnight using an automated HyBrite co-denaturation oven (Abbott Molecular, Rungis, France). Slides were then placed in washing buffer 1 (2x SSC/0.1% NP40) at room temperature for 2 min to remove the coverslips, immersed in 73°C washing buffer 2 (0.4x SSC/0.3% NP40) for 2 min, then placed in washing buffer 1 at room temperature for 1 min, passed through graded ethanol, air-dried and mounted with Vectashield mounting medium with DAPI (Vector laboratories, CliniSciences, Nanterre, France). Slides were analyzed with an epifluorescence microscope equipped with appropriate single band-pass filters (Abbott Molecular). FISH analysis was performed by 2 trained lecturers. The slide was considered as interpretable if the following requirements were fulfilled: (1) bright and distinct signals (2) low background relatively free of fluorescent particles or haziness.

Quantitative fluorescence in situ hybridization

Quantitative fluorescence in situ hybridization (Q-FISH) was performed as previously mentioned [34]. Briefly, cell nuclei were prepared and spread onto microscope slides. The hybridization was performed according to manufacturer’s instructions (DakoCytomation, Trappes, France) with minor modifications. Thus, after pre-treatment steps nuclei were denatured in formamide 70% in 2xSSC at 70°C, rinsed in cold ethanol series and air dried. Peptide nucleic acid (PNA) probe for telomere conjugated with Cyanine 3 (Cy3) dye, was denatured at 80°C for 5 min. Nuclei were hybridized at 30°C for 2 hours into a ThermoBrite system (Leica Microsystemes, Nanterre, France). Then, slides were rinsed, washed, air dried and mounted with Vectashield (Vector Laboratories) containing 4′,6-diamidino-2-phenylindole dihydrochloride (DAPI).

Statistical analysis

Data were analyzed according to the number of dependent variables and to the nature of dependent and independent variables. In order to test the normality of our data the Shapiro-Wilk test was used. For non-parametric statistical analysis the Wilcoxon-Mann Whitney or the Fisher’s exact test were used. This test is suitable for sample number below 30. To remove outlier values, Dixon’s test was applied. Survival was evaluated by Cox's proportional hazards analysis. The significance level was taken at P ≤ 0.05. The statistical tests were performed using the MedCalc software (Ostend, Belgium).

Immunohistochemical for PD1, P53 and MDM2 expression

Immunohistochemical analysis was performed on fixed skin sections. The following primary antibodies and dilutions were used: Anti-p53 diluted 1:100 (Clone DO-7, DakoCytomation, Denmark), anti-PD1 diluted 1:50 (Clone MRQ-22 Zytomed, Diagomics, Blagnac, France) and anti-MDM2 diluted 1:50 (Clone IF2, ThermoFisher Scientific). Sections were incubated in 90°C 30mn/ER1 for P53, 20mn/ER2 for PD1 and 20mn/ER1 for MDM2. The analysis was performed using Bond-Max (Leica Microsystemes) according to the manufacturers’ recommendations. Both percentage of positive tumor cells and staining intensity (1+ to 3+) were recorded according to guidelines in an attempt to determine an H-score [36]. The predominant pattern (percentage x intensity) was finally recorded.

Results

Description of the cohort

Characteristics of the patients are summarized in Table 1 and S1 Table.

thumbnail
Table 1. Clinical and biological features of patients with primary Sézary Syndrome (SS) or with past-history of Mycosis Fungoides (MF) at secondary SS stage.

https://doi.org/10.1371/journal.pone.0173171.t001

After review of initial and current stage, 26 patients presented with primary SS and 9 patients with secondary SS had a past-history of MF, mostly erythrodermic (n = 7). Patients with primary SS had a median age of 75.5 years at diagnosis and gender ratio (M/F) was 1. For patients with a past-history of MF or secondary SS, median age at MF diagnosis was 62 years and the gender ratio was 0.5. Mean and median follow-up of patients with primary SS were 52.0 and 45.5 months, respectively. Mean and median follow-up of patients with secondary SS were 70.0 and 50.0 months, respectively. In all patients, one blood sample at SS diagnosis was analyzed for TP53 mutation as well as for TCRG gene rearrangement showing a monoclonal rearrangement (S1 Table). In the control group, median age was 74 and the gender ratio was 3.

Detection of somatic TP53 mutations in peripheral blood mononuclear cells of patients with SS

We performed targeted deep sequencing of all coding and intron-exon flanking sequences of the TP53 gene and studied 63 blood samples of 35 patients with SS using the Personal Genome Machine (PGM) and a commercially available NGS TP53 panel kit. Somatic TP53 mutations were identified in the blood of 10 out of 35 patients (29%) corresponding to 4 nonsense SNVs, 3 missense SNVs, 2 in-frame INDELs, and 1 frameshift INDEL. Interestingly, each case exhibited a single and unique mutation (Fig 1).

thumbnail
Fig 1. Schematic representation of somatic mutations identified in TP53 by targeted deep sequencing (n = 35).

Mutation sites were marked and amino acid changes were indicated. Colors and shapes indicated the kind of mutation and the presence in the COSMIC database.

https://doi.org/10.1371/journal.pone.0173171.g001

Among the 10 TP53 mutations, two of them were present within the domain encoding for tetramerization property, and seven in the DNA-binding domain. Prediction analysis of Polyphen-2 (Polymorphism Phenotyping v2, http://genetics.bwh.harvard.edu/pph2/) showed a score greater than 0,96 for the 3 missense mutations. The three variants with INDELs were found absent in the COSMIC v77 database. The 2 in-frame deletions of 5 and 8 Amino Acids involved the DNA-binding domain and the frameshift mutation resulted in a stop codon at the beginning of this domain.

TP53 mutation is present in clonal T-cells

For patient P3, we first analyzed TCR Vβ expression on PBMC using the panel of 24 monoclonal antibodies to TCR Vβ families. Such analysis identified the expansion of a restricted TCR Vβ22+ CD4+CD8- population and these cells were isolated from other PBMC using flow cytometry cell sorting (Figure A in S1 Fig). After DNA extraction, we performed both clonality analysis of TCRG gene rearrangement and targeted deep sequencing of the TP53 gene (Figure B in S1 Fig). The same TP53 mutation as in the original unselected PBMC sample was detected in the TCR Vβ22+ CD4+CD8- sorted population with a 99% allelic frequency while a wild-type status was observed in the other sorted PBMC cells (Figure C in S1 Fig). This was correlated with the detection of a monoclonal TCRG rearrangement in the TCR Vβ22+ CD4+CD8- sorted cells and a polyclonal profile in the others (Figures A, B and C in S1 Fig).

TP53 mutations can be reliably detected in sequential blood specimens and in the skin of patients with SS

We analyzed 28 sequential blood samples from patients with TP53 mutation in the blood. We also analyzed 7 skin samples of patients with primary SS and TP53 mutation as well as 3 skin samples in patients with secondary SS. All BAM files of wild-type samples of these patients were checked using Alamut Software (Interactive Biosoftware, Rouen, France) in order to detect very low frequency mutations.

The same and unique mutation was found both in skin and blood samples in the 10 cases with TP53 mutation. In these patients, sequential blood samples collected after treatment showed the persistence of the same TP53 mutated clone. In six of these 10 patients and also in nine patients without mutation, sequential blood samples were collected after extracorporeal photophoresis and no additional TP53 mutation was observed. Altogether, the incidence of TP53 mutation was equivalent in patients who ever received UVA irradiation or not.

Several whole exome sequencing analysis reported a high rate of C>T transition in SS point mutations suggesting a causal role for UV exposure [25, 37]. In our cohort, 6 out of 7 point mutations corresponded to a C>T transition. However, such C>T transitions were not observed at NpCpC sites which correspond to canonical motifs modified after UVB exposure but occurred at NpCpG sites in 3 out of these 6 cases which has been attributed to a senescence-related process [38].

Polymorphism P72R or P72P of the TP53 gene

A specific TP53 polymorphism Pro72-to-Arg (P72R) was suggested to be associated with cancer susceptibility or aggressiveness [39]. Indeed, this SNP was found homozygous or heterozygous in 5 out of 6 SS patients with wild-type TP53 status [29]. Alternatively, Mc Girt et al observed a significant increase in the frequency of Pro72Pro (P72P) in a cohort of patients with MF and suggesting that this polymorphism could contribute to MF susceptibility in patients at a younger age than those with SS [40]. We therefore evaluated the presence of these single nucleotide polymorphisms (SNP) both at the diagnosis and during follow-up in SS samples. In the 35 patients, 22 were homozygous Arg72Arg (63%), 12 were heterozygous Pro72Arg (34%), and one was homozygous Pro72Pro (3%). There was no difference between our patients and the general population for the allelic frequency for this SNP: Arg72Arg (56%), Pro72Arg (37%), Pro72Pro (7%). When focusing on the 9 patients with MF and then secondary SS, we observed a slight increase of the homozygous Arg72Arg (78%) but such difference may result from the small size of this subgroup. In the control group of 8 inflammatory dermatoses, the distribution was similar with Arg72Arg (n = 6, 67%) and Pro72Arg (n = 2, 33%).

TP53 Copy Number Variation (CNV)

Blood samples of 32 patients with SS were analyzed by FISH for TP53 locus. No material was available in three patients including two patients with TP53 mutation. FISH analysis showed a recurrent deletion of TP53 in 27 out of 32 patients (84%) (Fig 2). In 8 of 10 patients with TP53 mutation, a deletion of the TP53 locus of at least one allele was observed. Two of the 8 cases exhibited 3 or 4 centromeric signals for chromosome 17 with loss of one or two TP53 alleles, respectively. Six of these 8 patients showed a diploid pattern with loss of one TP53 allele associated with mutation of the remaining allele suggesting a double knock-out.

thumbnail
Fig 2. Clinical data and somatic alterations of TP53 gene in Sézary syndrome.

Age at diagnosis (years), Clinical status means survival time in months after diagnosis of Sézary syndrome until the death or last clinical status, ⱡ data determined by fluorescence in situ hybridization, ¥ data determined by quantitative fluorescence in situ hybridization.

https://doi.org/10.1371/journal.pone.0173171.g002

Altogether, an imbalance and/or mutation of TP53 were detected in 27 out of 32 patients who were analyzed both by NGS and FISH analyses (84%).

Clinical and biological features associated with TP53 alterations

Patients with TP53 alterations (mean age 68) were younger than patients without TP53 alteration (mean age 82) (p = 0.0075, Mann-Whitney test). We also analyzed survival according to TP53 status but found no statistical difference. Survival was however shorter in patients with previous history of MF than in patients with primary SS whatever the TP53 status (28 months versus 37 months, respectively, p = 0.0038) (S2 Fig).

When studying the number of circulating Sézary cells (less or more than 1000/μL, phenotypic criterion used for B2 staging (loss of CD7 and/or CD26), large-cell histological transformation, no specific feature was observed according to TP53 status (Fig 2). However, large cell transformation was associated with worse prognosis (p = 0.0075) (S2 Fig).

When studying telomere length by quantitative FISH and q-PCR as reported previously [34], TP53 alterations were more frequently observed in cases with shorter telomeres than in cases with longer telomeres (p<0.05, Fisher exact test) (Fig 2 and S3 Fig).

P53 expression

We analyzed p53 and mdm-2 expression in tumoral cells in 17 skin samples of patients with SS by visual comparison with PD1 staining of SS cells to minimize counting of non-tumoral cells (data not shown). Indeed, p53 protein level was found increased by abnormal stabilization of the mutated form and sometimes associated with mdm-2 protein overexpression in CTCL cases [41, 42]. The cutoff used for IHC positivity was an immunoreactivity of more than 30% of tumoral cells but we also evaluated staining intensity as reported for individual cases (Fig 2). In patients with TP53 mutation, p53 expression was observed in 89% of cases (8/9). Five of these cases presented a double TP53 alteration (one mutated allele and one deleted allele) but were stained for p53 in 30–80% of tumour cells. Only one sample was negative corresponding to the patient with TP53 frameshift mutation and deletion. We also observed p53 immunostaining in 44% of the patients (4/9) with TP53 wild-type status. To explore whether p53 stabilization could be associated with the expression of the key negative regulator mdm-2, we analyzed immunostaining for mdm-2 in 14 skin biopsies. None of them was found positive (data not shown).

TP53 alterations are not detected in the skin and blood of 8 patients with chronic erythrodermic inflammatory dermatoses

Five pairs of blood and skin samples and three skin samples of 8 patients with reactive erythroderma were studied by Next Generation Sequencing (NGS) using the Ion Ampliseq TP53 panel. No TP53 mutation was detected.

Clonality analysis of TCRG gene rearrangements showed a polyclonal profile in samples of 4 patients (skin analysis in one patient, parallel blood and skin analysis in 3 patients). A monoclonal TCRG gene rearrangement was detected in samples of 4 patients (only at the skin level in 3 patients, both at the skin and blood level in one patient but with a different monoclonal band).

Skin sections of 7 of these 8 patients with reactive erythroderma were subjected to FISH analysis along with skin sections of 2 patients with SS and TP53 monoallelic loss at the blood level. Interestingly, no TP53 loss was observed in the dermal lymphoid infiltrate of patients with reactive erythroderma (Figure A in S4 Fig). In the two patients with SS, the upper-dermal infiltrate contained more than 20% of cells with an allelic loss pattern corresponding to the one observed at the level of PBMC of the same patient (Figure B in S4 Fig).

Discussion

Next generation sequencing (NGS) techniques provide novel tools to determine the genetic landscape of cancer subtypes and to examine in depth the presence of single gene mutations. Using DNA extracted from unselected PBMC for routine TCR gene clonality analysis, a TP53 mutation was detected in 10 out of 35 of SS cases (29%). Such rate is in the range of those reported in the different whole exome sequencing series using purified CD4+ blood cells (for a review see [22]). A unique mutation per patient was detected underscoring the need of the genetic coverage achieved by NGS. Here, we also demonstrate the presence of TP53 mutations before treatment in 9 of these 10 patients. Targeted deep sequencing (>3000X) also detected the same mutation at the skin level in all patients with SS. While TP53 mutations have not been detected by NGS in patients with classical MF [23, 40], they have been observed in some patients with erythrodermic MF and a blood T-cell clone [41]. The presence of an identical TP53 mutation at the skin and blood level in patients at erythrodermic MF stage evolving into secondary SS establishes a biological link with SS, as suggested by others [25, 27].

In SS, the pattern of TP53 gene mutations was suggestive of a role for UV radiation as point mutations were C>T transition occurring at di-pyrimidine sites [25, 41]. This would support that SS cells may originate from skin resident cells while others suggested that they derive from central memory T-cells with circulating capacities [43]. Although our data are limited to a single gene analysis, we have not observed that C>T transitions occurred at NpCpC sites for the TP53 gene. TP53 driving mutations in such elderly patients with SS may rather result from aging than from UV exposure, as in chronic lymphoid leukemia (CLL) [44]. However, as reported in two patients with a double inactivated TP53 clone [26], we observed a relative younger age for patients with TP53 alterations (mutation and/or deletion) than without.

Strikingly, TP53 mutations have been detected at the skin level not only in B and T–cell cutaneous lymphomas but also in indolent lymphoproliferations such as lymphomatoid papulosis [41]. This strongly suggests that other cells such as keratinocytes may account for a mutated status when skin specimens only are analyzed. It also underscores the reliability of our parallel testing of blood and skin specimens for TP53 status assessment.

At the biological level, the TP53 mutations reported here were found to mainly affect the DNA binding or tetramerization domain of TP53 [45]. By prediction analysis, they were found deleterious with a high probability score. As recently shown in breast cancers, NGS represents an unequivocal tool for the detection of mutations confirmed as mutant by the functional assay [46].

Using a simple FISH analysis on either unselected PBMC or skin sections, we also found the presence of at least one allelic alteration of TP53 in 83% of SS cases. By comparison with patients with inflammatory erythroderma, we provide here evidence that TP53 alterations may be used to assess malignancy of circulating cells in patients with either erythroderma or a history of MF. Since the submission of our manuscript, another group reported the use of a panel of 11 FISH probes for the assessment of leukemic CTCL including SS, Transformed MF, Follicular-MF and MF [14]. Such customized panel was designed according to a recent exome study of selected CD4+ cells and evaluated against criteria established by the same group for leukemic evaluation in CTCL [27, 47]. Interestingly, TP53 loss was a constant finding in their 10 SS cases [14]. Using international criteria to define primary and secondary SS, our study strongly supports that TP53 assessment by FISH may be a robust tool for leukemic assessment of CTCL. This approach is also applicable to skin specimens for the differential diagnosis with inflammatory erythroderma.

The TP53 tumor suppressor gene is involved in maintaining genomic stability and mediating response of cancer cells to therapy. In CLL, TP53 inactivation is associated with complex karyotypes at onset of the disease providing a fitness advantage under therapy [44]. Similarly, most of our SS cases with TP53 alteration displayed complex karyotypes (our personal data) and copy number alterations [17]. No correlation between TP53 status and prognosis or large cell transformation was observed which may depend on other genetic events such as MYC gain or CDKN2A deletion [21, 48]. In SS cell lines such as SeAx, Hut78 and HH, TP53 mutations were previously found associated with functional resistance to nutlin-3, an mdm-2 antagonist [29]. Whether TP53 alteration may impact therapeutic response should be evaluated along with other genomic features in prospective trials.

SS mutational landscape has been reported to be complex with high inter-individual variability, especially in the generation of original fusion genes per patient [25, 26, 28]. TP53 alteration in SS not only represents the most common genetic feature but may have a pivotal role in disease initiation as suggested by the early detection at diagnosis. Moreover, TP53 alterations were found to correlate with telomere shortening in our patients. Such association may allow Sézary cells to escape senescence or checkpoint control after telomere crisis generating diverse chromosomal breaks and imbalances.

Altogether, TP53 alterations were detected by targeted NGS or FISH analysis in 29 out of 35 patients with SS (83%) and this incidence was slightly higher in the 32 cases analyzed by both techniques (84%). Therefore, either FISH alone or a combined analysis with NGS may be a valuable tool for the diagnosis of patients with erythroderma and a circulating monoclonal T-cell population, especially in cases lacking the phenotypic criteria proposed by the ISCL/EORTC [3]. Both phenotypic and genotypic criteria will probably replace the morphological criterion of Sézary cell count which lacks reproducibility [47]. Although less recurrent genetic imbalances such as MYC gain (40%) and MNT loss (30%) together with variable genomic or epigenetic criteria have been proposed to differentiate SS from erythrodermic inflammatory dermatoses [4, 49], our study is original in supporting the introduction of a simple and robust TP53 testing in algorithms for the diagnosis of patients with erythroderma as well as for leukemic assessment in patients with CTCL.

Supporting information

S1 Table. Patient features and histopathological data.

* for patients with large cell transformation.

https://doi.org/10.1371/journal.pone.0173171.s001

(DOCX)

S1 Fig. Flow cytometry, clonality analysis and TP53 assessment demonstrate the presence of TP53 mutation in TCRVβ22+ clonal T-cells.

A: Flow cytometry of patient P3 peripheral blood mononuclear cells (PBMC) (left column) followed by cell sorting of TCRVβ22+ tumor cells (median column) and non-TCRVβ22+ cells (right column). B: TCRG gene rearrangement analysis by the BIOMED-2 gene scanning assay C: Results of targeted deep sequencing of the TP53 gene (The presence of an identical TCRG monoclonal rearrangement in total PBMC and sorted TCRVβ22+ tumor cells correlates with the presence of an identical TP53 c.991C>T nonsense mutation while non-TCRVβ22+ sorted cells exhibit a polyclonal profile and a wild-type TP53 status.

https://doi.org/10.1371/journal.pone.0173171.s002

(TIF)

S2 Fig. Survival analysis according to Cox's proportional hazards analysis between patients with primary or secondary Sézary syndrome (A) and patients with or without large cell transformation (B).

https://doi.org/10.1371/journal.pone.0173171.s003

(TIF)

S3 Fig. Correlation between telomeres length and TP53 status according to Fisher’s exact test.

The length of telomeres was determined by quantitative fluorescence in situ hybridization (FISH). TP53 status was assessed both by FISH and targeted deep sequencing.

https://doi.org/10.1371/journal.pone.0173171.s004

(TIF)

S4 Fig. Fluorescence in Situ Hybridization (FISH) analysis of skin sections of patients with either inflammatory dermatoses or Sézary syndrome.

A:. A balanced and diploid FISH pattern after hybridization with the locus specific TP53 probe (red signals) and the control chromosome 17 probe (green signals) is observed in inflammatory cells of the upper dermis in one patient with inflammatory erythroderma. B: A monoallelic TP53 deletion with loss of one red signal corresponding to the TP53 locus is detected in Sézary cells (arrows) present in the upper dermis in one patient with Sézary syndrome.

https://doi.org/10.1371/journal.pone.0173171.s005

(TIF)

Acknowledgments

This work is part of the MATULILA project funding to Anne Pham-Ledard supported by the University Hospital of Bordeaux, the Ligue contre le Cancer Comité de Gironde grant to Sandrine Poglio and the Association pour la Recherche Contre le Cancer grant to Edith Chevret. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. There was no additional external funding received for this study.

Author Contributions

  1. Conceptualization: AG APL SP MBB JPM.
  2. Data curation: AG EL MPC TB SP SB.
  3. Formal analysis: AG EL MPC.
  4. Funding acquisition: APL SP EC.
  5. Investigation: APL MBB.
  6. Methodology: AG EL MV MPC SP SB JPV BV MBB.
  7. Project administration: JPM.
  8. Resources: BV JPM.
  9. Software: AG TB.
  10. Supervision: JPM.
  11. Validation: AG EL MPC APL SP BV MBB JPM.
  12. Visualization: AG EL MPC SP JPM.
  13. Writing – original draft: AG MPC APL EC MBB JPM.
  14. Writing – review & editing: AG MPC APL EC MBB JPM.

References

  1. 1. Willemze R, Kerl H, Sterry W, Berti E, Cerroni L, Chimenti S, et al. EORTC classification for primary cutaneous lymphomas: a proposal from the Cutaneous Lymphoma Study Group of the European Organization for Research and Treatment of Cancer. Blood. 1997;90(1):354–71. pmid:9207472
  2. 2. Willemze R, Jaffe ES, Burg G, Cerroni L, Berti E, Swerdlow SH, et al. WHO-EORTC classification for cutaneous lymphomas. Blood. 2005 May 15;105(10):3768–85. pmid:15692063
  3. 3. Olsen E, Vonderheid E, Pimpinelli N, Willemze R, Kim Y, Knobler R, et al. Revisions to the staging and classification of mycosis fungoides and Sezary syndrome: a proposal of the International Society for Cutaneous Lymphomas (ISCL) and the cutaneous lymphoma task force of the European Organization of Research and Treatment of Cancer (EORTC). Blood. 2007 Sep 15;110(6):1713–22. pmid:17540844
  4. 4. Boonk SE, Zoutman WH, Marie-Cardine A, van der Fits L, Out-Luiting JJ, Mitchell TJ, et al. Evaluation of immunophenotypic and molecular biomarkers for Sezary syndrome using standard operating procedures: multicenter study of 59 cases. The Journal of investigative dermatology. 2016 Feb 27.
  5. 5. Bagot M, Moretta A, Sivori S, Biassoni R, Cantoni C, Bottino C, et al. CD4(+) cutaneous T-cell lymphoma cells express the p140-killer cell immunoglobulin-like receptor. Blood. 2001 Mar 1;97(5):1388–91. pmid:11222384
  6. 6. Delfau-Larue MH, Laroche L, Wechsler J, Lepage E, Lahet C, Asso-Bonnet M, et al. Diagnostic value of dominant T-cell clones in peripheral blood in 363 patients presenting consecutively with a clinical suspicion of cutaneous lymphoma [In Process Citation]. Blood. 2000;96(9):2987–92. pmid:11049975
  7. 7. Vega F, Luthra R, Medeiros LJ, Dunmire V, Lee SJ, Duvic M, et al. Clonal heterogeneity in mycosis fungoides and its relationship to clinical course. Blood. 2002 Nov 1;100(9):3369–73. pmid:12384439
  8. 8. Muche JM, Lukowsky A, Heim J, Friedrich M, Audring H, Sterry W. Demonstration of frequent occurrence of clonal T cells in the peripheral blood but not in the skin of patients with small plaque parapsoriasis. Blood. 1999 Aug 15;94(4):1409–17. pmid:10438729
  9. 9. Wong HK. STAT Assays with a TWIST: Differentiating Sezary Syndrome from Erythrodermic Inflammatory Dermatitis. The Journal of investigative dermatology. 2016 Jul;136(7):1313–5. pmid:27342033
  10. 10. Fraser-Andrews E, Woolford A, Russel-Jones R, Seed P, Whittaker S. Detection of a peripheral blood T-cell clone is an independant prognostic marker in mycosis fungoides. The Journal of investigative dermatology. 2000;114:117–21. pmid:10620126
  11. 11. Scarisbrick J, Whittaker S, Evans A, Fraser-Andrews E, Child F, Dean A, et al. Prognostic significance of tumor burden in the blood of patients with erythrodermic primary cutaneous T-cell lymphoma. Blood. 2001;97:624–30. pmid:11157477
  12. 12. Beylot-Barry M, Sibaud V, Thiébaut R, Vergier B, Beylot C, Delaunay M, et al. Evidence that an identical skin and blood T-cell clone is an independent prognostic factor in primary cutaneous T-cell lymphomas. The Journal of investigative dermatology. 2001;117:920–6. pmid:11676833
  13. 13. Marcus Muche J, Karenko L, Gellrich S, Karhu R, Kytola S, Kahkonen M, et al. Cellular coincidence of clonal T cell receptor rearrangements and complex clonal chromosomal aberrations-a hallmark of malignancy in cutaneous T cell lymphoma. The Journal of investigative dermatology. 2004 Mar;122(3):574–8. pmid:15086537
  14. 14. Weed J, Gibson J, Lewis J, Carlson K, Foss F, Choi J, et al. FISH panel for leukemic CTCL. The Journal of investigative dermatology. 2016 Nov 8.
  15. 15. Vermeer MH, van Doorn R, Dijkman R, Mao X, Whittaker S, van Voorst Vader PC, et al. Novel and highly recurrent chromosomal alterations in Sezary syndrome. Cancer Res. 2008 Apr 15;68(8):2689–98. pmid:18413736
  16. 16. Karenko L, Kahkonen M, Hyytinen ER, Lindlof M, Ranki A. Notable losses at specific regions of chromosomes 10q and 13q in the Sezary syndrome detected by comparative genomic hybridization. The Journal of investigative dermatology. 1999 Mar;112(3):392–5. pmid:10084322
  17. 17. Laharanne E, Oumouhou N, Bonnet F, Carlotti M, Gentil C, Chevret E, et al. Genome-Wide Analysis of Cutaneous T-Cell Lymphomas Identifies Three Clinically Relevant Classes. The Journal of investigative dermatology. 2010 Feb 4.
  18. 18. Mao X, Lillington D, Scarisbrick JJ, Mitchell T, Czepulkowski B, Russell-Jones R, et al. Molecular cytogenetic analysis of cutaneous T-cell lymphomas: identification of common genetic alterations in Sezary syndrome and mycosis fungoides. Br J Dermatol. 2002 Sep;147(3):464–75. pmid:12207585
  19. 19. Mao X, Lillington DM, Czepulkowski B, Russell-Jones R, Young BD, Whittaker S. Molecular cytogenetic characterization of Sezary syndrome. Genes Chromosomes Cancer. 2003 Mar;36(3):250–60. pmid:12557225
  20. 20. Salgado R, Servitje O, Gallardo F, Vermeer MH, Ortiz-Romero PL, Karpova MB, et al. Oligonucleotide Array-CGH Identifies Genomic Subgroups and Prognostic Markers for Tumor Stage Mycosis Fungoides. The Journal of investigative dermatology. 2009 Sep 17.
  21. 21. Laharanne E, Chevret E, Idrissi Y, Gentil C, Longy M, Ferrer J, et al. CDKN2A-CDKN2B deletion defines an aggressive subset of cutaneous T-cell lymphoma. Mod Pathol. 2010 Apr;23(4):547–58.
  22. 22. Chevret E, Merlio JP. Sezary Syndrome: Translating Genetic Diversity into Personalized Medicine. The Journal of investigative dermatology. 2016 Jul;136(7):1319–24. pmid:27342034
  23. 23. da Silva Almeida AC, Abate F, Khiabanian H, Martinez-Escala E, Guitart J, Tensen CP, et al. The mutational landscape of cutaneous T cell lymphoma and Sezary syndrome. Nature genetics. 2015 Dec;47(12):1465–70. pmid:26551667
  24. 24. Kiel MJ, Sahasrabuddhe AA, Rolland DC, Velusamy T, Chung F, Schaller M, et al. Genomic analyses reveal recurrent mutations in epigenetic modifiers and the JAK-STAT pathway in Sezary syndrome. Nature communications. 2015;6:8470. pmid:26415585
  25. 25. Wang L, Ni X, Covington KR, Yang BY, Shiu J, Zhang X, et al. Genomic profiling of Sezary syndrome identifies alterations of key T cell signaling and differentiation genes. Nature genetics. 2015 Dec;47(12):1426–34. pmid:26551670
  26. 26. Prasad A, Rabionet R, Espinet B, Zapata L, Puiggros A, Melero C, et al. Identification of gene mutations and fusion genes in patients with Sezary Syndrome. The Journal of investigative dermatology. 2016 Mar 30.
  27. 27. Choi J, Goh G, Walradt T, Hong BS, Bunick CG, Chen K, et al. Genomic landscape of cutaneous T cell lymphoma. Nature genetics. 2015 Sep;47(9):1011–9. pmid:26192916
  28. 28. Woollard WJ, Pullabhatla V, Lorenc A, Patel VM, Butler RM, Bayega A, et al. Candidate driver genes involved in genome maintenance and DNA repair in Sezary syndrome. Blood. 2016 Jun 30;127(26):3387–97. pmid:27121473
  29. 29. Lamprecht B, Kreher S, Mobs M, Sterry W, Dorken B, Janz M, et al. The tumour suppressor p53 is frequently nonfunctional in Sezary syndrome. Br J Dermatol. 2012 Aug;167(2):240–6. pmid:22384858
  30. 30. Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic acids research. 2010 Sep;38(16):e164. pmid:20601685
  31. 31. Theodorou I, Delfau-Larue M, Bigorne C, Lahet C, Cochet G, Bagot M, et al. Cutaneous T-cell infiltrates: analysis of T-cell receptor γ gene rearrangement by polymerase chain reaction and denaturing gradient gel electrophoresis. Blood. 1995;86:305–10. pmid:7795237
  32. 32. Langerak AW, Groenen PJ, Bruggemann M, Beldjord K, Bellan C, Bonello L, et al. EuroClonality/BIOMED-2 guidelines for interpretation and reporting of Ig/TCR clonality testing in suspected lymphoproliferations. Leukemia. 2012 Oct;26(10):2159–71. pmid:22918122
  33. 33. Marty M, Prochazkova M, Laharanne E, Chevret E, Longy M, Jouary T, et al. Primary cutaneous T-cell lymphomas do not show specific NAV3 gene deletion or translocation. The Journal of investigative dermatology. 2008 Oct;128(10):2458–66. pmid:18509363
  34. 34. Chevret E, Andrique L, Prochazkova-Carlotti M, Ferrer J, Cappellen D, Laharanne E, et al. Telomerase functions beyond telomere maintenance in primary cutaneous T-cell lymphoma. Blood. 2014 Mar 20;123(12):1850–9. pmid:24470590
  35. 35. Pham-Ledard A, Prochazkova-Carlotti M, Andrique L, Cappellen D, Vergier B, Martinez F, et al. Multiple genetic alterations in primary cutaneous large B-cell lymphoma, leg type support a common lymphomagenesis with activated B-cell-like diffuse large B-cell lymphoma. Mod Pathol. 2014 Mar;27(3):402–11. pmid:24030746
  36. 36. Hirsch FR, Varella-Garcia M, Bunn PA Jr., Di Maria MV, Veve R, Bremmes RM, et al. Epidermal growth factor receptor in non-small-cell lung carcinomas: correlation between gene copy number and protein expression and impact on prognosis. J Clin Oncol. 2003 Oct 15;21(20):3798–807. pmid:12953099
  37. 37. Woollard WJ, Kalaivani NP, Jones CL, Roper C, Tung L, Lee JJ, et al. Independent Loss of Methylthioadenosine Phosphorylase (MTAP) in Primary Cutaneous T-Cell Lymphoma. The Journal of investigative dermatology. 2016 Feb 9.
  38. 38. Alexandrov LB, Nik-Zainal S, Wedge DC, Aparicio SA, Behjati S, Biankin AV, et al. Signatures of mutational processes in human cancer. Nature. 2015 Aug 22;500(7463):415–21.
  39. 39. Whibley C, Pharoah PD, Hollstein M. p53 polymorphisms: cancer implications. Nature reviews. 2009 Feb;9(2):95–107. pmid:19165225
  40. 40. McGirt LY, Jia P, Baerenwald DA, Duszynski RJ, Dahlman KB, Zic JA, et al. Whole-genome sequencing reveals oncogenic mutations in mycosis fungoides. Blood. 2015 Jul 23;126(4):508–19. pmid:26082451
  41. 41. McGregor JM, Crook T, Fraser-Andrews EA, Rozycka M, Crossland S, Brooks L, et al. Spectrum of p53 gene mutations suggests a possible role for ultraviolet radiation in the pathogenesis of advanced cutaneous lymphomas. The Journal of investigative dermatology. 1999 Mar;112(3):317–21. pmid:10084308
  42. 42. Marks DI, Vonderheid EC, Kurz BW, Bigler RD, Sinha K, Morgan DA, et al. Analysis of p53 and mdm-2 expression in 18 patients with Sezary syndrome. British journal of haematology. 1996 Mar;92(4):890–9. pmid:8616082
  43. 43. Campbell JJ, Clark RA, Watanabe R, Kupper TS. Sezary syndrome and mycosis fungoides arise from distinct T-cell subsets: a biologic rationale for their distinct clinical behaviors. Blood. 2010 Aug 5;116(5):767–71. pmid:20484084
  44. 44. Landau DA, Tausch E, Taylor-Weiner AN, Stewart C, Reiter JG, Bahlo J, et al. Mutations driving CLL and their evolution in progression and relapse. Nature. 2015 Oct 22;526(7574):525–30. pmid:26466571
  45. 45. Leroy B, Anderson M, Soussi T. TP53 mutations in human cancer: database reassessment and prospects for the next decade. Human mutation. 2014 Jun;35(6):672–88. pmid:24665023
  46. 46. Iggo R, Rudewicz J, Monceau E, Sevenet N, Bergh J, Sjoblom T, et al. Validation of a yeast functional assay for p53 mutations using clonal sequencing. The Journal of pathology. 2013 Dec;231(4):441–8. pmid:23897043
  47. 47. Gibson JF, Huang J, Liu KJ, Carlson KR, Foss F, Choi J, et al. Cutaneous T-cell lymphoma (CTCL): Current practices in blood assessment and the utility of T-cell receptor (TCR)-Vbeta chain restriction. Journal of the American Academy of Dermatology. 2016 May;74(5):870–7. pmid:26874819
  48. 48. van Doorn R, van Kester MS, Dijkman R, Vermeer MH, Mulder AA, Szuhai K, et al. Oncogenomic analysis of mycosis fungoides reveals major differences with Sezary syndrome. Blood. 2009 Jan 1;113(1):127–36. pmid:18832135
  49. 49. van Doorn R, Slieker RC, Boonk SE, Zoutman WH, Goeman JJ, Bagot M, et al. Epigenomic Analysis of Sezary Syndrome Defines Patterns of Aberrant DNA Methylation and Identifies Diagnostic Markers. The Journal of investigative dermatology. 2016 Apr 23.