Abstract
Background
Loss-of-function mutations in the filaggrin gene (FLG), which encodes an epidermal protein crucial for the formation of a functional skin barrier, have been identified as a major predisposing factor in the etiopathogenesis of atopic dermatitis (AD). Recent reports of relatively low frequencies of FLG-null mutations among specific ethnic groups with AD necessitated analysis of the epigenetic regulation which may control FLG expression without altering its DNA sequence.
Methods
Quantitative polymerase chain reaction was performed to determine the restoration of FLG mRNA expression in normal human epidermal keratinocyte (NHEK) cells after treatment with epigenetic modulating agents. Bisulfite genomic sequencing and pyrosequencing analyses of the FLG promoter region were conducted to identify the citical CpG sites relevant to FLG expression. We performed small-scale pilot study for epidermal tissues obtained from Korean patients with severe AD.
Results
We here show that DNA methylation in the FLG with non-CpG island promoter is responsible for the transcriptional regulation of FLG in undifferentiated NHEK cells. The methylation frequencies in a single CpG site of the FLG promoter were significantly higher in lesional epidermis than those in matched nonlesional epidermis of subjects with severe AD.
Profilaggrin, filaggrin (filament-aggregating protein), and their proteolyzed amino acids are multifunctional protein components contributing to formation of the intact stratum corneum barrier1. According to the “outside–inside” view, atopic dermatitis (AD) is a common chronic skin disease that is characterized by abnormal permeability barrier owing to the deficiency of filaggrin2. Since two prevalent loss-of-function (LOF) mutations (p.R501X and c.2282del4) of the filaggrin gene (FLG) were identified in the European populations with ichthyosis vulgaris (IV) and atopic eczema34, a comprehensive analytical strategy for screening of FLG variants has been implemented in the cohorts with not only IV or AD but also secondary allergic diseases5. Although various FLG-null mutations have been revealed in the selected Asian populations (including Japanese67, Chinese89, Singaporean Chinese10, Taiwanese11, and Korean1213), the prevalence and mutation rates of FLG-null alleles are different from European cohorts with IV or AD. Additionally, FLG LOF variants were found to be less common in patients of African ancestry with AD compared with European and Asian patients1415. A recent study revealed FLG mutations in 22.2% of African Americans with IV/AD16, whereas no FLG LOF mutations were detected in South African amaXhosa patients with AD17. Besides the genetic background predisposing to skin barrier dysfunction, non-genetic alterations of FLG seem to be correlated with AD pathogenesis in non-European populations, which have shown relatively low frequencies of FLG-null mutations.
An epigenetic trait is defined as a stably heritable phenotype resulting from changes in a chromosome without alterations in the DNA sequence18. Epigenetic mechanisms including DNA methylation, histone modifications, and non-coding RNAs have been implicated in a huge variety of human diseases correlating with gene–environment interactions1920. Methylation of the cytosine base within a CpG dinucleotide context in gene promoters is generally associated with transcriptional repression2122. Aberrant DNA methylation and subsequent silencing of a tumor suppressor gene have also been demonstrated in skin tumors, such as basal cell carcinoma, and in autoimmunityrelated skin disorders23. Recently, epigenetic alterations have been proposed as an alternative mechanism of predisposition to AD in people with wild-type FLG. Ziyab et al.24, for example, reported that DNA methylation at the single CpG site in the FLG gene body showed a significant interaction with FLG LOF variants on the risk for eczema. Nonetheless, epigenetic factors that may account for the deficiency of filaggrin and its degradation products in AD remain elusive. We hypothesized that epigenetic changes at the FLG gene promoter might contribute to modulation of FLG expression in human epidermal keratinocytes and might be associated with one of the etiopathogenic mechanisms that trigger the deficit of filaggrin in the affected epidermis of AD patients.
This study was conducted according to the Declaration of Helsinki Principles and was approved by the Institutional Review Board of the Chung-Ang University Hospital (IRB no. C2015258 [1716]). All patients gave their written informed consent to participate in the study.
Normal human epidermal keratinocyte (NHEK) cells were obtained from the PromoCell (Heidelberg, Germany). NHEK cells were grown in Keratinocyte Growth Medium 2 (PromoCell) supplemented with Keratinocyte Growth Medium 2 SupplementMix (PromoCell). All cells were maintained at 37℃ in a humidified incubator at 5% CO2. We treated the undifferentiated NHEK cells with 2′-deoxy-5-azacytidine (DAC) for 72 hours and with trichostatin A (TSA) for 24 hours. DAC and TSA were obtained from Sigma-Aldrich (St. Louis, MO, USA).
Punch biopsy specimens were collected from the paired non-lesional and lesional skin of patients with severe AD (n=10). Epidermal samples were isolated under a dissecting microscope. Genomic DNA extraction was performed using the QIAamp DNA Mini Kit (Qiagen, Hilden, Germany).
As described in detail previously25, total RNA was isolated from the cell lines using the RNeasy Mini Kit (Qiagen) following the manufacturer's instructions. cDNA synthesis for real-time two-step reverse transcription polymerase chain reaction (PCR) was performed on 1 µg of total RNA using the QuantiTect Reverse Transcription Kit (Qiagen), and then 1 µl of diluted cDNA was utilized for cycling reactions using the Rotor-Gene SYBR Green PCR Kit (Qiagen). The amplification and quantitative analysis were performed on the Rotor-Gene Q 5plex HRM system (Qiagen). Thermal cycling was conducted using the default conditions of the Rotor-Gene Q Series Software (Qiagen), which consisted of 5 minutes at 95℃ followed by 40 rounds of 5 seconds at 95℃ and 10 seconds at 60℃. The transcript level measured was normalized to the results of a QuantiTect Primer Assay (Qiagen) for ACTB (encoding β-actin), which served as an internal control.
As described in detail previously25, a 928-bp fragment harboring four CpG sites located at the 5′-end region of FLG (NC_ 000001.11), a 5′ NsiI site, and a 3′ HindIII site was synthesized and cloned into the pCpGfree-basic-Lucia reporter plasmid (InvivoGen, San Diego, CA, USA) that codes for a secreted coelenterazine-utilizing variant of luciferase. Briefly, 5 µg of the Lucia reporter plasmid was methylated using 12 U of M.SssI CpG methyltransferase (New England BioLabs, Ipswich, MA, USA) in vitro at 37℃ for 20 hours. After purification with the DNA Clean & Concentrato™-5 Kit (Zymo Research, Irvine, CA, USA), methylation of the plasmid was verified by bisulfite pyrosequencing of the selected CpG units in the FLG promoter (Supplementary Fig. 1). Next, HEK293T cells were transiently transfected with unmethylated or methylated reporter plasmids by means of the jetPEI reagent (Polyplustransfection, New York, NY, USA). Luciferase activity of the Lucia reporter in the supernatant of the cells at 72 hours after the transfection was measured using the QUANTI-Luc assay reagent (InvivoGen). In each experiment, the cells were cotransfected with the pCMV-CLuc 2 Control Plasmid (New England BioLabs), which encodes a secreted variant of Cypridina luciferase, as a normalization control. Luminescence induced by Cypridina luciferase in the cell supernatant was determined using the BioLux Cypridina Luciferase Assay Kit (New England BioLabs). Reporter activity was normalized by calculating the ratio of Lucia to CLuc activities.
For bisulfite genomic sequencing analysis, as described in detail previously25, genomic DNA was extracted from the cell lines using the QIAamp DNA Mini Kit. Bisulfite modification of the genomic DNA was performed using the EZ DNA Methylation-Lightning™ Kit (Zymo Research). The bisulfite-converted genomic DNA was amplified using primer sets specific to the proximal promoter region of FLG (positions −738 to −162 upstream of the transcription start site (TSS), which is denoted as +1 in Supplementary Fig. 2) containing four CpG units. The sequences of the PCR primers used were as follows: forward (5′-AGAAGGAAGAGTATGTGGAATATG-3′) and reverse (5′- CAACAACCTATATTTACTTCCCAAC-3′). Cycling conditions were as follows: initial denaturation at 94℃ for 10 minutes; 45~50 cycles of denaturation at 94℃ for 30 seconds, annealing at 56℃ for 30 seconds, and extension at 72℃ for 30 seconds; and the final extension at 72℃ for 10 minutes. PCR products were purified and subcloned into the pGEM-T Easy Vector (Promega, Madison, WI, USA) for subsequent sequencing analysis. The nucleotide sequences of 12~15 independent clones were analyzed.
Furthermore, the bisulfite-modified genomic DNA obtained from normal individuals and AD patients was amplified using primer sets specific to two key CpG sites within the FLG promoter. The sequences of the PCR primers used were as follows: CpG 4 site, forward (5′-GAAGGAAGAGTATGTGGAATATGT-3′), 5′-biotinylated reverse (5′-CACTAAAAACATAAATTTAATTAAACAAAACTC-3′), and sequencing (5′- TGTGGAATATGTTTTTGATGTTA-3′); CpG 2 site, forward (5′-AGATGGAATATATAGATTAAAAAAGAATA-3′), 5′-biotinylated reverse (5′-TAAATACACTTCCTAATCCTTATCTCT-3′), and sequencing (5′-TTTTTATTATAAATTGAATTTTAAG-3′). Pyrosequencing reactions were performed using the PyroMark Gold Q96 Reagents (Qia gen) and quantitative analysis was conducted on the PyroMark Q96 ID platform (Qiagen).
Data were presented as mean±standard error of the mean and analyzed by Student's t-test. To determine significant differences in the quantitative methylation data generated by pyrosequencing measurements, comparisons between nonlesional and lesional tissue samples were made using the Wilcoxon signed-rank test. Differences with a p<0.05 were considered statistically significant. Statistical analyses were conducted using the IBM SPSS Statistics ver. 20 software (IBM Corp., Armonk, NY, USA).
To test whether epigenetic regulation of FLG takes place in NHEK cells, we first examined the restoration of FLG expression in undifferentiated NHEK cells after treatment with either the DNA methyltransferase inhibitor DAC or the histone deacetylase inhibitor TSA (Fig. 1A). The DAC-induced DNA demethylation and TSA-mediated histone acetylation triggered a significant increase in FLG mRNA levels in undifferentiated NHEK cells. These data implied that the loss of FLG expression in undifferentiated NHEK cells correlated with epigenetic silencing involving both DNA methylation and histone modifications.
Next, to determine whether DNA methylation in the 5′-end region of FLG was important for modulating the transcriptional activity of the core promoter, we conducted an in vitro methylated reporter assay (Fig. 1B). We synthesized an FLG fragment containing the non-CpG island (CGI) promoter and four CpG dinucleotide units therein (Fig. 1C). In vitro methylation was performed after cloning this fragment into a CpG-free luciferase plasmid to assess the role of CpG methylation in FLG transcription. HEK293T cells transfected with the methylated FLG reporter plasmid showed markedly decreased luciferase activity as compared with the cells transfected with the unmethylated FLG plasmid. The magnitude of the reduction between unmethylated control and methylated plasmid was approximately 90%, suggesting that DNA methylation directly led to transcriptional silencing of FLG with non-CGI promoter.
To identify the critical CpG sites important for regulating FLG expression in undifferentiated NHEK cells, we performed bisulfite genomic sequencing of the proximal promoter region of FLG. We determined the methylation profiles of the four CpG dinucleotides (i.e., CpG 1 to 4) located in the FLG promoter, assigning the first CpG dinucleotide to a site 185 bp upstream of FLG TSS (CpG 1). As shown in Fig. 1C, the methylation profile of the FLG promoter revealed a substantial difference in methylation frequency at both the CpG 2 and CpG 4 loci (−410 and −702 bp upstream of the FLG TSS, respectively) in the DAC-treated NHEK cells, when compared to control cells. These findings indicate that the site-specific CpG dinucleotides within the FLG promoter might be preferentially affected by a DNA demethylating agent in undifferentiated NHEK cells.
To evaluate whether the epigenotyping for FLG by using of the putative DNA methylation marker (i.e., CpG 2 and CpG 4 dinucleotide units within the FLG promoter) was applicable to the epidermal tissues from AD patients, we performed a small-scale pilot study comparing methylation differences of these CpG loci between non-lesional and lesional epidermis obtained from Korean patients with severe AD (Eczema Area Severity Index>15; n=10; Table 1). After the triplicate experiments of bisulfite pyrosequencing for each pairs of matched samples, we observed that CpG 2 sites were characterized by overall higher methylation frequencies in lesional tissues than in matched nonlesional tissues (Fig. 2A). Although the differential methylation status at the CpG 4 unit was detected during the bisulfite genomic sequencing of FLG in undifferentiated NHEK cells, no significant correlation was observed in the quantitative methylation analysis between nonlesional and lesional epidermal tissue samples of AD subjects (Fig. 2B). Nonetheless, a statistically significant association of methylation frequencies was shown by a single CpG unit, i.e., CpG 2, within the FLG promoter (p=0.0343) in the skin biopsy samples categorized into non-lesional and lesional samples from the AD patients (Fig. 2B).
Epigenetic alterations in AD, including histone modifications and miRNAs, have been discovered in a variety of genes except FLG2627. The 5′-end region of FLG is devoid of any defined CGI (length ≥200 bp, Guanine-Cytosine content ≥50%, and observed-to-expected CpG ratio [ObsCpG/ExpCpG] ≥0.6) (Supplementary Fig. 2)28. The potential involvement of variable methylation in the non-CGI promoter of FLG has been mostly overlooked. In contrast to genes with CGIs at their TSSs, substantial fluctuations occur in the promoter methylation levels of genes that are CpG-poor at the TSS22. Because CpG-poor regulatory regions tend to acquire a low methylation state when occupied by transcription factors2930, DNA methylation at non-CGI promoters have high information content about the ongoing activity of transcription factors21. Several studies have shown that no significant promoter methylation is associated with FLG expression in AD cases2431 and in buccal cells32. Herein, we focused on the CpG sites that were not examined in previous epigenetic analyses, because the methylation array-based studies on FLG in AD subjects could not have included all the differentially methylated CpG loci in the non-CGI promoter of FLG. Noticeably, non-CGI promoters, which are generally hypermethylated, remain transcriptionally active regardless of their methylation state3334, but our findings demonstrate that DNA methylation of only four CpG sites directly suppress the non-CGI promoter of FLG. Moreover, significant differences were observed in the methylation frequency of the single CpG site between non-lesional and lesional epidermis obtained from Korean patients with severe AD, suggesting that the possible epigenetic alterations in vivo might depend on the methylation status at this site-specific CpG unit within the FLG promoter.
It has been previously found that the levels of filaggrin breakdown products, such as trans-urocanic acid and pyrrolidone-5-carboxylic acid, are significantly reduced in patients with AD without FLG-null mutations1735. Our results support the insight that epigenetic alterations may account for discordant phenotypes of patients with a skin disease who are wild-type on the prevalent mutations of the disease36. Because a given cytosine can either be completely methylated or unmethylated, ‘variable methylation’ is the outcome of averaging these binary states22. Hence, what is measured for a given epidermal sample at a given CpG site, i.e., the CpG 2 unit within the FLG promoter, is the percentage of epidermal cells that are methylated. Considering the significant differences in the methylation frequencies among the dermatologist-diagnosed AD samples, this key CpG 2 unit in the FLG promoter can serve as an AD-associated DNA methylation biomarker. Indeed, a single CpG site showing differential methylation between different disease states is known as a methylation variable position (MVP), which can be considered as the epigenetic equivalent of a single nucleotide polymorphism3738. This might suggest that the methylation frequency at any given MVP may serve as a novel epigenetic signature for physiological and pathologic status3940. Thus, our observations also suggest that the CpG 2 unit within the FLG promoter might serve as an MVP candidate of filaggrindeficient skin phenotype.
In conclusion, we provided evidence for one of the mechanisms via which significant contributors to filaggrin deficiency work in the lesional skin of patients with AD. To extrapolate the plausible correlation between epigenotype of FLG and genetic background of affected individuals and controls, future studies should involve a larger AD cohort where the genotyping for the FLG-null mutations is performed beforehand. Furthermore, our finding may have more promising clinical relevance in conjunction with FLG-genotyping studies, which include diseases of various severity.
ACKNOWLEDGMENT
The authors thank Chung-Ang University Hospital Biomedical Research Institute for kindly offering research equipment.
This work was supported by a grant (714001-07) from the Center for Industrialization of Natural Nutraceuticals through the Agriculture, Food and Rural Affairs Research Center Support Program, Ministry of Agriculture, Food and Rural Affairs, Republic of Korea.
References
1. Brown SJ, McLean WH. One remarkable molecule: filaggrin. J Invest Dermatol. 2012; 132(3 Pt 2):751–762. PMID: 22158554.
2. Elias PM, Steinhoff M. “Outside-to-inside” (and now back to “outside”) pathogenic mechanisms in atopic dermatitis. J Invest Dermatol. 2008; 128:1067–1070. PMID: 18408746.
3. Smith FJ, Irvine AD, Terron-Kwiatkowski A, Sandilands A, Campbell LE, Zhao Y, et al. Loss-of-function mutations in the gene encoding filaggrin cause ichthyosis vulgaris. Nat Genet. 2006; 38:337–342. PMID: 16444271.
4. Palmer CN, Irvine AD, Terron-Kwiatkowski A, Zhao Y, Liao H, Lee SP, et al. Common loss-of-function variants of the epidermal barrier protein filaggrin are a major predisposing factor for atopic dermatitis. Nat Genet. 2006; 38:441–446. PMID: 16550169.
5. Sandilands A, Terron-Kwiatkowski A, Hull PR, O'Regan GM, Clayton TH, Watson RM, et al. Comprehensive analysis of the gene encoding filaggrin uncovers prevalent and rare mutations in ichthyosis vulgaris and atopic eczema. Nat Genet. 2007; 39:650–654. PMID: 17417636.
6. Nomura T, Sandilands A, Akiyama M, Liao H, Evans AT, Sakai K, et al. Unique mutations in the filaggrin gene in Japanese patients with ichthyosis vulgaris and atopic dermatitis. J Allergy Clin Immunol. 2007; 119:434–440. PMID: 17291859.
7. Hamada T, Sandilands A, Fukuda S, Sakaguchi S, Ohyama B, Yasumoto S, et al. De novo occurrence of the filaggrin mutation p.R501X with prevalent mutation c.3321delA in a Japanese family with ichthyosis vulgaris complicated by atopic dermatitis. J Invest Dermatol. 2008; 128:1323–1325. PMID: 18007582.
8. Zhang X, Liu S, Chen X, Zhou B, Liu D, Lei G, et al. Novel and recurrent mutations in the filaggrin gene in Chinese patients with ichthyosis vulgaris. Br J Dermatol. 2010; 163:63–69. PMID: 20222934.
9. Zhang H, Guo Y, Wang W, Shi M, Chen X, Yao Z. Mutations in the filaggrin gene in Han Chinese patients with atopic dermatitis. Allergy. 2011; 66:420–427. PMID: 21039602.
10. Chen H, Common JE, Haines RL, Balakrishnan A, Brown SJ, Goh CS, et al. Wide spectrum of filaggrin-null mutations in atopic dermatitis highlights differences between Singaporean Chinese and European populations. Br J Dermatol. 2011; 165:106–114. PMID: 21428977.
11. Hsu CK, Akiyama M, Nemoto-Hasebe I, Nomura T, Sandilands A, Chao SC, et al. Analysis of Taiwanese ichthyosis vulgaris families further demonstrates differences in FLG mutations between European and Asian populations. Br J Dermatol. 2009; 161:448–451. PMID: 19416262.
12. Kim EJ, Jeong MS, Li K, Park MK, Lee MK, Yoon Y, et al. Genetic polymorphism of FLG in Korean ichthyosis vulgaris patients. Ann Dermatol. 2011; 23:170–176. PMID: 21747615.
13. Park J, Jekarl DW, Kim Y, Kim J, Kim M, Park YM. Novel FLG null mutations in Korean patients with atopic dermatitis and comparison of the mutational spectra in Asian populations. J Dermatol. 2015; 42:867–873. PMID: 25997159.
14. Winge MC, Bilcha KD, Liedén A, Shibeshi D, Sandilands A, Wahlgren CF, et al. Novel filaggrin mutation but no other loss-of-function variants found in Ethiopian patients with atopic dermatitis. Br J Dermatol. 2011; 165:1074–1080. PMID: 21692775.
15. Margolis DJ, Gupta J, Apter AJ, Hoffstad O, Papadopoulos M, Rebbeck TR, et al. Exome sequencing of filaggrin and related genes in African-American children with atopic dermatitis. J Invest Dermatol. 2014; 134:2272–2274. PMID: 24608987.
16. Polcari I, Becker L, Stein SL, Smith MS, Paller AS. Filaggrin gene mutations in African Americans with both ichthyosis vulgaris and atopic dermatitis. Pediatr Dermatol. 2014; 31:489–492. PMID: 24920311.
17. Thawer-Esmail F, Jakasa I, Todd G, Wen Y, Brown SJ, Kroboth K, et al. South African amaXhosa patients with atopic dermatitis have decreased levels of filaggrin breakdown products but no loss-of-function mutations in filaggrin. J Allergy Clin Immunol. 2014; 133:280–282. 282.e1–282.e2. PMID: 24369804.
18. Berger SL, Kouzarides T, Shiekhattar R, Shilatifard A. An operational definition of epigenetics. Genes Dev. 2009; 23:781–783. PMID: 19339683.
19. Jaenisch R, Bird A. Epigenetic regulation of gene expression: how the genome integrates intrinsic and environmental signals. Nat Genet. 2003; 33:245–254. PMID: 12610534.
20. Botchkarev VA, Gdula MR, Mardaryev AN, Sharov AA, Fessing MY. Epigenetic regulation of gene expression in keratinocytes. J Invest Dermatol. 2012; 132:2505–2521. PMID: 22763788.
21. Schübeler D. Function and information content of DNA methylation. Nature. 2015; 517:321–326. PMID: 25592537.
22. Jones PA. Functions of DNA methylation: islands, start sites, gene bodies and beyond. Nat Rev Genet. 2012; 13:484–492. PMID: 22641018.
23. Li Y, Sawalha AH, Lu Q. Aberrant DNA methylation in skin diseases. J Dermatol Sci. 2009; 54:143–149. PMID: 19395242.
24. Ziyab AH, Karmaus W, Holloway JW, Zhang H, Ewart S, Arshad SH. DNA methylation of the filaggrin gene adds to the risk of eczema associated with loss-of-function variants. J Eur Acad Dermatol Venereol. 2013; 27:e420–e423. PMID: 23003573.
25. Lee J, Han JH, Jang A, Kim JW, Hong SA, Myung SC. DNA methylation-mediated downregulation of DEFB1 in prostate cancer cells. PLoS One. 2016; 11:e0166664. PMID: 27835705.
26. Alaskhar Alhamwe B, Khalaila R, Wolf J, von Bülow V, Harb H, Alhamdan F, et al. Histone modifications and their role in epigenetics of atopy and allergic diseases. Allergy Asthma Clin Immunol. 2018; 14:39. PMID: 29796022.
28. Gardiner-Garden M, Frommer M. CpG islands in vertebrate genomes. J Mol Biol. 1987; 196:261–282. PMID: 3656447.
29. Stadler MB, Murr R, Burger L, Ivanek R, Lienert F, Schöler A, et al. DNA-binding factors shape the mouse methylome at distal regulatory regions. Nature. 2011; 480:490–495. PMID: 22170606.
30. Ziller MJ, Gu H, Müller F, Donaghey J, Tsai LT, Kohlbacher O, et al. Charting a dynamic DNA methylation landscape of the human genome. Nature. 2013; 500:477–481. PMID: 23925113.
31. Rodríguez E, Baurecht H, Wahn AF, Kretschmer A, Hotze M, Zeilinger S, et al. An integrated epigenetic and transcriptomic analysis reveals distinct tissue-specific patterns of DNA methylation associated with atopic dermatitis. J Invest Dermatol. 2014; 134:1873–1883. PMID: 24739813.
32. Tan HT, Ellis JA, Koplin JJ, Martino D, Dang TD, Suaini N, et al. HealthNuts Study Investigators. Methylation of the filaggrin gene promoter does not affect gene expression and allergy. Pediatr Allergy Immunol. 2014; 25:608–610. PMID: 24912553.
33. Weber M, Hellmann I, Stadler MB, Ramos L, Pääbo S, Rebhan M, et al. Distribution, silencing potential and evolutionary impact of promoter DNA methylation in the human genome. Nat Genet. 2007; 39:457–466. PMID: 17334365.
34. Meissner A, Mikkelsen TS, Gu H, Wernig M, Hanna J, Sivachenko A, et al. Genome-scale DNA methylation maps of pluripotent and differentiated cells. Nature. 2008; 454:766–770. PMID: 18600261.
35. Kezic S, O'Regan GM, Yau N, Sandilands A, Chen H, Campbell LE, et al. Levels of filaggrin degradation products are influenced by both filaggrin genotype and atopic dermatitis severity. Allergy. 2011; 66:934–940. PMID: 21261659.
36. Bataille V, Lens M, Spector TD. The use of the twin model to investigate the genetics and epigenetics of skin diseases with genomic, transcriptomic and methylation data. J Eur Acad Dermatol Venereol. 2012; 26:1067–1073. PMID: 22243446.
37. Tost J, Dunker J, Gut IG. Analysis and quantification of multiple methylation variable positions in CpG islands by Pyrosequencing. Biotechniques. 2003; 35:152–156. PMID: 12866415.
38. Rakyan VK, Down TA, Balding DJ, Beck S. Epigenome-wide association studies for common human diseases. Nat Rev Genet. 2011; 12:529–541. PMID: 21747404.
39. Rakyan VK, Hildmann T, Novik KL, Lewin J, Tost J, Cox AV, et al. DNA methylation profiling of the human major histocompatibility complex: a pilot study for the human epigenome project. PLoS Biol. 2004; 2:e405. PMID: 15550986.
40. Huang Q, Huang JF, Zhang B, Baum L, Fu WL. Methylation variable position profiles of hMLH1 promoter CpG islands in human sporadic colorectal carcinoma. Diagn Mol Pathol. 2012; 21:24–33. PMID: 22306672.
SUPPLEMENTARY MATERIALS
Supplementary data can be found via http://anndermatol.org/src/sm/ad-32-122-s001.pdf.