The Genetic Basis of Morphological Diversity in Domesticated Goldfish

A Comparative Transcriptome and Proteome Analysis of the Molecular Mechanism Underlying Anterior to Dorsal Eye Rotation in the Celestial-Eye Goldfish (Carassius auratus)

Goldfish (Carassius auratus), subjected to millennia of artificial selection and breeding, have diversified into numerous ornamental varieties, such as the celestial-eye (CE) goldfish, noted for its unique dorsal eye rotation. Previous studies have primarily focused on anatomical modifications in CE goldfish eyes, yet the molecular underpinnings of their distinctive eye orientation remain poorly understood. This study employed high-throughput transcriptome and proteome sequencing on 110-day-old full-sibling CE goldfish, which displayed either anterior or upward eye rotations. Verification of these findings was conducted using quantitative PCR (qPCR) for transcriptomic data and parallel reaction monitoring (PRM) for proteomic analysis. Our research identified 73,685 genes and 7717 proteins, pinpointing 8 common differentially expressed genes (DEGs) and proteins (DEPs) implicated in cytoskeleton remodeling, cell adhesion, apoptosis, and optic nerve regeneration. Enrichment analyses further delineated pathways associated with apoptosis, necroptosis, and cell adhesion molecules. The results indicated a significant role for genes involved in cytoskeletal dynamics, nervous system function, and apoptotic processes in the dorsal eye rotation of CE goldfish. Analyses of abnormalities in ocular membrane structures, along with disturbances in lipid and protein synthesis metabolism and energy metabolism during developmental stages, provided compelling evidence for the potential use of CE goldfish as a model organism in studying human eye-related disorders. This investigation provided the first comprehensive transcriptomic and proteomic overview of eye rotation in CE goldfish, offering insights crucial for the genetic breeding of new ornamental fish varieties.

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.

An improved transcriptome annotation reveals asymmetric expression and distinct regulation patterns in allotetraploid common carp

In allotetraploid common carp, protein-coding homoeologs presented divergent expression levels between the two subgenomes. However, whether subgenome dominance occurs in other transcriptional and post-transcriptional events remains unknown. Using Illumina RNA sequencing and PacBio full-length sequencing, we refined the common carp transcriptome annotation and explored differences in four transcriptional and post-transcriptional events between the two subgenomes. The results revealed that the B subgenome presented more alternative splicing events, as did lncRNAs and circRNAs. However, the expression levels, tissue specificity, sequence features, and functions of lncRNAs and circRNAs did not significantly differ between the two subgenomes, suggesting a common regulatory mechanism shared by the two subgenomes. Furthermore, both the number and base substitution frequency of RNA editing events were greater in the B subgenome. Functional analyses of these transcriptional events also revealed subgenome bias. Genes that undergo alternative splicing in the A subgenome participate in more biological processes, and lncRNA targets show a preference between subgenomes. CircRNA host genes in the B subgenome were associated with more biological functions, and RNA editing preferentially occurred in noncoding regions or led to nonsynonymous mutations in the B subgenome. Taken together, the refined transcriptome annotation revealed complicated and imbalanced expression strategies in allotetraploid common carp.

A comparative analysis of planarian genomes reveals regulatory conservation in the face of rapid structural divergence

The planarian Schmidtea mediterranea is being studied as a model species for regeneration, but the assembly of planarian genomes remains challenging. Here, we report a high-quality haplotype-phased, chromosome-scale genome assembly of the sexual S2 strain of S. mediterranea and high-quality chromosome-scale assemblies of its three close relatives, S. polychroa, S. nova, and S. lugubris. Using hybrid gene annotations and optimized ATAC-seq and ChIP-seq protocols for regulatory element annotation, we provide valuable genome resources for the planarian research community and a first comparative perspective on planarian genome evolution. Our analyses reveal substantial divergence in protein-coding sequences and regulatory regions but considerable conservation within promoter and enhancer annotations. We also find frequent retrotransposon-associated chromosomal inversions and interchromosomal translocations within the genus Schmidtea and, remarkably, independent and nearly complete losses of ancestral metazoan synteny in Schmidtea and two other flatworm groups. Overall, our results suggest that platyhelminth genomes can evolve without syntenic constraints.

Phylogenomic analyses of all species of swordtail fishes (genus Xiphophorus) show that hybridization preceded speciation

Hybridization has been recognized to play important roles in evolution, however studies of the genetic consequence are still lagging behind in vertebrates due to the lack of appropriate experimental systems. Fish of the genus Xiphophorus are proposed to have evolved with multiple ancient and ongoing hybridization events. They have served as an informative research model in evolutionary biology and in biomedical research on human disease for more than a century. Here, we provide the complete genomic resource including annotations for all described 26 Xiphophorus species and three undescribed taxa and resolve all uncertain phylogenetic relationships. We investigate the molecular evolution of genes related to cancers such as melanoma and for the genetic control of puberty timing, focusing on genes that are predicted to be involved in pre-and postzygotic isolation and thus affect hybridization. We discovered dramatic size-variation of some gene families. These persisted despite reticulate evolution, rapid speciation and short divergence time. Finally, we clarify the hybridization history in the entire genus settling disputed hybridization history of two Southern swordtails. Our comparative genomic analyses revealed hybridization ancestries that are manifested in the mosaic fused genomes and show that hybridization often preceded speciation.

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.

Hagfish genome elucidates vertebrate whole-genome duplication events and their evolutionary consequences

Polyploidy or whole-genome duplication (WGD) is a major event that drastically reshapes genome architecture and is often assumed to be causally associated with organismal innovations and radiations. The 2R hypothesis suggests that two WGD events (1R and 2R) occurred during early vertebrate evolution. However, the timing of the 2R event relative to the divergence of gnathostomes (jawed vertebrates) and cyclostomes (jawless hagfishes and lampreys) is unresolved and whether these WGD events underlie vertebrate phenotypic diversification remains elusive. Here we present the genome of the inshore hagfish, Eptatretus burgeri. Through comparative analysis with lamprey and gnathostome genomes, we reconstruct the early events in cyclostome genome evolution, leveraging insights into the ancestral vertebrate genome. Genome-wide synteny and phylogenetic analyses support a scenario in which 1R occurred in the vertebrate stem-lineage during the early Cambrian, and 2R occurred in the gnathostome stem-lineage, maximally in the late Cambrian-earliest Ordovician, after its divergence from cyclostomes. We find that the genome of stem-cyclostomes experienced an additional independent genome triplication. Functional genomic and morphospace analyses demonstrate that WGD events generally contribute to developmental evolution with similar changes in the regulatory genome of both vertebrate groups. However, appreciable morphological diversification occurred only in the gnathostome but not in the cyclostome lineage, calling into question the general expectation that WGDs lead to leaps of bodyplan complexity.

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.

Identification of Candidate Genes for Red-Eyed (Albinism) Domestic Guppies Using Genomic and Transcriptomic Analyses

Guppies are small tropical fish with brightly colored bodies and variable tail shapes. There are two phenotypes of domestic guppy eye color: red and black. The wild type is black-eyed. The main object of this study was to identify candidate genes for the red-eyed phenotype in domestic guppies. We hope to provide molecular genetic information for the development of new domestic guppy strains. Additionally, the results also contribute to basic research concerning guppies. In this study, 121 domestic guppies were used for genomic analysis (GWAS), and 44 genes were identified. Furthermore, 21 domestic guppies were used for transcriptomic analysis, and 874 differentially expressed genes (DEGs) were identified, including 357 upregulated and 517 downregulated genes. Through GO and KEGG enrichment, we identified some important terms or pathways mainly related to melanin biosynthesis and ion transport. qRT-PCR was also performed to verify the differential expression levels of four important candidate genes (TYR, OCA2, SLC45A2, and SLC24A5) between red-eyed and black-eyed guppies. Based on the results of genomic and transcriptomic analyses, we propose that OCA2 is the most important candidate gene for the red-eyed phenotype in guppies.

Phylogenomics analyses of all species of Swordtails (Genus Xiphophorus ) highlights hybridization precedes speciation

Hybridization has been recognized as an important driving force for evolution, however studies of the genetic consequence and its cause are still lagging behind in vertebrates due to the lack of appropriate experimental systems. Fish of the central American genus Xiphophorus were proposed to have evolved with multiple ancient and ongoing hybridization events, and served as a valuable research model in evolutionary biology and in biomedical research on human disease for more than a century. Here, we provide the complete genome resource and its annotation of all 26 Xiphophorus species. On this dataset we resolved the so far conflicting phylogeny. Through comparative genomic analyses we investigated the molecular evolution of genes related to melanoma, for a main sexually selected trait and for the genetic control of puberty timing, which are predicted to be involved in pre-and postzygotic isolation and thus to influence the probability of interspecific hybridization in Xiphophorus . We demonstrate dramatic size-variation of some gene families across species, despite the reticulate evolution and short divergence time. Finally, we clarify the hybridization history in the genus Xiphophorus genus, settle the long dispute on the hybridization origin of two Southern swordtails, highlight hybridizations precedes speciation, and reveal the distribution of hybridization ancestry remaining in the fused genome.

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.

Maternal dominance contributes to subgenome differentiation in allopolyploid fishes

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

Analysis of the genetic loci of pigment pattern evolution in vertebrates

Vertebrate pigmentation patterns are amongst the best characterised model systems for studying the genetic basis of adaptive evolution. The wealth of available data on the genetic basis for pigmentation evolution allows for analysis of trends and quantitative testing of evolutionary hypotheses. We employed Gephebase, a database of genetic variants associated with natural and domesticated trait variation, to examine trends in how cis-regulatory and coding mutations contribute to vertebrate pigmentation phenotypes, as well as factors that favour one mutation type over the other. We found that studies with lower ascertainment bias identified higher proportions of cis-regulatory mutations, and that cis-regulatory mutations were more common amongst animals harbouring a higher number of pigment cell classes. We classified pigmentation traits firstly according to their physiological basis and secondly according to whether they affect colour or pattern, and identified that carotenoid-based pigmentation and variation in pattern boundaries are preferentially associated with cis-regulatory change. We also classified genes according to their developmental, cellular, and molecular functions. We found a greater proportion of cis-regulatory mutations in genes implicated in upstream developmental processes compared to those involved in downstream cellular functions, and that ligands were associated with a higher proportion of cis-regulatory mutations than their respective receptors. Based on these trends, we discuss future directions for research in vertebrate pigmentation evolution.

A median fin derived from the lateral plate mesoderm and the origin of paired fins

The development of paired appendages was a key innovation during evolution and facilitated the aquatic to terrestrial transition of vertebrates. Largely derived from the lateral plate mesoderm (LPM), one hypothesis for the evolution of paired fins invokes derivation from unpaired median fins via a pair of lateral fin folds located between pectoral and pelvic fin territories1. Whilst unpaired and paired fins exhibit similar structural and molecular characteristics, no definitive evidence exists for paired lateral fin folds in larvae or adults of any extant or extinct species. As unpaired fin core components are regarded as exclusively derived from paraxial mesoderm, any transition presumes both co-option of a fin developmental programme to the LPM and bilateral duplication2. Here, we identify that the larval zebrafish unpaired pre-anal fin fold (PAFF) is derived from the LPM and thus may represent a developmental intermediate between median and paired fins. We trace the contribution of LPM to the PAFF in both cyclostomes and gnathostomes, supporting the notion that this is an ancient trait of vertebrates. Finally, we observe that the PAFF can be bifurcated by increasing bone morphogenetic protein signalling, generating LPM-derived paired fin folds. Our work provides evidence that lateral fin folds may have existed as embryonic anlage for elaboration to paired fins.

Transposon signatures of allopolyploid genome evolution

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

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

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

Genomic and Transcriptional Profiles of Kelch-like (klhl) Gene Family in Polyploid Carassius Complex

Genome duplication supplies raw genetic materials and has been thought to be essential for evolutionary innovation and ecological adaptation. Here, we select Kelch-like (klhl) genes to study the evolution of the duplicated genes in the polyploid Carassius complex, including amphidiploid C. auratus and amphitriploid C. gibelio. Phylogenetic, chromosomal location and read coverage analyses indicate that most of Carassius klhl genes exhibit a 2:1 relationship with zebrafish orthologs and confirm two rounds of polyploidy, an allotetraploidy followed by an autotriploidy, occurred during Carassius evolution. The lineage-specific expansion and biased retention/loss of klhl genes are also found in Carassius. Transcriptome analyses across eight adult tissues and seven embryogenesis stages reveal varied expression dominance and divergence between the two species. The expression of klhls in response to Carassius herpesvirus 2 infection shows different expression changes corresponding to distinct herpesvirus resistances in three C. gibelio gynogenetic clones. Finally, we find that most C. gibelio klhl genes possess three alleles except eight genes that have lost one or two alleles due to genome rearrangement. The allele expression bias is prosperous for Cgklhl genes and varies during embryogenesis owning to the sequential expression manner of the alleles. The current study provides global insights into the genomic and transcriptional evolution of duplicated genes in a given superfamily resulting from multiple rounds of polyploidization.

Bioelectricity in Developmental Patterning and Size Control: Evidence and Genetically Encoded Tools in the Zebrafish Model

Developmental patterning is essential for regulating cellular events such as axial patterning, segmentation, tissue formation, and organ size determination during embryogenesis. Understanding the patterning mechanisms remains a central challenge and fundamental interest in developmental biology. Ion-channel-regulated bioelectric signals have emerged as a player of the patterning mechanism, which may interact with morphogens. Evidence from multiple model organisms reveals the roles of bioelectricity in embryonic development, regeneration, and cancers. The Zebrafish model is the second most used vertebrate model, next to the mouse model. The zebrafish model has great potential for elucidating the functions of bioelectricity due to many advantages such as external development, transparent early embryogenesis, and tractable genetics. Here, we review genetic evidence from zebrafish mutants with fin-size and pigment changes related to ion channels and bioelectricity. In addition, we review the cell membrane voltage reporting and chemogenetic tools that have already been used or have great potential to be implemented in zebrafish models. Finally, new perspectives and opportunities for bioelectricity research with zebrafish are discussed.

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

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

The genetic basis and potential molecular mechanism of yellow-albino northern snakehead (Channa argus)

Body colour is an important economic trait for commercial fishes. Recently, a new colour morph displaying market-favoured yellow skin (termed as yellow-mutant, YM) of northern snakehead (Channa argus) was discovered in China. We confirmed that YM snakehead is an albino with complete loss of melanin in the skin and eyes by histological and ultrastructural observations, and inherited as a recessive Mendelian trait. By applying genomic analysis approaches, in combination with gene knockdown and rescue experiments, we suggested a non-sense mutation in slc45a2 (c.383G > A) is the causation for the YM snakehead. Notably, significantly higher levels of key melanogenesis genes (tyr, tyrp1, dct and pmel) and phospho-MITF protein were detected in YM snakehead than those in wild-type individuals, and the underlying mechanism was further investigated by comparative transcriptomic analysis. Results revealed that differential expressed genes involved in pathways like MAPK, WNT and calcium signalling were significantly induced in YM snakehead, which might account for the increased amount of melanogenesis elements, and presumably be stimulated by fibroblast-derived melanogenic factors in a paracrine manner. Our study clarified the genetic basis of colour variation in C. argus and provided the preliminary clue indicating the potential involvement of fibroblasts in pigmentation in fish.

Single-cell transcriptomics of the goldfish retina reveals genetic divergence in the asymmetrically evolved subgenomes after allotetraploidization

The recent whole-genome duplication (WGD) in goldfish (Carassius auratus) approximately 14 million years ago makes it a valuable model for studying gene evolution during the early stages after WGD. We analyzed the transcriptome of the goldfish retina at the level of single-cell (scRNA-seq) and open chromatin regions (scATAC-seq). We identified a group of genes that have undergone dosage selection, accounting for 5% of the total 11,444 ohnolog pairs. We also identified 306 putative sub/neo-functionalized ohnolog pairs that are likely to be under cell-type-specific genetic variation at single-cell resolution. Diversification in the expression patterns of several ohnolog pairs was observed in the retinal cell subpopulations. The single-cell level transcriptome analysis in this study uncovered the early stages of evolution in retinal cell of goldfish after WGD. Our results provide clues for understanding the relationship between the early stages of gene evolution after WGD and the evolution of diverse vertebrate retinal functions.

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.

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

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

Symmetric subgenomes and balanced homoeolog expression stabilize the establishment of allopolyploidy in cyprinid fish

BACKGROUND

Interspecific postzygotic reproduction isolation results from large genetic divergence between the subgenomes of established hybrids. Polyploidization immediately after hybridization may reset patterns of homologous chromosome pairing and ameliorate deleterious genomic incompatibility between the subgenomes of distinct parental species in plants and animals. However, the observation that polyploidy is less common in vertebrates raises the question of which factors restrict its emergence. Here, we perform analyses of the genome, epigenome, and gene expression in the nascent allotetraploid lineage (2.95 Gb) derived from the intergeneric hybridization of female goldfish (Carassius auratus, 1.49 Gb) and male common carp (Cyprinus carpio, 1.42 Gb), to shed light on the changes leading to the stabilization of hybrids.

RESULTS

We firstly identify the two subgenomes derived from the parental lineages of goldfish and common carp. We find variable unequal homoeologous recombination in somatic and germ cells of the intergeneric F1 and allotetraploid (F22 and F24) populations, reflecting high plasticity between the subgenomes, and rapidly varying copy numbers between the homoeolog genes. We also find dynamic changes in transposable elements accompanied by genome merger and duplication in the allotetraploid lineage. Finally, we observe the gradual decreases in cis-regulatory effects and increases in trans-regulatory effects along with the allotetraploidization, which contribute to increases in the symmetrical homoeologous expression in different tissues and developmental stages, especially in early embryogenesis.

CONCLUSIONS

Our results reveal a series of changes in transposable elements, unequal homoeologous recombination, cis- and trans-regulations (e.g. DNA methylation), and homoeologous expression, suggesting their potential roles in mediating adaptive stabilization of regulatory systems of the nascent allotetraploid lineage. The symmetrical subgenomes and homoeologous expression provide a novel way of balancing genetic incompatibilities, providing a new insight into the early stages of allopolyploidization in vertebrate evolution.

Coordinated patterning of zebrafish caudal fin symmetry by a central and two peripheral organizers

BACKGROUND

Caudal fin symmetry characterizes teleosts and likely contributes to their evolutionary success. However, the coordinated development and patterning of skeletal elements establishing external symmetry remains incompletely understood. We explore the spatiotemporal emergence of caudal skeletal elements in zebrafish to consider evolutionary and developmental origins of caudal fin symmetry.

RESULTS

Transgenic reporters and skeletal staining reveal that the hypural diastema-defining gap between hypurals 2 and 3 forms early and separates progenitors of two plates of connective tissue. Two sets of central principal rays (CPRs) synchronously, sequentially, and symmetrically emerge around the diastema. The two dorsal- and ventral-most rays (peripheral principal rays, PPRs) arise independently and earlier than adjacent CPRs. Muscle and tendon markers reveal that different muscles attach to CPR and PPR sets.

CONCLUSIONS

We propose that caudal fin symmetry originates from a central organizer that establishes the hypural diastema and bidirectionally patterns surrounding tissue into two plates of connective tissue and two mirrored sets of CPRs. Further, two peripheral organizers unidirectionally specify PPRs, forming a symmetric "composite" fin derived from three fields. Distinct CPR and PPR ontogenies may represent developmental modules conferring ray identities, muscle connections, and biomechanical properties. Our model contextualizes mechanistic studies of teleost fin morphological variation.

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.

Rodlet cells in kidney of goldfish (Carassius auratus, Linnaeus 1758): A light and confocal microscopy study

Rodlet cells (RCs) have always been an enigma for scientists. RCs have been given a variety of activities over the years, including ion transport, osmoregulation, and sensory function. These cells, presumably as members of the granulocyte line, are present only in teleosts and play a role in the innate immune response. RCs are migratory cells found in a variety of organs, including skin, vascular, digestive, uropoietic, reproductive, and respiratory systems, and present distinct physical properties that make them easily recognizable in tissues and organs. The development of RCs can be divided into four stages: granular, transitional, mature, and ruptured, having different morphological characteristics. Our study aims to characterize the different stages of these cells by histomorphological and histochemical techniques. Furthermore, we characterized these cells at all stages with peroxidase and fluorescence immunohistochemical techniques using different antibodies: S100, tubulin, α-SMA, piscidin, and for the first time TLR-2. From our results, the immunoreactivity of these cells to the antibodies performed may confirm that RCs play a role in fish defense mechanisms, helping to expand the state of the art on immunology and immune cells of teleosts.

Modulation of bioelectric cues in the evolution of flying fishes

Changes to allometry, or the relative proportions of organs and tissues within organisms, is a common means for adaptive character change in evolution. However, little is understood about how relative size is specified during development and shaped during evolution. Here, through a phylogenomic analysis of genome-wide variation in 35 species of flying fishes and relatives, we identify genetic signatures in both coding and regulatory regions underlying the convergent evolution of increased paired fin size and aerial gliding behaviors. To refine our analysis, we intersected convergent phylogenomic signatures with mutants with altered fin size identified in distantly related zebrafish. Through these paired approaches, we identify a surprising role for an L-type amino acid transporter, lat4a, and the potassium channel, kcnh2a, in the regulation of fin proportion. We show that interaction between these genetic loci in zebrafish closely phenocopies the observed fin proportions of flying fishes. The congruence of experimental and phylogenomic findings point to conserved, non-canonical signaling integrating bioelectric cues and amino acid transport in the establishment of relative size in development and evolution.

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.

Genomic Basis of Striking Fin Shapes and Colors in the Fighting Fish

Resolving the genomic basis underlying phenotypic variations is a question of great importance in evolutionary biology. However, understanding how genotypes determine the phenotypes is still challenging. Centuries of artificial selective breeding for beauty and aggression resulted in a plethora of colors, long-fin varieties, and hyper-aggressive behavior in the air-breathing Siamese fighting fish (Betta splendens), supplying an excellent system for studying the genomic basis of phenotypic variations. Combining whole-genome sequencing, quantitative trait loci mapping, genome-wide association studies, and genome editing, we investigated the genomic basis of huge morphological variation in fins and striking differences in coloration in the fighting fish. Results revealed that the double tail, elephant ear, albino, and fin spot mutants each were determined by single major-effect loci. The elephant ear phenotype was likely related to differential expression of a potassium ion channel gene, kcnh8. The albinotic phenotype was likely linked to a cis-regulatory element acting on the mitfa gene and the double-tail mutant was suggested to be caused by a deletion in a zic1/zic4 coenhancer. Our data highlight that major loci and cis-regulatory elements play important roles in bringing about phenotypic innovations and establish Bettas as new powerful model to study the genomic basis of evolved changes.

Loss of αklotho causes reduced motor ability and short lifespan in zebrafish

The klotho gene encodes a transmembrane protein αKlotho that interacts with a fibroblast growth factor (FGF) receptor in renal tubular epithelial cells and functions as a co-receptor for FGF23, which is an osteocytes-derived hormone. This bone-to-kidney signal promotes urinary phosphate excretion. Interestingly, αKlotho knockout mice show an accelerated aging and a shortened life span. Similarly, C. elegans lacking the αklotho homologue showed a short life span. However, the physiological basis of aging-related function of αklotho remain unclear. The αklotho-deficient vertebrate animals other than mice have been awaited as an alternative model of premature aging. We here employed zebrafish in our study and revealed that αklotho mutant zebrafish appeared to be normal at 3 months postfertilization (mpf) but eventually underwent premature death by 9 mpf, while normal zebrafish is known to survive for 42 months. We also assessed the motor ability of zebrafish in a forced swimming assay and found that αklotho mutant zebrafish displayed reduced swimming performance before their survival declined. A recent study also reported a similar finding that αklotho-deficient zebrafish exhibited a short life span and reduced spontaneous movements. Taken together, these results suggest that αKlotho mutant zebrafish show premature aging and are useful to investigate aging in vertebrates.

A novel application of bubble-eye strain of Carassius auratus for ex vivo fish immunological studies

In this study, we investigated a new application of bubble-eye goldfish (commercially available strain with large bubble-shaped eye sacs) for immunological studies in fishes utilizing the technical advantage of examining immune cells in the eye sac fluid ex vivo without sacrificing animals. As known in many aquatic species, the common goldfish strain showed an increased infection sensitivity at elevated temperature, which we demonstrate may be due to an immune impairment using the bubble-eye goldfish model. Injection of heat-killed bacterial cells into the eye sac resulted in an inflammatory symptom (surface reddening) and increased gene expression of pro-inflammatory cytokines observed in vivo, and elevated rearing temperature suppressed the induction of pro-inflammatory gene expressions. We further conducted ex vivo experiments using the immune cells harvested from the eye sac and found that the induced expression of pro-inflammatory cytokines was suppressed when we increased the temperature of ex vivo culture, suggesting that the temperature response of the eye-sac immune cells is a cell autonomous function. These results indicate that the bubble-eye goldfish is a suitable model for ex vivo investigation of fish immune cells and that the temperature-induced infection susceptibility in the goldfish may be due to functional impairments of immune cells.

Genome Evolution: Domestication of the Allopolyploid Goldfish

Goldfish are popular ornamental animals with morphologically highly diverse strains generated by artificial selection over the past millennium. New genome analyses reveal the genetics underlying some of the most iconic goldfish phenotypes and illuminate the domestication of these diverse strains following genome duplication.

Igf signaling couples retina growth with body growth by modulating progenitor cell division

How the body and organs balance their relative growth is of key importance for coordinating size and function. This is of particular relevance in organisms, which continue to grow over their entire life span. We addressed this issue in the neuroretina of medaka fish (Oryzias latipes), a well-studied system with which to address vertebrate organ growth. We reveal that a central growth regulator, Igf1 receptor (Igf1r), is necessary and sufficient for proliferation control in the postembryonic retinal stem cell niche: the ciliary marginal zone (CMZ). Targeted activation of Igf1r signaling in the CMZ uncouples neuroretina growth from body size control, and we demonstrate that Igf1r operates on progenitor cells, stimulating their proliferation. Activation of Igf1r signaling increases retinal size while preserving its structural integrity, revealing a modular organization in which progenitor differentiation and neurogenesis are self-organized and highly regulated. Our findings position Igf signaling as a key module for controlling retinal size and composition, with important evolutionary implications.

Highlighted Article: Targeted activation of Igf1r signaling in the retinal stem cell niche increases retina size through expanding the progenitor but not stem cell population.

The Developmental and Genetic Architecture of the Sexually Selected Male Ornament of Swordtails

Sexual selection results in sex-specific characters like the conspicuously pigmented extension of the ventral tip of the caudal fin-the "sword"-in males of several species of Xiphophorus fishes. To uncover the genetic architecture underlying sword formation and to identify genes that are associated with its development, we characterized the sword transcriptional profile and combined it with genetic mapping approaches. Results showed that the male ornament of swordtails develops from a sexually non-dimorphic prepattern of transcription factors in the caudal fin. Among genes that constitute the exclusive sword transcriptome and are located in the genomic region associated with this trait we identify the potassium channel, Kcnh8, as a sword development gene. In addition to its neural function kcnh8 performs a known role in fin growth. These findings indicate that during evolution of swordtails a brain gene has been co-opted for an additional novel function in establishing a male ornament.

Global tissue transcriptomic analysis to improve genome annotation and unravel skin pigmentation in goldfish

Goldfish is an ornamental fish with diverse phenotypes. However, the limited genomic resources of goldfish hamper our understanding of the genetic basis for its phenotypic diversity. To provide enriched genomic resources and infer possible mechanisms underlying skin pigmentation, we performed a large-scale transcriptomic sequencing on 13 adult goldfish tissues, larvae at one- and three-days post hatch, and skin tissues with four different color pigmentation. A total of 25.52 Gb and 149.80 Gb clean data were obtained using the PacBio and Illumina platforms, respectively. Onto the goldfish reference genome, we mapped 137,674 non-redundant transcripts, of which 5.54% was known isoforms and 78.53% was novel isoforms of the reference genes, and the remaining 21,926 isoforms are novel isoforms of additional new genes. Both skin-specific and color-specific transcriptomic analyses showed that several significantly enriched genes were known to be involved in melanogenesis, tyrosine metabolism, PPAR signaling pathway, folate biosynthesis metabolism and so on. Thirteen differentially expressed genes across different color skins were associated with melanogenesis and pteridine synthesis including mitf, ednrb, mc1r, tyr, mlph and gch1, and xanthophore differentiation such as pax7, slc2a11 and slc2a15. These transcriptomic data revealed pathways involved in goldfish pigmentation and improved the gene annotation of the reference genome.

Upregulation of the PPAR signaling pathway and accumulation of lipids are related to the morphological and structural transformation of the dragon-eye goldfish eye

Goldfish comprise around 300 different strains with drastically altered and aesthetical morphologies making them suitable models for evolutionary developmental biology. The dragon-eye strain is characterized by protruding eyes (analogous to those of Chinese dragons). Although the strain has been selected for about 400 years, the mechanism of its eye development remains unclear. In this study, a stable dragon-eye goldfish strain with a clear genetic background was rapidly established and studied. We found that upregulation of the PPAR signaling pathway accompanied by an increase in lipid accumulation might trigger the morphological and structural transformation of the eye in dragon-eye goldfish. At the developmental stage of proptosis (eye protrusion), downregulation of the phototransduction pathway was consistent with the structural defects and myopia of the dragon-eye strain. With the impairment of retinal development, cytokine-induced inflammation was activated, especially after proptosis, similar to the pathologic symptoms of many human ocular diseases. In addition, differentially expressed transcription factors were significantly enriched in the PAX and homeobox families, two well-known transcription factor families involved in eye development. Therefore, our findings reveal the dynamic changes in key pathways during eye development in dragon-eye goldfish, and provide insights into the molecular mechanisms underlying drastically altered eyes in goldfish and human ocular disease.

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

Goldfish have been subjected to over 1,000 y of intensive domestication and selective breeding. In this report, we describe a high-quality goldfish genome (2n = 100), anchoring 95.75% of contigs into 50 pseudochromosomes. Comparative genomics enabled us to disentangle the two subgenomes that resulted from an ancient hybridization event. Resequencing 185 representative goldfish variants and 16 wild crucian carp revealed the origin of goldfish and identified genomic regions that have been shaped by selective sweeps linked to its domestication. Our comprehensive collection of goldfish varieties enabled us to associate genetic variations with a number of well-known anatomical features, including features that distinguish traditional goldfish clades. Additionally, we identified a tyrosine-protein kinase receptor as a candidate causal gene for the first well-known case of Mendelian inheritance in goldfish-the transparent mutant. The goldfish genome and diversity data offer unique resources to make goldfish a promising model for functional genomics, as well as domestication.

Excitation-Contraction Coupling in the Goldfish (Carassius auratus) Intact Heart

Cardiac physiology of fish models is an emerging field given the ease of genome editing and the development of transgenic models. Several studies have described the cardiac properties of zebrafish (Denio rerio). The goldfish (Carassius auratus) belongs to the same family as the zebrafish and has emerged as an alternative model with which to study cardiac function. Here, we propose to acutely study electrophysiological and systolic Ca²⁺ signaling in intact goldfish hearts. We assessed the Ca²⁺ dynamics and the electrophysiological cardiac function of goldfish, zebrafish, and mice models, using pulsed local field fluorescence microscopy, intracellular microelectrodes, and flash photolysis in perfused hearts. We observed goldfish ventricular action potentials (APs) and Ca²⁺ transients to be significantly longer when compared to the zebrafish. The action potential half duration at 50% (APD₅₀) of goldfish was 370.38 ± 8.8 ms long, and in the zebrafish they were observed to be only 83.9 ± 9.4 ms. Additionally, the half duration of the Ca²⁺ transients was also longer for goldfish (402.1 ± 4.4 ms) compared to the zebrafish (99.1 ± 2.7 ms). Also, blocking of the L-type Ca²⁺ channels with nifedipine revealed this current has a major role in defining the amplitude and the duration of goldfish Ca²⁺ transients. Interestingly, nifedipine flash photolysis experiments in the intact heart identified whether or not the decrease in the amplitude of Ca²⁺ transients was due to shorter APs. Moreover, an increase in temperature and heart rate had a strong shortening effect on the AP and Ca²⁺ transients of goldfish hearts. Furthermore, ryanodine (Ry) and thapsigargin (Tg) significantly reduced the amplitude of the Ca²⁺ transients, induced a prolongation in the APs, and altogether exhibited the degree to which the Ca²⁺ release from the sarcoplasmic reticulum contributed to the Ca²⁺ transients. We conclude that the electrophysiological properties and Ca²⁺ signaling in intact goldfish hearts strongly resembles the endocardial layer of larger mammals.