Identification of candidate genes for prostate cancer-risk SNPs utilizing a normal prostate tissue eQTL data set

Molecular analysis of archival diagnostic prostate cancer biopsies identifies genomic similarities in cases with progression post-radiotherapy, and those with de novo metastatic disease

BACKGROUND

It is important to identify molecular features that improve prostate cancer (PCa) risk stratification before radical treatment with curative intent. Molecular analysis of historical diagnostic formalin-fixed paraffin-embedded (FFPE) prostate biopsies from cohorts with post-radiotherapy (RT) long-term clinical follow-up has been limited. Utilizing parallel sequencing modalities, we performed a proof-of-principle sequencing analysis of historical diagnostic FFPE prostate biopsies. We compared patients with (i) stable PCa (sPCa) postprimary or salvage RT, (ii) progressing PCa (pPCa) post-RT, and (iii) de novo metastatic PCa (mPCa).

METHODS

A cohort of 19 patients with diagnostic prostate biopsies (n = 6 sPCa, n = 5 pPCa, n = 8 mPCa) and mean 4 years 10 months follow-up (diagnosed 2009-2016) underwent nucleic acid extraction from demarcated malignancy. Samples underwent 3'RNA sequencing (3'RNAseq) (n = 19), nanoString analysis (n = 12), and Illumina 850k methylation (n = 8) sequencing. Bioinformatic analysis was performed to coherently identify differentially expressed genes and methylated genomic regions (MGRs).

RESULTS

Eighteen of 19 samples provided useable 3'RNAseq data. Principal component analysis (PCA) demonstrated similar expression profiles between pPCa and mPCa cases, versus sPCa. Coherently differentially methylated probes between these groups identified ~600 differentially MGRs. The top 50 genes with increased expression in pPCa patients were associated with reduced progression-free survival post-RT (p < 0.0001) in an external cohort.

CONCLUSIONS

3'RNAseq, nanoString and 850k-methylation analyses are each achievable from historical FFPE diagnostic pretreatment prostate biopsies, unlocking the potential to utilize large cohorts of historic clinical samples. Profiling similarities between individuals with pPCa and mPCa suggests biological similarities and historical radiological staging limitations, which warrant further investigation.

Identify Regulatory eQTLs by Multiome Sequencing in Prostate Single Cells

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 targeted 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.

Integrative cross-omics and cross-context analysis elucidates molecular links underlying genetic effects on complex traits

Genetic effects on functionally related 'omic' traits often co-occur in relevant cellular contexts, such as tissues. Motivated by the multi-tissue methylation quantitative trait loci (mQTLs) and expression QTLs (eQTLs) analysis, we propose X-ING (Cross-INtegrative Genomics) for cross-omics and cross-context integrative analysis. X-ING takes as input multiple matrices of association statistics, each obtained from different omics data types across multiple cellular contexts. It models the latent binary association status of each statistic, captures the major association patterns among omics data types and contexts, and outputs the posterior mean and probability for each input statistic. X-ING enables the integration of effects from different omics data with varying effect distributions. In the multi-tissue cis-association analysis, X-ING shows improved detection and replication of mQTLs by integrating eQTL maps. In the trans-association analysis, X-ING reveals an enrichment of trans-associations in many disease/trait-relevant tissues.

Characterizing prostate cancer risk through multi-ancestry genome-wide discovery of 187 novel risk variants

The transferability and clinical value of genetic risk scores (GRSs) across populations remain limited due to an imbalance in genetic studies across ancestrally diverse populations. Here we conducted a multi-ancestry genome-wide association study of 156,319 prostate cancer cases and 788,443 controls of European, African, Asian and Hispanic men, reflecting a 57% increase in the number of non-European cases over previous prostate cancer genome-wide association studies. We identified 187 novel risk variants for prostate cancer, increasing the total number of risk variants to 451. An externally replicated multi-ancestry GRS was associated with risk that ranged from 1.8 (per standard deviation) in African ancestry men to 2.2 in European ancestry men. The GRS was associated with a greater risk of aggressive versus non-aggressive disease in men of African ancestry (P = 0.03). Our study presents novel prostate cancer susceptibility loci and a GRS with effective risk stratification across ancestry groups.

Epigenome editing in cancer: Advances and challenges for potential therapeutic options

Cancers are diseases caused by genetic and non-genetic environmental factors. Epigenetic alterations, some attributed to non-genetic factors, can lead to cancer development. Epigenetic changes can occur in tumor suppressors or oncogenes, or they may contribute to global cell state changes, making cells abnormal. Recent advances in gene editing technology show potential for cancer treatment. Herein, we will discuss our current knowledge of epigenetic alterations occurring in cancer and epigenetic editing technologies that can be applied to developing therapeutic options.

Combined CRISPRi and proteomics screening reveal a cohesin-CTCF-bound allele contributing to increased expression of RUVBL1 and prostate cancer progression

Genome-wide association studies along with expression quantitative trait locus (eQTL) mapping have identified hundreds of single-nucleotide polymorphisms (SNPs) and their target genes in prostate cancer (PCa), yet functional characterization of these risk loci remains challenging. To screen for potential regulatory SNPs, we designed a CRISPRi library containing 9,133 guide RNAs (gRNAs) to cover 2,166 candidate SNP loci implicated in PCa and identified 117 SNPs that could regulate 90 genes for PCa cell growth advantage. Among these, rs60464856 was covered by multiple gRNAs significantly depleted in screening (FDR < 0.05). Pooled SNP association analysis in the PRACTICAL and FinnGen cohorts showed significantly higher PCa risk for the rs60464856 G allele (p value = 1.2 × 10-16 and 3.2 × 10-7, respectively). Subsequent eQTL analysis revealed that the G allele is associated with increased RUVBL1 expression in multiple datasets. Further CRISPRi and xCas9 base editing confirmed that the rs60464856 G allele leads to elevated RUVBL1 expression. Furthermore, SILAC-based proteomic analysis demonstrated allelic binding of cohesin subunits at the rs60464856 region, where the HiC dataset showed consistent chromatin interactions in prostate cell lines. RUVBL1 depletion inhibited PCa cell proliferation and tumor growth in a xenograft mouse model. Gene-set enrichment analysis suggested an association of RUVBL1 expression with cell-cycle-related pathways. Increased expression of RUVBL1 and activation of cell-cycle pathways were correlated with poor PCa survival in TCGA datasets. Our CRISPRi screening prioritized about one hundred regulatory SNPs essential for prostate cell proliferation. In combination with proteomics and functional studies, we characterized the mechanistic role of rs60464856 and RUVBL1 in PCa progression.

Combined CRISPRi and proteomics screening reveal a cohesin-CTCF-bound allele contributing to increased expression of RUVBL1 and prostate cancer progression

Genome-wide association studies along with expression quantitative trait loci (eQTL) mapping have identified hundreds of single nucleotide polymorphisms (SNPs) and their target genes in prostate cancer (PCa), yet functional characterization of these risk loci remains challenging. To screen for potential regulatory SNPs, we designed a CRISPRi library containing 9133 guide RNAs (gRNAs) to target 2,166 candidate SNP sites implicated in PCa and identified 117 SNPs that could regulate 90 genes for PCa cell growth advantage. Among these, rs60464856 was covered by multiple gRNAs significantly depleted in the screening (FDR<0.05). Pooled SNP association analysis in the PRACTICAL and FinnGen cohorts showed significantly higher PCa risk for the rs60464856 G allele (pvalue=1.2E-16 and 3.2E-7). Subsequent eQTL analysis revealed that the G allele is associated with increased RUVBL1 expression in multiple datasets. Further CRISPRi and xCas9 base editing proved the rs60464856 G allele leading to an elevated RUVBL1 expression. Furthermore, SILAC-based proteomic analysis demonstrated allelic binding of cohesin subunits at the rs60464856 region, where HiC dataset showed consistent chromatin interactions in prostate cell lines. RUVBL1 depletion inhibited PCa cell proliferation and tumor growth in xenograft mouse model. Gene set enrichment analysis suggested an association of RUVBL1 expression with cell-cycle-related pathways. An increased expression of RUVBL1 and activations of cell-cycle pathways were correlated with poor PCa survival in TCGA datasets. Together, our CRISPRi screening prioritized about one hundred regulatory SNPs essential for prostate cell proliferation. In combination with proteomics and functional studies, we characterized the mechanistic role of rs60464856 and RUVBL1 in PCa progression.

Association between interleuleukin-1β polymorphism (rs16944) and biomarkers levels in Iraqi patients with prostate cancer

OBJECTIVE

Prostate cancer (PCa) is the second-leading cause of mortality in men and the most commonly diagnosed non-cutaneous male malignancy. Host genetic factors, and inflammation-induced cytokines, play a key role in prostate oncogenesis. Single Nucleotide Polymorphisms (SNP) in cytokine genes were suggested to increase the susceptibility for PCa development and progression. This study aimed to investigate the association between the SNP (rs16944) in the interleukin-1 β (IL-1β) gene and the serum levels of Prostate Specific Antigen (PSA) Prolactine (PRL), testosterone, and IL-1β in Iraqi PCa patients versus healthy controls.

METHODS

Taqman Real Time-PCR, was performed to investigate the IL-1β (rs16944) polymorphism in 100 Iraqi PCa patients and 50 age-matched healthy controls in a case-control study. Serum levels of PSA, PRL, and testosterone were determined by ELISA and FIA, and associated with the IL-1β serum level as well as with the SNP (rs 16944). The association between the clinico-pathological parameters and the genotype distribution of PCa patients was also studied.

RESULTS

There level of IL-1β was significant increased in the serum of PCa patients compared to controls (P = 8.19 10^-7). Serum levels for other biomarkers such as PSA, PRL and testosterone were also significantly elevated in patients compared to controls (P < 0.0001). No differences were seen for genotype and allele distribution between PCa patients and controls. Nevertheless, in the group of controls, we found that 36% carried the GG genotype against only 26% in the patients group.This suggests that this could be a protective genotype (OR 0.62, P = 0.254). In addition, we found that the GA genotype is slightly more frequent in patients as compared to controls (OR 1.22, P = 0.605). Interestingly, serum levels of IL-1β, PSA, PRL and testosterone were significantly higher in PCa patients carrying the GA genotype, and the GA and AA genotypes are strongly associated with the aggressive behavior of the disease such as advanced TNM, and high Gleason score.

CONCLUSION

Our data suggest that both serum IL-1β level and IL-1β SNP (rs16944) may be considered as candidate biomarkers for PCa. Moreover, the GA, and AA genotypes carriers along with high sera levels of IL-1β, PSA and PRL, have an increased risk for PCa with aggressive behavior in Iraqi men.

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

The vast majority of disease-associated single nucleotide polymorphisms (SNP) identified from genome-wide association studies (GWAS) are localized in non-coding regions. A significant fraction of these variants impact transcription factors binding to enhancer elements and alter gene expression. To functionally interrogate the activity of such variants we developed snpSTARRseq, a high-throughput experimental method that can interrogate the functional impact of hundreds to thousands of non-coding variants on enhancer activity. snpSTARRseq dramatically improves signal-to-noise by utilizing a novel sequencing and bioinformatic approach that increases both insert size and the number of variants tested per loci. Using this strategy, we interrogated known prostate cancer (PCa) risk-associated loci and demonstrated that 35% of them harbor SNPs that significantly altered enhancer activity. Combining these results with chromosomal looping data we could identify interacting genes and provide a mechanism of action for 20 PCa GWAS risk regions. When benchmarked to orthogonal methods, snpSTARRseq showed a strong correlation with in vivo experimental allelic-imbalance studies whereas there was no correlation with predictive in silico approaches. Overall, snpSTARRseq provides an integrated experimental and computational framework to functionally test non-coding genetic variants.

Prostate Cancer Transcriptomic Regulation by the Interplay of Germline Risk Alleles, Somatic Mutations, and 3D Genomic Architecture

Prostate cancer is one of the most heritable human cancers. Genome-wide association studies have identified at least 185 prostate cancer germline risk alleles, most noncoding. We used integrative three-dimensional (3D) spatial genomics to identify the chromatin interaction targets of 45 prostate cancer risk alleles, 31 of which were associated with the transcriptional regulation of target genes in 565 localized prostate tumors. To supplement these 31, we verified transcriptional targets for 56 additional risk alleles using linear proximity and linkage disequilibrium analysis in localized prostate tumors. Some individual risk alleles influenced multiple target genes; others specifically influenced only distal genes while leaving proximal ones unaffected. Several risk alleles exhibited widespread germline-somatic interactions in transcriptional regulation, having different effects in tumors with loss of PTEN or RB1 relative to those without. These data clarify functional prostate cancer risk alleles in large linkage blocks and outline a strategy to model multidimensional transcriptional regulation.

SIGNIFICANCE

Many prostate cancer germline risk alleles are enriched in the noncoding regions of the genome and are hypothesized to regulate transcription. We present a 3D genomics framework to unravel risk SNP function and describe the widespread germline-somatic interplay in transcription control. This article is highlighted in the In This Issue feature, p. 2711.

Extensive germline-somatic interplay contributes to prostate cancer progression through HNF1B co-option of TMPRSS2-ERG

Genome-wide association studies have identified 270 loci conferring risk for prostate cancer (PCa), yet the underlying biology and clinical impact remain to be investigated. Here we observe an enrichment of transcription factor genes including HNF1B within PCa risk-associated regions. While focused on the 17q12/HNF1B locus, we find a strong eQTL for HNF1B and multiple potential causal variants involved in the regulation of HNF1B expression in PCa. An unbiased genome-wide co-expression analysis reveals PCa-specific somatic TMPRSS2-ERG fusion as a transcriptional mediator of this locus and the HNF1B eQTL signal is ERG fusion status dependent. We investigate the role of HNF1B and find its involvement in several pathways related to cell cycle progression and PCa severity. Furthermore, HNF1B interacts with TMPRSS2-ERG to co-occupy large proportion of genomic regions with a remarkable enrichment of additional PCa risk alleles. We finally show that HNF1B co-opts ERG fusion to mediate mechanistic and biological effects of the PCa risk-associated locus 17p13.3/VPS53/FAM57A/GEMIN4. Taken together, we report an extensive germline-somatic interaction between TMPRSS2-ERG fusion and genetic variations underpinning PCa risk association and progression.

Analysis of the potential role of photocurable hydrogel in patient-derived glioblastoma organoid culture through RNA sequencing

Patient-derived glioblastoma organoid (GBO) growth in hydrogels recapitulates key features of parental tumors, making GBOs a useful tool for fundamental research on cancer biology and offer deeper insight into the development of innovative therapeutic strategies for cancer treatment. Matrigel as a natural hydrogel has been widely used for 3D culture in most tumor organoid studies, but the volatility in its biochemical and biophysical properties makes it difficult to be further applied in GBO cultures. Thus, several kinds of biomimetic hydrogels from synthetic or biological polymers have been developed for tumor organoid growth. Here, we innovatively utilize a photocurable hydrogel-based biomimetic instructive system containing gelatin methacryloyl (GelMA) mixed with a hyaluronic acid (HA) hydrogel as a scaffold for generating GBOs. Furthermore, we evaluated the GBO biological properties at the transcriptome level, which showed that GBOs cultured with this hydrogel retain the expression profile of key neurodevelopmental markers, driving mutations and alternative splicing of parental tumors. Notably, GBOs cultured with the photocurable hydrogel may provide a platform for precision cancer medicine, bridging the gap between basic research and clinical application. Although significant challenges remain, biomimetic hydrogels can provide an exceptional window for the construction of tumor organoids to ensure the accuracy of the research and clinical data.

A microRNA Transcriptome-wide Association Study of Prostate Cancer Risk

Large genome-wide association studies have identified hundreds of single-nucleotide polymorphisms associated with increased risk of prostate cancer (PrCa), and many of these risk loci is presumed to confer regulatory effects on gene expression. While eQTL studies of long RNAs has yielded many potential risk genes, the relationship between PrCa risk genetics and microRNA expression dysregulation is understudied. We performed an microRNA transcriptome-wide association study of PrCa risk using small RNA sequencing and genome-wide genotyping data from N = 441 normal prostate epithelium tissue samples along with N = 411 prostate adenocarcinoma tumor samples from the Cancer Genome Atlas (TCGA). Genetically regulated expression prediction models were trained for all expressed microRNAs using the FUSION TWAS software. TWAS for PrCa risk was performed with both sets of models using single-SNP summary statistics from the recent PRACTICAL consortium PrCa case-control OncoArray GWAS meta-analysis. A total of 613 and 571 distinct expressed microRNAs were identified in the normal and tumor tissue datasets, respectively (overlap: 480). Among these, 79 (13%) normal tissue microRNAs demonstrated significant cis-heritability (median cis-h2 = 0.15, range: 0.03–0.79) for model training. Similar results were obtained from TCGA tumor samples, with 48 (9%) microRNA expression models successfully trained (median cis-h2 = 0.14, range: 0.06–0.60). Using normal tissue models, we identified two significant TWAS microRNA associations with PrCa risk: over-expression of mir-941 family microRNAs (PTWAS = 2.9E-04) and reduced expression of miR-3617-5p (PTWAS = 1.0E-03). The TCGA tumor TWAS also identified a significant association with miR-941 overexpression (PTWAS = 9.7E-04). Subsequent finemapping of the TWAS results using a multi-tissue database indicated limited evidence of causal status for each microRNA with PrCa risk (posterior inclusion probabilities &lt;0.05). Future work will examine downstream regulatory effects of microRNA dysregulation as well as microRNA-mediated risk mechanisms via competing endogenous RNA relationships.

Identification of novel prostate cancer genes in patients stratified by Gleason classification: role of antitumoral genes

Prostate cancer (PCa) is a tumor with a great heterogeneity, both at a molecular and clinical level. Despite its global good prognosis, cases can vary from indolent to lethal metastatic and scientific efforts are aimed to discern those with worse outcomes. Current prognostic markers, as Gleason score, fall short when it comes to distinguishing these cases. Identification of new early biomarkers to enable a better PCa distinction and classification remains a challenge. In order to identify new genes implicated in PCa progression we conducted several differential gene expression analyses over paired samples comparing primary PCa tissue against healthy prostatic tissue of PCa patients. The results obtained show that this approach is a serious alternative to overcome patient heterogeneity. We were able to identify 250 genes whose expression varies along with tissue differentiation—healthy to tumor tissue, 161 of these genes are described here for the first time to be related to PCa. The further manual curation of these genes allowed to annotate 39 genes with antitumoral activity, 22 of them described for the first time to be related to PCa proliferation and metastasis. These findings could be replicated in different cohorts for most genes. Results obtained considering paired differential expression, functional annotation and replication results point to: CGREF1, UNC5A, C16orf74, LGR6, IGSF1, QPRT and CA14 as possible new early markers in PCa. These genes may prevent the progression of the disease and their expression should be studied in patients with different outcomes.

Interpreting Functional Impact of Genetic Variations by Network QTL for Genotype–Phenotype Association Study

An enormous challenge in the post-genome era is to annotate and resolve the consequences of genetic variation on diverse phenotypes. The genome-wide association study (GWAS) is a well-known method to identify potential genetic loci for complex traits from huge genetic variations, following which it is crucial to identify expression quantitative trait loci (eQTL). However, the conventional eQTL methods usually disregard the systematical role of single-nucleotide polymorphisms (SNPs) or genes, thereby overlooking many network-associated phenotypic determinates. Such a problem motivates us to recognize the network-based quantitative trait loci (QTL), i.e., network QTL (nQTL), which is to detect the cascade association as genotype → network → phenotype rather than conventional genotype → expression → phenotype in eQTL. Specifically, we develop the nQTL framework on the theory and approach of single-sample networks, which can identify not only network traits (e.g., the gene subnetwork associated with genotype) for analyzing complex biological processes but also network signatures (e.g., the interactive gene biomarker candidates screened from network traits) for characterizing targeted phenotype and corresponding subtypes. Our results show that the nQTL framework can efficiently capture associations between SNPs and network traits (i.e., edge traits) in various simulated data scenarios, compared with traditional eQTL methods. Furthermore, we have carried out nQTL analysis on diverse biological and biomedical datasets. Our analysis is effective in detecting network traits for various biological problems and can discover many network signatures for discriminating phenotypes, which can help interpret the influence of nQTL on disease subtyping, disease prognosis, drug response, and pathogen factor association. Particularly, in contrast to the conventional approaches, the nQTL framework could also identify many network traits from human bulk expression data, validated by matched single-cell RNA-seq data in an independent or unsupervised manner. All these results strongly support that nQTL and its detection framework can simultaneously explore the global genotype–network–phenotype associations and the underlying network traits or network signatures with functional impact and importance.

A transcriptome-wide association study identifies novel candidate susceptibility genes for prostate cancer risk

A large proportion of heritability for prostate cancer risk remains unknown. Transcriptome-wide association study combined with validation comparing overall levels will help to identify candidate genes potentially playing a role in prostate cancer development. Using data from the Genotype-Tissue Expression Project, we built genetic models to predict normal prostate tissue gene expression using the statistical framework PrediXcan, a modified version of the unified test for molecular signatures and Joint-Tissue Imputation. We applied these prediction models to the genetic data of 79 194 prostate cancer cases and 61 112 controls to investigate the associations of genetically determined gene expression with prostate cancer risk. Focusing on associated genes, we compared their expression in prostate tumor vs normal prostate tissue, compared methylation of CpG sites located at these loci in prostate tumor vs normal tissue, and assessed the correlations between the differentiated genes' expression and the methylation of corresponding CpG sites, by analyzing The Cancer Genome Atlas (TCGA) data. We identified 573 genes showing an association with prostate cancer risk at a false discovery rate (FDR) ≤ 0.05, including 451 novel genes and 122 previously reported genes. Of the 573 genes, 152 showed differential expression in prostate tumor vs normal tissue samples. At loci of 57 genes, 151 CpG sites showed differential methylation in prostate tumor vs normal tissue samples. Of these, 20 CpG sites were correlated with expression of 11 corresponding genes. In this TWAS, we identified novel candidate susceptibility genes for prostate cancer risk, providing new insights into prostate cancer genetics and biology.

H3K27ac HiChIP in prostate cell lines identifies risk genes for prostate cancer susceptibility

Genome-wide association studies (GWASs) have identified more than 200 prostate cancer (PrCa) risk regions, which provide potential insights into causal mechanisms. Multiple lines of evidence show that a significant proportion of PrCa risk can be explained by germline causal variants that dysregulate nearby target genes in prostate-relevant tissues, thus altering disease risk. The traditional approach to explore this hypothesis has been correlating GWAS variants with steady-state transcript levels, referred to as expression quantitative trait loci (eQTLs). In this work, we assess the utility of chromosome conformation capture (3C) coupled with immunoprecipitation (HiChIP) to identify target genes for PrCa GWAS risk loci. We find that interactome data confirm previously reported PrCa target genes identified through GWAS/eQTL overlap (e.g., MLPH). Interestingly, HiChIP identifies links between PrCa GWAS variants and genes well-known to play a role in prostate cancer biology (e.g., AR) that are not detected by eQTL-based methods. HiChIP predicted enhancer elements at the AR and NKX3-1 prostate cancer risk loci, and both were experimentally confirmed to regulate expression of the corresponding genes through CRISPR interference (CRISPRi) perturbation in LNCaP cells. Our results demonstrate that looping data harbor additional information beyond eQTLs and expand the number of PrCa GWAS loci that can be linked to candidate susceptibility genes.

Mechanistic insights into genetic susceptibility to prostate cancer

Prostate cancer (PCa) is the second most common cancer in men and is a highly heritable disease that affects millions of individuals worldwide. Genome-wide association studies have to date discovered nearly 270 genetic loci harboring hundreds of single nucleotide polymorphisms (SNPs) that are associated with PCa susceptibility. In contrast, the functional characterization of the mechanisms underlying PCa risk association is still growing. Given that PCa risk-associated SNPs are highly enriched in noncoding cis-regulatory genomic regions, accumulating evidence suggests a widespread modulation of transcription factor chromatin binding and allelic enhancer activity by these noncoding SNPs, thereby dysregulating gene expression. Emerging studies have shown that a proportion of noncoding variants can modulate the formation of transcription factor complexes at enhancers and CTCF-mediated 3D genome architecture. Interestingly, DNA methylation-regulated CTCF binding could orchestrate a long-range chromatin interaction between PCa risk enhancer and causative genes. Additionally, one-causal-variant-two-risk genes or multiple-risk-variant-multiple-genes are prevalent in some PCa risk-associated loci. In this review, we will discuss the current understanding of the general principles of SNP-mediated gene regulation, experimental advances, and functional evidence supporting the mechanistic roles of several PCa genetic loci with potential clinical impact on disease prevention and treatment.

Meta-Analyses of Splicing and Expression Quantitative Trait Loci Identified Susceptibility Genes of Glioma

Background

The functions of most glioma risk alleles are unknown. Very few studies had evaluated expression quantitative trait loci (eQTL), and insights of susceptibility genes were limited due to scarcity of available brain tissues. Moreover, no prior study had examined the effect of glioma risk alleles on alternative RNA splicing.

Objective

This study explored splicing quantitative trait loci (sQTL) as molecular QTL and improved the power of QTL mapping through meta-analyses of both cis eQTL and sQTL.

Methods

We first evaluated eQTLs and sQTLs of the CommonMind Consortium (CMC) and Genotype-Tissue Expression Project (GTEx) using genotyping, or whole-genome sequencing and RNA-seq data. Alternative splicing events were characterized using an annotation-free method that detected intron excision events. Then, we conducted meta-analyses by pooling the eQTL and sQTL results of CMC and GTEx using the inverse variance-weighted model. Afterward, we integrated QTL meta-analysis results (Q < 0.05) with the Glioma International Case Control Study (GICC) GWAS meta-analysis (case:12,496, control:18,190), using a summary statistics-based mendelian randomization (SMR) method.

Results

Between CMC and GTEx, we combined the QTL data of 354 unique individuals of European ancestry. SMR analyses revealed 15 eQTLs in 11 loci and 32 sQTLs in 9 loci relevant to glioma risk. Two loci only harbored sQTLs (1q44 and 16p13.3). In seven loci, both eQTL and sQTL coexisted (2q33.3, 7p11.2, 11q23.3 15q24.2, 16p12.1, 20q13.33, and 22q13.1), but the target genes were different for five of these seven loci. Three eQTL loci (9p21.3, 20q13.33, and 22q13.1) and 4 sQTL loci (11q23.3, 16p13.3, 16q12.1, and 20q13.33) harbored multiple target genes. Eight target genes of sQTLs (C2orf80, SEC61G, TMEM25, PHLDB1, RP11-161M6.2, HEATR3, RTEL1-TNFRSF6B, and LIME1) had multiple alternatively spliced transcripts.

Conclusion

Our study revealed that the regulation of transcriptome by glioma risk alleles is complex, with the potential for eQTL and sQTL jointly affecting gliomagenesis in risk loci. QTLs of many loci involved multiple target genes, some of which were specific to alternative splicing. Therefore, quantitative trait loci that evaluate only total gene expression will miss many important target genes.

Single nucleotide polymorphisms (SNPs) in prostate cancer: its implications in diagnostics and therapeutics

Prostate cancer is one of the most frequently diagnosed malignancies in developed countries and approximately 248,530 new cases of prostate cancer are likely to be diagnosed in the United States in 2021. During the late 1990s and 2000s, the prostate cancer-related death rate has decreased by 4% per year on average because of advancements in prostate-specific antigen (PSA) testing. However, the non-specificity of PSA to distinguish between benign and malignant forms of cancer is a major concern in the management of prostate cancer. Despite other risk factors in the pathogenesis of prostate cancer, recent advancement in molecular genetics suggests that genetic heredity plays a crucial role in prostate carcinogenesis. Approximately, 60% of heritability and more than 100 well-recognized single-nucleotide-polymorphisms (SNPs) have been found to be associated with prostate cancer and constitute a major risk factor in the development of prostate cancer. Recent findings revealed that a low to moderate effect on the progression of prostate cancer of individual SNPs was observed compared to a strong progressive effect when SNPs were in combination. Here, in this review, we made an attempt to critically analyze the role of SNPs and associated genes in the development of prostate cancer and their implications in diagnostics and therapeutics. A better understanding of the role of SNPs in prostate cancer susceptibility may improve risk prediction, enhance fine-mapping, and furnish new insights into the underlying pathophysiology of prostate cancer.

A Large-Scale Association Study Detects Novel Rare Variants, Risk Genes, Functional Elements, and Polygenic Architecture of Prostate Cancer Susceptibility

To identify rare variants associated with prostate cancer susceptibility and better characterize the mechanisms and cumulative disease risk associated with common risk variants, we conducted an integrated study of prostate cancer genetic etiology in two cohorts using custom genotyping microarrays, large imputation reference panels, and functional annotation approaches. Specifically, 11,984 men (6,196 prostate cancer cases and 5,788 controls) of European ancestry from Northern California Kaiser Permanente were genotyped and meta-analyzed with 196,269 men of European ancestry (7,917 prostate cancer cases and 188,352 controls) from the UK Biobank. Three novel loci, including two rare variants (European ancestry minor allele frequency < 0.01, at 3p21.31 and 8p12), were significant genome wide in a meta-analysis. Gene-based rare variant tests implicated a known prostate cancer gene (HOXB13), as well as a novel candidate gene (ILDR1), which encodes a receptor highly expressed in prostate tissue and is related to the B7/CD28 family of T-cell immune checkpoint markers. Haplotypic patterns of long-range linkage disequilibrium were observed for rare genetic variants at HOXB13 and other loci, reflecting their evolutionary history. In addition, a polygenic risk score (PRS) of 188 prostate cancer variants was strongly associated with risk (90th vs. 40th-60th percentile OR = 2.62, P = 2.55 × 10-191). Many of the 188 variants exhibited functional signatures of gene expression regulation or transcription factor binding, including a 6-fold difference in log-probability of androgen receptor binding at the variant rs2680708 (17q22). Rare variant and PRS associations, with concomitant functional interpretation of risk mechanisms, can help clarify the full genetic architecture of prostate cancer and other complex traits. SIGNIFICANCE: This study maps the biological relationships between diverse risk factors for prostate cancer, integrating different functional datasets to interpret and model genome-wide data from over 200,000 men with and without prostate cancer.See related commentary by Lachance, p. 1637.

Genome-wide association analysis reveals regulation of at-risk loci by DNA methylation in prostate cancer

Epigenetic changes are potentially important for the ontogeny and progression of tumors but are not usually studied because of the complexity of analyzing transcript regulation resulting from epigenetic alterations. Prostate cancer (PCa) is characterized by variable clinical manifestations and frequently unpredictable outcomes. We performed an expression quantitative trait loci (eQTL) analysis to identify the genomic regions that regulate gene expression in PCa and identified a relationship between DNA methylation and clinical information. Using multi-level information published in The Cancer Genome Atlas, we performed eQTL-based analyses on DNA methylation and gene expression. To better interpret these data, we correlated loci and clinical indexes to identify the important loci for both PCa development and progression. Our data demonstrated that although only a small proportion of genes are regulated via DNA methylation in PCa, these genes are enriched in important cancer-related groups. In addition, single nucleotide polymorphism analysis identified the locations of CpG sites and genes within at-risk loci, including the 19q13.2–q13.43 and 16q22.2–q23.1 loci. Further, an epigenetic association study of clinical indexes detected risk loci and pyrosequencing for site validation. Although DNA methylation-regulated genes across PCa samples are a small proportion, the associated genes play important roles in PCa carcinogenesis.

CRISPRi screens reveal a DNA methylation-mediated 3D genome dependent causal mechanism in prostate cancer

Prostate cancer (PCa) risk-associated SNPs are enriched in noncoding cis-regulatory elements (rCREs), yet their modi operandi and clinical impact remain elusive. Here, we perform CRISPRi screens of 260 rCREs in PCa cell lines. We find that rCREs harboring high risk SNPs are more essential for cell proliferation and H3K27ac occupancy is a strong indicator of essentiality. We also show that cell-line-specific essential rCREs are enriched in the 8q24.21 region, with the rs11986220-containing rCRE regulating MYC and PVT1 expression, cell proliferation and tumorigenesis in a cell-line-specific manner, depending on DNA methylation-orchestrated occupancy of a CTCF binding site in between this rCRE and the MYC promoter. We demonstrate that CTCF deposition at this site as measured by DNA methylation level is highly variable in prostate specimens, and observe the MYC eQTL in the 8q24.21 locus in individuals with low CTCF binding. Together our findings highlight a causal mechanism synergistically driven by a risk SNP and DNA methylation-mediated 3D genome architecture, advocating for the integration of genetics and epigenetics in assessing risks conferred by genetic predispositions.

Trans-ancestry genome-wide association meta-analysis of prostate cancer identifies new susceptibility loci and informs genetic risk prediction

Prostate cancer is a highly heritable disease with large disparities in incidence rates across ancestry populations. We conducted a multiancestry meta-analysis of prostate cancer genome-wide association studies (107,247 cases and 127,006 controls) and identified 86 new genetic risk variants independently associated with prostate cancer risk, bringing the total to 269 known risk variants. The top genetic risk score (GRS) decile was associated with odds ratios that ranged from 5.06 (95% confidence interval (CI), 4.84-5.29) for men of European ancestry to 3.74 (95% CI, 3.36-4.17) for men of African ancestry. Men of African ancestry were estimated to have a mean GRS that was 2.18-times higher (95% CI, 2.14-2.22), and men of East Asian ancestry 0.73-times lower (95% CI, 0.71-0.76), than men of European ancestry. These findings support the role of germline variation contributing to population differences in prostate cancer risk, with the GRS offering an approach for personalized risk prediction.

An integrative multi-omics analysis to identify candidate DNA methylation biomarkers related to prostate cancer risk

It remains elusive whether some of the associations identified in genome-wide association studies of prostate cancer (PrCa) may be due to regulatory effects of genetic variants on CpG sites, which may further influence expression of PrCa target genes. To search for CpG sites associated with PrCa risk, here we establish genetic models to predict methylation (N = 1,595) and conduct association analyses with PrCa risk (79,194 cases and 61,112 controls). We identify 759 CpG sites showing an association, including 15 located at novel loci. Among those 759 CpG sites, methylation of 42 is associated with expression of 28 adjacent genes. Among 22 genes, 18 show an association with PrCa risk. Overall, 25 CpG sites show consistent association directions for the methylation-gene expression-PrCa pathway. We identify DNA methylation biomarkers associated with PrCa, and our findings suggest that specific CpG sites may influence PrCa via regulating expression of candidate PrCa target genes.

Prostate cancer risk SNP rs10993994 is a trans-eQTL for SNHG11 mediated through MSMB

How genome-wide association studies-identified single-nucleotide polymorphisms (SNPs) affect remote genes remains unknown. Expression quantitative trait locus (eQTL) association meta-analysis on 496 prostate tumor and 602 normal prostate samples with 117 SNPs revealed novel cis-eQTLs and trans-eQTLs. Mediation testing and colocalization analysis demonstrate that MSMB is a cis-acting mediator for SNHG11 (P < 0.01). Removing rs10993994 in LNCaP cell lines by CRISPR/Cas9 editing shows that the C-allele corresponds with an over 100-fold increase in MSMB expression and 5-fold increase in SNHG11 compared with the T-allele. Colocalization analysis confirmed that the same set of SNPs associated with MSMB expression is associated with SNHG11 expression (posterior probability of shared variants is 66.6% in tumor and 91.4% in benign). These analyses further demonstrate variants driving MSMB expression differ in tumor and normal, suggesting regulatory network rewiring during tumorigenesis.

Pathway Analysis of Genes Identified through Post-GWAS to Underpin Prostate Cancer Aetiology

Understanding the functional role of risk regions identified by genome-wide association studies (GWAS) has made considerable recent progress and is referred to as the post-GWAS era. Annotation of functional variants to the genes, including cis or trans and understanding their biological pathway/gene network enrichments, is expected to give rich dividends by elucidating the mechanisms underlying prostate cancer. To this aim, we compiled and analysed currently available post-GWAS data that is validated through further studies in prostate cancer, to investigate molecular biological pathways enriched for assigned functional genes. In total, about 100 canonical pathways were significantly, at false discovery rate (FDR) < 0.05), enriched in assigned genes using different algorithms. The results have highlighted some well-known cancer signalling pathways, antigen presentation processes and enrichment in cell growth and development gene networks, suggesting risk loci may exert their functional effect on prostate cancer by acting through multiple gene sets and pathways. Additional upstream analysis of the involved genes identified critical transcription factors such as HDAC1 and STAT5A. We also investigated the common genes between post-GWAS and three well-annotated gene expression datasets to endeavour to uncover the main genes involved in prostate cancer development/progression. Post-GWAS generated knowledge of gene networks and pathways, although continuously evolving, if analysed further and targeted appropriately, will have an important impact on clinical management of the disease.

DNA methylation and cis-regulation of gene expression by prostate cancer risk SNPs

Genome-wide association studies have identified more than 100 SNPs that increase the risk of prostate cancer (PrCa). We identify and compare expression quantitative trait loci (eQTLs) and CpG methylation quantitative trait loci (meQTLs) among 147 established PrCa risk SNPs in primary prostate tumors (n = 355 from a Seattle-based study and n = 495 from The Cancer Genome Atlas, TCGA) and tumor-adjacent, histologically benign samples (n = 471 from a Mayo Clinic study). The role of DNA methylation in eQTL regulation of gene expression was investigated by data triangulation using several causal inference approaches, including a proposed adaptation of the Causal Inference Test (CIT) for causal direction. Comparing eQTLs between tumors and benign samples, we show that 98 of the 147 risk SNPs were identified as eQTLs in the tumor-adjacent benign samples, and almost all 34 eQTL identified in tumor sets were also eQTLs in the benign samples. Three lines of results support the causal role of DNA methylation. First, nearly 100 of the 147 risk SNPs were identified as meQTLs in one tumor set, and almost all eQTLs in tumors were meQTLs. Second, the loss of eQTLs in tumors relative to benign samples was associated with altered DNA methylation. Third, among risk SNPs identified as both eQTLs and meQTLs, mediation analyses suggest that over two-thirds have evidence of a causal role for DNA methylation, mostly mediating genetic influence on gene expression. In summary, we provide a comprehensive catalog of eQTLs, meQTLs and putative cancer genes for known PrCa risk SNPs. We observe that a substantial portion of germline eQTL regulatory mechanisms are maintained in the tumor development, despite somatic alterations in tumor genome. Finally, our mediation analyses illuminate the likely intermediary role of CpG methylation in eQTL regulation of gene expression.

Integrative comparison of the genomic and transcriptomic landscape between prostate cancer patients of predominantly African or European genetic ancestry

Men of predominantly African Ancestry (AA) have higher prostate cancer (CaP) incidence and worse survival than men of predominantly European Ancestry (EA). While socioeconomic factors drive this disparity, genomic factors may also contribute to differences in the incidence and mortality rates. To compare the prevalence of prostate tumor genomic alterations and transcriptomic profiles by patient genetic ancestry, we evaluated genomic profiles from The Cancer Genome Atlas (TCGA) CaP cohort (n = 498). Patient global and local genetic ancestry were estimated by computational algorithms using genotyping data; 414 (83.1%) were EA, 61 (12.2%) were AA, 11 (2.2%) were East Asian Ancestry (EAA), 10 (2.0%) were Native American (NA), and 2 (0.4%) were other ancestry. Genetic ancestry was highly concordant with self-identified race/ethnicity. Subsequent analyses were limited to 61 AA and 414 EA cases. Significant differences were observed by ancestry in the frequency of SPOP mutations (20.3% AA vs. 10.0% EA; p = 5.6×10-03), TMPRSS2-ERG fusions (29.3% AA vs. 39.6% EA; p = 4.4×10-02), and PTEN deletions/losses (11.5% AA vs. 30.2% EA; p = 3.5×10-03). Differentially expressed genes (DEGs) between AAs and EAs showed significant enrichment for prostate eQTL target genes (p = 8.09×10-48). Enrichment of highly expressed DEGs for immune-related pathways was observed in AAs, and for PTEN/PI3K signaling in EAs. Nearly one-third of DEGs (31.3%) were long non-coding RNAs (DE-lncRNAs). The proportion of DE-lncRNAs with higher expression in AAs greatly exceeded that with lower expression in AAs (p = 1.2×10-125). Both ChIP-seq and RNA-seq data suggested a stronger regulatory role for AR signaling pathways in DE-lncRNAs vs. non-DE-lncRNAs. CaP-related oncogenic lncRNAs, such as PVT1, PCAT1 and PCAT10/CTBP1-AS, were found to be more highly expressed in AAs. We report substantial heterogeneity in the prostate tumor genome and transcriptome between EA and AA. These differences may be biological contributors to racial disparities in CaP incidence and outcomes.

Discovery of rare coding variants in OGDHL and BRCA2 in relation to breast cancer risk in Chinese women

The missing heritability of breast cancer could be partially attributed to rare variants (MAF < 0.5%). To identify breast cancer-associated rare coding variants, we conducted whole-exome sequencing (~50×) in genomic DNA samples obtained from 831 breast cancer cases and 839 controls of Chinese females. Using burden tests for each gene that included rare missense or predicted deleterious variants, we identified 29 genes showing promising associations with breast cancer risk. We replicated the association for two genes, OGDHL and BRCA2, at a Bonferroni-corrected p < 0.05, by genotyping an independent set of samples from 1,628 breast cancer cases and 1,943 controls. The association for OGDHL was primarily driven by three predicted deleterious variants (p.Val827Met, p.Pro839Leu, p.Phe836Ser; p < 0.01 for all). For BRCA2, we characterized a total of 27 disruptive variants, including 18 nonsense, six frameshift and three splicing variants, whereas they were only detected in cases, but none of the controls. All of these variants were either very rare (AF < 0.1%) or not detected in >4,500 East Asian women from the genome Aggregation database (gnomAD), providing additional support to our findings. Our study revealed a potential novel gene and multiple disruptive variants of BRCA2 for breast cancer risk, which may identify high-risk women in Chinese populations.

Single-nucleotide polymorphism rs13426236 contributes to an increased prostate cancer risk via regulating MLPH splicing variant 4

A prostate cancer risk single-nucleotide polymorphism (SNP), rs13426236, is significantly associated with melanophilin (MLPH) expression. To functionally characterize role of the rs13426236 in prostate cancer, we first performed splicing-specific expression quantitative trait loci analysis and refined the significant association of rs13426236 allele G with an increased expression of MLPH splicing transcript variant 4 (V4) (P = 7.61E-5) but not other protein-coding variants (V1-V3) (P > .05). We then performed an allele-specific reporter assay to determine if SNP-containing sequences functioned as an active enhancer. Compared to allele A, allele G of rs13426236 showed significantly higher luciferase activity on the promoter of the splicing transcript V4 (P < .03) but not on the promoter of transcript V1 (P > .05) in two prostate cancer cell lines (DU145 and 22Rv1). Cell transfection assays showed stronger effect of transcript V4 than V1 on promoting cell proliferation, invasion, and antiapoptotic activities. RNA profiling analysis demonstrated that transcript V4 overexpression caused significant expression changes in glycosylation/glycoprotein and metal-binding gene ontology pathways (FDR < 0.01). We also found that both transcripts V4 and V1 were significantly upregulated in prostate adenocarcinoma (P ≤ 2.49E-6) but only transcript V4 upregulation was associated with poor recurrence-free survival (P = .028, hazard ratio = 1.63, 95% confidence interval = 1.05-2.42) in The Cancer Genome Atlas data. This study provides strong evidence showing that prostate cancer risk SNP rs13426236 upregulates expression of MLPH transcript V4, which may function as a candidate oncogene in prostate cancer.

Integrative genomic analyses of APOBEC-mutational signature, expression and germline deletion of APOBEC3 genes, and immunogenicity in multiple cancer types

BACKGROUND

Although APOBEC-mutational signature is found in tumor tissues of multiple cancers, how a common germline APOBEC3A/B deletion affects the mutational signature remains unclear.

METHODS

Using data from 10 cancer types generated as part of TCGA, we performed integrative genomic and association analyses to assess inter-relationship of expressions for isoforms APOBEC3A and APOBEC3B, APOBEC-mutational signature, germline APOBEC3A/B deletions, neoantigen loads, and tumor infiltration lymphocytes (TILs).

RESULTS

We found that expression level of the isoform uc011aoc transcribed from the APOBEC3A/B chimera was associated with a greater burden of APOBEC-mutational signature only in breast cancer, while germline APOBEC3A/B deletion led to an increased expression level of uc011aoc in multiple cancer types. Furthermore, we found that the deletion was associated with elevated APOBEC-mutational signature, neoantigen loads and relative composition of T cells (CD8+) in TILs only in breast cancer. Additionally, we also found that APOBEC-mutational signature significantly contributed to neoantigen loads and certain immune cell abundances in TILs across cancer types.

CONCLUSIONS

These findings reveal new insights into understanding the genetic, biological and immunological mechanisms through which APOBEC genes may be involved in carcinogenesis, and provide potential genetic biomarker for the development of disease prevention and cancer immunotherapy.

Identifying Putative Susceptibility Genes and Evaluating Their Associations with Somatic Mutations in Human Cancers

Genome-wide association studies (GWASs) have identified hundreds of genetic risk variants for human cancers. However, target genes for the majority of risk loci remain largely unexplored. It is also unclear whether GWAS risk-loci-associated genes contribute to mutational signatures and tumor mutational burden (TMB) in cancer tissues. We systematically conducted cis-expression quantitative trait loci (cis-eQTL) analyses for 294 GWAS-identified variants for six major types of cancer-colorectal, lung, ovary, prostate, pancreas, and melanoma-by using transcriptome data from the Genotype-Tissue Expression (GTEx) Project, the Cancer Genome Atlas (TCGA), and other public data sources. By using integrative analysis strategies, we identified 270 candidate target genes, including 99 with previously unreported associations, for six cancer types. By analyzing functional genomic data, our results indicate that 180 genes (66.7% of 270) had evidence of cis-regulation by putative functional variants via proximal promoter or distal enhancer-promoter interactions. Together with our previously reported associations for breast cancer risk, our results show that 24 genes are shared by at least two cancer types, including four genes for both breast and ovarian cancer. By integrating mutation data from TCGA, we found that expression levels of 33 and 66 putative susceptibility genes were associated with specific mutational signatures and TMB of cancer-driver genes, respectively, at a Bonferroni-corrected p < 0.05. Together, these findings provide further insight into our understanding of how genetic risk variants might contribute to carcinogenesis through the regulation of susceptibility genes that are related to the biogenesis of somatic mutations.

Analysis of Over 140,000 European Descendants Identifies Genetically Predicted Blood Protein Biomarkers Associated with Prostate Cancer Risk

Several blood protein biomarkers have been associated with prostate cancer risk. However, most studies assessed only a small number of biomarkers and/or included a small sample size. To identify novel protein biomarkers of prostate cancer risk, we studied 79,194 cases and 61,112 controls of European ancestry, included in the PRACTICAL/ELLIPSE consortia, using genetic instruments of protein quantitative trait loci for 1,478 plasma proteins. A total of 31 proteins were associated with prostate cancer risk including proteins encoded by GSTP1, whose methylation level was shown previously to be associated with prostate cancer risk, and MSMB, SPINT2, IGF2R, and CTSS, which were previously implicated as potential target genes of prostate cancer risk variants identified in genome-wide association studies. A total of 18 proteins inversely correlated and 13 positively correlated with prostate cancer risk. For 28 of the identified proteins, gene somatic changes of short indels, splice site, nonsense, or missense mutations were detected in patients with prostate cancer in The Cancer Genome Atlas. Pathway enrichment analysis showed that relevant genes were significantly enriched in cancer-related pathways. In conclusion, this study identifies 31 candidates of protein biomarkers for prostate cancer risk and provides new insights into the biology and genetics of prostate tumorigenesis. SIGNIFICANCE: Integration of genomics and proteomics data identifies biomarkers associated with prostate cancer risk.

Single-Nucleotide Polymorphisms Sequencing Identifies Candidate Functional Variants at Prostate Cancer Risk Loci

Genome-wide association studies have identified over 150 risk loci that increase prostate cancer risk. However, few causal variants and their regulatory mechanisms have been characterized. In this study, we utilized our previously developed single-nucleotide polymorphisms sequencing (SNPs-seq) technology to test allele-dependent protein binding at 903 SNP sites covering 28 genomic regions. All selected SNPs have shown significant cis-association with at least one nearby gene. After preparing nuclear extract using LNCaP cell line, we first mixed the extract with dsDNA oligo pool for protein–DNA binding incubation. We then performed sequencing analysis on protein-bound oligos. SNPs-seq analysis showed protein-binding differences (>1.5-fold) between reference and variant alleles in 380 (42%) of 903 SNPs with androgen treatment and 403 (45%) of 903 SNPs without treatment. From these significant SNPs, we performed a database search and further narrowed down to 74 promising SNPs. To validate this initial finding, we performed electrophoretic mobility shift assay in two SNPs (rs12246440 and rs7077275) at CTBP2 locus and one SNP (rs113082846) at NCOA4 locus. This analysis showed that all three SNPs demonstrated allele-dependent protein-binding differences that were consistent with the SNPs-seq. Finally, clinical association analysis of the two candidate genes showed that CTBP2 was upregulated, while NCOA4 was downregulated in prostate cancer (p < 0.02). Lower expression of CTBP2 was associated with poor recurrence-free survival in prostate cancer. Utilizing our experimental data along with bioinformatic tools provides a strategy for identifying candidate functional elements at prostate cancer susceptibility loci to help guide subsequent laboratory studies.

Association of imputed prostate cancer transcriptome with disease risk reveals novel mechanisms

Here we train cis-regulatory models of prostate tissue gene expression and impute expression transcriptome-wide for 233,955 European ancestry men (14,616 prostate cancer (PrCa) cases, 219,339 controls) from two large cohorts. Among 12,014 genes evaluated in the UK Biobank, we identify 38 associated with PrCa, many replicating in the Kaiser Permanente RPGEH. We report the association of elevated TMPRSS2 expression with increased PrCa risk (independent of a previously-reported risk variant) and with increased tumoral expression of the TMPRSS2:ERG fusion-oncogene in The Cancer Genome Atlas, suggesting a novel germline-somatic interaction mechanism. Three novel genes, HOXA4, KLK1, and TIMM23, additionally replicate in the RPGEH cohort. Furthermore, 4 genes, MSMB, NCOA4, PCAT1, and PPP1R14A, are associated with PrCa in a trans-ethnic meta-analysis (N = 9117). Many genes exhibit evidence for allele-specific transcriptional activation by PrCa master-regulators (including androgen receptor) in Position Weight Matrix, Chip-Seq, and Hi-C experimental data, suggesting common regulatory mechanisms for the associated genes.

Identification of Novel Susceptibility Loci and Genes for Prostate Cancer Risk: A Transcriptome-Wide Association Study in Over 140,000 European Descendants

Genome-wide association study-identified prostate cancer risk variants explain only a relatively small fraction of its familial relative risk, and the genes responsible for many of these identified associations remain unknown. To discover novel prostate cancer genetic loci and possible causal genes at previously identified risk loci, we performed a transcriptome-wide association study in 79,194 cases and 61,112 controls of European ancestry. Using data from the Genotype-Tissue Expression Project, we established genetic models to predict gene expression across the transcriptome for both prostate models and cross-tissue models and evaluated model performance using two independent datasets. We identified significant associations for 137 genes at P < 2.61 × 10-6, a Bonferroni-corrected threshold, including nine genes that remained significant at P < 2.61 × 10-6 after adjusting for all known prostate cancer risk variants in nearby regions. Of the 128 remaining associated genes, 94 have not yet been reported as potential target genes at known loci. We silenced 14 genes and many showed a consistent effect on viability and colony-forming efficiency in three cell lines. Our study provides substantial new information to advance our understanding of prostate cancer genetics and biology. SIGNIFICANCE: This study identifies novel prostate cancer genetic loci and possible causal genes, advancing our understanding of the molecular mechanisms that drive prostate cancer.

The prostate cancer risk variant rs55958994 regulates multiple gene expression through extreme long-range chromatin interaction to control tumor progression

Genome-wide association studies identified single-nucleotide polymorphism (SNP) rs55958994 as a significant variant associated with increased susceptibility to prostate cancer. However, the mechanisms by which this SNP mediates increased risk to cancer are still unknown. In this study, we show that this variant is located in an enhancer active in prostate cancer cells. Deletion of this enhancer from prostate tumor cells resulted in decreased tumor initiation, tumor growth, and invasive migration, as well as a loss of stem-like cells. Using a combination of capture chromosome conformation capture (Capture-C) and RNA sequencing, we identified genes on the same and different chromosomes as targets regulated by the enhancer. Furthermore, we show that expression of individual candidate target genes in an enhancer-deleted cell line rescued different aspects of tumorigenesis. Our data suggest that the rs55958994-associated enhancer affects prostate cancer progression by influencing expression of multiple genes via long-range chromatin interactions.

Identification of novel genetic variants for type 2 diabetes, childhood obesity, and their pleiotropic loci

Obesity has result in increased prevalence of type 2 diabetes (T2D) in children. The genetic mechanisms underlying their relationship, however, are not fully understood. Here, we aim to identify novel SNPs associated with T2D and childhood obesity (CO), especially their pleiotropic loci. We integrated the summary statistics for two independent GWASs of T2D (n = 149,821) and childhood body mass index (CBMI) (n = 35,668) using the pleiotropy-informed conditional false discovery rate (cFDR) method. By leveraging the information of different levels of association for CBMI, we observed a strong enrichment of genetic variants associated with T2D. We identified 139 T2D-associated SNPs with 125 novel ones (cFDR < 0.05). Conditioned on T2D, we identified 37 significant SNPs for CBMI (cFDR < 0.05), including 25 novel ones. The conjunctional cFDR (ccFDR) analysis showed ten novel pleiotropic loci for T2D and CBMI (ccFDR < 0.05). Interestingly, the novel SNP rs1996023 is located at protein coding gene GNPDA2 (ccFDR = 1.28E-02), which has been reported to influence the risk of T2D and CO through central nervous system. Our findings may help to explain a greater proportion of the heritability for human traits and advance the understanding of the common pathophysiology between T2D and CO.

An expanded variant list and assembly annotation identifies multiple novel coding and noncoding genes for prostate cancer risk using a normal prostate tissue eQTL data set

Prostate cancer (PrCa) is highly heritable; 284 variants have been identified to date that are associated with increased prostate cancer risk, yet few genes contributing to its development are known. Expression quantitative trait loci (eQTL) studies link variants with affected genes, helping to determine how these variants might regulate gene expression and may influence prostate cancer risk. In the current study, we performed eQTL analysis on 471 normal prostate epithelium samples and 249 PrCa-risk variants in 196 risk loci, utilizing RNA sequencing transcriptome data based on ENSEMBL gene definition and genome-wide variant data. We identified a total of 213 genes associated with known PrCa-risk variants, including 141 protein-coding genes, 16 lncRNAs, and 56 other non-coding RNA species with differential expression. Compared to our previous analysis, where RefSeq was used for gene annotation, we identified an additional 130 expressed genes associated with known PrCa-risk variants. We detected an eQTL signal for more than half (n = 102, 52%) of the 196 loci tested; 52 (51%) of which were a Group 1 signal, indicating high linkage disequilibrium (LD) between the peak eQTL variant and the PrCa-risk variant (r2>0.5) and may help explain how risk variants influence the development of prostate cancer.

An AR-ERG transcriptional signature defined by long-range chromatin interactomes in prostate cancer cells

The aberrant activities of transcription factors such as the androgen receptor (AR) underpin prostate cancer development. While the AR cis-regulation has been extensively studied in prostate cancer, information pertaining to the spatial architecture of the AR transcriptional circuitry remains limited. In this paper, we propose a novel framework to profile long-range chromatin interactions associated with AR and its collaborative transcription factor, erythroblast transformation-specific related gene (ERG), using chromatin interaction analysis by paired-end tag (ChIA-PET). We identified ERG-associated long-range chromatin interactions as a cooperative component in the AR-associated chromatin interactome, acting in concert to achieve coordinated regulation of a subset of AR target genes. Through multifaceted functional data analysis, we found that AR-ERG interaction hub regions are characterized by distinct functional signatures, including bidirectional transcription and cotranscription factor binding. In addition, cancer-associated long noncoding RNAs were found to be connected near protein-coding genes through AR-ERG looping. Finally, we found strong enrichment of prostate cancer genome-wide association study (GWAS) single nucleotide polymorphisms (SNPs) at AR-ERG co-binding sites participating in chromatin interactions and gene regulation, suggesting GWAS target genes identified from chromatin looping data provide more biologically relevant findings than using the nearest gene approach. Taken together, our results revealed an AR-ERG-centric higher-order chromatin structure that drives coordinated gene expression in prostate cancer progression and the identification of potential target genes for therapeutic intervention.

Post-GWAS in prostate cancer: from genetic association to biological contribution

Genome-wide association studies (GWAS) have been successful in deciphering the genetic component of predisposition to many human complex diseases including prostate cancer. Germline variants identified by GWAS progressively unravelled the substantial knowledge gap concerning prostate cancer heritability. With the beginning of the post-GWAS era, more and more studies reveal that, in addition to their value as risk markers, germline variants can exert active roles in prostate oncogenesis. Consequently, current research efforts focus on exploring the biological mechanisms underlying specific susceptibility loci known as causal variants by applying novel and precise analytical methods to available GWAS data. Results obtained from these post-GWAS analyses have highlighted the potential of exploiting prostate cancer risk-associated germline variants to identify new gene networks and signalling pathways involved in prostate tumorigenesis. In this Review, we describe the molecular basis of several important prostate cancer-causal variants with an emphasis on using post-GWAS analysis to gain insight into cancer aetiology. In addition to discussing the current status of post-GWAS studies, we also summarize the main molecular mechanisms of potential causal variants at prostate cancer risk loci and explore the major challenges in moving from association to functional studies and their implication in clinical translation.

Genetic variation in transcription factors and photosynthesis light-reaction genes regulates photosynthetic traits

Transcription factors (TFs) play crucial roles in regulating the production of the components required for photosynthesis; elucidating the mechanisms by which underlying genetic variation in TFs affects complex photosynthesis-related traits may improve our understanding of photosynthesis and identify ways to improve photosynthetic efficiency. Promoter analysis of 96 nuclear-encoded Populus tomentosa Carr. genes within this pathway revealed 47 motifs responsive to light, stress, hormones and organ-specific regulation, as well as 86 TFs that might bind these motifs. Using phenotype-genotype associations, we identified 244 single-nucleotide polymorphisms (SNPs) within 105 genes associated with 12 photosynthesis-related traits. Most (30.33%) of these SNPs were located in intronic regions and these SNPs explained 18.66% of the mean phenotypic variation in the photosynthesis-related traits. Additionally, expression quantitative trait loci (eQTL) mapping identified 216 eQTLs associated with 110 eGenes (genes regulated by eQTLs), explaining 14.12% of the variability of gene expression. The lead SNPs of 12.04% of the eQTLs also contributed to phenotypic variation. Among these, a SNP in zf-Dof 5.6 (G120_9287) affected photosynthesis by modulating the expression of a sub-regulatory network of eight other TFs, which in turn regulate 55 photosynthesis-related genes. Furthermore, epistasis analysis identified a large interacting network representing 732 SNP-SNP pairs, of which 354 were photosynthesis gene-TF pairs, emphasizing the important roles of TFs in affecting photosynthesis-related traits. We combined eQTL and epistasis analysis and found 32 TFs harboring eQTLs being epistatic to their targets (identified by eQTL analysis), of which 15 TFs were also associated with photosynthesis traits. We therefore constructed a schematic model of TFs involved in regulating the photosynthetic light reaction pathway. Taken together, our results provide insight into the genetic regulation of photosynthesis, and may drive progress in the marker-assisted selection of desirable P. tomentosa genotypes with more efficient photosynthesis.

Deciphering the Emerging Complexities of Molecular Mechanisms at GWAS Loci

Genome-wide association studies (GWASs) have identified thousands of loci associated with hundreds of complex diseases and traits, and progress is being made toward elucidating the causal variants and genes underlying these associations. Functional characterization of mechanisms at GWAS loci is a multi-faceted challenge. Challenges include linkage disequilibrium and allelic heterogeneity at each locus, the noncoding nature of most loci, and the time and cost needed for experimentally evaluating the potential mechanistic contributions of genes and variants. As GWAS sample sizes increase, more loci are identified, and the complexities of individual loci emerge. Loci can consist of multiple association signals, each of which can reflect the influence of multiple variants, inseparable by association analyses. Each signal within a locus can influence the same or different target genes. Experimental studies of genes and variants can differ on the basis of cell type, cellular environment, or other context-specific variables. In this review, we describe the complexity of mechanisms at GWAS loci-including multiple signals, multiple variants, and/or multiple genes-and the implications these complexities hold for experimental study design and interpretation of GWAS mechanisms.

Large-scale transcriptome-wide association study identifies new prostate cancer risk regions

Although genome-wide association studies (GWAS) for prostate cancer (PrCa) have identified more than 100 risk regions, most of the risk genes at these regions remain largely unknown. Here we integrate the largest PrCa GWAS (N = 142,392) with gene expression measured in 45 tissues (N = 4458), including normal and tumor prostate, to perform a multi-tissue transcriptome-wide association study (TWAS) for PrCa. We identify 217 genes at 84 independent 1 Mb regions associated with PrCa risk, 9 of which are regions with no genome-wide significant SNP within 2 Mb. 23 genes are significant in TWAS only for alternative splicing models in prostate tumor thus supporting the hypothesis of splicing driving risk for continued oncogenesis. Finally, we use a Bayesian probabilistic approach to estimate credible sets of genes containing the causal gene at a pre-defined level; this reduced the list of 217 associations to 109 genes in the 90% credible set. Overall, our findings highlight the power of integrating expression with PrCa GWAS to identify novel risk loci and prioritize putative causal genes at known risk loci.

Identification of expression quantitative trait loci associated with schizophrenia and affective disorders in normal brain tissue

Schizophrenia and the affective disorders, here comprising bipolar disorder and major depressive disorder, are psychiatric illnesses that lead to significant morbidity and mortality worldwide. Whilst understanding of their pathobiology remains limited, large case-control studies have recently identified single nucleotide polymorphisms (SNPs) associated with these disorders. However, discerning the functional effects of these SNPs has been difficult as the associated causal genes are unknown. Here we evaluated whether schizophrenia and affective disorder associated-SNPs are correlated with gene expression within human brain tissue. Specifically, to identify expression quantitative trait loci (eQTLs), we leveraged disorder-associated SNPs identified from 11 genome-wide association studies with gene expression levels in post-mortem, neurologically-normal tissue from two independent human brain tissue expression datasets (UK Brain Expression Consortium (UKBEC) and Genotype-Tissue Expression (GTEx)). Utilizing stringent multi-region meta-analyses, we identified 2,224 cis-eQTLs associated with expression of 40 genes, including 11 non-coding RNAs. One cis-eQTL, rs16969968, results in a functionally disruptive missense mutation in CHRNA5, a schizophrenia-implicated gene. Importantly, comparing across tissues, we find that blood eQTLs capture < 10% of brain cis-eQTLs. Contrastingly, > 30% of brain-associated eQTLs are significant in tibial nerve. This study identifies putatively causal genes whose expression in region-specific tissue may contribute to the risk of schizophrenia and affective disorders.

eQTL discovery and their association with severe equine asthma in European Warmblood horses

BACKGROUND

Severe equine asthma, also known as recurrent airway obstruction (RAO), is a debilitating, performance limiting, obstructive respiratory condition in horses that is phenotypically similar to human asthma. Past genome wide association studies (GWAS) have not discovered coding variants associated with RAO, leading to the hypothesis that causative variant(s) underlying the signals are likely non-coding, regulatory variant(s). Regions of the genome containing variants that influence the number of expressed RNA molecules are expression quantitative trait loci (eQTLs). Variation associated with RAO that also regulates a gene's expression in a disease relevant tissue could help identify candidate genes that influence RAO if that gene's expression is also associated with RAO disease status.

RESULTS

We searched for eQTLs by analyzing peripheral blood mononuclear cells (PBMCs) from two half-sib families and one unrelated cohort of 82 European Warmblood horses that were previously treated in vitro with: no stimulation (MCK), lipopolysaccharides (LPS), recombinant cyathostomin antigen (RCA), and hay-dust extract (HDE). We identified high confidence eQTLs that did not violate linear modeling assumptions and were not significant due to single outlier individuals. We identified a mean of 4347 high confidence eQTLs in four treatments of PBMCs, and discovered two trans regulatory hotspots regulating genes involved in related biological pathways. We corroborated previous RAO associated single nucleotide polymorphisms (SNPs), and increased the resolution of past GWAS by analyzing 1,056,195 SNPs in 361 individuals. We identified four RAO-associated SNPs that only regulate gene expression of dexamethasone-induced protein (DEXI), however we found no significant association between DEXI gene expression and presence of RAO.

CONCLUSIONS

Thousands of genetic variants regulate gene expression in PBMCs of European Warmblood horses in cis and trans. Most high confidence eSNPs are significantly enriched near the transcription start sites of their target genes. Two trans regulatory hotspots on chromosome 11 and 13 regulate many genes involved in transmembrane cell signaling and neurological development respectively when PBMCs are treated with HDE. None of the top fifteen RAO associated SNPs strongly influence disease status through gene expression regulation.

High-throughput screening of prostate cancer risk loci by single nucleotide polymorphisms sequencing

Functional characterization of disease-causing variants at risk loci has been a significant challenge. Here we report a high-throughput single-nucleotide polymorphisms sequencing (SNPs-seq) technology to simultaneously screen hundreds to thousands of SNPs for their allele-dependent protein-binding differences. This technology takes advantage of higher retention rate of protein-bound DNA oligos in protein purification column to quantitatively sequence these SNP-containing oligos. We apply this technology to test prostate cancer-risk loci and observe differential allelic protein binding in a significant number of selected SNPs. We also test a unique application of self-transcribing active regulatory region sequencing (STARR-seq) in characterizing allele-dependent transcriptional regulation and provide detailed functional analysis at two risk loci (RGS17 and ASCL2). Together, we introduce a powerful high-throughput pipeline for large-scale screening of functional SNPs at disease risk loci.

Functional characterization of disease-causing variants at risk loci in cancer is challenging. Here, in prostate cancer the authors report a pipeline for high-throughput single-nucleotide polymorphisms sequencing (SNPs-seq) for large scale screening of functional SNPs at disease risk loci.

Single Nucleotide Polymorphism Genotyping in Single-Molecule Electronic Circuits

Establishing low-cost, high-throughput, simple, and accurate single nucleotide polymorphism (SNP) genotyping techniques is beneficial for understanding the intrinsic relationship between individual genetic variations and their biological functions on a genomic scale. Here, a straightforward and reliable single-molecule approach is demonstrated for precise SNP authentication by directly measuring the fluctuations in electrical signals in an electronic circuit, which is fabricated from a high-gain field-effect silicon nanowire decorated with a single hairpin DNA, in the presence of different target DNAs. By simply comparing the proportion difference of a probe-target duplex structure throughout the process, this study implements allele-specific and accurate SNP detection. These results are supported by the statistical analyses of different dynamic parameters such as the mean lifetime and the unwinding probability of the duplex conformation. In comparison with conventional polymerase chain reaction and optical methods, this convenient and label-free method is complementary to existing optical methods and also shows several advantages, such as simple operation and no requirement for fluorescent labeling, thus promising a futuristic route toward the next-generation genotyping technique.