Next Article in Journal
Standardization of Data Analysis for RT-QuIC-Based Detection of Chronic Wasting Disease
Previous Article in Journal
Neoehrlichia mikurensis—A New Emerging Tick-Borne Pathogen in North-Eastern Poland?
Previous Article in Special Issue
Serotyping and Antimicrobial Susceptibility Profiling of Glaesserella parasuis Isolated from Diseased Swine in Brazil
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Molecular Characterization and Phylogenetic Analysis of Outer membrane protein P2 (OmpP2) of Glaesserella (Haemophilus) parasuis Isolates in Central State of Peninsular Malaysia

1
Department of Veterinary Clinical Studies, Faculty of Veterinary Medicine, Universiti Putra Malaysia, UPM, Serdang 43400, Malaysia
2
Department of Bioprocess Technology, Faculty of Biotechnology and Biomolecular Sciences, Universiti Putra Malaysia, UPM, Serdang 43400, Malaysia
3
Department of Veterinary Medicine, College of Veterinary Medicine, National Chiayi University, Chiayi City 60004, Taiwan
*
Authors to whom correspondence should be addressed.
Pathogens 2023, 12(2), 308; https://doi.org/10.3390/pathogens12020308
Submission received: 30 December 2022 / Revised: 6 February 2023 / Accepted: 9 February 2023 / Published: 12 February 2023
(This article belongs to the Special Issue About Glaesserella (Haemophilus) Parasuis and Glässer’s Disease)

Abstract

:
Glaesserella (Haemophilus) parasuis, the etiological agent of Glässer’s disease, is an economically significant pathogen commonly associated with serofibrinous polyserositis, arthritis, fibrinous bronchopneumonia and/or meningitis. This study is the first attempt to molecularly characterize and provide a detailed overview of the genetic variants of G. parasuis present in Malaysia, in reference to its serotype, virulence-associated trimeric autotransporters (vtaA) gene and outer membrane protein P2 (OmpP2) gene. The G. parasuis isolates (n = 11) from clinically sick field samples collected from two major pig producing states (Selangor and Perak) were selected for analysis. Upon multiplex PCR, the majority of the isolates (eight out of 11) were identified to be serotype 5 or 12, and interestingly, serotypes 3, 8 and 15 were also detected, which had never been reported in Malaysia prior to this. Generally, virulent vtaA was detected for all isolates, except for one, which displayed a nonvirulent vtaA. A phylogenetic analysis of the OmpP2 gene revealed that the majority of Malaysian isolates were clustered into genotype 1, which could be further divided into Ia and Ib, while only one isolate was clustered into genotype 2.

1. Introduction

Glässer’s disease is caused by the previously known Haemophilus parasuis. It has recently been reclassified to Glaesserella parasuis after a series of detailed phylogenomic analysis, in honor of K. Glässer, who first described the disease in pigs, caused by a fastidious, Gram-negative coccobacillus isolated from the respiratory tract or systemic sites of swine [1]. G. parasuis may be present in most commercial pig farms as a colonizer of the upper respiratory tract of healthy pigs. It can be detected in the nasal cavity as early as 2 days after birth in nonvaccinated herds, as part of their normal respiratory microbiota [2]. On the other hand, opportunistic G. parasuis strains can lead to the onset of lower respiratory tract disease (involved in porcine respiratory disease complex) after weaning, or other clinical presentations characterized by serofibrinous polyserositis, arthritis, and/or meningitis. This often sequels with high morbidity and mortality, especially in the naive swine population, causing significant economic impact to the pig production industry [3].
To outline effective surveillance, control and preventive measures, various efforts in the characterization and identification of G. parasuis virulence indicators have been described over the years, gaining a broader understanding of the disease’s nature. Typically, G. parasuis is characterized by its serotype (1 to 15) [4,5]. The correlation between serotype and virulence is doubtful; there are conflicting reports on different strains of the same serotype showing variable degrees of symptoms, or different pigs challenged by the same strain showing variable symptoms [6,7,8,9]. Having said that, molecular typing or sequencing methods present a major advancement for epidemiological and vaccine studies, where virulence genes, or even single gene point mutations, can be identified [6,10]. For example, molecular pathotype (i.e., the leader sequence (LS)-PCR classification of virulence and nonvirulence based on vta-locus) [11,12], and also by other molecular genotyping methods [6,8,9]. Furthermore, the molecular pathotyping method can be performed easily and correspond well to clinical background information [12].
The virulence-associated trimeric autotransporter (vtaA) is one important virulent gene involved in adherence to the host cell and extracellular protein during an infection [13,14]. The vtaA genes of G. parasuis are differentiated into three groups based on the YadA domains; vtaA1 and vtaA2 are associated with virulence and are a potential indicator of virulent isolates, while vtaA3 is highly conserved and is present in all strains [15]. A leader sequence (LS)-PCR based on the extended signal peptide region (ESPR) of the vtaA, reported by Galofré-Milà et al. (2017), distinguishes nonvirulent G. parasuis commensals from the potentially disease-causing virulent isolates, enabling the prediction of virulence potential of G. parasuis [11].
Aside from the trimeric autotransporters (vtaA) antigen, the outer membrane protein P2 (OmpP2) is amongst the protective antigen identified for G. parasuis subunit vaccine developments [16,17,18]. The OmpP2 gene is the most abundant protein amongst the outer membrane protein, highly conserved and abundant in Gram-negative bacteria [19,20]. It could be divided into two genotypes, which are genetic type-I and type-II [21]. The length of the OmpP2 gene of genetic type-I ranges from 1077 bp to 1095 bp, while the length of the OmpP2 gene of genetic type-II ranges from 1167 bp to 1203 bp [22]. The OmpP2 gene is believed to be associated with the virulence of G. parasuis, as it is able to induce proinflammatory cytokine mRNA expression in porcine alveolar macrophages (PAM) [23,24]. The OmpP2 is considered an immunodominant porin, where it can be targeted as a potential vaccine candidate, subjected to immunogenic characterization, pathogenesis studies, and also serve as a diagnostic antigen [25,26,27,28].
In this study, to characterize the G. parasuis isolates of diseased pigs from major pig producing states (Selangor and Perak) in Malaysia, we applied molecular PCR-based and sequencing methods, which were: molecular serotyping using serotype-specific primers, as described by Howell et al. (2015); molecular pathotyping based on the leader sequence (LS)-PCR of vtaA gene; as well as molecular genotyping via the phylogenetic analysis of the outer membrane protein P2 (OmpP2) gene.

2. Materials and Methods

2.1. Bacteria Isolates

The central states of Perak and Selangor contributes the highest pig production in Peninsular Malaysia. Farm sizes from both states have sow numbers from 100 to 4000. Due to the pandemic and the threat of African Swine Fever, access to farms, or farms submitting samples specifically for bacteria isolation, was limited. There were 28 Glaesserella parasuis isolates retrievable from archive samples (2018–2020; n = 24) and recent clinical cases (2021–2022; n = 4) presented to the swine unit of the Faculty of Veterinary Medicine, Universiti Putra, Malaysia. The isolation of G. parasuis from clinical samples was only performed upon special request by farmers and was not a common routine procedure. The isolation of G. parasuis was carried out by inoculating samples submitted on plated chocolate agar (Oxoid, USA), identified by conventional PCR, and positive isolates were preserved at 10% skim milk kept in −80 °C.
The skim milk preserved isolates were subcultured by inoculating 5 µL of skim milk containing G. parasuis on plated chocolate agar (Thermo Fisher Scientific, Inc., Waltham, MA, USA) and incubated in a candle jar (which provided a condition of approximate 5% CO2) at 37 °C for 24 to 48 h. A loopful of bacteria colonies from the pure culture was suspended into 100 µL deionized water reaching the turbidity of 0.5 McFarland standard. DNA was extracted from the bacteria suspension using a column-based extraction kit (DNeasy® Blood and Tissue Kit, Qiagen, Hilden, North Rhine-Westphalia, Germany) performed according to the manufacturer’s protocol specified for Gram-negative bacteria. The extracted DNA was used for the subsequent PCR assays. To confirm that the bacteria colonies were indeed G. parasuis, a conventional PCR assay was carried out using a published species-specific primer of G. parasuis [4]. Briefly, a 20 µL PCR reaction mixture of 12.5 µL 2× Mytaq™ Red Mix PCR master mix (Meridian Bioscience®, Cincinnati, Ohio, USA), 0.5 nM of forward and reverse primer each, 2 µL (50–100 ng/µL) template and 3.5 µL nuclease-free water was prepared for each sample. The PCR assay was performed at 95 °C for 1 min, followed by 30 cycles of 95 °C for 15 s, 57 °C for 30 s and 72 °C for 30 s, and a final extension of 72 °C for 5 min. All PCR products were subjected to 2% agarose gel electrophoresis in TAE buffer 80 V for 45 min, using a 100-bp molecular weight marker (Qiagen, Hilden, North Rhine-Westphalia, Germany) as a guide. Electrophoresed gels were visualized using a UV transilluminator and gel documentation system (Syngene, Frederick, MD, USA). A PCR band at 275 bp denoted sample was positive of G. parasuis.

2.2. Multiplex PCR Assay for Serotyping

The isolates were further serotyped by the molecular method published by Howell et al., 2015, following modifications by Schuwerk et al., 2020 [12], which divided the assay into two sets of multiplex PCR (mPCR). The PCR reactions were set up consisting of 12.5 µL 2× Mytaq™ HS Red Mix PCR master mix (Meridian Bioscience®, Cincinnati, OH, USA), 0.5 nM of forward and reverse primer each, 2 µL (50–100 ng/µL) template and nuclease-free water topped-up to a volume of 25 µL. The PCR assay was performed at 95 °C for 1 min, followed by 30 cycles of 95 °C for 15 s, 58 °C for 30 s and 72 °C for 30 s, and a final extension of 72 °C for 5 min. All PCR products were subjected to 2% agarose gel electrophoresis in TAE buffer 80 V for 45 min, using a 100-bp molecular weight marker (Qiagen, Hilden, North Rhine-Westphalia, Germany) as a guide. Electrophoresed gels were visualized using a UV transilluminator and gel documentation system (Syngene, Frederick, MD, USA). The PCR bands were measured using the Image Lab software v6.1 (Bio-rad Laboratories, Inc., Hercules, CA, USA).

2.3. Conventional PCR Assay for vtaA Gene Identification

One forward and two reverse primers, outlined by Galofré-Milà et al., 2017 [11], were used to distinguish the vtaA gene of G. parasuis into virulent and nonvirulent strains. The PCR reactions were set up comprising of 12.5 µL 2× Mytaq™ HS Red Mix PCR master mix (Meridian Bioscience®, Cincinnati, OH, USA), 0.5 nM of forward and reverse primer each, 2 µL (50–100 ng/µL) template and nuclease-free water topped-up to a volume of 25 µL. The PCR assay was performed at 95 °C for 1 min, followed by 30 cycles of 95 °C for 15 s, 54 °C for 30 s, 72 °C for 30 s and a final extension of 72 °C for 5 min. The PCR products were subjected to 2% agarose gel electrophoresis in TAE buffer 80 V for 45 min, using a 100-bp molecular weight marker (Qiagen, Hilden, North Rhine-Westphalia, Germany) as a guide. Electrophoresed gels were visualized using a UV transilluminator and gel documentation system (Syngene, Frederick, MD, USA). The PCR bands were measured using the Image Lab software v6.1 (Bio-rad Laboratories, Inc., Hercules, CA, USA).

2.4. The PCR Amplification and Bioinformatic Analysis of OmpP2 Gene

The extracted DNA were subjected to conventional PCR assay targeting the G. parasuis OmpP2 gene using primers described by Li et al., 2012 [27]. Briefly, the PCR reactions were set up consisting of 12.5 µL 2× Mytaq™ Red Mix PCR master mix (Meridian Bioscience®, Cincinnati, OH, USA), 0.5 nM of primer each, 2 µL template and nuclease-free water topped-up to a volume of 25 µL. The PCR assay was performed at 95 °C for 1 min, followed by 30 cycles of 95 °C for 15 s, 54 °C for 30 s, an extension at 72 °C for 30 s and a final extension of 72 °C for 5 min. The PCR products were subjected to 2% agarose gel electrophoresis in TAE buffer 80 V for 45 min, using a 100-bp molecular weight marker (Qiagen, Hilden, North Rhine-Westphalia, Germany) as a guide. Electrophoresed gels were visualized using a UV transilluminator and gel documentation system (Syngene, Hercules, CA, USA).
The PCR product displaying a single band at approximate 1100 bp was purified using MEGAquick-spinTM Plus total fragment DNA purification kit (iNtRON Biotechnology, Seongnam-si, Gyeonggi-do, South Korea) and outsourced to a sequencing company (Macrogen Asia Pacific Pte. Ltd., Singapore). Samples were sequenced for both forward and reverse strands via Sanger sequencing. Only good quality sequences with quality scores ≥40 were analyzed. Consensus sequences were obtained by aligning the forward sequence and the reverse complement of the reverse sequences using the Mega 11 software [29]. The sequences were analyzed using bioinformatic tools, such as the Mega 11 software for phylogenetic tree construction, and the Clustal Omega (EMBL-EBI) web-based tool for pairwise similarity comparison [29,30]. The sequences were deposited into the NCBI GenBank; their accession numbers are listed in Table 1.

3. Results

3.1. Bacteria Isolates

Among the 28 retrievable samples, 11 G. parasuis isolates were selected from seven larger producing farms with sow number ranges from 300 to 4000. One to three isolates were selected from eight pigs with different isolation sites, namely, lung, pleural swab, brain, peritoneal swab, joint and tonsil.

3.2. The Molecular Serotyping of G. Parasuis

The G. parasuis isolates were serotyped via mPCR with serotype-specific primers, as shown in Figure 1a,b. It is noteworthy that the mPCR used can distinguish between 14 out of 15 previously described serovars, except serovar 5 and serovar 12, which were detected by the same pair of primers [4].
Serotype 5 or 12 (73%) was found to be the most prevalent serotype, followed by serotypes 3 (9%), 8 (9%) and 15 (9%). These were isolated from cases with polyserositis (73%), arthritis (45%) and fibrinous pneumoniae (18%). There were also multiple G. parasuis serotypes isolated from one infected pig showing clinical signs of polyserositis from Selangor, namely, serotype 5 or 12 (GP/UPM/MY005) and serotype 15 (GP/UPM/MY006). Interestingly, serotype 15 was isolated from the tonsils, whereas serotype 5 or 12 was isolated from brain samples in that single animal.

3.3. Virulence-Associated Trimeric Autotransporters (vtaA) Gene

The gel electrophoresis of PCR assay revealed that only one isolate (GP/UPM/MY006) displayed a nonvirulent vtaA band (208 bp) as shown in Figure 2a and 10 out of 11 isolates displayed virulent vtaA bands (179 bp) as shown in Figure 2b. The isolates that displayed virulent vtaA were all sampled from clinical sites, e.g., lung, brain and joint, while the only nonvirulent vtaA isolate was isolated from a carrier site, which was the tonsil of a clinically ill pig.

3.4. Genotyping via OmpP2 Gene

3.4.1. Amplification of OmpP2 Gene

The primers used in this study worked well at 54 °C temperature to amplify the OmpP2 gene and produced a target band estimated of 1100 bp (Figure 3) for all 11 isolates, based on PCR optimization using gradient annealing temperature. The 1100 bp region of the OmpP2 gene amplified by the primers, as evidenced by the NCBI BLAST search results of the sequenced PCR product.

3.4.2. The Bioinformatic Analysis of G. parasuis OmpP2 Gene Sequence

The sequences of the OmpP2 gene of G. parasuis isolates of Malaysia vary in length from 1077 to 1191 nucleotides. Pairwise comparison using Clustal Omega revealed Malaysia isolates shared a 94.25% to 100% homologous identity among each other as listed in Table S1. In reference to the other 81 G. parasuis strains, Malaysian isolates shared a high homology of 92.98% to 99.87%, except for two strains, Hs-DY06 and Hs-DY13 of China (59.65–65.19%), which is unexpectedly different to the other 79 reference strains. In comparison to Actinobacillus pleuropneumonia and Haemophilus influenzae, the OmpP2 gene was indeed distinct to G. parasuis with only 43.83–55.94% homology.
In exclusion to the above-mentioned two distinct G. parasuis reference strains (Hs-DY06 and Hs-DY13), a total of 79 reference sequences of the G. parasuis OmpP2 gene were selected for the construction of a phylogenetic tree (Figure 4). The sequences were divided into two clusters, which corresponded with the grouping of genetic type-I (Cluster I) and genetic type-II (Cluster 2) from previous studies [21,31]. Within Cluster 1, it could further be divided into two groups, which we tentatively referred to as genetic type-Ia (Cluster Ia) and genetic type-Ib (Cluster Ib). The Malaysian isolates of GP/UPM/MY001, GP/UPM/MY002, GP/UPM/MY003 and GP/UPM/MY004 were closely related to the USA strain 685-99, the Germany strain 84-15995, China strains Hs-DY05 and F603, and clustered among genetic type Ia group. On the other hand, isolates GP/UPM/MY007, GP/UPM/MY005, GP/UPM/MY008, GP/UPM/MY009, GP/UPM/MY010 and GP/UPM/MY011 were more closely related to China strain SW124 and USA strain MN-H; clustered among genetic type Ib. There was only one isolate that was clustered among genetic type-II, that was isolate GP/UPM/MY006, which was most closely related to China strain LHDR_HPS_1_2.
The further comparison of subclusters (Cluster Ia and Cluster Ib) and clusters (Cluster I and II) revealed a homology of 95.2–98.9% and 92.98–97.32%. All sequences in Cluster II generally have a longer nucleotide length (1167–1203) as compared to Cluster I (1077–1095). Information on each isolate, and their results, were summarized in Table 1.

4. Discussion

Although serotyping provides only part of the information on the virulence of G. parasuis, it is important in terms of selecting the best vaccine for control. This is because current vaccine employs the inactivation/attenuation of bacteria whole cell from one or two serotypes, which might not produce desirable protection to heterologous challenge. Currently available commercial vaccines targets serotypes 4, 5 and 12; 1 and 6, or serotype 5 only [35], while the only approved commercial vaccine in Malaysia targets serotype 3, 4 and 5 [36]. When compared to the results of this study, vaccines which may be suitable would most likely be those against the highly virulent serotype 5 or 12 strains, and to a lesser extent those against serotype 8 and 15. Nonetheless, the majority of the isolates described in this study, which were obtained from two major pig producing states in Malaysia (Perak and Selangor), were detected to be serotype 5 or 12. This is consistent to the predominance of serotype 5 or 12 found in North America, Europe and Asia [37]. For a better overview of the G. parasuis serotypic distribution within the country, the current study may be expanded to include all pig producing states in Malaysia.
The upregulation of vtaA proteins during infection supports the role of the protein in infection, especially by virulent strains [14]. This is coherent to findings where virulent vtaA was detected in all strains isolated from clinical sites (lung, pleural wall joint, brain and peritoneal wall), and nonvirulent vtaA was detected in one strain that was isolated from the carrier site (tonsil). This demonstrated that the difference in vtaA might be due to isolation sites, whether G. parasuis was isolated from carrier or systemic sites, which was similar to findings in this study [12].
After the Nipah virus outbreak in Malaysia (1998–1999), there is a lack of national breeder units within the country [38]. Since then, Malaysia has been importing live pigs, for breeding purposes, from Canada, Denmark, Finland, France, South Korea, Sweden, the United Kingdom and the USA [39]. However, the team could not find an agreement to the point where the findings of closely related strains were due to direct transmission from imported breeder animals to the farm, as Malaysian isolates were found to be closely related to Germany and China. This is highly possible due to the lack of the OmpP2 gene sequence information deposited by the country that Malaysia is importing from, as well as from neighboring countries. Nevertheless, as only one gene was used for comparison, the information gained was limited.
In this study, strains in Cluster I (genetic type-I) had a shorter OmpP2 gene nucleotide length, which mainly displayed virulent vtaA, and the opposite was seen in strains of Cluster II. The shorter length of OmpP2 proteins was demonstrated to exhibit significantly increased resistance to complement killing, complementing the idea that genetic type I strains are likely to be causing disease [40]. The detection of virulent vtaA in Malaysia isolates belonging to OmpP2 genetic type I, and nonvirulent vtaA in isolates belonging to OmpP2 genetic type II, seems to be parallel. However, mapping information from previous literature [41] showed that strains from genetic type II, such as strain H465 and SW124 of Germany, also possess virulent vtaA, even though all strains in the genetic type I were found to possess virulent vtaA. This may indicate that the two genes are not directly related, but when combined, both may yield more severe clinical infections; this is worth being studied further.

5. Conclusions

In conclusion, this study reveals that more than 91% of the isolates isolated from pigs suffering from Glässer’s disease belong to the OmpP2 genotype 1 cluster, and that a substantial proportion of disease (91%) was caused by vtaA virulent strains.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/pathogens12020308/s1, Table S1: Percent identity matrix of pairwise comparison of 11 Malaysian G. parasuis OmpP2 gene nucleotide sequences, with 83 reference sequences obtained from the NCBI GenBank. This comprises reference gene sequences of 81 G. parasuis OmpP2 gene of different countries, one Actinobacillus pleuropneumoniae OmpP2 gene and one Haemophilus influenzae OmpP2 gene. Sequences were labeled in the format of NCBI GenBank Accession number_Strain ID_Country_Location in country_ serotype; depending on which information was available. The Malaysian gene sequences ID were highlighted in blue. Values in the figure were colored to scale, with the highest value in green, median value in yellow and lowest value in red.

Author Contributions

Conceptualization, C.Y.L., S.E.L., L.Y.P., J.L., P.T.O. and M.W.C.F.; methodology, C.Y.L. and H.X.O.; software, C.Y.L., C.Y.T. and H.X.O.; validation, C.Y.L., P.T.O. and M.W.C.F.; formal analysis, C.Y.L., H.X.O., P.T.O. and M.W.C.F.; writing—original draft preparation, C.Y.L. and M.W.C.F.; writing—review and editing, C.Y.L., H.X.O., P.T.O. and M.W.C.F.; visualization, C.Y.L.; supervision, M.W.C.F., J.L. and P.T.O.; project administration, C.Y.L. and S.E.L.; funding acquisition, P.T.O. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

This study was granted approval from the Institutional Animal Care and Use Committee (IACUC) under AUP Code UPM/IACUC/AUP- AUP-R064/2020 and was conducted adhering to the guidelines as stated in the Code of Practice for Care and Use of Animals for Scientific Purposes, as stipulated by Universiti Putra Malaysia.

Informed Consent Statement

Not applicable.

Data Availability Statement

The Malaysian Glaesserella parasuis OmpP2 sequences reported in this study have been deposited at GenBank (http://www.ncbi.nlm.nih.gov) under accession numbers OP651005–OP651015.

Acknowledgments

The research team is thankful to the Veterinary Clinical Laboratory and Virology Laboratory, Faculty of Veterinary Medicine, Universiti Putra Malaysia, for assistance in this research.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dickerman, A.; Bandara, A.B.; Inzana, T.J. Phylogenomic analysis of Haemophilus parasuis and proposed reclassification to Glaesserella parasuis, gen. nov., comb. nov. Int. J. Syst. Evol. Microbiol. 2020, 70, 180–186. [Google Scholar] [CrossRef] [PubMed]
  2. Cerdà-Cuéllar, M.; Naranjo, J.F.; Verge, A.; Nofrarías, M.; Cortey, M.; Olvera, A.; Segalés, J.; Aragon, V. Sow vaccination modulates the colonization of piglets by Haemophilus parasuis. Vet. Microbiol. 2010, 145, 315–320. [Google Scholar] [CrossRef]
  3. Zhang, B.; Tang, C.; Liao, M.; Yue, H. Update on the pathogenesis of Haemophilus parasuis infection and virulence factors. Vet. Microbiol. 2014, 168, 1–7. [Google Scholar] [CrossRef]
  4. Howell, K.J.; Peters, S.E.; Wang, J.; Hernandez-Garcia, J.; Weinert, L.A.; Luan, S.L.; Chaudhuri, R.R.; Angen, Ø.; Aragon, V.; Williamson, S.M.; et al. Development of a multiplex PCR assay for rapid molecular serotyping of Haemophilus parasuis. J. Clin. Microbiol. 2015, 53, 3812–3821. [Google Scholar] [CrossRef] [PubMed]
  5. Kielstein, P.; Rapp-Gabrielson, V.J. Designation of 15 serovars of Haemophilus parasuis on the basis of immunodiffusion using heat-stable antigen extracts. J. Clin. Microbiol. 1992, 30, 862–865. [Google Scholar] [CrossRef]
  6. Olvera, A.; Segalés, J.; Aragón, V. Update on the diagnosis of Haemophilus parasuis infection in pigs and novel genotyping methods. Vet. J. 2007, 174, 522–529. [Google Scholar] [CrossRef]
  7. Aragon, V.; Cerdà-Cuéllar, M.; Fraile, L.; Mombarg, M.; Nofrarías, M.; Olvera, A.; Sibila, M.; Solanes, D.; Segalés, J. Correlation between clinico-pathological outcome and typing of Haemophilus parasuis field strains. Vet. Microbiol. 2010, 142, 387–393. [Google Scholar] [CrossRef]
  8. Brockmeier, S.L.; Loving, C.L.; Mullins, M.A.; Register, K.B.; Nicholson, T.L.; Wiseman, B.S.; Baker, R.B.; Kehrli, M.E., Jr. Virulence, transmission, and heterologous protection of four isolates of Haemophilus parasuis. Clin. Vaccine. Immunol. 2013, 20, 1466–1472. [Google Scholar] [CrossRef]
  9. Howell, K.J.; Weinert, L.A.; Chaudhuri, R.R.; Luan, S.L.; Peters, S.E.; Corander, J.; Harris, D.; Angen, Ø.; Aragon, V.; Bensaid, A.; et al. The use of genome wide association methods to investigate pathogenicity, population structure and serovar in Haemophilus parasuis. BMC Genom. 2014, 15, 1179. [Google Scholar] [CrossRef] [PubMed]
  10. Ruiz, A.; Oliveira, S.; Torremorell, M.; Pijoan, C. Outer membrane proteins and DNA profiles in strains of Haemophilus parasuis recovered from systemic and respiratory sites. J. Clin. Microbiol. 2001, 39, 1757–1762. [Google Scholar] [CrossRef] [Green Version]
  11. Galofré-Milà, N.; Correa-Fiz, F.; Lacouture, S.; Gottschalk, M.; Strutzberg-Minder, K.; Bensaid, A.; Pina-Pedrero, S.; Aragon, V. A robust PCR for the differentiation of potential virulent strains of Haemophilus parasuis. BMC Vet. Res. 2017, 13, 124. [Google Scholar] [CrossRef]
  12. Schuwerk, L.; Hoeltig, D.; Waldmann, K.H.; Strutzberg-Minder, K.; Valentin-Weigand, P.; Rohde, J. Serotyping and pathotyping of Glaesserella parasuis isolated 2012–2019 in Germany comparing different PCR-based methods. Vet. Res. 2020, 51, 137. [Google Scholar] [CrossRef]
  13. Cotter, S.E.; Surana, N.K.; Geme, J.W.S., III. Trimeric autotransporters: A distinct subfamily of autotransporter proteins. Trends Microbiol. 2005, 13, 199–205. [Google Scholar] [CrossRef]
  14. Costa-Hurtado, M.; Barba-Vidal, E.; Maldonado, J.; Aragon, V. Update on Glässer’s disease: How to control the disease under restrictive use of antimicrobials. Vet. Microbiol. 2020, 242, 108595. [Google Scholar] [CrossRef] [PubMed]
  15. Pina, S.; Olvera, A.; Barceló, A.; Bensaid, A. Trimeric autotransporters of Haemophilus parasuis: Generation of an extensive passenger domain repertoire specific for pathogenic strains. J. Bacteriol. 2009, 191, 576–587. [Google Scholar] [CrossRef] [PubMed]
  16. Olvera, A.; Pina, S.; Pérez-Simó, M.; Aragón, V.; Segalés, J.; Bensaid, A. Immunogenicity and protection against Haemophilus parasuis infection after vaccination with recombinant virulence associated trimeric autotransporters (VtaA). Vaccine 2011, 29, 2797–2802. [Google Scholar] [CrossRef] [PubMed]
  17. Confer, A.W.; Ayalew, S. The OmpA family of proteins: Roles in bacterial pathogenesis and immunity. Vet. Microbiol. 2013, 163, 207–222. [Google Scholar] [CrossRef]
  18. Liu, H.; Xue, Q.; Zeng, Q.; Zhao, Z. Haemophilus parasuis vaccines. Vet. Immunol. Immunopathol. 2016, 180, 53–58. [Google Scholar] [CrossRef]
  19. Mullins, M.A.; Register, K.B.; Bayles, D.O.; Loving, C.L.; Nicholson, T.L.; Brockmeier, S.L.; Dyer, D.W.; Phillips, G.J. Characterization and comparative analysis of the genes encoding Haemophilus parasuis outer membrane proteins P2 and P5. J. Bacteriol. 2009, 191, 5988–6002. [Google Scholar] [CrossRef]
  20. Zhou, M.; Guo, Y.; Zhao, J.; Hu, Q.; Hu, Y.; Zhang, A.; Chen, H.; Jin, M. Identification and characterization of novel immunogenic outer membrane proteins of Haemophilus parasuis serovar 5. Vaccine 2009, 27, 5271–5277. [Google Scholar] [CrossRef]
  21. Li, P.; Bai, J.; Li, J.X.; Zhang, G.L.; Song, Y.H.; Li, Y.F.; Wang, X.W.; Jiang, P. Molecular cloning, sequencing, and expression of the outer membrane protein P2 gene of Haemophilus parasuis. Res. Vet. Sci. 2012, 93, 736–742. [Google Scholar] [CrossRef]
  22. Li, L.; Yang, Z.; Wang, Y.; Jin, Q.; Wu, X.; Zhang, R. The omp2 gene of HPS-type bacteria cloning and sequence analysis isolates from Sichuan Province. Engineering 2012, 5, 153–158. [Google Scholar] [CrossRef]
  23. Zhang, B.; Xu, C.; Liao, M. Outer membrane protein P2 of the Haemophilus parasuis SC096 strain contributes to adherence to porcine alveolar macrophages cells. Vet. Microbiol. 2012, 158, 226–227. [Google Scholar] [CrossRef] [PubMed]
  24. Zhou, S.; He, X.; Xu, C.; Zhang, B.; Feng, S.; Zou, Y.; Li, J.; Liao, M. The outer membrane protein P2 (OmpP2) of Haemophilus parasuis induces proinflammatory cytokine mRNA expression in porcine alveolar macrophages. Vet. J. 2014, 199, 461–464. [Google Scholar] [CrossRef]
  25. Davies, R.L.; Lee, I. Sequence diversity and molecular evolution of the heat-modifiable outer membrane protein gene (ompA) of Mannheimia (Pasteurella) haemolytica, Mannheimia glucosida, and Pasteurella trehalosi. J. Bacteriol. 2004, 186, 5741–5752. [Google Scholar] [CrossRef]
  26. Hiltke, T.J.; Schiffmacher, A.T.; Dagonese, A.J.; Sethi, S.; Murphy, T.F. Horizontal transfer of the gene encoding outer membrane protein P2 of nontypeable Haemophilus influenzae, in a patient with chronic obstructive pulmonary disease. J. Infect. Dis. 2003, 188, 114–117. [Google Scholar] [CrossRef] [PubMed]
  27. Oliveira, S.; Pijoan, C. Haemophilus parasuis: New trends on diagnosis, epidemiology and control. Vet. Microbiol. 2004, 99, 1–12. [Google Scholar] [CrossRef]
  28. Sikkema, D.J.; Murphy, T.F. Molecular analysis of the P2 porin protein of nontypeable Haemophilus influenzae. Infect. Immun. 1992, 60, 5204–5211. [Google Scholar] [CrossRef]
  29. Tamura, K.; Stecher, G.; Kumar, S. MEGA11: Molecular evolutionary genetics analysis version 11. Mol. Biol. Evol. 2021, 38, 3022–3027. [Google Scholar] [CrossRef] [PubMed]
  30. Sievers, F.; Wilm, A.; Dineen, D.; Gibson, T.J.; Karplus, K.; Li, W.; Lopez, R.; McWilliam, H.; Remmert, M.; Söding, J.; et al. Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Mol. Syst. Biol. 2011, 7, 539. [Google Scholar] [CrossRef]
  31. Wu, X.; Xiao, L.; Wang, Y.; Yao, X.; Yang, Z. Genetic variants and phylogenetic analysis of Haemophilus parasuis (HPS) OMPP2 detected in Sichuan, China from 2013 to 2015. J. Vet. Med. Sci. 2017, 79, 1648–1651. [Google Scholar] [CrossRef] [PubMed]
  32. Saitou, N.; Nei, M. The neighbor-joining method: A new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 1987, 4, 406–425. [Google Scholar]
  33. Tamura, K.; Nei, M.; Kumar, S. Prospects for inferring very large phylogenies by using the neighbor-joining method. Proc. Natl. Acad. Sci. USA 2004, 101, 11030–11035. [Google Scholar] [CrossRef] [PubMed]
  34. Stecher, G.; Tamura, K.; Kumar, S. Molecular evolutionary genetics analysis (MEGA) for macOS. Mol. Biol. Evol. 2020, 37, 1237–1239. [Google Scholar] [CrossRef] [PubMed]
  35. Zhao, Z.; Liu, H.; Xue, Y.; Chen, K.; Liu, Z.; Xue, Q.; Wang, C. Analysis of efficacy obtained with a trivalent inactivated Haemophilus parasuis serovars 4, 5, and 12 vaccine and commercial vaccines against Glässer’s disease in piglets. Can. J. Vet. Res. 2017, 81, 22–27. [Google Scholar]
  36. Department of Veterinary Services. List of Approved Veterinary Vaccines. 2018. Available online: https://www.dvs.gov.my/dvs/resources/user_1/DVS%20pdf/zooka/2018/03.05.2018/List_of_Approved_Veterinary_Vaccine_updated_May_2018.pdf (accessed on 27 January 2023).
  37. Macedo, N.; Gottschalk, M.; Strutzberg-Minder, K.; Van, C.N.; Zhang, L.; Zou, G.; Zhou, R.; Marostica, T.; Clavijo, M.J.; Tucker, A.; et al. Molecular characterization of Glaesserella parasuis strains isolated from North America, Europe and Asia by serotyping PCR and LS-PCR. Vet. Res. 2021, 52, 68. [Google Scholar] [CrossRef] [PubMed]
  38. Singh, M.S.G.; Fong, R.W.J. Swine breeding and production in Malaysia. In Proceedings of the International Symposium on Recent Progress in Swine Breeding and Raising Technologies, Tainan, Taiwan, 3–4 June 2014; pp. 153–166. [Google Scholar]
  39. Choo, L.C. Malaysia Country Report. In Proceedings of the OIE Second Regional Workshop on Swine Disease Control, Beijing, China, 27–29 June 2017. [Google Scholar]
  40. Zhang, B.; Feng, S.; Xu, C.; Zhou, S.; He, Y.; Zhang, L.; Zhang, J.; Guo, L.; Liao, M. Serum resistance in Haemophilus parasuis SC096 strain requires outer membrane protein P2 expression. FEMS Microbiol. Lett. 2012, 326, 109–115. [Google Scholar] [CrossRef]
  41. Scherrer, S.; Rademacher, F.; Stephan, R.; Peterhans, S. Identification of Glaesserella parasuis and Differentiation of Its 15 Serovars Using High-Resolution Melting Assays. Pathogens 2022, 11, 752. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Serotyping using primers by Howell et al. (2015) [4]: (a) with an adopted method of the set 2 recommendation by Schuwerk et al., (2020) [12]. M, molecular weight marker (100 bp); NTC, nontemplate control; 1, GP/UPM/MY001; 2, GP/UPM/MY002; 3, GP/UPM/MY003; 4, GP/UPM/MY004; 5, GP/UPM/MY005; 6, GP/UPM/MY006; 7, GP/UPM/MY007; 8, GP/UPM/MY008; 9, GP/UPM/MY008; 10, GP/UPM/MY010; 11, GP/UPM/MY011; (b) with the set 1 recommendation by Schuwerk et al., (2020) [12]. M, molecular weight marker (100 bp); NTC, nontemplate control; 1, GP/UPM/MY001; 2, GP/UPM/MY002; 3, GP/UPM/MY003; 4, GP/UPM/MY004.
Figure 1. Serotyping using primers by Howell et al. (2015) [4]: (a) with an adopted method of the set 2 recommendation by Schuwerk et al., (2020) [12]. M, molecular weight marker (100 bp); NTC, nontemplate control; 1, GP/UPM/MY001; 2, GP/UPM/MY002; 3, GP/UPM/MY003; 4, GP/UPM/MY004; 5, GP/UPM/MY005; 6, GP/UPM/MY006; 7, GP/UPM/MY007; 8, GP/UPM/MY008; 9, GP/UPM/MY008; 10, GP/UPM/MY010; 11, GP/UPM/MY011; (b) with the set 1 recommendation by Schuwerk et al., (2020) [12]. M, molecular weight marker (100 bp); NTC, nontemplate control; 1, GP/UPM/MY001; 2, GP/UPM/MY002; 3, GP/UPM/MY003; 4, GP/UPM/MY004.
Pathogens 12 00308 g001
Figure 2. A PCR gel photo of vtaA amplification. (a) M, molecular weight marker (100 bp); NTC, nontemplate control; NC, negative control (S. suis); 1–4, samples yielded poor result that were not including in this publication; 5, GP/UP/MY006; (b) M, molecular weight marker (100 bp); NTC, nontemplate control; NC, negative control (S. suis); 1, sample yielded negative result that was not including in this publication; 2, GP/UPM/MY001; 3, GP/UPM/MY002; 4, GP/UPM/MY003; 5, GP/UPM/MY004; 6, GP/UPM/MY005; V; virulent.
Figure 2. A PCR gel photo of vtaA amplification. (a) M, molecular weight marker (100 bp); NTC, nontemplate control; NC, negative control (S. suis); 1–4, samples yielded poor result that were not including in this publication; 5, GP/UP/MY006; (b) M, molecular weight marker (100 bp); NTC, nontemplate control; NC, negative control (S. suis); 1, sample yielded negative result that was not including in this publication; 2, GP/UPM/MY001; 3, GP/UPM/MY002; 4, GP/UPM/MY003; 5, GP/UPM/MY004; 6, GP/UPM/MY005; V; virulent.
Pathogens 12 00308 g002
Figure 3. An electrophoresed gel photo for the amplification of OmpP2 gene (GP/UPM/MY011 not shown here). M, molecular weight marker (100 bp); 1, nontemplate control; 2, GP/UPM/MY001; 3, negative control (Streptococcus suis); 4, GP/UPM/MY002; 5, GP/UPM/MY003; 6, GP/UPM/MY004; 7, GP/UPM/MY005; 8, GP/UPM/MY006; 9, GP/UPM/MY007; 10, GP/UPM/MY008; 11, GP/UPM/MY008; 12, GP/UPM/MY010.
Figure 3. An electrophoresed gel photo for the amplification of OmpP2 gene (GP/UPM/MY011 not shown here). M, molecular weight marker (100 bp); 1, nontemplate control; 2, GP/UPM/MY001; 3, negative control (Streptococcus suis); 4, GP/UPM/MY002; 5, GP/UPM/MY003; 6, GP/UPM/MY004; 7, GP/UPM/MY005; 8, GP/UPM/MY006; 9, GP/UPM/MY007; 10, GP/UPM/MY008; 11, GP/UPM/MY008; 12, GP/UPM/MY010.
Pathogens 12 00308 g003
Figure 4. A phylogenetic tree of the Glaesserella parasuis OmpP2 gene nucleotide sequences. Evolutionary history was inferred using the neighbor-joining method [32]. The optimal tree is shown. The tree is drawn to scale, with branch lengths in the same units as those of the evolutionary distances used to infer the phylogenetic tree. The evolutionary distances were computed using the maximum composite likelihood method [33] and are in the units of the number of base substitutions per site. This analysis involved 90 nucleotide sequences, i.e., 11 Malaysian G. parasuis OmpP2 gene nucleotide sequences and 79 reference sequences from the NCBI GenBank. All ambiguous positions were removed for each sequence pair (pairwise deletion option). There were a total of 1299 positions in the final dataset. Evolutionary analyses were conducted in the MEGA 11 software [29,34]. The GenBank accession numbers, strain, origin country/location and serotype (if available) are as indicated. Malaysian gene sequences were additionally highlighted with a box.
Figure 4. A phylogenetic tree of the Glaesserella parasuis OmpP2 gene nucleotide sequences. Evolutionary history was inferred using the neighbor-joining method [32]. The optimal tree is shown. The tree is drawn to scale, with branch lengths in the same units as those of the evolutionary distances used to infer the phylogenetic tree. The evolutionary distances were computed using the maximum composite likelihood method [33] and are in the units of the number of base substitutions per site. This analysis involved 90 nucleotide sequences, i.e., 11 Malaysian G. parasuis OmpP2 gene nucleotide sequences and 79 reference sequences from the NCBI GenBank. All ambiguous positions were removed for each sequence pair (pairwise deletion option). There were a total of 1299 positions in the final dataset. Evolutionary analyses were conducted in the MEGA 11 software [29,34]. The GenBank accession numbers, strain, origin country/location and serotype (if available) are as indicated. Malaysian gene sequences were additionally highlighted with a box.
Pathogens 12 00308 g004
Table 1. A summary of isolate information and results.
Table 1. A summary of isolate information and results.
Isolate IDStateFarmAnimalYearFarm Size (Sow Number)Postmortem LesionSite of IsolationGenBank Accession NumberSerotypevtaAGenetic Type
GP/UPM/MY001SelangorS1S100120181000Polyserositis, arthritislungOP6510055 or 12VirulentIa
GP/UPM/MY002SelangorS2S20012018300ArthritisjointOP6510063VirulentIa
GP/UPM/MY003SelangorS3S30012018300Fibrinous pneumoniaelungOP6510078VirulentIa
GP/UPM/MY004SelangorS3S30022018300Fibrinous pneumoniaelungOP6510085 or 12VirulentIa
GP/UPM/MY005SelangorS4S40012020900PolyserositisbrainOP6510095 or 12VirulentIb
GP/UPM/MY006SelangorS4S40012020900PolyserositistonsilOP65101015NonvirulentII
GP/UPM/MY007PerakP1P10012020500PolyserositisbrainOP6510115 or 12VirulentIb
GP/UPM/MY008SelangorS5S50012021400Polyserositispleural swabOP6510125 or 12VirulentIb
GP/UPM/MY009PerakP2P200120224000Polyserositis, arthritislungOP6510135 or 12VirulentIb
GP/UPM/MY010PerakP2P200120224000Polyserositis, arthritisperitoneal swabOP6510145 or 12VirulentIb
GP/UPM/MY011PerakP2P200120224000Polyserositis, arthritisjointOP6510155 or 12VirulentIb
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Lee, C.Y.; Ong, H.X.; Tan, C.Y.; Low, S.E.; Phang, L.Y.; Lai, J.; Ooi, P.T.; Fong, M.W.C. Molecular Characterization and Phylogenetic Analysis of Outer membrane protein P2 (OmpP2) of Glaesserella (Haemophilus) parasuis Isolates in Central State of Peninsular Malaysia. Pathogens 2023, 12, 308. https://doi.org/10.3390/pathogens12020308

AMA Style

Lee CY, Ong HX, Tan CY, Low SE, Phang LY, Lai J, Ooi PT, Fong MWC. Molecular Characterization and Phylogenetic Analysis of Outer membrane protein P2 (OmpP2) of Glaesserella (Haemophilus) parasuis Isolates in Central State of Peninsular Malaysia. Pathogens. 2023; 12(2):308. https://doi.org/10.3390/pathogens12020308

Chicago/Turabian Style

Lee, Chee Yien, Hui Xin Ong, Chew Yee Tan, Suet Ee Low, Lai Yee Phang, Jyhmirn Lai, Peck Toung Ooi, and Michelle Wai Cheng Fong. 2023. "Molecular Characterization and Phylogenetic Analysis of Outer membrane protein P2 (OmpP2) of Glaesserella (Haemophilus) parasuis Isolates in Central State of Peninsular Malaysia" Pathogens 12, no. 2: 308. https://doi.org/10.3390/pathogens12020308

APA Style

Lee, C. Y., Ong, H. X., Tan, C. Y., Low, S. E., Phang, L. Y., Lai, J., Ooi, P. T., & Fong, M. W. C. (2023). Molecular Characterization and Phylogenetic Analysis of Outer membrane protein P2 (OmpP2) of Glaesserella (Haemophilus) parasuis Isolates in Central State of Peninsular Malaysia. Pathogens, 12(2), 308. https://doi.org/10.3390/pathogens12020308

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop