When adaptation is slowed down: Genomic analysis of evolutionary stasis in thermal tolerance during biological invasion in a novel climate

Research conducted during the past two decades has demonstrated that biological invasions are excellent models of rapid evolution. Even so, characteristics of invasive populations such as a short time for recombination to assemble optimal combinations of alleles may occasionally limit adaptation to new environments. Here, we investigated such genetic constraints to adaptation in the invasive brown anole (Anolis sagrei)-a tropical ectotherm that was introduced to the southeastern United States, a region with a much colder climate than in its native Caribbean range. We examined thermal physiology for 30 invasive populations and tested for a climatic cline in cold tolerance. Also, we used genomics to identify mechanisms that may limit adaptation. We found no support for a climatic cline, indicating that thermal tolerance did not shift adaptively. Concomitantly, population genomic results were consistent with the occurrence of recombination cold spots that comprise more than half of the genome and maintain long-range associations among alleles in invasive populations. These genomic regions overlap with both candidate thermal tolerance loci that we identified using a standard genome-wide association test. Moreover, we found that recombination cold spots do not have a large contribution to population differentiation in the invasive range, contrary to observations in the native range. We suggest that limited recombination is constraining the contribution of large swaths of the genome to adaptation in invasive brown anoles. Our study provides an example of evolutionary stasis during invasion and highlights the possibility that reduced recombination occasionally slows down adaptation in invasive populations.

Why do sex chromosomes progressively lose recombination?

Progressive recombination loss is a common feature of sex chromosomes. Yet, the evolutionary drivers of this phenomenon remain a mystery. For decades, differences in trait optima between sexes (sexual antagonism) have been the favoured hypothesis, but convincing evidence is lacking. Recent years have seen a surge of alternative hypotheses to explain progressive extensions and maintenance of recombination suppression: neutral accumulation of sequence divergence, selection of nonrecombining fragments with fewer deleterious mutations than average, sheltering of recessive deleterious mutations by linkage to heterozygous alleles, early evolution of dosage compensation, and constraints on recombination restoration. Here, we explain these recent hypotheses and dissect their assumptions, mechanisms, and predictions. We also review empirical studies that have brought support to the various hypotheses.

Evolution, Inheritance, and Strata Formation of the W Chromosome in Duck (Anas platyrhynchos)

The nonrecombining female-limited W chromosome is predicted to experience unique evolutionary processes. Difficulties in assembling W chromosome sequences have hindered the identification of duck W-linked sequences and their evolutionary footprint. To address this, we conducted three initial contig-level genome assemblies and developed a rigorous pipeline by which to successfully expand the W-linked data set, including 11 known genes and 24 newly identified genes. Our results indicate that the W chromosome expression may not be subject to female-specific selection; a significant convergent pattern of upregulation associated with increased female-specific selection was not detected. The genetic stability of the W chromosome is also reflected in the strong evolutionary correlation between it and the mitochondria; the complete consistency of the cladogram topology constructed from their gene sequences proves the shared maternal coevolution. By detecting the evolutionary trajectories of W-linked sequences, we have found that recombination suppression started in four distinct strata, of which three were conserved across Neognathae. Taken together, our results have revealed a unique evolutionary pattern and an independent stratum evolutionary pattern for sex chromosomes.

The origin and evolution of sex chromosomes, revealed by sequencing of the Silene latifolia female genome

White campion (Silene latifolia, Caryophyllaceae) was the first vascular plant where sex chromosomes were discovered. This species is a classic model for studies on plant sex chromosomes due to presence of large, clearly distinguishable X and Y chromosomes that originated de novo about 11 million years ago (mya), but lack of genomic resources for this relatively large genome (∼2.8 Gb) remains a significant hurdle. Here we report S. latifolia female genome assembly integrated with sex-specific genetic maps of this species, focusing on sex chromosomes and their evolution. The analysis reveals a highly heterogeneous recombination landscape with strong reduction in recombination rate in the central parts of all chromosomes. Recombination on the X chromosome in female meiosis primarily occurs at the very ends, and over 85% of the X chromosome length is located in a massive (∼330 Mb) gene-poor, rarely recombining pericentromeric region (Xpr). The results indicate that the non-recombining region on the Y chromosome (NRY) initially evolved in a relatively small (∼15 Mb), actively recombining region at the end of the q-arm, possibly as a result of inversion on the nascent X chromosome. The NRY expanded about 6 mya via linkage between the Xpr and the sex-determining region, which may have been caused by expanding pericentromeric recombination suppression on the X chromosome. These findings shed light on the origin of sex chromosomes in S. latifolia and yield genomic resources to assist ongoing and future investigations into sex chromosome evolution.

Onset and stepwise extensions of recombination suppression are common in mating-type chromosomes of Microbotryum anther-smut fungi

Sex chromosomes and mating-type chromosomes can display large genomic regions without recombination. Recombination suppression often extended stepwise with time away from the sex- or mating-type-determining genes, generating evolutionary strata of differentiation between alternative sex or mating-type chromosomes. In anther-smut fungi of the Microbotryum genus, recombination suppression evolved repeatedly, linking the two mating-type loci and extended multiple times in regions distal to the mating-type genes. Here, we obtained high-quality genome assemblies of alternative mating types for four Microbotryum fungi. We found an additional event of independent chromosomal rearrangements bringing the two mating-type loci on the same chromosome followed by recombination suppression linking them. We also found, in a new clade analysed here, that recombination suppression between the two mating-type loci occurred in several steps, with first an ancestral recombination suppression between one of the mating-type locus and its centromere; later, completion of recombination suppression up to the second mating-type locus occurred independently in three species. The estimated dates of recombination suppression between the mating-type loci ranged from 0.15 to 3.58 million years ago. In total, this makes at least nine independent events of linkage between the mating-type loci across the Microbotryum genus. Several mating-type locus linkage events occurred through the same types of chromosomal rearrangements, where similar chromosome fissions at centromeres represent convergence in the genomic changes leading to the phenotypic convergence. These findings further highlight Microbotryum fungi as excellent models to study the evolution of recombination suppression.

New insights into Xenopus sex chromosome genomics from the Marsabit clawed frog X. borealis

In many groups, sex chromosomes change frequently but the drivers of their rapid evolution are varied and often poorly characterized. With an aim of further understanding sex chromosome turnover, we investigated the polymorphic sex chromosomes of the Marsabit clawed frog, Xenopus borealis, using genomic data and a new chromosome-scale genome assembly. We confirmed previous findings that 54.1 Mb of chromosome 8L is sex-linked in animals from east Kenya and a laboratory strain, but most (or all) of this region is not sex-linked in natural populations from west Kenya. Previous work suggests possible degeneration of the Z chromosomes in the east population because many sex-linked transcripts of this female heterogametic population have female-biased expression, and we therefore expected this chromosome to not be present in the west population. In contrast, our simulations support a model where most or all of the sex-linked portion of the Z chromosome from the east acquired autosomal segregation in the west, and where much genetic variation specific to the large sex-linked portion of the W chromosome from the east is not present in the west. These recent changes are consistent with the hot-potato model, wherein sex chromosome turnover is favoured by natural selection if it purges a (minimally) degenerate sex-specific sex chromosome, but counterintuitively suggest natural selection failed to purge a Z chromosome that has signs of more advanced and possibly more ancient regulatory degeneration. These findings highlight complex evolutionary dynamics of young, rapidly evolving Xenopus sex chromosomes and set the stage for mechanistic work aimed at pinpointing additional sex-determining genes in this group.

Genomic analyses of the Linum distyly supergene reveal convergent evolution at the molecular level

Supergenes govern multi-trait-balanced polymorphisms in a wide range of systems; however, our understanding of their origins and evolution remains incomplete. The reciprocal placement of stigmas and anthers in pin and thrum floral morphs of distylous species constitutes an iconic example of a balanced polymorphism governed by a supergene, the distyly S-locus. Recent studies have shown that the Primula and Turnera distyly supergenes are both hemizygous in thrums, but it remains unknown whether hemizygosity is pervasive among distyly S-loci. As hemizygosity has major consequences for supergene evolution and loss, clarifying whether this genetic architecture is shared among distylous species is critical. Here, we have characterized the genetic architecture and evolution of the distyly supergene in Linum by generating a chromosome-level genome assembly of Linum tenue, followed by the identification of the S-locus using population genomic data. We show that hemizygosity and thrum-specific expression of S-linked genes, including a pistil-expressed candidate gene for style length, are major features of the Linum S-locus. Structural variation is likely instrumental for recombination suppression, and although the non-recombining dominant haplotype has accumulated transposable elements, S-linked genes are not under relaxed purifying selection. Our findings reveal remarkable convergence in the genetic architecture and evolution of independently derived distyly supergenes, provide a counterexample to classic inversion-based supergenes, and shed new light on the origin and maintenance of an iconic floral polymorphism.

Stepwise evolution of a butterfly supergene via duplication and inversion

Supergenes maintain adaptive clusters of alleles in the face of genetic mixing. Although usually attributed to inversions, supergenes can be complex, and reconstructing the precise processes that led to recombination suppression and their timing is challenging. We investigated the origin of the BC supergene, which controls variation in warning coloration in the African monarch butterfly, Danaus chrysippus. By generating chromosome-scale assemblies for all three alleles, we identified multiple structural differences. Most strikingly, we find that a region of more than 1 million bp underwent several segmental duplications at least 7.5 Ma. The resulting duplicated fragments appear to have triggered four inversions in surrounding parts of the chromosome, resulting in stepwise growth of the region of suppressed recombination. Phylogenies for the inversions are incongruent with the species tree and suggest that structural polymorphisms have persisted for at least 4.1 Myr. In addition to the role of duplications in triggering inversions, our results suggest a previously undescribed mechanism of recombination suppression through independent losses of divergent duplicated tracts. Overall, our findings add support for a stepwise model of supergene evolution involving a variety of structural changes.

This article is part of the theme issue ‘Genomic architecture of supergenes: causes and evolutionary consequences’.

The rattlesnake W chromosome: A GC-rich retroelement refugium with retained gene function across ancient evolutionary strata

Sex chromosomes diverge after the establishment of recombination suppression, resulting in differential sex-linkage of genes involved in genetic sex determination and dimorphic traits. This process produces systems of male or female heterogamety wherein the Y and W chromosomes are only present in one sex and are often highly degenerated. Sex-limited Y and W chromosomes contain valuable information about the evolutionary transition from autosomes to sex chromosomes, yet detailed characterizations of the structure, composition, and gene content of sex-limited chromosomes are lacking for many species. In this study, we characterize the female-specific W chromosome of the prairie rattlesnake (Crotalus viridis) and evaluate how recombination suppression and other processes have shaped sex chromosome evolution in ZW snakes. Our analyses indicate that the rattlesnake W chromosome is over 80% repetitive and that an abundance of GC-rich mdg4 elements has driven an overall high degree of GC-richness despite a lack of recombination. The W chromosome is also highly enriched for repeat sequences derived from endogenous retroviruses and likely acts as a “refugium” for these and other retroelements. We annotated 219 putatively functional W-linked genes across at least two evolutionary strata identified based on estimates of sequence divergence between Z and W gametologs. The youngest of these strata is relatively gene-rich, however gene expression across strata suggests retained gene function amidst a greater degree of degeneration following ancient recombination suppression. Functional annotation of W-linked genes indicates a specialization of the W chromosome for reproductive and developmental function since recombination suppression from the Z chromosome.

Recent expansion of the non-recombining sex-linked region on Silene latifolia sex chromosomes

Evolution of a non-recombining sex-specific region on the Y (or W) chromosome (NRY) is a key step in sex chromosome evolution, but how recombination suppression evolves is not well understood. Studies in many different organisms indicated that NRY evolution often involves several expansion steps. Why such NRY expansions occur remains unclear, although it is though that they are likely driven by sexually antagonistic selection. This paper describes a recent NRY expansion due to shift of the pseudoautosomal boundary on the sex chromosomes of a dioecious plant Silene latifolia. The shift resulted in inclusion of at least 16 pseudoautosomal genes into the NRY. This region is pseudoautosomal in closely related Silene dioica and Silene diclinis, indicating that the NRY expansion occurred in S. latifolia after it speciated from the other species ~120 thousand years ago. As S. latifolia and S. dioica actively hybridise across Europe, interspecific gene flow could blur the PAR boundary in these species. The pseudoautosomal genes have significantly elevated genetic diversity (π ~ 3% at synonymous sites), which is consistent with balancing selection maintaining diversity in this region. The recent shift of the PAR boundary in S. latifolia offers an opportunity to study the process of on-going NRY expansion.

Consequences of partially recessive deleterious genetic variation for the evolution of inversions suppressing recombination between sex chromosomes

The evolution of suppressed recombination between sex chromosomes is widely hypothesized to be driven by sexually antagonistic selection (SA), where tighter linkage between the sex-determining gene(s) and nearby SA loci is favored when it couples male-beneficial alleles to the proto-Y chromosome, and female-beneficial alleles to the proto-X. Despite limited empirical evidence, the SA selection hypothesis overshadows several alternatives, including an incomplete but often-repeated "sheltering hypothesis" that suggests that expansion of the sex-linked region (SLR) reduces homozygous expression of partially recessive deleterious mutations at selected loci. Here, we use population genetic models to evaluate the consequences of deleterious mutational variation for the evolution of neutral chromosomal inversions expanding the SLR on proto-Y chromosomes. We find that SLR-expanding inversions face a race against time: lightly loaded inversions are initially beneficial, but eventually become deleterious as they accumulate new mutations, and must fix before this window of opportunity closes. The outcome of this race is strongly influenced by inversion size, the mutation rate, and the dominance coefficient of deleterious mutations. Yet, small inversions have elevated fixation probabilities relative to neutral expectations for biologically plausible parameter values. Our results demonstrate that deleterious genetic variation can plausibly drive recombination suppression in small steps and would be most consistent with empirical patterns of small evolutionary strata or gradual recombination arrest.

Some thoughts about the words we use for thinking about sex chromosome evolution

Sex chromosomes are familiar to most biologists since they first learned about genetics. However, research over the past 100 years has revealed that different organisms have evolved sex-determining systems independently. The differences in the ages of systems, and in how they evolved, both affect whether sex chromosomes have evolved. However, the diversity means that the terminology used tends to emphasize either the similarities or the differences, sometimes causing misunderstandings. In this article, I discuss some concepts where special care is needed with terminology. The following four terms regularly create problems: ‘sex chromosome’, ‘master sex-determining gene’, ‘evolutionary strata’ and ‘genetic degeneration’. There is no generally correct or wrong use of these words, but efforts are necessary to make clear how they are to be understood in specific situations. I briefly outline some widely accepted ideas about sex chromosomes, and then discuss these ‘problem terms’, highlighting some examples where careful use of the words helps bring to light current uncertainties and interesting questions for future work.

This article is part of the theme issue ‘Sex determination and sex chromosome evolution in land plants’.

Dissection of Structural Reorganization of Wheat 5B Chromosome Associated With Interspecies Recombination Suppression

Chromosomal rearrangements that lead to recombination suppression can have a significant impact on speciation, and they are also important for breeding. The regions of recombination suppression in wheat chromosome 5B were identified based on comparisons of the 5B map of a cross between the Chinese Spring (CS) variety of hexaploid wheat and CS-5Bdic (genotype CS with 5B substituted with its homologue from tetraploid Triticum dicoccoides) with several 5B maps of tetraploid and hexaploid wheat. In total, two regions were selected in which recombination suppression occurred in cross CS × CS-5Bdic when compared with other maps: one on the short arm, 5BS_RS, limited by markers BS00009810/BS00022336, and the second on the long arm, 5BL_RS, between markers Ra_c10633_2155 and BS00087043. The regions marked as 5BS_RS and 5BL_RS, with lengths of 5 Mb and 3.6 Mb, respectively, were mined from the 5B pseudomolecule of CS and compared to the homoeologous regions (7.6 and 3.8 Mb, respectively) of the 5B pseudomolecule of Zavitan (T. dicoccoides). It was shown that, in the case of 5BS_RS, the local heterochromatin islands determined by the satellite DNA (119.2) and transposable element arrays, as well as the dissimilarity caused by large insertions/deletions (chromosome rearrangements) between 5BSs aestivum/dicoccoides, are likely the key determinants of recombination suppression in the region. Two major and two minor segments with significant loss of similarity were recognized within the 5BL_RS region. It was shown that the loss of similarity, which can lead to suppression of recombination in the 5BL_RS region, is caused by chromosomal rearrangements, driven by the activity of mobile genetic elements (both DNA transposons and long terminal repeat retrotransposons) and their divergence during evolution. It was noted that the regions marked as 5BS_RS and 5BL_RS are associated with chromosomal rearrangements identified earlier by С-banding analysis of intraspecific polymorphism of tetraploid emmer wheat. The revealed divergence in 5BS_RS and 5BL_RS may be a consequence of interspecific hybridization, plant genetic adaptation, or both.

Tempo of degeneration across independently evolved non-recombining regions

Recombination is beneficial over the long term, allowing more effective selection. Despite long-term advantages of recombination, local recombination suppression can evolve and lead to genomic degeneration, in particular on sex chromosomes. Here, we investigated the tempo of degeneration in nonrecombining regions, that is, the function curve for the accumulation of deleterious mutations over time, leveraging on 22 independent events of recombination suppression identified on mating-type chromosomes of anther-smut fungi, including newly identified ones. Using previously available and newly generated high-quality genome assemblies of alternative mating types of 13 Microbotryum species, we estimated degeneration levels in terms of accumulation of nonoptimal codons and nonsynonymous substitutions in nonrecombining regions. We found a reduced frequency of optimal codons in the nonrecombining regions compared with autosomes, that was not due to less frequent GC-biased gene conversion or lower ancestral expression levels compared with recombining regions. The frequency of optimal codons rapidly decreased following recombination suppression and reached an asymptote after ca. 3 Ma. The strength of purifying selection remained virtually constant at d_N/d_S = 0.55, that is, at an intermediate level between purifying selection and neutral evolution. Accordingly, nonsynonymous differences between mating-type chromosomes increased linearly with stratum age, at a rate of 0.015 per My. We thus develop a method for disentangling effects of reduced selection efficacy from GC-biased gene conversion in the evolution of codon usage and we quantify the tempo of degeneration in nonrecombining regions, which is important for our knowledge on genomic evolution and on the maintenance of regions without recombination.

Evolutionary Genomics of Sex-Related Chromosomes at the Base of the Green Lineage

Although sex is now accepted as a ubiquitous and ancestral feature of eukaryotes, direct observation of sex is still lacking in most unicellular eukaryotic lineages. Evidence of sex is frequently indirect and inferred from the identification of genes involved in meiosis from whole genome data and/or the detection of recombination signatures from genetic diversity in natural populations. In haploid unicellular eukaryotes, sex-related chromosomes are named mating-type (MTs) chromosomes and generally carry large genomic regions where recombination is suppressed. These regions have been characterized in Fungi and Chlorophyta and determine gamete compatibility and fusion. Two candidate MT+ and MT− alleles, spanning 450–650 kb, have recently been described in Ostreococcus tauri, a marine phytoplanktonic alga from the Mamiellophyceae class, an early diverging branch in the green lineage. Here, we investigate the architecture and evolution of these candidate MT+ and MT− alleles. We analyzed the phylogenetic profile and GC content of MT gene families in eight different genomes whose divergence has been previously estimated at up to 640 Myr, and found evidence that the divergence of the two MT alleles predates speciation in the Ostreococcus genus. Phylogenetic profiles of MT trans-specific polymorphisms in gametologs disclosed candidate MTs in two additional species, and possibly a third. These Mamiellales MT candidates are likely to be the oldest mating-type loci described to date, which makes them fascinating models to investigate the evolutionary mechanisms of haploid sex determination in eukaryotes.

Compatibility between Solubility and Enhanced Crystallinity of Benzotriazole-Based Small Molecular Acceptors with Less Bulky Alkyl Chains for Organic Solar Cells

Optimizing the molecular structures of organic solar cell (OSC) materials and boosting the power conversion efficiencies are the eternal theme in the solar energy region. A series of fused benzotriazole (BTA)-based A-DA'D-A structures of nonfullerene acceptors (such as Y18) were developed for application in efficient OSCs, in which high quantum efficiencies and low voltage losses could be achieved because of the optimized electron-deficient core and specific molecular geometry. Here, based on the BTA core, the bulky alkyl chain on the BTA unit was further tailored to minimize the lateral alkyl chains and enhance the crystallinity while maintaining an adequate solubility. The resulting NFAs of BTA-C1, BTA-C5, and BTA-C6 are synthesized. Compared with the well-designed molecular Y18 (BTA-C8), we found that simply replacing the 2-ethylhexyl chain with a single methyl (BTA-C1) can easily improve the fill factor up to 77%, but its poor light absorption capacity and large domain size impeded further efficiency improvement. In particular, the BTA-C5, with a shortened branch alkyl chain of 2-methylbutyl, achieves suitable solubility and enhanced crystallinity. Significantly, owing to the balanced charge carrier mobility and suitable phase separation, the BTA-C5-based binary single-junction OSCs achieve a high efficiency of 17.11%, which is one of the top values in BTA-based OSCs.

Mating-Type Locus Organization and Mating-Type Chromosome Differentiation in the Bipolar Edible Button Mushroom Agaricus bisporus

In heterothallic basidiomycete fungi, sexual compatibility is restricted by mating types, typically controlled by two loci: PR, encoding pheromone precursors and pheromone receptors, and HD, encoding two types of homeodomain transcription factors. We analysed the single mating-type locus of the commercial button mushroom variety, Agaricus bisporus var. bisporus, and of the related variety burnettii. We identified the location of the mating-type locus using genetic map and genome information, corresponding to the HD locus, the PR locus having lost its mating-type role. We found the mip1 and β-fg genes flanking the HD genes as in several Agaricomycetes, two copies of the β-fg gene, an additional HD2 copy in the reference genome of A. bisporus var. bisporus and an additional HD1 copy in the reference genome of A. bisporus var. burnettii. We detected a 140 kb-long inversion between mating types in an A. bisporus var. burnettii heterokaryon, trapping the HD genes, the mip1 gene and fragments of additional genes. The two varieties had islands of transposable elements at the mating-type locus, spanning 35 kb in the A. bisporus var. burnettii reference genome. Linkage analyses showed a region with low recombination in the mating-type locus region in the A. bisporus var. burnettii variety. We found high differentiation between β-fg alleles in both varieties, indicating an ancient event of recombination suppression, followed more recently by a suppression of recombination at the mip1 gene through the inversion in A. bisporus var. burnettii and a suppression of recombination across whole chromosomes in A. bisporus var. bisporus, constituting stepwise recombination suppression as in many other mating-type chromosomes and sex chromosomes.

Sex chromosome degeneration, turnover, and sex-biased expression of sex-linked transcripts in African clawed frogs (Xenopus)

The tempo of sex chromosome evolution-how quickly, in what order, why and how their particular characteristics emerge during evolution-remains poorly understood. To understand this further, we studied three closely related species of African clawed frog (genus Xenopus), that each has independently evolved sex chromosomes. We identified population polymorphism in the extent of sex chromosome differentiation in wild-caught Xenopus borealis that corresponds to a large, previously identified region of recombination suppression. This large sex-linked region of X. borealis has an extreme concentration of genes that encode transcripts with sex-biased expression, and we recovered similar findings in the smaller sex-linked regions of Xenopus laevis and Xenopus tropicalis. In two of these species, strong skews in expression (mostly female-biased in X. borealis, mostly male-biased in X. tropicalis) are consistent with expectations associated with recombination suppression, and in X. borealis, we hypothesize that a degenerate ancestral Y-chromosome transitioned into its contemporary Z-chromosome. These findings indicate that Xenopus species are tolerant of differences between the sexes in dosage of the products of multiple genes, and offer insights into how evolutionary transformations of ancestral sex chromosomes carry forward to affect the function of new sex chromosomes. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part I)'.

Evidence of Recombination Suppression Blocks on the Y Chromosome of Date Palm (Phoenix dactylifera)

The genus Phoenix includes the fruit producing date palm tree among 14 species that are all dioecious. Females produce the fruit that are high in sugar content and used in multiple countries ranging from North Africa to South Asia, especially from the Phoenix dactylifera, Phoenix sylvestris, and Phoenix canariensis species. While females produce the fruit, understanding of the genetic basis of sex control only began recently. Through genus-wide sequencing of males and females we recently identified three genes that are conserved in all males and absent in all females of the genus and confirmed an XY sex chromosome system. While our previous study focused on conservation of male-specific sequences at the genus-level, it would be of interest to better understand the spread of male-specific sequences away from the core conserved male genes on the Y chromosome during speciation. To this end, we enumerated male-specific 16 bp sequences using three male/female pairs from the western subpopulation of date palm and documented the density of these sequences in contigs of a phased date palm genome assembly. Here we show that male specific sequences in the date palm Y chromosome have likely spread in defined events that appear as blocks of varying density with significant changes in density between them. Collinearity of genes in these blocks with oil palm shows high synteny with chromosome 10 between megabase 15 and 23 and reveals that large sections of the date palm Y chromosome have maintained the ancestral structure even as recombination has stopped between X and Y.

The Genomic Architecture and Evolutionary Fates of Supergenes

Supergenes are genomic regions containing sets of tightly linked loci that control multi-trait phenotypic polymorphisms under balancing selection. Recent advances in genomics have uncovered significant variation in both the genomic architecture as well as the mode of origin of supergenes across diverse organismal systems. Although the role of genomic architecture for the origin of supergenes has been much discussed, differences in the genomic architecture also subsequently affect the evolutionary trajectory of supergenes and the rate of degeneration of supergene haplotypes. In this review, we synthesize recent genomic work and historical models of supergene evolution, highlighting how the genomic architecture of supergenes affects their evolutionary fate. We discuss how recent findings on classic supergenes involved in governing ant colony social form, mimicry in butterflies, and heterostyly in flowering plants relate to theoretical expectations. Furthermore, we use forward simulations to demonstrate that differences in genomic architecture affect the degeneration of supergenes. Finally, we discuss implications of the evolution of supergene haplotypes for the long-term fate of balanced polymorphisms governed by supergenes.

Divergence and Remarkable Diversity of the Y Chromosome in Guppies

The guppy sex chromosomes show an extraordinary diversity in divergence across populations and closely related species. In order to understand the dynamics of the guppy Y chromosome, we used linked-read sequencing to assess Y chromosome evolution and diversity across upstream and downstream population pairs that vary in predator and food abundance in three replicate watersheds. Based on our population-specific genome assemblies, we first confirmed and extended earlier reports of two strata on the guppy sex chromosomes. Stratum I shows significant accumulation of male-specific sequence, consistent with Y divergence, and predates the colonization of Trinidad. In contrast, Stratum II shows divergence from the X, but no Y-specific sequence, and this divergence is greater in three replicate upstream populations compared with their downstream pair. Despite longstanding assumptions that sex chromosome recombination suppression is achieved through inversions, we find no evidence of inversions associated with either Stratum I or Stratum II. Instead, we observe a remarkable diversity in Y chromosome haplotypes within each population, even in the ancestral Stratum I. This diversity is likely due to gradual mechanisms of recombination suppression, which, unlike an inversion, allow for the maintenance of multiple haplotypes. In addition, we show that this Y diversity is dominated by low-frequency haplotypes segregating in the population, suggesting a link between haplotype diversity and female preference for rare Y-linked color variation. Our results reveal the complex interplay between recombination suppression and Y chromosome divergence at the earliest stages of sex chromosome divergence.

Evolutionary stasis of the pseudoautosomal boundary in strepsirrhine primates

Sex chromosomes are typically comprised of a non-recombining region and a recombining pseudoautosomal region. Accurately quantifying the relative size of these regions is critical for sex-chromosome biology both from a functional and evolutionary perspective. The evolution of the pseudoautosomal boundary (PAB) is well documented in haplorrhines (apes and monkeys) but not in strepsirrhines (lemurs and lorises). Here, we studied the PAB of seven species representing the main strepsirrhine lineages by sequencing a male and a female genome in each species and using sex differences in coverage to identify the PAB. We found that during primate evolution, the PAB has remained unchanged in strepsirrhines whereas several recombination suppression events moved the PAB and shortened the pseudoautosomal region in haplorrhines. Strepsirrhines are well known to have much lower sexual dimorphism than haplorrhines. We suggest that mutations with antagonistic effects between males and females have driven recombination suppression and PAB evolution in haplorrhines.

Genetic diversity, demographic history and neo-sex chromosomes in the Critically Endangered Raso lark

Small effective population sizes could expose island species to inbreeding and loss of genetic variation. Here, we investigate factors shaping genetic diversity in the Raso lark, which has been restricted to a single islet for approximately 500 years, with a population size of a few hundred. We assembled a reference genome for the related Eurasian skylark and then assessed diversity and demographic history using RAD-seq data (75 samples from Raso larks and two related mainland species). We first identify broad tracts of suppressed recombination in females, indicating enlarged neo-sex chromosomes. We then show that genetic diversity across autosomes in the Raso lark is lower than in its mainland relatives, but inconsistent with long-term persistence at its current population size. Finally, we find that genetic signatures of the recent population contraction are overshadowed by an ancient expansion and persistence of a very large population until the human settlement of Cape Verde. Our findings show how genome-wide approaches to study endangered species can help avoid confounding effects of genome architecture on diversity estimates, and how present-day diversity can be shaped by ancient demographic events.

The Origin and Evolution of Chromosomal Reciprocal Translocation in Quasipaa boulengeri (Anura, Dicroglossidae)

Chromosomal rearrangements have long fascinated evolutionary biologists for being widely implicated in causing genetic differentiation. Suppressed recombination has been demonstrated in various species with inversion; however, there is controversy over whether such recombination suppression would facilitate divergence in reciprocal translocation with reduced fitness. In this study, we used the spiny frog, Quasipaa boulengeri, whose western Sichuan Basin populations exhibit translocation polymorphisms, to test whether the genetic markers on translocated (rearranged) or normal chromosomes have driven this genetic differentiation. We also investigated its overall genetic structure and the possibility of chromosomal fixation. Whole-chromosome painting and genetic structure clustering suggested a single origin of the translocation polymorphisms, and high-throughput sequencing of rearranged chromosomes isolated many markers with known localizations on chromosomes. Using these markers, distinct patterns of gene flow were found between rearranged and normal chromosomes. Genetic differentiation was only found in the translocated chromosomes, not in normal chromosomes or the mitochondrial genome. Hybrid unfitness cannot explain the genetic differentiation, as then the differentiation would be observed throughout the whole genome. Our results suggest that suppressed recombination drives genetic differentiation into a balanced chromosomal polymorphism. Mapping to a reference genome, we found that the region of genetic differentiation covered a wide range of translocated chromosomes, not only in the vicinity of chromosomal breakpoints. Our results imply that the suppressed recombination region could be extended by accumulation of repetitive sequences or capture of alleles that are adapted to the local environment, following the spread and/or fixation of chromosomal rearrangement.

Multiple chromosomal rearrangements in a hybrid zone between Littorina saxatilis ecotypes

Both classical and recent studies suggest that chromosomal inversion polymorphisms are important in adaptation and speciation. However, biases in discovery and reporting of inversions make it difficult to assess their prevalence and biological importance. Here, we use an approach based on linkage disequilibrium among markers genotyped for samples collected across a transect between contrasting habitats to detect chromosomal rearrangements de novo. We report 17 polymorphic rearrangements in a single locality for the coastal marine snail, Littorina saxatilis. Patterns of diversity in the field and of recombination in controlled crosses provide strong evidence that at least the majority of these rearrangements are inversions. Most show clinal changes in frequency between habitats, suggestive of divergent selection, but only one appears to be fixed for different arrangements in the two habitats. Consistent with widespread evidence for balancing selection on inversion polymorphisms, we argue that a combination of heterosis and divergent selection can explain the observed patterns and should be considered in other systems spanning environmental gradients.

LTR Retrotransposons Show Low Levels of Unequal Recombination and High Rates of Intraelement Gene Conversion in Large Plant Genomes

The accumulation and removal of transposable elements (TEs) is a major driver of genome size evolution in eukaryotes. In plants, long terminal repeat (LTR) retrotransposons (LTR-RTs) represent the majority of TEs and form most of the nuclear DNA in large genomes. Unequal recombination (UR) between LTRs leads to removal of intervening sequence and formation of solo-LTRs. UR is a major mechanism of LTR-RT removal in many angiosperms, but our understanding of LTR-RT-associated recombination within the large, LTR-RT-rich genomes of conifers is quite limited. We employ a novel read-based methodology to estimate the relative rates of LTR-RT-associated UR within the genomes of four conifer and seven angiosperm species. We found the lowest rates of UR in the largest genomes studied, conifers and the angiosperm maize. Recombination may also resolve as gene conversion, which does not remove sequence, so we analyzed LTR-RT-associated gene conversion events (GCEs) in Norway spruce and six angiosperms. Opposite the trend for UR, we found the highest rates of GCEs in Norway spruce and maize. Unlike previous work in angiosperms, we found no evidence that rates of UR correlate with retroelement structural features in the conifers, suggesting that another process is suppressing UR in these species. Recent results from diverse eukaryotes indicate that heterochromatin affects the resolution of recombination, by favoring gene conversion over crossing-over, similar to our observation of opposed rates of UR and GCEs. Control of LTR-RT proliferation via formation of heterochromatin would be a likely step toward large genomes in eukaryotes carrying high LTR-RT content.

Recombination Suppression in PbS Quantum Dot Heterojunction Solar Cells by Energy-Level Alignment in the Quantum Dot Active Layers

Using spatial energy-level gradient engineering with quantum dots (QDs) of different sizes to increase the generated carrier collection at the junction of a QD heterojunction solar cell (QDHSC) is a hopeful route for improving the energy-conversion efficiency. However, the results of current related research have shown that a variable band-gap structure in a QDHSC will create an appreciable increase, not in the illumination current density, but rather in the fill factor. In addition, there are a lack of studies on the mechanism of the effect of these graded structures on the photovoltaic performance of QDHSCs. This study presents the development of air atmosphere solution-processed TiO₂/PbS QDs/Au QDHSCs by engineering the energy-level alignment (ELA) of the active layer via the use of a sorted order of differently sized QD layers (four QD sizes). In comparison to the ungraded device (without the ELA), the optimized graded architecture (containing the ELA) solar cells exhibited a great increase (21.4%) in short-circuit current density ( J_sc). As a result, a J_sc value greater than 30 mA/cm² has been realized in planar, thinner absorption layer (∼300 nm) PbS QDHSCs, and the open-circuit voltage ( V_oc) and power-conversion efficiency (PCE) were also improved. Through characterization by the light intensity dependences of the J_sc and V_oc and transient photovoltage decay, we find that (i) the ELA structure, serving as an electron-blocking layer, reduces the interfacial recombination at the PbS/anode interface, and (ii) the ELA structure can drive more carriers toward the desirable collection electrode, and the additional carriers can fill the trap states, reducing the trap-assisted recombination in the PbS QDHSCs. This work has clearly elucidated the mechanism of the recombination suppression in the graded QDHSCs and demonstrated the effects of ELA structure on the improvement of J_sc. The charge recombination mechanisms characterized in this work would be able to shed light on further improvements of QDHSCs, which could even benefit other types of solar cells.

Microsatellites mapping for non-model species with chromosomal rearrangement: a case study in the frog Quasipaa boulengeri (Anura: Dicroglossidae)

Gene mapping is an important resource for understanding the evolution of genes and cytogenetics. Model species with a known genetic map or genome sequence allow for the selection of genetic markers on a desired chromosome, while it is hard to locate these markers on chromosomes of non-model species without such references. A frog species, Quasipaa boulengeri, shows chromosomal rearrangement polymorphisms, making itself a fascinating model for chromosomal speciation mediated by suppressed recombination. However, no markers have been located on its rearranged chromosomes. We present a complete protocol to map microsatellites based on mechanical microdissection and chromosome amplification techniques. Following this protocol, we mapped 71 microsatellites of Q. boulengeri at the chromosome level. In total, eight loci were assigned to rearranged chromosomes, and the other 63 loci might attach to other chromosomes. These microsatellites could be used to compare the gene flow and verify the chromosomal suppressed recombination hypothesis in Q. boulengeri. This integrated protocol could be effectively used to map genes to chromosomes for non-model species.

Sequence diversity patterns suggesting balancing selection in partially sex-linked genes of the plant Silene latifolia are not generated by demographic history or gene flow

DNA sequence diversity in genes in the partially sex-linked pseudoautosomal region (PAR) of the sex chromosomes of the plant Silene latifolia is higher than expected from within-species diversity of other genes. This could be the footprint of sexually antagonistic (SA) alleles that are maintained by balancing selection in a PAR gene (or genes) and affect polymorphism in linked genome regions. SA selection is predicted to occur during sex chromosome evolution, but it is important to test whether the unexpectedly high sequence polymorphism could be explained without it, purely by the combined effects of partial linkage with the sex-determining region and the population's demographic history, including possible introgression from Silene dioica. To test this, we applied approximate Bayesian computation-based model choice to autosomal sequence diversity data, to find the most plausible scenario for the recent history of S. latifolia and then to estimate the posterior density of the most relevant parameters. We then used these densities to simulate variation to be expected at PAR genes. We conclude that an excess of variants at high frequencies at PAR genes should arise in S. latifolia populations only for genes with strong associations with fully sex-linked genes, which requires closer linkage with the fully sex-linked region than that estimated for the PAR genes where apparent deviations from neutrality were observed. These results support the need to invoke selection to explain the S. latifolia PAR gene diversity, and encourage further work to test the possibility of balancing selection due to sexual antagonism.

Plant Sex Chromosomes

Although individuals in most flowering plant species, and in many haploid plants, have both sex functions, dioecious species-in which individuals have either male or female functions only-are scattered across many taxonomic groups, and many species have genetic sex determination. Among these, some have visibly heteromorphic sex chromosomes, and molecular genetic studies are starting to uncover sex-linked markers in others, showing that they too have fully sex-linked regions that are either too small or are located in chromosomes that are too small to be cytologically detectable from lack of pairing, lack of visible crossovers, or accumulation of heterochromatin. Detailed study is revealing that, like animal sex chromosomes, plant sex-linked regions show evidence for accumulation of repetitive sequences and genetic degeneration. Estimating when recombination stopped confirms the view that many plants have young sex-linked regions, making plants of great interest for studying the timescale of these changes.

Plant contributions to our understanding of sex chromosome evolution

A minority of angiosperms have male and female flowers separated in distinct individuals (dioecy), and most dioecious plants do not have cytologically different (heteromorphic) sex chromosomes. Plants nevertheless have several advantages for the study of sex chromosome evolution, as genetic sex determination has evolved repeatedly and is often absent in close relatives. I review sex-determining regions in non-model plant species, which may help us to understand when and how (and, potentially, test hypotheses about why) recombination suppression evolves within young sex chromosomes. I emphasize high-throughput sequencing approaches that are increasingly being applied to plants to test for non-recombining regions. These data are particularly illuminating when combined with sequence data that allow phylogenetic analyses, and estimates of when these regions evolved. Together with comparative genetic mapping, this has revealed that sex-determining loci and sex-linked regions evolved independently in many plant lineages, sometimes in closely related dioecious species, and often within the past few million years. In reviewing recent progress, I suggest areas for future work, such as the use of phylogenies to allow the informed choice of outgroup species suitable for inferring the directions of changes, including testing whether Y chromosome-like regions are undergoing genetic degeneration, a predicted consequence of losing recombination.

The status of supergenes in the 21st century: recombination suppression in Batesian mimicry and sex chromosomes and other complex adaptations

I review theoretical models for the evolution of supergenes in the cases of Batesian mimicry in butterflies, distylous plants and sex chromosomes. For each of these systems, I outline the genetic evidence that led to the proposal that they involve multiple genes that interact during ‘complex adaptations’, and at which the mutations involved are not unconditionally advantageous, but show advantages that trade‐off against some disadvantages. I describe recent molecular genetic studies of these systems and questions they raise about the evolution of suppressed recombination. Nonrecombining regions of sex chromosomes have long been known, but it is not yet fully understood why recombination suppression repeatedly evolved in systems in distantly related taxa, but does not always evolve. Recent studies of distylous plants are tending to support the existence of recombination‐suppressed genome regions, which may include modest numbers of genes and resemble recently evolved sex‐linked regions. For Batesian mimicry, however, molecular genetic work in two butterfly species suggests a new supergene scenario, with a single gene mutating to produce initial adaptive phenotypes, perhaps followed by modifiers specifically refining and perfecting the new phenotype.

RAD mapping reveals an evolving, polymorphic and fuzzy boundary of a plant pseudoautosomal region

How loss of genetic exchanges (recombination) evolves between sex chromosomes is a long-standing question. Suppressed recombination may evolve when a sexually antagonistic (SA) polymorphism occurs in a partially sex-linked 'pseudoautosomal' region (or 'PAR'), maintaining allele frequency differences between the two sexes, and creating selection for closer linkage with the fully sex-linked region of the Y chromosome in XY systems, or the W in ZW sex chromosome systems. Most evidence consistent with the SA polymorphism hypothesis is currently indirect, and more studies of the genetics and population genetics of PAR genes are clearly needed. The sex chromosomes of the plant Silene latifolia are suitable for such studies, as they evolved recently and the loss of recombination could still be ongoing. Here, we used RAD sequencing to genetically map sequences in this plant, which has a large genome (c. 3 gigabases) and no available whole-genome sequence. We mapped 83 genes on the sex chromosomes, and comparative mapping in the related species S. vulgaris supports previous evidence for additions to an ancestral PAR and identified at least 12 PAR genes. We describe evidence that recombination rates have been reduced in meiosis of both sexes, and differences in recombination between S. latifolia families suggest ongoing recombination suppression. Large allele frequency differences between the sexes were found at several loci closely linked to the PAR boundary, and genes in different regions of the PAR showed striking sequence diversity patterns that help illuminate the evolution of the PAR.

Independent stratum formation on the avian sex chromosomes reveals inter-chromosomal gene conversion and predominance of purifying selection on the W chromosome

We used a comparative approach spanning three species and 90 million years to study the evolutionary history of the avian sex chromosomes. Using whole transcriptomes, we assembled the largest cross-species dataset of W-linked coding content to date. Our results show that recombination suppression in large portions of the avian sex chromosomes has evolved independently, and that long-term sex chromosome divergence is consistent with repeated and independent inversions spreading progressively to restrict recombination. In contrast, over short-term periods we observe heterogeneous and locus-specific divergence. We also uncover four instances of gene conversion between both highly diverged and recently evolved gametologs, suggesting a complex mosaic of recombination suppression across the sex chromosomes. Lastly, evidence from 16 gametologs reveal that the W chromosome is evolving with a significant contribution of purifying selection, consistent with previous findings that W-linked genes play an important role in encoding sex-specific fitness.

Flooding tolerance in interspecific introgression lines containing chromosome segments from teosinte (Zea nicaraguensis) in maize (Zea mays subsp. mays)

BACKGROUND AND AIMS

Nicaraguan teosinte (Zea nicaraguensis), a species found in frequently flooded areas, provides useful germplasm for breeding flooding-tolerant maize (Z. mays subsp. mays). The objective of this study was to select flooding-tolerant lines using a library of introgression lines (ILs), each containing a chromosome segment from Z. nicaraguensis in the maize inbred line Mi29.

METHODS

To produce the ILs, a single F1 plant derived from a cross between maize Mi29 and Z. nicaraguensis was backcrossed to Mi29 three times, self-pollinated four times and genotyped using simple sequence repeat markers. Flooding tolerance was evaluated at the seedling stage under reducing soil conditions.

KEY RESULTS

By backcrossing and selfing, a series of 45 ILs were developed covering nearly the entire maize genome. Five flooding-tolerant lines were identified from among the ILs by evaluating leaf injury. Among these, line IL#18, containing a Z. nicaraguensis chromosome segment on the long arm of chromosome 4, showed the greatest tolerance to flooding, suggesting the presence of a major quantitative trait locus (QTL) in that region. The presence of the QTL was verified by examining flooding tolerance in a population segregating for the candidate region of chromosome 4. There was no significant relationship between the capacity to form constitutive aerenchyma and flooding tolerance in the ILs, indicating the presence of other factors related to flooding tolerance under reducing soil conditions.

CONCLUSIONS

A flooding-tolerant genotype, IL#18, was identified; this genotype should be useful for maize breeding. In addition, because the chromosome segments of Z. nicaraguensis in the ILs cover nearly the entire genome and Z. nicaraguensis possesses several unique traits related to flooding tolerance, the ILs should be valuable material for additional QTL detection and the development of flooding-tolerant maize lines.

Expansion of the pseudo-autosomal region and ongoing recombination suppression in the Silene latifolia sex chromosomes

There are two very interesting aspects to the evolution of sex chromosomes: what happens after recombination between these chromosome pairs stops and why suppressed recombination evolves. The former question has been intensively studied in a diversity of organisms, but the latter has been studied largely theoretically. To obtain empirical data, we used codominant genic markers in genetic mapping of the dioecious plant Silene latifolia, together with comparative mapping of S. latifolia sex-linked genes in S. vulgaris (a related hermaphrodite species without sex chromosomes). We mapped 29 S. latifolia fully sex-linked genes (including 21 newly discovered from transcriptome sequencing), plus 6 genes in a recombining pseudo-autosomal region (PAR) whose genetic map length is ∼25 cM in both male and female meiosis, suggesting that the PAR may contain many genes. Our comparative mapping shows that most fully sex-linked genes in S. latifolia are located on a single S. vulgaris linkage group and were probably inherited from a single autosome of an ancestor. However, unexpectedly, our maps suggest that the S. latifolia PAR region expanded through translocation events. Some genes in these regions still recombine in S. latifolia, but some genes from both addition events are now fully sex-linked. Recombination suppression is therefore still ongoing in S. latifolia, and multiple recombination suppression events have occurred in a timescale of few million years, much shorter than the timescale of formation of the most recent evolutionary strata of mammal and bird sex chromosomes.