Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Dating the Origin of Language Using Phonemic Diversity

  • Charles Perreault ,

    Contributed equally to this work with: Charles Perreault, Sarah Mathew

    Affiliation Santa Fe Institute, Santa Fe, New Mexico, United States of America

  • Sarah Mathew

    Contributed equally to this work with: Charles Perreault, Sarah Mathew

    sarah.mathew@arklab.su.se

    Affiliation Centre for the Study of Cultural Evolution, Stockholm University, Stockholm, Sweden

Abstract

Language is a key adaptation of our species, yet we do not know when it evolved. Here, we use data on language phonemic diversity to estimate a minimum date for the origin of language. We take advantage of the fact that phonemic diversity evolves slowly and use it as a clock to calculate how long the oldest African languages would have to have been around in order to accumulate the number of phonemes they possess today. We use a natural experiment, the colonization of Southeast Asia and Andaman Islands, to estimate the rate at which phonemic diversity increases through time. Using this rate, we estimate that present-day languages date back to the Middle Stone Age in Africa. Our analysis is consistent with the archaeological evidence suggesting that complex human behavior evolved during the Middle Stone Age in Africa, and does not support the view that language is a recent adaptation that has sparked the dispersal of humans out of Africa. While some of our assumptions require testing and our results rely at present on a single case-study, our analysis constitutes the first estimate of when language evolved that is directly based on linguistic data.

Introduction

A capacity for language is a hallmark of our species [1], [2], yet we know little about the timing of its appearance. Language appears in the archaeological record only recently, with the advent of lexicographic writing around 5,400 years ago [3]. Therefore, investigators have addressed the origin of language by studying the evolutionary history of anatomical features [4][7] and genes [8][15] that are associated with speech production. This research suggests that other Homo species had the ability to produce speech sounds that overlap with the range of speech sounds of modern humans, and that species such as Neanderthals possessed genes that, in humans, play a role in language. But we do not know whether these archaic hominins actually produced speech, and if so, to which extent it was similar to our capacity for language. As of now, the anatomical and genetic data lack the resolution necessary to differentiate proto-language from modern human language. Until this resolution is improved, we need alternative lines of evidence in order to better understand the timing of language origin.

Here, we use phonemic diversity data to date the origin of language. Phonemic diversity denotes the number of perceptually distinct units of sound–consonants, vowels and tones–in a language. The worldwide pattern of phonemic diversity potentially contains the statistical signal of the expansion of modern humans on the planet [16]. As human populations left Africa, 60–70 kya, and expanded into the rest of the world [1], [17], they underwent a series of bottlenecks. This serial founder effect has led to a clinal loss of genetic [18][20], phenotypic [21][23] and phonemic diversity [16] that can be observed in present-day human populations. African languages today have some of the largest phonemic inventories in the world, while the smallest inventories are found in South America and Oceania, some of the last regions of the globe to be colonized. The loss of phonemes through serial founder effect is consistent with other lines of evidence that indicate that phonemic diversity is determined by cultural transmission forces, rather than cognitive or functional constraints. First, phonemic diversity varies considerably among languages, and several languages function with a restricted number of phonemes. Rotokas, a language of New Guinea, and Pirahã, spoken in South-America, both have 11 phonemes [24], [25], while !Xun, a language spoken in Southern Africa has 141 phonemes. Second, as predicted by theoretical models linking cultural transmission and demography [26][28], phonemic diversity correlates positively with speaker population size [16], [29]. And finally, phonemic diversity also correlates positively with the number of surrounding languages [16], suggesting that phonemes, like other cultural traits, can be borrowed. Phonemic diversity not only evolves culturally, but it also evolves slowly [16]. That the languages outside of Africa might have not recovered their original phonemic diversity, despite thousands of years of history in their respective continent, and despite all the historical, linguistic and social factors that lead to linguistic change [30][36], suggests that phonemic diversity changes over long time scales. Here, we take advantage of the fact that phonemic diversity evolves culturally and slowly, and use it as a slow-clock to date the origin of language.

By focusing on phonemes rather than cognates–words that share a common ancestry–we are able to circumvent problems that prevent current historical linguistic approaches from tackling the problem of dating the origin of language. Glottochronology uses the number of cognates that languages share to estimate when they diverged [37][39]. However, because cognates change over short time scales, the time-depth resolution of glottochronology is limited to a few thousand years [8]. Several historical, social and demographic factors influence cognate evolution, [31], [32], [34], [40], [41], a main one being frequency of word use. Common words evolve more slowly than rare ones [42]. Frequency of word-use alone predicts 50% of the variation in rates of cognate change, and can generate cognate half-lives that range from 750 years to more than 10,000 years [42]. Such variation in rates of cognate change is problematic for glottochronology, because glottochronology assumes a constant rate of cognate change [43], [44]. The assumption of a constant rate of change can be relaxed by applying phylogenetic methods to cognate datasets. These methods are powerful tools for estimating the date of divergence of language families [45][47]. Nonetheless, the temporal scope of this method is, at least in its current state, too limited to address questions about the origin of language. For instance, the average word half-life among Indo-European languages is about 2,530 years [42]. Here we circumvent the problem of variation in rates by averaging rates of phoneme accumulation over a large spatial and temporal scale.

Given that languages accumulate phonemes over long time scales, we ask how long African languages had to have been around in order to reach their current phonemic diversity. We start by building two related mathematical models that describe two ways by which phonemic diversity can rise through time. In the first model, phonemic inventory increases linearly with time, while in the second model phonemic inventory increases exponentially. Then, we parametrize the two models with empirical data. Finally, we use rewritten forms of the models to estimate the time span over which phonemes would have had to accumulate in Africa.

We do not attempt to capture all the factors that influence phonemic inventory size. The state of our knowledge does not allow us to formalize the specific mechanisms by which phonemic diversity increases and decreases. Therefore, our models are agnostic about the particular mechanisms of change in phonemic diversity, and capture only the net effect of these mechanisms on phonemic diversity. We summarize this net effect as a single number, a rate of phoneme accumulation through time. Note that phonemic changes that occur within a language and that do not lead to a net change in the size of the phonemic inventory are not relevant to our analysis. The crucial assumption underlying our models is that the net effect of the factors leading to phonemic gain is greater than the net effect of those leading to loss. When this assumption is met, all other things being equals, phonemic diversity increases through time.

The method used in this paper to date the origin of language is built upon various assumptions that require further testing. An assumption underlying the empirical parametrization of the model is that human populations have lost phonemes through a drift-loss process during their expansion across the world [16]. However, this hypothesis is not widely accepted among linguists. Problems with the drift-loss hypothesis are discussed in a collection of commentaries published in Linguistic Typology [48][60] and Science [61][63]. Overall, these commentaries highlight the fact that, while Atkinson’s hypothesis remains viable, alternative hypothesis to the worldwide pattern of phonemic diversity have yet to be satisfyingly rejected [64], [65]. As we describe our method and material below, we specify the other assumptions that we have made and that also need further investigation to be validated. Despite these caveats, our approach constitutes a novel solution to the difficult question of dating the origin of language.

Analysis

We start by estimating the rate at which languages accumulate phonemes. Controlling for distance from Africa, the phonemic diversity of a language depends on the speaker population size, the geographic area over which the language is spoken, and local linguistic diversity [16]. This suggests that new phonemes are more likely to appear in large populations. It also suggests that phonemes can be borrowed through contact between groups and languages [16].

With that in mind, consider the hypothetical case of two small populations, B and C, that dispersed from the same parent population, A, t years ago (Figure 1). Suppose that B and C are similar in size so that they both experience approximately the same loss in phonemic diversity due to the founder effect. Now, suppose that population B colonizes a large continental territory and subsequently expands and diversifies linguistically [66], [67]. In contrast, population C settles on a small island that does not allow for population expansion and language diversification. Because of the differences between the regions colonized by B and C, population B will accumulate phonemes at a faster rate than population C. Furthermore, if population C evolves on a sufficiently small island and remains isolated for most of its history, then the rate of phoneme accumulation in C will be low, and its phonemic diversity will remain approximately stable through time. Consequently, the present-day difference between the phonemic diversity of B and C can be attributed to the new phonemes accumulated within population B. Thus, the current phonemic diversity of population C has remained through time a good approximation of the original phonemic diversity of population B. When this is true, and if the date of colonization, t, is known, then it is possible to estimate the phoneme accumulation rate in a large population as(1)assuming that phonemic inventories increase linearly, and(2)assuming that phonemic inventories increase exponentially. and are the current phonemic diversity of populations B and C, and is the time elapsed between divergence of B and C, and the moment when their present phonemic inventories were measured. The linear model (Equation 1) is appropriate when phonemes increase independently of a language’s phonemic diversity. The exponential model (Equation 2) captures the alternative situation where the rate at which phonemes accumulate increases with a language’s phonemic diversity. Such dependence would arise, for instance, if each phoneme has the potential to give rise to new phonemes.

thumbnail
Figure 1. A model of change in phonemic diversity through drift and recovery.

At time two small populations, B and C, emigrate from population and colonize two different regions. Population B settles on a large landmass, and subsequently grows and diversifies linguistically. As a result, the average phonemic diversity of population B increases with time. Conversely, the phonemic diversity of population C remains stable through time because it occupies a small, isolated island. Therefore, the phonemic diversity of population C can be used to approximate what the phonemic diversity of population B would have been at time Large dots denote high phonemic diversity and small dots denote low phonemic diversity.

https://doi.org/10.1371/journal.pone.0035289.g001

To estimate and empirically, we take advantage of a natural experiment that approximates the scenario outlined in Figure 1, the migration history of humans in mainland Southeast Asia and the Andaman Islands. Both Southeast Asia and the Andaman Islands were colonized during the Pleistocene dispersal of modern humans out of Africa, a process that started 70–60 kya [71]. Genetic data indicate that humans dispersed in Asia following a coastal route, from India to Australia [17], [68][70], and that both Southeast Asia and Andaman Islands were colonized from a population that occupied the region spanning from southern India to the Malay Peninsula [69], [71], [72]. This dispersal was rapid. Genetic analyses estimate that it occurred approximately 65 kya [69], [71], and the archaeological record puts humans both in Southeast Asia [73] and Australia [74] at least 45 kya. Relative to the long temporal scale over which phonemes accumulate, we expect that the Andaman Islands and Mainland Southeast Asia were colonized simultaneously.

Populations in Southeast Asia and Andaman Islands differed demographically and linguistically. Like population B above, human groups expanded considerably after their arrival in Southeast Asia. By 40–20 kya, more than half of the total human population is estimated to have lived in South and Southeast Asia [75]. Today, about 160 million people live in Mainland Southeast Asia, and speak more than 60 languages. Conversely, we expect the Andaman population to have mirrored population C in the example above, and to have gained few novel phonemes, because of their low population size and remarkable degree of isolation. The Andaman Islands constitute a fragmented landscape of about 200 small islands, with a carrying capacity estimated to about 5000 individuals before contact with Europeans [76]. Genetic analyses suggest that the inhabitants of Andaman Islands have remained isolated since their arrival during the Pleistocene, up until the mid-19th century [70], [72], [77]. The 13 languages spoken on the islands at that time period are linguistic isolates, with no clear relationship to other Asian languages [78][81].

We estimate the parameters and in Equations 1 and 2 as follows. Assuming that Mainland Southeast Asia and Andaman Islands were colonized at some point in time between 45 kya and 65 kya, we use 45 and 65 k as lower and upper bounds of We obtained the phonemic diversity of languages of Mainland Southeast Asia and Andaman Islands using data from the UCLA Phonological Segment Inventory Database (UPSID) [24], [25]. While the categorical scaled measurements of phonemic diversity of the World Atlas of Language Structures (WALS) [82] were sufficient to detect a potential global serial founder effect [16], they are inadequate for the calculation of a phoneme accumulation rate. The UPSID contains the number of phonemic segments of a global sample of 451 languages. We estimate by taking the average phonemic inventory size of the languages in Mainland Southeast Asia. Assuming an eastward, coastal migration route, we have excluded the Asian languages that are located west of Andaman Islands (such as the languages from India and Nepal), as well as those spoken in Myanmar and the Malay Peninsula, because they could have served as departure points for the colonization of Andaman Islands (Figure 2). The 20 languages retained in our sample are thus those spoken in Cambodia, Vietnam, Laos and Southwest China (Table 1). The average phonemic diversity of the resulting sample is (errors represent one standard error). Great Andamanese (ISO 639-2: apq) is the only Andamanese language to appear in UPSID. Its phonemic diversity, 24, serves as our estimate of

thumbnail
Figure 2. Approximate location of the languages included in the Mainland Southeast Asia sample.

The languages located inside the shaded area were excluded from the sample because the region is a potential departure point for the colonization of Andaman Islands or Mainland Southeast Asia.

https://doi.org/10.1371/journal.pone.0035289.g002

Setting to 41.21 and to 24, we obtain range estimates for the phoneme accumulation parameters and for a large, linguistically diverse population (Table 2). Note that, in the real world, we expect and to vary through time and space, both within and between languages, as a result of various linguistic forces and historical contingencies. In contrast, our estimates of and are averaged over 20 languages, that are dispersed over a vast spatial area, and that have been evolving in the region for perhaps as long as 60 ky. By using a time and space-averaged value, we are attempting to eliminate the effect of local contingencies and estimate the expected value of the rate of phoneme accumulation of human languages. We need such time and space-averaged value especially since we are dating an event that happened thousands of years ago, by using the average present-day phonemic diversity of multiple African languages.

Using the rates of phoneme accumulation and we calculate the time it would take for a language to acquire the phonemic diversity observed today in African languages, (3)or(4)where is the number of phonemes the first human languages started with. Phonemic diversity is assumed to have increased linearly in Equation 3, and exponentially in Equation 4.

To estimate we use the average phonemic diversity of African languages that natively possess clicks [83], [84]. We do so because they comprise the African languages that have had the longest continuous history, and as a result are the ones that have lost phonemes due to founder effect the least recently. The largest language groups in Africa–Afro-Asiatic, Niger-Congo and Nilo-Saharan–underwent recently considerable geographic expansion [85], which could have decreased their phonemic diversity through serial founder effect. This idea is consistent with the fact that the average phonemic diversity of Afro-Asiatic, Niger-Congo and Nilo-Saharan languages is 36, 33, and 29 respectively, while the average phonemic diversity of African languages outside these families is 75. The African languages in UPSID outside of these three families are Hadza, Khoekhoe, Sandawe and !Xun. All of these languages use click consonants. Genetic analyses suggests that the speakers of these languages may have had the longest continuous population history [85][89], with mitochondrial DNA and Y chromosome variation indicating that the divergence between the click language speakers is at least as old as the divergence between any other pair of human populations [85], [86]. The main click language branches–Hadza, Sandawe and South African Khoisan (the last one includes Khoekhoe and !Xun) are estimated to have diverged as early as 55–35 kya [85], [86], with Hadza and Sandawe splitting 20–15 kya [85]. We have also included the Dahalo language in our sample. Dahalo is an Afro-Asiatic language, but the occurrence of click sounds in its core vocabulary suggests that it natively may have had clicks [90]. Using the five African click languages present in UPSID, we estimate to be (Table 3).

We cannot know what the initial number of phonemes of the first human language, was. A reasonable assumption is that it is equal to the smallest phonemic inventory ever observed Therefore, we have set to 11 phonemes. On the other hand, it is possible that the languages with the lowest phonemic diversity today are outliers, and that a central value of the world’s phonemic diversity better approximates the initial phonemic diversity of human languages. We show how changing to the median phonemic diversity of the languages in the UPSID sample affects the result.

Results

When is 45–65 kya, the linear and the exponential growth models yield values of 232–159 kya and 225–156 kya, respectively. Setting to the median phonemic diversity, 29, decreases our estimate to 163–112 kya and 75–108 kya for the linear and exponential growth models respectively. We have also estimated intervals around using one standard error around and the rates of accumulation and The value of is minimized when phonemic diversity in Africa is low and phoneme accumulation rate is high. Conversely, is maximized when phonemic diversity in Africa is high and phoneme accumulation rate is low. Therefore, the upper bound for under linear growth is obtained by setting Equation 2 to and its lower bound is obtained by setting Equation 2 to Similarly, under exponential growth, the upper bound of is and the lower bound is The resulting date ranges are shown in Figure 3.

thumbnail
Figure 3. Date estimates for the origin of present-day languages.

Range estimates are calculated assuming an initial phonemic diversity, of 11 and 29, and a colonization time for Mainland Southeast Asia and Andaman Islands of 45 kya and 65 kya (parameter ). The lozenges denote central values, and the error bars represent calculations made with one standard deviation of the rate of phoneme accumulation (parameters and ), and the current phonemic diversity of African languages as described in the Results section. White and black lozenges represent results for the linear and exponential models of phonemic gain, respectively.

https://doi.org/10.1371/journal.pone.0035289.g003

These estimates are fairly insensitive to changes in model assumptions. We have considered the possibility that we are overestimating the phonemic diversity of African languages by restricting our sample to click languages. Click sounds may be evolving independently of non-click sounds. This would mean that a language could accumulate non-click phonemes at a certain rate, such as and while simultaneously accumulating click-sounds at another rate If this is true, then we cannot compare the African languages in our sample which have been accumulating clicks and non-click phonemes simultaneously, to Mainland Southeast Asian and Andaman languages which do not contain click sounds. To account for this possibility, we excluded click sounds from the phoneme inventory counts of African languages. The average non-click phonemic diversity of our sample of African click languages is Using this value decreases our estimate to 158–108 kya and 187–129 kya for linear and exponential growth respectively. We have also tested the robustness of our results by excluding the Dahalo language from our sample of African languages. While Dahalo is thought to natively possess clicks, it is an Afro-Asiatic language [83] and as such might bias our sample of African languages towards lower phonemic diversity. Removing it from the sample increases our estimate to 244–167 kya and 230–159 kya for linear and exponential growth respectively. Finally, we also have increased our sample of Mainland Southeast Asian languages. Previously, we had excluded the languages spoken in Thailand, Malaysia, and Myanmar, as the colonizers of the Andaman Islands could have departed from one of these regions. By relaxing this assumption and including in our sample all the Mainland Southeast Asian languages contained in UPSID (the languages spoken in Myanmar, Thailand, Malaysia, Laos, Cambodia, Vietnam and Southwest China), we find that the average phonemic diversity in the region, is which increases our estimate to 242–168 kya and 236–163 kya for linear and exponential growth, respectively.

Discussion

Our analysis suggests that language appears early in the history of our species. It does not support the idea that language is a recent adaptation that could have sparked the colonization of the globe by our species about 50 kya [1], [91]. Rather, our result is consistent with the archaeological evidence suggesting that human behavior became increasingly complex during the Middle Stone Age (MSA) in Africa, sometime between 350–150 kya [92][100]. However, we cannot rule out the possibility that other linguistic adaptations, that are independent of phonemic evolution, arose later and triggered the out-of-Africa expansion.

Our date estimate for the origin of language roughly coincides with the date range for the emergence of modern humans. Fossil evidence suggests that anatomically modern humans were present by 195–160 kya [101][104], and fossils classified as Homo helmei, that may be anatomically modern or nearly modern, are dated to 300–250 kya [95], [100]. Coalescence times from genetic data suggest that a genetic population bottleneck, possibly associated with a speciation event, occurred 200–100 kya [85], [105][108].

A population bottleneck causing a loss of phonemes would push back, or even reset the phonemic clock. As a result, our date estimates should be treated as minimum ages for the origin of language. It is thus possible that language arose before the last speciation event in our lineage, or even before the appearance of behavioral modernity.

Our date estimates should be treated with caution. Our results hinge on a series of assumptions in addition to the ones laid out in the Material and Method section. We assume that the rate of phoneme accumulation of Southeast Asia and Africa were similar. We assume that the Andaman languages did not accumulate new phonemes following the colonization of the Andaman Islands, or lose phonemes when their populations crashed upon contact with Europeans. We assume that the founding populations that settled Andaman Islands and Mainland Southeast Asia have lost an equivalent number of phonemes due to drift. Also, the UPSID phoneme counts do not include tonal distinctions. The absence of tonal distinctions in our data could add noise to our analysis, and bias it if it leads us to underestimate the phonemic diversity of one of the continental regions, Africa and Mainland Southeast Asia, more so than the other. We assume that the rate of accumulation of phonemes does not decrease as phonemic inventory size increases. An accumulation rate that decreases with phonemic diversity would lead us to underestimate the antiquity of present-day phonemic inventories. A similar bias would also occur if the phoneme accumulation rate changed through time as our species evolved. Furthermore, our estimate of the rate of phoneme accumulation is based on a single historical case. We are not aware of other colonization sequences that resembles the one outlined in Figure 1 that would also be ancient enough to allow for phonemic inventories to increase. However, despite the caveats we have highlighted here, this analysis constitutes the first appraisal of when language evolved to be based directly on linguistic data.

Acknowledgments

We thank Tanmoy Bhattacharya, Robert Boyd, P. Jeffrey Brantingham, Juergen Neubauer and three anonymous reviewers for insightful comments and criticisms that helped improve this paper.

Author Contributions

Conceived and designed the experiments: SM CP. Performed the experiments: SM CP. Analyzed the data: SM CP. Contributed reagents/materials/analysis tools: SM CP. Wrote the paper: SM CP.

References

  1. 1. Klein RG (2009) The human career. Chicago: University of Chicago Press, 989 pp.
  2. 2. Tattersall I (2009) Human origins: Out of Africa. Proc Natl Acad Sci USA 106: 16018–16021.
  3. 3. Powell BB (2009) Writing: Theory and history of the technology of civilization. Malden, MA: Wiley-Blackwell.
  4. 4. Fitch W (2000) The evolution of speech: a comparative review. Trends in Cognitive Sciences 4: 258–267.
  5. 5. Lieberman DE (2008) Speculations about the selective basis for modern human craniofacial form. Evolutionary Anthropology 17: 55–68.
  6. 6. Lieberman P (2007) The evolution of human speech: Its anatomical and neural bases. Current Anthropology 48: 39–66.
  7. 7. Houghton P (1993) Neandertal supralaryngeal vocal tract. American Journal of Physical Anthropology 90: 139–46.
  8. 8. Campbell L (2000) Time perspective in linguistics. Renfrew C, McMahon A, Trask L, editors, Time Depth in Historical Linguistics, Volume 1, Cambridge: McDonald Institute for Archaeological Research, chapter 1. pp. 3–31.
  9. 9. Enard W, Przeworski M, Fisher SE, Lai CSL, Wiebe V, et al. (2002) Molecular evolution of FOXP2, a gene involved in speech and language. Nature 418: 869–72.
  10. 10. Fisher SE, Marcus GF (2006) The eloquent ape: genes, brains and the evolution of language. Nature Reviews Genetics 7: 9–20.
  11. 11. Fisher SE, Scharff C (2009) FOXP2 as a molecular window into speech and language. Trends in Genetics 25: 166–77.
  12. 12. Zhang J, Webb DM, Podlaha O (2002) Accelerated protein evolution and origins of human-specific features: Foxp2 as an example. Genetics 162: 1825–35.
  13. 13. Lieberman P (2006) The FOXP2 gene, human cognition and language. International congress Series 1296: 115–126.
  14. 14. Coop G, Bullaughey K, Luca F, Przeworski M (2008) The timing of selection at the human FOXP2 gene. Molecular biology and evolution 25: 1257–9.
  15. 15. Krause J, Lalueza-Fox C, Orlando L, Enard W, Green RE, et al. (2007) The derived FOXP2 variant of modern humans was shared with Neandertals. Current Biology 17: 1908–12.
  16. 16. Atkinson QD (2011) Phonemic diversity supports a serial founder effect model of language expansion from Africa. Science 332: 346–349.
  17. 17. Mellars P (2006) Why did modern human populations disperse from Africa ca. 60,000 years ago? A new model. Proc Natl Acad Sci USA 103: 9381–9387.
  18. 18. Li JZ, Absher DM, Tang H, Southwick AM, Casto AM, et al. (2008) Worldwide human relationships inferred from genome-wide patterns of variation. Science (New York, NY) 319: 1100–4.
  19. 19. Prugnolle F, Manica A, Balloux F (2005) Geography predicts neutral genetic diversity of human populations. Current Biology 15: R159–R160.
  20. 20. Ramachandran S, Deshpande O, Roseman CC, Rosenberg NA, Feldman MW, et al. (2005) Support from the relationship of genetic and geographic distance in human populations for a serial founder effect originating in Africa. Proc Natl Acad Sci USA 102: 15942–7.
  21. 21. Betti L, Balloux F, Amos W, Hanihara T, Manica A (2009) Distance from Africa, not climate, explains within-population phenotypic diversity in humans. Proceedings of the Royal Society B Biological Sciences 276: 809–814.
  22. 22. Von Cramon-Taubadel N, Lycett SJ (2008) Brief communication: human cranial variation fits iterative founder effect model with African origin. American Journal of Physical Anthropology 136: 108–113.
  23. 23. Manica A, Amos W, Balloux F, Hanihara T (2007) The effect of ancient population bottlenecks on human phenotypic variation. Nature 448: 346–348.
  24. 24. Maddieson I (1984) Patterns of sounds. Cambridge: Cambridge University Press.
  25. 25. Maddieson I, Precoda K (1990) Updating UPSID. UCLA Working Papers in Phonetics 74: 104–111.
  26. 26. Henrich J (2004) Demography and cultural evolution: How adaptive cultural processes can produce maladaptive losses: The Tasmanian case. American Antiquity 69: 197–214.
  27. 27. Powell A, Shennan S, Thomas MG (2009) Late Pleistocene demography and the appearance of modern human behavior. Science 324: 1298–1301.
  28. 28. Shennan S (2001) Demography and cultural innovation: a model and its implications for the emergence of modern human culture. Cambridge Archaeological Journal 11: 5–16.
  29. 29. Hay J, Bauer L (2007) Phoneme inventory size and population size. Language 83: 388–400.
  30. 30. Croft W (2000) Explaining language change: an evolutionary approach. Oxford, UK: Oxford university Press.
  31. 31. Labov W (1994) Principles of linguistic change: Internal factors. Oxford, UK: Blackwell.
  32. 32. Labov W (2001) Principles of linguistic change: Social factors. Oxford, UK: Blackwell.
  33. 33. Labov W (2010) Principles of linguistic change: Cognitive and cultural factors.
  34. 34. Thomason SG, Kaufman T (1988) Language contact, creolization, and genetic linguistics. Berkeley, CA: University of California Press.
  35. 35. Yang CD (2003) Knowledge and learning in natural language. Oxford: Oxford University Press.
  36. 36. Trudgill P (2004) Linguistic and social typology: The Austronesian migrations and phoneme inventories. Linguistic Typology 8: 305–320.
  37. 37. Swadesh M (1952) Lexicostatistic dating of prehistoric ethnic contacts. Proceedings American Philosophical Society 96: 452–463.
  38. 38. Swadesh M (1955) Towards greater accuracy in lexicostatistic dating. International Journal of American Linguistics 21: 121–137.
  39. 39. Lees RB (1953) The basis of glottochronology. Language 29: 113–127.
  40. 40. Milroy J, Milroy L (1985) Linguistic change, social network and speaker innovation. Journal of Linguistics 21: 339–384.
  41. 41. Nettle D (1999) Is the rate of linguistic change constant? Lingua 108: 119–136.
  42. 42. Pagel M, Atkinson QD, Meade A (2007) Frequency of word-use predicts rates of lexical evolution throughout Indo-European history. Nature 449: 717–20.
  43. 43. Bergsland K, Vogt H (1962) On the validity of glottochronology. Current Anthropology 3: 115.
  44. 44. Blust R, Renfrew C, McMahon A, Trask RL (2001) Why lexicostatistics doesn’t work: the “universal constant” hypothesis and the Austronesian languages. Renfrew C, McMahon A, Trask L, editors, Time depth in historical linguistics, McDonald Institute for Archaeological Research, volume 2, chapter 13. pp. 311–331.
  45. 45. Gray RD, Atkinson QD (2003) Language-tree divergence times support the Anatolian theory of Indo-European origin. Nature 426: 435–439.
  46. 46. Gray RD, Drummond AJ, Greenhill SJ (2009) Language phylogenies reveal expansion pulses and pauses in Pacific settlement. Science 323: 479–483.
  47. 47. Gray RD, Atkinson QD, Greenhill SJ (2011) Language evolution and human history: what a difference a date makes. Philosophical Transactions of the Royal Society B: Biological Sciences 366: 1090–1100.
  48. 48. Bybee J (2011) How plausible is the hypothesis that population size and dispersal are related to phoneme inventory size? Introducing and commenting on a debate. Linguistic Typology 15: 147–153.
  49. 49. Trudgill P (2011) Social structure and phoneme inventories. Linguistic Typology 15: 155–160.
  50. 50. Donohue M, Nichols J (2011) Does phoneme inventory size correlate with population size? Linguistic Typology 15: 161–170.
  51. 51. Dahl O (2011) Are small languages more or less complex than big ones? Linguistic Typology 15: 171–175.
  52. 52. Wichmann , Rama T, Holman , Eric W (2011) Phonological diversity, word length, and population sizes across languages: The ASJP evidence. Linguistic Typology 15: 177–197.
  53. 53. Sproat R (2011) Phonemic diversity and the out-of-Africa theory. Linguistic Typology 15: 199–206.
  54. 54. Bowern C (2011) Out of Africa? The logic of phoneme inventories and founder effects. Linguistic Typology 15: 207–216.
  55. 55. Pericliev V (2011) On phonemic diversity and the origin of language in Africa. Linguistic Typology 15: 217–221.
  56. 56. Ringe D (2011) A pilot study for an investigation into Atkinson’s hypothesis. Linguistic Typology 15: 223–231.
  57. 57. Rice K (2011) Athabaskan languages and serial founder effects. Linguistic Typology 15: 233–250.
  58. 58. Ross B, Donahue M (2011) The many origins of diversity and complexity in phonology. Linguistic Typology 15: 251–265.
  59. 59. Maddieson I, Bhattacharya T, Smith DE, Croft W (2011) Geographical distribution of phonological complexity. Linguistic Typology 15: 267–279.
  60. 60. Jaeger TF, Graff P, Croft W, Pontillo D (2011) Mixed e_ect models for genetic and areal dependencies in linguistic typology. 2011 15: 281–319.
  61. 61. Cysouw M, Dediu D, Moran S (2012) Comment on “Phonemic diversity supports a serial founder effect model of language expansion from Africa”. Science 335: 657.
  62. 62. Wang C, Ding Q, Tao H, Li H (2012) Comment on “Phonemic Diversity Supports a Serial Founder Effect Model”. Science 335: 657.
  63. 63. Tuyl RV, Pereltsvaig A (2012) Comment on “Phonemic Diversit’y Supports a Serial Founder Effect Model”. Science 335: 657.
  64. 64. Atkinson QD (2011) Linking spatial patterns of language variation to ancient demography and population migrations. Linguistic Typology 15: 321–332.
  65. 65. Atkinson QD (2012) Response to Comments on \Phonemic Diversity Supports a Serial Founder Effect Model of Expansion from Africa". Science 335: 657.
  66. 66. Lightfoot D (2006) How new languages emerge. Cambridge, UK: Cambridge University Press, 199 pp.
  67. 67. Nettle D (1998) Explaining global patterns of language diversity. Journal of Anthropological Archaeology 17: 354–374.
  68. 68. Forster P, Matsumura S (2006) Did early humans go north or south? Science 308: 965–966.
  69. 69. Macaulay V, Hill C, Achilli A, Rengo C, Clarke D, et al. (2005) Single, rapid coastal settlement of Asia revealed by analysis of complete mitochondrial genomes. Science 308: 1034–6.
  70. 70. Thangaraj K, Chaubey G, Kivisild T, Reddy AG, Singh VK, et al. (2005) Reconstructing the origin of Andaman Islanders. Science 308: 996.
  71. 71. Palanichamy MG, Sun C, Agrawal S, Bandelt HJ, Kong QP, et al. (2004) Phylogeny of mitochondrial DNA macrohaplogroup N in India, based on complete sequencing: implications for the peopling of South Asia. American Journal of Human Genetics 75: 966–78.
  72. 72. Endicott P, Gilbert MTP, Stringer C, Lalueza-Fox C, Willerslev E, et al. (2003) The genetic origins of the Andaman Islanders. American Journal of Human Genetics 72: 178–84.
  73. 73. Barker G (2005) The archaeology of foraging and farming at Niah Cave, Sarawak. Asian Perspectives 44: 90–106.
  74. 74. O’Connell JF, Allen J (2004) Dating the colonization of Sahul (Pleistocene Australia-New Guinea): a review of recent research. Journal of Archaeological Science 31: 835–853.
  75. 75. Atkinson QD, Gray RD, Drummond AJ (2008) mtDNA variation predicts population size in humans and reveals a major Southern Asian chapter in human prehistory. Molecular biology and evolution 25: 468–74.
  76. 76. Ericksen P, Beckerman S (1975) Population determinants in the Andaman Islands. Mankind 10: 105–107.
  77. 77. Thangaraj K, Singh L, Reddy AG, Rao VR, Sehgal SC, et al. (2003) Genetic affinities of the Andaman Islanders, a vanishing human population. Current Biology 13: 86–93.
  78. 78. Greenberg JH (1971) The Indo-Pacific hypothesis. Sebeok TA, editor, Current Trends in Linguistics vol 8: Linguistics in Oceania, Mouton de Gruyter, volume 8. pp. 809–871.
  79. 79. Wurm SA (1971) Classifications of Australian languages, including Tasmanian. Sebeok TA, editor, Current trends in linguistics, vol. 8: Linguistics in Oceania, Mouton de Gruyter. pp. 721–778.
  80. 80. Wurm SA, McElhanan K (1975) New Guinea area languages and language study, vol 1: Papuan languages and the New Guinea linguistic scene. Canberra: Australian National University.
  81. 81. Ruhlen M (1987) A guide to the world’s languages. Stanford, CA: Stanford University Press.
  82. 82. Dryer MS, Haspelmath M (2011) The World Atlas of Language Structures online. Munich: Max Planck Digital Library.
  83. 83. Güldemann T, Stoneking M (2008) A historical appraisal of clicks: A linguistic and genetic population perspective. Annual Review of Anthropology 37: 93–109.
  84. 84. Greenberg JH (1963) The languages of Africa. Bloomington: Indiana University Press.
  85. 85. Tishkoff SA, Gonder MK, Henn BM, Mortensen H, Knight A, et al. (2007) History of clickspeaking populations of Africa inferred from mtDNA and Y chromosome genetic variation. Molecular biology and evolution 24: 2180–95.
  86. 86. Knight A, Underhill PA, Mortensen HM, Zhivotovsky LA, Lin AA, et al. (2003) African Y chromosome and mtDNA divergence provides insight into the history of click languages. Current Biology 13: 464–473.
  87. 87. Semino O, Santachiara-Benerecetti AS, Falaschi F, Cavalli-Sforza LL, Underhill PA (2002) Ethiopians and Khoisan share the deepest clades of the human Y-chromosome phylogeny. American Journal of Human G’enetics 70: 265–8.
  88. 88. Chen YS, Olckers A, Schurr TG, Kogelnik AM, Huoponen K, et al. (2000) mtDNA variation in the South African Kung and Khwe-and their genetic relationships to other African populations. American journal of human genetics 66: 1362–83.
  89. 89. Watson E, Bauer K, Aman R, Weiss G, von Haeseler A, et al. (1996) mtDNA sequence diversity in Africa. American Journal of Human Genetics 59: 437–44.
  90. 90. Sands B, Güldemann T (2009) What click languages can and can’t tell us about language origins. Botha RP, Knight C, editors, The cradle of language (studies in the evolution of language), USA: Oxford University Press, chapter 11. pp. 204–218.
  91. 91. Diamond J (1991) The rise and fall of the third chimpanzee. Hutchinson Radius (Vintage Edition 1992), 360 pp.
  92. 92. Brown KS, Marean CW, Herries AIR, Jacobs Z, Tribolo C, et al. (2009) Fire as an engineering tool of early modern humans. Science 325: 859–62.
  93. 93. D’Errico F, Henshilwood C, Vanhaeren M, Van Niekerk K (2005) Nassarius kraussianus shell beads from Blombos Cave: evidence for symbolic behaviour in the Middle Stone Age. Journal of Human Evolution 48: 3–24.
  94. 94. D’Errico F, Henshilwood CS (2007) Additional evidence for bone technology in the southern African Middle Stone Age. Journal of Human Evolution 52: 142–163.
  95. 95. Foley R, Lahr MM (2003) On stony ground: Lithic technology, human evolution, and the emergence of culture. Evolutionary Anthropology 12: 109–122.
  96. 96. Henshilwood CS, Marean CW (2003) The origin of modern human behavior: Critique of the models and their test implications. Current Anthropology 44: 627–651.
  97. 97. Marean CA, Assefa Z (2005) The Middle and Upper Pleistocene African record for the biological and behavioral origins of modern humans. Stahl A, editor, African archaeology: A critical introduction, Malden, MA: Blackwell Publishing, chapter 4. pp. 93–129.
  98. 98. Marean CW, Bar-Matthews M, Bernatchez J, Fisher E, Goldberg P, et al. (2007) Early human use of marine resources and pigment in South Africa during the Middle Pleistocene. Nature 449: 905–908.
  99. 99. Marean CW (2011) Coastal South Africa and the coevolution of the modern human lineage and the Coastal adaptation. Bicho NF, Haws JA, Davis LG, editors, Trekking the shore: Changing coastlines and the antiquity of coastal settlement, New York, NY: Springer New York, Interdisciplinary Contributions to Archaeology. pp. 421–440. https://doi.org/10.1007/978–1-4419–8219–3
  100. 100. Mcbrearty S, Brooks AS (2000) The revolution that wasn’t: a new interpretation of the origin of modern human behavior. Journal of Human Evolution 39: 453–563.
  101. 101. Clark JD, Beyene Y, WoldeGabriel G, Hart WK, Renne PR, et al. (2003) Stratigraphic, chronological and behavioural contexts of Pleistocene Homo sapiens from Middle Awash, Ethiopia. Nature 423: 747–752.
  102. 102. McDougall I, Brown FH, Fleagle JG (2005) Stratigraphic placement and age of modern humans from Kibish, Ethiopia. Nature 433: 733–736.
  103. 103. Smith TM, Ta_oreau P, Reid DJ, Grün R, Eggins S, et al. (2007) Earliest evidence of modern human life history in North African early Homo sapiens. Proc Natl Acad Sci USA 104: 6128–6133.
  104. 104. White TD, Asfaw B, DeGusta D, Gilbert H, Richards GD, et al. (2003) Pleistocene Homo sapiens from Middle Awash, Ethiopia. Nature 423: 742–747.
  105. 105. Behar DM, Villems R, Soodyall H, Blue-Smith J, Pereira L, et al. (2008) The dawn of human matrilineal diversity. Journal of Human Genetics 82: 1130–1140.
  106. 106. Fagundes NJR, Ray N, Beaumont M, Neuenschwander S, Salzano FM, et al. (2007) Statistical evaluation of alternative models of human evolution. Proc Natl Acad Sci USA 104: 17614–17619.
  107. 107. Ingman M, Kaessmann H, Pääbo S, Gyllensten U (2000) Mitochondrial genome variation and the origin of modern humans. Nature 408: 708–713.
  108. 108. Gonder MK, Mortensen HM, Reed FA, De Sousa A, Tishkoff SA (2007) Whole-mtDNA genome sequence analysis of ancient African lineages. Molecular Biology and Evolution 24: 757–768.