Partial Clonality Expands the Opportunity for Spatial Adaptation.

AbstractReproductive mode may strongly impact adaptation in spatially varying populations linked by dispersal, especially when sexual and clonal offspring differ in dispersal. We determined how spatial structure affects adaptation in populations with mixed clonal and sexual reproduction. In a source-sink quantitative genetic deterministic model (with stabilizing selection around different optima), greater clonal reproduction or parent-offspring association (a measure of the part of the parent's phenotype other than the additive genetic component inherited by clonal offspring) increased the selective difference (difference between phenotypic optima) allowing sink populations to adapt. Given dispersal differences between clonally and sexually produced juveniles, adaptation increased with an increasing fraction of clonal dispersers. When considering migrational meltdown, partially clonal reproduction reduced cases where dispersal caused habitat loss. Stochastic individual-based simulations support these results, although the effect of differential dispersal was reversed, with decreased clonal dispersal allowing greater adaptation. These results parallel earlier findings that for an instantaneous shift in phenotypic optimum, increasing clonality allowed population persistence for a greater shift; here, selective change is spatial rather than temporal. These results may help explain the success of many partially clonal organisms in invading new habitats, complementing traditional explanations based on avoiding Allee effects.

Environmental fluctuations dampen the effects of clonal reproduction on evolutionary rescue.

Evolutionary rescue occurs when genetic change allows a population to persist in response to an environmental change that would otherwise have led to extinction. Most studies of evolutionary rescue assume that species have either fully clonal or fully sexual reproduction; however, many species have partially clonal reproductive strategies in which they reproduce both clonally and sexually. Furthermore, the few evolutionary rescue studies that have evaluated partially clonal reproduction did not consider fluctuations in the environment, which are nearly ubiquitous in nature. Here, we use individual-based simulations to investigate how environmental fluctuations (either uncorrelated or positively autocorrelated) influence the effect of clonality on evolutionary rescue. We show that, for moderate magnitudes of environmental fluctuations, as was found in the absence of fluctuations, increasing the degree of clonality increases the probability of population persistence in response to an abrupt environmental change, but decreases persistence in response to a continuous, directional environmental change. However, with large magnitudes of fluctuations, both the benefits of clonality following a step change and the detrimental effects of clonality following a continuous, directional change are generally reduced; in fact, in the latter scenario, increasing clonality can even become beneficial if environmental fluctuations are autocorrelated. We also show that increased generational overlap dampens the effects of environmental fluctuations. Overall, we demonstrate that understanding the evolutionary rescue of partially clonal organisms requires not only knowledge of the species life history and the type of environmental change, but also an understanding of the magnitude and autocorrelation of environmental fluctuations.

Sex and Asex: A Clonal Lexicon.

Organisms across the tree of life have complex life cycles that include both sexual and asexual reproduction or that are obligately asexual. These organisms include ecologically dominant species that structure many terrestrial and marine ecosystems, as well as many pathogens, pests, and invasive species. We must consider both the evolution and maintenance of these various reproductive modes and how these modes shape the genetic diversity, adaptive evolution, and ability to persist in the species that exhibit them. Thus, having a common framework is a key aspect of understanding the biodiversity that shapes our planet. In the 2019 AGA President's Symposium, Sex and Asex: The genetics of complex life cycles, researchers investigating a wide range of taxonomic models and using a variety of modes of investigation coalesced around a common theme-understanding not only how such complex life cycles may evolve, but how they are shaped by the evolutionary and ecological forces around them. In this introduction to the Special Issue from the symposium, we give an overview of some of the key ideas and areas of investigation (a common clonal lexicon, we might say) and introduce the breadth of work submitted by symposium participants.

Evolutionary Rescue in a Linearly Changing Environment: Limits on Predictability.

Populations subject to substantial environmental change that decreases absolute fitness (expected number of offspring per individual) to less than one must adapt to persist. The probability of adaptive evolutionary rescue may be influenced by factors intrinsic to the organism itself, or by features specific to the individual population and its environment. An important question (given the increasing prevalence of environmental change) is the predictability of evolutionary rescue. We used an individual-based simulation model and a related analytic model to examine population persistence, given a continuously changing environment that leads to a linear change in the optimum for a phenotypic trait under selection. Population persistence was not well predicted by the population genetics at the start of environmental change, which contrasts strongly with the results shown in prior work for persistence after a sudden environmental change. Larger populations, which had a greater scope for the generation and maintenance of beneficial genetic variation, showed a clear advantage, but increasing the rate of environmental change always decreased the probability of persistence. Extinctions occurred throughout the period of continuous change, and populations that went extinct showed little sign of their eventual fate until shortly before extinction. Partially clonal populations showed less predictability and greater vulnerability to extinction when impacted by continuous change than did fully sexual populations-any advantage gained by the initial transmission of well-adapted phenotypes via clonal reproduction is lost as the phenotypic optimum continues to shift and the generation of novel variation is required for continuous adaptation.

Effects of Clonal Reproduction on Evolutionary Lag and Evolutionary Rescue.

Evolutionary lag-the difference between mean and optimal phenotype in the current environment-is of keen interest in light of rapid environmental change. Many ecologically important organisms have life histories that include stage structure and both sexual and clonal reproduction, yet how stage structure and clonality interplay to govern a population's rate of evolution and evolutionary lag is unknown. Effects of clonal reproduction on mean phenotype partition into two portions: one that is phenotype dependent, and another that is genotype dependent. This partitioning is governed by the association between the nonadditive genetic plus random environmental component of phenotype of clonal offspring and their parents. While clonality slows phenotypic evolution toward an optimum, it can dramatically increase population survival after a sudden step change in optimal phenotype. Increased adult survival slows phenotypic evolution but facilitates population survival after a step change; this positive effect can, however, be lost given survival-fecundity trade-offs. Simulations indicate that the benefits of increased clonality under environmental change greatly depend on the nature of that change: increasing population persistence under a step change while decreasing population persistence under a continuous linear change requiring de novo variation. The impact of clonality on the probability of persistence for species in a changing world is thus inexorably linked to the temporal texture of the change they experience.

Warmest extreme year in U.S. history alters thermal requirements for tree phenology.

The frequency of extreme warm years is increasing across the majority of the planet. Shifts in plant phenology in response to extreme years can influence plant survival, productivity, and synchrony with pollinators/herbivores. Despite extensive work on plant phenological responses to climate change, little is known about responses to extreme warm years, particularly at the intraspecific level. Here we investigate 43 populations of white ash trees (Fraxinus americana) from throughout the species range that were all grown in a common garden. We compared the timing of leaf emergence during the warmest year in U.S. history (2012) with relatively non-extreme years. We show that (a) leaf emergence among white ash populations was accelerated by 21 days on average during the extreme warm year of 2012 relative to non-extreme years; (b) rank order for the timing of leaf emergence was maintained among populations across extreme and non-extreme years, with southern populations emerging earlier than northern populations; (c) greater amounts of warming units accumulated prior to leaf emergence during the extreme warm year relative to non-extreme years, and this constrained the potential for even earlier leaf emergence by an average of 9 days among populations; and (d) the extreme warm year reduced the reliability of a relevant phenological model for white ash by producing a consistent bias toward earlier predicted leaf emergence relative to observations. These results demonstrate a critical need to better understand how extreme warm years will impact tree phenology, particularly at the intraspecific level.

Mutation-selection balance and mixed mating with asexual reproduction.

The effects of asexual reproduction on both the number of deleterious mutations per gamete and the mean fitness under mutation-selection balance are investigated. We use two simulation models, considering both finite and infinite populations. The two models incorporate asexual reproduction with varying levels of outcrossing and selfing, degrees of dominance and selection coefficients. The values for mean fitness and number of deleterious mutations per gamete are compared within and among finite and infinite populations to identify the effect of asexual reproduction on levels of load, and how asexual reproduction may interact with genetic drift (population size). Increasing asexual reproduction resulted in an increase in mean fitness and a decrease in the average number of deleterious mutations per gamete for both nearly recessive and additive alleles in both the infinite and finite simulations. Increased mean fitness with increasing asexuality is possibly due to two interacting forces: a greater opportunity for selection to act on heterozygous versus homozygous mutations and the shielding of a proportion of the population from meiotic mutations due to asexual reproduction. The results found here highlight the need to consider asexual reproduction along with mixed mating in models of genetic load and mutation-selection balance.

Who are we now? A demographic assessment of three evolution societies.

Scientific societies have the potential to catalyze support for communities that have been historically excluded from science. Many of these societies have formed committees to propose and administer initiatives to promote the career and well-being of their members, with a special emphasis on racial and ethnic minorities. Yet, these societies are rarely armed with data to inform their proposals. Three of the evolution societies (American Society of Naturalists, "ASN"; Society of Systematic Biologists, "SSB"; Society for the Study of Evolution, "SSE") have also formed Diversity, Equity, and Inclusion committees in the last few years. As a first step in determining the needs of the societies, these committees collected data on the demographic characteristics of the societies' constituents by surveying the attendants of the Evolution 2019 meeting. Here, we report the proportions for different demographic groups in attendance at the meeting and compare these proportions to the demographics of recipients of Ph.D. degrees either in evolutionary biology or in the broader life sciences, as well as population demographics of the USA. Our results indicate that historically excluded groups are still underrepresented across US-based evolutionary biology professional societies. We explore whether demographic composition differs at different professional stages and find that representation for women and LGBTQ+ members decreases as the career stage progresses. We also find some evidence for heterogeneity across societies in terms of racial composition. Finally, we discuss the caveats and limitations of our procedures. Our results will serve to inform future efforts to collect demographic data at the society levels, which should in turn be used to design and implement evidence-based initiatives for inclusion and equity. This report should be a starting point for systematic efforts to characterize the ever-changing representation in evolutionary biology and to work toward the inclusion of all groups.

Nonclonal coloniality: Genetically chimeric colonies through fusion of sexually produced polyps in the hydrozoan Ectopleura larynx.

Hydrozoans typically develop colonies through asexual budding of polyps. Although colonies of Ectopleura are similar to other hydrozoans in that they consist of multiple polyps physically connected through continuous epithelia and shared gastrovascular cavity, Ectopleura larynx does not asexually bud polyps indeterminately. Instead, after an initial phase of limited budding in a young colony, E. larynx achieves its large colony size through the aggregation and fusion of sexually (nonclonally) produced polyps. The apparent chimerism within a physiologically integrated colony presents a potential source of conflict between distinct genetic lineages, which may vary in their ability to access the germline. To determine the extent to which the potential for genetic conflict exists, we characterized the types of genetic relationships between polyps within colonies, using a RAD-Seq approach. Our results indicate that E. larynx colonies are indeed comprised of polyps that are clones and sexually reproduced siblings and offspring, consistent with their life history. In addition, we found that colonies also contain polyps that are genetically unrelated, and that estimates of genome-wide relatedness suggests a potential for conflict within a colony. Taken together, our data suggest that there are distinct categories of relationships in colonies of E. larynx, likely achieved through a range of processes including budding, regeneration, and fusion of progeny and unrelated polyps, with the possibility for a genetic conflict resolution mechanism. Together these processes contribute to the reevolution of the ecologically important trait of coloniality in E. larynx.

Associations between cytoplasmic and nuclear Loci in hybridizing populations.

We extend current multilocus models to describe the effects of migration, recombination, selection, and nonrandom mating on sets of genes in diploids with varied modes of inheritance, allowing us to consider the patterns of nuclear and cytonuclear associations (disequilibria) under various models of migration. We show the relationship between the multilocus notation recently presented by Kirkpatrick, Johnson, and Barton (developed from previous work by Barton and Turelli) and the cytonuclear parameterization of Asmussen, Arnold, and Avise and extend this notation to describe associations between cytoplasmic elements and multiple nuclear genes. Under models with sexual symmetry, both nuclear-nuclear and cytonuclear disequilibria are equivalent. They differ, however, in cases involving some type of sexual asymmetry, which is then reflected in the asymmetric inheritance of cytoplasmic markers. An example given is the case of different migration rates in males and females; simulations using 2, 3, 4, or 5 unlinked autosomal markers with a maternally inherited cytoplasmic marker illustrate how nuclear-nuclear and cytonuclear associations can be used to separately estimate female and male migration rates. The general framework developed here allows us to investigate conditions where associations between loci with different modes of inheritance are not equivalent and to use this nonequivalence to test for deviations from simple models of admixture.

The role of pathogen shedding in linking within- and between-host pathogen dynamics.

A model linking within- and between-host pathogen dynamics via pathogen shedding (emission of pathogens throughout the course of infection) is developed, and several aspects of host availability and co-infection are considered. In this model, the rate of pathogen shedding affects both the pathogen population size within a host (also affecting host mortality) and the rate of infection of new hosts. Our goal is to ascertain how the rate of shedding is likely to evolve, and what factors permit coexistence of alternative shedding rates in a pathogen population. For a constant host population size (where an increase in infected hosts necessarily decreases susceptible hosts), important differences arise depending on whether pathogens compete only for susceptible (uninfected) hosts, or whether co-infection allows for competition for infected hosts. With no co-infection, the pathogen type that can persist with the lowest number of susceptible hosts will outcompete any other, which under the assumptions of the model is the pathogen with the highest basic reproduction number. This is often a pathogen with a relatively high shedding rate (s). If within-host competition is allowed, a trade-off develops due to the conflicting effects of shedding on within- and between-host pathogen dynamics, with within-host competition favoring clones with low shedding rates while between-host competition benefits clones with higher shedding rates. With within-host competition for the same host cells, low shedding rate clones should eliminate high-s clones in a co-infected host, if equilibrium is reached. With co-infection, but no within-host competition, pathogen clones still interact by affecting the mortality of co-infected hosts; here, coexistence is more likely. With co-infection, two clones can coexist if one is the superior competitor for uninfected hosts and the other for co-infected hosts.

A statistical characterization of consistent patterns of human immunodeficiency virus evolution within infected patients.

Within-patient HIV populations evolve rapidly because of a high mutation rate, short generation time, and strong positive selection pressures. Previous studies have identified "consistent patterns" of viral sequence evolution. Just before HIV infection progresses to AIDS, evolution seems to slow markedly, and the genetic diversity of the viral population drops. This evolutionary slowdown could be caused either by a reduction in the average viral replication rate or because selection pressures weaken with the collapse of the immune system. The former hypothesis (which we denote "cellular exhaustion") predicts a simultaneous reduction in both synonymous and nonsynonymous evolution, whereas the latter hypothesis (denoted "immune relaxation") predicts that only nonsynonymous evolution will slow. In this paper, we present a set of statistical procedures for distinguishing between these alternative hypotheses using DNA sequences sampled over the course of infection. The first component is a new method for estimating evolutionary rates that takes advantage of the temporal information in longitudinal DNA sequence samples. Second, we develop a set of probability models for the analysis of evolutionary rates in HIV populations in vivo. Application of these models to both synonymous and nonsynonymous evolution affords a comparison of the cellular-exhaustion and immune-relaxation hypotheses. We apply the procedures to longitudinal data sets in which sequences of the env gene were sampled over the entire course of infection. Our analyses (1) statistically confirm that an evolutionary slowdown occurs late in infection, (2) strongly support the immune-relaxation hypothesis, and (3) indicate that the cessation of nonsynonymous evolution is associated with disease progression.

Viral infection in internally structured hosts. I. Conditions for persistent infection.

For a virus population within its host, two important levels of structure can be considered: multiple cell types which can be infected, and tissue types or body compartments which may be coupled via movement. We develop a model with both types of structure. Migration between compartments can create "sources" and "sinks" within the virus population, where realized viral growth rate and abundance is lowered in some compartments compared to what would be observed in isolation. Using both analytical and numerical methods, we investigate how this within-host spatial structure affects the conditions for persistent viral infection. We find that migration between compartments makes the establishment of infection more difficult than it would be in the absence of migration, implying that within-host spatial structure combined with viral movement decreases the likelihood of viral establishment. If migration is symmetrical and compartments are heterogeneous, an increase in migration rates between compartments generally makes establishment less likely. This may help to explain the tissue specificity observed for many viruses. There are, however, important exceptions to this result. These include circumstances where the virus initially invades a compartment that is unfavorable to population growth and migration is necessary to infect other parts of the host body. Stochastic aspects of viral establishment may also favor increased migration as it tends to dampen the amplitude of fluctuations in population size during the initial transient phase of establishment.

The genealogy of a sequence subject to purifying selection at multiple sites.

We investigate the effect of purifying selection at multiple sites on both the shape of the genealogy and the distribution of mutations on the tree. We find that the primary effect of purifying selection on a genealogy is to shift the distribution of mutations on the tree, whereas the shape of the tree remains largely unchanged. This result is relevant to the large number of coalescent estimation procedures, which generally assume neutrality for segregating polymorphisms--applying these estimators to evolutionarily constrained sequences could lead to a significant degree of bias. We also estimate the statistical power of several neutrality tests in detecting weak to moderate purifying selection and find that the power is quite good for some parameter combinations. This result contrasts with previous studies, which predicted low statistical power because of the minor effect that weak purifying selection has on the shape of a genealogy. Finally, we investigate the effect of Hill-Robertson interference among linked deleterious mutations on patterns of molecular variation. We find that dependence among selected loci can substantially reduce the efficacy of even fairly strong purifying selection.

Linking dynamical and population genetic models of persistent viral infection.

This article develops a theoretical framework to link dynamical and population genetic models of persistent viral infection. This linkage is useful because, while the dynamical and population genetic theories have developed independently, the biological processes they describe are completely interrelated. Parameters of the dynamical models are important determinants of evolutionary processes such as natural selection and genetic drift. We develop analytical methods, based on coupled differential equations and Markov chain theory, to predict the accumulation of genetic diversity within the viral population as a function of dynamical parameters. These methods are first applied to the standard model of viral dynamics and then generalized to consider the infection of multiple host cell types by the viral population. Each cell type is characterized by specific parameter values. Inclusion of multiple cell types increases the likelihood of persistent infection and can increase the amount of genetic diversity within the viral population. However, the overall rate of gene sequence evolution may actually be reduced.