WO2005098035A2

WO2005098035A2 - Method for the quantification of methylated dna

Info

Publication number: WO2005098035A2
Application number: PCT/EP2005/003793
Authority: WO
Inventors: David GÜTIG; Dirk Habighorst; Antje Kluth; Armin Schmitt; Matthias Schuster; Ina Schwope
Original assignee: Epigenomics Ag
Priority date: 2004-04-06
Filing date: 2005-04-06
Publication date: 2005-10-20
Also published as: EP1733054A2; US20050287553A1; CA2559426A1; AU2005231971A1; WO2005098035A3

Abstract

The method according to the invention concerns a method for the quantification of two different variations of a DNA sequence. Particularly, the invention relates to a quantification of methylated DNA. For this purpose, the DNA to be investigated is first converted so that cytosine is converted to uracil, while 5-methylcytosine remains unchanged. Then the converted DNA is amplified by means of a real-time PCR. However, probes are utilized, one of which is specific for the methylated status and one for the unmethylated status of the DNA. The degree of methylation of the investigated DNA can be calculated from the ratio of the signal intensities of the probes or from the Ct values. The method according to the invention is particularly suitable for the diagnosis and prognosis of cancer and other disorders associated with a change in the methylation status, as well as for the prediction of drug interactions.

Description

Method for the quantification of methylated DNA

The present invention concerns a method for the quantification of methylated cytosine positions in DNA. Background of the invention

5-Methylcytosine is the most frequent covalently modified base in the DNA of eukaryotic cells. It plays an important biological role, among other things, in the regulation of transcription, in genetic imprinting and in tumorigenesis (for review: Millar et al.: Five not four: History and significance of the fifth base. In: The Epigenome, S. Beck and A. Olek (eds.), Wiley-VCH Publishers, Weinheim 2003, pp. 3-20). The identification of

5-methylcytosine is particularly of considerable interest for cancer diagnosis. It is difficult to detect methylcytosine, of course, since cytosine and 5-methylcytosine have the same base-pairing behavior. The conventional DNA analysis methods based on hybridization are thus not applicable. Correspondingly, current methods for methylation analysis operate according to two different principles. In the first one, methylation-specific restriction enzymes are utilized, and in the second one, a selective chemical conversion of unmethylated cytosines to uracil is conducted (so-called bisulfite treatment; see e.g., PCT/EP2004/011715). The enzymatically or chemically pretreated DNA is then amplified for the most part and can be analyzed in different ways (for review: WO 02/072880 pp. 1 ff); Fraga and Estella: DNA methylation: a profile of methods and applications. Biotechniques. 2002 Sep;33(3):632, 634, 636-49). For a sensitive analysis, the chemically pretreated DNA is usually amplified by means of a PCR method. A selective amplification only of methylated (or with the reverse approach, unmethylated) DNA can be assured by employing methylation-specific primers or blockers (so-called methylation-sensitive PCR/MSP or "Heavy Methyl" methods; see: Herman et al.:

Methylation-specific PCR: a novel PCR assay for methylation status of CpG islands. Proc Natl Acad Sci U S A. 1996 Sep 3;93(18):9821-6 Cottrell et al.: A real-time PCR assay for DNA-methylation using methylation-specific blockers. Nucl. Acids Res. 2004 32: e10). On the other hand, it is also possible to amplify the DNA first in a non-methylation-specific manner and then to analyze the amplificates by means of methylation-specific probes (for review: Trinh et al.: DNA methylation analysis by MethyLight technology. Methods. 2001 Dec;25(4):456-62). The named PCR methods are also applicable as real-time PCR variants. These make it possible to detect the methylation status directly in the course of the PCR without the necessity of a subsequent analysis of the products ("Methyϋght" - WO 00/70090; US 6,331 ,393; Trinh et al. 2001 , loc. cit.)

A quantification of the degree of methylation is necessary for different applications, e.g., for classifications of tumors, for prognostic information or for the prediction of drug effects. Different methods are known for the quantification of the degree of methylation. First, an amplification of the DNA is produced, in part, e.g., with Ms-SNuPE, with hybridizations on microarrays, with hybridization assays in solution or with direct bisulfite sequencing (for review: Fraga and Estella 2002, loc. cit.). A problem with these "end point analyses" consists of the fact that the amplification can occur non uniformly, among other things, due to obstruction of product, enzyme instability and a decrease in concentration of the reaction components. A correlation between the quantity of amplificate and the quantity of DNA utilized is thus not always suitable. Quantification is thus sensitive to error (see: Kains: The PCR plateau phase - towards an understanding of its limitations. Biochem. Biophys. Acta 1494 (2000) 23-27). Threshold-value analysis, which is based on a real-time PCR determines the quantity of amplificate, in contrast, not at the end of the amplification, but in the exponential phase of the amplification. This method presumes that the amplification efficiency is constant in the exponential phase. The so-called threshold value Ct is a measure for that PCR cycle, in which the signal in the exponential phase of the amplification is greater for the first time than the background signal. Absolute quantification then results by means of a comparison of the Ct value of the investigated DNA with the Ct value of a standard (see: Trinh et al. 2001 , loc. cit.; Lehmann et al.: Quantitative assessment of promoter hypermethylation during breast cancer development. Am J Pathol. 2002 Feb;160(2):605-12). A problem of Ct analysis consists of the fact that with high DNA concentrations only a small resolution can be achieved. This applies also when high degrees of methylation are to be determined via PMR values (for PMR values, see: Eads et al., CANCER RESEARCH 61 , 3410-3418, April 15, 2001.) The amplification of a reference gene, e.g., the β-actin gene, is also necessary for this type of Ct analysis (see: Trinh et al. 2001 , loc cit.). In the following, a new real-time PCR method is described for quantitative methylation analysis. Here, a non-methylation-specific, conversion-specific amplification of the target DNA is produced. The amplificates are detected by means of the hybridization of two different methylation-specific real-time PCR probes. Here, one of the probes is specific for the methylated state, while the other probe is specific for the unmethylated state. The two probes bear different fluorescent dyes. A quantification of the degree of methylation can be produced within specific PCR cycles employing the ratio of signal intensities of the two probes. Alternatively, the Ct values of two fluorescent channels can also be drawn on for the quantification of the methylation. In both cases, a quantification of the degree of methylation is possible without the necessity of determining the absolute DNA quantity. A simultaneous amplification of a reference gene or a determination of the PMR values is thus not necessary. In addition, the method according to the invention supplies reliable values for both large and small DNA quantities as well as for high and low degrees of methylation.

Based on the special importance of cytosine methylation, there is a great technical need for quantitative methods of methylation analysis. The method according to the invention overcomes the above-described disadvantages of the prior art and thus represents an important technical advance.

Description

The method according to the invention for the quantification of methylated DNA is characterized in that the following steps are conducted: a) the DNA to be investigated is reacted in such a way that 5-methylcytosine remains unchanged, while unmethylated cytosine is converted into uracil or into another base which is distinguished from cytosine in its base-pairing behavior; b) the converted DNA is amplified in the presence of two real-time probes, wherein one of the probes is specific for the methylated state, and the other probe is specific for the unmethylated state of the DNA; c) it is determined at different time points how far the amplification has proceeded by detecting the hybridization of the probes to the amplificates, d) the degree of methylation of the investigated DNA is determined.

In the following the present invention is also referred to as "QM" (quantitative methylation) assay.

In the first step of this embodiment, the DNA to be investigated is reacted with a chemical or with an enzyme in such a way that 5-methylcytosine remains unchanged, while unmethylated cytosine is converted into uracil or into another base which is distinguished from cytosine in its base-pairing behavior. In this case, the DNA to be investigated can originate from different sources, each time depending on the diagnostic or scientific objective. For diagnostic objectives, tissue samples are preferably used as the initial material, but body fluids, particularly serum, can also be used. It is also possible to use DNA from sputum, stool, urine, or cerebrospinal fluid. Preferably, the DNA is first isolated from the biological sample. The DNA is extracted according to standard methods from blood, e.g., with the use of the Qiagen UltraSens DNA extraction kit. The isolated DNA can then be fragmented, e.g., by reaction with restriction enzymes. The reaction conditions and the enzymes employed are known to the person skilled in the art and result, e.g., from the protocols supplied by the manufacturers. Then the DNA is converted chemically or by means of enzymes. A chemical conversion by means of bisulfite is preferably conducted. The bisulfite conversion is known to the person skilled in the art in different variants (see, e.g.: Frommer et al.: A genomic sequencing protocol that yields a positive display of 5-methylcytosine residues in individual DNA strands. Proc Natl Acad Sci U S A. 1992 Mar 1 ; 89(5): 1827-31 ; Olek, A modified and improved method for bisulphite based cytosine methylation analysis. Nucleic Acids Res. 1996 Dec 15;24(24):5064-6; PCT/EP2004/011715). It is particularly preferred that the bisulfite conversion is conducted in the presence of denaturing solvents, e.g., dioxane and a radical trap (see: PCT/EP2004/011715). In another preferred embodiment, the DNA is not chemically converted, but rather is converted by enzymes. This is conceivable, e.g., by the use of cytidine deaminases, which convert unmethylated cytidine more rapidly than methylated cytidine. An appropriate enzyme has recently been identified (Bransteitter et al.: Activation-induced cytidine deaminase deaminates deoxycytidine on single-stranded DNA but requires the action of RNase. Proc Natl. Acad Sci U S A. 2003 Apr 1 ;100(7):4102-7).

In the second step of the method according to the invention, the converted DNA is amplified in the presence of two real-time probes, wherein one of the probes is specific for the methylated state, and the other probe is specific for the unmethylated state. Preferably, an amplification is conducted by means of an exponential amplification process, particularly preferably by means of a PCR. Primers are used for the amplification, which are specific for the chemically or enzymatically converted DNA. In this case, non-methylation-specific primers are preferably utilized, i.e., primers which do not make available CG or methylation-specific TG or CA dinucleotide. A uniform amplification of methylated and unmethylated DNA is conducted with these primers. It is also possible, however, to amplify a larger sequence region in a methylation-specific manner and thus to quantify specific cytosine positions within this sequence by means of the method according to the invention. The design of methylation-specific and non-methylation-specific primers and the PCR reaction conditions belong to the prior art (see: e.g.: US Patent 6,331,393; Trinh et al 2001 , loc. cit). The primers are preferably located close to the probe. The length of the amplicon should not exceed 200 bp. The melting temperature ™ should be from 52 to 60 °C (depending on probe-Tm, approx. 5-7°C below probe-Tm).

The amplification is conducted in the presence of two different probes, wherein one of the probes is specific for the methylated state of the DNA , while the other probe is specific for the unmethylated state of the DNA. The methylation-specific probes correspondingly bear at least one CpG dinucleotide, while the non-methylation-specific probes make available at least one specific TG or CA dinucleotide. Preferably the probes bear three specific dinucleotides. Both probes cover the same CpG-positions. Melting temperatures of the probes should be similar. Moreover, the probes should cover positions representing converted C-positions in order to ensure conversion-specific detection. The probes involve real-time probes. These real-time probes are understood in the following to be probes which permit the amplificates to be detected during the amplification. Different real-time PCR variants are familiar to the person skilled in the art, e.g., Lightcycler, Taqman, Sunrise, Molecular Beacon or Eclipse probes. The particulars on constructing and detecting these probes belong to the prior art (see: US Patent 6,331 ,393 with additional citations). Thus the design of the probes is carried out manually or by means of the "PrimerExpress" software of Applied Biosystems (for Taqman probes) or via the MGB Eclipse design software of Epoch Biosciences (for Eclipse probes). Preferably, Taqman probes are used, which are utilized most preferably in combination with Minor Groove Binders (MGB).

Taqman probe design preferably follows the design guidelines given by Applied Biosystems for the "Taqman Allelic Discriminiation" assay . So both probes preferably have the same 5'-end, which has impact on the 5'-exonuclease activity of the polymerase. Runs of identical nucleotides (> 4 bases, esp. G) should be avoided. Preferably, there is no G at 5'-end (quenching). The probes should contain more Cs than Gs and the polymorphic site should preferably be located approx. in the middle third of the sequence. Preferred reporter dyes are FAM and VIC.

The amplification is preferably conducted together with both probes in one vessel, so that the reaction conditions for both probes are identical. This embodiment also leads to an increased specificity, since the probes compete for binding sites. It is necessary, of course, that the two probes bear different labels. On the other hand, it is also possible to conduct the amplifications in different vessels. In this way, disruptive interactions between the fluorescent dyes can be avoided. When performing amplifications and detection with 2 probes in 2 vessels a competing, unlabeled oligonucleotide can be used in order to increase specificity of probe binding.

In the third step of the method according to the invention, it is determined at different time points how far the amplification has proceeded. This is done by detecting hybridizations during the individual amplification cycles. Depending on the probes utilized, detection is made according to the prior art.

In the fourth step of the method according to the invention, the degree of methylation of the investigated DNA is determined. This can be done by means of different embodiments. In a preferred embodiment, the degree of methylation of the investigated DNA is determined from the ratio of the signal intensities of the two probes. This can be accomplished by means of the following formula:

The notation I_CG indicates the signal intensity of the probe specific for the methylated state and ITG indicates the signal intensity of the probe specific for the unmethylated state.

The signal intensities during a PCR cycle in the exponential amplification phase of the PCR are particularly preferably placed in a ratio to one another. A calculation is preferably carried out close to the cycle, in which the amplification reaches its maximal increase. This corresponds to the point of inflection of the fluorescent intensity curve or the maximum of its first derivative.

The calculation is thus conducted at a time point which preferably lies at up to five cycles before or after the inflection point, particularly preferably up to two cycles before or after the inflection point, and most particularly preferred up to one cycle before or after the inflection point. In the optimal embodiment, the calculation occurs directly at the inflection point.

For the case when the inflection points of the two curves lie in different cycles, the calculation is preferably conducted at the inflection point of the curve which has the highest signal at this time point.

The determination of the inflection point is preferably made by means of the first derivative of the fluorescent intensity curves. The derivatives are preferably first subjected to a smoothing ("Spline", see: Press, W. H., Teukolsky, S. A., Vetterling, W. T., Flannery, B. P. (2002). Numerical Recipes in C. Cambridge: University Press; Chapter 3.3.).

In another preferred embodiment, the calculation of the degree of methylation is conducted not by means of the ratio of the fluorescent intensities, but by means of the ratio of threshold values at which a certain signal intensity will be exceeded, e.g., at the Ct values (see above). The determination of Ct values is found in the prior art (see: Trinh et al., loc.cit, 2002). The degree of methylation can then be determined via the following formula: degree of methylation = 100/ (1+2^ΔCt).

In addition, it is conceivable to use other criteria for calculating the degree of methylation, e.g., the area under the fluorescent curve (area under the curve) or the maximal slope of the curves or the maximum of the second derivative of amplification.

A quantification by means of the above-described method is very well possible if the assay conditions have been previously optimized in this respect. An optimization is conducted with different methylation standards (e.g., with 0%, 5%, 10%, 25%, 50%, 75% and 100% degree of methylation). DNA, which covers the entire genomic DNA or a representative portion thereof, is preferably used as the standard. The different degrees of methylation are obtained by appropriate mixtures of methylated and unmethylated DNA. The production of methylated DNA is relatively simple with the use of Sssl methylase. This enzyme converts all unmethylated cytosines in the sequence context CG to 5-methylcytosine. Sperm DNA, which provides only a small degree of methylation, can be used as completely unmethylated DNA (see: Trinh et al. 2001 , loc.cit.). The preparation of unmethylated DNA is preferably conducted by means of a so-called genome-wide amplification (WGA - whole genome amplification; for review: Hawkins et al.: Whole genome amplification-applications and advances. Curr Opin Biotechnol. 2002 Feb; 13(1): 65-7) WGA) . Here, wide parts of the genome will be amplified by means of "random" or degenerate primers. Since only unmethylated cytosine nucleotides will be provided in the amplification, a completely unmethylated DNA results after several amplification cycles. Thus a "Multiple Displacement Amplification" by means of the φ29 polymerase is preferably produced (MDA, see: Dean et al. 2002 loc.cit.; US Patent 6,124,120). Correspondingly produced DNA is available from different commercial suppliers ("GenomiPhi" of Amersham Biosciences, www4.amershambiosciences.com; "Repli-g" of Molecular Staging, www.moIecularstaging.com). The production of methylation standards is described in great detail in the European Patent Application 04 090 037.5 (date filed: Feb. 5, 2004; Applicant: Epigenomics AG). By calculating the quotient of the signals which are detected for the methylated state and the sum of the signals which are detected for the methylated and the unmethylated state, the measured methylation rate is obtained. If this is plotted against the theoretical methylation rates (corresponding to the proportion of methylated DNA in the defined mixtures) and the regression line which passes through the measured points is determined, a calibration curve is obtained. A calibration is conducted preferably with different quantities of DNA, e.g., with 0.1 , 1 and 10 ng of DNA per batch.

Assays are particularly suitable for quantification by means of the method according to the invention, if the calibration curves for the time point of the exponential amplification provide a y-axis crossing as close as possible to zero. Methylation states that are adjacent should be distinguished by a high Fisher score (preferably greater than 1 , most preferably greater than 3).

In addition, it is advantageous if a y-axis intercept is provided that is as small as possible and a Fisher score is provided that is as high as possible (preferably greater than 1 , particularly preferred greater than 3). In addition, it is advantageous if the curves have a slope and a regression close to the value 1.

The assays can be optimized in this respect by means of varying the primers, the probes, the temperature program and the other reaction parameters using standard tests.

As has already been mentioned above, the methylation rate can be determined with the method according to the invention, independently from a standard curve. If a standard curve is prepared, however, then the absolute content of methylated DNA can be determined also very simply by means of the method according to the invention. A particularly preferred use of the method according to the invention lies in the diagnosis or prognosis of cancer diseases or other disorders associated with a change of the methylation status. These include, among others: CNS malfunctions; symptoms of aggression or behavioral disturbances; clinical, psychological and social consequences of brain damage; psychotic disturbances and personality disorders; dementia and/or associated syndromes; cardiovascular disease, malfunction and damage; malfunction, damage or disease of the gastrointestinal tract; malfunction, damage or disease of the respiratory system; lesion, inflammation, infection, immunity and/or convalescence; malfunction, damage or disease of the body as a consequence of an abnormality in the development process; malfunction, damage or disorder of the skin, the muscles, the connective tissue or the bones; endocrine and metabolic malfunction, damage or disease; headaches or sexual malfunction. The method according to the invention is also suitable for predicting undesired drug interactions and for the differentiation of cell types or tissues or for the investigation of cell differentiation.

A kit is also included according to the invention, and this comprises two primers, a polymerase, as well as a probe specific for the methylated state and a probe specific for the unmethylated state as well as, optionally, additional reagents necessary for a PCR and/or a bisulfite reagent.

It is known to the person skilled in the art that all of the above-named embodiments of the method according to the invention can be used not only for the methylation analysis, but also for the quantification of sequence differences in RNA or in DNA. For this, the first step of the described method — the chemical or enzymatic conversion — is not conducted. Thus it is possible to investigate an allele-specific gene expression or a gene duplication by means of the method according to the invention. In addition, it is possible to investigate single nucleotide polymorphisms (SNPs) from the pooled samples. Another application is the quantification of different strains of microorganisms. Therefore a general principle of the present invention is a quantification of two different variation of a DNA sequence, characterized in that the following steps are conducted: a) the DNA is amplified in the presence of two real-time probes, whereby one of the probes is specific for one variation of the DNA sequence, and the other probe is specific for the other variation of the DNA sequence; b) it is determined at different time points how far the amplification has proceeded by detecting the hybridization of the probes to the amplificates, c) the proportions of the two sequence variations is determined.

All of the above-described preferred embodiments can be transferred correspondingly to these applications outside the sphere of methylation analysis. The resulting applications are thus also a part of this invention. The person skilled in the art knows how he must modify the above-described methods.

The same applies to the indicated uses and the kit.

Therefore a part of the present invention is also a kit comprising two primers, a polymerase, a probe specific for one variation of the DNA sequence and a probe specific for the other variation of the DNA sequence to be investigated. The kit may optionally contain additional reagents necessary for a PCR.

One application of the method according to the invention is the investigation of allele-specific gene expression (see for allele-specific gene expression: Lo et al.: Allelic variation in gene expression is common in the human genome. Genome Res. 2003 Aug;13(8):1855-62; Weber et al.: A real-time polymerase chain reaction assay for quantification of allele ratios and correction of amplification bias. Anal Biochem 2003 Sep 15;320(2):252-8). First, a reverse transcription is necessary for the application of the method according to the invention. In particular, a method for the quantification of allele-specific gene expression is offered according to the invention, which is characterized in that: a) the RNA to be investigated is reverse-transcribed, b) the cDNA is amplified in the presence of two real-rime probes, whereby one of the probes is specific for one of the alleles and the other probe is specific for the other allele, c) at different time points it is determined how far the amplification has proceeded by detecting the hybridizations of the probes to the amplificates, d) the allele-specific gene expression is quantified.

In the first step of this embodiment, the RNA to be investigated is reverse-transcribed. Appropriate methods are found in the prior art (see: Lo et al. 2003, loc.cit.). Usually the RNA is isolated first. Various commercially available kits can be used for this purpose (e.g. Micro-Fast Track, Invitrogen; RNAzol B, Tel-Test). The cDNA is then produced by means of a commercially available reverse transcriptase (e.g., from Invitrogen).

In the second step of the method according to the invention, the cDNA is amplified in the presence of two real-time probes, whereby one of the probes is specific for the sequence of one allele, and the other probe is specific for the sequence of the other allele. The probes involve real-time probes, thus, e.g., Lightcycler, Taqman, Sunrise, Molecular Beacon or Eclipse probes. The particulars on constructing and detecting these probes belong to the prior art (see above).

Preferably, the amplification is conducted by means of an exponential amplification process, most preferably by means of a PCR. Primers are used for the amplification, which amplify the DNA of both alleles in a uniform manner. The design of primers and probes as well as the PCR reaction conditions belong to the prior art (see above). The amplification is preferably conducted together with both probes in one vessel, so that the reaction conditions for both probes are identical (see above).

In the third step of the method according to the invention, it is determined at different time points how far the amplification has proceeded. This is done by detecting the hybridizations of the probes to the amplificates during the individual amplification cycles Depending on the probes utilized, detection is made according to the prior art (see above).

In the fourth step of the method according to the invention, the allele-specific gene expression is quantified. This can be done — as is described above in detail for the methylation analysis — by means of different embodiments. In a preferred embodiment, quantification is made by means of the ratio of signal intensities of the two probes. However, it is also possible to draw on the area under the fluorescent curves or the maximal slope of the curves for quantifying the ratio of the threshold values (see above).

As was described in detail for the methylation analysis, a quantification is very well possible if the assay conditions have been previously optimized in this respect. For this purpose, a calibration curve is plotted by means of a standard series which contains different proportions of the two allele sequences. The quality criteria (y-axis intercept, Fisher score, slope regression) described in detail for the methylation analysis also apply here in a general sense.

Another application of the method according to the invention outside the sphere of methylation analysis is the investigation of single nucleotide polymorphisms (SNPs) from pooled samples. A pool of samples is meaningful for different objectives, e.g., for identifying genes which take part in the emergence of complex disorders (see: Shifman et al.: Quantitative technologies for allele frequency estimation of SNPs in DNA pools. Mol Cell Probes 2002 Dec;16(6):429-34). Accordingly, a method for investigating SNPs from pooled samples is included in the invention, which is characterized in that: a) the sample to be investigated is amplified in the presence of two real-time probes, whereby one of the probes is specific for the sequence of one SNP, and the other probe is specific for the sequence of the other SNP, b) at different time points it is determined how far the amplification has proceeded, by detecting the hybridizations of the probed to the amplificates , c) it is concluded from this which SNP at what fraction is represented in the pool. A gene duplication can also be investigated according to the same principle (see: Pielberg et al.: A sensitive method for detecting variation in copy numbers of duplicated genes. Genome Res 2003 Sep;13(9):2171-7).

Another application of the method according to the invention is the investigation of mutations in microorganisms. Thus the proportion of wild type and the proportion of mutant strain can be determined in a sample. Such an application can be of importance for therapeutic decisions (see, e.g.: Nelson et al.: Detection of all single-base mismatches in solution by chemiluminescence. Nucleic Acids Res 1996 Dec 15;24(24):4998-5003). This embodiment of the method according to the invention is accordingly characterized in that a) the sample to be investigated is amplified in the presence of two real-time probes, whereby one of the probes is specific for the sequence of the wild type, and the other probe is specific for the sequence of the mutant strain, d) at different time points it is determined how far the amplification has proceeded, by detecting the hybridizations of the probes to the amplificates e) it is concluded from this which strain at that fraction is represented in the sample.

Examples Example 1

It shall be pointed out that a reliable quantification of DNA methylation is possible with the method according to the invention. For this purpose, the degree of methylation of the two genes S100A2 and TFF1 will be analyzed. Calibration curves with several DNA mixtures of different degrees of methylation were plotted. A series of DNA mixtures of known degrees of methylation were used as the standard (0, 5, 10, 25, 50, 75 and 100% methylated DNA). For the production of this "gold standard", completely methylated and completely unmethylated DNA were mixed together in different ratios. The completely unmethylated DNA was obtained from Molecular Staging. It was prepared there by means of a multiple displacement amplification of human genomic DNA from whole blood. The completely methylated DNA was produced by means of an Sssl treatment of the completely unmethylated DNA according to the manufacturer's instructions. The DNA was then bisulfite-converted (see PCT/EP2004/011715). For the real-time PCR assays, primer pairs were used which were specific for the bisulfite conversion. The primers, however, were nonspecific for methylation, i.e., they did not contain CpG positions. Two bisulfite-specific MGB-Taqman probes (Applied Biosystems) were also utilized. These probes comprised 2 CpG positions. One probe was specific for the methylated state and was labeled with FAM. The second probe was specific for the unmethylated state and bore a VIC label (see Fig. 1). The following primers and probes were used for TFF1 : methylation-specific probe: 6FAM-ACACCGTTCGTaaaa- MGBNFQ (Seq ID1), non-methylation-specific probe VIC-ACACCATTCATaaaaT-MGBNFQ (Seq ID 2), Forward Primer: AGtTGGTGATGtTGATtAGAGtt (Seq ID 3), Reverse Primer CCCTCCCAaTaTaCAAATAAaaaCTa (Seq ID 4). The following oligonucleotides were utilized for S100A2: methylation-specific probe: 6FAM- tTCGTGTAtATAtATGCGttTG- MGBNFQ (Seq ID 5), non-methylation-specific probe VIC- tTTGTGTAtATAtATGTGttTGTG-MGBNFQ (Seq ID 6), Forward Primer TttTGTGTGAGAGGtTGTGAGtAt (Seq ID 7), Reverse Primer CCTCCTaATaTCCCCCAaCT (Seq ID 8). The real-time PCR was carried out in an

ABI7700 Sequence Detection System (Applied Biosystems) in a 20 μ\ reaction volume. The final concentrations in the reaction mixtures amounted to: IxTaqMan Buffer A (Applied Biosystems) containing ROX as a passive reference dye, 2.5 mmol/l MgCI₂ (Applied Biosystems), 1 U of AmpliTaq Gold DNA polymerase (Applied Biosystems), 625 nmol/l primers, 200 nmol/l probes, 200 μmol/l dNTPs. The temperature profile for the TFF1 assay was conducted as follows: 10 min activation at 94°C, followed by 45 cycles of 15 s at 94°C denaturing and 60 s at 60°C annealing + elongation. The fluorescence was measured during the 60°C step (Fig. 2). The annealing was conducted at 62°C for the S100A2 assay. The data analysis was conducted according to the recommendations of Applied Biosystems. The degrees of methylation were determined according to the following formula: methylation rate: = delta Rn CG probe/ (delta Rn CG probe + delta Rn TG probe). By plotting the measured methylation rates against the theoretical methylation rates, a calibration curve was prepared for each PCR cycle (Fig. 3). The suitability of the individual curves for the quantification was determined by means of the following curve parameters: slope, R², y-axis intercept as well as Fisher scores for the classification of adjacent methylation levels each time (Fig. 3). From the same experiments, calibration curves were plotted on the threshold cycles (Ct), wherein the methylation rate was calculated with the following formula: methylation rate = 100/ (1 + 2 ^{delta C} ) (Fig. 4). If the suitability of the different cycles (optimal cycle, in which the slope of the amplification curve is maximal, vs. final cycle) is compared with the suitability of the calibration based on Ct values, it can be seen that overall the calibration by means of the optimal cycle produces the best curve parameters (Fig. 5): slope close to 1 , R² close to 1 , y-axis intercept close to 0, Fisher scores >1.

Example 2 It will be shown that the method according to the invention makes possible a reliable quantification of the methylation of different types of samples. For this purpose, a portion of the biological sample material was fresh frozen, and the remainder was embedded in paraffin. Then the DNA was isolated from the sample first according to the standard techniques and after this, it was bisulfited (see, e.g. German Patent Application 10347400.5). Then the DNA was amplified by means of two non-methylation-specific primers in the presence of two Taqman oligonucleotide probes. One of the oligonucleotide probes was specific for the methylated state, and the other for the unmethylated state of the investigated gene. Both probes had a reporter fluorescent dye at the 5'-end and a quencher at the 3'-end. The reactions were calibrated with DNA standards of a defined methylation status as described above. The β-actin gene (ACTB) was investigated for determining the quantity of sample DNA. The primers and probes utilized here did not provide CpG dinucleotides, so that the amplification was produced here independently of the methylation status. Thus only one probe was necessary here. The following oligonucleotides were used: Primer 1 : TGGTGATGGAGGAGGTTTAGTAAGT (SEQ ID NO: 9); Primer 2: AACCAATAAAACCTACTCCTCCCTTAA (SEQ ID NO: 10); probe: 6FAM-ACCACCACCCAACACACAATAACAAACACA-TAMRA or Dabcyl (SEQ ID NO. 1 ). The following reaction components were utilized: 3 mmol/l MgCl buffer, 10x buffer, Hotstart TAQ. The following temperature program was used: 95°C for 10 minutes, then 45 cycles: 95 °C, 15 sec; 62 °C, 1 min. The fluorescent signals were recorded with a Lightcycler device. The degree of methylation of a specific locus was determined by the following formula: degree of methylation = 100 * I CG / (ICG + ITG) (I = fluorescent intensity of the CG or TG probe) Table 1 shows the results of Example 2. "Fresh" denotes fresh frozen tissue, "PET" stands for paraffin-embedded tissue. In all, 18 sample pairs were investigated. It was shown that the method according to the invention makes possible a quantification from both types of samples.

Table 1

Example 3 Reliability of the QM assay within a broad range of input DNA. Different amounts of bisulfite DNA ( 50, 10, 5, 1 ng) derived from nine different samples (fresh frozen tissue samples and paraffin embedded tissue samples) were analysed by the QM assay. The results are illustrated in Figure 6. It is shown that the QM assays perform well in a wide range of input DNA. The determined methylation degree is independent of the DNA input amount. The standard deviation does not exceed a value of ± 5 percentage points around the mean of measured methylation rate. This value of the standard deviation is caused by the interplate variability (see Example 4) . Example 4

Reproducibility of QM assay. In order to investigate the reproducibility of the QM assay 12 different QM assays were conducted in five separate runs. As indicated in Fig. 7, the assays showed a low intra- and inter-plate variability. The confidence interval is around ± 5 percentage points of the mean of the methylation rate (Fig. 8). Example 5 In order to provide a comparison of methylation analysis by means of array ("chip") analysis to the assay of the present invention, methylation of the gene PITX2 was analysed in patients with breast cancer.

The following study was based on samples from 236 breast cancer patients, wherein all patients were NO (nodal status negative), and older than 35 years. In all cases surgery was performed before 1998. All patients were ER+ (estrogen receptor positive), and the tumors were graded to be T1-3, G1-3. In this study all patients received Tamoxifen directly after surgery, and the outcome was assessed according to the length of disease-free survival.

The DNA samples were extracted using the Wizzard Kit (Promega). Total genomic DNA of all samples was bisulfite treated converting unmethylated cytosines to uracil.

Methylated cytosines remained conserved. Bisulfite treatment was performed with minor modifications according to the protocol described in Olek et al. (1996). After bisulfitation 10 ng of each DNA sample was used in subsequent mPCR reactions containing 6-8 primer pairs. Each reaction contained the following: 2.5 pmol each primer; 11.25 ng DNA (bisulfite treated); Multiplex PCR Master mix (Qiagen); The primer oligonucleotides used to generate the amplificate, were: GTAGGGGAGGGAAGTAGATGT (SEQ ID NO: 12); TCCTCAACTCTACAAACCTAAAA (SEQ ID NO: 13). Initial denaturation was carried out at 95°C for 15 min. Forty cycles were carried out as follows: Denaturation at 95°C for 30 sec, followed by annealing at 57°C for 90 sec, primer elongation at 72°C for 90 sec. A final elongation at 72°C was carried out for 10 min. All PCR products from each individual sample were then hybridised to glass slides carrying a pair of immobilised oligonucleotides for each CpG position under analysis. Each of these detection oligonucleotides was designed to hybridise to the bisulphite converted sequence around one CpG site which was either originally unmethylated (TG) or methylated (CG). Hybridisation conditions were selected to allow the detection of the single nucleotide differences between the TG and CG variants. 5 μl volume of each multiplex PCR product was diluted in 10 x Ssarc buffer . The reaction mixture was then hybridised to the detection oligonucleotides as follows. Denaturation at 95°C, cooling down to 10 °C, hybridisation at 42°C overnight followed by washing with 10 x Ssarc and dH2O at 42°C. The sequences of the oligonucleotides used were the following: AGTCGGGAGAGCGAAA (SEQ ID NO 14); AGTTGGGAGAGTGAAA (SEQ ID NO 15). Fluorescent signals from each hybridised oligonucleotide were detected using genepix scanner and software. Ratios for the two signals (from the CG oligonucleotide and the TG oligonucleotide used to analyse each CpG position) were calculated based on comparison of intensity of the fluorescent signals.

The log methylation ratio (log(CG/TG)) at each CpG position is determined according to a standardised pre-processing pipeline that includes the following steps: For each spot the median background pixel intensity is subtracted from the median foreground pixel intensity (this gives a good estimate of background corrected hybridisation intensities): For both CG and TG detection oligonucleotides of each CpG position the background corrected median of 4 redundant spot intensities is taken; For each chip and each CpG position the log(CG/TG) ratio is calculated; For each sample the median of log(CG/TG) intensities over the redundant chip repetitions is taken. This ratio has the property that the hybridisation noise has approximately constant variance over the full range of possible methylation rates (Huber et al., 2002).

The same samples were then analysed by means of the assay of the present invention.

The amount of sample DNA amplified was quantified by reference to the gene (β-actin

(ACTB)) to normalize for input DNA. For standardization the primers and the probe for analysis of the ACTB gene lacked CpG dinucleotides so that amplification is possible regardless of methylation levels. As there are no methylation variable positions, only one probe oligonucleotide is required.

The following oligonucleotides were used in the reaction to amplify the control amplificate: Control Primerl : TGGTGATGGAGGAGGTTTAGTAAGT (SEQ ID NO: 16) Control Primer2: AACCAATAAAACCTACTCCTCCCTTAA (SEQ ID NO: 17) Control Probe: 6FAM-ACCACCACCCAACACACAATAACAAACACA-TAMRA or Dabcyl (SEQ ID NO: 18); The following primers are used to generate an amplificate within the PITX2 sequence comprising the CpG sites of interest: Primers for PITX bisulfite amplificate length: 144 bp PITX2: GTAGGGGAGGGAAGTAGATGTT (SEQ ID NO: 19); PITX2: TTCTAATCCTCCTTTCCACAATAA (SEQ ID NO: 20); Probes: PITX2cg1 : FAM-AGTCGGAGTCGGGAGAGCGA-Darquencher (SEQ ID NO: 21); As an alternative quencher TAMRA was also used in additional experiments: FAM-AGTCGGAGTCGGGAGAGCGA-TAMRA; PITX2tg1 : YAKIMA YELLOW-AG TTGGAGTTGGGAGAGTGAAAGGAGA-Darquencher (SEQ ID NO: 22).

The extent of methylation at a specific locus was determined by the following formula:

methylation rate = 100 ^* I (CG) / (l(CG) + l(TG)), (I = Intensity of the fluorescence of CG-probe or TG-probe).

PCR components were ordered from Eurogentec: 3 mM MgCI2 buffer, 10x buffer, Hotstart TAQ; Program (45 cycles): 95 °C, 10 min; 95 °C, 15 sec; 62 °C, 1 min

Results

For each assay the methylation (and where relevant mean methylation over multiple oligo-pairs) for each amplificate was calculated and the population split into groups according to their mean methylation values, wherein one group was composed of individuals with a methylation score higher than the median and a second group composed of individuals with a methylation score lower than the median. Results are shown in figures 9 to 11. The survival curves generated by microarray analysis could be well confirmed by the new QM assay (fig 9 and 10). The correlation plot between microarray and QM assay is shown in fig 11 , indicating a co-efficient of 0.87. Therefore, methylation markers pre-validated by microarray methylation analysis are well transferable to the QM-assay format.

Brief description of the figures

Figure 1 shows the principle of the QM assay. Primers are used for the amplification, which are bisulfite-specific, but contain no CpG positions (shown as black circles); the probes are specific for the methylated or the unmethylated state of the covered CpG positions; if one uses both probes in the same reaction, then they are labeled with different fluorescent dyes (R1 , R2; Q = quencher).

Figure 2 shows the results of Example 1. Represented is the detection of the amplification products of TFF1 in each cycle (x-axis) by means of fluorescent signals of the hybridized probes (y-axis: fluorescent intensity); A: Amplification curves of DNA mixtures of known methylation levels detected with the FAM-labeled probe for the methylated state; B: corresponding detection with the VIC-labeled probe for the unmethylated state.

Figure 3 shows the results of Example 1. Represented are the calibration curves based on fluorescent intensities in the optimal cycle (maximum of the first derivative of the amplification curve) and corresponding curve parameters; A, B: Cycle 36 of the amplification of TFF1 , 1 ng of initial DNA; A: slope, R², y-axis intercept; B: whisker plots of Fisher scores; C, D: Cycle 35 of the amplification of S100A2, ng of initial DNA; C : slope, R², y-axis intercept; D: whisker plots of Fisher scores

Figure 4 shows the results of Example 1. Represented are the calibration curves based on Ct values and corresponding curve parameters, amplification of TFF1 on 1 ng of DNA; A: slope, R², y-axis intercept; B: whisker plots of Fisher scores Figure 5 shows the results of Example 1. Represented is a comparison of the curve parameters (slope, R², y-axis intercept, Fisher scores for differentiating adjacent methylation levels) of the calibration curves, which are obtained in different techniques for evaluation (based on fluorescent intensities in the optimal cycle or at the end point or based on Ct values) of amplification curves; A: Amplification of S100A2 on 10 ng of initial DNA; B: Amplification of TFF1 on 10 ng of initial DNA.The y-axis shows the values of the different quality parameters which are presented along the x-axis: a= linearity, b=slope, c=y-intercept, d=Fischer 0:5; e=Fischer5:10; f=Fischer 10:25; g=Fischer25:50; h=Fischer 50:75; l=Fischer75:100. The black columns represent the present invention calculating the methylation rate by the optimal amplification cycle. The white columns represent determination by end point analysis, and the grey columns represent the Ct-value analysis.

Figure 6 shows the results of Example 3. The y-axis shows the methylation rate in percent. Nine different samples were investigated, each with 50 ng (left), 10 ng (second from the left), 5ng (second from the right) and 1 ng (right) input of bisulfite DNA. In any case the standard deviation does not exceed 5%.

Figure 7 shows the result of Example 4. 12 different QM assays were conducted in five separate runs. The y-axis shows the methylation rate in percent. The different runs showed a low intra- and inter-plate variability.

Figure 8 shows the results of Example 4. 12 different QM assays were conducted in five separate runs. The y-axis shows the methylation rate in percent, the x axis the number of repetitions. The calculated confidence interval is around ± 5 percentage points of the mean of the methylation rate.

Figure 9 shows the results of example 5 (chip). The X axis shows the metastasis free survival times of the patients in years, and the Y axis shows the proportion of recurrence free survival patients in %. The lower curve shows the proportion of metastasis free patients in the population with above median methylation levels, and the upper curve shows the proportion of metastasis free patients in the population with below median methylation levels.

Figure 10 shows the results of example 5 (QM assay). The X axis shows the metastasis free survival times of the patients in years, and the Y axis shows the proportion of recurrence free survival patients in %. The lower curve shows the proportion of metastasis free patients in the population with above median methylation levels, and the upper curve shows the proportion of metastasis free patients in the population with below median methylation levels.

Figure 11 shows the correlation of measured methylation values using the chip platform (Y axis) and the assay of the present invention (Y-axis) of each patient. The correlation co-efficient is 0.87.

Claims

Patent Claims

1) Method for the quantification of methylated DNA, hereby characterized in that a) the DNA to be investigated is reacted in such a way that 5-methylcytosine remains unchanged, while unmethylated cytosine is converted into uracil or into another base, which is distinguished from cytosine in its base-pairing behavior, b) the converted DNA is amplified in the presence of two real-time probes, whereby one of the probes is specific for the methylated state, and the other probe is specific for the unmethylated state of the DNA; c) it is determined at different time points how far the amplification has proceeded by detecting the hybridization of the probes to the amplificates, d) the degree of methylation of the investigated DNA is determined.

2) The method according to claim 1, further characterized in that an exponential amplification method is used in step b).

3) The method according to claim 2, further characterized in that a PCR is used.

4) The method according to at least one of claims 1 to 3, further characterized in that the amplification is conducted by means of primers that are not methylation-specific.

5) The method according to at least one of claims 1 to 4, further characterized in that Lightcycler, Taqman, Sunrise, Molecular Beacon or Eclipse probes are used as real-time probes.

6) The method according to claim 4, further characterized in that Taqman probes are utilized in combination with minor groove binders as real-time probes.

7) The method according to at least one of claims 1 to 6, further characterized in that the amplification is conducted in the presence of both probes in the same vessel.

8) The method according to at least one of claims 1 to 7, further characterized in that the degree of methylation is calculated from the ratio of the signal intensities of the two probes at a specific time point.

9) The method according to claim 8, further characterized in that the degree of methylation is determined from the ratio of the signal intensities at a time point during the exponential amplification phase.

10) The method according to claim 9, further characterized in that the degree of methylation is determined from the ratio of the signal intensities at a time point which lies 5 cycles before or after the time point at which the amplification reaches its maximal slope (inflection point of the fluorescent intensity curves).

11) The method according to claim 10, further characterized in that the degree of methylation is determined from the ratio of the signal intensities at a time point which lies 2 cycles before or after the time point at which the amplification reaches its maximal slope.

12) The method according to claim 11 , further characterized in that the degree of methylation is determined from the ratio of the signal intensities at a time point which lies 1 cycle before or after the time point at which the amplification reaches its maximal slope.

13) The method according to claim 12, further characterized in that the degree of methylation is determined from the ratio of the signal intensities at a time point at which the amplification reaches its maximal slope (inflection point of the fluorescent intensity curves).

14) The method according to at least one of claims 1 to 7, further characterized in that the calculation of the degree of methylation is made by means of the ratio of threshold values, at which a certain signal intensity is exceeded.

15) The method according to claim 14, further characterized in that the calculation is made by means of the ratio of Ct values.

16) The method according to claim 15, further characterized in that the calculation is made by means of the following formula: degree of methylation = 100/(1 +2^ΔCt).

17) The method according to at least one of claims 1 to 7, further characterized in that the calculation of the degree of methylation is made by means of the ratio of the area under the fluorescent intensity curves or by means of the maximal slope of the curves.

18) The method according to at least one of claims 1 to 17, further characterized in that prior to step a), the assay conditions are optimized such that the fluorescent intensity curves provide a y-axis intercept that is as small as possible and a Fisher score that is as high as possible for the time point of the exponential amplification.

19) The method according to at least one of claims 1 to 18, further characterized in that prior to step a), the assay conditions are optimized such that the fluorescent intensity curves have a slope and a regression close to the value 1 for the time point of the exponential amplification.

20) A method for the absolute determination of the content of methylated DNA, hereby characterized in that a) a method is conducted according to one of claims 1-19, b) the result is compared with that of a standard curve.

21) The method according to at least one of claims 1 to 20, further characterized in that the quantification is carried out for the diagnosis of cancer diseases or other disorders associated with a change in the methylation status.

22) The method according to at least one of claims 1 to 20, further characterized in that the quantification is carried out for predicting undesired drug interactions and for the differentiation of cell types or tissues, or for the investigation of cell differentiation.

23) A kit which comprises two primers, a polymerase, as well as a probe specific for the methylated state and a probe specific for the unmethylated state, as well as, optionally, contains additional reagents necessary for a PCR and/or a bisulfite reagent.

24) A method for the quantification of two different variation of a DNA sequence characterized in that the following steps are conducted: a) the DNA is amplified in the presence of two real-time probes, whereby one of the probes is specific for one variation of the DNA sequence, and the other probe is specific for the other variation of the DNA sequence; b) it is determined at different time points how far the amplification has proceeded by detecting the hybridization of the probes to the amplificates, c) the proportions of the two sequence variations are determined.

25. A kit comprising two primers, a polymerase, a probe specific for one variation of the DNA sequence and a probe specific for the other variation of the DNA sequence to be investigated, and optionally additional reagents necessary for a PCR.

26) A method for the quantification of allele-specific gene expression, hereby characterized in that: a) the RNA to be investigated is reverse-transcribed, b) the cDNA is amplified in the presence of two real-time probes, whereby one of the probes is specific for one of the alleles and the other probe is specific for the other allele, c) it is determined at different time points how far the amplification has proceeded by detecting the hybridization of the probes to the amplificates,

d) the allele-specific gene expression is quantified. 27) A method for the investigation of SNPs from pooled samples, hereby characterized in that: a) the sample to be investigated is amplified in the presence of two real-time probes, whereby one of the probes is specific for the sequence of one SNP, and the other probe is specific for the sequence of the other SNP, b) it is determined at different time points how far the amplification has proceeded by detecting the hybridization of the probes to the amplificates, c) a conclusion as to the SNPs is made.

28) A method for determining the fraction of a wild-type microorganism and the fraction of a mutant strain in a sample, hereby characterized in that: a) the sample to be investigated is amplified in the presence of two real-time probes, whereby one of the probes is specific for the sequence of the wild type, and the other probe is specific for the sequence of the mutant strain, b) it is determined at different time points how far the amplification has proceeded by detecting the hybridization of the probes to the amplificates, c) it is concluded from this which strain at what fraction is represented in the sample.