Transgene-mediated skeletal phenotypic variation in zebrafish.

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, harbours a huge amount of foreign DNA. Such DNA, in the form of transposable elements, can affect genome function in a major way, and transgene biology needs to be included in our understanding of the genome. Here we examine an unexpected phenotypic effect of the chromosomally integrated transgene fli1a-F-hsp70l:Gal4VP16 that serves as a model for transgene function generally. We examine larval fras1 mutant zebrafish (Danio rerio). Gal4VP16 is a potent transcriptional activator that is 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 widespread 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, for example Gal4VP16, are responsible for the enhancements, rather than the effect on neighbouring 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.

Epigenetic factors Dnmt1 and Uhrf1 coordinate intestinal 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.

SCPP Genes and Their Relatives in Gar: Rapid Expansion of Mineralization Genes in Osteichthyans.

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.

Evolution of gene expression after whole-genome duplication: New insights from the spotted gar genome.

Whole-genome duplications (WGDs) 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-related 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), offers 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 divergence 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.

Transcriptomes of post-mitotic neurons identify the usage of alternative pathways during adult and embryonic neuronal differentiation.

BACKGROUND: Understanding the mechanisms by which neurons are generated and specified, and how they integrate into functional circuits is key to being able to treat disorders of the nervous system and acute brain trauma. Much of what we know about neuronal differentiation has been studied in developing embryos, but differentiation steps may be very different during adult neurogenesis. For this reason, we compared the transcriptomes of newly differentiated neurons in zebrafish embryos and adults. RESULTS: Using a 4tU RNA labeling method, we isolated and sequenced mRNA specifically from cells of one day old embryos and adults expressing the transgene HA-uprt-mcherry under control of the neuronal marker elavl3. By categorizing transcript products into different protein classes, we identified similarities and differences of gene usage between adult and embryonic neuronal differentiation. We found that neurons in the adult brain and in the nervous system of one day old embryos commonly use transcription factors - some of them identical - during the differentiation process. When we directly compared adult differentiating neurons to embryonic differentiating neurons, however, we found that during adult neuronal differentiation, the expression of neuropeptides and neurotransmitter pathway genes is more common, whereas classical developmental signaling through secreted molecules like Hedgehog or Wnt are less enriched, as compared to embryonic stages. CONCLUSIONS: We conclude that both adult and embryonic differentiating neurons show enriched use of transcription factors compared to surrounding cells. However, adult and embryonic developing neurons use alternative pathways to differentiate. Our study provides evidence that adult neuronal differentiation is distinct from the better characterized embryonic neuronal differentiation process. This important insight and the lists of enriched genes we have identified will now help pave the way to a better understanding of the mechanisms of embryonic and adult neuronal differentiation and how to manipulate these processes.

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.

Direct male development in chromosomally ZZ 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.

Direct Male Development in Chromosomally ZZ 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 strain 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.

Evolution and developmental expression of the sodium-iodide symporter (NIS, slc5a5) gene family: Implications for perchlorate toxicology.

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.

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.

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.

A Search for Sex-Linked Loci in the Agamid Lizard, Calotes versicolor.

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.

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.

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.

Endocrine disruption and differential gene expression in sentinel fish on St. Lawrence Island, Alaska: Health implications for indigenous residents.

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.

MicroRNA Profiling during Craniofacial Development: Potential Roles for Mir23b and Mir133b.

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.

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

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

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

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

Wild sex in zebrafish: loss of the natural sex determinant in domesticated strains.

Sex determination can be robustly genetic, strongly environmental, or genetic subject to environmental perturbation. The genetic basis of sex determination is unknown for zebrafish (Danio rerio), a model for development and human health. We used RAD-tag population genomics to identify sex-linked polymorphisms. After verifying this "RAD-sex" method on medaka (Oryzias latipes), we studied two domesticated zebrafish strains (AB and TU), two natural laboratory strains (WIK and EKW), and two recent isolates from nature (NA and CB). All four natural strains had a single sex-linked region at the right tip of chromosome 4, enabling sex genotyping by PCR. Genotypes for the single nucleotide polymorphism (SNP) with the strongest statistical association to sex suggested that wild zebrafish have WZ/ZZ sex chromosomes. In natural strains, "male genotypes" became males and some "female genotypes" also became males, suggesting that the environment or genetic background can cause female-to-male sex reversal. Surprisingly, TU and AB lacked detectable sex-linked loci. Phylogenomics rooted on D. nigrofasciatus verified that all strains are monophyletic. Because AB and TU branched as a monophyletic clade, we could not rule out shared loss of the wild sex locus in a common ancestor despite their independent domestication. Mitochondrial DNA sequences showed that investigated strains represent only one of the three identified zebrafish haplogroups. Results suggest that zebrafish in nature possess a WZ/ZZ sex-determination mechanism with a major determinant lying near the right telomere of chromosome 4 that was modified during domestication. Strains providing the zebrafish reference genome lack key components of the natural sex-determination system but may have evolved variant sex-determining mechanisms during two decades in laboratory culture.