MBE Advance Access originally published online on March 21, 2008
Molecular Biology and Evolution 2008 25(6):1209-1218; doi:10.1093/molbev/msn068
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Research Articles |
Climate Change and Postglacial Human Dispersals in Southeast Asia
,1
* Institute of Integrative and Comparative Biology, Faculty of Biological Sciences, University of Leeds, Leeds, UK
Transfusion Medicine and Anthropology Laboratory, Mackay Memorial Hospital, Tamsui, Taiwan
Leverhulme Centre for Human Evolutionary Studies, University of Cambridge, Cambridge, UK
Department of Health, Tainan Hospital, Tainan City, Taiwan
|| The Estonian Biocentre, Tartu, Estonia
¶ Department of Forensic and Chemical Sciences, Anglia Ruskin University, Cambridge, UK
# Department of Statistics, University of Glasgow, Glasgow, UK
** School of Archaeology and Anthropology, The Australian National University, Canberra, Australia
School of Anthropology, University of Oxford, Oxford, UK
E-mail: m.b.richards{at}leeds.ac.uk
 |
Abstract |
|---|
 TOP Abstract IntroductionMaterial and MethodsResultsDiscussionSupplementary MaterialAcknowledgementsReferences |
|---|
Modern humans have been living in Island Southeast Asia (ISEA) for at least 50,000 years. Largely because of the influence of linguistic studies, however, which have a shallow time depth, the attention of archaeologists and geneticists has usually been focused on the last 6,000 years—in particular, on a proposed Neolithic dispersal from China and Taiwan. Here we use complete mitochondrial DNA (mtDNA) genome sequencing to spotlight some earlier processes that clearly had a major role in the demographic history of the region but have hitherto been unrecognized. We show that haplogroup E, an important component of mtDNA diversity in the region, evolved in situ over the last 35,000 years and expanded dramatically throughout ISEA around the beginning of the Holocene, at the time when the ancient continent of Sundaland was being broken up into the present-day archipelago by rising sea levels. It reached Taiwan and Near Oceania more recently, within the last
8,000 years. This suggests that global warming and sea-level rises at the end of the Ice Age, 15,000–7,000 years ago, were the main forces shaping modern human diversity in the region.
Key Words: complete mtDNA genomes • Island Southeast Asia • Neolithic • postglacial • late glacial
|
Introduction |
|---|
|
TOPAbstractIntroductionMaterial and MethodsResultsDiscussionSupplementary MaterialAcknowledgementsReferences |
|---|
Between the 1970s and the 1990s, archaeologists and linguists developed a model for the settlement of Island Southeast Asia (ISEA) that became widely accepted. The "out of Taiwan" model suggested that settlement occurred in 2 waves of dispersal. The first of these, a late Pleistocene arrival of so-called "Australo-Melanesian" people
50,000 years ago (Barker et al. 2007
), led to the initial peopling of both Sahul (Hudjashov et al. 2007) and the ancient continent of Sundaland, which comprised Indo-China, Malaysia, and western Indonesia before the sea-level rises that began 19,000 years ago (fig. 1). These people were then replaced or assimilated in the mid-Holocene by the "Austronesian expansion"—a maritime dispersal of rice agriculturalists from southern China via Taiwan into ISEA between 5,500 and 4,000 years ago, who brought with them the Austronesian languages now spoken throughout the region (Bellwood 1997; Blust 1999; Diamond and Bellwood 2003). This view implies that these were the 2 most important demographic events shaping modern human diversity in the region, with the latter generally accorded far greater attention than the former.
|
The benchmark status of a mid-Holocene out of Taiwan dispersal has been questioned by genetic studies, especially those focusing on variation in the nonrecombining mitochondrial DNA (mtDNA) and Y chromosome (Richards et al. 1998
; Kayser et al. 2000, 2001; Su et al. 2000; Capelli et al. 2001; Trejaut et al. 2005). In particular, from a detailed consideration of more than 1,000 mtDNA control region sequences of samples from throughout ISEA and Taiwan, we have recently shown that, at most, only 20% of modern mtDNAs in ISEA were introduced at the time of the Neolithic (Hill et al. 2007).
These results could be compatible with a version of out of Taiwan that allows for a greater level of assimilation of indigenous Pleistocene lineages during the Austronesian expansion than had previously been anticipated (Hurles et al. 2002; Paz 2002). However, some scholars have suggested a more radical view in which the key demographic processes involved not (or not only) the spread of farming but late glacial and postglacial dispersals as a result of climate change and sea-level rises (Meacham 1984–1985; Oppenheimer 1998; Solheim 2006). We have recently shown that the majority of the mtDNAs in the region indeed coalesce to founder ages between 5,000 and 25,000 years ago, suggesting independent genetic evidence for such a view (Hill et al. 2007). Climate change at the end of the last Ice Age has been shown to have had an important effect on human genetic diversity in western Europe, where increased glaciation drove hunter-gatherer groups into southern refuge areas (Torroni et al. 1998; Torroni, Bandelt, et al. 2001; Achilli et al. 2004; Gamble et al. 2005; Pereira et al. 2005). The potential impact in Southeast Asia was at least as profound: around half of the land area of the continent of Sundaland was lost to the sea between 15,000 and 7,000 years ago, with a concomitant doubling of the length of coastline as the resulting archipelago was formed (Bird et al. 2005; Sathiamurthy and Voris 2006). Moreover, it seems likely that the sea level rose in a series of rapid steps that could have had a catastrophic impact on human populations (Blanchon and Shaw 1995; Pelejero et al. 1999; Hanebuth et al. 2000).
Here we use the sequence variation in complete mtDNA genomes to investigate the possibility that a signal of late glacial and postglacial dispersals exists in ISEA (Richards and Macaulay 2001). Although many mtDNA clusters coalesce in the time window of interest (Hill et al. 2007), we have focused particularly on mtDNA haplogroup E. This haplogroup accounts for 15% of mtDNA lineages in ISEA and Taiwan and is almost entirely restricted to this region; it therefore presents a particularly clear and highly informative distribution (Hill et al. 2007). Although rarely found elsewhere, it nests within a larger haplogroup, M9, that includes a sister clade to E, M9a'b, much of which is restricted to Japan and the East Asian mainland.
|
Material and Methods |
|---|
|
TOPAbstractIntroductionMaterial and MethodsResultsDiscussionSupplementary MaterialAcknowledgementsReferences |
|---|
Subjects
We completely sequenced 50 haplogroup E mtDNAs: 25 from 9 locations throughout Indonesia/East Malaysia, 9 from 4 locations in the Philippines, 14 from 8 aboriginal groups in Taiwan (including several who no longer speak Austronesian languages), and 2 from the north coast of Papua New Guinea (PNG). Samples were selected on the basis of hypervariable segment I (HVS-I) control region variation and the results of restriction typing for diagnostic sites. Haplogroup E includes an HVS-I motif with transitions at positions 16223, 16362, and 16390 with respect to the reference sequence (Andrews et al. 1999
); in addition, we tested for the diagnostic loss of restriction site –7598 HhaI in the coding region. We also sequenced 3 "pre-M9a" lineages identified on the basis of control region polymorphisms at nucleotide positions (nps) 16234 and 16362 shared with M9a. The human DNA samples used in the UK were anonymized unlinked samples in accordance with UK ethical guidelines, and the work was approved by the University of Leeds, Faculty of Biological Sciences Ethics Committee. The Taiwan samples were collected with the informed consent of the subjects and the work approved by the Human Experiment Committee of the Mackay Memorial Hospital of Taipei.
We included one published E sequence from the Philippines in the analysis (Ingman and Gyllensten 2003) and a further published coding region sequence for the age estimates (Herrnstadt et al. 2002); we also included a further 13 published M9a sequences (Ingman et al. 2000; Tanaka et al. 2004), an unclassified M9* sequence (diverging independently from the root of M9) previously identified in Malaysia (Macaulay et al. 2005), and an M9b sequence (Kong et al. 2006).
Sequence Analysis
We performed the complete sequencing as previously described (Torroni, Rengo, et al. 2001), either at the University of Leeds using an ABI 16-capillary 3130XL DNA Analyzer or in Taipei using an ABI 48-capillary 3730 DNA Analyzer.
Phylogenetic and Statistical Analyses
For the phylogenetic analysis, we used (unweighted) the reduced median (RM) network algorithm of the package Network 4.1 (Bandelt et al. 1995), with all characters included with the exception of size differences in the polyC tracts at positions 309 and 315. A previously published M9* sequence from Malaysia (Macaulay et al. 2005) served as a suitable outgroup, and the RM network equates to the single most parsimonious tree in this case with the exception of a handful of characters. These are 3 fast-evolving control region characters, at positions 153, 16390, and 16519, which generate reticulations in the RM network. In the E2a clade, the positions 16185 and 5460 also generated a reticulation. The mutation at position 16185 appears independently far less frequently than that at 5460 in a database of 2,000 complete sequences, so we assumed a single event at this site in E2a.
Time depths were estimated from the more slowly evolving coding region of the molecule between nps 577 and 16023. We obtained the time estimates by converting the averaged distance (
) of the haplotypes of a clade from the respective root haplotype, accompanied by a heuristic estimate of standard error, 
, following Saillard et al. (2000). We also obtained maximum likelihood (ML) estimates from the coding region where possible (i.e., for clades defined by coding region mutations), using PAML 3.13 (Yang 1997), assuming the HKY85 mutation model with gamma-distributed rates (approximated by a discrete distribution with 32 categories). We converted in the coding region and mutational distance in ML to time using a mutation rate estimate of 1.26 x 10–8 base substitutions per nucleotide per year, corresponding to one substitution per 5,140 years in the coding region (Mishmar et al. 2003).
We used HVS-I information to augment the age estimates in some cases, using a transition rate of 1 in 20,180 years for the region between nps 16090 and 16365 (Forster et al. 1996). Table 1 shows a comparison of overall time depths of nodes in the gene tree estimated from both coding-region sequences and HVS-I data, indicating that the age estimates are comparable. Additionally, we compared the coding region ML age estimates with the estimates (table 2). The ML age estimates did not differ substantially from the ones calculated with .
|
|
The reconstructed tree was scaled against the
coding region point age estimates of the major nodes. The estimates for 2 of the clades, E1a1 and E1a1a, were unfeasible because the point estimate for the ancestor node E1a1 (8,900 ± 2,150) was lower than the derived node E1a1a (10,300 ± 2,500). Both E1a and E1a1 present HVS-I motifs corresponding to the root type of E (with transitions at positions 16223, 16362, and 16390). Observed HVS-I variation in these lineages was completely restricted, albeit rarely, to ISEA. The samples presented in the complete sequence tree for the node E1a1 are from Taiwan, most likely the sink (or descendant) population, and almost lack HVS-I variation. This founder effect in Taiwan has therefore probably reduced the overall age estimate of the node E1a1. So the more frequent and star-like E1a1a clade should present a more reliable age estimate than its ancestor node and is therefore the one scaled in the tree.
Haplotype diversity was calculated as 1 – 
ixi2, where xi is the relative frequency of the ith haplotype in the sample (Torroni, Bandelt, et al. 2001).
Phylogeographic Distributions
We used the database of published mtDNA HVS-I sequences and our own unpublished data in order to draw more robust conclusions about frequency distributions than is possible using the numerically more limited complete genome data. Haplogroup E includes an HVS-I motif with transitions at positions 16223, 16362, and 16390. E1a1a has an additional transition at position 16291, E1b is distinguished in HVS-I by a transition at position 16261, E2 by a transition at position 16051, and E2b by a transition at position 16185. M9a carries transitions at positions 16234, 16316, and 16362; 16234 and 16362 occur on the "pre-M9a" branch that diverges from the ancestor of M9 50,000 years ago. Given these motifs, and the size of the HVS-I database, we can have considerable confidence that our frequency estimates are broadly accurate. M9b is more difficult because of its scarcity: it has only been identified twice: once as a complete sequence from western China (Kong et al. 2006) and once as the HVS-I sequence (motif: 16051, 16209, 16223, and 16362) in a Yao speaker from southeast China (Wen, Li, Gao, et al. 2004).
We displayed frequency distributions using the Kriging algorithm of Surfer 8 using the geographic points indicated in figure 1. The haplogroup E distributions were based on a total of 1,704 samples from throughout ISEA (Hill et al. 2006, 2007; authors’ unpublished data), 184 from West Papua (authors’ unpublished data), and 962 from Taiwan (Tajima et al. 2003; Trejaut et al. 2005; authors’ unpublished data). The frequencies displayed in Taiwan correspond solely to Austronesian-speaking aboriginals. The distribution of the haplogroup E sequences is shown in supplementary table S1 (Supplementary Material online). All the haplogroup E sequences identified on the basis of their HVS-I motif in ISEA were further checked by sequencing for position 10834 (defining E1). Samples with the control region transition at position 16051 were checked by sequencing for the position 8730 (defining E2a) and 8 Indonesian and 4 Taiwanese from the E1 branch that lacked either one of the informative control region variants (at nps 16291 or 16261) were checked for np 6620 position (defining E1a). Published data from China, Thailand, Malaysia, West Papua, and Australia were also included in the analyses (van Holst Pellekaan et al. 1998; Fucharoen et al. 2001; Kivisild et al. 2002; Tommaseo-Ponzetta et al. 2002; Yao et al. 2002a; Yao et al. 2002b; Yao and Zhang 2002; Wen, Li, Gao, et al. 2004; Wen, Li, Lu, et al. 2004; Wen, Xie, et al. 2004). There are 3 potential haplogroup E sequences in a Chinese data set of 2,618 sequences (i.e., 0.1%), 2 in southern China (Kivisild et al. 2002; Yao et al. 2002b) and 1 in western China (Yao et al. 2000) which may be "accidentals" from ISEA/Taiwan. The East Asian data sets were additionally used to build a frequency distribution map for M9a.
|
Results |
|---|
|
TOPAbstractIntroductionMaterial and MethodsResultsDiscussionSupplementary MaterialAcknowledgementsReferences |
|---|
Haplogroup M9
The phylogenetic tree of M9 with reconstructed character evolution is shown in figure 2. Although M9a is restricted to mainland Asia and Japan, the deepest branches in the clade encompassing it (pre-M9a) are found in Indo-China, China, and Taiwan and the next deepest in ISEA. M9b is the result of the earliest split from pre-M9a within M9 but its phylogeographic usefulness is minimal due to its apparent extreme rarity, with only 2 instances to date in southern and western China (Wen, Li, Gao, et al. 2004
; Kong et al. 2006). Along with the finding of an unclassified M9* lineage in the aboriginal inhabitants of Malaysia (Macaulay et al. 2005) and the insular distribution of haplogroup E itself this distribution suggests that M9 as a whole is most likely to have had a Southeast Asian origin, 50,000 years ago. M9a itself appears to have spread north into the East Asian mainland 15,000 years ago, after the last glacial maximum (LGM). This is supported by the phylogeography of HVS-I variation: to date almost all the deeper branching lineages within the published pre-M9a HVS-I data are found in southern China and Southeast Asia and all northeast Asian lineages belong to M9a, itself as characterized by a transition at np 16316.
|
It is worth noting that M9 is only defined by a single nucleotide change, at position 4491, and its status as a clade may therefore be open to some doubt. A database of
2,000 published mtDNA coding region sequences suggests that position 4491 evolves slowly in the global mtDNA tree, where we just detected 6 occurrences in more than 10,000 mutations. Nevertheless, it is possible that the Malaysian outgroup labeled M9* may in fact be the result of an independent mutation at position 4491. The link between haplogroups E and pre-M9a is slightly stronger as it also includes a mutation at the control region site 16362, although this position is fast evolving (Bandelt et al. 2006). The inferred ancestry of pre-M9a in Southeast Asia provides an additional phylogeographic argument for a common ancestry. None of our principal arguments to follow rely on the inferred reality of M9 as a clade, but the phylogenetic link would strengthen the case for their common ancestry in the original Sunda population.
Haplogroup E
Haplogroup E itself has a likely time depth of more than 30,000 years but falls into 2 major subclades, E1 and E2, which are 17,000 and 9,500 years old, respectively (fig. 2). They are each subdivided into 2 further principal subclades, E1a, E1b, E2a, and E2b, which range from 4,500 to 11,500 years old. These subclades have a very distinctive geographic distribution, which is highly informative about the demographic history of the region. Although all 4 subclades are found in ISEA, so that ISEA lineages occur throughout the tree, only 2 of the 4 (E1a [fig. 3a] and E2b [fig. 3d]) are found in Taiwan. The remaining subclades, E1b (fig. 3b) and E2a (fig. 3c), are both largely restricted to ISEA, albeit extending occasionally to New Guinea and the Malay Peninsula, suggesting that both arose in ISEA and dispersed fairly recently east and west. Of the 2 subclades shared with Taiwan, E2b, which is only 4,300 (±2,300) years old, falls into 2 branches, both present in ISEA, with the Taiwanese lineages forming a clear, shallower founder subclade. In other words, the diversity of E2b seen in Taiwan is a subset of the clade found in ISEA. This implies an origin in ISEA and subsequent migration to Taiwan.
|
Therefore, it is clear that 3 of the 4 subclades of haplogroup E arose in ISEA. The fourth and oldest subclade, E1a, is 11,700 (±3,450) years old, and Taiwanese and ISEA lineages are interleaved within it so that a direction of dispersal is less obvious. However, the overall diversity translates into an age estimate that is somewhat higher in ISEA (13,100 ± 3,700 years) than in Taiwan (10,750 ± 4,200 years). Furthermore, 1 of the 2 major subclades of E1a, E1a2 (at 9,400 ± 2,850 years), is only present in ISEA, suggesting an origin in this region. It is also worth noting that there is very low diversity within the control region in the E1a* sequences (i.e., excluding E1a1a) in Taiwan (haplotype diversity = 0.131 and an estimate of 1,400 ± 1,000 years using
), suggesting a recent founder effect. The derived and most common clade E1a1a presents control region diversity in Taiwan that approximates the diversity seen in what would have been Sundaland, using both haplotype diversity (table 3) and statistics. It should be noted that the value in Taiwan is elevated by the presence of 2 high–frequency derived haplotypes that, together with the root type, account for more than 80% of the samples of this clade in Taiwan (supplementary table S1, Supplementary Material online), again strongly suggesting a founder effect.
|
This problem is also amenable to dissection by an exploratory founder analysis that assumes either Taiwan or ISEA as the source in turn. This procedure estimates the time of arrival of lineages in the putative sink population by subtracting the diversity already present in the putative source at the time of the dispersal (Richards et al. 2000
). Applying such a reciprocal founder analysis, the founder age for E1a in Taiwan assuming an ISEA source would be 6,550 (±1,850) years, whereas a dispersal in the opposite direction would yield an age of 10,000 (±2,400) years in ISEA, reflecting the higher proportion of private variants in the latter compared with Taiwan. This again suggests that ISEA is the more likely source. We further assessed the lower diversity in Taiwan by simply calculating age estimates for each of the major clades or paragroups in the tree separately for Taiwan and ISEA (the latter including PNG) (table 4). The only case in which the age estimate is higher in Taiwan is haplogroup E2b. However, in that case, the diversity in Taiwan is restricted to the E2b1 node, corresponding clearly to a shallow founder clade. Moreover, the HVS-I haplotype diversity of E2b in Taiwan is extremely low (table 3).
|
We then explored the diversity statistics from the more numerous control region data in further detail, in order to further test this conclusion. We used HVS-I haplotype diversity within each subclade (and the paragroup E1a*, which includes the phylogenetically deeper parts of E1a) to assess variation between regions and the possible impact of recent founder effects. Table 3 shows the diversity of each subclade in both ISEA (divided into the former Sunda area, the Philippines, and the islands of Wallacea to the east, which were never connected to the mainland) and Taiwan. This shows that the higher diversity for each of the clades shared between Taiwan and ISEA is always found in one of the regions of ISEA and never in Taiwan.
Spatial frequency distributions suggest that the deeper branches of E1a and E2a may have diversified in Wallacea, where these haplogroups now occur at the highest frequencies (fig. 3a and c). However, caution is warranted as the sample sizes of these haplogroups are low. The 2 subclades with higher sample sizes (E1a1a and E1b) have a higher diversity in Sundaland than in Wallacea or the Philippines. This pattern is particularly striking in the most common clade, E1a1a. Furthermore, all subclades of haplogroup E are present in Borneo but not in Wallacea.
This suggests that the most likely location for the "homeland" of haplogroup E as a whole—that is to say, the place in which the defining mutations arose—is somewhere in the vicinity of northeast Sundaland. However, the diversity and frequency patterns suggest that it is quite possible that the subclades E2a and E1a* may have arisen in northwest Wallacea. In this case, the area of diversification of haplogroup E lineages would be distributed around the coastlines of the Sulawesi and Sulu Seas, as might be expected for lineages carried by a mobile, maritime-oriented population.
Wherever the precise origin of haplogroup E may lie, spatial frequency distribution and diversity strongly suggest that the haplogroup arose in ISEA and some of its subclades spread subsequently to Taiwan. Furthermore, we can infer a timescale for these processes using the coding region molecular clock. Early subclades are largely restricted to ISEA and date to the early Holocene. These are E1a1a1 (age 7,700 ± 3,150 years), E1a2 (age 9,400 ± 2,850 years), E1b (age 6,850 ± 2,150 years), and E2a (age 6,700 ± 1,800 years). This implies that the diversity of haplogroup E has been accumulating within ISEA for a minimum of 7,000–10,000 years.
Applying a founder analysis to the subclades found in Taiwan, assuming ISEA as the source, yielded arrival time estimates of 6,550 ± 1,850 years ago for E1a and 5,150 ± 2,950 years ago for E2b. As these values were quite similar, we also applied a founder analysis to the entire haplogroup, combining data from both subclades, effectively assuming a single dispersal. This overall founder age for Taiwan came to 6,250 ± 1,600 years. These results indicate that haplogroup E evolved in situ in ISEA and was involved in range expansions from 12,000 years ago, reaching Taiwan roughly 4,000–8,000 years ago.
We can also consider the control region database from a phylogeographic point of view to add additional detail to this picture. Control region sequence patterns suggest range expansions across eastern Indonesia in the mid-Holocene. The derived E1a1a lineages (which are distinguished by a diagnostic transition at np 16291 in the control region) present 2 clear zones of diversity across ISEA: western Indonesia (13,450 ± 5,600 years) and Wallacea (5,000 ± 1,700 years). The diversity seen in the Philippines amounts to 6,600 (±3,800) years. Therefore, it seems likely that much of the Philippines and Wallacea were settled at least in part during the mid-Holocene following diversification in the Sulu Sea/Sulawesi Sea region.
Several haplogroup E (mainly E1b) lineages in Near Oceania suggest subsequent gene flow further east (Friedlaender et al. 2007), perhaps along a sphere of interaction or "voyaging corridor" between Southeast Asia and the Solomon Islands that may have been established by 6,000 years ago (Terrell and Welsch 1997) or earlier (O'Connor and Veth 2005). A few putative haplogroup E lineages observed in the control region database for China also suggest possible very minor (and probably very recent) gene flow toward the mainland.
|
Discussion |
|---|
|
TOPAbstractIntroductionMaterial and MethodsResultsDiscussionSupplementary MaterialAcknowledgementsReferences |
|---|
Haplogroup E allows us to illuminate demographic processes in ISEA prior to the proposed Austronesian expansion from China/Taiwan. This haplogroup probably evolved within the descendants of the first settlers of Sundaland, who arrived carrying haplogroup M >50,000 years ago (Macaulay et al. 2005
). Haplogroup E evolved, most likely from haplogroup M9, 35,000 years ago and was caught up in a dramatic series of dispersals and expansions that began in eastern Sundaland/northwest Wallacea from 12,000 years ago. This may well have been part of the same set of processes that also led to the dispersal of its probable sister haplogroup, M9a, north along the Pacific rim toward Korea and Japan from 15,000 years ago (fig. 3e). This pattern confirms the suggestion of postglacial expansion and dispersal proposed for haplogroup E on the basis of HVS-I sequence variation by Hill et al. (2007). Indeed, the HVS-I data for ISEA suggested expansions and dispersals carrying other mtDNA haplogroups in addition to haplogroup E at the same time (Hill et al. 2007). These will need to be similarly tested at the level of complete genome sequencing in order to resolve the signals clearly. The present data clearly demonstrate that these demographic processes took place and that they had a significant impact on the mtDNA pool of modern Southeast Asians.
The most plausible explanation for the spread of these genetic signatures is the impact on coastal-dwelling populations of the rapid global warming and sea-level rises that led to the inundation of the Sunda shelf by meltwater at the end of the last Ice Age (Oppenheimer 1998; Lin et al. 2005). Although sea levels began to rise gradually after the end of the LGM 19,000 years ago, it is likely that there were 3 major episodes of accelerated sea-level rise and flooding, probably due to ice sheet collapse and sometimes referred to as Catastrophic Rise Events 1–3, at 14,500, 11,500, and probably also 7,500 years ago (Blanchon and Shaw 1995; Oppenheimer 1998; Pelejero et al. 1999; Hanebuth et al. 2000; Voris 2000; Lambeck and Chappell 2001; Bird et al. 2005). It has been suggested that these episodes triggered major displacements of human groups living on the Sunda coastline and had an important role in shaping subsequent life in the region, in particular its maritime orientation and the development of sailing technology (Oppenheimer 1998; Solheim 2006).
With the opening up of approximately twice as much coastline as there had previously been, the coastal populations of Sundaland (and perhaps also northwestern Wallacea) would have been well placed to expand into this rapidly growing habitat. Conversely, the hinterland-based populations would have faced severe land pressures from both the drowning of their ancestral lands and the encroachment of rainforest over the savannas and monsoon forests to which they had adapted (Bird et al. 2005). A scenario of overall population growth, combined with intense competition at optimal resource locations on the growing coastline, would have provided the same motivations for dispersal that are conventionally associated with agriculturalist expansion. The fact that haplogroup E presents only 3 branches beyond 12,000 years ago, and only 2 long branches from 17,000 to 33,000 years ago, may however provide a glimpse of the scale of the preceding catastrophe: the heavy lineage losses may perhaps suggest a population decline before the more recent reexpansions.
These results open up a barely explored perspective on the population history of Homo sapiens in Southeast Asia. They provide a window onto major, environmentally related demographic changes that took place before the proposed Neolithic/Austronesian dispersal that is usually the main point of reference for the prehistory of this region. Moreover, archaeological evidence provides independent support for a large-scale dispersal across ISEA (or central/eastern Sundaland and Wallacea, as the region then was) during the 15,000 years after the LGM.
Prior to the LGM, there is no very clear distinction between the stone artifact assemblages of mainland Southeast Asia (MSEA—then, northwestern Sundaland) and ISEA. However, following the LGM, 2 geographically discrete traditions or "technocomplexes" emerged. The Hoabinhian appeared across MSEA (including the Malay Peninsula) and the northern two-thirds of Sumatra. It was based on pebble tools flaked along their edges, usually on one face but occasionally bifacially (Solheim 1994; Bellwood 1997; Bulbeck 2003; Forestier 2007). In contrast, the "flake–blade technocomplex" was based on flakes detached from rotated multiplatform cores. This appears to have emerged as a distinctive stone tool technology approximately 25, 000–30,000 years ago restricted to the islands and coastlines of the Sulu Sea region. After it had spread to northern Borneo 18,000 years ago, it evidently dispersed throughout ISEA following the LGM (Bellwood 1997; Tanudirjo 2001; Mahirta 2003; Arifin 2004; Veth et al. 2005; Mijares 2006; Forestier 2007; Mijares 2007).
The expansion of the Hoabinhian from South China/Vietnam has already been related to the arrival of the R9b and N9a mtDNA haplogroups into the Malay Peninsula (Hill et al. 2006). Here, we now find a very close relationship between the geographic extent of post-LGM flake–blade industries and the haplogroup E lineages, again suggesting a close match between archaeological and genetic reconstructions.
Although Pleistocene assemblages assigned to the flake–blade technocomplex lack distinct stone tool types, post-LGM dispersal of Hoabinhian populations into ISEA, or flake–blade technocomplex populations into MSEA/Sumatra, would be expected to leave the mark of anomalous assemblages within their regional context—which are yet to be recorded archaeologically. The archaeology of Taiwan, however, does reveal the incursion of a distinct stone tool technology, which is particularly similar to North Luzon Holocene industries (cf., Mijares 2007). Although most Taiwan flaked stone industries of the last 15,000 years are assigned to the Changpanian culture, which is related to the Hoabinhian (Bellwood 1997; Rolett et al. 2002), Zhang (2000) describes the appearance of a flake–blade industry, based on flakes struck from multiplatform cores, at 4 coastal sites in east Taiwan, all dated to 5,000–6,000 years ago. This would fit well with our estimated arrival from the south of E1a lineages in Taiwan.
Both the genetic and archaeological evidence therefore suggest that a maritime-oriented culture was in operation around the coastlines of what are now the Sulu and Sulawesi Seas by the LGM, which would have been preadapted to benefit from rising sea levels (Meacham 1984–1985; Oppenheimer 1998; Solheim 2006). The signature of dispersals from this region at the end of the Ice Age is evident both in the distribution of mtDNA haplogroup E lineages and the expansion of the flake–blade technocomplex.
These results only account for 15% of ISEA mtDNA lineages, but preliminary analyses based on HVS-I diversity suggest that many others may share a similar history (Hill et al. 2007). Complete genome analysis, such as we have carried out here, will be needed to test this further. It also remains to be seen whether a similar signal will be detected in the Y chromosome variation of Southeast Asians; the distribution of the main Y chromosome haplogroups in ISEA and Taiwan (Capelli et al. 2001) would be consistent with this scenario, but greater phylogenetic resolution is necessary to identify the directions of dispersal. A recent study of autosomal short tandem repeat variation in Taiwan and the western Pacific (Friedlaender et al. 2008) indicates the possibility of an ISEA component in the diversity of remote Pacific islanders that needs to be explored further.
The discovery of a signature of climate change in the human genome of Southeast Asians echoes similar findings for mtDNA haplogroups V and H in Europe. These haplogroups were carried from a glacial refuge in southwest Europe into western, central, and northern Europe following the retreat of the glaciers from 15,000 years ago (Richards et al. 2000; Torroni, Bandelt, et al. 2001; Gamble et al. 2004; Pereira et al. 2005). Our results therefore reemphasize the critical role of climate change—in particular global warming after the LGM—in shaping the evolution of the modern human gene pool.
|
Supplementary Material |
|---|
|
TOPAbstractIntroductionMaterial and MethodsResultsDiscussionSupplementary MaterialAcknowledgementsReferences |
|---|
Supplementary table S1 is available at Molecular Biology and Evolution online (http://www.mbe.oxfordjournals.org/). The sequence data for this paper appear in GenBank under accession numbers: EF093535 [GenBank] –EF093558 [GenBank] and EF185793 [GenBank] –EF185816.
|
Acknowledgements |
|---|
|
TOPAbstractIntroductionMaterial and MethodsResultsDiscussionSupplementary MaterialAcknowledgementsReferences |
|---|
We thank the Bradshaw Foundation, the British Academy, and the National Science Council of Taiwan (grant 95-2627-H-195-001) for financial support. P.S. was supported by a Marie Curie Early Stage Training Grant. We also thank John Clegg, A. S. M Sofro, and the Taipei County Labor Department for the samples, Harold Voris and Clara Simpson for the map used in figure 1, Toomas Kivisild for critical comments, and Alessandro Achilli for advice on the use of Surfer.
|
Footnotes |
|---|
1 These authors contributed equally to this work.
Connie Mulligan, Associate Editor
|
References |
|---|
|
TOPAbstractIntroductionMaterial and MethodsResultsDiscussionSupplementary MaterialAcknowledgementsReferences |
|---|
Achilli A, Rengo C, Magri C. (21 co-authors). The molecular dissection of mtDNA haplogroup H confirms that the Franco-Cantabrian glacial refuge was a major source for the European gene pool. Am J Hum Genet (2004) 75:910–918.[CrossRef][Web of Science][Medline]
Andrews RM, Kubacka I, Chinnery PF, Lightowlers R, Turnbull D, Howell N. Reanalysis and revision of the Cambridge reference sequence for human mitochondrial DNA. Nat Genet (1999) 23:147.[CrossRef][Web of Science][Medline]
Arifin K. Early human occupation of the East Kalimantan rainforest (the upper Birang river region, Berau) (2004) Canberra (Australia): The Australian National University. 277–278.
Bandelt H-J, Forster P, Sykes BC, Richards MB. Mitochondrial portraits of human populations using median networks. Genetics (1995) 141:743–753.[Abstract]
Bandelt HJ, Kong QP, Yao YG, Richards M, Macaulay V. Estimation of mutation rates and coalescence times: some caveats. In: Mitochondrial DNA and the evolution of Homo sapiens—Bandelt HJ, Macaulay V, Richards M, eds. (2006) Berlin (Germany): Springer. 47–90.
Barker G, Barton H, Bird M. (27 co-authors). The ‘human revolution’ in lowland tropical Southeast Asia: the antiquity and behavior of anatomically modern humans at Niah Cave (Sarawak, Borneo). J Hum Evol (2007) 52:243–261.[CrossRef][Medline]
Bellwood P. Prehistory of the Indo-Malaysian Archipelago (1997) 2nd ed. Honolulu (Hawaii): University of Hawaii Press.
Bird MI, Taylor D, Hunt C. Palaeoenvironments of insular Southeast Asia during the last glacial period: a savanna corridor in Sundaland? Q Sci Rev (2005) 24:2228–2242.[CrossRef]
Blanchon P, Shaw J. Reef drowning during the last glaciation: evidence for catastrophic sea-level rise and ice-sheet collapse. Geology (1995) 23:4–8.
Blust R. Subgrouping, circularity and extinction: some issues in comparative Austronesian linguistics. In: Symposium Series of the Institute of Linguistics (1999) 1. Academica Sinica. 31–94.
Bulbeck D. Hunter-gatherer occupation of the Malay Peninsula from the Ice Age to the Iron Age. In: Under the canopy: the archaeology of tropical rain forests—Mercader J, ed. (2003) New Brunswick (NJ): Rutgers University Press. 119–160.
Capelli C, Wilson JF, Richards M, Stumpf MPH, Gratrix F, Oppenheimer S, Underhill P, Pascali VL, Ko T-M, Goldstein DB. A predominantly indigenous paternal heritage for the Austronesian-speaking peoples of insular Southeast Asia and Oceania. Am J Hum Genet (2001) 68:432–443.[CrossRef][Web of Science][Medline]
Diamond J, Bellwood P. Farmers and their languages: the first expansions. Science (2003) 300:597–603.
Forestier H. Ribuan Gunung, Ribuan Alat Batu: Prasejarah Song Keplek, Gunung Sewu, Jawa Timur (2007) Jakarta (Indonesia): Kepustakaan Populer Gramedia.
Forster P, Harding R, Torroni A, Bandelt H-J. Origin and evolution of Native American mtDNA variation: a reappraisal. Am J Hum Genet (1996) 59:935–945.[Web of Science][Medline]
Friedlaender JS, Friedlaender FR, Hodgson JA, Stolz M, Koki G, Horvat G, Zhadanov S, Schurr T, Merriwether DA. Melanesian mtDNA complexity. PLoS One (2007) 2:e248.[CrossRef]
Friedlaender JS, Friedlaender FR, Reed FA. (12 co-authors). The genetic structure of Pacific Islanders. PLoS Genet (2008) 4:e19.[CrossRef][Medline]
Fucharoen G, Fucharoen S, Horai S. Mitochondrial DNA polymorphisms in Thailand. J Hum Genet (2001) 46:115–125.[CrossRef][Web of Science][Medline]
Gamble C, Davies W, Pettitt P, Hazelwood L, Richards M. The archaeological and genetic foundations of the European population during the Lateglacial: implications for ‘agricultural thinking’. Cambridge Arch J (2005) 15:193–223.[CrossRef]
Gamble C, Davies W, Pettitt P, Richards M. Climate change and evolving human diversity in Europe during the last glacial. Philos Trans R Soc Lond B Biol Sci (2004) 359:243–254.[CrossRef][Web of Science]
Hanebuth T, Stattegger K, Grootes PM. Rapid flooding of the Sunda Shelf: a late-glacial sea-level record. Science (2000) 288:1033–1035.
Herrnstadt C, Elson JL, Fahy E. (11 co-authors). Reduced-median-network analysis of complete mitochondrial DNA coding-region sequences from the major African, Asian, and European haplogroups. Am J Hum Genet (2002) 70:1152–1171. (erratum 1171:1448–1449).[CrossRef][Web of Science][Medline]
Hill C, Soares P, Mormina M. (12 co-authors). Phylogeography and ethnogenesis of aboriginal Southeast Asians. Mol Biol Evol (2006) 23:2480–2491.
Hill C, Soares P, Mormina M. (11 co-authors). A mitochondrial stratigraphy for Island Southeast Asia. Am J Hum Genet (2007) 80:29–43.[CrossRef][Web of Science][Medline]
Hudjashov G, Kivisild T, Underhill PA. (12 co-authors). Revealing the prehistoric settlement of Australia by Y chromosome and mtDNA analysis. Proc Natl Acad Sci USA (2007) 104:8726–8730.
Hurles ME, Nicholson J, Bosch E, Renfrew C, Sykes BC, Jobling MA. Y chromosomal evidence for the origins of Oceanic-speaking peoples. Genetics (2002) 160:289–303.
Ingman M, Gyllensten U. Mitochondrial genome variation and evolutionary history of Australian and New Guinean aborigines. Genome Res (2003) 13:1600–1606.
Ingman M, Kaessmann H, Pääbo S, Gyllensten U. Mitochondrial genome variation and the origin of modern humans. Nature (2000) 408:708–713.[CrossRef][Medline]
Kayser M, Brauer S, Weiss G, Schiefenhövel W, Underhill PA, Stoneking M. Independent histories of human Y chromosomes from Melanesia and Australia. Am J Hum Genet (2001) 68:173–190.[CrossRef][Web of Science][Medline]
Kayser M, Brauer S, Weiss G, Underhill PA, Roewer L, Schiefenhövel W, Stoneking M. Melanesian origin of Polynesian Y chromosomes. Curr Biol (2000) 10:1237–1246.[CrossRef][Web of Science][Medline]
Kivisild T, Tolk H-V, Parik J, Wang Y, Papiha SS, Bandelt H-J, Villems R. The emerging limbs and twigs of the East Asian mtDNA tree. Mol Biol Evol (2002) 19:1737–1751. (erratum 1720:1162).
Kong Q-P, Bandelt H-J, Sun C. (12 co-authors). Updating the East Asian mtDNA phylogeny: a prerequisite for the identification of pathogenic mutations. Hum Mol Genet (2006) 15:2076–2086.
Lambeck K, Chappell J. Sea level change through the last glacial cycle. Science (2001) 292:679–686.
Lin M, Chu C-C, Broadberry RE, Yu L-C, Loo JH, Trejaut JA. Genetic diversity of Taiwan's indigenous peoples: possible relationship with insular Southeast Asia. In: The peopling of East Asia—Sagart L, Blench R, Sanchesz-Mazas A, eds. (2005) Abingdon (VA): RoutledgeCurzon. 230–247.
Macaulay V, Hill C, Achilli A. (21 co-authors). Single, rapid coastal settlement of Asia revealed by analysis of complete mitochondrial genomes. Science (2005) 308:1034–1036.
Mahirta. Human Occupation on Rote and Sawu Islands. Nusa Tenggara Timur (2003) Canberra (Australia): The Australian National University. 166–168.
Meacham W. On the improbability of Austronesian origins in South China. Asian Perspect (1984–1985) 26:89–106.
Mijares ASB. Unravelling prehistory: the archaeology of north-eastern Luzon. (2006) Canberra (Australia): The Australian National University. 170–175.
Mijares ASB. Unearthing Prehistory: the Archaeology of Northeastern Luzon, Philippine Islands (2007) Oxford (UK): John and Erica Hedges Ltd.
Mishmar D, Ruiz-Pesini E, Golik P. (13 co-authors). Natural selection shaped regional mtDNA variation in humans. Proc Natl Acad Sci USA (2003) 100:171–176.
O'Connor S, Veth P. Early Holocene shell hooks from Lene Hara Cave, East Timor establish complex fishing technology was in use in Island Southeast Asia 5000 years before Austronesian settlement. Antiquity (2005) 81:532–535.
Oppenheimer S. Eden in the East (1998) London: Wiedenfeld and Nicholson.
Paz V. Island southeast Asia—spread or friction zone? In: Examining the farming/language dispersal hypothesis—Renfrew C, ed. (2002) Cambridge: MacDonald Institute for Archaeological Research. 275–285.
Pelejero C, Kienast M, Wang L, Grimalt JO. The flooding of Sundaland during the last deglaciation: imprints in hemipelagic sediments from the southern South China Sea. Earth Planet Sci Lett (1999) 171:661–671.[CrossRef]
Pereira L, Richards M, Goios A. (13 co-authors). High-resolution mtDNA evidence for the late-glacial resettlement of Europe from an Iberian refugium. Genome Res (2005) 15:19–24.
Richards M, Macaulay V. The mitochondrial gene tree comes of age. Am J Hum Genet (2001) 68:1315–1320.[CrossRef][Web of Science][Medline]
Richards M, Macaulay V, Hickey E. (37 co-authors). Tracing European founder lineages in the Near Eastern mitochondrial gene pool. Am J Hum Genet (2000) 67:1251–1276.[Web of Science][Medline]
Richards M, Oppenheimer S, Sykes B. mtDNA suggests Polynesian origins in eastern Indonesia. Am J Hum Genet (1998) 63:1234–1236.[CrossRef][Web of Science][Medline]
Rolett BV, Tianlong J, Gongwu L. Early seafaring in the Taiwan Strait and the search for Austronesian origins. J East Asian Arch (2002) 4:307–319.[CrossRef]
Saillard J, Forster P, Lynnerup N, Bandelt H-J, Nørby SS. mtDNA variation among Greenland Eskimos: the edge of the Beringian expansion. Am J Hum Genet (2000) 67:718–726.[CrossRef][Web of Science][Medline]
Sathiamurthy E, Voris HK. Maps of Holocene sea level transgression and submerged lakes on the Sunda Shelf. The Nat Hist J Chulalongkorn Univ (2006) (Suppl 2):1–43.
Solheim WG. South-east Asia and Korea from the beginnings of food production to the first states. In: The history of humanity—De Laet SJ, ed. (1994) London: Routledge. 468–481.
Solheim WG. Archaeology and culture in Southeast Asia: unraveling the Nusantao (2006) Quezon City (Philippines): University of Philippines Press.
Su B, Jin L, Underhill P. (11 co-authors). Polynesian origins: insights from the Y chromosome. Proc Natl Acad Sci USA (2000) 97:8225–8228.
Tajima A, Sun C-H, Pan I-H, Ishida T, Saitou N, Horai S. Mitochondrial DNA polymorphisms in nine aboriginal groups of Taiwan: implications for the population history of aboriginal Taiwanese. Hum Genet (2003) 113:24–33.[CrossRef][Web of Science][Medline]
Tanaka M, Cabrera VM, González AM, et al, (28 co-authors). Mitochondrial genome variation in Eastern Asia and the peopling of Japan. Genome Res (2004) 14:1832–1850.
Tanudirjo DA. Islands in Between: prehistory of the Northeastern Indonesian Archipelago. (2001) Canberra (Australia): The Australian National University.
Terrell JE, Welsch RL. Lapita and the temporal geography of prehistory. Antiquity (1997) 71:548–572.[Web of Science]
Tommaseo-Ponzetta M, Attimonelli M, De Robertis M, Tanzariello F, Saccone C. Mitochondrial DNA variability of West New Guinea populations. Am J Phys Anthropol (2002) 117:49–67.[CrossRef][Web of Science][Medline]
Torroni A, Bandelt H-J, D'Urbano L. (11 co-authors). mtDNA analysis reveals a major late Paleolithic population expansion from southwestern to northeastern Europe. Am J Hum Genet (1998) 62:1137–1152.[CrossRef][Web of Science][Medline]
Torroni A, Bandelt H-J, Macaulay V. (33 co-authors). A signal, from human mtDNA, of post-glacial recolonization in Europe. Am J Hum Genet (2001a) 69:844–852.[CrossRef][Web of Science][Medline]
Torroni A, Rengo C, Guida V. (12 co-authors). Do the four clades of the mtDNA haplogroup L2 evolve at different rates? Am J Hum Genet (2001b) 69:1348–1356.[CrossRef][Web of Science][Medline]
Trejaut JA, Kivisild T, Loo JH, Lee CL, He CL, Xi JR, Li ZY, Lin M. Traces of archaic mitochondrial lineages persist in Austronesian speaking Formosan populations. Unambiguous linkage between Aboriginal Taiwanese and Polynesians to the exclusion of Mainland Asians. PLoS Biol (2005) 3:e247.[CrossRef][Medline]
van Holst Pellekaan SM, Frommer M, Sved JA, Boettcher B. Mitochondrial control-region sequence variation in aboriginal Australians. Am J Hum Genet (1998) 62:435–449.[CrossRef][Web of Science][Medline]
Veth P, Spriggs M, O'Connor S. Continuity in tropical cave use: examples from East Timor and the Aru Islands. Asian Perspect (2005) 44:180–192.[CrossRef]
Voris HK. Maps of Pleistocene sea levels in Southeast Asia: shorelines, river systems and time durations. J Biogeography (2000) 27:1153–1167.[CrossRef]
Wen B, Li H, Gao S. (18 co-authors). Genetic structure of Hmong-Mien speaking populations in East Asia as revealed by mtDNA lineages. Mol Biol Evol (2004a) 22:725–734.[CrossRef][Web of Science][Medline]
Wen B, Li H, Lu D. (18 co-authors). Genetic evidence supports demic diffusion of Han culture. Nature (2004b) 431:302–305.[CrossRef][Medline]
Wen B, Xie X, Gao S. (13 co-authors). Analyses of genetic structure of Tibeto-Burman populations reveals sex-biased admixture in southern Tibeto-Burmans. Am J Hum Genet (2004c) 74:856–865.[CrossRef][Web of Science][Medline]
Yang Z. PAML: a program package for phylogenetic analysis by maximum likelihood. CABIOS (1997) 13:555–556.[Medline]
Yao Y-G, Kong Q-P, Bandelt H-J, Kivisild T, Zhang Y-P. Phylogeographic differentiation of mitochondrial DNA in Han Chinese. American Journal of Human Genetics (2002a) 70:635–651.[CrossRef][Web of Science][Medline]
Yao Y-G, Nie L, Harpending H, Fu Y-X, Yuan Z-G, Zhang Y-P. Genetic relationship of Chinese ethnic populations revealed by mtDNA sequence diversity. American Journal of Physical Anthropology (2002b) 118:63–76.[CrossRef][Web of Science][Medline]
Yao Y-G, Zhang Y-P. Phylogeographic analysis of mtDNA variation in four ethnic populations from Yunnan Province: new data and a reappraisal. Journal of Human Genetics (2002) 47:311–318.[CrossRef][Web of Science][Medline]
Yao YG, Lu XM, Luo HR, Li WH, Zhang YP. Gene admixture in the Silk Road region of China: evidence from mtDNA and melanocortin1 receptor polymorphism. Genes and Genetic Systems (2000) 75:173.[CrossRef]
Zhang S. The epipaleolithic in China. J East Asian Arch (2000) 2:51–66.[CrossRef]
CiteULike 
Connotea 
Del.icio.us What's this?
This article has been cited by other articles:
|
|
 |
|
  D. H. Parks, M. Porter, S. Churcher, S. Wang, C. Blouin, J. Whalley, S. Brooks, and R. G. Beiko GenGIS: A geospatial information system for genomic data Genome Res., October 1, 2009; 19(10): 1896 - 1904. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
J. Xing, W. S. Watkins, D. J. Witherspoon, Y. Zhang, S. L. Guthery, R. Thara, B. J. Mowry, K. Bulayeva, R. B. Weiss, and L. B. Jorde Fine-scaled human genetic structure revealed by SNP microarrays Genome Res., May 1, 2009; 19(5): 815 - 825. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 R. D. Gray, A. J. Drummond, and S. J. Greenhill Language Phylogenies Reveal Expansion Pulses and Pauses in Pacific Settlement Science, January 23, 2009; 323(5913): 479 - 483. [Abstract] [Full Text] [PDF]
|
|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||