Epigenetic modifications, including DNA methylation, are heritable changes essential for normal development via the control of gene expression (1). These heritable changes are established during the cell cycle and cell division, leading to distinct identities of gene expression profiles while the primary sequences of DNA are the same. Failure in the proper maintenance of epigenetic machinery can lead to an alteration of gene expression and subsequent carcinogenesis (2). Recent evidence in the field of epigenetic analysis showed that cancer cells exhibit global epigenetic alterations in addition to various genetic changes observed during malignant transformation (3). DNA methylation is one of the most extensively studied epigenetic modifications characterized as the methylation of C5 position of cytosine, resulting in its modification into 5-methylcytosine. DNA methylation occurs mainly at the cytosine-phosphate-guanine (CpG) sites, which are localized in the promoter region of tumor suppressor genes. The transcription of genes initiates from the promoter regions characterized by relatively high frequencies of CpG dinucleotides. For this reason, aberrantly methylated CpG islands are currently defined as promising prognostic and predictive markers associated with cancer development (4).
Meis homeobox 1 (
Selecting the appropriate methods for the evaluation of gene promoter methylation status based on their capacity to distinguish 5-methylcytosine and nonmethylated CpG sites is the most crucial step toward subsequent implementation of the detecting method into clinical practice, particularly without biases regarding false positive and false negative results. Currently, several innovative methods and high-throughput molecular biology and genomics technologies are available to detect the degree of DNA methylation status (10). PCR-based methods have been widely used to evaluate DNA methylation status for decades. The first generation of PCR used methylation-specific-restricted enzymes to distinguish methylated and unmethylated alleles. These methods were followed by methylation-specific PCR (MSP) based on the amplification of sodium bisulfite-treated DNA and consecutively analyzing the presence or absence of converted unmethylated cytosine to uracil after the process of deamination (11). Droplet digital PCR (ddPCR) is a relatively novel method for estimating the absolute quantification of the target sequence of DNA. Recent data show that ddPCR manifests higher sensitivity and specificity than routinely used qPCR (12). In the present article, we describe the optimization of in-house developed ddPCR assay for the identification of the
Designing of ddPCR was conducted by the employment of free on-line application: The Eukaryotic Promoter Database (EPD), MethPrimer, and Primer3.
Firstly, EPD application was used (
Consequently, we used the commercially available software MethPrimer (
The final step included usage of free available primer design software Primer3 (
Primer and probe sequences for
Primers | |
---|---|
Name | Sequence (5´→3´) |
TGGGGAGAGAGTTTGTAGG | |
ACACAAACACCACACACC | |
Name | Sequence |
CGGTCGCGGGTTATTGTTTGC | |
TGGTGGTTGTGGGTTATTGTTTGTGT |
The ddPCR platform (Bio-Rad Laboratories, Hercules, CA, USA) was used for the optimization of methylated (M-Probe, FAM labeled) and unmethylated (UnM-Probe, HEX labeled) probes. As a template sequence, we used commercially available fully methylated and fully unmethylated EpiTect DNA controls (Qiagen, Hilden, Germany). The reaction mix for ddPCR contained 10 μL of Supermix for Probes (No dUTP) (Bio-Rad Laboratories), 0.45 μL of each primer, 0.45 μL of each probe, 1 μL of methylated and unmethylated control DNA adjusted with 7.75 μL of water up to a final volume of 20 μL. Subsequently, 70 μL of Droplet Generation Oil for Probes, along with the reaction mixture, was loaded into a DG8 cartridge and inserted into the QX200 Droplet Generator (Bio-Rad Laboratories) to divide each sample into 20,000 droplets. After generating oil emulsion, approximately 40 μL of the sample was transferred into 96-well PCR plates, covered with a pierceable foil, and heat-sealed by Bio-Rad’s PX1 system. The first optimization step of the ddPCR assay covers estimating optimal PCR conditions. The thermal PCR consisted of enzyme activation for 10 min at 95 °C, followed by 40 cycles of denaturation for 30 sec at 94 °C, annealing/extension step for 1 min with temperature gradient, finished by one cycle of enzyme deactivation for 10 min at 98 °C. To define the optimal conditions of annealing, we performed the above-mentioned temperature gradient in the range of 50 – 62 °C. After finishing the PCR program, we loaded the plate onto the QX200 Droplet Reader (Bio-Rad Laboratories) for a final analysis.
The second optimization step was performing a concentration gradient of eight different dilutions of DNA controls to confirm the dynamic range and linearity of the method. The reaction mix for ddPCR contained 10 μL of Supermix for Probes (No dUTP) (Bio-Rad Laboratories), 0.45 μL of each primer, 0.45 μL of each probe, diluted methylated and unmethylated DNA controls in water with 8000, 4000, 2000, 1000, 500, 250, 125, and 62 copies per reaction adjusted with a variable volume (regarding to DNA input) of water up to a final volume of 20 μL. Further workflow and analysis were the same as described above.
Using commercial methylated and unmethylated DNA, we determined the best annealing temperature to be 57.1 °C for both methylated and unmethylated probes. Both DNA controls were combined with a methylated probe and an unmethylated probe. The methylated probe (Fig. 3A) manifested a sufficient sensitivity in the detection of methylated DNA and, at the same time, a perfect selectivity capacity. The unmethylated probe showed a high sensitivity for detecting unmethylated control DNA, but the same probes also manifested a cross-reactivity resulting in the generation of positive droplets detected in methylated control DNA (approx. 9%) (Fig. 3B).
Eight different dilutions of DNA controls were included in the concentration gradient to confirm the dynamic range and linearity of the method. A methylated (Fig. 4A) and an unmethylated DNA control (Fig. 4B) were used for the dilution. Similar to the concentration gradient, both controls were combined with a methylated probe and an unmethylated probe. The optimal concentration was determined by considering signal amplitude, the overall count of positive droplets, and binding specificity. In addition, thresholds were set at the level of 1500 for FAM and 1400 for HEX.
DNA methylation plays an essential role in cancer initiation, progression, and development of its metastatic form (13). Development and optimization of novel methods for the analysis of specific methylation patterns provides new opportunities and suggests a potential concept for clinical implementation (10). The ddPCR platform allows the precise detection and absolute quantification of targeted DNA sequences. Usage of bisulfite-treated DNA and subsequent application of specific probes for methylated and unmethylated sequences is also suitable for quantifying methylated DNA at single-base resolution utilizing the ddPCR platform (11,14).
In current study, we designed the probes with high specificity and sensitivity targeting the promoter region of the
In conclusion, we described the design and optimization of the ddPCR assay recognizing methylated CpG dinucleotides in the