Diversity and Convergence of Sex-Determination Mechanisms in Teleost Fish

Sexual reproduction is prevalent across diverse taxa. However, sex-determination mechanisms are so diverse that even closely related species often differ in sex-determination systems. Teleost fish is a taxonomic group with frequent turnovers of sex-determining mechanisms and thus provides us with great opportunities to investigate the molecular and evolutionary mechanisms underlying the turnover of sex-determining systems. Here, we compile recent studies on the diversity of sex-determination mechanisms in fish. We demonstrate that genes in the TGF-β signaling pathway are frequently used for master sex-determining (MSD) genes. MSD genes arise via two main mechanisms, duplication-and-transposition and allelic mutations, with a few exceptions. We also demonstrate that temperature influences sex determination in many fish species, even those with sex chromosomes, with higher temperatures inducing differentiation into males in most cases. Finally, we review theoretical models for the turnover of sex-determining mechanisms and discuss what questions remain elusive.

Genomic Approaches Are Improving Taxonomic Representation in Genetic Studies of Speciation

Until recently, our understanding of the genetics of speciation was limited to a narrow group of model species with a specific set of characteristics that made genetic analysis feasible. Rapidly advancing genomic technologies are eliminating many of the distinctions between laboratory and natural systems. In light of these genomic developments, we review the history of speciation genetics, advances that have been gleaned from model and non-model organisms, the current state of the field, and prospects for broadening the diversity of taxa included in future studies. Responses to a survey of speciation scientists across the world reveal the ongoing division between the types of questions that are addressed in model and non-model organisms. To bridge this gap, we suggest integrating genetic studies from model systems that can be reared in the laboratory or greenhouse with genomic studies in related non-models where extensive ecological knowledge exists.

Repositioning of centromere-associated repeats during karyotype evolution in Oryzias fishes

The karyotype, which is the number and shape of chromosomes, is a fundamental characteristic of all eukaryotes. Karyotypic changes play an important role in many aspects of evolutionary processes, including speciation. In organisms with monocentric chromosomes, it was previously thought that chromosome number changes were mainly caused by centric fusions and fissions, whereas chromosome shape changes, that is, changes in arm numbers, were mainly due to pericentric inversions. However, recent genomic and cytogenetic studies have revealed examples of alternative cases, such as tandem fusions and centromere repositioning, found in the karyotypic changes within and between species. Here, we employed comparative genomic approaches to investigate whether centromere repositioning occurred during karyotype evolution in medaka fishes. In the medaka family (Adrianichthyidae), the three phylogenetic groups differed substantially in their karyotypes. The Oryzias latipes species group has larger numbers of chromosome arms than the other groups, with most chromosomes being metacentric. The O. javanicus species group has similar numbers of chromosomes to the O. latipes species group, but smaller arm numbers, with most chromosomes being acrocentric. The O. celebensis species group has fewer chromosomes than the other two groups and several large metacentric chromosomes that were likely formed by chromosomal fusions. By comparing the genome assemblies of O. latipes, O. javanicus, and O. celebensis, we found that repositioning of centromere-associated repeats might be more common than simple pericentric inversion. Our results demonstrated that centromere repositioning may play a more important role in karyotype evolution than previously appreciated.

Mechanisms of Intrinsic Postzygotic Isolation: From Traditional Genic and Chromosomal Views to Genomic and Epigenetic Perspectives

Intrinsic postzygotic isolation typically appears as reduced viability or fertility of interspecific hybrids caused by genetic incompatibilities between diverged parental genomes. Dobzhansky-Muller interactions among individual genes, and chromosomal rearrangements causing problems with chromosome synapsis and recombination in meiosis, have both long been considered as major mechanisms behind intrinsic postzygotic isolation. Recent research has, however, suggested that the genetic basis of intrinsic postzygotic isolation can be more complex and involves, for example, overall divergence of the DNA sequence or epigenetic changes. Here, we review the mechanisms of intrinsic postzygotic isolation from genic, chromosomal, genomic, and epigenetic perspectives across diverse taxa. We provide empirical evidence for these mechanisms, discuss their importance in the speciation process, and highlight questions that remain unanswered.

Strontium isotope analysis of otoliths reveals differences in the habitat salinity among three sympatric stickleback species of the genus Pungitius

The analysis of otolith Sr isotope ratios (87Sr/86Sr) is a powerful method to study fish migration in freshwater areas. However, few studies have applied this method to study fish movement in brackish-water environments. Furthermore, despite the fact that habitat differentiation has been shown to drive genetic differentiation and reproductive isolation among stickleback fish, no studies have used the otolith 87Sr/86Sr ratios to analyze habitat differentiation between stickleback ecotypes and species. In this study, we analyzed the otolith 87Sr/86Sr ratios of three sympatric stickleback species of the genus Pungitius in the Shiomi River on Hokkaido Island, Japan: P. tymensis, the brackish-water type of the P. pungitius-P. sinensis complex, and the freshwater type of the P. pungitius-P. sinensis complex. First, we created a mixing equation to depict the relationship between habitat salinity and the 87Sr/86Sr ratios of river water. We found that the otolith 87Sr/86Sr ratios differed significantly among the three species, indicating that the three species utilize habitats with different salinities: P. tymensis and the brackish-water type inhabit freshwater and brackish-water environments, respectively, with the freshwater type using intermediate habitats. In addition, we found that some freshwater individuals moved to habitats with higher salinities as they grew. Our study demonstrates that the analysis of otolith 87Sr/86Sr ratios is a useful method for studying the habitat use of fish in brackish-water environments and habitat differentiation among closely related sympatric and parapatric species.

A Cryptic Sex-Linked Locus Revealed by the Elimination of a Master Sex-Determining Locus in Medaka Fish

AbstractSex chromosomes rapidly turn over in several taxonomic groups. Sex chromosome turnover is generally thought to start with the appearance of a new sex-determining gene on an autosome while an old sex-determining gene still exists, followed by the fixation of the new one. However, we do not know how prevalent the transient state is, where multiple sex-determining loci coexist within natural populations. Here, we removed a Y chromosome with a master male-determining gene DMY from medaka fish using high temperature-induced sex-reversed males. After four generations, the genomic characteristics of a sex chromosome were found on one chromosome, which was an autosome in the original population. Thus, the elimination of a master sex-determining locus can reveal a cryptic locus with a possible sex-determining effect, which can be the seed for sex chromosome turnover. Our results suggest that populations that seem to have a single-locus XY system may have other chromosomal regions with sex-determining effects. In conclusion, the coexistence of multiple sex-determining genes in a natural population may be more prevalent than previously thought. Experimental elimination of a master sex-determining locus may serve as a promising method for finding a locus that can be a protosex chromosome.

Do sex-linked male meiotic drivers contribute to intrinsic hybrid incompatibilities? Recent empirical studies from flies and rodents

Intrinsic hybrid incompatibility is one of the important isolating barriers between species. In organisms with sex chromosomes, intrinsic hybrid incompatibility often follows two rules: Haldane's rule and large-X effects. One explanation for these two rules is that sex chromosomes are hotspots for meiotic drivers that can cause intrinsic hybrid incompatibility between geographically isolated populations. Although this hypothesis seems plausible and several empirical data are consistent with it, we are still unsure whether such mechanisms occur in nature, particularly with respect to speciation with gene flow. Here, we review empirical studies that have investigated the roles of meiotic drive in sex-chromosome evolution and speciation and propose future studies necessary for testing this hypothesis.

Plate reduction in southern Japanese freshwater populations of threespine stickleback (Gasterosteus aculeatus)

Adaptation to similar environments can lead to the evolution of similar phenotypes in phylogenetically independent lineages. However, the extent of parallel evolution often varies. Because such variations can be due to environmental heterogeneity among seemingly similar habitats, identification of the environmental factors that cause non-parallel patterns can provide valuable insight into the ecological factors associated with phenotypic diversification. Armor plate reduction in replicate freshwater populations of the threespine stickleback (Gasterosteus aculeatus) represents a well-known example of parallel evolution. Many freshwater populations in multiple regions of the Northern Hemisphere have reduced plate numbers, but not all freshwater populations exhibit plate reduction. In this study, we characterized plate number variation in Japanese freshwater populations and investigated the association between plate number and several abiotic environmental factors. We found that most freshwater populations have not reduced plate numbers in Japan. Plate reduction tends to occur in habitats with warmer winter temperatures at lower latitudes in Japan. In contrast, low dissolved calcium levels or water turbidity had no significant effects on plate reduction, although these were reported to be associated with plate reduction in Europe. Although our data are consistent with the hypothesis that winter temperatures are associated with plate reduction, further studies on the relationship between temperatures and fitness using sticklebacks with varying plate numbers are necessary to confirm this hypothesis and understand the factors causing variations in the extent of parallel evolution.

Multiple colonizations and hybridization of a freshwater fish group on a satellite island of Sulawesi

Repeated colonizations and resultant hybridization may increase lineage diversity on an island if introgression occurs only in a portion of the indigenous island lineage. Therefore, to precisely understand how island biodiversity was shaped, it is essential to reconstruct the history of secondary colonization and resultant hybridization both in time and space. In this study, we reconstructed the history of multiple colonizations of the Oryzias woworae species group, a freshwater fish group of the family Adrianichthyidae, from Sulawesi Island to its southeast satellite island, Muna Island. Phylogenetic and species tree analyses using genome-wide single-nucleotide polymorphisms revealed that all local populations on Muna Island were monophyletic, but that there were several genetically distinct lineages within the island. Population structure and phylogenetic network analyses demonstrated that colonization of this island occurred more than once, and that secondary colonization and resultant introgressive hybridization occurred only in one local population on the island. The spatially heterogeneous introgression induced by the multiple colonizations were also supported by differential admixture analyses. In addition, the differential admixture analyses detected reverse colonization from Muna Island to the Sulawesi mainland. Coalescence-based demographic inference estimated that these mutual colonizations occurred during the middle to late Quaternary period, during which sea level repeatedly declined; this indicates that the colonizations occurred via land bridges. We conclude that these mutual colonizations between Muna Island and the Sulawesi mainland, and the resultant spatially heterogeneous introgression shaped the current biodiversity of this species group in this area.

Convergent copy number increase of genes associated with freshwater colonization in fishes

Copy number variation (CNV) can cause phenotypic changes. However, in contrast to amino acid substitutions and cis-regulatory changes, little is known about the functional categories of genes in which CNV is important for adaptation to novel environments. It is also unclear whether the same genes repeatedly change the copy numbers for adapting to similar environments. Here, we investigate CNV associated with freshwater colonization in fishes, which was observed multiple times across different lineages. Using 48 ray-finned fishes across diverse orders, we identified 23 genes whose copy number increases were associated with freshwater colonization. These genes showed enrichment for peptide receptor activity, hexosyltransferase activity and unsaturated fatty acid metabolism. We further revealed that three of the genes showed copy number increases in freshwater populations compared to marine ancestral populations of the stickleback genus Gasterosteus. These results indicate that copy number increases of genes involved in fatty acid metabolism (FADS2), immune function (PSMB8a) and thyroid hormone metabolism (UGT2) may be important for freshwater colonization at both the inter-order macroevolutionary scale and at the intra-genus microevolutionary scale. Further analysis across diverse taxa will help to understand the role of CNV in the adaptation to novel environments. This article is part of the theme issue 'Genetic basis of adaptation and speciation: from loci to causative mutations'.

Speciation and adaptation research meets genome editing

Understanding the genetic basis of reproductive isolation and adaptive traits in natural populations is one of the fundamental goals in evolutionary biology. Genome editing technologies based on CRISPR-Cas systems and site-specific recombinases have enabled us to modify a targeted genomic region as desired and thus to conduct functional analyses of target loci, genes and mutations even in non-conventional model organisms. Here, we review the technical properties of genome editing techniques by classifying them into the following applications: targeted gene knock-out for investigating causative gene functions, targeted gene knock-in of marker genes for visualizing expression patterns and protein functions, precise gene replacement for identifying causative alleles and mutations, and targeted chromosomal rearrangement for investigating the functional roles of chromosomal structural variations. We describe examples of their application to demonstrate functional analysis of naturally occurring genetic variations and discuss how these technologies can be applied to speciation and adaptation research. This article is part of the theme issue 'Genetic basis of adaptation and speciation: from loci to causative mutations'.

Genetic basis for the evolution of pelvic-fin brooding, a new mode of reproduction, in a Sulawesian fish

Modes of reproduction in animals are diverse, with different modes having evolved independently in multiple lineages across a variety of taxa. However, an understanding of the genomic change driving the transition between different modes of reproduction is limited. Several ricefishes (Adrianichthyidae) on the island of Sulawesi have a unique mode of reproduction called "pelvic-fin brooding," wherein females carry externally fertilized eggs until hatching using their pelvic fins. Phylogenomic analysis demonstrated pelvic-fin brooders to have evolved at least twice in two distant clades of the Adrianichthyidae. We investigated the genetic architecture of the evolution of this unique mode of reproduction. Morphological analyses and laboratory observations revealed that females of pelvic-fin brooders have longer pelvic fins and a deeper abdominal concavity, and that they can carry an egg clutch for longer than non-brooding adrianichthyids, suggesting that these traits play important roles in this reproductive mode. Quantitative trait locus mapping using a cross between a pelvic-fin brooder Oryzias eversi and a non-brooding O. dopingdopingensis reveals different traits involved in pelvic-fin brooding to be controlled by different loci on different chromosomes. Genomic analyses of admixture detected no signatures of introgression between two lineages with pelvic-fin brooders, indicating that introgression is unlikely to be responsible for repeated evolution of pelvic-fin brooding. These findings suggest that multiple independent mutations may have contributed to the convergent evolution of this novel mode of reproduction.

Local Adaptation and the Evolution of Genome Architecture in Threespine Stickleback

Theory predicts that local adaptation should favor the evolution of a concentrated genetic architecture, where the alleles driving adaptive divergence are tightly clustered on chromosomes. Adaptation to marine versus freshwater environments in threespine stickleback has resulted in an architecture that seems consistent with this prediction: divergence among populations is mainly driven by a few genomic regions harboring multiple quantitative trait loci for environmentally adapted traits, as well as candidate genes with well-established phenotypic effects. One theory for the evolution of these "genomic islands" is that rearrangements remodel the genome to bring causal loci into tight proximity, but this has not been studied explicitly. We tested this theory using synteny analysis to identify micro- and macro-rearrangements in the stickleback genome and assess their potential involvement in the evolution of genomic islands. To identify rearrangements, we conducted a de novo assembly of the closely related tubesnout (Aulorhyncus flavidus) genome and compared this to the genomes of threespine stickleback and two other closely related species. We found that small rearrangements, within-chromosome duplications, and lineage-specific genes (LSGs) were enriched around genomic islands, and that all three chromosomes harboring large genomic islands have experienced macro-rearrangements. We also found that duplicates and micro-rearrangements are 9.9× and 2.9× more likely to involve genes differentially expressed between marine and freshwater genotypes. While not conclusive, these results are consistent with the explanation that strong divergent selection on candidate genes drove the recruitment of rearrangements to yield clusters of locally adaptive loci.

Genetic basis of speciation and adaptation: from loci to causative mutations

Does evolution proceed in small steps or large leaps? How repeatable is evolution? How constrained is the evolutionary process? Answering these long-standing questions in evolutionary biology is indispensable for both understanding how extant biodiversity has evolved and predicting how organisms and ecosystems will respond to changing environments in the future. Understanding the genetic basis of phenotypic diversification and speciation in natural populations is key to properly answering these questions. The leap forward in genome sequencing technologies has made it increasingly easier to not only investigate the genetic architecture but also identify the variant sites underlying adaptation and speciation in natural populations. Furthermore, recent advances in genome editing technologies are making it possible to investigate the functions of each candidate gene in organisms from natural populations. In this article, we discuss how these recent technological advances enable the analysis of causative genes and mutations and how such analysis can help answer long-standing evolutionary biology questions.

This article is part of the theme issue ‘Genetic basis of adaptation and speciation: from loci to causative mutations’.

Construction of a chromosome-level Japanese stickleback species genome using ultra-dense linkage analysis with single-cell sperm sequencing

It is still difficult to construct the genomes of higher organisms as their genome sequences must be extended to the length of the chromosome by linkage analysis. In this study, we attempted to provide an innovative alternative to conventional linkage analysis by devising a method to genotype sperm using 10× Genomics single-cell genome sequencing libraries to generate a linkage map without interbreeding individuals. A genome was assembled using sperm from the Japanese stickleback Gasterosteus nipponicus, with single-cell genotyping yielding 1 864 430 very dense hetero-SNPs and an average coverage per sperm cell of 0.13×. In total, 1665 sperm were used, which is an order of magnitude higher than the number of recombinations used for conventional linkage analysis. We then improved the linkage analysis tool scaffold extender with low depth linkage analysis (SELDLA) to analyze the data according to the characteristics of the single-cell genotyping data. Finally, we were able to determine the chromosomal location (97.1%) and orientation (64.4%) of the contigs in the 456 Mb genome of G. nipponicus, sequenced using nanopores. This method promises to be a useful tool for determining the genomes of non-model organisms for which breeding systems have not yet been established by linkage analysis.

Species divergence and repeated ancient hybridization in a Sulawesian lake system

An increasing volume of empirical studies demonstrated that hybridization between distant lineages may have promoted speciation in various taxa. However, the timing, extent and direction of introgressive hybridization remain unknown in many cases. Here, we report a possible case in which repeated hybridization promoted divergence of Oryzias ricefishes (Adrianichthyidae) on Sulawesi, an island of Wallacea. Four Oryzias species are endemic to the Malili Lake system in central Sulawesi, which is composed of five tectonic lakes; of these, one lake is inhabited by two species. Morphological and population genomic analyses of genome-wide single-nucleotide polymorphisms revealed that these two sympatric species are phylogenetically sister to but substantially reproductively isolated from each other. Analyses of admixture and comparison of demographic models revealed that the two sympatric species experienced several substantial introgressions from outgroup populations that probably occurred soon after they had secondary contact with each other in the lake. However, the ratio of migrants from the outgroups was estimated to be different between the two species, which is consistent with the hypothesis that these introgressions aided their divergence or prevented them from forming a hybrid swarm. Repeated lake fragmentations and fusions may have promoted diversification of this freshwater fish species complex that is endemic to this ancient lake system.

Copy number variation of a fatty acid desaturase gene Fads2 associated with ecological divergence in freshwater stickleback populations

Fitness of aquatic animals can be limited by the scarcity of nutrients such as long-chain polyunsaturated fatty acids, especially docosahexaenoic acid (DHA). DHA availability from diet varies among aquatic habitats, imposing different selective pressures on resident animals to optimize DHA acquisition and synthesis. For example, DHA is generally poor in freshwater ecosystems compared to marine ecosystems. Our previous work revealed that, relative to marine fishes, several freshwater fishes evolved higher copy numbers of the fatty acid desaturase2 (Fads2) gene, which encodes essential enzymes for DHA biosynthesis, likely compensating for the limited availability of DHA in freshwater. Here, we demonstrate that Fads2 copy number also varies between freshwater sticklebacks inhabiting lakes and streams with stream fish having higher Fads2 copy number. Additionally, populations with benthic-like morphology possessed higher Fads2 copy number than those with planktivore-like morphology. This may be because benthic-like fish mainly feed on DHA-deficient prey such as macroinvertebrates whereas planktivore-like fish forage more regularly on DHA-rich prey, like copepods. Our results suggest that Fads2 copy number variation arises from ecological divergence not only between organisms exploiting marine and freshwater habitats but also between freshwater organisms exploiting divergent resources.

Mesozoic origin and 'out-of-India' radiation of ricefishes (Adrianichthyidae)

The Indian subcontinent has an origin geologically different from Eurasia, but many terrestrial animal and plant species on it have congeneric or sister species in other parts of Asia, especially in the Southeast. This faunal and floral similarity between India and Southeast Asia is explained by either of the two biogeographic scenarios, 'into-India' or 'out-of-India'. Phylogenies based on complete mitochondrial genomes and five nuclear genes were undertaken for ricefishes (Adrianichthyidae) to examine which of these two biogeographic scenarios fits better. We found that Oryzias setnai, the only adrianichthyid distributed in and endemic to the Western Ghats, a mountain range running parallel to the western coast of the Indian subcontinent, is sister to all other adrianichthyids from eastern India and Southeast-East Asia. Divergence time estimates and ancestral area reconstructions reveal that this western Indian species diverged in the late Mesozoic during the northward drift of the Indian subcontinent. These findings indicate that adrianichthyids dispersed eastward 'out-of-India' after the collision of the Indian subcontinent with Eurasia, and subsequently diversified in Southeast-East Asia. A review of geographic distributions of 'out-of-India' taxa reveals that they may have largely fuelled or modified the biodiversity of Eurasia.

The evolutionary ecology of fatty-acid variation: Implications for consumer adaptation and diversification

The nutritional diversity of resources can affect the adaptive evolution of consumer metabolism and consumer diversification. The omega-3 long-chain polyunsaturated fatty acids eicosapentaenoic acid (EPA; 20:5n-3) and docosahexaenoic acid (DHA; 22:6n-3) have a high potential to affect consumer fitness, through their widespread effects on reproduction, growth and survival. However, few studies consider the evolution of fatty acid metabolism within an ecological context. In this review, we first document the extensive diversity in both primary producer and consumer fatty acid distributions amongst major ecosystems, between habitats and amongst species within habitats. We highlight some of the key nutritional contrasts that can shape behavioural and/or metabolic adaptation in consumers, discussing how consumers can evolve in response to the spatial, seasonal and community-level variation of resource quality. We propose a hierarchical trait-based approach for studying the evolution of consumers' metabolic networks and review the evolutionary genetic mechanisms underpinning consumer adaptation to EPA and DHA distributions. In doing so, we consider how the metabolic traits of consumers are hierarchically structured, from cell membrane function to maternal investment, and have strongly environment-dependent expression. Finally, we conclude with an outlook on how studying the metabolic adaptation of consumers within the context of nutritional landscapes can open up new opportunities for understanding evolutionary diversification.

Multiple paths to the same destination: Influence of gene flow on convergent evolution

Elucidation of the genetic mechanisms of convergent evolution, the evolution of similar or the same phenotypes in phylogenetically independent lineages, helps predict how populations will respond to the same selective pressures. Convergent evolution can be caused by either the fixation of identical-by-descent alleles, independent mutations at the same gene, or mutations in different genes controlling the same trait. To what extent does the fixation of identical-by-descent alleles lead to convergent evolution in isolated populations where inflow of adaptive alleles from other populations is limited? In a From the Cover article in this issue of Molecular Ecology, Kemppainen et al. (2021) compared the genetic basis for the reduction of pelvic structures in three isolated freshwater populations of nine-spined stickleback (Pungitius pungitius) from Northern Europe. The authors used quantitative trait loci (QTL) mapping to reveal that the pelvic reduction in these three populations was caused by mutations at different genetic loci. In contrast to studies in three-spined stickleback (Gasterosteus aculeatus), where independently derived Pitx1 mutations were shown to be responsible for plate reduction across multiple freshwater populations, Kemppainen et al. (2021) found Pitx1 to be the candidate causative gene for only one population of P. pungitius. This study highlights the importance of genetic studies of convergent evolution, not only in the presence of gene flow but also in its absence for a better understanding of the genetic architecture of convergent evolution.

Genetic basis for variation in the number of cephalic pores in a hybrid zone between closely related species of goby, Gymnogobius breunigii and Gymnogobius castaneus

Populations or species exploiting different habitats can differ in sensory perception as a result of divergent adaptation. In bony fish, the water current is perceived via neuromasts, the end organ of the lateral line system. Although fish in different habitats are known to vary in neuromasts, we know little about the genetic basis for such variation. Here, we investigate the genetic basis for variation in supraorbital neuromasts in a hybrid zone between the Japanese gobies Gymnogobius breunigii and Gymnogobius castaneus. The former has supraorbital canal neuromasts with six cephalic pores, whereas the latter has only superficial neuromasts with no canals or pores in the supraorbital region. Our genomic analysis showed that G. breunigii and G. castaneus occur mainly in the lower and mid/upper reaches, respectively. In a river in northern Japan, hybrids were found at the sites between the habitats of the two species. These hybrids exhibited anomalies of cephalic pores. Using this hybrid zone, we conducted genome-wide association studies and identified one locus significantly associated with the number of pores. Genomic cline analysis in the hybrid zone demonstrated that this locus exhibited a higher introgression rate compared with the genomic background, indicating the possibility of adaptive introgression.

Tempo and mode in karyotype evolution revealed by a probabilistic model incorporating both chromosome number and morphology

Karyotype, including the chromosome and arm numbers, is a fundamental genetic characteristic of all organisms and has long been used as a species-diagnostic character. Additionally, karyotype evolution plays an important role in divergent adaptation and speciation. Centric fusion and fission change chromosome numbers, whereas the intra-chromosomal movement of the centromere, such as pericentric inversion, changes arm numbers. A probabilistic model simultaneously incorporating both chromosome and arm numbers has not been established. Here, we built a probabilistic model of karyotype evolution based on the "karyograph", which treats karyotype evolution as a walk on the two-dimensional space representing the chromosome and arm numbers. This model enables analysis of the stationary distribution with a stable karyotype for any given parameter. After evaluating their performance using simulated data, we applied our model to two large taxonomic groups of fish, Eurypterygii and series Otophysi, to perform maximum likelihood estimation of the transition rates and reconstruct the evolutionary history of karyotypes. The two taxa significantly differed in the evolution of arm number. The inclusion of speciation and extinction rates demonstrated possibly high extinction rates in species with karyotypes other than the most typical karyotype in both groups. Finally, we made a model including polyploidization rates and applied it to a small plant group. Thus, the use of this probabilistic model can contribute to a better understanding of tempo and mode in karyotype evolution and its possible role in speciation and extinction.

Genome editing reveals fitness effects of a gene for sexual dichromatism in Sulawesian fishes

Sexual selection drives rapid phenotypic diversification of mating traits. However, we know little about the causative genes underlying divergence in sexually selected traits. Here, we investigate the genetic basis of male mating trait diversification in the medaka fishes (genus Oryzias) from Sulawesi, Indonesia. Using linkage mapping, transcriptome analysis, and genome editing, we identify csf1 as a causative gene for red pectoral fins that are unique to male Oryzias woworae. A cis-regulatory mutation enables androgen-induced expression of csf1 in male fins. csf1-knockout males have reduced red coloration and require longer for mating, suggesting that coloration can contribute to male reproductive success. Contrary to expectations, non-red males are more attractive to a predatory fish than are red males. Our results demonstrate that integrating genomics with genome editing enables us to identify causative genes underlying sexually selected traits and provides a new avenue for testing theories of sexual selection.

Multiple waves of freshwater colonization of the three-spined stickleback in the Japanese Archipelago

Background

The three-spined stickleback (Gasterosteus aculeatus) is a remarkable system to study the genetic mechanisms underlying parallel evolution during the transition from marine to freshwater habitats. Although the majority of previous studies on the parallel evolution of sticklebacks have mainly focused on postglacial freshwater populations in the Pacific Northwest of North America and northern Europe, we recently use Japanese stickleback populations for investigating shared and unique features of adaptation and speciation between geographically distant populations. However, we currently lack a comprehensive phylogeny of the Japanese three-spined sticklebacks, despite the fact that a good phylogeny is essential for any evolutionary and ecological studies. Here, we conducted a phylogenomic analysis of the three-spined stickleback in the Japanese Archipelago.

Results

We found that freshwater colonization occurred in multiple waves, each of which may reflect different interglacial isolations. Some of the oldest freshwater populations from the central regions of the mainland of Japan (hariyo populations) were estimated to colonize freshwater approximately 170,000 years ago. The next wave of colonization likely occurred approximately 100,000 years ago. The inferred origins of several human-introduced populations showed that introduction occurred mainly from nearby habitats. We also found a new habitat of the three-spined stickleback sympatric with the Japan Sea stickleback (Gasterosteus nipponicus).

Conclusions

These Japanese stickleback systems differ from those in the Pacific Northwest of North America and northern Europe in terms of divergence time and history. Stickleback populations in the Japanese Archipelago offer valuable opportunities to study diverse evolutionary processes in historical and contemporary timescales.

Differences in the contributions of sex linkage and androgen regulation to sex-biased gene expression in juvenile and adult sticklebacks

Different evolutionary interests between males and females can lead to the evolution of sexual dimorphism. However, intersex genetic correlations due to the shared genome can constrain the evolution of sexual dimorphism, resulting in intra-locus sexual conflict. One of the mechanisms resolving this conflict is sex linkage, which allows males and females to carry different alleles on sex chromosomes. Another is a regulatory mutation causing sex-biased gene expression, which is often mediated by gonadal steroids in vertebrates. How do these two mechanisms differ in the contributions to the resolution of intra-locus sexual conflict? The magnitude of sexual conflict often varies between the juvenile and adult stages. Because gonadal steroids change in titre during development, we hypothesized that gonadal steroids play a role in sexual dimorphism expression only at certain developmental stages, whereas sex linkage is more important for sexual dimorphism expressed throughout life. Our brain transcriptome analysis of juvenile and adult threespine sticklebacks showed that the majority of genes that were sex-biased in both stages were sex-linked. The relative contribution of androgen-dependent regulation to the sex-biased transcriptome increased and that of sex linkage declined in adults compared to juveniles. The magnitude of the sex differences was greater in sex-linked genes than androgen-responsive genes, suggesting that sex linkage is more effective than androgen regulation in the production of large sex differences in gene expression. Overall, our data are consistent with the hypothesis that sex linkage is effective in resolving sexual conflict throughout life, whereas androgen-dependent regulation can contribute to temporary resolution of sexual conflict.

Genome-wide patterns of divergence and introgression after secondary contact between Pungitius sticklebacks

Speciation is a continuous process. Although it is known that differential adaptation can initiate divergence even in the face of gene flow, we know relatively little about the mechanisms driving complete reproductive isolation and the genomic patterns of divergence and introgression at the later stages of speciation. Sticklebacks contain many pairs of sympatric species differing in levels of reproductive isolation and divergence history. Nevertheless, most previous studies have focused on young species pairs. Here, we investigated two sympatric stickleback species, Pungitius pungitius and P. sinensis, whose habitats overlap in eastern Hokkaido; these species show hybrid male sterility, suggesting that they may be at a late stage of speciation. Our demographic analysis using whole-genome sequence data showed that these species split 1.73 Ma and came into secondary contact 37 200 years ago after a period of allopatry. This long period of allopatry might have promoted the evolution of intrinsic incompatibility. Although we detected on-going gene flow and signatures of introgression, overall genomic divergence was high, with considerable heterogeneity across the genome. The heterogeneity was significantly associated with variation in recombination rate. This sympatric pair provides new avenues to investigate the late stages of the stickleback speciation continuum. This article is part of the theme issue 'Towards the completion of speciation: the evolution of reproductive isolation beyond the first barriers'.

Patterns of genomic divergence and introgression between Japanese stickleback species with overlapping breeding habitats

With only a few absolute geographic barriers in marine environments, the factors maintaining reproductive isolation among marine organisms remain elusive. However, spatial structuring in breeding habitat can contribute to reproductive isolation. This is particularly important for marine organisms that migrate to use fresh- or brackish water environments to breed. The Japanese Gasterosteus stickleback species, the Pacific Ocean three-spined stickleback (G. aculeatus) and the Japan Sea stickleback (G. nipponicus) overwinter in the sea, but migrate to rivers for spawning. Although they co-occur at several locations across the Japanese islands, they are reproductively isolated. Our previous studies in Bekanbeushi River showed that the Japan Sea stickleback spawns in the estuary, while the Pacific Ocean stickleback mainly spawns further upstream in freshwater. Overall genomic divergence was very high with many interspersed regions of introgression. Here, we investigated genomic divergence and introgression between the sympatric species in the much shorter Tokotan River, where they share spawning sites. The levels of genome-wide divergence were reduced and introgression was increased, suggesting that habitat isolation substantially contributes to a reduction in gene flow. We also found that genomic regions of introgression were largely shared between the two systems. Furthermore, some regions of introgression were located near loci with a heterozygote advantage for juvenile survival. Taken together, introgression may be partially driven by adaptation in this system. Although, the two species remain clearly genetically differentiated. Regions with low recombination rates showed especially low introgression. Speciation reversal is therefore likely prevented by barriers other than habitat isolation.

Accumulation of Deleterious Mutations in Landlocked Threespine Stickleback Populations

Colonization of new habitats often reduces population sizes and may result in the accumulation of deleterious mutations by genetic drift. Compared with the genomic basis for adaptation to new environments, genome-wide analysis of deleterious mutations in isolated populations remains limited. In the present study, we investigated the accumulation of deleterious mutations in five endangered freshwater populations of threespine stickleback (Gasterosteus aculeatus) in the central part of the mainland of Japan. Using whole-genome resequencing data, we first conducted phylogenomic analysis and confirmed at least two independent freshwater colonization events in the central mainland from ancestral marine ecotypes. Next, analyses of single nucleotide polymorphisms showed a substantial reduction of heterozygosity in freshwater populations compared with marine populations. Reduction in heterozygosity was more apparent at the center of each chromosome than the peripheries and on X chromosomes compared with autosomes. Third, bioinformatic analysis of deleterious mutations showed increased accumulation of putatively deleterious mutations in the landlocked freshwater populations compared with marine populations. For the majority of populations examined, the frequencies of putatively deleterious mutations were higher on X chromosomes than on autosomes. The interpopulation comparison indicated that the majority of putatively deleterious mutations may have accumulated independently. Thus, whole-genome resequencing of endangered populations can help to estimate the accumulation of deleterious mutations and inform us of which populations are the most severely endangered. Furthermore, analysis of variation among chromosomes can give insights into whether any particular chromosomes are likely to accumulate deleterious mutations.

Diversity in reproductive seasonality in the three-spined stickleback, Gasterosteus aculeatus

The annual timing of reproduction is a key life history trait with a large effect on fitness. Populations often vary in the timing and duration of reproduction to adapt to different seasonality of ecological and environmental variables between habitats. However, little is known about the molecular genetic mechanisms underlying interpopulation variation in reproductive seasonality. Here, we demonstrate that the three-spined stickleback (Gasterosteus aculeatus) is a good model for molecular genetic analysis of variations in reproductive seasonality. We first compiled data on reproductive seasons of diverse ecotypes, covering marine-anadromous, lake and stream ecotypes, of three-spined stickleback inhabiting a wide range of latitudes. Our analysis showed that both ecotype and latitude significantly contribute to variation in reproductive seasons. Stream ecotypes tend to start breeding earlier and end later than other ecotypes. Populations from lower latitudes tend to start breeding earlier than those from higher latitudes in all three ecotypes. Additionally, stream ecotypes tend to have extended breeding seasons at lower latitudes than at higher latitudes, leading to nearly year-round reproduction in the most southern stream populations. A review of recent progress in our understanding of the physiological mechanisms underlying seasonal reproduction in the three-spined stickleback indicates that photoperiod is an important external cue that stimulates and/or suppresses reproduction in this species. Taking advantage of genomic tools available for this species, the three-spined stickleback will be a good model to investigate what kinds of genes and mutations underlie variations in the physiological signalling pathways that regulate reproduction in response to photoperiod.

A key metabolic gene for recurrent freshwater colonization and radiation in fishes

Colonization of new ecological niches has triggered large adaptive radiations. Although some lineages have made use of such opportunities, not all do so. The factors causing this variation among lineages are largely unknown. Here, we show that deficiency in docosahexaenoic acid (DHA), an essential ω-3 fatty acid, can constrain freshwater colonization by marine fishes. Our genomic analyses revealed multiple independent duplications of the fatty acid desaturase gene Fads2 in stickleback lineages that subsequently colonized and radiated in freshwater habitats, but not in close relatives that failed to colonize. Transgenic manipulation of Fads2 in marine stickleback increased their ability to synthesize DHA and survive on DHA-deficient diets. Multiple freshwater ray-finned fishes also show a convergent increase in Fads2 copies, indicating its key role in freshwater colonization.

The Origin of a New Sex Chromosome by Introgression between Two Stickleback Fishes

Introgression is increasingly recognized as a source of genetic diversity that fuels adaptation. Its role in the evolution of sex chromosomes, however, is not well known. Here, we confirm the hypothesis that the Y chromosome in the ninespine stickleback, Pungitius pungitius, was established by introgression from the Amur stickleback, P. sinensis. Using whole genome resequencing, we identified a large region of Chr 12 in P. pungitius that is diverged between males and females. Within but not outside of this region, several lines of evidence show that the Y chromosome of P. pungitius shares a most recent common ancestor not with the X chromosome, but with the homologous chromosome in P. sinensis. Accumulation of repetitive elements and gene expression changes on the new Y are consistent with a young sex chromosome in early stages of degeneration, but other hallmarks of Y chromosomes have not yet appeared. Our findings indicate that porous species boundaries can trigger rapid sex chromosome evolution.

Convergent evolution of body color between sympatric freshwater fishes via different visual sensory evolution

Although there are many examples of color evolution potentially driven by sensory drive, only few studies have examined whether distinct species inhabiting the same environments evolve similar body colors via shared sensory mechanisms. In this study, we tested whether two sympatric freshwater fish taxa, halfbeaks of the genus Nomorhamphus and ricefishes of the genus Oryzias in Sulawesi Island, converge in both body color and visual sensitivity. After reconstructing the phylogeny separately for Nomorhamphus and Oryzias using transcriptome-wide sequences, we demonstrated positive correlations of body redness between these two taxa across environments, even after phylogenetic corrections, which support convergent evolution. However, substantial differences were observed in the expression profiles of opsin genes in the eyes between Nomorhamphus and Oryzias. Particularly, the expression levels of the long wavelength-sensitive genes were negatively correlated between the taxa, indicating that they have different visual sensitivities despite living in similar light environments. Thus, the convergence of body colorations between these two freshwater fish taxa was not accompanied by convergence in opsin sensitivities. This system presents a case in which body color convergence can occur between sympatric species via different mechanisms.

A four-questions perspective on public information use in sticklebacks (Gasterosteidae)

Whether learning primarily reflects general processes or species-specific challenges is a long-standing matter of dispute. Here, we present a comprehensive analysis of public information use (PI-use) in sticklebacks (Gasterosteidae). PI-use is a form of social learning by which animals are able to assess the relative quality of resources, here prey patches, by observing the behaviour of others. PI-use was highly specific with only Pungitius and their closest relative Culaea inconstans showing evidence of PI-use. We saw no effects of ontogenetic experience upon PI-use in Pungitius pungitius. Experiments with live demonstrators and animated fish revealed that heightened activity and feeding strikes by foraging conspecifics are important cues in the transmission of PI. Finally, PI-use was the only form of learning in which P. pungitius and another stickleback, Gasterosteus aculeatus differed. PI-use in sticklebacks is species-specific and may represent an 'ecological specialization' for social foraging. Whether this reflects selection on perception, attentional or cognitive processes remains to be determined.

The genomic landscape at a late stage of stickleback speciation: High genomic divergence interspersed by small localized regions of introgression

Speciation is a continuous process and analysis of species pairs at different stages of divergence provides insight into how it unfolds. Previous genomic studies on young species pairs have revealed peaks of divergence and heterogeneous genomic differentiation. Yet less known is how localised peaks of differentiation progress to genome-wide divergence during the later stages of speciation in the presence of persistent gene flow. Spanning the speciation continuum, stickleback species pairs are ideal for investigating how genomic divergence builds up during speciation. However, attention has largely focused on young postglacial species pairs, with little knowledge of the genomic signatures of divergence and introgression in older stickleback systems. The Japanese stickleback species pair, composed of the Pacific Ocean three-spined stickleback (Gasterosteus aculeatus) and the Japan Sea stickleback (G. nipponicus), which co-occur in the Japanese islands, is at a late stage of speciation. Divergence likely started well before the end of the last glacial period and crosses between Japan Sea females and Pacific Ocean males result in hybrid male sterility. Here we use coalescent analyses and Approximate Bayesian Computation to show that the two species split approximately 0.68-1 million years ago but that they have continued to exchange genes at a low rate throughout divergence. Population genomic data revealed that, despite gene flow, a high level of genomic differentiation is maintained across the majority of the genome. However, we identified multiple, small regions of introgression, occurring mainly in areas of low recombination rate. Our results demonstrate that a high level of genome-wide divergence can establish in the face of persistent introgression and that gene flow can be localized to small genomic regions at the later stages of speciation with gene flow.

Impact of the huge 2011 Tohoku-oki tsunami on the phenotypes and genotypes of Japanese coastal threespine stickleback populations

On March 11, 2011, a large earthquake occurred, causing a tsunami which struck the Pacific coast of northeast Japan. We investigated the ecological and genetic effects of the large tsunami on the threespine stickleback (genus Gasterosteus) populations in Otsuchi Town, which was one of the most severely damaged areas after the tsunami. Our environmental surveys showed that spring water may have contributed to the habitat recovery. Morphological analysis of the stickleback before and after the tsunami showed morphological shifts in the gill raker number, which is a foraging trait. Genetic analyses revealed that the allelic richness of one population was maintained after the tsunami, whereas that of another decreased in 2012 and then started to recover in 2013. Additionally, we found that the large tsunami and ground subsidence created new spring water-fed pools with sticklebacks, suggesting that the tsunami brought sticklebacks into these pools. Genetic analysis of this population showed that this population might be derived from hybridization between freshwater Gasterosteus aculeatus and anadromous G. nipponicus. Overall, our data indicate that tsunamis can influence morphologies and genetic structures of freshwater fishes. Furthermore, spring water may play important roles in the maintenance and creation of fish habitats, faced with environmental disturbance.

Contribution of gene flow to the evolution of recombination suppression in sex chromosomes

Polymorphism of alleles that benefit one sex but harm the other (sexually antagonistic alleles) generates selective pressures for reduced recombination between themselves and sex-determination loci. Such polymorphism can be maintained within a population when selection coefficients are sufficiently balanced between males and females. However, if regulatory mutations restrict gene expression only to one sex, these alleles become neutral in the other sex and easily fixed within a population, removing the selective pressures for recombination suppression in sex chromosomes. When there is spatial variation in selection regimes, however, alleles that are deleterious in one sex and neutral in the other can be maintained in other neighboring populations and gene flow may continuously supply deleterious alleles. We hypothesized that this maintenance of genetic variation may promote the establishment of recombination suppression in sex chromosomes even in cases where selection is limited to one sex. Using individual-based simulations, we show that spatial variation in male-limited selection and gene flow can promote the establishment of Y-autosome fusions, a special case of recombination suppression in sex chromosomes. This can be explained by the fact that fused Y-chromosomes that capture alleles that are beneficial for local males have a higher mean fitness compared to unfused Y chromosomes in the presence of deleterious gene flow. We also simulated the case of sex-concordant selection and found that gene flow of alleles that are deleterious in both sexes did not substantially increase the establishment rates of Y-autosome fusions across the parameter space examined. This can be accounted for by the fact that foreign alleles that are deleterious in both sexes can be efficiently removed from the population compared to alleles that are neutral in females. These results indicate that how gene flow affects the establishment rates of Y-autosome fusions depends largely on selection regimes. Spatial variation in sex-specific selection and gene flow should be appreciated as a factor affecting sex chromosome evolution.

Different contributions of local- and distant-regulatory changes to transcriptome divergence between stickleback ecotypes

Differential gene expression can play an important role in phenotypic evolution and divergent adaptation. Although differential gene expression can be caused by both local- and distant-regulatory changes, we know little about their relative contribution to transcriptome evolution in natural populations. Here, we conducted expression quantitative trait loci (eQTL) analysis to investigate the genetic architecture underlying transcriptome divergence between marine and stream ecotypes of threespine sticklebacks (Gasterosteus aculeatus). We identified both local and distant eQTLs, some of which constitute hotspots, regions with a disproportionate number of significant eQTLs relative to the genomic background. The majority of local eQTLs including those in the hotspots caused expression changes consistent with the direction of transcriptomic divergence between ecotypes. Genome scan analysis showed that many local eQTLs overlapped with genomic regions of high differentiation. In contrast, nearly half of the distant eQTLs including those in the hotspots caused opposite expression changes, and few overlapped with regions of high differentiation, indicating that distant eQTLs may act as a constraint of transcriptome evolution. Finally, a comparison between two salinity conditions revealed that nearly half of eQTL hotspots were environment specific, suggesting that analysis of genetic architecture in multiple conditions is essential for predicting response to selection.

Whole-genome sequencing reveals small genomic regions of introgression in an introduced crater lake population of threespine stickleback

Invasive species pose a major threat to biological diversity. Although introduced populations often experience population bottlenecks, some invasive species are thought to be originated from hybridization between multiple populations or species, which can contribute to the maintenance of high genetic diversity. Recent advances in genome sequencing enable us to trace the evolutionary history of invasive species even at whole‐genome level and may help to identify the history of past hybridization that may be overlooked by traditional marker‐based analysis. Here, we conducted whole‐genome sequencing of eight threespine stickleback (Gasterosteus aculeatus) individuals, four from a recently introduced crater lake population and four of the putative source population. We found that both populations have several small genomic regions with high genetic diversity, which resulted from introgression from a closely related species (Gasterosteus nipponicus). The sizes of the regions were too small to be detected with traditional marker‐based analysis or even some reduced‐representation sequencing methods. Further amplicon sequencing revealed linkage disequilibrium around an introgression site, which suggests the possibility of selective sweep at the introgression site. Thus, interspecies introgression might predate introduction and increase genetic variation in the source population. Whole‐genome sequencing of even a small number of individuals can therefore provide higher resolution inference of history of introduced populations.

Purifying Selection Maintains Dosage-Sensitive Genes during Degeneration of the Threespine Stickleback Y Chromosome

Sex chromosomes are subject to unique evolutionary forces that cause suppression of recombination, leading to sequence degeneration and the formation of heteromorphic chromosome pairs (i.e., XY or ZW). Although progress has been made in characterizing the outcomes of these evolutionary processes on vertebrate sex chromosomes, it is still unclear how recombination suppression and sequence divergence typically occur and how gene dosage imbalances are resolved in the heterogametic sex. The threespine stickleback fish (Gasterosteus aculeatus) is a powerful model system to explore vertebrate sex chromosome evolution, as it possesses an XY sex chromosome pair at relatively early stages of differentiation. Using a combination of whole-genome and transcriptome sequencing, we characterized sequence evolution and gene expression across the sex chromosomes. We uncovered two distinct evolutionary strata that correspond with known structural rearrangements on the Y chromosome. In the oldest stratum, only a handful of genes remain, and these genes are under strong purifying selection. By comparing sex-linked gene expression with expression of autosomal orthologs in an outgroup, we show that dosage compensation has not evolved in threespine sticklebacks through upregulation of the X chromosome in males. Instead, in the oldest stratum, the genes that still possess a Y chromosome allele are enriched for genes predicted to be dosage sensitive in mammals and yeast. Our results suggest that dosage imbalances may have been avoided at haploinsufficient genes by retaining function of the Y chromosome allele through strong purifying selection.

Comparative analysis of Japanese three-spined stickleback clades reveals the Pacific Ocean lineage has adapted to freshwater environments while the Japan Sea has not

Divergent selection and adaptive divergence can increase phenotypic diversification amongst populations and lineages. Yet adaptive divergence between different environments, habitats or niches does not occur in all lineages. For example, the colonization of freshwater environments by ancestral marine species has triggered adaptive radiation and phenotypic diversification in some taxa but not in others. Studying closely related lineages differing in their ability to diversify is an excellent means of understanding the factors promoting and constraining adaptive evolution. A well-known example of the evolution of increased phenotypic diversification following freshwater colonization is the three-spined stickleback. Two closely related stickleback lineages, the Pacific Ocean and the Japan Sea occur in Japan. However, Japanese freshwater stickleback populations are derived from the Pacific Ocean lineage only, suggesting the Japan Sea lineage is unable to colonize freshwater. Using stable isotope data and trophic morphology, we first show higher rates of phenotypic and ecological diversification between marine and freshwater populations within the Pacific Ocean lineage, confirming adaptive divergence has occurred between the two lineages and within the Pacific Ocean lineage but not in the Japan Sea lineage. We further identified consistent divergence in diet and foraging behaviour between marine forms from each lineage, confirming Pacific Ocean marine sticklebacks, from which all Japanese freshwater populations are derived, are better adapted to freshwater environments than Japan Sea sticklebacks. We suggest adaptive divergence between ancestral marine populations may have played a role in constraining phenotypic diversification and adaptive evolution in Japanese sticklebacks.