Supplementary Material for: Integrative Network Analysis of Multi-Omics Data in the Link between Plasma Carotenoid Concentrations and Lipid Profile

Introduction: Carotenoids, which are a reliable biomarker of fruit and vegetable consumption, are positively associated with the lipid profile. Circulating carotenoid concentrations may interact with several omics profiles including genome, transcriptome, and epigenome. Few studies have used multi-omics approaches, and they rarely include environmental factors, such as diet. Objective: The objective of this observational study was to examine the potential role of multi-omics data in the interconnection between diet, represented by total carotenoids, and lipid profile using weighted gene correlation network analysis (WGCNA). Methods: Blood leukocyte DNA methylation levels of 472,245 CpG sites and whole blood gene expression levels of 18,160 transcripts were tested for associations with total carotenoid concentrations using regressions in 48 healthy subjects. WGCNA was used to identify co-omics modules and hub genes related to the lipid profile. Results: Among genes associated with total carotenoid concentrations, a total of 236 genes were identified at both DNA methylation and gene expression levels. Using WGCNA, six modules, consisting of groups of highly correlated genes represented by colors, were identified and linked to the lipid profile. Probes clustered in the turquoise and green modules correlated with plasma lipid concentrations. A total of 28 hub genes were identified. Conclusions: Genome-wide DNA methylation and gene expression levels were both associated with plasma total carotenoid concentrations. Several hub genes, mostly involved in lipid metabolism and inflammatory response with several genetic variants associated with plasma lipid concentrations, came out of the integrative analysis. This provides a comprehensive understanding of the interactive molecular system between carotenoids, omics, and plasma lipid profile.