Journal of Heredity 2004:95(2):158-164
© 2004 The American Genetic Association
Brief Communication |
Isolation of Y Chromosome-specific Microsatellites in the Horse and Cross-species Amplification in the Genus Equus
From the Institut für Tierzucht und Genetik, Veterinärmedizinische Universität Wien Veterinärplatz 1, A-1210 Wien, Austria (Wallner, Brem, and Müller); INRA, Centre de Recherches de Jouy, Departement de Genetique Animale, Jouy-en-Josas, France (Piumi); and Ludwig Boltzmann-Institut für immuno-, zyto- und molekulargenetische Forschung, Veterinärplatz 1, A-1210 Wien, Austria (Brem and Achmann). We thank W. Zimmermann, C. Waltzer, H. Korb, E. Dvorak, H. Siegismund, the Zoological Society of San Diego, and the Zoological Garden Munich for samples. We are grateful to M. Wigler for a detailed RDA protocol. Thanks to G. Muir for comments on the manuscript. This study was financially supported by the Hochschuljubiläumsstiftung der Stadt Wien (H-234/99) and Agrobiogen GmbH.
Address correspondence to Roland Achmann, genteQ GmbH, Falkenried 88, D-20251 Hamburg, Germany, or e-mail: achmann{at}genteQ.de.
| Abstract |
|---|
|
|
|---|
Y chromosome polymorphisms such as microsatellites or single nucleotide polymorphisms represent a paternal counterpart to mitochondrial DNA (mtDNA) for evolutionary and phylogeographic studies. The use of Y chromosome haplotyping in natural populations of species other than humans is still hindered by the lack of sequence information necessary for polymorphism screening. Here we used representational difference analysis (RDA) followed by a screen of a bacterial artificial chromosome (BAC) library for repetitive sequences to obtain polymorphic Y-chromosomal markers. The procedure was performed for the domestic horse (Equus caballus) and we report the first six Y-chromosomal microsatellite markers for this species. Three markers were also useful for haplotyping taxa of the zebra/ass lineage. Y-chromosomal microsatellite markers show a single haplotype in the domestic horse, whereas notable variation has been observed in the other members of the genus Equus.
Recent advances in genome analysis and population genetics in humans suggest that Y chromosome haplotyping may develop into an important tool for studying natural populations (Hurles and Jobling 2001). With the exception of the pseudoautosomal region, the entire Y chromosome acts as a single nonrecombining unit, which is male specific and effectively haploid. This ensures that original combinations of mutational events along male lineages are preserved as single-linked unambiguous haplotypes. Y-chromosomal markers are required as a paternal analogue to mitochondrial DNA (mtDNA). Polymorphic markers on the nonrecombining portion of the Y chromosome range from rarely occurring biallelic point mutations (single nucleotide polymorphisms [SNPs]) to faster-mutating multiallelic polymorphisms found at minisatellite or microsatellite (short tandem repeat [STR]) loci. Y-chromosomal SNPs were often found to be population specific (Hammer et al. 1997). Therefore, for useful measures of diversity of recently diverged populations (e.g., breeds in domestic animals), more polymorphic microsatellite markers are necessary.
In contrast to sequence data on human and mouse, Y chromosome information is scarce for other mammals (Hurles and Jobling 2001). For equids, where mtDNA and autosomal microsatellites have been used to investigate the phylogenetic relationship within the genus (Oakenfull et al. 2000) and horse domestication (e.g., Jansen et al. 2002; Vila et al. 2001), polymorphic Y-chromosomal markers are not available at present. The construction of a haploid chromosome-specific library requires the preselection of a particular chromosome by flow cytometry (Bergstrom et al. 1998) or microdissection (Ponce de Leon et al. 1996). Both methods are technically demanding and suffer from limitations such as separation efficiency and minute amounts of DNA as starting material (Claussen et al. 1997). Alternatively, haploid markers can be isolated by the use of sequence tags (STSs) as probes for selecting chromosome-specific clones from genomic libraries containing large inserts. The inserts of identified clones can be further screened for microsatellite markers.
Here we used representational difference analysis (RDA) (Lisitsyn et al. 1993) for the isolation of Y chromosome-specific STSs. RDA is a subtractive approach that permits the cloning of differences between similar genomes. The principle behind RDA is that unique target sequences present in one DNA population (tester, heterogametic XY) can be purified through repetitive steps of difference enrichment by hybridization with an excess of a second DNA population (driver, homogametic XX). RDA-enriched Y-specific STSs were used as probes for screening a bacterial artificial chromosome (BAC) library to subsequently isolate microsatellites by subcloning. We demonstrate the practicability of this strategy to isolate Y-specific STSs and microsatellite markers for the domestic horse (Equus caballus). Moreover, we present data on the variability of the first six equine Y-chromosomal microsatellite markers in the domestic horse and other equine species.
| Materials and Methods |
|---|
|
|
|---|
Representational Difference Analysis
The RDA procedure was applied essentially as previously described (Lisitsyn et al. 1993), using the restriction enzyme BglII and corresponding amplification adapter oligonucleotides. Genomic DNA from a Lipizzan stallion and Lipizzan mares was used as tester and driver, respectively. In order to avoid a significant reduction in the kinetic enrichment of the target due to the "Cot effect" (Mathieu-Daudé 1996), the number of polymerase chain reaction (PCR) cycles used to reamplify the product after each subtractive round was kept to a minimum (first round: 12 cycles; third round: 22 cycles). Fragments of sizes between 500 and 1000 bp of the amplicon pools after two and three rounds of subtraction were cloned into pUC19 (MBI Fermentas). Inserts from 30 clones, 7 from the library after two subtractive rounds and 23 after three subtractive rounds, were sequenced with BigDye terminator cycle sequencing chemistry (Applied Biosystems) using standard M13 primers. Sequences were classified by their similarity into different sequence families using Autoassembler software (Applied Biosystems). For each family a randomly selected sequence was designated as the reference sequence. Sequences were compared for homologies to known sequences by using BLASTN. Repetitive elements were identified by the program Repeat Masker (Smit and Green 2003; program available at http://ftp.genome.washington.edu/RM/RepeatMasker.html).
Primers were designed for each reference sequence (Table 1) and Y chromosome specificity was tested by PCR amplification of genomic DNA from males and females. Amplification was performed in 15 µl reactions containing 20 ng DNA, 0.5 µM of each primer, 0.25 mM of each dNTP, 2 mM MgCl2, 1x AmpliTaq buffer, and 0.5 U AmpliTaq Gold (Applied Biosystems). PCR conditions were 10 min at 95°C, followed by 35 cycles of 30 s at 95°C, 40 s at annealing temperature (see Table 1), 90 s at 72°C, and a final extension of 10 min at 72°C. Ethidium bromide-stained PCR products were visualized on 2% agarose gels. DNA samples (one per breed) from 14 different breeds (Akhal-Teke, Andalusian, Arabian, Austrian Warmblood, Icelandic Horse, Kladruber, Lipizzan, Mongolian Domestic Horse, Noric, Quarter Horse, Shetland Pony, Tarpan, Thoroughbred, Trakehner) were directly sequenced for male-specific STSs.
|
Isolation of Microsatellites from Selected BAC Clones
A horse BAC library was screened for Y-chromosomal clones as described (Godard et al. 1998) with primers specific for eight Y-chromosomal STSs isolated by RDA (we excluded locus SH3-B-14 because of difficulties in PCR amplification) and a horse SRY-specific primer pair (Hasegawa et al. 1999). We isolated BAC DNA as described at http://dga.jouy.inra.fr/grafra/ (last visited December 2003). BAC DNA was prescreened for the presence of (CA/GT)n microsatellites by dot-blot hybridization.
Microsatellite enriched libraries were established as described by Armour et al. (1994). The enrichment of (CA/GT)n-containing elements was performed as described in Traxler et al. (2000), with the modification that adapters RBgl24 and RBgl12 (Lisitsyn et al. 1993) were used for reamplification. Twenty randomly chosen microsatellite-enriched subclones originating from each BAC clone were subjected to colony PCR using standard M13 primers. Amplicons were controlled for the presence of (CA/GT)n sequences by dot-blot hybridization (see above). For positive clones, plasmid DNA was extracted and sequenced.
Microsatellite Genotyping
We designed specific PCR primers for identified microsatellite sequences and checked their Y-chromosomal location as described above. One primer of each pair was labeled with a fluorescent dye to allow simultaneous analysis of fragments on an ABI310 instrument (Applied Biosystems). PCR was performed in a 15 µl volume containing primers (for concentration see Table 2) and the PCR reaction mix as described above. The DNA was initially denatured at 95°C for 10 min followed by 35 cycles of 30 s at 95°C, 40 s at annealing temperature (see Table 2), 90 s at 72°C, and a final extension of 30 min at 72°C. Alleles were sized relative to the internal size standard Tamra GS 500 using Genescan version 2.1 (Applied Biosystems).
|
Microsatellite variability was evaluated with DNA from 49 male horses of 32 different breeds representing the domestic horse (Akhal-Teke, Andalusian, Appaloosa, Arabian, Austrian Warmblood, Barb, Connemara, Icelandic Horse, Irish Tinker, Kladruber, Lipizzan, Mangolarga Marchoder, Minishetty, Missourie Foxtrotter, Mongolian Domestic Horse, New Forest, Noric, Norwegian Fjord, Oldenburger, Paint, Pinto, Quarter Horse, Shagya Arabian, Shetland Pony, Saddlebred, Shire Horse, Tarpan, Thoroughbred, Trakehner, Trotter, Welsh Pony, Wurttemberger). In addition DNA samples from other equids were analyzed: 10 Przewalski's Horses (Equus przewalskii), 3 Donkeys (Equus asinus), 1 Onager (Equus hemionus onager), 3 Kiangs (Equus kiang), 3 Grevy's Zebras (Equus grevyi), 14 Hartmann's Zebras (Equus zebra hartmannae), 5 Damara Zebras (Equus burchelli antiquorum), and 2 Grant's Zebras (Equus burchelli boehmi).
| Results and Discussion |
|---|
|
|
|---|
Representational Difference Analysis
In order to isolate specific markers from the equine Y chromosome we performed RDA with male equine genomic DNA as tester and female genomic DNA as driver. Figure 1 shows amplicons after the first, second, and third round of RDA. The second and third round products show an enrichment of discrete bands and a reduction in the intensity of smearing background DNA. Second and third round amplicons were size fractionated and subcloned. Thirty randomly selected clones were sequenced and sequences clustered into 14 families. Eight of the 14 sequence families showed similarity to common genome-wide repeats and one clone showed homology to a human Y-chromosomal sequence (Table 1).
|
For 9 of the 14 reference sequences, PCR on male samples amplified a single band, whereas female samples could not be amplified. The exclusive amplification of male DNA strongly supports a Y-chromosomal origin of the isolated sequences (Table 1).
Because four of the nine Y-specific sequences show a high homology to known repetitive elements (LINEs or ERVs), we tested for multiple copies on the Y chromosome by direct sequencing of PCR products from 14 male horses from different breeds. Two STSs with partial homologies to endogenous retroviral sequences showed heterozygous sequencing signals at some nucleotide positions (Eca SH3-B-6 at 4 positions, Eca SH3-B-14 at 10 positions), indicating multiple, divergent copies on the horse Y chromosome. Because of the repetitive nature of the Y chromosome (Foote et al. 1992) and the preferential isolation of repeats due to kinetic enrichment in the amplification steps (e.g., Panaud et al. 2002), the isolation of Y-specific repeats is not surprising. Sequencing of the remaining seven Y-chromosomal sequences showed no variable positions over a total length of 3033 nucleotides among the 14 individuals from different breeds. This observation is not unexpected, because the diversity on the Y chromosome is known to be extremely low (Shen et al. 2000).
Screening of a Horse BAC Library
We screened a male horse BAC library (Godard et al. 1998) with primers specific for nine Y-chromosomal STSs and selected 15 BAC clones in total (Table 1). Three markers exclusively identified a single clone and four markers simultaneously amplified in two to four clones. STSs SH3-B-12 and SH3-B-19 gave positive signals for the same two clones, indicating that they are located very closely on the Y chromosome. Hybridization of a (CA/GT)n probe to BAC DNA dot blots identified 11 BAC clones containing simple tandem repeat sequences. The application of an enrichment method (Armour et al. 1994, Traxler et al. 2000) permitted a quick isolation of dinucleotide microsatellites from BACs. The screening of microsatellite-enriched BAC DNA (20 random subclones per original BAC clone) by dot blot hybridization revealed that out of 220 subclones tested for microsatellite sequences, 56 showed a positive signal. We sequenced the insert of all positive subclones and identified 11 different microsatellite loci. Six of eight loci were successfully amplified from horse DNA (Table 2) and proved to be male specific. Electropherograms of Y-specific markers displayed the typical peak pattern (main peak preceded by one or two stutter peaks) of dinucleotide microsatellites (Figure 2a). However, marker Eca.YH12 produced a peak pattern with three main peaks (Figure 2a). Eca.YH12 could not be amplified from female DNA and amplification from the corresponding BAC clone 68B8 revealed a single main peak (Figure 2b). We assume that this microsatellite has multiple copies on the horse Y chromosome. The duplication or triplication of a larger Y-chromosomal region, followed by a change in the number of repeats within the microsatellite is a frequently observed phenomenon in humans (Kayser et al. 2000).
|
We assessed the variability of the six equine Y-specific microsatellites in a sample of 49 male domestic horses. Although males from 32 distinct breeds from different continents were investigated, we detected only a single haplotype (see Table 3).
|
The Przewalski's Horse (n = 10) haplotypes differed from that of the domestic horse at two loci (Eca.YH12 and Eca.YA16; Table 3). In contrast to the domestic horse, we identified two haplotypes, which were distinguished from the domestic horse haplotype by the same two loci. Eca.YH12 displayed either one or two detectable alleles (Figure 2c,d), which can be explained by the fact that Eca.YH12 is not triplicated in the Przewalski's Horse. Alternatively, the locus could be triplicated, but shows identical alleles at two or three loci. Finally, mutations in the primer binding sites may affect the amplification of this locus, so that one or two copies are not amplified at all (nonamplifiable alleles; Achmann et al. 2001). Our data show that although E. przewalskii has undergone a severe bottleneck (Boyd and Houpt 1994), at least two stallions must have contributed to the present-day population.
We further extended the analysis to several other equine species (Table 3). Three markers (Eca.YP9, Eca.YM2, and Eca.YE1) could be amplified under less stringent PCR conditions (Table 2). For E. asinus, E. hemionus onager, E. kiang, E. zebra hartmannae, E. burchelli antiquorum, E. burchelli boehmi, and E. grevyi unique haplotypes were observed, which allow the discrimination between each species/subspecies (Table 3). We discovered interspecific haplotypic variation in E. asinus, E. grevyi, and E. zebra hartmannae (Table 3), although we analyzed only a small number of animals (n = 314) in these species.
The six newly characterized Y-specific microsatellites displayed no variability in the domestic horse. This is in contrast to the genetic variation normally observed for equine autosomal markers such as microsatellites (e.g., Vila et al. 2001) and polymorphic genes (Shubitowski et al. 2001) as well as for mtDNA (Jansen et al. 2001; Vila et al. 2001). Although some of the isolated Y-chromosomal microsatellite loci for the horse are imperfect and have a low total number of repeats, they show variability in other equine species. Thus a low mutation rate alone cannot account for the observed Y chromosome uniformity in the domestic horse.
We suggest that the observation of a single Y chromosome haplotype in the horse compared to the high levels of mtDNA variation (Jansen et al. 2002; Vila et. al 2001) may reflect a strong bias toward females in breeding and trade. The effective population size of the Y chromosome corresponds to the number of breeding males in a population. Intensive breeding strategies in which selected breeding stallions are used for many mares (Levine 1999), together with upgrading many breeds by crossing in Arabian and thoroughbred stallions, have probably resulted in the fixation of a single Y chromosome haplotype. Further investigations are necessary to confirm this hypothesis, because a global selective sweep, due to advantageous mutations in a Y-specific gene, could also have led to the fixation of a single haplotype. Whatever the explanation for the apparently monomorphic Y chromosome in domestic horses is, it seems, that today's horse Y chromosomes can be traced back to a recent common ancestor.
In summary, we show that RDA in combination with the screen of a BAC library is a quick and effective strategy to isolate Y-chromosomal microsatellite markers. Since BAC libraries are available for several domestic animals, our approach for the isolation of Y-specific markers can be easily performed in other species. An increasing number of Y-specific markers could help to extend our knowledge about sex-specific differences in ecology, behavior, or migration not only in humans, but also in many other mammalian species.
| Footnotes |
|---|
Corresponding Editor: Oliver A. Ryder
Received May 29, 2003
Accepted January 15, 2004
| References |
|---|
|
|
|---|
-
Achmann R, Huber T, Wallner B, Dovc P, Muller M, Brem G, 2001. Base substitutions in the sequences flanking microsatellite markers HMS3 and ASB2 interfere with parentage testing in the Lipizzan horse. Anim Genet. 32:52.[CrossRef][Web of Science][Medline]
Armour J, Neumann R, Gobey S, Jeffreys A, 1994. Isolation of human simple repeat loci by hybridization selection. Hum Mol Genet. 3:599-604.
Boyd L, Houpt K, 1994. Przewalski's horse, the history and biology of an endangered species. Albany: State University of New York Press.
Bergstrom D, Grieco D, Sonti M, Fawcett J, Bell-Prince C, Cram LS, Narayanswami S, Simpson EM, 1998. The mouse Y-chromosome: enrichment, sizing, and cloning by bivariate flow cytometry. Genomics. 48:304-313.[CrossRef][Web of Science][Medline]
Claussen U, Senger G, Chuboba I, 1997. The use of microdissected chromosomes in genome mapping. In: Genome mapping. A practical approach (Dear PH, ed). Oxford: IRL Press; 185197.
Foote S, Vollrath D, Hilton A, Page DC, 1992. The human Y chromosome: overlapping DNA clones spanning the euchromatic region. Science. 258:60-66.
Godard S, Schibler L, Oustry A, Cribiu E, Guerin G, 1998. Construction of a horse BAC library and cytogenetical assignment of 20 type I and type II markers. Mamm Genome. 9:633-637.[CrossRef][Web of Science][Medline]
Hasegawa T, Ishida M, Harigaya T, Sato F, Ishida N, Mukoyama H, 1999. Linear SRY transcript in equine testis. J Vet Med Sci. 61:97-100.[CrossRef][Web of Science][Medline]
Hammer M, Spurdle A, Karafet T, Bonner M, Wood ET, Novelletto A, Malaspina P, Mitchell RJ, Horai S, Jenkins T, Zegura SL, 1997. The geographic distribution of human Y chromosome variation. Genetics. 145:787-805.[Abstract]
Hurles M, Jobling M, 2001. Haploid chromosomes in molecular ecology: lessons from the human Y. Mol Ecol. 10:1599-1613.[CrossRef][Medline]
Jansen T, Forster P, Levine MA, Oelke H, Hurles M, Renfrew C, Weber J, Olek K, 2002. Mitochondrial DNA and the origins of the domestic horse. Proc Natl Acad Sci USA. 99:10905-10910.
Kayser M, Roewer L, Hedman M, Henke L, Henke J, Brauer S, Kruger C, Krawczak M, Nagy M, Dobosz T, Szibor R, de Knijff P, Stoneking M, Sajantila A, 2000. Characteristics and frequency of germline mutations at microsatellite loci from the human Y chromosome, as revealed by direct observation in father/son pairs. Am J Hum Genet. 66:1580-1588.[CrossRef][Web of Science][Medline]
Levine M, 1999. Botai and the origins of horse domestication. J Anthropol Archaeol. 18:29-78.
Lisitsyn N, Lisitsyn N, Wigler M, 1993. Cloning the differences between two complex genomes. Science. 259:946-951.[Abstract]
Oakenfull A, Lim A, Ryder O, 2000. A survey of equid mitochondrial DNA: implications for the evolution, genetic diversity and conservation of Equus. Conserv Genet. 1:341-355.[CrossRef]
Mathieu-Daudé F, Welsh J, Vogt T, McClelland M, 1996. DNA rehybridization during PCR: the "Cot effect" and its consequences. Nucleic Acids Res. 24:2080-2086.
Panaud O, Vitte C, Hivert J, Muzlak S, Talag J, Brar D, Sarr A, 2002. Characterization of transposable elements in the genome of rice (Oryza sativa L.) using representational difference analysis (RDA). Mol Genet Genomics. 268:113-121.[CrossRef][Web of Science][Medline]
Ponce de Leon F, Ambady S, Hawkins G, Kappes S, Bishop MD, Robl JM, Beattie CW, 1996. Development of a bovine X chromosome linkage group and painting probes to assess cattle, sheep, and goat X chromosome segment homologies. Proc Natl Acad Sci USA. 93:3450-3454.
Shen P, Wang F, Underhill P, Franco C, Yang W, Roxas A, Sung R, Lin AA, Hyman RW, Vollrath D, Davis RW, Cavalli-Sforza LL, Oefner PJ, 2000. Population genetic implications from sequence variation in four Y chromosome genes. Proc Natl Acad Sci USA. 97:7354-7359.
Shubitowski DM, Venta PJ, Douglass CL, Zhou RX, Ewart SL, 2001. Polymorphism identification within 50 equine gene-specific sequence tagged sites. Anim Genet. 32:78-88.[CrossRef][Web of Science][Medline]
Smit AFA, Green P, 2003. RepeatMasker (visited December 2003) http://ftp.genome.washington.edu/RM/RepeatMasker.html).
Traxler B, Brem G, Muller M, Achmann R, 2000. Polymorphic DNA microsatellites in the domestic pigeon, Columba livia var. domestica. Mol Ecol. 9:366-368.[CrossRef][Medline]
Vila C, Leonard J, Götherström A, Marklund S, Sandberg K, Liden K, Wayne RK, Ellegren H, 2001. Widespread origins of domestic horse lineages. Science. 291:474-477.
This article has been cited by other articles:
![]() |
S.-J. Luo, W. E. Johnson, V. A. David, M. Menotti-Raymond, R. Stanyon, Q. X. Cai, T. Beck, N. Yuhki, J. Pecon-Slattery, J. L. D. Smith, et al. Development of Y Chromosome Intraspecific Polymorphic Markers in the Felidae J. Hered., July 31, 2007; (2007) esm063v2. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. J. Royo, I. Alvarez, A. Beja-Pereira, A. Molina, I. Fernandez, J. Jordana, E. Gomez, J. P. Gutierrez, and F. Goyache The Origins of Iberian Horses Assessed via Mitochondrial DNA J. Hered., November 1, 2005; 96(6): 663 - 669. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Raudsepp, A. Santani, B. Wallner, S. R. Kata, C. Ren, H.-B. Zhang, J. E. Womack, L. C. Skow, and B. P. Chowdhary A detailed physical map of the horse Y chromosome PNAS, June 22, 2004; 101(25): 9321 - 9326. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||



