Genome duplication, a trait shared by 22000 species of ray-finned fish

Decoding the fish genome opens a new era in important trait research and molecular breeding in China

Aquaculture represents the fastest-growing global food production sector, as it has become an essential component of the global food supply. China has the world's largest aquaculture industry in terms of production volume. However, the sustainable development of fish culture is hindered by several concerns, including germplasm degradation and disease outbreaks. The practice of genomic breeding, which relies heavily on genome information and genotypephenotype relationships, has significant potential for increasing the efficiency of aquaculture production. In 2014, the completion of the genome sequencing and annotation of the Chinese tongue sole signified the beginning of the fish genomics era in China. Since then, domestic researchers have made dramatic progress in functional genomic studies. To date, the genomes of more than 60 species of fish in China have been assembled and annotated. Based on these reference genomes, evolutionary, comparative, and functional genomic studies have revolutionized our understanding of a wide range of biologically and economically important traits of fishes, including growth and development, sex determination, disease resistance, metamorphosis, and pigmentation. Furthermore, genomic tools and breeding techniques such as SNP arrays, genomic selection, and genome editing have greatly accelerated genetic improvement through the incorporation of functional genomic information into breeding activities. This review aims to summarize the current status, advances, and perspectives of the genome resources, genomic study of important traits, and genomic breeding techniques of fish in China. The review will provide aquaculture researchers, fish breeders, and farmers with updated information concerning fish genomic research and breeding technology. The summary will help to promote the genetic improvement of production traits and thus will support the sustainable development of fish aquaculture.

Comparative analysis of Krüppel-like factors expression in the retinas of zebrafish and mice during development and after injury

The Krüppel-like factors (KLFs) have emerged as important transcriptional regulators of various cellular processes, including neural development. Some of them have been described as intrinsic factors involved in axon regeneration in the central nervous system (CNS) of vertebrates. Zebrafish are known for their ability to regenerate several tissues in adulthood, including the CNS, a capability lost during vertebrate evolution and absent in adult mammals. The role that KLFs could play in this differential ability remains unknown. Therefore, in this study, we analyzed the endogenous response of certain KLFs implicated in axon regeneration (KLFs 6, 7, 9, and 13) during retina development and after axon injury. The results showed that the expression of Klfs 6, 7, and 13 decreases in the developing retina of mice but not in zebrafish, while the mRNA levels of Klf9 strongly increase in both species. The response to injury was further analyzed using optic nerve crush (ONC) as a model of lesion. Our analysis during the acute phase (hours) demonstrated an induction of Klfs 6 and 7 expression exclusively in the zebrafish retina, while Klfs 9 and 13 mRNA levels increased in both species. Further analysis of the chronic response (days) showed that mRNA levels of Klf6 transiently increase in the retinas of both zebrafish and mice, whereas those of Klf7 decrease later after optic nerve injury. In addition, the analysis revealed that the expression of Klf9 decreases, while that of Klf13 increases in the retinas of zebrafish in response to optic nerve injury but remains unaltered in mice. Altogether, these findings support the hypothesis that KLFs may play a role in the differential axon regeneration abilities exhibited by fish and mice.

The AMPK and AKT/GSK3β pathways are involved in recombinant proteins fibroblast growth factor 1 (rFGF1 and rFGF1a) improving glycolipid metabolism in rainbow trout (Oncorhynchus mykiss) fed a high carbohydrate diet

Fibroblast growth factor 1 (FGF1) regulates vertebrate cell growth, proliferation and differentiation, and energy metabolism. In this study, we cloned rainbow trout (Oncorhynchus mykiss) fgf1 and fgf1a, prepared their recombinant proteins (rFGF1 and rFGF1a), and described the molecular mechanisms by which they improve glycolipid metabolism in carnivorous fish. A 31-d feeding trial was conducted to investigate whether they could enhance glycolipid metabolism in rainbow trout on high-carbohydrate diets (HCD). A total of 720 rainbow trout (8.9 ± 0.5 g) were equally divided into 4 groups: the chow diet (CD) group injected with PBS, the HCD group injected with PBS, the HCD group injected with rFGF1 (400 ng/g body weight), and the HCD group injected with rFGF1a (400 ng/g body weight). The results showed that short-term HCD had a significant positive effect on the specific growth rate (SGR) of rainbow trout (P < 0.05). However, it led to an increase in crude fat, serum triglyceride (TG) and glucose content, as well as serum glutamic pyruvic transaminase (GPT) and glutamic oxalacetic transaminase (GOT) contents (P < 0.05), suggesting a negative health effect of HCD. Nevertheless, rFGF1 and rFGF1a showed beneficial therapeutic effects. They significantly reduced the crude fat content of the liver, serum TG, GOT, and GPT contents caused by HCD (P < 0.05). The upregulation in atgl, hsl, and acc2 mRNAs implied the promotion of TG catabolism. Moreover, rFGF1 and rFGF1a contributed to promoting lipolysis by activating the AMPK pathway and reducing lipid accumulation in the liver caused by HCD. In addition, the rFGF1 and rFGF1a-treated groups significantly reduced serum glucose levels and elevated hepatic glycogen content under HCD, and increased glucose uptake by hepatocytes. We observed a decrease in mRNA levels for pepck, g6pase, and pygl, along with an increase in mRNA levels for gys, glut2, and gk in the liver. Furthermore, these proteins regulated hepatic gluconeogenesis and glycogen synthesis by increasing the phosphorylation level of AKT, ultimately leading to an increase in GSK3β phosphorylation. In conclusion, this study demonstrates that rFGF1 and rFGF1a can enhance lipolysis and glucose utilization in rainbow trout by activating the AMPK pathway and AKT/GSK3β axis.

The Peroxisome Proliferator-Activated Receptors of Ray-Finned Fish: Unique Structures, Elusive Functions

The Actinopterygian and specifically the Teleostean peroxisome proliferator-activated receptors (PPARs) present an impressive variability and complexity in their structures, both at the gene and protein levels. These structural differences may also reflect functional divergence from their mammalian homologs, or even between fish species. This review, taking advantage of the data generated from the whole-genome sequencing of several fish species, highlights the differences in the primary structure of the receptors, while discussing results from the literature pertaining to the functions of fish PPARs and their activation by natural and synthetic compounds.

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.

Characteristics of CD4-1 gene and its immune responses against Aeromonas veronii infection by activating NF-κB signaling in Qihe crucian carp Carassius auratus

CD4-1 found in bony fish contains four extracellular immunoglobulin (Ig)-like domains similar to that of mammalian CD4, which is crucial for the activation of CD4+ helper T-cell. However, there is limited knowledge regarding the molecular markers, immune functions and regulation mechanism of CD4-1 in teleosts due to their vast diversity. In this study, we cloned and characterized two isoforms of Qihe crucian carp CD4-1, designated as CaCD4-1.1 and CaCD4-1.2. We further explored their expression responses upon stimulation with Aeromonas veronii, and the regulation of their immune responses against A. veronii by NF-κB. The ORF of CaCD4-1.1 and CaCD4-1.2 cDNA encoded 477 and 466 amino acids, respectively. Both proteins contained seven conserved cysteine residues in the extracellular domain, and a CCC motif in their cytoplasm, respectively. However, CaCD4-1.1 exhibited a relatively limited similarity with CaCD4-1.2 in the ectodomain. The quantitative real-time polymerase chain reaction (qRT-PCR) analysis revealed that the mRNA expression of CaCD4-1.1 and CaCD4-1.2 exhibited differential constitutive expression across all examined tissues. Furthermore, the expression level of CD4-1.2 was higher than that of CD4-1.1 in the gills, head kidney, and spleen of Qihe crucian carp subjected to A. veronii challenge, while it was lower in the trunk kidney. Inhibition of NF-κB activity resulted in a decrease in the expression levels of CD4-1.1 and CD4-1.2 mRNA in the gill, while inducing an increase in expression levels in the spleen, in accordance with the observed ultrastructural changes in both organs. Interestingly, the impact of NF-κB on the mRNA expression level of CD4-1.1 appears to be stronger than that of CD4-1.2. Our results suggest that CaCD4-1.1 and CaCD4-1.2 could be expressed on T cells and antigen-sampling cells that exhibit similar characteristics to mammalian M cells, respectively, and differentially regulated by NF-κB in adaptive immune responses against bacterial infection. This research contributes to a better understanding of the crucial role of CD4-1 in the immune response of Qihe crucian carp and provide novel insights for the prevention and treatment of fish diseases in aquaculture.

CRK and NCK adaptors may functionally overlap in zebrafish neurodevelopment, as indicated by common binding partners and overlapping expression patterns

CRK adaptor proteins are important for signal transduction mechanisms driving cell proliferation and positioning during vertebrate central nervous system development. Zebrafish lacking both CRK family members exhibit small, disorganized retinas with 50% penetrance. The goal of this study was to determine whether another adaptor protein might functionally compensate for the loss of CRK adaptors. Expression patterns in developing zebrafish, and bioinformatic analyses of the motifs recognized by their SH2 and SH3 domains, suggest NCK adaptors are well-positioned to compensate for loss of CRK adaptors. In support of this hypothesis, proteomic analyses found CRK and NCK adaptors share overlapping interacting partners including known regulators of cell adhesion and migration, suggesting their functional intersection in neurodevelopment.

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.

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.

Assembly of the poorly differentiated Verasper variegatus W chromosome by different sequencing technologies

The assembly of W and Y chromosomes poses significant challenges in vertebrate genome sequencing and assembly. Here, we successfully assembled the W chromosome of Verasper variegatus with a length of 20.48 Mb by combining population and PacBio HiFi sequencing data. It was identified as a young sex chromosome and showed signs of expansion in repetitive sequences. The major component of the expansion was Ty3/Gypsy. The ancestral Osteichthyes karyotype consists of 24 protochromosomes. The sex chromosomes in four Pleuronectiformes species derived from a pair of homologous protochromosomes resulting from a whole-genome duplication event in teleost fish, yet with different sex-determination systems. V. variegatus and Cynoglossus semilaevis adhere to the ZZ/ZW system, while Hippoglossus stenolepis and H. hippoglossus follow the XX/XY system. Interestingly, V. variegatus and H. hippoglossus derived from one protochromosome, while C. semilaevis and H. stenolepis derived from another protochromosome. Our study provides valuable insights into the evolution of sex chromosomes in flatfish and sheds light on the important role of whole-genome duplication in shaping the evolution of sex chromosomes.

The zebrafish paralog six2b is required for early proximal pronephros morphogenesis

The transcription factor Six2 plays a crucial role in maintaining self-renewing nephron progenitor cap mesenchyme (CM) during metanephric kidney development. In mouse and human, expression at single-cell resolution has detected Six2 in cells as they leave the CM pool and differentiate. The role Six2 may play in these cells as they differentiate remains unknown. Here, we took advantage of the zebrafish pronephric kidney which forms directly from intermediate mesoderm to test six2b function during pronephric tubule development and differentiation. Expression of six2b during early zebrafish development was consistent with a role in pronephros formation. Using morpholino knock-down and CRISPR/Cas9 mutagenesis, we show a functional role for six2b in the development of proximal elements of the pronephros. By 48 h post-fertilization, six2b morphants and mutants showed disrupted pronephric tubule morphogenesis. We observed a lower-than-expected frequency of phenotypes in six2b stable genetic mutants suggesting compensation. Supporting this, we detected increased expression of six2a in six2b stable mutant embryos. To further confirm six2b function, F0 crispant embryos were analyzed and displayed similar phenotypes as morphants and stable mutants. Together our data suggests a conserved role for Six2 during nephrogenesis and a role in the morphogenesis of the proximal tubule.

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.

Gene expression in notochord and nuclei pulposi: a study of gene families across the chordate phylum

The transition from notochord to vertebral column is a crucial milestone in chordate evolution and in prenatal development of all vertebrates. As ossification of the vertebral bodies proceeds, involutions of residual notochord cells into the intervertebral discs form the nuclei pulposi, shock-absorbing structures that confer flexibility to the spine. Numerous studies have outlined the developmental and evolutionary relationship between notochord and nuclei pulposi. However, the knowledge of the similarities and differences in the genetic repertoires of these two structures remains limited, also because comparative studies of notochord and nuclei pulposi across chordates are complicated by the gene/genome duplication events that led to extant vertebrates. Here we show the results of a pilot study aimed at bridging the information on these two structures. We have followed in different vertebrates the evolutionary trajectory of notochord genes identified in the invertebrate chordate Ciona, and we have evaluated the extent of conservation of their expression in notochord cells. Our results have uncovered evolutionarily conserved markers of both notochord development and aging/degeneration of the nuclei pulposi.

Stanniocalcin 1a regulates organismal calcium balance and survival by suppressing Trpv6 expression and inhibiting IGF signaling in zebrafish

Stanniocalcin 1 (Stc1) is well known for its role in regulating calcium uptake in fish by acting on ionocytes or NaR cells. A hallmark of NaR cells is the expression of Trpv6, a constitutively open calcium channel. Recent studies in zebrafish suggest that genetical deletion of Stc1a and Trpv6 individually both increases IGF signaling and NaR cell proliferation. While trpv6-/- fish suffered from calcium deficiency and died prematurely, stc1a-/- fish had elevated body calcium levels but also died prematurely. The relationship between Stc1a, Trpv6, and IGF signaling in regulating calcium homeostasis and organismal survival is unclear. Here we report that loss of Stc1a increases Trpv6 expression in NaR cells in an IGF signaling-dependent manner. Treatment with CdCl2, a Trpv6 inhibitor, reduced NaR cell number in stc1a -/- fish to the sibling levels. Genetic and biochemical analysis results suggest that Stc1a and Trpv6 regulate NaR cell proliferation via the same IGF pathway. Alizarin red staining detected abnormal calcium deposits in the yolk sac region and kidney stone-like structures in stc1a -/- fish. Double knockout or pharmacological inhibition of Trpv6 alleviated these phenotypes, suggesting that Stc1a inhibit epithelial Ca2+ uptake by regulating Trpv6 expression and activity. stc1a-/- mutant fish developed cardiac edema, body swelling, and died prematurely. Treatment of stc1a-/- fish with CdCl2 or double knockout of Trpv6 alleviated these phenotypes. These results provide evidence that Stc1a regulates calcium homeostasis and organismal survival by suppressing Trpv6 expression and inhibiting IGF signaling in ionocytes.

The Himalayan uplift and evolution of aquatic biodiversity across Asia: Snowtrout (Cyprininae: Schizothorax) as a test case

Global biodiversity hotspots are often remote, tectonically active areas undergoing climatic fluctuations, such as the Himalaya Mountains and neighboring Qinghai-Tibetan Plateau (QTP). They provide biogeographic templates upon which endemic biodiversity can be mapped to infer diversification scenarios. Yet, this process can be somewhat opaque for the Himalaya, given substantial data gaps separating eastern and western regions. To help clarify, we evaluated phylogeographic and phylogenetic hypotheses for a widespread fish (Snowtrout: Cyprininae; Schizothorax) by sequencing 1,140 base pair of mtDNA cytochrome-b (cytb) from Central Himalaya samples (Nepal: N = 53; Bhutan: N = 19), augmented with 68 GenBank sequences (N = 60 Schizothorax/N = 8 outgroups). Genealogical relationships (N = 132) were analyzed via maximum likelihood (ML), Bayesian (BA), and haplotype network clustering, with clade divergence estimated via TimeTree. Snowtrout seemingly originated in Central Asia, dispersed across the QTP, then into Bhutan via southward-flowing tributaries of the east-flowing Yarlung-Tsangpo River (YLTR). Headwaters of five large Asian rivers provided dispersal corridors from Central into eastern/southeastern Asia. South of the Himalaya, the YLTR transitions into the Brahmaputra River, facilitating successive westward colonization of Himalayan drainages first in Bhutan, then Nepal, followed by far-western drainages subsequently captured by the (now) westward-flowing Indus River. Two distinct Bhutanese phylogenetic groups were recovered: Bhutan-1 (with three subclades) seemingly represents southward dispersal from the QTP; Bhutan-2 apparently illustrates northward colonization from the Lower Brahmaputra. The close phylogenetic/phylogeographic relationships between the Indus River (Pakistan) and western tributaries of the Upper Ganges (India/Nepal) potentially implicate an historic, now disjunct connection. Greater species-divergences occurred across rather than within-basins, suggesting vicariance as a driver. The Himalaya is a component of the Earth's largest glacial reservoir (i.e., the "third-pole") separate from the Arctic/Antarctic. Its unique aquatic biodiversity must be defined and conserved through broad, trans-national collaborations. Our study provides an initial baseline for this process.

Characterising the Physiological Responses of Chinook Salmon (Oncorhynchus tshawytscha) Subjected to Heat and Oxygen Stress

In New Zealand, during the hottest periods of the year, some salmon farms in the Marlborough Sounds reach water temperatures above the optimal range for Chinook salmon. High levels of mortality are recorded during these periods, emphasising the importance of understanding thermal stress in this species. In this study, the responses of Chinook salmon (Oncorhynchus tshawytscha) to chronic, long-term changes in temperature and dissolved oxygen were investigated. This is a unique investigation due to the duration of the stress events the fish were exposed to. Health and haematological parameters were analysed alongside gene expression results to determine the effects of thermal stress on Chinook salmon. Six copies of heat shock protein 90 (HSP90) were discovered and characterised: HSP90AA1.1a, HSP90AA1.2a, HSP90AA1.1b, HSP90AA1.2b, HSP90AB1a and HSP90AB1b, as well as two copies of SOD1, named SOD1a and SOD1b. The amino acid sequences contained features similar to those found in other vertebrate HSP90 and SOD1 sequences, and the phylogenetic tree and synteny analysis provided conclusive evidence of their relationship to other vertebrate HSP90 and SOD1 genes. Primers were designed for qPCR to enable the expression of all copies of HSP90 and SOD1 to be analysed. The expression studies showed that HSP90 and SOD1 were downregulated in the liver and spleen in response to longer term exposure to high temperatures and lower dissolved oxygen. HSP90 was also downregulated in the gill; however, the results for SOD1 expression in the gill were not conclusive. This study provides important insights into the physiological and genetic responses of Chinook salmon to temperature and oxygen stress, which are critical for developing sustainable fish aquaculture in an era of changing global climates.

The duplication of genomes and genetic networks and its potential for evolutionary adaptation and survival during environmental turmoil

The importance of whole-genome duplication (WGD) for evolution is controversial. Whereas some view WGD mainly as detrimental and an evolutionary dead end, there is growing evidence that polyploidization can help overcome environmental change, stressful conditions, or periods of extinction. However, despite much research, the mechanistic underpinnings of why and how polyploids might be able to outcompete or outlive nonpolyploids at times of environmental upheaval remain elusive, especially for autopolyploids, in which heterosis effects are limited. On the longer term, WGD might increase both mutational and environmental robustness due to redundancy and increased genetic variation, but on the short-or even immediate-term, selective advantages of WGDs are harder to explain. Here, by duplicating artificially generated Gene Regulatory Networks (GRNs), we show that duplicated GRNs-and thus duplicated genomes-show higher signal output variation than nonduplicated GRNs. This increased variation leads to niche expansion and can provide polyploid populations with substantial advantages to survive environmental turmoil. In contrast, under stable environments, GRNs might be maladaptive to changes, a phenomenon that is exacerbated in duplicated GRNs. We believe that these results provide insights into how genome duplication and (auto)polyploidy might help organisms to adapt quickly to novel conditions and to survive ecological uproar or even cataclysmic events.

Chromosome-level genome of butterflyfish unveils genomic features of unique colour patterns and morphological traits

Chaetodontidae, known as butterflyfishes, are typical fish in coral ecosystems, exhibiting remarkable interspecific differences including body colour patterns and feeding ecology. In this study, we report genomes of three butterflyfish species (Chelmon rostratus, Chaetodon trifasciatus and Chaetodon auriga) and a closely related species from the Pomacanthidae family, Centropyge bicolour, with an average genome size of 65,611 Mb. Chelmon rostratus, comprising 24 chromosomes assembled to the chromosome level, could be served as a reference genome for butterflyfish. By conducting a collinearity analysis between butterflyfishes and several fishes, we elucidated the specific and conserved genomic features of butterflyfish, with particular emphasis on novel genes arising from tandem duplications and their potential functions. In addition to the two melanocyte-specific tyr genes commonly found in fish, we found the gene tyrp3, a new tyrosinase-related proteins gene in the reef fish, including butterflyfish and clownfish, implicating their involvement in the pigmentation diversity of fish. Additionally, we observed a tandem duplication expansion of three copies of nell1 gene in C. rostratus genome, which likely contribute to its unique jaw development and distinctive morphology of its sharp mouth. These results provided valuable genomic resources for further investigations into the genetic diversity and evolutionary adaptations of reef fish.

Molecular Evolution of Aryl Hydrocarbon Receptor Signaling Pathway Genes

The Aryl hydrocarbon receptor is an ancient transcriptional factor originally discovered as a sensor of dioxin. In addition to its function as a receptor of environmental toxicants, it plays an important role in development. Although a significant amount of research has been carried out to understand the AHR signal transduction pathway and its involvement in species' susceptibility to environmental toxicants, none of them to date has comprehensively studied its evolutionary origins. Studying the evolutionary origins of molecules can inform ancestral relationships of genes. The vertebrate genome has been shaped by two rounds of whole-genome duplications (WGD) at the base of vertebrate evolution approximately 600 million years ago, followed by lineage-specific gene losses, which often complicate the assignment of orthology. It is crucial to understand the evolutionary origins of this transcription factor and its partners, to distinguish orthologs from ancient non-orthologous homologs. In this study, we have investigated the evolutionary origins of proteins involved in the AHR pathway. Our results provide evidence of gene loss and duplications, crucial for understanding the functional connectivity of humans and model species. Multiple studies have shown that 2R-ohnologs (genes and proteins that have survived from the 2R-WGD) are enriched in signaling components relevant to developmental disorders and cancer. Our findings provide a link between the AHR pathway's evolutionary trajectory and its potential mechanistic involvement in pathogenesis.

Evolution and divergence of teleost adrenergic receptors: why sometimes 'the drugs don't work' in fish

Adrenaline and noradrenaline, released as hormones and/or neurotransmitters, exert diverse physiological functions in vertebrates, and teleost fishes are widely used as model organisms to study adrenergic regulation; however, such investigations often rely on receptor subtype-specific pharmacological agents (agonists and antagonists; see Glossary) developed and validated in mammals. Meanwhile, evolutionary (phylogenetic and comparative genomic) studies have begun to unravel the diversification of adrenergic receptors (ARs) and reveal that whole-genome duplications and pseudogenization events in fishes results in notable distinctions from mammals in their genomic repertoire of ARs, while lineage-specific gene losses within teleosts have generated significant interspecific variability. In this Review, we visit the evolutionary history of ARs (including α1-, α2- and β-ARs) to highlight the prominent interspecific differences in teleosts, as well as between teleosts and other vertebrates. We also show that structural modelling of teleost ARs predicts differences in ligand binding affinity compared with mammalian orthologs. To emphasize the difficulty of studying the roles of different AR subtypes in fish, we collate examples from the literature of fish ARs behaving atypically compared with standard mammalian pharmacology. Thereafter, we focus on specific case studies of the liver, heart and red blood cells, where our understanding of AR expression has benefited from combining pharmacological approaches with molecular genetics. Finally, we briefly discuss the ongoing advances in 'omics' technologies that, alongside classical pharmacology, will provide abundant opportunities to further explore adrenergic signalling in teleosts.

Emergence of the Synucleins

This study establishes the origin and evolutionary history of the synuclein genes. A combination of phylogenetic analyses of the synucleins from twenty-two model species, characterization of local synteny similarities among humans, sharks and lampreys, and statistical comparisons among lamprey and human chromosomes, provides conclusive evidence for the current diversity of synuclein genes arising from the whole-genome duplications (WGDs) that occurred in vertebrates. An ancestral synuclein gene was duplicated in a first WGD, predating the diversification of all living vertebrates. The two resulting genes are still present in agnathan vertebrates. The second WGD, specific to the gnathostome lineage, led to the emergence of the three classical synuclein genes, SNCA, SNCB and SNCG, which are present in all jawed vertebrate lineages. Additional WGDs have added new genes in both agnathans and gnathostomes, while some gene losses have occurred in particular species. The emergence of synucleins through WGDs prevented these genes from experiencing dosage effects, thus avoiding the potential detrimental effects associated with individual duplications of genes that encode proteins prone to aggregation. Additional insights into the structural and functional features of synucleins are gained through the analysis of the highly divergent synuclein proteins present in chondrichthyans and agnathans.

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.

Draft Genome of White-blotched River Stingray Provides Novel Clues for Niche Adaptation and Skeleton Formation

The white-blotched river stingray (Potamotrygon leopoldi) is a cartilaginous fish native to the Xingu River, a tributary of the Amazon River system. As a rare freshwater-dwelling cartilaginous fish in the Potamotrygonidae family in which no member has the genome sequencing information available, P. leopoldi provides the evolutionary details in fish phylogeny, niche adaptation, and skeleton formation. In this study, we present its draft genome of 4.11 Gb comprising 16,227 contigs and 13,238 scaffolds, with contig N50 of 3937 kb and scaffold N50 of 5675 kb in size. Our analysis shows that P. leopoldi is a slow-evolving fish that diverged from elephant sharks about 96 million years ago. Moreover, two gene families related to the immune system (immunoglobulin heavy constant delta genes and T-cell receptor alpha/delta variable genes) exhibit expansion in P. leopoldi only. We also identified the Hox gene clusters in P. leopoldi and discovered that seven Hox genes shared by five representative fish species are missing in P. leopoldi. The RNA sequencing data from P. leopoldi and other three fish species demonstrate that fishes have a more diversified tissue expression spectrum when compared to mammals. Our functional studies suggest that lack of the gc gene encoding vitamin D-binding protein in cartilaginous fishes (both P. leopoldi and Callorhinchus milii) could partly explain the absence of hard bone in their endoskeleton. Overall, this genome resource provides new insights into the niche adaptation, body plan, and skeleton formation of P. leopoldi, as well as the genome evolution in cartilaginous fishes.

phox2ba: The Potential Genetic Link behind the Overlap in the Symptomatology between CHARGE and Central Congenital Hypoventilation Syndromes

CHARGE syndrome typically results from mutations in the gene encoding chromodomain helicase DNA-binding protein 7 (CHD7). CHD7 is involved in regulating neural crest development, which gives rise to tissues of the skull/face and the autonomic nervous system (ANS). Individuals with CHARGE syndrome are frequently born with anomalies requiring multiple surgeries and often experience adverse events post-anesthesia, including oxygen desaturations, decreased respiratory rates, and heart rate abnormalities. Central congenital hypoventilation syndrome (CCHS) affects ANS components that regulate breathing. Its hallmark feature is hypoventilation during sleep, clinically resembling observations in anesthetized CHARGE patients. Loss of PHOX2B (paired-like homeobox 2b) underlies CCHS. Employing a chd7-null zebrafish model, we investigated physiologic responses to anesthesia and compared these to loss of phox2b. Heart rates were lower in chd7 mutants compared to the wild-type. Exposure to tricaine, a zebrafish anesthetic/muscle relaxant, revealed that chd7 mutants took longer to become anesthetized, with higher respiratory rates during recovery. chd7 mutant larvae demonstrated unique phox2ba expression patterns. The knockdown of phox2ba reduced larval heart rates similar to chd7 mutants. chd7 mutant fish are a valuable preclinical model to investigate anesthesia in CHARGE syndrome and reveal a novel functional link between CHARGE syndrome and CCHS.

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.

Duox is the primary NADPH oxidase responsible for ROS production during adult caudal fin regeneration in zebrafish

Sustained elevated levels of reactive oxygen species (ROS) have been shown to be essential for regeneration in many organisms. This has been shown primarily via the use of pharmacological inhibitors targeting the family of NADPH oxidases (NOXes). To identify the specific NOXes involved in ROS production during adult caudal fin regeneration in zebrafish, we generated nox mutants for duox, nox5 and cyba (a key subunit of NOXes 1-4) and crossed these lines with a transgenic line ubiquitously expressing HyPer, which permits the measurement of ROS levels. Homozygous duox mutants had the greatest effect on ROS levels and rate of fin regeneration among the single mutants. However, duox:cyba double mutants showed a greater effect on fin regeneration than the single duox mutants, suggesting that Nox1-4 also play a role during regeneration. This work also serendipitously found that ROS levels in amputated adult zebrafish fins oscillate with a circadian rhythm.

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.

Comparative genomics of Hox and ParaHox genes among major lineages of Branchiopoda with emphasis on tadpole shrimps

Hox and ParaHox genes (HPHGs) are key developmental genes that pattern regional identity along the anterior–posterior body axis of most animals. Here, we identified HPHGs in tadpole shrimps (Pancrustacea, Branchiopoda, Notostraca), an iconic example of the so-called “living fossils” and performed a comparative genomics analysis of HPHGs and the Hox cluster among major branchiopod lineages. Notostraca possess the entire Hox complement, and the Hox cluster seems to be split into two different subclusters, although we were not able to support this finding with chromosome-level assemblies. However, the genomic structure of Hox genes in Notostraca appears more derived than that of Daphnia spp., which instead retains the plesiomorphic condition of a single compact cluster. Spinicaudata and Artemia franciscana show instead a Hox cluster subdivided across two or more genomic scaffolds with some orthologs either duplicated or missing. Yet, branchiopod HPHGs are similar among the various clades in terms of both intron length and number, as well as in their pattern of molecular evolution. Sequence substitution rates are in fact generally similar for most of the branchiopod Hox genes and the few differences we found cannot be traced back to natural selection, as they are not associated with any signals of diversifying selection or substantial switches in selective modes. Altogether, these findings do not support a significant stasis in the Notostraca Hox cluster and further confirm how morphological evolution is not tightly associated with genome dynamics.

Structure and evolution of neuronal wiring receptors and ligands

One of the fundamental properties of a neuronal circuit is the map of its connections. The cellular and developmental processes that allow for the growth of axons and dendrites, selection of synaptic targets, and formation of functional synapses use neuronal surface receptors and their interactions with other surface receptors, secreted ligands, and matrix molecules. Spatiotemporal regulation of the expression of these receptors and cues allows for specificity in the developmental pathways that wire stereotyped circuits. The families of molecules controlling axon guidance and synapse formation are generally conserved across animals, with some important exceptions, which have consequences for neuronal connectivity. Here, we summarize the distribution of such molecules across multiple taxa, with a focus on model organisms, evolutionary processes that led to the multitude of such molecules, and functional consequences for the diversification or loss of these receptors.

The first crystal structure of CD8αα from a cartilaginous fish

Introduction

Cartilaginous fishes are the most evolutionary-distant vertebrates from mammals and possess an immunoglobulin (Ig)- and T cell-mediated adaptive immunity. CD8 is the hallmark receptor of cytotoxic T cells and is required for the formation of T cell receptor-major histocompatibility complex (TCR-MHC) class I complexes.

Methods

RACE PCR was used to obtain gene sequences. Direct dilution was applied for the refolding of denatured recombinant CD8 protein. Hanging-drop vapor diffusion method was performed for protein crystallization.

Results

In this study, CD8α and CD8β orthologues (termed ScCD8α and ScCD8β) were identified in small-spotted catshark (Scyliorhinus canicula). Both ScCD8α and ScCD8β possess an extracellular immunoglobulin superfamily (IgSF) V domain as in previously identified CD8 proteins. The genes encoding CD8α and CD8β are tandemly linked in the genomes of all jawed vertebrates studied, suggesting that they were duplicated from a common ancestral gene before the divergence of cartilaginous fishes and other vertebrates. We determined the crystal structure of the ScCD8α ectodomain homodimer at a resolution of 1.35 Å and show that it exhibits the typical topological structure of CD8α from endotherms. As in mammals, the homodimer formation of ScCD8αα relies upon interactions within a hydrophobic core although this differs in position and amino acid composition. Importantly, ScCD8αα shares the canonical cavity required for interaction with peptide-loaded MHC I in mammals. Furthermore, it was found that ScCD8α can co-immunoprecipitate with ScCD8β, indicating that it can form both homodimeric and heterodimeric complexes.

Conclusion

Our results expand the current knowledge of vertebrate CD8 dimerization and the interaction between CD8α with p/MHC I from an evolutionary perspective.

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.

The evolution of vitamin C biosynthesis and transport in animals

Background

Vitamin C (VC) is an indispensable antioxidant and co-factor for optimal function and development of eukaryotic cells. In animals, VC can be synthesized by the organism, acquired through the diet, or both. In the single VC synthesis pathway described in animals, the penultimate step is catalysed by Regucalcin, and the last step by l-gulonolactone oxidase (GULO). The GULO gene has been implicated in VC synthesis only, while Regucalcin has been shown to have multiple functions in mammals.

Results

Both GULO and Regucalcin can be found in non-bilaterian, protostome and deuterostome species. Regucalcin, as here shown, is involved in multiple functions such as VC synthesis, calcium homeostasis, and the oxidative stress response in both Deuterostomes and Protostomes, and in insects in receptor-mediated uptake of hexamerin storage proteins from haemolymph. In Insecta and Nematoda, however, there is no GULO gene, and in the latter no Regucalcin gene, but species from these lineages are still able to synthesize VC, implying at least one novel synthesis pathway. In vertebrates, SVCT1, a gene that belongs to a family with up to five members, as here shown, is the only gene involved in the uptake of VC in the gut. This specificity is likely the result of a subfunctionalization event that happened at the base of the Craniata subphylum. SVCT-like genes present in non-Vertebrate animals are likely involved in both VC and nucleobase transport. It is also shown that in lineages where GULO has been lost, SVCT1 is now an essential gene, while in lineages where SVCT1 gene has been lost, GULO is now an essential gene.

Conclusions

The simultaneous study, for the first time, of GULO, Regucalcin and SVCTs evolution provides a clear picture of VC synthesis/acquisition and reveals very different selective pressures in different animal taxonomic groups.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02040-7.

Commentary: Zebrafish as a Model for Osteoporosis-An Approach to Accelerating Progress in Drug and Exercise-Based Treatment

Osteoporosis (OP) is a degenerative disease characterized by reduced bone strength and increased fracture risk. As the global population continues to age, the prevalence and economic burden of osteoporosis can be expected to rise substantially, but there remain various gaps in the field of OP care. For instance, there is a lack of anti-fracture drugs with proven long-term efficacy. Likewise, though exercise remains widely recommended in OP prevention and management, data regarding the safety and efficacy for patients after vertebral fracture remain limited. This lack of evidence may be due to the cost and inherent difficulties associated with exercise-based OP research. Thus, the current research landscape highlights the need for novel research strategies that accelerate OP drug discovery and allow for the low-cost study of exercise interventions. Here, we outline an example of one strategy, the use of zebrafish, which has emerged as a potential model for the discovery of anti-osteoporosis therapeutics and study of exercise interventions. The strengths, limitations, and potential applications of zebrafish in OP research will be outlined.

Mechanisms of Na+ uptake from freshwater habitats in animals

Life in fresh water is osmotically and energetically challenging for living organisms, requiring increases in ion uptake from dilute environments. However, mechanisms of ion uptake from freshwater environments are still poorly understood and controversial, especially in arthropods, for which several hypothetical models have been proposed based on incomplete data. One compelling model involves the proton pump V-type H+ ATPase (VHA), which energizes the apical membrane, enabling the uptake of Na+ (and other cations) via an unknown Na+ transporter (referred to as the "Wieczorek Exchanger" in insects). What evidence exists for this model of ion uptake and what is this mystery exchanger or channel that cooperates with VHA? We present results from studies that explore this question in crustaceans, insects, and teleost fish. We argue that the Na+/H+ antiporter (NHA) is a likely candidate for the Wieczorek Exchanger in many crustaceans and insects; although, there is no evidence that this is the case for fish. NHA was discovered relatively recently in animals and its functions have not been well characterized. Teleost fish exhibit redundancy of Na+ uptake pathways at the gill level, performed by different ion transporter paralogs in diverse cell types, apparently enabling tolerance of low environmental salinity and various pH levels. We argue that much more research is needed on overall mechanisms of ion uptake from freshwater habitats, especially on NHA and other potential Wieczorek Exchangers. Such insights gained would contribute greatly to our general understanding of ionic regulation in diverse species across habitats.

An atlas of fish genome evolution reveals delayed rediploidization following the teleost whole-genome duplication

Teleost fishes are ancient tetraploids descended from an ancestral whole-genome duplication that may have contributed to the impressive diversification of this clade. Whole-genome duplications can occur via self-doubling (autopolyploidy) or via hybridization between different species (allopolyploidy). The mode of tetraploidization conditions evolutionary processes by which duplicated genomes return to diploid meiotic pairing, and subsequent genetic divergence of duplicated genes (cytological and genetic rediploidization). How teleosts became tetraploid remains unresolved, leaving a fundamental gap in the interpretation of their functional evolution. As a result of the whole-genome duplication, identifying orthologous and paralogous genomic regions across teleosts is challenging, hindering genome-wide investigations into their polyploid history. Here, we combine tailored gene phylogeny methodology together with a state-of-the-art ancestral karyotype reconstruction to establish the first high-resolution comparative atlas of paleopolyploid regions across 74 teleost genomes. We then leverage this atlas to investigate how rediploidization occurred in teleosts at the genome-wide level. We uncover that some duplicated regions maintained tetraploidy for more than 60 million years, with three chromosome pairs diverging genetically only after the separation of major teleost families. This evidence suggests that the teleost ancestor was an autopolyploid. Further, we find evidence for biased gene retention along several duplicated chromosomes, contradicting current paradigms that asymmetrical evolution is specific to allopolyploids. Altogether, our results offer novel insights into genome evolutionary dynamics following ancient polyploidizations in vertebrates.

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.

The evolution of the testis transcriptome in pregnant male pipefishes and seahorses

In many animals, sperm competition and sexual conflict are thought to drive the rapid evolution of male-specific genes, especially those expressed in the testes. A potential exception occurs in the male pregnant pipefishes, where females transfer eggs to the males, eliminating testes from participating in these processes. Here, we show that testis-related genes differ dramatically in their rates of molecular evolution and expression patterns in pipefishes and seahorses (Syngnathidae) compared to other fish. Genes involved in testis or sperm function within syngnathids experience weaker selection in comparison to their orthologs in spawning and livebearing fishes. An assessment of gene turnover and expression in the testis transcriptome suggests that syngnathids have lost (or significantly reduced expression of) important classes of genes from their testis transcriptomes compared to other fish. Our results indicate that more than 50 million years of male pregnancy have removed syngnathid testes from the molecular arms race that drives the rapid evolution of male reproductive genes in other taxa.

Zebrafish optineurin: genomic organization and transcription regulation

Optineurin (OPTN) is involved in a variety of mechanisms such as autophagy, vesicle trafficking, and NF-κB signaling. Mutations in the OPTN gene have been associated with different pathologies including glaucoma, amyotrophic lateral sclerosis and Paget's disease of bone. Since the relationship between fish and mammalian OPTN is not well understood the objective of the present work was to characterize the zebrafish optn gene and protein structure and to investigate its transcriptional regulation. Through a comparative in silico analysis, we observed that zebrafish optn presents genomic features similar to those of its human counterpart, including its neighboring genes and structure. A comparison of OPTN protein from different species revealed a high degree of conservation in its functional domains and 3D structure. Furthermore, our in vitro transient-reporter analysis identified a functional promoter in the upstream region of the zebrafish optn gene, along with a region important for its transcription regulation. Site-directed mutagenesis revealed that the NF-κB motif is responsible for the activation of this region. In conclusion, with this study, we characterize zebrafish optn and our results indicate that zebrafish can be considered as an alternative model to study OPTN's biological role in bone-related diseases.

Granulin as an important immune molecule involved in lamprey tissue repair and regeneration by promoting cell proliferation and migration

Progranulin (PGRN) is an autocrine growth factor that regulates cell proliferation, migration, wound healing, and tissue repair in mammals. Lamprey is the most primitive of the extant vertebrates and is regarded as the survivor of a once flourishing group of paleozoic vertebrates, with a history of more than 500 million years. To date, the evolutionary dynamics and the underlying function of the PGRNs remain largely unclear in lamprey. Here, we screened four genes encoding PGRNs from the genomes of Lethenteron reissneri and Petromyzon marinus, including one long form (named Lr-PGRN-L) and three short forms (named Lr-PGRN-S1, Lr-PGRN-S2, and Lr-PGRN-S3), and performed phylogenetic tree, functional domain, and synteny analyses to identify the evolutionary history of the four Lr-PGRNs. In addition, the expressions of the four Lr-pgrn family genes and the immune response against various pathogenic challenges were also investigated. We found that these genes were widely distributed in various tissues of lamprey and performed a variety of functions. Moreover, our results suggest that Lr-PGRN-S1 induces cell migration and proliferation, and is involved in repair after skin and spinal cord injury under appropriate conditions. Our findings are valuable because they improve the understanding of the evolutionary relationship of vertebrate pgrn genes, as well as providing new insights into the diverse and important roles of Lr-PGRNs.

Temporal Dynamic Analysis of Alternative Splicing During Embryonic Development in Zebrafish

Alternative splicing is pervasive in mammalian genomes and involved in embryo development, whereas research on crosstalk of alternative splicing and embryo development was largely restricted to mouse and human and the alternative splicing regulation during embryogenesis in zebrafish remained unclear. We constructed the alternative splicing atlas at 18 time-course stages covering maternal-to-zygotic transition, gastrulation, somitogenesis, pharyngula stages, and post-fertilization in zebrafish. The differential alternative splicing events between different developmental stages were detected. The results indicated that abundance alternative splicing and differential alternative splicing events are dynamically changed and remarkably abundant during the maternal-to-zygotic transition process. Based on gene expression profiles, we found splicing factors are expressed with specificity of developmental stage and largely expressed during the maternal-to-zygotic transition process. The better performance of cluster analysis was achieved based on the inclusion level of alternative splicing. The biological function analysis uncovered the important roles of alternative splicing during embryogenesis. The identification of isoform switches of alternative splicing provided a new insight into mining the regulated mechanism of transcript isoforms, which always is hidden by gene expression. In conclusion, we inferred that alternative splicing activation is synchronized with zygotic genome activation and discovered that alternative splicing is coupled with transcription during embryo development in zebrafish. We also unveiled that the temporal expression dynamics of splicing factors during embryo development, especially co-orthologous splicing factors. Furthermore, we proposed that the inclusion level of alternative splicing events can be employed for cluster analysis as a novel parameter. This work will provide a deeper insight into the regulation of alternative splicing during embryogenesis in zebrafish.

Comparing the folding landscapes of evolutionarily divergent procaspase-3

All caspases evolved from a common ancestor and subsequently developed into two general classes, inflammatory or apoptotic caspases. The caspase-hemoglobinase fold has been conserved throughout nearly one billion years of evolution and is utilized for both the monomeric and dimeric subfamilies of apoptotic caspases, called initiator and effector caspases, respectively. We compared the folding and assembly of procaspase-3b from zebrafish to that of human effector procaspases in order to examine the conservation of the folding landscape. Urea-induced equilibrium folding/unfolding of procaspase-3b showed a minimum three-state folding pathway, where the native dimer isomerizes to a partially folded dimeric intermediate, which then unfolds. A partially folded monomeric intermediate observed in the folding landscape of human procaspase-3 is not well-populated in zebrafish procaspase-3b. By comparing effector caspases from different species, we show that the effector procaspase dimer undergoes a pH-dependent conformational change, and that the conformational species in the folding landscape exhibit similar free energies. Together, the data show that the landscape for the caspase-hemoglobinase fold is conserved, yet it provides flexibility for species-specific stabilization or destabilization of folding intermediates resulting in changes in stability. The common pH-dependent conformational change in the native dimer, which yields an enzymatically inactive species, may provide an additional, albeit reversible, mechanism for controlling caspase activity in the cell.

The Zebrafish model in dermatology: an update for clinicians

Recently, the zebrafish has been established as one of the most important model organisms for medical research. Several studies have proved that there is a high level of similarity between human and zebrafish genomes, which encourages the use of zebrafish as a model for understanding human genetic disorders, including cancer. Interestingly, zebrafish skin shows several similarities to human skin, suggesting that this model organism is particularly suitable for the study of neoplastic and inflammatory skin disorders. This paper appraises the specific characteristics of zebrafish skin and describes the major applications of the zebrafish model in dermatological research.

Zebrafish mutants in vegfab can affect endothelial cell proliferation without altering ERK phosphorylation and are phenocopied by loss of PI3K signaling

The formation of appropriately patterned blood vessel networks requires endothelial cell migration and proliferation. Signaling through the Vascular Endothelial Growth Factor A (VEGFA) pathway is instrumental in coordinating these processes. mRNA splicing generates short (diffusible) and long (extracellular matrix bound) Vegfa isoforms. The differences between these isoforms in controlling cellular functions are not understood. In zebrafish, vegfaa generates short and long isoforms, while vegfab only generates long isoforms. We found that mutations in vegfaa had an impact on endothelial cell (EC) migration and proliferation. Surprisingly, mutations in vegfab more strongly affected EC proliferation in distinct blood vessels, such as intersegmental blood vessels in the zebrafish trunk and central arteries in the head. Analysis of downstream signaling pathways revealed no change in MAPK (ERK) activation, while inhibiting PI3 kinase signaling phenocopied vegfab mutant phenotypes in affected blood vessels. Together, these results suggest that extracellular matrix bound Vegfa might act through PI3K signaling to control EC proliferation in a distinct set of blood vessels during angiogenesis.

Evolutionary Genomics Reveals Multiple Functions of Arylalkylamine N-Acetyltransferase in Fish

As an important hormone, melatonin participates in endocrine regulation of diverse functions in vertebrates. Its biosynthesis is catalyzed by four cascaded enzymes, among them, arylalkylamine N-acetyltransferase (AANAT) is the most critical one. Although only single aanat gene has been identified in most groups of vertebrates, researchers including us have determined that fish have the most diverse of aanat genes (aanat1a, aanat1b, and aanat2), playing various potential roles such as seasonal migration, amphibious aerial vision, and cave or deep-sea adaptation. With the rapid development of genome and transcriptome sequencing, more and more putative sequences of fish aanat genes are going to be available. Related phylogeny and functional investigations will enrich our understanding of AANAT functions in various fish species.

Spironolactone affects cardiovascular and craniofacial development in zebrafish embryos (Danio rerio)

Spironolactone, a potassium-sparing diuretic and aldosterone antagonist, is a mineralocorticoid hormone commonly prescribed to patients suffering from heart failure, hirsutism, dermatological afflictions, and hypertension. Interestingly, relatively little work has been done on the development of vertebrate embryos after exposure to this compound. Here, we treat zebrafish embryos with spironolactone at 10^-6 M, 10^-7 M, or 10^-8 M, and observe them after three to seven days of exposure. While no effect was observed in mortality, we did detect differences in cardiovascular development at 3 dpf and craniofacial development at 5 dpf. At 10^-6 M, smaller atria, ventricles, and blood vessels were observed. The highest concentrations also caused a longer ceratohyal/Meckel's distance, longer palatoquadrate, and smaller angles between the palatoquadrate and both the ceratohyal and Meckel's. Further research of spironolactone's effects on embryonic development could lead to a better understanding of the compound resulting in improved public and environmental health.

M-CSFR/CSF1R signaling regulates myeloid fates in zebrafish via distinct action of its receptors and ligands

csf1ra and csf1rb are indispensable for macrophage differentiation and, together with csf1a, regulate embryonic macrophage fates.

Il34 regulates zebrafish granulocyte development through csf1rb.

Macrophage colony-stimulating factor receptor (M-CSFR/CSF1R) signaling is crucial for the differentiation, proliferation, and survival of myeloid cells. The CSF1R pathway is a promising therapeutic target in many human diseases, including neurological disorders and cancer. Zebrafish are commonly used for human disease modeling and preclinical therapeutic screening. Therefore, it is necessary to understand the proper function of cytokine signaling in zebrafish to reliably model human-related diseases. Here, we investigate the roles of zebrafish Csf1rs and their ligands (Csf1a, Csf1b, and Il34) in embryonic and adult myelopoiesis. The proliferative effect of exogenous Csf1a on embryonic macrophages is connected to both receptors, Csf1ra and Csf1rb, however there is no evident effect of Csf1b in zebrafish embryonic myelopoiesis. Furthermore, we uncover an unknown role of Csf1rb in zebrafish granulopoiesis. Deregulation of Csf1rb signaling leads to failure in myeloid differentiation, resulting in neutropenia throughout the whole lifespan. Surprisingly, Il34 signaling through Csf1rb seems to be of high importance as both csf1rb^Δ4bp-deficient and il34^Δ5bp-deficient zebrafish larvae lack granulocytes. Our single-cell RNA sequencing analysis of adult whole kidney marrow (WKM) hematopoietic cells suggests that csf1rb is expressed mainly by blood and myeloid progenitors, and the expression of csf1ra and csf1rb is nonoverlapping. We point out differentially expressed genes important in hematopoietic cell differentiation and immune response in selected WKM populations. Our findings could improve the understanding of myeloid cell function and lead to the further study of CSF1R pathway deregulation in disease, mostly in cancerogenesis.

Long-read sequencing of the zebrafish genome reorganizes genomic architecture

BACKGROUND

Nanopore sequencing technology has revolutionized the field of genome biology with its ability to generate extra-long reads that can resolve regions of the genome that were previously inaccessible to short-read sequencing platforms. Over 50% of the zebrafish genome consists of difficult to map, highly repetitive, low complexity elements that pose inherent problems for short-read sequencers and assemblers.

RESULTS

We used long-read nanopore sequencing to generate a de novo assembly of the zebrafish genome and compared our assembly to the current reference genome, GRCz11. The new assembly identified 1697 novel insertions and deletions over one kilobase in length and placed 106 previously unlocalized scaffolds. We also discovered additional sites of retrotransposon integration previously unreported in GRCz11 and observed the expression of these transposable elements in adult zebrafish under physiologic conditions, implying they have active mobility in the zebrafish genome and contribute to the ever-changing genomic landscape.

CONCLUSIONS

We used nanopore sequencing to improve upon and resolve the issues plaguing the current zebrafish reference assembly, GRCz11. Zebrafish is a prominent model of human disease, and our corrected assembly will be useful for studies relying on interspecies comparisons and precise linkage of genetic events to disease phenotypes.

Genome-wide characterization of gap junction (connexins and pannexins) genes in turbot (Scophthalmus maximus L.): evolution and immune response following Vibrio anguillarum infection

Gap junction (GJ), a special intercellular junction between different cell types, directly connects the cytoplasm of adjacent cells, allows various molecules, ions and electrical impulses to pass through the intercellular regulatory gate, and plays vital roles in response to bacterial infection. Up to date, the information about the GJ in turbot (Scophthalmus maximus L.) is still limited. In current study, 43 gap junction genes were identified in turbot, phylogeny analysis suggested that gap junctions from turbot and other species were clustered into six groups, GJA, GJB, GJC, GJD, GJE and PANX, and turbot GJs together with respective GJs from Japanese flounder, half-smooth tongue sole and large yellow croaker, sharing same ancestors. In addition, these 43 GJ genes distributed in different chromosomes unevenly. According to gene structure and domain analysis, these genes (in GJA-GJE group) were highly conserved in that most of them contain the transmembrane area, connexin domain (CNX) and cysteine-rich domain (connexin CCC), while PANXs contain Pfam Innexin. Although only one tandem duplication was identified in turbot gap junction gene, 235 pairs of segmental duplications were identified in the turbot genome. To further investigate their evolutionary relationships, Ka/Ks was calculated, and results showed that most ratios were lower than 1, indicating they had undergone negative selection. Finally, expression analysis showed that gap junction genes were widely distributed in turbot tissues and significantly regulated after Vibrio anguillarum infection. Taken together, our research could provide valuable information for further exploration of the function of gap junction genes in teleost.