Next Article in Journal
Prospective Study of Low-Radiation and Low-Iodine Dose Aortic CT Angiography in Obese and Non-Obese Patients: Image Quality and Impact of Patient Characteristics
Next Article in Special Issue
The Use of Chest Magnetic Resonance Imaging in Malignant Pleural Mesothelioma Diagnosis
Previous Article in Journal
Applications of Optical Coherence Tomography in the Ocular Diagnosis: From the Tear Film to the Sclera
Previous Article in Special Issue
Unusual Histology in Mesothelioma: A Report of Two Cases with a Brief Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Gene Expression Analysis of Biphasic Pleural Mesothelioma: New Potential Diagnostic and Prognostic Markers

1
Unit of Pathological Anatomy, University Hospital of Pisa, 56100 Pisa, Italy
2
Department of Surgical, Medical, Molecular Pathology and Critical Area, University of Pisa, 56100 Pisa, Italy
3
Unit of Pneumology, University Hospital of Pisa, 56100 Pisa, Italy
4
Division of Thoracic Surgery, Department of Surgical, Medical and Molecular Pathology and Critical Care Medicine, University of Pisa, 56100 Pisa, Italy
5
Minimally Invasive and Robotic Thoracic Surgery, Robotic Multispecialty Center of Surgery, University Hospital of Pisa, 56100 Pisa, Italy
6
Pathology Unit, Department of Mental and Physical Health and Preventive Medicine, Università degli Studi della Campania “L Vanvitelli”, 80138 Naples, Italy
*
Author to whom correspondence should be addressed.
Diagnostics 2022, 12(3), 674; https://doi.org/10.3390/diagnostics12030674
Submission received: 4 February 2022 / Revised: 2 March 2022 / Accepted: 8 March 2022 / Published: 10 March 2022
(This article belongs to the Special Issue Recent Advances in Diagnosis of Malignant Mesothelioma)

Abstract

:
Biphasic is the second most common histotype of pleural mesothelioma (PM). It shares epithelioid and sarcomatoid features and is challenging to diagnose. The aim of this study was to identify biphasic PM markers to improve subtyping and prognosis definition. The expression levels of 117 cancer genes, evaluated using the nanoString system, were compared between the three major histotypes (epithelioid, sarcomatoid, and biphasic), and expression differences within biphasic PM were evaluated in relation to the percentage of epithelioid components. Biphasic PM overexpressed CTNNA1 and TIMP3 in comparison to sarcomatoid, and COL16A1 and SDC1 in comparison to epithelioid PM. CFB, MSLN, CLDN15, SERPINE1, and PAK4 were deregulated among all histotypes, leading to the hypothesis of a gradual expression from epithelioid to sarcomatoid PM. According to gene expression, biphasic PM samples were divided in two clusters with a significant difference in the epithelioid component. ADCY4, COL1A1, and COL4A2 were overexpressed in the biphasic group with a low percentage of epithelioid component. Survival analysis using TCGA data showed that high COL1A1 and COL4A2 expression levels correlate with poor survival in PM patients. Herein, we identified markers with the potential to improve diagnosis and prognostic stratification of biphasic PM, which is still an orphan tumor.

1. Introduction

Pleural mesothelioma (PM) is a rare and highly aggressive tumor, mostly associated with asbestos exposure [1,2]. PM is characterized by a long latency period (10–30 years) from asbestos exposure to the clinical onset of disease, with a median overall survival ranging from 12–18 months and a 5-year survival rate of about 5% [3,4]. According to the 5th World Health Organization (WHO) Thoracic Tumors classification, diffuse PM includes three major histotypes: epithelioid (70–80% of cases), sarcomatoid (10%), and biphasic or mixed (10–20%) [5]. An accurate histological classification is mandatory to inform treatment decision and constitutes the most reliable prognostic factor. In fact, the reported median overall survival rate for epithelioid, biphasic and sarcomatoid PM is 19, 12, and 4 months, respectively [6,7]. Patients with epithelioid mesotheliomas can be eligible for surgery and multimodality approach, while sarcomatoid tumors do not achieve benefit from surgery [1]. The management of patients suffering from biphasic mesothelioma is challenging since this histotype displays a large degree of heterogeneity. A PM is classified as biphasic if both epithelioid and sarcomatoid components are present, with each comprising at least 10% of the tumor, and the proportion of sarcomatoid component is the major prognostic indicator and treatment driver [8,9]. Currently, there is no recommendation on grading sarcomatoid features, and percentage cut-offs of sarcomatoid component clearly associated with poor prognosis have not been identified [8]. Furthermore, the identification of a spindle cell component as malignant is not easy, and a low inter-observer correlation has been reported in recognizing biphasic PM [8]. In addition, the histological classification of PM is even more difficult on small biopsies, with a concordance with surgical resection ranging from 72 to 83% [10]. Recently, Chirieac and collaborators compared pre-surgery biopsies and surgical resection specimens from 759 PM patients: 112 patients (19%), initially diagnosed as epithelioid, were reclassified as biphasic after surgical resection; 15 (11%) and 4 (3%) cases initially diagnosed as biphasic were reclassified as epithelioid and sarcomatoid, respectively; 7 (19%) cases initially diagnosed as sarcomatoid were reclassified as biphasic [10].
The identification of histotype biomarkers to refine the histological classification and patients’ risk stratification is warranted. In this context, gene expression analysis proved to be a useful approach to better subtype PM [11]. Bueno and collaborators, by performing a transcriptome analysis, identified four distinct molecular subtypes: sarcomatoid, epithelioid, biphasic–sarcomatoid and biphasic–epithelioid, with the epithelioid subtype showing the best overall survival. The molecular subtyping was highly concordant with the histological diagnosis for sarcomatoid PM. A total of 62% of cases histologically diagnosed as epithelioid were differently reclassified, mostly as biphasic–epithelioid cases. A total of 21% of histologically biphasic samples were reclassified as sarcomatoid, 33.8% as biphasic–epitheliod and 40% as biphasic–sarcomatoid [11]. In the same way, Alcala and collaborators, using RNA-seq data, found that molecular profiles explain the prognosis of PM better than histological classification, with significant differences in the expression of proangiogenic and immune checkpoint genes [12].
The aim of our study was to identify histotype-specific biomarkers by evaluating the expression levels of 117 cancer genes, previously reported as deregulated in PM [13]. We mainly focused on biphasic PM, in order to underline differences with epithelioid and sarcomatoid PM and to investigate new potential prognostic biomarkers correlated with the amount of epithelioid component.

2. Materials and Methods

2.1. Study Population

In this study, PM samples (biopsies and surgical resections) were retrospectively collected from the archives of the Pathology Unit, University Hospital of Pisa (Pisa, Italy) and of the Pathology Unit, “Luigi Vanvitelli” University Hospital of Naples (Naples, Italy). The study was conducted in accordance with the principles of Helsinki declaration of 1975 and approved by our institutional ethics committee (protocol code 9989, 2019). Informed written consent for publication was not required because no sensitive data that could identify patients were used in this study.
In detail, 25 biphasic (13 biopsies and 12 surgical resections), 32 epithelioid (10 biopsies and 22 surgical resections) and 18 sarcomatoid (11 biopsies and 7 surgical resections) samples were analyzed. All tumor samples were formalin-fixed and paraffin-embedded (FFPE); histological diagnosis and pathological features were independently reviewed by two pathologists (GA and GF) according to the WHO 5th edition criteria [5] (Figure 1). In particular, epithelioid mesotheliomas are usually composed of round, epithelioid cells, typically displaying eosinophilic cytoplasm and round nuclei with small nucleoli, although cytological atypia can occur. Epithelioid mesothelioma displays a wide range of histological patterns; common patterns are tubulopapillary, solid and trabecular, while adenomatoid is less common. Sarcomatoid mesotheliomas are instead composed of elongated/spindle cells arranged in solid sheets or within a fibrous stroma. Biphasic mesothelioma is a highly characteristic pattern containing a mixture of epithelioid and sarcomatoid areas within the same tumor. According to WHO criteria, in definitive resection specimens each pattern should constitute at least 10% of the neoplasm, while in small biopsy specimens any tumor can be diagnosed as biphasic PM, regardless of percentages of each component. However, reporting the percentage of sarcomatoid component is recommended because of potential implications for prognosis and therapeutic management. In the present study, for biphasic PM the percentage of epithelioid component was independently determined by two pathologists. All cases were discussed and the median value was considered for statistical evaluation. Since sarcomatoid areas may sometimes be difficult to distinguish from reactive stroma, we evaluated the fibrous proliferation using immunohistochemistry tests. In detail, immunohistochemical expression of mesothelial markers (i.e., calretinin, WT1, D2-40, and CK-PAN) or BAP-1 loss were helpful for diagnosis. Finally, the most representative paraffin block of each tumor was selected for gene expression analysis.

2.2. Gene Expression Analysis

Gene expression analysis was performed using a previously defined nCounter custom codeset including 117 PM target genes and 6 housekeeping genes, synthesized by nanoString Technologies (NanoString Technologies, Seattle, WA, USA), as previously described [13]. The 117 genes included in the panel belong to pathways with a crucial role in cancer development and progression (i.e., epithelial to mesenchimal transition (EMT), cell matrix, angiogenesis, DNA methylation, focal adhesion) [13].
For each sample, three FFPE unstained sections underwent a standard deparaffinization procedure; enrichment of tumor cells was performed by manual macrodissection, and RNA was purified using the Qiagen RNeasy FFPE kit (Qiagen, Hilden, Germany) according to the manufacturers’ instructions.
Quantity and quality of total RNA were assessed by a Xpose spectrophotometer (Trinean, Gentbrugge, Belgium). A total of 150 ng of RNA was used for nanoString analysis; hybridization was performed for 18 h at 65 °C in a SensoQuest thermal cycler (SensoQuest, Gottingen, Germany). The clean-up of samples and counts of digital reports were executed as described by the manufacturer (NanoString Technologies).

2.3. Data Analysis

Raw data were normalized using the nSolver v.4.0 (NanoString Technologies), and batch effects were removed by using the ComBat function of the Surrogate Variable Analysis (sva) Bioconductor package v.3.42.0. To observe gross similarity patterns of expression, hierarchical clustering and principal component analysis were performed using the packages heatmap3 v.1.1.9 and factoextra v.1.0.7 respectively. Differentially expressed genes (DEG) were computed by the moderate t-statistics and the Benjamini–Hochberg correction following the procedures of the limma package v.3.50.0. A false discovery rate (FDR) of 0.05 was considered significant. Common deregulated genes among contrasts were presented using a Venn diagram. To highlight gene expression differences within biphasic PM, non-negative matrix factorization algorithm was applied comparing the results with a random expression matrix. The best fitting clustering was used to predict the biphasic PM group based on gene expression; the analyses were performed using the NMF package v.0.23.0. The Mann–Whitney U test was used for testing differences according to the epithelial component.
For validation, TCGA data on Mesothelioma [11] were downloaded from cBioPortal (https://www.cbioportal.org/, last accessed on 10 December 2021) [14,15]. Patients were dichotomized based on median expression value of the selected genes. Survival curves were built using the Kaplan–Meier method and compared by the log-rank test following the procedures of the survival v.3.2-13 and survminer v.0.4.9 packages. All analyses and plots were generated in R v.4.1.2 environment (https://www.r-project.org/, last accessed on 15 December 2021).

3. Results

3.1. Differentially Expressed Genes between Histotypes

All RNA samples were adequate for gene expression analysis and none failed the nanoString test.
Principal component analysis indicated the existence of specific expression patterns for epithelioid and sarcomatoid PM, whereas biphasic PM showed a partial overlap with the other histotypes (Figure 2). The analysis of differentially expressed genes revealed 34 genes deregulated between biphasic and epithelioid PM, 68 between sarcomatoid and epithelioid PM, and 34 between sarcomatoid and biphasic PM (Tables S1–S3; Figures S1–S3). The overlapping of deregulated genes among contrasts is reported in the Venn diagram (Figure 3). Only four genes were specifically up-regulated in biphasic in comparison to sarcomatoid (CTNNA1 and TIMP3) and epithelioid PM (COL16A1 and SDC1). Twelve genes were down-regulated only in sarcomatoid in comparison to epithelioid PM (ACSL1, CDK7, EEF2, EGR3, EIF4G1, ESR2, HEG1, LGALS3BP, PPARA, PTGS2, SOD1, and TPPP) and 5 were up-regulated only in sarcomatoid compared to epithelioid PM (AURKA, MMP7, PDGFRB, PLK2, UBE2T). Five out of one-hundred-and-seventeen genes were deregulated among all groups (CFB, MSLN, CLDN15, SERPINE1, and PAK4) (Figure 4).

3.2. Gene Expression Differences within Biphasic PM

The best fitting clustering by the non-negative matrix factorization algorithm (Figures S4 and S5) predicted two biphasic PM groups based on gene expression. The percentage of the epithelioid component, independently evaluated, was significantly different between the two identified biphasic PM clusters, p = 0.03 (Figure 5). The comparison of gene expression levels between the two identified biphasic PM clusters revealed three genes significantly up-regulated in the cluster with a low percentage of epithelioid component: ADCY4, COL1A1, and COL4A2. Survival analysis using TCGA data showed that high COL1A1 and COL4A2 expression levels are correlated with a worse overall survival in PM patients (Figure 6).

4. Discussion

The histological subtyping of diffuse pleural mesothelioma is still the major prognostic factor, which is crucial for treatment planning [16]. Patients with epithelioid PM can be eligible for surgery, while those with sarcomatoid PM do not benefit from surgery [3,17]. Since biphasic PM shares both epithelioid and sarcomatoid features, it is hard both to diagnose and to treat [8,18]. Several studies have pointed out the diagnostic difficulties related to the pathological identification of biphasic PM, and the determination of sarcomatoid and epithelioid components, which greatly impact prognosis [8,10,18]. However, only few studies have been specifically focused on this histotype, mainly due to its rarity, and the identification of diagnostic and prognostic biomarkers is urgently needed.
Herein, we compared expression levels of 117 genes with a key role in cancer development and progression between the three major PM histotypes. As expected, expression levels mainly differ between sarcomatoid and epithelioid PM. Biphasic PM showed a heterogeneous gene expression profile, which overlaps to that of epitheliod PM in the majority of cases. Only four genes were specifically up-regulated in biphasic in comparison to sarcomatoid (CTNNA1 and TIMP3) and epithelioid PM (COL16A1 and SDC1). These deregulated genes impact on EMT process, leading from epithelial to mesenchymal phenotype [19]. In particular, CTNNA1 and TIMP3 are both involved in the inhibition of EMT [20,21]. COL16A1 encodes the alpha chain of type XVI collagen and plays a role in maintaining the integrity of the extracellular matrix [22]. SDC1 mediates cell binding, cell signaling and cytoskeletal organizations, and it can have a different impact according to cancer type. For instance, in gastric and colorectal cancer a reduced expression of SDC1 is associated with advanced stages, whereas in pancreatic and breast cancer its overexpression was found to promote cell growth and proliferation [23].
Interestingly, 5 out of 117 genes were deregulated among all groups, thus suggesting that their expression may constitute a gradient from epithelioid to sarcomatoid PM. Again, almost all of them are involved in EMT, which can give intermediate cell states between the epithelial and mesenchymal ones [24].
SERPINE1 was up-regulated in sarcomatoid in comparison to both epithelioid and biphasic, and in biphasic in comparison to epithelioid PM. This protein belongs to the serine proteinase inhibitor (serpin) superfamily; it is the main inhibitor of tissue plasminogen activator, essential for extracellular matrix remodeling, and urokinase. SERPINE1 can induce epithelial to mesenchimal transition and the acquisition of stem cell properties, thus promoting the transition from a proliferative to an invasive tumor phenotype. SERPINE1 expression is correlated also to the activation of hypoxia related factors; it is able to promote tumor cell migration and metastasis, and it is usually associated with a poor prognosis [25,26,27].
CFB, MSLN and CLDN15 were all up-regulated in epithelioid compared to sarcomatoid and biphasic PM, and in biphasic in comparison to sarcomatoid PM. CLDN15 and MSLN are well known epithelial markers expressed in epithelioid PM. CLDN15 is an integral membrane protein; it is a component of tight junctions and plays a crucial role in maintaining cell polarity and signal transductions [28,29]. It has already been reported as the most significantly up regulated gene in epithelioid in comparison to other PM histotypes [11]; indeed, it is usually suppressed in cells undergoing EMT [30]. MSLN gene encodes a cell-surface glicoprotein, expressed at low levels on normal mesothelial cells. MSLN can promote tumor invasion and malignant transformation, and it is highly expressed in epithelioid and biphasic PM, but not in sarcomatoid histotype [31,32,33]. To date, soluble mesothelin is the only tumor biomarker to receive US Food and Drug Administration approval for clinical use in mesothelioma, and it is under evaluation also as therapeutic target [31,33,34]. CFB encodes the complement factor B (CFB) involved in the alternative pathway of complement activation. Complement-related proteins are involved in cell-cell and cell-stroma interactions [35] and influence tumor microenvironment and tumor progression. Shimazaki et al. found a correlation between the expression of stromal CFB and an enrichment in immunosuppressive regulatory T-cells, myeloid-derived suppressor cells and tumor-associated macrophages in pancreatic ductal adenocarcinoma [35]. A low CFB expression has been associated with a poor overall survival and a more aggressive disease in different solid tumors, including lung adenocarcinoma [36] and pancreatic ductal adenocarcinoma [35]. In this study a CFB low expression was found in sarcomatoid PM, the most aggressive histotype.
PAK4 was up-regulated in biphasic PM in comparison to both sarcomatoid and epithelioid histotypes and in epithelioid in comparison to sarcomatoid PM. It encodes a serine threonine kinase with a pro-oncogenic function when overexpressed in human cells. PAK4 is involved in carcinogenesis by promoting cell proliferation, invasion, migration and anti-apoptotic activity. It has a critical role in actin assembly and in the regulation of cytoskeleton dynamics [37]. It has been demonstrated that tumors overexpressing PAK4 have an activated Wnt/b-catetin signaling pathway, and activated PAK4 promotes anchorage independent cell-growth [38,39]. In addition, PAK4 is involved in the regulation of the anti-tumor immune response and is usually associated with a lack of response to PD-1 blockade. Interestingly, PAK4 and PD-L1 over-expression is mediated by the same signaling pathways (mainly PI3K/AKT, and Wnt/β-catenin), and PAK4 inhibition can improve response to immunotherapy [38,39].
The discrimination between biphasic and the other PM histotypes is not always feasible, particularly on small biopsies [10], and the identification of a spindle cell component as malignant can be difficult [8,18]. In this context, the up-mentioned genes can work as diagnostic biormarkes for PM subtyping.
In addition, it has already been demonstrated that gene expression profiles are able to improve also prognostic definition, and that biphasic PM cannot be considered a unique entity, from a clinical nor from a biological point of view [11,12].
The non-negative matrix factorization algorithm identified two clusters of biphasic PM on the basis of gene expression. These two groups had a different percentage of epithelioid component and express differently three genes, namely ADCY4, COL1A1 and COL4A2, whose up-regulation is associated with poor survival of PM patients as showed by TCGA data analysis. In our study, all these genes were up-regulated in the cluster of biphasic PM with a lower percentage of epithelioid component, which is usually associated with a worse prognosis [6,8]. From a biological point of view, COL1A1 and COL4A2 have a key role in the organization of extracellular structure, angiogenesis, metallopeptidase activity and EMT promotion [40,41]. Recently, it has been proved, in a murine model, that inhibition of collagen production can delay malignant mesothelioma tumor growth [42]. In addition, COL1A1 has already been described as a potential prognostic marker for mesothelioma, with a high expression correlated with a poorer overall survival [41]. Furthermore, Zhang and collaborators found that in PM tissues COL1A1 expression is negatively and positively correlated with neutrophil and macrophages infiltration, respectively. High levels of macrophages infiltration ere indicative of a worse clinical outcome for PM patients [41].
ADCY4 is a member of the family of adenylate cyclases responsible for the formation of the secondary messenger cyclic adenosine (cAMP). ADCY4 expression levels have been correlated with overall survival in some cancers, such as breast cancer and in lung adenocarcinoma [43,44]. ADCY4 has also a role in inflammation: cAMP is involved in the regulation of caspase 11 inflammasome activation and this impacts on perpetuation or resolution of inflammation itself [43]. For these reasons, ADCY4, COL1A1 and COL4A2 can be valuable candidate as prognostic markers for biphasic PM.
In this study the sample size is limited, but it is relatively high considering the rarity of these tumors. In addition, it was not possible to perform survival analysis on our retrospective cohort because clinical data were not available; TCGA cases include mostly epithelioid PM, with few biphasic and sarcomatoid cases. Consequently, a prospective validation on biphasic PM should be performed. Nevertheless, the biological role, previous reports and the impact on survival of PM patients support our findings.
Although gene expression analysis by the nanoString system proved to be suitable for FFPE samples, it requires specific equipment and expertise. In this context, a prospective validation of identified biomarkers using also more common techniques, such as immunohistochemistry, is warranted both to confirm obtained results and to favor their eventual introduction in clinical practice.
In conclusion, we identified potential diagnostic markers, such as CTNNA1, TIMP3, COL16A1, SDC1, CFB, MSLN, CLDN15, SERPINE1, and PAK4, which could be helpful in the discrimination of biphasic PM from both epithelioid and sarcomatoid subtypes. We propose also ADCY4, COL1A1 and COL4A2 as potential prognostic markers in biphasic PM and probably in all PM irrespective of the histotype. The validation of these markers could offer additional tools to improve the management of this lethal and orphan disease.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/diagnostics12030674/s1. Figure S1: Volcano plot. Biphasic versus epithelioid mesothelioma; Figure S2: Volcano plot. Sarcomatoid versus biphasic mesothelioma; Figure S3: Volcano plot. Sarcomatoid versus epithelioid mesothelioma; Figure S4: Non-negative matrix factorization algorithm performance using the 117 genes compared to a random matrix; Figure S5: Consensus plot. Best fitting clustering by the non-negative matrix factorization algorithm. Table S1: Differentially expressed genes between biphasic and epithelioid mesothelioma; Table S2: Differentially expressed genes between sarcomatoid and biphasic mesothelioma; Table S3: Differentially expressed genes between sarcomatoid and epithelioid mesothelioma.

Author Contributions

Conceptualization, R.B., A.M.P., G.A. and G.F.; selection of cases, A.C., F.M., M.L. and R.F.; pathological revision, G.A., C.D., A.P., R.F. and G.F.; molecular analyses, R.B.; data analysis, A.M.P.; original draft preparation, R.B. and A.M.P.; review and supervision, G.A. and G.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded in part by The Kazan McClain Partners’ Foundation (micro research grant IMIG 2018 attendees).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by our Institutional Review Board (protocol code 9989, date of approval 20 February 2019).

Informed Consent Statement

Informed written consent for publication was not required because no sensitive and clinical data that can identify patients were used in this study.

Data Availability Statement

The data presented in this study are contained within the article or Supplementary Material.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Opitz, I. Management of Malignant Pleural Mesothelioma-The European Experience. J. Thorac. Dis. 2014, 6 (Suppl. 2), S238–S252. [Google Scholar] [PubMed]
  2. Robinson, B.M. Malignant Pleural Mesothelioma: An Epidemiological Perspective. Ann. Cardiothorac. Surg. 2012, 1, 491–496. [Google Scholar] [PubMed]
  3. Opitz, I.; Friess, M.; Kestenholz, P.; Schneiter, D.; Frauenfelder, T.; Nguyen-Kim, T.D.L.; Seifert, B.; Hoda, M.A.; Klepetko, W.; Stahel, R.A.; et al. A New Prognostic Score Supporting Treatment Allocation for Multimodality Therapy for Malignant Pleural Mesothelioma. J. Thorac. Oncol. 2015, 10, 1634–1641. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Sinn, K.; Mosleh, B.; Hoda, M.A. Malignant Pleural Mesothelioma: Recent Developments. Curr. Opin. Oncol. 2021, 33, 80–86. [Google Scholar] [CrossRef] [PubMed]
  5. WHO Classification of Tumours Editorial Board. Thoracic Tumours, 5th ed.; International Agency for Research on Cancer: Lyon, France, 2021. [Google Scholar]
  6. Hung, Y.P.; Chirieac, L.R. Pathology of Malignant Pleural Mesothelioma. Thorac. Surg. Clin. 2020, 30, 367–382. [Google Scholar] [CrossRef] [PubMed]
  7. Meyerhoff, R.R.; Yang, C.-F.J.; Speicher, P.J.; Gulack, B.C.; Hartwig, M.G.; D’Amico, T.A.; Harpole, D.H.; Berry, M.F. Impact of Mesothelioma Histologic Subtype on Outcomes in the Surveillance, Epidemiology, and End Results Database. J. Surg. Res. 2015, 196, 23–32. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Galateau Salle, F.; Le Stang, N.; Nicholson, A.G.; Pissaloux, D.; Churg, A.; Klebe, S.; Roggli, V.L.; Tazelaar, H.D.; Vignaud, J.M.; Attanoos, R.; et al. New Insights on Diagnostic Reproducibility of Biphasic Mesotheliomas: A Multi-Institutional Evaluation by the International Mesothelioma Panel From the MESOPATH Reference Center. J. Thorac. Oncol. 2018, 13, 1189–1203. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Travis, W.D. The 2015 WHO Classification of Lung Tumors. Pathologe 2014, 35 (Suppl. 2), 188. [Google Scholar] [CrossRef] [PubMed]
  10. Chirieac, L.R.; Hung, Y.P.; Foo, W.C.; Hofer, M.D.; VanderLaan, P.A.; Richards, W.G.; Sugarbaker, D.J.; Bueno, R. Diagnostic Value of Biopsy Sampling in Predicting Histology in Patients with Diffuse Malignant Pleural Mesothelioma. Cancer 2019, 125, 4164–4171. [Google Scholar] [CrossRef] [PubMed]
  11. Bueno, R.; Stawiski, E.W.; Goldstein, L.D.; Durinck, S.; De Rienzo, A.; Modrusan, Z.; Gnad, F.; Nguyen, T.T.; Jaiswal, B.S.; Chirieac, L.R.; et al. Comprehensive Genomic Analysis of Malignant Pleural Mesothelioma Identifies Recurrent Mutations, Gene Fusions and Splicing Alterations. Nat. Genet. 2016, 48, 407–416. [Google Scholar] [CrossRef] [PubMed]
  12. Alcala, N.; Mangiante, L.; Le-Stang, N.; Gustafson, C.E.; Boyault, S.; Damiola, F.; Alcala, K.; Brevet, M.; Thivolet-Bejui, F.; Blanc-Fournier, C.; et al. Redefining Malignant Pleural Mesothelioma Types as a Continuum Uncovers Immune-Vascular Interactions. EBioMedicine 2019, 48, 191–202. [Google Scholar] [CrossRef] [Green Version]
  13. Bruno, R.; Alì, G.; Giannini, R.; Proietti, A.; Lucchi, M.; Chella, A.; Melfi, F.; Mussi, A.; Fontanini, G. Malignant Pleural Mesothelioma and Mesothelial Hyperplasia: A New Molecular Tool for the Differential Diagnosis. Oncotarget 2017, 8, 2758–2770. [Google Scholar] [CrossRef] [Green Version]
  14. Cerami, E.; Gao, J.; Dogrusoz, U.; Gross, B.E.; Sumer, S.O.; Aksoy, B.A.; Jacobsen, A.; Byrne, C.J.; Heuer, M.L.; Larsson, E.; et al. The CBio Cancer Genomics Portal: An Open Platform for Exploring Multidimensional Cancer Genomics Data. Cancer Discov. 2012, 2, 401–404. [Google Scholar] [CrossRef] [Green Version]
  15. Gao, J.; Aksoy, B.A.; Dogrusoz, U.; Dresdner, G.; Gross, B.; Sumer, S.O.; Sun, Y.; Jacobsen, A.; Sinha, R.; Larsson, E.; et al. Integrative Analysis of Complex Cancer Genomics and Clinical Profiles Using the CBioPortal. Sci. Signal. 2013, 6, pl1. [Google Scholar] [CrossRef] [Green Version]
  16. Wadowski, B.; De Rienzo, A.; Bueno, R. The Molecular Basis of Malignant Pleural Mesothelioma. Thorac. Surg. Clin. 2020, 30, 383–393. [Google Scholar] [CrossRef]
  17. Scherpereel, A.; Opitz, I.; Berghmans, T.; Psallidas, I.; Glatzer, M.; Rigau, D.; Astoul, P.; Bölükbas, S.; Boyd, J.; Coolen, J.; et al. ERS/ESTS/EACTS/ESTRO Guidelines for the Management of Malignant Pleural Mesothelioma. Eur. Respir. J. 2020, 55, 1900953. [Google Scholar] [CrossRef]
  18. Brcic, L.; Kern, I. Clinical Significance of Histologic Subtyping of Malignant Pleural Mesothelioma. Transl. Lung Cancer Res. 2020, 9, 924–933. [Google Scholar] [CrossRef]
  19. Ramundo, V.; Zanirato, G.; Aldieri, E. The Epithelial-to-Mesenchymal Transition (EMT) in the Development and Metastasis of Malignant Pleural Mesothelioma. Int. J. Mol. Sci. 2021, 22, 12216. [Google Scholar] [CrossRef]
  20. Chi, Q.; Xu, H.; Song, D.; Wang, Z.; Wang, Z.; Ma, G. α-E-Catenin (CTNNA1) Inhibits Cell Proliferation, Invasion and EMT of Bladder Cancer. Cancer Manag. Res. 2020, 12, 12747–12758. [Google Scholar] [CrossRef]
  21. Guan, Z.; Zhang, J.; Song, S.; Dai, D. Promoter Methylation and Expression of TIMP3 Gene in Gastric Cancer. Diagn. Pathol. 2013, 8, 110. [Google Scholar] [CrossRef] [Green Version]
  22. Grässel, S.; Bauer, R.J. Collagen XVI in Health and Disease. Matrix Biol. 2013, 32, 64–73. [Google Scholar] [CrossRef] [PubMed]
  23. Liao, S.; Liu, C.; Zhu, G.; Wang, K.; Yang, Y.; Wang, C. Relationship between SDC1 and Cadherin Signalling Activation in Cancer. Pathol. Res. Pract. 2020, 216, 152756. [Google Scholar] [CrossRef] [PubMed]
  24. Dongre, A.; Weinberg, R.A. New Insights into the Mechanisms of Epithelial–Mesenchymal Transition and Implications for Cancer. Nat. Rev. Mol. Cell Biol. 2019, 20, 69–84. [Google Scholar] [CrossRef] [PubMed]
  25. Duffy, M.; McGowan, P.; Gallagher, W. Cancer Invasion and Metastasis: Changing Views. J. Pathol. 2008, 214, 283–293. [Google Scholar] [CrossRef] [PubMed]
  26. Pavón, M.A.; Arroyo-Solera, I.; Céspedes, M.V.; Casanova, I.; León, X.; Mangues, R. UPA/UPAR and SERPINE1 in Head and Neck Cancer: Role in Tumor Resistance, Metastasis, Prognosis and Therapy. Oncotarget 2016, 7, 57351–57366. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Sun, Z.-J.; Yu, G.-T.; Huang, C.-F.; Bu, L.-L.; Liu, J.-F.; Ma, S.-R.; Zhang, W.-F.; Liu, B.; Zhang, L. Hypoxia Induces TFE3 Expression in Head and Neck Squamous Cell Carcinoma. Oncotarget 2016, 7, 11651–11663. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Oliveira, S.S.; Morgado-Díaz, J.A. Claudins: Multifunctional Players in Epithelial Tight Junctions and Their Role in Cancer. Cell. Mol. Life Sci. 2007, 64, 17–28. [Google Scholar] [CrossRef] [PubMed]
  29. Osanai, M.; Takasawa, A.; Murata, M.; Sawada, N. Claudins in Cancer: Bench to Bedside. Pflüg. Arch.—Eur. J. Physiol. 2017, 469, 55–67. [Google Scholar] [CrossRef] [PubMed]
  30. Lamouille, S.; Xu, J.; Derynck, R. Molecular Mechanisms of Epithelial–Mesenchymal Transition. Nat. Rev. Mol. Cell Biol. 2014, 15, 178–196. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Hu, Z.I.; Ghafoor, A.; Sengupta, M.; Hassan, R. Malignant Mesothelioma: Advances in Immune Checkpoint Inhibitor and Mesothelin-targeted Therapies. Cancer 2021, 127, 1010–1020. [Google Scholar] [CrossRef] [PubMed]
  32. Pastan, I.; Hassan, R. Discovery of Mesothelin and Exploiting It as a Target for Immunotherapy. Cancer Res. 2014, 74, 2907–2912. [Google Scholar] [CrossRef] [Green Version]
  33. Creaney, J.; Robinson, B.W.S. Malignant Mesothelioma Biomarkers. Chest 2017, 152, 143–149. [Google Scholar] [CrossRef] [Green Version]
  34. Hassan, R.; Thomas, A.; Alewine, C.; Le, D.T.; Jaffee, E.M.; Pastan, I. Mesothelin Immunotherapy for Cancer: Ready for Prime Time? J. Clin. Oncol. 2016, 34, 4171–4179. [Google Scholar] [CrossRef] [Green Version]
  35. Shimazaki, R.; Takano, S.; Satoh, M.; Takada, M.; Miyahara, Y.; Sasaki, K.; Yoshitomi, H.; Kagawa, S.; Furukawa, K.; Takayashiki, T.; et al. Complement Factor B Regulates Cellular Senescence and Is Associated with Poor Prognosis in Pancreatic Cancer. Cell. Oncol. 2021, 44, 937–950. [Google Scholar] [CrossRef]
  36. He, C.; Li, Y.; Zhang, R.; Chen, J.; Feng, X.; Duan, Y. Low CFB Expression Is Independently Associated with Poor Overall and Disease-free Survival in Patients with Lung Adenocarcinoma. Oncol. Lett. 2021, 21, 478. [Google Scholar] [CrossRef]
  37. Callow, M.G.; Clairvoyant, F.; Zhu, S.; Schryver, B.; Whyte, D.B.; Bischoff, J.R.; Jallal, B.; Smeal, T. Requirement for PAK4 in the Anchorage-Independent Growth of Human Cancer Cell Lines. J. Biol. Chem. 2002, 277, 550–558. [Google Scholar] [CrossRef] [Green Version]
  38. Abril-Rodriguez, G.; Torrejon, D.Y.; Liu, W.; Zaretsky, J.M.; Nowicki, T.S.; Tsoi, J.; Puig-Saus, C.; Baselga-Carretero, I.; Medina, E.; Quist, M.J.; et al. PAK4 Inhibition Improves PD-1 Blockade Immunotherapy. Nat. Cancer 2020, 1, 46–58. [Google Scholar] [CrossRef] [Green Version]
  39. Naïja, A.; Merhi, M.; Inchakalody, V.; Fernandes, Q.; Mestiri, S.; Prabhu, K.S.; Uddin, S.; Dermime, S. The Role of PAK4 in the Immune System and Its Potential Implication in Cancer Immunotherapy. Cell. Immunol. 2021, 367, 104408. [Google Scholar] [CrossRef]
  40. Ryu, D.; Lee, C. Expression Quantitative Trait Loci for PI3K/AKT Pathway. Medicine 2017, 96, e5817. [Google Scholar] [CrossRef]
  41. Zhang, C.; Liu, S.; Wang, X.; Liu, H.; Zhou, X.; Liu, H. COL1A1 Is a Potential Prognostic Biomarker and Correlated with Immune Infiltration in Mesothelioma. BioMed Res. Int. 2021, 2021, 5320941. [Google Scholar] [CrossRef]
  42. Abayasiriwardana, K.S.; Wood, M.K.; Prêle, C.M.; Birnie, K.A.; Robinson, B.W.; Laurent, G.J.; McAnulty, R.J.; Mutsaers, S.E. Inhibition of Collagen Production Delays Malignant Mesothelioma Tumor Growth in a Murine Model. Biochem. Biophys. Res. Commun. 2019, 510, 198–204. [Google Scholar] [CrossRef] [PubMed]
  43. Chen, R.; Zeng, L.; Zhu, S.; Liu, J.; Zeh, H.J.; Kroemer, G.; Wang, H.; Billiar, T.R.; Jiang, J.; Tang, D.; et al. CAMP Metabolism Controls Caspase-11 Inflammasome Activation and Pyroptosis in Sepsis. Sci. Adv. 2019, 5, eaav5562. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Yu, Y.; Tian, X. Analysis of Genes Associated with Prognosis of Lung Adenocarcinoma Based on GEO and TCGA Databases. Medicine 2020, 99, e20183. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Histological images of diffuse pleural mesothelioma. (A) Epithelioid mesothelioma characterized by a tubular and solid pattern; (B) biphasic mesothelioma showing both epithelioid and sarcomatoid malignant areas; (C) sarcomatoid mesothelioma characterized by malignant spindle cells within a fibrous stroma. (Hematoxylin and Eosin staining, magnification 10×).
Figure 1. Histological images of diffuse pleural mesothelioma. (A) Epithelioid mesothelioma characterized by a tubular and solid pattern; (B) biphasic mesothelioma showing both epithelioid and sarcomatoid malignant areas; (C) sarcomatoid mesothelioma characterized by malignant spindle cells within a fibrous stroma. (Hematoxylin and Eosin staining, magnification 10×).
Diagnostics 12 00674 g001
Figure 2. Principal component analysis. Epithelioid (blue) and sarcomatoid (magenta) PM are clearly separated, biphasic PM seems more heterogeneous and overlaps with the two other histotypes, especially epitheliod. MB, biphasic pleural mesothelioma; ME, epithelioid pleural mesothelioma; MM, sarcomatoid pleural mesothelioma.
Figure 2. Principal component analysis. Epithelioid (blue) and sarcomatoid (magenta) PM are clearly separated, biphasic PM seems more heterogeneous and overlaps with the two other histotypes, especially epitheliod. MB, biphasic pleural mesothelioma; ME, epithelioid pleural mesothelioma; MM, sarcomatoid pleural mesothelioma.
Diagnostics 12 00674 g002
Figure 3. Venn diagram. Overlapping of differentially expressed genes among epithelioid, biphasic and sarcomatoid PM. Color codes identify subgroups of deregulated genes. In bold are genes, consistently upregulated in histotypes with higher sarcomatoid component, while plain font is for genes consistently upregulated in histotypes with higher epithelioid component. Asterisks identify genes that are highly expressed in biphasic compared to both epithelioid and sarcomatoid tumors. B, biphasic pleural mesothelioma; E, epithelioid pleural mesothelioma; S, sarcomatoid pleural mesothelioma.
Figure 3. Venn diagram. Overlapping of differentially expressed genes among epithelioid, biphasic and sarcomatoid PM. Color codes identify subgroups of deregulated genes. In bold are genes, consistently upregulated in histotypes with higher sarcomatoid component, while plain font is for genes consistently upregulated in histotypes with higher epithelioid component. Asterisks identify genes that are highly expressed in biphasic compared to both epithelioid and sarcomatoid tumors. B, biphasic pleural mesothelioma; E, epithelioid pleural mesothelioma; S, sarcomatoid pleural mesothelioma.
Diagnostics 12 00674 g003
Figure 4. Violin plots. Genes deregulated among the three PM subtypes. * FDR < 0.05; ** FDR < 0.01; *** FDR < 0.001.
Figure 4. Violin plots. Genes deregulated among the three PM subtypes. * FDR < 0.05; ** FDR < 0.01; *** FDR < 0.001.
Diagnostics 12 00674 g004
Figure 5. Violin plot of epithelial component according to biphasic groups identified by the non-negative matrix factorization algorithm on the basis of gene expression levels. The percentage of epithelioid component was separately and independently determined and it was statistically different among the two identified clusters. Median percentage of epithelioid component in group 1: 40% (interquartile range: 20–52.5%). Median percentage of epithelioid component in group 2: 55% (interquartile range: 50–72.5%).
Figure 5. Violin plot of epithelial component according to biphasic groups identified by the non-negative matrix factorization algorithm on the basis of gene expression levels. The percentage of epithelioid component was separately and independently determined and it was statistically different among the two identified clusters. Median percentage of epithelioid component in group 1: 40% (interquartile range: 20–52.5%). Median percentage of epithelioid component in group 2: 55% (interquartile range: 50–72.5%).
Diagnostics 12 00674 g005
Figure 6. Overall survival curves of mesothelioma cases from TCGA cohort.
Figure 6. Overall survival curves of mesothelioma cases from TCGA cohort.
Diagnostics 12 00674 g006
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bruno, R.; Poma, A.M.; Alì, G.; Distefano, C.; Proietti, A.; Chella, A.; Lucchi, M.; Melfi, F.; Franco, R.; Fontanini, G. Gene Expression Analysis of Biphasic Pleural Mesothelioma: New Potential Diagnostic and Prognostic Markers. Diagnostics 2022, 12, 674. https://doi.org/10.3390/diagnostics12030674

AMA Style

Bruno R, Poma AM, Alì G, Distefano C, Proietti A, Chella A, Lucchi M, Melfi F, Franco R, Fontanini G. Gene Expression Analysis of Biphasic Pleural Mesothelioma: New Potential Diagnostic and Prognostic Markers. Diagnostics. 2022; 12(3):674. https://doi.org/10.3390/diagnostics12030674

Chicago/Turabian Style

Bruno, Rossella, Anello Marcello Poma, Greta Alì, Claudia Distefano, Agnese Proietti, Antonio Chella, Marco Lucchi, Franca Melfi, Renato Franco, and Gabriella Fontanini. 2022. "Gene Expression Analysis of Biphasic Pleural Mesothelioma: New Potential Diagnostic and Prognostic Markers" Diagnostics 12, no. 3: 674. https://doi.org/10.3390/diagnostics12030674

APA Style

Bruno, R., Poma, A. M., Alì, G., Distefano, C., Proietti, A., Chella, A., Lucchi, M., Melfi, F., Franco, R., & Fontanini, G. (2022). Gene Expression Analysis of Biphasic Pleural Mesothelioma: New Potential Diagnostic and Prognostic Markers. Diagnostics, 12(3), 674. https://doi.org/10.3390/diagnostics12030674

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop