Key Points
Richter syndrome has genomic complexity intermediate between chronic lymphocytic leukemia and diffuse large B-cell lymphoma.
Inactivation of TP53 and of CDKN2A is a main mechanism in the transformation to Richter syndrome.
Abstract
Richter syndrome (RS) occurs in up to 15% of patients with chronic lymphocytic leukemia (CLL). Although RS, usually represented by the histologic transformation to a diffuse large B-cell lymphoma (DLBCL), is associated with a very poor outcome, especially when clonally related to the preexisting CLL, the mechanisms leading to RS have not been clarified. To better understand the pathogenesis of RS, we analyzed a series of cases including 59 RS, 28 CLL phase of RS, 315 CLL, and 127 de novo DLBCL. RS demonstrated a genomic complexity intermediate between CLL and DLBCL. Cell-cycle deregulation via inactivation of TP53 and of CDKN2A was a main mechanism in the histologic transformation from CLL phase, being present in approximately one half of the cases, and affected the outcome of the RS patients. A second major subgroup was characterized by the presence of trisomy 12 and comprised one third of the cases. Although RS shared some of the lesions seen in de novo DLBCL, its genomic profile was clearly separate. The CLL phase preceding RS had not a generalized increase in genomic complexity compared with untransformed CLL, but it presented clear differences in the frequency of specific genetic lesions.
Introduction
Richter syndrome (RS), which is defined by the histologic transformation to a diffuse large B-cell lymphoma (DLBCL), occurs in up to 15% patients with chronic lymphocytic leukemia (CLL).1-3 RS is associated with a very poor outcome, especially when clonally related.1-3
Mechanisms leading to RS have not been clarified, and strong molecular and clinical predictors of transformation at the time of the diagnosis of CLL are needed. In the vast majority of the studies, authors have assessed RS for the genetic events known to occur in CLL and in DLBCL.2-9 Only a few, unbiased, genome-wide studies have been reported, which were performed on a limited number of samples.7,10 So far, the most important events recently elucidated regarding RS pathogenesis are TP53 inactivation and the activation of MYC pathway,3 which also could be mediated by somatic mutations activating NOTCH1.9
Better characterization of the genetics of RS could serve diverse purposes. The assessment of specific copy number (CN) changes underlying the CLL phase of patients with RS could help in an early identification of CLL patients with an enhanced transforming risk. The evaluation of genomic aberrations acquired during the transition from CLL phase to RS eventually would lead to a better understanding of the implicated biologic processes and potentially introduce better-tailored therapeutic approaches.
With these aims, we performed genome-wide DNA profiling on a large series of samples derived from the RS and CLL phases of patients with RS and compared them with untransformed CLL (CLL-U) and de novo DLBCL.
Materials and methods
Tumor panel
All clinical specimens were derived from pathological sites and obtained in the course of routine diagnostic procedures, before the initiation of therapy. CLL-phase and CLL-U samples were collected from peripheral blood at diagnosis. RS and DLBCL samples were collected from nodal or extranodal biopsies. Diagnoses were made following the recommended criteria.1 All RS cases represented DLBCL transformed from a typical CLL. The distinction between germinal center B-cell (GCB) and non-GCB type, on the basis of gene expression profiling, was available for 108 de novo DLBCL samples as previously published.11-13 Cases were selected on the basis of the availability of frozen material, with the fraction of neoplastic cells in the specimen representing >70% of overall cellularity as determined by morphologic and/or immunophenotypic studies. Informed consent was obtained in accordance to the Declaration of Helsinki following the procedures approved by the local ethical committees and institutional review boards of each participating institution. The study was approved by the Bellinzona Ethical Committee.
Genome-wide DNA profiling
Genome-wide DNA profiles were obtained from high-molecular-weight genomic DNA by use of the Genome-Wide Human SNP Array 6.0 (Affymetrix, Santa Clara, CA), as previously reported.14 The raw data reported in this article have been deposited in the Gene Expression Omnibus database (accession number GSE50252). For 73 of 127 DLBCL cases, CEL files were kindly provided by Laura Pasqualucci.11 The raw CN was extracted from CEL files as previously reported.14,15 Genomic profiles were segmented with the Fast first-derivative segmentation algorithm.15 Profiles were considered as poor quality and discarded from further analyses if they showed severe oversegmentation or no aberrations at all as evaluated by two independent investigators (I.K., F.B.). Minimal common regions (MCRs)16 were calculated. MCRs containing genes encoding the immunoglobulin heavy chain (IGHV) genes and the κ and λ light chains were discarded.13 Focal aberrations additionally were identified with the Genomic Identification of Significant Targets in Cancer 2.0 algorithm17 following the default criteria on the filtered and merged segments.
Mutational analysis
Rearrangements of the IGHV-D-J were amplified and sequenced as previously described.3 We documented the clonal relationship between the RS phase and the CLL phase by analyzing the IGHV-D-J rearrangements.3 Mutational analysis of the TP53 and NOTCH1 genes was performed as previously described.3,18 Mutational analysis of the CDKN2A gene was performed by the use of 5′-ACC GGA GGA AGA AAG AGG AG-3′ and 5′-AGA ATC GAA GCG CTA CCT GA-3′ as forward and reverse primers for the p16INK4a isoform (NM_000077) and 5′-TGG GTC CCA GTC TGC AGT TA-3′ and 5′-TAG CCT GGG CTA GAG ACG AA-3′ as forward and reverse primers for the p14ARF isoform (NM_058195).
Analysis of clinical data
Overall survival (OS) was calculated from the time of histologic transformation to the last follow-up or death of any cause. Log-rank test was used to investigate the impact on OS of categorical variables. The cumulative probability of OS was plotted as curve according to the Kaplan-Meier method. The χ2 test was used to compare differences in clinical parameters among subgroups. A P value of < 0.05 was considered as statistically significant. Statistical analyses were performed with Stata/SE v.12.1 (StataCorp, College Station, TX).
Results
RS has a genomic complexity intermediate between CLL and DLBCL
To better understand the pathogenesis of RS, we have analyzed, by using the same high-resolution genome wide-DNA profiling assay, tumor samples derived from: (1) RS, that is, DLBCL in patients with a previous documented history of CLL (n = 59, Table 1 19 ); (2) CLL phase of RS, that is, CLL samples from patients who have then undergone histologic transformation to DLBCL (n = 28, 15 with paired RS); (3) CLL, from patients who have not undergone histologic transformation to an aggressive lymphoma (n = 315); and (4) de novo DLBCL, that is, from patients without having a previous history of CLL (n = 127).
The median number of CN changes was three among CLL-U, 3.5 for CLL-phase, 8.5 for RS, and 16 for DLBCL. The differences were statistically significant between RS and CLL phase (P = .0035) or CLL-U (P < .0001) and between RS and DLBCL (P = .003). Similarly, the frequency of recurrent lesions also showed a progressive increase across the spectrum of the four diseases (Figure 1A-D).
RS presents a series of recurrent lesions
In RS, the four most prevalent MCRs, likely containing the candidate cancer genes, were deletions at 17p (affecting TP53), at 13q14.3 (DLEU2/MIR15A/MIR16-1), trisomy 12, and losses at 9p21 (CDKN2A; Table 2). Recurrent homozygous losses were detected at 9p21 (7%; CDKN2A) and 13q14.3 (3%; DLEU2/MIR15A/MIR16). The targets of the 17p and 9p deletions also were sequenced. TP53 was found mutated in 12 of 21 (57%) cases, and the mutations were concomitant to 17p loss in 5 of 12. No CDKN2A mutations were found among 16 cases analyzed (12 without 9p21 loss). There were amplifications targeting known oncogenes: JAK2 (two cases, in one also KDM4C/JMJD2C), MIR17HG (one case), MYC (one case), BCL2 (one case).
To identify significantly recurrently affected regions, we applied the GISTIC algorithm (Figure 1E, Table 3). The MIR17HG locus, coding for the mir17-92 miRNA cluster, at 13q31.3 was the only significant region of gain. On the converse, multiple regions of losses were identified, with, among the most significant, again deletions at 13q14.2-q14.3 (DLEU2/MIR15A/MIR16-1) and 9p21 (CDKN2A), but also at 14q24.2-q32.33.
CLL-phase presents a limited number of recurrent lesions
In CLL-phase samples, the most common MCRs were deletions at 17p (TP53) and at 13q14.3 (DLEU2/MIR15A/MIR16), followed by trisomy 12, 11q, and 8p losses and 2p gains (Table 2). The only recurrent homozygous deletion targeted the 13q14.3 locus. The GISTIC algorithm did not identify any significant gain and only two losses, affecting the DLEU2/MIR15A/MIR16 and the MGA loci, respectively (data not shown).
CDKN2A loss is the most common acquired lesions at transformation
To identify aberrations acquired at histologic transformation, we compared RS with CLL phase (Table 4). Losses at 9p21 (CDKN2A) were the most statistically significant difference, being absent in CLL phase and detected in up to 20% of RS samples. In three cases with paired CLL phase and RS, the loss was clearly acquired at transformation. As also evident in Figure 1, RS presented a series of additional lesions more commonly than CLL phase, although the differences did not reach the statistical significance: losses at 6q and 9q or gains at 8q, 13q, 15q, and 18q. Conversely, losses at 13q14.3 (DLEU2/MIR15A/MIR16), 17p (TP53) and 11q (ATM), 8p, and also trisomy 12 and gains at 2p (REL) were observed at an overall similar frequency in both CLL phase and RS. Among 15 paired samples, TP53 loss was present in six cases, and in one case the lesion was acquired at the time of transformation. Figure 2 shows the acquisition of multiple genomic aberrations at the transformation in three paired genomic profiles.
CLL in patients who develop RS exhibits a greater frequency of certain genomic aberrations compared with untransformed CLL
To identify lesions capable of identifying CLL patients at greater risk of RS transformation, we then compared the genomic profiles observed in the CLL-phase samples in patients who developed RS with the samples of CLL cases who with a median follow-up of 6.25 years have not transformed to RS (Table 4). The most significant differences between the two groups concerned a greater frequency of losses at 17p (TP53), 15q (MGA, including homozygous losses), and gains at 2p (MYCN, REL) in samples derived from CLL phase, which, conversely, presented less deletions targeting MIR15/MIR16 (13q14.3). Although not reaching statistical significance, also +12 and ATM losses were observed more frequently in CLL phase of RS-transformed cases than in CLL-U: 29% vs 18% and 11% vs 7%, respectively.
Because CLL genomic lesions show a differential distribution on the basis of the IGHV gene mutational status20 and the CLL phase of RS patients is mostly associated with unmutated IGHV,2,8 the genomic profiles of CLL phase were compared with CLL-U both with unmutated IGHV (uCLL-U, n = 111) and with mutated IGHV genes (mCLL-U, n = 199). Even when compared with the uCLL-U samples, CLL-phase samples more commonly displayed the loss of TP53 (32% vs 11%; P = .0046), but other lesions (−8p, RS 11% vs uCLL-U 4%; +2p, 14% vs 8%; −13q14.3, 32% vs 44%), although differentially affected, did not reach statistical significance. Similar frequencies were observed for 15q15 (MGA, RS 11% vs uCLL-U 7%; homozygous losses, 4% vs 1%) and 11q21-q23 (ATM, 11% vs 14%) losses, and for trisomy 12 (29% vs 24%).
Differences in distribution of genomic aberrations were largely more evident when CLL phase was compared with 199 mCLL-U cases. CLL phase presented more commonly 2p25.3-p12 gains (REL/MYCN, 14%% vs 1%; P < .00001), losses affecting 8p (11% vs 2%; P = .028), 11q21-q23.3 (ATM; 11% vs 2%; P = .028), 15q15 (MGA, 14% vs 2%, P = .0002), and 17p (TP53, 32% vs 4%; P < .00001); in contrast, they exhibited less frequently 13q14.3 deletions (DLEU2/MIR15A/MIR16B, 32% vs 58%; P = .0077).
In CLL, different patterns of 13q14.3 deletions have been described, mainly determined by concomitant loss or otherwise of both the MIR15/MIR16 and of the more centromeric RB1 locus.14,21-23 Here, we observed that the 13q14 losses also encompassed the RB1 locus more commonly in CLL phase and RS than in CLL-U: 56%, 63%, and 42%, respectively, but the differences were not statistically significant (P = .064).
RS shares some of the DLBCL lesions, but its profile remains separate
Finally, we compared RS with de novo DLBCL (Table 4). A large series of statistically significant differences were observed, largely attributable to a lower frequency in RS of frequently associated lesions of de novo DLBCL. The only lesions that appeared more commonly in RS were lesions typical for CLL (13q14.3 and 11q22.3 losses and +12) and the 7q31-q36 deletion. Lesions that were the most significantly overrepresented in RS vs CLL-phase or in the CLL-phase vs CLL-U, namely losses of CDKN2A and TP53, showed similar frequency in both RS and DLBCL: 21% vs 17% and 33% vs 24%, respectively.
Because two main DLBCL subtypes are recognized on the basis of gene expression profiling, GCB, and non-GCB DLBCL, segregating unique recurrent genomic lesions,16 we compared the RS vs the two individual DLBCL subgroups (supplemental Table 1; see the Blood Web site). The levels of statistical significance of the observed differences were largely maintained, showing differences with both GCB- and non−GCB-DLBCL.
Two main genetic pathways lead to RS
Considering the most frequent lesions, we found that the 59 RS cases could be separated in three main subgroups that were mutually exclusive (Figure 3A). A first group, accounting for one half of the RS cases (30/60), was characterized by TP53 inactivation (by loss or by somatic mutations) and/or CDKN2A loss, alongside MYC activation (five gains and two translocations), 13q14.3 loss (n = 9) and additional lesions. A second group presented almost exclusively +12 (17/60, 28%). A third group was heterogeneous in terms of genomic lesions (13/60, 22%). Seven had genomic aberrations other than TP53/CDKN2A inactivation or +12, often including some of the most recurrent lesions; six cases did not apparently present any aberration.
NOTCH1, which can be activated by somatic mutations in CLL and RS,9,24-27 was sequenced in 28 cases: seven mutations were detected: 5 of 7 (71%) mutations were detected in cases bearing +12 and were mutually exclusive with MYC activation.
The vast majority of patients with a high proliferation rate (Ki-67+ cells > 70%) presented TP53 and/or CDKN2 inactivation (79%, 11/14), and the latter group showed a trend for a greater proliferation rate than the remaining cases (P = .057). A high proliferation rate was observed in 61% of the cases (11/18) with inactivation of TP53 and/or CDKN2, in 33% (2/6) with +12 and in 17% (1/6) of the genetically heterogeneous group. Among the 24 cases with available follow-up data, the cases with TP53 and/or CDKN2A inactivation had poorer OS than the remaining patients (P = .0105; hazard ratio; 3.6, 95% confidence interval, 1.3-10.1; Figure 3B). No associations were found for any of the three RS subgroups with IGHV gene mutational status, B-cell receptor stereotypy, expression of BCL6, MUM1, CD5, CD23, or the clonal relatedness of the CLL phase with the RS phase.
Discussion
We analyzed a large series of RS and its preexisting CLL phase, untransformed CLL, and de novo DLBCL cases by high-density genome-wide DNA microarray to better define the genetic aberrations underlying RS and to understand its relationship with the other three conditions. Our results are indicating that: (1) RS has a genomic complexity intermediate between CLL and DLBCL, as suggested by the CN aberration rates; (2) being present in approximately one half of the cases, inactivation of TP53 and of CDKN2A is a major mechanism in the histologic transformation from CLL phase and can affect the outcome of patients with RS; (3) a second major subgroup defined by genetic lesions is characterized by the presence of trisomy 12 and comprises approximately one third of the cases; (4) RS shares some of the DLBCL lesions, but its profile remains clearly separate; (5) the CLL phase preceding RS does not have an overall increase in genomic complexity compared with CLL-U, but it presents clear differences in the frequency of specific genetic lesions.
In RS cases in the present study, both the median number of lesions per sample and the recurrent lesions were intermediate between CLL and DLBCL, whereas the profiles of the latter two diseases demonstrated findings similar to what is reported in the literature.9,11 RS bore more lesions than CLL-U or even CLL phase, which is in accordance with the greater biologic and clinical aggressiveness of the disease. However, the observation that RS showed less CN changes than DLBCL did not seem in keeping with the more aggressive clinical course in RS. In the immunodeficiency-related DLBCL, the presence of Epstein-Barr virus infection is associated with a lower number of genomic aberrations.28,29 Although RS arises in the preexisting immunodeficiency associated with CLL,1 Epstein-Barr virus infection does not play a major role in RS2,3 and is an unlikely explanation for our data. Epigenetic changes, such as specific methylation patterns,30 or still-unidentified somatic mutations, might explain the observed relative low genomic complexity in RS vs DLBCL.
The most common lesions in RS concerned TP53 losses, trisomy 12, 13q14.3 losses (DLEU2/MIR15A/MIR16B), and CDKN2A losses. All these lesions occurred in >20% of the cases and were accompanied by a series of additional aberrations, such as gains of the MIR17HG locus, losses at 14q32 (including TRAF3), at 15q15 (MGA), at 1p (including TNFRSF14), 11q (including ATM, BIRC3, MIR34B, MIR34C), and on the long arm of chromosome 7.
Apart from TP53 inactivation, trisomy 12 and 13q14.3 losses, already present in the CLL phase, the vast majority of the RS genomic aberrations are acquired at the time of transformation. The most frequent acquired lesions were deletions targeting CDKN2A, a known tumor suppressor gene involved in cell-cycle regulation. CDKN2A inactivation, by DNA losses and not by somatic mutations, was exclusively observed in RS and in DLBCL at similar frequencies, but neither in CLL-U nor in CLL phase. Importantly, in DLBCL, CDKN2A losses have been associated with a poor outcome and with a non−GCB-phenotype.16,31,32 In addition, there are previous suggestions of the importance of cell-cycle deregulation in RS.30,33-35 In our RS series, CDKN2A losses also were associated with poor outcome and occurred mostly concomitantly with TP53 inactivation.
The inactivation of either CDKN2A or TP53 was accompanied by additional genomic lesions (RB1 losses, MYC activation, gains of MIR17HG) and defined a subgroup of patients characterized by high proliferation rate and deregulation of cell cycle, apoptosis, senescence, and cellular metabolism. Of clinical relevance, these patients had a poorer outcome. Similarly to what has been very recently suggested for a subset of de novo DLBCL,35 it might be speculated that these patients might benefit by treatment with compounds targeting the dysregulated apoptotic machinery, inhibiting the cyclin-dependent kinases or reactivating TP53 activity.
TP53 inactivation and/or CDKN2A loss were mutually exclusive with trisomy 12, indicating that RS might develop via two main genetic pathways. The most common mechanism, used by at least one half of the cases, entails the accumulation of a series of genomic lesions, resulting in TP53 inactivation, MYC activation, and, finally, loss of CDKN2A. A second mechanism, occurring in one third of the patients, is driven by the early acquisition of trisomy 12. This lesion is followed by the acquisition of NOTCH1 mutations, in agreement with recent publications,9,24,25 and/or of a limited number of additional genomic aberrations not exclusively observed in this group of patients.
Approximately one fifth of the cases had neither TP53/CDKN2A inactivation nor trisomy 12 and exhibited a heterogeneous pattern of genomic aberrations, with individual patients bearing lesions (11q losses, 14q losses, MYC activation, NOTCH1 mutation, MIR17HG gain) also detected in the two major RS subgroups defined by TP53/CDKN2A losses or trisomy 12.
In agreement with the literature,2,3 the histology of all our cases of RS was that of DLBCL, and the vast majority of them presented a non-GCB phenotype. Overall, the RS genomic profile appeared statistically distant from both non−GCB-DLBCL and GCB-DLBCL but RS did not fully recapitulate the genomic lesions observed in DLBCL, neither in the GCB- nor the non−GCB-subgroups. RS shared some of the lesions commonly observed in DLBCL (such as gains of 1q, 3q, 13q, and 18q and losses at 9p and 17p), but other lesions, recurrent in DLBCL (such as gains of chromosomes 5, 7, or 6q losses) were clearly underrepresented in RS. Indeed, we had also previously reported that RS has an extremely low rate of deletions and mutations inactivating the two known tumor suppressor genes mapped at 6q, coding for PRDM1/BLIMP1 (6q21) and TNFAIP3/A20 (6q23),3 which prevent plasmacytic differentiation and activate the nuclear factor κB pathway, respectively.36
The lesions specifically detected in RS and not, or at lower frequency, in DLBCL comprised mainly the aberrations typically observed in CLL (+12, and 13q14.3 and 11q23 losses).20 The deletions targeting chromosome 13 often encompassed the RB1 gene, alongside DLEU2/MIR15A/MIR16B. This pattern is similar to the one recently described in DLBCL, in which the lesion seems to contribute to cell-cycle deregulation and also to the tumor immune escape,36 and compatible with the notion that that a wider 13q14 lesion, comprising also RB1, is associated with a more aggressive form of disease.14,21,23,38 Losses affecting the long arm of 7 losses were observed exclusively in RS; these losses were more telomeric than the deletions present in splenic marginal zone lymphomas,39 and their extension did not allow the identification of putative tumor suppressor genes.
RS cells usually carry unmutated IGHV genes; thus, it was not surprising that the genomic profiles of the paired CLL-phase samples were similar to that seen in uCLL-U. CLL phase presented an elevated frequency of specific genomic aberrations (such as 2p gains, losses at 8p, 15q, and TP53 inactivation) but not a general increase of genomic complexity. In our previous study of 323 unselected CLL cases, gains occurring at 2p, possibly affecting REL and/or MYCN, were associated with a poor outcome and also with a high risk of RS transformation20 : this observation is herein confirmed on a larger number of CLL-phase samples. Losses at 8p have been reported to be associated with TP53 inactivation and with a poorer outcome in lymphoid tumors, including CLL and DLBCL, often together with TP53 inactivation.20,40 Deletions at 15q affected MGA, a negative MYC regulator, which, thus, might be directly involved in CLL and RS pathogenesis.41,42
Finally, indications on new therapeutic targets for RS patients come from the observed genomic amplifications of 9p24.1 (JAK2, JMJD2C), 13q31.3 (MIR17HG), and 18q21.33 (BCL2), which are all directly involved in the positive regulation of MYC expression and activity,43,44 and, notably, amplifications of 8q24.21 (MYC) also were observed underlining the importance of the gene in RS pathogenesis. Moreover, MIR17HG also contributes to the activation of the PI3K/AKT/mTOR pathway.45 Our data provide the rational for the evaluation of new targeted therapies such JAK2 inhibitors, MYC inhibitors,46 BCL2 inhibitors, phosphatidylinositol 3-kinase inhibitors, AKT inhibitors, and mammalian target of rapamycin inhibitors.47
In conclusion, our data establish specific patterns of genomic lesions in RS patients with prognostic and potentially therapeutic relevance. Also, lesions differentiating RS from DLBCL and the RS CLL phase from untransformed CLL have been identified.
Partially presented as an oral communication at the 53rd annual meeting of the American Society of Hematology, San Diego, CA, December 10-13, 2011.
The data reported in this article have been deposited in the Gene Expression Omnibus database (accession number GSE50252).
The online version of this article contains a data supplement.
The publication costs of this article were defrayed in part by page charge payment. Therefore, and solely to indicate this fact, this article is hereby marked “advertisement” in accordance with 18 USC section 1734.
Acknowledgments
This work was supported by Oncosuisse OCS 02296-08-2008 (to F. Bertoni), the Helmut Horten Foundation (to F. Bertoni), the San Salvatore Foundation (to F. Bertoni), the Nelia et Amadeo Barletta Foundation (to F. Bertoni), Computational Life Sciences/Ticino in Rete (to F. Bertoni), and DFG-SFB1074, subproject B2 (to S.S.). E.C. was a recipient of an European Society for Medical Oncology (ESMO) Fellowship Grant.
Authorship
Contribution: E.C. interpreted data and cowrote the manuscript; A.R. performed genomic profiling; I.K. performed statistical analysis and interpreted data; P.M.V.R. performed statistical analysis; D.R. and G.G. designed research, collected and characterized tumor samples, and provided advice; J.C.S, D.O., K.S., T.P., F. Berger, K.H.Y., F.M., R.R., T.C.G., W.C.C., E.M.O., M. Lucioni, R.M., G.I., M. Ladetto, F.F., S.C., H.V., S.H.S., S.S., M.A.P, A.M., D.S., E.N., A.Z., V.G., G.B., and G.O. collected and characterized tumor samples; E.Z. provided advice and cowrote the manuscript; F. Bertoni designed research, performed statistical analysis, interpreted data, and cowrote the manuscript; and all authors approved the final manuscript.
Conflict-of-interest disclosure: The authors declare no competing financial interests.
Correspondence: Francesco Bertoni, MD, Lymphoma and Genomics Research Program, Institute of Oncology Research, via Vincenzo Vela 6, 6500 Bellinzona, Switzerland; e-mail: frbertoni@me.com.
References
Author notes
E.C. and A.R contributed equally to this study.