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

Uncovering somatic genetic drivers in prostate cancer through comprehensive genome-wide analysis

Given that hereditary prostate cancer (PCa) accounts for only a small fraction of PCa phenotypes, there is still a substantial journey ahead in exploring the somatic genetic drivers contributing to sporadic PCa. The expression quantitative trait loci (eQTLs) data were sourced from the GTEx dataset for prostate-specific genes, and the summary statistic information was collected for 5854 genes. Genetic associations with PCa were extracted from three well-established consortiums: the UK Biobank (9131 cases and 173,493 controls), the PRACTICAL study (79,148 cases and 61,106 controls), and the FinnGen cohort (13,216 cases and 119,948 controls). To prioritize potential causal targets, additional analysis, including the protein-protein interaction (PPI), The Cancer Genome Atlas (TCGA) dataset, and the single-cell-type expression analysis, was performed. Generally, a total of 150 common significant genes with the same causal association with PCa were identified. Out of the 150 genes examined, 67.33% (101/150) were found to have protein-coding functions, while only 30.67% (46/150) of these genes had prior mentions in the scientific literature. Notably, the analysis of the TCGA dataset showed that only 44.67% (67/150) of the genes produced consistent results with the Mendelian randomization (MR) analysis. Furthermore, the evaluation of single-cell RNA-seq data and colocalization analysis singled out MSMB as a critical gene associated with the occurrence of PCa. We pinpointed a range of prostate-specific genes that display causal associations with the onset of PCa. Among these, the MSMB gene emerged as a pivotal factor linked to PCa, demonstrating robust consistency across all four assessments, including the MR, TCGA dataset, single-cell RNA-seq data, and colocalization analysis. These findings provided fresh perspectives on the pathogenesis of PCa and presented potential targets for drug development.

Enhancer transcription profiling reveals an enhancer RNA-driven ferroptosis and new therapeutic opportunities in prostate cancer

Enhancer RNAs (eRNAs), a subclass of non-coding RNAs transcribed from enhancer regions, have emerged as critical regulators of gene expression; however, their functional roles in prostate cancer remain largely unexplored. In this study, we performed integrated chromatin accessibility and transcriptomic analyses using ATAC-seq and RNA-seq on twenty pairs of prostate cancer and matched benign tissues. By incorporating chromatin immunoprecipitation sequencing data, we identified a subset of differentially expressed eRNAs significantly associated with genes involved in prostate development and oncogenic signaling pathways. Among these, lactotransferrin-eRNA (LTFe) was markedly downregulated in prostate cancer tissues, with functional analyses revealing its tumor-suppressive role. Mechanistically, LTFe promotes the transcription of its target gene, lactotransferrin (LTF), by interacting with heterogeneous nuclear ribonucleoprotein F (HNRNPF) and facilitating enhancer-promoter chromatin interactions. Furthermore, we demonstrate that the LTFe-LTF axis facilitates ferroptosis by modulating iron transport. Notably, androgen receptor (AR) signaling disrupts LTFe-associated chromatin looping, leading to ferroptosis resistance. Therapeutically, co- administration of the AR inhibitor enzalutamide and the ferroptosis inducer RSL3 significantly suppressed tumor growth, offering a promising strategy for castration-resistant prostate cancer. Collectively, this study provides novel insights into the mechanistic role of eRNAs in prostate cancer, highlighting the LTFe-LTF axis as a critical epigenetic regulator and potential therapeutic target for improved treatment outcomes.

The role and mechanism of hepatocyte nuclear factor 1β in the occurrence and development of different human tumors: A pan-cancer analysis

Carcinomatosis is one of the leading threats to human public fitness. HNF1B is a critical transcription element in vertebrate proliferation and oncogenesis, which has been shown to play roles in reactive oxygen species (ROS) metabolism. Our previous results have identified HNF1B as a tumor suppressor that could inhibit the malignant progression of prostate cancer. Yet there is no pan-carcinomatosis analysis of HNF1B, which could help us better understand common and unique underlying mechanisms in mankind knubs to enhance novel and competent treatment. Here, in our research, we evaluated the utterance pattern and explored the function of HNF1B in 33 knub categories using the data from the Cancer Genome Atlas Program (TCGA), Gene Expression Omnibus (GEO), and CLNICAL PROTEOMICTUMOR ANALYSIS CONSORTIUM (CPTAC) dataset. We found different HNF1B roles in various cancer types. HNF1B was upregulated in CHOL, STAD, KIRP, and THCA, and was downregulated in GBM, KICH, COAD, KIRC, LUSC, SARC, PAAD, and TGCT. Prognostic analyses indicated that higher HNF1B displayed better illness outcomes in BLCA, READ, and PRAD, while poorer outcomes in LUSC and THYM. HNF1B mutation was most frequent in endometrial cancer but was not associated with disease prognosis. It was discovered that HNF1B utterance relevant to endothelial cell penetration status in BLCA, ESCA, LUAD, LUSC, and TGCT, and carcinomatosis-associated fibroblast infiltration was observed in ESCA, KIRC, LIHC, and TGCT. Moreover, functional enrichment analysis disclosed that metabolism-related functions were implicated in the function of HNF1B. Taken together, our pan- carcinomatosis analysis showed the complicated roles of HNF1B in a variety of carcinomatoses, being able to improve the extensive comprehension of HNF1B's role in tumorigenesis.

HNF1B polymorphisms and endometrial cancer risk: validation of identified loci and evaluation of novel variants

BACKGROUND

This study aimed to validate HNF1B single nucleotide polymorphisms (SNPs) associated with endometrial cancer risk in a Chinese Han population and explore novel SNPs. Our findings enhance the understanding of genetic components and are crucial for detection strategies and personalized medicine.

METHODS

We genotyped four HNF1B SNPs in 637 patients and 667 controls using Agena MassARRAY. Logistic regression calculated odds ratios (ORs) and 95% CI. Forest plots visualize stratified analyses. Multiple comparisons tested genetic loci-clinical indicator associations.

RESULTS

The study confirmed that rs4430796 (A>G) reduced endometrial cancer risk (OR = 0.83, 95% CI: 0.70-0.99, p = 0.041). Additionally, novel HNF1B mutations were associated with endometrial cancer risk: rs7405776 in individuals under the age of 55 (OR = 047, 95% CI: 0.25-0.91, p = 0.025) and nonsmokers (OR = 0.42, 95% CI: 0.23-0.75, p = 0.004), and rs11651755 in drinkers (OR = 0.39, 95% CI: 0.17-0.90, p = 0.027) and nonsmokers (OR = 0.43, 95% CI: 0.23-0.81, p = 0.009). The SNP rs4430796 was also associated with the CA125 level (p < 0.05).

CONCLUSION

HNF1B polymorphisms influence endometrial cancer risk in the Chinese Han population. Further studies are needed to explore the functional roles and clinical practicality of these SNPs.

Integrated analysis of single-cell RNA-seq and bulk RNA-seq revealed key genes for bone metastasis and chemoresistance in prostate cancer

BACKGROUND

Prostate cancer (PCa) is a serious malignancy. The main causes of PCa aggravation and death are unexplained resistance to chemotherapy and bone metastases.

OBJECTIVE

This study aimed to investigate the molecular mechanisms associated with the dynamic processes of progression, bone metastasis, and chemoresistance in PCa.

METHODS

Through comprehensive analysis of single-cell RNA sequencing (scRNA-seq) data, Gene Expression Omnibus (GEO) tumor progression and metastasis-related genes were identified. These genes were subjected to lasso regression modeling using the Cancer Genome Atlas (TCGA) database. Tartrate-resistant acid phosphatase (TRAP) staining and real-time quantitative PCR (RT-qPCR) were used to evaluate osteoclast differentiation. CellMiner was used to confirm the effect of LDHA on chemoresistance. Finally, the relationship between LDHA and chemoresistance was verified using doxorubicin-resistant PCa cell lines.

RESULTS

7928 genes were identified as genes related to tumor progression and metastasis. Of these, 7 genes were found to be associated with PCa prognosis. The scRNA-seq and TCGA data showed that the expression of LDHA was higher in tumors and associated with poor prognosis of PCa. In addition, upregulation of LDHA in PCa cells induces osteoclast differentiation. Additionally, high LDHA expression was associated with resistance to Epirubicin, Elliptinium acetate, and doxorubicin. Cellular experiments demonstrated that LDHA knockdown inhibited doxorubicin resistance in PCa cells.

CONCLUSIONS

LDHA may play a potential contributory role in PCa initiation and development, bone metastasis, and chemoresistance. LDHA is a key target for the treatment of PCa.

Oncogene goosecoid is transcriptionally regulated by E2F1 and correlates with disease progression in prostate cancer

BACKGROUND

Although some well-established oncogenes are involved in cancer initiation and progression such as prostate cancer (PCa), the long tail of cancer genes remains to be defined. Goosecoid ( GSC ) has been implicated in cancer development. However, the comprehensive biological role of GSC in pan-cancer, specifically in PCa, remains unexplored. The aim of this study was to investigate the role of GSC in PCa development.

METHODS

We performed a systematic bioinformatics exploration of GSC using datasets from The Cancer Genome Atlas, Genotype-Tissue Expression, Gene Expression Omnibus, German Cancer Research Center, and our in-house cohorts. First, we evaluated the expression of GSC and its association with patient prognosis, and identified GSC -relevant genetic alterations in cancers. Further, we focused on the clinical characterization and prognostic analysis of GSC in PCa. To understand the transcriptional regulation of GSC by E2F transcription factor 1 ( E2F1 ), we performed chromatin immunoprecipitation quantitative polymerase chain reaction (qPCR). Functional experiments were conducted to validate the effect of GSC on the tumor cellular phenotype and sensitivity to trametinib.

RESULTS

GSC expression was elevated in various tumors and significantly correlated with patient prognosis. The alterations of GSC contribute to the progression of various tumors especially in PCa. Patients with PCa and high GSC expression exhibited worse progression-free survival and biochemical recurrence outcomes. Further, GSC upregulation in patients with PCa was mostly accompanied with higher Gleason score, advanced tumor stage, lymph node metastasis, and elevated prostate-specific antigen (PSA) levels. Mechanistically, the transcription factor, E2F1 , stimulates GSC by binding to its promoter region. Detailed experiments further demonstrated that GSC acted as an oncogene and influenced the response of PCa cells to trametinib treatment.

CONCLUSIONS

GSC was highly overexpressed and strongly correlated with patient prognosis in PCa. We found that GSC , regulated by E2F1 , acted as an oncogene and impeded the therapeutic efficacy of trametinib in PCa.

Combined SNPs sequencing and allele specific proteomics capture reveal functional causality underpinning the 2p25 prostate cancer susceptibility locus

Genome wide association studies (GWASs) have identified numerous risk loci associated with prostate cancer, yet unraveling their functional significance remains elusive. Leveraging our high-throughput SNPs-seq method, we pinpointed rs4519489 within the multi-ancestry GWAS-discovered 2p25 locus as a potential functional SNP due to its significant allelic differences in protein binding. Here, we conduct a comprehensive analysis of rs4519489 and its associated gene, NOL10, employing diverse cohort data and experimental models. Clinical findings reveal a synergistic effect between rs4519489 genotype and NOL10 expression on prostate cancer prognosis and severity. Through unbiased proteomics screening, we reveal that the risk allele A of rs4519489 exhibits enhanced binding to USF1, a novel oncogenic transcription factor (TF) implicated in prostate cancer progression and prognosis, resulting in elevated NOL10 expression. Furthermore, we elucidate that NOL10 regulates cell cycle pathways, fostering prostate cancer progression. The concurrent expression of NOL10 and USF1 correlates with aggressive prostate cancer characteristics and poorer prognosis. Collectively, our study offers a robust strategy for functional SNP screening and TF identification through high-throughput SNPs-seq and unbiased proteomics, highlighting the rs4519489-USF1-NOL10 regulatory axis as a promising biomarker or therapeutic target for clinical diagnosis and treatment of prostate cancer.

LncRNA VPS9D1-AS1 regulates miR-187-3p/fibroblast growth factor receptor-like 1 axis to promote proliferation, migration, and invasion of prostate cancer cells

The morbidity and mortality of prostate cancer are increasing year by year, and the survival rate of prostate cancer patients after treatment is low. Therefore, investigating the molecular mechanism underlying prostate cancer is crucial for developing effective treatments. Recent studies have shown the important role of long-chain non-coding RNAs (lncRNAs) in tumorigenesis. VPS9D1-AS1 can modulate the progression of multiple cancers, but its molecular action mechanism in prostate cancer remains unknown. This study, therefore, intended to investigate the regulatory mechanism of VPS9D1-AS1 in prostate cancer. First, differentially expressed lncRNAs in prostate cancer were identified through bioinformatics approaches. The target lncRNA for the study was determined by reviewing the relevant literature and its downstream miRNA/mRNA axis was uncovered. Then, quantitative reverse transcription polymerase chain reaction was introduced to assess the expression of VPS9D1-AS1, miR-187-3p, and fibroblast growth factor receptor-like 1 (FGFRL1) at a cellular level, and Western blot was conducted to assess the protein level of FGFRL1 in cells. The results indicated that VPS9D1-AS1 and FGFRL1 were highly expressed in prostate cancer while miR-187-3p was less expressed. Besides, MTT, colony formation, wound healing, and cell invasion assays showed that silencing VPS9D1-AS1 inhibited the viability, migration ability, and invasion ability of prostate cancer cells. Dual-luciferase assay and RNA binding protein immunoprecipitation assay were performed to explore the interplay of miR-187-3p and VPS9D1-AS1 or FGFRL1. The results showed that VPS9D1-AS1 could sponge miR-187-3p, and FGFRL1 could serve as a direct target of miR-187-3p. Moreover, combined with the results of the rescue experiment, VPS9D1-AS1 was found to upregulate FGFRL1 by competitively sponging miR-187-3p to accelerate the malignant behaviors of prostate cancer cells. In conclusion, VPS9D1-AS1 could promote the phenotype progression of prostate cancer cells through targeting the miR-187-3p/FGFRL1 axis, and it has the potential to be a target for prostate cancer patients.

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.

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.

Fusion sequencing via terminator-assisted synthesis (FTAS-seq) identifies TMPRSS2 fusion partners in prostate cancer

Genetic rearrangements that fuse an androgen-regulated promoter area with a protein-coding portion of an originally androgen-unaffected gene are frequent in prostate cancer, with the fusion between transmembrane serine protease 2 (TMPRSS2) and ETS transcription factor ERG (ERG) (TMPRSS2-ERG fusion) being the most prevalent. Conventional hybridization- or amplification-based methods can test for the presence of expected gene fusions, but the exploratory analysis of currently unknown fusion partners is often cost-prohibitive. Here, we developed an innovative next-generation sequencing (NGS)-based approach for gene fusion analysis termed fusion sequencing via terminator-assisted synthesis (FTAS-seq). FTAS-seq can be used to enrich the gene of interest while simultaneously profiling the whole spectrum of its 3'-terminal fusion partners. Using this novel semi-targeted RNA-sequencing technique, we were able to identify 11 previously uncharacterized TMPRSS2 fusion partners and capture a range of TMPRSS2-ERG isoforms. We tested the performance of FTAS-seq with well-characterized prostate cancer cell lines and utilized the technique for the analysis of patient RNA samples. FTAS-seq chemistry combined with appropriate primer panels holds great potential as a tool for biomarker discovery that can support the development of personalized cancer therapies.

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.