Integrative splicing-quantitative-trait-locus analysis reveals risk loci for non-small-cell lung cancer

Splicing quantitative trait loci (sQTLs) have been demonstrated to contribute to disease etiology by affecting alternative splicing. However, the role of sQTLs in the development of non-small-cell lung cancer (NSCLC) remains unknown. Thus, we performed a genome-wide sQTL study to identify genetic variants that affect alternative splicing in lung tissues from 116 individuals of Chinese ancestry, which resulted in the identification of 1,385 sQTL-harboring genes (sGenes) containing 378,210 significant variant-intron pairs. A comprehensive characterization of these sQTLs showed that they were enriched in actively transcribed regions, genetic regulatory elements, and splicing-factor-binding sites. Moreover, sQTLs were largely distinct from expression quantitative trait loci (eQTLs) and showed significant enrichment in potential risk loci of NSCLC. We also integrated sQTLs into NSCLC GWAS datasets (13,327 affected individuals and 13,328 control individuals) by using splice-transcriptome-wide association study (spTWAS) and identified alternative splicing events in 19 genes that were significantly associated with NSCLC risk. By using functional annotation and experiments, we confirmed an sQTL variant, rs35861926, that reduced the risk of lung adenocarcinoma (rs35861926-T, OR = 0.88, 95% confidence interval [CI]: 0.82-0.93, p = 1.87 × 10^-5) by promoting FARP1 exon 20 skipping to downregulate the expression level of the long transcript FARP1-011. Transcript FARP1-011 promoted the migration and proliferation of lung adenocarcinoma cells. Overall, our study provided informative lung sQTL resources and insights into the molecular mechanisms linking sQTL variants to NSCLC risk.