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

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.

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.

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.

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.

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.

Genome resequencing clarifies phylogeny and reveals patterns of selection in the toxicogenomics model Pimephales promelas

Background

The fathead minnow (Pimephales promelas) is a model species for toxicological research. A high-quality genome reference sequence is available, and genomic methods are increasingly used in toxicological studies of the species. However, phylogenetic relationships within the genus remain incompletely known and little population-genomic data are available for fathead minnow despite the potential effects of genetic background on toxicological responses. On the other hand, a wealth of extant samples is stored in museum collections that in principle allow fine-scale analysis of contemporary and historical genetic variation.

Methods

Here we use short-read shotgun resequencing to investigate sequence variation among and within Pimephales species. At the genus level, our objectives were to resolve phylogenetic relationships and identify genes with signatures of positive diversifying selection. At the species level, our objective was to evaluate the utility of archived-sample resequencing for detecting selective sweeps within fathead minnow, applied to a population introduced to the San Juan River of the southwestern United States sometime prior to 1950.

Results

We recovered well-supported but discordant phylogenetic topologies for nuclear and mitochondrial sequences that we hypothesize arose from mitochondrial transfer among species. The nuclear tree supported bluntnose minnow (P. notatus) as sister to fathead minnow, with the slim minnow (P. tenellus) and bullhead minnow (P. vigilax) more closely related to each other. Using multiple methods, we identified 11 genes that have diversified under positive selection within the genus. Within the San Juan River population, we identified selective-sweep regions overlapping several sets of related genes, including both genes that encode the giant sarcomere protein titin and the two genes encoding the MTORC1 complex, a key metabolic regulator. We also observed elevated polymorphism and reduced differentation among populations (F_ST) in genomic regions containing certain immune-gene clusters, similar to what has been reported in other taxa. Collectively, our data clarify evolutionary relationships and selective pressures within the genus and establish museum archives as a fruitful resource for characterizing genomic variation. We anticipate that large-scale resequencing will enable the detection of genetic variants associated with environmental toxicants such as heavy metals, high salinity, estrogens, and agrichemicals, which could be exploited as efficient biomarkers of exposure in natural populations.

Potential mechanisms of Na+/K+-ATPase attenuation by heat and pesticides co-exposure in goldfish: role of cellular apoptosis, oxidative/nitrative stress, and antioxidants in gills

In this study, we examined the dose-dependent effects of an environmentally relevant pesticide cocktail (metalachlor, linuron, isoproturon, tebucanazole, aclonifen, atrazine, pendimethalin, and azinphos-methyl) and temperature change (22 vs. 32 °C for 4-week exposure) on Na⁺/K⁺-ATPase, 3-nitrotyrosine protein (NTP), dinitrophenyl protein (DNP), catalase (CAT), and superoxide dismutase (SOD) expressions in gills of goldfish (Carassius auratus). Histopathological analysis showed widespread damage to gill in elevated temperature (32 °C) and pesticide co-exposure groups, including fusion of secondary lamellae, club-shaped primary lamellae, rupture of epithelial layer, loss of normal architecture, and hemorrhaging. Immunohistochemical and qRT-PCR analyses showed significant decreases in Na⁺/K⁺-ATPase protein and mRNA expressions in gills exposed to higher temperature and pesticides; however, combined exposure to heat and pesticides significantly increases NTP, DNP, CAT, and SOD expressions. In situ TUNEL assay revealed elevated levels of apoptotic cells in response to combined exposure. Collectively, our results suggest the combined effects of heat and pesticide stress cause cellular damage, upregulate oxidative/nitrative stress biomarkers, and increase apoptotic cells, downregulate Na⁺/K⁺-ATPase expression in gills. This provides new evidence for oxidant/antioxidant-dependent mechanisms for downregulation of Na⁺/K⁺-ATPase expression in gills during combined exposure.

Equilibrated evolution of the mixed auto-/allopolyploid haplotype-resolved genome of the invasive hexaploid Prussian carp

Understanding genome evolution of polyploids requires dissection of their often highly similar subgenomes and haplotypes. Polyploid animal genome assemblies so far restricted homologous chromosomes to a ‘collapsed’ representation. Here, we sequenced the genome of the asexual Prussian carp, which is a close relative of the goldfish, and present a haplotype-resolved chromosome-scale assembly of a hexaploid animal. Genome-wide comparisons of the 150 chromosomes with those of two ancestral diploid cyprinids and the allotetraploid goldfish and common carp revealed the genomic structure, phylogeny and genome duplication history of its genome. It consists of 25 syntenic, homeologous chromosome groups and evolved by a recent autoploid addition to an allotetraploid ancestor. We show that de-polyploidization of the alloploid subgenomes on the individual gene level occurred in an equilibrated fashion. Analysis of the highly conserved actinopterygian gene set uncovered a subgenome dominance in duplicate gene loss of one ancestral chromosome set.

The haplotype-resolved assembly of the asexual invasive Prussian carp shows six genome copies (AAABBB), evolved from two ancestral species by a recent self-addition (AB) to its hybrid-tetraploid (AABB) goldfish ancestor. Equilibrated gene loss led to subgenome dominance.

Comparative genome anatomy reveals evolutionary insights into a unique amphitriploid fish

Triploids are rare in nature because of difficulties in meiotic and gametogenic processes, especially in vertebrates. The Carassius complex of cyprinid teleosts contains sexual tetraploid crucian carp/goldfish (C. auratus) and unisexual hexaploid gibel carp/Prussian carp (C. gibelio) lineages, providing a valuable model for studying the evolution and maintenance mechanism of unisexual polyploids in vertebrates. Here we sequence the genomes of the two species and assemble their haplotypes, which contain two subgenomes (A and B), to the chromosome level. Sequencing coverage analysis reveals that C. gibelio is an amphitriploid (AAABBB) with two triploid sets of chromosomes; each set is derived from a different ancestor. Resequencing data from different strains of C. gibelio show that unisexual reproduction has been maintained for over 0.82 million years. Comparative genomics show intensive expansion and alterations of meiotic cell cycle-related genes and an oocyte-specific histone variant. Cytological assays indicate that C. gibelio produces unreduced oocytes by an alternative ameiotic pathway; however, sporadic homologous recombination and a high rate of gene conversion also exist in C. gibelio. These genomic changes might have facilitated purging deleterious mutations and maintaining genome stability in this unisexual amphitriploid fish. Overall, the current results provide novel insights into the evolutionary mechanisms of the reproductive success in unisexual polyploid vertebrates.

Genome sequencing and haplotype assembly of two cyprinid teleosts, a sexual tetraploid and an unisexual hexaploid, reveal insights into the evolutionary mechanisms underpinning the reproductive success of unisexual polyploid vertebrates.

Multiomic atlas with functional stratification and developmental dynamics of zebrafish cis-regulatory elements

Zebrafish, a popular organism for studying embryonic development and for modeling human diseases, has so far lacked a systematic functional annotation program akin to those in other animal models. To address this, we formed the international DANIO-CODE consortium and created a central repository to store and process zebrafish developmental functional genomic data. Our data coordination center ( https://danio-code.zfin.org ) combines a total of 1,802 sets of unpublished and re-analyzed published genomic data, which we used to improve existing annotations and show its utility in experimental design. We identified over 140,000 cis-regulatory elements throughout development, including classes with distinct features dependent on their activity in time and space. We delineated the distinct distance topology and chromatin features between regulatory elements active during zygotic genome activation and those active during organogenesis. Finally, we matched regulatory elements and epigenomic landscapes between zebrafish and mouse and predicted functional relationships between them beyond sequence similarity, thus extending the utility of zebrafish developmental genomics to mammals.

Two duplicated gsdf homeologs cooperatively regulate male differentiation by inhibiting cyp19a1a transcription in a hexaploid fish

Although evolutionary fates and expression patterns of duplicated genes have been extensively investigated, how duplicated genes co-regulate a biological process in polyploids remains largely unknown. Here, we identified two gsdf (gonadal somatic cell-derived factor) homeologous genes (gsdf-A and gsdf-B) in hexaploid gibel carp (Carassius gibelio), wherein each homeolog contained three highly conserved alleles. Interestingly, gsdf-A and gsdf-B transcription were mainly activated by dmrt1-A (dsx- and mab-3-related transcription factor 1) and dmrt1-B, respectively. Loss of either gsdf-A or gsdf-B alone resulted in partial male-to-female sex reversal and loss of both caused complete sex reversal, which could be rescued by a nonsteroidal aromatase inhibitor. Compensatory expression of gsdf-A and gsdf-B was observed in gsdf-B and gsdf-A mutants, respectively. Subsequently, we determined that in tissue culture cells, Gsdf-A and Gsdf-B both interacted with Ncoa5 (nuclear receptor coactivator 5) and blocked Ncoa5 interaction with Rora (retinoic acid-related orphan receptor-alpha) to repress Rora/Ncoa5-induced activation of cyp19a1a (cytochrome P450, family 19, subfamily A, polypeptide 1a). These findings illustrate that Gsdf-A and Gsdf-B can regulate male differentiation by inhibiting cyp19a1a transcription in hexaploid gibel carp and also reveal that Gsdf-A and Gsdf-B can interact with Ncoa5 to suppress cyp19a1a transcription in vitro. This study provides a typical case of cooperative mechanism of duplicated genes in polyploids and also sheds light on the conserved evolution of sex differentiation.

Aquaculture omics: An update on the current status of research and data analysis

Aquaculture is the fast-growing agricultural sector and has the ability to meet the growing demand for protein nutritional security for future population. In future aquaculture is going to be the major source of fish proteins as capture fisheries reached at its maximum. However, several challenges need to overcome such as lack of genetically improved strains/varieties, lack of species-specific feed/functional feed, round the year availability of quality fish seed, pollution of ecosystems and increased frequencies of disease occurrence etc. In recent years, the continuous development of high throughput sequencing technology has revolutionized the biological sciences and provided necessary tools. Application of 'omics' in aquaculture research have been successfully used to resolve several productive and reproductive issues and thus ensure its sustainability and profitability. To date, high quality draft genomes of over fifty fish species have been generated and successfully used to develop large number of single nucleotide polymorphism markers (SNPs), marker panels and other genomic resources etc in several aquaculture species. Similarly, transcriptome profiling and miRNAs analysis have been used in aquaculture research to identify key transcripts and expression analysis of candidate genes/miRNAs involved in reproduction, immunity, growth, development, stress toxicology and disease. Metagenome analysis emerged as a promising scientific tool to analyze the complex genomes contained within microbial communities. Metagenomics has been successfully used in the aquaculture sector to identify novel and potential pathogens, antibiotic resistance genes, microbial roles in microcosms, microbial communities forming biofloc, probiotics etc. In the current review, we discussed application of high-throughput technologies (NGS) in the aquaculture sector.

A Novel TRIM9 Protein Promotes NF-κB Activation Through Interacting With LvIMD in Shrimp During WSSV Infection

The TRIpartite Motif (TRIM) proteins play key roles in cell differentiation, apoptosis, development, autophagy, and innate immunity in vertebrates. In the present study, a novel TRIM9 homolog (designated as LvTRIM9-1) specifically expressed in the lymphoid organ of shrimp was identified from the Pacific whiteleg shrimp Litopenaeus vannamei. Its deduced amino acid sequence possesses the typical features of TRIM proteins, including a RING domain, two B-boxes, a coiled-coil domain, a FN3 domain, and a SPRY domain. The transcripts of LvTRIM9-1 were mainly located in the lymphoid tubules of the lymphoid organ. Knockdown of LvTRIM9-1 could apparently inhibit the transcriptions of some genes from white spot syndrome virus (WSSV) and reduce the viral propagation in the lymphoid organ. Overexpression of LvTRIM9-1 in mammalian cells could activate the promoter activity of NF-κB, and an in vivo experiment in shrimp showed that knockdown of LvTRIM9-1 reduced the expression of LvRelish in the lymphoid organ. Yeast two-hybridization and co-immunoprecipitation (Co-IP) assays confirmed that LvTRIM9-1 could directly interact with LvIMD, a key component of the IMD pathway, through its SPRY domain. These data suggest that LvTRIM9-1 could activate the IMD pathway in shrimp via interaction with LvIMD. This is the first evidence to show the regulation of a TRIM9 protein on the IMD pathway through its direct interaction with IMD, which will enrich our knowledge on the role of TRIM proteins in innate immunity of invertebrates.