How Important Is Variation in Extrinsic Reproductive Isolation to the Process of Speciation?

The strength of reproductive isolation (RI) between two or more lineages during the process of speciation can vary by the ecological conditions. However, most speciation research has been limited to studying how ecologically dependent RI varies among a handful of broadly categorized environments. Very few studies consider the variability of RI and its effects on speciation at finer scales-that is, within each environment due to spatial or temporal environmental heterogeneity. Such variation in RI across time and/or space may inhibit speciation through leaky reproductive barriers or promote speciation by facilitating reinforcement. To investigate this overlooked aspect of speciation research, we conducted a literature review of existing studies of variation in RI in the field and then conducted individual-based simulations to examine how variation in hybrid fitness across time and space affects the degree of gene flow. Our simulations indicate that the presence of variation in hybrid fitness across space and time often leads to an increase in gene flow compared to scenarios where hybrid fitness remains static. This observation can be attributed to the convex relationship between the degree of gene flow and the strength of selection on hybrids. Our simulations also show that the effect of variation in RI on facilitating gene flow is most pronounced when RI, on average, is relatively low. This finding suggests that it could serve as an important mechanism to explain why the completion of speciation is often challenging. While direct empirical evidence documenting variation in extrinsic RI is limited, we contend that it is a prevalent yet underexplored phenomenon. We support this argument by proposing common scenarios in which RI is likely to exhibit variability and thus influence the process of speciation.

Effects of environment and genotype on dispersal differ across departure, transfer and settlement in a butterfly metapopulation.

Active dispersal is driven by extrinsic and intrinsic factors at the three stages of departure, transfer and settlement. Most empirical studies capture only one stage of this complex process, and knowledge of how much can be generalized from one stage to another remains unknown. Here we use genetic assignment tests to reconstruct dispersal across 5 years and 232 habitat patches of a Glanville fritillary butterfly (Melitaea cinxia) metapopulation. We link individual dispersal events to weather, landscape structure, size and quality of habitat patches, and individual genotype to identify the factors that influence the three stages of dispersal and post-settlement survival. We found that nearly all tested factors strongly affected departure probabilities, but that the same factors explained very little variation in realized dispersal distances. Surprisingly, we found no effect of dispersal distance on post-settlement survival. Rather, survival was influenced by weather conditions, quality of the natal habitat patch, and a strong interaction between genotype and occupancy status of the settled habitat patch, with more mobile genotypes having higher survival as colonists rather than as immigrants. Our work highlights the multi-causality of dispersal and that some dispersal costs can only be understood by considering extrinsic and intrinsic factors and their interaction across the entire dispersal process.

A modified niche model for generating food webs with stage-structured consumers: The stabilizing effects of life-history stages on complex food webs.

Almost all organisms grow in size during their lifetime and switch diets, trophic positions, and interacting partners as they grow. Such ontogenetic development introduces life-history stages and flows of biomass between the stages through growth and reproduction. However, current research on complex food webs rarely considers life-history stages. The few previously proposed methods do not take full advantage of the existing food web structural models that can produce realistic food web topologies.We extended the niche model developed by Williams and Martinez (Nature, 2000, 404, 180-183) to generate food webs that included trophic species with a life-history stage structure. Our method aggregated trophic species based on niche overlap to form a life-history structured population; therefore, it largely preserved the topological structure of food webs generated by the niche model. We applied the theory of allometric predator-prey body mass ratio and parameterized an allometric bioenergetic model augmented with biomass flow between stages via growth and reproduction to study the effects of a stage structure on the stability of food webs.When life-history stages were linked via growth and reproduction, more food webs persisted, and persisting food webs tended to retain more trophic species. Topological differences between persisting linked and unlinked food webs were small to modest. The slopes of biomass spectra were lower, and weak interaction links were more prevalent in the linked food webs than the unlinked ones, suggesting that a life-history stage structure promotes characteristics that can enhance stability of complex food webs.Our results suggest a positive relationship between the complexity and stability of complex food webs. A life-history stage structure in food webs may play important roles in dynamics of and diversity in food webs.

The effects of functional response and host abundance fluctuations on genetic rescue in parasitoids with single-locus sex determination.

Many parasitoids have single-locus complementary sex determination (sl-CSD), which produces sterile or inviable males when homozygous at the sex determining locus. A previous study theoretically showed that small populations have elevated risks of extinction due to the positive feedback between inbreeding and small population size, referred to as the diploid male vortex. A few modeling studies have suggested that the diploid male vortex may not be as common because balancing selection at sex determining loci tends to maintain high allelic diversity in spatially structured populations. However, the generality of the conclusion is yet uncertain, as they were drawn either from models developed for particular systems or from a general-purpose competition model. To attest the conclusion, we study several well-studied host-parasitoid models that incorporate functional response specifying the number of attacked hosts given a host density and derive the conditions for a diploid male vortex in a single population. Then, we develop spatially structured individual-based versions of the models to include female behavior, diploid male fertility, and temporal fluctuations. The results show that producing a handful of successful offspring per female parasitoid could enable parasitoid persistence when a typical number of CSD alleles are present. The effect of functional response depends on the levels of fluctuations in host abundance, and inviable or partially fertile diploid males and a small increase in dispersal can alleviate the risk of a diploid male vortex. Our work supports the generality of effective genetic rescue in spatially connected parasitoid populations with sl-CSD. However, under more variable climate, the efficacy of the CSD mechanism may substantially decline.

Scaling up the effects of inbreeding depression from individuals to metapopulations.

Inbreeding is common in nature, and many laboratory studies have documented that inbreeding depression can reduce the fitness of individuals. Demonstrating the consequences of inbreeding depression on the growth and persistence of populations is more challenging because populations are often regulated by density- or frequency-dependent selection and influenced by demographic and environmental stochasticity. A few empirical studies have shown that inbreeding depression can increase extinction risk of local populations. The importance of inbreeding depression at the metapopulation level has been conjectured based on population-level studies but has not been evaluated. We quantified the impact of inbreeding depression affecting the fitness of individuals on metapopulation persistence in heterogeneous habitat networks of different sizes and habitat configuration in a context of natural butterfly metapopulations. We developed a spatial individual-based simulation model of metapopulations with explicit genetics. We used Approximate Bayesian Computation to fit the model to extensive demographic, genetic and life-history data available for the well-studied Glanville fritillary butterfly (Melitaea cinxia) metapopulations in the Åland islands in SW Finland. We compared 18 semi-independent habitat networks differing in size and fragmentation. The results show that inbreeding is more frequent in small habitat networks, and consequently, inbreeding depression elevates extinction risks in small metapopulations. Metapopulation persistence and neutral genetic diversity maintained in the metapopulations increase with the total habitat amount in and mean patch size of habitat networks. Dispersal and mating behaviour interact with landscape structure to determine how likely it is to encounter kin while looking for mates. Inbreeding depression can decrease the viability of small metapopulations even when they are strongly influenced by stochastic extinction-colonization dynamics and density-dependent selection. The findings from this study support that genetic factors, in addition to demographic factors, can contribute to extinctions of small local populations and also of metapopulations.

Increased fluctuation in a butterfly metapopulation leads to diploid males and decline of a hyperparasitoid.

Climate change can increase spatial synchrony of population dynamics, leading to large-scale fluctuation that destabilizes communities. High trophic level species such as parasitoids are disproportionally affected because they depend on unstable resources. Most parasitoid wasps have complementary sex determination, producing sterile males when inbred, which can theoretically lead to population extinction via the diploid male vortex (DMV). We examined this process empirically using a hyperparasitoid population inhabiting a spatially structured host population in a large fragmented landscape. Over four years of high host butterfly metapopulation fluctuation, diploid male production by the wasp increased, and effective population size declined precipitously. Our multitrophic spatially structured model shows that host population fluctuation can cause local extinctions of the hyperparasitoid because of the DMV. However, regionally it persists because spatial structure allows for efficient local genetic rescue via balancing selection for rare alleles carried by immigrants. This is, to our knowledge, the first empirically based study of the possibility of the DMV in a natural host-parasitoid system.

Genetics of dispersal.

Dispersal is a process of central importance for the ecological and evolutionary dynamics of populations and communities, because of its diverse consequences for gene flow and demography. It is subject to evolutionary change, which begs the question, what is the genetic basis of this potentially complex trait? To address this question, we (i) review the empirical literature on the genetic basis of dispersal, (ii) explore how theoretical investigations of the evolution of dispersal have represented the genetics of dispersal, and (iii) discuss how the genetic basis of dispersal influences theoretical predictions of the evolution of dispersal and potential consequences. Dispersal has a detectable genetic basis in many organisms, from bacteria to plants and animals. Generally, there is evidence for significant genetic variation for dispersal or dispersal-related phenotypes or evidence for the micro-evolution of dispersal in natural populations. Dispersal is typically the outcome of several interacting traits, and this complexity is reflected in its genetic architecture: while some genes of moderate to large effect can influence certain aspects of dispersal, dispersal traits are typically polygenic. Correlations among dispersal traits as well as between dispersal traits and other traits under selection are common, and the genetic basis of dispersal can be highly environment-dependent. By contrast, models have historically considered a highly simplified genetic architecture of dispersal. It is only recently that models have started to consider multiple loci influencing dispersal, as well as non-additive effects such as dominance and epistasis, showing that the genetic basis of dispersal can influence evolutionary rates and outcomes, especially under non-equilibrium conditions. For example, the number of loci controlling dispersal can influence projected rates of dispersal evolution during range shifts and corresponding demographic impacts. Incorporating more realism in the genetic architecture of dispersal is thus necessary to enable models to move beyond the purely theoretical towards making more useful predictions of evolutionary and ecological dynamics under current and future environmental conditions. To inform these advances, empirical studies need to answer outstanding questions concerning whether specific genes underlie dispersal variation, the genetic architecture of context-dependent dispersal phenotypes and behaviours, and correlations among dispersal and other traits.

Expansion of cooperatively growing populations: Optimal migration rates and habitat network structures.

Range expansion of species is driven by the interactions among individual- and population-level processes and the spatial pattern of habitats. In this work we study how cooperatively growing populations spread on networks representing the skeleton of complex landscapes. By separating the slow and fast variables of the expansion process, we are able to give analytical predictions for the critical conditions that divide the dynamic behaviors into different phases (extinction, localized survival, and global expansion). We observe a resonance phenomenon in how the critical condition depends on the expansion rate, indicating the existence of an optimal strategy for global expansion. We derive the conditions for such optimal migration in locally treelike graphs and numerically study other structured networks. Our results highlight the importance of both the underlying interaction pattern and migration rate of the expanding populations for range expansion. We also discuss potential applications of the results to biological control and conservation.

Mechanisms by Which Phenotypic Plasticity Affects Adaptive Divergence and Ecological Speciation.

Phenotypic plasticity is the ability of one genotype to produce different phenotypes depending on environmental conditions. Several conceptual models emphasize the role of plasticity in promoting reproductive isolation and, ultimately, speciation in populations that forage on two or more resources. These models predict that plasticity plays a critical role in the early stages of speciation, prior to genetic divergence, by facilitating fast phenotypic divergence. The ability to plastically express alternative phenotypes may, however, interfere with the early phase of the formation of reproductive barriers, especially in the absence of geographic barriers. Here, we quantitatively investigate mechanisms under which plasticity can influence progress toward adaptive genetic diversification and ecological speciation. We use a stochastic, individual-based model of a predator-prey system incorporating sexual reproduction and mate choice in the predator. Our results show that evolving plasticity promotes the evolution of reproductive isolation under diversifying environments when individuals are able to correctly select a more profitable habitat with respect to their phenotypes (i.e., adaptive habitat choice) and to assortatively mate with relatively similar phenotypes. On the other hand, plasticity facilitates the evolution of plastic generalists when individuals have a limited capacity for adaptive habitat choice. We conclude that plasticity can accelerate the evolution of a reproductive barrier toward adaptive diversification and ecological speciation through enhanced phenotypic differentiation between diverging phenotypes.

Genomics and the origin of species.

Speciation is a fundamental evolutionary process, the knowledge of which is crucial for understanding the origins of biodiversity. Genomic approaches are an increasingly important aspect of this research field. We review current understanding of genome-wide effects of accumulating reproductive isolation and of genomic properties that influence the process of speciation. Building on this work, we identify emergent trends and gaps in our understanding, propose new approaches to more fully integrate genomics into speciation research, translate speciation theory into hypotheses that are testable using genomic tools and provide an integrative definition of the field of speciation genomics.

Evolutionary suicide as a consequence of runaway selection for greater aggregation tendency.

Aggregation of individuals is a common phenomenon in nature. By aggregating, individuals can reap benefits but may also be subject to associated costs from increased competition. The benefits of aggregation can depend on population density, which in turn can be affected by aggregation when it determines reproductive success of individuals. The Allee effect is often considered to be one of the factors that can explain the evolution of aggregation behavior. We investigated this hypothesis with a mathematical model which integrates population dynamics and evolution. Individuals gain synergistically from aggregation but suffer from scramble competition with aggregation tendency as an evolving trait. We found that aggregation behavior can stabilize the population dynamics and reduce population growth. The results show that the Allee effect alone is not sufficient for aggregative behavior to evolve as an evolutionarily stable strategy. We also found that weak local competition does not promote aggregation due to feedback from the population level: under low competition, the population can achieve high density such that aggregation becomes costly rather than beneficial. Our model instead exhibits an escalation of aggregation tendency, leading to the extinction of the population in a process known as evolutionary suicide. We conclude that for aggregation to evolve as an evolutionarily stable strategy we need to consider other factors such as inter-patch dispersal to new patches and avoidance of excessively large groups.

Toward an integration of evolutionary biology and ecosystem science.

At present, the disciplines of evolutionary biology and ecosystem science are weakly integrated. As a result, we have a poor understanding of how the ecological and evolutionary processes that create, maintain, and change biological diversity affect the flux of energy and materials in global biogeochemical cycles. The goal of this article was to review several research fields at the interfaces between ecosystem science, community ecology and evolutionary biology, and suggest new ways to integrate evolutionary biology and ecosystem science. In particular, we focus on how phenotypic evolution by natural selection can influence ecosystem functions by affecting processes at the environmental, population and community scale of ecosystem organization. We develop an eco-evolutionary model to illustrate linkages between evolutionary change (e.g. phenotypic evolution of producer), ecological interactions (e.g. consumer grazing) and ecosystem processes (e.g. nutrient cycling). We conclude by proposing experiments to test the ecosystem consequences of evolutionary changes.

Persistence of pathogens with short infectious periods in seasonal tick populations: the relative importance of three transmission routes.

BACKGROUND: The flaviviruses causing tick-borne encephalitis (TBE) persist at low but consistent levels in tick populations, despite short infectious periods in their mammalian hosts and transmission periods constrained by distinctly seasonal tick life cycles. In addition to systemic and vertical transmission, cofeeding transmission has been proposed as an important route for the persistence of TBE-causing viruses. Because cofeeding transmission requires ticks to feed simultaneously, the timing of tick activity may be critical to pathogen persistence. Existing models of tick-borne diseases do not incorporate all transmission routes and tick seasonality. Our aim is to evaluate the influence of seasonality on the relative importance of different transmission routes by using a comprehensive mathematical model. METHODOLOGY/PRINCIPAL FINDINGS: We developed a stage-structured population model that includes tick seasonality and evaluated the relative importance of the transmission routes for pathogens with short infectious periods, in particular Powassan virus (POWV) and the related "deer tick virus," emergent encephalitis-causing flaviviruses in North America. We used the next generation matrix method to calculate the basic reproductive ratio and performed elasticity analyses. We confirmed that cofeeding transmission is critically important for such pathogens to persist in seasonal tick populations over the reasonable range of parameter values. At higher but still plausible rates of vertical transmission, our model suggests that vertical transmission can strongly enhance pathogen prevalence when it operates in combination with cofeeding transmission. CONCLUSIONS/SIGNIFICANCE: Our results demonstrate that the consistent prevalence of POWV observed in tick populations could be maintained by a combination of low vertical, intermediate cofeeding and high systemic transmission rates. When vertical transmission is weak, nymphal ticks support integral parts of the transmission cycle that are critical for maintaining the pathogen. We also extended the model to pathogens that cause chronic infections in hosts and found that cofeeding transmission could contribute to elevating prevalence even in these systems. Therefore, the common assumption that cofeeding transmission is not relevant in models of chronic host infection, such as Lyme disease, could lead to underestimating pathogen prevalence.