The evolutionary origin and domestication history of goldfish (Carassius auratus)

Chronic hypoxia induces alternative splicing of transcripts in the goldfish brain

Several species evolved mechanisms to tolerate periods of severe environmental hypoxia and anoxia. Among them, goldfish are unique as they do not enter a comatose state under such conditions. Taking advantage of the recently published and annotated goldfish genome, we had previously profiled the transcriptomic response of the goldfish brain under normoxic (21 kPa oxygen saturation, N) and hypoxic conditions (2.1 kPa oxygen saturation) after 1 and 4 weeks (1WH, 4WH). Using the RNA-Seq data, we report the occurrence of alternative mRNA splicing (skipped exon, retained intron, alternative 3' or 5' splice sites, and mutually exclusive exons). At 1WH/N, 1004 significant alternative splicing events on 769 gene loci were identified, increasing to 1187 on 963 loci at 4WH/N. There were 305 loci with alternatively spliced transcripts common to both 1WH/N and 4WH/N, 221 of which exhibited the same precise location and splicing mechanism. Specific gene transcripts affected by alternative splicing events were almost entirely different from previously identified differentially expressed genes under chronic hypoxia. GO-term enrichment analyses of gene loci of alternatively spliced transcripts, however, did include similar pathways as previously identified for DEGs. These include epigenetic machinery, ion channel activity (1WH/N), glutamate signaling (4WH/N), endothelial cell function, and ATP hydrolyzation pathways (1WH/N + 4WH/N). We describe selected examples of alternatively spliced transcripts to discuss possible functional relevance in the goldfish brain response to chronic hypoxia. Together, our data identified an additional layer of regulation in brain pathways relevant to hypoxia tolerance in goldfish, which complement previously reported gene expression changes.

Enhancing genome editing efficiency in goldfish (Carassius auratus) through utilization of CRISPR-Cas12a (Cpf1) temperature dependency

The CRISPR/Cas technology has demonstrated revolutionary potential across various fields, including agriculture, medicine, and food safety detection. However, the utility of CRISPR/Cas12a, a particularly promising gene-editing tool, is constrained by its temperature sensitivity, limiting its application in low-temperature environments. In this study, we developed a gene-editing technique based on the CRISPR/Cas12a system in the poikilothermic species goldfish Carassius auratus. We systematically evaluated the editing efficiencies of LbCas12a and AsCas12a on the tyrosinase (tyr) gene under varying temperature conditions. Our results revealed a pronounced temperature dependence of Cas12a, with elevated temperatures markedly enhancing its editing activity, particularly for AsCas12a. A brief one-hour high-temperature treatment was sufficient to achieve effective gene disruption, underscoring CRISPR/Cas12a as a rapid and efficient gene-editing tool. Temperature was utilized as a conditional trigger for Cas12a-mediated gene knockout, enabling precise modulation of gene disruption at specific embryonic developmental stages. Whole-genome resequencing of the mutants confirmed the absence of off-target effects, further emphasizing the precision of this editing process. These findings indicated that CRISPR/Cas12a represented a viable alternative to the widely utilized CRISPR/Cas9 system and could be applied in conjunction, thereby expanding the potential applications of gene-editing technologies.

Fish models to explore epigenetic determinants of hypoxia-tolerance

The occurrence of environmental hypoxia in freshwater and marine aquatic systems has increased over the last century and is predicted to further increase with climate change. As members of the largest extant vertebrate group, freshwater fishes, and to a much lesser extent marine fishes, are vulnerable to increased occurrence of hypoxia. This is important as fishes render important ecosystem services and have important cultural and economic roles. Evolutionarily successful, fishes have adapted to diverse aquatic freshwater and marine habitats with different oxygen conditions. While some fishes exhibit genetic adaptions to tolerate hypoxia and even anoxia, others are limited to oxygen-rich habitats. Recent advances in molecular epigenetics have shown that some epigenetic machinery, especially histone- and DNA demethylases, is directly dependent on oxygen and modulates important transcription-regulating epigenetic marks in the process. At the post-transcriptional level, hypoxia has been shown to affect non-coding microRNA abundance. Together, this evidence adds a new molecular epigenetic basis to study hypoxia tolerance in fishes. Here, we review the documented and predicted changes in environmental hypoxia in aquatic systems and discuss the diversity and comparative physiology of hypoxia tolerance in fishes, including molecular and physiological adaptations. We then discuss how recent mechanistic advances in environmental epigenetics can inform future work probing the role of oxygen-dependent epigenetic marks in shaping organismal hypoxia-tolerance in fishes with a focus on within- and between-species variation, acclimation, inter- and multigenerational plasticity, and multiple climate-change stressors. We conclude by describing the translational potential of this approach for conservation physiology, ecotoxicology, and aquaculture.

Variation and Interaction of Distinct Subgenomes Contribute to Growth Diversity in Intergeneric Hybrid Fish

Intergeneric hybridization greatly reshapes regulatory interactions among allelic and non-allelic genes. However, their effects on growth diversity remain poorly understood in animals. In this study, we conducted whole-genome sequencing and RNA sequencing analyses in diverse hybrid varieties resulting from the intergeneric hybridization of goldfish (Carassius auratus red var.) and common carp (Cyprinus carpio). These hybrid individuals were characterized by distinct mitochondrial genomes and copy number variations. Through a weighted gene correlation network analysis, we identified 3693 genes as candidate growth-regulating genes. Among them, the expression of 3672 genes in subgenome R (originating from goldfish) displayed negative correlations with body weight, whereas 20 genes in subgenome C (originating from common carp) exhibited positive correlations. Notably, we observed intriguing expression patterns of solute carrier family 2 member 12 (slc2a12) in subgenome C, showing opposite correlations with body weight that changed with water temperatures, suggesting differential interactions between feeding activity and weight gain in response to seasonal changes for hybrid animals. In 40.30% of alleles, we observed dominant trans-regulatory effects in the regulatory interactions between distinct alleles from subgenomes R and C. Integrating analyses of allele-specific expression and DNA methylation data revealed that DNA methylation on both subgenomes shaped the relative contribution of allelic expression to the growth rate. These findings provide novel insights into the interactions of distinct subgenomes that underlie heterosis in growth traits and contribute to a better understanding of multiple allelic traits in animals.

Identification and functional analysis of six melanocortin-4-receptor-like (MC4R-like) mutations in goldfish (Carassius auratus)

Melanocortin receptor-4 (MC4R) belongs to the G protein-coupled receptor family, characterized by a classical structure of seven transmembrane domains (7TMD). They play an important role in food intake and weight regulation. In the present study, we identified melanocortin-4-receptor-like (caMC4RL) mutants of goldfish from the Qian River in the Qin Ling region and characterized their functional properties, including the constitutive activities of the mutants, ligand-induced cAMP and ERK1/2 accumulation, and AMPK activation. The results show that six caMC4RL mutants were identified in goldfish from the Qian River in the Qin Ling region, and are located in the conserved position of the Cyprinidae MC4Rs. The mutations (E57K, P296S, and R302T/K) result in the loss of Gs signaling function. The mutations (P296 and R302T/K) exhibited biased signaling in response to ACTH stimulation in the MAPK/ERK pathway. In addition, the E57K mutant may play a role in weight regulation and could serve as molecular markers for molecular breeding. These data will provide fundamental information for functional studies of teleost GPCR mutants and MC4R isoforms.

Goldfish phenomics reveals commonalities and a lack of universality in the domestication process for ornamentation

Domestication process effects are manifold, affecting genotype and phenotype, and assumed to be universal in animals by part of the scientific community. While mammals and birds have been thoroughly investigated, from taming to intensive selective breeding, fish domestication remains comparatively unstudied. The most widely bred and traded ornamental fish species worldwide, the goldfish, underwent the effect of long-term artificial selection on differing skeletal and soft tissue modules through ornamental domestication. Here, we provide a global morphological analysis in this emblematic ornamental domesticated fish. We demonstrate that goldfish exhibit unique morphological innovations in whole-body, cranial, and sensory (Weberian ossicles and brain) anatomy compared to their evolutionary clade, highlighting a remarkable morphological disparity within a single species comparable to that of a macroevolutionary radiation. In goldfish, as in the case of dogs and pigeons in their respective evolutionary contexts, the most ornamented varieties are extremes in the occupied morphological space, emphasizing the power of artificial selection for nonadaptive traits. Using 21st century tools on a dataset comprising the 16 main goldfish breeds, 23 wild close relatives, and 39 cypriniform species, we show that Charles Darwin's expressed wonder at the goldfish is justified. There is a commonality of overall pattern in the morphological differentiation of domesticated forms selected for ornamental purposes, but the singularity of goldfish occupation and extension within (phylo)morphospaces, speaks against a universality in the domestication process.

Chromosome-level genome assembly of the mud carp (Cirrhinus molitorella) using PacBio HiFi and Hi-C sequencing

The mud carp (Cirrhinus molitorella) is an important economic farmed fish, mainly distributed in South China and Southeast Asia due to its strong adaptability and high yield. Despite its economic importance, the paucity of genomic information has constrained detailed genetic research and breeding efforts. In this study, we utilized PacBio HiFi long-read sequencing and Hi-C technologies to generate a meticulously assembled chromosome-level genome of the mud carp. This assembly spans 1,033.41 Mb, with an impressive 99.82% distributed across 25 chromosomes. The contig N50 and scaffold N50 are 33.29 Mb and 39.86 Mb, respectively. The completeness of the mud carp genome assembly is highlighted by a BUSCO score of 98.05%. We predict 25,865 protein-coding genes, with a BUSCO score of 96.54%, and functional annotations for 91.83% of these genes. Approximately 52.21% of the genome consists of repeat elements. This high-fidelity genome assembly is a vital resource for advancing molecular breeding, comparative genomics, and evolutionary studies of the mud carp and related species.

A novel goldfish orthotopic xenograft model of hepatocellular carcinoma developed to evaluate antitumor drug efficacy

Tumor xenograft animal models play a crucial role in hepatocellular carcinoma (HCC) research. Mice xenograft models are time consuming, laborious and expensive while zebrafish tumor xenograft models are cost-effective and effortless. However, the development of orthotopic xenograft models for HCC in zebrafish embryos has been challenging due to the small size of zebrafish livers. In this study, we utilized 7-day-old goldfish embryos as hosts and successfully established an orthotopic xenograft model of HCC in goldfish livers. Through injecting fluorescence labeled HCC cells into the liver of goldfish, we could visualize the proliferation and migration of tumor cells in vivo. In addition, we found that the temperature of 36 °C was better for tumor cell survival in goldfish larvae compared to 28 °C, assessed by EdU and TUNEL assays. Moreover, macrophage infiltration in the goldfish liver could be evaluated by neutral red staining. Finally, we evaluated the efficacy of the targeted therapy drug Sorafenib and the traditional Chinese medicine, Huaier granules, alone or in combination in the goldfish HCC orthotopic xenograft model. We found that the combination therapy showed the best efficacy against HCC cells in terms of macrophage infiltration, polarization as well as tumor cells proliferation, metastasis and apoptosis. In conclusion, the proposed goldfish HCC orthotopic xenograft model opens new avenues for HCC related research, including evaluation of tumor progression, cell interactions in the immune microenvironment, drug efficacy, and screening of anti-tumor drugs.

The Identification of Proteomic Signatures Associated with Alkaline Tolerance in the Skin Mucus of Crucian Carp (Carassius auratus)

The skin is covered by a protective mucus layer, which is essential to the innate defense mechanism of fish. Investigating the response of skin mucus to various toxic stresses is crucial for enhancing its ability to tackle environmental challenges and developing strategies to mitigate toxic effects. Alkalinity stress assays (50 mmol/L NaHCO3) were conducted on crucian carp (Carassius auratus) from Lake Dali Nur (pH = 9.6) and Ping Xiang red crucian carp from freshwater (pH = 7) over 7 days. The expression of skin mucous proteins was analyzed using the liquid chromatography (LC)-spectrometry (MS)/MS Analysis-Data-independent acquisition (DIA) mode. A total of 12,537 proteins were identified across 20 samples from four groups, with 12,025 quantified. In the alkaline water population, high alkali stress resulted in the up-regulation of 139 proteins and the down-regulation of 500 proteins. In contrast, the freshwater population showed an increase in 112 proteins and a decrease in 120; both populations had a total of 23 genes up-regulated and 21 down-regulated. The protein regulatory network for the alkaline water group included 3146 pairwise interactions among 464 nodes, with only 20 being differentially expressed proteins. Conversely, the freshwater group's network comprised just 1027 specific interactions across 337 nodes, with 6 corresponding to differentially expressed proteins. A common protein regulatory network responding to high alkali stress was extracted and visualized for both populations. Based on their regulatory relationships and expression levels, these proteins are hypothesized to play similar roles under high alkali stress. Notably, the alpha-globin fragment and keratin type I cytoskeletal 13-like proteins showed markedly up-regulated expression, with the alpha-globin fragment increasing nearly a thousandfold from an extremely low level. This suggests it could serve as a potential biomarker for alkali tolerance, warranting further investigation.

Genetic improvement and genomic resources of important cyprinid species: status and future perspectives for sustainable production

Cyprinid species are the most cultured aquatic species around the world in terms of quantity and total value. They account for 25% of global aquaculture production and significantly contribute to fulfilling the demand for fish food. The aquaculture of these species is facing severe concerns in terms of seed quality, rising feed costs, disease outbreaks, introgression of exotic species, environmental impacts, and anthropogenic activities. Numerous researchers have explored biological issues and potential methods to enhance cyprinid aquaculture. Selective breeding is extensively employed in cyprinid species to enhance specific traits like growth and disease resistance. In this context, we have discussed the efforts made to improve important cyprinid aquaculture practices through genetic and genomic approaches. The recent advances in DNA sequencing technologies and genomic tools have revolutionized the understanding of biological research. The generation of a complete genome and other genomic resources in cyprinid species has significantly strengthened molecular-level investigations into disease resistance, growth, reproduction, and adaptation to changing environments. We conducted a comprehensive review of genomic research in important cyprinid species, encompassing genome, transcriptome, proteome, metagenome, epigenome, etc. This review reveals that considerable data has been generated for cyprinid species. However, the seamless integration of this valuable data into genetic selection programs has yet to be achieved. In the upcoming years, genomic techniques, gene transfer, genome editing tools are expected to bring a paradigm shift in sustainable cyprinid aquaculture production. The comprehensive information presented here will offer insights for the cyprinid aquaculture research community.

Allopolyploid subgenome identification and implications for evolutionary analysis

Whole-genome duplications (WGDs) are widespread genomic events in eukaryotes that are hypothesized to contribute to the evolutionary success of many lineages, including flowering plants, Saccharomyces yeast, and vertebrates. WGDs generally can be classified into autopolyploids (ploidy increase descended from one species) or allopolyploids (ploidy increase descended from multiple species). Assignment of allopolyploid progenitor species (called subgenomes in the polyploid) is important to understanding the biology and evolution of polyploids, including the asymmetric subgenome evolution following hybridization (biased fractionation). Here, I review the different methodologies used to identify the ancestors of allopolyploid subgenomes, discuss the advantages and disadvantages of these methods, and outline the implications of how these methods affect the subsequent evolutionary analysis of these genomes.

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.

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

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

Exploring the origin of a unique mutant allele in twin-tail goldfish using CRISPR/Cas9 mutants

Artificial selection has been widely applied to genetically fix rare phenotypic features in ornamental domesticated animals. For many of these animals, the mutated loci and alleles underlying rare phenotypes are known. However, few studies have explored whether these rare genetic mutations might have been fixed due to competition among related mutated alleles or if the fixation occurred due to contingent stochastic events. Here, we performed genetic crossing with twin-tail ornamental goldfish and CRISPR/Cas9-mutated goldfish to investigate why only a single mutated allele-chdS with a E127X stop codon (also called chdAE127X)-gives rise to the twin-tail phenotype in the modern domesticated goldfish population. Two closely related chdS mutants were generated with CRISPR/Cas9 and compared with the E127X allele in F2 and F3 generations. Both of the CRISPR/Cas9-generated alleles were equivalent to the E127X allele in terms of penetrance/expressivity of the twin-tail phenotype and viability of carriers. These findings indicate that multiple truncating mutations could have produced viable twin-tail goldfish. Therefore, the absence of polymorphic alleles for the twin-tail phenotype in modern goldfish likely stems from stochastic elimination or a lack of competing alleles in the common ancestor. Our study is the first experimental comparison of a singular domestication-derived allele with CRISPR/Cas9-generated alleles to understand how genetic fixation of a unique genotype and phenotype may have occurred. Thus, our work may provide a conceptual framework for future investigations of rare evolutionary events in domesticated animals.

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.

Population structure and adaptability analysis of Schizothorax o'connori based on whole-genome resequencing

BACKGROUND

Schizothorax o'connori is an endemic fish distributed in the upper and lower reaches of the Yarlung Zangbo River in China. It has experienced a fourth round of whole gene replication events and is a good model for exploring the genetic differentiation and environmental adaptability of fish in the Qinghai-Tibet Plateau. The uplift of the Qinghai-Tibet Plateau has led to changes in the river system, thereby affecting gene exchange and population differentiation between fish populations. With the release of fish whole genome data, whole genome resequencing has been widely used in genetic evolutionary analysis and screening of selected genes in fish, which can better elucidate the genetic basis and molecular environmental adaptation mechanisms of fish. Therefore, our purpose of this study was to understand the population structure and adaptive characteristics of S. o'connori using the whole-genome resequencing method.

RESULTS

The results showed that 23,602,746 SNPs were identified from seven populations, mostly distributed on chromosomes 2 and 23. There was no significant genetic differentiation between the populations, and the genetic diversity was relatively low. However, the Zangga population could be separated from the Bomi, Linzhi, and Milin populations in the cluster analysis. Based on historical dynamics analysis of the population, the size of the ancestral population of S. o'connori was affected by the late accelerated uplift of the Qinghai Tibet Plateau and the Fourth Glacial Age. The selected sites were mostly enriched in pathways related to DNA repair and energy metabolism.

CONCLUSION

Overall, the whole-genome resequencing analysis provides valuable insights into the population structure and adaptive characteristics of S. o'connori. There was no obvious genetic differentiation at the genome level between the S. o'connori populations upstream and downstream of the Yarlung Zangbo River. The current distribution pattern and genetic diversity are influenced by the late accelerated uplift of the Qinghai Tibet Plateau and the Fourth Ice Age. The selected sites of S. o'connori are enriched in the energy metabolism and DNA repair pathways to adapt to the low temperature and strong ultraviolet radiation environment at high altitude.

Effects of Water Temperature on the Body Color and Expression of the Genes Related to Body Color Regulation in the Goldfish

Melanin-concentrating hormone (MCH), melanocyte-stimulating hormone (MSH), and somatolactin (SL) in the hypothalamus-pituitary axis are associated with body color regulation in teleost fish. Although these hormones' production and secretion respond well to light environments, such as background color, little is known about the effects of different water temperatures. We investigated the effects of water temperature, 10°C, 20°C, and 30°C, on body color and the expression of these genes and corresponding receptor genes in goldfish. The body color in white background (WBG) becomes paler at the higher water temperature, although no difference was observed in black background (BBG). Brain mRNA contents of proMCH genes (pmch1 and pmch2) increased at 30°C and 20°C compared to 10°C in WBG, respectively. Apparent effects of background color and temperature on the pituitary mRNA contents of a POMC gene (pomc) were not observed. The pituitary mRNA contents of the SLα gene were almost double those on a WBG at any temperature, while those of the SLβ gene (slb) at 30°C tended to be higher than those at 10°C and 20°C on WBG and BBG. The scale mRNA contents of the MCH receptor gene (mchr2) in WBG were higher than those in BBG at 30°C. The highest scale mRNA contents of MSH receptor (mc1r and mc5r) on BBG were observed at 20°C, while the lowest respective mRNA levels were observed at 30°C on WBG. These results highlight the importance of temperature for the endocrinological regulation of body color, and darker background color may stabilize those endocrine functions.

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.

Genetic Diversity and Signatures of Selection in the Roughskin Sculpin (Trachidermus fasciatus) Revealed by Whole Genome Sequencing

The roughskin sculpin (Trachidermus fasciatus) is an endangered fish species in China. In recent years, artificial breeding technology has made significant progress, and the population of roughskin sculpin has recovered in the natural environment through enhancement programs and the release of juveniles. However, the effects of released roughskin sculpin on the genetic structure and diversity of wild populations remain unclear. Studies on genetic diversity analysis based on different types and numbers of molecular markers have yielded inconsistent results. In this study, we obtained 2,610,157 high-quality SNPs and 494,698 InDels through whole-genome resequencing of two farmed populations and one wild population. Both farmed populations showed consistent levels of genomic polymorphism and a slight increase in linkage compared with wild populations. The population structure of the two farmed populations was distinct from that of the wild population, but the degree of genetic differentiation was low (overall average Fst = 0.015). Selective sweep analysis showed that 523,529 genes were selected in the two farmed populations, and KEGG enrichment analysis showed that the selected genes were related to amino acid metabolism, which might be caused by artificial feeding. The findings of this study provide valuable additions to the existing genomic resources to help conserve roughskin sculpin populations.

In Vivo Imaging Sheds Light on the Susceptibility and Permissivity of Carassius auratus to Cyprinid Herpesvirus 2 According to Developmental Stage

Cyprinid herpesvirus 2 (CyHV-2) is a virus that causes mass mortality in economically important Carassius spp. However, there have been no comprehensive studies into host susceptibility or permissivity with respect to developmental stage, and the major portal of viral entry into the host is still unclear. To help bridge these knowledge gaps, we developed the first ever recombinant strain of CyHV-2 expressing bioluminescent and fluorescent reporter genes. Infection of Carassius auratus hosts with this recombinant by immersion facilitated the exploitation of various in vivo imaging techniques to establish the spatiotemporal aspects of CyHV-2 replication at larval, juvenile, and adult developmental stages. While less susceptible than later developmental stages, larvae were most permissive to CyHV-2 replication, leading to rapid systemic infection and high mortality. Permissivity to CyHV-2 decreased with advancing development, with adults being the least permissive and, thus, also exhibiting the least mortality. Across all developmental stages, the skin was the most susceptible and permissive organ to infection at the earliest sampling points post-infection, indicating that it represents the major portal of entry into these hosts. Collectively these findings provide important fundamental insights into CyHV-2 pathogenesis and epidemiology in Carassius auratus with high relevance to other related economically important virus-host models.

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.

Nigericin treatment activates endoplasmic reticulum apoptosis pathway in goldfish kidney leukocytes

Nigericin has been reported to induce apoptosis and pyroptosis in mammalian models. However, the effects and mechanism underlying the immune responses of teleost HKLs induced by nigericin remain enigmatic. To decipher the mechanism after nigericin treatment, the transcriptomic profile of goldfish HKLs was analyzed. The results demonstrated that a total of 465 differently expressed genes (DEGs) with 275 up-regulated and 190 down-regulated genes were identified between the control and nigericin treated groups. Among them, the top 20 DEG KEGG enrichment pathways were observed including apoptosis pathways. In addition, the expression level of selected genes (ADP4, ADP5, IRE1, MARCC, ALR1, DDX58) by quantitative real-time PCR showed a significant change after treatment with nigericin, which was generally identical to the expression patterns of the transcriptomic data. Furthermore, the treatment could induce cell death of HKLs, which was confirmed by LDH release and annexin V-FITC/PI assays. Taken together, our results support the idea that nigericin treatment might activate the IRE1-JNK apoptosis pathway in goldfish HKLs, which will provide insights into the mechanisms underlying HKLs immunity towards apoptosis or pyroptosis regulation in teleosts.

Efficiently Editing Multiple Duplicated Homeologs and Alleles for Recurrent Polyploids

Research on the evolutionary fate of duplicated genes in recurrent polyploids is scarce due to the difficulties in disentangling the different homeologs and alleles of duplicated genes. This chapter describes the detailed procedures to identify different homeologs and alleles of duplicated genes, to analyze their molecular characteristics, and to reveal their functional divergence by gene editing with CRISPR/Cas9 (clustered regularly interspaced short palindromic repeats/CRISPR-associated system 9). Using the gene editing approach, we efficiently constructed multiple knockout mutant lines with single or simultaneously disrupted different homeologs or alleles in a recurrent polyploid fish, demonstrating its usability for targeting and mutating multiple divergent homeologs and alleles in recurrent duplicated genomes.

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.

Sequencing and Assembly of Polyploid Genomes

Polyploidy has been observed throughout major eukaryotic clades and has played a vital role in the evolution of angiosperms. Recent polyploidizations often result in highly complex genome structures, posing challenges to genome assembly and phasing. Recent advances in sequencing technologies and genome assembly algorithms have enabled high-quality, near-complete chromosome-level assemblies of polyploid genomes. Advances in novel sequencing technologies include highly accurate single-molecule sequencing with HiFi reads, chromosome conformation capture with Hi-C technique, and linked reads sequencing. Additionally, new computational approaches have also significantly improved the precision and reliability of polyploid genome assembly and phasing, such as HiCanu, hifiasm, ALLHiC, and PolyGembler. Herein, we review recently published polyploid genomes and compare the various sequencing, assembly, and phasing approaches that are utilized in these genome studies. Finally, we anticipate that accurate and telomere-to-telomere chromosome-level assembly of polyploid genomes could ultimately become a routine procedure in the near future.

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.

Pleiotropic functions of chordin gene causing drastic morphological changes in ornamental goldfish

Breeders and fanciers have established many peculiar morphological phenotypes in ornamental goldfish. Among them, the twin-tail and dorsal-finless phenotypes have particularly intrigued early and recent researchers, as equivalent morphologies are extremely rare in nature. These two mutated phenotypes appeared almost simultaneously within a short time frame and were fixed in several strains. However, little is known about how these two different mutations could have co-occurred during such a short time period. Here, we demonstrate that the chordin gene, a key factor in dorsal-ventral patterning, is responsible not only for the twin-tail phenotype but also for the dorsal-finless phenotype. Our F2 backcrossing and functional analyses revealed that the penetrance/expressivity of the dorsal-finless phenotype can be suppressed by the wild-type allele of chdS. Based on these findings, we propose that chdSwt may have masked the expression of the dorsal-finless phenotype, acting as a capacitor buffering gene to allow accumulation of genetic mutations. Once this gene lost its original function in the twin-tail goldfish lineages, the dorsal-finless phenotype could be highly expressed. Thus, this study experimentally demonstrates that the rapid genetic fixation of morphological mutations during a short domestication time period may be related to the robustness of embryonic developmental mechanisms.

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.

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.

Midnight/midday-synchronized expression of cryptochrome genes in the eyes of three teleost species, zebrafish, goldfish, and medaka

Photoperiodic responses are observed in many organisms living in the temperate zones. The circadian clock is involved in photoperiodic time measurement; however, the underlying molecular mechanism for detection of the day length remains unknown. We previously compared the expression profiles of the Cryptochrome(Cry) genes in the zebrafish eye and reported that Cry1ab has a double peak with variable expression duration depending on the photoperiod. In this study, to understand commonalities and differences in the photoperiodic responses of ocular Cry genes, we identified Cryptochrome genes in two other teleost species, goldfish and medaka, living in temperate zones, and measured ocular Cry mRNA levels in all of the three species, under different photoperiods (long-day [14 h light: 10 h dark] and short-day [10 h light: 14 h dark] and in constant darkness. Cry1ab mRNA levels did not show dual peaks in goldfish or medaka under the examined conditions; however, the mRNA expression profiles of many Crys were altered in all three species, depending on the day length and light condition. Based on their expression profiles, Cry mRNA peaks were classified into three groups that better synchronize to sunrise (light-on), midnight/midday (middle points of the dark/light periods), or sunset (light-off). These results suggest the presence of multiple oscillators that oscillate independently or a complex oscillator in which Cry expression cycles change in a photoperiod-dependent manner in the eye.

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.

Threats at home? Assessing the potential ecological impacts and risks of commonly traded pet fishes

Invasive alien species (IAS) are major drivers of global biodiversity loss, and the poorly regulated international pet trade is a source of emerging and future invaders. Predictions of the likely ecological impacts and risks of such IAS have been significantly enhanced in recent years with new metrics, which require application to many more actual and potential IAS. Hence, this study assesses the potential ecological impacts and risks of two readily available pet trade species: goldfish, Carassius auratus, a species with non-native populations worldwide; and white cloud mountain minnow, Tanichthys albonubes, a species with a limited invasion history to date. First, we compared the per capita feeding rates of these non-native species with two European trophically analogous natives – the stone loach, Barbatula barbatula, and the common minnow, Phoxinus phoxinus – using the Comparative Functional Response method. Second, we used foraging experiments in conspecific pairs to determine synergistic, neutral or antagonistic intraspecific interactions. Third, we performed novel object experiments using the two pet trade species to assess boldness, a known “dispersal enhancing trait”. Goldfish had the highest maximum feeding rates of the four species, while white cloud mountain minnows had the lowest. Neutral interactions were observed for all four species in the paired foraging experiments, with goldfish having the highest consumption and white cloud mountain minnows having the lowest. Goldfish demonstrated greater boldness, being more active during the experimental trials and more likely to approach a novel object than white cloud mountain minnows. Further, combining maximum feeding rates, boldness and species availabilities from our survey of pet shops, we assessed the relative invasion risks (RIR) of the two non-natives. This highlighted goldfish as the higher risk and most worthy of management prioritisation, mirroring its more extensive invasion history. We propose that such metrics have potential to direct future IAS policy decisions and management towards the ever-increasing rates of biological invasions worldwide.

Lactobacillus casei (IBRC-M 10,711) ameliorates the growth retardation, oxidative stress, and immunosuppression induced by malathion toxicity in goldfish (Carassius auratus)

Probiotics can functionality improve fish wellbeing and are suggested as antioxidative agents to protect fish from xenobiotics toxicity. Herein, dietary Lactobacillus casei (IBRC-M 10,711) was included in the diets of goldfish (Carassius auratus) to protect against malathion toxicity. Fish (12.47 ± 0.06 g) were randomly allocated to six groups (triplicates), as follows: T1) control; T2) fish exposed to 50% of malathion 96 h LC50; T3) L. casei at 106 CFU/g diet; T4) L. casei at 107 CFU/g diet; T5) fish exposed to 50% of malathion 96 h LC50 + L. casei at 106 CFU/g diet; T6) fish exposed to 50% of malathion 96 h LC50 + L. casei at 107 CFU/g diet. After 60 days, goldfish fed T4 had the highest final body weight (FBW), weight gain (WG), and specific growth rate (SGR), and the lowest feed conversion ratio (FCR) among the groups (P < 0.05). However, the T2 group showed lower FBW, WG, and SGR and higher FCR than fish in T1 (P < 0.05). Fish in the T4 group had the highest blood total proteins, albumin, and globulin, while fish in T2 had the lowest levels (P < 0.05). Fish in the group T2 had the highest triglycerides, cholesterol, cortisol, lactate dehydrogenase (LDH), alanine aminotransferase (ALT), aspartate aminotransferase (AST), and alkaline phosphatase (ALP) levels in the blood, while fish fed T4 had the lowest values (P < 0.05). The superoxide dismutase (SOD) and catalase (CAT) showed the highest activities in T3 and T4 groups, and the lowest SOD was seen in the T2 group, whereas the lowest CAT was seen in the T2, T5, and T6 groups (P < 0.05). Fish in the T5 and T6 groups had higher glutathione peroxidase (GSH-Px) activities than fish in T1 and T2 groups but T3 and T4 groups showed the highest values (P < 0.05). T2 group had the highest malondialdehyde (MDA) level, while T3 and T4 groups had the lowest MDA level (P < 0.05). Blood immunoglobulin (Ig) and lysozyme activity were significantly higher in T3 and T4 groups and lower in the T2 group than in the control (P < 0.05). The alternative complement pathway (ACH50) was significantly higher in T2, T3, T4, T5, and T6 groups than in the T1 group (P < 0.05). Skin mucus Ig was significantly higher in T3 and T4 groups and lower in the T2 group than in the control (P < 0.05). The highest lysozyme activity, protease, and ACH50 in the skin mucus samples were in the T4 group, while the lowest values were in the T2 group (P < 0.05). In conclusion, dietary L. casei protects goldfish from malathion-induced growth retardation, oxidative stress, and immunosuppression.

Genomic consequences of domestication of the Siamese fighting fish

Siamese fighting (betta) fish are among the most popular and morphologically diverse pet fish, but the genetic bases of their domestication and phenotypic diversification are largely unknown. We assembled de novo the genome of a wild Betta splendens and whole-genome sequenced 98 individuals across five closely related species. We find evidence of bidirectional hybridization between domesticated ornamental betta and other wild Betta species. We discover dmrt1 as the main sex determination gene in ornamental betta and that it has lower penetrance in wild B. splendens. Furthermore, we find genes with signatures of recent, strong selection that have large effects on color in specific parts of the body or on the shape of individual fins and that most are unlinked. Our results demonstrate how simple genetic architectures paired with anatomical modularity can lead to vast phenotypic diversity generated during animal domestication and launch betta as a powerful new system for evolutionary genetics.

Domestication of farmed fish via the attenuation of stress responses mediated by the hypothalamus-pituitary-inter-renal endocrine axis

Human-directed domestication of terrestrial animals traditionally requires thousands of years for breeding. The most prominent behavioral features of domesticated animals include reduced aggression and enhanced tameness relative to their wild forebears, and such behaviors improve the social tolerance of domestic animals toward both humans and crowds of their own species. These behavioral responses are primarily mediated by the hypothalamic-pituitary-adrenal (inter-renal in fish) (HPA/I) endocrine axis, which is involved in the rapid conversion of neuronal-derived perceptual information into hormonal signals. Over recent decades, growing evidence implicating the attenuation of the HPA/I axis during the domestication of animals have been identified through comprehensive genomic analyses of the paleogenomic datasets of wild progenitors and their domestic congeners. Compared with that of terrestrial animals, domestication of most farmed fish species remains at early stages. The present review focuses on the application of HPI signaling attenuation to accelerate the domestication and genetic breeding of farmed fish. We anticipate that deeper understanding of HPI signaling and its implementation in the domestication of farmed fish will benefit genetic breeding to meet the global demands of the aquaculture industry.

Parallel subgenome structure and divergent expression evolution of allo-tetraploid common carp and goldfish

How two subgenomes in allo-tetraploids adapt to coexistence and coordinate through structure and expression evolution requires extensive studies. In the present study, we report an improved genome assembly of allo-tetraploid common carp, an updated genome annotation of allo-tetraploid goldfish and the chromosome-scale assemblies of a progenitor-like diploid Puntius tetrazona and an outgroup diploid Paracanthobrama guichenoti. Parallel subgenome structure evolution in the allo-tetraploids was featured with equivalent chromosome components, higher protein identities, similar transposon divergence and contents, homoeologous exchanges, better synteny level, strong sequence compensation and symmetric purifying selection. Furthermore, we observed subgenome expression divergence processes in the allo-tetraploids, including inter-/intrasubgenome trans-splicing events, expression dominance, decreased expression levels, dosage compensation, stronger expression correlation, dynamic functionalization and balancing of differential expression. The potential disorders introduced by different progenitors in the allo-tetraploids were hypothesized to be alleviated by increasing structural homogeneity and performing versatile expression processes. Resequencing three common carp strains revealed two major ecotypes and uncovered candidate genes relevant to growth and survival rate.

Susceptibility of Goldfish to Cyprinid Herpesvirus 2 (CyHV-2) SH01 Isolated from Cultured Crucian Carp

Cyprinid herpesvirus 2 (CyHV-2), a member of the Alloherpesviridae family belonging to the genus Cyprinivirus, is a fatal contagious aquatic pathogen that affects goldfish (Carassius auratus) and crucian carp (Carassius carassius). Although crucian carp and goldfish belong to the genus Carassius, it is unclear whether they are susceptible to the same CyHV-2 isolate. In addition, the origin of the crucian carp-derived CyHV-2 virus isolate remains unclear. CyHV-2 SH01 was isolated during herpesviral hematopoietic necrosis disease (HVHN) outbreaks in crucian carp at a local fish farm near Shanghai. CyHV-2 SH01 was confirmed by PCR and Western blot analysis of kidney, spleen, muscle, and blood tissue from the diseased crucian carp. Moreover, histopathological and ultra-pathological analyses revealed pathological changes characteristic of CyHV-2 SH01 infection in the tissues of the diseased crucian carp. In the present study, goldfish and crucian carp were challenged with CyHV-2 SH01 to elucidate viral virulence. We found that CyHV-2 SH01 could cause rapid and fatal disease progression in goldfish and crucian carp 24 h post-injection at 28 °C. Experimental infection of goldfish by injection indicated that the average virus titer in the kidney of the goldfish was 10^3.47 to 10^3.59 copies/mg. In addition, tissues exhibited the most prominent histopathological changes (cellular wrinkling and shrinkage, cytoplasmic vacuolation, fusion of the gill lamellae, and hepatic congestion) in CyHV-2 SH01-infected goldfish and crucian carp. Thus, crucian carp and goldfish showed a high sensitivity, with typical symptoms, to HVHN disease caused by CyHV-2 SH01.

Beyond fish eDNA metabarcoding: Field replicates disproportionately improve the detection of stream associated vertebrate species

Fast, reliable, and comprehensive biodiversity monitoring data are needed for environmental decision making and management. Recent work on fish environmental DNA (eDNA) metabarcoding shows that aquatic diversity can be captured fast, reliably, and non-invasively at moderate costs. Because water in a catchment flows to the lowest point in the landscape, often a stream, it can collect traces of terrestrial species via surface or subsurface runoff along its way or when specimens come into direct contact with water (e.g., when drinking). Thus, fish eDNA metabarcoding data can provide information on fish but also on other vertebrate species that live in riparian habitats. This additional data may offer a much more comprehensive approach for assessing vertebrate diversity at no additional costs. Studies on how the sampling strategy affects species detection especially of stream-associated communities, however, are scarce. We therefore performed an analysis on the effects of biological replication on both fish as well as (semi-)terrestrial species detection. Along a 2 km stretch of the river Mulde (Germany), we collected 18 1-L water samples and analyzed the relation of detected species richness and quantity of biological replicates taken. We detected 58 vertebrate species, of which 25 were fish and lamprey, 18 mammals, and 15 birds, which account for 50%, 22.2%, and 7.4% of all native species to the German federal state of Saxony-Anhalt. However, while increasing the number of biological replicates resulted in only 24.8% more detected fish and lamprey species, mammal, and bird species richness increased disproportionately by 68.9% and 77.3%, respectively. Contrary, PCR replicates showed little stochasticity. We thus emphasize to increase the number of biological replicates when the aim is to improve general species detections. This holds especially true when the focus is on rare aquatic taxa or on (semi-)terrestrial species, the so-called ‘bycatch’. As a clear advantage, this information can be obtained without any additional sampling or laboratory effort when the sampling strategy is chosen carefully. With the increased use of eDNA metabarcoding as part of national fish bioassessment and monitoring programs, the complimentary information provided on bycatch can be used for biodiversity monitoring and conservation on a much broader scale.

The Metabolomic Response of Crucian Carp (Carassius carassius) to Anoxia and Reoxygenation Differs between Tissues and Hints at Uncharacterized Survival Strategies

The anoxia-tolerant crucian carp (Carassius carassius) has been studied in detail for numerous years, with particular focus on unravelling the underlying physiological mechanisms of anoxia tolerance. However, relatively little work has been focused on what occurs beyond anoxia, and often the focus is a single organ or tissue type. In this study, we quantified more than 100 metabolites by capillary electrophoresis-mass spectrometry (CE-MS) in brain, heart, liver, and blood plasma from four experimental groups, being normoxic (control) fish, anoxia-exposed fish, and two groups that had been exposed to anoxia followed by reoxygenation for either 3 h or 24 h. The heart, which maintains cardiac output during anoxia, unexpectedly, was slower to recover compared to the brain and liver, mainly due to a slower return to control concentrations of the energy-carrying compounds ATP, GTP, and phosphocreatine. Crucian carp accumulated amino acids in most tissues, and also surprisingly high levels of succinate in all tissues investigated during anoxia. Purine catabolism was enhanced, leading to accumulation of uric acid during anoxia and increasing urea formation that continued into 24 h of reoxygenation. These tissue-specific differences in accumulation and distribution of the metabolites may indicate an intricate system of transport between tissues, opening for new avenues of investigation of possible mechanisms aimed at reducing the generation of reactive oxygen species (ROS) and resultant tissue damage during reoxygenation.

Abnormal meiosis in fertile and sterile triploid cyprinid fish

Meiosis is the key process for producing mature gametes. A natural fertile triploid Carassius auratus population (3nDTCC) and an artificially derived sterile triploid crucian carp (3nCC) have been previously observed, providing suitable model organisms for investigating meiosis characteristics in triploid fish. In the present study, the microstructures and ultrastructures of spermatogenesis were studied in these fishes. TdT-mediated dUTP nick end labeling detection was performed to investigate the apoptosis of spermatocytes. Fluorescence in situ hybridization was employed to trace chromatin pairing. In addition, the mRNA expressions of cell cycle-related genes (i.e., cell division control 2 and cell cycle protein B) were determined by quantitative realtime polymerase chain reaction to illustrate the molecular mechanism of abnormal meiosis in the 3nCC. The results showed that the 3nCC undergoes an irregular prophase I, with the chromosomes distributed in a unipolar radial manner and exhibiting partial pairing, hindered metaphase I, and degenerated cells in the subsequent stages. Meanwhile, the 3nDTCC presented a relatively regular meiotic prophase I with complete conjugate chromosome pairs and chromosomes distributed along the karyotheca, which were presented as a ring structure by slicing. Only the spreads with 130-150 irregular chromosomes can be easily detected in the 3nDTCC, suggesting that it may undergo an abnormal metaphase I. This study provides new insights into the meiosis of fertile and sterile triploid cyprinid fish.