51
|
Cancer Biomarkers from Genome-Scale DNA Methylation: Comparison of Evolutionary and Semantic Analysis Methods. MICROARRAYS 2015; 4:647-70. [PMID: 27600245 PMCID: PMC4996413 DOI: 10.3390/microarrays4040647] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/27/2015] [Revised: 11/09/2015] [Accepted: 11/18/2015] [Indexed: 11/16/2022]
Abstract
DNA methylation profiling exploits microarray technologies, thus yielding a wealth of high-volume data. Here, an intelligent framework is applied, encompassing epidemiological genome-scale DNA methylation data produced from the Illumina’s Infinium Human Methylation 450K Bead Chip platform, in an effort to correlate interesting methylation patterns with cancer predisposition and, in particular, breast cancer and B-cell lymphoma. Feature selection and classification are employed in order to select, from an initial set of ~480,000 methylation measurements at CpG sites, predictive cancer epigenetic biomarkers and assess their classification power for discriminating healthy versus cancer related classes. Feature selection exploits evolutionary algorithms or a graph-theoretic methodology which makes use of the semantics information included in the Gene Ontology (GO) tree. The selected features, corresponding to methylation of CpG sites, attained moderate-to-high classification accuracies when imported to a series of classifiers evaluated by resampling or blindfold validation. The semantics-driven selection revealed sets of CpG sites performing similarly with evolutionary selection in the classification tasks. However, gene enrichment and pathway analysis showed that it additionally provides more descriptive sets of GO terms and KEGG pathways regarding the cancer phenotypes studied here. Results support the expediency of this methodology regarding its application in epidemiological studies.
Collapse
|
52
|
Sun Z, Cunningham J, Slager S, Kocher JP. Base resolution methylome profiling: considerations in platform selection, data preprocessing and analysis. Epigenomics 2015; 7:813-28. [PMID: 26366945 PMCID: PMC4790440 DOI: 10.2217/epi.15.21] [Citation(s) in RCA: 71] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
Bisulfite treatment-based methylation microarray (mainly Illumina 450K Infinium array) and next-generation sequencing (reduced representation bisulfite sequencing, Agilent SureSelect Human Methyl-Seq, NimbleGen SeqCap Epi CpGiant or whole-genome bisulfite sequencing) are commonly used for base resolution DNA methylome research. Although multiple tools and methods have been developed and used for the data preprocessing and analysis, confusions remains for these platforms including how and whether the 450k array should be normalized; which platform should be used to better fit researchers' needs; and which statistical models would be more appropriate for differential methylation analysis. This review presents the commonly used platforms and compares the pros and cons of each in methylome profiling. We then discuss approaches to study design, data normalization, bias correction and model selection for differentially methylated individual CpGs and regions.
Collapse
Affiliation(s)
- Zhifu Sun
- Division of Biomedical Statistics & Informatics, Mayo Clinic, Rochester, MN 55905, USA
| | | | - Susan Slager
- Division of Biomedical Statistics & Informatics, Mayo Clinic, Rochester, MN 55905, USA
| | - Jean-Pierre Kocher
- Division of Biomedical Statistics & Informatics, Mayo Clinic, Rochester, MN 55905, USA
| |
Collapse
|
53
|
Kling T, Johansson P, Sanchez J, Marinescu VD, Jörnsten R, Nelander S. Efficient exploration of pan-cancer networks by generalized covariance selection and interactive web content. Nucleic Acids Res 2015; 43:e98. [PMID: 25953855 PMCID: PMC4551906 DOI: 10.1093/nar/gkv413] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2014] [Accepted: 04/17/2015] [Indexed: 12/25/2022] Open
Abstract
Statistical network modeling techniques are increasingly important tools to analyze cancer genomics data. However, current tools and resources are not designed to work across multiple diagnoses and technical platforms, thus limiting their applicability to comprehensive pan-cancer datasets such as The Cancer Genome Atlas (TCGA). To address this, we describe a new data driven modeling method, based on generalized Sparse Inverse Covariance Selection (SICS). The method integrates genetic, epigenetic and transcriptional data from multiple cancers, to define links that are present in multiple cancers, a subset of cancers, or a single cancer. It is shown to be statistically robust and effective at detecting direct pathway links in data from TCGA. To facilitate interpretation of the results, we introduce a publicly accessible tool (cancerlandscapes.org), in which the derived networks are explored as interactive web content, linked to several pathway and pharmacological databases. To evaluate the performance of the method, we constructed a model for eight TCGA cancers, using data from 3900 patients. The model rediscovered known mechanisms and contained interesting predictions. Possible applications include prediction of regulatory relationships, comparison of network modules across multiple forms of cancer and identification of drug targets.
Collapse
Affiliation(s)
- Teresia Kling
- Sahlgrenska Cancer Center and Dept of Molecular and Clinical Medicine, University of Gothenburg, Box 425, SE-405 30 Gothenburg, Sweden
| | - Patrik Johansson
- Department of Immunology, Genetics and Pathology (IGP) and Science for Life Laboratory, Uppsala University, Rudbecklaboratoriet, SE-751 85 Uppsala, Sweden
| | - José Sanchez
- Mathematical Sciences, University of Gothenburg and Chalmers University of Technology, SE-412 96 Gothenburg, Sweden
| | - Voichita D Marinescu
- Department of Immunology, Genetics and Pathology (IGP) and Science for Life Laboratory, Uppsala University, Rudbecklaboratoriet, SE-751 85 Uppsala, Sweden
| | - Rebecka Jörnsten
- Mathematical Sciences, University of Gothenburg and Chalmers University of Technology, SE-412 96 Gothenburg, Sweden
| | - Sven Nelander
- Department of Immunology, Genetics and Pathology (IGP) and Science for Life Laboratory, Uppsala University, Rudbecklaboratoriet, SE-751 85 Uppsala, Sweden
| |
Collapse
|
54
|
Demerath EW, Guan W, Grove ML, Aslibekyan S, Mendelson M, Zhou YH, Hedman ÅK, Sandling JK, Li LA, Irvin MR, Zhi D, Deloukas P, Liang L, Liu C, Bressler J, Spector TD, North K, Li Y, Absher DM, Levy D, Arnett DK, Fornage M, Pankow JS, Boerwinkle E. Epigenome-wide association study (EWAS) of BMI, BMI change and waist circumference in African American adults identifies multiple replicated loci. Hum Mol Genet 2015; 24:4464-79. [PMID: 25935004 DOI: 10.1093/hmg/ddv161] [Citation(s) in RCA: 230] [Impact Index Per Article: 25.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2015] [Accepted: 04/13/2015] [Indexed: 02/06/2023] Open
Abstract
Obesity is an important component of the pathophysiology of chronic diseases. Identifying epigenetic modifications associated with elevated adiposity, including DNA methylation variation, may point to genomic pathways that are dysregulated in numerous conditions. The Illumina 450K Bead Chip array was used to assay DNA methylation in leukocyte DNA obtained from 2097 African American adults in the Atherosclerosis Risk in Communities (ARIC) study. Mixed-effects regression models were used to test the association of methylation beta value with concurrent body mass index (BMI) and waist circumference (WC), and BMI change, adjusting for batch effects and potential confounders. Replication using whole-blood DNA from 2377 White adults in the Framingham Heart Study and CD4+ T cell DNA from 991 Whites in the Genetics of Lipid Lowering Drugs and Diet Network Study was followed by testing using adipose tissue DNA from 648 women in the Multiple Tissue Human Expression Resource cohort. Seventy-six BMI-related probes, 164 WC-related probes and 8 BMI change-related probes passed the threshold for significance in ARIC (P < 1 × 10(-7); Bonferroni), including probes in the recently reported HIF3A, CPT1A and ABCG1 regions. Replication using blood DNA was achieved for 37 BMI probes and 1 additional WC probe. Sixteen of these also replicated in adipose tissue, including 15 novel methylation findings near genes involved in lipid metabolism, immune response/cytokine signaling and other diverse pathways, including LGALS3BP, KDM2B, PBX1 and BBS2, among others. Adiposity traits are associated with DNA methylation at numerous CpG sites that replicate across studies despite variation in tissue type, ethnicity and analytic approaches.
Collapse
Affiliation(s)
- Ellen W Demerath
- Division of Epidemiology and Community Health, School of Public Health, University of Minnesota, Minneapolis, MN 55454, USA,
| | - Weihua Guan
- Division of Biostatistics, School of Public Health, University of Minnesota, Minneapolis, MN 55454, USA
| | - Megan L Grove
- Human Genetics Center, School of Public Health, University of Texas Health Sciences Center at Houston, Houston, TX 77030, USA
| | | | - Michael Mendelson
- Population Sciences Branch, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, MD 20824, USA, Framingham Heart Study, Framingham, MA 01702, USA, Department of Cardiology, Boston Children's Hospital, Boston, MA 02215, USA
| | - Yi-Hui Zhou
- Department of Statistics, North Carolina State University, Raleigh, NC 27695, USA
| | - Åsa K Hedman
- Wellcome Trust Centre for Human Genetics, University of Oxford, Oxford, UK, Department of Medical Sciences, Molecular Medicine and Science for Life Laboratory, Uppsala University, Uppsala, Sweden
| | - Johanna K Sandling
- Department of Medical Sciences, Molecular Medicine and Science for Life Laboratory, Uppsala University, Uppsala, Sweden, Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, UK
| | - Li-An Li
- Human Genome Sequencing Center, Baylor College of Medicine, Houston, TX 77030, USA
| | | | - Degui Zhi
- Department of Biostatistics, School of Public Health, University of Alabama at Birmingham, Birmingham, AL 35294, USA
| | - Panos Deloukas
- Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, UK, William Harvey Research Institute, Barts and The London School of Medicine and Dentistry, Queen Mary University of London, London, UK, Princess Al-Jawhara Al-Brahim Centre of Excellence in Research of Hereditary Disorders (PACER-HD), King Abdulaziz University, Jeddah 21589, Saudi Arabia
| | - Liming Liang
- Population Sciences Branch, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, MD 20824, USA, Framingham Heart Study, Framingham, MA 01702, USA, Departments of Epidemiology and Biostatistics, School of Public Health, Harvard University, Boston, MA 02115, USA
| | - Chunyu Liu
- Framingham Heart Study, Framingham, MA 01702, USA, Department of Biostatistics, Boston University, Boston, MA 02118, USA
| | - Jan Bressler
- Human Genetics Center, School of Public Health, University of Texas Health Sciences Center at Houston, Houston, TX 77030, USA
| | - Tim D Spector
- Department of Twin Research and Genetic Epidemiology, King's College London, London SE1 7EH, UK
| | - Kari North
- Department of Epidemiology, University of North Carolina at Chapel Hill, Chapel Hill, NC 27514, USA
| | - Yun Li
- Department of Genetics, Department of Biostatistics and Department of Computer Science, University of North Carolina, Chapel Hill, NC 27599, USA
| | - Devin M Absher
- Hudson Alpha Institute for Biotechnology, Huntsville, AL 34806, USA
| | - Daniel Levy
- Population Sciences Branch, National Heart, Lung, and Blood Institute, National Institutes of Health, Bethesda, MD 20824, USA, Framingham Heart Study, Framingham, MA 01702, USA
| | | | - Myriam Fornage
- Human Genetics Center, School of Public Health, University of Texas Health Sciences Center at Houston, Houston, TX 77030, USA, Human Genome Sequencing Center, Baylor College of Medicine, Houston, TX 77030, USA
| | - James S Pankow
- Division of Epidemiology and Community Health, School of Public Health, University of Minnesota, Minneapolis, MN 55454, USA
| | - Eric Boerwinkle
- Human Genetics Center, School of Public Health, University of Texas Health Sciences Center at Houston, Houston, TX 77030, USA, Human Genome Sequencing Center, Baylor College of Medicine, Houston, TX 77030, USA
| |
Collapse
|
55
|
Ma Z, Teschendorff AE, Leijon A, Qiao Y, Zhang H, Guo J. Variational Bayesian Matrix Factorization for Bounded Support Data. IEEE TRANSACTIONS ON PATTERN ANALYSIS AND MACHINE INTELLIGENCE 2015; 37:876-889. [PMID: 26353300 DOI: 10.1109/tpami.2014.2353639] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
A novel Bayesian matrix factorization method for bounded support data is presented. Each entry in the observation matrix is assumed to be beta distributed. As the beta distribution has two parameters, two parameter matrices can be obtained, which matrices contain only nonnegative values. In order to provide low-rank matrix factorization, the nonnegative matrix factorization (NMF) technique is applied. Furthermore, each entry in the factorized matrices, i.e., the basis and excitation matrices, is assigned with gamma prior. Therefore, we name this method as beta-gamma NMF (BG-NMF). Due to the integral expression of the gamma function, estimation of the posterior distribution in the BG-NMF model can not be presented by an analytically tractable solution. With the variational inference framework and the relative convexity property of the log-inverse-beta function, we propose a new lower-bound to approximate the objective function. With this new lower-bound, we derive an analytically tractable solution to approximately calculate the posterior distributions. Each of the approximated posterior distributions is also gamma distributed, which retains the conjugacy of the Bayesian estimation. In addition, a sparse BG-NMF can be obtained by including a sparseness constraint to the gamma prior. Evaluations with synthetic data and real life data demonstrate the good performance of the proposed method.
Collapse
|
56
|
Ali O, Cerjak D, Kent JW, James R, Blangero J, Carless MA, Zhang Y. An epigenetic map of age-associated autosomal loci in northern European families at high risk for the metabolic syndrome. Clin Epigenetics 2015; 7:12. [PMID: 25806089 PMCID: PMC4372177 DOI: 10.1186/s13148-015-0048-6] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2014] [Accepted: 01/16/2015] [Indexed: 12/20/2022] Open
Abstract
BACKGROUND The prevalence of chronic diseases such as cancer, type 2 diabetes, metabolic syndrome (MetS), and cardiovascular disease increases with age in all populations. Epigenetic features are hypothesized to play important roles in the pathophysiology of age-associated diseases, but a map of these markers is lacking. We searched for genome-wide age-associated methylation signatures in peripheral blood of individuals at high risks for MetS by profiling 485,000 CpG sites in 192 individuals of Northern European ancestry using the Illumina HM450 array. Subjects (ages 6-85 years) were part of seven extended families, and 73% of adults and 32% of children were overweight or obese. RESULTS We found 22,122 genome-wide significant age-associated CpG sites (P α=0.05 = 3.65 × 10(-7) after correction for multiple testing) of which 14,155 are positively associated with age while 7,967 are negatively associated. By applying a positional density-based clustering algorithm, we generated a map of epigenetic 'hot-spots' of age-associated genomic segments, which include 290 age-associated differentially methylated CpG clusters (aDMCs), of which 207 are positively associated with age. Gene/pathway enrichment analyses were performed on these clusters using FatiGO. Genes localized to both the positively (n = 241) and negatively (n = 16) age-associated clusters are significantly enriched in specific KEGG pathways and GO terms. The most significantly enriched pathways are the hedgehog signaling pathway (adjusted P = 3.96 × 10(-3)) and maturity-onset diabetes of the young (MODY) (adjusted P = 6.26 × 10(-3)) in the positive aDMCs and type I diabetes mellitus (adjusted P = 3.69 × 10(-7)) in the negative aDMCs. We also identified several epigenetic loci whose age-associated change rates differ between subjects diagnosed with MetS and those without. CONCLUSION We conclude that in a family cohort at high risk for MetS, age-associated epigenetic features enrich in biological pathways important for determining the fate of fat cells and for insulin production. We also observe that several genes known to be related to MetS show differential epigenetic response to age in individuals with and without MetS.
Collapse
Affiliation(s)
- Omar Ali
- />Department of Pediatrics, Medical College of Wisconsin, Milwaukee, Wisconsin USA
| | - Diana Cerjak
- />TOPS Obesity and Metabolic Research Center, Department of Medicine, Medical College of Wisconsin, Milwaukee, Wisconsin USA
- />Human and Molecular Genetics Center, Medical College of Wisconsin, 8701 Watertown Plank Road, Milwaukee, WI 53226, Wisconsin USA
| | - Jack W Kent
- />Department of Genetics, Texas Biomedical Research Institute, San Antonio, Texas USA
| | - Roland James
- />TOPS Obesity and Metabolic Research Center, Department of Medicine, Medical College of Wisconsin, Milwaukee, Wisconsin USA
- />Human and Molecular Genetics Center, Medical College of Wisconsin, 8701 Watertown Plank Road, Milwaukee, WI 53226, Wisconsin USA
| | - John Blangero
- />Department of Genetics, Texas Biomedical Research Institute, San Antonio, Texas USA
| | - Melanie A Carless
- />Department of Genetics, Texas Biomedical Research Institute, San Antonio, Texas USA
| | - Yi Zhang
- />TOPS Obesity and Metabolic Research Center, Department of Medicine, Medical College of Wisconsin, Milwaukee, Wisconsin USA
- />Human and Molecular Genetics Center, Medical College of Wisconsin, 8701 Watertown Plank Road, Milwaukee, WI 53226, Wisconsin USA
| |
Collapse
|
57
|
Ng JMK, Yu J. Promoter hypermethylation of tumour suppressor genes as potential biomarkers in colorectal cancer. Int J Mol Sci 2015; 16:2472-96. [PMID: 25622259 PMCID: PMC4346847 DOI: 10.3390/ijms16022472] [Citation(s) in RCA: 132] [Impact Index Per Article: 14.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2014] [Revised: 12/02/2014] [Accepted: 01/14/2015] [Indexed: 02/06/2023] Open
Abstract
Colorectal cancer (CRC) is a common malignancy and the fourth leading cause of cancer deaths worldwide. It results from the accumulation of multiple genetic and epigenetic changes leading to the transformation of colon epithelial cells into invasive adenocarcinomas. In CRC, epigenetic changes, in particular promoter CpG island methylation, occur more frequently than genetic mutations. Hypermethylation contributes to carcinogenesis by inducing transcriptional silencing or downregulation of tumour suppressor genes and currently, over 600 candidate hypermethylated genes have been identified. Over the past decade, a deeper understanding of epigenetics coupled with technological advances have hinted at the potential of translating benchtop research into biomarkers for clinical use. DNA methylation represents one of the largest bodies of literature in epigenetics, and hence has the highest potential for minimally invasive biomarker development. Most progress has been made in the development of diagnostic markers and there are currently two, one stool-based and one blood-based, biomarkers that are commercially available for diagnostics. Prognostic and predictive methylation markers are still at their infantile stages.
Collapse
Affiliation(s)
- Jennifer Mun-Kar Ng
- Institute of Digestive Disease and Department of Medicine and Therapeutics, State Key Laboratory of Digestive Disease, Li Ka Shing Institute of Health Sciences, CUHK Shenzhen Research Institute, The Chinese University of Hong Kong, Hong Kong, China
| | - Jun Yu
- Institute of Digestive Disease and Department of Medicine and Therapeutics, State Key Laboratory of Digestive Disease, Li Ka Shing Institute of Health Sciences, CUHK Shenzhen Research Institute, The Chinese University of Hong Kong, Hong Kong, China.
| |
Collapse
|
58
|
|
59
|
Abstract
Just as genome-wide association studies (GWAS) grew from the field of genetic epidemiology, so too do epigenome-wide association studies (EWAS) derive from the burgeoning field of epigenetic epidemiology, with both aiming to understand the molecular basis for disease risk. While genetic risk of disease is currently unmodifiable, there is hope that epigenetic risk may be reversible and or modifiable. This review will take a look back at the origins of this field and revisit the past early efforts to conduct EWAS using the 27k Illumina methylation beadarrays, to the present where most investigators are using the 450k Illumina beadarrays and finally to the future where next generation sequencing based methods beckon. There have been numerous diseases, exposures and lifestyle factors investigated with EWAS, with several significant associations now identified. However, much like the GWAS studies, EWAS are likely to require large international consortium-based approaches to reach the numbers of subjects, and statistical and scientific rigor, required for robust findings.
Collapse
Affiliation(s)
- James M Flanagan
- Epigenetics Unit, Division of Cancer, Department of Surgery and Cancer, Faculty of Medicine, Imperial College London, 4th Floor IRDB, Hammersmith Campus, Du Cane Road, London, W12 0NN, UK,
| |
Collapse
|
60
|
Wahl S, Fenske N, Zeilinger S, Suhre K, Gieger C, Waldenberger M, Grallert H, Schmid M. On the potential of models for location and scale for genome-wide DNA methylation data. BMC Bioinformatics 2014; 15:232. [PMID: 24994026 PMCID: PMC4227139 DOI: 10.1186/1471-2105-15-232] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2014] [Accepted: 06/26/2014] [Indexed: 01/05/2023] Open
Abstract
Background With the help of epigenome-wide association studies (EWAS), increasing knowledge on the role of epigenetic mechanisms such as DNA methylation in disease processes is obtained. In addition, EWAS aid the understanding of behavioral and environmental effects on DNA methylation. In terms of statistical analysis, specific challenges arise from the characteristics of methylation data. First, methylation β-values represent proportions with skewed and heteroscedastic distributions. Thus, traditional modeling strategies assuming a normally distributed response might not be appropriate. Second, recent evidence suggests that not only mean differences but also variability in site-specific DNA methylation associates with diseases, including cancer. The purpose of this study was to compare different modeling strategies for methylation data in terms of model performance and performance of downstream hypothesis tests. Specifically, we used the generalized additive models for location, scale and shape (GAMLSS) framework to compare beta regression with Gaussian regression on raw, binary logit and arcsine square root transformed methylation data, with and without modeling a covariate effect on the scale parameter. Results Using simulated and real data from a large population-based study and an independent sample of cancer patients and healthy controls, we show that beta regression does not outperform competing strategies in terms of model performance. In addition, Gaussian models for location and scale showed an improved performance as compared to models for location only. The best performance was observed for the Gaussian model on binary logit transformed β-values, referred to as M-values. Our results further suggest that models for location and scale are specifically sensitive towards violations of the distribution assumption and towards outliers in the methylation data. Therefore, a resampling procedure is proposed as a mode of inference and shown to diminish type I error rate in practically relevant settings. We apply the proposed method in an EWAS of BMI and age and reveal strong associations of age with methylation variability that are validated in an independent sample. Conclusions Models for location and scale are promising tools for EWAS that may help to understand the influence of environmental factors and disease-related phenotypes on methylation variability and its role during disease development.
Collapse
Affiliation(s)
- Simone Wahl
- Research Unit of Molecular Epidemiology, Helmholtz Zentrum München - German Research Center for Environmental Health, Ingolstädter Landstrasse 1, 85764 Neuherberg, Germany.
| | | | | | | | | | | | | | | |
Collapse
|
61
|
Ching T, Song MA, Tiirikainen M, Molnar J, Berry M, Towner D, Garmire LX. Genome-wide hypermethylation coupled with promoter hypomethylation in the chorioamniotic membranes of early onset pre-eclampsia. Mol Hum Reprod 2014; 20:885-904. [PMID: 24944161 DOI: 10.1093/molehr/gau046] [Citation(s) in RCA: 42] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
Pre-eclampsia is the leading cause of fetal and maternal morbidity and mortality. Early onset pre-eclampsia (EOPE) is a disorder that has severe maternal and fetal outcomes, whilst its etiology is poorly understood. We hypothesize that epigenetics plays an important role to mediate the development of EOPE and conducted a case-control study to compare the genome-wide methylome difference between chorioamniotic membranes from 30 EOPE and 17 full-term pregnancies using the Infinium Human Methylation 450 BeadChip arrays. Bioinformatics analysis tested differential methylation (DM) at CpG site level, gene level, and pathway and network level. A striking genome-wide hypermethylation pattern coupled with hypomethylation in promoters was observed. Out of 385 184 CpG sites, 9995 showed DM (2.6%). Of those DM sites, 91.9% showed hypermethylation (9186 of 9995). Over 900 genes had DM associated with promoters. Promoter-based DM analysis revealed that genes in canonical cancer-related pathways such as Rac, Ras, PI3K/Akt, NFκB and ErBB4 were enriched, and represented biological functional alterations that involve cell cycle, apoptosis, cancer signaling and inflammation. A group of genes previously found to be up-regulated in pre-eclampsia, including GRB2, ATF3, NFKB2, as well as genes in proteasome subunits (PSMA1, PMSE1, PSMD1 and PMSD8), harbored hypomethylated promoters. Contrarily, a cluster of microRNAs, including mir-519a1, mir-301a, mir-487a, mir-185, mir-329, mir-194, mir-376a1, mir-486 and mir-744 were all hypermethylated in their promoters in the EOPE samples. These findings collectively reveal new avenues of research regarding the vast epigenetic modifications in EOPE.
Collapse
Affiliation(s)
- Travers Ching
- Molecular Bioscience and Bioengineering Graduate Program, University of Hawaii at Manoa, Honolulu, HI 96822, USA Epidemiology Program, University of Hawaii Cancer Center, Honolulu, HI 96813, USA
| | - Min-Ae Song
- Molecular Bioscience and Bioengineering Graduate Program, University of Hawaii at Manoa, Honolulu, HI 96822, USA Genomics Shared Resources Core, University of Hawaii Cancer Center, Honolulu, HI 96813, USA
| | - Maarit Tiirikainen
- Genomics Shared Resources Core, University of Hawaii Cancer Center, Honolulu, HI 96813, USA
| | - Janos Molnar
- Genomics Shared Resources Core, University of Hawaii Cancer Center, Honolulu, HI 96813, USA
| | - Marla Berry
- Department of Cell and Molecular Biology, John A. Burns School of Medicine, University of Hawaii, Honolulu, HI 96813, USA
| | - Dena Towner
- Department of Obstetrics, Gynecology and Women's Health, John A. Burns School of Medicine, University of Hawaii, Honolulu, HI 96826, USA
| | - Lana X Garmire
- Molecular Bioscience and Bioengineering Graduate Program, University of Hawaii at Manoa, Honolulu, HI 96822, USA Epidemiology Program, University of Hawaii Cancer Center, Honolulu, HI 96813, USA
| |
Collapse
|
62
|
Ma Z, Teschendorff AE, Yu H, Taghia J, Guo J. Comparisons of non-Gaussian statistical models in DNA methylation analysis. Int J Mol Sci 2014; 15:10835-54. [PMID: 24937687 PMCID: PMC4100184 DOI: 10.3390/ijms150610835] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2014] [Revised: 05/12/2014] [Accepted: 06/10/2014] [Indexed: 12/25/2022] Open
Abstract
As a key regulatory mechanism of gene expression, DNA methylation patterns are widely altered in many complex genetic diseases, including cancer. DNA methylation is naturally quantified by bounded support data; therefore, it is non-Gaussian distributed. In order to capture such properties, we introduce some non-Gaussian statistical models to perform dimension reduction on DNA methylation data. Afterwards, non-Gaussian statistical model-based unsupervised clustering strategies are applied to cluster the data. Comparisons and analysis of different dimension reduction strategies and unsupervised clustering methods are presented. Experimental results show that the non-Gaussian statistical model-based methods are superior to the conventional Gaussian distribution-based method. They are meaningful tools for DNA methylation analysis. Moreover, among several non-Gaussian methods, the one that captures the bounded nature of DNA methylation data reveals the best clustering performance.
Collapse
Affiliation(s)
- Zhanyu Ma
- Pattern Recognition and Intelligent System Lab.,Beijing University of Posts and Telecommunications, No. 10 Xitucheng Road,Beijing 100876, China.
| | - Andrew E Teschendorff
- Computational Systems Genomics, CAS-MPG Partner Institute for Computational Biology, Shanghai Institute for Biological Sciences, Chinese Academy of Sciences, 320 Yue Yang Road, Shanghai 200031, China.
| | - Hong Yu
- Pattern Recognition and Intelligent System Lab.,Beijing University of Posts and Telecommunications, No. 10 Xitucheng Road,Beijing 100876, China.
| | - Jalil Taghia
- Communication Theory Lab., KTH - Royal Institute of Technology, Osquldas väg 10,10044 Stockholm, Sweden.
| | - Jun Guo
- Pattern Recognition and Intelligent System Lab.,Beijing University of Posts and Telecommunications, No. 10 Xitucheng Road,Beijing 100876, China.
| |
Collapse
|
63
|
Hsiao CL, Hsieh AR, Lian IB, Lin YC, Wang HM, Fann CSJ. A novel method for identification and quantification of consistently differentially methylated regions. PLoS One 2014; 9:e97513. [PMID: 24818602 PMCID: PMC4018258 DOI: 10.1371/journal.pone.0097513] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2014] [Accepted: 04/16/2014] [Indexed: 12/28/2022] Open
Abstract
Advances in biotechnology have resulted in large-scale studies of DNA methylation. A differentially methylated region (DMR) is a genomic region with multiple adjacent CpG sites that exhibit different methylation statuses among multiple samples. Many so-called “supervised” methods have been established to identify DMRs between two or more comparison groups. Methods for the identification of DMRs without reference to phenotypic information are, however, less well studied. An alternative “unsupervised” approach was proposed, in which DMRs in studied samples were identified with consideration of nature dependence structure of methylation measurements between neighboring probes from tiling arrays. Through simulation study, we investigated effects of dependencies between neighboring probes on determining DMRs where a lot of spurious signals would be produced if the methylation data were analyzed independently of the probe. In contrast, our newly proposed method could successfully correct for this effect with a well-controlled false positive rate and a comparable sensitivity. By applying to two real datasets, we demonstrated that our method could provide a global picture of methylation variation in studied samples. R source codes to implement the proposed method were freely available at http://www.csjfann.ibms.sinica.edu.tw/eag/programlist/ICDMR/ICDMR.html.
Collapse
Affiliation(s)
- Ching-Lin Hsiao
- Institute of Biomedical Sciences, Academia Sinica, Taipei, Taiwan
| | - Ai-Ru Hsieh
- Institute of Biomedical Sciences, Academia Sinica, Taipei, Taiwan
| | - Ie-Bin Lian
- Department of Mathematics, National Changhua University of Education, Changhua, Taiwan
| | - Ying-Chao Lin
- Institute of Biomedical Sciences, Academia Sinica, Taipei, Taiwan
| | - Hui-Min Wang
- Institute of Biomedical Sciences, Academia Sinica, Taipei, Taiwan
| | - Cathy S. J. Fann
- Institute of Biomedical Sciences, Academia Sinica, Taipei, Taiwan
- * E-mail:
| |
Collapse
|
64
|
Jiao Y, Widschwendter M, Teschendorff AE. A systems-level integrative framework for genome-wide DNA methylation and gene expression data identifies differential gene expression modules under epigenetic control. Bioinformatics 2014; 30:2360-6. [DOI: 10.1093/bioinformatics/btu316] [Citation(s) in RCA: 116] [Impact Index Per Article: 11.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
|
65
|
Vucic EA, Chari R, Thu KL, Wilson IM, Cotton AM, Kennett JY, Zhang M, Lonergan KM, Steiling K, Brown CJ, McWilliams A, Ohtani K, Lenburg ME, Sin DD, Spira A, MacAulay CE, Lam S, Lam WL. DNA methylation is globally disrupted and associated with expression changes in chronic obstructive pulmonary disease small airways. Am J Respir Cell Mol Biol 2014; 50:912-22. [PMID: 24298892 PMCID: PMC4068945 DOI: 10.1165/rcmb.2013-0304oc] [Citation(s) in RCA: 96] [Impact Index Per Article: 9.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2013] [Accepted: 12/03/2013] [Indexed: 01/06/2023] Open
Abstract
DNA methylation is an epigenetic modification that is highly disrupted in response to cigarette smoke and involved in a wide spectrum of malignant and nonmalignant diseases, but surprisingly not previously assessed in small airways of patients with chronic obstructive pulmonary disease (COPD). Small airways are the primary sites of airflow obstruction in COPD. We sought to determine whether DNA methylation patterns are disrupted in small airway epithelia of patients with COPD, and evaluate whether changes in gene expression are associated with these disruptions. Genome-wide methylation and gene expression analysis were performed on small airway epithelial DNA and RNA obtained from the same patient during bronchoscopy, using Illumina's Infinium HM27 and Affymetrix's Genechip Human Gene 1.0 ST arrays. To control for known effects of cigarette smoking on DNA methylation, methylation and gene expression profiles were compared between former smokers with and without COPD matched for age, pack-years, and years of smoking cessation. Our results indicate that aberrant DNA methylation is (1) a genome-wide phenomenon in small airways of patients with COPD, and (2) associated with altered expression of genes and pathways important to COPD, such as the NF-E2-related factor 2 oxidative response pathway. DNA methylation is likely an important mechanism contributing to modulation of genes important to COPD pathology. Because these methylation events may underlie disease-specific gene expression changes, their characterization is a critical first step toward the development of epigenetic markers and an opportunity for developing novel epigenetic therapeutic interventions for COPD.
Collapse
Affiliation(s)
- Emily A. Vucic
- Department of Integrative Oncology, British Columbia Cancer Research Centre, Vancouver, British Columbia, Canada
| | - Raj Chari
- Department of Genetics, Harvard Medical School, Boston, Massachusetts
| | - Kelsie L. Thu
- Department of Integrative Oncology, British Columbia Cancer Research Centre, Vancouver, British Columbia, Canada
| | - Ian M. Wilson
- Department of Integrative Oncology, British Columbia Cancer Research Centre, Vancouver, British Columbia, Canada
| | - Allison M. Cotton
- Department of Medical Genetics, University of British Columbia, Vancouver, British Columbia, Canada
| | - Jennifer Y. Kennett
- Department of Integrative Oncology, British Columbia Cancer Research Centre, Vancouver, British Columbia, Canada
| | - May Zhang
- Department of Integrative Oncology, British Columbia Cancer Research Centre, Vancouver, British Columbia, Canada
| | - Kim M. Lonergan
- Department of Integrative Oncology, British Columbia Cancer Research Centre, Vancouver, British Columbia, Canada
| | - Katrina Steiling
- Division of Computational Biomedicine, Department of Medicine, Boston University Medical Center, Boston, Massachusetts; and
| | - Carolyn J. Brown
- Department of Medical Genetics, University of British Columbia, Vancouver, British Columbia, Canada
| | - Annette McWilliams
- Department of Integrative Oncology, British Columbia Cancer Research Centre, Vancouver, British Columbia, Canada
| | - Keishi Ohtani
- Department of Integrative Oncology, British Columbia Cancer Research Centre, Vancouver, British Columbia, Canada
| | - Marc E. Lenburg
- Division of Computational Biomedicine, Department of Medicine, Boston University Medical Center, Boston, Massachusetts; and
| | - Don D. Sin
- University of British Columbia James Hogg Research Centre and the Institute of Heart and Lung Health, St. Paul’s Hospital, Vancouver, British Columbia, Canada
| | - Avrum Spira
- Division of Computational Biomedicine, Department of Medicine, Boston University Medical Center, Boston, Massachusetts; and
| | - Calum E. MacAulay
- Department of Integrative Oncology, British Columbia Cancer Research Centre, Vancouver, British Columbia, Canada
| | - Stephen Lam
- Department of Integrative Oncology, British Columbia Cancer Research Centre, Vancouver, British Columbia, Canada
| | - Wan L. Lam
- Department of Integrative Oncology, British Columbia Cancer Research Centre, Vancouver, British Columbia, Canada
| |
Collapse
|
66
|
Chou WC, Cheng AL, Brotto M, Chuang CY. Visual gene-network analysis reveals the cancer gene co-expression in human endometrial cancer. BMC Genomics 2014; 15:300. [PMID: 24758163 PMCID: PMC4234489 DOI: 10.1186/1471-2164-15-300] [Citation(s) in RCA: 78] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2013] [Accepted: 04/04/2014] [Indexed: 11/10/2022] Open
Abstract
Background Endometrial cancers (ECs) are the most common form of gynecologic malignancy. Recent studies have reported that ECs reveal distinct markers for molecular pathogenesis, which in turn is linked to the various histological types of ECs. To understand further the molecular events contributing to ECs and endometrial tumorigenesis in general, a more precise identification of cancer-associated molecules and signaling networks would be useful for the detection and monitoring of malignancy, improving clinical cancer therapy, and personalization of treatments. Results ECs-specific gene co-expression networks were constructed by differential expression analysis and weighted gene co-expression network analysis (WGCNA). Important pathways and putative cancer hub genes contribution to tumorigenesis of ECs were identified. An elastic-net regularized classification model was built using the cancer hub gene signatures to predict the phenotypic characteristics of ECs. The 19 cancer hub gene signatures had high predictive power to distinguish among three key principal features of ECs: grade, type, and stage. Intriguingly, these hub gene networks seem to contribute to ECs progression and malignancy via cell-cycle regulation, antigen processing and the citric acid (TCA) cycle. Conclusions The results of this study provide a powerful biomarker discovery platform to better understand the progression of ECs and to uncover potential therapeutic targets in the treatment of ECs. This information might lead to improved monitoring of ECs and resulting improvement of treatment of ECs, the 4th most common of cancer in women.
Collapse
Affiliation(s)
| | | | | | - Chun-Yu Chuang
- Department of Biomedical Engineering and Environmental Sciences, National Tsing Hua University, Hsinchu 30013, Taiwan.
| |
Collapse
|
67
|
Aphinyanaphongs Y, Fu LD, Li Z, Peskin ER, Efstathiadis E, Aliferis CF, Statnikov A. A comprehensive empirical comparison of modern supervised classification and feature selection methods for text categorization. J Assoc Inf Sci Technol 2014. [DOI: 10.1002/asi.23110] [Citation(s) in RCA: 40] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Affiliation(s)
- Yindalon Aphinyanaphongs
- Center for Health Informatics and Bioinformatics; New York University Langone Medical Center; 227 East 30th Street New York NY 10016
- Department of Medicine; New York University School of Medicine; 550 First Avenue New York NY 10016
| | - Lawrence D. Fu
- Center for Health Informatics and Bioinformatics; New York University Langone Medical Center; 227 East 30th Street New York NY 10016
- Department of Medicine; New York University School of Medicine; 550 First Avenue New York NY 10016
| | - Zhiguo Li
- Center for Health Informatics and Bioinformatics; New York University Langone Medical Center; 227 East 30th Street New York NY 10016
| | - Eric R. Peskin
- Center for Health Informatics and Bioinformatics; New York University Langone Medical Center; 227 East 30th Street New York NY 10016
| | - Efstratios Efstathiadis
- Center for Health Informatics and Bioinformatics; New York University Langone Medical Center; 227 East 30th Street New York NY 10016
| | - Constantin F. Aliferis
- Center for Health Informatics and Bioinformatics; New York University Langone Medical Center; 227 East 30th Street New York NY 10016
- Department of Pathology; New York University School of Medicine; 550 First Avenue New York NY 10016
- Department of Biostatistics; Vanderbilt University; 1211 Medical Center Drive Nashville TN 37232
| | - Alexander Statnikov
- Center for Health Informatics and Bioinformatics; New York University Langone Medical Center; 227 East 30th Street New York NY 10016
- Department of Medicine; New York University School of Medicine; 550 First Avenue New York NY 10016
| |
Collapse
|
68
|
Morris TJ, Butcher LM, Feber A, Teschendorff AE, Chakravarthy AR, Wojdacz TK, Beck S. ChAMP: 450k Chip Analysis Methylation Pipeline. Bioinformatics 2014; 30:428-30. [PMID: 24336642 PMCID: PMC3904520 DOI: 10.1093/bioinformatics/btt684] [Citation(s) in RCA: 632] [Impact Index Per Article: 63.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2013] [Revised: 11/03/2013] [Accepted: 11/19/2013] [Indexed: 12/18/2022] Open
Abstract
UNLABELLED The Illumina Infinium HumanMethylation450 BeadChip is a new platform for high-throughput DNA methylation analysis. Several methods for normalization and processing of these data have been published recently. Here we present an integrated analysis pipeline offering a choice of the most popular normalization methods while also introducing new methods for calling differentially methylated regions and detecting copy number aberrations. AVAILABILITY AND IMPLEMENTATION ChAMP is implemented as a Bioconductor package in R. The package and the vignette can be downloaded at bioconductor.org
Collapse
Affiliation(s)
- Tiffany J Morris
- Medical Genomics, and Statistical Genomics, UCL Cancer Institute, University College London, London WC1E 6BT, UK, CAS-MPG Partner Institute for Computational Biology, Shanghai Institute for Biological Sciences, Shanghai 200031, China and Institute of Environmental Medicine, Karolinska Instytutet, 17177 Stockholm, Sweden
| | | | | | | | | | | | | |
Collapse
|
69
|
Epigenetic analysis of neurocognitive development at 1 year of age in a community-based pregnancy cohort. Behav Genet 2014; 44:113-25. [PMID: 24452678 DOI: 10.1007/s10519-014-9641-2] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2013] [Accepted: 01/09/2014] [Indexed: 10/25/2022]
Abstract
Multiple studies show that molecular genetic changes and epigenetic modifications affect the risk of cognitive disability or impairment. However, the role of epigenetic variation in cognitive development of neurotypical young children remains largely unknown. Using data from a prospective, community-based study of mother-infant pairs, we investigated the association of DNA methylation patterns in neonatal umbilical cord blood with cognitive and language development at 1 year of age. No CpG loci achieved genome-wide significance, although a small number of weakly suggestive associations with Bayley-III Receptive Communication scales were noted. While umbilical cord blood is a convenient resource for genetic analyses of birth outcomes, our results do not provide conclusive evidence that its use for DNA methylation profiling yields epigenetic markers that are directly related to postnatal neurocognitive outcomes at 1 year of age.
Collapse
|
70
|
Harlaar N, Hutchison KE. Alcohol and the methylome: design and analysis considerations for research using human samples. Drug Alcohol Depend 2013; 133:305-16. [PMID: 23968814 DOI: 10.1016/j.drugalcdep.2013.07.026] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 01/22/2013] [Revised: 07/26/2013] [Accepted: 07/26/2013] [Indexed: 12/14/2022]
Abstract
BACKGROUND A growing number of studies in human samples have sought to determine whether chronic alcohol use and alcohol use disorders (AUDs) may be associated with epigenetic factors, such as DNA methylation. We review the extant literature in light of some of the challenges that currently affect the design and interpretation of epigenetic research in human samples. METHOD A literature search was used to identify studies that have examined DNA methylation in relation to alcohol use or AUDs in human samples (through July 2013). A total of 22 studies were identified. RESULTS Associations with quantitative or diagnostic phenotypes of alcohol use or AUDs have been reported for several genes. However, all studies to date have relied on relatively small samples and cross-sectional study designs. Additionally, attempts to replicate results have been rare. More generally, research progress is hampered by several issues, including limitations of the technologies used to assess DNA methylation, tissue- and cell-specificity of methylation patterns, the difficulties of relating observed methylation differences at a given locus to a functional effect, and limited knowledge about the molecular mechanisms underlying the effects of alcohol on DNA methylation. CONCLUSIONS Although we share the optimism that epigenetics may lead to new insights into the etiology and pathophysiology of AUDs, the methodological and scientific challenges associated with conducting methylomic research in human samples need to be carefully considered when designing and evaluating such studies.
Collapse
Affiliation(s)
- Nicole Harlaar
- University of Colorado Boulder, Boulder, CO 80309-0345, USA.
| | | |
Collapse
|
71
|
Wilhelm-Benartzi CS, Koestler DC, Karagas MR, Flanagan JM, Christensen BC, Kelsey KT, Marsit CJ, Houseman EA, Brown R. Review of processing and analysis methods for DNA methylation array data. Br J Cancer 2013; 109:1394-402. [PMID: 23982603 PMCID: PMC3777004 DOI: 10.1038/bjc.2013.496] [Citation(s) in RCA: 129] [Impact Index Per Article: 11.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/24/2013] [Revised: 07/23/2013] [Accepted: 07/30/2013] [Indexed: 12/21/2022] Open
Abstract
The promise of epigenome-wide association studies and cancer-specific somatic DNA methylation changes in improving our understanding of cancer, coupled with the decreasing cost and increasing coverage of DNA methylation microarrays, has brought about a surge in the use of these technologies. Here, we aim to provide both a review of issues encountered in the processing and analysis of array-based DNA methylation data and a summary of the advantages of recent approaches proposed for handling those issues, focusing on approaches publicly available in open-source environments such as R and Bioconductor. We hope that the processing tools and analysis flowchart described herein will facilitate researchers to effectively use these powerful DNA methylation array-based platforms, thereby advancing our understanding of human health and disease.
Collapse
Affiliation(s)
- C S Wilhelm-Benartzi
- Epigenetics Unit, Division of Cancer, Department of Surgery and Cancer, Faculty of Medicine, Ovarian Cancer Action Research Centre, Imperial College London, 4th floor IRDB, Hammersmith Campus, Du Cane Road, London W12 0NN, UK
| | - D C Koestler
- Section of Biostatistics and Epidemiology, Geisel School of Medicine at Dartmouth College, Hanover, NH 03755, USA
| | - M R Karagas
- Section of Biostatistics and Epidemiology, Geisel School of Medicine at Dartmouth College, Hanover, NH 03755, USA
| | - J M Flanagan
- Epigenetics Unit, Division of Cancer, Department of Surgery and Cancer, Faculty of Medicine, Ovarian Cancer Action Research Centre, Imperial College London, 4th floor IRDB, Hammersmith Campus, Du Cane Road, London W12 0NN, UK
| | - B C Christensen
- Section of Biostatistics and Epidemiology, Geisel School of Medicine at Dartmouth College, Hanover, NH 03755, USA
- Department of Pharmacology and Toxicology, Geisel School of Medicine at Dartmouth College, Hanover, NH 03755, USA
| | - K T Kelsey
- Department of Pathology and Laboratory Medicine, Brown University, Providence, RI, USA
- Department of Epidemiology, Brown University, Providence, RI, USA
| | - C J Marsit
- Section of Biostatistics and Epidemiology, Geisel School of Medicine at Dartmouth College, Hanover, NH 03755, USA
- Department of Pharmacology and Toxicology, Geisel School of Medicine at Dartmouth College, Hanover, NH 03755, USA
| | - E A Houseman
- Department of Public Health, Oregon State University, Corvallis, OR, USA
| | - R Brown
- Epigenetics Unit, Division of Cancer, Department of Surgery and Cancer, Faculty of Medicine, Ovarian Cancer Action Research Centre, Imperial College London, 4th floor IRDB, Hammersmith Campus, Du Cane Road, London W12 0NN, UK
- Section of Molecular Pathology, Institute for Cancer Research, Sutton, UK
| |
Collapse
|
72
|
Li P, Demirci F, Mahalingam G, Demirci C, Nakano M, Meyers BC. An integrated workflow for DNA methylation analysis. J Genet Genomics 2013; 40:249-60. [PMID: 23706300 DOI: 10.1016/j.jgg.2013.03.010] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2013] [Revised: 03/25/2013] [Accepted: 03/25/2013] [Indexed: 10/27/2022]
Abstract
The analysis of cytosine methylation provides a new way to assess and describe epigenetic regulation at a whole-genome level in many eukaryotes. DNA methylation has a demonstrated role in the genome stability and protection, regulation of gene expression and many other aspects of genome function and maintenance. BS-seq is a relatively unbiased method for profiling the DNA methylation, with a resolution capable of measuring methylation at individual cytosines. Here we describe, as an example, a workflow to handle DNA methylation analysis, from BS-seq library preparation to the data visualization. We describe some applications for the analysis and interpretation of these data. Our laboratory provides public access to plant DNA methylation data via visualization tools available at our "Next-Gen Sequence" websites (http://mpss.udel.edu), along with small RNA, RNA-seq and other data types.
Collapse
Affiliation(s)
- Pingchuan Li
- Department of Plant & Soil Sciences, Delaware Biotechnology Institute, University of Delaware, Newark, DE 19711, USA
| | | | | | | | | | | |
Collapse
|
73
|
Ma Z, Teschendorff AE. A variational Bayes beta mixture model for feature selection in DNA methylation studies. J Bioinform Comput Biol 2013; 11:1350005. [PMID: 23859269 DOI: 10.1142/s0219720013500054] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
An increasing number of studies are using beadarrays to measure DNA methylation on a genome-wide basis. The purpose is to identify novel biomarkers in a wide range of complex genetic diseases including cancer. A common difficulty encountered in these studies is distinguishing true biomarkers from false positives. While statistical methods aimed at improving the feature selection step have been developed for gene expression, relatively few methods have been adapted to DNA methylation data, which is naturally beta-distributed. Here we explore and propose an innovative application of a recently developed variational Bayesian beta-mixture model (VBBMM) to the feature selection problem in the context of DNA methylation data generated from a highly popular beadarray technology. We demonstrate that VBBMM offers significant improvements in inference and feature selection in this type of data compared to an Expectation-Maximization (EM) algorithm, at a significantly reduced computational cost. We further demonstrate the added value of VBBMM as a feature selection and prioritization step in the context of identifying prognostic markers in breast cancer. A variational Bayesian approach to feature selection of DNA methylation profiles should thus be of value to any study undergoing large-scale DNA methylation profiling in search of novel biomarkers.
Collapse
Affiliation(s)
- Zhanyu Ma
- KTH-Royal Institute of Technology, School of Electrical Engineering, SE-100 44, Stockholm, Sweden.
| | | |
Collapse
|
74
|
Abstract
DNA methylation is an epigenetic mark that has suspected regulatory roles in a broad range of biological processes and diseases. The technology is now available for studying DNA methylation genome-wide, at a high resolution and in a large number of samples. This Review discusses relevant concepts, computational methods and software tools for analysing and interpreting DNA methylation data. It focuses not only on the bioinformatic challenges of large epigenome-mapping projects and epigenome-wide association studies but also highlights software tools that make genome-wide DNA methylation mapping more accessible for laboratories with limited bioinformatics experience.
Collapse
Affiliation(s)
- Christoph Bock
- CeMM Research Center for Molecular Medicine of the Austrian Academy of Sciences, 1090 Vienna, Austria.
| |
Collapse
|
75
|
Teschendorff AE, Marabita F, Lechner M, Bartlett T, Tegner J, Gomez-Cabrero D, Beck S. A beta-mixture quantile normalization method for correcting probe design bias in Illumina Infinium 450 k DNA methylation data. ACTA ACUST UNITED AC 2012; 29:189-96. [PMID: 23175756 PMCID: PMC3546795 DOI: 10.1093/bioinformatics/bts680] [Citation(s) in RCA: 1088] [Impact Index Per Article: 90.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
Motivation: The Illumina Infinium 450 k DNA Methylation Beadchip is a prime candidate technology for Epigenome-Wide Association Studies (EWAS). However, a difficulty associated with these beadarrays is that probes come in two different designs, characterized by widely different DNA methylation distributions and dynamic range, which may bias downstream analyses. A key statistical issue is therefore how best to adjust for the two different probe designs. Results: Here we propose a novel model-based intra-array normalization strategy for 450 k data, called BMIQ (Beta MIxture Quantile dilation), to adjust the beta-values of type2 design probes into a statistical distribution characteristic of type1 probes. The strategy involves application of a three-state beta-mixture model to assign probes to methylation states, subsequent transformation of probabilities into quantiles and finally a methylation-dependent dilation transformation to preserve the monotonicity and continuity of the data. We validate our method on cell-line data, fresh frozen and paraffin-embedded tumour tissue samples and demonstrate that BMIQ compares favourably with two competing methods. Specifically, we show that BMIQ improves the robustness of the normalization procedure, reduces the technical variation and bias of type2 probe values and successfully eliminates the type1 enrichment bias caused by the lower dynamic range of type2 probes. BMIQ will be useful as a preprocessing step for any study using the Illumina Infinium 450 k platform. Availability: BMIQ is freely available from http://code.google.com/p/bmiq/. Contact:a.teschendorff@ucl.ac.uk Supplementary information:Supplementary data are available at Bioinformatics online
Collapse
Affiliation(s)
- Andrew E Teschendorff
- Statistical Genomics Group, UCL Cancer Institute, University College London, London WC1E 6BT, UK.
| | | | | | | | | | | | | |
Collapse
|