Mandarin fish (Sinipercidae) genomes provide insights into innate predatory feeding

Assembly and analysis of Sinipercidae fish sex chromosomes reveals that a supergene drives sex chromosome origin and turnover

The sex chromosomes of fish are usually homomorphic and recently derived, making them an ideal model to understand the origin and evolution of sex chromosomes in vertebrates. Here, combined Pacbio, Hi-C, and Illumina sequencing of the male mandarin fish (Siniperca chuatsi) genome enabled the assembly of a high-quality chromosome-level genome, including highly resolved X and Y chromosome assemblies. Genome wide association study (GWAS) and coverage information analysis revealed a 2.0 Mbp sex determining region (SDR). A Y-specific duplication of the anti-Müllerian hormone (amh) gene, amhy, was identified that encoded a truncated AMH protein. Loss of function and overexpression experiments demonstrated that amhy may act as a male sex-determining gene (SDG). GWAS led to the identification of a common sex chromosome (Chr24) and SDG (amhy) in S. scherzeri, while another sex chromosome (Chr11) was identified in Coreoperca whiteheadi. Interestingly, completely conserved protein coding sequences of amhy were commonly identified in male S. chuatsi, S. scherzeri, and S. knerii. These results support a classic model for the origin and evolution of early sex chromosomes and suggest that amhy is a super SDG that can drive the origin and turnover of sex chromosome in vertebrates.

WASTER: Practical de novo phylogenomics from low-coverage short reads

The advent of affordable whole-genome sequencing has spurred numerous large-scale projects aimed at inferring the tree of life, yet achieving a complete species-level phylogeny remains a distant goal due to significant costs and computational demands. Traditional species tree inference methods, though effective, are hampered by the need for high-coverage sequencing, high-quality genomic alignments, and extensive computational resources. To address these challenges, this study introduces WASTER, a novel de novo tool for inferring species trees directly from short-read sequences. WASTER employs a k-mer based approach for identifying variable sites, circumventing the need for genome assembly and alignment. Using simulations, we demonstrate that WASTER achieves accuracy comparable to that of traditional alignment-based methods, even for low sequencing depth, and has substantially higher accuracy than other alignment-free methods. We validate WASTER's efficacy on real data, where it accurately reconstructs phylogenies of eukaryotic species with as low depth as 1.5X. WASTER provides a fast and efficient solution for phylogeny estimation in cases where genome assembly and/or alignment may bias analyses or is challenging, for example due to low sequencing depth. It also provides a method for generating guide trees for tree-based alignment algorithms. WASTER's ability to accurately estimate trees from low-coverage sequencing data without relying on assembly and alignment will lead to substantially reduced sequencing and computational costs in phylogenomic projects.

Dynamic Outlier Slicing Allows Broader Exploration of Adaptive Divergence: A Comparison of Individual Genome and Pool-Seq Data Linked to Humic Adaptation in Perch

How genetic variation contributes to adaptation at different environments is a central focus in evolutionary biology. However, most free-living species still lack a comprehensive understanding of the primary molecular mechanisms of adaptation. Here, we characterised the targets of selection associated with drastically different aquatic environments-humic and clear water-in the common freshwater fish, Eurasian perch (Perca fluviatilis). By using whole-genome sequencing (WGS) on a large population dataset (n = 42 populations) and analysing 873,788 SNPs, our primary aim was to uncover novel and confirm known footprints of selection. We compared individual and pooled WGS, and developed a novel approach, termed dynamic outlier slicing, to assess how the choice of outlier-calling stringency influences functional and Gene Ontology (GO) enrichment. By integrating genome-environment association (GEA) analysis with allele frequency-based approaches, we estimated composite selection signals (CSS) and identified 2679 outlier SNPs distributed across 324 genomic regions, involving 468 genes. Dynamic outlier slicing identified robust enrichment signals in five annotation categories (upstream, downstream, synonymous, 5'UTR and 3'UTR) highlighting the crucial role of regulatory elements in adaptive evolution. Furthermore, GO analyses revealed strong enrichment of molecular functions associated with gated channel activity, transmembrane transporter activity and ion channel activity, emphasising the importance of osmoregulation and ion balance maintenance. Our findings demonstrate that despite substantial random drift and divergence, WGS of high number of population pools enabled the identification of strong selection signals associated with adaptation to both humic and clear water environments, providing robust evidence of widespread adaptation. We anticipate that the dynamic outlier slicing method we developed will enable a more thorough exploration of adaptive divergence across a diverse range of species.

Genomes reveal pervasive distant hybridization in nature among cyprinid fishes

BACKGROUND

Genomic data have unveiled a fascinating aspect of the evolutionary past, showing that the mingling of different species through hybridization has left its mark on the histories of numerous life forms. However, the relationship between hybridization events and the origins of cyprinid fishes remains unclear.

RESULTS

In this study, we generated de novo assembled genomes of 8 cyprinid fishes and conducted phylogenetic analyses on 24 species. Widespread allele sharing across species boundaries was observed within 7 subfamilies of cyprinid fishes. Based on a systematic analysis of multiple tissues, we found that the testis exhibited a conserved pattern of divergence between the herbivorous Megalobrama amblycephala and the carnivorous Culter alburnus, suggesting a potential link to incomplete reproductive isolation. Significant differences in the expression of 4 genes (dpp2, ctrl, psb7, and ppce) in the liver and intestine, accompanied by variations in enzyme activities, indicated swift divergence in digestive enzyme secretion. Moreover, we identified introgressed genes linked to organ development in sympatric fishes with analogous feeding habits within the Cultrinae and Leuciscinae subfamilies.

CONCLUSIONS

Our findings highlight the significant role played by incomplete reproductive isolation and frequent gene flow events, particularly those associated with the development of digestive organs, in driving speciation among cyprinid fishes in diverse freshwater ecosystems.

Ecological diversification of sea catfishes is accompanied by genome-wide signatures of positive selection

Habitat transitions have shaped the evolutionary trajectory of many clades. Sea catfishes (Ariidae) have repeatedly undergone ecological transitions, including colonizing freshwaters from marine environments, leading to an adaptive radiation in Australia and New Guinea alongside non-radiating freshwater lineages elsewhere. Here, we generate and analyze one long-read reference genome and 66 short-read whole genome assemblies, in conjunction with genomic data for 54 additional species. We investigate how three major ecological transitions have shaped genomic variation among ariids over their ~ 50 million-year evolutionary history. Our results show that relatively younger freshwater lineages exhibit a higher incidence of positive selection than their more ancient marine counterparts. They also display a larger disparity in body shapes, a trend that correlates with a heightened occurrence of positive selection on genes associated with body size and elongation. Although positive selection in the Australia and New Guinea radiation does not stand out compared to non-radiating lineages overall, selection across the prolactin gene family during the marine-to-freshwater transition suggests that strong osmoregulatory adaptations may have facilitated their colonization and radiation. Our findings underscore the significant role of selection in shaping the genome and organismal traits in response to habitat shifts across macroevolutionary scales.

Transcriptome analysis reveals the importance of phototransduction during the first-feeding in mandarin fish (Siniperca chuatsi)

The mandarin fish (Siniperca chuatsi), as a typical freshwater carnivorous fish, has high economic value. Mandarin fish have a peculiar feeding habit of feeding on other live fry during the first-feeding period, while rejecting zooplankton or particulate feed, which may be attributed to the low expression of zooplankton-associated gene sws1 in mandarin fish. The domesticated strain of mandarin fish could feed on Artemia at 3 days post hatching (dph). However, the mechanism of mandarin fish larvae recognize and forage Artemia as food is still unclear. In this study, we employed transcriptional analysis to identify the representative differential pathways between mandarin fish larvae unfed and fed with Artemia at 3 dph. The comparative transcriptome analysis has unveiled a tapestry of genetic expression, highlighting 403 genes that have been up-regulated and 259 that have been down-regulated, all of which constitute the differentially expressed genes (DEGs). KEGG pathway analysis revealed that the number of differentially expressed genes in the photoconductive signaling pathway was the largest. Next, the Vorinostat (suberoylanilide hydroxamic acid, SAHA) was used to assess whether sws1 induced ingestion of Artemia in mandarin fish larvae. We discovered that SAHA-treated larvae had more food intake of Artemia and up-regulated the transcription level of npy, which might have been associated with the up-regulated of sws1 opsin. Additionally, exposure to 0.5 µM SAHA increased the expression of genes involved in phototransduction pathway. These findings would provide insights on the molecular processes involved in mandarin fish larvae feeding on Artemia at the first-feeding stage.

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

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

Mandarin fish von Hippel-Lindau protein regulates the NF-κB signaling pathway via interaction with IκB to promote fish ranavirus replication

The von Hippel-Lindau tumor suppressor protein (VHL), an E3 ubiquitin ligase, functions as a critical regulator of the oxygen-sensing pathway for targeting hypoxia-inducible factors. Recent evidence suggests that mammalian VHL may also be critical to the NF-κB signaling pathway, although the specific molecular mechanisms remain unclear. Herein, the roles of mandarin fish ( Siniperca chuatsi) VHL ( scVHL) in the NF-κB signaling pathway and mandarin fish ranavirus (MRV) replication were explored. The transcription of scVHL was induced by immune stimulation and MRV infection, indicating a potential role in innate immunity. Dual-luciferase reporter gene assays and reverse transcription quantitative PCR (RT-qPCR) results demonstrated that scVHL evoked and positively regulated the NF-κB signaling pathway. Treatment with NF-κB signaling pathway inhibitors indicated that the role of scVHL may be mediated through scIKKα, scIKKβ, scIκBα, or scp65. Co-immunoprecipitation (Co-IP) analysis identified scIκBα as a novel target protein of scVHL. Moreover, scVHL targeted scIκBα to catalyze the formation of K63-linked polyubiquitin chains to activate the NF-κB signaling pathway. Following MRV infection, NF-κB signaling remained activated, which, in turn, promoted MRV replication. These findings suggest that scVHL not only positively regulates NF-κB but also significantly enhances MRV replication. This study reveals a novel function of scVHL in NF-κB signaling and viral infection in fish.

The Chinese longsnout catfish genome provides novel insights into the feeding preference and corresponding metabolic strategy of carnivores

Fish show variation in feeding habits to adapt to complex environments. However, the genetic basis of feeding preference and the corresponding metabolic strategies that differentiate feeding habits remain elusive. Here, by comparing the whole genome of a typical carnivorous fish (Leiocassis longirostris Günther) with that of herbivorous fish, we identify 250 genes through both positive selection and rapid evolution, including taste receptor taste receptor type 1 member 3 (tas1r3) and trypsin We demonstrate that tas1r3 is required for carnivore preference in tas1r3-deficient zebrafish and in a diet-shifted grass carp model. We confirm that trypsin correlates with the metabolic strategies of fish with distinct feeding habits. Furthermore, marked alterations in trypsin activity and metabolic profiles are accompanied by a transition of feeding preference in tas1r3-deficient zebrafish and diet-shifted grass carp. Our results reveal a conserved adaptation between feeding preference and corresponding metabolic strategies in fish, and provide novel insights into the adaptation of feeding habits over the evolution course.

Effect of different feeding regimens on physiological indicators, intestinal transcriptome, and bacterial flora of mandarin fish (Siniperca chuatsi)

Mandarin fish (Siniperca chuatsi) represents a typical carnivorous freshwater economic fish in China. Recently, the study of their feeding behavior to acclimate formulated diets has become a research focus. This study evaluated the effects of various diets on the body composition, nutritional content, digestive enzyme activity, gene expression, and gut microbiota of mandarin fish. Firstly, no significant differences were found in the muscle's basic nutritional components (moisture, crude protein, crude fat, and crude ash), as well as in the fatty acid and amino acid content, between the live feed group (LFSC) and the compound feed group (CFSC). However, mandarin fish in the LFSC group exhibited significantly higher lipase activity in the liver and intestine compared to the CFSC group, while amylase activity in the intestine showed an opposite pattern. Additionally, intestinal transcriptome analysis revealed 6238 differentially expressed genes and identified several differentially expressed clock genes associated with diet type. Furthermore, gut microbiota analysis indicated that different feeding regimens influenced microbial composition, revealing correlations between bacterial genera and intestinal gene expression levels. These findings provided novel insights into the gut microbiota and transcriptomic responses of mandarin fish to different dietary types.

Multi-genome comparisons reveal gain-and-loss evolution of anti-Mullerian hormone receptor type 2 as a candidate master sex-determining gene in Percidae

BACKGROUND

The Percidae family comprises many fish species of major importance for aquaculture and fisheries. Based on three new chromosome-scale assemblies in Perca fluviatilis, Perca schrenkii, and Sander vitreus along with additional percid fish reference genomes, we provide an evolutionary and comparative genomic analysis of their sex-determination systems.

RESULTS

We explored the fate of a duplicated anti-Mullerian hormone receptor type-2 gene (amhr2bY), previously suggested to be the master sex-determining (MSD) gene in P. flavescens. Phylogenetically related and structurally similar amhr2 duplicates (amhr2b) were found in P. schrenkii and Sander lucioperca, potentially dating this duplication event to their last common ancestor around 19-27 Mya. In P. fluviatilis and S. vitreus, this amhr2b duplicate has been likely lost while it was subject to amplification in S. lucioperca. Analyses of the amhr2b locus in P. schrenkii suggest that this duplication could be also male-specific as it is in P. flavescens. In P. fluviatilis, a relatively small (100 kb) non-recombinant sex-determining region (SDR) was characterized on chromosome 18 using population-genomics approaches. This SDR is characterized by many male-specific single-nucleotide variations (SNVs) and no large duplication/insertion event, suggesting that P. fluviatilis has a male heterogametic sex-determination system (XX/XY), generated by allelic diversification. This SDR contains six annotated genes, including three (c18h1orf198, hsdl1, tbc1d32) with higher expression in the testis than in the ovary.

CONCLUSIONS

Together, our results provide a new example of the highly dynamic sex chromosome turnover in teleosts and provide new genomic resources for Percidae, including sex-genotyping tools for all three known Perca species.

Zebrafish Feed Intake: A Systematic Review for Standardizing Feeding Management in Laboratory Conditions

Zebrafish are one of the most used animal models in biological research and a cost-effective alternative to rodents. Despite this, nutritional requirements and standardized feeding protocols have not yet been established for this species. This is important to avoid nutritional effects on experimental outcomes, and especially when zebrafish models are used in preclinical studies, as many diseases have nutritional confounding factors. A key aspect of zebrafish nutrition is related to feed intake, the amount of feed ingested by each fish daily. With the goal of standardizing feeding protocols among the zebrafish community, this paper systematically reviews the available data from 73 studies on zebrafish feed intake, feeding regimes (levels), and diet composition. Great variability was observed regarding diet composition, especially regarding crude protein (mean 44.98 ± 9.87%) and lipid content (9.91 ± 5.40%). Interestingly, the gross energy levels of the zebrafish diets were similar across the reviewed studies (20.39 ± 2.10 kilojoules/g of feed). In most of the reviewed papers, fish received a predetermined quantity of feed (feed supplied). The authors fed the fish according to the voluntary intake and then calculated feed intake (FI) in only 17 papers. From a quantitative point of view, FI was higher than when a fixed quantity (pre-defined) of feed was supplied. Also, the literature showed that many biotic and abiotic factors may affect zebrafish FI. Finally, based on the FI data gathered from the literature, a new feeding protocol is proposed. In summary, a daily feeding rate of 9-10% of body weight is proposed for larvae, whereas these values are equal to 6-8% for juveniles and 5% for adults when a dry feed with a proper protein and energy content is used.

De Novo Genome Assembly of the Whitespot Parrotfish (Scarus forsteni): A Valuable Scaridae Genomic Resource

Scarus forsteni, a whitespot parrotfish from the Scaridae family, is a herbivorous fish inhabiting coral reef ecosystems. The deterioration of coral reefs has highly affected the habitats of the parrotfish. The decline in genetic diversity of parrotfish emphasizes the critical importance of conserving their genetic variability to ensure the resilience and sustainability of marine ecosystems for future generations. In this study, a genome of S. forsteni was assembled de novo through using Illumina and Nanopore sequencing. The 1.71-Gb genome of S. forsteni, was assembled into 544 contigs (assembly level: contig). It exhibited an N50 length of 17.97 Mb and a GC content percentage of 39.32%. Our BUSCO analysis revealed that the complete protein of the S. forsteni genome had 98.10% integrity. Combined with structure annotation data, 34,140 (74.81%) genes were functionally annotated out of 45,638 predicted protein-coding genes. Upon comparing the genome size and TE content of teleost fishes, a roughly linear relationship was observed between these two parameters. However, TE content is not a decisive factor in determining the genome size of S. forsteni. Population history analysis results indicate that S. forsteni experienced two major population expansions, both of which occurred before the last interglacial period. In addition, through a comparative genomic analysis of the evolutionary relationship of other species, it was found that S. forsteni had the closest relationship with Cheilinus undulatus, another member of the Labridae family. Our expansion and contraction analysis of the gene family showed that the expansion genes were mainly associated with immune diseases, organismal systems, and cellular processes. At the same time, cell transcription and translation, sex hormone regulation, and other related pathways were also more prominent in the positive selection genes. The genomic sequence of S. forsteni offers valuable resources for future investigations on the conservation, evolution, and behavior of fish species.

Regulatory effect of NK homeobox 1 (NKX2.1) on melanocortin 4 receptor (Mc4r) promoter in Mandarin fish

The melanocortin 4 receptor (MC4R) is a G protein-coupled transporter that mediates the regulation of thyroid hormones and leptin on energy balance and food intake. However, the mechanisms of transcriptional regulation of Mc4r by thyroid hormone and leptin in fish have been rarely reported. The messenger RNA expression of Mc4r gene was significantly higher in brain than those in other tissues of mandarin fish. We analyzed the structure and function of a 2029 bp sequence of Mc4r promoter. Meanwhile, overexpression of NKX2.1 and incubation with leptin significantly increased Mc4r promoter activity, but triiodothyronine showed the opposite effect. In addition, mutations in the NKX2.1 binding site abolished not only the activation of Mc4r promoter activity by leptin but also the inhibitory effect of thyroid hormones on Mc4r promoter activity. In summary, these results suggested that thyroid hormones and leptin might regulate the transcriptional expression of Mc4r through NKX2.1.

Multi-genome comparisons reveal gain-and-loss evolution of the anti-Mullerian hormone receptor type 2 gene, an old master sex determining gene, in Percidae

The Percidae family comprises many fish species of major importance for aquaculture and fisheries. Based on three new chromosome-scale assemblies in Perca fluviatilis, Perca schrenkii and Sander vitreus along with additional percid fish reference genomes, we provide an evolutionary and comparative genomic analysis of their sex-determination systems. We explored the fate of a duplicated anti-Mullerian hormone receptor type-2 gene (amhr2bY), previously suggested to be the master sex determining (MSD) gene in P. flavescens. Phylogenetically related and structurally similar amhr2 duplications (amhr2b) were found in P. schrenkii and Sander lucioperca, potentially dating this duplication event to their last common ancestor around 19-27 Mya. In P. fluviatilis and S. vitreus, this amhr2b duplicate has been lost while it was subject to amplification in S. lucioperca. Analyses of the amhr2b locus in P. schrenkii suggest that this duplication could be also male-specific as it is in P. flavescens. In P. fluviatilis, a relatively small (100 kb) non-recombinant sex-determining region (SDR) was characterized on chromosome-18 using population-genomics approaches. This SDR is characterized by many male-specific single-nucleotide variants (SNVs) and no large duplication/insertion event, suggesting that P. fluviatilis has a male heterogametic sex determination system (XX/XY), generated by allelic diversification. This SDR contains six annotated genes, including three (c18h1orf198, hsdl1, tbc1d32) with higher expression in testis than ovary. Together, our results provide a new example of the highly dynamic sex chromosome turnover in teleosts and provide new genomic resources for Percidae, including sex-genotyping tools for all three known Perca species.

Characterization of myosin heavy chain (MYH) genes and their differential expression in white and red muscles of Chinese perch, Siniperca chuatsi

Fish skeletal muscle is composed of two anatomically and functionally different fiber layers, white or fast and red or slow muscles. Myosin, the major structural protein of fish skeletal muscle, contains multiple myosin heavy chain (MYH) isoforms involved in the high plasticity of muscle in response to varying functional demands and/or environmental changes. In this study, we comparatively assayed the cellular and ultrastructural feature of white and red skeletal muscles. Then, a total of 28 class II myosin heavy chain genes were identified in by searching the Chinese perch genome database. Among them, 14 genes code for the fast-muscle-type myosin heavy chain, and 7 genes code for the slow-muscle-type myosin heavy chain. Further, the different isoform gene structures, function domains, phylogenetic relations, and muscle-fiber type-specific expression were characterized. This is the first systematic work on the molecular characterization of class II myosin heavy chain isoforms and the differential analysis of their expression in red and white muscle tissues in Chinese perch Siniperca chuatsi. Our work provided valuable information for a better understanding of myh genes and their molecular characteristics, and the correlations of multiple myosin isoforms with potential functions in response to varying functional demands and/or environmental changes.

Effects and regulatory pathway of proopinmelanocortin on feeding habit domestication in mandarin fish

Proopiomelanocortin (POMC) is a hormone precursor, and has been reported to participate in domestication. However, its effects on feeding habit domestication in fish are poorly understood. Mandarin fish (Siniperca chuatsi) feeds solely on live prey fish since first-feeding. In the present study, the high expression of pomc in mandarin fish, both the pomc siRNA and MC4R inhibitor treatments increased the success rate of domestication from live prey fish to dead prey fish and food intake of dead prey fish, suggesting the role of pomc on the special feeding habit of live prey fish in mandarin fish. In addition, one c-fos binding site was identified in the region that from -1053 bp to -931 bp upstream of the transcription start site of pomc, and this region exhibited positive promoter activity. The mandarin fish brain cells treated with c-fos siRNA displayed suppressed pomc mRNA expression, indicating that c-fos positively regulated pomc expression. Furthermore, the mRNA expression of c-fos was higher in the mandarin fish which were more difficult to domesticate. The results of ChIP assay and inhibitor treatment confirmed that the activation of c-fos gene by histone H3K4me3 was catalyzed by Setd1b in mandarin fish. Three open peaks were found at the upstream regulatory region of setd1b by ATAC-seq, and the mRNA expression of setd1b was higher in the mandarin fish which were more difficult to domesticate. These results indicated that Setd1b could methylate histone H3K4 to activate the c-fos transcription, maintaining the high expression of pomc, which might contribute to the special feeding habit of mandarin fish.

Molecular Mechanisms Underlying Vertebrate Adaptive Evolution: A Systematic Review

Adaptive evolution is a process in which variation that confers an evolutionary advantage in a specific environmental context arises and is propagated through a population. When investigating this process, researchers have mainly focused on describing advantageous phenotypes or putative advantageous genotypes. A recent increase in molecular data accessibility and technological advances has allowed researchers to go beyond description and to make inferences about the mechanisms underlying adaptive evolution. In this systematic review, we discuss articles from 2016 to 2022 that investigated or reviewed the molecular mechanisms underlying adaptive evolution in vertebrates in response to environmental variation. Regulatory elements within the genome and regulatory proteins involved in either gene expression or cellular pathways have been shown to play key roles in adaptive evolution in response to most of the discussed environmental factors. Gene losses were suggested to be associated with an adaptive response in some contexts. Future adaptive evolution research could benefit from more investigations focused on noncoding regions of the genome, gene regulation mechanisms, and gene losses potentially yielding advantageous phenotypes. Investigating how novel advantageous genotypes are conserved could also contribute to our knowledge of adaptive evolution.

Whole-genome resequencing and bisulfite sequencing provide new insights into the feeding habit domestication in mandarin fish (Siniperca chuatsi)

Mandarin fish (Siniperca chuatsi) is one of the most economically important fish in China. However, it has the peculiar feeding habit that it feeds solely on live prey fish since first-feeding, while refuses dead prey fish or artificial diets. After the specific training procedure, partial individuals could accept dead prey fish and artificial diets. The genetic basis of individual difference in artificial diet feeding habit is still unknown. In the present study, the resequencing was performed between 10 individuals which could be domesticated to accept artificial diets and 10 individuals which could not. Through the selective sweep analysis based on heterozygosity (Hp) and population differentiation coefficient (Fst), 57 candidate windows were identified as the putative selected regions for feeding habit domestication of mandarin fish, involved in 149 genes. These genes were related to memory, vision and olfaction function, which could be potential targets of molecular marker assistant breeding of artificial diet feeding trait. Beside of the DNA sequence, we also explored the potential role of DNA methylation in feeding habit domestication in mandarin fish. Whole-genome bisulfite sequencing was performed between the individuals which could be domesticated to accept artificial diets and those could not. 5,976 differentially methylated regions were identified, referring to 3,522 genes, such as the genes involved in cAMP signaling pathway. The DNA methylation changes of these genes might contribute to the adaption of artificial diets in mandarin fish. In conclusion, the putative selected regions and the differentially methylated regions were identified in the whole genome, providing new insights into the feeding habit domestication from live prey fish to artificial diets in mandarin fish. And the involved genes were identified as the candidate genes for molecular breeding of artificial diet utilization in mandarin fish.

Comparative Study of the Molecular Characterization, Evolution, and Structure Modeling of Digestive Lipase Genes Reveals the Different Evolutionary Selection Between Mammals and Fishes

Vertebrates need suitable lipases to digest lipids for the requirement of energy and essential nutrients; however, the main digestive lipase genes of fishes have certain controversies. In this study, two types of digestive lipase genes (pancreatic lipase (pl) and bile salt-activated lipase (bsal)) were identified in mammals and fishes. The neighborhood genes and key active sites of the two lipase genes were conserved in mammals and fishes. Three copies of PL genes were found in mammals, but only one copy of the pl gene was found in most of the fish species, and the pl gene was even completely absent in some fish species (e.g., zebrafish, medaka, and common carp). Additionally, the hydrophobic amino acid residues (Ile and Leu) which are important to pancreatic lipase activity were also absent in most of the fish species. The PL was the main digestive lipase gene in mammals, but the pl gene seemed not to be the main digestive lipase gene in fish due to the absence of the pl gene sequence and the important amino acid residues. In contrast, the bsal gene existed in all fish species, even two to five copies of bsal genes were found in most of the fishes, but only one copy of the BSAL gene was found in mammals. The amino acid residues of bile salt-binding sites and the three-dimensional (3D) structure modeling of Bsal proteins were conserved in most of the fish species, so bsal might be the main digestive lipase gene in fish. The phylogenetic analysis also indicated that pl or bsal showed an independent evolution between mammals and fishes. Therefore, we inferred that the evolutionary selection of the main digestive lipase genes diverged into two types between mammals and fishes. These findings will provide valuable evidence for the study of lipid digestion in fish.

Genomic architecture of adaptive radiation and hybridization in Alpine whitefish

Adaptive radiations represent some of the most remarkable explosions of diversification across the tree of life. However, the constraints to rapid diversification and how they are sometimes overcome, particularly the relative roles of genetic architecture and hybridization, remain unclear. Here, we address these questions in the Alpine whitefish radiation, using a whole-genome dataset that includes multiple individuals of each of the 22 species belonging to six ecologically distinct ecomorph classes across several lake-systems. We reveal that repeated ecological and morphological diversification along a common environmental axis is associated with both genome-wide allele frequency shifts and a specific, larger effect, locus, associated with the gene edar. Additionally, we highlight the possible role of introgression between species from different lake-systems in facilitating the evolution and persistence of species with unique trait combinations and ecology. These results highlight the importance of both genome architecture and secondary contact with hybridization in fuelling adaptive radiation.

Screening of Genes Related to Sex Determination and Differentiation in Mandarin Fish (Siniperca chuatsi)

Mandarin fish has an XX/XY sex-determination system. The female mandarin fish is typically larger than the male. Sex identification and the discovery of genes related to sex determination in mandarin fish have important theoretical significance in the elucidation of the regulation and evolutionary mechanism of animal reproductive development. In this study, the chromosome-level genome of a female mandarin fish was assembled, and we found that LG24 of the genome was an X chromosome. A total of 61 genes on the X chromosome showed sex-biased expression. Only six gonadal genes (LG24G00426, LG24G003280, LG24G003300, LG24G003730, LG24G004200, and LG24G004770) were expressed in the testes, and the expression of the other gene LG24G003870 isoform 1 in the ovaries was significantly higher than that in the testes (p < 0.01). Five (except LG24G003280 and LG24G003300) of the seven aforementioned genes were expressed at the embryonic development stage, suggesting their involvement in early sex determination. The expression of LG24G004770 (encoding HS6ST 3-B-like) was also significantly higher in female muscles than in male muscles (p < 0.01), indicating other functions related to female growth. ZP3 encoded by LG24G003870 isoform 1 increased the C-terminal transmembrane domain, compared with that encoded by other fish zp3 isoforms, indicating their different functions in sex determination or differentiation. This study provides a foundation for the identification of sex-determining genes in mandarin fish.

Genome-wide identification and expression patterns of opsin genes during larval development in Chinese perch (Siniperca chuatsi)

Vision is important for fish to forage food and fishes express opsin genes to receive visual signals. Chinese perch (Siniperca chuatsi) larvae prey on other fish species larvae at firstfeeding but donoteat any zooplankton, the expression of opsin genes in S. chuatsilarvae is unknown. In this study, we conducted a whole-genome analysis and demonstrated that S. chuatsihave5cone opsin genes (sws1, sws2Aα, sws2Aβ, rh2and lws)and 2 rod opsin genes (rh1and rh1-exorh). The syntenicanalysisshowedthe flanking genes ofall opsin genes were conserved during fish evolution, but the ancestorof S. chuatsimightlost some opsin gene copies duringtheevolution.The phylogeneticanalysisshowed sws1of S. chuatsiwas closest to those of Lates calcariferwhich had a truncated sws1gene; the sws2Aα, sws2Aβ,lws,rh2,rh1 andrh1-exorh of S. chuatsihad a closer relationship with those of Percomorpha fishes.Importantly, results of in situhybridization showed the sws1 opsingene,which is related to forage zooplankton,had extremely low levelexpression in retinaat early stages.Surprisingly, the rh2 opsin gene had a high level expression at firstfeeding stage. The sws2Aα, sws2Aβand lwshad a little expression at early stages but the lwsshowed a increasing trend with larval development, rh1 opsin gene expression appeared at15 dph. In thisstudy, we found a specialpattern of visual opsin genes expression in S. chuatsi, it might influence the larval first feeding and feeding habit.

Improved contiguity of the threespine stickleback genome using long-read sequencing

While the cost and time for assembling a genome has drastically decreased, it still remains a challenge to assemble a highly contiguous genome. These challenges are rapidly being overcome by the integration of long-read sequencing technologies. Here, we use long-read sequencing to improve the contiguity of the threespine stickleback fish (Gasterosteus aculeatus) genome, a prominent genetic model species. Using Pacific Biosciences sequencing, we assembled a highly contiguous genome of a freshwater fish from Paxton Lake. Using contigs from this genome, we were able to fill over 76.7% of the gaps in the existing reference genome assembly, improving contiguity over fivefold. Our gap filling approach was highly accurate, validated by 10X Genomics long-distance linked-reads. In addition to closing a majority of gaps, we were able to assemble segments of telomeres and centromeres throughout the genome. This highlights the power of using long sequencing reads to assemble highly repetitive and difficult to assemble regions of genomes. This latest genome build has been released through a newly designed community genome browser that aims to consolidate the growing number of genomics datasets available for the threespine stickleback fish.

Comparative transcriptomic analysis reveals an association of gibel carp fatty liver with ferroptosis pathway

BACKGROUND

Fatty liver has become a main problem that causes huge economic losses in many aquaculture modes. It is a common physiological or pathological phenomenon in aquaculture, but the causes and occurring mechanism are remaining enigmatic.

METHODS

Each three liver samples from the control group of allogynogenetic gibel carp with normal liver and the overfeeding group with fatty liver were collected randomly for the detailed comparison of histological structure, lipid accumulation, transcriptomic profile, latent pathway identification analysis (LPIA), marker gene expression, and hepatocyte mitochondria analyses.

RESULTS

Compared to normal liver, larger hepatocytes and more lipid accumulation were observed in fatty liver. Transcriptomic analysis between fatty liver and normal liver showed a totally different transcriptional trajectory. GO terms and KEGG pathways analyses revealed several enriched pathways in fatty liver, such as lipid biosynthesis, degradation accumulation, peroxidation, or metabolism and redox balance activities. LPIA identified an activated ferroptosis pathway in the fatty liver. qPCR analysis confirmed that gpx4, a negative regulator of ferroptosis, was significantly downregulated while the other three positively regulated marker genes, such as acsl4, tfr1 and gcl, were upregulated in fatty liver. Moreover, the hepatocytes of fatty liver had more condensed mitochondria and some of their outer membranes were almost ruptured.

CONCLUSIONS

We reveal an association between ferroptosis and fish fatty liver for the first time, suggesting that ferroptosis might be activated in liver fatty. Therefore, the current study provides a clue for future studies on fish fatty liver problems.