The Landscape of Circular RNA in Cancer.

Circular RNAs (circRNAs) are an intriguing class of RNA due to their covalently closed structure, high stability, and implicated roles in gene regulation. Here, we used an exome capture RNA sequencing protocol to detect and characterize circRNAs across >2,000 cancer samples. When compared against Ribo-Zero and RNase R, capture sequencing significantly enhanced the enrichment of circRNAs and preserved accurate circular-to-linear ratios. Using capture sequencing, we built the most comprehensive catalog of circRNA species to date: MiOncoCirc, the first database to be composed primarily of circRNAs directly detected in tumor tissues. Using MiOncoCirc, we identified candidate circRNAs to serve as biomarkers for prostate cancer and were able to detect circRNAs in urine. We further detected a novel class of circular transcripts, termed read-through circRNAs, that involved exons originating from different genes. MiOncoCirc will serve as a valuable resource for the development of circRNAs as diagnostic or therapeutic targets across cancer types.

Inactivation of CDK12 Delineates a Distinct Immunogenic Class of Advanced Prostate Cancer.

Using integrative genomic analysis of 360 metastatic castration-resistant prostate cancer (mCRPC) samples, we identified a novel subtype of prostate cancer typified by biallelic loss of CDK12 that is mutually exclusive with tumors driven by DNA repair deficiency, ETS fusions, and SPOP mutations. CDK12 loss is enriched in mCRPC relative to clinically localized disease and characterized by focal tandem duplications (FTDs) that lead to increased gene fusions and marked differential gene expression. FTDs associated with CDK12 loss result in highly recurrent gains at loci of genes involved in the cell cycle and DNA replication. CDK12 mutant cases are baseline diploid and do not exhibit DNA mutational signatures linked to defects in homologous recombination. CDK12 mutant cases are associated with elevated neoantigen burden ensuing from fusion-induced chimeric open reading frames and increased tumor T cell infiltration/clonal expansion. CDK12 inactivation thereby defines a distinct class of mCRPC that may benefit from immune checkpoint immunotherapy.

Oncogenic Role of THOR, a Conserved Cancer/Testis Long Non-coding RNA.

Large-scale transcriptome sequencing efforts have vastly expanded the catalog of long non-coding RNAs (lncRNAs) with varying evolutionary conservation, lineage expression, and cancer specificity. Here, we functionally characterize a novel ultraconserved lncRNA, THOR (ENSG00000226856), which exhibits expression exclusively in testis and a broad range of human cancers. THOR knockdown and overexpression in multiple cell lines and animal models alters cell or tumor growth supporting an oncogenic role. We discovered a conserved interaction of THOR with IGF2BP1 and show that THOR contributes to the mRNA stabilization activities of IGF2BP1. Notably, transgenic THOR knockout produced fertilization defects in zebrafish and also conferred a resistance to melanoma onset. Likewise, ectopic expression of human THOR in zebrafish accelerated the onset of melanoma. THOR represents a novel class of functionally important cancer/testis lncRNAs whose structure and function have undergone positive evolutionary selection.

Integrative clinical genomics of advanced prostate cancer.

Toward development of a precision medicine framework for metastatic, castration-resistant prostate cancer (mCRPC), we established a multi-institutional clinical sequencing infrastructure to conduct prospective whole-exome and transcriptome sequencing of bone or soft tissue tumor biopsies from a cohort of 150 mCRPC affected individuals. Aberrations of AR, ETS genes, TP53, and PTEN were frequent (40%-60% of cases), with TP53 and AR alterations enriched in mCRPC compared to primary prostate cancer. We identified new genomic alterations in PIK3CA/B, R-spondin, BRAF/RAF1, APC, ß-catenin, and ZBTB16/PLZF. Moreover, aberrations of BRCA2, BRCA1, and ATM were observed at substantially higher frequencies (19.3% overall) compared to those in primary prostate cancers. 89% of affected individuals harbored a clinically actionable aberration, including 62.7% with aberrations in AR, 65% in other cancer-related genes, and 8% with actionable pathogenic germline alterations. This cohort study provides clinically actionable information that could impact treatment decisions for these affected individuals.

Targeting SWI/SNF ATPases in enhancer-addicted prostate cancer.

The switch/sucrose non-fermentable (SWI/SNF) complex has a crucial role in chromatin remodelling1 and is altered in over 20% of cancers2,3. Here we developed a proteolysis-targeting chimera (PROTAC) degrader of the SWI/SNF ATPase subunits, SMARCA2 and SMARCA4, called AU-15330. Androgen receptor (AR)+ forkhead box A1 (FOXA1)+ prostate cancer cells are exquisitely sensitive to dual SMARCA2 and SMARCA4 degradation relative to normal and other cancer cell lines. SWI/SNF ATPase degradation rapidly compacts cis-regulatory elements bound by transcription factors that drive prostate cancer cell proliferation, namely AR, FOXA1, ERG and MYC, which dislodges them from chromatin, disables their core enhancer circuitry, and abolishes the downstream oncogenic gene programs. SWI/SNF ATPase degradation also disrupts super-enhancer and promoter looping interactions that wire supra-physiologic expression of the AR, FOXA1 and MYC oncogenes themselves. AU-15330 induces potent inhibition of tumour growth in xenograft models of prostate cancer and synergizes with the AR antagonist enzalutamide, even inducing disease remission in castration-resistant prostate cancer (CRPC) models without toxicity. Thus, impeding SWI/SNF-mediated enhancer accessibility represents a promising therapeutic approach for enhancer-addicted cancers.

Expressed pseudogenes in the transcriptional landscape of human cancers.

Pseudogene transcripts can provide a novel tier of gene regulation through generation of endogenous siRNAs or miRNA-binding sites. Characterization of pseudogene expression, however, has remained confined to anecdotal observations due to analytical challenges posed by the extremely close sequence similarity with their counterpart coding genes. Here, we describe a systematic analysis of pseudogene "transcription" from an RNA-Seq resource of 293 samples, representing 13 cancer and normal tissue types, and observe a surprisingly prevalent, genome-wide expression of pseudogenes that could be categorized as ubiquitously expressed or lineage and/or cancer specific. Further, we explore disease subtype specificity and functions of selected expressed pseudogenes. Taken together, we provide evidence that transcribed pseudogenes are a significant contributor to the transcriptional landscape of cells and are positioned to play significant roles in cellular differentiation and cancer progression, especially in light of the recently described ceRNA networks. Our work provides a transcriptome resource that enables high-throughput analyses of pseudogene expression.

Development of an orally bioavailable mSWI/SNF ATPase degrader and acquired mechanisms of resistance in prostate cancer.

Mammalian switch/sucrose nonfermentable (mSWI/SNF) ATPase degraders have been shown to be effective in enhancer-driven cancers by functioning to impede oncogenic transcription factor chromatin accessibility. Here, we developed AU-24118, an orally bioavailable proteolysis-targeting chimera (PROTAC) degrader of mSWI/SNF ATPases (SMARCA2 and SMARCA4) and PBRM1. AU-24118 demonstrated tumor regression in a model of castration-resistant prostate cancer (CRPC) which was further enhanced with combination enzalutamide treatment, a standard of care androgen receptor (AR) antagonist used in CRPC patients. Importantly, AU-24118 exhibited favorable pharmacokinetic profiles in preclinical analyses in mice and rats, and further toxicity testing in mice showed a favorable safety profile. As acquired resistance is common with targeted cancer therapeutics, experiments were designed to explore potential mechanisms of resistance that may arise with long-term mSWI/SNF ATPase PROTAC treatment. Prostate cancer cell lines exposed to long-term treatment with high doses of a mSWI/SNF ATPase degrader developed SMARCA4 bromodomain mutations and ABCB1 (ATP binding cassette subfamily B member 1) overexpression as acquired mechanisms of resistance. Intriguingly, while SMARCA4 mutations provided specific resistance to mSWI/SNF degraders, ABCB1 overexpression provided broader resistance to other potent PROTAC degraders targeting bromodomain-containing protein 4 and AR. The ABCB1 inhibitor, zosuquidar, reversed resistance to all three PROTAC degraders tested. Combined, these findings position mSWI/SNF degraders for clinical translation for patients with enhancer-driven cancers and define strategies to overcome resistance mechanisms that may arise.

Targeting SWI/SNF ATPases in H3.3K27M diffuse intrinsic pontine gliomas.

Diffuse midline gliomas (DMGs) including diffuse intrinsic pontine gliomas (DIPGs) bearing lysine-to-methionine mutations in histone H3 at lysine 27 (H3K27M) are lethal childhood brain cancers. These tumors harbor a global reduction in the transcriptional repressive mark H3K27me3 accompanied by an increase in the transcriptional activation mark H3K27ac. We postulated that H3K27M mutations, in addition to altering H3K27 modifications, reprogram the master chromatin remodeling switch/sucrose nonfermentable (SWI/SNF) complex. The SWI/SNF complex can exist in two main forms termed BAF and PBAF that play central roles in neurodevelopment and cancer. Moreover, BAF antagonizes PRC2, the main enzyme catalyzing H3K27me3. We demonstrate that H3K27M gliomas show increased protein levels of the SWI/SNF complex ATPase subunits SMARCA4 and SMARCA2, and the PBAF component PBRM1. Additionally, knockdown of mutant H3K27M lowered SMARCA4 protein levels. The proteolysis targeting chimera (PROTAC) AU-15330 that simultaneously targets SMARCA4, SMARCA2, and PBRM1 for degradation exhibits cytotoxicity in H3.3K27M but not H3 wild-type cells. AU-15330 lowered chromatin accessibility measured by ATAC-Seq at nonpromoter regions and reduced global H3K27ac levels. Integrated analysis of gene expression, proteomics, and chromatin accessibility in AU-15330-treated cells demonstrated reduction in the levels of FOXO1, a key member of the forkhead family of transcription factors. Moreover, genetic or pharmacologic targeting of FOXO1 resulted in cell death in H3K27M cells. Overall, our results suggest that H3K27M up-regulates SMARCA4 levels and combined targeting of SWI/SNF ATPases in H3.3K27M can serve as a potent therapeutic strategy for these deadly childhood brain tumors.

Targeting the lipid kinase PIKfyve upregulates surface expression of MHC class I to augment cancer immunotherapy.

Despite the remarkable clinical success of immunotherapies in a subset of cancer patients, many fail to respond to treatment and exhibit resistance. Here, we found that genetic or pharmacologic inhibition of the lipid kinase PIKfyve, a regulator of autophagic flux and lysosomal biogenesis, upregulated surface expression of major histocompatibility complex class I (MHC-I) in cancer cells via impairing autophagic flux, resulting in enhanced cancer cell killing mediated by CD8+ T cells. Genetic depletion or pharmacologic inhibition of PIKfyve elevated tumor-specific MHC-I surface expression, increased intratumoral functional CD8+ T cells, and slowed tumor progression in multiple syngeneic mouse models. Importantly, enhanced antitumor responses by Pikfyve-depletion were CD8+ T cell- and MHC-I-dependent, as CD8+ T cell depletion or B2m knockout rescued tumor growth. Furthermore, PIKfyve inhibition improved response to immune checkpoint blockade (ICB), adoptive cell therapy, and a therapeutic vaccine. High expression of PIKFYVE was also predictive of poor response to ICB and prognostic of poor survival in ICB-treated cohorts. Collectively, our findings show that targeting PIKfyve enhances immunotherapies by elevating surface expression of MHC-I in cancer cells, and PIKfyve inhibitors have potential as agents to increase immunotherapy response in cancer patients.

Proxalutamide reduces SARS-CoV-2 infection and associated inflammatory response.

Early in the COVID-19 pandemic, data suggested that males had a higher risk of developing severe disease and that androgen deprivation therapy might be associated with protection. Combined with the fact that TMPRSS2 (transmembrane serine protease 2), a host entry factor for the SARS-CoV-2 virus, was a well-known androgen-regulated gene, this led to an upsurge of research investigating androgen receptor (AR)-targeting drugs. Proxalutamide, an AR antagonist, was shown in initial clinical studies to benefit COVID-19 patients; however, further validation is needed as one study was retracted. Due to continued interest in proxalutamide, which is in phase 3 trials, we examined its ability to impact SARS-CoV-2 infection and downstream inflammatory responses. Proxalutamide exerted similar effects as enzalutamide, an AR antagonist prescribed for advanced prostate cancer, in decreasing AR signaling and expression of TMPRSS2 and angiotensin-converting enzyme 2 (ACE2), the SARS-CoV-2 receptor. However, proxalutamide led to degradation of AR protein, which was not observed with enzalutamide. Proxalutamide inhibited SARS-CoV-2 infection with an IC50 value of 97 nM, compared to 281 nM for enzalutamide. Importantly, proxalutamide inhibited infection by multiple SARS-CoV-2 variants and synergized with remdesivir. Proxalutamide protected against cell death in response to tumor necrosis factor alpha and interferon gamma, and overall survival of mice was increased with proxalutamide treatment prior to cytokine exposure. Mechanistically, we found that proxalutamide increased levels of NRF2, an essential transcription factor that mediates antioxidant responses, and decreased lung inflammation. These data provide compelling evidence that proxalutamide can prevent SARS-CoV-2 infection and cytokine-induced lung damage, suggesting that promising clinical data may emerge from ongoing phase 3 trials.

Distinct structural classes of activating FOXA1 alterations in advanced prostate cancer.

ABTRACT: Forkhead box A1 (FOXA1) is a pioneer transcription factor that is essential for the normal development of several endoderm-derived organs, including the prostate gland1,2. FOXA1 is frequently mutated in hormone-receptor-driven prostate, breast, bladder and salivary-gland tumours3-8. However, it is unclear how FOXA1 alterations affect the development of cancer, and FOXA1 has previously been ascribed both tumour-suppressive9-11 and oncogenic12-14 roles. Here we assemble an aggregate cohort of 1,546 prostate cancers and show that FOXA1 alterations fall into three structural classes that diverge in clinical incidence and genetic co-alteration profiles, with a collective prevalence of 35%. Class-1 activating mutations originate in early prostate cancer without alterations in ETS or SPOP, selectively recur within the wing-2 region of the DNA-binding forkhead domain, enable enhanced chromatin mobility and binding frequency, and strongly transactivate a luminal androgen-receptor program of prostate oncogenesis. By contrast, class-2 activating mutations are acquired in metastatic prostate cancers, truncate the C-terminal domain of FOXA1, enable dominant chromatin binding by increasing DNA affinity and-through TLE3 inactivation-promote metastasis driven by the WNT pathway. Finally, class-3 genomic rearrangements are enriched in metastatic prostate cancers, consist of duplications and translocations within the FOXA1 locus, and structurally reposition a conserved regulatory element-herein denoted FOXA1 mastermind (FOXMIND)-to drive overexpression of FOXA1 or other oncogenes. Our study reaffirms the central role of FOXA1 in mediating oncogenesis driven by the androgen receptor, and provides mechanistic insights into how the classes of FOXA1 alteration promote the initiation and/or metastatic progression of prostate cancer. These results have direct implications for understanding the pathobiology of other hormone-receptor-driven cancers and rationalize the co-targeting of FOXA1 activity in therapeutic strategies.

Expression of L1 Cell Adhesion Molecule, a Nephronal Principal Cell Marker, in Nephrogenic Adenoma.

Nephrogenic adenoma (NA) is a benign, reactive lesion seen predominantly in the urinary bladder and often associated with antecedent inflammation, instrumentation, or an operative history. Its histopathologic diversity can create diagnostic dilemmas and pathologists use morphologic evaluation along with available immunohistochemical (IHC) markers to navigate these challenges. IHC assays currently do not designate or specify NA's potential putative cell of origin. Leveraging single-cell RNA-sequencing technology, we nominated a principal (P) cell-collecting duct marker, L1 cell adhesion molecule (L1CAM), as a potential biomarker for NA. IHC characterization revealed L1CAM to be positive in all 35 (100%) patient samples of NA; negative expression was seen in the benign urothelium, benign prostatic glands, urothelial carcinoma (UCA) in situ, prostatic adenocarcinoma, majority of high-grade UCA, and metastatic UCA. In the study, we also used single-cell RNA sequencing to nominate a novel compendium of biomarkers specific for the proximal tubule, loop of Henle, and distal tubule (DT) (including P and intercalated cells), which can be used to perform nephronal mapping using RNA in situ hybridization and IHC technology. Employing this technique on NA we found enrichment of both the P-cell marker L1CAM and, the proximal tubule type-A and -B cell markers, PDZKI1P1 and PIGR, respectively. The cell-type markers for the intercalated cell of DTs (LINC01187 and FOXI1), and the loop of Henle (UMOD and IRX5), were found to be uniformly absent in NA. Overall, our findings show that based on cell type-specific implications of L1CAM expression, the shared expression pattern of L1CAM between DT P cells and NA. L1CAM expression will be of potential value in assisting surgical pathologists toward a diagnosis of NA in challenging patient samples.

Single-cell analyses of renal cell cancers reveal insights into tumor microenvironment, cell of origin, and therapy response.

Diverse subtypes of renal cell carcinomas (RCCs) display a wide spectrum of histomorphologies, proteogenomic alterations, immune cell infiltration patterns, and clinical behavior. Delineating the cells of origin for different RCC subtypes will provide mechanistic insights into their diverse pathobiology. Here, we employed single-cell RNA sequencing (scRNA-seq) to develop benign and malignant renal cell atlases. Using a random forest model trained on this cell atlas, we predicted the putative cell of origin for more than 10 RCC subtypes. scRNA-seq also revealed several attributes of the tumor microenvironment in the most common subtype of kidney cancer, clear cell RCC (ccRCC). We elucidated an active role for tumor epithelia in promoting immune cell infiltration, potentially explaining why ccRCC responds to immune checkpoint inhibitors, despite having a low neoantigen burden. In addition, we characterized an association between high endothelial cell types and lack of response to immunotherapy in ccRCC. Taken together, these single-cell analyses of benign kidney and RCC provide insight into the putative cell of origin for RCC subtypes and highlight the important role of the tumor microenvironment in influencing ccRCC biology and response to therapy.

AGO2 promotes tumor progression in KRAS-driven mouse models of non-small cell lung cancer.

Lung cancer is the deadliest malignancy in the United States. Non-small cell lung cancer (NSCLC) accounts for 85% of cases and is frequently driven by activating mutations in the gene encoding the KRAS GTPase (e.g., KRASG12D). Our previous work demonstrated that Argonaute 2 (AGO2)-a component of the RNA-induced silencing complex (RISC)-physically interacts with RAS and promotes its downstream signaling. We therefore hypothesized that AGO2 could promote KRASG12D-dependent NSCLC in vivo. To test the hypothesis, we evaluated the impact of Ago2 knockout in the KPC (LSL-KrasG12D/+;p53f/f;Cre) mouse model of NSCLC. In KPC mice, intratracheal delivery of adenoviral Cre drives lung-specific expression of a stop-floxed KRASG12D allele and biallelic ablation of p53 Simultaneous biallelic ablation of floxed Ago2 inhibited KPC lung nodule growth while reducing proliferative index and improving pathological grade. We next applied the KPHetC model, in which the Clara cell-specific CCSP-driven Cre activates KRASG12D and ablates a single p53 allele. In these mice, Ago2 ablation also reduced tumor size and grade. In both models, Ago2 knockout inhibited ERK phosphorylation (pERK) in tumor cells, indicating impaired KRAS signaling. RNA sequencing (RNA-seq) of KPC nodules and nodule-derived organoids demonstrated impaired canonical KRAS signaling with Ago2 ablation. Strikingly, accumulation of pERK in KPC organoids depended on physical interaction of AGO2 and KRAS. Taken together, our data demonstrate a pathogenic role for AGO2 in KRAS-dependent NSCLC. Given the prevalence of this malignancy and current difficulties in therapeutically targeting KRAS signaling, our work may have future translational relevance.

Targeting transcriptional regulation of SARS-CoV-2 entry factors ACE2 and TMPRSS2.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus responsible for COVID-19, employs two key host proteins to gain entry and replicate within cells, angiotensin-converting enzyme 2 (ACE2) and the cell surface transmembrane protease serine 2 (TMPRSS2). TMPRSS2 was first characterized as an androgen-regulated gene in the prostate. Supporting a role for sex hormones, males relative to females are disproportionately affected by COVID-19 in terms of mortality and morbidity. Several studies, including one employing a large epidemiological cohort, suggested that blocking androgen signaling is protective against COVID-19. Here, we demonstrate that androgens regulate the expression of ACE2, TMPRSS2, and androgen receptor (AR) in subsets of lung epithelial cells. AR levels are markedly elevated in males relative to females greater than 70 y of age. In males greater than 70 y old, smoking was associated with elevated levels of AR and ACE2 in lung epithelial cells. Transcriptional repression of the AR enhanceosome with AR or bromodomain and extraterminal domain (BET) antagonists inhibited SARS-CoV-2 infection in vitro. Taken together, these studies support further investigation of transcriptional inhibition of critical host factors in the treatment or prevention of COVID-19.

Highly Recurrent IDH1 Mutations in Prostate Cancer With Psammomatous Calcification.

Prostate cancer is a heterogeneous disease with several well-recognized morphologic subtypes and histologic variants-subsets of which are enriched for or associated with specific genomic alterations. Herein, we report a cohort of 4 unique prostate cancers characterized by intratumoral psammomatous calcification-which we have termed prostate cancer with psammomatous calcification (PCWPC). Clinicopathologic review demonstrates that PCWPCs are high-grade (grade group ≥3) tumors that involve the anterior prostate, and integrative targeted next-generation sequencing reveals recurrent hotspot IDH1 mutations. This morphology-molecular correlation is independently confirmed in The Cancer Genome Atlas prostatic adenocarcinoma cohort, with 3 of the 5 IDH1-mutant prostate cancers showing psammomatous calcification (rφ = 0.67; Fisher exact test, P < .0001). Overall, these findings suggest that PCWPC represents a novel subtype of prostate cancer enriched for an anterior location and the presence of hotspot IDH1 mutations. Recognition of these unique morphologic features could help identify IDH1-mutant prostate cancer cases retrospectively and prospectively-facilitating future large research studies and enabling clinical trial enrollment and precision medicine approaches for patients with advanced and/or aggressive disease.

Integrative clinical genomics of metastatic cancer.

Metastasis is the primary cause of cancer-related deaths. Although The Cancer Genome Atlas has sequenced primary tumour types obtained from surgical resections, much less comprehensive molecular analysis is available from clinically acquired metastatic cancers. Here we perform whole-exome and -transcriptome sequencing of 500 adult patients with metastatic solid tumours of diverse lineage and biopsy site. The most prevalent genes somatically altered in metastatic cancer included TP53, CDKN2A, PTEN, PIK3CA, and RB1. Putative pathogenic germline variants were present in 12.2% of cases of which 75% were related to defects in DNA repair. RNA sequencing complemented DNA sequencing to identify gene fusions, pathway activation, and immune profiling. Our results show that integrative sequence analysis provides a clinically relevant, multi-dimensional view of the complex molecular landscape and microenvironment of metastatic cancers.