Abstract:While genome-wide association studies and expression quantitative trait loci (eQTL) analysis have made significant progress in identifying noncoding variants associated with prostate cancer risk and bulk tissue transcriptome changes, the regulatory effect of these genetic elements on gene expression remains largely unknown. Recent developments in single-cell sequencing have made it possible to perform ATAC-seq and RNA-seq profiling simultaneously to capture functional associations between chromatin accessibility and gene expression. In this study, we tested our hypothesis that this multiome single-cell approach allows for mapping regulatory elements and their target genes at prostate cancer risk loci. We applied a 10X Multiome ATAC + Gene Expression platform to encapsulate Tn5 transposase-tagged nuclei from multiple prostate cell lines for a total of 65,501 high quality single cells from RWPE1, RWPE2, PrEC, BPH1, DU145, PC3, 22Rv1 and LNCaP cell lines. To address data sparsity commonly seen in the single-cell sequencing, we performed targeted sequencing to enrich sequencing data at prostate cancer risk loci involving 2,730 candidate germline variants and 273 associated genes. Although not increasing the number of captured cells, the targetd multiome data did improve eQTL gene expression abundance by about 20% and chromatin accessibility abundance by about 5%. Based on this multiomic profiling, we further associated RNA expression alterations with chromatin accessibility of germline variants at single cell levels. Cross validation analysis showed high overlaps between the multiome associations and the bulk eQTL findings from GTEx prostate cohort. We found that about 20% of GTEx eQTLs were covered within the significant multiome associations (p-value ≤ 0.05, gene abundance percentage ≥ 5%), and roughly 10% of the multiome associations could be identified by significant GTEx eQTLs. We also analyzed accessible regions with available heterozygous SNP reads and observed more frequent association in genomic regions with allelically accessible variants (p = 0.0055). Among these findings were previously reported regulatory variants including rs60464856-RUVBL1 (multiome p-value = 0.0099 in BPH1) and rs7247241-SPINT2 (multiome p-value = 0.0002- 0.0004 in 22Rv1). We also functionally validated a new regulatory SNP and its target gene rs2474694-VPS53 (multiome p-value = 0.00956 in BPH1 and 0.00625 in DU145) by reporter assay and SILAC proteomics sequencing. Taken together, our data demonstrated the feasibility of the multiome single-cell approach for identifying regulatory SNPs and their regulated genes.

Regsnps-Asb: A Computational Framework for Identifying Allele-Specific Transcription Factor Binding from ATAC-seq Data

atSNP: transcription factor binding affinity testing for regulatory SNP detection

Regsnps: a Strategy for Prioritizing Regulatory Single Nucleotide Substitutions.

Allele-specific transcription factor binding across human brain regions offers mechanistic insight into eQTLs

Systematic identification of regulatory variants associated with cancer risk

A New Approach for Discovering Functional Links Connecting Non-Coding Regulatory Variants to Gene Targets

Identify Regulatory eQTLs by Multiome Sequencing in Prostate Single Cells

Comparative and Functional Genomic Resource for Mechanistic Studies of Human Blood Pressure–Associated Single Nucleotide Polymorphisms

Integrative modeling of eQTLs and cis-regulatory elements suggest mechanisms underlying cell type specificity of eQTLs

Haplotype-aware modeling of cis-regulatory effects highlights the gaps remaining in eQTL data

Allele-specific binding variants causing ChIP-seq peak height of histone modification are not enriched in expression QTL annotations

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

Sasquatch: predicting the impact of regulatory SNPs on transcription factor binding from cell- and tissue-specific DNase footprints

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

Quantitative analysis of cis-regulatory elements in transcription with KAS-ATAC-seq

Single nucleotide polymorphisms that modulate microRNA regulation of gene expression in tumors

Systematic Analysis of the Effects of Genetic Variants on Chromatin Accessibility to Decipher Functional Variants in Non-Coding Regions

Prediction of Functional Regulatory SNPs in Monogenic and Complex Disease.

Distant regulatory effects of genetic variation in multiple human tissues

A Novel Method for in Silico Identification of Regulatory SNPs in Human Genome.

Systematic Analysis of Binding of Transcription Factors to Noncoding Variants.