Print version ISSN 1415-4757
Genet. Mol. Biol. vol.25 no.1 São Paulo 2002
Evaluation of polymorphism in ten microsatellite loci in Uruguayan sheep flocks
Ivanna Tomasco, G. Wlasiuk and E.P. Lessa
Laboratorio de Evolución, Facultad de Ciencias, Universidad de la República, Montevideo, Uruguay.
Send correspondence to Ivanna Tomasco. Facultad de Ciencias, Universidad de la República, Iguá 4225, 11400 Montevideo, Uruguay. E-mail: email@example.com.
The allele frequencies of 10 microsatellite loci previously described for sheep as BM1314, BM6526, OarFCB128, OarHH64, OarCP20, OarHH47, OarFCB48, OarHH35, OarHH72 and BM2508 were estimated for the Uruguayan flocks. A representative sample of 101 individuals composed by the two predominant breeds (76% Corriedale and 24% Australian Merino) was used. The sample did not show a significant tendency towards substructuring, in spite of presenting some significantly different allele frequencies between races. The Corriedale sample presents three loci in which the presence of null alleles is possible. The markers were highly variable, showing between 7 and 15 alleles each. The Polymorphism Information Content Index ranged from 0.63 to 0.87 and the Exclusion Probability from 0.39 to 0.75 for a cumulative Exclusion Probability of 99.98%. These results suggest the effectiveness of this set of loci for testing genetic relatedness. This is the first report of microsatellite variation in Corriedale.
Key words: sheep, polymorphism, microsatellites.
Received: March 14, 2002; accepted: March 25, 2002.
For every program of breeding and selection, it is desirable to rely on a precise assessment of parentage to assist, for example, in the elimination of undesirable recessive alleles, in early selection for future breeding, or to ensure accurate pedigrees and registration. Microsatellites, segments of the nuclear genome composed of tandem repeat of short-sequence motifs, have become excellent candidate markers for this kind of studies (Queller et al., 1993), since they are numerous, highly variable and easy to score.
A large number of highly polymorphic microsatellites have been characterized and mapped in domestic animals, including sheep, cattle and other ruminants (de Gortari et al., 1997, 1998; Hayes et al., 1996; Jenkings et al., 1997), facilitating the use of these markers in parentage testing. For this kind of test to be acceptable, it must provide a high degree of certainty in the assignment of parentage, which can be evaluated with the Exclusion Probability (Weir, 1996). This requires a previous knowledge of allele frequencies in the population of interest and certainty about the assumptions of the model, like independence of the markers utilized and Hardy-Weinberg equilibrium.
The goals of this study were to analyze the polymorphism of a set of microsatellites previously described as polymorphic in sheep and to evaluate their usefulness for relatedness testing in Uruguayan sheep flocks.
MATERIAL AND METHODS
We analyzed ten microsatellite loci in a sample of the Uruguayan flocks considering only pure-breed individuals. The analysis were carried out separately for two subsamples: one of 77 Corriedale individuals from a fair, and another of 24 Australian Merinos from an experimental field (all of them are expected to be unrelated ewes), as well as on the combined dataset. Corriedale, the predominant race in the country, has not been previously characterized using microsatellite markers.
Total DNA extractions were made with SDS/proteinase K/NaCl/alcohol precipitation (modified from Miller et al., 1988) from skin preserved in 95% ethyl-alcohol. The microsatellite loci OarHH35, OarHH64, OarHH47 (Henry et al., 1993), OarVH72 (Pierson et al., 1993), OarFCB128 (Buchanan and Crawford 1993), OarFCB48 (Buchanan et al., 1994), OarCP20, BM1314, BM6526, (Crawford et al., 1995), and BM2508 (Mulsant et al., 1998) were amplified by PCR (Polymerase Chain Reaction) for all animals; excluding OarHH35 and OarHH64 all of them map on different chromosomes. Amplification was carried out in a total volume of 12.5 mL containing the following constituents: 6.25 mL of DNA (0.4 mg/mL) used as a template, 1X Taq Polymerase Buffer, 240 mM of each dNTP, 240 nM of each primer, 0.5 units of Taq Polymerase and variable concentration of MgCl2 depending on the locus (from 1.5 to 3.0 mM). PCR amplifications were performed in a RapidcyclerTM (Idaho Technology) by an initial denaturation of 3 min at 94 °C, followed by 35 cycles of 45 s of denaturation at 94 °C, 45 s of annealing at 50 °C and 45 s of extension at 72 °C, and a final extension of 7 min at 72 °C.
The amplified products were electrophoresed in 7%, 0.80 mm thick denaturing polyacrylamide gels (1000-1200 V, 1-2 h), and the DNA bands were visualized by silver staining (Sanguinetti et al., 1994). Genotypes were scored by at least two persons independently and the absolute sizes of the alleles were determined in relation to a 10bp DNA size standard (GIBCO BRL).
On the basis of allele and genotype frequencies, Hardy-Weinberg equilibrium and linkage disequilibrium were tested using the software GenePop, version 3.1b (Raymond and Rousset, 1995). A Bonferroni correction was used to adjust significance levels in all tests to an overall a = 0.05. Observed heterozygosity (Hobs), Polymorphism Information Content (PIC, Botstein et al., 1980) and Exclusion Probability assuming that the genotype of the lamb and one of its parents is known (Q, according to Weir, 1996) were calculated.
The software Structure (Pritchard et al., 2000) was used in order to assess levels of subdivision between two races (10,000 Burn-in period, 10,000 MCMC repetitions, and default parameters).
Allele length in base pairs, allele frequencies, Hobs, PIC, Q and the number of individuals analyzed per locus are given in Table I. All loci were highly polymorphic, showing between 7 and 15 alleles, PICs between 0.63 and 0.87 and Hobs in the range of 0.62 to 0.84 for the whole sample. The two loci that map on the same chromosome (OarHH35 and OarHH64) did not show a significant linkage (p > 0.373, SE = 0.052).
The loci OarHH35, OarHH64 and OarFCB128 showed significant departures from Hardy-Weinberg equilibrium in the whole population and in the Corriedale sample (Table II). With the exception of OarFCB128, OarVH72 and BM6526, all loci presented significant differences in allele frequencies between Merino and Corriedale samples. Some alleles were much more frequent in one race than in the other (e.g., allele 123 in OarHH35 and 146 in OarHH47, which are 14 and 8 times more frequent, respectively in Merino than in Corriedale, or allele 75 in OarCP20 which is seven times more frequent in Corriedale than in Merino). Other alleles were found only in one sample, and sometimes at high frequency (e.g., alleles 134 and 142 of OarHH47 in Corriedale, and allele 85 of OarCP20 in Merino). In several cases, these frequencies differed from those reported in the literature (see Henry et al., 1993, for OarHH47, OarHH35, OarHH64; Buchanan and Crawford 1993, for OarFCB128; Ede et al., 1995, for OarCP20; and Buchanan et al., 1994, for OarFCB48). There were also differences in one base pair between absolute allele sizes reported for OarFCB48 and OarFCB128 and our observations but it was easy to resolve the correspondences on the basis of allele frequencies. Direct comparisons would be needed to assess the reason for this discrepancy.
Analyses of subdivision using the software Structure (Pritchard et al., 2000) did not find a significant tendency towards substructuring in the sample (data not shown), in spite of the presence of two breeds.
The loci used in the analysis are effectively independent, since the only two that map on the same chromosome (OarHH35 and OarHH64) did not show significant linkage disequilibrium. The departures from Hardy-Weinberg equilibrium shown by OarHH35, OarHH64 and OarFCB128 in the Corriedale sample were probably the cause of the departure in the whole population. This could be due to the presence of null alleles in these loci in Corriedale or be the consequence of several years of intensive selection on these breeds (Usha et al., 1995). To avoid this problem, these loci can be excluded from parentage analysis and the cumulative Q would vary from 99.99% to 99.77% and from 99.98% to 99.75% for the whole population and the Corriedale stock, respectively. For Merinos, the cumulative Q with ten loci was 99.93%.
In spite of the presence of several differences in allele frequencies between races, these were insufficient to result in substantial separation between these breeds. These two races are known to be related, since Corriedale was originated from crosses between pure Merino and Lincoln sheep (Majala, 1997), so this result is not entirely surprising. Coupled with the fact that most of Uruguayan sheep are not completely pure, these observations suggest that allele frequency estimations derived from the whole sample may be sufficient for future applications in relatedness tests. Also, such tests may exclude the loci suspected to have null alleles without a marked reduction in the probability of exclusion.
The fact that estimated frequencies differed from those reported in the literature in most loci analyzed, stresses the importance of independent estimates for different regions and races.
In sum, these results indicate that the selected system of markers is highly effective for relatedness studies, with a Probability of Exclusion greater than 99.9% considering all loci. The significance of this report, the first of this kind in the country, is that it offers interesting perspectives for the incorporation of molecular genetic techniques to animal breeding in Uruguay. In addition, our results represent an original contribution of allele frequencies for regional Corriedale and Australian Merino races.
We are grateful to Jaime Mendoza, Roberto Cardellino and Mario Azzarini, from SUL (Secretariado Uruguayo de la Lana) for encouragement and collaboration in the field work, and to CSIC (Comisión Sectorial de Investigación Científica de la Universidad de la República) for financial support of this work.
Botstein D, White RL, Skolnick M and David R (1980) Construction of a genetic linkage map in man using restriction fragment length polymorphisms. Am. J. Hum. Genet. 32:314-331. [ Links ]
Buchanan F and Crawford AM (1993) Ovine microsatellites at the OarFCB11, OarFCB128, OarFCB193, OarFCB266 and OarFCB304 loci. Anim. Genet. 24:145. [ Links ]
Buchanan FC, Galloway SM and Crawford AM (1994) Ovine microsatellites at the OarFCB5, OarFCB19, OarFCB20, OarFCB48, OarFCB129 and OarFCB226 loci. Anim. Genet. 25:60. [ Links ]
Crawford AM, Dodds KG, Ede AJ, Pierson CA, Montgomery GW, Garmonsway HG, Beattie AE, Davies K, Maddox JF, Kappes SW, et al. (1995). An autosomal genetic linkage map of the sheep genome. Genetics. 140:703-724. [ Links ]
de Gortari MJ, Freking BA, Kappes SM, Leymaster KA, Crawford AM, Stone RT and Beattie CW (1997) Extensive genomic conservation of cattle microsatellite heterozygosity in sheep. Animal. Genetics. 28:274-290. [ Links ]
de Gortari MJ, Freking BA, Cuthbertson RP, Kappes SM, Keele JW, Stone RT, Leymaster KA, Dodds KG, Crawford AM and Beattie CW (1998) A second-generation linkage map of the sheep genome. Mamm. Genome. 9:204-209. [ Links ]
Ede AJ, Pierson CA and Crawford AM (1995) Ovine microsatellites at the OarCP9, OarCP16, OarCP20, OarCP21, OarCP23 and OarCP26 loci. Anim. Genet. 26:129-130. [ Links ]
Henry HM, Penty JM, Pierson CA and Crawford AM (1993) Ovine microsatellites at the OarHH35, OarHH41, OarHH44, OarHH47and OarHH64 loci. Anim. Genet. 24:222. [ Links ]
Hayes H, Le Chalony C, Goubin G, Mercier D, Payen E, Bignon C and Kohno K (1996) Localization of ZNF164, ZNF146, GGTA1, SOX2, PRLR and EEF2 on homologous cattle, sheep and goat chromosomes by fluorescent in situ hybridization and comparison with the human gene map. Cytogenet. Cell Genet. 72:342-346. [ Links ]
Jenkins ZA, Henry HM, Galloway SM, Dodds KG and Montgomery GW (1997) Comparative linkage mapping of genes on sheep chromosome 3 provides evidence of chromosomal rearrangements in the evolution of Bovidae. Cytogenet. Cell Genet. 78:272-274. [ Links ]
Majala K (1997) Genetics aspects of domestication, common breeds and their origin. In: Piper L and Rubinsky A (eds) The genetics of sheep. CAB International, U.K , pp. 13-50. [ Links ]
Miller SA, Dikes DD and Polesky HH (1988) A simple salting procedure for extracting DNA from human nucleated cells. Nucleic Acids Res. 16:215. [ Links ]
Mulsant P, Schibler L, Lecerf F, Riquet J, Chitour N, Eggen A, Cribiu EP, Lanneluc I and Elsen JM (1998) Regional mapping of the FecB (Booroola) region of sheep chromosome 6. Anim. Genet. 29:36-37. [ Links ]
Pierson CA, Hanrahan V, Ede AJ and Crawford AM (1993) Ovine microsatellites at the OarVH34, OarVH58, OarVH61 and OarVH72 loci. Anim. Genet. 24:224. [ Links ]
Pritchard JK, Stephens M and Donnelly P (2000) Inference of population structure using multilocus genotype data. Genetics. 155:945-959. [ Links ]
Queller DC, Strassmann JE and Hughes CR (1993) Microsatellites and Kinship. Tree. 8:285-288. [ Links ]
Raymond M and Rousset F (1995) GenePop (version 1.2). Population genetics software for exact test and ecumenisism. J. Hered. 86:248-249. [ Links ]
Sanguinetti CJ, Neto ED and Simpson AJG (1994) Rapid silver staining and recovery of PCR products separated on polyacrilamide gels. Biotechniques. 17:915-918. [ Links ]
Usha AP, Simpson SO and Williams JL (1995) Probability of random sire exclusion using microsatellite markers for parentage verification. Anim. Genet. 26:155-161. [ Links ]
Weir BS (1996) Genetic data analysis II. Sinauer Associates, Inc. Publishers. Sunderland, Massachusetts, 445 pp. [ Links ]