Comprehensive genomic profiling on metastatic Melanoma: results from a network screening from 7 Italian Cancer Centres.

BACKGROUND: The current therapeutic algorithm for Advanced Stage Melanoma comprises of alternating lines of Targeted and Immuno-therapy, mostly via Immune-Checkpoint blockade. While Comprehensive Genomic Profiling of solid tumours has been approved as a companion diagnostic, still no approved predictive biomarkers are available for Melanoma aside from BRAF mutations and the controversial Tumor Mutational Burden. This study presents the results of a Multi-Centre Observational Clinical Trial of Comprehensive Genomic Profiling on Target and Immuno-therapy treated advanced Melanoma. METHODS: 82 samples, collected from 7 Italian Cancer Centres of FFPE-archived Metastatic Melanoma and matched blood were sequenced via a custom-made 184-gene amplicon-based NGS panel. Sequencing and bioinformatics analysis was performed at a central hub. Primary analysis was carried out via the Ion Reporter framework. Secondary analysis and Machine Learning modelling comprising of uni and multivariate, COX/Lasso combination, and Random Forest, was implemented via custom R/Python scripting. RESULTS: The genomics landscape of the ACC-mela cohort is comparable at the somatic level for Single Nucleotide Variants and INDELs aside a few gene targets. All the clinically relevant targets such as BRAF and NRAS have a comparable distribution thus suggesting the value of larger scale sequencing in melanoma. No comparability is reached at the CNV level due to biotechnological biases and cohort numerosity. Tumour Mutational Burden is slightly higher in median for Complete Responders but fails to achieve statistical significance in Kaplan-Meier survival analysis via several thresholding strategies. Mutations on PDGFRB, NOTCH3 and RET were shown to have a positive effect on Immune-checkpoint treatment Overall and Disease-Free Survival, while variants in NOTCH4 were found to be detrimental for both endpoints. CONCLUSIONS: The results presented in this study show the value and the challenge of a genomics-driven network trial. The data can be also a valuable resource as a validation cohort for Immunotherapy and Target therapy genomic biomarker research.

NR2F1 shapes mitochondria in the mouse brain, providing new insights into Bosch-Boonstra-Schaaf optic atrophy syndrome.

The nuclear receptor NR2F1 acts as a strong transcriptional regulator in embryonic and postnatal neural cells. In humans, mutations in the NR2F1 gene cause Bosch-Boonstra-Schaaf optic atrophy syndrome (BBSOAS), a rare neurodevelopmental disorder characterized by multiple clinical features including vision impairment, intellectual disability and autistic traits. In this study, we identified, by genome-wide and in silico analyses, a set of nuclear-encoded mitochondrial genes as potential genomic targets under direct NR2F1 transcriptional control in neurons. By combining mouse genetic, neuroanatomical and imaging approaches, we demonstrated that conditional NR2F1 loss of function within the adult mouse hippocampal neurogenic niche results in a reduced mitochondrial mass associated with mitochondrial fragmentation and downregulation of key mitochondrial proteins in newborn neurons, the genesis, survival and functional integration of which are impaired. Importantly, we also found dysregulation of several nuclear-encoded mitochondrial genes and downregulation of key mitochondrial proteins in the brain of Nr2f1-heterozygous mice, a validated BBSOAS model. Our data point to an active role for NR2F1 in the mitochondrial gene expression regulatory network in neurons and support the involvement of mitochondrial dysfunction in BBSOAS pathogenesis.

3plex enables deep computational investigation of triplex forming lncRNAs.

Long non-coding RNAs (lncRNAs) regulate gene expression through different molecular mechanisms, including DNA binding via the formation of RNA:DNA:DNA triple helices (TPXs). Despite the increasing amount of experimental evidence, TPXs investigation remains challenging. Here we present 3plex, a software able to predict TPX interactions in silico. Given an RNA sequence and a set of DNA sequences, 3plex integrates 1) Hoogsteen pairing rules that describe the biochemical interactions between RNA and DNA nucleotides, 2) RNA secondary structure prediction and 3) determination of the TPX thermal stability derived from a collection of TPX experimental evidences. We systematically collected and uniformly re-analysed published experimental lncRNA binding sites on human and mouse genomes. We used these data to evaluate 3plex performance and showed that its specific features allow a reliable identification of TPX interactions. We compared 3plex with the other available software and obtained comparable or even better accuracy at a fraction of the computation time. Interestingly, by inspecting collected data with 3plex we found that TPXs tend to be shorter and more degenerated than previously expected and that the majority of analysed lncRNAs can directly bind to the genome by TPX formation. Those results suggest that an important fraction of lncRNAs can exert its biological function through this mechanism. The software is available at https://github.com/molinerisLab/3plex.

DNMT3B supports meso-endoderm differentiation from mouse embryonic stem cells.

The correct establishment of DNA methylation patterns during mouse early development is essential for cell fate specification. However, the molecular targets as well as the mechanisms that determine the specificity of the de novo methylation machinery during differentiation are not completely elucidated. Here we show that the DNMT3B-dependent DNA methylation of key developmental regulatory regions at epiblast-like cells (EpiLCs) provides an epigenetic priming that ensures flawless commitment at later stages. Using in vitro stem cell differentiation and loss of function experiments combined with high-throughput genome-wide bisulfite-, bulk-, and single cell RNA-sequencing we dissected the specific role of DNMT3B in cell fate. We identify DNMT3B-dependent regulatory elements on the genome which, in Dnmt3b knockout (3BKO), impair the differentiation into meso-endodermal (ME) progenitors and redirect EpiLCs towards the neuro-ectodermal lineages. Moreover, ectopic expression of DNMT3B in 3BKO re-establishes the DNA methylation of the master regulator Sox2 super-enhancer, downmodulates its expression, and restores the expression of ME markers. Taken together, our data reveal that DNMT3B-dependent methylation at the epiblast stage is essential for the priming of the meso-endodermal lineages and provide functional characterization of the de novo DNMTs during EpiLCs lineage determination.

Targeting tumor-stroma communication by blocking endothelin-1 receptors sensitizes high-grade serous ovarian cancer to PARP inhibition.

PARP inhibitors (PARPi) have changed the treatment paradigm of high-grade serous ovarian cancer (HG-SOC). However, the impact of this class of inhibitors in HG-SOC patients with a high rate of TP53 mutations is limited, highlighting the need to develop combinatorial therapeutic strategies to improve responses to PARPi. Here, we unveil how the endothelin-1/ET-1 receptor (ET-1/ET-1R) axis, which is overexpressed in human HG-SOC and associated with poor prognosis, instructs HG-SOC/tumor microenvironment (TME) communication via key pro-malignant factors and restricts the DNA damage response induced by the PARPi olaparib. Mechanistically, the ET-1 axis promotes the p53/YAP/hypoxia inducible factor-1α (HIF-1α) transcription hub connecting HG-SOC cells, endothelial cells and activated fibroblasts, hence fueling persistent DNA damage signal escape. The ET-1R antagonist macitentan, which dismantles the ET-1R-mediated p53/YAP/HIF-1α network, interferes with HG-SOC/stroma interactions that blunt PARPi efficacy. Pharmacological ET-1R inhibition by macitentan in orthotopic HG-SOC patient-derived xenografts synergizes with olaparib to suppress metastatic progression, enhancing PARPi survival benefit. These findings reveal ET-1R as a mechanistic determinant in the regulation of HG-SOC/TME crosstalk and DNA damage response, indicating the use of macitentan in combinatorial treatments with PARPi as a promising and emerging therapy.

Circulating miR-1246 and miR-485-3p as Promising Biomarkers of Clinical Response and Outcome in Melanoma Patients Treated with Targeted Therapy.

Despite the significant improvements in advanced melanoma therapy, there is still a pressing need for biomarkers that can predict patient response and prognosis, and therefore support rational treatment decisions. Here, we investigated whether circulating miRNAs could be biomarkers of clinical outcomes in patients treated with targeted therapy. Using next-generation sequencing, we profiled plasma miRNAs at baseline and at progression in patients treated with BRAF inhibitors (BRAFi) or BRAFi + MEKi. Selected miRNAs associated with response to therapy were subjected to validation by real-time quantitative RT-PCR. Receiver Operating Characteristics (ROC), Kaplan-Meier and univariate and multivariate Cox regression analyses were performed on the validated miR-1246 and miR-485-3p baseline levels. The median baseline levels of miR-1246 and miR-485-3p were significantly higher and lower, respectively, in the group of patients not responding to therapy (NRs) as compared with the group of responding patients (Rs). In Rs, a trend toward an increase in miR-1246 and a decrease in miR-485-3p was observed at progression. Baseline miR-1246 level and the miR-1246/miR-485-3p ratio showed a good ability to discriminate between Rs and NRs. Poorer PFS and OS were observed in patients with unfavorable levels of at least one miRNA. In multivariate analysis, a low level of miR-485-3p and a high miR-1246/miR-485-3p ratio remained independent negative prognostic factors for PFS, while a high miR-1246/miR-485-3p ratio was associated with an increased risk of mortality, although statistical significance was not reached. Evaluation of miR-1246 and miR-485-3p baseline plasma levels might help clinicians to identify melanoma patients most likely to be unresponsive to targeted therapy or at higher risk for short-term PFS and mortality, thus improving their management.

Dihydroorotate dehydrogenase inhibition reveals metabolic vulnerability in chronic myeloid leukemia.

The development of different generations of BCR-ABL1 tyrosine kinase inhibitors (TKIs) has led to the high overall survival of chronic myeloid leukemia (CML) patients. However, there are CML patients who show resistance to TKI therapy and are prone to progress to more advanced phases of the disease. So, implementing an alternative approach for targeting TKIs insensitive cells would be of the essence. Dihydroorotate dehydrogenase (DHODH) is an enzyme in the de novo pyrimidine biosynthesis pathway that is located in the inner membrane of mitochondria. Here, we found that CML cells are vulnerable to DHODH inhibition mediated by Meds433, a new and potent DHODH inhibitor recently developed by our group. Meds433 significantly activates the apoptotic pathway and leads to the reduction of amino acids and induction of huge metabolic stress in CML CD34+ cells. Altogether, our study shows that DHODH inhibition is a promising approach for targeting CML stem/progenitor cells and may help more patients discontinue the therapy.

The Emerging Roles of Long Non-Coding RNAs in Intellectual Disability and Related Neurodevelopmental Disorders.

In the human brain, long non-coding RNAs (lncRNAs) are widely expressed in an exquisitely temporally and spatially regulated manner, thus suggesting their contribution to normal brain development and their probable involvement in the molecular pathology of neurodevelopmental disorders (NDD). Bypassing the classic protein-centric conception of disease mechanisms, some studies have been conducted to identify and characterize the putative roles of non-coding sequences in the genetic pathogenesis and diagnosis of complex diseases. However, their involvement in NDD, and more specifically in intellectual disability (ID), is still poorly documented and only a few genomic alterations affecting the lncRNAs function and/or expression have been causally linked to the disease endophenotype. Considering that a significant fraction of patients still lacks a genetic or molecular explanation, we expect that a deeper investigation of the non-coding genome will unravel novel pathogenic mechanisms, opening new translational opportunities. Here, we present evidence of the possible involvement of many lncRNAs in the etiology of different forms of ID and NDD, grouping the candidate disease-genes in the most frequently affected cellular processes in which ID-risk genes were previously collected. We also illustrate new approaches for the identification and prioritization of NDD-risk lncRNAs, together with the current strategies to exploit them in diagnosis.

The acetyltransferase p300 is recruited in trans to multiple enhancer sites by lncSmad7.

The histone acetyltransferase p300 (also known as KAT3B) is a general transcriptional coactivator that introduces the H3K27ac mark on enhancers triggering their activation and gene transcription. Genome-wide screenings demonstrated that a large fraction of long non-coding RNAs (lncRNAs) plays a role in cellular processes and organ development although the underlying molecular mechanisms remain largely unclear (1,2). We found 122 lncRNAs that interacts directly with p300. In depth analysis of one of these, lncSmad7, is required to maintain ESC self-renewal and it interacts to the C-terminal domain of p300. lncSmad7 also contains predicted RNA-DNA Hoogsteen forming base pairing. Combined Chromatin Isolation by RNA precipitation followed by sequencing (ChIRP-seq) together with CRISPR/Cas9 mutagenesis of the target sites demonstrate that lncSmad7 binds and recruits p300 to enhancers in trans, to trigger enhancer acetylation and transcriptional activation of its target genes. Thus, these results unveil a new mechanism by which p300 is recruited to the genome.

Laniakea@ReCaS: exploring the potential of customisable Galaxy on-demand instances as a cloud-based service.

BACKGROUND: Improving the availability and usability of data and analytical tools is a critical precondition for further advancing modern biological and biomedical research. For instance, one of the many ramifications of the COVID-19 global pandemic has been to make even more evident the importance of having bioinformatics tools and data readily actionable by researchers through convenient access points and supported by adequate IT infrastructures. One of the most successful efforts in improving the availability and usability of bioinformatics tools and data is represented by the Galaxy workflow manager and its thriving community. In 2020 we introduced Laniakea, a software platform conceived to streamline the configuration and deployment of "on-demand" Galaxy instances over the cloud. By facilitating the set-up and configuration of Galaxy web servers, Laniakea provides researchers with a powerful and highly customisable platform for executing complex bioinformatics analyses. The system can be accessed through a dedicated and user-friendly web interface that allows the Galaxy web server's initial configuration and deployment. RESULTS: "Laniakea@ReCaS", the first instance of a Laniakea-based service, is managed by ELIXIR-IT and was officially launched in February 2020, after about one year of development and testing that involved several users. Researchers can request access to Laniakea@ReCaS through an open-ended call for use-cases. Ten project proposals have been accepted since then, totalling 18 Galaxy on-demand virtual servers that employ ~ 100 CPUs, ~ 250 GB of RAM and ~ 5 TB of storage and serve several different communities and purposes. Herein, we present eight use cases demonstrating the versatility of the platform. CONCLUSIONS: During this first year of activity, the Laniakea-based service emerged as a flexible platform that facilitated the rapid development of bioinformatics tools, the efficient delivery of training activities, and the provision of public bioinformatics services in different settings, including food safety and clinical research. Laniakea@ReCaS provides a proof of concept of how enabling access to appropriate, reliable IT resources and ready-to-use bioinformatics tools can considerably streamline researchers' work.

Machine learning for RNA sequencing-based intrinsic subtyping of breast cancer.

Stratification of breast cancer (BC) into molecular subtypes by multigene expression assays is of demonstrated clinical utility. In principle, global RNA-sequencing (RNA-seq) should enable reconstructing existing transcriptional classifications of BC samples. Yet, it is not clear whether adaptation to RNA-seq of classifiers originally developed using PCR or microarrays, or reconstruction through machine learning (ML) is preferable. Hence, we focused on robustness and portability of PAM50, a nearest-centroid classifier developed on microarray data to identify five BC "intrinsic subtypes". We found that standard PAM50 is profoundly affected by the composition of the sample cohort used for reference construction, and we propose a strategy, named AWCA, to mitigate this issue, improving classification robustness, with over 90% of concordance, and prognostic ability; we also show that AWCA-based PAM50 can even be applied as single-sample method. Furthermore, we explored five supervised learners to build robust, single-sample intrinsic subtype callers via RNA-seq. From our ML-based survey, regularized multiclass logistic regression (mLR) displayed the best performance, further increased by ad-hoc gene selection on the global transcriptome. On external test sets, mLR classifications reached 90% concordance with PAM50-based calls, without need of reference sample; mLR proven robustness and prognostic ability make it an equally valuable single-sample method to strengthen BC subtyping.

B lymphocytes directly contribute to tissue fibrosis in patients with IgG4-related disease.

BACKGROUND: IgG4-related disease (IgG4-RD) is a fibroinflammatory condition marked by rapid clinical improvement after selective depletion of B lymphocytes with rituximab. This feature suggests that B cells might participate in fibrogenesis and wound healing. OBJECTIVE: In the present work we aimed to demonstrate that B lymphocytes contribute directly to tissue fibrosis in patients with IgG4-RD. METHODS: Total circulating CD19+ B lymphocytes, naive B cells, memory B cells, or plasmablasts from patients with IgG4-RD were cultivated with human fibroblasts. Profibrotic soluble factors and collagen production in cocultures were assessed by using ELISAs and Luminex assays. RNA sequencing and quantitative RT-PCR were used to assess fibroblast activation in the presence of B cells, as well as induction of profibrotic pathways in B-cell subsets. Relevant profibrotic and inflammatory molecules were confirmed in vitro by using functional experiments and on IgG4-RD tissue sections by using multicolor immunofluorescence studies. RESULTS: B cells from patients with IgG4-RD (1) produced the profibrotic molecule platelet-derived growth factor B and stimulated collagen production by fibroblasts; (2) expressed enzymes implicated in extracellular matrix remodeling, such as lysyl oxidase homolog 2; (3) produced the chemotactic factors CCL4, CCL5, and CCL11; and (4) induced production of these same chemokines by activated fibroblasts. Plasmablasts expressed sets of genes implicated in fibroblast activation and proliferation and therefore represent cells with intrinsic profibrotic properties. CONCLUSION: We have demonstrated that B cells contribute directly to tissue fibrosis in patients with IgG4-RD. These unanticipated profibrotic properties of B lymphocytes, particularly plasmablasts, might be relevant for fibrogenesis in patients with other fibroinflammatory disorders and for wound-healing processes in physiologic conditions.

Expanding the spectrum of genes responsible for hereditary motor neuropathies.

BACKGROUND: Inherited peripheral neuropathies (IPNs) represent a broad group of genetically and clinically heterogeneous disorders, including axonal Charcot-Marie-Tooth type 2 (CMT2) and hereditary motor neuropathy (HMN). Approximately 60%-70% of cases with HMN/CMT2 still remain without a genetic diagnosis. Interestingly, mutations in HMN/CMT2 genes may also be responsible for motor neuron disorders or other neuromuscular diseases, suggesting a broad phenotypic spectrum of clinically and genetically related conditions. Thus, it is of paramount importance to identify novel causative variants in HMN/CMT2 patients to better predict clinical outcome and progression. METHODS: We designed a collaborative study for the identification of variants responsible for HMN/CMT2. We collected 15 HMN/CMT2 families with evidence for autosomal recessive inheritance, who had tested negative for mutations in 94 known IPN genes, who underwent whole-exome sequencing (WES) analyses. Candidate genes identified by WES were sequenced in an additional cohort of 167 familial or sporadic HMN/CMT2 patients using next-generation sequencing (NGS) panel analysis. RESULTS: Bioinformatic analyses led to the identification of novel or very rare variants in genes, which have not been previously associated with HMN/CMT2 (ARHGEF28, KBTBD13, AGRN and GNE); in genes previously associated with HMN/CMT2 but in combination with different clinical phenotypes (VRK1 and PNKP), and in the SIGMAR1 gene, which has been linked to HMN/CMT2 in only a few cases. These findings were further validated by Sanger sequencing, segregation analyses and functional studies. CONCLUSIONS: These results demonstrate the broad spectrum of clinical phenotypes that can be associated with a specific disease gene, as well as the complexity of the pathogenesis of neuromuscular disorders.

The Suv420h histone methyltransferases regulate PPAR-γ and energy expenditure in response to environmental stimuli.

Obesity and its associated metabolic abnormalities have become a global emergency with considerable morbidity and mortality. Epidemiologic and animal model data suggest an epigenetic contribution to obesity. Nevertheless, the cellular and molecular mechanisms through which epigenetics contributes to the development of obesity remain to be elucidated. Suv420h1 and Suv420h2 are histone methyltransferases responsible for chromatin compaction and gene repression. Through in vivo, ex vivo, and in vitro studies, we found that Suv420h1 and Suv420h2 respond to environmental stimuli and regulate metabolism by down-regulating peroxisome proliferator-activated receptor gamma (PPAR-γ), a master transcriptional regulator of lipid storage and glucose metabolism. Accordingly, mice lacking Suv420h proteins activate PPAR-γ target genes in brown adipose tissue to increase mitochondria respiration, improve glucose tolerance, and reduce adipose tissue to fight obesity. We conclude that Suv420h proteins are key epigenetic regulators of PPAR-γ and the pathways controlling metabolism and weight balance in response to environmental stimuli.

SP1 and STAT3 Functionally Synergize to Induce the RhoU Small GTPase and a Subclass of Non-canonical WNT Responsive Genes Correlating with Poor Prognosis in Breast Cancer.

Breast cancer is a heterogeneous disease whose clinical management is very challenging. Although specific molecular features characterize breast cancer subtypes with different prognosis, the identification of specific markers predicting disease outcome within the single subtypes still lags behind. Both the non-canonical Wingless-type MMTV Integration site (WNT) and the Signal Transducer and Activator of Transcription (STAT)3 pathways are often constitutively activated in breast tumors, and both can induce the small GTPase Ras Homolog Family Member U RhoU. Here we show that RhoU transcription can be triggered by both canonical and non-canonical WNT ligands via the activation of c-JUN N-terminal kinase (JNK) and the recruitment of the Specificity Protein 1 (SP1) transcription factor to the RhoU promoter, identifying for the first time SP1 as a JNK-dependent mediator of WNT signaling. RhoU down-regulation by silencing or treatment with JNK, SP1 or STAT3 inhibitors leads to impaired migration and invasion in basal-like MDA-MB-231 and BT-549 cells, suggesting that STAT3 and SP1 can cooperate to induce high RhoU expression and enhance breast cancer cells migration. Moreover, in vivo concomitant binding of STAT3 and SP1 defines a subclass of genes belonging to the non-canonical WNT and the Interleukin (IL)-6/STAT3 pathways and contributing to breast cancer aggressiveness, suggesting the relevance of developing novel targeted therapies combining inhibitors of the STAT3 and WNT pathways or of their downstream mediators.

The Trithorax protein Ash1L promotes myoblast fusion by activating Cdon expression.

Myoblast fusion (MF) is required for muscle growth and repair, and its alteration contributes to muscle diseases. The mechanisms governing this process are incompletely understood, and no epigenetic regulator has been previously described. Ash1L is an epigenetic activator belonging to the Trithorax group of proteins and is involved in FSHD muscular dystrophy, autism and cancer. Its physiological role in skeletal muscle is unknown. Here we report that Ash1L expression is positively correlated with MF and reduced in Duchenne muscular dystrophy. In vivo, ex vivo and in vitro experiments support a selective and evolutionary conserved requirement for Ash1L in MF. RNA- and ChIP-sequencing indicate that Ash1L is required to counteract Polycomb repressive activity to allow activation of selected myogenesis genes, in particular the key MF gene Cdon. Our results promote Ash1L as an important epigenetic regulator of MF and suggest that its activity could be targeted to improve cell therapy for muscle diseases.

Correction to: ZIKA virus elicits P53 activation and genotoxic stress in human neural progenitors similar to mutations involved in severe forms of genetic microcephaly.

The authors wish to point out that the name of the first author is appearing incorrectly on Pubmed: it should be El Ghouzzi V (and not Ghouzzi VE). In addition, the words "and p53" appear at the end of the title in the original publication ( https://www.nature.com/articles/cddis2016266 ) and in the previous erratum version ( https://www.nature.com/articles/cddis2016446 ). This is not correct.

Evolutionary Rewiring of Human Regulatory Networks by Waves of Genome Expansion.

Genome expansion is believed to be an important driver of the evolution of gene regulation. To investigate the role of a newly arising sequence in rewiring regulatory networks, we estimated the age of each region of the human genome by applying maximum parsimony to genome-wide alignments with 100 vertebrates. We then studied the age distribution of several types of functional regions, with a focus on regulatory elements. The age distribution of regulatory elements reveals the extensive use of newly formed genomic sequence in the evolution of regulatory interactions. Many transcription factors have expanded their repertoire of targets through waves of genomic expansions that can be traced to specific evolutionary times. Repeated elements contributed a major part of such expansion: many classes of such elements are enriched in binding sites of one or a few specific transcription factors, whose binding sites are localized in specific portions of the element and characterized by distinctive motif words. These features suggest that the binding sites were available as soon as the new sequence entered the genome, rather than being created later by accumulation of point mutations. By comparing the age of regulatory regions to the evolutionary shift in expression of nearby genes, we show that rewiring through genome expansion played an important role in shaping human regulatory networks.

A functional strategy to characterize expression Quantitative Trait Loci.

The study of genetic variation has been revolutionized by the advent of high-throughput technologies able to determine the complete genomic sequence of thousands of individuals. Understanding the functional relevance of variants is, however, still a difficult task, especially when focusing on non-coding variants. Most of the variants associated with disease by Genome-Wide Association Studies (GWAS) are indeed non-coding, and presumably exert their effects by altering gene regulation. Expression Quantitative Trait Loci (eQTL) studies represent an important step in understanding the functional relevance of regulatory variants. We propose a new strategy to detect and characterize eQTLs, based on the effect of variants on the Total Binding Affinity (TBA) profiles of regulatory regions. Using a large dataset of coupled genome and expression data, we show that TBA-based inference allows the identification of eQTLs not revealed by traditional methods and helps in their interpretation in terms of altered transcription factor binding.

In silico prediction of lncRNA function using tissue specific and evolutionary conserved expression.

BACKGROUND: In recent years long non coding RNAs (lncRNAs) have been the subject of increasing interest. Thanks to many recent functional studies, the existence of a large class of lncRNAs with potential regulatory functions is now widely accepted. Although an increasing number of lncRNAs is being characterized and shown to be involved in many biological processes, the functions of the vast majority lncRNA genes is still unknown. Therefore computational methods able to take advantage of the increasing amount of publicly available data to predict lncRNA functions could be very useful. RESULTS: Since coding genes are much better annotated than lncRNAs, we attempted to project known functional information regarding proteins onto non coding genes using the guilt by association principle: if a gene shows an expression profile that correlates with those of a set of coding genes involved in a given function, that gene is probably involved in the same function. We computed gene coexpression for 30 human tissues and 9 vertebrates and mined the resulting networks with a methodology inspired by the rank product algorithm used to identify differentially expressed genes. Using different types of reference data we can predict putative new annotations for thousands of lncRNAs and proteins, ranging from cellular localization to relevance for disease and cancer. CONCLUSIONS: New function of coding genes and lncRNA can be profitably predicted using tissue specific coexpression, as well as expression of orthologous genes in different species. The data are available for download and through a user-friendly web interface at www.funcpred.com .