- Research article
Genetic variation analysis of the Bali street dog using microsatellites
BMC Geneticsvolume 6, Article number: 6 (2005)
Approximately 800,000 primarily feral dogs live on the small island of Bali. To analyze the genetic diversity in this population, forty samples were collected at random from dogs in the Denpasar, Bali region and tested using 31 polymorphic microsatellites. Australian dingoes and 28 American Kennel Club breeds were compared to the Bali Street Dog (BSD) for allelic diversity, heterozygosities, F-statistics, GST estimates, Nei's DA distance and phylogenetic relationships.
The BSD proved to be the most heterogeneous, exhibiting 239 of the 366 total alleles observed across all groups and breeds and had an observed heterozygosity of 0.692. Thirteen private alleles were observed in the BSD with an additional three alleles observed only in the BSD and the Australian dingo. The BSD was related most closely to the Chow Chow with a FST of 0.088 and also with high bootstrap support to the Australian dingo and Akita in the phylogenetic analysis.
This preliminary study into the diversity and relationship of the BSD to other domestic and feral dog populations shows the BSD to be highly heterogeneous and related to populations of East Asian origin. These results indicate that a viable and diverse population of dogs existed on the island of Bali prior to its geographic isolation approximately 12,000 years ago and has been little influenced by domesticated European dogs since that time.
Bali, a province of the Republic of Indonesia, is an island just 87 km from north to south and 142 km from east to west and home to more than 2.9 million people . Approximately 800,000 stray dogs (Fig. 1) also live on the island based on a survey conducted by the Bali Street Dog Foundation (personal communication). Only a small percentage of these dogs live in homes or are provided routine veterinary care .
More than 90% of the residents of Bali are Hindu  with myth and ritual playing a vital part of daily life . The dog is also an important part of Balinese life and mythology. A popular tale from the Mahabharata  describes King Yudisthira's journey to Heaven's Gate, and his love for a dog that befriended him on his arduous and tragic journey (Fig. 2). As a direct result of such mythology, BSDs are treated with a degree of reverence and are often provided ceremonial food offerings . The deliberate killing of street dogs is not typically practiced, because Balinese people believe that all things should be allowed to die naturally . These cultural mores have contributed to the current overpopulation of dogs on the island.
As a result of overpopulation, many BSDs suffer from chronic skin diseases, internal parasites, parvo- and distemper-virus infections, and malnutrition. In an effort to reduce the dog population and to care for their medical needs, the Bali Street Dog Foundation (Yayasan Yudisthira Swarga) was founded in 1998 . They provide emergency care, treatment for skin disease and parasites, sterilization, public education on the plight of feral dogs, and improved veterinarian training. Twenty to 30 dogs are sterilized each day, with more than 9,000 dogs sterilized to date.
The BSD population is of interest for both its genetic diversity and historical relationships. It is also a population that has bred more or less randomly for thousands of years with limited genetic influx, due mainly to geographic barriers and a strict rabies control program in effect since 1926. The present study is concerned with the genetic diversity of this unique canine population and its relationship to other canine subpopulations in Asia and throughout the world. Data presented herein was derived from the DNA testing of 40 BSD samples from the Denpasar city region of Bali with 31 polymorphic microsatellite loci. The genetic diversity of the BSD was compared to that of the Australian dingo and 28 American Kennel Club (AKC) breeds.
Analysis of locus diversity across all 30 subpopulations revealed that the number of observed alleles ranged from six to 20 with a total of 366 for all loci (Table 1). Overall heterozygosity of the loci was high, with an average of 0.779, and all but four loci having HT values greater than 0.700. Average HS was 0.577 for the 30 subpopulations, with all but three loci having HS values greater than 0.500. The HS and HT values were closest for C23.123 and farthest for C22.279 and C10.404. HWE analysis revealed that all but one locus had at least one population out of equilibrium for the 30 populations sampled. C01.424, C31.646 and CPH16 had 7 populations out of HWE and AHT130 did not have any populations with p values below 0.05. The level of locus diversity attributable to subpopulation structure was evaluated with two statistics – RST and FST. Both statistics gave similar average values at 0.230 and 0.236 respectively. However, RST ranged from 0.098 to 0.486 while FST ranged from 0.179 to 0.328.
Bali street dog diversity
Overall, the BSD was the most genetically diverse population surveyed here, displaying 239 total alleles out of the 366 seen in all 30 subpopulations or 65.3% of the total observed alleles (Table 2). The Australian dingo displayed 144 alleles and the AKC breeds displayed 138.8 on average. Analysis of expected (HE) and observed (HO) heterozygosities (Table 2) revealed that the BSD had a 44.0% higher HE than the Australian dingo (0.736 vs. 0.511) and a 28.4% higher HE than the average AKC breed (0.573). HO was also highest in the BSD at 0.692, versus 0.426 in the Australian dingo and 0.563 in the average AKC breed.
In order to evaluate the bias of sampling twice the number of BSDs, twenty samples were taken at random from the total pool of 40 for bootstrap determinations, and the number of observed alleles, HE, and HO were calculated. This process was repeated for 10,000 iterations and the average value for each measurement was determined. The average bootstrap value for the number of observed alleles for 20 BSDs was 214.6. The average bootstrapped HE and HO values for the BSD were 0.727 and 0.692 respectively.
To understand the loss of approximately 24 observed alleles after bootstrapping, the allele frequencies for the BSD at each locus was examined. While the BSD had the highest number of alleles, they also had the highest number of alleles with a frequency below 5% (67 out of 239, data not shown).
FIS estimates were calculated to assess the level of inbreeding for each subpopulation (Table 2). The BSD had the lowest value at 0.097 and the Australian dingo had the highest at 0.194 with the average of AKC breeds at 0.137.
Allele frequency analysis also revealed that 10 private alleles were observed in the BSD as well as three alleles shared only with the Australian dingo (Table 3). The majority of private alleles in the BSD were below 5% in frequency with the exception of AHT121 where the private alleles had a combined frequency of 18.75%. The BSD and the Australian dingo also shared three alleles not seen in any AKC breed at locus C10.404 with a combined frequency of 16.25% in the BSD and 75% in the Australian dingo.
Several additional unique alleles were found only in the BSD, Australian dingo, Chow Chow and Akita; demonstrating a closer relationship of the BSD to Asian versus non-Asian dogs (data not shown). Appearing at the highest frequency was allele 201 of locus CPH08 in the BSD, Australian dingo and Chow Chow at frequencies of 21.3%, 10% and 20%, respectively. Further, the BSD, Australian dingo, Chow Chow and Akita share allele 113 of locus C22.279 at frequencies of 23.8%, 67.5%, 15% and 5%, respectively. Results for locus PEZ08 demonstrate a lack of influence of European alleles where a high frequency of deviations from n+4 alleles were observed in the AKC breeds sampled yet no deviation from n+4 alleles was observed in the BSD or the Australian dingo.
A pairwise FST analysis was performed for each locus between the BSD and the two closest subpopulations: the Australian dingo and Chow Chow (Table 1). The BSD was most similar to the Australian dingo at locus CPH16 with a FST of 0.005 and to the Chow Chow at locus C23.123 with a FST of 0.002. The BSD was most dissimilar to the Australian dingo and Chow Chow at locus RVC1 with a FST of 0.355 and 0.216 respectively. Several loci had similar distances for both population pairs, such as locus C01.424 or locus C20.446 and may indicate areas of the genome that are neutral to either environmental or human selection.
Genetic distance relationships
Further distance analysis was performed for all 31 loci between all 30 subpopulations using both Nei's DA distance and pairwise FST estimates (Table 4). Across all loci, the BSD shared allele frequencies most closely with the Chow Chow (DA = 0.242, FST = 0.088) and the Australian dingo (DA = 0.242, FST = 0.126), and least closely with the Airedale Terrier (DA = 0.454, FST = 0.258).
Genetic distance relationships amongst the five Asian subpopulations were further explored using neighbor-joining dendograms with four non-Asian subpopulations for comparison (Fig. 3). The BSD, Chow Chow, Australian dingo and Akita clustered together in 90% of the trees. The BSD, Chow Chow and Australian dingo further clustered in 87% of the trees. The BSD and Australian dingo maintained their relationship within the larger cluster in 84% of the trees. In the remainder of the tree, the Rhodesian Ridgeback, Greyhound, Airedale Terrier and Borzoi maintained a relationship in 51% of the trees and the Airedale Terrier/Borzoi cluster was seen in 63% of the trees. The Pug did not maintain a relationship with any other breed in this analysis, but was intermediate to the Asian and non-Asian subpopulations.
Microsatellites have been previously used to assess genetic diversity and relationships in feral dog subpopulations [6, 7]. Kim et al.  found that HO was high in three feral dog subpopulations of Korea, Sakhalin and Taiwan, ranging from 0.539 in the Taiwanese to 0.717 in the Korean dogs. Given that the loci used in that study had an average allele number of 7.75, these values are similar to the HO of 0.692 observed in the BSD. Wilton et al.  surveyed a population of Australian dingoes and found an average HO of 0.387 using microsatellites with an average allele number of 6.93, similar to the HO of 0.426 for the Australian dingoes reported herein with an average allele number of 11.8.
Given the size of the island of Bali, it is extraordinary that 800,000 feral dogs can thrive and maintain such high levels of genetic diversity. Of all the subpopulations surveyed here, the BSD has the highest number of observed alleles, the highest heterozygosity, the fewest number of loci out of Hardy-Weinberg equilibrium and the lowest FIS. Even after adjusting for sample size, the BSD maintains their status as the most heterogeneous population in the study. Unlike the Australian dingo which exhibits a much lower level of diversity, the BSD findings suggest either a large founding population on Bali and/or a consistent genetic influx since the geographic isolation of ~12,000 years ago. This data also supports that the BSD appears to approximate a randomly breeding population with little selection pressure.
When comparing the heterogeneity of the BSD to that observed within the AKC breeds some caveats should be addressed. One may initially expect long established, well-defined dog breeds to be much less heterogeneous than reported here. While some breeds do have a low HE, such as the Boxer with a HE of 0.320, breeds like the Jack Russell Terrier have a high HE of 0.713 and overall their HE is higher than that of the dingo. Of first note, the selection of the dogs that contribute to a breed composition mostly occurs prior to official breed recognition primarily by genetic drift due to geographic isolation and selection for particular working or physical characteristics. After official breed recognition future breeding choices are based primarily on the availability of sires and dams that approximate the breed standard. As a result, there is a founding population that proceeds to breed mostly by convenience. Also, many breed standards have changed considerably over the years resulting in retention of a certain level of diversity within each breed, some breeds retaining much more than others. Finally, dogs comprising the comparison AKC breeds were sampled from across the United States, removing any geographical bias of the genotypes observed and slightly elevating the heterozygosities.
The average allelic diversity of the loci used in the present study was 11.8 alleles per locus, versus 7.75 in the Kim et al work . However, the average number of alleles observed is 4.6 among the subpopulations in the present study and the average HT is 0.577. The average values for the 11 subpopulations surveyed in the Kim et al  work were 4.34 and 0.547, respectively. The higher total allelic diversity in the present study is likely due to the fact that nearly three times more subpopulations were studied.
RST and FST values were nearly identical across all subpopulations and all loci, indicating that approximately 23% of the differences observed in allele frequencies can be attributed to differences between subpopulations. FST provides an unbiased estimate of genetic drift between subpopulations by comparing alleles identical by state. RST takes advantage of the stepwise mutation model, which assumes that mutations most often occur as whole repeat unit losses or gains from the original allele size. As a result, the number of mutations provides an estimate of time from divergence. It is interesting, therefore, to compare RST and FST values by locus. Eighteen of the 31 loci studied have an RST to FST ratio greater than 1.1 (Table 1) indicating that the populations have been separated for a sufficient amount of time for mutations to impact genetic structure. An interesting exception is observed at CPH16 where the ratio is 0.420. CPH16 may have a mutation pattern where both stepwise additions and subtractions occur at equal and high frequency. Of note, the average pairwise RST value between the BSD and each of the 29 comparison subpopulations is 0.056 at locus CPH16. The highest RST to FST ratio occurs at locus CPH03 with a value of 1.724. Interestingly, the BSD and the Australian dingo have a pairwise RST value of 0.017 at CPH03, whereas the average value of the BSD compared to the other 28 subpopulations has a value of 0.254. The distance between the BSD and the Australian dingo at CPH03 may support that those two populations were isolated most recently from each other relative to the other 28 subpopulations.
Bali street dog origin
The origin of the people of Bali is clouded by myth and a scarcity of archeological findings. Therefore, the origin of the dog on Bali is also speculative. Nonetheless, a hypothesis can be formed based on known human and dog histories. Current evidence points to an early migration of humans from Africa through Indonesia and into Australia approximately 60,000 to 70,000 years ago [8, 9]. Recent excavations have also revealed that there was a great expansion into Indonesia from China between 4,000 and 5,000 years ago that could have contributed to a population pre-existing on Bali . Supportive evidence that Indonesia was populated prior to 5,000 years ago is a higher degree of heterogeneity in the Indonesian population than seen in the North Asian population, suggesting that the Indonesia was populated earlier than regions to the North . The "Slow Boat Model" for the peopling of Polynesia also suggests a prolonged mixing of Southeast Asians with Indonesians, which predated migration to the East . In short, Indonesia appears to be a human genetic melting pot with genetic influences over tens of thousands of years.
The dog on the island of Bali may also be a parallel "canine genetic melting pot." While the domestication date of the dog is in much dispute , approximately 14,000 years ago is accepted as a late date. During the earliest human migrations through Indonesia however, it is highly possible that wolf packs or feral dogs traveled the same routes, establishing a feral population on Bali in the process. Even if humans were not capable of taming the dog at that time, dogs could still have benefited from close proximity to humans. Figure 4 shows a superimposition of the proposed geographic origin for five Asian and four non-Asian dog subpopulations presented herein and the major theorized human migration routes. It is noteworthy that the BSD, Chow Chow and Australian dingo, related breeds by genetic analysis, all share one proposed human migration route.
If a feral dog population was established on the island of Bali more than 14,000 years ago, then that population became isolated approximately 10,000 years ago when the sea levels drastically rose, submerging the land bridges of the Indonesia archipelago . Geographic isolation was unlikely to have been absolute; genetic diversity of the BSD was invariably enhanced at various times by the influx of new dogs. At the time humans migrated to Indonesia from China, dogs were known to be domesticated and undoubtedly accompanied people as companions . Mitochondrial DNA sequencing evidence suggests that the dingo was introduced into Australia about that time from the Indonesian archipelago [15, 8, 9]. Bali's documented history of repeated war and trade spanning the last 2,000 years [1, 16, 17] represents actions that are often associated with the introduction of new animals. Indeed, a somewhat free movement of dogs probably occurred up to 1926, when the import of dogs to Bali was greatly curtailed as a means to prevent the introduction of rabies . This policy greatly reduced, though not eliminated, new outside introductions of new dogs to the island. In contrast to the Australian dingo population, which appears to have undergone a severe population bottleneck or founder effect based on microsatellite alleles and mtDNA , the BSD population maintains a high level of genetic variation. There is no evidence for a genetic bottleneck or small founding population for the BSD.
The relatedness of the BSD to the Australian dingo and the Chow Chow is evidenced by common unique alleles and allele frequencies despite the very different levels of genetic diversity between the subpopulations. According to the hypothesis presented herein, one could imagine that feral dog subpopulations were established throughout Indonesia with much mixing until ~12,000 years ago. At that time, each population became closed with little influx of new genetic material until humans migrated south from Asia between 4,000 and 5,000 years ago. The degree of influx since that period would have been influenced by the frequency of trade and conflict, factors determined by accessibility, available natural resources, and political structure. The island of Bali is historically a less visited island than it's neighbor Java and therefore the indigenous dog population would have been subjected to less influence.
This study into the diversity and relationship of the BSD to other domestic and feral dog populations shows the BSD to be highly diverse and related to populations of East Asian origin. These results indicate that a viable and diverse population of dogs existed on the island of Bali prior to its geographic isolation approximately 12,000 years ago and has been little influenced by domesticated European dogs since that time. It would be of interest to study feral subpopulations on other islands in the archipelago to determine if the same level of diversity is observed elsewhere, or if the situation on Bali is truly unique. Y-chromosome, mitochondrial and MHC marker typing on the BSD, as well as feral dogs from other regions, would help to determine if indeed dogs followed the same migration routes as their likely human companions.
BSDs were randomly captured and taken to a BSD Foundation field clinic for treatment or sterilization and simultaneously sampled for DNA collection with buccal swabs. Familial relationships of the BSDs sampled could not be easily determined; therefore the sample population was doubled (40 vs. 19–20 samples) over that of other study groups. Blood samples from the Australian dingo were taken from captive animals in Australia. Australian dingoes were known to be unrelated by at least one generation.
Dogs from 28 American Kennel Club (AKC) breeds, equally representing the AKC group designations, were sampled with buccal swabs for a previous study . Twenty dogs were tested for each breed, with the exception of two breeds (Doberman Pinscher and the Border Collie) that comprised 19 individuals. The 28 breeds included were: Airedale Terrier, Akita, American Eskimo, Australian Shepherd, Belgian Tervuren, Bernese Mountain Dog, Border Collie, Borzoi, Boxer, Brittany, Bull Terrier, Bulldog, Chow Chow, Doberman Pinscher, Golden Retriever, Greyhound, Jack Russell Terrier, Keeshond, Labrador Retriever, Miniature Bull Terrier, Norwegian Elkhound, Papillon, Pembroke Welsh Corgi, Pomeranian, Pug, Rhodesian Ridgeback, Weimaraner, and Yorkshire Terrier. Dogs within each breed were unrelated by at least one generation.
Thirty-one of the 100 microsatellites multiplexed into 12 PCRs by the Veterinary Genetics Laboratory  had been previously used to evaluate the Australian dingo samples (unpublished data). For comparison purposes, those same 31 microsatellites were selected for use in the present study. All markers but one (PEZ02) were mapped on either the 1999 canine genetic linkage map  or the Radiation hybrid map . Loci selected for study represented 25 of the 38 autosomes of the dog, with five autosomes represented by two loci. The average distance for the markers on chromosomes CFA06, CFA11, CFA20 and CFA23 is 23.5 cM and 23.4 Mb between AHT139 and RVC1 on CFA15. As a result, only 25 loci are known to be unlinked. PEZ02 has not been mapped and may be linked to a marker in the study.
Forward primers were synthesized and dye labeled with either Fam, Hex or Vic, or Tamra or Ned (Applied Biosystems, Inc. (ABI), Foster City, CA). Reverse primers were synthesized by Operon (Alameda, CA). Primer sequences and concentrations for all markers are available as Additional file 1.
Sample preparation and PCR conditions
BSD and AKC breed DNA was derived from buccal cells harvested from the inside of the cheek with nylon bristle cytology brushes (Medical Packaging Corp., Camarillo, CA). Samples were collected by owners or field volunteers and submitted directly to the laboratory. DNA was extracted by heating a single swab for 10 min at 95°C in 400 μl 50 mM NaOH and then neutralized with 140 μl 1 M Tris-HCl, pH 8.0. Australian dingo DNA was extracted from blood using a standard sodium hydroxide digest.
A 2 μl aliquot of extract was used in each PCR which equates to approximately 50 ng DNA. All markers and DNAs were amplified with a PCR reagent mix of 1X PCR buffer (ABI), 4.17 mM MgCl2, 200 μM of each dNTP (Hoffmann-La Roche Inc, Nutley, NJ), 0.6 unit AmpliTaq (ABI), and 2% DMSO (Sigma) then covered with 15 ul Chill-out™ Liquid Wax (MJ Research, Inc., Waltham, MA) to prevent evaporation. One of five thermal cycler programs was used for each primer mix ranging from 56° to 64° degrees for the annealing temperature. All PCR work was done in polycarbonate 96-well v-bottom microtiter plates (USA Scientific, Ocala, FL) on MJ Research PTC-100 thermal cyclers (MJ Research, Inc., Waltham, MA). Protocols are also available in Additional 1.
Gel electrophoresis conditions and DNA fragment analysis
One μl aliquots of PCR product were mixed with 2 μl Fluorescent Ladder (CXR) 60–400 (Promega 400) or Internal Lane Standard 600 (Promega 600) (Promega, Madison, WI) fluorescent size standard, denatured on MJ Research PTC-100 thermal cyclers for three minutes at 95°C, then held at 5°C or placed on ice for at least one minute before gel loading. Two μl aliquots were then loaded onto a 6% denaturing polyacrylamide gel and run on an ABI 377 Automated Sequencer using ABI 10" × 7 1/8" short plates (12 cm). Gels were run at 1.10 kV (constant) voltage, 60.0 mA current, 200 W power, 51°C and 40.0 mW (constant) laser power for up to 2 hours when using Promega 400, and up to 3 hours using Promega 600. DNA fragment analysis was performed with in-house designed STRand software , which replaces ABI Genotyper and Genescan software. This data was then transferred to an in-house database compatible with the STRand software.
Allelic diversity and observed heterozygosities (HO) were determined by direct counting for each of the 30 subpopulations. Hardy-Weinberg equilibrium (HWE) tests were performed using Genepop version 3.4 . Pairwise FST estimates and subpopulation expected heterozygosities (HE) for the 30 breeds or dog groups were performed using Genepop version 3.4 . FIS estimates (inbreeding coefficient of each subpopulation) for each allele following Weir and Cockerham  were calculated using Genepop version 3.4 and are presented as averages across all loci.
Gene diversity or total population heterozygosity (HT) and its associated parameters, HS (average heterozygosity among subpopulations) and GST (coefficient of genetic differentiation), were calculated across all loci using the public domain software, DISPAN . Two additional measures of variance, FST  and RST [27, 28] were calculated using Genepop version 3.4. A pairwise genetic distance matrix using Nei's DA distance was also created using DISPAN with bootstrapping. Genotype data for all populations is available in Additional file 2.
Phylogenetic tree construction
Allele frequencies were used to compute a matrix of genetic distances , which were then used to construct a phylogenetic tree of relationships among 5 Asian and 4 non-Asian dog subpopulations. Takezaki's  POPTREE program was used to create a neighbor joining tree using DA distances with 1000 bootstrap replications. The output of POPTREE was then converted to the New Hampshire format for editing in the stand alone program TREEVIEW version 1.6.6  and bootstrap values were added.
Bali Street Dog
FST RST, GST: F-statistics indices
HT, HE, HO: Heterozygosity indices
Number of alleles
American Eskimo Dog
Bernese Mountain Dog
Miniature Bull Terrier
Jack Russell Terrier
Pembroke Welsh Corgi
Hobart A, Ramseyer U, Leemann A: The peoples of Bali. 1996, Oxford: Blackwell Publishers
Yayasan Yudisthira Swarga (Bali Street Dog Foundation). [http://www.yamp.com/balidogs/]
Vishva Hindu Parishad. [http://www.vhp.org/englishsite/d.Dimensions_of_VHP/qVishwa%20Samanvya/vishvahinduparishadabroad.htm]
Buck W: Mahabhrata. 1973, Berkeley: University of California Press
Department of Foreign Affairs of the Republic of Indonesia. State Gazette. 1912, , No. 452 and Ministry of Agriculture Decree No. 892/Kpts/The Nation.56/1997
Kim KS, Tanabe Y, Park CK, Ha JH: Genetic variability in East Asian dogs using microsatellite loci analysis. Journal of Heredity. 2001, 92 (5): 398-403. 10.1093/jhered/92.5.398.
Wilton AN, Steward DJ, Zafiris K: Microsatellite variation in the Australian dingo. Journal of Heredity. 1999, 90 (1): 108-11. 10.1093/jhered/90.1.108.
Dayton L: Modern human origins meeting: "Tracing the road down under". Science. 2003, 302 (5645): 555-Oct 24
Dayton L: Modern human origins meeting: "On the trail of the first Australian dingo". Science. 2003, 302 (5645): 555-6. Oct 24
Faradz SM, Pattiiha MZ, Leigh DA, Jenkins M, Leggo J, Buckley MF, Holden JJ: Genetic diversity at the FMR1 locus in the Indonesian population. Annals of Human Genetics. 2000, 64 (Pt 4): 329-39. 10.1046/j.1469-1809.2000.6440329.x.
Gibbons A: The peopling of the Pacific. Science. 2001, 291 (5509): 1735-7. 10.1126/science.291.5509.1735. Mar 2
Pennisi E: Canine evolution: A shaggy dog history. Science. 2002, 298 (5598): 1540-2. 10.1126/science.298.5598.1540. Nov 22
Hanebuth T, Stattegger K, Grootes PM: Rapid flooding of the sunda shelf: A late-glacial sea-level record. Science. 2000, 288 (5468): 1033-5. 10.1126/science.288.5468.1033. May 12
Savolainen P, Zhang YP, Luo J, Lundeberg J, Leitner T: Genetic evidence for an East Asian origin of domestic dogs. Science. 2002, 298 (5598): 1610-3. 10.1126/science.1073906. Nov 22
Savolainen P, Leitner T, Wilton AN, Matisoo-Smith E, Lundeberg J: A detailed picture of the Origin of the Australian dingo, obtained from the study of Mitochondrial DNA. Proceedings of the National Academy of Sciences of the USA. 2004, 101 (33): 12387-12390. 10.1073/pnas.0401814101.
Ardika IW: Archaeological research in northeastern Bali, Indonesia. PhD thesis. 1991, Australian National University
Hartaningsih N, Dharma DMN, Rudiyanto MD: Anjing Bali: Pemuliabiakan dan pelestarian (Bali Street Dogs: breeding and conservation). 1999, Yogyakarta: Kanisius Press
Wilton AN: DNA Methods of Assessing Australian dingo Purity. A Symposium on the Australian dingo. Edited by: Dickman CR, Lunney D. 2001, Mosman NSW: Royal Zoological Society of New South Wales, 49-55.
Irion DN, Schaffer AL, Famula TR, Eggleston ML, Hughes SS, Pedersen NC: Analysis of genetic variation in 28 dog breed subpopulations with 100 microsatellite markers. Journal of Heredity. 2003, 94 (1): 81-7. 10.1093/jhered/esg004.
Eggleston ML, Irion DN, Schaffer AL, Hughes SS, Draper JE, Robertson KR, Millon LV, Pedersen NC: PCR multiplexed microsatellite panels to expedite canine genetic disease linkage analysis. Animal Biotechnology. 2002, 13 (2): 223-235. 10.1081/ABIO-120016191.
Neff MW, Broman KW, Mellersh CS, Ray K, Acland GM, Aguirre GD, Ziegle JS, Ostrander EA, Rine J: A second-generation genetic linkage map of the domestic dog, Canis familiaris. Genetics. 1999, 151 (2): 803-20.
Guyon R, Lorentzen TD, Hitte C, Kim L, Cadieu E, Parker HG, Quignon P, Lowe JK, Renier C, Gelfenbeyn B, Vignaux F, DeFrance HB, Gloux S, Mahairas GG, Andre C, Galibert F, Ostrander EA: A 1-Mb Resolution Radiation Hybrid Map of the Canine Genome. Proceedings of the National Academy of Sciences of the USA. 2003, 100: 5286-5291. 10.1073/pnas.0831002100.
STRand Nucleic Acid Analysis Software. [http://www.vgl.ucdavis.edu/STRand]
Raymond M, Rousset F: GENEPOP (version 1.2): population genetics software for exact tests and ecumenicism. Journal of Heredity. 1995, 86: 248-249.
Weir BS, Cockerham CC: Estimating F-statistics for the analysis of population structure. Evolution. 1984, 38: 1358-1370.
Ota T: Program DISPAN: Genetic distance and phylogenetic analysis. 1993, Copyright of the Pennsylvania State University, University Park, PA, [http://iubio.bio.indiana.edu:7780/archive/00000022/]
Rousset F: Equilibrium values of measure of population subdivision for stepwise mutation processes. Genetics. 1996, 142: 1357-1362.
Michalakis Y, Excoffier L: A generic estimation of population subdivision using distances between alleles with special interest to microsatellite loci. Genetics. 1996, 142: 1061-1064.
Nei M: Molecular evolutionary genetics. 1987, New York: Columbia University Press
Takezaki N: POPTREE version for DOS. [http://www.bio.psu.edu/People/Faculty/Nei/Lab/Programs.html]
Page RDM: TreeView version 1.6.6. Distributed by the Division of Environmental and Evolutionary Biology, Institute of Biomedical and Life Sciences. 2001, University of Glasgow, Glasgow G12 8QQ, Scotland, UK, [http://taxonomy.zoology.gla.ac.uk/rod/treeview.html]
We thank the Yayasan Yudisthira Swarga (Bali Street Dog Foundation) and director Sherry Grant for bringing this interesting dog population to our attention and their effort in sampling feral dogs.
DNI performed the majority of data acquisition and analysis, wrote first draft of the manuscript and prepared the final draft for submission. ALS performed the majority of sample processing, assisted in data acquisition and the writing of the subsequent drafts of the manuscript as well as final draft preparation. SG sampled the dogs tested, provided background for the manuscript and assisted in the final draft preparation. ANW provided the Australian dingo data for comparison and assisted in the subsequent drafts of the manuscript. NCP directed the research and assisted in the writing of the manuscript. All authors read and approved the final manuscript.