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

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

Lack of Polymorphism on the Drosophila Fourth Chromosome Resulting from Selection

Evolutionary processes can be inferred from comparisons of intraspecific polymorphism and interspecific divergence. We sequenced a 1.1-kb fragment of the cubitus interruptus Dominant (ciD) locus located on the nonrecombining fourth chromosome for ten natural lines of Drosophila melanogaster and nine of Drosophila simulans. We found no polymorphism within D. melanogaster and a single polymorphism within D. simulans; divergence between the species was about 5%. Comparison with the alcohol dehydrogenase gene and its two flanking regions in D. melanogaster, for which comparable data are available, revealed a statistically significant departure from neutrality in all three tests. This lack of polymorphism in the ciD locus may reflect recent positive selective sweeps on the fourth chromosome with extreme hitchhiking generated by the lack of recombination. By simulation, we estimate there to be a 50% chance that the selective sweeps occurred within the past 30,000 years in D. melanogaster and 75,000 in D. simulans.     

Neutral and Nonneutral Mitochondrial Genetic Variation in Deep-Sea Clams from the Family Vesicomyidae

Nucleotide sequences at two mitochondrial genes from 57 individuals representing eight species of deep-sea clams (Vesicomyidae) were examined for variation consistent with the neutral model of molecular evolution. One gene, cytochrome oxidase subunit I (COI), deviated from the expectations of neutrality by containing an excess of intraspecific nonsynonymous polymorphism. Additionally, one species, Calyptogena kilmeri, showed a significant excess of rare polymorphism specifically at the COI locus. In contrast, a second mitochondrial gene, the large-subunit 16S ribosomal RNA gene (16S), showed little deviation from neutrality either between or within species. Together, COI and 16S show no deviation from neutral expectations by the HKA test, produce congruent phylogenetic relationships between species, and show correlated numbers of fixed differences between species and polymorphism within species. These patterns of both neutral and nonneutral evolution within the mitochondrial genome are most consistent with a model where intraspecific nonsynonymous polymorphism at COI is near neutrality. In addition to examining the forces of molecular evolution, we extend hypotheses about interspecific relationships within this family for geographical locations previously unexamined by molecular methods including habitats near the Middle Atlantic, the Aleutian Trench, and Costa Rica. Received: 10 March 1999 / Accepted: 13 September 1999     

The role of gene conversion in determining sequence variation and divergence in the Est-5 gene family in Drosophila pseudoobscura.

Nucleotide sequences of eight Est-5A and Est-5C genes corresponding to previously sequenced Est-5B genes in Drosophila pseudoobscura were determined to compare patterns of polymorphism and divergence among members of this small gene family. The three esterase genes were also sequenced from D. persimilis and D. miranda for interspecific comparisons. The data provide evidence that gene conversion between loci contributes to polymorphism and to the homogenization of the Est5 genes. For Est-5B, which encodes one of the most highly polymorphic proteins in Drosophila, 12% of the segregating amino acid variants appear to have been introduced via gene conversion from other members of the gene family. Interlocus gene conversion can also explain high sequence similarity, especially at synonymous sites, between Est-5B and Est-5A. Tests of neutrality using interspecific comparisons show that levels of polymorphism conform to neutral expectations at each Est-5 locus. However, McDonald-Kreitman tests based on intraspecific gene comparisons indicate that positive selection on amino acids has accompanied Est-5 gene duplication and divergence in D. pseudoobscura.     

DNA水平自然选择作用的检测

上个世纪60年代,Kimura提出的“中性进化”假说使经典的达尔文自然选择学说遭遇了前所未有的挑战。但新近的研究表明：在DNA水平,越来越多的证据支持“自然选择”的进化理论。这些研究成果得益于近年来大量群体和基因组DNA数据的积累,以及理论群体遗传学的发展。在DNA水平检测选择作用是否存在的方法包括两大类：种内多态性检验和种间差异度检验。前者以Tajima(1989)提出的D检验为代表,后者大都基于“中性条件下,种内与种间进化速率一致”的原理。这些方法以中性假说作为零假设,结合统计检验方法分析DNA数据,被称为“中性检验”。这些方法对于解决一些有关进化的基础理论问题和人类遗传学及生物信息学的深入研究都具有重要意义。本文介绍几个应用广泛的检测方法,以使国内的读者了解它们的基本思路和操作方法。     

A Nonparametric Test Reveals Selection for Rapid Flowering in the Arabidopsis Genome

The detection of footprints of natural selection in genetic polymorphism data is fundamental to understanding the genetic basis of adaptation, and has important implications for human health. The standard approach has been to reject neutrality in favor of selection if the pattern of variation at a candidate locus was significantly different from the predictions of the standard neutral model. The problem is that the standard neutral model assumes more than just neutrality, and it is almost always possible to explain the data using an alternative neutral model with more complex demography. Today's wealth of genomic polymorphism data, however, makes it possible to dispense with models altogether by simply comparing the pattern observed at a candidate locus to the genomic pattern, and rejecting neutrality if the pattern is extreme. Here, we utilize this approach on a truly genomic scale, comparing a candidate locus to thousands of alleles throughout the Arabidopsis thaliana genome. We demonstrate that selection has acted to increase the frequency of early-flowering alleles at the vernalization requirement locus FRIGIDA. Selection seems to have occurred during the last several thousand years, possibly in response to the spread of agriculture. We introduce a novel test statistic based on haplotype sharing that embraces the problem of population structure, and so should be widely applicable.     

A Nonparametric Test Reveals Selection for Rapid Flowering in the Arabidopsis Genome

The detection of footprints of natural selection in genetic polymorphism data is fundamental to understanding the genetic basis of adaptation, and has important implications for human health. The standard approach has been to reject neutrality in favor of selection if the pattern of variation at a candidate locus was significantly different from the predictions of the standard neutral model. The problem is that the standard neutral model assumes more than just neutrality, and it is almost always possible to explain the data using an alternative neutral model with more complex demography. Today's wealth of genomic polymorphism data, however, makes it possible to dispense with models altogether by simply comparing the pattern observed at a candidate locus to the genomic pattern, and rejecting neutrality if the pattern is extreme. Here, we utilize this approach on a truly genomic scale, comparing a candidate locus to thousands of alleles throughout the Arabidopsis thaliana genome. We demonstrate that selection has acted to increase the frequency of early-flowering alleles at the vernalization requirement locus FRIGIDA. Selection seems to have occurred during the last several thousand years, possibly in response to the spread of agriculture. We introduce a novel test statistic based on haplotype sharing that embraces the problem of population structure, and so should be widely applicable.     

HLA DNA sequence variation among human populations: molecular signatures of demographic and selective events

Molecular differences between HLA alleles vary up to 57 nucleotides within the peptide binding coding region of human Major Histocompatibility Complex (MHC) genes, but it is still unclear whether this variation results from a stochastic process or from selective constraints related to functional differences among HLA molecules. Although HLA alleles are generally treated as equidistant molecular units in population genetic studies, DNA sequence diversity among populations is also crucial to interpret the observed HLA polymorphism. In this study, we used a large dataset of 2,062 DNA sequences defined for the different HLA alleles to analyze nucleotide diversity of seven HLA genes in 23,500 individuals of about 200 populations spread worldwide. We first analyzed the HLA molecular structure and diversity of these populations in relation to geographic variation and we further investigated possible departures from selective neutrality through Tajima's tests and mismatch distributions. All results were compared to those obtained by classical approaches applied to HLA allele frequencies.Our study shows that the global patterns of HLA nucleotide diversity among populations are significantly correlated to geography, although in some specific cases the molecular information reveals unexpected genetic relationships. At all loci except HLA-DPB1, populations have accumulated a high proportion of very divergent alleles, suggesting an advantage of heterozygotes expressing molecularly distant HLA molecules (asymmetric overdominant selection model). However, both different intensities of selection and unequal levels of gene conversion may explain the heterogeneous mismatch distributions observed among the loci. Also, distinctive patterns of sequence divergence observed at the HLA-DPB1 locus suggest current neutrality but old selective pressures on this gene. We conclude that HLA DNA sequences advantageously complement HLA allele frequencies as a source of data used to explore the genetic history of human populations, and that their analysis allows a more thorough investigation of human MHC molecular evolution.     

Sequence variation and haplotype structure at the putative flowering-time locus COL1 of Brassica nigra

Motivated by a previous study indicating that polymorphism at an indel, Ind2, within the Brassica nigra COL1 gene is significantly associated with flowering time, we searched for evidence of selection in a sample of 41 complete sequences of B. nigra COL1. The within-gene population recombination rate is moderate, and all neutrality tests used in the present study failed to detect departure from the standard neutral model or evidence of selection. The haplotype structure of the 5'-half of the gene is primarily associated with the demographic history of the species and more specifically with the split between European and Ethiopian populations, whereas the structure of the 3'-half reflects the polymorphism at Ind2. This could be the result of selection or a combination of recombination and migration during the history of the sample of sequences. Without additional information on polymorphism in flanking areas, these two alternatives are difficult to tell apart. If selection acted on the gene, we suggest that if the indel itself is not the target of selection, among the polymorphic sites cosegregating with the polymorphism at Ind2, replacement polymorphisms around sites 890 and 1260 are the most likely quantitative trait nucleotides within the gene.     

Definition and Estimation of Higher-Order Gene Fixation Indices

Fixation indices summarize the associations between genes that arise from the joint effects of inbreeding and selection. In this paper, fixation indices are derived for pairs, triplets and quadruplets of genes at a single multiallelic locus. The fixation indices are obtained by dividing cumulants by constants; the cumulants describe the statistical distribution of alleles and the constants are functions of gene frequency. The use of cumulants instead of moments is necessary only for four-gene indices, when the fourth cumulant is used. A second type of four-gene index is also required, and this index is based upon the covariation of second-order cumulants. At multiallelic loci, a large number of indices is possible. If alleles are selectively neutral, the number of indices is reduced and the relationship between gene identity and gene cumulants is shown.--Two-gene indices can always be estimated from genotypic frequency data at a single polymorphic locus. Three-gene indices are also estimable except when allele frequency equals one-half. Four-gene indices are not estimable unless selection is assumed to have an equal effect upon each allele (such as under selective neutrality) and the locus contains at least three alleles of unequal frequency. For diallelic or selected loci, an alternative four-gene fixation index is proposed. This index incorporates both types of four-gene associations but cannot be related to gene identity.     

Insights into genome differentiation: pheromone-binding protein variation and population history in the European corn borer (Ostrinia nubilalis)

Examination of sequence variation at nuclear loci can give insights into population history and gene flow that cannot be derived from other commonly used molecular markers, such as allozymes. Here, we report on sequence variation at a single nuclear locus, the pheromone-binding protein (PBP) locus, in the European corn borer (Ostrinia nubilalis). The European corn borer has been divided into three races in New York State on the basis of differences in pheromone communication and life history. Previous allozyme data have suggested that there is a small but significant amount of genetic differentiation between these races. The PBP does not appear to be involved in the pheromone differences between these races. Examination of variation at the PBP locus in the three races reveals no fixed differences between races despite high levels of polymorphism. There also appears to have been considerable recombination in the history of the pheromone-binding protein alleles. Observation of both recombination between alleles and lack of significant nucleotide or insertion/deletion divergence between races leads us to suggest that these populations are either recently diverged or have continued to exchange genetic material subsequent to divergence in pheromone communication and life history.     

Genealogical footprints of speciation processes in wild tomatoes: demography and evidence for historical gene flow

Multilocus studies assessing patterns of nucleotide polymorphism within and among closely related species provide access to genealogical information bearing on demographic and geographic aspects of their speciation history. However, the technical difficulties in obtaining sufficient sequence data have severely limited this approach thus far, especially in outbred plant taxa. We employ the analytical framework of divergence population genetics in testing the isolation model of speciation in three self-incompatible species of wild tomatoes (clade Lycopersicon), in particular the assumption of divergence without gene flow. Based on DNA sequence data for 13 nuclear loci, average levels of silent polymorphism vary more than three-fold among species. We estimate a large effective population size for the ancestral species, quite similar to that of the highly polymorphic L. peruvianum. The other two species, however, exhibit concordant signatures of population-size reduction. These demographic inferences are biologically plausible and consistent with results obtained from standard neutrality tests. While the isolation model cannot be rejected by goodness-of-fit criteria, patterns of intragenic linkage disequilibrium in L. peruvianum are indicative of historical introgression at least in some regions of the genome. Considered jointly with the geographic pattern of postzygotic reproductive isolation, our results suggest that speciation occurred under residual gene flow, implying natural selection as one of the evolutionary forces driving the diversification of tomato lineages.     

Molecular Population Genetics of the Distal Portion of the X Chromosome in Drosophila: Evidence for Genetic Hitchhiking of the Yellow-Achaete Region

We have estimated DNA sequence variation and differentiation within and between Drosophila melanogaster and its sibling species, Drosophila simulans, using six-cutter restriction site variation at yellow-achaete (y-ac), phosphogluconate dehydrogenase (Pgd), and period (per). These three gene regions are of varying distance from the telomere of the X chromosome and range from very low to moderate rates of recombination in D. melanogaster. According to Tajima's test of neutrality, the Pgd region has been influenced by balancing selection in D. melanogaster. This is consistent with previous data suggesting the allozyme polymorphism at this locus is visible to selection. The Hudson, Kreitman, Aguadé test of neutrality reveals a significant departure from neutrality for the y-ac region compared to the per or rosy regions in D. simulans. There is also a significant departure for the y-ac region compared to the Adh 5' flanking region in D. melanogaster. In both species the departure appears to be due to reduced variation at y-ac compared to that expected from divergence between D. simulans and D. melanogaster. We conclude that recent hitchhiking associated with the selective fixation of one or more advantageous mutants in the y-ac region is the best explanation for reduced variation at y-ac.     

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 duplicated amylase gene regions in Drosophila melanogaster: evidence of positive selection in the coding regions and selective constraints in the cis-regulatory regions

In this study, we randomly sampled Drosophila melanogaster from Japanese and Kenyan natural populations. We sequenced duplicated (proximal and distal) Amy gene regions to test whether the patterns of polymorphism were consistent with neutral molecular evolution. F(st) between the two geographically distant populations, estimated from Amy gene regions, was 0.084, smaller than reported values for other loci, comparing African and Asian populations. Furthermore, little genetic differentiation was found at a microsatellite locus (DROYANETSB) in these samples (G'st = -0.018). The results of several tests (Tajima's, Fu and Li's, and Wall's tests) were not significantly different from neutrality. However, a significantly higher level of fixed replacement substitutions was detected by a modified McDonald and Kreitman test for both populations. This indicates that positive selection occurred during or immediately after the speciation of D. melanogaster. Sliding-window analysis showed that the proximal region 1, a part of the proximal 5' flanking region, was conserved between D. melanogaster and its sibling species, D. simulans. An HKA test was significant when the proximal region 1 was compared with the 5' flanking region of Alcohol dehydrogenase (Adh), indicating a severe selective constraint on the Amy proximal region 1. These results suggest that natural selection has played an important role in the molecular evolution of Amy gene regions in D. melanogaster.     

Testing the generalized neutrality hypothesis

Various tests of the hypothesis of selective neutrality based on gene frequency are now available. These tests take as null hypothesis the concept of “strict neutrality”: all new mutants are required to be selectively identical to each other. For evolutionary questions, however, (as opposed to those of genetic polymorphism), a wider null hypothesis might be of interest. Since deleterious alleles have essentially no evolutionary importance, one might wish to test the null hypothesis that only neutral or deleterious mutations occur. The principal alternative to this hypothesis is that there exists heterotic selection of some form for some alleles tending to maintain a level of genetic polymorphism higher than that under neutrality. In this paper an assessment is made of the usefulness of a test of strict neutrality first proposed by this author (Ewens, 1972) as a test of null hypothesis of “generalized neutrality,” i.e. that only neutral or deleterious alleles occur. At the same time some remarks will be made about estimation of the fundamental parameter θ defining these processes.     

Statistical power when testing for genetic differentiation

A variety of statistical procedures are commonly employed when testing for genetic differentiation. In a typical situation two or more samples of individuals have been genotyped at several gene loci by molecular or biochemical means, and in a first step a statistical test for allele frequency homogeneity is performed at each locus separately, using, e.g. the contingency chi-square test, Fisher's exact test, or some modification thereof. In a second step the results from the separate tests are combined for evaluation of the joint null hypothesis that there is no allele frequency difference at any locus, corresponding to the important case where the samples would be regarded as drawn from the same statistical and, hence, biological population. Presently, there are two conceptually different strategies in use for testing the joint null hypothesis of no difference at any locus. One approach is based on the summation of chi-square statistics over loci. Another method is employed by investigators applying the Bonferroni technique (adjusting the P-value required for rejection to account for the elevated alpha errors when performing multiple tests simultaneously) to test if the heterogeneity observed at any particular locus can be regarded significant when considered separately. Under this approach the joint null hypothesis is rejected if one or more of the component single locus tests is considered significant under the Bonferroni criterion. We used computer simulations to evaluate the statistical power and realized alpha errors of these strategies when evaluating the joint hypothesis after scoring multiple loci. We find that the 'extended' Bonferroni approach generally is associated with low statistical power and should not be applied in the current setting. Further, and contrary to what might be expected, we find that 'exact' tests typically behave poorly when combined in existing procedures for joint hypothesis testing. Thus, while exact tests are generally to be preferred over approximate ones when testing each particular locus, approximate tests such as the traditional chi-square seem preferable when addressing the joint hypothesis.     

EVIDENCE FOR SELECTION ON A CHORDATE HISTOCOMPATIBILITY LOCUS

Allorecognition is the ability of an organism to differentiate self or close relatives from unrelated individuals. The best known applications of allorecognition are the prevention of inbreeding in hermaphroditic species (e.g., the self‐incompatibility [SI] systems in plants), the vertebrate immune response to foreign antigens mediated by MHC loci, and somatic fusion, where two genetically independent individuals physically join to become a chimera. In the few model systems where the loci governing allorecognition outcomes have been identified, the corresponding proteins have exhibited exceptional polymorphism. But information about the evolution of this polymorphism outside MHC is limited. We address this subject in the ascidian Botryllus schlosseri, where allorecognition outcomes are determined by a single locus, called FuHC (Fusion/HistoCompatibility). Molecular variation in FuHC is distributed almost entirely within populations, with very little evidence for differentiation among different populations. Mutation plays a larger role than recombination in the creation of FuHC polymorphism. A selection statistic, neutrality tests, and distribution of variation within and among different populations all provide evidence for selection acting on FuHC, but are not in agreement as to whether the selection is balancing or directional.     

Nucleotide polymorphism in colicin E2 gene clusters: evidence for nonneutral evolution

To explore the molecular mechanisms behind the diversification of colicin gene clusters, we examined DNA sequence polymorphism for the colicin gene clusters of 14 colicin E2 (ColE2) plasmids obtained from natural isolates of Escherichia coli. Two types of ColE2 plasmids are revealed, with type II gene clusters generated by recombination between type I ColE2 and ColE7 gene clusters. The levels and patterns of DNA polymorphism are different between the two types. Type I polymorphism is distributed evenly along the gene cluster, while type II accumulates polymorphism at an elevated rate in the 5' end of the colicin gene. These differences may be explained by recombinational origins of type II gene clusters. The pattern of divergence between the ColE2 gene cluster and its close relative ColE9 is not correlated with the pattern of polymorphism within ColE2, suggesting that this gene cluster is not evolving in a neutral fashion. A statistical test confirms significant departures from the predictions of neutrality. These data lend further support to the hypothesis that colicin gene clusters may evolve under the influence of nonneutral forces.