Phenotypic and functional translation of IL33 genetics in asthma.

BACKGROUND: Asthma is a complex disease with multiple phenotypes that may differ in disease pathobiology and treatment response. IL33 single nucleotide polymorphisms (SNPs) have been reproducibly associated with asthma. IL33 levels are elevated in sputum and bronchial biopsies of patients with asthma. The functional consequences of IL33 asthma SNPs remain unknown. OBJECTIVE: This study sought to determine whether IL33 SNPs associate with asthma-related phenotypes and with IL33 expression in lung or bronchial epithelium. This study investigated the effect of increased IL33 expression on human bronchial epithelial cell (HBEC) function. METHODS: Association between IL33 SNPs (Chr9: 5,815,786-6,657,983) and asthma phenotypes (Lifelines/DAG [Dutch Asthma GWAS]/GASP [Genetics of Asthma Severity & Phenotypes] cohorts) and between SNPs and expression (lung tissue, bronchial brushes, HBECs) was done using regression modeling. Lentiviral overexpression was used to study IL33 effects on HBECs. RESULTS: We found that 161 SNPs spanning the IL33 region associated with 1 or more asthma phenotypes after correction for multiple testing. We report a main independent signal tagged by rs992969 associating with blood eosinophil levels, asthma, and eosinophilic asthma. A second, independent signal tagged by rs4008366 presented modest association with eosinophilic asthma. Neither signal associated with FEV1, FEV1/forced vital capacity, atopy, and age of asthma onset. The 2 IL33 signals are expression quantitative loci in bronchial brushes and cultured HBECs, but not in lung tissue. IL33 overexpression in vitro resulted in reduced viability and reactive oxygen species-capturing of HBECs, without influencing epithelial cell count, metabolic activity, or barrier function. CONCLUSIONS: We identify IL33 as an epithelial susceptibility gene for eosinophilia and asthma, provide mechanistic insight, and implicate targeting of the IL33 pathway specifically in eosinophilic asthma.

Phenotypic and functional translation of IL1RL1 locus polymorphisms in lung tissue and asthmatic airway epithelium.

The IL1RL1 (ST2) gene locus is robustly associated with asthma; however, the contribution of single nucleotide polymorphisms (SNPs) in this locus to specific asthma subtypes and the functional mechanisms underlying these associations remain to be defined. We tested for association between IL1RL1 region SNPs and characteristics of asthma as defined by clinical and immunological measures and addressed functional effects of these genetic variants in lung tissue and airway epithelium. Utilizing 4 independent cohorts (Lifelines, Dutch Asthma GWAS [DAG], Genetics of Asthma Severity and Phenotypes [GASP], and Manchester Asthma and Allergy Study [MAAS]) and resequencing data, we identified 3 key signals associated with asthma features. Investigations in lung tissue and primary bronchial epithelial cells identified context-dependent relationships between the signals and IL1RL1 mRNA and soluble protein expression. This was also observed for asthma-associated IL1RL1 nonsynonymous coding TIR domain SNPs. Bronchial epithelial cell cultures from asthma patients, exposed to exacerbation-relevant stimulations, revealed modulatory effects for all 4 signals on IL1RL1 mRNA and/or protein expression, suggesting SNP-environment interactions. The IL1RL1 TIR signaling domain haplotype affected IL-33-driven NF-κB signaling, while not interfering with TLR signaling. In summary, we identify that IL1RL1 genetic signals potentially contribute to severe and eosinophilic phenotypes in asthma, as well as provide initial mechanistic insight, including genetic regulation of IL1RL1 isoform expression and receptor signaling.

Duration of breastfeeding is associated with leptin (LEP) DNA methylation profiles and BMI in 10-year-old children.

BACKGROUND: Breastfeeding is protective against many long-term diseases, yet the mechanisms involved are unknown. Leptin gene (LEP) is reported to be associated with body mass index (BMI). On the other hand, breastfeeding duration has been found to be associated with DNA methylation (DNAm) of the LEP gene. Therefore, epigenetic regulation of LEP may represent the mechanism underlying the protective effect of breastfeeding duration against obesity. METHODS: In the Isle of Wight Birth Cohort, peripheral blood DNAm at 23 cytosine-phosphate-guanine sites (CpGs) in the LEP locus in 10-year-old (n = 297) samples and 16 CpGs in 18-year-old (n = 305) samples, were generated using the Illumina Infinium MethylationEPIC and HumanMethylation450 Beadchips respectively and tested for association with breastfeeding duration (total and exclusive) using linear regression. To explore the association between breastfeeding durations and genome-wide DNAm, epigenome-wide association studies (EWASs) and differential methylation region (DMR) analyses were performed. BMI trajectories spanning the first 18 years of life were used as the outcome to test the association with breastfeeding duration (exposure) using multi-nominal logistic regression. Mediation analysis was performed for significant CpG sites. RESULTS: Both total and exclusive breastfeeding duration were associated with DNAm at four LEP CpG sites at 10 years (P value < 0.05), and not at 18 years. Though no association was observed between breastfeeding duration and genome-wide DNAm, DMR analyses identified five significant differentially methylated regions (Sidak adjusted P value < 0.05). Breastfeeding duration was also associated with the early transient overweight trajectory. Furthermore, DNAm of LEP was associated with this trajectory at one CpG site and early persistent obesity at another, though mediation analysis was not significant. CONCLUSIONS: Breastfeeding duration is associated with LEP methylation at age 10 years and BMI trajectory. LEP DNAm is also significantly associated with BMI trajectories throughout childhood, though sample sizes were small. However, mediation analysis did not demonstrate that DNAm of LEP explained the protective effect of breastfeeding against childhood obesity.

Epigenome-wide association study of asthma and wheeze characterizes loci within HK1.

BACKGROUND: To identify novel epigenetic markers of adolescent asthma and replicate findings in an independent cohort, then explore whether such markers are detectable at birth, predictive of early-life wheeze, and associated with gene expression in cord blood. METHODS: We performed epigenome-wide screening with recursive random forest feature selection and internal validation in the IOW birth cohort. We then tested whether we could replicate these findings in the independent cohort ALSPAC and followed-up our top finding with children of the IOW cohort. RESULTS: We identified 10 CpG sites associated with adolescent asthma at a 5% false discovery rate (IOW, n = 370), five of which exhibited evidence of associations in the replication study (ALSPAC, n = 720). One site, cg16658191, within HK1 displayed particularly strong associations after cellular heterogeneity adjustments in both cohorts (ORIOW = 0.17, 95% CI 0.04-0.57) (ORALSPAC = 0.57, 95% CI 0.38-0.87). Additionally, higher expression of HK1 (OR = 3.81, 95% CI 1.41-11.77) in cord blood was predictive of wheezing in infancy (n = 82). CONCLUSION: We identified novel associations between asthma and wheeze with methylation at cg16658191 and the expression of HK1, which may serve as markers of, predictors of, and potentially etiologic factors involved in asthma and early life wheeze.

Moderate-to-severe asthma in individuals of European ancestry: a genome-wide association study.

BACKGROUND: Few genetic studies that focus on moderate-to-severe asthma exist. We aimed to identity novel genetic variants associated with moderate-to-severe asthma, see whether previously identified genetic variants for all types of asthma contribute to moderate-to-severe asthma, and provide novel mechanistic insights using expression analyses in patients with asthma. METHODS: In this genome-wide association study, we used a two-stage case-control design. In stage 1, we genotyped patient-level data from two UK cohorts (the Genetics of Asthma Severity and Phenotypes [GASP] initiative and the Unbiased BIOmarkers in PREDiction of respiratory disease outcomes [U-BIOPRED] project) and used data from the UK Biobank to collect patient-level genomic data for cases and controls of European ancestry in a 1:5 ratio. Cases were defined as having moderate-to-severe asthma if they were taking appropriate medication or had been diagnosed by a doctor. Controls were defined as not having asthma, rhinitis, eczema, allergy, emphysema, or chronic bronchitis as diagnosed by a doctor. For stage 2, an independent cohort of cases and controls (1:5) was selected from the UK Biobank only, with no overlap with stage 1 samples. In stage 1 we undertook a genome-wide association study of moderate-to-severe asthma, and in stage 2 we followed up independent variants that reached the significance threshold of p less than 1 × 10-6 in stage 1. We set genome-wide significance at p less than 5 × 10-8. For novel signals, we investigated their effect on all types of asthma (mild, moderate, and severe). For all signals meeting genome-wide significance, we investigated their effect on gene expression in patients with asthma and controls. FINDINGS: We included 5135 cases and 25 675 controls for stage 1, and 5414 cases and 21 471 controls for stage 2. We identified 24 genome-wide significant signals of association with moderate-to-severe asthma, including several signals in innate or adaptive immune-response genes. Three novel signals were identified: rs10905284 in GATA3 (coded allele A, odds ratio [OR] 0·90, 95% CI 0·88-0·93; p=1·76 × 10-10), rs11603634 in the MUC5AC region (coded allele G, OR 1·09, 1·06-1·12; p=2·32 × 10-8), and rs560026225 near KIAA1109 (coded allele GATT, OR 1·12, 1·08-1·16; p=3·06 × 10-9). The MUC5AC signal was not associated with asthma when analyses included mild asthma. The rs11603634 G allele was associated with increased expression of MUC5AC mRNA in bronchial epithelial brush samples via proxy SNP rs11602802; (p=2·50 × 10-5) and MUC5AC mRNA was increased in bronchial epithelial samples from patients with severe asthma (in two independent analyses, p=0·039 and p=0·022). INTERPRETATION: We found substantial shared genetic architecture between mild and moderate-to-severe asthma. We also report for the first time genetic variants associated with the risk of developing moderate-to-severe asthma that regulate mucin production. Finally, we identify candidate causal genes in these loci and provide increased insight into this difficult to treat population. FUNDING: Asthma UK, AirPROM, U-BIOPRED, UK Medical Research Council, and Rosetrees Trust.

Changes in DNA Methylation from Age 18 to Pregnancy in Type 1, 2, and 17 T Helper and Regulatory T-Cells Pathway Genes.

To succeed, pregnancies need to initiate immune biases towards T helper 2 (Th2) responses, yet little is known about what establishes this bias. Using the Illumina 450 K platform, we explored changes in DNA methylation (DNAm) of Th1, Th2, Th17, and regulatory T cell pathway genes before and during pregnancy. Female participants were recruited at birth (1989), and followed through age 18 years and their pregnancy (2011-2015). Peripheral blood DNAm was measured in 245 girls at 18 years; from among these girls, the DNAm of 54 women was repeatedly measured in the first (weeks 8-21, n = 39) and second (weeks 22-38, n = 35) halves of pregnancy, respectively. M-values (logit-transformed ß-values of DNAm) were analyzed: First, with repeated measurement models, cytosine-phosphate-guanine sites (CpGs) of pathway genes in pregnancy and at age 18 (nonpregnant) were compared for changes (p ≤ 0.05). Second, we tested how many of the 348 pathway-related CpGs changed compared to 10 randomly selected subsets of all other CpGs and compared to 10 randomly selected subsets of other CD4+-related CpGs (348 in each subset). Contrasted to the nonpregnant state, 27.7% of Th1-related CpGs changed in the first and 36.1% in the second half of pregnancy. Among the Th2 pathway CpGs, proportions of changes were 35.1% (first) and 33.8% (second half). The methylation changes suggest involvement of both Th1 and Th2 pathway CpGs in the immune bias during pregnancy. Changes in regulatory T cell and Th17 pathways need further exploration.

Subclonal Evolution of Cancer-Related Gene Mutations in p53 Immunopositive Patches in Human Skin.

Normal sun-exposed skin contains numerous epidermal patches that stain positive for p53 protein (p53 immunopositive patches, PIPs), which are considered potential early precursors of skin cancer. Although the TP53 gene is mutated in many PIPs, it is unclear whether PIPs contain any other cancer-related mutations. Here we report that PIPs, predominantly <3,000 p53 immunopositive cells in size, within normal chronically exposed skin contain mutations in multiple genes that are mutated in cutaneous squamous cell cancers. These mutations in the PIPs were not detected within the non-PIP epidermis of corresponding normal chronically exposed skin. Although some of these genetic alterations are clonal in the PIPs, many of the mutations are subclonal within these lesions. Similar mutations are seen in later precancers (actinic keratoses and Bowen's disease). Our results demonstrate that PIPs in chronically exposed skin contain multiple mutations in cancer-related genes. In addition, the results indicate that the clonal evolution of mutations that are seen within later precancerous lesions and in established malignancy can also occur in PIPs within normal human skin.

Role of DNA Methylation in Type 2 Diabetes Etiology: Using Genotype as a Causal Anchor.

Several studies have investigated the relationship between genetic variation and DNA methylation with respect to type 2 diabetes, but it is unknown if DNA methylation is a mediator in the disease pathway or if it is altered in response to disease state. This study uses genotypic information as a causal anchor to help decipher the likely role of DNA methylation measured in peripheral blood in the etiology of type 2 diabetes. Illumina HumanMethylation450 BeadChip data were generated on 1,018 young individuals from the Avon Longitudinal Study of Parents and Children (ALSPAC) cohort. In stage 1, 118 unique associations between published type 2 diabetes single nucleotide polymorphisms (SNPs) and genome-wide methylation (methylation quantitative trait loci [mQTLs]) were identified. In stage 2, a further 226 mQTLs were identified between 202 additional independent non-type 2 diabetes SNPs and CpGs identified in stage 1. Where possible, associations were replicated in independent cohorts of similar age. We discovered that around half of known type 2 diabetes SNPs are associated with variation in DNA methylation and postulated that methylation could either be on a causal pathway to future disease or could be a noncausal biomarker. For one locus (KCNQ1), we were able to provide further evidence that methylation is likely to be on the causal pathway to disease in later life.

Tetanus vaccination is associated with differential DNA-methylation: Reduces the risk of asthma in adolescence.

BACKGROUND: Vaccinations have been suggested to be associated with increased risk of allergic diseases. Tetanus vaccination is one of the most frequently administered vaccines as a part of wound management and was also found to be associated with increased serum IgE levels. We hypothesized that the vaccination modifies the risk of allergic diseases through epigenetic changes such as DNA methylation. METHOD: Data on tetanus vaccination between 10 and 18years of age was collected from a birth cohort established on the Isle of Wight UK in 1989. DNA methylation data were collected from individuals at different ages (at birth [n=30], age 10 [n=34], age 18 [n=245] and during pregnancy [n=121]) using the Illumina Infinium HumanMethylation450K array. Firstly, we performed an epigenome-wide screening to identify cytosine-phosphate-guanine sites (CpGs) associated with tetanus vaccination in 18-year-olds. Secondly, we tested their association with asthma, allergic sensitization, eczema, serum IgE and pulmonary lung function (FVC, FEV1, FEV1/FVC, and FEF25-75%). We then described changes in the methylation of the selected CpG sites over age, and by vaccination status. RESULTS: Tetanus vaccination was found to be associated with decreased methylation of cg14472551 (p value 0.5×10-5, FDR-adjusted p value 2.1×10-4) and increased methylation of cg01669161 (p value 0.0007, FDR-adjusted p value 0.014). Both CpGs, in turn, were associated with decreased risk of asthma at 18years of age. Cg14472551 is located in an intron of KIAA1549L, whose protein binds to a B-cell commitment transcription factor; cg01669161 is located between an antisense regulator of the proteasome assembly chaperone PSMG3, and TFAMP1, a pseudogene. Increased methylation of cg01669161 was also associated with decreased serum IgE levels. CONCLUSION: DNA methylation changes following tetanus vaccination may offer a novel prospect to explain a differential occurrence of asthma in adolescence.

DNA methylation and genetic polymorphisms of the Leptin gene interact to influence lung function outcomes and asthma at 18 years of age.

The leptin gene (LEP) plays a regulatory role in satiety, inflammation, and allergy. Prior findings linking leptin to asthma motivated us to investigate whether DNA methylation (DNA-M) of CpG (cytosine-phosphate-guanine) sites in concert with single nucleotide polymorphisms (SNPs) of LEP can explain the risk of asthma and lung function. Methylation of CpG sites was assessed using the Illumina Infinium Human Methylation 450 beadchip in blood samples collected from 10- and 18-year-old boys and girls from the Isle of Wight (IOW) birth cohort (UK). Four LEP SNPs were genotyped. Linear and log linear models were used for the analysis, adjusting for false discovery rate (FDR). The analyses were repeated in the BAMSE cohort (Sweden). In the IOW study, the interaction of cg00666422 and rs11763517 (CT vs TT and CC) was associated with FEV1 (FDR-adjusted p-value: 0.03), FEV1/FVC ratio (FDR-adjusted p-value: 0.0096), and FEF25-75% (FDR-adjusted p-value: 0.00048) such that they decreased with increasing DNA-M. The interaction of the same CpG-SNP pair was also associated with increased risk of asthma at age 18. We replicated the findings for FEV1/FVC and FEF25-75% in a smaller sample of 34 participants at age 10. Regarding the BAMSE cohort, although, the interaction of cg00666422 and rs11763517 on lung function were not significant, the direction of the effect was the same as in IOW cohort. Thus, penetrance of LEP genotype seems to be modified by methylation at cg00666422 and is linked to airway obstruction and asthma.

An Efficient Approach to Screening Epigenome-Wide Data.

Screening cytosine-phosphate-guanine dinucleotide (CpG) DNA methylation sites in association with some covariate(s) is desired due to high dimensionality. We incorporate surrogate variable analyses (SVAs) into (ordinary or robust) linear regressions and utilize training and testing samples for nested validation to screen CpG sites. SVA is to account for variations in the methylation not explained by the specified covariate(s) and adjust for confounding effects. To make it easier to users, this screening method is built into a user-friendly R package, ttScreening, with efficient algorithms implemented. Various simulations were implemented to examine the robustness and sensitivity of the method compared to the classical approaches controlling for multiple testing: the false discovery rates-based (FDR-based) and the Bonferroni-based methods. The proposed approach in general performs better and has the potential to control both types I and II errors. We applied ttScreening to 383,998 CpG sites in association with maternal smoking, one of the leading factors for cancer risk.

Gene expression differences in relation to age and social environment in queen and worker bumble bees.

Eusocial insects provide special insights into the genetic pathways influencing aging because of their long-lived queens and flexible aging schedules. Using qRT-PCR in the primitively eusocial bumble bee Bombus terrestris (Linnaeus), we investigated expression levels of four candidate genes associated with taxonomically widespread age-related pathways (coenzyme Q biosynthesis protein 7, COQ7; DNA methyltransferase 3, Dnmt3; foraging, for; and vitellogenin, vg). In Experiment 1, we tested how expression changes with queen relative age and productivity. We found a significant age-related increase in COQ7 expression in queen ovary. In brain, all four genes showed higher expression with increasing female (queen plus worker) production, with this relationship strengthening as queen age increased, suggesting a link with the positive association of fecundity and longevity found in eusocial insect queens. In Experiment 2, we tested effects of relative age and social environment (worker removal) in foundress queens and effects of age and reproductive status in workers. In this experiment, workerless queens showed significantly higher for expression in brain, as predicted if downregulation of for is associated with the cessation of foraging by foundress queens following worker emergence. Workers showed a significant age-related increase in Dnmt3 expression in fat body, suggesting a novel association between aging and methylation in B. terrestris. Ovary activation was associated with significantly higher vg expression in fat body and, in younger workers, in brain, consistent with vitellogenin's ancestral role in regulating egg production. Overall, our findings reveal a mixture of novel and conserved features in age-related genetic pathways under primitive eusociality.

TSLP polymorphisms, allergen exposures, and the risk of atopic disorders in children.

BACKGROUND: Thymic stromal lymphopoietin (TSLP) polymorphisms influence atopy risk. TSLP might constitute a key interface between the environment and the allergic immune response. However, whether the effects of TSLP polymorphisms on atopic dermatitis (AD) are modified by allergic sensitization is not clear. OBJECTIVE: To evaluate the joint effect of allergic sensitization and TSLP polymorphisms on AD and to test whether TSLP polymorphisms increase the risk of asthma in children with AD. METHODS: A total of 1,520 kindergarten children (375 with AD and 1,145 controls) selected from the Childhood Environment and Allergic Diseases Study cohort in 2010 were enrolled. Information about allergic diseases and environmental exposures was collected by questionnaire. Skin prick tests were performed to measure allergic sensitization. TSLP polymorphisms were genotyped by TaqMan assay. Logistic regressions were conducted to estimate the association among TSLP polymorphisms, allergic sensitization, and AD. For replication, a subsample of the British Isle of Wight birth cohort was used. RESULTS: The TSLP rs2289278 CC genotype increased the risk of AD (odds ratio 1.90, 95% confidence interval 1.12-3.22). In children sensitized to certain allergens, a genetic predisposition (rs2289278 genotype CC) significantly increased the risk of AD. These findings were replicated in the British subsample using rs2289276 genotypes TT and TC, which are in linkage disequilibrium with rs2289278. In subjects with AD, the rs2289278 C allele also significantly increased the risk of developing asthma (odds ratio 8.31, 95% confidence interval 1.08-64.13). CONCLUSION: The association of rs2289278 with AD was stronger in children with allergic sensitization than in children without atopy. TSLP polymorphisms also increased the risk of asthma in children with AD.

Multi-ancestry genome-wide association study of 21,000 cases and 95,000 controls identifies new risk loci for atopic dermatitis.

Genetic association studies have identified 21 loci associated with atopic dermatitis risk predominantly in populations of European ancestry. To identify further susceptibility loci for this common, complex skin disease, we performed a meta-analysis of >15 million genetic variants in 21,399 cases and 95,464 controls from populations of European, African, Japanese and Latino ancestry, followed by replication in 32,059 cases and 228,628 controls from 18 studies. We identified ten new risk loci, bringing the total number of known atopic dermatitis risk loci to 31 (with new secondary signals at four of these loci). Notably, the new loci include candidate genes with roles in the regulation of innate host defenses and T cell function, underscoring the important contribution of (auto)immune mechanisms to atopic dermatitis pathogenesis.

DNA methylation loci associated with atopy and high serum IgE: a genome-wide application of recursive Random Forest feature selection.

BACKGROUND: The prevalence of allergic diseases are increasing worldwide, emphasizing the need to elucidate their pathogeneses. The aims of this study were to use a two-stage design to identify DNA methylation levels at cytosine-phosphate-guanine (CpG) sites across the genome associated with atopy and high serum immunoglobulin E (IgE), then to replicate our findings in an independent cohort. METHODS: Atopy was assessed via skin prick tests and high serum IgE. Methylation levels were measured from whole blood using the Illumina Infinium HumanMethylation450 BeadChip from 18-year-old women (n = 245) and men (n = 122) in the Isle of Wight birth cohort. After data cleaning and processing, and removing probes with possible single nucleotide polymorphisms, DNA methylation levels from 254,460 CpG sites from the 245 women were subjected to recursive Random Forest feature selection for stage 1. The sites selected from stage 1 were tested in stage 2 for associations with atopy and high IgE levels (>200 kU/L) via logistic regression adjusted for predicted cell-type proportions and sex. Sites significantly associated with atopy in stage 2 underwent replication tests in the independent Swedish birth cohort BAMSE (n = 464). RESULTS: In stage 1, 62 sites were selected, of which 22 were associated with atopy in stage 2 (P-value range 6.5E-9 to 1.4E-5) and 12 associated with high IgE levels (P-value range 1.1E-5 to 7.1E-4) at the Bonferroni adjusted alpha (0.05/62 = 0.0008). Of the 19 available sites, 13 were replicated. CONCLUSIONS: We identified 13 novel epigenetic loci associated with atopy and high IgE that could serve as candidate loci for future studies; four were within genes with known roles in the immune response (cg04983687 in the body of ZFPM1, cg18219873 in the 5'UTR of PRG2, cg27469152 in the 3'UTR of EPX, and cg09332506 in the body of COPA).

Does allergy begin in utero?

It has been recognized for centuries that allergic disease runs in families, implying a role for genetic factors in determining individual susceptibility. More recently, a range of evidence shows that many of these genetic factors, together with in utero environmental exposures, lead to the development of allergic disease through altered immune and organ development. Environmental exposures during pregnancy including diet, nutrient intake and toxin exposures can alter the epigenome and interact with inherited genetic and epigenetic risk factors to directly and indirectly influence organ development and immune programming. Understanding of these factors will be essential in identifying at-risk individuals and possible development of therapeutic interventions for the primary prevention of allergic disease. In this review, we summarize the evidence that suggests allergic disease begins in utero, together with possible mechanisms for the effect of environmental exposures during pregnancy on allergic disease risk, including epigenetics.

A statistical method for single sample analysis of HumanMethylation450 array data: genome-wide methylation analysis of patients with imprinting disorders.

BACKGROUND: The Illumina Infinium HumanMethylation450 BeadChip is an array-based technology for analysing DNA methylation at approximately 475,000 differentially methylated cytosines across the human genome. Hitherto, the array has been used for case-control studies, where sample numbers can be sufficient to yield statistically robust data on a genome-wide basis. We recently reported an informatic pipeline capable of yielding statistically and biologically significant results using only five cases, which expanded the use of this technology to rare disease studies. However, the clinical application of these technologies requires the ability to perform robust analysis of individual patients. RESULTS: Here we report a novel informatic approach for methylation array analysis of single samples, using the Crawford-Howell t-test. We tested our approach on patients with ultra-rare imprinting disorders with aberrant DNA methylation at multiple locations across the genome, which was previously detected by targeted testing. However, array analysis outperformed targeted assays in three ways: it detected loci not normally analysed by targeted testing, detected methylation changes too subtle to detect by the targeted testing and reported broad and consistent methylation changes across genetic loci not captured by point testing. CONCLUSIONS: This method has potential clinical utility for human disorders where DNA methylation change may be a biomarker of disease.

Oral contraceptives modify the effect of GATA3 polymorphisms on the risk of asthma at the age of 18 years via DNA methylation.

BACKGROUND: The prevalence of asthma in girls increases after puberty. Previous studies have detected associations between sex hormones and asthma, as well as between sex hormones and T helper 2 (Th2) asthma-typical immune responses. Therefore, we hypothesized that exogenous or endogenous sex hormone exposure (represented by oral contraceptive pill (OCP) use and early menarche, respectively) are associated with DNA methylation (DNA-M) of the Th2 transcription factor gene, GATA3, in turn affecting the risk of asthma in girls, possibly in interaction with genetic variants. Blood samples were collected from 245 female participants aged 18 years randomly selected for methylation analysis from the Isle of Wight birth cohort, UK. Information on use of OCPs, age at menarche, and concurrent asthma were assessed by questionnaire. Genome-wide DNA-M was determined using the Illumina Infinium HumanMethylation450 beadchip. In a first stage, we tested the interaction between sex hormone exposure and genetic variants on DNA-M of specific cytosine-phosphate-guanine (CpG) sites. In a second stage, we determined whether these CpG sites interact with genetic variants in GATA3 to explain the risk of asthma. RESULTS: Interactions between OCP use and seven single nucleotide polymorphisms (SNPs) of GATA3 were analyzed for 14 CpG sites (stage 1). The interaction between OCP use and SNP rs1269486 was found to be associated with the methylation level of cg17124583 (P = 0.002, false discovery rate (FDR) adjusted P = 0.04). DNA-M of this same CpG site was also influenced by the interaction between age at menarche and rs1269486 (P = 0.0017). In stage 2, we found that cg17124583 modified the association of SNP rs422628 with asthma risk at the age of 18 years (P = 0.006, FDR adjusted P = 0.04). Subjects with genotype AG showed an increase in average risk ratio (RR) from 0.31 (95% CI: 0.10 to 0.8) to 11.65 (95% CI: 1.71 to 79.5) when methylation level increased from 0.02 to 0.12, relative to genotype AA. CONCLUSION: A two-stage model consisting of genetic variants in the GATA3 gene, OCP use, age at menarche, and DNA-M may explain how sex hormones in women can increase the asthma prevalence after puberty.

The interplay of DNA methylation over time with Th2 pathway genetic variants on asthma risk and temporal asthma transition.

BACKGROUND: Genetic effects on asthma of genes in the T-helper 2 (Th2) pathway may interact with epigenetic factors including DNA methylation. We hypothesized that interactions between genetic variants and methylation in genes in this pathway (IL4, IL4R, IL13, GATA3, and STAT6) influence asthma risk, that such influences are age-dependent, and that methylation of some CpG sites changes over time in accordance with asthma transition. We tested these hypotheses in subsamples of girls from a population-based birth cohort established on the Isle of Wight, UK, in 1989. RESULTS: Logistic regression models were applied to test the interaction effect of DNA methylation and SNP on asthma within each of the five genes. Bootstrapping was used to assess the models identified. From 1,361 models fitted at each age of 10 and 18 years, 8 models, including 4 CpGs and 8 SNPs, showed potential associations with asthma risk. Of the 4 CpGs, methylation of cg26937798 (IL4R) and cg23943829 (IL4) changes between ages 10 and 18 (both higher at 10; P = 9.14 × 10(-6) and 1.07 × 10(-5), respectively). At age 10, the odds of asthma tended to decrease as cg12405139 (GATA3) methylation increased (log-OR = -12.15; P = 0.049); this effect disappeared by age 18. At age 18, methylation of cg09791102 (IL4R) was associated with higher risk of asthma among subjects with genotype GG compared to AG (P = 0.003), increased cg26937798 methylation among subjects with rs3024685 (IL4R) genotype AA (P = 0.003) or rs8832 (IL4R) genotype GG (P = 0.01) was associated with a lower asthma risk; these CpGs had no effect at age 10. Increasing cg26937798 methylation over time possibly reduced the risk of positive asthma transition (asthma-free at age 10 → asthma at age 18; log-OR = -3.11; P = 0.069) and increased the likelihood of negative transition (asthma at age 10 → asthma-free at age 18; log-OR = 3.97; P = 0.074). CONCLUSIONS: The interaction of DNA methylation and SNPs in Th2 pathway genes is likely to contribute to asthma risk. This effect may vary with age. Methylation of some CpGs changed over time, which may influence asthma transition.

Genome-wide DNA methylation analysis of patients with imprinting disorders identifies differentially methylated regions associated with novel candidate imprinted genes.

BACKGROUND: Genomic imprinting is allelic restriction of gene expression potential depending on parent of origin, maintained by epigenetic mechanisms including parent of origin-specific DNA methylation. Among approximately 70 known imprinted genes are some causing disorders affecting growth, metabolism and cancer predisposition. Some imprinting disorder patients have hypomethylation of several imprinted loci (HIL) throughout the genome and may have atypically severe clinical features. Here we used array analysis in HIL patients to define patterns of aberrant methylation throughout the genome. DESIGN: We developed a novel informatic pipeline capable of small sample number analysis, and profiled 10 HIL patients with two clinical presentations (Beckwith-Wiedemann syndrome and neonatal diabetes) using the Illumina Infinium Human Methylation450 BeadChip array to identify candidate imprinted regions. We used robust statistical criteria to quantify DNA methylation. RESULTS: We detected hypomethylation at known imprinted loci, and 25 further candidate imprinted regions (nine shared between patient groups) including one in the Down syndrome critical region (WRB) and another previously associated with bipolar disorder (PPIEL). Targeted analysis of three candidate regions (NHP2L1, WRB and PPIEL) showed allelic expression, methylation patterns consistent with allelic maternal methylation and frequent hypomethylation among an additional cohort of HIL patients, including six with Silver-Russell syndrome presentations and one with pseudohypoparathyroidism 1B. CONCLUSIONS: This study identified novel candidate imprinted genes, revealed remarkable epigenetic convergence among clinically divergent patients, and highlights the potential of epigenomic profiling to expand our understanding of the normal methylome and its disruption in human disease.