Abstract:Background The disease-associated non-coding variants identified by genome-wide association studies (GWASs) were enriched in open chromatin regions (OCRs) and implicated in gene regulation. Genetic variants in OCRs thus may exert regulatory functions and contribute to non-small cell lung cancer (NSCLC) susceptibility. Objective To fine map potential functional variants in GWAS loci that contribute to NSCLC predisposition using chromatin accessibility and histone modification data and explore their functions by population study and biochemical experimental analyses. Methods We mapped the chromatin accessible regions of lung tissues using data of assay for transposase-accessible chromatin using sequencing (ATAC-seq) in The Cancer Genome Atlas (TCGA) and prioritized potential regulatory variants within lung cancer GWAS loci by aligning with histone signatures using data of chromatin immunoprecipitation assays followed by sequencing (ChIP-seq) in the Encyclopedia of DNA Elements (ENCODE). A two-stage case–control study with 1,830 cases and 2,001 controls was conducted to explore the associations between candidate variants and NSCLC risk in Chinese population. Bioinformatic annotations and biochemical experiments were performed to further reveal the potential functions of significant variants. Results Sixteen potential functional single-nucleotide polymorphisms (SNPs) were selected as candidates from bioinformatics analyses. Three variants out of the 16 candidate SNPs survived after genotyping in stage 1 case–control study, and only the results of SNP rs13064999 were successfully validated in the analyses of stage 2 case–control study. In combined analyses, rs13064999 was significantly associated with NSCLC risk [additive model; odds ratio (OR) = 1.17; 95%CI, 1.07–1.29; p = 0.001]. Functional annotations indicated its potential enhancer bioactivity, and dual-luciferase reporter assays revealed a significant increase in luciferase activity for the reconstructed plasmid with rs13064999 A allele, when compared to the one with wild-type G allele ( p A549 < 0.001, p SK-MES-1 = 0.004). Further electrophoretic mobility shift assays (EMSA) and super-shift assays confirmed a stronger affinity of HP1γ for the binding motif containing SNP rs13064999 A allele. Conclusion These findings suggested that the functional variant rs13064999, identified by the integration of ATAC-seq and ChIP-seq data, contributes to the susceptibility of NSCLC by affecting HP1γ binding, while the exact biological mechanism awaits further exploration.

Integrating apaQTL and eQTL analysis identifies a potential causal variant associated with lung adenocarcinoma risk in the Chinese population

Association Between the LINC00673 Rs11655237 C> T Polymorphisms with Cancer Risk in the Chinese Population: A Meta-Analysis

SNP Rs17079281 Decreases Lung Cancer Risk Through Creating an YY1-Binding Site to Suppress DCBLD1 Expression

Colocalization analysis of 3′ UTR alternative polyadenylation quantitative trait loci reveals novel mechanisms underlying associations with lung function

Crosstalk among Alternative Polyadenylation, Genetic Variants and Ubiquitin Modification Contribute to Lung Adenocarcinoma Risk

Genetic variants that impact alternative polyadenylation in cancer represent candidate causal risk loci

Alternative polyadenylation-related genetic variants contribute to bladder cancer risk

The Association of Rs710886 in Lncrna PCAT1 with Bladder Cancer Risk in a Chinese Population.

Identification of novel pQTL‐SNPs associated with lung adenocarcinoma risk: A multi‐stage study

Large-scale alternative polyadenylation-wide association studies to identify putative cancer susceptibility genes

Large-scale alternative polyadenylation (APA)-wide association studies to identify putative susceptibility genes in human common cancers

A Functional Polymorphism in Accessible Chromatin Region Confers Risk of Non-Small Cell Lung Cancer in Chinese Population

Genetic Basis of Alternative Polyadenylation is an Emerging Molecular Phenotype for Human Traits and Diseases

Alternative polyadenylation shapes the molecular and clinical features of lung adenocarcinoma

A novel APA-based prognostic signature may predict the prognosis of lung adenocarcinoma in an East Asian population

An Atlas of Alternative Polyadenylation Quantitative Trait Loci Contributing to Complex Trait and Disease Heritability

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

A genome-wide association study identifies two new lung cancer susceptibility loci at 13q12.12 and 22q12.2 in Han Chinese

rs401681 and rs402710 confer lung cancer susceptibility by regulating TERT expression instead of CLPTM1L in East Asian populations

Genome-wide association scan of tag SNPs identifies a susceptibility locus for lung cancer at 15q25.1

Genome‐Wide 3′‐UTR Single Nucleotide Polymorphism Association Study Identifies Significant Prostate Cancer Risk‐Associated Functional Loci at 8p21.2 in Chinese Population