Genetic variability of 10 microsatellite markers in the characterization of Brazilian Nellore cattle (Bos indicus)

ANIMAL GENETICS

RESEARCH ARTICLE

Genetic variability of 10 microsatellite markers in the characterization of Brazilian Nellore cattle (Bos indicus)

Marcelo Cervini^I; Flávio Henrique-Silva^II; Norma Mortari^I; Euclides Matheucci Jr^I

^ILaboratório de Imunogenética, Departamento de Genética e Evolução, Universidade Federal de São Carlos, São Carlos, SP, Brazil

^IILaboratório de Biologia Molecular, Departamento de Genética e Evolução, Universidade Federal de São Carlos, São Carlos, SP, Brazil

^{Send correspondence to}

ABSTRACT

We assessed the polymorphism of 10 microsatellites in Brazilian Nellore cattle (Bos indicus) using a commercial multiplex system. Allele frequencies, polymorphism information content, heterozygosity and exclusion probability were calculated. Allele frequencies revealed that in the sample analyzed the markers were not equally polymorphic. The exclusion probabilities and the polymorphism information content of some loci in Nellore cattle were lower than in Bos taurus breeds. When all the microsatellites were considered the combined exclusion probability was 0.9989. This multiplex analysis can contribute toward pedigree information, adequate genetic improvements and breeding programs.

Key words: alleles, frequencies, microsatellite, Nellore, polymorphism, zebu.

INTRODUCTION

The Bos indicus Nellore herd is one of the largest commercial beef herds in the world and is well-adapted to tropical regions. According to the Brazilian Ministry of Agriculture, Livestock and Supply (http://www.agricultura. gov.br), the Nellore herd is the most important beef herd in Brazil, where the total number of both purebred and crossbreed Nellore cattle totals over 140 million head. Due to genetic improvement programs and adequate international sanitary standards, Brazil has ranked as the top beef exporter since 2003, with an export volume exceeding 1.18 million tons (Brazilian Association of Meat Export Industries http://www.abiec.com.br). Accurate pedigree information is essential to maintaining the quality of breed improvement programs and molecular markers have become an important genetic tool in animal genetics studies, allowing the analysis of genetic variability within and between herds. Many of the current molecular marker techniques are based on variations of the polymerase chain reaction (PCR) such as random amplification of polymorphic DNA (RAPD).

Microsatellites markers have been widely used as a genetic markers in bovine population studies and pedigree verification (Visscher et al. 2002, Hansen et al, 2002, Ibeagha-Awemu and Erhardt, 2005), mainly because of their large polymorphism information content, widespread distribution in the eukaryotic genome (Tautz and Renz, 1984) and robust methodology. Microsatellites have been effective in evaluating differences within cattle breeds and in determining population substructures (MacHugh et al, 1998; Ciampolini et al., 1995). More than 1400 microsatellites have been mapped in the cattle genome (Luikart et al., 1999) and some of them have been employed in population genetics studies and parentage verification. Many microsatellite loci have been used in Nellore improvement programs but, to date, there have been no reports of pedigree verification studies using microsatellite markers, pedigree verification in Brazilian livestock currently being based on blood groups and biochemical polymorphism analyses.

The aim of the study described in this paper was to characterize Brazilian Nellore cattle through the analysis of the genetic variability of ten microsatellite markers and to evaluate if these markers are informative in parentage tests.

Materials and Methods

Sample collection and DNA extraction

We sampled 200 unrelated adult Nellore cattle (150 dams and 50 bulls) registered in their breeding associations and randomly selected from private and research herds belonging to 43 farms located in various regions of Brazil. Blood samples were collected in heparinized glass tubes and total genomic DNA isolated as described by Debomoy et al. (1991) and stored at -20 °C.

Microsatellite amplification

As recommended by the International Society of Animal Genetics (ISAG), ten microsatellites (Table 1) were selected for the analysis, using the Stockmarks for Cattle Bovine Genotyping Kit (Applied Biosystems Division, Perkin-Elmer, Foster City, CA). Multiplex amplification was carried out in a final volume of 15 µL containing 50 ng of template DNA, 0.5 units of AmpliTaq Gold^TM polymerase (PE Applied Biosystems, Foster City, CA), 3.0 µL Stockmarks Buffer, 400 µM of each dNTP and 5.5 µL of primer mix (Table 1). The reactions were carried out using a Programmable Thermal Controller PTC-100^TM (MJ Research, INC) in an initial denaturation phase of 15 min at 95 °C, followed by 31 cycles of 45 s at 94 °C, 45 s at 61 °C and 1 min at 72 °C. A final extension was carried out at 72 °C for 1 h and then at 25 °C for 2 h. After amplification, 90 µL of water was added to the tubes and 0.4 µL of this solution was mixed with 2 µL loading mix (DI formamide: dye: GS350Rox - 6:1:1) and analyzed in a 6% (w/v) denaturing gel using an ABI PRISM^TM 377 DNA Sequencer. The fluorescence data was collected by GeneScan^TM Analysis 2.0 and analyzed using Genotyper^TM 2.0 software.

Data analysis

The GENEPOP package Version 3.4 (Raymond & Rousset, 1995) was used to calculate an exact test for deviation from Hardy-Weinberg equilibrium (HWE), allele frequencies and heterozygotic deficiency. Since the microsatellite loci have more than four alleles, an unbiased estimate of the exact HWE probability was calculated using the Markov chain method of Guo & Thompson (1992). The gene diversity (D) was calculated with FSTAT 2.9.3.2 (Goudet, 2001). Exclusion probability (EP), combined exclusion probability (CEP), expected heterozygosity (He) and observed heterozygosity (Ho), and polymorphism information content (PIC) were calculated using Cervus 2.0 software (Marshall et al., 1998).

RESULTS

Ninety-four alleles were detected from the 10 loci surveyed, yielding a mean value of 9.4 alleles per locus. The allele frequencies of 10 microsatellites are listed in Table 2. Allele frequencies revealed that not all markers were equally informative. The TGLA227, BM1824 and TGLA53 loci each had one allele with a much higher frequency than the other alleles (75 bp, 180 bp and 160 bp respectively). The loci ETH10 and ETH3 each had two alleles with high frequencies (209 bp -207 bp and 115 bp -117 bp, respectively). The number of alleles per locus ranged from six for TGLA227 to 16 for TGLA122. The TGLA122 locus showed the highest allele polymorphism, while the INRA023 locus displayed the highest exclusion probability. Six loci (TGLA 53, ETH10, ETH3, ETH225, TGLA122 and INRA023) deviated significantly (p < 0.05) from HWE. A significant deficit of heterozygosity (p < 0.01) was detected in the TGLA53, ETH10 ETH225, TGLA122, and INRA023 loci. The ETH3 locus did not show heterozygote deficiency, although its P value was close to p < 0.01. The mean PIC value was 0.640 and the mean expected heterozygosity value was 0.679. Expected and observed heterozygosity, probability of exclusion and PIC values are shown in Table 3. The combined probability of parentage exclusion for the 10 microsatellites was 0.9989.

DISCUSSION

Accurate cattle pedigree information is essential for the optimal development of breed and selection programs, improving productivity in the animal industry. Misidentification of parentage can lead to breeding inaccuracy, causing great financial losses in herd management and in the beef industry. Geldermann et al. (1986) estimated misidentification rates of 13% using blood group factors and biochemical polymorphisms in cattle. Ron et al. (1996) found a 5% misidentification rate using microsatellite analysis in Israeli dairy cattle. Rosa (1997) reported a misidentification rate of 15% in Brazilian livestock, based on restriction fragment length polymorphism (RFLP) and microsatellite analysis. Microsatellites are the most widely used molecular markers in pedigree control. The use of microsatellites with high polymorphism information content would help to correctly identify individual cattle, allowing for the better operation of cattle breeding programs.

Little information is available regarding the allele frequencies of the ten microsatellites studied in this work and the other variability estimates for Nellore cattle and, to date, there are no estimates of the multiplex variability in Brazilian cattle breeds. Since the evaluation of polymorphism is strictly dependent on the allele number and the frequency distribution of the alleles, estimates of allele frequencies are essential.

A comparison of the results obtained for B. indicus Nellore cattle with those of B. taurus breeds indicated a difference in variability for some loci, which are highly informative in B. taurus but less informative in B. indicus (Nellore). The exclusion probability values for the TGLA227, ETH10 and TGLA53 loci of taurine cattle described by Peelman et al. (1998) and Heyen et al. (1997) are much higher than that of Nellore cattle.

According to Peelman et al. (1998), who analyzed Belgium cattle, the number of TGLA53 locus alleles in Holstein Friesian (13 alleles), Belgian Red Pied (12 alleles), East Flemish (12 alleles) and Belgian Blue (10 alleles) cattle were very similar to that found in Nellore cattle (13 alleles). However, we found that the exclusion probability for the TGLA53 locus in Brazilian Nellore (EP = 0.256) cattle is much lower than in the four Belgian breeds (Holstein Friesian = 0.742, Belgian Red Pied = 0.711, East Flemish = 0.698 and Belgian Blue = 0.682). We obtained similar results for the TGLA227 locus (EP 0.230), much lower than the values described by Heyen (1997) for Holstein (0.69), Red Angus (0.63) and Gelbvieh (0.68) cattle. Thus, the effectiveness of these markers in European B. taurus cattle is not always the same for Indian B. indicus zebu (Brahman, Nellore) cattle. The substitution of the markers with low variability values for others with improved EP values could render this multiplex more efficient for pedigree verification and individual identification in Nellore. Characterization of Brazilian cattle breeds with microsatellite loci is useful to identify informative markers for each breed, optimizing parentage tests along with the variability values of each marker, thus using the least number of markers with higher levels of information while simultaneously facilitating genotypic identification.

The combined exclusion probability value for the 10 loci was 0.9989, an acceptable value more than ideal for parentage tests (Baron et al, 2002). Jia et al (2004) showed that the CEP value was 0.9957 for Holstein Friesian cattle using six microsatellite markers, while Radko et al (2002) obtained a CEP value of 0.9999 using 11 microsatellites and it is known that the CEP values found in Nellore cattle is lower than that of other taurine breeds (Heyen et al., 1997).

In our study we found significant (p < 0.01) deviations from HWE for six loci (TGLA122, INRA023, TGLA53, ETH10, ETH225 and ETH3). Machado et al (2003) also found significant deviations from HWE for Nellore, Gyr and Guzerat cattle breeds using microsatellites markers. Almeida et al (2000) found that the TGLA122 locus was in HWE in Brazilian hybrid bovine breed (5/8 Aberdeen Angus x 3/8 Nellore). We found deviations from HWE caused by heterozygote deficiency at the TGLA122, INRA023, TGLA53, ETH10 and ETH225 loci. Beja et al. (2003) and Loftus et al. (1999) found deviations from HWE in other European bovine populations, also caused by a heterozygosity deficit, and similar results have been reported by Loftus (1999) in six populations, including Indian Nellore cattle.

Several factors can lead to heterozygote deficiency, including null alleles, assortive mating, the Wahlund effect, selection against heterozygotes, inbreeding, or a combination of these. Null alleles are alleles that are not amplified (usually due to a mutation in one of the primer binding sites) and are commonly reported in microsatellite studies as being the source of heterozygosity deficit (Pemberton et al., 1995). The frequency of microsatellite loci containing null alleles has proved to be as high as 30% in humans (Callen et al., 1993). In paternity tests, an undetected null allele may have profound consequences, since it may cause rejection of an otherwise correctly assigned parent (Holm et al., 2001).

To date, there are no reports of studies on Nellore cattle indicating the presence of null alleles for the markers analyzed, although the presence of null alleles has previously been observed in segregation analyses using other microsatellite loci in Nellore cattle (Tambasco et al., 2000). This hypothesis cannot be excluded because segregation analysis using the loci evaluated in this study has not yet been undertaken for Nellore cattle.

Despite the paucity of information provided by some of the loci analyzed in this study, the use of this multiplex analysis proved efficient in Nellore characterization and can be used in pedigree verification.

Acknowledgements

Marcelo Cervini has a fellowship from the Brazilian agency Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES). This work was financially supported by Fundação de Apoio Institucional Desenvolvimento Científico e Tecnológico - Universidade Federal de São Carlos (FAI-UFSCar).

Rosa AJM (1997) Caracterização da raça Nelore e testes de paternidade por marcadores moleculares. Master's Thesis, Escola Superior de Agricultura "Luiz de Queiroz", Universidade de São Paulo, Piracicaba.

Internet Resources

Goudet J (2001) FSTAT, a program to estimate and test gene diversities and exation indices (version 2.9.3). Available from http://www.unil.ch/izea/softwares/fstat.html.

Received: May 11, 2005; Accepted: November 16, 2005.

Associate Editor: Pedro Franklin Barbosa

• Baron EE, Martinez ML, Verneque RS and Coutinho LL (2002) Parentage testing and effect of misidentification on the estimation of breeding value in Gir cattle. Genet Mol Biol 25:389-394.
• Barendse W, Armitage SM, Kossarek LM, Shalom A, Kirkpatrick BW, Ryan AM, Clayton D, Li L, Neibergs HL, Zhang N, Grosse WM, Weiss J, Creighton P, McCarthy F, Ron M, Teale AJ, Fries R, McGraw RA, Moore SS, Georges M, Soller M, Womack JE and Hetzel DJS (1994) A genetic linkage map of the bovine genome. Nature Genet 6:227.
• Beja-Pereira A, Alexandrino P, Bessa I, Carretero Y, Dunner S, Ferrand N, Jordana J, Laloe D, Moazami-Goudarzi K, Sanches A and Cañon J (2003) Genetic characterization of Southwestern European bovine breeds: A historical and biogeographical reassessment with a set of 16 microsatellites. J Hered 94:243-50
• Bishop MD, Kappes SM and Keele JW et al (1994) A genetic linkage map for cattle. Genetics 136:619-39.
• Botstein D, White RL, Skolnick M and Davis RW (1980) Construction of a genetic linkage map in man using restriction fragment length polymorphisms. Amer J Human Genet 32:314-31.
• Callen DF, Thompson AD, Shen Y, Phillips HA, Richards RI, Mulley JC and Sutherland GR (1993) Incidence and origin of null alleles in the (AC)n microsatellite markers. Amer J Human Genet 52:922-27.
• Ciampolini R, Moazami-Goudarzi K, Vaiman D, Dillman C, Mazzanti E, Foulley JL, Leveziel H and Cianci D (1995) Individual multilocus genotypes using microsatellite polymorphisms to permit the analysis of the genetic variability within and between Italian beef cattle breeds. J Animal Sci 73:3259-68.
• Curi RA and Lopes CR (2002) Evaluation of nine microsatellite loci and misidentification paternity frequency in a population of Gyr breed bovines. Braz J Veterin Res Animal Sci 39:129-35.
• Debomoy KL and Nurnberger JI (1991) A rapid non-enzymatic method for preparation of HMW DNA from blood for RFLP studies. Nucleic Acids Res 19:5444.
• Geldermann H, Pieper U and Weber WE (1986) Effect of misidentification on the estimation of breeding value and heritability in cattle. J Animal Sci 63:1759-68.
• Guo S and Thompson EA (1992) Performing exact test of Hardy-Weinberg proportions for multiple alleles. Biometrics 48:361-72.
• Hansen C, Shrestha JN, Parker RJ, Crow GH, McAlpine PJ and Derr JN (2002) Genetic diversity among Canadienne, Brown Swiss, Holstein, and Jersey cattle of Canada based on 15 bovine microsatellite markers. Genome 45:897-904.
• Heyen DW, Beever JE, Da Y, Evert RE, Green C, Bates SRE, Ziegle JS and Lewin HA (1997) Exclusion probabilities of 22 bovine microsatellite markers in fluorescsnt multiplexes for semi-automated parentage testing. Animal Genet 28:21-27.
• Holm LE, Loeschcke V and Bendixen C (2001) Elucidation of the molecular basis of a null allele in a rainbow trout microsatellite. Marine Biotechnol 3:555-60.
• Ibeagha-Awemu EM and Erhardt G (2005) Genetic structure and differentiation of 12 African Bos indicus and Bos taurus cattle breeds, inferred from protein and microsatellite polymorphisms. J Anim Breed Genet 122:12-20
• Jamieson A and Taylor St CS (1997) Comparisons of three probability formulae for parentage exclusion. Animal Genet 28:397-400.
• Loftus R, Ertrugrul O, Harba A, El-Barodys M, MacHugh D, Park S and Bradley D (1999) A microsatellite survey of cattle from a centre of origin: The Near East. Mol Ecol 8:2015-22.
• Luikart G, Biju-Duval M-P, Ertugrul O, Zagdsuren Y, Maudet C and Taberlet P (1999) Power of 22 microsatellite markers in fluorescent multiplexes for parentage testing in goats (Capra hircus). Animal Genet 30:431-38.
• MacHugh DE, Loftus RT, Bradley DG, Sharp PM and Cunningham P (1994) Microsatellite DNA variation within and among European cattle breeds. Proc Royal Soc 256:25-31.
• MacHugh DE, Loftus RT, Cunningham P and Bradley DG (1998) Genetic structure of seven European cattle breeds assessed using 20 microsatellite markers. Animal Genet 29:333-40.
• Marshall TC, Slate J, Kruuk LEB and Pemberton JM (1998) Statistical confidence for likelihood-based paternity inference in natural population. Mol Ecol 7:639-55.
• Peelman LJ, Mortiaux F, Van Zeveren A, Dansercoer A, Mommens G, Coopman F, Bouquet Y, Burny A, Renaville R and Portelle D (1998) Evaluation of the genetic variability of 23 bovine microsatellite markers in four Belgian cattle breeds. Animal Genet 29:161-67.
• Pemberton JM, Slate J, Bancroft DR and Barret JA (1995) Nonamplifying alleles at microsatellite loci: A caution for parentage and population studies. Mol Ecol 4:249-52.
• Radko A, Duniec M, Zabek T, Janik A and Natonek M (2002) Polymorphism of 11 microsatellite DNA sequence and their usefulness for paternity control in cattle. Polish Soc Veterin Sci 58:708-710.
• Raymond M and Rousset F (1995) GENEPOP (version 1.2): Population genetics software for exact tests and ecumenicism. J Hered 86:248-49.
• Ron M, Blanc Y, Band M, Ezra E and Weller JI (1996) Misidentification rate in the Israeli Dairy Cattle population and its implications for genetic improvement. J Dairy Sci 79:676-81.
• Steffen P, Eggen A, Dietz AB, Womack JE and Stranzinger G (1993) Isolation and mapping of polymorphic microsatellites in cattle. Animal Genet 24:121.
• Tambasco DD, de Alencar MM, Coutinho LL, Tambasco AJ, Tambasco MD and Regitano LCA (2000) Caracterização molecular de animais da raça Nelore utilizando microssatélites e genes candidatos. Rev Bras Zootecnia 29:1044-49.
• Tautz D and Renz M (1984) Simple sequences are ubiquitous components of eukaryotics genomes. Nucleic Acids Res 12:4127-38.
• Toldo S, Fries SR, Steffen P, Neibergs HL and Barendse W (1993) Physically mapped cosmid-derived microsatellite markers as anchor loci on the bovine chromosome. Mammalian Genome 4:720-27.
• Vaiman D, Mercier D, Moazami-Goudarzi K, Eggen A, Ciampolini R, Lepingle A, Velmala R, Kaukinen J, Varvio SL, Martin P, Levéziel H and Guérin G (1994) A set of 99 cattle microsatellite: characterization, synteny mapping and polymorphism. Mammalian Genome 5:288-297.
• Visscher PM, Woolliams JA, Smith D and Williams JL (2002) Estimation of pedigree errors in the UK dairy population using microsatellite markers and the impact on selection. J Dairy Sci 85:2368-75.

Publication Dates

• Publication in this collection
  01 Sept 2006
• Date of issue
  2006

History

• Received
  11 May 2005
• Accepted
  16 Nov 2005