Metastatic castration-resistant prostate cancer (mCRPC) is the lethal form of the disease. Several recent studies have identified genomic alterations in mCRPC, but the clinical implications of these ...genomic alterations have not been fully elucidated.
To use whole-genome sequencing (WGS) to assess the association between key driver gene alterations and overall survival (OS), and to use whole-transcriptome RNA sequencing to identify genomic drivers of enzalutamide resistance.
We performed survival analyses and gene set enrichment analysis (GSEA) on WGS and RNA sequencing results for a cohort of 101 mCRPC patients.
OS was the clinical endpoint for all univariate and multivariable survival analyses. Candidate drivers of enzalutamide resistance were identified in an unbiased manner, and mutations of the top candidate were further assessed for enrichment among enzalutamide-resistant patients using Fisher's exact test.
Harboring two DNA alterations in RB1 was independently predictive of poor OS (median 14.1 vs 42.0mo; p=0.007) for men with mCRPC. GSEA identified the Wnt/β-catenin pathway as the top differentially modulated pathway among enzalutamide-resistant patients. Furthermore, β-catenin mutations were exclusive to enzalutamide-resistant patients (p=0.01) and independently predictive of poor OS (median 13.6 vs 41.7mo; p=0.025).
The presence of two RB1 DNA alterations identified in our WGS analysis was independently associated with poor OS among men with mCRPC. The Wnt/β-catenin pathway plays an important role in enzalutamide resistance, with differential pathway expression and enrichment of β-catenin mutations in enzalutamide-resistant patients. Moreover, β-catenin mutations were predictive of poor OS in our cohort.
We observed a correlation between genomic findings for biopsy samples from metastases from men with metastatic castration-resistant prostate cancer (mCRPC) and clinical outcomes. This work sheds new light on clinically relevant genomic alterations in mCRPC and provides a roadmap for the development of new personalized treatment regimens in mCRPC.
The presence of two DNA alterations in RB1 is associated with poor overall survival independently of other clinicopathologic factors in metastatic castration-resistant prostate cancer. In addition, Wnt/β-catenin pathway activation and β-catenin mutations are associated with enzalutamide resistance and poor overall survival.
Prostate cancer is the second leading cause of male cancer deaths. Castration-resistant prostate cancer (CRPC) is a lethal stage of the disease that emerges when endocrine therapies are no longer ...effective at suppressing activity of the androgen receptor (AR) transcription factor. The purpose of this study was to identify genomic mechanisms that contribute to the development and progression of CRPC.
We used whole-genome and targeted DNA-sequencing approaches to identify mechanisms underlying CRPC in an aggregate cohort of 272 prostate cancer patients. We analyzed structural rearrangements at the genome-wide level and carried out a detailed structural rearrangement analysis of the
locus. We used genome engineering to perform experimental modeling of
gene rearrangements and long-read RNA sequencing to analyze effects on expression of AR and truncated AR variants (AR-V).
was among the most frequently rearranged genes in CRPC tumors.
gene rearrangements promoted expression of diverse AR-V species.
gene rearrangements occurring in the context of
amplification correlated with AR overexpression. Cell lines with experimentally derived
gene rearrangements displayed high expression of tumor-specific AR-Vs and were resistant to endocrine therapies, including the AR antagonist enzalutamide.
gene rearrangements are an important mechanism of resistance to endocrine therapies in CRPC.
MicroRNAs (miRNAs) are small, non-coding RNAs that play essential roles in plant growth, development, and stress response. We conducted a genome-wide survey of maize miRNA genes, characterizing their ...structure, expression, and evolution. Computational approaches based on homology and secondary structure modeling identified 150 high-confidence genes within 26 miRNA families. For 25 families, expression was verified by deep-sequencing of small RNA libraries that were prepared from an assortment of maize tissues. PCR-RACE amplification of 68 miRNA transcript precursors, representing 18 families conserved across several plant species, showed that splice variation and the use of alternative transcriptional start and stop sites is common within this class of genes. Comparison of sequence variation data from diverse maize inbred lines versus teosinte accessions suggest that the mature miRNAs are under strong purifying selection while the flanking sequences evolve equivalently to other genes. Since maize is derived from an ancient tetraploid, the effect of whole-genome duplication on miRNA evolution was examined. We found that, like protein-coding genes, duplicated miRNA genes underwent extensive gene-loss, with approximately 35% of ancestral sites retained as duplicate homoeologous miRNA genes. This number is higher than that observed with protein-coding genes. A search for putative miRNA targets indicated bias towards genes in regulatory and metabolic pathways. As maize is one of the principal models for plant growth and development, this study will serve as a foundation for future research into the functional roles of miRNA genes.
Noncoding RNAs (ncRNAs) are emerging as key molecules in human cancer, with the potential to serve as novel markers of disease and to reveal uncharacterized aspects of tumor biology. Here we discover ...121 unannotated prostate cancer-associated ncRNA transcripts (PCATs) by ab initio assembly of high-throughput sequencing of polyA(+) RNA (RNA-Seq) from a cohort of 102 prostate tissues and cells lines. We characterized one ncRNA, PCAT-1, as a prostate-specific regulator of cell proliferation and show that it is a target of the Polycomb Repressive Complex 2 (PRC2). We further found that patterns of PCAT-1 and PRC2 expression stratified patient tissues into molecular subtypes distinguished by expression signatures of PCAT-1-repressed target genes. Taken together, our findings suggest that PCAT-1 is a transcriptional repressor implicated in a subset of prostate cancer patients. These findings establish the utility of RNA-Seq to identify disease-associated ncRNAs that may improve the stratification of cancer subtypes.
Although recurrent gene fusions involving erythroblastosis virus E26 transformation-specific (ETS) family transcription factors are common in prostate cancer, their products are considered ...'undruggable' by conventional approaches. Recently, rare targetable gene fusions involving the anaplastic lymphoma receptor tyrosine kinase (ALK) gene, have been identified in 1-5% of lung cancers, suggesting that similar rare gene fusions may occur in other common epithelial cancers, including prostate cancer. Here we used paired-end transcriptome sequencing to screen ETS rearrangement-negative prostate cancers for targetable gene fusions and identified the SLC45A3-BRAF (solute carrier family 45, member 3-v-raf murine sarcoma viral oncogene homolog B1) and ESRP1-RAF1 (epithelial splicing regulatory protein-1-v-raf-1 murine leukemia viral oncogene homolog-1) gene fusions. Expression of SLC45A3-BRAF or ESRP1-RAF1 in prostate cells induced a neoplastic phenotype that was sensitive to RAF and mitogen-activated protein kinase kinase (MAP2K1) inhibitors. Screening a large cohort of patients, we found that, although rare, recurrent rearrangements in the RAF pathway tend to occur in advanced prostate cancers, gastric cancers and melanoma. Taken together, our results emphasize the key role of RAF family gene rearrangements in cancer, suggest that RAF and MEK inhibitors may be useful in a subset of gene fusion-harboring solid tumors and demonstrate that sequencing of tumor transcriptomes and genomes may lead to the identification of rare targetable fusions across cancer types.
While next-generation sequencing (NGS) has become the primary technology for discovering gene fusions, we are still faced with the challenge of ensuring that causative mutations are not missed while ...minimizing false positives. Currently, there are many computational tools that predict structural variations (SV) and gene fusions using whole genome (WGS) and transcriptome sequencing (RNA-seq) data separately. However, as both WGS and RNA-seq have their limitations when used independently, we hypothesize that the orthogonal validation from integrating both data could generate a sensitive and specific approach for detecting high-confidence gene fusion predictions. Fortunately, decreasing NGS costs have resulted in a growing quantity of patients with both data available. Therefore, we developed a gene fusion discovery tool, INTEGRATE, that leverages both RNA-seq and WGS data to reconstruct gene fusion junctions and genomic breakpoints by split-read mapping. To evaluate INTEGRATE, we compared it with eight additional gene fusion discovery tools using the well-characterized breast cell line HCC1395 and peripheral blood lymphocytes derived from the same patient (HCC1395BL). The predictions subsequently underwent a targeted validation leading to the discovery of 131 novel fusions in addition to the seven previously reported fusions. Overall, INTEGRATE only missed six out of the 138 validated fusions and had the highest accuracy of the nine tools evaluated. Additionally, we applied INTEGRATE to 62 breast cancer patients from The Cancer Genome Atlas (TCGA) and found multiple recurrent gene fusions including a subset involving estrogen receptor. Taken together, INTEGRATE is a highly sensitive and accurate tool that is freely available for academic use.
Snail1 is a zinc finger transcriptional repressor whose pathological expression has been linked to cancer cell epithelial-mesenchymal transition (EMT) programs and the induction of tissue-invasive ...activity, but pro-oncogenic events capable of regulating Snail1 activity remain largely uncharacterized. Herein, we demonstrate that p53 loss-of-function or mutation promotes cancer cell EMT by de-repressing Snail1 protein expression and activity. In the absence of wild-type p53 function, Snail1-dependent EMT is activated in colon, breast, and lung carcinoma cells as a consequence of a decrease in miRNA-34 levels, which suppress Snail1 activity by binding to highly conserved 3' untranslated regions in Snail1 itself as well as those of key Snail1 regulatory molecules, including β-catenin, LEF1, and Axin2. Although p53 activity can impact cell cycle regulation, apoptosis, and DNA repair pathways, the EMT and invasion programs initiated by p53 loss of function or mutation are completely dependent on Snail1 expression. These results identify a new link between p53, miR-34, and Snail1 in the regulation of cancer cell EMT programs.
Polycomb Repressive Complexes (PRC1 and PRC2)-mediated epigenetic regulation is critical for maintaining cellular homeostasis. Members of Polycomb Group (PcG) proteins including EZH2, a PRC2 ...component, are upregulated in various cancer types, implicating their role in tumorigenesis. Here, we have identified several microRNAs (miRNAs) that are repressed by EZH2. These miRNAs, in turn, regulate the expression of PRC1 proteins BMI1 and RING2. We found that ectopic overexpression of EZH2-regulated miRNAs attenuated cancer cell growth and invasiveness, and abrogated cancer stem cell properties. Importantly, expression analysis revealed an inverse correlation between miRNA and PRC protein levels in cell culture and prostate cancer tissues. Taken together, our data have uncovered a coordinate regulation of PRC1 and PRC2 activities that is mediated by miRNAs.
► We identified miRNAs that are negatively regulated by PRC2 protein EZH2 ► EZH2-regulated miRNAs, in turn, regulate PRC1 proteins BMI1 and RING2 ► miRNAs and PRC protein levels are inversely correlated in prostate cancer ► PRC1 and PRC2 activities are coordinately regulated via miRNAs
Significance Paramyxoviruses, the cause of many important human and animal diseases, constitute a large family of enveloped negative-stranded RNA viruses including parainfluenza virus 5 (PIV5). The ...virion RNA is associated with ∼2,600 protomers of N-protein in the form of a helical ribonucleoprotein (RNP) (nucleocapsid). The RNP serves as the template for the viral polymerase in vivo. When expressed, N forms a 13 member N-ring that resembles the building block of the RNP. We have determined the atomic structure of the N-ring from PIV5 with 78 bound RNA residues. Precisely, six nucleotides of RNA are associated with each N protomer. Modeling PIV5-N in an “open conformation” in the ring structure indicates how transcription/replication could occur with minimal changes to the nucleocapsid structure.
Parainfluenza virus 5 (PIV5) is a member of the Paramyxoviridae family of membrane-enveloped viruses with a negative-sense RNA genome that is packaged and protected by long filamentous nucleocapsid-helix structures (RNPs). These RNPs, consisting of ∼2,600 protomers of nucleocapsid (N) protein, form the template for viral transcription and replication. We have determined the 3D X-ray crystal structure of the nucleoprotein (N)-RNA complex from PIV5 to 3.11-Å resolution. The structure reveals a 13-mer nucleocapsid ring whose diameter, cavity, and pitch/height dimensions agree with EM data from early studies on the Paramyxovirinae subfamily of native RNPs, indicating that it closely represents one-turn in the building block of the RNP helices. The PIV5-N nucleocapsid ring encapsidates a nuclease resistant 78-nt RNA strand in its positively charged groove formed between the N-terminal (NTD) and C-terminal (CTD) domains of its successive N protomers. Six nucleotides precisely are associated with each N protomer, with alternating three-base-in three-base-out conformation. The binding of six nucleotides per protomer is consistent with the “rule of six” that governs the genome packaging of the Paramyxovirinae subfamily of viruses. PIV5-N protomer subdomains are very similar in structure to the previously solved Nipah-N structure, but with a difference in the angle between NTD/CTD at the RNA hinge region. Based on the Nipah-N structure we modeled a PIV5-N open conformation in which the CTD rotates away from the RNA strand into the inner spacious nucleocapsid-ring cavity. This rotation would expose the RNA for the viral polymerase activity without major disruption of the nucleocapsid structure.
Enhancer of zeste homolog 2 (EZH2) is a mammalian histone methyltransferase that contributes to the epigenetic silencing of target genes and regulates the survival and metastasis of cancer cells. ...EZH2 is overexpressed in aggressive solid tumors by mechanisms that remain unclear. Here we show that the expression and function of EZH2 in cancer cell lines are inhibited by microRNA-101 (miR-101). Analysis of human prostate tumors revealed that miR-101 expression decreases during cancer progression, paralleling an increase in EZH2 expression. One or both of the two genomic loci encoding miR-101 were somatically lost in 37.5% of clinically localized prostate cancer cells (6 of 16) and 66.7% of metastatic disease cells (22 of 33). We propose that the genomic loss of miR-101 in cancer leads to overexpression of EZH2 and concomitant dysregulation of epigenetic pathways, resulting in cancer progression.