RISK STRATIFICATION METHOD FOR A PATIENT HAVING A POLYMORPHISM

Abstract

A risk stratification method for a patient in a disease state and specifically patients presenting a tumor, includes determining if the patient is a homozygote or heterozygote and further determining the allelic expression for the patient, CC, T/C, or C/T. For patients having the cytosine methylated, they have a T/C allelic expression and patients without a methylated cytosine have a C/T allelic expression. A patient with a TT allelic expression is classified as a highest risk patient, a patient with a T/C allelic expression is classified as a second highest risk patient, a patient with a C/T allelic expression is classified as a third highest risk patients and a patient with a CC allelic expression is classified as a lowest risk patient. The risk stratification method may further include identification of an abnormal expression or mutation/function of a gene product produced by CTCF biding site 6.

Claims

1. A risk stratification method for a patient in a disease state comprising the steps of: a) determining whether said patient is a homozygote patient or a heterozygote patient with respect to CTCF biding site 6: b) determining whether the homozygote patient has a Thymine/Thymine, (TT) allelic expression or a Cytosine/Cytosine (CC) allelic expression for said CTCF biding site 6; c) determining for the heterozygote patient whether the cytosine is methylated; whereby when the cytosine is not methylated the patient has a Thymine/Cytosine, (TC) allelic expression for said CTCF biding and whereby when the cytosine is methylated the patient has a Cytosine/Thymine (CT) allelic expression for said CTCF biding site 6; wherein the TC allelic expression has a paternal cytosine and the CT allelic expression has a maternal cytosine; d) risk stratifying the patient according to their CTCF biding site 6 allelic expression wherein; i) the patient with a TT allelic expression is classified as highest risk, patient; ii) the patient with a TC allelic expression classified as a second highest risk patient: iii) the patient with; a CT allelic expression is classified as a third highest risk patient; iv) the patient with a CC allelic expression is classified as a lowest risk patient.

2. The risk stratification method of claim 1, wherein the step of determining for the heterozygote patient whether the cytosine is methylated comprises bisulfite conversion and quantitative methylation sensitive pyrosequencing.

3. The risk stratification method of claim 1, wherein the step of determining for the heterozygote patient whether the cytosine is methylated comprises directly sequencing parental DNA.

4. The risk stratification method of claim 1, wherein the method further comprises the step of determining the disease state of said patient by identification of a patient tumor.

5. The risk stratification method of claim 4, wherein the tumor is an infantile hemangioma.

6. The risk stratification method of claim 4, wherein the tumor is breast cancer.

7. The risk stratification method of claim 4, wherein the tumor is ovarian cancer.

8. The risk stratification method of claim 4, wherein the tumor is testicular cancer.

9. The risk stratification method of claim 4, wherein the tumor is liver cancer.

10. The risk stratification method of claim 4, wherein the tumor is lung cancer.

11. The risk stratification method of claim 4, wherein the tumor is brain cancer.

12. The risk stratification method of claim 1, wherein the method further comprises the step of determining the disease state of said patient by identification of an abnormal expression or mutation/function of a gene product produced by CTCF biding site 6.

13. The risk stratification method of claim 12, wherein the gene product is IGF2.

14. The risk stratification method of claim 12, wherein the gene product is H19.

15. The risk stratification method of claim 12, wherein the gene product is H19 antisense.

16. The risk stratification method of claim 12, herein the gene product is IGF2 antisense.

17. The risk stratification method of claim 12, wherein the gene product is a micro rna within the gene locus.

18. The risk stratification method of claim 12, wherein the gene product is an isoform.

19. The risk stratification method of claim 1, wherein the method further comprises the step of determining the disease state of said patient by identification of an abnormal expression or mutation/function of a binding gene product that binds to the CTCF biding site 6.

20. The risk stratification method of claim 19, wherein the binding gene product is CTCF.

21. The risk stratification method of claim 19, wherein the binding gene product is BORIS.

22. The risk stratification method of claim 19, wherein the binding gene product is a biding isoform.

23. The risk stratification method of claim 1, wherein the method further comprises the step of determining the disease state of said patient by identification of an abnormal expression or mutation/function of a binding partner of CTCF biding site 6.

24. The risk stratification method of claim 1, further comprising the steps of; a) determining an expression level of BORIS and an expression level of CTCF with respect to CTCF biding site 6; b) determining a % CTCF that is the expression level of CTCF as a function of the sum of the expression level of CTCF and the expression level of BORIS; c) risk stratifying the patient according to said % CTCF wherein; i) wherein when the % CTCF is less than 20%, the patient is classified as a highest % CTCF risk patient; ii) wherein when the % CTCF is less than 50% and greater than 20%, the patient is classified as a second highest % CTCF risk patient; iii) wherein when the % CTCF is less than 80% and greater than 50%, the patient is classified as a third highest % CTCF risk patient; iv) wherein when the % CTCF is 100% or less and greater than 8 the patient is classified as a lowest % CTCF risk patient;

25. The risk stratification method of claim 1, further comprising the steps of: a) determining an expression level of BORIS and an expression level of CTCF with respect to CTCF biding site 6; b) determining a CTCF-BORIS difference which is a difference in the expression level of CTCF to the expression level of BORIS; c) risk stratifying the patient according to said BORIS/CTCF difference wherein; i) when the CTCF-BORIS difference is less than zero, the patient is classified as a highest CTCF-BORIS risk patient; and ii) when the CTCF-BORIS difference is greater than zero the patient is classified as a lowest CTCF-BORIS risk patient.

26. The risk stratification method of claim 1, further comprising the steps of: a) determining an expression level of BORIS and an expression level of CTCF with respect to CTCF biding site 6; b) determining a % CTCF that is the expression level of CTCF as a function of the sum of the expression level of CTCF and the expression level of BORIS; c) determining a rate of % CTCF; d) risk stratifying the patient according to said rate of % CTCF wherein; i) a patient with a negative rate of % CTCF is a highest rate of % CTCF risk patient; ii) a patient with a neutral or positive rate of % CTCF is a lowest rate of % CTCF risk patient.

27. A risk stratification method for a patient in a disease state comprising the steps of: a) determining an expression level of BORIS and an expression level of CTCF with respect to CTCF biding site 6; b) determining a % CTCF that is the expression level of CTCF as a function of the sum of the expression level of CTCF and the expression level of BORlS; c) determining a rate of % CTCF; d) risk stratifying the patient according to said rate of CTCF wherein; i) a patient with a negative rate of % CTCF is a highest rate of % CTCF risk patient; ii) a patient with a neutral or positive rate of % CTCF is a lowest rate of % CTCF risk patient.

28. The risk stratification method of claim 27, wherein the method further comprises the step of determining the disease state of said patient by identification of a patient tumor.

29. The risk stratification method of claim 27, wherein the tumor is an infantile hemangioma:

30. A medication response prediction method for a patient in a disease state comprising the steps of: a) determining whether said patient is a homozygote patient or a heterozygote patient with respect to CTCF biding site 6; b) determining whether the homozygote patient has a Thymine/Thymine, (TT) allelic expression or a Cytosine/Cytosine (CC) allelic expression for said CTCF biding site 6; c) determining for the heterozygote patient whether the cytosine is methylated; whereby when the cytosine is not methylated the patient has a Thymine/Cytosine, (TC) allelic expression for said CTCF biding site 6; and whereby when the cytosine is methylated the patient has a Cytosine/Thymine (CT) allelic expression for said CTCF biding site 6; wherein the TC allelic expression has a paternal cytosine and the CT allelic expression has a maternal cytosine; wherein a non-TT allelic expression includes a TC allelic expression, a CT allelic expression and a CC allelic expression; d) predicting a medication response for said patient according to their CTCF biding site 6 allelic expression wherein; i) for the patient with a non-TT allelic expression, the medication response is predicted to be most effective; ii) for the patient with a TT allelic expression, the medication response is predicted to be least effective.

31. The medication response prediction method of claim 30 wherein the method further comprises the step of determining the disease state of said patient by identification of a patient tumor.

32. The medication response prediction method of claim 31 wherein the medication is a cortical steroid.

33. The medication response prediction method of claim 31, wherein the tumor is an infantile hemangioma.

34. The medication response prediction method of claim 31, wherein the tumor is a cancerous tumor.

35. The medication response prediction method of claim 30, wherein the method further comprises the step of determining the disease state of said patient by identification of an abnormal expression or mutation/function of a gene product produced by CTCF biding site 6.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0022] The accompanying drawings are included to provide a further understanding of the invention and are incorporated in and constitute a part of this specification, illustrate embodiments of the invention, and together with the description serve to explain the principles of the invention.

[0023] FIG. 1 shows a drawing of an infantile hemangioma.

[0024] FIG. 2A shows a diagram representing the genotyping via direct sequencing of blood samples.

[0025] FIG. 2B shows a chart of bisulfite conversion and quantitative methylation sensitive pyrosequencing.

[0026] FIGS. 3A and 3B show a chart of IGF2 transcription by clinical stage wherein FIG. 3A shows IGF2 Transcript and FIG. 3B shows % CTCF.

[0027] FIG. 4A shows a graph of the percent CTCF with respect to CTCF and BORIS versus the age of the lesion.

[0028] FIG. 4B shows the CSUM of % CTCF versus the age of the lesion.

[0029] FIG. 5 shows a graph of IGF2 Transcript versus the % CTCF.

[0030] FIG. 6 shows Western Analysis of 24 IH Samples Via 5 Independent Western Blots.

[0031] FIG. 7 shows a schematic integrating CTCF and BORIS expression via both transcript and protein with clinical stage.

[0032] FIGS. 8A and 8B show a graph of the parent of origin, paternal and maternal respectively, specific effects of CTCF BS6 on IGF2 and H19 transcription.

[0033] FIG. 9A shows a graphs of percent methylation versus CpG number.

[0034] FIG. 9B shows a graph of percent methylation of the H19 promoter by pyrosequencing versus CTCF transcript.

[0035] FIG. 9C shows methylation of the H19 promoter.

[0036] FIG. 9D shows a sequence showing the H19 promoter, CpG#4 of the bisulfite sequencing test corresponds to the CCCGGG Pst1 digestion, site of the Southern analysis.

[0037] FIG. 10A shows a graph of the diameter of an infantile hemangioma versus the age of the lesion in days for TT and non-TT patients.

[0038] FIG. 10B shows a graph of the diameter of an infantile hemangioma versus the age of the lesion in days for TC, C/T and C/C patients.

DETAILED DESCRIPTION OF THE ILLUSTRATED EMBODIMENTS

[0039] Corresponding reference characters indicate corresponding parts throughout the several views of the figures. The figures represent an illustration of some of the embodiments of the present invention and are not to be construed as limiting the scope of the invention in any manner. Further, the figures are not necessarily to scale, some features may be exaggerated to show details of particular components. Therefore, specific structural and functional details disclosed herein are not to be interpreted as limiting, but merely as a representative basis for teaching one skilled in the art to variously employ the present invention.

[0040] As used herein, the terms comprises, comprising includes, including, has, having or any other variation thereof, are intended to cover a non-exclusive inclusion. For example, a process, method, article, or apparatus that comprises a list of elements is not necessarily limited to only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Also, use of a or an are employed to describe elements and components described herein. This is done merely for convenience and to give a general sense of the scope of the invention. This description should be read to include one or at least one and the singular also includes the plural unless it is obvious that it is meant otherwise.

[0041] Certain exemplary embodiments of the present invention are described herein and are illustrated in the accompanying figures. The embodiments described are only for purposes of illustrating the present invention and should not be interpreted as limiting the scope of the invention. Other embodiments of the invention, and certain modifications, combinations and improvements of the described embodiments, will occur to those skilled in the art and all such alternate embodiments, combinations, modifications and improvements are within the scope of the present invention.

[0042] This application incorporates by reference in its entirety, the following publication; [0043] Brent Schultz, Xiaopan Yao, Yanhong Deng, Milton Waner, Christopher Spock, Laura Tom, John Persing, Deepak Narayan, (2015) A Common Polymorphism Within The IGF2 Imprinting Control Region Is Associated with Parent of Origin Specific Effects in Infantile Hemangiomas PLOS ONE|DOI:10.1371/journal.pone.0113168

Materials and Methods

[0044] This, article describes a study that involved analyzing post excisional tissue from surgical candidates. The decision to operate was in no way influenced by the study. Clinical data was gathered retrospectively from this same group. All surgical candidates had clinical measurements available for analysis. All samples and clinical data were collected in accordance with the approved HIC protocol (#0507000430) as reviewed by the Yale University Medical School IRB. This protocol was approved specifically for this study. Written consent was obtained from each patient's legal guardian prior to surgery. All data obtained including clinical measurements were stored in a de-identified format.

Specimen Collection:

[0045] Please refer to the Master Data Table for details. Those specimens later confirmed to be hemangioma tissue, via Glut-1 positive histology, were considered for this project. Further, only discrete solitary lesions that were not found in the setting of a syndrome were considered. Those patients where prior surgical resections of the lesion were performed were also deemed ineligible. Of note, those lesions previously treated with laser were not excluded, as the effects of laser treatment are relatively superficial. However, during specimen collection, all areas that appeared grossly to be affected by laser treatment were excised before further processing. Briefly, forty-two samples were collected (See Table 1 for details.) Of these, two samples (#41 and #42), were excluded from all analyses; sample #41 was Glut-1 negative on histology and #42 was subject to prior resections. Nineteen IH samples were selected at random for methylation analysis of the H19 promoter by southern. These 19 samples were also analyzed for methylation specific pyrosequencing of the same region. An additional two samples (numbers 1 and 14) were also subjected to H19 methylation specific pyrosequencing to bolster the number of samples with both H19 methylation and transcriptional data. Regarding transcriptional analysis, nineteen samples were found to have suitably intact RNA for quantitative RT PCR. Specimens for transcriptional analysis were separated into three categories: 1) Proliferative, 2) Quiescent, and 3) Involuting phases. As a lesion's stage is, by definition, clinical, an experienced physician staged the IH at the time of surgery. Data regarding clinical stage was gathered prospectively. The determination of the clinical stage was made by one of three highly experienced surgeons regarding vascular anomalies, using interval growth, patient age and the color/turgor of the lesion at the time of resection as criteria. General characteristics of these categories are as follows: 1) proliferative hemangiomas were generally less than 1.5 years of age with interval growth between the last two clinical visits preceding surgery, no lightening of lesion color was noted. 2) Quiescent, hemangiomas: no interval growth between the last two clinic visits preceding surgery, lightening of color also played a factor in these determinations. Involuting hemangiomas: interval regression by measurement between the last two clinic visits preceding surgery, further color changes were often but not always noted. Table 1:

[0046] In total 34 samples were genotyped for a polymorphism within CTCF BS 6 and parental contributions were determined for heterozygotes (see FIG. 2 and bisulfite sequencing methods for details.) Only lesions of the head and neck were included. Of these samples, 3 were excluded because they were not on the head or neck. Two other samples were excluded because one was not Glut-1 positive on histology and one had a previous resection of the same lesion prior to evaluation, (See Table 1 for details.) Thus 29 individuals were included in this analysis. Charts were retrospectively reviewed from patients treated either at the Yale University Plastic Surgery Center (New Haven Conn.), or the Vascular Birthmarks Institute (New York, N.Y.). The age of the lesion was then compared to the size of the lesion as determined below. These data were plotted and separated by CTCF BS6 genotype and parental contribution in heterozygotes. ANCOVA analyses were then performed on putative growth curves. The age at the time of resection, with corresponding size, was used only for proliferating lesions. For those lesions that were resected at the time of involution, or quiescence, the size of the lesion at the clinical visit where quiescence was first noted was used. In the case of medical interventions, the size of the lesion before a response was noted was used. This information was used to create a clinical table of results were factors such as ulceration, steroid/chemotherapeutic, and laser treatment were also noted (note that no beta blockers were used in the sample population.) Thus, different ages are associated with most individuals when comparing the clinical data table and the master data table. For instance: Sample #7 has an age of 2304 days assigned in the aster Data Table 1. However, in the Clinical Table the age assigned to sample #7 is 1050 days. The difference in age assignments represent the age of the patient when the lesion was excised (this age was used for the molecular analyses) versus the age of the patient either before the first response to medical intervention was noted or when the lesion first entered the quiescent phase. Thus, for sample number 7: The lesion entered the quiescent phase (as determined retrospectively) at age 1050 days but was then excised at age 2304 days. Regarding the assessment of lesion size, if multiple dimensions were given, the largest was used. In some cases, only one dimension was given so volumetric estimations could not be calculated for every patient. Thus, patients' lesions were standardized to a greatest diameter equivalent. This measurement was correlated with clinical photographs when available. All data utilized varied by less than 10% between stated measurement and photographic estimation when available. Lesions were classified into one of three growth phases: proliferative, quiescent and involuting. Only sizes of lesions that were in the proliferative or quiescent phase were used in this study. All data was stored in a de-identified format with a unique accession number for each patient.

DDNA Preservation and Extraction

[0047] Immediately following tissue resection, DNA was isolated using the Qiagen DNeasy Tissue Mini Kit according to the manufacturer's protocol. Only samples with an A260/A280 measurement of 1.8 or above that ran as a single band on the gel were further analyzed,

[0048] FIG. 2. Deducing Parental Contributions From Direct Sequencing and Bisulfite Pyrosequencing. FIG. 2A: 29 patients were genotyped via direct sequencing of blood samples for a known polymorphism within the core CTCF BS 6 sequence (rs10732516.) All homozygous genotypes could be deduced from this information alone. FIG. 2B: All samples (heterozygotes and homozygotes) were subjected to bisulfite conversion and quantitative methylation sensitive pyrosequencing Methylation occurs only on the paternal chromosome for CTCF BS6. In normal tissue, such as patient matched control blood, this assay is capable of isolating the genotype of the paternal chromosome. As thymidine cannot be methylated, those individuals with a paternal T at rs10732516 were not methylated at CpG#5, Paternal C carrying individuals were methylated at CpG#5. Thus, the maternal and paternal contribution to CTCFBS6 can be deduced. This assay sidesteps the need for directly sequencing parents' DNA and eliminates the potential ambiguity ensuing from heterozygous parents. Note: The methylation values of this assay are subject to primer bias, Tost et al (25.) This is evident by the 3 distinct groupings of methylation levels, which are artifactual.

RNA Preservation and Extraction

[0049] Immediately following tissue resection, 100-500 mg of tissue was stored in Quiagen RNA Later solution according to the manufacturer's protocol. RNA was extracted via liquid nitrogen powder homogenization using Invitrogen Trizol reagent according to the standard protocol, 10 g of total RNA from each sample was then treated with DNase Qiagen mini-elute columns according to manufacturer's specifications. RNA integrity was then assessed using Agilent bioanalyzer 2100 (provided as a service of the Keck Center at Yale University.) Those samples with 18s/28s ratio of 1.8 or greater were converted into cDNA using the ABI 4368813 cDNA: archive kit. All samples were then stored at 80 degrees C.

Quantitative rtPCR for CTCF, BORIS, H19 and IGF2

[0050] Nineteen IH samples with suitable RNA, as, previously specified, were subjected to fluorescent quantitative RT-PCR using ABI Taqman primers that were previously validated by the manufacturer and spanned intron exon boundaries. For reasons of sample scarcity, not all samples were subjected to each assay (See FIG. 1 for details.) The assays were: IGF2assay number Hs00171254_m1, H19assay number Hs00399294_g1, CTCFassay number Hs00198081_m1, and BORISassay number Hs00540744_m1. Gene quantification was performed using the standard curve method via pooled sample cDNA (equal contributions from each sample) and successive two fold dilutions, beginning from 50 ng and ending with 0.39 ng. All reactions were performed on the ABI 79005 thermocycler using default cycling conditions previously optimized for these assays. Reactions were performed in duplicate and average CT values, if they agreed within 0.4 cycles, were used to calculate absolute quantity. Three runs of RT PCR were performed with overlapping samples in each run to allow normalization of the data. Not all samples were subjected to every assay depending upon sample quantity. Of Note: Sample #4 does not have an H19 transcription value, as on duplicate plating for rtPCR, the CT values did not agree within 0.4 cycles. Furthermore, samples 15-19 were the final rtPCR of the three runs performed and due to sample scarcity and the need to construct standard curves from pooled samples, only CTCF and BORIS rtPCR's were performed.

Western Analysis

[0051] 24 samples were subjected to Western analysis. As this process is tissue intensive, younger samples such as #21 and #23 could only be used for this analysis as insufficient tissue was left for further processing. Other samples were selected biased toward analyzing those samples with transcriptional data in order to compare transcriptional phenomena to translational. However, as the analysis proceeded, presentation gels were constructed to demonstrate key transition points in CTCF and BORIS translation in samples that had not been treated with steroids. Briefly, 50 mg of each sample were processed with a rotary homogenizer in 200 ml of RIPA lysis buffer. After centrifugation lysates were created, using a standard beta-mercapto-ethanol with SDS. PAGE was performed with 36 g of protein per well in NuPage 10% Bis-Tris precast gels inMOPS buffer at 100 volts. PAGE separated proteins were then transferred for two hours to a PVDF membrane (Bio-Rad) in a standard transfer buffer at 100 mAmps, Anti BORIS antibody (Abcam 18337) was used at 1/5000 dilution in TBST with 5% cows milk overnight. Two concentrations of anti-CTCF were used1:10,000 and 1:5,000to better visualize late and early rises in CTCF protein (see FIG. 3 legend for details.) As anti-CTCF and anti-BORIS were both rabbit polyclonal antibodies they could be visualized simultaneously on the same film following incubation with the anti-rabbit secondary conjugated to horseradish peroxidase and ECL treatment. Images were then scanned and adjusted for brightness and contrast in Adobe Photoshop.

Bisulfite Methylation Analysis Using Quantitative Pyrosequencing

[0052] This method was first described by Grunau et al and Dupont et al. Protocols specific for each assay in this study are available in S1 Supplementary Information or upon request. Incomplete bisulfite conversion was detected by designing amplicons that contained at least 1 unmethylated cytosine. Primer bias was controlled for by establishing methylation curves of 100% methylated DNA titrated against known amounts of whole genome amplified PCR products that, by definition, are unmethylated. These methylation curves allow experimental samples to be calibrated against known standards. Primers for Exon 9, the H19 promoter and CTCFBS6 as well as the bisulfite-converted sequences they amplify are available in S1 Supplementary Information. The presence of an A/G polymorphism, approximately 130 base pairs downstream of CTCF BS6 leads to primer bias and distorts the absolute methylation values of CTCF BS6, Tost et al.

Deducing Parental Contributions of Alleles at CTCF BS6

[0053] All samples were subjected to direct sequencing of CTCF BS6 containing the polymorphism rs10732516. Primer design and reaction conditions are available in S1 Supplementary Information. The DNA samples were subjected in parallel to methylation sensitive pyrosequencing of the same polymorphism; please see the section titled Specimen Collection for further details. Comparing these results allows each parental contribution to be deduced, see FIG. 1 for full details.

[0054] As shown in FIG. 3, IGF2 transcription by clinical stage shows an inverse relationship to the % CTCF of identical stage. FIG. 3A: IGF 2 mRNA levels were approximately 6 lower in the involuting samples compared to their proliferating counterparts. Proliferating vs. involuted p=0.02, Proliferating vs. quiescent p=0.0081, Quiescent vs involuting p=0.01. Wilcoxian Rank Sum Test. Error bars represent standard deviation. FIG. 3B: % CTCF changes significantly according to clinical stage. Prolif vs. quiescent p=0.01, quiescent vs. invol p=0.006, prolif vs. invol p=0.63, Wilcoxian two sample test. Error bars represent standard deviation. Note: All samples in the IGF2 analysis were included in the % CTCF analysis, with additional samples.

Genomic Southern Analysis for the H19Promoter

[0055] Eighteen were analyzed at a CLIA certified molecular diagnostic laboratory where this assay is performed as a clinical test for Beckwith-Wiedemann Syndrome. The assay is originally described by Debaun et at Norms for this test were previously established with 30 normal controls at 55% methylation with a standard deviation of 5%. All samples were run with a normal and Beckwith-Wiedemann control. The assay exploits a CCCGGG site in the H19 promoter that is cut by the methylation sensitive restriction enzyme Pst1.

Statistical Analysis

[0056] Descriptive statistics were used to present patient characteristics. The difference in expression of IGF2 transcript across the three developmental stages of IH was evaluated using a Kruskal-Wallis test. To evaluate if the relative amount of CTCF compared to BORIS transcript changes predictably over time, change point analysis was performed. Change point analysis indicates the likelihood that a change in transcript expression occurred in the sample population by confidence level and a confidence interval regarding when those changes occur. The % CTCF [CTCF/(CTCF+BORIS)100] was used to develop a change point model that was then compared against clinical staging and IGF2 expression in the sample population. A full explanation of the methods used, as well as a shareware change-point analyzer is presented as an online resource: Taylor, Wayne A. (2000), Change-Point Analysis: A Powerful New Tool For Detecting Changes, (http://www.variation.com/cpa/tech/changepoint.html.) To evaluate the association of IGF2 transcript and the relative amounts of CTCF, a linear regression model was fitted, with the % CTCF and age as covariates. The correlation and partial correlation were, also calculated. Partial correlations indicate what percentage of variance in IGF2 can be explained by CTCF % alone. Analysis of the covariance model (ANCOVA) was fitted to examine if the correlation between the IGF2 transcript and the difference between CTCF and BORIS varied by the paternal genotype at CTCF CS6, once adjusted by age.

Results Master Data Table

[0057] In total, 40 samples were analyzed on a molecular basis. A description of basic demographics, with genotypes at CTCF BS 6 and transcript expression values for IGF2, H19, CTCF and BORIS with correlative methylation data was compiled. Please see FIG. 1: Master Data Table for details. This table can be utilized to confirm any statistical analysis presented in this study.

Expression of IGF2, CTCF and BORIS

[0058] FIGS. 4A and 4B Analyzing the Percentage of CTCF Compared to Total CTCF and BORIS Transcript.

[0059] FIG. 4A Using the two samples with the lowest % CTCF (420 and 418 days as common points (Purple)) two curves with high correlation to age can be appreciated. % CTCF steadily decreases as lesions age until approximately 400 days (red and purple points), then CTCF once again increases compared to BORIS (purple and blue points). This roughly correlates with the transition from proliferative to quiescent lesions. FIG. 4B CSUM of % CTCF demonstrates statistically significant variation about the mean of % CTCF according to age. For a full explanation of the CSUM data and commensurate change point analysis.

[0060] IGF2 transcription differed significantly by clinical stages (p<0.0001, Kruskal-Wallis test). Plateau stage lesions expressed significantly higher levels of IGF2 than proliferating (p=0.0081, Wilcoxon rank sum test) and involuted samples (p=0.02). Involuted hemangiomas expressed the lowest levels of IGF2, approximately 6 lower than their proliferating counterparts (p=0.01).

[0061] To potentially explain the changes in IGF2 transcription, quantitative RT-PCR was performed for CTCF and BORIS. CTCF and BORIS are co-expressed in all samples. However, the percentage of CTCF transcript compared to total CTCF and BORIS in a given sample [CTCF/(CTCF+BORIS)100] varied significantly over developmental time (FIGS. 3B and 4A). This was confirmed by a change point model (FIG. 4B): the Y axis is the cumulative sum (CUSUM) of the differences between % CTCF and the average value of % CTCF. A segment of the CUSUM chart with an upward slope indicates a period where the values tend to be above the overall average. Likewise, a segment with a downward slope indicates a period of time where the values tend to be below the overall average. Based on this analysis, two change points, one estimated at 418 days and the other at 1277 days, were detected. Prior to approximately 418 days (90% Cl: 368-547 days), the value of % CTCF tends to maintain a nigher level with an average value in this stage equal to 59%. In the second stage (418-1277 days), the level of % CTCF is low with an average 28%. After the second change point at approximately 1277days (90% Cl: 760-1500 days), % CTCF has recovered to a high level again with the average 66%. See S2 Supplementary Information for a bar graph analysis of the change point model of % CTCF expression. This result is highly similar to the results obtained by separating samples according to clinical stage (S3 Supplementary Information) % CTCF transcription in IH varies according to clinical stage (FIG. 3B) Furthermore, the % CTCF varied inversely with IGF2 transcription (compare FIG. 3A and 3B) These two graphs clearly demonstrate that a higher % CTCF corresponds with lower levels of IGF2 expression. Moreover, a strong positive correlation was detected between BORIS and IGF2 transcription (p=0.0028). Though CTCF alone does not significantly correlate with IGF2 transcript levels, taking both CTCF and BORIS into account using % CTCF is the strongest predictor of IGF2 mRNA expression p=0.0004 (FIG. 5).

[0062] As described herein, a reduction, especially a rapid reduction, in the % CTCF, as expressed in FIG. 4A, 4B, and 5, may be an indication of an aggressive tumor and may be used, to risk stratify the patient. A % CTCF may be determined as a proportion of the sum of the expression level of CTCF and expression level of BORIS. It has been found that a reduction in the % CTCF overtime, or rate of % CTCF, in indicative of a more aggressive tumor, thereby putting the patient in an higher risk classification. For example, as shown in FIG. 4A, if the line fit of % CTCF versus age in days has a negative slope, then a patient may be classified in as a highest risk patient and if the slope is positive, then the patient may be classified as a lowest risk patient.

[0063] As shown in FIG. 5, IGF2 transcript levels correlate inversely with the percentage of CTCF compared to a total of CTCF+BORIS. This data represents the first demonstration of the potentially antagonistic effects of CTCF and BORIS on a target gene's transcription through a continuous curve. P=0.0004 ANCOVA Model. Age effect was not significant in the model p=241. The % CTCF is a stronger statistical predictor of IGF2 expression than BORIS alone 0.0004 vs. 0.0028 respectively. (N=15)

[0064] FIG. 6, shows Western Analysis of 24 IH samples via 5 independent Western Blots. Twenty samples with 4 blots depicted, demonstrates 4 stages of CTCF and BORIS expression. 6-1: A low concentration of anti-CTCF (1:10,000) demonstrates the complete spectrum of CTCF expression with increases early (210 to 418) and late (760 to 987) in protein expression. (Note, a testes negative control was included as well as a venous malformation denoted as VM.) 6-2 through 6-4 were probed with 1:5000 concentration of anti-CTCF that more clearly demonstrates the early rise in CTCF that occurs after 367 days. 6-2 suggests an early increase in BORIS with precipitous downregulation after 244 days. 6-3 and 6-4 expand this critical age range demonstrating a period from 244 days to 367 where BORIS is downregulated but CTCF is not yet upregulated. Note, samples marked with an asterisk were also subjected to CTCF and BORIS rtPCR.

[0065] Western analysis of CTCF and BORIS confirms and expands upon the transcript data (FIG. 6) As expected in proliferating lesions, BORIS transcript and protein levels steadily rise in early stage samples (FIG. 4A transcript data, FIG. 6-2 Western Analysis.) Furthermore, during the transition from quiescent to involuting samples, CTCF mRNA and protein increase compared to BORIS (FIGS. 3B, 4A and 4B transcript, and FIG. 6-1 protein.) Thus, the western and transcriptional data globally confirm one another at the endpoints of IH development. However, the protein data suggests a third change in CTCF and BORIS levels that the transcript change point analysis does not. This third proteomic change appears to take place at the late proliferating to early quiescent phase. It coincides with the so-called late proliferative stage in IH that is suggested by clinicians but not universally accepted. These data provide the first molecular support for what was previously a clinical category: the late proliferative stage of IH growth. The relative expression of CTCF and BORIS via both transcript and protein levels, is predictive of clinical stage and IGF2 expression. See FIG. 7 for details.

[0066] FIG. 7. shows a schematic integrating CTCF and BORIS expression via both transcript and protein with clinical stage. The Western analysis suggests 4 stages of CTCF and BORIS (see FIG. 3, panels 1 to 4) each stage leading to higher levels of CTCF expression relative to BORIS. These interval changes in protein expression closely correlate with clinical stage. Furthermore, bars above the stages represent the 95% confidence intervals of the two change points identified by quantitative rtPCR. These data show remarkable agreement reinforcing the idea that relative CTCF and BORIS expression levels closely mirror the clinical stage of the lesions tested. Of note, the CTCF and BORIS protein data also suggest a molecular basis for a late proliferative stage.

[0067] The CTCF to BORIS transcript difference (C-B) predicts IGF2 transcription according to the paternal allele at CTCF BS6. This study utilizes existing technologies: Direct sequencing of the known polymorphism of CTCF BS6 (rs10732516) with a previously described methylation assay for CTCF BS6. Applying these two assays in a novel manner (see FIG. 2) allowed us to deduce both maternal and paternal contributions to CTCF BS6which will be referred to as the maternal and paternal contribution (FIG. 8.) This paternal contribution likely has significant effects on IGF2 production as it relates to CTCF and BORIS.

[0068] IGF2 mRNA is demonstrated to be inversely related to CTCF and positively correlated to BORIS transcripts when plotted against % CTCF (FIG. 5).

[0069] By using the difference between CTCF and BORIS (C-B) rather than the % CTCF, this relationship can be differentiated by the paternally contributed allele at CTCF BS6 (FIG. 8A). The paternal allele of a common C/T polymorphism within CTCF BS6 (rs10732516) corresponds with two strikingly different CTCF-BORIS vs. IGF2 curves. After age adjustment, the effect of CTCF-BORIS on IGF2 transcription was significantly different between patients bearing different paternal alleles (p=0.05, ANCOVA model) and there is a strong correlation between IGF2 expression and CTCF-BORIS (p=0.0007). The samples bearing a paternal C allele, appear to demonstrate a six fold steeper slope of IGF2 mRNA relative to the CTCF-BORIS difference, compared to their paternal T bearing counterparts. This allele was identified in both tissue and patient matched control blood. Of note, heterozygote analysis revealed no clear relationship regarding these factors according to the maternal allele (p=0.95, ANCOVA model). It remains a possibility that the paternal allele effect may be steroid treatment driven as more samples with the paternal T allele were treated with steroids than the paternal C allele. This potential bias was investigated with an odds ratio calculation sorting steroid treatment according to paternal genotype. The odds ratio suggested that paternal T samples were more likely to be, treated with steroids but this result did not reach statistical significance. As the allele specific analysis was done on only proliferative samples the odds ratio calculation was performed twice, once including involuted samples and once to their exclusion. However, it is acknowledged that the odds ratio suggested a potential steroid treatment bias in paternal T samples, which may have become significant n a larger patient cohort. Lastly, IGF2 expression in IH was mono-allelic in all 5 informative heterozygotes tested for a known IGF2 polymorphism in exon 9. IGF2 imprinting status appears to be maintained despite BORIS expression.

[0070] H19 transcript levels correlate positively with CTCF mRNA according to the maternally contributed allele at CTCFBS6. After age adjustment, CTCF transcript levels alone correlated positively with H19 transcription but only when separated by maternal genotype ((p=0.0150, FIG. 8B) Moreover, this positive correlation is significantly different among patients with different maternal alleles (p=0.0129, ANCOVA). The correlation between CTCF and H19 transcription is stronger in patients bearing a maternal C allele compared to their maternal T counterparts. There were no identifiable relationships between H19 transcription and either the paternal genotype at CTCF BS6 or BORIS mRNA, p=0.8 ANCOVA (S4A Supplementary Information.) There also appeared to be no relationship between H19 expression and clinical stage of the, hemangioma (S4C Supplementary Information).

[0071] It remains a possibility that the maternal allele effect may be steroid treatment driven as more samples with the maternal T allele were treated with steroids than the maternal C allele. This potential bias was investigated with an odds ratio calculation sorting steroid treatment according to maternal genotype. The odds ratio suggested that maternal T samples were more likely to be treated with steroids but this result did not reach statistical significance (See S4F Supplementary Information.) However, it is acknowledged that this odds ratio may have become statistically in a larger sample size; the effect was not great enough to significantly bias the sample size.

[0072] CTCF transcript levels alone correlate with demethylation of the H19 Promoter. All IH samples tested demonstrated significant hypomethylation at the H19 promoter compared to matched patient blood controls (FIG. 9A). This was demonstrated with bisulfite specific pyrosequencing (FIG. 9B) and confirmed with methylation sensitive enzyme digest and southern analysis (FIG. 9C). As pyrosequencing is quantitative, these data were correlated with CTCF and BORIS expression in matched samples. Decreased promoter methylation correlated with higher levels of CTCF (p=0.015; simple regression. FIG. 9B). Yet, this finding may be subject to confounders as it is not statistically significant after age adjustment (p=0.17). As the H19 promoter is paternally methylated and maternally demethylated, a methylation level below 50% would entail demethylation of the paternal allele in at least a subset of cells within a given sample. However, no bi-allelic expression of H19 could be detected in 5 heterozygote samples (S5 Supplementary Information.) Furthermore, H19 promoter demethylation did not strongly correlate with H19 expression (S6 Supplementary Information.)

[0073] FIG. 9. shows CTCF expression and H19 promoter methylation. FIG. 9A: Increased CTCF transcript level correlates with demethylation of the H19 Promoter. Those samples with the highest CTCF expression were the least methylated ranging from 34% to 14%. However, demethylation of the H19 promoter did not correlate strictly with H19 transcript expression (S6 Supplementary Information). FIG. 9B and 9CThe H19 promoter (see FIG. 9D) is hypomethylated, demonstrated by bisultife converted pyrosequencing (9B) and methylation sensitive restriction digest with southern hybridization (9C.) 25 IH samples, and 13 matched blood controls were subjected to bisulfite converted pyrosequencing. 13 IH samples and 13 matched blood controls were subjected to southern analysis with methylation sensitive Pst1 digestion. Two representative gels show, 5 IH samples, a Beckwith-Weidman positive control and a 50% methylated normal control. FIG. 9D: sequence showing the H19 promoterCpG#4 of the bisulfite sequencing test corresponds to the CCCGGG Pst1 digestion site of the Southern analysis. Other CpG's tested are in bold. This CpG is in close proximity to the transcription start site of H19 (blue arrow) and an overlapping putative CTCF binding site identified by positional weight matrix analysis.

[0074] Multiple imprinted sites within the IGF2/H19 locus are abnormally methylated in IH compared to matched control blood.

[0075] It remains a formal possibility that the normalization of CTCF to BORIS ratios, as well as decreased IGF2 transcription, in involuting and involuted samples is not due to an intracellular phenomenon but rather the incremental replacement of abnormal IH tissue (vascular stroma) with normal tissue (fat.) Thus, the results presented are the product of tissue heterogeneity. It is acknowledged that IH lesions transform from a vascular tumor into a fibrofatty residuum; therefore, the transitional phases are by definition composed of heterogeneous cell populations. However; no evidence was observed that the fibrofatty residuum of an involuted IH represents normal tissue. To the contrary, many of the methylation abnormalities discovered by this study are either stable or progressive from early to late clinical stages. For instance, the H19 promoter is significantly demethylated in all IH samples (see FIG. 9A-9C). However, the demethylation is progressive with age (S7 Supplementary Information.) Furthermore, it was also found focal demethylation at Exon 9 and hypermethylation at DMR0, deviating from the expected 50% for these known imprinted sites. Both findings remained consistent in all IH clinical types and were age independent (S8 and S9 Supplementary Information for details.) If IH tissues were being replaced by normal fat it would be expect the methylation abnormalities demonstrated in this work to normalize, not remain constant or even progress with age. Given this argument, the simplest explanation for the methylation data is that IH tissue begins as epigenetically abnormal vascular stroma and transforms into epigenetically abnormal fat or is at least replaced by the like.

CTCF BS6 Genotypes Correlate With Clinical Outcomes

[0076] Mechanisms aside, parent of origin specific effects are demonstrated at the molecular level regarding expression patterns of both IGF2 and H19. However, the question remains whether these molecular phenotypes may translate into clinically significant growth patterns. Table 2 is a complete table of all, patients participating in this retrospective clinical study. For details of subject inclusion please see methods section. Size of the lesion as well as the date of examination was included with relevant clinical information such as sex, medical treatments utilized and presence of ulceration during clinical course. Each patient was sorted according to CTCF Binding Site Six Genotype and paternal contribution for heterozygotes. By plotting the size of IH lesions against the CTCF BS6 genotypes, four distinct growth curves emerge (FIG. 10A and 10B).

[0077] As shown in Table 2, many of the patients were treated with a cortical steroid-methyl pregnisone, either by injection or systemic oral intake. Table 3 shows that the patients with a non-TT allelic expression had a much better response to this medication, wherein only 25% had a failure of response to the medication. In contrast, the patients with aa TT allelic expression had the least effective response to the Cortical steroid treatment.

TABLE-US-00001 TABLE 3 Allelic Expression Failure of Steroids Non-Failure of Steroids TT 7 2 Non-TT 5 15

[0078] The results of this analysis shows that the type of allelic expression can be used to predict the effectiveness of Cortical steroid treatment. The odds ratio from this stud is 10.5 that the TT group will fail steroids, with a sensitivity of 58.33% (95% Cl 27-84%) and a specificity of 88.33% (95% Cl 63-98%), The failure of steroids is defined in this analysis as the patient requiring surgical treatment despite the treatment with the medication, cortical steroids.

[0079] FIG. 10. Clinical Correlation of Hemangioma Growth Rates with Parental Contributions to CTCF BS 6. FIG. 10A: This retrospective analysis of 29 samples, 9 TT, 20 non TT, demonstrates significantly distinct growth curves over a large age range. The ANCOVA model has identified age as a predictor of size p=0.0007. The association between tumor size and age is significantly different among the genotypes of TT, C/T, T/C and CC p<0.0001. Of Note the paternal contribution is presented first and the maternal is second. The interaction terms of parentally specific genotypes allowed us to test if the slopes of the curves between tumor size and age are different among the genotypes. This analysis indicted that an increase in 1 day of age is associated with 0.016 cm of growth in the TT group. This is significantly higher than the non TT group p=0.0019. FIG. 10B; Growth analysis focusing on the non TT group. Each non TT growth curve varied independently and significantly from the TT samples (CC vs. TT: P<0.0001, CT vs. TT: P<0.0008, TC vs. TT: P=0.0025). Furthermore, these data suggest parent of origin specific effects as those samples with identical genotypes but opposite parental contributions displayed statistically significant differences in growth curves. The paternal T/maternal C genotype grew at approximately twice the rate as their paternal C/maternal T carrying counterparts (p=0.05). The homozygous C group appeared to have a roughly flat growth rate between the heterozygotes and did not significantly vary with either heterozygote group (CC vs. C/T p=0.99,CC vs. T/C p=0.74)

[0080] The association between tumor size and age (days) are significantly different among these four genotypes (separating heterozygotes by their respective parental contributions) CC, C/T, T/C, TT (p=0.0162, ACOVA.) The most impressive growth phenotype was exhibited by homozygous T samples reaching an average of 7.8 cm before excision (FIG. 11A) Comparing the TT group against all non-TT subjects, the difference in lesion size, increased significantly with age (p=0.0001, ANCOVA). Thus in this study, TT lesions grew more rapidly than non-TT genotypes. In fact, each growth curveseparated by maternal and paternal genotypevaried independently and significantly from the TT samples (CC vs. TT: P<0.0001, CT vs. TT: P<0.0008, TC vs. TT: P=0.0025). Although genotype appears to have no effect at approximately 100 days, after three months, lesions begin to distinguish themselves suggesting distinct growth velocities. Furthermore, these data suggest parent of origin specific effects as those samples with identical genotypes but opposite parental contributions displayed statistically significant differences in growth curves. Namely, those lesions carrying the paternal T/maternal C genotype grew at approximately twice the rate as their paternal C/maternal T carrying counterparts p=0.05 (FIG. 10B) Furthermore, each heterozygote growth curve varied by age with high correlation of r squared above 0.85. However, it must be emphasized that sample size is relatively low (particularly in the C paternal/T maternal group) and the p-value just reached the threshold of significance.

TABLE-US-00002 TABLE 4 CLINICAL TABLE OF ULCERATION BY TT AND NON-TT GENOTYPE Retrospectively collected results with associated descriptive statistics: TT Lesions have a significantly higher associated odds ratio for ulceration. This proposed clinical test may be most useful in ruling out the chance of ulceration early in the disease course as sensitivity and negative predictive value are high. A larger prospective study is warranted. Ulcerated Not Ulcerated TT 6 3 Non TT 0 20 Sensitivity = 100% Odds Ratio = 76.1 Specificity = 86.96% 95% CI = 3.4-1676 Positive Predictive Value = 66.67% P = .006 Negative Predictive Value = 100%

[0081] TT Lesions Have a Significantly Higher Associated Odds Ratio for Ulceration. This proposed clinical test could be most useful in ruling out the chance of ulceration early in the disease course as sensitivity and negative predictive value are high. A larger prospective study is warranted.

[0082] Lastly, size is a highly significant clinical outcome when studying IH. However, of similar importance is ulceration. Once an IH ulcerates, it is usually painful for the patient and is given to bleeding which can be clinically significant. Ulceration is usually a marker of rapid disease progression and heralds an escalation of care. This can entail the institution of laser therapy, pharmacologic intervention or surgical excision. Not surprisingly, ulceration alone can prompt surgical treatment regardless of size or location of the lesion. To study the risk of ulceration an odds ratio calculation was performed comparing TT and non-TT lesions. The TT lesions had an odds ratio of 76.1 for ulceration p=0.006 (Table 4). Although this is a small sample cohort preliminary specificity was performed, sensitivity and positive and negative predictive value calculations (Table 4). These early results suggest the highest clinical usefulness of the proposed test in ruling out potential future ulceration. Although encouraging, these results will need to be corroborated prospectively in a larger cohort.

[0083] It will be apparent to those skilled in the art that various modifications, combinations and variations can be made in the present invention without departing from the spirit or scope of the invention. Specific embodiments, features and elements described herein may be modified, and/or combined in any suitable manner. Thus, it is intended that the present invention cover the modifications, combinations and variations of this invention provided they come within the scope of the appended claims and their equivalents.

[0084] The following references are hereby incorporated by reference in their entirety:

REFERENCES

[0085] 1. Mulliken J B, Glowacki J (1982) Hemangiomas and vascular malformations in infants and children: a classification based on endothelial characteristics. Plastic and reconstructive surgery 69: 412-422. PMlD: 7063565

[0086] 2. Mulliken J B, Glowacki J (1982) Classification of pediatric vascular lesions. Plastic and reconstructive surgery 70: 120-121. PMID: 7089103

[0087] 3. Glowacki J, Mulliken J B (1982) Mast cells in hemangiomas and vascular malformations. Pediatrics 70: 48-51.;PMID: 7088632

[0088] 4. North P E, Waner M, Brodsky M C (2002) Are infantile hemangioma of placental origin? Ophthalmology 109: 223-224, PMID, 11825799

[0089] 5. North P E, Waner M, Mizeracki A, Mrak R E, Nicholas R, et al. (2001) A unique microvascular phenotype shared by juvenile hemangiomas and human placenta. Archives of dermatology 137: 559-570. PMID: 11346333

[0090] 6. North P E, Waner M, Mizeracki A, Mihm M C J (2000) GLUT1: a newly discovered immunohistochemical marker for juvenile hemangiomas. Human pathology 31: 11-22. PMlD: 10665907

[0091] 7. Chang L C, Haggstrom A N, Drolet B A, Baselga E, Chamlin S L, et al. (2008) Growth Characteristics of Infantile Hemangiomas: Implications for Management. Pediatrics 122: 360-367. doi: 10.1542/peds. 2007-2767 PMID: 18676554

[0092] 8. Ritter M R, Dorrell M I, Edmonds J, Friedlander S F, Friedlander M (2002) Insulin-like growth factor 2 and potential regulators of hemangioma growth and involution identified by large-scale expression analysis. Proceedings of the National Academy of Sciences of the United States of America 99: 7455-7460, PMID: 12032304

[0093] 9. Yu Y, Wylie-Sears J, Boscolo E, Mulliken J B, Bischoff J (2004) Genomic imprinting of IGF2 is maintained in infantile hemangioma despite its high level of expression. Molecular medicine (Cambridge, Mass 10: 117-123,

[0094] 10. Weksberg R (2003) Beckwith-Wiedemann syndrome demonstrates a role for epigenetic control of normal development. Human molecular genetics 12: 61R-68. doi:10.1093/hmg/ddg067

[0095] 11. DeBaun M R, Niemitz E L, McNeil D E, Brandenburg S A, Lee M P, et al. (2002) Epigenetic alterations of H19 and LIT1 distinguish patients with Beckwith-Wiedemann syndrome with cancer and birth defects. American journal of human genetics 70: 604-611. PMID: 11813134

[0096] 12. Ferguson-Smith A C, Surani M A (2001) Imprinting and the epigenetic asymmetry between parental genomes. Science 293:1086-1089. PMID: 11498578

[0097] 13. Reik W, Walter J (2001) Genomic imprinting: parental influence on the genome. Nature reviews 2: 21-32. PMID: 11253064

[0098] 14. Chernukhin I, Shamsuddin S, Kang S Y, Bergstrom R, Kwon Y W, et al. (2007) Ctcf Interacts with and Recruits the Largest Subunit of Rna Polymerase li to Ctcf Target Sites Genome-Wide. Molecular and cellular biology.

[0099] 15. Thorvaldsen J L, Duran K L, Bartolomei M S (1998) Deletion of the H19 differentially methylated domain results in loss of imprinted expression of H19 and Igf2. Genes & development 12: 3693-3702. PMID: 9851976

[0100] 16. Schoenherr C J, Levorse J M, Tilghman S M (2003) CTCF maintains differential methylation at the Igf2/H19 locus. Nature genetics 33: 66-69. PMID: 12461525

[0101] 17. Florentino F P, Giordano A (2011) The tumor suppressor role of CTCF, J Cell Physiol. doi: 10.1002/jcp. 22780

[0102] 18. Kurukuti S, Tiwari V K, Tavoosidana G, Pugacheva E, Murrell A, et al. (2006) CTCF binding at the H19 imprinting control region mediates maternally inherited higher-order chromatin conformation to restrict enhancer access to Igf2. Proceedings of the National Academy of Sciences of the United States of America 103: 10684-10689. PMID: 16815976

[0103] 19. Grunau C, Clark S J, Rosenthal A (2001) Bisulfite genomic sequencing: systematic investigation of critical experimental parameters. Nucleic Acids Research 29: e65-e65. doi: 10.1093/nar/29.13.e65 PMID: 11433041

[0104] 20. Dupont J-M, Tost J, Jammes H, Gut IG (2004) De novo quantitative bisulfite sequencing using the pyrosequencing technology. Analytical Biochemistry 333: 119-127. doi: 10.1016/j.ab.2004.05.007 PMID: 15351288

[0105] 21. Tost J, Jammes H, Dupont J M, Buffat C, Robert B, et at (2006) Non-random, individual-specific methylation profiles are present at the sixth CTCF binding site in the human H19/IGF2 imprinting control region. Nucleic Acids Research 34: 5438-5448. PMID: 17012269

[0106] 22. Klenova E M, Morse H C, Ohlsson R, Lobanenkov V V (2002) The novel BORIS+CTCF gene family is uniquely involved in the epigenetics of normal biology and cancer. Seminars in cancer biology 12: 399-414. PMID: 12191639

[0107] 23. Loukinov D I, Pugacheva E, Vatolin S, Pack S D, Moon H, et al. (2002) BORIS, a novel male germ-linespecific protein associated with epigenetic reprogramming events, shares the same 11-zinc-finger domain with CTCF, the insulator protein involved in reading imprinting marks in the soma. Proceedings of the National Academy of Sciences of the United States of America 99: 6806-6811. PMID: 12011441

[0108] 24. Loukinov D, Ghochikyan A, Mkrtichyan M, Ichim T E, Lobanenkov V V, et al. (2006) Antitumor efficacy of DNA vaccination to the epigenetically acting tumor promoting transcription factor BORIS and CD80 molecular adjuvant. J Cell Biochem 98: 1037-1043, doi: 10.1002/jcb.20953 PMID: 16741971

[0109] 25. Jones T A, Ogunkolade B W, Szary J, Aarum J, Mumin M A, et al. (2011) Widespread expression of BORIS/CTCFL in normal and cancer cells. PLoS ONE 6: e22399. doi: 10.1371/journal.pone.0022399 PMID: 21811597

[0110] 26. Woloszynska-Read A, Zhang W, Yu J, Link P A, Mhawech-Fauceglia P, et al. (2011) Coordinated cancer germline antigen promoter and global DNA hypomethylation in ovarian cancer: association with the BORIS/CTCF expression ratio and advanced stage. Clin Cancer Res 17: 2170-2180. doi: 10.1158/1078-0432.CCR-10-2315 PMID: 21296871

[0111] 27. Renaud S, Pugacheva E M, Delgado M D, Braunschweig R, Abdullaev Z, et al. (2007) Expression of the CTCF-paralogous cancer-testis, gene, brother of the regulator of imprinted sites (BORIS), is regulated by three alternative promoters modulated by CpG methylation and by CTCF and p53 transcription factors, Nucleic Acids Research 35: 7372-7388, PMID: 17962299

[0112] 28. Burke L J, Zhang R, Bartkuhn M, Tiwari V K, Tavoosidana G, et al. (2005) CTCF binding and higher order chromatin structure of the H19 locus are maintained in mitotic chromatin. The EMBO journal 24: 3291-3300. PMID: 16107875

[0113] 29. Filippova G N, Qi C F, Ulmer J E, Moore J M, Ward M D, et al. (2002) Tumor-associated zinc finger mutations in the CTCF transcription factor selectively alter tts DNA-binding specificity. Cancer Research 62: 48-52. PMID: 11782357

[0114] 30. Zhang Y, Shields T, Crenshaw T, Hao Y, Moulton T, et al. (1993) Imprinting of human H19: allele-specific CpG methylation, loss of the active allele in Wilms tumor, and potential for somatic allele switching. American journal of human genetics 53: 113-124, PMID: 8391213

[0115] 31. Bell AC, Felsenfeld G (2000) Methylation of a CTCF-dependent boundary controls imprinted expression of the Igf2 gene. Nature 405: 482-485. PMlD: 10839546

[0116] 32. Kanduri C, Pant V, Loukinov D, Pugacheva E, Qi C F, et al. (2000) Functional association of CTCF with the insulator upstream of the H19 gene is parent of origin-specific and methylation-sensitive. Curr Biol 10: 853-856. PMID: 10899010

[0117] 33. Feinberg A P, Cui H, Ohlsson R (2002) DNA methylation and genomic imprinting: insights from cancer into epigenetic mechanisms. Seminars in cancel biology 12: 389-398. PMID: 12191638

[0118] 34. Bell A C, West A G, Felsenfeld C (1999) The protein CTCF is required for the enhancer blocking activity of vertebrate insulators. Cell 98: 387-396. PMID: 10458813

[0119] 35.Wolffe A P (2000) Transcriptional control: imprinting insulation. Curr Biol 10: R463-R465, PMID: 10873799

[0120] 36. Ulaner G A, Vu T H, Li T, Hu J F, Yao X M, et al. (2003) Loss of imprinting of IGF2 and H19 in osteosarcoma is accompanied by reciprocal methylation changes of a CTCF-binding site. Human molecular genetics 12: 535-549. PMID: 12588801

[0121] 37. Pant V, Kurukuti S, Pugacheva E, Shamsuddin S, Mariano P, et al. (2004) Mutation of a single CTCF target site within the H19 imprinting control region leads to loss of Igf2 imprinting and complex patterns of de novo methylation upon maternal inheritance. Molecular and cellular biology 24: 3497-3504.PMlD: 15060168

[0122] 38. Jelinic P, Stehle J C, Shaw P (2006) The testis-specific factor CTCFL cooperates with the protein methyltransferase PRMT7 in H19 imprinting control region methylation. PLoS biology 4: e355. PMlD: 17048991

[0123] 39. Adkins R M, Spines C, Morrison J C, Hill J B, Watson E M, et al. (2010) Association of birth weight with polymorphisms in the IGF2, H19, and IGF2R genes. Pediatr Res 68: 429-434. doi: 10.1203/PDR. 0b013e3181f1ca99 PMID: 20639793

[0124] 40. Pearce C L, Doherty J A, Van Den Berg D J, Moysich K, Hsu C, et al. (2011) Genetic variation in insulinlike growth factor 2 may play a role in ovarian cancer risk. Human molecular genetics 20: 2263-2272. doi: 10.1093/hmg/ddr087 PMID: 21422097

[0125] 41. Neuhausen S L, Brummel S, Ding V C, Steele L, Nathanson K L, et al. (2011) Genetic Variation in IGF2 and HTRA1 and Breast Cancer Risk among BRCA1 and BRCA2 Carriers. Cancer Epidemiology Biomarkers & Prevention 20: 1690-1702. doi: 10.1158/1055-9985,EPI-10-1338 PMID: 21708937

[0126] 42. Coolen M W, Statham A L, Qu W, Campbell M J, Henders A K, et al. (2011) Impact of the genome on the epigenome is manifested in DNA methylation patterns of imprinted regions in monozygotic and dizygotic twins. PLoS ONE 6: e25590. doi: 10.1371/journal.pone.0025590 PMID: 21991322

[0127] 43. Sleutels F, Soochit W, Bartkuhn M, Heath H, Dienstbach S. et al. (2012) The male germ cell gene regulator CTCFL is functionally different from CTCF and binds CTCF-like consensus sites in a nucleosome composition-dependent manner. Epigenetics Chromatin 5: 8. doi: 10.1188/1756-8935-5-8 PMID: 22709888