Widespread and Functional RNA Circularization in Localized Prostate Cancer.

The cancer transcriptome is remarkably complex, including low-abundance transcripts, many not polyadenylated. To fully characterize the transcriptome of localized prostate cancer, we performed ultra-deep total RNA-seq on 144 tumors with rich clinical annotation. This revealed a linear transcriptomic subtype associated with the aggressive intraductal carcinoma sub-histology and a fusion profile that differentiates localized from metastatic disease. Analysis of back-splicing events showed widespread RNA circularization, with the average tumor expressing 7,232 circular RNAs (circRNAs). The degree of circRNA production was correlated to disease progression in multiple patient cohorts. Loss-of-function screening identified 11.3% of highly abundant circRNAs as essential for cell proliferation; for ∼90% of these, their parental linear transcripts were not essential. Individual circRNAs can have distinct functions, with circCSNK1G3 promoting cell growth by interacting with miR-181. These data advocate for adoption of ultra-deep RNA-seq without poly-A selection to interrogate both linear and circular transcriptomes.

Phenotype prediction from single-cell RNA-seq data using attention-based neural networks.

MOTIVATION: A patient's disease phenotype can be driven and determined by specific groups of cells whose marker genes are either unknown or can only be detected at late-stage using conventional bulk assays such as RNA-Seq technology. Recent advances in single-cell RNA sequencing (scRNA-seq) enable gene expression profiling in cell-level resolution, and therefore have the potential to identify those cells driving the disease phenotype even while the number of these cells is small. However, most existing methods rely heavily on accurate cell type detection, and the number of available annotated samples is usually too small for training deep learning predictive models. RESULTS: Here, we propose the method ScRAT for phenotype prediction using scRNA-seq data. To train ScRAT with a limited number of samples of different phenotypes, such as coronavirus disease (COVID) and non-COVID, ScRAT first applies a mixup module to increase the number of training samples. A multi-head attention mechanism is employed to learn the most informative cells for each phenotype without relying on a given cell type annotation. Using three public COVID datasets, we show that ScRAT outperforms other phenotype prediction methods. The performance edge of ScRAT over its competitors increases as the number of training samples decreases, indicating the efficacy of our sample mixup. Critical cell types detected based on high-attention cells also support novel findings in the original papers and the recent literature. This suggests that ScRAT overcomes the challenge of missing marker genes and limited sample number with great potential revealing novel molecular mechanisms and/or therapies. AVAILABILITY AND IMPLEMENTATION: The code of our proposed method ScRAT is published at https://github.com/yuzhenmao/ScRAT.

Establishment and validation of preclinical models of SMARCA4-inactivated and ARID1A/ARID1B co-inactivated dedifferentiated endometrial carcinoma.

OBJECTIVE: Dedifferentiated endometrial cancer (DDEC) is an uncommon and clinically highly aggressive subtype of endometrial cancer characterized by genomic inactivation of SWItch/Sucrose Non-Fermentable (SWI/SNF) complex protein. It responds poorly to conventional systemic treatment and its rapidly progressive clinical course limits the therapeutic windows to trial additional lines of therapies. This underscores a pressing need for biologically accurate preclinical tumor models to accelerate therapeutic development. METHODS: DDEC tumor from surgical samples were implanted into immunocompromised mice for patient-derived xenograft (PDX) and cell line development. The histologic, immunophenotypic, genetic and epigenetic features of the patient tumors and the established PDX models were characterized. The SMARCA4-deficienct DDEC model was evaluated for its sensitivity toward a KDM6A/B inhibitor (GSK-J4) that was previously reported to be effective therapy for other SMARCA4-deficient cancer types. RESULTS: All three DDEC models exhibited rapid growth in vitro and in vivo, with two PDX models showing spontaneous development of metastases in vivo. The PDX tumors maintained the same undifferentiated histology and immunophenotype, and exhibited identical genomic and methylation profiles as seen in the respective parental tumors, including a mismatch repair (MMR)-deficient DDEC with genomic inactivation of SMARCA4, and two MMR-deficient DDECs with genomic inactivation of both ARID1A and ARID1B. Although the SMARCA4-deficient cell line showed low micromolecular sensitivity to GSK-J4, no significant tumor growth inhibition was observed in the corresponding PDX model. CONCLUSIONS: These established patient tumor-derived models accurately depict DDEC and represent valuable preclinical tools to gain therapeutic insights into this aggressive tumor type.

GRB10 sustains AR activity by interacting with PP2A in prostate cancer cells.

Prostate cancer (PCa) progression is driven by androgen receptor (AR) signaling. Unfortunately, androgen-deprivation therapy and the use of even more potent AR pathway inhibitors (ARPIs) cannot bring about a cure. ARPI resistance (ie, castration-resistant PCa, CRPC) will inevitably develop. Previously, we demonstrated that GRB10 is an AR transcriptionally repressed gene that functionally contributes to CRPC development and ARPI resistance. GRB10 expression is elevated prior to CRPC development in our patient-derived xenograft models and is significantly upregulated in clinical CRPC samples. Here, we analyzed transcriptomic data from GRB10 knockdown in PCa cells and found that AR signaling is downregulated. While the mRNA expression of AR target genes decreased upon GRB10 knockdown, AR expression was not affected at the mRNA or protein level. We further found that phosphorylation of AR serine 81 (S81), which is critical for AR transcriptional activity, is decreased by GRB10 knockdown and increased by its overexpression. Luciferase assay using GRB10-knockdown cells also indicate reduced AR activity. Immunoprecipitation coupled with mass spectrometry revealed an interaction between GRB10 and the PP2A complex, which is a known phosphatase of AR. Further validations and analyses showed that GRB10 binds to the PP2Ac catalytic subunit with its PH domain. Mechanistically, GRB10 knockdown increased PP2Ac protein stability, which in turn decreased AR S81 phosphorylation and reduced AR activity. Our findings indicate a reciprocal feedback between GRB10 and AR signaling, implying the importance of GRB10 in PCa progression.

Establishment and characterization of a novel treatment-related neuroendocrine prostate cancer cell line KUCaP13.

The prevalence of neuroendocrine prostate cancer (NEPC) arising from adenocarcinoma (AC) upon potent androgen receptor (AR) pathway inhibition is increasing. Deeper understanding of NEPC biology and development of novel therapeutic agents are needed. However, research is hindered by the paucity of research models, especially cell lines developed from NEPC patients. We established a novel NEPC cell line, KUCaP13, from tissue of a patient initially diagnosed with AC which later recurred as NEPC. The cell line has been maintained permanently in vitro under regular cell culture conditions and is amenable to gene engineering with lentivirus. KUCaP13 cells lack the expression of AR and overexpress NEPC-associated genes, including SOX2, EZH2, AURKA, PEG10, POU3F2, ENO2, and FOXA2. Importantly, the cell line maintains the homozygous deletion of CHD1, which was confirmed in the primary AC of the index patient. Loss of heterozygosity of TP53 and PTEN, and an allelic loss of RB1 with a transcriptomic signature compatible with Rb pathway aberration were revealed. Knockdown of PEG10 using shRNA significantly suppressed growth in vivo. Introduction of luciferase allowed serial monitoring of cells implanted orthotopically or in the renal subcapsule. Although H3K27me was reduced by EZH2 inhibition, reversion to AC was not observed. KUCaP13 is the first patient-derived, treatment-related NEPC cell line with triple loss of tumor suppressors critical for NEPC development through lineage plasticity. It could be valuable in research to deepen the understanding of NEPC.

AITL: Adversarial Inductive Transfer Learning with input and output space adaptation for pharmacogenomics.

MOTIVATION: The goal of pharmacogenomics is to predict drug response in patients using their single- or multi-omics data. A major challenge is that clinical data (i.e. patients) with drug response outcome is very limited, creating a need for transfer learning to bridge the gap between large pre-clinical pharmacogenomics datasets (e.g. cancer cell lines), as a source domain, and clinical datasets as a target domain. Two major discrepancies exist between pre-clinical and clinical datasets: (i) in the input space, the gene expression data due to difference in the basic biology, and (ii) in the output space, the different measures of the drug response. Therefore, training a computational model on cell lines and testing it on patients violates the i.i.d assumption that train and test data are from the same distribution. RESULTS: We propose Adversarial Inductive Transfer Learning (AITL), a deep neural network method for addressing discrepancies in input and output space between the pre-clinical and clinical datasets. AITL takes gene expression of patients and cell lines as the input, employs adversarial domain adaptation and multi-task learning to address these discrepancies, and predicts the drug response as the output. To the best of our knowledge, AITL is the first adversarial inductive transfer learning method to address both input and output discrepancies. Experimental results indicate that AITL outperforms state-of-the-art pharmacogenomics and transfer learning baselines and may guide precision oncology more accurately. AVAILABILITY AND IMPLEMENTATION: https://github.com/hosseinshn/AITL. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

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.

CircMiner: accurate and rapid detection of circular RNA through splice-aware pseudo-alignment scheme.

MOTIVATION: The ubiquitous abundance of circular RNAs (circRNAs) has been revealed by performing high-throughput sequencing in a variety of eukaryotes. circRNAs are related to some diseases, such as cancer in which they act as oncogenes or tumor-suppressors and, therefore, have the potential to be used as biomarkers or therapeutic targets. Accurate and rapid detection of circRNAs from short reads remains computationally challenging. This is due to the fact that identifying chimeric reads, which is essential for finding back-splice junctions, is a complex process. The sensitivity of discovery methods, to a high degree, relies on the underlying mapper that is used for finding chimeric reads. Furthermore, all the available circRNA discovery pipelines are resource intensive. RESULTS: We introduce CircMiner, a novel stand-alone circRNA detection method that rapidly identifies and filters out linear RNA sequencing reads and detects back-splice junctions. CircMiner employs a rapid pseudo-alignment technique to identify linear reads that originate from transcripts, genes or the genome. CircMiner further processes the remaining reads to identify the back-splice junctions and detect circRNAs with single-nucleotide resolution. We evaluated the efficacy of CircMiner using simulated datasets generated from known back-splice junctions and showed that CircMiner has superior accuracy and speed compared to the existing circRNA detection tools. Additionally, on two RNase R treated cell line datasets, CircMiner was able to detect most of consistent, high confidence circRNAs compared to untreated samples of the same cell line. AVAILABILITY AND IMPLEMENTATION: CircMiner is implemented in C++ and is available online at https://github.com/vpc-ccg/circminer. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Structural variation and fusion detection using targeted sequencing data from circulating cell free DNA.

MOTIVATION: Cancer is a complex disease that involves rapidly evolving cells, often forming multiple distinct clones. In order to effectively understand progression of a patient-specific tumor, one needs to comprehensively sample tumor DNA at multiple time points, ideally obtained through inexpensive and minimally invasive techniques. Current sequencing technologies make the 'liquid biopsy' possible, which involves sampling a patient's blood or urine and sequencing the circulating cell free DNA (cfDNA). A certain percentage of this DNA originates from the tumor, known as circulating tumor DNA (ctDNA). The ratio of ctDNA may be extremely low in the sample, and the ctDNA may originate from multiple tumors or clones. These factors present unique challenges for applying existing tools and workflows to the analysis of ctDNA, especially in the detection of structural variations which rely on sufficient read coverage to be detectable. RESULTS: Here we introduce SViCT , a structural variation (SV) detection tool designed to handle the challenges associated with cfDNA analysis. SViCT can detect breakpoints and sequences of various structural variations including deletions, insertions, inversions, duplications and translocations. SViCT extracts discordant read pairs, one-end anchors and soft-clipped/split reads, assembles them into contigs, and re-maps contig intervals to a reference genome using an efficient k-mer indexing approach. The intervals are then joined using a combination of graph and greedy algorithms to identify specific structural variant signatures. We assessed the performance of SViCT and compared it to state-of-the-art tools using simulated cfDNA datasets with properties matching those of real cfDNA samples. The positive predictive value and sensitivity of our tool was superior to all the tested tools and reasonable performance was maintained down to the lowest dilution of 0.01% tumor DNA in simulated datasets. Additionally, SViCT was able to detect all known SVs in two real cfDNA reference datasets (at 0.6-5% ctDNA) and predict a novel structural variant in a prostate cancer cohort. AVAILABILITY: SViCT is available at https://github.com/vpc-ccg/svict. Contact:faraz.hach@ubc.ca.

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.

MOLI: multi-omics late integration with deep neural networks for drug response prediction.

MOTIVATION: Historically, gene expression has been shown to be the most informative data for drug response prediction. Recent evidence suggests that integrating additional omics can improve the prediction accuracy which raises the question of how to integrate the additional omics. Regardless of the integration strategy, clinical utility and translatability are crucial. Thus, we reasoned a multi-omics approach combined with clinical datasets would improve drug response prediction and clinical relevance. RESULTS: We propose MOLI, a multi-omics late integration method based on deep neural networks. MOLI takes somatic mutation, copy number aberration and gene expression data as input, and integrates them for drug response prediction. MOLI uses type-specific encoding sub-networks to learn features for each omics type, concatenates them into one representation and optimizes this representation via a combined cost function consisting of a triplet loss and a binary cross-entropy loss. The former makes the representations of responder samples more similar to each other and different from the non-responders, and the latter makes this representation predictive of the response values. We validate MOLI on in vitro and in vivo datasets for five chemotherapy agents and two targeted therapeutics. Compared to state-of-the-art single-omics and early integration multi-omics methods, MOLI achieves higher prediction accuracy in external validations. Moreover, a significant improvement in MOLI's performance is observed for targeted drugs when training on a pan-drug input, i.e. using all the drugs with the same target compared to training only on drug-specific inputs. MOLI's high predictive power suggests it may have utility in precision oncology. AVAILABILITY AND IMPLEMENTATION: https://github.com/hosseinshn/MOLI. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Alignment-free clustering of UMI tagged DNA molecules.

MOTIVATION: Next-Generation Sequencing has led to the availability of massive genomic datasets whose processing raises many challenges, including the handling of sequencing errors. This is especially pertinent in cancer genomics, e.g. for detecting low allele frequency variations from circulating tumor DNA. Barcode tagging of DNA molecules with unique molecular identifiers (UMI) attempts to mitigate sequencing errors; UMI tagged molecules are polymerase chain reaction (PCR) amplified, and the PCR copies of UMI tagged molecules are sequenced independently. However, the PCR and sequencing steps can generate errors in the sequenced reads that can be located in the barcode and/or the DNA sequence. Analyzing UMI tagged sequencing data requires an initial clustering step, with the aim of grouping reads sequenced from PCR duplicates of the same UMI tagged molecule into a single cluster, and the size of the current datasets requires this clustering process to be resource-efficient. RESULTS: We introduce Calib, a computational tool that clusters paired-end reads from UMI tagged sequencing experiments generated by substitution-error-dominant sequencing platforms such as Illumina. Calib clusters are defined as connected components of a graph whose edges are defined in terms of both barcode similarity and read sequence similarity. The graph is constructed efficiently using locality sensitive hashing and MinHashing techniques. Calib's default clustering parameters are optimized empirically, for different UMI and read lengths, using a simulation module that is packaged with Calib. Compared to other tools, Calib has the best accuracy on simulated data, while maintaining reasonable runtime and memory footprint. On a real dataset, Calib runs with far less resources than alignment-based methods, and its clusters reduce the number of tentative false positive in downstream variation calling. AVAILABILITY AND IMPLEMENTATION: Calib is implemented in C++ and its simulation module is implemented in Python. Calib is available at https://github.com/vpc-ccg/calib. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Modeling human prostate cancer progression in vitro.

Detailed mechanisms involved in prostate cancer (CaP) development and progression are not well understood. Current experimental models used to study CaP are not well suited to address this issue. Previously, we have described the hormonal progression of non-tumorigenic human prostate epithelial cells (BPH1) into malignant cells via tissue recombination. Here, we describe a method to derive human cell lines from distinct stages of CaP that parallel cellular, genetic and epigenetic changes found in patients with cancers. This BPH1-derived Cancer Progression (BCaP) model represents different stages of cancer. Using diverse analytical strategies, we show that the BCaP model reproduces molecular characteristics of CaP in human patients. Furthermore, we demonstrate that BCaP cells have altered gene expression of shared pathways with human and transgenic mouse CaP data, as well as, increasing genomic instability with TMPRSS2-ERG fusion in advanced tumor cells. Together, these cell lines represent a unique model of human CaP progression providing a novel tool that will allow the discovery and experimental validation of mechanisms regulating human CaP development and progression. This BPH1-derived Cancer Progression (BCaP) model represents different stages of cancer. The BCaP model reproduces molecular characteristics of prostate cancer. The cells have altered gene expression with TMPRSS2-ERG fusion representing a unique model for prostate cancer progression.

Metagenomic and metatranscriptomic analysis of human prostate microbiota from patients with prostate cancer.

BACKGROUND: Prostate cancer (PCa) is the most common malignant neoplasm among men in many countries. Since most precancerous and cancerous tissues show signs of inflammation, chronic bacterial prostatitis has been hypothesized to be a possible etiology. However, establishing a causal relationship between microbial inflammation and PCa requires a comprehensive analysis of the prostate microbiome. The aim of this study was to characterize the microbiome in prostate tissue of PCa patients and investigate its association with tumour clinical characteristics as well as host expression profiles. RESULTS: The metagenome and metatranscriptome of tumour and the adjacent benign tissues were assessed in 65 Chinese radical prostatectomy specimens. Escherichia, Propionibacterium, Acinetobacter and Pseudomonas were abundant in both metagenome and metatranscriptome, thus constituting the core of the prostate microbiome. The biodiversity of the microbiomes could not be differentiated between the matched tumour/benign specimens or between the tumour specimens of low and high Gleason Scores. The expression profile of ten Pseudomonas genes was strongly correlated with that of eight host small RNA genes; three of the RNA genes may negatively associate with metastasis. Few viruses could be identified from the prostate microbiomes. CONCLUSIONS: This is the first study of the human prostate microbiome employing an integrated metagenomics and metatranscriptomics approach. In this Chinese cohort, both metagenome and metatranscriptome analyses showed a non-sterile microenvironment in the prostate of PCa patients, but we did not find links between the microbiome and local progression of PCa. However, the correlated expression of Pseudomonas genes and human small RNA genes may provide tantalizing preliminary evidence that Pseudomonas infection may impede metastasis.

Characterization of transcriptomic signature of primary prostate cancer analogous to prostatic small cell neuroendocrine carcinoma.

Prostatic small cell neuroendocrine carcinoma (SC/NE) is well studied in metastatic castration-resistant prostate cancer; however, it is not well characterized in the primary setting. Herein, we used gene expression profiling of SC/NE prostate cancer (PCa) to develop a 212 gene signature to identify treatment-naïve primary prostatic tumors that are molecularly analogous to SC/NE (SC/NE-like PCa). The 212 gene signature was tested in several cohorts confirming similar molecular profile between prostatic SC/NE and small cell lung carcinoma. The signature was then translated into a genomic score (SCGScore) using modularized logistic regression modeling and validated in four independent cohorts achieving an average AUC >0.95. The signature was evaluated in more than 25,000 primary adenocarcinomas to characterize the biology, prognosis and potential therapeutic response of predicted SC/NE-like tumors. Assessing SCGScore in a prospective cohort of 17,967 RP and 6,697 biopsy treatment-naïve primary tumors from the Decipher Genomic Resource Information Database registry, approximately 1% of the patients were found to have a SC/NE-like transcriptional profile, whereas 0.5 and 3% of GG1 and GG5 patients respectively showed to be SC/NE-like. More than 80% of these patients are genomically high-risk based on Decipher score. Interrogating in vitro drug sensitivity analyses, SC/NE-like prostatic tumors showed higher response to PARP and HDAC inhibitors.

Metagenomic analysis reveals a rich bacterial content in high-risk prostate tumors from African men.

BACKGROUND: Inflammation is a hallmark of prostate cancer (PCa), yet no pathogenic agent has been identified. Men from Africa are at increased risk for both aggressive prostate disease and infection. We hypothesize that pathogenic microbes may be contributing, at least in part, to high-risk PCa presentation within Africa and in turn the observed ethnic disparity. METHODS: Here we reveal through metagenomic analysis of host-derived whole-genome sequencing data, the microbial content within prostate tumor tissue from 22 men. What is unique about this study is that patients were separated by ethnicity, African vs European, and environments, Africa vs Australia. RESULTS: We identified 23 common bacterial genera between the African, Australian, and Chinese prostate tumor samples, while nonbacterial microbes were notably absent. While the most abundant genera across all samples included: Escherichia, Propionibacterium, and Pseudomonas, the core prostate tumor microbiota was enriched for Proteobacteria. We observed a significant increase in the richness of the bacterial communities within the African vs Australian samples (t = 4.6-5.5; P = .0004-.001), largely driven by eight predominant genera. Considering core human gut microbiota, African prostate tissue samples appear enriched for Escherichia and Acidovorax, with an abundance of Eubacterium associated with host tumor hypermutation. CONCLUSIONS: Our study provides suggestive evidence for the presence of a core, bacteria-rich, prostate microbiome. While unable to exclude for fecal contamination, the observed increased bacterial content and richness within the African vs non-African samples, together with elevated tumor mutational burden, suggests the possibility that bacterially-driven oncogenic transformation within the prostate microenvironment may be contributing to aggressive disease presentation in Africa.

Computational identification of micro-structural variations and their proteogenomic consequences in cancer.

Motivation: Rapid advancement in high throughput genome and transcriptome sequencing (HTS) and mass spectrometry (MS) technologies has enabled the acquisition of the genomic, transcriptomic and proteomic data from the same tissue sample. We introduce a computational framework, ProTIE, to integratively analyze all three types of omics data for a complete molecular profile of a tissue sample. Our framework features MiStrVar, a novel algorithmic method to identify micro structural variants (microSVs) on genomic HTS data. Coupled with deFuse, a popular gene fusion detection method we developed earlier, MiStrVar can accurately profile structurally aberrant transcripts in tumors. Given the breakpoints obtained by MiStrVar and deFuse, our framework can then identify all relevant peptides that span the breakpoint junctions and match them with unique proteomic signatures. Observing structural aberrations in all three types of omics data validates their presence in the tumor samples. Results: We have applied our framework to all The Cancer Genome Atlas (TCGA) breast cancer Whole Genome Sequencing (WGS) and/or RNA-Seq datasets, spanning all four major subtypes, for which proteomics data from Clinical Proteomic Tumor Analysis Consortium (CPTAC) have been released. A recent study on this dataset focusing on SNVs has reported many that lead to novel peptides. Complementing and significantly broadening this study, we detected 244 novel peptides from 432 candidate genomic or transcriptomic sequence aberrations. Many of the fusions and microSVs we discovered have not been reported in the literature. Interestingly, the vast majority of these translated aberrations, fusions in particular, were private, demonstrating the extensive inter-genomic heterogeneity present in breast cancer. Many of these aberrations also have matching out-of-frame downstream peptides, potentially indicating novel protein sequence and structure. Availability and implementation: MiStrVar is available for download at https://bitbucket.org/compbio/mistrvar, and ProTIE is available at https://bitbucket.org/compbio/protie. Contact: cenksahi@indiana.edu. Supplementary information: Supplementary data are available at Bioinformatics online.

MechRNA: prediction of lncRNA mechanisms from RNA-RNA and RNA-protein interactions.

Motivation: Long non-coding RNAs (lncRNAs) are defined as transcripts longer than 200 nt that do not get translated into proteins. Often these transcripts are processed (spliced, capped and polyadenylated) and some are known to have important biological functions. However, most lncRNAs have unknown or poorly understood functions. Nevertheless, because of their potential role in cancer, lncRNAs are receiving a lot of attention, and the need for computational tools to predict their possible mechanisms of action is more than ever. Fundamentally, most of the known lncRNA mechanisms involve RNA-RNA and/or RNA-protein interactions. Through accurate predictions of each kind of interaction and integration of these predictions, it is possible to elucidate potential mechanisms for a given lncRNA. Results: Here, we introduce MechRNA, a pipeline for corroborating RNA-RNA interaction prediction and protein binding prediction for identifying possible lncRNA mechanisms involving specific targets or on a transcriptome-wide scale. The first stage uses a version of IntaRNA2 with added functionality for efficient prediction of RNA-RNA interactions with very long input sequences, allowing for large-scale analysis of lncRNA interactions with little or no loss of optimality. The second stage integrates protein binding information pre-computed by GraphProt, for both the lncRNA and the target. The final stage involves inferring the most likely mechanism for each lncRNA/target pair. This is achieved by generating candidate mechanisms from the predicted interactions, the relative locations of these interactions and correlation data, followed by selection of the most likely mechanistic explanation using a combined P-value. We applied MechRNA on a number of recently identified cancer-related lncRNAs (PCAT1, PCAT29 and ARLnc1) and also on two well-studied lncRNAs (PCA3 and 7SL). This led to the identification of hundreds of high confidence potential targets for each lncRNA and corresponding mechanisms. These predictions include the known competitive mechanism of 7SL with HuR for binding on the tumor suppressor TP53, as well as mechanisms expanding what is known about PCAT1 and ARLn1 and their targets BRCA2 and AR, respectively. For PCAT1-BRCA2, the mechanism involves competitive binding with HuR, which we confirmed using HuR immunoprecipitation assays. Availability and implementation: MechRNA is available for download at https://bitbucket.org/compbio/mechrna. Supplementary information: Supplementary data are available at Bioinformatics online.

Markers of MEK inhibitor resistance in low-grade serous ovarian cancer: EGFR is a potential therapeutic target.

BACKGROUND: Although low-grade serous ovarian cancer (LGSC) is rare, case-fatality rates are high as most patients present with advanced disease and current cytotoxic therapies are not overly effective. Recognizing that these cancers may be driven by MAPK pathway activation, MEK inhibitors (MEKi) are being tested in clinical trials. LGSC respond to MEKi only in a subgroup of patients, so predictive biomarkers and better therapies will be needed. METHODS: We evaluated a number of patient-derived LGSC cell lines, previously classified according to their MEKi sensitivity. Two cell lines were genomically compared against their matching tumors samples. MEKi-sensitive and MEKi-resistant lines were compared using whole exome sequencing and reverse phase protein array. Two treatment combinations targeting MEKi resistance markers were also evaluated using cell proliferation, cell viability, cell signaling, and drug synergism assays. RESULTS: Low-grade serous ovarian cancer cell lines recapitulated the genomic aberrations from their matching tumor samples. We identified three potential predictive biomarkers that distinguish MEKi sensitive and resistant lines: KRAS mutation status, and EGFR and PKC-alpha protein expression. The biomarkers were validated in three newly developed LGSC cell lines. Sub-lethal combination of MEK and EGFR inhibition showed drug synergy and caused complete cell death in two of four MEKi-resistant cell lines tested. CONCLUSIONS: KRAS mutations and the protein expression of EGFR and PKC-alpha should be evaluated as predictive biomarkers in patients with LGSC treated with MEKi. Combination therapy using a MEKi with EGFR inhibition may represent a promising new therapy for patients with MEKi-resistant LGSC.

Aneustat (OMN54) has aerobic glycolysis-inhibitory activity and also immunomodulatory activity as indicated by a first-generation PDX prostate cancer model.

Aneustat (OMN54) is a multivalent, botanical anticancer candidate therapeutic. A recent Phase-I clinical trial has indicated that it is well tolerated by patients and has immunomodulatory activity. In our study, using in vitro and in vivo prostate cancer models, we investigated Aneustat with regard to effects on (i) cancer-generated immunosuppression based on aerobic glycolysis leading to acidification of the tumor microenvironment, and (ii) immune-related processes such as macrophage differentiation and shifts in the intratumoral levels of host immune cells. Aneustat markedly reduced glucose consumption, lactic acid secretion, glycolysis-related gene expressions and proliferation of human LNCaP prostate cancer cells. In addition, Aneustat induced differentiation of RAW264.7 macrophages to the M1 anticancer phenotype. Treatment of LNCaP xenografts and first-generation patient-derived prostate cancer tissue xenografts with Aneustat in both cases led to a marked shift in intratumoral host (mouse/patient) immune cell levels: a higher ratio of cytotoxic CD8+ T/Treg cells, higher numbers of NK cells, lower numbers of Treg cells and MDSCs, i.e. changes favoring the host anticancer immune response. Taken together, the data indicate that Aneustat has immunomodulatory activity based on inhibition of aerobic glycolysis which in patients may lead to reduction of cancer-induced immunosuppression. Furthermore, first-generation patient-derived cancer tissue xenograft models may be used for screening compounds for immunomodulatory activity.