Norms of reaction and diversifying selection

The numbers of progeny produced by comparable numbers of female Drosophila melanogaster of 26 geographic strains on nine different culture media are examined in the context of norms of reaction. Having emphasized that diversifying selection is seldom discussed simultaneously with its seemingly related topic, norms of reaction, I present the following argument: diversifying selection has generally been viewed as involving sub-populations inhabiting separate localities and subject to different patterns of selection, norms of reaction as variation whose weighted average determines the relative fitnesses of different genotypes within individual sub-populations. Should environmental challenges frequently involve life or death (including sterility) outcomes, norms of reaction involving components of fitness engender diversifying selection within local populations (demes).     

Geographical variation in selection, from phenotypes to molecules

Molecular technologies now allow researchers to isolate quantitative trait loci (QTLs) and measure patterns of gene sequence variation within chromosomal regions containing important polymorphisms. I develop a simulation model to investigate gene sequence evolution within genomic regions that harbor QTLs. The QTLs influence a trait experiencing geographical variation in selection, which is common in nature and produces obvious differentiation at the phenotypic level. Counter to expectations, the simulations suggest that selection can substantially affect quantitative genetic variation without altering the amount and pattern of molecular variation at sites closely linked to the QTLs. Even with large samples of gene sequences, the likelihood of rejecting neutrality is often low. The exception is situations where strong selection is combined with low migration among demes, conditions that may be common in many plant species. The results have implications for gene sequence surveys and, perhaps more generally, for interpreting the apparently weak connection between levels of molecular and quantitative trait variation within species.     

The stolbur phytoplasma antigenic membrane protein gene stamp is submitted to diversifying positive selection

Surface proteins play an important role in phytoplasma life cycle. The antigenic membrane protein (AMP) of 'Candidatus Phytoplasma asteris' has been shown to interact with the insect microfilaments. Due to the synteny of the groL-amp-nadE genetic locus between phytoplasma genomes, the gene stamp that encodes the antigenic membrane protein of stolbur phytoplasma has been cloned and characterized. It encodes a 157 amino acid-long protein with a predicted signal peptide and a C-terminal hydrophobic alpha-helix. STAMP was 26-40% identical to AMP of 'Candidatus Phytoplasma asteris' strains and 40% identical to AMP of 'Ca. P. japonicum'. The expression of STAMP in Escherichia. coli produced a 16 kDa peptide recognized by an anti-stolbur monoclonal antibody. Stamp was more variable than the house-keeping gene secY and the ratio between non-synonymous over synonymous mutations (dN/dS ) was 2.78 for stamp as compared to 0.64 for secY. This indicates that stamp is submitted to a positive diversifying selection pressure.     

Positive diversifying selection in avian Mx genes

Mx proteins are interferon-induced GTPases that confer antiviral activities against RNA viruses. We analysed the molecular evolution of the Mx gene in birds using data on interspecific divergence in anseriform and galliform birds, and on intraspecific diversity in commercial chicken lines, local Chinese chicken breeds as well as in the mallard. The overall ratio of non-synonymous to synonymous substitution was unusually high, 0.80, indicating relaxed constraint or positive selection. Evidence for the latter was provided by that a total of 11–18 codons were found to have evolved under positive selection. The great majority of these codons are located in a region unique to birds at the N-terminal end of the Mx protein. We found an excess of non-synonymous polymorphisms relative to synonymous variants in all comparisons. This, together with positive Tajima’s D values in the local Chinese chicken breeds and in the mallard suggests that balancing selection is acting in avian Mx genes. As such, Mx mimics the major histocompatibility complex system, indicating that heterozygous individuals are better off withstanding pathogen attack. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. S.B. and L.Q. contributed equally to this work.     

Estimating diversifying selection and functional constraint in the presence of recombination

Models of molecular evolution that incorporate the ratio of nonsynonymous to synonymous polymorphism (dN/dS ratio) as a parameter can be used to identify sites that are under diversifying selection or functional constraint in a sample of gene sequences. However, when there has been recombination in the evolutionary history of the sequences, reconstructing a single phylogenetic tree is not appropriate, and inference based on a single tree can give misleading results. In the presence of high levels of recombination, the identification of sites experiencing diversifying selection can suffer from a false-positive rate as high as 90%. We present a model that uses a population genetics approximation to the coalescent with recombination and use reversible-jump MCMC to perform Bayesian inference on both the dN/dS ratio and the recombination rate, allowing each to vary along the sequence. We demonstrate that the method has the power to detect variation in the dN/dS ratio and the recombination rate and does not suffer from a high false-positive rate. We use the method to analyze the porB gene of Neisseria meningitidis and verify the inferences using prior sensitivity analysis and model criticism techniques.     

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.     

Genome-wide signatures of population bottlenecks and diversifying selection in European wolves

Genomic resources developed for domesticated species provide powerful tools for studying the evolutionary history of their wild relatives. Here we use 61K single-nucleotide polymorphisms (SNPs) evenly spaced throughout the canine nuclear genome to analyse evolutionary relationships among the three largest European populations of grey wolves in comparison with other populations worldwide, and investigate genome-wide effects of demographic bottlenecks and signatures of selection. European wolves have a discontinuous range, with large and connected populations in Eastern Europe and relatively smaller, isolated populations in Italy and the Iberian Peninsula. Our results suggest a continuous decline in wolf numbers in Europe since the Late Pleistocene, and long-term isolation and bottlenecks in the Italian and Iberian populations following their divergence from the Eastern European population. The Italian and Iberian populations have low genetic variability and high linkage disequilibrium, but relatively few autozygous segments across the genome. This last characteristic clearly distinguishes them from populations that underwent recent drastic demographic declines or founder events, and implies long-term bottlenecks in these two populations. Although genetic drift due to spatial isolation and bottlenecks seems to be a major evolutionary force diversifying the European populations, we detected 35 loci that are putatively under diversifying selection. Two of these loci flank the canine platelet-derived growth factor gene, which affects bone growth and may influence differences in body size between wolf populations. This study demonstrates the power of population genomics for identifying genetic signals of demographic bottlenecks and detecting signatures of directional selection in bottlenecked populations, despite their low background variability.     

The signature of positive selection on standing genetic variation

Considerable interest is focused on the use of polymorphism data to identify regions of the genome that underlie recent adaptations. These searches are guided by a simple model of positive selection, in which a mutation is favored as soon as it arises. This assumption may not be realistic, as environmental changes and range expansions may lead previously neutral or deleterious alleles to become beneficial. We examine what effect this mode of selection has on patterns of variation at linked neutral sites by implementing a new coalescent model of positive directional selection on standing variation. In this model, a neutral allele arises and drifts in the population, then at frequency f becomes beneficial, and eventually reaches fixation. Depending on the value of f, this scenario can lead to a large variance in allele frequency spectra and in levels of linkage disequilibrium at linked, neutral sites. In particular, for intermediate f, the beneficial substitution often leads to a loss of rare alleles--a pattern that differs markedly from the signature of directional selection currently relied on by researchers. These findings highlight the importance of an accurate characterization of the effects of positive selection, if we are to reliably identify recent adaptations from polymorphism data.     

Purification and structure of the host-specific toxin from Helminthosporium carbonum race 1

The host-specific toxin from Helminthosporiumcarbonum race 1 was purified from culture filtrates by solvent extraction, gel filtration, and high pressure liquid chromatography. High resolution mass spectrometry of the purified toxin gave a MW of 436.2318 and an elemental composition C₂₁H₃₂N₄O₆. Amino acid analysis and proton and¹³C-NMR indicated a peptide containing four amino acids. Their sequence was determined by gas chromatography mass spectrometry. Finally, digestion of the amino acids with D- and L-amino acid oxidases gave the complete structure cyclo[(L-2-amino-9, 10-epoxy-8-oxodecanoyl)-D-prolyl-L-alanyl-L-alanyl].     

Rapid detection of positive selection in genes and genomes through variation clusters

Positive selection in genes and genomes can point to the evolutionary basis for differences among species and among races within a species. The detection of positive selection can also help identify functionally important protein regions and thus guide protein engineering. Many existing tests for positive selection are excessively conservative, vulnerable to artifacts caused by demographic population history, or computationally very intensive. I here propose a simple and rapid test that is complementary to existing tests and that overcomes some of these problems. It relies on the null hypothesis that neutrally evolving DNA regions should show a Poisson distribution of nucleotide substitutions. The test detects significant deviations from this expectation in the form of variation clusters, highly localized groups of amino acid changes in a coding region. In applying this test to several thousand human-chimpanzee gene orthologs, I show that such variation clusters are not generally caused by relaxed selection. They occur in well-defined domains of a protein's tertiary structure and show a large excess of amino acid replacement over silent substitutions. I also identify multiple new human-chimpanzee orthologs subject to positive selection, among them genes that are involved in reproductive functions, immune defense, and the nervous system.     

Comparative genomics and diversifying selection of the clustered vertebrate protocadherin genes

To explain the mechanism for specifying diverse neuronal connections in the brain, Sperry proposed that individual cells carry chemoaffinity tags on their surfaces. The enormous complexity of these connections requires a tremendous diversity of cell-surface proteins. A large number of neural transmembrane protocadherin (Pcdh) proteins is encoded by three closely linked human and mouse gene clusters (alpha, beta, and gamma). To gain insight into Pcdh evolution, I performed comprehensive comparative cDNA and genomic DNA analyses for the three clusters in the chimpanzee, rat, and zebrafish genomes. I found that there are species-specific duplications in vertebrate Pcdh genes and that additional diversity is generated through alternative splicing within the zebrafish "variable" and "constant" regions. Moreover, different codons (sites) in the mammalian Pcdh ectodomains (ECs) are under diversifying selection, with some under diversity-enhancing positive Darwinian selection and others, including calcium-binding sites, under strong purifying selection. Interestingly, almost all positively selected codon positions are located on the surface of ECs 2 and 3. These diversified residues likely play an important role in combinatorial interactions of Pcdh proteins, which could provide the staggering diversity required for neuronal connections in the brain. These results also suggest that adaptive selection is an additional evolutionary factor for increasing Pcdh diversity.     

Preservation of a pseudogene by gene conversion and diversifying selection

Interlocus gene conversion is considered a crucial mechanism for generating novel combinations of polymorphisms in duplicated genes. The importance of gene conversion between duplicated genes has been recognized in the major histocompatibility complex and self-incompatibility genes, which are likely subject to diversifying selection. To theoretically understand the potential role of gene conversion in such situations, forward simulations are performed in various two-locus models. The results show that gene conversion could significantly increase the number of haplotypes when diversifying selection works on both loci. We find that the tract length of gene conversion is an important factor to determine the efficacy of gene conversion: shorter tract lengths can more effectively generate novel haplotypes given the gene conversion rate per site is the same. Similar results are also obtained when one of the duplicated genes is assumed to be a pseudogene. It is suggested that a duplicated gene, even after being silenced, will contribute to increasing the variability in the other locus through gene conversion. Consequently, the fixation probability and longevity of duplicated genes increase under the presence of gene conversion. On the basis of these findings, we propose a new scenario for the preservation of a duplicated gene: when the original donor gene is under diversifying selection, a duplicated copy can be preserved by gene conversion even after it is pseudogenized.     

Evidence for diversifying selection at the pyoverdine locus of Pseudomonas aeruginosa

Pyoverdine is the primary siderophore of the gram-negative bacterium Pseudomonas aeruginosa. The pyoverdine region was recently identified as the most divergent locus alignable between strains in the P. aeruginosa genome. Here we report the nucleotide sequence and analysis of more than 50 kb in the pyoverdine region from nine strains of P. aeruginosa. There are three divergent sequence types in the pyoverdine region, which correspond to the three structural types of pyoverdine. The pyoverdine outer membrane receptor fpvA may be driving diversity at the locus: it is the most divergent alignable gene in the region, is the only gene that showed substantial intratype variation that did not appear to be generated by recombination, and shows evidence of positive selection. The hypothetical membrane protein PA2403 also shows evidence of positive selection; residues on one side of the membrane after protein folding are under positive selection. R', previously identified as a type IV strain, is clearly derived from a type III strain via a 3.4-kb deletion which removes one amino acid from the pyoverdine side chain peptide. This deletion represents a natural modification of the product of a nonribosomal peptide synthetase enzyme, whose consequences are predictive from the DNA sequence. There is also linkage disequilibrium between the pyoverdine region and pvdY, a pyoverdine gene separated by 30 kb from the pyoverdine region. The pyoverdine region shows evidence of horizontal transfer; we propose that some alleles in the region were introduced from other soil bacteria and have been subsequently maintained by diversifying selection.     

A host-specific factor is necessary for efficient folding of the autotransporter plasmid-encoded toxin

Autotransporters are the most common virulence factors secreted from Gram-negative pathogens. Until recently, autotransporter folding and outer membrane translocation were thought to be self-mediated events that did not require accessory factors. Here, we report that two variants of the autotransporter plasmid-encoded toxin are secreted by a lab strain of Escherichia coli. Biophysical analysis and cell-based toxicity assays demonstrated that only one of the two variants was in a folded, active conformation. The misfolded variant was not produced by a pathogenic strain of enteroaggregative E. coli and did not result from protein overproduction in the lab strain of E. coli. Our data suggest a host-specific factor is required for efficient folding of plasmid-encoded toxin.     

Calreticulin molecular evolution: a strong purifying and episodic diversifying selection result

Calreticulin (CRT) is a low molecular weight protein present in vertebrates, invertebrates and higher plants. Its multiple functions have been demonstrated. It plays an important role as a chaperone and Ca²⁺ buffer inside sarcoplasmic/endoplasmic reticulum (SR/ER), and outside the ER in many physiological/pathological processes. Recently it has been observed that CRT over-expression or its absence is linked to various pathological conditions, such as malignant evolution and progression, and these facts really increased its study interests. Using an evolution approach CRT was further characterized. Several Bayesian phylogenetic analyses were performed using coding and amino acid sequences. CRT molecular evolution was investigated for the presence of negative or/and positive selection using HyPhy package. The results indicated that the purifying selection might have operated over the whole CRT primary structure. Although, an episodic diversifying selection was also found on the analyzed CRT sequences.     

Patterns of diversifying selection in the phytotoxin-like scr74 gene family of Phytophthora infestans

Phytophthora infestans, the organism responsible for the Irish famine, causes late blight, a re-emerging disease of potato and tomato. Little is known about the molecular evolution of P. infestans genes. To identify candidate effector genes (virulence or avirulence genes) that may have co-evolved with the host, we mined expressed sequence tag (EST) data from infection stages of P. infestans for secreted and potentially polymorphic genes. This led to the identification of scr74, a gene that encodes a predicted 74-amino acid secreted cysteine-rich protein with similarity to the Phytophthora cactorum phytotoxin PcF. The expression of scr74 was upregulated approximately 60-fold 2 to 4 days after inoculation of tomato and was also significantly induced during early stages of colonization of potato. The scr74 gene was found to belong to a highly polymorphic gene family within P. infestans with 21 different sequences identified. Using the approximate and maximum likelihood (ML) methods, we found that diversifying selection likely caused the extensive polymorphism observed within the scr74 gene family. Pairwise comparisons of 17 scr74 sequences revealed elevated ratios of nonsynonymous to synonymous nucleotide-substitution rates, particularly in the mature region of the proteins. Using ML, all 21 polymorphic amino acid sites were identified to be under diversifying selection. Of these 21 amino acids, 19 are located in the mature protein region, suggesting that selection may have acted on the functional portions of the proteins. Further investigation of gene copy number and organization revealed that the scr74 gene family comprises at least three copies located in a region of no more than 300 kb of the P. infestans genome. We found evidence that recombination contributed to sequence divergence within at least one gene locus. These results led us to propose an evolutionary model that involves gene duplication and recombination, followed by functional divergence of scr74 genes. This study provides support for using diversifying selection as a criterion for identifying candidate effector genes from sequence databases.