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.

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.