Open Access Highly Accessed Research article

Dog Y chromosomal DNA sequence: identification, sequencing and SNP discovery

Christian Natanaelsson1, Mattias CR Oskarsson1, Helen Angleby1, Joakim Lundeberg1, Ewen Kirkness2 and Peter Savolainen1*

Author Affiliations

1 School of Biotechnology, KTH, Royal Institute of Technology, AlbaNova University Center, 10691 Stockholm, Sweden

2 The Institute for Genomic Research (TIGR), Rockville, MD 20850, USA

For all author emails, please log on.

BMC Genetics 2006, 7:45  doi:10.1186/1471-2156-7-45

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2156/7/45


Received:10 July 2006
Accepted:6 October 2006
Published:6 October 2006

© 2006 Natanaelsson et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Population genetic studies of dogs have so far mainly been based on analysis of mitochondrial DNA, describing only the history of female dogs. To get a picture of the male history, as well as a second independent marker, there is a need for studies of biallelic Y-chromosome polymorphisms.

However, there are no biallelic polymorphisms reported, and only 3200 bp of non-repetitive dog Y-chromosome sequence deposited in GenBank, necessitating the identification of dog Y chromosome sequence and the search for polymorphisms therein. The genome has been only partially sequenced for one male dog, disallowing mapping of the sequence into specific chromosomes. However, by comparing the male genome sequence to the complete female dog genome sequence, candidate Y-chromosome sequence may be identified by exclusion.

Results

The male dog genome sequence was analysed by Blast search against the human genome to identify sequences with a best match to the human Y chromosome and to the female dog genome to identify those absent in the female genome. Candidate sequences were then tested for male specificity by PCR of five male and five female dogs.

32 sequences from the male genome, with a total length of 24 kbp, were identified as male specific, based on a match to the human Y chromosome, absence in the female dog genome and male specific PCR results. 14437 bp were then sequenced for 10 male dogs originating from Europe, Southwest Asia, Siberia, East Asia, Africa and America. Nine haplotypes were found, which were defined by 14 substitutions. The genetic distance between the haplotypes indicates that they originate from at least five wolf haplotypes. There was no obvious trend in the geographic distribution of the haplotypes.

Conclusion

We have identified 24159 bp of dog Y-chromosome sequence to be used for population genetic studies. We sequenced 14437 bp in a worldwide collection of dogs, identifying 14 SNPs for future SNP analyses, and giving a first description of the dog Y-chromosome phylogeny.

Background

Population genetic studies of the earliest history of the domestic dog have so far mainly been based on the analysis of mitochondrial DNA (mtDNA), which represents a single genetic marker and, since it is maternally inherited, can only describe the history of females [1-4]. In order to have a second marker to corroborate the results based on mtDNA and to obtain a picture also for the history of male dogs, studies of the dog Y chromosome would be valuable. The Y-chromosomal genetic variation has been analysed in a few studies so far, but only for microsatellite variation [5,6]. Microsatellites have the disadvantage that mutations may be recurrent, rendering the phylogenetic relations between haplotypes not entirely clear. Single nucleotide polymorphisms (SNPs) on the other hand are an, in practice, digital form of information, since multiple mutations in the same nucleotide position is unlikely. It has been shown for humans that individuals may have Y-chromosomal haplotypes that are identical based on a large number of microsatellite markers even though the haplotypes as defined by biallelic markers are different, a contradiction which is due to recurrent mutations of the microsatellites [7]. This shows that Y-chromosomal biallelic markers are imperative for phylogenetic studies where they can serve as a backbone for the Y-chromosome phylogeny upon which the more detailed microsatellite variation can be imposed. Furthermore, for the dating of ancient population genetic events, analysis of DNA sequence variation offers a better tool than do microsatellites. There is thus a need for DNA sequence information for the dog Y chromosome.

However, so far only 3200 bp of nonrepetetive Y-chromosomal dog sequence has been reported in GenBank [8,9]. With the recent sequencing of the dog genome [10] the sequence for the autosomal and the X chromosomes were identified but, since the sequence was from a female dog, the Y chromosome sequence was not obtained. The genome of a male dog has earlier been surveyed by the sequencing with 1.5× coverage of a male dog [11], but because of the low coverage of this genome sequence an assembly into the different chromosomes is not possible and a definite identification of Y chromosomal sequences is therefore not possible.

In order to identify Y-chromosomal sequence we therefore compared the two genome datasets by Blast analysis, identifying sequences present in the male genome and absent in the female genome, and also compared the sequences with human Y chromosome sequence for further screening. The candidate sequences were then checked for male specificity by PCR analysis of male and female dogs. In order to identify SNPs to be used for population genetic studies we thereafter sequenced 14437 bp of the identified Y-chromosome sequence in 10 individuals representing dogs from all major populations around the world, giving a phylogenetic backbone for future population genetic studies of the Y-chromosome variation in dogs.

Results and discussion

Identification of Y chromosomal sequence

In order to identify Y-chromosome sequences in the male dog genome sequence we performed a combination of Blast and PCR analysis. For a sequence to be accepted as Y specific three criteria had to be fulfilled: (i) it should, in a Blast search against the human genome, have a best hit or be linked to a read with a best hit, to the human Y chromosome, (ii) it should in a Blast search against the reference female dog genome sequence have a negative result and (iii) it should in a PCR against samples from five male and five female dogs yield PCR product only for the male samples.

Identification of candidate Y-chromosome sequences by Blast analysis

The Blast analysis was performed in two steps; to reject the thousands of male sequences that lack identity with the female genome sequence owing to gaps in the female genome assembly or polymorphisms between the two genomes, a screen of the male dog genome sequence against the human Y chromosome sequence was done before comparing with the female dog genome sequence. Assembled survey-sequence data, based on approximately two million contigs and singleton reads, from the genome of a male poodle [11] was compared to the reference human genome sequence using blastn. For 32 contigs, the highest-scoring alignments were with a single region on the human Y-chromosome (1; rows 1–32). Component reads from these contigs were linked physically to mate reads from the other end of the insert in the same clones and these were, in turn, often part of distinct contigs consisting of several other reads. These sequences were also considered as possible Y-chromosome sequence (1; rows 33–93) giving totally 93 potential Y-chromosome sequences which were further compared to the reference female boxer genome sequence [11]. Sequences sharing >95% nucleotide identity for more than 90% of their lengths with regions of the boxer genome sequence were discounted, resulting in 43 remaining sequences considered as potential Y-chromosome sequences (1; Boxer Hit:"-").

Additional file 1. Identification of Y-chromosome sequences; results from blast analysis and PCR screen.

a Name of shotgun sequence.

b GenBank accession number of shotgun sequence used as template for primer design.

c Best hit versus Human genome. Shown here is the highest-scoring alignment with a single region on the human genome.

d Female dog genome hit. A match versus the female dog genome (sequences sharing >95% nucleotide identity for more than 90% of their lengths with regions of the female dog genome sequence) is indicated with a "+" and a non-match is indicated with a "-".

e,f Result from the first (two males and two females) and second (three males and three females) PCR screen. "Candidate" indicates a sequence that gave a PCR fragment for male dogs only, "Not male specific" a sequence that gave a PCR fragment for both male and female dogs, a dash indicates a sequence that was not tested, "Approved" indicates a sequence that gave a PCR fragment for all five male dogs tested and no band for the five female dogs, and "Failed" indicates a sequence that did not give three male bands and/or gave female bands in the second test.

Format: XLS Size: 28KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Confirmation of male specificity by PCR analysis

The 43 sequences considered potential Y-chromosomal sequence were tested for male specificity by PCR in two steps, first against a panel of two male and two female dogs (DNA sample panel 1), and the sequences passing this screen were further tested against a panel of three male and three female dogs (DNA sample panel 2). To investigate the correctness of excluding sequences with match to the female genome sequence, even though they had a best hit to the Y chromosome in comparison to the human genome, 29 of the sequences excluded based on this criterion were also tested. Upon agarose gel electrophoresis, sequences giving bands or smears of the same size for both males and females were excluded from further analysis, as were sequences giving no bands at all. Sequences giving bands only for males were kept as candidates.

35 sequences, of the 43 considered potentially Y-chromosome specific, passed the first PCR screen (1). Of the eight sequences which did not pass, seven gave bands or smears for both males and females and one gave no bands at all. Of the 29 sequences excluded based on match to the female genome all failed, 28 giving both male and female bands and one giving no bands at all. The 35 sequences that passed the first screen were selected for the second PCR screen against three more male and female dogs to confirm the male specificity. 32 of the 35 candidates were approved (giving bands of correct size for all five males, and no band at all for all five females) and thus considered male specific. Of the three failed sequences two had a band for only one of the male dogs and the third one had a weak band for one female dog. The GenBank Accession number for these sequences, totally 24159 bp in length, can be found in 1.

Sequence variation

In order to study the genetic variation for the Y-chromosome among dogs, and to identify nucleotide diversity to be used in population genetic studies, we performed DNA sequencing of 10 male dogs for 24 of the sequences identified as male specific, for a total of 14437 bp. Twelve of these sequences formed six pairs obtained from the same shotgun clone, with intervening unknown sequence. This sequence was determined through primer walking. The sequence data generated for the ten samples, including the here defined intervening sequence, have been deposited in GenBank, under accession numbers DQ973626DQ973805. In order to avoid geographical bias for the identified polymorphisms and to cover a large amount of the genetic variation among dogs globally, the 10 dogs were chosen to represent several different populations and types of domestic dogs from all continents except Australia (DNA sample panel 3).

Among the 10 dogs, 9 haplotypes were found which were defined by a total of 14 polymorphic positions, 13 substitutions and 1 indel (Figure 1). As expected for haploid sequence, none of the individuals were heterozygous for any nucleotide position. The average number of differences between two individuals was 4.2, giving a π value of 2.91 × 10-4. A π value for the wolf Y chromosome of 0.4 × 10-4 has earlier been reported [12], seven times lower than the here reported value for dogs. The reason for this discrepancy is not clear. One explanation might be that the here analysed dog sequences probably mostly derive from non-genic regions while the wolf sequences were mostly intronic sequence, possibly giving different substitution rates for the two data sets. Another explanation might be a geographical bias for the wolf sample (30/36 animals were from Russia (unspecified region)/the Baltics/Scandinavia, which because of the great mobility of wolves may be quite closely related populations). It can be noted that for the human Y chromosome π = 1.51 × 10-4 (2.3 Mb of the NRY sequenced) [13], which is in best agreement with the value for the dogs in this study.

thumbnailFigure 1. Sequence alignment showing base substitutions and indels for the 10 analysed samples. Only variable sites are shown, with sequence positions given above. Identity with the reference sequence (Yorkshire Terrier) is denoted by a stop, substitution by a different base letter, and deletion by a dash. Positions are denoted relative to the sequences submitted to GenBank (DQ973626DQ973805).

Phylogenetic relationships

The fact that analysis of ten individuals resulted in nine haplotypes indicates that there are numerous more Y-chromosome haplotypes to be found in the global dog population. A phylogenetic tree displays the relationship between the nine haplotypes (Figure 2). There is no clear geographic pattern for the distribution of haplotypes; possibly the sample is too small to show any such trends. The nucleotide substitution rate for the human Y chromosome is approximately 1.65 × 10-9 per bp and year [14], and it has been shown that the average genome mutation rate is approximately constant among mammals [15]. Assuming a mutation rate of 1.65 × 10-9 substitutions per year and bp also for the dog Y chromosome, in the lack of a species specific rate, we get a rate of one substitution per 42000 years for the analysed region. Assuming that the domestic dog originated around 15,000 years ago [2,16] this suggests that the nine dog Y-chromosomal haplotypes, several of which are more than one mutation step from its nearest neighbours, would originate from at least five wolf haplotypes. This agrees with an origin from a minimum of six female wolves according to mtDNA studies [2] and of totally a minimum of 21 individuals according to a study of dog MHC-variation [17].

thumbnailFigure 2. Maximum parsimony phylogenetic tree of the nine haplotypes found in the sequence analysis. Sequence types (circles) and empty nodes (solid dots) are separated by one mutational step (substitutions; indels are not shown).

Conclusion

We have in this study, using Blast analysis and PCR screening among male and female dogs, identified 24159 bp of dog Y-chromosome sequence to be used for population genetic studies. We sequenced 14437 bp of this sequence in a worldwide collection of 10 dogs, giving a first glimpse into the dog Y-chromosome phylogeny. Nine haplotypes were found among the ten dogs, indicating that there are numerous more Y-chromosome haplotypes to be found in the global dog population. The 14 identified SNPs can be used for SNP analyses, and the SNP based phylogeny will serve as a backbone for studies of Y chromosome microsatellite variation.

Methods

Blast analysis of genome shotgun sequences

Assembled survey-sequence data from the genome of a male poodle [11] was searched against the reference human genome sequence [18] to identify sequences with a best match to the Y chromosome, using NCBI-blastn with the following parameters; -e 1e-3 -q -11 -r 10 -G 20 -E 10 -X 150 -F T -b 1 -v 1 -m 8. Sequences with several matches were excluded as repetitive sequences. The male poodle sequences were also searched against the reference female boxer genome sequence [19] using NCBI-blastn with the following parameters; -W 40 -b 1 -v 1 -m 8. A match was defined as an alignment with >95% nucleotide identity for >90% of the length of the poodle sequence, and sequences without a match were considered to be potential Y-chromosome sequences (1).

PCR analysis

DNA samples

Samples used for the first screening analysis for male specificity (Panel 1, blood samples): Males; German Sheppard (region of origin: Europe) and Doberman Pincher (Europe), Females; Toy Poodle (Europe) and Springer Spaniel (Europe). Samples used for the second screening analysis for male specificity (Panel 2, buccal epithelial cell samples): Males; Leonberger (Europe), Tibetan Spaniel (East Asia) and Samoyed (Siberia), Females; Leonberger (Europe), Chow-Chow (East Asia) and Thai Ridgeback (Southeast Asia). Samples used for sequence analysis (Panel 3, all males and buccal epithelial cell samples): Yorkshire Terrier (Europe), Belgian Groenendael (Europe), Caanan dog (Southwest Asia), Afghan Hound (Southwest Asia), East Siberian Laika (Siberia), Tibetan Spaniel (East Asia), Thai Ridgeback (Southeast Asia), Akita Inu (Japan), Basenji (Africa) and Mexican Hairless Dog/Xoloitzquintle (America).

DNA extraction

Blood samples were treated with heparin or EDTA. The DNA extractions from blood were performed using the protocol #PT 3628-1, version #PR 22673, for the "NucleoSpin Blood Kit" (Biosciences Clontech). Prior to PCR, extracted DNA from the blood samples was diluted by a factor of eight. Buccal epithelial cell samples were collected using Whatman FTA-indicating cards according to the manufacturer's specifications [20].

DNA amplification for test of male specificity

PCR amplification of target sequences was performed using a forward and reverse primer (2). The PCR mixture consisted of 2 mM MgCl2 (Invitrogen), 20 mM TRIS-HCl pH 8.4 (Invitrogen), 50 mM KCl (Invitrogen), 200 mM of each dNTP (GE Healthcare), 1 unit PlatinumTaq DNA Polymerase (Invitrogen) and 200 mM of each primer in a total volume of 50 μl. The DNA template was 1 μl DNA extract from blood samples or a ∅ 2 mm piece from an FTA indication card. The reaction was run in a ThermoHybaid MBS 0.2 S (Thermo Electron Corporation) and the PCR program consisted of a predenaturation step (94°C, 2 min), followed by 40 cycles of denaturation (94°C, 30 sec), primer annealing (55°C, 30 sec) and extension (72°C, 3 min) followed by a final extension step (72°C, 10 min). The presence of amplification products was checked by analyzing five μl PCR product with agarose gel electrophoresis and staining with ethidium bromide.

Additional file 2. Primers used for PCR and sequence analysis and positions analysed in DNA sequencing.

a Name of analysed sequence. GenBank accession number given in 1. The shotgun sequencing approach yields two sequences for each clone, one from each end of the insert. Such a pair can be analyzed as a single unit in the PCR screen (if the intervening sequence is short enough to allow a PCR). Twenty sequences, among the 72 selected for the first PCR screening, were such paired sequences and they are listed together and the upstream sequence is listed first. For these ten pairs one PCR primer was placed in each sequence, giving amplification over both sequences. For the 24 sequences that were sequenced in ten dogs, the name used in GenBank is given in parentheses. "$" indicates sequences that here are treated reverse complementary compared to the shotgun sequence found in GenBank.

b PCR primers: P = PCR primer used in PCR screens; O = Outer primer and I = Inner primer in nested amplification; f = forward primer and r = reverse primer.

c Sequencing primers: f = forward primer and r = reverse primer in sequence analysis.

d Analysed positions. Here the region of the original shotgun sequence that was analysed is given. In four cases two separate parts of the shotgun sequence was analysed, because of intervening sequence giving bad quality reads.

Format: XLS Size: 71KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

DNA amplification for sequence analysis

PCR amplification of target sequences was performed in a nested configuration to improve specificity of the amplification (2). Both outer and inner PCRs had PCR mixtures and PCR programs identical to the amplification for test of male specificity, with the exceptions that the template for the inner PCR was 0.5 μl of PCR product from the outer PCR and the number of cycles for the outer and inner PCR was, respectively, 15 and 35.

DNA sequence analysis

All nucleotide positions were sequenced by at least one forward and one reverse sequence read, using the primers given in 2. For the cycle sequencing reaction, 1 μl of the inner amplification product was mixed with 13 μl 1× cycle sequencing buffer (26 mM Tris, pH 9.0, 6.5 mM MgCl2), 1 μl big dye terminator (ABI Prism Big Dye Terminator Cycle Sequencing Ready Reaction Kit v2.1, Applied Biosystems) and 0.25 μM primer, for a total reaction volume of 20 μl. The reaction was run in a ThermoHybaid MBS 02 S. The cycle sequencing program consisted of 30 cycles of denaturation (96°C, 10 s), primer annealing (55°C, 15 s) and extension (60°C, 4 min). The cycle sequencing products were ethanol precipitated and analysed on an ABI 3700 according to the manufacturers' directions (Applied Biosystems). The DNA sequences were edited using Sequencing Analysis 2.1.1 (Applied Biosystems) and assembled into contigs and further edited in Sequencher 4.1 (Gene Codes Corporation). Alignment and comparison of sequences were performed in Sequencher 4.1.

Authors' contributions

CN designed and carried out PCR screening, performed phylogenetic analysis and drafted the manuscript. MO designed sequence analysis and performed DNA sequencing. HA designed and carried out PCR screening and performed DNA sequencing. JL participated in co-ordinating and planning the study and helped drafting the manuscript. EK designed and performed Blast analysis and helped to draft the manuscript. PS conceived the study, and participated in its design and coordination and drafted the manuscript. All authors read and approved the final manuscript.

Acknowledgements

CN was supported by the Swedish Research Council and the Knut and Alice Wallenberg Foundation. MRCO was supported by the Swedish Research Council. HA was supported by VINNOVA (contract number PL4043-4). JL was supported by the faculty of the Royal Institute of Technology. EK was supported by The Institute for Genomic Research. PS was supported by the Swedish Research Council. This work was supported by the Swedish Kennel Club, the Magn. Bergvall Foundation, the Carl Trygger Foundation and the Swedish Research Council.

References

  1. Vila C, Savolainen P, Maldonado JE, Amorim IR, Rice JE, Honeycutt RL, Crandall KA, Lundeberg J, Wayne RK: Multiple and Ancient Origins of the Domestic Dog.

    Science 1997, 276:1687-1689. PubMed Abstract | Publisher Full Text OpenURL

  2. Savolainen P, Zhang YP, Luo J, Lundeberg J, Leitner T: Genetic evidence for an East Asian origin of domestic dogs.

    Science 2002, 298:1610-1613. PubMed Abstract | Publisher Full Text OpenURL

  3. 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.

    PNAS 2004, 101:12387-12390. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  4. Leonard JA, Wayne RK, Wheeler J, Valadez R, Guillen S, Vila C: Ancient DNA Evidence for Old World Origin of New World Dogs.

    Science 2002, 298:1613-1616. PubMed Abstract | Publisher Full Text OpenURL

  5. Bannasch DL, Bannasch MJ, Ryun JR, Famula TR, Pedersen NC: Y chromosome haplotype analysis in purebred dogs.

    Mammalian genome 2005, 16:273-280. PubMed Abstract | Publisher Full Text OpenURL

  6. Sundqvist AK, Björnerfeldt S, Leonard JA, Hailer F, Hedhammar Å, Ellegren H, Vila C: Unequal Contribution of Sexes in the Origin of Dog Breeds.

    Genetics 2006, 172:1121-1128. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  7. Jobling MA, Tyler-Smith C: New uses for new haplotypes the human Y chromosome, disease and selection.

    Trends in Genetics 2000, 16:356-362. PubMed Abstract | Publisher Full Text OpenURL

  8. Olivier M, Lust G: Two DNA sequences specific for the canine Y chromosome.

    Animal Genetics 1998, 29:146-149. PubMed Abstract | Publisher Full Text OpenURL

  9. Pujar S, Kothapalli KS, Kirkness E, Van Wormer RH, Meyers-Wallen VN: Exclusion of Lhx9 as a Candidate Gene for Sry-Negative XX Sex Reversal in the American Cocker Spaniel Model.

    Journal of Heredity 2005, 96:452-454. PubMed Abstract | Publisher Full Text OpenURL

  10. Lindblad-Toh K, Wade CM, Mikkelsen TS, Karlsson EK, Jaffe DB, Kamal M, Clamp M, Chang JL, Kulbokas EJ 3rd, Zody MC, Mauceli E, Xie X, Breen M, Wayne RK, Ostrander EA, Ponting CP, Galibert F, Smith DR, DeJong PJ, Kirkness E, Alvarez P, Biagi T, Brockman W, Butler J, Chin CW, Cook A, Cuff J, Daly MJ, DeCaprio D, Gnerre S, Grabherr M, Kellis M, Kleber M, Bardeleben C, Goodstadt L, Heger A, Hitte C, Kim L, Koepfli KP, Parker HG, Pollinger JP, Searle SM, Sutter NB, Thomas R, Webber C, Baldwin J, Abebe A, Abouelleil A, Aftuck L, Ait-Zahra M, Aldredge T, Allen N, An P, Anderson S, Antoine C, Arachchi H, Aslam A, Ayotte L, Bachantsang P, Barry A, Bayul T, Benamara M, Berlin A, Bessette D, Blitshteyn B, Bloom T, Blye J, Boguslavskiy L, Bonnet C, Boukhgalter B, Brown A, Cahill P, Calixte N, Camarata J, Cheshatsang Y, Chu J, Citroen M, Collymore A, Cooke P, Dawoe T, Daza R, Decktor K, DeGray S, Dhargay N, Dooley K, Dooley K, Dorje P, Dorjee K, Dorris L, Duffey N, Dupes A, Egbiremolen O, Elong R, Falk J, Farina A, Faro S, Ferguson D, Ferreira P, Fisher S, FitzGerald M, Foley K, Foley C, Franke A, Friedrich D, Gage D, Garber M, Gearin G, Giannoukos G, Goode T, Goyette A, Graham J, Grandbois E, Gyaltsen K, Hafez N, Hagopian D, Hagos B, Hall J, Healy C, Hegarty R, Honan T, Horn A, Houde N, Hughes L, Hunnicutt L, Husby M, Jester B, Jones C, Kamat A, Kanga B, Kells C, Khazanovich D, Kieu AC, Kisner P, Kumar M, Lance K, Landers T, Lara M, Lee W, Leger JP, Lennon N, Leuper L, LeVine S, Liu J, Liu X, Lokyitsang Y, Lokyitsang T, Lui A, Macdonald J, Major J, Marabella R, Maru K, Matthews C, McDonough S, Mehta T, Meldrim J, Melnikov A, Meneus L, Mihalev A, Mihova T, Miller K, Mittelman R, Mlenga V, Mulrain L, Munson G, Navidi A, Naylor J, Nguyen T, Nguyen N, Nguyen C, Nguyen T, Nicol R, Norbu N, Norbu C, Novod N, Nyima T, Olandt P, O'Neill B, O'Neill K, Osman S, Oyono L, Patti C, Perrin D, Phunkhang P, Pierre F, Priest M, Rachupka A, Raghuraman S, Rameau R, Ray V, Raymond C, Rege F, Rise C, Rogers J, Rogov P, Sahalie J, Settipalli S, Sharpe T, Shea T, Sheehan M, Sherpa N, Shi J, Shih D, Sloan J, Smith C, Sparrow T, Stalker J, Stange-Thomann N, Stavropoulos S, Stone C, Stone S, Sykes S, Tchuinga P, Tenzing P, Tesfaye S, Thoulutsang D, Thoulutsang Y, Topham K, Topping I, Tsamla T, Vassiliev H, Venkataraman V, Vo A, Wangchuk T, Wangdi T, Weiand M, Wilkinson J, Wilson A, Yadav S, Yang S, Yang X, Young G, Yu Q, Zainoun J, Zembek L, Zimmer A, Lander ES: Genome sequence, comparative analysis and haplotype structure of the domestic dog.

    Nature 2005, 438:803-819. PubMed Abstract | Publisher Full Text OpenURL

  11. Kirkness E, Bafna V, Halpern AL, Levy S, Remington K, Rusch DB, Delcher AL, Pop M, Wang W, Fraser CM, Venter JC: The Dog Genome: Survey Sequencing and Comparative Analysis.

    Science 2003, 301:1898-1903. PubMed Abstract | Publisher Full Text OpenURL

  12. Hellborg L, Ellegren H: Low Levels of Nucleotide Diversity in Mammalian Y Chromosomes.

    Mol Biol Evol 2004, 21:158-163. PubMed Abstract | Publisher Full Text OpenURL

  13. Sachidanandam R, Weissman D, Schmidt SC, Kakol JM, Stein LD, Marth G, Sherry S, Mullikin JC, Mortimore BJ, Willey DL, Hunt SE, Cole CG, Coggill PC, Rice CM, Ning Z, Rogers J, Bentley DR, Kwok PY, Mardis ER, Yeh RT, Schultz B, Cook L, Davenport R, Dante M, Fulton L, Hillier L, Waterston RH, McPherson JD, Gilman B, Schaffner S, Van Etten WJ, Reich D, Higgins J, Daly MJ, Blumenstiel B, Baldwin J, Stange-Thomann N, Zody MC, Linton L, Lander ES, Altshuler D: A map of human genome sequence variation containing 1.42 million single nucleotide polymorphisms.

    Nature 2001, 409:928-933. PubMed Abstract | Publisher Full Text OpenURL

  14. Schaffner SF: The X chromosome in population genetics.

    Nat Rev Genet 2004, 5:43-51. PubMed Abstract | Publisher Full Text OpenURL

  15. Kumar S, Subramanian S: Mutation rates in mammalian genomes.

    PNAS 2002, 99:803-808. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  16. Clutton-Brock J: Origins of the dog: domestication and early history. In The domestic dog: its Evolution, Behaviour and Interactions with People. Edited by Serpell J. UK: Cambridge University Press; 1995:7-20. OpenURL

  17. Vila C, Seddon J, Ellegren H: Genes of domestic mammals augmented by backcrossing with wild ancestors.

    Trends in Genetics 2005, 21:214-218. PubMed Abstract | Publisher Full Text OpenURL

  18. UCSC Genome Bioinformatics, Human Genome, Jul. 2003 (hg16) [http://hgdownload.cse.ucsc.edu/goldenPath/hg16] webcite

  19. UCSC Genome Bioinformatics, Dog Genome, Jul. 2004 (canFam1) [http://hgdownload.cse.ucsc.edu/goldenPath/canFam1] webcite

  20. Whatman FTA protocol BD08 [http://www.whatman.com/repository/documents/s7/51613.pdf] webcite