Chromosome inversions, local adaptation and speciation

A Chromosome Inversion Creates a Supergene for Sex and Colour in Lake Malawi Cichlids

Cichlid fishes have the highest rates of evolutionary turnover of sex chromosomes among vertebrates. Many large structural polymorphisms in the radiation of cichlids in Lake Malawi are associated with sex chromosomes and may also carry adaptive variation. Here, we investigate the structure and evolutionary history of an inversion polymorphism that includes both a ZW sex locus and an orange-blotch colour polymorphism in the rock-dwelling cichlid fishes of Lake Malawi. We use long-read sequencing to characterise the sequence and breakpoints of the inversion. We quantify allele frequency differences across the inversion in population samples of the genera Metriaclima and Labeotropheus. We also examine expression differences of genes in the inversion. The simple inversion spans 7 Mb and is flanked by CACTA transposons that may have catalysed the rearrangement. The region includes ~600 genes, several of which show large differences in expression. Some of these genes are candidates for the sex and colour phenotypes. This inversion is an accessible model system for studying the role of structural polymorphisms and sex chromosome turnover in the adaptive radiation of cichlids in the lakes of East Africa.

Discovery of major QTL and a massive haplotype associated with cannabinoid biosynthesis in drug-type Cannabis

Cannabis (Cannabis sativa L.), once sidelined by decades of prohibition, has now gained recognition as a multifaceted and promising plant in both medical research and commercial applications following its recent legalization. This study leverages a genome-wide association study (GWAS) on 174 drug-type Cannabis accessions from the legal Canadian market, focusing on identifying quantitative trait loci (QTL) and candidate genes associated with eleven cannabinoid traits using 282K common single-nucleotide polymorphisms. This approach aims to transform our understanding of Cannabis genetics. We have pinpointed 33 significant markers that significantly influence cannabinoid production, promising to drive the development of Cannabis varieties with specific cannabinoid profiles. Among the notable findings is a massive haplotype of ∼60 Mb on chromosome 7 in Type I (i.e., tetrahydrocannabinol [THC]-dominant) accessions, highlighting a major genetic influence on cannabinoid profiles. These insights offer valuable guidance for Cannabis breeding programs, enabling the use of precise genetic markers to select and refine promising Cannabis varieties. This approach promises to speed up the breeding process, reduce costs significantly compared to traditional methods, and ensure that the resulting Cannabis varieties are optimized for specific medical and recreational needs. This study marks a significant stride toward fully integrating Cannabis into modern agricultural practices and genetic research, paving the way for future innovations.

3D Genome Constrains Breakpoints of Inversions That Can Act as Barriers to Gene Flow in the Stickleback

DNA within the nucleus is organised into a well-regulated three-dimensional (3D) structure. However, how such 3D genome structures influence speciation processes remains largely elusive. Recent studies have shown that 3D genome structures influence mutation rates, including the occurrence of chromosomal rearrangement. For example, breakpoints of chromosomal rearrangements tend to be located at topologically associating domain (TAD) boundaries. Here, we hypothesised that TAD structures may constrain the location of chromosomal inversions and thereby shape the genomic landscape of divergence between species with ongoing gene flow, given that inversions can act as barriers to gene flow. To test this hypothesis, we used a pair of Japanese stickleback species, Gasterosteus nipponicus (Japan Sea stickleback) and G. aculeatus (three-spined stickleback). We first constructed chromosome-scale genome assemblies of both species using high fidelity long reads and high-resolution proximity ligation data and identified several chromosomal inversions. Second, via population genomic analyses, we revealed higher genetic differentiation in inverted regions than in colinear regions and no gene flow within inversions, which contrasts with the significant gene flow in colinear regions. Third, using Hi-C data, we revealed 3D genome structures of sticklebacks, delineated by A/B compartments and TADs. Finally, we found that inversion breakpoints tend to be located at TAD boundaries. Thus, our study demonstrates that the 3D genome constrains breakpoints of inversions that can act as barriers to gene flow in the stickleback. Further integration of 3D genome analyses with population genomics could provide novel insights into how the 3D genome influences speciation.

The Genetics of Host Plant Acceptance in Pea Aphids

The evolution of host-associated sympatric populations in phytophagous insects (so called "host races") connects adaptive divergence to barriers to gene flow. Pea aphid (Acyrthosiphon pisum) host races specialise on legume species, and because host plant choice leads to assortative mating, the genetic basis of host plant acceptance is key to understanding speciation. Aphids use smell and taste in their host plant selection. While chemosensory genes frequently emerge as "outliers" in genome scans, their link to plant acceptance behaviour remains unclear. We examined the genetic basis of host-associated phenotypes using an F2 cross between two pea aphid host-associated races (specialised on alfalfa-Medicago sativa- and pea-Pisum sativum), assaying behaviour on both host plants and conducting QTL and regional heritability analyses based on a high-resolution linkage map. We identified five regions of moderate effect associated with acceptance of alfalfa, two with pea acceptance and two with survival on alfalfa. Two QTLs, one for alfalfa and one for pea acceptance, are located within a large rearranged region on chromosome 1, while other QTLs linked to alfalfa acceptance and survival are in the same region on chromosome 3-linking host plant choice to fitness. These findings highlight the polygenic basis of acceptance behaviour and the role of gene clustering and chromosomal rearrangements in promoting coupling among barrier loci. We identified 60 chemosensory genes within regions connected to acceptance, 24 of which were divergent among pea aphid races in previous genome scan or gene expression analyses. Evidence linking these genes to acceptance phenotypes supports their role in determining host plant specificity and as barrier loci contributing to pea aphid speciation.

Putting Structural Variants Into Practice: The Role of Chromosomal Inversions in the Management of Marine Environments

Major threats to marine species and ecosystems include overfishing, invasive species, pollution and climate change. The changing climate not only imposes direct threats through the impacts of severe marine heatwaves, cyclones and ocean acidification but also complicates fisheries and invasive species management by driving species range shifts. The dynamic nature of these threats means that the future of our oceans will depend on the ability of species to adapt. This has led to calls for genetic interventions focussed on enhancing species' adaptive capacity, including translocations, restocking and selective breeding. Assessing the benefits and risks of such approaches requires an improved understanding of the genetic architecture of adaptive variation, not only in relation to climate-resilient phenotypes but also locally adapted populations and the fitness of hybrids. Large structural genetic variants such as chromosomal inversions play an important role in local adaptation by linking multiple adaptive loci. Consequently, inversions are likely to be particularly important when managing for adaptive capacity. However, under some circumstances, they also accumulate deleterious mutations, potentially increasing the risk of inbreeding depression. Genetic management that takes account of these dual roles on fitness is likely to be more effective at ensuring population persistence. We summarise evolutionary factors influencing adaptive and deleterious variation of inversions, review inversions found in marine taxa, and provide a framework to predict the consequences of ignoring inversions in key management scenarios. We conclude by describing practical methods to bridge the gap between evolutionary theory and practical application of inversions in conservation.

Breakpoint-chiasma interference in pericentric inversion heterokaryotypes

Heterozygous inversion breakpoints inhibit the formation of chiasmata in their vicinity, and it has been suggested that they do so through the same mechanism that also causes interference between chiasmata. In this paper, I therefore extend my earlier model of chiasma interference to account for interference between breakpoints and chiasmata in pericentric inversion heterokaryotypes. Using this model to analyze recombination and sterility datasets for Drosophila melanogaster, I find support for the hypothesis that inversion breakpoints interfere with chiasmata in the same way and to the same degree that other chiasmata do. I also find that breakpoints, like chiasmata, appear to show negative interference in the pericentromeric region, and positive interference elsewhere. I discuss the implications of these findings in light of the recent HEI10 coarsening interference hypothesis, and conclude with some remarks about the evolutionary origin of chiasma interference.

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.

Pan-genomic analysis highlights genes associated with agronomic traits and enhances genomics-assisted breeding in alfalfa

Alfalfa (Medicago sativa L.), a globally important forage crop, is valued for its high nutritional quality and nitrogen-fixing capacity. Here, we present a high-quality pan-genome constructed from 24 diverse alfalfa accessions, encompassing a wide range of genetic backgrounds. This comprehensive analysis identified 433,765 structural variations and characterized 54,002 pan-gene families, highlighting the pivotal role of genomic diversity in alfalfa domestication and adaptation. Key structural variations associated with salt tolerance and quality traits were discovered, with functional analysis implicating genes such as MsMAP65 and MsGA3ox1. Notably, overexpression of MsGA3ox1 led to a reduced stem-leaf ratio and enhanced forage quality. The integration of genomic selection and marker-assisted breeding strategies improved genomic estimated breeding values across multiple traits, offering valuable genomic resources for advancing alfalfa breeding. These findings provide insights into the genetic basis of important agronomic traits and establish a solid foundation for future crop improvement.

Genetic Differentiation is Constrained to Chromosomal Inversions and Putative Centromeres in Locally Adapted Populations With Higher Gene Flow

The impact of genome structure on adaptation is a growing focus in evolutionary biology, revealing an important role for structural variation and recombination landscapes in shaping genetic diversity across genomes and among populations. This is particularly relevant when local adaptation occurs despite gene flow, where clustering of differentiated loci can maintain locally adapted variants by reducing recombination between them. However, the limited genomic resources for nonmodel species, including reference genomes and recombination maps, have constrained our understanding of these patterns. In this study, we leverage the Atlantic silverside-a nonmodel fish with extensive local adaptation across a steep latitudinal gradient-as an ideal system to explore how genome structure influences adaptation under varying levels of gene flow, using a newly available reference genome and multiple recombination maps. Analyzing 168 genomes from four populations, we found a continuum of genome-wide differentiation increasing from south to north, reflecting higher connectivity among southern populations and reduced gene flow at northern latitudes. With increasing gene flow, the number and clustering of FST outlier loci also increased, with differentiated loci found exclusively within large haploblocks harboring inversions and smaller peaks overlapping putative centromeric regions. Notably, sequence divergence was only evident in inversions, supporting their role in adaptive divergence with gene flow, whereas centromeric regions appeared differentiated because of low recombination and diversity, with no indication of elevated divergence. Our results support the hypothesis that clustered genomic architectures evolve with high gene flow and enhance our understanding of how inversions and centromeres are linked to different evolutionary processes.

Are you with me? Co-occurrence tests from community ecology can identify positive and negative epistasis between inversions in Mimulus guttatus

Chromosomal inversions are structural genetic variants that can play a crucial role in adaptive evolution and speciation. Patterns of attraction and repulsion among unlinked inversions - whether they tend to be carried by the same or different individuals- can indicate how selection is acting on these polymorphisms. In this study, we compare analytical techniques using data from 64 inversions that segregate among 1373 F2 plants of the yellow monkeyflower Mimulus guttatus. Mendelian assortment provides a strong null hypothesis for [Formula: see text] contingency tests. Here, we show how co-occurrence metrics used in community ecology can provide additional insight regarding coupling and repulsion of inversions at genotypic level. The centered Jaccard/Tanimoto index and the affinity score describe the specific way that inversions interact to generate epistasis for plant survival. We further explore the use of network analysis to visualize inversion interactions and to identify essential third and fourth order interactions, which expand the traditional pairwise scope of the co-occurrence metrics. We suggest that a combination of different statistical approaches will provide the most complete characterization of the fitness effects, both for inversions and other polymorphisms essential to adaptation and speciation.

Karyotype evolution and speciation in Orthoptera

Karyotype evolution might fuel speciation and can thereby contribute to species diversity. To test the hypothesis that speciation and karyotype change are linked, we estimated anagenetic and cladogenetic rates of karyotype evolution as well as speciation rates in Orthoptera. We compiled the male diploid chromosome number and the number of visible chromosome arms (the fundamental number) from published sources for 1,541 species. Chromosome-associated speciation rates were estimated by jointly modelling cladogenetic and anagenetic character evolution and the phylogenetic birth-death process in a Bayesian statistical framework using a subset of 516 species from 14 families. Our findings unveiled heterogeneity among orthopteran families in the pace of karyotype evolution and whether it was linked to speciation. In 6/14 clades, we found evidence supporting speciation-associated (cladogenetic) karyotype changes, while in 6/14 clades karyotype evolution was primarily anagenetic. The remaining clades (2/14) showed uncertainty in favour of either model. We further analyzed whether flightless phenotype, and thus less mobile species, showed higher rates of karyotype evolution. We showed that the flightless phenotype is associated with the rate of chromosome loss. The finding indicates contrasting patterns of karyotype evolution within specific orthopteran lineages, thus emphasizing substantial diversity in the pace of this evolutionary process. It also implies that substantial changes in chromosome number, arising from instances of chromosomal gains and losses, are recurring events in orthopterans that are associated with reproductive isolation and speciation, at least in some groups.

Unexpected absence of a multiple-queen supergene haplotype from supercolonial populations of Formica ants

Ants exhibit many complex social organization strategies. One particularly elaborate strategy is supercoloniality, in which a colony consists of many interconnected nests (= polydomy) with many queens (= polygyny). In many species of Formica ants, an ancient queen number supergene determines whether a colony is monogyne (= headed by single queen) or polygyne. The presence of the rearranged P haplotype typically leads colonies to be polygyne. However, the presence and function of this supergene polymorphism have not been examined in supercolonial populations. Here, we use genomic data from species in the Formica rufa group to determine whether the P haplotype leads to supercoloniality. In a Formica paralugubris population, we find that nests are polygyne despite the absence of the P haplotype in workers. We find spatial genetic ancestry patterns in nests consistent with supercolonial organization. Additionally, we find that the P haplotype is also absent in workers from supercolonial Formica aquilonia and Formica aquilonia × polyctena hybrid populations but is present in some Formica polyctena workers. We conclude that the P haplotype is not necessary for supercoloniality in the Formica rufa group, despite its long-standing association with non-supercolonial polygyny across the Formica genus.

Adaptation repeatedly uses complex structural genomic variation

Structural elements are widespread across genomes, but their complexity and role in repeatedly driving local adaptation remain unclear. In this work, we use phased genome assemblies to show that adaptive divergence in cryptic color pattern in a stick insect is repeatedly underlain by structural variation, but not a simple chromosomal inversion. We found that color pattern in populations of stick insects on two mountains is associated with translocations that have also been inverted. These translocations differ in size and origin on each mountain, but they overlap partially and involve some of the same gene regions. Moreover, this structural variation is subject to divergent selection and arose without introgression between species. Our results show how the origin of structural variation provides a mechanism for repeated bouts of adaptation.

The potential of inversions to accumulate balanced sexual antagonism is supported by simulations and Drosophila experiments

Chromosomal inversion polymorphisms can be common, but the causes of their persistence are often unclear. We propose a model for the maintenance of inversion polymorphism, which requires that some variants contribute antagonistically to two phenotypes, one of which has negative frequency-dependent fitness. These conditions yield a form of frequency-dependent disruptive selection, favoring two predominant haplotypes segregating alleles that favor opposing antagonistic phenotypes. An inversion associated with one haplotype can reduce the fitness load incurred by generating recombinant offspring, reinforcing its linkage to the haplotype and enabling both haplotypes to accumulate more antagonistic variants than expected otherwise. We develop and apply a forward simulator to examine these dynamics under a tradeoff between survival and male display. These simulations indeed generate inversion-associated haplotypes with opposing sex-specific fitness effects. Antagonism strengthens with time, and can ultimately yield karyotypes at surprisingly predictable frequencies, with striking genotype frequency differences between sexes and between developmental stages. To test whether this model may contribute to well-studied yet enigmatic inversion polymorphisms in Drosophila melanogaster, we track inversion frequencies in laboratory crosses to test whether they influence male reproductive success or survival. We find that two of the four tested inversions show significant evidence for the tradeoff examined, with In(3 R)K favoring survival and In(3 L)Ok favoring male reproduction. In line with the apparent sex-specific fitness effects implied for both of those inversions, In(3 L)Ok was also found to be less costly to the viability and/or longevity of males than females, whereas In(3 R)K was more beneficial to female survival. Based on this work, we expect that balancing selection on antagonistically pleiotropic traits may provide a significant and underappreciated contribution to the maintenance of natural inversion polymorphism.

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.

The long and short of hyperdivergent regions

The increasing prevalence of genome sequencing and assembly has uncovered evidence of hyperdivergent genomic regions - loci with excess genetic diversity - in species across the tree of life. Hyperdivergent regions are often enriched for genes that mediate environmental responses, such as immunity, parasitism, and sensory perception. Especially in self-fertilizing species where the majority of the genome is homozygous, the existence of hyperdivergent regions might imply the historical action of evolutionary forces such as introgression and/or balancing selection. We anticipate that the application of new sequencing technologies, broader taxonomic sampling, and evolutionary modeling of hyperdivergent regions will provide insights into the mechanisms that generate and maintain genetic diversity within and between species.

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.

Comparative genomics provides insights into the biogeographic and biochemical diversity of meliaceous species

Meliaceous plants such as Azadirachta indica (neem) and Melia azedarach (chinaberry) contain large amounts of limonoids with unique anti-insect activities. However, genes responsible for downstream modifications of limonoids are not well known. Here, we improve the genome assemblies of neem and chinaberry to the telomere-to-telomere (T2T) level. Allopatric speciation of the two plants is confirmed by the lineage-specific inversion of chromosome 12 in the neem lineage. We further identify two BAHD-acetyltransferases (ATs) in chinaberry (MaAT8824 and MaAT1704) that catalyse acetylation at both the C-12 and C-3 hydroxyl groups of limonoids, whereas the syntenic neem copy (AiAT0635) does not possess this activity. A critical N-terminal region (SAGAVP) is crucial for the acetylation of AiAT0635, and swapping it into the MaAT8824 version (CHRSSG) can endow it with acetylation activity. Our improved genome assemblies provide insights into allopatric speciation of neem, as well as limonoid biosynthesis and chemical diversity in meliaceous plants.

The Distribution and Dispersal of Large Haploblocks in a Superspecies

Haploblocks are regions of the genome that coalesce to an ancestor as a single unit. Differentiated haplotypes in these regions can result from the accumulation of mutational differences in low-recombination chromosomal regions, especially when selective sweeps occur within geographically structured populations. We introduce a method to identify large well-differentiated haploblock regions (LHBRs), based on the variance in standardised heterozygosity (ViSHet) of single nucleotide polymorphism (SNP) genotypes among individuals, calculated across a genomic region (500 SNPs in our case). We apply this method to the greenish warbler (Phylloscopus trochiloides) ring species, using a newly assembled reference genome and genotypes at more than 1 million SNPs among 257 individuals. Most chromosomes carry a single distinctive LHBR, containing 4-6 distinct haplotypes that are associated with geography, enabling detection of hybridisation events and transition zones between differentiated populations. LHBRs have exceptionally low within-haplotype nucleotide variation and moderately low between-haplotype nucleotide distance, suggesting their establishment through recurrent selective sweeps at varying geographic scales. Meiotic drive is potentially a powerful mechanism of producing such selective sweeps, and the LHBRs are likely to often represent centromeric regions where recombination is restricted. Links between populations enable introgression of favoured haplotypes and we identify one haploblock showing a highly discordant distribution compared to most of the genome, being present in two distantly separated geographic regions that are at similar latitudes in both east and central Asia. Our results set the stage for detailed studies of haploblocks, including their genomic location, gene content and contribution to reproductive isolation.

Structural variation, selection, and diversification of the NPIP gene family from the human pangenome

The NPIP (nuclear pore interacting protein) gene family has expanded to high copy number in humans and African apes where it has been subject to an excess of amino acid replacement consistent with positive selection (1). Due to the limitations of short-read sequencing, NPIP human genetic diversity has been poorly understood. Using highly accurate assemblies generated from long-read sequencing as part of the human pangenome, we completely characterize 169 human haplotypes (4,665 NPIP paralogs and alleles). Of the 28 NPIP paralogs, just three (NPIPB2, B11, and B14) are fixed at a single copy, and only a single locus, B2, shows no structural variation. Four NPIP paralogs map to large segmental duplication blocks that mediate polymorphic inversions (355 kbp-1.6 Mbp) corresponding to microdeletions associated with developmental delay and autism. Haplotype-based tests of positive selection and selective sweeps identify two paralogs, B9 and B15, within the top percentile for both tests. Using full-length cDNA data from 101 tissue/cell types, we construct paralog-specific gene models and show that 56% (31/55 most abundant isoforms) have not been previously described in RefSeq. We define six distinct translation start sites and other protein structural features that distinguish paralogs, including a variable number tandem repeat that encodes a beta helix of variable size that emerged ~3.1 million years ago in human evolution. Among the 28 NPIP paralogs, we identify distinct tissue and developmental patterns of expression with only a few maintaining the ancestral testis-enriched expression. A subset of paralogs (NPIPA1, A5, A6-9, B3-5, and B12/B13) show increased brain expression. Our results suggest ongoing positive selection in the human population and rapid diversification of NPIP gene models.

Evolution of Hybrid Inviability Associated With Chromosome Fusions

Chromosomal rearrangements, such as inversions, have received considerable attention in the speciation literature due to their hampering effects on recombination. Less is known about how other rearrangements, such as chromosome fissions and fusions, can affect the evolution of reproductive isolation. Here, we use crosses between populations of the wood white butterfly (Leptidea sinapis) with different karyotypes to identify genomic regions associated with hybrid inviability. We map hybrid inviability candidate loci by contrasting allele frequencies between F2 hybrids that survived until the adult stage with individuals of the same cohort that succumbed to hybrid incompatibilities. Hybrid inviability candidate regions have high genetic differentiation between parental populations, reduced recombination rates, and are enriched near chromosome fusions. By analysing sequencing coverage, we exclude aneuploidies as a direct link between hybrid inviability and chromosome fusions. Instead, our results point to an indirect relationship between hybrid inviability and chromosome fusions, possibly related to reduced recombination in fused chromosomes. Thus, we map postzygotic isolation to chromosomal rearrangements, providing crucial empirical evidence for the idea that chromosome number differences between taxa can contribute to speciation.

Incomplete recombination suppression fuels extensive haplotype diversity in a butterfly colour pattern supergene

Supergenes can evolve when recombination-suppressing mechanisms like inversions promote co-inheritance of alleles at two or more polymorphic loci that affect a complex trait. Theory shows that such genetic architectures can be favoured under balancing selection or local adaptation in the face of gene flow, but they can also bring costs associated with reduced opportunities for recombination. These costs may in turn be offset by rare 'gene flux' between inverted and ancestral haplotypes, with a range of possible outcomes. We aimed to shed light on these processes by investigating the 'BC supergene', a large genomic region comprising multiple rearrangements associated with three distinct wing colour morphs in Danaus chrysippus, a butterfly known as the African monarch, African queen and plain tiger. Using whole-genome resequencing data from 174 individuals, we first confirm the effects of BC on wing colour pattern: background melanism is associated with SNPs in the promoter region of yellow, within an inverted subregion of the supergene, while forewing tip pattern is most likely associated with copy-number variation in a separate subregion of the supergene. We then show that haplotype diversity within the supergene is surprisingly extensive: there are at least six divergent haplotype groups that experience suppressed recombination with respect to each other. Despite high divergence between these haplotype groups, we identify an unexpectedly large number of natural recombinant haplotypes. Several of the inferred crossovers occurred between adjacent inversion 'modules', while others occurred within inversions. Furthermore, we show that new haplotype groups have arisen through recombination between two pre-existing ones. Specifically, an allele for dark colouration in the promoter of yellow has recombined into distinct haplotype backgrounds on at least two separate occasions. Overall, our findings paint a picture of dynamic evolution of supergene haplotypes, fuelled by incomplete recombination suppression.

Inversions contribute disproportionately to parallel genomic divergence in dune sunflowers

The probability of parallel genetic evolution is a function of the strength of selection and constraints imposed by genetic architecture. Inversions capture locally adapted alleles and suppress recombination between them, which limits the range of adaptive responses. In addition, the combined phenotypic effect of alleles within inversions is likely to be greater than that of individual alleles; this should further increase the contributions of inversions to parallel evolution. We tested the hypothesis that inversions contribute disproportionately to parallel genetic evolution in independent dune ecotypes of Helianthus petiolaris. We analysed habitat data and identified variables underlying parallel habitat shifts. Genotype-environment association analyses of these variables indicated parallel responses of inversions to shared selective pressures. We also confirmed larger seed size across the dunes and performed quantitative trait locus mapping with multiple crosses. Quantitative trait loci shared between locations fell into inversions more than expected by chance. We used whole-genome sequencing data to identify selective sweeps in the dune ecotypes and found that the majority of shared swept regions were found within inversions. Phylogenetic analyses of shared regions indicated that within inversions, the same allele typically was found in the dune habitat at both sites. These results confirm predictions that inversions drive parallel divergence in the dune ecotypes.

Common Ancestry of the Id Locus: Chromosomal Rearrangement and Polygenic Possibilities

The diversity in dermal pigmentation and plumage color among domestic chickens is striking, with Black Bone Chickens (BBC) particularly notable for their intense melanin hyperpigmentation. This unique trait is driven by a complex chromosomal rearrangement on chromosome 20 at the Fm locus, resulting in the overexpression of the EDN3 (a gene central to melanocyte regulation). In contrast, the inhibition of dermal pigmentation is regulated by the Id locus. Although prior studies using genetic crosses, GWAS, and gene expression analysis have investigated the genetic underpinnings of the Id locus, its precise location and functional details remain elusive. Our study aims to precisely locate the Id locus, identify associated chromosomal rearrangements and candidate genes influencing dermal pigmentation, and examine the ancestral status of the Id locus in BBC breeds. Using public genomic data from BBC and non-BBC breeds, we refined the Id locus to a ~1.6 Mb region that co-localizes with Z amplicon repeat units at the distal end of the q-arm of chromosome Z within a 10.36 Mb inversion in Silkie BBC. Phylogenetic and population structure analyses reveal that the Id locus shares a common ancestry across all BBC breeds, much like the Fm locus. Selection signatures and highly differentiated BBC-specific SNPs within the MTAP gene position it as the prime candidate for the Id locus with CCDC112 and additional genes, suggesting a possible polygenic nature. Our results suggest that the Id locus is shared among BBC breeds and may function as a supergene cluster in shank and dermal pigmentation variation.

Inversions Can Accumulate Balanced Sexual Antagonism: Evidence from Simulations and Drosophila Experiments

Chromosomal inversion polymorphisms can be common, but the causes of their persistence are often unclear. We propose a model for the maintenance of inversion polymorphism, which requires that some variants contribute antagonistically to two phenotypes, one of which has negative frequency-dependent fitness. These conditions yield a form of frequency-dependent disruptive selection, favoring two predominant haplotypes segregating alleles that favor opposing antagonistic phenotypes. An inversion associated with one haplotype can reduce the fitness load incurred by generating recombinant offspring, reinforcing its linkage to the haplotype and enabling both haplotypes to accumulate more antagonistic variants than expected otherwise. We develop and apply a forward simulator to examine these dynamics under a tradeoff between survival and male display. These simulations indeed generate inversion-associated haplotypes with opposing sex-specific fitness effects. Antagonism strengthens with time, and can ultimately yield karyotypes at surprisingly predictable frequencies, with striking genotype frequency differences between sexes and between developmental stages. To test whether this model may contribute to well-studied yet enigmatic inversion polymorphisms in Drosophila melanogaster, we track inversion frequencies in laboratory crosses to test whether they influence male reproductive success or survival. We find that two of the four tested inversions show significant evidence for the tradeoff examined, with ln(3R)K favoring survival and ln(3L)Ok favoring male reproduction. In line with the apparent sex-specific fitness effects implied for both of those inversions, ln(3L)Ok was also found to be less costly to the viability and/or longevity of males than females, whereas ln(3R)K was more beneficial to female survival. Based on this work, we expect that balancing selection on antagonistically pleiotropic traits may provide a significant and underappreciated contribution to the maintenance of natural inversion polymorphism.

Evolutionary adaptation under climate change: Aedes sp. demonstrates potential to adapt to warming

Climate warming is expected to shift the distributions of mosquitoes and mosquito-borne diseases, promoting expansions at cool range edges and contractions at warm range edges. However, whether mosquito populations could maintain their warm edges through evolutionary adaptation remains unknown. Here, we investigate the potential for thermal adaptation in Aedes sierrensis, a congener of the major disease vector species that experiences large thermal gradients in its native range, by assaying tolerance to prolonged and acute heat exposure, and its genetic basis in a diverse, field-derived population. We found pervasive evidence of heritable genetic variation in mosquito heat tolerance, and phenotypic trade-offs in tolerance to prolonged versus acute heat exposure. Further, we found genomic variation associated with prolonged heat tolerance was clustered in several regions of the genome, suggesting the presence of larger structural variants such as chromosomal inversions. A simple evolutionary model based on our data estimates that the maximum rate of evolutionary adaptation in mosquito heat tolerance will exceed the projected rate of climate warming, implying the potential for mosquitoes to track warming via genetic adaptation.

Rapid Parallel Adaptation in Distinct Invasions of Ambrosia Artemisiifolia Is Driven by Large-Effect Structural Variants

When introduced to multiple distinct ranges, invasive species provide a compelling natural experiment for understanding the repeatability of adaptation. Ambrosia artemisiifolia is an invasive, noxious weed, and chief cause of hay fever. Leveraging over 400 whole-genome sequences spanning the native-range in North America and 2 invasions in Europe and Australia, we inferred demographically distinct invasion histories on each continent. Despite substantial differences in genetic source and effective population size changes during introduction, scans of both local climate adaptation and divergence from the native-range revealed genomic signatures of parallel adaptation between invasions. Disproportionately represented among these parallel signatures are 37 large haploblocks-indicators of structural variation-that cover almost 20% of the genome and exist as standing genetic variation in the native-range. Many of these haploblocks are associated with traits important for adaptation to local climate, like size and the timing of flowering, and have rapidly reformed native-range clines in invaded ranges. Others show extreme frequency divergence between ranges, consistent with a response to divergent selection on different continents. Our results demonstrate the key role of large-effect standing variants in rapid adaptation during range expansion, a pattern that is robust to diverse invasion histories.

Population Genomics of Adaptive Radiation

Adaptive radiations are rich laboratories for exploring, testing, and understanding key theories in evolution and ecology because they offer spectacular displays of speciation and ecological adaptation. Particular challenges to the study of adaptive radiation include high levels of species richness, rapid speciation, and gene flow between species. Over the last decade, high-throughput sequencing technologies and access to population genomic data have lessened these challenges by enabling the analysis of samples from many individual organisms at whole-genome scales. Here we review how population genomic data have facilitated our knowledge of adaptive radiation in five key areas: (1) phylogenetics, (2) hybridization, (3) timing and rates of diversification, (4) the genomic basis of trait evolution, and (5) the role of genome structure in divergence. We review current knowledge in each area, highlight outstanding questions, and focus on methods that facilitate detection of complex patterns in the divergence and demography of populations through time. It is clear that population genomic data are revolutionising the ability to reconstruct evolutionary history in rapidly diversifying clades. Additionally, studies are increasingly emphasising the central role of gene flow, re-use of standing genetic variation during adaptation, and structural genomic elements as facilitators of the speciation process in adaptive radiations. We highlight hybridization-and the hypothesized processes by which it shapes diversification-and questions seeking to bridge the divide between microevolutionary and macroevolutionary processes as rich areas for future study. Overall, access to population genomic data has facilitated an exciting era in adaptive radiation research, with implications for deeper understanding of fundamental evolutionary processes across the tree of life.

Opening the species box: what parsimonious microscopic models of speciation have to say about macroevolution

In the last two decades, lineage-based models of diversification, where species are viewed as particles that can divide (speciate) or die (become extinct) at rates depending on some evolving trait, have been very popular tools to study macroevolutionary processes. Here, we argue that this approach cannot be used to break down the inner workings of species diversification and that "opening the species box" is necessary to understand the causes of macroevolution, but that too detailed speciation models also fail to make robust macroevolutionary predictions. We set up a general framework for parsimonious models of speciation that rely on a minimal number of mechanistic principles: (a) reproductive isolation is caused by excessive dissimilarity between genotypes; (b) dissimilarity results from a balance between differentiation processes and homogenizing processes; and (c) dissimilarity can feed back on these processes by decelerating homogenization. We classify such models according to the main homogenizing process: (a) clonal evolution models (ecological drift), (b) models of genetic isolation (gene flow), and (c) models of isolation by distance (spatial drift). We review these models and their specific predictions on macroscopic variables such as species abundances, speciation rates, interfertility relationships, or phylogenetic tree structure. We propose new avenues of research by displaying conceptual questions remaining to be solved and new models to address them: the failure of speciation at secondary contact, the feedback of dissimilarity on homogenization, and the emergence in space of breeding barriers.

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.

Divergence and gene flow history at two large chromosomal inversions underlying ecotype differentiation in the long-snouted seahorse

Chromosomal inversions can play an important role in divergence and reproductive isolation by building and maintaining distinct allelic combinations between evolutionary lineages. Alternatively, they can take the form of balanced polymorphisms that segregate within populations until one arrangement becomes fixed. Many questions remain about how inversion polymorphisms arise, how they are maintained over the long term, and ultimately, whether and how they contribute to speciation. The long-snouted seahorse (Hippocampus guttulatus) is genetically subdivided into geographic lineages and marine-lagoon ecotypes, with shared structural variation underlying lineage and ecotype divergence. Here, we aim to characterize structural variants and to reconstruct their history and suspected role in ecotype formation. We generated a near chromosome-level genome assembly and described genome-wide patterns of diversity and divergence through the analysis of 112 whole-genome sequences from Atlantic, Mediterranean, and Black Sea populations. By also analysing linked-read sequencing data, we found evidence for two chromosomal inversions that were several megabases in length and showed contrasting allele frequency patterns between lineages and ecotypes across the species range. We reveal that these inversions represent ancient intraspecific polymorphisms, one likely being maintained by divergent selection and the other by pseudo-overdominance. A possible selective coupling between the two inversions was further supported by the absence of specific haplotype combinations and a putative functional interaction between the two inversions in reproduction. Lastly, we detected gene flux eroding divergence between inverted alleles at varying levels for the two inversions, with a likely impact on their dynamics and contribution to divergence and speciation.

Evolutionary history of an island endemic, the Azorean common quail

Oceanic islands are characterized by conditions that favour diversification into endemic lineages that can be very different from their mainland counterparts. This can be the result of fast phenotypic divergence due to drift or the result of slower adaptation to local conditions. This uniqueness can obscure their evolutionary history. Here we used morphological, stable isotope, genetic and genomic data to characterize common quails (Coturnix coturnix) in the Azores archipelago and assess the divergence from neighbouring common quail populations. Historical documents suggested that these quails could have a recent origin associated with the arrival of humans in the last centuries. Our results show that Azorean quails constitute a well-differentiated lineage with small size and dark throat pigmentation that has lost the migratory ability and that diverged from mainland quail lineages more than 0.8 mya, contrary to the notion of a recent human-mediated arrival. Even though some Azorean quails carry an inversion that affects 115 Mbp of chromosome 1 and that has been associated with the loss of the migratory behaviour in other common quail populations, half of the analysed individuals do not have that inversion and still do not migrate. The long coexistence and evolution in isolation in the Azores of two chromosomal variants (with and without the inversion) is best explained by balancing selection. Thus, a unique and long evolutionary history led to the island endemic that we know today, C. c. conturbans.

Chromosomal inversion polymorphisms are widespread across the species ranges of rough periwinkles (Littorina saxatilis and L. arcana)

Inversions are thought to play a key role in adaptation and speciation, suppressing recombination between diverging populations. Genes influencing adaptive traits cluster in inversions, and changes in inversion frequencies are associated with environmental differences. However, in many organisms, it is unclear if inversions are geographically and taxonomically widespread. The intertidal snail, Littorina saxatilis, is one such example. Strong associations between putative polymorphic inversions and phenotypic differences have been demonstrated between two ecotypes of L. saxatilis in Sweden and inferred elsewhere, but no direct evidence for inversion polymorphism currently exists across the species range. Using whole genome data from 107 snails, most inversion polymorphisms were found to be widespread across the species range. The frequencies of some inversion arrangements were significantly different among ecotypes, suggesting a parallel adaptive role. Many inversions were also polymorphic in the sister species, L. arcana, hinting at an ancient origin.

Orthopteran Neo-Sex Chromosomes Reveal Dynamics of Recombination Suppression and Evolution of Supergenes

The early evolution of sex chromosomes has remained obscure for more than a century. The Vandiemenella viatica species group of morabine grasshoppers is highly suited for studying the early stages of sex chromosome divergence and degeneration of the Y chromosome. This stems from the fact that neo-XY sex chromosomes have independently evolved multiple times by X-autosome fusions with different autosomes. Here, we generated new chromosome-level assemblies for two chromosomal races representing karyotypes with and without neo-sex chromosomes (P24XY and P24X0), and sequence data of a third chromosomal race with a different neo-XY chromosome system (P25XY). Interestingly, these two neo-XY chromosomal races are formed by different X-autosome fusions (involving chr1 and chrB, respectively), and we found that both neo-Y chromosomes have partly ceased to recombine with their neo-X counterpart. We show that the neo-XY chromosomes have diverged through accumulation of SNPs and structural mutations, and that many neo-Y-linked genes have degenerated since recombination ceased. However, the non-recombining regions of neo-Y chromosomes host non-degenerated genes crucial for sex determination, such as sex-lethal and transformer, alongside genes associated with spermatogenesis, fertility, and reproduction, illustrating their integrative role as a masculinizing supergene. Contrary to expectations, the neo-Y chromosomes showed (slightly) lower density of transposable elements (TEs) compared to other genomic regions. The study reveals the unique dynamics of young sex chromosomes, with evolution of recombination suppression and pronounced decay of (some) neo-sex chromosome genes, and provides a compelling case illustrating how chromosomal fusions and post-fusion mutational processes contribute to the evolution of supergenes.

Conserved islands of divergence associated with adaptive variation in sockeye salmon are maintained by multiple mechanisms

Local adaptation is facilitated by loci clustered in relatively few regions of the genome, termed genomic islands of divergence. The mechanisms that create and maintain these islands and how they contribute to adaptive divergence is an active research topic. Here, we use sockeye salmon as a model to investigate both the mechanisms responsible for creating islands of divergence and the patterns of differentiation at these islands. Previous research suggested that multiple islands contributed to adaptive radiation of sockeye salmon. However, the low-density genomic methods used by these studies made it difficult to fully elucidate the mechanisms responsible for islands and connect genotypes to adaptive variation. We used whole genome resequencing to genotype millions of loci to investigate patterns of genetic variation at islands and the mechanisms that potentially created them. We discovered 64 islands, including 16 clustered in four genomic regions shared between two isolated populations. Characterisation of these four regions suggested that three were likely created by structural variation, while one was created by processes not involving structural variation. All four regions were small (< 600 kb), suggesting low recombination regions do not have to span megabases to be important for adaptive divergence. Differentiation at islands was not consistently associated with established population attributes. In sum, the landscape of adaptive divergence and the mechanisms that create it are complex; this complexity likely helps to facilitate fine-scale local adaptation unique to each population.

Genomic insights into variation in thermotolerance between hybridizing swordtail fishes

Understanding how organisms adapt to changing environments is a core focus of research in evolutionary biology. One common mechanism is adaptive introgression, which has received increasing attention as a potential route to rapid adaptation in populations struggling in the face of ecological change, particularly global climate change. However, hybridization can also result in deleterious genetic interactions that may limit the benefits of adaptive introgression. Here, we used a combination of genome-wide quantitative trait locus mapping and differential gene expression analyses between the swordtail fish species Xiphophorus malinche and X. birchmanni to study the consequences of hybridization on thermotolerance. While these two species are adapted to different thermal environments, we document a complicated architecture of thermotolerance in hybrids. We identify a region of the genome that contributes to reduced thermotolerance in individuals heterozygous for X. malinche and X. birchmanni ancestry, as well as widespread misexpression in hybrids of genes that respond to thermal stress in the parental species, particularly in the circadian clock pathway. We also show that a previously mapped hybrid incompatibility between X. malinche and X. birchmanni contributes to reduced thermotolerance in hybrids. Together, our results highlight the challenges of understanding the impact of hybridization on complex ecological traits and its potential impact on adaptive introgression.

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.

Human Impacts on Great Lakes Walleye Sander vitreus Structure, Diversity and Local Adaptation

Artificial propagation and wild release may influence the genetic integrity of wild populations. This practice has been prevalent in fisheries for centuries and is often termed 'stocking'. In the Laurentian Great Lakes (Great Lakes here-on), walleye populations faced declines from the 1950s to the 1970s, prompting extensive stocking efforts for restoration. By the mid-2010s, walleye populations showed signs of recovery, but the genetic legacy of stocking on population structure at the genomic level remains unclear. Using a dataset of 45,600 genome-aligned SNP loci genotyped in 1075 walleye individuals, we investigated the genetic impacts of over 50 years of stocking across the Great Lakes. Population structure was associated with both natural geographic barriers and stocking from non-native sources. Admixture between Lake Erie walleye and walleye from the re-populated Tittabawassee River indicate that stocking may have re-distributed putatively adaptive alleles around the Great Lakes. Genome scans identified FST outliers and evidence of selective sweeps, indicating local adaptation of spawning populations is likely. Notably, one genomic region showed strong differentiation between Muskegon River and walleye from the Tittabawassee River, which was re-populated by Muskegon strain walleye, suggesting admixture and selection both impact the observed genetic diversity. Overall, our study underscores how artificial propagation and translocations can significantly alter the evolutionary trajectory of populations. The findings highlight the complex interplay between stocking practices and population genetic diversity, emphasising the need for careful management strategies to preserve the genetic integrity of wild populations amidst conservation efforts.

Admixture With Cultivated Sunflower Likely Facilitated Establishment and Spread of Wild Sunflower (Helianthus annuus) in Argentina

A better understanding of the genetic and ecological factors underlying successful invasions is critical to mitigate the negative impacts of invasive species. Here, we study the invasion history of Helianthus annuus populations from Argentina, with particular emphasis on the role of post-introduction admixture with cultivated sunflower (also H. annuus) and climate adaptation driven by large haploblocks. We conducted genotyping-by-sequencing of samples of wild populations as well as Argentinian cultivars and compared them with wild (including related annual Helianthus species) and cultivated samples from the native range. We also characterised samples for 11 known haploblocks associated with environmental variation in native populations to test whether haploblocks contributed to invasion success. Population genomics analyses supported two independent geographic sources for Argentinian populations, the central United States and Texas, but no significant contribution of related annual Helianthus species. We found pervasive admixture with cultivated sunflower, likely as result of post-introduction hybridization. Genomic scans between invasive populations and their native sources identified multiple genomic regions of divergence, possibly indicative of selection, in the invaded range. These regions significantly overlapped between the two native-invasive comparisons and showed disproportionally high crop ancestry, suggesting that crop alleles contributed to invasion success. We did not find evidence of climate adaptation mediated by haploblocks, yet outliers of genome scans were enriched in haploblock regions and, for at least two haploblocks, the cultivar haplotype was favoured in Argentina. Our results show that admixture with cultivated sunflower played a major role in the establishment and spread of H. annuus populations in Argentina.

Inversions encounter relaxed genetic constraints and balance birth and death of TPS genes in Curcuma

Evolutionary dynamics of inversion and its impact on biochemical traits are a puzzling question. Here, we show abundance of inversions in three Curcuma species (turmeric, hidden ginger and Siam tulip). Genes within inversions display higher long terminal repeat content and lower expression level compared with genomic background, suggesting inversions in Curcuma experience relaxed genetic constraints. It is corroborated by depletion of selected SNPs and enrichment of deleterious mutations in inversions detected among 56 Siam tulip cultivars. Functional verification of tandem duplicated terpene synthase (TPS) genes reveals that genes within inversions become pseudogenes, while genes outside retain catalytic function. Our findings suggest that inversions act as a counteracting force against tandem duplication in balancing birth and death of TPS genes and modulating terpenoid contents in Curcuma. This study provides an empirical example that inversions are likely not adaptive but affect biochemical traits.

De novo genome assemblies of two cryptodiran turtles with ZZ/ZW and XX/XY sex chromosomes provide insights into patterns of genome reshuffling and uncover novel 3D genome folding in amniotes

Understanding the evolution of chromatin conformation among species is fundamental to elucidate the architecture and plasticity of genomes. Nonrandom interactions of linearly distant loci regulate gene function in species-specific patterns, affecting genome function, evolution, and, ultimately, speciation. Yet, data from nonmodel organisms are scarce. To capture the macroevolutionary diversity of vertebrate chromatin conformation, here we generate de novo genome assemblies for two cryptodiran (hidden-neck) turtles via Illumina sequencing, chromosome conformation capture, and RNA-seq: Apalone spinifera (ZZ/ZW, 2n = 66) and Staurotypus triporcatus (XX/XY, 2n = 54). We detected differences in the three-dimensional (3D) chromatin structure in turtles compared to other amniotes beyond the fusion/fission events detected in the linear genomes. Namely, whole-genome comparisons revealed distinct trends of chromosome rearrangements in turtles: (1) a low rate of genome reshuffling in Apalone (Trionychidae) whose karyotype is highly conserved when compared to chicken (likely ancestral for turtles), and (2) a moderate rate of fusions/fissions in Staurotypus (Kinosternidae) and Trachemys scripta (Emydidae). Furthermore, we identified a chromosome folding pattern that enables "centromere-telomere interactions" previously undetected in turtles. The combined turtle pattern of "centromere-telomere interactions" (discovered here) plus "centromere clustering" (previously reported in sauropsids) is novel for amniotes and it counters previous hypotheses about amniote 3D chromatin structure. We hypothesize that the divergent pattern found in turtles originated from an amniote ancestral state defined by a nuclear configuration with extensive associations among microchromosomes that were preserved upon the reshuffling of the linear genome.

Pervasive fitness trade-offs revealed by rapid adaptation in large experimental populations of Drosophila melanogaster

Life-history trade-offs are an inherent feature of organismal biology that evolutionary theory posits play a key role in patterns of divergence within and between species. Efforts to quantify trade-offs are largely confined to phenotypic measurements and the identification of negative genetic-correlations among fitness-relevant traits. Here, we use time-series genomic data collected during experimental evolution in large, genetically diverse populations of Drosophila melanogaster to directly measure the manifestation of trade-offs in response to temporally fluctuating selection pressures on ecological timescales. Specifically, we quantify the genome-wide signal of antagonistic pleiotropy suggestive of trade-offs between reproduction and stress tolerance. We further identify a putative role of two cosmopolitan inversions in these trade-offs, and show that loci responding to selection during lab-based, reproduction selection exhibit signals of fluctuating selection in an outdoor mesocosm exposed to natural environmental conditions. Our results demonstrate the utility of time-series genomic data in revealing the presence and genomic architecture underlying fitness trade-offs, and add credence to models positing a role of generic life history trade-offs in the maintenance of variation in natural populations.

The origin and maintenance of supergenes contributing to ecological adaptation in Atlantic herring

Chromosomal inversions are associated with local adaptation in many species. However, questions regarding how they are formed, maintained and impact various other evolutionary processes remain elusive. Here, using a large genomic dataset of long-read and short-read sequencing, we ask these questions in one of the most abundant vertebrates on Earth, the Atlantic herring. This species has four megabase-sized inversions associated with ecological adaptation that correlate with water temperature. The S and N inversion alleles at these four loci dominate in the southern and northern parts, respectively, of the species distribution in the North Atlantic Ocean. By determining breakpoint coordinates of the four inversions and the structural variations surrounding them, we hypothesize that these inversions are formed by ectopic recombination between duplicated sequences immediately outside of the inversions. We show that these are old inversions (>1 MY), albeit formed after the split between the Atlantic herring and its sister species, the Pacific herring. There is evidence for extensive gene flux between inversion alleles at all four loci. The large Ne of herring combined with the common occurrence of opposite homozygotes across the species distribution has allowed effective purifying selection to prevent the accumulation of genetic load and repeats within the inversions.

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.

Genomics of natural populations: gene conversion events reveal selected genes within the inversions of Drosophila pseudoobscura

When adaptive phenotypic variation or quantitative trait loci map within an inverted segment of a chromosome, researchers often despair because the suppression of crossing over will prevent the discovery of selective target genes that established the rearrangement. If an inversion polymorphism is old enough, then the accumulation of gene conversion tracts offers the promise that quantitative trait loci or selected loci within inversions can be mapped. The inversion polymorphism of Drosophila pseudoobscura is a model system to show that gene conversion analysis is a useful tool for mapping selected loci within inversions. D. pseudoobscura has over 30 different chromosomal arrangements on the third chromosome (Muller C) in natural populations and their frequencies vary with changes in environmental habitats. Statistical tests of five D. pseudoobscura gene arrangements identified outlier genes within inverted regions that had potentially heritable variation, either fixed amino acid differences or differential expression patterns. We use genome sequences of the inverted third chromosome (Muller C) to infer 98,443 gene conversion tracts for a total coverage of 142 Mb or 7.2× coverage of the 19.7 Mb chromosome. We estimated gene conversion tract coverage in the 2,668 genes on Muller C and tested whether gene conversion coverage was similar among arrangements for outlier vs non-outlier loci. Outlier genes had lower gene conversion tract coverage among arrangements than the non-outlier genes suggesting that selection removes exchanged DNA in the outlier genes. These data support the hypothesis that the third chromosome in D. pseudoobscura captured locally adapted combinations of alleles prior to inversion mutation events.