Three genomes in the algal genus Volvox reveal the fate of a haploid sex-determining region after a transition to homothallism.

Transitions between separate sexes (dioecy) and other mating systems are common across eukaryotes. Here, we study a change in a haploid dioecious green algal species with male- and female-determining chromosomes (U and V). The genus Volvox is an oogamous (with large, immotile female gametes and small, motile male gametes) and includes both heterothallic species (with distinct male and female genotypes, associated with a mating-type system that prevents fusion of gametes of the same sex) and homothallic species (bisexual, with the ability to self-fertilize). We date the origin of an expanded sex-determining region (SDR) in Volvox to at least 75 Mya, suggesting that homothallism represents a breakdown of dioecy (heterothallism). We investigated the involvement of the SDR of the U and V chromosomes in this transition. Using de novo whole-genome sequences, we identified a heteromorphic SDR of ca 1 Mbp in male and female genotypes of the heterothallic species Volvox reticuliferus and a homologous region (SDLR) in the closely related homothallic species Volvox africanus, which retained several different hallmark features of an SDR. The V. africanus SDLR includes a large region resembling the female SDR of the presumptive heterothallic ancestor, whereas most genes from the male SDR are absent. However, we found a multicopy array of the male-determining gene, MID, in a different genomic location from the SDLR. Thus, in V. africanus, an ancestrally female genotype may have acquired MID and thereby gained male traits.

Genomic and Secretomic Analyses of the Newly Isolated Fungus Perenniporia fraxinea SS3 Identified CAZymes Potentially Related to a Serious Pathogenesis of Hardwood Trees.

Perenniporia fraxinea can colonize living trees and cause severe damage to standing hardwoods by secreting a number of carbohydrate-activate enzymes (CAZymes), unlike other well-studied Polyporales. However, significant knowledge gaps exist in understanding the detailed mechanisms for this hardwood-pathogenic fungus. To address this issue, five monokaryotic P. fraxinea strains, SS1 to SS5, were isolated from the tree species Robinia pseudoacacia, and high polysaccharide-degrading activities and the fastest growth were found for P. fraxinea SS3 among the isolates. The whole genome of P. fraxinea SS3 was sequenced, and its unique CAZyme potential for tree pathogenicity was determined in comparison to the genomes of other nonpathogenic Polyporales. These CAZyme features are well conserved in a distantly related tree pathogen, Heterobasidion annosum. Furthermore, the carbon source-dependent CAZyme secretions of P. fraxinea SS3 and a nonpathogenic and strong white-rot Polyporales member, Phanerochaete chrysosporium RP78, were compared by activity measurements and proteomic analyses. As seen in the genome comparisons, P. fraxinea SS3 exhibited higher pectin-degrading activities and higher laccase activities than P. chrysosporium RP78, which were attributed to the secretion of abundant glycoside hydrolase family 28 (GH28) pectinases and auxiliary activity family 1_1 (AA1_1) laccases, respectively. These enzymes are possibly related to fungal invasion into the tree lumens and the detoxification of tree defense substances. Additionally, P. fraxinea SS3 showed secondary cell wall degradation capabilities at the same level as that of P. chrysosporium RP78. Overall, this study suggested mechanisms for how this fungus can attack the cell walls of living trees as a serious pathogen and differs from other nonpathogenic white-rot fungi. IMPORTANCE Many studies have been done to understand the mechanisms underlying the degradation of plant cell walls of dead trees by wood decay fungi. However, little is known about how some of these fungi weaken living trees as pathogens. P. fraxinea belongs to the Polyporales, a group of strong wood decayers, and is known to aggressively attack and fell standing hardwood trees all over the world. Here, we report CAZymes potentially related to plant cell wall degradation and pathogenesis factors in a newly isolated fungus, P. fraxinea SS3, by genome sequencing in conjunction with comparative genomic and secretomic analyses. The present study provides insights into the mechanisms of the degradation of standing hardwood trees by the tree pathogen, which will contribute to the prevention of this serious tree disease.

MetaPlatanus: a metagenome assembler that combines long-range sequence links and species-specific features.

De novo metagenome assembly is effective in assembling multiple draft genomes, including those of uncultured organisms. However, heterogeneity in the metagenome hinders assembly and introduces interspecies misassembly deleterious for downstream analysis. For this purpose, we developed a hybrid metagenome assembler, MetaPlatanus. First, as a characteristic function, it assembles the basic contigs from accurate short reads and then iteratively utilizes long-range sequence links, species-specific sequence compositions, and coverage depth. The binning information was also used to improve contiguity. Benchmarking using mock datasets consisting of known bacteria with long reads or mate pairs revealed the high contiguity MetaPlatanus with a few interspecies misassemblies. For published human gut data with nanopore reads from potable sequencers, MetaPlatanus assembled many biologically important elements, such as coding genes, gene clusters, viral sequences, and over-half bacterial genomes. In the benchmark with published human saliva data with high-throughput nanopore reads, the superiority of MetaPlatanus was considerably more evident. We found that some high-abundance bacterial genomes were assembled only by MetaPlatanus as near-complete. Furthermore, MetaPlatanus can circumvent the limitations of highly fragmented assemblies and frequent interspecies misassembles obtained by the other tools. Overall, the study demonstrates that MetaPlatanus could be an effective approach for exploring large-scale structures in metagenomes.

MoG+: a database of genomic variations across three mouse subspecies for biomedical research.

Laboratory mouse strains have mosaic genomes derived from at least three major subspecies that are distributed in Eurasia. Here, we describe genomic variations in ten inbred strains: Mus musculus musculus-derived BLG2/Ms, NJL/Ms, CHD/Ms, SWN/Ms, and KJR/Ms; M. m. domesticus-derived PGN2/Ms and BFM/Ms; M. m. castaneus-derived HMI/Ms; and JF1/Ms and MSM/Ms, which were derived from a hybrid between M. m. musculus and M. m. castaneus. These strains were established by Prof. Moriwaki in the 1980s and are collectively named the "Mishima Battery". These strains show large phenotypic variations in body size and in many physiological traits. We resequenced the genomes of the Mishima Battery strains and performed a comparative genomic analysis with dbSNP data. More than 81 million nucleotide coordinates were identified as variant sites due to the large genetic distances among the mouse subspecies; 8,062,070 new SNP sites were detected in this study, and these may underlie the large phenotypic diversity observed in the Mishima Battery. The new information was collected in a reconstructed genome database, termed MoG+ that includes new application software and viewers. MoG+ intuitively visualizes nucleotide variants in genes and intergenic regions, and amino acid substitutions across the three mouse subspecies. We report statistical data from the resequencing and comparative genomic analyses and newly collected phenotype data of the Mishima Battery, and provide a brief description of the functions of MoG+, which provides a searchable and unique data resource of the numerous genomic variations across the three mouse subspecies. The data in MoG+ will be invaluable for research into phenotype-genotype links in diverse mouse strains.

Efficient de novo assembly of highly heterozygous genomes from whole-genome shotgun short reads.

Although many de novo genome assembly projects have recently been conducted using high-throughput sequencers, assembling highly heterozygous diploid genomes is a substantial challenge due to the increased complexity of the de Bruijn graph structure predominantly used. To address the increasing demand for sequencing of nonmodel and/or wild-type samples, in most cases inbred lines or fosmid-based hierarchical sequencing methods are used to overcome such problems. However, these methods are costly and time consuming, forfeiting the advantages of massive parallel sequencing. Here, we describe a novel de novo assembler, Platanus, that can effectively manage high-throughput data from heterozygous samples. Platanus assembles DNA fragments (reads) into contigs by constructing de Bruijn graphs with automatically optimized k-mer sizes followed by the scaffolding of contigs based on paired-end information. The complicated graph structures that result from the heterozygosity are simplified during not only the contig assembly step but also the scaffolding step. We evaluated the assembly results on eukaryotic samples with various levels of heterozygosity. Compared with other assemblers, Platanus yields assembly results that have a larger scaffold NG50 length without any accompanying loss of accuracy in both simulated and real data. In addition, Platanus recorded the largest scaffold NG50 values for two of the three low-heterozygosity species used in the de novo assembly contest, Assemblathon 2. Platanus therefore provides a novel and efficient approach for the assembly of gigabase-sized highly heterozygous genomes and is an attractive alternative to the existing assemblers designed for genomes of lower heterozygosity.

Chromosomal Aneuploidy Improves the Brewing Characteristics of Sake Yeast.

The effect of chromosomal aneuploidy on the brewing characteristics of brewery yeasts has not been studied. Here we report that chromosomal aneuploidy in sake brewery yeast (Saccharomyces cerevisiae) leads to the development of favorable brewing characteristics. We found that pyruvate-underproducing sake yeast, which produces less off-flavor diacetyl, is aneuploid and trisomic for chromosomes XI and XIV. To confirm that this phenotype is due to aneuploidy, we obtained 45 haploids with various chromosomal additions and investigated their brewing profiles. A greater number of chromosomes correlated with a decrease in pyruvate production. Especially, sake yeast haploids with extra chromosomes in addition to chromosome XI produced less pyruvate than euploids. Mitochondrion-related metabolites and intracellular oxygen species in chromosome XI aneuploids were higher than those in euploids, and this effect was canceled in their "petite" strains, suggesting that an increase in chromosomes upregulated mitochondrial activity and decreased pyruvate levels. These findings suggested that an increase in chromosome number, including chromosome XI, in sake yeast haploids leads to pyruvate underproduction through the augmentation of mitochondrial activity. This is the first report proposing that aneuploidy in brewery yeasts improves their brewing profile.IMPORTANCE Chromosomal aneuploidy has not been evaluated in development of sake brewing yeast strains. This study shows the relationship between chromosomal aneuploidy and brewing characteristics of brewery yeast strains. High concentrations of pyruvate during sake storage give rise to α-acetolactate and, in turn, to high concentrations of diacetyl, which is considered an off-flavor. It was demonstrated that pyruvate-underproducing sake yeast is trisomic for chromosome XI and XIV. Furthermore, sake yeast haploids with extra chromosomes produced reduced levels of pyruvate and showed metabolic processes characteristic of increased mitochondrial activity. This novel discovery will enable the selection of favorable brewery yeasts by monitoring the copy numbers of specific chromosomes through a process that does not involve generation/use of genetically modified organisms.

The ancestor of extant Japanese fancy mice contributed to the mosaic genomes of classical inbred strains.

Commonly used classical inbred mouse strains have mosaic genomes with sequences from different subspecific origins. Their genomes are derived predominantly from the Western European subspecies Mus musculus domesticus, with the remaining sequences derived mostly from the Japanese subspecies Mus musculus molossinus. However, it remains unknown how this intersubspecific genome introgression occurred during the establishment of classical inbred strains. In this study, we resequenced the genomes of two M. m. molossinus-derived inbred strains, MSM/Ms and JF1/Ms. MSM/Ms originated from Japanese wild mice, and the ancestry of JF1/Ms was originally found in Europe and then transferred to Japan. We compared the characteristics of these sequences to those of the C57BL/6J reference sequence and the recent data sets from the resequencing of 17 inbred strains in the Mouse Genome Project (MGP), and the results unequivocally show that genome introgression from M. m. molossinus into M. m. domesticus provided the primary framework for the mosaic genomes of classical inbred strains. Furthermore, the genomes of C57BL/6J and other classical inbred strains have long consecutive segments with extremely high similarity (>99.998%) to the JF1/Ms strain. In the early 20th century, Japanese waltzing mice with a morphological phenotype resembling that of JF1/Ms mice were often crossed with European fancy mice for early studies of "Mendelism," which suggests that the ancestor of the extant JF1/Ms strain provided the origin of the M. m. molossinus genome in classical inbred strains and largely contributed to its intersubspecific genome diversity.

Coelacanth genomes reveal signatures for evolutionary transition from water to land.

Coelacanths are known as "living fossils," as they show remarkable morphological resemblance to the fossil record and belong to the most primitive lineage of living Sarcopterygii (lobe-finned fishes and tetrapods). Coelacanths may be key to elucidating the tempo and mode of evolution from fish to tetrapods. Here, we report the genome sequences of five coelacanths, including four Latimeria chalumnae individuals (three specimens from Tanzania and one from Comoros) and one L. menadoensis individual from Indonesia. These sequences cover two African breeding populations and two known extant coelacanth species. The genome is ∼2.74 Gbp and contains a high proportion (∼60%) of repetitive elements. The genetic diversity among the individuals was extremely low, suggesting a small population size and/or a slow rate of evolution. We found a substantial number of genes that encode olfactory and pheromone receptors with features characteristic of tetrapod receptors for the detection of airborne ligands. We also found that limb enhancers of bmp7 and gli3, both of which are essential for limb formation, are conserved between coelacanth and tetrapods, but not ray-finned fishes. We expect that some tetrapod-like genes may have existed early in the evolution of primitive Sarcopterygii and were later co-opted to adapt to terrestrial environments. These coelacanth genomes will provide a cornerstone for studies to elucidate how ancestral aquatic vertebrates evolved into terrestrial animals.

Comparative genome and transcriptome analyses of the social amoeba Acytostelium subglobosum that accomplishes multicellular development without germ-soma differentiation.

BACKGROUND: Social amoebae are lower eukaryotes that inhabit the soil. They are characterized by the construction of a starvation-induced multicellular fruiting body with a spore ball and supportive stalk. In most species, the stalk is filled with motile stalk cells, as represented by the model organism Dictyostelium discoideum, whose developmental mechanisms have been well characterized. However, in the genus Acytostelium, the stalk is acellular and all aggregated cells become spores. Phylogenetic analyses have shown that it is not an ancestral genus but has lost the ability to undergo cell differentiation. RESULTS: We performed genome and transcriptome analyses of Acytostelium subglobosum and compared our findings to other available dictyostelid genome data. Although A. subglobosum adopts a qualitatively different developmental program from other dictyostelids, its gene repertoire was largely conserved. Yet, families of polyketide synthase and extracellular matrix proteins have not expanded and a serine protease and ABC transporter B family gene, tagA, and a few other developmental genes are missing in the A. subglobosum lineage. Temporal gene expression patterns are astonishingly dissimilar from those of D. discoideum, and only a limited fraction of the ortholog pairs shared the same expression patterns, so that some signaling cascades for development seem to be disabled in A. subglobosum. CONCLUSIONS: The absence of the ability to undergo cell differentiation in Acytostelium is accompanied by a small change in coding potential and extensive alterations in gene expression patterns.

Comparative analysis of chimpanzee and human Y chromosomes unveils complex evolutionary pathway.

The mammalian Y chromosome has unique characteristics compared with the autosomes or X chromosomes. Here we report the finished sequence of the chimpanzee Y chromosome (PTRY), including 271 kb of the Y-specific pseudoautosomal region 1 and 12.7 Mb of the male-specific region of the Y chromosome. Greater sequence divergence between the human Y chromosome (HSAY) and PTRY (1.78%) than between their respective whole genomes (1.23%) confirmed the accelerated evolutionary rate of the Y chromosome. Each of the 19 PTRY protein-coding genes analyzed had at least one nonsynonymous substitution, and 11 genes had higher nonsynonymous substitution rates than synonymous ones, suggesting relaxation of selective constraint, positive selection or both. We also identified lineage-specific changes, including deletion of a 200-kb fragment from the pericentromeric region of HSAY, expansion of young Alu families in HSAY and accumulation of young L1 elements and long terminal repeat retrotransposons in PTRY. Reconstruction of the common ancestral Y chromosome reflects the dynamic changes in our genomes in the 5-6 million years since speciation.

Transcriptomic comparison between populations selected for higher and lower mobility in the red flour beetle Tribolium castaneum.

Movement is an important behavior observed in a wide range of taxa. Previous studies have examined genes controlling movement using wing polymorphic insects and genes controlling wing size. However, few studies have investigated genes controlling movement activity rather than morphological traits. In the present study, we conducted RNA sequencing using populations with higher (WL) and lower (WS) mobility established by artificial selection in the red flour beetle Tribolium castaneum and compared gene expression levels between selected populations with two replicate lines. As a result, we found significant differences between the selected populations in 677 genes expressed in one replicate line and 1198 genes expressed in another replicate line, of which 311 genes were common to the two replicate lines. Furthermore, quantitative PCR focusing on 6 of these genes revealed that neuropeptide F receptor gene (NpF) was significantly more highly expressed in the WL population than in the WS population, which was common to the two replicate lines. We discuss differences in genes controlling movement between walking activity and wing polymorphism.

Genome assembly of Genji firefly (Nipponoluciola cruciata) reveals novel luciferase-like luminescent proteins without peroxisome targeting signal.

The Genji firefly, Nipponoluciola cruciata, is an aquatic firefly endemic to Japan, inhabiting a wide area of the Japanese archipelago. The luminescence of fireflies is a scientifically interesting phenomenon, and many studies have evaluated this species in Japan. In this study, we sequenced the whole genome of male N. cruciata and constructed a high-quality genome assembly of 662 Mb with a BUSCO completeness of 99.1% in the genome mode. Using the detected set of 15,169 protein-coding genes, the genomic structures and genetic background of luminescence-related genes were also investigated. We found four new firefly luciferase-like genes in the genome. The highest bioluminescent activity was observed for LLa2, which originated from ancestral PDGY, a mitochondrial acyl-CoA synthetase. A thioesterase candidate, NcruACOT1, which is involved in d-luciferin biosynthesis, was expressed in the lantern. Two opsins were also detected and the absorption wavelength of the UV-type opsin candidate shifted from UV to blue. These findings provide an important resource for unravelling the adaptive evolution of fireflies in terms of luminescence and vision.

BRD9 determines the cell fate of hematopoietic stem cells by regulating chromatin state.

ATP-dependent chromatin remodeling SWI/SNF complexes exist in three subcomplexes: canonical BAF (cBAF), polybromo BAF (PBAF), and a newly described non-canonical BAF (ncBAF). While cBAF and PBAF regulate fates of multiple cell types, roles for ncBAF in hematopoietic stem cells (HSCs) have not been investigated. Motivated by recent discovery of disrupted expression of BRD9, an essential component of ncBAF, in multiple cancers, including clonal hematopoietic disorders, we evaluate here the role of BRD9 in normal and malignant HSCs. BRD9 loss enhances chromatin accessibility, promoting myeloid lineage skewing while impairing B cell development. BRD9 significantly colocalizes with CTCF, whose chromatin recruitment is augmented by BRD9 loss, leading to altered chromatin state and expression of myeloid-related genes within intact topologically associating domains. These data uncover ncBAF as critical for cell fate specification in HSCs via three-dimensional regulation of gene expression and illuminate roles for ncBAF in normal and malignant hematopoiesis.

Human chromosome 11 DNA sequence and analysis including novel gene identification.

Chromosome 11, although average in size, is one of the most gene- and disease-rich chromosomes in the human genome. Initial gene annotation indicates an average gene density of 11.6 genes per megabase, including 1,524 protein-coding genes, some of which were identified using novel methods, and 765 pseudogenes. One-quarter of the protein-coding genes shows overlap with other genes. Of the 856 olfactory receptor genes in the human genome, more than 40% are located in 28 single- and multi-gene clusters along this chromosome. Out of the 171 disorders currently attributed to the chromosome, 86 remain for which the underlying molecular basis is not yet known, including several mendelian traits, cancer and susceptibility loci. The high-quality data presented here--nearly 134.5 million base pairs representing 99.8% coverage of the euchromatic sequence--provide scientists with a solid foundation for understanding the genetic basis of these disorders and other biological phenomena.

Profiling a single-stranded DNA region within an rDNA segment that affects the loading of bacterial condensin.

Bacterial condensin preferentially loads onto single-stranded DNA (ssDNA) in vitro and onto rDNA in vivo to support proper chromosome compaction. Thus, the actively transcribing rDNA would provide the ssDNA region for the loading of bacterial condensin. We attempted to detect the ssDNA region in the rrnI gene in situ. Non-denaturing sodium bisulfite treatment catalyzed the conversion of cytosines to thymines (CT-conversion) at the melted DNA of a genome. Using next-generation sequencing, we generated an average of 11,000 reads covering each cytosine on the rDNA segment. In principle, the CT-conversion rate is an accurate guide to detect ssDNA segment. We detected multiple ssDNA segments throughout the rDNA. The deletion mutations of the rDNA hindered the ssDNA formation at the 100-500 bp segment downstream of the promoter. These data support the idea that the ssDNA segment plays a crucial role as the condensin-loading site and suggest the mechanism of condensin loading onto rDNA.

Relationship Between the Structure of Alkylpsoralens and Their Antiproliferative Activity in HL60 Cells.

BACKGROUND/AIM: As part of our continuing investigation in coumarin derivatives as potential anticancer substances, a series of alkylpsoralens were synthesized, and their antiproliferative activity was evaluated in leukemic HL60 cells. MATERIALS AND METHODS: Alkylpsoralens were systematically synthesized from the combination of several chloroketones and 7-hydroxycoumarin derivatives. RESULTS: Among the compounds synthesized, 4,4',8-trimethylpsoralen demonstrated the most potent activity (IC50=6.6 µM). CONCLUSION: The correlation between the alkylation pattern and antiproliferative activity showed the importance of the C4-methyl and C8-methyl moieties in the psoralen nucleus as well as the importance of lipophilicity for their antiproliferative activity.

Complete Genome Sequences of Three Phocaeicola vulgatus Strains Isolated from a Healthy Japanese Individual.

Phocaeicola vulgatus (formerly Bacteroides vulgatus) is a pathogenic anaerobic bacterium frequently involved in human infections. We present the complete genome sequences of three Phocaeicola vulgatus strains isolated from the same healthy person, determined by hybrid assembly using Nanopore long-read sequencing and DNBseq short-read sequencing.

Horizontal transfer of a plasmid possessing mcr-1 marked with a single nucleotide mutation between Escherichia coli isolates from community residents.

OBJECTIVES: The widespread dissemination of phenotypic colistin-resistant (COR) bacteria in the community threatens public health. The horizontal gene transfer of the mobile colistin resistance gene via plasmids is thought to be one of the main mechanisms for dissemination. However, genotypic evidence to prove this in community settings is limited. This study used genome analysis to demonstrate the direct horizontal colistin resistance gene transfer via plasmids in isolates from the community. RESULTS: A total of 19 isolates of COR Escherichia coli from stool specimens of 23 residents from seven households in the Vietnamese community were assessed in this study. The whole-genome sequence data of isolates were acquired using a combination of DNBSEQ short-reads and Nanopore long-read sequencing. Analysis of genomic data was performed using online tools such as Geneious. Analysis of the genomic information of COR E. coli isolates revealed that the isolates from two residents of different households had a similar IncP1 plasmid possessing mcr-1.1, marked with a single nucleotide mutation at the same position. The study provided direct evidence to prove that mcr was horizontally transmitted among bacteria in community residents.

Nintedanib induces gene expression changes in the lung of induced-rheumatoid arthritis-associated interstitial lung disease mice.

Nintedanib is a multi-tyrosine kinase inhibitor widely used to treat progressive fibrosing interstitial lung diseases because it slows the reduction in forced vital capacity. However, the prognosis for patients treated with nintedanib remains poor. To improve nintedanib treatment, we examined the effects of nintedanib on gene expression in the lungs of induced-rheumatoid arthritis-associated interstitial lung disease model mice, which develop rheumatoid arthritis and subsequent pulmonary fibrosis. Using next-generation sequencing, we identified 27 upregulated and 130 downregulated genes in the lungs of these mice after treatment with nintedanib. The differentially expressed genes included mucin 5B and heat shock protein 70 family genes, which are related to interstitial lung diseases, as well as genes associated with extracellular components, particularly the myocardial architecture, suggesting unanticipated effects of nintedanib. Of the genes upregulated in the nintedanib-treated lung, expression of regulatory factor X2, which is suspected to be involved in cilia movement, and bone morphogenetic protein receptor type 2, which is involved in the pathology of pulmonary hypertension, was detected by immunohistochemistry and RNA in situ hybridization in peripheral airway epithelium and alveolar cells. Thus, the present findings indicate a set of genes whose expression alteration potentially underlies the effects of nintedanib on pulmonary fibrosis. It is expected that these findings will contribute to the development of improved nintedanib strategies for the treatment of progressive fibrosing interstitial lung diseases.

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.