Journal of Heredity Advance Access originally published online on April 2, 2007
Journal of Heredity 2007 98(3):232-237; doi:10.1093/jhered/esm009
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Inbreeding and Genetic Structure in the Endangered Sorraia Horse Breed: Implications for its Conservation and Management
From the Departamento de Biologia Animal, Faculdade de Ciências, Centro de Biologia Ambiental, Universidade de Lisboa, Edifício C2-Piso 3, Campo Grande, 1749-016 Lisboa, Portugal (Luís and Oom); and the Department of Veterinary Integrative Biosciences, College of Veterinary Medicine and Biomedical Sciences, Texas A&M University, College Station, TX 77843 (Cothran)
Address correspondence to C. Luís at the address above, or e-mail: cmluis{at}fc.ul.pt.
The Sorraia horse is a closed breed with reduced effective population size and considered in critical maintained risk status. The breed exists in 2 main breeding populations, one in Portugal and one in Germany, with a smaller population size. A set of 22 microsatellite loci was used to examine genetic diversity and structure of the Sorraia horse breed and to compare individual inbreeding coefficient F, estimated from pedigree data, with individual heterozygosity and mean d2. The Sorraia horse shows lower levels of microsatellite diversity when compared with other horse breeds. Due to management strategies, there are clear differences in the genetic structure of the 2 main Sorraia horse populations. Individual heterozygosity was shown to be a good estimator, used together with or as an alternative to inbreeding coefficient, in predicting fitness and evaluating the inbreeding level of the Sorraia horse. The information gathered in this study, combined with information available from previous studies, offers an important and wide information base for the future development of an effective breeding management of the Sorraia horse in order to preserve this endangered breed.
The Sorraia horse is a breed with unusual characteristics. It was recovered and maintained with the purpose of preserving an equine type considered to be representative of the ancestral Iberian horse, from which several domestic horse breeds may have been derived and that probably existed in the Iberian Peninsula since the Palaeolithic period (d'Andrade 1926, 1937, 1945, 1954). In 1920, a small herd of horses with evident primitive characteristics (grullo and dun coat colors, dark muzzle area, black dorsal stripe, black-tipped ears, and striped legs) was accidentally discovered by Ruy d'Andrade in the Sorraia valley. By 1937, in order to preserve a horse that he considered to be of an ancestral type, Ruy d'Andrade selected 10 individuals for the founder population (3 stallions and 7 mares1 of them pregnant) and started a breeding program trying to fix their primitive characteristics in the progeny. The unknown stallion that sired the pregnant mare was also considered a founder, and in 1948, a very striped dun Argentine Criollo stallion was added to the founding population. Since then, no further introductions were made to the population.
In 1975, 2 other farms started breeding the Sorraia horse, and in 1976, a subpopulation was created in Germany with 3 stallions and 3 mares imported from Portugal with no further immigration until the recent past (Oom 1992; Oom and Luís 2001). Although there are now other breeders in Germany, these animals are the founders of the majority of the population bred outside Portugal.
The Sorraia horse is a closed breed, with complete pedigree for all individuals tracing back to the founders. In 2004, the first volume of the Sorraia Studbook was published (Oom et al. 2004) containing 564 animals from the foundation of the breed until 31 December 2002. Less than 200 individuals are alive, and only 10 of the 12 founders are still represented in the living population (Oom et al. 2004; Kjöllerström 2005).
The Sorraia breed has a reduced effective population size and, according to Food and Agriculture Organization (FAO) criteria, is considered in critical maintained risk status (FAO 1998; DAD-IS Database http://dad.fao.org). The low population size has led to high levels of mean kinship (mk) coefficient (0.282
mk
0.423) and average inbreeding (F = 0.363), as estimated from genealogical data analysis (Kjöllerström 2005; Oom and Serrano 2006).
Compelling reasons for conserving genetic resources of domestic livestock, including horses, can be classified into a few categories that include genetic insurance for the species and for agriculture, historical and cultural factors, and concerns for the extinction of variants that might have scientific value (Sponenberg 2000). The main objective of conservation genetics is to preserve variability within populations under the hypothesis of correlation between genetic variation and population viability (Hall and Bradley 1995). Studying homozygosity as an index for inbreeding is a valuable management tool and may be critical for the Sorraia horse.
Assessment of genetic variability in populations has been greatly facilitated by the advent of high-resolution molecular genetic techniques for population genetic studies, thus complementing the evaluation of genetic diversity by pedigree data analysis. Because of their high variability, microsatellites are often used to determine the amount and pattern of genetic variation in populations (Goldstein and Schlotterer 1999). Using microsatellite loci information, some studies have estimated individual heterozygosity and mean d2 to measure the extent of inbreeding in populations (e.g., Coltman et al. 1998; Coulson et al. 1998; Ellegreen 1999; Hedrick et al. 2001; Curik et al. 2003; Markert et al. 2004).
With only 12 founders and no further genetic input into the population since its foundation, the Sorraia horse shows limited genetic variability across several genetic markers, most likely due to founder effect and genetic drift as detected in several previous works (Oom and Cothran 1994; Luís, Bastos-Silveira, et al. 2002; Luís, Cothran, et al. 2002; Luís et al. 2005). There is also a significant indication that inbreeding levels in the Sorraia horse are negatively correlated with both male and female fertility rates and with juvenile survival (Kjöllerström J, do Amaral JP, do Mar Oom M, in preparation). Those problems may be of particular concern as such a small effective population size makes the population even more sensitive to demographic stochasticity.
The aim of this study is to use a set of microsatellite data to examine genetic diversity and structure of the Sorraia horse breed and to compare individual inbreeding coefficient F, estimated from pedigree data, with individual heterozygosity and mean d2. With this study, we intend to gather genetic information that may complement the pedigree information and be useful for managing the genetic diversity of the endangered Sorraia breed.
| Materials and Methods |
|---|
|
|
|---|
Sampling and DNA Extraction
A total of 131 Sorraia horses, belonging to the 2 main populations of this breed, were used as source material for DNA extraction: 99 individuals from Portugal (including individuals from 3 main breeding farms) and 32 individuals from Germany (including the 5 founders of this population). DNA was extracted from fresh whole-blood samples using the high-salt extraction procedure (Montgomery and Sise 1990) and from frozen samples with the QIAamp DNA Blood Mini Kit (Qiagen Inc., Valencia, CA).
Microsatellite Analysis
Twenty-two dinucleotide horse microsatellites, distributed over 16 different autosomes (Table 1), were analyzed: AHT4, AHT5 (Binns et al. 1995), ASB2 (Breen et al. 1997), HMS1, HMS2, HMS3, HMS6, HMS7 (Guérin et al. 1994), HTG4, HTG6, HTG8 (Ellegren et al. 1992), HTG10, HTG14 (Marklund et al. 1994), LEX20, LEX23 (Coogle et al.1996), LEX36, LEX41 (Coogle et al. 1997), NVHEQ18 (Roed et al. 1997), UCDEQ5 (Eggleston-Stott et al. 1996), UCDEQ405, UCDEQ425 (Eggleston-Stott et al. 1997), and VHL20 (Van Haeringen et al. 1994). Polymerase chain reactions (PCRs) and amplification were performed as described by Luís, Cothran, et al. (2002).
|
PCR products were separated by electrophoresis in a 6% Long Ranger gel solution (Cambrex Bio Science, Rockland Inc., Rockland, ME) using an automated fluorescence Li-Cor 4200S sequencer (Li-Cor, Lincoln, NE). Alleles were scored according to PCR product size with RFLPscan 3.1 software (Scanalytics CPS Inc., Rockville, MD).
Data Analysis
The number of individuals typed, mean number of alleles (MNA), the observed heterozygosity (Ho), and the unbiased expected heterozygosity (He) were estimated with the software package GENETIX v. 4.04 (Belkhir et al. 2002). Polymorphic information content (PIC), and the probability of exclusion (PE) were assessed using the computer program CERVUS (Marshall et al. 1998).
The software GENEPOP v. 3.3 (Raymond and Rousset 1995) was used to test for deviations from HardyWeinberg equilibrium (HWE) across all loci and subpopulations applying the exact test with default settings of the Markov chain Monte Carlo methodology (Guo and Thompson 1992). P values were adjusted by applying the Hochberg's step-up Bonferroni method (Hochberg 1988).
F statistics of Wright (1951) were used to analyze genetic structuring of the Sorraia breed. The fixation coefficients (FIS and FST) were estimated according to Weir and Cockerham (1984) using the FSTAT 2.9.3.2 [EC] computer program (Goudet 2002).
The existence of genetic bottlenecks was tested using the program BOTTLENECK (Cornuet and Luikart 1996). In recently bottlenecked populations, the majority of loci will exhibit an excess of heterozygotes, thus exceeding the heterozygosity expected in a population at mutation drift equilibrium. The excess of such heterozygosity was estimated under the assumption of a two-phase microsatellite mutation (TPM) model (Luikart et al. 1998) and under application of the Wilcoxon signed-rank test. BOTTLENECK was further employed to test for an L-shaped distribution of allele frequencies of each marker locus because such a distribution is not typical for populations that have undergone short but severe bottleneck (Luikart et al. 1998).
Individual heterozygosity was calculated as the number of loci at which an individual was heterozygous, divided by the total number of loci at which an individual was scored. Mean d2 was calculated as the squared distance in repeat units between the 2 alleles an individual had at a microsatellite locus, averaged over all loci at which an individual was scored, according to equation:
|
|
Because not all individuals were successfully scored at all loci, individuals that had been scored at fewer than 11 loci were omitted from individual heterozygosity and mean d2 analyses.
The inbreeding coefficients values used are included in the Sorraia Studbook (Oom et al. 2004) and were estimated according to Oom (1992) with all genealogies tracing back to the founders.
Descriptive statistics and the Pearson's correlation between individual heterozygosity, mean d2, and inbreeding coefficient variables were calculated using SAS® software package (RELEASE 8.2; SAS Institute Inc., Cary, NC, 2003). Pearson's correlation was used because the variables have a normal distribution.
| Results |
|---|
|
|
|---|
Genetic Diversity
The total number of alleles found for the 22 microsatellites analyzed was 73, and the MNA per locus was 3.3, with a range of 25 alleles. Ho ranged between 0.088 (HMS2) and 0.705 (LEX36). The highest values for He, PIC, and PE were found for microsatellite HTG4, whereas the lowest ones were found for microsatellite HMS2 (Table 1).
After Hochberg's step-up Bonferroni correction for multiple comparisons (Hochberg 1988), a significant deviation from HWE was only observed for microsatellite HTG10. Allele frequencies are available from the authors on request. Diversity statistics across the main 2 Sorraia populations (Portugal and Germany) are presented in Table 1. Microsatellite HTG10 is monomorphic for the German population. The MNA was higher for the Portuguese population, but the Ho and He were higher for the German population. A positive nonsignificant FIS estimation (P > 0.05) was obtained for the Portuguese population, whereas a highly significant negative estimation (P < 0.001) was obtained for Germany.
For the Portuguese population, under the assumption of the TPM, the Wilcoxon test gave a significant result (P = 0.0165) indicative of a bottleneck. However, the presence of an L-shaped distribution of the allele frequencies does not indicate a recent bottleneck in this population. For the German population, the Wilcoxon test revealed the existence of a highly significant bottleneck (P < 0.001) that was confirmed by the absence of an L-shaped distribution of the allele frequencies.
Population Structure
The pairwise FST coefficient, representing the differentiation of the 2 main Sorraia populations, was highly significant with a value of 0.09 (P < 0.001). Thus, 9% of the microsatellite variability is explained by the subdivision of populations, whereas the remaining variability is explained by the variation within populations.
Inbreeding, Individual Heterozygosity, and Mean d2
Frequency distributions of the mean d2, individual heterozygosity, and inbreeding coefficient are shown in Figure 1. Mean individual heterozygosity and mean d2 were 46.4 ± 11.3% and 36.9 ± 16.7% (in bp2), respectively, whereas mean inbreeding coefficient for the analyzed samples was 32.5 ± 7.2%. There was a highly positive correlation (r = 0.395, P < 0.0001) between individual heterozygosity and mean d2. A nonsignificant negative correlation was observed between inbreeding coefficient and mean d2 (r = 0.116, P = 0.198), whereas a very significant negative correlation was observed between inbreeding coefficient and individual heterozygosity (r = 0.268, P < 0.005).
|
| Discussion |
|---|
|
|
|---|
The Sorraia horse shows lower levels of microsatellite diversity when compared with the majority of other horse breeds (e.g., He: 0.6770.770, Ho: 0.6940.765, MNA: 5.27.8, Cañon et al. 2000; He: 0.6460.732, Ho: 0.6280.671, MNA: 4.77.5, Cunningham 2001; He: 0.751, Ho: 0.732, MNA: 4.5, Morais et al. 2005; He: 0.4420.770, Ho: 0.4520.785, MNA: 2.04.7, Juras et al. 2003; He: 0.675, Ho: 0.663, MNA: 7.1, Achmann et al. 2004; He: 0.6500.765, Ho: 0.6330.777, MNA: 4.338.08, Solis et al. 2005). The German and Portuguese subpopulations differ from each other in variability. As expected, due to founder effect, the German population has a reduced number of alleles for most loci analyzed compared with the Portuguese one; however, no loss of founders' alleles at tested loci has occurred since the foundation of this population due to genetic drift. Both Ho and He were higher in the German population, despite statistical evidence for the genetic bottleneck that took place in 1976.The lower heterozygosity found in Portugal is most probably caused by the breeding management used in most Portuguese breeding farms with breeding done in pasture by free live-cover, and the use of only 1 stallion per year per herd, with this same stallion often used during several breeding years. This strategy is completely different from the one used in Germany, where breeding is mostly done by hand-cover with different males used almost each year and more then one per herd, a management strategy that appears to be more suited for genetic conservation because in this population the genetic variability parameters are higher.
As expected, the significant FST coefficient estimated indicates a relatively low gene flow between the Portuguese and German subpopulations sampled due to the reproductive isolation that existed between them until very recently when some Sorraia from Portugal were taken to Germany by new German breeders. Thus, the management strategies implemented in the 2 populations seem to be the main reason for the existence of different genetic structures in both populations.
Founder effect, genetic drift, and artificial selection/breeding management are among the factors that likely contributed to the differentiation within the Sorraia horse breed, owing to the low size of the analyzed populations. Migration can be considered null because there has been no migration between populations since fragmentation, at least until very recently and after sampling collection for this study was completed. Also, in most domestic species, the drift process can be accelerated by the unequal representation of both sexes, which is the case of the Sorraia where there are far fewer males in the breeding population.
A major conclusion from this study related to genetic conservation in the Sorraia breed is that the Portuguese breeders should try to focus the breeding management strategy in increasing heterozygosity. Some efforts have been made, over the last 3 breeding years, by the Sorraia Breeder's Association to advise Portuguese breeders of the best suitable stallions to cover the mares from different herds, based on pedigree information and subsequent genetic analysis performed with GENES software v. 11.8 (Lacy 1998). It will be interesting to analyze the results of this practice in the future in terms of molecular genetic variation within the Portuguese population and the Sorraia total population.
Following a strategy that is fully recommended in conservation plans, in recent years, and after the sampling for this study was performed, the Portuguese population has been further fragmented into different and smaller subpopulations. The intention is to decrease the loss of the genetic variability per generation of the overall population by promoting a more effective exchange of stallions between subpopulations (which will undergo different local effects of genetic drift) and to reduce kinship between animals of different subpopulations.
Pedigrees are essential to provide information for breeding management plans, on the assumption that data registered is correct. A panel of 15 microsatellites showing a PE of 99.8% is now used for parentage testing in the Sorraia horse. The set of microsatellites used in this study would provide superior power of exclusion and thus would provide a superior panel for paternity testing of the Sorraia horse.
For the Sorraia horse, a population that is subject to close inbreeding, individual heterozygosity can more effectively be used to evaluate inbreeding level of the population and, so, inbreeding depression. Therefore, within the Sorraia breed, this molecular estimator could be used, together with, or alternatively to inbreeding coefficient, in predicting fitness. The results from this work support the findings by Tsitrone et al. (2001), Slate and Pemberton (2002), and Markert et al. (2004) stating that heterozygosity is a more robust measure of inbreeding than is mean d2 and that multilocus heterozygosity data can substitute or provide a useful complement to data on inbreeding depression obtained through pedigree analysis.
The results of the present study are in accordance with the known history of the Sorraia horse breed that is a highly inbred, small, and closed population. As stated by Fernández et al. (2004), information obtained from the analysis of neutral molecular markers can be used to make decisions about individuals' contributions to the next generation in order to maximize gene or allelic diversity. Therefore, the information gathered in this study, combined with information available from pedigree and reproductive fitness analysis (Oom 1992; Matos 1996; Oom and Luís 2001; Kjöllerström 2005), from morphological measures (Oom 1992), from mitochondrial DNA markers (Luís, Cothran, et al. 2002; Luís, Bastos-Silveira, Cothran, and Oom 2006; Luís, Bastos-Silveira, Cothran, Costa-Ferreira, et al. 2006), and immunogenetics (Luís et al. 2005), offers an important and wide information base for the future development of effective breeding management of the Sorraia horse in order to preserve this endangered breed.
| Acknowledgments |
|---|
We would like to thank Joana Morais, Ana Rita Grosso, and Joana Cavaco Silva for help provided in the laboratory. C.L. was supported by a PhD grant (SFRH/BD/3318/2000) from the Portuguese Foundation for Science and Technology (FCT/MCT).
| Footnotes |
|---|
Corresponding Editor: Ernest Bailey
Received August 7, 2006
Accepted February 12, 2007
| References |
|---|
|
|
|---|
-
Achmann R, Curik I, Dovc P, Kavar T, Bodo I, Habe F, Marti E, Sölkner J, Brem G. Microsatellite diversity, population subdivision and gene flow in the Lipizzan horse. Anim Genet (2004) 35:285292.[CrossRef][ISI][Medline]
Belkhir K, Borsa P, Chikhi L, Raufaste N, Bonhomme F. GENETIX version 4.04 logiciel sous WindowsTM pour la génétique des populations. (2002) Montpellier (France): Laboratoire Génome, Populations, Interactions CNRS UMR 5000, Université Montpellier II.
Binns MM, Holmes NG, Holliman A, Scott AM. The identification of polymorphic microsatellite loci in the horse and their use in thoroughbred parentage testing. Br Vet J (1995) 151:916.[ISI][Medline]
Breen M, Lindgren G, Binns MM, Norman J, Irvin Z, Bell K, Sandberg K, Ellegren H. Genetical and physical assignments of equine microsatellitesfirst integration of anchored markers in horse genome mapping. Mamm Genome (1997) 8:267273.[CrossRef][ISI][Medline]
Cañon J, Checa ML, Carleos C, Vega-Pla JL, Vallejo M, Dunner S. The genetic structure of Spanish Celtic horse breeds inferred from microsatellite data. Anim Genet (2000) 31:3948.[CrossRef][ISI][Medline]
Coltman DW, Bowen WD, Wright JM. Birth weight and neonatal survival of harbour seal pups are positively correlated with genetic variation measure by microsatellites. Proc R Soc Lond B (1998) 265:803809.[Medline]
Coogle L, Reid R, Bailey E. Equine dinucleotide repeat loci LEX015-LEX024. Anim Genet (1996) 27:217218.[ISI][Medline]
Coogle L, Reid R, Bailey E. Equine dinucleotide repeat loci LEX034-LEX048. Anim Genet (1997) 28:309.[ISI][Medline]
Cornuet JM, Luikart G. Description and power analysis of two tests for detecting recent population bottlenecks from allele frequency data. Genetics (1996) 144:20012014.[Abstract]
Coulson TN, Pemberton JM, Albon SD, Beaumont M, Marshall TC, Slate J, Guinness FE, Clutton-Brock TH. Microsatellites reveal heterosis in red deer. Proc R Soc Lond B (1998) 265:489495.[Medline]
Cunningham EP, Dooley JJ, Splan RK, Bradley DG. Microsatellite diversity, pedigree relatedness and the contributions of founder lineages to thoroughbred horses. Anim Genet (2001) 32:360364.[CrossRef][ISI][Medline]
Curik I, Zechner P, Sölkner J, Achman R, Bodo I, Dovc P, Kavar T, Marti E, Brem G. Inbreeding, microsatellite heterozygosity, and morphological traits in Lipizzan horses. J Hered (2003) 94:125132.
d'Andrade R. Apontamentos para um estudo sobre a origem e domesticação do cavalo na Península Ibérica (1926) Lisboa (Portugal): Aproximações Centro Tipográfico Colonial.
d'Andrade R. Les chevaux du Sorraia. C R XII Congre Int Zool (1937) 236:82370.
d'Andrade R. O Cavalo do Sorraia. Bol Pecu (1945) 13:113.
d'Andrade R. Alredor del Caballo Español. (1954) Lisboa (Portugal): Colección de Estúdios.
Eggleston-Stott ML, DelValle A, Bautista M, Dileanis S, Wictum E, Bowling AT. Nine equine dinucleotide repeats at microsatellite loci UCDEQ136, UCDEQ405, UCDEQ412, UCDEQ425, UCDEQ437, UCDEQ467, UCDEQ487, UCDEQ502 and UCDEQ505. Anim Genet (1997) 28:370371.[Medline]
Eggleston-Stott ML, Valle A, Bowling AT, Bautista M, Malyi W. Four equine dinucleotide repeats at microsatellite loci UCDEQ5, UCDEQ14, UCDEQ46 and UCDEQ62. Anim Genet (1996) 27:129.[ISI][Medline]
Ellegren H, Johansson M, Sandberg K, Andersson L. Cloning of highly polymorphic microsatellites in the horse. Anim Genet (1992) 23:133142.[ISI][Medline]
Ellegren H. Inbreeding and relatedness in Scandinavian grey wolves Canis lupus. Hereditas (1999) 130:239244.[CrossRef][ISI][Medline]
FAO. Secondary Guidelines for Development of National Farm Animal Genetic Resources Management Plans (1998) FAO (Rome): Management of Small Populations at Risk. Initiative for Domestic Animal Diversity.
Fernández J, Toro MA, Caballero A. Managing individuals' contributions to maximize the allelic diversity maintained in small, conserved populations. Conserv Biol (2004) 18:13581367.[CrossRef]
Goldstein DB, Schlotterer C. Microsatellites: evolution and application. (1999) Oxford: Oxford University Press.
GoudetJ. 2002. FSTAT, a program to estimate and test gene diversities and fixation indices (version 2.9.3.2) [Internet]. Available from: http://www2.unil.ch/popgen/softwares/fstat.htm.
Guérin G, Bertaud M, Amigues Y. Characterization of seven new horse microsatellites: HMS1, HMS2, HMS3, HMS5, HMS6, HMS7, and HMS8. Anim Genet (1994) 25:62.[ISI][Medline]
Guo SW, Thompson EA. Performing the exact test of Hardy-Weinberg proportions for multiple alleles. Biometrics (1992) 48:361372.[CrossRef][ISI][Medline]
Hall SJG, Bradley DG. Conserving livestock breed biodiversity. TREE (1995) 10:267270.
Hedrick P, Fredrickson R, Ellegren H. Evaluation of d2, a microsatellite measure of inbreeding and outbreeding, in wolves with a known pedigree. Evolution (2001) 55:12561260.[CrossRef][ISI][Medline]
Hochberg Y. A sharper Bonferroni procedure for multiple tests of significance. Biometrika (1988) 75:800802.
Juras R, Cothran EG, Klimas R. Genetic analysis of three Lithuanian horse breeds. Acta Agric Scand (2003) 53:180185.
Kjöllerström HJ. O Studbook como instrumento na análise genética e demográfica de populações ameaçadas: o caso do Cavalo do Sorraia (2005) Faculdade de Ciências, Universidade de Lisboa: Lisboa (Portugal): Relatório de Estágio.
Lacy R. GENES: a computer program for pedigree analysis and genetic management. (1998) Chicago (IL): Chicago Zoological Society.
Luikart G, Sherwin WB, Steele BM, Allendorf FW. Usefulness of molecular markers for detecting population bottlenecks and monitoring genetic change. Mol Ecol (1998) 7:963974.[CrossRef][Medline]
Luís C, Bastos-Silveira C, Cothran EG, Oom MM. Variation in the mitochondrial control region sequence between the two maternal lines of the Sorraia horse breed. Genet Mol Biol (2002) 25:309311.
Luís C, Bastos-Silveira C, Cothran EG, Oom MM. Iberian origins of New World horse breeds. J Hered (2006) 97:107113.
Luís C, Bastos-Silveira C, Cothran EG, Costa-Ferreira J, Oom MM. A lost Sorraia maternal lineage found in the Lusitano horse breed. J Anim Breed Genet (2006) 123:399402.[CrossRef][ISI][Medline]
Luís C, Cothran EG, Oom MM. Microsatellites in Portuguese autochthonous horse breeds: usefulness for parentage testing. Genet Mol Biol (2002) 25:131134.
Luís C, Cothran EG, Oom MM, Bailey E. Major histocompatibility complex locus DRA polymorphism in the endangered Sorraia horse and related breeds. J Anim Breed Genet (2005) 122:6972.[CrossRef][ISI][Medline]
Markert JA, Grant PR, Grant BR, Keller LF, Coombs JL, Petren K. Neutral locus heterozygosity, inbreeding, and survival in Darwin's ground finches (Geospiza fortis and G. scandens). Heredity (2004) 92:306315.[CrossRef][ISI][Medline]
Marklund S, Ellegren H, Eriksson S, Sandberg K, Andersson L. Parentage testing and linkage analysis in the horse using a set of highly polymorphic microsatellites. Anim Genet (1994) 25:1923.[ISI][Medline]
Marshall TC, Slate J, Kruuk LEB, Pemberton JM. Statistical confidence for likelihood-based paternity inference in natural populations. Mol Ecol (1998) 7:639655.[CrossRef][Medline]
Matos S. Contributo para a gestão genética do cavalo do Sorraia (1996) Lisbon (Portugal): Relatório de Estágio, Faculdade de Ciências, Universidade de Lisboa.
Montgomery GW, Sise JA. Extraction of DNA from sheep white blood cells. N Z J Agric Res (1990) 33:437441.
Morais J, Oom MM, Malta-Vacas J, Luís C. Genetic structure of an endangered Portuguese semiferal pony breed, the Garrano. Biochem Genet (2005) 43:347364.[CrossRef][ISI][Medline]
Oom MM. O Cavalo Lusitano. Uma raça em recuperação [PhD dissertation]. (1992) Lisbon (Portugal): University of Lisbon.
Oom MM, Andrade JLd', Costa-Ferreira J. Stud book da raça Sorraia. Portugal: Associação Internacional de Criadores Cavalo Sorraia (2004).
Oom MM, Cothran EG. The genetic variation of an endangered breed: the Sorraia horse. Anim Genet (1994) 25:15.
Oom MM, Luís C. The Sorraia horse (Portugal): genetic variation and genetic management of an endangered breed. In: Conservation Genetics Meeting, September 1214 (2001) Switzerland: Lausanne.
Oom MM, Serrano J. O cavalo do Sorraia. Um património histórico e genético a preservar. (2006) Lisboa (Portugal): Ed. Cosmos.
Raymond M, Rousset F. GENEPOP (version 3.3): population genetics software for exact tests and ecumenicism. J Hered (1995) 86:248249.
Roed KH, Midthjell L, Bjornstad G, Olsaker I. Equine dinucleotide repeat microsatellites at the NVHEQ5, NVHEQ7, NVHEQ11, NVHEQ18 and NVHEQ24 loci. Anim Genet (1997) 28:381382.[ISI][Medline]
Slate J, Pemberton JM. Comparing molecular measures for detecting inbreeding depression. J Evol Biol (2002) 15:2031.[CrossRef][ISI]
Solis A, Jugo BM, Mériaux JC, Iriondo M, Mazón LI, Aguirre AI, Vicario A, Estomba A. Genetic diversity within and among four south European native horse breeds based on microsatellite DNA analysis: implications for conservation J Hered. 9 (2005) 6:670678.
Sponenberg DP. Genetic resources and their conservation. In: The genetics of the horseBowling AT, Ruvinsky A, eds. (2000) Wallingford: CABI Publishing. 387410.
Tsitrone A, Rousset F, David P. Heterosis, marker mutational processes and population inbreeding history. Genetics (2001) 159:18451859.
Van Haeringen H, Bowling AT, Stott ML, Lenstra JA, Zwaagstra KA. A highly polymorphic horse microsatellite locus: VHL20. Anim Genet (1994) 25:207.[Medline]
Weir BS, Cockerham CC. Estimating F-statistics for the analysis of population structure. Evolution (1984) 38:13581370.[CrossRef][ISI]
Wright S. The genetical structure of populations. Ann Eugen (1951) 15:323354.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
