Molecular Classification of Colorectal Cancer by microRNA Profiling: Correlation with the Consensus Molecular Subtypes (CMS) and Validation of miR-30b Targets.

Colorectal cancer consensus molecular subtypes (CMSs) are widely accepted and constitutes the basis for patient stratification to improve clinical practice. We aimed to find whether miRNAs could reproduce molecular subtypes, and to identify miRNA targets associated to the High-stroma/CMS4 subtype. The expression of 939 miRNAs was analyzed in tumors classified in CMS. TALASSO was used to find gene-miRNA interactions. A miR-mRNA regulatory network was constructed using Cytoscape. Candidate gene-miR interactions were validated in 293T cells. Hierarchical-Clustering identified three miRNA tumor subtypes (miR-LS; miR-MI; and miR-HS) which were significantly associated (p < 0.001) to the reported mRNA subtypes. miR-LS correlated with the low-stroma/CMS2; miR-MI with the mucinous-MSI/CMS1 and miR-HS with high-stroma/CMS4. MicroRNA tumor subtypes and association to CMSs were validated with TCGA datasets. TALASSO identified 1462 interactions (p < 0.05) out of 21,615 found between 176 miRs and 788 genes. Based on the regulatory network, 88 miR-mRNA interactions were selected as candidates. This network was functionally validated for the pair miR-30b/SLC6A6. We found that miR-30b overexpression silenced 3'-UTR-SLC6A6-driven luciferase expression in 293T-cells; mutation of the target sequence in the 3'-UTR-SLC6A6 prevented the miR-30b inhibitory effect. In conclusion CRC subtype classification using a miR-signature might facilitate a real-time analysis of the disease course and treatment response.

Functional Enrichment Analysis of Regulatory Elements.

Statistical methods for enrichment analysis are important tools to extract biological information from omics experiments. Although these methods have been widely used for the analysis of gene and protein lists, the development of high-throughput technologies for regulatory elements demands dedicated statistical and bioinformatics tools. Here, we present a set of enrichment analysis methods for regulatory elements, including CpG sites, miRNAs, and transcription factors. Statistical significance is determined via a power weighting function for target genes and tested by the Wallenius noncentral hypergeometric distribution model to avoid selection bias. These new methodologies have been applied to the analysis of a set of miRNAs associated with arrhythmia, showing the potential of this tool to extract biological information from a list of regulatory elements. These new methods are available in GeneCodis 4, a web tool able to perform singular and modular enrichment analysis that allows the integration of heterogeneous information.

Genome-wide signatures of flowering adaptation to climate temperature: Regional analyses in a highly diverse native range of Arabidopsis thaliana.

Current global change is fueling an interest to understand the genetic and molecular mechanisms of plant adaptation to climate. In particular, altered flowering time is a common strategy for escape from unfavourable climate temperature. In order to determine the genomic bases underlying flowering time adaptation to this climatic factor, we have systematically analysed a collection of 174 highly diverse Arabidopsis thaliana accessions from the Iberian Peninsula. Analyses of 1.88 million single nucleotide polymorphisms provide evidence for a spatially heterogeneous contribution of demographic and adaptive processes to geographic patterns of genetic variation. Mountains appear to be allele dispersal barriers, whereas the relationship between flowering time and temperature depended on the precise temperature range. Environmental genome-wide associations supported an overall genome adaptation to temperature, with 9.4% of the genes showing significant associations. Furthermore, phenotypic genome-wide associations provided a catalogue of candidate genes underlying flowering time variation. Finally, comparison of environmental and phenotypic genome-wide associations identified known (Twin Sister of FT, FRIGIDA-like 1, and Casein Kinase II Beta chain 1) and new (Epithiospecifer Modifier 1 and Voltage-Dependent Anion Channel 5) genes as candidates for adaptation to climate temperature by altered flowering time. Thus, this regional collection provides an excellent resource to address the spatial complexity of climate adaptation in annual plants.

miRNA profiling during antigen-dependent T cell activation: A role for miR-132-3p.

microRNAs (miRNAs) are tightly regulated during T lymphocyte activation to enable the establishment of precise immune responses. Here, we analyzed the changes of the miRNA profiles of T cells in response to activation by cognate interaction with dendritic cells. We also studied mRNA targets common to miRNAs regulated in T cell activation. pik3r1 gene, which encodes the regulatory subunits of PI3K p50, p55 and p85, was identified as target of miRNAs upregulated after T cell activation. Using 3'UTR luciferase reporter-based and biochemical assays, we showed the inhibitory relationship between miR-132-3p upregulation and expression of the pik3r1 gene. Our results indicate that specific miRNAs whose expression is modulated during T cell activation might regulate PI3K signaling in T cells.

Gene signature associated with benign neurofibroma transformation to malignant peripheral nerve sheath tumors.

Benign neurofibromas, the main phenotypic manifestations of the rare neurological disorder neurofibromatosis type 1, degenerate to malignant tumors associated to poor prognosis in about 10% of patients. Despite efforts in the field of (epi)genomics, the lack of prognostic biomarkers with which to predict disease evolution frustrates the adoption of appropriate early therapeutic measures. To identify potential biomarkers of malignant neurofibroma transformation, we integrated four human experimental studies and one for mouse, using a gene score-based meta-analysis method, from which we obtained a score-ranked signature of 579 genes. Genes with the highest absolute scores were classified as promising disease biomarkers. By grouping genes with similar neurofibromatosis-related profiles, we derived panels of potential biomarkers. The addition of promoter methylation data to gene profiles indicated a panel of genes probably silenced by hypermethylation. To identify possible therapeutic treatments, we used the gene signature to query drug expression databases. Trichostatin A and other histone deacetylase inhibitors, as well as cantharidin and tamoxifen, were retrieved as putative therapeutic means to reverse the aberrant regulation that drives to malignant cell proliferation and metastasis. This in silico prediction corroborated reported experimental results that suggested the inclusion of these compounds in clinical trials. This experimental validation supported the suitability of the meta-analysis method used to integrate several sources of public genomic information, and the reliability of the gene signature associated to the malignant evolution of neurofibromas to generate working hypotheses for prognostic and drug-responsive biomarkers or therapeutic measures, thus showing the potential of this in silico approach for biomarker discovery.

Distinct Molecular Signature of Murine Fetal Liver and Adult Hematopoietic Stem Cells Identify Novel Regulators of Hematopoietic Stem Cell Function.

During ontogeny, fetal liver (FL) acts as a major site for hematopoietic stem cell (HSC) maturation and expansion, whereas HSCs in the adult bone marrow (ABM) are largely quiescent. HSCs in the FL possess faster repopulation capacity as compared with ABM HSCs. However, the molecular mechanism regulating the greater self-renewal potential of FL HSCs has not yet extensively been assessed. Recently, we published RNA sequencing-based gene expression analysis on FL HSCs from 14.5-day mouse embryo (E14.5) in comparison to the ABM HSCs. We reanalyzed these data to identify key transcriptional regulators that play important roles in the expansion of HSCs during development. The comparison of FL E14.5 with ABM HSCs identified more than 1,400 differentially expressed genes. More than 200 genes were shortlisted based on the gene ontology (GO) annotation term "transcription." By morpholino-based knockdown studies in zebrafish, we assessed the function of 18 of these regulators, previously not associated with HSC proliferation. Our studies identified a previously unknown role for tdg, uhrf1, uchl5, and ncoa1 in the emergence of definitive hematopoiesis in zebrafish. In conclusion, we demonstrate that identification of genes involved in transcriptional regulation differentially expressed between expanding FL HSCs and quiescent ABM HSCs, uncovers novel regulators of HSC function.

Notch-regulated miR-223 targets the aryl hydrocarbon receptor pathway and increases cytokine production in macrophages from rheumatoid arthritis patients.

Evidence links aryl hydrocarbon receptor (AHR) activation to rheumatoid arthritis (RA) pathogenesis, although results are inconsistent. AHR agonists inhibit pro-inflammatory cytokine expression in macrophages, pivotal cells in RA aetiopathogenesis, which hints at specific circuits that regulate the AHR pathway in RA macrophages. We compared microRNA (miR) expression in CD14(+) cells from patients with active RA or with osteoarthritis (OA). Seven miR were downregulated and one (miR-223) upregulated in RA compared to OA cells. miR-223 upregulation correlated with reduced Notch3 and Notch effector expression in RA patients. Overexpression of the Notch-induced repressor HEY-1 and co-culture of healthy donor monocytes with Notch ligand-expressing cells showed direct Notch-mediated downregulation of miR-223. Bioinformatics predicted the AHR regulator ARNT (AHR nuclear translocator) as a miR-223 target. Pre-miR-223 overexpression silenced ARNT 3'UTR-driven reporter expression, reduced ARNT (but not AHR) protein levels and prevented AHR/ARNT-mediated inhibition of pro-inflammatory cytokine expression. miR-223 counteracted AHR/ARNT-induced Notch3 upregulation in monocytes. Levels of ARNT and of CYP1B1, an AHR/ARNT signalling effector, were reduced in RA compared to OA synovial tissue, which correlated with miR-223 levels. Our results associate Notch signalling to miR-223 downregulation in RA macrophages, and identify miR-223 as a negative regulator of the AHR/ARNT pathway through ARNT targeting.

Highly proliferative primitive fetal liver hematopoietic stem cells are fueled by oxidative metabolic pathways.

Hematopoietic stem cells (HSCs) in the fetal liver (FL) unlike adult bone marrow (BM) proliferate extensively, posing different metabolic demands. However, metabolic pathways responsible for the production of energy and cellular building blocks in FL HSCs have not been described. Here, we report that FL HSCs use oxygen dependent energy generating pathways significantly more than their BM counterparts. RNA-Seq analysis of E14.5 FL versus BM derived HSCs identified increased expression levels of genes involved in oxidative phosphorylation (OxPhos) and the citric acid cycle (TCA). We demonstrated that FL HSCs contain more mitochondria than BM HSCs, which resulted in increased levels of oxygen consumption and reactive oxygen species (ROS) production. Higher levels of DNA repair and antioxidant pathway gene expression may prevent ROS-mediated (geno)toxicity in FL HSCs. Thus, we here for the first time highlight the underestimated importance of oxygen dependent pathways for generating energy and building blocks in FL HSCs.

Proteogenomics Dashboard for the Human Proteome Project.

dasHPPboard is a novel proteomics-based dashboard that collects and reports the experiments produced by the Spanish Human Proteome Project consortium (SpHPP) and aims to help HPP to map the entire human proteome. We have followed the strategy of analog genomics projects like the Encyclopedia of DNA Elements (ENCODE), which provides a vast amount of data on human cell lines experiments. The dashboard includes results of shotgun and selected reaction monitoring proteomics experiments, post-translational modifications information, as well as proteogenomics studies. We have also processed the transcriptomics data from the ENCODE and Human Body Map (HBM) projects for the identification of specific gene expression patterns in different cell lines and tissues, taking special interest in those genes having little proteomic evidence available (missing proteins). Peptide databases have been built using single nucleotide variants and novel junctions derived from RNA-Seq data that can be used in search engines for sample-specific protein identifications on the same cell lines or tissues. The dasHPPboard has been designed as a tool that can be used to share and visualize a combination of proteomic and transcriptomic data, providing at the same time easy access to resources for proteogenomics analyses. The dasHPPboard can be freely accessed at: http://sphppdashboard.cnb.csic.es.

Quest for Missing Proteins: Update 2015 on Chromosome-Centric Human Proteome Project.

This paper summarizes the recent activities of the Chromosome-Centric Human Proteome Project (C-HPP) consortium, which develops new technologies to identify yet-to-be annotated proteins (termed "missing proteins") in biological samples that lack sufficient experimental evidence at the protein level for confident protein identification. The C-HPP also aims to identify new protein forms that may be caused by genetic variability, post-translational modifications, and alternative splicing. Proteogenomic data integration forms the basis of the C-HPP's activities; therefore, we have summarized some of the key approaches and their roles in the project. We present new analytical technologies that improve the chemical space and lower detection limits coupled to bioinformatics tools and some publicly available resources that can be used to improve data analysis or support the development of analytical assays. Most of this paper's content has been compiled from posters, slides, and discussions presented in the series of C-HPP workshops held during 2014. All data (posters, presentations) used are available at the C-HPP Wiki (http://c-hpp.webhosting.rug.nl/) and in the Supporting Information.

NFFinder: an online bioinformatics tool for searching similar transcriptomics experiments in the context of drug repositioning.

Drug repositioning, using known drugs for treating conditions different from those the drug was originally designed to treat, is an important drug discovery tool that allows for a faster and cheaper development process by using drugs that are already approved or in an advanced trial stage for another purpose. This is especially relevant for orphan diseases because they affect too few people to make drug research de novo economically viable. In this paper we present NFFinder, a bioinformatics tool for identifying potential useful drugs in the context of orphan diseases. NFFinder uses transcriptomic data to find relationships between drugs, diseases and a phenotype of interest, as well as identifying experts having published on that domain. The application shows in a dashboard a series of graphics and tables designed to help researchers formulate repositioning hypotheses and identify potential biological relationships between drugs and diseases. NFFinder is freely available at http://nffinder.cnb.csic.es.

NMF-mGPU: non-negative matrix factorization on multi-GPU systems.

BACKGROUND: In the last few years, the Non-negative Matrix Factorization ( NMF ) technique has gained a great interest among the Bioinformatics community, since it is able to extract interpretable parts from high-dimensional datasets. However, the computing time required to process large data matrices may become impractical, even for a parallel application running on a multiprocessors cluster. In this paper, we present NMF-mGPU, an efficient and easy-to-use implementation of the NMF algorithm that takes advantage of the high computing performance delivered by Graphics-Processing Units ( GPUs ). Driven by the ever-growing demands from the video-games industry, graphics cards usually provided in PCs and laptops have evolved from simple graphics-drawing platforms into high-performance programmable systems that can be used as coprocessors for linear-algebra operations. However, these devices may have a limited amount of on-board memory, which is not considered by other NMF implementations on GPU. RESULTS: NMF-mGPU is based on CUDA ( Compute Unified Device Architecture ), the NVIDIA's framework for GPU computing. On devices with low memory available, large input matrices are blockwise transferred from the system's main memory to the GPU's memory, and processed accordingly. In addition, NMF-mGPU has been explicitly optimized for the different CUDA architectures. Finally, platforms with multiple GPUs can be synchronized through MPI ( Message Passing Interface ). In a four-GPU system, this implementation is about 120 times faster than a single conventional processor, and more than four times faster than a single GPU device (i.e., a super-linear speedup). CONCLUSIONS: Applications of GPUs in Bioinformatics are getting more and more attention due to their outstanding performance when compared to traditional processors. In addition, their relatively low price represents a highly cost-effective alternative to conventional clusters. In life sciences, this results in an excellent opportunity to facilitate the daily work of bioinformaticians that are trying to extract biological meaning out of hundreds of gigabytes of experimental information. NMF-mGPU can be used "out of the box" by researchers with little or no expertise in GPU programming in a variety of platforms, such as PCs, laptops, or high-end GPU clusters. NMF-mGPU is freely available at https://github.com/bioinfo-cnb/bionmf-gpu .

A subpopulation of CD103(pos) ICOS(pos) Treg cells occurs at high frequency in lymphopenic mice and represents a lymph node specific differentiation stage.

Regulatory T (Treg) cells are pivotal for the maintenance of peripheral tolerance by controlling self-reactive, chronic, and homeostatic T-cell responses. Here, we report that the increase in Treg-cell suppressive function observed in lymphopenic mice correlates with the degree of lymphopenia and is caused by a higher frequency of a novel subpopulation of CD103(pos) ICOS(pos) Treg cells. Though present in the thymus, CD103(pos) ICOS(pos) Treg cells are not generated there but recirculate from the periphery to that site. The acquisition and maintenance of this distinctive phenotype requires the LN microenvironment and the in situ availability of antigen. Contrary to conventional effector and other Treg cells, the cellularity of CD103(pos) ICOS(pos) Treg cells is not affected by the absence of IL-7 and thymic stroma lymphopoetin. Given their increased frequency in lymphopenia, the absolute number of CD103(pos) ICOS(pos) Treg cells remains unchanged in the periphery irrespective of a paucity of total Treg cells. We furthermore demonstrate, with cell transfers in mice, that the CD103(pos) ICOS(pos) phenotype represents a LN-specific differentiation stage arrived at by several other Treg-cell subsets. Thus, tissue-specific cues determine the overall potency of the peripheral Treg-cell pool by shaping its subset composition.

miR-217 is an oncogene that enhances the germinal center reaction.

microRNAs are a class of regulators of gene expression that have been shown critical for a great number of biological processes; however, little is known of their role in germinal center (GC) B cells. Although the GC reaction is crucial to ensure a competent immune response, GC B cells are also the origin of most human lymphomas, presumably due to bystander effects of the immunoglobulin gene remodeling that takes place at these sites. Here we report that miR-217 is specifically upregulated in GC B cells. Gain- and loss-of-function mouse models reveal that miR-217 is a positive modulator of the GC response that increases the generation of class-switched antibodies and the frequency of somatic hypermutation. We find that miR-217 down-regulates the expression of a DNA damage response and repair gene network and in turn stabilizes Bcl-6 expression in GC B cells. Importantly, miR-217 overexpression also promotes mature B-cell lymphomagenesis; this is physiologically relevant as we find that miR-217 is overexpressed in aggressive human B-cell lymphomas. Therefore, miR-217 provides a novel molecular link between the normal GC response and B-cell transformation.

Improving miRNA-mRNA interaction predictions.

BACKGROUND: MicroRNAs are short RNA molecules that post-transcriptionally regulate gene expression. Today, microRNA target prediction remains challenging since very few have been experimentally validated and sequence-based predictions have large numbers of false positives. Furthermore, due to the different measuring rules used in each database of predicted interactions, the selection of the most reliable ones requires extensive knowledge about each algorithm. RESULTS: Here we propose two methods to measure the confidence of predicted interactions based on experimentally validated information. The output of the methods is a combined database where new scores and statistical confidences are re-assigned to each predicted interaction. The new scores allow the robust combination of several databases without the effect of low-performing algorithms dragging down good-performing ones. The combined databases obtained using both algorithms described in this paper outperform each of the existing predictive algorithms that were considered for the combination. CONCLUSIONS: Our approaches are a useful way to integrate predicted interactions from different databases. They reduce the selection of interactions to a unique database based on an intuitive score and allow comparing databases between them.

CentrosomeDB: a new generation of the centrosomal proteins database for Human and Drosophila melanogaster.

We present the second generation of centrosomeDB, available online at http://centrosome.cnb.csic.es, with a significant expansion of 1357 human and drosophila centrosomal genes and their corresponding information. The centrosome of animal cells takes part in important biological processes such as the organization of the interphase microtubule cytoskeleton and the assembly of the mitotic spindle. The active research done during the past decades has produced lots of data related to centrosomal proteins. Unfortunately, the accumulated data are dispersed among diverse and heterogeneous sources of information. We believe that the availability of a repository collecting curated evidences of centrosomal proteins would constitute a key resource for the scientific community. This was our first motivation to introduce CentrosomeDB in NAR database issue in 2009, collecting a set of human centrosomal proteins that were reported in the literature and other sources. The intensive use of this resource during these years has encouraged us to present this new expanded version. Using our database, the researcher is offered the possibility to study the evolution, function and structure of the centrosome. We have compiled information from many sources, including Gene Ontology, disease-association, single nucleotide polymorphisms and associated gene expression experiments. Special interest has been paid to protein-protein interaction.

Surfing transcriptomic landscapes. A step beyond the annotation of chromosome 16 proteome.

The Spanish team of the Human Proteome Project (SpHPP) marked the annotation of Chr16 and data analysis as one of its priorities. Precise annotation of Chromosome 16 proteins according to C-HPP criteria is presented. Moreover, Human Body Map 2.0 RNA-Seq and Encyclopedia of DNA Elements (ENCODE) data sets were used to obtain further information relative to cell/tissue specific chromosome 16 coding gene expression patterns and to infer the presence of missing proteins. Twenty-four shotgun 2D-LC-MS/MS and gel/LC-MS/MS MIAPE compliant experiments, representing 41% coverage of chromosome 16 proteins, were performed. Furthermore, mapping of large-scale multicenter mass spectrometry data sets from CCD18, MCF7, Jurkat, and Ramos cell lines into RNA-Seq data allowed further insights relative to correlation of chromosome 16 transcripts and proteins. Detection and quantification of chromosome 16 proteins in biological matrices by SRM procedures are also primary goals of the SpHPP. Two strategies were undertaken: one focused on known proteins, taking advantage of MS data already available, and the second, aimed at the detection of the missing proteins, is based on the expression of recombinant proteins to gather MS information and optimize SRM methods that will be used in real biological samples. SRM methods for 49 known proteins and for recombinant forms of 24 missing proteins are reported in this study.

Sumoylated hnRNPA2B1 controls the sorting of miRNAs into exosomes through binding to specific motifs.

Exosomes are released by most cells to the extracellular environment and are involved in cell-to-cell communication. Exosomes contain specific repertoires of mRNAs, microRNAs (miRNAs) and other non-coding RNAs that can be functionally transferred to recipient cells. However, the mechanisms that control the specific loading of RNA species into exosomes remain unknown. Here we describe sequence motifs present in miRNAs that control their localization into exosomes. The protein heterogeneous nuclear ribonucleoprotein A2B1 (hnRNPA2B1) specifically binds exosomal miRNAs through the recognition of these motifs and controls their loading into exosomes. Moreover, hnRNPA2B1 in exosomes is sumoylated, and sumoylation controls the binding of hnRNPA2B1 to miRNAs. The loading of miRNAs into exosomes can be modulated by mutagenesis of the identified motifs or changes in hnRNPA2B1 expression levels. These findings identify hnRNPA2B1 as a key player in miRNA sorting into exosomes and provide potential tools for the packaging of selected regulatory RNAs into exosomes and their use in biomedical applications.

Unraveling a novel transcription factor code determining the human arterial-specific endothelial cell signature.

Endothelial cells (ECs) lining arteries and veins have distinct molecular/functional signatures. The underlying regulatory mechanisms are incompletely understood. Here, we established a specific fingerprint of freshly isolated arterial and venous ECs from human umbilical cord comprising 64 arterial and 12 venous genes, representing distinct functions/pathways. Among the arterial genes were 8 transcription factors (TFs), including Notch target HEY2, the current "gold standard" determinant for arterial EC (aEC) specification. Culture abrogated differential gene expression in part due to gradual loss of canonical Notch activity and HEY2 expression. Notably, restoring HEY2 expression or Delta-like4-induced Notch signaling in cultured ECs only partially reinstated the aEC gene signature, whereas combined overexpression of the 8 TFs restored this fingerprint more robustly. Whereas some TFs stimulated few genes, others boosted a large proportion of arterial genes. Although there was some overlap and cross-regulation, the TFs largely complemented each other in regulating the aEC gene profile. Finally, overexpression of the 8 TFs in human umbilical vein ECs conveyed an arterial-like behavior upon their implantation in a Matrigel plug in vivo. Thus, our study shows that Notch signaling determines only part of the aEC signature and identifies additional novel and complementary transcriptional players in the complex regulation of human arteriovenous EC identity.

Chromothripsis: breakage-fusion-bridge over and over again.

The acquisition of massive but localized chromosome translocations, a phenomenon termed chromothripsis, has received widespread attention since its discovery over a year ago. Until recently, chromothripsis was believed to originate from a single catastrophic event, but the molecular mechanisms leading to this event are yet to be uncovered. Because a thorough interpretation of the data are missing, the phenomenon itself has wrongly acquired the status of a mechanism used to justify many kinds of complex rearrangements. Although the assumption that all translocations in chromothripsis originate from a single event has met with criticism, satisfactory explanations for the intense but localized nature of this phenomenon are still missing. Here, we show why the data used to describe massive catastrophic rearrangements are incompatible with a model comprising a single event only and propose a molecular mechanism in which a combination of known cellular pathways accounts for chromothripsis. Instead of a single traumatic event, the protection of undamaged chromosomes by telomeres can limit repetitive breakage-fusion-bridge events to a single chromosome arm. Ultimately, common properties of chromosomal instability, such as aneuploidy and centromere fission, might establish the complex genetic pattern observed in this genomic state.