- 1Key Laboratory of Human Genetics and Environmental Medicine, School of Public Health, Xuzhou Medical University, Xuzhou, China
- 2Laboratory of Zoonoses, China Animal Health and Epidemiology Center, Qingdao, China
- 3Testing Laboratory, Inner Mongolia Animal Disease Control Center, Hohhot, China
- 4College of Animal Science and Technology, Shandong Agriculture University, Taian, China
Brucella BP26 proves to be a highly immunogenic antigen with excellent specificity in brucellosis detection. In China, the authorized use of the Bp26-deleted vaccine M5ΔBP26 for preventing small ruminant brucellosis highlights the importance of developing accurate detection methods targeting BP26, particularly for the diagnosis of differentiation between infected and vaccinated animals (DIVA). Using the traditional mouse hybridoma technique, we successfully obtained 12 monoclonal antibodies (mAbs) targeting BP26. The efficacy of these mAbs in detecting various animal brucellosis cases using the competitive ELISA method was evaluated. Among them, only the E10 mAb exhibited significant efficiency, being inhibited by 100, 97.62, and 100% of brucellosis-positive sera from cattle, small ruminants, and canines, respectively. The E10-based competitive enzyme-linked immunosorbent assay (cELISA) outperformed the BP26-based indirect enzyme-linked immunosorbent assay (iELISA) in accuracy, particularly for cattle and small ruminant brucellosis, with cELISA sensitivity reaching 97.62% compared to 64.29% for iELISA for small ruminants. Although cELISA showed slightly lower specificity than iELISA, it still maintained high accuracy in canine brucellosis detection. The epitope of mAb E10 was identified in the amino acid sequence QPIYVYPDDKNNLKEPTITGY, suggesting its potential as a diagnostic antigen for brucellosis. In conclusion, the E10-based cELISA presents an effective means of detecting animal brucellosis, particularly significant for DIVA diagnosis in China, where the BP26-mutant vaccine is widely used.
1 Introduction
The 26-kDa periplasmic protein of Brucella, BP26, serves as a robust immunogenic antigen capable of eliciting an antibody response in Brucella-infected animals (1–8). Despite a delayed and weaker antibody response compared to Brucella O-polysaccharide (OPS) antigen (3), BP26 exhibits higher specificity in brucellosis detection (9), with minimal cross-reactivity with sera infected with other bacterial pathogens. Additionally, BP26 presents an advantage in detecting brucellosis caused by rough Brucella strains, which lack OPS antigen on their surface (4, 10).
Furthermore, studies have demonstrated that BP26 gene-deleted vaccines maintain their protective efficacy against wild-type Brucella strains, making BP26-deficient vaccines, such as S19 and Rev1, preferable for achieving serological differentiation between infected and vaccinated animals (DIVA) (11–14). In China, the authorization of the BP26-mutant vaccine M5ΔBP26 for preventing small ruminant brucellosis highlights the importance of developing detection methods targeting BP26 in facilitating DIVA strategies with the widespread use of M5ΔBP26 vaccines (15).
We and other researchers have evaluated the efficacy of the BP26-based indirect enzyme-linked immunosorbent assay (iELISA) in serologically detecting various forms of animal brucellosis (7, 16). Fewer studies have explored the competitive enzyme-linked immunosorbent assay (cELISA) using anti-BP26 monoclonal antibodies (mAbs) for brucellosis detection (17, 18). In this study, we used the traditional mouse hybridoma technique to generate a panel of mAbs targeting BP26. Using a collection of brucellosis-positive sera, we assessed the ability of these mAbs to detect different forms of animal brucellosis, identifying one mAb with high efficiency in cELISA for this purpose. We anticipate that our findings will contribute to the development of novel diagnostic methods for brucellosis.
2 Materials and methods
2.1 Ethics statement
The animal study was reviewed and approved by the Experimental Animal Ethics Committee of Xuzhou Medical University (approval number: 202208W051).
2.2 Production of anti-BP26 mAbs
The BP26 protein (reference strain: 16M/ATCC 23456/NCTC 10094) was expressed using a prokaryotic expression system established in our laboratory (16). The immunization of BALB/c mice and hybridoma screening procedures followed methods outlined in a previously published study (18). Anti-BP26 mAbs were generated by intraperitoneal inoculation of mice with hybridoma cells. Mouse ascites containing mAbs were collected and purified using a commercial Protein G column (Sangon Biotech, Shanghai, China). The purity of the mAbs was assessed using SDS-PAGE.
2.3 Sera
All sera utilized in this experiment were archived at the Chinese Animal Epidemiology and Health Center (CAHEC) with well-documented backgrounds. Positive sera were sourced from animals from which wild-type Brucella had been isolated, while negative sera originated from animals raised in brucellosis-free zones. A total of 245 sera were included in the study, comprising 96 from small ruminants (naturally infected: 42; negative: 54), 96 from dairy cattle (naturally infected: 44; negative: 52), and 53 from canines (naturally infected: 7; negative: 46). Additionally, 231 sera from ruminants immunized with the M5ΔBP26 vaccine were collected. These sera were used to assess the efficacy of both the iELISA and cELISA methods.
Furthermore, rabbit sera (purchased from Tianjin Biochip Corporation, Tianjin, China) artificially infected with Y. enterocolitica O9, E. coli (O157:H7, O116), Salmonella urban, Ochrobactrum anthropi, and Vibrio cholerae were utilized to assess the specificity of both the iELISA and cELISA methods. Rabbit sera underwent the same testing procedures similar to those for small ruminants, cattle, and canines, including dilutions. HRP-conjugated goat anti-rabbit IgG (diluted at 1:10,000, Bioworld, MN, United States) was used as the secondary antibody.
2.4 The characterization of BP26 epitopes recognized by mAbs
Six linear polypeptides of BP26 were synthesized, and iELISA was conducted as previously described (19). In brief, 100 μL of peptide-KLH conjugate (10 μg/mL) was coated onto microplates (Corning, NY, United States) in 0.1 M of carbonate buffer (pH 9.6) and incubated overnight at 4°C. The plates were then washed three times with phosphate-buffered saline containing 0.05% Tween-20 (PBST) and blocked with 5% skimmed milk for 2 h at 37°C. After washing three times with PBST, 100 μL of purified mAb at 1 μg/mL was added to each well and incubated at 37°C for 1 h with the same volume of normal mouse serum used as a negative control. The plates were washed again three times and then incubated with 100 μL of horseradish peroxidase (HRP)-conjugated recombinant protein G (diluted at 1:10,000) (Bersee, Beijing, China) for 45 min at 37°C. TMB was used as the colorimetric substrate, and the optical density was measured at 450 nm (OD450) using an ELISA plate reader (BioTek, United States). A ratio of OD450 value of mAb to normal mouse sera of above 1.5 was considered a positive result.
2.5 iELISA and cELISA
The iELISA was conducted following previously established procedures with slight modifications (16). One change was the concentration of BP26 protein used to coat the microplates, which was adjusted to 1 μg/mL. Additionally, commercial blocking buffer (BioFX, United States) was substituted for 5% skimmed milk.
For cELISA, critical conditions, such as the coating concentration of BP26 and the volume of serum and mAbs, were optimized using checkerboard titration. HRP-conjugated rabbit anti-mouse IgG (diluted at 1:10,000) (Sangon, China) served as the secondary antibody to detect the mAb. The remaining steps were identical to those of the iELISA. All samples were repeated three times.
2.6 Statistical analysis
The overall sensitivity and specificity of iELISA and cELISA in detecting animal brucellosis were calculated using receiver operating characteristic (ROC) curves. The optimal cutoff values were defined using the highest sum of sensitivity and specificity. For each optimal cutoff value, the main parameters, such as sensitivity, specificity, positive predictive value (PPV), and negative predictive value (NPV), were calculated. The dot plot was performed using GraphPad Prism version 6.05.
3 Results
3.1 The production of anti-BP26 mAbs and epitope recognition analysis
After screening, a total of 12 hybridoma cell lines were found to secrete mAbs, which were reactive to BP26. Then, mAbs were prepared from mouse ascites and coded as E1–E12. Based on the results of BP26 polypeptides iELISA, the mAb E10 reacted to the polypeptide of QPIYVYPDDKNNLKEPTITGY and three mAbs coded as E4, E5, and E8 reacted to polypeptide AAAPDNSVPIAAGENSYNVSVNVVFE, while the remaining eight mAbs did not react to any of these polypeptides (Table 1).
3.2 The inhibition of animal brucellosis-positive sera to mAbs
cELISA was conducted individually for cattle, small ruminants, and canine brucellosis-positive sera using these mAbs. The inhibition percentage (IP) for each mAb was calculated as [100 − (OD450 value of positive serum/OD450 value of negative serum) × 100] %. Based on the IP values, the positive sera were categorized into five groups: <20%, 20–39%, 40–59%, 60–79%, and ≥80%.
If a serum with an IP value of more than 20% was positive, E10 inhibited 100% of cattle brucellosis sera, 97.62% of small ruminant sera, and 100% of canine sera (Table 2). With the same brucellosis-positive sera, the other mAbs exhibited significantly lower IPs compared to E10. Similarly, when the threshold value of IP was raised to 40%, E10 still displayed the highest IPs for cattle, small ruminants, and canine brucellosis-positive sera, with values of 38.64, 57.14, and 100%, respectively. Based on these results, E10 was selected as the preferred choice for assembling a BP26-based cELISA kit.
3.3 The efficiency of the E10-based cELISA in detecting animal brucellosis
To evaluate the efficacy of mAb E10 in detecting animal brucellosis, cELISA was conducted using positive and negative sera from various domestic animals, including cattle, small ruminants, and canines. The optimal coating concentration of BP26 protein was set at 10 μg/mL, and the optimal serum volume used was 30 μL, which was mixed with 70 μL of pre-diluted mAb (0.875 μg/mL) before transferring to each well of the microplate.
The largest area under the ROC curve was observed for small ruminant sera (AUC = 0.9974), followed by cattle (AUC = 0.9934) and canine sera (AUC = 0.9814) (Figure 1). Using the optimal cutoff values, the cELISA accuracy for small ruminants was the highest (97.94%), with PPV and NPV recorded as 97.62 and 98.18%, respectively (Table 3). Cattle sera showed similar accuracy, although the PPV was slightly lower than that of small ruminant sera. Canine sera exhibited the highest PPV (100%), but due to the limited number of samples, the NPV was the lowest, resulting in the lowest overall accuracy among the three groups of sera.
3.4 The efficiency of the BP26-based iELISA in detecting animal brucellosis
Compared to cELISA, the same collections of brucellosis-positive and brucellosis-negative sera were also tested by the BP26-based iELISA. The optimal coating concentration of BP26 protein for iELISA was 0.3 μg/mL, much lower than cELISA.
According to the result of iELISA, the largest area under the ROC curve was obtained for canine sera (AUC = 1.0000), followed by cattle (AUC = 0.9578) and small ruminant sera (AUC = 0.8126) (Figure 2). Using the optimal cutoff values, the accuracy of cELISA for canine sera was the highest (100%), meaning that all the positive and negative sera were correctly distinguished. Lower accuracy values were obtained for cattle and small ruminant sera (89.58 and 82.30%, respectively) (Table 4).
3.5 Cross-reaction to other serum
If a ratio of S/N (OD450, sample/negative) >2.0 in iELISA or N/S >2.0 in cELISA was considered to be positive, neither E10-based cELISA nor BP26-based iELISA identified the rabbit sera-infected by Y. enterocolitica O9, E. coli (O157:H7, O116), Salmonella urban, Ochrobactrum anthropi, and Vibrio cholerae.
4 Discussion
Among the protein antigens utilized for brucellosis detection, BP26 has been extensively studied. Previous research has shown that the BP26-based iELISA can effectively detect various cases of animal brucellosis. However, our experiment suggests that cELISA using anti-BP26 mAbs may be more efficient. Overall, the cELISA method established in this study outperformed iELISA, particularly in terms of accuracy for detecting cattle and small ruminant sera (Tables 3, 4). Notably, there was a significant difference in small ruminant detection between cELISA (97.62%) and iELISA (64.29%). Conversely, iELISA showed higher accuracy in detecting canine brucellosis, achieving 100% accuracy as compared to 94.34% for cELISA. This difference primarily stemmed from slightly lower specificity in cELISA, although sensitivity remained comparable to iELISA. Unfortunately, the limited availability of brucellosis-positive sera in canine in our experiment reflects the insufficient attention given to canine brucellosis in China, resulting in a shortage of qualified serum samples. Thus, further evaluation is necessary to assess the efficacy of anti-BP26 mAb-based cELISA in detecting rough Brucella-induced brucellosis.
BP26 is an immunogenic protein capable of stimulating various mAbs targeting different epitopes within the mouse model (17, 18). In our experiment, mAb E10 demonstrated significant potential in detecting animal brucellosis, while mAbs E6 and E12 showed some ability but did not meet the necessary sensitivity and specificity requirements for diagnostics. The remaining nine mAbs were largely ineffective as diagnostic reagents. As mAb E10 reacted to the polypeptide QPIYVYPDDKNNLKEPTITGY, its recognized epitope likely lies within this sequence. Interestingly, Qiu et al. (18) reported two epitopes recognized by anti-BP26 mAbs, one of which, with the sequence QPIYVYPD, overlaps with the polypeptide used in our experiment. Although we cannot confirm their exact identity, the antigenicity of the QPIYVYPDDKNNLKEPTITGY sequence suggests that it may be a useful diagnostic antigen for brucellosis. Furthermore, this proposed epitope falls within the BP26 region as described by Seco-Mediavilla et al. (20), between amino acids 55 and 152, which showed superior specificity in detecting small ruminant sera. This finding aligns with our data, indicating that the E10-based cELISA also exhibits high specificity for brucellosis in small ruminants.
Currently, widespread vaccination with attenuated Brucella live vaccines is a key strategy for controlling animal brucellosis in high-prevalence areas, with vaccines possessing DIVA competency being particularly favored by authorities. In China, the commercialization and increasing utilization of the M5ΔBP26 vaccine, which contains a Bp26 deletion, for immunizing small ruminants are notable. The corresponding iELISA utilizing BP26 as an antigen has been proposed to differentiate between M5ΔBP26 vaccinated sera and naturally infected ones (21). Given that the mAb E10-based cELISA established in this study demonstrated better sensitivity and specificity than iELISA, our unpublished data showed that all of the small ruminant sera collected from 7 days to 207 days after vaccination by M5ΔBP26 were not recognized by this cELISA method, demonstrating that the cELISA might be more suitable for the diagnosis of DIVA.
In conclusion, brucellosis remains a significant zoonotic disease, particularly in China, where it is highly endemic in certain regions. The cELISA and iELISA methods outlined in this study not only facilitate the detection of brucellosis in cattle, small ruminants, and canines but also serve as DIVA tests when the Bp26-deleted vaccine is administered. It is hoped that the vaccination and DIVA test strategy will expedite the eradication of animal brucellosis in China.
4.1 Limitations
Although our constructed cELISA effectively distinguishes animal brucellosis sera from non-brucellosis sera and can differentiate between immune sera and naturally infected sera in ruminants, our study has certain limitations. First, the sample size is limited, and expanding it will require verifying the sensitivity and specificity of the constructed method. Second, the M5ΔBP26 vaccine strain is specific to China, limiting the applicability of our method to distinguish vaccine-immune sera, which is used worldwide. Finally, our study only included small ruminants immunized sera, and further investigation is needed to assess its ability to differentiate between immunized and naturally infected sera in other animals.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding authors.
Ethics statement
The animal study was approved by Experimental Animal Ethics Committee of Xuzhou Medical University. The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
XG: Conceptualization, Data curation, Writing – original draft. MS: Conceptualization, Writing – original draft. YG: Conceptualization, Writing – original draft. YW: Writing – review & editing. XY: Writing – review & editing. ML: Writing – review & editing. JL: Writing – review & editing. XS: Writing – review & editing. XF: Writing – review & editing. HZ: Funding acquisition, Writing – review & editing. SS: Conceptualization, Writing – review & editing. JW: Conceptualization, Writing – review & editing. DY: Funding acquisition, Methodology, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This study was supported by the Xuzhou Science and Technology Bureau (Grant number KC23306), the Medical Research Program of Jiangsu Commission of Health (Grant number Z2023080), the National Key Research and Development Program of China (Grant number 2021YFD1800402), and the Postgraduate Research and Practice Innovation Program of Jiangsu Province (Grant number KYCX23-2963).
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.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2024.1389728/full#supplementary-material
Supplementary TABLE 1 Raw data for Figure 1.
Supplementary TABLE 2 Raw data for Figure 2.
Supplementary TABLE 3 Details of the production of BP26 antigen (epitope) and monoclonal antibodies.
References
1. Rossetti, OL, Arese, AI, Boschiroli, ML, and Cravero, SL. Cloning of Brucella abortus gene and characterization of expressed 26-kilodalton periplasmic protein: potential use for diagnosis. J Clin Microbiol. (1996) 34:165–9. doi: 10.1128/jcm.34.1.165-169.1996
2. Cloeckaert, A, Debbarh, HS, Vizcaíno, N, Saman, E, Dubray, G, and Zygmunt, MS. Cloning, nucleotide sequence, and expression of the Brucella melitensis bp26 gene coding for a protein immunogenic in infected sheep. FEMS Microbiol Lett. (1996) 140:139–44. doi: 10.1111/j.1574-6968.1996.tb08327.x
3. Cloeckaert, A, Baucheron, S, Vizcaino, N, and Zygmunt, MS. Use of recombinant BP26 protein in serological diagnosis of Brucella melitensis infection in sheep. Clin Diagn Lab Immunol. (2001) 8:772–5. doi: 10.1128/CDLI.8.4.772-775.2001
4. Zygmunt, MS, Baucheron, S, Vizcaino, N, Bowden, RA, and Cloeckaert, A. Single-step purification and evaluation of recombinant BP26 protein for serological diagnosis of Brucella ovis infection in rams. Vet Microbiol. (2002) 87:213–20. doi: 10.1016/S0378-1135(02)00052-4
5. Kumar, S, Tuteja, U, Kumar, A, and Batra, HV. Expression and purification of the 26 kDa periplasmic protein of Brucella abortus: a reagent for the diagnosis of bovine brucellosis. Biotechnol Appl Biochem. (2008) 49:213–8. doi: 10.1042/BA20070111
6. Tiwari, AK, Kumar, S, Pal, V, Bhardwaj, B, and Rai, GP. Evaluation of the recombinant 10-kilodalton immunodominant region of the BP26 protein of Brucella abortus for specific diagnosis of bovine brucellosis. Clin Vaccine Immunol. (2011) 18:1760–4. doi: 10.1128/CVI.05159-11
7. Tian, M, Song, M, Yin, Y, Lian, Z, Li, Z, Hu, H, et al. Characterization of the main immunogenic proteins in Brucella infectio n for their application in diagnosis of brucellosis. Comp Immunol Microbiol Infect Dis. (2020) 70:101462. doi: 10.1016/j.cimid.2020.101462
8. Nagalingam, M, Basheer, TJ, Balamurugan, V, Shome, R, Kumari, SS, Reddy, GBM, et al. Comparative evaluation of the immunodominant proteins of Brucella abortus for the diagnosis of cattle brucellosis. Vet World. (2021) 14:803–12. doi: 10.14202/vetworld.2021.803-812
9. Chart, H, Okubadejo, OA, and Rowe, B. The serological relationship between Escherichia coli O157 and Yersinia enterocolitica O9 using sera from patients with brucellosis. Epidemiol Infect. (1992) 108:77–85. doi: 10.1017/S0950268800049529
10. Yao, M, Liu, M, Chen, X, Li, J, Li, Y, Wei, YR, et al. Comparison of BP26, Omp25 and Omp31 and a multiepitope-based fusion protein in the serological detection of canine brucellosis. Infect Drug Resist. (2022) 15:5301–8. doi: 10.2147/IDR.S374432
11. Boschiroli, ML, Cravero, SL, Arese, AI, Campos, E, and Rossetti, OL. Protection against infection in mice vaccinated with a Brucella abortus mutant. Infect Immun. (1997) 65:798–800. doi: 10.1128/iai.65.2.798-800.1997
12. Guilloteau, LA, Laroucau, K, Olivier, M, Grillo, MJ, Marin, CM, Verger, J-M, et al. Residual virulence and immunogenicity of CGV26 and CGV2631 B. melitens is Rev. 1 deletion mutant strains in sheep after subcutaneous or conju nctival vaccination. Vaccine. (2006) 24:3461–8. doi: 10.1016/j.vaccine.2006.02.007
13. Jacques, I, Verger, JM, Laroucau, K, Grayon, M, Vizcaino, N, Peix, A, et al. Immunological responses and protective efficacy against Brucella melitensis induced by bp26 and omp31 B. melitensis Rev. 1 deletion mutants in sheep. Vaccine. (2007) 25:794–805. doi: 10.1016/j.vaccine.2006.09.051
14. Grilló, MJ, Marín, CM, Barberán, M, de Miguel, MJ, Laroucau, K, Jacques, I, et al. Efficacy of bp26 and bp26/omp31 B. melitensis Rev. 1 deletion mutants against Brucella ovis in rams. Vaccine. (2009) 27:187–91. doi: 10.1016/j.vaccine.2008.10.065
15. Li, T, Tong, Z, Huang, M, Tang, L, Zhang, H, and Chen, C. Brucella melitensis M5-90Δbp26 as a potential live vaccine that allows for the distinction between natural infection and immunization. Can J Microbiol. (2017) 63:719–29. doi: 10.1139/cjm-2017-0179
16. Bai, Q, Li, H, Wu, X, Shao, J, Sun, M, and Yin, D. Comparative analysis of the main outer membrane proteins of Brucella in the diagnosis of brucellosis. Biochem Biophys Res Commun. (2021) 560:126–31. doi: 10.1016/j.bbrc.2021.04.127
17. Debbarh, HS, Zygmunt, MS, Dubray, G, and Cloeckaert, A. Competitive enzyme-linked immunosorbent assay using monoclonal antibodies to the Brucella melitensis BP26 protein to evaluate antibody responses in infected and B. melitensis Rev. 1 vaccinated sheep. Vet Microbiol. (1996) 53:325–37. doi: 10.1016/S0378-1135(96)01265-5
18. Qiu, J, Wang, W, Wu, J, Zhang, H, Wang, Y, Qiao, J, et al. Characterization of periplasmic protein BP26 epitopes of Brucella melitensis reacting with murine monoclonal and sheep antibodies. PLoS One. (2012) 7:e34246. doi: 10.1371/journal.pone.0034246
19. Yin, D, Bai, Q, Wu, X, Li, H, Shao, J, Sun, M, et al. Correction: paper-based ELISA diagnosis technology for human brucellosis based on a multiepitope fusion protein. PLoS Negl Trop Dis. (2021) 17:e0011079. doi: 10.1371/journal.pntd.0011079
20. Seco-Mediavilla, P, Verger, JM, Grayon, M, Cloeckaert, A, Marín, CM, Zygmunt, MS, et al. Epitope mapping of the Brucella melitensis BP26 immunogenic protein: usefulness for diagnosis of sheep brucellosis. Clin Diagn Lab Immunol. (2003) 10:647–51. doi: 10.1128/CDLI.10.4.647-651.2003
21. Liu, WX, Hu, S, Qiao, ZJ, Chen, WY, Liu, LT, Wang, FK, et al. Expression, purification, and improved antigenic specificity of a truncated recombinant bp26 protein of Brucella melitensis M5-90: a potential antigen for differential serodiagnosis of brucellosis in sheep and goats. Biotechnol Appl Biochem. (2011) 58:32–8. doi: 10.1002/bab.11
Keywords: brucellosis, competitive ELISA, BP26, diagnosis, monoclonal antibodies
Citation: Guo X, Sun M, Guo Y, Wu Y, Yan X, Liu M, Li J, Sun X, Fan X, Zhang H, Sun S, Wang J and Yin D (2024) Production and evaluation of anti-BP26 monoclonal antibodies for the serological detection of animal brucellosis. Front. Vet. Sci. 11:1389728. doi: 10.3389/fvets.2024.1389728
Edited by:
Nitin Vasantrao Kurkure, Maharashtra Animal and Fishery Sciences University, IndiaReviewed by:
Jeanni Fehrsen, Onderstepoort Veterinary Institute (ARC-SA), South AfricaLaxmi Narayan Sarangi, National Dairy Development Board, India
Copyright © 2024 Guo, Sun, Guo, Wu, Yan, Liu, Li, Sun, Fan, Zhang, Sun, Wang and Yin. 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.
*Correspondence: Shufang Sun, c3Vuc2h1ZmFuZ0BjYWhlYy5jbg==; Jianlong Wang, bm1neG1jd2psQDE2My5jb20=; Dehui Yin, eWluZGgxNkB4emhtdS5lZHUuY24=
†These authors have contributed equally to this work