US20110151438A9 - Methods of Analysis of Methylation - Google Patents

Methods of Analysis of Methylation Download PDF

Info

Publication number
US20110151438A9
US20110151438A9 US11/923,649 US92364907A US2011151438A9 US 20110151438 A9 US20110151438 A9 US 20110151438A9 US 92364907 A US92364907 A US 92364907A US 2011151438 A9 US2011151438 A9 US 2011151438A9
Authority
US
United States
Prior art keywords
sequence
methylation
probes
probe
complementary
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.)
Abandoned
Application number
US11/923,649
Other versions
US20080108073A1 (en
Inventor
Shivani Nautiyal
Malek Faham
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.)
Affymetrix Inc
Original Assignee
Affymetrix Inc
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
Priority claimed from US10/300,311 external-priority patent/US7208295B2/en
Application filed by Affymetrix Inc filed Critical Affymetrix Inc
Priority to US11/923,649 priority Critical patent/US20110151438A9/en
Assigned to AFFYMETRIX, INC. reassignment AFFYMETRIX, INC. ASSIGNMENT OF ASSIGNORS INTEREST (SEE DOCUMENT FOR DETAILS). Assignors: FAHAM, MALEK, NAUTIYAL, SHIVANI
Publication of US20080108073A1 publication Critical patent/US20080108073A1/en
Publication of US20110151438A9 publication Critical patent/US20110151438A9/en
Assigned to GENERAL ELECTRIC CAPITAL CORPORATION, AS AGENT reassignment GENERAL ELECTRIC CAPITAL CORPORATION, AS AGENT SECURITY AGREEMENT Assignors: AFFYMETRIX, INC.
Assigned to AFFYMETRIX, INC. reassignment AFFYMETRIX, INC. RELEASE BY SECURED PARTY (SEE DOCUMENT FOR DETAILS). Assignors: GENERAL ELECTRIC CAPITAL CORPORATION, AS AGENT
Priority to US15/281,721 priority patent/US10407717B2/en
Priority to US16/528,222 priority patent/US10822642B2/en
Priority to US17/080,813 priority patent/US20210147903A1/en
Abandoned legal-status Critical Current

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/6813Hybridisation assays
    • C12Q1/6827Hybridisation assays for detection of mutation or polymorphism
    • C12Q1/683Hybridisation assays for detection of mutation or polymorphism involving restriction enzymes, e.g. restriction fragment length polymorphism [RFLP]
    • 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/6813Hybridisation assays
    • 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/6844Nucleic acid amplification reactions
    • C12Q1/6848Nucleic acid amplification reactions characterised by the means for preventing contamination or increasing the specificity or sensitivity of an amplification reaction
    • 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/6844Nucleic acid amplification reactions
    • C12Q1/6853Nucleic acid amplification reactions using modified primers or templates
    • 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/6844Nucleic acid amplification reactions
    • C12Q1/6853Nucleic acid amplification reactions using modified primers or templates
    • C12Q1/6855Ligating adaptors
    • 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
    • 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
    • C12Q2521/00Reaction characterised by the enzymatic activity
    • C12Q2521/30Phosphoric diester hydrolysing, i.e. nuclease
    • C12Q2521/331Methylation site specific nuclease
    • 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
    • C12Q2523/00Reactions characterised by treatment of reaction samples
    • C12Q2523/10Characterised by chemical treatment
    • C12Q2523/125Bisulfite(s)
    • 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
    • C12Q2525/00Reactions involving modified oligonucleotides, nucleic acids, or nucleotides
    • C12Q2525/10Modifications characterised by
    • C12Q2525/155Modifications characterised by incorporating/generating a new priming site
    • 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
    • C12Q2565/00Nucleic acid analysis characterised by mode or means of detection
    • C12Q2565/50Detection characterised by immobilisation to a surface
    • C12Q2565/501Detection characterised by immobilisation to a surface being an array of oligonucleotides
    • 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
    • C12Q2600/00Oligonucleotides characterized by their use
    • C12Q2600/154Methylation markers

Definitions

  • genomic DNA is modified by bisulfite. Fragments of the modified DNA are generated with defined ends using locus specific primer extension.
  • FIG. 1 is a schematic representation illustrating methods for appending known sequences to a single-stranded nucleic acid template at specific positions.
  • a template with defined ends is obtained by primer extension.
  • FIG. 1B common sequences are ligated to the ends of the template.
  • FIG. 1C shows an alternate embodiment for digestion of the template probe.
  • FIG. 7 shows a method for determining methylation state of a restriction site using dU probes.
  • FIG. 7A shows differential cleavage and end filling steps for determining methylation using dU probes that distinguish between methylated and unmethylated restriction sites.
  • FIG. 8 shows 384 multiplex amplification using dU probes.
  • Computer software products of the invention typically include computer readable medium having computer-executable instructions for performing the logic steps of the method of the invention.
  • Suitable computer readable medium include floppy disk, CD-ROM/DVD/DVD-ROM, hard-disk drive, flash memory, ROM/RAM, magnetic tapes and etc.
  • the computer executable instructions may be written in a suitable computer language or combination of several languages.
  • coli DNA ligase which catalyzes the formation of a phosphodiester bond between juxtaposed 5′-phosphate and 3′-hydroxyl termini in duplex DNA containing cohesive ends; and T4 RNA ligase which catalyzes ligation of a 5′ phosphoryl-terminated nucleic acid donor to a 3′ hydroxyl-terminated nucleic acid acceptor through the formation of a 3′->5′ phosphodiester bond
  • substrates include single-stranded RNA and DNA as well as dinucleoside pyrophosphates; or any other methods described in the art.
  • Fragmented DNA may be treated with one or more enzymes, for example, an endonuclease, prior to ligation of adaptors to one or both ends to facilitate ligation by generating ends that are compatible with ligation.
  • genomic is all the genetic material in the chromosomes of an organism.
  • DNA derived from the genetic material in the chromosomes of a particular organism is genomic DNA.
  • a genomic library is a collection of clones made from a set of randomly generated overlapping DNA fragments representing the entire genome of an organism.
  • Probe 130 includes at least first target complementarity region 120 and oligonucleotide positioning regions 119 and 121 which are directly adjacent to 120 .
  • template region 113 hybridizes to template complementarity region 120 of probe 130 and oligonucleotide regions 15 and 17 concurrently hybridize to oligonucleotide positioning regions 119 and 121 of probe 130 .
  • Bisulfite treatment allows the methylation status of cytosines to be detected by a variety of methods.
  • any method that may be used to detect a SNP may be used, for examples, see Syvanen, Nature Rev. Gen. 2:930-942 (2001).
  • Methods such as single base extension (SBE) may be used or hybridization of sequence specific probes similar to allele specific hybridization methods.
  • SBE single base extension
  • MIP Molecular Inversion Probe
  • molecular inversion probes may be used to determine methylation status after methylation dependent modification.
  • a MIP may be designed for each cytosine to be interrogated.
  • the MIP includes a locus specific region that hybridizes upstream and one that hybridizes downstream of an interrogation site and can be extended through the interrogation site, incorporating a base that is complementary to the interrogation position.
  • the interrogation position may be the cytosine of interest after bisulfite modification and amplification of the region and the detection can be similar to detection of a polymorphism.
  • Separate reactions may be performed for each NTP so extension only takes place in the reaction containing the base corresponding to the interrogation base or the different products may be differentially labeled.
  • methods for identifying methylation based on differential cleavage by restriction enzymes are used.
  • Methylation-sensitive restriction analysis followed by PCR amplification or Southern analysis have been disclosed, for example, in Huang et al. Cancer Res. 57:1030-1034 (1997), Zuccotti et al, Methods in Enzym. 225:557-567 (1993), Carrel et al. Am Jour. Med. Genet. 64:27-30 (1996) and Chang et al. Plant MoL Bio. Rep. 10:362-366 (1992).
  • Enzymes with a recognition site that includes the sequence GCGC include, for example, Hin P1 I, Hha I, Afe I, Kas I, Nar I, Sfo I, Bbe I, and Fsp I.
  • Enzymes with a recognition site that includes the sequence TCGA include, for example, Taq I, Cla I (MSRE), BspD I (MSRE), PaeR7 I, Tli I, Xho I, Sal I, and BstB I.
  • FIG. 3 illustrates an embodiment in which priming sites containing methylated cytosines are added to genomic DNA using uracil containing target specific template probes 309 .
  • the uracil containing template probes each have a target specific region 313 and 314 that is flanked by common priming sites 311 and 312 .
  • the common priming sites may be the same or different.
  • Oligonucleotide 315 is complementary to 311 and 312 and contains methylated cytosines that are resistant to bisulfite conversion.
  • the oligonucleotide 315 and the genomic targets 305 and 307 hybridize to the dU probes to form duplexes 317 .
  • FIG. 4 illustrates the alteration of DNA sequence resulting from treatment of DNA with bisulfite.
  • SEQ ID NOs. 1 and 2 are the complementary starting strands. Methylated cytosines are indicated by “m”.
  • SEQ ID NO. 3 corresponds to SEQ ID NO. 1 after treatment and SEQ ID NO. 4 corresponds to SEQ ID NO. 2 after treatment. Unmethylated cytosines have been converted to “U” and methylated cytosines remain “C”.
  • SEQ ID NO. 3 gives SEQ ID NOs. 5 and 6, while SEQ ID NO. 4 results in SEQ ID NOs. 7 and 8. Many of the G-C base pairs from the starting sequence have been converted to A-T base pairs; resulting in a relatively low GC content for the final sequence.
  • FIG. 5 illustrates a method for methylation analysis of a selected group of target sequences using dU probe technology.
  • Genomic DNA [ 501 ] is digested with one or more restriction enzymes to generate restriction fragments [ 503 ].
  • the restriction fragments are mixed with a collection of dU probes for specific targets [ 504 ] and common primer sequences [ 505 ] and the nucleic acids are denatured and allowed to anneal.
  • the annealing step results in the formation of complexes between the dU probes, target fragments and primers [ 506 ].
  • FIG. 6 illustrates a method of determining the methylation state of a cytosine in a methylation sensitive restriction site.
  • Fragment 601 has an unmethylated MspI/HpaII site while fragment 602 has a methylated site.
  • MspI and HpaII are isoschizomers with MspI being methylation insensitive and HpaII being methylation sensitive.
  • the DNA is cleaved with HpaII so that unmethylated sites are cleaved while methylated sites remain uncleaved.
  • the fragments are then end filled using Klenow.
  • the DNA is then cleaved with the MspI so that the methylated sites are cleaved.
  • Probes 604 and 605 can hybridize to either the product generated from the methylated or unmethylated sites but ligation between the 603 a and 606 is blocked and ligation of 603 b and 607 is blocked.
  • the reactions on the left have ligase while the reactions on the right are without ligase. As expected, amplification products are only observed in the presence of ligase.
  • Each of the different buffer conditions gave specific amplification products in the expected size range. Buffer A gives more amplification product than Buffer B at lower dU probe concentrations. The reactions have increasing amounts of dU probe as indicated graphically above each gel image.
  • one or more reaction conditions may be varied.
  • the basic starting conditions are 6 ⁇ SSPE with 7.5% PEG
  • the annealing conditions are 94° C. for 1 min, 70° C. for 40 min, 65° C. for 40 min, 60° C. for 40 min, 55° C. for 10 min, and 50° C. for 10 min.
  • ranges of conditions that may be used include 0-15 mM MgCl 2 , 0-200 mM KCl, variation of the pH, addition of RecA and SSB, including the use of thermophilic SSB, and variation of the reaction volume.
  • AP sites may then be cleaved enzymatically using an AP endonuclease or, under certain conditions, an AP lyase.
  • the AP site may be cleaved using ApeI, an Endo IV or Fpg (formamidopyrimidine [fapy]-DNA glycosylase; also known as 8-oxoguanine DNA glycosylase).
  • Fpg cleaves both 3′ and 5′ to the AP site, removing the AP site and leaving a 1 base gap.
  • the AP sites may be cleaved chemically, such as by treatment with 1,4 diaminobutane and heat.
  • AP sites may also be cleaved by exposure to high pH.
  • the probe may be amplified after circularization. Circularized probes that have been linearized may be amplified, for example by PCR using primers to regions 502 and 506 . Other methods of amplification may also be used. For example, rolling circle amplification may be used to generate multiple copies of the circularized probes. See, for example, U.S. Pat. Nos. 5,648,245 and 5,854,033 and Fire and Xu, Proc Natl Acad Sci USA 92:4641-5 (1995). As used herein, the term amplification includes the production of RNA transcripts by polymerization driven from a phage promoter. For example, a T7 RNA polymerase promoter sequence may be incorporated into the probe and used to generate multiple RNA copies using a T7 RNA polymerase.
  • Nucleic acid amplification methods useful in the methods of the present invention are well known in the art and include, e.g., polymerase chain reaction (PCR), nucleic acid sequence-based amplification (NASBA), self-sustained sequence recognition (3 SR), ligase chain reaction (LCR), transcription-mediated amplification (TMA), rolling circle amplification (RCA), and strand displacement amplification (SDA).
  • PCR polymerase chain reaction
  • NASBA nucleic acid sequence-based amplification
  • SR self-sustained sequence recognition
  • LCR ligase chain reaction
  • TMA transcription-mediated amplification
  • RCA rolling circle amplification
  • SDA strand displacement amplification

Abstract

Methods for determining the methylation status of a plurality of cytosines are disclosed. In some aspects genomic DNA target sequences containing CpGs are targeted for analysis by multiplex amplification using target specific probes that can be specifically degraded prior to amplification. The targets may be modified with bisulfite prior to amplification. In another aspect targets are cut with methylation sensitive or insensitive restriction enzymes and marked with a tag using the target specific probes. The presence or absence of methylation may be determined using methylation sensitive restriction enzyme or bisulfite treatment. Detection in many embodiments employs hybridization to tag arrays, genotyping arrays or resequencing arrays.

Description

    RELATED APPLICATIONS
  • The present application claims priority to U.S. application No. 60/862,735 filed Oct. 24, 2006, the disclosure of which is incorporated herein by reference in its entirety.
  • BACKGROUND OF THE INVENTION
  • The genomes of higher eukaryotes contain the modified nucleoside 5-methyl cytosine (5-meC). This modification is usually found as part of the dinucleotide CpG.
  • Cytosine is converted to 5-methylcytosine in a reaction that involves flipping a target cytosine out of an intact double helix and transfer of a methyl group from S-adenosylmethionine by a methyltransferase enzyme (Klimasauskas et al. Cell 76:357-369, 1994). This enzymatic conversion is the only epigenetic modification of DNA known to exist in vertebrates and is essential for normal embryonic development (Bird, Cell 70:5-8, 1992; Laird and Jaenisch, Human Mol. Genet. 3:1487-1495, 1994; and Li et al. Cell 69:915-926, 1992).
  • The frequency of the CpG dinucleotide in the human genome is only about 20% of the statistically expected frequency, possibly because of spontaneous deamination of 5-meC to T (Schoreret et al. Proc. Natl. Acad Sci. USA 89:957-961, 1992). There are about 28 million CpG doublets in a haploid copy of the human genome and it is estimated that about 70-80% of the cytosines at CpGs are methylated. Regions where CpG is present at levels that are approximately the expected frequency are referred to as “CpG islands” (Bird, A. P., Nature 321:209-213, 1986). These regions have been estimated to comprise about 1% of vertebrate genomes and account for about 15% of the total number of CpG dinucleotides. CpG islands are typically between 0.2 and 1 kb in length and are often located upstream of housekeeping and tissue-specific genes. CpG islands are often located upstream of transcribed regions, but may also extend into transcribed regions. About 2-4% of cytosines are methylated and probably the majority of cytosines that are 5′ of Gs are methylated. Most of the randomly distributed CpGs are methylated, but only about 20% of the CpGs in CpG islands are methylated. Recent studies on CpG islands suggest that promoters segregate into two classes by CpG content. See, Saxonov et al. PNAS 103(5):1412-7 (2006).
  • DNA methylation is an epigenetic determinant of gene expression. Patterns of CpG methylation are heritable, tissue specific, and correlate with gene expression.
  • The consequence of methylation is usually gene silencing. DNA methylation also correlates with other cellular processes including embryonic development, chromatin structure, genomic imprinting, somatic X-chromosome inactivation in females, inhibition of transcription and transposition of foreign DNA and timing of DNA replication. When a gene is highly methylated it is less likely to be expressed, possibly because CpG methylation prevents transcription factors from recognizing their cognate binding sites. Proteins that bind methylated DNA may also recruit histone deacetylase to condense adjacent chromatin. Such “closed” chromatin structures prevent binding of transcription factors. Thus the identification of sites in the genome containing 5-meC is important in understanding cell-type specific programs of gene expression and how gene expression profiles are altered during both normal development and diseases such as cancer. Precise mapping of DNA methylation patterns in CpG islands has become essential for understanding diverse biological processes such as the regulation of imprinted genes, X chromosome inactivation, and tumor suppressor gene silencing in human cancer caused by increase methylation.
  • Methylation of cytosine residues in DNA plays an important role in gene regulation. Methylation of cytosine may lead to decreased gene expression by, for example, disruption of local chromatin structure, inhibition of transcription factor-DNA binding, or by recruitment of proteins which interact specifically with methylated sequences and prevent transcription factor binding. DNA methylation is required for normal embryonic development and changes in methylation are often associated with disease. Genomic imprinting, X chromosome inactivation, chromatin modification, and silencing of endogenous retroviruses all depend on establishing and maintaining proper methylation patterns. Abnormal methylation is a hallmark of cancer cells and silencing of tumor suppressor genes is thought to contribute to carcinogenesis. Methylation mapping using microarray-based approaches may be used, for example, to profile cancer cells revealing a pattern of DNA methylation that may be used, for example, to diagnose a malignancy, predict treatment outcome or monitor progression of disease. Methylation in eukaryotes can also function to inhibit the activity of viruses and transposons, see Jones et al. EMBO J. 17:6385-6393 (1998). Alterations in the normal methylation process have also been shown to be associated with genomic instability (Lengauer et al. Proc. Natl. Acad. Sci. USA 94:2545-2550, 1997). Such abnormal epigenetic changes may be found in many types of cancer and can serve as potential markers for oncogenic transformation.
  • SUMMARY OF THE INVENTION
  • Methods for analyzing the methylation status of cytosines in genomic DNA are disclosed.
  • In some aspects the methods include a step of multiplex amplification of a plurality of regions of interest. The methods provide for the addition of known priming sequences to the 5′ and 3′ ends of the sequences to be amplified so that subsequent amplification may be performed using primers to the known priming sequences. Such multiplexed amplification reactions provide high specificity and uniform amplification of templates.
  • In a first aspect, the invention provides a method for multiplex locus specific amplification of a plurality of templates to provide a plurality of templates with known 5′ and 3′ ends.
  • The template may be derived from cDNA or genomic DNA, from a single individual or from a plurality of individuals. The template may, for example, be genomic DNA derived from a eukaryote, such as a human being.
  • The multiplex methods of the present invention may include at least 10 templates of distinct sequence, at least 100 templates of distinct sequence, at least 1000 templates of distinct sequence, or more. Usefully, at least one of the first and second oligonucleotides comprises a bar code sequence, thus allowing concurrent detection of all amplified templates.
  • In one embodiment, genomic DNA is modified by bisulfite. Fragments of the modified DNA are generated with defined ends using locus specific primer extension.
  • The extension products have defined ends and are then hybridized to a dU probe and adaptor sequences are ligated to the ends. The dU probe is degraded and the adaptor ligated fragments are amplified.
  • In some aspects the methods are used to classify a tissue into a class, for example, a known tumor class. The hybridization pattern obtained from the tissue sample, using the disclosed methods, is compared to hybridization patterns from samples from tissues of known tumor class, obtained using the disclosed methods.
  • In one aspect a method for analyzing the methylation of a plurality of cytosines in a plurality of target sequences is disclosed. A genomic DNA sample is fragmented to generate fragments that include a mixture of target fragments and non-target fragments. The fragments are mixed with a common primer sequence and a collection of dU probes that are complementary to different target sequences to be analyzed. Each dU probe has a sequence that is complementary to a different target fragment flanked at both ends by the complement of the common primer sequence. Target fragments and common primer sequences hybridize to dU probes to form ligation complexes and ligase is added to ligate the common primer sequences to the target fragments in the ligation complexes. The dU probes are digested using UDG and the ligated products are treated with bisulfite and amplified. The amplification product is analyzed, for example, by hybridization to an array to determine the methylation state of cytosines in the starting sample by detecting sequence changes corresponding to bisulfite modification.
  • In another aspect DNA is fragmented with a methylation sensitive enzyme so that only unmethylated DNA is fragmented. The overhang created by cleavage is filled in with a DNA polymerase, marking the unmethylated fragments with an additional sequence. The DNA is then fragmented with an isoschizomer of the first enzyme that is methylation insensitive so it will cleave the methylated sites. The fragments are then hybridized to dU probes that are designed to hybridize to either the fragment generated by cleavage then filling or cleavage alone. Amplification products are generated that are differentially detectable, for example, by being marked with different tag sequences.
  • BRIEF DESCRIPTION OF THE DRAWINGS
  • The above and other objects and advantages of the present invention will be apparent upon consideration of the following detailed description taken in conjunction with the accompanying drawings, in which like characters refer to like parts throughout, and in which:
  • FIG. 1 is a schematic representation illustrating methods for appending known sequences to a single-stranded nucleic acid template at specific positions. In FIG. 1A a template with defined ends is obtained by primer extension. In FIG. 1B common sequences are ligated to the ends of the template. FIG. 1C shows an alternate embodiment for digestion of the template probe.
  • FIG. 2A illustrates one embodiment of a dU probe.
  • FIG. 2B illustrates the use of dU probes for multiplex amplification of specific target sequences from a complex sample.
  • FIG. 3 is a schematic representation illustrating a method for using dU probes in conjunction with bisulfite modification to analyze methylation of a collection of pre-selected targets.
  • FIG. 4 shows the changes resulting from treatment of DNA with bisulfite followed by PCR amplification.
  • FIG. 5 illustrates a method of using dU probes in combination with bisulfite treatment for methylation analysis.
  • FIG. 6 is a schematic representation illustrating a method for determining if a restriction site is methylated. FIG. 6A illustrates methyl sensitive digestion and hybridization of methylation specific dU probes. FIG. 6B shows how differentially tagged fill-in and non fill specific dU probes are used to detect presence or absence of methylation.
  • FIG. 7 shows a method for determining methylation state of a restriction site using dU probes. FIG. 7A shows differential cleavage and end filling steps for determining methylation using dU probes that distinguish between methylated and unmethylated restriction sites.
  • FIG. 7B shows differential dU probes as template for ligation of tags and primers to different cleavage products from FIG. 7A.
  • FIG. 7C shows a schematic and a gel image of fragments resulting from restriction enzyme cleavage, dU probe hybridization and template directed ligation, followed by PCR amplification of 9 targets.
  • FIG. 8 shows 384 multiplex amplification using dU probes.
  • FIG. 9 shows a gel image of dU capture following cleavage using fill-in and no fill specific probes as illustrated schematically in FIG. 7
  • FIG. 10 is a schematic representation illustrating embodiments for determining the epihaplotype of two neighboring cytosines.
  • FIG. 11 is a schematic representation illustrating the use of molecular inversion probes.
  • FIG. 12 is a schematic illustrating probe inversion, amplification, and detection of a tagged fragment.
  • DETAILED DESCRIPTION OF THE INVENTION
  • General
  • The present invention has many preferred embodiments and relies on many patents, applications and other references for details known to those of the art. Therefore, when a patent, application, or other reference is cited or repeated below, it should be understood that it is incorporated by reference in its entirety for all purposes as well as for the proposition that is recited.
  • As used in this application, the singular form “a,” “an,” and “the” include plural references unless the context clearly dictates otherwise. For example, the term “an agent” includes a plurality of agents, including mixtures thereof.
  • An individual is not limited to a human being, but may also include other organisms including but not limited to mammals, plants, fungi, bacteria or cells derived from any of the above.
  • Throughout this disclosure, various aspects of this invention can be presented in a range format. It should be understood that the description in range format is merely for convenience and brevity and should not be construed as an inflexible limitation on the scope of the invention. Accordingly, the description of a range should be considered to have specifically disclosed all the possible subranges as well as individual numerical values within that range. For example, description of a range such as from 1 to 6 should be considered to have specifically disclosed subranges such as from 1 to 3, from 1 to 4, from 1 to 5, from 2 to 4, from 2 to 6, from 3 to 6 etc., as well as individual numbers within that range, for example, 1, 2, 3, 4, 5, and 6. This applies regardless of the breadth of the range.
  • The practice of the present invention may employ, unless otherwise indicated, conventional techniques and descriptions of organic chemistry, polymer technology, molecular biology (including recombinant techniques), cell biology, biochemistry, and immunology, which are within the skill of the art. Such conventional techniques include polymer array synthesis, hybridization, ligation, and detection of hybridization using a label. Specific illustrations of suitable techniques can be had by reference to the example herein below. However, other equivalent conventional procedures can, of course, also be used. Such conventional techniques and descriptions can be found in standard laboratory manuals such as Genome Analysis: A Laboratory Manual Series (Vols. I-IV), Using Antibodies: A Laboratory Manual, Cells: A Laboratory Manual, PCR Primer: A Laboratory Manual, and Molecular Cloning: A Laboratory Manual (all from Cold Spring Harbor Laboratory Press), Stryer, L. (1995) Biochemistry (4th Ed.) Freeman, New York, Gait, “Oligonucleotide Synthesis: A Practical Approach” 1984, IRL Press, London, Nelson and Cox (2000), Lehninger, Principles of Biochemistry 3rd Ed., W.H. Freeman Pub., New York, N.Y. and Berg et al. (2002) Biochemistry, 5th Ed., W.H. Freeman Pub., New York, N.Y., all of which are herein incorporated in their entirety by reference for all purposes.
  • The present invention can employ solid substrates, including arrays in some preferred embodiments. Methods and techniques applicable to polymer (including protein) array synthesis have been described in U.S. Ser. No. 09/536,841, WO 00/58516, U.S. Pat. Nos. 5,143,854, 5,242,974, 5,252,743, 5,324,633, 5,384,261, 5,405,783, 5,424,186, 5,451,683, 5,482,867, 5,491,074, 5,527,681, 5,550,215, 5,571,639, 5,578,832, 5,593,839, 5,599,695, 5,624,711, 5,631,734, 5,795,716, 5,831,070, 5,837,832, 5,856,101, 5,858,659, 5,936,324, 5,968,740, 5,974,164, 5,981,185, 5,981,956, 6,025,601, 6,033,860, 6,040,193, 6,090,555, 6,136,269, 6,269,846 and 6,428,752, in PCT Applications Nos. PCT/US99/00730 (International Publication No. WO 99/36760) and PCT/US01/04285 (International Publication No. WO 01/58593), which are all incorporated herein by reference in their entirety for all purposes.
  • Patents that describe synthesis techniques in specific embodiments include U.S. Pat. Nos. 5,412,087, 6,147,205, 6,262,216, 6,310,189, 5,889,165, and 5,959,098. Nucleic acid arrays are described in many of the above patents, but the same techniques are applied to polypeptide arrays.
  • Nucleic acid arrays that are useful in the present invention include those that are commercially available from Affymetrix (Santa Clara, Calif.) under the brand name GeneChip®. Example arrays are shown on the website at affymetrix.com.
  • The present invention also contemplates many uses for polymers attached to solid substrates. These uses include gene expression monitoring, profiling, library screening, genotyping and diagnostics. Gene expression monitoring and profiling methods can be shown in U.S. Pat. Nos. 5,800,992, 6,013,449, 6,020,135, 6,033,860, 6,040,138, 6,177,248 and 6,309,822. Genotyping and uses therefore are shown in U.S. Ser. Nos. 10/442,021, 10/013,598 (U.S. Patent Application Publication 20030036069), and U.S. Pat. Nos. 5,856,092, 6,300,063, 5,858,659, 6,284,460, 6,361,947, 6,368,799 and 6,333,179. Other uses are embodied in U.S. Pat. Nos. 5,871,928, 5,902,723, 6,045,996, 5,541,061, and 6,197,506.
  • The present invention also contemplates sample preparation methods in certain preferred embodiments. Prior to or concurrent with hybridization to an array, the sample may be amplified by a variety of mechanisms, some of which may employ PCR. See, for example, PCR Technology: Principles and Applications for DNA Amplification (Ed. H.A. Erlich, Freeman Press, NY, N.Y., 1992); PCR Protocols: A Guide to Methods and Applications (Eds. Innis, et al., Academic Press, San Diego, Calif., 1990); Mattila et al., Nucleic Acids Res. 19, 4967 (1991); Eckert et al., PCR Methods and Applications 1, 17 (1991); PCR (Eds. McPherson et al., IRL Press, Oxford); and U.S. Pat. Nos. 4,683,202, 4,683,195, 4,800,159, 4,965,188, and 5,333,675. The sample may be amplified on the array. See, for example, U.S. Pat. No. 6,300,070 which is incorporated herein by reference.
  • Other suitable amplification methods include the ligase chain reaction (LCR) (for example, Wu and Wallace, Genomics 4, 560 (1989), Landegren et al., Science 241, 1077 (1988) and Barringer et al. Gene 89:117 (1990)), transcription amplification (Kwoh et al., Proc. Natl. Acad. Sci. USA 86, 1173 (1989) and WO88/10315), self-sustained sequence replication (Guatelli et al., Proc. Nat. Acad. Sci. USA, 87, 1874 (1990) and WO90/06995), selective amplification of target polynucleotide sequences (U.S. Pat. No. 6,410,276), consensus sequence primed polymerase chain reaction (CP-PCR) (U.S. Pat. No. 4,437,975), arbitrarily primed polymerase chain reaction (AP-PCR) (U.S. Pat. Nos. 5,413,909, 5,861,245), rolling circle amplification (RCA) (for example, Fire and Xu, PNAS 92:4641 (1995) and Liu et al. J. Am. Chem. Soc. 118:1587 (1996)) and nucleic acid based sequence amplification (NABSA), (See, U.S. Pat. Nos. 5,409,818, 5,554,517, and 6,063,603).
  • Other amplification methods that may be used are described in, U.S. Pat. Nos. 5,242,794, 5,494,810, 4,988,617 and in U.S. Ser. No. 09/854,317. Other amplification methods are also disclosed in Dahl et al., Nuc. Acids Res. 33(8):e71 (2005) and circle to circle amplification (C2CA) Dahl et al., PNAS 101:4548 (2004). Locus specific amplification and representative genome amplification methods may also be used.
  • Additional methods of sample preparation and techniques for reducing the complexity of a nucleic sample are described in Dong et al. Genome Research 11, 1418 (2001), in U.S. Pat. Nos. 6,872,529, 6,361,947, 6,391,592 and 6,107,023, US Patent Publication Nos. 20030096235 and 20030082543 and U.S. patent application Ser. No. 09/916,135.
  • Methods for conducting polynucleotide hybridization assays have been well developed in the art. Hybridization assay procedures and conditions will vary depending on the application and are selected in accordance with the general binding methods known including those referred to in: Maniatis et al. Molecular Cloning: A Laboratory Manual (2nd Ed. Cold Spring Harbor, N.Y, 1989); Berger and Kimmel Methods in Enzymology, Vol. 152, Guide to Molecular Cloning Techniques (Academic Press, Inc., San Diego, Calif., 1987); Young and Davism, P.N.A.S, 80: 1194 (1983). Methods and apparatus for carrying out repeated and controlled hybridization reactions have been described in U.S. Pat. Nos. 5,871,928, 5,874,219, 6,045,996 and 6,386,749, 6,391,623 each of which are incorporated herein by reference.
  • The present invention also contemplates signal detection of hybridization between ligands in certain preferred embodiments. See U.S. Pat. Nos. 5,143,854, 5,578,832; 5,631,734; 5,834,758; 5,936,324; 5,981,956; 6,025,601; 6,141,096; 6,185,030; 6,201,639; 6,218,803; and 6,225,625, in U.S. Ser. No. 10/389,194 and in PCT Application PCT/US99/06097 (published as WO99/47964), each of which also is hereby incorporated by reference in its entirety for all purposes.
  • Methods and apparatus for signal detection and processing of intensity data are disclosed in, for example, U.S. Pat. Nos. 5,143,854, 5,547,839, 5,578,832, 5,631,734, 5,800,992, 5,834,758; 5,856,092, 5,902,723, 5,936,324, 5,981,956, 6,025,601, 6,090,555, 6,141,096, 6,185,030, 6,201,639; 6,218,803; and 6,225,625, in U.S. Ser. Nos. 10/389,194, 60/493,495 and in PCT Application PCT/US99/06097 (published as WO 99/47964), each of which also is hereby incorporated by reference in its entirety for all purposes. Instruments and software may also be purchased commercially from various sources, including Affymetrix.
  • The practice of the present invention may also employ conventional biology methods, software and systems. Computer software products of the invention typically include computer readable medium having computer-executable instructions for performing the logic steps of the method of the invention. Suitable computer readable medium include floppy disk, CD-ROM/DVD/DVD-ROM, hard-disk drive, flash memory, ROM/RAM, magnetic tapes and etc. The computer executable instructions may be written in a suitable computer language or combination of several languages. Basic computational biology methods are described in, for example Setubal and Meidanis et al., Introduction to Computational Biology Methods (PWS Publishing Company, Boston, 1997); Salzberg, Searles, Kasif, (Ed.), Computational Methods in Molecular Biology, (Elsevier, Amsterdam, 1998); Rashidi and Buehler, Bioinformatics Basics: Application in Biological Science and Medicine (CRC Press, London, 2000) and Ouelette and Bzevanis Bioinformatics: A Practical Guide for Analysis of Gene and Proteins (Wiley & Sons, Inc., 2nd ed., 2001). See U.S. Pat. No. 6,420,108.
  • Methods for detection of methylation status are disclosed, for example, in Fraga and Esteller, BioTechniques 33:632-649 (2002) and Dahl and Guldberg Biogerontology 4:233-250 (2003). Methylation detection using bisulfite modification and target specific PCR have been disclosed, for example, in U.S. Pat. Nos. 5,786,146, 6,200,756, 6,143,504, 6,265,171, 6,251,594, 6,331,393, and 6,596,493. U.S. Pat. No. 6,884,586 disclosed methods for methylation analysis using nicking agents and isothermal amplification.
  • The present invention may also make use of various computer program products and software for a variety of purposes, such as probe design, management of data, analysis, and instrument operation. See, U.S. Pat. Nos. 5,593,839, 5,795,716, 5,733,729, 5,974,164, 6,066,454, 6,090,555, 6,185,561, 6,188,783, 6,223,127, 6,229,911 and 6,308,170.
  • Additionally, the present invention may have preferred embodiments that include methods for providing genetic information over networks such as the Internet as shown in U.S. Ser. Nos. 10/197,621, 10/063,559 (United States Publication No. 20020183936), 10/065,856, 10/065,868, 10/328,818, 10/328,872, 10/423,403, and 60/482,389.
  • All documents, i.e., publications and patent applications, cited in this disclosure, including the foregoing, are incorporated by reference herein in their entireties for all purposes to the same extent as if each of the individual documents were specifically and individually indicated to be so incorporated by reference herein in its entirety.
  • DEFINITIONS
  • “Adaptor sequences” or “adaptors” are generally oligonucleotides of at least 5, 10, or 15 bases and preferably no more than 50 or 60 bases in length; however, they may be even longer, up to 100 or 200 bases. Adaptor sequences may be synthesized using any methods known to those of skill in the art. For the purposes of this invention they may, as options, comprise primer binding sites, recognition sites for endonucleases, common sequences and promoters. The adaptor may be entirely or substantially double stranded or entirely single stranded. A double stranded adaptor may comprise two oligonucleotides that are at least partially complementary. The adaptor may be phosphorylated or unphosphorylated on one or both strands.
  • Adaptors may be more efficiently ligated to fragments if they comprise a substantially double stranded region and a short single stranded region which is complementary to the single stranded region created by digestion with a restriction enzyme. For example, when DNA is digested with the restriction enzyme EcoRI the resulting double stranded fragments are flanked at either end by the single stranded overhang 5′-AATT-3′, an adaptor that carries a single stranded overhang 5′-AATT-3′ will hybridize to the fragment through complementarity between the overhanging regions. This “sticky end” hybridization of the adaptor to the fragment may facilitate ligation of the adaptor to the fragment but blunt ended ligation is also possible. Blunt ends can be converted to sticky ends using the exonuclease activity of the Klenow fragment. For example when DNA is digested with PvuII the blunt ends can be converted to a two base pair overhang by incubating the fragments with Klenow in the presence of dTTP and dCTP. Overhangs may also be converted to blunt ends by filling in an overhang or removing an overhang.
  • Methods of ligation will be known to those of skill in the art and are described, for example in Sambrook et at. (2001) and the New England BioLabs catalog both of which are incorporated herein by reference for all purposes. Methods include using T4 DNA Ligase which catalyzes the formation of a phosphodiester bond between juxtaposed 5′ phosphate and 3′ hydroxyl termini in duplex DNA or RNA with blunt and sticky ends; Taq DNA Ligase which catalyzes the formation of a phosphodiester bond between juxtaposed 5′ phosphate and 3′ hydroxyl termini of two adjacent oligonucleotides which are hybridized to a complementary target DNA; E. coli DNA ligase which catalyzes the formation of a phosphodiester bond between juxtaposed 5′-phosphate and 3′-hydroxyl termini in duplex DNA containing cohesive ends; and T4 RNA ligase which catalyzes ligation of a 5′ phosphoryl-terminated nucleic acid donor to a 3′ hydroxyl-terminated nucleic acid acceptor through the formation of a 3′->5′ phosphodiester bond, substrates include single-stranded RNA and DNA as well as dinucleoside pyrophosphates; or any other methods described in the art. Fragmented DNA may be treated with one or more enzymes, for example, an endonuclease, prior to ligation of adaptors to one or both ends to facilitate ligation by generating ends that are compatible with ligation.
  • Adaptors may also incorporate modified nucleotides that modify the properties of the adaptor sequence. For example, phosphorothioate groups may be incorporated in one of the adaptor strands. A phosphorothioate group is a modified phosphate group with one of the oxygen atoms replaced by a sulfur atom. In a phosphorothioated oligo (often called an “S-Oligo”), some or all of the internucleotide phosphate groups are replaced by phosphorothioate groups. The modified backbone of an S-Oligo is resistant to the action of most exonucleases and endonucleases. Phosphorothioates may be incorporated between all residues of an adaptor strand, or at specified locations within a sequence. A useful option is to sulfurize only the last few residues at each end of the oligo. This results in an oligo that is resistant to exonucleases, but has a natural DNA center.
  • The term “array” as used herein refers to an intentionally created collection of molecules which can be prepared either synthetically or biosynthetically. The molecules in the array can be identical or different from each other. The array can assume a variety of formats, for example, libraries of soluble molecules; libraries of compounds tethered to resin beads, silica chips, or other solid supports.
  • The term “epigenetic” as used herein refers to factors other than the primary sequence of the genome that affect the development or function of an organism, they can affect the phenotype of an organism without changing the genotype. Epigenetic factors include modifications in gene expression that are controlled by heritable but potentially reversible changes in DNA methylation and chromatin structure. Methylation patterns are known to correlate with gene expression and in general highly methylated sequences are poorly expressed.
  • The term “genome” as used herein is all the genetic material in the chromosomes of an organism. DNA derived from the genetic material in the chromosomes of a particular organism is genomic DNA. A genomic library is a collection of clones made from a set of randomly generated overlapping DNA fragments representing the entire genome of an organism.
  • Restriction enzymes or restriction endonucleases and their properties are well known in the art. A wide variety of restriction enzymes are commercially available, from, for example, New England Biolabs. Restriction enzymes recognize a sequence specific sites (recognition site) in DNA. Typically the recognition site varies from enzyme to enzyme and may also vary in length. Isoschizomers are enzymes that share the same recognition site. Restriction enzymes may cleave close to or within their recognition site or outside of the recognition site. Often the recognition site is symmetric because the enzyme binds the double stranded DNA as homodimers.
  • Recognition sequences may be continuous or may be discontinuous, for example, two half sites separated by a variable region. Cleavage can generate blunt ends or short single stranded overhangs.
  • In preferred aspects of the present invention enzymes that include at least one CpG dinucleotide in the recognition site may be used. Enzymes with a recognition site that includes the sequence CCGG include, for example, Msp I, Hpa II, Age I, Xma I, Sma I, NgoM IV, Nae I, and BspE I. Enzymes with a recognition site that includes the sequence CGCG include, for example, BstU I, Mlu I, Sac II, BssH II and Nru I. Enzymes with a recognition site that includes the sequence GCGC include, for example, Hin P1 I, Hha I, Afe I, Kas I, Nar I, Sfo I, Bbe I, and Fsp I. Enzymes with a recognition site that includes the sequence TCGA include, for example, Taq I, Cla I, BspD I, PaeR7 I, Tli I, Xho I, Sal I, and BstB I. For additional enzymes that contain CpG in the recognition sequence. See, for example, the New England Biolabs catalog and web site. In some aspects two restriction enzymes may have a different recognition sequence but generate identical overhangs or compatible cohesive ends. For example, the overhangs generated by cleavage with Hpa II or Msp I can be ligated to the overhang generated by cleavage with Taq I. Some restriction enzymes that include CpG in the recognition site are unable to cleave if the site is methylated, these are methylation sensitive. Other enzymes that contain CpG in their recognition site can cleave regardless of the presence of methylation, these are methylation insensitive. Examples of methylation insensitive enzymes, that include a CpG in the recognition site, include BsaW I (WCCGGW), BsoB I, BssS I, Msp I, and Taq I. Examples of methylation sensitive enzymes, that include a CpG in the recognition site, include Aat II, Aci I, Acl I, Afe I, Age I, Asc I, Ava I, BmgB I, BsaA I, BsaH I, BspD I, Eag I, Fse I, Fau I, Hpa II, HinP1 I, Nar I, and SnaB I.
  • The terms “solid support”, “support”, and “substrate” as used herein are used interchangeably and refer to a material or group of materials having a rigid or semi-rigid surface or surfaces. In many embodiments, at least one surface of the solid support will be substantially flat, although in some embodiments it may be desirable to physically separate synthesis regions for different compounds with, for example, wells, raised regions, pins, etched trenches, or the like. According to other embodiments, the solid support(s) will take the form of beads, resins, gels, microspheres, or other geometric configurations. See U.S. Pat. No. 5,744,305 for exemplary substrates.
  • “Specific binding” refers to the ability of two molecular species concurrently present in a heterogeneous (inhomogeneous) sample to bind to one another in preference to binding to other molecular species in the sample. Typically, a specific binding interaction will discriminate over adventitious binding interactions in the reaction by at least two-fold, more typically by at least 10-fold, often at least 100-fold. Typically, the affinity or avidity of a specific binding reaction is least about 107 M−1, using at least 108 M−1 to at least about 109 M−1, and often greater, including affinities or avidities up to 1010 M−1 to 1012 M−1.
  • Methods for Multiplexed Capture of Genomic DNA
  • Many of the embodiments described below for analysis of methylation employ methods of capturing a population of pre-selected target sequences from a genomic DNA sample. In some aspects the methods rely on dU probe technology as previously disclosed in US Patent Application No. 20030096291 and U.S. Pat. No. 7,208,295. Additional methods for making and using dU probes are also disclosed in 60/887,546 filed Jan. 31, 2007. Briefly, dU probes may be generated by performing PCR using a pair of primers for each target sequence. Each primer contains a target specific region in the 3′ portion and a 5′ common sequence. The PCR is performed in the presence of dUTP so uracil is incorporated into the amplification product. The amplicons from different reactions can be pooled and amplified as a single reaction using primers to the common regions in the primers. This results in a pool of dU probes for a collection of target sequences. The pool need only be generated once and aliquots can be used for amplification of the targets from different samples. To use the dU probes an aliquot is mixed with the target sample which has preferably been fragmented to generate ends that are compatible with the target specific region of the dU probes. Sequences complementary to the common regions are also added. The dU probe acts like a splint and the target sequence from the sample is ligated to the common regions and can then be amplified using common primers. These methods allow for amplification of a limited number of specific target sequences from a complex background, for example, 100 to 2,000 different exons of interest from genes of interest or promoter regions of interest can be amplified from human genomic DNA. The sequence of the dU probes determines what sequences will be amplified and variation in the target, for example, polymorphisms that are present in the sample but were not present in the nucleic acid used to generate the dU probes, still allow for amplification of the sequence in the target.
  • DNA captured by these methods can be used for a variety of purposes, including, for example, methylation analysis, genotyping analysis, resequencing analysis, copy number analysis, haplotype analysis, and variant detection. The captured DNA is ligated to common priming sequences at the 5′ and 3′ ends of the captured genomic DNA to facilitate amplification. The captured DNA can be treated prior to amplification with the common primers, for example, the DNA can be bisulfite modified to preserve an indication of methylation status. Other treatments are also possible, for example, the captured DNA may be subjected to an affinity separation prior to amplification. For example, 5 meC containing captured DNA may be isolated using an antibody to 5 meC and one or both fractions may be subsequently amplified.
  • In preferred aspects, a collection of template probes corresponding to a collection of targets of interest are used as template to facilitate the ligation of common priming sequences to the ends of the target sequences in a nucleic acid sample to be analyzed. After the ligation, the template probes are digested or separated so they are not targets for subsequence amplification. The target sequences can then be amplified from the nucleic acid sample being analyzed and the amplification product can be interrogated.
  • Each template probe is complementary to a different target of interest flanked by a priming sequence at the 5′ end and a second priming sequence at the 3′ end. The targets are allowed to hybridize to the corresponding template probe and common priming sequences are ligated to the ends of the target in a subsequent step. The template probes are removed, for example, by digestion and the targets are amplified using primers to the common priming sequences. Non-targets do not have the common priming sequences so they are not amplified. This allows for multiplex amplification of a large number of target sequences, for example, 200 to more than 20,000 selected target sequences from a complex sample such as a genome. The length of each target may be, for example, about 100 to 1000, about 200 to 1000, about 200 to 500, about 200-2000 or about 100 to 5,000 bases. In one embodiment the ends of the targets may be defined by restriction sites in the genomic DNA sequence.
  • The template probes may also be used to mark the targets with one or more additional sequences. In a particularly preferred embodiment each template probe has a unique barcode sequence between one of the priming sequences and the target complementary region. The target is hybridized to the template probe and oligonucleotides that are complementary to the tags and to the common priming sequences are added and allowed to hybridize to the template. The pieces are hybridized to the template probe so that the ends are juxtaposed and can be ligated to form a contiguous sequence. In some aspects template probes have more than one barcode sequence. The barcode sequence can be used as a unique identifier of subsequent products.
  • The template probes are synthesized so that the target complementary region has defined ends. The genomic DNA sample from which targets are to be amplified is treated so that the targets have defined ends that correspond to the template probes. This can be accomplished in a variety of ways, for example, the ends may be generated by restriction enzyme digestion or by PCR.
  • One method for obtaining target sequences with defined ends is illustrated in FIG. 1A, in which a selected target sequence 105 with defined end bases at 101 and 103 is amplified from a larger sequence using locus specific primer extension. The resulting template 120 has ends that are defined by the sequence of oligonucleotide primers 109 and 113. The 5′ end of 109 corresponds to position 103 and the 5′ end of primer 113 corresponds to position 101. In some embodiments the starting DNA is bisulfite modified genomic DNA.
  • In the first step of the method, illustrated in FIG. 1A, primer 109 is mixed with the target sequence and extended to form extension product 107. Primer 113 is hybridized to extension product 107 and extended to generate extension product 115 which has the sequence of target 105. A second copy of primer 109 can hybridize to extension product 115 to generate a double stranded fragment.
  • In FIG. 1B single stranded fragment 115 Probe 130 includes at least first target complementarity region 120 and oligonucleotide positioning regions 119 and 121 which are directly adjacent to 120. In the annealing step, template region 113 hybridizes to template complementarity region 120 of probe 130 and oligonucleotide regions 15 and 17 concurrently hybridize to oligonucleotide positioning regions 119 and 121 of probe 130.
  • The nucleotide of template complementarity region 120 and the nucleotide of oligonucleotide positioning regions 119 and 121 that are directly adjacent within probe 130 define junctions within probe 130, and may be referred to as junctional nucleotides.
  • In preferred aspects the template probes 130 are synthesized in the presence of dUTP so that U is incorporated into the probe in at least one position and preferably multiple positions. The resulting “dU probes” can be degraded by UDG treatment after they are used as template probes.
  • In an alternative embodiment to the dU probe, the template probes need not contain uracil and are degraded by exonuclease treatment instead of by UDG treatment. This embodiment is illustrated in FIG. 1C. In this embodiment the 5′ common priming sequence [121] that is ligated to the 5′ end of the genomic target DNA sequence [119] is resistant to 5′ exonuclease activity so the template probe [130] is digested by the exonuclease but the ligated target sequence is not. The targeting region [120] of the probe and the common priming regions [125] of the probe need not have dU incorporated. The sequence may be made exonuclease resistant for example, by inclusion of 4 or 5 phosphorothioate linkages at the 5′ end of the oligonucleotide. This has been shown to render oligonucleotides resistant to treatment by T7 exonuclease (Nikiforov, PCR Methods and Applications, 3:285-291, 1994). In another aspect, the 3′ oligonucleotide [122] may be resistant to a 3′ to 5′ exonuclease activity and the template probe [130] may be digested with a 3′ to 5′ exonuclease activity after capture.
  • In another aspect STRANDASE λ Exonuclease, available from NOVAGEN is used to digest the template probe [130]. This enzyme digests one strand of a DNA duplex from a 5′ phosphorylated end. The enzyme can be used to prepare single stranded DNA from PCR products synthesized using one phosphorylated primer and one non-phosphorylated primer. If the 5′ end of the template probe [130] is phosphorylated and the 5′ end of oligonucleotide 121 is non-phosphorylated STRANDASE λ Exonuclease may be used to degrade the template probe after ligation and prior to amplification.
  • Methods for Detection of Methylation
  • Mammalian methylation patterns are complex and change during development, see van Steensel and Henikoff BioTechniques 35: 346-357 (2003).
  • Methylation in promoter regions is generally accompanied by gene silencing and loss of methylation or loss of the proteins that bind to the methylated CpG can lead to diseases in humans, for example, Immunodeficiency Craniofacial Syndrome and Rett Syndrome, Bestor (2000) Hum. Mol. Genet. 9:2395-2402. DNA methylation may be gene-specific and occurs genome-wide.
  • Methods for detecting methylation status have been described in, for example U.S. Pat. Nos. 6,214,556, 5,786,146, 6,017,704, 6,265,171, 6,200,756, 6,251,594, 5,912,147, 6,331,393, 6,605,432, and 6,300,071 and US Patent Application publication Nos. 20030148327, 20030148326, 20030143606, 20030082609, 20050153347 and 20050009059, each of which are incorporated herein by reference. Other array based methods of methylation analysis are disclosed in U.S. Patent Publication No. 20050196792 and U.S. patent application Ser. No. 11/213,273 filed Aug. 26, 2005. For a review of selected methylation detection methods, see, Oakeley, E. J., Pharmacology & Therapeutics 84:389-400 (1999). Available methods include, but are not limited to: reverse-phase HPLC, thin-layer chromatography, SssI methyltransferases with incorporation of labeled methyl groups, the chloracetaldehyde reaction, differentially sensitive restriction enzymes, hydrazine or permanganate treatment (m5C is cleaved by permanganate treatment but not by hydrazine treatment), sodium bisulfite, combined bisulphate-restriction analysis, and methylation sensitive single nucleotide primer extension.
  • Other related methods of methylation analysis and arrays that may be useful in conjunction with the disclosed methods include those disclosed in U.S. Provisional Application Nos. 60/744,705 and 60/788,520.
  • In some embodiments the methods include treatment of the sample with bisulfite. Unmethylated cytosine is converted to uracil through a three-step process during sodium bisulfite modification. The steps are sulphonation to convert cytosine to cytosine sulphonate, deamination to convert cytosine sulphonate to uracil sulphonate and alkali desulphonation to convert uracil sulphonate to uracil. Conversion on methylated cytosine is much slower and is not observed at significant levels in a 4-16 hour reaction. See Clark et al, Nucleic Acids Res., 22(15):2990-7 (1994). If the cytosine is methylated it will remain a cytosine (see FIG. 4). If the cytosine is unmethylated it will be converted to uracil. When the modified strand is copied, through, for example, extension of a locus specific primer, a random or degenerate primer or a primer to an adaptor, a G will be incorporated in the interrogation position (opposite the C being interrogated) if the C was methylated and an A will be incorporated in the interrogation position if the C was unmethylated. When the double stranded extension product is amplified those Cs that were converted to U's and resulted in incorporation of A in the extended primer will be replaced by Ts during amplification. Those Cs that were not modified and resulted in the incorporation of G will remain as C. Bisulfite treatment can degrade the DNA making it difficult to amplify. Also the sequence degeneracy resulting from the treatment complicates primer design. The treatment may also result in incomplete desulfonation, depurination and other as yet uncharacterized DNA damage, making downstream processing more challenging. The treatment can also result in preferential amplification of unmethylated DNA relative to methylated DNA. This may be mitigated by increasing the PCR extension time.
  • Kits for DNA bisulfite modification are commercially available from, for example, Human Genetic Signatures' Methyleasy and Chemicon's CpGenome Modification Kit. See also, WO04096825A1, which describes bisulfite modification methods and Olek et al. Nuc. Acids Res. 24:5064-6 (1994), which discloses methods of performing bisulfite treatment and subsequent amplification on material embedded in agarose beads. In one aspect a catalyst such as diethylenetriamine may be used in conjunction with bisulfite treatment, see Komiyama and Oshima, Tetrahedron Letters 35:8185-8188 (1994). Diethylenetriamine has been shown to catalyze bisulfite ion-induced deamination of 2′-deoxycytidine to 2′-deoxyuridine at pH 5 efficiently. Other catalysts include ammonia, ethylene-diamine, 3,3′-diaminodipropylamine, and spermine. In some aspects deamination is performed using sodium bisulfite solutions of 3-5 M with an incubation period of 12-16 hours at about 50° C. A faster procedure has also been reported using 9-10 M bisulfite pH 5.4 for about 10 minutes at 90° C., see Hayatsu et al, Proc. Jpn. Acad. Ser. B 80:189-194 (2004).
  • Bisulfite treatment allows the methylation status of cytosines to be detected by a variety of methods. For example, any method that may be used to detect a SNP may be used, for examples, see Syvanen, Nature Rev. Gen. 2:930-942 (2001). Methods such as single base extension (SBE) may be used or hybridization of sequence specific probes similar to allele specific hybridization methods. In another aspect the Molecular Inversion Probe (MIP) assay may be used.
  • In a preferred aspect, molecular inversion probes, described in Hardenbol et al. Genome Res. 15:269-275 (2005) and in U.S. Pat. No. 6,858,412, may be used to determine methylation status after methylation dependent modification. A MIP may be designed for each cytosine to be interrogated. In a preferred aspect the MIP includes a locus specific region that hybridizes upstream and one that hybridizes downstream of an interrogation site and can be extended through the interrogation site, incorporating a base that is complementary to the interrogation position. The interrogation position may be the cytosine of interest after bisulfite modification and amplification of the region and the detection can be similar to detection of a polymorphism. Separate reactions may be performed for each NTP so extension only takes place in the reaction containing the base corresponding to the interrogation base or the different products may be differentially labeled.
  • In one embodiment methods for identifying methylation based on differential cleavage by restriction enzymes are used. Methylation-sensitive restriction analysis followed by PCR amplification or Southern analysis have been disclosed, for example, in Huang et al. Cancer Res. 57:1030-1034 (1997), Zuccotti et al, Methods in Enzym. 225:557-567 (1993), Carrel et al. Am Jour. Med. Genet. 64:27-30 (1996) and Chang et al. Plant MoL Bio. Rep. 10:362-366 (1992).
  • In some embodiments of the present invention, at least one oligonucleotide is appended to at least 2 templates of distinct sequence, typically at least 5 templates of distinct sequence, even at least 10, 20, 30, 40, or even at least 50 templates of distinct sequence, and may be appended to 100, 500, 1000, even 5000 or more templates of distinct sequence.
  • Barcode tags are short nucleic acids having sequence that is designed algorithmically to maximize discrimination on a microarray displaying complements of the respective tags; a 1:1 correspondence as between tag sequence and nucleic acid to which it is appended permits each such nucleic acid to be identified by detection of the bar code uniquely associated therewith. See, e.g. Shoemaker et al. Nature Genet. 14(4):450-6 (1996); EP 0799897; Fan et al. Genome Res. 10:853-60 (2000); and U.S. Pat. No. 6,150,516, the disclosures of which are incorporated herein by reference in their entireties.
  • In the methods of the present invention, a distinct bar code sequence may be included in each species. In these embodiments, the terminal region of each species of oligonucleotide is distinct in sequence, and can anneal only to a single species of probe. The 1:1 correspondence as between tag sequence and template-appended oligonucleotide thus permits each template or template amplification product to be identified by detection of the barcode uniquely associated therewith. Appending common first and second priming sites directly to each of the plurality of templates of distinct sequence—without prior amplification of the template—facilitates the subsequent stoichiometric amplification and analysis of a wide variety of templates of distinct sequence, obviating the problems of unequal amplification observed with many multiplex PCR approaches. By permitting the de novo design of the priming sites, independent of considerations of template sequence, the methods of the present invention also permit amplification with primers having optimal hybridization characteristics, decreasing artifacts such as primer dimer formation.
  • In preferred aspects of the present methods enzymes that include at least one CpG dinucleotide in the recognition site may be used. Enzymes with a recognition site that includes the sequence CCGG include, for example, Msp I, Hpa II, Age I, Xma I, Sma I, NgoM IV, Nae I, and BspE I. Enzymes with a recognition site that includes the sequence CGCG include, for example, BstU I (CGCG, MSRE), Mlu I (ACGCGT, MSRE), Sac II (CCGCGG, MSRE), BssH II (GCGCGC, MSRE) and Nru I (TCGCGA, MSRE). Not I, BstZ I, Csp I and Eag I have two CpG's in their recognition sites and cleavage is blocked by CpG methylation. Enzymes with a recognition site that includes the sequence GCGC include, for example, Hin P1 I, Hha I, Afe I, Kas I, Nar I, Sfo I, Bbe I, and Fsp I. Enzymes with a recognition site that includes the sequence TCGA include, for example, Taq I, Cla I (MSRE), BspD I (MSRE), PaeR7 I, Tli I, Xho I, Sal I, and BstB I. For additional enzymes that contain CpG in the recognition sequence and for information about the enzyme's sensitivity to methylation see, for example, the New England Biolabs catalog and web site. In some aspects two restriction enzymes may have a different recognition sequence but generate identical overhangs or compatible cohesive ends. For example, the overhangs generated by cleavage with Hpa II or Msp I can be ligated to the overhang generated by cleavage with Taq I. Some restriction enzymes that include CpG in the recognition site are unable to cleave if the site is methylated, these are methylation sensitive restriction enzymes (MSRE). Other enzymes that contain CpG in their recognition site can cleave regardless of the presence of methylation, these are methylation insensitive restriction enzymes (MIRE). A third type of enzyme cleaves only when the recognition site is methylated, and are referred to herein as methylation dependent restriction enzymes (MDRE). Examples of MIREss that have a CpG in the recognition sequence include, for example, BsaW I (WCCGGW), BsoB I, BssS I, Msp I, and Taq I. Examples of MSREs, that include a CpG in the recognition site, include Aat II, Aci I, Acl I, Afe I, Age I, Asc I, Ava I, BmgB I, BsaA I, BsaH I, BspD I, Cla I, Eag I, Fse I, Fau I, Hae III, Hpa II, HinP1 I, Mlu I, Nar I, Not I, Nru I, Pvu I, Sac II, Sal I, Sma I and SnaB I. In preferred aspects a pair of enzymes that have differential sensitivity to methylation and cleave at the same recognition sequence. with one member of the pair being a MSRE and the other member being a MIRE is used. Such pairs include, for example, Hpa II (MSRE) and Msp I (MIRE). Both have recognition sequence CCGG. Another pair that may be used is Sma I (MSRE) and Xma I (MIRE). In another aspect the pair is Bis I a MDRE and Fnu4HI a MSRE. Both enzymes recognize the site GCNGC but Bis I cleaves specifically at methylated C within the sequence G(5mC) ˆNGC [Degtyarev et al., Izv. Sib. Otd. Akad. Nauk SSSR 15:25-26 (1989) and Chmuzh et al. Biotekhnologia 3:22-26 (2005)], (ˆ indicates the position of cleavage within the recognition sequence) while cleavage by Fnu4HI at GCˆNGC is blocked by methylation. BthCI also cuts at GCNGˆC. GlaI is another MDRE recognizing G(5m)CˆGC, see Chernukhin et al., SibEnzyme Scientific Library, October 2005. MSREs to pair with GlaI may include, for example, HpaI (GCGˆC) and HinP1I (GˆCGC). DpnI is another MDRE recognizing G(m6A)ˆTC while MboI (ˆGATC), ChaI (GˆTCA) and BstKTI (GATˆC) recognize unmethylated GATC sites.
  • FIG. 2A and FIG. 2B illustrate the use of the dU probe for multiplex amplification. FIG. 2A shows the components of the dU probe [200]. The dU probe has a central targeting region [201] composed of dU, dG, dC and dA, a tag region [203] and two common priming regions, one at the 5′ end of the probe [205] and the second at the 3′ end of the probe [207]. There is a dU probe for each target to be amplified. Each probe in the collection of dU probes preferably has the same common regions as other probes in the collection, facilitating universal priming for amplification by, for example, PCR. The dU probes differ in the sequence of the targeting region and the tag region, although groups of dU probes may share the same tag region in some aspects.
  • In FIG. 2B, the dU probes [200] are mixed with the sample [209] from which the targets [211] are to be amplified along with sequences that are complementary to the tag region [213] and to the common regions [215, 217]. The target, tag complements, and common region complements hybridize to the dU probe to form a complex [219]. A ligase is added to ligate the ends of the annealed fragments to form a ligation product [221]. UDG is added to cleave the dU probes and the remaining ligated target [221] containing fragment is amplified. The presence of the fragment may then be detected, for example, by hybridization to an array of probes. In some aspects the array of probes may include probes that are complementary to the tag sequence. In another aspect the array is a genotyping array. In another aspect the array is a resequencing array. In another aspect the array is a copy number array. In another aspect the array is a gene expression array.
  • Detection of bisulfite converted DNA on arrays is complicated for several reasons. First, the converted sequences are generally AT rich and often have low complexity. Second, because there may be multiple CpGs in a single probe target and each may be methylated or unmethylated there is a combinatorial problem with the number of different possible perfect match probes that may be necessary. For example, if the starting sequence is GACGAACGCGGCT (SEQ ID NO. 9), there are three cytosines that may be methylated or unmethylated. This generates 8 possible outcomes after bisulfite treatment:
    GACGAACGCGGCT (SEQ ID NO.9)
    GACGAACGTGGCT (SEQ ID NO.10)
    GACGAATGTGGCT (SEQ ID NO.11)
    GATGAATGTGGCT (SEQ ID NO.12)
    GATGAATGCGGCT (SEQ ID NO.13)
    GATGAACGCGGCT (SEQ ID NO.14)
    GACGAATGCGGCT (SEQ ID NO.15)
    GATGAACGTGGCT (SEQ ID NO.16)

    Each different possible outcome could be targeted by a different perfectly complementary probe. For the 8 possible outcomes there would be 8 different perfect match probes on the array. The probe could be complementary to either strand.
  • A pilot experiment was performed using the following sequences as target loci (CGs are underlined):
    (SEQ ID NO.17)
    AGCTGGTGATGCTGATCAGAGCCTCTGTAGTCTTAAATGACTTTTCTAAC
    TAATTCTAAATCTTCAGAACCCATCGTATAAAAAGGCCATACCTTCTGGA
    GGGACGTCGATGGTATTAGGATAGAAGCACCAGGGGACCCCACGAACGGT
    GTCGTCGAAACAGCAGCCCTTATTTGCACACTGGGAGGG
    (SEQ ID NO.18)
    GGGACCACCCTTATAAGGCTCGGAGGCCGCGAGGCCTTCGCTGGAGTTTC
    GCCGCCGCAGTCTTCGCCACCAGTGAGTACGCGCGGCCCGCGTCCCCGGG
    GATGGGGCTCAGAGCTCCCAGCATGGG
    (SEQ ID NO.19)
    GCAGGGGAGGGAAGCAGATGCCAGCGGGCCGAAGAGTCGGGAGCCGGAGC
    CGGGAGAGCGAAAGGAGAGGGGACCTGGCGGGGCACTTAGGAGCCAACCG
    AGGAGCAGGAGCACGGACTCCCACTGTGGAAAGGAGGACCAGAA
  • An array was designed to interrogate both forward and reverse strands of the three target loci. For each CpG 21, 23 and 25 base probes were included. For the 21 mers and 23 mers the probes were tiled from −1 to +5 with the central CpG at the central or 0 position. For the 25 mer probes they were tiled from −3 to +3 with the central CpG at 0. If multiple CpGs were present in a probe all combinatorial possibilities were represented on the array. All probes were included with 9 fold redundancy. Completely unmethylated or completely methylated DNA was bisulfite converted (Zymo kit), subjected to locus specific amplification, labeled and hybridized to the array (2 pM or 20 pM). The hybridization conditions were 3M TMAC, 10 mM Tris pH 7.8, 0.01% Tween-20, 0.5 mg/ml BSA, 0.1 mg/ml HS DNA at 49° C. overnight. Low stringency wash conditions were 6×SSPE and 0.1% Tween-20. High stringency wash conditions were 0.6×SSPE and 0.1% Tween 20.
  • For data analysis only the probes that detect either completely methylated (“C probe signal”) or completely unmethylated (“T probe signal”) were used since the starting DNA was completely methylated or unmethylated. Two values were plotted: the log signal sum=log (C probe signal+T probe signal) and the signal contrast=(T probe signal-C probe signal)/(C probe signal+T probe signal). For methylated the signal contrast is expect to be −1 (all signals from C probes) and for unmethylated the signal contrast is expected to be 1 (all signal from T probes). Plotting “log sig sum” versus “sig contrast” allows clean distinction for most signal points. Analysis of the data to determine effects of probe length shows no dramatic distinction between 21, 23 and 25 mer probes, although there is a cluster of 21 mer probes around 0 which indicates that the C and T probes are behaving similarly and not distinguishing well between methylated and unmethylated.
  • Mixtures of methylated and unmethylated DNA at varying percentages were tested including 0, 25, 50, 75 and 100% methylated. The points cluster in the expected order, with 100% being to the left of 75% which is to the left of 50% which is to the left of 25% which is to the left of 0%, but 75, 50 and 25 are shifted slightly to the left of expected, with 25% clustering around a sig contrast of 0, 50 at a sig contrast of about −0.25 and 75 at a sig contrast of about −0.75.
  • To determine how well the bisulfite read out works in the context of a more complex background (1) unmethylated genomic DNA (2) methylated genomic DNA or (3) unmethylated 2000-plex DNA generated from dU capture probes was treated with bisulfite, (1) and (2) were subjected to locus specific amplification, (3) was amplified using the common primers from the dU probes. The amplification products were either pooled prior to fragmentation and hybridization or fragmented and hybridized without pooling. The results indicate that even in a complex background of bisulfite converted unmethylated DNA the methylation states of the spike ins are clearly distinguishable.
  • FIG. 3 illustrates an embodiment in which priming sites containing methylated cytosines are added to genomic DNA using uracil containing target specific template probes 309. The uracil containing template probes each have a target specific region 313 and 314 that is flanked by common priming sites 311 and 312. The common priming sites may be the same or different. Oligonucleotide 315 is complementary to 311 and 312 and contains methylated cytosines that are resistant to bisulfite conversion. The oligonucleotide 315 and the genomic targets 305 and 307 hybridize to the dU probes to form duplexes 317. Ligation between the ends of targets 305 and 307 and oligonucleotides 315 results in product complexes [319] that include the targets with common primers ligated to either end in a complex with the complementary dU probe. The duplexes in 319 are subjected to UDG cleavage to digest 309 leaving singles started genomic DNA flanked by common priming sites [321]. The target sequences 305 and 307 hybridize to the target specific regions 313 and 314, respectively, of the template probes. The uracil containing probe is then digested with UDG to leave target genomic DNAs flanked by common priming sites [321] which are then treated with bisulfite to generate bisulfite modified genomic fragments flanked by unmodified common priming sites [323] which are then amplified using primers to the common sequences to obtain multiplex amplified targets [325] that can be analyzed to determine if the cytosines from the genomic DNA were methylated or unmethylated. Analysis may be, for example, by hybridization to an array of probes with sequences that are specific for the sequence resulting from presence or absence of methylation, but other methods such as sequence specific PCR may also be used. See, for example, U.S. Pat. Nos. 6,265,171 and 6,200,756.
  • FIG. 4 illustrates the alteration of DNA sequence resulting from treatment of DNA with bisulfite. SEQ ID NOs. 1 and 2 are the complementary starting strands. Methylated cytosines are indicated by “m”. SEQ ID NO. 3 corresponds to SEQ ID NO. 1 after treatment and SEQ ID NO. 4 corresponds to SEQ ID NO. 2 after treatment. Unmethylated cytosines have been converted to “U” and methylated cytosines remain “C”. After PCR amplification SEQ ID NO. 3 gives SEQ ID NOs. 5 and 6, while SEQ ID NO. 4 results in SEQ ID NOs. 7 and 8. Many of the G-C base pairs from the starting sequence have been converted to A-T base pairs; resulting in a relatively low GC content for the final sequence.
  • FIG. 5 illustrates a method for methylation analysis of a selected group of target sequences using dU probe technology. Genomic DNA [501] is digested with one or more restriction enzymes to generate restriction fragments [503]. The restriction fragments are mixed with a collection of dU probes for specific targets [504] and common primer sequences [505] and the nucleic acids are denatured and allowed to anneal. The annealing step results in the formation of complexes between the dU probes, target fragments and primers [506]. There is a ligation step and a cleanup to enrich for the complexes with ligated targets, common primers and dU probes [507] and then the complexes are treated to digest the dU probes, leaving the targets ligated with primer sequences at either end [509]. These are treated with bisulfite to generate bisulfite modified targets flanked by common priming sequences [511] which may then be amplified by PCR to generate amplification products of the selected targets [513]. In this embodiment the primer sites are subjected to modification by bisulfite prior to amplification so steps should be taken to maintain the sequence of the primers. In one aspect the primers that are ligated to the target sequences contain 5-methyl cytosine so they are not changed in sequence by the bisulfite treatment. The amplification product is enriched for the targets of the dU probes. The amplification product is analyzed by hybridization to an array [515].
  • FIG. 6 illustrates a method of determining the methylation state of a cytosine in a methylation sensitive restriction site. Fragment 601 has an unmethylated MspI/HpaII site while fragment 602 has a methylated site. MspI and HpaII are isoschizomers with MspI being methylation insensitive and HpaII being methylation sensitive. The DNA is cleaved with HpaII so that unmethylated sites are cleaved while methylated sites remain uncleaved. The fragments are then end filled using Klenow. The DNA is then cleaved with the MspI so that the methylated sites are cleaved. The resulting fragments have one strand 603 that has a variable end sequence, the variation is the difference between end filling and no in filling. In the example illustrated the unmethylated fragment results in 603 a with CCG at the 3′ end and the methylated fragment results in 603 b with C at the 3′ end. The next step is the addition of dU probes that are specific for each type of fragment. Probe 604 is specific for 603 a and probe 605 is specific for 603 b. The dU probes have different tag sequences 606 and 607. The probes are shown hybridized to both fragments. The dU probes have common priming sequences 608 and target complementary region 609. The tags are optional and allow for differential detection. In FIG. 6B sequences that are complementary to the common priming sequences 608 and the tag sequences 606 and 607 are hybridized to the dU probes. Probes 604 and 605 can hybridize to either the product generated from the methylated or unmethylated sites but ligation between the 603 a and 606 is blocked and ligation of 603 b and 607 is blocked.
  • Fragment 603 a ligates with 607 and 608 to generate target 610 and fragment 603 b ligates with 606 and 608 to generate target 611. The dU probes can then be digested to leave targets 610 and 611 which are tagged with either 606 or 607, depending on methylation, and are flanked by common priming sites 608. Targets 610 and 611 are amplified by PCR using a primer complementary to 608. The methylation state of the site can be determined by detecting which tag sequence or sequences are present in the amplified sample. If 607 is present the site was unmethylated and if 606 is present the site was methylated. The common priming sequences and the tag are ligated to the fragments in a methylation specific manner. The dU probe is digested and the products are amplified and the tags detected.
  • In FIG. 7A a fragment with an unmethylated restriction site [701] is shown on the left and a fragment where the site is methylated [703] is shown on the right. The fragments are digested with a methylation sensitive restriction enzyme so that the fragment on the left is cleaved but the fragment on the right is not. Klenow is added and the ends of the fragments on the left are end filled to add complementary sequences [705]. The products are then digested with a methylation insensitive enzyme so that the site in the fragment on the right is cleaved but the ends are not filled. The difference between strand 707 and strand 709 is that 707 has additional sequence 705—this is filled in and will be complementary to the fill-in specific probe while the sequence 709 is not filled in and will be targeted by the no-fill specific probe.
  • In FIG. 7B the different fragments 707 and 709 are differentially detected using dU probes. The dU probes include a probe [711] specific for 709 (no fill specific) and a probe [713] specific for 707 (fill-in specific). Probe 713 has an extra region [737] corresponding to the filled in region resulting from cleavage with the methylation sensitive enzyme. Probes 711 and 713 also vary in the sequence present at the tag complement regions, where 731 c is complementary to tag 731 and 733 c is complementary to tag 733. The dU probes can form complexes with the non-cognate targets to form complexes [715 and 713] but only the correctly formed complexes [719 and 721] provide the proper template for ligation of all the elements to form the complexes with the correctly ligated products [723 and 725]. After UDG treatment the ligated fragment products [727 and 729] are differentially labeled with different tags [731 and 733] and flanked by common primers [735] for amplification. The different tag sequences allow detection of the initial methylation event. If the site is methylated only 729 is detected (tag 731). If the site is unmethylated only 727 is detected (tag 733).
  • FIG. 7C shows a schematic of a dU probe capture experiment using dU probes for 9 different targets. The results are shown in the gel image on the right. The amount of dU probe added was varied in each reaction. Lane 1 is a standard ladder, lane 2 has no dU probe, lane 3 has 2 amol, lane 4 has 6, lane 5 has 18, lane 6 has 54, and lane 7 has 162. Two different dU probes are shown schematically. There is 0.1 amol of the target present.
  • FIG. 8 shows the results of 384-plex capture using a mixture of 384 target specific dU probe methodology. Stained agarose gel images are shown. The reactions run on the upper gel used either Buffer A (15 mM MgCl2, 33 mM Tris-HCl pH 7.5, 0.1 mg/ml BSA) or Buffer B (33 mM Tris-acetate pH 7.9, 10 mM MgCl2, 66 mM K-acetate, 0.1 mg/ml BSA) as indicated and the reactions run on the lower gel used either Buffer C (33 mM Tris-acetate pH 7.9, 10 mM MgCl2, 66 mM K-acetate, 0.1 mg/ml BSA, 0.2 U Primer Navigator) or Buffer D (6×SSPE, 7.5% PEG-8000) as indicated. The reactions on the left have ligase while the reactions on the right are without ligase. As expected, amplification products are only observed in the presence of ligase. Each of the different buffer conditions gave specific amplification products in the expected size range. Buffer A gives more amplification product than Buffer B at lower dU probe concentrations. The reactions have increasing amounts of dU probe as indicated graphically above each gel image. In many embodiments one or more reaction conditions may be varied. In some aspects the basic starting conditions are 6×SSPE with 7.5% PEG, the annealing conditions are 94° C. for 1 min, 70° C. for 40 min, 65° C. for 40 min, 60° C. for 40 min, 55° C. for 10 min, and 50° C. for 10 min. Other ranges of conditions that may be used include 0-15 mM MgCl2, 0-200 mM KCl, variation of the pH, addition of RecA and SSB, including the use of thermophilic SSB, and variation of the reaction volume.
  • FIG. 9 shows an image of a gel separating the products of amplification in the presence of fill-in specific probes (center) or no fill specific probes (right side of gel) or both (left side of gel). HpaII and MspI are the enzymes used. The genomic DNA is from liver.
  • In preferred aspects there are 200-1000, 1000-2500, 2000-5000, 5000-10,000, 10,000-20,000 or more than 20,000 different target sequences analyzed. There may be a dU probe for each target sequence to be amplified.
  • In another aspect, common priming sites are added to each of a plurality of templates of distinct sequence by first hybridizing locus specific primers to the targets and extending the locus specific primers.
  • Many methods of determining the genotype or methylation status of two positions that are near to one another do not provide information about chromosome. For example, if SNP1 and SNP2 are in the same gene and both are interrogated and found to be heterozygotes (A/B) there are 4 possibilities for the each copy of the gene-it could be A for SNP1 and A for SNP2, B for SNP1 and B for SNP2, A for SNP1 and B for SNP2 or B for SNP1 and A for SNP2. The genotypes of the individual SNPs don't provide information about the genotype of the neighboring SNPs unless there is additional information about linkage. Similarly, if the methylation state of a cytosine is determined by bisulfite modification based analysis it may be useful to determine if both C's are methylated on the same chromosome or on different chromosomes.
  • Analysis of determining “epihaplotypes” using inversion probes are shown in FIG. 10. In one embodiment molecular inversion probe methods are used to analyze the methylation status of neighboring cytosines. When determining the genotype of two neighboring or linked SNPs in a diploid organism, many methods do not provide information about what the haplotype is, for example, if SNP1 is heterozygous A/G and SNP2 is heterozygous C/T, there is no information about how the SNPs are combined on each of the chromosomes. It would be useful to know if the A allele of SNP 1 is linked to the C or T allele of SNP 2 or if the G allele of SNP 1 is linked to the C or T allele of SNP 2. Similarly, with methylation in a diploid it may be desirable to know if two neighboring cytosines are both methylated in one copy of the gene and not methylated in the other copy. In one aspect genomic DNA is treated with bisulfite to differentially modify methylated and unmethylated cytosines. Each MIP has an interrogation position corresponding to a first cytosine position (1019) and the second cytosine position (1020) is positioned between the target complementary arms (1016 and 1018) of the MIP. In the figure, the interrogation base (1017) of the MIP is shown as a G corresponding to methylation at 1019. One of skill in the art will recognize that a MIP probe could be designed with a T at the interrogation position which would circularize only if the position was unmethylated. Similarly, the MIP probe could be designed so that the interrogation position was complementary to position 1023 with or without methylation and position 1019 could be positioned in the gap between 1026 and 1018.
  • As shown in the figure, the G at 1017 hybridizes to the C at 1019 and can be extended to fill the gap 1020, including the complement of position 1023. In the chromosome (1022) where position 1019 was not methylated and was converted to a U by bisulfite, the G does not hybridize efficiently and is not extended to fill in the gap. After gap filling and ligation step (1030) the MIP on the left is circular while the MIP on the right is linear. Linear MIPs may be digested by exonuclease in step 1032 and the remaining circular probes can be linearized at the cleavage site 1004 and amplified using primers complementary to the primer sites 1002 and 1006. The amplified fragments 1036 can be interrogated for both the presence of tag sequences 1010 and for the sequence present at the second cytosine position. The presence of the tag indicates the methylation status of position 1019, if tag 1010 is present position 1019 was methylated. The methylation status of position 1023 is interrogated by allele specific hybridization. MIP probes for different methylation states of position 1019 have different tags and are separated into separate reactions.
  • The array (1037) has probe 1039 to tag 1010 and probe 1041 and 1043 that are complementary to portions of 1016 and 1018 and to gap 1020 with 1041 being perfectly complementary to the gap generated if position 1023 was not methylated and probe 1043 being perfectly complementary to the gap generated if position 1023 was methylated. This method may also be used to detect the genotypes of two neighboring SNPs.
  • FIG. 11 shows a schematic of a method of forming a closed circular probe from a linear MIP. For additional details on features of molecular inversion probes and methods of using molecular inversion probes see U.S. Pat. No. 6,858,412 which is incorporated herein in its entirety for disclosure of methods related to MIP. The MIP has sequences 1116 and 1118 that are complementary to template 1100 and hybridize to 1100 leaving a gap 1120. The gap may be 1 or more bases or it may be simply a break in the phosphate backbone that may be closed by ligation. The MIP also has barcode sequence 1110, priming sites 1102 and 1106, first cleavage site 1104 and optional spacer sequences 1108 and 1114 which can contain other functional sequences such as priming sites, restriction sites, promoter sites, additional barcode sequences or detectable labels.
  • An optional extension of either 1116 or 1118 (depending on which has a free 3′ end) may be used to fill in the gap at 1120 followed by ligation of the free ends of the linear probe at 1124. Linear probes may be digested by exonuclease treatment. Circular probes may then be linearized, for example by cleavage at 1104. In one aspect the cleavage site 1104 comprises one or more uracil bases. The probe may be treated after the ligation step with uracil-DNA glycosylase (“UDG”), which catalyzes the release of free uracil from uracil-containing DNA, creating apurinic (“AP”) sites.
  • AP sites may then be cleaved enzymatically using an AP endonuclease or, under certain conditions, an AP lyase. For example, the AP site may be cleaved using ApeI, an Endo IV or Fpg (formamidopyrimidine [fapy]-DNA glycosylase; also known as 8-oxoguanine DNA glycosylase). Fpg cleaves both 3′ and 5′ to the AP site, removing the AP site and leaving a 1 base gap. Alternatively, the AP sites may be cleaved chemically, such as by treatment with 1,4 diaminobutane and heat. AP sites may also be cleaved by exposure to high pH.
  • FIG. 12 shows a MIP before and after probe inversion. Before inversion the ends are 1201 and 1213. After circularization and cleavage at 1205 (probe inversion) the ends are 1207 and 1203. Primers for 1207 and 1203 are used to amplify the inverted probe. The amplification products are typically double stranded while the probe is single stranded. The amplification products can be cleaved at 1211 to separate 1207 and 1209 from the remainder of the probe. The tag sequence [1209] may be detected by hybridization to an array [1217] of tag probes [1219]. A label [1221] may be used to detect a hybridized tag sequence.
  • In some aspects the probes can include purines at a cleavage site such as 8-oxoguanine, 8-oxoadenine, fapy-guanine, methyl-fapy-guanine, fapy-adenine, aflatoxin B-fapy-guanine, 5-hydroxy-cytosine, and 5-hydroxy-uracil, that mimic damaged purines. Fpg glycosylase will release these residues from DNA and remove the resulting AP site, leaving a 1 nucleotide gap.
  • The probe may be amplified after circularization. Circularized probes that have been linearized may be amplified, for example by PCR using primers to regions 502 and 506. Other methods of amplification may also be used. For example, rolling circle amplification may be used to generate multiple copies of the circularized probes. See, for example, U.S. Pat. Nos. 5,648,245 and 5,854,033 and Fire and Xu, Proc Natl Acad Sci USA 92:4641-5 (1995). As used herein, the term amplification includes the production of RNA transcripts by polymerization driven from a phage promoter. For example, a T7 RNA polymerase promoter sequence may be incorporated into the probe and used to generate multiple RNA copies using a T7 RNA polymerase.
  • In preferred aspects, the amplification product is DNA produced by polymerization primed using one or more oligonucleotides (“primers”) that are capable of hybridizing to one or more priming sites within one or more of the oligonucleotides appended to the template. For example, a first primer capable of binding to a first priming site present in the first oligonucleotide may be used to prime unidirectional amplification. A second primer capable of binding to the complement of the second priming site present in the second oligonucleotide may be used concurrently to prime bidirectional amplification. In embodiments in which first and second priming sites are reverse complements of one another, the first and second primers may be the same. Amplification may be isothermal or thermal cycling.
  • Nucleic acid amplification methods useful in the methods of the present invention are well known in the art and include, e.g., polymerase chain reaction (PCR), nucleic acid sequence-based amplification (NASBA), self-sustained sequence recognition (3 SR), ligase chain reaction (LCR), transcription-mediated amplification (TMA), rolling circle amplification (RCA), and strand displacement amplification (SDA).
  • In another aspect the enzyme CIRCLIGASE (Epicenter) which ligates single stranded DNAs into circles may be used. DNA fragments may be denatured after digestion and CIRCLIGASE may be used to circularize the fragments after bisulfite treatment. The circles may then be amplified, for example by rolling circle amplification. The enzyme is thermostable and ATP-dependent, but will ligate ends of ssDNA in the absence of a complementary sequence without the requirement that the ends be annealed adjacently on one another.
  • In another aspect one or more of the methods disclosed herein is combined with the use of affinity bases pull-downs of 5mC containing DNA. Antibodies that recognize 5mC are commercially available. For example, ABCAM sells a 5-methyl cytosine antibody (ab1884).
  • All patents, patent publications, and other published references mentioned herein are hereby incorporated by reference in their entireties as if each had been individually and specifically incorporated by reference herein. While preferred illustrative embodiments of the present invention are described, one skilled in the art will appreciate that the present invention may be practiced by other than the described embodiments, which are presented for purposes of illustration only and not by way of limitation. The present invention is limited only by the claims that follow.

Claims (25)

1. A method for identifying sites of methylation in a genomic DNA sample said method comprising:
(a) amplifying a plurality of target sequences by a method comprising:
(i) hybridizing a plurality of first locus specific primers to the genomic DNA sample and extending said first locus specific primers,
(ii) hybridizing a plurality of second locus specific primers to the product of (i) and extending said second locus specific primers to obtain a plurality of templates comprising 5′ and 3′ ends defined by the first and second locus specific primers,
(iii) hybridizing a plurality of uracil containing probes to the product of (ii), wherein the uracil containing probes comprise a 5′ first common sequence, a 3′ second common sequence and a template complementary sequence, a first oligonucleotide that is complementary to the 5′ first common sequence and a second oligonucleotide that is complementary to the 3′ second common sequence, wherein the first and second oligonucleotides comprise methyl cytosine in place of cytosine,
(iv) ligating the first oligonucleotide to the 3′ end of the template and ligating the second oligonucleotide to the 5′ end of the template to generate a plurality of template fragments of known sequence with common sequences at the 5′ and 3′ ends;
(b) treat with uracil DNA glycosidase to fragment the uracil containing probes;
(c) treat with bisulfite to obtain bisulfite modified targets;
(d) amplify the bisulfite modified targets using primers to said first and second common sequences; and
(e) analyze cytosine positions in the target sequence for methylation by determining the sequence of the position in the amplified bisulfite modified targets obtained in step (d).
2. The method of claim 1 wherein step (e) comprises: hybridization of the amplified target to a pair of sequence specific probes for each cytosine position to be analyzed wherein the pair of methylation specific probes comprises a first probe that is perfectly complementary to the cytosine position and surrounding bases with bisulfite modification blocked by methylation and a second probe that is perfectly complementary to the cytosine position and surrounding bases with bisulfite modification, wherein the presence or absence of methylation is determined by analyzing the pattern of hybridization of the amplified bisulfite modified target to the first and second probes.
3. The method of claim 1 wherein sequence specific probes are attached to a solid support, selected from the group consisting of plurality of beads, one or more silica chip, one or more glass slides and one or more membranes, and wherein there are at least 1,000 pairs of sequence specific probes attached to said solid support.
4. A method for determining the methylation status of a plurality of cytosines in a plurality of restriction sites for a methylation sensitive restriction enzyme (MSRE) in a genomic DNA sample comprising:
(a) fragmenting the genomic DNA sample with the MSRE;
(b) filling the ends generated by cleavage with the MSRE using a DNA polymerase;
(c) fragmenting the product of (b) with a methylation insensitive restriction enzyme (MIRE) that recognizes the same restriction site;
(d) hybridize the products of (c) with methylation specific dU probes, wherein each comprises a barcode tag sequence, a 5′ first common priming sequence and a 3′ second common priming sequence;
(e) add a ligase and oligonucleotides that are complementary to the barcode tag sequences, the first common priming sequence and the second common priming sequence under conditions to allow ligation of the target to the oligonucleotides;
(f) digest the products of (e) with uracil DNA glycosidase;
(g) amplify the product of (f) by PCR using primers to the first and second common priming sequence; and,
(h) detect the presence or absence of amplified barcode tag sequences from (g) by hybridization to determine methylation status of selected cytosines.
5. The method of claim 4 wherein for each restriction site being interrogated there is a first dU probe and a second dU probe, wherein the first dU probe is complementary to the fragment generated if the MSRE cut at the restriction site and the second dU probe is complementary to the fragment generated if the MSRE did not cut at the restriction site and wherein each dU probe has a different barcode tag sequence.
6. The method of claim 5 wherein the presence of the barcode sequence of the first dU probe indicates that the restriction site was unmethylated in at least some fragments and presence of the barcode tag sequence of the second dU probe indicates that the restriction site was methylated in at least some fragments.
7. The method of claim 4 wherein the barcode tag sequences are detected by hybridization to an array of barcode tag complements in known or determinable positions of the array.
8. The method of claim 7 wherein the array comprises more than 10,000 different barcode tag complements.
9. The method of claim 4 wherein the barcode tag sequence is between 20 and 30 bases in length.
10. The method of claim 4 wherein the DNA polymerase is Klenow.
11. The method of claim 4 wherein at least 500 different restriction sites are analyzed.
12. The method of claim 4 wherein the MSRE is Hpa II or Sma I and the MIRE is Msp I or Xma I.
13. The method of claim 7 wherein the barcode tag complements are attached to a solid support selected from the group consisting of beads, chips, and membranes.
14. A method of analyzing methylation status of a plurality of restriction sites comprising:
(a) cutting the genomic DNA with a MDRE;
(b) end filling the restriction sites that were cut with the MDRE;
(c) cutting the product of (b) with a MSRE that is an isoschizomers of the MDRE;
(d) adding first and second dU probes to the products of (c) along with ligase, tag oligonucleotides and common sequence oligonucleotides;
(e) treating the products of (d) with UDG to fragment dU probes;
(f) amplifying the products of (e) by PCR with common primers; and
(g) detecting the tag sequences present in (f).
15. The method of claim 14 wherein the MDRE is BisI and the MSRE is Fnu4HI.
16. A method of detecting the presence or absence of methylation at a restriction site in a genomic DNA sample, said method comprising:
(a) treating said genomic DNA sample with a methylation sensitive restriction enzyme;
(b) treating the products of (a) with Klenow;
(c) treating the products of (b) with a methylation insensitive isoschizomers of the methylation sensitive restriction enzyme used in (a);
(d) adding to the products of (c) the following nucleic acids:
(i) a first dU probe and a second dU probe wherein said first dU probe is complementary to a strand of genomic DNA immediately adjacent to the restriction site and including the region filled in by Klenow and said second dU probe is complementary to the same strand of genomic DNA but does not include the region filled in by Klenow and wherein said first and second dU probes further comprise first and second common priming regions wherein said first dU probe has a first tag complement region and said second dU probe has a second tag complement region,
(ii) a first tag oligonucleotide that is complementary to said first tag complement region and a second tag complement that is complementary to said second tag complement;
(iii) oligonucleotides complementary to the first and second common priming regions; and
(iv) ligase;
(d) incubating to allow ligation;
(e) adding a UDG activity to the product of (d);
(f) amplifying the product of (e) by PCR using primers to the common priming regions; and
(g) detecting the presence of said first tag sequence or said second tag sequence, wherein presence of said first tag sequence indicates that the restriction site was unmethylated and presence of said second tag sequence indicates that the restriction site was methylated.
17. The method of claim 16 wherein said step of detecting comprises hybridization to an array of probes complementary to tag sequences.
18. A method of detecting the presence or absence of methylation at a plurality of restriction sites according to the method of claim 16.
19. The method of claim 18 wherein the plurality of restriction sites includes at least 200 different restriction sites.
20. A method for analyzing the methylation of a plurality of cytosines in a plurality of target sequences, said method comprising:
(a) obtaining a genomic DNA sample;
(b) fragmenting the genomic DNA sample to obtain fragments, wherein said fragments comprise a mixture of target fragments and non-target fragments;
(c) mixing the fragments with (i) a plurality of template probes, each template probe comprising a target complementarity region, a first common priming sequence and a second common priming sequence, wherein said first common priming sequence is 3′ of the target complementarity region and said second common priming sequences is 5′ of the target complementarity region, (ii) an oligonucleotide complementary to said first common priming sequence and an oligonucleotide that is complementary to said second common priming sequence, wherein at least one of the oligonucleotides is exonuclease resistant, and (iii) a ligase, to obtain ligation products;
(d) treating the products of (c) with an exonuclease to digest the template probes;
(e) treating the products of (d) with bisulfite;
(f) amplifying the products of (e) by PCR using primers to the common primer sequences to obtain an amplification product;
(g) hybridizing the amplification product from step (g) to an array of probes to obtain a hybridization pattern; and
(h) analyzing the hybridization pattern to detect the presence or absence of methylation at a plurality of cytosines in the target sequences.
21. The method of claim 20 wherein the exonuclease is a 5′ to 3′ exonuclease and the oligonucleotide complementary to said first common priming sequence is resistant to 5′ to 3′ exonuclease digestion.
22. The method of claim 20 wherein the oligonucleotide complementary to said first common priming sequence comprises a plurality of phosphorothioate linkages.
23. The method of claim 20 wherein the exonuclease is a 3′ to 5′ exonuclease and the oligonucleotide complementary to said second common priming sequence is resistant to 3′ to 5′ exonuclease digestion.
24. The method of claim 20 wherein the array comprises a first plurality of probes that are perfectly complementary to a genomic region containing a CpG after bisulfite treatment and amplification if the region was unmethylated and a second plurality of probes that are perfectly complementary to the same region after bisulfite modification and amplification if the region was fully methylated.
25. The method of claim 20 wherein the common primer sequences contain 5 methyl cytosine in place of cytosine.
US11/923,649 2001-11-19 2007-10-24 Methods of Analysis of Methylation Abandoned US20110151438A9 (en)

Priority Applications (4)

Application Number Priority Date Filing Date Title
US11/923,649 US20110151438A9 (en) 2001-11-19 2007-10-24 Methods of Analysis of Methylation
US15/281,721 US10407717B2 (en) 2001-11-19 2016-09-30 Methods of analysis of methylation
US16/528,222 US10822642B2 (en) 2001-11-19 2019-07-31 Methods of analysis of methylation
US17/080,813 US20210147903A1 (en) 2001-11-19 2020-10-26 Methods of analysis of methylation

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
US33169301P 2001-11-19 2001-11-19
US10/300,311 US7208295B2 (en) 2001-11-19 2002-11-19 Multiplex oligonucleotide addition and target amplification
US86273506P 2006-10-24 2006-10-24
US11/739,654 US7754451B2 (en) 2001-11-19 2007-04-24 Multiplex oligonucleotide addition and target amplification
US11/923,649 US20110151438A9 (en) 2001-11-19 2007-10-24 Methods of Analysis of Methylation

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
US11/739,654 Continuation-In-Part US7754451B2 (en) 2001-11-19 2007-04-24 Multiplex oligonucleotide addition and target amplification

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US15/281,721 Division US10407717B2 (en) 2001-11-19 2016-09-30 Methods of analysis of methylation

Publications (2)

Publication Number Publication Date
US20080108073A1 US20080108073A1 (en) 2008-05-08
US20110151438A9 true US20110151438A9 (en) 2011-06-23

Family

ID=39360154

Family Applications (4)

Application Number Title Priority Date Filing Date
US11/923,649 Abandoned US20110151438A9 (en) 2001-11-19 2007-10-24 Methods of Analysis of Methylation
US15/281,721 Expired - Lifetime US10407717B2 (en) 2001-11-19 2016-09-30 Methods of analysis of methylation
US16/528,222 Expired - Lifetime US10822642B2 (en) 2001-11-19 2019-07-31 Methods of analysis of methylation
US17/080,813 Abandoned US20210147903A1 (en) 2001-11-19 2020-10-26 Methods of analysis of methylation

Family Applications After (3)

Application Number Title Priority Date Filing Date
US15/281,721 Expired - Lifetime US10407717B2 (en) 2001-11-19 2016-09-30 Methods of analysis of methylation
US16/528,222 Expired - Lifetime US10822642B2 (en) 2001-11-19 2019-07-31 Methods of analysis of methylation
US17/080,813 Abandoned US20210147903A1 (en) 2001-11-19 2020-10-26 Methods of analysis of methylation

Country Status (1)

Country Link
US (4) US20110151438A9 (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017201102A1 (en) 2016-05-16 2017-11-23 Accuragen Holdings Limited Method of improved sequencing by strand identification
US10155980B2 (en) 2016-08-15 2018-12-18 Accuragen Holdings Limited Compositions and methods for detecting rare sequence variants
US10752942B2 (en) 2015-10-09 2020-08-25 Accuragen Holdings Limited Methods and compositions for enrichment of amplification products
US11203782B2 (en) 2018-03-29 2021-12-21 Accuragen Holdings Limited Compositions and methods comprising asymmetric barcoding
US11286519B2 (en) 2013-12-11 2022-03-29 Accuragen Holdings Limited Methods and compositions for enrichment of amplification products
US11597973B2 (en) 2013-12-11 2023-03-07 Accuragen Holdings Limited Compositions and methods for detecting rare sequence variants
US11859246B2 (en) 2013-12-11 2024-01-02 Accuragen Holdings Limited Methods and compositions for enrichment of amplification products

Families Citing this family (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110151438A9 (en) * 2001-11-19 2011-06-23 Affymetrix, Inc. Methods of Analysis of Methylation
JP5151167B2 (en) * 2007-01-31 2013-02-27 住友化学株式会社 DNA methylation measurement method
US20090233809A1 (en) * 2008-03-04 2009-09-17 Affymetrix, Inc. Resequencing methods for identification of sequence variants
US8541207B2 (en) 2008-10-22 2013-09-24 Illumina, Inc. Preservation of information related to genomic DNA methylation
EP2464748A4 (en) * 2009-08-12 2012-12-05 Commw Scient Ind Res Org Detection and analysis of methylation in nucleic acid sequences
US8759036B2 (en) 2011-03-21 2014-06-24 Affymetrix, Inc. Methods for synthesizing pools of probes
US10174383B2 (en) * 2014-08-13 2019-01-08 Vanadis Diagnostics Method of estimating the amount of a methylated locus in a sample
US10683534B2 (en) * 2015-01-27 2020-06-16 BioSpyder Technologies, Inc. Ligation assays in liquid phase
EP3209801B1 (en) * 2014-10-20 2023-06-28 Commonwealth Scientific and Industrial Research Organisation Genome methylation analysis
CN107868774A (en) * 2016-09-23 2018-04-03 上海紫众生物科技有限公司 A kind of thermal starting archaeal dna polymerase, preparation method and applications
KR102476709B1 (en) * 2016-11-21 2022-12-09 나노스트링 테크놀로지스, 인크. Chemical compositions and methods of using same
KR20210061962A (en) 2018-05-14 2021-05-28 나노스트링 테크놀로지스, 인크. Chemical composition and method of use thereof
CN111850101B (en) * 2020-06-29 2021-12-28 西安交通大学 Visual distinguishing method for single-cell DNA epigenetic modification space positioning and adjacent distribution
CN113308518B (en) * 2021-06-02 2022-10-11 重庆大学 DNA methylation hypersensitive detection method and application thereof
WO2023082240A1 (en) * 2021-11-15 2023-05-19 深圳华大智造科技股份有限公司 Method for detecting genome-wide dna methylation modification
CN114214408B (en) * 2021-12-22 2023-02-17 重庆大学附属肿瘤医院 Method, probe library and kit for detecting tumor ctDNA methylation with high throughput and high sensitivity

Citations (67)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5118604A (en) * 1990-06-15 1992-06-02 Yale University Restriction-fragment DNA probes and probe clusters
US5217889A (en) * 1990-10-19 1993-06-08 Roninson Igor B Methods and applications for efficient genetic suppressor elements
US5580730A (en) * 1994-08-19 1996-12-03 Olympus America, Inc. Enzyme digestion method for the detection of amplified DNA
US5710000A (en) * 1994-09-16 1998-01-20 Affymetrix, Inc. Capturing sequences adjacent to Type-IIs restriction sites for genomic library mapping
US5786146A (en) * 1996-06-03 1998-07-28 The Johns Hopkins University School Of Medicine Method of detection of methylated nucleic acid using agents which modify unmethylated cytosine and distinguishing modified methylated and non-methylated nucleic acids
US5871917A (en) * 1996-05-31 1999-02-16 North Shore University Hospital Research Corp. Identification of differentially methylated and mutated nucleic acids
US5912147A (en) * 1996-10-22 1999-06-15 Health Research, Inc. Rapid means of quantitating genomic instability
US5942609A (en) * 1998-11-12 1999-08-24 The Porkin-Elmer Corporation Ligation assembly and detection of polynucleotides on solid-support
US6017704A (en) * 1996-06-03 2000-01-25 The Johns Hopkins University School Of Medicine Method of detection of methylated nucleic acid using agents which modify unmethylated cytosine and distinguishing modified methylated and non-methylated nucleic acids
US6156501A (en) * 1993-10-26 2000-12-05 Affymetrix, Inc. Arrays of modified nucleic acid probes and methods of use
US6214556B1 (en) * 1997-11-27 2001-04-10 Epigenomics Ag Method for producing complex DNA methylation fingerprints
US6251594B1 (en) * 1997-06-09 2001-06-26 Usc/Norris Comprehensive Cancer Ctr. Cancer diagnostic method based upon DNA methylation differences
US6300071B1 (en) * 1998-07-29 2001-10-09 Keygene N.V. Method for detecting nucleic acid methylation using AFLP™
US6309823B1 (en) * 1993-10-26 2001-10-30 Affymetrix, Inc. Arrays of nucleic acid probes for analyzing biotransformation genes and methods of using the same
US20010046669A1 (en) * 1999-02-24 2001-11-29 Mccobmie William R. Genetically filtered shotgun sequencing of complex eukaryotic genomes
US6331393B1 (en) * 1999-05-14 2001-12-18 University Of Southern California Process for high-throughput DNA methylation analysis
US20020006623A1 (en) * 1998-05-04 2002-01-17 Allan Bradley Compositions and methods for array-based genomic nucleic acid analysis of biological molecules
US6383754B1 (en) * 1999-08-13 2002-05-07 Yale University Binary encoded sequence tags
US20020123053A1 (en) * 2000-06-23 2002-09-05 University Of Chicago Methods for isolating centromere DNA
US20020137086A1 (en) * 2001-03-01 2002-09-26 Alexander Olek Method for the development of gene panels for diagnostic and therapeutic purposes based on the expression and methylation status of the genes
US20020197639A1 (en) * 2001-06-08 2002-12-26 Shia Michael A. Methods and products for analyzing nucleic acids based on methylation status
US20030022215A1 (en) * 2001-04-23 2003-01-30 Dana-Farber Cancer Institute, Inc. Methods for rapid screening of polymorphisms, mutations and methylation
US6514698B1 (en) * 1997-08-29 2003-02-04 Osvaldo J. Lopez DNA methyltransferase genotyping
US20030032026A1 (en) * 1999-07-26 2003-02-13 Kurt Berlin Method for relative quantification of methylation of cytosine bases in DNA samples
US20030036081A1 (en) * 2001-07-02 2003-02-20 Epigenomics Ag Distributed system for epigenetic based prediction of complex phenotypes
US20030068620A1 (en) * 2000-09-20 2003-04-10 Sanford Markowitz Methods and compositions for detecting cancers associated with methylation of hMLH1 promoter DNA
US20030082600A1 (en) * 2001-03-09 2003-05-01 Alexander Olek Highly sensitive method for the detection of cytosine methylation patters
US20030082609A1 (en) * 2000-04-06 2003-05-01 Alexander Olek Diagnosis of diseases associated with gene regulation
US6562569B1 (en) * 1998-12-30 2003-05-13 Oligos Etc. Inc. Arrays with modified oligonucleotide and polynucleotide compositions
US20030096291A1 (en) * 2001-11-19 2003-05-22 Malek Faham Multiplex oligonucleotide addition and target amplification
US20030096289A1 (en) * 2001-11-16 2003-05-22 Osamu Suzuki Oligonucleotide-immobilized substrate for detecting methylation
US20030099997A1 (en) * 2001-10-24 2003-05-29 Bestor Timothy H. Method for gene identification based on differential DNA methylation
US20030104464A1 (en) * 2000-06-19 2003-06-05 Kurt Berlin Method for the high-parallel analysis of polymorphisms
US20030113737A1 (en) * 2001-01-24 2003-06-19 Genomic Expression Aps Assay and kit for analyzing gene expression
US20030119025A1 (en) * 2000-02-25 2003-06-26 Alexander Olek Ligase/polymerase method for detecting cytosine methylation in dna samples
US20030129620A1 (en) * 2000-02-25 2003-07-10 Alexander Olek Method for detecting cytosine methylation in dna samples
US20030129602A1 (en) * 1999-02-05 2003-07-10 Huang Tim Hui-Ming High-throughput methods for detecting DNA methylation
US20030143606A1 (en) * 2000-06-30 2003-07-31 Alexander Olek Diagnosis of diseases associated with the immune system by determining cytosine methylation
US20030148290A1 (en) * 2002-02-06 2003-08-07 Susan Cottrell Quantitative methylation detection in DNA samples
US6613511B1 (en) * 1997-04-21 2003-09-02 Xzillion Gmbh & Co. Characterizing DNA
US20030170684A1 (en) * 2000-02-07 2003-09-11 Jian-Bing Fan Multiplexed methylation detection methods
US20030170689A1 (en) * 2001-05-11 2003-09-11 Rexagen Corporation DNA microarrays comprising active chromatin elements and comprehensive profiling therewith
US20030180779A1 (en) * 2002-03-15 2003-09-25 Epigenomics Ag Discovery and diagnostic methods using 5-methylcytosine DNA glycosylase
US20030212455A1 (en) * 2000-03-16 2003-11-13 Bas Van Steensel Identification of in vivo dna binding loci of chromatin proteins using a tethered nucleotide modification enzyme
US20030215842A1 (en) * 2002-01-30 2003-11-20 Epigenomics Ag Method for the analysis of cytosine methylation patterns
US20030224040A1 (en) * 2002-03-07 2003-12-04 Baylin Stephen B. Genomic screen for epigenetically silenced genes associated with cancer
US6677121B2 (en) * 1999-04-06 2004-01-13 Agilix Corporation Fixed address analysis of sequence tags
US20040029128A1 (en) * 2002-08-08 2004-02-12 Epigenomics, Inc. Methods and nucleic acids for the analysis of CpG dinucleotide methylation status associated with the calcitonin gene
US20040038254A1 (en) * 2002-03-08 2004-02-26 Risa Peoples Compositions and methods for detecting nucleic acid methylation
US20040048275A1 (en) * 2000-10-23 2004-03-11 Per Guldberg Materials and methods relating to nucleic acid amplification and profiling
US6713258B2 (en) * 2000-03-29 2004-03-30 Center For The Application Of Molecular Biology To International Agriculture (Cambia) Methods for genotyping by hybridization analysis
US20040072197A1 (en) * 2001-11-08 2004-04-15 Jones Peter A. Assay for the detection and quantitation of hemimethylation
US20040081976A1 (en) * 2002-03-07 2004-04-29 David Sidransky Genomic screen for epigenetically silenced tumor suppressor genes
US20040086944A1 (en) * 2000-11-13 2004-05-06 Grigg Geoffrey Walter Detection of methylated dna molecules
US20040102905A1 (en) * 2001-03-26 2004-05-27 Epigenomics Ag Method for epigenetic feature selection
US20040115663A1 (en) * 2001-10-26 2004-06-17 Kurt Berlin Method for the detection of cytosine methylations in immobilized dna samples
US20040146868A1 (en) * 2003-01-24 2004-07-29 Epigenomics Ag Methods and nucleic acids for the analysis of CpG dinucleotide methylation status associated with the development of peripheral zone prostate cancer
US20040152080A1 (en) * 2000-06-19 2004-08-05 Kurt Berlin Method for detecting cytosine methylations
US20040248090A1 (en) * 1999-12-06 2004-12-09 Alexander Olek Method for the parallel detection of the degree of methylation of genomic dna
US20050009059A1 (en) * 2003-05-07 2005-01-13 Affymetrix, Inc. Analysis of methylation status using oligonucleotide arrays
US20050021240A1 (en) * 2000-11-02 2005-01-27 Epigenomics Ag Systems, methods and computer program products for guiding selection of a therapeutic treatment regimen based on the methylation status of the DNA
US20050064428A1 (en) * 2002-01-08 2005-03-24 Kurt Berlin Method for detecting cytosine-methylation patterns by exponential ligation of hybridised probe oligo-nucleotides (mla)
US6884586B2 (en) * 2001-07-15 2005-04-26 Keck Graduate Institute Methylation analysis using nicking agents
US20050196792A1 (en) * 2004-02-13 2005-09-08 Affymetrix, Inc. Analysis of methylation status using nucleic acid arrays
US20050227265A1 (en) * 2004-02-10 2005-10-13 Francis Barany Method for detection of promoter methylation status
US20050266417A1 (en) * 2003-09-12 2005-12-01 Francis Barany Methods for identifying target nucleic acid molecules
US20060292585A1 (en) * 2005-06-24 2006-12-28 Affymetrix, Inc. Analysis of methylation using nucleic acid arrays

Family Cites Families (59)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TW205070B (en) 1986-06-30 1993-05-01 Takeda Pharm Industry Co Ltd
US5185243A (en) 1988-08-25 1993-02-09 Syntex (U.S.A.) Inc. Method for detection of specific nucleic acid sequences
CA1339731C (en) 1988-10-12 1998-03-17 Charles T. Caskey Multiplex genomic dna amplification for deletion detection
US5508178A (en) 1989-01-19 1996-04-16 Rose; Samuel Nucleic acid amplification using single primer
GB8920097D0 (en) 1989-09-06 1989-10-18 Ici Plc Amplification processes
US5439793A (en) 1990-07-19 1995-08-08 Syntex (U.S.A.) Inc. Method for producing a polynucleotide having an intramolecularly base-paired structure
US5595891A (en) 1990-07-19 1997-01-21 Behringwerke Ag Method for producing a polynucleotide for use in single primer amplification
DE69230873T3 (en) 1991-09-24 2006-04-06 Keygene N.V. Selective restriction fragment amplification: general procedure for DNA fingerprinting
US5612199A (en) 1991-10-11 1997-03-18 Behringwerke Ag Method for producing a polynucleotide for use in single primer amplification
US5731171A (en) 1993-07-23 1998-03-24 Arch Development Corp. Sequence independent amplification of DNA
US5928905A (en) 1995-04-18 1999-07-27 Glaxo Group Limited End-complementary polymerase reaction
US6117679A (en) 1994-02-17 2000-09-12 Maxygen, Inc. Methods for generating polynucleotides having desired characteristics by iterative selection and recombination
US5834252A (en) 1995-04-18 1998-11-10 Glaxo Group Limited End-complementary polymerase reaction
US5851770A (en) 1994-04-25 1998-12-22 Variagenics, Inc. Detection of mismatches by resolvase cleavage using a magnetic bead support
ATE226983T1 (en) 1994-08-19 2002-11-15 Pe Corp Ny COUPLED AMPLICATION AND LIGATION PROCEDURE
US5491063A (en) 1994-09-01 1996-02-13 Hoffmann-La Roche Inc. Methods for in-solution quenching of fluorescently labeled oligonucleotide probes
US5846719A (en) 1994-10-13 1998-12-08 Lynx Therapeutics, Inc. Oligonucleotide tags for sorting and identification
US5565340A (en) 1995-01-27 1996-10-15 Clontech Laboratories, Inc. Method for suppressing DNA fragment amplification during PCR
US5882856A (en) 1995-06-07 1999-03-16 Genzyme Corporation Universal primer sequence for multiplex DNA amplification
US5728526A (en) 1995-06-07 1998-03-17 Oncor, Inc. Method for analyzing a nucleotide sequence
PT843735E (en) 1995-07-11 2002-11-29 Forfas Trading As Biores Irela DETECTION OF GLYCOSYLASE-NUCLEOTIDED NUCLEOTIDED SEQUENCES
US5888731A (en) 1995-08-30 1999-03-30 Visible Genetics Inc. Method for identification of mutations using ligation of multiple oligonucleotide probes
US6458530B1 (en) 1996-04-04 2002-10-01 Affymetrix Inc. Selecting tag nucleic acids
EP0912761A4 (en) 1996-05-29 2004-06-09 Cornell Res Foundation Inc Detection of nucleic acid sequence differences using coupled ligase detection and polymerase chain reactions
WO1999049293A2 (en) 1998-03-24 1999-09-30 Boston Probes, Inc. Methods, kits and compositions pertaining to detection complexes
GB9812768D0 (en) 1998-06-13 1998-08-12 Zeneca Ltd Methods
US6355422B1 (en) 1998-06-24 2002-03-12 City Of Hope Single tube PCR assay for detection of chromosomal mutations: application to the inversion hotspot in the factor VIII gene including optional use of subcycling PCR
US6743605B1 (en) 1998-06-24 2004-06-01 Enzo Life Sciences, Inc. Linear amplification of specific nucleic acid sequences
ATE374833T1 (en) 1998-11-09 2007-10-15 Eiken Chemical METHOD FOR SYNTHESIS OF NUCLEIC ACIDS
JP4087028B2 (en) 1999-06-14 2008-05-14 シャープ株式会社 External circuit mounting method and thermocompression bonding apparatus
US20020119448A1 (en) 1999-06-23 2002-08-29 Joseph A. Sorge Methods of enriching for and identifying polymorphisms
ATE409236T1 (en) 1999-07-14 2008-10-15 Packard Bioscience Company DERIVATIZED NUCLEIC ACIDS AND THEIR USE
US7118860B2 (en) 1999-10-29 2006-10-10 Stratagene California Methods for detection of a target nucleic acid by capture
US6528254B1 (en) 1999-10-29 2003-03-04 Stratagene Methods for detection of a target nucleic acid sequence
US6207424B1 (en) 1999-11-23 2001-03-27 Maxim Biotech, Inc. Self-primed amplification system
DE19957827C2 (en) 1999-11-25 2003-06-12 Epigenomics Ag Use of an oligomer array with PNA and / or DNA oligomers on a surface
CA2399733C (en) 2000-02-07 2011-09-20 Illumina, Inc. Nucleic acid detection methods using universal priming
JP2003535599A (en) 2000-06-06 2003-12-02 ティーエム バイオサイエンス コーポレイション Nucleic acid capture unit and its use
AU2001281311A1 (en) 2000-07-10 2002-01-21 Epigenx Pharmaceutical, Inc. Detecting methylated cytosine in polynucleotides
CA2421102A1 (en) 2000-08-25 2002-03-07 Lovelace Respiratory Research Institute Nested methylation-specific polymerase chain reaction cancer detection method
US7129044B2 (en) 2000-10-04 2006-10-31 The Board Of Trustees Of The Leland Stanford Junior University Renaturation, reassociation, association and hybridization of nucleic acid molecules
US6350580B1 (en) 2000-10-11 2002-02-26 Stratagene Methods for detection of a target nucleic acid using a probe comprising secondary structure
AU1070002A (en) 2000-10-23 2002-05-06 Cancer Res Ventures Ltd Materials and methods relating to nucleic acid amplification and profiling
IE20000887A1 (en) 2000-11-03 2002-12-11 Univ College Cork Nat Univ Ie Method for the amplification and optional characterisation of nucleic acids
EP2327794A3 (en) 2001-01-25 2011-09-14 TM Bioscience Corporation Polynucleotides for use as tags and tag complements, manufacture and use thereof
US7101663B2 (en) 2001-03-02 2006-09-05 University of Pittsburgh—of the Commonwealth System of Higher Education PCR method
US20060194275A1 (en) 2001-04-13 2006-08-31 Incyte Corporation Transporter and ion channels
WO2003023065A1 (en) 2001-09-06 2003-03-20 Syngenta Participations Ag Dna methylation patterns
EP1438428A4 (en) 2001-09-14 2005-11-16 Univ Queensland Detection of dna methylation
US20060183115A1 (en) 2001-09-26 2006-08-17 Epigenx Pharmaceuticals Assays for dna methylation changes
EP1442129B1 (en) 2001-10-18 2008-09-10 Lovelace Respiratory Research Institute Cancer monitoring by aberrant promotor methylation of the transcription factor genes pax5 alpha pax5 beta, novel loop helix loop protein, novel gene 2, and beta3 genes
US20110151438A9 (en) * 2001-11-19 2011-06-23 Affymetrix, Inc. Methods of Analysis of Methylation
CA2467814A1 (en) 2001-12-03 2003-06-12 Illumina, Inc. Multiplexed methylation detection methods
EP1344832A1 (en) 2002-03-15 2003-09-17 Epigenomics AG Methods and nucleic acids for the analysis of methylation within the gene melastatin
AU2002952944A0 (en) 2002-11-27 2002-12-12 Human Genetic Signatures Pty Ltd Restoration of methylation states in cells
US7288373B2 (en) 2003-05-02 2007-10-30 Human Genetic Signatures Pty Ltd. Treatment of methylated nucleic acid
US7368239B2 (en) 2003-08-29 2008-05-06 Applera Corporation Method and materials for polyamine catalyzed bisulfite conversion of cytosine to uracil
US7534873B2 (en) 2003-08-29 2009-05-19 Applied Biosystems, Llc Method and materials for quaternary amine catalyzed bisulfite conversion of cytosine to uracil
WO2020218649A1 (en) 2019-04-26 2020-10-29 주식회사 카이렌 Oxygen mask including temperature adjusting module

Patent Citations (77)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5118604A (en) * 1990-06-15 1992-06-02 Yale University Restriction-fragment DNA probes and probe clusters
US5217889A (en) * 1990-10-19 1993-06-08 Roninson Igor B Methods and applications for efficient genetic suppressor elements
US6309823B1 (en) * 1993-10-26 2001-10-30 Affymetrix, Inc. Arrays of nucleic acid probes for analyzing biotransformation genes and methods of using the same
US6156501A (en) * 1993-10-26 2000-12-05 Affymetrix, Inc. Arrays of modified nucleic acid probes and methods of use
US5580730A (en) * 1994-08-19 1996-12-03 Olympus America, Inc. Enzyme digestion method for the detection of amplified DNA
US6509160B1 (en) * 1994-09-16 2003-01-21 Affymetric, Inc. Methods for analyzing nucleic acids using a type IIs restriction endonuclease
US5710000A (en) * 1994-09-16 1998-01-20 Affymetrix, Inc. Capturing sequences adjacent to Type-IIs restriction sites for genomic library mapping
US5871917A (en) * 1996-05-31 1999-02-16 North Shore University Hospital Research Corp. Identification of differentially methylated and mutated nucleic acids
US6017704A (en) * 1996-06-03 2000-01-25 The Johns Hopkins University School Of Medicine Method of detection of methylated nucleic acid using agents which modify unmethylated cytosine and distinguishing modified methylated and non-methylated nucleic acids
US6200756B1 (en) * 1996-06-03 2001-03-13 The Johns Hopkins University School Of Medicine Methods for identifying methylation patterns in a CpG-containing nucleic acid
US5786146A (en) * 1996-06-03 1998-07-28 The Johns Hopkins University School Of Medicine Method of detection of methylated nucleic acid using agents which modify unmethylated cytosine and distinguishing modified methylated and non-methylated nucleic acids
US5912147A (en) * 1996-10-22 1999-06-15 Health Research, Inc. Rapid means of quantitating genomic instability
US6613511B1 (en) * 1997-04-21 2003-09-02 Xzillion Gmbh & Co. Characterizing DNA
US6251594B1 (en) * 1997-06-09 2001-06-26 Usc/Norris Comprehensive Cancer Ctr. Cancer diagnostic method based upon DNA methylation differences
US6514698B1 (en) * 1997-08-29 2003-02-04 Osvaldo J. Lopez DNA methyltransferase genotyping
US6214556B1 (en) * 1997-11-27 2001-04-10 Epigenomics Ag Method for producing complex DNA methylation fingerprints
US20020006623A1 (en) * 1998-05-04 2002-01-17 Allan Bradley Compositions and methods for array-based genomic nucleic acid analysis of biological molecules
US6300071B1 (en) * 1998-07-29 2001-10-09 Keygene N.V. Method for detecting nucleic acid methylation using AFLP™
US5942609A (en) * 1998-11-12 1999-08-24 The Porkin-Elmer Corporation Ligation assembly and detection of polynucleotides on solid-support
US6562569B1 (en) * 1998-12-30 2003-05-13 Oligos Etc. Inc. Arrays with modified oligonucleotide and polynucleotide compositions
US6605432B1 (en) * 1999-02-05 2003-08-12 Curators Of The University Of Missouri High-throughput methods for detecting DNA methylation
US20030129602A1 (en) * 1999-02-05 2003-07-10 Huang Tim Hui-Ming High-throughput methods for detecting DNA methylation
US20010046669A1 (en) * 1999-02-24 2001-11-29 Mccobmie William R. Genetically filtered shotgun sequencing of complex eukaryotic genomes
US6677121B2 (en) * 1999-04-06 2004-01-13 Agilix Corporation Fixed address analysis of sequence tags
US20020086324A1 (en) * 1999-05-14 2002-07-04 Laird Peter W. Process for high throughput DNA methylation analysis
US6331393B1 (en) * 1999-05-14 2001-12-18 University Of Southern California Process for high-throughput DNA methylation analysis
US20030032026A1 (en) * 1999-07-26 2003-02-13 Kurt Berlin Method for relative quantification of methylation of cytosine bases in DNA samples
US6383754B1 (en) * 1999-08-13 2002-05-07 Yale University Binary encoded sequence tags
US20040248090A1 (en) * 1999-12-06 2004-12-09 Alexander Olek Method for the parallel detection of the degree of methylation of genomic dna
US20030170684A1 (en) * 2000-02-07 2003-09-11 Jian-Bing Fan Multiplexed methylation detection methods
US20030119025A1 (en) * 2000-02-25 2003-06-26 Alexander Olek Ligase/polymerase method for detecting cytosine methylation in dna samples
US20030129620A1 (en) * 2000-02-25 2003-07-10 Alexander Olek Method for detecting cytosine methylation in dna samples
US20030212455A1 (en) * 2000-03-16 2003-11-13 Bas Van Steensel Identification of in vivo dna binding loci of chromatin proteins using a tethered nucleotide modification enzyme
US6713258B2 (en) * 2000-03-29 2004-03-30 Center For The Application Of Molecular Biology To International Agriculture (Cambia) Methods for genotyping by hybridization analysis
US20030148327A1 (en) * 2000-04-06 2003-08-07 Alexander Olek Diagnosis of diseases associated with metastasis
US20030148326A1 (en) * 2000-04-06 2003-08-07 Alexander Olek Diagnosis of diseases associated with dna transcription
US20030162194A1 (en) * 2000-04-06 2003-08-28 Alexander Olek Diagnosis of diseases associated with apoptosis
US20030082609A1 (en) * 2000-04-06 2003-05-01 Alexander Olek Diagnosis of diseases associated with gene regulation
US20040152080A1 (en) * 2000-06-19 2004-08-05 Kurt Berlin Method for detecting cytosine methylations
US20030104464A1 (en) * 2000-06-19 2003-06-05 Kurt Berlin Method for the high-parallel analysis of polymorphisms
US20040023282A1 (en) * 2000-06-23 2004-02-05 Song Luo Methods for isolating centromere DNA
US6649347B2 (en) * 2000-06-23 2003-11-18 Chromatin, Inc. Use of methylated nucleic acid segments for isolating centromere DNA
US20020123053A1 (en) * 2000-06-23 2002-09-05 University Of Chicago Methods for isolating centromere DNA
US20030143606A1 (en) * 2000-06-30 2003-07-31 Alexander Olek Diagnosis of diseases associated with the immune system by determining cytosine methylation
US20040023230A1 (en) * 2000-06-30 2004-02-05 Alexander Olek Method and nucleic acids for pharmacogenomic methylation analysis
US20030068620A1 (en) * 2000-09-20 2003-04-10 Sanford Markowitz Methods and compositions for detecting cancers associated with methylation of hMLH1 promoter DNA
US20040048275A1 (en) * 2000-10-23 2004-03-11 Per Guldberg Materials and methods relating to nucleic acid amplification and profiling
US20050021240A1 (en) * 2000-11-02 2005-01-27 Epigenomics Ag Systems, methods and computer program products for guiding selection of a therapeutic treatment regimen based on the methylation status of the DNA
US20040086944A1 (en) * 2000-11-13 2004-05-06 Grigg Geoffrey Walter Detection of methylated dna molecules
US20030113737A1 (en) * 2001-01-24 2003-06-19 Genomic Expression Aps Assay and kit for analyzing gene expression
US20020137086A1 (en) * 2001-03-01 2002-09-26 Alexander Olek Method for the development of gene panels for diagnostic and therapeutic purposes based on the expression and methylation status of the genes
US20030082600A1 (en) * 2001-03-09 2003-05-01 Alexander Olek Highly sensitive method for the detection of cytosine methylation patters
US20040102905A1 (en) * 2001-03-26 2004-05-27 Epigenomics Ag Method for epigenetic feature selection
US20030022215A1 (en) * 2001-04-23 2003-01-30 Dana-Farber Cancer Institute, Inc. Methods for rapid screening of polymorphisms, mutations and methylation
US20030170689A1 (en) * 2001-05-11 2003-09-11 Rexagen Corporation DNA microarrays comprising active chromatin elements and comprehensive profiling therewith
US20020197639A1 (en) * 2001-06-08 2002-12-26 Shia Michael A. Methods and products for analyzing nucleic acids based on methylation status
US20030036081A1 (en) * 2001-07-02 2003-02-20 Epigenomics Ag Distributed system for epigenetic based prediction of complex phenotypes
US6884586B2 (en) * 2001-07-15 2005-04-26 Keck Graduate Institute Methylation analysis using nicking agents
US20030099997A1 (en) * 2001-10-24 2003-05-29 Bestor Timothy H. Method for gene identification based on differential DNA methylation
US20040115663A1 (en) * 2001-10-26 2004-06-17 Kurt Berlin Method for the detection of cytosine methylations in immobilized dna samples
US20040072197A1 (en) * 2001-11-08 2004-04-15 Jones Peter A. Assay for the detection and quantitation of hemimethylation
US20030096289A1 (en) * 2001-11-16 2003-05-22 Osamu Suzuki Oligonucleotide-immobilized substrate for detecting methylation
US20030096291A1 (en) * 2001-11-19 2003-05-22 Malek Faham Multiplex oligonucleotide addition and target amplification
US20050064428A1 (en) * 2002-01-08 2005-03-24 Kurt Berlin Method for detecting cytosine-methylation patterns by exponential ligation of hybridised probe oligo-nucleotides (mla)
US20030215842A1 (en) * 2002-01-30 2003-11-20 Epigenomics Ag Method for the analysis of cytosine methylation patterns
US20030148290A1 (en) * 2002-02-06 2003-08-07 Susan Cottrell Quantitative methylation detection in DNA samples
US20040081976A1 (en) * 2002-03-07 2004-04-29 David Sidransky Genomic screen for epigenetically silenced tumor suppressor genes
US20030224040A1 (en) * 2002-03-07 2003-12-04 Baylin Stephen B. Genomic screen for epigenetically silenced genes associated with cancer
US20040038254A1 (en) * 2002-03-08 2004-02-26 Risa Peoples Compositions and methods for detecting nucleic acid methylation
US20030180779A1 (en) * 2002-03-15 2003-09-25 Epigenomics Ag Discovery and diagnostic methods using 5-methylcytosine DNA glycosylase
US20040029128A1 (en) * 2002-08-08 2004-02-12 Epigenomics, Inc. Methods and nucleic acids for the analysis of CpG dinucleotide methylation status associated with the calcitonin gene
US20040146868A1 (en) * 2003-01-24 2004-07-29 Epigenomics Ag Methods and nucleic acids for the analysis of CpG dinucleotide methylation status associated with the development of peripheral zone prostate cancer
US20050009059A1 (en) * 2003-05-07 2005-01-13 Affymetrix, Inc. Analysis of methylation status using oligonucleotide arrays
US20050266417A1 (en) * 2003-09-12 2005-12-01 Francis Barany Methods for identifying target nucleic acid molecules
US20050227265A1 (en) * 2004-02-10 2005-10-13 Francis Barany Method for detection of promoter methylation status
US20050196792A1 (en) * 2004-02-13 2005-09-08 Affymetrix, Inc. Analysis of methylation status using nucleic acid arrays
US20060292585A1 (en) * 2005-06-24 2006-12-28 Affymetrix, Inc. Analysis of methylation using nucleic acid arrays

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Gitan et al. Genome Research (2001) 12: 158-164. *

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11286519B2 (en) 2013-12-11 2022-03-29 Accuragen Holdings Limited Methods and compositions for enrichment of amplification products
US11859246B2 (en) 2013-12-11 2024-01-02 Accuragen Holdings Limited Methods and compositions for enrichment of amplification products
US11597973B2 (en) 2013-12-11 2023-03-07 Accuragen Holdings Limited Compositions and methods for detecting rare sequence variants
US11578359B2 (en) 2015-10-09 2023-02-14 Accuragen Holdings Limited Methods and compositions for enrichment of amplification products
US10752942B2 (en) 2015-10-09 2020-08-25 Accuragen Holdings Limited Methods and compositions for enrichment of amplification products
US11427866B2 (en) 2016-05-16 2022-08-30 Accuragen Holdings Limited Method of improved sequencing by strand identification
EP3458586B1 (en) 2016-05-16 2022-12-28 Accuragen Holdings Limited Method of improved sequencing by strand identification
WO2017201102A1 (en) 2016-05-16 2017-11-23 Accuragen Holdings Limited Method of improved sequencing by strand identification
CN109511265A (en) * 2016-05-16 2019-03-22 安可济控股有限公司 The method for improving sequencing is identified by chain
US10724088B2 (en) 2016-08-15 2020-07-28 Accuragen Holdings Limited Compositions and methods for detecting rare sequence variants
US11643683B2 (en) 2016-08-15 2023-05-09 Accuragen Holdings Limited Compositions and methods for detecting rare sequence variants
US10155980B2 (en) 2016-08-15 2018-12-18 Accuragen Holdings Limited Compositions and methods for detecting rare sequence variants
US11203782B2 (en) 2018-03-29 2021-12-21 Accuragen Holdings Limited Compositions and methods comprising asymmetric barcoding

Also Published As

Publication number Publication date
US20210147903A1 (en) 2021-05-20
US10822642B2 (en) 2020-11-03
US20200024640A1 (en) 2020-01-23
US10407717B2 (en) 2019-09-10
US20080108073A1 (en) 2008-05-08
US20170081703A1 (en) 2017-03-23

Similar Documents

Publication Publication Date Title
US10822642B2 (en) Methods of analysis of methylation
US11821030B2 (en) Multiplex targeted amplification using flap nuclease
CA2950295C (en) Method for identification and enumeration of nucleic acid sequence, expression, copy, or dna methylation changes, using combined nuclease, ligase, polymerase, and sequencing reactions
US20050009059A1 (en) Analysis of methylation status using oligonucleotide arrays
US9828640B2 (en) Analysis of methylation using nucleic acid arrays
US20060292585A1 (en) Analysis of methylation using nucleic acid arrays
US20140378340A1 (en) Methods for Genotyping
US20050196792A1 (en) Analysis of methylation status using nucleic acid arrays
EP1723256B1 (en) CpG-AMPLICON AND ARRAY PROTOCOL
US20080254453A1 (en) Analysis of methylation using selective adaptor ligation
WO2023107682A1 (en) Methods for the detection of a nucleic acid

Legal Events

Date Code Title Description
AS Assignment

Owner name: AFFYMETRIX, INC., CALIFORNIA

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:NAUTIYAL, SHIVANI;FAHAM, MALEK;REEL/FRAME:020356/0928;SIGNING DATES FROM 20071217 TO 20080103

Owner name: AFFYMETRIX, INC., CALIFORNIA

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:NAUTIYAL, SHIVANI;FAHAM, MALEK;SIGNING DATES FROM 20071217 TO 20080103;REEL/FRAME:020356/0928

AS Assignment

Owner name: GENERAL ELECTRIC CAPITAL CORPORATION, AS AGENT, MA

Free format text: SECURITY AGREEMENT;ASSIGNOR:AFFYMETRIX, INC.;REEL/FRAME:028465/0541

Effective date: 20120625

AS Assignment

Owner name: AFFYMETRIX, INC., CALIFORNIA

Free format text: RELEASE BY SECURED PARTY;ASSIGNOR:GENERAL ELECTRIC CAPITAL CORPORATION, AS AGENT;REEL/FRAME:037109/0132

Effective date: 20151028

STCB Information on status: application discontinuation

Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION