Virulence Gene Characterization and Serotyping of Major Bacterial Pathogens Isolated from Bovine Respiratory Disease in Ethiopia ()
1. Introduction
Ethiopia has the largest number of livestock in Africa and home for various animal species. The country is naturally endowed with different agro-ecological zones and suitable for livestock production. Cattle are the most important livestock subsector in the country and the main providers of nutrition and income for livelihoods. The total cattle population in the country is estimated to be 59.5 million and distributed evenly throughout the country, with higher density in the highland areas [1] [2]. However, the productivity potential remained very low to meet the growing demand of the population. This is due to various constraints such as high prevalence of veterinary disease, low animal health service coverage, and shortage of access to inputs. Moreover, animal health issues are the major barriers to productivity. Prevalence of veterinary important diseases like respiratory illness exert negative impact on productivity and contributes to economic losses in the country [3] [4].
BRD is considered as one of the potential disease in all age of cattle [5] [6]. BRD is caused by diverse bacterial pathogens including Pasteurella multocida (P. multocida), Mannheimia haemolytica (M. haemolytica), Histophilus somni (H. somni), Mycoplasma species, and Trueperella Pyogens (T. pyogens) [7] [8]. Acute infection can be caused by virulent strains of M. haemolytica, P. multocida, and H. somni. These pathogens are opportunistic, possess different virulence genes, and are capable of infecting the lower airway and lung parenchyma [7] [9] [10]. Moreover, M. haemolytica and H. somni do not constantly involve stressors (environmental and management) or other concomitant infections to cause fatal pneumonia [7].
Strains of M. haemolytica can be classified into 17 serotypes using indirect haemagglutination assay of the capsular polysaccharide antigen. 12 serotypes belong to biotype A (1, 2, 5, 6, 7, 8, 9, 12, 13, 14, 16, and 17) and are reclassified as M. haemolytica. 4 serotypes belong to biotype T (3, 4, 10, and 15) and are worth mentioned as B. trehalosi, while M. glucosida was belonged to serotype A: 11 [11] [12]. P. multocida was grouped into five capsular types (A, B, D, E, and F) based on cap loci variation and characterized using primers designed for polymerase chain reaction (PCR) assay [13] [14]. Besides, P. multocida was further classified based on its lipopolysaccharide (LPS) into eight LPS (L1 - L8) genotypes [15].
M. haemolytica and P. multocida possess a range of virulence factors that influence the host defense mechanism and enable colonization and infection [15]. M. haemolytica pathogenicity is associated with various virulence factors of outer membrane proteins, adhesins, leukotoxin, superoxide dismutase, protectins, and hyaluronidase. Virulence genes like transferring-binding protein A (TbpA), leukotoxinC (lktC), ferric-binding protein A (FbpA) which acquire similar properties to the transferrin, and several outer membrane proteins facilitate the development of the disease in the host [16] [17] [18]. Besides, M. haemolytica serotype-specific antigen (ssa) represents virulence-associated gene coding an important virulence marker responsible for the pathogenic potential of M. haemolytica [18] [19] [20]. Likewise, the Rpt2 locus (species-specific) modulates the type III restriction-modification system which acts as a barrier to the introduction of foreign DNA [21].
P. multocida possesses multiple virulence factors associated with fimbriae, capsule, polysaccharide, endotoxins, lipopolysaccharide (LPS), exotoxins, multocidins or siderophores, plasmids, outer membrane proteins (OMPs), and extracellular enzymes. The virulence-associated genes such as toxA, TbpA, sodA, sodC, ompA, ompH, and plpB play a significant role in the pathogenesis of the pathogen [15] [22] [23] [24]. P. multocida surface lipoproteins (pmSLP) are a class of soluble proteins that represent surface and are anchored to the outer membrane [25]. Hence, the identification of these virulence factors indicate the pathogenicity of bacterial pathogens.
Effective control of BRD is likely requires efficient intervention, definitive diagnosis, efficacious vaccines, therapeutic measures, and improved management practices [26]. Hence, BRD call for efficient controlling strategies. Thus, continuous outbreak surveillance and characterization of the major bacterial pathogens associated with BRD are essential for the development of an effective preventive approach. Moreover, molecular epidemiology is essential to determine the distribution of potential pathogenic serotypes circulating in the country. However, there are few studies on the characterization of M. haemolytica and P. multocida circulating in Ethiopia. Besides, information regarding serotypes and virulence genes distribution associated with BRD cases in cattle of Ethiopia is lacking. Therefore, the present study aimed to investigate BRD causing bacterial pathogens, distribution of capsular types, serotypes, and virulence-associated genes of M. haemolytica and P. multocida.
2. Material and Methods
2.1. Sample Collection and Study Area
Samples were collected from BRD outbreak cases, showing respiratory signs including depression, appetite loss, respiratory distress, high fever, nasal discharge, and lacrimation. 166 samples were collected from November 2020 to April 2021. Nasopharyngeal swab samples (n = 107) and a section of lung tissue samples (n = 59) were excised from the pneumonic lungs. Samples were transported to the National Veterinary Institute (NVI), research and development laboratory, Ethiopia. Samples were collected from Hawassa, Ziway, and Mojo. Hawassa is located at 7˚3'1.3464''N and 38˚29'43.8144''E with an elevation of 1712 meters above sea level (m.a.s.l). Ziway is located at 7˚55'39.792''N and 38˚43'15.9888''E at an altitude of 1645 m.a.s.l. Mojo is coordinated between 8˚35'13.2648''N and 39˚7'23.3256''E with an elevation of 1780 m.a.s.l.
2.2. Isolation of Bacterial Pathogens
Bacterial strains were isolated using the standard bacteriological assay. Briefly, pneumonic lung tissue samples were minced, reconstituted in 4 ml sterile physiological saline (pH 7.0 ± 0.2), centrifuged at 3000 × g for 3 min, and the supernatant was discarded. The remaining sediment was reconstituted with 100 µl sterile physiological saline. Ten µl of pneumonic lung sample suspension and nasopharyngeal swabs were streaked comparably onto blood agar base (HiMedia, India) with 5% sheep blood and MacConkey agar (HiMedia, India). Plates were incubated at 37˚C for 24 - 48 hrs aerobically.
Colonies characteristics of Pasteurellaceae were further assayed as per the standard bacteriological method. Presumptive isolates were further analyzed using biochemical test (catalase and oxidase reaction, ornithine decarboxylase (ODC), indole production, nitrate reduction and urease). Identification of BRD-associated bacteria pathogens to species level was conducted using carbohydrate fermentation reaction (glucose, lactose, sucrose, arabinose, trehalose, dulcitol, mannitol, sorbitol, and D-xylose).
2.3. Molecular Assay
A PCR assay was conducted on presumptive bacterial pathogens of BRD as a confirmatory assay. Species-specific universal genes Rpt2 and Kmt1 were used for molecular detection of M. haemolytica and P. multocida, respectively. Genomic DNA was extracted using the Qiagen DNeasy Blood and Tissue Kit (QIAGEN GmbH, Germany) following the manufacturer’s instruction.
2.4. Serotyping
Rabbit antisera were produced against M. haemolytica reference strains which were kindly provided by the NVI, Ethiopia. The assay was conducted using a rapid plate agglutination test against rabbit serum and classified to their respective serotypes [27]. Biotypes of P. multocida were assigned based on their sugar fermentation profiles [28].
2.5. M. haemolytica and P. multocida Virulence-Associated Genes
Virulence-associated genes detection was conducted using the PCR method on twenty-two randomly selected isolates of M. haemolytica (n = 11) and P. multocida (n = 11). PCR was conducted to investigate virulence-associated genes of M. haemolytica (ssa, sodA, FbpA, TbpA, and lktC) and P.multocida (toxA, FbpA, TbpA, and pmSLP) (Table 1 and Table 2).
2.6. Electrophoresis
Amplified PCR products were electrophoresed in 1.5% agarose gels. PCR products (10 µl) were mixed with loading buffer (6×) and loaded into each separate
Table 1. Primer pairs and PCR condition of virulence-associated genes.
F: forward primer; R: reverse primer; bp: base pair.
wells. A 100 bp/1kb plus DNA molecular marker was added into the first and last lane and run at 120 V for 60 min. Products were visualized for the desired size of DNA bands under a gel documentation system (Uvitec, UK).
2.7. Data Analysis
Data were coded and stored in Excel and analyzed using statistical software (STATA 11). Descriptive statistics were used and statistical result was considered at p < 0.05.
Table 2. PCR condition of the different virulent gene primers.
3. Results
3.1. Outbreak Investigation
Diseased animals showed high fever (>40.0˚C), coughing, anorexia, severe respiratory distress, marked depression, salivation, and absence of rales when the ventral lung auscultated. Infected animals appeared dull and respiratory grunts were observed in advanced cases. Cattle slaughtered at the abattoir were inspected for typical gross pathological lesions. The affected parts of the lung showed firm, friable, irregularity in shape, consolidation, and dark red coloration were frequently observed in pneumonic cases. In some advanced cases, pulmonary parenchymal consolidation and interstitial edema were evidenced.
3.2. Biochemical Characterization
Bacterial pathogens isolated from nasopharyngeal and pneumonic lung samples of cattle revealed the identification of M. haemolytica, P. multocida, B. trehalosi, and H. somni. Isolates were Gram-negative, coccobacilli, and pleomorphic (Table 3).
3.3. Isolation of Bacterial Pathogens
Bacteriological assay showed that 37 (22.29%) M. haemolytica isolates and 26 (15.66%) P. multocida isolates were identified. Besides, 12 (7.23%) B. trehalosi isolates and 15 (9.04%) H. somni isolates were phenotypically characterized from outbreak samples at p > 0.05 (Table 4).
Table 3. Biochemical profile of presumptive isolates.
*+: Positive; **-: negative.
Table 4. Prevalence of bacterial pathogens.
*The result is not significant at p < 0.05.
3.4. PCR Assay
PCR assay confirmed that all (100%) M. haemolytica isolates and 25 (96.15%) P. multocida isolates were positive for the Rpt2 and kmt1 universal genes, respectively. Besides, all B. trehalosi isolates were found positive for the sodA gene. Capsular typing of P. multocida isolates posive for kmt1 gene revealed that 21 (84.0%) isolates were belonged to capsular type A (hyaD-hyaC) and 4 (16.0%) were belonged capsular type D (dcbF) (Figure 1).
3.5. Serotyping and Biotying
Serotyping of Rpt2 PCR-positive M. haemolytica isolates revealed that 26 (70.27%) isolates were belonged to serotype A: 1, 4 (10.81%) isolates were serotype A: 2, and 7 (18.92%) were belonged to M. haemolytica serotype A: 6. P. multocida kmt1 and hyaD-hyaC PCR-positive isolates belonged to biotype A: 1, A: 2, and A: 3 from 3 (14.29%), 2 (9.52%), and 16 (76.19%) isolates, respectively.
3.6. Virulence Gene Detection
M. haemolytica strains harbored ssa in 7 (63.63%) and sodA in 10 (90.91%) strains, while FbpA, TbpA, and lktC genes were found in all isolates. P. multocida exhibited 100% toxA, FbpA, TbpA, and pmSLP genes (Table 5 and Figure 2).
Figure 1. Agarose gel electrophoresis (a) (lane M: molecular ladder; lane 1 - 5: M. haemolytica ssa gene (327 bp); lane 7 - 10: P. multocida kmt1 gene (460 bp); lane 11 - 14: P. multocida capsular type D (657 bp); lane N: negative control; lane P1 - P3: positive control); (b) P. multocida capsular type A (1044 bp); (c) B. trehalosi isolates sodA gene (144 bp).
Figure 2. Agarose gel electrophoresis of Virulence-associated genes (a) FbpA gene of M. haemolytica (1000 bp); (b) FbpA gene of P. multocida (500 bp); (c) TbpA gene of P. multocida (728 bp); (d) lktC gene of M. haemolytica (440 bp).
4. Discussion
BRD is a multiplex illness that causes a severe form of pneumonia in all age groups of cattle. It is caused by range of factors such as multiple bacterial and viral pathogens together with environmental stressors are responsible for the outbreak of this disease [29]. BRD infected animals are detected late or not detected
Table 5. Virulence genes detected in M. haemolytica and P. multocida.
at all during the disease development. Hence, early intervention is essential to effectively control the disease and reduce morbidity and mortality. Besides, compiling the prevalence of the major causative pathogens of BRD might help in early intervention and disease control [30] [31]. Therefore, outbreak monitoring, investigation, and identification of bacterial pathogens are significant to design effective controlling strategies.
BRD outbreak investigation aids in defining the frequency of the major infectious pathogens. The clinical signs observed in the current study include coughing, depression, loss of appetite, respiratory distress, and elevated temperature (DART). These findings are consistent with previous studies [7] [8]. Pneumonic lung gross pathological lesion inspection showed irregularity in shape, dark red coloration, marbling, non-friable foci, or fibrinous pleuritis, pulmonary parenchymal consolidation, and interstitial edema in advanced cases were consistent with previous findings [7] [30] [32]. However, diagnosis of BRD infection based on clinical symptoms is difficult. Consequently, diagnosis has to be supported with the identification of the exact pathogens.
In the current outbreak investigation, bacteriological assay revealed higher prevalence of M. haemolytica 37 (22.29%) followed by P. multocida 25 (15.06%). Besides, H. somni and B. trehalosi are identified from 15 (9.04%) and 12 (7.23%) outbreak samples, respectively. This study showed the potential impact of M. haemolytica and P. multocida in the study areas. The current finding indicated a lower prevalence of M. haemolytica as compared to previous reports of 50.59% [3], 46.4% [4] and 29.2% [33]. However, it is higher than other reports of 10.67% and 10.13% of M. haemolytica from the central part of the country [34] [35]. Alike, 14.3% and 12.67% B. trehalosi incidence was reported in previous studies [4] [33]. Therefore, this study highlights the significant role of M. haemolytica, B. trehalosi, P. multocida, and H. somni in causing BRD outbreaks.
PCR assay designed to amplify the rpt2 locus (species-specific) universal gene of M. haemolytica used as a confirmatory assay [18] [21]. Besides, primers targeting fragment of the kmt1 gene of P. multocida produce unique amplification product to all strains [12] [14]. The PCR result of M. haemolytica and P. multocida isolates using the rpt2 and kmt1 universal genes confirmed 37 (100%) and 25 (96.15%) of M. haemolytica and P. multocida isolates, respectively. PCR results of M. haemolytica and P. multocida using the universal primers rpt2 and kmt1 are in agreement with the previous findings in Ethiopia from pneumonic cattle and sheep [3] [12]. Besides, the PCR assay designed to detect P. multocida capsular type A and D strains using specific cap loci (hyaD-hyaC) and (dcbF) genes, respectively [13]. The capsular typing revealed 21 (84.0%) capsular type A (hyaD-hyaC) and 4 (16.0%) capsular type D (dcbF). These findings are in accordance with previous reports [36] [37] and these studies have proved that P. multocida (capsular type A) isolates are the principal causative agent of BRD [38] [39].
Serotyping of PCR positive M. haemolytica (n = 37) showed that 70.27% isolates belonged to serotype A: 1, 10.81% serotype A: 2, and 18.92% belonged to serotype A: 6. Higher prevalence of M. haemolytica serotype A: 1 and distribution of serotype A: 2 and A: 6 in the present study was agreed with other reports who identified these serotypes from cattle suffering from BRD [3] [16] [26]. Besides, biotyping of P. multocida (n = 21) strains indicated that 76.19% isolates were A: 3, while 14.29% and 9.52% were belonged to A: 1 and A: 2, respectively. Higher prevalence of P. multocida A: 3 in the present study was in agreement with previous reports [37] [39].
In the present study, PCR based assay confirmed the distribution of virulence-associated genes. Potential virulence factors of M. haemolytica include ssa, sodA, FbpA, TbpA, and lktC. Among these virulence-associated genes, Lkt is important in inducing pneumonia. The findings revealed that M. haemolytica isolates harbored 100% FbpA, TbpA, and lktC genes while the ssa gene was identified from 63.64% and sodA gene from 90.91% strains. PCR assay results of M. haemolytica are comparable to previous findings, who reported lktC as species-specific in most M. haemolytica isolates [16] [17]. Besides, ssa gene of serotype-specific antigen of outer membrane protein serine protease is considered as the most important virulence factor in M. haemolytica isolates [29].
P. multocida isolates exhibited toxA, FbpA, TbpA, and pmSLP in all isolates. These virulence genes play important roles in the pathogenesis of BRD [15] [23] [26]. The gene encoding dermonecrotoxin (toxA) defines the pathogenic potential of P. multocida strains. This virulence gene is responsible for the clinical symptoms associated with pneumonia [15] [23] [39]. Higher prevalence of TbpA was reported in isolates from diseased cattle [20] [36]. Besides, the surface lipoprotein (SLP) of P. multocida is essential for host colonization [40] and the FbpA gene possess similar properties with transferrins and has a great role in transporting and removing iron across the outer membrane [41]. Hence, among these virulence genes the toxA and tpbA genes of P. multocida are important epidemiological gene markers [26] [42]. Identification of virulence genes implies significant information about the pathogenicity of M. haemolytica and P. multocida isolates.
5. Conclusion
In conclusion, remarkable evidence was identified in the distribution of capsular type, serotypes, and virulence-associated genes of M. haemolytica and P. multocida. PCR assay indicated that M. haemolytica (A: 1) and P. multocida (A: 3) were the most prevalent isolates to cause BRD. Detection of virulence genes of M. haemolytica (ssa, sodA, FbpA, TbpA, and lktC) and P. multocida (toxA, FbpA, TbpA, and pmSLP) in most of the strains implies the pathogenic potential of both pathogens to cause disease outbreak. Hence, the current findings provide relevant information to understand the capsular types, serotypes, and virulence gene incidence associated with M. haemolytica and P. multocida. Therefore, continuous outbreak surveillance from wider areas of the country and molecular epidemiology are indispensable in designing efficient prevention strategies.
Acknowledgements
The authors are grateful to acknowledge the Research and Development Laboratories (NVI, Ethiopia) for allowing the laboratory facility to conduct this research. Great appreciation was forwarded to Veterinary clinics located in the study areas, animal owners, and abattoir staff for their cooperation during sample collection. The authors also thank the Department of Biotechnology, Koneru Lakshmaiah Education Foundation (KLEF) for supporting the study.
Data Availability
All data supporting the findings of this study can be obtained from the corresponding author upon formal request.
Ethical Consideration
Consent was first obtained from animal owners. Samples collection followed scientific protocols and animal handling employed with basic animal welfare ethics. Laboratory assay was executed following the standard bacteriological and molecular methods.
List of Abbreviations
BRD: Bovine Respiratory Disease; FbpA: Ferric-binding protein A; LPS: Lipopolysaccharide; lktC: leukotoxinC; pmSLP: P. multocida surface lipoproteins; PCR: Polymerase Chain Reaction; TbpA: Transferring-binding protein A; ssa: M. haemolytica serotype-specific antigen.