Identification of the hormetic dose-response and regulatory network of multiple metals co-exposure-related hypertension via integration of metallomics and adverse outcome pathways

Environmental stressors, including heavy metals, can be associated with hypertension development. However, little information regarding the dose-response relationship and toxicity mechanisms of metal mixtures with hypertension development is currently available. Therefore, we recruited 940 participants from six factories in northeastern China and measured the urinary concentrations of 19 metals. Then, we used Bayesian kernel machine regression (BKMR) to explore associations between metals co-exposure and hypertension. The BKMR model indicated a hermetic dose-response relationship between eight urinary metals (Co, Cr, Ni, Cd, As, Fe, Zn, and Pb) and hypertension risk. Moreover, heterogeneous and non-linear association patterns were detected across different metals/metalloids concentrations. Next, for the first time, we analyzed data of chemicals containing specific metal elements in the Comparative Toxicogenomics Database (CTD) from a disease perspective and provided insights from various biological levels to explain heavy metal co-exposure-related hypertension. On the molecular scale, 43 chemical components and 112 potential target genes were detected for metal exposure-related hypertension. Further, the network topology analysis indicated that target genes such as insulin (INS, degree = 78), albumin (ALB, degree = 74), renin (REN, degree = 71), interleukin-6 (IL6, degree = 70), endothelin 1 (EDN1, degree = 70), and endothelial nitric oxide synthase (NOS3, degree = 69) have a strong correlation with heavy metals co-exposure. Finally, we used integrative analyses in the adverse outcome pathway (AOP) wiki to analyze the co-exposure of heavy metals and hypertension and support an integrated metallomics approach. We selected the AOP 149 as the framework and found that the molecular initiating events (MIEs) of hypertension stems from the oxidation of AA residues on critical peptides of the NO pathway. The NOS3 was particularly promising since its subunit has three metal ion cross-linking domains with Zn