Integration of apaQTL and eQTL analysis reveals novel SNPs associated with occupational pulmonary fibrosis risk

To explore the association between apaQTL/eQTL-SNPs and the susceptibility to silicosis. A silicosis-related GWAS was initially conducted to screen for single nucleotide polymorphisms (SNPs) associated with the risk of silicosis. Candidate SNPs with apaQTL and eQTL functions were then obtained from the 3'aQTL-atlas and GTEx databases. Subsequently, additional case-control studies were performed to validate the relationship between the candidate apaQTL/eQTL-SNPs and the risk of silicosis. Finally, experiments were conducted to illustrate APA events occurring at different alleles of the identified apaQTL/eQTL-SNPs. The combined results of the GWAS and iMLDR validations indicate that the variant T allele of the rs2974341 located on SMIM19 (additive model: OR = 0.66, the 95% CI = 0.53-0.84, P = 0.001) and the variant T allele of the rs2390488 located on TMTC4 (additive model: OR = 0.72, 95% CI = 0.57-0.90, P = 0.005) were significantly associated with decreased risk of developing silicosis susceptibility. Furthermore, 3'RACE experiments verified the presence of two poly (A) sites (proximal and distal) in SMIM19, rs2974341 may remotely regulate the binding between miRNA-3646 and SMIM19 with its high LD locus rs2974353 to affect the expression level of SMIM19. The rs2974341 variant T allele may contribute to the generation of the shorter 3'UTR transcript of SMIM19 and affect the binding of miRNA-3646 to the target gene SMIM19. The apaQTL/eQTL-SNPs may provide new perspectives for evaluating the regulatory function of SNPs in the development of silicosis.