- Research article
Effectiveness of 10 polymorphic microsatellite markers for parentage and pedigree analysis in plateau pika (Ochotona curzoniae)
BMC Geneticsvolume 11, Article number: 101 (2010)
The plateau pika (Ochotona curzoniae) is an underground-dwelling mammal, native to the Tibetan plateau of China. A set of 10 polymorphic microsatellite loci has been developed earlier. Its reliability for parentage assignment has been tested in a plateau pika population. Two family groups with a known pedigree were used to validate the power of this set of markers.
The error in parentage assignment using a combination of these 10 loci was very low as indicated by their power of discrimination (0.803 - 0.932), power of exclusion (0.351 - 0.887), and an effectiveness of the combined probability of exclusion in parentage assignment of 99.999%.
All the offspring of a family could be assigned to their biological mother; and their father or relatives could also be identified. This set of markers therefore provides a powerful and efficient tool for parentage assignment and other population analyses in the plateau pika.
Plateau pikas (Ochotona curzoniae) are small lagomorphs that inhabit the high alpine grasslands of the Tibetan plateau of China. They live in cohesive families and occupy burrow systems. Plateau pikas exhibit monogamy, polygyny, polyandry and promiscuous mating systems . Approximately 57.8% of pikas exhibit philopatry, and dispersal movements are extremely restricted, although some dispersal may occur to ensure spatial separation of kin that may otherwise mate . Inbreeding would be expected to occur under these circumstances. Dominant males monopolize mating in order to maximise reproductive fitness and minimise inbreeding depression. Previous methods to determine the level of inbreeding and how it affects the population depended mainly on direct observation due to the lack of molecular tools. Although family group behaviors have been described through observation in the plateau pika [2–5], details of family structures lack corroborative molecular evidence. In some breeding systems such as lekking, polygyny, polyandry and cooperative breeding, it may be impossible to determine parentage from direct observations . Therefore, molecular tools such as microsatellites markers are necessary to obtain genetic information about family structure, social behavior and dispersal. Microsatellite markers, also called short tandem repeats (STRs), are ideal molecular markers for various genetic studies because they are highly polymorphic, codominant in the manner of inheritance, and easy to genotype by using the polymerase chain reaction (PCR) . Previous studies have demonstrated the ability of microsatellite markers to determine parentage [8, 9], analyze pedigree  and identify species  using species-specific STRs. The success rates of parentage assignment in some studies are high, for example, assignment success for a hatchery population of brown sole  and Chinese shrimp  reached 92.2% and 92.9%, respectively.
In this study, we evaluated the reliability of our previously developed 10 microsatellite loci  for parentage determination and discuss their use for future population genetic analyses of the plateau pika populations.
The mean proportion of individuals genotyped was 0.9556 and expected heterozygosity (He) ranges from 0.784 to 0.905, while PIC ranged from 0.623 to 0.857. The number of alleles (k) per locus ranged from 5 to 11. The PE for individual loci was from 0.351 to 0.887, and the CPE for the 10 loci used in this study was as high as 99.999% (Table1). Moreover, we find that the CPE for a core set of 5 loci in parentage determination reached 99.99% (Table2) and the CPE increased much slowly with each additional locus (Figure1). Combined non-exclusion probabilities for the first parent, second parent and parent pairs are 2.28E-03, 4.9E-05 and 5E-08, respectively.
For the family-8, LOD values for maternity were from 2.21 to 3.86, which indicates that all offspring are correctly assigned to their mother, pF8, and these results are consistent with pedigree information. LOD values for the candidate paternity of family-8 ranges from 2.26 to 5.88 and their assignment to their father, pM8, is consistent with observations in the field (Table2). Nevertheless, there is a mismatched locus between p8-2 and pM8, namely EU518189. Similarity of two individuals based on allele sharing of the 10 loci and the proportion of shared alleles between individuals was calculated (Table3).
For the family-10, LOD values for maternity were 4.29 and 6.27 and the offspring could be correctly assigned to their biological mother, pF10 (Table2). The proportion of shared alleles between pF10 and offspring was 0.60 and 0.65 (Table3). None of the offspring from either family could be assigned to pF3. However, LOD values for paternity of family-10 are below 0. Proportion of shared alleles between p10-1, p10-2 and pM12 are 0.25 and 0.31, respectively. Allele sharing proportion between offspring from family-10 with pM10 is nearly equal to that with pM12 (Table3), with an average of 0.28 versus 0.25, but both are below 0.5. This indicates that the genetic distance between family-10 and pM10 is smaller than that between family-10 and pM12. Thus, the paternity for the embryos in family-10 could not be conclusively attributed to either pM12 or pM10.
The use of microsatellite genotyping for parentage assignment and population genetics studies is a common procedure in wild animals. He and PIC values presented here for a set of microsatellites are high compared to similar studies in other species . Paternity exclusion is a commonly used method for parentage determination. The CPE for our set of 10 microsatellite loci is over 99.99% (Table1). The CPE value for this set of microsatellite markers is higher than that of 15 microsatellites used in the horse for pedigree verification . Figure 1 shows that the combined probability of exclusion is already quite high for a small number of markers. Although it is not yet established whether all the 10 loci used here are autosomal, their application for the purpose of parentage assignment appears to be satisfactory, especially when data from all 10 is analysed together, as was done here.
As demonstrated in this study, all the offspring could be assigned to their true mother and the paternity for family-8 was satisfactorily established, except for a mismatch between embryo, p8-2 and the male pM8. This is probably due to a mutation or a null allele in locus EU518189. For the family-10, a father could not be identified from the candidate fathers. The members from family-10 were genotyped three times to insure veracity of the methods used. There are still 4 loci with genotype mismatches in the paternity-offspring pairs. A negative LOD score means that the candidate parent is not likely to be the true parent. This could be because the true father is probably not present in the candidate pool. The exclusion analysis shows that pM12 is closely related to family-10 with half of the total loci matching, which makes him more likely to be their father than pM10, and this is inconsistent with the observation that pM10 was captured while moving within the family-10 burrow. Furthermore, the allele sharing matrix between embryos from family-10 and pM10 shows that pM10 and pM12 are both closely related with family-10. Taken together, these results demonstrate that pM10 and pM12 may be close relatives of p10-1 and p10-2. This indicates that the set of markers show great potential in not only parentage assignment, but also in determining close relatives. Although pM10 was captured along with family-10, it does not appear to be the biological father of the offspring analyzed here. We conclude that there are other males living in a family, and that philopatry exists in social plateau pikas which is consistent with former observation. In this study, pM12 showed significant genetic similarity to family-10, but wasn't captured within the family-10 home range. This supports Dobson's observation that there exists dispersal in this species .
In conclusion, we have developed a set of 10 microsatellite markers that has proved to be a powerful tool for parentage verification and individual identification of plateau pikas from the Qinghai-Xizang (Tibetan) plateau of China.
Sample collections and molecular techniques
Plateau pikas were trapped by the string-noose method  in the southeastern region of the Tibetan Plateau, about 17 km east of Dawu town (34°24′N, 100°21′E, at an elevation of 3846 m), Maqin County, Qinghai Province, People's Republic of China. This study was approved by China Zoological Society. The captured adult pikas were euthanized humanely by experienced laboratory assistant by cervical dislocation as per the recommendations by the China Wildlife Conservation Association. The animals were dissected and muscle tissue samples were obtained. Pregnant females were anesthetized before dissection. The muscle and embryo tissues were collected and stored in 95% ethanol. Animals from two families, designated family-8 and family-10 were analyzed. Tissue samples obtained from three embryos (p8-1, p8-2, p8-3) from family-8; and two embryos (p10-1, p10-2) from family-10, were from the females, pF8 and pF10, respectively. Based on direct observations, the corresponding probable fathers of these embryos were the males, pM8, pM10. Eight other males from the area were randomly collected and used in the analysis (Figure2. Another female, pF3 was also included as an additional candidate mother.
All the sample tissues were washed with distilled water to remove the ethanol from the surface before use in DNA extraction. Genomic DNA was isolated from the muscle tissues of the parents and their progeny using the standard phenol-chloroform extraction protocol. Negative controls were included to confirm that no contamination had occurred during the extraction procedure. The 10 pairs of microsatellite primers employed here [EU518185, EU518186, EU518196, EU518194, EU518192, EU518193, EU518189, EU518184, EU518187, EU518191] have been described in our earlier work [13, 17]. These 10 loci selected for this study were chosen based on their ease of amplification, particularly when analyzing impure or trace samples that occur frequently in wildlife studies, while the other 3 loci [EU518195, EU518188, EU51819] described in our previous study require a high quality of genomic DNA that might be difficult to obtain in some studies. PCR was carried out in a total volume of 10 μL containing approximately 10 ng of genomic DNA, 1 pmol unlabelled reverse primer, 1 pmol fluorescently labeled forward primer and Taq DNA polymerase (TakaRa), in an appropriate buffer. The other reactions conditions were as follows: an initial denaturation at 94°C for 5 min, followed by 10 cycles of 94°C for 30 s and 72°C for 30 s, and then 20 cycles of 89°C for 30 s and 72°C for 30 s, with a final extension at 72°C for 20 min [12, 14]. Appropriate controls were included for each set of PCR reactions to ascertain that the template and the primers were not contaminated. PCR products were purified using ethanol-ammonium acetate. Genotyping was based on fragment length polymorphism of the fluorescently tagged DNA fragments using the 3730 DNA Analyzer with GeneScan™-500 LIZ™ Size Standard and GeneMarkerV1.75 software (all from Applied Biosystems, Inc.).
Genetic parentage analysis
The allelic size data set was checked for numeric errors and null alleles at 95% confidence interval using MICRO-CHECKER . Hardy-Weinberg equilibrium was tested in GENEPOP version 3.4 . Cervus 3.0  was used for the analysis of all the following parameters: parentage analysis by calculating logarithm of odds (LOD) scores given by the natural logarithm of the overall likelihood ratio, homozygotes (Ho), heterozygotes (He), power of discrimination (PD), polymorphism information content (PIC), average non-exclusion probability for one candidate parent (NE-1P), average non-exclusion probability for one candidate parent given the genotype of a known parent of the opposite sex (NE-2P); average non-exclusion probability for a candidate parent pair (NE-PP), estimated null allele frequency F (Null), average non-exclusion probability for identity of two unrelated individuals (NE-I), average non-exclusion probability for identity of two siblings (NE-SI) . Power of exclusion and combined power of exclusion of the 10 loci were calculated from allele frequencies by PowerStats V1.2 (Promega Corporation, U.S.A.). Furthermore, distances between pairs of individuals based on allele sharing were calculated using Microsatellite Tookit http://www.animalgenomics.ucd.ie/sdepark/ms-toolkit/.
Dobson FS, Smith AT, Gao Xue W: The mating system and gene dynamics of plateau pikas. Behavioural processes. 2000, 51 (1-3): 101-110. 10.1016/S0376-6357(00)00122-4.
Dobson FS, Smith AT, Gao WX: Social and ecological influences on dispersal and philopatry in the plateau pika (Ochotona curzoniae). Behavioral Ecology. 1998, 9 (6): 622-10.1093/beheco/9.6.622.
Qu Jiapeng YM, Wenjing LI, Kexin LI, Zhang Yanming, Smith Andrew: Seasonal variation of family group structure of plateau pikas(Ochotona curzonlae). Acta Theriologica Sinica. 2008, 28 (2): 144-150.
Smith AT, Gao WX: Social relationships of adult black-lipped pikas (Ochotona curzoniae). Journal of Mammalogy. 1991, 231-247. 10.2307/1382094.
Jiapeng Q, Kexin L, Min Y, Wenjing L, Yanming Z, Smith A: Seasonal dynamic pattern of spacial territory in social groups of plateau pikas (Ochotona curzoniae). ACTA THERIOLOGICA SINICA. 2007, 27 (3): 215-220.
Richardson DS, Jury FL, Blaakmeer K, Komdeur J, Burke T: Parentage assignment and extra-group paternity in a cooperative breeder: the Seychelles warbler (Acrocephalus sechellensis). Molecular Ecology. 2001, 10: 2263-2273. 10.1046/j.0962-1083.2001.01355.x.
Estoup A, Gharbi K, SanCristobal M, Chevalet C, Haffray P, Guyomard R: Parentage assignment using microsatellites in turbot (Scophthalmus maximus) and rainbow trout (Oncorhynchus mykiss) hatchery populations. Canadian Journal of Fisheries and Aquatic Sciences. 1998, 55 (3): 715-723. 10.1139/cjfas-55-3-715.
O'Reilly PT, Herbinger C, Wright JM: Analysis of parentage determination in Atlantic salmon (Salmo salar) using microsatellites. Animal Genetics. 1998, 29 (5): 363-370. 10.1046/j.1365-2052.1998.295359.x.
Dong S, Kong J, Zhang T, Meng X, Wang R: Parentage determination of Chinese shrimp(Fenneropenaeus chinensis) based on microsatellite DNA markers. Aquaculture(Amsterdam). 2006, 258 (1-4): 283-288.
Sefc KM, Steinkellner H, Glssl J, Kampfer S, Regner F: Reconstruction of a grapevine pedigree by microsatellite analysis. TAG Theoretical and Applied Genetics. 1998, 97 (1): 227-231. 10.1007/s001220050889.
Li CD, Rossnagel BG, Scoles GJ: The development of oat microsatellite markers and their use in identifying relationships among Avena species and oat cultivars. TAG Theoretical and Applied Genetics. 2000, 101 (8): 1259-1268. 10.1007/s001220051605.
Kim SG, Morishima K, Satoh N, Fujioka T, Saito S, Arai K: Parentage assignment in hatchery population of brown sole Pleuronectes herzensteini by microsatellite DNA markers. Fisheries Science. 2007, 73 (5): 1087-1093. 10.1111/j.1444-2906.2007.01440.x.
Li K, Geng J, Yang J, Zhang Y, Hu S: Isolation and characterization of 13 microsatellite loci in the plateau pika (Ochotonacurzoniae). Conservation Genetics. 2009, 10 (3): 785-787. 10.1007/s10592-008-9662-6.
Mukesh M, Sodhi M, Bhatia S, Mishra B: Genetic diversity of Indian native cattle breeds as analysed with 20 microsatellite loci. Journal of Animal Breeding and Genetics. 2004, 121 (6): 416-424. 10.1111/j.1439-0388.2004.00468.x.
Tozaki T, Kakoi H, Mashima S, Hirota K, Hasegawa T, Ishida N, Miura N, Choi-Miura NH, Tomita M: Population study and validation of paternity testing for Thoroughbred horses by 15 microsatellite loci. The Journal of veterinary medical science/the Japanese Society of Veterinary Science. 2001, 63 (11): 1191-
Dobson F, Smith A, Gao W: Social and ecological influences on dispersal and philopatry in the plateau pika (Ochotona curzoniae). Behavioral Ecology. 1998, 9 (6): 622-10.1093/beheco/9.6.622.
Geng J, Li K, Zhang Y, Hu S: A Modified Enrichment Method to Construct Microsatellite Library from Plateau Pika Genome (Ochotona curzoniae). Genomics, Proteomics & Bioinformatics. 2010, 8 (1): 72-76.
Van Oosterhout C, Hutchinson WF, Wills DPM, Shipley P: MICRO-CHECKER: software for identifying and correcting genotyping errors in microsatellite data. 2004, 4: 535-538.
Raymond M, Rousset F: GENEPOP (version 1.2): population genetics software for exact tests and ecumenicism. Am Genetic Assoc. 1995, 86: 248-
Marshall TC, Slate J, Kruuk LEB, Pemberton JM: Statistical confidence for likelihood-based paternity inference in natural populations. Molecular Ecology. 1998, 7 (5): 639-655. 10.1046/j.1365-294x.1998.00374.x.
We would like to thank Dr. Ji Weihong of Massey University for her comments and careful editing of this manuscript. This work was supported by the Key Innovation Research Program, Chinese Academy of Sciences (Grant KSCX2-YW-N-06; KSCX2-YW-N-90), and '973' Project of China Ministry of Science and Technology (Grant 2007CB109100).
KXL performed all the molecular experiments, analyzed all the data set, drafted manuscript, make the Figure 1, Figure 2, Table 1, made revisions for this manuscript. JNG made table 2 and table 3 and described them. JPQ collected the samples. YMZ took response for all this paper. SNH directed all the experiment and revised this manuscript. All authors read and approved the final manuscript.
Kexin Li, Jianing Geng contributed equally to this work.