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

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.

In this genomic study on Alpine whitefish radiations, the authors reveal details on the genetic architecture underlying the repeated eco-morphological diversification and the role of hybridization in the evolution of endemic whitefish species.

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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07621-2.