Next Article in Journal
Resection Margins in Head and Neck Cancer Surgery: An Update of Residual Disease and Field Cancerization
Next Article in Special Issue
Ancillary Diagnostic Investigations in Malignant Pleural Mesothelioma
Previous Article in Journal
The Migration and Invasion of Oral Squamous Carcinoma Cells: Matrix, Growth Factor and Signalling Involvement
Previous Article in Special Issue
Pathological Characterization of Tumor Immune Microenvironment (TIME) in Malignant Pleural Mesothelioma
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

New DNA Methylation Signals for Malignant Pleural Mesothelioma Risk Assessment

1
Department of Medical Sciences, University of Turin, 10126 Turin, Italy
2
Italian Institute for Genomic Medicine, IIGM, 10060 Candiolo, Italy
3
Candiolo Cancer Institute, FPO-IRCCS, 10060 Candiolo, Italy
4
Mesothelioma Unit, Azienda Ospedaliera SS. Antonio e Biagio e Cesare Arrigo, 15121 Alessandria, Italy
5
Medical Statistics, Department of Translational Medicine, University of Eastern Piedmont, 28100 Novara, Italy
6
Cancer Epidemiology Unit, CPO-Piemonte, 28100 Novara, Italy
7
Department of Health Sciences, University of Eastern Piedmont, 28100 Novara, Italy
8
IRCCS Istituto delle Scienze Neurologiche di Bologna, 40126 Bologna, Italy
9
Department of Integrated Activities Research and Innovation–Azienda Ospedaliera SS. Antonio e Biagio e Cesare Arrigo, 15122 Alessandria, Italy
10
Cancer Epidemiology Unit, Department of Medical Sciences, University of Turin, 10126 Turin, Italy
11
Interdepartmental Center for Studies on Asbestos and Other Toxic Particulates “G. Scansetti”, University of Turin, 10126 Turin, Italy
12
Medical Genetics Unit, AOU Città della Salute e della Scienza, 10126 Turin, Italy
*
Authors to whom correspondence should be addressed.
Cancers 2021, 13(11), 2636; https://doi.org/10.3390/cancers13112636
Submission received: 30 March 2021 / Revised: 24 May 2021 / Accepted: 24 May 2021 / Published: 27 May 2021
(This article belongs to the Special Issue Malignant Mesothelioma)

Abstract

:

Simple Summary

Our study investigated DNA methylation differences in easily accessible white blood cells (WBCs) between malignant pleural mesothelioma (MPM) cases and asbestos-exposed cancer-free controls. A multiple regression model highlighted that the methylation level of two single CpGs (cg03546163 in FKBP5 and cg06633438 in MLLT1) are independent MPM markers. The epigenetic changes at the FKBP5 and MLLT1 genes were robustly associated with MPM in asbestos-exposed subjects. Interaction analyses showed that MPM cases and cancer-free controls showed DNAm differences which may be linked to asbestos exposure.

Abstract

Malignant pleural mesothelioma (MPM) is a rare and aggressive neoplasm. Patients are usually diagnosed when current treatments have limited benefits, highlighting the need for noninvasive tests aimed at an MPM risk assessment tool that might improve life expectancy. Three hundred asbestos-exposed subjects (163 MPM cases and 137 cancer-free controls), from the same geographical region in Italy, were recruited. The evaluation of asbestos exposure was conducted considering the frequency, the duration and the intensity of occupational, environmental and domestic exposure. A genome-wide methylation array was performed to identify novel blood DNA methylation (DNAm) markers of MPM. Multiple regression analyses adjusting for potential confounding factors and interaction between asbestos exposure and DNAm on the MPM odds ratio were applied. Epigenome-wide analysis (EWAS) revealed 12 single-CpGs associated with the disease. Two of these showed high statistical power (99%) and effect size (>0.05) after false discovery rate (FDR) multiple comparison corrections: (i) cg03546163 in FKBP5, significantly hypomethylated in cases (Mean Difference in beta values (MD) = −0.09, 95% CI = −0.12|−0.06, p = 1.2 × 10−7), and (ii) cg06633438 in MLLT1, statistically hypermethylated in cases (MD = 0.07, 95% CI = 0.04|0.10, p = 1.0 × 10−6). Based on the interaction analysis, asbestos exposure and epigenetic profile together may improve MPM risk assessment. Above-median asbestos exposure and hypomethylation of cg03546163 in FKBP5 (OR = 20.84, 95% CI = 8.71|53.96, p = 5.5 × 10−11) and hypermethylation of cg06633438 in MLLT1 (OR = 11.71, 95% CI = 4.97|29.64, p = 5.9 × 10−8) genes compared to below-median asbestos exposure and hyper/hypomethylation of single-CpG DNAm, respectively. Receiver Operation Characteristics (ROC) for Case-Control Discrimination showed a significant increase in MPM discrimination when DNAm information was added in the model (baseline model, BM: asbestos exposure, age, gender and white blood cells); area under the curve, AUC = 0.75; BM + cg03546163 at FKBP5. AUC = 0.89, 2.1 × 10−7; BM + cg06633438 at MLLT1. AUC = 0.89, 6.3 × 10−8. Validation and replication procedures, considering independent sample size and a different DNAm analysis technique, confirmed the observed associations. Our results suggest the potential application of DNAm profiles in blood to develop noninvasive tests for MPM risk assessment in asbestos-exposed subjects.

1. Introduction

Mesothelioma has a long latency period, usually emerging 20–40 years after asbestos exposure [1]. Malignant pleural mesothelioma (MPM) is rare (prevalence 1–9/100,000), but the corresponding annual death toll worldwide is still estimated at about 40,000 [2,3]. Each year, 125 million people are exposed to asbestos, according to a World Health Organization report [4]. The International Agency for Research on Cancer confirmed that all fibrous forms of asbestos are carcinogenic to humans. The main outcome of exposure is mesothelioma, but cancer at other sites, such as respiratory-tract tumors, are moderately frequent [4]. Previous in vitro studies have demonstrated the cytotoxic effects of asbestos fibers [5,6].
A significant association between MPM and asbestos exposure has been reported, showing a clear, increasing trend in the odds ratio (OR) with increasing cumulative exposure among subjects exposed to over 10 fiber/mL-years [7]. Another study reported that the incidence of malignant mesothelioma (MM) was strongly associated with the proximity of one’s residence to an asbestos exposure source [8].
DNA methylation (DNAm) is an epigenetic mechanism involved in gene expression regulation. In particular, dysregulation of promoter DNAm and histone modification are epigenetic mechanisms involved in human malignancies [9].
According to recent papers, both DNAm and genetic alterations may contribute to MPM tumorigenesis [10,11,12,13,14,15]. Whereas the genome remains consistent throughout one’s lifetime, factors like ageing, lifestyle, environmental exposures and diseases can modify DNAm. The adaptive nature of DNAm means that it can be used to link environmental factors to the development of pathologic phenotypes such as cancers. Fasanelli et al. observed an association between exposure to tobacco and site-specific CpG methylation. They also used peripheral blood DNA to evaluate the importance of these epigenetic alterations in the aetiology of lung cancer [16].
There is less information on the mechanisms and clinical outcomes of epigenetic derangements in MPM [17,18,19]. Several studies have evaluated DNAm alterations in MM samples [20,21,22], but few of them focused on DNAm alteration in blood as a circulating marker. Fischer et al. examined serum DNAm of nine gene-specific promoters from MM cases [23]. A more recent paper identified hypomethylation of a single CpG in FKBP5 in whole blood cells as a predictor of overall survival in MPM cases [13]. Guarrera et al. evaluated methylation levels in DNA from whole blood leukocytes as potential diagnostic markers for MPM and found a differential methylation between asbestos-exposed MPM cases and controls, mainly in genes related to the immune system [11]. The identification of reliable DNAm biomarkers with high sensitivity and specificity for MPM risk assessment would be a major advancement.
This study was undertaken with the goal to identify new biomarkers for MPM risk assessment and to determine if peripheral blood DNAm profiles have any predictive value. The second goal was to evaluate the interaction effect of asbestos exposure with DNAm on MPM risk. Currently, there are no sensitive testing methods available for the screening of asbestos-exposed individuals who are at high risk of developing MPM. Thus, the identification of reliable MPM diagnostic biomarkers in peripheral blood might provide a tool for detecting the disease at an early stage.

2. Results

2.1. Epigenome-Wide Association Study (EWAS)

CpGs (445,254) passed quality control procedures and were considered for statistical analyses. EWAS revealed two statistically significant differentially methylated single-CpGs between case and control groups: cg03546163 in the FKBP5 gene (Mean Difference in beta values (MD) = 0.09, 95% CI = −0.12|−0.06, p = 1.2 × 10−7, p = 0.028) and cg06633438 in the MLLT1 gene (MD = 0.07, 95% CI = 0.04|0.10, p = 1.0 × 10−6, p = 0.049) after False Discovery Rate (FDR) post hoc correction (Figure 1; Table 1).
Another 10 CpGs showed hypo/hypermethylation in MPM considering FDR < 0.05 but not effect size (MD) cut off ≥ |0.05| (Table 1).
Bootstrap was computed to estimate measures of accuracy using random sampling methods. The bias-corrected and accelerated (BCa) bootstrap interval was calculated for cg03546163 in FKBP5 (95% CIBCa = −0.16|−0.10, z0 = −0.008, a = 0.002) and cg06633438 in MLLT1 (95% CIBCa = −0.06|−0.1, z0 = −0.011, a = 0.0004) genes, confirming the robustness of the results considering the sample under study.
Statistically significant differences in MD between cases and controls were found in the WBCs estimated (monocytes, p = 6.0 × 10−3; granulocytes, p = 2.2 × 10−16; B cells, p = 1.1 × 10−12; NK cells, p = 3.6 × 10−4; CD4+ T, p = 2.2 × 10−16; CD8+ T, p = 6.8 × 10−11; Naïve CD4T, p = 0.012; Naïve CD8T, p = 7.0 × 10−3).
In order to assess if smoking status, classified as current, former and never-smokers, could modify DNAm profiles, we performed a multivariate regression analysis with the same model used for the discovery phase. No evidence of methylation differences linked to different smoking levels was found for any of the twelve statistically significant CpGs.

2.2. Receiver Operation Characteristics (ROC) for Case-Control Discrimination

The baseline model (BM) including age, gender, asbestos exposure and WBCs was compared with BM adding the DNAm levels of cg03546163 or cg06633438. Receiver Operation Characteristics 8ROC9 curves showed a significant increase in MPM discrimination when DNAm information was added in the model (Table 2).

2.3. Interaction Analysis

CpG sites and asbestos exposure were considered as predictors of MPM risk in the interaction model. Categorical variables (quantile information) were used considering median values.
We tested the interaction between asbestos exposure and DNAm levels at cg03546163 in FKBP5 and cg06633438 in MLLT1.
Considering cg03546163 in FKBP5, DNA hypermethylation and low asbestos exposure levels were used as references, while for cg06633438 in MLLT1, DNA hypomethylation and low asbestos exposure levels were set as references (Table 3).
The OR was estimated as the relationship between the combination of single-CpGs DNAm levels and asbestos exposure quantile, and the reference (low median asbestos exposure and hypermethylation status for cg03546163, or hypomethylation status for cg06633438). Age, gender, population stratification, and WBCs were included in the GLM (family = binomial) to adjust the interaction effect.
The relationship between asbestos exposures and single-CpG DNAm levels was evaluated. An increase of one unit of asbestos exposure (rank transformed fibers/mL years) was related to the FKBP5 gene (β = −0.016, 95% CI = −0.031|−0.001, p = 0.044) and MLLT1 gene (β = −0.014, 95% CI = 0.001|0.026, p = 0.035) methylation level variations.
Strong association between asbestos exposure and MPM risk, considering dichotomous distribution of asbestos exposure, was found (OR = 6.11, 95% CI = 3.73|10.20, p = 1.8 × 10−12). Quartile distribution of asbestos exposure was evaluated to estimate the potential incremental association with MPM risk (1st quartile: used as reference; 2nd quartile: OR = 1.83, 95% CI = 0.93|3.69, p = 0.09; 3rd quartile: OR = 6.63, 95% CI = 3.30|13.81, p = 2.1 × 10−7; 4th quartile: OR = 11.00, 95% CI = 5.26|24.30, p = 7.3 × 10−10).

2.4. Validation and Replication

For the replication and validation approaches, an independent sample of 140 MPM cases and 104 cancer-free asbestos-exposed controls from the same areas were considered, using a targeted DNAm analysis technique.
The direction and magnitude of the association was consistent for cg03546163 and cg06633438 DNAm. Patients showed significantly lower DNAm at cg03546163 (MD = −0.061, 95% CI = −0.087|−0.036, p = 4.5 × 10−6) and higher DNAm at cg06633438 (MD = 0.024, 95% CI = 0.061|0.013, p = 4.0 × 10−2) compared with controls. A multivariate analysis confirmed that DNAm at cg03546163 in FKBP5 and cg06633438 in MLLT1 were independently associated with MPM detection.

3. Discussion

In the present study, we used a whole genome microarray approach to investigate DNAm in WBCs from MPM cases and asbestos-exposed cancer-free controls from a region with a history of asbestos exposure (Piedmont, Italy) [10] in order to identify new noninvasive epigenetic markers related to MPM. The identification of reliable MPM diagnostic biomarkers in peripheral blood might improve risk assessment.
We observed hypomethylation of CpG cg03546163, located in the 5′ UTR region of FKBP5 gene, in MPM cases compared to controls.
Epigenetic activation of the FKBP Prolyl Isomerase 5 (FKBP5) gene has been shown to be associated with increased stress sensitivity and the risk of psychiatric disorders [24]. FKBP5 is an immunophilin and has an important role in immunoregulation and in protein folding and trafficking. It plays a role in transcriptional complexes and acts as a cotranscription factor, along with other proteins in the FKBP family [25]. The suggestion of a possible role of FKBP5 in the development and progression of different types of cancer has stemmed from several studies. In particular, high protein expression has been linked to either suppression or promotion of tumour growth, depending on tumour type and microenvironment [26,27].
FKBP5 is involved in the NF-kB and AKT signaling pathways, both of which are implicated in tumorigenesis [28]. Notably, NF-kB appears to be frequently constitutively activated in malignant tumours and involved in the modulation of genes linked to cell motility, neoangiogenesis, proliferation and programmed cell death [29]. An epigenetic upregulation of FKBP5 could promote NF-kB activation [30]. STAT3-NFkB activity is involved in chemoresistance in MM cells [31], and NFkB was shown to be constitutively active as a result of asbestos-induced chronic inflammation [32].
CpG cg06633438 located in the body region of the MLLT1 gene was hypermethylated in cases compared to controls.
The MLLT1 gene encodes the ENL protein, a histone acetylation reader component of the super elongation complex (SEC), which promotes transcription at the elongation stage by suppressing transient pausing by the polymerase at multiple sites along the DNA. In acute myeloid leukemia, MLLT1 regulates chromatin remodeling and gene expression of many important proto-oncogenes [31]. Yoshikawa and colleagues suggested that mesothelioma may be the consequence of the somatic inactivation of chromatin-remodeling complexes and/or histone modifiers, including MLLT1 [30].
In mesothelioma patients with short-term recurrence after surgery, frequent 19p13.2 loss was reported. This region encompasses several putative tumor suppressors or oncogenes, including MLLT1 [32].
Interestingly, MLLT1 and FKBP5 showed opposite behavior, increasing and decreasing DNAm levels, respectively, in relation to MPM. This finding could reflect the opposite expression profiles of the two genes among all the different subtypes of white blood cells in normal human hematopoiesis, as reported in the Blood Spot database (http://servers.binf.ku.dk/bloodspot/, accessed on 26 May 2021) (Figure 2) [33].
Our interaction analysis showed that considering DNAm levels at FKBP5 and MLLT1 genes together with asbestos exposure levels may help to better define MPM risk for asbestos-exposed subjects.
Six single-CpGs showed differential methylation in patients, including those located in C13orf30, CDC42BPB, ZNF629 and ALG9 genes; the other six were not annotated to named genes. ALG9 is a glycogene whose reduced expression has been described during the epithelial-to-mesenchymal transition, an essential process also involved in cancer progression [34]. The CDC42BPB gene is ubiquitously expressed in mammals and encodes a serine/threonine protein kinase, a member of the MRCK family [35]. The role of MRCKs in cytoskeletal reorganization during cell migration and invasion has been characterized [36]. The biological function of C13orf30 and ZNF629, a DNA-binding transcription factor, is still to be established.
MPM cases and asbestos-exposed controls showed different proportions of estimated WBCs, which may denote the crucial implication of the immune system. It is known that in cancer, including mesothelioma, the immune system is affected [37], and there is evidence that asbestos directs antigen overstimulation, and that reactive oxygen species production induces functional changes in WBCs [38]. Indeed, in MPM cases, we showed a reduction of estimated CD4+ and CD8+ T lymphocytes, suggesting a weaker adaptive immune system [39]. This may reflect the possible occurrence of functional changes in WBC subtypes in MPM [40,41].
The need for reliable biomarkers is of extreme relevance for a disease such as MPM, which is characterized by the accumulation and persistence of asbestos fibers in the lungs, leading to a long latency period before clear clinical signs of the tumor are detectable. Several biomarkers for early MPM detection (e.g., mesothelin, osteopontin and fibulin-3) have been proposed so far; however, some of them are still under investigation [42]. In this context, DNAm changes in easily-accessible WBCs may provide a useful tool to better assess MPM risk in asbestos-exposed subjects.
Our findings that DNAm levels in single-CpGs in FKBP5 and MLLT1 genes are independent markers of MPM in asbestos-exposed subjects suggest the potential use of blood DNAm analysis as a noninvasive test for MPM detection.
Some somatic gene alterations in lung cancer have been linked to tobacco smoke, but few data are available on the role of asbestos fibers: Andujar and colleagues investigate the mechanism of P16/CDKN2A alterations in lung cancer including asbestos-exposed patients. P16/CDKN2A gene inactivation in asbestos-exposed non-small-cell lung carcinoma (NSCLC) cases, a tumor independent of tobacco smoking but associated with asbestos exposure, mainly occurs via promoter hypermethylation, loss of heterozygosity and homozygous deletion, suggesting a possible relationship with an effect of asbestos fibers [43].We observed epigenetic deregulations in the blood of MPM patients compared to that of cancer-free controls, suggesting the potential use of DNAm for risk stratification among asbestos-exposed individuals.
If this observation can be verified in prospectively collected samples, it may be possible to use CpGs methylation to further improve MPM risk estimation for subjects with occupational and/or environmental asbestos exposure.

Limitation of the Study

Leukocyte DNAm may be a nonspecific marker related to a general, tumour-induced inflammatory status rather than a specific MPM biomarker. Further studies are therefore needed to support our findings.
One main limitation of the functional interpretation of our results is that all our cases had already developed MPM at recruitment: thus, our findings likely reflect disease status rather than being markers of the dynamic processes leading to MPM onset. The lack of MPM tissue from the same subjects also poses major constraints to the functional interpretation of our findings.
Notwithstanding the above limitations, the discrimination between MPM cases and asbestos-exposed cancer-free controls improved when DNAm levels were considered together with asbestos exposure levels.

4. Material and Methods

4.1. Study Population

Study subjects were part of a wider, ongoing collaborative study, which is actively enrolling MPM cases and cancer-free controls in the municipality of Casale Monferrato (Piedmont Region, Italy). This area was chosen due to its exceptionally high incidence of mesothelioma, caused by widespread occupational and environmental asbestos exposure originating from the Eternit asbestos-cement plant, which was operational until 1986 [44]. Additional MPM cases and cancer-free controls were recruited from other main hospitals of the Piedmont Region (in the municipalities of Turin, Novara and Alessandria). The ongoing collaborative study includes MPM cases diagnosed between incident MPM cases diagnosed between 2000 and 2010 after histological and/or cytological confirmation, and matched controls [45].
The present study included 159 MPM cases and 137 cancer-free controls from a larger case-control study, all of whom had genetic and blood DNAm data [46], good quality DNA at the time of the analyses, and information on asbestos exposure above the background level, as defined in Ferrante et al. [47]. MPM cases and asbestos-exposed cancer-free controls were matched by date of birth (±18 months) and gender. An additional 244 (140 MPM cases and 104 cancer-free controls) independent samples from the same case-control study were included for validation/replication analyses.
Table 4 and Table 5 shows the descriptive characteristics of controls and cases (Min, 1st Q, Median, Mean, 3rd Q and Max) that were considered in the statistical analysis (gender, age, asbestos exposure and WBC estimates: monocytes, granulocytes, natural killer, B cells, CD4+ T and CD8+ T). Asbestos exposure (occupational, environmental and domestic) was normalized considering frequency, duration and intensity. Smoking status (current, former and never smokers) is also explained in Table 6.
Our study complied with the Declaration of Helsinki principles and conformed to ethical requirements. All volunteers signed an informed consent form at enrollment. The study protocol was approved by the Ethics Committee of the Italian Institute for Genomic Medicine (IIGM, Candiolo, Italy).

4.2. Exposure Assessment

Information on occupational history and lifestyle habits were collected from all subjects through interviewer-administered questionnaires, which were completed during face-to-face interviews at enrollment. Job titles were coded in two ways according to the International Standard Classification of Occupations [47] and the Statistical Classification of Economic Activities in the European Community.
A cumulative exposure index was computed considering frequency, duration and intensity of asbestos exposure. Occupational, environmental and domestic asbestos exposure were evaluated by an experienced occupational epidemiologist [47], and exposure doses (fibers/mL years) were rank-transformed to remove skewness.

4.3. Blood DNAm Analysis and Beta-Value Extraction

DNAm levels were measured in DNA from whole blood collected at enrollment using the Infinium HumanMethylation450 BeadChip (Illumina, San Diego, CA, USA). For blood DNAm analysis (including quality control) please refer to the previous work of the same group [11].
We used the R package ‘methylumi’ to analyze DNAm data. The average methylation value at each locus was computed as the ratio of the intensity of the methylated signal over the total signal (unmethylated + methylated) [48]. Beta-values ranging from 0 (no methylation) to 1 (full methylation) represent the percentage of methylation at each individual CpG locus.
We excluded the following from the analyses: (i) single beta-values with a p-value for detection ≥ 0.01; (ii) CpG loci that had missing beta-values in more than 20% of the assayed samples; (iii) CpG loci detected by probes containing single nucleotide polymorphisms (SNPs) with MAF ≥ 0.05 in the CEPH (Utah residents with ancestry from northern and western Europe, CEU) population; and iv) samples with a global call rate ≤ 95%. We also excluded CpGs on chromosomes X and Y.

4.4. Batch Effect, Population Stratification and White Blood Cells Estimations

All differential methylation analyses were corrected for “control probes” Principal Components (PCs) to account for variability and batch effects in methylation assays. We used PCs assessed by principal component analysis of the BeadChip’s built-in control probes as a correction factor for statistical analyses of microarray data. This method allows researchers to account for the technical variability in the different steps in DNAm analysis, from bisulfite conversion to BeadChip processing [49].
An individual’s geographic origins may influence DNAm profiles, which could potentially introduce bias. To take this into consideration, we took advantage of the available data from our previous study, which includes a genome-wide genotyping dataset from the same study subjects [50]. When genome-wide genotyping was used to calculate the first PCs, they were shown to correlate with different geographic origins [51].
For each subject, we extracted WBC subtype percentages, estimated based on genome-wide methylation data. This method provides quantification of the composition of leukocytes than can be achieved by simple histological or flow cytometric assessments, with an admissible range of variability [52].

4.5. Statistical Analyses

4.5.1. Epigenome-Wide Association Study

An association test was used to analyze the mean differences (MD) in single-CpG methylation between MPM cases and asbestos-exposed cancer-free controls. We performed multiple regression analysis adjusted for age, gender, estimated WBCs (monocytes, granulocytes, natural killer, B cells, CD4+ T and CD8+ T), population stratification (first 2 PCs) and technical variability (first 10 PCs). For multiple comparison tests, a FDR p value ≤ 0.05 was considered statistically significant.
Bootstrapping was performed using random sampling methods to estimate the measures of accuracy defined in terms of bias, variance, confidence intervals and prediction error. Bootstrapping can also be applied to control and check the results for stability. The bias-corrected and accelerated (BCa) bootstrap interval was calculated with regard to single CpGs.
ROC for Case-Control Discrimination was implemented, and the AUC metric was applied to estimate the predictive performance of a binary classification (cases/controls). The baseline model (BM) included age, gender, asbestos exposure and WBCs, and was compared with the BM after adding the DNAm levels of statistically significant, single-CpGs at EWAS. AUC differences between BMs before and after the addition of DNAm levels were estimated with DeLong’s test.

4.5.2. Statistical Power

To ensure a study power greater than 99% (two-tailed test at α = 0.05 and β = 0.01), only CpGs with a MD between cases and controls ≥ |0.05| were selected.
Covariates were included step-by-step in a sensitivity analysis to validate the association output considering effect size, standard error, 95% confidence interval and p value variations.
Gene set enrichment analyses were carried out on CpGs with a False Discovery Rate p value (PFDR) ≤ 0.05 to identify pathways that may be affected by MPM-related changes in methylation.
All statistical analyses were conducted using the open source software R (4.0.2).

4.5.3. Interaction Analysis

Logistic regression was used to analyze the relationship between CpGs and asbestos exposure in MPM risk (odds ratio), adjusting for age, gender, SNP PCs and WBCs estimates. Asbestos exposure was classified as above-median or below-median, and CpG methylation was categorized as above-median or below-median.
MPM risk for a given CpG level and asbestos exposure was expressed by ORij, where i indicates the asbestos exposure (below-median or above-median) and j indicates the CpG (above-median or below-median). Considering the direction of the effect, the same approach was used: for hypomethylated CpGs, above-median was used as the reference level, whereas below-median was used for hypermethylated CpGs.
Subjects with below-median asbestos exposure and reference-level CpG DNAm were considered the baseline group, and their MPM risk was coded as OR00 = 1. Interaction was analyzed with respect to both additive and multiplicative models based on the ORs obtained by logistic regression.
Synergistic interaction (positive interaction) implies that the combined action of two factors in an additive model is greater than the sum of their individual effects. Antagonistic interaction, on the other hand, means that when two factors are present in an additive model, the action of one reduces the effect of the other.
Multivariable logistic regression models were used to explore any deviations from a multiplicative model, including asbestos exposure, CpG and the corresponding interaction term (CpG × exposure). All models were adjusted for age, gender, SNP PCs, technical covariates and WBCs estimates. p-values < 0.05 were considered statistically significant.

4.6. Validation and Replication

DNAm signal validation and replication was done by the EpiTYPER MassARRAY assay (Agena Bioscience). This assay uses a MALDI-TOF mass spectrometry-based method to quantitatively assess the DNA methylation state of the CpG sites of interest [53]. DNA (500 ng) was bisulfite-converted as indicated in Section 4.3.
PCR amplification, treatment with SAP solution and Transcription/RNase A cocktails were performed according to the manufacturer’s instructions, and the mass spectra were analyzed by an EpiTYPER analyzer. The MassARRAY assay cannot discriminate between CpGs located in close proximity in the sequence, so instead, the close neighboring CpGs are analyzed as “Units”, i.e., the measured methylation level is the average of the methylation levels of the CpGs cumulatively analyzed within the Unit. In the case of cg03546163, the measured methylation level is the average between two CpG sites located in very close proximity (Figure S1). For cg06633438, the two adjacent signals were considered, since the results for the model did not differ for effect size, standard error, 95% CI or p value (Figure S2).

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/cancers13112636/s1, Figure S1: Location of cg03546163 investigated by EpiTYPER MassARRAY, Figure S2: Location of cg06633438 investigated by EpiTYPER MassARRAY.

Author Contributions

Conceptualization, G.C., G.M., C.M., D.M., I.D. and F.G.; methodology, G.C.; software, G.C.; validation, G.C., A.A. and C.P.; formal analysis, G.C.; investigation, G.C., G.M., S.G., A.A., A.R., E.C. and C.C.; resources, G.C., G.M., C.M., D.M., I.D., F.G., D.F., M.S., M.L.V. and R.L.; data curation, G.C., C.M., D.M., I.D., F.G., C.P., D.F., M.S., M.L.V. and R.L.; writing—original draft preparation, G.C., S.G., A.A., A.R., E.C. and C.C.; writing—review and editing, G.C., G.M., C.M., D.M., I.D., F.G., S.G., A.A., A.R., E.C., C.C., C.P., D.F., M.S., M.L.V. and R.L.; visualization, G.C., G.M., S.G., A.A., A.R., E.C., C.C., C.P., D.F., M.S., M.L.V. and R.L.; supervision, G.C. and G.M.; project administration, G.C. and G.M.; funding acquisition, G.M., C.M., D.M., I.D. and F.G. All authors have read and agreed to the published version of the manuscript.

Funding

The research leading to these results received funding from AIRC under IG 2018-ID. 21390 project (to G.M.), from the HERMES (Hereditary Risk in MESothelioma) Project funded by the offer of compensation to the inhabitants of Casale Monferrato deceased or affected by mesothelioma (to I.D. and C.M.), and from the Ministero dell’Istruzione, dell’Università e della Ricerca–MIUR project “Dipartimenti di Eccellenza 2018–2022” (n° D15D18000410001, to G.M.) to the Department of Medical Sciences, University of Torino.

Institutional Review Board Statement

Ethics approval and consent to participate: The study protocol was approved by the Ethics Committee of the Italian Institute for Genomic Medicine (IIGM, Candiolo, Italy).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The methylation of single individuals cannot be published due to informed consent limitations.

Acknowledgments

The authors wish to thank all the patients contributing to this research, Caterina Casadio and Ezio Piccolini for their past clinical contribution and Paolo Garagnani for the Sequenom analysis supervision.

Conflicts of Interest

The authors declare that they have no competing interest.

References

  1. Sekido, Y. Molecular pathogenesis of malignant mesothelioma. Carcinogenesis 2013, 34, 1413–1419. [Google Scholar] [CrossRef]
  2. Rossini, M.; Rizzo, P.; Bononi, I.; Clementz, A.; Ferrari, R.; Martini, F.; Tognon, M.G. New Perspectives on Diagnosis and Therapy of Malignant Pleural Mesothelioma. Front. Oncol. 2018, 8, 91. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Furuya, S.; Chimed-Ochir, O.; Takahashi, K.; David, A.; Takala, J. Global Asbestos Disaster. Int. J. Environ. Res. Public Health 2018, 15, 1000. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Straif, K.; Benbrahim-Tallaa, L.; Baan, R.; Grosse, Y.; Secretan, B.; El Ghissassi, F.; Bouvard, V.; Guha, N.; Freeman, C.; Galichet, L.; et al. A review of human carcinogens—Part C: Metals, arsenic, dusts, and fibres. Lancet Oncol. 2009, 10, 453–454. [Google Scholar] [CrossRef]
  5. Jaurand, M.C. Mechanisms of fiber-induced genotoxicity. Environ. Health Perspect. 1997, 105, 1073–1084. [Google Scholar]
  6. Kelsey, K.T.; Yano, E.; Liber, H.L.; Little, J.B. The in vitro genetic effects of fibrous erionite and crocidolite asbestos. Br. J. Cancer 1986, 54, 107–114. [Google Scholar] [CrossRef] [Green Version]
  7. Iwatsubo, Y.; Pairon, J.C.; Boutin, C.; Ménard, O.; Massin, N.; Caillaud, D.; Orlowski, E.; Galateau-Salle, F.; Bignon, J.; Brochard, P. Pleural mesothelioma: Dose-response relation at low levels of asbestos exposure in a French population-based case-control study. Am. J. Epidemiol. 1998, 148, 133–142. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Howel, D.; Arblaster, L.; Swinburne, L.; Schweiger, M.; Renvoize, E.; Hatton, P. Routes of asbestos exposure and the development of mesothelioma in an English region. Occup. Environ. Med. 1997, 54, 403–409. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Kanherkar, R.R.; Bhatia-Dey, N.; Csoka, A.B. Epigenetics across the human lifespan. Front. Cell Dev. Biol. 2014, 2, 49. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Ferrante, D.; Bertolotti, M.; Todesco, A.; Mirabelli, D.; Terracini, B.; Magnani, C. Cancer mortality and incidence of mesothelioma in a cohort of wives of asbestos workers in Casale Monferrato, Italy. Environ. Health Perspect. 2007, 115, 1401–1405. [Google Scholar] [CrossRef]
  11. Guarrera, S.; Viberti, C.; Cugliari, G.; Allione, A.; Casalone, E.; Betti, M.; Ferrante, D.; Aspesi, A.; Casadio, C.; Grosso, F.; et al. Peripheral Blood DNA Methylation as Potential Biomarker of Malignant Pleural Mesothelioma in Asbestos-Exposed Subjects. J. Thorac. Oncol. 2019, 14, 527–539. [Google Scholar] [CrossRef] [Green Version]
  12. Matullo, G.; Guarrera, S.; Betti, M.; Fiorito, G.; Ferrante, D.; Voglino, F.; Cadby, G.; Di Gaetano, C.; Rosa, F.; Russo, A.; et al. Genetic variants associated with increased risk of malignant pleural mesothelioma: A genome-wide association study. PLoS ONE 2013, 8, e61253. [Google Scholar] [CrossRef] [Green Version]
  13. Cugliari, G.; Catalano, C.; Guarrera, S.; Allione, A.; Casalone, E.; Russo, A.; Grosso, F.; Ferrante, D.; Viberti, C.; Aspesi, A.; et al. DNA Methylation of FKBP5 as Predictor of Overall Survival in Malignant Pleural Mesothelioma. Cancers 2020, 12, 3470. [Google Scholar] [CrossRef]
  14. Betti, M.; Aspesi, A.; Sculco, M.; Matullo, G.; Magnani, C.; Dianzani, I. Genetic predisposition for malignant mesothelioma: A concise review. Mutat. Res. 2019, 781, 1–10. [Google Scholar] [CrossRef]
  15. Guo, G.; Chmielecki, J.; Goparaju, C.; Heguy, A.; Dolgalev, I.; Carbone, M.; Seepo, S.; Meyerson, M.; Pass, H.I. Whole-exome sequencing reveals frequent genetic alterations in BAP1, NF2, CDKN2A, and CUL1 in malignant pleural mesothelioma. Cancer Res. 2015, 75, 264–269. [Google Scholar] [CrossRef] [Green Version]
  16. Fasanelli, F.; Baglietto, L.; Ponzi, E.; Guida, F.; Campanella, G.; Johansson, M.; Grankvist, K.; Johansson, M.; Assumma, M.B.; Naccarati, A.; et al. Hypomethylation of smoking-related genes is associated with future lung cancer in four prospective cohorts. Nat. Commun. 2015, 15, 10192. [Google Scholar] [CrossRef] [PubMed]
  17. Moore, L.D.; Le, T.; Fan, G. DNA methylation and its basic function. Neuropsychopharmacology 2013, 38, 23–38. [Google Scholar]
  18. Bononi, A.; Napolitano, A.; Pass, H.I.; Yang, H.; Carbone, M. Latest developments in our understanding of the pathogenesis of mesothelioma and the design of targeted therapies. Expert. Rev. Respir. Med. 2015, 9, 633–654. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Vandermeers, F.; Neelature, S.; Costa, C.; Hubaux, R.; Cosse, J.P.; Willems, L. The role of epigenetics in malignant pleural mesothelioma. Lung Cancer 2013, 81, 311–318. [Google Scholar] [CrossRef] [PubMed]
  20. Zhang, X.; Tang, N.; Rishi, A.K.; Pass, H.I.; Wali, A. Methylation profile landscape in mesothelioma: Possible implications in early detection, disease progression, and therapeutic options. Methods Mol. Biol. 2015, 1238, 235–247. [Google Scholar] [PubMed]
  21. Goto, Y.; Shinjo, K.; Kondo, Y.; Shen, L.; Toyota, M.; Suzuki, H.; Gao, W.; An, N.; Fujii, M.; Murakami, H.; et al. Epigenetic profiles distinguish malignant pleural mesothelioma from lung adenocarcinoma. Cancer Res. 2009, 69, 9073–9082. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Christensen, B.C.; Houseman, E.A.; Poage, G.M.; Godleski, J.J.; Bueno, R.; Sugarbaker, D.J.; Wiencke, J.K.; Nelson, H.H.; Marsit, C.J.; Kelsey, K.T. Integrated profiling reveals a global correlation between epigenetic and genetic alterations in mesothelioma. Cancer Res. 2010, 70, 5686–5694. [Google Scholar] [CrossRef] [Green Version]
  23. Fischer, J.R.; Ohnmacht, U.; Rieger, N.; Zemaitis, M.; Stoffregen, C.; Kostrzewa, M.; Buchholz, E.; Manegold, C.; Lahm, H. Promoter methylation of RASSF1A, RARbeta and DAPK predict poor prognosis of patients with malignant mesothelioma. Lung Cancer 2006, 54, 109–116. [Google Scholar] [CrossRef]
  24. Matosin, N.; Halldorsdottir, T.; Binder, E.B. Understanding the molecular mechanisms underpinning gene by environment interactions in psychiatric disorders: The FKBP5 model. Biol. Psychiatry 2018, 83, 821–830. [Google Scholar] [CrossRef] [PubMed]
  25. Kang, C.B.; Hong, Y.; Dhe-Paganon, S.; Sup Yoon, H. FKBP family proteins: Immunophilins with versatile biological functions. Neurosignals 2008, 16, 318–325. [Google Scholar] [CrossRef]
  26. Li, L.; Lou, Z.; Wang, L. The role of FKBP5 in cancer aetiology and chemoresistance. Br. J. Cancer 2011, 104, 19–23. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Romano, S.; D’Angelillo, A.; Romano, M.F. Pleiotropic roles in cancer biology for multifaceted proteins FKBPs. Biochim. Biophys. Acta 2015, 1850, 2061–2068. [Google Scholar] [CrossRef]
  28. Staibano, S.; Mascolo, M.; Ilardi, G.; Siano, M.; De Rosa, G. Immunohistochemical analysis of FKBP51 in human cancers. Curr. Opin. Pharmacol. 2011, 11, 338–347. [Google Scholar] [CrossRef] [PubMed]
  29. Zannas, A.S.; Jia, M.; Hafner, K.; Baumert, J.; Wiechmann, T.; Pape, J.C.; Arloth, J.; Ködel, M.; Martinelli, S.; Roitman, M.; et al. Epigenetic upregulation of FKBP5 by aging and stress contributes to NF-kappaB-driven inflammation and cardiovascular risk. Proc. Natl. Acad. Sci. USA 2019, 116, 11370–11379. [Google Scholar] [CrossRef] [Green Version]
  30. Yoshikawa, Y.; Sato, A. Biallelic germline and somatic mutations in malignant mesothelioma: Multiple mutations in transcription regulators including mSWI/SNF genes. Int. J. Cancer 2015, 136, 560–571. [Google Scholar] [CrossRef]
  31. Zhou, J.; Ng, Y.; Chng, W.J. ENL: Structure, function, and roles in hematopoiesis and acute myeloid leukemia. Cell Mol. Life Sci. 2018, 75, 3931–3941. [Google Scholar] [CrossRef] [PubMed]
  32. Canino, C.; Luo, Y.; Marcato, P.; Blandino, G.; Pass, H.I.; Cioce, M. STAT3-NFkB/DDIT3/CEBPbeta axis modulates ALDH1A3 expression in chemoresistant cell subpopulations. Oncotarget 2015, 6, 12637–12653. [Google Scholar] [CrossRef] [Green Version]
  33. Bagger, F.O.; Kinalis, S.; Rapin, N. BloodSpot: A database of healthy and malignant haematopoiesis updated with purified and single cell mRNA sequencing profiles. Nucleic Acids Res. 2019, 47, 881–885. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Schins, R.P.; Donaldson, K. Nuclear Factor Kappa-B Activation by Particles and Fibers. Inhal. Toxicol. 2000, 12, 317–326. [Google Scholar] [CrossRef]
  35. Tan, Z.; Lu, W.; Li, X.; Yang, G.; Guo, J.; Yu, H.; Li, Z.; Guan, F. Altered N-Glycan expression profile in epithelial-to-mesenchymal transition of NMuMG cells revealed by an integrated strategy using mass spectrometry and glycogene and lectin microarray analysis. J. Proteome Res. 2014, 13, 2783–2795. [Google Scholar] [CrossRef] [PubMed]
  36. Kale, V.P.; Hengst, J.A.; Desai, D.H.; Dick, T.E.; Choe, K.N.; Colledge, A.L.; Takahashi, Y.; Sung, S.S.; Amin, S.G.; Yun, J.K. A novel selective multikinase inhibitor of ROCK and MRCK effectively blocks cancer cell migration and invasion. Cancer Lett. 2014, 354, 299–310. [Google Scholar] [CrossRef] [Green Version]
  37. Tanrikulu, A.C.; Abakay, A.; Komek, H.; Abakay, O. Prognostic value of the lymphocyte-to-monocyte ratio and other inflammatory markers in malignant pleural mesothelioma. Environ. Health Prev. Med. 2016, 21, 304–311. [Google Scholar] [CrossRef] [Green Version]
  38. Zhao, Z.S.; Manser, E. PAK and other Rho-associated kinases—Effectors with surprisingly diverse mechanisms of regulation. Biochem. J. 2005, 386 Pt 2, 201–214. [Google Scholar] [CrossRef]
  39. Nishimura, Y.; Kumagai-Takei, N.; Matsuzaki, H.; Lee, S.; Maeda, M.; Kishimoto, T.; Fukuoka, K.; Nakano, T.; Otsuki, T. Functional alteration of natural killer cells and cytotoxic T lymphocytes upon asbestos exposure and in malignant mesothelioma patients. Biomed. Res. Int. 2015, 2015, 238431. [Google Scholar] [CrossRef]
  40. Maeda, M.; Nishimura, Y.; Kumagai, N.; Hayashi, H.; Hatayama, T.; Katoh, M.; Miyahara, N.; Yamamoto, S.; Hirastuka, J.; Otsuki, T. Dysregulation of the immune system caused by silica and asbestos. J. Immunotoxicol. 2010, 7, 268–278. [Google Scholar] [CrossRef]
  41. Miura, Y.; Nishimura, Y.; Katsuyama, H.; Maeda, M.; Hayashi, H.; Dong, M.; Hyodoh, F.; Tomita, M.; Matsuo, Y.; Uesaka, A.; et al. Involvement of IL-10 and Bcl-2 in resistance against an asbestos-induced apoptosis of T cells. Apoptosis 2006, 11, 1825–1835. [Google Scholar] [CrossRef]
  42. Cristaudo, A.; Bonotti, A.; Guglielmi, G.; Fallahi, P.; Foddis, R. Serum mesothelin and other biomarkers: What have we learned in the last decade? J. Thorac. Dis. 2018, 10, 353–359. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Andujar, P.; Wang, J.; Descatha, A.; Galateau-Sallé, F.; Abd-Alsamad, I.; Billon-Galland, M.A.; Blons, H.; Clin, B.; Danel, C.; Housset, B.; et al. p16INK4A inactivation mechanisms in non-small-cell lung cancer patients occupationally exposed to asbestos. Lung Cancer 2010, 67, 23–30. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Dianzani, I.; Gibello, L.; Biava, A.; Giordano, M.; Bertolotti, M.; Betti, M.; Ferrante, D.; Guarrera, S.; Betta, G.P.; Mirabelli, D.; et al. Polymorphisms in DNA repair genes as risk factors for asbestos-related malignant mesothelioma in a general population study. Mutat. Res. 2006, 599, 124–134. [Google Scholar] [CrossRef]
  45. Betti, M.; Ferrante, D.; Padoan, M.; Guarrera, S.; Giordano, M.; Aspesi, A.; Mirabelli, D.; Casadio, C.; Ardissone, F.; Ruffini, E.; et al. XRCC1 and ERCC1 variants modify malignant mesothelioma risk: A case-control study. Mutat. Res. 2011, 708, 11–20. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Betti, M.; Casalone, E.; Ferrante, D.; Aspesi, A.; Morleo, G.; Biasi, A.; Sculco, M.; Mancuso, G.; Guarrera, S.; Righi, L.; et al. Germline mutations in DNA repair genes predispose asbestos-exposed patients to malignant pleural mesothelioma. Cancer Lett. 2017, 405, 38–45. [Google Scholar] [CrossRef] [PubMed]
  47. Ferrante, D.; Mirabelli, D.; Tunesi, S.; Terracini, B.; Magnani, C. Pleural mesothelioma and occupational and non-occupational asbestos exposure: A case-control study with quantitative risk assessment. Occup. Environ. Med. 2016, 73, 147–153. [Google Scholar] [CrossRef] [Green Version]
  48. Du, P.; Zhang, X.; Huang, C.C.; Jafari, N.; Kibbe, W.A.; Hou, L.; Lin, S.M. Comparison of Beta-value and M-value methods for quantifying methylation levels by microarray analysis. BMC Bioinform. 2010, 11, 587. [Google Scholar] [CrossRef] [Green Version]
  49. Lehne, B.; Drong, A.W.; Loh, M.; Zhang, W.; Scott, W.R.; Tan, S.T.; Afzal, U.; Scott, J.; Jarvelin, M.R.; Elliott, P.; et al. A coherent approach for analysis of the Illumina HumanMethylation450 BeadChip improves data quality and performance in epigenome-wide association studies. Genome Biol. 2015, 16, 37. [Google Scholar] [CrossRef] [Green Version]
  50. Campanella, G.; Polidoro, S.; Di Gaetano, C.; Fiorito, F.; Guarrera, S.; Krogh, V.; Palli, D.; Panico, S.; Sacerdote, C.; Tumino, R.; et al. Epigenetic signatures of internal migration in Italy. Int. J. Epidemiol. 2015, 44, 1442–1449. [Google Scholar] [CrossRef] [Green Version]
  51. Di Gaetano, C.; Voglino, F.; Guarrera, S.; Fiorito, G.; Rosa, F.; Di Blasio, A.M.; Manzini, P.; Dianzani, I.; Betti, M.; Cusi, D.; et al. An overview of the genetic structure within the Italian population from genome-wide data. PLoS ONE 2012, 7, e43759. [Google Scholar] [CrossRef] [PubMed]
  52. Houseman, E.A.; Accomando, W.P.; Koestler, D.C.; Christensen, B.C.; Marsit, C.J.; Nelson, H.H.; Wiencke, J.K.; Kelsey, K.T. DNA methylation arrays as surrogate measures of cell mixture distribution. BMC Bioinform. 2012, 13, 86. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Ehrich, M.; Nelson, M.R.; Stanssens, P.; Zabeau, M.; Liloglou, T.; Xinarianos, G.; Cantor, C.R.; Field, J.K.; van den Boom, D. Quantitative high-throughput analysis of DNA methylation patterns by base-specific cleavage and mass spectrometry. Proc. Natl. Acad. Sci. USA 2005, 102, 15785–15790. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Manhattan plot for EWAS test on 450 k single CpGs. Single-CpG DNAm was used as dependent variable adjusting for age, gender, White blood cells (WBCs: monocytes, granulocytes, natural killer, B cells, CD4+ T and CD8+ T) estimation, population stratification and technical variability. FDR post hoc line highlights statistically significant differences between cases and controls at single CpG level.
Figure 1. Manhattan plot for EWAS test on 450 k single CpGs. Single-CpG DNAm was used as dependent variable adjusting for age, gender, White blood cells (WBCs: monocytes, granulocytes, natural killer, B cells, CD4+ T and CD8+ T) estimation, population stratification and technical variability. FDR post hoc line highlights statistically significant differences between cases and controls at single CpG level.
Cancers 13 02636 g001
Figure 2. Expression profiles in normal human haematopoiesis. MLLT1 (A) and FKBP5 (B) expression profiles in normal human haematopoiesis as reported in the Blood Spot database (http://servers.binf.ku.dk/bloodspot/, accessed on 26 May 2021).
Figure 2. Expression profiles in normal human haematopoiesis. MLLT1 (A) and FKBP5 (B) expression profiles in normal human haematopoiesis as reported in the Blood Spot database (http://servers.binf.ku.dk/bloodspot/, accessed on 26 May 2021).
Cancers 13 02636 g002
Table 1. Differential DNAm analyses ordered by effect size. Information about single-CpGs, including location-related values and model outputs (effect size, standard error, p values).
Table 1. Differential DNAm analyses ordered by effect size. Information about single-CpGs, including location-related values and model outputs (effect size, standard error, p values).
Probe IDChrMap PositionGene SymbolUcsc Refgene GroupSnp ProbeEffect SizeStandard Errorp ValueFdrSignificance
cg0286923512124726864 rs732235270.0580.0111.3 × 10−70.028
cg03546163635654363FKBP55′UTR −0.0890.0161.3 × 10−70.028*§₼
cg0235304812124718401 0.0330.0062.2 × 10−70.032
cg06633438196272158MLLT1Body 0.0690.0141.0 × 10−60.049*§₼
cg188603291343354421C13orf30TSS1500 0.0500.0101.3 × 10−60.049
cg1978219014103487004CDC42BPBBody 0.0430.0091.2 × 10−60.049
cg06834916595610 0.0370.0081.4 × 10−60.049
cg094796501685578516 rs48434490.0370.0071.2 × 10−60.049
cg266809891685560739 rs803326600.0360.0077.6 × 10−70.049
cg254095541234871422 0.0340.0071.1 × 10−60.049
cg012013991630793389ZNF629Body 0.0300.0066.1 × 10−70.049
cg1728326611111717611ALG9Body −0.0300.0061.1 × 10−60.049
Control group was set as reference. Adjustment covariates: age, gender, population stratification, WBCs (monocytes, granulocytes, natural killer, B cells, CD4+ T and CD8+ T) estimation and technical variability. *: statistically significant at p value < 0.05; §: statistically significant at FDR post hoc adjustments. ₼: statistically significant at beta = 0.01.
Table 2. Disease discrimination test considering (AUC) comparison between baseline model and models additionally including single-CpG.
Table 2. Disease discrimination test considering (AUC) comparison between baseline model and models additionally including single-CpG.
ModelAUCDeLong’s Test
BM (asbestos exposure, age, gender and WBCs)0.75Reference
BM + cg03546163 (FKBP5)0.892.1 × 10−7
BM + cg06633438 (MLLT1)0.896.3 × 10−8
Models are shown as baseline model (BM) or BM + Single CpG DNAm. AUC Differences between considered model and BM were estimated with the DeLong’s test.
Table 3. Interaction between asbestos exposure and single CpG DNAm on the MPM Odds ratios.
Table 3. Interaction between asbestos exposure and single CpG DNAm on the MPM Odds ratios.
DNAmAsbestos ExposureORStd. Error95% CIp Value
cg03546163 (FKBP5)
HypoLow2.791.511.26|6.330.013
HyperHigh7.211.543.17|17.274.6 × 10−6
HypoHigh20.841.598.71|53.965.5 × 10−11
cg06633438 (MLLT1)
HyperLow1.291.630.70|3.810.258
HypoHigh7.271.553.17|17.655.3 × 10−6
HyperHigh11.711.574.97|29.645.9 × 10−8
Reference for cg03546163 in FKBP5: hypermethylation and low asbestos exposure levels; Reference for cg06633438 in MLLT1: hypomethylation and low asbestos exposure levels.
Table 4. Descriptive characteristics of cancer-free control group.
Table 4. Descriptive characteristics of cancer-free control group.
VariableControls (Male 100, Female 37)
Min1st QMedianMean3rd QMax
Age41.6057.4165.6564.5972.6390.94
Asbestos exposure−2.71−0.97−0.48−0.440.091.73
Monocytes0.000.050.060.070.080.26
Granulocytes0.360.540.600.620.680.99
Natural Killer0.000.040.070.080.110.29
B cells0.000.070.090.090.110.19
CD4+ T0.000.100.140.140.190.35
CD8+ T0.000.030.060.070.100.23
Minimum (Min), First Quartile (1st Q), Median, Mean, Third Quartile (3rt Q) and Maximum (Max) of variables related to cancer-free controls.
Table 5. Descriptive characteristics of MPM group.
Table 5. Descriptive characteristics of MPM group.
VariableCases (Male 113, Female 50)
Min1st QMedianMean3rd QMax
Age33.9061.1968.6867.5975.1790.80
Asbestos exposure−2.71−0.210.390.370.982.94
Monocytes0.000.050.070.080.100.20
Granulocytes0.370.670.740.740.811.03
Natural Killer0.000.020.050.060.080.23
B cells0.000.050.060.060.080.16
CD4+ T0.000.030.070.080.110.22
CD8+ T0.000.000.020.030.040.22
Minimum (Min), First Quartile (1st Q), Median, Mean, Third Quartile (3rt Q) and Maximum (Max) of variables related to MPM cases.
Table 6. Descriptive characteristics of smoking status stratified by disease.
Table 6. Descriptive characteristics of smoking status stratified by disease.
Smoking HabitsCases (163)Controls (137)
n%n%
Current smokers2917.793021.90
Former smokers5433.136043.80
Never smokers7546.014734.31
n and % of the three levels of smoking status stratified by disease.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cugliari, G.; Allione, A.; Russo, A.; Catalano, C.; Casalone, E.; Guarrera, S.; Grosso, F.; Ferrante, D.; Sculco, M.; La Vecchia, M.; et al. New DNA Methylation Signals for Malignant Pleural Mesothelioma Risk Assessment. Cancers 2021, 13, 2636. https://doi.org/10.3390/cancers13112636

AMA Style

Cugliari G, Allione A, Russo A, Catalano C, Casalone E, Guarrera S, Grosso F, Ferrante D, Sculco M, La Vecchia M, et al. New DNA Methylation Signals for Malignant Pleural Mesothelioma Risk Assessment. Cancers. 2021; 13(11):2636. https://doi.org/10.3390/cancers13112636

Chicago/Turabian Style

Cugliari, Giovanni, Alessandra Allione, Alessia Russo, Chiara Catalano, Elisabetta Casalone, Simonetta Guarrera, Federica Grosso, Daniela Ferrante, Marika Sculco, Marta La Vecchia, and et al. 2021. "New DNA Methylation Signals for Malignant Pleural Mesothelioma Risk Assessment" Cancers 13, no. 11: 2636. https://doi.org/10.3390/cancers13112636

APA Style

Cugliari, G., Allione, A., Russo, A., Catalano, C., Casalone, E., Guarrera, S., Grosso, F., Ferrante, D., Sculco, M., La Vecchia, M., Pirazzini, C., Libener, R., Mirabelli, D., Magnani, C., Dianzani, I., & Matullo, G. (2021). New DNA Methylation Signals for Malignant Pleural Mesothelioma Risk Assessment. Cancers, 13(11), 2636. https://doi.org/10.3390/cancers13112636

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