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

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.

Temporal and habitat adaptations in Drosophila subobscura populations: changes in chromosomal inversions

In insects, chromosomal inversion polymorphism has been related with different adaptations, including global warming. Regarding this environmental change, Drosophila subobscura stands out as a useful model species due to its rich inversion polymorphism covering the whole karyotype. The main aims of this research were to analyze the differentiation of this polymorphism in Jastrebac Mt. (Serbia) depending on the different habitats (beech and oak forests) and over time. These latter changes were studied in relation to climatic variables (mean, minimum and maximum temperatures, humidity and rainfall). Microdifferentiation was observed between beech and oak forests, mainly for the A and O chromosomes. However, the changes over time turned out to be larger than those due to habitat. In Jastrebac Mt., temperatures increased over time, with this increase being significant for mean and minimum one. The Multidimensional Scaling analysis showed a relation between chromosomal inversions and temperatures (mainly minimum) in Jastrebac Mt. and other Serbian populations of D. subobscura. In beech forest of Jastrebac Mt., the Chromosomal Thermal Index increased over time from 1990 to 1994, but showing a possible stabilization in 2023. This result was observed in other studied Serbian populations. Although those are preliminary results, it might hypothesize that there may be a threshold for the action of natural selection, increasing 'warm' adaptive inversions and decreasing 'cold' ones. The possible reasons for this hypothesis are also discussed.

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.

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.

Chromosomal structural variation loci HSS1 and HSS6 lead to hybrid sterility in rice

KEY MESSAGE

Chromosomal structural variation leads to rice hybrid sterility, and this study has established a new inheritance model and offers a novel insight into hybrid sterility in rice. The utilization of heterosis among different species significantly enhances rice yield. However, the direct application of this advantage is hindered by hybrid sterility (HS). Here we identify a stable semi-sterile inbred line (SSIL) in the hybridization between rice species O. sativa indica and wild rice O. longistaminata. Both pollen and spikelet fertility in heterozygous SSIL plants are semi-sterile, and the homozygous plants show normal fertility. Interestingly, unlike previously reported hybrid sterility loci, SSIL does not induce segregation distortion in progeny. The genotypes of chromosome 1 (Chr.1) exhibit tight pseudo-linkage on Chr.6 in SSIL population. HSS1 and HSS6 are named as the two loci of Chr.1 and Chr.6. Cytological observations revealed abnormalities in male and female gametes during meiosis, ultimately resulting in semi-sterile pollen and spikelet. HSS1 was narrowed down to a 190.6-kb interval and HSS6 to a 1391.6-kb region using an SSIL population of 10,393 plants. Through high-throughput sequencing and observation of chromosomal behavior during meiosis, a reciprocal translocation between the short arm of Chr.1 and the end of Chr.6 in the translocated RD23-type (RT-type) gamete in SSIL was discovered. The presence of RT-type gametes in SSIL is the direct causative factor for semi-sterility in both male and female gametes. This translocation led to abnormal SSIL-RT/L synapses and the formation of tetravalent ring structures in chromosomes during pachynema. Our findings have uncovered a pair of reciprocal chromosome translocations that control reproductive isolation in rice. These insights offer valuable guidance for optimizing hybrid breeding applications, ultimately enhancing the benefits of heterosis.

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.

Advances in Understanding the Karyotype Evolution of Tetrapulmonata and Two Other Arachnid Taxa, Ricinulei and Solifugae

BACKGROUND/OBJECTIVES

Arachnids are a megadiverse arthropod group. The present study investigated the chromosomes of pedipalpid tetrapulmonates (orders Amblypygi, Thelyphonida, Schizomida) and two arachnid orders of uncertain phylogenetic placement, Ricinulei and Solifugae, to reconstruct their karyotype evolution. Except for amblypygids, the cytogenetics of these arachnid orders was almost unknown prior to the present study.

METHODS

Chromosomes were investigated using methods of standard (Giemsa-stained preparations, banding techniques) and molecular cytogenetics (fluorescence in situ hybridization, comparative genomic hybridization).

RESULTS AND CONCLUSIONS

New data for 38 species, combined with previously published data, suggest that ancestral arachnids possessed low to moderate 2n (22-40), monocentric chromosomes, one nucleolus organizer region (NOR), low levels of heterochromatin and recombinations, and no or homomorphic sex chromosomes. Karyotypes of Pedipalpi and Solifugae diversified via centric fusions, pericentric inversions, and changes in the pattern of NORs and, in solifuges, also through tandem fusions. Some solifuges display an enormous amount of constitutive heterochromatin and high NOR number. It is hypothesized that the common ancestor of amblypygids, thelyphonids, and spiders exhibited a homomorphic XY system, and that telomeric heterochromatin and NORs were involved in the evolution of amblypygid sex chromosomes. The new findings support the Cephalosomata clade (acariforms, palpigrades, and solifuges). Hypotheses concerning the origin of acariform holocentric chromosomes are presented. Unlike current phylogenetic hypotheses, the results suggest a sister relationship between Schizomida and a clade comprising other tetrapulmonates as well as a polyploidization in the common ancestor of the clade comprising Araneae, Amblypygi, and Thelyphonida.

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.

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.

Short-term evolutionary implications of an introgressed size-determining supergene in a vulnerable population

The Thorny Skate (Amblyraja radiata) is a vulnerable species displaying a discrete size-polymorphism in the northwest Atlantic Ocean (NWA). We conducted whole genome sequencing of samples collected across its range. Genetic diversity was similar at all sampled sites, but we discovered a ~ 31 megabase bi-allelic supergene associated with the size polymorphism, with the larger size allele having introgressed in the last ~160,000 years B.P. While both Gulf of Maine (GoM) and Canadian (CAN) populations exhibit the size polymorphism, we detected a significant deficit of heterozygotes at the supergene and longer stretches of homozygosity in GoM population. This suggests inbreeding driven by assortative mating for size in GoM but not in CAN. Coalescent-based demographic modelling reveals strong migration between regions maintaining genetic variability in the recombining genome, preventing speciation between morphs. This study highlights short-term context-dependent evolutionary consequences of a size-determining supergene providing new insights for the management of vulnerable species.

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.

Divergent Population Structure in Five Common Rockfish Species of Puget Sound, WA Suggests the Need for Species-Specific Management

Quantifying connectivity between endangered or threatened marine populations is critical information for management and conservation, especially where abundance and productivity differ among such populations. Spatial patterns of such connectivity depend not only on extrinsic factors such as oceanography and bathymetry but also on intrinsic species-specific factors such as life history, demography and the location of glacial refugia. Nevertheless, population structure is often inferred from related or ecologically similar species. For example, the population structure in most rockfish species (Sebastes spp.) in the Salish Sea and the US West Coast is currently inferred from genetic data of three species that are known to hybridise in Puget Sound. Here, we determined the population structure and connectivity in five Puget Sound Rockfish species (Black [Sebastes melanops], Yellowtail [S. flavidus], Redstripe [S. proriger], Greenstriped [S. elongatus], and Puget Sound Rockfish [S. emphaeus]) from over 12,000 restriction-site associated DNA sequencing (RADseq) loci. We found species-specific patterns of genetic differentiation, attributable to both extrinsic and intrinsic factors. Specifically, Black and Puget Sound rockfishes showed no genetic differentiation; Yellowtail and Greenstriped rockfishes were structured according to known geographic barriers; and Redstripe Rockfish revealed evidence for temporal genetic differentiation, suggesting irregular recruitment influences population structure. Only Yellowtail Rockfish followed the federal DPS boundaries generally assumed for rockfish, further emphasizing the importance of species-specific management for the effective recovery and management of these rockfish populations and of marine species 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.

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.

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.

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.

Genomics-Informed Range Predictions Under Global Warming Reveal Reduced Adaptive Diversity Whilst Buffering Range Shifts for a Marine Snail

Understanding the genetic basis of local adaptation in thermal performance is useful for predicting species distribution shifts under anthropogenic climate change. Many species are distributed across multiple biogeographic regions, and the uniquely adapted populations in each region may respond to future ocean warming with distinct distribution changes. In the present study, we investigated phylogeographic patterns, thermal sensitivity, and genetic differentiation in the intertidal snail Littorina brevicula along China's coast. Whole-genome sequencing results based on a newly assembled chromosome-level genome revealed two genetic lineages, with a north-south divergence that is linked to the thermal environment. Within each lineage, individuals could be further subdivided into genetic subgroups that differ at key genomic loci underpinning differences in upper heat tolerance. Heat stress drives adaptive divergence across multiple levels of organization, from the individual to the biogeographic level. Taking into account genetic diversity associated with variation in heat tolerance, a physiological species distribution model (pSDM) was applied to predict the distributions of the different genetic subgroups in response to climate change. Both northern and southern lineages were predicted to experience declines in habitat suitability under a 4°C future warming scenario, and that a genotypic subset of snails from the southern lineage may even be driven to extinction. These findings illustrate that even when a species' range is maintained, it can nonetheless experience a significant decrease in adaptive diversity as a result of climate change. The integrated approach presented here, which considered both physiological and adaptive genetic variation at the level of individuals within a biogeographical context, provided new insights into how marine species can respond to global warming.

Evolution and Plasticity of Gene Expression Under Progressive Warming in Drosophila subobscura

Understanding the molecular mechanisms of thermal adaptation is crucial to predict the impacts of global warming. However, there is still a lack of research on the effects of rising temperatures over time and of studies involving different populations from the same species. The present study focuses on these two aspects, which are of great importance in understanding how organisms cope and adapt to ongoing changes in their environment. This study investigates the impact of global warming on the gene expression patterns of Drosophila subobscura populations from two different latitudinal locations after 23 generations of evolution. Our results indicate that evolutionary changes depend on the genetic background of the populations, with different starting points for thermal evolution, and that high-latitude populations show more pronounced evolutionary changes, with some evidence of convergence towards low-latitude populations. We found an interplay between plasticity and selection, with the high-latitude population showing fewer initial plastic genes and lower levels of adaptive plasticity, but a greater magnitude of change in both plastic and selective responses during evolution under warming conditions compared with its low-latitude counterpart. A substantial proportion of the transcriptome was observed to be evolving, despite the lack of observable response at higher-order phenotypic traits. The interplay between plasticity and selection may prove to be an essential component in shaping species' evolutionary responses to climate change. Furthermore, the value of conducting studies on multiple populations of the same species is emphasised, given the identification of differences between populations with different backgrounds in several contexts.

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.

Transposable elements in Drosophila montana from harsh cold environments

BACKGROUND

Substantial discoveries during the past century have revealed that transposable elements (TEs) can play a crucial role in genome evolution by affecting gene expression and inducing genetic rearrangements, among other molecular and structural effects. Yet, our knowledge on the role of TEs in adaptation to extreme climates is still at its infancy. The availability of long-read sequencing has opened up the possibility to identify and study potential functional effects of TEs with higher precision. In this work, we used Drosophila montana as a model for cold-adapted organisms to study the association between TEs and adaptation to harsh climates.

RESULTS

Using the PacBio long-read sequencing technique, we de novo identified and manually curated TE sequences in five Drosophila montana genomes from eco-geographically distinct populations. We identified 489 new TE consensus sequences which represented 92% of the total TE consensus in D. montana. Overall, 11-13% of the D. montana genome is occupied by TEs, which as expected are non-randomly distributed across the genome. We identified five potentially active TE families, most of them from the retrotransposon class of TEs. Additionally, we found TEs present in the five analyzed genomes that were located nearby previously identified cold tolerant genes. Some of these TEs contain promoter elements and transcription binding sites. Finally, we detected TEs nearby fixed and polymorphic inversion breakpoints.

CONCLUSIONS

Our research revealed a significant number of newly identified TE consensus sequences in the genome of D. montana, suggesting that non-model species should be studied to get a comprehensive view of the TE repertoire in Drosophila species and beyond. Genome annotations with the new D. montana library allowed us to identify TEs located nearby cold tolerant genes, and present at high population frequencies, that contain regulatory regions and are thus good candidates to play a role in D. montana cold stress response. Finally, our annotations also allow us to identify for the first time TEs present in the breakpoints of three D. montana inversions.

Identification of Genomic Regions Associated with Differences in Flowering Time and Inflorescence Architecture between Melastoma candidum and M. normale

Understanding the genetic basis of species differences in flowering time and inflorescence architecture can shed light on speciation and molecular breeding. Melastoma shows rapid speciation, with about 100 species formed in the past few million years, and, meanwhile, possesses high ornamental values. Two largely sympatric and closely related species of this genus, M. candidum and M. normale, differ markedly in flowering time and flower number per inflorescence. Here, we constructed an F2 population between M. candidum and M. normale, and used extreme bulks for flowering time and flower number per inflorescence in this population to identify genomic regions underlying the two traits. We found high differentiation on nearly the whole chromosome 7 plus a few regions on other chromosomes between the two extreme bulks for flowering time. Large chromosomal inversions on chromosome 7 between the two species, which contain flowering-related genes, can explain recombinational suppression on the chromosome. We identified 1872 genes with one or more highly differentiated SNPs between the two bulks for flowering time, including CSTF77, FY, SPA3, CDF3, AGL8, AGL15, FHY1, COL9, CIB1, FKF1 and FAR1, known to be related to flowering. We also identified 680 genes with one or more highly differentiated SNPs between the two bulks for flower number per inflorescence, including PNF, FIL and LAS, knows to play important roles in inflorescence development. These large inversions on chromosome 7 prevent us from narrowing down the genomic region(s) associated with flowering time differences between the two species. Flower number per inflorescence in Melastoma appears to be controlled by multiple genes, without any gene of major effect. Our study indicates that large chromosomal inversions can hamper the identification of the genetic basis of important traits, and the inflorescence architecture of Melastoma species may have a complex genetic basis.

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, facilitating 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 acute heat tolerance, which phenotypically trades off with tolerance to prolonged heat exposure. A simple evolutionary model based on our data shows that the estimated maximum rate of evolutionary adaptation in mosquito heat tolerance typically exceeds that of projected climate warming under idealized conditions. Our findings indicate that natural mosquito populations may have the potential to track projected warming via genetic adaptation. Prior climate-based projections may thus underestimate the range of mosquito and mosquito-borne disease distributions under future climate conditions.

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.

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.

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.

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.

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.

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.

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.