@article{14742,
  abstract     = {Chromosomal rearrangements (CRs) have been known since almost the beginning of genetics.
While an important role for CRs in speciation has been suggested, evidence primarily stems
from theoretical and empirical studies focusing on the microevolutionary level (i.e., on taxon
pairs where speciation is often incomplete). Although the role of CRs in eukaryotic speciation at
a macroevolutionary level has been supported by associations between species diversity and
rates of evolution of CRs across phylogenies, these findings are limited to a restricted range of
CRs and taxa. Now that more broadly applicable and precise CR detection approaches have
become available, we address the challenges in filling some of the conceptual and empirical
gaps between micro- and macroevolutionary studies on the role of CRs in speciation. We
synthesize what is known about the macroevolutionary impact of CRs and suggest new research avenues to overcome the pitfalls of previous studies to gain a more comprehensive understanding of the evolutionary significance of CRs in speciation across the tree of life.},
  author       = {Lucek, Kay and Giménez, Mabel D. and Joron, Mathieu and Rafajlović, Marina and Searle, Jeremy B. and Walden, Nora and Westram, Anja M and Faria, Rui},
  issn         = {1943-0264},
  journal      = {Cold Spring Harbor Perspectives in Biology},
  keywords     = {General Biochemistry, Genetics and Molecular Biology},
  number       = {11},
  publisher    = {Cold Spring Harbor Laboratory},
  title        = {{The impact of chromosomal rearrangements in speciation: From micro- to macroevolution}},
  doi          = {10.1101/cshperspect.a041447},
  volume       = {15},
  year         = {2023},
}

@article{11479,
  abstract     = {Understanding population divergence that eventually leads to speciation is essential for evolutionary biology. High species diversity in the sea was regarded as a paradox when strict allopatry was considered necessary for most speciation events because geographical barriers seemed largely absent in the sea, and many marine species have high dispersal capacities. Combining genome-wide data with demographic modelling to infer the demographic history of divergence has introduced new ways to address this classical issue. These models assume an ancestral population that splits into two subpopulations diverging according to different scenarios that allow tests for periods of gene flow. Models can also test for heterogeneities in population sizes and migration rates along the genome to account, respectively, for background selection and selection against introgressed ancestry. To investigate how barriers to gene flow arise in the sea, we compiled studies modelling the demographic history of divergence in marine organisms and extracted preferred demographic scenarios together with estimates of demographic parameters. These studies show that geographical barriers to gene flow do exist in the sea but that divergence can also occur without strict isolation. Heterogeneity of gene flow was detected in most population pairs suggesting the predominance of semipermeable barriers during divergence. We found a weak positive relationship between the fraction of the genome experiencing reduced gene flow and levels of genome-wide differentiation. Furthermore, we found that the upper bound of the ‘grey zone of speciation’ for our dataset extended beyond that found before, implying that gene flow between diverging taxa is possible at higher levels of divergence than previously thought. Finally, we list recommendations for further strengthening the use of demographic modelling in speciation research. These include a more balanced representation of taxa, more consistent and comprehensive modelling, clear reporting of results and simulation studies to rule out nonbiological explanations for general results.},
  author       = {De Jode, Aurélien and Le Moan, Alan and Johannesson, Kerstin and Faria, Rui and Stankowski, Sean and Westram, Anja M and Butlin, Roger K. and Rafajlović, Marina and Fraisse, Christelle},
  issn         = {1752-4571},
  journal      = {Evolutionary Applications},
  number       = {2},
  pages        = {542--559},
  publisher    = {Wiley},
  title        = {{Ten years of demographic modelling of divergence and speciation in the sea}},
  doi          = {10.1111/eva.13428},
  volume       = {16},
  year         = {2023},
}

@article{12157,
  abstract     = {Polygenic adaptation is thought to be ubiquitous, yet remains poorly understood. Here, we model this process analytically, in the plausible setting of a highly polygenic, quantitative trait that experiences a sudden shift in the fitness optimum. We show how the mean phenotype changes over time, depending on the effect sizes of loci that contribute to variance in the trait, and characterize the allele dynamics at these loci. Notably, we describe the two phases of the allele dynamics: The first is a rapid phase, in which directional selection introduces small frequency differences between alleles whose effects are aligned with or opposed to the shift, ultimately leading to small differences in their probability of fixation during a second, longer phase, governed by stabilizing selection. As we discuss, key results should hold in more general settings and have important implications for efforts to identify the genetic basis of adaptation in humans and other species.},
  author       = {Hayward, Laura and Sella, Guy},
  issn         = {2050-084X},
  journal      = {eLife},
  keywords     = {General Immunology and Microbiology, General Biochemistry, Genetics and Molecular Biology, General Medicine, General Neuroscience},
  publisher    = {eLife Sciences Publications},
  title        = {{Polygenic adaptation after a sudden change in environment}},
  doi          = {10.7554/elife.66697},
  volume       = {11},
  year         = {2022},
}

@article{12166,
  abstract     = {Kerstin Johannesson is a marine ecologist and evolutionary biologist based at the Tjärnö Marine Laboratory of the University of Gothenburg, which is situated in the beautiful Kosterhavet National Park on the Swedish west coast. Her work, using marine periwinkles (especially Littorina saxatilis and L. fabalis) as main model systems, has made a remarkable contribution to marine evolutionary biology and our understanding of local adaptation and its genetic underpinnings.},
  author       = {Westram, Anja M and Butlin, Roger},
  issn         = {1365-294X},
  journal      = {Molecular Ecology},
  keywords     = {Genetics, Ecology, Evolution, Behavior and Systematics},
  number       = {1},
  pages        = {26--29},
  publisher    = {Wiley},
  title        = {{Professor Kerstin Johannesson–winner of the 2022 Molecular Ecology Prize}},
  doi          = {10.1111/mec.16779},
  volume       = {32},
  year         = {2022},
}

@article{12234,
  abstract     = {Hybrid speciation—the origin of new species resulting from the hybridization of genetically divergent lineages—was once considered rare, but genomic data suggest that it may occur more often than once thought. In this study, Noguerales and Ortego found genomic evidence supporting the hybrid origin of a grasshopper that is able to exploit a broader range of host plants than either of its putative parents.},
  author       = {Stankowski, Sean},
  issn         = {1558-5646},
  journal      = {Evolution},
  keywords     = {General Agricultural and Biological Sciences, Genetics, Ecology, Evolution, Behavior and Systematics},
  number       = {11},
  pages        = {2784--2785},
  publisher    = {Wiley},
  title        = {{Digest: On the origin of a possible hybrid species}},
  doi          = {10.1111/evo.14632},
  volume       = {76},
  year         = {2022},
}

@article{12247,
  abstract     = {Chromosomal inversions have been shown to play a major role in a local adaptation by suppressing recombination between alternative arrangements and maintaining beneficial allele combinations. However, so far, their importance relative to the remaining genome remains largely unknown. Understanding the genetic architecture of adaptation requires better estimates of how loci of different effect sizes contribute to phenotypic variation. Here, we used three Swedish islands where the marine snail Littorina saxatilis has repeatedly evolved into two distinct ecotypes along a habitat transition. We estimated the contribution of inversion polymorphisms to phenotypic divergence while controlling for polygenic effects in the remaining genome using a quantitative genetics framework. We confirmed the importance of inversions but showed that contributions of loci outside inversions are of similar magnitude, with variable proportions dependent on the trait and the population. Some inversions showed consistent effects across all sites, whereas others exhibited site-specific effects, indicating that the genomic basis for replicated phenotypic divergence is only partly shared. The contributions of sexual dimorphism as well as environmental factors to phenotypic variation were significant but minor compared to inversions and polygenic background. Overall, this integrated approach provides insight into the multiple mechanisms contributing to parallel phenotypic divergence.},
  author       = {Koch, Eva L. and Ravinet, Mark and Westram, Anja M and Johannesson, Kerstin and Butlin, Roger K.},
  issn         = {1558-5646},
  journal      = {Evolution},
  keywords     = {General Agricultural and Biological Sciences, Genetics, Ecology, Evolution, Behavior and Systematics},
  number       = {10},
  pages        = {2332--2346},
  publisher    = {Wiley},
  title        = {{Genetic architecture of repeated phenotypic divergence in Littorina saxatilis evolution}},
  doi          = {10.1111/evo.14602},
  volume       = {76},
  year         = {2022},
}

@article{12264,
  abstract     = {Reproductive isolation (RI) is a core concept in evolutionary biology. It has been the central focus of speciation research since the modern synthesis and is the basis by which biological species are defined. Despite this, the term is used in seemingly different ways, and attempts to quantify RI have used very different approaches. After showing that the field lacks a clear definition of the term, we attempt to clarify key issues, including what RI is, how it can be quantified in principle, and how it can be measured in practice. Following other definitions with a genetic focus, we propose that RI is a quantitative measure of the effect that genetic differences between populations have on gene flow. Specifically, RI compares the flow of neutral alleles in the presence of these genetic differences to the flow without any such differences. RI is thus greater than zero when genetic differences between populations reduce the flow of neutral alleles between populations. We show how RI can be quantified in a range of scenarios. A key conclusion is that RI depends strongly on circumstances—including the spatial, temporal and genomic context—making it difficult to compare across systems. After reviewing methods for estimating RI from data, we conclude that it is difficult to measure in practice. We discuss our findings in light of the goals of speciation research and encourage the use of methods for estimating RI that integrate organismal and genetic approaches.},
  author       = {Westram, Anja M and Stankowski, Sean and Surendranadh, Parvathy and Barton, Nicholas H},
  issn         = {1420-9101},
  journal      = {Journal of Evolutionary Biology},
  keywords     = {Ecology, Evolution, Behavior and Systematics},
  number       = {9},
  pages        = {1143--1164},
  publisher    = {Wiley},
  title        = {{What is reproductive isolation?}},
  doi          = {10.1111/jeb.14005},
  volume       = {35},
  year         = {2022},
}

@article{12265,
  author       = {Westram, Anja M and Stankowski, Sean and Surendranadh, Parvathy and Barton, Nicholas H},
  issn         = {1420-9101},
  journal      = {Journal of Evolutionary Biology},
  keywords     = {Ecology, Evolution, Behavior and Systematics},
  number       = {9},
  pages        = {1200--1205},
  publisher    = {Wiley},
  title        = {{Reproductive isolation, speciation, and the value of disagreement: A reply to the commentaries on ‘What is reproductive isolation?’}},
  doi          = {10.1111/jeb.14082},
  volume       = {35},
  year         = {2022},
}

@misc{13066,
  abstract     = {Chromosomal inversions have been shown to play a major role in local adaptation by suppressing recombination between alternative arrangements and maintaining beneficial allele combinations. However, so far, their importance relative to the remaining genome remains largely unknown. Understanding the genetic architecture of adaptation requires better estimates of how loci of different effect sizes contribute to phenotypic variation. Here, we used three Swedish islands where the marine snail Littorina saxatilis has repeatedly evolved into two distinct ecotypes along a habitat transition. We estimated the contribution of inversion polymorphisms to phenotypic divergence while controlling for polygenic effects in the remaining genome using a quantitative genetics framework. We confirmed the importance of inversions but showed that contributions of loci outside inversions are of similar magnitude, with variable proportions dependent on the trait and the population. Some inversions showed consistent effects across all sites, whereas others exhibited site-specific effects, indicating that the genomic basis for replicated phenotypic divergence is only partly shared. The contributions of sexual dimorphism as well as environmental factors to phenotypic variation were significant but minor compared to inversions and polygenic background. Overall, this integrated approach provides insight into the multiple mechanisms contributing to parallel phenotypic divergence.},
  author       = {Koch, Eva and Ravinet, Mark and Westram, Anja M and Jonannesson, Kerstin and Butlin, Roger},
  publisher    = {Dryad},
  title        = {{Data from: Genetic architecture of repeated phenotypic divergence in Littorina saxatilis ecotype evolution}},
  doi          = {10.5061/DRYAD.M905QFV4B},
  year         = {2022},
}

@article{10604,
  abstract     = {Maternally inherited Wolbachia transinfections are being introduced into natural mosquito populations to reduce the transmission of dengue, Zika, and other arboviruses. Wolbachia-induced cytoplasmic incompatibility provides a frequency-dependent reproductive advantage to infected females that can spread transinfections within and among populations. However, because transinfections generally reduce host fitness, they tend to spread within populations only after their frequency exceeds a critical threshold. This produces bistability with stable equilibrium frequencies at both 0 and 1, analogous to the bistability produced by underdominance between alleles or karyotypes and by population dynamics under Allee effects. Here, we analyze how stochastic frequency variation produced by finite population size can facilitate the local spread of variants with bistable dynamics into areas where invasion is unexpected from deterministic models. Our exemplar is the establishment of wMel Wolbachia in the Aedes aegypti population of Pyramid Estates (PE), a small community in far north Queensland, Australia. In 2011, wMel was stably introduced into Gordonvale, separated from PE by barriers to A. aegypti dispersal. After nearly 6 years during which wMel was observed only at low frequencies in PE, corresponding to an apparent equilibrium between immigration and selection, wMel rose to fixation by 2018. Using analytic approximations and statistical analyses, we demonstrate that the observed fixation of wMel at PE is consistent with both stochastic transition past an unstable threshold frequency and deterministic transformation produced by steady immigration at a rate just above the threshold required for deterministic invasion. The indeterminacy results from a delicate balance of parameters needed to produce the delayed transition observed. Our analyses suggest that once Wolbachia transinfections are established locally through systematic introductions, stochastic “threshold crossing” is likely to only minimally enhance spatial spread, providing a local ratchet that slightly—but systematically—aids area-wide transformation of disease-vector populations in heterogeneous landscapes.},
  author       = {Turelli, Michael and Barton, Nicholas H},
  issn         = {2056-3744},
  journal      = {Evolution Letters},
  keywords     = {genetics, ecology, evolution, behavior and systematics},
  number       = {1},
  pages        = {92--105},
  publisher    = {Wiley},
  title        = {{Why did the Wolbachia transinfection cross the road? Drift, deterministic dynamics, and disease control}},
  doi          = {10.1002/evl3.270},
  volume       = {6},
  year         = {2022},
}

@article{10658,
  abstract     = {We analyse how migration from a large mainland influences genetic load and population numbers on an island, in a scenario where fitness-affecting variants are unconditionally deleterious, and where numbers decline with increasing load. Our analysis shows that migration can have qualitatively different effects, depending on the total mutation target and fitness effects of deleterious variants. In particular, we find that populations exhibit a genetic Allee effect across a wide range of parameter combinations, when variants are partially recessive, cycling between low-load (large-population) and high-load (sink) states. Increased migration reduces load in the sink state (by increasing heterozygosity) but further inflates load in the large-population state (by hindering purging). We identify various critical parameter thresholds at which one or other stable state collapses, and discuss how these thresholds are influenced by the genetic versus demographic effects of migration. Our analysis is based on a ‘semi-deterministic’ analysis, which accounts for genetic drift but neglects demographic stochasticity. We also compare against simulations which account for both demographic stochasticity and drift. Our results clarify the importance of gene flow as a key determinant of extinction risk in peripheral populations, even in the absence of ecological gradients. This article is part of the theme issue ‘Species’ ranges in the face of changing environments (part I)’.},
  author       = {Sachdeva, Himani and Olusanya, Oluwafunmilola O and Barton, Nicholas H},
  issn         = {1471-2970},
  journal      = {Philosophical Transactions of the Royal Society B},
  number       = {1846},
  publisher    = {The Royal Society},
  title        = {{Genetic load and extinction in peripheral populations: The roles of migration, drift and demographic stochasticity}},
  doi          = {10.1098/rstb.2021.0010},
  volume       = {377},
  year         = {2022},
}

@article{10736,
  abstract     = {Predicting function from sequence is a central problem of biology. Currently, this is possible only locally in a narrow mutational neighborhood around a wildtype sequence rather than globally from any sequence. Using random mutant libraries, we developed a biophysical model that accounts for multiple features of σ70 binding bacterial promoters to predict constitutive gene expression levels from any sequence. We experimentally and theoretically estimated that 10–20% of random sequences lead to expression and ~80% of non-expressing sequences are one mutation away from a functional promoter. The potential for generating expression from random sequences is so pervasive that selection acts against σ70-RNA polymerase binding sites even within inter-genic, promoter-containing regions. This pervasiveness of σ70-binding sites implies that emergence of promoters is not the limiting step in gene regulatory evolution. Ultimately, the inclusion of novel features of promoter function into a mechanistic model enabled not only more accurate predictions of gene expression levels, but also identified that promoters evolve more rapidly than previously thought.},
  author       = {Lagator, Mato and Sarikas, Srdjan and Steinrueck, Magdalena and Toledo-Aparicio, David and Bollback, Jonathan P and Guet, Calin C and Tkačik, Gašper},
  issn         = {2050-084X},
  journal      = {eLife},
  publisher    = {eLife Sciences Publications},
  title        = {{Predicting bacterial promoter function and evolution from random sequences}},
  doi          = {10.7554/eLife.64543},
  volume       = {11},
  year         = {2022},
}

@article{10787,
  abstract     = {A species distributed across diverse environments may adapt to local conditions. We ask how quickly such a species changes its range in response to changed conditions. Szép et al. (Szép E, Sachdeva H, Barton NH. 2021 Polygenic local adaptation in metapopulations: a stochastic eco-evolutionary model. Evolution75, 1030–1045 (doi:10.1111/evo.14210)) used the infinite island model to find the stationary distribution of allele frequencies and deme sizes. We extend this to find how a metapopulation responds to changes in carrying capacity, selection strength, or migration rate when deme sizes are fixed. We further develop a ‘fixed-state’ approximation. Under this approximation, polymorphism is only possible for a narrow range of habitat proportions when selection is weak compared to drift, but for a much wider range otherwise. When rates of selection or migration relative to drift change in a single deme of the metapopulation, the population takes a time of order m−1 to reach the new equilibrium. However, even with many loci, there can be substantial fluctuations in net adaptation, because at each locus, alleles randomly get lost or fixed. Thus, in a finite metapopulation, variation may gradually be lost by chance, even if it would persist in an infinite metapopulation. When conditions change across the whole metapopulation, there can be rapid change, which is predicted well by the fixed-state approximation. This work helps towards an understanding of how metapopulations extend their range across diverse environments.
This article is part of the theme issue ‘Species’ ranges in the face of changing environments (Part II)’.},
  author       = {Barton, Nicholas H and Olusanya, Oluwafunmilola O},
  issn         = {1471-2970},
  journal      = {Philosophical Transactions of the Royal Society B: Biological Sciences},
  keywords     = {General Agricultural and Biological Sciences, General Biochemistry, Genetics and Molecular Biology},
  number       = {1848},
  publisher    = {The Royal Society},
  title        = {{The response of a metapopulation to a changing environment}},
  doi          = {10.1098/rstb.2021.0009},
  volume       = {377},
  year         = {2022},
}

@phdthesis{11128,
  abstract     = {Although we often see studies focusing on simple or even discrete traits in studies of colouration,
the variation of “appearance” phenotypes found in nature is often more complex, continuous
and high-dimensional. Therefore, we developed automated methods suitable for large datasets
of genomes and images, striving to account for their complex nature, while minimising human
bias. We used these methods on a dataset of more than 20, 000 plant SNP genomes and
corresponding fower images from a hybrid zone of two subspecies of Antirrhinum majus with
distinctly coloured fowers to improve our understanding of the genetic nature of the fower
colour in our study system.
Firstly, we use the advantage of large numbers of genotyped plants to estimate the haplotypes in
the main fower colour regulating region. We study colour- and geography-related characteristics
of the estimated haplotypes and how they connect to their relatedness. We show discrepancies
from the expected fower colour distributions given the genotype and identify particular
haplotypes leading to unexpected phenotypes. We also confrm a signifcant defcit of the
double recessive recombinant and quite surprisingly, we show that haplotypes of the most
frequent parental type are much less variable than others.
Secondly, we introduce our pipeline capable of processing tens of thousands of full fower
images without human interaction and summarising each image into a set of informative scores.
We show the compatibility of these machine-measured fower colour scores with the previously
used manual scores and study impact of external efect on the resulting scores. Finally, we use
the machine-measured fower colour scores to ft and examine a phenotype cline across the
hybrid zone in Planoles using full fower images as opposed to discrete, manual scores and
compare it with the genotypic cline.},
  author       = {Matejovicova, Lenka},
  isbn         = {978-3-99078-016-9},
  issn         = {2663-337X},
  pages        = {112},
  publisher    = {Institute of Science and Technology Austria},
  title        = {{Genetic basis of flower colour as a model for adaptive evolution}},
  doi          = {10.15479/at:ista:11128},
  year         = {2022},
}

@misc{11321,
  abstract     = {Here are the research data underlying the publication "Effects of fine-scale population structure on the distribution of heterozygosity in a long-term study of Antirrhinum majus" Further information are summed up in the README document. },
  author       = {Surendranadh, Parvathy and Arathoon, Louise S and Baskett, Carina and Field, David and Pickup, Melinda and Barton, Nicholas H},
  publisher    = {Institute of Science and Technology Austria},
  title        = {{Effects of fine-scale population structure on the distribution of heterozygosity in a long-term study of Antirrhinum majus}},
  doi          = {10.15479/at:ista:11321},
  year         = {2022},
}

@article{11334,
  abstract     = {Hybridization is a common evolutionary process with multiple possible outcomes. In vertebrates, interspecific hybridization has repeatedly generated parthenogenetic hybrid species. However, it is unknown whether the generation of parthenogenetic hybrids is a rare outcome of frequent hybridization between sexual species within a genus or the typical outcome of rare hybridization events. Darevskia is a genus of rock lizards with both hybrid parthenogenetic and sexual species. Using capture sequencing, we estimate phylogenetic relationships and gene flow among the sexual species, to determine how introgressive hybridization relates to the origins of parthenogenetic hybrids. We find evidence for widespread hybridization with gene flow, both between recently diverged species and deep branches. Surprisingly, we find no signal of gene flow between parental species of the parthenogenetic hybrids, suggesting that the parental pairs were either reproductively or geographically isolated early in their divergence. The generation of parthenogenetic hybrids in Darevskia is, then, a rare outcome of the total occurrence of hybridization within the genus, but the typical outcome when specific species pairs hybridize. Our results question the conventional view that parthenogenetic lineages are generated by hybridization in a window of divergence. Instead, they suggest that some lineages possess specific properties that underpin successful parthenogenetic reproduction.},
  author       = {Freitas, Susana and Westram, Anja M and Schwander, Tanja and Arakelyan, Marine and Ilgaz, Çetin and Kumlutas, Yusuf and Harris, David James and Carretero, Miguel A. and Butlin, Roger K.},
  issn         = {1558-5646},
  journal      = {Evolution},
  number       = {5},
  pages        = {899--914},
  publisher    = {Wiley},
  title        = {{Parthenogenesis in Darevskia lizards: A rare outcome of common hybridization, not a common outcome of rare hybridization}},
  doi          = {10.1111/evo.14462},
  volume       = {76},
  year         = {2022},
}

@phdthesis{11388,
  abstract     = {In evolve and resequence experiments, a population is sequenced, subjected to selection and
then sequenced again, so that genetic changes before and after selection can be observed at
the genetic level. Here, I use these studies to better understand the genetic basis of complex
traits - traits which depend on more than a few genes.
In the first chapter, I discuss the first evolve and resequence experiment, in which a population
of mice, the so-called "Longshanks" mice, were selected for tibia length while their body mass
was kept constant. The full pedigree is known. We observed a selection response on all
chromosomes and used the infinitesimal model with linkage, a model which assumes an infinite
number of genes with infinitesimally small effect sizes, as a null model. Results implied a very
polygenic basis with a few loci of major effect standing out and changing in parallel. There
was large variability between the different chromosomes in this study, probably due to LD.
In chapter two, I go on to discuss the impact of LD, on the variability in an allele-frequency
based summary statistic, giving an equation based on the initial allele frequencies, average
pairwise LD, and the first four moments of the haplotype block copy number distribution. I
describe this distribution by referring back to the founder generation. I then demonstrate
how to infer selection via a maximum likelihood scheme on the example of a single locus and
discuss how to extend this to more realistic scenarios.
In chapter three, I discuss the second evolve and resequence experiment, in which a small
population of Drosophila melanogaster was selected for increased pupal case size over 6
generations. The experiment was highly replicated with 27 lines selected within family and a
known pedigree. We observed a phenotypic selection response of over one standard deviation.
I describe the patterns in allele frequency data, including allele frequency changes and patterns
of heterozygosity, and give ideas for future work.},
  author       = {Belohlavy, Stefanie},
  isbn         = {978-3-99078-018-3},
  pages        = {98},
  publisher    = {Institute of Science and Technology Austria},
  title        = {{The genetic basis of complex traits studied via analysis of evolve and resequence experiments}},
  doi          = {10.15479/at:ista:11388},
  year         = {2022},
}

@article{11411,
  abstract     = {Many studies have quantified the distribution of heterozygosity and relatedness in natural populations, but few have examined the demographic processes driving these patterns. In this study, we take a novel approach by studying how population structure affects both pairwise identity and the distribution of heterozygosity in a natural population of the self-incompatible plant Antirrhinum majus. Excess variance in heterozygosity between individuals is due to identity disequilibrium, which reflects the variance in inbreeding between individuals; it is measured by the statistic g2. We calculated g2 together with FST and pairwise relatedness (Fij) using 91 SNPs in 22,353 individuals collected over 11 years. We find that pairwise Fij declines rapidly over short spatial scales, and the excess variance in heterozygosity between individuals reflects significant variation in inbreeding. Additionally, we detect an excess of individuals with around half the average heterozygosity, indicating either selfing or matings between close relatives. We use 2 types of simulation to ask whether variation in heterozygosity is consistent with fine-scale spatial population structure. First, by simulating offspring using parents drawn from a range of spatial scales, we show that the known pollen dispersal kernel explains g2. Second, we simulate a 1,000-generation pedigree using the known dispersal and spatial distribution and find that the resulting g2 is consistent with that observed from the field data. In contrast, a simulated population with uniform density underestimates g2, indicating that heterogeneous density promotes identity disequilibrium. Our study shows that heterogeneous density and leptokurtic dispersal can together explain the distribution of heterozygosity.},
  author       = {Surendranadh, Parvathy and Arathoon, Louise S and Baskett, Carina and Field, David and Pickup, Melinda and Barton, Nicholas H},
  issn         = {1943-2631},
  journal      = {Genetics},
  number       = {3},
  publisher    = {Oxford University Press},
  title        = {{Effects of fine-scale population structure on the distribution of heterozygosity in a long-term study of Antirrhinum majus}},
  doi          = {10.1093/genetics/iyac083},
  volume       = {221},
  year         = {2022},
}

@article{11447,
  abstract     = {Empirical essays of fitness landscapes suggest that they may be rugged, that is having multiple fitness peaks. Such fitness landscapes, those that have multiple peaks, necessarily have special local structures, called reciprocal sign epistasis (Poelwijk et al. in J Theor Biol 272:141–144, 2011). Here, we investigate the quantitative relationship between the number of fitness peaks and the number of reciprocal sign epistatic interactions. Previously, it has been shown (Poelwijk et al. in J Theor Biol 272:141–144, 2011) that pairwise reciprocal sign epistasis is a necessary but not sufficient condition for the existence of multiple peaks. Applying discrete Morse theory, which to our knowledge has never been used in this context, we extend this result by giving the minimal number of reciprocal sign epistatic interactions required to create a given number of peaks.},
  author       = {Saona Urmeneta, Raimundo J and Kondrashov, Fyodor and Khudiakova, Kseniia},
  issn         = {1522-9602},
  journal      = {Bulletin of Mathematical Biology},
  keywords     = {Computational Theory and Mathematics, General Agricultural and Biological Sciences, Pharmacology, General Environmental Science, General Biochemistry, Genetics and Molecular Biology, General Mathematics, Immunology, General Neuroscience},
  number       = {8},
  publisher    = {Springer Nature},
  title        = {{Relation between the number of peaks and the number of reciprocal sign epistatic interactions}},
  doi          = {10.1007/s11538-022-01029-z},
  volume       = {84},
  year         = {2022},
}

@article{11546,
  abstract     = {Local adaptation leads to differences between populations within a species. In many systems, similar environmental contrasts occur repeatedly, sometimes driving parallel phenotypic evolution. Understanding the genomic basis of local adaptation and parallel evolution is a major goal of evolutionary genomics. It is now known that by preventing the break-up of favourable combinations of alleles across multiple loci, genetic architectures that reduce recombination, like chromosomal inversions, can make an important contribution to local adaptation. However, little is known about whether inversions also contribute disproportionately to parallel evolution. Our aim here is to highlight this knowledge gap, to showcase existing studies, and to illustrate the differences between genomic architectures with and without inversions using simple models. We predict that by generating stronger effective selection, inversions can sometimes speed up the parallel adaptive process or enable parallel adaptation where it would be impossible otherwise, but this is highly dependent on the spatial setting. We highlight that further empirical work is needed, in particular to cover a broader taxonomic range and to understand the relative importance of inversions compared to genomic regions without inversions.},
  author       = {Westram, Anja M and Faria, Rui and Johannesson, Kerstin and Butlin, Roger and Barton, Nicholas H},
  issn         = {1471-2970},
  journal      = {Philosophical Transactions of the Royal Society B: Biological Sciences},
  keywords     = {General Agricultural and Biological Sciences, General Biochemistry, Genetics and Molecular Biology},
  number       = {1856},
  publisher    = {Royal Society of London},
  title        = {{Inversions and parallel evolution}},
  doi          = {10.1098/rstb.2021.0203},
  volume       = {377},
  year         = {2022},
}

