Below are the first 10 and last 10 pages of uncorrected machine-read text (when available) of this chapter, followed by the top 30 algorithmically extracted key phrases from the chapter as a whole.
Intended to provide our own search engines and external engines with highly rich, chapter-representative searchable text on the opening pages of each chapter. Because it is UNCORRECTED material, please consider the following text as a useful but insufficient proxy for the authoritative book pages.
Do not use for reproduction, copying, pasting, or reading; exclusively for search engines.
OCR for page 205
In the Light of Evolution, Volume I: Adaptation and Complex Design 11 Plant Domestication, a Unique Opportunity to Identify the Genetic Basis of Adaptation JEFFREY ROSS-IBARRA, PETER L. MORRELL, and BRANDON S. GAUT Despite the fundamental role of plant domestication in human history and the critical importance of a relatively small number of crop plants to modern societies, we still know little about adaptation under domestication. Here we focus on efforts to identify the genes responsible for adaptation to domestication. We start from a historical perspective, arguing that Darwin’s conceptualization of domestication and unconscious selection provides valuable insight into the evolutionary history of crops and also provides a framework to evaluate modern methods used to decipher the genetic mechanisms underlying phenotypic change. We then review these methods, framing the discussion in terms of the phenotype–genotype hierarchy. Top-down approaches, such as quantitative trait locus and linkage disequilibrium mapping, start with a phenotype of interest and use genetic analysis to identify candidate genes. Bottom-up approaches, alternatively, use population genetic analyses to identify potentially adaptive genes and then rely on standard bioinformatics and reverse genetic tools to connect selected genes to a phenotype. We discuss the successes, advantages, and challenges of each, but we conclude that bottom-up approaches to understanding domestication as an adaptive process hold greater promise both for the study of Department of Ecology and Evolutionary Biology, University of California, Irvine, CA 92697-2525.
OCR for page 206
In the Light of Evolution, Volume I: Adaptation and Complex Design adaptation and as a means to identify genes that contribute to agronomically important traits. Plant domestication fundamentally altered the course of human history. The adaptation of plants to cultivation was vital to the shift from hunter–gatherer to agricultural societies, and it stimulated the rise of cities and modern civilization. Humans still rely on crops that were domesticated >10,000 years ago in such diverse places as Central America, New Guinea, and the Fertile Crescent. Nonetheless, modern humans are reliant on a surprisingly small number of crops: Nearly 70% of the calories consumed by humans are supplied by only 15 crops (Table 11.1). The cereals are particularly important, with five crops (rice, wheat, maize, sugarcane, and barley) contributing more than half of the calories consumed. Despite the critical importance of these crops, in most cases little is known about their domestication. Some obvious questions pertain to the domesticators: Who were they? How did they identify the incipient crop? What were their cultivation methods? Other questions concern crop history: What was the wild progenitor of the modern crop? Did domesti- TABLE 11.1 Major World Crops Ranked by Metric Tonnage Rank by Tonnage* Common Name Rank by Calories Consumed* Ploidy Propagation Life History 1 Sugarcane 4 8× O,V P 2 Maize 3 2× O A 3 Wheat 2 6× O A 4 Rice 1 2× S A 5 Potatoes 6 4× O,V AP 6 Sugar beet 8 2× O A 7 Soybeans 5 2× S A 8 Cassava 9 4× O,V AP 9 Palm kernel 7 2× O P 10 Barley 11 2× S A 11 Sweet potatoes 15 4–6× O,V AP 12 Tomatoes 30 2× S A 13 Watermelons 38 2× O A 14 Bananas 19 3× V P 15 Brassicas 37 2× O A O, outcrossing; S, selfing; V, vegetative; P, perennial; A, annual; AP, perennial species generally cultivated as annuals. *Data are from the Food and Agriculture Organization of the United Nations (www.fao.org, 2004).
OCR for page 207
In the Light of Evolution, Volume I: Adaptation and Complex Design cation occur more than once? If so, where? The application of phylogeographic methods is beginning to inform the answers to this latter set of questions (Heun et al., 1997; Harter et al., 2004), but the picture for any one crop remains far from complete. In this chapter we focus on a third set of questions that revolve around the phenotypic changes associated with domestication. The first question is whether phenotypic changes driven by artificial selection are an apt analogy for adaptation in nature. We take a historical perspective on this issue, arguing that Darwin’s conceptualization of domestication provides valuable insight into our view of crop history and provides a framework for evaluating methods used to decipher the genetic mechanisms underlying phenotypic change. We then review these methods, framing the discussion in terms of the phenotype–genotype hierarchy, evaluating both “top-down” and “bottom-up” approaches. We conclude by arguing that an appreciation of domestication as an adaptive process has the potential to reveal far more about the genes contributing to agronomic traits than has been learned to date. IS ARTIFICIAL SELECTION ANALOGOUS TO NATURAL SELECTION? A HISTORICAL PERSPECTIVE In the opening chapter of the Origin of Species, Charles Darwin introduced the idea of natural selection with an analogy to domestication (Darwin, 1859b). The importance of domestication to Darwin’s thinking is evident even in early sketches of his work (Darwin, 1909), and Darwin himself claimed that the example of domestication was fundamental to the formulation of his theory (C. Darwin, 1859b, 1868, 1958; for a different perspective see Rheinberger and McLaughlin, 1984). But for Darwin domestication was more than a useful analogy: he saw it as a model of adaptation from which inferences about the nature of variation and selection in natural systems could be drawn (Cornell, 1984; Rheinberger and McLaughlin, 1984). Darwin’s assertion of the importance of domestication in understanding the evolutionary process was not universally accepted, however. One of the most vocal critics of Darwin’s views on domestication was Alfred Russell Wallace. Even in their joint publication announcing the theory of natural selection, Wallace denies the relevance of domestication: “We see, then, that no inferences as to varieties in a state of nature can be deduced from the observation of those occurring among domestic animals” (Darwin and Wallace, 1858). “It has always been considered a weakness in Darwin’s work,” he later writes, “that he based his theory, primarily, on the evidence of variation in domesticated animals and cultivated plants” (Wallace, 1889).
OCR for page 208
In the Light of Evolution, Volume I: Adaptation and Complex Design Wallace found fault with two aspects of domestication as a heuristic for understanding adaptation in nature. He argued first that the analogy was flawed: artificial selection requires an intelligent selector, whereas no such force acts in natural systems. Additionally, he insisted that the selection itself was fundamentally different, leading to intrinsically different kinds of variation. Domesticated species, he wrote, “are abnormal, irregular, artificial; they are subject to varieties which never occur and never can occur in a state of nature: their very existence depends altogether on human care; so far are many of them removed from that just proportion of faculties, that true balance of organization, by means of which alone an animal left to its own resources can preserve its existence and continue its race” (Darwin and Wallace, 1858). Both of Wallace’s lines of argument find modern audiences, from those who see a fundamental difference between the conscious selection of humans and natural processes (Stebbins, 1950) to those who argue that variation in domesticated species differs from that in nature (Coyne and Lande, 1985). Yet Darwin directly addressed these ideas with his explicit recognition of unconscious selection (Darwin, 1909). Darwin (1868) divided human-mediated selection into two components: methodical selection, “systematically endeavor[ing] to modify a breed according to some predetermined standard,” and unconscious selection, “that which follows from men naturally preserving the most valued and destroying the less valued individuals, without any thought of altering the breed.” Unconscious selection, he posited, was no different from natural selection. Humans change the conditions in which cultivated species live and reproduce, and this change exerts selection on the population even in the absence of a choice or predetermined goal by the cultivator. Although the term “unconscious selection” fell out of use for many years (Darlington, 1963), students of crop evolution nonetheless recognized its fundamental role in domestication. Both Vavilov (1926) and Engelbrecht (Zeven, 1973) viewed the initial stages of domestication as determined entirely by unconscious selection, and modern workers widely cite the central role of unconscious selection, sometimes referred to as automatic (Harlan, 1992) or unintentional (Rindos, 1984), in effecting observed variation in domesticated species. These authors argue that many of the phenotypic changes associated with domestication are likely to have arisen via unconscious selection and, like Darwin, view unconscious selection, and much of the process of domestication, as illustrative of the process and effects of natural selection. The phenotypic changes associated with adaptation under domestication are substantial. Many of these changes are shared across a broad array of domesticated plants; this suite of changes is commonly referred to as the “domestication syndrome” (Hammer, 1984). Common features
OCR for page 209
In the Light of Evolution, Volume I: Adaptation and Complex Design of the domestication syndrome are larger fruit or grain, reduced branching, gigantism, the loss or reduction of seed dispersal, the loss of seed dormancy, changes in photoperiod sensitivity, and the loss or reduction of toxic compounds (Hammer, 1984; Gepts, 2004). For example, the major cereal crops in Table 11.1 (rice, wheat, maize, and barley) all experienced a series of parallel phenotypic shifts brought about by domestication, including reduced seed dispersal, reduced branching or tillering, decreased seed dormancy, synchronized seed maturation, an increase in grain size, and larger inflorescences. We emphasize three salient points about the phenotypic changes associated with domestication. First, with the possible exception of characteristics such as color or fruit size that were clearly desirable by humans, most features of the domestication syndrome are likely the result of unconscious selection (Hammer, 1984; Heiser, 1988; Harlan, 1992; for a more inclusive view see Darwin, 1868). Second, the traits most clearly resulting from unconscious selection are those that would have been difficult for early cultivators to notice or that would have changed without any direct effort. Seed dormancy, for example, would be selected against by almost any method of cultivation, even without a conscious decision to plant only nondormant individuals. Finally, like its natural counterpart, unconscious selection is not limited to visible phenotypes; much of the adaptation under domestication may have involved physiological or developmental changes corresponding to the new edaphic, photosynthetic, hydrological, and competitive regimes associated with cultivation. TWO APPROACHES TO FINDING ADAPTIVE GENES Dramatic shifts in phenotype associated with domestication are not only important as evolutionary examples; they have broad economic and societal consequences. There is substantial interest in discovering the genes and genetic mechanisms that contribute to phenotypic changes associated with domestication, because their identification may facilitate trait manipulation through modified breeding strategies (McCouch, 2004). We discuss two approaches to this goal, starting at opposite ends of the phenotype–genotype hierarchy. To date, most research has followed what we call a top-down approach, which begins with a phenotype of interest and then identifies causative genomic regions via genetic analyses such as quantitative trait locus (QTL) and linkage disequilibrium (LD) mapping (Fig. 11.1 Left). An alternative approach is to build on Darwin’s view of domestication, starting with the concept of adaptation and moving from the bottom up. In this approach, population genetic methods are used to search for the signal of adaptation in a set of genes, and traditional molecular methods are used to move from gene to phenotype (Fig. 11.1 Right). Here we introduce
OCR for page 210
In the Light of Evolution, Volume I: Adaptation and Complex Design FIGURE 11.1 Schematic of the phenotype–genotype hierarchy as represented by top-down and bottom-up approaches. each of these approaches, discuss the methodologies available for their implementation, and assess their strengths and weaknesses. FROM THE TOP DOWN: QTL AND LD MAPPING To date, all of the successes at identifying genes underlying the adaptive changes during domestication have originated from top-down approaches, beginning with the phenotype and using genetic analyses to uncover genomic regions and eventually candidate genes responsible for the phenotype of interest. The most successful method for finding these genes has been QTL mapping, but association or LD methods are rapidly gaining favor in the plant genomics community. While it is beyond the scope of this article to provide a comprehensive review of QTL and LD mapping, we review some empirical findings and highlight some of the challenges of spanning the gap between phenotype and genotype. QTL Mapping Given a trait of interest, QTL mapping was the first (and is still the most widely used) method available for localizing the genetic basis of a trait (e.g., Sax, 1923). QTL mapping has led to all of the major successes in the identification and cloning of genes underlying domestication traits (Doebley et al., 2006). The best-known examples come from tomato and maize. In the mid-1980s Tanksley and coworkers (Paterson et al., 1988) initiated QTL analysis of fruit mass in a cross between wild and domesticated tomato, localizing six QTLs. With extensive mapping efforts, they were able to isolate a region encompassing the major QTL fruitweight2.2 (fw2.2). They also demonstrated the phenotypic effect of fw2.2 with transgenic analysis (Frary et al., 2000). At about the same time Doebley and coworkers
OCR for page 211
In the Light of Evolution, Volume I: Adaptation and Complex Design (Doebley et al., 1990; Doebley and Stec, 1991) mapped differences in plant architecture and plant yield between maize and its wild ancestor, teosinte. Subsequent mapping and mutation analyses led to the isolation of major genes that govern phenotypic differences between maize and teosinte, including teosinte branched1 (tb1), a gene controlling lateral branching (Doebley et al., 1995), and teosinte glume architecture (tga), which contributes to differences in inflorescence architecture (Wang et al., 2005). These successes highlight the value of the QTL approach, but the method is not without its limitations. It can, for example, be difficult to develop mapping populations for perennial, inbreeding, and vegetatively propagated crops. Thus, some of the 15 crops in Table 11.1, such as bananas and palm trees, are intractable for study by QTL approaches. It is also important to remember that the results of QTL analysis often depend on the environment (Paterson et al., 1988) as well as the parental lines used in the cross (Doebley and Stec, 1991; Li et al., 2006a). Caution is therefore warranted in interpreting the generality of QTLs, especially in cases of multiple domestication or local adaptation. There are also numerous statistical issues, the most important of which is the limited power to accurately estimate the number and size of QTLs, an observation that has become known as the Beavis effect (Beavis, 1994, 1998). Although this limitation has not proven problematic for cloning genes of large phenotypic effect, statistical power poses a major concern for more classically quantitative traits like size, weight, or yield that are likely to be determined by a larger number of QTLs of smaller phenotypic effect, and statistical concerns become even more problematic for the estimation of complex phenomena such as epistasis (Carlborg and Haley, 2004). QTL studies have provided and will continue to provide considerable utility for identifying genes and genomic regions that contribute to phenotypes of interest. Moreover, the rate at which such genes are identified will continue to increase as genomic data become available for more species; this increase is already evident in the 2006 publication year, which witnessed an explosion of the isolation of genes contributing to major phenotypic differences between domesticates and their wild ancestors. Although not solely attributable to QTL approaches, genes isolated in 2006 included two rice shattering genes (Konishi et al., 2006; Li et al., 2006b), a rice kernel color gene (Sweeney et al., 2006), a wheat shattering gene (Simons et al., 2006), and a wheat senescence gene affecting nutritional content (Uauy et al., 2006). Even so, only a handful of genes have been isolated by these approaches (Doebley et al., 2006), and the total output has been surprisingly small given both the large amount of money and human capital invested in QTL studies and the economic and societal importance of a relatively small number of plants (Table 11.1). Furthermore, the genes isolated to date are genes of very large effect, i.e., the “low hanging fruit”
OCR for page 212
In the Light of Evolution, Volume I: Adaptation and Complex Design (Doebley et al., 2006). Substantially more effort will likely be required to identify and clone genes of smaller effect. LD Mapping In the hope of overcoming some of the limitations of QTL analysis, plant researchers have moved toward LD mapping as an additional means to identify genomic regions that contribute to phenotypes. In practice, LD mapping can be separated into two types, each focusing on a different level of genetic analysis. The first, like most QTL approaches, aims to identify genome-wide variation that associates with phenotypic variation. This requires measures of genetic variability in markers representing most of the genome and tests of phenotype–genotype association for each marker. The second type of association analysis attempts to pinpoint the causative genetic mutation(s) that effect phenotype; these latter studies typically focus on variation in one or few candidate genes rather than whole genomes. The primary advantage of LD mapping is that it can rely on population samples; there is no need for crosses and the production of large numbers of progeny. This is an obvious benefit for the study of bananas, palms, or other long-lived perennial species (Table 11.1) and in general allows studies to proceed more rapidly. In addition, the population sample may contain many more informative meioses (i.e., all those that have occurred in the evolutionary history of the sample) than a traditional QTL mapping population. As a result, the phenotype of interest may be associated with a much smaller chromosomal segment than in a QTL population, in theory providing greater mapping resolution. Like QTL methods, however, there are several features of experimental design that need to be carefully considered when undertaking LD mapping. First, distinguishing true associations from statistical noise requires large sample sizes, both for statistical power and to correct for multiple tests (Long and Langley, 1999; Macdonald and Long, 2004). Even with large sample sizes, researchers may have to assume that the effects of individual mutations are additive; testing for epistatic interactions between hundreds of markers further exacerbates the problem of multiple tests (Macdonald et al., 2005). One way to reduce this problem is to test for associations between phenotypes and haplotypes (or “haplotype blocks”) rather than individual markers (Clark, 2004). But unless haplotypes can be inferred experimentally (Morrell et al., 2006), as in selfing taxa such as barley and rice (Table 11.1), the necessary computational inference of haplotypes can prove an impediment to this approach. Another design challenge is sample origin. Geographic structure or other departures from panmixis can result in spurious associations in
OCR for page 213
In the Light of Evolution, Volume I: Adaptation and Complex Design which a genotype is associated with a geographic region rather than a phenotype. This will become especially problematic for phenotypes that vary by geographic region, such as flowering time or photoperiod sensitivity. Many important crops (e.g., barley, rice, and soybean) are derived from wild populations with extensive geographic structure (Lin et al., 2001; Morrell et al., 2003; Kuroda et al., 2006). This structure is often reflected in the domesticate as well, especially in cases involving multiple independent domestications (Londo et al., 2006; Morrell and Clegg, 2007). Unfortunately, for many of the crops in Table 11.1 we have little information about the location of domestication or population structure in wild populations. The conspicuous exceptions are rice, maize, barley, and wheat, whose domestication histories are becoming more clear (Matsuoka et al., 2002; Salamini et al., 2002; Londo et al., 2006; Morrell and Clegg, 2007). Studies of human diseases (Weiss and Clark, 2002) suggest that basic research on demographic history and population structure will be crucial to the success of LD mapping in plants. The final design challenge that we will consider here is marker (usually SNP) density. LD mapping studies are very powerful when the causative mutation is genotyped (Risch and Merikangas, 1996; Long and Langley, 1999). If the causative mutation is not genotyped, it is still possible to identify association via markers that are in LD with the causative mutation. However, the extent of LD can vary dramatically among plant species (Flint-Garcia et al., 2003; Morrell et al., 2005), among genomic regions within plants (Gaut et al., 2007), and among population samples (Tenaillon et al., 2001; Ching et al., 2002; Caldwell et al., 2006). The distribution of LD is also affected by homologous gene conversion, which predominantly disrupts short-range LD patterns (Jeffreys and May, 2004; Padhukasahasram et al., 2004; Morrell et al., 2006). Study design, statistical analysis, and controlling for biological challenges such as population structure are very active areas of research (Rosenberg and Nordborg, 2006; Yu and Buckler, 2006), but several large-scale plant LD mapping studies are currently underway despite having little background information about the extent of LD and geographic structure in the populations being studied. The difficulties inherent to LD mapping are reflected in the literature. In a genome-wide association study, Aranzana et al. (2005) confirmed several Arabidopsis QTLs for flowering time and pathogen resistance but also noted a high rate of false positive associations. Workers using large wild-caught populations of Drosophila have been unable to verify associations identified in lab populations, suggesting that some results may not be replicable regardless of sample size, the number of SNPs genotyped, or the care taken in study design (Macdonald and Long, 2004; Macdonald et al., 2005). Furthermore, failure to identify an association between a candidate gene and a phenotype of interest is likely underreported.
OCR for page 214
In the Light of Evolution, Volume I: Adaptation and Complex Design Despite these drawbacks, LD mapping has had some successes. One early example successfully linked phenotypic variation in malting quality in barley to haplotype variation at the β-amylase2 gene, a locus involved in starch hydrolysis. Differences in the coding region of barley β-amylase2 affect thermostability of the enzyme (Ma et al., 2001; Clark et al., 2003), and SNP genotyping confirmed that cultivars with high malting quality and the high-thermostability enzyme share a common haplotype (Polakova et al., 2003; Malysheva-Otto and Roder, 2006). Resequencing of candidate genes has also been used in foxtail millet and rice to determine the genetic basis of waxy or sticky grains. Mutations at the waxy (granule-bound starch synthase) locus result in changes in amylose content in the endosperm, resulting in the sticky grains popular in eastern and southern Asia (Domon et al., 2002; Kawase et al., 2005; Olsen et al., 2006). LD mapping has also been used to verify associations inferred from QTL or other approaches (Thornsberry et al., 2001; Szalma et al., 2005; Balasubramanian et al., 2006; Breseghello and Sorrells, 2006b). A promising future direction for LD mapping is the use of synthetic populations derived from a relatively small number of founders (Breseghello and Sorrells, 2006a), facilitating QTL and LD mapping in a single population while minimizing complications due to population structure (Flint-Garcia et al., 2005; Breseghello and Sorrells, 2006a; Yu and Buckler, 2006). History, Adaptation, and Population Genetics Extensive work is required to isolate a candidate gene for a particular trait, but a phenotype–genotype association is no guarantee that the trait or its candidate gene has been historically important or is an adaptation. It is tempting to conclude that observable phenotypic differences are adaptive, particularly in domesticated organisms where selection is strong and the direction of selection can be surmised. However, many of the differences between domesticates and their progenitors may not be adaptive, at least from a human perspective; for example, QTLs decreasing protein content in wheat (Uauy et al., 2006) and seed size in sunflower (Burke et al., 2002) are unlikely to have been directly selected during domestication. A number of alternative processes can explain observed phenotype–genotype associations, including genetic drift, selection on a correlated trait, pleiotropy, or even natural selection working in opposition to anthropogenic selection. It therefore behooves us to endeavor to test adaptive hypotheses rather than assume them to be true (Gould and Lewontin, 1979). To understand the process of adaptation during domestication, one must first consider the genetic history associated with domestication.
OCR for page 215
In the Light of Evolution, Volume I: Adaptation and Complex Design Domestication of all plants and animals led to a reduction in genetic diversity (Ladizinsky, 1985; Doebley, 1989; Gepts, 2004), and thus all genes in any domesticated plant necessarily have a history that includes a recent demographic event, the bottleneck associated with domestication (Fig. 11.2). Population subdivision in the wild ancestor, ongoing introgression between the crop and wild relatives, and multiple domestication events can also have demographic impacts. Genes important for domestication were also subjected to conscious or unconscious directional selection, experiencing a reduction in variation over and above that associated with any demographic events (Fig. 11.2). The level of diversity remaining at a given locus in a domesticate is thus expected to be inversely proportional to the locus’s adaptive importance during domestication. Thus, the major genes contributing to agronomically important traits may lack variation entirely (Whitt et al., 2002). With a candidate gene in hand, molecular population genetic methods can be used to test adaptive hypotheses. Conceptually, the approach is simple: under the selection scenario described in Fig. 11.2, one expects that genes contributing to adaptive traits will have low genetic variation relative to nonselected genes. In addition, a strongly selected gene may have other discriminating features, such as an excess of low frequency polymorphisms or high intralocus LD (Przeworski, 2002). It is thus essen- FIGURE 11.2 Schematic representation of a population bottleneck and its effect on a neutral gene and a selected gene. In Upper, shaded circles representgenetic diversity. The bottleneck reduces diversity in neutral genes, but selection decreases diversity beyond that caused by the bottleneck alone. Lower illustrates sequence haplotypes of these two hypothetical genes. The neutral gene lost several haplotypes through the domestication bottleneck, but the selected gene is left with only one haplotype containing the selected site.
OCR for page 216
In the Light of Evolution, Volume I: Adaptation and Complex Design tial to assay genetic polymorphism in a number of randomly selected reference genes and compare them to the candidate gene. However, this is often difficult to do well, requiring many randomly sampled genes (see below) and computationally intensive simulation methods to estimate the underlying demographic model. FROM THE BOTTOM UP: MOLECULAR POPULATION GENETICS In marked contrast to the top-down or phenotype-first approaches already discussed, bottom-up approaches start by identifying genes with the signature of adaptation using population genetics and then make use of a broad array of genetic tools to identify the phenotypes to which these genes contribute. Bottom-up approaches are relatively new, and many of the methodologies are still being developed, but we believe that they have the potential to revolutionize crop genetics. Here we briefly introduce some of the methods and outline the challenges involved in identifying candidate genes using population genetics. Fitting a Demographic Model Ideally, bottom-up approaches begin by assaying genetic diversity in hundreds of loci, preferably from a sample of ≈100 individuals representing both the domesticate and its wild ancestor. Given sequence polymorphism data, several factors will affect the ability to detect the signal of adaptation, including the strength and history of selection, rates of mutation and recombination, and the demographic history of the population (Wright and Gaut, 2005). As mentioned above, demographic considerations are particularly important for crop plants, likely invalidating standard population genetic tests designed to detect the signal of selection. The standard tests typically assume that populations evolve according to the idealized Wright–Fisher model, with panmictic populations of constant population size. When these assumptions are inaccurate, as they certainly are for most domesticated species (Fig. 11.2), tests to detect selection can be wildly inaccurate. For example, computer simulations show that Tajima’s D, a commonly used test statistic for selection, identifies up to 25% of loci as selected after a change in population size due to a bottleneck, even when there has been no selection (Wright and Gaut, 2005). Another recent method incorrectly infers selection up to 90% of the time when Wright–Fisher assumptions do not hold (Jensen et al., 2005). Thus, departures from standard assumptions dramatically decrease the reliability of tests for selection and can distort the signature of selection beyond recognition (Slatkin and Wiehe, 1998; Przeworski et al., 2005). Clearly, one should view with skepticism studies of domesticated crops
OCR for page 217
In the Light of Evolution, Volume I: Adaptation and Complex Design that employ standard population genetics tests to infer selection and thus the historical importance of a gene. How, then, does one address the problem of demography? One way is to develop a demographic model that provides a reasonable fit to available data and then apply statistical tests of selection under that demographic model. The estimation of demographic history from DNA sequence data was first applied to maize (Eyre-Walker et al., 1998; Hilton and Gaut, 1998). In these early studies, sequence variation was assessed at a handful of loci from maize and its wild ancestor, teosinte. It was explicitly assumed that there had been no artificial selection on these genes. Observed polymorphism data were compared with data simulated under a historical coalescent model that included a population bottleneck. The size and duration of the bottleneck were varied via simulation, and bottleneck parameters that best fit the observed data were determined. For example, over a time frame of 2,800 years (an estimate of the duration of domestication based on archaeological data), the effective size of the population maintained through the bottleneck was estimated to be ≈2,900 individuals (Hilton and Gaut, 1998). This indicates that high genetic diversity in maize is not necessarily due to a large founding population. Taken further, such inferences can be applied to better understand the agricultural practices of early domesticators (Hillman and Davies, 1990). Since these initial studies, it has become possible to gather sequence polymorphism data from hundreds of loci. With many loci, it is no longer necessary or appropriate to assume that none of the genes has been targeted by selection, and it becomes possible both to infer the proportion of genes under selection and to identify those genes. At the same time, coalescent models have improved and can now include demographic factors such as recombination, population growth, and introgression (Hudson, 2002). The bottom-up approach should be especially powerful when applied to domesticated species, for three reasons. First, archaeological remains provide independent information about the timing of the domestication bottleneck, and its effects are relatively well understood. Second, artificial selection is strong and domestication is recent on an evolutionary time scale, so that the signature of selection should be highly detectable in patterns of genetic diversity (Przeworski, 2002; Olsen et al., 2006; Teshima et al., 2006). Third, polymorphism can be compared between a crop and its wild ancestor, greatly increasing inferential power (Voight et al., 2005) and helping to discriminate among evolutionary events before, during, or after domestication (Wright and Gaut, 2005). Examples of this demographic approach have appeared in the literature with increasing frequency (95) and have been incorporated into testing for selection in humans and Drosophila as well (Tenaillon et al., 2004). At present, however, the process
OCR for page 218
In the Light of Evolution, Volume I: Adaptation and Complex Design of estimating a demographic model is time-consuming and computationally intensive and requires substantial population genetic expertise. Thus far, the bottom-up approach to studying domestication has been applied only to maize. Wright et al. (2005) formulated a demographic model of maize domestication using sequence polymorphism data from 793 genes in 14 maize inbred lines and 16 haploid plants from its wild progenitor, teosinte. With these data, Wright et al. (2005) first sought to estimate a plausible demographic model and then asked whether the data were more likely if directional selection on a subset of loci was included in the model. Applying a novel likelihood ratio approach to this problem, they estimated that 2–4% of their loci were linked to a target of artificial selection during domestication. Their approach also allowed them to rank loci in terms of evidence for selection. The list of selected genes is enriched for functions related to transcription factors, genes implicated in plant growth, and genes involved in amino acid biosynthesis. Moreover, genes identified as targets of selection clustered nonrandomly around previously identified QTLs for domestication traits (Wright et al., 2005) and are more highly expressed than random genes only in the maize ear (K. M. Hufford and B.S.G., unpublished results), an organ expected a priori to be the target of selection. Empirical Ranking The demographic approach for finding candidate “adaptive” genes is model-intensive. As an alternative to estimating the demographic model, several studies have simply ranked genes empirically (Toomajian et al., 2006; Voight et al., 2006). This is an acceptable, but not optimal, solution based on a straightforward idea. Under the selection scenario described in Fig. 11.2, one expects that genes contributing to adaptive traits should have low genetic variation or skewed allele frequencies compared with nonselected genes (Tajima, 1989). Without knowing the exact demographic model, it thus makes sense to assay genetic polymorphism in a number of genes, compare them, and rank them by summary statistics. The candidate gene, if selected, should fall into the extreme tail of the distribution of summary statistics like S, the number of SNPs in the gene, or Tajima’s D, a measure of the allele frequency spectrum. If the gene is extreme, then the polymorphism data are consistent with an adaptive hypothesis. This idea can be applied to a genome-wide sample of genes to identify candidate genes de novo via bottom-up methods, or, alternatively, to compare a candidate gene identified by top-down approaches to a sample of reference loci. Although empirical ranking is a suitable approach, its efficacy depends greatly on the particulars of individual evolutionary histories
OCR for page 219
In the Light of Evolution, Volume I: Adaptation and Complex Design and the number of sampled loci used. Simulations show that the false discovery rate of this method may be high for recessive genes, for genes selected from standing variation, and for populations that have undergone demographic change (Teshima et al., 2006), all factors likely to have played a role in adaptive events under domestication. Similarly, although statistical methods are available to explicitly test for selection using this approach, recent results suggest that the false positive rate of these statistics is also high (but can be controlled, e.g., Thornton and Jensen, 2006). To illustrate these potential problems empirically, we describe results from resampling the data of Wright et al. (2005). We simulated scenarios in which a researcher tests for selection in a candidate gene and ranks the candidate relative to a sample of reference loci. To do this, we chose both a “candidate” locus and a “reference” sample of loci, sampling with replacement from the complete data set. For each candidate locus we asked whether it was extreme in terms of low genetic diversity (measured by S) or in its frequency distribution (measured by Tajima’s D). We made two comparisons for each locus, first using the observed measure in maize, then using the difference between maize and teosinte. We treated each of the ≈800 loci as a candidate locus and compared it to a random set of reference loci. We initially set the size of the reference sample to five loci, then repeated the experiment with samples of 10, 20, and 50. The entire process was repeated 1,000 times, giving an estimate of the distribution of numbers of candidate genes rejected for each reference sample size. Results of this resampling are shown in Fig. 11.3, presented for both summary statistics and a scenario in which the neutrality of a candidate locus is rejected if either of the statistics is extreme. Two valuable insights emerge from this exercise. The first is that using small samples of reference loci gives very poor results. With a sample of five loci for comparison, our simulations rejected >12% of the candidate loci for each of the statistics and rejected ≈35% of the loci if the most extreme of the two statistics was used. Samples of size 10 improve the situation, but our simulations still rejected more than twice as many loci as the model-based approach (Wright et al., 2005). The second insight is that including ancestral individuals improves the test greatly. Because loci in domesticated plants have undergone a demographic bottleneck, a substantial number of loci have lost much or all of their variation, even in the absence of selection: of the 793 loci in the data set, 90 have no diversity. Using only observed diversity in maize for comparison, these 90 loci are always rejected as extreme, regardless of the sample size. When comparative data from teosinte are included it becomes clear that some of the loci with zero diversity in maize have low levels of diversity in teosinte as well, suggesting that these are low-diversity genes rather than genes disproportionately affected by domestication.
OCR for page 220
In the Light of Evolution, Volume I: Adaptation and Complex Design FIGURE 11.3 Resampling tests to examine empirical ranking methods for finding candidate genes. Left represents sequence data from maize alone; Right demonstrates the difference in statistics between maize and teosinte. The statistics are the number of SNPs (S) (Top), Tajima’s D (Middle), and a combination of both (Bottom) (see text). For each graph, the heavy line represents the median number of genes, of ≈800, that are inferred to be under selection. Boxes represent the central 50% of the data, and lines extend out to 3/2 of the interquantile range. From Gene to Phenotype via Molecular Genetics The biggest drawback to bottom-up approaches is that the candidate genes are not associated with a phenotype. In theory this can be rectified by using the array of genetic tools available for many model species. For the first five species listed in Table 11.1, for example, a number of genomics tools exist to aid in connecting a candidate gene to a phenotype: databases for ESTs, microarray, and gene expression data; targeted mutagenesis lines; genetic maps; and partial or complete genome sequences. For model crops such as maize and rice, reverse genetics methodologies have been transformed into high-throughput pipelines suitable for the analysis of large numbers of genes (McCallum et al., 2000). The link from gene to phenotype for nonmodel species, however, may be daunting. Although high-throughput analysis of phenotype is a distant possibility for many species in Table 11.1, bioinformatics and standard reverse genetics techniques can still provide a bounty of information regarding possible phenotypes. For many crops without extensive genetic resources, valuable information can nonetheless be gleaned from comparative
OCR for page 221
In the Light of Evolution, Volume I: Adaptation and Complex Design genomic analysis of gene function or expression in related species. At worst, comparative bioinformatics can provide insight into the general class of gene, potentially offering information about the role the gene has played during domestication. Similarly, many standard reverse genetics approaches such as RNA interference and transgenic methods can lead to significant clues as to gene function. PERSPECTIVE: TOP-DOWN VS. BOTTOM-UP Both top-down and bottom-up approaches will continue to prove useful for the study of adaptation to domestication. With the current rate of increase of genomic information for many crop species, we expect that the dramatic increase in top-down success stories seen in 2006 will continue for some time, identifying some of the genes of large effect that contribute to the phenotypes associated with domestication. Anytime the goal is to identify genes underlying a specific phenotype of interest, these top-down approaches will continue to be the best choice. We argue, however, that top-down approaches have a severe and insuperable limitation for the study of adaptation: the requirement of identifying a phenotype a priori. It is plausible (and even likely) that alleles influencing fruit size in tomato (Nesbitt and Tanksley, 2002) or inflorescence structure in maize (Bomblies and Doebley, 2005) have evolved as adaptations to domestication; the available genetic evidence does not speak to a history of selection, and without comparative genetic or experimental evidence of selection, adaptive hypotheses for these genes must remain, at best, hypotheses. We must be careful not to assume adaptation simply because a gene correlates with a trait of agronomic importance, and the converse is equally true: there are likely many genes that, although not responsible for obvious morphological change, will nonetheless show evidence of selection and adaptation under domestication. Furthermore, the genetic history of crops creates a dilemma for QTL and LD mapping approaches: mapping requires segregating genetic diversity at the gene of interest, but genes governing historically important phenotypes are expected to have low genetic diversity in the domesticate (Fig. 11.2). QTL and LD approaches that do not include wild populations are likely to miss many of the genes contributing to agronomic traits that were important during early domestication. In contrast to QTL studies that can use wild × domesticate crosses, LD mapping is faced with a Catch-22: including both wild and domesticated individuals will lead to spurious associations due to sample origin, but purely wild populations will be depauperate for the domesticated phenotype (and genotype) of interest. In addition to their freedom from the constraints of a priori phenotypic choices, bottom-up approaches have several advantages for finding genes
OCR for page 222
In the Light of Evolution, Volume I: Adaptation and Complex Design that contribute to adaptive traits and that will be useful in an agronomic context. These advantages include the following: (i) segregating variation is not required to identify genes of interest; (ii) far fewer plant samples are needed than for LD mapping, with only tens (<100) of samples (Teshima et al., 2006) often sufficing as opposed to hundreds or thousands (Long and Langley, 1999); (iii) like LD mapping, bottom-up approaches can be applied to species that reproduce slowly and lack genetic tools; and (iv) they allow inferences about demographic history, providing historical insights into the process of domestication. We should note, of course, that bottom-up and top-down approaches are not mutually exclusive; for example, bottom-up approaches in maize are also being used to identify candidate genes for LD mapping (Yu and Buckler, 2006). Although they have advantages, bottom-up approaches are also not a panacea, for at least four reasons. First, their success will vary among species, depending on levels and distribution of genetic diversity. For example, initial surveys of genetic diversity in sorghum have failed to identify selected genes (Hamblin et al., 2006). This failure is in part a limitation of the study system, because sorghum has low genetic diversity, but it may also reflect inefficient sample design. Simulation studies suggest that these methods should be quite powerful with moderate (<100) sample sizes, even with diversity levels as low as those found in sorghum (Teshima et al., 2006, and M. Przeworski, personal communication), but empirical studies relied on a sample of 17 domesticated individuals and only one wild plant (Hamblin et al., 2006). Second, genes identified as selected may have been targets of selection or may be linked to a target of selection (through “hitchhiking”). For example, selection on the rice gene waxy appears to have affected patterns of sequence diversity in 29 additional genes. This lack of resolution is, however, a shortcoming shared with QTL and LD methods, because in all cases it is difficult to differentiate between a target (or “causal”) marker and linkage effects (Weigel and Nordborg, 2005). In fact, when the genomic locations of genes are available, the expected chromosomal resolution of bottom-up approaches is at worst similar to QTL and LD mapping. Third, like top-down approaches, bottom-up approaches may not be feasible for all crops. The limitation here is not generation time (as in QTL studies), but rather levels of genetic diversity, polyploidy, and population structure. Polyploidy makes population genetic analysis difficult, requiring careful separation of homeologs and their independent evolutionary histories. As with LD mapping, unrecognized population structure can be problematic for population genetic analyses, producing patterns that can be mistaken for selection. Finally, bottom-up approaches share a major limitation with both association and QTL mapping. All three methods identify candidate genes or regions, but verification requires additional functional characterization (Weigel and Nordborg, 2005). It is worth noting that in many cases
OCR for page 223
In the Light of Evolution, Volume I: Adaptation and Complex Design this last step, connecting a candidate gene to a phenotype via functional studies, is often not much easier for top-down approaches than for candidate genes identified by using population genetics. Over the past 25 years, top-down approaches have yielded a list of ≈30 genes with well characterized phenotypic effects in plants (Doebley et al., 2006). It is known that these 30 are genes of major effect, i.e., either Mendelian factors or major QTLs, but for most it has not been determined whether they have played an important adaptive role historically. In contrast, limited application of bottom-up approaches in maize have identified ≈50 genes with a signature of adaptation. It is a statistical certainty that some of these will prove to be false positives, but it is also likely that some of these genes contribute to phenotypes that would not or could not be studied via QTL or LD mapping. In the last year, the number of published, large-scale studies seeking to identify selected genes has exploded. Screens for selection have been applied to polymorphism data from humans (Bustamante et al., 2005; Voight et al., 2006) and Arabidopsis (Toomajian et al., 2006) as well as maize (Wright et al., 2005; Yamasaki et al., 2005). To a much more limited extent bottom-up approaches are being applied to other domesticated species, i.e., rice (Olsen et al., 2006; Zhu et al., 2007), flax (Allaby et al., 2005), sorghum (Hamblin et al., 2006), and dogs (Pollinger et al., 2005). We argue that there is an opportunity, in fact, a pressing need, for a broad-based initiative to implement bottom-up approaches in 15–20 important crops, not unlike a multispecies HapMap project. Such an initiative would be relatively inexpensive given new sequencing technologies and would have far-reaching consequences beyond identifying candidate genes. Important side benefits would include broader-based information on LD, SNP discovery on a panel of sufficient size to limit ascertainment biases, and evolutionary analyses of polymorphism in a genomic context. Data compared across species may also provide insights into the process of adaptation. For example, such data could inform the age-old question as to whether parallel phenotypic changes, such as the domestication syndrome, evolve via parallel genetic mechanisms (Paterson et al., 1995). Wide-scale implementation of bottom-up approaches across species would be of potential agronomic benefit, but would also provide a unique opportunity to identify the genetic basis of adaptation. ACKNOWLEDGMENTS We thank S. J. Macdonald and J. G. Waines for helpful discussion, and we thank two anonymous reviewers. This work was supported by National Science Foundation Grants DEB-0426166 (to B.S.G.), DBI-0321467 (to B.S.G.), and DEB-0129247 (to M. T. Clegg).
OCR for page 224
In the Light of Evolution, Volume I: Adaptation and Complex Design This page intentionally left blank.