Polymorphism of locus DRB3.2 in populations of Creole Cattle from Northern MexicoAbstract
The polymorphism of locus BoLA-DRB3.2 of the Major Histocompatibility Complex was evaluated in two northern Mexican Creole cattle populations, Chihuahua (n = 47) and Tamaulipas (n = 51). The BoLA-DRB3.2 locus was typed by amplification and digestion with restriction endonuclease enzymes (PCR-RFLP). Fifty-two alleles were detected (28 previously reported and 24 new ones). In the Chihuahua population, 18 alleles and 5.5 effective alleles were found, while in the Tamaulipas population there were 34 and 10.8, respectively. The allele frequencies ranged from 0.011 to 0.383 in Chihuahua and from 0.010 to 0.206 in Tamaulipas. The frequencies of the new alleles in both cattle populations were low (0.010 to 0.053). The expected heterozygosity was 0.827 and 0.916, respectively, for the Chihuahua and Tamaulipas populations. Both populations presented a heterozygote deficit: [Chihuahua F IS = 0.1 (p = 0.019) and Tamaulipas F IS = 0.317 (p < 0.001)]. In conclusion, this study showed that the Mexican Creole cattle have many low-frequency alleles, several of which are exclusive to these populations. Genetic distances obtained show that the Mexican Creole cattle population is composed of independent populations, far apart from other South American Creole populations.
BoLA-DRB3.2; polymorphism; Creole cattle; PCR-RFLP
GENETICS
Polymorphism of locus DRB3.2 in populations of Creole Cattle from Northern Mexico
Ilda G. Fernández1,2
José Gonzalo Ríos Ramírez2
Amanda Gayosso Vázquez3
Raúl Ulloa Arvizu3
Rogelio A. Alonso Morales3
1Departamento de Ciencias Médico Veterinarias, Universidad Autónoma Agraria Antonio Narro, Torreón, Coahuila, México
2Departamento de Reproducción, Facultad de Zootecnia, Universidad Autónoma de Chihuahua, Chihuahua, Chihuahua, México
3Departamento de Genética y Bioestadística, Facultad de Medicina Veterinaria y Zootecnia, Universidad Nacional Autónoma de México, México, D.F., México
The polymorphism of locus BoLA-DRB3.2 of the Major Histocompatibility Complex was evaluated in two northern Mexican Creole cattle populations, Chihuahua (n = 47) and Tamaulipas (n = 51). The BoLA-DRB3.2 locus was typed by amplification and digestion with restriction endonuclease enzymes (PCR-RFLP). Fifty-two alleles were detected (28 previously reported and 24 new ones). In the Chihuahua population, 18 alleles and 5.5 effective alleles were found, while in the Tamaulipas population there were 34 and 10.8, respectively. The allele frequencies ranged from 0.011 to 0.383 in Chihuahua and from 0.010 to 0.206 in Tamaulipas. The frequencies of the new alleles in both cattle populations were low (0.010 to 0.053). The expected heterozygosity was 0.827 and 0.916, respectively, for the Chihuahua and Tamaulipas populations. Both populations presented a heterozygote deficit: [Chihuahua FIS = 0.1 (p = 0.019) and Tamaulipas FIS = 0.317 (p < 0.001)]. In conclusion, this study showed that the Mexican Creole cattle have many low-frequency alleles, several of which are exclusive to these populations. Genetic distances obtained show that the Mexican Creole cattle population is composed of independent populations, far apart from other South American Creole populations.
The Mexican Creole cattle originated from Bos taurus populations brought to Mexico by the Spaniards in the XVI century, mainly from the Iberian peninsula. Currently, the Mexican Creole cattle population consists of unselected populations with low productivity features, such as low body weight and low milk production. On the other hand, they seem to have the ability of surviving in extreme climates and of resisting to diseases. These cattle are seriously endangered by the rapid displacement due to improved modern breeds, but some populations still remain in remote areas in South, Central and North Mexico. In northern Mexico, Creole cattle are raised by small communities under precarious extensive conditions (low-input and low-management) in a semi-arid climate. In the state of Chihuahua, they are reared by peasants and by indigenous breeders of the Sierra Tarahumara, ever since their original introduction by Jesuit priests in 1667 (Ríos, 2004). In the state of Tamaulipas, similarly, the Creole cattle are raised by small farmers, although the breed was introduced from different regions of the country, mainly by a genetic rescue program that the Mexican government started around 1980.
Certain populations of Mexican Creole cattle have become a biotype of commercial interest for rodeo sporting events in the USA and Canada, because they possess features required for this activity. These features are hardiness, docility, speed, as well as tolerance to physical punishment and orientation, shape, length, thickness, and strength of the horns (Félix Portillo et al., 2006).
One of the characteristics that have probably evolved due to natural selection is the resistance to infectious diseases, because the Creole cattle are usually not immunized. The Major Histocompatibility Complex (MHC) is one of the main components of the immune system. Also called bovine lymphocyte antigen (BoLA) complex in cattle, it is located on chromosome 23 and contains some genes involved in antigen presentation necessary to trigger the immune response (Hess et al., 1999). The class II BoLA-DR region is composed of one DRA locus and at least three DRB loci (DP, DQ, and DR) (Muggli-Cockett and Stone, 1989). The MHC allele diversity is associated with the ability to recognize a large number of antigens, resulting in a more efficient immune response (Behl et al., 2007). So it may be interesting to assess the level of allele diversity in the BoLA-DRB3.2 gene in a population that is under great pressure for survival under tough conditions, such as the Mexican Creole cattle populations. In addition, this polymorphism can be used to study the genetic relationships between populations and to assess their levels of genetic differentiation. Presently, extensive information is available on the levels of genetic diversity of exon 2 of the DRB3 gene in different populations of cattle, obtained by amplification of this segment by PCR, and subsequent digestion with endonucleases (Gilliespie et al., 1999; Kelly et al., 2003; Martinez et al., 2005; Behl et al., 2007). The objective of this study was to investigate the level of genetic diversity present in the MHC BoLA-DRB3.2 locus in two Creole cattle populations from northern Mexico.
Dendrogram depicting the genetic distances among Chihuahua and Tamaulipas Creole cattle, South American Creole cattle and taurine and Zebu breeds.
Samples were obtained from 98 Creole Bos taurus cattle from two northern states of Mexico: Chihuahua State samples were obtained from three localities (Morelos: n = 26, Chinipas: n = 7, and Guachochi: n = 14), all located in Sierra Tarahumara. All of the Tamaulipas samples came from a particular region (Soto la Marina; n = 51).
DNA was purified from blood samples following standard procedures (Sambrook et al., 1989). A 284-bp segment of DNA containing the BoLA-DRB3.2 gene was amplified by PCR. For this purpose, we used a pair of oligonucleotides described previously by van Eijk et al. (1992) (HL030: 5'-ATCCTCTCTCTGCAGCACATTTC C-3' and HL032: 5'-TCGCCGCTGCACAGTGAAACTC TC-3'). The amplification of the BoLA-DRB3.2 gene was performed in a single PCR reaction rather than in two subsequent hemi-nested reactions as originally proposed by van Eijk et al. (1992). PCR reactions were carried out in a total volume of 50 μL, containing 150 ng of genomic DNA, 1.5 mM MgCl2, 10 mM Tris – HCl pH 8.4, 50 mM KCl, 1 mg/mL gelatin, 0.15 mg/mL bovine serum albumin, 0.1% Triton-X100, 0.2 mM deoxynucleotide triphosphates (dNTPs; Biogénica S.A. de C.V. Mexico), 1 μM of each oligonucleotide HL030 and HL032, and 2.5 U Taq DNA polymerase (Biogénica, S. A. de C. V. Mexico). Amplification was carried out in a Thermal Cycler (Omn-E HYBAID, UK), with an initial denaturation step of 3 min at 94 °C, followed by 30 cycles of 30 s at 94 °C, 30 s at 68 °C and 30 s at 72 °C, and a final extension step of 5 min at 72 °C.
Three aliquots of 60 μL were digested independently, each containing 10 μL of the PCR-amplified products, 1X of the recommended digestion buffer, and 5 U of each of the restriction enzymes RsaI, BstYI, and HaeIII (New England Biolabs Inc., Mass., USA), and incubated at the recommended temperature for 14 h. As a control of BstYI complete digestion, 10 μL of the digestion mix were added to 160 ng of pBR322 DNA and incubated in parallel.
The digested products were resolved by vertical gel electrophoresis (8 x 10 cm) on an 8% polyacrylamide gel (Sambrook et al., 1989), and stained with silver nitrate using a commercial kit (BIO-RAD Silver Stain Plus 161-0448). As a molecular weight marker, pBR322 DNA digested with MspI was used. The size of the digestion products was estimated with the RFLP Molecular Weight Program, version 2.17 (UVP, Inc., Upland, CA, USA). Allele identification from the restriction patterns and its nomenclature have been previously reported by van Eijk et al. (1992), Gelhaus et al. (1995), Maillard et al. (1999), and Gilliespie et al. (1999).
The PCR products were cloned into a T vector (pUC19 DNA was digested with EcoRV and T tails added at the 3' end, as described by Sambrook et al., 1989). The cloned products were subjected to PCR and the amplified products digested with restriction enzymes.
Statistical analysis was carried out using the GENEPOP v.3.4 program (Raymond and Rousset, 1995a) to calculate allele frequencies, the value FIS (Weir and Cockerham, 1984) that measures the heterozygote deficit (Rousset and Raymond, 1995), and the exact test for population differentiation (Raymond and Rousset, 1995b). Statistical significance was measured by the Markov Chains method (Guo and Thompson, 1992). The allele frequencies used for comparison and genetic analysis were those published for Argentine Creole cattle (Giovambattista et al., 1996), Uruguay Creole cattle (Kelly et al., 2003), Blancorejinegro Creole of Colombia (Martinez et al., 2005), Saavedreño cattle of Bolivia (Ripoli et al., 2004), Holstein and Jersey (Sharif et al., 1998), and Zebu Gir (da Mota et al., 2002). The genetic distances were calculated according to the method described by Reynolds et al. (1983), using the Gendist program, implemented in the PHYLIP package (Felsenstein, 1993). Genetic distances were used to construct a neighbor-joining phylogenetic tree (Huson and Bryant, 2006). The unbiased expected heterozygosity or gene diversity was calculated from the previously reported frequencies (Nei and Chesser, 1983), and the effective allele number (ne) was obtained as the inverse of the homozygosity (Nei, 1987).
In the 98 animals tested, we identified 52 BoLA-DRB3.2 alleles; 28 had already been reported (van Eijk et al., 1992; Gelhaus et al., 1995; Maillard et al., 1999; Gilliespie et al., 1999; Behl et al., 2007) and 24 were new patterns. The identified alleles and their frequencies are shown in Table 1. In the Chihuahua Creole cattle, we found 18 alleles of which seven were new, and in the Tamaulipas Creole cattle we found 34 alleles, of which 17 were new; these two populations shared seven alleles, of which two were new ones (nbd and sba).
Genetic diversity was 0.83 and 0.92 for Chihuahua and Tamaulipas Creole cattle, respectively. The effective allele number was 5.5 for Chihuahua and 10.8 for Tamaulipas. As shown in Table 1, in the Chihuahua population, the cumulative frequency of six alleles out of 18 (*24, *19, ibb, *15, xbb, and *23) accounted for 75% of the allele frequency. For the population of Tamaulipas, eight alleles out of 34 (*23, *24, *39, *17, *18, *15, tbb, and *22) accounted for 66.6% of the whole allele frequency.
The FIS parameter assesses the deviation between expected and observed heterozygosity (Weir and Cockerham, 1984). In our study, the FIS values obtained for the Chihuahua Creole and Tamaulipas populations were 0.1 (p = 0.0196) and 0.317 (p < 0.001), respectively. Positive values mean that both groups present a heterozygote deficit. The FIS values for the three subpopulations of the Chihuahua population (Chinipas, Guachochi and Morelos) were also estimated (-0.063, 0.233 and 0.075, respectively), and only the Morelos population presented heterozygote deficiency at a significant level (p = 0.034).
The FST value measures the degree of genetic differentiation between populations. This value was FST = 0.043 between the Chihuahua and Tamaulipas groups, which means that 4% of the genetic variation corresponded to differences between populations, whereas the remaining 96% were due to differences among individuals within each population. These results indicate little genetic differentiation between these two populations. We also performed a genetic differentiation test among the subpopulations of Chihuahua (Morelos, Chinipas and Guachochi) and Tamaulipas (FST = 0.033). The allele distribution was similar among the three Chihuahua subpopulations (p > 0.05) and between Chinipas and Tamaulipas (p > 0.05); however, it was different between Guachochi and Morelos vs. Tamaulipas (p < 0.001). In order to analyze the genetic relationships of the Mexican Creole populations with other American Creole cattle and with other commercial breed cattle, we evaluated their genetic distances by the Reynolds method, using the allele frequencies reported for Argentina (Giovambattista et al., 1996), Uruguay (Kelly et al., 2003), Colombia (Martinez et al., 2005), Bolivia (Ripoli et al., 2004), Holstein and Jersey breeds (Sharif et al., 1998), and for the non-taurine population Zebu Gir (da Mota et al., 2002). The genetic relationships among the different populations are shown in Figure 1.
The results of this study show that Chihuahua and Tamaulipas Creole cattle have a high degree of BoLA-DRB3.2 polymorphism, in concordance with findings in other populations, such as the American Creole cattle (Giovambattista et al., 1996; Kelly et al., 2003; Ripoli et al., 2004; Martinez et al., 2005), European taurine breeds (Sharif et al., 1998), Asian taurine breeds (Takeshima et al., 2003), and the Zebu Gir breed (da Mota et al., 2002; Behl et al., 2007). In this study, allele DRB3.2*24 showed the highest frequency (0.38) in Chihuahua Creole cattle, and it has also been reported to be present at frequencies higher than 0.1 in all taurine populations (Giovambattista et al., 1996; Sharif et al., 1998; Gilliespie et al., 1999; Kelly et al., 2003). This allele was further reported to be associated with a higher antibody-mediated immune response and with a lower cell-mediated immune response (Rupp et al., 2007). In the Tamaulipas population, the most frequently found allele was DRB3.2*23 (0.206). This allele has been detected in almost all populations of taurine cattle (Giovambattista et al., 1996; Sharif et al., 1998; Gilliespie et al., 1999), but not in Uruguayan Creole cattle (Kelly et al., 2003). The high frequency of alleles DRB3.2*23 and DRB3.2*24 found in Tamaulipas and Chihuahua Creole cattle may be the result of positive natural selection, but this remains to be demonstrated.
Allele DRB3.2*10 was observed at low frequencies (< 0.1) in most of the groups except in Tamaulipas Creole and Blancorejinegro Creole (Martinez et al., 2005). Allele DRB3.2*16 was more frequent in Uruguayan Creole (Kelly et al., 2003) and Saavedreño (Ripoli et al., 2004), but in our study it was not detected in Tamaulipas Creole and was present in Chihuahua Creole at a low frequency (1.1%). Gelhaus et al. (1995) demonstrated the presence of eight alleles in African breeds (N'Dama and Boran); we found some of these alleles in the Tamaulipas population (DRB3.2*16 and DRB3.2*32) and in the Chihuahua Creole cattle (DRB3.2*6 and DRB3.2*10).
It is known that the MHC genes are highly polymorphic in the populations studied so far (Dietz et al., 1997). Although we found a high number of alleles in Mexican populations, the Chihuahua Creole cattle had the lowest effective allele number (5.5), even lower than that of Holstein (7.4) and Jersey (8.2) populations, which are known to be undergoing intense selection and inbreeding. Actually, in the Mexican Creole cattle many alleles at low frequencies were found, but presenting some exclusiveness. In fact, considering the ratio between the effective allele number and the total number of alleles, the Argentine and Uruguayan Creole showed a high relationship (0.8 and 0.72, respectively), whereas in the Chihuahua and Tamaulipas Creole cattle this relationship was very low (0.3 and 0.317, respectively). This may be the result of natural selection on a particular BoLA allele pool, leading to the development of resistance to local infectious diseases.
Of the novel alleles reported in this study, sequence analysis was made in 10 of them by Félix Portillo et al. (2006), who found that two of these corresponded to the previously reported alleles DRB3*1602 and DRB3*1501. The remaining eight alleles rendered unique nucleic acid sequences after a search in the official BoLA database. Even though the nucleotide sequencing of BoLA genes is a more precise way of allele identification and comparison, for laboratories with limited resources the PCR-RFLP technique is a more useful and effective tool for BoLA-DRB3.2 genotyping, because it is easier to perform, cheaper and faster.
In the Chihuahua and Tamaulipas Creole cattle, we observed an excess of homozygosity that may be due to the combined effect of inbreeding and the Wahlund effect. The Chihuahua Creole cattle animals belonged to three subpopulations from the Sierra Tarahumara, which are isolated and have only few sires available. In turn, the Tamaulipas Creole cattle, although the population descended from animals from different regions of Mexico, are currently confined and consist of small, closed herds.
Finally, Figure 1 presents a dendrogram that summarizes the genetic distances and shows that the Mexican Creole cattle is an independent branch, far apart from the South American Creole populations. This may indicate that this population has a different origin from those populations, or that it has experienced extreme genetic differentiation, probably due to a founder effect. It is also evident that there is extensive differentiation among the different populations of Creole cattle in America, although the Argentinean and Uruguayan populations are closer related, being placed in a common node. Furthermore, although the Chihuahua and Tamaulipas Creole populations are in the same branch, there is a great genetic distance between them. This may be the result of the ancient isolation and independent origin of the Chihuahua group. This information is important for the knowledge of this endangered genetic resource and for the establishment of strategic conservation programs. Furthermore, our findings suggest that these local cattle populations may be a valuable reservoir of genetic diversity for the development of defined local breeds or in the improvement of commercial cattle lines.
Acknowledgments
The authors wish to express their gratitude to all members of the Laboratorio de Genética Molecular, Facultad de Medicina Veterinaria y Zootecnia, Universidad Nacional Autónoma de México, for their help, and particularly to MC María del Carmen Castro Méndez and to Dr. Refugio Cortés Fernández.
References
Behl J.D., Verma N.K., Behl R., Mukesh M., Ahlawat S.P.S. (2007) Characterization of genetic polymorphism of the bovine lymphocyte antigen DRB3.2 locus in Kankrej cattle (Bos indicus). J Dairy Sci 90:2997-3001.
da Mota A.F., Gabriel J.E., Martinez M.L., Coutinho L.L. (2002) Distribution of bovine lymphocyte antigen (BoLA-DRB3) alleles in Brazilian dairy Gir cattle (Bos indicus). Eur J Immunogenet 29:223-227.
Dietz A.B., Cohen N.D., Timms L., Kehrli Jr M.E. (1997) Bovine lymphocyte antigen class II alleles as risk factors for high somatic cell counts in milk of lactating dairy cows. J Dairy Sci 80:406-412.
Félix Portillo M., Ríos Ramírez J.G., Erosa de la Vega G.E., Rodríguez Almeida F. (2006) Sequencing of new BoLA-DRB3.2 alleles detected in Mexican Creole cattle. Téc Pecu Méx 44:15-25.
Felsenstein J. (1993) PHYLIP: Phylogeny Inference Package, v.3.55Seattle.University of Washington
Gelhaus A., Schnittger L., Mehlitz D., Horstmann R.D., Meyer C.G. (1995) Sequence and PCR-RFLP analysis of 14 novel BoLA-DRB3 alleles. Anim Genet 26:147-153.
Gilliespie B.E., Jayarao B.M., Dowlen H.H., Oliver S.P. (1999) Analysis and frequency of bovine lymphocyte antigen DRB3.2 alleles in Jersey cows. J Dairy Sci 82:2049-2053.
Giovambattista G., Golijow C.D., Dulout F.N., Lojo M.M. (1996) Gene frequencies of DRB3.2 locus of Argentine Creole cattle. Anim Genet 27:55
Guo S.W., Thompson E.A. (1992) Performing the exact test of Hardy-Weinberg proportions for multiple alleles. Biometrics 48:361-372.
Hess M., Goldammer T., Gelhaus A., Ried K., Rappold G., Eggen A., Bishop M.D., Schwerin M., Horstmann R.D. (1999) Physical assignment of the bovine MHC class IIa and IIb genes. Cytogenet Cell Genet 85:244-247.
Huson D., Bryant D. (2006) Application of phylogenetic networks in evolutionary studies. Mol Biol Evol 23:254-267.
Kelly L., Nicolini P., D'Angelo M., Nimo A., Rincón G., Piagio J., Postiglioni A. (2003) Polimorfismo de gen DRB3.2 en bovinos criollos del Uruguay. Arch Zootec 52:77-80.
Maillard J.C., Renard C., Chardon P., Chantal I., Bensaid A. (1999) Characterization of 18 new BoLA-DRB3 alleles. Anim Genet 30:200-203.
Martínez R., Toro R., Montoya F., Burbano M., Tobón J., Gallego J., Ariza F. (2005) Caracterización del locus BoLa-DRB3 en ganado criollo colombiano y asociación con resistencia a enfermedades. Arch Zootec 54:349-356.
Muggli-Cockett N.E., Stone R.T. (1989) Partial nucleotide sequence of a bovine major histocompatibility class II DR beta-like gene. Anim Genet 20:361-367.
Nei M. (1987) Molecular Evolutionary GeneticsNew YorkColumbia University Press512
Nei M., Chesser K. (1983) Estimation of fixation indices and gene diversities. Ann Hum Genet 47:253
Raymond M., Rousset F. (1995a) GENEPOP, v. 1.2: Population genetics software for exact tests and ecumenicism. J Hered 86:248-249.
Raymond M., Rousset F. (1995b) An exact test for population differentiation. Evolution 49:1280-1283.
Reynolds J., Weir B.S., Cockerham C.C. (1983) Estimation of the coancestry coefficient: Basis for a short-term genetic distance. Genetics 105:767-779.
Ríos J.G. (2004) Perspectivas de la producción de ganado de rodeo en Chihuahua. Chihuahua Ganadero 29:5
Ripoli M.V., Lirón V.J.P., De Luca J.C., Rojas F., Giovambattista G. (2004) Gene frequency distribution of the BoLA-DRB3 locus in Saavedreño Creole dairy cattle. Biochem Genet 42:231-240.
Rousset F., Raymond M. (1995) Testing heterozygote excess and deficiency. Genetics 140:1413-1419.
Rupp R., Hernandez A., Mallard B.A. (2007) Association of bovine leukocyte antigen (BoLA) DRB3.2 with immune response mastitis, and production and type traits in Canadian Holstein. J Dairy Sci 90:1029-1038.
Sambrook J., Fritsch E.F., Maniatis T. (1989) Molecular Cloning: A Laboratory Manual2nd editionNew York.Cold Spring Harbor Press
Sharif S., Mallard B.A., Wilkie B.N., Sargeant J.M., Scott H.M., Dekkers J.C.M., Leslie K.E. (1998) Associations of the bovine major histocompatibility complex DRB3 (BoLA-DRB3) alleles with occurrence of diseases and milk somatic cell score in Canadian dairy cattle. Anim Genet 29:185-193.
Takeshima S., Saitou N., Morita M., Inoko H., Aida Y. (2003) The diversity of bovine MHC class II DRB3 genes in Japanese Black, Japanese Shorthorn, Jersey and Holstein cattle in Japan. Gene 316:111-118.
van Eijk M.J.T., Stewart-Haynes J.A., Lewin H.A. (1992) Extensive polymorphisms of the BoLA-DRB3 gene distinguished by PCR-RFLP. Anim Genet 23:483-496.
Weir B., Cockerham C. (1984) Estimating F-statistics for the analysis of population structure. Evolution 38:1358-1370.
Internet Resources
Report of the 8th Workshop of the ISAG BoLA Nomenclature Committee (2002) Available at June2,2007 from: http://www.projects.roslin.ac.uk/bola/wk98.html
- Allele distribution and gene frequency of locus DRB3.2 in Creole cattle and in improved breeds.
ne = Effective number of alleles. Alleles 1-53 ISAG (2002). n = Sample size. aChihuahua, bTamaulipas, cBlancorejinegro; Martínez et al. (2005), dGiovambattista et al. (1996), eKelly et al. (2003), fRipoli et al. (2004), gda Mota et al. (2002), hSharif et al. (1998), i*Gilliespie et al. (1999), **Behl et al. (2007).
Received: November 30, 2007; Accepted: June 30, 2008
Ilda G. Fernández. Departamento de Ciencias Médico Veterinarias, Universidad Autónoma Agraria Antonio Narro, Periférico Raúl López Sánchez y Carretera a Santa Fe, 27054 Torreón, Coahuila, México. E-mail: ilda_fernandez_garcia@yahoo.com.mx.
Publication Dates
-
Publication in this collection
19 Nov 2008 -
Date of issue
2008
History
-
Received
3011 -
Accepted
3006
