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

Expansion and Functional Diversification of Long-Wavelength-Sensitive Opsin in Anabantoid Fishes

Gene duplication is one of the most important sources of novel genotypic diversity and the subsequent evolution of phenotypic diversity. Determining the evolutionary history and functional changes of duplicated genes is crucial for a comprehensive understanding of adaptive evolution. The evolutionary history of visual opsin genes is very dynamic, with repeated duplication events followed by sub- or neofunctionalization. While duplication of the green-sensitive opsins rh2 is common in teleost fish, fewer cases of multiple duplication events of the red-sensitive opsin lws are known. In this study, we investigate the visual opsin gene repertoire of the anabantoid fishes, focusing on the five lws opsin genes found in the genus Betta. We determine the evolutionary history of the lws opsin gene by taking advantage of whole-genome sequences of nine anabantoid species, including the newly assembled genome of Betta imbellis. Our results show that at least two independent duplications of lws occurred in the Betta lineage. The analysis of amino acid sequences of the lws paralogs of Betta revealed high levels of diversification in four of the seven transmembrane regions of the lws protein. Amino acid substitutions at two key-tuning sites are predicted to lead to differentiation of absorption maxima (λmax) between the paralogs within Betta. Finally, eye transcriptomics of B. splendens at different developmental stages revealed expression shifts between paralogs for all cone opsin classes. The lws genes are expressed according to their relative position in the lws opsin cluster throughout ontogeny. We conclude that temporal collinearity of lws expression might have facilitated subfunctionalization of lws in Betta and teleost opsins in general.

Identification of pigmentation genes in skin, muscle and tail of a Thai-flag variety of Siamese fighting fish Betta splendens

Pigmentation genes expressed in skin, body muscle and tail of Thai-flag compared with Blue, White and Red varieties of Siamese fighting fish Betta splendens were identified. In total, 22,919 new unigenes were found. Pearson correlation and PCA analysis revealed that expression profiles of genes in muscle, skin and tail across solid color variety were similar. In contrast, those in skin and red tail part of Thai-flag were closely related but they showed different expression profiles with the white tail part. Moreover, 21,347-64,965 SNPs were identified in exonic regions of identified genes. In total, 28,899 genes were differentially expressed between paired comparisons of libraries where 13,907 genes (48.12 %) were upregulated and 14,992 genes (51.88 %) were downregulated. DEGs between paired libraries were 106-5775 genes relative to the compared libraries (56-2982 and 50-2782 for upregulated and downregulated DEGs). Interestingly, 432 pigmentation genes of B. splendens were found. Of these, 297 DEGs showed differential expression between varieties. Many DEGs in melanogenesis (Bsmcr1r, Bsmcr5r, and Bsslc2a15b), tyrosine metabolism (Bstyr, Bstyrp1b and Bsdct), stripe repressor (BsAsip1 and BsAsip2b), pteridine (Bsgch2) and carotenoid (BsBco2) biosynthesis were downregulated in the Thai-flag compared with solid color varieties. Expression of Bsbco1l, Bsfrem2b, Bskcnj13, Bszic2a and Bspah in skin, muscle and tail of Thai-flag, Blue, Red and White varieties was analyzed by qRT-PCR and revealed differential expression between fish varieties and showed anatomical tissue-preferred expression patterns in the same fish variety. The information could be applied to assist genetic-based development of new B. splendens varieties in the future.

The reference genome of Macropodus opercularis (the paradise fish)

Amongst fishes, zebrafish (Danio rerio) has gained popularity as a model system over most other species and while their value as a model is well documented, their usefulness is limited in certain fields of research such as behavior. By embracing other, less conventional experimental organisms, opportunities arise to gain broader insights into evolution and development, as well as studying behavioral aspects not available in current popular model systems. The anabantoid paradise fish (Macropodus opercularis), an "air-breather" species has a highly complex behavioral repertoire and has been the subject of many ethological investigations but lacks genomic resources. Here we report the reference genome assembly of M. opercularis using long-read sequences at 150-fold coverage. The final assembly consisted of 483,077,705 base pairs (~483 Mb) on 152 contigs. Within the assembled genome we identified and annotated 20,157 protein coding genes and assigned ~90% of them to orthogroups.

The paradise fish, an advanced animal model for behavioral genetics and evolutionary developmental biology

Paradise fish (Macropodus opercularis) is an air-breathing freshwater fish species with a signature labyrinth organ capable of extracting oxygen from the air that helps these fish to survive in hypoxic environments. The appearance of this evolutionary innovation in anabantoids resulted in a rewired circulatory system, but also in the emergence of species-specific behaviors, such as territorial display, courtship and parental care in the case of the paradise fish. Early zoologists were intrigued by the structure and function of the labyrinth apparatus and a series of detailed descriptive histological studies at the beginning of the 20th century revealed the ontogenesis and function of this specialized system. A few decades later, these fish became the subject of numerous ethological studies, and detailed ethograms of their behavior were constructed. These latter studies also demonstrated a strong genetic component underlying their behavior, but due to lack of adequate molecular tools, the fine genetic dissection of the behavior was not possible at the time. The technological breakthroughs that transformed developmental biology and behavioral genetics in the past decades, however, give us now a unique opportunity to revisit these old questions. Building on the classic descriptive studies, the new methodologies will allow us to follow the development of the labyrinth apparatus at a cellular resolution, reveal the genes involved in this process and also the genetic architecture behind the complex behaviors that we can observe in this species.

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.

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.

Genetic manipulation of betta fish

Betta splendens, also known as Siamese fighting fish or "betta," is a freshwater fish species renowned for its astonishing morphological diversity and extreme aggressive behavior. Despite recent advances in our understanding of the genetics and neurobiology of betta, the lack of tools to manipulate their genome has hindered progress at functional and mechanistic levels. In this study, we outline the use of three genetic manipulation technologies, which we have optimized for use in betta: CRISPR/Cas9-mediated knockout, CRISPR/Cas9-mediated knockin, and Tol2-mediated transgenesis. We knocked out three genes: alkal2l, bco1l, and mitfa, and analyzed their effects on viability and pigmentation. Furthermore, we knocked in a fluorescent protein into the mitfa locus, a proof-of-principle experiment of this powerful technology in betta. Finally, we used Tol2-mediated transgenesis to create fish with ubiquitous expression of GFP, and then developed a bicistronic plasmid with heart-specific expression of a red fluorescent protein to serve as a visible marker of successful transgenesis. Our work highlights the potential for the genetic manipulation of betta, providing valuable resources for the effective use of genetic tools in this animal model.

Combination of GWAS and FST-based approaches identified loci associated with economic traits in sugarcane

Sugarcane is a globally important plant for both sugar and biofuel production. Although conventional breeding has played an important role in increasing the productivity of sugarcane, it takes a long time to achieve breeding goals such as high yield and resistant to diseases. Molecular breeding, including marker-assisted breeding and genomic selection, can accelerate genetic improvement by selecting elites at the seedling stage with DNA markers. However, only a few DNA markers associated with important traits were identified in sugarcane. The purpose of this study was to identify DNA markers associated with sugar content, stalk diameter, and sugarcane top borer resistance. The sugarcane samples with trait records were genotyped using the restriction site-associated DNA sequencing (RADseq) technology. Using F_ST analysis and genome-wide association study (GWAS), a total of 9, 23 and 9 DNA variants (single nucleotide polymorphisms (SNPs)/insertions and deletions (indels)) were associated with sugar content, stalk diameter, and sugarcane top borer resistance, respectively. The identified genetic variants were on different chromosomes, suggesting that these traits are complex and determined by multiple genetic factors. These DNA markers identified by both approaches have the potential to be used in selecting elite clones at the seeding stage in our sugarcane breeding program to accelerate genetic improvement. Certainly, it is essential to verify the reliability of the identified DNA markers associated with traits before they are used in molecular breeding in other populations.

Transcriptome analyses of betta fish (Betta splendens) provide novel insights into fin regeneration and color-related genes

The betta fish (Betta splendens), an important ornamental fish, haswell-developed and colorful fins.After fin amputation, betta fish can easily regenerate finssimilar to the originalsin terms of structureand color. The powerful fin regeneration ability and a variety of colors in the betta fish are fascinating. However, the underlying molecular mechanisms are still not fully understood. In this study, tail fin amputation and regeneration experiments were performed on two kinds of betta fish: red and white color betta fish. Then, transcriptome analyseswere conducted to screen out fin regeneration and color-relatedgenes in betta fish. Through enrichment analyses of differentially expressed genes (DEGs), we founda series of enrichment pathways and genes related to finregeneration, including cell cycle (i.e. plcg2), TGF-beta signaling pathway (i.e. bmp6), PI3K-Akt signaling pathway (i.e. loxl2aand loxl2b), Wnt signaling pathway(i.e. lef1), gap junctions (i.e. cx43), angiogenesis (i.e. foxp1), and interferon regulatory factor (i.e. irf8). Meanwhile, some fin color-related pathways and genes were identified in betta fish, especially melanogenesis (i.e. tyr, tyrp1a, tyrp1b, and mc1r) and carotenoid color genes (i.e. pax3, pax7, sox10, and ednrba). In conclusion, this studycan not only enrich the research onfish tissue regeneration, but also has a potential significance for the aquaculture and breeding of the betta fish.

The GRN concept as a guide for evolutionary developmental biology

Organismal phenotypes result largely from inherited developmental programs, usually executed during embryonic and juvenile life stages. These programs are not blank slates onto which natural selection can draw arbitrary forms. Rather, the mechanisms of development play an integral role in shaping phenotypic diversity and help determine the evolutionary trajectories of species. Modern evolutionary biology must, therefore, account for these mechanisms in both theory and in practice. The gene regulatory network (GRN) concept represents a potent tool for achieving this goal whose utility has grown in tandem with advances in "omic" technologies and experimental techniques. However, while the GRN concept is widely utilized, it is often less clear what practical implications it has for conducting research in evolutionary developmental biology. In this Perspective, we attempt to provide clarity by discussing how experiments and projects can be designed in light of the GRN concept. We first map familiar biological notions onto the more abstract components of GRN models. We then review how diverse functional genomic approaches can be directed toward the goal of constructing such models and discuss current methods for functionally testing evolutionary hypotheses that arise from them. Finally, we show how the major steps of GRN model construction and experimental validation suggest generalizable workflows that can serve as a scaffold for project design. Taken together, the practical implications that we draw from the GRN concept provide a set of guideposts for studies aiming at unraveling the molecular basis of phenotypic diversity.

Genetic manipulation of betta fish

Betta splendens , also known as Siamese fighting fish or 'betta', are renowned for their astonishing morphological diversity and extreme aggressive behavior. Despite recent advances in our understanding of the genetics and neurobiology of betta, the lack of tools to manipulate their genome has hindered progress at functional and mechanistic levels. In this study, we outline the use of three genetic manipulation technologies, which we have optimized for use in betta: CRISPR/Cas9-mediated knockout, CRISPR/Cas9-mediated knockin, and Tol2-mediated transgenesis. We knocked out three genes: alkal2l, bco1l , and mitfa , and analyzed their effects on viability and pigmentation. Furthermore, we successfully knocked in a fluorescent protein into the mitfa locus, a proof-of-principle experiment of this powerful technology in betta. Finally, we used Tol2-mediated transgenesis to create fish with ubiquitous expression of GFP, and then developed a bicistronic plasmid with heart-specific expression of a red fluorescent protein to serve as a visible marker of successful transgenesis. Our work highlights the potential for the genetic manipulation of betta, providing valuable resources for the effective use of genetic tools in this animal model.

Redundancy analysis, genome-wide association studies and the pigmentation of brown trout (Salmo trutta L.)

The association of molecular variants with phenotypic variation is a main issue in biology, often tackled with genome-wide association studies (GWAS). GWAS are challenging, with increasing, but still limited, use in evolutionary biology. We used redundancy analysis (RDA) as a complimentary ordination approach to single- and multitrait GWAS to explore the molecular basis of pigmentation variation in brown trout (Salmo trutta) belonging to wild populations impacted by hatchery fish. Based on 75,684 single nucleotide polymorphic (SNP) markers, RDA, single- and multitrait GWAS allowed the extraction of 337 independent colour patterning loci (CPLs) associated with trout pigmentation traits, such as the number of red and black spots on flanks. Collectively, these CPLs (i) mapped onto 35 out of 40 brown trout linkage groups indicating a polygenic genomic architecture of pigmentation, (ii) were found to be associated with 218 candidate genes, including 197 genes formerly mentioned in the literature associated to skin pigmentation, skin patterning, differentiation or structure notably in a close relative, the rainbow trout (Onchorhynchus mykiss), and (iii) related to functions relevant to pigmentation variation (e.g., calcium- and ion-binding, cell adhesion). Annotated CPLs include genes with well-known pigmentation effects (e.g., PMEL, SLC45A2, SOX10), but also markers associated with genes formerly found expressed in rainbow or brown trout skins. RDA was also shown to be useful to investigate management issues, especially the dynamics of trout pigmentation submitted to several generations of hatchery introgression.

Platyfish bypass the constraint of the caudal fin ventral identity in teleosts

BACKGROUND

The caudal fin of teleosts is characterized by dorsoventral symmetry. Despite this external morphology, the principal rays of this appendage connect to bones below the notochord, indicating the ventral (hypochordal) identity of this organ.

RESULTS

Here, we report that this typical architecture of the caudal fin is not fully conserved in the platyfish (Xiphophorus maculatus) and the guppy (Poecilia reticulata), representatives of the Poeciliidae family. We show that in these species, 3-4 principal rays connect to bones above the notochord, suggesting an epichordal contribution. Consistently, as examined in platyfish, dorsal identity genes zic1/4 were highly expressed in these rays, providing molecular evidence of their epichordal origin. Developmental analysis revealed that the earliest rays above the notochord emerge at the 10-ray stage of fin morphogenesis. In contrast to zebrafish and medaka, platyfish and guppies display a mirrored shape of dorsal and ventral processes of the caudal endoskeleton. Our study suggests that an ancestral bauplan expanded in poeciliids by advancing its symmetrical pattern.

CONCLUSION

The platyfish evolved a fin architecture with the epichordal origin of its upper principal rays and a high level of symmetry in the caudal endoskeleton. This innovative architecture highlights the adaptation of the teleost skeleton.

Linked selection, differential introgression and recombination rate variation promote heterogeneous divergence in a pair of yellow croakers

Understanding the mechanisms underlying heterogeneous genomic divergence is of particular interest in evolutionary biology. Highly differentiated genomic regions, known as genomic islands, often evolve between diverging lineages. These genomic islands may be related to selection promoting adaptation or reproductive isolation. Based on whole genome assembly and genome-wide RAD sequencing in a pair of yellow croakers (genus: Larimichthys), we investigated the evolutionary processes shaping genomic landscapes of divergence. Demographic modelling indicated that the two species diverged following a secondary contact scenario, where differential introgression and linked selection were suggested to be involved in heterogeneous genomic divergence. We identified reduced recombination rate in genomic islands and a relatively good conservation of both genetic diversity and recombination landscapes between species, which highlight the roles of linked selection and recombination rate variation in promoting heterogeneous divergence in the common ancestral lineage of the two species. In addition, we found a positive correlation between differentiation (F_ST ) and absolute sequence divergence (D_xy ), and elevated D_xy in genomic islands, indicating that the genomic landscape of divergence was not shaped by linked selection alone. Restricted gene flow in highly differentiated regions has probably remodelled the landscape of heterogeneous genomic divergence. This study highlights that highly differentiated genomic regions can also arise from a combination of linked selection and differential gene flow in interaction with varying recombination rates.

The genetic architecture of phenotypic diversity in the Betta fish (Betta splendens)

The Betta fish displays a remarkable variety of phenotypes selected during domestication. However, the genetic basis underlying these traits remains largely unexplored. Here, we report a high-quality genome assembly and resequencing of 727 individuals representing diverse morphotypes of the Betta fish. We show that current breeds have a complex domestication history with extensive introgression with wild species. Using a genome-wide association study, we identify the genetic basis of multiple traits, including coloration patterns, the "Dumbo" phenotype with pectoral fin outgrowth, extraordinary enlargement of body size that we map to a major locus on chromosome 8, the sex determination locus that we map to dmrt1, and the long-fin phenotype that maps to the locus containing kcnj15. We also identify a polygenic signal related to aggression, involving multiple neural system-related genes such as esyt2, apbb2, and pank2. Our study provides a resource for developing the Betta fish as a genetic model for morphological and behavioral research in vertebrates.

Transcriptome comparison for identification of pigmentation-related genes in different color varieties of Siamese fighting fish Betta splendens

Transcriptome comparison was performed to identify genes expressed in skin, muscle and tails of mono-color (Red, Blue, Black, White and Yellow), bi-color (Cambodian) and multi-color (Marble) varieties of Siamese fighting fish Betta splendens. In total, 163,140 unigenes covering 26.348 Gb were found. Of these, 93,899 (57.55 %) unigenes significantly matched at least one database. In total, 5039 differentially expressed genes (DEGs) were found where 2415 genes (47.93 %) showed higher expression and 2624 genes (52.07 %) showed lower expression for all pairwise comparisons. DEGs between paired color varieties were 133-443. Of these, 38-220 genes were more highly expressed while 37-280 genes were more lowly expressed relative to the compared varieties. A total of 897 sequences (148 genes) significantly matched pigmentation-related genes of Danio rerio (E-value < 1e-06). Of these, 19 DEGs were identified. Examples are tyrosinase-related protein 1a (BsTyrp1a), epidermal growth factor receptor (BsEgfr) and neurofibronin 1a (BsNf1a). Moreover, 711,123 SNPs were identified and 1365 of these were located in pigmentation-related genes. Interestingly, an A > C₄₇₄ SNP in the gene BsTrpm7 and an indel (position 3571) in the BsItgb1a gene were found only in Cambodian. A C > T₂₅₂₀ SNP in BsFzd4 and 10 of 11 SNPs in BsTyrp1a were found only in Black. Different expression levels (P < 0.05) were found for tyrosinase (BsTyr), BsTyrp1a, BsNf1a and BsEgf1 among skin, body muscle and tails of the same variety and among the same tissues of different varieties (Red, Green, Blue, Black, Cambodian and Multi-colors, N = 5 each).

A chromosome-level genome assembly of chia provides insights into high omega-3 content and coat color variation of its seeds

Chia (Salvia hispanica) is a functional food crop for humans. Although its seeds contain high omega-3 fatty acids, the seed yield of chia is still low. Genomic resources available for this plant are limited. We report the first high-quality chromosome-level genome sequence of chia. The assembled genome size was 347.6 Mb and covered 98.1% of the estimated genome size. A total of 31 069 protein-coding genes were predicted. The absence of recent whole-genome duplication and the relatively low intensity of transposable element expansion in chia compared to its sister species contribute to its small genome size. Transcriptome sequencing and gene duplication analysis reveal that the expansion of the fab2 gene family is likely to be related to the high content of omega-3 in seeds. The white seed coat color is determined by a single locus on chromosome 4. This study provides novel insights into the evolution of Salvia species and high omega-3 content, as well as valuable genomic resources for genetic improvement of important commercial traits of chia and its related species.

Care and Use of Siamese Fighting Fish (Betta Splendens) for Research

Betta splendens, also called Siamese fighting fish or 'betta,' are a popular species in the fishkeeping hobby. Native to South- east Asia, betta have been selectively bred for their fighting ability for hundreds of years, which has resulted in the species' characteristic male aggression. More recently, betta have been bred for a number of ornamental traits such as coloration, fin morphology, and body size. Betta have unique characteristics and an evolutionary history that make them a useful model for studies in the fields of behavior, endocrinology, neurobiology, genetics, development, and evolution. However, standard laboratory procedures for raising and keeping these fish are not well established, which has limited their use. Here, we briefly review the past and present use of betta in research, with a focus on their utility in behavioral, neurobiological, and evolutionary studies. We then describe effective husbandry practices for maintaining betta as a research colony.

The Insertion in the 3' UTR of Pmel17 Is the Causal Variant for Golden Skin Color in Tilapia

Understanding of the relationships between genotypes and phenotypes is a central problem in biology. Although teleosts have colorful phenotypes, not much is known about their underlying mechanisms. Our previous study showed that golden skin color in Mozambique tilapia was mapped in the major locus containing the Pmel gene, and an insertion in 3' UTR of Pmel17 was fully correlated with the golden color. However, the molecular mechanism of how Pmel17 determines the golden skin color is unknown. In this study, knockout of Pmel17 with CRISPR/Cas9 in blackish tilapias resulted in golden coloration, and rescue of Pmel17 in golden tilapias recovered the wild-type blackish color, indicating that Pmel17 is the gene determining the golden and blackish color. Functional analysis in vitro showed that the insertion in the 3' UTR of Pmel17 reduced the transcripts of Pmel17. Our data supplies more evidence to support that Pmel17 is the gene for blackish and golden colors, and highlights that the insertion in the 3' UTR of Pmel17 is the causative mutation for the golden coloration.

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.

Transposon-induced epigenetic silencing in the X chromosome as a novel form of dmrt1 expression regulation during sex determination in the fighting fish

Background

Fishes are the one of the most diverse groups of animals with respect to their modes of sex determination, providing unique models for uncovering the evolutionary and molecular mechanisms underlying sex determination and reversal. Here, we have investigated how sex is determined in a species of both commercial and ecological importance, the Siamese fighting fish Betta splendens.

Results

We conducted association mapping on four commercial and two wild populations of B. splendens. In three of the four commercial populations, the master sex determining (MSD) locus was found to be located in a region of ~ 80 kb on LG2 which harbours five protein coding genes, including dmrt1, a gene involved in male sex determination in different animal taxa. In these fish, dmrt1 shows a male-biased gonadal expression from undifferentiated stages to adult organs and the knockout of this gene resulted in ovarian development in XY genotypes. Genome sequencing of XX and YY genotypes identified a transposon, drbx1, inserted into the fourth intron of the X-linked dmrt1 allele. Methylation assays revealed that epigenetic changes induced by drbx1 spread out to the promoter region of dmrt1. In addition, drbx1 being inserted between two closely linked cis-regulatory elements reduced their enhancer activities. Thus, epigenetic changes, induced by drbx1, contribute to the reduced expression of the X-linked dmrt1 allele, leading to female development. This represents a previously undescribed solution in animals relying on dmrt1 function for sex determination. Differentiation between the X and Y chromosomes is limited to a small region of ~ 200 kb surrounding the MSD gene. Recombination suppression spread slightly out of the SD locus. However, this mechanism was not found in the fourth commercial stock we studied, or in the two wild populations analysed, suggesting that it originated recently during domestication.

Conclusions

Taken together, our data provide novel insights into the role of epigenetic regulation of dmrt1 in sex determination and turnover of SD systems and suggest that fighting fish are a suitable model to study the initial stages of sex chromosome evolution.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01205-y.

Masculinizing Effects of Chrysin-Loaded Poloxamer Micelles on Siamese Fighting Fish

Siamese fighting fish (Betta splendens) are freshwater fish that are commonly found in Thailand and other Southeast Asian countries. In the present study, chrysin-loaded polymeric micelles (CPs) were developed and investigated for the masculinizing effects, survival rate, growth indices, and toxicity on Siamese fighting fish. CPs were prepared using a poloxamer. The micelle system of CPs that were formed at a chrysin-to-polymer ratio of 1:2 was found to be the most suitable monodispersed system and exhibited a nanosized diameter (74.2 ± 1.6 nm) with a narrow size distribution (0.288 ± 0.012). In vivo studies were performed using Siamese fighting fish larvae as animal models. In the in vivo toxicity study, the fish larvae were immersed in aqueous systems containing CPs that had five different chrysin concentrations of 1, 10, 100, 1000, and 10,000 ng/mL for 24, 48, and 72 h. Blank polymeric micelles and water were used as controls. The in vivo masculinization effect of CPs with different chrysin concentrations on the fish larvae was evaluated after 5 weeks of exposure. The results demonstrated that CPs with a chrysin concentration of 1000 ng/mL showed a masculinization effect of 94.59 ± 2.76% with a high fish larvae survival rate of 72.45 ± 5.09% and low toxicity. It was concluded that the developed CPs had a significant effect on the sex reversal of Siamese fighting fish larvae with a high survival rate.