Unique pattern of <Emphasis Type="Italic">R-</Emphasis>gene variation within populations in Arabidopsis

An understanding of the variation pattern in disease resistance (R) genes is essential for its use in breeding programs aimed at neutralizing the threat of pathogens. Although the variation between populations is well known, there is little research about R-gene variation patterns within populations. Here, we investigate the polymorphism at three R-gene loci of 39 individual plants from nine populations of Arabidopsis thaliana. Our data suggest that alleles of each locus from individuals within a local population were either nearly identical, or highly diverse as ones between populations. The vast majority (92.5%) of within-population variation was shared globally, with high levels of allelic diversity (up to 11.7%) and abundant diverse-alleles. This unique pattern of within-population variation at R-loci suggests that individual plants within a population had the great potential to maintain a high level of globally-shared polymorphisms, and that the diversifying selection was the major force maintaining such polymorphisms. Consequently, the shared-polymorphism became recyclable for new R-genes, as the corresponding avirulence re-emerges in pathogen populations. Nucleotide sequence data reported are available in the GenBank databases under the accession numbers EF368054-EF368158.     

Comparing bivariate reaction norms among species: time and size at metamorphosis in three species of Hyla (Anura: Hylidae)

Summary Univariate and bivariate methods for comparing norms of reaction among species are discussed and illustrated with an example using North American hylid treefrogs. Norms of reaction for size at metamorphosis (SM) and length of larval period (LP) were compared among treefrog species raised at different food levels (Hyla cinerea vs H. gratiosa) and at different temperatures (H. cinerea vs H. gratiosa vs H. squirella). Hyla cinerea and Hyla gratiosa show parallel norms of reaction across food levels and temperatures. Across temperatures, H. squirella shows a much smaller change in SM relative to change in LP than do H. cinerea and H. gratiosa. This difference in shape of reaction norms may reflect different histories of selection resulting from these species' use of different larval habitats.     

SEX‐SPECIFIC PATTERNS OF MORPHOLOGICAL DIVERSIFICATION: EVOLUTION OF REACTION NORMS AND STATIC ALLOMETRIES IN NERIID FLIES

The consequences of sex‐specific selection for patterns of diversification remain poorly known. Because male secondary sexual traits are typically costly to express, and both costs and benefits are likely to depend on ambient environment and individual condition, such traits may be expected to diversify via changes in reaction norms as well as the scaling of trait size with body size (static allometry). We investigated morphological diversification within two species of Australian neriid flies (Telostylinus angusticollis, Telostylinus lineolatus) by rearing larvae from several populations on larval diets varying sixfold in nutrient concentration. Mean body size varied among populations of T. angusticollis, but body size reaction norms did not vary within either species. However, we detected diversification of reaction norms for body shape in males and females within both species. Moreover, unlike females, males also diversified in static allometry slope and reaction norms for static allometry slope of sexual and nonsexual traits. Our findings reveal qualitative sex differences in patterns of morphological diversification, whereby shape–size relationships diversify extensively in males, but remain conserved in females despite extensive evolution of trait means. Our results highlight the importance of incorporating plasticity and allometry in studies of adaptation and diversification.     

Patterns of divergence of a large family of nodule cysteine‐rich peptides in accessions of Medicago truncatula

The nodule cysteine‐rich (NCR) groups of defensin‐like (DEFL) genes are one of the largest gene families expressed in the nodules of some legume plants. They have only been observed in the inverted repeat loss clade (IRLC) of legumes, which includes the model legume Medicago truncatula. NCRs are reported to play an important role in plant–microbe interactions. To understand their diversity we analyzed their expression and sequence polymorphisms among four accessions of M. truncatula. A significant expression and nucleotide variation was observed among the genes. We then used 26 accessions to estimate the selection pressures shaping evolution among the accessions by calculating the nucleotide diversity at non‐synonymous and synonymous sites in the coding region. The mature peptides of the orthologous NCRs had signatures of both purifying and diversifying selection pressures, unlike the seed DEFLs, which predominantly exhibited purifying selection. The expression, sequence variation and apparent diversifying selection in NCRs within the Medicago species indicates rapid and recent evolution, and suggests that this family of genes is actively evolving to adapt to different environments and is acquiring new functions.     

Tertiary Windowing to Detect Positive Diversifying Selection

As a protein-encoding gene evolves, different selective pressures act on the gene temporally and spatially. An examination of the ratio of nonsynonymous-to-synonymous nucleotide substitution rate ratios (K_a/K_s) has proven to be a valuable method to examine selective pressures on protein encoding genes, including detecting positive diversifying selection. To gain power over averaging all sites in a gene together, examination of sites in primary sequence windows has frequently been employed. However, selection acts on folded proteins and sites that are close in tertiary space may not be close in primary sequence. A new method for the examination of K_a/K_s ratios based upon windows in tertiary structure is introduced and applied to the leptin gene family in mammals. Tertiary sequence windowing detects new sites under positive diversifying selection and detects positive diversifying selection with a more significant signal along various branches of the leptin gene family tree.Reviewing Editor: Dr. Rasmus Nielsen     

The evolution of tolerance to damage in Gentianella campestris: natural selection and the quantitative genetics of tolerance

In the framework of phenotypic plasticity, tolerance to browsing can be operationally defined as a norm of reaction comparing plant performance in undamaged and damaged conditions. Genetic variation in tolerance is then indicated by heterogeneity in the slopes of norms of reaction from a population. We investigated field gentian (Gentianella campestris) tolerance to damage in the framework of phenotypic plasticity using a sample of maternal lines from natural populations grown under common garden conditions and randomly split into either a control or an artificial clipping treatment. We found a diversity of tolerance norms of reaction at both the population and family level: the impacts of clipping ranged from poor tolerance (negative slope) to overcompensation (positive slope). We detected heterogeneity in tolerance norms of reaction in four populations. Similarly, we found a variety of plastic architectural responses to clipping and genetic variation in these responses in several populations. Overall, we found that the most tolerant populations were late flowering and also exhibit the greatest plastic increases in node (meristem) production in response to damage. We studied damage-imposed natural selection on plasticity in plant architecture in 10 of the sampled populations. In general, there was strong positive direct selection on final number of nodes for both control and clipped plants. However, the total selection on nodes (direct + indirect selection) within each treatment category depended heavily on the frequency of damage and cross-treatment genetic correlations in node production. In some cases, strong correlated responses to selection across the damage treatment led to total selection against nodes in the more rare environment. This could ultimately lead to the evolution of maladaptive phenotypes in one or both of the treatment categories. These results suggest that tolerance and a variety of architectural responses to damage may evolve by both direct and indirect responses to natural selection. While the present study demonstrates the potential importance of cross-treatment genetic correlations in directing the evolution of tolerance traits, such as branch or node production, we did not find any strong evidence of genetic trade-offs in candidate tolerance traits between undamaged and damaged conditions. This revised version was published online in July 2006 with corrections to the Cover Date.     

The evolution of reaction norms: simple models for age and size at maturity

This paper presents a simple model for the evolution of reaction norms for age and size at maturity that predicts reaction norms with a variety of shapes. Using realistic parameter values the model predicts reaction norms close to those observed in Drosophila. The major assumptions of the model are: 1) that net reproductive rate is maximized, 2) that growth is determinate, and 3) that mortality rates are independent of age and size at maturity. If, additionally, juvenile mortality is uncorrelated with a growth coefficient, k, the model predicts that selection favors maturation later at a smaller size when k is reduced by environmental factors and that decreased juvenile mortality leads to delayed maturity. These two predictions conform with those found by previous models using other measures of fitness. Correlations between k and juvenile mortality can change the shape of the predicted reaction norm. Depending on the precise form of the correlation, the model can predict done- or bowl-shaped reaction norms and can predict delayed or earlier maturity as k decreases. These shapes are qualitatively different from those predicted by previous models that used different fitness measures. Systematic estimates of the parameter values for this and for related models are required to determine the appropriate fitness measure for models of reaction norms.     

THE CONTRIBUTION OF NEW MUTATIONS TO GENOTYPE-ENVIRONMENT INTERACTION FOR FITNESS IN DROSOPHILA MELANOGASTER

Many studies have documented the existence of genotype-environment interaction (GEI) for traits closely related to fitness in natural populations. A type of GEI that is commonly observed is changes in the fitness ranking of genetic groups (families, clones, or inbred lines) in different environments. We refer to such changes in ranking as crossing of reaction norms for fitness. A common interpretation of crossing of reaction norms for fitness is that selection favors different alleles in the different environments (i.e., that “trade-offs” exist). If this is the case, selection could maintain genetic variation, and even lead to reproductive isolation between subpopulations using different environments. Even if the same alleles are favored in every environment, however, deleterious mutations that vary in the magnitude of their effect depending on environment could cause reaction norms for fitness to cross. If deleterious mutations with environment-dependent effects are responsible for maintaining much of the variation leading to crossing of reaction norms for fitness in natural populations, it should be possible to observe crossing of reaction norms for fitness among otherwise genetically identical lines bearing newly arisen spontaneous mutations. We examined the contribution of new mutations to GEI for fitness in Drosophila melanogaster. Eighteen lines were derived from a common, highly inbred base stock, and maintained at a population size of 10 pairs for over 200 generations, to allow them to accumulate spontaneous mutations. Because of the small population size of the lines, selection against mildly deleterious mutations should have been relatively ineffective. The lines were tested for productivity (number of surviving adult progeny from a standard number of parents) in five different environmental treatments, comprising different food media, temperatures, and levels of competition. The lines showed highly significant GEI for productivity, owing largely to considerable changes in ranking in the different environments. We conclude that mutations that are deleterious on average, but whose quantitative effects depend on environment, could be responsible for maintaining much of the variation leading to crossing of reaction norms for fitness that has been observed in samples of D. melanogaster from the wild.     

Seasonal migration promoting assortative mating in Littorina brevicula on a boulder shore in Japan

Littorina brevicula Philippi is one of the most common snails found in the upper intertidal zone of Japan. In Amakusa, some of the population of L. brevicula migrate to the lower zone in the winter, while the rest stay in the upper zone. Thus, during the winter, which is its reproductive season, the population of L. brevicula divides into two sub-populations. This leads to a hypothesis that the migration pattern in winter is genetically controlled and this behavioural dimorphism is maintained by reproductive isolation between the two sub-populations. In order to test this hypothesis, the following three points were investigated: (1) whether the same snails migrate in a similar way every winter, (2) whether there is a significant tidal level preference in snails, and (3) whether reproductive isolation occurs between the two sub-populations. The results showed (1) the migration behaviour of each snail was consistent over two successive winters, i.e. the same group of snails migrated downward every winter and the same group of snails stayed in the upper zone every winter, (2) transplanted snails moved toward the original zones where they were caught, suggesting that the snails actively selected their tidal zone in winter, and (3) most of the snails copulated within each sub-population. Therefore, reproductive isolation between the two sub-populations was considered to be established to some extent by the dimorphic migration behaviour. In conclusion, the migratory behaviour of L. brevicula is determined separately for each individual and might be genetically controlled, and the behavioural dimorphism may be maintained by partial reproductive isolation between the two sub-populations.     

Latitudinal variation in thermal reaction norms of post‐winter pupal development in two butterflies differing in phenological specialization

Latitudinal clines in thermal reaction norms of development are a common phenomenon in temperate insects. Populations from higher latitudes often develop faster throughout the range of relevant temperatures (i.e countergradient variation) because they must be able to complete their life cycle within a shorter seasonal time window compared to populations at lower latitudes. In the present study, we experimentally demonstrate that two species of butterflies Anthocharis cardamines (L.) and Pieris napi (L.) instead show a cogradient variation in thermal reaction norms of post‐winter pupal development so that lower latitude populations develop faster than higher latitude populations. The two species share host plants but differ in the degree of phenological specialization, as well as in the patterns of voltinism. We suggest that the pattern in A. cardamines, a univoltine phenological specialist feeding exclusively on flowers and seedpods, is the result of selection for matching to the phenological pattern of its local host plants. The other species, P. napi, is a phenological generalist feeding on the leaves of the hosts and it shows a latitudinal cline in voltinism. Because the latitudinal pattern in P. napi was an effect of slow development in a fraction of the pupae from the most northern population, we hypothesize that this population may include both bivoltine and univoltine genotypes. Consequently, although the two species both showed cogradient patterns in thermal reaction norms, it appears likely that this was for different reasons. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 113 , 981–991.     

Contrasting patterns of evolution between allelic groups at a single locus in <Emphasis Type="Italic">Arabidopsis</Emphasis>

Heterogeneities in evolutionary pattern among different loci are commonly observed. To see whether the heterogeneity can also be observed among allelic groups in a single locus, we investigated the coding sequence and the flanking regions of Rpp13, a disease resistance gene in up to 60 accession lines from worldwide populations in Arabidopsis thaliana. An extraordinarily high level of polymorphism (π=0.098) and four distinct clades were found in the leucine-rich repeat (LRR) region in this gene. No obvious geographic relationship with the clades was observed, and such clades were not observed in the other regions in and around this gene. The average genetic diversity among the clades ranged from 10 to 14.6% in the LRR. The levels of polymorphism within each clade varied largely, and significant heterogeneity in evolutionary rates among clades was detected. A statistically significant departure from neutrality was also detected by Fu & Li’s tests. These results suggest that both directional and diversifying selection are working on this locus, and that natural selection can cause heterogeneity in evolutionary rate, even among allele groups in a locus. Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

Inversion and allozyme polymorphism show contrasting patterns of microgeographical population structure in a natural population of Drosophila buzzatii from Argentina

Second chromosome inversions and genotypic frequencies at seven allozyme loci were determined in a natural population of the cactophilic species Drosophila buzzatii that uses as breeding sites the necrotic cladodes of the prickly pear Opuntia quimilo and the rotting stems of cardón, Trichocereus terschekii. Different processes govern the evolutionary fate of inversion and allozyme polymorphisms. A pattern of heterotic balance for inversions seems to be acting uniformly in each breeding site and could depend on different regimes of density‐dependent selection within cactus hosts. Patterns of variation of allozymes revealed significant heterogeneity in allele frequencies for Esterase‐1 (Est‐1) among O. quimilo rots and Aldehyde oxidase (Aldox) and Xanthine dehydrogenase (Xdh) among T. terschekii substrates and showed gene‐cactus effects only for Esterase‐2 (Est‐2). Consistent and significant excesses of homozygotes were detected at both the within‐rot and in the total population levels that could be accounted for by diversifying selection among individual breeding sites.     

VARIATION IN REACTION NORMS AMONG POPULATIONS OF THE GRASS BOUTELOUA RIGIDISETA

Recent research has emphasized the importance of investigating the reaction norms of quantitative traits to understand evolution in natural environments. In this study, genetic differences in reaction norms among eight populations of the grass Bouteloua rigidiseta were examined using clonal replicates of genotypes planted in a common garden with two levels of competition (single B. rigidiseta without competition and single B. rigidiseta surrounded by four Erioneuron pilosum). The populations were found to be genetically differentiated for a variety of traits. Differences in reaction norms of size-specific fecundity (spikelet clusters per tiller number) were detected among the populations: some showed little response to competition; in others size-specific fecundity was much greater in the absence of competition. This divergence in reaction norms among these populations may be the result of past selection (including the cost of plasticity), or genetic drift.     

Gene duplications,divergence and recombination shape adaptive evolution of the fish ectoparasite Gyrodactylus bullatarudis

Determining the molecular basis of parasite adaptation to its host is an important component in understanding host–parasite coevolution and the epidemiology of parasitic infections. Here, we investigate short‐ and long‐term adaptive evolution in the eukaryotic parasite Gyrodactylus bullatarudis infecting Caribbean guppies (Poecilia reticulata), by comparing the reference genome of Tobagonian G. bullatarudis with other Platyhelminthes, and by analysing resequenced samples from local Trinidadian populations. At the macroevolutionary timescale, we observed duplication of G‐protein and serine proteases genes, which are probably important in host–parasite arms races. Serine protease also showed strong evidence of ongoing, diversifying selection at the microevolutionary timescale. Furthermore, our analyses revealed that a hybridization event, involving two divergent genomes, followed by recombination has dramatically affected the genetic composition of Trinidadian populations. The recombinant genotypes invaded Trinidad and replaced local parasites in all populations. We localized more than 300 genes in regions fixed in local populations for variants of different origin, possibly due to diversifying selection pressure from local host populations. In addition, around 70 genes were localized in regions identified as heterozygous in some, but not all, individuals. This pattern is consistent with a very recent spread of recombinant parasites. Overall, our results are consistent with the idea that recombination between divergent genomes can result in particularly successful parasites.     

A comparative assessment of purification techniques for mesophyll protoplasts of Brassica napus L.

The effects of three different general purification protocols have been assessed quantitatively using mesophyll protoplasts of Brassica napus. Within the initial sample two distinct sub-populations were determined. The methods used influenced the ratio of the vacuolated to chloroplastic type protoplast sub-populations. Overall recovery rates of the initial sample varied according to the method used from 38% to 27%, but the relative recovery of the sub-populations varied considerably with a purified ratio of between 1.0:0.78 to 1.0:7.0. Size distribution profiles of the initial and purified populations are also presented.     

The Interaction of Protein Structure, Selection, and Recombination on the Evolution of the Type-1 Fimbrial Major Subunit (fimA) from Escherichia coli

Fimbrial adhesins allow bacteria to interact with and attach to their environment. The bacteria possibly benefit from these interactions, but all external structures including adhesins also allow bacteria to be identified by other organisms. Thus adhesion molecules might be under multiple forms of selection including selection to constrain functional interactions or evolve novel epitopes to avoid recognition. We address these issues by studying genetic diversity in the Escherichia coli type-1 fimbrial major subunit, fimA. Overall, sequence diversity in fimA is high (π= 0.07) relative to that in other E. coli genes. High diversity is a function of positive diversifying selection, as detected by d _N/d _S ratios higher than 1.0, and amino acid residuces subject to diversifying selection are nonrandomly clustered on the exterior surface of the peptide. In addition, McDonald and Kreitman tests suggest that there has been historical but not current directional selection at fimA between E. coli and Salmonella. Finally, some regions of the fimA peptide appear to be under strong structural constraint within E. coli, particularly the interior regions of the molecule that is involved in subunit to subunit interaction. Recombination also plays a major role contributing to E. coli fimA allelic variation and estimates of recombination (2N _e c) and mutation (2N _eμ) are about the same. Recombination may act to separate the diverse evolutionary forces in different regions of the fimA peptide. Received: 13 April 2000 / Accepted: 28 October 2000     

Rapid Evolution by Positive Selection and Gene Gain and Loss: PLA<Subscript>2</Subscript> Venom Genes in Closely Related <Emphasis Type="Italic">Sistrurus</Emphasis> Rattlesnakes with Divergent Diets

Rapid evolution of snake venom genes by positive selection has been reported previously but key features of this process such as the targets of selection, rates of gene turnover, and functional diversity of toxins generated remain unclear. This is especially true for closely related species with divergent diets. We describe the evolution of PLA₂ gene sequences isolated from genomic DNA from four taxa of Sistrurus rattlesnakes which feed on different prey. We identified four to seven distinct PLA₂ sequences in each taxon and phylogenetic analyses suggest that these sequences represent a rapidly evolving gene family consisting of both paralogous and homologous loci with high rates of gene gain and loss. Strong positive selection was implicated as a driving force in the evolution of these protein coding sequences. Exons coding for amino acids that make up mature proteins have levels of variation two to three times greater than those of the surrounding noncoding intronic sequences. Maximum likelihood models of coding sequence evolution reveal that a high proportion (∼30%) of all codons in the mature protein fall into a class of codons with an estimated d _N /d _S (ω) ratio of at least 2.8. An analysis of selection on individual codons identified nine residues as being under strong (p < 0.01) positive selection, with a disproportionately high proportion of these residues found in two functional regions of the PLA₂ protein (surface residues and putative anticoagulant region). This is direct evidence that diversifying selection has led to high levels of functional diversity due to structural differences in proteins among these snakes. Overall, our results demonstrate that both gene gain and loss and protein sequence evolution via positive selection are important evolutionary forces driving adaptive divergence in venom proteins in closely related species of venomous snakes.     

Adaptive molecular evolution of MC1R gene reveals the evidence for positive diversifying selection in indigenous goat populations

Detecting signatures of selection can provide a new insight into the mechanism of contemporary breeding and artificial selection and further reveal the causal genes associated to the phenotypic variation. However, the signatures of selection on genes entailing for profitable traits between Chinese commercial and indigenous goats have been poorly interpreted. We noticed footprints of positive selection at MC 1R gene containing SNP s genotyped in five Chinese native goat breeds. An experimental distribution of F _ST was built based on approximations of F _ST for each SNP across five breeds. We identified selection using the high F _ST outlier method and found that MC 1R candidate gene show evidence of positive selection. Furthermore, adaptive selection pressure on specific codons was determined using different codon based on maximum‐likelihood methods; signature of positive selection in mammalian MC 1R was explored in individual codons. Evolutionary analyses were inferred under maximum likelihood models, the HyPhy package implemented in the DATAMONKEY Web Server. The results of codon selection displayed positive diversifying selection at the sites were mainly involved in development of genetic variations in coat color in various mammalian species. Positive diversifying selection inferred with recent evolutionary changes in domesticated goat MC 1R provides new insights that the gene evolution may have been modulated by domestication events in goats.     

The effect of small-scale environmental changes on disease incidence and severity in a natural plant-pathogen interaction

Summary The incidence and severity of Rhynchosporium secalis infections were assessed in a large population of Hordeum leporinum. Transects were set out in four directions from five trees to determine the effect of shading. Under the tree canopy 60.3% of H. leporinum plants were infected while only 11.2% were infected away from the canopy. Disease severity, on those plants which were infected, was higher under the canopy (mean 12.4% and 13.0% leaf area diseased for the flag and first leaves, respectively) than away from the canopy (means of 7.8% and 5.0% for the flag and first leaves respectively). Plants under the tree canopy contained on average 23% more nitrogen, raising the possibility that the susceptibility of the host changed in response to nitrogen levels. However, the observed pattern is also consistent with the hypothesis that shade-associated changes in the environment enhanced the ability of the pathogen to infect and develop on the host. The data clearly demonstrate the importance of small-scale environmental factors on natural host-pathogen interactions. These environmental factors may cause differential selection for disease resistance within a host population, which may ultimately lead to the formation of sub-populations with differing levels of resistance.