A review of statistical strategies to integrate biomarkers of chemical exposure with biomarkers of effect applied in omic-scale environmental epidemiology

Humans are exposed to a growing list of synthetic chemicals, some of them becoming a major public health concern due to their capacity to impact multiple biological endpoints and contribute to a range of chronic diseases. The integration of endogenous (omic) biomarkers of effect in environmental health studies has been growing during the last decade, aiming to gain insight into potential mechanisms linking the exposures and the clinical conditions. The emergence of high-throughput omic platforms has raised a list of statistical challenges posed by the large dimension and complexity of data generated. Thus, the aim of the present study was to critically review the current state-of-the-science about statistical approaches used to integrate endogenous biomarkers in environmental-health studies linking chemical exposures with health outcomes. The present review specifically focused on internal exposure to environmental chemical pollutants, involving both persistent organic pollutants (POPs) and non-persistent pollutants like phthalates or bisphenols, and metals. We identified 42 eligible articles published since 2016, reporting 48 different statistical workflows, mostly focused on POPs and using metabolomic profiling in the intermediate layer. The outcomes were mainly binary and focused on metabolic disorders. A large diversity of statistical strategies were reported to integrate chemical mixtures and endogenous biomarkers to characterize their associations with health conditions. Multivariate regression models were the most predominant statistical method reported in the published workflows, however some studies applied latent based methods or multipollutant models to overcome the specific constraints of omic or exposure data. A minority of studies used formal mediation analysis to characterize the indirect effects mediated by the endogenous biomarkers. The principles of each specific statistical method and overall workflow set-up are summarized in the light of highlighting their applicability, strengths and weaknesses or interpretability to gain insight into the causal structures underlying the triad: exposure, effect-biomarker and outcome.