Leveraging biological and statistical covariates improves the detection power in epigenome-wide association testing

Food Allergy Genetics and Epigenetics: A Review of Genome-Wide Association Studies

In this review, we provide an overview of food allergy genetics and epigenetics aimed at clinicians and researchers. This includes a brief review of the current understanding of genetic and epigenetic mechanisms, inheritance of food allergy, as well as a discussion of advantages and limitations of the different types of studies in genetic research. We specifically focus on the results of genome-wide association studies in food allergy, which have identified 16 genetic variants that reach genome-wide significance, many of which overlap with other allergic diseases, including asthma, atopic dermatitis, and allergic rhinitis. Identified genes for food allergy are mainly involved in epithelial barrier function (e.g., FLG, SERPINB7) and immune function (e.g., HLA, IL4). Epigenome-wide significant findings at 32 loci are also summarized as well as 14 additional loci with significance at a false discovery of < 1 × 10-4. Integration of epigenetic and genetic data is discussed in the context of disease mechanisms, many of which are shared with other allergic diseases. The potential utility of genetic and epigenetic discoveries is deliberated. In the future, genetic and epigenetic markers may offer ways to predict the presence or absence of clinical IgE-mediated food allergy among sensitized individuals, likelihood of development of natural tolerance, and response to immunotherapy.

Early B cell transcriptomic markers of measles-specific humoral immunity following a 3rd dose of MMR vaccine

B cell transcriptomic signatures hold promise for the early prediction of vaccine-induced humoral immunity and vaccine protective efficacy. We performed a longitudinal study in 232 healthy adult participants before/after a 3rd dose of MMR (MMR3) vaccine. We assessed baseline and early transcriptional patterns in purified B cells and their association with measles-specific humoral immunity after MMR vaccination using two analytical methods ("per gene" linear models and joint analysis). Our study identified distinct early transcriptional signatures/genes following MMR3 that were associated with measles-specific neutralizing antibody titer and/or binding antibody titer. The most significant genes included: the interleukin 20 receptor subunit beta/IL20RB gene (a subunit receptor for IL-24, a cytokine involved in the germinal center B cell maturation/response); the phorbol-12-myristate-13-acetate-induced protein 1/PMAIP1, the brain expressed X-linked 2/BEX2 gene and the B cell Fas apoptotic inhibitory molecule/FAIM, involved in the selection of high-affinity B cell clones and apoptosis/regulation of apoptosis; as well as IL16 (encoding the B lymphocyte-derived IL-16 ligand of CD4), involved in the crosstalk between B cells, dendritic cells and helper T cells. Significantly enriched pathways included B cell signaling, apoptosis/regulation of apoptosis, metabolic pathways, cell cycle-related pathways, and pathways associated with viral infections, among others. In conclusion, our study identified genes/pathways linked to antigen-induced B cell proliferation, differentiation, apoptosis, and clonal selection, that are associated with, and impact measles virus-specific humoral immunity after MMR vaccination.

2dGBH: Two-dimensional group Benjamini-Hochberg procedure for false discovery rate control in two-way multiple testing of genomic data

MOTIVATION

Emerging omics technologies have introduced a two-way grouping structure in multiple testing, as seen in single-cell omics data, where the features can be grouped by either genes or cell types. Traditional multiple testing methods have limited ability to exploit such two-way grouping structure, leading to potential power loss.

RESULTS

We propose a new 2D Group Benjamini-Hochberg (2dGBH) procedure to harness the two-way grouping structure in omics data, extending the traditional one-way adaptive GBH procedure. Using both simulated and real datasets, we show that 2dGBH effectively controls the false discovery rate across biologically relevant settings, and it is more powerful than the BH or q-value procedure and more robust than the one-way adaptive GBH procedure.

AVAILABILITY AND IMPLEMENTATION

2dGBH is available as an R package at: https://github.com/chloelulu/tdGBH. The analysis code and data are available at: https://github.com/chloelulu/tdGBH-paper.

Linking Prenatal Environmental Exposures to Lifetime Health with Epigenome-Wide Association Studies: State-of-the-Science Review and Future Recommendations

BACKGROUND

The prenatal environment influences lifetime health; epigenetic mechanisms likely predominate. In 2016, the first international consortium paper on cigarette smoking during pregnancy and offspring DNA methylation identified extensive, reproducible exposure signals. This finding raised expectations for epigenome-wide association studies (EWAS) of other exposures.

OBJECTIVE

We review the current state-of-the-science for DNA methylation associations across prenatal exposures in humans and provide future recommendations.

METHODS

We reviewed 134 prenatal environmental EWAS of DNA methylation in newborns, focusing on 51 epidemiological studies with meta-analysis or replication testing. Exposures spanned cigarette smoking, alcohol consumption, air pollution, dietary factors, psychosocial stress, metals, other chemicals, and other exogenous factors. Of the reproducible DNA methylation signatures, we examined implementation as exposure biomarkers.

RESULTS

Only 19 (14%) of these prenatal EWAS were conducted in cohorts of 1,000 or more individuals, reflecting the still early stage of the field. To date, the largest perinatal EWAS sample size was 6,685 participants. For comparison, the most recent genome-wide association study for birth weight included more than 300,000 individuals. Replication, at some level, was successful with exposures to cigarette smoking, folate, dietary glycemic index, particulate matter with aerodynamic diameter < 10 μ m and < 2.5 μ m , nitrogen dioxide, mercury, cadmium, arsenic, electronic waste, PFAS, and DDT. Reproducible effects of a more limited set of prenatal exposures (smoking, folate) enabled robust methylation biomarker creation.

DISCUSSION

Current evidence demonstrates the scientific premise for reproducible DNA methylation exposure signatures. Better powered EWAS could identify signatures across many exposures and enable comprehensive biomarker development. Whether methylation biomarkers of exposures themselves cause health effects remains unclear. We expect that larger EWAS with enhanced coverage of epigenome and exposome, along with improved single-cell technologies and evolving methods for integrative multi-omics analyses and causal inference, will expand mechanistic understanding of causal links between environmental exposures, the epigenome, and health outcomes throughout the life course. https://doi.org/10.1289/EHP12956.

Designing Studies for Epigenetic Biomarker Development in Autoimmune Rheumatic Diseases

In just a few years, the number of epigenetic studies in autoimmune rheumatic and inflammatory diseases has greatly increased. This is in part due to the need of identifying additional determinants to genetics to explain the pathogenesis and development of these disorders. In this regard, epigenetics provides potential mechanisms that determine gene function, are linked to environmental factors, and could explain a wide range of phenotypic variability among patients with these diseases. Despite the high interest and number of studies describing epigenetic alterations under these conditions and exploring their relationship to various clinical aspects, few of the proposed biomarkers have yet reached clinical practice. The potential of epigenetic markers is high, as these alterations link measurable features with a number of biological traits. In the present article, we present published studies in the field, discuss some frequent limitations in the existing research, and propose a number of considerations that should be taken into account by those starting new projects in the field, with an aim to generate biomarkers that could make it into the clinics.

Challenges in Analyzing Functional Epigenetic Data in Perspective of Adolescent Psychiatric Health

The formative period of adolescence plays a crucial role in the development of skills and abilities for adulthood. Adolescents who are affected by mental health conditions are at risk of suicide and social and academic impairments. Gene-environment complementary contributions to the molecular mechanisms involved in psychiatric disorders have emphasized the need to analyze epigenetic marks such as DNA methylation (DNAm) and non-coding RNAs. However, the large and diverse bioinformatic and statistical methods, referring to the confounders of the statistical models, application of multiple-testing adjustment methods, questions regarding the correlation of DNAm across tissues, and sex-dependent differences in results, have raised challenges regarding the interpretation of the results. Based on the example of generalized anxiety disorder (GAD) and depressive disorder (MDD), we shed light on the current knowledge and usage of methodological tools in analyzing epigenetics. Statistical robustness is an essential prerequisite for a better understanding and interpretation of epigenetic modifications and helps to find novel targets for personalized therapeutics in psychiatric diseases.

Leveraging Gene-Level Prediction as Informative Covariate in Hypothesis Weighting Improves Power for Rare Variant Association Studies

Gene-based rare variant association studies (RVASs) have low power due to the infrequency of rare variants and the large multiple testing burden. To correct for multiple testing, traditional false discovery rate (FDR) procedures which depend solely on P-values are often used. Recently, Independent Hypothesis Weighting (IHW) was developed to improve the detection power while maintaining FDR control by leveraging prior information for each hypothesis. Here, we present a framework to increase power of gene-based RVASs by incorporating prior information using IHW. We first build supervised machine learning models to assign each gene a prediction score that measures its disease risk, using the input of multiple biological features, fed with high-confidence risk genes and local background genes selected near GWAS significant loci as the training set. Then we use the prediction scores as covariates to prioritize RVAS results via IHW. We demonstrate the effectiveness of this framework through applications to RVASs in schizophrenia and autism spectrum disorder. We found sizeable improvements in the number of significant associations compared to traditional FDR approaches, and independent evidence supporting the relevance of the genes identified by our framework but not traditional FDR, demonstrating the potential of our framework to improve power of gene-based RVASs.

2dFDR: a new approach to confounder adjustment substantially increases detection power in omics association studies

One challenge facing omics association studies is the loss of statistical power when adjusting for confounders and multiple testing. The traditional statistical procedure involves fitting a confounder-adjusted regression model for each omics feature, followed by multiple testing correction. Here we show that the traditional procedure is not optimal and present a new approach, 2dFDR, a two-dimensional false discovery rate control procedure, for powerful confounder adjustment in multiple testing. Through extensive evaluation, we demonstrate that 2dFDR is more powerful than the traditional procedure, and in the presence of strong confounding and weak signals, the power improvement could be more than 100%.

Leveraging biological and statistical covariates improves the detection power in epigenome-wide association testing

BACKGROUND

Epigenome-wide association studies (EWAS), which seek the association between epigenetic marks and an outcome or exposure, involve multiple hypothesis testing. False discovery rate (FDR) control has been widely used for multiple testing correction. However, traditional FDR control methods do not use auxiliary covariates, and they could be less powerful if the covariates could inform the likelihood of the null hypothesis. Recently, many covariate-adaptive FDR control methods have been developed, but application of these methods to EWAS data has not yet been explored. It is not clear whether these methods can significantly improve detection power, and if so, which covariates are more relevant for EWAS data.

RESULTS

In this study, we evaluate the performance of five covariate-adaptive FDR control methods with EWAS-related covariates using simulated as well as real EWAS datasets. We develop an omnibus test to assess the informativeness of the covariates. We find that statistical covariates are generally more informative than biological covariates, and the covariates of methylation mean and variance are almost universally informative. In contrast, the informativeness of biological covariates depends on specific datasets. We show that the independent hypothesis weighting (IHW) and covariate adaptive multiple testing (CAMT) method are overall more powerful, especially for sparse signals, and could improve the detection power by a median of 25% and 68% on real datasets, compared to the ST procedure. We further validate the findings in various biological contexts.

CONCLUSIONS

Covariate-adaptive FDR control methods with informative covariates can significantly increase the detection power for EWAS. For sparse signals, IHW and CAMT are recommended.