Comethyl: A network-based methylome approach to investigate the multivariate nature of health and disease

Health outcomes are frequently shaped by difficult to dissect inter-relationships between biological, behavioral, social, and environmental factors. DNA methylation patterns reflect such multi-variate intersections, providing a rich source of novel biomarkers and insight into disease etiologies. Recent advances in whole-genome bisulfite sequencing (WGBS) enable investigation of DNA methylation over all genomic CpGs, but existing bioinformatic approaches lack accessible system-level tools. Here, we develop the R package Comethyl, for weighted gene correlation network analysis (WGCNA) of user-defined genomic regions that generates modules of comethylated regions, which are then tested for correlations with sample traits. First, regions are defined by CpG genomic location or regulatory annotation and filtered based on CpG count, sequencing depth, and variability. Next, correlation networks are used to find modules of interconnected nodes using methylation values within the selected regions. Each module containing multiple comethylated regions is reduced in complexity to a single eigennode value, which is then tested for correlations with experimental metadata. Comethyl has the ability to cover the noncoding regulatory regions of the genome with high relevance to interpretation of genome-wide association studies and integration with other types of epigenomic data. We demonstrate the utility of Comethyl on a dataset of male cord blood samples from newborns later diagnosed with autism spectrum disorder (ASD) versus typical development. Comethyl successfully identified an ASD-associated module containing gene regions with brain glial functions. Comethyl is expected to be useful in uncovering the multi-variate nature of health disparities for a variety of common disorders. Comethyl is available at github.com/cemordaunt/comethyl. Description of the AuthorsCharles E. Mordaunt, Ph.D. developed Comethyl while a postdoctoral fellow in the department of Medical Microbiology and Immunology at UC Davis. He is currently a Computational Biologist at GSK. Julia S. Mouat is a doctoral student in the Integrative Genetics and Genomics graduate group at UC Davis with interests in health disparities and intergenerational epigenetic risk factors for autism spectrum disorders. Rebecca J. Schmidt, Ph.D. is an Associate Professor of Public Health Sciences at UC Davis, with expertise in the use of epigenetics in epidemiology and neurodevelopmental disorders. Janine M. LaSalle, Ph.D. is a Professor of Medical Microbiology and Immunology, Co-Director of the Perinatal Origins of Disparities Center, and Deputy Director of the Environmental Health Sciences Center at UC Davis, with expertise in epigenomics and neurodevelopmental disorders.