Balancing selection at the tomato RCR3 Guardee gene family maintains variation in strength of pathogen defense

Coevolution between hosts and pathogens is thought to occur between interacting molecules of both species. This results in the maintenance of genetic diversity at pathogen antigens (or so-called effectors) and host resistance genes such as the major histocompatibility complex (MHC) in mammals or resistance (R) genes in plants. In plant-pathogen interactions, the current paradigm posits that a specific defense response is activated upon recognition of pathogen effectors via interaction with their corresponding R proteins. According to the "Guard-Hypothesis," R proteins (the "guards") can sense modification of target molecules in the host (the "guardees") by pathogen effectors and subsequently trigger the defense response. Multiple studies have reported high genetic diversity at R genes maintained by balancing selection. In contrast, little is known about the evolutionary mechanisms shaping the guardee, which may be subject to contrasting evolutionary forces. Here we show that the evolution of the guardee RCR3 is characterized by gene duplication, frequent gene conversion, and balancing selection in the wild tomato species Solanum peruvianum. Investigating the functional characteristics of 54 natural variants through in vitro and in planta assays, we detected differences in recognition of the pathogen effector through interaction with the guardee, as well as substantial variation in the strength of the defense response. This variation is maintained by balancing selection at each copy of the RCR3 gene. Our analyses pinpoint three amino acid polymorphisms with key functional consequences for the coevolution between the guardee (RCR3) and its guard (Cf-2). We conclude that, in addition to coevolution at the "guardee-effector" interface for pathogen recognition, natural selection acts on the "guard-guardee" interface. Guardee evolution may be governed by a counterbalance between improved activation in the presence and prevention of auto-immune responses in the absence of the corresponding pathogen.     

Evolutionary Analyses Suggest a Function of MxB Immunity Proteins Beyond Lentivirus Restriction

Viruses impose diverse and dynamic challenges on host defenses. Diversifying selection of codons and gene copy number variation are two hallmarks of genetic innovation in antiviral genes engaged in host-virus genetic conflicts. The myxovirus resistance (Mx) genes encode interferon-inducible GTPases that constitute a major arm of the cell-autonomous defense against viral infection. Unlike the broad antiviral activity of MxA, primate MxB was recently shown to specifically inhibit lentiviruses including HIV-1. We carried out detailed evolutionary analyses to investigate whether genetic conflict with lentiviruses has shaped MxB evolution in primates. We found strong evidence for diversifying selection in the MxB N-terminal tail, which contains molecular determinants of MxB anti-lentivirus specificity. However, we found no overlap between previously-mapped residues that dictate lentiviral restriction and those that have evolved under diversifying selection. Instead, our findings are consistent with MxB having a long-standing and important role in the interferon response to viral infection against a broader range of pathogens than is currently appreciated. Despite its critical role in host innate immunity, we also uncovered multiple functional losses of MxB during mammalian evolution, either by pseudogenization or by gene conversion from MxA genes. Thus, although the majority of mammalian genomes encode two Mx genes, this apparent stasis masks the dramatic effects that recombination and diversifying selection have played in shaping the evolutionary history of Mx genes. Discrepancies between our study and previous publications highlight the need to account for recombination in analyses of positive selection, as well as the importance of using sequence datasets with appropriate depth of divergence. Our study also illustrates that evolutionary analyses of antiviral gene families are critical towards understanding molecular principles that govern host-virus interactions and species-specific susceptibility to viral infection.     

Rapid genome‐wide evolution in Brassica rapa populations following drought revealed by sequencing of ancestral and descendant gene pools

There is increasing evidence that evolution can occur rapidly in response to selection. Recent advances in sequencing suggest the possibility of documenting genetic changes as they occur in populations, thus uncovering the genetic basis of evolution, particularly if samples are available from both before and after selection. Here, we had a unique opportunity to directly assess genetic changes in natural populations following an evolutionary response to a fluctuation in climate. We analysed genome‐wide differences between ancestors and descendants of natural populations of Brassica rapa plants from two locations that rapidly evolved changes in multiple phenotypic traits, including flowering time, following a multiyear late‐season drought in California. These ancestor‐descendant comparisons revealed evolutionary shifts in allele frequencies in many genes. Some genes showing evolutionary shifts have functions related to drought stress and flowering time, consistent with an adaptive response to selection. Loci differentiated between ancestors and descendants (F_ST outliers) were generally different from those showing signatures of selection based on site frequency spectrum analysis (Tajima's D), indicating that the loci that evolved in response to the recent drought and those under historical selection were generally distinct. Very few genes showed similar evolutionary responses between two geographically distinct populations, suggesting independent genetic trajectories of evolution yielding parallel phenotypic changes. The results show that selection can result in rapid genome‐wide evolutionary shifts in allele frequencies in natural populations, and highlight the usefulness of combining resurrection experiments in natural populations with genomics for studying the genetic basis of adaptive evolution.     

Selection analysis on the rapid evolution of a secondary sexual trait

Evolutionary analyses of population translocations (experimental or accidental) have been important in demonstrating speed of evolution because they subject organisms to abrupt environmental changes that create an episode of selection. However, the strength of selection in such studies is rarely measured, limiting our understanding of the evolutionary process. This contrasts with long-term, mark–recapture studies of unmanipulated populations that measure selection directly, yet rarely reveal evolutionary change. Here, we present a study of experimental evolution of male colour in Trinidadian guppies where we tracked both evolutionary change and individual-based measures of selection. Guppies were translocated from a predator-rich to a low-predation environment within the same stream system. We used a combination of common garden experiments and monthly sampling of individuals to measure the phenotypic and genetic divergence of male coloration between ancestral and derived fish. Results show rapid evolutionary increases in orange coloration in both populations (1 year or three generations), replicating the results of previous studies. Unlike previous studies, we linked this evolution to an individual-based analysis of selection. By quantifying individual reproductive success and survival, we show, for the first time, that males with more orange and black pigment have higher reproductive success, but males with more black pigment also have higher risk of mortality. The net effect of selection is thus an advantage of orange but not black coloration, as reflected in the evolutionary response. This highlights the importance of considering all components of fitness when understanding the evolution of sexually selected traits in the wild.     

The genetic basis of traits regulating sperm competition and polyandry: can selection favour the evolution of good- and sexy-sperm?

The good-sperm and sexy-sperm (GS-SS) hypotheses predict that female multiple mating (polyandry) can fuel sexual selection for heritable male traits that promote success in sperm competition. A major prediction generated by these models, therefore, is that polyandry will benefit females indirectly via their sons' enhanced fertilization success. Furthermore, like classic 'good genes' and 'sexy son' models for the evolution of female preferences, GS-SS processes predict a genetic correlation between genes for female mating frequency (analogous to the female preference) and those for traits influencing fertilization success (the sexually selected traits). We examine the premise for these predictions by exploring the genetic basis of traits thought to influence fertilization success and female mating frequency. We also highlight recent debates that stress the possible genetic constraints to evolution of traits influencing fertilization success via GS-SS processes, including sex-linked inheritance, nonadditive effects, interacting parental genotypes, and trade-offs between integrated ejaculate components. Despite these possible constraints, the available data suggest that male traits involved in sperm competition typically exhibit substantial additive genetic variance and rapid evolutionary responses to selection. Nevertheless, the limited data on the genetic variation in female mating frequency implicate strong genetic maternal effects, including X-linkage, which is inconsistent with GS-SS processes. Although the relative paucity of studies on the genetic basis of polyandry does not allow us to draw firm conclusions about the evolutionary origins of this trait, the emerging pattern of sex linkage in genes for polyandry is more consistent with an evolutionary history of antagonistic selection over mating frequency. We advocate further development of GS-SS theory to take account of the complex evolutionary dynamics imposed by sexual conflict over mating frequency.     

Genetic diversity and the evolutionary history of plant immunity genes in two species of Zea

Plant pathogenesis-related genes (PR genes) code for enzymes, enzyme inhibitors, and other peptides that confer resistance to pathogens and herbivores. Although several PR genes have been the subject of molecular population genetic analyses, a general understanding of their long-term evolutionary dynamics remains incomplete. Here we analyze sequence data from 17 PR genes from two closely related teosinte species of central Mexico. In addition to testing whether patterns of diversity at individual loci depart from expectations under a neutral model, we compared patterns of diversity at defense genes, as a class, to nondefense genes. In Zea diploperennis, the majority of defense genes have patterns of diversity consistent with neutral expectations while at least two genes showed evidence of recent positive selection consistent with arms-race models of antagonistic coevolution. In Zea mays ssp. parviglumis, by contrast, analyses of both defense and nondefense genes revealed strong and consistent departures from the neutral model, suggestive of nonequilibrium population dynamics or population structure. Nevertheless, we found a significant excess of replacement polymorphism in defense genes compared to nondefense genes. Although we cannot exclude relaxed selective constraint as an explanation, our results are consistent with temporally variable (transient or episodic) selection or geographically variable selection acting on parviglumis defense genes. The different patterns of diversity found in the two Zea species may be explained by parviglumis' greater distribution and population structure together with geographic variation in selection.     

Predicting evolutionary responses to selection on polyandry in the wild: additive genetic covariances with female extra-pair reproduction

The evolutionary forces that underlie polyandry, including extra-pair reproduction (EPR) by socially monogamous females, remain unclear. Selection on EPR and resulting evolution have rarely been explicitly estimated or predicted in wild populations, and evolutionary predictions are vulnerable to bias due to environmental covariances and correlated selection through unmeasured traits. However, evolutionary responses to (correlated) selection on any trait can be directly predicted as additive genetic covariances (cov_A) with appropriate components of relative fitness. I used comprehensive life-history, paternity and pedigree data from song sparrows (Melospiza melodia) to estimate cov_A between a female''s liability to produce extra-pair offspring and two specific fitness components: relative annual reproductive success (ARS) and survival to recruitment. All three traits showed non-zero additive genetic variance. Estimates of cov_A were positive, predicting evolution towards increased EPR, but 95% credible intervals overlapped zero. There was therefore no conclusive prediction of evolutionary change in EPR due to (correlated) selection through female ARS or recruitment. Negative environmental covariance between EPR and ARS would have impeded evolutionary prediction from phenotypic selection differentials. These analyses demonstrate an explicit quantitative genetic approach to predicting evolutionary responses to components of (correlated) selection on EPR that should be unbiased by environmental covariances and unmeasured traits.     

An information‐theoretic approach to estimating the composite genetic effects contributing to variation among generation means: Moving beyond the joint‐scaling test for line cross analysis

The pace and direction of evolution in response to selection, drift, and mutation are governed by the genetic architecture that underlies trait variation. Consequently, much of evolutionary theory is predicated on assumptions about whether genes can be considered to act in isolation, or in the context of their genetic background. Evolutionary biologists have disagreed, sometimes heatedly, over which assumptions best describe evolution in nature. Methods for estimating genetic architectures that favor simpler (i.e., additive) models contribute to this debate. Here we address one important source of bias, model selection in line cross analysis (LCA). LCA estimates genetic parameters conditional on the best model chosen from a vast model space using relatively few line means. Current LCA approaches often favor simple models and ignore uncertainty in model choice. To address these issues we introduce Software for Analysis of Genetic Architecture (SAGA), which comprehensively assesses the potential model space, quantifies model selection uncertainty, and uses model weighted averaging to accurately estimate composite genetic effects. Using simulated data and previously published LCA studies, we demonstrate the utility of SAGA to more accurately define the components of complex genetic architectures, and show that traditional approaches have underestimated the importance of epistasis.     

Genetic constraints and sexual dimorphism in immune defense

The absence of continued evolutionary change despite the presence of genetic variation and directional selection is very common. Genetic correlations between traits can reduce the evolvability of traits. One intriguing example might be found in a sexual conflict over sexually dimorphic traits: a common genetic architecture constrains the response to selection on a trait subjected to sexually asymmetric selection pressures. Here we show that males and females of the mealworm beetle Tenebrio molitor differ in the quantitative genetic architecture of four traits related to immune defense and condition. Moreover, high genetic correlations between the sexes constitute a genetic constraint to the evolution of sexual dimorphism in immune defense. Our results suggest a general mechanism by which sexual conflict can promote evolutionary stasis. We furthermore show negative genetic correlations, strong indications of trade-offs, between immune traits for two pairs of traits in females.     

The probability of parallel genetic evolution from standing genetic variation

Parallel evolution is often assumed to result from repeated adaptation to novel, yet ecologically similar, environments. Here, we develop and analyse a mathematical model that predicts the probability of parallel genetic evolution from standing genetic variation as a function of the strength of phenotypic selection and constraints imposed by genetic architecture. Our results show that the probability of parallel genetic evolution increases with the strength of natural selection and effective population size and is particularly likely to occur for genes with large phenotypic effects. Building on these results, we develop a Bayesian framework for estimating the strength of parallel phenotypic selection from genetic data. Using extensive individual‐based simulations, we show that our estimator is robust across a wide range of genetic and evolutionary scenarios and provides a useful tool for rigorously testing the hypothesis that parallel genetic evolution is the result of adaptive evolution. An important result that emerges from our analyses is that existing studies of parallel genetic evolution frequently rely on data that is insufficient for distinguishing between adaptive evolution and neutral evolution driven by random genetic drift. Overcoming this challenge will require sampling more populations and the inclusion of larger numbers of loci.     

PLEIOTROPY CAUSES LONG-TERM GENETIC CONSTRAINTS ON LIFE-HISTORY EVOLUTION IN BRASSICA RAPA

Fundamental, long-term genetic trade-offs constrain life-history evolution in wild crucifer populations. I studied patterns of genetic constraint in Brassica rapa by estimating genetic correlations among life-history components by quantitative genetic analyses among ten wild populations, and within four populations. Genetic correlations between age and size at first reproduction were always greater than +0.8 within and among all populations studied. Although quantitative genetics might provide insight about genetic constraints if genetic parameters remain approximately constant, little evidence has been available to determine the constancy of genetic correlations. I found strong and consistent estimates of genetic correlations between life-history components, which were very similar within four natural populations. Population differentiation also showed these same trade-offs, resulting from long-term genetic constraint. For some traits, evolutionary changes among populations were incompatible with a model of genetic drift. Historical patterns of natural selection were inferred from population differentiation, suggesting that correlated response to selection has caused some traits to evolve opposite to the direct forces of natural selection. Comparison with Arabidopsis suggests that these life-history trade-offs are caused by genes that regulate patterns of resource allocation to different components of fitness. Ecological and energetic models may correctly predict these trade-offs because there is little additive genetic variation for rates of resource acquisition, but resource allocation is genetically variable.     

Population genetics of duplicated disease-defense genes,hm1 and hm2, in maize (Zea mays ssp. mays L.) and its wild ancestor (Zea mays ssp. parviglumis)

Plant defense genes are subject to nonneutral evolutionary dynamics. Here we investigate the evolutionary dynamics of the duplicated defense genes hm1 and hm2 in maize and its wild ancestor Zea mays ssp. parviglumis. Both genes have been shown to confer resistance to the fungal pathogen Cochliobolus carbonum race 1, but the effectiveness of resistance differs between loci. The genes also display different population histories. The hm1 locus has the highest nucleotide diversity of any gene yet sampled in the wild ancestor of maize, and it contains a large number of indel polymorphisms. There is no evidence, however, that high diversity in hm1 is a product of nonneutral evolution. In contrast, hm2 has very low nucleotide diversity in the wild ancestor of maize. The distribution of hm2 polymorphic sites is consistent with nonneutral evolution, as indicated by Tajima's D and other neutrality tests. In addition, one hm2 haplotype is more frequent than expected under the equilibrium neutral model, suggesting hitchhiking selection. Both defense genes retain >80% of the level of genetic variation in maize relative to the wild ancestor, and this level is similar to other maize genes that were not subject to artificial selection during domestication.     

Unstable DNA Repair Genes Shaped by Their Own Sequence Modifying Phenotypes

The question of whether natural selection favors genetic stability or genetic variability is a fundamental problem in evolutionary biology. Bioinformatic analyses demonstrate that selection favors genetic stability by avoiding unstable nucleotide sequences in protein encoding DNA. Yet, such unstable sequences are maintained in several DNA repair genes, thereby promoting breakdown of repair and destabilizing the genome. Several studies have therefore argued that selection favors genetic variability at the expense of stability. Here we propose a new evolutionary mechanism, with supporting bioinformatic evidence, that resolves this paradox. Combining the concepts of gene-dependent mutation biases and meiotic recombination, we argue that unstable sequences in the DNA mismatch repair (MMR) genes are maintained by their own phenotype. In particular, we predict that human MMR maintains an overrepresentation of mononucleotide repeats (monorepeats) within and around the MMR genes. In support of this hypothesis, we report a 31% excess in monorepeats in 250 kb regions surrounding the seven MMR genes compared to all other RefSeq genes (1.75 vs. 1.34%, P = 0.0047), with a particularly high content in PMS2 (2.41%, P = 0.0047) and MSH6 (2.07%, P = 0.043). Based on a mathematical model of monorepeat frequency, we argue that the proposed mechanism may suffice to explain the observed excess of repeats around MMR genes. Our findings thus indicate that unstable sequences in MMR genes are maintained through evolution by the MMR mechanism. The evolutionary paradox of genetically unstable DNA repair genes may thus be explained by an equilibrium in which the phenotype acts back on its own genotype.     

EXPERIMENTAL EVOLUTION OF RESISTANCE IN BRASSICA RAPA: CORRELATED RESPONSE OF TOLERANCE IN LINES SELECTED FOR GLUCOSINOLATE CONTENT

The evolutionary response of plant populations to selection for increased defense may be constrained by costs of defense. The purpose of this study was to investigate such constraints on the evolution of defense due to a cost of defense manifested as a trade-off between defense and tolerance. Variation in the response to artificial damage (tolerance) among lines of Brassica rapa that had been artificially selected for foliar glucosinolate content (defense) was examined. Leaf area was removed from replicates of three selection lines (high glucosinolates, control, and low glucosinolates) at three damage levels (0%, 20%, and 60% damage). An external cost of defense would result in a statistically significant selection line by damage treatment interaction, with those selected for high defense expressing less tolerance than those selected for low defense. Damage treatment had a significant overall effect on estimated total fitness, with fitness declining with increasing damage level. Further, selection line also had a significant overall effect on estimated total fitness, with low-defense selection lines having higher fitness compared to both control and high-defense selection lines. More importantly, a cost of defense in terms of tolerance was demonstrated by a significant selection line-by-damage treatment interaction. This interaction was in the direction to demonstrate a genetic trade-off between defense and tolerance, with low-defense selection lines decreasing estimated total fitness in response to damage less than both control and high-defense selection lines. Variation in tolerance among selection lines was due to the greater ability of low-defense lines to maintain fruit and seed production despite the presence of damage. In terms of tolerance, this cost of glucosinolate production in B. rapa could constrain the evolution of increased defense and, in so doing, maintain individuals within the population that are poorly defended yet tolerant.     

The molecular population genetics of regulatory genes

Regulatory loci, which may encode both trans acting proteins as well as cis acting promoter regions, are crucial components of an organism's genetic architecture. Although evolution of these regulatory loci is believed to underlie the evolution of numerous adaptive traits, there is little information on natural variation of these genes. Recent molecular population genetic studies, however, have provided insights into the extent of natural variation at regulatory genes, the evolutionary forces that shape them and the phenotypic effects of molecular regulatory variants. These recent analyses suggest that it may be possible to study the molecular evolutionary ecology of regulatory diversification by examining both the extent and patterning of regulatory gene diversity, the phenotypic effects of molecular variation at these loci and their ecological consequences.     

On indirect genetic effects in structured populations

Indirect genetic effects (IGEs) occur when the phenotype of an individual, and possibly its fitness, depends, at least in part, on the genes of its social partners. The effective result is that environmental sources of phenotypic variance can themselves evolve. Simple models have shown that IGEs can alter the rate and direction of evolution for traits involved in interactions. Here we expand the applicability of the theory of IGEs to evolution in metapopulations by including nonlinear interactions between individuals and population genetic structure. Although population subdivision alone generates some dramatic and nonintuitive evolutionary dynamics for interacting phenotypes, the combination of nonlinear interactions with subdivision reveals an even greater importance of IGEs. The presence of genetic structure links the evolution of interacting phenotypes and the traits that influence their expression ("effector traits") even in the absence of genetic correlations. When nonlinear social effects occur in subdivided populations, evolutionary response is altered and can even oppose the direction expected due to direct selection. Because population genetic structure allows for multilevel selection, we also investigate the role of IGEs in determining the response to individual and group selection. We find that nonlinear social effects can cause interference between levels of selection even when they act in the same direction. In some cases, interference can be so extreme that the actual evolutionary response to multilevel selection is opposite in direction to that predicted by summing selection at each level. This theoretical result confirms empirical data that show higher levels of selection cannot be ignored even when selection acts in the same direction at all levels.     

Detecting lineage-specific adaptive evolution of brain-expressed genes in human using rhesus macaque as outgroup

Comparative genetic analysis between human and chimpanzee may detect genetic divergences responsible for human-specific characteristics. Previous studies have identified a series of genes that potentially underwent Darwinian positive selection during human evolution. However, without a closely related species as outgroup, it is difficult to identify human-lineage-specific changes, which is critical in delineating the biological uniqueness of humans. In this study, we conducted phylogeny-based analyses of 2633 human brain-expressed genes using rhesus macaque as the outgroup. We identified 47 candidate genes showing strong evidence of positive selection in the human lineage. Genes with maximal expression in the brain showed a higher evolutionary rate in human than in chimpanzee. We observed that many immune-defense-related genes were under strong positive selection, and this trend was more prominent in chimpanzee than in human. We also demonstrated that rhesus macaque performed much better than mouse as an outgroup in identifying lineage-specific selection in humans.