doi:10.1111/j.1420-9101.2007.01410.x MINI REVIEW Neutral theory: a historical perspective E. G. LEIGH JR Smithsonian Tropical Research Institute, Balboa, Panama Keywords: forest ecology; neutral theory; null hypotheses; population genetics. Abstract To resolve a panselectionist paradox, the population geneticist Kimura invented a neutral theory, where each gene is equally likely to enter the next generation whatever its allelic type. To learn what could be explained without invoking Darwinian adaptive divergence, Hubbell devised a similar neutral theory for forest ecology, assuming each tree is equally likely to reproduce whatever its species. In both theories, some predictions worked; neither theory proved universally true. Simple assumptions allow neutral theorists to treat many subjects still immune to more realistic theory. Ecologists exploit far fewer of these possibilities than population geneticists, focussing instead on species abundance distributions, where their predictions work best, but most closely match non-neutral predictions. Neutral theory cannot explain adaptive divergence or ecosystem function, which ecologists must understand. By addressing new topics and predicting changes in time, however, ecological neutral theory can provide probing null hypotheses and stimulate more realistic theory. Introduction To resolve what he thought was a paradox in panselec- tionism, Kimura (1968a) proposed a neutral theory of population genetics. Here, all genotypes have equal fitness: the only causes of change in al?ele frequencies are mutation, migration and demographic stochasticity. This theory provided a unifying frame for earlier studies of Wright (1931), Mal?cot (1948), Kimura (1955), Karlin & McGregor (1964, 1967) and many others. It evoked hostility (Ohta & Gillespie, 1996), but also stimulated elegant and creative mathematics (Ewens, 1972; Karlin & McGregor, 1972; Nagylaki, 1974, 1976; Watterson, I974a,b, 1976,1984;Kingman, 1977, 1978,1982; Sawyer, I977a,b, 1979; Sawyer & Fleischmann, 1979; Donnelly, 1986; Donnelly &Tavar?, 1987, 1995; Hoppe, 1987). To learn whether species-area curves and species abundance distributions can be explained without invok- ing differences among species, Hubbell (1979, 1997, 2001), following an earlier attempt by Caswell (1976), constructed a neutral theory of forest dynamics and diversity. Here, a tree's species does not affect its prospects of death or reproduction. Hubbell modelled Correspondence: Egbert Giles Leigh Jr, Smithsonian Tropical Research Institute, Unit 0948, APO AA 34002-0948, USA. Tel.: -^507 212 8940; fax: 4-507 212 8937; e-mail: leighl@si.edu this theory on the neutral population genetics of a single multi-allelic locus in a haploid population. Neutral theory rewarded Hubbell, as it had Kimura, with abun- dant controversy, and with a variety of quantitative predictions, some in accord with observation. Theories that included more biology seemed unable to match these neutral predictions in an unforced manner. Hubbell (2001), however, was a best-seller, whereas Kimura (1983) was not, even though Kimura's theory may prove more useful (Nee, 2005). Why were their receptions so different? What does this difference reveal about ecology compared with population genetics, and about the prospective contribution of neutral theory to each? To answer, I will: 1 review selected questions posed about neutral pro- cesses in population genetics, including Kimura's neutral theory, and the reasons for posing them; 2 outline those techniques for studying neutral processes in population genetics that allow the clearest and simplest presentation of Hubbell's neutral theory; 3 outline the motives and achievements of Hubbell's (2001) neutral theory of forest ecology, with special attention to its successful predictions; 4 discuss criticisms of Hubbell's neutral theory, evaluate its weaknesses and strengths, and suggest explanations for the contrasting reception of neutral theory in population genetics and ecology; ? 2007 THE AUTHOR. J. EVOL. BIOL. 20 (2007) 2075-2091 JOURNAL COMPILATION ? 2007 EUROPEAN SOCIETY FOR EVOLUTIONARY BIOLOGY 2075 2076 E. G. LEIGH JR 5 suggest ways to expand this theory's scope and effectiveness, and outiine its potential contributions to forest ecology. Neutral theory in population genetics: questions the This section outlines investigations concerning neutral processes in genetics from 1922 to their culmination in Kimura's (1968a, 1983) neutral theory. What questions were asked? Why were they asked? The presentation will emphasize questions relevant to the neutral theory of forest ecology; so, it focuses on inflnite-allele models for a single locus in a haploid population. Despite the extensive study of neutral processes and the enormous body of techniques developed for this study, no population genet- icist, not even Kimura, sought to deny the importance of adaptive evolution. Instead, all major workers were inter- ested, at least to some degree, in how neutral processes affected adaptive evolution. Wright (1932) and Kimura (1983) even asked how neutral processes might rearrange genetic variation in ways that facilitate adaptive evolution. Ronald Fisher asked questions, and developed tech- niques, relevant to the neutral theory to understand the origin of adaptation more clearly. Fisher (1922, 1930:73- 76) used generating functions to calculate the prospects of survival and spread of a neutral mutation, because chance is the decisive influence on an allele's prospects when its numbers are few. Fisher (1958:84-87) also used gener- ating functions to show that if a neutral mutant with a single bearer at generation 0 survived to generation n, it would have probabOity exp{-2jin) of being represented then by more than j individuals. He concluded that chance would take unrealistically long to spread a neutral mutant through a large population. Fisher (1922, 1930) also showed that, absent mutation and selection, genetic variation in large populations declines so slowly that mutation is too rare to dictate the direction of evolution. Finafly, Fisher (1930) argued that sexual reproduction evolved to enable the simultaneous fixation of new favourable mutations that first occur in different individ- uals of a population. This explanation presupposes finite populations and the role of chance in the spread of new mutations, but it focuses on adaptation. Fisher's studies of neutral processes were essential parts of his theory of adaptation by natural selection. Sewafl Wright also investigated neutral processes to understand the origin of adaptation. Wright (1932) believed that a genotype's fitness could not be accurately estimated from the average fitnesses of its component genes. He concluded that many populations would have to cross an 'adaptive vaUey' of lower fitness to attain higher fitness. Such crossings could only occur by chance. He therefore believed that adaptive evolution occurred most readily in large populations divided into many subpopu- lations that exchanged occasional migrants. In each subpopulation, afiele frequencies would undergo 'genetic drift', leading to the formation of new genotypes. When an especially favourable genotype spread through a subpopu- lation, migrants would export it to the rest of the species (Wright, 1932). Wright (1931) therefore studied neutral processes to learn how they affect the occurrence of new gene combinations. To this end, he calculated the prob- ability density of possible frequencies of a neutral al?ele in a subpopulation where immigration and mutation are in balance with genetic drift. This calculation enabled him to learn how low the rate m of exchange of migrants must be to allow significant differentiation among subpopulations (Wright, 1931:127-128). Although population geneticists including Kimura (1955) had already treated a variety of neutral processes, Kimura only proposed his neutral theory in 1968 (Kimura, 1968a). He did not deny adaptation. Rather, he invented the neutral theory to circumvent an appar- ent paradox. Haldane (1957) had calculated that the spread by natural selection of an al?ele of initial fre- quency q(0) through a population of constant size N would cost N\n[Hij{0)] selective deaths that would not have occurred were this al?ele substituted instantly. Kimura (1968a, 1983:25-26) proposed his neutral theory when he learned that a population experiences a nucle- otide substitution every 2 years: he thought no popula- tion could support the mortality that these substitutions would impose were they driven by selection. Inventing the neutral theory, however, brought surprising rewards. In 1969, Kimura told me that the neutral theory was the first theory to allow him to make interesting, testable, seemingly valid, predictions in population genetics. Lewontin ( 1974) also noticed the difficulties of extracting useful predictions from non-neutral population genetics, but he did not become a neutral theorist. Kimura's primary achievement in neutral theory was to show that a population's rate of neutral evolution is independent of its size. In a population of JV adults, a new mutation carried by a single adult has frequency UN. Its 'expected' future frequency is UN: it either disappears, with probability 1 - 1/JV, or takes over the population, with probabOity UN. If u new mutations occur per individual, Nu occur per generation of N individuals, of which u wOl be fixed. Therefore, the mutation rate per individual is the substitution rate per population (Kimura, 1968a ). This achievement was the theoretical basis for the 'molecular clock' (Zuckerkandl & Pauling, 1965). Where the rate of molecular evolution is constant (clock-like), an easOy understood, easily appUed technique for inferring molecular phylogenies, the unweighted pair group meth- od with arithmetic mean (UPGMA) is trustworthy (Nei, 1987:293-298; Felsenstein, 2004:161-166). The 'molecu- lar clock' was a principal target for critics of the neutral theory (GiUespie, 1986). Kimura & Crow (1964) considered the balance between mutation and random extinction at a single locus in a population of N haploids, where successive generations are distinct, each adult has probability u of carrying a new ? 2007 THE AUTHOR. J. EVOL. BIOL. 20 (2007) 2075-2091 JOURNAL COMPILATION ? 2007 EUROPEAN SOCIETY FOR EVOLUTIONARY BIOLOGY The neutral theory: a historical perspective 2077 mutation and each new mutation is a new al?ele. Here, the probability F that two adults, sampled with replacement, carry the same al?ele is 1/(1 + 2Nu). Moreover, the expected number of al?eles with frequency between c? and q + dq is 4>(q)?q = 2Nu(l - qf'^"~''dq/q (Kimura & Crow, 1964:731, Kimura, 1968b:252). Finally, Kimura & Ohta (1969:765) showed that in a population of N haploid adults, a new immutable mutant now carried by a single adult would take an average of 2JV generations to spread through the population, if it is lucky enough to do so. If it is now carried by ; adults, it takes an average of 2{N?j)(N - j)\n[N?{N - j)] genera- tions to spread through the population, if it is lucky enough to do so. Since, Kimura's neutral theory has proved valid for 'third position' nucleotides with no effect on amino acid sequences (Ohta & Gillespie, 1996). It has also become a source of increasingly subtle null hypotheses by which to test for the presence of natural selection (Hartl & Clark, 1997:316; Nielsen, 2001). Subjects of these tests range from genome-wide phenomena (Bernardi, 2007) to the interplay of linkage, mutation and random drift at nearby sites along a DNA sequence (Sabeti et al, 2002). Neutral theory in population genetics: useful techniques This section outlines techniques that simplify the analysis of neutral processes at single loci in haploid populations, first studied by other methods. I illustrate techniques inspired by Moran (1962) and Mal?cot (1948) that allow clearer presentation of the basic theses of the neutral theory of forest ecology and make its further development easier. Models of populations with overlapping generations, especially Moran's (Moran, 1962:78-81; Ewens, 2004:104-109), are much easier to solve than models where successive generations are distinct. To see how such models simplify neutral theory, consider a multi-allelic locus in a panmictic population with N haploid adults. There are ?V time-steps per generation (the time required for iVadults to die). Every time-step, let an adult be chosen at random to die, and another, which might (in contrast to the model of Hubbell, 1979) be the same one, to provide the gene at our locus for the instantly maturing replace- ment (hereafter, we say, 'provide the replacement young') for the dying adult. Finally, let each adult have probability M < 1 of bearing a new al?ele. probability that a new al?ele is born at that time-step is u. At steady state, u ^ 4)(\)(N - \)IN^, 4)(\)^uN^I (N - 1 ). The probability PijJ - 1 ) that one of the j bearers of an al?ele dies at the next time-step and is replaced by the bearer of another al?ele is nearly j(N - j)?N^; the proba- bility Pij - \,j) that one of an allele's j - 1 bearers produces a nonmutant young then, replacing a bearer of another al?ele, is (j - 1){N -j + 1)(1 - u)/N^. At steady state, the number 0(y - l)P(j - l,j) of al?eles with mem- bership increasing from 7 - 1 to j balances the number (j) (J)P(JJ ~ 1 ) with membership decreasing from j to j - 1. As^(j- l)Pij-l,i) = 4>ij)Pij,i- I), ^(j) = p(j-i,M(j-i)/p(j,j-i) ^ 0?- 1)?- l)(W-i+ l){l-u)/[j{N-j)], 0(2)^(j) = Ne(l - uy/[j(N-j)]. Let j < JV (almost always true if ? > 1) and set j/N = x, UN = q, 4>(j)IN = :>{q)dq = e(l - uf'dq/[qil - q)] ^ 6(1 - qf-'dq/q, = x= NI{N + 6) (Kimura & Crow, 1964:731). If I - u and 9 = a, then ^(/?)~0(i-M)Vi = =?Vi, (1) the log-series of Fisher et al. (1943). In the Moran model, 6 = uNI(l - u) corresponds to Hubbell's (2001:121) fundamental biodiversity number. The probability F that two genes are of the same al?ele, and its uses Mal?cot (1948, 1969) showed that the behaviour of the probability F{t) that two adults sampled from the same haploid population at generation t carry the same al?ele at a given locus is easily predicted. First, he showed how F{t) changes in a population where al?eles do not mutate (u = 0), successive generations are distinct and each generation has N reproductives. With probability 1/JV, two adults sampled from generation ? -1- 1 inherit their gene from the same parent, so F = 1. With probability I - \IN their genes come from different parents, so they have probability F{t) of being the same al?ele. Therefore, lUloran and the steady-state distribution of genes over al?eles Let 4)ij) al?eles be represented by j genes apiece. What is the value of (?!)(/) at steady state? The probability P( 1,0) that the only carrier of its al?ele dies at the next time-step, and either does not reproduce, or replaces itself with a mutant young, is \(N - 1 -1- u)IN^ ^ (N - 1)?N^, whereas the F{t+l) = l/N+{l-l/N)F{t), 1-F{t+1) = {l-l/N)[l-F{t)]. l-F{t) = {l-l/N)[l-F{t-l)] = (1 - l/N)'[l - F{0)] ^[l- F{0)]exp{-t/N). (cf. Mal?cot, 1948:32; Crow & Kimura, 1970:101-102). 'Genetic drift' diminishes genetic variation, as measured ? 2007 THE AUTHOR. J. EVOL. BIOL. 20 (2007) 2075-2091 JOURNAL COMPILATION ? 2007 EUROPEAN SOCIETY FOR EVOLUTIONARY BIOLOGY 2078 E. G. LEIGH JR by 1 - -F, but does so very slowly, at a rate inversely proportional to population size, as Fisher (1922, 1930) first showed. Now let generations overlap according to the Moran model. Here, one adult dies at each time-step (each UN generation) and is immediately replaced by the instantly maturing young of an adult randomly chosen from the N alive just before the aforementioned adult died. Sample two adults at time t. With probability 1 - 2IN, neither dies by time ? + \IN. If one dies, and if al?eles do not mutate, its replacement has probability 1 IN of inheriting its gene from the survivor. Here, F(t + l/N) = F(i)(l - l/N) + (2/N)[\/N + [l-l/N)F{t)]=F{t) + 2[l-F{t)]/N\ l-F(t+ l/N) = [1 - F(t)][\ - l/N^] = [\-F{\/N)][\-2/N^f' Ri [1 -?'(l/]V)]exp(-2?/]V) (cf. Leigh et al, 1993, 2004a; Gilbert et al, 2006). Genetic drift diminishes 1 - _F twice as fast in populations where generations overlap. Now let each adult have probability u of bearing a new mutation at our locus. If generations overlap according to the Moran model, what is the equilibrium probability F that two adults sampled from a population of N bear the same al?ele? These adults' genes must have a common ancestor, and F is the probability that no mutation occurred in the line of descent of either gene since their latest common ancestor. Sample two adults at time t. The probability that one descended unmutated from the other at the last time-step is (1 - u)/2N^, and the probability that one arose as a mutant then is 2u/N. If neither happened then, the probabilities of these two events are the same for the time-step next before, and so forth. The probability -F that these genes' latest common ancestor lived after the last mutation in their ancestry is thus (Ewens, 2004:329, 339-340) [2{l-u)/N^]/[2u/N + 2{l-u)/N^] = {l-u)/{Nu+l-u) = l/[l+Nu/{l-u)]. If we set Nul(l - u) = e (Ewens, 2004:294), then F = 1/(1 -I- 6). If successive generations are distinct, the probability that one of our adults is a new mutant is nearly 2M, the probability that their al?eles are inherited unmutated from the same parent is (1 - u)^IN, and F^{l/N)/{2u+l/N) = l/{l+2Nu) (2) (cf. Mal?cot, 1948:34-35; Kimura & Crow, 1964; Gillespie, 2004:46-47). Here, as in Hubbell (2001:121), the funda- mental biodiversity number ? is 2Nu and F = 1/(1 -i- ?). Next, what is the average number k(n) of al?eles in a sample of n adults? Again, let generations overlap. Moran style. Then k{2) = 1-1-6/(1-1- 6), for 0/(1 -i- 6) = 1 - F is the probability that the second adult's al?ele differs from the first's. If n adults are sampled at one time-step, the probability that, at the last time-step, one inherited its al?ele unmutated from one of the n - I others is n{n - 1)(1 - u)IN^, whereas the probability that one was born as a mutant then is nu/N. Reasoning as before, the probability that a new mutant arose since the most recent common ancestor of any pair of our al?eles, which is k(n) - k(n-l), is (nulN)l[nulN + n{n - \)(\ - u)l JV^] = e/[l -I-(H - 1)0] (Ewens, 2004:306-307). Thus k{n), the expected number of al?eles in a sample of n adults, is 0/0 + 0/(0 + 1) + 0/(0 + 2) + ? ? ? + 0/(0 + n-l) (Ewens, 2004:338-340). If n is large, k(n) R; 0 ln[(0 -i- n)l 0] = 0 ln[l -I- K/0]. If we set 0 = ot, the relation between the number n of adults sampled and the expected number k of al?eles among them is predicted by the log-series (Fisher ei a/., 1943; Watterson, 1974a). Few species consist of one panmictic population. What happens in populations that exchange migrants with one or more others? First, consider a population of N adults where u = 0. Let each adult have probability m of being an immigrant from a large, panmictic source pool (as in Hubbell, 2001:83-86). Let the probability F* that two members of the source pool carry the same al?ele be constant. F* is also the probability that an immigrant has the same al?ele as a random adult in our population (Leigh et al, 1993). At steady state, where the probability F that two adults sampled from our population have the same al?ele does not change, F{t + UN) = F(t) is F= (1 -2/N)F + {2/N){mF* + (1 - m)[llN + (1 - 1/W)F]}, F = mF* + {\- m)/N + (1 - m){\ - \/N)F, NmF = NmF* + 1 - m - (1 - m)F, F = {NmF* + 1 - m)l(Nm +\ -m), \-F = Nm{l - F*)/{Nm +l-m). (3) Populations with fewer immigrants (lower Nm) are more homozygous. The probability that in a widespread, uniform popula- tion, two adults r km apart carry the same al?ele has long fascinated theorists (Wright, 1943, 1946; Mal?cot, 1948; Kimura & Weiss, 1964). To illustrate the method as simply as possible, consider a linear series of populations with N haploid adults apiece, which exchange migrants only with adjacent populations. Let each adult have probability u of carrying a new al?ele. How does the probability F(?) that two adults from different popula- tions separated by ? - 1 intervening populations decline as n increases? Let generations overlap. Let each adult ? 2007 THE AUTHOR. J. EVOL. BIOL. 20 (2007) 2075-2091 JOURNAL COMPILATION ? 2007 EUROPEAN SOCIETY FOR EVOLUTIONARY BIOLOGY The neutral theory: a historical perspective 2079 have probability 1 - m of inheriting its gene from its own population, probability mil of inheriting it from the population to its right, and equal probability of inheriting it from the population to its left. At steady state, let F{Q) be the probability that two adults from the same population carry the same al?ele, and let the probability that any two adults carry the same al?ele depend only on the distance between them. If n > 1, the probability that two adults n populations apart carry the same al?ele is F{n) = (1 - 2di)P(?) + 2di(l - u){m[F{n + 1) + F(M- l)]/2 + F(n)(l -m)], F{n) = (1 - u){F{n) + [F{n + 1) - F{n)]m/2 -[F{n)-F{n- l)]m/2}. Set F{n+1] ? F(n) ^ dF/dn, [F{n + 1 ) - F{n)] ?F{n-l)]^d^F/dn\ [F{n) Then d^F/dn^ ^ 2uF/m, F{n) ^ F{0)ex^\-ny'{2ulm)]. (6) The negative square root was chosen because F must decrease as n increases. Showing that F(0)=\l [l+N^(2um)} (Chave & Leigh, 2002) is too difficult to do here. This problem, or a close analogue, was solved for diploid populations where successive generations are distinct by Mal?cot (1969:81-84), Nagylaki (1974, 1976) and Sawyer (1977b). The two-dimensional case is more difficult to solve, but the methods for deriving the equation for F are the same (Chave & Leigh, 2002). /d,(s) = di+(l -2di)s + dte^ =s + dt{l-sf, ht\ft{s)\ =f,+?t{s) =m+dt[i -m]\f,+i,{s) -m = dm = dt[i - m\, di = d{i/[i -mf}, ? = i/[i -m\ -1/(1 -s); ft{s) = [t + s{l-t)\l[l+t-st]. Here, /,(0) = i/(l + t) is the probability that A has died out by time ? (Feller, 1968:480). Set po(i) =Et = i/(l + t). Then the probability that A still survives at time t is 1-E, = l/(l+i) We may now set /,(s) equal to (4) \Ft + {l- 2Et)s\/[l - EtS] =E, + {1- 2Et)s + E^s + E,(1 -2Ey+E?s^ + --- = ?/(l + ?) + [1/(1 + t)f\s + teV(l +t) + iV/(l + tf + iV/(l+?)' + ??? For every t, the expected number of A's is 1, so the expected number of A, given that the al?ele survives, is 1/[1 - Po(i)] = 1 + t- The probability Pk{t) that there are k A's at time t is i*"'/(l + i)*""- If A's stiU survive, the probability is i*~'/(l + i)* that there are k A's then, and the probability that there are at least k A's then is [i/(l + ?)]*"'- When i> 1, [?/(l+?)r'^exp(-V?)- (5) If a neutral mutant survives so long, an average of k generations elapse before there are k A-bearers. If successive generations are distinct, however, the proba- bifity is exp( - 2klt) that a neutral mutant surviving to generation t has over k bearers then (Fisher, 1930). Branching processes, Moran style Fisher (1930, 1958) used branching processes to study the initial spread of new mutations. He assumed that successive generations are distinct. Branching processes are simpler where an allele's bearers die one at a time. Consider a large population where an adult living at time ? has probability di of dying by time t + dt and equal probability of bearing an instantly maturing young by then (so time is measured in generations). If one adult carries an al?ele A at time 0 and A's do not mutate, what is the probabifity Pk{t) that there are k A- bearers at time i? To answer, form the generating function ??S) = Y,Pk{t)s As different adults reproduce independently, /,+?,(s)= ft\fw(s)]- As/o(s) = s (there is one A when t = 0), Neutral theory of forest ecology: achievements Having discussed the motives, techniques and achieve- ments for studying neutral processes in population gene- tics, let us turn to HubbeU's (1979, 1997, 2001) neutral theory of forest ecology. The questions motivating this theory are: must we know how different species coexist to predict species abundance distributions or species-area curves (Hubbefi, 2001: 6)? Are a plot's tree species present because they are the ones most suited to its microhabitats, or because they are the ones whose seeds happened to get there (niche assembly vs. dispersal assembly: Hubbell, 2001:8)? Here, there are two possibOities. Do too few seeds of a given species disperse from a large panmictic source pool to this plot to reach the microhabitats that best suit them, even though any seed could disperse across the source pool (Etienne & Alonso, 2005)? Or are seeds unable to reach the best microhabitats because they disperse only a limited distance from their parents? ? 2007 THE AUTHOR. J. EVOL. BIOL. 20 (2007) 2075-2091 JOURNAL COMPILATION ? 2007 EUROPEAN SOCIETY FOR EVOLUTIONARY BIOLOGY 2080 E. G. LEIGH JR Hubbell (2001:6-7) also considered the neutral theory as a 'teething ring' for theory with more biological content. Only recently did Hubbell (2006) suggest using the neutral theory as a null hypothesis. Hubbell (2001) thought that the neutral theory explained species-area curves. His chapter on species- area curves, however, provided a quantitative prediction only for areas so small that panmixia applies. The log- series implies that the number S of species in a sample of JV trees is 5 = a ln(l + JV/ot), where a is Fisher's alpha (Fisher et al, 1943). On some tropical forest plots < 50 ha in area, this is nearly true (Condit et al, 1996a, 2004). Caswell (1976) was the first ecologist to predict this relation from the neutral theory. Bramson et al. (1996, 1998) essayed a neutral theory of species-area curves on a square lattice of trees, each a unit distance from its four nearest neighbours. Here, a tree alive at time t has probability di of dying by time t + at. If it dies, each of its four nearest neighbours has probability 1/4 of being seed-parent of the replacement young, which has probability u of being a new species. Bramson et al. (1998) scale N relative to \lu, setting N= (Huf, In N = y\n{l/u). If M < 1, W > 1, JV? < 1, then the number of species in a square of N contiguous trees is given by S/N = u\y\ri{\lu)f/2it = ?(In Nf/ln; S = Nu{\nNfl2n, (Bramson et al, 1998, equation 1.3). If the mean square dispersal distance of offspring from parent is a^ rather than 1, then S = Nu(\n Nftlna^ (Chave et al, 2002:20). Bramson ?tal. (1996) scale JV relative to L^, setting JV = L^^ where L = 1 /VM, the distance a random-walking line of descent with unit step-length is likely to propagate across a lattice in I lu time-steps before it mutates (Hubbell, 2001:174). If JV > \lu, then 5 = 2JVM[ln(l/v^M)]Vji = JVM[ln(l/M)]V27i (Bramson et al, 1996, equation 1.2). For L so large that a species rarely spreads so far before dying out, S increases linearly with JV. No one, however, has yet derived the familiar species-individual curve, 5 = cJV^, where z < 0.4, (MacArthur & Wilson, 1967), from the neutral theory. Although 50-ha plots have stimulated analyses of species-area curves (Condit et al, 1996a, 2004; Plotkin et al, 2000a,b), species-area curves have played a minor role in empirical discussions of the neutral theory, presumably because they offer much less support for it than, say, species abundance distributions. Accordingly, species-area curves will not be discussed further in this paper. The most famous achievement of Hubbell's (1997, 2001) neutral theory is accounting for the distribution of trees > 10 cm dbh over species on a 50-ha plot at Barro Colorado Island, Panama, assuming that a fraction m of its seed-parents were chosen from a large source pool. and assuming panmixia and neutral dynamics for both the plot and the source pool. This example will be examined in detail. With neutral dynamics, the number of species repre- sented by k trees apiece in a panmictic source pool with M trees and speciation rate Mis (^(?;) ?a {aLlk)(l - M)* (eqn 1). In the Moran model, a = Mul(l - u). Hubbell assumed that the plot's most common tree species were equally dense in the source pool. In the log-series best fitting these common trees, a = 50. Setting a = 50 for the source pool, and setting immigration rate w = 0.1 fit the species abundance distribution of trees > 10 cm dbh on the plot. Volkov et al (2003) and McKane et al (2004) provide an easily understood way to calculate the species abun- dance distribution of trees on a plot, given m and the source pool's a. The number \?/(j) of species with j trees apiece on the plot is roughly M k=j where M is the number of trees in the source pool, 4>{k) is the number of species with k trees apiece in this pool, and Puii) is the probability that a species with k trees in the source pool has j of the plot's JV trees. If a species has k trees in the source pool, the probability that one of these trees is seed-parent to a tree on the plot at the next time-step is mklM. Thus, if this species now has no trees on the plot, the probability Pk{0,l) that it has one after the next time-step is mklM. If this species now has one tree on the plot, the probability PA:(1,0) that it dies at the next time-step and is replaced by some other species is (1/JV)[(1 - w)(JV - 1)/JV-i-OT(M -/i)/M)], where 1/JV is the probability that the one tree on the plot of this species dies, (1 - m){N - 1)/JV is the probability that it is replaced from the plot by another species, and m(M - k)lM is the probabifity that it is replaced from the source pool by another species. At steady state, P*(0)P*(0,1) = P*(1)P*(1,0), so Pfc(l) =P,(0)Pfc(0,l)/P,(l,0). Similarly, P^(l)P^{l,2) = P,{2)P,{2,1), so P,(2)=Pt(0)P,(0,l)Pi(l,2)/P,(l,0)Pi(2,l). The calculation continues similarly for larger j, and is completed by setting '^ZjLoPkij) = 1- The calculation is inexact, because it assumes that the memberships of species on the plot change independently, whereas they must actually sum to JV. Etienne (2005) provides more precise maximum-fikelihood methods of calculating m and the source pool's a and the most probable distribu- tion of species abundances on the plot. Is the migration rate m estimated by this method more than an arbitrary fitted parameter? If each tree's seeds are distributed about their parent according to the same radially symmetric Gaussian, then m = O.l implies a root mean square dispersal distance from parent to offspring ? 2007 THE AUTHOR. J. EVOL. BIOL. 20 (2007) 2075-2091 JOURNAL COMPILATION ? 2007 EUROPEAN SOCIETY FOR EVOLUTIONARY BIOLOGY The neutral theory: a historical perspective 2081 near 59 m (Leigh et al, 2004a). H. C. MuUer-Landau fit radially symmetric Gaussians to distributions of seeds about parent trees for 65 different species. Root mean square dispersal distance varied greatly from species to species, in violation of the neutral theory's assumptions, but the average over all species was 55 m (Condit et al, 2002). Is this match of observation to prediction mean- ingful? There is no such thing, however, as a panmictic source pool. To develop the neutral theory more reasonably, one must let seeds disperse only to a limited distance from their parents. Following Mal?cot (1948) and Nag- ylaki (1974), Chave & Leigh (2002) constructed a neutral theory of species turnover, assuming limited dispersal of seeds from parents. Consider a limitless forest with p trees per m^, where each tree has probability u of being a new species, a tree alive at time t has probability di of dying by time ? + di, and the probability is exp(-r^/o-^) that a dead tree is replaced by the immediately maturing young of a tree < r m away, where a is the root mean square distance from a tree to its young. Then, at steady state, the probability F(r) that two trees r km apart belong to the same species is nearly 2Ko[{r/a)^u\/[np(j^ +\n(\/u% if r > CT (Chave & Leigh, 2002). KQ is the modified Bessel function of order 0: when x > 3, KO(X)K?-^(nl2x)eTq){-x) (Olver, 1965). Here, x = rIL, where L = al^u is the distance a species is likely to spread before it dies out (Hubbell, 2001:174). Barro Colorado's plot has one tree > 10 cm dbh per 23.3 m^. If we set cr = 57 m, M = 4.8 X 10"** and 1/p = 23.3 m^, the predicted F(r) fits the trend of the data for trees > 10 cm dbh on 1-ha sample plots over 100 m apart in central Panama (Condit et al, 2002). To be sure, there was much scatter about this trend: species composition of plots on similar soils were more similar, and species composition of plots on very different soils were more different, than the neutral theory predicted (Condit et al, 2002). The root mean square dispersal distance of 57 m calculated by Condit et al. (2002) to fit the species turnover data is close both to 55 m, the average of the root mean square dispersal distances calculated for tree species on Barro Colorado's plot by H. C. Muller- Landau (Condit et al, 2002) and to the 59 m calculated by assuming that 10% of the trees on this plot have seed-parents outside it (Leigh et al, 2004a). Does this mean that it is safe, for certain predictions at least, to assume panmictic source pools? Evaluating Hubbell's neutral theory Hubbell's (1979, 1997, 2001) neutral theory of forest ecology has aroused abundant controversy, as one might expect, for it ignores Darwin's (1859:110-116) concept of adaptive divergence of character (Alonso et al, 2006). This concept was crucial to Darwin's (1859) explanation of diversity, as it was to Fisher's (1930:125-131), Lack's (1947), MacArthur's (1961, 1972), Mayr's (1963), Grant's (1986), Schluter's (2000), Mayr & Diamond's (2001) and many others. Neutral theorists' defence of a simple theory capable of useful quantitative prediction clashes with other biologists' pursuit of Darwin's approach. This section considers empirical objections to the neutral theory, and evaluates the theory's weak- nesses and strengths. Empirical objections to neutral theory Empirical objections to the neutral theory fall under four heads. First, its fundamental axiom, that all trees are alike before natural selection regardless of their species, is false (Hubbell, 2001:6). Different species respond differ- ently to environmental change (Condit et al, 1996b). Different species adjust differently to the master trade-off between growing fast in bright light and surviving in shade (Brokaw, 1987; Pacala et al, 1996; Wright et al, 2003). Different species are adapted to different habitats (Fine et al, 2004). Seedlings and saplings grow more slowly, and die more quickly, where conspecifics are more common (Harms et al, 2000; Hubbell et al, 2001). Most, if not all, of these differences reflect adaptive divergence. No one denies these facts, although some, such as Chave (2004:250-251), question their impact on the neutral theory's predictions. The second category of criticisms stems from tests of 'common-sense' deductions from the neutral theory. Terborgh ei a/, (1996:563) assumed that floodplain forests separated by upland forest bluffs developed 'in the absence of spatial correlation' with each other, so that, in a neutral world, their species compositions would develop independently. Mature floodplain forests 30 km apart on the Rio Manu, isolated from each other by 2 km of riverfront bluffs separating their floodplains, are unexpectedly similar in species composition and relative abundance of common species (Terborgh et al, 1996). Patterns of succession from newly formed riverbank are also remarkably similar on these floodplains (Foster et al, 1986; Salo et al, 1986; Terborgh & Petren, 1991). These similarities are especially remarkable because the species composition of floodplain forest differ greatly from that of the upland forest separating one floodplain from another. Hubbell (2001:331) countered that seed dispersal main- tained the homogeneity of the floodplain forests. None- theless, the similar patterns of successional development in these forests, and the striking differences between floodplain forest and adjoining upland forest, imply severe restrictions to the generality of the neutral theory. Clark & McLachlan (2003) assumed that the neutral theory implies that the species composition of forests 50 km apart in a uniform landscape drift independently, and should therefore diverge progressively. Examining pollen records for the last ten thousand years in lakes ? 2007 THE AUTHOR. J. EVOL. BIOL. 20 (2007) 2075-2091 JOURNAL COMPILATION ? 2007 EUROPEAN SOCIETY FOR EVOLUTIONARY BIOLOGY 2082 E. G. LEIGH JR scattered over 10 000 km , they found no evidence for progressive divergence in tree species composition, and declared the neutral theory falsified. Clark and McLachlan also cited Tsuga canadensis, which suddenly became rare over a million square kilometres over 5000 years ago, and recovered its former abundance after 1500 years, as if wasted by a new pathogen to which it later acquired immunity (Allison et ai, 1986). Hubbell & Foster (1986:322) had said that return to previous abundance after the cessation of an abundance-depress- ing disturbance is evidence against the neutral theory. Tsuga provides that evidence. In reply, Volkov et al. (2004), ignoring Tsuga, argued that dispersal can main- tain the homogeneity of the forest over some (unspec- ified) distance. Neither Hubbell nor Volkov used the neutral theory to predict the homogenizing power of dispersal: that requires reckoning with the limited dispersal of seeds from their parents. Instead, both relied on the fictitious concept of a panmictic source pool, where seeds can disperse from one end to the other. The third category of objections concerns neutral theorists' emphasis on species abundance distributions. Most neutral theory tests compare predicted species abundance distributions with those observed on 50-ha plots, especially Barro Colorado Island's plot. McGiU's (2003) claim that a log-normal fit Barro Colorado's species abundance distribution better than Hubbell's simulations led to progressively more exact calculations (Volkov etal, 2003; Etienne & Olff, 2004, 2005). The sampfing theory of Etienne & Olff (2005) calculated most likely values of source pool a and immigration rate m from the source pool from a plot's species abundance distribution. These methods, however, could not judge decisively between Etienne and Olff's neutral theory prediction and the log-normal. A simpler calculation predicts source pool a, given immigration rate m and the probability F that two trees > 10 cm dbh in a plot with AT trees are conspecific. If trees in the source pool are distributed over species in a log- series with a = ? (eqn 1), the probabifity that two trees sampled from the source pool are conspecific is 1/(1 -i- ?). If a fraction m of the plot's JV trees have seed-parents in the source pool, then, by eqn 3, \Nml{\ + e) + 1 - m)]/{Nm + 1 ? m On Barro Colorado's plot, ^7= 1/(38.0) (calculated from 1995 census data suppfied by Suzanne Lao from R. Condit and the Center for Tropical Forest Science). Here, if w = 0.1, then 6, the source pool's a, is 37.6; if m = 0.093, then 6 = 37.65. The most likely values of m and source pool cc calculated from this 1995 census by Etienne's sampling theory are 0.093 and 47.7 respec- tively (Chave et al., 2006). I do not know why these two methods give such different expectations for 6. The parameter m in Etienne's species abundance distribution is, moreover, seldom measured, nor have neutral theorists compared estimates made by different methods. A rough calculation (Leigh et al, 2004a:260- 261) suggests that iim = 0.093 for a 50-ha plot, as Chave et al. (2006) report for Barro Colorado, then root mean square dispersal distance cr is 55 m, close to the 57 m that gives the best fit for species turnover near Barro Colorado (Condit et al, 2002). The same rough calculation shows that if CT = 77.5 m, the value that best fits species turnover around Yasuni in Amazonian Ecuador (Condit et al, 2002), then, for Yasuni's 25-ha plot, m = 0.175. Applying Etienne's estimator to the trees > 10 cm dbh on the 25-ha plot at Yasuni yields a most fikely value of 0.43 for m (Chave etal, 2006). Is the fit for Barro Colorado fortuitous? As Ricklefs (2006:1426) remarks, proper test of a mechanistic theory involves more than curve-fitting. Furthermore, a log-series distribution in the source pool does not necessarOy imply neutral dynamics, as Gillespie (1991) pointed out for the neutral theory of population genetics. Gillespie's (1991) non-neutral explanation for the log-series has a parallel in ecology - tree species with identical mortality which coexist by temporally partition- ing tree fall gaps: in different years, different species are most successful in occupying these gaps (Chesson & Warner, 1981; Hatfleld& Chesson, 1989). Purves & Pa cala (2005) found that a log-series distribution in the source pool does not imply the absence of strong non-neutral processes. Regrettably, MacArthur's (1966) skepticism concerning the usefulness of species abundance distribu- tions has not been much heeded. The fourth group of criticisms centres on tests of neutral predictions of changes in time. For example, Leigh (1981) showed that, according to the neutral theory, the average time to extinction of a population now containing N adults was more than N generations. Without assembling supporting data, he claimed that very common species do not last so long in the fossil record, and concluded that extinctions of such species must be caused by some change in their environment. Ricklefs (2003), who assembled supporting data, inde- pendently found that common species are shorter lived than the neutral theory predicts. Neutral predictions of changes in time often make good nuO hypotheses because falsifying them, unlike falsifying most neutral predictions of species abundance distributions (McGOl et al, 2006), often suggests alternate hypotheses which shed light on the ecological organization of tropical forests. Two such tests will be discussed. The neutral theory predicts limits to how fast a tree population changes. If a tree species forms a smaU proportion of the forest, and if its trees die one by one, the number of its deaths should not differ more from the number of its recruitments than the number of heads from the number of tails in a corresponding number of tosses of a fair coin (cf. Gilbert et al, 2006:307). Census data from Barro Colorado's 50-ha plot was provided by Suzanne Lao on behalf of R. Condit and the Smithsonian's Center for Tropical Forest Science. Barro ? 2007 THE AUTHOR. J. EVOL. BIOL. 20 (2007) 2075-2091 JOURNAL COMPILATION ? 2007 EUROPEAN SOCIETY FOR EVOLUTIONARY BIOLOGY The neutral theory: a historical perspective 2083 Colorado's 1985 census counted 1900 Trkhilia tuberculata > 10 cm dbh, of which 239 died by the 1990 census, which counted 1783 T. tuberculata > 10 cm dbh. There- fore, 122 T. tuberculata recruited between 1985 and 1990 (1900 - 239 + 122 = 1783). If deaths and recruitments were equally likely, the expected number of deaths in this species between 1985 and 1990 would be (239 + 122)/2 = 180.5, and its standard deviation would be V(361/4) = 9.5. The number dying is over six standard deviations higher than predicted: the number of trees of this species declined significantly between 1985 and 1990. Similar reasoning shows that the numbers of trees > 10 cm dbh of three of the plot's eight most common species (omitting the clonal palm Oenocarpus mapora, where deaths are difficult to ascertain) increased signif- icantly during all four 5-year census intervals between 1985 and 2005, and two others decreased during every interval. For eight of the plot's 11 next most common species, the number of trees > 10 cm dbh dying between 1985 and 2005 was in significant imbalance with the number of recruitments. During any one census interval, deaths were in significant imbalance with recruitment in at least 11 of these 19 species (Table 1). Changes differed significantly among these tree species. Hubbell & Foster (1990) emphasized the prevalent disequOibrium among trees on this plot, without comparing the magnitude of the change in populations of different tree species with neutral predictions. The neutral theory also predicts how fast chance can spread a new tree species or clade, if it survives so long. Ocotea (Lauraceae) dispersed southward across the sea then separating North and South America, and Symphonia globulifera (Guttiferae) dispersed across the sea from Africa to South America, about 20 mifiion years ago (Chanderbali et al, 2001; Dick et al, 2003). 5. globulifera now averages more than one tree > 10 cm dbh per hectare throughout Amazonia: South America must contain over 10 mifiion reproductive adults of this species (Leigh et al, 2004b:449). If annual mortality averages 2%, less than 500 000 tree generations have elapsed since this species invaded South America. If so few invaders reached adulthood that only one tree's chloroplasts has surviving descendants, eqn 5 implies that the probability that Symphonia multiplied so quickly by chance is about e~^?, or 2x10"^ (Leigh et al, 2004b:448-449). The same applies to Ocotea. Other plants have also invaded Amazonia successfully. At least 20% of the 1000 free-standing woody species on a 25-ha plot in Amazonian Ecuador descend from invaders that crossed the sea after South America became an island continent (Pennington & Dick, 2004). Table 1 Balance between death and recruitment in selected tree species on the 50-ha forest dynamics plot of Barro Colorado Island, Panama. Species N85 D N90 D N95 D N 00 D N05 Trichina tuberculata 1900 t239 1783 t166 1681 t207 1572 t204 t1429 Faramea occidentaiis 1400 *146 1649 *196 1717 *202 1828 *260 ?1909 Pouisenia armata 857 t161 801 t134 755 t168 671 t128 t630 Aiseis blacldana 854 *19 936 *15 981 *20 1013 *35 *1046 Quararibea asterolepis 694 38 699 *30 724 37 723 47 714 Gustavia superba 641 11 649 8 644 22 638 t30 619 Virola sebifera 587 46 605 46 615 t70 589 t84 559 Hirtella triandra 552 *33 618 *22 681 *37 717 *39 *765 Ccrdia lasiocalyx 445 47 444 t92 364 t91 295 t65 t248 Guarea guidonia 362 *16 376 18 376 29 370 24 359 Socratea excrrhiza 357 t64 336 41 346 t56 325 t56 t297 Protium tenuifololium 350 39 358 *21 381 23 392 30 *406 Tetragastris panamensis 330 *12 362 *21 379 *15 397 21 *399 Prioria ccpaifera 308 *7 335 13 344 12 352 7 *357 Beilschmiedia p?ndula 304 19 296 15 295 t20 282 t22 t270 Tabernaemontana arb?rea 303 24 304 *14 323 *15 341 *20 *362 Cordia bicolor 259 *16 285 *23 326 29 328 31 323 IHeisteria concinna 255 *13 274 *7 288 16 288 22 *292 Cecropia insignis 246 40 257 32 263 *29 281 *41 *342 Numbers, ?V, of trees of each species in the census indicated, and the number D dying before the next census. Deaths are considered in significant imbalance with recruitment if, letting D be the number of deaths between censuses t and ? + 1, \N{t+ 1) - N{t)\ > 2^{N{t -HI) - N{t) + 2D}. Deaths are considered in significant imbalance with recruitment if \N{05) - N(85)| > 2V{JV(05) - ?V(85) + 2 ^ D}, where Y^D is the sum of the deaths in this species occurring during each of the four census intervals. * before the number of deaths denotes a significant increase within the census interval, and t denotes a significant decrease. * before N 05 denotes a significant increase over the whole 20-year period, and t before N 05 denotes a significant decrease during this period. ? 2007 THE AUTHOR. J. EVOL. BIOL. 20 (2007) 2075-2091 JOURNAL COMPILATION ? 2007 EUROPEAN SOCIETY FOR EVOLUTIONARY BIOLOGY 2084 E. G. LEIGH JR Although their advantages spread these species all across Amazonia, no species makes up more than a small proportion of the forest: the invaders could not replace their predecessors. Tropical tree species must therefore differ in ways that allow a great many of them to coexist, with none able to replace the others, contra Huston (1994). Understanding tropical tree diversity requires learning what these differences are. This finding, which is crucial for future research directions, was revealed by using a null hypothesis from the neutral theory. Weaknesses and strengths of neutral theory Now that the neutral theory has been evaluated in empirical and historical context, it is time to sum up this theory's current weaknesses and strengths. One weak- ness of the neutral theory is the narrow range of its concerns - an extraordinary happenstance, considering the theory's justly vaunted potential for synthesis. Unfike population geneticists, ecological neutral theorists have ignored changes in time, failing to test whether neutral patterns are generated by neutral processes. They have abandoned the task of testing whether many pairs of species diverged simultaneously in response to a new barrier to their neutral counterparts in population genetics (Hickerson et al, 2006). Most analytical neutral theorists have abjured spatially explicit theory. Chave & Leigh's (2002) spatially explicit theory of species turn- over assumes a most implausible steady state: unlike the population geneticist Nagylaki (1974, 1976), no ecologist has considered the manner, or the speed, of approach to this steady state. With one exception (Chave et al, 2002), analytical predictions of species-area curves, which require spatially explicit mathematics, have been left to professional mathematicians (Bramson et al., 1996, 1998), whose derivations few biologists can understand. Instead, ecological neutral theorists have focussed almost exclusively on species abundance distributions in plots that receive migrants from a panmictic source pool. They have fit predicted formulae to observed distributions of abundance (often very successfully) without trying to eliminate, or at least test the impor- tance of, unrealistic auxifiary assumptions such as pan- mictic source pools. Only professional mathematicians, not read by biologists, have tried to derive species abundance distributions where young disperse only a limited distance from their parents (Bramson et al., 1998). The probability u that a tree belongs to a new species, a fitted parameter, takes the place of a proper neutral theory of speciation. The values of u used to fit the 'species turnover' formula - the probability F(r) that two trees r km apart are conspecific - to different data sets are unsettfingly different. The value of u fitted for central Panama is 1300 times the value fitted for Yasuni in Ecuadorian Amazonia and 2 600 000 times that for Mantj in south-eastern Amazonian Peru (Condit et al.. 2002, legend to Fig. 2). The speciation rate u also appears in the fundamental biodiversity number ? = Nu, another fitted parameter. As N is the number of trees in a fictitious panmictic source pool, the presence of u in the parameter ? sheds no additional light on speciation rate. Much work must be done to construct a truly testable neutral theory, free of confounding auxifiary assump- tions (Dayton, 1973), that can generate trustworthy nuU hypotheses. These weaknesses may not ruin all usefulness of current neutral models. Latimer et al. (2005) used the neutral theory to compare dispersal limitation and speci- ation rate among woody plants of the Cape fynbos with those among Amazonian trees. Unfortunately, this anal- ysis treated a set of scattered plots as a single panmictic plot to calculate Nu for, and immigration rate m from, their source pool. Nonetheless, the authors' conclusion that speciation was far higher, and dispersal distances much shorter, in the fynbos than among Amazonian trees, is supported by empirical data. Surely, analysing these sets of plots by the methods of Condit et al. (2002) would have yielded the same conclusion with fewer questionable assumptions. Again, nothing appears more ridiculous than assum- ing that each tree has probability u of being the first member of a new species. Hubbell (2001) explored the effects of replacing mutational speciation by the random splitting of populations. Etienne et al. (2007) assumed that (mutational) speciation rate within a species was inversely proportional to its population size, so that each species produced daughter species at the same rate. The 'standard model' of HubbeO (2001) and Etienne (2005) outperformed both alternatives. Indeed, it appears that most plant species descend from small, local populations (Leigh et al, 2004b). If so, assuming that species start from single trees may not alter predictions of ecological interest. Equation 4 impfies that a new species with a single tree has probability I In of having descendants n generations later. \i k <^n, the probability that a new species starting with k trees has no descendants n generations later is (1 - 1/H)* Ki exp{-kln) ^ I - kin. This species has probability 1 - exp(-?:/n) K? kin of surviving to that time. The probability that only one of the k founders has descen- dants in the female fine then is {kln){\-kln) {k/n)exp{-k/n) The probability that, if this species survives to generation n, all members alive then descend in the female line from a single ancestor is [(fc/K)exp(-/i/M)]/[l - exp(-/i/?)].lf this is nearly 1, then the distribution of the number of descendants at generation n of a species starting with k trees, given that it still survives then, is the same as if the species had started with a single tree. Thus, if Nu adults belonging to new species are born into a source pool each generation, a plot's species abundance distribution should ? 2007 THE AUTHOR. J. EVOL. BIOL. 20 (2007) 2075-2091 JOURNAL COMPILATION ? 2007 EUROPEAN SOCIETY FOR EVOLUTIONARY BIOLOGY The neutral theory: a historical perspective 2085 be the same whether these Nu adults belong to one species, Nu species, or any number in between, so long as every species on the plot has at least n^ Nu adults in the source pool. Some weaknesses, however, will be harder to cure. Deriving spatially explicit theory of species-area curves or species abundance distributions for different scales, assuming that offspring disperse only a limited distance from their parents will not be easy. Despite Hubbell's breadth of vision, the neutral theory of forest ecology has yet to attract theoreticians with the nose for crucial problems and the mathematical skills to solve them of a Fisher, a Kimura, a Mal?cot or a Nagylaki. Nonetheless, if the interest and the promise of the neutral theory can be made clear without confusing it with truth, or even the truth of ideal-gas models, the theory should attract skilful workers. The neutral theory's fundamental, unavoidable weak- ness is, of course, its false basic assumption. Differences between tree species matter. Ignoring these differences makes the neutral theory unable to treat ecosystem properties (Westoby & Wright, 2006). Moreover, found- ing the neutral theory on the irrelevance of species differences means that we cannot know in advance which of its predictions will work. Why does the neutral theory predict (some!) species abundance dis- tributions so well, and species-area curves so poorly? Not knowing why this theory works where it does, we can never be sure to what new topics it can safely be applied. The neutral theory's false assumptions restrict its usefulness to: 1 domains where neutrality prevails, such as extinction from a plot of species represented there by one or two trees, or short-term fluctuations about equilibrium of a plot's tree species composition (Gilbert et al, 2006); 2 ?rst stabs at otherwise impenetrable topics, especially those requiring spatially explicit theory, such as spe- cies-area curves or species turnover; 3 predictions trancending the validity of the neutral theory's assumptions, that need more realistic deriva- tions, such as species abundance distributions; and 4 employment as a null hypothesis. The neutral theory wins adherents through its suc- cessful predictions, especially of species abundance dis- tributions and, to a lesser extent, species turnover. These successes, however, depend on unmeasured parameters such as the speciation rate u and unrealistic assumptions such as panmictic source pools. Nonetheless, the neutral theory may enable economical summary, and perhaps even a workable method of interpolation, for many types of diversity data from extensive, species-rich forests such as Amazonia. The neutral theory's fundamental and abiding strength, however, is the very simplicity of its assump- tions. These simple assumptions allow the development of theory on topics such as species-area curves, phylo- genetic patterns and species turnover, that are not otherwise amenable to precise, mechanistic theory. The neutral theory of species turnover has drawn attention to how the local origin and finite spread of different plant species contributes to species turnover (Leigh et al, 2004a). The neutral theory may well make other such contributions in future. At the very least, ecological neutral theory, like that in population genetics, provides natural, yet probing null hypotheses with which one can assess the significance of a variety of biological processes (Leigh et al, 1993, 2004b; Nee, 2005; Gilbert et al, 2006). Why was neutral theory more popular among ecologists? Finally, why the contrasting reception of neutral theory in population genetics vs. ecology? Population genetics was founded by evolutionary biologists (Fisher, 1918, 1930; Haldane, 1924, 1932; Wright, 1931, 1932). Darwin (1859) established the origin of adaptation as the central question of evolutionary biology, and population genet- icists helped show how effectively natural selection shaped adaptation (Fisher, 1930; Haldane, 1932; Wright, 1932). Early on, Wright (1931:153) thought that differ- ences between related species or genera were 'nonadap- tive'. Gause's (1934) competitive exclusion principle, and careful studies such as that of Pittendrigh (1950) on how related species coexist, convinced most evolutionary biologists that most species differences were adaptive (Hutchinson, 1959; Mayr, 1963). Although Kimura never denied the importance of natural selection in adaptive evolution, many disliked his idea that most gene substitutions were neutral. Surely, as knowledge in- creased, the idea of neutral genes would melt away (Gifiespie, 1991), as have the idea that species differences are nonadaptive, and the creationists' hope in unbridge- able missing links. Since 1930, Darwin's (1859, 1871) synthesis has acquired immense authority as the soundness of his judgment and the prescience of his vision have become progressively more evident. Nowa- days, most 'anti-adaptationists' in evolutionary biology merely seek to replace a nearly automatic presumption of adaptation by shared standards for testing adaptive explanations (Orzack & Sober, 2001). In population genetics, no neutral theorist seeks to overthrow or marginalize adaptation by natural selection: here, the scope of the neutral theory controversy is comparatively limited. Explaining genetic diversity triggered a bitter contro- versy that raged for many years. Was genetic diversity maintained by mutation-selection balance or by heter- osis (Lewontin, 1974)? Balance between mutation and genetic drift presented a third possibility. Genetic diver- sity, however, increased far more slowly with population size than the neutral theory predicted (Lewontin, 1974:209). A 'nearly neutral' theory, based on a prob- ability distribution of selective disadvantages of nearly ? 2007 THE AUTHOR. J. EVOL. BIOL. 20 (2007) 2075-2091 JOURNAL COMPILATION ? 2007 EUROPEAN SOCIETY FOR EVOLUTIONARY BIOLOGY 2086 E. G. LEIGH JR neutral al?eles, fit the relation between a population's size and its genetic diversity far better (Kimura, 1979), confirming the roles of both mutation-selection balance and neutral process in maintaining genetic diversity. The neutral theory reinforced the very fruitful (even if inexact) idea of clocked molecular phylogenies, contin- ues to inspire elegant, creative mathematics (Ewens, 2004), provides useful nufi hypotheses (Donnelly & lavar?, 1995; Kreitman & Akashi, 1995; Sabeti et al, 2002), suggests the surprising importance of genetic drift in evolutionary change and accounts for patterns of nucleotide substitutions at sOent sites (Ohta & Gillespie, 1996). Non-neutral population genetics, however, is stOl contributing significantly to evolutionary theory, not least by analyses of interactions between different levels of selection and other aspects of the evolution of cooperation (Maynard Smith & Szathm?ry, 1995; see also Leigh & Rowell, 1995). Accordingly, few evolution- ary biologists now view the neutral theory as the salvation of population genetics. Unlike evolutionary theory, ecology lacks a founder who established a magisterial synthesis. Darwin (1859:110-116) sketched the outlines of an evolutionary ecology, implicitly based on a principle of competitive exclusion and explicitly based on his principle of 'diver- gence of character'. He also used island biology to outfine an ecosystem ecology showing why larger land masses support more diversity and more intense competition, and are more resistant to introduced invaders (Darwin, 1859:104-108; Leigh et al, 2007). Darwin's ideas, how- ever, formed a few pages in his argument for the role of natural selection in evolution. Ecologists ignored them. Later, these ideas were rediscovered, in part, by Elton (1927), Lack (1947), Hutchinson (1959) and MacArthur (1961, 1972) among others. Many ecologists, however, denied the existence or the importance of density regulation and/or competitive exclusion (Andrewartha & Birch, 1954; Huston, 1994), even though these processes were central to Darwin's case for natural selection. The confusion was probably aggravated because theory played a less secure, less weU-defined role in ecology than in evolution, where population genetics was based on Mendel's laws. Ecologists divide sharply over what factors maintain species diversity, and why tropical settings are more diverse (Willig et al, 2003; Leigh ?tal, 2004b). Ecolo- gists disagree on whether trade-offs and avoiding competitive exclusion are the keys to understanding species diversity (Huston, 1994; Leigh et al, 2004b): no consensus is in sight. Some propose explicit mecha- nisms to explain diversity, such as the limitation of plant populations by specialized pests (Janzen, 1970), temporal segregation in recruitment (Chesson & Warner, 1981) and trade-offs such as growing fast in bright light vs. surviving in shade (Pacala et al, 1996), or growing fast vs. investing in anti-herbivore defence (Fine et al, 2006). Other explanations of diversity gradients, however, invoke factors such as mutation rate (Rohde, 1992) or aspects of climate (Hawkins et al, 2003); here, the distinction between cause and correlation, and between causal processes and the environmental conditions that modulate their effects, is less clear. The neutral theory held out the hope of replacing a chaos of competing qualitative explanations of diversity that no argument, observation or experi- ment seemed able to resolve, by a theory capable of successful quantitative prediction. No wonder that some greeted the neutral theory as ecology's salvation. There is a more fundamental issue. Everyone wants a quantitative theory that accounts for both pattern and the processes that generate it. No general theory in community ecology, however, has accomplished this. Neutral theorists have chosen to explain pattern, and ignore the processes that generate it: their critics prefer qualitative theory that deals with process. The contrast- ing attitudes to process of neutral theorists and their critics are iUustrated by a neutral theorist's reaction to Leigh et al's (1993) demonstration that a set of four tree species is spreading nonrandomly quickly on 80-Year-old islets without resident mammals in central Panama's Gatun Lake. Chave (2004:250) responded that 'only a few non-neutral species may have confused [this] anal- ysis'. Where Chave saw confusion, Leigh et al (1993) found a clue to why species diversity on these islets dropped nonrandomly quickly after they were isolated from the mainland. The four spreading species had large seeds, some of which escaped insect attack even if not buried, whereas seeds of most other large-seeded trees escape insect attack only if buried by agoutis. On these mammal-free islets, these four species have a great advantage over other large-seeded trees. The neutral theory could not explain the changes on these islets, but it provided the null hypothesis that helped to reveal what was happening there. More generally, understanding a community's response to any kind of change requires knowledge of how species differ, and how these differences affect their responses to environmental change (Paine, 1974, 1992, 2002; Foster, 1990; Leigh et al, 1993; King, 1994). Here, a qualitative theory of process is far more useful than a quantitative theory of pattern that cannot account for process. This circumstance does not lack irony. Some decades ago, the ecological theory of competition-driven adaptive divergence drew its evidence from pattern, not process (see Hutchinson, 1959; MacArthur, 1960, 1969, 1972). This theory became focused on process and based on experiment (see, for example, Kitajima, 1994; Schl?ter, 1994; Fine et al, 2004, 2006; Grant & Grant, 2006), largely thanks to prodding from critics such as Dayton (1973), Connor & Simberloff (1979) and Council (1980). Thanks to their critics, proponents of adaptive divergence can now argue the importance of process to pattern- oriented neutral theorists. ? 2007 THE AUTHOR. J. EVOL. BIOL. 20 (2007) 2075-2091 JOURNAL COMPILATION ? 2007 EUROPEAN SOCIETY FOR EVOLUTIONARY BIOLOGY The neutral theory: a historical perspective 2087 The usefulness of the useless A Chinese sage, Chuang Tzu, warned against ignoring the usefulness of the useless, remarking that 'only those who already know the value of the useless can be talked to about the useful' (Waley, 1939:3). His warning is especially apt for neutral theory. Neutral theory could rejuvenate forest ecology by uniting ecology and bioge- ography, especially if it includes evolution too. Its simple assumptions should enable neutral theory to base this synthesis on the limited dispersal of offspring from parents (Bramson et al, 1996, 1998; Chave & Leigh, 2002), and predict: 1 species abundance distributions on different-sized plots, species-area curves over different scales and species turnover (the decreased similarity of species composition between two plots with increased distance between them) where extinction balances speciation and dispersal; 2 the manner, and the speed, of approach to the above equilibrium; 3 the changes, short and long term, of species composi- tion on different-sized plots; 4 the distribution of different species' range sizes (Bell, 2001); 5 the rate of spread of new species lucky enough to survive; and 6 the average further lifetime of a species that now has N adults. This synthesis should be supplemented by a submodel embodying a neutral theory of speciation. Hubbell (2001) foresaw many elements of this synthesis, but much of what he foresaw has yet to be done. Such a synthesis would have many uses. First, it would serve as a model of what a more realistic theory should aim for. Second, there will be much to learn from why some neutral predictions are wrong. For example, understanding why the neutral theory of species-area curves fails for areas A in the range where the number S of species in area A varies as S = cA' might help us to decipher the biological basis of this empirical law. Third, this synthesis will provide null hypotheses by which the effectiveness of different biological processes can be assessed (Nee, 2005). There is no point in disproving the neutral theory for the pleasure thereof: everyone, even its advocates, knows that its assumptions are wrong (Hubbell, 2001:6). Falsifying the neutral theory is useful only when its predictions are compared with those of other mechanistic hypotheses (McGill et al, 2006:1421). A proper, spatially expUcit, neutral theory would allow us to detect whether a species, or a collection of species, is spreading, or a plot's species composition changing, nonrandomly quickly (cf. Leigh et al, 1993; Gilbert et al, 2006). It might also enable us to decide whether the plant species growing together on a small plot, or those that have successfully colonized an oceanic island, are more distantly, or more closely, related than we would expect by chance. Does a plot's environment favour a set of related species (Webb, 2000), or do differences in physiological function or the specialized pests attracted allow plants to grow and survive better if their nearest neighbours are more distantly related (Cavender-Bares et al, 2004; Webb et al, 2006)? Finally, neutral theory can serve as a stepping stone to a more realistic theory. Neutral theory may identify 'natural measures' of quantities, such as species diversity or similarity in species composition, which facilitate the development of more realistic theory (Leigh et al, 2004a). Neutral theory provides methods for tackling more realistic problems. The variety of 'diffusion meth- ods' in population genetics (Watterson, 1996) first used to study neutral processes (Fisher, 1922) has helped to solve a far wider range of problems (Kimura, 1994; Watterson, 1996; Ewens, 2004). Finafiy, Polya (1954) has emphasized the importance of solving related, simpler problems as stepping stones to solving more difficult ones. The neutral synthesis may often serve as this simpler stepping stone to developing a more realistic theory. For example, Zilfio et al (2005) modified neutral theory by assuming that each species was density regulated, as a first attempt at incorporating the effects of pest pressure. Neutral theory may serve as a particu- larly essential stepping stone to the spatially explicit theory ecology so sorely needs (Tilman & Kareiva, 1997). Often, the first neutral derivation will not permit gener- alization: only later ones will open practicable routes to a more realistic theory. The neutral theory will accordingly have little future unless some of its proponents revel in the beauty of its mathematics. An elegant, general, spatially explicit neutral theory, however, could provide an essential basis for both empirical and theoretical ecology. Acknowledgments I am most grateful to my Institute's librarians, Vielka Chang-Yau and Angel Aguirre, for procuring innumer- able 'pdfs,' to Suzanne Lao for providing data on Barro Colorado's 50-ha plot in usable form, and to Richard Condit, of the Smithsonian's Center for Tropical Forest Science, for permitting its use in this paper. I am also most grateful to J?r?me Chave and Rampal Etienne for their kindness in explaining various points of neutral theory to one of its more persistent critics. Finally, I am thankful to the plants and animals of Barro Colorado Island, and to Claudio Monteverdi's Selva Morale e Spirituale, for helping to set the controversies over the neutral theory in perspective. References Allison, T.D., Moeller, R.E. & Davis, M.B. 1986. Pollen in laminated sediments provides evidence for a mid-Holocene forest pathogen outbreak. Ecology 67: 1101-1105. ? 2007 THE AUTHOR. J. EVOL. BIOL. 20 (2007) 2075-2091 JOURNAL COMPILATION ? 2007 EUROPEAN SOCIETY FOR EVOLUTIONARY BIOLOGY 2088 E. G. LEIGH JR Alonso, D., Etienne, R.S. & McKane, A.J. 2006. The merits of neutral theory. TREE 8: 451-457. Andrewartha, H.G. & Birch, L.C. 1954. The Distribution and Abundance of Animals. University of Chicago Press, Chicago, IL. Bell, G. 2001. Neutral macroecology. Science 293: 2413-2418. Bernardi, G. 2007. The neoselectionist theory of genome evolution. PNAS 104: 8385-8390. Bramson, M., Cox, J.T. & Durrett, R. 1996. Spatial models for species area curves. Ann. Probab. 24: 1727-1751. Bramson, M., Cox, J.T. & Durrett, R. 1998. A spatial model for the abundance of species. Ann. Probab. 28: 658-709. Brokaw, N.V.L. 1987. Gap-phase regeneration of three pioneer tree species in a tropical forest. J. Ecol. 75: 9-19. Caswell, H. 1976. Community structure: a neutral model analysis. Ecol. Monogr. 46: 327-354. Cavender-Bares, J., Ackerly, D.D., Baum, D.A. & Bazzaz, F.A. 2004. Phylogenetic overdispersion in Floridian oak commu- nities. Am. Nat. 163: 823-843. Chanderbali, A.S., van der Werff, H. & Renner, S.S. 2001. Phylogeny and historical biogeography of Lauraceae: evidence from the chloroplast and nuclear genomes. Ann. MO Bot. Gard. 88: 104-134. Chave, J. 2004. Neutral theory and community ecology. Ecol. Lett. 7: 241-253. Chave, J. & Leigh, E.G. Jr 2002. A spatially explicit neutral model of /?-diversity in tropical forests. Theor. Popul. Biol. 62: 153-166. Chave, J., Muller-Landau, H.C. & Levin, S.A. 2002. Comparing classical community models: theoretical consequences for patterns of diversity. Am. Nat. 159: 1-23. Chave, J., Alonso, D. & Etienne, R.S. 2006. Comparing models of species abundance. Nature 441: (7089): El. doi: 10.1038/ nature04826. Chesson, P.L. & Warner, R.R. 1981. Environmental variability promotes coexistence in lottery competitive systems. Am. Nat. 117: 923-943. Clark, J.S. & McLachlan, J.S. 2003. Stability of forest biodiver- sity. Nature 423: 635-638. Condit, R., Hubbell, S.P., LaFrankie, J.V., Sukumar, R., Manok- aran, N., Foster, R.B. & Ashton, P.S. 1996a. Species-area and species-individual relationships for tropical trees: a compari- son of three 50-ha plots. J. Ecol. 84: 549-562. Condit, R., Hubbell, S.P. & Foster, R.B. 1996b. Changes in tree species abundance in a neotropical forest: impact of climate change. J. Trop. Ecol. 12: 231-256. Condit, R., Pitman, N., Leigh, E.G. Jr, Chave, J., Terborgh, J., Foster, R.B., et al. 2002. Beta-diversity in tropical forest trees. Science 295: 666-669. Condit, R., Leigh, E.G. Jr, Lao, S. & CTFS Working Group. 2004. Species-area relationships and diversity measures in the forest dynamics plots. In: Tropical Forest Diversity and Dynamism (E. C. Losos & E. G. Leigh Jr, eds), pp. 79-89. University of Chicago Press, Chicago, IL. Connell, J.H. 1980. Diversity and the coevolution of competi- tors, or the ghost of competition past. Oikos 35: 131-138. Connor, F.F. & Simberloff, D. 1979. The assembly of species communities: chance or competition? Ecology 60: 1132-1140. Crow, J.F. & Kimura, M. 1970. An Introduction to Population Genetics Theory. Harper and Row, New York. Darwin, C. 1859. On the Origin of Species. John Murray, London. Darwin, C. I87I. The Descent of Man and Selection in Relation to Sex. John Murray, London. Dayton, P.K. 1973. Two cases of resource partitioning in an intertidal community: making the right predictions for the wrong reason. Am. Nat. 107: 662-670. Dick, C.W., Abdul-Salim, K. & Bermingham, E. 2003. Molecular systematic analysis reveals cryptic Tertiary diversification of a widespread tropical rain forest tree. Am. Nat. 162: 691-703. Donnelly, P. 1986. Partition structures, Polya urns, the Ewens sampling formula, and the ages of al?eles. Theor. Popul. Biol. 30: 271-288. Donnelly, P. & Tavar?, S. 1987. The population genetics of the infinitely-manyneutralallelesmodel. J.A????/?.B/o/. 25: 381-391. Donnelly, P. & Tavar?, S. 1995. Coalescents and genealogical structure under neutrality. Annu. Rev. Genet. 29: 401-421. Elton, C. 1927. Animal Ecology. Sidgwick and Jackson, London. Etienne, R.S. 2005. A new sampling formula for neutral biodiversity. Ecol. Lett 8: 253-260. Etienne, R.S. & Alonso, D. 2005. A dispersal-limited sampling theory for species and al?eles. Ecol. Lett. 8: 1147-1156. Etienne, R.S. & Olff, H. 2004. A novel genealogical approach to neutral biodiversity theory. Ecol. Lett. 7: 170-175. Etienne, R.S. & Olff, H. 2005. Confronting different models of community structure to species-abundance data: a Bayesian model comparison. Ecol. Lett. 8: 493-504. Etienne, R.S., Apol, M.E.F., Olff, H. & Weissing, F.J. 2007. Modes of speciation and the neutral theory of biodiversity. Otkos 116: 241-258. Ewens, W.J. 1972. The sampling theory of selectively neutral al?eles. Theor. Popul. Biol. 3: 87-112. Ewens, W.J. 2004. Mathematical Population Genetics. L Theoretical Introduction. Springer, New York. Feller, W. 1968. An Introduction to Probability Theory and its Applications, Vol. I, 3rd edn. John Wiley, New York. Felsenstein, J. 2004. Inferring Phytogenies. Sinauer Associates, Sunderland, MA. Fine, P.V.A., Mesones, I. & Coley, P.D. 2004. Herbivores promote habitat specialization by trees in Amazonian forests. Science 305: 663-665. Fine, P.V.A., Miller, Z.J., Mesones, I., Irazuzta, S., Appel, H.M., Stevens, M.H.H., S??ksj?rvi, I., Schultz, J.C. & Coley, P.D. 2006. The growth-defense trade-off and habitat specialization by plants in Amazonian forests. Ecology 87 (Supplement): S150-SI62. Fisher, R.A. 1918. The correlation between relatives on the supposition of Mendelian inheritance. Trans. R. Soc. Edlnb. 52: 399-433. Fisher, R.A. 1922. On the dominance ratio. Proc. R. Soc. Edinb. 42: 321-341. Fisher, R.A. 1930. The Genetical Theory of Natural Selection. Clarendon Press, Oxford. Fisher, R.A. 1958. The Genetical Theory of Natural Selection, 2nd edn. Dover, New York. Fisher, R.A., Corbet, A.S. & Williams, C.B. 1943. The relation between the number of species and the number of individuals in a random sample of an animal population. J. Anim. Ecol. 12: 42-58. Foster, R.B. 1990. Long-term change in the successional forest community of the Rio Manu floodplain. In: Pour Neotropical Rainforests (A. H. Gentry, ed.), pp. 565-572. Yale University Press, New Haven, CT. ? 2007 THE AUTHOR. J. EVOL. BIOL. 20 (2007) 2075-2091 JOURNAL COMPILATION ? 2007 EUROPEAN SOCIETY FOR EVOLUTIONARY BIOLOGY The neutral theory: a historical perspective 2089 Foster, R.B., Arce, B.J. & Wacliter, T.S. 1986. Dispersal and the sequential plant communities in Amazonian Peru flood plain. In: Frugivores and Seed Dispersal (A. Estrada & T. H. Fleming, eds), pp. 357-370. W. Junk, Dordrecht. Cause, G.F. 1934. The Struggle for Existence. Wilhams & Wilkins, Baltimore, MD. Gilbert, B., Laurance, W.F., Leigh, E.G. Jr & Nascimento, H.E.M. 2006. Can neutral theory predict the responses of Amazonian tree communities to forest fragmentation? Am. Nat. 168: 304-317. Gillespie, J.H. 1986. Natural selection and the molecular clock. Mol. Biol. Evol. 3: 138-155. Gillespie, J.H. 1991. The Causes of Molecular Evolution. Oxford University Press, New York. Gillespie, J.H. 2004. Population Genetics, 2nd edn. Johns Hopkins University Press, Baltimore, MD. Grant, P.R. 1986. Ecology and Evolution of Darwin's Finches. Princeton University Press, Princeton, NJ. Grant, P.R. & Grant, B.R. 2006. Evolution of character displace- ment in Darwin's finches. Science 313: 224-226. Haldane, J.B.S. 1924. A mathematical theory of natural and arti- ficial selection. Part I. Trans. Camh. Philos. Soc. 23: 19-41. Haldane, J.B.S. 1932. The Causes of Evolution. Longmans Green, London. Haldane, J.B.S. 1957. The cost of natural selection. J. Genet. 55: 511-524. Harms, K.E., Wright, S.J., Calder?n, O., Hern?ndez, A. & Herr?, E.A. 2000. Pervasive density-dependent recruitment enhances seedling diversity in a tropical forest. Nature 404: 493-495. Hard, D.L. & Clark, A.G. 1997. Principles of Population Genetics, 3rd edn. Sinauer Associates, Sunderland, MA. Hatfield, J.S. &? Chesson, P.L. 1989. Diffusion analysis and spatial distribution of the two-species lottery competition model. Theor. Popul. Biol. 36: 251-266. Hawkins, B.A., Field, R., Cornell, H.V., Currie, D.J., Gu?gan, J.-F., Kaufman, D.M., Kerr, J.T., Mittelbach, G.G., Oberdorff, T., O'Brien, E.M., Porter, E.E. & Turner, J.R.G. 2003. Energy, water, and broad-scale geographic patterns of species richness. Ecology 8i: 3105-3117. Hickerson, M.J., Stahl, E.A. & Lessios, H.A. 2006. Test for simultaneous divergence using approximate Bayesian com- putation. Evolution 60: 2435-2453. Hoppe, F.M. 1987. The samphng theory of neutral al?eles and an urn model in population genetics. J. Math. Biol. 25: 123-159. Hubbell, S.P. 1979. Tree dispersion, abundance and diversity in a tropical dry forest. Science 203: 1299-1309. Hubbell, S.P. 1997. A unified theory of biogeography and relative species abundance and its application to tropical rain forests and coral reefs. Coral Reefs 16 (Suppl.): S9-S21. Hubbell, S.P. 2001. The Unified Neutral Theory of Biodiversity and Biogeography. Princeton University Press, Princeton, NJ. Hubbell, S.P. 2006. Neutral theory and the evolution of ecological equivalence. Ecology 87: 1387-1398. Hubbell, S.P. & Foster, R.B. 1986. Biology, chance, and history and the structure of tropical rain forest tree communities. In: Community Ecology (J. Diamond & T. J. Case, eds), pp. 314-329. Harper & Row, New York. Hubbell, S.P. & Foster, R.B. 1990. Structure, dynamics and equilibrium status of old-growth forest on Barro Colorado Island. In: Four Neotropical Rainforests (A. H. Gentry, ed.), pp. 522-541. Yale University Press, New Haven, CT. Hubbell, S.P., Ahumada, J.A., Condit, R. & Foster, R.B. 2001. Local neighborhood effects on long-term survival of individual trees in a neotropical forest. Ecol. Res. 16: 859-875. Huston, M.A. 1994. Biological Diversity. Cambridge University Press, Cambridge. Hutchinson, G.E. 1959. Homage to Santa Rosalia, or why are there so many kinds of animals? Am. Nat. 93: 145-159. Janzen, D.H. 1970. Herbivores and the number of tree species in tropical forests. Am. Nat 104: 501-528. Karlin, S. & McGregor, J. 1964. Direct product branching processes and related Markov chains. PNAS 51: 598-602. Karlin, S. & McGregor, J. 1967. The number of mutant forms maintained in a population. Proc. Fifth Berkeley Symp. Math. Stat. Probab. 4: 415-438. Karlin, S. & McGregor, J. 1972. Addendum to a paper by W. Ewens. Theor. Popul. Biol. 3: 113-116. Kimura, M. 1955. Solution of a process of random genetic drift with a continuous model. PNAS 41: 144-150. Kimura, M. 1968a. Evolutionary rate at the molecular level. Nature 217: 624-626. Kimura, M. 1968b. Genetic variability maintained in a finite population due to mutational production of neutral and nearly neutral isoalleles. Genet. Res. 11: 247-269. Kimura, M. 1979. Model of effectively neutral mutations where selective constraint is incorporated. PNAS 76: 3440-3444. Kimura, M. 1983. The Neutral Theory of Molecular Evolution. Cambridge University Press, Cambridge. Kimura, M. 1994. Population Genetics, Molecular Evolution, and the Neutral Theory: Selected Papers. University of Chicago Press, Chicago, IL. Kimura, M. & Crow, J.F. 1964. The number of al?eles that can be maintained in a finite population. Genetics 49: 725-738. Kimura, M. & Ohta, T. 1969. The average number of generations until fixation of a mutant gene in a finite population. Genetics 61: 763-771. Kimura, M. & Weiss, G.H. 1964. The stepping stone model of population structure and the decrease of genetic correlation with distance. Genetics 4:9: 561-576. King, D.A. 1994. Influence of light level on the growth and morphology of saplings in a Panamanian forest. Am. J. Bot. 81: 948-957. Kingman, J.F.C. 1977. The population structure associated with the Ewens sampling formula. Theor. Popul. Biol. 11: 274-283. Kingman, J.F.C. 1978. Random partitions in population genet- ics. Proc. R. Soc. Und. A 361: 1-20. Kingman, J.F.C. 1982. The coalescent. Stochastic Process. Appl. 13: 235-248. Kitajima, K. 1994. Relative importance of photosynthetic traits and allocation patterns as correlates of seedling shade toler- ance of 13 tropical trees. Oecologia 98: 419-428. Kreitman, M. & Akashi, H. 1995. Molecular evidence for natural selection. Annu. Rev. Ecol. Syst. 26: 403-422. Lack, D. 1947. Darwin's Finches. Cambridge University Press, Cambridge. Latimer, A.M., Silander, J.A. Jr & Cowling, R.M. 2005. Neutral ecological theory reveals isolation and rapid speciation in a biodiversity hot spot. Science 309: I722-I725. Leigh, E.G. Jr 1981. The average lifetime of a population in a varying environment. J. Theor. Biol. 90: 213-239. Leigh, E.G. Jr & Rowell, T.E. 1995. The evolution of mutuahsm and other forms of harmony at various levels of biological organization. ?cologie 26: I3I-I58. ? 2007 THE AUTHOR. J. EVOL. BIOL. 20 (2007) 2075-2091 JOURNAL COMPILATION ? 2007 EUROPEAN SOCIETY FOR EVOLUTIONARY BIOLOGY 2090 E. G. LEIGH JR Leigh, E.G. Jr, Wright, S.J., Putz, F.E. & Herre, E.A. 1993. The decline of tree diversity on newly isolated tropical islands: a test of a null hypothesis and some implications. Evol. Ecol. 7: 76-102. Leigh, E.G. Jr, Condit, R. & Loo de Lao, S. 2004a. The neutral theory of forest ecology. In: Tropical Forest Diversity and Dynamism (E. C. Losos & E. G. Leigh Jr, eds), pp. 244-263. University of Chicago Press, Chicago, IL. Leigh, E.G. Jr, Davidar, P., Dick, C.W., Puyravaud, J.-P., Terborgh, J., ter Steege, H. & Wright, S.J. 2004b. Why do some tropical forests have so many species of trees? Biotropica 36: 447-473. Leigh, E.G. Jr, Hladik, A., Hladik, CM. & Jolly, A. 2007. The biogeography of large islands, or, how does the size of the ecological theater affect the evolutionary play? Revue d'?cologie (La Terre et la Vie), in press. Lewontin, R.C. 1974. The Genetic Basis of Evolutionary Change. Columbia University Press, New York. MacArthur, R.H. I960. On the relative abundance of species. Am. Nat. 94: 25-36. MacArthur, R.H. 1961. Population effects of natural selection. Am. Nat. 95: 195-199. MacArthur, R.H. 1966. Note on Mrs. Pielou's comments. Ecology 47: 1074. MacArthur, R.H. 1969. Patterns of communities in the tropics. Biol. J. Linn. Soc. 1: 19-30. MacArthur, R.H. 1972. Geographical Ecology. Harper and Row, New York. MacArthur, R.H. & Wilson, E.O. 1967. The Theory of Island Biogeography. Princeton University Press, Princeton, NJ. Mal?cot, G. 1948. Les math?matiques de l'h?r?dit?. Masson et Cie, Paris. Mal?cot, G. 1969. The Mathematics of Heredity. W. H. Freeman, San Francisco, CA. Maynard Smith, J. & Szathm?ry, E. 1995. The Major Transitions in Evolution. Freeman/Spektrum, Oxford. Mayr, E. 1963. Animal Species and Evolution. Harvard University Press, Cambridge, MA. Mayr, E. & Diamond, J. 2001. The Birds of Northern Melanesia: Speciation, Ecology and Biogeography. Oxford University Press, New York. McGill, B.J. 2003. A test of the unified neutral theory of biodiversity. Nature 422: 881-885. McGill, B.J., Maurer, B.A. & Weiser, M.D. 2006. Empirical evaluation of neutral theory. Ecology 87: 1411-1423. McKane, A.J., Alonso, D. & Sol?, R.V. 2004. Analytic solution of Hubbell's model of local community dynamics. Theor. Popul. Biol. 65: 67-73. Moran, P.A.P. 1962. The Statistical Processes of Evolutionary Theory. Clarendon Press, Oxford. Nagylaki, T. 1974. The decay of genetic variability in geograph- ically structured populations. PNAS 71: 2932-2936. Nagylaki, T. 1976. The decay of genetic variability in geographically structured populations. II. Theor. Popul. Biol. 10: 70-82. Nee, S. 2005. The neutral theory of biodiversity: do the numbers add up? Fund. Ecol. 19: 173-176. Nei, M. 1987. Molecular Evolutionary Genetics. Columbia Univer- sity Press, New York. Nielsen, R. 2001. Statistical tests of neutrality in the age of genomics. Heredity 86: 641-647. Ohta, T. & Gillespie, J.H. 1996. Development of neutral and nearly neutral theories. Theor. Popul. Biol. 49: 128-142. Olver, F.W.J. 1965. Bessel functions of integer order. In: Handbook of Mathematical Functions (M. Abramowitz &? I. Stegun, eds), pp. 355-433. Dover, New York. Orzack, S.H. & Sober, E. (eds) 2001. Adaptationism and Optimality. Cambridge University Press, Cambridge. Pacala, S.W., Canham, CD., Saponara, J., Silander, J.A. Jr, Kobe, R.K. & Ribbens, E. 1996. Forest models defined by field measurements: estimation, error analysis and dynamics. Ecol. Monogr. 66: 1-43. Paine, R.T. 1974. Intertidal community structure. Oecologia 15: 93-120. Paine, R.T. 1992. Food-web analysis through field measure- ments of per capita interaction strength. Nature 355: 73-75. Paine, R.T. 2002. Trophic control of production in a rocky intertidal community. Science 296: 736-739. Pennington, R.T. & Dick, CW. 2004. The role of immigrants in the assembly of the South American rainforest tree fiora. Philos. Trans. R. Soc. Lond. B 359: I6II-I622. Pittendrigh, CS. 1950. The ecotypic specialization of Anopheles homunculus; and its relation to competition with A. bellator. Evolution 4: 64-78. Plotkin, J.B., Potts, M.D., Leslie, N., Manokaran, N., LaFrankie, J. & Ashton, P.S. 2000a. Species-area curves, spatial aggrega- tion, and habitat specialization in tropical forests. J. Theor. Biol. 207: 81-99. Plotkin, J.B., Potts, M.D., Yu, D.W., Bunyavejchewin, S., Condit, R., Foster, R., Hubbefi, S., LaFrankie, J., Manokaran, N., Seng, L.H., Sukumar, R., Novak, M.A. & Ashton, P.S. 2000b. Predicting species diversity in tropical forests. PNAS 97: I0850-I0854. Polya, G. 1954. Induction and Analogy in Mathematics. Princeton University Press, Princeton, NJ. Purves, D.W. & Pacala, S.W. 2005. Ecological drift in niche- structured communities: neutral pattern does not imply neutral process. In: Biotic Interactions in the Tropics (D. Burslem, M. Pinard & S. Hartley, eds) pp. 107-138. Cambridge Univer- sity Press, Cambridge. Ricklefs, R.E. 2003. A comment on Hubbell's zero-sum ecolog- ical drift model. Oikos 100: 185-192. Ricklefs, R.E. 2006. The unified neutral theory of biodiversity: do the numbers add up? Ecology 87: 1424-1431. Rohde, K. 1992. Latitudinal gradients in species diversity: the search for the primary cause. Oikos 65: 514-527. Sabeti, P.C., Reich, D.E., Higgins, J.M., Levine, H.Z.P., Richter, D.J., Schaffner, S.F., Gabriel, S.B., Platko, J.V., Patterson, N.J., McDonald, G.J., Ackerman, H.C, Campbell, S.J., Altshuler, D., Cooper, R., Kwiatkowski, D., Ward, R. & Lander, E.S. 2002. Detecting recent positive selection in the human genome from haplotype structure. Nature ^19: 832-837. Salo, J., Kalliola, R., H?kkinen, I., M?kinen, Y., Niemel?, P., Puhakka, M. & Coley, P.D. 1986. River dynamics and the diversity of Amazon lowland forest. Nature 322: 254?258. Sawyer, S. 1977a. Rates of consolidation in a selectively neutral migration model. Ann. Probab. 5: 486-493. Sawyer, S. 1977b. Asymptotic properties of the equilibrium probability of identity in a geographically structured popula- tion. Adv. Appl. Probab. 9: 268-282. Sawyer, S. 1979. A fimit theorem for patch sizes in a selectively- neutral migration model. J. Appl. Probab. 16: 482-495. ? 2007 THE AUTHOR. J. EVOL. BIOL. 20 (2007) 2075-2091 JOURNAL COMPILATION ? 2007 EUROPEAN SOCIETY FOR EVOLUTIONARY BIOLOGY The neutral theory: a historical perspective 2091 Sawyer, S. & Fleischmann, J. 1979. Maximum geographic range o? a mutant al?ele considered as a subtype of a Brownian branching random field. PNAS 76: 872-875. Schl?ter, D. 1994. Experimental evidence that competition promotes divergence in adaptive radiation. Science 266: 798-801. Schl?ter, D. 2000. The Ecology of Adaptive Radiation. Oxford University Press, New York. Terborgh, J. & Petren, K. 1991. Development of habitat structure through succession in an Amazonian floodplain forest. In: Habitat Structure: The Physical Arrangement of Objects in Space (S. S. Bell, E. D. McCoy & H. R. Mushinsky, eds), pp. 28-46. Chapman & Hall, New York. Terborgh, J., Foster, R.B. & Percy Nunez, V. 1996. Tropical tree communities: a test of the non-equilibrium hypothesis. Ecology 77: 561-567. Tilman, D. & Kareiva, P. 1997. Spatial Ecology. Princeton University Press, Princeton, NJ. Volkov, I., Banavar, J.R., Hubbell, S.P. & Maritan, A. 2003. Neutral theory and relative species abundance in ecology. Nature m-. 1035-1037. Volkov, I., Banavar, J.R., Maritan, S. & Hubbell, S.P. 2004. Neutral theory: the stability of forest biodiversity. Nature 427: 696. Waley, A. 1939. Three Ways of Thought in Ancient China. George Allen and Unwin, London. Watterson, G.A. 1974a. Models for the logarithmic species abundance distribution. Theor. Popul. Biol. 6: 217-250. Watterson, G.A. 1974b. The sampling theory of selectively neutral al?eles. Adv. Appl. Probab. 6: 463-488. Watterson, G.A. 1976. The stationary distribution of the infinitely-many neutral al?eles diffusion model. J. Appl. Probab. 13: 639-651. Watterson, G.A. 1984. Lines of descent and the coalescent. Theor. Popul. Biol. 26: 77-92. Watterson, G.A. 1996. Motoo Kimura's use of diffusion theory in population genetics. Theor. Popul. Biol. 49: 154-188. Webb, CO. 2000. Exploring the phylogenetic structure of ecological communities: an example for rainforest trees. Am. Nat. 156: 145-155. Webb, CO., Gilbert, G.S. & Donoghue, M.J. 2006. Phylodiver- sity-dependent seedling mortality, size structure and disease in a Bornean rain forest. Ecology 87 (Supplement): S123-S131. Westoby, M. & Wright, I.J. 2006. Land-plant ecology on the basis of functional traits. TREE 21: 261-268. Willig, M.R., Kaufman, D.M. & Stevens, R.D. 2003. Latitudinal gradients in biodiversity: pattern, process, scale and synthesis. Annu. Rev. Ecol. Evol. Syst. 34: 273-309. Wright, S. 1931. Evolution in Mendelian populations. Genetics 16: 97-159. Wright, S. 1932. The roles of mutation, inbreeding, crossbreed- ing and selection in evolution. Proc. Sixth Int. Congr. Genet. 1: 356-366. Wright, S. 1943. Isolation by distance. Genetics 28: 114-138. Wright, S. 1946. Isolation by distance under diverse systems of mating. Genetics 31: 39-59. Wright, S.J., Muller-Landau, H.C, Condit, R. & Hubbell, S.P. 2003. Gap-dependent recruitment, realized vital rates, and size distributions of tropical trees. Ecology 84: 3174-3185. Zilho, T., Volkov, I., Banavar, J.R., Hubbell, S.P. & Maritan, A. 2005. Spatial scaling in model plant communities. Phys. Rev. Lett. 95: 098I0I-I-098I0I-4. Zuckerkandl, E. & Pauling, L. 1965. Evolutionary divergence and convergence in proteins. In: Evolving Genes and Proteins (V. Bryson & H. J. Vogel, eds), pp. 97-116. Academic Press, New York. Received 6 February 2007; revised 14 June 2007; accepted 18 June 2007 ? 2007 THE AUTHOR. J. EVOL. BIOL. 20 (2007) 2075-2091 JOURNAL COMPILATION ? 2007 EUROPEAN SOCIETY FOR EVOLUTIONARY BIOLOGY