Abstract:
A high-throughput sequencing method for detecting methylated CpG islands includes: processing a DNA sample by using a modifier, and converting cytosine in the DNA sample into uracil, and keeping 5′methylcytosine unchanged; amplifying the obtained segment by using a primer A and DNA polymerase, to obtain a segment having one end being capable of anchoring a junction primer C; amplifying the obtained segment by using a primer B and DNA polymerase, to obtain a segment gathering methylated CpG islands and having two ends being capable of separately anchoring junction primers C and D; amplifying the obtained segment at a PCR exponent by using the junction primers C and D and the DNA polymerase, to obtain the amplified product; and separating and purifying the amplified product, to form a high-throughput sequencing library and perform computer sequencing, and data analysis.

Description:
BACKGROUND OF THE INVENTION 
       [0001]    1. Field of the Invention 
         [0002]    The present invention relates to detection of methylated CpG islands in a genome using high-throughput sequencing. 
         [0003]    2. Description of Related Art 
         [0004]    DNA methylation is the modification of the cytosine (C) to the 5′-methylated-cytosine (5mC) by adding a methyl group to the C5 position of the cytosine. It plays important roles in many biological processes, including regulation of transcription, repression of transposon, genomic imprinting and X-chromosome inactivation, and is a hot topic in the field of molecular biology. In vertebrates including human, DNA methylation mainly occurs at the CpG site (CpG indicates dinucleotide of which the guanine (G) base immediately follows the cytosine base along the DNA strand). The average content of CpG dinucleotides in the vertebrate genome is lower than their expected frequency. However, in some region of the genome, CpG dinucleotides show the expected or even higher frequency and these regions are referred to as CpG islands. CpG islands are mainly found in gene promoters. In the human genome, there are about 30,000 CpG islands, with more than 50% of CpG islands being located in promotes and more than 60% of promoters contain a CpG island. Methylation of the promoter CpG islands leads to silencing of gene expression; this mechanism participates in many biological processes, including X-chromosome inactivation, genomic imprinting, differentiation of embryonic stem cells, development of germ cells, as well as initiation and progressing of cancers. Intragenic and intergenic CpG islands may be unidentified promoters. Comprehensive understanding the biological functions of CpG islands methylation requires systematic and high-efficient techniques. 
         [0005]    Traditional techniques for detecting DNA methylation, including restriction enzyme digestion, restriction enzyme digestion and PCR, methylation-specific PCR, can only detect single or a small number of sites. With the progresses in the high-throughput sequencing technologies in recent years, researchers begin to systematically profile DNA methylation at the level of whole genome. The current techniques for detecting DNA methylation using high-throughput sequencing including: 1) methylated DNA immuoprecipitation; 2) methylated CpG immuoprecipitation, and 3) bisulfite sequencing. The former two techniques capture methylated DNA by using an antibody or a recombinant methylated CpG binding protein, followed by high-throughput sequencing; they can only measure the DNA methylation status in a semi-quantitative manner and has a resolution of about 100 base-pair (bp). The bisulfite sequencing is based on the fact that sodium bisulfite treatment converts unmethylated cytosines (C) to uracils (U) whereas methylated cytosines (5mC) are not affected; DNA methylation status thus can be acquired by subsequent high-throughput sequencing. This technique, which is the gold standard of DNA methylation analysis, has a resolution as accurate as a single base pair. In 2009, scientists reported the first human whole-genome DNA methylation map at single-base-pair resolution by using the bisulfite sequencing. However, since this technique sequences the whole genome, the cost is very high, which hampers its application to a large number of samples. The reduced representation bisulfite sequencing technique has been developed (Gu H, et al. Preparation of reduced representation bisulfite sequencing libraries for genome-scale DNA methylation profiling. Nat Protoc. 2011 6(4):468-81.). This method enriches genomic regions of promoters and the CpG islands by means of MspI enzyme digestion and gel purification, followed by end-repair, A tailing, adapter ligation, size selection and PCR amplification to obtain the library. Though it is more cost-effective than the whole-genome bisulfite sequencing, the process of library establishment is complex and requires five to six days. In addition, the enrichment process is not able to distinguish between the methylated and unmethylated CpG islands, which increase sequencing cost. Another method (patent: methods and application for establishing library of high-throughput sequencing, CN103103624A) captures genomic regions including the CpG islands by using specific probes and then performs bisulfite sequencing. However, the process for capturing and library establishment is rather time-consuming. 
         [0006]    Therefore, techniques for high-efficient detection of the methylated CpG islands remain limited. 
       SUMMARY OF THE INVENTION 
       [0007]    The present invention enriches and amplifies the methylated CpG islands from a bisulfite-converted genomic DNA sample by means of a high-CpG-density oligonucleotide primer, simultaneously anchoring the adapter sequences to the amplicons for amplification and then sequencing via a three-steps PCR reaction. This technique allows for highly efficient enrichment of the methylated CpG islands and rapid establishment of the high-throughput sequencing library, representing a novel, high-efficient and cost-effective method for detecting the methylated CpG islands using high-throughput sequencing. 
         [0008]    In a first aspect, the present invention provides a method for detecting the methylated CpG islands using the high-throughput sequencing, comprising the following steps: 
         [0009]    (a) Treating a DNA sample with a modifying agent to form the modified DNA wherein cytosine bases but not 5′-methyl-cytosine bases of the DNA sample are modified to uracil bases. 
         [0010]    The DNA sample can be any suitable DNA polymer. Preferably, the DNA sample is the genomic DNA of animals, plants or microorganisms; more preferably the DNA sample is the genomic DNA of human or mouse. 
         [0011]    The DNA sample can derive from human cells, tissues, blood, body fluid, urine, excrement or their combination. In one preferred embodiment, the DNA sample is the cell-free DNA of human plasma or serum. Optionally, the DNA sample can also be the genomic DNA of human whole blood cells or human cancer cell lines. 
         [0012]    The agent for treating the DNA sample modifies cytosine bases but not 5′-methyl-cytosine bases under conditions to form single-stranded DNA. The modifying agent can be selected from bisulfite, acetate or citrate. Preferably, the agent is bisulfite. Optionally, bisulfite treatment of the DNA sample can be achieved by using commercial kits such as MethylCode Bisulfite Conversion Kit (Invitrogen) EZ DNA methylation-Gold Kit (ZYMO) or EpiTect Bisulfite Kit (Qiagen). 
         [0013]    (b) Providing Primer A and DNA polymerase to the modified DNA to allow at least one round of linear amplification to form the semi-amplicon capable of anchoring Adapter Primer C at one end. 
         [0014]    The primer A is composed of a 3′ portion and a 5′ portion, wherein the 3′ portion contains 4 or more nucleotides capable of binding to the modified DNA and allowing amplification. Preferably, the 3′ portion contains only C, A and T except for the CpG dinucleotide. More preferably, the second nucleotide at the 3′ terminal is C. Preferably, the 3′ portion allows preliminary enrichment of the methylated CpG islands, thereby forming the semi-amplicon enriched with the methylated CpG islands and capable of anchoring Adapter Primer C at one end in step (b). The preliminary enriching process, while still need further enriching process in step (c), can increase the net efficiency for enrichment of the methylated CpG islands. To achieve the preliminary enrichment, the 3′ portion of primer A needs to contain a certain density of CpG or C. Preferably, the 3′ portion of primer A is a moderate-CpG-density sequence of which the 3′-terminal 7 nucleotides contain 1 CpG dinucleotide; alternatively, the 3′ portion of primer A is a high-C-density sequence of which the 3′-terminal 5 nucleotides contain 3 or more C, preferably all the 3′-terminal 3 nucleotides are C. Optionally the 3′ portion of primer A is a high-CpG-density sequence of which the 3′-terminal 7 nucleotides contain 2 or 3 CpG dinucleotides, here further enriching process in step (c) is not indispensible (see below). 
         [0015]    The 5′ portion of primer A is used to anchor Adapter Primer C; it allows Adapter Primer C to bind to its reverse complementary sequence for PCR amplification. The word “anchor” is to describe the function of the 5′ portion to join Adapter Primer C to the amplicon via PCR. Preferably, the sequences of the 5′-terminal 15 to 40 nucleotides of Prime A are the same as that of the 3′-terminal 15 to 40 nucleotides of Adapter Primer C. 
         [0016]    The DNA polymerase can be any suitable polymerase, such as Taq polymerase, ExTaq polymerase, LATaq polymerase, AmpliTaq polymerase, Amplitaq Gold polymerase, Titanium Taq polymerase, KlenTaq polymerase, Platinum Taq polymerase, Accuprime Taq polymerase, Pyrobest DNA polymerase, Pfu polymerase, Pfu polymerase turbo, Phusion polymerase, Pwo polymerase, Vent polymerase, Vent Exo-polymerase, Sequenase™ polymerase, 9° Nm DNA polymerase, Therminator DNA polymerase, Expand DNA polymerase, rTth DNA polymerase, DyNazyme™ EXT polymerase, DNA polymerase I, T7 DNA polymerase, T4 DNA polymerase, Bst DNA polymerase, phi-29 DNA polymerase, and Klenow fragment. 
         [0017]    Preferably, the DNA polymerase is capable of strand displacement. It can be any suitable polymerase with the strand displacement activity, including but not restricted to DNA polymerase I (Klenow) large fragment (New England Biolabs (NEB) cateloge: M0210S), Klenow fragment (exo-) (NEB catalogue: M0212S), Bst DNA polymerase large fragment (NEB cateloge: M0275S), Vent(exo-) (NEB cateloge: M0257S), Deep Vent (exo-) (NEB cateloge: M0259S), M-MulV reverse transcriptase (NEB cateloge: M0253S), 9 Nm DNA polymerase (NEB cateloge: M0260S) and Phi-29 DNA polymerase (NEB cateloge: M0269S). In one preferred embodiment, the DNA polymerase is Klenow fragment (exo-). 
         [0018]    Preferably, the DNA polymerase is deficient of exonuclease activity. 
         [0019]    Linear amplification refers to that the amount of amplified product increases in a linear instead of an exponential relationship to the amplification times. The step needs at least one round of linear amplification; preferably, 2-30 rounds of linear amplification are performed. 
         [0020]    (c) Amplifying the semi-amplicon by using Primer B and DNA polymerase to form the full-amplicon enriched with methylated CpG islands and capable of anchoring Adapter Primer C at one end and Adapter Primer D at the other end; 
         [0021]    Primer B is composed of a 3′ portion and a 5′ portion. The 3′ portion contains 7 or more nucleotides and is a high-CpG-density sequence capable of binding to the semi-amplicon and allowing amplification and enrichment of the methylated CpG islands, wherein the high-CpG-density sequence is a sequence of which the 3′-terminal 7 nucleotides contain 2 or 3 CpG dinucleotides. Preferably, the 3′ portion of Primer B contains only A and T except for the CpG dinucleotides. 
         [0022]    Enrichment of the methylated CpG islands is due to the preference of Primer B to binding to the methylated CpG islands, while not binding to the regions outside the methylated CpG islands, in a genome. This preference is determined by the high-CpG-density nature of the 3′ portion of the primer. Our bioinformatics analysis demonstrates that the high-CpG-density sequences are not randomly distributed in the genome, but accumulate at the CpG islands to a certain degree. Hence, a high-CpG-density primer is capable of enriching the methylated CpG islands. Preferably, the 3′ portion is also a high-GC-content sequence, which indicates that the total number of C and G bases within the 3′-terminal 10 nucleotides of the primer is equal to or greater than 7. A high-GC-content sequence will not only more accumulate in the CpG islands, but also help additionally increasing the annealing temperature of the primer. Technical personnel in this field can understand that the annealing temperature of a primer is positively correlated with its GC content. The according high-GC-content feature of the primer makes its annealing temperature reaching about 40 to 60 centigrade, thus facilitating binding of the primer to the template and increasing the amplification efficiency. Preferably, the 3′ portion of Primer B is an extremely-high-CpG-density sequence, of which the 3′-terminal 7 nucleotides of the primer contains 3 CpG dinucleotides. 
         [0023]    The 5′ portion of Primer B is used to anchor Primer D. it allows Adapter Primer D to bind to its reverse complementary sequence for PCR amplification. The word “anchor” is to describe the function of the 5′ portion to join Adapter Primer D to the amplicon via PCR. Preferably, the sequence of the 5′ portion of Primer B is the same as the 3′-terminal 15 to 40 nucleotides of Primer D. 
         [0024]    The polymerase can be any suitable polymerase mentioned before. Preferably, the polymerase is hot-start. Preferably, the annealing temperature of the reaction is between about 40 and 60 centigrade. 
         [0025]    When the 3′ portion of Primer A has a character of high-CpG-density, the semi-amplicons in the step (c) have already enriched the methylated CpG islands (and been able to anchor Adapter Primer C at one end). Thus, the full-amplicons, which enrich the methylated CpG islands and are able to anchor Adapter primer C at one end and Adapter primer D at the other ends, can be obtained by means of Primer B whose 3′ portion is just able to bind and amplify the semi-amplicons without a further function of CpG island enrichment. 
         [0026]    (d) Amplifying the full-amplicon by using Adapter Primer C, Adapter Primer D and DNA polymerase to form the final-amplicon via PCR exponential amplification. 
         [0027]    The word “Adapter Primer” herein refers to that the function of the primer is similar to the “adapter” used in conventional methods for establishing the high-throughput sequencing library (such as Illumnia&#39;s TruSeq DNA Sample Prep Kit and Applied Biosystems (ABI)&#39;s The SOLiD™ Fragment Library Construction Kit), which allows binding of the DNA fragments to the flow cell for subsequent amplification and sequencing. Different from the usage of ligation reaction to add “adapter” in conventional methods for establishing high-throughput sequencing library, the present invention adds “adapters” to each end of the amplicons by means of PCR reaction via the anchor sequences at the 5′ portion of Primers A and B. Technical personnel in this field can understand that either Primer C or Primer D can contain “barcode” sequences, which facilitate to simultaneously sequence multiple libraries in one flow cell. The Primer C and D correspond to the adapter sequences of a given high-throughput sequencing platform, which includes, but does not limited to, the Illumina&#39;s Genome Analyzer IIx, HiSeq and MiSeq platforms, ABI&#39;s SoLiD, 5500 W Series Genetic Analyzer and Ion Torrent PGM platforms, and Roche454&#39;s GS Junior and GS FLX+ platforms. 
         [0028]    (e) Separating and purifying the final-amplicon to form the library for high-throughput sequencing, sequencing the library and analyzing the data. 
         [0029]    The approach to separate and purify the final amplicons can be any suitable method, including but not restricted to magnetic beads-based, column-based and gel electrophoresis-based purification. Preferably, the purification method is able to make size-selection for the amplicons. Preferably, the amplicons between 160 and 400 bp are selected. Preferably, the amplicons are separated by using 2% agarose gel electrophoresis and the fragments between 160 and 400 bp are excised and then purified. 
         [0030]    The high-throughput sequencing platform for analysis includes, but does not limited to, the Illumina&#39;s Genome Analyzer IIx, HiSeq and MiSeq platforms, ABI&#39;s SoLiD, 5500 W Series Genetic Analyzer and Ion Torrent PGM platforms, and Roche454&#39;s GS Junior and GS FLX+ platforms. 
         [0031]    The method for data analysis is not limited and can be any suitable software for data analysis and sequence alignment, which includes, but does not limited to Bismark, BSMAP, Bowtie and SOAP. 
         [0032]    In a second aspect, the present invention provides a kit for detecting methylated CpG islands using high-throughput sequencing. The kit comprises a set of primers including the Primer A, Primer B, Adapter Primer C, Adapter Primer D, and the DNA polymerases according to the first aspect, as well as instructions for the kit. 
     
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         [0033]      FIG. 1  Schematic of the present invention for detecting methylated CpG islands based on the high-throughput sequencing (SEQ ID NOS: 10-20). 
           [0034]      FIG. 2  Agarose gel electrophoresis results of the present invention for amplification of the genomic DNAs extracted from HeLa cells and human peripheral blood mononuclear cells (hPBMCs). 
           [0035]      FIG. 3  Bioanalyzer 2100 results showing the size distribution of the libraries established from the genomic DNAs of HeLa cells and hPBMCs by the present invention. 
           [0036]      FIG. 4  IGV views of the present invention in comparison with the whole-genome bisulfite sequencing for the genomic regions containing MAEL, ILDR2 or CDKN2A. 
           [0037]      FIG. 5  The percentage of a number of short nucleotide sequences containing 1 to 3 CpG dinucleotides located within versus outside the CpG islands in the human genome. 
           [0038]      FIG. 6  Enrichment of the CpG islands by means of Primer A and Primer B with various 3′ portion sequences. 
       
    
    
     DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
       [0039]    In order that the present invention may be more clearly understood, preferred forms will described with reference to the following examples and drawings. It will be appreciated by persons skilled in the art that numerous variations and/or modification may be made to the invention as shown in the specific embodiments without departing from the spirit or scope of the invention as broadly described. 
         [0040]    The present invention enriches and amplifies the methylated CpG islands from a bisulfite-converted genomic DNA sample by means of a high-CpG-density oligonucleotide primer, simultaneously anchoring the adapter sequences to the amplicons for high-throughput sequencing via a three-steps PCR reaction. An example of the present invention is shown below ( FIG. 1 ). 
         [0041]    1. Bisulfite Conversion of the DNA Sample. 
         [0042]    The bisulfite conversion is performed by using the MethylCode™ Bisulfite Conversion Kit (Invitrogen) according to the protocol provided by the manufacturer. Detailed steps are as follows: 
         [0000]    1.1 Preparing the CT conversion reagent: Add 900 μl of sterile distilled water, 50 μl of resuspension buffer, and 300 μl of dilution buffer directly to one tube of CT conversion reagent; mix by shaking or intermittent brief vortexing for 10 minutes for dissolvation, keep protected from light at room temperature;
 
1.2 Add 20 μl of the DNA sample ranging from 500 pg to 500 ng to a PCR tube;
 
1.3 Add 130 μl of CT conversion reagent to the DNA sample, and mix by flicking the tube or pipetting up and down;
 
1.4 Place the tube in a thermal cycler and run the following program: 98° C. for 10 minutes, 64° C. for 2.5 hours and 4° C. storage for up to 20 hours;
 
1.5 Place a spin column in a collection tube and add 600 μl of binding buffer to the column;
 
1.6 Add the sample from step 1.4 to the binding buffer in the column; close the cap and mix by inverting several times;
 
1.7 Centrifuge at full speed (&gt;10,000×g) for 30 seconds, then discard the flow-through;
 
1.8 Add 100 μl of wash buffer prepared with ethanol to the column, and centrifuge at full speed for 30 seconds; discard the flow-through;
 
1.9 Add 200 μl of desulphonation buffer to the column and let the column stand at room temperature for 15-20 minutes;
 
1.10 Centrifuge at full speed for 30 seconds; discard the flow-through;
 
1.11 Add 100 μl of wash buffer with ethanol to the column and centrifuge at full speed for 30 seconds; discard the flow-through;
 
1.12 Repeat the wash in step 1.11 one more time, then transfer the spin column to a new, clean 1.5-ml microcentrifuge tube;
 
1.13 Add 10 μl of elution buffer directly to the column matrix; centrifuge at full speed for 30 seconds to elute the DNA.
 
         [0043]    2. Lineal Amplification by Means of Primer a and DNA Polymerase. 
         [0000]    2.1 Prepare the following reaction mixture in a PCR tube for the DNA sample obtained from step 1.13: 
         [0000]                                        Contents   Volume                       The DNA sample (from   4 μl           step 1.13)               NEBuffer 2   1.5 μl             dNTP (2.5 mM)   1.5 μl             Primer A (5 uM) *   1 μl           Water   6 μl           Klenow (exo-) **   0.3 μl × 5 (times, see                step 2.5) = 1.5 μl           Total   15 μl                                          
A/T/C), the 3′ portion of Primer A is underlined by a wave line ’ and the 5′ portion is underlined by a straight line ’ showed as SEQ ID No: 1 in the SEQUENCE LISTING       **: add at step 2.3            
2.2 Place the tube in a thermal cycler and run the following program: 95° C. for 2 minutes and 4° C. for storage;
 
2.3 Add 0.3 μl of Klenow Fragment (exo-) (NEB catalog: M0212S), mix and briefly centrifuge;
 
2.4 Place the tube in a thermal cycler and run the following program: 4° C. for 50 seconds, 10° C. for 50 seconds, 20° C. for 50 seconds, 30° C. for 50 seconds, 37° C. for 5 minutes, 95° C. for 10 seconds, 4° C. for pause;
 
2.5 Repeat step 2.3 and 2.4 for a total of 4 times, omit the step of 4° C. for pause at the last time;
 
2.6 Run the following program to inactivate the Klenow Fragment: 75° C. for 20 minutes and 50° C. for pause.
 
         [0044]    3. Amplification by Means of Primer B and DNA Polymerase. 
         [0000]    3.1 Prepare the following reaction mixture in a new PCR tube 
         [0000]                                Contents   Volume                   Ex Taq Buffer   0.5 μl       Primer B (10 uM) *   0.5 μl       HS Ex Taq (Takara catalog: DRR006B)   0.3 μl       Water   3.7 μl       Total     5 μl                                         
A/T/G); the 3′ portion is underlined by a wave line, and the 5′ portion is underlined by a straight line ’showed as SEQ ID No: 2 in the SEQUENCE LISTING            
3.2 Place the tube in a thermal cycler and preheat to 50° C.;
 
3.3 Add the preheated mixture to the products (the semi-amplicons) from step 2.6, mix by pipetting up and down 5-6 times;
 
3.4 Place the tube in a thermal cycler and run the following program: 95° C. for 3 minutes, 50° C. for 1 minutes, 72° C. for 1 minutes, 50° C. for pause.
 
         [0045]    4. Exponential Amplification Using Adapter Primer C, Adapter Primer D and DNA Polymerase. 
         [0000]    4.1 Prepare the following reaction mixture in a new PCR tube: 
         [0000]                                Contents   Volume                   Ex Taq Buffer    12 μl               Adapter Primer C (100 uM) *     2.4 μl               Adapter Primer D (100 uM) **     2.4 μl               dNTP (2.5 mM)    12 μl               HS Ex Taq (Takara catalog: DRR006B)     1.2 μl               Water    90.8 μl               Total   120 μl               *: Adapter Primer C:       AATGATACGGCGACCACCGAGATCTACAC TCTTTCCCTACACGACGCTCTTCCGATCT ; the underlined oligonucleotide sequence is the same as the 5′ portion of Primer A ′ showed as SEQ ID No: 3 in the SEQUENCE LISTING;       **: Adapter Primer D:       CAAGCAGAAGACGGCATACGAGAT CTGATC GTGAC TGGAGTTCAGACGTGTGCTCTTCCGATCT , the underlined (by a straight line) oligonucleotide sequence is the same as the 5′ portion of Primer B and the barcode sequence (corresponding to the Illumina index 9) is       underlined by a double-line ′ showed as SEQ ID No: 4 in the SEQUENCE LISTING;            
4.2 Preheat the mixture at 50° C. in a thermal cycler;
 
4.3 Add the preheated mixture to the product (full-amplicons) from step 3.4, mix by pipetting up and down 5-6 times;
 
4.4 Place the tube in a thermal cycler and run the following program: 95° C. for 3 minutes;
 
4.5 Run the following program: 20 cycles of 95° C. for 30 seconds, 67° C. for 30 seconds, 72° C. for 1 minute, and then 4° C. for storage.
 
         [0046]    5. Size-Selection, Purification, High Throughput Sequencing and Data Analysis. 
         [0000]    5.1 Prepare 2% agarose gel, add 1×SYBR Safe (Invitrogen);
 
5.2 Run the products (final-amplicons) from step 4.5 in the agarose gel by electrophosis;
 
5.3 Imagine the gel;
 
5.4 Excise the gel for DNA fragments ranging from 160 to 400 bp;
 
5.5 Purify the DNA from the gel by using the QIAquick Gel Extraction Kit (Qiagen) according to the manufacturer&#39;s protocol, obtaining the library for high throughput sequencing;
 
5.6 Analyze the size distribution of the library by using the Bioanalyzer 2100 (Agilent) and quantify the concentration of the library by using QPCR;
 
5.7 Sequence the library on an Illumina HiSeq2000 sequencer at pair-end 100 bp reads to obtain the raw sequencing data;
 
5.8 Data analysis: firstly discard whole or any subsets of adapter sequences and low quality reads, then map the reads to the human genome reference (hg19), and then perform subsequent bioinformatics analysis.
 
         [0047]      FIG. 2  shows 2% agarose gel electrophosis results of the products amplified from the genomic DNAs of HeLa cells and hPBMCs by using the present invention. Lane 1 and 4 are DNA markers; lane 2 is the result of the HeLa cells from 15 ng of bisulfite-converted DNA; lane 3 is the result of the hPBMCs (from an adult male) from 30 ng of bisulfite-converted DNA; lane 5 is the non-template control.  FIG. 2  shows that the genomic DNAs of both the HeLa cells and hPBMCs are positively amplified. These results demonstrate that the present method, which is able to work with as low as 15 ng of genomic DNA as the starting material, is highly efficient and sensitive. 
         [0048]      FIG. 3  shows the Bioanalyzer 2100 results of the libraries established from the genomic DNAs of HeLa cells and hPBMCs by using the present invention. The libraries are made from the amplified products mentioned in  FIG. 2  by excision and purification of the 2% agarose gel after electrophosis.  FIG. 3A  show the result of the HeLa cells and  FIG. 3B  show the result of the hPBMCs. The results show that the fragment size of the libraries ranges from 160 to 280 bp. 
         [0049]    The two libraries are sequenced on the HiSeq2000 platform, obtaining 1.3 and 1.5 gigabase (Gb) raw data for the HeLa cells and the hPBMCs, respectively. The statistical analysis of the data is shown in Table 1. The CpG islands comprise about 0.7% of the whole human genome; in contrast, in the sequencing data of the present invention for the HeLa cells and the hPBMCs, 39% and 20% of the reads map to the CpG islands, respectively. Thus, the present invention enriches the CpG islands for 56 and 29 times for the HeLa cells and the hPBMCs, respectively. The present invention achieves an average depth of 20-30× for the genomic regions including the methylated CpG islands by sequencing only 1˜2 Gb raw data. In contrast, the whole genome bisulfate sequencing requires 150˜200 Gb raw data to achieve a similar depth. This demonstrates that the present invention greatly increases the efficiency for detecting the methylated CpG islands based on the high-throughput sequencing. 
         [0000]    
       
         
               
             
               
               
               
               
               
               
               
               
               
               
               
               
             
               
               
               
               
               
               
               
               
               
               
               
               
             
           
               
                 TABLE 1 
               
             
             
               
                   
               
               
                 Data analysis of the present invention for the genomic DNAs of the HeLa cells and the hPBMCs. 
               
             
          
           
               
                   
                   
                   
                   
                   
                   
                   
                 Reads 
                   
                   
                   
                   
               
               
                   
                   
                   
                   
                 After 
                 Uniquely 
                   
                 overlapping 
                 Percentage 
                 Enrichment 
               
               
                   
                 Raw 
                 Raw 
                 Clean 
                 adapter 
                 mapped 
                   
                 the CpG 
                 of reads 
                 of the CpG 
                   
                 Average 
               
               
                   
                 data 
                 data 
                 data 
                 trimming 
                 reads 
                   
                 islands 
                 overlapping 
                 islands 
                 Number 
                 sequencing 
               
               
                   
                 (reads 
                 (Bases 
                 (Reads 
                 (Reads 
                 (Reads 
                 Mapping 
                 (Reads 
                 the CpG 
                 (observed/ 
                 of CpG 
                 depth for 
               
               
                   
                 (M)) 
                 (Gb)) 
                 (M)) 
                 (M)) 
                 (M)) 
                 efficiency 
                 (M)) 
                 islands 
                 expected) 
                 (depth &gt;= 3) 
                 each CpG 
               
               
                   
                   
               
             
          
           
               
                 HeLa 
                 12.9 
                 1.3 
                 7.5 
                 6.4 
                 3.6 
                 57% 
                 1.4 
                 39% 
                 56 
                 495559 
                 27× 
               
               
                 hPBMC 
                 14.7 
                 1.5 
                 8.1 
                 5.9 
                 3.1 
                 53% 
                 0.6 
                 20% 
                 29 
                 435131 
                 21× 
               
               
                   
               
             
          
         
       
     
         [0050]      FIG. 4  shows IGV views of the present invention in comparison with the whole-genome bisulfite sequencing for the genomic regions containing MAEL, ILDR2 or CDKN2A.  FIG. 4A  shows the genomic region containing MAEL and ILDR2. MAEL (maelstrom spermatogenic transposon silencer) is a gene specially expressed in the testis. The promoter CpG islands of this gene are demethylated in the germ cell while highly methylated in the somatic cell. On the contrary, the promoter CpG island of ILDR2 (immunoglobulin-like domain containing receptor 2) is demethylated in the somatic cell. For the genomic DNAs of both the HeLa cells and the hPBMCs, the promoter CpG island of MAEL is amplified by the present invention; the sequencing results shows that it is highly methylated in both genomes. In contrast, the promoter CpG island of ILDR2 is not amplified for both genomic DNAs. These results demonstrate that the present invention is able to selectively amplify the methylated CpG islands while does not amplify the demethylated CpG islands. The whole-genome bisulfite sequencing (of the human brain tissue) verifies that, in the normal somatic cell, the promoter CpG island of MAEL is highly methylated while that of MAEL is not methylated. On the other hand, it demonstrates that most isolated CpG sites are highly methylated, and they are selectively discarded during the process of the present invention for enriching the methylated CpG islands and establishing the library.  FIG. 4B  shows the genomic region containing CDKN2A (cyclin-dependent kinase inhibitor 2A). The CDKN2A genomic region contains several CpG islands being demethylated in the somatic cell while some of them are methylated in cancer cells. The present invention shows that 2 of 4 CpG islands are amplified in the genomic DNA of the HeLa cells and the sequencing results shows that they are highly methylated, while all 4 CpG islands are not amplified in the genomic DNA of the hPBMCs. These results further demonstrate that the present invention is able to accurately and high-efficiently enrich and amplify the methylated CpG islands for the high throughput sequencing. 
         [0051]      FIG. 5  shows the bioinformatics analysis for the percentage of a number of short nucleotide sequences containing 1 to 3 CpG dinucleotides located within versus outside the CpG islands in the human genome. The results demonstrate that the enrichment of a short nucleotide sequence in the CpG islands is positively correlated with the number of CpG dinucleotide it contains. 
         [0052]      FIG. 6  (SEQ ID NOS: 1, 2, 5-9) shows enrichment of the CpG islands by means of Primer A and Primer B with a variety of 3′ portion sequences for the genomic DNA of the HeLa cells. The results show that the CpG islands enrichment is positively correlated to the CpG density of the primers. H=C/A/T/D=G/A/T R=G/A N=A/T/C/G. 
         [0053]    The above description of exemplary embodiments of the invention has been presented for the purposes of illustration and description. It is not intended to be exhaustive or to limit the invention to the precise form described, and many modifications and variations are possible in light of the teaching above. The embodiments were chosen and described in order to best explain the principles of the invention and its practical applications to thereby enable others skilled in the art to best utilize the invention in various embodiments and with various modifications as are suited to the particular use contemplated.