Hox genes control homocercal caudal fin development and evolution

Ancient bony fishes had heterocercal tails, like modern sharks and sturgeons, with asymmetric caudal fins and a vertebral column extending into an elongated upper lobe. Teleost fishes, in contrast, developed a homocercal tail characterized by two separate equal-sized fin lobes and the body axis not extending into the caudal fin. A similar heterocercal-to-homocercal transition occurs during teleost ontogeny, although the underlying genetic and developmental mechanisms for either transition remain unresolved. Here, we investigated the role of hox13 genes in caudal fin formation as these genes control posterior identity in animals. Analysis of expression profiles of zebrafish hox13 paralogs and phenotypes of CRISPR/Cas9-induced mutants showed that double hoxb13a and hoxc13a mutants fail to form a caudal fin. Furthermore, single mutants display heterocercal-like morphologies not seen since Mesozoic fossil teleosteomorphs. Relaxation of functional constraints after the teleost genome duplication may have allowed hox13 duplicates to neo- or subfunctionalize, ultimately contributing to the evolution of a homocercal tail in teleost fishes.

Cold-Driven Hemoglobin Evolution in Antarctic Notothenioid Fishes Prior to Hemoglobin Gene Loss in White-Blooded Icefishes

Expression of multiple hemoglobin isoforms with differing physiochemical properties likely helps species adapt to different environmental and physiological conditions. Antarctic notothenioid fishes inhabit the icy Southern Ocean and display fewer hemoglobin isoforms, each with less affinity for oxygen than temperate relatives. Reduced hemoglobin multiplicity was proposed to result from relaxed selective pressure in the cold, thermally stable, and highly oxygenated Antarctic waters. These conditions also permitted the survival and diversification of white-blooded icefishes, the only vertebrates living without hemoglobin. To understand hemoglobin evolution during adaptation to freezing water, we analyzed hemoglobin genes from 36 notothenioid genome assemblies. Results showed that adaptation to frigid conditions shaped hemoglobin gene evolution by episodic diversifying selection concomitant with cold adaptation and by pervasive evolution in Antarctic notothenioids compared to temperate relatives, likely a continuing adaptation to Antarctic conditions. Analysis of hemoglobin gene expression in adult hematopoietic organs in various temperate and Antarctic species further revealed a switch in hemoglobin gene expression underlying hemoglobin multiplicity reduction in Antarctic fish, leading to a single hemoglobin isoform in adult plunderfishes and dragonfishes, the sister groups to icefishes. The predicted high hemoglobin multiplicity in Antarctic fish embryos based on transcriptomic data, however, raises questions about the molecular bases and physiological implications of diverse hemoglobin isoforms in embryos compared to adults. This analysis supports the hypothesis that the last common icefish ancestor was vulnerable to detrimental mutations affecting the single ancestral expressed alpha- and beta-globin gene pair, potentially predisposing their subsequent loss.

Genomics of cold adaptations in the Antarctic notothenioid fish radiation

Numerous novel adaptations characterise the radiation of notothenioids, the dominant fish group in the freezing seas of the Southern Ocean. To improve understanding of the evolution of this iconic fish group, here we generate and analyse new genome assemblies for 24 species covering all major subgroups of the radiation, including five long-read assemblies. We present a new estimate for the onset of the radiation at 10.7 million years ago, based on a time-calibrated phylogeny derived from genome-wide sequence data. We identify a two-fold variation in genome size, driven by expansion of multiple transposable element families, and use the long-read data to reconstruct two evolutionarily important, highly repetitive gene family loci. First, we present the most complete reconstruction to date of the antifreeze glycoprotein gene family, whose emergence enabled survival in sub-zero temperatures, showing the expansion of the antifreeze gene locus from the ancestral to the derived state. Second, we trace the loss of haemoglobin genes in icefishes, the only vertebrates lacking functional haemoglobins, through complete reconstruction of the two haemoglobin gene clusters across notothenioid families. Both the haemoglobin and antifreeze genomic loci are characterised by multiple transposon expansions that may have driven the evolutionary history of these genes.

Pathways to polar adaptation in fishes revealed by long-read sequencing

Long-read sequencing is driving a new reality for genome science in which highly contiguous assemblies can be produced efficiently with modest resources. Genome assemblies from long-read sequences are particularly exciting for understanding the evolution of complex genomic regions that are often difficult to assemble. In this study, we utilized long-read sequencing data to generate a high-quality genome assembly for an Antarctic eelpout, Ophthalmolycus amberensis, the first for the globally distributed family Zoarcidae. We used this assembly to understand how O. amberensis has adapted to the harsh Southern Ocean and compared it to another group of Antarctic fishes: the notothenioids. We showed that selection has largely acted on different targets in eelpouts relative to notothenioids. However, we did find some overlap; in both groups, genes involved in membrane structure, thermal tolerance and vision have evidence of positive selection. We found evidence for historical shifts of transposable element activity in O. amberensis and other polar fishes, perhaps reflecting a response to environmental change. We were specifically interested in the evolution of two complex genomic loci known to underlie key adaptations to polar seas: haemoglobin and antifreeze proteins (AFPs). We observed unique evolution of the haemoglobin MN cluster in eelpouts and related fishes in the suborder Zoarcoidei relative to other Perciformes. For AFPs, we identified the first species in the suborder with no evidence of afpIII sequences (Cebidichthys violaceus) in the genomic region where they are found in all other Zoarcoidei, potentially reflecting a lineage-specific loss of this cluster. Beyond polar fishes, our results highlight the power of long-read sequencing to understand genome evolution.

Genome structures resolve the early diversification of teleost fishes

Accurate species phylogenies are a prerequisite for all evolutionary research. Teleosts are the largest and most diversified group of extant vertebrates, but relationships among their three oldest extant lineages remain unresolved. On the basis of seven high-quality new genome assemblies in Elopomorpha (tarpons, eels), we revisited the topology of the deepest branches of the teleost phylogeny using independent gene sequence and chromosomal rearrangement phylogenomic approaches. These analyses converged to a single scenario that unambiguously places the Elopomorpha and Osteoglossomorpha (arapaima, elephantnose fish) in a monophyletic sister group to all other teleosts, i.e., the Clupeocephala lineage (zebrafish, medaka). This finding resolves more than 50 years of controversy on the evolutionary relationships of these lineages and highlights the power of combining different levels of genome-wide information to solve complex phylogenies.

Coordinated patterning of zebrafish caudal fin symmetry by a central and two peripheral organizers

BACKGROUND

Caudal fin symmetry characterizes teleosts and likely contributes to their evolutionary success. However, the coordinated development and patterning of skeletal elements establishing external symmetry remains incompletely understood. We explore the spatiotemporal emergence of caudal skeletal elements in zebrafish to consider evolutionary and developmental origins of caudal fin symmetry.

RESULTS

Transgenic reporters and skeletal staining reveal that the hypural diastema-defining gap between hypurals 2 and 3 forms early and separates progenitors of two plates of connective tissue. Two sets of central principal rays (CPRs) synchronously, sequentially, and symmetrically emerge around the diastema. The two dorsal- and ventral-most rays (peripheral principal rays, PPRs) arise independently and earlier than adjacent CPRs. Muscle and tendon markers reveal that different muscles attach to CPR and PPR sets.

CONCLUSIONS

We propose that caudal fin symmetry originates from a central organizer that establishes the hypural diastema and bidirectionally patterns surrounding tissue into two plates of connective tissue and two mirrored sets of CPRs. Further, two peripheral organizers unidirectionally specify PPRs, forming a symmetric "composite" fin derived from three fields. Distinct CPR and PPR ontogenies may represent developmental modules conferring ray identities, muscle connections, and biomechanical properties. Our model contextualizes mechanistic studies of teleost fin morphological variation.

A parasite outbreak in notothenioid fish in an Antarctic fjord

Climate changes can promote disease outbreaks, but their nature and potential impacts in remote areas have received little attention. In a hot spot of biodiversity on the West Antarctic Peninsula, which faces among the fastest changing climates on Earth, we captured specimens of two notothenioid fish species affected by large skin tumors at an incidence never before observed in the Southern Ocean. Molecular and histopathological analyses revealed that X-cell parasitic alveolates, members of a genus we call Notoxcellia, are the etiological agent of these tumors. Parasite-specific molecular probes showed that xenomas remained within the skin but largely outgrew host cells in the dermis. We further observed that tumors induced neovascularization in underlying tissue and detrimentally affected host growth and condition. Although many knowledge gaps persist about X-cell disease, including its mode of transmission and life cycle, these findings reveal potentially active biotic threats to vulnerable Antarctic ecosystems.

An ancient truncated duplication of the anti‐Mullerian hormone receptor type 2 gene is a potential conserved master sex determinant in the Pangasiidae catfish family

The evolution of sex determination (SD) in teleosts is amazingly dynamic, as reflected by the variety of different master sex-determining genes identified. Pangasiids are economically important catfishes in South Asian countries, but little is known about their SD system. Here, we generated novel genomic resources for 12 Pangasiids and characterized their SD system. Based on a Pangasianodon hypophthalmus chromosome-scale genome assembly, we identified an anti-Müllerian hormone receptor type Ⅱ gene (amhr2) duplication, which was further characterized as being sex-linked in males and expressed only in testes. These results point to a Y chromosome male-specific duplication (amhr2by) of the autosomal amhr2a. Sequence annotation revealed that the P. hypophthalmus Amhr2by is truncated in its N-terminal domain, lacking the cysteine-rich extracellular part of the receptor that is crucial for ligand binding, suggesting a potential route for its neofunctionalization. Reference-guided assembly of 11 additional Pangasiids, along with sex-linkage studies, revealed that this truncated amhr2by duplication is a male-specific conserved gene in Pangasiids. Reconstructions of the amhr2 phylogeny suggested that amhr2by arose from an ancient duplication/insertion event at the root of the Siluroidei radiation that is dated to ~100 million years ago. Together these results bring multiple lines of evidence supporting that amhr2by is an ancient and conserved master sex-determining gene in Pangasiids, a finding that highlights the recurrent use of the transforming growth factor β pathway, which is often used for the recruitment of teleost master SD genes, and provides another empirical case towards firther understanding of dynamics of SD systems.

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.

Discovery of novel fish papillomaviruses: From the Antarctic to the commercial fish market

Fish papillomaviruses form a newly discovered group broadly recognized as the Secondpapillomavirinae subfamily. This study expands the documented genomes of the fish papillomaviruses from six to 16, including one from the Antarctic emerald notothen, seven from commercial market fishes, one from data mining of sea bream sequence data, and one from a western gull cloacal swab that is likely diet derived. The genomes of secondpapillomaviruses are ∼6 kilobasepairs (kb), which is substantially smaller than the ∼8 kb of terrestrial vertebrate papillomaviruses. Each genome encodes a clear homolog of the four canonical papillomavirus genes, E1, E2, L1, and L2. In addition, we identified open reading frames (ORFs) with short linear peptide motifs reminiscent of E6/E7 oncoproteins. Fish papillomaviruses are extremely diverse and phylogenetically distant from other papillomaviruses suggesting a model in which terrestrial vertebrate-infecting papillomaviruses arose after an evolutionary bottleneck event, possibly during the water-to-land transition.

Advancing human disease research with fish evolutionary mutant models

Model organism research is essential to understand disease mechanisms. However, laboratory-induced genetic models can lack genetic variation and often fail to mimic the spectrum of disease severity. Evolutionary mutant models (EMMs) are species with evolved phenotypes that mimic human disease. EMMs complement traditional laboratory models by providing unique avenues to study gene-by-environment interactions, modular mutations in noncoding regions, and their evolved compensations. EMMs have improved our understanding of complex diseases, including cancer, diabetes, and aging, and illuminated mechanisms in many organs. Rapid advancements of sequencing and genome-editing technologies have catapulted the utility of EMMs, particularly in fish. Fish are the most diverse group of vertebrates, exhibiting a kaleidoscope of specialized phenotypes, many that would be pathogenic in humans but are adaptive in the species' specialized habitat. Importantly, evolved compensations can suggest avenues for novel disease therapies. This review summarizes current research using fish EMMs to advance our understanding of human disease.

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.

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.

Bone microstructure and bone mineral density are not systemically different in Antarctic icefishes and related Antarctic notothenioids

Ancestors of the Antarctic icefishes (family Channichthyidae) were benthic and had no swim bladder, making it energetically expensive to rise from the ocean floor. To exploit the water column, benthopelagic icefishes were hypothesized to have evolved a skeleton with “reduced bone,” which gross anatomical data supported. Here, we tested the hypothesis that changes to icefish bones also occurred below the level of gross anatomy. Histology and micro‐CT imaging of representative craniofacial bones (i.e., ceratohyal, frontal, dentary, and articular) of extant Antarctic fish species specifically evaluated two features that might cause the appearance of “reduced bone”: bone microstructure (e.g., bone volume fraction and structure linear density) and bone mineral density (BMD, or mass of mineral per volume of bone). Measures of bone microstructure were not consistently different in bones from the icefishes Chaenocephalus aceratus and Champsocephalus gunnari, compared to the related benthic notothenioids Notothenia coriiceps and Gobionotothen gibberifrons. Some quantitative measures, such as bone volume fraction and structure linear density, were significantly increased in some icefish bones compared to homologous bones of non‐icefish. However, such differences were rare, and no microstructural measures were consistently different in icefishes across all bones and species analyzed. Furthermore, BMD was similar among homologous bones of icefish and non‐icefish Antarctic notothenioids. In summary, “reduced bone” in icefishes was not due to systemic changes in bone microstructure or BMD, raising the prospect that “reduced bone” in icefish occurs only at the gross anatomic level (i.e., smaller or fewer bones). Given that icefishes exhibit delayed skeletal development compared to non‐icefish Antarctic fishes, combining these phenotypic data with genomic data might clarify genetic changes driving skeletal heterochrony.

Productivity and Change in Fish and Squid in the Southern Ocean

Southern Ocean ecosystems are globally important and vulnerable to global drivers of change, yet they remain challenging to study. Fish and squid make up a significant portion of the biomass within the Southern Ocean, filling key roles in food webs from forage to mid-trophic species and top predators. They comprise a diverse array of species uniquely adapted to the extreme habitats of the region. Adaptations such as antifreeze glycoproteins, lipid-retention, extended larval phases, delayed senescence, and energy-conserving life strategies equip Antarctic fish and squid to withstand the dark winters and yearlong subzero temperatures experienced in much of the Southern Ocean. In addition to krill exploitation, the comparatively high commercial value of Antarctic fish, particularly the lucrative toothfish, drives fisheries interests, which has included illegal fishing. Uncertainty about the population dynamics of target species and ecosystem structure and function more broadly has necessitated a precautionary, ecosystem approach to managing these stocks and enabling the recovery of depleted species. Fisheries currently remain the major local driver of change in Southern Ocean fish productivity, but global climate change presents an even greater challenge to assessing future changes. Parts of the Southern Ocean are experiencing ocean-warming, such as the West Antarctic Peninsula, while other areas, such as the Ross Sea shelf, have undergone cooling in recent years. These trends are expected to result in a redistribution of species based on their tolerances to different temperature regimes. Climate variability may impair the migratory response of these species to environmental change, while imposing increased pressures on recruitment. Fisheries and climate change, coupled with related local and global drivers such as pollution and sea ice change, have the potential to produce synergistic impacts that compound the risks to Antarctic fish and squid species. The uncertainty surrounding how different species will respond to these challenges, given their varying life histories, environmental dependencies, and resiliencies, necessitates regular assessment to inform conservation and management decisions. Urgent attention is needed to determine whether the current management strategies are suitably precautionary to achieve conservation objectives in light of the impending changes to the ecosystem.

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.

RADSex: A computational workflow to study sex determination using restriction site-associated DNA sequencing data

The study of sex determination and sex chromosome organization in nonmodel species has long been technically challenging, but new sequencing methodologies now enable precise and high-throughput identification of sex-specific genomic sequences. In particular, restriction site-associated DNA sequencing (RAD-Seq) is being extensively applied to explore sex determination systems in many plant and animal species. However, software specifically designed to search for and visualize sex-biased markers using RAD-Seq data is lacking. Here, we present RADSex, a computational analysis workflow designed to study the genetic basis of sex determination using RAD-Seq data. RADSex is simple to use, requires few computational resources, makes no prior assumptions about the type of sex-determination system or structure of the sex locus, and offers convenient visualization through a dedicated R package. To demonstrate the functionality of RADSex, we re-analysed a published data set of Japanese medaka, Oryzias latipes, where we uncovered a previously unknown Y chromosome polymorphism. We then used RADSex to analyse new RAD-Seq data sets from 15 fish species spanning multiple taxonomic orders. We identified the sex determination system and sex-specific markers in six of these species, five of which had no known sex-markers prior to this study. We show that RADSex greatly facilitates the study of sex determination systems in nonmodel species thanks to its speed of analyses, low resource usage, ease of application and visualization options. Furthermore, our analysis of new data sets from 15 species provides new insights on sex determination in fish.

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.

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.

The genome sequence of the channel bull blenny, Cottoperca gobio (Günther, 1861)

We present a genome assembly for Cottoperca gobio (channel bull blenny, (Günther, 1861)); Chordata; Actinopterygii (ray-finned fishes), a temperate water outgroup for Antarctic Notothenioids. The size of the genome assembly is 609 megabases, with the majority of the assembly scaffolded into 24 chromosomal pseudomolecules. Gene annotation on Ensembl of this assembly has identified 21,662 coding genes.

Retinal oxygen supply shaped the functional evolution of the vertebrate eye

The retina has a very high energy demand but lacks an internal blood supply in most vertebrates. Here we explore the hypothesis that oxygen diffusion limited the evolution of retinal morphology by reconstructing the evolution of retinal thickness and the various mechanisms for retinal oxygen supply, including capillarization and acid-induced haemoglobin oxygen unloading. We show that a common ancestor of bony fishes likely had a thin retina without additional retinal oxygen supply mechanisms and that three different types of retinal capillaries were gained and lost independently multiple times during the radiation of vertebrates, and that these were invariably associated with parallel changes in retinal thickness. Since retinal thickness confers multiple advantages to vision, we propose that insufficient retinal oxygen supply constrained the functional evolution of the eye in early vertebrates, and that recurrent origins of additional retinal oxygen supply mechanisms facilitated the phenotypic evolution of improved functional eye morphology.

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.

Developmental tuning of mineralization drives morphological diversity of gill cover bones in sculpins and their relatives

The role of osteoblast placement in skeletal morphological variation is relatively well understood, but alternative developmental mechanisms affecting bone shape remain largely unknown. Specifically, very little attention has been paid to variation in later mineralization stages of intramembranous ossification as a driver of morphological diversity. We discover the occurrence of specific, sometimes large, regions of nonmineralized osteoid within bones that also contain mineralized tissue. We show through a variety of histological, molecular, and tomographic tests that this “extended” osteoid material is most likely nonmineralized bone matrix. This tissue type is a significant determinant of gill cover bone shape in the teleostean suborder Cottoidei. We demonstrate repeated evolution of extended osteoid in Cottoidei through ancestral state reconstruction and test for an association between extended osteoid variation and habitat differences among species. Through measurement of extended osteoid at various stages of gill cover development in species across the phylogeny, we gain insight into possible evolutionary developmental origins of the trait. We conclude that this fine‐tuned developmental regulation of bone matrix mineralization reflects heterochrony at multiple biological levels and is a novel mechanism for the evolution of diversity in skeletal morphology. This research lays the groundwork for a new model in which to study bone mineralization and evolutionary developmental processes, particularly as they may relate to adaptation during a prominent evolutionary radiation of 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.

Adaptation of Proteins to the Cold in Antarctic Fish: A Role for Methionine?

The evolution of antifreeze glycoproteins has enabled notothenioid fish to flourish in the freezing waters of the Southern Ocean. Whereas successful at the biodiversity level to life in the cold, paradoxically at the cellular level these stenothermal animals have problems producing, folding, and degrading proteins at their ambient temperatures of -1.86 °C. In this first multi-species transcriptome comparison of the amino acid composition of notothenioid proteins with temperate teleost proteins, we show that, unlike psychrophilic bacteria, Antarctic fish provide little evidence for the mass alteration of protein amino acid composition to enhance protein folding and reduce protein denaturation in the cold. The exception was the significant overrepresentation of positions where leucine in temperate fish proteins was replaced by methionine in the notothenioid orthologues. We hypothesize that these extra methionines have been preferentially assimilated into the genome to act as redox sensors in the highly oxygenated waters of the Southern Ocean. This redox hypothesis is supported by analyses of notothenioids showing enrichment of genes associated with responses to environmental stress, particularly reactive oxygen species. So overall, although notothenioid fish show cold-associated problems with protein homeostasis, they may have modified only a selected number of biochemical pathways to work efficiently below 0 °C. Even a slight warming of the Southern Ocean might disrupt the critical functions of this handful of key pathways with considerable impacts for the functioning of this ecosystem in the future.

Evolution of caudal fin ray development and caudal fin hypural diastema complex in spotted gar, teleosts, and other neopterygian fishes

BACKGROUND

The caudal fin of actinopterygians transitioned from a heterocercal dorsoventrally asymmetrical fin to a homocercal externally symmetrical fin in teleosts through poorly understood evolutionary developmental mechanisms. We studied the caudal skeleton of major living actinopterygian lineages, including polypteriformes, acipenseriformes, Holostei (gars and bowfin), and teleosts, compared with reports of extinct neopterygians and basal teleosteans. We focused on the hypural diastema complex, which includes (1) a gap between hypurals 2 and 3, that (2) separates two plates of connective tissue at (3) the branching of caudal vasculature; these features had been considered as a shared, derived trait of teleosts, a synapomorphy.

RESULTS

These studies revealed that gars and teleosts share all three features of the hypural diastema complex. Absence of a complex with these features from bowfin, fossil Holostei, and stem Teleostei argues in favor of repetitive, independent emergence in several neopterygian and basal Teleostei lineages, or less likely, many independent losses. We further observed that, in gars and teleosts, the earliest developing lepidotrichia align with the horizontal adult body axis, thus participating in external symmetry.

CONCLUSIONS

These results suggest that the hypural diastema complex in teleosts and gars represents a homoplasy among neopterygians and that it emerged repeatedly by parallel evolution due to shared inherited underlying genetic and developmental programs (latent homology). Because the hypural diastema complex exists in gars with heterocercal tails, this complex is independent of homocercality. Developmental Dynamics 247:832-853, 2018. © 2018 Wiley Periodicals, Inc.

Skeletal development in the heterocercal caudal fin of spotted gar (lepisosteus oculatus) and other lepisosteiformes

BACKGROUND

The caudal fin of actinopterygians experienced substantial morphological changes during evolution. In basal actinopterygians, the caudal fin skeleton supports an asymmetrical heterocercal caudal fin, while most teleosts have a symmetrical homocercal caudal fin. The transition from the ancestral heterocercal form to the derived homocercal caudal fin remains poorly understood. Few developmental studies provide an understanding of derived and ancestral characters among basal actinopterygians. To fill this gap, we examined the development of the caudal fin of spotted gar Lepisosteus oculatus, one of only eight living species of Holostei, the sister group to the teleosts.

RESULTS

Our observations of animals from fertilization to more than a year old provide the most detailed description of the development of caudal fin skeletal elements in any Holostean species. We observed two different types of distal caudal radials replacing two transient plates of connective tissue, identifying two hypaxial ensembles separated by a space between hypurals 2 and 3. These features have not been described in any gar species, but can be observed in other gar species, and thus represent anatomical structures common to lepisosteiformes.

CONCLUSIONS

The present work highlights the power and importance of ontogenic studies and provides bases for future evolutionary and morphological investigations on actinopterygians fins. Developmental Dynamics 247:724-740, 2018. © 2018 Wiley Periodicals, Inc.

[Reproductive behavior of the green birdmouth wrasse Gomphosus caeruleus on a Reunion Island reef: Mode of reproduction, environmental factors and reproductive strategy alternation]

The green birdmouth wrasse Gomphosus caeruleus is present all year round on the coral reefs of Reunion Island (Indian Ocean). A group of individuals was followed on one of these reefs with the objective of studying the reproduction mode of the species, the influence of environmental factors, and social behaviors on the control of reproduction. Our observations revealed that G. caeruleus is, like many Labridae, a protogynous hermaphrodite species, probably diandric, that the reproduction of G. caeruleus is, like in other reef fish species, influenced by the lunar cycle with a peak of reproductive activity during waxing gibbous phase, and that G. caeruleus displays social behavior leading to alternating haremic mating system on a single territory and lek-like mating systems without aggressions between males. These observations enhanced our knowledge of the reproduction of Labridae and reef species.

Gonadal soma controls ovarian follicle proliferation through Gsdf in zebrafish

BACKGROUND

Aberrant signaling between germ cells and somatic cells can lead to reproductive disease and depends on diffusible signals, including transforming growth factor-beta (TGFB) -family proteins. The TGFB-family protein Gsdf (gonadal soma derived factor) controls sex determination in some fish and is a candidate for mediating germ cell/soma signaling.

RESULTS

Zebrafish expressed gsdf in somatic cells of bipotential gonads and expression continued in ovarian granulosa cells and testicular Sertoli cells. Homozygous gsdf knockout mutants delayed leaving the bipotential gonad state, but then became a male or a female. Mutant females ovulated a few oocytes, then became sterile, accumulating immature follicles. Female mutants stored excess lipid and down-regulated aromatase, gata4, insulin receptor, estrogen receptor, and genes for lipid metabolism, vitellogenin, and steroid biosynthesis. Mutant females contained less estrogen and more androgen than wild-types. Mutant males were fertile. Genomic analysis suggests that Gsdf, Bmp15, and Gdf9, originated as paralogs in vertebrate genome duplication events.

CONCLUSIONS

In zebrafish, gsdf regulates ovarian follicle maturation and expression of genes for steroid biosynthesis, obesity, diabetes, and female fertility, leading to ovarian and extra-ovarian phenotypes that mimic human polycystic ovarian syndrome (PCOS), suggesting a role for a related TGFB signaling molecule in the etiology of PCOS. Developmental Dynamics 246:925-945, 2017. © 2017 Wiley Periodicals, Inc.

The advantage of channeling nucleotides for very processive functions

Nucleoside triphosphate (NTP)s, like ATP (adenosine 5'-triphosphate) and GTP (guanosine 5'-triphosphate), have long been considered sufficiently concentrated and diffusible to fuel all cellular ATPases (adenosine triphosphatases) and GTPases (guanosine triphosphatases) in an energetically healthy cell without becoming limiting for function. However, increasing evidence for the importance of local ATP and GTP pools, synthesised in close proximity to ATP- or GTP-consuming reactions, has fundamentally challenged our view of energy metabolism. It has become evident that cellular energy metabolism occurs in many specialised 'microcompartments', where energy in the form of NTPs is transferred preferentially from NTP-generating modules directly to NTP-consuming modules. Such energy channeling occurs when diffusion through the cytosol is limited, where these modules are physically close and, in particular, if the NTP-consuming reaction has a very high turnover, i.e. is very processive. Here, we summarise the evidence for these conclusions and describe new insights into the physiological importance and molecular mechanisms of energy channeling gained from recent studies. In particular, we describe the role of glycolytic enzymes for axonal vesicle transport and nucleoside diphosphate kinases for the functions of dynamins and dynamin-related GTPases.

The advantage of channeling nucleotides for very processive functions

Nucleoside triphosphate (NTP)s, like ATP (adenosine 5'-triphosphate) and GTP (guanosine 5'-triphosphate), have long been considered sufficiently concentrated and diffusible to fuel all cellular ATPases (adenosine triphosphatases) and GTPases (guanosine triphosphatases) in an energetically healthy cell without becoming limiting for function. However, increasing evidence for the importance of local ATP and GTP pools, synthesised in close proximity to ATP- or GTP-consuming reactions, has fundamentally challenged our view of energy metabolism. It has become evident that cellular energy metabolism occurs in many specialised 'microcompartments', where energy in the form of NTPs is transferred preferentially from NTP-generating modules directly to NTP-consuming modules. Such energy channeling occurs when diffusion through the cytosol is limited, where these modules are physically close and, in particular, if the NTP-consuming reaction has a very high turnover, i.e. is very processive. Here, we summarise the evidence for these conclusions and describe new insights into the physiological importance and molecular mechanisms of energy channeling gained from recent studies. In particular, we describe the role of glycolytic enzymes for axonal vesicle transport and nucleoside diphosphate kinases for the functions of dynamins and dynamin-related GTPases.

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.

Embryogenesis and early skeletogenesis in the antarctic bullhead notothen, Notothenia coriiceps

BACKGROUND

Environmental temperature influences rates of embryonic development, but a detailed staging series for vertebrate embryos developing in the subzero cold of Antarctic waters is not yet available from fertilization to hatching. Given projected warming of the Southern Ocean, it is imperative to establish a baseline to evaluate potential effects of changing climate on fish developmental dynamics.

RESULTS

We studied the Bullhead notothen (Notothenia coriiceps), a notothenioid fish inhabiting waters between -1.9 and +2 °C. In vitro fertilization produced embryos that progressed through cleavage, epiboly, gastrulation, segmentation, organogenesis, and hatching. We compared morphogenesis spatially and temporally to Zebrafish and medaka. Experimental animals hatched after about 6 months to early larval stages. To help understand skeletogenesis, we analyzed late embryos for expression of sox9 and runx2, which regulate chondrogenesis, osteogenesis, and eye development. Results revealed that, despite their prolonged developmental time course, N. coriiceps embryos developed similarly to those of other teleosts with large yolk cells.

CONCLUSIONS

Our studies set the stage for future molecular analyses of development in these extremophile fish. Results provide a foundation for understanding the impact of ocean warming on embryonic development and larval recruitment of notothenioid fish, which are key factors in the marine trophic system. Developmental Dynamics 245:1066-1080, 2016. © 2016 Wiley Periodicals, Inc.

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.

X-Linked Retinitis Pigmentosa 2 Is a Novel Maternal-Effect Gene Required for Left-Right Asymmetry in Zebrafish

Retinitis pigmentosa 2 (RP2) gene is responsible for up to 20% of X-linked retinitis pigmentosa, a severe heterogeneous genetic disorder resulting in progressive retinal degeneration in humans. In vertebrates, several bodies of evidence have clearly established the role of Rp2 protein in cilia genesis and/or function. Unexpectedly, some observations in zebrafish have suggested the oocyte-predominant expression of the rp2 gene, a typical feature of maternal-effect genes. In the present study, we investigate the maternal inheritance of rp2 gene products in zebrafish eggs in order to address whether rp2 could be a novel maternal-effect gene required for normal development. Although both rp2 mRNA and corresponding protein are expressed during oogenesis, rp2 mRNA is maternally inherited, in contrast to Rp2 protein. A knockdown of the protein transcribed from both rp2 maternal and zygotic mRNA results in delayed epiboly and severe developmental defects, including eye malformations, that were not observed when only the protein from zygotic origin was knocked down. Moreover, the knockdown of maternal and zygotic Rp2 revealed a high incidence of left-right asymmetry establishment defects compared to only zygotic knockdown. Here we show that rp2 is a novel maternal-effect gene exclusively expressed in oocytes within the zebrafish ovary and demonstrate that maternal rp2 mRNA is essential for successful embryonic development and thus contributes to egg developmental competence. Our observations also reveal that Rp2 protein translated from maternal mRNA is important to allow normal heart loop formation, thus providing evidence of a direct maternal contribution to left-right asymmetry establishment.

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.

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.

Maternally inherited npm2 mRNA is crucial for egg developmental competence in zebrafish

The molecular mechanisms underlying and determining egg developmental competence remain poorly understood in vertebrates. Nucleoplasmin (Npm2) is one of the few known maternal effect genes in mammals, but this maternal effect has never been demonstrated in nonmammalian species. A link between developmental competence and the abundance of npm2 maternal mRNA in the egg was previously established using a teleost fish model for egg quality. The importance of maternal npm2 mRNA for egg developmental competence remains unknown in any vertebrate species. In the present study, we aimed to characterize the contribution of npm2 maternal mRNA to early developmental success in zebrafish using a knockdown strategy. We report here the oocyte-specific expression of npm2 and maternal inheritance of npm2 mRNA in zebrafish eggs. The knockdown of the protein translated from this maternal mRNA results in developmental arrest before the onset of epiboly and subsequent embryonic death, a phenotype also observed in embryos lacking zygotic transcription. Npm2 knockdown also results in impaired transcription of the first-wave zygotic genes. Our results show that npm2 is also a maternal effect gene in a nonmammalian vertebrate species and that maternally inherited npm2 mRNA is crucial for egg developmental competence. We also show that de novo protein synthesis from npm2 maternal mRNA is critical for developmental success beyond the blastula stage and required for zygotic genome activation. Finally, our results suggest that npm2 maternal mRNA is an important molecular factor of egg quality in fish and possibly in all vertebrates.

Evolution of the miR199-214 cluster and vertebrate skeletal development

MicroRNA (miRs) are short non-coding RNAs that fine-tune the regulation of gene expression to coordinate a wide range of biological processes. MicroRNAs are transcribed from miR genes and primary miR transcripts are processed to approximately 22 nucleotide single strand mature forms that function as repressors of transcript translation when bound to the 3'UTR of protein coding transcripts in association with the RISC. Because of their role in the regulation of gene expression, miRs are essential players in development by acting on cell fate determination and progression toward cell differentiation. The miR199 and miR214 genes occupy an intronic cluster located on the opposite strand of the Dynamin3 gene. These miRNAs play major roles in a broad variety of developmental processes and diseases, including skeletal development and several types of cancer. In the work reported here, we first deciphered the origin of the miR199 and miR214 families by following evolution of miR paralogs and their host Dynamin paralogs. We then examined the expression patterns of miR199 and miR214 in developing zebrafish embryos and demonstrated their regulation through a common primary transcript. Results suggest an evolutionarily conserved regulation across vertebrate lineages. Our expression study showed predominant expression patterns for both miR in tissues surrounding developing craniofacial skeletal elements consistent with expression data in mouse and human, thus indicating a conserved role of miR199 and miR214 in vertebrate skeletogenesis.

The Nme gene family in fish

The Nme gene family, also known as Nm23 or NDPK, is a very ancient gene family that can be found in all kingdoms of life. In the late eighties, a gene of the Nme family, NME1, was identified as the first metastatic suppressor gene, resulting in a major interest for this family. Due to the complexity of the family, the need for a unified and evolutionary-supported gene nomenclature was recently stressed by the scientific community. Based on a complete evolutionary history study of the gene family in metazoans and vertebrates, a unified nomenclature was recently proposed and accepted by gene nomenclature consortia. In addition to its well-documented role in tumor metastasis, members of the Nme family are also involved in a wide variety of cellular and physiological processes. Available data in non-mammalian species remain, however, scarce with the noticeable exception of Drosophila in which a major role in development was reported. In fish, very few studies have specifically investigated the role of nme genes. Several transcriptomic and proteomic studies have, however, revealed the expression of nme genes in various fish organs and tissues, in mature oocytes, and during embryonic development. Altogether, interest for the Nme gene family in fish is growing and new functions/roles in fish biology are expected to be discovered in the forthcoming years. Here, we briefly review the current knowledge of the Nme family in fish.

The NME gene family in zebrafish oogenesis and early development

After the recent report of the expression of several nme genes in the zebrafish gonads, the present study aimed at further analyzing the expression of nme genes in the ovary with special attention for the nme transcripts that are maternally inherited and could thus participate in the determination of oocyte developmental competence. The expression levels of all groups I and II nme genes were characterized by QPCR in a panel of zebrafish tissues. The nme genes exhibiting an ovarian expression were subsequently monitored throughout oogenesis and early development, and their expression sites characterized using in situ hybridization. Here, we show that nme2b1, nme3, nme4, and nme6 are highly expressed in the ovary and present in the zebrafish oocyte throughout oogenesis. While the four transcripts are maternally inherited, nme3 and nme6 display a typical maternal profile and are detected in the zebrafish early embryo. In contrast to nme3, nme6, abundance exhibits a sharp decrease during early embryogenesis. After zygotic genome activation, we observed an increased expression of nme2b1, nme2b2, nme3, and nme6. The present study provides a comprehensive overview of the expression of nme family members during zebrafish oogenesis and early development. In addition, the maternal origin of two nme transcripts in the early embryo is reported here for the first time in any vertebrate species. Together, our observations suggest an important role of the nme family in oocyte and embryo development in vertebrates.

Nme gene family evolutionary history reveals pre-metazoan origins and high conservation between humans and the sea anemone, Nematostella vectensis

Background

The Nme gene family is involved in multiple physiological and pathological processes such as cellular differentiation, development, metastatic dissemination, and cilia functions. Despite the known importance of Nme genes and their use as clinical markers of tumor aggressiveness, the associated cellular mechanisms remain poorly understood. Over the last 20 years, several non-vertebrate model species have been used to investigate Nme functions. However, the evolutionary history of the family remains poorly understood outside the vertebrate lineage. The aim of the study was thus to elucidate the evolutionary history of the Nme gene family in Metazoans.

Methodology/Principal Findings

Using a total of 21 eukaryote species including 14 metazoans, the evolutionary history of Nme genes was reconstructed in the metazoan lineage. We demonstrated that the complexity of the Nme gene family, initially thought to be restricted to chordates, was also shared by the metazoan ancestor. We also provide evidence suggesting that the complexity of the family is mainly a eukaryotic innovation, with the exception of Nme8 that is likely to be a choanoflagellate/metazoan innovation. Highly conserved gene structure, genomic linkage, and protein domains were identified among metazoans, some features being also conserved in eukaryotes. When considering the entire Nme family, the starlet sea anemone is the studied metazoan species exhibiting the most conserved gene and protein sequence features with humans. In addition, we were able to show that most of the proteins known to interact with human NME proteins were also found in starlet sea anemone.

Conclusion/Significance

Together, our observations further support the association of Nme genes with key cellular functions that have been conserved throughout metazoan evolution. Future investigations of evolutionarily conserved Nme gene functions using the starlet sea anemone could shed new light on a wide variety of key developmental and cellular processes.

Nme protein family evolutionary history, a vertebrate perspective

BACKGROUND

The Nme family, previously known as Nm23 or NDPK, is involved in various molecular processes including tumor metastasis and some members of the family, but not all, exhibit a Nucleoside Diphosphate Kinase (NDPK) activity. Ten genes are known in humans, in which some members have been extensively studied. In non-mammalian species, the Nme protein family has received, in contrast, far less attention. The picture of the vertebrate Nme family remains thus incomplete and orthology relationships with mammalian counterparts were only partially characterized. The present study therefore aimed at characterizing the Nme gene repertoire in vertebrates with special interest for teleosts, and providing a comprehensive overview of the Nme gene family evolutionary history in vertebrates.

RESULTS

In the present study, we present the evolutionary history of the Nme family in vertebrates and characterize the gene family repertoire for the first time in several non-mammalian species. Our observations show that vertebrate Nme genes can be separated in two evolutionary distinct groups. Nme1, Nme2, Nme3, and Nme4 belong to Group I while vertebrate Nme5, Nme6, Nme7, Nme8, and Nme9 belong to Group II. The position of Nme10 is in contrast more debatable due to its very specific evolutionary history. The present study clearly indicates that Nme5, Nme6, Nme7, and Nme8 originate from duplication events that occurred before the chordate radiation. In contrast, Nme genes of the Group I have a very different evolutionary history as our results suggest that they all arise from a common gene present in the chordate ancestor. In addition, expression patterns of all zebrafish nme transcripts were studied in a broad range of tissues by quantitative PCR and discussed in the light of the function of their mammalian counterparts.

CONCLUSION

This work offers an evolutionary framework that will pave the way for future studies on vertebrate Nme proteins and provides a unified vertebrate Nme nomenclature that is consistent with the nomenclature in use in mammals. Based on protein structure and expression data, we also provide new insight into molecular functions of Nme proteins among vertebrates and raise intriguing questions on the roles of Nme proteins in gonads.