Paternal diet defines offspring chromatin state and intergenerational obesity.

The global rise in obesity has revitalized a search for genetic and epigenetic factors underlying the disease. We present a Drosophila model of paternal-diet-induced intergenerational metabolic reprogramming (IGMR) and identify genes required for its encoding in offspring. Intriguingly, we find that as little as 2 days of dietary intervention in fathers elicits obesity in offspring. Paternal sugar acts as a physiological suppressor of variegation, desilencing chromatin-state-defined domains in both mature sperm and in offspring embryos. We identify requirements for H3K9/K27me3-dependent reprogramming of metabolic genes in two distinct germline and zygotic windows. Critically, we find evidence that a similar system may regulate obesity susceptibility and phenotype variation in mice and humans. The findings provide insight into the mechanisms underlying intergenerational metabolic reprogramming and carry profound implications for our understanding of phenotypic variation and evolution.

Design and testing of a humanized porcine donor for xenotransplantation.

Recent human decedent model studies1,2 and compassionate xenograft use3 have explored the promise of porcine organs for human transplantation. To proceed to human studies, a clinically ready porcine donor must be engineered and its xenograft successfully tested in nonhuman primates. Here we describe the design, creation and long-term life-supporting function of kidney grafts from a genetically engineered porcine donor transplanted into a cynomolgus monkey model. The porcine donor was engineered to carry 69 genomic edits, eliminating glycan antigens, overexpressing human transgenes and inactivating porcine endogenous retroviruses. In vitro functional analyses showed that the edited kidney endothelial cells modulated inflammation to an extent that was indistinguishable from that of human endothelial cells, suggesting that these edited cells acquired a high level of human immune compatibility. When transplanted into cynomolgus monkeys, the kidneys with three glycan antigen knockouts alone experienced poor graft survival, whereas those with glycan antigen knockouts and human transgene expression demonstrated significantly longer survival time, suggesting the benefit of human transgene expression in vivo. These results show that preclinical studies of renal xenotransplantation could be successfully conducted in nonhuman primates and bring us closer to clinical trials of genetically engineered porcine renal grafts.

Multiomics of Tissue Extracellular Vesicles Identifies Unique Modulators of Atherosclerosis and Calcific Aortic Valve Stenosis.

BACKGROUND: Fewer than 50% of patients who develop aortic valve calcification have concomitant atherosclerosis, implying differential pathogenesis. Although circulating extracellular vesicles (EVs) act as biomarkers of cardiovascular diseases, tissue-entrapped EVs are associated with early mineralization, but their cargoes, functions, and contributions to disease remain unknown. METHODS: Disease stage-specific proteomics was performed on human carotid endarterectomy specimens (n=16) and stenotic aortic valves (n=18). Tissue EVs were isolated from human carotid arteries (normal, n=6; diseased, n=4) and aortic valves (normal, n=6; diseased, n=4) by enzymatic digestion, (ultra)centrifugation, and a 15-fraction density gradient validated by proteomics, CD63-immunogold electron microscopy, and nanoparticle tracking analysis. Vesiculomics, comprising vesicular proteomics and small RNA-sequencing, was conducted on tissue EVs. TargetScan identified microRNA targets. Pathway network analyses prioritized genes for validation in primary human carotid artery smooth muscle cells and aortic valvular interstitial cells. RESULTS: Disease progression drove significant convergence (P<0.0001) of carotid artery plaque and calcified aortic valve proteomes (2318 proteins). Each tissue also retained a unique subset of differentially enriched proteins (381 in plaques; 226 in valves; q<0.05). Vesicular gene ontology terms increased 2.9-fold (P<0.0001) among proteins modulated by disease in both tissues. Proteomics identified 22 EV markers in tissue digest fractions. Networks of proteins and microRNA targets changed by disease progression in both artery and valve EVs revealed shared involvement in intracellular signaling and cell cycle regulation. Vesiculomics identified 773 proteins and 80 microRNAs differentially enriched by disease exclusively in artery or valve EVs (q<0.05); multiomics integration found tissue-specific EV cargoes associated with procalcific Notch and Wnt signaling in carotid arteries and aortic valves, respectively. Knockdown of tissue-specific EV-derived molecules FGFR2, PPP2CA, and ADAM17 in human carotid artery smooth muscle cells and WNT5A, APP, and APC in human aortic valvular interstitial cells significantly modulated calcification. CONCLUSIONS: The first comparative proteomics study of human carotid artery plaques and calcified aortic valves identifies unique drivers of atherosclerosis versus aortic valve stenosis and implicates EVs in advanced cardiovascular calcification. We delineate a vesiculomics strategy to isolate, purify, and study protein and RNA cargoes from EVs entrapped in fibrocalcific tissues. Integration of vesicular proteomics and transcriptomics by network approaches revealed novel roles for tissue EVs in modulating cardiovascular disease.

Screening for Inhibitors of YAP Nuclear Localization Identifies Aurora Kinase A as a Modulator of Lung Fibrosis.

Idiopathic pulmonary fibrosis is a progressive lung disease with limited therapeutic options that is characterized by pathological fibroblast activation and aberrant lung remodeling with scar formation. YAP (Yes-associated protein) is a transcriptional coactivator that mediates mechanical and biochemical signals controlling fibroblast activation. We previously identified HMG-CoA (3-hydroxy-3-methylglutaryl coenzyme A) reductase inhibitors (statins) as YAP inhibitors based on a high-throughput small-molecule screen in primary human lung fibroblasts. Here we report that several Aurora kinase inhibitors were also identified from the top hits of this screen. MK-5108, a highly selective inhibitor for AURKA (Aurora kinase A), induced YAP phosphorylation and cytoplasmic retention and significantly reduced profibrotic gene expression in human lung fibroblasts. The inhibitory effect on YAP nuclear translocation and profibrotic gene expression is specific to inhibition of AURKA, but not Aurora kinase B or C, and is independent of the Hippo pathway kinases LATS1 and LATS2 (Large Tumor Suppressor 1 and 2). Further characterization of the effects of MK-5108 demonstrate that it inhibits YAP nuclear localization indirectly via effects on actin polymerization and TGFß (Transforming Growth Factor ß) signaling. In addition, MK-5108 treatment reduced lung collagen deposition in the bleomycin mouse model of pulmonary fibrosis. Our results reveal a novel role for AURKA in YAP-mediated profibrotic activity in fibroblasts and highlight the potential of small-molecule screens for YAP inhibitors for identification of novel agents with antifibrotic activity.

Genomic analyses identify molecular subtypes of pancreatic cancer.

Integrated genomic analysis of 456 pancreatic ductal adenocarcinomas identified 32 recurrently mutated genes that aggregate into 10 pathways: KRAS, TGF-ß, WNT, NOTCH, ROBO/SLIT signalling, G1/S transition, SWI-SNF, chromatin modification, DNA repair and RNA processing. Expression analysis defined 4 subtypes: (1) squamous; (2) pancreatic progenitor; (3) immunogenic; and (4) aberrantly differentiated endocrine exocrine (ADEX) that correlate with histopathological characteristics. Squamous tumours are enriched for TP53 and KDM6A mutations, upregulation of the TP63∆N transcriptional network, hypermethylation of pancreatic endodermal cell-fate determining genes and have a poor prognosis. Pancreatic progenitor tumours preferentially express genes involved in early pancreatic development (FOXA2/3, PDX1 and MNX1). ADEX tumours displayed upregulation of genes that regulate networks involved in KRAS activation, exocrine (NR5A2 and RBPJL), and endocrine differentiation (NEUROD1 and NKX2-2). Immunogenic tumours contained upregulated immune networks including pathways involved in acquired immune suppression. These data infer differences in the molecular evolution of pancreatic cancer subtypes and identify opportunities for therapeutic development.

Plasma MicroRNA Profiling of Plasmodium falciparum Biomass and Association with Severity of Malaria Disease.

Severe malaria (SM) is a major public health problem in malaria-endemic countries. Sequestration of Plasmodium falciparum-infected erythrocytes in vital organs and the associated inflammation leads to organ dysfunction. MicroRNAs (miRNAs), which are rapidly released from damaged tissues into the host fluids, constitute a promising biomarker for the prognosis of SM. We applied next-generation sequencing to evaluate the differential expression of miRNAs in SM and in uncomplicated malaria (UM) in children in Mozambique. Six miRNAs were associated with in vitro P. falciparum cytoadhesion, severity in children, and P. falciparum biomass. Relative expression of hsa-miR-4497 quantified by TaqMan-quantitative reverse transcription PCR was higher in plasma of children with SM than those with UM (p<0.048) and again correlated with P. falciparum biomass (p = 0.033). These findings suggest that different physiopathological processes in SM and UM lead to differential expression of miRNAs and suggest a pathway for assessing their prognostic value malaria.

Unification of miRNA and isomiR research: the mirGFF3 format and the mirtop API.

MOTIVATION: MicroRNAs (miRNAs) are small RNA molecules (∼22 nucleotide long) involved in post-transcriptional gene regulation. Advances in high-throughput sequencing technologies led to the discovery of isomiRs, which are miRNA sequence variants. While many miRNA-seq analysis tools exist, the diversity of output formats hinders accurate comparisons between tools and precludes data sharing and the development of common downstream analysis methods. RESULTS: To overcome this situation, we present here a community-based project, miRNA Transcriptomic Open Project (miRTOP) working towards the optimization of miRNA analyses. The aim of miRTOP is to promote the development of downstream isomiR analysis tools that are compatible with existing detection and quantification tools. Based on the existing GFF3 format, we first created a new standard format, mirGFF3, for the output of miRNA/isomiR detection and quantification results from small RNA-seq data. Additionally, we developed a command line Python tool, mirtop, to create and manage the mirGFF3 format. Currently, mirtop can convert into mirGFF3 the outputs of commonly used pipelines, such as seqbuster, isomiR-SEA, sRNAbench, Prost! as well as BAM files. Some tools have also incorporated the mirGFF3 format directly into their code, such as, miRge2.0, IsoMIRmap and OptimiR. Its open architecture enables any tool or pipeline to output or convert results into mirGFF3. Collectively, this isomiR categorization system, along with the accompanying mirGFF3 and mirtop API, provide a comprehensive solution for the standardization of miRNA and isomiR annotation, enabling data sharing, reporting, comparative analyses and benchmarking, while promoting the development of common miRNA methods focusing on downstream steps of miRNA detection, annotation and quantification. AVAILABILITY AND IMPLEMENTATION: https://github.com/miRTop/mirGFF3/ and https://github.com/miRTop/mirtop. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Huntington's disease brain-derived small RNAs recapitulate associated neuropathology in mice.

Progressive motor alterations and selective death of striatal medium spiny neurons (MSNs) are key pathological hallmarks of Huntington's disease (HD), a neurodegenerative condition caused by a CAG trinucleotide repeat expansion in the coding region of the huntingtin (HTT) gene. Most research has focused on the pathogenic effects of the resultant protein product(s); however, growing evidence indicates that expanded CAG repeats within mutant HTT mRNA and derived small CAG repeat RNAs (sCAG) participate in HD pathophysiology. The individual contribution of protein versus RNA toxicity to HD pathophysiology remains largely uncharacterized and the role of other classes of small RNAs (sRNA) that are strongly perturbed in HD is uncertain. Here, we demonstrate that sRNA produced in the putamen of HD patients (HD-sRNA-PT) are sufficient to induce HD pathology in vivo. Mice injected with HD-sRNA-PT show motor abnormalities, decreased levels of striatal HD-related proteins, disruption of the indirect pathway, and strong transcriptional abnormalities, paralleling human HD pathology. Importantly, we show that the specific blockage of sCAG mitigates HD-sRNA-PT neurotoxicity only to a limited extent. This observation prompted us to identify other sRNA species enriched in HD putamen with neurotoxic potential. We detected high levels of tRNA fragments (tRFs) in HD putamen, and we validated the neurotoxic potential of an Alanine derived tRF in vitro. These results highlight that HD-sRNA-PT are neurotoxic, and suggest that multiple sRNA species contribute to striatal dysfunction and general transcriptomic changes, favoring therapeutic strategies based on the blockage of sRNA-mediated toxicity.

Screening for YAP Inhibitors Identifies Statins as Modulators of Fibrosis.

Idiopathic pulmonary fibrosis is a lung disease with limited therapeutic options that is characterized by pathological fibroblast activation and aberrant lung remodeling with scar formation. YAP (Yes-associated protein) is a transcriptional coactivator that mediates mechanical and biochemical signals controlling fibroblast activation. In this study, we developed a high-throughput small-molecule screen for YAP inhibitors in primary human lung fibroblasts. Multiple HMG-CoA (hydroxymethylglutaryl-coenzyme A) reductase inhibitors (statins) were found to inhibit YAP nuclear localization via induction of YAP phosphorylation, cytoplasmic retention, and degradation. We further show that the mevalonate pathway regulates YAP activation, and that simvastatin treatment reduces fibrosis markers in activated human lung fibroblasts and in the bleomycin mouse model of pulmonary fibrosis. Finally, we show that simvastatin modulates YAP in vivo in mouse lung fibroblasts. Our results highlight the potential of small-molecule screens for YAP inhibitors and provide a mechanism for the antifibrotic activity of statins in idiopathic pulmonary fibrosis.

Functional Impact and Evolution of a Novel Human Polymorphic Inversion That Disrupts a Gene and Creates a Fusion Transcript.

Despite many years of study into inversions, very little is known about their functional consequences, especially in humans. A common hypothesis is that the selective value of inversions stems in part from their effects on nearby genes, although evidence of this in natural populations is almost nonexistent. Here we present a global analysis of a new 415-kb polymorphic inversion that is among the longest ones found in humans and is the first with clear position effects. This inversion is located in chromosome 19 and has been generated by non-homologous end joining between blocks of transposable elements with low identity. PCR genotyping in 541 individuals from eight different human populations allowed the detection of tag SNPs and inversion genotyping in multiple populations worldwide, showing that the inverted allele is mainly found in East Asia with an average frequency of 4.7%. Interestingly, one of the breakpoints disrupts the transcription factor gene ZNF257, causing a significant reduction in the total expression level of this gene in lymphoblastoid cell lines. RNA-Seq analysis of the effects of this expression change in standard homozygotes and inversion heterozygotes revealed distinct expression patterns that were validated by quantitative RT-PCR. Moreover, we have found a new fusion transcript that is generated exclusively from inverted chromosomes around one of the breakpoints. Finally, by the analysis of the associated nucleotide variation, we have estimated that the inversion was generated ~40,000-50,000 years ago and, while a neutral evolution cannot be ruled out, its current frequencies are more consistent with those expected for a deleterious variant, although no significant association with phenotypic traits has been found so far.

The small RNA content of human sperm reveals pseudogene-derived piRNAs complementary to protein-coding genes.

At the end of mammalian sperm development, sperm cells expel most of their cytoplasm and dispose of the majority of their RNA. Yet, hundreds of RNA molecules remain in mature sperm. The biological significance of the vast majority of these molecules is unclear. To better understand the processes that generate sperm small RNAs and what roles they may have, we sequenced and characterized the small RNA content of sperm samples from two human fertile individuals. We detected 182 microRNAs, some of which are highly abundant. The most abundant microRNA in sperm is miR-1246 with predicted targets among sperm-specific genes. The most abundant class of small noncoding RNAs in sperm are PIWI-interacting RNAs (piRNAs). Surprisingly, we found that human sperm cells contain piRNAs processed from pseudogenes. Clusters of piRNAs from human testes contain pseudogenes transcribed in the antisense strand and processed into small RNAs. Several human protein-coding genes contain antisense predicted targets of pseudogene-derived piRNAs in the male germline and these piRNAs are still found in mature sperm. Our study provides the most extensive data set and annotation of human sperm small RNAs to date and is a resource for further functional studies on the roles of sperm small RNAs. In addition, we propose that some of the pseudogene-derived human piRNAs may regulate expression of their parent gene in the male germline.

Specific small-RNA signatures in the amygdala at premotor and motor stages of Parkinson's disease revealed by deep sequencing analysis.

MOTIVATION: Most computational tools for small non-coding RNAs (sRNA) sequencing data analysis focus in microRNAs (miRNAs), overlooking other types of sRNAs that show multi-mapping hits. Here, we have developed a pipeline to non-redundantly quantify all types of sRNAs, and extract patterns of expression in biologically defined groups. We have used our tool to characterize and profile sRNAs in post-mortem brain samples of control individuals and Parkinson's disease (PD) cases at early-premotor and late-symptomatic stages. RESULTS: Clusters of co-expressed sRNAs mapping onto tRNAs significantly separated premotor and motor cases from controls. A similar result was obtained using a matrix of miRNAs slightly varying in sequence (isomiRs). The present framework revealed sRNA alterations at premotor stages of PD, which might reflect initial pathogenic perturbations. This tool may be useful to discover sRNA expression patterns linked to different biological conditions. AVAILABILITY AND IMPLEMENTATION: The full code is available at http://github.com/lpantano/seqbuster CONTACT: lpantano@hsph.harvard.edu or eulalia.marti@crg.eu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Comparative analysis of LIN28-RNA binding sites identified at single nucleotide resolution.

It remains a formidable challenge to characterize the diverse complexes of RNA binding proteins and their targets. While crosslink and immunoprecipitation (CLIP) methods are powerful techniques that identify RNA targets on a global scale, the resolution and consistency of these methods is a matter of debate. Here we present a comparative analysis of LIN28-pre-let-7 UV-induced crosslinking using a tandem mass spectrometry (MS/MS) and deep sequencing interrogation of in vitro crosslinked complexes. Interestingly, analyses by the two methods diverge in their identification of crosslinked nucleotide identity - whereas bioinformatics and sequencing analyses suggest guanine in mammalian cells, MS/MS identifies uridine. This work suggests the need for comprehensive analysis and validation of crosslinking methodologies.

InvFEST, a database integrating information of polymorphic inversions in the human genome.

The newest genomic advances have uncovered an unprecedented degree of structural variation throughout genomes, with great amounts of data accumulating rapidly. Here we introduce InvFEST (http://invfestdb.uab.cat), a database combining multiple sources of information to generate a complete catalogue of non-redundant human polymorphic inversions. Due to the complexity of this type of changes and the underlying high false-positive discovery rate, it is necessary to integrate all the available data to get a reliable estimate of the real number of inversions. InvFEST automatically merges predictions into different inversions, refines the breakpoint locations, and finds associations with genes and segmental duplications. In addition, it includes data on experimental validation, population frequency, functional effects and evolutionary history. All this information is readily accessible through a complete and user-friendly web report for each inversion. In its current version, InvFEST combines information from 34 different studies and contains 1092 candidate inversions, which are categorized based on internal scores and manual curation. Therefore, InvFEST aims to represent the most reliable set of human inversions and become a central repository to share information, guide future studies and contribute to the analysis of the functional and evolutionary impact of inversions on the human genome.

A pathogenic mechanism in Huntington's disease involves small CAG-repeated RNAs with neurotoxic activity.

Huntington's disease (HD) is an autosomal dominantly inherited disorder caused by the expansion of CAG repeats in the Huntingtin (HTT) gene. The abnormally extended polyglutamine in the HTT protein encoded by the CAG repeats has toxic effects. Here, we provide evidence to support that the mutant HTT CAG repeats interfere with cell viability at the RNA level. In human neuronal cells, expanded HTT exon-1 mRNA with CAG repeat lengths above the threshold for complete penetrance (40 or greater) induced cell death and increased levels of small CAG-repeated RNAs (sCAGs), of ≈21 nucleotides in a Dicer-dependent manner. The severity of the toxic effect of HTT mRNA and sCAG generation correlated with CAG expansion length. Small RNAs obtained from cells expressing mutant HTT and from HD human brains significantly decreased neuronal viability, in an Ago2-dependent mechanism. In both cases, the use of anti-miRs specific for sCAGs efficiently blocked the toxic effect, supporting a key role of sCAGs in HTT-mediated toxicity. Luciferase-reporter assays showed that expanded HTT silences the expression of CTG-containing genes that are down-regulated in HD. These results suggest a possible link between HD and sCAG expression with an aberrant activation of the siRNA/miRNA gene silencing machinery, which may trigger a detrimental response. The identification of the specific cellular processes affected by sCAGs may provide insights into the pathogenic mechanisms underlying HD, offering opportunities to develop new therapeutic approaches.

Genome assembly comparison identifies structural variants in the human genome.

Numerous types of DNA variation exist, ranging from SNPs to larger structural alterations such as copy number variants (CNVs) and inversions. Alignment of DNA sequence from different sources has been used to identify SNPs and intermediate-sized variants (ISVs). However, only a small proportion of total heterogeneity is characterized, and little is known of the characteristics of most smaller-sized (<50 kb) variants. Here we show that genome assembly comparison is a robust approach for identification of all classes of genetic variation. Through comparison of two human assemblies (Celera's R27c compilation and the Build 35 reference sequence), we identified megabases of sequence (in the form of 13,534 putative non-SNP events) that were absent, inverted or polymorphic in one assembly. Database comparison and laboratory experimentation further demonstrated overlap or validation for 240 variable regions and confirmed >1.5 million SNPs. Some differences were simple insertions and deletions, but in regions containing CNVs, segmental duplication and repetitive DNA, they were more complex. Our results uncover substantial undescribed variation in humans, highlighting the need for comprehensive annotation strategies to fully interpret genome scanning and personalized sequencing projects.

A highly expressed miR-101 isomiR is a functional silencing small RNA.

BACKGROUND: MicroRNAs (miRNAs) are short non-coding regulatory RNAs that control gene expression usually producing translational repression and gene silencing. High-throughput sequencing technologies have revealed heterogeneity at length and sequence level for the majority of mature miRNAs (IsomiRs). Most isomiRs can be explained by variability in either Dicer1 or Drosha cleavage during miRNA biogenesis at 5' or 3' of the miRNA (trimming variants). Although isomiRs have been described in different tissues and organisms, their functional validation as modulators of gene expression remains elusive. Here we have characterized the expression and function of a highly abundant miR-101 5'-trimming variant (5'-isomiR-101). RESULTS: The analysis of small RNA sequencing data in several human tissues and cell lines indicates that 5'-isomiR-101 is ubiquitously detected and a highly abundant, especially in the brain. 5'-isomiR-101 was found in Ago-2 immunocomplexes and complementary approaches showed that 5'-isomiR-101 interacted with different members of the silencing (RISC) complex. In addition, 5'-isomiR-101 decreased the expression of five validated miR-101 targets, suggesting that it is a functional variant. Both the binding to RISC members and the degree of silencing were less efficient for 5'-isomiR-101 compared with miR-101. For some targets, both miR-101 and 5'-isomiR-101 significantly decreased protein expression with no changes in the respective mRNA levels. Although a high number of overlapping predicted targets suggest similar targeted biological pathways, a correlation analysis of the expression profiles of miR-101 variants and predicted mRNA targets in human brains at different ages, suggest specific functions for miR-101- and 5'-isomiR-101. CONCLUSIONS: These results suggest that isomiRs are functional variants and further indicate that for a given miRNA, the different isomiRs may contribute to the overall effect as quantitative and qualitative fine-tuners of gene expression.

Microarray and deep sequencing cross-platform analysis of the mirRNome and isomiR variation in response to epidermal growth factor.

BACKGROUND: Epidermal Growth Factor (EGF) plays an important function in the regulation of cell growth, proliferation, and differentiation by binding to its receptor (EGFR) and providing cancer cells with increased survival responsiveness. Signal transduction carried out by EGF has been extensively studied at both transcriptional and post-transcriptional levels. Little is known about the involvement of microRNAs (miRNAs) in the EGF signaling pathway. miRNAs have emerged as major players in the complex networks of gene regulation, and cancer miRNA expression studies have evidenced a direct involvement of miRNAs in cancer progression. RESULTS: In this study, we have used an integrative high content analysis approach to identify the specific miRNAs implicated in EGF signaling in HeLa cells as potential mediators of cancer mediated functions. We have used microarray and deep-sequencing technologies in order to obtain a global view of the EGF miRNA transcriptome with a robust experimental cross-validation. By applying a procedure based on Rankprod tests, we have delimited a solid set of EGF-regulated miRNAs. After validating regulated miRNAs by reverse transcription quantitative PCR, we have derived protein networks and biological functions from the predicted targets of the regulated miRNAs to gain insight into the potential role of miRNAs in EGF-treated cells. In addition, we have analyzed sequence heterogeneity due to editing relative to the reference sequence (isomiRs) among regulated miRNAs. CONCLUSIONS: We propose that the use of global genomic miRNA cross-validation derived from high throughput technologies can be used to generate more reliable datasets inferring more robust networks of co-regulated predicted miRNA target genes.

A non-biased framework for the annotation and classification of the non-miRNA small RNA transcriptome.

MOTIVATION: Recent progress in high-throughput sequencing technologies has largely contributed to reveal a highly complex landscape of small non-coding RNAs (sRNAs), including novel non-canonical sRNAs derived from long non-coding RNA, repeated elements, transcription start sites and splicing site regions among others. The published frameworks for sRNA data analysis are focused on miRNA detection and prediction, ignoring further information in the dataset. As a consequence, tools for the identification and classification of the sRNAs not belonging to miRNA family are currently lacking. RESULTS: Here, we present, SeqCluster, an extension of the currently available SeqBuster tool to identify and analyze at different levels the sRNAs not annotated or predicted as miRNAs. This new module deals with sequences mapping onto multiple locations and permits a highly versatile and user-friendly interaction with the data in order to easily classify sRNA sequences with a putative functional importance. We were able to detect all known classes of sRNAs described to date using SeqCluster with different sRNA datasets.

SeqBuster, a bioinformatic tool for the processing and analysis of small RNAs datasets, reveals ubiquitous miRNA modifications in human embryonic cells.

High-throughput sequencing technologies enable direct approaches to catalog and analyze snapshots of the total small RNA content of living cells. Characterization of high-throughput sequencing data requires bioinformatic tools offering a wide perspective of the small RNA transcriptome. Here we present SeqBuster, a highly versatile and reliable web-based toolkit to process and analyze large-scale small RNA datasets. The high flexibility of this tool is illustrated by the multiple choices offered in the pre-analysis for mapping purposes and in the different analysis modules for data manipulation. To overcome the storage capacity limitations of the web-based tool, SeqBuster offers a stand-alone version that permits the annotation against any custom database. SeqBuster integrates multiple analyses modules in a unique platform and constitutes the first bioinformatic tool offering a deep characterization of miRNA variants (isomiRs). The application of SeqBuster to small-RNA datasets of human embryonic stem cells revealed that most miRNAs present different types of isomiRs, some of them being associated to stem cell differentiation. The exhaustive description of the isomiRs provided by SeqBuster could help to identify miRNA-variants that are relevant in physiological and pathological processes. SeqBuster is available at http://estivill_lab.crg.es/seqbuster.