Epigenetics meets metabolomics: an epigenome-wide association study with blood serum metabolic traits

Genome-wide DNA methylation analysis of cannabis use disorder in a veteran cohort enriched for posttraumatic stress disorder

Cannabis use has been increasing over the past decade, not only in the general US population, but particularly among military veterans. With this rise in use has come a concomitant increase in cannabis use disorder (CUD) among veterans. Here, we performed an epigenome-wide association study for lifetime CUD in an Iraq/Afghanistan era veteran cohort enriched for posttraumatic stress disorder (PTSD) comprising 2,310 total subjects (1,109 non-Hispanic black and 1,201 non-Hispanic white). We also investigated CUD interactions with current PTSD status and examined potential indirect effects of DNA methylation (DNAm) on the relationship between CUD and psychiatric diagnoses. Four CpGs were associated with lifetime CUD, even after controlling for the effects of current smoking (AHRR cg05575921, LINC00299 cg23079012, VWA7 cg22112841, and FAM70A cg08760398). Importantly, cg05575921, a CpG strongly linked to smoking, remained associated with lifetime CUD even when restricting the analysis to veterans who reported never smoking cigarettes. Moreover, CUD interacted with current PTSD to affect cg05575921 and cg23079012 such that those with both CUD and PTSD displayed significantly lower DNAm compared to the other groups. Finally, we provide preliminary evidence that AHRR cg05575921 helps explain the association between CUD and any psychiatric diagnoses, specifically mood disorders.

Differential DNA methylation and metabolite profiling of Atlantic killifish (Fundulus heteroclitus) from the New Bedford Harbor Superfund site

Atlantic killifish (Fundulus heteroclitus) is a valuable model in evolutionary toxicology to study how the interactions between genetic and environmental factors serve the adaptive ability of organisms to resist chemical pollution. Killifish populations inhabiting environmental toxicant-contaminated New Bedford Harbor (NBH) show phenotypes tolerant to polychlorinated biphenyls (PCBs) and differences at the transcriptional and genomic levels. However, limited research has explored epigenetic alterations and metabolic effects in NBH killifish. To identify the involvement of epigenetic and metabolic regulation in the adaptive response of killifish, we investigated tissue- and sex-specific differences in global DNA methylation and metabolomic profiles of NBH killifish populations, compared to sensitive populations from a non-polluted site, Scorton Creek (SC). The results revealed that liver-specific global DNA hypomethylation and differential metabolites were evident in fish from NBH compared with those from SC. The sex-specific differences were not greater than the tissue-specific differences. We demonstrated liver-specific enriched metabolic pathways (e.g., amino acid metabolic pathways converged into the urea cycle and glutathione metabolism), suggesting possible crosstalk between differential metabolites and DNA hypomethylation in the livers of NBH killifish. Additional investigation of methylated gene regions is necessary to understand the functional role of DNA hypomethylation in the regulation of enzyme-encoding genes associated with metabolic processes and physiological changes in NBH populations.

Novel Time-dependent Multi-omics Integration in Sepsis-associated Liver Dysfunction

The recently developed technologies that allow the analysis of each single omics have provided an unbiased insight into ongoing disease processes. However, it remains challenging to specify the study design for the subsequent integration strategies that can associate sepsis pathophysiology and clinical outcomes. Here, we conducted a time-dependent multi-omics integration (TDMI) in a sepsis-associated liver dysfunction (SALD) model. We successfully deduced the relation of the Toll-like receptor 4 (TLR4) pathway with SALD. Although TLR4 is a critical factor in sepsis progression, it is not specified in single-omics analyses but only in the TDMI analysis. This finding indicates that the TDMI-based approach is more advantageous than single-omics analyses in terms of exploring the underlying pathophysiological mechanism of SALD. Furthermore, TDMI-based approach can be an ideal paradigm for insightful biological interpretations of multi-omics datasets that will potentially reveal novel insights into basic biology, health, and diseases, thus allowing the identification of promising candidates for therapeutic strategies.

Metabolome subtyping reveals multi-omics characteristics and biological heterogeneity in major psychiatric disorders

Growing evidence suggests that major psychiatric disorders (MPDs) share common etiologies and pathological processes. However, the diagnosis is currently based on descriptive symptoms, which ignores the underlying pathogenesis and hinders the development of clinical treatments. This highlights the urgency of characterizing molecular biomarkers and establishing objective diagnoses of MPDs. Here, we collected untargeted metabolomics, proteomics and DNA methylation data of 327 patients with MPDs, 131 individuals with genetic high risk and 146 healthy controls to explore the multi-omics characteristics of MPDs. First, differential metabolites (DMs) were identified and we classified MPD patients into 3 subtypes based on DMs. The subtypes showed distinct metabolomics, proteomics and DNA methylation signatures. Specifically, one subtype showed dysregulation of complement and coagulation proteins, while the DNA methylation showed abnormalities in chemical synapses and autophagy. Integrative analysis in metabolic pathways identified the important roles of the citrate cycle, sphingolipid metabolism and amino acid metabolism. Finally, we constructed prediction models based on the metabolites and proteomics that successfully captured the risks of MPD patients. Our study established molecular subtypes of MPDs and elucidated their biological heterogeneity through a multi-omics investigation. These results facilitate the understanding of pathological mechanisms and promote the diagnosis and prevention of MPDs.

Metabolomic biomarkers of habitual B vitamin intakes unveil novel differentially methylated positions in the human epigenome

BACKGROUND

B vitamins such as folate (B9), B6, and B12 are key in one carbon metabolism, which generates methyl donors for DNA methylation. Several studies have linked differential methylation to self-reported intakes of folate and B12, but these estimates can be imprecise, while metabolomic biomarkers can offer an objective assessment of dietary intakes. We explored blood metabolomic biomarkers of folate and vitamins B6 and B12, to carry out epigenome-wide analyses across up to three European cohorts. Associations between self-reported habitual daily B vitamin intakes and 756 metabolites (Metabolon Inc.) were assessed in serum samples from 1064 UK participants from the TwinsUK cohort. The identified B vitamin metabolomic biomarkers were then used in epigenome-wide association tests with fasting blood DNA methylation levels at 430,768 sites from the Infinium HumanMethylation450 BeadChip in blood samples from 2182 European participants from the TwinsUK and KORA cohorts. Candidate signals were explored for metabolite associations with gene expression levels in a subset of the TwinsUK sample (n = 297). Metabolomic biomarker epigenetic associations were also compared with epigenetic associations of self-reported habitual B vitamin intakes in samples from 2294 European participants.

RESULTS

Eighteen metabolites were associated with B vitamin intakes after correction for multiple testing (Bonferroni-adj. p < 0.05), of which 7 metabolites were available in both cohorts and tested for epigenome-wide association. Three metabolites - pipecolate (metabolomic biomarker of B6 and folate intakes), pyridoxate (marker of B6 and folate) and docosahexaenoate (DHA, marker of B6) - were associated with 10, 3 and 1 differentially methylated positions (DMPs), respectively. The strongest association was observed between DHA and DMP cg03440556 in the SCD gene (effect = 0.093 ± 0.016, p = 4.07E-09). Pyridoxate, a catabolic product of vitamin B6, was inversely associated with CpG methylation near the SLC1A5 gene promoter region (cg02711608 and cg22304262) and with SLC7A11 (cg06690548), but not with corresponding changes in gene expression levels. The self-reported intake of folate and vitamin B6 had consistent but non-significant associations with the epigenetic signals.

CONCLUSION

Metabolomic biomarkers are a valuable approach to investigate the effects of dietary B vitamin intake on the human epigenome.

DNA Methylation Near CPT1A and Changes in Triglyceride-rich Lipoproteins in Response to Weight-loss Diet Interventions

CONTEXT

Carnitine palmitoyltransferase-1A, encoded by the CPT1A gene, plays a key role in the oxidation of long-chain fatty acids in the mitochondria and may be important in triglyceride metabolism. Previous work has shown that high fat intake was negatively associated with CPT1A methylation and positively associated with CPT1A expression.

OBJECTIVE

We aim to investigate the association of DNA methylation (DNAm) at the CPT1A gene with reductions in triglycerides and triglyceride-rich lipoproteins (TRLs) in response to weight-loss diet interventions.

METHODS

The current study included 538 White participants, who were randomly assigned to 1 of 4 diets varying in macronutrient components. We defined the regional DNAm at CPT1A as the average methylation level over CpGs within 500 bp of the 3 triglyceride-related DNAm sites.

RESULTS

Dietary fat intake significantly modified the association between baseline DNAm at CPT1A and 2-year changes in total plasma triglycerides, independent of concurrent weight loss. Among participants assigned to a low-fat diet, a higher regional DNAm level at CPT1A was associated with a greater reduction in total plasma triglycerides at 2 years (P = .01), compared with those assigned to a high-fat diet (P = .64) (P interaction = .018). Further investigation on lipids and apolipoproteins in very low-density lipoprotein (VLDL) revealed similar interaction patterns for 2-year changes in VLDL-triglycerides, VLDL-cholesterol, and VLDL-apolipoprotein B (P interaction = .009, .002, and .016, respectively), but not for VLDL-apoC-III (P interaction = .36).

CONCLUSION

Participants with a higher regional DNAm level at CPT1A benefit more in long-term improvement in triglycerides, particularly in the TRLs and related apolipoproteins when consuming a low-fat weight-loss diet.

Identifying therapeutic biomarkers of zoledronic acid by metabolomics

Zoledronic acid (ZOL) is a potent antiresorptive agent that increases bone mineral density (BMD) and reduces fracture risk in postmenopausal osteoporosis (PMOP). The anti-osteoporotic effect of ZOL is determined by annual BMD measurement. In most cases, bone turnover markers function as early indicators of therapeutic response, but they fail to reflect long-term effects. We used untargeted metabolomics to characterize time-dependent metabolic shifts in response to ZOL and to screen potential therapeutic markers. In addition, bone marrow RNA-seq was performed to support plasma metabolic profiling. Sixty rats were assigned to sham-operated group (SHAM, n = 21) and ovariectomy group (OVX, n = 39) and received sham operation or bilateral ovariectomy, respectively. After modeling and verification, rats in the OVX group were further divided into normal saline group (NS, n = 15) and ZOL group (ZA, n = 18). Three doses of 100 μg/kg ZOL were administrated to the ZA group every 2 weeks to simulate 3-year ZOL therapy in PMOP. An equal volume of saline was administered to the SHAM and NS groups. Plasma samples were collected at five time points for metabolic profiling. At the end of the study, selected rats were euthanatized for bone marrow RNA-seq. A total number of 163 compound were identified as differential metabolites between the ZA and NS groups, including mevalonate, a critical molecule in target pathway of ZOL. In addition, prolyl hydroxyproline (PHP), leucyl hydroxyproline (LHP), 4-vinylphenol sulfate (4-VPS) were identified as differential metabolites throughout the study. Moreover, 4-VPS negatively correlated with increased vertebral BMD after ZOL administration as time-series analysis revealed. Bone marrow RNA-seq showed that the PI3K-AKT signaling pathway was significantly associated with ZOL-mediated changes in expression (adjusted-p = 0.018). In conclusion, mevalonate, PHP, LHP, and 4-VPS are candidate therapeutic markers of ZOL. The pharmacological effect of ZOL likely occurs through inhibition of the PI3K-AKT signaling pathway.

DNA methylation in peripheral blood leukocytes for the association with glucose metabolism and invasive breast cancer

BACKGROUND

Insulin resistance (IR) is a well-established factor for breast cancer (BC) risk in postmenopausal women, but the interrelated molecular pathways on the methylome are not explicitly described. We conducted a population-level epigenome-wide association (EWA) study for DNA methylation (DNAm) probes that are associated with IR and prospectively correlated with BC development, both overall and in BC subtypes among postmenopausal women.

METHODS

We used data from Women's Health Initiative (WHI) ancillary studies for our EWA analyses and evaluated the associations of site-specific DNAm across the genome with IR phenotypes by multiple regressions adjusting for age and leukocyte heterogeneities. For our analysis of the top 20 IR-CpGs with BC risk, we used the WHI and the Cancer Genomic Atlas (TCGA), using multiple Cox proportional hazards and logit regressions, respectively, accounting for age, diabetes, obesity, leukocyte heterogeneities, and tumor purity (for TCGA). We further conducted a Gene Set Enrichment Analysis.

RESULTS

We detected several EWA-CpGs in TXNIP, CPT1A, PHGDH, and ABCG1. In particular, cg19693031 in TXNIP was replicated in all IR phenotypes, measured by fasting levels of glucose, insulin, and homeostatic model assessment-IR. Of those replicated IR-genes, 3 genes (CPT1A, PHGDH, and ABCG1) were further correlated with BC risk; and 1 individual CpG (cg01676795 in POR) was commonly detected across the 2 cohorts.

CONCLUSIONS

Our study contributes to better understanding of the interconnected molecular pathways on the methylome between IR and BC carcinogenesis and suggests potential use of DNAm markers in the peripheral blood cells as preventive targets to detect an at-risk group for IR and BC in postmenopausal women.

Epigenome-wide DNA methylation in leukocytes and toenail metals: The normative aging study

BACKGROUND

Environmental metal exposures have been associated with multiple deleterious health endpoints. DNA methylation (DNAm) may provide insight into the mechanisms underlying these relationships. Toenail metals are non-invasive biomarkers, reflecting a medium-term time exposure window.

OBJECTIVES

This study examined variation in leukocyte DNAm and toenail arsenic (As), cadmium (Cd), lead (Pb), manganese (Mn), and mercury (Hg) among elderly men in the Normative Aging Study, a longitudinal cohort.

METHODS

We repeatedly collected samples of blood and toenail clippings. We measured DNAm in leukocytes with the Illumina HumanMethylation450 K BeadChip. We first performed median regression to evaluate the effects of each individual toenail metal on DNAm at three levels: individual cytosine-phosphate-guanine (CpG) sites, regions, and pathways. Then, we applied a Bayesian kernel machine regression (BKMR) to assess the joint and individual effects of metal mixtures on DNAm. Significant CpGs were identified using a multiple testing correction based on the independent degrees of freedom approach for correlated outcomes. The approach considers the effective degrees of freedom in the DNAm data using the principal components that explain >95% variation of the data.

RESULTS

We included 564 subjects (754 visits) between 1999 and 2013. The numbers of significantly differentially methylated CpG sites, regions, and pathways varied by metals. For example, we found six significant pathways for As, three for Cd, and one for Mn. The As-associated pathways were associated with cancer (e.g., skin cancer) and cardiovascular disease, whereas the Cd-associated pathways were related to lung cancer. Metal mixtures were also associated with 47 significant CpG sites, as well as pathways, mainly related to cancer and cardiovascular disease.

CONCLUSIONS

This study provides an approach to understanding the potential epigenetic mechanisms underlying observed relations between toenail metals and adverse health endpoints.

DNA methylation changes and inflammaging in aging-associated diseases

Aging as an inevitable phenomenon is associated with pervasive changes in physiological functions. There is a relationship between aging and the increase of several chronic diseases. Most age-related disorders are accompanied by an underlying chronic inflammatory state, as demonstrated by local infiltration of inflammatory cells and greater levels of proinflammatory cytokines in the bloodstream. Within inflammaging, many epigenetic events, especially DNA methylation, change. During the aging process, due to aberrations of DNA methylation, biological processes are disrupted, leading to the emergence or progression of a variety of human diseases, including cancer, neurodegenerative disorders, cardiovascular disease and diabetes. The focus of this review is on DNA methylation, which is involved in inflammaging-related activities, and how its dysregulation leads to human disorders.

DNA Methylation Profiles of Ovarian Granular Cells from Fluorosis Female Patients Suffering Reproductive Dysfunctions

Fluorosis often causes female reproductive dysfunction. A rapid turnover of DNA methylation is a pathological change in many human diseases, including female reproductive dysfunction. The role of DNA methylation in fluorosis was unknown and investigated in this experiment. Fifty fluorosis women patients were selected as High F group and forty-six healthy women were recruited as Control group were enrolled. In addition, ovarian granulosa cells were obtained from five women in High F group and five women in Control group. All ten women went through in vitro fertilization (IVF) process with DNA methylation sequencing. KGN cells (human granulosa cell line) were cultured with different concentrations of sodium fluoride (0-8 mM NaF) for 24 h for the in vitro study. The level of DNA methylation in blood samples was significantly higher in High F group than that in Control group. The level of serum estradiol (E₂) was significantly lower in women from High F group, while the levels of serum luteinizing hormone (LH) and follicle-stimulating hormone (FSH) in High F group were significantly higher than that in Control group. The methylation sequences of KGN cells relating to autophagy were significantly changed by NaF treatment dose-dependently. Based on these results, it is concluded that DNA methylation and autophagy may play a significant role in the pathophysiology of reproductive dysfunction caused by fluorosis.

Metabolic characteristics related to the hazardous effects of environmental arsenic on humans: A metabolomic review

Arsenic (As) is a toxic metalloid exist ubiquitously in environment. Epidemiological studies and laboratory animal studies have verified that As damages multiple organs or tissues in the body and is associated with a variety of diseases. Changes in metabolites usually indicate disturbances in metabolic pathways and specific metabolites are considered as biomarkers of diseases or drugs/toxins or environmental effects. Metabolomics is the quantitative measurement of the dynamic multi-parameter metabolic responses of biological systems due to pathophysiological or genetic changes. Current years, some metabolomic studies on the hazardous effect of environmental As on humans have been reported. In this paper, we first overviewed the metabolomics studies of environmental As exposure in humans since 2011, emphasizing on the data mining process of metabolic characteristics related to the hazardous effects of environmental As on humans. Then, the relationship between metabolic characteristics and the toxic mechanism of environmental As exposure in humans were discussed, and finally, the prospects of metabolomics studies on populations exposed to environmental As were put forward. Our paper may shed light on the study of mechanisms, prevention and individualized treatment of As poisoning.

Maternal hyperandrogenism is associated with a higher risk of type 2 diabetes mellitus and overweight in adolescent and adult female offspring: a long-term population-based follow-up study

PURPOSE

Adverse intrauterine environment may predispose offspring to cardio-metabolic dysfunction in later life. In this study, we aimed to investigate the effects of maternal hyperandrogenism (MH) on cardio-metabolic risk factors in female offspring in later life.

METHODS

This prospective population-based study included 211 female offspring with MH and 757 female offspring without MH (controls). Both groups were followed from baseline to the date of incidence of events, censoring, or end of the study period, whichever came first. Age scaled unadjusted and adjusted cox regression models were applied to assess the hazard ratios (HR) and 95% confidence intervals (CIs) for the association of MH with pre-diabetes (pre-DM), type 2 diabetes mellitus (T2DM), overweight and obesity in offspring of both groups. Statistical analysis was performed using the software package STATA; significance level was set at P < 0.05.

RESULTS

This study revealed a higher risk of T2DM (unadjusted HR 2.67, 95% CI 1.33-5.36) and overweight (unadjusted HR 1.41, 95% CI 1.06-1.88) in female offspring with MH, compared to controls. Results remained unchanged after adjustment for potential confounders including body mass index, education, physical activity, mother's age at delivery, birth weight, and childhood obesity. However, no significant difference was observed in the risk of pre-DM and obesity in females with MH, compared to controls in both unadjusted and adjusted models.

CONCLUSION

This pioneer study with a long-term follow-up demonstrated that MH increases the risk of developing T2DM and being overweight in female offspring in later life. Further long-term population-based studies are needed to confirm these findings.

Additive and Interactive Genetically Contextual Effects of HbA1c on cg19693031 Methylation in Type 2 Diabetes

Type 2 diabetes mellitus (T2D) has a complex genetic and environmental architecture that underlies its development and clinical presentation. Despite the identification of well over a hundred genetic variants and CpG sites that associate with T2D, a robust biosignature that could be used to prevent or forestall clinical disease has not been developed. Based on the premise that underlying genetic variation influences DNA methylation (DNAm) independently of or in combination with environmental exposures, we assessed the ability of local and distal gene x methylation (GxMeth) interactive effects to improve cg19693031 models for predicting T2D status in an African American cohort. Using genome-wide genetic data from 506 subjects, we identified a total of 1476 GxMeth terms associated with HbA1c values. The GxMeth SNPs map to biological pathways associated with the development and complications of T2D, with genetically contextual differences in methylation observed only in diabetic subjects for two GxMeth SNPs (rs2390998 AG vs. GG, p = 4.63 × 10−11, Δβ = 13%, effect size = 0.16 [95% CI = 0.05, 0.32]; rs1074390 AA vs. GG, p = 3.93 × 10−4, Δβ = 9%, effect size = 0.38 [95% CI = 0.12, 0.56]. Using a repeated stratified k-fold cross-validation approach, a series of balanced random forest classifiers with random under-sampling were built to evaluate the addition of GxMeth terms to cg19693031 models to discriminate between normoglycemic controls versus T2D subjects. The results were compared to those obtained from models incorporating only the covariates (age, sex and BMI) and the addition of cg19693031. We found a post-pruned classifier incorporating 10 GxMeth SNPs and cg19693031 adjusted for covariates predicted the T2D status, with the AUC, sensitivity, specificity and precision of the positive target class being 0.76, 0.81, 0.70 and 0.63, respectively. Comparatively, the AUC, sensitivity, specificity and precision using the covariates and cg19693031 were only 0.71, 0.74, 0.67 and 0.59, respectively. Collectively, we demonstrate correcting for genetic confounding of cg19693031 improves its ability to detect type 2 diabetes. We conclude that an integrated genetic–epigenetic approach could inform personalized medicine programming for more effective prevention and treatment of T2D.

DNA methylation level of the gene encoding thioredoxin-interacting protein in peripheral blood cells is associated with metabolic syndrome in the Japanese general population

Metabolic syndrome (MetS) is cluster of metabolic diseases, including abdominal obesity, hyperglycemia, high blood pressure, and dyslipidemia, that directly escalate the risk of type 2 diabetes, heart disease, and stroke. Thioredoxin-interacting protein (TXNIP) is a binding protein for thioredoxin, a molecule that is a key inhibitor of cellular oxidation, and thus regulates the cellular redox state. Epigenetic alteration of the TXNIP-encoding locus has been associated with components of MetS. In the present study, we sought to determine whether the level of TXNIP methylation in blood is associated with MetS in the general Japanese population. DNA was extracted from the peripheral blood cells of 37 subjects with and 392 subjects without MetS. The level of TXNIP methylation at cg19693031 was assessed by the bisulfite-pyrosequencing method. We observed that TXNIP methylation levels were lower in MetS subjects (median 74.9%, range 71.7-78.4%) than in non-MetS subjects (median 77.7%, range 74.4-80.5%; p = 0.0024). Calculation of the confounding factor-adjusted odds ratios (ORs) and 95% confidence intervals (CIs) for hypomethylation revealed that subjects with MetS exhibited significantly higher ORs for hypomethylation than did those without MetS (OR, 2.92; 95% CI, 1.33-6.62; p = 0.009). Our findings indicated that lower levels of TXNIP methylation are associated with MetS in the general Japanese population. Altered levels of DNA methylation in TXNIP at cg19693031 might play an important role in the pathogenesis of MetS.

Epigenetic scores for the circulating proteome as tools for disease prediction

Protein biomarkers have been identified across many age-related morbidities. However, characterising epigenetic influences could further inform disease predictions. Here, we leverage epigenome-wide data to study links between the DNA methylation (DNAm) signatures of the circulating proteome and incident diseases. Using data from four cohorts, we trained and tested epigenetic scores (EpiScores) for 953 plasma proteins, identifying 109 scores that explained between 1% and 58% of the variance in protein levels after adjusting for known protein quantitative trait loci (pQTL) genetic effects. By projecting these EpiScores into an independent sample (Generation Scotland; n = 9537) and relating them to incident morbidities over a follow-up of 14 years, we uncovered 137 EpiScore-disease associations. These associations were largely independent of immune cell proportions, common lifestyle and health factors, and biological aging. Notably, we found that our diabetes-associated EpiScores highlighted previous top biomarker associations from proteome-wide assessments of diabetes. These EpiScores for protein levels can therefore be a valuable resource for disease prediction and risk stratification.

Epigenome-wide association study of serum urate reveals insights into urate co-regulation and the SLC2A9 locus

Elevated serum urate levels, a complex trait and major risk factor for incident gout, are correlated with cardiometabolic traits via incompletely understood mechanisms. DNA methylation in whole blood captures genetic and environmental influences and is assessed in transethnic meta-analysis of epigenome-wide association studies (EWAS) of serum urate (discovery, n = 12,474, replication, n = 5522). The 100 replicated, epigenome-wide significant (p < 1.1E-7) CpGs explain 11.6% of the serum urate variance. At SLC2A9, the serum urate locus with the largest effect in genome-wide association studies (GWAS), five CpGs are associated with SLC2A9 gene expression. Four CpGs at SLC2A9 have significant causal effects on serum urate levels and/or gout, and two of these partly mediate the effects of urate-associated GWAS variants. In other genes, including SLC7A11 and PHGDH, 17 urate-associated CpGs are associated with conditions defining metabolic syndrome, suggesting that these CpGs may represent a blood DNA methylation signature of cardiometabolic risk factors. This study demonstrates that EWAS can provide new insights into GWAS loci and the correlation of serum urate with other complex traits.

Blood DNA methylation markers associated with type 2 diabetes, fasting glucose, and HbA1c levels: An epigenome-wide association study in 316 adult twin pairs

DNA methylation plays an important role in the development and etiology of type 2 diabetes; however, few epigenomic studies have been conducted on twins. Herein, a two-stage study was performed to explore the associations between DNA methylation and type 2 diabetes, fasting plasma glucose, and HbA1c. DNA methylation in 316 twin pairs from the Chinese National Twin Registry (CNTR) was measured using Illumina Infinium BeadChips. In the discovery sample, the results revealed that 63 CpG sites and 6 CpG sites were significantly associated with fasting plasma glucose and HbA1c, respectively. In the replication sample, cg19690313 in TXNIP was associated with both fasting plasma glucose (P = 1.23 × 10^-17, FDR < 0.001) and HbA1c (P = 2.29 × 10^-18, FDR < 0.001). Furthermore, cg04816311, cg08309687, and cg09249494 may provide new insight in the metabolic mechanism of HbA1c. Our study provides solid evidence that cg19690313 on TXNIP correlates with HbA1c and fasting plasma glucose levels.

Overview of Metabolomic Analysis and the Integration with Multi-Omics for Economic Traits in Cattle

Metabolomics has been applied to measure the dynamic metabolic responses, to understand the systematic biological networks, to reveal the potential genetic architecture, etc., for human diseases and livestock traits. For example, the current published results include the detected relevant candidate metabolites, identified metabolic pathways, potential systematic networks, etc., for different cattle traits that can be applied for further metabolomic and integrated omics studies. Therefore, summarizing the applications of metabolomics for economic traits is required in cattle. We here provide a comprehensive review about metabolomic analysis and its integration with other omics in five aspects: (1) characterization of the metabolomic profile of cattle; (2) metabolomic applications in cattle; (3) integrated metabolomic analysis with other omics; (4) methods and tools in metabolomic analysis; and (5) further potentialities. The review aims to investigate the existing metabolomic studies by highlighting the results in cattle, integrated with other omics studies, to understand the metabolic mechanisms underlying the economic traits and to provide useful information for further research and practical breeding programs in cattle.

Using a Two-Sample Mendelian Randomization Method in Assessing the Causal Relationships Between Human Blood Metabolites and Heart Failure

Background: Heart failure (HF) is the main cause of morbidity and mortality worldwide, and metabolic dysfunction is an important factor related to HF pathogenesis and development. However, the causal effect of blood metabolites on HF remains unclear.

Objectives: Our chief aim is to investigate the causal relationships between human blood metabolites and HF risk.

Methods: We used an unbiased two-sample Mendelian randomization (MR) approach to assess the causal relationships between 486 human blood metabolites and HF risk. Exposure information was obtained from Sample 1, which is the largest metabolome-based genome-wide association study (mGWAS) data containing 7,824 Europeans. Outcome information was obtained from Sample 2, which is based on the results of a large-scale GWAS meta-analysis of HF and contains 47,309 cases and 930,014 controls of Europeans. The inverse variance weighted (IVW) model was used as the primary two-sample MR analysis method and followed the sensitivity analyses, including heterogeneity test, horizontal pleiotropy test, and leave-one-out analysis.

Results: We observed that 11 known metabolites were potentially related to the risk of HF after using the IVW method (P < 0.05). After adding another four MR models and performing sensitivity analyses, we found a 1-SD increase in the xenobiotics 4-vinylphenol sulfate was associated with ~22% higher risk of HF (OR [95%CI], 1.22 [1.07–1.38]).

Conclusions: We revealed that the 4-vinylphenol sulfate may nominally increase the risk of HF by 22% after using a two-sample MR approach. Our findings may provide novel insights into the pathogenesis underlying HF and novel strategies for HF prevention.

Connecting the epigenome, metabolome and proteome for a deeper understanding of disease

Epigenome-wide association studies (EWAS) identify genes that are dysregulated by the studied clinical endpoints, thereby indicating potential new diagnostic biomarkers, drug targets and therapy options. Combining EWAS with deep molecular phenotyping, such as approaches enabled by metabolomics and proteomics, allows further probing of the underlying disease-associated pathways. For instance, methylation of the TXNIP gene is associated robustly with prevalent type 2 diabetes and further with metabolites that are short-term markers of glycaemic control. These associations reflect TXNIP's function as a glucose uptake regulator by interaction with the major glucose transporter GLUT1 and suggest that TXNIP methylation can be used as a read-out for the organism's exposure to glucose stress. Another case is the association between DNA methylation of the AHRR and F2RL3 genes with smoking and a protein that is involved in the reprogramming of the bronchial epithelium. These examples show that associations between DNA methylation and intermediate molecular traits can open new windows into how the body copes with physiological challenges. This knowledge, if carefully interpreted, may indicate novel therapy options and, together with monitoring of the methylation state of specific methylation sites, may in the future allow the early diagnosis of impending disease. It is essential for medical practitioners to recognize the potential that this field holds in translating basic research findings to clinical practice. In this review, we present recent advances in the field of EWAS with metabolomics and proteomics and discuss both the potential and the challenges of translating epigenetic associations, with deep molecular phenotypes, to biomedical applications.

Low THRB (thyroid hormone receptor beta) Promoter Methylation Levels in Peripheral Blood Leukocytes Induced By Systematic Inflammation Are Involved in Low Thyroid Hormone Function in Metabolic Syndrome

[Figure: see text].

Characterizing the properties of bisulfite sequencing data: maximizing power and sensitivity to identify between-group differences in DNA methylation

BACKGROUND

The combination of sodium bisulfite treatment with highly-parallel sequencing is a common method for quantifying DNA methylation across the genome. The power to detect between-group differences in DNA methylation using bisulfite-sequencing approaches is influenced by both experimental (e.g. read depth, missing data and sample size) and biological (e.g. mean level of DNA methylation and difference between groups) parameters. There is, however, no consensus about the optimal thresholds for filtering bisulfite sequencing data with implications for the reproducibility of findings in epigenetic epidemiology.

RESULTS

We used a large reduced representation bisulfite sequencing (RRBS) dataset to assess the distribution of read depth across DNA methylation sites and the extent of missing data. To investigate how various study variables influence power to identify DNA methylation differences between groups, we developed a framework for simulating bisulfite sequencing data. As expected, sequencing read depth, group size, and the magnitude of DNA methylation difference between groups all impacted upon statistical power. The influence on power was not dependent on one specific parameter, but reflected the combination of study-specific variables. As a resource to the community, we have developed a tool, POWEREDBiSeq, which utilizes our simulation framework to predict study-specific power for the identification of DNAm differences between groups, taking into account user-defined read depth filtering parameters and the minimum sample size per group.

CONCLUSIONS

Our data-driven approach highlights the importance of filtering bisulfite-sequencing data by minimum read depth and illustrates how the choice of threshold is influenced by the specific study design and the expected differences between groups being compared. The POWEREDBiSeq tool, which can be applied to different types of bisulfite sequencing data (e.g. RRBS, whole genome bisulfite sequencing (WGBS), targeted bisulfite sequencing and amplicon-based bisulfite sequencing), can help users identify the level of data filtering needed to optimize power and aims to improve the reproducibility of bisulfite sequencing studies.

Epigenome-wide association meta-analysis of DNA methylation with coffee and tea consumption

Coffee and tea are extensively consumed beverages worldwide which have received considerable attention regarding health. Intake of these beverages is consistently linked to, among others, reduced risk of diabetes and liver diseases; however, the mechanisms of action remain elusive. Epigenetics is suggested as a mechanism mediating the effects of dietary and lifestyle factors on disease onset. Here we report the results from epigenome-wide association studies (EWAS) on coffee and tea consumption in 15,789 participants of European and African-American ancestries from 15 cohorts. EWAS meta-analysis of coffee consumption reveals 11 CpGs surpassing the epigenome-wide significance threshold (P-value <1.1×10-7), which annotated to the AHRR, F2RL3, FLJ43663, HDAC4, GFI1 and PHGDH genes. Among them, cg14476101 is significantly associated with expression of the PHGDH and risk of fatty liver disease. Knockdown of PHGDH expression in liver cells shows a correlation with expression levels of genes associated with circulating lipids, suggesting a role of PHGDH in hepatic-lipid metabolism. EWAS meta-analysis on tea consumption reveals no significant association, only two CpGs annotated to CACNA1A and PRDM16 genes show suggestive association (P-value <5.0×10-6). These findings indicate that coffee-associated changes in DNA methylation levels may explain the mechanism of action of coffee consumption in conferring risk of diseases.

The remodel of the "central dogma": a metabolomics interaction perspective

BACKGROUND

In 1957, Francis Crick drew a linear diagram on a blackboard. This diagram is often called the "central dogma." Subsequently, the relationships between different steps of the "central dogma" have been shown to be considerably complex, mostly because of the emerging world of small molecules. It is noteworthy that metabolites can be generated from the diet through gut microbiome metabolism, serve as substrates for epigenetic modifications, destabilize DNA quadruplexes, and follow Lamarckian inheritance. Small molecules were once considered the missing link in the "central dogma"; however, recently they have acquired a central role, and their general perception as downstream products has become reductionist. Metabolomics is a large-scale analysis of metabolites, and this emerging field has been shown to be the closest omics associated with the phenotype and concomitantly, the basis for all omics.

AIM OF REVIEW

Herein, we propose a broad updated perspective for the flux of information diagram centered in metabolomics, including the influence of other factors, such as epigenomics, diet, nutrition, and the gut- microbiome.

KEY SCIENTIFIC CONCEPTS OF REVIEW

Metabolites are the beginning and the end of the flux of information.

Global biochemical analysis of plasma, serum and whole blood collected using various anticoagulant additives

Introduction

Analysis of blood for the evaluation of clinically relevant biomarkers requires precise collection and sample handling by phlebotomists and laboratory staff. An important consideration for the clinical application of metabolomics are the different anticoagulants utilized for sample collection. Most studies that have characterized differences in metabolite levels in various blood collection tubes have focused on single analytes. We define analyte levels on a global metabolomics platform following blood sampling using five different, but commonly used, clinical laboratory blood collection tubes (i.e., plasma anticoagulated with either EDTA, lithium heparin or sodium citrate, along with no additive (serum), and EDTA anticoagulated whole blood).

Methods

Using an untargeted metabolomics platform we analyzed five sample types after all had been collected and stored at -80°C. The biochemical composition was determined and differences between the samples established using matched-pair t-tests.

Results

We identified 1,117 biochemicals across all samples and detected a mean of 1,036 in the sample groups. Compared to the levels of metabolites in EDTA plasma, the number of biochemicals present at statistically significant different levels (p<0.05) ranged from 452 (serum) to 917 (whole blood). Several metabolites linked to screening assays for rare diseases including acylcarnitines, bilirubin and heme metabolites, nucleosides, and redox balance metabolites varied significantly across the sample collection types.

Conclusions

Our study highlights the widespread effects and importance of using consistent additives for assessing small molecule levels in clinical metabolomics. The biochemistry that occurs during the blood collection process creates a reproducible signal that can identify specimens collected with different anticoagulants in metabolomic studies.

Impact statement

In this manuscript, normal/healthy donors had peripheral blood collected using multiple anticoagulants as well as serum during a fasted blood draw. Global metabolomics is a new technology being utilized to draw clinical conclusions and we interrogated the effects of different anticoagulants on the levels of biochemicals from each of the donors. Characterizing the effects of the anticoagulants on biochemical levels will help researchers leverage the information using global metabolomics in order to make conclusions regarding important disease biomarkers.

DNA methylation of blood cells is associated with prevalent type 2 diabetes in a meta-analysis of four European cohorts

BACKGROUND

Type 2 diabetes (T2D) is a heterogeneous disease with well-known genetic and environmental risk factors contributing to its prevalence. Epigenetic mechanisms related to changes in DNA methylation (DNAm), may also contribute to T2D risk, but larger studies are required to discover novel markers, and to confirm existing ones.

RESULTS

We performed a large meta-analysis of individual epigenome-wide association studies (EWAS) of prevalent T2D conducted in four European studies using peripheral blood DNAm. Analysis of differentially methylated regions (DMR) was also undertaken, based on the meta-analysis results. We found three novel CpGs associated with prevalent T2D in Europeans at cg00144180 (HDAC4), cg16765088 (near SYNM) and cg24704287 (near MIR23A) and confirmed three CpGs previously identified (mapping to TXNIP, ABCG1 and CPT1A). We also identified 77 T2D associated DMRs, most of them hypomethylated in T2D cases versus controls. In adjusted regressions among diabetic-free participants in ALSPAC, we found that all six CpGs identified in the meta-EWAS were associated with white cell-types. We estimated that these six CpGs captured 11% of the variation in T2D, which was similar to the variation explained by the model including only the common risk factors of BMI, sex, age and smoking (R2 = 10.6%).

CONCLUSIONS

This study identifies novel loci associated with T2D in Europeans. We also demonstrate associations of the same loci with other traits. Future studies should investigate if our findings are generalizable in non-European populations, and potential roles of these epigenetic markers in T2D etiology or in determining long term consequences of T2D.

DNA methylation and lipid metabolism: an EWAS of 226 metabolic measures

BACKGROUND

The discovery of robust and trans-ethnically replicated DNA methylation markers of metabolic phenotypes, has hinted at a potential role of epigenetic mechanisms in lipid metabolism. However, DNA methylation and the lipid compositions and lipid concentrations of lipoprotein sizes have been scarcely studied. Here, we present an epigenome-wide association study (EWAS) (N = 5414 total) of mostly lipid-related metabolic measures, including a fine profiling of lipoproteins. As lipoproteins are the main players in the different stages of lipid metabolism, examination of epigenetic markers of detailed lipoprotein features might improve the diagnosis, prognosis, and treatment of metabolic disturbances.

RESULTS

We conducted an EWAS of leukocyte DNA methylation and 226 metabolic measurements determined by nuclear magnetic resonance spectroscopy in the population-based KORA F4 study (N = 1662) and replicated the results in the LOLIPOP, NFBC1966, and YFS cohorts (N = 3752). Follow-up analyses in the discovery cohort included investigations into gene transcripts, metabolic-measure ratios for pathway analysis, and disease endpoints. We identified 161 associations (p value < 4.7 × 10-10), covering 16 CpG sites at 11 loci and 57 metabolic measures. Identified metabolic measures were primarily medium and small lipoproteins, and fatty acids. For apolipoprotein B-containing lipoproteins, the associations mainly involved triglyceride composition and concentrations of cholesterol esters, triglycerides, free cholesterol, and phospholipids. All associations for HDL lipoproteins involved triglyceride measures only. Associated metabolic measure ratios, proxies of enzymatic activity, highlight amino acid, glucose, and lipid pathways as being potentially epigenetically implicated. Five CpG sites in four genes were associated with differential expression of transcripts in blood or adipose tissue. CpG sites in ABCG1 and PHGDH showed associations with metabolic measures, gene transcription, and metabolic measure ratios and were additionally linked to obesity or previous myocardial infarction, extending previously reported observations.

CONCLUSION

Our study provides evidence of a link between DNA methylation and the lipid compositions and lipid concentrations of different lipoprotein size subclasses, thus offering in-depth insights into well-known associations of DNA methylation with total serum lipids. The results support detailed profiling of lipid metabolism to improve the molecular understanding of dyslipidemia and related disease mechanisms.

Multi-omics integration in biomedical research - A metabolomics-centric review

Recent advances in high-throughput technologies have enabled the profiling of multiple layers of a biological system, including DNA sequence data (genomics), RNA expression levels (transcriptomics), and metabolite levels (metabolomics). This has led to the generation of vast amounts of biological data that can be integrated in so-called multi-omics studies to examine the complex molecular underpinnings of health and disease. Integrative analysis of such datasets is not straightforward and is particularly complicated by the high dimensionality and heterogeneity of the data and by the lack of universal analysis protocols. Previous reviews have discussed various strategies to address the challenges of data integration, elaborating on specific aspects, such as network inference or feature selection techniques. Thereby, the main focus has been on the integration of two omics layers in their relation to a phenotype of interest. In this review we provide an overview over a typical multi-omics workflow, focusing on integration methods that have the potential to combine metabolomics data with two or more omics. We discuss multiple integration concepts including data-driven, knowledge-based, simultaneous and step-wise approaches. We highlight the application of these methods in recent multi-omics studies, including large-scale integration efforts aiming at a global depiction of the complex relationships within and between different biological layers without focusing on a particular phenotype.

Integrative Analysis of Gene-Specific DNA Methylation and Untargeted Metabolomics Data from the ELEMENT Cohort

Epigenetic modifications, such as DNA methylation, influence gene expression and cardiometabolic phenotypes that are manifest in developmental periods in later life, including adolescence. Untargeted metabolomics analysis provide a comprehensive snapshot of physiological processes and metabolism and have been related to DNA methylation in adults, offering insights into the regulatory networks that influence cellular processes. We analyzed the cross-sectional correlation of blood leukocyte DNA methylation with 3758 serum metabolite features (574 of which are identifiable) in 238 children (ages 8-14 years) from the Early Life Exposures in Mexico to Environmental Toxicants (ELEMENT) study. Associations between these features and percent DNA methylation in adolescent blood leukocytes at LINE-1 repetitive elements and genes that regulate early life growth (IGF2, H19, HSD11B2) were assessed by mixed effects models, adjusting for sex, age, and puberty status. After false discovery rate correction (FDR q < 0.05), 76 metabolites were significantly associated with LINE-1 DNA methylation, 27 with HSD11B2, 103 with H19, and 4 with IGF2. The ten identifiable metabolites included dicarboxylic fatty acids (five associated with LINE-1 or H19 methylation at q < 0.05) and 1-octadecanoyl-rac-glycerol (q < 0.0001 for association with H19 and q = 0.04 for association with LINE-1). We then assessed the association between these ten known metabolites and adiposity 3 years later. Two metabolites, dicarboxylic fatty acid 17:3 and 5-oxo-7-octenoic acid, were inversely associated with measures of adiposity (P < .05) assessed approximately 3 years later in adolescence. In stratified analyses, sex-specific and puberty-stage specific (Tanner stage = 2 to 5 vs Tanner stage = 1) associations were observed. Most notably, hundreds of statistically significant associations were observed between H19 and LINE-1 DNA methylation and metabolites among children who had initiated puberty. Understanding relationships between subclinical molecular biomarkers (DNA methylation and metabolites) may increase our understanding of genes and biological pathways contributing to metabolic changes that underlie the development of adiposity during adolescence.

Integrative Metabolomics to Identify Molecular Signatures of Responses to Vaccines and Infections

Approaches to the identification of metabolites have progressed from early biochemical pathway evaluation to modern high-dimensional metabolomics, a powerful tool to identify and characterize biomarkers of health and disease. In addition to its relevance to classic metabolic diseases, metabolomics has been key to the emergence of immunometabolism, an important area of study, as leukocytes generate and are impacted by key metabolites important to innate and adaptive immunity. Herein, we discuss the metabolomic signatures and pathways perturbed by the activation of the human immune system during infection and vaccination. For example, infection induces changes in lipid (e.g., free fatty acids, sphingolipids, and lysophosphatidylcholines) and amino acid pathways (e.g., tryptophan, serine, and threonine), while vaccination can trigger changes in carbohydrate and bile acid pathways. Amino acid, carbohydrate, lipid, and nucleotide metabolism is relevant to immunity and is perturbed by both infections and vaccinations. Metabolomics holds substantial promise to provide fresh insight into the molecular mechanisms underlying the host immune response. Its integration with other systems biology platforms will enhance studies of human health and disease.

An epigenome-wide association study of metabolic syndrome and its components

The role of metabolic syndrome (MetS) as a preceding metabolic state for type 2 diabetes and cardiovascular disease is widely recognised. To accumulate knowledge of the pathological mechanisms behind the condition at the methylation level, we conducted an epigenome-wide association study (EWAS) of MetS and its components, testing 1187 individuals of European ancestry for approximately 470 000 methylation sites throughout the genome. Methylation site cg19693031 in gene TXNIP —previously associated with type 2 diabetes, glucose and lipid metabolism, associated with fasting glucose level (P = 1.80 × 10⁻⁸). Cg06500161 in gene ABCG1 associated both with serum triglycerides (P = 5.36 × 10⁻⁹) and waist circumference (P = 5.21 × 10⁻⁹). The previously identified type 2 diabetes–associated locus cg08309687 in chromosome 21 associated with waist circumference for the first time (P = 2.24 × 10⁻⁷). Furthermore, a novel HDL association with cg17901584 in chromosome 1 was identified (P = 7.81 × 10⁻⁸). Our study supports previous genetic studies of MetS, finding that lipid metabolism plays a key role in pathology of the syndrome. We provide evidence regarding a close interplay with glucose metabolism. Finally, we suggest that in attempts to identify methylation loci linking separate MetS components, cg19693031 appears to represent a strong candidate.

Epigenetic loci for blood pressure are associated with hypertensive target organ damage in older African Americans from the genetic epidemiology network of Arteriopathy (GENOA) study

BACKGROUND

Hypertension is a major modifiable risk factor for arteriosclerosis that can lead to target organ damage (TOD) of heart, kidneys, and peripheral arteries. A recent epigenome-wide association study for blood pressure (BP) identified 13 CpG sites, but it is not known whether DNA methylation at these sites is also associated with TOD.

METHODS

In 1218 African Americans from the Genetic Epidemiology Network of Arteriopathy (GENOA) study, a cohort of hypertensive sibships, we evaluated the associations between methylation at these 13 CpG sites measured in peripheral blood leukocytes and five TOD traits assessed approximately 5 years later.

RESULTS

Ten significant associations were found after adjustment for age, sex, blood cell counts, time difference between CpG and TOD measurement, and 10 genetic principal components (FDR q < 0.1): two with estimated glomerular filtration rate (eGFR, cg06690548, cg10601624), six with urinary albumin-to-creatinine ratio (UACR, cg16246545, cg14476101, cg19693031, cg06690548, cg00574958, cg22304262), and two with left ventricular mass indexed to height (LVMI, cg19693031, cg00574958). All associations with eGFR and four associations with UACR remained significant after further adjustment for body mass index (BMI), smoking status, and diabetes. We also found significant interactions between cg06690548 and BMI on UACR, and between 3 CpG sites (cg19693031, cg14476101, and cg06690548) and diabetes on UACR (FDR q < 0.1). Mediation analysis showed that 4.7% to 38.1% of the relationship between two CpG sites (cg19693031 and cg00574958) and two TOD measures (UACR and LVMI) was mediated by blood pressure (Bonferroni-corrected P < 0.05). Mendelian randomization analysis suggests that methylation at two sites (cg16246545 and cg14476101) in PHGDH may causally influence UACR.

CONCLUSIONS

In conclusion, we found compelling evidence for associations between arteriosclerotic traits of kidney and heart and previously identified blood pressure-associated DNA methylation sites. This study may lend insight into the role of DNA methylation in pathological mechanisms underlying target organ damage from hypertension.

Methylome-wide association study of central adiposity implicates genes involved in immune and endocrine systems

Aim: We conducted a methylome-wide association study to examine associations between DNA methylation in whole blood and central adiposity and body fat distribution, measured as waist circumference, waist-to-hip ratio and waist-to-height ratio adjusted for body mass index, in 2684 African-American adults in the Atherosclerosis Risk in Communities study. Materials & methods: We validated significantly associated cytosine-phosphate-guanine methylation sites (CpGs) among adults using the Women's Health Initiative and Framingham Heart Study participants (combined n = 5743) and generalized associations in adolescents from The Raine Study (n = 820). Results & conclusion: We identified 11 CpGs that were robustly associated with one or more central adiposity trait in adults and two in adolescents, including CpG site associations near TXNIP, ADCY7, SREBF1 and RAP1GAP2 that had not previously been associated with obesity-related traits.

A multi-omic analysis of birthweight in newborn cord blood reveals new underlying mechanisms related to cholesterol metabolism

BACKGROUND

Birthweight reflects in utero exposures and later health evolution. Despite existing studies employing high-dimensional molecular measurements, the understanding of underlying mechanisms of birthweight remains limited.

METHODS

To investigate the systems biology of birthweight, we cross-sectionally integrated the methylome, the transcriptome, the metabolome and a set of inflammatory proteins measured in cord blood samples, collected from four birth-cohorts (n = 489). We focused on two sets of 68 metabolites and 903 CpGs previously related to birthweight and investigated the correlation structures existing between these two sets and all other omic features via bipartite Pearson correlations.

RESULTS

This dataset revealed that the set of metabolome and methylome signatures of birthweight have seven signals in common, including three metabolites [PC(34:2), plasmalogen PC(36:4)/PC(O-36:5), and a compound with m/z of 781.0545], two CpGs (on the DHCR24 and SC4MOL gene), and two proteins (periostin and CCL22). CCL22, a macrophage-derived chemokine has not been previously identified in relation to birthweight. Since the results of the omics integration indicated the central role of cholesterol metabolism, we explored the association of cholesterol levels in cord blood with birthweight in the ENVIRONAGE cohort (n = 1097), finding that higher birthweight was associated with increased high-density lipoprotein cholesterol and that high-density lipoprotein cholesterol was lower in small versus large for gestational age newborns.

CONCLUSIONS

Our data suggests that an integration of different omic-layers in addition to single omics studies is a useful approach to generate new hypotheses regarding biological mechanisms. CCL22 and cholesterol metabolism in cord blood play a mechanistic role in birthweight.

NRF2 loss recapitulates heritable impacts of paternal cigarette smoke exposure

Paternal cigarette smoke (CS) exposure is associated with increased risk of behavioral disorders and cancer in offspring, but the mechanism has not been identified. Here we use mouse models to investigate mechanisms and impacts of paternal CS exposure. We demonstrate that CS exposure induces sperm DNAme changes that are partially corrected within 28 days of removal from CS exposure. Additionally, paternal smoking is associated with changes in prefrontal cortex DNAme and gene expression patterns in offspring. Remarkably, the epigenetic and transcriptional effects of CS exposure that we observed in wild type mice are partially recapitulated in Nrf2-/- mice and their offspring, independent of smoking status. Nrf2 is a central regulator of antioxidant gene transcription, and mice lacking Nrf2 consequently display elevated oxidative stress, suggesting that oxidative stress may underlie CS-induced heritable epigenetic changes. Importantly, paternal sperm DNAme changes do not overlap with DNAme changes measured in offspring prefrontal cortex, indicating that the observed DNAme changes in sperm are not directly inherited. Additionally, the changes in sperm DNAme associated with CS exposure were not observed in sperm of unexposed offspring, suggesting the effects are likely not maintained across multiple generations.

Notch1 signaling mediated dysfunction of bone marrow mesenchymal stem cells derived from cyanotic congenital heart disease

Bone marrow mesenchymal stem cells (BMSCs) derived from cyanotic congenital heart disease (CCHD) exhibit deficient multi-lineage differentiation potential due to the abnormal accumulation of D-galactose. However, the underlying mechanisms have not yet been explored. Here, the multi-lineage differentiation potential of the BMSCs from CCHD and non-CCHD (NCHD) patients were assessed. BMSCs from CCHD patients exhibited inferior multi-lineage differentiation potential with reduced Notch1 protein and mRNA level. Bisulfite sequencing PCR results showed the methylation level of Notch1 promoter was raised, which inhibited the binding of NF-Ya. Exposure BMSCs from NCHD patients with D-galactose under hypoxia (4% O₂) decreased the expression of Notch1. And activating Notch1 partially restored the deficient BMSCs of CCHD patients. In conclusion, the impaired multi-lineage differentiation potential of BMSCs from CCHD patients is owing to the decreased Notch1 level with a remarkable hypermethylation in its promoter region. Activated Notch1 signaling pathway could partially restore the deficient BMSCs in the CCHD patients, which may provide a new method on cell therapy in patients with CCHD.

Famine Exposure during the Fetal Period Increased the Risk of Dyslipidemia in Female Adults

Chronic diseases are gradually becoming public health challenges around the world. This study was designed to explore the relationship between early life exposure to famine in China and the risk of dyslipidemia during adulthood. A total of 604 subjects born from 1955 to 1965 received a health checkup and completed a questionnaire survey at the health management center of the tertiary hospital in Hefei, China, in 2013. A logistic regression model was used to assess dyslipidemia in famine-exposed and nonexposed individuals. Overall, exposure to famine during the fetal (OR 1.37; 95% CI: 0.90-2.10; and p = 0.14) and childhood (OR 0.97; 95% CI: 0.67-1.41; and p = 0.89) periods did not significantly increase the risk of dyslipidemia in adulthood compared with no exposure group. For females, exposure to famine during the fetal period (OR 2.00; 95% CI: 1.03-3.86; and p = 0.04) significantly increased the risk of dyslipidemia in adulthood compared to no exposure; however, this difference was not found in males. Exposure to famine in early life leads to altered lipid distribution in adulthood, and the risk of dyslipidemia significantly increased in adult women who were exposed to famine during the fetal period. Our study further validated the relationship between famine exposure during pregnancy and increased risk of dyslipidemia in female adult offspring. This study provides a scientific basis for the prevention and control of abnormal blood lipid levels in adults.

Epigenetics meets proteomics in an epigenome-wide association study with circulating blood plasma protein traits

DNA methylation and blood circulating proteins have been associated with many complex disorders, but the underlying disease-causing mechanisms often remain unclear. Here, we report an epigenome-wide association study of 1123 proteins from 944 participants of the KORA population study and replication in a multi-ethnic cohort of 344 individuals. We identify 98 CpG-protein associations (pQTMs) at a stringent Bonferroni level of significance. Overlapping associations with transcriptomics, metabolomics, and clinical endpoints suggest implication of processes related to chronic low-grade inflammation, including a network involving methylation of NLRC5, a regulator of the inflammasome, and associated pQTMs implicating key proteins of the immune system, such as CD48, CD163, CXCL10, CXCL11, LAG3, FCGR3B, and B2M. Our study links DNA methylation to disease endpoints via intermediate proteomics phenotypes and identifies correlative networks that may eventually be targeted in a personalized approach of chronic low-grade inflammation.

Clinical and lifestyle related factors influencing whole blood metabolite levels - A comparative analysis of three large cohorts

Objective

Human blood metabolites are influenced by a number of lifestyle and environmental factors. Identification of these factors and the proper quantification of their relevance provides insights into human biological and metabolic disease processes, is key for standardized translation of metabolite biomarkers into clinical applications, and is a prerequisite for comparability of data between studies. However, so far only limited data exist from large and well-phenotyped human cohorts and current methods for analysis do not fully account for the characteristics of these data. The primary aim of this study was to identify, quantify and compare the impact of a comprehensive set of clinical and lifestyle related factors on metabolite levels in three large human cohorts. To achieve this goal, we improve current methodology by developing a principled analysis approach, which could be translated to other cohorts and metabolite panels.

Methods

63 Metabolites (amino acids, acylcarnitines) were quantified by liquid chromatography tandem mass spectrometry in three cohorts (total N = 16,222). Supported by a simulation study evaluating various analytical approaches, we developed an analysis pipeline including preprocessing, identification, and quantification of factors affecting metabolite levels. We comprehensively identified uni- and multivariable metabolite associations considering 29 environmental and clinical factors and performed metabolic pathway enrichment and network analyses.

Results

Inverse normal transformation of batch corrected and outlier removed metabolite levels accompanied by linear regression analysis proved to be the best suited method to deal with the metabolite data. Association analyses revealed numerous uni- and multivariable significant associations. 15 of the analyzed 29 factors explained >1% of variance for at least one of the metabolites. Strongest factors are application of steroid hormones, reticulocytes, waist-to-hip ratio, sex, haematocrit, and age. Effect sizes of factors are comparable across studies.

Conclusions

We introduced a principled approach for the analysis of MS data allowing identification, and quantification of effects of clinical and lifestyle factors with metabolite levels. We detected a number of known and novel associations broadening our understanding of the regulation of the human metabolome. The large heterogeneity observed between cohorts could almost completely be explained by differences in the distribution of influencing factors emphasizing the necessity of a proper confounder analysis when interpreting metabolite associations.

•

Amino-acids and acylcarnitines analyzed in three studies with >16,000 individuals.

•

Develop a generic and adaptable bioinformatics workflow.

•

Analysis of the impact of 29 clinical and life-style factors on blood metabolites.

•

Analysis of network between factors and metabolites.

•

Comparison of results between studies.

Personalized Epigenome Remodeling Under Biochemical and Psychological Changes During Long-Term Isolation Environment

It has been reported that several aspects of human health could be disturbed during a long-term isolated environment (for instance, the Mars-500 mission), including psychiatric disorders, circadian disruption, temporal dynamics of gut microbiota, immune responses, and physical-activity-related neuromuscular performance. Nevertheless, the mechanisms underlying these disturbances and the interactions among different aspects of human adaptation to extreme environments remain to be elucidated. Epigenetic features, like DNA methylation, might be a linking mechanism that explains the involvement of environmental factors between the human genome and the outcome of health. We conducted an exploration of personalized longitudinal DNA methylation patterns of the peripheral whole blood cells, profiling six subjects across six sampling points in the Mars-500 mission. Specifically, we developed a Personalized Epigenetic-Phenotype Synchronization Analysis (PeSa) algorithm to explore glucose- and mood-state-synchronized DNA methylation sites, focusing on finding the dynamic associations between epigenetic patterns and phenotypes in each individual, and exploring the underling epigenetic connections between glucose and mood-state disturbance. Results showed that DMPs (differentially methylated-probes) were significantly enriched in pathways related to glucose metabolism (Type II diabetes mellitus pathway), mood state (Long-term depression) and circadian rhythm (Circadian entrainment pathway) during the mission. Furthermore, our data revealed individualized glucose-synchronized and mood-state-synchronized DNA methylation sites, and PTPRN2 was found to be associated with both glucose and mood state disturbances across all six subjects. Our findings suggest that personalized phenotype-synchronized epigenetic features could reflect the effects on the human body, including the disturbances of glucose and mood-states. The association analysis of DNA methylation and phenotypes, like the PeSa analysis, could provide new possibilities in understanding the intrinsic relationship between phenotypic changes of the human body adapting to long-term isolation environmental factors.

Metabolomic Pattern Predicts Incident Coronary Heart Disease

Objective- Alterations in the serum metabolome may be detectable in at-risk individuals before the onset of coronary heart disease (CHD). Identifying metabolomic signatures associated with CHD may provide insight into disease pathophysiology and prevention. Approach and Results- Metabolomic profiling (chromatography-mass spectrometry) was performed in 2232 African Americans and 1366 European Americans from the ARIC study (Atherosclerosis Risk in Communities). We applied Cox regression with least absolute shrinkage and selection operator to select metabolites associated with incident CHD. A metabolite risk score was constructed to evaluate whether the metabolite risk score predicted CHD risk beyond traditional risk factors. After 30 years of follow-up, we observed 633 incident CHD cases. Thirty-two metabolites were selected by least absolute shrinkage and selection operator to be associated with CHD, and 19 of the 32 showed significant individual associations with CHD, including a sugar substitute, erythritol. Theophylline (hazard ratio [95% CI] =1.16 [1.09-1.25]) and gamma-linolenic acid (hazard ratio [95% CI] =0.89 [0.81-0.97]) showed the greatest positive and negative associations with CHD, respectively. A 1 SD greater standardized metabolite risk score was associated with a 1.37-fold higher risk of CHD (hazard ratio [95% CI] =1.37 [1.27-1.47]). Adding the metabolite risk score to the traditional risk factors significantly improved model predictive performance (30-year risk prediction: Δ C statistics [95% CI] =0.016 [0.008-0.024], continuous net reclassification index [95% CI] =0.522 [0.480-0.556], integrated discrimination index [95% CI] =0.038 [0.019-0.065]). Conclusions- We identified 19 metabolites from known and novel metabolic pathways that collectively improved CHD risk prediction. Visual Overview- An online visual overview is available for this article.

Fine-Mapping of the Human Blood Plasma N-Glycome onto Its Proteome

Most human proteins are glycosylated. Attachment of complex oligosaccharides to the polypeptide part of these proteins is an integral part of their structure and function and plays a central role in many complex disorders. One approach towards deciphering this human glycan code is to study natural variation in experimentally well characterized samples and cohorts. High-throughput capable large-scale methods that allow for the comprehensive determination of blood circulating proteins and their glycans have been recently developed, but so far, no study has investigated the link between both traits. Here we map for the first time the blood plasma proteome to its matching N-glycome by correlating the levels of 1116 blood circulating proteins with 113 N-glycan traits, determined in 344 samples from individuals of Arab, South-Asian, and Filipino descent, and then replicate our findings in 46 subjects of European ancestry. We report protein-specific N-glycosylation patterns, including a correlation of core fucosylated structures with immunoglobulin G (IgG) levels, and of trisialylated, trigalactosylated, and triantennary structures with heparin cofactor 2 (SERPIND2). Our study reveals a detailed picture of protein N-glycosylation and suggests new avenues for the investigation of its role and function in the associated complex disorders.

Gestational route to healthy birth (GaRBH): protocol for an Indian prospective cohort study

INTRODUCTION

Pregnancy is characterised by a high rate of metabolic shifts from early to late phases of gestation in order to meet the raised physiological and metabolic needs. This change in levels of metabolites is influenced by gestational weight gain (GWG), which is an important characteristic of healthy pregnancy. Inadequate/excessive GWG has short-term and long-term implications on maternal and child health. Exploration of gestational metabolism is required for understanding the quantitative changes in metabolite levels during the course of pregnancy. Therefore, our aim is to study trimester-specific variation in levels of metabolites in relation to GWG and its influence on fetal growth and newborn anthropometric traits at birth.

METHODS AND ANALYSIS

A prospective longitudinal study is planned (start date: February 2018; end date: March 2023) on pregnant women that are being recruited in the first trimester and followed in subsequent trimesters and at the time of delivery (total 3 follow-ups). The study is being conducted in a hospital located in Bikaner district (66% rural population), Rajasthan, India. The estimated sample size is of 1000 mother-offspring pairs. Information on gynaecological and obstetric history, socioeconomic position, diet, physical activity, tobacco and alcohol consumption, depression, anthropometric measurements and blood samples is being collected for metabolic assays in each trimester using standardised methods. Mixed effects regression models will be used to assess the role of gestational weight in influencing metabolite levels in each trimester. The association of maternal levels of metabolites with fetal growth, offspring's weight and body composition at birth will be investigated using regression modelling.

ETHICS AND DISSEMINATION

The study has been approved by the ethics committees of the Department of Anthropology, University of Delhi and Sardar Patel Medical College, Rajasthan. We are taking written informed consent after discussing the various aspects of the study with the participants in the local language.

Network Medicine in Pathobiology

The past decade has witnessed exponential growth in the generation of high-throughput human data across almost all known dimensions of biological systems. The discipline of network medicine has rapidly evolved in parallel, providing an unbiased, comprehensive biological framework through which to interrogate and integrate systematically these large-scale, multi-omic data to enhance our understanding of disease mechanisms and to design drugs that reflect a deep knowledge of molecular pathobiology. In this review, we discuss the key principles of network medicine and the human disease network and explore the latest applications of network medicine in this multi-omic era. We also highlight the current conceptual and technological challenges, which serve as exciting opportunities by which to improve and expand the network-based applications beyond the artificial boundaries of the current state of human pathobiology.

Investigating the effects of lycopene and green tea on the metabolome of men at risk of prostate cancer: The ProDiet randomised controlled trial

Lycopene and green tea consumption have been observationally associated with reduced prostate cancer risk, but the underlying mechanisms have not been fully elucidated. We investigated the effect of factorial randomisation to a 6-month lycopene and green tea dietary advice or supplementation intervention on 159 serum metabolite measures in 128 men with raised PSA levels (but prostate cancer-free), analysed by intention-to-treat. The causal effects of metabolites modified by the intervention on prostate cancer risk were then assessed by Mendelian randomisation, using summary statistics from 44,825 prostate cancer cases and 27,904 controls. The systemic effects of lycopene and green tea supplementation on serum metabolic profile were comparable to the effects of the respective dietary advice interventions (R2 = 0.65 and 0.76 for lycopene and green tea respectively). Metabolites which were altered in response to lycopene supplementation were acetate [β (standard deviation difference vs. placebo): 0.69; 95% CI = 0.24, 1.15; p = 0.003], valine (β: -0.62; -1.03, -0.02; p = 0.004), pyruvate (β: -0.56; -0.95, -0.16; p = 0.006) and docosahexaenoic acid (β: -0.50; -085, -0.14; p = 0.006). Valine and diacylglycerol were lower in the lycopene dietary advice group (β: -0.65; -1.04, -0.26; p = 0.001 and β: -0.59; -1.01, -0.18; p = 0.006). A genetically instrumented SD increase in pyruvate increased the odds of prostate cancer by 1.29 (1.03, 1.62; p = 0.027). An intervention to increase lycopene intake altered the serum metabolome of men at risk of prostate cancer. Lycopene lowered levels of pyruvate, which our Mendelian randomisation analysis suggests may be causally related to reduced prostate cancer risk.

Paving the Way to Precision Nutrition Through Metabolomics

Nutrition is an interdisciplinary science that studies the interactions of nutrients with the body in relation to maintenance of health and well-being. Nutrition is highly complex due to the underlying various internal and external factors that could model it. Thus, hacking this complexity requires more holistic and network-based strategies that could unveil these dynamic system interactions at both time and space scales. The ongoing omics era with its high-throughput molecular data generation is paving the way to embrace this complexity and is deeply reshaping the whole field of nutrition. Understanding the future paths of nutrition science is of importance from both translational and clinical perspectives. Basic nutrients which might include metabolites are important in nutrition science. Moreover, metabolites are key biological communication channels and represent an appealing functional readout at the interface of different major influential factors that define health and disease. Metabolomics is the technology that enables holistic and systematic analyses of metabolites in a biological system. Hence, given its intrinsic functionality, its tight connection to metabolism and its high clinical actionability potential, metabolomics is a very appealing technology for nutrition science. The ultimate goal is to deliver a tailored and clinically relevant nutritional recommendations and interventions to achieve precision nutrition. This work intends to present an update on the applications of metabolomics to personalize nutrition in translational and clinical settings. It also discusses the current conceptual shifts that are remodeling clinical nutrition practices in this Precision Medicine era. Finally, perspectives of clinical nutrition in the ever-growing, data-driven healthcare landscape are presented.