Inversions are bigger on the X chromosome

The evolution of suppressed recombination between sex chromosomes and the lengths of evolutionary strata

The idea that sex-differences in selection drive the evolution of suppressed recombination between sex chromosomes is well-developed in population genetics. Yet, despite a now classic body of theory, empirical evidence that sexually antagonistic (SA) selection drives the evolution of recombination arrest remains equivocal and alternative hypotheses underdeveloped. Here, we investigate whether the length of "evolutionary strata" formed by chromosomal inversions (or other large-effect recombination modifiers) expanding the nonrecombining sex-linked region (SLR) on sex chromosomes can be informative of how selection influenced their fixation. We develop population genetic models to show how the length of an SLR-expanding inversion and the presence of partially recessive deleterious mutational variation affect the fixation probability of three different classes of inversions: (i) intrinsically neutral, (ii) directly beneficial (i.e., due to breakpoint or positional effects), and (iii) those capturing SA loci. Our models indicate that inversions capturing an SA locus initially in linkage disequilibrium with the ancestral SLR exhibit a strong fixation bias toward small inversions, while neutral, beneficial, and inversions capturing a genetically unlinked SA locus tend to favor larger inversions and exhibit similar distributions of fixed inversion lengths. The footprint of evolutionary stratum size left behind by different selection regimes is strongly influenced by parameters affecting the deleterious mutation load, the physical position of the ancestral SLR, and the distribution of new inversion lengths.

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

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

Model-based inference of a dual role for HOPS in regulating guard cell vacuole fusion

Guard cell movements depend, in part, on the remodelling of vacuoles from a highly fragmented state to a fused morphology during stomata opening. Indeed, full opening of plant stomata requires vacuole fusion to occur. Fusion of vacuole membranes is a highly conserved process in eukaryotes, with key roles played by two multi-subunit complexes: HOPS (homotypic fusion and vacuolar protein sorting) and SNARE (soluble NSF attachment protein receptor). HOPS is a vacuole tethering factor that is thought to chaperone SNAREs from apposing vacuole membranes into a fusion-competent complex capable of rearranging membranes. In plants, recruitment of HOPS subunits to the tonoplast has been shown to require the presence of the phosphoinositide phosphatidylinositol 3-phosphate. However, chemically depleting this lipid induces vacuole fusion. To resolve this counter-intuitive observation regarding the role of HOPS in regulating plant vacuole morphology, we defined a quantitative model of vacuole fusion dynamics and used it to generate testable predictions about HOPS-SNARE interactions. We derived our model by using simulation-based inference to integrate prior knowledge about molecular interactions with limited, qualitative observations of emergent vacuole phenotypes. By constraining the model parameters to yield the emergent outcomes observed for stoma opening-as induced by two distinct chemical treatments-we predicted a dual role for HOPS and identified a stalled form of the SNARE complex that differs from phenomena reported in yeast. We predict that HOPS has contradictory actions at different points in the fusion signalling pathway, promoting the formation of SNARE complexes, but limiting their activity.

Two Nested Inversions in the X Chromosome Differentiate the Dominant Malaria Vectors in Europe, Anopheles atroparvus and Anopheles messeae

The Maculipennis subgroup of malaria mosquitoes includes both dominant malaria vectors and non-vectors in Eurasia. Understanding the genetic factors, particularly chromosomal inversions, that differentiate Anopheles species can provide valuable insights for vector control strategies. Although autosomal inversions between the species in this subgroup have been characterized based on the chromosomal banding patterns, the number and positions of rearrangements in the X chromosome remain unclear due to the divergent banding patterns. Here, we identified two large X chromosomal inversions, approximately 13 Mb and 10 Mb in size, using fluorescence in situ hybridization. The inversion breakpoint regions were mapped by hybridizing 53 gene markers with polytene chromosomes of An. messeae. The DNA probes were designed based on gene sequences from the annotated An. atroparvus genome. The two nested inversions resulted in five syntenic blocks. Only two small syntenic blocks, which encompass 181 annotated genes in the An. atroparvus genome, changed their position and orientation in the X chromosome. The analysis of the An. atroparvus genome revealed an enrichment of gene ontology terms associated with immune system and mating behavior in the rearranged syntenic blocks. Additionally, the enrichment of DNA transposons was found in sequences homologous to three of the four breakpoint regions. This study demonstrates the successful application of the physical genome mapping approach to identify rearrangements that differentiate species in insects with polytene chromosomes.

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.

Model-based inference of a plant-specific dual role for HOPS in regulating guard cell vacuole fusion

Stomata are the pores on a leaf surface that regulate gas exchange. Each stoma consists of two guard cells whose movements regulate pore opening and thereby control CO2 fixation and water loss. Guard cell movements depend in part on the remodeling of vacuoles, which have been observed to change from a highly fragmented state to a fused morphology during stomata opening. This change in morphology requires a membrane fusion mechanism that responds rapidly to environmental signals, allowing plants to respond to diurnal and stress cues. With guard cell vacuoles being both large and responsive to external signals, stomata represent a unique system in which to delineate mechanisms of membrane fusion. Fusion of vacuole membranes is a highly conserved process in eukaryotes, with key roles played by two multi-subunit complexes: HOPS (homotypic fusion and vacuolar protein sorting) and SNARE (soluble NSF attachment protein receptor). HOPS is a vacuole tethering factor that is thought to chaperone SNAREs from apposing vacuole membranes into a fusion-competent complex capable of rearranging membranes. To resolve a counter-intuitive observation regarding the role of HOPS in regulating plant vacuole morphology, we derived a quantitative model of vacuole fusion dynamics and used it to generate testable predictions about HOPS-SNARE interactions. We derived our model by applying simulation-based inference to integrate prior knowledge about molecular interactions with limited, qualitative observations of emergent vacuole phenotypes. By constraining the model parameters to yield the emergent outcomes observed for stoma opening - as induced by two distinct chemical treatments - we predicted a dual role for HOPS and identified a stalled form of the SNARE complex that differs from phenomena reported in yeast. We predict that HOPS has contradictory actions at different points in the fusion signaling pathway, promoting the formation of SNARE complexes, but limiting their activity.

Single-haplotype comparative genomics provides insights into lineage-specific structural variation during cat evolution

The role of structurally dynamic genomic regions in speciation is poorly understood due to challenges inherent in diploid genome assembly. Here we reconstructed the evolutionary dynamics of structural variation in five cat species by phasing the genomes of three interspecies F1 hybrids to generate near-gapless single-haplotype assemblies. We discerned that cat genomes have a paucity of segmental duplications relative to great apes, explaining their remarkable karyotypic stability. X chromosomes were hotspots of structural variation, including enrichment with inversions in a large recombination desert with characteristics of a supergene. The X-linked macrosatellite DXZ4 evolves more rapidly than 99.5% of the genome clarifying its role in felid hybrid incompatibility. Resolved sensory gene repertoires revealed functional copy number changes associated with ecomorphological adaptations, sociality and domestication. This study highlights the value of gapless genomes to reveal structural mechanisms underpinning karyotypic evolution, reproductive isolation and ecological niche adaptation.

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.

The evolution of suppressed recombination between sex chromosomes and the lengths of evolutionary strata

The idea that sex differences in selection drive the evolution of suppressed recombination between sex chromosomes is well developed in population genetics. Yet, despite a now classic body of theory, empirical evidence that sexually antagonistic selection drives the evolution of recombination arrest remains equivocal and alternative hypotheses underdeveloped. Here, we investigate whether the length of "evolutionary strata" formed by chromosomal inversions (or other large-effect recombination modifiers) expanding the non-recombining sex-linked region (SLR) on sex chromosomes can be informative of how selection influenced their fixation. We develop population genetic models to show how the length of an SLR-expanding inversion, and the presence of partially recessive deleterious mutational variation, affect the fixation probability of three different classes of inversions: (1) intrinsically neutral, (2) directly beneficial (i.e., due to breakpoint or positional effects), and (3) those capturing sexually antagonistic (SA) loci. Our models indicate that neutral inversions, and those capturing an SA locus in linkage disequilibrium with the ancestral SLR, will exhibit a strong fixation bias toward small inversions; while unconditionally beneficial inversions, and those capturing a genetically unlinked SA locus, will favor fixation of larger inversions. The footprint of evolutionary stratum size left behind by different selection regimes is strongly influenced by parameters affecting the deleterious mutation load, the physical position of the ancestral SLR, and the distribution of new inversion lengths.

What inversion lengths can tell us about their evolution

Genome structural differences, such as inversions, are extremely common between species and within populations. Theoretical models of how and why such inversions evolve have affirmed that they are able to evolve under both adaptive and non-adaptive scenarios (reviewed in Kirkpatrick 2010). What has remained difficult, however, is distinguishing these scenarios from each other. In this issue of Molecular Ecology, Connallon and Olito (2022) present a model that examines how adaptive and non-adaptive scenarios lead to different distributions of inversion sizes. The authors present several important predictions including an expectation that larger inversions should evolve under local adaptation scenarios and much smaller inversions should evolve when they are either underdominant or directly beneficial. Finally, the authors ask how the presence of deleterious mutations within populations affects the probability of fixing inversions of different types. The paper is therefore an important step in synthesizing decades of inversion theory.

Solving the puzzle of Fe homeostasis by integrating molecular, mathematical, and societal models

To ensure optimal utilization and bioavailability, iron uptake, transport, subcellular localization, and assimilation are tightly regulated in plants. Herein, we examine recent advances in our understanding of cellular responses to Fe deficiency. We then use intracellular mechanisms of Fe homeostasis to discuss how formalizing cell biology knowledge via a mathematical model can advance discovery even when quantitative data is limited. Using simulation-based inference to identify plausible systems mechanisms that conform to known emergent phenotypes can yield novel, testable hypotheses to guide targeted experiments. However, this approach relies on the accurate encoding of domain-expert knowledge in exploratory mathematical models. We argue that this would be facilitated by fostering more "systems thinking" life scientists and that diversifying your research team may be a practical path to achieve that goal.

Natural selection and the distribution of chromosomal inversion lengths

Chromosomal inversions contribute substantially to genome evolution, yet the processes governing their evolutionary dynamics remain poorly understood. Theory suggests that a readily measurable property of inversions-their length-can potentially affect their evolutionary fates. Emerging data on the lengths of polymorphic and fixed inversions may therefore provide clues to the evolutionary processes promoting inversion establishment. However, formal predictions for the distribution of inversion lengths remain incomplete, making empirical patterns difficult to interpret. We model the relation between inversion length and establishment probability for four inversion types: (1) neutral, (2) underdominant, (3) directly beneficial, and (4) indirectly beneficial, with selection favouring the latter because they capture locally adapted alleles at migration-selection balance and suppress recombination between them. We also consider how deleterious mutations affect the lengths of established inversions. We show that length distributions of common polymorphic and fixed inversions systematically differ among inversion types. Small rearrangements contribute the most to genome evolution under neutral and underdominant scenarios of selection, with the lengths of neutral inversion substitutions increasing, and those of underdominant substitutions decreasing, with effective population size. Among directly beneficial inversions, small rearrangements are preferentially fixed, whereas intermediate-to-large inversions are maintained as balanced polymorphisms via associative overdominance. Finally, inversions established under the local adaptation scenario are predominantly intermediate-to-large. Such inversions remain polymorphic or approach fixation within the local populations where they are favoured. Our models clarify how inversion length distributions relate to processes of inversion establishment, providing a platform for testing how natural selection shapes the evolution of genome structure.