Comprehensive Molecular Characterization of Muscle-Invasive Bladder Cancer.

We report a comprehensive analysis of 412 muscle-invasive bladder cancers characterized by multiple TCGA analytical platforms. Fifty-eight genes were significantly mutated, and the overall mutational load was associated with APOBEC-signature mutagenesis. Clustering by mutation signature identified a high-mutation subset with 75% 5-year survival. mRNA expression clustering refined prior clustering analyses and identified a poor-survival "neuronal" subtype in which the majority of tumors lacked small cell or neuroendocrine histology. Clustering by mRNA, long non-coding RNA (lncRNA), and miRNA expression converged to identify subsets with differential epithelial-mesenchymal transition status, carcinoma in situ scores, histologic features, and survival. Our analyses identified 5 expression subtypes that may stratify response to different treatments.

Bridge: A New Algorithm for Rooting Orthologous Genes in Large-Scale Evolutionary Analyses.

Orthology information has been used for searching patterns in high-dimensional data, allowing transferring functional information between species. The key concept behind this strategy is that orthologous genes share ancestry to some extent. While reconstructing the history of a single gene is feasible with the existing computational resources, the reconstruction of entire biological systems remains challenging. In this study, we present Bridge, a new algorithm designed to infer the evolutionary root of orthologous genes in large-scale evolutionary analyses. The Bridge algorithm infers the evolutionary root of a given gene based on the distribution of its orthologs in a species tree. The Bridge algorithm is implemented in R and can be used either to assess genetic changes across the evolutionary history of orthologous groups or to infer the onset of specific traits in a biological system.

spongEffects: ceRNA modules offer patient-specific insights into the miRNA regulatory landscape.

MOTIVATION: Cancer is one of the leading causes of death worldwide. Despite significant improvements in prevention and treatment, mortality remains high for many cancer types. Hence, innovative methods that use molecular data to stratify patients and identify biomarkers are needed. Promising biomarkers can also be inferred from competing endogenous RNA (ceRNA) networks that capture the gene-miRNA gene regulatory landscape. Thus far, the role of these biomarkers could only be studied globally but not in a sample-specific manner. To mitigate this, we introduce spongEffects, a novel method that infers subnetworks (or modules) from ceRNA networks and calculates patient- or sample-specific scores related to their regulatory activity. RESULTS: We show how spongEffects can be used for downstream interpretation and machine learning tasks such as tumor classification and for identifying subtype-specific regulatory interactions. In a concrete example of breast cancer subtype classification, we prioritize modules impacting the biology of the different subtypes. In summary, spongEffects prioritizes ceRNA modules as biomarkers and offers insights into the miRNA regulatory landscape. Notably, these module scores can be inferred from gene expression data alone and can thus be applied to cohorts where miRNA expression information is lacking. AVAILABILITY AND IMPLEMENTATION: https://bioconductor.org/packages/devel/bioc/html/SPONGE.html.

TreeAndLeaf: an R/Bioconductor package for graphs and trees with focus on the leaves.

MOTIVATION: Dendrogram is a classical diagram for visualizing binary trees. Although efficient to represent hierarchical relations, it provides limited space for displaying information on the leaf elements, especially for large trees. RESULTS: Here, we present TreeAndLeaf, an R/Bioconductor package that implements a hybrid layout strategy to represent tree diagrams with focus on the leaves. The TreeAndLeaf package combines force-directed graph and tree layout algorithms using a single visualization system, allowing projection of multiple layers of information onto a graph-tree diagram. The Supplementary Information provides two case studies that use breast cancer data from epidemiological and experimental studies. AVAILABILITY AND IMPLEMENTATION: TreeAndLeaf is written in the R language, and is available from the Bioconductor project at http://bioconductor.org/packages/TreeAndLeaf/ (version≥1.4.2). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Novel lncRNAs Co-Expression Networks Identifies LINC00504 with Oncogenic Role in Luminal A Breast Cancer Cells.

Long non-coding RNAs (lncRNAs) are functional transcripts with more than 200 nucleotides. These molecules exhibit great regulatory capacity and may act at different levels of gene expression regulation. Despite this regulatory versatility, the biology of these molecules is still poorly understood. Computational approaches are being increasingly used to elucidate biological mechanisms in which these lncRNAs may be involved. Co-expression networks can serve as great allies in elucidating the possible regulatory contexts in which these molecules are involved. Herein, we propose the use of the pipeline deposited in the RTN package to build lncRNAs co-expression networks using TCGA breast cancer (BC) cohort data. Worldwide, BC is the most common cancer in women and has great molecular heterogeneity. We identified an enriched co-expression network for the validation of relevant cell processes in the context of BC, including LINC00504. This lncRNA has increased expression in luminal subtype A samples, and is associated with prognosis in basal-like subtype. Silencing this lncRNA in luminal A cell lines resulted in decreased cell viability and colony formation. These results highlight the relevance of the proposed method for the identification of lncRNAs in specific biological contexts.

RTNsurvival: an R/Bioconductor package for regulatory network survival analysis.

MOTIVATION: Transcriptional networks are models that allow the biological state of cells or tumours to be described. Such networks consist of connected regulatory units known as regulons, each comprised of a regulator and its targets. Inferring a transcriptional network can be a helpful initial step in characterizing the different phenotypes within a cohort. While the network itself provides no information on molecular differences between samples, the per-sample state of each regulon, i.e. the regulon activity, can be used for describing subtypes in a cohort. Integrating regulon activities with clinical data and outcomes would extend this characterization of differences between subtypes. RESULTS: We describe RTNsurvival, an R/Bioconductor package that calculates regulon activity profiles using transcriptional networks reconstructed by the RTN package, gene expression data, and a two-tailed Gene Set Enrichment Analysis. Given regulon activity profiles across a cohort, RTNsurvival can perform Kaplan-Meier analyses and Cox Proportional Hazards regressions, while also considering confounding variables. The Supplementary Information provides two case studies that use data from breast and liver cancer cohorts and features uni- and multivariate regulon survival analysis. AVAILABILITY AND IMPLEMENTATION: RTNsurvival is written in the R language, and is available from the Bioconductor project at http://bioconductor.org/packages/RTNsurvival/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

RTNduals: an R/Bioconductor package for analysis of co-regulation and inference of dual regulons.

MOTIVATION: Transcription factors (TFs) are key regulators of gene expression, and can activate or repress multiple target genes, forming regulatory units, or regulons. Understanding downstream effects of these regulators includes evaluating how TFs cooperate or compete within regulatory networks. Here we present RTNduals, an R/Bioconductor package that implements a general method for analyzing pairs of regulons. RESULTS: RTNduals identifies a dual regulon when the number of targets shared between a pair of regulators is statistically significant. The package extends the RTN (Reconstruction of Transcriptional Networks) package, and uses RTN transcriptional networks to identify significant co-regulatory associations between regulons. The Supplementary Information reports two case studies for TFs using the METABRIC and TCGA breast cancer cohorts. AVAILABILITY AND IMPLEMENTATION: RTNduals is written in the R language, and is available from the Bioconductor project at http://bioconductor.org/packages/RTNduals/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

FGFR2 risk SNPs confer breast cancer risk by augmenting oestrogen responsiveness.

The fibroblast growth factor receptor 2 (FGFR2) locus is consistently the top hit in genome-wide association studies for oestrogen receptor-positive (ER(+)) breast cancer. Yet, its mode of action continues to be controversial. Here, we employ a systems biology approach to demonstrate that signalling via FGFR2 counteracts cell activation by oestrogen. In the presence of oestrogen, the oestrogen receptor (ESR1) regulon (set of ESR1 target genes) is in an active state. However, signalling by FGFR2 is able to reverse the activity of the ESR1 regulon. This effect is seen in multiple distinct FGFR2 signalling model systems, across multiple cells lines and is dependent on the presence of FGFR2. Increased oestrogen exposure has long been associated with an increased risk of breast cancer. We therefore hypothesized that risk variants should reduce FGFR2 expression and subsequent signalling. Indeed, transient transfection experiments assaying the three independent variants of the FGFR2 risk locus (rs2981578, rs35054928 and rs45631563) in their normal chromosomal context show that these single-nucleotide polymorphisms (SNPs) map to transcriptional silencer elements and that, compared with wild type, the risk alleles augment silencer activity. The presence of risk variants results in lower FGFR2 expression and increased oestrogen responsiveness. We thus propose a molecular mechanism by which FGFR2 can confer increased breast cancer risk that is consistent with oestrogen exposure as a major driver of breast cancer risk. Our findings may have implications for the clinical use of FGFR2 inhibitors.

Reduced Neurite Density in Neuronal Cell Cultures Exposed to Serum of Patients with Bipolar Disorder.

BACKGROUND: Increased inflammatory markers and oxidative stress have been reported in serum among patients with bipolar disorder (BD). The aim of this study is to assess whether biochemical changes in the serum of patients induces neurotoxicity in neuronal cell cultures. METHODS: We challenged the retinoic acid-differentiated human neuroblastoma SH-SY5Y cells with the serum of BD patients at early and late stages of illness and assessed neurite density and cell viability as neurotoxic endpoints. RESULTS: Decreased neurite density was found in neurons treated with the serum of patients, mostly patients at late stages of illness. Also, neurons challenged with the serum of late-stage patients showed a significant decrease in cell viability. CONCLUSIONS: Our findings showed that the serum of patients with bipolar disorder induced a decrease in neurite density and cell viability in neuronal cultures.

GALANT: a Cytoscape plugin for visualizing data as functional landscapes projected onto biological networks.

SUMMARY: Network-level visualization of functional data is a key aspect of both analysis and understanding of biological systems. In a continuing effort to create clear and integrated visualizations that facilitate the gathering of novel biological insights despite the overwhelming complexity of data, we present here the GrAph LANdscape VisualizaTion (GALANT), a Cytoscape plugin that builds functional landscapes onto biological networks. By using GALANT, it is possible to project any type of numerical data onto a network to create a smoothed data map resembling the network layout. As a Cytoscape plugin, GALANT is further improved by the functionalities of Cytoscape, the popular bioinformatics package for biological network visualization and data integration. AVAILABILITY: http://www.lbbc.ibb.unesp.br/galant.

Towards a genome-wide transcriptogram: the Saccharomyces cerevisiae case.

Analysis of genome-wide expression data poses a challenge to extract relevant information. The usual approaches compare cellular expression levels relative to a pre-established control and genes are clustered based on the correlation of their expression levels. This implies that cluster definitions are dependent on the cellular metabolic state, eventually varying from one experiment to another. We present here a computational method that order genes on a line and clusters genes by the probability that their products interact. Protein-protein association information can be obtained from large data bases as STRING. The genome organization obtained this way is independent from specific experiments, and defines functional modules that are associated with gene ontology terms. The starting point is a gene list and a matrix specifying interactions. Considering the Saccharomyces cerevisiae genome, we projected on the ordering gene expression data, producing plots of transcription levels for two different experiments, whose data are available at Gene Expression Omnibus database. These plots discriminate metabolic cellular states, point to additional conclusions, and may be regarded as the first versions of 'transcriptograms'. This method is useful for extracting information from cell stimuli/responses experiments, and may be applied with diagnostic purposes to different organisms.

Expression-based subtypes define pathologic response to neoadjuvant immune-checkpoint inhibitors in muscle-invasive bladder cancer.

Checkpoint immunotherapy (CPI) has increased survival for some patients with advanced-stage bladder cancer (BCa). However, most patients do not respond. Here, we characterized the tumor and immune microenvironment in pre- and post-treatment tumors from the PURE01 neoadjuvant pembrolizumab immunotherapy trial, using a consolidative approach that combined transcriptional and genetic profiling with digital spatial profiling. We identify five distinctive genetic and transcriptomic programs and validate these in an independent neoadjuvant CPI trial to identify the features of response or resistance to CPI. By modeling the regulatory network, we identify the histone demethylase KDM5B as a repressor of tumor immune signaling pathways in one resistant subtype (S1, Luminal-excluded) and demonstrate that inhibition of KDM5B enhances immunogenicity in FGFR3-mutated BCa cells. Our study identifies signatures associated with response to CPI that can be used to molecularly stratify patients and suggests therapeutic alternatives for subtypes with poor response to neoadjuvant immunotherapy.

A transcriptional network of cell cycle dysregulation in noninvasive papillary urothelial carcinoma.

Human cancers display a restricted set of expression profiles, despite diverse mutational drivers. This has led to the hypothesis that select sets of transcription factors act on similar target genes as an integrated network, buffering a tumor's transcriptional state. Noninvasive papillary urothelial carcinoma (NIPUC) with higher cell cycle activity has higher risk of recurrence and progression. In this paper, we describe a transcriptional network of cell cycle dysregulation in NIPUC, which was delineated using the ARACNe algorithm applied to expression data from a new cohort (n = 81, RNA sequencing), and two previously published cohorts. The transcriptional network comprised 121 transcription factors, including the pluripotency factors SOX2 and SALL4, the sex hormone binding receptors ESR1 and PGR, and multiple homeobox factors. Of these 121 transcription factors, 65 and 56 were more active in tumors with greater and less cell cycle activity, respectively. When clustered by activity of these transcription factors, tumors divided into High Cell Cycle versus Low Cell Cycle groups. Tumors in the High Cell Cycle group demonstrated greater mutational burden and copy number instability. A putative mutational driver of cell cycle dysregulation, such as homozygous loss of CDKN2A, was found in only 50% of High Cell Cycle NIPUC, suggesting a prominent role of transcription factor activity in driving cell cycle dysregulation. Activity of the 121 transcription factors strongly associated with expression of EZH2 and other members of the PRC2 complex, suggesting regulation by this complex influences expression of the transcription factors in this network. Activity of transcription factors in this network also associated with signatures of pluripotency and epithelial-to-mesenchymal transition (EMT), suggesting they play a role in driving evolution to invasive carcinoma. Consistent with this, these transcription factors differed in activity between NIPUC and invasive urothelial carcinoma.

Comprehensive Characterization of the Regulatory Landscape of Adrenocortical Carcinoma: Novel Transcription Factors and Targets Associated with Prognosis.

We reconstructed a transcriptional regulatory network for adrenocortical carcinoma (ACC) using transcriptomic and clinical data from The Cancer Genome Atlas (TCGA)-ACC cohort. We investigated the association of transcriptional regulatory units (regulons) with overall survival, molecular phenotypes, and immune signatures. We annotated the ACC regulons with cancer hallmarks and assessed single sample regulon activities in the European Network for the Study of Adrenal Tumors (ENSAT) cohort. We found 369 regulons associated with overall survival and subdivided them into four clusters: RC1 and RC2, associated with good prognosis, and RC3 and RC4, associated with worse outcomes. The RC1 and RC3 regulons were highly correlated with the 'Steroid Phenotype,' while the RC2 and RC4 regulons were highly correlated with a molecular proliferation signature. We selected two regulons, NR5A1 (steroidogenic factor 1, SF-1) and CENPA (Centromeric Protein A), that were consistently associated with overall survival for further downstream analyses. The CENPA regulon was the primary regulator of MKI-67 (a marker of proliferation KI-67), while the NR5A1 regulon is a well-described transcription factor (TF) in ACC tumorigenesis. We also found that the ZBTB4 (Zinc finger and BTB domain-containing protein 4) regulon, which is negatively associated with CENPA in our transcriptional regulatory network, is also a druggable anti-tumorigenic TF. We anticipate that the ACC regulons may be used as a reference for further investigations concerning the complex molecular interactions in ACC tumors.

FOXA1 repression drives lineage plasticity and immune heterogeneity in bladder cancers with squamous differentiation.

Cancers arising from the bladder urothelium often exhibit lineage plasticity with regions of urothelial carcinoma adjacent to or admixed with regions of divergent histomorphology, most commonly squamous differentiation. To define the biologic basis for and clinical significance of this morphologic heterogeneity, here we perform integrated genomic analyses of mixed histology bladder cancers with separable regions of urothelial and squamous differentiation. We find that squamous differentiation is a marker of intratumoral genomic and immunologic heterogeneity in patients with bladder cancer and a biomarker of intrinsic immunotherapy resistance. Phylogenetic analysis confirms that in all cases the urothelial and squamous regions are derived from a common shared precursor. Despite the presence of marked genomic heterogeneity between co-existent urothelial and squamous differentiated regions, no recurrent genomic alteration exclusive to the urothelial or squamous morphologies is identified. Rather, lineage plasticity in bladder cancers with squamous differentiation is associated with loss of expression of FOXA1, GATA3, and PPARG, transcription factors critical for maintenance of urothelial cell identity. Of clinical significance, lineage plasticity and PD-L1 expression is coordinately dysregulated via FOXA1, with patients exhibiting morphologic heterogeneity pre-treatment significantly less likely to respond to immune checkpoint inhibitors.

Unraveling Immune-Related lncRNAs in Breast Cancer Molecular Subtypes.

Breast cancer (BRCA) is the most leading cause of cancer worldwide. It is a heterogeneous disease with at least five molecular subtypes including luminal A, luminal B, basal-like, HER2-enriched, and normal-like. These five molecular subtypes are usually stratified according to their mRNA profile patterns; however, ncRNAs are increasingly being used for this purpose. Among the ncRNAs class, the long non-coding RNAs (lncRNAs) are molecules with more than 200 nucleotides with versatile regulatory roles; and high tissue-specific expression profiles. The heterogeneity of BRCA can also be reflected regarding tumor microenvironment immune cells composition, which can directly impact a patient's prognosis and therapy response. Using BRCA immunogenomics data from a previous study, we propose here a bioinformatics approach to include lncRNAs complexity in BRCA molecular and immune subtype. RNA-seq data from The Cancer Genome Atlas (TCGA) BRCA cohort was analyzed, and signal-to-noise ratio metrics were applied to create these subtype-specific signatures. Five immune-related signatures were generated with approximately ten specific lncRNAs, which were then functionally analyzed using GSEA enrichment and survival analysis. We highlighted here some lncRNAs in each subtype. LINC01871 is related to immune response activation and favorable overall survival in basal-like samples; EBLN3P is related to immune response suppression and progression in luminal B, MEG3, XXYLT1-AS2, and LINC02613 were related with immune response activation in luminal A, HER2-enriched and normal-like subtypes, respectively. In this way, we emphasize the need to know better the role of lncRNAs as regulators of immune response to provide new perspectives regarding diagnosis, prognosis and therapeutical targets in BRCA molecular subtypes.

Adrenocortical Carcinoma Steroid Profiles: In Silico Pan-Cancer Analysis of TCGA Data Uncovers Immunotherapy Targets for Potential Improved Outcomes.

Despite progress in understanding the biology of adrenocortical carcinoma (ACC), treatment options have not dramatically changed in the last three decades, nor have we learned how to avoid some of its long-term side effects. Our goal was to improve the understanding of immune pathways that may include druggable targets to enhance immune responses of patients with ACC, focusing on immune evasion and the activation of immune cells against ACC. Our strategy was aimed at improving insight regarding gene expression without steroid interference. Using approaches based on high and low steroid phenotypes (HSP and LSP, respectively), we characterized immune pathways using The Cancer Genome Atlas (TCGA) ACC cohort data. Although previous studies have suggested that patients with ACC receive minimal benefit from immunotherapy, high expression of immune modulators was noted in patients with LSP, suggesting the activation of these biomarkers may be an important adjuvant therapy target after clearance of excess glucocorticoids. In addition, patients with LSP ACC had higher immune cell infiltration than patients with HSP ACC and other cancer subtypes. Our findings can be summarized as follows (1): we confirmed and improved the definition of two immune response pathways to ACC (HSP and LSP) based on in silico transcriptome analysis (2), we demonstrated the steroid profile should be considered, otherwise analyses of ACC immune characteristics can generate confounding results (3), among the overexpressed immunotherapy targets, we demonstrated that LSP was rich in PDCD1LG2 (PD-L2) and both HSP and LSP overexpressed CD276 (B7-H3), which was associated with resistance to anti-PD1 therapy and may have accounted for the modest results of previous clinical trials, and (4) identification of patients with LSP or HSP ACC can be used to help determine whether immunotherapy should be used. In conclusion, we highlighted the differences between LSP and HSP, drawing attention to potential therapeutic targets (CD276, PDCD1, and PDCD1LG2). Treatments to reduce immune evasion, as well as the use of other natural and pharmacological immune activators, should include prior pharmacological inhibition of steroidogenesis. Attempts to combine these with tumor cell proliferation inhibitors, if they do not affect cells of the immune system, may produce interesting results.

AKT3 Expression in Mesenchymal Colorectal Cancer Cells Drives Growth and Is Associated with Epithelial-Mesenchymal Transition.

Colorectal cancer (CRC) is a heterogeneous disease that can currently be subdivided into four distinct consensus molecular subtypes (CMS) based on gene expression profiling. The CMS4 subtype is marked by high expression of mesenchymal genes and is associated with a worse overall prognosis compared to other CMSs. Importantly, this subtype responds poorly to the standard therapies currently used to treat CRC. We set out to explore what regulatory signalling networks underlie the CMS4 phenotype of cancer cells, specifically, by analysing which kinases were more highly expressed in this subtype compared to others. We found AKT3 to be expressed in the cancer cell epithelium of CRC specimens, patient derived xenograft (PDX) models and in (primary) cell cultures representing CMS4. Importantly, chemical inhibition or knockout of this gene hampers outgrowth of this subtype, as AKT3 controls expression of the cell cycle regulator p27KIP1. Furthermore, high AKT3 expression was associated with high expression of epithelial-mesenchymal transition (EMT) genes, and this observation could be expanded to cell lines representing other carcinoma types. More importantly, this association allowed for the identification of CRC patients with a high propensity to metastasise and an associated poor prognosis. High AKT3 expression in the tumour epithelial compartment may thus be used as a surrogate marker for EMT and may allow for a selection of CRC patients that could benefit from AKT3-targeted therapy.

Stage-stratified molecular profiling of non-muscle-invasive bladder cancer enhances biological, clinical, and therapeutic insight.

Understanding the molecular determinants that underpin the clinical heterogeneity of non-muscle-invasive bladder cancer (NMIBC) is essential for prognostication and therapy development. Stage T1 disease in particular presents a high risk of progression and requires improved understanding. We present a detailed multi-omics study containing gene expression, copy number, and mutational profiles that show relationships to immune infiltration, disease recurrence, and progression to muscle invasion. We compare expression and genomic subtypes derived from all NMIBCs with those derived from the individual disease stages Ta and T1. We show that sufficient molecular heterogeneity exists within the separate stages to allow subclassification and that this is more clinically meaningful for stage T1 disease than that derived from all NMIBCs. This provides improved biological understanding and identifies subtypes of T1 tumors that may benefit from chemo- or immunotherapy.