Caribbean Journal of Science, Vol. 39, No. 1, 11-22, 2003 Copyright 2003 College of Arts and Sciences University of Puerto Rico, Mayagu¨ez

Mitochondrial DNA Analysis in : Strong Maternal Ancestry of Closely Related Amerindians and Implications for the Peopling of Northwestern Venezuela

GLADYS TORO-LABRADOR1,OSWALD R. WEVER2, AND JUAN C. MARTI´NEZ-CRUZADO1,3

1Department of Biology, University of Puerto Rico, P.O. Box 9012, Mayagu¨ez, Puerto Rico 00681-9012 2Department of Internal Medicine, Dr. Horacio Oduber Hospital, Oranjestad, Aruba 3Corresponding author: [email protected]

ABSTRACT.—The continental origin of human mtDNA can be identified by its haplogroup determination through restriction fragment length polymorphism (RFLP) analysis. Hair root samples from 16 Aruban informed volunteers were analyzed by this method. Thirteen samples had mtDNAs of Amerindian origin and the remainder had their origin in sub-Saharan Africa; thus suggesting a substantial Amerindian maternal ancestry in Arubans, and helping explain the high incidence of health-related conditions common in Am- erindian populations. Nine of the thirteen Amerindian mtDNAs belonged to haplogroup D, suggesting that despite intense Indian population movements through historical times, most of the mtDNAs shared a recent common origin. Our results, combined with the pre- and post-Columbian history of Aruba and northwestern Venezuela, lead us to hypothesize that the frequency of haplogroup D mtDNAs should be high in the Caquetı´o tribe of coastal Falco´n in Venezuela and at least moderate in the Guajiro Indians of eastern Colombia. In addition, if the native Aruban Indians contributed substantially to the current mtDNA pool of Arubans despite the repeated historical deportation events that presumably removed the complete popula- tions from the island, they must also have possessed haplogroup D mtDNAs in high frequency. Our results do not rule out an origin for the Amerinds of northwestern Venezuela in ancient Amazonian ceramic cultures.

Because of its strict maternal inheritance, INTRODUCTION mtDNA studies only measure female in- High-resolution restriction endonuclease put. Thus, the results of such studies in analysis of the mitochondrial DNAs (mtDNAs) populations that have emanated from a of people from various continents, fol- colonization process are frequently not in lowed by phylogenetic reconstruction concordance with those obtained from analysis, has identified many groups of studies based on the Y chromosome, which evolutionarily related haplotypes (hap- displays strict paternal inheritance pat- logroups) that are mostly continent- terns. The results of combining data from specific. Thus, haplogroups can be defined mtDNA and Y chromosome studies typi- by sets of evolutionarily stable restriction cally reveal a predominance of asymmetric motifs that are common to all of their hap- matings between men from colonizing lotypes. Furthermore, the continental ori- populations and women from colonized gin of mtDNA obtained from a mixed ones (Pinto et al., 1996; Batista Dos Santos population can be determined by its hap- et al., 1999; Carvajal-Carmona et al., 2000). logroup identity. Hence, restriction frag- Aruban preceramic archaeological sites ment length polymorphism (RFLP) analy- belong to the Meso-Indian I period (Rouse ses have been used, sometimes in and Cruxent, 1963). The oldest dated site is combination with control region sequences, no more than 4500 years old. The Archaic to estimate the contribution of peoples settlers were apparently fisher-hunter- from different continents to the mtDNA gatherers organized in bands of 10 to 15 pool of mixed populations (Rando et al., people who sustained themselves prima- 1998; Alves-Silva et al., 2000; Green et al., rily on seafood (Versteeg, 1991). Polished 2000). axe-heads made of greenish-gray volcanic 11 12 G. TORO LABRADOR ET AL. material associated with the Meso-Indian large region that included the San Cristo´bal industry of Venezuela, and a sherd of a pot- area of the Andes, south of Lake Mara- tery disc found at the Canashitu preceramic caibo, and coastal areas from Lake Mara- site (500 BC-500 AD), suggest that the Ar- caibo to Cumana´as well as the islands of chaics traded with the ceramic cultures of Aruba, Curac¸ao, and Bonaire. Four subtra- Venezuela (Dahn, 1970). Because of their ditions developed from this expansion (Ol- proximity, we believe that the most likely iver, 1989), of which the Dabajuran subtra- candidates for such trade lived in the Para- dition took over coastal Falco´n, including guana´peninsula, located approximately 30 the Paraguana´peninsula and the islands of km south of Aruba. Aruba, Curac¸ao, and Bonaire (Fig. 1). The The most recent skeletal remains of Aru- Indians occupying this region were known ban Archaics have been found in Malmok as Caquetı´os by the Spaniards, and their (600 AD-900 AD). These remains bear clear language (Caquetı´o) belongs to the Arawa- distinctions when compared with those of kan family of languages (Oliver, 1989). the ceramic culture that arrived approxi- The dominant hypothesis is that the mately at 1000 AD. The Archaics possessed Arawakan family of languages originated high, narrow, and long dolichocephalic approximately 5500 years ago in the junc- skulls, while those from the ceramic culture tion between the Negro and Amazon rivers were low, wide, and short (brachyce- (Rouse, 1986). Proto-Arawakan speakers phalic). Furthermore, all Archaic skulls had migrated north, upstream along the Negro shovel-shaped incisors while such teeth are River, through the Casiquiare Canal, and rare in skulls from the ceramic culture (Ver- downstream along the Orinoco River. steeg, 1993; Tacoma, 1991). Splits occurred at the Meta-Orinoco and the The ceramic culture of Aruba belonged Apure-Orinoco junctions. Some of these an- to the Dabajuran subtradition of the cient people migrated further downstream Dabajuroid tradition. Because of similar el- along the Orinoco River, all the way to the ements in their pottery, the Dabajuroid and coast, and to the Lesser and Greater An- Tierroid traditions of western Venezuela tilles, while others migrated west, up- are believed to share a common ancestral stream along the Meta or the Apure River. culture that may have existed somewhere Those that migrated up the Meta River between the Turbio, Upper Cojedes, and reached the Andes. Some of those migrat- Upper Yaracuy rivers (Oliver, 1995). The ing up the Apure River diverged up the Dabajuroid tradition expanded in the Mar- Cojedes River and eventually reached the acaibo basin after 445 BC; it occupied a coast migrating downstream along the Ya-

FIG. 1. Geographic distribution of the Dabajuran subtradition (shaded area). MITOCHONDRIAL DNA ANALYSIS IN ARUBA 13 racuy River. The latter may have given rise the presence of Indians on the islands (Ver- to the Macro-Dabajuroid tradition (Oliver, steeg and Ruiz, 1995). Because of the com- 1989). plexity of the Aruba cave labyrinths, it is The Dabajurans of Aruba developed an possible that the Indians were mostly na- agrarian lifestyle. Thus, whereas prece- tives who had escaped deportation, but ramic sites have large shell middens, sug- they could have been recent migrants from gesting that the Archaics consumed many the mainland. In addition, substantial bivalves and oysters (Versteeg and Ruiz, mainland to Aruba migrations of Indian es- 1995), ceramic sites have comparatively capees occured from 1529 to 1556 during fewer shells (Versteeg and Rostain, 1997) the development of the Venezuelan colony among which those of the queen conch, (Haviser, 1991). Strombus gigas, predominate (Boerstra, Aruba was neglected by the Spaniards 1983). from 1533 until the Dutch conquest of 1636, There is strong evidence that Aruban when Spanish and Amerindian languages Dabajurans maintained close contacts with (especially Caquetı´o) were widely spoken. mainland cultures. For example, close to Upon the Dutch conquest the Spaniards 700 chert flakes were found at the Tanki fled, and all the Indians were deported to Flip ceramic site (Heidecker and Siegel, the mainland because they were regarded 1969), and since no chert source is found in as sympathetic to the Spaniards (Hartog, Aruba, these cherts must have been im- 1961). However, in that same year the ported. In addition, the avemorphic motifs Dutch West Indian Company assigned painted in pottery found in Aruba seem to Aruba the duty of breeding horses and have originated in non-Dabajuroid cul- cattle, and Indians were chosen for these tures. These motifs likely came from the La endeavors because they had a good repu- Guajira peninsula of Colombia, where the tation as wild-horse hunters. Thus, Indian oldest site containing artifacts with these migration to Aruba was stimulated. Also, motifs occurs in the Rancherı´a Valley. Ap- some Indians in war with Spaniards west of parently, only later did the avemorphic Maracaibo fled to Aruba (Hartog, 1961). motifs made it to the Venezuelan The importance of Aruba diminished af- Dabajuran site of Los Me´danos (Oliver, ter the 1648 Dutch Peace Treaty with Spain 1989), suggesting that the Aruban Indians and the island was neglected again. In 1655, traded not only with their Dabajuran kins the Dutch West Indian Company recog- in Venezuela, but with other cultures as nized free Indians in Aruba as trade part- well. ners. These Indians were assigned a piece It is estimated that some 450 to 600 Indi- of land on which to maintain themselves ans lived in Aruba at the time of the Span- through cultivation; they also cut and sold ish discovery (Dijkhoff, 1997) in 1499. To- wood and exploited marine resources (Ver- gether with Curac¸ao and Bonaire, Aruba steeg and Ruiz, 1995). In his description of was declared an island without use in 1513, the Aruban way of life during the second and two years later some 2000 Caquetı´o In- half of the 17th century, A.O. Exquemelin dians from the three islands combined points out that the Indians spoke Spanish, were transported to Hispaniola to work in were Catholic, and were visited frequently mines (Hartog, 1961). These Indians pre- by Spanish priests from the mainland (Har- sumably comprised the entire population tog, 1961). As an example of their strong of the islands (Versteeg and Ruiz, 1995), but links with the mainland, some 200 Indians 150 to 200 were returned to Aruba and Cu- agreed to leave Aruba in 1723 to raise the rac¸ao starting in 1526 to work on the expor- Venezuelan town of El Carrizal under the tation of brasilwood, kwihi, and divi-divi. ecclesiastic jurisdiction of the city of Coro. The Indians returned to Aruba and Cura- These Indians were described by the Span- c¸ao were mainly Caquetı´os, but some iards as Caquetı´os (Fortique, 1989). Arawaks from other islands The migration of non-Indians to Aruba were included in the group (Hartog, 1961). was scant. From 1640 to 1754, Aruba was to Upon their return, the Spaniards noted be colonized only with the permission of 14 G. TORO LABRADOR ET AL. the Dutch West Indian Company director. consistently lower (i.e. 7 % of 2536 in 1957) A.O. Exquemelin described the Lieutenant (De Ruyter, 1957) and even 5 % of 474 in as the only white man 1968 and 8 % of 1027 in 2000 (O.R. Wever, living on the island, and he was presum- unpubl. data). Blood-group O prevalences ably the only person who owned black were never intermediate between African slaves. African slaves are mentioned for the (47 %) and European (44 %) values but first time in the history of Aruba in 1750, were instead consistently higher (52 % of when they revolted and killed four resi- 474 in 1968, 54 % of 1577 in 1996 and 49 % dents (Alofs, 1996). Aruban Indians had of 1027 in 2000: O.R. Wever, unpubl. data). their own slaves. They raided or bought Hence, Aruban prevalences for Rhesus- captive mainland Indian minors of both negative (cde) and for blood-group O are sexes called red slaves (Hartog, 1961). intermediate between African values and The migration of non-Indians increased Amerindian values (cde 0 %; blood-group O substantially under a 1785 land tax that 100 %). made it easier to get a license to settle (Har- Although the 1816 census scores 33 % of tog, 1961). Settlers from various European the Aruba inhabitants as full-blooded Indi- countries moved to Aruba with their black ans, after 1820 only whites, blacks and col- slaves, but these migrations were stymied ored people are mentioned as ethnic com- by a famine in 1800 that occurred as a con- ponents (Hartog, 1961; Alofs, 1996; sequence of drought and war. Aruba Dijkhoff, 1997). We hope to show that this shifted between English and Dutch hands poorly defined category of colored people several times from 1803 to 1816, when the not only covered Afro-European mulattos English recapitulated. but mostly included mestizos and Amerin- Migration and admixture increased dra- dians. The persistence of a substantial pro- matically during the 19th and 20th centu- portion of Amerindian genetic ancestry in ries. The Caquetı´o language disappeared Arubans may explain the current high from Aruba during the first half of the 19th prevalences of health-related conditions century, and oral tradition holds that the that are common in Amerindians, such as last full-blooded Aruban Indian died diabetes (Muneta et al., 1993; Ramachan- around 1862 (Hartog, 1961). However, the dran, 1994; Lee et al., 1995), pterygium red slave trade also increased during the (Hilgers, 1959), and lactase deficiency (Cas- 19th century until the abolition of in key et al., 1977; Newcomer et al., 1977). 1863. Baptism records show that red slaves We provide further evidence of a strong during this period were mainly Guajiro In- Amerindian ancestry among Arubans, dians (Nooyen, 1965). The 1816 census of showing that most Arubans have mtDNAs Aruba reports 1732 inhabitants classified of Amerindian origin. We also show that into 584 colored free people (probably the Amerindian mtDNAs have low genetic mixed), 564 full-blooded Indians, 133 col- diversity and that their haplogroup distri- ored (red) slaves, 240 blacks, and 211 bution is highly structured, suggesting that whites. Women in these categories were the Amerindian mtDNA pool of Aruba is 313, 307, 48, 137, and 101, respectively homogeneous despite the intense migra- (Hartog, 1961). tion events that occurred during historical This historic data shows that Arubans times. Finally, we compare our mtDNA re- are a mixed population, much like those of sults with those of various South American many Latin American nations. Rhesus and tribes, as well as with linguistic and ar- ABO blood group findings suggested that a chaeological evidence, in an attempt to ad- trihybrid admixture pattern including Am- vance our understanding on the pre- erindian genes explains better the observed Columbian events that gave rise to the Aruban blood-group distributions than the people found by the Spaniards in Aruba. bihybrid Afro-European model of the Car- ibbean. Rhesus-negative (cde) prevalences MATERIALS AND METHODS were never intermediate between African Samples.—Hair follicle samples were ob- (9 %) and European (16 %) values but were tained with informed consent from 16 un- MITOCHONDRIAL DNA ANALYSIS IN ARUBA 15 related Aruban patients at the medical of- were heated at 95°C for 2.5 min before be- fice of Dr. Oswald R. Wever in Oranjestad, ing subjected to 32 cycles of 30 s at 94 °C, 1 Aruba. The mothers and maternal grand- min at 54°C, and 70 s at 72°C. Three restric- mothers of 15 participants were born in tion sites not tested by Martı´nez-Cruzado Aruba, while those of the other participant et al. (2001) were tested in this study; these were born in Cuba. The mothers of the par- define the sub-Saharan African hap- ticipants were all born between 1903 and logroups L3d (−8616 DpnII) and L3e (+2349 1940 with a mean at 1923. It was not pos- DpnII) (Rando et al., 1998), which are sible to date with precision the birthdays of equivalent to haplogroups L3b and L3a of the maternal grandmothers. Chen et al. (2000), respectively, and the Experimental strategy.—The haplogroups Caucasian haplogroup I (−4529 HaeII). The found here were described initially as 8616 and 2349 DpnII sites were tested from monomorphic at the 10394 DdeI and 10397 a 140 bp fragment amplified using primers AluI sites (Chen et al., 1995; Wallace, 1995). L8558 (5’ TCTGTTCGCTTCATTCATTG 3’) This fact combined with the proximity of and H8657 (5’ TGATTAGTCATTGT- these two restriction sites, which allows TGGGTG 3’) and from a 187 bp fragment analysis of both sites produced from a amplified using primers L2272 (5’ TCCT- single amplicon, makes their testing a con- CACACCCAATTGGAC 3’) and H2420 (5’ venient start in the haplogroup identifica- ATGCCTGTGTTGGGTTGAC 3’), respec- tion of mtDNAs obtained from a mixed tively. The 4529 HaeII site was tested from a population. By identifying the 10394 DdeI/ 199 bp fragment amplified using primers 10397 AluI motif of an unknown mtDNA, L4462 (5’ AAAATGTTGGTTATACCCTTC the number of candidate haplogroups is 3’) and H4620 (5’ ATGGCAGCTTCTGTG- quickly reduced. The major haplogroups GAAC 3’). classified by their 10394 DdeI/10397 AluI To confirm the haplogroup of the motifs and their specific markers and geo- mtDNA identified as belonging to hap- graphic origins are presented in Table 1. logroup A, the sequence of its HV-I region We grouped our samples by their 10394 was determined. An 879 bp fragment con- DdeI/10397 AluI motif and tested all the taining the HV-I region was amplified us- samples within each group for particular ing the primers L15766 (5’ ATTCTAACCT- haplogroup-specific markers. Because GAATCGGAGG 3’) and H34 (5’ some haplogroup I mtDNAs possess the ACCAAATGCATGGAGAGCTCC 3’). The haplogroup L3d-specific marker, the iden- amplicon was purified using the High Pure tification of mtDNAs as belonging to hap- PCR Product Purification Kit (Roche Mo- logroup L3d was confirmed by showing lecular Biochemicals) as per manufacturer’s that they did not possess the haplogroup instructions. Six hundred nanograms of the I-specific marker. Since one Italian hap- product and 10 uL of each sequencing logroup U mtDNA possessing the hap- primer (1 ␮M concentration) were sent to logroup A-specific marker has been found the University of Medicine and Dentistry of (Torroni et al., 1997), the identification of a New Jersey, New Jersey Medical School mtDNA as belonging to haplogroup A was Molecular Resource Facility, for automated confirmed through the sequence determi- sequencing. An Applied Biosystems (ABI) nation of its HV-I region. model 3100 capillary sequencer was used Experimental procedure.—Total nucleic ac- after cycle sequencing with Dye Terminator ids were released from the hair follicles and mix version 2.0. The primers used for specific mtDNA fragments were amplified sequencing were L15854 (5’ CCTAATCC- and subjected to restriction analyses. Pro- TAATACCAACTATC 3’) and H16526 (5’ cedures were performed as in Martı´nez- GGGAACGTGTGGGCTATTTAGG 3’). Di- Cruzado et al. (2001) with the following ex- versity for the Amerindian mtDNAs was ceptions: 1. The amplification reaction calculated using the method of Tajima mixtures contained 7 uL of sample DNA, (1989), h =[1−∑xi2]n/(n − 1), where xi is the bovine serum albumin at 0.3 mg/mL, and 2 frequency of each haplogroup and n is the units of Taq DNA polymerase. 2. Mixtures population size. 16

TABLE 1. Major RFLP-defined haplogroups, their origins, 10394 DdeI/10397 AluI motifs, and haplogroup-specific markers.

10394 DdeI/10397 Haplogroup-specific 10394 DdeI/ Haplogroup-specific AluI motif Haplogroup Origin1 marker2 10397 AluI motif Haplogroup Origin marker

(−/−) A Native American +663 HaeIII (+/+) C Native American +13262 AluI AL. ET LABRADOR TORO G. (−/−) B Native American 9 bp deletion (+/+) D Native American −5176 AluI Region V (−/−) F Asian −12406 HpaI (+/+) E Asian −7598 HhaI (−/−) H Caucasian −7025 AluI (+/+) G Asian +4830 HaeII (−/−) K3 Caucasian +12308 HinfI/ (+/−) I Caucasian −4529 HaeII −9052 HaeII (−/−) T Caucasian +4216 NlaIII (+/−) J Caucasian +4216 NlaIII (−/−) U Caucasian + African +12308 HinfI/ (+/−) K Caucasian +12308 HinfI/−9052 HaeII +9052 HaeII (−/−) V Caucasian −4577 NlaIII (+/−) L1 & L2 African +3592 HpaI (−/−) W Caucasian +8249 AvaII (+/−) L3b African +10084 TaqI (−/−) X Caucasian +14465 AccI (+/−) L3d African −8616 DpnII + Native (+/−) L3e African +2349 DpnII American 1The Native American haplogroups have their origin in Asia. 2Numbers refer to the nucleotide position according to the Cambridge Reference Sequence (Anderson et al., 1981). 3Haplogroup K is the only one with two 10394 DdeI/AluI motifs in moderate frequencies. MITOCHONDRIAL DNA ANALYSIS IN ARUBA 17

RESULTS 2349 DpnII site showed the −8616 DpnII motif and thus belonged to haplogroup A hierarchical strategy (Fig. 2) was used L3d. Because several of the mtDNAs be- to identify the mtDNA haplogroups of the longing to the Caucasian haplogroup I that 16 samples. Each sample was first tested for have been analyzed for their entire coding the 10394 DdeI/10397 AluI motif. Two sequences show a transversion at position samples lacked both restriction sites and 8616 that eliminates the 8616 DpnII site each was tested for the markers specific for (Finnila¨et al., 2001; Herrnstadt et al., 2002), haplogroups A and B; one sample had the it was necessary to confirm the identity of 663 HaeIII site but not the 9 bp deletion at the mtDNA belonging to haplogroup L3d. region V, and thus belonged to haplogroup This was done by testing and showing that A, while the other had the inverse profile it lacked the −4529 HaeII haplogroup I- and thus belonged to haplogroup B. The specific marker. identity of the haplogroup A mtDNA was Eleven samples possessed the 10394 DdeI confirmed with the DNA sequence of its and the 10397 AluI sites and all were tested HV-I region. The sequence presented the for the 13262 AluI site that identifies hap- transitions at positions 16111, 16223, 16290, logroup C and for the lack of the 5176 AluI and 16319 that characterize virtually all site specific for haplogroup D. Nine haplogroup A mtDNAs of Native Ameri- samples lacked both AluI sites and thus be- can origin (Ward et al., 1991; Horai et al., longed to haplogroup D; the other two pos- 1993; Torroni et al., 1993; Santos et al., 1994; sessed both sites and thus belonged to hap- Starikovskaya et al., 1998). logroup C. The test for the 5176 AluI site Three samples had the 10394 DdeI site that identified eight of the haplogroup D but lacked the 10397 AluI site. These samples is shown in Fig. 3. All results are samples were tested first for the 3592 HpaI summarized in Table 2. site that characterizes haplogroups L1 and L2 (which represent over 70 % of all sub- DISCUSSION Saharan Africans) but the three samples lacked the site. They were then tested for The abundant literature available on the the 2349 DpnII site that is specific for hap- human mtDNA variation and continent- logroup L3e. Two of the samples, including specific polymorphisms has allowed us to the one of Cuban ancestry, had the site and use a hierarchical strategy to quickly iden- thus belonged to haplogroup L3e. The tify the haplogroup and the continental ori- three samples were also tested for the gin of each of the mtDNA samples. This −8616 DpnII motif that characterizes hap- exempted us from the need to perform tests logroup L3d. Only the sample lacking the on all samples for the presence of Cauca-

FIG. 2. Hierarchical strategy used to identify the haplogroups of the 16 samples. The 10394 DdeI/10397 AluI motif was tested first. Samples (−/−) for the motif were tested for the markers specific for haplogroups A (+663 HaeIII) and B (Region V 9 bp deletion). Those (+/−) for the motif were tested for the 3592 HpaI site shared by haplogroups L1 and L2, and for the markers specific for haplogroups L3e (+2349 DpnII) and L3d (−8616 DpnII). The samples (+/+) for the motif were tested for the markers specific for haplogroups C (+13262 AluI) and D (−5176 AluI). 18 G. TORO LABRADOR ET AL.

FIG. 3. Test for the 5176 AluI site. Fragments 149 bp long were amplified using primers L5121 and H5229 (Martı´nez-Cruzado et al., 2001), subjected to AluI digestion and fractionated on 3% agarose gels. Lanes 1 and 15: phiX174-HaeIII molecular weight marker (New England Biolabs). Lanes 2 to 4: 149 bp control fragments not subjected to AluI digestion. Lane 5: blank. Lanes 6 to 14: experimental fragments. Only the sample in lane 7 was cut by AluI at the 5176 site, producing two undistinguishable fragments of approximately equal size. All other samples lacked the AluI site and thus belong to haplogroup D.

TABLE 2. Haplogroup distribution those possessing two of the haplogroup- specific markers reported here. MtDNAs Haplogroup Ancestry Frequency with haplogroup-specific markers of two A Native American 1 (6.3)1 haplogroups which shared the same 10394 B Native American 1 (6.3) DdeI/10397 AluI motif were found in only C Native American 2 (12.5) three cases. Specifically, 9 haplogroup I D Native American 9 (56.3) mtDNAs were found to possess the −8616 L3d sub-Saharan African 1 (6.3) L3e sub-Saharan African 2 (12.5) DpnII motif (Torroni et al., 1994c, 1997; Herrnstadt et al., 2002) that defines hap- 1Numbers in parenthesis represent percent frequen- logroup L3d. In addition, 4 haplogroup A cies. mtDNAs had the 9 bp deletion (Torroni et al., 1993, 1994b) that defines haplogroup sian haplogroups HV, I, JT, and W, the sub- B. Finally, 1 haplogroup U mtDNA con- Saharan African haplogroup L3b, the Cau- tained the +663 HaeIII haplogroup A- casian-Native American haplogroup X, and specific marker (Torroni et al., 1997). Thus, the Caucasian-African haplogroup U. our results were confirmed by testing our It could be argued that because of the haplogroup L3d mtDNA for the hap- relatively high mutation rate of the coding logroup I-specific marker, testing the two region of the human mtDNA, estimated at mtDNAs with a (−/−) motif for the mark- 2.2 to 2.9 % per million years (Torroni et al., ers specific for haplogroups A and B, and 1994a), all markers used to define any hap- sequencing the HV-I region of the hap- logroup could arise by an independent mu- logroup A mtDNA. tation in mtDNAs belonging to other From the interviews, we learned that 15 haplogroups, thus producing the of the 16 mtDNAs studied were present in misidentification of the mtDNAs. How- Aruba by the first decades of the 20th cen- ever, such events are very rare. In combi- tury (when the mothers of the participants nation, Finnila¨et al. (2001) and Herrnstadt were born) and that 13 (86.7 %) of these et al. (2002) sequenced the complete coding were of Amerindian origin, suggesting a region of 627 Caucasian, 69 Native Ameri- substantial Amerindian maternal ancestry can or Asian, and 56 African mtDNAs. In in Arubans. This is not surprising because addition, 418 mtDNAs of Native American the 1816 census category of colored free origin, 397 of Caucasian, and 214 of African people probably consisted of people of origin have been subjected to high- mixed ancestry. Typically, admixtures in resolution restriction endonuclease analy- the New World have consisted of men from sis (Torroni et al., 1993, 1994b, 1994c, 1996, colonizing and women from colonized 1997; Chen et al., 1995, 2000; Huoponen et populations, thus producing mixed popu- al., 1997). We searched these mtDNAs for lations with Amerindian mtDNAs in very MITOCHONDRIAL DNA ANALYSIS IN ARUBA 19 high frequencies (Green et al., 2000; Carva- versities are common in tribes of the New jal-Carmona et al., 2000; Alves-Silva et al., World, suggesting that many of them were 2000; Martı´nez-Cruzado et al., 2001). Thus, founded by small groups of closely related it is reasonable to speculate that the num- women. However, only 11 of the 55 New ber of women with Amerindian mtDNAs World tribes listed by Merriwether et al. in 1816 was equal or higher than the sum of (1995) have a haplogroup distribution (as the colored free people, the red slaves, and measured by the frequency of the most the full-blooded Indians. This adds to 668 common haplogroup) as highly structured out of 906 (73.7 %) women. The remainder as that of the Aruba population. Further- was composed of 137 (15.1 %) “blacks” and more, only 9 of the 25 Colombian tribes 101 (11.1 %) “whites”. studied by Keyeux et al. (2002) were calcu- The absence of mtDNAs of Caucasian lated to have genetic diversities lower than origin in our samples probably stems from that displayed in the Aruba study. These sampling error. However, it is noteworthy observations suggest that despite the in- that 90 of the 101 women counted in the tense movement of Indians in and out of 1816 census as whites were born in Aruba, Aruba during historical times, most Aruba leaving open the possibility that many of Amerindian mtDNAs originate from one or them were mixed. For comparison, while from a few sibling tribes. more than 80 % of the population of Puerto The most likely candidate is the Caquetı´o Rico was classified as “white” in the 2000 tribe. During the first years of colonization, census, more than half and over a quarter the native Indians of Aruba were described of this same population has mtDNAs of by the Spaniards as Caquetı´os (Hartog, Amerindian and sub-Saharan African ori- 1961). In addition, the Caquetı´os in the gin, respectively (Martı´nez-Cruzado et al., mainland were the tribe geographically 2001, and unpubl. data). Thus, a high Am- closest to Aruba, and archaeological evi- erindian or sub-Saharan African mtDNA dence points towards close ties between frequency among those women classified both groups during pre-Columbian times as white in the 1816 Aruba census cannot (Oliver, 1989). Thus, the most likely sce- be ruled out. nario is that most of the Aruban Indians Although our results show only a small that returned to Aruba from Hispaniola window of the universe that existed in were Caquetı´os, as were also most of those Aruba at the beginning of the 20th century, that were found in Aruba at this time. The the observed frequency is very similar to latter were likely either deportation escap- what may have been expected from an ees or mainland migrants. This scenario is equivalent study done in 1816. Hence, our supported by the fact that the Indians that results suggest that most women who mi- migrated from Aruba to the mainland in grated to Aruba during the 19th century 1723 to found the town of El Carrizal were had Amerindian maternal ancestry. This is described as Caquetı´os (Fortique, 1989). supported by experiments showing that the Some of the mtDNAs studied could have Rhesus-negative and blood-group O preva- a Guajiro Indian origin because many lences in Aruba are consistently out of the Guajiros were sent as slaves to Aruba dur- ranges imposed by African and European ing the 19th century (Nooyen, 1965). We prevalences, and more towards those im- may hypothesize that haplogroup D is pre- posed by Amerindian prevalences (O.R. dominant among Caquetı´os and at least Wever, unpubl. data). Thus, the mtDNA re- very common among Guajiros, suggesting sults presented here should be substanti- that both tribes are genetically closely re- ated through a study using a sample set lated. We expect that similar studies in Ca- better representing the Aruba population. quetı´os and Guajiros, and in ancient re- The haplogroup distribution of the Am- mains from these two tribes and Aruba will erindian mtDNAs of Aruba (Table 2) is reveal high haplogroup D frequencies. highly structured, and thus their genetic di- Similar haplogroup distributions do not versity (calculated at 0.5256) is low. Highly necessarily imply strong genetic relation- structured distributions and low genetic di- ships. Strictly speaking, the finding that 20 G. TORO LABRADOR ET AL. two Amerindians share a haplogroup D nario, the Caquetı´os would be located close mtDNA only means that their most recent to the northwestern edge of this expansion, common maternal ancestor cannot be older as the Wayuu tribe in the Guajira peninsula than the age of the haplogroup D woman of Colombia displays clear genetic differ- who crossed the Bering Strait some 29 000 ences. The latter tribe has haplogroups A years ago (Forster et al., 1996; Silva et al., and B at a 53 % combined frequency and 2002). However, we could expect a closer lacks haplogroup D (Mesa et al., 2000; relationship between these tribes due to Keyeux et al., 2002). their geographic proximity and because, al- It is intriguing that almost 70 % of the though distinct, the Caquetı´o and Guajiran Amerindian mtDNAs in Aruba belong to languages belong to the Arawakan family the same haplogroup, despite strong ar- (Rouse, 1986). It is noteworthy that among chaeological evidence indicating that the South American Indians, genetic variation ceramic people who came to Aruba ca. 1000 is more strongly correlated with linguistic AD were physically distinct to the Archaics than with geographic variation (Fagundes (Versteeg, 1991, 1993; Tacoma, 1991). A ge- et al., 2002). Median network analyses netically homogeneous population would (Bandelt et al., 1995, 1999) of mitochondrial suggest the absence of any admixture be- HV-I and HV-II region sequences are nec- tween the ceramic and preceramic peoples, essary to ascertain a common origin for the and thus the complete prehistoric replace- tribes. ment of one population by another. De- Because of their language kinship, it tailed analyses of haplogroup D mtDNAs could be postulated that all Arawakan from Caquetı´os might reveal two or more speaking tribes north of the Amazon River groups within haplogroup D that could share a recent genetic ancestry, and that re- represent independent populations whose peated founder events stemming from a admixture gave rise to the Caquetı´os. In common ancestor population generated de- combination with mtDNA studies on an- rived populations with the markedly dif- cient remains, these studies should lead us ferent mtDNA haplogroup distributions towards a better understanding of the peo- found in some Arawakan speaking tribes. pling and culture development of the re- For example, haplogroup D is absent in gion. Puerto Rico, where haplogroups A and C constitute over 91 % of the Amerindian mtDNAs (Martı´nez-Cruzado et al., 2001). Presumably, these mtDNAs derive mainly Acknowlegments.—We are grateful to do- from the extinct Taı´nos, whose Taı´no lan- nors of hair follicles used in this project. Dr. guage belonged to the Arawakan family of Anna Maria K.E. Wever-Merks and Dr. languages. The Amazonian origin of most John Kock (dermatologist) took most hair of the haplogroup C mtDNAs in Puerto samples. We thank Raymundo Dijkhoff Rico is supported by their −7013 RsaI motif (Aruban Museum of Archaeology) for pro- (Martı´nez-Cruzado, unpubl. data), which viding valuable information and revising had been found only in Amazonian tribes our statements about the pre- and post- (Torroni et al., 1993a). Although 69.2 % of Columbian history of Aruba and north- the Amerindian mtDNAs of Aruba belong western Venezuela. We also thank Dr. John to haplogroup D, the common origin of Uscian and two anonymous reviewers for both populations in the ancient proto- their useful comments on the manuscript, Arawakan populations of the Amazon is and Robert Donnelly (University of Medi- supported by the high frequencies of hap- cine and Dentistry of New Jersey, New Jer- logroups C (45.5 %) and D (34.8 %) in tribes sey Medical School Molecular Resource Fa- north of the Amazon River, such as the Ma- cility) for the oligonucleotide primers used kiritare, the Yanomama, the Macushi, the in this work and for the automated DNA Marubo, and the Wapishana (Torroni et al., sequencing services. This study was sup- 1993), even when not all of these tribes ported by a National Science Foundation speak Arawakan languages. Under this sce- grant to JCMC. MITOCHONDRIAL DNA ANALYSIS IN ARUBA 21

LITERATURE CITED logenetic network for European mtDNA. Am. J. Hum. Genet. 68: 1475-1484. Alofs, L. 1996. Slaven zonder plantage. Slavernij en Forster, P., R. Harding, A. Torroni, and H.J. Bandelt. emancipatie op Aruba, 1750-1863. Charuba, Edi- 1996. Origin and evolution of Native American cion Educativo no.2, Aruba. mtDNA variation: a reappraisal. Am. J. Hum. Alves-Silva, J. et al. 2000. The ancestry of Brazilian Genet. 59: 935-945. mtDNA lineages. Am. J. Hum. Genet. 67: 444-461. Fortique, J.R. 1989. Noticias sobre el Carrizal, pueblo Anderson, S. et al. 1981. Sequence and organization of de indios caquetı´os. Boletı´n de la Academia Nacio- the human mitochondrial genome. Nature 290: nal de la Historia, 72: 57-80, Caracas, Venezuela. 457-465. Green, L.D., J.N. Deer, and A. Knight. 2000. mtDNA Bandelt, H.J., P. Forster, B.C. Sykes, and M.B. Rich- affinities of the peoples of North-Central Mexico. ards. 1995. Mitochondrial portraits of human Am. J. Hum. Genet. 66: 989-998. populations using median networks. Genetics 141: Hartog, J. 1961. Aruba past and present, from the time 743-753. of the Indians until today. De Wit, Aruba, 451 pp. Bandelt, H.J., P. Forster, and A. Rohl. 1999. Median- Haviser, J.B. 1991. The first Bonaireans. Reports of the joining networks for inferring intraspecific phylog- Archaeological-Anthropological Institute of the enies. Mol. Biol. Evol. 16: 37-48. Netherland Antilles no.10, Curac¸ao. Batista Dos Santos, S.E., J.D. Rodrigues, A.K.C. Heidecker, L. and M.I. Siegel. 1969. Preliminary report Ribeiro-Dos-Santos, and M.A. Zago. 1999. Differ- on the excavation Henrı´quez-I Site, Tanki Flip, ential contribution of indigenous men and women Aruba, Netherland Antilles. Fla. Anthropol. 22: 12- to the formation of an urban population in the 26. Amazon region as revealed by mtDNA and Y- Herrnstadt, C. et al. 2002. Reduced-median-network DNA. Am. J. Phys. Anthropol. 109: 175-180. analysis of complete mitochondrial DNA coding- Boerstra, E.H.J. 1983. Some of the soil marks in the region sequences for the major African, Asian, and Tanki Flip excavation, Aruba, N.A.: the ditches. European haplogroups. Am. J. Hum. Genet. 70: Proceedings of the Ninth International Congress 1152-1171. for the Study of the Pre-Columbian Cultures of the Hilgers, J.H.C. 1959. Pterygium on the island of Lesser Antilles: 173-183. Aruba. Ph.D. thesis, Amsterdam. Carvajal-Carmona, L.G. et al. 2000. Strong Amerind/ Horai, S. et al. 1993. Peopling of the Americas, white sex bias and a possible Sephardic contribu- founded by four major lineages of mitochondrial tion among the founders of a population in north- DNA. Mol. Biol. Evol. 10: 23-47. west Colombia. Am. J. Hum. Genet. 67: 1287-1295. Huoponen, K. et al. 1997. Mitochondrial DNA and Y Caskey, D.A. et al. 1977. Effects of age on lactose mal- chromosome-specific polymorphisms in the Semi- absorption in Oklahoma Native Americans as de- nole tribe of Florida. Eur. J. Hum. Genet. 5: 25-34. termined by breath H2 analysis. Digest. Diseases Keyeux, G., C. Rodas, N. Gelvez, and D. Carter. 2002. 22: 113. Possible migration routes into de- Chen, Y.-S., A. Torroni, L. Excoffier, A.S. Santachiara- duced from mitochondrial DNA studies in Colom- Benerecetti, and D.C. Wallace. 1995. Analysis of bian Amerindian populations. Hum. Biol. 74: 211- mtDNA variation in African populations reveals 233. the most ancient of all human continent-specific Lee, E.T. et al. 1995. Diabetes and impared glucose haplogroups. Am. J. Hum. Genet. 57: 133-149. tolerance in three American Indian populations Chen, Y.-S. et al. 2000. mtDNA variation in the South aged 45-74 years: the strong heart study. Diabetes African Kung and Khwe - and their genetic rela- Care 18: 599-610. tionship to other African populations. Am. J. Hum. Martı´nez-Cruzado, J.C. et al. 2001. Mitochondrial Genet. 66: 1362-1383. DNA analysis reveals substantial Native American Dahn, J. 1970. Ecological and archaeological presenta- ancestry in Puerto Rico. Hum. Biol. 73: 491-511. tion of the Canashitu site in Aruba, Netherland Merriwether, D.A., F. Rothhammer, and R.E. Ferrell. Antilles. Rapport 14-9-1970. 1995. Distribution of the four founding lineage De Ruyter, J.A.M. 1957. Analysis of some results in haplotypes in Native Americans suggests a single hospital obstetrics in Aruba. Ph.D. thesis, Leiden wave of migration for the New World. Am. J. Phys. University, Holland. Anthropol. 98: 411-430. Dijkhoff, R.A.C.F. 1997. Tanki Flip / Henrı´quez-I: an Mesa, N.R. et al. 2000. Autosomal, mtDNA, and Y- early Urumaco site in Aruba. Master thesis, Leiden chromosome diversity in Amerinds: pre- and post- University, Holland. Columbian patterns of gene flow in South Fagundes, N.J.R., S.L. Bonatto, S.M. Callegari-Jacques, America. Am. J. Hum. Genet. 67: 1277-1286. and F.M. Salzano. 2002. Genetic, geographic, and Muneta, B., J. Newman, S. Wetterall, and J. Stevenson. linguistic variation among South American Indi- 1993. Diabetes and associated risk factors among ans: possible sex influence. Am. J. Phys. Anthropol. native Americans. Diabetes Care 16: 1619-1620. 117: 68-78. Newcomer, A.D., P.J. Thomas, D.B. McGill, and A.F. Finnila¨, S., M.S. Lehtonen, and K. Majamaa. 2001. Phy- Hofmann. 1977. Lactase deficiency: a common ge- 22 G. TORO LABRADOR ET AL.

netic trait of the American Indian. Gastroenterol- tral mutation hypothesis by DNA polymorphism. ogy 72: 234. Genetics 123: 585-595. Nooyen, R.H. 1965. Millefiori di Aruba. Stichting Torroni A. et al. 1993. Asian affinities and continental Aruba Nostra, 115 pp. radiation of the four founding Native American Oliver, J.R. 1989. The archaeological, linguistic and mtDNAs. Am. J. Hum. Genet. 53: 563-590. ethno-historical evidence for the expansion of Torroni A., J.V. Neel, R. Barrantes, T.G. Schurr, and Arawakan into northwestern Venezuela and D.C. Wallace. 1994a. Mitochondrial DNA “clock” northeastern Colombia. Ph.D. thesis, University of for the Amerinds and its implications for timing Illinois, Urbana-Champaign. their entry into North America. Proc. Natl. Acad. Oliver, J.R. 1995. The archaeology of Maticora, west- Sci. USA 91: 1158-1162. ern Venezuela: a prehistoric longhouse structure in Torroni A. et al. 1994b. mtDNA and Y-chromosome site Fal-7. Proceedings of the Fifteenth Interna- polymorphisms in four Native American popula- tional Congress for Caribbean Archaeology: 1-52. tions from southern Mexico. Am. J. Hum. Genet. Pinto, F., A.M. Gonza´lez, M. Herna´ndez, J.M. Larruga, 54: 303-318. and V.M. Cabrera. 1996. Genetic relationship be- Torroni A. et al. 1994c. mtDNA and the origin of Cau- tween the Canary Islanders and their African and casians: identification of ancient Caucasian-specific Spanish ancestors inferred from mitochondrial haplogroups, one of which is prone to a recurrent DNA sequences. Ann. Hum. Genet. 60: 321-330. somatic duplication in the D-loop region. Am. J. Ramachandran, A. 1994. Increasing prevalence of Hum. Genet. 55: 760-776. NIDDM. A public health problem of developing Torroni A. et al. 1996. Classification of European nations. Dialogue 11-15. mtDNAs from an analysis of three European popu- Rando, J.C. et al. 1998. Mitochondrial DNA analysis of lations. Genetics 144: 1835-1850. northwest African populations reveals genetic ex- Torroni A. et al. 1997. Haplotype and phylogenetic changes with European, Near-Eastern, and sub- analyses suggest that one European-specific Saharan populations. Ann. Hum. Genet. 62: 531- mtDNA background plays a role in the expression 550. of Leber Hereditary Optic Neuropathy by increas- Rouse, I. 1986. Migrations in prehistory: inferring ing the penetrance of the primary mutations 11778 population movement from cultural remains. Yale and 14484. Am. J. Hum. Genet. 60: 1107-1121. University Press. New Haven, Connecticut, 202 pp. Versteeg, A.H. 1991. The Indians of Aruba. Publica- Rouse, I. and J.M. Cruxent. 1963. Venezuelan archae- tions of the Archaeological Museum of Aruba no.4, ology. Yale University Press, New Haven, Con- Aruba. necticut, 179 pp. Versteeg, A.H. 1993. Malmok revisited: geometric lay- Santos, M., R.H. Ward, and R. Barrantes. 1994. out and patterns in a preceramic cemetery. Pro- mtDNA variation in the Chibcha Amerindian ceedings of the Fourteenth International Congress Hue´tar from Costa Rica. Hum. Biol. 66: 963-977. for Caribbean Archaeology: 115-124. Silva, W.A. et al. 2002. Mitochondrial genome diver- Versteeg, A.H. and S. Rostain. 1997. The archaeology sity of Native Americans supports a single early of Aruba: the Tanki Flip site. Publication of the entry of founder populations into America. Am. J. Archaeological Museum of Aruba no.8, Aruba, 519 Hum. Genet. 71: 187-192. pp. Starikovskaya, Y.B., R.I. Sukernik, T.G. Schurr, A.M. Versteeg, A.H. and A.C. Ruiz. 1995. Reconstructing Kogelnik, and D.C. Wallace. 1998. mtDNA diver- brasilwood island: the archaeology and landscape sity in Chukchi and Siberian Eskimos: implications of Indian Aruba. Publications of the Archaeologi- for the genetic history of ancient Beringia and the cal Museum of Aruba no.6, Aruba. peopling of the New World. Am. J. Hum. Genet. Wallace, D.C. 1995. Mitochondrial DNA variation in 63: 1473-1491. human evolution, degenerative disease and aging. Tacoma, J. 1991. Precolumbian human skeletal re- Am. J. Hum. Genet. 57: 201-223. mains from Curac¸ao, Aruba and Bonaire. Proceed- Ward, R.H., B.L. Frazier, K. Dew-Jager, and S. Pa¨a¨bo. ings of the Thirteenth International Congress for 1991. Extensive mitochondrial diversity within a Caribbean Archaeology, Part II: 802-813. single Amerindian tribe. Proc. Natl. Acad. Sci. USA Tajima, F. 1989. Statistical method for testing the neu- 88: 8720-8724.