Revisiting the Impact of Inversions in Evolution: From Population Genetic Markers to Drivers of Adaptive Shifts and Speciation?

Comparative karyotype analysis of the red brocket deer (M. americana sensu lato and M. rufa) complex: evidence of drastic chromosomal evolution and implications on speciation process

Chromosomal rearrangements are often associated with playing a role in the speciation process. However, the underlying mechanism that favors the genetic isolation associated with chromosomal changes remains elusive. In this sense, the genus Mazama is recognized by its high level of karyotype diversity among species with similar morphology. A cryptic species complex has been identified within the genus, with the red brocket deer (Mazama americana and Mazama rufa) being the most impressive example. The chromosome variation was clustered in cytotypes with diploid numbers ranging from 42 to 53 and was correlated with geographical location. We conducted an analysis of chromosome evolution of the red brocket deer complex using comparative chromosome painting and Bacterial Artificial Chromosome (BAC) clones among different cytotypes. The aim was to deepen our understanding of the karyotypic relationships within the red brocket, thereby elucidating the significant chromosome variation among closely related species. This underscores the significance of chromosome changes as a key evolutionary process shaping their genomes. The results revealed the presence of three distinct cytogenetic lineages characterized by significant karyotypic divergence, suggesting the existence of efficient post-zygotic barriers. Tandem fusions constitute the main mechanism driving karyotype evolution, following a few centric fusions, inversion X-autosomal fusions. The BAC mapping has improved our comprehension of the karyotypic relationships within the red brocket deer complex, prompting questions regarding the role of these changes in the speciation process. We propose the red brocket as a model group to investigate how chromosomal changes contribute to isolation and explore the implications of these changes in taxonomy and conservation.

Chromosomal inversions harbour excess mutational load in the coral, Acropora kenti, on the Great Barrier Reef

The future survival of coral reefs in the Anthropocene depends on the capacity of corals to adapt as oceans warm and extreme weather events become more frequent. Targeted interventions designed to assist evolutionary processes in corals require a comprehensive understanding of the distribution and structure of standing variation, however, efforts to map genomic variation in corals have so far focussed almost exclusively on SNPs, overlooking structural variants that have been shown to drive adaptive processes in other taxa. Here, we show that the reef-building coral, Acropora kenti, harbours at least five large, highly polymorphic structural variants, all of which exhibit signatures of strongly suppressed recombination in heterokaryotypes, a feature commonly associated with chromosomal inversions. Based on their high minor allele frequency, uniform distribution across habitats and elevated genetic load, we propose that these inversions in A. kenti are likely to be under balancing selection. An excess of SNPs with high impact on protein-coding genes within these loci elevates their importance both as potential targets for adaptive selection and as contributors to genetic decline if coral populations become fragmented or inbred in future.

Why do sex chromosomes progressively lose recombination?

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

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 In(3R)K favoring survival and In(3L)Ok favoring male reproduction. In line with the apparent sex-specific fitness effects implied for both of those inversions, In(3L)Ok was also found to be less costly to the viability and/or longevity of males than females, whereas In(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.

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.

From sub-Saharan Africa to China: Evolutionary history and adaptation of Drosophila melanogaster revealed by population genomics

Drosophila melanogaster is a widely used model organism for studying environmental adaptation. However, the genetic diversity of populations in Asia is poorly understood, leaving a notable gap in our knowledge of the global evolution and adaptation of this species. We sequenced genomes of 292 D. melanogaster strains from various ecological settings in China and analyzed them along with previously published genome sequences. We have identified six global genetic ancestry groups, despite the presence of widespread genetic admixture. The strains from China represent a unique ancestry group, although detectable differentiation exists among populations within China. We deciphered the global migration and demography of D. melanogaster, and identified widespread signals of adaptation, including genetic changes in response to insecticides. We validated the effects of insecticide resistance variants using population cage trials and deep sequencing. This work highlights the importance of population genomics in understanding the genetic underpinnings of adaptation, an effort that is particularly relevant given the deterioration of ecosystems.

Speciation across biomes: Rapid diversification with reproductive isolation in the Australian delicate mice

Phylogeographic studies of continental clades, especially when combined with palaeoclimate modelling, provide powerful insight into how environment drives speciation across climatic contexts. Australia, a continent characterized by disparate modern biomes and dynamic climate change, provides diverse opportunity to reconstruct the impact of past and present environments on diversification. Here, we use genomic-scale data (1310 exons and whole mitogenomes from 111 samples) to investigate Pleistocene diversification, cryptic diversity, and secondary contact in the Australian delicate mice (Hydromyini: Pseudomys), a recent radiation spanning almost all Australian environments. Across northern Australia, we find no evidence for introgression between cryptic lineages within Pseudomys delicatulus sensu lato, with palaeoclimate models supporting contraction and expansion of suitable habitat since the last glacial maximum. Despite multiple contact zones, we also find little evidence of introgression at a continental scale, with the exception of a potential hybrid zone in the mesic biome. In the arid zone, combined insights from genetic data and palaeomodels support a recent expansion in the arid specialist P. hermannsburgensis and contraction in the semi-arid P. bolami. In the face of repeated secondary contact, differences in sperm morphology and chromosomal rearrangements are potential mechanisms that maintain species boundaries in these recently diverged species. Additionally, we describe the western delicate mouse as a new species and recommend taxonomic reinstatement of the eastern delicate mouse. Overall, we show that speciation in an evolutionarily young and widespread clade has been driven by environmental change, and potentially maintained by divergence in reproductive morphology and chromosome rearrangements.

Diverse pathways to speciation revealed by marine snails

Speciation is a key evolutionary process that is not yet fully understood. Combining population genomic and ecological data from multiple diverging pairs of marine snails (Littorina) supports the search for speciation mechanisms. Placing pairs on a one-dimensional speciation continuum, from undifferentiated populations to species, obscured the complexity of speciation. Adding multiple axes helped to describe either speciation routes or reproductive isolation in the snails. Divergent ecological selection repeatedly generated barriers between ecotypes, but appeared less important in completing speciation while genetic incompatibilities played a key role. Chromosomal inversions contributed to genomic barriers, but with variable impact. A multidimensional (hypercube) approach supported framing of questions and identification of knowledge gaps and can be useful to understand speciation in many other systems.

The emerging importance of cross-ploidy hybridisation and introgression

Natural hybridisation is now recognised as pervasive in its occurrence across the Tree of Life. Resurgent interest in natural hybridisation fuelled by developments in genomics has led to an improved understanding of the genetic factors that promote or prevent species cross-mating. Despite this body of work overturning many widely held assumptions about the genetic barriers to hybridisation, it is still widely thought that ploidy differences between species will be an absolute barrier to hybridisation and introgression. Here, we revisit this assumption, reviewing findings from surveys of polyploidy and hybridisation in the wild. In a case study in the British flora, 203 hybrids representing 35% of hybrids with suitable data have formed via cross-ploidy matings, while a wider literature search revealed 59 studies (56 in plants and 3 in animals) in which cross-ploidy hybridisation has been confirmed with genetic data. These results show cross-ploidy hybridisation is readily overlooked, and potentially common in some groups. General findings from these studies include strong directionality of hybridisation, with introgression usually towards the higher ploidy parent, and cross-ploidy hybridisation being more likely to involve allopolyploids than autopolyploids. Evidence for adaptive introgression across a ploidy barrier and cases of cross-ploidy hybrid speciation shows the potential for important evolutionary outcomes.

Chromosomal Inversions and the Demography of Speciation in Drosophila montana and Drosophila flavomontana

Chromosomal inversions may play a central role in speciation given their ability to locally reduce recombination and therefore genetic exchange between diverging populations. We analyzed long- and short-read whole-genome data from sympatric and allopatric populations of 2 Drosophila virilis group species, Drosophila montana and Drosophila flavomontana, to understand if inversions have contributed to their divergence. We identified 3 large alternatively fixed inversions on the X chromosome and one on each of the autosomes 4 and 5. A comparison of demographic models estimated for inverted and noninverted (colinear) chromosomal regions suggests that these inversions arose before the time of the species split. We detected a low rate of interspecific gene flow (introgression) from D. montana to D. flavomontana, which was further reduced inside inversions and was lower in allopatric than in sympatric populations. Together, these results suggest that the inversions were already present in the common ancestral population and that gene exchange between the sister taxa was reduced within inversions both before and after the onset of species divergence. Such ancestrally polymorphic inversions may foster speciation by allowing the accumulation of genetic divergence in loci involved in adaptation and reproductive isolation inside inversions early in the speciation process, while gene exchange at colinear regions continues until the evolving reproductive barriers complete speciation. The overlapping X inversions are particularly good candidates for driving the speciation process of D. montana and D. flavomontana, since they harbor strong genetic incompatibilities that were detected in a recent study of experimental introgression.

The Complex Landscape of Structural Divergence Between the Drosophila pseudoobscura and D. persimilis Genomes

Structural genomic variants are key drivers of phenotypic evolution. They can span hundreds to millions of base pairs and can thus affect large numbers of genetic elements. Although structural variation is quite common within and between species, its characterization depends upon the quality of genome assemblies and the proportion of repetitive elements. Using new high-quality genome assemblies, we report a complex and previously hidden landscape of structural divergence between the genomes of Drosophila persimilis and D. pseudoobscura, two classic species in speciation research, and study the relationships among structural variants, transposable elements, and gene expression divergence. The new assemblies confirm the already known fixed inversion differences between these species. Consistent with previous studies showing higher levels of nucleotide divergence between fixed inversions relative to collinear regions of the genome, we also find a significant overrepresentation of INDELs inside the inversions. We find that transposable elements accumulate in regions with low levels of recombination, and spatial correlation analyses reveal a strong association between transposable elements and structural variants. We also report a strong association between differentially expressed (DE) genes and structural variants and an overrepresentation of DE genes inside the fixed chromosomal inversions that separate this species pair. Interestingly, species-specific structural variants are overrepresented in DE genes involved in neural development, spermatogenesis, and oocyte-to-embryo transition. Overall, our results highlight the association of transposable elements with structural variants and their importance in driving evolutionary divergence.

Rapid and Repeated Climate Adaptation Involving Chromosome Inversions following Invasion of an Insect

Following invasion, insects can become adapted to conditions experienced in their invasive range, but there are few studies on the speed of adaptation and its genomic basis. Here, we examine a small insect pest, Thrips palmi, following its contemporary range expansion across a sharp climate gradient from the subtropics to temperate areas. We first found a geographically associated population genetic structure and inferred a stepping-stone dispersal pattern in this pest from the open fields of southern China to greenhouse environments of northern regions, with limited gene flow after colonization. In common garden experiments, both the field and greenhouse groups exhibited clinal patterns in thermal tolerance as measured by critical thermal maximum (CTmax) closely linked with latitude and temperature variables. A selection experiment reinforced the evolutionary potential of CTmax with an estimated h2 of 6.8% for the trait. We identified 3 inversions in the genome that were closely associated with CTmax, accounting for 49.9%, 19.6%, and 8.6% of the variance in CTmax among populations. Other genomic variations in CTmax outside the inversion region were specific to certain populations but functionally conserved. These findings highlight rapid adaptation to CTmax in both open field and greenhouse populations and reiterate the importance of inversions behaving as large-effect alleles in climate adaptation.

Unravelling inversions: Technological advances, challenges, and potential impact on crop breeding

Inversions, a type of chromosomal structural variation, significantly influence plant adaptation and gene functions by impacting gene expression and recombination rates. However, compared with other structural variations, their roles in functional biology and crop improvement remain largely unexplored. In this review, we highlight technological and methodological advancements that have allowed a comprehensive understanding of inversion variants through the pangenome framework and machine learning algorithms. Genome editing is an efficient method for inducing or reversing inversion mutations in plants, providing an effective mechanism to modify local recombination rates. Given the potential of inversions in crop breeding, we anticipate increasing attention on inversions from the scientific community in future research and breeding applications.

Structural Variants and Speciation: Multiple Processes at Play

Research on the genomic architecture of speciation has increasingly revealed the importance of structural variants (SVs) that affect the presence, abundance, position, and/or direction of a nucleotide sequence. SVs include large chromosomal rearrangements such as fusion/fissions and inversions and translocations, as well as smaller variants such as duplications, insertions, and deletions (CNVs). Although we have ample evidence that SVs play a key role in speciation, the underlying mechanisms differ depending on the type and length of the SV, as well as the ecological, demographic, and historical context. We review predictions and empirical evidence for classic processes such as underdominance due to meiotic aberrations and the coupling effect of recombination suppression before exploring how recent sequencing methodologies illuminate the prevalence and diversity of SVs. We discuss specific properties of SVs and their impact throughout the genome, highlighting that multiple processes are at play, and possibly interacting, in the relationship between SVs and speciation.

Adaptation to seasonal reproduction and environment-associated factors drive temporal and spatial differentiation in northwest Atlantic herring despite gene flow

Understanding how marine organisms adapt to local environments is crucial for predicting how populations will respond to global climate change. The genomic basis, environmental factors and evolutionary processes involved in local adaptation are however not well understood. Here we use Atlantic herring, an abundant, migratory and widely distributed marine fish with substantial genomic resources, as a model organism to evaluate local adaptation. We examined genomic variation and its correlation with environmental variables across a broad environmental gradient, for 15 spawning aggregations in Atlantic Canada and the United States. We then compared our results with available genomic data of northeast Atlantic populations. We confirmed that population structure lies in a fraction of the genome including likely adaptive genetic variants of functional importance. We discovered 10 highly differentiated genomic regions distributed across four chromosomes. Nine regions show strong association with seasonal reproduction. One region, corresponding to a known inversion on chromosome 12, underlies a latitudinal pattern discriminating populations north and south of a biogeographic transition zone on the Scotian Shelf. Genome-environment associations indicate that winter seawater temperature best correlates with the latitudinal pattern of this inversion. The variation at two so-called 'islands of divergence' related to seasonal reproduction appear to be private to the northwest Atlantic. Populations in the northwest and northeast Atlantic share variation at four of these divergent regions, simultaneously displaying significant diversity in haplotype composition at another four regions, which includes an undescribed structural variant approximately 7.7 Mb long on chromosome 8. Our results suggest that the timing and geographic location of spawning and early development may be under diverse selective pressures related to allelic fitness across environments. Our study highlights the role of genomic architecture, ancestral haplotypes and selection in maintaining adaptive divergence in species with large population sizes and presumably high gene flow.

Genetic architecture of ecological divergence between Oryza rufipogon and Oryza nivara

Ecological divergence due to habitat difference plays a prominent role in the formation of new species, but the genetic architecture during ecological speciation and the mechanism underlying phenotypic divergence remain less understood. Two wild ancestors of rice (Oryza rufipogon and Oryza nivara) are a progenitor-derivative species pair with ecological divergence and provide a unique system for studying ecological adaptation/speciation. Here, we constructed a high-resolution linkage map and conducted a quantitative trait locus (QTL) analysis of 19 phenotypic traits using an F2 population generated from a cross between the two Oryza species. We identified 113 QTLs associated with interspecific divergence of 16 quantitative traits, with effect sizes ranging from 1.61% to 34.1% in terms of the percentage of variation explained (PVE). The distribution of effect sizes of QTLs followed a negative exponential, suggesting that a few genes of large effect and many genes of small effect were responsible for the phenotypic divergence. We observed 18 clusters of QTLs (QTL hotspots) on 11 chromosomes, significantly more than that expected by chance, demonstrating the importance of coinheritance of loci/genes in ecological adaptation/speciation. Analysis of effect direction and v-test statistics revealed that interspecific differentiation of most traits was driven by divergent natural selection, supporting the argument that ecological adaptation/speciation would proceed rapidly under coordinated selection on multiple traits. Our findings provide new insights into the understanding of genetic architecture of ecological adaptation and speciation in plants and help effective manipulation of specific genes or gene cluster in rice breeding.

Transposable Element Insertions Are Associated with Batesian Mimicry in the Pantropical Butterfly Hypolimnas misippus

Hypolimnas misippus is a Batesian mimic of the toxic African Queen butterfly (Danaus chrysippus). Female H. misippus butterflies use two major wing patterning loci (M and A) to imitate three color morphs of D. chrysippus found in different regions of Africa. In this study, we examine the evolution of the M locus and identify it as an example of adaptive atavism. This phenomenon involves a morphological reversion to an ancestral character that results in an adaptive phenotype. We show that H. misippus has re-evolved an ancestral wing pattern present in other Hypolimnas species, repurposing it for Batesian mimicry of a D. chrysippus morph. Using haplotagging, a linked-read sequencing technology, and our new analytical tool, Wrath, we discover two large transposable element insertions located at the M locus and establish that these insertions are present in the dominant allele responsible for producing mimetic phenotype. By conducting a comparative analysis involving additional Hypolimnas species, we demonstrate that the dominant allele is derived. This suggests that, in the derived allele, the transposable elements disrupt a cis-regulatory element, leading to the reversion to an ancestral phenotype that is then utilized for Batesian mimicry of a distinct model, a different morph of D. chrysippus. Our findings present a compelling instance of convergent evolution and adaptive atavism, in which the same pattern element has independently evolved multiple times in Hypolimnas butterflies, repeatedly playing a role in Batesian mimicry of diverse model species.

Chromosomes of Asian cyprinid fishes: Novel insight into the chromosomal evolution of Labeoninae (Teleostei, Cyprinidae)

The Labeoninae subfamily is a highly diversified but demonstrably monophyletic lineage of cyprinid fishes comprising five tribes and six incertae sedis genera. This widely distributed assemblage contains some 48 genera and around 480 recognized species distributed in freshwaters of Africa and Asia. In this study, the karyotypes and other chromosomal properties of five Labeoninae species found in Thailand Labeo chrysophekadion (Labeonini) and Epalzeorhynchos bicolor, Epalzeorhynchos munense, Henicorhynchus siamensis, Thynnichthys thynnoides (´Osteochilini´) were examined using conventional and molecular cytogenetic protocols. Our results confirmed a diploid chromosome number (2n) invariably 2n = 50, but the ratio of uni- and bi-armed chromosomes was highly variable among their karyotypes, indicating extensive structural chromosomal rearrangements. Karyotype of L. chrysophekadion contained 10m+6sm+20st+14a, 32m+10sm+8st for H. siamensis, 20m+12sm+10st+8a in E. bicolor, 20m+8sm+8st+14a in E. munense, and 18m+24sm+8st in T. thynnoides. Except for H. siamensis, which had four sites of 5S rDNA sites, other species under study had only one chromosome pair with those sites. In contrast, only one pair containing 18S rDNA sites were found in the karyotypes of three species, whereas two sites were found in that of E. bicolor. These cytogenetic patterns indicated that the cytogenomic divergence patterns of these labeonine species largely corresponded to the inferred phylogenetic tree. In spite of the 2n stability, diverse patterns of rDNA and microsatellite distribution as well as their various karyotype structures demonstrated significant evolutionary differentiation of Labeoninae genomes as exemplified in examined species. Labeoninae offers a traditional point of view on the evolutionary forces fostering biological diversity, and the recent findings add new pieces to comprehend the function of structural chromosomal rearrangements in adaption and speciation.

Postzygotic barriers persist despite ongoing introgression in hybridizing Mimulus species

The evolution of postzygotic isolation is thought to be a key step in maintaining species boundaries upon secondary contact, yet the dynamics and persistence of hybrid incompatibilities in naturally hybridizing species are not well understood. Here, we explore these issues using genetic mapping in three independent populations of recombinant inbred lines between naturally hybridizing monkeyflowers, Mimulus guttatus and Mimulus nasutus, from the sympatric Catherine Creek population. We discover that the three M. guttatus founders differ dramatically in admixture history, with nearly a quarter of one founder's genome introgressed from M. nasutus. Comparative genetic mapping in the three RIL populations reveals three new putative inversions, each one segregating among the M. guttatus founders, two due to admixture. We find strong, genome-wide transmission ratio distortion in all RILs, but patterns are highly variable among the three populations. At least some of this distortion appears to be explained by epistatic selection favouring parental genotypes, but tests of inter-chromosomal linkage disequilibrium also reveal multiple candidate Dobzhansky-Muller incompatibilities. We also map several genetic loci for hybrid pollen viability, including two interacting pairs that coincide with peaks of distortion. Remarkably, even with this limited sample of three M. guttatus lines, we discover abundant segregating variation for hybrid incompatibilities with M. nasutus, suggesting this population harbours diverse contributors to postzygotic isolation. Moreover, even with substantial admixture, hybrid incompatibilities between Mimulus species persist, suggesting postzygotic isolation might be a potent force in maintaining species barriers in this system.

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.

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.

Alloteropsis semialata as a study system for C4 evolution in grasses

BACKGROUND

Numerous groups of plants have adapted to CO2 limitations by independently evolving C4 photosynthesis. This trait relies on concerted changes in anatomy and biochemistry to concentrate CO2 within the leaf and thereby boost productivity in tropical conditions. The ecological and economic importance of C4 photosynthesis has motivated intense research, often relying on comparisons between distantly related C4 and non-C4 plants. The photosynthetic type is fixed in most species, with the notable exception of the grass Alloteropsis semialata. This species includes populations exhibiting the ancestral C3 state in southern Africa, intermediate populations in the Zambezian region and C4 populations spread around the palaeotropics.

SCOPE

We compile here the knowledge on the distribution and evolutionary history of the Alloteropsis genus as a whole and discuss how this has furthered our understanding of C4 evolution. We then present a chromosome-level reference genome for a C3 individual and compare the genomic architecture with that of a C4 accession of A. semialata.

CONCLUSIONS

Alloteropsis semialata is one of the best systems in which to investigate the evolution of C4 photosynthesis because the genetic and phenotypic variation provides a fertile ground for comparative and population-level studies. Preliminary comparative genomic investigations show that the C3 and C4 genomes are highly syntenic and have undergone a modest amount of gene duplication and translocation since the different photosynthetic groups diverged. The background knowledge and publicly available genomic resources make A. semialata a great model for further comparative analyses of photosynthetic diversification.

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.

The holocentric chromosome microevolution: From phylogeographic patterns to genomic associations with environmental gradients

Geographic isolation and chromosome evolution are two of the major drivers of diversification in eukaryotes in general, and specifically, in plants. On one hand, range shifts induced by Pleistocene glacial oscillations deeply shaped the evolutionary trajectories of species in the Northern Hemisphere. On the other hand, karyotype variability within species or species complexes may have adaptive potential as different karyotypes may represent different recombination rates and linkage groups that may be associated with locally adapted genes or supergenes. Organisms with holocentric chromosomes are ideal to study the link between local adaptation and chromosome evolution, due to their high cytogenetic variability, especially when it seems to be related to environmental variation. Here, we integrate the study of the phylogeography, chromosomal evolution and ecological requirements of a plant species complex distributed in the Western Euro-Mediterranean region (Carex gr. laevigata, Cyperaceae). We aim to clarify the relative influence of these factors on population differentiation and ultimately on speciation. We obtained a well-resolved RADseq phylogeny that sheds light on the phylogeographic patterns of molecular and chromosome number variation, which are compatible with south-to-north postglacial migration. In addition, landscape genomics analyses identified candidate loci for local adaptation, and also strong significant associations between the karyotype and the environment. We conclude that karyotype distribution in C. gr. laevigata has been constrained by both range shift dynamics and local adaptation. Our study demonstrates that chromosome evolution may be responsible, at least partially, for microevolutionary patterns of population differentiation and adaptation in Carex.

Female meiotic drive shapes the distribution of rare inversion polymorphisms in Drosophila melanogaster

In all species, new chromosomal inversions are constantly being formed by spontaneous rearrangement and then stochastically eliminated from natural populations. In Drosophila, when new chromosomal inversions overlap with a preexisting inversion in the population, their rate of elimination becomes a function of the relative size, position, and linkage phase of the gene rearrangements. These altered dynamics result from complex meiotic behavior wherein overlapping inversions generate asymmetric dyads that cause both meiotic drive/drag and segmental aneuploidy. In this context, patterns in rare inversion polymorphisms of a natural population can be modeled from the fundamental genetic processes of forming asymmetric dyads via crossing-over in meiosis I and preferential segregation from asymmetric dyads in meiosis II. Here, a mathematical model of crossover-dependent female meiotic drive is developed and parameterized with published experimental data from Drosophila melanogaster laboratory constructs. This mechanism is demonstrated to favor smaller, distal inversions and accelerate the elimination of larger, proximal inversions. Simulated sampling experiments indicate that the paracentric inversions directly observed in natural population surveys of D. melanogaster are a biased subset that both maximizes meiotic drive and minimizes the frequency of lethal zygotes caused by this cytogenetic mechanism. Incorporating this form of selection into a population genetic model accurately predicts the shift in relative size, position, and linkage phase for rare inversions found in this species. The model and analysis presented here suggest that this weak form of female meiotic drive is an important process influencing the genomic distribution of rare inversion polymorphisms.

Chromosomal Radiation: A model to explain karyotypic diversity in cryptic species

We present a concept that explains the pattern of occurrence of widely distributed organisms with large chromosomal diversity, large or small molecular divergence, and the insufficiency or absence of morphological identity. Our model is based on cytogenetic studies associated with molecular and biological data and can be applied to any lineage of sister species, chronospecies, or cryptic species. Through the evaluation of the karyotypic macrostructure, as the physical location of genes e satellites DNAs, in addition to phylogenetic reconstructions from mitochondrial and nuclear genes, per example, we have observed morphologically indistinguishable individuals presenting different locally fixed karyomorphs with phylogeographic discontinuity. The biological process behind this pattern is seen in many groups of cryptic species, in which variation lies mainly in the organization of their genomes but not necessarily in the ecosystems they inhabit or in their external morphology. It's similar to the processes behind other events observed in the distribution of lineages. In this work, we explore the hypothesis of a process analogous to ecological-evolutionary radiation, which we called Chromosomal Radiation. Chromosomal Radiation can be adaptive or non-adaptive and applied to different groups of organisms.

Ecology, the pace-of-life, epistatic selection and the maintenance of genetic variation in life-history genes

Evolutionary genetics has long struggled with understanding how functional genes under selection remain polymorphic in natural populations. Taking as a starting point that natural selection is ultimately a manifestation of ecological processes, we spotlight an underemphasized and potentially ubiquitous ecological effect that may have fundamental effects on the maintenance of genetic variation. Negative frequency dependency is a well-established emergent property of density dependence in ecology, because the relative profitability of different modes of exploiting or utilizing limiting resources tends to be inversely proportional to their frequency in a population. We suggest that this may often generate negative frequency-dependent selection (NFDS) on major effect loci that affect rate-dependent physiological processes, such as metabolic rate, that are phenotypically manifested as polymorphism in pace-of-life syndromes. When such a locus under NFDS shows stable intermediate frequency polymorphism, this should generate epistatic selection potentially involving large numbers of loci with more minor effects on life-history (LH) traits. When alternative alleles at such loci show sign epistasis with a major effect locus, this associative NFDS will promote the maintenance of polygenic variation in LH genes. We provide examples of the kind of major effect loci that could be involved and suggest empirical avenues that may better inform us on the importance and reach of this process.

Genomic analyses reveal dead-end hybridization between two deeply divergent kiwifruit species rather than homoploid hybrid speciation

Despite the importance of hybridization in evolution, the evolutionary consequence of homoploid hybridizations in plants remains poorly understood. Specially, homoploid hybridization events have been rarely documented due to a lack of genomic resources and methodological limitations. Actinidia zhejiangensis was suspected to have arisen from hybridization of Actinidia eriantha and Actinidia hemsleyana or Actinidia rufa. However, this species was very rare in nature and exhibited sympatric distribution with its potential parent species, which implied it might be a spontaneous hybrid of ongoing homoploid hybridization. Here, we illustrate the dead-end homoploid hybridization and genomic basis of isolating barriers between A. eriantha and A. hemsleyana through whole genome sequencing and population genomic analyses. Chromosome-scale genome assemblies of A. zhejiangensis and A. hemsleyana were generated. The chromosomes of A. zhejiangensis are confidently assigned to the two haplomes, and one of them originates from A. eriantha and the other originates from A. hemsleyana. Whole genome resequencing data reveal that A. zhejiangensis are mainly F1 hybrids of A. hemsleyana and A. eriantha and gene flow initiated about 0.98 million years ago, implying both strong genetic barriers and ongoing hybridization between these two deeply divergent kiwifruit species. Five inversions containing genes involved in pollen germination and pollen tube growth might account for the fertility breakdown of hybrids between A. hemsleyana and A. eriantha. Despite its distinct morphological traits and long recurrent hybrid origination, A. zhejiangensis does not initiate speciation. Collectively, our study provides new insights into homoploid hybridization in plants and provides genomic resources for evolutionary and functional genomic studies of kiwifruit.

Postzygotic barriers persist despite ongoing introgression in hybridizing Mimulus species

The evolution of postzygotic isolation is thought to be a key step in maintaining species boundaries upon secondary contact, yet the dynamics and persistence of hybrid incompatibilities in sympatric species are not well understood.Here, we explore these issues using genetic mapping in three populations of recombinant inbred lines between naturally hybridizing monkeyflowers Mimulus guttatus and M. nasutus from the sympatric Catherine Creek population.The three M. guttatus founders differ dramatically in admixture history. Comparative genetic mapping also reveals three putative inversions segregating among the M. guttatus founders, two due to admixture. We observe strong, genome-wide transmission ratio distortion, but patterns are highly variable among populations. Some distortion is explained by epistatic selection favoring parental genotypes, but tests of inter-chromosomal linkage disequilibrium also reveal multiple candidate Dobzhansky-Muller incompatibilities. We also map several genetic loci for hybrid fertility, including two interacting pairs coinciding with peaks of distortion.Remarkably, in this limited sample of M. guttatus, we discover abundant segregating variation for hybrid incompatibilities with M. nasutus, suggesting this population harbors diverse contributors to postzygotic isolation. Moreover, even with substantial admixture, hybrid incompatibilities between Mimulus species persist, suggesting postzygotic isolation might be a potent force in maintaining species barriers in this system.

The Adaptive Value of Chromosomal Inversions and Climatic Change-Studies on the Natural Populations of Drosophila subobscura from the Balkans

The adaptive value of the Drosophila subobscura chromosomal inversion polymorphism with regard to environmental effects is well-known. However, the specific details of the inversion adaptations to the global warming scenario deserve to be analyzed. Toward this aim, polymorphism and karyotypes were studied in 574 individuals from Petnica (Serbia) in annual samples taken in June for the period 2019-2022. Comparing the results of Petnica (Cfa: humid subtropical climate) with those from Avala (Serbia: Cfb, temperate oceanic climate) and Font Groga (Barcelona, Spain; Csa: hot-summer Mediterranean climate), significant differences were observed for their chromosomal polymorphism. In Petnica, inversions from U and E chromosomes mainly reacted significantly with regard to temperature, humidity, and rainfall. Moreover, the inversion polymorphism from Petnica (2019-2022) was compared with that from 1995. In this period, a significant increase in mean and maximum temperature was observed. However, to properly explain the observed variations of inversions over time, it was necessary to carefully analyze annual seasonal changes and particular heat wave episodes. Interestingly, yearly fluctuations of U chromosome 'warm'-adapted inversions corresponded with opposite changes in 'non-thermal' inversions. Perhaps these types of inversions were not correctly defined with regard to thermal adaptation, or these fluctuations were also due to adaptations to other physical and/or biological variables. Finally, a joint study of chromosomal inversion polymorphism from many Balkan populations of D. subobscura indicated that different climatic regions presented distinct composition, including thermal-adapted inversions.

The genome of Lactuca saligna, a wild relative of lettuce, provides insight into non-host resistance to the downy mildew Bremia lactucae

Lactuca saligna L. is a wild relative of cultivated lettuce (Lactuca sativa L.), with which it is partially interfertile. Hybrid progeny suffer from hybrid incompatibility (HI), resulting in reduced fertility and distorted transmission ratios. Lactuca saligna displays broad-spectrum resistance against lettuce downy mildew caused by Bremia lactucae Regel and is considered a non-host species. This phenomenon of resistance in L. saligna is called non-host resistance (NHR). One possible mechanism behind this NHR is through the plant-pathogen interaction triggered by pathogen recognition receptors, including nucleotide-binding leucine-rich repeat (NLR) proteins and receptor-like kinases (RLKs). We report a chromosome-level genome assembly of L. saligna (accession CGN05327), leading to the identification of two large paracentric inversions (>50 Mb) between L. saligna and L. sativa. Genome-wide searches delineated the major resistance clusters as regions enriched in NLRs and RLKs. Three of the enriched regions co-locate with previously identified NHR intervals. RNA-seq analysis of Bremia-infected lettuce identified several differentially expressed RLKs in NHR regions. Three tandem wall-associated kinase-encoding genes (WAKs) in the NHR8 interval display particularly high expression changes at an early stage of infection. We propose RLKs as strong candidates for determinants of the NHR phenotype of L. saligna.

Are scattered microsatellites weak chromosomal markers? Guided mapping reveals new insights into Trachelyopterus (Siluriformes: Auchenipteridae) diversity

The scattered distribution pattern of microsatellites is a challenging problem in fish cytogenetics. This type of array hinders the identification of useful patterns and the comparison between species, often resulting in over-limited interpretations that only label it as "scattered" or "widely distributed". However, several studies have shown that the distribution pattern of microsatellites is non-random. Thus, here we tested whether a scattered microsatellite could have distinct distribution patterns on homeologous chromosomes of closely related species. The clustered sites of 18S and 5S rDNA, U2 snRNA and H3/H4 histone genes were used as a guide to compare the (GATA)n microsatellite distribution pattern on the homeologous chromosomes of six Trachelyopterus species: T. coriaceus and Trachelyopterus aff. galeatus from the Araguaia River basin; T. striatulus, T. galeatus and T. porosus from the Amazonas River basin; and Trachelyopterus aff. coriaceus from the Paraguay River basin. Most species had similar patterns of the (GATA)n microsatellite in the histone genes and 5S rDNA carriers. However, we have found a chromosomal polymorphism of the (GATA)n sequence in the 18S rDNA carriers of Trachelyopterus galeatus, which is in Hard-Weinberg equilibrium and possibly originated through amplification events; and a chromosome polymorphism in Trachelyopterus aff. galeatus, which combined with an inversion polymorphism of the U2 snRNA in the same chromosome pair resulted in six possible cytotypes, which are in Hardy-Weinberg disequilibrium. Therefore, comparing the distribution pattern on homeologous chromosomes across the species, using gene clusters as a guide to identify it, seems to be an effective way to further the analysis of scattered microsatellites in fish cytogenetics.

Population History Shapes Responses to Different Temperature Regimes in Drosophila subobscura

Drosophila subobscura is considered a good model species for investigation of a population's ability to adapt and cope with climate changes. Decade long research has shown that inversion frequencies change in response to environmental factors indicating their role in adaptation to novel environments. The mechanisms behind organisms' responses to temperature are complex, involving changes in physiology, behavior, gene expression and regulation. On the other hand, a population's ability to respond to suboptimal conditions depends on standing genetic variation and population history. In order to elucidate the role of local adaptation in population response to the changing temperature, we investigated the response to temperature in D. subobscura individuals originating from two different altitudes by combining traditional cytogenetic techniques with assessing the levels of Hsp70 protein expression. Inversion polymorphism was assessed in the flies sampled from natural populations and in flies reared in laboratory conditions at three different temperatures after five and sixteen generations and Hsp70 protein expression profile in 12th generation flies at the basal level and after heat shock induction. Our results indicate that local adaptation and population history influence population response to the changing temperature.

An Ancestral Balanced Inversion Polymorphism Confers Global Adaptation

Since the pioneering work of Dobzhansky in the 1930s and 1940s, many chromosomal inversions have been identified, but how they contribute to adaptation remains poorly understood. In Drosophila melanogaster, the widespread inversion polymorphism In(3R)Payne underpins latitudinal clines in fitness traits on multiple continents. Here, we use single-individual whole-genome sequencing, transcriptomics, and published sequencing data to study the population genomics of this inversion on four continents: in its ancestral African range and in derived populations in Europe, North America, and Australia. Our results confirm that this inversion originated in sub-Saharan Africa and subsequently became cosmopolitan; we observe marked monophyletic divergence of inverted and noninverted karyotypes, with some substructure among inverted chromosomes between continents. Despite divergent evolution of this inversion since its out-of-Africa migration, derived non-African populations exhibit similar patterns of long-range linkage disequilibrium between the inversion breakpoints and major peaks of divergence in its center, consistent with balancing selection and suggesting that the inversion harbors alleles that are maintained by selection on several continents. Using RNA-sequencing, we identify overlap between inversion-linked single-nucleotide polymorphisms and loci that are differentially expressed between inverted and noninverted chromosomes. Expression levels are higher for inverted chromosomes at low temperature, suggesting loss of buffering or compensatory plasticity and consistent with higher inversion frequency in warm climates. Our results suggest that this ancestrally tropical balanced polymorphism spread around the world and became latitudinally assorted along similar but independent climatic gradients, always being frequent in subtropical/tropical areas but rare or absent in temperate climates.

Heterozygous inversion breakpoints suppress meiotic crossovers by altering recombination repair outcomes

Heterozygous chromosome inversions suppress meiotic crossover (CO) formation within an inversion, potentially because they lead to gross chromosome rearrangements that produce inviable gametes. Interestingly, COs are also severely reduced in regions nearby but outside of inversion breakpoints even though COs in these regions do not result in rearrangements. Our mechanistic understanding of why COs are suppressed outside of inversion breakpoints is limited by a lack of data on the frequency of noncrossover gene conversions (NCOGCs) in these regions. To address this critical gap, we mapped the location and frequency of rare CO and NCOGC events that occurred outside of the dl-49 chrX inversion in D. melanogaster. We created full-sibling wildtype and inversion stocks and recovered COs and NCOGCs in the syntenic regions of both stocks, allowing us to directly compare rates and distributions of recombination events. We show that COs outside of the proximal inversion breakpoint are distributed in a distance-dependent manner, with strongest suppression near the inversion breakpoint. We find that NCOGCs occur evenly throughout the chromosome and, importantly, are not suppressed near inversion breakpoints. We propose a model in which COs are suppressed by inversion breakpoints in a distance-dependent manner through mechanisms that influence DNA double-strand break repair outcome but not double-strand break formation. We suggest that subtle changes in the synaptonemal complex and chromosome pairing might lead to unstable interhomolog interactions during recombination that permits NCOGC formation but not CO formation.

Pan-genome inversion index reveals evolutionary insights into the subpopulation structure of Asian rice

Understanding and exploiting genetic diversity is a key factor for the productive and stable production of rice. Here, we utilize 73 high-quality genomes that encompass the subpopulation structure of Asian rice (Oryza sativa), plus the genomes of two wild relatives (O. rufipogon and O. punctata), to build a pan-genome inversion index of 1769 non-redundant inversions that span an average of ~29% of the O. sativa cv. Nipponbare reference genome sequence. Using this index, we estimate an inversion rate of ~700 inversions per million years in Asian rice, which is 16 to 50 times higher than previously estimated for plants. Detailed analyses of these inversions show evidence of their effects on gene expression, recombination rate, and linkage disequilibrium. Our study uncovers the prevalence and scale of large inversions (≥100 bp) across the pan-genome of Asian rice and hints at their largely unexplored role in functional biology and crop performance.

Standing genetic variation and chromosome differences drove rapid ecotype formation in a major malaria mosquito

Species distributed across heterogeneous environments often evolve locally adapted ecotypes, but understanding of the genetic mechanisms involved in their formation and maintenance in the face of gene flow is incomplete. In Burkina Faso, the major African malaria mosquito Anopheles funestus comprises two strictly sympatric and morphologically indistinguishable yet karyotypically differentiated forms reported to differ in ecology and behavior. However, knowledge of the genetic basis and environmental determinants of An. funestus diversification was impeded by lack of modern genomic resources. Here, we applied deep whole-genome sequencing and analysis to test the hypothesis that these two forms are ecotypes differentially adapted to breeding in natural swamps versus irrigated rice fields. We demonstrate genome-wide differentiation despite extensive microsympatry, synchronicity, and ongoing hybridization. Demographic inference supports a split only ~1,300 y ago, closely following the massive expansion of domesticated African rice cultivation ~1,850 y ago. Regions of highest divergence, concentrated in chromosomal inversions, were under selection during lineage splitting, consistent with local adaptation. The origin of nearly all variations implicated in adaptation, including chromosomal inversions, substantially predates the ecotype split, suggesting that rapid adaptation was fueled mainly by standing genetic variation. Sharp inversion frequency differences likely facilitated adaptive divergence between ecotypes by suppressing recombination between opposing chromosomal orientations of the two ecotypes, while permitting free recombination within the structurally monomorphic rice ecotype. Our results align with growing evidence from diverse taxa that rapid ecological diversification can arise from evolutionarily old structural genetic variants that modify genetic recombination.

Genome assembly, structural variants, and genetic differentiation between lake whitefish young species pairs (Coregonus sp.) with long and short reads

Nascent pairs of ecologically differentiated species offer an opportunity to get a better glimpse at the genetic architecture of speciation. Of particular interest is our recent ability to consider a wider range of genomic variants, not only single-nucleotide polymorphisms (SNPs), thanks to long-read sequencing technology. We can now identify structural variants (SVs) such as insertions, deletions and other rearrangements, allowing further insights into the genetic architecture of speciation and how different types of variants are involved in species differentiation. Here, we investigated genomic patterns of differentiation between sympatric species pairs (Dwarf and Normal) belonging to the lake whitefish (Coregonus clupeaformis) species complex. We assembled the first reference genomes for both C. clupeaformis sp. Normal and C. clupeaformis sp. Dwarf, annotated the transposable elements and analysed the genomes in the light of related coregonid species. Next, we used a combination of long- and short-read sequencing to characterize SVs and genotype them at the population scale using genome-graph approaches, showing that SVs cover five times more of the genome than SNPs. We then integrated both SNPs and SVs to investigate the genetic architecture of species differentiation in two different lakes and highlighted an excess of shared outliers of differentiation. In particular, a large fraction of SVs differentiating the two species correspond to insertions or deletions of transposable elements (TEs), suggesting that TE accumulation may represent a key component of genetic divergence between the Dwarf and Normal species. Together, our results suggest that SVs may play an important role in speciation and that, by combining second- and third-generation sequencing, we now have the ability to integrate SVs into speciation genomics.

Population genomic evidence of selection on structural variants in a natural hybrid zone

Structural variants (SVs) can promote speciation by directly causing reproductive isolation or by suppressing recombination across large genomic regions. Whereas examples of each mechanism have been documented, systematic tests of the role of SVs in speciation are lacking. Here, we take advantage of long-read (Oxford nanopore) whole-genome sequencing and a hybrid zone between two Lycaeides butterfly taxa (L. melissa and Jackson Hole Lycaeides) to comprehensively evaluate genome-wide patterns of introgression for SVs and relate these patterns to hypotheses about speciation. We found >100,000 SVs segregating within or between the two hybridizing species. SVs and SNPs exhibited similar levels of genetic differentiation between species, with the exception of inversions, which were more differentiated. We detected credible variation in patterns of introgression among SV loci in the hybrid zone, with 562 of 1419 ancestry-informative SVs exhibiting genomic clines that deviated from null expectations based on genome-average ancestry. Overall, hybrids exhibited a directional shift towards Jackson Hole Lycaeides ancestry at SV loci, consistent with the hypothesis that these loci experienced more selection on average than SNP loci. Surprisingly, we found that deletions, rather than inversions, showed the highest skew towards excess ancestry from Jackson Hole Lycaeides. Excess Jackson Hole Lycaeides ancestry in hybrids was also especially pronounced for Z-linked SVs and inversions containing many genes. In conclusion, our results show that SVs are ubiquitous and suggest that SVs in general, but especially deletions, might disproportionately affect hybrid fitness and thus contribute to reproductive isolation.

Experimental introgression in Drosophila: Asymmetric postzygotic isolation associated with chromosomal inversions and an incompatibility locus on the X chromosome

Interspecific gene flow (introgression) is an important source of new genetic variation, but selection against it can reinforce reproductive barriers between interbreeding species. We used an experimental approach to trace the role of chromosomal inversions and incompatibility genes in preventing introgression between two partly sympatric Drosophila virilis group species, D. flavomontana and D. montana. We backcrossed F₁ hybrid females from a cross between D. flavomontana female and D. montana male with the males of the parental species for two generations and sequenced pools of parental strains and their reciprocal second generation backcross (BC₂ mon and BC₂ fla) females. Contrasting the observed amount of introgression (mean hybrid index, HI) in BC₂ female pools along the genome to simulations under different scenarios allowed us to identify chromosomal regions of restricted and increased introgression. We found no deviation from the HI expected under a neutral null model for any chromosome for the BC₂ mon pool, suggesting no evidence for genetic incompatibilities in backcrosses towards D. montana. In contrast, the BC₂ fla pool showed high variation in the observed HI between different chromosomes, and massive reduction of introgression on the X chromosome (large X-effect). This observation is compatible with reduced recombination combined with at least one dominant incompatibility locus residing within the X inversion(s). Overall, our study suggests that genetic incompatibilities arising within chromosomal inversions can play an important role in speciation.

Karyotype differentiation in the Nothobranchiusugandensis species group (Teleostei, Cyprinodontiformes), seasonal fishes from the east African inland plateau, in the context of phylogeny and biogeography

The karyotype differentiation of the twelve known members of the Nothobranchiusugandensis Wildekamp, 1994 species group is reviewed and the karyotype composition of seven of its species is described herein for the first time using a conventional cytogenetic protocol. Changes in the architecture of eukaryotic genomes often have a major impact on processes underlying reproductive isolation, adaptation and diversification. African annual killifishes of the genus Nothobranchius Peters, 1868 (Teleostei: Nothobranchiidae), which are adapted to an extreme environment of ephemeral wetland pools in African savannahs, feature extensive karyotype evolution in small, isolated populations and thus are suitable models for studying the interplay between karyotype change and species evolution. The present investigation reveals a highly conserved diploid chromosome number (2n = 36) but a variable number of chromosomal arms (46-64) among members of the N.ugandensis species group, implying a significant role of pericentric inversions and/or other types of centromeric shift in the karyotype evolution of the group. When superimposed onto a phylogenetic tree based on molecular analyses of two mitochondrial genes the cytogenetic characteristics did not show any correlation with the phylogenetic relationships within the lineage. While karyotypes of many other Nothobranchius spp. studied to date diversified mainly via chromosome fusions and fissions, the N.ugandensis species group maintains stable 2n and the karyotype differentiation seems to be constrained to intrachromosomal rearrangements. Possible reasons for this difference in the trajectory of karyotype differentiation are discussed. While genetic drift seems to be a major factor in the fixation of chromosome rearrangements in Nothobranchius, future studies are needed to assess the impact of predicted multiple inversions on the genome evolution and species diversification within the N.ugandensis species group.

Genome architecture and diverged selection shaping pattern of genomic differentiation in wild barley

Divergent selection of populations in contrasting environments leads to functional genomic divergence. However, the genomic architecture underlying heterogeneous genomic differentiation remains poorly understood. Here, we de novo assembled two high-quality wild barley (Hordeum spontaneum K. Koch) genomes and examined genomic differentiation and gene expression patterns under abiotic stress in two populations. These two populations had a shared ancestry and originated in close geographic proximity but experienced different selective pressures due to their contrasting micro-environments. We identified structural variants that may have played significant roles in affecting genes potentially associated with well-differentiated phenotypes such as flowering time and drought response between two wild barley genomes. Among them, a 29-bp insertion into the promoter region formed a cis-regulatory element in the HvWRKY45 gene, which may contribute to enhanced tolerance to drought. A single SNP mutation in the promoter region may influence HvCO5 expression and be putatively linked to local flowering time adaptation. We also revealed significant genomic differentiation between the two populations with ongoing gene flow. Our results indicate that SNPs and small SVs link to genetic differentiation at the gene level through local adaptation and are maintained through divergent selection. In contrast, large chromosome inversions may have shaped the heterogeneous pattern of genomic differentiation along the chromosomes by suppressing chromosome recombination and gene flow. Our research offers novel insights into the genomic basis underlying local adaptation and provides valuable resources for the genetic improvement of cultivated barley.

Widespread chromosomal rearrangements preceded genetic divergence in a monitor lizard, Varanus acanthurus (Varanidae)

Chromosomal rearrangements are often associated with local adaptation and speciation because they suppress recombination, and as a result, rearrangements have been implicated in disrupting gene flow. Although there is strong evidence to suggest that chromosome rearrangements are a factor in genetic isolation of divergent populations, the underlying mechanism remains elusive. Here, we applied an integrative cytogenetics and genomics approach testing whether chromosomal rearrangements are the initial process, or a consequence, of population divergence in the dwarf goanna, Varanus acanthurus. Specifically, we tested whether chromosome rearrangements are indicators of genetic barriers that can be used to identify divergent populations by looking at gene flow within and between populations with rearrangements. We found that gene flow was present between individuals with chromosome rearrangements within populations, but there was no gene flow between populations that had similar chromosome rearrangements. Moreover, we identified a correlation between reduced genetic variation in populations with a higher frequency of homozygous submetacentric individuals. These findings suggest that chromosomal rearrangements were widespread prior to divergence, and because we found populations with higher frequencies of submetacentric chromosomes were associated with lower genetic diversity, this could indicate that polymorphisms within populations are early indicators of genetic drift.

Sympatric or micro-allopatric speciation in a glacial lake? Genomic islands support neither

Apparent cases of sympatric speciation may actually be due to micro-allopatric or micro-parapatric speciation. One way to distinguish between these models is to examine the existence and nature of genomic islands of divergence, wherein divergent DNA segments are interspersed with low-divergence segments. Such islands should be rare or absent under micro-allopatric speciation but common in cases of speciation with gene flow. Sympatric divergence of endemic fishes is known from isolated saline, crater, postglacial, and ancient lakes. Two morphologically distinct cyprinid fishes, Gymnocypris eckloni scoliostomus (GS) and G. eckloni eckloni (GE), in a small glacial lake on the Qinghai-Tibet Plateau, Lake Sunmcuo, match the biogeographic criteria of sympatric speciation. In this study, we examined genome-wide variation in 46 individuals from these two groups. The divergence time between the GS and GE lineages was estimated to be 20-60 Kya. We identified 54 large genomic islands (≥100 kb) of speciation, which accounted for 89.4% of the total length of all genomic islands. These islands harboured divergent genes related to olfactory receptors and olfaction signals that may play important roles in food selection and assortative mating in fishes. Although the genomic islands clearly indicated speciation with gene flow and rejected micro-allopatric speciation, they were too large to support the hypothesis of sympatric speciation. Theoretical and recent empirical studies suggested that continual gene flow in sympatry should give rise to many small genomic islands (as small as a few kilobases in size). Thus, the observed pattern is consistent with the extensive evidence on parapatric speciation, in which adjacent habitats facilitate divergent selection but also permit gene flow during speciation. We suggest that many, if not most, of the reported cases of sympatric speciation are likely to be micro-parapatric speciation.

Construction of ultra-high-density genetic linkage map of a sorghum-sudangrass hybrid using whole genome resequencing

The sorghum-sudangrass hybrid is a vital annual gramineous herbage. Few reports exist on its ultra-high-density genetic map. In this study, we sought to create an ultra-high-density genetic linkage map for this hybrid to strengthen its functional genomics research and genetic breeding. We used 150 sorghum-sudangrass hybrid F2 individuals and their parents (scattered ear sorghum and red hull sudangrass) for high-throughput sequencing on the basis of whole genome resequencing. In total, 1,180.66 Gb of data were collected. After identification, filtration for integrity, and partial segregation, over 5,656 single nucleotide polymorphism markers of high quality were detected. An ultra-high-density genetic linkage map was constructed using these data. The markers covered approximately 2,192.84 cM of the map with average marker intervals of 0.39 cM. The length ranged from 115.39 cM to 264.04 cM for the 10 linkage groups. Currently, this represents the first genetic linkage map of this size, number of molecular markers, density, and coverage for sorghum-sudangrass hybrid. The findings of this study provide valuable genome-level information on species evolution and comparative genomics analysis and lay the foundation for further research on quantitative trait loci fine mapping and gene cloning and marker-assisted breeding of important traits in sorghum-sudangrass hybrids.

Analysis of the Whole-Genome Sequences from an Equus Parent-Offspring Trio Provides Insight into the Genomic Incompatibilities in the Hybrid Mule

Interspecific hybridization often shows negative effects on hybrids. However, only a few multicellular species, limited to a handful of plants and animals, have shown partial genetic mechanisms by which hybridization leads to low fitness in hybrids. Here, to explore the outcome of combining the two genomes of a horse and donkey, we analyzed the whole-genome sequences from an Equus parent-offspring trio using Illumina platforms. We generated 41.39× and 46.21× coverage sequences for the horse and mule, respectively. For the donkey, a 40.38× coverage sequence was generated and stored in our laboratory. Approximately 24.86 million alleles were discovered that varied from the reference genome. Single nucleotide polymorphisms were used as polymorphic markers for assigning alleles to their parental genomic inheritance. We identified 25,703 Mendelian inheritance error single nucleotide polymorphisms in the mule genome that were not inherited from the parents through Mendelian inheritance. A total of 555 de novo single nucleotide polymorphisms were also identified. The rate of de novo single nucleotide polymorphisms was 2.21 × 10^-7 in the mule from the Equus parent-offspring trio. This rate is obviously higher than the natural mutation rate for Equus, which is also consistent with the previous hypothesis that interracial crosses may have a high mutation rate. The genes associated with these single nucleotide polymorphisms are mainly involved in immune processes, DNA repair, and cancer processes. The results of the analysis of three genomes from an Equus parent-offspring trio improved our knowledge of the consequences of the integration of parental genomes in mules.

Getting higher on rugged landscapes: Inversion mutations open access to fitter adaptive peaks in NK fitness landscapes

Molecular evolution is often conceptualised as adaptive walks on rugged fitness landscapes, driven by mutations and constrained by incremental fitness selection. It is well known that epistasis shapes the ruggedness of the landscape’s surface, outlining their topography (with high-fitness peaks separated by valleys of lower fitness genotypes). However, within the strong selection weak mutation (SSWM) limit, once an adaptive walk reaches a local peak, natural selection restricts passage through downstream paths and hampers any possibility of reaching higher fitness values. Here, in addition to the widely used point mutations, we introduce a minimal model of sequence inversions to simulate adaptive walks. We use the well known NK model to instantiate rugged landscapes. We show that adaptive walks can reach higher fitness values through inversion mutations, which, compared to point mutations, allows the evolutionary process to escape local fitness peaks. To elucidate the effects of this chromosomal rearrangement, we use a graph-theoretical representation of accessible mutants and show how new evolutionary paths are uncovered. The present model suggests a simple mechanistic rationale to analyse escapes from local fitness peaks in molecular evolution driven by (intragenic) structural inversions and reveals some consequences of the limits of point mutations for simulations of molecular evolution.

Ninety years ago, Wright translated Darwin’s core idea of survival of the fittest into rugged landscapes—a highly influential metaphor—with peaks representing high values of fitness separated by valleys of lower fitness. In this picture, once a population has reached a local peak, the adaptive dynamics may stall as further adaptation requires crossing a valley. At the DNA level, adaptation is often modelled as a space of genotypes that is explored through point mutations. Therefore, once a local peak is reached, any genotype fitter than that of the peak will be away from the neighbourhood of genotypes accessible through point mutations. Here we present a simple computational model for inversion mutations, one of the most frequent structural variations, and show that adaptive processes in rugged landscapes can escape from local peaks through intragenic inversion mutations. This new escape mechanism reveals the innovative role of inversions at the DNA level and provides a step towards more realistic models of adaptive dynamics, beyond the dominance of point mutations in theories of molecular evolution.

Genomic architecture and functional effects of potential human inversion supergenes

Supergenes are involved in adaptation in multiple organisms, but they are little known in humans. Genomic inversions are the most common mechanism of supergene generation and maintenance. Here, we review the information about two large inversions that are the best examples of potential human supergenes. In addition, we do an integrative analysis of the newest data to understand better their functional effects and underlying genetic changes. We have found that the highly divergent haplotypes of the 17q21.31 inversion of approximately 1.5 Mb have multiple phenotypic associations, with consistent effects in brain-related traits, red and white blood cells, lung function, male and female characteristics and disease risk. By combining gene expression and nucleotide variation data, we also analysed the molecular differences between haplotypes, including gene duplications, amino acid substitutions and regulatory changes, and identify CRHR1, KANLS1 and MAPT as good candidates to be responsible for these phenotypes. The situation is more complex for the 8p23.1 inversion, where there is no clear genetic differentiation. However, the inversion is associated with several related phenotypes and gene expression differences that could be linked to haplotypes specific of one orientation. Our work, therefore, contributes to the characterization of both exceptional variants and illustrates the important role of inversions. This article is part of the theme issue 'Genomic architecture of supergenes: causes and evolutionary consequences'.

Inversions and parallel evolution

Local adaptation leads to differences between populations within a species. In many systems, similar environmental contrasts occur repeatedly, sometimes driving parallel phenotypic evolution. Understanding the genomic basis of local adaptation and parallel evolution is a major goal of evolutionary genomics. It is now known that by preventing the break-up of favourable combinations of alleles across multiple loci, genetic architectures that reduce recombination, like chromosomal inversions, can make an important contribution to local adaptation. However, little is known about whether inversions also contribute disproportionately to parallel evolution. Our aim here is to highlight this knowledge gap, to showcase existing studies, and to illustrate the differences between genomic architectures with and without inversions using simple models. We predict that by generating stronger effective selection, inversions can sometimes speed up the parallel adaptive process or enable parallel adaptation where it would be impossible otherwise, but this is highly dependent on the spatial setting. We highlight that further empirical work is needed, in particular to cover a broader taxonomic range and to understand the relative importance of inversions compared to genomic regions without inversions.

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