Estimating a nucleotide substitution rate for maize from polymorphism at a major domestication locus

To estimate a rate for single nucleotide substitutions for maize (Zea mays ssp. mays), we have taken advantage of data from genetic and archaeological studies of the domestication of maize from its wild ancestor, teosinte (Z. mays ssp. parviglumis). Genetic studies have shown that the teosinte branched1 (tb1) gene was a major target of human selection during maize domestication, and sequence diversity in the intergenic region 5' to the tb1-coding sequence is extraordinarily low. We show that polymorphism in this region is consistent with new mutation following fixation for a small number of tb1 haplotypes during domestication. Archeological studies suggest that maize was domesticated approximately 6,250-10,000 years ago and subsequently the size of the maize population is thought to have expanded rapidly. Using the observed number of mutations within the region of selection at tb1, the approximate age of maize domestication, and approximations for the maize genealogy, we have derived estimates for the nucleotide substitution rate for the tb1 intergenic region. Using two approaches, one of which is a coalescent approach, we obtain rate estimates of approximately 2.9 x 10(-8) and 3.3 x 10(-8) substitutions per site per year. We also show that the pattern of polymorphism in the tb1 intergenic region appears to have been strongly affected by the mutagenic effect of DNA methylation. Excluding target sites of symmetric DNA methylation (CG and CNG sites) from analysis, the mutation rate estimates are reduced by approximately 50%-60%, while the rates for CG and CNG sites are nearly an order of magnitude higher. We use rate estimates from the tb1 region to estimate the timing of expansion of transposable elements in the maize genome and suggest that this expansion occurred primarily within the last million years.     

Allele frequency distribution under recurrent selective sweeps

The allele frequency of a neutral variant in a population is pushed either upward or downward by directional selection on a linked beneficial mutation ("selective sweeps"). DNA sequences sampled after the fixation of the beneficial allele thus contain an excess of rare neutral alleles. This study investigates the allele frequency distribution under selective sweep models using analytic approximation and simulation. First, given a single selective sweep at a fixed time, I derive an expression for the sampling probabilities of neutral mutants. This solution can be used to estimate the time of the fixation of a beneficial allele from sequence data. Next, I obtain an approximation to mean allele frequencies under recurrent selective sweeps. Under recurrent sweeps, the frequency spectrum is skewed toward rare alleles. However, the excess of high-frequency derived alleles, previously shown to be a signature of single selective sweeps, disappears with recurrent sweeps. It is shown that, using this approximation and multilocus polymorphism data, genomewide parameters of directional selection can be estimated.     

The signature of positive selection on standing genetic variation

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

Overexpression of the maize Teosinte Branched1 gene in wheat suppresses tiller development

The number of viable shoots influences the overall architecture and productivity of wheat (Triticum aestivum L.). The development of lateral branches, or tillers, largely determines the resultant canopy. Tillers develop from the outgrowth of axillary buds, which form in leaf axils at the crown of the plant. Tiller number can be reduced if axillary buds are not formed or if the outgrowth of these buds is restricted. The teosinte branched1 (tb1) gene in maize, and homologs in rice and Arabidopsis, genetically regulate vegetative branching. In maize, increased expression of the tb1 gene restricts the outgrowth of axillary buds into lateral branches. In this study, the maize tb1 gene was introduced through transformation into the wheat cultivar "Bobwhite" to determine the effect of tb1 overexpression on wheat shoot architecture. Examination of multiple generations of plants reveals that tb1 overexpression in wheat results in reduced tiller and spike number. In addition, the number of spikelets on the spike and leaf number were significantly greater in tb1-expressing plants, and the height of these plants was also reduced. These data reveal that the function of the tb1 gene and genetic regulation of lateral branching via the tb1 mode of action is conserved between wheat, rice, maize and Arabidopsis. Thus, the tb1 gene can be used to alter plant architecture in agriculturally important crops like wheat.     

Historical population size change of bowhead whales inferred from DNA sequence polymorphism data

Nucleotide sequence data from the mitochondrial control region were used from a phylogenetic context to investigate the long-term history of a population of bowhead whales (Balaena mysticetus). In addition, the coalescence time of these sequences was used to estimate the age of the inferred patterns of population size change. The results indicate that mitochondrial genetic polymorphism was not affected by a recent bottleneck that occurred near the turn of the 20th century, thereby preserving the signature of historical population size change in the mitochondrial genome. Further analysis showed that this population underwent an expansion initiated in the Middle to Late Pleistocene. As such, early Holocene changes in Arctic sea ice distribution appear to have had little influence on patterns of genetic variability in this population.     

Expression patterns and mutant phenotype of teosinte branched1 correlate with growth suppression in maize and teosinte

The evolution of domesticated maize from its wild ancestor teosinte is a dramatic example of the effect of human selection on agricultural crops. Maize has one dominant axis of growth, whereas teosinte is highly branched. The axillary branches in maize are short and feminized whereas the axillary branches of teosinte are long and end in a male inflorescence under normal growth conditions. Previous QTL and molecular analysis suggested that the teosinte branched1 (tb1) gene of maize contributed to the architectural difference between maize and teosinte. tb1 mutants of maize resemble teosinte in their overall architecture. We analyzed the tb1 mutant phenotype in more detail and showed that the highly branched phenotype was due to the presence of secondary and tertiary axillary branching, as well as to an increase in the length of each node, rather than to an increase in the number of nodes. Double-mutant analysis with anther ear1 and tassel seed2 revealed that the sex of the axillary inflorescence was not correlated with its length. RNA in situ hybridization showed that tb1 was expressed in maize axillary meristems and in stamens of ear primordia, consistent with a function of suppressing growth of these tissues. Expression in teosinte inflorescence development suggests a role in pedicellate spikelet suppression. Our results provide support for a role for tb1 in growth suppression and reveal the specific tissues where suppression may occur.     

Tillering in the sugary1 sweet corn is maintained by overriding the teosinte branched1 repressive signal

The evolution of apical dominance in maize during domestication from teosinte is associated with higher expression from the teosinte branched1 (tb1) gene that inhibits tiller bud outgrowth. Unlike many standard maize varieties, the sweet corn inbred P39 that carries a mutation in a starch biosynthesis gene sugary1 produces multiple tillers and providing an opportunity to explore the diversification of the tb1 signal in maize. Through gene expression analysis, we show that tiller buds in P39 continue to grow by overriding the high expression level of tb1 that arrests bud outgrowth in maize inbred B73. In addition, we demonstrate that while B73 is largely non-responsive to shade, both P39 and teosinte respond through tb1-independent and tb1-dependent molecular mechanisms, respectively, leading to inhibition of tiller bud outgrowth.     

Identification of a locus under complex positive selection in Drosophila simulans by haplotype mapping and composite-likelihood estimation

The recent action of positive selection is expected to influence patterns of intraspecific DNA sequence variation in chromosomal regions linked to the selected locus. These effects include decreased polymorphism, increased linkage disequilibrium, and an increased frequency of derived variants. These effects are all expected to dissipate with distance from the selected locus due to recombination. Therefore, in regions of high recombination, it should be possible to localize a target of selection to a relatively small interval. Previously described patterns of intraspecific variation in three tandemly arranged, testes-expressed genes (janusA, janusB, and ocnus) in Drosophila simulans included all three of these features. Here we expand the original sample and also survey nucleotide polymorphism at three neighboring loci. On the basis of recombination events between derived and ancestral alleles, we localize the target of selection to a 1.5-kb region surrounding janusB. A composite-likelihood-ratio test based on the spatial distribution and frequency of derived polymorphic variants corroborates this result and provides an estimate of the strength of selection. However, the data are difficult to reconcile with the simplest model of positive selection, whereas a new composite-likelihood method suggests that the data are better described by a model in which the selected allele has not yet gone to fixation.     

The coalescent with selection on copy number variants

We develop a coalescent-based simulation tool to generate patterns of single nucleotide polymorphisms (SNPs) in a wide region encompassing both the original and duplicated genes. Selection on the new duplicated copy and interlocus gene conversion between the two copies are incorporated. This simulation enables us to explore how selection on duplicated copies affects the pattern of SNPs. The fixation of an advantageous duplicated copy causes a strong reduction in polymorphism not only in the duplicated copy but also in its flanking regions, which is a typical signature of a selective sweep by positive selection. After fixation, polymorphism gradually increases by accumulating neutral mutations and eventually reaches the equilibrium value if there is no gene conversion. When gene conversion is active, the number of SNPs in the duplicated copy quickly increases by transferring SNPs from the original copy; therefore, the time when we can recognize the signature of selection is decreased. Because this effect of gene conversion is restricted only to the duplicated region, more power to detect selection is expected if a flanking region to the duplicated copy is used.     

Evidence for a natural allelic series at the Maize domestication Locus teosinte branched1

Despite numerous quantitative trait loci and association mapping studies, our understanding of the extent to which natural allelic series contribute to the variation for complex traits is limited. In this study, we investigate the occurrence of a natural allelic series for complex traits at the teosinte branched1 (tb1) gene in natural populations of teosinte (Zea mays ssp. parviglumis, Z. mays ssp. mexicana, and Z. diploperennis). Previously, tb1 was shown to confer large effects on both plant architecture and ear morphology between domesticated maize and teosinte; however, the effect of tb1 on trait variation in natural populations of teosinte has not been investigated. We compare the effects of nine teosinte alleles of tb1 that were introgressed into an isogenic maize inbred background. Our results provide evidence for a natural allelic series at tb1 for several complex morphological traits. The teosinte introgressions separate into three distinct phenotypic classes, which correspond to the taxonomic origin of the alleles. The effects of the three allelic classes also correspond to known morphological differences between the teosinte taxa. Our results suggest that tb1 contributed to the morphological diversification of teosinte taxa as well as to the domestication of maize.     

No association between the angiotensin-converting enzyme ID polymorphism and elite endurance athlete status.

Several studies have reported that the insertion (I) allele of the angiotensin-converting enzyme (ACE) I/deletion (D) polymorphism is associated with enhanced responsiveness to endurance training and is more common in endurance athletes than in sedentary controls. We tested the latter hypothesis in a cohort of 192 male endurance athletes with maximal oxygen uptake >/=75 ml. kg(-1). min(-1) and 189 sedentary male controls. The ACE ID polymorphism in intron 16 was typed with the three-primer polymerase chain reaction method. Both the genotype (P = 0.214) and allele (P = 0.095) frequencies were similar in the athletes and the controls. Further analyses in the athletes revealed no excess of the I allele among the athletes within the highest quartile (> 80 ml. kg(-1). min(-1)) or decile (>83 ml. kg(-1). min(-1)) of maximal oxygen uptake. These data from the GENATHLETE cohort do not support the hypothesis that the ACE ID polymorphism is associated with a higher cardiorespiratory endurance performance level.     

Molecular differentiation at nuclear loci in French host races of the European corn borer (Ostrinia nubilalis)

French populations of the European corn borer consist of two sympatric and genetically differentiated host races. As such, they are well suited to study processes that could be involved in sympatric speciation, but the initial conditions of host-race divergence need to be elucidated. Gene genealogies can provide insight into the processes involved in speciation. We used DNA sequences of four nuclear genes to (1) document the genetic structure of the two French host races previously delineated with allozyme markers, (2) find genes directly or indirectly involved in reproductive isolation between host races, and (3) estimate the time since divergence of the two taxa and see whether this estimate is compatible with this divergence being the result of a host shift onto maize after its introduction into Europe approximately 500 years ago. Gene genealogies revealed extensive shared polymorphism, but confirmed the previously observed genetic differentiation between the two host races. Significant departures from the predictions of neutral molecular evolution models were detected at three loci but were apparently unrelated to reproductive isolation between host races. Estimates of time since divergence between French host races varied from approximately 75,000 to approximately 150,000 years, suggesting that the two taxa diverged recently but probably long before the introduction of maize into Europe.     

A framework for estimating the fixation time of an advantageous allele in stepping-stone models

Determining how population subdivision increases the fixation time of an advantageous allele is an important problem in evolutionary genetics as this influences many processes. Here, I lay out a framework for calculating the fixation time of a positively selected allele in a subdivided population, as a function of the number of demes present, the migration rate between them and the manner in which they are connected. Using this framework, it becomes clear that a beneficial allele's fixation time is significantly reduced through migration continuously introducing copies of the allele into a newly colonized subpopulation, increasing its frequency within these demes. The effect that migration has on allele frequency needs to be explicitly taken into account to produce a realistic estimate of fixation time. This behaviour is most prominent when demes are arranged on a two-dimensional torus, in comparison with populations where demes are arranged in a circle. This is because each subpopulation is connected to several neighbours over a torus, so that there are multiple paths that an allele can take in order to fix. As a consequence, some demes experience a greater influx and efflux of migrants than others. Analytical results are found to be very accurate when compared to stochastic simulations, and are generally robust if there are a large number of demes, or if the allele is weakly selected for.     

Patterns of molecular evolution associated with two selective sweeps in the Tb1-Dwarf8 region in maize

We focused on a region encompassing a major maize domestication locus, Tb1, and a locus involved in the flowering time variation, Dwarf8 (D8), to investigate the consequences of two closely linked selective sweeps on nucleotide variation and gain some insights into maize geographical diffusion, through climate adaptation. First, we physically mapped D8 at approximately 300 kb 3' of Tb1. Second, we analyzed patterns of nucleotide variation at Tb1, D8, and seven short regions (400-700 bp) located in the Tb1-D8 region sequenced on a 40 maize inbred lines panel encompassing early-flowering temperate and late-flowering tropical lines. The pattern of polymorphism along the region is characterized by two valleys of depleted polymorphism while the region in between exhibits an appreciable amount of diversity. Our results reveal that a region approximately 100 kb upstream of the D8 gene exhibits hallmarks of divergent selection between temperate and tropical lines and is likely closer than the D8 gene to the target of selection for climate adaptation. Selection in the tropical lines appears more recent than in the temperate lines, suggesting an initial domestication of early-flowering maize. Simulation results indicate that the polymorphism pattern is consistent with two interfering selective sweeps at Tb1 and D8.     

A Coalescent Model for a Sweep of a Unique Standing Variant

The use of genetic polymorphism data to understand the dynamics of adaptation and identify the loci that are involved has become a major pursuit of modern evolutionary genetics. In addition to the classical “hard sweep” hitchhiking model, recent research has drawn attention to the fact that the dynamics of adaptation can play out in a variety of different ways and that the specific signatures left behind in population genetic data may depend somewhat strongly on these dynamics. One particular model for which a large number of empirical examples are already known is that in which a single derived mutation arises and drifts to some low frequency before an environmental change causes the allele to become beneficial and sweeps to fixation. Here, we pursue an analytical investigation of this model, bolstered and extended via simulation study. We use coalescent theory to develop an analytical approximation for the effect of a sweep from standing variation on the genealogy at the locus of the selected allele and sites tightly linked to it. We show that the distribution of haplotypes that the selected allele is present on at the time of the environmental change can be approximated by considering recombinant haplotypes as alleles in the infinite-alleles model. We show that this approximation can be leveraged to make accurate predictions regarding patterns of genetic polymorphism following such a sweep. We then use simulations to highlight which sources of haplotypic information are likely to be most useful in distinguishing this model from neutrality, as well as from other sweep models, such as the classic hard sweep and multiple-mutation soft sweeps. We find that in general, adaptation from a unique standing variant will likely be difficult to detect on the basis of genetic polymorphism data from a single population time point alone, and when it can be detected, it will be difficult to distinguish from other varieties of selective sweeps. Samples from multiple populations and/or time points have the potential to ease this difficulty.     

Recent and recurrent selective sweeps of the antiviral RNAi gene Argonaute-2 in three species of Drosophila

Antagonistic host-parasite interactions can drive rapid adaptive evolution in genes of the immune system, and such arms races may be an important force shaping polymorphism in the genome. The RNA interference pathway gene Argonaute-2 (AGO2) is a key component of antiviral defense in Drosophila, and we have previously shown that genes in this pathway experience unusually high rates of adaptive substitution. Here we study patterns of genetic variation in a 100-kbp region around AGO2 in three different species of Drosophila. Our data suggest that recent independent selective sweeps in AGO2 have reduced genetic variation across a region of more than 50 kbp in Drosophila melanogaster, D. simulans, and D. yakuba, and we estimate that selection has fixed adaptive substitutions in this gene every 30-100 thousand years. The strongest signal of recent selection is evident in D. simulans, where we estimate that the most recent selective sweep involved an allele with a selective advantage of the order of 0.5-1% and occurred roughly 13-60 Kya. To evaluate the potential consequences of the recent substitutions on the structure and function of AGO2, we used fold-recognition and homology-based modeling to derive a structural model for the Drosophila protein, and this suggests that recent substitutions in D. simulans are overrepresented at the protein surface. In summary, our results show that selection by parasites can consistently target the same genes in multiple species, resulting in areas of the genome that have markedly reduced genetic diversity.     

Natural selection for the Duffy-null allele in the recently admixed people of Madagascar

While gene flow between distantly related populations is increasingly recognized as a potentially important source of adaptive genetic variation for humans, fully characterized examples are rare. In addition, the role that natural selection for resistance to vivax malaria may have played in the extreme distribution of the protective Duffy-null allele, which is nearly completely fixed in mainland sub-Saharan Africa and absent elsewhere, is controversial. We address both these issues by investigating the evolution of the Duffy-null allele in the Malagasy, a recently admixed population with major ancestry components from both East Asia and mainland sub-Saharan Africa. We used genome-wide genetic data and extensive computer simulations to show that the high frequency of the Duffy-null allele in Madagascar can only be explained in the absence of positive natural selection under extreme demographic scenarios involving high genetic drift. However, the observed genomic single nucleotide polymorphism diversity in the Malagasy is incompatible with such extreme demographic scenarios, indicating that positive selection for the Duffy-null allele best explains the high frequency of the allele in Madagascar. We estimate the selection coefficient to be 0.066. Because vivax malaria is endemic to Madagascar, this result supports the hypothesis that malaria resistance drove fixation of the Duffy-null allele in mainland sub-Saharan Africa.     

Mechanisms of reinforcement in natural and simulated polymorphic populations

Reinforcement speciation is the process whereby selection against hybrids drives the evolution of enhanced pre‐mating reproductive isolation. Work has focused on divergent mating preferences (assortative mating) but pre‐mating isolation can also arise via various migration modification behaviours, such as divergent habitat preferences. The relative importance of these two different mechanisms of reinforcement remains unclear. A recent theoretical model (Yukilevich–True model) found that relative fixation probabilities between these mechanisms can vary. Additionally, natural populations of Timema cristinae walking‐sticks exhibit variation (polymorphism) in both mechanisms, generating questions about the patterns expected for allele frequencies prior to fixation, during the early stages of the speciation process. In the present study, we report: (1) new analyses examining the correlation between fixation probabilities for assortative mating and migration modification in the Yukilevich–True model; (2) novel simulations examining allele frequencies in polymorphic populations; and (3) empirical patterns of reinforcement in T. cristinae in the context of theoretical predictions. Simulations of both types yielded congruent results, revealing that the outcome of reinforcement was dependent on the strength of selection. Under weak selection, reinforcement by either mechanism is unlikely. Under intermediate selection, the conditions favoring the rise and fixation of one mechanism favored the rise and fixation of the other. However, assortative mating evolved somewhat more readily than migration modification. Populations of T. cristinae, which experience such intermediate selection, supported these predictions. Under strong selection, the evolution of migration modification generally interfered with the evolution of assortative mating by decreasing migration between populations, thereby reducing selection for assortative mating. Congruence of the results for allele frequencies versus fixation probabilities suggests that similar patterns of reinforcement are expected during different stages of the speciation process. © 2008 The Linnean Society of London, Biological Journal of the Linnean Society, 2008, 95 , 305–319.     

European ACP1*C allele has recessive deleterious effects on early life viability

The acid phosphatase locus (ACP1) is a classical polymorphism that has been surveyed in hundreds of human populations worldwide. Among individuals of European ancestry, the ACP1*C allele occurs with an average frequency of approximately 0.05, whereas it is nearly absent in all other human populations. It has been hypothesized that this allele is maintained by overdominant selection among European populations. Here, we analyze ACP1 protein polymorphism data from more than 50,000 individuals previously surveyed in 67 populations across Europe as well as inheritance data from more than 6,000 European parent-offspring pairs to assess the signature of natural selection currently acting on this allele. Although we see a significant excess of ACP1*C heterozygotes relative to Hardy-Weinberg expectations, we find no evidence that natural selection favors ACP1*C heterozygotes. Instead, ACP1*C appears to have a strongly deleterious and recessive fitness effect. We observed only 48.9% of expected homozygous offspring from heterozygous parents and significantly fewer homozygotes than expected within populations. Because parent-offspring pairs indicate a significant deficiency of ACP1*C homozygotes, we infer that viability selection is acting on ACP1*C homozygotes very early in life, perhaps before birth. We estimate that approximately 1.2% of all couples of European ancestry are composed of individuals who both carry the APC1*C allele. As such, selection against ACP1*C homozygotes may represent a nonnegligible contribution to the overall number of spontaneous abortions among women of European ancestry and may cause substantial fertility reductions among some combinations of parental genotypes.