Do Interactions among Microbial Symbionts Cause Selection for Greater Pathogen Virulence?

AbstractThe ecological and evolutionary consequences of microbiome treatments aimed at protecting plants and animals against infectious disease are not well understood, even as such biological control measures become more common in agriculture and medicine. Notably, we lack information on the impacts of symbionts on pathogen fitness with which to project the consequences of competition for the evolution of virulence. To address this gap, we estimated fitness consequences for a common plant pathogen, Ustilago maydis, over differing virulence levels and when the host plant (Zea mays) is coinfected with a defensive symbiont (Fusarium verticillioides) and compared these fitness estimates to those obtained when the symbiont is absent. Here, virulence is measured as the reduction in the growth of the host caused by pathogen infection. Results of aster statistical models demonstrate that the defensive symbiont most negatively affects pathogen infection and that these effects propagate through subsequent stages of disease development to cause lower pathogen fitness across all virulence levels. Moreover, the virulence level at which pathogen fitness is maximal is higher in the presence of the defensive symbiont than in its absence. Thus, as expected from theory for multiple parasites, competition from the defensive symbiont may cause selection for increased pathogen virulence. More broadly, we consider that the evolutionary impacts of interactions between pathogens and microbial symbionts will depend critically on biological context and environment and that interactions among diverse microbial symbionts in spatially heterogeneous communities contribute to the maintenance of the highly diverse symbiotic functions observed in these communities.

An integrated analysis of phenotypic selection on insect body size and development time.

Most studies of phenotypic selection do not estimate selection or fitness surfaces for multiple components of fitness within a unified statistical framework. This makes it difficult or impossible to assess how selection operates on traits through variation in multiple components of fitness. We describe a new generation of aster models that can evaluate phenotypic selection by accounting for timing of life-history transitions and their effect on population growth rate, in addition to survival and reproductive output. We use this approach to estimate selection on body size and development time for a field population of the herbivorous insect, Manduca sexta (Lepidoptera: Sphingidae). Estimated fitness surfaces revealed strong and significant directional selection favoring both larger adult size (via effects on egg counts) and more rapid rates of early larval development (via effects on larval survival). Incorporating the timing of reproduction and its influence on population growth rate into the analysis resulted in larger values for size in early larval development at which fitness is maximized, and weaker selection on size in early larval development. These results illustrate how the interplay of different components of fitness can influence selection on size and development time. This integrated modeling framework can be readily applied to studies of phenotypic selection via multiple fitness components in other systems.

Inferring fitness landscapes.

Since 1983, study of natural selection has relied heavily on multiple regression of fitness on the values for a set of traits via ordinary least squares (OLSs), as proposed by Lande and Arnold, to obtain an estimate of the quadratic relationship between fitness and the traits, the fitness surface. However, well-known statistical problems with this approach can affect inferences about selection. One key concern is that measures of lifetime fitness do not conform to a normal or any other standard sampling distribution, as needed to justify the usual statistical tests. Another is that OLS may yield an estimate of the sign of the fitness function's curvature that is opposite to the truth. We here show that the recently developed aster modeling approach, which explicitly models the components of fitness as the basis for inferences about lifetime fitness, eliminates these problems. We illustrate selection analysis via aster using simulated datasets involving five fitness components expressed in each of four years. We demonstrate that aster analysis yields accurate estimates of the fitness function in cases in which OLS misleads, as well as accurate confidence regions for directional selection gradients. Further, to evaluate selection when many traits are under consideration, we recommend model selection by information criteria and frequentist model averaging.

Unifying life-history analyses for inference of fitness and population growth.

The lifetime fitnesses of individuals comprising a population determine its numerical dynamics, and genetic variation in fitness results in evolutionary change. This dual importance of individual fitness is well understood, but empirical fitness records generally violate the assumptions of standard statistical approaches. This problem has undermined comprehensive study of fitness and impeded empirical synthesis of the numerical and genetic dynamics of populations. Recently developed aster models remedy this problem by explicitly modeling the dependence of later-expressed components of fitness (e.g., fecundity) on those expressed earlier (e.g., survival to reproduce). Moreover, aster models employ different sampling distributions for different components of fitness (e.g., binomial for survival over a given interval and Poisson for fecundity). Analysis is done by maximum likelihood, and the resulting distributions for lifetime fitness closely approximate observed data. We illustrate the breadth of aster models' utility with three examples demonstrating estimation of the finite rate of increase, comparison of mean fitness among genotypic groups, and analysis of phenotypic selection. Aster models offer a unified approach to addressing the breadth of questions in evolution and ecology for which life-history data are gathered.

A comprehensive model of mutations affecting fitness and inferences for Arabidopsis thaliana.

As the ultimate source of genetic variation, spontaneous mutation is essential to evolutionary change. Theoretical studies over several decades have revealed the dependence of evolutionary consequences of mutation on specific mutational properties, including genomic mutation rates, U, and the effects of newly arising mutations on individual fitness, s. The recent resurgence of empirical effort to infer these properties for diverse organisms has not achieved consensus. Estimates, which have been obtained by methods that assume mutations are unidirectional in their effects on fitness, are imprecise. Both because a general approach must allow for occurrence of fitness-enhancing mutations, even if these are rare, and because recent evidence demands it, we present a new method for inferring mutational parameters. For the distribution of mutational effects, we retain Keightley's assumption of the gamma distribution, to take advantage of the flexibility of its shape. Because the conventional gamma is one sided, restricting it to unidirectional effects, we include an additional parameter, rho, as an amount it is displaced from zero. Estimation is accomplished by Markov chain Monte Carlo maximum likelihood. Through a limited set of simulations, we verify the accuracy of this approach. We apply it to analyze data on two reproductive fitness components from a 17-generation mutation-accumulation study of a Columbia accession of Arabidopsis thaliana in which 40 lines sampled in three generations were assayed simultaneously. For these traits, U approximately/= 0.1-0.2, with distributions of mutational effects broadly spanning zero, such that roughly half the mutations reduce reproductive fitness. One evolutionary consequence of these results is lower extinction risks of small populations of A. thaliana than expected from the process of mutational meltdown. A comprehensive view of the evolutionary consequences of mutation will depend on quantitatively accounting for fitness-enhancing, as well as fitness-reducing, mutations.