Front. Vet. Sci.Frontiers in Veterinary ScienceFront. Vet. Sci.2297-1769Frontiers Media S.A.10.3389/fvets.2023.1256928Veterinary ScienceOriginal ResearchIdentification of Anaplasma marginale, Babesia bovis and Babesia bigemina resistance alleles in Crioulo Lageano cattle using PCR-SBT and BoLA-DRB3 gene sequencingCasaMariana da Silva1*GiovambattistaGuillermo2FontequeGraziela Vieira1MiguelEllen Lara3VogelCarla Ivane Ganz4MilettiLuiz Claudio4TakeshimaShin-Nosuke5FontequeJoandes Henrique61Graduate Program in Animal Science (PPGCA), State University of Santa Catarina (UDESC), Lages, Brazil2Facultad de Ciencias Veterinarias UNLP, IGEVET–Institute of Veterinary Genetics, La Plata, Argentina3Scientific Initiation Program, State University of Santa Catarina (UDESC), Lages, Brazil4Department of Animal Production and Food, State University of Santa Catarina (UDESC), Lages, Brazil5Department of Food and Nutrition, Faculty of Human Life, Jumonji University, Niiza, Saitama, Japan6Department of Veterinary Medicine, State University of Santa Catarina (UDESC), Lages, Brazil
Edited by: Albano Beja-Pereira, University of Porto, Portugal
Reviewed by: Shamik Polley, West Bengal University of Animal and Fishery Sciences, India; Teshale Sori Tolera, Addis Ababa University, Ethiopia; Abdulaziz Alouffi, King Abdulaziz City for Science and Technology, Saudi Arabia
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
Introduction
The BoLA-DRB3 gene in cattle is associated with tolerance to several infectious diseases, such as neosporosis, dermatophilosis, leukosis, and mastitis.
Methods
This study used PCR-SBT and BoLA-DRB3 gene sequencing to determine the association between the presence or absence of Anaplasma marginale, Babesia bovis, and Babesia bigemina infections in 208 Crioulo Lageano cattle and alleles present in the population. The chi-square test and odds ratio analysis were employed to establish the association.
Results
Of the BoLA-DRB3 gene alleles present in the population, two alleles were significantly associated with resistance to A. marginale infections: BoLA-DRB3001:01 (p < 0.001; OR = 0.224), which had a frequency of 7.93%, and BoLA-DRB3024:06 (p = 0.007; OR < 0.00001), which had a frequency of 0.72%. Regarding B. bovis infection, the BoLA-DRB3*011:01 allele (p = 0.002; OR = 0.271) had a frequency of 6% in the population and was associated with resistance to the infection. None of the alleles was associated with resistance to infection by B. bigemina.
Discussion
The Crioulo Lageano breed has alleles that may confer resistance against infection by A. marginale and B. bovis.
Locally adapted cattle breeds are crucial for sustainability and food security. Their genetic diversity is especially important in light of environmental challenges like climate change, as it can help improve production characteristics and meet consumer demands (1). Additionally, thanks to such genetic variability, these breeds offer greater resistance to various infectious diseases (2).
Anaplasmosis, caused by Anaplasma marginale, is among the most prevalent diseases in cattle and causes enormous economic losses due to the high morbidity and mortality and costs of its control (3). Similarly, babesiosis, caused by Babesia sp., leads to high morbidity and mortality when clinically manifested, as well as loss of productive efficiency, generating great economic losses (4). It is estimated that in Brazil these two diseases caused by Rhipicephalus microplus (Canestrini, 1888) ticks causes losses of around 3.24 million dollars (5).
The high costs of treatments and the loss of animals affected by anaplasmosis and babesiosis require new alternatives to control these diseases, and one of the most promising alternatives is the selection of animals with genotypes of resistance to clinical diseases (6, 7). These animals are important sources of genetic material for the improvement of herds susceptible to hemoparasites. In this regard, several studies have proven the importance of the BoLA-DRB3 gene with resistance to various infectious diseases (8–10).
The major histocompatibility complex (MHC) is characterized by an intimate relationship between genes that play a role in resistance/susceptibility to infectious diseases, in addition to acting on the innate and adaptive immune response (11). In cattle, MHC (BoLA – bovine leukocyte antigen) is located on chromosome 23 (12) and its genes are grouped into three classes, I, II, and III, according to their function (13). Class II genes are the most studied and are divided into two separate regions of the chromosome. These two parts are known as IIa and IIb (12, 14).
Class IIa has two sub-regions, one composed of the DQ family genes (DQA and DQB) and the other of the DR family genes (DRA and DRB) (15). At least three genes of the DRB family are described in the bovine species: DRB1, which is considered a pseudogene; DRB2, which is little expressed; and DRB3, which is described as highly expressed and polymorphic (16) and associated with tolerance to several diseases, such as neosporosis (17), dermatophilosis (18), leukosis (19, 20), and mastitis (21).
The Crioulo Lageano breed has already accumulated almost five centuries of natural selection and presents characteristics related to resistance to diseases and parasites (22). The challenge of extreme climate conditions during the rigorous winters of the South region of Brazil, the lack of food, predators, and abandonment in vast rural areas are believed to have contributed to the formation of a breed adapted to the environmental conditions of regions of altitude and the subtropical climate (23).
Thus, the presence of resistance alleles to babesiosis and anaplasmosis in locally adapted breeds, such as the Crioulo Lageano breed, associated with the productive potential, demonstrates their importance as a source of genetic material for other herds. Thus, if the resistance of animals to hemoparasites can be increased, reducing treatment costs, loss of animals, and the presence of drug residues in products of animal origin, ensuring sustainable production from locally adapted breeds. For this reason, the aim of this work was to determine resistance BoLA-DRB3 alelles for anaplasmosis and babesiosis infections.
Materials and methodsAnimals and sample acquisition
A total of 208 DNA samples from Crioulo Lageano cattle, including young and adult animals, males and females, from all available categories (bulls, cows, calves, and heifers), from in situ breeding properties, all located in the mountainous region of Santa Catarina, South Brazil, with animals exposed to the same climatic and pasture conditions were used. This breed originated from Bos taurus lineages have high meat quality and is very rustic with great adaptability to hostily conditions (24).
Blood samples from which DNA was extracted were collected for a previous study approved by CETEA Protocol No. 2461171115.
Physical examination
A physical examination was performed to verify clinical signs compatible with clinical diseases. Heart rate (HR), respiratory rate (RR), ruminal movements (RM), rectal temperature, and mucosal color were measured.
Molecular analysis
First a single multiplex-PCR techniques were used to amplify DNA from A. marginale, B. bovis, and B. bigemina. This reaction included for A. marginale, a set of primers based on the MSP5 gene (MSP5 F: 5’-CGC AGA TCT AGC AAA ATC GGC GAG AGG TTT ACC ACT TC-3’e MSP5 R: 5’-GCG CTG CAG TGG CGC AAA ATG CCC GAC ATA CC -3′), whereas for B. bovis and B. bigemina were those described by Figueroa et al. (25) (BoF: 5 -CAC GAG GAA GGA ACT ACC GAT GTT GA-3 and BoR: 5 - CCA AGG AGC TTC AAC GTA CGA GGT CA-3, for B. bovis and BiIA: 5’-CAT CTA ATT TCT CTC CAT ACC CCT CC-3′ e BiIB: 5’-CCT CGG CTT CAA CTC TGA TGC CAA AG −3′). Then a second nested-PCR (n-PCR) was performed for B. bovis and B. bigemina according to Figueroa et al. (25) BoFN 5 - TCA ACA AGG TAC TCT ATA TGG CTA CC-3 e BoRN 5 -CTA CCG AGC AGA ACC TTC TTC ACC AT-3 for B. bovis and (BiIAN: 5’-CGC AAG CCC AGC ACG CCC CGG TGC-3′ e BiIBN: 5’-CCG ACC TGG ATA GGC TGT GTG ATG-3′). The Multiplex-PCR reaction was performed in a final volume of 25 μL of solution with 1 U of Taq Polymerase enzyme GoTaq™ Hot Start Polymerase (Promega, WI, United States), 8.5 pmoles of each primer for each agent, 0.2 mM nucleotides (dNTPs), 3.5 mM magnesium chloride, 5 μL 5X Green GoTaq™ Flexi Buffer (Promega, WI, USA), and 3 μL of DNA (concentration between 20 and 100 ng/μL). Positive controls, using DNA samples donated by Empresa Brasileira de Pesquisa Agropecuária (Embrapa) from experimentally infected animals for each agent and negative controls (that uses water instead of DNA), were used for each reaction.
The temperature conditions applied in the thermocycler (BIOCYCLER™, USA) for the two reactions involved initial denaturation at 95°C for 2 min, followed by 30 cycles of 94°C for 1 min, 54.2°C for 1 min, and 73°C for 1 min, and also a final extension at 73°C for 7 min.
The electrophoresis of amplification products was carried out in a horizontal tank, on a 2% agarose gel added with Unisafe Dye 20,000x (Uniscience, United States). A 100-bp molecular weight marker was used as a standard to determine the size of the sample bands. The electrical source conditions were 140 V and 400 mA for 1 h, with visualization by exposure to ultraviolet light. Bands with a size close to 458 bp for A. marginale, 350 bp for B. bovis, and 278 bp for B. bigemina were considered positive in the first reaction, and bands with a size of approximately 290 bp for B. bovis and 170 bp for B. bigemina were considered positive in the second reaction (Table 1).
Sizes of bands in agarose gel 2% for each agent in first and second PCR reactions.
Size of bands on 2% agarose gel
Agent
Anaplasma marginale
Bebisia bovis
Bebisia bigemina
First Mulptiplex-PCR reaction
458pb
350pb
278pb
Nested-PCR reaction
-
290pb
170pb
BoLA-DRB-3.2 gene genotyping
Amplification of the BoLA-DRB3 gene second exon was performed through single-step PCR, with primers HLO30 (5′-ATC CTC TCT CTG CAG CAC ATT TCC-3′) and HLO32 (5′-TCG CCG CTG CAC AGT GAA ACT CTC-3′) (26, 27).
The reaction conditions consisted of initial denaturation at 95°C for 5 min, followed by 35 cycles of 95°C for 1 min, 65°C for 1 min, and 72°C for 1 min, and a final extension at 72°C for 10 min. The electrophoresis of amplification products was carried out on a 1.5% agarose gel added with Unisafe Dye 20,000x (Uniscience, United States). A 100 bp molecular weight marker was used as a standard to determine the size of the sample bands. The electrical source conditions were 100 V and 400 mA for 40 min and visualization by exposure to ultraviolet light. The obtained fragment had a size of 284 bp from exon 2 of the BoLA-DRB3 gene. The products of this PCR reaction were then sent to a specialized company (ACTGene – Rio Grande do Sul, Brasil) for sequencing by Sanger method, for the exact identification of alleles. Raw sequence data were analyzed using Assign 400ATF ver. 1.0.2.41 software (Conexio Genomics, Fremantle, Australia) as described by Takeshima et al. (28), and using a database of all BoLA-DRB3 alleles described so far in the IPD-MHC.1 Using this information, the software identify the genotype that explain the multiple polymorphic sites within the BoLA-DRB3 second exon.
Statistical analysis
Allele frequency and number of alleles for BoLA-DRB3 gene was obtained by direct counting. For association analysis, animals were grouped in positive and negative for each pathogen, considering that all animals were equaly exposed. Then, the association between each allele with the infection by each pathogen was determined using the chi-square test (p < 0.05), followed by the odds ratio analysis, with SigmaPlot12 software (SigmaPlot version 12.0 for Windows, Systat Software Inc., San Jose, CA, United States). Alleles with OR < 1 are considered resistance alleles and those with OR > 1 are considered susceptibility alleles.
Results
The clinical examination of the samples showed that no Crioulo Lageano animal had clinical signs compatible with anaplasmosis or babesiosis. Table 2 shows the proportions of occurrence of infections for the 208 Crioulo Lageano animals based on the molecular diagnosis performed for the agents A. marginale, B. bovis, and B. bigemina.
Absolute and relative proportion of Crioulo Lageano animals positive and negative for the agents Anaplasma marginale, Babesia bovis, and Babesia bigemina.
Anaplasma marginale
Babesia bovis
Babesia bigemina
Proportion
Positive
Negative
Positive
Negative
Positive
Negative
Absolute
165
43
144
64
127
81
Relative
79.3%
20.7%
69.2%
30.8%
61.0%
39.0%
The association analysis results based on the chi-square test and the odds ratio were summarized in Table 3. These analysis showed that among the alleles of the BoLA-DRB3 gene found in the population, BoLA-DRB3*001:01 (p < 0.001; OR = 0.224), with a frequency of 7.93% in the population, and BoLA-DRB3*024:06 (p = 0.007; OR < 0.00001), with a frequency of 0.72%, are significantly associated with resistance to A. marginale infection in the Crioulo Lageano breed.
Chi-square (Χ2) and Odds ratio analysis between BoLA-DRB3 alleles and presence/absence of the agents Anaplasma marginale, Babesia bovis, and Babesia bigemina, to determine resistance/susceptibility alleles for each agent in the Crioulo Lageano cattle.
BoLA-DRB3 allele
Anaplasma marginale
Babesia bovis
Babesia bigemina
Χ2 (value of p)
Odds ratio
Χ2 (value of p)
Odds ratio
Χ2 (value of p)
Odds ratio
DRB3*1801
0.152 (0.697)
1.269
3.739 (0.053)
0.513
0.017 (0.895)
0.991
DRB3*2801
0.477 (0.490)
0.705
0.088 (0.766)
0.841
0.889 (0.346)
0.684
DRB3*2201
0.325 (0.569)
1.432
0.653 (0.419)
1.480
3.419 (0.064)
2.191
DRB3*0101
12.115 (<0.001)
0.224
2.919 (0.088)
0.502
0.017 (0.896)
0.980
DRB3*1101
0.460 (0.498)
0.651
9.259 (0.002)
0.271
1.368 (0.242)
0.568
DRB3*1501
0.442 (0.506)
1.785
1.435 (0.231)
2.190
0.411 (0.522)
1.489
DRB3*5702
0.061 (0.804)
0.976
0.016 (0.898)
0.950
0.230 (0.632)
1.390
DRB3*0501
0.156 (0.693)
0.724
0.120 (0.729)
1.324
0.461 (0.497)
0.681
DRB3*1601
0.007 (0.930)
0.825
0.906 (0.341)
1.945
0.369 (0.544)
0.688
DRB3*0201
1.474 (0.225)
0.460
0.862 (0.353)
2.129
0.003 (0.951)
0.907
DRB3*0701
3.125 (0.077)
-
1.792 (0.181)
3.219
0.019 (0.888)
0.813
DRB3*4401
3.125 (0.077)
-
0.618 (0.432)
1.970
0.128 (0.720)
1.403
DRB3*2601
1.081 (0.298)
3.766
0.004 (0.948)
1.231
0.033 (0.854)
0.955
DRB3*1701
1.792 (0.181)
-
0.005 (0.939)
1.190
0.581 (0.446)
0.522
DRB3*2406
7.235 (0.007)
<0.00001
0.525 (0.469)
0.220
0.631 (0.427)
-
The bold values means that is a significant association.
The BoLA-DRB3*011:01 allele (p = 0.002; OR = 0.271), with a frequency of 6% in the population, is associated with resistance to B. bovis infection. None of the alleles was associated with resistance or susceptibility to infection by B. bigemina.
Discussion
This is the first work on the association of BoLA-DRB3 alleles with infection by A. marginale, B. bovis, and B. bigemina in cattle in Brazil, as well as in the Crioulo Lageano breed.
In Creole cattle breeds, association studies between BoLA-DRB3 alleles and infectious deaseases were only carried out in the Hartón del Valle breed from Colombia using the PCR-SBT method, reporting a possible resistance alleles to leukosis (29, 30).
In the state of Santa Catarina, the highlands region has already been classified as being of enzootic instability for anaplasmosis, with a prevalence of 27.24% of positive samples in the general herd (31), including the occurrence of an outbreak of the disease in this region (32). The prevalence obtained in breeding properties of the Crioulo Lageano breed was 79.74% and the enzootic instability was not confirmed (33).
Variations in the prevalence of the disease may be associated with the time of sample collection, as the period in which the blood collections were carried out in this study, between summer and autumn, presents a higher hatching of tick larvae and, consequently, more parasitized cattle (34), favoring infection by A. marginale, B. bovis and B. bigemina. For the present study, the proportion of animals infected with A. marginale was 79.3%, with no animal showing clinical signs of the disease (Table 2).
The absence of regular control of vectors and the extensive rearing of animals favors constant contact with hemoparasites, reinforcing their immunity against the causative agent of the clinical disease, which would not necessarily be linked to the presence of certain alleles of the BoLA-DRB3 gene in the population since none of the sampled animals showed clinical signs.
Purebred B. taurus animals are more susceptible to clinical disease caused by A. marginale compared to animals derived from crosses between B. taurus and B. indicus and purebred B. indicus animals (35). This fact was not confirmed for Crioulo Lageano cattle, belonging to the B. taurus lineage. In addition to greater genetic diversity when compared to commercial breeds, native breeds are adapted to the specific environmental conditions of their region after than more than five centuries of natural selection and, in general, there are reports of great resistance to local pathogens (22). In Argentina, Argentine Creole cattle, also Bos taurus, demonstrated phenotypic resistance to Riphicephalus (Boophilus) microplus tick infestation in experimental infection (36). Tick resistance may also be a factor that can influence resistance to hemoparasites.
Regarding babesiosis, the data obtained for the state of Santa Catarina indicated that the prevalence of B. bovis in the general herd is 29.57%, while B. bigemina presented a prevalence of 16.73% (31). The occurrence obtained for these agents in this study was much higher for the Crioulo Lageano breed (69.2% for B. bovis and 61% for B. bigemina) but none of the animals showed clinical signs of the diseases. Despite the use of the same methodology of this study, this difference may have been due to the heterogeneity of the herd analyzed in the work by Vieira et al. (31), which, when analyzing the general herd in the state, obtained samples of B. taurus and B. indicus animals. Furthermore, this variation may also be associated with the time of sample collection, since between summer and autumn, the period in which blood samples were collected for this study, there is a greater hatching of tick larvae and, consequently, more parasitized cattle (34).
Importantly, not all the genes involved in the expression of resistance to babesiosis are known, but it is known that it is a hereditary characteristic (37) and that B. taurus indicus animals tend to be more resistant to babesiosis and their vectors than B. taurus taurus animals (38, 39). However, Crioulo Lageano animals, even those belonging to the B. taurus taurus lineage, do not develop the clinical disease although a large number of animals carry the agent.
In terms of possible genetic resistance to these diseases, the BoLA-DRB3*011:01 allele is frequently found in Creole breeds and has already been associated with resistance to B. bigemina infection (40), which was not confirmed in the Crioulo Lageano breed although it is an allele with high frequency in the population. The difference in the methodology used to identify the alleles, as well as the proportion of infected animals in this study, may justify the association or not of the allele with the disease.
Duangjinda et al. (26) identified the following resistance alleles after genotyping by the PCR-RFLP method: the *14 and *41 alleles for A. marginale; the *14 allele for B. bovis; and the *10 and *51 alleles for B. bigemina. We reiterate the fact that PCR-RFLP has lower precision than the PCR-SBT method and some patterns found may not correspond exactly to the expected allele. In addition, different BoLA-DRB3 alleles share the same PCR-RFLP defined variants. For this reason, results can not be compared.
An association could be established between the alleles BoLA-DRB3*001:01 and BoLA-DRB3*024:06 and A. marginale infection in the Crioulo Lageano breed. Similarly, the BoLA-DRB3*011:01 allele was associated with resistance to B. bovis infection. None of the alleles was associated with resistance to infection by the agent B. bigemina. The BoLA-DRB3*001:01 allele has been associated with mastitis resistance in Holstein cattle (41), while the BoLA-DRB3*011:01 allele has been associated with resistance to B. bigemina and also to resistance to bovine leukosis in the Hartón del Valle breed (29, 40), in addition to being associated with resistance to mastitis (41). The BoLA-DRB3*024:06 allele has previously been described in African breeds, including Sudanese breeds (42) but has not yet been associated with any type of disease.
Research has shown that Crioulo Lageano cattle are more resistant to botflies and ticks than Aberdeen Angus cattle (43). These findings suggest that factors beyond genetics play a role in determining susceptibility to diseases like anaplasmosis and babesiosis. Other factors, such as lineage, contact with vectors, and environmental adaptation, must also be considered alongside genetics.
In addition to the resistance of animals to the infection itself by the agents of anaplasmosis and babesiosis, resistance to the disease vector may also be occurring. BoLA-DRB3 alleles have already been associated with resistance to the tick Rhipicephalus (Boophilus) microplus (44), which would directly influence the acquisition of the infection by the animals.
Data referring to resistance alleles can help animal breeding programs. However, regarding the hemoparasites explored in this study, it is important to maintain susceptible individuals in herds to keep the stability of agents in the population and the heterozygosity of individuals, which is essential to increase the immune response to different types of infection.
Conclusion
The presence of alleles associated with resistance to A. marginale and B. bovis in Crioulo Lageano cattle highlights the importance of maintaining herd variability and heterozygosity. Further research should investigate the relationship between genetic findings and environmental and management factors that contribute to resistance to infectious disease in Crioulo Lageano animals.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The animal studies were approved by CETEA - Comitê de Ética em Experimentação Animal Universidade do Estado de Santa Catarina–UDESC. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent was obtained from the owners for the participation of their animals in this study.
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was supported by the Fundação de Amparo à Pesquisa e Inovação do Estado de Santa Catarina (Foundation for Research and Innovation Support of Santa Catarina State; FAPESC – TERM OF GRANT SPONTANEOUS DEMAND – PESQUISA 2015; No.: 2015TR1543; PROCESSO No.: FAPESC1827/2015) and the Graduate Support Program, partially financed by the Coordination for the Improvement of Higher Education Personnel – Brazil (CAPES – Finance Code 001) for providing postgraduate fellowships and financing the translation service by the Graduate Support Program (PROAP) and Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), chamada CNPq Nº 4/2021 - Bolsas de Produtividade em Pesquisa - PQ - 2021 (315577/2021-5).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
We thank Luciana Gatto Brito from Embrapa Eastern Amazon for providing positive controls for B. bovis and primers for PCR and the Associação Brasileira dos Criadores da Raça Crioula Lageana (ABCCL) for providing the animals for this experiment. We would like to thank Dra. Maria do Socorro Maués Albuquerque from Empresa Brasileira de Pesquisa Agropecuária - Embrapa Recursos Genéticos e Biotecnologia (Embrapa - DF).
ReferencesFioravantiMSSilvaTMMouraMICostaGLDiasJMGuimarãesLKP. Resistance and resilience to diseases in local ruminant breeds: a focus on South America. Arch. Zootec. (2020) 69:338–52.SpringbettAJMacKenzieKWoolliamsJABishopSC. The contribution of genetic diversity to the spread of infectious diseases in livestock populations. Genetics. (2003) 165:1465–74. doi: 10.1093/genetics/165.3.1465, PMID: 14668395BrownCD. Dynamics and impact of tick-borne diseases of cattle. Trop Anim Health Prod. (1997) 29:1S–3S. doi: 10.1007/BF02632905, PMID: 9512735MahmmodYS. Molecular detection of natural Babesia bovis infection from clinically infected and apparently healthy water buffaloes (Bubalus bubalis) and crossbred cattle. J Buffalo Sci. (2012) 1:55–60. doi: 10.6000/1927-520X.2012.01.01.10GrisiLLeiteRCMartinsJRDBarrosATMAndreottiRCançadoPHD. Reassessment of the potential economic impact of cattle parasites in Brazil. Ver Bras de Parasitol Vet. (2014) 23:150–6. doi: 10.1590/S1984-29612014042, PMID: 25054492SchrootenCBovenhuisHVan ArendonkJAMBijmaP. Genetic progress in multistage dairy cattle breeding schemes using genetic markers. J Dairy Sci. (2005) 88:1569–81. doi: 10.3168/jds.S0022-0302(05)72826-5, PMID: 15778327StellaALohuisMMPagnaccoGJansenGB. Strategies for continual application of marker-assisted selection in an open nucleus population. J Dairy Sci. (2002) 85:2358–67. doi: 10.3168/jds.S0022-0302(02)74316-6, PMID: 12362469KulbergSHeringstadBGuttersrudOAOlsakerI. Study on the association of BoLA-DRB3. 2 alleles with clinical mastitis in Norwegian red cows. J Anim Breed Genet. (2007) 124:201–7. doi: 10.1111/j.1439-0388.2007.00662.x, PMID: 17651322RuppRHernandezAMallardBA. Association of bovine leukocyte antigen (BoLA) DRB3.2 with immune response, mastitis, and production and type traits in Canadian Holsteins. J Dairy Sci. (2007) 90:1029–38. doi: 10.3168/jds.S0022-0302(07)71589-8, PMID: 17235182SharifSMallardBAWilkieBNSargeantJMScottHMDekkersJCM. Associations of the bovine major histocompatibility complex DRB3 (BoLA-DRB3) alleles with occurrence of disease and milk somatic cell score in Canadian dairy cattle. Anim Genet. (1998) 29:185–93. doi: 10.1111/j.1365-2052.1998.00318.x, PMID: 9720177PennDJ. Major histocompatibility complex (MHC). Encyclopedia of Life Sciences, pp. 1–7. (2002).StafuzzaNBGrecoAJGrantJRAbbeyCAGillCARaudseppT. A high-resolution radiation hybrid map of the river buffalo major histocompatibility complex and comparison with BoLA. Anim Genet. (2013) 44:369–76. doi: 10.1111/age.12015, PMID: 23216319MosaferJHeydarpourMManshadERussellGSulimovaGE. Distribution of BoLA-DRB3 allelic frequencies and identification of two new alleles in Iranian buffalo breed. Sci World J. (2012) 2012:1–6. doi: 10.1100/2012/863024, PMID: 22454612IannuzziLGallagherDSWomackJEMeoGDSkowLCFerraraL. Chromosomal localization of the major histocompatibility complex in cattle and river buffalo by fluorescent in situ hybridization. Hereditas. (1993) 118:187–90. doi: 10.1111/j.1601-5223.1993.t01-1-00187.x, PMID: 8335474DongxiaoSYuanZ. Note: polymorphisms of the second exon of MHC-DRB gene in Chinese local sheep and goat. Biochem Genet. (2004) 42:385–90. doi: 10.1023/B:BIGI.0000039812.46203.89, PMID: 15524315KumarSSangwanMLAhlawatSBarwarA. Polymorphism in DRB3 exon 2 by PCR-RFLP and its association with mastitis in Murrah buffaloes. Indian J Biotechnol. (2011) 10:232–4.SchwabAEGearyTGBaillargeonPSchwabAJFecteauG. Association of BoLA DRB3 and DQA1 alleles with susceptibly to Neospora caninum and reproductive outcome in Quebec Holstein cattle. Vet Parasitol. (2009) 165:136–40. doi: 10.1016/j.vetpar.2009.07.004, PMID: 19646818MaillardJCMartinezDBensaidA. An amino acid sequence coded by the exon 2 of the BoLA DRB3 gene associated with a BoLA class I specificity constitutes a likely genetic marker of resistance to dermatophilosis in Brahman zebu cattle of Martinique (FWI). Annals of the New York Academy of Sciences. (1996) 791:185–197.LoCWBorjiginLSaitoSFukunagaKSaitouEOkazakiK. BoLA-DRB3 polymorphism is associated with differential susceptibility to bovine leukemia virus-induced lymphoma and proviral load. Viruses. (2020) 12:352. doi: 10.3390/v12030352, PMID: 32235771LoCWTakeshimaSNOkadaKSaitouEFujitaTMatsumotoY. Association of bovine leukemia virus-induced lymphoma with BoLA-DRB3 polymorphisms at DNA, amino acid, and binding pocket property levels. Pathogens. (2021) 10:437. doi: 10.3390/pathogens10040437, PMID: 33917549IbrahimEAAllamNAKotbEEEl-RafeyGAEl-DeenMMAFadlallahMG. Sequence-based typing-study on the relationship between subclinical mastitis and BoLA-DRB3. 2* allelic polymorphism in Egyptian cows. Glob Vet. (2012) 9:8–22.MarianteADSAlbuquerqueMEgitoAAPaivaSRCastroSTR. Conservação de raças brasileiras ameaçadas de extinção e a importância de sua inserção em sistemas de produção. Agrociencia. (2005) 9:459–64.CamargoMARMartinsVMV. Raça bovina Crioula Lageana, um patrimônio genético. Hora Vet. (2005) 24:61–4.MartinsVMV. Raça Crioula Lageana: o esteio do ontem, o labor do hoje a oportunidade do amanhã. Associação Brasileira de Criadores da Raça Crioula Lageana. (2009).FigueroaJVChievesLPJohnsonGSBueningGM. Multiplex polymerase chain reaction based assay for the detection of Babesia bigemina, Babesia bovis and Anaplasma marginale DNA in bovine blood. Vet Parasitol. (1993) 50:69–81. doi: 10.1016/0304-4017(93)90008-B, PMID: 8291198DuangjindaMJindatajakYTipvongWSriwarothaiJPattarajindaVKatawatinS. Association of BoLA-DRB3 alleles with tick-borne disease tolerance in dairy cattle in a tropical environment. Vet Parasitol. (2013) 196:314–20. doi: 10.1016/j.vetpar.2013.03.005, PMID: 23571121SigurdardottirSBorschCGustafssonKAnderssonL. Cloning and sequence analysis of 14 DRB alleles of the bovine major histocompatibility complex by using the polymerase chain reaction. Anim Genet. (1991) 22:199–209. doi: 10.1111/j.1365-2052.1991.tb00670.x, PMID: 1928826TakeshimaSNMatsumotoYMiyasakaTArainga-RamirezMSaitoHOnumaM. A new method for typing bovine major histocompatibility complex class II DRB3 alleles by combining two established PCR sequence-based techniques. Tissue Antigens. (2011) 78:208–13. doi: 10.1111/j.1399-0039.2011.01708.x, PMID: 21623735HernándezDYÁlvarezLAMuñozJE. Evaluación de la resistencia genética del ganado criollo Hartón del Valle al virus de la leucosis bovina. AICA. (2011) 2011:169–72.HernándezHDYMuñozFJEÁlvarezFLA. Asociación del locus BOLA-DRB3. 2 con el virus de la leucosis bovina en el ganado criollo colombiano. Rev Colombiana Cienc Anim. (2014) 6:319–26. doi: 10.24188/recia.v6.n2.2014.435VieiraLLCaneverMFCardozoLLCardosoCPHerkenhoffMENetoAT. Prevalence of Anaplasma marginale, Babesia bovis, and Babesia bigemina in cattle in the Campos de Lages region, Santa Catarina state, Brazil, estimated by multiplex-PCR. Parasite Epidemiol Control. (2019) 6:e00114. doi: 10.1016/j.parepi.2019.e00114, PMID: 31528739CaneverMFVieiraLLReckCRichterLMilettiLC. First evaluation of an outbreak of bovine babesiosis and anaplasmosis in southern Brazil using multiplex PCR. Korean J Parasitol. (2014) 52:507–11. doi: 10.3347/kjp.2014.52.5.507, PMID: 25352699CasaMSVettoriJMde SouzaKMMilettiLCVogelCIGLimaALF. High prevalence of Anaplasma marginale in the Crioula Lageana cattle. J Infect Dev Ctries. (2020) 14:623–30. doi: 10.3855/jidc.11691, PMID: 32683353SouzaAPGonzalesJCRamosCIPaloschiCGde MoraesAN. Variação sazonal de Boophilus microplus no Planalto Catarinense. Pesq Agrop Bras. (1988) 23:627–30.BockREDe VosAJKingstonTGMcLellanDJ. Effect of breed of cattle on innate resistance to infection with Babesia bovis, B bigemina and Anaplasma marginale. Aust Vet J. (1997) 75:337–40. doi: 10.1111/j.1751-0813.1997.tb15706.x, PMID: 9196820OrtegaMFGiovambattistaGCutulléCSantosDDNavaSBonamyM. Phenotypic evaluation of genetic resistance to the tick Rhipicephalus (Boophilus) microplus in argentine creole cattle. Ticks Tick Borne Dis. (2023) 14:102223–7. doi: 10.1016/j.ttbdis.2023.102223BilhassiTBOliveiraHNIbelliAMGigliotiRRegitanoLCOliveira-SequeiraTC. Quantitative study of Babesia bovis infection in beef cattle from São Paulo state. Brazil Ticks Tick Borne Dis. (2014) 5:234–8. doi: 10.1016/j.ttbdis.2013.11.002, PMID: 24522252JonssonNN. The productivity effects of cattle tick (Boophilus microplus) infestation on cattle, with particular reference to Bos indicus cattle and their crosses. Vet Parasitol. (2006) 137:1–10. doi: 10.1016/j.vetpar.2006.01.010, PMID: 16472920PiperEKJacksonLABielefeldt-OhmannHGondroCLew-TaborAEJonssonNN. Tick-susceptible Bos taurus cattle display an increased cellular response at the site of larval Rhipicephalus (Boophilus) microplus attachment, compared with tick-resistant Bos indicus cattle. Int J Parasitol. (2010) 40:431–41. doi: 10.1016/j.ijpara.2009.09.009, PMID: 19852965BolañosIHernandezDÁlvarezL. Asociación de los alelos del gen BoLA-DRB3 con la infección natural de Babesia spp en el ganado criollo Hartón del Valle. Arch Zootec. (2017) 66:113–20. doi: 10.21071/az.v66i253.2133YoshidaTFurutaHKondoYMukoyamaH. Association of BoLA-DRB3 alleles with mastitis resistance and susceptibility in Japanese Holstein cows. Anim Sci J. (2012) 83:359–66. doi: 10.1111/j.1740-0929.2011.00972.x, PMID: 22574787SalimBTakeshimaSNNakaoRMoustafaMAAhmedMKAKambalS. BoLA-DRB3 gene haplotypes show divergence in native Sudanese cattle from taurine and indicine breeds. Sci Rep. (2021) 11:17202. doi: 10.1038/s41598-021-96330-7, PMID: 34433838CardosoCPSilvaBFGonçalvesDSTagliariNJSaitoMEAmaranteAF. Resistence against ectoparasites in CriouloLageano and crossbred Angus cattle in southern Brazil under natural conditions. Pesqui Vet Bras. (2014) 34:141–6. doi: 10.1590/S0100-736X2014000200008MartinezMLMachadoMANascimentoCSSilvaMVGTeodoroRLFurlongJ. Association of BoLA-DRB3. 2 alleles with tick (Boophilus microplus) resistance in cattle. Genet Mol Res. (2006) 5:513–24. PMID: 17117367