1. Introduction
Pleural mesothelioma (PM) is a rare and severe tumor affecting the pleura, strictly correlated with asbestos exposure [
1]. PM usually has a lag time equal to 30-40 years between asbestos exposure and the first clinical manifestations, with limited therapeutic options and an overall survival (OS) ranging from 4 to 14 months [
2,
3]. There are three main PM histotypes: epithelioid (70-80% of cases), biphasic (10-20%), and sarcomatoid (10%) [
4,
5,
6], with different therapeutic and prognostic implications. In fact, epithelioid PM patients are potentially eligible for surgery and for a multimodality approach, and have a better OS compared to other histotypes [
2,
7,
8].
The diagnosis of PM can be extremely challenging both for small tissue biopsies, if tissue invasion is not clearly assessable, and for pleural effusions [
9,
10]. In particular, the morphologic differentiation of PM from reactive mesothelial hyperplasia (MH) is not always possible for effusions so that a tissue sample can be required first to confirm malignancy and then to subtype PM [
5,
11]. According to the International Mesothelioma Interest Group (IMIG) recommendations, a cytological diagnosis of PM is possible by coupling morphological examination with ancillary tests [
12,
13]. BRCA1-associated protein 1 (BAP1) and cyclin-dependent kinase inhibitor 2A (
CDKN2A, better known as
p16) are the most valuable markers for the differential diagnosis of malignant and benign pleural lesions [
14]. Loss of the expression of BAP1, evaluated by immunohistochemistry (IHC), and the presence of
p16 homozygous deletion evaluated by Fluorescent In Situ Hybridization (FISH), are highly specific markers for PM, since these alterations have never been described in benign lesions [
14,
15,
16,
17,
18]. Furthermore, methylthioadenosine phosphorylase (MTAP) has been recently suggested as a valuable surrogate marker for
p16, more easily evaluable by IHC rather than by FISH. Indeed,
MTAP and
CDKN2A genes map closely on the same chromosome (9p21) and if a deletion is present both genes are usually involved. A high concordance has been reported between MTAP loss of expression and
p16 deletion [
19]. Taken separately, BAP1,
p16 and MTAP, have a low diagnostic sensitivity (40-60%) for PM, which can be greatly improved by their combination (70-90%) [
2,
5,
13,
18]. However, the absence of BAP1,
p16 or MTAP alterations does not rule out a possibility of PM [
2,
5,
14]. In addition, the cytological diagnosis of biphasic PM is even more challenging than that of epithelioid PM, considering that cells from the sarcomatoid component do not shed into the effusions [
14] and the efficacy of available markers has not been completely investigated for this subtype [
4,
9,
20].
In this context, our group has previously developed a 117-gene expression panel and a related classification model able to discriminate epithelioid PM from MH [
21]. This panel has already been tested and successfully compared with BAP1 and
p16 on two independent retrospective series: one including malignant (epithelioid PM) and benign pleural tissues, giving a sensitivity of 97% and a specificity of 100% [
22], and the other including only pleural effusions (both cell-blocks and stained smears) with a sensitivity and specificity both equal to 100% [
23]. We have already demonstrated that among the 117 genes included in the panel there are some histotype specific markers, whose expression can help discriminate between biphasic and epithelioid PM [
24].
In this study we performed a prospective validation of the 117-gene expression panel as a diagnostic tool for PM cytological diagnosis and we compared its performance with that of BAP1, p16 and MTAP. Our purpose was to improve the diagnostic algorithm for the cytological differential diagnosis not only of epithelioid PM and MH, but also of biphasic PM.
3. Results
3.1. Study population
In this study, a total of 85 patients with PM or MH were consecutively enrolled. In detail, the epithelioid PM cases were 36 (27 cell-blocks and 9 Papanicolaou stained smears; 10 females and 26 males, with a median age of 71 years); biphasic PM were 15 (12 cell-blocks and 3 Papanicolaou stained smears; 2 females and 13 males, with a median age of 71 years); MH were 34 (33 cell-blocks and 1 Papanicolaou stained smear; 17 females and 17 males, with a median age of 70 years). Eight samples (3 Papanicolaou stained smears, including 1 MH and 2 epithelioid PMs, and 5 cell-blocks, including 1 MH and 4 epithelioid PM cases) were not adequate (less than 10% of atypical cells, poor quantity and quality RNA or consumed cell-blocks) neither for IHC, FISH, nor for gene expression analyses.
Finally, 77 cases successfully underwent gene expression analysis, 68 were cell-blocks and were also screened for BAP1, MTAP and p16.
3.2. Gene expression analysis
All the 77 samples passed the quality checks and were considered adequate for statistical and bioinformatics analysis.
After PCA, the first two principal components (PC) accounted for more than 50% of variability and were used for plotting. MH and PM were separated on PC1 with some overlapping (
Figure 2A). Unsupervised clustering revealed different gene expression profiles for malignant and benign cases; in particular, it is possible to identify two main clusters, one including only malignant and the one including mainly benign cases (
Figure 2B).
Considering all the 77 samples (both cell-blocks and stained smears), the classification model had an area under the curve (AUC) equal to 0.90 (95% CI 0.81-0.97), a sensitivity of 0.84 (95% CI 0.71-0.93), a specificity of 1 (95% CI 0.94-1), an accuracy of 0.91 (95% CI 0.83-0.96), a positive predictive value (PPV) of 1 (95% CI 0.95-1), and a negative predictive value (NPV) of 0.82 (95% CI 0.71-0.91). Moreover, it showed a similar performance in the differential diagnosis of both epithelioid (27 out of 30 cases correctly classified as malignant) and biphasic PM (13 out of 15 cases correctly classified) from benign effusions.
Deregulated genes between epithelioid and biphasic PM effusions were determined and compared with the already reported deregulated genes between epithelioid and biphasic PM pleural tissues [
24]. The top consistently deregulated genes between epithelioid and biphasic PM cases – both on pleural tissues and effusions – were
MSLN, CLDN15 and
CFB (Supplementary
Table 1). The score built by averaging the expression levels of these genes showed an accuracy of 0.80 (95% CI 0.63-0.90) in discriminating epithelioid and biphasic PM directly on pleural effusions.
3.3. BAP1, MTAP, p16 tests
BAP1 IHC was successfully executed on all 68 cell-blocks; MTAP IHC provided assessable results in 66 out of 68 cases, while p16 FISH in 65 out of 68 cases. In detail, two cases (both MHs) were indeterminate for MTAP because of inadequate staining in internal control positive cells and/or heterogeneous staining in mesothelial cells. Three cases (1 epithelioid PM and 2 biphasic PMs) were not evaluable for p16, probably due to their suboptimal pre-analytical conditions (i.e., fixation time, storage).
Overall, sixty-five cases provided clear results for BAP1, MTAP,
p16 and for the 117-gene expression panel. These cases were used to compare the performance of single biomarkers and their combinations, as reported in
Table 1 and in
Figure 3. Interestingly, 8 PM cases showed neither BAP1 nor MTAP/
p16 alterations. Seven out of 8 misdiagnosed PM cases were correctly classified as malignant using the 117-gene expression panel.
4. Discussion
Pleural effusions are among the first clinical manifestations of PM and can provide important diagnostic information [
4,
11,
13,
14]. However, the diagnosis of PM on pleural effusions is still a controversial issue and several factors must be taken into account: cytomorphologic features overlap between PM and MH cells and this fact, along with the lack of standardization [
9] in the collection and preparation procedures, can result in variable interpretations across cytopathologists. The introduction of BAP1, MTAP and
p16 tests has greatly improved the diagnostic accuracy of pleural effusions, but it is necessary to perform at least two of these tests (i.e. BAP1 and MTAP) to reach a good sensitivity [
4,
5,
14,
17,
25,
26]. Although theoretically evaluable both on cell-blocks and smears, IHC tests are standardized mainly on cell-blocks, and stained smears are not always sufficient to perform more than one molecular test. In addition, interpretation challenges are well-known (i.e. imperfect internal control staining, mix of nuclear and cytoplasmic staining), and heterogeneous staining is frequently observed in overlapping cells, typical of cytology samples [
19]. BAP1 and MTAP IHC have been reported to provide ambiguous results in 5% of tissue samples, and this rate is equal to 10% for cell-blocks [
5,
19,
27]. In these cases it is necessary to perform a biopsy, which delays the start of treatment. Moreover, not all patients tolerate invasive diagnostic procedures [
9,
10].
In this study, we validated the 117-gene panel as a diagnostic tool for the cytological diagnosis of PM: as a single test, this panel reached a sensitivity and a diagnostic accuracy comparable to those of combined BAP1 and MTAP evaluation. Data interpretation is more objective and the analysis can be executed directly on stained-smears, just after pathological evaluation, thus avoiding cell-block preparation and saving time. It has been reported that when cytological diagnosis is reliable and treatment can be immediately initiated median OS can be extended [
5,
28]. Moreover, the nCounter NanoString procedure is based on direct counting of the mRNA target molecules, with no amplification step [
28,
29], so that RNA quantity and quality, usually scarce from FFPE or stained smears, are not limiting factors and the failure rate is extremely low [
30]. In our study none of the analyzed 77 cases failed the gene expression analysis.
In this prospective cohort, we also confirmed the usefulness of BAP1 and MTAP IHC tests on pleural effusions, whose combination reached an excellent diagnostic accuracy. On the contrary, in our series the p16 FISH test did not improve diagnostic performance. However, 8 out of 35 PM cell-blocks (22.8%) would not have been correctly diagnosed by ancillary tests. Seven out of eight were correctly classified as malignant by the NanoString panel, thus reducing the rate of false negative cases to 2.8% and increasing both diagnostic accuracy and sensitivity.
A limit of this approach is that specificity may be slightly reduced, since the gene panel misclassified some MH samples. Consistently with our previous studies, we applied the same classification algorithm that uses all the 117 genes. There is still room for improvement both in terms of gene selection and testing of different algorithms.
We are aware that not all routine laboratories are equipped with a NanoString system and IHC is a more affordable first-level test. However, gene expression analysis can be considered before planning a biopsy in non-decisive pleural fluids, possibly referring to specialized centers.
Another advantage of the 117-gene panel is the possibility to evaluate more genes simultaneously and to obtain much more information alongside the differential diagnosis between malignant and benign pleural effusions. PM subtyping drives patients’ therapy and prognosis definition, but biphasic PM can be difficult to diagnose and to discriminate from other histotypes [
31]. It is composed by both epithelioid and sarcomatoid morphology [
4,
9], and the higher is the proportion of sarcomatoid component the worse is the prognosis [
31]. However, the identification of a spindle-cell component is not easy and a considerable interobserver variation in recognizing biphasic subtypes has been reported [
31,
32]. On pleural effusions it is impossible to morphologically discriminate between biphasic and epithelioid PMs, considering that sarcomatoid cells do not shed in the effusions [
14]. It has already been demonstrated that gene expression profiles on tumor tissues can improve PM subtyping, particularly for biphasic PM [
33,
34]; however, to the best of our knowledge, no data are available concerning the role of gene expression profiles in PM subtyping on pleural effusions. In a previous study, performed on a retrospective series of pleural tissues from patients with epithelioid, biphasic and sarcomatoid PM, we found that our panel includes some genes expressed in a histotype specific manner [
24]. In this work, we confirmed the deregulation of
CFB,
MSLN, and
CLDN15 on pleural effusions from epithelioid and biphasic PM. The score built by averaging the expression levels of these genes showed an accuracy of 0.80 in discriminating cytological specimens from the two most common histotypes. In this context, once malignancy has been assessed, the possibility of subtyping epithelioid and biphasic PM directly on pleural effusions can suggest another important clinical application of this gene expression panel.
It has to be acknowledged that the malignant rate of our study is probably higher than most consecutive series. This can be due both to the lack of sufficient material from some benign effusions and to the elevated number of patients with suspicion of malignancy referring to our institution.
Finally, pleural effusions are precious diagnostic material for mesothelioma and the 117-gene panel could improve the clinical management of PM patients both as a single test or in a stepwise diagnostic approach.