Map of benzoapyrene metabolites-DNA adducts in human bronchial epithelial-like cells: Based on chromatin immunoprecipitation followed by sequencing technology

Background The active metabolite of benzoapyrene (BaP), 7,8-dihydroxy-9,10-epoxybenzoapyrene (BPDE), can form adducts with DNA, but the spectrum of BPDE-DNA adducts is unclear.

Objective To identify the distribution of BPDE adduct sites and associated genes at the whole-genome level by chromatin immunoprecipitation followed by sequencing (ChIP-Seq), and serve as a basis for further exploring the toxicological mechanisms of BaP.

Methods Human bronchial epithelial-like cells (16HBE) were cultured to the fourth generation inthe logarithmic growth phase. Cells were harvested and added to chromatin immunoprecipitation lysis buffer. The lysate was divided into experimental and control groups. The experimental group received a final concentration of 20 μmol·L⁻¹ BPDE solution, while the control group received an equivalent volume of dimethyl sulfoxide solution. The cells were then incubated at 37 °C for 24 h. Chromatin fragments of 100-500 bp were obtained through sonication. BPDE-specific antibody (anti-BPDE 8E11) was used to enrich DNA fragments with BPDE adducts. High-throughput sequencing was conducted to detect BPDE adduct sites. The top 1000 peak sequences were subjected to motif analysis using MEME and DREME software. BPDE adduct target genes at the whole-genome level were annotated, and Gene Ontology (GO) functional analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of BPDE adduct target genes were conducted using bioinformatics techniques.

Results The high-throughput sequencing detected a total of 842 BPDE binding sites, distributed across various chromosomes. BPDE covalently bound to both coding and non-coding regions of genes, with 73.9% binding sites located in intergenic regions, 19.6% in intronic regions, and smaller proportions in upstream 2 kilobase, exonic, downstream 2 kilobase, and 5' untranslated regions. Regarding the top 1000 peak sequences, four reliable motifs were identified, revealing that sites rich in adenine (A) and guanine (G) were prone to binding. Through the enrichment analysis of binding sites, a total of 199 BPDE-adduct target genes were identified, with the majority located on chromosomes 1, 5, 7, 12, 17, and X. The GO analysis indicated that these target genes were mainly enriched in nucleic acid and protein binding, participating in the regulation of catalytic activity, transport activity, translation elongation factor activity, and playing important roles in cell division, differentiation, motility, substance transport, and information transfer. The KEGG analysis revealed that these target genes were primarily enriched in pathways related to cardiovascular diseases, cancer, and immune-inflammatory responses.

Conclusion Using ChIP-Seq, 199 BPDE adduct target genes at genome-wide level are identified, impacting biological functions such as cell division, differentiation, motility, substance transport, and information transfer. These genes are closely associated with cardiovascular diseases, tumors, and immune-inflammatory responses.