Short Evolutionary Timescales Research Articles

Dynamics and stage-specificity of between-population gene expression divergence in the Drosophila melanogaster larval fat body

Gene expression variation is pervasive across all levels of organismal organization, including development. Few studies, however, have examined variation in developmental transcriptional dynamics among populations, or how it contributes to phenotypic divergence. Indeed, the evolution of gene expression dynamics when both the evolutionary and temporal timescale are comparatively short remains relatively uncharacterized. Here, we examined coding and non-coding gene expression in the fat body of an ancestral African and a derived European Drosophila melanogaster population across three developmental stages spanning ten hours of larval development. Between populations, expression divergence was largely stage-specific. We detected higher expression variation during the late wandering stage, which may be a general feature of this stage. During this stage, we also detected higher and more extensive lncRNA expression in Europe, suggesting that lncRNA expression may be more important in derived populations. Interestingly, the temporal breadth of protein-coding and lncRNA expression became more restricted in the derived population. Taken together with the signatures of potential local adaptation that we detected at the sequence level in 9-25% of candidate genes (those showing evidence of expression divergence between populations), this finding suggests that gene expression becomes more developmental stage-specific during adaptation to new environments. We further used RNAi to identify several candidate genes that likely contribute to known phenotypic divergence between these populations. Our results shed light on the evolution and dynamics of expression variation over short developmental and evolutionary timescales, and how this variation contributes to population and phenotypic divergence.

The Evolution of Mitochondrial Genomes between Two Cymbidium Sister Species: Dozens of Circular Chromosomes and the Maintenance and Deterioration of Genome Synteny.

Plant mitochondrial genomes (mitogenomes) exhibit fluid genome architectures, which could lead to the rapid erosion of genome synteny over a short evolutionary time scale. Among the species-rich orchid family, the leafy Cymbidium lancifolium and leafless Cymbidium macrorhizon are sister species with remarkable differences in morphology and nutritional physiology. Although our understanding of the evolution of mitochondria is incomplete, these sister taxa are ideal for examining this subject. In this study, the complete mitogenomes of C. lancifolium and C. macrorhizon, totaling 704,244 bp and 650,751 bp, respectively, were assembled. In the 2 mitogenomes, 38 protein-coding genes, 18 cis- and 6 trans-spliced introns, and approximately 611 Kb of homologous sequences are identical; overall, they have 99.4% genome-wide similarity. Slight variations in the mitogenomes of C. lancifolium and C. macrorhizon in repeat content (21.0 Kb and 21.6 Kb, respectively) and mitochondrial DNA of plastid origin (MIPT; 38.2 Kb and 37.5 Kb, respectively) were observed. The mitogenome architectures of C. lancifolium and C. macrorhizon are complex and comprise 23 and 22 mini-circular chromosomes, respectively. Pairwise comparisons indicate that the two mitogenomes are largely syntenic, and the disparity in chromosome numbers is likely due to repeat-mediated rearrangements among different chromosomes. Notably, approximately 93.2 Kb C. lancifolium mitochondrial sequences lack any homology in the C. macrorhizon mitogenome, indicating frequent DNA gains and losses, which accounts mainly for the size variation. Our findings provide unique insights into mitogenome evolution in leafy and leafless plants of sister species and shed light on mitogenome dynamics during the transition from mixotrophy to mycoheterotrophy.

Sexually antagonistic selection maintains genetic variance when sexual dimorphism evolves.

Genetic variance (VG) in fitness related traits is often unexpectedly high, evoking the question how VG can be maintained in the face of selection. Sexually antagonistic (SA) selection favouring alternative alleles in the sexes is common and predicted to maintain VG, while directional selection should erode it. Both SA and sex-limited directional selection can lead to sex-specific adaptations but how each affect VG when sexual dimorphism evolves remain experimentally untested. Using replicated artificial selection on the seed beetle Callosobruchus maculatus body size we recently demonstrated an increase in size dimorphism under SA and male-limited (ML) selection by 50% and 32%, respectively. Here we test their consequences on genetic variation. We show that SA selection maintained significantly more ancestral, autosomal additive genetic variance than ML selection, while both eroded sex-linked additive variation equally. Ancestral female-specific dominance variance was completely lost under ML, while SA selection consistently sustained it. Further, both forms of selection preserved a high genetic correlation between the sexes (rm,f). These results demonstrate the potential for sexual antagonism to maintain more genetic variance while fuelling sex-specific adaptation in a short evolutionary time scale, and are in line with predicted importance of sex-specific dominance reducing sexual conflict over alternative alleles.

Evolution of female song and duetting in the chaffinch (Fringilla) species complex

Female song is ancestral to songbirds and shows considerable phylogenetic signal, but its presence also appears to be labile and correlated with life‐history and ecology. While previous studies have examined the evolution of female song across species‐rich families, here we studied female song in island populations of a recently diverged species‐complex, the chaffinches (genus Fringilla). We show this behaviour has recently evolved in these populations, and probably on two independent occasions. In F. teydea, female song is performed in loose duets with males, while in F. coelebs, female song is produced solo. Populations with singing females showed year‐round territoriality and were found in regions with low seasonality – both factors previously connected with high rates of female singing. To determine the relative saliency of female songs to conspecifics, and the degree to which they instigate territorial defence behaviours, we performed a series of speaker playback experiments. In F. c. canariensis, female song could induce comparable responses to male song, whilst duetting F. teydea playback produced similar responses to solo songs, and therefore may relate to within‐pair communication instead. Our results suggest female song can be a highly labile trait that can evolve over short evolutionary timescales.

Evolutionary conservation and diversification of auditory neural circuits that process courtship songs in Drosophila

Acoustic communication signals diversify even on short evolutionary time scales. To understand how the auditory system underlying acoustic communication could evolve, we conducted a systematic comparison of the early stages of the auditory neural circuit involved in song information processing between closely-related fruit-fly species. Male Drosophila melanogaster and D. simulans produce different sound signals during mating rituals, known as courtship songs. Female flies from these species selectively increase their receptivity when they hear songs with conspecific temporal patterns. Here, we firstly confirmed interspecific differences in temporal pattern preferences; D. simulans preferred pulse songs with longer intervals than D. melanogaster. Primary and secondary song-relay neurons, JO neurons and AMMC-B1 neurons, shared similar morphology and neurotransmitters between species. The temporal pattern preferences of AMMC-B1 neurons were also relatively similar between species, with slight but significant differences in their band-pass properties. Although the shift direction of the response property matched that of the behavior, these differences are not large enough to explain behavioral differences in song preferences. This study enhances our understanding of the conservation and diversification of the architecture of the early-stage neural circuit which processes acoustic communication signals.

Testing the success of palaeontological methods in the delimitation of clam shrimp (Crustacea, Branchiopoda) on extant species

AbstractFossil spinicaudatan taxonomy heavily relies on carapace features (size, shape, ornamentation) and palaeontologists have greatly refined methods to study and describe carapace variability. Whether carapace features alone are sufficient for distinguishing between species of a single genus has remained untested. In our study, we tested common palaeontological methods on 481 individuals of the extant Australian genus Ozestheria that have been previously assigned to ten species based on genetic analysis. All species are morphologically distinct based on geometric morphometrics (p ≤ 0.001), but they occupy overlapping regions in Ozestheria morphospace. Linear discriminant analysis of Fourier shape coefficients reaches a mean model performance of 93.8% correctly classified individuals over all possible 45 pairwise species comparisons. This can be further increased by combining the size and shape datasets. Nine of the ten examined species are clearly sexually dimorphic but male and female morphologies strongly overlap within species with little influence on model performance. Ornamentation is commonly species‐diagnostic; seven ornamentation types are distinguished of which six are species‐specific while one is shared by four species. A transformation of main ornamental features (e.g. from punctate to smooth) can occur among closely related species suggesting short evolutionary timescales. Our overall results support the taxonomic value of carapace features, which should also receive greater attention in the taxonomy of extant species. The extensive variation in carapace shape and ornamentation is noteworthy and several species would probably have been assigned to different genera or families if these had been fossils, bearing implications for the systematics of fossil Spinicaudata.

Neural networks: solving the chemistry of the interstellar medium

ABSTRACT Non-equilibrium chemistry is a key process in the study of the interstellar medium (ISM), in particular the formation of molecular clouds and thus stars. However, computationally, it is among the most difficult tasks to include in astrophysical simulations, because of the typically high (>40) number of reactions, the short evolutionary time-scales (about 104 times less than the ISM dynamical time), and the characteristic non-linearity and stiffness of the associated ordinary differential equations system (ODEs). In this proof of concept work, we show that Physics Informed Neural Networks (PINN) are a viable alternative to traditional ODE time integrators for stiff thermochemical systems, i.e. up to molecular hydrogen formation (9 species and 46 reactions). Testing different chemical networks in a wide range of densities (−2 < log n/cm−3 < 3) and temperatures (1 < log T/K < 5), we find that a basic architecture can give a comfortable convergence only for simplified chemical systems: to properly capture the sudden chemical and thermal variations, a Deep Galerkin Method is needed. Once trained (∼103 GPUhr), the PINN well reproduces the strong non-linear nature of the solutions (errors $\lesssim 10{{\ \rm per\ cent}}$) and can give speed-ups up to a factor of ∼200 with respect to traditional ODE solvers. Further, the latter have completion times that vary by about $\sim 30{{\ \rm per\ cent}}$ for different initial n and T, while the PINN method gives negligible variations. Both the speed-up and the potential improvement in load balancing imply that PINN-powered simulations are a very palatable way to solve complex chemical calculation in astrophysical and cosmological problems.

Polygenic plague resistance in the great gerbil uncovered by population sequencing.

Pathogens can elicit high selective pressure on hosts, potentially altering genetic diversity over short evolutionary timescales. Intraspecific variation in immune response is observable as variable survivability from specific infections. The great gerbil (Rhombomys opimus) is a rodent plague host with a heterogenic but highly resistant phenotype. Here, we investigate the genomic basis for plague-resistant phenotypes by exposing wild-caught great gerbils to plague (Yersinia pestis). Whole genome sequencing of 10 survivors and 10 moribund individuals revealed a subset of genomic regions showing elevated differentiation. Gene ontology analysis of candidate genes in these regions demonstrated enrichment of genes directly involved in immune functions, cellular metabolism and the regulation of apoptosis as well as pathways involved in transcription, translation, and gene regulation. Transcriptomic analysis revealed that the early activated great gerbil immune response to plague consisted of classical components of the innate immune system. Our approach combining challenge experiments with transcriptomics and population level sequencing, provides new insight into the genetic background of plague-resistance and confirms its complex nature, most likely involving multiple genes and pathways of both the immune system and regulation of basic cellular functions.

A butterfly pan-genome reveals that a large amount of structural variation underlies the evolution of chromatin accessibility.

Despite insertions and deletions being the most common structural variants (SVs) found across genomes, not much is known about how much these SVs vary within populations and between closely related species, nor their significance in evolution. To address these questions, we characterized the evolution of indel SVs using genome assemblies of three closely related Heliconius butterfly species. Over the relatively short evolutionary timescales investigated, up to 18.0% of the genome was composed of indels between two haplotypes of an individual Heliconius charithonia butterfly and up to 62.7% included lineage-specific SVs between the genomes of the most distant species (11 Mya). Lineage-specific sequences were mostly characterized as transposable elements (TEs) inserted at random throughout the genome and their overall distribution was similarly affected by linked selection as single nucleotide substitutions. Using chromatin accessibility profiles (i.e., ATAC-seq) of head tissue in caterpillars to identify sequences with potential cis-regulatory function, we found that out of the 31,066 identified differences in chromatin accessibility between species, 30.4% were within lineage-specific SVs and 9.4% were characterized as TE insertions. These TE insertions were localized closer to gene transcription start sites than expected at random and were enriched for sites with significant resemblance to several transcription factor binding sites with known function in neuron development in Drosophila We also identified 24 TE insertions with head-specific chromatin accessibility. Our results show high rates of structural genome evolution that were previously overlooked in comparative genomic studies and suggest a high potential for structural variation to serve as raw material for adaptive evolution.

Low protein expression enhances phenotypic evolvability by intensifying selection on folding stability.

Protein abundance affects the evolution of protein genotypes, but we do not know how it affects the evolution of protein phenotypes. Here we investigate the role of protein abundance in the evolvability of green fluorescent protein (GFP) towards the novel phenotype of cyan fluorescence. We evolve GFP in E. coli through multiple cycles of mutation and selection and show that low GFP expression facilitates the evolution of cyan fluorescence. A computational model whose predictions we test experimentally helps explain why: lowly expressed proteins are under stronger selection for proper folding, which facilitates their evolvability on short evolutionary time scales. The reason is that high fluorescence can be achieved by either few proteins that fold well or by many proteins that fold less well. In other words, we observe a synergy between a protein's scarcity and its stability. Because many proteins meet the essential requirements for this scarcity-stability synergy, it may be a widespread mechanism by which low expression helps proteins evolve new phenotypes and functions.

Widespread stop-codon recoding in bacteriophages may regulate translation of lytic genes.

Bacteriophages (phages) are obligate parasites that use host bacterial translation machinery to produce viral proteins. However, some phages have alternative genetic codes with reassigned stop codons that are predicted to be incompatible with bacterial translation systems. We analysed 9,422 phage genomes and found that stop-codon recoding has evolved in diverse clades of phages that infect bacteria present in both human and animal gut microbiota. Recoded stop codons are particularly over-represented in phage structural and lysis genes. We propose that recoded stop codons might function to prevent premature production of late-stage proteins. Stop-codon recoding has evolved several times in closely related lineages, which suggests that adaptive recoding can occur over very short evolutionary timescales.

Adaptive Divergence without Distinct Species Relationships Indicate Early Stage Ecological Speciation in Species of the Rhododendronpseudochrysanthum Complex Endemic to Taiwan.

The testing association of environmental variables with genetic and epigenetic variation could be crucial to deciphering the effects of environmental factors playing roles as selective drivers in ecological speciation. Although ecological speciation may occur in closely related species, species boundaries may not be established over a short evolutionary timescale. Here, we investigated the genetic and epigenetic variations using amplified fragment length polymorphism (AFLP) and methylation-sensitive amplification polymorphism (MSAP), respectively, and tested their associations with environmental variables in populations of four closely related species in the R. pseudochrysanthum complex. No distinctive species relationships were found using genetic clustering analyses, neighbor-joining tree, and neighbor-net tree based on the total AFLP variation, which is suggestive of the incomplete lineage sorting of ancestral variation. Nonetheless, strong isolation-by-environment and adaptive divergence were revealed, despite the significant isolation-by-distance. Annual mean temperature, elevation, normalized difference vegetation index, and annual total potential evapotranspiration were found to be the most important environmental variables explaining outlier genetic and epigenetic variations. Our results suggest that the four closely related species of the R. pseudochrysanthum complex share the polymorphism of their ancestor, but reproductive isolation due to ecological speciation can occur if local environmental divergence persists over time.

Fitness of Isidorella newcombi Following Multi-generational Cu Exposures: Mortality, Cellular Biomarkers and Life History Responses

The effects of multigenerational Cu exposure on the freshwater gastropod Isidorella newcombi were investigated. Snails were exposed to a range of treatment-specific Cu concentrations in the parental to F2 generations, and a common Cu concentration in the F3 generation. In the parental to F2 generations, some general responses to 3 days Cu exposures included reduced survival and feeding in snails exposed to higher Cu concentrations. This suggested that the snails exposed to the high Cu concentration were experiencing Cu-induced stress that may apply selection pressure. In the F3 generation, when all treatments were exposed to a common Cu concentration, increased survival was correlated with the pre-exposure Cu concentration history. Snails that had been pre-exposed to Cu also displayed reduced stress at a sub-lethal level, indicated by lower lysosomal destabilisation (LD). Mortality and LD responses in the F3 generation were not related to Cu tissue concentrations, indicating increased tolerance and reduced stress were not related to changes in Cu bioaccumulation. Total antioxidant capacity increased in the higher Cu concentration pre-exposure treatments which could be associated with lower Cu-induced stress, however, this is not supported by the oxidative damage marker lipid peroxidation, which also increased. While Cu tissue concentrations and oxidative stress markers were assessed to determine underlying reasons for increased tolerance in snails from a population with a multi-generational exposure history to Cu, the results were not conclusive. Despite this, it was demonstrated through increased survival and reduced LD that Cu tolerance can develop over a short evolutionary time scale.

Ancestry adjustment improves genome-wide estimates of regional intolerance.

Genomic regions subject to purifying selection are more likely to carry disease-causing mutations than regions not under selection. Cross species conservation is often used to identify such regions but with limited resolution to detect selection on short evolutionary timescales such as that occurring in only one species. In contrast, genetic intolerance looks for depletion of variation relative to expectation within a species, allowing species-specific features to be identified. When estimating the intolerance of noncoding sequence, methods strongly leverage variant frequency distributions. As the expected distributions depend on ancestry, if not properly controlled for, ancestral population source may obfuscate signals of selection. We demonstrate that properly incorporating ancestry in intolerance estimation greatly improved variant classification. We provide a genome-wide intolerance map that is conditional on ancestry and likely to be particularly valuable for variant prioritization.

Rapid Cis-Trans Coevolution Driven by a Novel Gene Retroposed from a Eukaryotic Conserved CCR4-NOT Component in Drosophila.

Young, or newly evolved, genes arise ubiquitously across the tree of life, and they can rapidly acquire novel functions that influence a diverse array of biological processes. Previous work identified a young regulatory duplicate gene in Drosophila, Zeus that unexpectedly diverged rapidly from its parent, Caf40, an extremely conserved component in the CCR4–NOT machinery in post-transcriptional and post-translational regulation of eukaryotic cells, and took on roles in the male reproductive system. This neofunctionalization was accompanied by differential binding of the Zeus protein to loci throughout the Drosophila melanogaster genome. However, the way in which new DNA-binding proteins acquire and coevolve with their targets in the genome is not understood. Here, by comparing Zeus ChIP-Seq data from D. melanogaster and D. simulans to the ancestral Caf40 binding events from D. yakuba, a species that diverged before the duplication event, we found a dynamic pattern in which Zeus binding rapidly coevolved with a previously unknown DNA motif, which we term Caf40 and Zeus-Associated Motif (CAZAM), under the influence of positive selection. Interestingly, while both copies of Zeus acquired targets at male-biased and testis-specific genes, D. melanogaster and D. simulans proteins have specialized binding on different chromosomes, a pattern echoed in the evolution of the associated motif. Using CRISPR-Cas9-mediated gene knockout of Zeus and RNA-Seq, we found that Zeus regulated the expression of 661 differentially expressed genes (DEGs). Our results suggest that the evolution of young regulatory genes can be coupled to substantial rewiring of the transcriptional networks into which they integrate, even over short evolutionary timescales. Our results thus uncover dynamic genome-wide evolutionary processes associated with new genes.

Quantifying bacterial evolution in the wild: A birthday problem for Campylobacter lineages.

Measuring molecular evolution in bacteria typically requires estimation of the rate at which nucleotide changes accumulate in strains sampled at different times that share a common ancestor. This approach has been useful for dating ecological and evolutionary events that coincide with the emergence of important lineages, such as outbreak strains and obligate human pathogens. However, in multi-host (niche) transmission scenarios, where the pathogen is essentially an opportunistic environmental organism, sampling is often sporadic and rarely reflects the overall population, particularly when concentrated on clinical isolates. This means that approaches that assume recent common ancestry are not applicable. Here we present a new approach to estimate the molecular clock rate in Campylobacter that draws on the popular probability conundrum known as the ‘birthday problem’. Using large genomic datasets and comparative genomic approaches, we use isolate pairs that share recent common ancestry to estimate the rate of nucleotide change for the population. Identifying synonymous and non-synonymous nucleotide changes, both within and outside of recombined regions of the genome, we quantify clock-like diversification to estimate synonymous rates of nucleotide change for the common pathogenic bacteria Campylobacter coli (2.4 x 10−6 s/s/y) and Campylobacter jejuni (3.4 x 10−6 s/s/y). Finally, using estimated total rates of nucleotide change, we infer the number of effective lineages within the sample time frame–analogous to a shared birthday–and assess the rate of turnover of lineages in our sample set over short evolutionary timescales. This provides a generalizable approach to calibrating rates in populations of environmental bacteria and shows that multiple lineages are maintained, implying that large-scale clonal sweeps may take hundreds of years or more in these species.

Autophagy as a Mechanism for Adaptive Prediction-Mediated Emergence of Drug Resistance.

Drug resistance is a major problem in treatment of microbial infections and cancers. There is growing evidence that a transient drug tolerant state may precede and potentiate the emergence of drug resistance. Therefore, understanding the mechanisms leading to tolerance is critical for combating drug resistance and for the development of effective therapeutic strategy. Through laboratory evolution of yeast, we recently demonstrated that adaptive prediction (AP), a strategy employed by organisms to anticipate and prepare for a future stressful environment, can emerge within 100 generations by linking the response triggered by a neutral cue (caffeine) to a mechanism of protection against a lethal agent (5-fluoroorotic acid, 5-FOA). Here, we demonstrate that mutations selected across multiple laboratory-evolved lines had linked the neutral cue response to core genes of autophagy. Across these evolved lines, conditional activation of autophagy through AP conferred tolerance, and potentiated subsequent selection of mutations in genes specific to overcoming the toxicity of 5-FOA. These results offer a new perspective on how extensive genome-wide genetic interactions of autophagy could have facilitated the emergence of AP over short evolutionary timescales to potentiate selection of 5-FOA resistance-conferring mutations.

Exposure to foreign gut microbiota can facilitate rapid dietary shifts

Dietary niche is fundamental for determining species ecology; thus, a detailed understanding of what drives variation in dietary niche is vital for predicting ecological shifts and could have implications for species management. Gut microbiota can be important for determining an organism’s dietary preference, and therefore which food resources they are likely to exploit. Evidence for whether the composition of the gut microbiota is plastic in response to changes in diet is mixed. Also, the extent to which dietary preference can be changed following colonisation by new gut microbiota from different species is unknown. Here, we use Drosophila spp. to show that: (1) the composition of an individual’s gut microbiota can change in response to dietary changes, and (2) ingestion of foreign gut microbes can cause individuals to be attracted to food types they previously had a strong aversion to. Thus, we expose a mechanism for facilitating rapid shifts in dietary niche over short evolutionary timescales.

Mobile Gene Sequence Evolution within Individual Human Gut Microbiomes Is Better Explained by Gene-Specific Than Host-Specific Selective Pressures.

Pangenomes—the cumulative set of genes encoded by a population or species—arise from the interplay of horizontal gene transfer, drift, and selection. The balance of these forces in shaping pangenomes has been debated, and studies to date focused on ancient evolutionary time scales have suggested that pangenomes generally confer niche adaptation to their bacterial hosts. To shed light on pangenome evolution on shorter evolutionary time scales, we inferred the selective pressures acting on mobile genes within individual human microbiomes from 176 Fiji islanders. We mapped metagenomic sequence reads to a set of known mobile genes to identify single nucleotide variants (SNVs) and calculated population genetic metrics to infer deviations from a neutral evolutionary model. We found that mobile gene sequence evolution varied more by gene family than by human social attributes, such as household or village. Patterns of mobile gene sequence evolution could be qualitatively recapitulated with a simple evolutionary simulation without the need to invoke the adaptive value of mobile genes to either bacterial or human hosts. These results stand in contrast with the apparent adaptive value of pangenomes over longer evolutionary time scales. In general, the most highly mobile genes (i.e., those present in more distinct bacterial host genomes) tend to have higher metagenomic read coverage and an excess of low-frequency SNVs, consistent with their rapid spread across multiple bacterial species in the gut. However, a subset of mobile genes—including those involved in defense mechanisms and secondary metabolism—showed a contrasting signature of intermediate-frequency SNVs, indicating species-specific selective pressures or negative frequency-dependent selection on these genes. Together, our evolutionary models and population genetic data show that gene-specific selective pressures predominate over human or bacterial host-specific pressures during the relatively short time scales of a human lifetime.

Heterogeneous Histories of Recombination Suppression on Stickleback Sex Chromosomes.

How consistent are the evolutionary trajectories of sex chromosomes shortly after they form? Insights into the evolution of recombination, differentiation, and degeneration can be provided by comparing closely related species with homologous sex chromosomes. The sex chromosomes of the threespine stickleback (Gasterosteus aculeatus) and its sister species, the Japan Sea stickleback (G. nipponicus), have been well characterized. Little is known, however, about the sex chromosomes of their congener, the blackspotted stickleback (G. wheatlandi). We used pedigrees to obtain experimentally phased whole genome sequences from blackspotted stickleback X and Y chromosomes. Using multispecies gene trees and analysis of shared duplications, we demonstrate that Chromosome 19 is the ancestral sex chromosome and that its oldest stratum evolved in the common ancestor of the genus. After the blackspotted lineage diverged, its sex chromosomes experienced independent and more extensive recombination suppression, greater X–Y differentiation, and a much higher rate of Y degeneration than the other two species. These patterns may result from a smaller effective population size in the blackspotted stickleback. A recent fusion between the ancestral blackspotted stickleback Y chromosome and Chromosome 12, which produced a neo-X and neo-Y, may have been favored by the very small size of the recombining region on the ancestral sex chromosome. We identify six strata on the ancestral and neo-sex chromosomes where recombination between the X and Y ceased at different times. These results confirm that sex chromosomes can evolve large differences within and between species over short evolutionary timescales.