Interrogation of functional variants in COPD GWAS loci by massively parallel reporter assays

Genome-wide association study (GWAS) loci often contain many linked variants, making it difficult to determine which variant is functionally relevant. Massively parallel reporter assays (MPRA) allow experimental testing of candidate variants to identify those with regulatory activity. Prior chronic obstructive pulmonary disease (COPD) MPRA studies have largely focused on individual loci, whereas broader multi-locus, multi-cell-type interrogation remains limited.

We aim to identify functional variants in five COPD GWAS loci across three lung-relevant cell types.

We screened 1120 variants using MPRA in epithelial (16HBE), fibroblast (MRC5), and endothelial (HUVEC) cells followed by reporter assay validation. Public Hi-C, ChIP-seq and ATAC-seq datasets were analyzed to evaluate chromatin context near candidate variants. We further performed CRISPR interference (CRISPRi) targeting variant-containing regions and measured gene expression by RT-qPCR in primary normal human bronchial epithelial (NHBE) cells using two gRNAs per variant. Co-immunoprecipitation was performed to test interaction between selected candidate genes.

In MPRA, we identified 25 variants with allele-specific effects (∼2% of tested variants). Enrichment of H3K27Ac and open chromatin near rs35421223 was detected in 16HBE cells. CRISPRi identified two SNP-gene pairs, RUVBL1 and RAB7A regulated by rs35421223 in both the 16HBE cell line and primary NHBE cells. We detected interaction between RUVBL1 and the known COPD gene product FAM13A.

Screening COPD loci across three cell types identified functional regulatory variants and linked them to candidate target genes for future mechanistic studies.