Expanding the annotation of zebrafish microRNAs based on small RNA sequencing

Mapping the m1A, m5C, m6A and m7G methylation atlas in zebrafish brain under hypoxic conditions by MeRIP-seq

BACKGROUND

The epigenetic modifications play important regulatory roles in tissue development, maintenance of physiological functions and pathological process. RNA methylations, including newly identified m1A, m5C, m6A and m7G, are important epigenetic modifications. However, how these modifications are distributed in the transcriptome of vertebrate brains and whether their abundance is altered under pathological conditions are still poorly understood. In this study, we chose the model animal of zebrafish to conduct a systematic study to investigate the mRNA methylation atlas in the brain.

RESULTS

By performing unbiased analyses of the m1A, m5C, m6A and m7G methylation of mRNA, we found that within the whole brain transcriptome, with the increase of the gene expression levels, the overall level of each of these four modifications on the related genes was also progressively increased. Further bioinformatics analysis indicated that the zebrafish brain has an abundance of m1A modifications. In the hypoxia-treated zebrafish brains, the proportion of m1A is decreased, affecting the RNA splicing and zebrafish endogenous retroviruses.

CONCLUSIONS

Our study presents the first comprehensive atlas of m1A, m5C, m6A and m7G in the epitranscriptome of the zebrafish brain and reveals the distribution of these modifications in mRNA under hypoxic conditions. These data provide an invaluable resource for further research on the involvement of m1A, m5C, m6A and m7G in the regulation of miRNA and repeat elements in vertebrates, and provide new thoughts to study the brain hypoxic injury on the aspect of epitranscriptome.

FishmiRNA: An evolutionarily supported microRNA annotation and expression database for ray-finned fishes

microRNAs are important post-transcriptional regulators of gene expression involved in countless biological processes and are widely studied across metazoans. While miRNA research continues to grow, the large community of fish miRNA researchers lacks exhaustive resources consistent among species. To fill this gap, we developed FishmiRNA, an evolutionarily supported microRNA annotation and expression database for ray-finned fishes: www.fishmirna.org. The self-explanatory database contains detailed, manually-curated miRNA annotations with orthology relationships rigorously established by sequence similarity and conserved syntenies, and expression data provided for each detected mature miRNA. In just few clicks, users can download the annotation and expression database in several convenient formats either in its entirety or a subset. Simple filters and BLAST search options also permit the simultaneous exploration and visual comparison of expression data for up to any ten mature miRNAs across species and organs. FishmiRNA was specifically designed for ease of use to reach a wide audience.

Circulating miRNA repertoire as a biomarker of metabolic and reproductive states in rainbow trout

Background

Circulating miRNAs (c-miRNAs) are found in most, if not all, biological fluids and are becoming well-established non-invasive biomarkers of many human pathologies. However, their features in non-pathological contexts and whether their expression profiles reflect normal life history events have received little attention, especially in non-mammalian species. The aim of the present study was to investigate the potential of c-miRNAs to serve as biomarkers of reproductive and metabolic states in fish.

Results

The blood plasma was sampled throughout the reproductive cycle of female rainbow trout subjected to two different feeding regimes that triggered contrasting metabolic states. In addition, ovarian fluid was sampled at ovulation, and all samples were subjected to small RNA-seq analysis, leading to the establishment of a comprehensive miRNA repertoire (i.e., miRNAome) and enabling subsequent comparative analyses to a panel of RNA-seq libraries from a wide variety of tissues and organs. We showed that biological fluid miRNAomes are complex and encompass a high proportion of the overall rainbow trout miRNAome. While sharing a high proportion of common miRNAs, the blood plasma and ovarian fluid miRNAomes exhibited strong fluid-specific signatures. We further revealed that the blood plasma miRNAome significantly changed depending on metabolic and reproductive states. We subsequently identified three evolutionarily conserved muscle-specific miRNAs or myomiRs (miR-1-1/2-3p, miR-133a-1/2-3p, and miR-206-3p) that accumulated in the blood plasma in response to high feeding rates, making these myomiRs strong candidate biomarkers of active myogenesis. We also identified miR-202-5p as a candidate biomarker for reproductive success that could be used to predict ovulation and/or egg quality.

Conclusions

Together, these promising results reveal the high potential of c-miRNAs, including evolutionarily conserved myomiRs, as physiologically relevant biomarker candidates and pave the way for the use of c-miRNAs for non-invasive phenotyping in various fish species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01163-5.

Evolution after Whole-Genome Duplication: Teleost MicroRNAs

MicroRNAs (miRNAs) are important gene expression regulators implicated in many biological processes, but we lack a global understanding of how miRNA genes evolve and contribute to developmental canalization and phenotypic diversification. Whole-genome duplication events likely provide a substrate for species divergence and phenotypic change by increasing gene numbers and relaxing evolutionary pressures. To understand the consequences of genome duplication on miRNA evolution, we studied miRNA genes following the teleost genome duplication (TGD). Analysis of miRNA genes in four teleosts and in spotted gar, whose lineage diverged before the TGD, revealed that miRNA genes were retained in ohnologous pairs more frequently than protein-coding genes, and that gene losses occurred rapidly after the TGD. Genomic context influenced retention rates, with clustered miRNA genes retained more often than nonclustered miRNA genes and intergenic miRNA genes retained more frequently than intragenic miRNA genes, which often shared the evolutionary fate of their protein-coding host. Expression analyses revealed both conserved and divergent expression patterns across species in line with miRNA functions in phenotypic canalization and diversification, respectively. Finally, major strands of miRNA genes experienced stronger purifying selection, especially in their seeds and 3'-complementary regions, compared with minor strands, which nonetheless also displayed evolutionary features compatible with constrained function. This study provides the first genome-wide, multispecies analysis of the mechanisms influencing metazoan miRNA evolution after whole-genome duplication.

In Silico Target Prediction of Overexpressed microRNAs from LPS-Challenged Zebrafish (Danio rerio) Treated with the Novel Anti-Inflammatory Peptide TnP

miRNAs regulate gene expression post-transcriptionally in various processes, e.g., immunity, development, and diseases. Since their experimental analysis is complex, in silico target prediction is important for directing investigations. TnP is a candidate peptide for anti-inflammatory therapy, first discovered in the venom of Thalassophryne nattereri, which led to miRNAs overexpression in LPS-inflamed zebrafish post-treatment. This work aimed to predict miR-21, miR-122, miR-731, and miR-26 targets using overlapped results of DIANA microT-CDS and TargetScanFish software. This study described 513 miRNAs targets using highly specific thresholds. Using Gene Ontology over-representation analysis, we identified their main roles in regulating gene expression, neurogenesis, DNA-binding, transcription regulation, immune system process, and inflammatory response. miRNAs act in post-transcriptional regulation, but we revealed that their targets are strongly related to expression regulation at the transcriptional level, e.g., transcription factors proteins. A few predicted genes participated concomitantly in many biological processes and molecular functions, such as foxo3a, rbpjb, rxrbb, tyrobp, hes6, zic5, smad1, e2f7, and npas4a. Others were particularly involved in innate immunity regulation: il17a/f2, pik3r3b, and nlrc6. Together, these findings not only provide new insights into the miRNAs mode of action but also raise hope for TnP therapy and may direct future experimental investigations.

Altered non-coding RNA expression profile in F1 progeny 1 year after parental irradiation is linked to adverse effects in zebrafish

Gamma radiation produces DNA instability and impaired phenotype. Previously, we observed negative effects on phenotype, DNA methylation, and gene expression profiles, in offspring of zebrafish exposed to gamma radiation during gametogenesis. We hypothesize that previously observed effects are accompanied with changes in the expression profile of non-coding RNAs, inherited by next generations. Non-coding RNA expression profile was analysed in F₁ offspring (5.5 h post-fertilization) by high-throughput sequencing 1 year after parental irradiation (8.7 mGy/h, 5.2 Gy total dose). Using our previous F₁-γ genome-wide gene expression data (GSE98539), hundreds of mRNAs were predicted as targets of differentially expressed (DE) miRNAs, involved in pathways such as insulin receptor, NFkB and PTEN signalling, linking to apoptosis and cancer. snRNAs belonging to the five major spliceosomal snRNAs were down-regulated in the F₁-γ group, Indicating transcriptional and post-transcriptional alterations. In addition, DEpiRNA clusters were associated to 9 transposable elements (TEs) (LTR, LINE, and TIR) (p = 0.0024), probable as a response to the activation of these TEs. Moreover, the expression of the lincRNAs malat-1, and several others was altered in the offspring F₁, in concordance with previously observed phenotypical alterations. In conclusion, our results demonstrate diverse gamma radiation-induced alterations in the ncRNA profiles of F₁ offspring observable 1 year after parental irradiation.

microRNA expression variation as a potential molecular mechanism contributing to adaptation to hydrogen sulphide

microRNAs (miRNAs) are post-transcriptional regulators of gene expression and can play an important role in modulating organismal development and physiology in response to environmental stress. However, the role of miRNAs in mediating adaptation to diverse environments in natural study systems remains largely unexplored. Here, we characterized miRNAs and their expression in Poecilia mexicana, a species of small fish that inhabits both normal streams and extreme environments in the form of springs rich in toxic hydrogen sulphide (H₂ S). We found that P. mexicana has a similar number of miRNA genes as other teleosts. In addition, we identified a large population of mature miRNAs that were differentially expressed between locally adapted populations in contrasting habitats, indicating that miRNAs may contribute to P. mexicana adaptation to sulphidic environments. In silico identification of differentially expressed miRNA-mRNA pairs revealed, in the sulphidic environment, the downregulation of miRNAs predicted to target mRNAs involved in sulphide detoxification and cellular homeostasis, which are pathways essential for life in H₂ S-rich springs. In addition, we found that predicted targets of upregulated miRNAs act in the mitochondria (16.6% of predicted annotated targets), which is the main site of H₂ S toxicity and detoxification, possibly modulating mitochondrial function. Together, the differential regulation of miRNAs between these natural populations suggests that miRNAs may be involved in H₂ S adaptation by promoting functions needed for survival and reducing functions affected by H₂ S. This study lays the groundwork for further research to directly demonstrate the role of miRNAs in adaptation to H₂ S. Overall, this study provides a critical stepping-stone towards a comprehensive understanding of the regulatory mechanisms underlying the adaptive variation in gene expression in a natural system.

Identification and characterization of a specific 13-miRNA expression signature during follicle activation in the zebrafish ovary

Ovarian folliculogenesis is always of great interest in reproductive biology. However, the molecular mechanisms that control follicle development, particularly the early phase of follicle activation or recruitment, still remain poorly understood. In an attempt to decipher the gene networks and signaling pathways involved in such transition, we conducted a transcriptomic analysis (RNA-seq) on zebrafish primary growth (PG, stage I; inactive) and previtellogenic (PV, stage II; activated) follicles. A total of 118 unique microRNAs (miRNAs) (11 downregulated and 83 upregulated during PG/PV transition) and 56711 unique messenger RNAs (mRNAs) (1839 downregulated and 7243 upregulated during PG/PV transition) were identified. Real-time quantitative polymerase chain reaction analysis confirmed differential expression of 46 miRNAs from 66 candidates (66.67%). Among which, we chose to focus on 13 miRNAs (let-7a, -7b, -7c-5p, -7d-5p, -7h, -7i; miR-21, -23a-3p, -27c-3p, -107a-3p, -125b-5p, -145-3p, and -202-5p) that exhibited significant differential expression between PG and PV follicles (P ≤ 0.045*). With this 13-miRNA expression signature alone, PG follicles can be well differentiated from PV follicles by hierarchical clustering, suggesting their functional relevance during PG-to-PV transition. By overlaying predicted target genes and the differentially expressed mRNAs revealed by the RNA-seq analysis, especially those showing reciprocal miRNA-mRNA expression patterns, we shortlisted a panel of miRNA downstream targets for luciferase reporter validation. The reporter assay confirmed the interactions of let-7i:: atg4a (P = 0.01*), miR-202-5p::c23h20orf24 (P = 0.0004***), and miR-144-5p::ybx1 (P = 0.003**), implicating these potential miRNA-mRNA gene pairs in follicle activation during folliculogenesis. Our transcriptomic data analyses suggest that miRNA-mediated post-transcriptional control may represent an important mechanism underlying follicle activation.

miRNA analysis with Prost! reveals evolutionary conservation of organ-enriched expression and post-transcriptional modifications in three-spined stickleback and zebrafish

MicroRNAs (miRNAs) can have organ-specific expression and functions; they can originate from dedicated miRNA genes, from non-canonical miRNA genes, or from mirror-miRNA genes and can also experience post-transcriptional variation. It remains unclear, however, which mechanisms of miRNA production or modification are organ-specific and the extent of their evolutionary conservation. To address these issues, we developed the software Prost! (PRocessing Of Short Transcripts), which, among other features, helps quantify mature miRNAs, accounts for post-transcriptional processing, such as nucleotide editing, and identifies mirror-miRNAs. Here, we applied Prost! to annotate and analyze miRNAs in three-spined stickleback (Gasterosteus aculeatus), a model fish for evolutionary biology reported to have a miRNome larger than most teleost fish. Zebrafish (Danio rerio), a distantly related teleost with a well-known miRNome, served as comparator. Our results provided evidence for the existence of 286 miRNA genes and 382 unique mature miRNAs (excluding mir430 gene duplicates and the vaultRNA-derived mir733), which doesn't represent a miRNAome larger than other teleost miRNomes. In addition, small RNA sequencing data from brain, heart, testis, and ovary in both stickleback and zebrafish identified suites of mature miRNAs that display organ-specific enrichment, many of which are evolutionarily-conserved in the brain and heart in both species. These data also supported the hypothesis that evolutionarily-conserved, organ-specific mechanisms may regulate post-transcriptional variations in miRNA sequence. In both stickleback and zebrafish, miR2188-5p was edited frequently with similar nucleotide changes in the seed sequence with organ specific editing rates, highest in the brain. In summary, Prost! is a new tool to identify and understand small RNAs, to help clarify a species' miRNA biology as shown here for an important model for the evolution of developmental mechanisms, and to provide insight into organ-enriched expression and the evolutionary conservation of miRNA post-transcriptional modifications.

Antarctic blackfin icefish genome reveals adaptations to extreme environments

Icefishes (suborder Notothenioidei; family Channichthyidae) are the only vertebrates that lack functional haemoglobin genes and red blood cells. Here, we report a high-quality genome assembly and linkage map for the Antarctic blackfin icefish Chaenocephalus aceratus, highlighting evolved genomic features for its unique physiology. Phylogenomic analysis revealed that Antarctic fish of the teleost suborder Notothenioidei, including icefishes, diverged from the stickleback lineage about 77 million years ago and subsequently evolved cold-adapted phenotypes as the Southern Ocean cooled to sub-zero temperatures. Our results show that genes involved in protection from ice damage, including genes encoding antifreeze glycoprotein and zona pellucida proteins, are highly expanded in the icefish genome. Furthermore, genes that encode enzymes that help to control cellular redox state, including members of the sod3 and nqo1 gene families, are expanded, probably as evolutionary adaptations to the relatively high concentration of oxygen dissolved in cold Antarctic waters. In contrast, some crucial regulators of circadian homeostasis (cry and per genes) are absent from the icefish genome, suggesting compromised control of biological rhythms in the polar light environment. The availability of the icefish genome sequence will accelerate our understanding of adaptation to extreme Antarctic environments.

MiRNAs: Biology, Biogenesis, their Web-based Tools, and Databases

INTRODUCTION

MicroRNAs (miRNAs), which are evolutionarily conserved, and endogenous non-coding RNAs, participate in the post-transcriptional regulation of eukaryotic genes. The biogenesis of miRNAs occurs in the nucleus. Then, in the cytoplasm, they are assembled along with some proteins in a ribonucleoprotein complex called RISC. miRNA component of the RISC complex binds to the complementary sequence of mRNA target depending on the degree of complementarity, and leads to mRNA degradation and/or inhibition of protein synthesis. miRNAs have been found in eukaryotes and some viruses play a role in development, metabolism, cell proliferation, growth, differentiation, and death.

OBJECTIVE

A large number of miRNAs and their targets were identified by different experimental techniques and computational approaches. The principal aim of this paper is to gather information about some miRNA databases and web-based tools for better and quicker access to relevant data.

RESULTS

Accordingly, in this paper, we collected and introduced miRNA databases and some webbased tools that have been developed by various research groups. We have categorized them into different classes including databases for viral miRNAs, and plant miRNAs, miRNAs in human beings, mice and other vertebrates, miRNAs related to human diseases, and target prediction, and miRNA expression. Also, we have presented relevant statistical information about these databases.

Identification of Novel MicroRNAs and Characterization of MicroRNA Expression Profiles in Zebrafish Ovarian Follicular Cells

MicroRNAs (miRNAs) are small non-coding RNAs that regulate gene expression primarily at the post-transcriptional levels and thereby play important roles in regulating many physiological and developmental processes. Oocyte maturation in fish is induced by hormones produced from the hypothalamus, pituitary, and ovary. Gonadotropin-releasing hormone (GnRH) stimulates the secretion of luteinizing hormone (LH), which in turn, induces the secretion of maturation-inducing hormone (MIH) from the ovary. It is documented that small early vitellogenic (or stage IIIa) follicles are unable to undergo oocyte maturation whereas oocytes in mid- to late vitellogenic (stage IIIb) follicles can be induced by LH and MIH to become mature. To determine whether miRNAs may be involved in the growth and acquisition of maturational competency of ovarian follicles, we determined the miRNA expression profiles in follicular cells collected from stage IIIa and IIIb follicles using next-generation sequencing. It was found that miRNAs are abundantly expressed in the follicular cells from both stages IIIa and IIIb follicles. Furthermore, bioinformatics analysis revealed the presence of 214 known, 31 conserved novel and 44 novel miRNAs in zebrafish vitellogenic ovarian follicular cells. Most mature miRNAs in follicular cells were found to be in the length of 22 nucleotides. Differential expression analysis revealed that 11 miRNAs were significantly up-regulated, and 13 miRNAs were significantly down-regulated in the stage IIIb follicular cells as compared with stage IIIa follicular cells. The expression of four of the significantly regulated miRNAs, dre-miR-22a-3p, dre-miR-16a, dre-miR-181a-3p, and dre-miR-29a, was validated by real-time PCR. Finally, gene enrichment and pathway analyses of the predicted targets of the significantly regulated miRNAs supported the involvement of several key signaling pathways in regulating ovarian function, including oocyte maturation. Taken together, this study identifies novel zebrafish miRNAs and characterizes miRNA expression profiles in somatic cells within the zebrafish ovarian follicles. The differential expression of miRNAs between stage IIIa and IIIb follicular cells suggests that these miRNAs are important regulators of zebrafish ovarian follicle development and/or oocyte maturation.

Omics and cytokine discovery in fish: Presenting the Yellowtail kingfish (Seriola lalandi) as a case study

A continued programme of research is essential to overcome production bottlenecks in any aquacultured fish species. Since the introduction of genetic and molecular techniques, the quality of immune research undertaken in fish has greatly improved. Thousands of species specific cytokine genes have been discovered, which can be used to conduct more sensitive studies to understand how fish physiology is affected by aquaculture environments or disease. Newly available transcriptomic technologies, make it increasingly easier to study the immunogenetics of farmed species for which little data exists. This paper reviews how the application of transcriptomic procedures such as RNA Sequencing (RNA-Seq) can advance fish research. As a case study, we present some preliminary findings using RNA-Seq to identify cytokine related genes in Seriola lalandi. These will allow in-depth investigations to understand the immune responses of these fish in response to environmental change or disease and help in the development of therapeutic approaches.

Genomic conservation of erythropoietic microRNAs (erythromiRs) in white-blooded Antarctic icefish

White-blooded Antarctic crocodile icefish are the only vertebrates known to lack functional hemoglobin genes and red blood cells throughout their lives. We do not yet know, however, whether extinction of hemoglobin genes preceded loss of red blood cells or vice versa, nor whether erythropoiesis regulators disappeared along with hemoglobin genes in this erythrocyte-null clade. Several microRNAs, which we here call erythromiRs, are expressed primarily in developing red blood cells in zebrafish, mouse, and humans. Abrogating some erythromiRs, like mir144 and mir451a, leads to profound anemia, demonstrating a functional role in erythropoiesis. Here, we tested two not mutually exclusive hypotheses: 1) that the loss of one or more erythromiR genes extinguished the erythropoietic program of icefish and/or led to the loss of globin gene expression through pseudogenization; and 2) that some erythromiR genes were secondarily lost after the loss of functional hemoglobin and red blood cells in icefish. We explored small RNA transcriptomes generated from the hematopoietic kidney marrow of four Antarctic notothenioids: two red-blooded species (bullhead notothen Notothenia coriiceps and emerald notothen Trematomus bernacchii) and two white-blooded icefish (blackfin icefish Chaenocephalus aceratus and hooknose icefish Chionodraco hamatus). The N. coriiceps genome assembly anchored analyses. Results showed that, like the two red-blooded species, the blackfin icefish genome possessed and the marrow expressed all known erythromiRs. This result indicates that loss of hemoglobin and red blood cells in icefish was not caused by loss of known erythromiR genes. Furthermore, expression of only one erythromiR, mir96, appears to have been lost after the loss of red blood cells and hemoglobin-expression was not detected in the erythropoietic organ of hooknose icefish but was present in blackfin icefish. All other erythromiRs investigated, including mir144 and mir451a, were expressed by all four species and thus are present in the genomes of at least the two white-blooded icefish. Our results rule out the hypothesis that genomic loss of any known erythromiRs extinguished erythropoiesis in icefish, and suggest that after the loss of red blood cells, few erythromiRs experienced secondary loss. Results suggest that functions independent of erythropoiesis maintained erythromiRs, thereby highlighting the evolutionary resilience of miRNA genes in vertebrate genomes.

The spotted gar genome illuminates vertebrate evolution and facilitates human-teleost comparisons

To connect human biology to fish biomedical models, we sequenced the genome of spotted gar (Lepisosteus oculatus), whose lineage diverged from teleosts before teleost genome duplication (TGD). The slowly evolving gar genome has conserved in content and size many entire chromosomes from bony vertebrate ancestors. Gar bridges teleosts to tetrapods by illuminating the evolution of immunity, mineralization and development (mediated, for example, by Hox, ParaHox and microRNA genes). Numerous conserved noncoding elements (CNEs; often cis regulatory) undetectable in direct human-teleost comparisons become apparent using gar: functional studies uncovered conserved roles for such cryptic CNEs, facilitating annotation of sequences identified in human genome-wide association studies. Transcriptomic analyses showed that the sums of expression domains and expression levels for duplicated teleost genes often approximate the patterns and levels of expression for gar genes, consistent with subfunctionalization. The gar genome provides a resource for understanding evolution after genome duplication, the origin of vertebrate genomes and the function of human regulatory sequences.

In silico genome wide mining of conserved and novel miRNAs in the brain and pineal gland of Danio rerio using small RNA sequencing data

MicroRNAs (miRNAs) are small, non-coding RNA molecules that bind to the mRNA of the target genes and regulate the expression of the gene at the post-transcriptional level. Zebrafish is an economically important freshwater fish species globally considered as a good predictive model for studying human diseases and development. The present study focused on uncovering known as well as novel miRNAs, target prediction of the novel miRNAs and the differential expression of the known miRNA using the small RNA sequencing data of the brain and pineal gland (dark and light treatments) obtained from NCBI SRA. A total of 165, 151 and 145 known zebrafish miRNAs were found in the brain, pineal gland (dark treatment) and pineal gland (light treatment), respectively. Chromosomes 4 and 5 of zebrafish reference assembly GRCz10 were found to contain maximum number of miR genes. The miR-181a and miR-182 were found to be highly expressed in terms of number of reads in the brain and pineal gland, respectively. Other ncRNAs, such as tRNA, rRNA and snoRNA, were curated against Rfam. Using GRCz10 as reference, the subsequent bioinformatic analyses identified 25, 19 and 9 novel miRNAs from the brain, pineal gland (dark treatment) and pineal gland (light treatment), respectively. Targets of the novel miRNAs were identified, based on sequence complementarity between miRNAs and mRNA, by searching for antisense hits in the 3'-UTR of reference RNA sequences of the zebrafish. The discovery of novel miRNAs and their targets in the zebrafish genome can be a valuable scientific resource for further functional studies not only in zebrafish but also in other economically important fishes.

miRNA Nomenclature: A View Incorporating Genetic Origins, Biosynthetic Pathways, and Sequence Variants

High-throughput sequencing of miRNAs has revealed the diversity and variability of mature and functional short noncoding RNAs, including their genomic origins, biogenesis pathways, sequence variability, and newly identified products such as miRNA-offset RNAs (moRs). Here we review known cases of alternative mature miRNA-like RNA fragments and propose a revised definition of miRNAs to encompass this diversity. We then review nomenclature guidelines for miRNAs and propose to extend nomenclature conventions to align with those for protein-coding genes established by international consortia. Finally, we suggest a system to encompass the full complexity of sequence variations (i.e., isomiRs) in the analysis of small RNA sequencing experiments.

A new model army: Emerging fish models to study the genomics of vertebrate Evo-Devo

Many fields of biology--including vertebrate Evo-Devo research--are facing an explosion of genomic and transcriptomic sequence information and a multitude of fish species are now swimming in this "genomic tsunami." Here, we first give an overview of recent developments in sequencing fish genomes and transcriptomes that identify properties of fish genomes requiring particular attention and propose strategies to overcome common challenges in fish genomics. We suggest that the generation of chromosome-level genome assemblies--for which we introduce the term "chromonome"--should be a key component of genomic investigations in fish because they enable large-scale conserved synteny analyses that inform orthology detection, a process critical for connectivity of genomes. Orthology calls in vertebrates, especially in teleost fish, are complicated by divergent evolution of gene repertoires and functions following two rounds of genome duplication in the ancestor of vertebrates and a third round at the base of teleost fish. Second, using examples of spotted gar, basal teleosts, zebrafish-related cyprinids, cavefish, livebearers, icefish, and lobefin fish, we illustrate how next generation sequencing technologies liberate emerging fish systems from genomic ignorance and transform them into a new model army to answer longstanding questions on the genomic and developmental basis of their biodiversity. Finally, we discuss recent progress in the genetic toolbox for the major fish models for functional analysis, zebrafish, and medaka, that can be transferred to many other fish species to study in vivo the functional effect of evolutionary genomic change as Evo-Devo research enters the postgenomic era.