CpG site-specific methylation as epi-biomarkers for the prediction of health risk in PAHs-exposed populations

DNMT1 regulates hypermethylation and silences hsa_circ_401351 in hydroquinone-induced malignant TK6 cells

BACKGROUND

Benzene and its metabolite hydroquinone (HQ) are widely used in daily life, and long-term exposure to benzene or HQ can induce acute myeloid leukemia (AML). Circular RNAs (circRNAs) are mostly produced by reverse splicing of gene exon mRNA precursors. The modulation of circRNA expression is connected to leukemia progression; however, the molecular mechanism is still unknown.

MATERIALS AND METHODS

In this study, the cells were divided into four groups: PBS control group (PBS-TK6), TK6 malignantly transformed cells induced by 10.0 μmol/L HQ (HQ-TK6), and HQ-TK6 cells treated with 5 μmol/L 5-AzaC (DNA methyltransferase inhibitor) for 24 h (HQ + 5-AzaC). HQ-TK6 cells were treated with 200 nmol/L TSA (histone deacetylation inhibitor) for 24 h (HQ + TSA). qRT-PCR was used to identify the differential hsa_circ_401351 expression between the four groups. We further determined the hsa_circ_401351 promoter methylation level with methylation-specific PCR. DNMT1 and DNMT3b were knocked down by CRISPR/Cas9 to elucidate the specific molecular mechanism of hsa_circ_401351 in HQ-TK6 cells. CCK-8 and flow cytometry detected cell proliferation and apoptosis, respectively, after hsa_circ_401351 was overexpressed in HQ-TK6 cells.

RESULTS

Compared with the PBS-TK6 group, the expression of hsa_circ_401351 was found to be lower in the HQ-TK6 group. Nevertheless, treatment with 5-AzaC or TSA increased hsa_circ_401351 expression, with the upregulation being more pronounced in the TSA group. The expression of hsa_circ_401351 in the DNMT1 knockdown group was dramatically increased by 50% compared to that in the control group, and the DNA methylation level of the hsa_circ_401351 promoter region was decreased. When hsa_circ_401351 was overexpressed, HQ-TK6 cell proliferation was significantly slowed after 48 h compared with the control group. Flow cytometry showed that cells were mainly arrested in G1 phase, and apoptosis was significantly enhanced. Similarly, qRT-PCR and Western blot data showed significant reductions in Caspase-3 mRNA and protein production, and Bcl-2 mRNA levels were also elevated.

CONCLUSIONS

Overall, our research showed that elevated DNMT1 expression in HQ-TK6 cells increased methylation levels and decreased expression of the hsa_circ_401351 promoter region, limiting its ability to suppress HQ-TK6 cell growth and enhance apoptosis.

Specific CpG sites methylation is associated with hematotoxicity in low-dose benzene-exposed workers

Benzene is a broadly used industrial chemicals which causes various hematologic abnormalities in human. Altered DNA methylation has been proposed as epigenetic biomarkers in health risk evaluation of benzene exposure, yet the role of methylation at specific CpG sites in predicting hematological effects remains unclear. In this study, we recruited 120 low-level benzene-exposed and 101 control male workers from a petrochemical factory in Maoming City, Guangdong Province, China. Urinary S-phenylmercapturic acid (SPMA) in benzene-exposed workers was 3.40-fold higher than that in control workers (P < 0.001). Benzene-induced hematotoxicity was characterized by reduced white blood cells counts and nuclear division index (NDI), along with an increased DNA damage and urinary 8-hydroxy-2'-deoxyguanosine (all P < 0.05). Methylation levels of TRIM36, MGMT and RASSF1a genes in peripheral blood lymphocytes (PBLCs) were quantified by pyrosequencing. CpG site 6 of TRIM36, CpG site 2, 4, 6 of RASSF1a and CpG site 1, 3 of MGMT methylation were recognized as hot CpG sites due to a strong correlation with both internal exposure and hematological effects. Notably, integrating hot CpG sites methylation of multiple genes reveal a higher efficiency in prediction of integrative damage compared to individual genes at hot CpG sites. The negative dose-response relationship between the combined methylation of hot CpG sites in three genes and integrative damage enabled the classification of benzene-exposed individuals into high-risk or low-risk groups using the median cut-off value of the integrative index. Subsequently, a prediction model for integrative damage in benzene-exposed populations was built based on the methylation status of the identified hot CpG sites in the three genes. Taken together, these findings provide a novel insight into application prospect of specific CpG site methylation as epi-biomarkers for health risk assessment of environmental pollutants.

Quantitative relationships of FAM50B and PTCHD3 methylation with reduced intelligence quotients in school aged children exposed to lead: Evidence from epidemiological and in vitro studies

At present, the application of DNA methylation (DNAm) biomarkers in environmental health risk assessment (EHRA) is more challenging due to the unclearly quantitative relationship between them. We aimed to explore the role of FAM50B and PTCHD3 at the level of signaling pathways, and establish the quantitative relationship between them and children's intelligence quotients (IQs). DNAm of target regions was measured in multiple cell models and was compared with the human population data. Then the dose-response relationships of lead exposure with neurotoxicity and DNAm were established by benchmark dose (BMD) model, followed by potential signaling pathway screening. Results showed that there was a quantitative linear relationship between children's IQs and FAM50B/PTCHD3 DNAm (DNAm between 51.40 % - 78.78 % and 31.41 % - 74.19 % for FAM50B and PTCHD3, respectively), and this relationship was more significant when children's IQs > 90. The receiver operating characteristic (ROC) and calibration curves showed that FAM50B/PTCHD3 DNAm had a satisfying accuracy and consistency in predicting children's IQs, which was confirmed by sensitivity analysis of gender and CpG site grouping data. In cell experiments, there was also a quantitative linear relationship between FAM50B DNAm and reactive oxygen species (ROS) production, which was mediated by PI3K-AKT signaling pathway. In addition, the lead BMD of ROS was close to that of FAM50B DNAm, suggesting that FAM50B DNAm was a suitable biomarker for the risk assessments of adverse outcomes induced by lead. Taken collectively, these results suggest that FAM50B/PTCHD3 can be applied to EHRA and the prevention/intervention of adverse effects of lead on children's IQs.

Benzo[b]fluoranthene induces male reproductive toxicity and apoptosis via Akt-Mdm2-p53 signaling axis in mouse Leydig cells: Integrating computational toxicology and experimental approaches

This study aims to explore the male reproductive toxicity of Benzo[b]fluoranthene (BbF) and related mechanisms. The results of computational toxicology analysis indicated male reproductive toxicity of BbF was related to apoptosis of Leydig cells and that Akt/p53 pathway might play a key role. In experiments, BbF induced testosterone decline, decreased concentration and motility of sperm and aggravated testicular pathological injury in mice. Besides, BbF led to apoptosis in Leydig cells, and decreased expressions of p-Akt and Bcl2, while improving the expressions of p53, Bax and Cleaved Caspase-3 in vivo and in vitro. Further, compared with BbF group, Akt activator SC79 significantly reduced cell apoptosis rate, improved cell viability, promoted the expressions of p-Akt and p-Mdm2, and reversed the above molecular expressions. Similarly, p53 inhibitor Pifithrin-α also significantly enhanced the cell vitality, alleviated the apoptosis of TM3 cells induced by BbF, and decreased the expressions of Bax and Cleaved Caspase-3, with the up-regulation of Bcl2. To sum up, by inhibiting Akt-Mdm2 signaling, BbF activated the p53-mediated mitochondrial apoptosis pathway, further inducing the apoptosis of Leydig cells, therefore resulting in testosterone decline and male reproductive damage. Besides, this study provided a valid mode integrating computational toxicology and experimental approaches in toxicity testing.

Diesel exhaust particles exposure induces liver dysfunction: Exploring predictive potential of human circulating microRNAs signature relevant to liver injury risk

Diesel exhaust particles (DEP) pollution should be taken seriously because it is an extensive environmental and occupational health concern. Exploring early effect biomarkers is crucial for monitoring and managing DEP-associated health risk assessment. Here, we found that serum levels of 67 miRNAs were dysregulated in DEP exposure group. Notably, 20 miRNAs were identified as each having a significant dose-response relationship with the internal exposure level of DEP. Further, we revealed that the DEP exposure could affect the liver function of subjects and that 7 miRNAs (including the well-known liver injury indicator, miR-122-5p) could serve as the novel epigenetic-biomarkers (epi-biomarkers) to reflect the liver-specific response to the DEP exposure. Importantly, an unprecedented prediction model using these 7 miRNAs was established for the assessment of DEP-induced liver injury risk. Finally, bioinformatic analysis indicated that the unique set of miRNA panel in serum might also contribute to the molecular mechanism of DEP exposure-induced liver damage. These results broaden our understanding of the adverse health outcomes of DEP exposure. Noteworthy, we believe this study could shed light on roles and functions of epigenetic biomarkers from environmental exposure to health outcomes by revealing the full chain of exposure-miRNAs-molecular pathways-disease evidence.

Human health risk assessment of 6:2 Cl-PFESA through quantitative in vitro to in vivo extrapolation by integrating cell-based assays, an epigenetic key event, and physiologically based pharmacokinetic modeling

Human health risk assessment of chemicals is essential but often relies on time-consuming and animal and labor-extensive procedures. Here, we develop a population-based, quantitative in vitro to in vivo extrapolation (QIVIVE) approach which depended on cellular effects monitored by in vitro assays, considered chemical internal concentration determined by LC-MS/MS, extrapolated into in vivo target tissue concentration through physiologically based pharmacokinetic (PBPK) modelling, and assessed populational health risk using in silico modelling. By applying this QIVIVE approach to 6:2 chlorinated polyfluorinated ether sulfonate (6:2 Cl-PFESA), as a representative of the emerging pollutants, we find that 6:2 Cl-PFESA disturbed lipid homeostasis in HepG2 cells through enhancement of lipid accumulation and fatty acid β-oxidation, during which miR-93-5p served as a key event towards toxicity and thus, could serve as an efficient toxicity marker for risk assessment; further, the disruption potency of lipid homeostasis of 6:2 Cl-PFESA for the most of studied populations in China might be of moderate concern. Together, our approach improved the reliability of QIVIVE during human health risk assessment, which can readily be used for other chemicals.

Genotoxicity of organic contaminants in the soil: A review based on bibliometric analysis and methodological progress

Organic contaminants (OCs) are ubiquitous in the environment, posing severe threats to human health and ecological balance. In particular, OCs and their metabolites could interact with genetic materials to induce genotoxicity, which has attracted considerable attention. In this review, bibliometric analysis was executed to analyze the publications on the genotoxicity of OCs in soil from 1992 to 2021. The result indicated that significant contributions were made by China and the United States in this field and the research hotspots were biological risks, damage mechanisms, and testing methods. Based on this, in this review, we summarized the manifestations and influencing factors of genotoxicity of OCs to soil organisms, the main damage mechanisms, and the most commonly utilized testing methods. OCs can induce genotoxicity and the hierarchical response of soil organisms, which could be influenced by the physicochemical properties of OCs and the properties of soil. Specific mechanisms of genotoxicity can be classified into DNA damage, epigenetic toxicity, and chromosomal aberrations. OCs with different molecular weights lead to genetic material damage by inducing the generation of ROS or forming adducts with DNA, respectively. The micronucleus test and the comet test are the most commonly used testing methods. Moreover, this review also pointed out that future studies should focus on the relationships between bioaccessibilities and genotoxicities, transcriptional regulatory factors, and potential metabolites of OCs to elaborate on the biological risks and mechanisms of genotoxicity from an overall perspective.

Effect and interaction of TNKS genetic polymorphisms and environmental factors on telomere damage in COEs-exposure workers

Coke oven emissions (COEs) contain many carcinogenic polycyclic aromatic hydrocarbons (PAHs). Telomere damage is an early biological marker reflecting long-term COEs-exposure. Whereas, whether the genetic variations of telomere-regulated gene TNKS have an effect on the COEs-induced telomere damage is unknown. So we detected the environmental exposure levels, relative telomere length (RTL), and TNKS genetic polymorphisms among 544 COEs-exposure workers and 238 healthy participants. We found that the RTL of the wild homozygous GG genotype in rs1055328 locus was statistically shorter compared with the CG+CC genotype for the healthy participants using covariance analysis(P = 0.008). In the Generalized linear model (GLM) analysis, TNKS rs1055328 GG could accelerate telomere shortening (P = 0.011); and the interaction between TNKS rs1055328 GG and COEs-exposure had an effect on RTL (P = 0.002). In conclusion, this study was the first to discover the role of TNKS rs1055328 locus in COEs-induced telomere damage, and proved that chromosomal damage was a combined consequence of environmental and genetic factors.

The hypermethylation of FOXP3 gene as an epigenetic marker for the identification of arsenic poisoning risk

Background and Purpose: Arsenic exposure can lead to skin lesions and multiple organ damage, which are not easily reversible and for which there is no effective therapeutics. Identification of reliable epigenetic markers is essential for early recognition of arsenic poisoning risk. Anomalous DNA methylation of immune homeostasis regulator FOXP3 is a critical mechanism for triggering arsenic poisoning. This study aims to explore the value of FOXP3 methylation in the identification of arsenic poisoning risk.Methods: 88 arsenic poisoning subjects and 41 references were recruited. Urinary arsenic contents and FOXP3 methylation in PBLCs was measured by ICP-MS and pyrosequencing, respectively.Results: The results showed that the elevated FOXP3 methylation in PBLCs were associated with the increased levels of urinary arsenic and were positively associated with the increased risk of arsenic poisoning and its progression. The result of mediation analysis revealed that 24.3% of the effect of arsenic exposure on the risk of arsenic poisoning was mediated by increased FOXP3 methylation. Additionally, we constructed a nomogram model with FOXP3 methylation as an epigenetic predictor to assess the probability of individual arsenic poisoning. The model showed a robust ability in the discrimination of arsenic poisoning risk, with an area under receiver operating characteristics curve of 0.897(0.845-0.949) and more than 70% accuracy. The calibration curves and the Harrell concordance index showed that the consistency rate between the probability predicted by the nomogram model and the actual probability is 89.7%.Conclusions: Taken together, we found the great potential of FOXP3 methylation for the identification of arsenic poisoning risk and provided a new approach to the application of epigenetic markers in accurately quantifying the risk of adverse outcomes.