An Atlas of Accessible Chromatin in Advanced Prostate Cancer Reveals the Epigenetic Evolution during Tumor Progression.

Metastatic castration-resistant prostate cancer (mCRPC) is a lethal disease that resists therapy targeting androgen signaling, the primary driver of prostate cancer. mCRPC resists androgen receptor (AR) inhibitors by amplifying AR signaling or by evolving into therapy-resistant subtypes that do not depend on AR. Elucidation of the epigenetic underpinnings of these subtypes could provide important insights into the drivers of therapy resistance. In this study, we produced chromatin accessibility maps linked to the binding of lineage-specific transcription factors (TF) by performing ATAC sequencing on 70 mCRPC tissue biopsies integrated with transcriptome and whole genome sequencing. mCRPC had a distinct global chromatin accessibility profile linked to AR function. Analysis of TF occupancy across accessible chromatin revealed 203 TFs associated with mCRPC subtypes. Notably, ZNF263 was identified as a putative prostate cancer TF with a significant impact on gene activity in the double-negative (AR- neuroendocrine-) subtype, potentially activating MYC targets. Overall, this analysis of chromatin accessibility in mCRPC provides valuable insights into epigenetic changes that occur during progression to mCRPC.

Integrated analyses highlight interactions between the three-dimensional genome and DNA, RNA and epigenomic alterations in metastatic prostate cancer.

The impact of variations in the three-dimensional structure of the genome has been recognized, but solid cancer tissue studies are limited. Here, we performed integrated deep Hi-C sequencing with matched whole-genome sequencing, whole-genome bisulfite sequencing, 5-hydroxymethylcytosine (5hmC) sequencing and RNA sequencing across a cohort of 80 biopsy samples from patients with metastatic castration-resistant prostate cancer. Dramatic differences were present in gene expression, 5-methylcytosine/5hmC methylation and in structural variation versus mutation rate between A and B (open and closed) chromatin compartments. A subset of tumors exhibited depleted regional chromatin contacts at the AR locus, linked to extrachromosomal circular DNA (ecDNA) and worse response to AR signaling inhibitors. We also identified topological subtypes associated with stark differences in methylation structure, gene expression and prognosis. Our data suggested that DNA interactions may predispose to structural variant formation, exemplified by the recurrent TMPRSS2-ERG fusion. This comprehensive integrated sequencing effort represents a unique clinical tumor resource.

Integrative analysis of ultra-deep RNA-seq reveals alternative promoter usage as a mechanism of activating oncogenic programmes during prostate cancer progression.

Transcription factor (TF) proteins regulate gene activity by binding to regulatory regions, most importantly at gene promoters. Many genes have alternative promoters (APs) bound by distinct TFs. The role of differential TF activity at APs during tumour development is poorly understood. Here we show, using deep RNA sequencing in 274 biopsies of benign prostate tissue, localized prostate tumours and metastatic castration-resistant prostate cancer, that AP usage increases as tumours progress and APs are responsible for a disproportionate amount of tumour transcriptional activity. Expression of the androgen receptor (AR), the key driver of prostate tumour activity, is correlated with elevated AP usage. We identified AR, FOXA1 and MYC as potential drivers of AP activation. DNA methylation is a likely mechanism for AP activation during tumour progression and lineage plasticity. Our data suggest that prostate tumours activate APs to magnify the transcriptional impact of tumour drivers, including AR and MYC.

A prostate cancer gastrointestinal transcriptional phenotype may be associated with diminished response to AR-targeted therapy.

Background: Prostate cancer is a heterogenous disease, but once it becomes metastatic it eventually becomes treatment resistant. One mechanism of resistance to AR-targeting therapy is lineage plasticity, where the tumor undergoes a transformation to an AR-indifferent phenotype, most studied in the context of neuroendocrine prostate cancer (NEPC). However, activation of additional de- or trans-differentiation programs, including a gastrointestinal (GI) gene expression program, has been suggested as an alternative method of resistance. In this study, we explored the previously identified GI prostate cancer phenotype (PCa-GI) in a large cohort of metastatic castration-resistant prostate cancer (mCRPC) patient biopsy samples. Methods: We analyzed a dataset of 634 mCRPC samples with batch effect corrected gene expression data from the West Coast Dream Team (WCDT), the East Coast Dream Team (ECDT), the Fred Hutchinson Cancer Research Center (FHCRC) and the Weill Cornell Medical center (WCM). Survival data was available from the WCDT and ECDT cohorts. We calculated a gene expression GI score using the sum of z-scores of genes from a published set of PCa-GI-defining genes (N=38). Survival analysis was performed using the Kaplan-Meier method and Cox proportional hazards regression with endpoint overall survival from time of biopsy to death of any cause. Results: We found that the PCa-GI score had a bimodal distribution, identifying a distinct set of tumors with an activated GI expression pattern. Approximately 35% of samples were classified as PCa-GI high, which was concordant with prior reports. Liver metastases had the highest median score but after excluding liver samples, 29% of the remaining samples were still classified as PCa-GI high, suggesting a distinct phenotype not exclusive to liver metastases. No correlation was observed between GI score and proliferation, AR signaling, or NEPC scores. Furthermore, the PCa-GI score was not associated with genomic alterations in AR, FOXA1, RB1, TP53 or PTEN. However, tumors with MYC amplifications showed significantly higher GI scores (p=0.0001). Patients with PCa-GI tumors had a shorter survival (HR=1.5 [1.1-2.1], p=0.02), but this result was not significant after adjusting for the liver as metastatic site (HR=1.2 [0.82-1.7], p=0.35). Patients with PCa-GI low samples had a better outcome after androgen receptor signaling inhibitors (ASI, abiraterone or enzalutamide) than other therapies (HR=0.37 [0.22-0.61], p=0.0001) while the benefit of ASI was smaller and non-significant for PCa-GI high samples (HR=0.55 [0.29-1.1], p=0.07). A differential pathway analysis identified FOXA2 signaling to be upregulated PCa-GI high tumors (FDR = 3.7 × 10-13). Conclusions: The PCa-GI phenotype is prevalent in clinical mCRPC samples and may represent a distinct biological entity. PCa-GI tumors may respond less to ASI and could offer a strategy to study novel therapeutic targets.

The Genomic and Epigenomic Landscape of Double-Negative Metastatic Prostate Cancer.

Systemic targeted therapy in prostate cancer is primarily focused on ablating androgen signaling. Androgen deprivation therapy and second-generation androgen receptor (AR)-targeted therapy selectively favor the development of treatment-resistant subtypes of metastatic castration-resistant prostate cancer (mCRPC), defined by AR and neuroendocrine (NE) markers. Molecular drivers of double-negative (AR-/NE-) mCRPC are poorly defined. In this study, we comprehensively characterized treatment-emergent mCRPC by integrating matched RNA sequencing, whole-genome sequencing, and whole-genome bisulfite sequencing from 210 tumors. AR-/NE- tumors were clinically and molecularly distinct from other mCRPC subtypes, with the shortest survival, amplification of the chromatin remodeler CHD7, and PTEN loss. Methylation changes in CHD7 candidate enhancers were linked to elevated CHD7 expression in AR-/NE+ tumors. Genome-wide methylation analysis nominated Krüppel-like factor 5 (KLF5) as a driver of the AR-/NE- phenotype, and KLF5 activity was linked to RB1 loss. These observations reveal the aggressiveness of AR-/NE- mCRPC and could facilitate the identification of therapeutic targets in this highly aggressive disease. SIGNIFICANCE: Comprehensive characterization of the five subtypes of metastatic castration-resistant prostate cancer identified transcription factors that drive each subtype and showed that the double-negative subtype has the worst prognosis.

The 5-Hydroxymethylcytosine Landscape of Prostate Cancer.

Analysis of DNA methylation is a valuable tool to understand disease progression and is increasingly being used to create diagnostic and prognostic clinical biomarkers. While conversion of cytosine to 5-methylcytosine (5mC) commonly results in transcriptional repression, further conversion to 5-hydroxymethylcytosine (5hmC) is associated with transcriptional activation. Here we perform the first study integrating whole-genome 5hmC with DNA, 5mC, and transcriptome sequencing in clinical samples of benign, localized, and advanced prostate cancer. 5hmC is shown to mark activation of cancer drivers and downstream targets. Furthermore, 5hmC sequencing revealed profoundly altered cell states throughout the disease course, characterized by increased proliferation, oncogenic signaling, dedifferentiation, and lineage plasticity to neuroendocrine and gastrointestinal lineages. Finally, 5hmC sequencing of cell-free DNA from patients with metastatic disease proved useful as a prognostic biomarker able to identify an aggressive subtype of prostate cancer using the genes TOP2A and EZH2, previously only detectable by transcriptomic analysis of solid tumor biopsies. Overall, these findings reveal that 5hmC marks epigenomic activation in prostate cancer and identify hallmarks of prostate cancer progression with potential as biomarkers of aggressive disease. SIGNIFICANCE: In prostate cancer, 5-hydroxymethylcytosine delineates oncogene activation and stage-specific cell states and can be analyzed in liquid biopsies to detect cancer phenotypes. See related article by Wu and Attard, p. 3880.

Predicting cancer drug TARGETS - TreAtment Response Generalized Elastic-neT Signatures.

We are now in an era of molecular medicine, where specific DNA alterations can be used to identify patients who will respond to specific drugs. However, there are only a handful of clinically used predictive biomarkers in oncology. Herein, we describe an approach utilizing in vitro DNA and RNA sequencing and drug response data to create TreAtment Response Generalized Elastic-neT Signatures (TARGETS). We trained TARGETS drug response models using Elastic-Net regression in the publicly available Genomics of Drug Sensitivity in Cancer (GDSC) database. Models were then validated on additional in-vitro data from the Cancer Cell Line Encyclopedia (CCLE), and on clinical samples from The Cancer Genome Atlas (TCGA) and Stand Up to Cancer/Prostate Cancer Foundation West Coast Prostate Cancer Dream Team (WCDT). First, we demonstrated that all TARGETS models successfully predicted treatment response in the separate in-vitro CCLE treatment response dataset. Next, we evaluated all FDA-approved biomarker-based cancer drug indications in TCGA and demonstrated that TARGETS predictions were concordant with established clinical indications. Finally, we performed independent clinical validation in the WCDT and found that the TARGETS AR signaling inhibitors (ARSI) signature successfully predicted clinical treatment response in metastatic castration-resistant prostate cancer with a statistically significant interaction between the TARGETS score and PSA response (p = 0.0252). TARGETS represents a pan-cancer, platform-independent approach to predict response to oncologic therapies and could be used as a tool to better select patients for existing therapies as well as identify new indications for testing in prospective clinical trials.

An integrated functional and clinical genomics approach reveals genes driving aggressive metastatic prostate cancer.

Genomic sequencing of thousands of tumors has revealed many genes associated with specific types of cancer. Similarly, large scale CRISPR functional genomics efforts have mapped genes required for cancer cell proliferation or survival in hundreds of cell lines. Despite this, for specific disease subtypes, such as metastatic prostate cancer, there are likely a number of undiscovered tumor specific driver genes that may represent potential drug targets. To identify such genetic dependencies, we performed genome-scale CRISPRi screens in metastatic prostate cancer models. We then created a pipeline in which we integrated pan-cancer functional genomics data with our metastatic prostate cancer functional and clinical genomics data to identify genes that can drive aggressive prostate cancer phenotypes. Our integrative analysis of these data reveals known prostate cancer specific driver genes, such as AR and HOXB13, as well as a number of top hits that are poorly characterized. In this study we highlight the strength of an integrated clinical and functional genomics pipeline and focus on two top hit genes, KIF4A and WDR62. We demonstrate that both KIF4A and WDR62 drive aggressive prostate cancer phenotypes in vitro and in vivo in multiple models, irrespective of AR-status, and are also associated with poor patient outcome.

Multiomics Characterization of Low-Grade Serous Ovarian Carcinoma Identifies Potential Biomarkers of MEK Inhibitor Sensitivity and Therapeutic Vulnerability.

Low-grade serous ovarian carcinoma (LGSOC) is a rare tumor subtype with high case fatality rates in patients with metastatic disease. There is a pressing need to develop effective treatments using newly available preclinical models for therapeutic discovery and drug evaluation. Here, we use multiomics integration of whole-exome sequencing, RNA sequencing, and mass spectrometry-based proteomics on 14 LGSOC cell lines to elucidate novel biomarkers and therapeutic vulnerabilities. Comparison of LGSOC cell line data with LGSOC tumor data enabled predictive biomarker identification of MEK inhibitor (MEKi) efficacy, with KRAS mutations found exclusively in MEKi-sensitive cell lines and NRAS mutations found mostly in MEKi-resistant cell lines. Distinct patterns of Catalogue of Somatic Mutations in Cancer mutational signatures were identified in MEKi-sensitive and MEKi-resistant cell lines. Deletions of CDKN2A/B and MTAP genes were more frequent in cell lines than tumor samples and possibly represent key driver events in the absence of KRAS/NRAS/BRAF mutations. These LGSOC cell lines were representative models of the molecular aberrations found in LGSOC tumors. For prediction of in vitro MEKi efficacy, proteomic data provided better discrimination than gene expression data. Condensin, minichromosome maintenance, and replication factor C protein complexes were identified as potential treatment targets in MEKi-resistant cell lines. This study suggests that CDKN2A/B or MTAP deficiency may be exploited using synthetically lethal treatment strategies, highlighting the importance of using proteomic data as a tool for molecular drug prediction. Multiomics approaches are crucial to improving our understanding of the molecular underpinnings of LGSOC and applying this information to develop new therapies. SIGNIFICANCE: These findings highlight the utility of global multiomics to characterize LGSOC cell lines as research models, to determine biomarkers of MEKi resistance, and to identify potential novel therapeutic targets.

Autoantibody Landscape in Patients with Advanced Prostate Cancer.

PURPOSE: Autoantibody responses in cancer are of great interest, as they may be concordant with T-cell responses to cancer antigens or predictive of response to cancer immunotherapies. Thus, we sought to characterize the antibody landscape of metastatic castration-resistant prostate cancer (mCRPC). EXPERIMENTAL DESIGN: Serum antibody epitope repertoire analysis (SERA) was performed on patient serum to identify tumor-specific neoepitopes. Somatic mutation-specific neoepitopes were investigated by associating serum epitope enrichment scores with whole-genome sequencing results from paired solid tumor metastasis biopsies and germline blood samples. A protein-based immunome-wide association study (PIWAS) was performed to identify significantly enriched epitopes, and candidate serum antibodies enriched in select patients were validated by ELISA profiling. A distinct cohort of patients with melanoma was evaluated to validate the top cancer-specific epitopes. RESULTS: SERA was performed on 1,229 serum samples obtained from 72 men with mCRPC and 1,157 healthy control patients. Twenty-nine of 6,636 somatic mutations (0.44%) were associated with an antibody response specific to the mutated peptide. PIWAS analyses identified motifs in 11 proteins, including NY-ESO-1 and HERVK-113, as immunogenic in mCRPC, and ELISA confirmed serum antibody enrichment in candidate patients. Confirmatory PIWAS, Identifying Motifs Using Next-generation sequencing Experiments (IMUNE), and ELISA analyses performed on serum samples from 106 patients with melanoma similarly revealed enriched cancer-specific antibody responses to NY-ESO-1. CONCLUSIONS: We present the first large-scale profiling of autoantibodies in advanced prostate cancer, utilizing a new antibody profiling approach to reveal novel cancer-specific antigens and epitopes. Our study recovers antigens of known importance and identifies novel tumor-specific epitopes of translational interest.

Identification of conserved evolutionary trajectories in tumors.

MOTIVATION: As multi-region, time-series and single-cell sequencing data become more widely available; it is becoming clear that certain tumors share evolutionary characteristics with others. In the last few years, several computational methods have been developed with the goal of inferring the subclonal composition and evolutionary history of tumors from tumor biopsy sequencing data. However, the phylogenetic trees that they report differ significantly between tumors (even those with similar characteristics). RESULTS: In this article, we present a novel combinatorial optimization method, CONETT, for detection of recurrent tumor evolution trajectories. Our method constructs a consensus tree of conserved evolutionary trajectories based on the information about temporal order of alteration events in a set of tumors. We apply our method to previously published datasets of 100 clear-cell renal cell carcinoma and 99 non-small-cell lung cancer patients and identify both conserved trajectories that were reported in the original studies, as well as new trajectories. AVAILABILITY AND IMPLEMENTATION: CONETT is implemented in C++ and available at https://github.com/ehodzic/CONETT. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Well-Differentiated Papillary Mesothelioma of the Peritoneum Is Genetically Distinct from Malignant Mesothelioma.

Well-differentiated papillary mesothelioma (WDPM) is an uncommon mesothelial proliferation that is most commonly encountered as an incidental finding in the peritoneal cavity. There is controversy in the literature about whether WDPM is a neoplasm or a reactive process and, if neoplastic, whether it is a variant or precursor of epithelial malignant mesothelioma or is a different entity. Using whole exome sequencing of five WDPMs of the peritoneum, we have identified distinct mutations in EHD1, ATM, FBXO10, SH2D2A, CDH5, MAGED1, and TP73 shared by WDPM cases but not reported in malignant mesotheliomas. Furthermore, we show that WDPM is strongly enriched with C > A transversion substitution mutations, a pattern that is also not found in malignant mesotheliomas. The WDPMs lacked the alterations involving BAP1, SETD2, NF2, CDKN2A/B, LASTS1/2, PBRM1, and SMARCC1 that are frequently found in malignant mesotheliomas. We conclude that WDPMs are neoplasms that are genetically distinct from malignant mesotheliomas and, based on observed mutations, do not appear to be precursors of malignant mesotheliomas.

Y-box binding protein-1 is crucial in acquired drug resistance development in metastatic clear-cell renal cell carcinoma.

BACKGROUND: Renal cell carcinoma (RCC) is a highly vascular tumor and patients with low risk metastatic RCC of clear-cell histological sub-type (mccRCC) are treated with tyrosine-kinase inhibitors (TKIs), sunitinib, as the first-line of treatment. Unfortunately, TKI resistance eventually develops, and the underlying molecular mechanism is not well understood. METHODS: RCC cell-line with metastatic clear-cell histology (Caki-1), and patient samples were analysed to identify the role of Y-box binding protein 1 (YB-1) and ATP-binding cassette sub-family B member 1 (ABCB-1) in acquired sunitinib-resistance development. Caki-1 was conditioned with increasing sunitinib doses to recapitulate acquired resistance development in clinics. Sunitinib-conditioned and wild-type Caki-1 were subjected to cell viability assay, scratch assay, chicken embryo chorioallantoic membrane engraftment and proteomics analysis. Classical biochemical assays like flow cytometry, immunofluorescent staining, immunohistochemical staining, optical coherence tomography imaging, Western Blot and RT-PCR assays were applied to determine the possible mechanism of sunitinib-resistance development and the effect of drug treatments. Publicly available data was also used to determine the role of YB-1 upregulation in ccRCC and the patients' overall survival. RESULTS: We demonstrate that YB-1 and ABCB-1 are upregulated in sunitinib-resistant in vitro, ex vivo, in vivo and patient samples compared to the sensitive samples. This provides evidence to a mechanism of acquired sunitinib-resistance development in mccRCC. Furthermore, our results establish that inhibiting ABCB-1 with elacridar, in addition to sunitinib, has a positive impact on reverting sunitinib-resistance development in in vitro, ex vivo and in vivo models. CONCLUSION: This work proposes a targeted therapy (elacridar and sunitinib) to re-sensitize sunitinib-resistant mccRCC and, possibly, slow disease progression.

Pathway and network analysis of more than 2500 whole cancer genomes.

The catalog of cancer driver mutations in protein-coding genes has greatly expanded in the past decade. However, non-coding cancer driver mutations are less well-characterized and only a handful of recurrent non-coding mutations, most notably TERT promoter mutations, have been reported. Here, as part of the ICGC/TCGA Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium, which aggregated whole genome sequencing data from 2658 cancer across 38 tumor types, we perform multi-faceted pathway and network analyses of non-coding mutations across 2583 whole cancer genomes from 27 tumor types compiled by the ICGC/TCGA PCAWG project that was motivated by the success of pathway and network analyses in prioritizing rare mutations in protein-coding genes. While few non-coding genomic elements are recurrently mutated in this cohort, we identify 93 genes harboring non-coding mutations that cluster into several modules of interacting proteins. Among these are promoter mutations associated with reduced mRNA expression in TP53, TLE4, and TCF4. We find that biological processes had variable proportions of coding and non-coding mutations, with chromatin remodeling and proliferation pathways altered primarily by coding mutations, while developmental pathways, including Wnt and Notch, altered by both coding and non-coding mutations. RNA splicing is primarily altered by non-coding mutations in this cohort, and samples containing non-coding mutations in well-known RNA splicing factors exhibit similar gene expression signatures as samples with coding mutations in these genes. These analyses contribute a new repertoire of possible cancer genes and mechanisms that are altered by non-coding mutations and offer insights into additional cancer vulnerabilities that can be investigated for potential therapeutic treatments.

Identification of gene signature for treatment response to guide precision oncology in clear-cell renal cell carcinoma.

Clear-cell renal cell carcinoma (ccRCC) is a common therapy resistant disease with aberrant angiogenic and immunosuppressive features. Patients with metastatic disease are treated with targeted therapies based on clinical features: low-risk patients are usually treated with anti-angiogenic drugs and intermediate/high-risk patients with immune therapy. However, there are no biomarkers available to guide treatment choice for these patients. A recently published phase II clinical trial observed a correlation between ccRCC patients' clustering and their response to targeted therapy. However, the clustering of these groups was not distinct. Here, we analyzed the gene expression profile of 469 ccRCC patients, using featured selection technique, and have developed a refined 66-gene signature for improved sub-classification of patients. Moreover, we have identified a novel comprehensive expression profile to distinguish between migratory stromal and immune cells. Furthermore, the proposed 66-gene signature was validated using a different cohort of 64 ccRCC patients. These findings are foundational for the development of reliable biomarkers that may guide treatment decision-making and improve therapy response in ccRCC patients.

Combinatorial Detection of Conserved Alteration Patterns for Identifying Cancer Subnetworks.

BACKGROUND: Advances in large-scale tumor sequencing have led to an understanding that there are combinations of genomic and transcriptomic alterations specific to tumor types, shared across many patients. Unfortunately, computational identification of functionally meaningful and recurrent alteration patterns within gene/protein interaction networks has proven to be challenging. FINDINGS: We introduce a novel combinatorial method, cd-CAP (combinatorial detection of conserved alteration patterns), for simultaneous detection of connected subnetworks of an interaction network where genes exhibit conserved alteration patterns across tumor samples. Our method differentiates distinct alteration types associated with each gene (rather than relying on binary information of a gene being altered or not) and simultaneously detects multiple alteration profile conserved subnetworks. CONCLUSIONS: In a number of The Cancer Genome Atlas datasets, cd-CAP identified large biologically significant subnetworks with conserved alteration patterns, shared across many tumor samples.

BAP1 haploinsufficiency predicts a distinct immunogenic class of malignant peritoneal mesothelioma.

BACKGROUND: Malignant peritoneal mesothelioma (PeM) is a rare and fatal cancer that originates from the peritoneal lining of the abdomen. Standard treatment of PeM is limited to cytoreductive surgery and/or chemotherapy, and no effective targeted therapies for PeM exist. Some immune checkpoint inhibitor studies of mesothelioma have found positivity to be associated with a worse prognosis. METHODS: To search for novel therapeutic targets for PeM, we performed a comprehensive integrative multi-omics analysis of the genome, transcriptome, and proteome of 19 treatment-naïve PeM, and in particular, we examined BAP1 mutation and copy number status and its relationship to immune checkpoint inhibitor activation. RESULTS: We found that PeM could be divided into tumors with an inflammatory tumor microenvironment and those without and that this distinction correlated with haploinsufficiency of BAP1. To further investigate the role of BAP1, we used our recently developed cancer driver gene prioritization algorithm, HIT'nDRIVE, and observed that PeM with BAP1 haploinsufficiency form a distinct molecular subtype characterized by distinct gene expression patterns of chromatin remodeling, DNA repair pathways, and immune checkpoint receptor activation. We demonstrate that this subtype is correlated with an inflammatory tumor microenvironment and thus is a candidate for immune checkpoint blockade therapies. CONCLUSIONS: Our findings reveal BAP1 to be a potential, easily trackable prognostic and predictive biomarker for PeM immunotherapy that refines PeM disease classification. BAP1 stratification may improve drug response rates in ongoing phases I and II clinical trials exploring the use of immune checkpoint blockade therapies in PeM in which BAP1 status is not considered. This integrated molecular characterization provides a comprehensive foundation for improved management of a subset of PeM patients.

High-throughput detection of RNA processing in bacteria.

BACKGROUND: Understanding the RNA processing of an organism's transcriptome is an essential but challenging step in understanding its biology. Here we investigate with unprecedented detail the transcriptome of Pseudomonas aeruginosa PAO1, a medically important and innately multi-drug resistant bacterium. We systematically mapped RNA cleavage and dephosphorylation sites that result in 5'-monophosphate terminated RNA (pRNA) using monophosphate RNA-Seq (pRNA-Seq). Transcriptional start sites (TSS) were also mapped using differential RNA-Seq (dRNA-Seq) and both datasets were compared to conventional RNA-Seq performed in a variety of growth conditions. RESULTS: The pRNA-Seq library revealed known tRNA, rRNA and transfer-messenger RNA (tmRNA) processing sites, together with previously uncharacterized RNA cleavage events that were found disproportionately near the 5' ends of transcripts associated with basic bacterial functions such as oxidative phosphorylation and purine metabolism. The majority (97%) of the processed mRNAs were cleaved at precise codon positions within defined sequence motifs indicative of distinct endonucleolytic activities. The most abundant of these motifs corresponded closely to an E. coli RNase E site previously established in vitro. Using the dRNA-Seq library, we performed an operon analysis and predicted 3159 potential TSS. A correlation analysis uncovered 105 antiparallel pairs of TSS that were separated by 18 bp from each other and were centered on single palindromic TAT(A/T)ATA motifs (likely - 10 promoter elements), suggesting that, consistent with previous in vitro experimentation, these sites can initiate transcription bi-directionally and may thus provide a novel form of transcriptional regulation. TSS and RNA-Seq analysis allowed us to confirm expression of small non-coding RNAs (ncRNAs), many of which are differentially expressed in swarming and biofilm formation conditions. CONCLUSIONS: This study uses pRNA-Seq, a method that provides a genome-wide survey of RNA processing, to study the bacterium Pseudomonas aeruginosa and discover extensive transcript processing not previously appreciated. We have also gained novel insight into RNA maturation and turnover as well as a potential novel form of transcription regulation. NOTE: All sequence data has been submitted to the NCBI sequence read archive. Accession numbers are as follows: [NCBI sequence read archive: SRX156386, SRX157659, SRX157660, SRX157661, SRX157683 and SRX158075]. The sequence data is viewable using Jbrowse on www.pseudomonas.com .

HIT'nDRIVE: patient-specific multidriver gene prioritization for precision oncology.

Prioritizing molecular alterations that act as drivers of cancer remains a crucial bottleneck in therapeutic development. Here we introduce HIT'nDRIVE, a computational method that integrates genomic and transcriptomic data to identify a set of patient-specific, sequence-altered genes, with sufficient collective influence over dysregulated transcripts. HIT'nDRIVE aims to solve the "random walk facility location" (RWFL) problem in a gene (or protein) interaction network, which differs from the standard facility location problem by its use of an alternative distance measure: "multihitting time," the expected length of the shortest random walk from any one of the set of sequence-altered genes to an expression-altered target gene. When applied to 2200 tumors from four major cancer types, HIT'nDRIVE revealed many potentially clinically actionable driver genes. We also demonstrated that it is possible to perform accurate phenotype prediction for tumor samples by only using HIT'nDRIVE-seeded driver gene modules from gene interaction networks. In addition, we identified a number of breast cancer subtype-specific driver modules that are associated with patients' survival outcome. Furthermore, HIT'nDRIVE, when applied to a large panel of pan-cancer cell lines, accurately predicted drug efficacy using the driver genes and their seeded gene modules. Overall, HIT'nDRIVE may help clinicians contextualize massive multiomics data in therapeutic decision making, enabling widespread implementation of precision oncology.