Abstract:Abstract Background. By SNP genotyping and RNA sequencing of 471 normal prostate samples, we recently created a prostate tissue-based eQTL dataset and identified significant eQTL signals at 51 prostate cancer risk loci. To functionally characterize these risk SNPs, we developed a massively parallel sequencing technology to screen SNPs for their allele-dependent protein binding differences. We combined this technology (called SNPs-seq) with another high throughput assay (called STARR-seq) to screen the risk loci with significant prostate-specific eQTL signals. Methods. To select candidate functional SNPs in eQTL regions, we took advantage of existing epigenomic datasets and available tools including ENCODE, HaploReg, and Regulome. For all selected SNPs, we first made allele-specific double-strand oligos and performed DNA-protein binding assays. We then performed sequencing analysis on the protein-bound DNA oligos and determined allele-specific protein binding differences. To evaluate reproducibility of SNPs-seq, we performed each assay in duplicates. We cloned SNPs-seq screened SNP regions showing allele-specific protein binding differences into the STARR-seq vector to further determine allele-specific enhancer activities. Finally, we performed EMSA and luciferase reporter assays to validate a set of promising candidate SNPs. Results. From 51 risk loci with strong eQTL signals, we selected 374 SNPs with strong indication of regulatory potential, as evidenced by overlapping with epigenomic marks. When comparing technical duplicates, sequence read counts from the SNPs-seq showed significant correlation with r2>=0.99. By normalizing input controls, we found 101 of the 374 SNPs showing significant allelic protein binding differences (>=1.5-fold binding difference between variant and reference alleles). Interestingly, three published functional SNPs (rs12769019, rs10993994, and rs4907792) were also among the significant SNPs, validating SNPs-seq as functional SNP screening tool. To further validate the candidate SNPs from SNPs-seq, we applied STARR-seq and tested the 101 SNPs-containing sequences (371-686bp) in LNCaP cell line under androgen treatment. This analysis revealed 11 SNPs that not only demonstrated enhancer/repressor activity but also functioned with allelic differences. EMSA and luciferase reporter assays confirmed 6 SNPs with allele-dependent enhancer/repressor activity. Conclusions. We developed a high throughput sequencing-based technology to screen large number candidate SNPs for their allelic protein binding differences. The SNPs-seq coupled with STARR-seq will provide a powerful strategy for functionally characterizing risk loci in prostate cancer and other common diseases. Further understanding genetic role of prostate cancer etiology may facilitate the translation of population-based discovery into biological mechanisms and eventually benefit clinical practice. Citation Format: Peng Zhang, Jing Zhu, Sufyan Suleman, Yong-Chen Guo, Mei-Jun Du, Li-Dong Wang, Gong-Hong Wei, Liang Wang. Functional characterization of prostate cancer risk loci by SNPs-seq and STARR-seq [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 1280. doi:10.1158/1538-7445.AM2017-1280

Optimized high-throughput screening of non-coding variants identified from genome-wide association studies

High-throughput Screening of Prostate Cancer Risk Loci by Single Nucleotide Polymorphisms Sequencing.

Abstract 1280: Functional Characterization of Prostate Cancer Risk Loci by SNPs-seq and STARR-seq

Systematic identification of regulatory variants associated with cancer risk

Identification and validation of regulatory SNPs that modulate transcription factor chromatin binding and gene expression in prostate cancer

Proteogenomic and Clinical Implications of Unique Recurrent Splice Variants in Clear Cell Renal Cell Carcinoma.

Post genome-wide association studies functional characterization of prostate cancer risk loci

High-throughput characterization of genetic effects on DNA-protein binding and gene transcription

DrGaP: a powerful tool for identifying driver genes and pathways in cancer sequencing studies.

Combined SNPs Sequencing and Allele Specific Proteomics Capture Reveal Functional Causality Underpinning the 2p25 Prostate Cancer Susceptibility Locus

Escan: Scan Regulatory Regions for Aggregate Association Testing Using Whole Genome Sequencing Data

scSNV-seq: high-throughput phenotyping of single nucleotide variants by coupled single-cell genotyping and transcriptomics

Association of Prostate Cancer Risk with SNPs in Regions Containing Androgen Receptor Binding Sites Captured by ChIP-On-chip Analyses.

A Sequential Methodology for the Rapid Identification and Characterization of Breast Cancer-Associated Functional SNPs.

Regsnps: a Strategy for Prioritizing Regulatory Single Nucleotide Substitutions.

Functional assessment of human enhancer activities using whole-genome STARR-sequencing

Systematic Enrichment Analysis of Potentially Functional Regions for 103 Prostate Cancer Risk-Associated Loci

SNP2APA: a Database for Evaluating Effects of Genetic Variants on Alternative Polyadenylation in Human Cancers

Chromatin Interactions and Candidate Genes at Ten Prostate Cancer Risk Loci.

Family-Based Association Tests for Genomewide Association Scans

SpRY-mediated screens facilitate functional dissection of non-coding sequences at single-base resolution