Divergence and Functional Degradation of a Sex Chromosome-like Supergene

Students Don't Learn the Way They Think They Do in a Large, Active-Learning Genetics Course

Recently, our course team transformed a large-enrollment introductory genetics course from being predominantly lecture based to active learning based. During class sessions, students engaged in problem solving, which occurs when a student attempts to solve a problem without knowing the path to complete it. We designed class activities incorporating three distinct pedagogies from cognitive psychology: inquiry-based prediction, tell-then-practice case studies, and worked examples. We used a within-subjects design to compare students' attitudes toward these activities and their learning gains. Postsurvey results indicated that students felt worked examples helped them perform well on exams (a performance goal) and understand the information (a mastery goal) significantly better than the other activity types. However, students reported that all activity types required similar effort. Interestingly, students exhibited larger learning gains from prediction activities compared with worked examples or tell-then-practice activities, as evidenced by a course pretest/posttest. We discuss potential reasons for this misalignment between perceived helpfulness and actual learning gains. Additionally, we evaluate the strengths and weaknesses of each pedagogical approach.

Converging on long and short: The genetics, molecular biology and evolution of heterostyly

Heterostyly is a fascinating floral polymorphism that enhances outcrossing. In heterostylous species the flowers of the two or three morphs differ in multiple traits, including reciprocal reproductive-organ placement and self-incompatibility. These traits are controlled by individual genes within an S-locus supergene, whose suppressed recombination ensures the coordinated inheritance of the morph phenotypes. Recent breakthroughs about the genetic and molecular basis of heterostyly have resulted from studies on many independently evolved instances and include the following: The S-locus is a hemizygous region comprising several individual genes in multiple heterostylous taxa. In many systems, a single gene within the S-locus plays dual roles in regulating both female traits of style length and self-incompatibility type, often involving brassinosteroid signalling. The S-loci have evolved through stepwise or segmental duplication in different lineages. The frequent breakdown of heterostyly generally results from individual mutations at the S-locus and leads to a genomic selfing syndrome. These discoveries suggest convergent and genetically constrained evolution of heterostyly at the molecular level.

Large Inversions Shape Diversification and Genome Evolution in Common Quails

Chromosomal inversions, by suppressing recombination, can profoundly shape genome evolution and drive adaptation. In the common quail (Coturnix coturnix), a highly mobile bird with a vast Palearctic breeding range, we previously identified a massive inversion on chromosome 1 associated with distinct phenotypes and restricted geographic distribution. Here, using a new de novo genome assembly, we characterise this inversion and uncover additional, ancient structural variation on chromosome 2 that segregates across the species' range: either two putatively linked inversions or a single, large inversion that appears as two due to scaffolding limitations. Together, the inversions encompass a remarkable 15.6% of the quail genome (153.6 Mbp), creating highly divergent haplotypes that diverged over a million years ago. While the chromosome 1 inversion is linked to phenotypic differences, including morphology and migratory behaviour, the chromosome 2 inversion(s) show no such association. Notably, all inversion regions exhibit reduced effective population size and a relaxation of purifying selection, evidenced by elevated nonsynonymous-to-synonymous substitution ratios (N/S). This suggests that inversions, particularly the geographically restricted one on chromosome 1, may act as engines of diversification, accelerating the accumulation of functional variation and potentially contributing to local adaptation, especially within isolated island populations. Our findings demonstrate how large-scale chromosomal rearrangements can compartmentalise a genome, fostering distinct evolutionary trajectories within a single, highly mobile species.

Phylogenomics resolves key relationships in Rumex and uncovers a dynamic history of independently evolving sex chromosomes

Sex chromosomes have evolved independently many times across eukaryotes. Despite a considerable body of literature on sex chromosome evolution, the causes and consequences of variation in their formation, degeneration, and turnover remain poorly understood. Chromosomal rearrangements are thought to play an important role in these processes by promoting or extending the suppression of recombination on sex chromosomes. Sex chromosome variation may also contribute to barriers to gene flow, limiting introgression among species. Comparative approaches in groups with sexual system variation can be valuable for understanding these questions. Rumex is a diverse genus of flowering plants harboring significant sexual system and karyotypic variation, including hermaphroditic and dioecious clades with XY (and XYY) sex chromosomes. Previous disagreement in the phylogenetic relationships among key species has rendered the history of sex chromosome evolution uncertain. Resolving this history is important for investigating the interplay of chromosomal rearrangements, introgression, and sex chromosome evolution in the genus. Here, we use new transcriptome assemblies from 11 species representing major clades in the genus, along with a whole-genome assembly generated for a key hermaphroditic species. Using phylogenomic approaches, we find evidence for the independent evolution of sex chromosomes across two major clades, and introgression from unsampled lineages likely predating the formation of sex chromosomes in the genus. Comparative genomic approaches revealed high rates of chromosomal rearrangement, especially in dioecious species, with evidence for a complex origin of the sex chromosomes through multiple chromosomal fusions. However, we found no evidence of elevated rates of fusion on the sex chromosomes in comparison with autosomes, providing no support for an adaptive hypothesis of sex chromosome expansion due to sexually antagonistic selection. Overall, our results highlight a complex history of karyotypic evolution in Rumex, raising questions about the role that chromosomal rearrangements might play in the evolution of large heteromorphic sex chromosomes.

Large inversions in Lake Malawi cichlids are associated with habitat preference, lineage, and sex determination

Chromosomal inversions are an important class of genetic variation that link multiple alleles together into a single inherited block that can have important effects on fitness. To study the role of large inversions in the massive evolutionary radiation of Lake Malawi cichlids, we used long-read technologies to identify four single and two tandem inversions that span half of each respective chromosome, and which together encompass over 10% of the genome. Each inversion is fixed in one of the two states within the seven major ecogroups, suggesting they played a role in the separation of the major lake lineages into specific lake habitats. One exception is within the benthic sub-radiation, where both inverted and non-inverted alleles continue to segregate within the group. The evolutionary histories of three of the six inversions suggest they transferred from the pelagic Diplotaxodon group into benthic ancestors at the time the benthic sub-radiation was seeded. The remaining three inversions are found in a subset of benthic species living in deep waters. We show that some of these inversions are used as XY sex-determination systems but are also likely limited to a subset of total lake species. Our work suggests that inversions have been under both sexual and natural selection in Lake Malawi cichlids and that they will be important to understanding how this adaptive radiation evolved.

The interplay of local adaptation and gene flow may lead to the formation of supergenes

Supergenes are genetic architectures resulting in the segregation of alternative combinations of alleles underlying complex phenotypes. The co-segregation of alleles at linked loci is often facilitated by polymorphic chromosomal rearrangements suppressing recombination locally. Supergenes are involved in many complex polymorphisms, including sexual, colour or behavioural polymorphisms in numerous plants, fungi, mammals, fish, and insects. Despite a long history of empirical and theoretical research, the formation of supergenes remains poorly understood. Here, using a two-island population genetic model, we explore how gene flow and the evolution of overdominant chromosomal inversions may jointly lead to the formation of supergenes. We show that the evolution of inversions in differentiated populations, both under disruptive selection, leads to an increase in frequency of poorly adapted, immigrant haplotypes. Indeed, rare allelic combinations, such as immigrant haplotypes, are more frequently reshuffled by recombination than common allelic combinations, and therefore benefit from the recombination suppression generated by inversions. When an inversion capturing a locally adapted haplotype spreads but is associated with a fitness cost hampering its fixation (e.g. a recessive mutation load), the maintenance of a non-inverted haplotype in the population is enhanced; under certain conditions, the immigrant haplotype persists alongside the inverted local haplotype, while the standard local haplotype disappears. This establishes a stable, local polymorphism with two non-recombining haplotypes encoding alternative adaptive strategies, that is, a supergene. These results bring new light to the importance of local adaptation, overdominance, and gene flow in the formation of supergenes and inversion polymorphisms in general.

Assessing the Association Between Animal Color and Behavior: A Meta-Analysis of Experimental Studies

Color varies in pattern and degree across the tree of life. In animals, genetic variation in color is hypothesized to have pleiotropic effects on a variety of behaviors due to shared dependence on underlying biochemical pathways. Such pleiotropy can constrain the independent evolution of color and behavior. Although associations between color and behavior have been reported, this relationship has not yet been addressed across a broad taxonomic scale with a formal meta-analysis. We used a phylogenetic meta-analytic approach to examine the relationship between individual variation in aggressive behavior and variation in multiple colors. Seventy-four studies met our inclusion criteria (vertebrates = 70; invertebrates = 4). After accounting for phylogeny and correcting for publication bias, there was a positive association between measures of aggression and degree or area of coloration (mean = 0.248, 95% CI = (0.044, 0.477)). Because this positive association was not restricted to melanin-based coloration, we conclude that this pattern does not strongly support the melanin pleiotropy hypothesis. Because the association was also not affected by moderators accounting for individual condition, social rank, or age, the results do not strongly support hypotheses that condition dependence accounts for relationships between color and aggressive behavior. The badge of status hypothesis predicts that arbitrary traits can evolve to signal aggression or social dominance. We propose that this is the most parsimonious explanation for the patterns we observe. Because of the lack of evidence for condition dependence in the association between color and aggression, we further propose that the genetic covariation between traits contributes to the evolution of the badges of status.

Testing the evolutionary theory of inversion polymorphisms in the yellow monkeyflower (Mimulus guttatus)

Chromosomal inversions have been implicated in a remarkable range of natural phenomena, but it remains unclear how much they contribute to standing genetic variation. Here, we evaluate 64 inversions that segregate within a single natural population of the yellow monkeyflower (Mimulus guttatus). Nucleotide diversity patterns confirm low internal variation for the derived orientation (predicted by recent origin), elevated diversity between orientations (predicted by natural selection), and localized fluctuations (predicted by gene flux). Sequence divergence between orientations varies idiosyncratically by position, not following the suspension bridge pattern predicted if the breakpoints are the targets of selection. Genetic variation in gene expression is not inflated close to inversion breakpoints but is clearly partitioned between orientations. Like sequence variation, the pattern of expression variation suggests that the capture of coadapted alleles is more important than the breakpoints for the fitness effects of inversions. This work confirms several evolutionary predictions for inversion polymorphisms, but clarity emerges only by synthesizing estimates across many loci.

A supergene controls facultative diapause in the crop pest Helicoverpa armigera

Many insect species, including the economically important pest Helicoverpa armigera, avoid unfavorable conditions by suspending development. This form of phenotypic plasticity-facultative diapause-is a complex trait, though its evolution and intricate genetic architecture remain poorly understood. To investigate how such a polygenic trait could be locally adapted, we explore its genetic architecture. We map a large-effect diapause-associated locus to the Z chromosome by crossing high- and low-latitude populations. By generating multiple chromosome-scale assemblies, we identify an ∼5.93-Mb chromosomal inversion that constitutes the locus. Within this inversion, 33 genes harbor divergent non-synonymous mutations, notably including three circadian rhythm genes: Period, Clock, and Cycle. CRISPR-Cas9 knockout experiments confirm that each gene is independently essential for pupal diapause. Thus, a diapause supergene arose within H. armigera via a Z chromosome inversion, enabling local climatic adaptation in this economically important crop pest.

Fitness consequences of structural variation inferred from a House Finch pangenome

Genomic structural variants (SVs) play a crucial role in adaptive evolution, yet their average fitness effects and characterization with pangenome tools are understudied in wild animal populations. We constructed a pangenome for House Finches (Haemorhous mexicanus), a model for studies of host-pathogen coevolution, using long-read sequence data on 16 individuals (32 de novo-assembled haplotypes) and one outgroup. We identified 887,118 SVs larger than 50 base pairs, mostly (60%) involving repetitive elements, with reduced SV diversity in the eastern US as a result of its introduction by humans. The distribution of fitness effects of genome-wide SVs was estimated using maximum likelihood approaches and revealed that SVs in both coding and noncoding regions were on average more deleterious than smaller indels or single nucleotide polymorphisms. The reference-free pangenome facilitated identification of a > 10-My-old, 11-megabase-long pericentric inversion on chromosome 1. We found that the genotype frequencies of the inversion, estimated from 135 birds widely sampled temporally and geographically, increased steadily over the 25 y since House Finches were first exposed to the bacterial pathogen Mycoplasma gallisepticum and showed signatures of balancing selection, capturing genes related to immunity and telomerase activity. We also observed shorter telomeres in populations with a greater number of years exposure to Mycoplasma. Our study illustrates the utility of long-read sequencing and pangenome methods for understanding wild animal populations, estimating fitness effects of genome-wide SVs, and advancing our understanding of adaptive evolution through structural variation.

One-Step Leaping Evolution from an Autosomal Pair to the Heteromorphic Sex Chromosomes

BACKGROUND

Sex chromosomes evolve from an autosomal pair after the acquisition of a sex-determining gene. The primary sex chromosomes are homomorphic in both sexes and often undergo heteromorphism in either sex (XY in males or ZW in females) in association with chromosome rearrangements such as inversion, which creates a non-recombining region, called a stratum. Then, multiple strata may form by sequential inversions and extend the non-recombining region, where gene divergence accelerates, and degeneration of the Y or W chromosome progressively occurs.

SUMMARY

In contrast to the conventional theory, we propose a shortcut in heteromorphic sex chromosome evolution, where an autosomal pair directly evolves into a heteromorphic sex chromosome pair. We illustrate this with two frog cases where Y chromosome or autosome, which is morphologically inverted, was introgressed from another species through interspecific hybridization, instantly forming a new heteromorphic sex chromosome pair. This event resulted in a distinct non-recombining region immediately after hybridization.

KEY MESSAGES

The introduction of an inverted chromosome from a different species may be associated with benefits in morphology, breeding behavior, hybrid viability, sex determination, and recovery of the sex ratio of the hybrids. We discuss the molecular mechanisms driving preferential mutations in the introduced, inverted chromosome through interspecific hybridization.

Repeated evolution on oceanic islands: comparative genomics reveals species-specific processes in birds

Understanding the interplay between genetic drift, natural selection, gene flow, and demographic history in driving phenotypic and genomic differentiation of insular populations can help us gain insight into the speciation process. Comparing patterns across different insular taxa subjected to similar selective pressures upon colonizing oceanic islands provides the opportunity to study repeated evolution and identify shared patterns in their genomic landscapes of differentiation. We selected four species of passerine birds (Common Chaffinch Fringilla coelebs/canariensis, Red-billed Chough Pyrrhocorax pyrrhocorax, House Finch  Haemorhous mexicanus and Dark-eyed/island Junco Junco hyemalis/insularis) that have both mainland and insular populations. Changes in body size between island and mainland populations were consistent with the island rule. For each species, we sequenced whole genomes from mainland and insular individuals to infer their demographic history, characterize their genomic differentiation, and identify the factors shaping them. We estimated the relative (Fst) and absolute (dxy) differentiation, nucleotide diversity (π), Tajima's D, gene density and recombination rate. We also searched for selective sweeps and chromosomal inversions along the genome. All species shared a marked reduction in effective population size (Ne) upon island colonization. We found diverse patterns of differentiated genomic regions relative to the genome average in all four species, suggesting the role of selection in island-mainland differentiation, yet the lack of congruence in the location of these regions indicates that each species evolved differently in insular environments. Our results suggest that the genomic mechanisms involved in the divergence upon island colonization-such as chromosomal inversions, and historical factors like recurrent selection-differ in each species, despite the highly conserved structure of avian genomes and the similar selective factors involved. These differences are likely influenced by factors such as genetic drift, the polygenic nature of fitness traits and the action of case-specific selective pressures.

Morph-linked variation in female pheromone signalling and male response in a polymorphic moth

Understanding the maintenance of genetic variation in reproductive strategies and polymorphisms in the wild requires a comprehensive examination of the complex interactions between genetic basis, behaviour and environmental factors. We tested the association between three colour genotypes and variation in female pheromone signalling and male antennal morphology in the wood tiger moth (Arctia plantaginis). These moths have genetically determined white (WW, Wy) and yellow (yy) hindwings that are linked to mating success and fitness, with heterozygotes (Wy) having an advantage. We hypothesized that attractiveness and reproductive success are correlated, with Wy females being more attractive than the other two genotypes which could contribute to maintaining the polymorphism. Female attractiveness was tested by baiting traps with females of the three colour genotypes both in low- (i.e. field set-up) and in high-population density (i.e. large enclosure set-up). Male's ability to reach females was correlated to their own colour genotype and antennal morphology (length, area and lamellae count). Contrary to our prediction, morph-related reproductive success and attractiveness were not correlated. Heavier Wy females attracted a lower proportion of males compared to WW and yy females. Specifically, an increase in weight corresponded to a decreased Wy but increased yy female attractiveness. yy females were generally more attractive than others likely due to earlier pheromone release. In males, lamellae count and genetic colour morph were linked to the male's ability to locate females. Furthermore, male traits affected their ability to reach females in a context-specific way. Males with denser antennae (i.e. higher lamellae count) and white males reached the females faster than yellows in the enclosure, while yellow males located females faster than whites in the field. Our results indicate that higher yy female attractiveness was likely affected by the combined effect of early pheromone release, female weight and higher population density. Males' searching success was affected by morph-specific behavioural strategies and local population density. Ultimately, the combined effect of genotype-related pheromone signalling strategies of females together with environment-dependent male behaviour affect male response and potentially contribute to maintaining variation in fitness-related traits.

Sexually selected differences in warbler plumage are related to a putative inversion on the Z chromosome

Large structural variants in the genome, such as inversions, may play an important role in producing population structure and local adaptation to the environment through suppression of recombination. However, relatively few studies have linked inversions to phenotypic traits that are sexually selected and may play a role in reproductive isolation. Here, we found that geographic differences in the sexually selected plumage of a warbler, the common yellowthroat (Geothlypis trichas), are largely due to differences in the Z (sex) chromosome (males are ZZ), which contains at least one putative inversion spanning 40% (31/77 Mb) of its length. The inversions on the Z chromosome vary dramatically east and west of the Appalachian Mountains, which provides evidence of cryptic population structure within the range of the most widespread eastern subspecies (G. t. trichas). In an eastern (New York) and western (Wisconsin) population of this subspecies, female prefer different male ornaments; larger black facial masks are preferred in Wisconsin and larger yellow breasts are preferred in New York. The putative inversion also contains genes related to vision, which could influence mating preferences. Thus, structural variants on the Z chromosome are associated with geographic differences in male ornaments and female choice, which may provide a mechanism for maintaining different patterns of sexual selection in spite of gene flow between populations of the same subspecies.

Structural genomic variation and behavioral interactions underpin a balanced sexual mimicry polymorphism

How phenotypic diversity originates and persists within populations are classic puzzles in evolutionary biology. While balanced polymorphisms segregate within many species, it remains rare for both the genetic basis and the selective forces to be known, leading to an incomplete understanding of many classes of traits under balancing selection. Here, we uncover the genetic architecture of a balanced sexual mimicry polymorphism and identify behavioral mechanisms that may be involved in its maintenance in the swordtail fish Xiphophorus birchmanni. We find that ∼40% of X. birchmanni males develop a "false gravid spot," a melanic pigmentation pattern that mimics the "pregnancy spot" associated with sexual maturity in female live-bearing fish. Using genome-wide association mapping, we detect a single intergenic region associated with variation in the false gravid spot phenotype, which is upstream of kitlga, a melanophore patterning gene. By performing long-read sequencing within and across populations, we identify complex structural rearrangements between alternate alleles at this locus. The false gravid spot haplotype drives increased allele-specific expression of kitlga, which provides a mechanistic explanation for the increased melanophore abundance that causes the spot. By studying social interactions in the laboratory and in nature, we find that males with the false gravid spot experience less aggression; however, they also receive increased attention from other males and are disdained by females. These behavioral interactions may contribute to the maintenance of this phenotypic polymorphism in natural populations. We speculate that structural variants affecting gene regulation may be an underappreciated driver of balanced polymorphisms across diverse species.

A sex chromosome polymorphism maintains divergent plumage phenotypes between extensively hybridizing yellowhammers (Emberiza citrinella) and pine buntings (E. leucocephalos)

Under allopatric speciation, populations of a species become isolated by a geographic barrier and develop reproductive isolation through genetic differentiation. When populations meet in secondary contact, the strength of evolved reproductive barriers determines the extent of hybridization and whether the populations will continue to diverge or merge together. The yellowhammer (Emberiza citrinella) and pine bunting (E. leucocephalos) are avian sister species that diverged in allopatry on either side of Eurasia during the Pleistocene glaciations. Though they differ greatly in plumage and form distinct genetic clusters in allopatry, these taxa show negligible mitochondrial DNA differentiation and hybridize extensively where they overlap in central Siberia, lending uncertainty to the state of reproductive isolation in the system. To assess the strength of reproductive barriers between taxa, we examined genomic differentiation across the system. We found that extensive admixture has occurred in sympatry, indicating that reproductive barriers between taxa are weak. We also identified a putative Z chromosome inversion region that underlies plumage variation in the system, with the 'pine bunting' haplotype showing dominance over the 'yellowhammer' haplotype. Our results suggest that yellowhammers and pine buntings are currently at a crossroads and that evolutionary forces may push this system towards either continued differentiation or population merging. However, even if these taxa merge, recombination suppression between putative chromosome Z inversion haplotypes may maintain divergent plumage phenotypes within the system. In this way, our findings highlight the important role hybridization plays in increasing the genetic and phenotypic variation as well as the evolvability of a system.

Supergenes are not necessary to explain the maintenance of complex alternative phenotypes

Evolutionary biology aims to explain the diversity seen in nature. Evolutionary theory provides frameworks to understand how simple polymorphisms or continuous variation are maintained, but phenotypes inherited as discrete suites of quantitative traits are difficult to fit into this framework. Supergenes have been proposed as a solution to this problem-if causal genes are co-located, they can be inherited as if a single gene, thus bridging the gap between simple polymorphisms and continuous traits. We develop models to ask: how are critical supergenes for maintaining phenotypic diversity? In our simplest model, without explicit genetic architectures, three alternative reproductive morphs are maintained in many of the parameter combinations we evaluated. For these same parameter values, models with demographic stochasticity, recombination and mutation (but without explicit genetic architecture) maintained only two of these three morphs, with stochasticity determining which morphs persisted. With explicit genetic architectures, regardless of whether causal loci were co-located in a supergene or distributed randomly, this stochasticity in which morphs are maintained was reduced. Even when phenotypic variation was lost, genetic diversity was maintained. Altogether, categorical traits with polygenic bases exhibited similar evolutionary dynamics to those determined by supergenes. Our work suggests that supergenes are not the only answer to the puzzle of how discrete polygenic phenotypic variation is maintained.

The Primula edelbergii S-locus is an example of a jumping supergene

Research on supergenes, non-recombining genomic regions housing tightly linked genes that control complex phenotypes, has recently gained prominence in genomics. Heterostyly, a floral heteromorphism promoting outcrossing in several angiosperm families, is controlled by the S-locus supergene. The S-locus has been studied primarily in closely related Primula species and, more recently, in other groups that independently evolved heterostyly. However, it remains unknown whether genetic architecture and composition of the S-locus are maintained among species that share a common origin of heterostyly and subsequently diverged across larger time scales. To address this research gap, we present a chromosome-scale genome assembly of Primula edelbergii, a species that shares the same origin of heterostyly with Primula veris (whose S-locus has been characterized) but diverged from it 18 million years ago. Comparative genomic analyses between these two species allowed us to show, for the first time, that the S-locus can 'jump' (i.e. translocate) between chromosomes maintaining its function in controlling heterostyly. Additionally, we found that four S-locus genes were conserved but reshuffled within the supergene, seemingly without affecting their expression, thus we could not detect changes explaining the lack of self-incompatibility in P. edelbergii. Furthermore, we confirmed that the S-locus is not undergoing genetic degeneration. Finally, we investigated P. edelbergii evolutionary history within Ericales in terms of whole genome duplications and transposable element accumulation. In summary, our work provides a valuable resource for comparative analyses aimed at investigating the genetics of heterostyly and the pivotal role of supergenes in shaping the evolution of complex phenotypes.

Drivers of genomic differentiation landscapes in populations of disparate ecological and geographical settings within mainland Apis cerana

Elucidating the evolutionary processes that drive population divergence can enhance our understanding of the early stages of speciation and inform conservation management decisions. The honeybee Apis cerana displays extensive population divergence, providing an informative natural system for exploring these processes. The mainland lineage A. cerana includes several peripheral subspecies with disparate ecological and geographical settings radiated from a central ancestor. Under this evolutionary framework, we can explore the patterns of genome differentiation and the evolutionary models that explain them. We can also elucidate the contribution of non-genomic spatiotemporal mechanisms (extrinsic features) and genomic mechanisms (intrinsic features) that influence these genomic differentiation landscapes. Based on 293 whole genomes, a small part of the genome is highly differentiated between central-peripheral subspecies pairs, while low and partial parallelism partly reflects idiosyncratic responses to environmental differences. Combined elements of recurrent selection and speciation-with-gene-flow models generate the heterogeneous genome landscapes. These elements weight differently between central-island and other central-peripheral subspecies pairs, influenced by glacial cycles superimposed on different geomorphologies. Although local recombination rates exert a significant influence on patterns of genomic differentiation, it is unlikely that low-recombination rates regions were generated by structural variation. In conclusion, complex factors including geographical isolation, divergent ecological selection and non-uniform genome features have acted concertedly in the evolution of reproductive barriers that could reduce gene flow in part of the genome and facilitate the persistence of distinct populations within mainland lineage of A. cerana.

Evolution of Chromosomal Inversions across an Avian Radiation

Chromosomal inversions are structural mutations that can play a prominent role in adaptation and speciation. Inversions segregating across species boundaries (trans-species inversions) are often taken as evidence for ancient balancing selection or adaptive introgression, but can also be due to incomplete lineage sorting. Using whole-genome resequencing data from 18 populations of 11 recognized munia species in the genus Lonchura (N = 176 individuals), we identify four large para- and pericentric inversions ranging in size from 4 to 20 Mb. All four inversions cosegregate across multiple species and predate the numerous speciation events associated with the rapid radiation of this clade across the prehistoric Sahul (Australia, New Guinea) and Bismarck Archipelago. Using coalescent theory, we infer that trans-specificity is improbable for neutrally segregating variation despite substantial incomplete lineage sorting characterizing this young radiation. Instead, the maintenance of all three autosomal inversions (chr1, chr5, and chr6) is best explained by selection acting along ecogeographic clines not observed for the collinear parts of the genome. In addition, the sex chromosome inversion largely aligns with species boundaries and shows signatures of repeated positive selection for both alleles. This study provides evidence for trans-species inversion polymorphisms involved in both adaptation and speciation. It further highlights the importance of informing selection inference using a null model of neutral evolution derived from the collinear part of the genome.

How structural variants shape avian phenotypes: Lessons from model systems

Despite receiving significant recent attention, the relevance of structural variation (SV) in driving phenotypic diversity remains understudied, although recent advances in long-read sequencing, bioinformatics and pangenomic approaches have enhanced SV detection. We review the role of SVs in shaping phenotypes in avian model systems, and identify some general patterns in SV type, length and their associated traits. We found that most of the avian SVs so far identified are short indels in chickens, which are frequently associated with changes in body weight and plumage colouration. Overall, we found that relatively short SVs are more frequently detected, likely due to a combination of their prevalence compared to large SVs, and a detection bias, stemming primarily from the widespread use of short-read sequencing and associated analytical methods. SVs most commonly involve non-coding regions, especially introns, and when patterns of inheritance were reported, SVs associated primarily with dominant discrete traits. We summarise several examples of phenotypic convergence across different species, mediated by different SVs in the same or different genes and different types of changes in the same gene that can lead to various phenotypes. Complex rearrangements and supergenes, which can simultaneously affect and link several genes, tend to have pleiotropic phenotypic effects. Additionally, SVs commonly co-occur with single-nucleotide polymorphisms, highlighting the need to consider all types of genetic changes to understand the basis of phenotypic traits. We end by summarising expectations for when long-read technologies become commonly implemented in non-model birds, likely leading to an increase in SV discovery and characterisation. The growing interest in this subject suggests an increase in our understanding of the phenotypic effects of SVs in upcoming years.

Genomic landscapes of divergence among island bird populations: Evidence of parallel adaptation but at different loci?

When populations colonise new environments, they may be exposed to novel selection pressures but also suffer from extensive genetic drift due to founder effects, small population sizes and limited interpopulation gene flow. Genomic approaches enable us to study how these factors drive divergence, and disentangle neutral effects from differentiation at specific loci due to selection. Here, we investigate patterns of genetic diversity and divergence using whole-genome resequencing (>22× coverage) in Berthelot's pipit (Anthus berthelotii), a passerine endemic to the islands of three north Atlantic archipelagos. Strong environmental gradients, including in pathogen pressure, across populations in the species range, make it an excellent system in which to explore traits important in adaptation and/or incipient speciation. First, we quantify how genomic divergence accumulates across the speciation continuum, that is, among Berthelot's pipit populations, between sub species across archipelagos, and between Berthelot's pipit and its mainland ancestor, the tawny pipit (Anthus campestris). Across these colonisation timeframes (2.1 million-ca. 8000 years ago), we identify highly differentiated loci within genomic islands of divergence and conclude that the observed distributions align with expectations for non-neutral divergence. Characteristic signatures of selection are identified in loci associated with craniofacial/bone and eye development, metabolism and immune response between population comparisons. Interestingly, we find limited evidence for repeated divergence of the same loci across the colonisation range but do identify different loci putatively associated with the same biological traits in different populations, likely due to parallel adaptation. Incipient speciation across these island populations, in which founder effects and selective pressures are strong, may therefore be repeatedly associated with morphology, metabolism and immune defence.

A supergene-controlling social structure in Alpine ants also affects the dispersal ability and fecundity of each sex

Social organization, dispersal and fecundity coevolve, but whether they are genetically linked remains little known. Supergenes are prime candidates for coupling adaptive traits and mediating sex-specific trade-offs. Here, we test whether a supergene that controls social structure in Formica selysi also influences dispersal-related traits and fecundity within each sex. In this ant species, single-queen colonies contain only the ancestral supergene haplotype M and produce MM queens and M males, while multi-queen colonies contain the derived haplotype P and produce MP queens, PP queens and P males. By combining multiple experiments, we show that the M haplotype induces phenotypes with higher dispersal potential and higher fecundity in both sexes. Specifically, MM queens, MP queens and M males are more aerodynamic and more fecund than PP queens and P males, respectively. Differences between MP and PP queens from the same colonies reveal a direct genetic effect of the supergene on dispersal-related traits and fecundity. The derived haplotype P, associated with multi-queen colonies, produces queens and males with reduced dispersal abilities and lower fecundity. More broadly, similarities between the Formica and Solenopsis systems reveal that supergenes play a major role in linking behavioural, morphological and physiological traits associated with intraspecific social polymorphisms.

Genes and gene networks underlying spatial cognition in food-caching chickadees

Substantial progress has been made in understanding the genetic architecture of phenotypes involved in a variety of evolutionary processes. Behavioral genetics remains, however, among the least understood. We explore the genetic architecture of spatial cognitive abilities in a wild passerine bird, the mountain chickadee (Poecile gambeli). Mountain chickadees cache thousands of seeds in the fall and require specialized spatial memory to recover these caches throughout the winter. We previously showed that variation in spatial cognition has a direct effect on fitness and has a genetic basis. It remains unknown which specific genes and developmental pathways are particularly important for shaping spatial cognition. To further dissect the genetic basis of spatial cognitive abilities, we combine experimental quantification of spatial cognition in wild chickadees with whole-genome sequencing of 162 individuals, a new chromosome-scale reference genome, and species-specific gene annotation. We have identified a set of genes and developmental pathways that play a key role in creating variation in spatial cognition and found that the mechanism shaping cognitive variation is consistent with selection against mildly deleterious non-coding mutations. Although some candidate genes were organized into connected gene networks, about half do not have shared regulation, highlighting that multiple independent developmental or physiological mechanisms contribute to variation in spatial cognitive abilities. A large proportion of the candidate genes we found are associated with synaptic plasticity, an intriguing result that leads to the hypothesis that certain genetic variants create antagonism between behavioral plasticity and long-term memory, each providing distinct benefits depending on ecological context.

Evolution and genetic architecture of sex-limited polymorphism in cuckoos

Sex-limited polymorphism has evolved in many species including our own. Yet, we lack a detailed understanding of the underlying genetic variation and evolutionary processes at work. The brood parasitic common cuckoo (Cuculus canorus) is a prime example of female-limited color polymorphism, where adult males are monochromatic gray and females exhibit either gray or rufous plumage. This polymorphism has been hypothesized to be governed by negative frequency-dependent selection whereby the rarer female morph is protected against harassment by males or from mobbing by parasitized host species. Here, we show that female plumage dichromatism maps to the female-restricted genome. We further demonstrate that, consistent with balancing selection, ancestry of the rufous phenotype is shared with the likewise female dichromatic sister species, the oriental cuckoo (Cuculus optatus). This study shows that sex-specific polymorphism in trait variation can be resolved by genetic variation residing on a sex-limited chromosome and be maintained across species boundaries.

A region of suppressed recombination misleads neoavian phylogenomics

Genomes are typically mosaics of regions with different evolutionary histories. When speciation events are closely spaced in time, recombination makes the regions sharing the same history small, and the evolutionary history changes rapidly as we move along the genome. When examining rapid radiations such as the early diversification of Neoaves 66 Mya, typically no consistent history is observed across segments exceeding kilobases of the genome. Here, we report an exception. We found that a 21-Mb region in avian genomes, mapped to chicken chromosome 4, shows an extremely strong and discordance-free signal for a history different from that of the inferred species tree. Such a strong discordance-free signal, indicative of suppressed recombination across many millions of base pairs, is not observed elsewhere in the genome for any deep avian relationships. Although long regions with suppressed recombination have been documented in recently diverged species, our results pertain to relationships dating circa 65 Mya. We provide evidence that this strong signal may be due to an ancient rearrangement that blocked recombination and remained polymorphic for several million years prior to fixation. We show that the presence of this region has misled previous phylogenomic efforts with lower taxon sampling, showing the interplay between taxon and locus sampling. We predict that similar ancient rearrangements may confound phylogenetic analyses in other clades, pointing to a need for new analytical models that incorporate the possibility of such events.

Remarkably High Repeat Content in the Genomes of Sparrows: The Importance of Genome Assembly Completeness for Transposable Element Discovery

Transposable elements (TE) play critical roles in shaping genome evolution. Highly repetitive TE sequences are also a major source of assembly gaps making it difficult to fully understand the impact of these elements on host genomes. The increased capacity of long-read sequencing technologies to span highly repetitive regions promises to provide new insights into patterns of TE activity across diverse taxa. Here we report the generation of highly contiguous reference genomes using PacBio long-read and Omni-C technologies for three species of Passerellidae sparrow. We compared these assemblies to three chromosome-level sparrow assemblies and nine other sparrow assemblies generated using a variety of short- and long-read technologies. All long-read based assemblies were longer (range: 1.12 to 1.41 Gb) than short-read assemblies (0.91 to 1.08 Gb) and assembly length was strongly correlated with the amount of repeat content. Repeat content for Bell's sparrow (31.2% of genome) was the highest level ever reported within the order Passeriformes, which comprises over half of avian diversity. The highest levels of repeat content (79.2% to 93.7%) were found on the W chromosome relative to other regions of the genome. Finally, we show that proliferation of different TE classes varied even among species with similar levels of repeat content. These patterns support a dynamic model of TE expansion and contraction even in a clade where TEs were once thought to be fairly depauperate and static. Our work highlights how the resolution of difficult-to-assemble regions of the genome with new sequencing technologies promises to transform our understanding of avian genome evolution.

Multiple Transitions between Y Chromosome and Autosome in Tago's Brown Frog Species Complex

Sex chromosome turnover is the transition between sex chromosomes and autosomes. Although many cases have been reported in poikilothermic vertebrates, their evolutionary causes and genetic mechanisms remain unclear. In this study, we report multiple transitions between the Y chromosome and autosome in the Japanese Tago's brown frog complex. Using chromosome banding and molecular analyses (sex-linked and autosomal single nucleotide polymorphisms, SNPs, from the nuclear genome), we investigated the frogs of geographic populations ranging from northern to southern Japan of two species, Rana tagoi and Rana sakuraii (2n = 26). Particularly, the Chiba populations of East Japan and Akita populations of North Japan in R. tagoi have been, for the first time, investigated here. As a result, we identified three different sex chromosomes, namely chromosomes 3, 7, and 13, in the populations of the two species. Furthermore, we found that the transition between the Y chromosome (chromosome 7) and autosome was repeated through hybridization between two or three different populations belonging to the two species, followed by restricted chromosome introgression. These dynamic sex chromosome turnovers represent the first such findings in vertebrates and imply that speciation associated with inter- or intraspecific hybridization plays an important role in sex chromosome turnover in frogs.

Spatial variation in population genomic responses to over a century of anthropogenic change within a tidal marsh songbird

Combating the current biodiversity crisis requires the accurate documentation of population responses to human-induced ecological change. However, our ability to pinpoint population responses to human activities is often limited to the analysis of populations studied well after the fact. Museum collections preserve a record of population responses to anthropogenic change that can provide critical baseline data on patterns of genetic diversity, connectivity, and population structure prior to the onset of human perturbation. Here, we leverage a spatially replicated time series of specimens to document population genomic responses to the destruction of nearly 90% of coastal habitats occupied by the Savannah sparrow (Passerculus sandwichensis) in California. We sequenced 219 sparrows collected from 1889 to 2017 across the state of California using an exome capture approach. Spatial-temporal analyses of genetic diversity found that the amount of habitat lost was not predictive of genetic diversity loss. Sparrow populations from southern California historically exhibited lower levels of genetic diversity and experienced the most significant temporal declines in genetic diversity. Despite experiencing the greatest levels of habitat loss, we found that genetic diversity in the San Francisco Bay area remained relatively high. This was potentially related to an observed increase in gene flow into the Bay Area from other populations. While gene flow may have minimized genetic diversity declines, we also found that immigration from inland freshwater-adapted populations into tidal marsh populations led to the erosion of divergence at loci associated with tidal marsh adaptation. Shifting patterns of gene flow through time in response to habitat loss may thus contribute to negative fitness consequences and outbreeding depression. Together, our results underscore the importance of tracing the genomic trajectories of multiple populations over time to address issues of fundamental conservation concern.

The genomics and evolution of inter-sexual mimicry and female-limited polymorphisms in damselflies

Sex-limited morphs can provide profound insights into the evolution and genomic architecture of complex phenotypes. Inter-sexual mimicry is one particular type of sex-limited polymorphism in which a novel morph resembles the opposite sex. While inter-sexual mimics are known in both sexes and a diverse range of animals, their evolutionary origin is poorly understood. Here, we investigated the genomic basis of female-limited morphs and male mimicry in the common bluetail damselfly. Differential gene expression between morphs has been documented in damselflies, but no causal locus has been previously identified. We found that male mimicry originated in an ancestrally sexually dimorphic lineage in association with multiple structural changes, probably driven by transposable element activity. These changes resulted in ~900 kb of novel genomic content that is partly shared by male mimics in a close relative, indicating that male mimicry is a trans-species polymorphism. More recently, a third morph originated following the translocation of part of the male-mimicry sequence into a genomic position ~3.5 mb apart. We provide evidence of balancing selection maintaining male mimicry, in line with previous field population studies. Our results underscore how structural variants affecting a handful of potentially regulatory genes and morph-specific genes can give rise to novel and complex phenotypic polymorphisms.

Social antagonism facilitates supergene expansion in ants

Antagonistic selection has long been considered a major driver of the formation and expansion of sex chromosomes. For example, sexually antagonistic variation on an autosome can select for suppressed recombination between that autosome and the sex chromosome, leading to a neo-sex chromosome. Autosomal supergenes, chromosomal regions containing tightly linked variants affecting the same complex trait, share similarities with sex chromosomes, raising the possibility that sex chromosome evolution models can explain the evolution of genome structure and recombination in other contexts. We tested this premise in a Formica ant species, wherein we identified four supergene haplotypes on chromosome 3 underlying colony social organization and sex ratio. We discovered a novel rearranged supergene variant (9r) on chromosome 9 underlying queen miniaturization. The 9r is in strong linkage disequilibrium with one chromosome 3 haplotype (P2) found in multi-queen (polygyne) colonies. We suggest that queen miniaturization is strongly disfavored in the single-queen (monogyne) background and is thus socially antagonistic. As such, divergent selection experienced by ants living in alternative social "environments" (monogyne and polygyne) may have contributed to the emergence of a genetic polymorphism on chromosome 9 and associated queen-size dimorphism. Consequently, an ancestral polygyne-associated haplotype may have expanded to include the polymorphism on chromosome 9, resulting in a larger region of suppressed recombination spanning two chromosomes. This process is analogous to the formation of neo-sex chromosomes and consistent with models of expanding regions of suppressed recombination. We propose that miniaturized queens, 16%-20% smaller than queens without 9r, could be incipient intraspecific social parasites.

How chromosomal inversions reorient the evolutionary process

Inversions are structural mutations that reverse the sequence of a chromosome segment and reduce the effective rate of recombination in the heterozygous state. They play a major role in adaptation, as well as in other evolutionary processes such as speciation. Although inversions have been studied since the 1920s, they remain difficult to investigate because the reduced recombination conferred by them strengthens the effects of drift and hitchhiking, which in turn can obscure signatures of selection. Nonetheless, numerous inversions have been found to be under selection. Given recent advances in population genetic theory and empirical study, here we review how different mechanisms of selection affect the evolution of inversions. A key difference between inversions and other mutations, such as single nucleotide variants, is that the fitness of an inversion may be affected by a larger number of frequently interacting processes. This considerably complicates the analysis of the causes underlying the evolution of inversions. We discuss the extent to which these mechanisms can be disentangled, and by which approach.

Overcoming the congenitally disadvantageous mutation through adaptation to environmental UV exposure in land snails

Congenital fitness-disadvantageous mutations are not maintained in the population; they are purged from the population through processes such as purifying selection. However, these mutations could persist in the population as polymorphisms when it is advantageous for the individuals carrying them in adapting to a specific external environment. We tested this hypothesis using the dimorphic land snail Euhadra peliomphala simodae in Japan; these snails have dark or bright coloured shells. The survival rate of dark snails at hatching was lower than that of the bright ones, as observed in the F1 progenies produced through crossing. Dark snails have a congenital fitness-disadvantageous mutation; however, they also have protection against ultraviolet radiation. They have a higher survival rate than the bright snails in a UV environment, as observed using the UV exposure experiments and UV transmittance measurements. This is a good example of a congenitally disadvantageous mutation that is advantageous for adapting to the external environment. These results explain the maintenance of polymorphism and highlight the genotypic and phenotypic diversity in the wild population.

A supergene in seaweed flies modulates male traits and female perception

Supergenes, tightly linked sets of alleles, offer some of the most spectacular examples of polymorphism persisting under long-term balancing selection. However, we still do not understand their evolution and persistence, especially in the face of accumulation of deleterious elements. Here, we show that an overdominant supergene in seaweed flies, Coelopa frigida, modulates male traits, potentially facilitating disassortative mating and promoting intraspecific polymorphism. Across two continents, the Cf-Inv(1) supergene strongly affected the composition of male cuticular hydrocarbons (CHCs) but only weakly affected CHC composition in females. Using gas chromatography-electroantennographic detection, we show that females can sense male CHCs and that there may be differential perception between genotypes. Combining our phenotypic results with RNA-seq data, we show that candidate genes for CHC biosynthesis primarily show differential expression for Cf-Inv(1) in males but not females. Conversely, candidate genes for odorant detection were differentially expressed in both sexes but showed high levels of divergence between supergene haplotypes. We suggest that the reduced recombination between supergene haplotypes may have led to rapid divergence in mate preferences as well as increasing linkage between male traits, and overdominant loci. Together this probably helped to maintain the polymorphism despite deleterious effects in homozygotes.

Multivariate Models of Animal Sex: Breaking Binaries Leads to a Better Understanding of Ecology and Evolution

"Sex" is often used to describe a suite of phenotypic and genotypic traits of an organism related to reproduction. However, these traits-gamete type, chromosomal inheritance, physiology, morphology, behavior, etc.-are not necessarily coupled, and the rhetorical collapse of variation into a single term elides much of the complexity inherent in sexual phenotypes. We argue that consideration of "sex" as a constructed category operating at multiple biological levels opens up new avenues for inquiry in our study of biological variation. We apply this framework to three case studies that illustrate the diversity of sex variation, from decoupling sexual phenotypes to the evolutionary and ecological consequences of intrasexual polymorphisms. We argue that instead of assuming binary sex in these systems, some may be better categorized as multivariate and nonbinary. Finally, we conduct a meta-analysis of terms used to describe diversity in sexual phenotypes in the scientific literature to highlight how a multivariate model of sex can clarify, rather than cloud, studies of sexual diversity within and across species. We argue that such an expanded framework of "sex" better equips us to understand evolutionary processes, and that as biologists, it is incumbent upon us to push back against misunderstandings of the biology of sexual phenotypes that enact harm on marginalized communities.

Whole genome assessment of a declining game bird reveals cryptic genetic structure and insights for population management

Population genomics applied to game species conservation can help delineate management units, ensure appropriate harvest levels and identify populations needing genetic rescue to safeguard their adaptive potential. The ruffed grouse (Bonasa umbellus) is rapidly declining in much of the eastern USA due to a combination of forest maturation and habitat fragmentation. More recently, mortality from West Nile Virus may have affected connectivity of local populations; however, genetic approaches have never explicitly investigated this issue. In this study, we sequenced 54 individual low-coverage (~5X) grouse genomes to characterize population structure, assess migration rates across the landscape to detect potential barriers to gene flow and identify genomic regions with high differentiation. We identified two genomic clusters with no clear geographic correlation, with large blocks of genomic differentiation associated with chromosomes 4 and 20, likely due to chromosomal inversions. After excluding these putative inversions from the data set, we found weak but nonsignificant signals of population subdivision. Estimated gene flow revealed reduced rates of migration in areas with extensive habitat fragmentation and increased genetic connectivity in areas with less habitat fragmentation. Our findings provide a benchmark for wildlife managers to compare and scale the genetic diversity and structure of ruffed grouse populations in Pennsylvania and across the eastern USA, and we also reveal structural variation in the grouse genome that requires further study to understand its possible effects on individual fitness and population distribution.

Born without a Silver Spoon: A Review of the Causes and Consequences of Adversity during Early Life

A huge amount of research attention has focused on the evolution of life histories, but most research focuses on dominant individuals that acquire a disproportionate level of reproductive success, while the life histories and reproductive tactics of subordinate individuals have received less attention. Here, we review the links between early life adversity and performance during adulthood in birds, and highlight instances in which subordinate individuals outperform dominant conspecifics. Subordinate individuals are those from broods raised under high risk of predation, with low availability of food, and/or with many parasites. Meanwhile, the broods of many species hatch or are born asynchronously and mitigation of the asynchrony is generally lacking from variation in maternal effects such as egg size and hormone deposition or genetic effects such as offspring sex or parentage. Subordinate individuals employ patterns of differential growth to attempt to mitigate the adversity they experience during early life, yet they overwhelmingly fail to overcome their initial handicap. In terms of surviving through to adulthood, subordinate individuals employ other "suboptimal" tactics, such as adaptively timing foraging behaviors to avoid dominant individuals. During adulthood, meanwhile, subordinate individuals rely on "suboptimal" tactics, such as adaptive dispersal behaviors and competing for partners at optimal times, because they represent the best options available to them to acquire copulations whenever possible. We conclude that there is a gap in knowledge for direct links between early life adversity and subordination during adulthood, meaning that further research should test for links. There are instances, however, where subordinate individuals employ "suboptimal" tactics that allow them to outperform dominant conspecifics during adulthood.

Evidence for a single, ancient origin of a genus-wide alternative life history strategy

Understanding the evolutionary origins and factors maintaining alternative life history strategies (ALHS) within species is a major goal of evolutionary research. While alternative alleles causing discrete ALHS are expected to purge or fix over time, one-third of the ~90 species of Colias butterflies are polymorphic for a female-limited ALHS called Alba. Whether Alba arose once, evolved in parallel, or has been exchanged among taxa is currently unknown. Using comparative genome-wide association study (GWAS) and population genomic analyses, we placed the genetic basis of Alba in time-calibrated phylogenomic framework, revealing that Alba evolved once near the base of the genus and has been subsequently maintained via introgression and balancing selection. CRISPR-Cas9 mutagenesis was then used to verify a putative cis-regulatory region of Alba, which we identified using phylogenetic foot printing. We hypothesize that this cis-regulatory region acts as a modular enhancer for the induction of the Alba ALHS, which has likely facilitated its long evolutionary persistence.

Comparative phylogeography reveals widespread cryptic diversity driven by ecology in Panamanian birds

Widespread species often harbor unrecognized genetic diversity, and investigating the factors associated with such cryptic variation can help us better understand the forces driving diversification. Here, we identify potential cryptic species based on a comprehensive dataset of COI mitochondrial DNA barcodes from 2,333 individual Panamanian birds across 429 species, representing 391 (59%) of the 659 resident landbird species of the country, as well as opportunistically sampled waterbirds. We complement this dataset with additional publicly available mitochondrial loci, such as ND2 and cytochrome b, obtained from whole mitochondrial genomes from 20 taxa. Using barcode identification numbers (BINs), we find putative cryptic species in 19% of landbird species, highlighting hidden diversity in the relatively well-described avifauna of Panama. Whereas some of these mitochondrial divergence events corresponded with recognized geographic features that likely isolated populations, such as the Cordillera Central highlands, the majority (74%) of lowland splits were between eastern and western populations. The timing of these splits are not temporally coincident across taxa, suggesting that historical events, such as the formation of the Isthmus of Panama and Pleistocene climatic cycles, were not the primary drivers of cryptic diversification. Rather, we observed that forest species, understory species, insectivores, and strongly territorial species-all traits associated with lower dispersal ability-were all more likely to have multiple BINs in Panama, suggesting strong ecological associations with cryptic divergence. Additionally, hand-wing index, a proxy for dispersal capability, was significantly lower in species with multiple BINs, indicating that dispersal ability plays an important role in generating diversity in Neotropical birds. Together, these results underscore the need for evolutionary studies of tropical bird communities to consider ecological factors along with geographic explanations, and that even in areas with well-known avifauna, avian diversity may be substantially underestimated.

LAY SUMMARY

- What factors are common among bird species with cryptic diversity in Panama? What role do geography, ecology, phylogeographic history, and other factors play in generating bird diversity?- 19% of widely-sampled bird species form two or more distinct DNA barcode clades, suggesting widespread unrecognized diversity.- Traits associated with reduced dispersal ability, such as use of forest understory, high territoriality, low hand-wing index, and insectivory, were more common in taxa with cryptic diversity. Filogeografía comparada revela amplia diversidad críptica causada por la ecología en las aves de Panamá.

RESUMEN

Especies extendidas frecuentemente tiene diversidad genética no reconocida, y investigando los factores asociados con esta variación críptica puede ayudarnos a entender las fuerzas que impulsan la diversificación. Aquí, identificamos especies crípticas potenciales basadas en un conjunto de datos de códigos de barras de ADN mitocondrial de 2,333 individuos de aves de Panama en 429 especies, representando 391 (59%) de las 659 especies de aves terrestres residentes del país, además de algunas aves acuáticas muestreada de manera oportunista. Adicionalmente, complementamos estos datos con secuencias mitocondriales disponibles públicamente de otros loci, tal como ND2 o citocroma b, obtenidos de los genomas mitocondriales completos de 20 taxones. Utilizando los números de identificación de código de barras (en ingles: BINs), un sistema taxonómico numérico que proporcina una estimación imparcial de la diversidad potencial a nivel de especie, encontramos especies crípticas putativas en 19% de las especies de aves terrestres, lo que destaca la diversidad oculta en la avifauna bien descrita de Panamá. Aunque algunos de estos eventos de divergencia conciden con características geográficas que probablemente aislaron las poblaciones, la mayoría (74%) de la divergencia en las tierras bajas se encuentra entre las poblaciones orientales y occidentales. El tiempo de esta divergencia no coincidió entre los taxones, sugiriendo que eventos históricos tales como la formación del Istmo de Panamá y los ciclos climáticos del pleistoceno, no fueron los principales impulsores de la especiación. En cambio, observamos asociaciones fuertes entre las características ecológicas y la divergencia mitocondriale: las especies del bosque, sotobosque, con una dieta insectívora, y con territorialidad fuerte mostraton múltiple BINs probables. Adicionalmente, el índice mano-ala, que está asociado a la capacidad de dispersión, fue significativamente menor en las especies con BINs multiples, sugiriendo que la capacidad de dispersión tiene un rol importamente en la generación de la diversidad de las aves neotropicales. Estos resultos demonstran la necesidad de que estudios evolutivos de las comunidades de aves tropicales consideren los factores ecológicos en conjunto con las explicaciones geográficos. Palabras clave: biodiversidad tropical, biogeografía, códigos de barras, dispersión, especies crípticas.

Inversions maintain differences between migratory phenotypes of a songbird

Structural rearrangements have been shown to be important in local adaptation and speciation, but have been difficult to reliably identify and characterize in non-model species. Here we combine long reads, linked reads and optical mapping to characterize three divergent chromosome regions in the willow warbler Phylloscopus trochilus, of which two are associated with differences in migration and one with an environmental gradient. We show that there are inversions (0.4-13 Mb) in each of the regions and that the divergence times between inverted and non-inverted haplotypes are similar across the regions (~1.2 Myrs), which is compatible with a scenario where inversions arose in either of two allopatric populations that subsequently hybridized. The improved genomes allow us to detect additional functional differences in the divergent regions, providing candidate genes for migration and adaptations to environmental gradients.

Chromosomal inversion polymorphisms shape the genomic landscape of deer mice

Chromosomal inversions are an important form of structural variation that can affect recombination, chromosome structure and fitness. However, because inversions can be challenging to detect, the prevalence and hence the significance of inversions segregating within species remains largely unknown, especially in natural populations of mammals. Here, by combining population-genomic and long-read sequencing analyses in a single, widespread species of deer mouse (Peromyscus maniculatus), we identified 21 polymorphic inversions that are large (1.5-43.8 Mb) and cause near-complete suppression of recombination when heterozygous (0-0.03 cM Mb-1). We found that inversion breakpoints frequently occur in centromeric and telomeric regions and are often flanked by long inverted repeats (0.5-50 kb), suggesting that they probably arose via ectopic recombination. By genotyping inversions in populations across the species' range, we found that the inversions are often widespread and do not harbour deleterious mutational loads, and many are likely to be maintained as polymorphisms by divergent selection. Comparisons of forest and prairie ecotypes of deer mice revealed 13 inversions that contribute to differentiation between populations, of which five exhibit significant associations with traits implicated in local adaptation. Taken together, these results show that inversion polymorphisms have a significant impact on recombination, genome structure and genetic diversity in deer mice and likely facilitate local adaptation across the widespread range of this species.

Origin and persistence of polymorphism in loci targeted by disassortative preference: a general model

The emergence and persistence of polymorphism within populations generally requires specific regimes of natural or sexual selection. Here, we develop a unified theoretical framework to explore how polymorphism at targeted loci can be generated and maintained by either disassortative mating choice or balancing selection due to, for example, heterozygote advantage. To this aim, we model the dynamics of alleles at a single locus A in a population of haploid individuals, where reproductive success depends on the combination of alleles carried by the parents at locus A. Our theoretical study of the model confirms that the conditions for the persistence of a given level of allelic polymorphism depend on the relative reproductive advantages among pairs of individuals. Interestingly, equilibria with unbalanced allelic frequencies were shown to emerge from successive introduction of mutants. We then investigate the role of the function linking allelic divergence to reproductive advantage on the evolutionary fate of alleles within the population. Our results highlight the significance of the shape of this function for both the number of alleles maintained and their level of genetic divergence. Large number of alleles are maintained with substantial replacement of alleles, when disassortative advantage slowly increases with allelic differentiation . In contrast, few highly differentiated alleles are predicted to be maintained when genetic differentiation has a strong effect on disassortative advantage. These opposite effects predicted by our model explain how disassortative mate choice may lead to various levels of allelic differentiation and polymorphism, and shed light on the effect of mate preferences on the persistence of balanced and unbalanced polymorphism in natural population.

Chromosome-scale assemblies of the male and female Populus euphratica genomes reveal the molecular basis of sex determination and sexual dimorphism

Reference-quality genomes of both sexes are essential for studying sex determination and sex-chromosome evolution, as their gene contents and expression profiles differ. Here, we present independent chromosome-level genome assemblies for the female (XX) and male (XY) genomes of desert poplar (Populus euphratica), resolving a 22.7-Mb X and 24.8-Mb Y chromosome. We also identified a relatively complete 761-kb sex-linked region (SLR) in the peritelomeric region on chromosome 14 (Y). Within the SLR, recombination around the partial repeats for the feminizing factor ARR17 (ARABIDOPSIS RESPONSE REGULATOR 17) was potentially suppressed by flanking palindromic arms and the dense accumulation of retrotransposons. The inverted small segments S1 and S2 of ARR17 exhibited relaxed selective pressure and triggered sex determination by generating 24-nt small interfering RNAs that induce male-specific hyper-methylation at the promoter of the autosomal targeted ARR17. We also detected two male-specific fusion genes encoding proteins with NB-ARC domains at the breakpoint region of an inversion in the SLR that may be responsible for the observed sexual dimorphism in immune responses. Our results show that the SLR appears to follow proposed evolutionary dynamics for sex chromosomes and advance our understanding of sex determination and the evolution of sex chromosomes in Populus.

Reference-quality genomes of both sexes of the dioecious tree species, Populus euphratica, provide further insight into the evolution of Populus sex chromosomes and highlight male-specific fusion genes that may contribute to sexual dimorphism.

Allele-specific cis-regulatory methylation of the gene for vasoactive intestinal peptide in white-throated sparrows

White-throated sparrows (Zonotrichia albicollis) offer a unique opportunity to connect genotype with behavioral phenotype. In this species, a rearrangement of the second chromosome is linked with territorial aggression; birds with a copy of this "supergene" rearrangement are more aggressive than those without it. The supergene has captured the gene VIP, which encodes vasoactive intestinal peptide, a neuromodulator that drives aggression in other songbirds. In white-throated sparrows, VIP expression is higher in the anterior hypothalamus of birds with the supergene than those without it, and expression of VIP in this region predicts the level of territorial aggression regardless of genotype. Here, we aimed to identify epigenetic mechanisms that could contribute to differential expression of VIP both in breeding adults, which exhibit morph differences in territorial aggression, and in nestlings, before territorial behavior develops. We extracted and bisulfite-converted DNA from samples of the hypothalamus in wild-caught adults and nestlings and used high-throughput sequencing to measure DNA methylation of a region upstream of the VIP start site. We found that the allele inside the supergene was less methylated than the alternative allele in both adults and nestlings. The differential methylation was attributed primarily to CpG sites that were shared between the alleles, not to polymorphic sites, which suggests that epigenetic regulation is occurring independently of the genetic differentiation within the supergene. This work represents an initial step toward understanding how epigenetic differentiation inside chromosomal inversions leads to the development of alternative behavioral phenotypes.

Temperature-associated selection linked to putative chromosomal inversions in king scallop (Pecten maximus)

The genomic landscape of divergence—the distribution of differences among populations or species across the genome—is increasingly characterized to understand the role that microevolutionary forces such as natural selection and recombination play in causing and maintaining genetic divergence. This line of inquiry has also revealed chromosome structure variation to be an important factor shaping the landscape of adaptive genetic variation. Owing to a high prevalence of chromosome structure variation and the strong pressure for local adaptation necessitated by their sessile nature, bivalve molluscs are an ideal taxon for exploring the relationship between chromosome structure variation and local adaptation. Here, we report a population genomic survey of king scallop (Pecten maximus) across its natural range in the northeastern Atlantic Ocean, using a recent chromosome-level genome assembly. We report the presence of at least three large (12–22 Mb), putative chromosomal inversions associated with sea surface temperature and whose frequencies are in contrast to neutral population structure. These results highlight a potentially large role for recombination-suppressing chromosomal inversions in local adaptation and suggest a hypothesis to explain the maintenance of differences in reproductive timing found at relatively small spatial scales across king scallop populations.

Dynamic molecular evolution of a supergene with suppressed recombination in white-throated sparrows

In white-throated sparrows, two alternative morphs differing in plumage and behavior segregate with a large chromosomal rearrangement. As with sex chromosomes such as the mammalian Y, the rearranged version of chromosome two (ZAL2m) is in a near-constant state of heterozygosity, offering opportunities to investigate both degenerative and selective processes during the early evolutionary stages of 'supergenes.' Here, we generated, synthesized, and analyzed extensive genome-scale data to better understand the forces shaping the evolution of the ZAL2 and ZAL2m chromosomes in this species. We found that features of ZAL2m are consistent with substantially reduced recombination and low levels of degeneration. We also found evidence that selective sweeps took place both on ZAL2m and its standard counterpart, ZAL2, after the rearrangement event. Signatures of positive selection were associated with allelic bias in gene expression, suggesting that antagonistic selection has operated on gene regulation. Finally, we discovered a region exhibiting long-range haplotypes inside the rearrangement on ZAL2m. These haplotypes appear to have been maintained by balancing selection, retaining genetic diversity within the supergene. Together, our analyses illuminate mechanisms contributing to the evolution of a young chromosomal polymorphism, revealing complex selective processes acting concurrently with genetic degeneration to drive the evolution of supergenes.

Radiation and hybridization underpin the spread of the fire ant social supergene

Supergenes are clusters of tightly linked genes that jointly produce complex phenotypes. Although widespread in nature, how such genomic elements are formed and how they spread are in most cases unclear. In the fire ant Solenopsis invicta and closely related species, a "social supergene controls whether a colony maintains one or multiple queens. Here, we show that the three inversions constituting the Social b (Sb) supergene emerged sequentially during the separation of the ancestral lineages of S. invicta and Solenopsis richteri. The two first inversions arose in the ancestral population of both species, while the third one arose in the S. richteri lineage. Once completely assembled in the S. richteri lineage, the supergene first introgressed into S. invicta, and from there into the other species of the socially polymorphic group of South American fire ant species. Surprisingly, the introgression of this large and important genomic element occurred despite recent hybridization being uncommon between several of the species. These results highlight how supergenes can readily move across species boundaries, possibly because of fitness benefits they provide and/or expression of selfish properties favoring their transmission.

Genomic architecture of supergenes: connecting form and function

Supergenes are tightly linked sets of loci that are inherited together and control complex phenotypes. While classical supergenes-governing traits such as wing patterns in Heliconius butterflies or heterostyly in Primula-have been studied since the Modern Synthesis, we still understand very little about how they evolve and persist in nature. The genetic architecture of supergenes is a critical factor affecting their evolutionary fate, as it can change key parameters such as recombination rate and effective population size, potentially redirecting molecular evolution of the supergene in addition to the surrounding genomic region. To understand supergene evolution, we must link genomic architecture with evolutionary patterns and processes. This is now becoming possible with recent advances in sequencing technology and powerful forward computer simulations. The present theme issue brings together theoretical and empirical papers, as well as opinion and synthesis papers, which showcase the architectural diversity of supergenes and connect this to critical processes in supergene evolution, such as polymorphism maintenance and mutation accumulation. Here, we summarize those insights to highlight new ideas and methods that illuminate the path forward for the study of supergenes in nature. This article is part of the theme issue 'Genomic architecture of supergenes: causes and evolutionary consequences'.

Iterative evolution of supergene-based social polymorphism in ants

Species commonly exhibit alternative morphs, with individual fate being determined during development by either genetic factors, environmental cues or a combination thereof. Ants offer an interesting case study because many species are polymorphic in their social structure. Some colonies contain one queen while others contain many queens. This variation in queen number is generally associated with a suite of phenotypic and life-history traits, including mode of colony founding, queen lifespan, queen-worker dimorphism and colony size. The basis of this social polymorphism has been studied in five ant lineages, and remarkably social morph seems to be determined by a supergene in all cases. These 'social supergenes' tend to be large, having formed through serial inversions, and to comprise hundreds of linked genes. They have persisted over long evolutionary timescales, in multiple lineages following speciation events, and have spread between closely related species via introgression. Their evolutionary dynamics are unusually complex, combining recessive lethality, spatially variable selection, selfish genetic elements and non-random mating. Here, we synthesize the five cases of supergene-based social polymorphism in ants, highlighting interesting commonalities, idiosyncrasies and implications for the evolution of polymorphisms in general. This article is part of the theme issue 'Genomic architecture of supergenes: causes and evolutionary consequences'.