Gene regulatory mechanisms underpinning prostate cancer susceptibility

Cancer-associated fibroblasts (CAFs) gene signatures predict outcomes in breast and prostate tumor patients

BACKGROUND

Over the last two decades, tumor-derived RNA expression signatures have been developed for the two most commonly diagnosed tumors worldwide, namely prostate and breast tumors, in order to improve both outcome prediction and treatment decision-making. In this context, molecular signatures gained by main components of the tumor microenvironment, such as cancer-associated fibroblasts (CAFs), have been explored as prognostic and therapeutic tools. Nevertheless, a deeper understanding of the significance of CAFs-related gene signatures in breast and prostate cancers still remains to be disclosed.

METHODS

RNA sequencing technology (RNA-seq) was employed to profile and compare the transcriptome of CAFs isolated from patients affected by breast and prostate tumors. The differentially expressed genes (DEGs) characterizing breast and prostate CAFs were intersected with data from public datasets derived from bulk RNA-seq profiles of breast and prostate tumor patients. Pathway enrichment analyses allowed us to appreciate the biological significance of the DEGs. K-means clustering was applied to construct CAFs-related gene signatures specific for breast and prostate cancer and to stratify independent cohorts of patients into high and low gene expression clusters. Kaplan-Meier survival curves and log-rank tests were employed to predict differences in the outcome parameters of the clusters of patients. Decision-tree analysis was used to validate the clustering results and boosting calculations were then employed to improve the results obtained by the decision-tree algorithm.

RESULTS

Data obtained in breast CAFs allowed us to assess a signature that includes 8 genes (ITGA11, THBS1, FN1, EMP1, ITGA2, FYN, SPP1, and EMP2) belonging to pro-metastatic signaling routes, such as the focal adhesion pathway. Survival analyses indicated that the cluster of breast cancer patients showing a high expression of the aforementioned genes displays worse clinical outcomes. Next, we identified a prostate CAFs-related signature that includes 11 genes (IL13RA2, GDF7, IL33, CXCL1, TNFRSF19, CXCL6, LIFR, CXCL5, IL7, TSLP, and TNFSF15) associated with immune responses. A low expression of these genes was predictive of poor survival rates in prostate cancer patients. The results obtained were significantly validated through a two-step approach, based on unsupervised (clustering) and supervised (classification) learning techniques, showing a high prediction accuracy (≥ 90%) in independent RNA-seq cohorts.

CONCLUSION

We identified a huge heterogeneity in the transcriptional profile of CAFs derived from breast and prostate tumors. Of note, the two novel CAFs-related gene signatures might be considered as reliable prognostic indicators and valuable biomarkers for a better management of breast and prostate cancer patients.

An Inherited Allele Confers Prostate Cancer Progression and Drug Resistance via RFX6/HOXA10-Orchestrated TGFβ Signaling

Genetic and epigenetic alterations are cancer hallmark characteristics. However, the role of inherited cancer predisposition alleles in co-opting lineage factor epigenetic reprogramming and tumor progression remains elusive. Here the FinnGen cohort phenome-wide analysis, along with multiple genome-wide association studies, has consistently identified the rs339331-RFX6/6q22 locus associated with prostate cancer (PCa) risk across diverse populations. It is uncovered that rs339331 resides in a reprogrammed androgen receptor (AR) binding site in PCa tumors, with the T risk allele enhancing AR chromatin occupancy. RFX6, an AR-regulated gene linked to rs339331, exhibits synergistic prognostic value for PCa recurrence and metastasis. This comprehensive in vitro and in vivo studies demonstrate the oncogenic functions of RFX6 in promoting PCa cell proliferation and metastasis. Mechanistically, RFX6 upregulates HOXA10 that profoundly correlates with adverse PCa outcomes and is pivotal in RFX6-mediated PCa progression, facilitating the epithelial-mesenchymal transition (EMT) and modulating the TGFβ/SMAD signaling axis. Clinically, HOXA10 elevation is associated with increased EMT scores, tumor advancement and PCa recurrence. Remarkably, reducing RFX6 expression restores enzalutamide sensitivity in resistant PCa cells and tumors. This findings reveal a complex interplay of genetic and epigenetic mechanisms in PCa pathogenesis and drug resistance, centered around disrupted prostate lineage AR signaling and abnormal RFX6 expression.

CRISPR genome-wide screening identifies PAK1 as a critical driver of ARSI cross-resistance in prostate cancer progression

Next-generation androgen receptor signaling inhibitors (ARSIs), such as enzalutamide (Enza) and darolutamide (Daro), are initially effective for the treatment of advanced prostate cancer (PCa) and castration-resistant prostate cancer (CRPC). However, patients often relapse and develop cross-resistance, which consequently makes drug resistance an inevitable cause of CRPC-related mortality. By conducting a comprehensive analysis of GEO datasets, CRISPR genome-wide screening results, ATAC-seq data, and RNA-seq data, we systemically identified PAK1 as a significant contributor to ARSI cross-resistance due to the activation of the PAK1/RELA/hnRNPA1/AR-V7 axis. Inhibition of PAK1 followed by suppression of NF-κB pathways and AR-V7 expression effectively overcomes ARSI cross-resistance. Our findings indicate that PAK1 represents a promising therapeutic target gene for the treatment of ARSI cross-resistant PCa patients in the clinic. STATEMENT OF SIGNIFICANCE: PAK1 drives ARSI cross-resistance in prostate cancer progression.

Dampened Regulatory Circuitry of TEAD1/ITGA1/ITGA2 Promotes TGFβ1 Signaling to Orchestrate Prostate Cancer Progression

The extracellular matrix (ECM) undergoes substantial changes during prostate cancer (PCa) progression, thereby regulating PCa growth and invasion. Herein, a meta-analysis of multiple PCa cohorts is performed which revealed that downregulation or genomic loss of ITGA1 and ITGA2 integrin genes is associated with tumor progression and worse prognosis. Genomic deletion of both ITGA1 and ITGA2 activated epithelial-to-mesenchymal transition (EMT) in benign prostate epithelial cells, thereby enhancing their invasive potential in vitro and converting them into tumorigenic cells in vivo. Mechanistically, EMT is induced by enhanced secretion and autocrine activation of TGFβ1 and nuclear targeting of YAP1. An unbiased genome-wide co-expression analysis of large PCa cohort datasets identified the transcription factor TEAD1 as a key regulator of ITGA1 and ITGA2 expression in PCa cells while TEAD1 loss phenocopied the dual loss of α1- and α2-integrins in vitro and in vivo. Remarkably, clinical data analysis revealed that TEAD1 downregulation or genomic loss is associated with aggressive PCa and together with low ITGA1 and ITGA2 expression synergistically impacted PCa prognosis and progression. This study thus demonstrated that loss of α1- and α2-integrins, either via deletion/inactivation of the ITGA1/ITGA2 locus or via loss of TEAD1, contributes to PCa progression by inducing TGFβ1-driven EMT.

Genome-wide meta-analysis, functional genomics and integrative analyses implicate new risk genes and therapeutic targets for anxiety disorders

Anxiety disorders are the most prevalent mental disorders. However, the genetic etiology of anxiety disorders remains largely unknown. Here we conducted a genome-wide meta-analysis on anxiety disorders by including 74,973 (28,392 proxy) cases and 400,243 (146,771 proxy) controls. We identified 14 risk loci, including 10 new associations near CNTNAP5, MAP2, RAB9BP1, BTN1A1, PRR16, PCLO, PTPRD, FARP1, CDH2 and RAB27B. Functional genomics and fine-mapping pinpointed the potential causal variants, and expression quantitative trait loci analysis revealed the potential target genes regulated by the risk variants. Integrative analyses, including transcriptome-wide association study, proteome-wide association study and colocalization analyses, prioritized potential causal genes (including CTNND1 and RAB27B). Evidence from multiple analyses revealed possibly causal genes, including RAB27B, BTN3A2, PCLO and CTNND1. Finally, we showed that Ctnnd1 knockdown affected dendritic spine density and resulted in anxiety-like behaviours in mice, revealing the potential role of CTNND1 in anxiety disorders. Our study identified new risk loci, potential causal variants and genes for anxiety disorders, providing insights into the genetic architecture of anxiety disorders and potential therapeutic targets.

The diagnostic or prognostic values of FADD in cancers based on pan‑cancer analysis

Previous studies have determined that aberrant expression of the fas-associated death domain (FADD) contributes to the development of cancer. However, no pan-cancer analysis has been reported to explore the relationship between FADD and various cancers. Multiple databases were screened to identify cancer datasets for the present study and to validate the expression of FADD in various tumors. The association of FADD alteration with cancer prognosis, clinical features and tumor immunity was also evaluated. Reverse transcription-quantitative PCR (RT-qPCR) was utilized to confirm the expression of FADD in breast, colon, liver and gastric cancer cells. Analysis of Gene Expression Omnibus database and The Cancer Genome Atlas database indicated that FADD was highly expressed in breast invasive carcinoma (BRCA), cervical squamous cell carcinoma and endocervical adenocarcinoma, cholangiocarcinoma, colon adenocarcinoma (COAD), esophageal carcinoma (ESCA), kidney renal clear cell carcinoma, kidney renal papillary cell carcinoma, liver hepatocellular carcinoma (LIHC), lung adenocarcinoma (LUAD) and prostate adenocarcinoma, whereas RT-qPCR results revealed that FADD was highly expressed in breast cancer and colon cancer. Further analyses demonstrated that FADD expression was significantly altered in ESCA, head and neck squamous cell carcinoma (HNSC), lung squamous cell carcinoma and BRCA. FADD expression was observed to be a risk factor of the overall survival in patients with HNSC, LIHC and LUAD as demonstrated by Kaplan-Meier and Cox regression analyses. The results of the present study demonstrated that FADD is highly expressed in numerous malignancies and can be utilized as a biomarker for the diagnosis of BRCA, COAD, LIHC and stomach adenocarcinoma. Moreover, FADD expression is a predictive risk factor for the development of HNSC, LIHC and LUAD and can potentially be used as a prognostic marker for these cancers.

Human-genome single nucleotide polymorphisms affecting transcription factor binding and their role in pathogenesis

Single nucleotide polymorphisms (SNPs) are the most common type of variation in the human genome. The vast majority of SNPs identified in the human genome do not have any effect on the phenotype; however, some can lead to changes in the function of a gene or the level of its expression. Most SNPs associated with certain traits or pathologies are mapped to regulatory regions of the genome and affect gene expression by changing transcription factor binding sites. In recent decades, substantial effort has been invested in searching for such regulatory SNPs (rSNPs) and understanding the mechanisms by which they lead to phenotypic differences, primarily to individual differences in susceptibility to diseases and in sensitivity to drugs. The development of the NGS (next-generation sequencing) technology has contributed not only to the identification of a huge number of SNPs and to the search for their association (genome-wide association studies, GWASs) with certain diseases or phenotypic manifestations, but also to the development of more productive approaches to their functional annotation. It should be noted that the presence of an association does not allow one to identify a functional, truly disease-associated DNA sequence variant among multiple marker SNPs that are detected due to linkage disequilibrium. Moreover, determination of associations of genetic variants with a disease does not provide information about the functionality of these variants, which is necessary to elucidate the molecular mechanisms of the development of pathology and to design effective methods for its treatment and prevention. In this regard, the functional analysis of SNPs annotated in the GWAS catalog, both at the genome-wide level and at the level of individual SNPs, became especially relevant in recent years. A genome-wide search for potential rSNPs is possible without any prior knowledge of their association with a trait. Thus, mapping expression quantitative trait loci (eQTLs) makes it possible to identify an SNP for which - among transcriptomes of homozygotes and heterozygotes for its various alleles - there are differences in the expression level of certain genes, which can be located at various distances from the SNP. To predict rSNPs, approaches based on searches for allele-specific events in RNA-seq, ChIP-seq, DNase-seq, ATAC-seq, MPRA, and other data are also used. Nonetheless, for a more complete functional annotation of such rSNPs, it is necessary to establish their association with a trait, in particular, with a predisposition to a certain pathology or sensitivity to drugs. Thus, approaches to finding SNPs important for the development of a trait can be categorized into two groups: (1) starting from data on an association of SNPs with a certain trait, (2) starting from the determination of allele-specific changes at the molecular level (in a transcriptome or regulome). Only comprehensive use of strategically different approaches can considerably enrich our knowledge about the role of genetic determinants in the molecular mechanisms of trait formation, including predisposition to multifactorial diseases.

GATA2 co-opts TGFβ1/SMAD4 oncogenic signaling and inherited variants at 6q22 to modulate prostate cancer progression

BACKGROUND

Aberrant somatic genomic alteration including copy number amplification is a hallmark of cancer genomes. We previously profiled genomic landscapes of prostate cancer (PCa), yet the underlying causal genes with prognostic potential has not been defined. It remains unclear how a somatic genomic event cooperates with inherited germline variants contribute to cancer predisposition and progression.

METHODS

We applied integrated genomic and clinical data, experimental models and bioinformatic analysis to identify GATA2 as a highly prevalent metastasis-associated genomic amplification in PCa. Biological roles of GATA2 in PCa metastasis was determined in vitro and in vivo. Global chromatin co-occupancy and co-regulation of GATA2 and SMAD4 was investigated by coimmunoprecipitation, ChIP-seq and RNA-seq assays. Tumor cellular assays, qRT-PCR, western blot, ChIP, luciferase assays and CRISPR-Cas9 editing methods were performed to mechanistically understand the cooperation of GATA2 with SMAD4 in promoting TGFβ1 and AR signaling and mediating inherited PCa risk and progression.

RESULTS

In this study, by integrated genomics and experimental analysis, we identified GATA2 as a prevalent metastasis-associated genomic amplification to transcriptionally augment its own expression in PCa. Functional experiments demonstrated that GATA2 physically interacted and cooperated with SMAD4 for genome-wide chromatin co-occupancy and co-regulation of PCa genes and metastasis pathways like TGFβ signaling. Mechanistically, GATA2 was cooperative with SMAD4 to enhance TGFβ and AR signaling pathways, and activated the expression of TGFβ1 via directly binding to a distal enhancer of TGFβ1. Strinkingly, GATA2 and SMAD4 globally mediated inherited PCa risk and formed a transcriptional complex with HOXB13 at the PCa risk-associated rs339331/6q22 enhancer, leading to increased expression of the PCa susceptibility gene RFX6.

CONCLUSIONS

Our study prioritizes causal genomic amplification genes with prognostic values in PCa and reveals the pivotal roles of GATA2 in transcriptionally activating the expression of its own and TGFβ1, thereby co-opting to TGFβ1/SMAD4 signaling and RFX6 at 6q22 to modulate PCa predisposition and progression.

Multiple Germline Events Contribute to Cancer Development in Patients with Li-Fraumeni Syndrome

Li-Fraumeni syndrome (LFS) is an autosomal dominant cancer-predisposition disorder. Approximately 70% of individuals who fit the clinical definition of LFS harbor a pathogenic germline variant in the TP53 tumor suppressor gene. However, the remaining 30% of patients lack a TP53 variant and even among variant TP53 carriers, approximately 20% remain cancer-free. Understanding the variable cancer penetrance and phenotypic variability in LFS is critical to developing rational approaches to accurate, early tumor detection and risk-reduction strategies. We leveraged family-based whole-genome sequencing and DNA methylation to evaluate the germline genomes of a large, multi-institutional cohort of patients with LFS (n = 396) with variant (n = 374) or wildtype TP53 (n = 22). We identified alternative cancer-associated genetic aberrations in 8/14 wildtype TP53 carriers who developed cancer. Among variant TP53 carriers, 19/49 who developed cancer harbored a pathogenic variant in another cancer gene. Modifier variants in the WNT signaling pathway were associated with decreased cancer incidence. Furthermore, we leveraged the noncoding genome and methylome to identify inherited epimutations in genes including ASXL1, ETV6, and LEF1 that confer increased cancer risk. Using these epimutations, we built a machine learning model that can predict cancer risk in patients with LFS with an area under the receiver operator characteristic curve (AUROC) of 0.725 (0.633-0.810).

Significance

Our study clarifies the genomic basis for the phenotypic variability in LFS and highlights the immense benefits of expanding genetic and epigenetic testing of patients with LFS beyond TP53. More broadly, it necessitates the dissociation of hereditary cancer syndromes as single gene disorders and emphasizes the importance of understanding these diseases in a holistic manner as opposed to through the lens of a single gene.

Associations of CTCF and FOXA1 with androgen and IGF pathways in men with localized prostate cancer

AIMS

To examine associations between the transcription factors CCCTC-binding factor (CTCF) and forkhead box protein A1 (FOXA1) and the androgen receptor (AR) and their association with components of the insulin-like growth factor (IGF)-pathway in a cohort of men with localized prostate cancer.

METHODS

Using prostate tissue samples collected during the Prostate cancer: Evidence of Exercise and Nutrition Trial (PrEvENT) trial (N = 70 to 92, depending on section availability), we assessed the abundance of CTCF, FOXA1, AR, IGFIR, p-mTOR, PTEN and IGFBP-2 proteins using a modified version of the Allred scoring system. Validation studies were performed using large, publicly available datasets (TCGA) (N = 489).

RESULTS

We identified a strong correlation between CTCF and AR staining with benign prostate tissue. CTCF also strongly associated with the IGFIR, with PTEN and with phospho-mTOR. FOXA1 was also correlated with staining for the IGF-IR, with IGFBP-2 and with staining for activated phosphor-mTOR. The staining for the IGF-IR was strongly correlated with the AR.

CONCLUSION

Our findings emphasise the close and complex links between the endocrine controls, well known to play an important role in prostate cancer, and the transcription factors implicated by the recent genetic evidence.

Dissecting the immune suppressive human prostate tumor microenvironment via integrated single-cell and spatial transcriptomic analyses

The treatment of low-risk primary prostate cancer entails active surveillance only, while high-risk disease requires multimodal treatment including surgery, radiation therapy, and hormonal therapy. Recurrence and development of metastatic disease remains a clinical problem, without a clear understanding of what drives immune escape and tumor progression. Here, we comprehensively describe the tumor microenvironment of localized prostate cancer in comparison with adjacent normal samples and healthy controls. Single-cell RNA sequencing and high-resolution spatial transcriptomic analyses reveal tumor context dependent changes in gene expression. Our data indicate that an immune suppressive tumor microenvironment associates with suppressive myeloid populations and exhausted T-cells, in addition to high stromal angiogenic activity. We infer cell-to-cell relationships from high throughput ligand-receptor interaction measurements within undissociated tissue sections. Our work thus provides a highly detailed and comprehensive resource of the prostate tumor microenvironment as well as tumor-stromal cell interactions.

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.

Exploiting Dependence of Castration-Resistant Prostate Cancer on the Arginine Vasopressin Signaling Axis by Repurposing Vaptans

Men with advanced prostate cancer are treated by androgen deprivation therapy but the disease recurs as incurable castration-resistant prostate cancer (CRPC), requiring new treatment options. We previously demonstrated that the G protein-coupled receptor (GPCR) arginine vasopressin receptor type1A (AVPR1A) is expressed in CRPC and promotes castration-resistant growth in vitro and in vivo. AVPR1A is part of a family of GPCR's including arginine vasopressin receptor type 2 (AVPR2). Interrogation of prostate cancer patient sample data revealed that coexpression of AVPR1A and AVPR2 is highly correlated with disease progression. Stimulation of AVPR2 with a selective agonist desmopressin promoted CRPC cell proliferation through cAMP/protein kinase A signaling, consistent with AVPR2 coupling to the G protein subunit alpha s. In contrast, blocking AVPR2 with a selective FDA-approved antagonist, tolvaptan, reduced cell growth. In CRPC xenografts, antagonizing AVPR2, AVPR1A, or both significantly reduced CRPC tumor growth as well as decreased on-target markers of tumor burden. Combinatorial use of AVPR1A and AVPR2 antagonists promoted apoptosis synergistically in CRPC cells. Furthermore, we found that castration-resistant cells produced AVP, the endogenous ligand for arginine vasopressin receptors, and knockout of AVP in CRPC cells significantly reduced proliferation suggesting possible AVP autocrine signaling. These data indicate that the AVP/arginine vasopressin receptor signaling axis represents a promising and clinically actionable target for CRPC.

IMPLICATIONS

The arginine vasopressin signaling axis in CRPC provides a therapeutic window that is targetable through repurposing safe and effective AVPR1A and AVPR2 antagonists.

Genome-Wide 3'-UTR Single Nucleotide Polymorphism Association Study Identifies Significant Prostate Cancer Risk-Associated Functional Loci at 8p21.2 in Chinese Population

MicroRNAs (miRNAs) are involved in the regulation of gene expression via incomplete base pairing to sequence motifs at the three prime untranslated regions (3'-UTRs) of mRNAs and play critical roles in the etiology of cancers. Single nucleotide polymorphisms (SNPs) in the 3'-UTR miRNA-binding regions may influence the miRNA affinity. However, this biological mechanism in prostate cancer (PCa) remains unclear. Here, a three-stage genome-wide association study of 3'-UTR SNPs (n=33 117) is performed in 5515 Chinese men. Three genome-wide significant variants are discovered at 8p21.2 (rs1567669, rs4872176, and rs4872177), which are all located in a linkage disequilibrium region of the NKX3-1 gene. Phenome-wide association analysis using the FinnGen data reveals a specific association of rs1567669 with PCa over 2,264 disease endpoints. Expression quantitative trait locus analyses based on both Chinese PCa cohort and the GTEx database show that risk alleles of these SNPs are significantly associated with low expression of NKX3-1. Based on the MirSNP database, dual-luciferase reporter assays show that risk alleles of these SNPs downregulate the expression of NKX3-1 via increased miRNA binding. These results indicate that the SNPs at the 3'-UTR of NKX3-1 significantly downregulate NKX3-1 expression by influencing the affinity of miRNA and increase the PCa risk.

Disassembly of α6β4-mediated hemidesmosomal adhesions promotes tumorigenesis in PTEN-negative prostate cancer by targeting plectin to focal adhesions

Loss of α6β4-dependent hemidesmosomal adhesions has been observed during prostate cancer progression. However, the significance and underlying mechanisms by which aberrant hemidesmosome assembly may modulate tumorigenesis remain elusive. Using an extensive CRISPR/Cas9-mediated genetic engineering approaches in different prostate cancer cell lines combined with in vivo tumorigenesis studies in mice, bone marrow-on-chip assays and bioinformatics, as well as histological analysis of prostate cancer patient cohorts, we demonstrated that simultaneous loss of PTEN and hemidesmosomal adhesions induced several tumorigenic properties including proliferation, migration, resistance to anoikis, apoptosis, and drug treatment in vitro, and increased metastatic capacity in vivo. These effects were plectin-depended and plectin was associated with actin-rich adhesions upon hemidesmosome disruption in PTEN-negative prostate cancer cells leading to activation of EGFR/PI3K/Akt- and FAK/Src-pathways. These results suggest that analysis of PTEN and hemidesmosomal proteins may have diagnostic value helping to stratify prostate cancer patients with high risk for development of aggressive disease and highlight actin-associated plectin as a potential therapeutic target specifically in PTEN/hemidesmosome dual-negative prostate cancer.

A Novel Angiogenesis-Related Gene Signature to Predict Biochemical Recurrence of Patients with Prostate Cancer following Radical Therapy

Background

Prostate cancer (PCa) is one of the most common malignancies in males; we aim to establish a novel angiogenesis-related gene signature for biochemical recurrence (BCR) prediction in PCa patients following radical therapy.

Methods

Gene expression profiles and corresponding clinicopathological data were acquired from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database. We quantified the levels of various cancer hallmarks and identified angiogenesis as the primary risk factor for BCR. Then machine learning methods combined with Cox regression analysis were used to screen prognostic genes and construct an angiogenesis-related gene signature, which was further verified in external cohorts. Furthermore, estimation of immune cell abundance and prediction of drug responses were also conducted to detect potential mechanisms.

Results

Angiogenesis was regarded as the leading risk factor for BCR in PCa patients (HR = 1.58, 95% CI: 1.38–1.81), and a novel prognostic signature based on three genes (NRP1, JAG2, and VCAN) was developed in the training cohort and successfully validated in another three independent cohorts. In all datasets, this signature was identified as a prognostic factor with promising and robust predictive values. Moreover, it also predicted higher infiltration of regulatory T cells and M2-polarized macrophages and increased drug sensitivity of angiogenesis inhibitors in high-risk patients.

Conclusions

The angiogenesis-related three-gene signature may serve as an independent prognostic factor for BCR, which would contribute to risk stratification and personalized management of PCa patients after radical therapy in clinical practice.

Functional genomic analysis delineates regulatory mechanisms of GWAS-identified bipolar disorder risk variants

Background

Genome-wide association studies (GWASs) have identified multiple risk loci for bipolar disorder (BD). However, pinpointing functional (or causal) variants in the reported risk loci and elucidating their regulatory mechanisms remain challenging.

Methods

We first integrated chromatin immunoprecipitation sequencing (ChIP-Seq) data from human brain tissues (or neuronal cell lines) and position weight matrix (PWM) data to identify functional single-nucleotide polymorphisms (SNPs). Then, we verified the regulatory effects of these transcription factor (TF) binding–disrupting SNPs (hereafter referred to as “functional SNPs”) through a series of experiments, including reporter gene assays, allele-specific expression (ASE) analysis, TF knockdown, CRISPR/Cas9-mediated genome editing, and expression quantitative trait loci (eQTL) analysis. Finally, we overexpressed PACS1 (whose expression was most significantly associated with the identified functional SNPs rs10896081 and rs3862386) in mouse primary cortical neurons to investigate if PACS1 affects dendritic spine density.

Results

We identified 16 functional SNPs (in 9 risk loci); these functional SNPs disrupted the binding of 7 TFs, for example, CTCF and REST binding was frequently disrupted. We then identified the potential target genes whose expression in the human brain was regulated by these functional SNPs through eQTL analysis. Of note, we showed dysregulation of some target genes of the identified TF binding–disrupting SNPs in BD patients compared with controls, and overexpression of PACS1 reduced the density of dendritic spines, revealing the possible biological mechanisms of these functional SNPs in BD.

Conclusions

Our study identifies functional SNPs in some reported risk loci and sheds light on the regulatory mechanisms of BD risk variants. Further functional characterization and mechanistic studies of these functional SNPs and candidate genes will help to elucidate BD pathogenesis and develop new therapeutic approaches and drugs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-022-01057-3.

HOXB13 suppresses de novo lipogenesis through HDAC3-mediated epigenetic reprogramming in prostate cancer

HOXB13, a homeodomain transcription factor, critically regulates androgen receptor (AR) activities and androgen-dependent prostate cancer (PCa) growth. However, its functions in AR-independent contexts remain elusive. Here we report HOXB13 interaction with histone deacetylase HDAC3, which is disrupted by the HOXB13 G84E mutation that has been associated with early-onset PCa. Independently of AR, HOXB13 recruits HDAC3 to lipogenic enhancers to catalyze histone deacetylation and suppress lipogenic regulators such as fatty acid synthase. Analysis of human tissues reveals that the HOXB13 gene is hypermethylated and downregulated in approximately 30% of metastatic castration-resistant PCa. HOXB13 loss or G84E mutation leads to lipid accumulation in PCa cells, thereby promoting cell motility and xenograft tumor metastasis, which is mitigated by pharmaceutical inhibition of fatty acid synthase. In summary, we present evidence that HOXB13 recruits HDAC3 to suppress de novo lipogenesis and inhibit tumor metastasis and that lipogenic pathway inhibitors may be useful to treat HOXB13-low PCa.

Regulatory variants at 2q33.1 confer schizophrenia risk by modulating distal gene TYW5 expression

Genome-wide association studies have shown that genetic variants at 2q33.1 are strongly associated with schizophrenia. However, potential causal variants in this locus and their roles in schizophrenia remain unknown. Here, we identified two functional variants (rs796364 and rs281759) that disrupt CTCF, RAD21 and FOXP2 binding at 2q33.1. We systematically investigated the regulatory mechanisms of these two variants with serial experiments, including reporter gene assays and electrophoretic mobility shift assay. Intriguingly, these two single nucleotide polymorphisms physically interacted with TYW5 and showed the most significant associations with TYW5 expression in human brain. Consistently, CRISPR-Cas9-mediated genome editing confirmed the regulatory effect of the two single nucleotide polymorphisms on TYW5 expression. Additionally, expression analysis indicated that TYW5 was significantly upregulated in brains of schizophrenia cases compared with controls, suggesting that rs796364 and rs281759 might confer schizophrenia risk by modulating TYW5 expression. We over-expressed TYW5 in mouse neural stem cells and rat primary neurons to mimic its upregulation in schizophrenia and found significant alterations in the proliferation and differentiation of neural stem cells, as well as dendritic spine density following TYW5 overexpression, indicating its important roles in neurodevelopment and spine morphogenesis. Furthermore, we independently confirmed the association between rs796364 and schizophrenia in a Chinese cohort of 8202 subjects. Finally, transcriptome analysis revealed that TYW5 affected schizophrenia-associated pathways. These lines of evidence consistently revealed that rs796364 and rs281759 might contribute to schizophrenia risk by regulating the expression of TYW5, a gene whose expression dysregulation affects two important schizophrenia pathophysiological processes (i.e. neurodevelopment and dendritic spine formation).

Functional genomics elucidates regulatory mechanisms of Parkinson's disease-associated variants

BACKGROUND

Genome-wide association studies (GWASs) have identified multiple risk loci for Parkinson's disease (PD). However, identifying the functional (or potential causal) variants in the reported risk loci and elucidating their roles in PD pathogenesis remain major challenges. To identify the potential causal (or functional) variants in the reported PD risk loci and to elucidate their regulatory mechanisms, we report a functional genomics study of PD.

METHODS

We first integrated chromatin immunoprecipitation sequencing (ChIP-Seq) (from neuronal cells and human brain tissues) data and GWAS-identified single-nucleotide polymorphisms (SNPs) in PD risk loci. We then conducted a series of experiments and analyses to validate the regulatory effects of these (i.e., functional) SNPs, including reporter gene assays, allele-specific expression (ASE), transcription factor (TF) knockdown, CRISPR-Cas9-mediated genome editing, and expression quantitative trait loci (eQTL) analysis.

RESULTS

We identified 44 SNPs (from 11 risk loci) affecting the binding of 12 TFs and we validated the regulatory effects of 15 TF binding-disrupting SNPs. In addition, we also identified the potential target genes regulated by these TF binding-disrupting SNPs through eQTL analysis. Finally, we showed that 4 eQTL genes of these TF binding-disrupting SNPs were dysregulated in PD cases compared with controls.

CONCLUSION

Our study systematically reveals the gene regulatory mechanisms of PD risk variants (including widespread disruption of CTCF binding), generates the landscape of potential PD causal variants, and pinpoints promising candidate genes for further functional characterization and drug development.

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.

Large Multicohort Study Reveals a Prostate Cancer Susceptibility Allele at 5p15 Regulating TERT via Androgen Signaling-Orchestrated Chromatin Binding of E2F1 and MYC

Aberrant telomerase reverse transcriptase (TERT) expression is crucial for tumor survival and cancer cells escaping apoptosis. Multiple TERT-locus variants at 5p15 have been discovered in association with cancer risk, yet the underlying mechanisms and clinical impacts remain unclear. Here, our association studies showed that the TERT promoter variant rs2853669 confers a risk of prostate cancer (PCa) in different ethnic groups. Further functional investigation revealed that the allele-specific binding of MYC and E2F1 at TERT promoter variant rs2853669 associates with elevated level of TERT in PCa. Mechanistically, androgen stimulations promoted the binding of MYC to allele T of rs2853669, thereby activating TERT, whereas hormone deprivations enhanced E2F1 binding at allele C of rs2853669, thus upregulating TERT expression. Notably, E2F1 could cooperate with AR signaling to regulate MYC expression. Clinical data demonstrated synergistic effects of MYC/E2F1/TERT expression or with the TT and CC genotype of rs2853669 on PCa prognosis and severity. Strikingly, single-nucleotide editing assays showed that the CC genotype of rs2853669 obviously promotes epithelial-mesenchymal transition (EMT) and the development of castration-resistant PCa (CRPC), confirmed by unbiased global transcriptome profiling. Our findings thus provided compelling evidence for understanding the roles of noncoding variations coordinated with androgen signaling and oncogenic transcription factors in mis-regulating TERT expression and driving PCa.

Allele-specific epigenetic activity in prostate cancer and normal prostate tissue implicates prostate cancer risk mechanisms

Genome-wide association studies (GWASs) of prostate cancer have identified >250 significant risk loci, but the causal variants and mechanisms for these loci remain largely unknown. Here, we sought to identify and characterize risk-harboring regulatory elements by integrating epigenomes from primary prostate tumor and normal tissues of 27 individuals across the H3K27ac, H3K4me3, and H3K4me2 histone marks and FOXA1 and HOXB13 transcription factors. We identified 7,371 peaks with significant allele specificity (allele-specific quantitative trait locus [asQTL] peaks). Showcasing their relevance to prostate cancer risk, H3K27ac T-asQTL peaks were the single annotation most enriched for prostate cancer GWAS heritability (40×), significantly higher than corresponding non-asQTL H3K27ac peaks (14×) or coding regions (14×). Surprisingly, fine-mapped GWAS risk variants were most significantly enriched for asQTL peaks observed in tumors, including asQTL peaks that were differentially imbalanced with respect to tumor-normal states. These data pinpointed putative causal regulatory elements at 20 GWAS loci, of which 11 were detected only in the tumor samples. More broadly, tumor-specific asQTLs were enriched for expression QTLs in benign tissues as well as accessible regions found in stem cells, supporting a hypothesis where some germline variants become reactivated during or after transformation and can be captured by epigenomic profiling of the tumor. Our study demonstrates the power of allele specificity in chromatin signals to uncover GWAS mechanisms, highlights the relevance of tumor-specific regulation in the context of cancer risk, and prioritizes multiple loci for experimental follow-up.

Functional Screenings Identify Regulatory Variants Associated with Breast Cancer Susceptibility

Genome-wide association studies (GWAS) have identified more than 2000 single nucleotide polymorphisms (SNPs) associated with breast cancer susceptibility, most of which are located in the non-coding region. However, the causal SNPs functioning as gene regulatory elements still remain largely undisclosed. Here, we applied a Dinucleotide Parallel Reporter sequencing (DiR-seq) assay to evaluate 288 breast cancer risk SNPs in nine different breast cancer cell lines. Further multi-omics analysis with the ATAC-seq (Assay for Transposase-Accessible Chromatin using sequencing), DNase-seq (DNase I hypersensitive sites sequencing) and histone modification ChIP-seq (Chromatin Immunoprecipitation sequencing) nominated seven functional SNPs in breast cancer cells. Functional investigations show that rs4808611 affects breast cancer progression by altering the gene expression of NR2F6. For the other site, rs2236007, the alteration promotes the binding of the suppressive transcription factor EGR1 and results in the downregulation of PAX9 expression. The downregulated expression of PAX9 causes cancer malignancies and is associated with the poor prognosis of breast cancer patients. Our findings contribute to defining the functional risk SNPs and the related genes for breast cancer risk prediction.

Parallel Reporter Assays Identify Altered Regulatory Role of rs684232 in Leading to Prostate Cancer Predisposition

Functional characterization of cancer risk-associated single nucleotide polymorphism (SNP) identified by genome-wide association studies (GWAS) has become a big challenge. To identify the regulatory risk SNPs that can lead to transcriptional misregulation, we performed parallel reporter gene assays with both alleles of 213 prostate cancer risk-associated GWAS SNPs in 22Rv1 cells. We disclosed 32 regulatory SNPs that exhibited different regulatory activities with two alleles. For one of the regulatory SNPs, rs684232, we found that the variation altered chromatin binding of transcription factor FOXA1 on the DNA region and led to aberrant gene expression of VPS53, FAM57A, and GEMIN4, which play vital roles in prostate cancer malignancy. Our findings reveal the roles and underlying mechanism of rs684232 in prostate cancer progression and hold great promise in benefiting prostate cancer patients with prognostic prediction and target therapies.

DCBLD1 is associated with the integrin signaling pathway and has prognostic value in non-small cell lung and invasive breast carcinoma

Germline single nucleotide polymorphisms in the promoter region of the DCBLD1 gene are associated with non-smoking cases of both non-small cell lung carcinoma (NSCLC) and human papillomavirus-negative head and neck cancer. However the clinical relevance and function of DCBLD1 remain unclear. This multicenter retrospective study was designed to evaluate the prognostic value and function of DCBLD1 in the four main solid cancers: NSCLC, invasive breast carcinoma, colorectal adenocarcinoma and prostate adenocarcinoma. We included the following cohorts: GSE81089 NSCLC, METABRIC invasive breast carcinoma, GSE14333 colorectal adenocarcinoma, GSE70770 prostate adenocarcinoma and The Cancer Genome Atlas (TCGA) Firehose Legacy cohorts of all four cancers. DCBLD1 gene expression was associated with a worse overall survival in multivariate analyses for both NSCLC cohorts (TCGA: P = 0.03 and GSE81089: P = 0.04) and both invasive breast carcinoma cohorts (TCGA: P = 0.02 and METABRIC: P < 0.001). Patients with high DCBLD1 expression showed an upregulation of the integrin signaling pathway in comparison to those with low DCBLD1 expression in the TCGA NSCLC cohort (FDR = 5.16 × 10^-14) and TCGA invasive breast carcinoma cohort (FDR = 1.94 × 10^-05).

Genetic Variants in Transcription Factor Binding Sites in Humans: Triggered by Natural Selection and Triggers of Diseases

Variants of transcription factor binding sites (TFBSs) constitute an important part of the human genome. Current evidence demonstrates close links between nucleotides within TFBSs and gene expression. There are multiple pathways through which genomic sequences located in TFBSs regulate gene expression, and recent genome-wide association studies have shown the biological significance of TFBS variation in human phenotypes. However, numerous challenges remain in the study of TFBS polymorphisms. This article aims to cover the current state of understanding as regards the genomic features of TFBSs and TFBS variants; the mechanisms through which TFBS variants regulate gene expression; the approaches to studying the effects of nucleotide changes that create or disrupt TFBSs; the challenges faced in studies of TFBS sequence variations; the effects of natural selection on collections of TFBSs; in addition to the insights gained from the study of TFBS alleles related to gout, its associated comorbidities (increased body mass index, chronic kidney disease, diabetes, dyslipidemia, coronary artery disease, ischemic heart disease, hypertension, hyperuricemia, osteoporosis, and prostate cancer), and the treatment responses of patients.

Identification of Germline Genetic Variants that Increase Prostate Cancer Risk and Influence Development of Aggressive Disease

Prostate cancer (PrCa) is a heterogeneous disease, which presents in individual patients across a diverse phenotypic spectrum ranging from indolent to fatal forms. No robust biomarkers are currently available to enable routine screening for PrCa or to distinguish clinically significant forms, therefore late stage identification of advanced disease and overdiagnosis plus overtreatment of insignificant disease both remain areas of concern in healthcare provision. PrCa has a substantial heritable component, and technological advances since the completion of the Human Genome Project have facilitated improved identification of inherited genetic factors influencing susceptibility to development of the disease within families and populations. These genetic markers hold promise to enable improved understanding of the biological mechanisms underpinning PrCa development, facilitate genetically informed PrCa screening programmes and guide appropriate treatment provision. However, insight remains largely lacking regarding many aspects of their manifestation; especially in relation to genes associated with aggressive phenotypes, risk factors in non-European populations and appropriate approaches to enable accurate stratification of higher and lower risk individuals. This review discusses the methodology used in the elucidation of genetic loci, genes and individual causal variants responsible for modulating PrCa susceptibility; the current state of understanding of the allelic spectrum contributing to PrCa risk; and prospective future translational applications of these discoveries in the developing eras of genomics and personalised medicine.

The landscape of RNA polymerase II-associated chromatin interactions in prostate cancer

Transcriptional dysregulation is a hallmark of prostate cancer (PCa). We mapped the RNA polymerase II-associated (RNA Pol II-associated) chromatin interactions in normal prostate cells and PCa cells. We discovered thousands of enhancer-promoter, enhancer-enhancer, as well as promoter-promoter chromatin interactions. These transcriptional hubs operate within the framework set by structural proteins - CTCF and cohesins - and are regulated by the cooperative action of master transcription factors, such as the androgen receptor (AR) and FOXA1. By combining analyses from metastatic castration-resistant PCa (mCRPC) specimens, we show that AR locus amplification contributes to the transcriptional upregulation of the AR gene by increasing the total number of chromatin interaction modules comprising the AR gene and its distal enhancer. We deconvoluted the transcription control modules of several PCa genes, notably the biomarker KLK3, lineage-restricted genes (KRT8, KRT18, HOXB13, FOXA1, ZBTB16), the drug target EZH2, and the oncogene MYC. By integrating clinical PCa data, we defined a germline-somatic interplay between the PCa risk allele rs684232 and the somatically acquired TMPRSS2-ERG gene fusion in the transcriptional regulation of multiple target genes - VPS53, FAM57A, and GEMIN4. Our studies implicate changes in genome organization as a critical determinant of aberrant transcriptional regulation in PCa.

Using machine learning method to identify MYLK as a novel marker to predict biochemical recurrence in prostate cancer

Aim: This study aims to identify novel marker to predict biochemical recurrence (BCR) in prostate cancer patients after radical prostatectomy with negative surgical margin. Materials & methods: The Cancer Genome Atlas database, Gene Expression Omnibus database and Cancer Cell Line Encyclopedia database were employed. The ensemble support vector machine-recursive feature elimination method was performed to select crucial gene for BCR. Results: We identified MYLK as a novel and independent biomarker for BCR in The Cancer Genome Atlas training cohort and confirmed in four independent Gene Expression Omnibus validation cohorts. Multi-omic analysis suggested that MYLK was a DNA methylation-driven gene. Additionally, MYLK had significant positive correlations with immune infiltrations. Conclusion:MYLK was identified and validated as a novel, robust and independent biomarker for BCR in prostate cancer.

Validation of prostate cancer risk variants rs10993994 and rs7098889 by CRISPR/Cas9 mediated genome editing

GWAS have identified numerous SNPs associated with prostate cancer risk. One such SNP is rs10993994. It is located in the β-microseminoprotein (MSMB) promoter region, mediates MSMB prostate secretion levels, and is linked to mRNA expression changes in both MSMB and the adjacent gene NCOA4. In addition, our previous work showed a second SNP, rs7098889, is in positive linkage disequilibrium with rs10993994 and associated with MSMB expression independent of rs10993994. Here, we generate a series of clones with single alleles removed by double guide RNA (gRNA) mediated CRISPR/Cas9 deletions, through which we demonstrate that each of these SNPs independently and greatly alters MSMB expression in an allele-specific manner. We further show that these SNPs have no substantial effect on the expression of NCOA4. These data demonstrate that a single SNP can have a large effect on gene expression and illustrate the importance of functional validation studies to deconvolute observed correlations. The method we have developed is generally applicable to test any SNP for which a relevant heterozygous cell line is available. AUTHOR SUMMARY: In pursuing the underlying biological mechanism of prostate cancer pathogenesis, scientists utilized the existence of common single nucleotide polymorphisms (SNPs) in the human genome as genetic markers to perform large scale genome wide association studies (GWAS) and have so far identified more than a hundred prostate cancer risk variants. Such variants provide an unbiased and systematic new venue to study the disease mechanism, and the next big challenge is to translate these genetic associations to the causal role of altered gene function in oncogenesis. The majority of these variants are waiting to be studied and lots of them may act in oncogenesis through gene expression regulation. To prove the concept, we took rs10993994 and its linked rs7098889 as an example and engineered single cell clones by allelic-specific CRISPR/Cas9 deletion to separate the effect of each allele. We observed that a single nucleotide difference would lead to surprisingly high level of MSMB gene expression change in a gene specific and cell-type specific manner. Our study strongly supports the notion that differential level of gene expression caused by risk variants and their associated genetic locus play a major role in oncogenesis and also highlights the importance of studying the function of MSMB encoded β-MSP in prostate cancer pathogenesis.

Complement C7 (C7), a Potential Tumor Suppressor, Is an Immune-Related Prognostic Biomarker in Prostate Cancer (PC)

Objectives: Prostate cancer (PC) is the second most frequent tumor in men, which has a high recurrence rate and poor prognosis. Therefore, this study aimed to identify novel prognostic biomarkers and therapeutic targets for immunotherapy and small molecule drugs for PC treatment.

Materials and Methods: The Estimation of Stromal and Immune cells in Malignant Tumor tissues using Expression data (ESTIMATE) algorithm was applied to calculate immune scores and stromal scores of TCGA-PRAD data. Differentially expressed genes (DEGs) were identified using R package “limma.” GO, KEGG, and DO analyses were performed to analyze DEGs. Overall survival and disease-free survival analyses were conducted for hub gene identification. To validate the hub gene at the mRNA and protein expression levels, genetic alterations were measured, and CCLE and Cox regression analyses were performed. Connectivity map (CMap) analysis and GSEA were performed for drug exploration and function analysis, respectively.

Results: Immune scores ranged from −1795.98 to 2339.39, and stomal scores ranged from −1877.60 to 1659.96. In total, 45 tumor microenvironment (TME)-related DEGs were identified, of which Complement C7 (C7) was selected and validated as a hub gene. CMap analysis identified six small molecule drugs as potential agents for PC treatment. Further analysis demonstrated that C7 expression was significantly correlated with clinical T, pathological N, and immune infiltration level.

Conclusions: In conclusion, of the 45 TME-related DEGs, C7 was shown to correlate with PC prognosis in patients, indicating it as a novel prognostic biomarker and immunotherapy target in PC. Additionally, six small molecule drugs showed strong therapeutic potential for PC treatment.

Pioneer of prostate cancer: past, present and the future of FOXA1

Prostate cancer is the most commonly diagnosed non-cutaneous cancers in North American men. While androgen deprivation has remained as the cornerstone of prostate cancer treatment, resistance ensues leading to lethal disease. Forkhead box A1 (FOXA1) encodes a pioneer factor that induces open chromatin conformation to allow the binding of other transcription factors. Through direct interactions with the Androgen Receptor (AR), FOXA1 helps to shape AR signaling that drives the growth and survival of normal prostate and prostate cancer cells. FOXA1 also possesses an AR-independent role of regulating epithelial-to-mesenchymal transition (EMT). In prostate cancer, mutations converge onto the coding sequence and cis-regulatory elements (CREs) of FOXA1, leading to functional alterations. In addition, FOXA1 activity in prostate cancer can be modulated post-translationally through various mechanisms such as LSD1-mediated protein demethylation. In this review, we describe the latest discoveries related to the function and regulation of FOXA1 in prostate cancer, pointing to their relevance to guide future clinical interventions.

Regulatory mechanisms of major depressive disorder risk variants

Major depressive disorder (MDD) is one of the most prevalent psychiatric disorders and a leading cause of disability worldwide. Though recent genome-wide association studies (GWAS) have identified multiple risk variants for MDD, how these variants confer MDD risk remains largely unknown. Here we systematically characterize the regulatory mechanism of MDD risk variants using a functional genomics approach. By integrating chromatin immunoprecipitation sequencing (ChIP-Seq) (from human brain tissues or neuronal cells) and position weight matrix (PWM) data, we identified 34 MDD risk SNPs that disrupt the binding of 15 transcription factors (TFs). We verified the regulatory effect of the TF binding-disrupting SNPs with reporter gene assays, allelic-specific expression analysis, and CRISPR-Cas9-mediated genome editing. Expression quantitative trait loci (eQTL) analysis identified the target genes that might be regulated by these regulatory risk SNPs. Finally, we found that NEGR1 (regulated by the TF binding-disrupting MDD risk SNP rs3101339) was dysregulated in the brains of MDD cases compared with controls, implying that rs3101339 may confer MDD risk by affecting NEGR1 expression. Our findings reveal how genetic variants contribute to MDD risk by affecting TF binding and gene regulation. More importantly, our study identifies the potential MDD causal variants and their target genes, thus providing pivotal candidates for future mechanistic study and drug development.

Hereditary Predisposition to Prostate Cancer: From Genetics to Clinical Implications

Prostate cancer (PrCa) ranks among the top five cancers for both incidence and mortality worldwide. A significant proportion of PrCa susceptibility has been attributed to inherited predisposition, with 10-20% of cases expected to occur in a hereditary/familial context. Advances in DNA sequencing technologies have uncovered several moderate- to high-penetrance PrCa susceptibility genes, most of which have previously been related to known hereditary cancer syndromes, namely the hereditary breast and ovarian cancer (BRCA1, BRCA2, ATM, CHEK2, and PALB2) and Lynch syndrome (MLH1, MSH2, MSH6, and PMS2) genes. Additional candidate genes have also been suggested, but further evidence is needed to include them in routine genetic testing. Recommendations based on clinical features, family history, and ethnicity have been established for more cost-efficient genetic testing of patients and families who may be at an increased risk of developing PrCa. The identification of alterations in PrCa predisposing genes may help to inform screening strategies, as well as treatment options, in the metastatic setting. This review provides an overview of the genetic basis underlying hereditary predisposition to PrCa, the current genetic screening recommendations, and the implications for clinical management of the disease.

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.

Annotations capturing cell type-specific TF binding explain a large fraction of disease heritability

Regulatory variation plays a major role in complex disease and that cell type-specific binding of transcription factors (TF) is critical to gene regulation. However, assessing the contribution of genetic variation in TF-binding sites to disease heritability is challenging, as binding is often cell type-specific and annotations from directly measured TF binding are not currently available for most cell type-TF pairs. We investigate approaches to annotate TF binding, including directly measured chromatin data and sequence-based predictions. We find that TF-binding annotations constructed by intersecting sequence-based TF-binding predictions with cell type-specific chromatin data explain a large fraction of heritability across a broad set of diseases and corresponding cell types; this strategy of constructing annotations addresses both the limitation that identical sequences may be bound or unbound depending on surrounding chromatin context and the limitation that sequence-based predictions are generally not cell type-specific. We partitioned the heritability of 49 diseases and complex traits using stratified linkage disequilibrium (LD) score regression with the baseline-LD model (which is not cell type-specific) plus the new annotations. We determined that 100 bp windows around MotifMap sequenced-based TF-binding predictions intersected with a union of six cell type-specific chromatin marks (imputed using ChromImpute) performed best, with an 58% increase in heritability enrichment compared to the chromatin marks alone (11.6× vs. 7.3×, P = 9 × 10-14 for difference) and a 20% increase in cell type-specific signal conditional on annotations from the baseline-LD model (P = 8 × 10-11 for difference). Our results show that TF-binding annotations explain substantial disease heritability and can help refine genome-wide association signals.

Androgen Receptor and Poly(ADP-ribose) Glycohydrolase Inhibition Increases Efficiency of Androgen Ablation in Prostate Cancer Cells

There is mounting evidence of androgen receptor signaling inducing genome instability and changing DNA repair capacity in prostate cancer cells. Expression of genes associated with base excision repair (BER) is increased with prostate cancer progression and correlates with poor prognosis. Poly(ADP-ribose) polymerase (PARP) and poly(ADP-ribose) glycohydrolase (PARG) are key enzymes in BER that elongate and degrade PAR polymers on target proteins. While PARP inhibitors have been tested in clinical trials and are a promising therapy for prostate cancer patients with TMPRSS2-ERG fusions and mutations in DNA repair genes, PARG inhibitors have not been evaluated. We show that PARG is a direct androgen receptor (AR) target gene. AR is recruited to the PARG locus and induces PARG expression. Androgen ablation combined with PARG inhibition synergistically reduces BER capacity in independently derived LNCaP and LAPC4 prostate cancer cell lines. A combination of PARG inhibition with androgen ablation or with the DNA damaging drug, temozolomide, significantly reduces cellular proliferation and increases DNA damage. PARG inhibition alters AR transcriptional output without changing AR protein levels. Thus, AR and PARG are engaged in reciprocal regulation suggesting that the success of androgen ablation therapy can be enhanced by PARG inhibition in prostate cancer patients.

DNA-PK, Nuclear mTOR, and the Androgen Pathway in Prostate Cancer

Androgen and its receptor (AR) are major drivers of prostate cancer (PCa), a leading cause of mortality in aging men. Thus, understanding the numerous mechanisms by which AR can promote the growth and proliferation of PCa cells and enable their escape from hormone-dependent therapies, eventually leading to metastasis and death of the patient, is essential to discover alternative therapeutic approaches. Recently, two structurally related members of the phosphatidylinositol 3-kinase-like protein kinase (PIKK) family, DNA-dependent protein kinase (DNA-PK) and mammalian target of rapamycin (mTOR), were shown to have a direct role in modulating AR activity on chromatin of PCa cells. In this review, the common features of DNA-PK and mTOR and the similarities in their noncanonical roles as transcription coregulators of the AR are highlighted. An outlook on how these findings could be translated into new approaches to manage and treat PCa is provided.

Androgen Receptor Dependence

Androgens and the androgen receptor (AR) play crucial roles in the biology of normal and diseased prostate tissue, including prostate cancer (PCa). This dependence is evidenced by the use of androgen depletion therapy (ADT) as the primary treatment for locally advanced, metastatic, or relapsed PCa. This dependence is further evidenced by the various mechanisms employed by PCa cells to re-activate the AR to circumvent the growth-inhibitory effects of ADT. Re-activation of the AR during ADT is central to the disease evolving into the lethal castration resistant PCa (CRPC) phenotype, which is responsible for nearly all PCa mortality. Thus, understanding the regulation of AR and AR signaling is important for understanding the development and progression of PCa. This understanding provides the foundation for development of newer approaches for targeting CRPC therapeutically.

Androgen receptor suppresses prostate cancer metastasis but promotes bladder cancer metastasis via differentially altering miRNA525-5p/SLPI-mediated vasculogenic mimicry formation

Early studies suggest that the androgen receptor (AR) may play differential roles in influencing prostate cancer (PCa) and bladder cancer (BCa) metastasis, but the underlying mechanisms remain unclear. Here, we found that the AR might function via differentially altering vasculogenic mimicry (VM) formation to either decrease PCa metastasis or increase BCa metastasis. Mechanism dissection showed that the AR could differentially alter the expression of the VM marker SLPI through miR-525-5p to regulate SLPI; moreover, it could either increase miR-525-5p transcription in PCa or decrease it in BCa via binding to different androgen-response-elements (AREs) located at different positions in the miR-525 precursor promoter. Further, results from liquid chromatography-mass spectrometry (LC-MS) showed that the co-factors of AR in PCa and BCa are NFIX and HDAC2, respectively. Together, these results provide the first detailed mechanism of how the AR can differentially alter PCa and BCa metastasis; thus, targeting the newly identified AR-miR-525-5p-SLPI axis may help suppress metastasis.

Expression of CCCTC-binding factor (CTCF) is linked to poor prognosis in prostate cancer

The chromatin-organizing factor CCCTC-binding factor (CTCF) is involved in transcriptional regulation, DNA-loop formation, and telomere maintenance. To evaluate the clinical impact of CTCF in prostate cancer, we analyzed CTCF expression by immunohistochemistry on a tissue microarray containing 17 747 prostate cancers. Normal prostate tissue showed negative to low CTCF expression, while in prostate cancers, CTCF expression was seen in 7726 of our 12 555 (61.5%) tumors and was considered low in 44.6% and high in 17% of cancers. Particularly, high CTCF expression was significantly associated with the presence of the transmembrane protease, serine 2:ETS-related gene fusion: Only 10% of ERG-negative cancers, but 30% of ERG-positive cancers had high-level CTCF expression (P < 0.0001). CTCF expression was significantly associated with advanced pathological tumor stage, high Gleason grade (P < 0.0001 each), nodal metastasis (P = 0.0122), and early biochemical recurrence (P < 0.0001). Multivariable modeling revealed that the prognostic impact of CTCF was independent from established presurgical parameters such as clinical stage and Gleason grade of the biopsy. Comparison with key molecular alterations showed strong associations with the expression of the Ki-67 proliferation marker and presence of phosphatase and tensin homolog deletions (P < 0.0001 each). The results of our study identify CTCF expression as a candidate biomarker for prognosis assessment in prostate cancer.

Genome-wide germline correlates of the epigenetic landscape of prostate cancer

Oncogenesis is driven by germline, environmental and stochastic factors. It is unknown how these interact to produce the molecular phenotypes of tumors. We therefore quantified the influence of germline polymorphisms on the somatic epigenome of 589 localized prostate tumors. Predisposition risk loci influence a tumor's epigenome, uncovering a mechanism for cancer susceptibility. We identified and validated 1,178 loci associated with altered methylation in tumoral but not nonmalignant tissue. These tumor methylation quantitative trait loci influence chromatin structure, as well as RNA and protein abundance. One prominent tumor methylation quantitative trait locus is associated with AKT1 expression and is predictive of relapse after definitive local therapy in both discovery and validation cohorts. These data reveal intricate crosstalk between the germ line and the epigenome of primary tumors, which may help identify germline biomarkers of aggressive disease to aid patient triage and optimize the use of more invasive or expensive diagnostic assays.

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.

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.

Androgen receptor enhancer usage and the chromatin regulatory landscape in human prostate cancers

The androgen receptor (AR) is commonly known as a key transcription factor in prostate cancer development, progression and therapy resistance. Genome-wide chromatin association studies revealed that transcriptional regulation by AR mainly depends on binding to distal regulatory enhancer elements that control gene expression through chromatin looping to gene promoters. Changes in the chromatin epigenetic landscape and DNA sequence can locally alter AR-DNA-binding capacity and consequently impact transcriptional output and disease outcome. The vast majority of reports describing AR chromatin interactions have been limited to cell lines, identifying numerous other factors and interacting transcription factors that impact AR chromatin interactions. Do these factors also impact AR cistromics - the genome-wide chromatin-binding landscape of AR - in vivo? Recent technological advances now enable researchers to identify AR chromatin-binding sites and their target genes in human specimens. In this review, we provide an overview of the different factors that influence AR chromatin binding in prostate cancer specimens, which is complemented with knowledge from cell line studies. Finally, we discuss novel perspectives on studying AR cistromics in clinical samples.