Evolution of the Human ASPM Gene, a Major Determinant of Brain Size

Corresponding author: Jianzhi Zhang, University of Michigan, 3003 Natural Science Bldg., 830 N. University Ave., Ann Arbor, MI 48109., jianzhi{at}umich.edu (E-mail)

Communicating editor: S. YOKOYAMA

	ABSTRACT

TOP ABSTRACT MATERIALS AND METHODS RESULTS DISCUSSION LITERATURE CITED

The size of human brain tripled over a period of 2 million years (MY) that ended 0.2–0.4 MY ago. This evolutionary expansion is believed to be important to the emergence of human language and other high-order cognitive functions, yet its genetic basis remains unknown. An evolutionary analysis of genes controlling brain development may shed light on it. ASPM (abnormal spindle-like microcephaly associated) is one of such genes, as nonsense mutations lead to primary microcephaly, a human disease characterized by a 70% reduction in brain size. Here I provide evidence suggesting that human ASPM went through an episode of accelerated sequence evolution by positive Darwinian selection after the split of humans and chimpanzees but before the separation of modern non-Africans from Africans. Because positive selection acts on a gene only when the gene function is altered and the organismal fitness is increased, my results suggest that adaptive functional modifications occurred in human ASPM and that it may be a major genetic component underlying the evolution of the human brain.

AMONG mammals, humans have an exceptionally big brain relative to their body size. For example, in comparison with chimpanzees, the brain weight of humans is 250% greater while the body is only 20% heavier (MCHENRY 1994 ). The dramatic evolutionary expansion of the human brain started from an average brain weight of 400–450 g 2–2.5 million years (MY) ago and ended with a weight of 1350–1450 g 0.2–0.4 MY ago (MCHENRY 1994 ; WOOD and COLLARD 1999 ). This process represents one of the most rapid morphological changes in evolution. It is generally believed that the brain expansion set the stage for the emergence of human language and other high-order cognitive functions and that it was caused by adaptive selection (DECAN 1992 ), yet the genetic basis of the expansion remains elusive. A study of human mutations that result in unusually small brains may help identify the genetic modifications that contributed to the human brain expansion. In this regard, primary microcephaly (small head) is of particular interest (MOCHIDA and WALSH 2001 ; BOND et al. 2002 ; KUMAR et al. 2002 ). Microcephaly is an autosomal recessive genetic disease with an incidence of 4–40 per million live births in western countries (MOCHIDA and WALSH 2001 ; KUMAR et al. 2002 ). It is defined as a head circumference >3 standard deviations below the population age-related mean, but with no associated malfunctions other than mild-to-moderate mental retardation (MOCHIDA and WALSH 2001 ; KUMAR et al. 2002 ). The reduction in head circumference correlates with a markedly reduced brain size. Microcephaly is genetically heterogeneous, associated with mutations in at least five loci (MOCHIDA and WALSH 2001 ; KUMAR et al. 2002 ), one of which was recently identified and named ASPM (abnormal spindle-like microcephaly associated; BOND et al. 2002 ). Four different homozygous mutations in ASPM introducing premature stop codons were found to cosegregate with the disease in four respective families, while none of these mutations were found in 200 normal human chromosomes (BOND et al. 2002 ). Because the brain size of a typical microcephaly patient (430 g; MOCHIDA and WALSH 2001 ; KUMAR et al. 2002 ) is comparable with those of early hominids such as the 2.3- to 3.0-MY-old Australopithecus africanus (420 g; MCHENRY 1994 ; WOOD and COLLARD 1999 ), I hypothesize that ASPM may be one of the genetic components underlying the human brain expansion. Signatures of accelerated evolution of ASPM under positive selection during human origins would strongly support my hypothesis, because the action of positive selection indicates a modification in gene function resulting in elevated organismal fitness (ZHANG et al. 2002 ). Below I provide population genetic and molecular evolutionary evidence for the operation of such adaptive selection on ASPM.

	MATERIALS AND METHODS

TOP ABSTRACT MATERIALS AND METHODS RESULTS DISCUSSION LITERATURE CITED

Sequencing of ASPM:
The human ASPM gene has 28 exons. All 28 exons were PCR amplified from genomic DNA samples of 14 human (Homo sapiens) individuals of different geographic origins (2 African Pygmies, 3 African Americans, 4 Europeans, 2 Southeast Asians, 1 Chinese, 1 Pacific islander, and 1 South American), using the high-fidelity Taq of Invitrogen (Carlsbad, CA). The PCR products were then purified and sequenced in both directions. Polymorphisms that appear only once (singletons) were confirmed by a second PCR-sequencing experiment. The human DNA samples were purchased from Coriell (Camden, NJ). I also amplified all 28 exons from one chimpanzee (Pan troglodytes) and one orangutan (Pongo pygmaeus) and sequenced the insert DNAs after the PCR products were cloned into the pCR4TOPO vector (Invitrogen). To trace the evolutionary origin of a large insertion/deletion in exon 18, two segments (Fig 1I and II) of exon 18 were also amplified and directly sequenced from genomic DNAs of hyrax (Procavia capensis), sea lion (Zalophus californianus), seal (Phoca vitulina), wolverine (Gulo gulo), fox (Alopex lagopus), dog (Canis familiaris), bear (Ursus maritimus), cat (Felis catus), pig (Sus scrofa), cow (Bos taurus), whale (Balaena mysticetus), rhesus monkey (Macaca mulatta), owl monkey (Aotus trivirgatus), and hamster (Cricetulus griseus).

View larger version (7K): [in this window] [in a new window] [Download PPT slide]   Figure 1. Nucleotide polymorphisms and substitutions in the human ASPM gene. The exons are shown by solid boxes and introns by lines. The sizes of exons, but not introns, are drawn to scale. The thick horizontal line above exon 18 shows the 867-nucleotide fragment that is not found in the mouse Aspm, with the two thin lines above it showing the two segments amplified in various mammals (see text and Fig 3). Fixed nucleotide substitutions in the human lineage after the human-chimpanzee split are shown by circles, whereas common and rare polymorphisms within humans are shown by squares and triangles, respectively. Nonsynonymous changes are solid symbols and synonymous changes are open symbols.

Data analysis:
The d_N/d_S ratios (LI 1997 ; NEI and KUMAR 2000 ) in the human, chimpanzee, and orangutan branches of the ASPM gene tree (Fig 2) were estimated using the maximum-likelihood method (YANG 1997 ), under the model of unequal codon frequencies (CodonFreq = 3) and unequal rates of transitions and transversions. This model fits the data significantly better than those with fixed codon frequencies (CodonFreq = 2) or equal rates of transitions and transversions. The hypothesis of d_N/d_S = 1 for a given branch and the hypothesis of equal d_N/d_S ratios between two branches were tested using the likelihood-ratio test, as well as a nonlikelihood method based on inferred ancestral sequences (ZHANG et al. 1997 , ZHANG et al. 1998 ). The ASPM sequence for the common ancestor of humans and chimpanzees was estimated using the Bayesian method (YANG 1997 ; ZHANG and NEI 1997 ). Population genetic analysis was conducted using DnaSP3.99 (ROZAS and ROZAS 1999 ). TAJIMA's (1989), FU and LI's (1993), and FAY and WU's (2000) tests were conducted using coalescent simulations under the assumption of no recombination across the gene. This assumption makes the tests more conservative. When recombination is considered, the P values become lower in Tajima's and Fu and Li's tests, but remain virtually the same in Fay and Wu's test. The recombination rate in the ASPM locus was estimated from KONG et al. 2002 to be 1.8 cM/10⁶ nucleotides. Since the coding sequences analyzed here span 62 x 10³ nucleotides of the genome, the recombination rate for the sequence is r = 1.8 x 10^-2 x 0.062 = 1.1 x 10^-3 recombination/generation. Assuming that the human effective population size is 10⁴ (TAKAHATA et al. 1995 ; HARPENDING et al. 1998 ), the population recombination parameter (4Nr) equals 4 x 10⁴ x 1.1 x 10^-3 = 44. Gene trees of segment I and II sequences from various mammals were reconstructed by the neighbor-joining method (SAITOU and NEI 1987 ) with KIMURA's (1980) two-parameter distances. Use of other distances or other methods (NEI and KUMAR 2000 ) does not change the result. Two thousand bootstrap replications (FELSENSTEIN 1985 ) were used to test the reliability of the trees. The MEGA2 software (KUMAR et al. 2001 ) was used for the phylogenetic analysis.

View larger version (12K): [in this window] [in a new window] [Download PPT slide]   Figure 2. The dN/dS ratios for the entire ASPM coding sequence in the human, chimpanzee, and orangutan lineages. The values were estimated from a likelihood analysis.

Computer simulation:
Computer simulation of ASPM evolution under pure neutrality was conducted following the procedure described in ZHANG and WEBB 2003 . The rate that an open reading frame (ORF) becomes disrupted is mainly determined by the sequence of the ORF, the rate of point mutations, and the rate of insertion/deletion (indel) mutations. Here I used a point mutation rate of 1 x 10^-9/site/year, which was estimated from a genomic comparison between humans and chimpanzees (YI et al. 2002 ). The relative mutational frequencies among the four nucleotides have only a negligible effect on the simulation result and I assumed that they are equal. I used an indel mutation rate of 1.1 x 10^-10/site/year, which was also estimated from a genomic comparison between the human and chimpanzee (BRITTEN 2002 ). I assumed that all indels with sizes that are multiples of three nucleotides do not disrupt an ORF. This simplifies the simulation but does not affect the results, because the majority of indels generated by mutations have small sizes (less than or equal to nine nucleotides; BRITTEN 2002 ). In the above genomic data analyzed, 19% of the total 1019 indels are of sizes that are multiples of three nucleotides. A simulation was then performed for 20,000 replications with the inferred ASPM coding sequence of the common ancestor of humans and chimpanzees. Under no functional constraints, the substitution rate is identical to the mutation rate and mutations are assumed to be random. An ORF is interrupted when an indel of a size that is not a multiple of three nucleotides or a nonsense point mutation occurs. I thus determined the t_1/2 for ASPM, or the time required for an intact ORF to be interrupted in half of the simulation replications. Under the above parameters, t_1/2 = 0.48 MY. When point and indel mutation rates are both halved, t_1/2 = 0.97 MY. The probability that ASPM retains its ORF after T MY of neutral evolution is ()^T/t_1/2. This estimation of the rate of pseudogenization is conservative, because other deleterious events such as insertions of transposable elements and null mutations at promoter regions and splicing sites are not considered here. The computer program for the simulation is from ZHANG and WEBB 2003 .

	RESULTS

TOP ABSTRACT MATERIALS AND METHODS RESULTS DISCUSSION LITERATURE CITED

Elevation of d_N/d_S in the human ASPM lineage:
Human ASPM has 28 coding exons, spanning 62 kb in chromosome 1p31 and encoding a huge protein of 3477 amino acids (Fig 1). I determined the entire coding sequences of ASPM from one human, one chimpanzee, and one orangutan, and compared them in the phylogenetic tree of the three species (Fig 2). The orangutan sequence is used as the outgroup for humans and chimpanzees so that nucleotide substitutions on the human and chimpanzee lineages can be separated. I did not sequence the gorilla because the gorilla sequence may not be appropriate as the outgroup due to incomplete lineage sorting (SATTA et al. 2000 ). Use of orangutan, a slightly more distant outgroup, solves this problem. A commonly used indicator of natural selection at the DNA sequence level is the ratio of the rate of nonsynonymous nucleotide substitution (d_N) to that of synonymous substitution (d_S). Most functional genes show d_N/d_S < 1, because a substantial proportion of nonsynonymous mutations are deleterious and are removed by purifying selection, whereas synonymous mutations are more or less neutral and are generally uninfluenced by selection. A gene may occasionally exhibit d_N/d_S > 1 when a large fraction of nonsynonymous mutations are advantageous and are driven to fixation by positive selection (LI 1997 ; NEI and KUMAR 2000 ). I estimated the d_N/d_S ratio for ASPM in each of the three tree branches (Fig 2), using a maximum-likelihood method, and found that d_N/d_S is lowest in the orangutan branch (0.43), higher in the chimpanzee branch (0.66), and highest in the human branch (1.03). The hypothesis of d_N/d_S = 1 is rejected for the orangutan branch (P < 0.001, likelihood-ratio test), but not for the other two branches, suggesting a difference in selection has occurred. Indeed, a test of the difference in d_N/d_S between the human and orangutan branches gives a marginally significant result (P = 0.064), but the difference between the chimpanzee and orangutan branches is not significant (P = 0.29), nor is the difference between human and chimpanzee branches (P = 0.45). Because the d_N/d_S ratio between the orangutan and mouse (Mus musculus) is also low (0.29), an increase of d_N/d_S in humans is more likely than a decrease in orangutans. The mouse sequence (GenBank accession no. AF533752) was not included in the phylogeny-based analysis as it is relatively distantly related to the ape sequences and contains multiple insertions and deletions, which would make the inference less reliable. Similar results are obtained when I first infer the ASPM sequence for the common ancestor of humans and chimpanzees and then estimate the d_N/d_S ratio by counting the numbers of synonymous and nonsynonymous nucleotide substitutions on each branch. For instance, this approach gives d_N/d_S = 1.13, 0.84, and 0.52, respectively, for the human, chimpanzee, and orangutan branches.

Complete functional relaxation does not adequately explain the elevation of d_N/d_S:
Two hypotheses may explain the increase in d_N/d_S to 1.03 during the evolution of human ASPM. First, the functional constraints and purifying selection on ASPM may have been completely relaxed and many deleterious nonsynonymous mutations were fixed by random genetic drift. Alternatively, advantageous nonsynonymous substitutions under positive selection occurred at some sites, while purifying selection acted at some other sites, resulting in an average d_N/d_S of 1. Under the first hypothesis, ASPM has been under pure neutral evolution since the human-chimpanzee separation 6–7 MY ago (BRUNET et al. 2002 ). Using rates of single-nucleotide mutations and insertion/deletion mutations estimated from human-chimpanzee genomic comparisons (BRITTEN 2002 ; YI et al. 2002 ), I conducted a computer simulation of neutral evolution of ASPM (see MATERIALS AND METHODS). I found that the probability that ASPM retains its open reading frame after 6 MY of neutral evolution is extremely low (1.7 x 10^-4). Even when the above two mutation rates are both halved, the probability is still very small (0.014), suggesting that ASPM must have been under purifying selection. The fact that nonsense mutations in ASPM lead to microcephaly also demonstrates the presence of functional constraints on the gene. Thus, the hypothesis of complete relaxation of functional constraints and lack of purifying selection for the past 6–7 MY of human evolution is inconsistent with the data, and some sites in ASPM must have been subject to purifying selection (d_N/d_S < 1). This result would imply, although not prove, that some other sites are under positive selection (d_N/d_S > 1), so that the average d_N/d_S across the entire protein is 1. However, it is difficult to rule out the possibility of an incomplete functional relaxation in human ASPM, which can lead to a d_N/d_S ratio of 1 when the number of substitutions is relatively small. A population genetic study may help resolve this question.

Signatures of purifying selection from population genetic data:
The entire coding sequence of ASPM is determined from 14 human individuals of different geographic origins. A total of 33 single-nucleotide polymorphisms are found (Table 1 and Table 2). The derived and ancestral alleles are inferred using the chimpanzee and orangutan sequences as outgroups. TAJIMA's (1989) and FU and LI's (1993) tests reveal slight departure of the data from the Wright-Fisher model of neutrality (D = -1.29, P = 0.081; F = -1.76, P = 0.074; Table 2). But FAY and WU's (2000) test, which is designed to detect recent selective sweeps, does not show a significant result (H = -2.08, P = 0.21). Thus, the negative D and F likely reflect recent population expansions and/or purifying background selection. A recent study suggested that negative D values may also be found under certain sampling schemes if there is fine-scale population differentiation (PTAK and PRZEWORSKI 2002 ). When the synonymous and nonsynonymous sites were analyzed separately, I detected significant negative D and F values at nonsynonymous sites (P < 0.05; Table 2), but not at synonymous sites. H is not significant at either type of site. These results suggest that the nonsynonymous sites in human ASPM are subject to purifying selection. It should be mentioned that the recombination rate in the ASPM region is 1.8 cM/10⁶ nucleotides (KONG et al. 2002 ), which translates into 1.1 x 10^-3 recombination/meiosis for the sequences analyzed here. This relatively high recombination rate localizes signatures of selection to a small region surrounding the selected sites. This might in part explain the above differences in the test results between synonymous and nonsynonymous sites.

Population genetic theory predicts that deleterious mutations do not reach high frequencies in populations, while neutral and advantageous mutations do. A comparison between rare and common polymorphisms may detect purifying selection of deleterious mutations (FAY et al. 2001 ). Fay et al. recommended a frequency of 10% for the derived allele as a cutoff between rare and common polymorphisms (FAY et al. 2001 , FAY et al. 2002 ). In the present sample of 28 chromosomes, derived alleles that appear one or two times are regarded as rare polymorphisms, and the rest are common. Because of the limited sample size, a truly rare allele may inadvertently appear more than twice in our sample and a truly common allele may inadvertently be regarded as rare. Using probability theory, I computed that the probability of the former error is <5% for an allele with frequency <3% and the probability of the latter error is <5% for an allele with frequency >20%. Thus, the present classification of rare and common alleles is expected to be relatively accurate. I observed that n_R = 15 nonsynonymous and s_R = 5 synonymous rare polymorphisms and n_C = 5 nonsynonymous and s_C = 8 synonymous common polymorphisms from the present data (Table 2; Fig 1). The ratio of n_C to n_R (5/15 = 0.333) is significantly lower than that of s_C to s_R (8/5 = 1.6; ² = 4.41, P < 0.05; Table 2). Since synonymous mutations are more or less neutral, the observed deficit of common nonsynonymous polymorphisms suggests that purifying selection has prevented the spread of nonsynonymous deleterious mutations. It is estimated by the likelihood method that there are N = 7459 and S = 2972 potentially nonsynonymous and synonymous sites in ASPM, respectively. Thus, for rare polymorphisms, there are n_R/N = 15/7459 = 2.01 x 10^-3 polymorphisms/nonsynonymous site and s_R/S = 5/2972 = 1.68 x 10^-3/synonymous site. Their difference is statistically insignificant (² = 0.09, P > 0.5). In contrast, for common polymorphisms, the number is significantly smaller per nonsynonymous site (n_C/N = 5/7459 = 0.67 x 10^-3) than per synonymous site (s_C/S = 8/2972 = 2.69 x 10^-3; ² = 6.98, P < 0.01), confirming that purifying selection has reduced the number of common nonsynonymous polymorphisms. This result also suggests the absence or rareness of advantageous nonsynonymous polymorphisms of ASPM that are currently segregating in humans, as such polymorphisms would predominantly show up as common polymorphisms and render n_C/N higher. This is consistent with the above result from Fay and Wu's test. The proportion of nonsynonymous polymorphisms not under purifying selection may be estimated by (n_C/N)/(s_C/S) = (0.67 x 10^-3)/(2.69 x 10^-3) = 0.25 or by (n_C/s_C)/(n_R/s_R) = (5/8)/(15/5) = 0.21. The two estimates are close to each other and to the d_N/d_S ratio between the mouse and orangutan (0.29). This indicates that human ASPM is currently under relatively strong purifying selection, and the strength of selection is comparable to or even greater than that in the long-term evolution of mammalian ASPM.

Comparison of polymorphism and divergence suggests past positive selection:
Because both the synonymous and nonsynonymous common polymorphisms are largely neutral, comparing them with the fixed substitutions between humans and chimpanzees can reveal the signature of selection that has influenced the substitution processes (MCDONALD and KREITMAN 1991 ; FAY et al. 2001 , FAY et al. 2002 ; SMITH and EYRE-WALKER 2002 ). This comparison shows a significant excess of fixed nonsynonymous substitutions (² = 3.88, P < 0.05, Table 2), suggesting that some nonsynonymous substitutions were fixed by positive selection. Because the expansion of brain size occurred in the human lineage after the human-chimpanzee split, it is more relevant to examine whether the human branch exhibits an excess of nonsynonymous substitutions. For this, the ASPM sequence of the common ancestor of humans and chimpanzees was inferred by the Bayesian method. Because the sequences considered are closely related, this inference is reliable, with the average posterior probability >0.999. Comparing the ancestral sequence with the polymorphic human sequences, I identified 16 nonsynonymous and 6 synonymous mutations that have been fixed in the human lineage (Table 2; Fig 1). Their ratio (16/6 = 2.67) is significantly greater than that for common polymorphisms (n_C/s_C = 5/8 = 0.63; ² = 4.00, P < 0.05). The number of neutral nonsynonymous substitutions may be estimated from the number of synonymous substitutions multiplied by n_C/s_C, which yielded 6 x (5/8) = 3.75 (FAY et al. 2001 , FAY et al. 2002 ; SMITH and EYRE-WALKER 2002 ). The number of nonsynonymous substitutions unexplainable by neutral evolution is 16 - 3.75 = 12, which may have been fixed by positive selection. It should be noted that a recent population expansion can cause an overestimate of the number of adaptive substitutions when slightly deleterious mutations are present. However, such overestimation is unlikely in the present case because the current effective population size of humans, even after the recent expansion, is still smaller than the long-term effective population size separating humans and chimpanzees and the effective population size of the common ancestor of humans and chimpanzees (TAKAHATA et al. 1995 ; CHEN and LI 2001 ; KAESSMANN et al. 2001 ; EYRE-WALKER 2002 ). It is interesting that there is no significant excess of nonsynonymous substitutions for either the chimpanzee or orangutan branches when the common polymorphisms and substitutions are compared (P > 0.05).

IQ repeats and brain size variation:
Human ASPM contains multiple calmodulin-binding IQ repeats (BOND et al. 2002 ). In a comparison of putative orthologous ASPM genes from the human, mouse, fruit fly (Drosophila melanogaster), and nematode (Caenorhabditis elegans), BOND et al. 2002 noticed that organisms with larger brains have more IQ repeats, implying a possible relation of IQ repeats and brain size. In particular, the predominant difference between the human and mouse ASPM genes is a large IQ-repeat-encoding insertion of 867 nucleotides at the end of exon 18. However, my data showed no difference in the number of IQ repeats between human and chimpanzee ASPM sequences. To trace the origin of the large insertion in human ASPM, I amplified and sequenced from several mammals two DNA segments that cover most of the insertion (Fig 1). Segment I is of 212 nucleotides and segment II is of 706 nucleotides. One or both segments were obtained from species belonging to primates, Cetartiodactyla, Carnivora, and Hyracoidea, but not from mouse or hamster (Fig 3). Phylogenetic analyses were conducted to confirm that the obtained sequences are orthologous to the human sequence (Fig 4). While nonamplification of a sequence does not prove its nonexistence, the amplification of the orthologous sequence indicates its presence. From the recently established mammalian phylogeny (MURPHY et al. 2001 ), it can be inferred that the large human insertion was already present in the common ancestor of most placental mammals, but was deleted in mouse and possibly in other rodents (Fig 3). Thus, this IQ-repeat-containing sequence does not explain the brain size variation among many nonrodent mammals.

View larger version (27K): [in this window] [in a new window] [Download PPT slide]   Figure 3. Amplification of segments I and II of the 867-nucleotide IQ-repeat-containing fragment in exon 18 of orthologous ASPM genes of mammals. The two segments are shown in Fig 1. "+" represents a successful amplification and sequencing and "-" represents no amplification. Note that the entire coding sequence of the gene is known for human, chimpanzee, orangutan, and mouse. The phylogenetic relationships of the mammalian orders are depicted following MURPHY et al. 2001 .

View larger version (20K): [in this window] [in a new window] [Download PPT slide]   Figure 4. Neighbor-joining trees of (A) segment I and (B) segment II sequences from ASPM exon 18 (see Fig 1). The aligned sequences of segment I have 155 nucleotide sites after the removal of primer-encoded regions, and those of segment II have 651 sites. Kimura's two-parameter distances are used. Bootstrap percentages from 2000 replications are shown on interior nodes.

	DISCUSSION

TOP ABSTRACT MATERIALS AND METHODS RESULTS DISCUSSION LITERATURE CITED

In the above, I provided evidence that advantageous amino acid substitutions unrelated to IQ repeats have been fixed by adaptive selection in human ASPM after the human-chimpanzee split, which strongly suggests that ASPM might be an important genetic component in the evolutionary expansion of human brain. The episode of positive selection on ASPM appears to have ended some time ago, as there is no evidence for positive selection on ASPM in current human populations; rather, relatively strong purifying selection is detected. Roughly, selective sweeps occurring in the past 0.5N generations may be detected (FAY and WU 2000 ), where N is the effective population size of humans and is thought to be 10,000 (TAKAHATA et al. 1995 ; HARPENDING et al. 1998 ). That is, the positive selection detected in ASPM occurred some time between 6–7 and 0.1 MY ago (0.5 x 10,000 generations x 20 years/generation). The latter date coincides with the suggested time of migration of modern humans out of Africa (reviewed in CAVALLI-SFORZA and FELDMAN 2003 ). It is also interesting to note that although the precise time when positive selection acted on ASPM is difficult to pinpoint, my estimate is consistent with the current understanding that the human brain expansion took place between 2–2.5 and 0.2–0.4 MY ago (MCHENRY 1994 ; WOOD and COLLARD 1999 ). Furthermore, a selective sweep in human FOXP2, a gene involved in speech and language development, has been detected (ENARD et al. 2002 ; ZHANG et al. 2002 ). This sweep was estimated to have occurred no earlier than 0.1–0.2 MY ago (ENARD et al. 2002 ; ZHANG et al. 2002 ). That is, the adaptive evolution of FOXP2 postdated that of ASPM, consistent with the common belief that a big brain may be a prerequisite for language (DECAN 1992 ).

Studies of ASPM in model organisms can help us understand how it impacts brain size. The mouse Aspm is highly expressed in the embryonic brain, particularly during cerebral cortical neurogenesis (BOND et al. 2002 ). The fruit fly ortholog asp is involved in organizing and binding together microtubules at the spindle poles and in forming the central mitotic spindle (GONZALEZ et al. 1990 ; WAKEFIELD et al. 2001 ). Mutations in asp cause dividing neuroblasts to arrest in metaphase, resulting in reduced central nervous system development (WAKEFIELD et al. 2001 ). The amino acid substitutions in human ASPM are located in exons 3, 18, 20, 21, and 22 (Fig 1), which encode a putative microtubule-binding domain and an IQ calmodulin-binding domain (BOND et al. 2002 ). These features suggest that the adaptive substitutions in human ASPM might be related to the regulation of mitosis in the nervous system, which can be tested in the future by functional assays of human ASPM as well as a laboratory-reconstructed ASPM protein of the common ancestor of humans and chimpanzees.

	FOOTNOTES

Sequence data from this article have been deposited with the EMBL/GenBank Data Libraries under accession nos. AY367065–87.

	ACKNOWLEDGMENTS

I thank David Webb for technical assistance and Douglas Futuyma, Priscilla Tucker, and David Webb for valuable comments on an earlier version of the manuscript. This work was supported by a start-up fund of the University of Michigan and a research grant from the National Institutes of Health (GM67030).

Manuscript received July 7, 2003; Accepted for publication August 20, 2003.

	LITERATURE CITED

TOP ABSTRACT MATERIALS AND METHODS RESULTS DISCUSSION LITERATURE CITED

BOND, J., E. ROBERTS, G. H. MOCHIDA, D. J. HAMPSHIRE, and S. SCOTT et al., 2002  ASPM is a major determinant of cerebral cortical size. Nat. Genet. 32:316-320.[Medline]

BRITTEN, R. J., 2002  Divergence between samples of chimpanzee and human DNA sequences is 5%, counting indels. Proc. Natl. Acad. Sci. USA 99:13633-13635.[Abstract/Free Full Text]

BRUNET, M., F. GUY, D. PILBEAM, H. T. MACKAYE, and A. LIKIUS et al., 2002  A new hominid from the Upper Miocene of Chad, Central Africa. Nature 418:145-151.

CAVALLI-SFORZA, L. L. and M. W. FELDMAN, 2003  The application of molecular genetic approaches to the study of human evolution. Nat. Genet. 33(Suppl.):266-275.

CHEN, F. C. and W. H. LI, 2001  Genomic divergences between humans and other hominoids and the effective population size of the common ancestor of humans and chimpanzees. Am. J. Hum. Genet. 68:444-456.[Medline]

DECAN, T. W., 1992 Biological aspects of language, pp. 128–133 in The Cambridge Encyclopedia of Human Evolution, edited by S. JONES, R. MARTIN and D. PILBEAM. Cambridge University Press, Cambridge, UK.

ENARD, W., M. PRZEWORSKI, S. E. FISHER, C. S. LAI, and V. WIEBE et al., 2002  Molecular evolution of FOXP2, a gene involved in speech and language. Nature 418:869-872.[Medline]

EYRE-WALKER, A., 2002  Changing effective population size and the McDonald-Kreitman test. Genetics 162:2017-2024.[Abstract/Free Full Text]

FAY, J. C. and C.-I WU, 2000  Hitchhiking under positive Darwinian selection. Genetics 155:1405-1413.[Abstract/Free Full Text]

FAY, J. C., G. J. WYCKOFF, and C.-I WU, 2001  Positive and negative selection on the human genome. Genetics 158:1227-1234.[Abstract/Free Full Text]

FAY, J. C., G. J. WYCKOFF, and C.-I WU, 2002  Testing the neutral theory of molecular evolution with genomic data from Drosophila. Nature 415:1024-1026.[Medline]

FELSENSTEIN, J., 1985  Confidence limits on phylogenies: an approach using the bootstrap. Evolution 39:783-791.

FU, Y. X. and W. H. LI, 1993  Statistical tests of neutrality of mutations. Genetics 133:693-709.[Abstract]

GONZALEZ, C., R. D. SAUNDERS, J. CASAL, I. MOLINA, and M. CARMENA et al., 1990  Mutations at the asp locus of Drosophila lead to multiple free centrosomes in syncytial embryos, but restrict centrosome duplication in larval neuroblasts. J. Cell Sci. 96:605-616.[Abstract/Free Full Text]

HARPENDING, H. C., M. A. BATZER, M. GURVEN, L. B. JORDE, and A. R. ROGERS et al., 1998  Genetic traces of ancient demography. Proc. Natl. Acad. Sci. USA 95:1961-1967.[Abstract/Free Full Text]

KIMURA, M., 1980  A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J. Mol. Evol. 16:111-120.[Medline]

KONG, A., D. F. GUDBJARTSSON, J. SAINZ, G. M. JONSDOTTIR, and S. A. GUDJONSSON et al., 2002  A high-resolution recombination map of the human genome. Nat. Genet. 31:241-247.[Medline]

KAESSMANN, H., V. WIEBE, G. WEISS, and S. PAABO, 2001  Great ape DNA sequences reveal a reduced diversity and an expansion in humans. Nat. Genet. 27:155-156.[Medline]

KUMAR, A., M. MARKANDAYA, and S. C. GIRIMAJI, 2002  Primary microcephaly: microcephalin and ASPM determine the size of the human brain. J. Biosci. 27:629-632.[Medline]

KUMAR, S., K. TAMURA, I. B. JAKOBSEN, and M. NEI, 2001  MEGA2: molecular evolutionary genetics analysis software. Bioinformatics 17:1244-1245.[Abstract/Free Full Text]

LI, W. H., 1997 Molecular Evolution. Sinauer, Sunderland, MA.

MCDONALD, J. H. and M. KREITMAN, 1991  Adaptive protein evolution at the Adh locus in Drosophila. Nature 351:652-654.[Medline]

MCHENRY, H. M., 1994  Tempo and mode in human evolution. Proc. Natl. Acad. Sci. USA 91:6780-6786.[Abstract/Free Full Text]

MOCHIDA, G. H. and C. A. WALSH, 2001  Molecular genetics of human microcephaly. Curr. Opin. Neurol. 14:151-156.[Medline]

MURPHY, W. J., E. EIZIRIK, S. J. O'BRIEN, O. MADSEN, and M. SCALLY et al., 2001  Resolution of the early placental mammal radiation using Bayesian phylogenetics. Science 294:2348-2351.[Abstract/Free Full Text]

NEI, M., and S. KUMAR, 2000 Molecular Evolution and Phylogenetics. Oxford University Press, New York.

PTAK, S. E. and M. PRZEWORSKI, 2002  Evidence for population growth in humans is confounded by fine-scale population structure. Trends Genet. 18:559-563.[Medline]

ROZAS, J. and R. ROZAS, 1999  DnaSP version 3: an integrated program for molecular population genetics and molecular evolution analysis. Bioinformatics 15:174-175.[Abstract/Free Full Text]

SAITOU, N. and M. NEI, 1987  The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 4:406-425.[Abstract]

SATTA, Y., J. KLEIN, and N. TAKAHATA, 2000  DNA archives and our nearest relative: the trichotomy problem revisited. Mol. Phylogenet. Evol. 14:259-275.[Medline]

SMITH, N. G. and A. EYRE-WALKER, 2002  Adaptive protein evolution in Drosophila. Nature 415:1022-1024.[Medline]

TAJIMA, F., 1989  Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics 123:585-595.[Abstract/Free Full Text]

TAKAHATA, N., Y. SATTA, and J. KLEIN, 1995  Divergence time and population size in the lineage leading to modern humans. Theor. Popul. Biol. 48:198-221.[Medline]

WAKEFIELD, J. G., S. BONACCORSI, and M. GATTI, 2001  The drosophila protein asp is involved in microtubule organization during spindle formation and cytokinesis. J. Cell Biol. 153:637-648.[Abstract/Free Full Text]

WOOD, B. and M. COLLARD, 1999  The human genus. Science 284:65-71.[Abstract/Free Full Text]

YANG, Z., 1997  PAML: a program package for phylogenetic analysis by maximum likelihood. Comput. Appl. Biosci. 13:555-556.[Free Full Text]

YI, S., D. L. ELLSWORTH, and W. H. LI, 2002  Slow molecular clocks in Old World monkeys, apes, and humans. Mol. Biol. Evol. 19:2191-2198.[Abstract/Free Full Text]

ZHANG, J. and M. NEI, 1997  Accuracies of ancestral amino acid sequences inferred by the parsimony, likelihood, and distance methods. J. Mol. Evol. 44(Suppl. 1):S139-S146.

ZHANG, J. and D. M. WEBB, 2003  Evolutionary deterioration of the vomeronasal pheromone transduction pathway in catarrhine primates. Proc. Natl. Acad. Sci. USA 100:8337-8341.[Abstract/Free Full Text]

ZHANG, J., S. KUMAR, and M. NEI, 1997  Small-sample tests of episodic adaptive evolution: a case study of primate lysozymes. Mol. Biol. Evol. 14:1335-1338.[Medline]

ZHANG, J., H. F. ROSENBERG, and M. NEI, 1998  Positive Darwinian selection after gene duplication in primate ribonuclease genes. Proc. Natl. Acad. Sci. USA 95:3708-3713.[Abstract/Free Full Text]

ZHANG, J., D. M. WEBB, and O. PODLAHA, 2002  Accelerated protein evolution and origins of human-specific features: Foxp2 as an example. Genetics 162:1825-1835.[Abstract/Free Full Text]

This article has been cited by other articles:

				Bin Liu Association of the Dusp6 (Mkp3) Gene With Mouse Brain Weight and Forebrain Structure J Child Neurol, June 1, 2008; 23(6): 624 - 627. [Abstract] [PDF]

				J.-k. Wang, Y. Li, and B. Su A common SNP of MCPH1 is associated with cranial volume variation in Chinese population Hum. Mol. Genet., May 1, 2008; 17(9): 1329 - 1335. [Abstract] [Full Text] [PDF]

				H. L. Kim and Y. Satta Population Genetic Analysis of the N-Acylsphingosine Amidohydrolase Gene Associated With Mental Activity in Humans Genetics, March 1, 2008; 178(3): 1505 - 1515. [Abstract] [Full Text] [PDF]

				N. Takahata Molecular Clock: An Anti-neo-Darwinian Legacy Genetics, May 1, 2007; 176(1): 1 - 6. [Full Text] [PDF]

				M. A. Bakewell, P. Shi, and J. Zhang More genes underwent positive selection in chimpanzee evolution than in human evolution PNAS, May 1, 2007; 104(18): 7489 - 7494. [Abstract] [Full Text] [PDF]

				N. Mekel-Bobrov, D. Posthuma, S. L. Gilbert, P. Lind, M. F. Gosso, M. Luciano, S. E. Harris, T. C. Bates, T. J.C. Polderman, L. J. Whalley, et al. The ongoing adaptive evolution of ASPM and Microcephalin is not explained by increased intelligence Hum. Mol. Genet., March 15, 2007; 16(6): 600 - 608. [Abstract] [Full Text] [PDF]

				D. Falk, C. Hildebolt, K. Smith, M. J. Morwood, T. Sutikna, Jatmiko, E. W. Saptomo, H. Imhof, H. Seidler, and F. Prior Brain shape in human microcephalics and Homo floresiensis PNAS, February 13, 2007; 104(7): 2513 - 2518. [Abstract] [Full Text] [PDF]

				E. J. Vallender and B. T. Lahn A primate-specific acceleration in the evolution of the caspase-dependent apoptosis pathway Hum. Mol. Genet., October 15, 2006; 15(20): 3034 - 3040. [Abstract] [Full Text] [PDF]

				I. Dupanloup and H. Kaessmann Evolutionary simulations to detect functional lineage-specific genes Bioinformatics, August 1, 2006; 22(15): 1815 - 1822. [Abstract] [Full Text] [PDF]

				J. L. Fish, Y. Kosodo, W. Enard, S. Paabo, and W. B. Huttner Aspm specifically maintains symmetric proliferative divisions of neuroepithelial cells PNAS, July 5, 2006; 103(27): 10438 - 10443. [Abstract] [Full Text] [PDF]

				R. P. Woods, N. B. Freimer, J. A. De Young, S. C. Fears, N. L. Sicotte, S. K. Service, D. J. Valentino, A. W. Toga, and J. C. Mazziotta Normal variants of Microcephalin and ASPM do not account for brain size variability Hum. Mol. Genet., June 15, 2006; 15(12): 2025 - 2029. [Abstract] [Full Text] [PDF]

				C. P. Ponting A novel domain suggests a ciliary function for ASPM, a brain size determining gene Bioinformatics, May 1, 2006; 22(9): 1031 - 1035. [Abstract] [Full Text] [PDF]

				A. Varki and T. K. Altheide Comparing the human and chimpanzee genomes: Searching for needles in a haystack Genome Res., December 1, 2005; 15(12): 1746 - 1758. [Abstract] [Full Text] [PDF]

				P. D. Evans, S. L. Gilbert, N. Mekel-Bobrov, E. J. Vallender, J. R. Anderson, L. M. Vaez-Azizi, S. A. Tishkoff, R. R. Hudson, and B. T. Lahn Microcephalin, a Gene Regulating Brain Size, Continues to Evolve Adaptively in Humans Science, September 9, 2005; 309(5741): 1717 - 1720. [Abstract] [Full Text] [PDF]

				N. Mekel-Bobrov, S. L. Gilbert, P. D. Evans, E. J. Vallender, J. R. Anderson, R. R. Hudson, S. A. Tishkoff, and B. T. Lahn Ongoing Adaptive Evolution of ASPM, a Brain Size Determinant in Homo sapiens Science, September 9, 2005; 309(5741): 1720 - 1722. [Abstract] [Full Text] [PDF]

				N. Kouprina, A. Pavlicek, N. K. Collins, M. Nakano, V. N. Noskov, J.-I. Ohzeki, G. H. Mochida, J. I. Risinger, P. Goldsmith, M. Gunsior, et al. The microcephaly ASPM gene is expressed in proliferating tissues and encodes for a mitotic spindle protein Hum. Mol. Genet., August 1, 2005; 14(15): 2155 - 2165. [Abstract] [Full Text] [PDF]

				Y.-q. Wang, Y.-p. Qian, S. Yang, H. Shi, C.-h. Liao, H.-K. Zheng, J. Wang, A. A. Lin, L. L. Cavalli-Sforza, P. A. Underhill, et al. Accelerated Evolution of the Pituitary Adenylate Cyclase-Activating Polypeptide Precursor Gene During Human Origin Genetics, June 1, 2005; 170(2): 801 - 806. [Abstract] [Full Text] [PDF]

				G. H. Mochida Topical Review: Cortical Malformation and Pediatric Epilepsy: A Molecular Genetic Approach J Child Neurol, April 1, 2005; 20(4): 300 - 303. [Abstract] [PDF]

				E. J. Vallender and B. T. Lahn Positive selection on the human genome Hum. Mol. Genet., October 1, 2004; 13(suppl_2): R245 - R254. [Abstract] [Full Text] [PDF]

P. D. Evans, J. R. Anderson, E. J. Vallender, S. S. Choi, and B. T. Lahn Reconstructing the evolutionary history of microcephalin, a gene controlling human brain size Hum. Mol. Genet., June 1, 2004; 13(11): 1139 - 1145