Research Article

Genetic Diversity of Pilin From Kazakh Isolates of Moraxella bovoculi

Marat Kuibagarov1, Assylbek Zhylkibayev1, Dinara Kamalova1, Anara Ryskeldina1, Nurdina Yerzhanova1, Yerlan Ramankulov1, Alexandr Shevtsov1*, John A Angelos2

1National Center for Biotechnology, 13/5, Kurgalzhynskoye road, Astana, 010000, Kazakhstan; 2Department of Medicine & Epidemiology, School of Veterinary Medicine, University of California, Davis, CA 95616, USA.

Received | September 20, 2022; Accepted | October 05, 2022; Published | October 15, 2022

*Correspondence | Alexandr Shevtsov, National Center for Biotechnology, 13/5, Kurgalzhynskoye road, Astana, 010000, Kazakhstan; Email: [email protected]

Citation | Kuibagarov M, Zhylkibayev A, Kamalova D, Ryskeldina A, Yerzhanova N, Ramankulov Y, Shevtsov A, Angelos JA (2022). Genetic diversity of pilin from kazakh isolates of moraxella bovoculi. Adv. Anim. Vet. Sci. 10(11): 2376-2383.

DOI | http://dx.doi.org/10.17582/journal.aavs/2022/10.11.2376.2383

ISSN (Online) | 2307-8316

 

 

Copyright: 2022 by the authors. Licensee ResearchersLinks Ltd, England, UK.

This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

Introduction

Infectious bovine keratoconjunctivitis (IBK), is a highly contagious and infectious disease of the bovine eye characterized by conjunctivitis and ulcerative keratitis that occurs worldwide (Holzhauer and Visser, 2004; Kneipp, 2021). Infectious keratoconjunctivitis also occurs in other domestic (Motha et al., 2003) and wild (Giacometti et al., 2002) species and is generally considered to be a multifactorial disease. The most commonly isolated pathogen from IBK is Moraxella bovis (M. bovis) (McConnel et al., 2007). The presence of M. bovoculi isolates in most eye samples taken from IBK-affected cattle suggests a role for M. bovoculi in pathogenesis (Angelos, 2010). The frequency of co-isolation of these two species may indicate that M. bovoculi serves as an opportunistic agent in IBK (Loy and Brodersen, 2014). In a recent study using the PCR method, M. bovoculi was detected in 75% of the samples of the total number of samples examined (Zheng et al., 2019). Other data on a positive correlation between the clinical signs of IBK and the presence of M. bovoculi also indirectly support that this organism is involved in the pathogenesis of IBK (Loy and Brodersen, 2014; Schnee et al., 2015).

Pilin is known to be one of the key pathogenic factors in M. bovis, as this protein allows moraxella to “stick” to the surface of the ocular epithelium (Greene et al., 2001; Lehr et al., 1985; Lepper and Hermans, 1986; Prieto et al., 2003; Ruehl et al., 1993). M.bovis pili (classified as type 4) show antigenic and phase variation and are classified into seven groups (A to G) (Moore and Lepper, 1991). These pili are composed of repeating, homologous polypeptide subunits called pilin. Type 4 pilins contain many variable and hypervariable domains that are a source of both structural and antigenic variation between species (Ruehl et al., 1993). Two different types of pilin molecules, called I and Q, have been identified as separate gene products. In vitro experiments have shown that phase variation allows M.bovis to switch from one type of pilin expression to another (Ruehl et al., 1993). The rather wide diversity of pilin genes among different M.bovis isolates may be the reason for the ineffectiveness of fortified pilin preparations proposed as vaccines against IKC. Sequence heterogeneity of a pilin gene (pilA) in M.bovoculi was also recently reported, however, diversity amongst deduced PilA sequences was much less in comparison to pilin amino acid sequences in M. bovis (Angelos et al., 2021).

That there are comparatively large differences between the amino acid sequences of PilA protein in M. bovoculi isolated from the nasopharynx of cattle without IBK symptoms versus isolates from the surface of the eyes from animals with IBK suggests a role for pilin in establishing M. bovoculi in different environments on the host.

To further expand our understanding of PilA in M. bovoculi and help guide potential development of vaccines to prevent IBK, we determined genetic diversity amongst structural pilin (PilA) sequences encoded by M. bovoculi from cattle with clinical signs of IBK from different regions of Kazakhstan.

Materials And Methods

For this study we used isolates of M. bovoculi that were isolated from cattle from various regions of Kazakhstan during May-August 2021.

Samples of biological material from eyes of cattle with clinical signs of IBK were taken using test tubes with Amies transport medium and plated on Columbia agar with 5% bovine serum and then incubated for 20-24 hours at 37°C. Colonies with characteristic growth for M. bovoculi were subcultured and used for primary identification by Matrix-assisted laser desorption-ionization time of flight mass spectrometry (MALDI-TOF MS). To do this, single colonies were suspended in 1 µl of a saturated solution of α-cyano-4-hydroxycinnamic acid (α-HCCA) with 50% acetonitrile and 2.5% trifluoroacetic acid (TFA), and dried in air. The chip with samples deposited on it was placed in a microflex LT MALDI-TOF mass spectrometer (Bruker Daltonics). The chip was positioned and calibrated against a calibration standard. The spectra were collected automatically using 40 laser pulses (60 Hz). The analyzed mass/charge range was 2000–20000 Da. The spectra were analyzed using the Bruker MALDI-TOF Biotyper v4.0 software; the minimum score for species identification was 1.8.

Isolates were stored at -80°C in trypticase soy broth with 40% glycerol.

Genomic DNA isolation was performed using QIAamp DNA Mini Kit, Qiagen.

The 16S rRNA gene fragment was amplified using universal primers 8f (5’-agagtttgatcctggctcag-3’) and 806R (5’- ggactaccagggtatctaat-3’) in a total volume of 30 µl (Srinivasan et al., 2015). The PCR reaction was performed with universal primers. PCR reaction included 2.5 mM of MgCl2, 300 nM of each primer, 200 nM each dNTP, 1X Taq Buffer with KCl (Thermo Scientific™), 1 U Taq DNA Polymerase (Thermo Scientific™, EP0401) and 3 ng DNA. The PCR program was 3 min at 95°C, 30 cycles at 95°C for 30 s, 55ºC 30 s and 72°C for 50 s, followed by 72°C for 10 min.

The pilA pilin gene was amplified using primers Mbovoc_Pilin_Dn (5-gtggggttacataaatataaaga-3’) and Mbovoc_Pilin_Up3 (5’-gattaatcaaaccttcaaacac-3’) as described (Angelos et al., 2021).

The efficiency of amplification was evaluated using gel electrophoresis. To determine the nucleotide sequence, PCR products were purified from unbound primers by an enzymatic method using Exonuclease I (Thermo Scientific™) and Shrimp Alkaline Phosphatase (Applied Biosystems™) as described (Werle et al., 1994). The sequencing reaction, using the above reaction primers, was performed using the BigDye® Terminator v3.1 Cycle Sequencing Kit (Applide Biosystems) according to the manufacturer’s instructions, followed by separation of the fragments on a 3730xlDNAAnalyzer automated genetic analyzer (Applide Biosystems). Nucleotide sequences were analyzed and combined into a common sequence in SeqMan software (DNAStar).

Nucleotide sequences of 16S rRNA were identified in BLAST (https://blast.ncbi.nlm.nih.gov/) by homologous sequences with the highest percentage of identity (Boratyn et al., 2013), and determination of percent identity with reference sequences of M. bovoculi, M. bovis, M. ovis, and M. lacunata (GeneBank acesion numbers: NR_043583, NR_028668, NR_0286701 and NR_036825, respectively) was calculated in MegAlign (DNAStar).

The nucleotide and amino acid sequence of the pilA gene was used to determine genetic diversity by building phylogenetic trees and percent identity. Percent identity was determined in MegAlign (DNAStar). Phylogenetic trees were built in MEGA X (Kumar et al., 2018). The evolutionary history was inferred using the Neighbor-Joining method (Saitou and Nei, 1987). The optimal tree is shown. The percentage of replicate trees in which the associated taxa clustered together in the bootstrap test (1000 replicates) are shown next to the branches (Felsenstein, 1985). The tree was 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 Kimura 2-parameter method (Kimura, 1980) and are in the units of the number of base substitutions per site. The rate variation among sites was modeled with a gamma distribution (shape parameter = 5). Codon positions included were 1st+2nd+3rd+Noncoding. All ambiguous positions were removed for each sequence pair (pairwise deletion option) (Kumar et al., 2018). The phylogenetic analysis of deduced amino acid sequences included PilA sequences of 10 previously identified PilA groups in M. bovoculi (Angelos et al., 2021).

Results

As a result of bacteriological studies, 35 bacterial isolates were identified on the basis of phenotypic characteristics as M. bovoculi. Using the MALDI-TOF method, all isolates were assigned to the genus Moraxella, however, in 33 isolates the protein profile corresponded to M. bovoculi, and in two isolates the protein profiles were identified as Moraxella_sg_Branhamella ovis, with a high Score Value (1.7-1.9). Genetic identification by nucleotide sequence analysis of the 16SrRNA gene identified all isolates as M. bovoculi (Table 1). The percentage of identity with the NR_043583 sequence of M. bovoculi was 99.8-100%, with M. bovis - 98.4-98.8%, with M. ovis - 95.7-96.1% and with M. lacunata - 99.4-99.6%. The isolates were from 4 farms in Kostanay, Akmola and Pavlodar regions (Figure 1).

 

Based on the nucleotide sequence of the pilA gene, the analyzed isolates were grouped into 10 genotypes (Figure 2) designated gNA1 through gNA10.

 

It was determined that 3 genotypes were represented by single isolates (gNA3, gNA5, gNA9) which were isolated from Akmola and Pavlodar regions. The gNA1 genotype contained 5 isolates from the Kostanay region. The largest gNA2 genotype represented 11 isolates: 7 were isolated on a farm in the city of Arkalyk, Kostanay region, 3 from the Bulandy district of the Akmola region, and 1 from the Akkuli district of the Pavlodar region. The gNA4 genotype includes 3 isolates from the Bulandy district of the Akmola region. The gNA6 genotype was represented by two isolates from the Bulandy district of the Akmola region, which, based on MALDI-TOF, had been identified as Branhamella ovis. The gNA7 genotype included 4 isolates from the Birzhan sal district of the Akmola region and 2 isolates from the Kostanay region. The gNA8 genotype contained 2 isolates from the Bulandy district of the

 

Table 1: Results of isolate identification

Sample number	Region	Date of sampling	Breed of cattle	Identification result
				MALDI-TOF)	16S rRNA
1	Kostanay region, Arkalyk	11.05.21	Kazakh white-headed	Moraxella_sg_Moraxella bovoculi	Moraxella bovoculi
2	Kostanay region, Arkalyk	11.05.21	Kazakh white-headed	Moraxella_sg_Moraxella bovoculi	Moraxella bovoculi
3	Kostanay region, Arkalyk	11.05.21	Kazakh white-headed	Moraxella_sg_Moraxella bovoculi	Moraxella bovoculi
4	Kostanay region, Arkalyk	11.05.21	Kazakh white-headed	Moraxella_sg_Moraxella bovoculi	Moraxella bovoculi
5	Kostanay region, Arkalyk	11.05.21	Kazakh white-headed	Moraxella_sg_Moraxella bovoculi	Moraxella bovoculi
6	Kostanay region, Arkalyk	11.05.21	Kazakh white-headed	Moraxella_sg_Moraxella bovoculi	Moraxella bovoculi
7	Kostanay region, Arkalyk	11.05.21	Kazakh white-headed	Moraxella_sg_Moraxella bovoculi	Moraxella bovoculi
8	Kostanay region, Arkalyk	11.05.21	Kazakh white-headed	Moraxella_sg_Moraxella bovoculi	Moraxella bovoculi
9	Kostanay region, Arkalyk	11.05.21	Kazakh white-headed	Moraxella_sg_Moraxella bovoculi	Moraxella bovoculi
10	Kostanay region, Arkalyk	11.05.21	Kazakh white-headed	Moraxella_sg_Moraxella bovoculi	Moraxella bovoculi
11	Kostanay region, Arkalyk	11.05.21	Kazakh white-headed	Moraxella_sg_Moraxella bovoculi	Moraxella bovoculi
12	Kostanay region, Arkalyk	11.05.21	Kazakh white-headed	Moraxella_sg_Moraxella bovoculi	Moraxella bovoculi
13	Kostanay region, Arkalyk	11.05.21	Kazakh white-headed	Moraxella_sg_Moraxella bovoculi	Moraxella bovoculi
14	Kostanay region, Arkalyk	11.05.21	Kazakh white-headed	Moraxella_sg_Moraxella bovoculi	Moraxella bovoculi
15	Akmola region, Bulandy district	31.05.21	Kazakh white-headed	Moraxella_sg_Moraxella bovoculi	Moraxella bovoculi
16	Akmola region, Bulandy district	31.05.21	Kazakh white-headed	Moraxella_sg_Moraxella bovoculi	Moraxella bovoculi
17	Akmola region, Bulandy district	31.05.21	Kazakh white-headed	Moraxella_sg_Moraxella bovoculi	Moraxella bovoculi
18	Akmola region, Bulandy district	31.05.21	Kazakh white-headed	Moraxella_sg_Moraxella bovoculi	Moraxella bovoculi
19	Akmola region, Bulandy district	31.05.21	Kazakh white-headed	Moraxella_sg_Moraxella bovoculi	Moraxella bovoculi
20	Akmola region, Bulandy district	31.05.21	Kazakh white-headed	Moraxella_sg_Branhamella ovis	Moraxella bovoculi
21	Akmola region, Bulandy district	31.05.21	Kazakh white-headed	Moraxella_sg_Moraxella bovoculi	Moraxella bovoculi
22	Akmola region, Bulandy district	31.05.21	Kazakh white-headed	Moraxella_sg_Moraxella bovoculi	Moraxella bovoculi
23	Akmola region, Bulandy district	31.05.21	Kazakh white-headed	Moraxella_sg_Branhamella ovis	Moraxella bovoculi
24	Akmola region, Bulandy district	31.05.21	Kazakh white-headed	Moraxella_sg_Branhamella ovis	Moraxella bovoculi
25	Akmola region, Bulandy district	31.05.21	Kazakh white-headed	Moraxella_sg_Moraxella bovoculi	Moraxella bovoculi
26	Akmola region, Birzhan sal district	9.08.21	Angus	Moraxella_sg_Moraxella bovoculi	Moraxella bovoculi
27	Akmola region, Birzhan sal district	9.08.21	Angus	Moraxella_sg_Moraxella bovoculi	Moraxella bovoculi
28	Akmola region, Birzhan sal district	9.08.21	Angus	Moraxella_sg_Moraxella bovoculi	Moraxella bovoculi
29	Akmola region, Birzhan sal district	9.08.21	Angus	Moraxella_sg_Moraxella bovoculi	Moraxella bovoculi
30	Akmola region, Birzhan sal district	9.08.21	Angus	Moraxella_sg_Moraxella bovoculi	Moraxella bovoculi
31	Pavlodar region, Akkuli district	19.08.21	Angus	Moraxella_sg_Moraxella bovoculi	Moraxella bovoculi
32	Pavlodar region, Akkuli district	19.08.21	Angus	Moraxella_sg_Moraxella bovoculi	Moraxella bovoculi
33	Pavlodar region, Akkuli district	19.08.21	Angus	Moraxella_sg_Moraxella bovoculi	Moraxella bovoculi
34	Pavlodar region, Akkuli district	19.08.21	Angus	Moraxella_sg_Moraxella bovoculi	Moraxella bovoculi
35	Pavlodar region, Akkuli district	19.08.21	Angus	Moraxella_sg_Moraxella bovoculi	Moraxella bovoculi

 

Table 2: Amino acid substitutions identified in the PilA protein

Isolate	Position
	33	61	62	81	98	130	131
10 Kost Arkalyk gNA1	R	S	N	V	I	G	Q
8 Kost Arkalyk gNA2	R	S	N	V	I	G	Q
18 Akmol Bulandy gNA3	R	S	N	V	I	G	Q
19 Akmol Bulandy gNA4	R	N	N	V	I	G	Q
28 Akmol Birzhan sal gNA5	R	S	N	V	I	G	Q
23 Akmol Bulandy gNA6	R	S	N	V	V	G	Q
30 Akmol Birzhan sal gNA7	R	S	N	V	I	G	Q
15 Akmol Bulandy gNA8	*	S	N	V	I	S	Q
33 Pavl Akkuli gNA9	R	S	N	V	I	D	K
31 Pavl Akkuli gNA10	R	S	S	P	I	G	Q

 

Akmola region. The gNA10 genotype was represented by 3 isolates from the Pavlodar region. Identity between genotypes varied from 99.8% to 96.1%, which, in nucleotide equivalents, was from 1 to 17 nucleotides (Figure 3).

Most of the variable nucleotides were silent mutations and did not lead to protein amino acid substitutions. A total of 8 nucleotides were missense mutations and resulted in amino acid substitutions (Table 2). The maximum difference between the genogroups was 2 amino acid residues out of 152.

The deduced PilA amino acid sequence from genotypes gNA1, gNA2, gNA3, gNA5 and gNA7 did not differ from each other. In the gNA4 genotype, serine was replaced by asparagine at position 61 (S61N). In the gNA6 genotype, at position 98 of the amino acid sequence, isoleucine was replaced by valine (I98V). In the gNA8 genotype, the substitution of cytosine for thymine at position 97 of the nucleotide sequence of the pilA gene led to the formation of a stop codon, which would result in the expression of a truncated protein in the 15_Akmol_Bulandy and 17_Akmol_Bulandy isolates. Also in this genotype, a substitution of glycine for serine at position 130 (G130S) was identified. In the gNA9 genotype, 2 amino acid substitutions were found at positions 130 and 131 of glycine for aspartic acid and glutamine for lysine, respectively (G130D and Q131K). In the gNA10 genotype, two amino acid substitutions of aspargine for serine at position 62 of the protein and valine for proline (N62S and V81P) were also found.

Phylogenetic analysis of the PilA deduced amino acid sequence of isolates from Kazakhstan with previously established groups made it possible to establish the group affiliation of the PilA gene genotypes (Figure 4). It was determined that more than 68% of the isolates (N=24) belonged to PilA group A. Three isolates each belonged to PilA groups C and F. Two isolates from the Akmola region under accession numbers 22 and 23 have a unique amino acid sequence and were assigned to a new group designated PilA group K. The two truncated PilA sequences (gNA8; isolates 15 and 17) were located in a separate clade and were assigned to an additional new group designated PilA group L. Isolate number 33 was also located separately and was designated as belonging to a new group designated PilA group M.

Nucleotide sequences of pilA gene has been deposited in GenBank under accession numbers OP352320 - OP352329.

 

Discussion

According to the conclusion of a number of researchers, IBK is one of the infections that cause the greatest economic damage to livestock (Brown et al., 1998; Prieto et al., 2013; Postma et al., 2008; Dennis et al., 2021). In the last decade in Kazakhstan, the number of cases of IBK has increased, which are primarily associated with the import of highly productive animals, the infection rate reaches 27% (Ivanov et al., 2021). Over the past few decades, numerous studies have been carried out to study the effectiveness of various types of vaccines for the prevention of IBK. In most cases, the emphasis was on whole-cell bacterial M bovis/M. bovoculi preparations and pilin antigens. The current lack of universally effective IBK vaccines suggests that further research is needed on pathogenicity and immunogenicity of Moraxella spp associated with IBK in cattle (Hille et al., 2022; Funk et al., 2009; O’Connor et al., 2019). The identification of conserved nucleotide sequences of genes encoding proteins with immunogenic potential is an important step towards creating improved vaccine preparations (Acquistapace et al., 2021). The data presented here should help guide future studies aimed at identifying more effective pilin-based vaccines to prevent IBK.

This is the first study in Kazakhstan to study the genetic homogeneity of M. bovoculi genes associated with IBK. Whether or not the differences that we identified in deduced PilA sequences correspond to differences in pathogenic potential of M. bovoculi remains to be determined. Further studies that examine temporal changes in PilA group circulating in herds of cattle may be needed to help answer this question.

Genetic identification showed that all 35 isolates from cattle in the farms of Kostanay, Akmola and Pavlodar regions were M. bovoculi. Some researchers suggest the lack of pigmentation around the eyes in animals of some breeds (Hereford and others) as one of the predisposing factors of IBK (Ward and Nielson, 1979). In this study, M. bovoculi isolates were from cattle with clinical signs of IBK of the Kazakh Whitehead (lack of pigmentation) as well as Angus breeds. Thus, in the cattle from which these isolates were derived, we did not observe any clear association between IBK and pigmentation.

As a result of the sequencing of the pilA gene sequences of the M. bovoculi isolates, it was determined that they belonged to 10 different genotypes, the identity between which varied from 96.1-99.8%. Due to the replacement of cytosine with thymine in the 97th position of the gNA8 genotype, a stop codon was formed, which would result in expression of a defective PilA protein in two isolates. PilA has been identified as a major virulence factor in some Gram-negative bacteria (Forslund et al., 2010), however, one previous study reported no effect on the pathogenicity of Pseudomonas aeruginosa isolates by turning off the pilA gene (Choi et al., 2002). Phylogenetic analysis of the amino acid sequence of PilA from M. bovoculi isolates from Kazakhstan, with previously established groups, made it possible to determine that 68% of Kazakh isolates belonged to group A, which is consistent with previous studies indicating a high prevalence of this group amongst strains isolated in the USA (Angelos et al., 2021). Groups C and F included 8.6% of isolates. Some of the studied isolates had a unique amino acid sequence and can be assigned to 3 additional PilA groups, K, L, and M.

In general, there is a high degree of homogeneity in the deduced amino acid sequence between the isolates included in the 13 established groups A-L.

More extensive and systematic studies will be needed to determine if any attempts by M.bovoculi to evade the host’s immune system is reflective of differences in PilA sequences.

Conclusion

The pilA gene from Kazakh isolates of M. bovoculi had a high degree of genetic homogeneity. PilA may be a worthwhile vaccine antigen to consider when creating vaccines to prevent IBK. Nevertheless, the identification of a putative truncated PilA in two isolates from animals with IBK requires further investigation of its role in pathogenesis.

Funding

The study was carried out with the financial support of the Ministry of Education and Science of the Republic of Kazakhstan within the framework of project AP08856254 for 2020-2022.

acknowledgements

The authors are grateful to Alexander Suminov for his great help in collecting samples.

Conflict of interest

The authors declare that there is no conflict of interest.

novelty statement

Sequencing of pilA in Kazakh isolates of M. bovoculi revealed 3 novel PilA deduced amino acid sequences. Additional research is needed to determine if PilA is involved in the biology of M. bovoculi as it relates to the presence of this organism in eyes of cattle with infectious bovine keratoconjunctivitis.

Authors contribution

This manuscript’s materials preparation, data analysis, and text writing were all done by all of the authors.

References

Acquistapace S., Umpiérrez A., Fernández-Ciganda S., Sotelo-Silveira J., Zunino P. (2021). Proteína de membrana externa CD de Moraxella bovis como antígeno potencial contra la queratoconjuntivitis infecciosa bovina. Vet. (Montevideo)., 57(216). https://doi.org/10.29155/vet.57.216.4.

Angelos JA. (2010). Moraxella bovoculi and Infectious Bovine Keratoconjunctivitis: Cause or Coincidence? Vet. Clin. North America-Food Anim. Pract. 26(1):73-+. Https://doi.org/10.1016/j.cvfa.2009.10.002

Angelos JA, Clothier KA, Agulto RL, Mandzyuk B, Tryland M. (2021). Relatedness of type IV pilin PilA amongst geographically diverse Moraxella bovoculi isolated from cattle with infectious bovine keratoconjunctivitis. J. Med. Microbiol. 70(2). Https://doi.org/ 10.1099/jmm.0.001293

Boratyn G. M., Camacho C., Cooper P. S., Coulouris G., Fong A., Ma N., Madden TL,  Matten WT,  McGinnis SD,  Merezhuk Y,  Raytselis Y,  Sayers EW,  Tao T,  Ye J,  Zaretskaya I. (2013). BLAST: a more efficient report with usability improvements. Nucl. Acids Res., 41(W1): W29-W33. Https://doi.org/ 10.1093/nar/gkt282

Brown MH, Brightman AH, Fenwick BW, Rider MA. (1998). Infectious bovine keratoconjunctivitis: a review. J. Vet. Intern Med. 12(4):259-66. Https://doi.org/10.1111/j.1939-1676.1998.tb02120.x

Choi J. Y., Sifri C. D., Goumnerov B. C., Rahme L. G., Ausubel F. M., Calderwood S. B. (2002). Identification of virulence genes in a pathogenic strain of Pseudomonas aeruginosa by representational difference analysis. J. Bacteriol., 184(4): 952-961. Https://doi.org/10.1128/jb.184.4.952-961.2002

Dennis EJ, Kneipp M. (2021). A Review of Global Prevalence and Economic Impacts of Infectious Bovine Keratoconjunctivitis. Vet. Clin. North Am. Food Anim. Pract. 37(2):355-69. Https://doi.org/10.1016/j.cvfa.2021.03.010

Felsenstein J. (1985). Confidence limits on phylogenies: an approach using the bootstrap. Evolution. 39(4):783-91. Https://doi.org/org/10.2307/2408678

Forslund A. L., Salomonsson E. N., Golovliov I., Kuoppa K., Michell S., Titball R., Oyston P., Noppa L., Sjöstedt A., Forsberg Å. (2010). The type IV pilin, PilA, is required for full virulence of Francisella tularensis subspecies tularensis. BMC Microbiol. 10(1): 1-8. Https://doi.org/ 10.1186/1471-2180-10-227

Funk L., O’Connor A. M., Maroney M., Engelken T., Cooper V. L., Kinyon J., Plummer P. (2009). A randomized and blinded field trial to assess the efficacy of an autogenous vaccine to prevent naturally occurring infectious bovine keratoconjunctivis (IBK) in beef calves. Vaccine., 27(34): 4585-4590. Https://doi.org/ 10.1016/j.vaccine.2009.05.082

Giacometti M, Janovsky M, Belloy L, Frey J. (2002). Infectious keratoconjunctivitis of ibex, chamois and other Caprinae. Revue Scientifique et Technique-Office International des Epizooties.; 21(1):335-46. Https://doi.org/10.20506/rst.21.2.1338

Greene WH, Grubbs ST, Potgieter LN. (2001). Detection of shared magnetic antigenic determinants on whole Moraxella bovis pili by use of antisera to cyanogen bromide-cleaved M. bovis pilus protein. Am. J. Vet. Res. 62(8):1279-84. Https://doi.org/ 10.2460/ajvr.2001.62.1279

Hille M. M., Spangler M. L., Clawson M. L., Heath K. D., Vu H. L., Rogers R. E., Loy J. D. (2022). A Five Year Randomized Controlled Trial to Assess the Efficacy and Antibody Responses to a Commercial and Autogenous Vaccine for the Prevention of Infectious Bovine Keratoconjunctivitis. Vaccines., 10(6): 916. Https://doi.org/10.3390/vaccines10060916

Holzhauer M, Visser I. (2004). Infectious bovine keratoconjunctivitis (IBK) in cows, clinical and lab review at four farms. Tijdschrift voor diergeneeskunde. 129(16):526-9.

Ivanov N. P., Bakiyeva F. A., Namet A. M., Sattarova R. S., Issakulova B. Z., Akmyrzayev N. Z. (2021). The epizootic situation of cattle moraxellosis in several economic entities of the Republic of Kazakhstan. Vet. World., 14(5): 1380. Https://doi.org/ 10.14202/vetworld.2021.1380-1388

Kimura M. (1980). A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J. Molecul. Evol. 16(2):111-20.  Https://doi.org/10.1007/BF01731581

Kneipp M. (2021). Defining and Diagnosing Infectious Bovine Keratoconjunctivitis. Vet. Clin. North Am. Food Anim. Pract. 37(2):237-52.  Https://doi.org/10.1016/j.cvfa.2021.03.001

Kumar S, Stecher G, Li M, Knyaz C, Tamura K. (2018). MEGA X: molecular evolutionary genetics analysis across computing platforms. Molecul. Biol, Evol. 35(6):1547.  Https://doi.org/ 10.1093/molbev/msy096

Loy JD, Brodersen BW. (2014). Moraxella spp. isolated from field outbreaks of infectious bovine keratoconjunctivitis: a retrospective study of case submissions from 2010 to 2013. J. Vet. Diagnost. Investigat. 26(6):761-8. Https://doi.org/10.1177/1040638714551403

Lehr C, Jayappa HG, Goodnow RA. (1985). Serologic and protective characterization of Moraxella bovis pili. Cornell Vet. 75(4):484-92.

Lepper A, Hermans I. (1986). Characterisation and quantitation of pilus antigens of Moraxella bovis by ELISA. Australian Vet. J. 63(12):401-5.

McConnel C, Shum L, House J. (2007). Infectious bovine keratoconjunctivitis antimicrobial therapy. Australian Vet. J..; 85(1‐2):65-9. Https://doi.org/10.1111/j.1751-0813.2006.00080.x

Moore LJ, Lepper AW. (1991). A unified serotyping scheme for Moraxella bovis. Vet. Microbiol. 29(1):75-83. Https://doi.org/10.1016/0378-1135(91)90111-r

Motha M, Frey J, Hansen M, Jamaludin R, Tham K. (2003). Detection of Mycoplasma conjunctivae in sheep affected with conjunctivitis and infectious keratoconjunctivitis. New Zealand Vet. J. 51(4):186-90. Https://doi.org/10.1080/00480169.2003.36362

O’Connor A., Cooper V., Censi L., Meyer E., Kneipp M. (2019). A 2-year randomized blinded controlled trial of a conditionally A 2 year randomized blinded controlled trial of a conditionally licensed Moraxella bovoculi vaccine to aid in prevention of infectious bovine keratoconjunctivitis in Angus beef calves. Medicine., 10: 12-2019. Https://doi.org/10.1111/jvim.15633

Postma GC, Carfagnini JC, Minatel L. (2008). Moraxella bovis pathogenicity: an update. Comp. Immunol. Microbiol. Infect. Dis. 31(6):449-58. Https://doi.org/10.1016/j.cimid.2008.04.001

Prieto CI, Rodriguez MaE, Bosch A, Chirdo FG, Yantorno OM. (2003). Whole-bacterial cell enzyme-linked immunosorbent assay for cell-bound Moraxella bovis pili. Vet. Microbiol. 91(2-3):157-68. Https://doi.org/10.1016/s0378-1135(02)00297-3

Prieto C, Serra DO, Martina P, Jacobs M, Bosch A, Yantorno OM. (2013). Evaluation of biofilm-forming capacity of Moraxella bovis, the primary causative agent of infectious bovine keratoconjunctivitis. Vet. Microbiol. 166(3-4):504-15. Https://doi.org/10.1016/j.vetmic.2013.06.003

Ruehl WW, Marrs C, Beard MK, Shokooki V, Hinojoza JR, Banks S, Bieber D,  Mattick J S (1993). Q pili enhance the attachment of Moraxella bovis to bovine corneas in vitro. Molecul. Microbiol. 7(2):285-8. Https://doi.org/10.1111/j.1365-2958.1993.tb01119.x

Saitou N, Nei M. (1987). The neighbor-joining method: a new method for reconstructing phylogenetic trees. Molecul. Biol. Evol. 4(4):406-25. Https://doi.org/10.1093/oxfordjournals.molbev.a040454

Schnee C, Heller M, Schubert E, Sachse K. (2015). Point prevalence of infection with Mycoplasma bovoculi and Moraxella spp. in cattle at different stages of infectious bovine keratoconjunctivitis. Vet. J. 203(1):92-6. Https://doi.org/10.1016/j.tvjl.2014.11.009

Srinivasan R, Karaoz U, Volegova M, MacKichan J, Kato-Maeda M, Miller S, NadarajanR,  Brodie E L,  Lynch  S V (2015). Use of 16S rRNA gene for identification of a broad range of clinically relevant bacterial pathogens. PLoS One. 10(2):e0117617. Https://doi.org/10.1371/journal.pone.0117617

Ward JK, Nielson MK. (1979). Pinkeye (Bovine Infectious Keratoconjunctivitis) in Beef-Cattle. J. Anim. Sci. 49(2):361-6. Https://doi.org/10.2527/jas1979.492361x

Werle E., Schneider C., Renner M., Völker M., Fiehn W. (1994). Convenient single-step, one tube purification of PCR products for direct sequencing. Nucl. Acids Res., 22(20): 4354. Https://doi.org/10.1093/nar/22.20.4354

Zheng W, Porter E, Noll L, Stoy C, Lu N, Wang Y, Liu X.,  Purvis T.,  Peddireddi L.,  Lubbers B.,  Hanzlicek G.,  Henningson J.,  Liu Z.,  Bai J. (2019). A multiplex real-time PCR assay for the detection and differentiation of five bovine pinkeye pathogens. J. Microbiol. Method.; 160:87-92. Https://doi.org/10.1016/j.mimet.2019.03.024