Dog Y chromosomal DNA sequence: identification, sequencing and SNP discovery
© Natanaelsson et al; licensee BioMed Central Ltd. 2006
Received: 10 July 2006
Accepted: 06 October 2006
Published: 06 October 2006
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.
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.
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.
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 . 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  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 , 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  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 (Additional file 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 (Additional file 1; rows 33–93) giving totally 93 potential Y-chromosome sequences which were further compared to the reference female boxer genome sequence . 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 (Additional file 1; Boxer Hit:"-").
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 (Additional file 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 Additional file 1.
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 DQ973626–DQ973805. 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).
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.
Blast analysis of genome shotgun sequences
Assembled survey-sequence data from the genome of a male poodle  was searched against the reference human genome sequence  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  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 (Additional file 1).
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).
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 .
DNA amplification for test of male specificity
PCR amplification of target sequences was performed using a forward and reverse primer (Additional file 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.
DNA amplification for sequence analysis
PCR amplification of target sequences was performed in a nested configuration to improve specificity of the amplification (Additional file 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 Additional file 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.
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.
- 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. 10.1126/science.276.5319.1687.PubMedView ArticleGoogle Scholar
- 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. 10.1126/science.1073906.PubMedView ArticleGoogle Scholar
- 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. 10.1073/pnas.0401814101.PubMedPubMed CentralView ArticleGoogle Scholar
- 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. 10.1126/science.1076980.PubMedView ArticleGoogle Scholar
- Bannasch DL, Bannasch MJ, Ryun JR, Famula TR, Pedersen NC: Y chromosome haplotype analysis in purebred dogs. Mammalian genome. 2005, 16: 273-280. 10.1007/s00335-004-2435-8.PubMedView ArticleGoogle Scholar
- 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. 10.1534/genetics.105.042358.PubMedPubMed CentralView ArticleGoogle Scholar
- Jobling MA, Tyler-Smith C: New uses for new haplotypes the human Y chromosome, disease and selection. Trends in Genetics. 2000, 16: 356-362. 10.1016/S0168-9525(00)02057-6.PubMedView ArticleGoogle Scholar
- Olivier M, Lust G: Two DNA sequences specific for the canine Y chromosome. Animal Genetics. 1998, 29: 146-149. 10.1046/j.1365-2052.1998.00299.x.PubMedView ArticleGoogle Scholar
- 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. 10.1093/jhered/esi058.PubMedView ArticleGoogle Scholar
- Lindblad-Toh K, Wade CM, Mikkelsen TS, Karlsson EK, Jaffe DB, Kamal M, Clamp M, Chang JL, Kulbokas EJ, 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. 10.1038/nature04338.PubMedView ArticleGoogle Scholar
- 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. 10.1126/science.1086432.PubMedView ArticleGoogle Scholar
- Hellborg L, Ellegren H: Low Levels of Nucleotide Diversity in Mammalian Y Chromosomes. Mol Biol Evol. 2004, 21: 158-163. 10.1093/molbev/msh008.PubMedView ArticleGoogle Scholar
- 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. 10.1038/35057149.PubMedView ArticleGoogle Scholar
- Schaffner SF: The X chromosome in population genetics. Nat Rev Genet. 2004, 5: 43-51. 10.1038/nrg1247.PubMedView ArticleGoogle Scholar
- Kumar S, Subramanian S: Mutation rates in mammalian genomes. PNAS. 2002, 99: 803-808. 10.1073/pnas.022629899.PubMedPubMed CentralView ArticleGoogle Scholar
- Clutton-Brock J: Origins of the dog: domestication and early history. The domestic dog: its Evolution, Behaviour and Interactions with People. Edited by: Serpell J. 1995, UK: Cambridge University Press, 7-20.Google Scholar
- Vila C, Seddon J, Ellegren H: Genes of domestic mammals augmented by backcrossing with wild ancestors. Trends in Genetics. 2005, 21: 214-218. 10.1016/j.tig.2005.02.004.PubMedView ArticleGoogle Scholar
- UCSC Genome Bioinformatics, Human Genome, Jul. 2003 (hg16). [http://hgdownload.cse.ucsc.edu/goldenPath/hg16]
- UCSC Genome Bioinformatics, Dog Genome, Jul. 2004 (canFam1). [http://hgdownload.cse.ucsc.edu/goldenPath/canFam1]
- Whatman FTA protocol BD08. [http://www.whatman.com/repository/documents/s7/51613.pdf]