Applications of selective neutrality tests to molecular ecology

Conceptual bases for quantifying the role of the environment on gene evolution: the participation of positive selection and neutral evolution

To demonstrate that a given change in the environment has contributed to the emergence of a given genotypic and phenotypic shift during the course of evolution, one should ask to what extent such shifts would have occurred without environmental change. Of course, such tests are rarely practical but phenotypic novelties can still be correlated to genomic shifts in response to environmental changes if enough information is available. We surveyed and re-evaluated the published data in order to estimate the role of environmental changes on the course of species and genomic evolution. Only a few published examples clearly demonstrate a causal link between a given environmental change and the fixation of a genomic variant resulting in functional modification (gain, loss or alteration of function). Many others suggested a link between a given phenotypic shift and a given environmental change but failed to identify the underlying genomic determinant(s) and/or the associated functional consequence(s). The proportion of genotypic and phenotypic variation that is fixed concomitantly with environmental changes is often considered adaptive and hence, the result of positive selection, even though alternative causes, such as genetic drift, are rarely investigated. Therefore, the second aim herein is to review evidence for the mechanisms leading to fixation.

Evidence for positive Darwinian selection on the hepcidin gene of Perciform and Pleuronectiform fishes

Hepcidin is a small cysteine-rich peptide that plays an important role in antimicrobial activity and in maintaining iron homeostasis in vertebrates. Here we report on the underlying mechanism that maintains high sequence diversities among the hepcidin-like variants of perciform and pleuronectiform fishes. In contrast to mammals, maximum likelihood-based codon substitution analyses revealed that positive Darwinian selection (nonsynonymous to synonymous substitution, omega > 1) is the likely cause of accelerated rate of amino acid substitutions in the hepcidin mature peptide region of these fishes. Comparison of models incorporating positive selection (omega > 1) at certain sites with models not incorporating positive selection (omega < 1) failed to reject (p = 0) the evidence of positive selection among the codon sites of percifom and pleuronectiform hepcidin. The adaptive evolution of this peptide in perciform and pleuronectiform fishes might be directed by pathogens when the host is exposed to new habitats/environments.

A spatial analysis method (SAM) to detect candidate loci for selection: towards a landscape genomics approach to adaptation

The detection of adaptive loci in the genome is essential as it gives the possibility of understanding what proportion of a genome or which genes are being shaped by natural selection. Several statistical methods have been developed which make use of molecular data to reveal genomic regions under selection. In this paper, we propose an approach to address this issue from the environmental angle, in order to complement results obtained by population genetics. We introduce a new method to detect signatures of natural selection based on the application of spatial analysis, with the contribution of geographical information systems (GIS), environmental variables and molecular data. Multiple univariate logistic regressions were carried out to test for association between allelic frequencies at marker loci and environmental variables. This spatial analysis method (SAM) is similar to current population genomics approaches since it is designed to scan hundreds of markers to assess a putative association with hundreds of environmental variables. Here, by application to studies of pine weevils and breeds of sheep we demonstrate a strong correspondence between SAM results and those obtained using population genetics approaches. Statistical signals were found that associate loci with environmental parameters, and these loci behave atypically in comparison with the theoretical distribution for neutral loci. The contribution of this new tool is not only to permit the identification of loci under selection but also to establish hypotheses about ecological factors that could exert the selection pressure responsible. In the future, such an approach may accelerate the process of hunting for functional genes at the population level.

Adaptive evolution of a stress response protein

BACKGROUND

Some cancers are mediated by an interplay between tissue damage, pathogens and localised innate immune responses, but the mechanisms that underlie these linkages are only beginning to be unravelled.

METHODS AND PRINCIPAL FINDINGS

Here we identify a strong signature of adaptive evolution on the DNA sequence of the mammalian stress response gene SEP53, a member of the epidermal differentiation complex fused-gene family known for its role in suppressing cancers. The SEP53 gene appears to have been subject to adaptive evolution of a type that is commonly (though not exclusively) associated with coevolutionary arms races. A similar pattern of molecular evolution was not evident in the p53 cancer-suppressing gene.

CONCLUSIONS

Our data thus raises the possibility that SEP53 is a component of the mucosal/epithelial innate immune response engaged in an ongoing interaction with a pathogen. Although the pathogenic stress mediating adaptive evolution of SEP53 is not known, there are a number of well-known candidates, in particular viruses with established links to carcinoma.

Adaptive evolution after duplication of penaeidin antimicrobial peptides

Penaeidin antimicrobial peptides in penaeid shrimps are an important component of their innate immune system that provides immunity against infection caused by several gram-positive bacteria and filamentous fungal species. Despite the knowledge on the identification and characterization of these peptides in penaeid shrimps, little is known about the evolutionary pattern of these peptides and the underlying genetic mechanisms that maintain high sequence diversities in the penaeidin gene family. Based on the phylogenetic analyses and maximum likelihood-based codon substitution analyses, here we present the convincing evidence that multiple copies of penaeidins have evolved by gene duplication, and positive Darwinian selection (adaptive evolution) is the likely cause of accelerated rate of amino acid substitutions among these duplicated genes. While the average ratio of non-synonymous to synonymous substitutions (omega) for the entire coding region of both active domains is 0.9805, few codon sites showed significantly higher omega (3.73). The likelihood ratio tests that compare models incorporating positive selection (omega>1) at certain codon sites with models not incorporating positive selection (omega<1), failed to reject (p=0) the evidence of positive Darwinian selection. The rapid adaptive evolution of this gene family might be directed by the pathogens and the faster rate of amino acid substitutions in the N-terminal proline-rich and C-terminal cysteine-rich domains could be due to their direct involvement in the protection against pathogens. When the host expose to different habitats/environment an accelerated rate of amino acid substitutions in both the active domains may also be expected.

Population structure and its effects on patterns of nucleotide polymorphism in teosinte (Zea mays ssp. parviglumis)

Surveys of nucleotide diversity in the wild ancestor of maize, Zea mays ssp. parviglumis, have revealed genomewide departures from the standard neutral equilibrium (NE) model. Here we investigate the degree to which population structure may account for the excess of rare polymorphisms frequently observed in species-wide samples. On the basis of sequence data from five nuclear and two chloroplast loci, we found significant population genetic structure among seven subpopulations from two geographic regions. Comparisons of estimates of population genetic parameters from species-wide samples and subpopulation-specific samples showed that population genetic subdivision influenced observed patterns of nucleotide polymorphism. In particular, Tajima's D was significantly higher (closer to zero) in subpopulation-specific samples relative to species-wide samples, and therefore more closely corresponded to NE expectations. In spite of these overall patterns, the extent to which levels and patterns of polymorphism within subpopulations differed from species-wide samples and NE expectations depended strongly on the geographic region (Jalisco vs. Balsas) from which subpopulations were sampled. This may be due to the demographic history of subpopulations in those regions. Overall, these results suggest that explicitly accounting for population structure may be important for studies examining the genetic basis of ecologically and agronomically important traits as well as for identifying loci that have been the targets of selection.

Molecular evidence for a species complex in the patagonian lizard Liolaemus bibronii and phylogeography of the closely related Liolaemus gracilis (Squamata: Liolaemini)

The lizard genus Liolaemus is endemic to temperate South America and includes 190 species. Liolaemus bibronii has a large geographic distribution and inhabits a great diversity of habitats, including the Monte, Steppe, and high Andean grassland environments. Liolaemus gracilis has a similar body size and shape to L. bibronii; the two are parapatrically distributed, and L. gracilis is also widely distributed. Here we use the mtDNA cytb sequence data of these two species to investigate lizard phylogeographic patterns in southern South America. L. bibronii is paraphyletic with respect to L. gracilis, Liolaemus ramirezae, Liolaemus robertmertensi and Liolaemus saxatilis; it is composed of many genetically different allopatric haploclades, some of which are reciprocally monophyletic. We also found evidence for introgression between L. bibronii and L. gracilis in the same area that introgression was hypothesized in the Liolaemus darwinii complex. We discuss the distribution of the major haploclades with inferences of their population histories, the concordance of these clades' distributions and histories with other lizard complexes studied with the same markers and methods, and taxonomic implications of these results.

Pattern of nucleotide substitution and divergence of prophenoloxidase in decapods

Despite the unprecedented development in identification and characterization of prophenoloxidase (proPO) in commercially important decapods, little is known about the evolutionary relationship, rate of amino acid replacement and differential selection pressures operating on proPO of different species of decapods. Here we report the evolutionary relationship among these nine decapod species based on proPO gene and types of selective pressures operating on proPO codon sites. Our analyses revealed that all the nine decapod species shared a common ancestor. The mean percentage sequence divergence at proPO gene was 34.4+/-0.6%. Pairwise estimates of nonsynonymous to synonymous ratio (omega) for Homarus americanus-H. gammarus is greater than one, therefore indicating adaptive evolution (functional diversification) of proPO in these two species. In contrast, strong purifying selection (omega<1) was observed in all other species pairs. However, phylogenetically closely related decapods revealed relatively higher omega value (omega=0.15+/-0.3) than the distantly related species pairs (omega=0.0075+/-0.005). These discrepancies could be due to higher fixation probability of beneficial mutation in closely related species. Maximum likelihood-based codon substitution analyses revealed a strong purifying selection operating on most of the codon sites, therefore suggesting proPO is functionally constrained (purifying selection). Codon substitution analyses have also revealed the evidence of strong purifying selection in haemocyanin subunits of decapods.

Combining population genomics and quantitative genetics: finding the genes underlying ecologically important traits

A central challenge in evolutionary biology is to identify genes underlying ecologically important traits and describe the fitness consequences of naturally occurring variation at these loci. To address this goal, several novel approaches have been developed, including 'population genomics,' where a large number of molecular markers are scored in individuals from different environments with the goal of identifying markers showing unusual patterns of variation, potentially due to selection at linked sites. Such approaches are appealing because of (1) the increasing ease of generating large numbers of genetic markers, (2) the ability to scan the genome without measuring phenotypes and (3) the simplicity of sampling individuals without knowledge of their breeding history. Although such approaches are inherently applicable to non-model systems, to date these studies have been limited in their ability to uncover functionally relevant genes. By contrast, quantitative genetics has a rich history, and more recently, quantitative trait locus (QTL) mapping has had some success in identifying genes underlying ecologically relevant variation even in novel systems. QTL mapping, however, requires (1) genetic markers that specifically differentiate parental forms, (2) a focus on a particular measurable phenotype and (3) controlled breeding and maintenance of large numbers of progeny. Here we present current advances and suggest future directions that take advantage of population genomics and quantitative genetic approaches - in both model and non-model systems. Specifically, we discuss advantages and limitations of each method and argue that a combination of the two provides a powerful approach to uncovering the molecular mechanisms responsible for adaptation.

Allelic diversity of metallothionein in Orchesella cincta (L.): traces of natural selection by environmental pollution

The advances made in statistical methods to detect selection from DNA sequence variation has resulted in an enormous increase in the number of studies reporting positive selection. However, a disadvantage of such statistical tests is that often no insight into the actual source of selection is obtained. Finer understanding of evolution can be obtained when those statistical tests are combined with field observations on allele frequencies. We assessed whether the metallothionein (mt) gene of Orchesella cincta (Collembola), which codes for a metal-binding protein, is subject to selection, by investigating alleles and allele frequencies among European metal-stressed and reference populations. Eight highly divergent alleles were resolved in Northwest Europe. At the nucleotide level, a total of 51 polymorphic sites (five of them implying amino-acid changes) were observed. Although statistical tests applied to the sequences alone showed no indication of selection, a G-test rejected the null hypothesis that alleles are homogeneously distributed over metal-stressed and reference populations. Analysis of molecular variance assigned a small, but significant amount of the total variance to differences between metal-stressed and non-stressed populations. In addition, it was shown that metal-stressed populations tend to be more genetically diversified at this locus than non-stressed ones. These results suggest that the mt gene and its surrounding DNA region are affected by environmental metal contamination. This study illustrates that, in addition to statistical tests, field observations on allele frequencies are needed to gain understanding of selection and adaptive evolution.

A multi-locus assessment of connectivity and historical demography in the bluehead wrasse (Thalassoma bifasciatum)

The pelagic larval stage of most coral reef fishes might allow extensive dispersal or, alternatively, some level of local recruitment might be important. Molecular markers can be used to obtain indirect estimates of dispersal to evaluate these alternatives, yet the extent of gene flow among populations is known for only a small number of species. The use of such markers must take into account the properties of the markers and the demographic history of the population when making inferences about current gene flow. In the Caribbean bluehead wrasse, Thalassoma bifasciatum, previous studies have found both substantial levels of local recruitment, in studies interpreting otolith microchemistry and, conversely, a lack of genetic differentiation inferred from mitochondrial DNA (mtDNA) restriction-fragment length polymorphism (RFLP) data and allozymes. However, if subtle differentiation exists, larger sample sizes and highly variable markers may be required to discern it. Here we present results from mitochondrial control region sequence and microsatellite data that indicate a lack of genetic differentiation at both small and large spatial scales. However, historical processes, such as changes in population size, may have affected the current distribution of genetic variation.

Selection on protein-coding genes of natural cyanobacterial populations

We examined the distribution of synonymous and non-synonymous changes in 12 protein-coding genes of natural populations of cyanobacteria to infer changes in gene functionality. By comparing mutation distributions within and across species using the McDonald-Kreitman test, we found data sets to contain evidence for purifying selection (hetR of Trichodesmium, nifH of Cylindrospermopsis raceborskii and rpoC1 of Anabaena lemmermannii) and positive selection (kaiC of Microcoleus chthonoplastes and rbcX of Anabaena and Aphanizomenon sp.). Other genes from the same set of clonal isolates (petB and rbcL in M. chthonoplastes and Anabaena/Aphanizomenon, respectively) did not harbour evidence for either form of selection. The results of branch models of codon evolution agreed fully with the results of the McDonald-Kreitman test in terms of significance and absolute value of the dN/dS estimates. The high frequency of gene-specific mutation patterns and their association with branches that separate closely related cyanobacterial genera suggest that evolutionary tests are suited to uncover gene-specific selective differentiation in cyanobacterial genomes. At the same time, given the lack of information about the history of cyanobacteria, analysis of larger numbers of protein-coding genes of clonal cyanobacterial isolates will produce more detailed pictures of the effects of natural selection.

An integrated approach for the comparative analysis of a multigene family: the nicotianamine synthase genes of barley

Recent genomic projects reveal that about half of the gene repertoire in plant genomes is made up by multigene families. In this paper, a set of structural and phylogenetic analyses have been applied to compare the differently sized nicotianamine synthase (NAS) gene families in barley and rice. Nicotianamine acts as a chelator of iron and other heavy metals and plays a key role in uptake, phloem transport and cytoplasmic distribution of iron, challenging efforts for the breeding of iron-efficient crop plants. Nine barley NAS genes have been mapped, and co-linearity of flanking genes in barley and rice was determined. The combined analyses reveal that the NAS multigene family members in barley originated through at least one duplication event that occurred before the divergence of rice and barley. Additional duplications appear to have occurred within each of the species. Although we detected no evidence for positive selection of recently duplicated genes within species, codon-based tests revealed evidence for positive selection having contributed to the divergence of some amino acids. The integrated comparative and phylogenetic analysis improved our current view of NAS gene family evolution, might facilitate the functional characterization of individual members and is applicable to other multigene families.

A genomic scan for divergent selection in a secondary contact zone between Atlantic and Gulf of Mexico oysters,Crassostrea virginica

The degree of population structure within species often varies considerably among loci. This makes it difficult to determine whether observed variance reflects neutral-drift stochasticity or locus-specific selection at one or more loci. This uncertainty is exacerbated when evolutionary equilibrium cannot be assumed and/or admixture potentially inflates genomic variance. Thus, the value of a 'genome scan', where locus-specific summary statistics are compared with a simulated neutral distribution among loci, may be limited in secondary contact zones if the null distribution is sensitive to the timing of secondary contact. Of particular interest here is the wide variance previously observed in locus-specific patterns of population structure between Atlantic and Gulf of Mexico populations of eastern oyster, Crassostrea virginica. To test the robustness of an equilibrium null model, we compared among-locus distributions of F(ST) simulated under migration-drift equilibrium and several nonequilibrium secondary contact histories. We then tested for evidence of divergent selection between two oyster populations on either side of a secondary contact zone using 215 amplified fragment length polymorphism (AFLP) loci. Constant-migration equilibrium and nonequilibrium secondary-contact simulations produced equivalent distributions of F(ST) when anchored by the global mean F(ST) observed in oysters, 0.0917. The 99th quantile of simulated neutral F(ST) encompassed most of the variation among oyster loci. Three AFLP loci exhibited F(ST) values higher than this threshold. Although no locus was significant after correcting for multiple tests, our results show in geographically clinal organisms: AFLPs can efficiently characterize the genomic distribution of F(ST); equilibrium models can be used to evaluate outliers; these procedures help focus research on genomic regions of interest.

The molecular population genetics of shoot development in Arabidopsis thaliana

Studies in Arabidopsis thaliana have provided us with a wealth of information about the genetic pathways that regulate plant morphogenesis. This developmental genetic treasure trove represents a fantastic resource for researchers interested in the microevolution of development. Several laboratories have begun using molecular population genetic analyses to investigate the evolutionary forces that act upon loci that regulate shoot morphogenesis. Much of this work has focused on coding sequence variation in transcription factors; however, recent studies have explored sequence variation in other types of proteins and in promoter regions. Several genes that regulate shoot development contain signatures of selective sweeps associated with positive selection or harbor putative balanced polymorphisms in coding and noncoding sequences. Other regulatory genes appear to be evolving neutrally, but have accumulated potentially deleterious replacement polymorphisms.

Comparative molecular population genetics of phycoerythrin locus in Prochlorococcus

As the only remainder type of phycobiliproteins in Prochlorococcus, the actual role of phycoerythrin still remains unknown. Previous studies revealed that two different forms of phycoerythrin gene were found in two ecotypes of Prochlorococcus that are specifically adapted to either high light (HL) or low light (LL) conditions. Here we analyze patterns of phycoerythrin nucleotide variation in the HL- and LL-Prochlorococcus populations. Our analyses reveal a significantly greater number of non-synonymous fixed substitutions in peB and peA than expected based on interspecific comparisons. This pattern of excess non-synonymous fixed substitutions is not seen in other five phycoerythrin-related genes (peZ/V/Y/T/S). Several neutrality statistical tests indicate an excess of rare frequency polymorphisms in the LL-Prochlorococcus data, but an excess of intermediate frequency polymorphisms in the HL-Prochlorococcus data. Distributions of the positively selected sites identified using the likelihood ratio test, when mapped onto the phycoerythrin tertiary structure, reveal that HL- and LL-phycoerythrin should be under different selective patterns. These findings may provide insights into the likely role of selection at the phycoerythrin locus and motivate further research to unveil the function of phycoerythrin in Prochlorococcus.

Molecular evolution of animal antimicrobial peptides: widespread moderate positive selection

An increasing number of studies in both vertebrates and invertebrates show that the evolution of antimicrobial peptides is driven by positive selection. Because these diverse molecules show potential for therapeutic applications, they are currently the targets of much structural and functional research, providing extensive background data for evolutionary studies. In this paper, patterns of molecular evolution in antimicrobial peptide genes are reviewed. Evidence for positive selection on antimicrobial peptides includes an excess of nonsynonymous nucleotide substitutions, an excess of charge-changing amino acid substitutions, nonneutral patterns of allelic variation, and functional assays in vivo and in vitro that show improved antimicrobial effects for derived sequence variants. Positive selection on antimicrobial peptides may be as common as, but perhaps weaker than, selection on the best-known example of adaptively evolving immunity genes, the major histocompatibility complex. Thus, antimicrobial peptides present a useful and underutilized model for the study of adaptive molecular evolution.

Microsatellites for ecologists: a practical guide to using and evaluating microsatellite markers

Recent improvements in genetic analysis and genotyping methods have resulted in a rapid expansion of the power of molecular markers to address ecological questions. Microsatellites have emerged as the most popular and versatile marker type for ecological applications. The rise of commercial services that can isolate microsatellites for new study species and genotype samples at reasonable prices presents ecologists with the unprecedented ability to employ genetic approaches without heavy investment in specialized equipment. Nevertheless, the lack of accessible, synthesized information on the practicalities and pitfalls of using genetic tools impedes ecologists' ability to make informed decisions on using molecular approaches and creates the risk that some will use microsatellites without understanding the steps needed to evaluate the quality of a genetic data set. The first goal of this synthesis is to provide an overview of the strengths and limitations of microsatellite markers and the risks, cost and time requirements of isolating and using microsatellites with the aid of commercial services. The second goal is to encourage the use and consistent reporting of thorough marker screening to ensure high quality data. To that end, we present a multistep screening process to evaluate candidate loci for inclusion in a genetic study that is broadly targeted to both novice and experienced geneticists alike.

Variation at two flowering time genes within and among populations of Arabidopsis thaliana: comparison with markers and traits

Flowering Locus C (FLC) and Frigida are two interacting genes controlling flowering time variation in Arabidopsis thaliana. Variation at these genes was surveyed in 12 A. thaliana populations sampled in France. These populations were also screened for variation at molecular markers [12 microsatellites and 19 cleaved amplified polymorphic sequence (CAPS) markers] and at seven quantitative traits measured with and without vernalization. Seven populations were highly polymorphic at markers (H(S) = 0.57 at microsatellites, 0.24 at CAPS) and showed heritable variation for bolting time and some other traits. Five populations were genetically fixed or nearly fixed. Q(ST) for bolting time without vernalization was significantly higher than F(ST), suggesting local divergent selection. One of the two haplotype groups at FLC (FLC(A)) was very predominant (frequency of 99%). The first exon of Frigida showed elevated nonsynonymous variation, and nine loss-of-function mutations were found throughout the gene. The association between loss-of-function and earlier bolting was confirmed. Overall, 18 Frigida haplotypes were detected. The pattern of variation at Frigida was largely similar to that found at markers and traits, with the same populations being fixed or highly diverse. Metapopulation dynamics is thus probably the main factor shaping genetic variation in A. thaliana. However, F(ST) for functional (FRI) vs. nonfunctional (FRI(Delta)) haplotypes was significantly higher than F(ST) at markers. This suggested that loss-of-function at Frigida is under local selection for flowering time.

Phylogeographic history and gene flow among giant Galápagos tortoises on southern Isabela Island

Volcanic islands represent excellent models with which to study the effect of vicariance on colonization and dispersal, particularly when the evolution of genetic diversity mirrors the sequence of geological events that led to island formation. Phylogeographic inference, however, can be particularly challenging for recent dispersal events within islands, where the antagonistic effects of land bridge formation and vicariance can affect movements of organisms with limited dispersal ability. We investigated levels of genetic divergence and recovered signatures of dispersal events for 631 Galápagos giant tortoises across the volcanoes of Sierra Negra and Cerro Azul on the island of Isabela. These volcanoes are among the most recent formations in the Galápagos (<0.7 million years), and previous studies based on genetic and morphological data could not recover a consistent pattern of lineage sorting. We integrated nested clade analysis of mitochondrial DNA control region sequences, to infer historical patterns of colonization, and a novel Bayesian multilocus genotyping method for recovering evidence of recent migration across volcanoes using eleven microsatellite loci. These genetic studies illuminate taxonomic distinctions as well as provide guidance to possible repatriation programs aimed at countering the rapid population declines of these spectacular animals.

Forest-tree population genomics and adaptive evolution

Forest trees have gained much attention in recent years as nonclassical model eukaryotes for population, evolutionary and ecological genomic studies. Because of low domestication, large open-pollinated native populations, and high levels of both genetic and phenotypic variation, they are ideal organisms to unveil the molecular basis of population adaptive divergence in nature. Population genomics, in its broad-sense definition, is an emerging discipline that combines genome-wide sampling with traditional population genetic approaches to understanding evolution. Here we briefly review traditional methods of studying adaptive genetic variation in forest trees, and describe a new, integrated population genomics approach. First, alleles (haplotypes) at candidate genes for adaptive traits and their effects on phenotypes need to be characterized via sequencing and association mapping. At this stage, functional genomics can assist in understanding gene action and regulation by providing detailed transcriptional profiles. Second, frequencies of alleles in native populations for causative single-nucleotide polymorphisms are estimated to identify patterns of adaptive variation across heterogeneous environments. Population genomics, through deciphering allelic effects on phenotypes and identifying patterns of adaptive variation at the landscape level, will in the future constitute a useful tool, if cost-effective, to design conservation strategies for forest trees.

DNA sequence variation and selection of tag single-nucleotide polymorphisms at candidate genes for drought-stress response in Pinus taeda L

Genetic association studies are rapidly becoming the experimental approach of choice to dissect complex traits, including tolerance to drought stress, which is the most common cause of mortality and yield losses in forest trees. Optimization of association mapping requires knowledge of the patterns of nucleotide diversity and linkage disequilibrium and the selection of suitable polymorphisms for genotyping. Moreover, standard neutrality tests applied to DNA sequence variation data can be used to select candidate genes or amino acid sites that are putatively under selection for association mapping. In this article, we study the pattern of polymorphism of 18 candidate genes for drought-stress response in Pinus taeda L., an important tree crop. Data analyses based on a set of 21 putatively neutral nuclear microsatellites did not show population genetic structure or genomewide departures from neutrality. Candidate genes had moderate average nucleotide diversity at silent sites (pi(sil) = 0.00853), varying 100-fold among single genes. The level of within-gene LD was low, with an average pairwise r2 of 0.30, decaying rapidly from approximately 0.50 to approximately 0.20 at 800 bp. No apparent LD among genes was found. A selective sweep may have occurred at the early-response-to-drought-3 (erd3) gene, although population expansion can also explain our results and evidence for selection was not conclusive. One other gene, ccoaomt-1, a methylating enzyme involved in lignification, showed dimorphism (i.e., two highly divergent haplotype lineages at equal frequency), which is commonly associated with the long-term action of balancing selection. Finally, a set of haplotype-tagging SNPs (htSNPs) was selected. Using htSNPs, a reduction of genotyping effort of approximately 30-40%, while sampling most common allelic variants, can be gained in our ongoing association studies for drought tolerance in pine.

Molecular evolution of the sheep prion protein gene

Transmissible spongiform encephalopathies (TSEs) are infectious, fatal neurodegenerative diseases characterized by aggregates of modified forms of the prion protein (PrP) in the central nervous system. Well known examples include variant Creutzfeldt-Jakob Disease (vCJD) in humans, BSE in cattle, chronic wasting disease in deer and scrapie in sheep and goats. In humans, sheep and deer, disease susceptibility is determined by host genotype at the prion protein gene (PRNP). Here I examine the molecular evolution of PRNP in ruminants and show that variation in sheep appears to have been maintained by balancing selection, a profoundly different process from that seen in other ruminants. Scrapie eradication programs such as those recently implemented in the UK, USA and elsewhere are based on the assumption that PRNP is under positive selection in response to scrapie. If, as these data suggest, that assumption is wrong, eradication programs will disrupt this balancing selection, and may have a negative impact on the fitness or scrapie resistance of national flocks.

Molecular Signatures of Natural Selection

There is an increasing interest in detecting genes, or genomic regions, that have been targeted by natural selection. The interest stems from a basic desire to learn more about evolutionary processes in humans and other organisms, and from the realization that inferences regarding selection may provide important functional information. This review provides a nonmathematical description of the issues involved in detecting selection from DNA sequences and SNP data and is intended for readers who are not familiar with population genetic theory. Particular attention is placed on issues relating to the analysis of large-scale genomic data sets.

The evolution of immune-related genes from disease carrying mosquitoes: diversity in a peptidoglycan- and a thioester-recognizing protein

Adaptive polymorphism may be common in immune system genes as co-evolutionary interactions foster diversity; either through ongoing positive selection (arms races), or balancing selection. DNA sequence diversity in two putative immune system genes was examined in species of the genus Anopheles and from Aedes aegypti. For one gene, encoding the peptidoglycan recognizing protein PGRPLB, there was evidence of purifying selection, suggesting that selection acts to eliminate sequence variation. For another gene, encoding the thioester-containing protein TEP3, higher levels of amino acid replacement were found than would be expected under neutral models of evolution - an indication that this gene has been subject to repeated bouts of positive selection.

Challenges for identifying functionally important genetic variation: the promise of combining complementary research strategies

Strategies for the identification of functional genetic variation underlying phenotypic traits of ecological and evolutionary importance have received considerable attention in the literature recently. This paper aims to bring together and compare the relative strengths and limitations of various potentially useful research strategies for dissecting functionally important genetic variation in a wide range of organisms. We briefly explore the relative strengths and limitations of traditional and emerging approaches and evaluate their potential use in free-living populations. While it is likely that much of the progress in functional genetic analyses will rely on progress in traditional model species, it is clear that with prudent choices of methods and appropriate sampling designs, much headway can be also made in a diverse range of species. We suggest that combining research approaches targeting different functional and biological levels can potentially increase understanding the genetic basis of ecological and evolutionary processes both in model and non-model organisms.

The evolutionary ecology of the major histocompatibility complex

The major histocompatibility complex (MHC) has become a paradigm for how selection can act to maintain adaptively important genetic diversity in natural populations. Here, we review the contribution of studies on the MHC in non-model species to our understanding of how selection affects MHC diversity, emphasising how ecological and ethological processes influence the tempo and mode of evolution at the MHC, and conversely, how variability at the MHC affects individual fitness, population dynamics and viability. We focus on three main areas: the types of information that have been used to detect the action of selection on MHC genes; the relative contributions of parasite-mediated and sexual selection on the maintenance of MHC diversity; and possible future lines of research that may help resolve some of the unanswered issues associated with MHC evolution.

How to use molecular marker data to measure evolutionary parameters in wild populations

Estimating the genetic basis of phenotypic traits and the selection pressures acting on them are central to our understanding of the evolution and conservation of wild populations. However, obtaining such evolutionary-related parameters is not an easy task as it requires accurate information on both relatedness among individuals and their breeding success. Polymorphic molecular markers are very useful in estimating relatedness between individuals and parentage analyses are now extensively used in most taxa. The next step in the application of molecular data to wild populations is to use them to derive estimates of evolutionary-related parameters for quantitative traits, such as quantitative genetic parameters (e.g. heritability, genetic correlations) and measures of selection (e.g. selection gradients). Despite their great appeal and potential, the optimal use of molecular tools is still debated and it remains unclear how they should best be used to obtain reliable estimates of evolutionary parameters in the wild. Here, we review the methods available for estimating quantitative genetic and selection parameters and discuss their merits and shortcomings, to provide a tool that summarizes the potential uses of molecular data to obtain such parameters in wild populations.

Comparative analyses of genetic diversities within tomato and pepper collections detected by retrotransposon-based SSAP, AFLP and SSR

The retrotransposon-based sequence-specific amplification polymorphism (SSAP) marker system was used to assess the genetic diversities of collections of tomato and pepper industrial lines. The utility of SSAP markers was compared to that of amplified fragment length polymorphism (AFLP) and simple sequence repeat (SSR) markers. On the basis of our results, SSAP is most informative of the three systems for studying genetic diversity in tomato and pepper, with a significant correlation of genetic relationships between different SSAP datasets and between SSAP, AFLP and SSR markers. SSAP showed about four- to ninefold more diversity than AFLP and had the highest number of polymorphic bands per assay ratio and the highest marker index. For tomato, SSAP is more suitable for inferring overall genetic variation and relationships, while SSR has the ability to detect specific genetic relationships. All three marker results for pepper showed general agreement with pepper types. Additionally, retrotransposon sequences isolated from one species can be used in related Solanaceae genera. These results suggest that different marker systems are suited for studying genetic diversity in different contexts depending on the group studied, where discordance between different marker systems can be very informative for understanding genetic relationships within the study group.

Molecular evolution of daphnia immunity genes: polymorphism in a gram-negative binding protein gene and an alpha-2-macroglobulin gene

Studies of DNA polymorphism have shown that some immune system genes of mammals and plants are exceptionally diverse, indicating that coevolution between these taxa and their parasites mediates positive selective sweeps and/or balancing selection. The genes of the arthropod immune system remain comparatively unstudied. We isolated two putative immune system genes from the cladoceran crustacean Daphnia and examined DNA sequence diversity. For one gene, encoding a putative gram-negative binding protein, we found evidence of only purifying selection, indicating that this gene is under strong functional constraint and that selection acts to eliminate amino acid variation. For another gene, encoding a putative alpha-2-macroglobulin, we found evidence of positive selection, indicating the possible involvement of this gene in a host-parasite arms race. We discuss the assumed function of these genes and offer speculation regarding which components of the arthropod immune system might experience diversifying adaptive evolution.

Statistical tests of the coalescent model based on the haplotype frequency distribution and the number of segregating sites

Several tests of neutral evolution employ the observed number of segregating sites and properties of the haplotype frequency distribution as summary statistics and use simulations to obtain rejection probabilities. Here we develop a "haplotype configuration test" of neutrality (HCT) based on the full haplotype frequency distribution. To enable exact computation of rejection probabilities for small samples, we derive a recursion under the standard coalescent model for the joint distribution of the haplotype frequencies and the number of segregating sites. For larger samples, we consider simulation-based approaches. The utility of the HCT is demonstrated in simulations of alternative models and in application to data from Drosophila melanogaster.

Functional divergence in tandemly duplicated Arabidopsis thaliana trypsin inhibitor genes

In multigene families, variation among loci and alleles can contribute to trait evolution. We explored patterns of functional and genetic variation in six duplicated Arabidopsis thaliana trypsin inhibitor (ATTI) loci. We demonstrate significant variation in constitutive and herbivore-induced transcription among ATTI loci that show, on average, 65% sequence divergence. Significant variation in ATTI expression was also found between two molecularly defined haplotype classes. Population genetic analyses for 17 accessions of A. thaliana showed that six ATTI loci arranged in tandem within 10 kb varied 10-fold in nucleotide diversity, from 0.0009 to 0.0110, and identified a minimum of six recombination events throughout the tandem array. We observed a significant peak in nucleotide and indel polymorphism spanning ATTI loci in the interior of the array, due primarily to divergence between the two haplotype classes. Significant deviation from the neutral equilibrium model for individual genes was interpreted within the context of intergene linkage disequilibrium and correlated patterns of functional differentiation. In contrast to the outcrosser Arabidopsis lyrata for which recombination is observed even within ATTI loci, our data suggest that response to selection was slowed in the inbreeding, annual A. thaliana because of interference among functionally divergent ATTI loci.

Mitochondrial gene diversity in the common vole Microtus arvalis shaped by historical divergence and local adaptations

The phylogeography of the common vole (Microtus arvalis) was examined by analysing mitochondrial DNA (mtDNA) sequence variation in 1044 base pairs (bp) of the cytochrome b (cytb) gene and in 322 bp of the control region (ctr) among 106 individuals from 58 locations. The geographical distribution of four previously recognized cytb evolutionary lineages in Europe was refined and a new lineage was found in southern Germany. All lineages were distributed allopatrically, except in one sample that was probably located in a contact zone. The occurrence of several lineages in the Alps is in keeping with their recent recolonization from distinct sources. The translation of 84 cytb DNA sequences produced 33 distinct proteins with relationships that differed from those of the DNA haplotypes, suggesting that the mtDNA lineages did not diverge in response to selection. In comparison with M. agrestis, a neutrality test detected no overall evidence for selection in the cytb gene, but a closer examination of a structural model showed that evolutionarily conserved and functionally important positions were often affected. A new phylogeographical test of random accumulation of nonsynonymous mutations generated significant results in three lineages. We therefore conclude that the molecular diversity of cytb in M. arvalis is overall the result of the demographic history of the populations, but that there have been several episodes of local adaptation to peculiar environments.

Molecular population genetic analysis differentiates two virulence mechanisms of the fungal avirulence gene NIP1

Deletion or alteration of an avirulence gene are two mechanisms that allow pathogens to escape recognition mediated by the corresponding resistance gene in the host. We studied these two mechanisms for the NIP1 avirulence gene in field populations of the fungal barley pathogen Rhynchosporium secalis. The product of the avirulence gene, NIP1, causes leaf necrosis and elicits a defense response on plants with the Rrs1 resistance gene. A high NIP1 deletion frequency (45%) was found among 614 isolates from different geographic populations on four continents. NIP1 was also sequenced for 196 isolates, to identify DNA polymorphisms and corresponding NIP1 types. Positive diversifying selection was found to act on NIP1. A total of 14 NIP1 types were found, 11 of which had not been described previously. The virulence of the NIP1 types was tested on Rrs1 and rrs1 barley lines. Isolates carrying three of these types were virulent on the Rrs1 cultivar. One type each was found in California, Western Europe, and Jordan. Additionally, a field experiment with one pair of near-isogenic lines was conducted to study the selection pressure imposed by Rrs1 on field populations of R. secalis. Deletion of NIP1 was the only mechanism used to infect the Rrs1 cultivar in the field experiment. In this first comprehensive study on the population genetics of a fungal avirulence gene, virulence to Rrs1 in R. secalis was commonly achieved through deletion of the NIP1 avirulence gene but rarely also through point mutations in NIP1.

On the effects of background selection in small populations on comparisons of molecular variation

Deleterious mutations affect genetic variation at linked neutral loci. Neutral variation can be reduced due to background selection, but in small population and with tight linkage such variation may increase due to associative overdominance. Here I report the results of computer simulations of diploid genotypes in small populations, where I look at the effect of deleterious mutations and linkage on comparisons of intra- and interspecific variation. Each chromosome consisted of 2000 loci where deleterious and neutral mutations occurred. The ratio of nonsynonymous to synonymous substitution rates (Ka/Ks) either increases with tight linkage or is unaffected, depending on the strength of selection. The ratio of the numbers of segregating mutations to the number of fixed mutations decreases under the conditions leading to background selection but can increase at tight linkage. Numbers of segregating sites (Sn) are less affected than nucleotide site diversity (pi), pi reduces more than Sn at intermediate linkage, but pi increases more than Sn when linkage is tight. Similar effects as found for Sn and pi are observed for heterozygosity and variance in allele size of tandem repeat loci.

Evidence for Positive Selection at the Pantophysin (Pan I) Locus in Walleye Pollock, Theragra chalcogramma

Nucleotide polymorphism at the pantophysin (Pan I) locus in walleye pollock, Theragra chalcogramma, was examined using DNA sequence data. Two distinct allelic lineages were detected in pollock, resulting from three amino acid replacement mutations in the first intravesicular domain of the protein. The common Pan I allelic group, comprising 94% of the samples, was less polymorphic (pi = 0.005) than the uncommon group (pi = 0.008), and nucleotide diversity in both was higher than for two allelic lineages in the related Atlantic cod, Gadus morhua. Phylogenetic analyses of Pan I sequences from these two species did not clearly resolve orthology among allelic groups, in part because of recombination that has occurred between the two pollock lineages. Conventional tests of neutrality comparing polymorphisms within and between homologous regions of the Pan I locus in walleye pollock and Atlantic cod did not detect the effects of selection. This result is likely attributed to low levels of synonymous divergence among allelic lineages and a lack of mutation-drift equilibrium inferred from nucleotide mismatch frequency distributions. However, the ratio of nonsynonymous to synonymous substitutions per site (dN/dS) exceeded unity in two intravesicular domains of the protein and the influence of positive selection at multiple codon sites was strongly inferred through the use of maximum-likelihood analyses. In addition, the frequency spectrum of linked neutral variation showed indirect effects of adaptive hitchhiking in pollock resulting from a selective sweep of the common allelic lineage. Recombination between the two allelic classes may have prevented complete loss of the older, more polymorphic lineage. The results suggest that recurrent sweeps driven by positive selection is the principle mode of evolution at the Pan I locus in gadid fishes.