An atlas of seven zebrafish hox cluster mutants provides insights into sub/neofunctionalization of vertebrate Hox clusters

Live Visualization of Calcified Bones in Zebrafish and Medaka Larvae and Juveniles Using Calcein and Alizarin Red S

Zebrafish and medaka are valuable model vertebrates for genetic studies. The advent of CRISPR-Cas9 technology has greatly enhanced our capability to produce specific gene mutants in zebrafish and medaka. Analyzing the phenotypes of these mutants is essential for elucidating gene function, though such analyses often yield unexpected results. Consequently, providing researchers with accessible and cost-effective phenotype analysis methods is crucial. A prevalent technique for investigating calcified bone development in these species involves using transgenic fish that express fluorescent proteins labeling calcified bones; however, acquiring these fish and isolating appropriate crosses can be time-consuming. We present a comprehensive protocol for visualizing ossified bones in zebrafish and medaka larvae and juveniles using calcein and alizarin red S staining, which is both economical and efficient. This method, applicable to live specimens during the ossification of bones, avoids apparent alterations in skeletal morphology and allows for the use of different fluorescent dyes in conjunction with transgenic labeling, thus enhancing the analysis of developmental processes in calcifying bones, such as vertebrae and fin rays. Key features • The calcified bones of alive zebrafish and medaka larvae and juveniles can be visualized repeatedly using simple and inexpensive calcein and alizarin red S. • No need to use transgenic fish to visualize ossified bones, allowing for rapid analysis of bone phenotypes in mutants. • Double staining is possible in transgenic fish with reporter genes such as GFP and DsRed using alizarin red S or calcein, which exhibit different fluorescence. • Ossification processes of bones such as vertebrae, ribs, and fin rays can be analyzed in mutants.

Zebrafish: unraveling genetic complexity through duplicated genes

The zebrafish is an invaluable model organism for genetic, developmental, and disease research. Although its high conservation with humans is often cited as justification for its use, the zebrafish harbors oft-ignored genetic characteristics that may provide unique insights into gene structure and function. Zebrafish, along with other teleost fish, underwent an additional round of whole genome duplication after their split from tetrapods-resulting in an abundance of duplicated genes when compared to other vertebrates. These duplicated genes have evolved in distinct ways over the ensuing 350 million years. Thus, each gene within a duplicated gene pair has nuanced differences that create a unique identity. By investigating both members of the gene pair together, we can elucidate the mechanisms that underly protein structure and function and drive the complex interplay within biological systems, such as signal transduction cascades, genetic regulatory networks, and evolution of tissue and organ function. It is crucial to leverage such studies to explore these molecular dynamics, which could have far-reaching implications for both basic science and therapeutic development. Here, we will review the role of gene duplications and the existing models for gene divergence and retention following these events. We will also highlight examples within each of these models where studies comparing duplicated genes in the zebrafish have yielded key insights into protein structure, function, and regulation.

Physical constraints on the positions and dimensions of the zebrafish swim bladder by surrounding bones

Precise regulation of organ size and position is crucial for optimal organ function. Since the swim bladder is primarily responsible for buoyancy in teleosts, early development and subsequent inflation of the swim bladder should be appropriately controlled with the body growth. However, the underlying mechanism remains unclear. In this study, we show that the size and position of the swim bladder are physically constrained by the surrounding bones in zebrafish. Non-invasive micro-CT scanning revealed that the anterior edge of the swim bladder is largely attached to the os suspensorium, which is an ossicle extending medioventrally from the 4th centrum. Additionally, we observed that hoxc6a mutants, which lack the os suspensorium, exhibited an anterior projection of the swim bladder beyond the 4th vertebra. During the swim bladder development, we found that the counterclockwise rotation of the os suspensorium correlates with posterior regression of the swim bladder, suggesting that the os suspensorium pushes the swim bladder posteriorly into its proper position. Furthermore, our results revealed a close association between the posterior region of the swim bladder and the pleural ribs. In hoxaa cluster mutants with additional ribs, the swim bladder expanded posteriorly, accompanied by an enlarged body cavity. Taken together, our results demonstrate the importance of the surrounding bones in the robust regulation of swim bladder size and position in zebrafish.

The functional roles of zebrafish HoxA- and HoxD-related clusters in the pectoral fin development

The paralogs 9-13 Hox genes in mouse HoxA and HoxD clusters are critical for limb development. When both HoxA and HoxD clusters are deleted in mice, significant limb truncation is observed compared to the phenotypes of single and compound mutants of Hox9-13 genes in these clusters. In zebrafish, mutations in hox13 genes in HoxA- and HoxD-related clusters result in abnormal morphology of pectoral fins, homologous to forelimbs. However, the effect of the simultaneous deletions of entire HoxA- and HoxD-related clusters on pectoral fin development remains unknown. Here, we generated mutants with several combinations of hoxaa, hoxab, and hoxda cluster deletions and analyzed the pectoral fin development. In hoxaa-/-;hoxab-/-;hoxda-/- larvae, the endoskeletal disc and the fin-fold are significantly shortened in developing pectoral fins. In addition, we show that this anomaly is due to defects in the pectoral fin growth after the fin bud formation. Furthermore, in the surviving adult mutants, micro-CT scanning reveals defects in the posterior portion of the pectoral fin which is thought to represent latent regions of the limb. Our results further support that the functional role of HoxA and HoxD clusters is conserved in the paired appendage formation in bony fishes.

Genome sequencing and assembly of Indian major carp, Cirrhinus mrigala (Hamilton, 1822)

The freshwater aquaculture species, Cirrhinus mrigala, commonly known as mrigal, holds significance in the carp-dominated aquaculture system, globally. Despite constituting 1.08% of the total freshwater aquaculture production, mrigal is the third most important Indian major carp. However, its genome and associated information is not available. This study aims to address this gap by generation high quality genome assembly using PacBio long reads, Illumina short reads and Hi-C scaffolding. The characterization of assembled highly contiguous genome, 1.057 Gb in size, revealed 39,091 genes with functional annotations. The orthology analysis based on direct orthologs and single copy ortholgs places C. mrigala in a distinct position within the Otophysi clade. Additionally, the study delves into Hox gene clusters, identifying 38 Hox genes distributed in seven clusters. The present genomic information offers potential applications for sustainable aquaculture management, including selection programs for economic traits.

The Hox code responsible for the patterning of the anterior vertebrae in zebrafish

The vertebral column is a characteristic structure of vertebrates. Genetic studies in mice have shown that Hox-mediated patterning plays a key role in specifying discrete anatomical regions of the vertebral column. Expression pattern analyses in several vertebrate embryos have provided correlative evidence that the anterior boundaries of Hox expression coincide with distinct anatomical vertebrae. However, because functional analyses have been limited to mice, it remains unclear which Hox genes actually function in vertebral patterning in other vertebrates. In this study, various zebrafish Hox mutants were generated for loss-of-function phenotypic analysis to functionally decipher the Hox code responsible for the zebrafish anterior vertebrae between the occipital and thoracic vertebrae. We found that Hox genes in HoxB- and HoxC-related clusters participate in regulating the morphology of the zebrafish anterior vertebrae. In addition, medaka hoxc6a was found to be responsible for anterior vertebral identity, as in zebrafish. Based on phenotypic similarities with Hoxc6 knockout mice, our results suggest that the Hox patterning system, including at least Hoxc6, may have been functionally established in the vertebral patterning of the common ancestor of ray-finned and lobe-finned fishes.

Teleost Hox code defines regional identities competent for the formation of dorsal and anal fins

The dorsal and anal fins can vary widely in position and length along the anterior-posterior axis in teleost fishes. However, the molecular mechanisms underlying the diversification of these fins remain unknown. Here, we used genetic approaches in zebrafish and medaka, in which the relative positions of the dorsal and anal fins are opposite, to demonstrate the crucial role of hox genes in the patterning of the teleost posterior body, including the dorsal and anal fins. By the CRISPR-Cas9-induced frameshift mutations and positional cloning of spontaneous dorsalfinless medaka, we show that various hox mutants exhibit the absence of dorsal or anal fins, or a stepwise posterior extension of these fins, with vertebral abnormalities. Our results indicate that multiple hox genes, primarily from hoxc-related clusters, encompass the regions responsible for the dorsal and anal fin formation along the anterior-posterior axis. These results further suggest that shifts in the anterior boundaries of hox expression which vary among fish species, lead to diversification in the position and size of the dorsal and anal fins, similar to how modulations in Hox expression can alter the number of anatomically distinct vertebrae in tetrapods. Furthermore, we show that hox genes responsible for dorsal fin formation are different between zebrafish and medaka. Our results suggest that a novel mechanism has occurred during teleost evolution, in which the gene network responsible for fin formation might have switched to the regulation downstream of other hox genes, leading to the remarkable diversity in the dorsal fin position.

In Search of a Target Gene for a Desirable Phenotype in Aquaculture: Genome Editing of Cyprinidae and Salmonidae Species

Aquaculture supplies the world food market with a significant amount of valuable protein. Highly productive aquaculture fishes can be derived by utilizing genome-editing methods, and the main problem is to choose a target gene to obtain the desirable phenotype. This paper presents a review of the studies of genome editing for genes controlling body development, growth, pigmentation and sex determination in five key aquaculture Salmonidae and Cyprinidae species, such as rainbow trout (Onchorhynchus mykiss), Atlantic salmon (Salmo salar), common carp (Cyprinus carpio), goldfish (Carassius auratus), Gibel carp (Carassius gibelio) and the model fish zebrafish (Danio rerio). Among the genes studied, the most applicable for aquaculture are mstnba, pomc, and acvr2, the knockout of which leads to enhanced muscle growth; runx2b, mutants of which do not form bones in myoseptae; lepr, whose lack of function makes fish fast-growing; fads2, Δ6abc/5Mt, and Δ6bcMt, affecting the composition of fatty acids in fish meat; dnd mettl3, and wnt4a, mutants of which are sterile; and disease-susceptibility genes prmt7, gab3, gcJAM-A, and cxcr3.2. Schemes for obtaining common carp populations consisting of only large females are promising for use in aquaculture. The immobilized and uncolored zebrafish line is of interest for laboratory use.

Hox genes and patterning the vertebrate body

The diversity of vertebrate body plans is dizzying, yet stunning for the many things they have in common. Vertebrates have inhabited virtually every part of the earth from its coldest to warmest climates. They locomote by swimming, flying, walking, slithering, or climbing, or combinations of these behaviors. And they exist in many different sizes, from the smallest of frogs, fish and lizards to giraffes, elephants, and blue whales. Despite these differences, vertebrates follow a remarkably similar blueprint for the establishment of their body plan. Within the relatively small amount of time required to complete gastrulation, the process through which the three germ layers, ectoderm, mesoderm, and endoderm are created, the embryo also generates its body axis and is simultaneously patterned. For the length of this axis, the genes that distinguish the neck from the rib cage or the trunk from the sacrum are the Hox genes. In vertebrates, there was evolutionary pressure to maintain this set of genes in the organism. Over the past decades, much has been learned regarding the regulatory mechanisms that ensure the appropriate expression of these genes along the main body axes. Genetic functions continue to be explored though much has been learned. Much less has been discerned on the identity of co-factors used by Hox proteins for the specificity of transcriptional regulation or what downstream targets and pathways are critical for patterning events, though there are notable exceptions. Current work in the field is demonstrating that Hox genes continue to function in many organs long after directing early patterning events. It is hopeful continued research will shed light on remaining questions regarding mechanisms used by this important and conserved set of transcriptional regulators.

Human split hand/foot variants are not as functional as wildtype human PRDM1 in the rescue of craniofacial defects

BACKGROUND

Split hand/foot malformation (SHFM) is a congenital limb disorder presenting with limb anomalies, such as missing, hypoplastic, or fused digits, and often craniofacial defects, including a cleft lip/palate, microdontia, micrognathia, or maxillary hypoplasia. We previously identified three novel variants in the transcription factor, PRDM1, that are associated with SHFM phenotypes. One individual also presented with a high arch palate. Studies in vertebrates indicate that PRDM1 is important for development of the skull; however, prior to our study, human variants in PRDM1 had not been associated with craniofacial anomalies.

METHODS

Using transient mRNA overexpression assays in prdm1a-/- mutant zebrafish, we tested whether the PRDM1 SHFM variants were functional and could lead to a rescue of the craniofacial defects observed in prdm1a-/- mutants. We also mined previously published CUT&RUN and RNA-seq datasets that sorted EGFP-positive cells from a Tg(Mmu:Prx1-EGFP) transgenic line that labels the pectoral fin, pharyngeal arches, and dorsal part of the head to examine Prdm1a binding and the effect of Prdm1a loss on craniofacial genes.

RESULTS

The prdm1a-/- mutants exhibit craniofacial defects including a hypoplastic neurocranium, a loss of posterior ceratobranchial arches, a shorter palatoquadrate, and an inverted ceratohyal. Injection of wildtype (WT) hPRDM1 in prdm1a-/- mutants partially rescues the palatoquadrate phenotype. However, injection of each of the three SHFM variants fails to rescue this skeletal defect. Loss of prdm1a leads to a decreased expression of important craniofacial genes by RNA-seq, including emilin3a, confirmed by hybridization chain reaction expression. Other genes including dlx5a/dlx6a, hand2, sox9b, col2a1a, and hoxb genes are also reduced. Validation by real-time quantitative PCR in the anterior half of zebrafish embryos failed to confirm the expression changes suggesting that the differences are enriched in prx1 expressing cells.

CONCLUSION

These data suggest that the three SHFM variants are likely not functional and may be associated with the craniofacial defects observed in the humans. Finally, they demonstrate how Prdm1a can directly bind and regulate genes involved in craniofacial development.

Morphological description of spontaneous pelvic fin loss in a neotropical cichlid fish

Pelvic fins are a characteristic structure of the vertebrate Bauplan. Yet, pelvic fin loss has occurred repeatedly across a wide diversity of other lineages of tetrapods and at least 48 times in teleost fishes. This pelvic finless condition is often associated with other morphological features such as body elongation, loss of additional structures, and bilateral asymmetry. However, despite the remarkable diversity in the several thousand cichlid fish species, none of them are characterized by the complete absence of pelvic fins. Here, we examined the musculoskeletal structure and associated bilateral asymmetry in Midas cichlids (Amphilophus cf. citrinellus) that lost their pelvic fins spontaneously in the laboratory. Due to this apparent mutational loss of the pelvic girdle and fins, the external and internal anatomy are described in a series of "normal" Midas individuals and their pelvic finless sibling tankmates. First, other traits associated with teleost pelvic fin loss, the genetic basis of pelvic fin loss, and the potential for pleiotropic effects of these genes on other traits in teleosts were all reviewed. Using these traits as a guide, we investigated whether other morphological differences were associated with the pelvic girdle/fin loss. The mean values of the masses of muscle of the pectoral fin, fin ray numbers in the unpaired fins, and oral jaw tooth numbers did not differ between the two pelvic fin morphotypes. However, significant differences in meristic values of the paired traits assessed were observed for the same side of the body between morphotypes. Notably, bilateral asymmetry was found exclusively for the posterior lateral line scales. Finally, we found limited evidence of pleiotropic effects, such as lateral line scale numbers and fluctuating asymmetry between the Midas pelvic fin morphotypes. The fast and relatively isolated changes in the Midas cichlids suggest minor but interesting pleiotropic effects could accompany loss of cichlid pelvic fins.

Testing biological actions of medicinal plants from northern Vietnam on zebrafish embryos and larvae: Developmental, behavioral, and putative therapeutical effects

Evaluating the risks and benefits of using traditional medicinal plants is of utmost importance for a huge fraction of the human population, in particular in Northern Vietnam. Zebrafish are increasingly used as a simple vertebrate model for testing toxic and physiological effects of compounds, especially on development. Here, we tested 12 ethanolic extracts from popular medicinal plants collected in northern Vietnam for their effects on zebrafish survival and development during the first 4 days after fertilization. We characterized more in detail their effects on epiboly, hatching, growth, necrosis, body curvature, angiogenesis, skeletal development and mostly increased movement behavior. Finally, we confirm the effect on epiboly caused by the Mahonia bealei extract by staining the actin filaments and performing whole genome gene expression analysis. Further, we show that this extract also inhibits cell migration of mouse embryo fibroblasts. Finally, we analyzed the chemical composition of the Mahonia bealei extract and test the effects of its major components. In conclusion, we show that traditional medicinal plant extracts are able to affect zebrafish early life stage development to various degrees. In addition, we show that an extract causing delay in epiboly also inhibits mammalian cell migration, suggesting that this effect may serve as a preliminary test for identifying extracts that inhibit cancer metastasis.

Zebrafish as a Model to Study Retinoic Acid Signaling in Development and Disease

Retinoic acid (RA) is a metabolite of vitamin A (retinol) that plays various roles in development to influence differentiation, patterning, and organogenesis. RA also serves as a crucial homeostatic regulator in adult tissues. The role of RA and its associated pathways are well conserved from zebrafish to humans in both development and disease. This makes the zebrafish a natural model for further interrogation into the functions of RA and RA-associated maladies for the sake of basic research, as well as human health. In this review, we explore both foundational and recent studies using zebrafish as a translational model for investigating RA from the molecular to the organismal scale.

Hox genes in development and beyond

Hox genes encode evolutionarily conserved transcription factors that are essential for the proper development of bilaterian organisms. Hox genes are unique because they are spatially and temporally regulated during development in a manner that is dictated by their tightly linked genomic organization. Although their genetic function during embryonic development has been interrogated, less is known about how these transcription factors regulate downstream genes to direct morphogenetic events. Moreover, the continued expression and function of Hox genes at postnatal and adult stages highlights crucial roles for these genes throughout the life of an organism. Here, we provide an overview of Hox genes, highlighting their evolutionary history, their unique genomic organization and how this impacts the regulation of their expression, what is known about their protein structure, and their deployment in development and beyond.

Genomic Characterization of hox Genes in Senegalese Sole (Solea senegalensis, Kaup 1858): Clues to Evolutionary Path in Pleuronectiformes

The Senegalese sole (Solea senegalensis, Kaup 1858), a marine flatfish, belongs to the Pleuronectiformes order. It is a commercially important species for fisheries and aquaculture. However, in aquaculture, several production bottlenecks have still to be resolved, including skeletal deformities and high mortality during the larval and juvenile phase. The study aims to characterize the hox gene clusters in S. senegalensis to understand better the developmental and metamorphosis process in this species. Using a BAC library, the clones that contain hox genes were isolated, sequenced by NGS and used as BAC-FISH probes. Subsequently the hox clusters were studied by sequence analysis, comparative genomics, and cytogenetic and phylogenetic analysis. Cytogenetic analysis demonstrated the localization of four BAC clones on chromosome pairs 4, 12, 13, and 16 of the Senegalese sole cytogenomic map. Comparative and phylogenetic analysis showed a highly conserved organization in each cluster and different phylogenetic clustering in each hox cluster. Analysis of structural and repetitive sequences revealed accumulations of polymorphisms mediated by repetitive elements in the hoxba cluster, mainly retroelements. Therefore, a possible loss of the hoxb7a gene can be established in the Pleuronectiformes lineage. This work allows the organization and regulation of hox clusters to be understood, and is a good base for further studies of expression patterns.

MicroRNAs as Indicators into the Causes and Consequences of Whole-Genome Duplication Events

Whole-genome duplications (WGDs) have long been considered the causal mechanism underlying dramatic increases to morphological complexity due to the neo-functionalization of paralogs generated during these events. Nonetheless, an alternative hypothesis suggests that behind the retention of most paralogs is not neo-functionalization, but instead the degree of the inter-connectivity of the intended gene product, as well as the mode of the WGD itself. Here, we explore both the causes and consequences of WGD by examining the distribution, expression, and molecular evolution of microRNAs (miRNAs) in both gnathostome vertebrates as well as chelicerate arthropods. We find that although the number of miRNA paralogs tracks the number of WGDs experienced within the lineage, few of these paralogs experienced changes to the seed sequence, and thus are functionally equivalent relative to their mRNA targets. Nonetheless, in gnathostomes, although the retention of paralogs following the 1R autotetraploidization event is similar across the two subgenomes, the paralogs generated by the gnathostome 2R allotetraploidization event are retained in higher numbers on one subgenome relative to the second, with the miRNAs found on the preferred subgenome showing both higher expression of mature miRNA transcripts and slower molecular evolution of the precursor miRNA sequences. Importantly, WGDs do not result in the creation of miRNA novelty, nor do WGDs correlate to increases in complexity. Instead, it is the number of miRNA seed sequences in the genome itself that not only better correlate to instances in complexification, but also mechanistically explain why complexity increases when new miRNA families are established.

Homeotic Genes: Clustering, Modularity, and Diversity

Hox genes code for transcription factors and are evolutionarily conserved. They regulate a plethora of downstream targets to define the anterior-posterior (AP) body axis of a developing bilaterian embryo. Early work suggested a possible role of clustering and ordering of Hox to regulate their expression in a spatially restricted manner along the AP axis. However, the recent availability of many genome assemblies for different organisms uncovered several examples that defy this constraint. With recent advancements in genomics, the current review discusses the arrangement of Hox in various organisms. Further, we revisit their discovery and regulation in Drosophila melanogaster. We also review their regulation in different arthropods and vertebrates, with a significant focus on Hox expression in the crustacean Parahyale hawaiensis. It is noteworthy that subtle changes in the levels of Hox gene expression can contribute to the development of novel features in an organism. We, therefore, delve into the distinct regulation of these genes during primary axis formation, segment identity, and extra-embryonic roles such as in the formation of hair follicles or misregulation leading to cancer. Toward the end of each section, we emphasize the possibilities of several experiments involving various organisms, owing to the advancements in the field of genomics and CRISPR-based genome engineering. Overall, we present a holistic view of the functioning of Hox in the animal world.