Hidden evolution: progress and limitations in detecting multifarious natural selection

From illustrative examples of research on the best-studied group of species to date, Drosophila melanogaster and its closest relatives, we argue that selection is multifarious, but often hidden. Selective fixation of new, highly advantageous alleles is the most parsimonious explanation for a typical pattern of molecular variation observed in genomic regions characterized by very low recombination: drastically reduced DNA sequence variation within species and typical levels of sequence divergence among species. At the same time, the identity of the gene (or genes) influenced by selection is not just difficult to discern; it may be impossible. Studies of the genetic basis of reproductive isolation demonstrate that, although the D. melanogaster complex species appear virtually identical, dozens of currently unidentified genes contribute to hybrid sterility. We argue that these findings are best explained by selectively-driven functional divergence and demonstrate the multifarious nature of selection. Although multifarious selection certainly occurs, the exact characters responsible for differences in survival and reproductive success are unknown. We do not see these inherent limits as a cause for despair or a problem for evolutionary biology. Instead, we hope to raise awareness of these complexities of evolution by highlighting both the progress and the limitations of characterizing multifarious natural selection.     

Recombination Detection Under Evolutionary Scenarios Relevant to Functional Divergence

Recombination can negatively impact methods designed to detect divergent gene function that rely on explicit knowledge of a gene tree. However, we know little about how recombination detection methods perform under evolutionary scenarios encountered in studies of functional molecular divergence. We use simulation to evaluate false positive rates for six recombination detection methods (GENECONV, MaxChi, Chimera, RDP, GARD-SBP, GARD-MBP) under evolutionary scenarios that might increase false positives. Broadly, these scenarios address: (i) asymmetric tree topology and sequence divergence, (ii) non-stationary codon bias and selection pressure, and (iii) positive selection. We also evaluate power to detect recombination under truly recombinant history. As with previous studies, we find that power increases with sequence divergence. However, we also find that accuracy to correctly infer the number of breakpoints is extremely low. When recombination is absent, increased sequence divergence leads to increased false positives. Furthermore, one method (GARD-SBP) is sensitive to tree shape, with higher false positive rates under an asymmetric tree topology. Somewhat surprisingly, all methods are robust to the simulated heterogeneity in codon bias, shifts in selection pressure and presence of positive selection. Based on these findings, we recommend that studies of functional divergence in systems where recombination is plausible can, and should, include a pre-test for recombination. Application of all methods to the core genome of Prochlorococcus reveals a substantial lack of concordance among results. Based on analysis of both real and simulated datasets we present some guidelines for the investigation of recombination in genes that may have experienced functional divergence.     

Evolution of RXLR-Class Effectors in the Oomycete Plant Pathogen Phytophthora ramorum

Phytophthora plant pathogens contain many hundreds of effectors potentially involved in infection of host plants. Comparative genomic analyses have shown that these effectors evolve rapidly and have been subject to recent expansions. We examined the recent sequence evolution of RXLR-class effector gene families in the sudden oak death pathogen, P. ramorum. We found that P. ramorum RXLR effectors have taken multiple evolutionary paths, including loss or gain of repeated domains, recombination or gene conversion among paralogs, and selection on point mutations. Sequencing of homologs from two subfamilies in P. ramorum’s closest known relatives revealed repeated gene duplication and divergence since speciation with P. lateralis. One family showed strong signatures of recombination while the other family has evolved primarily by point mutation. Comparison of a small number of the hundreds of RXLR-class effectors across three clonal lineages of P. ramorum shows striking divergence in alleles among lineages, suggesting the potential for functional differences between lineages. Our results suggest future avenues for examination of rapidly evolving effectors in P. ramorum, including investigation of the functional and coevolutionary significance of the patterns of sequence evolution that we observed.     

The influence of recombination on SNP diversity in chickens

The association between recombination rate and diversity, but not divergence is considered to be driven mainly by natural selection: fixation of positively selected variants and associated hitchhiking effects and/or background selection eliminating deleterious alleles. In the present study, we investigated the relationship between recombination rate, SNP diversity and interspecies divergence for 29 loci in chickens. We found that recombination rate is positively correlated with nucleotide diversity but is not correlated with interspecies divergence. It appears that variation in recombination rate explains over 30% of the variation in levels of diversity among 29 loci. Our data suggested that natural selection is a main factor in shaping SNP diversity in chickens. Since SNP diversity is significantly lower at Z-linked than at autosomal loci, we argued that genetic hitchhiking might be more important than background selection in producing the observed correlation.     

Evolutionary history of an MHC gene in two leporid species: characterisation of <Emphasis Type="Italic">Mhc-DQA</Emphasis> in the European brown hare and comparison with the European rabbit

We surveyed the genetic diversity of the expressed major histocompatibility complex class II DQA locus in natural populations of European brown hares, Lepus europaeus, from Austria and Belgium (267 individuals in total). Based on cDNA sequences, we designed hare-specific primers to amplify the highly variable second exon of the DQA gene. Using cloning–sequencing methodology and capillary electrophoresis single-strand conformation polymorphism, we found ten alleles of the DQA exon 2 locus across these two European regions, of which eight are described for the first time. To search for signals of selection and recombination in the evolution of the DQA gene within the leporids, we augmented our sample with orthologous DQA alleles from the European rabbit, Oryctolagus cuniculus, in order to carry out a species level, species pairwise comparison. We found evidence of recombination in the history of the DQA sequences in leporids with some recombinant alleles bridging the species divide. In both species, selection on peptide binding site codons can be detected, though stronger for the rabbit. This result suggests that there may be a differential selection pressure in the deeper evolutionary history of these two species due to differences in several demographic and ecological traits likely subjecting them to differential selection by parasites. Finally, evolutionary relationships show a widespread and statistically significant intermingling of alleles from the two species. The many macroparasites shared between hares and rabbits may explain this pattern of trans-species polymorphism. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. The nucleotide sequence data reported in this paper have been submitted to Genbank and have been assigned the accession numbers FJ225335–FJ225346.     

Non-neutral evolution of the major histocompatibility complex class II gene DRB1 in the sac-winged bat Saccopteryx bilineata

The immune genes of the major histocompatibility complex (MHC) are classical examples for high levels of genetic diversity and non-neutral evolution. This is particularly true for the regions containing the antigen-binding sites as, for instance, in the exon 2 of the MHC class II gene DRB. We surveyed, for the first time in the order Chiroptera, the genetic diversity within this exon in the sac-winged bat Saccopteryx bilineata. We detected 11 alleles among 85 bats, of which 79 were sampled in one population. Pairwise comparisons revealed that interallelic sequence differences ranged between 3 and 22%, although nucleotide substitutions were not evenly distributed along the exon sequence. This was most probably the result of intragenic recombination. High levels of sequence divergence and significantly more nonsynonymous than synonymous substitutions (d(N)/d(S)>1) suggest long-term balancing selection. Thus, the data are consistent with the hypothesis that recombination gives rise to new alleles at the DRB locus of the sac-winged bat, and these are maintained in the population through balancing selection. In this respect, the sac-winged bat closely resembles other mammalian species.     

Evolution of the phosphoglycerate mutase processed gene in human and chimpanzee revealing the origin of a new primate gene

Processed genes are created by retroposition from messenger RNA of expressed genes. The estimated amount of processed copies of genes in the human genome is 10,000-14,000. Some of these might be pseudogenes with the expected pattern for nonfunctional sequences, but some others might be an important source of new genes. We have studied the evolution of a Phosphoglycerate mutase processed gene (PGAM3) described in humans and believed to be a pseudogene. We sequenced PGAM3 in chimpanzee and macaque and obtained polymorphism data for human coding region. We found evidence that PGAM3 likely produces a functional protein, as an example of addressing functionality for human processed pseudogenes. First, the open reading frame was intact despite many deletions that occurred in the 3' untranslated region. Second, it appears that the gene is expressed. Finally, interspecies and intraspecies variation for PGAM3 was not consistent with the neutral model proposed for pseudogenes, suggesting that a new functional primate gene has originated. Amino acid divergence was significantly higher than synonymous divergence in PGAM3 lineage, supporting positive selection acting in this gene. This role of selection was further supported by the excess of rare alleles in a population genetic analysis. PGAM3 is located in a region of very low recombination; therefore, it is conceivable that the rapid fixation events in this newly arising gene may have contributed to a selective sweep of variation in the region.     

Selection and recombination drive the evolution of MHC class II DRB diversity in ungulates

Major histocompatibility complex (MHC) antigen-presenting genes are the most variable loci in vertebrate genomes. Host-parasite co-evolution is assumed to maintain the excessive polymorphism in the MHC loci. However, the molecular mechanisms underlying the striking diversity in the MHC remain contentious. The extent to which recombination contributes to the diversity at MHC loci in natural populations is still controversial, and there have been only few comparative studies that make quantitative estimates of recombination rates. In this study, we performed a comparative analysis for 15 different ungulates species to estimate the population recombination rate, and to quantify levels of selection. As expected for all species, we observed signatures of strong positive selection, and identified individual residues experiencing selection that were congruent with those constituting the peptide-binding region of the human DRB gene. However, in addition for each species, we also observed recombination rates that were significantly different from zero on the basis of likelihood-permutation tests, and in other non-quantitative analyses. Patterns of synonymous and non-synonymous sequence diversity were consistent with differing demographic histories between species, but recent simulation studies by other authors suggest inference of selection and recombination is likely to be robust to such deviations from standard models. If high rates of recombination are common in MHC genes of other taxa, re-evaluation of many inference-based phylogenetic analyses of MHC loci, such as estimates of the divergence time of alleles and trans-specific polymorphism, may be required.     

Effect of Divergence Time and Recombination Rate on Molecular Evolution of <Emphasis Type="Italic">Drosophila</Emphasis> INE-1 Transposable Elements and Other Candidates for Neutrally Evolving Sites

Interspecies divergence of orthologous transposable element remnants is often assumed to be simply due to genetic drift of neutral mutations that occurred after the divergence of the species. However, divergence may also be affected by other factors, such as variation in the mutation rate, ancestral polymorphisms, or selection. Here we attempt to determine the impact of these forces on divergence of three classes of sites that are often assumed to be selectively unconstrained (INE-1 TE remnants, sites within short introns, and fourfold degenerate sites) in two different pairwise comparisons of Drosophila (D. melanogaster vs. D. simulans and D. simulans vs. D. sechellia). We find that divergence of these three classes of sites is strongly influenced by the recombination environment in which they are located, and this is especially true for the closer D. simulans vs. D. sechellia comparison. We suggest that this is mainly a result of the contribution of ancestral polymorphisms in different recombination regions. We also find that intergenic INE-1 elements are significantly more diverged than intronic INE-1 in both pairwise comparisons, implying the presence of either negative selection or lower mutation rates in introns. Furthermore, we show that substitution rates in INE-1 elements are not associated with the length of the noncoding sequence in which they are located, suggesting that reduced divergence in long noncoding sequences is not due to reduced mutation rates in these regions. Finally, we show that GC content for each site within INE-1 sequences has evolved toward an equilibrium value (approximately 33%) since insertion.     

Molecular Evolution at the Self-Incompatibility Locus of <Emphasis Type="BoldItalic">Physalis longifolia</Emphasis> (Solanaceae)

The simple codominant expression of the alleles at the self-incompatibility locus ( S) of Solanaceae and their extraordinary spectrum of sequence diversity present an ideal case for understanding the molecular forces that shape the sequence polymorphism of genes involved in recognition reactions. Here, with unprecedented conspecific sequences of 33 Physalis longifolia S alleles, including 27 of their single introns, the tempo-spatial nucleotide substitution patterns were first detailed for the S locus and analyzed in genealogical time scales. Three major genealogical clades of the 33 P. longifolia S alleles were used to divide the genealogical time into the within-clade and the between-clade periods. During the within-clade period, the average nonsynonymous substitution rate was 50% higher than the intron substitution rate but the opposite trend emerged for the between-clade period. A new and simple method developed here was utilized to estimate the selection intensity in the coding regions. The magnitudes and the distribution of these estimates, in conjunction with the spatial substitution pattern among closely related S sequences, revealed an initial short-term action of strong positive selection and a continuous but weaker action of negative selection brought by functional/structural constraints on the S alleles. The two modes of selection can significantly modify the branch lengths of an S genealogy and may be ubiquitous to recognition systems.     

Recombination During In Vitro Evolution

Recombination, the swapping of large portions of genetic information between and among parental genotypes, can be applied to in vitro evolution experiments on functional nucleic acids. Both homologous and heterologous recombination can be achieved using standard laboratory techniques. In many cases, recombination can allow for the discovery of a ribozyme or DNAzyme phenotype that would not likely be encountered by reliance on point mutations alone. In addition, recombination can often aid in the discovery of global optima in sequence space and/or lessen the number of generations it would take to reach optima. Recombination is most efficiently used in combination with point mutations and applied after the first couple of rounds of selection but before high-fitness genotypes dominate the selection. The “recombination zone” describes that region of sequence space—defined by the residues that will ultimately participate in the function of the winning nucleic acid(s)—where recombination is expected to be the most beneficial in the search for high-fitness genotypes.[Reviewing Editor: Martin Kreitman]Author order determined by a single Bernoulli trial as implemented by RPS.     

Balancing selection and its effects on sequences in nearby genome regions

Our understanding of balancing selection is currently becoming greatly clarified by new sequence data being gathered from genes in which polymorphisms are known to be maintained by selection. The data can be interpreted in conjunction with results from population genetics models that include recombination between selected sites and nearby neutral marker variants. This understanding is making possible tests for balancing selection using molecular evolutionary approaches. Such tests do not necessarily require knowledge of the functional types of the different alleles at a locus, but such information, as well as information about the geographic distribution of alleles and markers near the genes, can potentially help towards understanding what form of balancing selection is acting, and how long alleles have been maintained.     

Rapid divergence of two classes of Haemophilus ducreyi

Haemophilus ducreyi, the etiologic agent of chancroid, expresses variants of several key virulence factors. While previous reports suggested that H. ducreyi strains formed two clonal populations, the differences between, and diversity within, these populations were unclear. To assess their variability, we examined sequence diversity at 11 H. ducreyi loci, including virulence and housekeeping genes, augmenting published data sets with PCR-amplified genes to acquire data for at least 10 strains at each locus. While sequences from all 11 loci place strains into two distinct groups, there was very little variation within each group. The difference between alleles of the two groups was variable and large at 3 loci encoding surface-exposed proteins (0.4 < K(S) < 1.3, where K(S) is divergence at synonymous sites) but consistently small at genes encoding cytoplasmic or periplasmic proteins (K(S) < 0.09). The data suggest that the two classes have recently diverged, that recombination has introduced variant alleles into at least 3 distinct loci, and that these alleles have been confined to one of the two classes. In addition, recombination is evident among alleles within, but not between, classes. Rather than clones of the same species, these properties indicate that the two classes may form distinct species.     

Unique nucleotide polymorphism of ankyrin gene cluster in <Emphasis Type="Italic">Arabidopsis</Emphasis>

The ankyrin (ANK) gene cluster is a part of a multigene family encoding ANK transmembrane proteins in Arabidopsis thaliana, and plays an important role in protein-protein interactions and in signal pathways. In contrast to other regions of a genome, the ANK gene cluster exhibits an extremely high level of DNA polymorphism in an ∼5-kb region, without apparent decay. Phylogenetic analysis detects two clear, deeply differentiated haplotypes (dimorphism). The divergence between haplotypes of accession Col-0 and Ler-0 (Hap-C and Hap-L) is estimated to be 10.7%, approximately equal to the 10.5% average divergence between A. thaliana and A. lyrata. Sequence comparisons for the ANK gene cluster homologues in Col-0 indicate that the members evolve independently, and that the similarity among paralogues is lower than between alleles. Very little intralocus recombination or gene conversion is detected in ANK regions. All these characteristics of the ANK gene cluster are consistent with a tandem gene duplication and birth-and-death process. The possible mechanisms for and implications of this elevated nucleotide variation are also discussed, including the suggestion of balancing selection.     

Patterns of Historical Balancing Selection on the Salmonid Major Histocompatibility Complex Class II β Gene

Allelic variation in the major histocompatibility class (MHC) IIB gene of salmonids is analyzed for patterns indicative of natural selection acting at the molecular level. Sequence data for the second exon of this MHC gene were generated for 11 species in three salmonid genera: Oncorhynchus, Salmo, and Salvelinus. Phylogenetic analysis of nucleotide sequences revealed: (1) monophyletic grouping of alleles from each genus, (2) transspecies evolution of alleles within Salmo and Salvelinus, and (3) differential patterns of transspecies evolution within the genus Oncorhynchus. Within Oncorhynchus, five of seven species had alleles that were species-specific or nearly so, while the remaining two, O. mykiss and O. clarkii, retained ancestral polymorphisms. The different patterns in Oncorhynchus and the other two genera could be due to historical demographic effects or functional differences in MHC molecules in the three genera, but the two hypotheses could not be distinguished with the current dataset. An analysis of recombination/gene conversion identified numerous recombinant alleles, which is consistent with what has been found in other vertebrate taxa. However, these gene conversion events could not account for the species-specific allelic lineages observed in five of the Oncorhynchus species. Analyses of the relative rates of nonsynonymous and synonymous substitutions revealed the signature of selection on the class IIB gene in all 11 of the salmonid species for both the ABS and the non-ABS codons. Codon-based analyses of selection identified seven codons that have experienced selection in the majority of the species. More than half of these sites were mammalian ABS codons, but several were not, suggesting subtle functional differences in the mammalian and teleost fish MHC molecules.     

Likelihood Analysis of the Chalcone Synthase Genes Suggests the Role of Positive Selection in Morning Glories (<Emphasis Type="Italic">Ipomoea</Emphasis>)

Chalcone synthase (CHS) is a key enzyme in the biosynthesis of flavonoides, which are important for the pigmentation of flowers and act as attractants to pollinators. Genes encoding CHS constitute a multigene family in which the copy number varies among plant species and functional divergence appears to have occurred repeatedly. In morning glories (Ipomoea), five functional CHS genes (A–E) have been described. Phylogenetic analysis of the Ipomoea CHS gene family revealed that CHS A, B, and C experienced accelerated rates of amino acid substitution relative to CHS D and E. To examine whether the CHS genes of the morning glories underwent adaptive evolution, maximum-likelihood models of codon substitution were used to analyze the functional sequences in the Ipomoea CHS gene family. These models used the nonsynonymous/synonymous rate ratio ( = d_N/d_S) as an indicator of selective pressure and allowed the ratio to vary among lineages or sites. Likelihood ratio test suggested significant variation in selection pressure among amino acid sites, with a small proportion of them detected to be under positive selection along the branches ancestral to CHS A, B, and C. Positive Darwinian selection appears to have promoted the divergence of subfamily ABC and subfamily DE and is at least partially responsible for a rate increase following gene duplication.     

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.     

MHC diversity in bottlenecked populations: a simulation model

The depletion of variation at MHC loci, which play a crucial role in pathogen recognition, has been postulated to be one of important extinction risk factors for endangered populations. Thus, it is important to understand how selection affects the level of polymorphism in these genes when populations undergo a reduction in size. We followed MHC diversity in computer simulations of population bottlenecks. The fates of MHC alleles in the simulations were determined either by drift, or by balancing selection resulting from host–parasite coevolution. We found that the impact of selection on MHC polymorphism in bottlenecked populations was dependent upon the timescales involved. Initially, selection maintained lower number of alleles than drift, but after ~40 generations of hosts selection maintained higher MHC diversity, as compared to drift. The adverse effects of decreased MHC polymorphism on population viability may be, to some extent, compensated for if selection helps to retain MHC alleles which show high functional diversity, which should allow protection against a broader range of pathogens. Our simulation shows, however, that the mean divergence of alleles retained under selection in bottlenecked populations is not, on average, significantly higher than the divergence due to drift.