GSE35829
Integrative functional genomics identifies an enhancer looping to the SOX9 gene disrupted by the 17q24.3 prostate cancer risk locus
Genome-wide association studies (GWAS) are identifying genetic predisposition to various diseases. The rs1859962 single nucleotide polymorphism (SNP) part of the 17q24.3 locus is a risk factor for prostate cancer (PCa). It defines a 130kb linkage disequilibrium (LD) block that lies in a ~2Mb gene desert area. Despite a role for the proximal SOX9 gene in PCa development, the functional biology driving the risk of this 17q24.3 risk locus is unknown. In the present study, we integrate genome-wide chromatin landscape datasets, namely epigenomes and chromatin openness from diverse cell-types to identify one PCa specific enhancer within the rs1859962 risk LD block. We reveal that this enhancer is part of a 1Mb chromatin loop with the SOX9 gene in PCa cells. The rs8072254 and rs1859961 SNPs part of this LD block map to this enhancer and impose allele-specific gene expression. The variant allele of rs1859961 directly decreases FoxA1 binding while increasing AP-1 binding compared to the reference allele. This latter is key in driving allele-specific gene expression. Together, our results demonstrate the risk associated with the PCa rs1859962 risk LD block is accounted for by multiple genetic variants mapping to a unique enhancer looping to the SOX9 oncogene. Allele-specific recruitment of the transcription factor AP-1 accounts in part for the increased enhancer activity ascribed to this PCa risk LD block. This further demonstrates that an integrative genomics approach can identify the functional biology disrupted by genetic risk-variants. Overall design: Examination of histone modification H3K36me3 in the prostate cancer LNCaP cell line under DHT treatment.