[go: up one dir, main page]

CN102796808B - Methylation high-flux detection method - Google Patents

Methylation high-flux detection method Download PDF

Info

Publication number
CN102796808B
CN102796808B CN201110133858.1A CN201110133858A CN102796808B CN 102796808 B CN102796808 B CN 102796808B CN 201110133858 A CN201110133858 A CN 201110133858A CN 102796808 B CN102796808 B CN 102796808B
Authority
CN
China
Prior art keywords
sequence
dna
methylation
library
sequencing
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201110133858.1A
Other languages
Chinese (zh)
Other versions
CN102796808A (en
Inventor
罗慧娟
吉冠玉
蒋慧
吴明枝
孙继华
吴仁花
王君文
高飞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
BGI Technology Solutions Co Ltd
Original Assignee
BGI Technology Solutions Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by BGI Technology Solutions Co Ltd filed Critical BGI Technology Solutions Co Ltd
Priority to CN201110133858.1A priority Critical patent/CN102796808B/en
Priority to US14/119,266 priority patent/US9133513B2/en
Priority to EP12789707.2A priority patent/EP2725125B1/en
Priority to DK12789707.2T priority patent/DK2725125T3/en
Priority to PCT/CN2012/075919 priority patent/WO2012159564A1/en
Publication of CN102796808A publication Critical patent/CN102796808A/en
Priority to HK13104998.9A priority patent/HK1178211B/en
Application granted granted Critical
Publication of CN102796808B publication Critical patent/CN102796808B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6869Methods for sequencing
    • C12Q1/6874Methods for sequencing involving nucleic acid arrays, e.g. sequencing by hybridisation
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12NMICROORGANISMS OR ENZYMES; COMPOSITIONS THEREOF; PROPAGATING, PRESERVING, OR MAINTAINING MICROORGANISMS; MUTATION OR GENETIC ENGINEERING; CULTURE MEDIA
    • C12N15/00Mutation or genetic engineering; DNA or RNA concerning genetic engineering, vectors, e.g. plasmids, or their isolation, preparation or purification; Use of hosts therefor
    • C12N15/09Recombinant DNA-technology
    • C12N15/10Processes for the isolation, preparation or purification of DNA or RNA
    • C12N15/1034Isolating an individual clone by screening libraries
    • C12N15/1065Preparation or screening of tagged libraries, e.g. tagged microorganisms by STM-mutagenesis, tagged polynucleotides, gene tags
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6806Preparing nucleic acids for analysis, e.g. for polymerase chain reaction [PCR] assay
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6869Methods for sequencing
    • CCHEMISTRY; METALLURGY
    • C40COMBINATORIAL TECHNOLOGY
    • C40BCOMBINATORIAL CHEMISTRY; LIBRARIES, e.g. CHEMICAL LIBRARIES
    • C40B50/00Methods of creating libraries, e.g. combinatorial synthesis
    • C40B50/06Biochemical methods, e.g. using enzymes or whole viable microorganisms

Landscapes

  • Chemical & Material Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Organic Chemistry (AREA)
  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Wood Science & Technology (AREA)
  • Zoology (AREA)
  • Genetics & Genomics (AREA)
  • Biochemistry (AREA)
  • General Engineering & Computer Science (AREA)
  • Biotechnology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Microbiology (AREA)
  • Molecular Biology (AREA)
  • General Health & Medical Sciences (AREA)
  • Biophysics (AREA)
  • Physics & Mathematics (AREA)
  • Analytical Chemistry (AREA)
  • Immunology (AREA)
  • Biomedical Technology (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Plant Pathology (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • General Chemical & Material Sciences (AREA)
  • Medicinal Chemistry (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)

Abstract

本发明提供一种甲基化高通量检测方法,特别提供了一种结合序列捕获和重亚硫酸盐测序的方法,该方法能够同时准确有效地分析几个样品中的目标区域甲基化状态,降低探针设计的难度,增加操作及应用的可行性,使得对全基因组内感兴趣的目标序列和区域进行高通量高精度的甲基化检测成为现实,并且该方法还具有针对性强以及节约成本和时间的特点。The present invention provides a high-throughput methylation detection method, especially a method combining sequence capture and bisulfite sequencing, which can accurately and effectively analyze the methylation status of target regions in several samples at the same time , reduce the difficulty of probe design, increase the feasibility of operation and application, and make high-throughput and high-precision methylation detection of target sequences and regions of interest in the whole genome a reality, and this method is also highly targeted As well as cost and time saving features.

Description

High-flux detection method methylates
Technical field
The present invention relates to high-throughput genomic methylation DNA detection field, a kind of method that binding sequence is caught and bisulfite (bisulfite) checks order is particularly provided.In addition, the invention still further relates to the order-checking of exon group and the methylation analysis sequencing technologies field based on bisulfite conversion.The invention still further relates to the label that methylates (Index) joint technique, can realize and on a chip, carry out catching of several samples simultaneously.Method of the present invention is specially adapted to s-generation sequencing technologies, especially solexa sequencing technologies.
Background technology
DNA methylation is epigenetics research field one of the direction of focus the most, also becomes just gradually that Mammals grows and the epigenetics mark of the various diseases such as cancer.DNA methylation is not only modified chromatin Structure, genome stability has vital role, and in eukaryote, DNA methylation participates in regulate several biological processes, as fetal development, genomic imprinting, x chromosome inactivation, the adjusting of genetic expression and silence, generation (the Brena RM et al.2006 of the various diseases such as silence and mammal tumor of retrotransposon; Egger G et al.2004; Gu H et al.2010).DNA methylation biomarker is that the early stage assessment of various diseases comprises the detection assessment to high-risk individuality, and a large amount of reference informations are provided.
At present, the order-checking of exon group and the methylation analysis order-checking based on bisulfite conversion are all ripe.Take Illumina Solexa, AB SOLiD and Roche 454 as the s-generation sequencing technologies of representative had been rapidly developed in recent years, become the important tool of genomics research.The feature of s-generation sequencing technologies maximum is high-throughput, it can check order to hundreds of millions of DNA fragmentations simultaneously, a current high-flux sequence instrument once can produce the data up to 300Gb, be equivalent to a people's genome sequencing 100 times, high-flux sequence is, by ultrasonic wave or additive method, genome is broken into a series of small segment, and add top connection (adaptor) in the both sides of small segment, then carry out the fundamental unit of bridge-type PCR or emulsion-based PCR (emulsion PCR) amplification formation order-checking by joint primer, again according to the partial sequence design common sequencing primer on joint, genomic dna is checked order.The method of current research DNA methylation mainly contains: the methylation analysis method of PCR-based, the methylation analysis method based on restriction enzyme, methylation analysis method and column chromatography etc. based on bisulfite.It is that to methylate the most frequently used be also method the most accurately in research that bisulfite is processed in conjunction with high-flux sequence, bisulfite sequencing (bisulfite genomics sequencing) is to utilize bisulfite that the unmethylated cytosine(Cyt) (C) in genomic dna is modified to uridylic (U), methylated cytosine(Cyt) (C) remains unchanged, in pcr amplification, uridylic (U) will become thymus pyrimidine (T), therefore the site that methylates just becomes a common C/T single base polymorphisms site, wherein the frequency of allelotrope cytosine(Cyt) (C) is the degree of gene methylation, it is large that the method has data volume, the shortcoming that cost is high.In conjunction with methylate DNA co-immunoprecipitation (the Methylated DNA Immunoprecipitation in addition of high-flux sequence, MeDIP) and methylated CpG in conjunction with territory (methyl-CpG binding domain, MBD) column chromatography, although they can well carry out enrichment to being rich in the region that methylates, but the relative homogeneous of coverage, lacks specific aim.
Sequence capturing technology is a kind of technology of genome specific region being carried out to selective enrichment, it separates interested region by suitable method from genome, and then checked order in this target area, this at low cost targetedly genomics research have very important meaning.
(Brenet F et al.2011 is found in research recently; Harder A et al.2010; Suzuki MM et al.2008): exon methylate recognize before density ratio and expect more general, and first exon and First Intron methylate and distribute greatlyr with the exon in downstream and the intron distributional difference that methylates, the methylated degree of downstream area is not closely related with genetic expression mostly.In a word, transcription initiation site around methylates, first exon region methylate with the relation of gene silencing than upstream promoter region methylate must be related to it more close.Analysis exon methylates the research of genetic expression is of great importance.Existing research is after genome bisulfite is processed mostly, according to C → T conversion, designs the probe of specific target area, the methylated fragment in acquisition target region.This method was carried out bisulfite processing by genome before sequence capturing, had increased the difficulty of probe design and had reduced the efficiency of sequence capturing, the coverage in object region, thereby having limited the popularity of applying.
Summary of the invention
For the problem existing in the above-mentioned research method that methylates based on bisulfite, the present invention proposes a kind of binding sequence and catch the method checking order with bisulfite, the method can be caught the target area in several samples by a sequence capturing test, evaluating objects region methylation state accurately and effectively, reduce the difficulty of probe design, increase the feasibility of operation and application, making that the high-precision detection that methylates of high-throughput is carried out in interested target sequence and region in full genome becomes a reality, and the method is pointed strong and the cost-saving and feature time also.
The object region the present invention relates to is exon region, exon group is caught with bisulfite order-checking and combined, the distribution situation that methylates that detects the full genome exon of mankind group with this, this finds that to research exon is with a wide range of applications to the effect of gene expression regulation.The main flow process of this technology: first genomic dna is interrupted at random, add after given joint and carry out sequence capturing by solution hybridization method, catch in the DNA that gets off and add foreign DNA and then carry out bisulfite processing, the DNA handling well is carried out to pcr amplification, TA clone detects bisulfite treatment effect, then library is carried out detection by quantitative and carried out high-flux sequence with new-generation sequencing instrument, finally carry out high-precision methylation status analysis within the scope of specific regions.This technical scheme is mainly made up of following five parts: selection, library preparation, sequence capturing, bisulfite processing, the order-checking of upper machine and the data analysis of sequence capturing probe.
Figure BDA0000062872370000031
the selection of sequence capturing probe
The technical program middle probe simplicity of design, can be based on current liquid phase or solid phase chip hybridization technique designing probe, probe length can be from 60mer-120mer not etc.As select probe SureSelect Human All Exon 38M kit (Agilent), this probe to cover the full exon region of people and part miRNA region, with a wherein chain complementation in genome object region, mean length is 120mer.
Figure BDA0000062872370000032
library preparation
The fragmentation of step 1 sample gene group DNA and foreign gene group DNA
Initial object research material and can be any species (for example people as the material of foreign gene group DNA, plant, insect) genomic dna, utilize physics or chemical process genomic dna to be interrupted to the fragment into 200-300bp as ultrasonic wave fragmentation method.Foreign gene group DNA preferably selects the λ DNA that does not methylate and modify; its effect be in the time that bisulfite is processed together with sample efficient co-processing; micro-DNA fragmentation is played a protective role, reduce to greatest extent the destruction of bisulfite to minim DNA.
Get and there is no the complete genome group DNA approximately 5 μ g of the pollution such as albumen, RNA and the λ DNA as foreign gene group, by ultrasonic fragmentation method, genomic dna is interrupted to the fragment into big or small 200-300bp.
Step 2 genomic dna end modified
DNA after interrupting at random reclaims after purifying, end reparation is carried out in effect by enzymes such as T4DNA polysaccharase (T4DNAPolymerase), Klenow fragment (Klenow Fragment) and T4 polynucleotide kinases (T4Polynucleotide Kinase), forms the DNA fragmentation of flat end.Then utilize Klenow fragment (Klenow Frgment) (3 '-5 ' exo-) polysaccharase and dATP to add " A " base at 3 ' end of the sequence filling.
Step 3 PEI (Paired-end Index) joint that methylates (Adapter), also referred to as the methylate connection of joint of two end labels
The sequence that 3 ' end adds " A " base under the effect of T4DNA ligase enzyme (T4DNA Ligase) with particular design and the joint modified of methylating (C site methylate modification) be connected, purifying reclaims, as reclaimed after the DNA in purification reaction system with MiniElute PCR Purification Kit (Qiagen), it is carried out quantitatively, as adopted Qubit (Invitrogen) method etc.
Figure BDA0000062872370000041
sequence capturing
Get 500-1000ng connection product library and carry out liquid phase or solid-phase hybridization, as adopt Agilent solution hybridization platform or Nimblegen solid phase or solution hybridization platform to carry out, after adding tab closure sequence to wait to hybridize in hybridization system simultaneously, collect sequence the purifying of catching by methods such as sex change, obtain the DNA molecular with hybridization probe regional complementarity.
Figure BDA0000062872370000042
bisulfite processing
By catch the DNA getting off add 200ng fragmentation λ DNA, then one reinstate bisulfite processing, as utilize EZ DNA Methylation-Gold Kit tM(ZYMO) make non-methylated cytosine be converted to uridylic.
Glue purification is cut in step 4 pcr amplification and library
DNA after changing with bisulfite is template, add PCR primer sequence, use for the warm start taq enzyme of the DNA after bisulfite conversion and carry out pcr amplification (r-taq of available routine or other polymeric enzymatic amplification), amplified production can carry out purifying by for example following three kinds of methods, is respectively magnetic beads for purifying, purification column purifying, 2% agarose gel electrophoresis purifying.It is quantitative that the product that purifying reclaims carries out QPCR, then treats machine order-checking.
Figure BDA0000062872370000051
upper machine order-checking and data analysis
By the sequence of catching after bisulfite is processed at s-generation order-checking platform, for example on Solexa order-checking platform, adopt the method checking order while synthesize to carry out sequencing.In order to make a distinction after order-checking deriving from DNA library prepared by different samples, the sequence label of 6bp or 8bp is incorporated into a side of fragment by joint or PCR primer, upper machine order-checking after can conveniently different libraries directly being mixed like this.When analytical data, the sequence of reference is hg18 (known whole genome sequence)
The analysis process of the raw data of order-checking gained is with reference to LI Y et al., Nature (2008) (JWang, et al., (2008) .The DNA Methylome of Human Peripheral Blood Mononuclear Cells.Nature, 456:60.), in actual applications, the selection meeting of reference data is along with studied genome source is different and different, fundamental analysis process comprises following key step: the C of the data normal chain that order-checking is obtained all changes into T, the G of complementary strand all changes into A, use SOAP program (Li R, Li Y., Kristiansen K. & Wang, J. (2008) .SOAP:short oligonucleotide alignment program.Bioinformatics, 24:713-714.) compare and be all converted into T to C respectively, G is all converted into the hg18 of A with reference on genome, wherein allow the mispairing of two bases.Then statistics comparison is in the section of reading (reads) of unique position, target area, and compare the above section of reading of reference sequences with these, filter as filtration parameter using Depth > 4, using the information that methylates in site by after filtering as the actual information that methylates measuring, and take carry out follow-up comparison information analysis for basic.
The sequence measurement that carries out again sequence capturing from the bisulfite processing postgenome of conventional study report is different, in the present invention, first carries out the sequence capturing of specific region, uses afterwards sulfiting again.Therefore, in the process of sequence capturing, the methylated state of sample does not change, and just does not need to consider the impact of object region methylation level on sequence capturing efficiency when design capture probe.Can catch all target areas and verify its situation that methylates, distribute to come design chips and needn't forethought methylate again, the data of acquisition are truer and fully.The present invention can make the probe design method in sequence capturing comparatively simple, and capture rate is high, so can be just analyzing gene group methylation level comprehensively and exactly.
The technical program middle probe simplicity of design, can be based on current liquid phase or solid phase chip hybridization technique designing probe, probe length can be from 60-120mer not etc.
This technology is in the time of library construction, the connection of label (Index) joint methylates before sequence capturing, can, by carry out sequence capturing after several sample mix simultaneously, greatly provide cost savings like this, thus realize large sample amount high-throughput methylation level detect.
Embodiment
Below in conjunction with embodiment, embodiment of the present invention are described in detail, but it will be understood to those of skill in the art that the following example is only for the present invention is described, and should not be considered as limiting scope of the present invention.
One aspect of the present invention provides the one high-flux detection method that methylates, it comprises the following steps: the selection of sequence capturing probe, library preparation, sequence capturing, bisulfite processing, the order-checking of upper machine and data analysis, wherein sequence capturing and bisulfite are processed the joint Connection Step prepared in library and pcr amplification and library and cut between glue purification step and carry out.
In a specific embodiment of the present invention, the described probe methylating in high-flux detection method can be the probe for liquid phase or solid phase chip hybridization, preferably the length of described probe is 60mer-120mer, and more preferably the mean length of described probe is 120mer.
In a specific embodiment of the present invention, the described library preparation process methylating in high-flux detection method comprises the connection of the label joint that methylates.
In a specific embodiment of the present invention, in the high-flux detection method that methylates, before sequence capturing, carry out following library preparation process:
The fragmentation of step 1 sample gene group DNA and foreign gene group DNA
Sample gene group DNA and can be any species as foreign gene group DNA, includes but not limited to the genomic dna of people, plant, insect; Utilize physics or chemical process, preferred ultrasonic wave fragmentation method interrupts preferred 5 μ g genomic dnas at random, preferably interrupts the fragment into 200-300bp; Preferably foreign gene group DNA is the λ DNA that does not methylate and modify;
Step 2 genomic dna end modified
DNA after interrupting at random reclaims after purifying, carries out end reparation by the effect of preferred T4DNA polysaccharase, Klenow fragment and T4 polynucleotide kinase enzyme, forms the DNA fragmentation of flat end; Then preferably utilize Klenow fragment (3 '-5 ' exo-) polysaccharase and dATP to add " A " base at 3 ' end of the sequence filling;
The methylate connection of joint of step 3 PEI
3 ' end adds that the sequence of " A " base is at ligase enzyme, under the effect of preferred T4DNA ligase enzyme with methylate modifications, the preferred C site joint of modifying that methylates connects, purifying recovery, carries out quantitatively it.
In a specific embodiment of the present invention, the PEI described in high-flux detection method that the methylates joint that methylates comprises:
The PE Index_ adapter 1 that methylates:
Phos/ tCAAGTaGATCGGAAGAGCACACGTCTGAACTCCAGTCAC (SEQ ID NO:1, all C site methylate modification), and
The PE Index_ adapter 2 that methylates:
TACACTCTTTCCCTACACGACGCTCTTCCGATCTA CTTGAT(SEQ ID NO:2)。
Can change in the present invention the base of underscore part in above-mentioned joint sequence, thereby produce more, different Tag primers, for multiple different sample.
In a specific embodiment of the present invention, the sequence capturing step methylating in high-flux detection method comprises the connection product library in the preparation process of library at solid phase or solution hybridization platform, preferably solution hybridization platform is hybridized, wherein preferably use 500-1000ng to connect the hybridization of product library, in hybridization system, add tab closure sequence simultaneously; After waiting to hybridize, collect sequence the purifying of catching by methods such as sex change, obtain and the DNA molecular of hybridization probe regional complementarity, wherein preferably externally aobvious subgroup is caught.
In a specific embodiment of the present invention, the tab closure sequence using in the sequence capturing step methylating in high-flux detection method is the PEI PEI that the uses complementary sequence of joint that methylates that methylates in the Connection Step of joint, and is selected from
The PE Index_ adapter 1 that methylates:
Phos/TCAAGTAGATCGGAAGAGCACACGTCTGAACTCCAGTCAC (SEQ ID NO:1), or
The PE Index_ adapter 2 that methylates:
The complementary sequence of TACACTCTTTCCCTACACGACGCTCTTCCGATCTACTTGAT (SEQ ID NO:2), all C site of wherein said joint is all through the modification that methylates.
In a specific embodiment of the present invention, the bisulfite treatment step methylating in high-flux detection method comprises outside the pale of civilization the DNA fragmentation obtaining in sequence capturing step source genomic dna, one reinstates bisulfite processes, and foreign gene group DNA preferably 200ng through the λ of fragmentation DNA.
The present invention further provides the one high-flux detection method that methylates, wherein after bisulfite is processed, carried out the step that glue purification is cut in pcr amplification in the preparation of following library and library:
DNA after changing with bisulfite is template, add PCR primer sequence to carry out pcr amplification, the polysaccharase that wherein used comprises for the warm start taq enzyme of the DNA after bisulfite conversion, conventional r-taq or other polymeric enzymatic amplification, the warm start taq enzyme of the DNA after preferably changing for bisulfite;
Amplified production is carried out to purifying, and purification process includes but not limited to magnetic beads for purifying, purification column purifying, 2% agarose gel electrophoresis purifying; It is quantitative that the product that purifying reclaims carries out QPCR, then treats machine order-checking.
In a specific embodiment of the present invention, pcr amplification step further comprises a side that the sequence label of preferred 6bp or 8bp is incorporated into DNA fragmentation by joint or PCR primer.
In a specific embodiment of the present invention, methylating, the PCR primer sequence using in root pcr amplification in high-flux detection method comprises:
The public primer of P1:
AATGATACGGCGACCACCGAGATCTACACTCTTTCCCTACACGACGCTCTTCCGATCT(SEQ ID NO:3),
Index1 primer:
CAAGCAGAAGACGGCATACGAGAT cTTGATgTGACTGGAGTTCAGACGTGTGCTCTTCCGATCT (SEQ ID NO:4), and
Index 2 primers:
CAAGCAGAAGAC GGCATAC GAGAT TCAAGTGTGACTGGAGTTCAGACGTGTGCTCTTCCGATCT(SEQ ID NO:5)。
Wherein can change the base of underscore part in above-mentioned PCR primer, thereby produce more, different Tag primers, for multiple different sample.
In a specific embodiment of the present invention, the upper machine order-checking and the data analysis that methylate in high-flux detection method comprise that the sequence that bisulfite treatment step is obtained is at order-checking platform, preferably s-generation order-checking platform, on preferred Solexa order-checking platform, check order, the line number of going forward side by side compares according to one's analysis.
The sequencing library that the present invention provides the high-flux detection method that methylates according to the present invention to build on the other hand.
The present invention provides the purposes for the high-flux sequence that methylates according to the sequencing library of the high-flux detection method structure that methylates according to the present invention on the other hand.
Further, the high-flux detection method that methylates according to the present invention is applicable to exon group order-checking, further preferably for detection of the distribution situation that methylates of the full genome exon of mankind group.
Accompanying drawing explanation
Fig. 1: MDC PCR 2100 detected results.
Fig. 2: YH PCR 2100 detected results.
Fig. 3: YH and the MDC coverage on every karyomit(e) distributes.Wherein red post YH represents YH Whole Blood Genomic DNA, green post mDC (Mature dendritic cells) represents ripe dendritic cell DNA, X-coordinate (chromosome) represents karyomit(e) number, and ordinate zou (coverage of ech chromosome) represents order-checking coverage.
Fig. 4: YH and MDC be depth profile on every karyomit(e).Its Green post YH represents YH Whole Blood Genomic DNA, and red post mDC represents ripe dendritic cell DNA, and X-coordinate (chromosome) represents karyomit(e) number, and ordinate zou (depth of ech chromosome) represents the order-checking degree of depth.
Fig. 5: YH MDC catches the distribution of fragment on gene.
Fig. 6: MDC catches target area and full genome bisulfite data comparison.X-coordinate and ordinate zou represent respectively methyl rate (methylation rate in target region) that target acquistion draws and the methyl rate of genome sequencing.Pearson Correlation represents Pearson correlation coefficient.
Fig. 7: YH catches target area and full genome bisulfite data comparison.X-coordinate and ordinate zou represent respectively methyl rate (methylation rate in target region) that target acquistion draws and the methyl rate of genome sequencing.Pearson Correlation represents Pearson correlation coefficient.
Embodiment
Main laboratory apparatus list
Figure BDA0000062872370000101
Figure BDA0000062872370000111
Related reagent list
Sequence list (5 '-> 3 ', what underscore part represented is label)
Figure BDA0000062872370000121
DNA segment:
Use covaris-S2 to interrupt instrument, YH Whole Blood Genomic DNA (Yan Di and Huang Di, two legendary rulers of remote antiquity's Whole Blood Genomic DNA, the genomic dna extracting from only one China man blood) and MDC DNA (people's immunocyte is DNA) sample respectively interrupt (external source λ DNA is also like this) take 5 μ g as initial amount.According to following parameter setting:
Treatment1 (processing 1) Duty/cycle (%) (duty ratio) 10
Intensity (intensity) 5
Cycle/burst (circulation/pulse) 200
Time (min) (time (second)) 50
Treatment2 (processing 2) Time (s) (time (second)) 0
Treatment (processing 3) Time (s) (time (second)) 0
Treatment4 (processing 4) Time (s) (time (second)) 0
Cycles (circulation) 3
After 2% agarose gel electrophoresis detection is qualified, (DNA fragmentation master tape concentrates between 200-300bp sample after interrupting, pollute without albumen, RNA), through QIAquick PCRPurification Kit (Qiagen) purifying, back dissolving is in the Elution Buffer (EB) of 32 μ l.
Genomic dna end modified:
DNA according to the form below obtained in the previous step is prepared to end in the centrifuge tube of 1.5ml and repairs reaction system:
Figure BDA0000062872370000131
Then centrifuge tube is put into the upper reaction of Thermomixer (Eppendrf) 30min that is adjusted to 20 ℃.React rear and carried out purifying with QIAquick PCR Purification Kit (Qiagen), finally sample has been dissolved in 34 μ l Elution Buffer (EB).
DNA fragmentation 3 ' end adds " A " base:
DNA according to the form below obtained in the previous step is prepared and is added " A " reaction system in the centrifuge tube of 1.5ml:
Figure BDA0000062872370000132
Then centrifuge tube is put into the upper reaction of Thermomixer (Eppendrf) 30min that is adjusted to 37 ℃.React rear and carried out purifying with MiniElute PCR Purification Kit (Qiagen), finally sample has been dissolved in to 16 μ L Elution Buffer (EB).
The methylate connection of joint of PEI:
DNA according to the form below obtained in the previous step is prepared to the joint ligation system that methylates:
Figure BDA0000062872370000133
Then be put into the upper reaction of Thermomixer (Eppendrf) 15min that is adjusted to 20 ℃.React rear and carried out purifying with MiniElute PCR Purification Kit (Qiagen), sample has been dissolved in to 12 μ L Elution Buffer (EB).Then adopt Qubit (Invitrogen) to carry out concentration detection, according to concentration detected result, the concentration of DNA is controlled to 147ng/ul (by concentrated or dilution mode).
Sequence capturing:
Take a hybridization (reagent is all from SureSelect Human All Exon 38Mkit) as example:
A. the preparation of hybridizing reagent A:
The preparation (C) of b.SureSelect Oligo Library Mix:
In PCR pipe, after packing 5 μ L Oligo Capture Library, add RNase Block (the wherein RNase Block: nuclease-free water=1: 3), be placed on ice after 2 μ L dilutions.
C. the configuration of library and block mixed liquid B:
Library 3.4μL 147ng/μL
SureSelect Block#1 2.5μL
SureSelect Block#2 2.5μL
Index Block 0.6μL
Amount to 9μL
The PEI joint that methylates is double-stranded DNA, while wherein using PE-index to methylate Adapter1 strand (SEQ ID NO:1), seal accordingly sequence (Index Block) and be the methylate complementary sequence of Adapter1 strand (SEQID NO:1) of PE-index in above-mentioned PEI methylates the Connection Step of joint; While using PE-index to methylate Adapter2 strand (SEQ ID NO:2) in above-mentioned PEI methylates the Connection Step of joint, the IndexBlock adding is the methylate complementary sequence of Adapter2 strand (SEQ ID NO:2) of PE-index.
In PCR pipe, configure mixed liquid B, and mixed solution mixed with pipettor, cover tightly pipe lid, be placed in PCR instrument by following condition operation thermal cycler program:
Figure BDA0000062872370000142
D. hybridizing reagent A is placed in thermal cycler and is at least kept 5 minutes in 65 ℃.
E. SureSelect Oligo Library Mix (C) is placed in thermal cycler and is at least kept 2 minutes in 65 ℃.
F. keep PCR pipe in 65 ℃, with 20 μ l pipettors, the Hybridization BufferA of 13 μ l is transferred in SureSelect Oligo Library Mix (C) pipe rapidly; Then with 20ul pipettor, the mixed solution in library and block mixed liquid B is all transferred in SureSelect Oligo Library Mix (C) pipe rapidly, with pipettor slowly up and down inhale blow 8 to 10 times.
G. cover tightly PCR pipe lid, in 65 ℃ of (thermal cycler heat lid is made as 105 ℃) hybridization 24 hours.
Sample wash-out (reagent is all from SureSelect Human All Exon 38M kit):
A. magnetic bead is prepared
Each hybridization is got 50 μ l DynabeadsM-280Streptavidin magnetic beads in a new 1.5mL centrifuge tube;
Use 200 μ L SureSelect Binding buffer to wash three times;
Finally add the resuspended magnetic bead of 200 μ l SureSelect Binding buffer.
B. elution samples
By be equipped with hybrid mixed liquid and magnetic bead centrifuge tube symmetry be fixed on BD Clay AdamsNutator Mixer or similarly upper 360 degree of device rotate and mix, under room temperature, hatch 30 minutes, instantaneous centrifugal 3s is placed on and on magnetic frame, removes supernatant;
With the resuspended magnetic bead of 500 μ l SureSelect Wash Buffer#1, and on whirlpool mixed instrument, vibrate and within 5 seconds, mix sample, under room temperature, hatch sample 15 minutes;
Wash according to the following steps magnetic bead and obtain object fragment:
A) centrifuge tube is shifted and is placed on magnetic frame, leave standstill 5-10min to clarification, and remove supernatant liquor as far as possible;
B) with the resuspended magnetic bead of 500 μ l SureSelect Wash Buffer#2, and on whirlpool mixed instrument, vibrate 5 seconds to mix sample;
C) sample is put in Thermomixer to 65 ℃ and hatched 10 minutes, put upside down centrifuge tube to mix sample, instantaneous centrifugal 3 seconds;
D) repeat " step a is to step c " 2 times, Wash Buffer#2 is removed clean as far as possible.
E) add 50 μ L SureSelect Elution Buffer, and on whirlpool mixed instrument, vibrate 5 seconds with resuspended magnetic bead, sample is put under room temperature and is hatched 10 minutes;
F) centrifuge tube is shifted and is placed on Dynal magnetic frame, leave standstill 5-10min to clarification;
G) with pipettor, the elutriant that contains sample is transferred in a new 1.5ml centrifuge tube;
H) in the DNA catching, add 50 μ l SureSelect Neutralization Buffer.
I) use 1.8 times of Ampure Beads to remove the salinity of catching in sample solution.
Catching product adds foreign DNA to carry out bisulfite co-processing:
Catch the external source λ DNA that adds 200ng fragmentation in product, then add bisulfite to process 2h.Bisulfite is processed and is adopted EZ DNA Methylation-Gold Kit tM(ZYMO), concrete steps are as follows:
A) preparation of CT Conversion Reagent: take out CTConversion Reagent (solid mixture) from test kit test kit, then add water, the M-Dissolving Buffer of 50 μ l and the M-Dilution Buffer of 300 μ l of 900 μ l in the CTConversion Reagent of a pipe.At room temperature dissolve and shake 10 minutes or on shaking table, shake 10 minutes.
B) preparation of M-WASH BUFFER: add the ethanol of 24ml 100% and prepare final operable M-Wash Buffer in M-WashBuffer, make hook (using marking pen mark) after adding ethanol on lid, show to use.
C) DNA to be converted is installed in PCR pipe to moisturizing to 20 μ l according to dividing.
D) in PCR pipe, add the CT Conversion Reagent of 130 μ l, by flicking test tube or pipettor operates biased sample.
E) sample hose is put into operation according to the following steps on PCR instrument:
Place 5 minutes for 98 ℃
Place 2.5 hours for 64 ℃
Carry out at once next step operation or storage (maximum 20 hours) at 4 ℃.
F) the M-Binding Buffer that adds 600 μ l is to Zymo-Spin IC tMin Column, and post is put in the Collection Tube providing as test kit.
G) filling sample is to Zymo-Spin IC tMcolumn contains M-Binding Buffer.Close the lid post is put upside down and carried out biased sample for several times.
H) (> 10,000xg) centrifugal 30 seconds at full speed, removes effluent liquid.
I) add the M-Wash Buffer of 200 μ l in post, centrifugal 30 seconds at full speed.
J) the M-Desulphonation Buffer that adds 200 μ l places 15 minutes in post and under room temperature (20 ℃-30 ℃), after cultivation, and centrifugal 30 seconds at full speed.
K) add the M-Wash Buffer of 200 μ l in post, centrifugal 30 seconds at full speed; Add again the M-Wash Buffer of 200 μ l and centrifugal 30 seconds.
L) the M-Elution Buffer that directly adds 10 μ l is in base for post matter.Post is placed in the pipe of 1.5ml, at full speed centrifugal come eluted dna.
Pcr amplification and library size are selected: DNA obtained in the previous step is prepared to PCR reaction system by following system:
The DNA of bisulfite processing 10μl
dNTP(2.5mM) 2μl
10x pcr buffer 2.5ul
JumpStartTM Taq DNA Polymerase 0.3μl
The public primer of P1 0.5μl
Index1 or 2 primers 0.5μl
dH 2O 9.2μl
Cumulative volume 25μl
PCR reaction conditions:
Figure BDA0000062872370000171
Amplified production is carried out to purifying with PCR Purification Kit (Qiagen), then carry out electrophoresis with 2% agarose gel, then glue being cut in object size library selects, adopt MiniElute PCR Purification Kit (Qiagen) to carry out glue purification recovery, finally library is dissolved in 20 μ l EB.
Library detection: TA clone detects sulphite to genomic efficiency of conversion; Agilent 2100Bioanalyzer detects library output; QPCR detection by quantitative library output.
Order-checking and data analysis:
Sample is carried out to two end sequencings in Solexa order-checking platform.By with the result of the full genome sulphite order-checking in same sample source and the data analysis comparison that full exon sequence is caught, analytical sequence is caught the feasibility of methylation analysis method.
Experimental result:
The Agilent 2100Bioanalyzer detected result of 1PCR product
The PCR product Agilent 2100Bioanalyzer detected result of Fig. 1 and Fig. 2 shows that the initial YH of 5 μ g (Yan Di and Huang Di, two legendary rulers of remote antiquity's Whole Blood Genomic DNA) and MDC DNA (ripe dendritic cell DNA) can build and utilize high-throughput new-generation sequencing instrument to carry out the library that methylates of high-flux sequence, show that in conjunction with actual sequencing data analytical results hereinafter described this inventive method is practical, can be applicable in practical study.
2TA clones result
Table 1YH and MDC DNA build library text storehouse quality examination result
Figure BDA0000062872370000181
2 libraries are selected respectively 36 and 41 clones and are carried out quality examination, and results conversion efficiency all, more than 99%, illustrates that bisulfite processing has realized efficient conversion.
3 data analyses
3.1 with full genome alignment rate comparison
Table 2YH, MDC and routine are built storehouse sequencing data result comparison result
As can be seen from the above table, the comparison rate of YH and MDC is respectively 91.80% and 87.17%.Comparison rate, the transformation efficiency of capture rate and BS is all in normal range.
The distribution plan of 3.2 coverages
Table 3YH, MDC coverage distribution plan
Figure BDA0000062872370000191
The distribution situation of table 3 display data under the different multiplier covering, in very concentrated being distributed in target area of current order-checking amount data.
3.2.1 the average coverage of each karyomit(e) distributes and checks order depth profile as can be seen from Figure 3, and sequencing data can cover the target area on every karyomit(e) substantially.
Fig. 4 is presented in target area, and the degree of depth on every karyomit(e) is all in 50 × left and right, and the comparatively accurate site information that methylates can be provided in the process of calculating methyl rate.
What Fig. 5 showed is the distributions on karyomit(e) of the data that capture on No. 12 karyomit(e)s, can find out that data only have enrichment in the place that has target area.
3.3 data dependence comparisons
Fig. 6 and Fig. 7 be respectively two samples (MDC and YH) use the sequencing data that obtains of the inventive method methylate correlation analysis with use full genome normally to build storehouse to obtain the methylate comparison of correlation analysis of corresponding sequencing data, data are filter depth 10 × the methyl rate distribution situation of above data; Can see that the methyl rate on same loci is basically identical, the relation conefficient of MDC sample is 0.94 (Pearson correlation coefficient), and the relation conefficient of YH sample is 0.93.
By being carried out to specific region, YH and MDC catch library construction, obtain sequencing data and conventional library order-checking corresponding data and carry out information analysis comparative result, from comparison efficiency, coverage, every chromosomal methyl rate and dependency, each side fraction of coverage is all good, methyl rate also has good consistence, and technique can reach very dark data volume among a small circle, these have all illustrated that the methylate high-flux sequence research of employing specific region is practicable, this method reduces the difficulty of probe design, increase the feasibility of operation and application, making that the high-precision detection that methylates of high-throughput is carried out in interested target sequence and region in full genome becomes a reality, and pointed strong and the cost-saving and feature time.
Reference:
[1]Brena RM,Huang TH,Plass C.Quantitative assessment of DNA methylation:Potential applications for disease diagnosis,classification,and prognosis in clinical settings.J Mol Med.2006May;84(5):365-77.
[2]Egger G,Liang G,Aparicio A et al.Epigenetics in human disease and prospects for epigenetic therapy.Nature 2004 May27;429(6990):457-63.
[3]Gu H,Bock C,Mikkelsen TS,et al.Genome-scale DNAmethylation mapping of clinical samples at single-nucleotide resolution.Nat Methods.2010 Feb;7(2):133-6.
[4]Hodges E,Smith AD,et al.High definition profiling of mammalian DNA methylation by array capture and single molecule bisulfite sequencing.Genome Res.2009Sep;19(9):1593-605.
[5]Brenet F,Moh M,et al.DNA methylation of the first exon is tightly linked to transcriptional silencing.PLoS One.2011 Jan18;6(1):e14524.
[6]Harder A,Titze S,et al.Monozygotic twins with neurofibromatosis type 1(NF1)display differences in methylation of NF1 gene promoter elements,5′untranslated region,exon and intron1.Twin Res Hum Genet.2010Dec;13(6):582-94.
[7]Suzuki MM,Bird A.DNA methylation landscapes:provocative insights from epigenomics.Nat Rev Genet.2008Jun;9(6):465-76.
[8]Li R,Li Y.,Kristiansen K.& Wang,J.(2008).SOAP:short oligonucleotide alignment program.Bioinformatics,24:713-714.
Figure IDA0000062872440000011
Figure IDA0000062872440000021

Claims (34)

1.一种甲基化高通量检测方法,其包括以下步骤:序列捕获探针的选择、文库制备、序列捕获、重亚硫酸盐处理、上机测序及数据分析,其特征在于序列捕获和重亚硫酸盐处理在文库制备的接头连接步骤和PCR扩增及文库切胶纯化步骤之间进行。1. A methylation high-throughput detection method, comprising the following steps: selection of sequence capture probes, library preparation, sequence capture, bisulfite treatment, computer sequencing and data analysis, characterized in that sequence capture and Bisulfite treatment is performed between the adapter ligation step of library preparation and the PCR amplification and library gel purification steps. 2.根据权利要求1的方法,其中所述探针是用于液相或者固相芯片杂交的探针。2. The method according to claim 1, wherein the probe is a probe for liquid-phase or solid-phase chip hybridization. 3.根据权利要求2的方法,其中所述探针的长度是60mer-120mer。3. The method according to claim 2, wherein the length of the probe is 60mer-120mer. 4.根据权利要求3的方法,其中所述探针的平均长度为120mer。4. The method according to claim 3, wherein the average length of the probes is 120mer. 5.根据权利要求1的方法,其中所述文库制备步骤包括甲基化标签接头的连接。5. The method according to claim 1, wherein said library preparation step comprises ligation of methylated tag adapters. 6.根据权利要求1的方法,其中所述文库制备步骤包括PEI甲基化接头的连接。6. The method according to claim 1, wherein said library preparation step comprises ligation of PEI methylated adapters. 7.根据权利要求1的方法,其中在序列捕获前进行如下文库制备步骤:7. The method according to claim 1, wherein the following library preparation steps are carried out before sequence capture: 步骤一样品基因组DNA以及外源基因组DNA的片段化Step 1 Fragmentation of sample genomic DNA and exogenous genomic DNA 样品基因组DNA和作为外源基因组DNA是任意物种,包括人、植物、昆虫的基因组DNA;利用物理或化学方法进行随机打断;The sample genomic DNA and the exogenous genomic DNA are any species, including human, plant, and insect genomic DNA; use physical or chemical methods for random interruption; 步骤二基因组DNA的末端修饰Step 2: Terminal modification of genomic DNA 随机打断后的DNA回收纯化后,进行末端修复,形成平末端的DNA片段;然后在补平的序列的3’末端加上“A”碱基;After the randomly interrupted DNA is recovered and purified, the end is repaired to form a blunt-ended DNA fragment; then an "A" base is added to the 3' end of the blunted sequence; 步骤三PEI甲基化接头的连接Step 3: Connection of PEI methylated linker 3’末端加上“A”碱基的序列在连接酶的作用下与甲基化修饰的接头进行连接,纯化回收,对其进行定量。The sequence with "A" base added at the 3' end is ligated with the methylated linker under the action of ligase, purified and recovered, and quantified. 8.根据权利要求7的方法,其中步骤一中,利用超声波片段化方法将5μg基因组DNA进行随机打断。8. The method according to claim 7, wherein in step 1, 5 μg of genomic DNA is randomly fragmented by ultrasonic fragmentation. 9.根据权利要求7的方法,其中步骤一中,将基因组DNA打断为200-300bp的片段。9. The method according to claim 7, wherein in step one, the genomic DNA is broken into fragments of 200-300bp. 10.根据权利要求7的方法,其中步骤一中,外源基因组DNA是没有甲基化修饰的λDNA。10. The method according to claim 7, wherein in step one, the exogenous genomic DNA is lambda DNA without methylation modification. 11.根据权利要求7的方法,其中步骤二中,通过T4DNA聚合酶、Klenow片段和T4聚核苷酸激酶酶的作用进行末端修复。11. The method according to claim 7, wherein in step 2, end repair is carried out by the action of T4 DNA polymerase, Klenow fragment and T4 polynucleotide kinase enzyme. 12.根据权利要求7的方法,其中步骤二中,利用Klenow片段(3’-5’exo-)聚合酶及dATP在补平的序列的3’末端加上“A”碱基。12. The method according to claim 7, wherein in step 2, an "A" base is added to the 3' end of the filled-in sequence using Klenow fragment (3'-5'exo-) polymerase and dATP. 13.根据权利要求7的方法,其中步骤三中,3’末端加上“A”碱基的序列在T4DNA连接酶的作用下与甲基化修饰的接头进行连接。13. The method according to claim 7, wherein in step 3, the sequence with the "A" base added to the 3' end is connected to the methylated linker under the action of T4 DNA ligase. 14.根据权利要求7的方法,其中步骤三中,3’末端加上“A”碱基的序列与C位点甲基化修饰的接头进行连接。14. The method according to claim 7, wherein in step 3, the sequence with the "A" base added at the 3' end is connected to the linker modified by methylation at the C site. 15.根据权利要求7的方法,其中所述PEI甲基化接头包括:15. The method according to claim 7, wherein said PEI methylated linker comprises: PE Index_甲基化adapter1:PE Index_methylation adapter1: Phos/TCAAGTAGATCGGAAGAGCACACGTCTGAACTCCAGTCAC,和Phos/ TCAAGT AGATCGGAAGAGCACACGTCTGAACTCCAGTCAC, and PE Index_甲基化adapter2:PE Index_methylation adapter2: TACACTCTTTCCCTACACGACGCTCTTCCGATCTACTTGAT,其中所述接头的所有C位点均经甲基化修饰。TACACTCTTTCCCTACACGACGCTCTTCCGATCTA CTTGAT , wherein all C sites of the linker are modified by methylation. 16.根据权利要求1的方法,其中序列捕获步骤包括将文库制备步骤中的连接产物文库在固相或液相杂交平台进行杂交,在杂交体系中同时加入接头封闭序列;待杂交完毕后,通过变性等方法收集捕获的序列并纯化,得到与杂交探针区域互补的DNA分子。16. The method according to claim 1, wherein the sequence capture step comprises hybridizing the ligation product library in the library preparation step on a solid-phase or liquid-phase hybridization platform, and adding an adapter blocking sequence to the hybridization system simultaneously; after the hybridization is completed, by Denaturation and other methods are used to collect the captured sequence and purify it to obtain a DNA molecule complementary to the region of the hybridization probe. 17.根据权利要求16的方法,其中序列捕获步骤包括将文库制备步骤中的连接产物文库在液相杂交平台进行杂交。17. The method according to claim 16, wherein the sequence capture step comprises hybridizing the ligation product library in the library preparation step on a liquid phase hybridization platform. 18.根据权利要求16的方法,其中使用500-1000ng连接产物文库杂交。18. The method according to claim 16, wherein 500-1000 ng of the ligation product library is used for hybridization. 19.根据权利要求16的方法,其中是对外显子组进行捕获。19. The method according to claim 16, wherein exome capture is performed. 20.根据权利要求16的方法,其中接头封闭序列是PEI甲基化接头的连接步骤中所使用PEI甲基化接头的互补序列,并且选自20. The method according to claim 16, wherein the adapter blocking sequence is the complementary sequence of the PEI methylated adapter used in the ligation step of the PEI methylated adapter, and is selected from PE Index_甲基化adapter1:PE Index_methylation adapter1: Phos/TCAAGTAGATCGGAAGAGCACACGTCTGAACTCCAGTCAC,或Phos/TCAAGTAGATCGGAAGAGCACACGTCTGAACTCCAGTCAC, or PE Index_甲基化adapter2:PE Index_methylation adapter2: TACACTCTTTCCCTACACGACGCTCTTCCGATCTACTTGAT的互补序列,其中所述接头的所有C位点均经甲基化修饰。The complementary sequence of TACACTCTTTCCCTACACGACGCTCTTCCGATCTACTTGAT, wherein all C sites of the linker are modified by methylation. 21.根据权利要求1的方法,其中重亚硫酸盐处理步骤包括将序列捕获步骤中得到的DNA与片段化外源基因组DNA,一起用重亚硫酸盐进行处理。21. The method according to claim 1, wherein the bisulfite treatment step comprises treating the DNA obtained in the sequence capture step together with the fragmented exogenous genomic DNA with bisulfite. 22.根据权利要求21的方法,其中外源基因组DNA是200ng经片段化的λDNA。22. The method according to claim 21, wherein the exogenous genomic DNA is 200 ng of fragmented lambda DNA. 23.根据权利要求1的方法,其中在重亚硫酸盐处理后进行如下文库制备中的PCR扩增及文库切胶纯化步骤:23. The method according to claim 1, wherein after bisulfite treatment, carry out the PCR amplification and library cutting gel purification steps in the following library preparation: 以重亚硫酸盐转换后的DNA为模板,加入PCR引物序列进行PCR扩增,其中所使用的聚合酶包括针对重亚硫酸盐转换后的DNA的热启动taq酶、常规的r-taq或其它聚合酶扩增;Using bisulfite-converted DNA as a template, add PCR primer sequences for PCR amplification, and the polymerase used includes hot-start taq enzyme for bisulfite-converted DNA, conventional r-taq or other Polymerase amplification; 对扩增产物进行纯化,纯化方法包括磁珠纯化、纯化柱纯化、2%的琼脂糖胶电泳纯化;纯化回收的产物进行QPCR定量,然后待上机测序。The amplified product was purified, and the purification methods included magnetic bead purification, purification column purification, and 2% agarose gel electrophoresis purification; the purified and recovered product was quantified by QPCR, and then sequenced on the machine. 24.根据权利要求23的方法,其中所使用的聚合酶是针对重亚硫酸盐转换后的DNA的热启动taq酶。24. The method according to claim 23, wherein the polymerase used is a hot-start taq enzyme directed against bisulfite-converted DNA. 25.根据权利要求23的方法,其中进一步包括将标签序列通过接头或者PCR引物引入到DNA片段的一侧。25. The method according to claim 23, further comprising introducing a tag sequence into one side of the DNA fragment through a linker or a PCR primer. 26.根据权利要求25的方法,其中标签序列是6bp或8bp的标签序列。26. The method according to claim 25, wherein the tag sequence is a 6bp or 8bp tag sequence. 27.根据权利要求23的方法,其中PCR引物序列包括:27. The method according to claim 23, wherein the PCR primer sequences comprise: P1公用引物:P1 public primers: AATGATACGGCGACCACCGAGATCTACACTCTTTCCCTACACGACGCTCTTCCGATCT,AATGATACGGCGACCACCGAGATCTACACTCTTTCCCTACACGACGCTCTTCCGATCT, Index1引物:Index1 primers: CAAGCAGAAGACGGCATACGAGATCTTGATGTGACTGGAGTTCAGACGTGTGCTCTTCCGATCT,和CAAGCAGAAGACGGCATACGAGAT CTTGAT GTGACTGGAGTTCAGACGTGTGCTCTTCCGATCT, and Index2引物:Index2 primers: CAAGCAGAAGACGGCATACGAGATTCAAGTGTGACTGGAGTTCAGACGTGTGCTCTTCCGATCT。 CAAGCAGAAGACGGCATACGAGATTCAAGTGTGACTGGAGTTCAGACGTGTGCTCTTCCGATCT . 28.根据权利要求1的方法,其中上机测序及数据分析包括将重亚硫酸盐处理步骤得到的序列在测序平台上进行测序,并进行数据分析比较。28. The method according to claim 1, wherein the on-machine sequencing and data analysis comprise sequencing the sequence obtained in the bisulfite treatment step on a sequencing platform, and performing data analysis and comparison. 29.根据权利要求28的方法,其中测序平台是第二代测序平台。29. The method according to claim 28, wherein the sequencing platform is a second generation sequencing platform. 30.根据权利要求28的方法,其中测序平台是Solexa测序平台。30. The method according to claim 28, wherein the sequencing platform is the Solexa sequencing platform. 31.根据权利要求1-30中任一项的方法构建的测序文库。31. A sequencing library constructed according to the method of any one of claims 1-30. 32.根据权利要求31的测序文库用于甲基化高通量测序的用途。32. Use of the sequencing library according to claim 31 for methylation high-throughput sequencing. 33.根据权利要求1-30中任一项的方法,适用于外显子组测序。33. The method according to any one of claims 1-30, adapted for exome sequencing. 34.根据权利要求1-30中任一项的方法,适用于检测人类全基因组外显子组的甲基化分布情况。34. The method according to any one of claims 1-30, which is suitable for detecting the methylation distribution of exomes in the whole human genome.
CN201110133858.1A 2011-05-23 2011-05-23 Methylation high-flux detection method Active CN102796808B (en)

Priority Applications (6)

Application Number Priority Date Filing Date Title
CN201110133858.1A CN102796808B (en) 2011-05-23 2011-05-23 Methylation high-flux detection method
US14/119,266 US9133513B2 (en) 2011-05-23 2012-05-22 High throughput methylation detection method
EP12789707.2A EP2725125B1 (en) 2011-05-23 2012-05-22 High throughput methylation detection method
DK12789707.2T DK2725125T3 (en) 2011-05-23 2012-05-22 HIGH-THROUGHPUT METHYLING DETECTION PROCEDURE
PCT/CN2012/075919 WO2012159564A1 (en) 2011-05-23 2012-05-22 High throughput methylation detection method
HK13104998.9A HK1178211B (en) 2013-04-24 High throughput methylation detection method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201110133858.1A CN102796808B (en) 2011-05-23 2011-05-23 Methylation high-flux detection method

Publications (2)

Publication Number Publication Date
CN102796808A CN102796808A (en) 2012-11-28
CN102796808B true CN102796808B (en) 2014-06-18

Family

ID=47196105

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201110133858.1A Active CN102796808B (en) 2011-05-23 2011-05-23 Methylation high-flux detection method

Country Status (5)

Country Link
US (1) US9133513B2 (en)
EP (1) EP2725125B1 (en)
CN (1) CN102796808B (en)
DK (1) DK2725125T3 (en)
WO (1) WO2012159564A1 (en)

Families Citing this family (34)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9365902B2 (en) * 2013-03-15 2016-06-14 Abbott Molecular Inc. Detection of bisulfite converted nucleotide sequences
CN104450872A (en) * 2013-09-25 2015-03-25 上海市肿瘤研究所 High-throughput multi-sample multi-target sing-base resolution methylation level detection method
CN105400776B (en) * 2014-09-12 2019-12-31 深圳华大智造科技有限公司 Oligonucleotide adapters and their application in the construction of single-strand circular libraries for nucleic acid sequencing
CN104611453B (en) * 2015-02-28 2016-08-24 基因科技(上海)有限公司 The methylated probe groups of diagnostic gene and application thereof
WO2016160965A1 (en) * 2015-03-30 2016-10-06 Rubicon Genomics, Inc. Methods and compositions for repair of dna ends by multiple enzymatic activities
CN104894233B (en) * 2015-04-22 2018-04-27 上海昂朴生物科技有限公司 A kind of multisample multiple clips DNA methylation high-flux sequence method
CN105002567B (en) * 2015-06-30 2017-10-13 北京百迈客生物科技有限公司 Construction method of high-throughput simplified methylation sequencing library without reference genome
CN105332063B (en) * 2015-08-13 2017-04-12 厦门飞朔生物技术有限公司 Construction method of single-tube and high-flux sequencing library
CN105506109A (en) * 2015-12-31 2016-04-20 安诺优达基因科技(北京)有限公司 Unicell simplified representative bisulfite sequencing method and kit
CN105462960B (en) * 2016-01-08 2019-02-12 杭州千基生物科技有限公司 A kind of method of DNA sulfite conversion and purification
CN105925562A (en) * 2016-05-10 2016-09-07 广州嘉检医学检测有限公司 Method and kit for enriching 4000 human pathogenic target genes
CN108504651B (en) * 2017-02-27 2020-09-08 深圳乐土生物科技有限公司 Library construction method and reagent for large sample size mixed library construction of PCR products based on high-throughput sequencing
CN108728903A (en) * 2017-04-21 2018-11-02 深圳市乐土精准医疗科技有限公司 The banking process of thalassemia large sample screening is used for based on high-flux sequence
CN107451419B (en) * 2017-07-14 2020-01-24 浙江大学 Method for generating simplified DNA methylation sequencing data by computer program simulation
CN107475779A (en) * 2017-09-22 2017-12-15 上海美吉医学检验有限公司 Library method for building up and its application suitable for unicellular RRBS sequencings
CN107937985A (en) * 2017-10-25 2018-04-20 人和未来生物科技(长沙)有限公司 A kind of construction method and detection method in micro fragmentation DNA methylation detection library
CN108300766A (en) * 2018-01-16 2018-07-20 四川大学 Methylate to chromatin open zone and mitochondria using transposase the method for research
CN113166809B (en) * 2018-12-29 2023-12-26 深圳华大生命科学研究院 DNA methylation detection method, kit, device and application
CN109609613B (en) * 2019-01-25 2022-03-11 艾吉泰康生物科技(北京)有限公司 DNA hydroxymethylation target region capture sequencing method
JP7203276B2 (en) * 2019-05-21 2023-01-12 深▲せん▼華大智造科技有限公司 Methods and kits for constructing sequencing libraries based on target regions of methylated DNA
CN112553298B (en) * 2019-09-10 2023-06-23 深圳华大生命科学研究院 A library construction method that can improve the bias of genome-wide methylation sequencing
CN110878298B (en) * 2019-11-22 2023-09-15 深圳市易基因科技有限公司 rRNA capture probe and application thereof
CN111020017A (en) * 2019-12-31 2020-04-17 北京优迅医学检验实验室有限公司 Method for evaluating DNA methylation conversion efficiency and application thereof
CN112662760A (en) * 2020-02-25 2021-04-16 博尔诚(北京)科技有限公司 Cancer gene methylation detection system and cancer in-vitro detection method implemented in cancer gene methylation detection system
CN112126987A (en) * 2020-08-19 2020-12-25 深圳思凝一云科技有限公司 Library construction method for sequencing methylation amplicon
CN111945232A (en) * 2020-08-27 2020-11-17 天津诺禾医学检验所有限公司 Y-shaped joint, kit and method for constructing RRBS library
CN112226486B (en) * 2020-09-24 2021-11-30 上海英基生物科技有限公司 Library construction method for targeted enrichment of genome methylated DNA and application thereof
CN112779320B (en) * 2020-12-04 2023-07-14 深圳市易基因科技有限公司 Multi-region DNA methylation detection probe design and detection method thereof
CN114974417A (en) * 2021-06-03 2022-08-30 广州燃石医学检验所有限公司 Methylation sequencing method and device
CN114058681A (en) * 2021-12-01 2022-02-18 大连晶泰生物技术有限公司 A methylation mutation detection method and kit based on target region capture
CN113948150B (en) * 2021-12-21 2022-04-19 北京迈基诺基因科技股份有限公司 JMML related gene methylation level evaluation method, model and construction method
CN114540497B (en) * 2022-02-25 2024-02-27 博尔诚(北京)科技有限公司 Marker for bladder cancer screening, probe composition and application thereof
AU2023245564A1 (en) * 2022-04-01 2024-11-14 Twist Bioscience Corporation Libraries for methylation analysis
CN118186057A (en) * 2022-12-13 2024-06-14 深圳湾实验室 Screening method of free DNA marker, DNA marker and application thereof

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2003035860A1 (en) * 2001-10-24 2003-05-01 The Trustees Of Columbia University In The City Of New York A method for gene identification based on differential dna methylation
WO2003087774A3 (en) * 2002-04-12 2003-12-31 Univ Missouri Ecist microarrays for dual screening of dna hypermethylation and gene silencing
WO2005090607A1 (en) * 2004-03-08 2005-09-29 Rubicon Genomics, Inc. Methods and compositions for generating and amplifying dna libraries for sensitive detection and analysis of dna methylation
CN101802223A (en) * 2007-08-15 2010-08-11 香港大学 Methods and compositions for high throughput bisulfite DNA-sequencing and uses thereof

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008096146A1 (en) * 2007-02-07 2008-08-14 Solexa Limited Preparation of templates for methylation analysis
CN102409408B (en) * 2010-09-21 2013-08-07 深圳华大基因科技服务有限公司 Method for accurate detection of whole genome methylation sites by utilizing trace genome DNA (deoxyribonucleic acid)
WO2012058634A2 (en) * 2010-10-28 2012-05-03 Salk Institute For Biological Studies Epigenomic induced pluripotent stem cell signatures

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2003035860A1 (en) * 2001-10-24 2003-05-01 The Trustees Of Columbia University In The City Of New York A method for gene identification based on differential dna methylation
WO2003087774A3 (en) * 2002-04-12 2003-12-31 Univ Missouri Ecist microarrays for dual screening of dna hypermethylation and gene silencing
WO2005090607A1 (en) * 2004-03-08 2005-09-29 Rubicon Genomics, Inc. Methods and compositions for generating and amplifying dna libraries for sensitive detection and analysis of dna methylation
CN101802223A (en) * 2007-08-15 2010-08-11 香港大学 Methods and compositions for high throughput bisulfite DNA-sequencing and uses thereof

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Direct selection of human genomic loci by microarray hybridization;Thomas J Albert et al.;《NATURE METHODS》;20071130;第4卷(第11期);903页右栏第2、3段,图1 *
Thomas J Albert et al..Direct selection of human genomic loci by microarray hybridization.《NATURE METHODS》.2007,第4卷(第11期),

Also Published As

Publication number Publication date
CN102796808A (en) 2012-11-28
WO2012159564A1 (en) 2012-11-29
DK2725125T3 (en) 2018-04-23
HK1178211A1 (en) 2013-09-06
EP2725125B1 (en) 2018-01-17
US20140256563A1 (en) 2014-09-11
EP2725125A1 (en) 2014-04-30
US9133513B2 (en) 2015-09-15
EP2725125A4 (en) 2014-12-24

Similar Documents

Publication Publication Date Title
CN102796808B (en) Methylation high-flux detection method
CN104250663B (en) The high-flux sequence detection method on methylated CpG island
CN103088433B (en) Construction method and application of genome-wide methylation high-throughput sequencing library
CN105463585B (en) Method and its application based on single strand dna structure sequencing library
CN105400776B (en) Oligonucleotide adapters and their application in the construction of single-strand circular libraries for nucleic acid sequencing
CN102409042B (en) High throughput genome methylated DNA enrichment method and indexes and index adapters used in same
CN102732629B (en) Method for concurrently determining gene expression level and polyadenylic acid tailing by using high-throughput sequencing
JP6899844B2 (en) Methods and kits for generating DNA libraries for massively parallel sequencing
CN104508144A (en) Methods and systems for determining and phasing haplotypes
CN103806111A (en) Construction method and application of high-throughout sequencing library
US10415083B2 (en) Long insert-based whole genome sequencing
CN105734679B (en) Nucleic acid target sequence captures the preparation method of sequencing library
CN102839168A (en) Nucleic acid probe, and preparation method and application thereof
CN110760936A (en) Method for constructing DNA methylation library and its application
CN111378720A (en) Construction method and application of sequencing library of long-chain non-coding RNA
US11261479B2 (en) Methods and compositions for enrichment of target nucleic acids
CN110741096A (en) Compositions and methods for detecting circulating tumor DNA
CN103374759B (en) A kind of detection of lung cancer shifts method and the application thereof of significant SNP
CN108130366A (en) A kind of method for building people miRNA sequencing libraries and carrying out high-flux sequence
CN111748621A (en) A probe library, kit and application for detecting lung cancer-related 41 genes
US20220325317A1 (en) Methods for generating a population of polynucleotide molecules
US20250163407A1 (en) Methods selectively depleting nucleic acid using rnase h
CN108103173A (en) A kind of method for building mouse miRNA sequencing libraries and carrying out high-flux sequence
HK1178211B (en) High throughput methylation detection method
HK1169458A (en) Ffpe sample nucleic acid library, its construction and the ffpe sample analysis method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
ASS Succession or assignment of patent right

Owner name: BGI TECHNOLOGY SOLUTIONS CO., LTD.

Free format text: FORMER OWNER: BGI-SHENZHEN CO., LTD.

Effective date: 20130715

Free format text: FORMER OWNER: BGI-SHENZHEN

Effective date: 20130715

C41 Transfer of patent application or patent right or utility model
TA01 Transfer of patent application right

Effective date of registration: 20130715

Address after: 518083 science and Technology Pioneer Park, comprehensive building, Beishan Industrial Zone, Yantian District, Guangdong, Shenzhen 201

Applicant after: BGI Technology Solutions Co., Ltd.

Address before: North Road No. 146, building 11F-3 Industrial Zone in Yantian District of Shenzhen city of Guangdong Province in 518083

Applicant before: BGI-Shenzhen Co., Ltd.

Applicant before: BGI-Shenzhen

REG Reference to a national code

Ref country code: HK

Ref legal event code: DE

Ref document number: 1178211

Country of ref document: HK

C14 Grant of patent or utility model
GR01 Patent grant
REG Reference to a national code

Ref country code: HK

Ref legal event code: GR

Ref document number: 1178211

Country of ref document: HK