RNA-seq profiling reveals different pathways between remodeled vessels and myocardium in hypertrophic cardiomyopathy.

OBJECTIVE: Coronary microvascular dysfunction (CMD) is a key pathophysiological feature of hypertrophic cardiomyopathy (HCM), contributing to myocardial ischemia and representing a critical determinant of patients' adverse outcome. The molecular mechanisms underlying the morphological and functional changes of CMD are still unknown. Aim of this study was to obtain insights on the molecular pathways associated with microvessel remodeling in HCM. METHODS: Interventricular septum myectomies from patients with obstructive HCM (n = 20) and donors' hearts (CTRL, discarded for technical reasons, n = 7) were collected. Remodeled intramyocardial arterioles and cardiomyocytes were microdissected by laser capture and next-generation sequencing was used to delineate the transcriptome profile. RESULTS: We identified 720 exclusive differentially expressed genes (DEGs) in cardiomyocytes and 1315 exclusive DEGs in remodeled arterioles of HCM. Performing gene ontology and pathway enrichment analyses, we identified selectively altered pathways between remodeled arterioles and cardiomyocytes in HCM patients and controls. CONCLUSIONS: We demonstrate the existence of distinctive pathways between remodeled arterioles and cardiomyocytes in HCM patients and controls at the transcriptome level.

Relative Information Gain: Shannon entropy-based measure of the relative structural conservation in RNA alignments.

Structural characterization of RNAs is a dynamic field, offering many modelling possibilities. RNA secondary structure models are usually characterized by an encoding that depicts structural information of the molecule through string representations or graphs. In this work, we provide a generalization of the BEAR encoding (a context-aware structural encoding we previously developed) by expanding the set of alignments used for the construction of substitution matrices and then applying it to secondary structure encodings ranging from fine-grained to more coarse-grained representations. We also introduce a re-interpretation of the Shannon Information applied on RNA alignments, proposing a new scoring metric, the Relative Information Gain (RIG). The RIG score is available for any position in an alignment, showing how different levels of detail encoded in the RNA representation can contribute differently to convey structural information. The approaches presented in this study can be used alongside state-of-the-art tools to synergistically gain insights into the structural elements that RNAs and RNA families are composed of. This additional information could potentially contribute to their improvement or increase the degree of confidence in the secondary structure of families and any set of aligned RNAs.

Over-Expression of UV-Damage DNA Repair Genes and Ribonucleic Acid Persistence Contribute to the Resilience of Dried Biofilms of the Desert Cyanobacetrium Chroococcidiopsis Exposed to Mars-Like UV Flux and Long-Term Desiccation.

The survival limits of the desert cyanobacterium Chroococcidiopsis were challenged by rewetting dried biofilms and dried biofilms exposed to 1.5 × 103 kJ/m2 of a Mars-like UV, after 7 years of air-dried storage. PCR-stop assays revealed the presence of DNA lesions in dried biofilms and an increased accumulation in dried-UV-irradiated biofilms. Different types and/or amounts of DNA lesions were highlighted by a different expression of uvrA, uvrB, uvrC, phrA, and uvsE genes in dried-rewetted biofilms and dried-UV-irradiated-rewetted biofilms, after rehydration for 30 and 60 min. The up-regulation in dried-rewetted biofilms of uvsE gene encoding an UV damage endonuclease, suggested that UV-damage DNA repair contributed to the repair of desiccation-induced damage. While the phrA gene encoding a photolyase was up-regulated only in dried-UV-irradiated-rewetted biofilms. Nucleotide excision repair genes were over-expressed in dried-rewetted biofilms and dried-UV-irradiated-rewetted biofilms, with uvrC gene showing the highest increase in dried-UV-irradiated-rewetted biofilms. Dried biofilms preserved intact mRNAs (at least of the investigated genes) and 16S ribosomal RNA that the persistence of the ribosome machinery and mRNAs might have played a key role in the early phase recovery. Results have implications for the search of extra-terrestrial life by contributing to the definition of habitability of astrobiologically relevant targets such as Mars or planets orbiting around other stars.

Modeling cancer drug response through drug-specific informative genes.

Recent advances in pharmacogenomics have generated a wealth of data of different types whose analysis have helped in the identification of signatures of different cellular sensitivity/resistance responses to hundreds of chemical compounds. Among the different data types, gene expression has proven to be the more successful for the inference of drug response in cancer cell lines. Although effective, the whole transcriptome can introduce noise in the predictive models, since specific mechanisms are required for different drugs and these realistically involve only part of the proteins encoded in the genome. We analyzed the pharmacogenomics data of 961 cell lines tested with 265 anti-cancer drugs and developed different machine learning approaches for dissecting the genome systematically and predict drug responses using both drug-unspecific and drug-specific genes. These methodologies reach better response predictions for the vast majority of the screened drugs using tens to few hundreds genes specific to each drug instead of the whole genome, thus allowing a better understanding and interpretation of drug-specific response mechanisms which are not necessarily restricted to the drug known targets.

Discovering sequence and structure landscapes in RNA interaction motifs.

RNA molecules are able to bind proteins, DNA and other small or long RNAs using information at primary, secondary or tertiary structure level. Recent techniques that use cross-linking and immunoprecipitation of RNAs can detect these interactions and, if followed by high-throughput sequencing, molecules can be analysed to find recurrent elements shared by interactors, such as sequence and/or structure motifs. Many tools are able to find sequence motifs from lists of target RNAs, while others focus on structure using different approaches to find specific interaction elements. In this work, we make a systematic analysis of RBP-RNA and RNA-RNA datasets to better characterize the interaction landscape with information about multi-motifs on the same RNAs. To achieve this goal, we updated our BEAM algorithm to combine both sequence and structure information to create pairs of patterns that model motifs of interaction. This algorithm was applied to several RNA binding proteins and ncRNAs interactors, confirming already known motifs and discovering new ones. This landscape analysis on interaction variability reflects the diversity of target recognition and underlines that often both primary and secondary structure are involved in molecular recognition.

Axitinib exposure triggers endothelial cells senescence through ROS accumulation and ATM activation.

Inhibitors of Vascular Endothelial Growth Factor target both tumor vasculature and cancer cells that have hijacked VEGF Receptors (VEGFRs) signaling for tumor growth-promoting activities. It is important to get precise insight in the specificity of cell responses to these antiangiogenic drugs to maximize their efficiency and minimize off-target systemic toxicity. Here we report that Axitinib, an inhibitor of VEGFRs currently in use as a second line treatment for advanced renal cell carcinoma, promotes senescence of human endothelial cells in vitro. A one-hour pulse of Axitinib is sufficient for triggering cell senescence. Mechanistically, this requires oxidative stress-dependent activation of the Ataxia Telangiectasia Mutated (ATM) kinase. Axitinib-mediated senescence promoting action is prevented by short-term treatment with antioxidants or ATM inhibitors, which conversely fail to prevent senescence induced by the DNA-damaging drug doxorubicin. Coherently, induction of oxidative stress-related genes distinguishes the response of endothelial cells to Axitinib from that to doxorubicin. Importantly, an Axitinib pulse causes cell senescence in glioblastoma cells. However, neither antioxidants nor ATM inhibitors can reverse this phenotype. Thus, antioxidants may selectively protect endothelial cells from Axitinib by decreasing systemic toxicity and maintaining a functional vascularization necessary for efficient delivery of chemotherapeutic drugs within the tumor mass.

The rs17084733 variant in the KIT 3' UTR disrupts a miR-221/222 binding site in gastrointestinal stromal tumour: a sponge-like mechanism conferring disease susceptibility.

Several miRNAs are dysregulated in gastrointestinal stromal tumours (GIST), and miR-221/222 appear to have a prominent role in GIST biology. Therefore, we investigated the role of DNA variants located in miR-221/222 precursor sequences and their target KIT 3'UTR. Ninety-five polymorphisms were analysed in 115 GIST cases and 88 healthy controls. KIT 3'UTR rs17084733 and pri-miR-222 rs75246947 were found significantly associated with GIST susceptibility. Specifically, KIT rs17084733 A allele was more common in GIST, particularly in KIT wild-type (WT) patients (Padj = 0.017). rs17084733 variant is located within one of the three miR-221/222 binding sites in the KIT 3'UTR, resulting in a mismatch in this seed region. Conversely, KIT mRNA levels were lower in patients carrying the variant allele, except for KIT mutant GIST. Luciferase assay data in GIST cells, generated using a construct containing all the three miR-221/222 binding sites, are consistent with KIT mRNA levels in GIST patients. Reporter assay data, generated using a construct containing only the site encompassing rs17084733, confirmed that this is a functional variant disrupting the miR-221/222 binding site. In conclusion, this is the first study investigating the role of SNPs on miR-221/222 precursor sequences and their binding region on KIT 3'UTR in GIST. We identified the KIT variant rs17084733 as a possible novel genetic biomarker for risk of developing KIT-WT GIST. Moreover, our findings suggest the role of one of the three miR-221/222 binding sites on KIT 3'UTR as endogenous sponge, soaking up and subtracting miR-221/222 to the other two sites characterized by a higher affinity.

BEAM web server: a tool for structural RNA motif discovery.

Motivation: RNA structural motif finding is a relevant problem that becomes computationally hard when working on high-throughput data (e.g. eCLIP, PAR-CLIP), often represented by thousands of RNA molecules. Currently, the BEAM server is the only web tool capable to handle tens of thousands of RNA in input with a motif discovery procedure that is only limited by the current secondary structure prediction accuracies. Results: The recently developed method BEAM (BEAr Motifs finder) can analyze tens of thousands of RNA molecules and identify RNA secondary structure motifs associated to a measure of their statistical significance. BEAM is extremely fast thanks to the BEAR encoding that transforms each RNA secondary structure in a string of characters. BEAM also exploits the evolutionary knowledge contained in a substitution matrix of secondary structure elements, extracted from the RFAM database of families of homologous RNAs. The BEAM web server has been designed to streamline data pre-processing by automatically handling folding and encoding of RNA sequences, giving users a choice for the preferred folding program. The server provides an intuitive and informative results page with the list of secondary structure motifs identified, the logo of each motif, its significance, graphic representation and information about its position in the RNA molecules sharing it. Availability and implementation: The web server is freely available at http://beam.uniroma2.it/ and it is implemented in NodeJS and Python with all major browsers supported. Contact: marco.pietrosanto@uniroma2.it. Supplementary information: Supplementary data are available at Bioinformatics online.

Genome assembly and transcriptome resource for river buffalo, Bubalus bubalis (2n = 50).

Water buffalo is a globally important species for agriculture and local economies. A de novo assembled, well-annotated reference sequence for the water buffalo is an important prerequisite for studying the biology of this species, and is necessary to manage genetic diversity and to use modern breeding and genomic selection techniques. However, no such genome assembly has been previously reported. There are 2 species of domestic water buffalo, the river (2 n = 50) and the swamp (2 n = 48) buffalo. Here we describe a draft quality reference sequence for the river buffalo created from Illumina GA and Roche 454 short read sequences using the MaSuRCA assembler. The assembled sequence is 2.83 Gb, consisting of 366 983 scaffolds with a scaffold N50 of 1.41 Mb and contig N50 of 21 398 bp. Annotation of the genome was supported by transcriptome data from 30 tissues and identified 21 711 predicted protein coding genes. Searches for complete mammalian BUSCO gene groups found 98.6% of curated single copy orthologs present among predicted genes, which suggests a high level of completeness of the genome. The annotated sequence is available from NCBI at accession GCA_000471725.1.

Caspase-8 contributes to angiogenesis and chemotherapy resistance in glioblastoma.

Caspase-8 is a key player in extrinsic apoptosis and its activity is often downregulated in cancer. However, human Caspase-8 expression is retained in some tumors, including glioblastoma (GBM), suggesting that it may support cancer growth in these contexts. GBM, the most aggressive of the gliomas, is characterized by extensive angiogenesis and by an inflammatory microenvironment that support its development and resistance to therapies. We have recently shown that Caspase-8 sustains neoplastic transformation in vitro in human GBM cell lines. Here, we demonstrate that Caspase-8, through activation of NF-kB, enhances the expression and secretion of VEGF, IL-6, IL-8, IL-1beta and MCP-1, leading to neovascularization and increased resistance to Temozolomide. Importantly, the bioinformatics analysis of microarray gene expression data derived from a set of high-grade human gliomas, shows that high Caspase-8 expression levels correlate with a worse prognosis.

RNA-Sequencing for profiling goat milk transcriptome in colostrum and mature milk.

BACKGROUND: In this work we aimed at sequencing and assembling the goat milk transcriptome corresponding at colostrum and 120 days of lactation. To reconstruct transcripts we used both the genome as reference, and a de novo assembly approach. Additionally, we aimed at identifying the differentially expressed genes (DEGs) between the two lactation stages and at analyzing the expression of genes involved in oligosaccharides metabolism. RESULTS: A total of 44,635 different transcripts, organized in 33,757 tentative genes, were obtained using the goat genome as reference. A significant sequence similarity match was found for 40,353 transcripts (90%) against the NCBI NT and for 35,701 (80%) against the NR databases. 68% and 69% of the de novo assembled transcripts, in colostrum and 120 days of lactation samples respectively, have a significant match with the merged transcriptome obtained using Cufflinks/Cuffmerge. CSN2, PAEP, CSN1S2, CSN3, LALBA, TPT1, FTH1, M-SAA3, SPP1, GLYCAM1, EEF1A1, CTSD, FASN, RPS29, CSN1S1, KRT19 and CHEK1 were found between the top fifteen highly expressed genes. 418 loci were differentially expressed between lactation stages, among which 207 and 122 were significantly up- and down-regulated in colostrum, respectively. Functional annotation and pathway enrichment analysis showed that in goat colostrum somatic cells predominate biological processes involved in glycolysis, carbohydrate metabolism, defense response, cytokine activity, regulation of cell proliferation and cell death, vasculature development, while in mature milk, biological process associated with positive regulation of lymphocyte activation and anatomical structure morphogenesis are enriched. The analysis of 144 different oligosaccharide metabolism-related genes showed that most of these (64%) were more expressed in colostrum than in mature milk, with eight expressed at very high levels (SLCA3, GMSD, NME2, SLC2A1, B4GALT1, B3GNT2, NANS, HEXB). CONCLUSIONS: To our knowledge, this is the first study comparing goat transcriptome of two lactation stages: colostrum and 120 days. Our findings suggest putative differences of expression between stages and can be envisioned as a base for further research in the topic. Moreover because a higher expression of genes involved in immune defense response, carbohydrate metabolism and related to oligosaccharide metabolism was identified in colostrum we here corroborate the potential of goat milk as a natural source of lactose-derived oligosaccharides and for the development of functional foods.

A novel method for the identification of conserved structural patterns in RNA: From small scale to high-throughput applications.

Functional RNA regions are often related to recurrent secondary structure patterns (or motifs), which can exert their role in several different ways, particularly in dictating the interaction with RNA-binding proteins, and acting in the regulation of a large number of cellular processes. Among the available motif-finding tools, the majority focuses on sequence patterns, sometimes including secondary structure as additional constraints to improve their performance. Nonetheless, secondary structures motifs may be concurrent to their sequence counterparts or even encode a stronger functional signal. Current methods for searching structural motifs generally require long pipelines and/or high computational efforts or previously aligned sequences. Here, we present BEAM (BEAr Motif finder), a novel method for structural motif discovery from a set of unaligned RNAs, taking advantage of a recently developed encoding for RNA secondary structure named BEAR (Brand nEw Alphabet for RNAs) and of evolutionary substitution rates of secondary structure elements. Tested in a varied set of scenarios, from small- to large-scale, BEAM is successful in retrieving structural motifs even in highly noisy data sets, such as those that can arise in CLIP-Seq or other high-throughput experiments.

c-MYC inhibition impairs hypoxia response in glioblastoma multiforme.

The c-MYC oncoprotein is a DNA binding transcription factor that enhances the expression of many active genes. c-MYC transcriptional signatures vary according to the transcriptional program defined in each cell type during differentiation. Little is known on the involvement of c-MYC in regulation of gene expression programs that are induced by extracellular cues such as a changing microenvironment. Here we demonstrate that inhibition of c-MYC in glioblastoma multiforme cells blunts hypoxia-dependent glycolytic reprogramming and mitochondria fragmentation in hypoxia. This happens because c-MYC inhibition alters the cell transcriptional response to hypoxia and finely tunes the expression of a subset of Hypoxia Inducible Factor 1-regulated genes. We also show that genes whose expression in hypoxia is affected by c-MYC inhibition are able to distinguish the Proneural subtype of glioblastoma multiforme, thus potentially providing a molecular signature for this class of tumors that are the least tractable among glioblastomas.

Genome-wide methylation analysis demonstrates that 5-aza-2-deoxycytidine treatment does not cause random DNA demethylation in fragile X syndrome cells.

BACKGROUND: Fragile X syndrome (FXS) is caused by CGG expansion over 200 repeats at the 5' UTR of the FMR1 gene and subsequent DNA methylation of both the expanded sequence and the CpGs of the promoter region. This epigenetic change causes transcriptional silencing of the gene. We have previously demonstrated that 5-aza-2-deoxycytidine (5-azadC) treatment of FXS lymphoblastoid cell lines reactivates the FMR1 gene, concomitant with CpG sites demethylation, increased acetylation of histones H3 and H4 and methylation of lysine 4 on histone 3. RESULTS: In order to check the specificity of the 5-azadC-induced DNA demethylation, now we performed bisulphite sequencing of the entire methylation boundary upstream the FMR1 promoter region, which is preserved in control wild-type cells. We did not observe any modification of the methylation boundary after treatment. Furthermore, methylation analysis by MS-MLPA of PWS/AS and BWS/SRS loci demonstrated that 5-azadC treatment has no demethylating effect on these regions. Genome-wide methylation analysis through Infinium 450K (Illumina) showed no significant enrichment of specific GO terms in differentially methylated regions after 5-azadC treatment. We also observed that reactivation of FMR1 transcription lasts up to a month after a 7-day treatment and that maximum levels of transcription are reached at 10-15 days after last administration of 5-azadC. CONCLUSIONS: Taken together, these data demonstrate that the demethylating effect of 5-azadC on genomic DNA is not random, but rather restricted to specific regions, if not exclusively to the FMR1 promoter. Moreover, we showed that 5-azadC has a long-lasting reactivating effect on the mutant FMR1 gene.

Next Generation Sequencing and Linkage Analysis for the Molecular Diagnosis of a Novel Overlapping Syndrome Characterized by Hypertrophic Cardiomyopathy and Typical Electrical Instability of Brugada Syndrome.

BACKGROUND: Familial hypertrophic cardiomyopathy (HCM) is an autosomal dominant inherited disorder; mutations in at least 20 genes have been associated. Brugada syndrome (BrS) is an autosomal dominant inherited disorder caused by mutations mainly in theSCN5Agene. A new clinical entity that consists of HCM, typical electrical instability of BrS and sudden death (SD), is described. METHODS AND RESULTS: The family was constituted by 7 members, 4 of who presented clinical features of HCM and electrical instability of BrS. The clinical presentation of proband was ventricular fibrillation. All members were clinically evaluated by physical examination, 12-lead electrocardiography, 2-dimensional echocardiography, stress test, electrocardiogram Holter, flecainide test, and electrophysiological study. An integrated linkage analysis and next generation sequencing (NGS) approach was used to identify the causative mutation. Linkage with the α-tropomyosin (TPM1) gene on chromosome 15q22 was identified. The NGS study identified a missense mutation within theTPM1gene (c.574G>A; p.E192K), exactly located in a binding domain with polycystin-2 protein. No other pathogenic mutations were identified. CONCLUSIONS: This is the first report of an association between HCM and BrS, and the first to use a combined approach of linkage and NGS to identify a causative mutation in SD. The present study expands the clinical spectrum of disorders associated with theTPM1gene and may be useful to report novel mechanisms of electrical instability in HCM and BrS.

The human rs1050286 polymorphism alters LOX-1 expression through modifying miR-24 binding.

The up-regulation of lectin-like oxidized low-density lipoprotein receptor-1 (LOX-1), encoded by the OLR1 gene, plays a fundamental role in the pathogenesis of atherosclerosis. Moreover, OLR1 polymorphisms were associated with increased susceptibility to acute myocardial infarction (AMI) and coronary artery diseases (CAD). In these pathologies, the identification of therapeutic approaches that can inhibit or reduce LOX-1 overexpression is crucial. Predictive analysis showed a putative hsa-miR-24 binding site in the 3'UTR of OLR1, 'naturally' mutated by the presence of the rs1050286 single nucleotide polymorphism (SNP). Luciferase assays revealed that miR-24 targets OLR1 3'UTR-G, but not 3'UTR-A (P < 0.0005). The functional relevance of miR-24 in regulating the expression of OLR1 was established by overexpressing miR-24 in human cell lines heterozygous (A/G, HeLa) and homozygous (A/A, HepG2) for rs1050286 SNP. Accordingly, HeLa (A/G), but not HepG2 (A/A), showed a significant down-regulation of OLR1 both at RNA and protein level. Our results indicate that rs1050286 SNP significantly affects miR-24 binding affinity to the 3'UTR of OLR1, causing a more efficient post-transcriptional gene repression in the presence of the G allele. On this basis, we considered that OLR1 rs1050286 SNP may contribute to modify OLR1 susceptibility to AMI and CAD, so ORL1 SNPs screening could help to stratify patients risk.

Revealing protein-lncRNA interaction.

Long non-coding RNAs (lncRNAs) are associated to a plethora of cellular functions, most of which require the interaction with one or more RNA-binding proteins (RBPs); similarly, RBPs are often able to bind a large number of different RNAs. The currently available knowledge is already drawing an intricate network of interactions, whose deregulation is frequently associated to pathological states. Several different techniques were developed in the past years to obtain protein-RNA binding data in a high-throughput fashion. In parallel, in silico inference methods were developed for the accurate computational prediction of the interaction of RBP-lncRNA pairs. The field is growing rapidly, and it is foreseeable that in the near future, the protein-lncRNA interaction network will rise, offering essential clues for a better understanding of lncRNA cellular mechanisms and their disease-associated perturbations.

Web-Beagle: a web server for the alignment of RNA secondary structures.

Web-Beagle (http://beagle.bio.uniroma2.it) is a web server for the pairwise global or local alignment of RNA secondary structures. The server exploits a new encoding for RNA secondary structure and a substitution matrix of RNA structural elements to perform RNA structural alignments. The web server allows the user to compute up to 10 000 alignments in a single run, taking as input sets of RNA sequences and structures or primary sequences alone. In the latter case, the server computes the secondary structure prediction for the RNAs on-the-fly using RNAfold (free energy minimization). The user can also compare a set of input RNAs to one of five pre-compiled RNA datasets including lncRNAs and 3' UTRs. All types of comparison produce in output the pairwise alignments along with structural similarity and statistical significance measures for each resulting alignment. A graphical color-coded representation of the alignments allows the user to easily identify structural similarities between RNAs. Web-Beagle can be used for finding structurally related regions in two or more RNAs, for the identification of homologous regions or for functional annotation. Benchmark tests show that Web-Beagle has lower computational complexity, running time and better performances than other available methods.

A simple protocol for the inference of RNA global pairwise alignments.

RNA alignment is an important step in the annotation and characterization of unknown RNAs, and several methods have been developed to meet the need of fast and accurate alignments. Being the performances of the aligning methods affected by the input RNA features, finding the most suitable method is not trivial. Indeed, no available method clearly outperforms the others. Here we present a simple workflow to help choosing the more suitable method for RNA pairwise alignment. We tested the performances of six algorithms, based on different approaches, on datasets created by merging publicly available datasets of known or curated RNA secondary structure annotations with datasets of curated RNA alignments. Then, we simulated the frequent case where the secondary structure is unknown by using the same alignment datasets but ignoring the known structure and instead predicting it. In conclusion, the proposed workflow for pairwise RNA alignment depends on the input RNA primary sequence identity and the availability of reliable secondary structures.

Exploiting holistic approaches to model specificity in protein phosphorylation.

Phosphate plays a chemically unique role in shaping cellular signaling of all current living systems, especially eukaryotes. Protein phosphorylation has been studied at several levels, from the near-site context, both in sequence and structure, to the crowded cellular environment, and ultimately to the systems-level perspective. Despite the tremendous advances in mass spectrometry and efforts dedicated to the development of ad hoc highly sophisticated methods, phosphorylation site inference and associated kinase identification are still unresolved problems in kinome biology. The sequence and structure of the substrate near-site context are not sufficient alone to model the in vivo phosphorylation rules, and they should be integrated with orthogonal information in all possible applications. Here we provide an overview of the different contexts that contribute to protein phosphorylation, discussing their potential impact in phosphorylation site annotation and in predicting kinase-substrate specificity.