De novo assembly of the goldfish (Carassius auratus) genome and the evolution of genes after whole-genome duplication

Olfactory toxicity of tetrabromobisphenol A to the goldfish Carassius auratus

Tetrabromobisphenol A (TBBPA) is one of the most extensively used brominated flame retardants and its increasing use in consumer products has raised concerns about its ecotoxicity. Given the ubiquity of TBBPA in aquatic environments, it is inevitable that these chemicals will enter the olfactory chambers of fish via water currents. Nevertheless, the olfactory toxicity of TBBPA to aquatic organisms and the underlying toxic mechanisms have yet to be elucidated. Therefore, we investigated the olfactory toxicity of TBBPA in the goldfish Carassius auratus, a model organism widely used in sensory biology. Results showed that exposure to TBBPA resulted in abnormal olfactory-mediated behaviors and diminished electro-olfactogram (EOG) responses, indicating reduced olfactory acuity. To uncover the underlying mechanisms of action, we examined the structural integrity of the olfactory epithelium (OE), expression levels of olfactory G protein-coupled receptors (GPCRs), enzymatic activities of ion transporters, and fluctuations in neurotransmitters. Additionally, comparative transcriptomic analysis was employed to investigate the molecular mechanisms further. Our study demonstrates for the first time that TBBPA at environmentally relevant levels can adversely affect the olfactory sensitivity of aquatic organisms by interfering with the transmission of aqueous stimuli to olfactory receptors, impeding the binding of odorants to their receptors, disrupting the olfactory signal transduction pathway, and ultimately affecting the generation of action potentials.

Recent and Recurrent Autopolyploidization Fueled Diversification of Snow Carp on the Tibetan Plateau

Whole-genome duplication (WGD), or polyploidization, is a major contributor to biodiversity. However, the establishment and survival of WGDs are often considered to be stochastic, since elucidating the processes of WGD establishment remains challenging. In the current study, we explored the processes leading to polyploidy establishment in snow carp (Cyprinidae: Schizothoracinae), a predominant component of the ichthyofauna of the Tibetan Plateau and its surrounding areas. Using large-scale genomic data from isoform sequencing, we analyzed ohnolog genealogies and divergence in hundreds to thousands of gene families across major snow carp lineages. Our findings demonstrated that independent autopolyploidization subsequent to speciation was prevalent, while autopolyploidization followed by speciation also occurred in the diversification of snow carp. This was further supported by matrilineal divergence and drainage evolution evidence. Contrary to the long-standing hypothesis that ancient polyploidization preceded the diversification of snow carp, we determined that polyploidy in extant snow carp was established by recurrent autopolyploidization events during the Pleistocene. These findings indicate that the diversification of extant snow carp resembles a coordinated duet: first, the uplift of the Tibetan Plateau orchestrated the biogeography and diversification of their diploid progenitors; then, the extensive Pliocene-Pleistocene climate changes acted as relay runners, further fueling diversification through recurrent autopolyploidization. Overall, this study not only reveals a hitherto unrecognized recent WGD lineage in vertebrates but also advances current understanding of WGD processes, emphasizing that WGD establishment is a nonstochastic event, emerging from numerous adaptations to environmental challenges and recurring throughout evolutionary history rather than merely in plants.

Is there direct photoentrainment in the goldfish liver? Wavelength-dependent regulation of clock genes and investigation of the opsin 7 family

Widespread direct photoentrainment in zebrafish peripheral tissues is linked to diverse non-visual opsins. To explore whether this broadly distributed photosensitivity is specific to zebrafish or is a general teleost feature, we investigated hepatic photosynchronization in goldfish. First, we focused on the opsin 7 family (OPN7, a key peripheral novel opsin in zebrafish), investigating its presence in the goldfish liver. Subsequently, we studied whether light can directly entrain the goldfish liver and retina clocks. Silico analysis revealed seven OPN7 paralogs from four gene families, suggesting expansion through whole-genome and tandem duplications. The paralogs of families OPN7a, OPN7b, and OPN7d were mainly localized in neural tissues, while OPN7c paralogs were more abundant in peripheral tissues-including the liver-suggesting divergent roles. Light (independently of the wavelength employed) directly induced the per2a clock gene in the retina both in vivo and in vitro, confirming expected photoentrainment. However, in the liver, photoinduction of per1a and cry1a only occurred in vivo, not in vitro. These results suggest an indirect light-entrainment mechanism of the goldfish hepatic clock, possibly mediated by other oscillators or photosensitive organs. Our findings challenge the assumption of widespread direct photosensitivity in the peripheral tissues of teleosts. Further research is needed to understand the role of tissue-specific photoentrainment and non-visual opsins in diverse teleost species.

Stem and progenitor cell proliferation are independently regulated by cell type-specific cyclinD genes

Regeneration and homeostatic turnover of solid tissues depend on the proliferation of symmetrically dividing adult stem cells, which either remain stem cells or differentiate based on their niche position. Here we demonstrate that in zebrafish lateral line sensory organs, stem and progenitor cell proliferation are independently regulated by two cyclinD genes. Loss of ccnd2a impairs stem cell proliferation during development, while loss of ccndx disrupts hair cell progenitor proliferation but allows normal differentiation. Notably, ccnd2a can functionally replace ccndx, indicating that the respective effects of these Cyclins on proliferation are due to cell type-specific expression. However, even though hair cell progenitors differentiate normally in ccndx mutants, they are mispolarized due to hes2 and Emx2 downregulation. Thus, regulated proliferation ensures that equal numbers of hair cells are polarized in opposite directions. Our study reveals cell type-specific roles for cyclinD genes in regulating the different populations of symmetrically dividing cells governing organ development and regeneration, with implications for regenerative medicine and disease.

Cardiac regeneration in goldfish (Carassius auratus) associated with increased expression of key extracellular matrix molecules

Cardiac regeneration is a natural phenomenon that occurs in many species outside of humans. The goldfish (Carassius auratus) is an understudied model of cardiac wound response, despite its ubiquity as pets as well as its relationship to the better-studied zebrafish. In this study, we examined the response of the goldfish heart to a resection injury. We found that by 70 days post-injury, goldfish scarlessly heal cardiac wounds under a certain size, with local cardiomyocyte proliferation driving the restoration of the myocardial layer. We also found the upregulation of extracellular matrix components related to cardiac regeneration in the injury site. This upregulation correlated with the level of cardiomyocyte proliferation occurring in the injury site, indicating an association between the two that warrants further exploration.

Pax7 is involved in leucophore formation in goldfish and gene knockout improves the transparency of transparent goldfish

Lines with few or no pigment cells have been established in fishes, and these lines are useful for bioimaging. The transparent goldfish (tra) line previously established by N-ethyl-N-nitrosourea (ENU) mutagenesis is also suitable for such experiments. However, in the case of tra, leucophores form in the adult fish, making it difficult to observe the organs inside body from outside the body. In this study, we attempted to create a knockout line of the pax7a and pax7b genes, which are thought to be involved in the formation of leucophores, to further improve the transparency of tra strain.Mutations were introduced by microinjection of the CRISPR/Cas9 mixture into single-cell embryos, mutant individuals were found in F0, and the next generation was generated to confirm the mutation patterns. As a result, multiple mutation patterns, including knockout, were obtained. The same pattern of knockout F1 with pax7a and pax7b mutations was crossed to generate a homozygous knockout in F2.In the resulting pax7b-/- (tra) fish but not in pax7a-/- (tra) fish, the number of leucophores was reduced compared to that in tra, and the transparency of the body was improved. It was suggested that pax7b plays an important role in leucophore formation in goldfish. The established transparent pax7b-/- (tra) goldfish line will be a useful model for bioimaging of the body interior.

Whole-genome duplication in an algal symbiont bolsters coral heat tolerance

The algal endosymbiont Durusdinium trenchii enhances the resilience of coral reefs under thermal stress. D. trenchii can live freely or in endosymbiosis, and the analysis of genetic markers suggests that this species has undergone whole-genome duplication (WGD). However, the evolutionary mechanisms that underpin the thermotolerance of this species are largely unknown. Here, we present genome assemblies for two D. trenchii isolates, confirm WGD in these taxa, and examine how selection has shaped the duplicated genome regions using gene expression data. We assess how the free-living versus endosymbiotic lifestyles have contributed to the retention and divergence of duplicated genes, and how these processes have enhanced the thermotolerance of D. trenchii. Our combined results suggest that lifestyle is the driver of post-WGD evolution in D. trenchii, with the free-living phase being the most important, followed by endosymbiosis. Adaptations to both lifestyles likely enabled D. trenchii to provide enhanced thermal stress protection to the host coral.

Mapping the chromatin accessibility landscape of zebrafish embryogenesis at single-cell resolution by SPATAC-seq

Currently, the dynamic accessible elements that determine regulatory programs responsible for the unique identity and function of each cell type during zebrafish embryogenesis lack detailed study. Here we present SPATAC-seq: a split-pool ligation-based assay for transposase-accessible chromatin using sequencing. Using SPATAC-seq, we profiled chromatin accessibility in more than 800,000 individual nuclei across 20 developmental stages spanning the sphere stage to the early larval protruding mouth stage. Using this chromatin accessibility map, we identified 604 cell states and inferred their developmental relationships. We also identified 959,040 candidate cis-regulatory elements (cCREs) and delineated development-specific cCREs, as well as transcription factors defining diverse cell identities. Importantly, enhancer reporter assays confirmed that the majority of tested cCREs exhibited robust enhanced green fluorescent protein expression in restricted cell types or tissues. Finally, we explored gene regulatory programs that drive pigment and notochord cell differentiation. Our work provides a valuable open resource for exploring driver regulators of cell fate decisions in zebrafish embryogenesis.

The dynamic transcriptomic response of the goldfish brain under chronic hypoxia

Oxygen is essential to fuel aerobic metabolism. Some species evolved mechanisms to tolerate periods of severe hypoxia and even anoxia in their environment. Among them, goldfish (Carassius auratus) are unique, in that they do not enter a comatose state under severely hypoxic conditions. There is thus significant interest in the field of comparative physiology to uncover the mechanistic basis underlying hypoxia tolerance in goldfish, with a particular focus on the brain. Taking advantage of the recently published and annotated goldfish genome, we profile the transcriptomic response of the goldfish brain under normoxic (21 kPa oxygen saturation) and, following gradual reduction, constant hypoxic conditions after 1 and 4 weeks (2.1 kPa oxygen saturation). In addition to analyzing differentially expressed protein-coding genes and enriched pathways, we also profile differentially expressed microRNAs (miRs). Using in silico approaches, we identify possible miR-mRNA relationships. Differentially expressed transcripts compared to normoxia were either common to both timepoints of hypoxia exposure (n = 174 mRNAs; n = 6 miRs), or exclusive to 1-week (n = 441 mRNAs; n = 23 miRs) or 4-week hypoxia exposure (n = 491 mRNAs; n = 34 miRs). Under chronic hypoxia, an increasing number of transcripts, including those of paralogous genes, was downregulated over time, suggesting a decrease in transcription. GO-terms related to the vascular system, oxidative stress, stress signalling, oxidoreductase activity, nucleotide- and intermediary metabolism, and mRNA posttranscriptional regulation were found to be enriched under chronic hypoxia. Known 'hypoxamiRs', such as miR-210-3p/5p, and miRs such as miR-29b-3p likely contribute to posttranscriptional regulation of these pathways under chronic hypoxia in the goldfish brain.

Genomic signatures associated with recurrent scale loss in cyprinid fish

Scale morphology represents a fundamental feature of fish and a key evolutionary trait underlying fish diversification. Despite frequent and recurrent scale loss throughout fish diversification, comprehensive genome-wide analyses of the genomic signatures associated with scale loss in divergent fish lineages remain scarce. In the current study, we investigated genome-wide signatures, specifically convergent protein-coding gene loss, amino acid substitutions, and cis-regulatory sequence changes, associated with recurrent scale loss in two divergent Cypriniformes lineages based on large-scale genomic, transcriptomic, and epigenetic data. Results demonstrated convergent changes in many genes related to scale formation in divergent scaleless fish lineages, including loss of P/Q-rich scpp genes (e.g. scpp6 and scpp7), accelerated evolution of non-coding elements adjacent to the fgf and fgfr genes, and convergent amino acid changes in genes (e.g. snap29) under relaxed selection. Collectively, these findings highlight the existence of a shared genetic architecture underlying recurrent scale loss in divergent fish lineages, suggesting that evolutionary outcomes may be genetically repeatable and predictable in the convergence of scale loss in fish.

Asymmetric evolution of ISG15 homologs and the immune adaptation to LBUSV infection in spotted seabass (Lateolabrax maculatus)

The battle between host and viral is ubiquitous across all ecosystems. Despite this, research is scarce on the antiviral characteristics of fish, particularly in those that primarily rely on innate immune responses. This study, comprehensively explored the genetic and antiviral features of ISG15 in spotted seabass, focusing on its response to largemouth bass ulcerative syndrome virus (LBUSV). Through whole-genome BLAST and PCR cloning, two ISG15 homologs, namely LmISG15a and LmISG15b, were identified in spotted seabass, both encoding highly conserved proteins. However, a distinctive contrast emerged in their expression patterns, with LmISG15a exhibiting high expression in immune organs while LmISG15b remained largely silent across various organs. Regulatory elements analysis indicated an asymmetric evolution of the two ISG15s, with the minimal expression of LmISG15b may attribute to the loss of a necessary ISRE and an additional instability "ATTTA" motif. Association analysis demonstrated a significant correlation between LmISG15a expression and LBUSV infection. Subsequent antiviral activity detection revealed that LmISG15a interacted with LBUSV, inhibiting its replication by activating ISGylation and downstream pro-inflammatory mediators. In summary, this study unveils a distinct evolutionary strategy of fish antiviral gene ISG15 and delineates its kinetic characteristics in response to LBUSV infection.

Evolutionary divergence of subgenomes in common carp provides insights into speciation and allopolyploid success

Hybridization and polyploidization have made great contributions to speciation, heterosis, and agricultural production within plants, but there is still limited understanding and utilization in animals. Subgenome structure and expression reorganization and cooperation post hybridization and polyploidization are essential for speciation and allopolyploid success. However, the mechanisms have not yet been comprehensively assessed in animals. Here, we produced a high-fidelity reference genome sequence for common carp, a typical allotetraploid fish species cultured worldwide. This genome enabled in-depth analysis of the evolution of subgenome architecture and expression responses. Most genes were expressed with subgenome biases, with a trend of transition from the expression of subgenome A during the early stages to that of subgenome B during the late stages of embryonic development. While subgenome A evolved more rapidly, subgenome B contributed to a greater level of expression during development and under stressful conditions. Stable dominant patterns for homoeologous gene pairs both during development and under thermal stress suggest a potential fixed heterosis in the allotetraploid genome. Preferentially expressing either copy of a homoeologous gene at higher levels to confer development and response to stress indicates the dominant effect of heterosis. The plasticity of subgenomes and their shifting of dominant expression during early development, and in response to stressful conditions, provide novel insights into the molecular basis of the successful speciation, evolution, and heterosis of the allotetraploid common carp.

Genome composition and GC content influence loci distribution in reduced representation genomic studies

BACKGROUND

Genomic architecture is a key evolutionary trait for living organisms. Due to multiple complex adaptive and neutral forces which impose evolutionary pressures on genomes, there is a huge variability of genomic features. However, their variability and the extent to which genomic content determines the distribution of recovered loci in reduced representation sequencing studies is largely unexplored.

RESULTS

Here, by using 80 genome assemblies, we observed that whereas plants primarily increase their genome size by expanding their intergenic regions, animals expand both intergenic and intronic regions, although the expansion patterns differ between deuterostomes and protostomes. Loci mapping in introns, exons, and intergenic categories obtained by in silico digestion using 2b-enzymes are positively correlated with the percentage of these regions in the corresponding genomes, suggesting that loci distribution mostly mirrors genomic architecture of the selected taxon. However, exonic regions showed a significant enrichment of loci in all groups regardless of the used enzyme. Moreover, when using selective adaptors to obtain a secondarily reduced loci dataset, the percentage and distribution of retained loci also varied. Adaptors with G/C terminals recovered a lower percentage of selected loci, with a further enrichment of exonic regions, while adaptors with A/T terminals retained a higher percentage of loci and slightly selected more intronic regions than expected.

CONCLUSIONS

Our results highlight how genome composition, genome GC content, RAD enzyme choice and use of base-selective adaptors influence reduced genome representation techniques. This is important to acknowledge in population and conservation genomic studies, as it determines the abundance and distribution of loci.

Identification and effective regulation of scarb1 gene involved in pigmentation change in autotetraploid Carassius auratus

The autotetraploid Carassius auratus (4nRR, 4 n=200, RRRR) is derived from whole-genome duplication of Carassius auratus red var. (RCC, 2 n=100, RR). In the current study, we demonstrated that chromatophores and pigment changes directly caused the coloration and variation of 4nRR skin (red in RCC, brownish-yellow in 4nRR). To further explore the molecular mechanisms underlying coloration formation and variation in 4nRR, we performed transcriptome profiling and molecular functional verification in RCC and 4nRR. Results revealed that scarb1, associated with carotenoid metabolism, underwent significant down-regulation in 4nRR. Efficient editing of this candidate pigment gene provided clear evidence of its significant role in RCC coloration. Subsequently, we identified four divergent scarb1 homeologs in 4nRR: two original scarb1 homeologs from RCC and two duplicated ones. Notably, three of these homeologs possessed two highly conserved alleles, exhibiting biased and allele-specific expression in the skin. Remarkably, after precise editing of both the original and duplicated scarb1 homeologs and/or alleles, 4nRR individuals, whether singly or multiply mutated, displayed a transition from brownish-yellow skin to a cyan-gray phenotype. Concurrently, the proportional areas of the cyan-gray regions displayed a gene-dose correlation. These findings illustrate the subfunctionalization of duplicated scarb1, with all scarb1 genes synergistically and equally contributing to the pigmentation of 4nRR. This is the first report concerning the functional differentiation of duplicated homeologs in an autopolyploid fish, substantially enriching our understanding of coloration formation and change within this group of organisms.

The hagfish genome and the evolution of vertebrates

As the only surviving lineages of jawless fishes, hagfishes and lampreys provide a crucial window into early vertebrate evolution1-3. Here we investigate the complex history, timing and functional role of genome-wide duplications4-7 and programmed DNA elimination8,9 in vertebrates in the light of a chromosome-scale genome sequence for the brown hagfish Eptatretus atami. Combining evidence from syntenic and phylogenetic analyses, we establish a comprehensive picture of vertebrate genome evolution, including an auto-tetraploidization (1RV) that predates the early Cambrian cyclostome-gnathostome split, followed by a mid-late Cambrian allo-tetraploidization (2RJV) in gnathostomes and a prolonged Cambrian-Ordovician hexaploidization (2RCY) in cyclostomes. Subsequently, hagfishes underwent extensive genomic changes, with chromosomal fusions accompanied by the loss of genes that are essential for organ systems (for example, genes involved in the development of eyes and in the proliferation of osteoclasts); these changes account, in part, for the simplification of the hagfish body plan1,2. Finally, we characterize programmed DNA elimination in hagfish, identifying protein-coding genes and repetitive elements that are deleted from somatic cell lineages during early development. The elimination of these germline-specific genes provides a mechanism for resolving genetic conflict between soma and germline by repressing germline and pluripotency functions, paralleling findings in lampreys10,11. Reconstruction of the early genomic history of vertebrates provides a framework for further investigations of the evolution of cyclostomes and jawed vertebrates.

Evolutionary mechanisms and practical significance of reproductive success and clonal diversity in unisexual vertebrate polyploids

Unisexual reproduction is generally relevant to polyploidy, and unisexual vertebrates are often considered an evolutionary "dead end" due to the accumulation of deleterious mutations and absence of genetic diversity. However, some unisexual polyploids have developed strategies to avoid genomic decay, and thus provide ideal models to unveil unexplored evolutionary mechanisms, from the reproductive success to clonal diversity creation. This article reviews the evolutionary mechanisms for overcoming meiotic barrier and generating genetic diversity in unisexual vertebrates, and summarizes recent research advancements in the polyploid Carassius complex. Gynogenetic gibel carp (Carassius gibelio) is a unique amphitriploid that has undergone a recurrent autotriploidy and has overcome the bottleneck of triploid sterility via gynogenesis. Recently, an efficient strategy in which ploidy changes, including from amphitriploid to amphitetraploid, then from amphitetraploid to novel amphitriploid, drive unisexual-sexual-unisexual reproduction transition and clonal diversity has been revealed. Based on this new discovery, multigenomic reconstruction biotechnology has been used to breed a novel strain with superior growth and stronger disease resistance. Moreover, a unique reproduction mode that combines both abilities of ameiotic oogenesis and sperm-egg fusion, termed as ameio-fusiongensis, has been discovered, and it provides an efficient approach to synthesize sterile allopolyploids. In order to avoid ecological risks upon escape and protect the sustainable property rights of the aquaculture seed industry, a controllable fertility biotechnology approach for precise breeding is being developed by integrating sterile allopolyploid synthesis and gene-editing techniques. This review provides novel insights into the origin and evolution of unisexual vertebrates and into the attempts being made to exploit new breeding biotechnologies in aquaculture.

A phylogenetics-based nomenclature system for steroid receptors in teleost fishes

Teleost fishes have emerged as tractable models for studying the neuroendocrine regulation of social behavior via molecular genetic techniques, such as CRISPR/Cas9 gene editing. Moreover, teleosts provide an opportunity to investigate the evolution of steroid receptors and their functions, as species within this lineage possess novel steroid receptor paralogs that resulted from a teleost-specific whole genome duplication. Although teleost fishes have grown in popularity as models for behavioral neuroendocrinology, there is not a consistent nomenclature system for steroid receptors and their genes, which may impede a clear understanding of steroid receptor paralogs and their functions. Here, we used a phylogenetic approach to assess the relatedness of protein sequences encoding steroid receptor paralogs in 18 species from 12 different orders of the Infraclass Teleostei. While most similarly named sequences grouped based on the established phylogeny of the teleost lineage, our analysis revealed several inconsistencies in the nomenclature of steroid receptor paralogs, particularly for sequences encoding estrogen receptor beta (ERβ). Based on our results, we propose a nomenclature system for teleosts in which Greek symbols refer to proteins and numbers refer to genes encoding different subtypes of steroid receptors within the five major groups of this nuclear receptor subfamily. Collectively, our results bridge a critical gap by providing a cohesive naming system for steroid receptors in teleost fishes, which will serve to improve communication, promote collaboration, and enhance our understanding of the evolution and function of steroid receptors across vertebrates.

The Goldfish Genome and Its Utility for Understanding Gene Regulation and Vertebrate Body Morphology

Goldfish, widely viewed as an ornamental fish, is a member of Cyprinidae family and has a very long history in research for both genetics and physiology studies. Among Cyprinidae, the chromosomal locations of orthologs and the amino acid sequences are usually highly conserved. Adult goldfish are 1000 times larger than adult zebrafish (who are in the same family of fishes), which can make it easier to perform several types of experiments compared to their zebrafish cousins. Comparing mutant phenotypes in orthologous genes between goldfish and zebrafish can often be very informative and provide a deeper insight into the gene function than studying the gene in either species alone. Comparative genomics and phenotypic comparisons between goldfish and zebrafish will provide new opportunities for understanding the development and evolution of body forms in the vertebrate lineage.

Gene Characterization of Nocturnin Paralogues in Goldfish: Full Coding Sequences, Structure, Phylogeny and Tissue Expression

The aim of this work is the full characterization of all the nocturnin (noc) paralogues expressed in a teleost, the goldfish. An in silico analysis of the evolutive origin of noc in Osteichthyes is performed, including the splicing variants and new paralogues appearing after teleostean 3R genomic duplication and the cyprinine 4Rc. After sequencing the full-length mRNA of goldfish, we obtained two isoforms for noc-a (noc-aa and noc-ab) with two splice variants (I and II), and only one for noc-b (noc-bb) with two transcripts (II and III). Using the splicing variant II, the prediction of the secondary and tertiary structures renders a well-conserved 3D distribution of four α-helices and nine β-sheets in the three noc isoforms. A synteny analysis based on the localization of noc genes in the patrilineal or matrilineal subgenomes and a phylogenetic tree of protein sequences were accomplished to stablish a classification and a long-lasting nomenclature of noc in goldfish, and valid to be extrapolated to allotetraploid Cyprininae. Finally, both goldfish and zebrafish showed a broad tissue expression of all the noc paralogues. Moreover, the enriched expression of specific paralogues in some tissues argues in favour of neo- or subfunctionalization.

Re-identification and characterization of grass carp Ctenopharyngodon idella TLR20

Toll-like receptors (TLRs) play a crucial role in the recognition of microbial-associated molecular patterns in the innate immune system. Fish TLRs have undergone significant gene expansion to adapt to complex aquatic environments. Among them, TLR20 from the TLR11 family actively responds to viral and bacterial invasions. Previous studies have reported two TLR20s in grass carp (Ctenopharyngodon idella), and in this study, we revised this conclusion. Based on the latest grass carp genome, we identified a new TLR20 member. These three TLR20s are arranged in tandem on chromosome 9, indicating that they are generated by gene duplication events. They were renamed CiTLR20.1 to CiTLR20.3 based on their chromosomal positions. The CiTLR20s in C. idella exhibit higher similarities with those in Danio rerio, Cyprinus carpio, and Megalobrama amblycephala, and lower similarities with those in other distantly related fish species. Selective pressure analysis revealed low conservation and negative evolution of TLR20s during evolution. The 3D structures of the three TLR20s showed significant differences, reflecting functional variations and different downstream adaptor molecule recruitment. Transcriptome data revealed tissue distribution differences of TLR20s, with TLR20.1 showing relatively low expression levels in all the tissues, while TLR20.2 and TLR20.3 showed higher expression in the head kidney, spleen, and gill. Additionally, TLR20.2 and TLR20.3 actively responded to GCRV-II infection, with higher upregulation of TLR20.2 in response to Aeromonas hydrophila challenge. In conclusion, this study corrected the number of grass carp TLR20 members and analyzed TLR20 from an evolutionary and structural perspective, exploring its role in antiviral and antibacterial defense. This study provides reference for future research on fish TLR20.

Maternal dominance contributes to subgenome differentiation in allopolyploid fishes

Teleost fishes, which are the largest and most diverse group of living vertebrates, have a rich history of ancient and recent polyploidy. Previous studies of allotetraploid common carp and goldfish (cyprinids) reported a dominant subgenome, which is more expressed and exhibits biased gene retention. However, the underlying mechanisms contributing to observed 'subgenome dominance' remains poorly understood. Here we report high-quality genomes of twenty-one cyprinids to investigate the origin and subsequent subgenome evolution patterns following three independent allopolyploidy events. We identify the closest extant relatives of the diploid progenitor species, investigate genetic and epigenetic differences among subgenomes, and conclude that observed subgenome dominance patterns are likely due to a combination of maternal dominance and transposable element densities in each polyploid. These findings provide an important foundation to understanding subgenome dominance patterns observed in teleost fishes, and ultimately the role of polyploidy in contributing to evolutionary innovations.

Source diversity of Artemia enrichment boosts goldfish (Carassius auratus) performance, β-carotene content, pigmentation, immune-physiological and transcriptomic responses

This study aimed to assess the impact of spirulina and/or canthaxanthin-enriched Artemia on the goldfish (Carassius auratus) growth, pigmentation, blood analysis, immunity, intestine and liver histomorphology, and expression of somatolactin (SL) and growth hormone (GH) genes. Artemia was enriched with spirulina and/or canthaxanthin for 24 h. Goldfish (N = 225, 1.10 ± 0.02 g) were tested in five experimental treatments, three replicates each: (T1) fish fed a commercial diet; (T2) fish fed a commercial diet and un-enriched Artemia (UEA); (T3) fish fed a commercial diet and spirulina-enriched Artemia (SEA); (T4) fish fed a commercial diet and canthaxanthin-enriched Artemia (CEA); and (T5) fish fed a commercial diet and spirulina and canthaxanthin-enriched Artemia (SCA) for 90 days. The results showed that the use of spirulina and/or canthaxanthin increased performance, β-carotene content and polyunsaturated fatty acids of Artemia. For goldfish, T5 showed the highest growth performance, β-carotene concentration and the lowest chromatic deformity. T5 also showed improved hematology profile, serum biochemical, and immunological parameters. Histomorphology of the intestine revealed an increase in villi length and goblet cells number in the anterior and middle intestine, with normal liver structure in T5. SL and GH gene expression in the liver and brain differed significantly among treatments with a significant increase in enriched Artemia treatments compared to T1 and T2. In conclusion, the use of spirulina and/or canthaxanthin improved performance of Artemia. Feeding goldfish spirulina and/or canthaxanthin-enriched Artemia improved performance, β-carotene content, pigmentation, health status and immune-physiological response.

New reports on iron related proteins: Molecular characterization of two ferroportin genes in common carp (Cyprinus carpio L.) and its expression pattern

Iron uptake, transport, and storage require the involvement of several proteins, including ferroportin (fpn), the sole known iron efflux transporter. Due to its critical function fpn has been studied, particularly in humans. Here, we characterized the ferroportin gene in common carp (Cyprinus carpio L.) and performed RNA-seq analysis to evaluate its constitutive transcription levels across different tissues. Our results indicate that C. carpio possesses two functional fpns with distinct expression patterns, highlighting the potential for functional divergence and expression differentiation among fpns in this species.

Symmetric expression of ohnologs encoding conserved antiviral responses in tetraploid common carp suggest absence of subgenome dominance after whole genome duplication

Allopolyploids often experience subgenome dominance, with one subgenome showing higher levels of gene expression and greater gene retention. Here, we address the functionality of both subgenomes of allotetraploid common carp (Cyprinus carpio) by analysing a functional network of interferon-stimulated genes (ISGs) crucial in anti-viral immune defence. As an indicator of subgenome dominance we investigated retainment of a core set of ohnologous ISGs. To facilitate our functional genomic analysis a high quality genome was assembled (WagV4.0). Transcriptome data from an in vitro experiment mimicking a viral infection was used to infer ISG expression. Transcriptome analysis confirmed induction of 88 ISG ohnologs on both subgenomes. In both control and infected states, average expression of ISG ohnologs was comparable between the two subgenomes. Also, the highest expressing and most inducible gene copies of an ohnolog pair could be derived from either subgenome. We found no strong evidence of subgenome dominance for common carp.

Identification and functional characterization of a novel TRPA1 gene from sea cucumber Apostichopus japonicus and interaction with miR-2013 in response to salt stress

Salinity is important abiotic factor influencing sea cucumber aquaculture. This study aimed to identify and functional study of a novel transient receptor potential cation channel subfamily A member 1 (TRPA1) involved in salinity stress through interaction with miR-2013 in the sea cucumber. The full-length cDNA sequence was 1369 bp in length and encoded 138 amino acids. The TRPA1 homolog protein was a hydrophilic protein without a signal peptide and was predicted to a spatial structure of seven helices and eight random coils and two major ANK functional domains. Bioinformatic analysis and luciferase reporter assays confirmed TRPA1 as a target gene of miR-2013. Quantitative PCR revealed that miR-2013 was induced upregulation after salinity stress, while TRPA1 showed upregulated expression with maximum expression at 24 h. The expression of miR-2013 and TRPA1 was negatively regulated. Transfection experiments were conducted to validate the role of miR-2013 and TRPA1 in salinity response. The results showed that miR-2013 was upregulated and TRPA1 was downregulated after transfection with miR-2013 mimics, while miR-2013 was downregulated and TRPA1 was upregulated after transfection with miR-2013 inhibitor. Transfection with si-TRPA1 homolog resulted in upregulation of miR-2013 and downregulation of TRPA1 homolog. These findings suggest that miR-2013 can regulate the expression of TRPA1 under salt stress, and highlight the importance of miR-2013 and TRPA1 in salt stress response. miR-2013 mimics improved the survival rate, while miR-2013 inhibitor and si-TRPA1 reduced it. These findings suggest that miR-2013 and TRPA1 play important roles in sea cucumbers adaptation to salinity changes.

Hybridization led to a rewired pluripotency network in the allotetraploid Xenopus laevis

After fertilization, maternally contributed factors to the egg initiate the transition to pluripotency to give rise to embryonic stem cells, in large part by activating de novo transcription from the embryonic genome. Diverse mechanisms coordinate this transition across animals, suggesting that pervasive regulatory remodeling has shaped the earliest stages of development. Here, we show that maternal homologs of mammalian pluripotency reprogramming factors OCT4 and SOX2 divergently activate the two subgenomes of Xenopus laevis, an allotetraploid that arose from hybridization of two diploid species ~18 million years ago. Although most genes have been retained as two homeologous copies, we find that a majority of them undergo asymmetric activation in the early embryo. Chromatin accessibility profiling and CUT&RUN for modified histones and transcription factor binding reveal extensive differences in predicted enhancer architecture between the subgenomes, which likely arose through genomic disruptions as a consequence of allotetraploidy. However, comparison with diploid X. tropicalis and zebrafish shows broad conservation of embryonic gene expression levels when divergent homeolog contributions are combined, implying strong selection to maintain dosage in the core vertebrate pluripotency transcriptional program, amid genomic instability following hybridization.

Zebrafish regulatory genomic resources for disease modelling and regeneration

In the past decades, the zebrafish has become a disease model with increasing popularity owing to its advantages that include fast development, easy genetic manipulation, simplicity for imaging, and sharing conserved disease-associated genes and pathways with those of human. In parallel, studies of disease mechanisms are increasingly focusing on non-coding mutations, which require genome annotation maps of regulatory elements, such as enhancers and promoters. In line with this, genomic resources for zebrafish research are expanding, producing a variety of genomic data that help in defining regulatory elements and their conservation between zebrafish and humans. Here, we discuss recent developments in generating functional annotation maps for regulatory elements of the zebrafish genome and how this can be applied to human diseases. We highlight community-driven developments, such as DANIO-CODE, in generating a centralised and standardised catalogue of zebrafish genomics data and functional annotations; consider the advantages and limitations of current annotation maps; and offer considerations for interpreting and integrating existing maps with comparative genomics tools. We also discuss the need for developing standardised genomics protocols and bioinformatic pipelines and provide suggestions for the development of analysis and visualisation tools that will integrate various multiomic bulk sequencing data together with fast-expanding data on single-cell methods, such as single-cell assay for transposase-accessible chromatin with sequencing. Such integration tools are essential to exploit the multiomic chromatin characterisation offered by bulk genomics together with the cell-type resolution offered by emerging single-cell methods. Together, these advances will build an expansive toolkit for interrogating the mechanisms of human disease in zebrafish.

Parallel evolution of fish bi-modal breathing and expansion of olfactory receptor (OR) genes: toward a universal ORs nomenclature

Olfactory receptors (ORs) play a key role in the prime sensorial perception, being highly relevant for intra/interspecific interactions. ORs are a subgroup of G-protein coupled receptors that exhibit highly complex subgenomes in vertebrates. However, OR repertoires remain poorly studied in fish lineages, precluding finely retracing their origin, evolution, and diversification, especially in the most basal groups. Here, we conduct an exhaustive gene screening upon 43 high-quality fish genomes exhibiting varied gene repertoires (2-583 genes). While the early vertebrates performed gas exchange through gills, we hypothesize that the emergence of new breathing structures (swim bladder and paired lungs) in early osteichthyans may be associated with expansions in the ORs gene families sensitive to airborne molecules. Additionally, we verify that the OR repertoire of moderns actinopterygians has not increased as expected following a whole genome duplication, likely due to regulatory mechanisms compensating the gene load excess. Finally, we identify 25 distinct OR families, allowing us to propose an updated universal nomenclature for the fish ORs.

ncOrtho: efficient and reliable identification of miRNA orthologs

MicroRNAs (miRNAs) are post-transcriptional regulators that finetune gene expression via translational repression or degradation of their target mRNAs. Despite their functional relevance, frameworks for the scalable and accurate detection of miRNA orthologs are missing. Consequently, there is still no comprehensive picture of how miRNAs and their associated regulatory networks have evolved. Here we present ncOrtho, a synteny informed pipeline for the targeted search of miRNA orthologs in unannotated genome sequences. ncOrtho matches miRNA annotations from multi-tissue transcriptomes in precision, while scaling to the analysis of hundreds of custom-selected species. The presence-absence pattern of orthologs to 266 human miRNA families across 402 vertebrate species reveals four bursts of miRNA acquisition, of which the most recent event occurred in the last common ancestor of higher primates. miRNA families are rarely modified or lost, but notable exceptions for both events exist. miRNA co-ortholog numbers faithfully indicate lineage-specific whole genome duplications, and miRNAs are powerful markers for phylogenomic analyses. Their exceptionally low genetic diversity makes them suitable to resolve clades where the phylogenetic signal is blurred by incomplete lineage sorting of ancestral alleles. In summary, ncOrtho allows to routinely consider miRNAs in evolutionary analyses that were thus far reserved to protein-coding genes.

Transposon signatures of allopolyploid genome evolution

Hybridization brings together chromosome sets from two or more distinct progenitor species. Genome duplication associated with hybridization, or allopolyploidy, allows these chromosome sets to persist as distinct subgenomes during subsequent meioses. Here, we present a general method for identifying the subgenomes of a polyploid based on shared ancestry as revealed by the genomic distribution of repetitive elements that were active in the progenitors. This subgenome-enriched transposable element signal is intrinsic to the polyploid, allowing broader applicability than other approaches that depend on the availability of sequenced diploid relatives. We develop the statistical basis of the method, demonstrate its applicability in the well-studied cases of tobacco, cotton, and Brassica napus, and apply it to several cases: allotetraploid cyprinids, allohexaploid false flax, and allooctoploid strawberry. These analyses provide insight into the origins of these polyploids, revise the subgenome identities of strawberry, and provide perspective on subgenome dominance in higher polyploids.

An efficient approach to synthesize sterile allopolyploids through the combined reproduction mode of ameiotic oogenesis and sperm-egg fusion in the polyploid Carassius complex

The association between polyploidy and reproduction transition, which is an intriguing issue in evolutionary genetics, can also be exploited as an approach for genetic improvement in agriculture. Recently, we generated novel amphitriploids (NA3n) by integrating the genomes of the gynogenetic Carassius gibelio and sexual C. auratus, and found gynogenesis was recovered in most NA3n females (NA3n♀I). Here, we discovered a unique reproduction mode, termed ameio-fusiongenesis, which combines the abilities of both ameiotic oogenesis and sperm-egg fusion, in a few NA3n females (NA3n♀II). These females inherited ameiotic oogenesis to produce unreduced eggs from gynogenetic C. gibelio and sperm-egg fusion from sexual C. auratus. Subsequently, we utilized this unique reproduction mode to generate a group of synthetic alloheptaploids by crossing NA3n♀II with Megalobrama amblycephala. They contained all chromosomes of maternal NA3n♀II and a chromosomal set of paternal M. amblycephala. Intergenomic chromosome translocations between NA3n♀II and M. amblycephala were also observed in a few somatic cells. Primary oocytes of the alloheptaploid underwent severe apoptosis owing to incomplete double-strand break repair at prophase I. Although spermatocytes displayed similar chromosome behavior at prophase I, they underwent apoptosis due to chromosome separation failure at metaphase I. Therefore, the alloheptaploid females and males were all sterile. Finally, we established a sustainable clone for the large-scale production of NA3n♀II and developed an efficient approach to synthesize diverse allopolyploids containing genomes of different cyprinid species. These findings not only broaden our understanding of reproduction transition but also offer a practical strategy for polyploidy breeding and heterosis fixing.

Characteristics and expression profiles of MHC class Ⅰ molecules in Carassius auratus

Major histocompatibility complex class Ⅰ (MHC Ⅰ) molecules play a vital role in adaptive immune systems in vertebrates by presenting antigens to effector T cells. Understanding the expression profiling of MHC Ⅰ molecules in fish is essential for improving our knowledge of the relationship between microbial infection and adaptive immunity. In this study, we conducted a comprehensive analysis of MHC Ⅰ gene characteristics in Carassius auratus, an important freshwater aquaculture fish in China that is susceptible to Cyprinid herpesvirus 2 (CyHV-2) infection. We identified approximately 20 MHC Ⅰ genes discussed, including U, Z, and L lineage genes. However, only U and Z lineage proteins were identified in the kidney of Carassius auratus using high pH reversed-phase chromatography and mass spectrometry. The L lineage proteins were either not expressed or present at an extremely low level in the kidneys of Carassius auratus. We also used targeted proteomics to analyze changes in protein MHC Ⅰ molecules abundance in healthy and CyHV-2-infected Carassius auratus. We observed that five MHC Ⅰ molecules were upregulated, and Caau-UFA was downregulated in the diseased group. This study is the first to reveal the expression of MHC Ⅰ molecules at a large scale in Cyprinids, which enhances our understanding of fish adaptive immune systems.

The hagfish genome and the evolution of vertebrates

As the only surviving lineages of jawless fishes, hagfishes and lampreys provide a critical window into early vertebrate evolution. Here, we investigate the complex history, timing, and functional role of genome-wide duplications in vertebrates in the light of a chromosome-scale genome of the brown hagfish Eptatretus atami. Using robust chromosome-scale (paralogon-based) phylogenetic methods, we confirm the monophyly of cyclostomes, document an auto-tetraploidization (1RV) that predated the origin of crown group vertebrates ~517 Mya, and establish the timing of subsequent independent duplications in the gnathostome and cyclostome lineages. Some 1RV gene duplications can be linked to key vertebrate innovations, suggesting that this early genomewide event contributed to the emergence of pan-vertebrate features such as neural crest. The hagfish karyotype is derived by numerous fusions relative to the ancestral cyclostome arrangement preserved by lampreys. These genomic changes were accompanied by the loss of genes essential for organ systems (eyes, osteoclast) that are absent in hagfish, accounting in part for the simplification of the hagfish body plan; other gene family expansions account for hagfishes' capacity to produce slime. Finally, we characterise programmed DNA elimination in somatic cells of hagfish, identifying protein-coding and repetitive elements that are deleted during development. As in lampreys, the elimination of these genes provides a mechanism for resolving genetic conflict between soma and germline by repressing germline/pluripotency functions. Reconstruction of the early genomic history of vertebrates provides a framework for further exploration of vertebrate novelties.

Triclosan and triclocarban weaken the olfactory capacity of goldfish by constraining odorant recognition, disrupting olfactory signal transduction, and disturbing olfactory information processing

Recently, increased production and consumption of disinfectants such as triclosan (TCS) and triclocarban (TCC) have led to massive pollution of the environment, which draws global concern over the potential risk to aquatic organisms. However, the olfactory toxicity of disinfectants in fish remains elusive to date. In the present study, the impact of TCS and TCC on the olfactory capacity of goldfish was assessed by neurophysiological and behavioral approaches. As shown by the reduced distribution shifts toward amino acid stimuli and hampered electro-olfactogram responses, our results demonstrated that TCS/TCC treatment would cause deterioration of the olfactory ability of goldfish. Our further analysis found that TCS/TCC exposure suppressed the expression of olfactory G protein-coupled receptors in the olfactory epithelium, restricted the transformation of odorant stimulation into electrical responses by disturbing the cAMP signaling pathway and ion transportation, and induced apoptosis and inflammation in the olfactory bulb. In conclusion, our results demonstrated that an environmentally realistic level of TCS/TCC would weaken the olfactory capacity of goldfish by constraining odorant recognition efficiency, disrupting olfactory signal generation and transduction, and disturbing olfactory information processing.

Promoter evolution of mammalian gene duplicates

BACKGROUND

Gene duplication is thought to be a central process in evolution to gain new functions. The factors that dictate gene retention following duplication as well paralog gene divergence in sequence, expression and function have been extensively studied. However, relatively little is known about the evolution of promoter regions of gene duplicates and how they influence gene duplicate divergence. Here, we focus on promoters of paralog genes, comparing their similarity in sequence, in the sets of transcription factors (TFs) that bind them, and in their overall promoter architecture.

RESULTS

We observe that promoters of recent duplications display higher sequence similarity between them and that sequence similarity rapidly declines between promoters of more ancient paralogs. In contrast, similarity in cis-regulation, as measured by the set of TFs that bind promoters of both paralogs, does not simply decrease with time from duplication and is instead related to promoter architecture-paralogs with CpG Islands (CGIs) in their promoters share a greater fraction of TFs, while CGI-less paralogs are more divergent in their TF binding set. Focusing on recent duplication events and partitioning them by their duplication mechanism enables us to uncover promoter properties associated with gene retention, as well as to characterize the evolution of promoters of newly born genes: In recent retrotransposition-mediated duplications, we observe asymmetry in cis-regulation of paralog pairs: Retrocopy genes are lowly expressed and their promoters are bound by fewer TFs and are depleted of CGIs, in comparison with the original gene copy. Furthermore, looking at recent segmental duplication regions in primates enable us to compare successful retentions versus loss of duplicates, showing that duplicate retention is associated with fewer TFs and with CGI-less promoter architecture.

CONCLUSIONS

In this work, we profiled promoters of gene duplicates and their inter-paralog divergence. We also studied how their characteristics are associated with duplication time and duplication mechanism, as well as with the fate of these duplicates. These results underline the importance of cis-regulatory mechanisms in shaping the evolution of new genes and their fate following duplication.

Disruption of T-box transcription factor eomesa results in abnormal development of median fins in Oujiang color common carp Cyprinus carpio

Median fins are thought to be ancestors of paired fins which in turn give rise to limbs in tetrapods. However, the developmental mechanisms of median fins remain largely unknown. Nonsense mutation of the T-box transcription factor eomesa in zebrafish results in a phenotype without dorsal fin. Compared to zebrafish, the common carp undergo an additional round of whole genome duplication, acquiring an extra copy of protein-coding genes. To verify the function of eomesa genes in common carp, we established a biallelic gene editing technology in this tetraploidy fish through simultaneous disruption of two homologous genes, eomesa1 and eomesa2. We targeted four sites located upstream or within the sequences encoding the T-box domain. Sanger sequencing data indicated the average knockout efficiency was around 40% at T1-T3 sites and 10% at T4 site in embryos at 24 hours post fertilization. The individual editing efficiency was high to about 80% at T1-T3 sites and low to 13.3% at T4 site in larvae at 7 days post fertilization. Among 145 mosaic F0 examined at four months old, three individuals (Mutant 1-3) showed varying degrees of maldevelopment in the dorsal fin and loss of anal fin. Genotyping showed the genomes of all three mutants were disrupted at T3 sites. The null mutation rates on the eomesa1 and eomesa2 loci were 0% and 60% in Mutant 1, 66.7% and 100% in Mutant 2, and 90% and 77.8% in Mutant 3, respectively. In conclusion, we demonstrated a role of eomesa in the formation and development of median fins in Oujiang color common carp and established an method that simultaneously disrupt two homologous genes with one gRNA, which would be useful in genome editing in other polyploidy fishes.

Full-Length RNA Sequencing Provides Insights into Goldfish Evolution under Artificial Selection

Goldfish Carassius auratus is an ideal model for exploring fish morphology evolution. Although genes underlying several ornamental traits have been identified, little is known about the effects of artificial selection on embryo gene expression. In the present study, hybrid transcriptome sequencing was conducted to reveal gene expression profiles of Celestial-Eye (CE) and Ryukin (RK) goldfish embryos. Full-length transcriptome sequencing on the PacBio platform identified 54,218 and 54,106 transcript isoforms in CE and RK goldfish, respectively. Of particular note was that thousands of alternative splicing (AS) and alternative polyadenylation (APA) events were identified in both goldfish breeds, and most of them were inter-breed specific. RT-PCR and Sanger sequencing showed that most of the predicted AS and APA were correct. Moreover, abundant long non-coding RNA and fusion genes were detected, and again most of them were inter-breed specific. Through RNA-seq, we detected thousands of differentially expressed genes (DEGs) in each embryonic stage between the two goldfish breeds. KEGG enrichment analysis on DEGs showed extensive differences between CE and RK goldfish in gene expression. Taken together, our results demonstrated that artificial selection has led to far-reaching influences on goldfish gene expression, which probably laid the genetic basis for hundreds of goldfish variations.

The miR-430 locus with extreme promoter density forms a transcription body during the minor wave of zygotic genome activation

In anamniote embryos, the major wave of zygotic genome activation starts during the mid-blastula transition. However, some genes escape global genome repression, are activated substantially earlier, and contribute to the minor wave of genome activation. The mechanisms underlying the minor wave of genome activation are little understood. We explored the genomic organization and cis-regulatory mechanisms of a transcription body, in which the minor wave of genome activation is first detected in zebrafish. We identified the miR-430 cluster as having excessive copy number and the highest density of Pol-II-transcribed promoters in the genome, and this is required for forming the transcription body. However, this transcription body is not essential for, nor does it encompasse, minor wave transcription globally. Instead, distinct minor-wave-specific promoter architecture suggests that promoter-autonomous mechanisms regulate the minor wave of genome activation. The minor-wave-specific features also suggest distinct transcription initiation mechanisms between the minor and major waves of genome activation.

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.

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.

Transcriptome profiling and differential expression analysis of the immune-related genes during the acute phase of infection with Photobacterium damselae subsp. damselae in silver pomfret (Pampus argenteus)

Silver pomfret has been widely cultured in China due to its high economic value. Photobacterium damselae subsp. damselae (PDD) is a Gram-negative bacterium that has been shown to infect many fish species. To increase knowledge of the molecular mechanisms of the host defense against PDD, we conducted transcriptome analysis of head kidney in silver pomfret at 24 h and 72 h post-infection (hpi) via Illumina sequencing. The de novo assembly resulted in the identification of 79,063 unigenes, with 59,386 (75.11%) successfully annotated in public databases (NR, NT, KO, Swiss-Prot, Pfam, GO, and KOG databases). Comparison of gene expression profiles between PBS-injected fish (sham control) and PDD-challenged fish revealed 329 and 570 differentially expressed genes (DEGs) were screened at 24 hpi and 72 hpi, respectively. The DEGs were enriched in multiple immune-related pathways such as Hepatitis C, Gastric acid secretion, CAMs and Leukocyte transendothelial migration pathways, Primary immunodeficieny, ECM-receptor interaction, PI3K-Akt signaling pathway. The data obtained in the present study offers valuable information for acute immune response of silver pomfret challenged with PDD, which will facilitate further investigations on strategies against Photobacterium spp. infection in teleosts.

Lactate-mediated neural plasticity genes emerged during the evolution of memory systems

The ability to record experiences and learning is present to different degrees in several species; however, the complexity and diversity of memory processes are cognitive function features that differentiate humans from other species. Lactate has recently been discovered to act as a signaling molecule for neuronal plasticity linked to long-term memory. Because lactate is not only an energy substrate for neurons but also a signaling molecule for plasticity (Magistretti and Allaman in Nat Rev Neurosci 19:235-249, 2018. https://doi.org/10.1038/nrn.2018.19 ), it is of particular interest to understand how and when memory-related genes and lactate-mediated neural plasticity (LMNP) genes emerged and evolved in humans. To understand the evolutionary origin and processes of memory and LMNP genes, we first collected information on genes related to memory and LMNP from the literature and then conducted a comparative analysis of these genes. We found that the memory and LMNP genes have different origins, suggesting that these genes may have become established gradually in evolutionarily and functional terms and not at the same time. We also found that memory and LMNP systems have a similar evolutionary history, having been formed with the gradual participation of newly emerging genes throughout their evolution. We propose that the function of LMNP as a signaling process may be evolutionarily associated with memory systems through an unidentified system that is linked by 13 common genes between memory and LMNP gene sets. This study provides evolutionary insight into the possible relationship between memory and the LMNP systems that deepens our understanding of the evolution of memory systems.

Genetic and epigenetic regulation of growth, reproduction, disease resistance and stress responses in aquaculture

Major progress has been made with genomic and genetic studies in aquaculture in the last decade. However, research on epigenetic regulation of aquaculture traits is still at an early stage. It is apparent that most, if not all, aquaculture traits are regulated at both genetic and epigenetic levels. This paper reviews recent progress in understanding of genetic and epigenetic regulation of important aquaculture traits such as growth, reproduction, disease resistance, and stress responses. Although it is challenging to make generalized statements, DNA methylation is mostly correlated with down-regulation of gene expression, especially when at promoters and enhancers. As such, methylation of growth factors and their receptors is negatively correlated with growth; hypomethylation of genes important for stress tolerance is correlated with increased stress tolerance; hypomethylation of genes important for male or female sex differentiation leads to sex differentiation into males or females, respectively. It is apparent that environmental regulation of aquaculture traits is mediated at the level of epigenetic regulation, and such environment-induced epigenetic changes appeared to be intergenerationally inherited, but evidences for transgenerational inheritance are still limited.

Cell Junction and Vesicle Trafficking-Mediated Melanosome/Melanin Transfer Are Involved in the Dynamic Transformation of Goldfish Carassius auratus Skin Color

Goldfish are one of the most popular models for studying the genetic diversity of skin color. Transcriptome sequencing (RNA-seq) and whole genome bisulfate sequencing (WGBS) of skin tissues from the third filial (F3) cyan (CN), black (BK), and white (WH) goldfish were conducted to analyze the molecular mechanism of color transformation in fish. The RNA-seq yielded 56 Gb of clean data and 56,627 transcripts from nine skin samples. The DEGs (differentially expressed genes) were enriched in cell junction cellular components and the tight junction pathway. Ninety-five homologs of the claudin family were predicted and 16 claudins were identified in correlation with skin color transformation. WGBS yielded 1079 Gb of clean data from 15 samples. Both the DEGs and the DMRs (differentially methylated regions) in the BK_CN group were found to be enriched in cytoskeleton reorganization and vesicle trafficking. Masson staining and TEM (transmission electron microscopy) confirmed the varied distribution and processes of melanosome/melanin in skin tissues. Our results suggested that cytoskeleton reorganization, cell junction, and the vesicle trafficking system played key roles in the transfer of the melanosome/melanin, and it was the extracellular translocation rather than the biosynthesis or metabolism of the melanin process that resulted in the color transformation of cyan goldfish. The data will facilitate the understanding of the molecular mechanisms underlying dynamic skin color transformation in goldfish.

Immunoglobulin heavy-chain loci in ancient allotetraploid goldfish

As an ancient allotetraploid species, goldfish (Carassius auratus) have two sets of subgenomes. In this study, immunoglobulin heavy-chain (IGH) genes were cloned from the red crucian carp (Carassius auratus red var.), and the corresponding loci were identified in the gynogenetic diploid red crucian carp (GRCC) genome as well as the genomes of three other goldfish strains (Wakin, G-12, and CaTCV-1). Examination showed that each goldfish strain possessed two sets of parallel IGH loci: a complete IGHA locus and a degenerated IGHB locus that was nearly 40 × smaller. In the IGHA locus, multiple τ chain loci were arranged in tandem between the μ&δ chain locus and the variable genes, but no τ-like genes were found in the IGHB locus.

Comparative genomics analysis and genome assembly integration with the recombination landscape contribute to Takifugu bimaculatus assembly refinement

Takifugu genus has been brought to the fore in scientific and practical research due to its compact genome, explosive speciation progress and economic value. Here we updated the chromosome-level genome of Takifugu bimaculatus by an ultra-high-density linkage map, a classic and accurate way of chromosome assembly. The map constituted a robust assembly frame, with 92.2% (372.77 Mb) of the draft genome cumulatively placed. With intraspecies and interspecies comparative genomic analysis, we developed a criterion to quantify the differences between assemblies and established a novel way to integrate information from multiple assemblies. The integrated assembly rectified potential mis-assemblies, greatly improving the genome contiguity and correctness. Our results rendered profound information on the genetic recombination of T. bimaculatus and provided new insights into effective genome assembly. The consolidated assembly will be a contributory tool of T. bimaculatus and broadly across the Takifugu by providing a convincing reference for genomic research.

Breaking the reproductive barrier of divergent species to explore the genomic landscape

Background: Climate change will have significant consequences for species. Species range shifts induce the emergence of new hybrid zones or the spatial displacement of pre-existing ones. These hybrid zones may become more porous as alleles are passed from one species to another. Currently, hybridization between highly divergent species living in sympatry seems extremely limited. Indeed, this phenomenon involves breaking two barriers. The first is the pre-mating barrier, related to the reproductive phenology of the two species. The second is the post-zygotic barrier, related to the genetic divergence between these species. Here, we were interested in identifying new hybridization patterns and potential implications, especially in the context of environmental modifications.

Methods: We sampled Telestes souffia and Parachondrostoma toxostoma wild specimens from different locations across France and genotyped them for SNP markers. We identified discriminant loci using F1-hybrid specimens and parental species and performed principal component analysis and Bayesian model-based clustering to analyze phylogenetic information. Furthermore, we assessed deviation in allele frequency from F1 to F2 and for Hardy–Weinberg equilibrium for F2 and assessed gene function associated with two F2 cohorts.

Results: We demonstrate that by breaking the ecological barrier, massive introgressive hybridization is possible between two endemic lineages of Cyprinidae belonging to two distinct genera. For both cohorts studied (=2 cm and >2 cm), a large majority of loci (>88%) presented no deviation in allele frequency and no departure from the Hardy–Weinberg equilibrium. For individuals beyond the 2 cm stage, two phenomena were observed. The first was an allelic imbalance in favor of P. toxostoma, for some genomic regions, with genes involved in developmental regulatory processes, cytoskeletal organization, and chromosome organization. The second was an excess of heterozygous loci coupled with an equilibrium of allelic frequencies for genes involved in immune response and kidney/liver development. Moreover, the 2 cm-sized specimens with high mortality yielded a particular genomic signature.

Conclusion: Our study displayed important results for understanding the early stages of hybridization between divergent lineages and predicting the emergence of future hybrid zones in the wild. Moreover, this hybridization generates a wide spectrum of hybrids that are a potential source of important evolutionary novelties.

A regulatory network of Sox and Six transcription factors initiate a cell fate transformation during hearing regeneration in adult zebrafish

Using adult zebrafish inner ears as a model for sensorineural regeneration, we ablated the mechanosensory receptors and characterized the single-cell epigenome and transcriptome at consecutive time points during hair cell regeneration. We utilized deep learning on the regeneration-induced open chromatin sequences and identified cell-specific transcription factor (TF) motif patterns. Enhancer activity correlated with gene expression and identified potential gene regulatory networks. A pattern of overlapping Sox- and Six-family TF gene expression and binding motifs was detected, suggesting a combinatorial program of TFs driving regeneration and cell identity. Pseudotime analysis of single-cell transcriptomic data suggested that support cells within the sensory epithelium changed cell identity to a “progenitor” cell population that could differentiate into hair cells. We identified a 2.6 kb DNA enhancer upstream of the sox2 promoter that, when deleted, showed a dominant phenotype that resulted in a hair-cell-regeneration-specific deficit in both the lateral line and adult inner ear.

Jimenez et al. interrogate the epigenomic and transcriptomic landscape of regenerating adult zebrafish inner-ear sensory epithelia. They show that the support-cell population transitions to an intermediate “progenitor” cell state that becomes new hair cells, and they demonstrate that the cell fate decisions may be driven by the coordinate regulation and spatial co-binding of Sox and Six transcription factors. By functionally validating a predicted regeneration-responsive enhancer upstream of sox2, they show that precise timing of sox2 expression is critical for hearing regeneration in zebrafish.

Changes in Ploidy Drive Reproduction Transition and Genomic Diversity in a Polyploid Fish Complex

Unisexual animals are commonly found in some polyploid species complexes, and most of these species have had a long evolutionary history. However, their method for avoiding genomic decay remains unclear. The polyploid Carassius complex naturally comprises the sexual amphidiploid C. auratus (crucian carp or goldfish) (AABB) and the gynogenetic amphitriploid C. gibelio (gibel carp) (AAABBB). Recently, we developed a fertile synthetic amphitetraploid (AAAABBBB) male from C. gibelio by incorporating a C. auratus genome. In this study, we generated novel amphitriploids (AAABBB) by backcrossing the amphitetraploid male with the amphidiploid C. auratus. Whole-genome resequencing revealed the genomic changes, including recombination and independent assortment between homologs of C. gibelio and C. auratus. The fertility, sex determination system, oocyte development, and fertilization behaviors of the novel amphitriploids were investigated. Approximately 80% of the novel amphitriploid females recovered the unisexual gynogenesis ability. Intriguingly, two types of primary oocyte (with and without homolog synapsis) were discovered, and their distinct development fates were observed. Type I oocytes entered apoptosis due to improper synaptonemal complex assembly and incomplete double-strand break repair, whereas subsequent type II oocytes bypassed meiosis through an alternative ameiotic pathway to develop into mature eggs. Moreover, gynogenesis was stabilized in their offspring, and a new array of diverse gynogenetic amphitriploid clones was produced. These revealed genomic changes and detailed cytological data provide comprehensive evidence that changes in ploidy drive unisexual and sexual reproduction transition, thereby resulting in genomic diversity and allowing C. gibelio avoid genomic decay.

DupScan: predicting and visualizing vertebrate genome duplication database

Duplicated genes prevail in vertebrates and are important in the acquisition of new genes and novelties. Whole genome duplication (WGD) is one of the sources of duplicated genes. It can provide raw materials for natural selection by increasing the flexibility and complexity of the genome. WGDs are the driving force for the evolution of vertebrates and contribute greatly to their species diversity, especially in fish species with complicated WGD patterns. Here, we constructed the DupScan database (https://dupscan.sysumeg.com/) by integrating 106 chromosomal-level genomes, which can analyze and visualize synteny at both the gene and genome scales, visualize the Ka, Ks, and 4DTV values, and browse genomes. DupScan was used to perform functional adaptation for the intricate WGD investigation based on synteny matching. DupScan supports the analysis of five WGD rounds (R): VGD2 (vertebrate genome duplication 2), Ars3R (Acipenser-ruthenus-specific 3R), Pss3R (Polyodon-spathula-specific 3R), Ts3R (teleost-specific duplication 3R), Ss4R (salmonid-specific 4R), and Cs4R (carp-specific 4R). DupScan serves as one-stop analysis platform for synteny and WGD research in which users can analyze and predict synteny and WGD patterns across 106 species of whole genome sequences. This further aided us in elucidating genome evolutionary patterns across over 60,000 vertebrate species with synteny and WGD events.