Search

Over the last few years, the number of microRNAs in the human genome has become a controversially debated issue. Several publications reported thousands of putative novel microRNAs not included in the curated microRNA gene database MirGeneDB and the repository miRBase. Recently, by using sequencing of ∼300 human tissues and cell lines, the human RNA atlas, an expanded inventory of human RNA annotations, was published, reporting thousands of putative microRNAs. We, the developers of established microRNA prediction tools and hosts of MirGeneDB, raise concerns about the frequently applied prediction and functional validation strategies, briefly discussing the drawbacks of false positive detections. By means of quantifying well-established biogenesis-derived features, we show that the reported novel microRNAs essentially represent false-positives and argue that the human microRNA complement, at about 550 microRNA genes, is already near complete. Output of available tools must be curated as false predictions will misguide scientists looking for biomarkers or therapeutic targets.

Whole genome duplications (WGDs) are major events that drastically reshape genome architecture and are causally associated with organismal innovations and radiations1. The 2R Hypothesis suggests that two WGD events (1R and 2R) occurred during early vertebrate evolution2,3. However, the veracity and timing of the 2R event relative to the divergence of gnathostomes (jawed vertebrates) and cyclostomes (jawless hagfishes and lampreys) is unresolved4-6 and whether these WGD events underlie vertebrate phenotypic diversification remains elusive7. Here we present the genome of the inshore hagfish, Eptatretus burgeri. Through comparative analysis with lamprey and gnathostome genomes, we reconstruct the early events in cyclostome genome evolution, leveraging insights into the ancestral vertebrate genome. Genome-wide synteny and phylogenetic analyses support a scenario in which 1R occurred in the vertebrate stem-lineage during the early Cambrian, and the 2R event occurred in the gnathostome stem-lineage in the late Cambrian after its divergence from cyclostomes. We find that the genome of stem-cyclostomes experienced two additional, independent genome duplications (herein CR1 and CR2). Functional genomic and morphospace analyses demonstrate that WGD events generally contribute to developmental evolution with similar changes in the regulatory genome of both vertebrate groups. However, appreciable morphological diversification occurred only after the 2R event, questioning the general expectation that WGDs lead to leaps of morphological complexity7.

Soft-bodied cephalopods such as the octopus are exceptionally intelligent invertebrates with a highly complex nervous system that evolved independently from vertebrates. Because of elevated RNA editing in their nervous tissues, we hypothesized that RNA regulation may play a major role in the cognitive success of this group. We thus profiled mRNAs and small RNAs in 18 tissues of the common octopus. We show that the major RNA innovation of soft-bodied cephalopods is a massive expansion of the miRNA gene repertoire. These novel miRNAs were primarily expressed in neuronal tissues, during development, and had conserved and thus likely functional target sites. The only comparable miRNA expansions happened, strikingly, in vertebrates. Thus, we propose that miRNAs are intimately linked to the evolution of complex animal brains.One-Sentence SummarymiRNAs are deeply linked to the emergence of complex brains.

Terrestrial mammals range over six orders of magnitude in size, while average cellular, or specific metabolic rate(sMR) varies only ∼5-fold, constrained by the strict matching of protein synthesis to cell size. Protein synthesis, the single most costly metabolic activity, accounts for 20% of ATP turnover, and the speed of translation can be increased only at the expense of increased energy dissipation due to misreading errors and stochastic translation noise. We hypothesized that global microRNA activity evolved under the same constraints assMRand would therefore exhibit a similar pattern of variation. This expectation was met by the number of microRNA families: based on the manually curated microRNA gene database MirGeneDB, the acquisition of microRNA families (mirFam) accelerated two-fold in late-branching mammals with body temperatures (Tb) > 36°C, relative to early-branching mammals (Tb< 36°C), independent of phylogenetic distance. WhensMRis fit tomirFamby maximum likelihood, the variation is homoscedastic, allowing us to compare models for the evolution ofsMRin relation tomirFam. WithmirFamas the predictor, an Ornstein-Uhlenbeck (OU) process of stabilizing selection accounted better for the biphasic evolution ofsMRthan did the corresponding Brownian motion model, and was optimized whensMRwas scaled toM0.75. OU simulations with an adaptive shift at the divergence of Boreoeutheria predicted 95% of the variation insMR. We infer that the advent of placentation led to the emergence of an adaptive optimum characterized by higher body temperatures, fasterMRs, and heightened global microRNA activity.

SummaryThe annotation of microRNAs, an important class of post-transcriptional regulators, depends on the availability of transcriptomics data and expert knowledge. This led to a large gap between novel genomes made available and high-quality microRNA complements. Using >16,000 microRNAs from the manually curated microRNA gene database MirGeneDB, we generated trained covariance models for all conserved microRNA families. These models are available in MirMachine, our new tool for the annotation of conserved microRNA complements from genomes only. We successfully applied MirMachine to a wide range of animal species, including those with very large genomes, additional genome duplications and extinct species, where smallRNA sequencing will be hard to achieve. We further describe a microRNA score of expected microRNAs that can be used to assess the completeness of genome assemblies. MirMachine closes a long-persisting gap in the microRNA field facilitating automated genome annotation pipelines and deeper studies on the evolution of genome regulation, even in extinct organisms.HighlightsAn annotation pipeline using trained covariance models of microRNA familiesEnables massive parallel annotation of microRNA complements of genomesMirMachine creates meaningful annotations for very large and extinct genomesmicroRNA score to assess genome assembly completenessGraphical abstract

The evolution of specialized cell-types is a long-standing interest of biologists, but given the deep time-scales very difficult to reconstruct or observe. microRNAs have been linked to the evolution of cellular complexity and may inform on specialization. The endothelium is a vertebrate specific specialization of the circulatory system that enabled a critical new level of vasoregulation. The evolutionary origin of these endothelial cells is unclear. We hypothesized that Mir-126, an endothelial cell-specific microRNA may be informative.We here reconstruct the evolutionary history of Mir-126. Mir-126 likely appeared in the last common ancestor of vertebrates and tunicates, a species without an endothelium, within an intron of the evolutionary much older EGF Like Domain Multiple (Egfl) locus. Mir-126 has a complex evolutionary history due to duplications and losses of both the host gene and the microRNA. Taking advantage of the strong evolutionary conservation of the microRNA among Olfactores, and using RNA in situ hybridization (RISH), we localized Mir-126 in the tunicateCiona robusta. We found exclusive expression of the mature Mir-126 in granular amebocytes, supporting a long-proposed scenario that endothelial cells arose from hemoblasts, a type of proto-endothelial amoebocyte found throughout invertebrates.This observed change of expression of Mir-126 from proto-endothelial amoebocytes in the tunicate to endothelial cells in vertebrates is the first direct observation of the evolution of a cell-type in relation to microRNA expression indicating that microRNAs can be a prerequisite of cell-type evolution.Research Highlightsdirect observation of cell-type evolutionhigh conservation of sequence enables for simple RISH experiment of expressionMir-126 follows the evolution of hematopoetic cells to endothelial cells

The evolutionary origins of Bilateria remain enigmatic. One of the more enduring proposals highlights similarities between a cnidarian-like planula larva and simple acoel-like flatworms. This idea is based in part on the view of the Xenacoelomorpha as an outgroup to all other bilaterians which are themselves designated the Nephrozoa (protostomes and deuterostomes). Genome data, which can help to elucidate phylogenetic relationships and provide important comparative data, remain sparse for early branching bilaterians. Here we assemble and analyse the genome of the simple, marine xenacoelomorphXenoturbella bocki, a key species for our understanding of early bilaterian and deuterostome evolution. Our highly contiguous genome assembly ofX. bockihas a size of ∼110 Mbp in 18 chromosome like scaffolds, with repeat content, and intron, exon and intergenic space comparable to other bilaterian invertebrates. We findX. bockito have a similar number of genes to other bilaterians and to have retained ancestral metazoan synteny. Key bilaterian signalling pathways are also largely complete and most bilaterian miRNAs are present. We conclude thatX. bockihas a complex genome typical of bilaterians, in contrast to the apparent simplicity of its body plan. Overall, our data do not provide evidence supporting the idea that Xenacoelomorpha are a primitively simple outgroup to other bilaterians and gene presence/absence data support a relationship with Ambulacraria.

The tropical liver fluke Fasciola gigantica affects livestock and humans in many Asian countries, large parts of Africa, and parts of Europe. Despite the public health and economic impacts of F. gigantica, understanding of F. gigantica biology and how the complex lifecycle of this liver fluke is transcriptionally regulated remain unknown. Here, we tested the hypothesis that the regulatory small non-coding RNAs (sncRNAs), microRNAs (miRNAs) and tRNA-derived fragments (tRFs) play roles in the adaptation of F. gigantica to its intermediate and definitive hosts. We sequenced sncRNAs of eight lifecycle stages of F. gigantica. In total, 56 miRNAs from 33 conserved families and four Fasciola-specific miRNAs were identified. Expression analysis of miRNAs suggested clear stage-related patterns. By leveraging the existing transcriptomic data, we predicted a miRNA-based regulation of metabolism, transport, growth and developmental processes. Also, by comparing miRNA complement of F. gigantica with that of Fasciola hepatica, we detected a high level of conservation and identified differences in some miRNAs, which can be used to distinguish the two species. Moreover, we found that tRFs at each lifecycle stage were predominantly derived by tRNA-Lys and tRNA-Gly at 5' half sites, but relatively high expression was related to the buffalo-infecting stages. Taken together, we provided a comprehensive overview of the dynamic transcriptional changes of small RNAs that occur during the developmental stages of F. gigantica. This global analysis of F. gigantica lifecycle stages revealed new roles of miRNAs and tRFs in parasite development and will facilitate future research into understanding of fasciolosis pathobiology.

DNA sequencing is the current key technology for historic or ancient biological samples and has led to many exciting discoveries in the field of paleogenomics. However, functional insights into tissue identity, cellular composition, or gene regulation cannot be gained from DNA. Recent analyses have shown that, under favorable conditions, RNA can also be sequenced from ancient samples, enabling studies at the transcriptomic and regulatory level. Analyzing ancient RNA data from a Pleistocene canid, we find hundreds of intact microRNAs that are taxonomically informative, show tissue specificity and have functionally predictive characteristics. With an extraordinary age of 14,300 yr, these microRNA sequences are by far the oldest ever reported. The authenticity of the sequences is further supported by (i) the presence of canid/Caniformia-specific sequences that never evolved outside of this clade, (ii) tissue-specific expression patterns (cartilage, liver, and muscle) that resemble those of modern dogs, and (iii) RNA damage patterns that are clearly distinct from those of fresh samples. By performing computational microRNA-target enrichment analyses on the ancient sequences, we predict microRNA functions consistent with their tissue pattern of expression. For instance, we find a liver-specific microRNA that regulates carbohydrate metabolism and starvation responses in canids. In summary, we show that straightforward paleotranscriptomic microRNA analyses can give functional glimpses into tissue identity, cellular composition, and gene regulatory activity of ancient samples and biological processes that took place in the Pleistocene, thus holding great promise for deeper insights into gene regulation in extinct animals based on ancient RNA sequencing.

MicroRNAs (miRNAs) regulate the levels of translation of messenger RNAs (mRNAs). At present, the major parameter that can explain the selection of the target mRNA and the efficiency of translation repression is the base pairing between the ‘seed’ region of the miRNA and its counterpart mRNA1. Here we use R1ρ relaxation-dispersion nuclear magnetic resonance2 and molecular simulations3 to reveal a dynamic switch—based on the rearrangement of a single base pair in the miRNA–mRNA duplex—that elongates a weak five-base-pair seed to a complete seven-base-pair seed. This switch also causes coaxial stacking of the seed and supplementary helix fitting into human Argonaute 2 protein (Ago2), reminiscent of an active state in prokaryotic Ago4,5. Stabilizing this transient state leads to enhanced repression of the target mRNA in cells, revealing the importance of this miRNA–mRNA structure. Our observations tie together previous findings regarding the stepwise miRNA targeting process from an initial ‘screening’ state to an ‘active’ state, and unveil the role of the RNA duplex beyond the seed in Ago2. Repression of a messenger RNA by a cognate microRNA depends not only on complementary base pairing, but also on the rearrangement of a single base pair, producing a conformation that fits better within the human Ago2 protein.

Although microRNAs (miRNAs) contribute to all hallmarks of cancer, miRNA dysregulation in metastasis remains poorly understood. The aim of this work was to reliably identify miRNAs associated with metastatic progression of colorectal cancer (CRC) using novel and previously published next-generation sequencing (NGS) datasets generated from 268 samples of primary (pCRC) and metastatic CRC (mCRC; liver, lung and peritoneal metastases) and tumor adjacent tissues. Differential expression analysis was performed using a meticulous bioinformatics pipeline, including only bona fide miRNAs, and utilizing miRNA-tailored quality control and processing. Five miRNAs were identified as up-regulated at multiple metastatic sites Mir-210_3p, Mir-191_5p, Mir-8-P1b_3p [mir-141–3p], Mir-1307_5p and Mir-155_5p. Several have previously been implicated in metastasis through involvement in epithelial-to-mesenchymal transition and hypoxia, while other identified miRNAs represent novel findings. The use of a publicly available pipeline facilitates reproducibility and allows new datasets to be added as they become available. The set of miRNAs identified here provides a reliable starting-point for further research into the role of miRNAs in metastatic progression.

Background Colorectal liver metastasis (CLM) is a leading cause of colorectal cancer mortality, and the response to immune checkpoint inhibition (ICI) in microsatellite-stable CRC has been disappointing. Administration of cytotoxic chemotherapy may cause increased density of tumor-infiltrating T cells, which has been associated with improved response to ICI. This study aimed to quantify and characterize T-cell infiltration in CLM using T-cell receptor (TCR) repertoire sequencing. Eighty-five resected CLMs from patients included in the Oslo CoMet study were subjected to TCR repertoire sequencing. Thirty-five and 15 patients had received neoadjuvant chemotherapy (NACT) within a short or long interval, respectively, prior to resection, while 35 patients had not been exposed to NACT. T-cell fractions were calculated, repertoire clonality was analyzed based on Hill evenness curves, and TCR sequence convergence was assessed using network analysis. Results Increased T-cell fractions (10.6% vs. 6.3%) were detected in CLMs exposed to NACT within a short interval prior to resection, while modestly increased clonality was observed in NACT-exposed tumors independently of the timing of NACT administration and surgery. While private clones made up >90% of detected clones, network connectivity analysis revealed that public clones contributed the majority of TCR sequence convergence. Conclusions TCR repertoire sequencing can be used to quantify T-cell infiltration and clonality in clinical samples. This study provides evidence to support chemotherapy-driven T-cell clonal expansion in CLM in a clinical context.

The telomere to telomere (T2T) genome project discovered and mapped ∼240 million additional base pairs of primarily telomeric and centromeric reads. Much of this sequence was comprised of satellite sequences and large segmental duplications. We evaluated the extent to which human bona fide microRNAs (miRNAs) may be found in additional paralogous genomic loci or if previously undescribed microRNAs are present in these newly sequenced regions of the human genome. New genomic regions of the T2T project spanning ∼240 million bp of sequence were obtained and evaluated by blastn for the human miRNAs contained in MirGeneDB2.0 (N=556) and miRBase (N = 1917) along with all species of MirGeneDB2.0 miRNAs (N=10,899). Additionally, bowtie was used to compare unmapped reads from >4,000 primary cell samples to the new T2T sequence. Based on sequence and structure, no bona fide miRNAs were identified. Ninety-seven miRNAs of questionable authenticity (frequently known repeat elements) were identified from the miRBase dataset across the newly described regions of the human genome. These 97 represent only 51 miRNA families due to paralogy of highly similar miRNAs such as 24 members of the hsa-mir-548 family. Altogether, this data strongly supports our having identified widely expressed bona fide miRNAs in the human genome and move us further toward the completion of human miRNA discovery.

Whole-genome duplications (WGDs) have long been considered the causal mechanism underlying dramatic increases to morphological complexity due to the neo-functionalization of paralogs generated during these events. Nonetheless, an alternative hypothesis suggests that behind the retention of most paralogs is not neo-functionalization, but instead the degree of the inter-connectivity of the intended gene product, as well as the mode of the WGD itself. Here, we explore both the causes and consequences of WGD by examining the distribution, expression, and molecular evolution of microRNAs (miRNAs) in both gnathostome vertebrates as well as chelicerate arthropods. We find that although the number of miRNA paralogs tracks the number of WGDs experienced within the lineage, few of these paralogs experienced changes to the seed sequence, and thus are functionally equivalent relative to their mRNA targets. Nonetheless, in gnathostomes, although the retention of paralogs following the 1R autotetraploidization event is similar across the two subgenomes, the paralogs generated by the gnathostome 2R allotetraploidization event are retained in higher numbers on one subgenome relative to the second, with the miRNAs found on the preferred subgenome showing both higher expression of mature miRNA transcripts and slower molecular evolution of the precursor miRNA sequences. Importantly, WGDs do not result in the creation of miRNA novelty, nor do WGDs correlate to increases in complexity. Instead, it is the number of miRNA seed sequences in the genome itself that not only better correlate to instances in complexification, but also mechanistically explain why complexity increases when new miRNA families are established.

DNA sequencing is the current key technology for historic or ancient biological samples and has led to many exciting discoveries in the field of paleogenomics. However, functional insights into tissue identity, cellular composition, or gene regulation cannot be gained from DNA. Recent analyses have shown that, under favorable conditions, RNA can also be sequenced from ancient samples, enabling studies at the transcriptomic and regulatory level. Analyzing ancient RNA data from a Pleistocene canid, we find hundreds of intact microRNAs that are taxonomically informative, show tissue specificity and have functionally predictive characteristics. With an extraordinary age of 14,300 yr, these microRNA sequences are by far the oldest ever reported. The authenticity of the sequences is further supported by (i) the presence of canid/Caniformia-specific sequences that never evolved outside of this clade, (ii) tissue-specific expression patterns (cartilage, liver, and muscle) that resemble those of modern dogs, and (iii) RNA damage patterns that are clearly distinct from those of fresh samples. By performing computational microRNA-target enrichment analyses on the ancient sequences, we predict microRNA functions consistent with their tissue pattern of expression. For instance, we find a liver-specific microRNA that regulates carbohydrate metabolism and starvation responses in canids. In summary, we show that straightforward paleotranscriptomic microRNA analyses can give functional glimpses into tissue identity, cellular composition, and gene regulatory activity of ancient samples and biological processes that took place in the Pleistocene, thus holding great promise for deeper insights into gene regulation in extinct animals based on ancient RNA sequencing.

Whole genome duplications (WGDs) have long been considered the causal mechanism underlying the dramatic increase in vertebrate morphological complexity relative to invertebrates. This is due to the retention and neo-functionalization of paralogues generated during these events, evolving new regulatory circuits, and ultimately morphological novelty. Nonetheless, an alternative hypothesis suggests that behind the retention of most paralogues is not neo-functionalization, but instead the degree of the inter-connectivity of the intended gene product, as well as the mode of the WGD itself. Here, we explore both the causes and consequences of WGD by examining the distribution, expression, and molecular evolution of microRNAs (miRNAs) in both gnathostome vertebrates as well as chelicerate arthropods. We find that although the number of miRNA paralogues tracks the number of WGDs experienced within the lineage, few of these paralogues experienced changes to the seed sequence, and thus are functionally equivalent relative to their mRNA targets. Nonetheless, the paralogues generated by the gnathostome 2R allotetraploidization event are retained in higher numbers on one sub-genome relative the second, with the miRNAs found on the preferred set of paralogons showing both higher expression of mature miRNA transcripts and slower molecular evolution of the precursor miRNA sequences. Importantly, WGDs do not result in the creation of miRNA novelty, nor do WGDs correlate to increases in complexity. Instead, it is the number of miRNA seed sequences in the genome itself that not only better correlate to instances in complexification, but also mechanistically explain why complexity increases when new miRNA families are established.

Since the original publication of sRNAtoolbox in 2015, small RNA research experienced notable advances in different directions. New protocols for small RNA sequencing have become available to address important issues such as adapter ligation bias, PCR amplification artefacts or to include internal controls such as spike-in sequences. New microRNA reference databases were developed with different foci, either prioritizing accuracy (low number of false positives) or completeness (low number of false negatives). Additionally, other small RNA molecules as well asmicroRNA sequence and length variants (isomiRs) have continued to gain importance. Finally, the number of microRNA sequencing studies deposited in GEO nearly triplicated from 2014 (280) to 2018 (764). These developments imply that fast and easy-to-use tools for expression profiling and subsequent downstream analysis of miRNAseq data are essential to many researchers. Key features in this sRNAtoolbox release include addition of all major RNA library preparation protocols to sRNAbench and improvements in sRNAde, a tool that summarizes several aspects of small RNA sequencing studies including the detection of consensus differential expression. A special emphasis was put on the user-friendliness of the tools, for instance sRNAbench now supports parallel launching of several jobs to improve reproducibility and user time efficiency., European Union [765492 to M.H.]; Spanish Government [AGL2017-88702-C2-2-R to M.H., J.L.O.]; Instituto de Salud Carlos III, FEDER funds [PIE16/00045 to J.A.M.]; Chair ‘Doctors Galera-Requena in cancer stem cell research’ to JMA and by the Ministry of Education of Spain [FPU13/05662 to R.L., IFI16/00041 to E.A.]; Strategic Research Area (SFO) program of the Swedish Research Council (to V.R.) through Stockholm University (to B.F.). Funding for open access charge: SpanishGovernment [AGL2017-88702-C2-2-R].

High-throughput sequencing for micro-RNAs (miRNAs) to obtain expression estimates is a central method of molecular biology. Surprisingly, there are a number of different approaches to converting sequencing output into micro-RNA counts. Each has their own strengths and biases that impact on the final data that can be obtained from a sequencing run. This chapter serves to make the reader aware of the trade-offs one must consider in analyzing small RNA sequencing data. It then compares two methods, miRge2.0 and the sRNAbench and the steps utilized to output data from their tools.

Made available in DSpace on 2021-06-25T10:55:23Z (GMT). No. of bitstreams: 0 Previous issue date: 2021-12-01 Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) Background: During vertebrate evolution, the heart has undergone remarkable changes that lead to morphophysiological differences in the fully formed heart of these species, such as chamber septation, heart rate frequency, blood pressure, and cardiac output volume. Despite these differences, the heart developmental process is guided by a core gene set conserved across vertebrates. Nonetheless, the regulatory mechanisms controlling the expression of genes involved in heart development and maintenance are largely uncharted. MicroRNAs (miRNAs) have been described as important regulatory elements in several biological processes, including heart biology. These small RNA molecules are broadly conserved in sequence and genomic context in metazoans. Mutations may occur in miRNAs and/or genes that contribute to the establishment of distinct repertoires of miRNA-target interactions, thereby favoring the differential control of gene expression and, consequently, the origin of novel phenotypes. In fact, several studies showed that miRNAs are integrated into genetic regulatory networks (GRNs) governing specific developmental programs and diseases. However, studies integrating miRNAs in vertebrate heart GRNs under an evolutionary perspective are still scarce. Results: We comprehensively examined and compared the heart miRNome of 20 species representatives of the five major vertebrate groups. We found 54 miRNA families with conserved expression and a variable number of miRNA families with group-specific expression in fishes, amphibians, reptiles, birds, and mammals. We also detected that conserved miRNAs present higher expression levels and a higher number of targets, whereas the group-specific miRNAs present lower expression levels and few targets. Conclusions: Both the conserved and group-specific miRNAs can be considered modulators orchestrating the core and peripheral genes of heart GRNs of vertebrates, which can be related to the morphophysiological differences and similarities existing in the heart of distinct vertebrate groups. We propose a hypothesis to explain evolutionary differences in the putative functional roles of miRNAs in the heart GRNs analyzed. Furthermore, we present new insights into the molecular mechanisms that could be helping modulate the diversity of morphophysiology in the heart organ of vertebrate species. Laboratório Especial de Toxinologia Aplicada (LETA) CeTICS Instituto Butantan Department of Chemical and Biological Sciences Institute of Biosciences of Botucatu São Paulo State University (UNESP) Department of Biophysics and Pharmacology Institute of Biosciences of Botucatu São Paulo State University (UNESP) Department of Biological Sciences Vanderbilt University Department of Molecular Biosciences The Wenner-Gren Institute (MBW) Stockholm University Department of Chemical and Biological Sciences Institute of Biosciences of Botucatu São Paulo State University (UNESP) Department of Biophysics and Pharmacology Institute of Biosciences of Botucatu São Paulo State University (UNESP) FAPESP: 2012/15589-7 FAPESP: 2013/06864-7

High-throughput sequencing for micro-RNAs (miRNAs) to obtain expression estimates is a central method of molecular biology. Surprisingly, there are a number of different approaches to converting sequencing output into micro-RNA counts. Each has their own strengths and biases that impact on the final data that can be obtained from a sequencing run. This chapter serves to make the reader aware of the trade-offs one must consider in analyzing small RNA sequencing data. It then compares two methods, miRge2.0 and the sRNAbench and the steps utilized to output data from their tools.

RNAcentral is a comprehensive database of non-coding RNA (ncRNA) sequences that provides a single access point to 44 RNA resources and >18 million ncRNA sequences from a wide range of organisms and RNA types. RNAcentral now also includes secondary (2D) structure information for >13 million sequences, making RNAcentral the world’s largest RNA 2D structure database. The 2D diagrams are displayed using R2DT, a new 2D structure visualization method that uses consistent, reproducible and recognizable layouts for related RNAs. The sequence similarity search has been updated with a faster interface featuring facets for filtering search results by RNA type, organism, source database or any keyword. This sequence search tool is available as a reusable web component, and has been integrated into several RNAcentral member databases, including Rfam, miRBase and snoDB. To allow for a more fine-grained assignment of RNA types and subtypes, all RNAcentral sequences have been annotated with Sequence Ontology terms. The RNAcentral database continues to grow and provide a central data resource for the RNA community. RNAcentral is freely available at https://rnacentral.org.

Aging is a key risk factor for chronic diseases of the elderly. MicroRNAs regulate post-transcriptional gene silencing through base-pair binding on their target mRNAs. We identified nonlinear changes in age-related microRNAs by analyzing whole blood from 1334 healthy individuals. We observed a larger influence of the age as compared to the sex and provide evidence for a shift to the 5’ mature form of miRNAs in healthy aging. The addition of 3059 diseased patients uncovered pan-disease and disease-specific alterations in aging profiles. Disease biomarker sets for all diseases were different between young and old patients. Computational deconvolution of whole-blood miRNAs into blood cell types suggests that cell intrinsic gene expression changes may impart greater significance than cell abundance changes to the whole blood miRNA profile. Altogether, these data provide a foundation for understanding the relationship between healthy aging and disease, and for the development of age-specific disease biomarkers., Aging is a key risk factor for chronic diseases of the elderly. Here the authors perform large-scale miRNA profiling of blood from individuals of a range of ages and show that common diseases alter the physiological age-related blood microRNA profile.

RNA sequencing data sets rapidly increase in quantity. For microRNAs (miRNAs), frequently dozens to hundreds of billion reads are generated per study. The quantification of annotated miRNAs and the prediction of new miRNAs are leading computational tasks. Now, the increased depth of coverage allows to gain deeper insights into the variability of miRNAs. The analysis of isoforms of miRNAs (isomiRs) is a trending topic, and a range of computational tools for the analysis of isomiRs has been developed. We provide an overview on 27 available computational solutions for the analysis of isomiRs. These include both stand-alone programs (17 tools) and web-based solutions (10 tools) and span a publication time range from 2010 to 2020. Seven of the tools were published in 2019 and 2020, confirming the rising importance of the topic. While most of the analyzed tools work for a broad range of organisms or are completely independent of a reference organism, several tools have been tailored for the analysis of human miRNA data or for plants. While 14 of the tools are general analysis tools of miRNAs, and isomiR analysis is one of their features, the remaining 13 tools have specifically been developed for isomiR analysis. A direct comparison on 20 deep sequencing data sets for selected tools provides insights into the heterogeneity of results. With our work, we provide users a comprehensive overview on the landscape of isomiR analysis tools and in that support the selection of the most appropriate tool for their respective research task.

Trematodes are widespread parasitic flatworms that significantly affect mankind either directly as human parasites, or indirectly via the infection of livestock and the related economic damage. The two most important trematode taxa are the blood flukes Schistosoma and the liver flukes Fasciola, but detection and differentiation of these parasites remains a challenge. Recently, microRNAs (miRNAs) were described from extracellular vesicles (EV) for both parasites secreted into respective hosts. These molecules have been proposed as mediators of parasite-host communication, and potential biomarkers for the detection of parasitic infections from host blood. Our aim here was to study similarities and differences in the miRNA complements of Schistosoma mansoni and Fasciola hepatica, EV-load in particular, to predict their targets and potential functions in the parasite-host interaction. We reanalyzed the known miRNA complements of S. mansoni and F. hepatica and found 16 and 4 previously overlooked, but deeply conserved miRNAs, respectively, further moving their complements closer together. We found distinct miRNA enrichment patterns in EVs both showing high levels of flatworm miRNAs with potential for the detection of an infection from blood. Two miRNAs of the protostome specific MIR-71 and MIR-277 families were highly expressed in EVs and could, therefore, have potential as biomarkers for trematode infection. Curiously, we identified nucleotide differences in the sequence of Mir-277-P2 between S. mansoni and F. hepatica that hold great promise for the distinction of both parasites. To test whether the EV-miRNAs of S. mansoni and F. hepatica could be modulating the expression of host genes, we predicted miRNA targets in 321 human and cattle messenger RNAs that overlapped between both hosts. Of several predicted targets, wnt signaling pathway genes stood out and their suppression likely leads to changes in the glucose concentration in host blood and the reduction of inflammatory and immune responses.

Although microRNAs (miRNA) are involved in all hallmarks of cancer, miRNA dysregulation in metastasis remains poorly understood and contradictory results have been published. The aim of this work was to identify miRNAs associated with metastatic progression of colorectal cancer (CRC). Novel and previously published next generation sequencing (NGS) datasets generated from 268 samples with primary (pCRC) and metastatic CRC (mCRC; liver, lung and peritoneal metastases) and tumor adjacent tissues were analyzed. Differential expression analysis was performed using a meticulous bioinformatics pipeline, including only bona fide miRNAs, utilizing miRNA-tailored quality control and processing, and applying a physiologically meaningful cut-off value (100 reads per million). The results were adjusted for host tissue background expression and samples from the different metastatic sites were independently analyzed. A metastatic signature containing five miRNAs up-regulated at multiple metastatic sites was identified (Mir-210_3p, Mir-191_5p, Mir-8-P1b_3p (mir-141-3p), Mir-1307_5p, and Mir-155_5p) along with a number of miRNAs that were differentially expressed at individual metastatic sites. Several of these have previously been implicated in metastasis through involvement in epithelial-to-mesenchymal transition and hypoxia, while other identified miRNAs represent novel findings. The identified differentially expressed miRNAs confirm known associations and contribute novel insights into miRNA involvement in the metastatic process. The use of open science practices facilitates reproducibility, and new datasets may easily be added to the publicly available pipeline to continuously improve the knowledge in the field. The identified set of miRNAs provides a reliable starting-point for further research into the role of miRNAs in metastatic progression.

BACKGROUND & AIMS: Chromosomal instability (CIN) is a carcinogenesis event that promotes metastasis and resistance to therapy by unclear mechanisms. Expression of the colon cancer-associated transcript 2 gene (CCAT2), which encodes a long noncoding RNA (lncRNA), associates with CIN, but little is known about how CCAT2 lncRNA regulates this cancer enabling characteristic.METHODS: We performed cytogenetic analysis of colorectal cancer (CRC) cell lines (HCT116, KM12C/SM, and HT29) overexpressing CCAT2 and colon organoids from C57BL/6N mice with the CCAT2 transgene and without (controls). CRC cells were also analyzed by immunofluorescence microscopy, gamma-H2AX, and senescence assays. CCAT2 transgene and control mice were given azoxymethane and dextran sulfate sodium to induce colon tumors. We performed gene expression array and mass spectrometry to detect downstream targets of CCAT2 lncRNA. We characterized interactions between CCAT2 with downstream proteins using MS2 pull-down, RNA immunoprecipitation, and selective 2'-hydroxyl acylation analyzed by primer extension analyses. Downstream proteins were overexpressed in CRC cells and analyzed for CIN. Gene expression levels were measured in CRC and non-tumor tissues from 5 cohorts, comprising more than 900 patients.RESULTS: High expression of CCAT2 induced CIN in CRC cell lines and increased resistance to 5-fluorouracil and oxaliplatin. Mice that expressed the CCAT2 transgene developed chromosome abnormalities, and colon organoids derived from crypt cells of these mice had a higher percentage of chromosome abnormalities compared with organoids from control mice. The transgenic mice given azoxymethane and dextran sulfate sodium developed more and larger colon polyps than control mice given these agents. Microarray analysis and mass spectrometry indicated that expression of CCAT2 increased expression of genes involved in ribosome biogenesis and protein synthesis. CCAT2 lncRNA interacted directly with and stabilized BOP1 ribosomal biogenesis factor (BOP1). CCAT2 also increased expression of MYC, which activated expression of BOP1. Overexpression of BOP1 in CRC cell lines resulted in chromosomal missegregation errors, and increased colony formation, and invasiveness, whereas BOP1 knockdown reduced viability. BOP1 promoted CIN by increasing the active form of aurora kinase B, which regulates chromosomal segregation. BOP1 was overexpressed in polyp tissues from CCAT2 transgenic mice compared with healthy tissue. CCAT2 lncRNA and BOP1 mRNA or protein were all increased in microsatellite stable tumors (characterized by CIN), but not in tumors with microsatellite instability compared with nontumor tissues. Increased levels of CCAT2 lncRNA and BOP1 mRNA correlated with each other and with shorter survival times of patients.CONCLUSIONS: We found that overexpression of CCAT2 in colon cells promotes CIN and carcinogenesis by stabilizing and inducing expression of BOP1 an activator of aurora kinase B. Strategies to target this pathway might be developed for treatment of patients with microsatellite stable colorectal tumors.

Background To date, mitochondrial genomes of more than one hundred flatworms (Platyhelminthes) have been sequenced. They show a high degree of similarity and a strong taxonomic bias towards parasitic lineages. The mitochondrial gene atp8 has not been confidently annotated in any flatworm sequenced to date. However, sampling of free-living flatworm lineages is incomplete. We addressed this by sequencing the mitochondrial genomes of the two small-bodied (about 1 mm in length) free-living flatworms Stenostomum sthenum and Macrostomum lignano as the first representatives of the earliest branching flatworm taxa Catenulida and Macrostomorpha respectively. Results We have used high-throughput DNA and RNA sequence data and PCR to establish the mitochondrial genome sequences and gene orders of S. sthenum and M. lignano. The mitochondrial genome of S. sthenum is 16,944 bp long and includes a 1,884 bp long inverted repeat region containing the complete sequences of nad3, rrnS, and nine tRNA genes. The model flatworm M. lignano has the smallest known mitochondrial genome among free-living flatworms, with a length of 14,193 bp. The mitochondrial genome of M. lignano lacks duplicated genes, however, tandem repeats were detected in a non-coding region. Mitochondrial gene order is poorly conserved in flatworms, only a single pair of adjacent ribosomal or protein-coding genes – nad4l-nad4 – was found in S. sthenum and M. lignano that also occurs in other published flatworm mitochondrial genomes. Unexpectedly, we unambiguously identified the full metazoan mitochondrial protein-coding gene complement including atp8 in S. sthenum and M. lignano. A subsequent search detected atp8 in all mitochondrial genomes of polyclad flatworms published to date, although the gene wasn’t previously annotated in these species. Conclusions Manual, but not automated genome annotation revealed the presence of atp8 in basally branching free-living flatworms, signifying both the importance of manual data curation and of diverse taxon sampling. We conclude that the loss of atp8 within flatworms is restricted to the parasitic taxon Neodermata. Electronic supplementary material The online version of this article (doi:10.1186/s12864-017-3807-2) contains supplementary material, which is available to authorized users.

A lack of knowledge of the cellular origin of microRNAs (miRNAs) has greatly confounded functional and biomarkers studies. Recently, three studies characterized miRNA expression patterns across >78 human cell types. These combined data expand our knowledge of miRNA expression localization and confirm that many miRNAs show cell-type specific expression patterns.

Ancient DNA sequencing is the key technology for paleogenomic studies and today a routine method in many laboratories. Recent analyses have shown that, under favoring conditions, also RNA can be sequenced from historical and even ancient samples. We have re-analyzed ancient RNA data from a Pleistocene canid and find - in addition to the previously described messenger RNA fragments - intact microRNAs, which are short transcripts with important gene regulatory functions. With an extraordinary age of 14,300 years, the canid microRNA profiles are the oldest ever reported. Despite their age, we show that the microRNA profiles are conclusive of taxonomic origin, tissue identity with organ- and cell-type specific signatures, and that they yield glimpses into gene regulatory activity and biological processes from the Pleistocene. In summary, we here show that straightforward microRNA analyses hold great promise for deeper insights into gene regulation in extinct animals.