Functional Interrogation of Enhancer Connectome Prioritizes Candidate Target Genes at Ovarian Cancer Susceptibility Loci

Identifying causal regulatory variants and their target genes from the majority of non-coding disease-associated genetic loci is the main challenge in post-Genome-Wide Association Studies (GWAS) functional studies. Although chromosome conformation capture (3C) and its derivative technologies have been successfully applied to nominate putative causal genes for non-coding variants, many GWAS target genes have not been identified yet. This study generated a high-resolution contact map from epithelial ovarian cancer (EOC) cells with two H3K27ac-HiChIP libraries and analyzed the underlying gene networks for 15 risk loci identified from the largest EOC GWAS. By combinatory analysis of 4,021 fine-mapped credible variants of EOC GWAS and high-resolution contact map, we obtained 162 target genes that mainly enriched in cancer related pathways. Compared with GTEx eQTL genes in ovarian tissue and annotated proximal genes, 132 HiChIP targets were first identified for EOC causal variants. More than half of the credible variants (CVs) involved interactions that were over 185 kb in distance, indicating that long-range transcriptional regulation is an important mechanism for the function of GWAS variants in EOC. We also found that many HiChIP gene targets showed significantly differential expressions between normal ovarian and EOC tumor samples. We validated one of these targets by manipulating the rs9303542 located region with CRISPR-Cas9 deletion and dCas9-VP64 activation experiments and found altered expression of HOXB7 and HOXB8 at 17q21.32. This study presents a systematic analysis to identify putative target genes for causal variants of EOC, providing an in-depth investigation of the mechanisms of non-coding regulatory variants in the etiology and pathogenesis of ovarian cancer.