Establishment and Characterization of a Cell Line (S-RMS1) Derived from an Infantile Spindle Cell Rhabdomyosarcoma with SRF-NCOA2 Fusion Transcript.

Background: Spindle cell rhabdomyosarcoma (S-RMS) is a rare tumor that was previously considered as an uncommon variant of embryonal RMS (ERMS) and recently reclassified as a distinct RMS subtype with NCOA2, NCOA1, and VGLL2 fusion genes. In this study, we established a cell line (S-RMS1) derived from a four-month-old boy with infantile spindle cell RMS harboring SRF-NCOA2 gene fusion. Methods: Morphological and molecular characteristics of S-RMS1 were analyzed and compared with two RMS cell lines, RH30 and RD18. Whole genome sequencing of S-RMS1 and clinical exome sequencing of genomic DNA were performed. Results: S-RMS1 showed cells small in size, with a fibroblast-like morphology and positivity for MyoD-1, myogenin, desmin, and smooth muscle actin. The population doubling time was 3.7 days. Whole genome sequencing demonstrated that S-RMS1 retained the same genetic profile of the tumor at diagnosis. A Western blot analysis showed downregulation of AKT-p and YAP-p while RT-qPCR showed upregulation of endoglin and GATA6 as well as downregulation of TGFßR1 and Mef2C transcripts. Conclusion: This is the first report of the establishment of a cell line from an infantile spindle cell RMS with SRF-NCOA2 gene fusion. S-RMS1 should represent a useful tool for the molecular characterization of this rare and almost unknown tumor.

Alzheimer's brains show inter-related changes in RNA and lipid metabolism.

Alzheimer's disease (AD) involves changes in both lipid and RNA metabolism, but it remained unknown if these differences associate with AD's cognition and/or post-mortem neuropathology indices. Here, we report RNA-sequencing evidence of inter-related associations between lipid processing, cognition level, and AD neuropathology. In two unrelated cohorts, we identified pathway-enriched facilitation of lipid processing and alternative splicing genes, including the neuronal-enriched NOVA1 and hnRNPA1. Specifically, this association emerged in temporal lobe tissue samples from donors where postmortem evidence demonstrated AD neuropathology, but who presented normal cognition proximate to death. The observed changes further associated with modified ATP synthesis and mitochondrial transcripts, indicating metabolic relevance; accordingly, mass-spectrometry-derived lipidomic profiles distinguished between individuals with and without cognitive impairment prior to death. In spite of the limited group sizes, tissues from persons with both cognitive impairment and AD pathology showed elevation in several drug-targeted genes of other brain, vascular and autoimmune disorders, accompanied by pathology-related increases in distinct lipid processing transcripts, and in the RNA metabolism genes hnRNPH2, TARDBP, CLP1 and EWSR1. To further detect 3'-polyadenylation variants, we employed multiple cDNA primer pairs. This identified variants that showed limited differences in scope and length between the tested cohorts, yet enabled superior clustering of demented and non-demented AD brains versus controls compared to total mRNA expression values. Our findings indicate inter-related cognition-associated differences in AD's lipid processing, alternative splicing and 3'-polyadenylation, calling for pursuing the underlying psychological and therapeutics implications.

A sketch of known and novel MYCN-associated miRNA networks in neuroblastoma.

Neuroblastoma (NB) originates from neural crest-derived precursors and represents the most common childhood extracranial solid tumour. MicroRNAs (miRNAs), a class of small non-coding RNAs that participate in a wide variety of biological processes by regulating gene expression, appear to play an essential role within the NB context. High-throughput next generation sequencing (NGS) was applied to study the miRNA transcriptome in a cohort of NB tumours with and without MYCN-amplification (MNA and MNnA, respectively) and in dorsal root ganglia (DRG), as a control. Out of the 128 miRNAs differentially expressed in the NB vs. DRG comparison, 47 were expressed at higher levels, while 81 were expressed at lower levels in the NB tumours. We also found that 23 miRNAs were differentially expressed in NB with or without MYCN-amplification, with 17 miRNAs being upregulated and 6 being downregulated in the MNA subtypes. Functional annotation analysis of the target genes of these differentially expressed miRNAs demonstrated that many mRNAs were involved in cancer-related pathways, such as DNA-repair and apoptosis as well as FGFR and EGFR signalling. In particular, we found that miR-628-3p negatively affects MYCN gene expression. Furthermore, we identified a novel miRNA candidate with variable expression in MNA vs. MNnA tumours, whose putative target genes are implicated in the mTOR pathway. The present study provides further insight into the molecular mechanisms that correlate miRNA dysregulation to NB development and progression.

A Review of Research-Grade Human Induced Pluripotent Stem Cells Qualification and Biobanking Processes.

Human induced pluripotent stem cell (hiPSC) biobanks are invaluable resources for basic and clinical research, since they provide a sustainable supply of accessible cell lines that meet high quality and safety standards. hiPSCs are particularly useful for understanding disease mechanisms, creating cell models for drug development, and generating novel clinical therapies. For clinical applications and drug discovery, it is fundamental that the acquired pluripotent cell lines never touch animal-derived products nor xenogeneic reagents (Good Manufacturing Practice-grade); whereas for research grade, it is sufficient to operate under Good Laboratory Practice conditions. However, regardless of the end use, it is important that every step in the whole process, starting from the original cells throughout expansion and manipulation, must be performed and recorded rigorously. Here, we describe our biobanking management system that is applied specifically to human pluripotent stem cells.

Histone Modifications in a Mouse Model of Early Adversities and Panic Disorder: Role for Asic1 and Neurodevelopmental Genes.

Hyperventilation following transient, CO2-induced acidosis is ubiquitous in mammals and heritable. In humans, respiratory and emotional hypersensitivity to CO2 marks separation anxiety and panic disorders, and is enhanced by early-life adversities. Mice exposed to the repeated cross-fostering paradigm (RCF) of interference with maternal environment show heightened separation anxiety and hyperventilation to 6% CO2-enriched air. Gene-environment interactions affect CO2 hypersensitivity in both humans and mice. We therefore hypothesised that epigenetic modifications and increased expression of genes involved in pH-detection could explain these relationships. Medullae oblongata of RCF- and normally-reared female outbred mice were assessed by ChIP-seq for H3Ac, H3K4me3, H3K27me3 histone modifications, and by SAGE for differential gene expression. Integration of multiple experiments by network analysis revealed an active component of 148 genes pointing to the mTOR signalling pathway and nociception. Among these genes, Asic1 showed heightened mRNA expression, coherent with RCF-mice's respiratory hypersensitivity to CO2 and altered nociception. Functional enrichment and mRNA transcript analyses yielded a consistent picture of enhancement for several genes affecting chemoception, neurodevelopment, and emotionality. Particularly, results with Asic1 support recent human findings with panic and CO2 responses, and provide new perspectives on how early adversities and genes interplay to affect key components of panic and related disorders.

BITS 2015: the annual meeting of the Italian Society of Bioinformatics.

This preface introduces the content of the BioMed Central journal Supplements related to the BITS 2015 meeting, held in Milan, Italy, from the 3th to the 5th of June, 2015.

In-depth characterization of breast cancer tumor-promoting cell transcriptome by RNA sequencing and microarrays.

Numerous studies have reported the existence of tumor-promoting cells (TPC) with self-renewal potential and a relevant role in drug resistance. However, pathways and modifications involved in the maintenance of such tumor subpopulations are still only partially understood. Sequencing-based approaches offer the opportunity for a detailed study of TPC including their transcriptome modulation. Using microarrays and RNA sequencing approaches, we compared the transcriptional profiles of parental MCF7 breast cancer cells with MCF7-derived TPC (i.e. MCFS). Data were explored using different bioinformatic approaches, and major findings were experimentally validated. The different analytical pipelines (Lifescope and Cufflinks based) yielded similar although not identical results. RNA sequencing data partially overlapped microarray results and displayed a higher dynamic range, although overall the two approaches concordantly predicted pathway modifications. Several biological functions were altered in TPC, ranging from production of inflammatory cytokines (i.e., IL-8 and MCP-1) to proliferation and response to steroid hormones. More than 300 non-coding RNAs were defined as differentially expressed, and 2,471 potential splicing events were identified. A consensus signature of genes up-regulated in TPC was derived and was found to be significantly associated with insensitivity to fulvestrant in a public breast cancer patient dataset. Overall, we obtained a detailed portrait of the transcriptome of a breast cancer TPC line, highlighted the role of non-coding RNAs and differential splicing, and identified a gene signature with a potential as a context-specific biomarker in patients receiving endocrine treatment.

Whole transcriptome RNA sequencing data from blood leukocytes derived from Parkinson's disease patients prior to and following deep brain stimulation treatment.

Recent evidence demonstrates the power of RNA sequencing (RNA-Seq) for identifying valuable and urgently needed blood biomarkers and advancing both early and accurate detection of neurological diseases, and in particular Parkinson's disease (PD). RNA sequencing technology enables non-biased, high throughput, probe-independent inspection of expression data and high coverage and both quantification of global transcript levels as well as the detection of expressed exons and junctions given a sufficient sequencing depth (coverage). However, the analysis of sequencing data frequently presents a bottleneck. Tools for quantification of alternative splicing from sequenced libraries hardly exist at the present time, and methods that support multiple sequencing platforms are especially lacking. Here, we describe in details a whole RNA-Seq transcriptome dataset produced from PD patient's blood leukocytes. The samples were taken prior to, and following deep brain stimulation (DBS) treatment while being on stimulation and following 1 h of complete electrical stimulation cessation and from healthy control volunteers. We describe in detail the methodology applied for analyzing the RNA-Seq data including differential expression of long noncoding RNAs (lncRNAs). We also provide details of the corresponding analysis of in-depth splice isoform data from junction and exon reads, with the use of the software AltAnalyze. Both the RNA-Seq raw (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE42608) and analyzed data (https://www.synapse.org/#!Synapse:syn2805267) may be found valuable towards detection of novel blood biomarkers for PD.

Preface: BITS2014, the annual meeting of the Italian Society of Bioinformatics.

This Preface introduces the content of the BioMed Central journal Supplements related to BITS2014 meeting, held in Rome, Italy, from the 26th to the 28th of February, 2014.

Impact of different ChIP-Seq protocols on DNA integrity and quality of bioinformatics analysis results.

Different ChIP-Seq protocols may have a significant impact on the final outcome in terms of quality, number and distribution of called peaks. Sample DNA undergoes a long procedure before the final sequencing step, and damaged DNA can result in excessive mismatches in the alignment with reference genome. In this letter, we present the effect of well-defined modifications (timing of formaldehyde crosslink reversal, brand of the sonicator) of standard ChIP-Seq protocol on parallel samples derived from the same cell line correlating the initial DNA quality control metrics to the final bioinformatics analysis results.

Deep Sequencing the microRNA profile in rhabdomyosarcoma reveals down-regulation of miR-378 family members.

BACKGROUND: Rhabdomyosarcoma (RMS) is a highly malignant tumour accounting for nearly half of soft tissue sarcomas in children. MicroRNAs (miRNAs) represent a class of short, non-coding, regulatory RNAs which play a critical role in different cellular processes. Altered miRNA levels have been reported in human cancers, including RMS. METHODS: Using deep sequencing technology, a total of 685 miRNAs were investigated in a group of alveolar RMSs (ARMSs), embryonal RMSs (ERMSs) as well as in normal skeletal muscle (NSM). Q-PCR, MTT, cytofluorimetry, migration assay, western blot and immunofluorescence experiments were carried out to determine the role of miR-378a-3p in cancer cell growth, apoptosis, migration and differentiation. Bioinformatics pipelines were used for miRNA target prediction and clustering analysis. RESULTS: Ninety-seven miRNAs were significantly deregulated in ARMS and ERMS when compared to NSM. MiR-378 family members were dramatically decreased in RMS tumour tissue and cell lines. Interestingly, members of the miR-378 family presented as a possible target the insulin-like growth factor receptor 1 (IGF1R), a key signalling molecule in RMS. MiR-378a-3p over-expression in an RMS-derived cell line suppressed IGF1R expression and affected phosphorylated-Akt protein levels. Ectopic expression of miR-378a-3p caused significant changes in apoptosis, cell migration, cytoskeleton organization as well as a modulation of the muscular markers MyoD1, MyoR, desmin and MyHC. In addition, DNA demethylation by 5-aza-2'-deoxycytidine (5-aza-dC) was able to up-regulate miR-378a-3p levels with a concomitant induction of apoptosis, decrease in cell viability and cell cycle arrest in G2-phase. Cells treated with 5-aza-dC clearly changed their morphology and expressed moderate levels of MyHC. CONCLUSIONS: MiR-378a-3p may function as a tumour suppressor in RMS and the restoration of its expression would be of therapeutic benefit in RMS. Furthermore, the role of epigenetic modifications in RMS deserves further investigations.

Emerging bioinformatics approaches for analysis of NGS-derived coding and non-coding RNAs in neurodegenerative diseases.

Neurodegenerative diseases in general and specifically late-onset Alzheimer's disease (LOAD) involve a genetically complex and largely obscure ensemble of causative and risk factors accompanied by complex feedback responses. The advent of "high-throughput" transcriptome investigation technologies such as microarray and deep sequencing is increasingly being combined with sophisticated statistical and bioinformatics analysis methods complemented by knowledge-based approaches such as Bayesian Networks or network and graph analyses. Together, such "integrative" studies are beginning to identify co-regulated gene networks linked with biological pathways and potentially modulating disease predisposition, outcome, and progression. Specifically, bioinformatics analyses of integrated microarray and genotyping data in cases and controls reveal changes in gene expression of both protein-coding and small and long regulatory RNAs; highlight relevant quantitative transcriptional differences between LOAD and non-demented control brains and demonstrate reconfiguration of functionally meaningful molecular interaction structures in LOAD. These may be measured as changes in connectivity in "hub nodes" of relevant gene networks (Zhang etal., 2013). We illustrate here the open analytical questions in the transcriptome investigation of neurodegenerative disease studies, proposing "ad hoc" strategies for the evaluation of differential gene expression and hints for a simple analysis of the non-coding RNA (ncRNA) part of such datasets. We then survey the emerging role of long ncRNAs (lncRNAs) in the healthy and diseased brain transcriptome and describe the main current methods for computational modeling of gene networks. We propose accessible modular and pathway-oriented methods and guidelines for bioinformatics investigations of whole transcriptome next generation sequencing datasets. We finally present methods and databases for functional interpretations of lncRNAs and propose a simple heuristic approach to visualize and represent physical and functional interactions of the coding and non-coding components of the transcriptome. Integrating in a functional and integrated vision coding and ncRNA analyses is of utmost importance for current and future analyses of neurodegenerative transcriptomes.

Long non-coding RNA and alternative splicing modulations in Parkinson's leukocytes identified by RNA sequencing.

The continuously prolonged human lifespan is accompanied by increase in neurodegenerative diseases incidence, calling for the development of inexpensive blood-based diagnostics. Analyzing blood cell transcripts by RNA-Seq is a robust means to identify novel biomarkers that rapidly becomes a commonplace. However, there is lack of tools to discover novel exons, junctions and splicing events and to precisely and sensitively assess differential splicing through RNA-Seq data analysis and across RNA-Seq platforms. Here, we present a new and comprehensive computational workflow for whole-transcriptome RNA-Seq analysis, using an updated version of the software AltAnalyze, to identify both known and novel high-confidence alternative splicing events, and to integrate them with both protein-domains and microRNA binding annotations. We applied the novel workflow on RNA-Seq data from Parkinson's disease (PD) patients' leukocytes pre- and post- Deep Brain Stimulation (DBS) treatment and compared to healthy controls. Disease-mediated changes included decreased usage of alternative promoters and N-termini, 5'-end variations and mutually-exclusive exons. The PD regulated FUS and HNRNP A/B included prion-like domains regulated regions. We also present here a workflow to identify and analyze long non-coding RNAs (lncRNAs) via RNA-Seq data. We identified reduced lncRNA expression and selective PD-induced changes in 13 of over 6,000 detected leukocyte lncRNAs, four of which were inversely altered post-DBS. These included the U1 spliceosomal lncRNA and RP11-462G22.1, each entailing sequence complementarity to numerous microRNAs. Analysis of RNA-Seq from PD and unaffected controls brains revealed over 7,000 brain-expressed lncRNAs, of which 3,495 were co-expressed in the leukocytes including U1, which showed both leukocyte and brain increases. Furthermore, qRT-PCR validations confirmed these co-increases in PD leukocytes and two brain regions, the amygdala and substantia-nigra, compared to controls. This novel workflow allows deep multi-level inspection of RNA-Seq datasets and provides a comprehensive new resource for understanding disease transcriptome modifications in PD and other neurodegenerative diseases.

Mutations in the FOG2/ZFPM2 gene are associated with anomalies of human testis determination.

In recent years, considerable advances have been made in our understanding of genetics of mammalian gonad development; however, the underlying genetic aetiology in the majority of patients with 46,XY disorders of sex development (DSD) still remains unknown. Based on mouse models, it has been hypothesized that haploinsufficiency of the Friend of GATA 2 (FOG2) gene could lead to 46,XY gonadal dysgenesis on specific inbred genetic backgrounds. Using whole exome sequencing, we identified independent missense mutations in FOG2 in two patients with 46,XY gonadal dysgenesis. One patient carried a non-synonymous heterozygous mutation (p.S402R), while the other patient carried a heterozygous p.R260Q mutation and a homozygous p.M544I mutation. Functional studies indicated that the failure of testis development in these cases could be explained by the impaired ability of the mutant FOG2 proteins to interact with a known regulator of early testis development, GATA4. This is the first example of mutations in the coding sequence of FOG2 associated with 46,XY DSD in human and adds to the list of genes in the human known to be associated with DSD.

Small noncoding RNAs in cells transformed by human T-cell leukemia virus type 1: a role for a tRNA fragment as a primer for reverse transcriptase.

UNLABELLED: The present study employed mass sequencing of small RNA libraries to identify the repertoire of small noncoding RNAs expressed in normal CD4(+) T cells compared to cells transformed with human T-cell leukemia virus type 1 (HTLV-1), the causative agent of adult T-cell leukemia/lymphoma (ATLL). The results revealed distinct patterns of microRNA expression in HTLV-1-infected CD4(+) T-cell lines with respect to their normal counterparts. In addition, a search for virus-encoded microRNAs yielded 2 sequences that originated from the plus strand of the HTLV-1 genome. Several sequences derived from tRNAs were expressed at substantial levels in both uninfected and infected cells. One of the most abundant tRNA fragments (tRF-3019) was derived from the 3' end of tRNA-proline. tRF-3019 exhibited perfect sequence complementarity to the primer binding site of HTLV-1. The results of an in vitro reverse transcriptase assay verified that tRF-3019 was capable of priming HTLV-1 reverse transcriptase. Both tRNA-proline and tRF-3019 were detected in virus particles isolated from HTLV-1-infected cells. These findings suggest that tRF-3019 may play an important role in priming HTLV-1 reverse transcription and could thus represent a novel target to control HTLV-1 infection. IMPORTANCE: Small noncoding RNAs, a growing family of regulatory RNAs that includes microRNAs and tRNA fragments, have recently emerged as key players in many biological processes, including viral infection and cancer. In the present study, we employed mass sequencing to identify the repertoire of small noncoding RNAs in normal T cells compared to T cells transformed with human T-cell leukemia virus type 1 (HTLV-1), a retrovirus that causes adult T-cell leukemia/lymphoma. The results revealed a distinct pattern of microRNA expression in HTLV-1-infected cells and a tRNA fragment (tRF-3019) that was packaged into virions and capable of priming HTLV-1 reverse transcription, a key event in the retroviral life cycle. These findings indicate tRF-3019 could represent a novel target for therapies aimed at controlling HTLV-1 infection.

Genome-wide analysis of histone marks identifying an epigenetic signature of promoters and enhancers underlying cardiac hypertrophy.

Cardiac hypertrophy, initially an adaptive response of the myocardium to stress, can progress to heart failure. The epigenetic signature underlying this phenomenon is poorly understood. Here, we report on the genome-wide distribution of seven histone modifications in adult mouse cardiomyocytes subjected to a prohypertrophy stimulus in vivo. We found a set of promoters with an epigenetic pattern that distinguishes specific functional classes of genes regulated in hypertrophy and identified 9,207 candidate active enhancers whose activity was modulated. We also analyzed the transcriptional network within which these genetic elements act to orchestrate hypertrophy gene expression, finding a role for myocyte enhancer factor (MEF)2C and MEF2A in regulating enhancers. We propose that the epigenetic landscape is a key determinant of gene expression reprogramming in cardiac hypertrophy and provide a basis for understanding the role of chromatin in regulating this phenomenon.

Transient exposure to low levels of insecticide affects metabolic networks of honeybee larvae.

The survival of a species depends on its capacity to adjust to changing environmental conditions, and new stressors. Such new, anthropogenic stressors include the neonicotinoid class of crop-protecting agents, which have been implicated in the population declines of pollinating insects, including honeybees (Apis mellifera). The low-dose effects of these compounds on larval development and physiological responses have remained largely unknown. Over a period of 15 days, we provided syrup tainted with low levels (2 µg/L(-1)) of the neonicotinoid insecticide imidacloprid to beehives located in the field. We measured transcript levels by RNA sequencing and established lipid profiles using liquid chromatography coupled with mass spectrometry from worker-bee larvae of imidacloprid-exposed (IE) and unexposed, control (C) hives. Within a catalogue of 300 differentially expressed transcripts in larvae from IE hives, we detect significant enrichment of genes functioning in lipid-carbohydrate-mitochondrial metabolic networks. Myc-involved transcriptional response to exposure of this neonicotinoid is indicated by overrepresentation of E-box elements in the promoter regions of genes with altered expression. RNA levels for a cluster of genes encoding detoxifying P450 enzymes are elevated, with coordinated downregulation of genes in glycolytic and sugar-metabolising pathways. Expression of the environmentally responsive Hsp90 gene is also reduced, suggesting diminished buffering and stability of the developmental program. The multifaceted, physiological response described here may be of importance to our general understanding of pollinator health. Muscles, for instance, work at high glycolytic rates and flight performance could be impacted should low levels of this evolutionarily novel stressor likewise induce downregulation of energy metabolising genes in adult pollinators.

Characterization of the genome of the dairy Lactobacillus helveticus bacteriophage {Phi}AQ113.

The complete genomic sequence of the dairy Lactobacillus helveticus bacteriophage ΦAQ113 was determined. Phage ΦAQ113 is a Myoviridae bacteriophage with an isometric capsid and a contractile tail. The final assembled consensus sequence revealed a linear, circularly permuted, double-stranded DNA genome with a size of 36,566 bp and a G+C content of 37%. Fifty-six open reading frames (ORFs) were predicted, and a putative function was assigned to approximately 90% of them. The ΦAQ113 genome shows functionally related genes clustered together in a genome structure composed of modules for DNA replication/regulation, DNA packaging, head and tail morphogenesis, cell lysis, and lysogeny. The identification of genes involved in the establishment of lysogeny indicates that it may have originated as a temperate phage, even if it was isolated from natural cheese whey starters as a virulent phage, because it is able to propagate in a sensitive host strain. Additionally, we discovered that the ΦAQ113 phage genome is closely related to Lactobacillus gasseri phage KC5a and Lactobacillus johnsonii phage Lj771 genomes. The phylogenetic similarities between L. helveticus phage ΦAQ113 and two phages that belong to gut species confirm a possible common ancestral origin and support the increasing consideration of L. helveticus as a health-promoting organism.

Estrogen-dependent dynamic profile of eNOS-DNA associations in prostate cancer.

In previous work we have documented the nuclear translocation of endothelial NOS (eNOS) and its participation in combinatorial complexes with Estrogen Receptor Beta (ERß) and Hypoxia Inducible Factors (HIFs) that determine localized chromatin remodeling in response to estrogen (E2) and hypoxia stimuli, resulting in transcriptional regulation of genes associated with adverse prognosis in prostate cancer (PCa). To explore the role of nuclear eNOS in the acquisition of aggressive phenotype in PCa, we performed ChIP-Sequencing on chromatin-associated eNOS from cells from a primary tumor with poor outcome and from metastatic LNCaP cells. We found that: 1. the eNOS-bound regions (peaks) are widely distributed across the genome encompassing multiple transcription factors binding sites, including Estrogen Response Elements. 2. E2 increased the number of peaks, indicating hormone-dependent eNOS re-localization. 3. Peak distribution was similar with/without E2 with ≈ 55% of them in extragenic DNA regions and an intriguing involvement of the 5' domain of several miRs deregulated in PCa. Numerous potentially novel eNOS-targeted genes have been identified suggesting that eNOS participates in the regulation of large gene sets. The parallel finding of downregulation of a cluster of miRs, including miR-34a, in PCa cells associated with poor outcome led us to unveil a molecular link between eNOS and SIRT1, an epigenetic regulator of aging and tumorigenicity, negatively regulated by miR-34a and in turn activating eNOS. E2 potentiates miR-34a downregulation thus enhancing SIRT1 expression, depicting a novel eNOS/SIRT1 interplay fine-tuned by E2-activated ER signaling, and suggesting that eNOS may play an important role in aggressive PCa.

Deep-sequencing of endothelial cells exposed to hypoxia reveals the complexity of known and novel microRNAs.

In order to understand the role of microRNAs (miRNAs) in vascular physiopathology, we took advantage of deep-sequencing techniques to accurately and comprehensively profile the entire miRNA population expressed by endothelial cells exposed to hypoxia. SOLiD sequencing of small RNAs derived from human umbilical vein endothelial cells (HUVECs) exposed to 1% O2 or normoxia for 24 h yielded more than 22 million reads per library. A customized bioinformatic pipeline identified more than 400 annotated microRNA/microRNA* species with a broad abundance range: miR-21 and miR-126 totaled almost 40% of all miRNAs. A complex repertoire of isomiRs was found, displaying also 5' variations, potentially affecting target recognition. High-stringency bioinformatic analysis identified microRNA candidates, whose predicted pre-miRNAs folded into a stable hairpin. Validation of a subset by qPCR identified 18 high-confidence novel miRNAs as detectable in independent HUVEC cultures and associated to the RISC complex. The expression of two novel miRNAs was significantly down-modulated by hypoxia, while miR-210 was significantly induced. Gene ontology analysis of their predicted targets revealed a significant association to hypoxia-inducible factor signaling, cardiovascular diseases, and cancer. Overexpression of the novel miRNAs in hypoxic endothelial cells affected cell growth and confirmed the biological relevance of their down-modulation. In conclusion, deep-sequencing accurately profiled known, variant, and novel microRNAs expressed by endothelial cells in normoxia and hypoxia.