Genetic variants with gene regulatory effects are associated with diisocyanate-induced asthma

Background Isocyanates are major causes of occupational asthma, but susceptibility and mechanisms of diisocyanate-induced asthma (DA) remain uncertain. Objective The aim of this study was to identify DA-associated functional genetic variants through next-generation sequencing (NGS), bioinformatics, and functional assays. Methods NGS was performed in 91 workers with DA. Fourteen loci with known DA-associated single nucleotide polymorphisms (SNPs) were sequenced and compared with data from 238 unexposed subjects. Ranking of DA-associated SNPs based on their likelihood to affect gene regulatory mechanisms in the lung yielded 21 prioritized SNPs. Risk and nonrisk oligonucleotides were tested for binding of nuclear extracts from A549, BEAS-2B, and IMR-90 lung cell lines by using electrophoretic mobility shift assays. DNA constructs were cloned into a pGL3 promoter vector for luciferase gene reporter assays. Results NGS detected 130 risk variants associated with DA (3.1 × 10 −6 to 6.21 × 10 −4 ), 129 of which were located in noncoding regions. The 21 SNPs prioritized by using functional genomic data sets were in or proximal to 5 genes: cadherin 17 ( CDH17 ; n = 10), activating transcription factor 3 ( ATF3 ; n = 7), family with sequence similarity, member A ( FAM71A ; n = 2), tachykinin receptor 1 ( TACR1 ; n = 1), and zinc finger and BTB domain-containing protein 16 ( ZBTB16 ; n = 1). Electrophoretic mobility shift assays detected allele-dependent nuclear protein binding in A549 cells for 8 of 21 variants. In the luciferase assay 4 of the 21 SNPs exhibited allele-dependent changes in gene expression. DNA affinity precipitation and mass spectroscopy of rs147978008 revealed allele-dependent binding of H1 histones, which was confirmed by using Western blotting. Conclusions We identified 5 DA-associated potential regulatory SNPs. Four variants exhibited effects on gene regulation ( ATF rs11571537, CDH17 rs2446824 and rs2513789, and TACR1 rs2287231). A fifth variant ( FAM71A rs147978008) showed nonrisk allele preferential binding to H1 histones. These results demonstrate that many DA-associated genetic variants likely act by modulating gene regulation.