Systematic identification of genotype-dependent enhancer variants in eosinophilic esophagitis.

Eosinophilic esophagitis (EoE) is a rare atopic disorder associated with esophageal dysfunction, including difficulty swallowing, food impaction, and inflammation, that develops in a small subset of people with food allergies. Genome-wide association studies (GWASs) have identified 9 independent EoE risk loci reaching genome-wide significance (p < 5 × 10-8) and 27 additional loci of suggestive significance (5 × 10-8 < p < 1 × 10-5). In the current study, we perform linkage disequilibrium (LD) expansion of these loci to nominate a set of 531 variants that are potentially causal. To systematically interrogate the gene regulatory activity of these variants, we designed a massively parallel reporter assay (MPRA) containing the alleles of each variant within their genomic sequence context cloned into a GFP reporter library. Analysis of reporter gene expression in TE-7, HaCaT, and Jurkat cells revealed cell-type-specific gene regulation. We identify 32 allelic enhancer variants, representing 6 genome-wide significant EoE loci and 7 suggestive EoE loci, that regulate reporter gene expression in a genotype-dependent manner in at least one cellular context. By annotating these variants with expression quantitative trait loci (eQTL) and chromatin looping data in related tissues and cell types, we identify putative target genes affected by genetic variation in individuals with EoE. Transcription factor enrichment analyses reveal possible roles for cell-type-specific regulators, including GATA3. Our approach reduces the large set of EoE-associated variants to a set of 32 with allelic regulatory activity, providing functional insights into the effects of genetic variation in this disease.

Very early onset eosinophilic esophagitis is common, responds to standard therapy, and demonstrates enrichment for CAPN14 genetic variants.

BACKGROUND: Eosinophilic esophagitis (EoE) is a chronic, food antigen-mediated disease characterized by esophageal dysfunction and intraepithelial eosinophil accumulation. OBJECTIVE: We hypothesized that very early onset EoE (V-EoE) would be enriched for early-life and genetic factors and have worse presentation and prognosis than later-onset pediatric EoE (L-EoE). METHODS: We conducted a single-site, retrospective review comparing patients diagnosed at age 12 months or less (V-EoE, n = 57) and age 14 to 18 years (L-EoE, n = 70). These patients underwent medical record, EoE Histology Scoring System, Endoscopic Reference Score, and EoE Diagnostic Panel assessment when sample availability permitted. Genetic association used 2 EoE genotype repositories. Data were analyzed using chi-square tests, t tests, Wilcoxon rank-sum tests, Spearman correlations, cluster analysis, and logistic regression. RESULTS: Among pediatric patients with EoE, diagnosis most commonly occurred within early life (0-24 months, 17%). V-EoE was more likely to attain histologic remission via dietary restriction (P < .0001). Basal zone hyperplasia and eosinophil inflammation were greater in V-EoE (P < .05). Esophageal strictures more commonly occurred in L-EoE (P = .03). V-EoE had lower endoscopic scores (P < .05). Molecular expression was very similar between groups. Cesarean delivery was more common in patients with V-EoE (P = .03). Patients with V-EoE demonstrated enrichment of CAPN14 common genetic variants. CONCLUSIONS: Early-life diagnosis of EoE is a common occurrence. V-EoE responds to standard therapy without early evidence for complications, suggesting a less severe prognosis than hypothesized. Molecular pathogenesis is preserved between V-EoE and L-EoE. Cesarean delivery and CAPN14 genetic variation likely promote earlier disease development.

Replication and meta-analyses nominate numerous eosinophilic esophagitis risk genes.

BACKGROUND: Eosinophilic esophagitis (EoE) is an emerging, chronic, rare allergic disease associated with marked eosinophil accumulation in the esophagus. Previous genome-wide association studies have provided strong evidence for 3 genome-wide susceptibility loci. OBJECTIVE: We sought to replicate known and suggestive EoE genetic risk loci and conduct a meta-analysis of previously reported data sets. METHODS: An EoE-Custom single-nucleotide polymophism (SNP) Chip containing 956 candidate EoE risk single-nucleotide polymorphisms was used to genotype 627 cases and 365 controls. Statistical power was enhanced by adding 1959 external controls and performing meta-analyses with 2 independent EoE genome-wide association studies. RESULTS: Meta-analysis identified replicated association and genome-wide significance at 6 loci: 2p23 (2 independent genetic effects) and 5q22, 10p14, 11q13, and 16p13. Seven additional loci were identified at suggestive significance (P < 10-6): 1q31, 5q23, 6q15, 6q21, 8p21, 17q12, and 22q13. From these risk loci, 13 protein-coding EoE candidate risk genes were expressed in a genotype-dependent manner. EoE risk genes were expressed in disease-relevant cell types, including esophageal epithelia, fibroblasts, and immune cells, with some expressed as a function of disease activity. The genetic risk burden of EoE-associated genetic variants was markedly larger in cases relative to controls (P < 10-38); individuals with the highest decile of genetic burden had greater than 12-fold risk of EoE compared with those within the lowest decile. CONCLUSIONS: This study extends the genetic underpinnings of EoE, highlighting 13 genes whose genotype-dependent expression expands our etiologic understanding of EoE and provides a framework for a polygenic risk score to be validated in future studies.

Genetic, Inflammatory, and Epithelial Cell Differentiation Factors Control Expression of Human Calpain-14.

Eosinophilic esophagitis (EoE) is a chronic, food-driven allergic disease resulting in eosinophilic esophageal inflammation. We recently found that EoE susceptibility is associated with genetic variants in the promoter of CAPN14, a gene with reported esophagus-specific expression. CAPN14 is dynamically up-regulated as a function of EoE disease activity and after exposure of epithelial cells to interleukin-13 (IL-13). Herein, we aimed to explore molecular modulation of CAPN14 expression. We identified three putative binding sites for the IL-13-activated transcription factor STAT6 in the promoter and first intron of CAPN14 Luciferase reporter assays revealed that the two most distal STAT6 elements were required for the ∼10-fold increase in promoter activity subsequent to stimulation with IL-13 or IL-4, and also for the genotype-dependent reduction in IL-13-induced promoter activity. One of the STAT6 elements in the promoter was necessary for IL-13-mediated induction of CAPN14 promoter activity while the other STAT6 promoter element was necessary for full induction. Chromatin immunoprecipitation in IL-13 stimulated esophageal epithelial cells was used to further support STAT6 binding to the promoter of CAPN14 at these STAT6 binding sites. The highest CAPN14 and calpain-14 expression occurred with IL-13 or IL-4 stimulation of esophageal epithelial cells under culture conditions that allow the cells to differentiate into a stratified epithelium. This work corroborates a candidate molecular mechanism for EoE disease etiology in which the risk variant at 2p23 dampens CAPN14 expression in differentiated esophageal epithelial cells following IL-13/STAT6 induction of CAPN14 promoter activity.

Eosinophilic oesophagitis endotype classification by molecular, clinical, and histopathological analyses: a cross-sectional study.

BACKGROUND: Eosinophilic oesophagitis is understood in terms of quantifiable histological, endoscopic, and molecular features. Data are scant for inter-relations of these features and their potential to identify distinct disease endotypes. We aimed to identify clinical-pathological correlations between endoscopic and histological disease variables by transcription profiling of the oesophagus of patients with eosinophilic oesophagitis of varying severity and disease activity states. METHODS: We did a cross-sectional study across ten hospital sites in the USA associated with the Consortium of Eosinophilic Gastrointestinal Disease Researchers. We analysed oesophageal biopsy specimens taken from paediatric and adult patients with eosinophilic oesophagitis (discovery cohort), using the eosinophilic oesophagitis diagnostic panel (EDP), a set of 96 informative transcripts. Histological and endoscopic features were assessed by quantification of oesophageal eosinophils and use of the eosinophilic oesophagitis histology scoring system (HSS) and the eosinophilic oesophagitis endoscopic reference score (EREFS). Associations among the various histological, endoscopic, and molecular features were analysed by Spearman correlation. Results were replicated in a biologically independent, single-centre, validation cohort of patients with active eosinophilic oesophagitis. FINDINGS: The discovery cohort contained 185 samples and the validation cohort comprised 100 specimens. In the discovery cohort, EDP showed intersite consistency, significant correlation with oesophageal eosinophils (p<0·0001), and similar findings between paediatric and adult patients. Of eight HSS domains, basal zone hyperplasia correlated with the EDP (median Spearman ρ 0·47 [IQR 0·36-0·60]). Of five EREFS features, distal furrows correlated with the EDP (median Spearman ρ 0·42 [0·32-0·50]). By analysing active eosinophilic oesophagitis in the discovery cohort, the EDP identified three clusters associated with distinct endotypes (termed EoEe1-3) despite similar eosinophil levels. EoEe1 was associated with a normal-appearing oesophagus (risk ratio [RR] 3·27, 95% CI 1·04-10·27; p=0·0443), an inverse association with a history of oesophageal dilation (0·27, 0·09-0·82; p=0·0105) and showed relatively mild histological, endoscopic, and molecular changes. EoEe2 showed an inflammatory and steroid-refractory phenotype (RR 2·77, 95% CI 1·11-6·95; p=0·0376) and had the highest expression of inflammatory cytokines and steroid-responding genes. EoEe3 was associated with a narrow-calibre oesophagus (RR 7·98, 95% CI 1·84-34·64; p=0·0013) and adult onset (2·22, 1·19-4·12; p=0·0155), and showed the highest degree of endoscopic and histological severity and the lowest expression of epithelial differentiation genes. These endotypes were replicated in the validation cohort by clustering and with an eosinophilic oesophagitis endotype-prediction algorithm. INTERPRETATION: Our new disease classification stratifies patients with eosinophilic oesophagitis into subgroups with potential clinical and therapeutic significance and provides a framework for a precision medicine approach to eosinophilic oesophagitis. FUNDING: National Institutes of Health.

Identification of endometrial cancer methylation features using combined methylation analysis methods.

BACKGROUND: DNA methylation is a stable epigenetic mark that is frequently altered in tumors. DNA methylation features are attractive biomarkers for disease states given the stability of DNA methylation in living cells and in biologic specimens typically available for analysis. Widespread accumulation of methylation in regulatory elements in some cancers (specifically the CpG island methylator phenotype, CIMP) can play an important role in tumorigenesis. High resolution assessment of CIMP for the entire genome, however, remains cost prohibitive and requires quantities of DNA not available for many tissue samples of interest. Genome-wide scans of methylation have been undertaken for large numbers of tumors, and higher resolution analyses for a limited number of cancer specimens. Methods for analyzing such large datasets and integrating findings from different studies continue to evolve. An approach for comparison of findings from a genome-wide assessment of the methylated component of tumor DNA and more widely applied methylation scans was developed. METHODS: Methylomes for 76 primary endometrial cancer and 12 normal endometrial samples were generated using methylated fragment capture and second generation sequencing, MethylCap-seq. Publically available Infinium HumanMethylation 450 data from The Cancer Genome Atlas (TCGA) were compared to MethylCap-seq data. RESULTS: Analysis of methylation in promoter CpG islands (CGIs) identified a subset of tumors with a methylator phenotype. We used a two-stage approach to develop a 13-region methylation signature associated with a "hypermethylator state." High level methylation for the 13-region methylation signatures was associated with mismatch repair deficiency, high mutation rate, and low somatic copy number alteration in the TCGA test set. In addition, the signature devised showed good agreement with previously described methylation clusters devised by TCGA. CONCLUSION: We identified a methylation signature for a "hypermethylator phenotype" in endometrial cancer and developed methods that may prove useful for identifying extreme methylation phenotypes in other cancers.

Methods for high-throughput MethylCap-Seq data analysis.

BACKGROUND: Advances in whole genome profiling have revolutionized the cancer research field, but at the same time have raised new bioinformatics challenges. For next generation sequencing (NGS), these include data storage, computational costs, sequence processing and alignment, delineating appropriate statistical measures, and data visualization. Currently there is a lack of workflows for efficient analysis of large, MethylCap-seq datasets containing multiple sample groups. METHODS: The NGS application MethylCap-seq involves the in vitro capture of methylated DNA and subsequent analysis of enriched fragments by massively parallel sequencing. The workflow we describe performs MethylCap-seq experimental Quality Control (QC), sequence file processing and alignment, differential methylation analysis of multiple biological groups, hierarchical clustering, assessment of genome-wide methylation patterns, and preparation of files for data visualization. RESULTS: Here, we present a scalable, flexible workflow for MethylCap-seq QC, secondary data analysis, tertiary analysis of multiple experimental groups, and data visualization. We demonstrate the experimental QC procedure with results from a large ovarian cancer study dataset and propose parameters which can identify problematic experiments. Promoter methylation profiling and hierarchical clustering analyses are demonstrated for four groups of acute myeloid leukemia (AML) patients. We propose a Global Methylation Indicator (GMI) function to assess genome-wide changes in methylation patterns between experimental groups. We also show how the workflow facilitates data visualization in a web browser with the application Anno-J. CONCLUSIONS: This workflow and its suite of features will assist biologists in conducting methylation profiling projects and facilitate meaningful biological interpretation.

Genome-wide methylation profiling in decitabine-treated patients with acute myeloid leukemia.

The outcome of older (≥ 60 years) acute myeloid leukemia (AML) patients is poor, and novel treatments are needed. In a phase 2 trial for older AML patients, low-dose (20 mg/m(2) per day for 10 days) decitabine, a DNA hypomethylating azanucleoside, produced 47% complete response rate with an excellent toxicity profile. To assess the genome-wide activity of decitabine, we profiled pretreatment and post treatment (day 25/course 1) methylomes of marrow samples from patients (n = 16) participating in the trial using deep-sequencing analysis of methylated DNA captured by methyl-binding protein (MBD2). Decitabine significantly reduced global methylation compared with pretreatment baseline (P = .001). Percent marrow blasts did not correlate with global methylation levels, suggesting that hypomethylation was related to the activity of decitabine rather than to a mere decrease in leukemia burden. Hypomethylation occurred predominantly in CpG islands and CpG island-associated regions (P ranged from .03 to .04) A significant concentration (P < .001) of the hypomehtylated CpG islands was found in chromosome subtelomeric regions, suggesting a differential activity of decitabine in distinct chromosome regions. Hypermethylation occurred much less frequently than hypomethylation and was associated with low CpG content regions. Decitabine-related methylation changes were concordant with those previously reported in distinct genes. In summary, our study supports the feasibility of methylome analyses as a pharmacodynamic endpoint for hypomethylating therapies.

Enrichment-based DNA methylation analysis using next-generation sequencing: sample exclusion, estimating changes in global methylation, and the contribution of replicate lanes.

BACKGROUND: DNA methylation is an important epigenetic mark and dysregulation of DNA methylation is associated with many diseases including cancer. Advances in next-generation sequencing now allow unbiased methylome profiling of entire patient cohorts, greatly facilitating biomarker discovery and presenting new opportunities to understand the biological mechanisms by which changes in methylation contribute to disease. Enrichment-based sequencing assays such as MethylCap-seq are a cost effective solution for genome-wide determination of methylation status, but the technical reliability of methylation reconstruction from raw sequencing data has not been well characterized. METHODS: We analyze three MethylCap-seq data sets and perform two different analyses to assess data quality. First, we investigate how data quality is affected by excluding samples that do not meet quality control cutoff requirements. Second, we consider the effect of additional reads on enrichment score, saturation, and coverage. Lastly, we verify a method for the determination of the global amount of methylation from MethylCap-seq data by comparing to a spiked-in control DNA of known methylation status. RESULTS: We show that rejection of samples based on our quality control parameters leads to a significant improvement of methylation calling. Additional reads beyond ~13 million unique aligned reads improved coverage, modestly improved saturation, and did not impact enrichment score. Lastly, we find that a global methylation indicator calculated from MethylCap-seq data correlates well with the global methylation level of a sample as obtained from a spike-in DNA of known methylation level. CONCLUSIONS: We show that with appropriate quality control MethylCap-seq is a reliable tool, suitable for cohorts of hundreds of patients, that provides reproducible methylation information on a feature by feature basis as well as information about the global level of methylation.

Cancer epigenetics: a perspective on the role of DNA methylation in acquired endocrine resistance.

Epigenetic mechanisms, including DNA methylation, are responsible for determining and maintaining cell fate, stably differentiating the various tissues in our bodies. Increasing evidence shows that DNA methylation plays a significant role in cancer, from the silencing of tumor suppressors to the activation of oncogenes and the promotion of metastasis. Recent studies also suggest a role for DNA methylation in drug resistance. This perspective article discusses how DNA methylation may contribute to the development of acquired endocrine resistance, with a focus on breast cancer. In addition, we discuss DNA methylome profiling and how recent developments in this field are shedding new light on the role of epigenetics in endocrine resistance. Hormone ablation is the therapy of choice for hormone-sensitive breast tumors, yet as many as 40% of patients inevitably relapse, and these hormone refractory tumors often have a poor prognosis. Epigenetic studies could provide DNA methylation biomarkers to predict and diagnose acquired resistance in response to treatment. Elucidation of epigenetic mechanisms may also lead to the development of new treatments that specifically target epigenetic abnormalities or vulnerabilities in cancer cells. Expectations must be tempered by the fact that epigenetic mechanisms of endocrine resistance remain poorly understood, and further study is required to better understand how altering epigenetic pathways with therapeutics can promote or inhibit endocrine resistance in different contexts. Going forward, DNA methylome profiling will become increasingly central to epigenetic research, heralding a network-based approach to epigenetics that promises to advance our understanding of the etiology of cancer in ways not previously possible.

Contactin 4 as an autism susceptibility locus.

Structural and sequence variation have been described in several members of the contactin (CNTN) and contactin-associated protein (CNTNAP) gene families in association with neurodevelopmental disorders, including autism. Using array comparative genome hybridization (CGH), we identified a maternally inherited ∼535 kb deletion at 3p26.3 encompassing the 5' end of the contactin 4 gene (CNTN4) in a patient with autism. Based on this finding and previous reports implicating genomic rearrangements of CNTN4 in autism spectrum disorders (ASDs) and 3p- microdeletion syndrome, we undertook sequencing of the coding regions of the gene in a local ASD cohort in comparison with a set of controls. Unique missense variants were identified in 4 of 75 unrelated individuals with ASD, as well as in 1 of 107 controls. All of the amino acid substitutions were nonsynonomous, occurred at evolutionarily conserved positions, and were, thus, felt likely to be deleterious. However, these data did not reach statistical significance, nor did the variants segregate with disease within all of the ASD families. Finally, there was no detectable difference in binding of two of the variants to the interacting protein PTPRG in vitro. Thus, additional larger studies will be necessary to determine whether CNTN4 functions as an autism susceptibility locus in combination with other genetic and/or environmental factors.

A Scalable, Flexible Workflow for MethylCap-Seq Data Analysis.

Advances in whole genome profiling have revolutionized the cancer research field, but at the same time have raised new bioinformatics challenges. For next generation sequencing (NGS), these include data storage, computational costs, sequence processing and alignment, delineating appropriate statistical measures, and data visualization. The NGS application MethylCap-seq involves the in vitro capture of methylated DNA and subsequent analysis of enriched fragments by massively parallel sequencing. Here, we present a scalable, flexible workflow for MethylCap-seq Quality Control, secondary data analysis, tertiary analysis of multiple experimental groups, and data visualization. This workflow and its suite of features will assist biologists in conducting methylation profiling projects and facilitate meaningful biological interpretation.