FishmiRNA: An Evolutionarily Supported MicroRNA Annotation and Expression Database for Ray-Finned Fishes.

MicroRNAs (miRNAs) are important post-transcriptional regulators of gene expression involved in countless biological processes and are widely studied across metazoans. Although 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 miRNA 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.

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.

A fish with no sex: gonadal and adrenal functions partition between zebrafish NR5A1 co-orthologs.

People with NR5A1 mutations experience testicular dysgenesis, ovotestes, or adrenal insufficiency, but we do not completely understand the origin of this phenotypic diversity. NR5A1 is expressed in gonadal soma precursor cells before expression of the sex-determining gene SRY. Many fish have two co-orthologs of NR5A1 that likely partitioned ancestral gene subfunctions between them. To explore ancestral roles of NR5A1, we knocked out nr5a1a and nr5a1b in zebrafish. Single-cell RNA-seq identified nr5a1a-expressing cells that co-expressed genes for steroid biosynthesis and the chemokine receptor Cxcl12a in 1-day postfertilization (dpf) embryos, as does the mammalian adrenal-gonadal (interrenal-gonadal) primordium. In 2dpf embryos, nr5a1a was expressed stronger in the interrenal-gonadal primordium than in the early hypothalamus but nr5a1b showed the reverse. Adult Leydig cells expressed both ohnologs and granulosa cells expressed nr5a1a stronger than nr5a1b. Mutants for nr5a1a lacked the interrenal, formed incompletely differentiated testes, had no Leydig cells, and grew far larger than normal fish. Mutants for nr5a1b formed a disorganized interrenal and their gonads completely disappeared. All homozygous mutant genotypes lacked secondary sex characteristics, including male breeding tubercles and female sex papillae, and had exceedingly low levels of estradiol, 11-ketotestosterone, and cortisol. RNA-seq showed that at 21dpf, some animals were developing as females and others were not, independent of nr5a1 genotype. By 35dpf, all mutant genotypes greatly under-expressed ovary-biased genes. Because adult nr5a1a mutants form gonads but lack an interrenal and conversely, adult nr5a1b mutants lack a gonad but have an interrenal, the adrenal, and gonadal functions of the ancestral nr5a1 gene partitioned between ohnologs after the teleost genome duplication, likely owing to reciprocal loss of ancestral tissue-specific regulatory elements. Identifying such elements could provide hints to otherwise unexplained cases of Differences in Sex Development.

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.

A Hormone That Lost Its Receptor: Anti-Müllerian Hormone (AMH) in Zebrafish Gonad Development and Sex Determination.

Fetal mammalian testes secrete Anti-Müllerian hormone (Amh), which inhibits female reproductive tract (Müllerian duct) development. Amh also derives from mature mammalian ovarian follicles, which marks oocyte reserve and characterizes polycystic ovarian syndrome. Zebrafish (Danio rerio) lacks Müllerian ducts and the Amh receptor gene amhr2 but, curiously, retains amh To discover the roles of Amh in the absence of Müllerian ducts and the ancestral receptor gene, we made amh null alleles in zebrafish. Results showed that normal amh prevents female-biased sex ratios. Adult male amh mutants had enormous testes, half of which contained immature oocytes, demonstrating that Amh regulates male germ cell accumulation and inhibits oocyte development or survival. Mutant males formed sperm ducts and some produced a few offspring. Young female mutants laid a few fertile eggs, so they also had functional sex ducts. Older amh mutants accumulated nonvitellogenic follicles in exceedingly large but sterile ovaries, showing that Amh helps control ovarian follicle maturation and proliferation. RNA-sequencing data partitioned juveniles at 21 days postfertilization (dpf) into two groups that each contained mutant and wild-type fish. Group21-1 upregulated ovary genes compared to Group21-2, which were likely developing as males. By 35 dpf, transcriptomes distinguished males from females and, within each sex, mutants from wild types. In adult mutants, ovaries greatly underexpressed granulosa and theca genes, and testes underexpressed Leydig cell genes. These results show that ancestral Amh functions included development of the gonadal soma in ovaries and testes and regulation of gamete proliferation and maturation. A major gap in our understanding is the identity of the gene encoding a zebrafish Amh receptor; we show here that the loss of amhr2 is associated with the breakpoint of a chromosome rearrangement shared among cyprinid fishes.

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.

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.

Expanding the annotation of zebrafish microRNAs based on small RNA sequencing.

MicroRNAs (miRs) are short non-coding RNAs that fine-tune the regulation of gene expression to coordinate a wide range of biological processes. Because of their role in the regulation of gene expression, miRs are essential players in development by acting on cell fate determination and progression towards cell differentiation and are increasingly relevant to human health and disease. Although the zebrafish Danio rerio is a major model for studies of development, genetics, physiology, evolution, and human biology, the annotation of zebrafish miR-producing genes remains limited. In the present work, we report deep sequencing data of zebrafish small RNAs from brain, heart, testis, and ovary. Results provide evidence for the expression of 56 un-annotated mir genes and 248 un-annotated mature strands, increasing the number of zebrafish mir genes over those already deposited in miRBase by 16% and the number of mature sequences by 63%. We also describe the existence of three pairs of mirror-mir genes and two mirtron genes, genetic features previously undescribed in non-mammalian vertebrates. This report provides information that substantially increases our knowledge of the zebrafish miRNome and will benefit the entire miR community.

Minimal Information for Neural Electromagnetic Ontologies (MINEMO): A standards-compliant method for analysis and integration of event-related potentials (ERP) data.

We present MINEMO (Minimal Information for Neural ElectroMagnetic Ontologies), a checklist for the description of event-related potentials (ERP) studies. MINEMO extends MINI (Minimal Information for Neuroscience Investigations)to the ERP domain. Checklist terms are explicated in NEMO, a formal ontology that is designed to support ERP data sharing and integration. MINEMO is also linked to an ERP database and web application (the NEMO portal). Users upload their data and enter MINEMO information through the portal. The database then stores these entries in RDF (Resource Description Framework), along with summary metrics, i.e., spatial and temporal metadata. Together these spatial, temporal, and functional metadata provide a complete description of ERP data and the context in which these data were acquired. The RDF files then serve as inputs to ontology-based labeling and meta-analysis. Our ultimate goal is to represent ERPs using a rich semantic structure, so results can be queried at multiple levels, to stimulate novel hypotheses and to promote a high-level, integrative account of ERP results across diverse study methods and paradigms.