Since the description of zebrafish (Danio rerio) in 1822, the identity of its closest living relative has been unclear. To address this problem, we sequenced the exomes of 10 species in genus Danio, using the closely related Devario aequipinnatus as outgroup, to infer relationships across the 25 chromosomes of the zebrafish genome. The majority of relationships within Danio were remarkably consistent across all chromosomes. Relationships of chromosome segments, however, depended systematically upon their genomic location within zebrafish chromosomes. Regions near chromosome centers identified Danio kyathit and/or Danio aesculapii as the closest relative of zebrafish, while segments near chromosome ends supported only D. aesculapii as the zebrafish sister species. Genome-wide comparisons of derived character states revealed that danio relationships are inconsistent with a simple bifurcating species history but support an ancient hybrid origin of the D. rerio lineage by homoploid hybrid speciation. We also found evidence of more recent gene flow limited to the high recombination ends of chromosomes and several megabases of chromosome 20 with a history distinct from the rest of the genome. Additional insights gained from incorporating genome structure into a phylogenomic study demonstrate the utility of such an approach for future studies in other taxa. The multiple genomic histories of species in the genus Danio have important implications for comparative studies in these morphologically varied and beautiful species and for our understanding of the hybrid evolutionary history of zebrafish.

The genetics of sex determination varies across taxa, sometimes even within a species. Major domesticated strains of zebrafish (Danio rerio), including AB and TU, lack a strong genetic sex determining locus, but strains more recently derived from nature, like Nadia (NA), possess a ZZ male/ZW female chromosomal sex-determination system. AB fish pass through a juvenile ovary stage, forming oocytes that survive in fish that become females but die in fish that become males. To understand mechanisms of gonad development in NA zebrafish, we studied histology and single cell transcriptomics in developing ZZ and ZW fish. ZW fish developed oocytes by 22 days post-fertilization (dpf) but ZZ fish directly formed testes, avoiding a juvenile ovary phase. Gonads of some ZW and WW fish, however, developed oocytes that died as the gonad became a testis, mimicking AB fish, suggesting that the gynogenetically derived AB strain is chromosomally WW. Single-cell RNA-seq of 19dpf gonads showed similar cell types in ZZ and ZW fish, including germ cells, precursors of gonadal support cells, steroidogenic cells, interstitial/stromal cells, and immune cells, consistent with a bipotential juvenile gonad. In contrast, scRNA-seq of 30dpf gonads revealed that cells in ZZ gonads had transcriptomes characteristic of testicular Sertoli, Leydig, and germ cells while ZW gonads had granulosa cells, theca cells, and developing oocytes. Hematopoietic and vascular cells were similar in both sex genotypes. These results show that juvenile NA zebrafish initially develop a bipotential gonad; that a factor on the NA W chromosome, or fewer than two Z chromosomes, is essential to initiate oocyte development; and without the W factor, or with two Z doses, NA gonads develop directly into testes without passing through the juvenile ovary stage. Sex determination in AB and TU strains mimics NA ZW and WW zebrafish, suggesting loss of the Z chromosome during domestication. Genetic analysis of the NA strain will facilitate our understanding of the evolution of sex determination mechanisms.

The vertebrate sodium–iodide symporter (NIS or SLC5A5) transports iodide into the thyroid follicular cells that synthesize thyroid hormone. The SLC5A protein family includes transporters of vitamins, minerals, and nutrients. Disruption of SLC5A5 function by perchlorate, a pervasive environmental contaminant, leads to human pathologies, especially hypothyroidism. Perchlorate also disrupts the sexual development of model animals, including threespine stickleback ( Gasterosteus aculeatus ) and zebrafish ( Danio rerio ), but the mechanism of action is unknown. To test the hypothesis that SLC5A5 paralogs are expressed in tissues necessary for the development of reproductive organs, and therefore are plausible candidates to mediate the effects of perchlorate on sexual development, we first investigated the evolutionary history of Slc5a paralogs to better understand potential functional trajectories of the gene family. We identified two clades of slc5a paralogs with respect to an outgroup of sodium/choline cotransporters ( slc5a7 ); these clades are the NIS clade of sodium/iodide and lactate cotransporters ( slc5a5 , slc5a6 , slc5a8 , slc5a8 , and slc5a12 ) and the SGLT clade of sodium/glucose cotransporters ( slc5a1 , slc5a2 , slc5a3 , slc5a4 , slc5a10 , and slc5a11 ). We also characterized expression patterns of slc5a genes during development. Stickleback embryos and early larvae expressed NIS clade genes in connective tissue, cartilage, teeth, and thyroid. Stickleback males and females expressed slc5a5 and its paralogs in gonads. Single‐cell transcriptomics (scRNA‐seq) on zebrafish sex‐genotyped gonads revealed that NIS clade‐expressing cells included germ cells ( slc5a5 , slc5a6a , and slc5a6b ) and gonadal soma cells ( slc5a8l ). These results are consistent with the hypothesis that perchlorate exerts its effects on sexual development by interacting with slc5a5 or its paralogs in reproductive tissues. These findings show novel expression domains of slc5 genes in stickleback and zebrafish, which suggest similar functions across vertebrates including humans, and provide candidates to mediate the effects of perchlorate on sexual development.

When considering relationships between genotype and phenotype we frequently ignore the fact that the genome of a typical animal, notably including that of a fish and a human, harbors a huge amount of foreign DNA. Some of it, including the DNA of “autonomous” transposable elements, can spontaneously mobilize to occupy new chromosomal sites and take on new functions, presenting a challenge to the host organism and also possibly introducing new fuel for evolutionary change. Transposable elements are useful for introducing transgenes, integrating them into host genomes with high efficiency. Transgenesis has become very widespread in biological research, and in our society at large. This year the governments of both Canada and the United States have approved the first use of ‘genetically engineered’ animals in food production, Atlantic salmon, Salmo salar . With the recent advent of amazing gene-editing technology, there is no doubt that the transgene industry will grow explosively in the coming years. The biology of transgenes needs to be included in our understanding of the genome. It is in this spirit that we have investigated an unexpected and novel phenotypic effect of the chromosomally integrated transgene fli1a-F-hsp70l:Gal4VP16 . We examine larval fras1 mutant zebrafish ( Danio rerio ). Gal4VP16 is a potent transcriptional activator, and already well known for toxicity and mediating unusual transcriptional effects. In the presence of the transgene, phenotypes in the neural crest-derived craniofacial skeleton, notably fusions and shape changes associated with loss of function fras1 mutations, are made more severe, as we quantify by scoring phenotypic penetrance, the fraction of mutants expressing the trait. A very interesting feature is that the enhancements are highly specific for fras1 mutant phenotypes – occurring in the apparent absence of more wide-spread changes. Except for the features due to the fras1 mutation, the transgene-bearing larvae appear generally healthy and to be developing normally. The transgene behaves as a genetic partial dominant: A single copy is sufficient for the enhancements, yet, for some traits, two copies may exert a stronger effect. We made new strains bearing independent insertions of the fli1a-F-hsp70l:Gal4VP16 transgene in new locations in the genome, and observed increased severities of the same phenotypes as observed for the original insertion. This finding suggests that sequences within the transgene, e.g. Gal4VP16, are responsible for the enhancements, rather than effect on neighboring host sequences (such as an insertional mutation). The specificity, and biological action underlying the traits, are subjects of considerable interest for further investigation, as we discuss. Our findings show that work with transgenes needs to be undertaken with caution and attention to detail.

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.

Intestinal tract development is a coordinated process involving signaling among the progenitors and developing cells from all three germ layers. Development of endoderm-derived intestinal epithelium has been shown to depend on epigenetic modifications, but whether that is also the case for intestinal tract cell types from other germ layers remains unclear. We found that functional loss of a DNA methylation machinery component, ubiquitin-like protein containing PHD and RING finger domains 1 (uhrf1), leads to reduced numbers of ectoderm-derived enteric neurons and severe disruption of mesoderm-derived intestinal smooth muscle. Genetic chimeras revealed that Uhrf1 functions both cell-autonomously in enteric neuron precursors and cell-non-autonomously in surrounding intestinal cells, consistent with what is known about signaling interactions between these cell types that promote one another’s development. Uhrf1 recruits the DNA methyltransferase Dnmt1 to unmethylated DNA during replication. Dnmt1 is also expressed in enteric neurons and smooth muscle progenitors. dnmt1 mutants have fewer enteric neurons and disrupted intestinal smooth muscle compared to wildtypes. Because dnmt1;uhrf1 double mutants have a similar phenotype to dnmt1 and uhrf1 single mutants, Dnmt1 and Uhrf1 must function together during enteric neuron and intestinal muscle development. This work shows that genes controlling epigenetic modifications are important to coordinate intestinal tract development, provides the first demonstration that these genes influence development of the ENS, and advances uhrf1 and dnmt1 as potential new Hirschsprung disease candidates.

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.

Fetal mammalian testes secrete Amh (Anti-Müllerian hormone), which inhibits female reproductive tract (Müllerian duct) development. Amh also derives from mature mammalian ovarian follicles, which marks oocyte reserve and characterizes PCOS (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 non-vitellogenic follicles in exceedingly large but sterile ovaries, showing that Amh helps control ovarian follicle maturation and proliferation. RNA-seq data partitioned 21-day post-fertilization (dpf) juveniles into two groups that each contained mutant and wild type fish. Group21-1 up-regulated ovary genes compared to Group21-2, which were likely developing as males. By 35dpf, transcriptomes distinguished males from females and, within each sex, mutants from wild types. In adult mutants, ovaries greatly under-expressed granulosa and theca genes and testes under-expressed 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.Article Summary Anti-Müllerian hormone (Amh) inhibits female reproductive duct development, signals oocyte reserve, and marks polycystic ovarian syndrome. Zebrafish lacks Müllerian ducts and the typical Amh receptor, questioning evolving roles of Amh. We made knockout mutations in zebrafish amh . Most mutants were female and the few males often had oocytes in their testes, showing that Amh promotes male development. Mutant reproductive ducts functioned, but testes were enormous and ovaries accumulated immature oocytes, showing that Amh regulates germ cell proliferation and maturation. Transcriptomics revealed that Amh controls development of steroid-producing gonad cells. Amh in zebrafish preserved ancestral roles despite losing Müllerian ducts and the Amh receptor.

The Indian garden lizard, Calotes versicolor, lacks cytologically recognizable sex chromosomes, and its mechanism of sex determination is unclear. We evaluated genotype-to-sex-phenotype association using RAD-seq in wild-caught males and females, 30 of each sex. Of 210,736 unique, 96-nt long RAD-tags, 48% contained polymorphisms, 23% of which were present in at least 40 of 60 individuals. Twenty one RAD-tags neared, but none achieved, the inclusion criteria for sex enrichment, as expected if C. versicolor lacks highly differentiated sex chromosomes. Three RAD-tags with alleles most strongly associated with sex tended to be heterozygous in females and to lack male-specific alleles, suggesting a ZW female/ZZ male system. Putative female alleles, however, were present in some males and lacking from some females, suggesting either recombination between these markers and the sex locus or sex reversal due to environmental or genetic factors. Paired-end, 250-nt reads from 1 male provided a fragmented draft genome assembly. Four sex-associated RAD-tags were identical to portions of 4 unique C. versicolor genomic contigs rather than linked to a single putative sex-linked region. The lack of strongly sex-linked loci coupled with weak evidence for temperature-associated sex determination intensifies the need for further investigation of the puzzling sex determination mechanism in C. versicolor.

Anti-Mullerian hormone (Amh) inhibits female reproductive duct development, signals oocyte reserve, and marks polycystic ovarian syndrome. Zebrafish lacks Mullerian ducts and the typical Amh receptor, questioning evolving roles of Amh. Yan et al. made knockout mutations in zebrafish... 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.

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.

People living a subsistence lifestyle in the Arctic are highly exposed to persistent organic pollutants, including polychlorinated biphenyls (PCBs). Formerly Used Defense (FUD) sites are point sources of PCB pollution; the Arctic contains thousands of FUD sites, many co-located with indigenous villages. We investigated PCB profiles and biological effects in freshwater fish (Alaska blackfish [Dallia pectoralis] and ninespine stickleback [Pungitius pungitius]) living upstream and downstream of the Northeast Cape FUD site on St. Lawrence Island in the Bering Sea. Despite extensive site remediation, fish remained contaminated with PCBs. Vitellogenin concentrations in males indicated exposure to estrogenic contaminants, and some fish were hypothyroid. Downstream fish showed altered DNA methylation in gonads and altered gene expression related to DNA replication, response to DNA damage, and cell signaling. This study demonstrates that, even after site remediation, contaminants from Cold War FUD sites in remote regions of the Arctic remain a potential health threat to local residents – in this case, Yupik people who had no influence over site selection and use by the United States military.

Intestinal epithelium development depends on epigenetic modifications, but whether that is also the case for other intestinal tract cell types remains unclear. We found that functional loss of a DNA methylation machinery component, ubiquitin-like protein containing PHD and RING finger domains 1 (uhrf1), leads to reduced enteric neuron number, changes in neuronal morphology, and severe intestinal smooth muscle disruption. Genetic chimeras revealed that Uhrf1 functions both cell-autonomously in enteric neuron progenitors and cell-non-autonomously in surrounding intestinal cells. Uhrf1 recruits the DNA methyltransferase Dnmt1 to unmethylated DNA during replication. Dnmt1 is also expressed in enteric neuron and smooth muscle progenitors. dnmt1 mutants show a strong reduction in enteric neuron number and disrupted intestinal smooth muscle. Because dnmt1 ; uhrf1 double mutants have a similar phenotype to dnmt1 and uhrf1 single mutants, Dnmt1 and Uhrf1 must function together during enteric neuron and intestinal muscle development. This work shows that genes controlling epigenetic modifications are important in coordinating intestinal tract development, provides the first demonstration that these genes are important in ENS development, and advances uhrf1 and dnmt1 as potential new Hirschsprung disease candidates. Summary This work provides evidence that DNA methylation factors are important in all cell types that contribute to development of a functional intestine.

MicroRNAs (miRNAs) are small RNA molecules (20-27 nt long) that are involved in eukaryotic gene regulation. They regulate targeted genes by RNA complementarity, generally between the miRNA seed region and the 3’ UTR of messenger RNAs. They have become relevant to the study of developmental stages and diseases. With the advance of sequencing technology, the cost of detecting the miRNA transcriptome has decreased substantially, and has led to the discovery of isomiRs (with slight sequence variations relative to the annotated miRNAs). As a consequence, the amount of miRNA data has increased exponentially in the last 5 years, together with the number of pipelines to analyze it. However, there is still a lack of consensus standards in storing and sharing such data, generating diverse output files; hindering comparison between tools, data sharing, and development of downstream analyses independent of the tools used for miRNA detection and quantification. Here, we present a community based, open source project to work toward the standardization of miRNA pipelines and encourage the development of downstream tools for visualization, differential expression, sample clustering, and model prediction analyses. This project is an international collaboration, with experts from different countries that have been developing miRNA pipelines and resources. We have described a standard format for the output of miRNA detection and quantification tools using small RNA-seq data. The format is based on the GFF3 standard in order to support reference coordinates and parent/child relationships between the features. The format contains information foreach sequence and its annotation to the miRNA precursor, the definition of reference or isomiR sequence, its quality, the isomiR type, and abundance in the data set. Moreover, we support a command line python tool to manage the miRNA GFF3 format ( miRTop ). Currently, miRTop can convert the output of commonly used small RNA-Seq pipelines, such as seqbuster, isomiR-SEA, sRNAbench, and Prost!, as well as BAM files to the miRNA GFF3 format. Importantly, the miRge pipeline has adapted the GFF3 format natively. miRTop can convert miRNA GFF3 files to a count matrix that can be easily imported to any downstream tool (i.e. for differential expression analysis). Furthermore, the tool can export the miRNA GFF file to the isomiRs package, improving the usability of this tool and rendering it independent of the upstream methods used for quantifying miRNA sequences. Together, we believe that the miRTop project will improve accessibility and uniformity of the results oft miRNA pipelines, the possibility to easily run benchmarking analyses, and promote the development of miRNA tools downstream of the detection and quantification steps.

Whole genome duplications (WGD) are important evolutionary events. Our understanding of underlying mechanisms, including the evolution of duplicated genes after WGD, however remains incomplete. Teleost fish experienced a common WGD (teleost-specific genome duplication, or TGD) followed by a dramatic adaptive radiation leading to more than half of all vertebrate species. The analysis of gene expression patterns following TGD at the genome level has been limited by the lack of suitable genomic resources. The recent concomitant release of the genome sequence of spotted gar (a representative of holosteans, the closest lineage of teleosts that lacks the TGD) and the tissue-specific gene expression repertoires of over 20 holostean and teleostean fish species, including spotted gar, zebrafish and medaka (the PhyloFish project), offered a unique opportunity to study the evolution of gene expression following TGD in teleosts. We show that most TGD duplicates gained their current status (loss of one duplicate gene or retention of both duplicates) relatively rapidly after TGD (i.e. prior to the divergence of medaka and zebrafish lineages). The loss of one duplicate is the most common fate after TGD with a probability of approximately 80%. In addition, the fate of duplicate genes after TGD, including subfunctionalization, neofunctionalization, or retention of two ‘similar’ copies occurred not only before, but also after the radiation of species tested, in consistency with a role of the TGD in speciation and/or evolution of gene function. Finally, we report novel cases of TGD ohnolog subfunctionalization and neofunctionalization that further illustrate the importance of these processes.

Gar is an actinopterygian that has bone, dentin, enameloid, and ganoin (enamel) in teeth and/or scales. Mineralization of these tissues involves genes encoding various secretory calcium-binding phosphoproteins (SCPPs) in osteichthyans, but no SCPP genes have been identified in chondrichthyans to date. In the gar genome, we identified 38 SCPP genes, seven of which encode acidic-residue-rich proteins and 31 encode Pro/Gln (P/Q) rich proteins. These gar SCPP genes constitute the largest known repertoire, including many newly identified P/Q-rich genes expressed in teeth and/or scales. Among gar SCPP genes, six acidic and three P/Q-rich genes were identified as orthologs of sarcopterygian genes. The sarcopterygian orthologs of most of these acidic genes are involved in bone and/or dentin formation, and sarcopterygian orthologs of all three P/Q-rich genes participate in enamel formation. The finding of these genes in gar suggests that an elaborate SCPP gene-based genetic system for tissue mineralization was already present in stem osteichthyans. While SCPP genes have been thought to originate from ancient SPARCL1, SPARCL1L1 appears to be more closely related to these genes, because it established a structure similar to acidic SCPP genes probably in stem gnathostomes, perhaps at about the same time with the origin of tissue mineralization. Assuming enamel evolved in stem osteichthyans, all P/Q-rich SCPP genes likely arose within the osteichthyan lineage. Furthermore, the absence of acidic SCPP genes in chondrichthyans might be explained by the secondary loss of earliest acidic genes. It appears that many SCPP genes expanded rapidly in stem osteichthyans and in basal actinopterygians.

Nat. Genet. 48, 427–437 (2016); published online 7 March 2016; corrected after print 25 April 2016 As we intended, other researchers have been able to use the draft spotted gar genome sequence available from the Broad Institute website since December 2011, the assembly LepOcu1 publicly available from NCBI since 13 January 2012 under accession code GCA000242695.

Defects in mid-facial development, including cleft lip/palate, account for a large number of human birth defects annually. In many cases, aberrant gene expression results in either a reduction in the number of neural crest cells (NCCs) that reach the frontonasal region and form much of the facial skeleton or subsequent failure of NCC patterning and differentiation into bone and cartilage. While loss of gene expression is often associated with developmental defects, aberrant upregulation of expression can also be detrimental. microRNAs (miRNAs) are a class of non-coding RNAs that normally repress gene expression by binding to recognition sequences located in the 3' UTR of target mRNAs. miRNAs play important roles in many developmental systems, including midfacial development. Here, we take advantage of high throughput RNA sequencing (RNA-seq) from different tissues of the developing mouse midface to interrogate the miRs that are expressed in the midface and select a subset for further expression analysis. Among those examined, we focused on four that showed the highest expression level in in situ hybridization analysis. Mir23b and Mir24.1 are specifically expressed in the developing mouse frontonasal region, in addition to areas in the perichondrium, tongue musculature and cranial ganglia. Mir23b is also expressed in the palatal shelves and in anterior epithelium of the palate. In contrast, Mir133b and Mir128.2 are mainly expressed in head and trunk musculature. Expression analysis of mir23b and mir133b in zebrafish suggests that mir23b is expressed in the pharyngeal arch, otic vesicle, and trunk muscle while mir133b is similarly expressed in head and trunk muscle. Functional analysis by overexpression of mir23b in zebrafish leads to broadening of the ethmoid plate and aberrant cartilage structures in the viscerocranium, while overexpression of mir133b causes a reduction in ethmoid plate size and a significant midfacial cleft. These data illustrate that miRs are expressed in the developing midface and that Mir23b and Mir133b may have roles in this developmental process.

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.