Next Article in Journal
Experimentally Determined Long Intrinsically Disordered Protein Regions Are Now Abundant in the Protein Data Bank
Next Article in Special Issue
Epigenetics in Lewy Body Diseases: Impact on Gene Expression, Utility as a Biomarker, and Possibilities for Therapy
Previous Article in Journal
Intratumoral Administration of a Novel Cytotoxic Formulation with Strong Tissue Dispersive Properties Regresses Tumor Growth and Elicits Systemic Adaptive Immunity in In Vivo Models
Previous Article in Special Issue
Pathogenomics of Uterine Fibroids Development
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 21
Issue 12
10.3390/ijms21124494
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Colorectal Cancer Early Detection in Stool Samples Tracing CpG Islands Methylation Alterations Affecting Gene Expression

by
Ana Florencia Vega-Benedetti
1,†,
Eleonora Loi
1,†,
Loredana Moi
1,
Sandra Orrù
2,
Pina Ziranu
3,
Andrea Pretta
3,
Eleonora Lai
3,
Marco Puzzoni
3,
Letizia Ciccone
4,
Andrea Casadei-Gardini
5,
Francesco Cabras
6,
Federica Fortunato
6,
Angelo Restivo
6,
Luigi Zorcolo
6,
Mario Scartozzi
3 and
Patrizia Zavattari
1,*
1
Department of Biomedical Sciences, Unit of Biology and Genetics, University of Cagliari, 09042 Cagliari, Italy
2
Department of Pathology, “A. Businco” Oncologic Hospital, ASL Cagliari, 09121 Cagliari, Italy
3
Department of Medical Oncology, University Hospital of Cagliari, 09042 Cagliari, Italy
4
Bio-Rad Laboratories srl, Southern Europe, 20090 Milan, Italy
5
Division of Oncology, Department of Oncology and Hematology, University Hospital of Modena, 41125 Modena, Italy
6
Department of Surgery, Colorectal Surgery Center, University of Cagliari, 09042 Cagliari, Italy
*
Author to whom correspondence should be addressed.
†
These authors contributed equally.
Int. J. Mol. Sci. 2020, 21(12), 4494; https://doi.org/10.3390/ijms21124494
Submission received: 12 May 2020 / Revised: 14 June 2020 / Accepted: 22 June 2020 / Published: 24 June 2020
(This article belongs to the Special Issue Epigenomics of Complex Traits and Diseases)
Browse Figures
Versions Notes

Abstract

:
Colorectal cancer (CRC) is a major cause of cancer mortality. Early diagnosis is relevant for its prevention and treatment. Since DNA methylation alterations are early events in tumourigenesis and can be detected in cell-free DNA, they represent promising biomarkers for early CRC diagnosis through non-invasive methods. In our previous work, we identified 74 early altered CpG islands (CGIs) associated with genes involved in cell cross-talking and cell signalling pathways. The aim of this work was to test whether methylation-based biomarkers could be detected in non-invasive matrices. Our results confirmed methylation alterations of GRIA4 and VIPR2 in CRC tissues, using MethyLight, as well as in stool samples, using a much more sensitive technique as droplet digital PCR. Furthermore, we analysed expression levels of selected genes whose promoter CGIs were hypermethylated in CRC, detecting downregulation at mRNA and protein levels in CRC tissue for GRIA4, VIPR2, SPOCK1 and SLC6A3. Most of these genes were already lowly expressed in colon normal tissues supporting the idea that cancer DNA methylation targets genes already barely expressed in the matched normal tissues. Our study suggests GRIA4 and VIPR2 as biomarkers for early CRC diagnosis using stool samples and confirms downregulation of genes hypermethylated in CRC.

1. Introduction

Colorectal cancer (CRC) is the third most common cancer worldwide and one of the most frequent cause of cancer death. Early diagnosis has a relevant importance for its prevention and treatment [1]. Therefore, to improve diagnosis in the early stage of colorectal cancer, researchers focused on the identification of biomarkers at genetic and epigenetic levels [2,3,4,5,6]. CRC diagnosis requires examination of tissue collected during colonoscopy. However, colonoscopy is a costly and invasive procedure, thus screening tests on stool samples have been developed. These tests include faecal occult blood test (FOBT) and faecal immunochemical test (FIT) but, although they can offer indications for possible CRC, their specificity and sensitivity are still limited [7].
Methylation alterations occur early in tumour progression [4,6,8] predating genetic mutations, and can be detected in non-invasive matrices, such as stool and plasma samples [9]. Therefore, they represent extremely useful biomarkers for early CRC diagnosis. In fact, several commercial kits analysing methylation-based biomarkers are available. However, these kits have shown low specificity and sensitivity for the detection of early stages CRC and advanced adenomas, generating false positive and negative results [10,11,12]. Recently, our research group identified 74 altered CpG Islands (CGI) in both CRCs and adenomas confirming that methylation alterations are early events during CRC tumorigenesis [4]. As highlighted in our previous work, CGIs that undergo early methylation alteration are associated with genes involved in cell cross-talking and cell signalling pathways, such as membrane receptors, secreted signalling proteins and cell adhesion molecules. These data would suggest that the first modifications of the cell towards a carcinogenesis path involve the interaction of the cell with the other cells and with the surrounding microenvironment.
In this work, we analysed selected genes that may be not only promising biomarkers for early CRC detection but are also representative of the aforementioned most altered biological processes in colon tumorigenesis. In particular, we studied methylation status of GRIA4 and VIPR2, gene and protein expression levels of GRIA4, VIPR2, SPOCK1 and SLC6A3 in 10 CRC tissue samples and their matched normal tissues. Moreover, we explored the possible usefulness of GRIA4 and VIPR2 associated CGIs as early CRC biomarkers in stool samples.
Among the selected biomarkers, the Glutamate Ionotropic Receptor AMPA Type Subunit 4 (GRIA4) encodes a homonym subunit of the AMPA tetrameric receptor complex [13]. The principal function of this receptor type, as a cationic ion channel, is mainly performed in the central nervous system, e.g., synaptic communication [14]. Although knockdown of GRIA4 has been associated with dysregulation of genes involved in invasion and metastasis [15], its functional role in cancer is not fully elucidated. Interestingly, GRIA4-associated CGI has been found hypermethylated in CRC and adenomas [3,4,16].
The vasoactive intestinal peptide receptor 2 (VIPR2; alias VPAC2) encodes a homonym transmembrane protein that is associated with a guanine nucleotide binding protein. VIPR2 is activated through the interaction with the ligands VIP, PACAP-38 and -27 [17]. Specifically, VIPR2 receptor is involved in smooth muscle relaxation, exocrine and endocrine secretion and it has been associated with physiological processes such as the circadian activity and the immune response [17]. This type 2 receptor is expressed in the central nervous system, smooth muscle and blood vessels [17,18]. VPAC2 has been poorly studied in cancer, but it has been detected in thyroid, gastric, lung and neuroendocrine tumours, among others [19,20,21].
The SPARC/Osteonectin, Cwcv and Kazal like domains proteoglycan 1 (SPOCK1) gene encodes a glycoprotein of the extracellular matrix whose principal role is in cell–cell and cell–matrix interactions. It has been reported that this gene is frequently altered in various tumour types such as lung, prostate and colon carcinomas [22,23,24]. Previous research demonstrated that SPOCK1 achieves cell cycle regulation through the PI3K/AKT pathway, but nevertheless its function is not well elucidated [24].
The Solute Carrier Family 6 Member 3 (SLC6A3) encodes a dopamine transporter that resorbs synaptic dopamine. Its function has been mainly investigated in psychiatric disorders [25] and Parkinson’s disease [26]. Concerning its role in human malignancies, SLC6A3 is also a target gene of Hypoxia Inducible Factor 1 (HIF-1), and thus its expression has recently been studied in glioblastomas in response to oxygen decline [27]. However, some evidence describes the participation of dopamine receptors in other types of tumours [28]. In fact, there is no proof of the role of SLC6A3 in colorectal cancer. Nevertheless, it has been reported that dopamine improves the efficacy of therapy in breast and colon cancer [29]. In a different context, hypermethylation was detected in dopamine pathway genes, including SLC6A3, as a consequence of a high fat and sugar diet suggesting a link between dopaminergic signalling and metabolic control [30].
Our work confirmed the methylation alterations of GRIA4- and VIPR2-associated CGIs in CRC tissues as well as non-invasive matrices such as stool samples. Moreover, our results showed a downregulation at mRNA and protein levels in tumour colon tissue for all the analysed biomarkers.

2. Results

2.1. Biomarkers’ Selection

Twenty-four CGIs fulfilled our selection criteria (see Materials and Methods), having an area under curve (AUC) higher than 0.95 in our discovery set [4] and colon adenocarcinoma validation set from The Cancer Genome Atlas (TCGA-COAD) (Table 1).
The four biomarkers’ selection was based on a combination of the above parameters but also on the fact that they were firstly discovered by our research group [4], and for their functional roles. Two CGIs (chr11:105481126-105481422 and chr7:158936507-158938492), respectively mapping on GRIA4 and VIPR2 gene promoters, were selected for the methylation analyses. VIPR2 was selected for its functional role and the involvement of vasoactive intestinal peptide receptors in cancer [21], while GRIA4 was chosen because it has previously shown methylation alterations in CRC tissues as well in plasma and stool samples [3,4].
These two CGIs were hypermethylated in CRC samples in our discovery set and TCGA validation dataset (Figure 1).

2.2. Methylation Analyses

Further methylation analyses of GRIA4 and VIPR2 were conducted on 10 CRC paired tissue samples using MethyLight qPCR.
Methylation levels of tumour samples were compared to those of their respective normal samples.
GRIA4 showed hypermethylation in 6/10 tumour samples, while VIPR2 was hypermethylated in 7/10 tumour samples. For some tumour samples methylation levels of GRIA4 (2/10 samples) or VIPR2 (2/10 samples) were similar to their respective normal samples, while in other cases, methylation of GRIA4 (2/10 samples) and VIPR2 (1/10 sample) was not detectable, probably due to the low content of tumour cells in these tissue samples. (Table 2).
Tissue specimens were examined by a histopathologist to investigate the possible causes of these methylation differences among the analysed samples. Tumour samples with undetectable methylation or methylation levels similar to the respective normal samples showed low content of tumour cells. On the other hand, tumour samples with high methylation levels showed a high content of tumour cells.
To investigate whether these alterations can also be detected through a non-invasive procedure as biomarkers, their methylation levels were evaluated in stool samples from the same cohort by MethyLight.
GRIA4 methylation was detected in 4/10 samples, while VIPR2 methylation was detected in 7/10 samples (Table 2).
Since MethyLight might not be sensitive enough for the detection of traces of methylated DNA, the same samples were analysed using a much more sensitive technique—droplet digital PCR (ddPCR). This method allowed to detect methylation of both GRIA4 and VIPR2 in 9/10 samples (Table 2).

2.3. mRNA Expression Study

The same 10 CRC paired tissue samples were analysed for gene expression of GRIA4, VIPR2 and two additional genes, SPOCK1 and SLC6A3, whose CGIs were hypermethylated in our CRC discovery set and successfully validated in TCGA dataset (Figure 2 and Table 1) [4]. Unfortunately, we did not have enough material to conduct further methylation analyses for SPOCK1 and SLC6A3.
The mean expression of each of the four genes was markedly lower in tumour than in normal tissue (Figure 3). Of note, in our experimental study SLC6A3 revealed the greatest average expression difference between tumour and normal tissues (Figure 3).
We investigated the gene expression of the four genes between tumour and normal tissues using TCGA expression data obtained by RNA-seq. This analysis confirmed the significant lower expression level of GRIA4 and VIPR2 in tumour tissues than in normal ones. In disagreement with our experimental data, SLC6A3 expression analysis revealed higher levels in tumours and SPOCK1 presented similar expression levels in both tissues (Figure 4).

2.4. Protein Expression Study

Protein expression level of the four biomarkers were analysed by Western Blot. SPOCK1, SLC6A3 and GluR4 were significantly low-expressed in tumour respect to normal tissues, following the same pattern of the mRNA expression level (Figure 5 and Figure S1). On the contrary, VIPR2 did not show statistically significant differences of expression between tumour and normal tissues.

3. Discussion

Colorectal cancer develops through the accumulation of genetic and epigenetic aberrations. Methylation alterations in CRC could be used as biomarkers for early diagnosis, and/or prognostic and predictive markers to improve therapy. In fact, our research group identified a panel of 74 altered CGIs that discriminates tumour and adenoma tissues from normal tissues [4]. In the current work, we focused on four genes, GRIA4, VIPR2, SPOCK1 and SLC6A3, whose CGIs were hypermethylated in CRC in our previous study. Further methylation analyses were performed in 10 paired samples for GRIA4 and VIPR2, confirming their alteration in CRC. The absence of hypermethylation of GRIA4 and VIPR2 associated CGIs in CRC tissue samples with low content of tumour cells highlights the high specificity of these two biomarkers for CRC detection. To note, sample collection during surgery is a rigorous step that must be carefully controlled and may lead to possible false results. In fact, in the current study a methodological error in tissue sample collection caused apparent contradictory results. However, this contradiction has been uncovered by an accurate histopathological analysis. We analysed stool samples from the same cohort of CRC patients. Stool samples from healthy individuals were not included in this study. Nevertheless, it is expected that they have an even lower number of methylated targets in a high background of unmethylated molecules compared to stool samples of CRC patients. To note, in CRC patients the concentration of exfoliated colonocytes, released in faeces, is 4.5-fold higher than in healthy individuals [31]. Importantly, our results showed that MethyLight gave false negatives in stool samples. Nevertheless, we were able to confirm that using a very sensitive method, such as ddPCR, these two methylation alterations were detectable in all except one stool sample. The lack of methylation positivity in this sample could be due to the loss of tumour DNA from the faeces suggesting the high specificity of our biomarkers for tumour DNA. In fact, stool samples are collected from the resected tumour tissues after their washout to clean the specimens for the pathology examination. Moreover, it must be considered that the patients follow a special diet including laxative drinks to clean the bowel before surgery [32].
Therefore, the detection of these methylation-based biomarkers using ddPCR would allow an early diagnosis of CRC and the follow-up of patients after tumour surgical resection and/or treated with chemotherapy and/or radiotherapy. In fact, ddPCR can overcome the technical challenges (i.e., poor DNA quality, presence of contaminating DNA and PCR inhibitors) responsible for a reduced performance of MethyLight for the detection of low copies of methylated DNA in such samples, being 25-fold more sensitive than conventional MethyLight [33]. Our results confirm that ddPCR is the gold standard method for detecting low copies of tumour DNA in non-invasive matrices.
Given that FOBT and other biomarkers screening tests recently introduced in the sanitary system give a high number of false positive results [7,9], our methylation biomarkers can potentially improve the current diagnosis system, reinforcing our previous results for GRIA4-associated CGI [3,4].
Moreover, we studied the transcript and protein expression levels of the four selected biomarkers, whose associated CGIs are located in promoter regions, demonstrating their downregulation in tumour tissues. Not all samples reflect a direct correlation between RNA and protein expression given that factors, such as post-transcriptional modification, translational regulation and protein half-lives, regulate gene expression. However, the mean values showed a correlation between transcript and protein expression levels. These results confirmed the association between promoter hypermethylation and gene downregulation already reported, considering methylation as a mechanism for gene transcriptional inactivation in cancer [34,35]. Moreover, the expression levels of most of the genes studied in this work, were already low in colon normal tissues (GTEx data, https://www.gtexportal.org) supporting the idea that methylation in cancer targets genes barely expressed in tissues where tumour arises [4,36,37,38]. Interestingly, our study shows that these genes are further repressed in tumour tissues that can be observed only by means of a targeted gene expression analysis as previously reported [4,36,39]. This further downregulation can be detected only by using methods such as qRT-PCR. In fact, the background levels of hybridization to probes in gene expression microarray or the low sequencing depth for low expressed genes in RNA-seq do not allow us to detect small gene expression changes of low-expressed genes [4,40,41].
These might be the possible reasons for the lack of validation of our gene expression results for SPOCK1 and SLC6A3 in TCGA RNA-seq data, unlike GRIA4 and VIPR2.
The current study confirmed GRIA4 methylation alteration in CRC [3,4] and gives further support of its possible role in tumour by means of the altered expression detected at the RNA and protein level.
There is little evidence of the role of the VPAC receptors in cancer although their expression was already reported in a variety of tumours. In particular, VPAC1 receptors are expressed in malignant epithelial neoplasms, while expression of VPAC2 receptors has mainly been found in some leiomyomas and gastrointestinal stromal tumours [18,42,43,44]. Our results demonstrated protein and transcript expression of VIPR2 in normal and tumour tissue enriching the knowledge of its possible involvement in CRC. However, we did not detect a statistically significant reduction of VIPR2 at protein level in tumour tissues, probably for the large variability among the samples, although a tendency towards downregulation can be observed.
Our finding that SPOCK1 is downregulated in CRC is in contrast with TCGA data and with previous research that demonstrated its overexpression in colon cancer [24,45] and in other tumour types [21,23]. In TCGA data, this gene was similarly expressed in tumour and normal samples. Moreover, Sanz-Pamplona et al. (2014) showed that SPOCK1 was upregulated in the normal mucosa adjacent to CRC tissue (minimum distance of 10 cm) compared to colon mucosa from healthy donors. However, its expression in the normal mucosa samples from CRC patients was similar to that of the tumour. In fact, at the molecular level peritumoural cells resembled more malignant cells than normal cells, i.e., the surrounding environment could already present alterations similar to the tumour tissue. Protein–protein network analysis showed that SPOCK1 is a protein secreted by adjacent mucosa and interacts with a receptor in tumour [45]. Of note, some of our tumour samples have shown higher SPOCK1 expression than their matched-normal sample, indicating that its expression is highly variable. Therefore, further studies to clarify these discordant results are needed.
The expression of the dopamine transporter SLC6A3 was significantly downregulated in our CRC samples, both at mRNA and protein levels. In contrast, this gene resulted upregulated in CRC in TCGA data. However, gene expression values of tumours markedly overlap those of normal samples in TCGA dataset. The discrepancy with our results might be related to the low expression of this gene in the colon (GTEx data: about 0.8 TPM) and even reduced expression in tumours that might be not detected by techniques such as RNA-seq and this result should be validated in a larger sample cohort by ultrasensitive methods such as ddPCR. SLC6A3 possible role in cancer might be due to the well-known protective effects of dopamine in tumours [28]. Previously, it was reported that Slc6a3 was hypermethylated in the promoter region in response to high-fat-sucrose diet and downregulated in the nervous system of prenatally stressed female rats and murine models of diet-induced obesity [46,47]. In addition, methylation aberrations at SLC6A3 gene were related to triglyceride level and obesity in humans [30,48]. Hence, the hypermethylation event at the promoter region of SLC6A3, detected previously by our group [4], and its low expression in CRC, reported in the current study, could suggest an alteration in metabolism, such as lipogenic pathways, and probably a oxidative stress environment leading to malignancy transformation [49]. This statement can also be supported by the involvement of SLC6A3 in response to hypoxia, in glioblastomas, as one of the target genes of HIF-1 [27].
In summary, our work was focused on the methylation and expression analyses of GRIA4, VIPR2, SPOCK1 and SLC6A3, belonging to gene families involved in the crosstalk between tumour cells and the environment. From our point of view, the identified hypermethylation and reduced expression of these normally low expressed genes but further downregulated in cancer, underline their functional involvement in a defined program of gene silencing during tumour transformation. This is in agreement with the “epigenetic switching” concept allowing the maintenance of a stable and permanent repression of genes important for cell proliferation and tumorigenesis leading to the restoring of a stem-cell-like state [50,51]. Of note, epigenome editing approaches may allow to re-establish the normal methylation and gene expression patterns of these genes, possibly reverting the tumour phenotype to that characteristic of a normal cell.

4. Materials and Methods

4.1. Tissue Samples

Tumour and matched-normal fresh-frozen (FF) tissue samples of 10 CRC patients were collected from the Department of General Surgery of the University of Cagliari (Italy). Normal samples were taken at a distance >10 cm from the neoplastic tissue. All the biological samples analysed were obtained with written informed consent signed from patients and ethical approval granted by the relative Ethics Committee.
Histopathological analysis was performed for all the tissue samples. Frozen section slides underwent standard haematoxylin and eosin (H&E) staining. Microscope images were acquired of each individual slide.
Patients’ clinical data are reported in Table 3.

4.2. Stool Samples

Stool samples were collected from the same cohort of ten CRC patients taken intraoperatively from the bowel resection specimen. All stools samples were immediately frozen after collection and stored at −80 °C until being processed.

4.3. Marker Selection

Biomarkers were selected from a panel of 74 altered CGIs, previously identified from a genome-wide methylation study of 18 primary CRCs and four matched peritumoural samples, 21 colorectal adenomas and three matched-normal intestinal mucosa samples by our research group [4].
Methylation data from this cohort (discovery set) and methylation data of TCGA-COAD were analysed for the biomarkers selection using the following pipeline summarized in Figure 6. In the first step, CGIs for which mean beta value was lower than 0.25 in peritumoural and normal samples were retained. Subsequently, only CGIs with a beta value higher than 0.45 in at least 75% of tumours and a beta value higher than 0.25 in not more than 25% of peritumoural and normal samples were selected. Receiver operating characteristic (ROC) curves for the 36 selected biomarkers were performed by R “ROCR” package. A group of 24 CGIs showed an AUC higher than 0.95 in our cohort and in TCGA-COAD.
Four biomarkers were finally selected for a combination of parameters (AUC, Δβ), functional roles and because they were originally identified by our group [4]. Additional filtering criteria included feasibility of the assay. Further methylation analyses were conducted only for two out of the four selected biomarkers (GRIA4 and VIPR2), because of the lack of biological material for the other assays.

4.4. Methylation Analyses

Genomic DNA was extracted from tumour and peritumoural tissue using the DNeasy Blood & Tissue Kit (Qiagen, Hilden, Germany). Stool DNA extraction was performed using the QIAamp DNA Stool Mini Kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions.
DNA samples were bisulfite converted using EZ DNA Methylation Gold Kit TM (Zymo Research, Irvine, CA, USA) according to the manufacturer’s instructions. Quality control and quantification of DNA were performed before and after bisulfite conversion.
GRIA4 and VIPR2 methylation was assessed by using MethyLight qPCR [52] in tissue and stool samples. The reference repetitive element Alu was used in a methylation-independent control reaction to normalize the amount of DNA input. Primers and probes were designed using Beacon Designer™ (Premier Biosoft, San Francisco, CA, USA) (Table 4). The probes were labelled with the 6-Carboxyfluorescein (6-FAM) fluorophore at the 5′ end. A primer-probe mix containing 300 nM of each primer and 100 nM of the probe was prepared. Each assay was performed in triplicate using: 15 μL of TaqMan Genotyping Master mix (Applied Biosystems, Foster City, CA, USA), 4.5 μL of primer-probe mix, 5 μL of bisulfite-converted DNA (10 ng/μL) and 5.5 μL of RNase-free water. A fully-methylated DNA (Human Methylated & Non-Methylated (WGA) DNA Set; Zymo Research, Irvine, CA, USA) was used as a positive control for the reaction. The experiment was conducted on a DNA Engine Opticon 2 Real-Time Cycler (Bio-Rad, Hercules, CA, USA) using the following thermal conditions: initial PCR activation step at 95 °C for 10 min, followed by 50 cycles of denaturation step at 95 °C for 15 s and annealing/extension step at 60 °C for 60 s.
A comparison between CRC and their matched-normal samples have been performed. CRC samples were considered as hypermethylated if delta Ct (mean Ct normal − mean Ct tumour) > 2 and as hypomethylated if delta Ct (mean Ct normal − mean Ct tumour) < −2.
Methylation of the two biomarkers was also evaluated by ddPCR in the stool samples from CRC patients using the same primers and probes used in the MethyLight assay (Table 4).
ddPCR reactions containing 2 × ddPCR Supermix for probes (Bio-Rad, Hercules, CA, USA), forward and reverse primer (900 nM), probe (250 nM) and 2 to 5 μL bisulfite-converted DNA in a final volume of 20 μL were partitioned into ~20,000 oil-emulsified droplets per well and replicated in three wells using a Bio-Rad QX200 droplet generator (Bio-Rad, Hercules, CA, USA). The droplets were transferred into 96-well plates and PCR was performed using the following conditions: 10 min at 95 °C, 40 cycles of 30 s at 95 °C followed by 60 s at 60 °C, then 5 min at 4 °C, 5 min at 95 °C. Plates were read on a Bio-Rad QX200 droplet reader (Bio-Rad, Hercules, CA, USA).
Data were analysed using the QuantaSoft 1.7.4 software (Bio-Rad). The droplet counts (positive or negative) from all replicated wells were combined to yield a ‘merged’ well. Concentration values (number of copies/μL) and Poisson confidence intervals were computed for each “merged” well.

4.5. mRNA Expression Analysis

Total RNA was extracted from tumoural and normal samples using the RNeasy Mini Kit (Qiagen, Hilden, Germany) following manufacturer’s protocol. An aliquot of 1μg RNA/sample was retrotranscribed using the High Capacity Kit (Applied Biosystems, Carlsbad, CA, USA). Gene expression was evaluated by qPCR using SsoAdvanced™ Universal SYBR® Green Supermix (Bio-Rad, Hercules, CA, USA) for each gene tested, GRIA4, SLC6A3, SPOCK1 and VIPR2, and for the endogenous gene TFRC. Primers sequences can be found in Table 5. PCR conditions were primary denaturation at 95 °C for 2 min followed by 50 cycles of denaturation at 95 °C for 15 s and annealing/extension at 60 °C for 1 min. After the amplification cycles, melting curves were produced by increasing the temperature from 65 °C to 95 °C holding each temperature for 5 sec and reading fluorescence every 0.5 °C.
The transcript levels were quantified using the ΔΔCT method. Statistical analyses were done by comparing the average ΔCt of the sample groups of interest using Welch’s t-test. We considered upregulation or downregulation when there existed statistical significant differences in the expression level of tumour tissues compared to normal tissues, whereas the expression between tissues was considered similar when the statistical analysis was not significant.

4.6. Protein Expression Analysis

The protein expression level of SLC6A3, SPOCK1, VIPR2 and GluR4 was studied by Western blot. Proteins from colon tissue were isolated using the Membrane Protein Extraction Kit (Mem-PERTM Plus, Thermo Fisher Scientific, Waltham, MA, USA) following the manufacturer’s instructions. Protein extracts were fractioned by 10% SDS-polyacrylamide gel and separated at 100 V during 90 min. Then, the proteins were transferred onto 0.45 µm nitrocellulose membrane, blocked in 5% milk overnight at 4 °C and then incubated 2 h with primary antibody GluR4 (0.5 µg/mL, PA5-18931, Thermo Fisher Scientific, Waltham, MA, USA), SPOCK1 (1 µg/mL, MA5-24039, Thermo Fisher Scientific, Waltham, MA, USA), VIPR2 (0.5 µg/mL, AB2266, Millipore, Burlington, MA, USA) and SLC6A3 (1 µg/mL, SAB2502027, Sigma-Aldrich, St. Louis, MO, USA). After hybridization with primary antibody, the membrane was washed with Tris-buffered saline containing Tween-20 three times, then incubated with horseradish peroxidase-labelled secondary antibody (Jackson ImmunoResearch, Ely, UK) for 1 h at room temperature and washed with Tris-buffered saline containing Tween-20 three times. Final detection was performed with ECL Chemiluminescent Western blotting reagents (Bio-Rad, Hercules, CA, USA). An antibody against NaK ATPase (1:40000, Ab76020, Abcam, Cambridge, UK) was used for gel-loading control. Unfortunately, there was not enough protein lysate of the subjects 19 and 21 to study the expression level of SPOCK1.
Western blot signals were quantified using ImageJ program. The intensity of each band was normalized with respect to that of NaK ATPase. The assay was performed in duplicate. Statistical differences between tumour and normal samples were calculated using Welch’s t-test.

4.7. Validation Analyses

Methylation and RNA-seq data from TCGA, including colorectal samples, were used to validate our methylation and gene expression results. Statistical differences between tumour and normal samples were calculated using Welch’s t-test.

5. Conclusions

In conclusion, our study suggests GRIA4 and VIPR2 as biomarkers for CRC early detection by non-invasive methods, in particular in stool samples. Further studies are needed to explore methylation detection of these two biomarkers in stool samples from population subjected to CRC screening and prove their specificity and sensitivity for early diagnosis. Moreover, we confirm the downregulation of genes hypermethylated in tumours and suggest the role of these genes in CRC although functional studies are needed to support this hypothesis.

Supplementary Materials

The following are available online at https://www.mdpi.com/1422-0067/21/12/4494/s1.

Author Contributions

Conceptualization, P.Z. (Patrizia Zavattari); Data curation, A.F.V.-B., E.L. (Eleonora Loi), L.M., P.Z. (Pina Ziranu), A.P., E.L. (Eleonora Lai), M.P., A.C.-G., M.S. and P.Z. (Patrizia Zavattari); Formal analysis, A.F.V.-B. and E.L. (Eleonora Loi); Funding acquisition, L.Z. and P.Z. (Patrizia Zavattari); Investigation, A.F.V.-B., E.L. (Eleonora Loi), L.M., S.O., P.Z. (Pina Ziranu), A.P., E.L. (Eleonora Lai), M.P., A.C.-G., M.S. and P.Z. (Patrizia Zavattari); Methodology, A.F.V.-B., E.L. (Eleonora Loi), L.M., S.O., L.C. and P.Z. (Patrizia Zavattari); Project administration, L.Z. and P.Z. (Patrizia Zavattari); Resources, L.C., F.C., F.F., A.R., L.Z. and P.Z. (Patrizia Zavattari); Software, E.L. (Eleonora Loi); Supervision, M.S. and P.Z. (Patrizia Zavattari); Visualization, A.F.V.-B., E.L. (Eleonora Loi) and S.O.; Writing—original draft, A.F.V.-B., E.L. (Eleonora Loi) and P.Z. (Patrizia Zavattari); Writing—review & editing, A.F.V.-B., E.L. (Eleonora Loi), M.S. and P.Z. (Patrizia Zavattari) All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by grants from Fondazione Banco di Sardegna (2012), Fondo per la Ricerca Locale (ex 60%), Università di Cagliari, POR FESR Sardegna (1C-177) to PZ, Regione Autonoma della Sardegna (CRP-79303) to PZ and LZ, Italian University Minister PON-AIM (F74I18000200001) to PZ.

Acknowledgments

We gratefully acknowledge S. Nanni from the “Università Cattolica del Sacro Cuore”, Italy, for the use of the ddPCR system (BioRad) for methylation assays in stool samples. We also thank A. Fadda for helping to draw some methylation assays.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Siegel, R.L.; Miller, K.D.; Jemal, A. Cancer statistics, 2016. CA Cancer J. Clin. 2016, 66, 7–30. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Bettegowda, C.; Sausen, M.; Leary, R.J.; Kinde, I.; Wang, Y.; Agrawal, N.; Bartlett, B.R.; Wang, H.; Luber, B.; Alani, R.M.; et al. Detection of circulating tumor DNA in early- and late-stage human malignancies. Sci. Transl. Med. 2014, 6, 224ra224. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Barault, L.; Amatu, A.; Siravegna, G.; Ponzetti, A.; Moran, S.; Cassingena, A.; Mussolin, B.; Falcomata, C.; Binder, A.M.; Cristiano, C.; et al. Discovery of methylated circulating DNA biomarkers for comprehensive noninvasive monitoring of treatment response in metastatic colorectal cancer. Gut 2018, 67, 1995–2005. [Google Scholar] [CrossRef] [PubMed]
  4. Fadda, A.; Gentilini, D.; Moi, L.; Barault, L.; Leoni, V.P.; Sulas, P.; Zorcolo, L.; Restivo, A.; Cabras, F.; Fortunato, F.; et al. Colorectal cancer early methylation alterations affect the crosstalk between cell and surrounding environment, tracing a biomarker signature specific for this tumor. Int. J. Cancer 2018, 143, 907–920. [Google Scholar] [CrossRef] [Green Version]
  5. Wright, M.; Beaty, J.S.; Ternent, C.A. Molecular Markers for Colorectal Cancer. Surg. Clin. North Am. 2017, 97, 683–701. [Google Scholar] [CrossRef]
  6. Vega-Benedetti, A.F.; Loi, E.; Moi, L.; Blois, S.; Fadda, A.; Antonelli, M.; Arcella, A.; Badiali, M.; Giangaspero, F.; Morra, I.; et al. Clustered protocadherins methylation alterations in cancer. Clin. Epigenetics 2019, 11, 100. [Google Scholar] [CrossRef] [Green Version]
  7. Morikawa, T.; Kato, J.; Yamaji, Y.; Wada, R.; Mitsushima, T.; Shiratori, Y. A comparison of the immunochemical fecal occult blood test and total colonoscopy in the asymptomatic population. Gastroenterology 2005, 129, 422–428. [Google Scholar] [CrossRef]
  8. Luo, Y.; Wong, C.J.; Kaz, A.M.; Dzieciatkowski, S.; Carter, K.T.; Morris, S.M.; Wang, J.; Willis, J.E.; Makar, K.W.; Ulrich, C.M.; et al. Differences in DNA methylation signatures reveal multiple pathways of progression from adenoma to colorectal cancer. Gastroenterology 2014, 147, 418–429.e8. [Google Scholar] [CrossRef] [Green Version]
  9. Toiyama, Y.; Okugawa, Y.; Goel, A. DNA methylation and microRNA biomarkers for noninvasive detection of gastric and colorectal cancer. Biochem. Biophys. Res. Commun. 2014, 455, 43–57. [Google Scholar] [CrossRef] [Green Version]
  10. Vedeld, H.M.; Goel, A.; Lind, G.E. Epigenetic biomarkers in gastrointestinal cancers: The current state and clinical perspectives. Semin. Cancer Biol. 2018, 51, 36–49. [Google Scholar] [CrossRef]
  11. Imperiale, T.F.; Ransohoff, D.F.; Itzkowitz, S.H.; Levin, T.R.; Lavin, P.; Lidgard, G.P.; Ahlquist, D.A.; Berger, B.M. Multitarget stool DNA testing for colorectal-cancer screening. New Engl. J. Med. 2014, 370, 1287–1297. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Epigenomics AG. EIFU Epi ProColon 2.0 CE- Rev 7, Instructions for Use; Epigenomics AG: Berlin, Germany, 2016. [Google Scholar]
  13. Traynelis, S.F.; Wollmuth, L.P.; McBain, C.J.; Menniti, F.S.; Vance, K.M.; Ogden, K.K.; Hansen, K.B.; Yuan, H.; Myers, S.J.; Dingledine, R. Glutamate receptor ion channels: Structure, regulation, and function. Pharmacol. Rev. 2010, 62, 405–496. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Martin, S.; Chamberlin, A.; Shinde, D.N.; Hempel, M.; Strom, T.M.; Schreiber, A.; Johannsen, J.; Ousager, L.B.; Larsen, M.J.; Hansen, L.K.; et al. De Novo Variants in GRIA4 Lead to Intellectual Disability with or without Seizures and Gait Abnormalities. Am. J. Hum. Genet. 2017, 101, 1013–1020. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Luksch, H.; Uckermann, O.; Stepulak, A.; Hendruschk, S.; Marzahn, J.; Bastian, S.; Staufner, C.; Temme, A.; Ikonomidou, C. Silencing of selected glutamate receptor subunits modulates cancer growth. Anticancer Res. 2011, 31, 3181–3192. [Google Scholar] [PubMed]
  16. Hauptman, N.; Jevsinek Skok, D.; Spasovska, E.; Bostjancic, E.; Glavac, D. Genes CEP55, FOXD3, FOXF2, GNAO1, GRIA4, and KCNA5 as potential diagnostic biomarkers in colorectal cancer. BMC Med Genom. 2019, 12, 54. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Dickson, L.; Finlayson, K. VPAC and PAC receptors: From ligands to function. Pharmacol. Ther. 2009, 121, 294–316. [Google Scholar] [CrossRef]
  18. Reubi, J.C.; Laderach, U.; Waser, B.; Gebbers, J.O.; Robberecht, P.; Laissue, J.A. Vasoactive intestinal peptide/pituitary adenylate cyclase-activating peptide receptor subtypes in human tumors and their tissues of origin. Cancer Res. 2000, 60, 3105–3112. [Google Scholar]
  19. Szilasi, M.; Buglyo, A.; Treszl, A.; Kiss, L.; Schally, A.V.; Halmos, G. Gene expression of vasoactive intestinal peptide receptors in human lung cancer. Int. J. Oncol. 2011, 39, 1019–1024. [Google Scholar] [CrossRef] [Green Version]
  20. Schulz, S.; Mann, A.; Novakhov, B.; Piggins, H.D.; Lupp, A. VPAC2 receptor expression in human normal and neoplastic tissues: Evaluation of the novel MAB SP235. Endocr. Connect. 2015, 4, 18–26. [Google Scholar] [CrossRef] [Green Version]
  21. Moody, T.W.; Nuche-Berenguer, B.; Jensen, R.T. VIP/PACAP, and their receptors and cancer. Curr. Opin. Endocrinol. Diabetes Obes. 2016, 23, 38–47. [Google Scholar] [CrossRef] [Green Version]
  22. Miao, L.; Wang, Y.; Xia, H.; Yao, C.; Cai, H.; Song, Y. SPOCK1 is a novel transforming growth factor-beta target gene that regulates lung cancer cell epithelial-mesenchymal transition. Biochem. Biophys. Res. Commun. 2013, 440, 792–797. [Google Scholar] [CrossRef] [PubMed]
  23. Chen, Q.; Yao, Y.T.; Xu, H.; Chen, Y.B.; Gu, M.; Cai, Z.K.; Wang, Z. SPOCK1 promotes tumor growth and metastasis in human prostate cancer. Drug Des. Dev. Ther. 2016, 10, 2311–2321. [Google Scholar] [CrossRef] [Green Version]
  24. Zhang, J.; Zhi, X.; Shi, S.; Tao, R.; Chen, P.; Sun, S.; Bian, L.; Xu, Z.; Ma, L. SPOCK1 is up-regulated and promotes tumor growth via the PI3K/AKT signaling pathway in colorectal cancer. Biochem. Biophys. Res. Commun. 2017, 482, 870–876. [Google Scholar] [CrossRef] [PubMed]
  25. Salatino-Oliveira, A.; Rohde, L.A.; Hutz, M.H. The dopamine transporter role in psychiatric phenotypes. Am. J. Med Genet. Part B Neuropsychiatr. Genet. 2018, 177, 211–231. [Google Scholar] [CrossRef]
  26. Robertson, B.D.; Al Jaja, A.S.; MacDonald, A.A.; Hiebert, N.M.; Tamjeedi, R.; Seergobin, K.N.; Schwarz, U.I.; Kim, R.B.; MacDonald, P.A. SLC6A3 Polymorphism Predisposes to Dopamine Overdose in Parkinson’s Disease. Front. Neurol. 2018, 9, 693. [Google Scholar] [CrossRef] [Green Version]
  27. Bao, L.; Chen, Y.; Lai, H.T.; Wu, S.Y.; Wang, J.E.; Hatanpaa, K.J.; Raisanen, J.M.; Fontenot, M.; Lega, B.; Chiang, C.M.; et al. Methylation of hypoxia-inducible factor (HIF)-1alpha by G9a/GLP inhibits HIF-1 transcriptional activity and cell migration. Nucleic Acids Res. 2018, 46, 6576–6591. [Google Scholar] [CrossRef] [Green Version]
  28. Rubi, B.; Maechler, P. Minireview: New roles for peripheral dopamine on metabolic control and tumor growth: Let’s seek the balance. Endocrinology 2010, 151, 5570–5581. [Google Scholar] [CrossRef] [Green Version]
  29. Sarkar, C.; Chakroborty, D.; Chowdhury, U.R.; Dasgupta, P.S.; Basu, S. Dopamine increases the efficacy of anticancer drugs in breast and colon cancer preclinical models. Clin. Cancer Res. Off. J. Am. Assoc. Cancer Res. 2008, 14, 2502–2510. [Google Scholar] [CrossRef] [Green Version]
  30. Ramos-Lopez, O.; Riezu-Boj, J.I.; Milagro, F.I.; Martinez, J.A.; Project, M. Dopamine gene methylation patterns are associated with obesity markers and carbohydrate intake. Brain Behav. 2018, 8, e01017. [Google Scholar] [CrossRef] [Green Version]
  31. De Maio, G.; Rengucci, C.; Zoli, W.; Calistri, D. Circulating and stool nucleic acid analysis for colorectal cancer diagnosis. World J. Gastroenterol. 2014, 20, 957–967. [Google Scholar] [CrossRef]
  32. Lynch, M.L.; Brand, M.I. Preoperative evaluation and oncologic principles of colon cancer surgery. Clin. Colon Rectal Surg. 2005, 18, 163–173. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Yu, M.; Carter, K.T.; Makar, K.W.; Vickers, K.; Ulrich, C.M.; Schoen, R.E.; Brenner, D.; Markowitz, S.D.; Grady, W.M. MethyLight droplet digital PCR for detection and absolute quantification of infrequently methylated alleles. Epigenetics 2015, 10, 803–809. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Baylin, S.B. DNA methylation and gene silencing in cancer. Nat. Clin. Pract. Oncol. 2005, 2, S4–S11. [Google Scholar] [CrossRef] [PubMed]
  35. Jones, P.A. Functions of DNA methylation: Islands, start sites, gene bodies and beyond. Nat. Rev. Genet. 2012, 13, 484–492. [Google Scholar] [CrossRef] [PubMed]
  36. Antonelli, M.; Fadda, A.; Loi, E.; Moi, L.; Zavattari, C.; Sulas, P.; Gentilini, D.; Cameli, C.; Bacchelli, E.; Badiali, M.; et al. Integrated DNA methylation analysis identifies topographical and tumoral biomarkers in pilocytic astrocytomas. Oncotarget 2018, 9, 13807–13821. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Sproul, D.; Nestor, C.; Culley, J.; Dickson, J.H.; Dixon, J.M.; Harrison, D.J.; Meehan, R.R.; Sims, A.H.; Ramsahoye, B.H. Transcriptionally repressed genes become aberrantly methylated and distinguish tumors of different lineages in breast cancer. Proc. Natl. Acad. Sci. USA 2011, 108, 4364–4369. [Google Scholar] [CrossRef] [Green Version]
  38. Moarii, M.; Boeva, V.; Vert, J.P.; Reyal, F. Changes in correlation between promoter methylation and gene expression in cancer. BMC Genom. 2015, 16, 873. [Google Scholar] [CrossRef] [Green Version]
  39. Loi, E.; Moi, L.; Fadda, A.; Satta, G.; Zucca, M.; Sanna, S.; Amini Nia, S.; Cabras, G.; Padoan, M.; Magnani, C.; et al. Methylation alteration of SHANK1 as a predictive, diagnostic and prognostic biomarker for chronic lymphocytic leukemia. Oncotarget 2019, 10, 4987–5002. [Google Scholar] [CrossRef] [Green Version]
  40. Sproul, D.; Meehan, R.R. Genomic insights into cancer-associated aberrant CpG island hypermethylation. Brief. Funct. Genom. 2013, 12, 174–190. [Google Scholar] [CrossRef]
  41. Sims, D.; Sudbery, I.; Ilott, N.E.; Heger, A.; Ponting, C.P. Sequencing depth and coverage: Key considerations in genomic analyses. Nat. Rev. Genet. 2014, 15, 121–132. [Google Scholar] [CrossRef]
  42. Oka, H.; Jin, L.; Reubi, J.C.; Qian, X.; Scheithauer, B.W.; Fujii, K.; Kameya, T.; Lloyd, R.V. Pituitary adenylate-cyclase-activating polypeptide (PACAP) binding sites and PACAP/vasoactive intestinal polypeptide receptor expression in human pituitary adenomas. Am. J. Pathol. 1998, 153, 1787–1796. [Google Scholar] [CrossRef] [Green Version]
  43. Fernandez-Martinez, A.B.; Carmena, M.J.; Arenas, M.I.; Bajo, A.M.; Prieto, J.C.; Sanchez-Chapado, M. Overexpression of vasoactive intestinal peptide receptors and cyclooxygenase-2 in human prostate cancer. Analysis of potential prognostic relevance. Histol. Histopathol. 2012, 27, 1093–1101. [Google Scholar] [CrossRef] [PubMed]
  44. Liu, S.; Zeng, Y.; Li, Y.; Guo, W.; Liu, J.; Ouyang, N. VPAC1 overexpression is associated with poor differentiation in colon cancer. Tumour Biol. J. Int. Soc. Oncodevelopmental Biol. Med. 2014, 35, 6397–6404. [Google Scholar] [CrossRef] [PubMed]
  45. Sanz-Pamplona, R.; Berenguer, A.; Cordero, D.; Mollevi, D.G.; Crous-Bou, M.; Sole, X.; Pare-Brunet, L.; Guino, E.; Salazar, R.; Santos, C.; et al. Aberrant gene expression in mucosa adjacent to tumor reveals a molecular crosstalk in colon cancer. Mol. Cancer 2014, 13, 46. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Paternain, L.; Batlle, M.A.; De la Garza, A.L.; Milagro, F.I.; Martinez, J.A.; Campion, J. Transcriptomic and epigenetic changes in the hypothalamus are involved in an increased susceptibility to a high-fat-sucrose diet in prenatally stressed female rats. Neuroendocrinology 2012, 96, 249–260. [Google Scholar] [CrossRef]
  47. Vucetic, Z.; Carlin, J.L.; Totoki, K.; Reyes, T.M. Epigenetic dysregulation of the dopamine system in diet-induced obesity. J. Neurochem. 2012, 120, 891–898. [Google Scholar] [CrossRef] [Green Version]
  48. Arpon, A.; Milagro, F.I.; Laja, A.; Segura, V.; de Pipaon, M.S.; Riezu-Boj, J.I.; Alfredo Martinez, J. Methylation changes and pathways affected in preterm birth: A role for SLC6A3 in neurodevelopment. Epigenomics 2018, 10, 91–103. [Google Scholar] [CrossRef]
  49. Tarrado-Castellarnau, M.; de Atauri, P.; Cascante, M. Oncogenic regulation of tumor metabolic reprogramming. Oncotarget 2016, 7, 62726–62753. [Google Scholar] [CrossRef] [Green Version]
  50. Widschwendter, M.; Fiegl, H.; Egle, D.; Mueller-Holzner, E.; Spizzo, G.; Marth, C.; Weisenberger, D.J.; Campan, M.; Young, J.; Jacobs, I.; et al. Epigenetic stem cell signature in cancer. Nat. Genet. 2007, 39, 157–158. [Google Scholar] [CrossRef]
  51. Gal-Yam, E.N.; Saito, Y.; Egger, G.; Jones, P.A. Cancer epigenetics: Modifications, screening, and therapy. Annu. Rev. Med. 2008, 59, 267–280. [Google Scholar] [CrossRef]
  52. Eads, C.A.; Danenberg, K.D.; Kawakami, K.; Saltz, L.B.; Blake, C.; Shibata, D.; Danenberg, P.V.; Laird, P.W. MethyLight: A high-throughput assay to measure DNA methylation. Nucleic Acids Res. 2000, 28, E32. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Ijms 21 04494 g001 550
Figure 1. Methylation values obtained from our colorectal (CRC) discovery set and The Cancer Genome Atlas (TCGA) validation set. Genomic organization of GRIA4 (A) and VIPR2 (B), including the localization of exons and CGIs. Mean β values, resulting from the average of the samples (normal and tumour) of each probe mapping on the altered CGIs associated with GRIA4 (A) and VIPR2 (B). When more than one CGI is shown, the altered one is enclosed in a yellow box.
Figure 1. Methylation values obtained from our colorectal (CRC) discovery set and The Cancer Genome Atlas (TCGA) validation set. Genomic organization of GRIA4 (A) and VIPR2 (B), including the localization of exons and CGIs. Mean β values, resulting from the average of the samples (normal and tumour) of each probe mapping on the altered CGIs associated with GRIA4 (A) and VIPR2 (B). When more than one CGI is shown, the altered one is enclosed in a yellow box.
Ijms 21 04494 g001
Ijms 21 04494 g002 550
Figure 2. Methylation values obtained from our CRC discovery set and TCGA validation set. Genomic organization of SPOCK1 (A) and SLC6A3 (B), including the localization of exons and CGIs. Mean β values, resulting from the average of the samples (normal and tumour) of each probe mapping on the altered CGIs associated with SPOCK1 (A) and SLC6A3 (B). When more than one CGI is shown, the altered one is enclosed in a yellow box.
Figure 2. Methylation values obtained from our CRC discovery set and TCGA validation set. Genomic organization of SPOCK1 (A) and SLC6A3 (B), including the localization of exons and CGIs. Mean β values, resulting from the average of the samples (normal and tumour) of each probe mapping on the altered CGIs associated with SPOCK1 (A) and SLC6A3 (B). When more than one CGI is shown, the altered one is enclosed in a yellow box.
Ijms 21 04494 g002
Ijms 21 04494 g003 550
Figure 3. Differential gene expression analysis between tumour and normal samples. Box plot fold change values of GRIA4 (A), VIPR2 (B), SPOCK1 (C) and SLC6A3 (D) for CRC and control samples. *** indicates p value < 0.0001.
Figure 3. Differential gene expression analysis between tumour and normal samples. Box plot fold change values of GRIA4 (A), VIPR2 (B), SPOCK1 (C) and SLC6A3 (D) for CRC and control samples. *** indicates p value < 0.0001.
Ijms 21 04494 g003
Ijms 21 04494 g004 550
Figure 4. RNA expression of GRIA4 (A), VIPR2 (B), SPOCK1 (C) and SLC6A3 (D) in TCGA dataset. Tumour and normal expression of the four biomarkers. *** indicates p value < 0.0001.
Figure 4. RNA expression of GRIA4 (A), VIPR2 (B), SPOCK1 (C) and SLC6A3 (D) in TCGA dataset. Tumour and normal expression of the four biomarkers. *** indicates p value < 0.0001.
Ijms 21 04494 g004
Ijms 21 04494 g005 550
Figure 5. Protein expression level for each paired sample. (A) Representative blots of GluR4, VIPR2, SPOCK1 and SLC6A3 protein expression in 10 CRC paired tissues samples. NaK ATPase was used as loading control. (B) Box plots of GluR4, VIPR2, SPOCK1 and SLC6A3 expression. Asterisks indicate statistically significant differences (* p value < 0.05, ** p value < 0.01, *** p value < 0.001). For the SPOCK1 expression study, there was not enough protein lysate for samples 19 and 21.
Figure 5. Protein expression level for each paired sample. (A) Representative blots of GluR4, VIPR2, SPOCK1 and SLC6A3 protein expression in 10 CRC paired tissues samples. NaK ATPase was used as loading control. (B) Box plots of GluR4, VIPR2, SPOCK1 and SLC6A3 expression. Asterisks indicate statistically significant differences (* p value < 0.05, ** p value < 0.01, *** p value < 0.001). For the SPOCK1 expression study, there was not enough protein lysate for samples 19 and 21.
Ijms 21 04494 g005
Ijms 21 04494 g006 550
Figure 6. Biomarkers’ selection pipeline. Pipeline for biomarkers selection using methylation data from our previous CRC discovery set and from TCGA-COAD dataset.
Figure 6. Biomarkers’ selection pipeline. Pipeline for biomarkers selection using methylation data from our previous CRC discovery set and from TCGA-COAD dataset.
Ijms 21 04494 g006
Table
Table 1. Selected CpG islands showing an area under curve (AUC) > 0.95 in our discovery set and in TCGA-COAD validation set.
Table 1. Selected CpG islands showing an area under curve (AUC) > 0.95 in our discovery set and in TCGA-COAD validation set.
CGIGeneΔβ
Discovery Set		AUC
Discovery Set		Δβ
Validation Set		AUC
Validation Set
chr2:182321761-182323029ITGA40.371.000.350.96
chr4:156129168-156130209NPY2R0.271.000.30.97
chr4:157997166-157997686GLRB0.301.000.360.97
chr4:107956555-107957453DKK20.320.970.330.97
chr5:136834016-136835146SPOCK10.291.000.330.98
chr5:140864527-140864748PCDHGA40.351.000.380.97
chr5:1444678-1446648SLC6A30.261.000.290.98
chr5:178016558-178017670COL23A10.290.980.320.96
chr5:159399004-159399928ADRA1B0.250.970.310.97
chr6:159589636-159591319FNDC10.331.000.330.97
chr6:73330942-73333109KCNQ50.361.000.330.96
chr7:28448716-28450028CREB50.271.000.270.96
chr7:158936507-158938492VIPR20.330.960.350.96
chr8:97505747-97507607SDC20.291.000.360.96
chr8:75896528-75897116CRISPLD10.210.960.270.97
chr10:15761423-15762101ITGA80.300.960.350.97
chr11:105481126-105481422GRIA40.401.000.410.96
chr11:133938850-133939681JAM30.301.000.290.97
chr12:117798076-117799448NOS10.250.970.270.96
chr13:110958891-110960590COL4A10.330.960.370.98
chr16:23846941-23848102PRKCB0.250.980.340.97
chr19:48918115-48918340GRIN2D0.331.000.380.96
chr21:28337856-28340237ADAMTS50.291.000.310.99
chr22:33453892-33454505SYN30.280.970.320.97
AUC: Area under curve; CGI: CpG island.
Table
Table 2. GRIA4 and VIPR2 methylation analyses results.
Table 2. GRIA4 and VIPR2 methylation analyses results.
GRIA4VIPR2
Tumour Tissue SampleStool Sample MethyLightStool Sample ddPCRTumour Tissue SampleStool Sample MethyLightStool Sample ddPCR
CRC_2Hyper methylatedMethylatedMethylatedHyper methylatedMethylatedMethylated
CRC_3Hyper methylatedMethylatedMethylatedHyper methylatedMethylatedMethylated
CRC_8Undetectable methylationUndetectable methylationMethylatedUndetectable methylationMethylatedMethylated
CRC_12Undetectable methylationUndetectable methylationMethylatedHyper methylatedUndetectable methylationMethylated
CRC_14Not differentially methylatedUndetectable methylationMethylatedHyper methylatedMethylatedMethylated
CRC_19Hyper methylatedMethylatedMethylatedNot differentially methylatedMethylatedMethylated
CRC_21Hyper methylatedUndetectable methylationUndetectable methylationHyper methylatedUndetectable methylationUndetectable methylation
CRC_29Not differentially methylatedMethylatedMethylatedHyper methylatedMethylatedMethylated
CRC_33Hyper methylatedUndetectable methylationMethylatedHyper methylatedMethylatedMethylated
CRC_34Hyper methylatedUndetectable methylationMethylatedNot differentially methylatedUndetectable methylationMethylated
CRC: Colorectal cancer.
Table
Table 3. Clinical characteristics of CRC patients.
Table 3. Clinical characteristics of CRC patients.
Sample IDTumour LocationStage at DiagnosisMucinous HistologyLymphovascular InvasionGradeUlcerative Neoplasia
CRC_2Left colonINONOG2NO
CRC_3Right colonIIINOYESG2YES
CRC_8RectumIVYESYESG2NO
CRC_12Right colon0NOYESG1NO
CRC_14RectumIIIYESYESG3NO
CRC_19Right colonIIYESYESG2YES
CRC_21Transversal colonIINOYESG2NO
CRC_29Right colonIINOYESG2NA
CRC_33Right colonIVNOYESG2NA
CRC_34RectumIIINOYESG2YES
CRC: Colorectal cancer.
Table
Table 4. Primers and probe sequences for MethyLight and ddPCR assay.
Table 4. Primers and probe sequences for MethyLight and ddPCR assay.
TargetForward Primer (5′–3′)Reverse Primer (5′–3′)Probe (5′–3′)
GRIA4GGGTTGGTGTAGGTTTGTTCTCCCCCCTTACTTTCTCACATACACACAAAACGCCGCGACCGCCACAC
VIPR2TCGGTTTCGAGTAGAGAGAATTGGAAACAAATACAAACGACCGCAAAACCCTTCCGAACGCACACCTAACCC
AluGGTTAGGTATAGTGGTTTATATTTGTAATTTTAGATATTAACTAAACTAATCTTAAACTCCTAACCTCACCTACCTTAACCTCCC
Table
Table 5. Primers for quantitative real-time polymerase chain reaction (qRT-PCR) assay.
Table 5. Primers for quantitative real-time polymerase chain reaction (qRT-PCR) assay.
GeneForward Primer (5′–3′)Reverse Primer (5′–3′)
CRC
gene expression assay		GRIA4TCATGTGGACAACATTGAGACAATCATAGAGTCCAAAAATGGCAAA
VIPR2GTCTCTTGCAACAGGAAGCATCTCAGGATGAAGGACAGGAA
SLC6A3CCATACTGAAAGGTGTGGGCTAGAAGAGATAGTGCAGCGCC
SPOCK1AGGTAAAATGCAGCCCTCACATTCCCCTTCTTTTGCCTGGG
TFRCGGCACAGCTCTCCTATTGAAACCAAAGTCTCCAGCACTCCAACT
CRC: Colorectal cancer.

Share and Cite

MDPI and ACS Style

Vega-Benedetti, A.F.; Loi, E.; Moi, L.; Orrù, S.; Ziranu, P.; Pretta, A.; Lai, E.; Puzzoni, M.; Ciccone, L.; Casadei-Gardini, A.; et al. Colorectal Cancer Early Detection in Stool Samples Tracing CpG Islands Methylation Alterations Affecting Gene Expression. Int. J. Mol. Sci. 2020, 21, 4494. https://doi.org/10.3390/ijms21124494

AMA Style

Vega-Benedetti AF, Loi E, Moi L, Orrù S, Ziranu P, Pretta A, Lai E, Puzzoni M, Ciccone L, Casadei-Gardini A, et al. Colorectal Cancer Early Detection in Stool Samples Tracing CpG Islands Methylation Alterations Affecting Gene Expression. International Journal of Molecular Sciences. 2020; 21(12):4494. https://doi.org/10.3390/ijms21124494

Chicago/Turabian Style

Vega-Benedetti, Ana Florencia, Eleonora Loi, Loredana Moi, Sandra Orrù, Pina Ziranu, Andrea Pretta, Eleonora Lai, Marco Puzzoni, Letizia Ciccone, Andrea Casadei-Gardini, and et al. 2020. "Colorectal Cancer Early Detection in Stool Samples Tracing CpG Islands Methylation Alterations Affecting Gene Expression" International Journal of Molecular Sciences 21, no. 12: 4494. https://doi.org/10.3390/ijms21124494

APA Style

Vega-Benedetti, A. F., Loi, E., Moi, L., Orrù, S., Ziranu, P., Pretta, A., Lai, E., Puzzoni, M., Ciccone, L., Casadei-Gardini, A., Cabras, F., Fortunato, F., Restivo, A., Zorcolo, L., Scartozzi, M., & Zavattari, P. (2020). Colorectal Cancer Early Detection in Stool Samples Tracing CpG Islands Methylation Alterations Affecting Gene Expression. International Journal of Molecular Sciences, 21(12), 4494. https://doi.org/10.3390/ijms21124494

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