Zebrafish hox clusters and vertebrate genome evolution

The myoarchitecture of the vertebrate cardiac ventricles: evolution and classification

The ventricle of the vertebrate heart is the main segment of the cardiac outflow region. Compared with other cardiac components, it shows remarkable histomorphological variation among different animal groups. This variation is especially apparent in the myocardium, which is generally classified into three main types: trabeculated, compact and mixed. The trabeculated or 'spongy' myocardium is characterized by the existence of trabeculae and deep recesses or intertrabecular spaces, lined by the endocardium. The compact type is composed of condensed myocardial fibers, with almost no trabeculated layer. The mixed type consists of an outer compact layer and an inner trabeculated layer. Among vertebrates, fishes show a great diversity of myocardial types. On this basis, the ventricular myoarchitecture has been categorized into four groups of varying complexity. This classification is made according to (i) the proportion of the two types of myocardium, trabeculated versus compact, and (ii) the vascularization of the heart wall. Here, we review the morphogenetic mechanisms that give rise to the different ventricular myoarchitecture in gnathostomes (i.e. jawed vertebrates) with special emphasis on the diversity of the ventricular myocardium throughout the phylogeny of ancient actinopterygians and teleosts. Finally, we propose that the classification of the ventricular myoarchitecture should be reconsidered, given that the degrees of myocardial compactness on which the current classification system is based do not constitute discrete states, but an anatomical continuum.

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.

A zebrafish gene with sequence similarities to human uromodulin and GP2 displays extensive evolutionary diversification among teleost and confers resistance to bacterial infection

In the process of investigating synaptic changes happening to mutants lacking postsynaptic receptors in the neuromuscular junction, we focused on a hitherto uncharacterized zebrafish gene zgc153932 whose expression was increased in the RNAseq and droplet digital PCR (ddPCR) analysis of a paralyzed mutant sofa potato. The zgc153932 gene which we named omcin5 (omc5) showed amino acid sequence similarity to human uromodulin and GP2, which are expressed in epithelial cells of the kidney and the gut respectively and bind to bacteria pili. omc5 had 14 paralogues in a ∼400 KB region on the chromosome 12 of the zebrafish genome. These genes, named omcin1 through 15, constitute a gene cluster which presumably arose from recent gene duplications in the zebrafish lineage. An antibody raised against the epitope common to 6-9 genes in the omcin family revealed expression in the cloaca of 1 day post fertilization (dpf) embryos which broadened to the urinary and digestive tracts by 5 dpf. Expression of omc5 was increased by exposure of embryos to Escherichia coli (E. coli). Survival of omc5 mutant embryos was shortened in the presence of E. coli, or when they were not maintained in germ-free conditions. Adults omc5 mutants also exhibited susceptibility to infection. Other teleost species which had omcin-like genes in their genomes showed a range of gene duplication, resulting in clusters of 1 to >15 omcin-like genes. We hereby identified a new gene family specific to teleost that include a microbial induced gene which confers resistance to bacterial infection.

Research progress of nanog gene in fish

Nanog is a crucial regulatory factor in maintaining the self-renewal and pluripotency of embryonic stem cells. It is involved in various biological processes, such as early embryonic development, cell reprogramming, cell cycle regulation, the proliferation and migration of primordial germ cells. While research on this gene has primarily focused on mammals, there has been a growing interest in studying nanog in fish. However, there is a notable lack of comprehensive reviews regarding this gene in fish, which is essential for guiding future research. This review aims to provide a thorough summary of the gene's structure, expression patterns, functions and regulatory mechanisms in fish. The findings suggest that nanog probably has both conserved and divergent functions in regulating cell pluripotency, early embryonic development, and germ cell development in teleosts compared to other species, including mammals. These insights lay the foundation for future research and applications of the nanog gene, providing a new perspective for understanding the evolution and conserved charactristics of teleost nanog.

Threonine dehydrogenase regulates neutrophil homeostasis but not H3K4me3 levels in zebrafish

l-threonine dehydrogenase (Tdh) is an enzyme that links threonine metabolism to epigenetic modifications and mitochondria biogenesis. In vitro studies show that it is critical for the regulation of trimethylation of histone H3 lysine 4 (H3K4me3) levels and cell fate determination of mouse embryonic stem cells (mESCs). However, whether Tdh regulates a developmental process in vivo and, if it does, whether it also primarily regulates H3K4me3 levels in this process as it does in mESCs, remains elusive. Here, we revealed that, in zebrafish hematopoiesis, tdh is preferentially expressed in neutrophils. Knockout of tdh causes a decrease in neutrophil number and slightly suppresses their acute injury-induced migration, but, unlike the mESCs, the level of H3K4me3 is not evidently reduced in neutrophils sorted from the kidney marrow of adult tdh-null zebrafish. These phenotypes are dependent on the enzymatic activity of Tdh. Importantly, a soluble supplement of nutrients that are able to fuel the acetyl-CoA pool, such as pyruvate, glucose and branched-chain amino acids, is sufficient to rescue the reduction in neutrophils caused by tdh deletion. In summary, our study presents evidence for the functional requirement of Tdh-mediated threonine metabolism in a developmental process in vivo. It also provides an animal model for investigating the nutritional regulation of myelopoiesis and immune response, as well as a useful tool for high-throughput drug/nutrition screening.

Expansion and Functional Diversification of Long-Wavelength-Sensitive Opsin in Anabantoid Fishes

Gene duplication is one of the most important sources of novel genotypic diversity and the subsequent evolution of phenotypic diversity. Determining the evolutionary history and functional changes of duplicated genes is crucial for a comprehensive understanding of adaptive evolution. The evolutionary history of visual opsin genes is very dynamic, with repeated duplication events followed by sub- or neofunctionalization. While duplication of the green-sensitive opsins rh2 is common in teleost fish, fewer cases of multiple duplication events of the red-sensitive opsin lws are known. In this study, we investigate the visual opsin gene repertoire of the anabantoid fishes, focusing on the five lws opsin genes found in the genus Betta. We determine the evolutionary history of the lws opsin gene by taking advantage of whole-genome sequences of nine anabantoid species, including the newly assembled genome of Betta imbellis. Our results show that at least two independent duplications of lws occurred in the Betta lineage. The analysis of amino acid sequences of the lws paralogs of Betta revealed high levels of diversification in four of the seven transmembrane regions of the lws protein. Amino acid substitutions at two key-tuning sites are predicted to lead to differentiation of absorption maxima (λmax) between the paralogs within Betta. Finally, eye transcriptomics of B. splendens at different developmental stages revealed expression shifts between paralogs for all cone opsin classes. The lws genes are expressed according to their relative position in the lws opsin cluster throughout ontogeny. We conclude that temporal collinearity of lws expression might have facilitated subfunctionalization of lws in Betta and teleost opsins in general.

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.

Disruption of the creb3l1 gene causes defects in caudal fin regeneration and patterning in zebrafish Danio rerio

BACKGROUND

The gene cAMP-Responsive Element Binding protein 3-like-1 (CREB3L1) has been implicated in bone development in mice, with CREB3L1 knock-out mice exhibiting fragile bones, and in humans, with CREB3L1 mutations linked to osteogenesis imperfecta. However, the mechanism through which Creb3l1 regulates bone development is not fully understood.

RESULTS

To probe the role of Creb3l1 in organismal physiology, we used CRISPR-Cas9 genome editing to generate a Danio rerio (zebrafish) model of Creb3l1 deficiency. In contrast to mammalian phenotypes, the Creb3l1 deficient fish do not display abnormalities in osteogenesis, except for a decrease in the bifurcation pattern of caudal fin. Both, skeletal morphology and overall bone density appear normal in the mutant fish. However, the regeneration of caudal fin postamputation is significantly affected, with decreased overall regenerate and mineralized bone area. Moreover, the mutant fish exhibit a severe patterning defect during regeneration, with a significant decrease in bifurcation complexity of the fin rays and distalization of the bifurcation sites. Analysis of genes implicated in bone development showed aberrant patterning of shha and ptch2 in Creb3l1 deficient fish, linking Creb3l1 with Sonic Hedgehog signaling during fin regeneration.

CONCLUSIONS

Our results uncover a novel role for Creb3l1 in regulating tissue growth and patterning during regeneration.

Molecular Dating of the Teleost Whole Genome Duplication (3R) Is Compatible With the Expectations of Delayed Rediploidization

Vertebrate evolution has been punctuated by three whole genome duplication events that have been implicated causally in phenotypic evolution, from the origin of phenotypic novelties to explosive diversification. Arguably, the most dramatic of these is the 3R whole genome duplication event associated with the origin of teleost fishes which comprise more than half of all living vertebrate species. However, tests of a causal relationship between whole genome duplication and teleost diversification have proven difficult due to the challenge of establishing the timing of these phenomena. Here we show, based on molecular clock dating of concatenated gene alignments, that the 3R whole genome duplication event occurred in the early-middle Permian (286.18 to 267.20 million years ago; Ma), 52.02 to 12.84 million years (Myr) before the divergence of crown-teleosts in the latest Permian-earliest Late Triassic (254.36 to 234.16 Ma) and long before the major pulses of teleost diversification in Ostariophysi and Percomorpha (56.37 to 100.17 Myr and at least 139.24 to 183.29 Myr later, respectively). The extent of this temporal gap between putative cause and effect precludes 3R as a deterministic driver of teleost diversification. However, these age constraints remain compatible with the expectations of a prolonged rediploidization process following whole genome duplication which, through the effects of chromosome rearrangement and gene loss, remains a viable mechanism to explain the evolution of teleost novelties and diversification.

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.

DLX genes and proteins in mammalian forebrain development

The vertebrate Dlx gene family encode homeobox transcription factors that are related to the Drosophila Distal-less (Dll) gene and are crucial for development. Over the last ∼35 years detailed information has accrued about the redundant and unique expression and function of the six mammalian Dlx family genes. DLX proteins interact with general transcriptional regulators, and co-bind with other transcription factors to enhancer elements with highly specific activity in the developing forebrain. Integration of the genetic and biochemical data has yielded a foundation for a gene regulatory network governing the differentiation of forebrain GABAergic neurons. In this Primer, we describe the discovery of vertebrate Dlx genes and their crucial roles in embryonic development. We largely focus on the role of Dlx family genes in mammalian forebrain development revealed through studies in mice. Finally, we highlight questions that remain unanswered regarding vertebrate Dlx genes despite over 30 years of research.

Zebrafish ApoB-Containing Lipoprotein Metabolism: A Closer Look

Zebrafish have become a powerful model of mammalian lipoprotein metabolism and lipid cell biology. Most key proteins involved in lipid metabolism, including cholesteryl ester transfer protein, are conserved in zebrafish. Consequently, zebrafish exhibit a human-like lipoprotein profile. Zebrafish with mutations in genes linked to human metabolic diseases often mimic the human phenotype. Zebrafish larvae develop rapidly and externally around the maternally deposited yolk. Recent work revealed that any disturbance of lipoprotein formation leads to the accumulation of cytoplasmic lipid droplets and an opaque yolk, providing a visible phenotype to investigate disturbances of the lipoprotein pathway, already leading to discoveries in MTTP (microsomal triglyceride transfer protein) and ApoB (apolipoprotein B). By 5 days of development, the digestive system is functional, making it possible to study fluorescently labeled lipid uptake in the transparent larvae. These and other approaches enabled the first in vivo description of the STAB (stabilin) receptors, showing lipoprotein uptake in endothelial cells. Various zebrafish models have been developed to mimic human diseases by mutating genes known to influence lipoproteins (eg, ldlra, apoC2). This review aims to discuss the most recent research in the zebrafish ApoB-containing lipoprotein and lipid metabolism field. We also summarize new insights into lipid processing within the yolk cell and how changes in lipid flux alter yolk opacity. This curious new finding, coupled with the development of several techniques, can be deployed to identify new players in lipoprotein research directly relevant to human disease.

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.

Embryonic expression patterns of TBL1 family in zebrafish

Except the addition of TBL1Y in human, transducing beta like 1 (TBL1) family mainly consists of two members TBL1X and TBL1XR1, taking part in multiple intracellular signaling pathways such as Wnt/β-catenin and NF-κB in cancer progression. However, the gene expression patterns of this family during embryonic development remain largely unknown. Here we took advantage of zebrafish model to characterize the spatial and temporal expression patterns of TBL1 family genes including tbl1x, tbl1xr1a and tbl1xr1b. The in situ hybridization studies of gene expression showed robust expressions of tbl1x and tbl1xr1b as maternal transcripts except tbl1xr1a. As the embryo develops, zygotic expressions of all TBL1 family members occur and have a redundant and broad pattern including in brain, neural retina, pharyngeal arches, otic vesicles, and pectoral fins. Ubiquitous expression of all family members were ranked from the strongest to the weakest: tbl1xr1a, tbl1x, and tbl1xr1b. In addition, one tbl1xr1a transcript tbl1xr1a202 showed unique and rich expression in the developing heart and lateral line neuromasts. Overall, all members of zebrafish TBL1 family shared numerous similarities and exhibited certain distinctions in the expression patterns, indicating that they might have redundant and exclusive functions to be further explored.

Hox genes control homocercal caudal fin development and evolution

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

Zebrafish cardiac repolarization does not functionally depend on the expression of the hERG1b-like transcript

Zebrafish provide a translational model of human cardiac function. Their similar cardiac electrophysiology enables screening of human cardiac repolarization disorders, drug arrhythmogenicity, and novel antiarrhythmic therapeutics. However, while zebrafish cardiac repolarization is driven by delayed rectifier potassium channel current (IKr), the relative role of alternate channel transcripts is uncertain. While human ether-a-go-go-related-gene-1a (hERG1a) is the dominant transcript in humans, expression of the functionally distinct alternate transcript, hERG1b, modifies the electrophysiological and pharmacologic IKr phenotype. Studies of zebrafish IKr are frequently translated without consideration for the presence and impact of hERG1b in humans. Here, we performed phylogenetic analyses of all available KCNH genes from Actinopterygii (ray-finned fishes). Our findings confirmed zebrafish cardiac zkcnh6a as the paralog of human hERG1a (hKCNH2a), but also revealed evidence of a hERG1b (hKCNH2b)-like N-terminally truncated gene, zkcnh6b, in zebrafish. zkcnh6b is a teleost-specific variant that resulted from the 3R genome duplication. qRT-PCR showed dominant expression of zkcnh6a in zebrafish atrial and ventricular tissue, with low levels of zkcnh6b. Functional evaluation of zkcnh6b in a heterologous system showed no discernable function under the conditions tested, and no influence on zkcnh6a function during the zebrafish ventricular action potential. Our findings provide the first descriptions of the zkcnh6b gene, and show that, unlike in humans, zebrafish cardiac repolarization does not rely upon co-assembly of zERG1a/zERG1b. Given that hERG1b modifies IKr function and drug binding in humans, our findings highlight the need for consideration when translating hERG variant effects and toxicological screens in zebrafish, which lack a functional hERG1b-equivalent gene.

Hox, homology, and parsimony: An organismal perspective

Hox genes are important regulators in animal development. They often show a mosaic of conserved (e.g., longitudinal axis patterning) and lineage-specific novel functions (e.g., development of skeletal, sensory, or locomotory systems). Despite extensive research over the past decades, it remains controversial at which node in the animal tree of life the Hox cluster evolved. Its presence already in the last common metazoan ancestor has been proposed, although the genomes of both putative earliest extant metazoan offshoots, the ctenophores and the poriferans, are devoid of Hox sequences. The lack of Hox genes in the supposedly "simple"-built poriferans and their low number in cnidarians and the basally branching bilaterians, the xenacoelomorphs, seems to support the classical notion that the number of Hox genes is correlated with the degree of animal complexity. However, the 4-fold increase of the Hox cluster in xiphosurans, a basally branching chelicerate clade, as well as the situation in some teleost fishes that show a multitude of Hox genes compared to, e.g., human, demonstrates, that there is no per se direct correlation between organismal complexity and Hox number. Traditional approaches have tried to base homology on the morphological level on shared expression profiles of individual genes, but recent data have shown that, in particular with respect to Hox and other regulatory genes, complex gene-gene interactions rather than expression signatures of individual genes alone are responsible for shaping morphological traits during ontogeny. Accordingly, for sound homology assessments and reconstructions of character evolution on organ system level, additional independent datasets (e.g., morphological, developmental) need to be included in any such analyses. If supported by solid data, proposed structural homology should be regarded as valid and not be rejected solely on the grounds of non-parsimonious distribution of the character over a given phylogenetic topology.

Behavioral Assays Dissecting NMDA Receptor Function in Zebrafish

Zebrafish are a powerful system to study brain development and to dissect the activity of complex circuits. One advantage is that they display complex behaviors, including prey capture, learning, responses to photic and acoustic stimuli, and social interaction (Dreosti et al., Front Neural Circuits 9:39, 2015; Bruckner et al., PLoS Biol 20:e3001838, 2022; Zoodsma et al., Mol Autism 13:38, 2022) that can be probed to assess brain function. Many of these behaviors are easily assayed at early larval stages, offering a noninvasive and high-throughput readout of nervous system function. Additionally, larval zebrafish readily uptake small molecules dissolved in water making them ideal for behavioral-based drug screens. Together, larval zebrafish and their behavioral repertoire offer a means to rapidly dissect brain circuitry and can serve as a template for high-throughput small molecule screens.NMDA receptor subunits are highly conserved in zebrafish compared to mammals (Zoodsma et al., Mol Autism 13:38, 2022; Cox et al., Dev Dyn 234:756-766, 2005; Zoodsma et al., J Neurosci 40:3631-3645, 2020). High amino acid and domain structure homology between humans and zebrafish underlie conserved functional similarities. Here we describe a set of behavioral assays that are useful to study the NMDA receptor activity in brain function.

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.

The evolution of NLRC3 subfamily genes in Sebastidae teleost fishes

BACKGROUND

With more than 36,000 valid fish species, teleost fishes constitute the most species-rich vertebrate clade and exhibit extensive genetic and phenotypic variation, including diverse immune defense strategies. NLRC3 subfamily genes, which are specific to fishes, play vital roles in the immune system of teleosts. The evolution of teleosts has been impacted by several whole-genome duplication (WGD) events, which might be a key reason for the expansions of the NLRC3 subfamily, but detailed knowledge of NLRC3 subfamily evolution in the family Sebastidae is still limited.

RESULTS

Phylogenetic inference of NLRC3 subfamily protein sequences were conducted to evaluate the orthology of NLRC3 subfamily genes in black rockfish (Sebastes schlegilii), 13 other fish species from the families Sebastidae, Serranidae, Gasterosteidae and Cyclopteridae, and three species of high vertebrates (bird, reptile and amphibian). WGD analyses were used to estimate expansions and contractions of the NLRC3 subfamily, and patterns of expression of NLRC3 subfamily genes in black rockfish following bacterial infections were used to investigate the functional roles of these genes in the traditional and mucosal immune system of the Sebastidae. Different patterns of gene expansions and contractions were observed in 17 fish and other species examined, and one and two whole-genome duplication events were observed in two members of family Sebastidae (black rockfish and honeycomb rockfish, Sebastes umbrosus), respectively. Subsequently, 179 copy numbers of NLRC3 genes were found in black rockfish and 166 in honeycomb rockfish. Phylogenetic analyses corroborated the conservation and evolution of NLRC3 orthologues between Sebastidae and other fish species. Finally, differential expression analyses provided evidence of the immune roles of NLRC3 genes in black rockfish during bacterial infections and gene ontology analysis also indicated other functional roles.

CONCLUSIONS

We hypothesize that NLRC3 genes have evolved a variety of different functions, in addition to their role in the immune response, as a result of whole genome duplication events during teleost diversification. Importantly, this study had underscored the importance of sampling across taxonomic groups, to better understand the evolutionary patterns of the innate immunity system on which complex immunological novelties arose. Moreover, the results in this study could extend current knowledge of the plasticity of the immune system.

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

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

Evolution and phylogenetic distribution of endo-α-mannosidase

While glycans underlie many biological processes, such as protein folding, cell adhesion, and cell-cell recognition, deep evolution of glycosylation machinery remains an understudied topic. N-linked glycosylation is a conserved process in which mannosidases are key trimming enzymes. One of them is the glycoprotein endo-α-1,2-mannosidase which participates in the initial trimming of mannose moieties from an N-linked glycan inside the cis-Golgi. It is unique as the only endo-acting mannosidase found in this organelle. Relatively little is known about its origins and evolutionary history; so far it was reported to occur only in vertebrates. In this work, a taxon-rich bioinformatic survey to unravel the evolutionary history of this enzyme, including all major eukaryotic clades and a wide representation of animals, is presented. The endomannosidase was found to be more widely distributed in animals and other eukaryotes. The protein motif changes in context of the canonical animal enzyme were tracked. Additionally, the data show the two canonical vertebrate endomannosidase genes, MANEA and MANEAL, arose at the second round of the two vertebrate genome duplications and one more vertebrate paralog, CMANEAL, is uncovered. Finally, a framework where N-glycosylation co-evolved with complex multicellularity is described. A better understanding of the evolution of core glycosylation pathways is pivotal to understanding biology of eukaryotes in general, and the Golgi apparatus in particular. This systematic analysis of the endomannosidase evolution is one step toward this goal.

A gene regulatory network combining Pax3/7, Sox10 and Mitf generates diverse pigment cell types in medaka and zebrafish

Neural crest cells generate numerous derivatives, including pigment cells, and are a model for studying how fate specification from multipotent progenitors is controlled. In mammals, the core gene regulatory network for melanocytes (their only pigment cell type) contains three transcription factors, Sox10, Pax3 and Mitf, with the latter considered a master regulator of melanocyte development. In teleosts, which have three to four pigment cell types (melanophores, iridophores and xanthophores, plus leucophores e.g. in medaka), gene regulatory networks governing fate specification are poorly understood, although Mitf function is considered conserved. Here, we show that the regulatory relationships between Sox10, Pax3 and Mitf are conserved in zebrafish, but the role for Mitf is more complex than previously emphasized, affecting xanthophore development too. Similarly, medaka Mitf is necessary for melanophore, xanthophore and leucophore formation. Furthermore, expression patterns and mutant phenotypes of pax3 and pax7 suggest that Pax3 and Pax7 act sequentially, activating mitf expression. Pax7 modulates Mitf function, driving co-expressing cells to differentiate as xanthophores and leucophores rather than melanophores. We propose that pigment cell fate specification should be considered to result from the combinatorial activity of Mitf with other transcription factors.

Temporospatial inhibition of Erk signaling is required for lymphatic valve formation

Intraluminal lymphatic valves (LVs) and lymphovenous valves (LVVs) are critical to ensure the unidirectional flow of lymphatic fluid. Morphological abnormalities in these valves always cause lymph or blood reflux, and result in lymphedema. However, the underlying molecular mechanism of valve development remains poorly understood. We here report the implication of Efnb2-Ephb4-Rasa1 regulated Erk signaling axis in lymphatic valve development with identification of two new valve structures. Dynamic monitoring of phospho-Erk activity indicated that Erk signaling is spatiotemporally inhibited in some lymphatic endothelial cells (LECs) during the valve cell specification. Inhibition of Erk signaling via simultaneous depletion of zygotic erk1 and erk2 or treatment with MEK inhibitor selumetinib causes lymphatic vessel hypoplasia and lymphatic valve hyperplasia, suggesting opposite roles of Erk signaling during these two processes. ephb4b mutants, efnb2a;efnb2b or rasa1a;rasa1b double mutants all have defective LVs and LVVs and exhibit blood reflux into lymphatic vessels with an edema phenotype. Importantly, the valve defects in ephb4b or rasa1a;rasa1b mutants are mitigated with high-level gata2 expression in the presence of MEK inhibitors. Therefore, Efnb2-Ephb4 signaling acts to suppress Erk activation in valve-forming cells to promote valve specification upstream of Rasa1. Not only do our findings reveal a molecular mechanism of lymphatic valve formation, but also provide a basis for the treatment of lymphatic disorders.

Loss of dlx5a/dlx6a Locus Alters Non-Canonical Wnt Signaling and Meckel's Cartilage Morphology

The dlx genes encode transcription factors that establish a proximal-distal polarity within neural crest cells to bestow a regional identity during craniofacial development. The expression regions of dlx paralogs are overlapping yet distinct within the zebrafish pharyngeal arches and may also be involved in progressive morphologic changes and organization of chondrocytes of the face. However, how each dlx paralog of dlx1a, dlx2a, dlx5a and dlx6a affects craniofacial development is still largely unknown. We report here that the average lengths of the Meckel's, palatoquadrate and ceratohyal cartilages in different dlx mutants were altered. Mutants for dlx5a-/- and dlx5i6-/-, where the entire dlx5a/dlx6a locus was deleted, have the shortest lengths for all three structures at 5 days post fertilization (dpf). This phenotype was also observed in 14 dpf larvae. Loss of dlx5i6 also resulted in increased proliferation of neural crest cells and expression of chondrogenic markers. Additionally, altered expression and function of non-canonical Wnt signaling were observed in these mutants suggesting a novel interaction between dlx5i6 locus and non-canonical Wnt pathway regulating ventral cartilage morphogenesis.

Transposon waves at the water-to-land transition

The major transitions in vertebrate evolution are associated with significant genomic reorganizations. In contrast to the evolutionary processes that occurred at the origin of vertebrates or prior to the radiation of teleost fishes, no whole-genome duplication events occurred during the water-to-land transition, and it remains an open question how did genome dynamics contribute to this prominent evolutionary event. Indeed, the recent sequencing of sarcopterygian and amphibian genomes has revealed that the extant lineages immediately preceding and succeeding this transition harbor an exceptional number of transposable elements and it is tempting to speculate that these sequences might have catalyzed the adaptations that enabled vertebrates to venture into land. Here, we review the genome dynamics associated with the major transitions in vertebrate evolution and discuss how the highly repetitive genomic landscapes revealed by recent efforts to characterize the genomes of amphibians and sarcopterygians argue for turbulent genome dynamics occurring before the water-to-land transition and possibly enabling it.

The vertebrate sialylation machinery: structure-function and molecular evolution of GT-29 sialyltransferases

Every eukaryotic cell is covered with a thick layer of complex carbohydrates with essential roles in their social life. In Deuterostoma, sialic acids present at the outermost positions of glycans of glycoconjugates are known to be key players in cellular interactions including host-pathogen interactions. Their negative charge and hydrophilic properties enable their roles in various normal and pathological states and their expression is altered in many diseases including cancers. Sialylation of glycoproteins and glycolipids is orchestrated by the regulated expression of twenty sialyltransferases in human tissues with distinct enzymatic characteristics and preferences for substrates and linkages formed. However, still very little is known on the functional organization of sialyltransferases in the Golgi apparatus and how the sialylation machinery is finely regulated to provide the ad hoc sialome to the cell. This review summarizes current knowledge on sialyltransferases, their structure-function relationships, molecular evolution, and their implications in human biology.

Whole-genome Duplications and the Long-term Evolution of Gene Regulatory Networks in Angiosperms

Angiosperms have a complex history of whole-genome duplications (WGDs), with varying numbers and ages of WGD events across clades. These WGDs have greatly affected the composition of plant genomes due to the biased retention of genes belonging to certain functional categories following their duplication. In particular, regulatory genes and genes encoding proteins that act in multiprotein complexes have been retained in excess following WGD. Here, we inferred protein-protein interaction (PPI) networks and gene regulatory networks (GRNs) for seven well-characterized angiosperm species and explored the impact of both WGD and small-scale duplications (SSDs) in network topology by analyzing changes in frequency of network motifs. We found that PPI networks are enriched in WGD-derived genes associated with dosage-sensitive intricate systems, and strong selection pressures constrain the divergence of WGD-derived genes at the sequence and PPI levels. WGD-derived genes in network motifs are mostly associated with dosage-sensitive processes, such as regulation of transcription and cell cycle, translation, photosynthesis, and carbon metabolism, whereas SSD-derived genes in motifs are associated with response to biotic and abiotic stress. Recent polyploids have higher motif frequencies than ancient polyploids, whereas WGD-derived network motifs tend to be disrupted on the longer term. Our findings demonstrate that both WGD and SSD have contributed to the evolution of angiosperm GRNs, but in different ways, with WGD events likely having a more significant impact on the short-term evolution of polyploids.

Genome-wide identification, evolution and expression analysis of tight junction gene family and the immune roles of claudin5 gene in turbot (Scophthalmus maximus L.)

Tight junction proteins (TJs) are important component proteins that maintaining the structure and function of TJs, connecting to each other to form a TJ complex between cells, maintaining the biological homeostasis of the internal environment. In this study, a total of 103 TJ genes were identified in turbot according to our whole-transcriptome database. Transmembrane TJs were divided into seven subfamilies, including claudin (CLDN), occludin (OCLD), tricellulin (MARVELD2), MARVEL domain containing 3 (MARVELD3), junctional adhesion molecules (JAM), immunoglobulin superfamily member 5 (IGSF5/JAM4), blood vessel epicardial substance (BVEs). Moreover, the majority of homologous pairs of TJ genes showed highly conserved alongside length, exon/intron number and motifs. As for phylogenetic analysis for 103 TJ genes, eight of them have undergone a positive selection and JAMB-like has undergone the most neutral evolution. The expression patterns of several TJ genes showed the lowest expression levels in blood, while the highest expression levels were detected in intestine, gill and skin, which all belong to mucosal tissues. Meanwhile, most examined TJ genes showed down-regulated expression patterns during bacterial infection, while several TJ genes exhibited up-regulated expression patterns at a later stage (24 h). At the same time, several potential candidate genes (such as CLDN-15, CLDN-3, CLDN-12, CLDN-5 and OCLD) were significantly down-regulated, which may indicate their important functions that involved in the regulation of bacterial infection. Currently, there is little research on CLDN5 in the intestine, but it is highly expressed in the intestine and has significant changes in intestinal expression after bacterial infection. Thus, we knocked down CLDN5 by the method of lentiviral infection. The result showed CLDN5 was related to cell migration (wound healing) and apoptosis, and the method of dualluciferasereporterassay showed that the functions of CLDN5 could be regulated by miR-24. The study of TJs may lead to a better understanding of the function of TJs in teleost.

Discovery of Teleost Plasma Kallikrein/Coagulation Factor XI-Like Gene from Channel Catfish (Ictalurus punctatus) and the Evidence that the Protein Encoded by it Acts as a Lectin

Mammalian plasma kallikrein (PK) and coagulation factor XI (fXI) are serine proteases that play in the kinin-kallikrein cascade and in the blood clotting pathway. These proteases share sequence homology and have four apple domains (APDs) and a serine protease domain (SPD) from their N-terminus to C-terminus. No homologs of these proteases are believed to be present in fish species, except for lobe-finned fish. Fish, however, have a unique lectin, named kalliklectin (KL), which is composed of APDs only. In the present study, we found genomic sequences encoding a protein with both APDs and SPD in a few cartilaginous and bony fishes, including the channel catfish Ictalurus punctatus, using bioinformatic analysis. Furthermore, we purified two ~ 70 kDa proteins from the blood plasma of the catfish using mannose-affinity and gel filtration chromatography sequentially. Using de novo sequencing with quadrupole time-of-flight tandem mass spectrometry, several internal amino acid sequences in these proteins were mapped onto possible PK/fXI-like sequences that are thought to be splicing variants. Exploration of APD-containing proteins in the hagfish genome database and phylogenetic analysis suggested that the PK/fXI-like gene originated from hepatocyte growth factor, and that the gene was acquired in a common ancestor of jawed fish. Synteny analysis provided evidence for chromosomal translocation around the PK/fXI-like locus that occurred in the common ancestor of holosteans and teleosts after separation from the lobe-finned fish lineage, or gene duplication into two chromosomes, followed by independent gene losses. This is the first identification of PK/fXI-like proteins in teleosts.

Genomic and Transcriptional Profiles of Kelch-like (klhl) Gene Family in Polyploid Carassius Complex

Genome duplication supplies raw genetic materials and has been thought to be essential for evolutionary innovation and ecological adaptation. Here, we select Kelch-like (klhl) genes to study the evolution of the duplicated genes in the polyploid Carassius complex, including amphidiploid C. auratus and amphitriploid C. gibelio. Phylogenetic, chromosomal location and read coverage analyses indicate that most of Carassius klhl genes exhibit a 2:1 relationship with zebrafish orthologs and confirm two rounds of polyploidy, an allotetraploidy followed by an autotriploidy, occurred during Carassius evolution. The lineage-specific expansion and biased retention/loss of klhl genes are also found in Carassius. Transcriptome analyses across eight adult tissues and seven embryogenesis stages reveal varied expression dominance and divergence between the two species. The expression of klhls in response to Carassius herpesvirus 2 infection shows different expression changes corresponding to distinct herpesvirus resistances in three C. gibelio gynogenetic clones. Finally, we find that most C. gibelio klhl genes possess three alleles except eight genes that have lost one or two alleles due to genome rearrangement. The allele expression bias is prosperous for Cgklhl genes and varies during embryogenesis owning to the sequential expression manner of the alleles. The current study provides global insights into the genomic and transcriptional evolution of duplicated genes in a given superfamily resulting from multiple rounds of polyploidization.

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.

Evolutionary loss of the ß1-adrenergic receptor in salmonids

Whole-genome duplications (WGDs) have been at the heart of the diversification of ß-adrenergic receptors (ß-ARs) in vertebrates. Non-teleost jawed vertebrates typically possess three ß-AR genes: adrb1 (ß1-AR), adrb2 (ß2-AR), and adrb3 (ß3-AR), originating from the ancient 2R (two rounds) WGDs. Teleost fishes, owing to the teleost-specific WGD, have five ancestral adrb paralogs (adrb1, adrb2a, adrb2b, adrb3a and adrb3b). Salmonids are particularly intriguing from an evolutionary perspective as they experienced an additional WGD after separating from other teleosts. Moreover, adrenergic regulation in salmonids, especially rainbow trout, has been intensively studied for decades. However, the repertoire of adrb genes in salmonids has not been yet characterized. An exhaustive genome survey of diverse salmonids, spanning five genera, complemented by phylogenetic sequence analysis, revealed each species has seven adrb paralogs: two adrb2a, two adrb2b, two adrb3a and one adrb3b. Surprisingly, salmonids emerge as the first known jawed vertebrate lineage to lack adrb1. adrb1 is nevertheless highly expressed in the hearts of non-salmonid teleosts, indicating that the wealth of data on adrenergic regulation in salmonids should be generalised to other teleost fishes with caution. It is hypothesised that the loss of adrb1 could have been viable because of the evolutionary radiation of adrb2 and adrb3 genes attributable to the salmonid WGD.

Bovine HOXA11 Gene Identified from RNA-Seq: mRNA Profile Analysis and Genetic Variation Detection Using ME Method and Their Associations with Carcass Traits

The Homeobox A11 (HOXA11) gene regulates limb skeletal development and muscle growth, thus, it was selected as a candidate gene for bovine carcass traits. In this study, we analyzed the mRNA expression level of HOXA11 in various tissues and cells, and determined the genetic variations in the HOXA11 gene, which might be used as molecular markers for cattle breeding. The mRNA expression profiles of HOXA11 in bovine different tissues showed that HOXA11 was highly expressed in both fat and muscle. The gene expression trend of HOXA11 in myoblasts and adipocytes indicated that HOXA11 might be involved in the differentiation of bovine myoblasts and adipocytes. The data in the Ensembl database showed that there are two putative insertion/deletion (InDel) polymorphisms in the bovine HOXA11 gene. The insertion site (rs515880802) was located in the upstream region (NC_037331.1: g. 68853364-68853365) and named as P1-Ins-4-bp, and the deletion site (rs517582703) was located in the intronic region (NC_037331.1: g. 68859510-68859517) and named as P2-Del-8-bp. These polymorphisms within the HOXA11 gene were identified and genotyped by PCR amplification, agarose gel electrophoresis and DNA sequencing in the 640 Shandong Black Cattle Genetic Resource (SDBCGR) population. Moreover, the mutation frequency was very low after detection, so the mathematical expectation (ME) method was used for detection. Statistical analysis demonstrated that P1-Ins-4-bp was significantly correlated with the beef shoulder (p = 0.012) and tongue root (p = 0.004). Meanwhile, P2-Del-8-bp displayed a significant correlation with the back tendon (p = 0.008), money tendon (p = 2.84 × 10^-4), thick flank (p = 0.034), beef shin (p = 9.09 × 10^-7), triangle thick flank (p = 0.04), triangle flank (p = 1.00 × 10^-6), rump (p = 0.018) and small tenderloin (p = 0.043) in the female SDBCGR population. In summary, these outcomes may provide a new perspective for accelerating the molecular breeding of cattle through marker-assisted selection (MAS) strategies.

Efficiently Editing Multiple Duplicated Homeologs and Alleles for Recurrent Polyploids

Research on the evolutionary fate of duplicated genes in recurrent polyploids is scarce due to the difficulties in disentangling the different homeologs and alleles of duplicated genes. This chapter describes the detailed procedures to identify different homeologs and alleles of duplicated genes, to analyze their molecular characteristics, and to reveal their functional divergence by gene editing with CRISPR/Cas9 (clustered regularly interspaced short palindromic repeats/CRISPR-associated system 9). Using the gene editing approach, we efficiently constructed multiple knockout mutant lines with single or simultaneously disrupted different homeologs or alleles in a recurrent polyploid fish, demonstrating its usability for targeting and mutating multiple divergent homeologs and alleles in recurrent duplicated genomes.

Genome-wide identification and characterization of toll-like receptor 5 (TLR5) in fishes

Toll-like receptors 5 (TLR5), a member of the toll-like receptors (TLRs) family, is a class of pattern recognition receptors (PRRs) that recognize pathogen-associated molecular patterns (PAMPs). It responds to vertebrate recognition of bacterial flagellin and participates in innate immune responses. However, genome-wide identification and characterization of TLR5 in fishes have not been investigated. Here, three TLR5M isotypes (TLR5Ma, TLR5Mb1, and TLR5Mb2) and a TLR5S are all extracted from fish genomes on the basis of phylogenetic and synteny analyses. We confirmed that the non-teleost fishes have one TLR5M gene, as well as additional TLR5 genes (TLR5M and TLR5S) in teleost fishes. In addition, some special teleost fishes possess two to three TLR5 genes, which have undergone the fourth whole-genome duplication (WGD). According to our results, we inferred that the diversity of TLR5 genes in fishes seems to be the result of combinations of WGD and gene loss. Furthermore, TLR5 isoforms displayed differences at the flagellin interaction sites and viral binding sites, and showed lineage-specific, which indicated that TLR5 duplicates may generate functional divergence. Bacterial experiments also supported the idea that CiTLR5Ma and CiTLR5Mb are subfunctionalized to sense bacterial flagellin. In summary, our present comparative genomic survey will benefit for further functional investigations of TLR5 genes in fish.

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.

Single-cell transcriptomics of the goldfish retina reveals genetic divergence in the asymmetrically evolved subgenomes after allotetraploidization

The recent whole-genome duplication (WGD) in goldfish (Carassius auratus) approximately 14 million years ago makes it a valuable model for studying gene evolution during the early stages after WGD. We analyzed the transcriptome of the goldfish retina at the level of single-cell (scRNA-seq) and open chromatin regions (scATAC-seq). We identified a group of genes that have undergone dosage selection, accounting for 5% of the total 11,444 ohnolog pairs. We also identified 306 putative sub/neo-functionalized ohnolog pairs that are likely to be under cell-type-specific genetic variation at single-cell resolution. Diversification in the expression patterns of several ohnolog pairs was observed in the retinal cell subpopulations. The single-cell level transcriptome analysis in this study uncovered the early stages of evolution in retinal cell of goldfish after WGD. Our results provide clues for understanding the relationship between the early stages of gene evolution after WGD and the evolution of diverse vertebrate retinal functions.

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.

A Chromosome-level assembly of the Japanese eel genome, insights into gene duplication and chromosomal reorganization

Japanese eels (Anguilla japonica) are commercially important species, harvested extensively for food. Currently, this and related species (American and European eels) are challenging to breed on a commercial basis. As a result, the wild stock is used for aquaculture. Moreover, climate change, habitat loss, water pollution, and altered ocean currents affect eel populations negatively. Accordingly, the International Union for Conservation of Nature lists Japanese eels as endangered and on its red list. Here we presented a high-quality genome assembly for Japanese eels and demonstrated that large chromosome reorganizations occurred in the events of third-round whole-genome duplications (3R-WRDs). Several chromosomal fusions and fissions have reduced the ancestral protochromosomal number of 25 to 19 in the Anguilla lineage. A phylogenetic analysis of the expanded gene families showed that the olfactory receptors (group δ and ζ genes) and voltage-gated Ca2+ channels expanded significantly. Both gene families are crucial for olfaction and neurophysiology. Additional tandem and proximal duplications occurred following 3R-WGD to acquire immune-related genes for an adaptive advantage against various pathogens. The Japanese eel assembly presented here can be used to study other Anguilla species relating to evolution and conservation.

Sex Lethal Gene Manipulates Gonadal Development of Medaka, Oryzias latipes, through Estrogenic Interventions

Germ cells are pivotal for gonadal sexuality maintenance and reproduction. Sex lethal (sxl), the somatic sex determining gene of Drosophila, is the known regulator and initiator of germ cell femininity in invertebrates. However, the role of the Sxl homologue has rarely been investigated in vertebrates. So, we used medaka to clarify the role of sxl in vertebrate gonadogenesis and sexuality and identified two Sxl homologues, i.e., Sxl1a and Sxl1b. We found that sxl1a specifically expresses in the primordial germ cells (PGC), ovary, (early gonia and oocytes), while sxl1b distributions are ubiquitous. An mRNA overexpression of sxl1a accelerated germ cell numbers in 10 DAH XY fish, and sxl1a knockdown (KD), on the other hand, induced PGC mis-migration, aberrant PGC structuring and ultimately caused significant germ cell reduction in XX fish. Using an in vitro promoter analysis and in vivo steroid treatment, we found a strong link between sxl1a and estrogenic germ cell-population maintenance. Further, using sxl1a-KD and erβ2-knockout fish, we determined that sxl1 acts through erβ2 and controls PGC sexuality. Cumulatively, our study highlights the novel role of sxl1a in germ cell maintenance and sexual identity assignment and thus might become a steppingstone to understanding the commonalities of animal sexual development.

A New Hope: A Hermaphroditic Nematode Enables Analysis of a Recent Whole Genome Duplication Event

Whole genome duplication (WGD) is often considered a major driver of evolution that leads to phenotypic novelties. However, the importance of WGD for evolution is still controversial because most documented WGD events occurred anciently and few experimental systems amenable to genetic analysis are available. Here, we report a recent WGD event in the hermaphroditic nematode Allodiplogaster sudhausi and present a comparison with a gonochoristic (male/female) sister species that did not undergo WGD. Self-fertilizing reproduction of A. sudhausi makes it amenable to functional analysis and an ideal system to study WGD events. We document WGD in A. sudhausi through karyotype analysis and whole genome sequencing, the latter of which allowed us to 1) identify functional bias in retention of protein domains and metabolic pathways, 2) show most duplicate genes are under evolutionary constraint, 3) show a link between sequence and expression divergence, and 4) characterize differentially expressed duplicates. We additionally show WGD is associated with increased body size and an abundance of repeat elements (36% of the genome), including a recent expansion of the DNA-hAT/Ac transposon family. Finally, we demonstrate the use of CRISPR/Cas9 to generate mutant knockouts, whereby two WGD-derived duplicate genes display functional redundancy in that they both need to be knocked out to generate a phenotype. Together, we present a novel experimental system that is convenient for examining and characterizing WGD-derived genes both computationally and functionally.

Three IRF4 paralogs act as negative regulators of type Ⅰ IFN responses in yellow catfish (Pelteobagrus fulvidraco)

IRF4 is a master member of the interferon regulatory factor (IRF) family playing vital regulatory roles in immune system development and function. Tetrapods have a single-copy IRF4 gene, while teleosts harbor duplicated IRF4 genes. This work describes three IRF4 paralogs from yellow catfish (Pelteobagrus fulvidraco), designated PfIRF4A, PfIRF4B and PfIRF4B-like. These genes all contain a typical IRF structural architecture. Phylogenic and synteny analyses indicate that they should arise from the teleost-specific whole-genome duplication. PfIRF4 genes are abundantly expressed in the immune-related tissues and upregulated by PolyI:C, LPS, and Edwardsiella ictaluri. Ectopic expression of these genes inhibits the activation of fish type Ⅰ IFN promoters and downregulates the transcription levels of IFN-responsive genes, thus allowing the efficient replication of a fish rhabdovirus, spring viremia of carp virus (SVCV). PfIRF4s possess a repressive effect on MyD88-mediated activation of IFN and NF-κB. Some differences are observed between each individual paralog. PfIRF4B is the main form expressed across the tissues and the most up-regulated one after pathogen induction. It exerts a stronger inhibitory effect on IFN antiviral response than the other two paralogs. PfIRF4A and PfIRF4B-like are primarily present in the nucleus, while PfIRF4B displays colocalization and direct associations with MyD88 in the cytoplasm. Overall, the data demonstrate that three PfIRF4 paralogs show shared and individual functional properties in the negative regulation of type Ⅰ IFN response.

Paralogous Genes Involved in Embryonic Development: Lessons from the Eye and Other Tissues

During evolution, gene duplications lead to a naturally increased gene dosage. Duplicated genes can be further retained or eliminated over time by purifying selection pressure. The retention probability is increased by functional diversification and by the acquisition of novel functions. Interestingly, functionally diverged paralogous genes can maintain a certain level of functional redundancy and at least a partial ability to replace each other. In such cases, diversification probably occurred at the level of transcriptional regulation. Nevertheless, some duplicated genes can maintain functional redundancy after duplication and the ability to functionally compensate for the loss of each other. Many of them are involved in proper embryonic development. The development of particular tissues/organs and developmental processes can be more or less sensitive to the overall gene dosage. Alterations in the gene dosage or a decrease below a threshold level may have dramatic phenotypic consequences or even lead to embryonic lethality. The number of functional alleles of particular paralogous genes and their mutual cooperation and interactions influence the gene dosage, and therefore, these factors play a crucial role in development. This review will discuss individual interactions between paralogous genes and gene dosage sensitivity during development. The eye was used as a model system, but other tissues are also included.

Sequential and additive expression of miR-9 precursors control timing of neurogenesis

MicroRNAs (miRs) have an important role in tuning dynamic gene expression. However, the mechanism by which they are quantitatively controlled is unknown. We show that the amount of mature miR-9, a key regulator of neuronal development, increases during zebrafish neurogenesis in a sharp stepwise manner. We characterize the spatiotemporal profile of seven distinct microRNA primary transcripts (pri-mir)-9s that produce the same mature miR-9 and show that they are sequentially expressed during hindbrain neurogenesis. Expression of late-onset pri-mir-9-1 is added on to, rather than replacing, the expression of early onset pri-mir-9-4 and -9-5 in single cells. CRISPR/Cas9 mutation of the late-onset pri-mir-9-1 prevents the developmental increase of mature miR-9, reduces late neuronal differentiation and fails to downregulate Her6 at late stages. Mathematical modelling shows that an adaptive network containing Her6 is insensitive to linear increases in miR-9 but responds to stepwise increases of miR-9. We suggest that a sharp stepwise increase of mature miR-9 is created by sequential and additive temporal activation of distinct loci. This may be a strategy to overcome adaptation and facilitate a transition of Her6 to a new dynamic regime or steady state.

Site pleiotropy of a stickleback Bmp6 enhancer

Development and regeneration are orchestrated by gene regulatory networks that operate in part through transcriptional enhancers. Although many enhancers are pleiotropic and are active in multiple tissues, little is known about whether enhancer pleiotropy is due to 1) site pleiotropy, in which individual transcription factor binding sites (TFBS) are required for activity in multiple tissues, or 2) multiple distinct sites that regulate expression in different tissues. Here, we investigated the pleiotropy of an intronic enhancer of the stickleback Bone morphogenetic protein 6 (Bmp6) gene. This enhancer was previously shown to regulate evolved changes in tooth number and tooth regeneration, and is highly pleiotropic, with robust activity in both fins and teeth throughout embryonic, larval, and adult life, and in the heart and kidney in adult fish. We tested the hypothesis that the pleiotropy of this enhancer is due to site pleiotropy of an evolutionarily conserved predicted Foxc1 TFBS. Transgenic analysis and site-directed mutagenesis experiments both deleting and scrambling this predicted Foxc1 TFBS revealed that the binding site is required for enhancer activity in both teeth and fins throughout embryonic, larval, and adult development, and in the heart and kidney in adult fish. Collectively these data support a model where the pleiotropy of this Bmp6 enhancer is due to site pleiotropy and this putative binding site is required for enhancer activity in multiple anatomical sites from the embryo to the adult.

Disruption of grin2B, an ASD-associated gene, produces social deficits in zebrafish

BACKGROUND

Autism spectrum disorder (ASD), like many neurodevelopmental disorders, has complex and varied etiologies. Advances in genome sequencing have identified multiple candidate genes associated with ASD, including dozens of missense and nonsense mutations in the NMDAR subunit GluN2B, encoded by GRIN2B. NMDARs are glutamate-gated ion channels with key synaptic functions in excitatory neurotransmission. How alterations in these proteins impact neurodevelopment is poorly understood, in part because knockouts of GluN2B in rodents are lethal.

METHODS

Here, we use CRISPR-Cas9 to generate zebrafish lacking GluN2B (grin2B-/-). Using these fish, we run an array of behavioral tests and perform whole-brain larval imaging to assay developmental roles and functions of GluN2B.

RESULTS

We demonstrate that zebrafish GluN2B displays similar structural and functional properties to human GluN2B. Zebrafish lacking GluN2B (grin2B-/-) surprisingly survive into adulthood. Given the prevalence of social deficits in ASD, we assayed social preference in the grin2B-/- fish. Wild-type fish develop a strong social preference by 3 weeks post fertilization. In contrast, grin2B-/- fish at this age exhibit significantly reduced social preference. Notably, the lack of GluN2B does not result in a broad disruption of neurodevelopment, as grin2B-/- larvae do not show alterations in spontaneous or photic-evoked movements, are capable of prey capture, and exhibit learning. Whole-brain imaging of grin2B-/- larvae revealed reduction of an inhibitory neuron marker in the subpallium, a region linked to ASD in humans, but showed that overall brain size and E/I balance in grin2B-/- is comparable to wild type.

LIMITATIONS

Zebrafish lacking GluN2B, while useful in studying developmental roles of GluN2B, are unlikely to model nuanced functional alterations of human missense mutations that are not complete loss of function. Additionally, detailed mammalian homologies for larval zebrafish brain subdivisions at the age of whole-brain imaging are not fully resolved.

CONCLUSIONS

We demonstrate that zebrafish completely lacking the GluN2B subunit of the NMDAR, unlike rodent models, are viable into adulthood. Notably, they exhibit a highly specific deficit in social behavior. As such, this zebrafish model affords a unique opportunity to study the roles of GluN2B in ASD etiologies and establish a disease-relevant in vivo model for future studies.

HOX-Gene Cluster Organization and Genome Duplications in Fishes and Mammals: Transcript Variant Distribution along the Anterior–Posterior Axis

Hox genes play a crucial role in morphogenesis, especially in anterior–posterior body axis patterning. The organization of Hox clusters in vertebrates is a result of several genome duplications: two rounds of duplication in the ancestors of all vertebrates and a third round that was specific for teleost fishes. Teleostei cluster structure has been significantly modified in the evolutionary processes by Hox gene losses and co-options, while mammals show no such tendency. In mammals, the Hox gene number in a single cluster is stable and generally large, and the numbers are similar to those in the Chondrichthyes. Hox gene alternative splicing activity slightly differs between fishes and mammals. Fishes and mammals have differences in their known alternative splicing activity for Hox gene distribution along the anterior–posterior body axis. The analyzed fish groups—the Coelacanthiformes, Chondrichthyes, and Teleostei—all have higher known alternative mRNA numbers from the anterior and posterior regions, whereas mammals have a more uniform Hox transcript distribution along this axis. In fishes, most Hox transcripts produce functioning proteins, whereas mammals have significantly more known transcripts that do not produce functioning proteins.

Identification, conservation, and expression of tiered pharmacogenes in zebrafish

The number of adverse drug events in the United States is critically high, with annual rates exceeding 1 million cases over the last nine years. One cause of adverse drug events is the underlying genetic variation that can alter drug responses. Pharmacogenomics is a growing field that seeks to better understand the relationship between a patient’s genetics and drug efficacy. Currently, pharmacogenomics relies largely on human trials, as there is not a well-developed animal model for studying preventative measures and alternative treatments. Here, we analyzed pharmacogene expression at two developmental time points in zebrafish to demonstrate the potential of using this model organism for high-throughput pharmacogenomics research. We found that 76% of tiered human pharmacogenes have a zebrafish ortholog, and of these, many have highly conserved amino acid sequences. Additional gene ontology analysis was used to classify pharmacogenes and identify candidate pathways for future modeling in zebrafish. As precision medicine burgeons, adopting a high-throughput in vivo model such as the zebrafish could greatly increase our understanding of the molecular pathology underlying adverse drug events.