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.

Powerful mapping of cis-genetic effects on gene expression across diverse populations reveals novel disease-critical genes

Leveraging trans-ethnic genetic risk scores to improve association power for complex traits in underrepresented populations

Inferring Gene-Disease Association by an Integrative Analysis of eQTL Genome-Wide Association Study and Protein-Protein Interaction Data.

Supplementary Material For: Inferring Gene-Disease Association by an Integrative Analysis of Eqtl Genome-Wide Association Study and Protein-Protein Interaction Data

MGAS: a Powerful Tool for Multivariate Gene-Based Genome-Wide Association Analysis

Leveraging functional genomic annotations and genome coverage to improve polygenic prediction of complex traits within and between ancestries

Quantifying Portable Genetic Effects and Improving Cross-Ancestry Genetic Prediction with GWAS Summary Statistics

Identify Regulatory eQTLs by Multiome Sequencing in Prostate Single Cells

Powerful and robust inference of complex phenotypes' causal genes with dependent expression quantitative loci by a median-based Mendelian randomization.

Optimization of Multi-Ancestry Polygenic Risk Score Disease Prediction Models

Leveraging Functional Annotation to Identify Genes Associated with Complex Diseases.

Cross-population enhancement of PrediXcan predictions with a gnomAD-based east Asian reference framework

Cis-SNPs Set Testing and PrediXcan Analysis for Gene Expression Data using Linear Mixed Models

-eQTL mapping in gene sets identifies network effects of genetic variants

Local genetic effects on gene expression across 44 human tissues

Exploring the phenotypic consequences of tissue specific gene expression variation inferred from GWAS summary statistics

Evaluating Marginal Genetic Correlation of Associated Loci for Complex Diseases and Traits Between European and East Asian Populations.

Integrating Predicted Transcriptome From Multiple Tissues Improves Association Detection

A covariance-enhanced approach to multi-tissue joint eQTL mapping with application to transcriptome-wide association studies

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

Improved multi-ancestry fine-mapping identifies -regulatory variants underlying molecular traits and disease risk