Qing Liu1,2†
Qingke Kong- 1College of Animal Science and Technology, Southwest University, Chongqing, China
- 2Institute of Preventive Veterinary Medicine, College of Veterinary Medicine, Sichuan Agricultural University, Chengdu, China
Pasteurella multocida is a pathogenic microorganism that causes a variety of serious diseases in humans and animals worldwide. The global regulator gene, fur, plays an important role in pathogenesis and regulates the virulence of many bacteria. Here, we identified a fur gene in P. multocida by complementing a Salmonella Choleraesuis Δfur mutant, and characterized a fur mutant strain of P. multocida. The P. multocida Δfur mutant strain exhibited no significant differences in growth and outer membrane protein (OMP) profiles when the complemented strain was compared to the parent. Ducks were used as the model organism to determine the virulence and protection efficacy induced by Δfur mutant strain. Animal experiments showed that colonization by the mutant was decreased by oral infection of live Δfur mutant strain. The LD50 of the ducks infected with the Δfur mutant was 146-fold higher than that of the ducks infected with the wild-type strain when administered through the oral route. Evaluation of the immunogenicity and protective efficacy of the Δfur mutant of P. multocida revealed strong serum IgY and bile IgA immune responses following oral inoculation with the Δfur strain. Ducks that were orally inoculated with the Δfur mutant strain demonstrated 62% protection efficacy against severe lethal challenge with the wild-type P. multocida. This study provides new insights into P. multocida virulence and the potential use of an attenuated vaccine against P. multocida.
Introduction
Iron is an essential limiting element in bacterial metabolism, and its availability plays a vital role in the pathogenesis of most bacteria (1). Many enzymes require iron as a necessary co-factor for a wide variety of physiological processes (2), such as respiration, DNA biosynthesis, gene regulation, and the tricarboxylic acid (TCA) cycle (3). Proteins such as iron-containing heme and iron-sulfur cluster proteins act not only as excellent electron carriers but also as environmental or intracellular sensors that regulate gene expression (4). However, excess iron is highly toxic to cells. In particular, iron is an enhancer of oxygen toxicity because this metal efficiently catalyzes the Fenton reaction, which converts hydrogen peroxide to the highly reactive hydroxyl radical (5). Therefore, during the process of iron uptake, bacteria maintain a dynamic balance of iron ions within an ideal range by adjusting rates of uptake. Precise regulation of the intracellular iron concentration is typically closely associated with iron acquisition genes, which are controlled via the ferric uptake regulator (Fur) protein. The Fur protein is a metal ion-dependent transcription factor that is ubiquitously expressed in bacteria and controls the transcription and expression of genes involved in a variety of cellular functions, including iron uptake (6). Fur binding to metal ions causes a conformational change in the regulator, which promotes and stabilizes its interaction with DNA (5). Fur-Fe2+ complexes bind to the DNA promoter regions (Fur boxes) of genes involved in iron acquisition, resulting in transcriptional repression or activation (7). In certain bacteria, Fur regulates biofilm formation (8–10) and plays an important role in disease pathogenesis, host-bacteria interactions that control virulence factor expression, and the in vitro regulation of toxin production. The absence or impairment of the fur gene induces decreased virulence of Staphylococcus aureus, Listeria monocytogenes, and Vibrio cholerae in mouse models (11–14). Moreover, in fish models, an Aeromonas salmonicida fur mutant exhibited significantly reduced virulence compared to the wild-type strain, demonstrating its potential for use as a live-attenuated vaccine (1). Pasteurella multocida is a Gram-negative, capsulated, non-motile, highly common bacterium that infects a range of animals through respiratory or oral routes (15), resulting in significant economic losses to industries worldwide (16). Currently, these bacteria are divided into 5 serotypes (A, B, D, E, and F) based on capsular typing and further classified into 16 serotypes (1–16) according to the type of lipopolysaccharide (LPS); different serotypes exhibit varying degrees of virulence in animals (17). Due to the lack of effective coverage in multi-serotype vaccines, antibiotics remain the most commonly employed treatment for avian cholera despite their correlation to increased drug resistance and food safety risks of contaminating bacteria (18). Furthermore, existing live-attenuated vaccines have been observed to revert to virulent strains (19). Therefore, the development of an efficient and safe live-attenuated vaccine for cholera control in fowl is urgently needed. In the present study, we characterized the P. multocida Fur protein. Specifically, the previously putative fur gene of P. multocida was identified through heterologous complementation of P. multocida fur in a Salmonella Choleraesuis Δfur mutant. A Δfur mutant of P. multocida, P. multocida 0818, isolated from a duck was produced. We investigated several phenotypes of the wild-type and Δfur mutant strains including growth curves, outer membrane protein (OMP) profiles, tolerance to polymyxin B and bile salt deoxycholate (DOC), and serum sensitivity. We also examined virulence colonization levels and immunogenicity of fur mutant in a natural host, the Sheldrake duck.
Materials and Methods
Bacterial Strains, Plasmids, Media, Growth Conditions, and Reagents.
The bacterial strains and plasmids used in the present study are listed in Table 1. The P. multocida 0818 strain was isolated from a pasteurellosis outbreak at a duck farm in Sichuan Province, China in 2012 (21). Because the genome of P. multocida 0818 has not been sequenced, all of the primers designed for P. multocida 0818 in Table 2 were based on the genomic sequence of P. multocida Pm70 (AE004439.1) (25). P. multocida was cultured in Bacto brain heart infusion (BHI) broth (BD Bioscience, San Jose, CA, USA), and Salmonella Choleraesuis and Escherichia coli were grown in Luria-Bertani (LB) broth. The antibiotics and reagents used in this study were purchased from Sigma-Aldrich (St. Louis, MO, USA). LB broth, BHI broth, and chromoazurol S (CAS) broth were routinely used (26). When required, the media were supplemented with 1.5% agar, kanamycin (Kan) (50 μg/ml), or chloramphenicol (Cm) (25 μg/ml for S. Choleraesuis and E. coli, 2.5 μg/ml for P. multocida) and ampicillin (Amp) (100 μg/ml). Diaminopimelic acid (DAP) was added to the LB broth (50 μg/ml) for E. coli growth because the χ7213 strain needs DAP to compensate for the Δasd mutation. Bacterial growth was spectrophotometrically detected at a wavelength of 600 nm.
Table 1. Bacterial strains and plasmids used in the present study.
Table 2. Primers used in the present study.
Molecular and Genetic Procedures
Restriction enzymes were obtained from New England BioLabs (Ipswich, MA, USA). PrimeSTAR Max DNA polymerase (Takara Bio Inc., Shiga, Japan) and Taq DNA polymerase (Vazyme Biotech Co., Ltd., Nanjing City, China) were used in all PCR assays. The primers employed in the present study are listed in Table 2. When required, a PCR purification kit and a Plasmid Mini Preparation Kit purchased from Tiangen Biotech (Beijing, China) were used, and the commercial sequencing was performed at BGI Tech (BGI Tech Solutions Co., Ltd., Shenzhen, China).
Detection of Secreted Siderophores
The fur gene of P. multocida 0818 was amplified and cloned into the low-copy-number plasmid pSS664 using two primers, Cfur-F/Cfur-R, which contained BamHI and HindIII cleavage sites, resulting in the recombinant plasmid pSS908. This plasmid was transformed into a Salmonella Choleraesuis Δfur mutant to generate the heterologous expression strain S984 (pSS908). Because the Fur protein is a global regulator, which negatively regulates iron uptake-related genes involved in siderophore production in Salmonella, the siderophore activity will reflect the status of Fur synthesis in Salmonella. Therefore, detecting siderophore production will be indicative of Fur synthesis in Salmonella (27). The production of siderophores was detected using chemical assays as previously described (26). The approach was based on the ability of a siderophore with a high affinity for iron (III). An indicator dye with chrome azurol S (CAS) and hexadecyltrimethylammonium bromide serves as an indicator. The CAS indicator solution consists of 60.5 mg of Chrome Azurol S, 27 mg of Ferric Chloride (FeCl3·6H2O), and 83.3 μl of concentrated hydrochloric acid (HCl) in 200 ml H2O. The color of this indicator dye is blue in the presence of iron (III); however, when a strong chelator, such as siderophore, removes the iron from this dye, it turns from blue to yellow/orange(26). Fifty milliliter CAS indicator solution was added to 1 L minimal basal media agar to make CAS indicator plates (26). The assays were performed by streaking the bacterial strain grown under iron-limiting culture on the CAS agar plate, after overnight inoculation (about 18 h) in 37°C condition, a yellow-orange around the bacterial colony is visualized, indicating bacterial siderophore production, which is negatively regulated by Fur, and the observed blue color around bacterial indicates absence of siderophore.
Expression and Purification of Fur Protein and Generation of a Polyclonal Mouse Anti-Fur Antibody
To obtain Fur protein, the fur gene was amplified from P. multocida 0818 genomic DNA using two primers, pET-fur-F/pET-fur-R, which contained BamHI and NdeI enzyme cleavage sites and six histidine tags in the C-terminus. The complete fur fragment was amplified, then digested using the aforementioned restriction enzymes, and subsequently inserted into the BamHI-NdeI-digested linear pET-32a expression vector, generating pET-32a-fur. The genetic elements of the recombinant plasmid, pET-32a-fur, were confirmed through sequencing and subsequently transformed into competent E. coli Rosetta cells with ampicillin selection. Fur protein expression was induced using different concentrations of isopropyl β-D-thioglucopyranoside (IPTG; Sigma-Aldrich, St. Louis, MO). The Fur protein was characterized by 12.5% SDS polyacrylamide gel electrophoresis (SDS-PAGE). Subsequently, 6 × His/Ni-NTA affinity chromatography was employed for Fur protein purification following the product procedure. To prepare a polyclonal mouse anti-Fur antibody, five female Kunming mice were immunized twice at a 2-week interval through intramuscular injection with the Fur protein (1 mg) and QuickAntibody-Mouse3 W adjuvant (Bio Dragon Immunotechnologies Co., Ltd., Beijing, China). Subsequently, a blood sample was collected from the eye socket, and mouse antisera were prepared 3 weeks after the second immunization. This antibody was used for detecting Fur expression in the wild-type strain, the Δfur mutant S985, and S985 (pSS911) using the western blotting procedure was performed as previously described (28).
Construction of a P. multocida fur Deletion Strain
To construct a fur deletion mutant strain, 500-bp upstream and downstream fragments of the fur gene were amplified using P. multocida 0818 genomic DNA as a template, with the following primers: Dfur-1F/Dfur-1R and Dfur-2F/Dfur-2R (Table 2), respectively. The primers Dfur-1R and Dfur-2F contained NotI and SbfI enzyme cleavage sites. The fused DNA fragment was obtained by amplifying the upstream and downstream fragments using a PCR with the Dfur-1F/Dfur-2R primers. The suicide plasmid T-vector pRE112 was digested with AhdI to generate a linear T-vector (23), and the fused fragment was subsequently ligated into this T-vector, followed by the terminal addition of a polyA tail using the Tailing-A Reaction Kit (Tiangen Biotech) to generate the plasmid pSS909. Subsequently, the KanR cassette sequence, which is 905 bp in length, was amplified by PCR using the primers Kan-F/Kan-R and the plasmid pET28a as a template, followed by a restriction enzyme digestion with NotI and SbfI. The resulting DNA fragment was inserted into pSS909 to construct the recombinant plasmid pSS910. pSS910 was subsequently transformed into E. coli χ7213 and mobilized into P. multocida 0818 through intergeneric transfer using a filter mating method (21). Single crossover integrates were selected after growth on BHI agar containing 50 μg/ml Kan without DAP (to select against E. coli). Successful integration was confirmed by PCR analysis of genomic DNA from the fur mutant strain (21). The fur gene of S. Choleraesuis was deleted by suicide plasmid pYA4807 following the standard method in our lab (29), the resulting mutant strain was named as S984.
To construct complementary fur mutant strains in P. multocida, the fur gene sequence (441 bp in length) was amplified from P. multocida 0818 genomic DNA using two primers, Pfur-F/Pfur-R, which contain BamHI and NotI cleavage sites. The fur fragment was cloned into the pMC-Express plasmid (kindly gifted from Paul R. Langford at Imperial College London) through enzymatic digestion to construct the recombinant plasmid pSS911, and the resulting plasmid was introduced into the Δfur P. multocida strain to generate complementary strains by electroporation.
Phenotype Characterization
The cell growth curve for P. multocida 0818 was analyzed after recording the OD600 value every 2 h over a 24 h period. The OM proteins of P. multocida and S. Choleraesuis were prepared as previously described (20, 21). Briefly, BHI growth medium was used to grow 200 mL of cells that were harvested after reaching an OD600 of 0.8. Following repeated freeze-thaw cycles, the cells were then centrifuged and washed once in HEPES buffer (10 mM, pH 7.4) and re-suspended in 20 ml of 4-(2-Hydroxyethyl) piperazine-1-ethanesulfonic acid buffer (HEPES buffer) (10 mM, pH 7.4) on ice. Subsequently, the cells were disrupted through sonication (six bursts, 10 s each), and the intact cells were removed through centrifugation (15,600 × g, 5 min, 4°C) (21). The supernatants containing OMPs were transferred to new tubes and centrifuged again (15,600 × g, 30 min, 4°C). The membrane pellets were re-suspended in 2 ml of HEPES buffer. To solubilize the cytoplasmic membrane, 2 ml of 2% sarcosyl was added followed by incubation at room temperature (RT) for 30 min with constant shaking. After centrifugation (15,600 × g, 30 min, 4°C), the pellets containing the OMPs were washed once with 0.5 ml of HEPES buffer and re-suspended in 50 μl of HEPES buffer. The OMP concentrations were measured using the BCA Protein Assay Kit (Thermo Scientific, Rockford, IL, USA). Protein samples were then supplemented with a sample buffer [50 mM Tris, 20% glycerol, 4% SDS, 0.005% bromophenol blue, and 5% β-mercaptoethanol] and boiled at 95°C for 5 min. Subsequently, the samples were assessed by 12.5% SDS-PAGE and stained with Coomassie Brilliant Blue R-250 (Sigma-Aldrich).
Polymyxin B and Deoxycholic Acid (DOC) Resistance Assays
Polymyxin B and DOC resistance assays were performed as previously described (30) with slight modifications. Cells were harvested from BHI broth and washed in PBS when P. multocida growth reached an OD600 of 0.8. The cells were subsequently suspended in BHI medium. The bacterial suspension was diluted to 1 × 106 colony-forming units (CFU) per 100 μl and mixed with polymyxin B (0.5 μg/ml), DOC (250 μg/ml) solution, or untreated, followed by an incubation at 37°C for 1 h. Subsequently, the bacteria were 10-fold serially diluted with BHI medium, and diluted suspensions were spread onto BHI agar plates to determine the number of CFU per 100 μl. The survival rate was calculated as mean CFU of polymyxin B-treated or DOC-treated culture/mean CFU of an untreated culture. All tests were performed in triplicate.
Determination of LD50 in Ducks
All animal experiments in this study were performed in strict accordance with the recommendations delineated in the Guide for the Care and Use of Laboratory Animals of the Ministry of Science and Technology of China. All animal procedures were approved by the Animal Ethics Committee of Sichuan Agricultural University (approval no. 2015-015). Ducklings used in the animal experiments were negative for P. multocida based on a serological assay. One-day-old Sheldrake ducks were purchased from a hatchery and acclimated for 7 days prior to the experiment. The P. multocida 0818 wild-type and Δfur strains were grown in BHI media, then harvested and re-suspended in PBS when the cultures reached an OD600 of 0.8–0.9. The suspensions were subsequently diluted into various concentration doses (1 × 105 CFU to 1 × 108 CFU per 100 μl of PBS). Groups of 1-week-old ducks (10–15 ducks per group) were orally infected with 100 μl of PBS containing different doses of P. multocida 0818 or S985 (Δfur), followed by observation for 2 weeks; the deaths were recorded daily.
Colonization of Sheldrake Tissues by P. multocida
Colonization of the liver, spleen, and lungs by P. multocida in ducks was evaluated as previously described (23). Two groups of 1-week-old Sheldrake ducks (5 ducks per group) were orally infected with 100 μl of PBS containing ~1 × 108 CFU of the P. multocida parent strain or the Δfur strain. All animals were euthanized 3 days later, and 0.1-g samples from each of the selected organs of the infected ducks were removed by dissection. The dissected organ samples were placed in sterile micro-sample bags containing 900 μl of BSG and ground into uniform particles. Subsequently, the homogenates were 10-fold serially diluted, and 100 μl of the diluted suspensions were plated onto agar plates containing BHI or BHI supplemented with Kan.
Serum Antibacterial Test
Duck sera were prepared from healthy adult ducks. The levels of serum antibacterial activity against the P. multocida strains were determined as previously described (31). P. multocida was grown in BHI media to an OD600 of 0.8, harvested and re-suspended in PBS, and the bacterial suspension was diluted to a concentration of ~1 × 104 CFU/ml. An aliquot of 100 μl containing 1 × 104 CFU/ml bacterial suspension was added to 900 μl of duck serum or heat-inactivated duck serum (heated at 56°C for 30 min). One hundred-microliter serum/bacterial mixtures were subsequently 10-fold serially diluted with PBS, and 100 μl of the diluted suspensions were spread onto BHI agar plates (supplemented with or without Kan). The remaining mixtures were cultured at 37°C for 3 h with shaking. Subsequently, the samples were diluted and spread using the same method. Viable counts were determined at 0 and 3 h. The growth rate was calculated as the average CFU/ml at 3 h divided by the average CFU/ml at 0 h. All tests were performed in triplicate.
Immunization and Challenge of Sheldrake Ducks
Sichuan Sheldrake ducks were challenged using a previously described method (21). One-day-old Sheldrake ducks were purchased from a hatchery and acclimated for 7 days prior to the experiment. The ducks were fasted for 6 h prior to oral inoculation with the P. multocida strains and were not fed until 30 min after inoculation. Groups of 65 1-week-old ducks were orally immunized with 100 μl of PBS containing 1 × 105 CFUs of the Δfur strain on day 8 and boosted with the same dose of Δfur on day 18. One set of controls was used: 30 1-week-old ducks that were orally inoculated with 100 μl of PBS. Those ducks that survived 28 days after oral inoculation were challenged with 1 × 108 CFUs of the P. multocida 0818 wild-type strain (500 times 50% lethal doses [LD50]). The ducks were fed twice daily with duckling feed without antibiotics and followed by observation, and deaths were recorded daily. Moribund ducks were euthanized and subsequently necropsied to evaluate the presence of P. multocida in the liver, spleen, and lungs.
Enzyme-Linked Immunosorbent Assay (ELISA)
Blood and bile samples were collected from each group on day 0 (prior to immunization), day 10 (before booster immunization), and day 20 (prior to challenge) as previously described (32). Five randomly selected ducks in the oral immunization group and control group were euthanized separately and then subjected to individual sampling. The P. multocida 0818 wild-type strain was harvested when the culture reached an OD600 of 0.8 in BHI media. A total of 1 × 1010 CFUs of P. multocida was obtained by centrifugation, washed, and suspended in PBS. Then, 1 × 1010 CFUs of whole-cell antigens of P. multocida 0818 were inactivated after heating at 80°C for 10 min. Subsequently, heat-inactivated whole-cell antigens of P. multocida or 0.25 μg/ml purified P. multocida OMP were independently applied to each well of a 96-well ELISA microtitre plate (Nunc-immuno MaxiSorp plate, Nunc, Roskilde, Denmark), both in a volume of 100 μl per well. The plates were then incubated at 4°C overnight, washed three times with 200 μl of PBST per well, and blocked with 2% BSA diluted in PBS (1 h at room temperature). The serum and bile samples from immunized and non-immunized ducks were then diluted to 1:160 and 1:40, respectively, and 100-μl aliquots of these dilutions were added to each well. After incubation for 1 h at 37°C, the samples were washed five times with PBST, and 100 μl of 1:5000-diluted horseradish peroxidase (HRP)-labeled goat anti-duck IgY or IgA was applied to the wells. The plates were incubated for 1 h at 37°C and washed five times with PBST.
Statistical Analysis
All statistics were performed using GraphPad Prism 5 software package (Graph Software, San Diego, CA). Numerical data were expressed as means ± SD and evaluated with a student t-test. Differences were considered significant at P < 0.05. The median lethal dose (LD50) was estimated using a profit analysis.
Results
Heterologous Expression of P. multocida fur Gene in S. Choleraesuis
Based on our bioinformatics analysis, the fur gene of P. multocida Pm70 was 441 bp in length and shared no highly similar sequences with either S. Typhimurium fur (Gene ID 1252213) or E. coli fur (Gene ID 945295); however, the 441-bp sequence of the Pm70 fur gene was predicted to encode a 146-amino acid polypeptide, which showed 67% identity to both S. Typhimurium Fur (Protein ID, AAL19637.1) and E. coli Fur (Protein ID, CDZ19538.1). Thus, we proposed that this protein might perform the same functions as other Fur proteins. To further evaluate the functionality of P. multocida Fur, Δfur mutants of S. Choleraesuis were complemented with the recombinant plasmid pSS908 carrying the P. multocida fur gene in the low-copy-number plasmid pSS664. S. Choleraesuis Δfur (S984) was utilized for complementation assays (Table 1). The Δfur mutant of S. Choleraesuis showed the synthesis and secretion of siderophores. The P. multocida fur gene complemented S. Choleraesuis Δfur, repressing the synthesis and secretion of siderophores in an iron-dependent manner. Moreover, the OMP profiles of S. Choleraesuis (S340), S. Choleraesuis Δfur (S984), and S. Choleraesuis Δfur expressing P. multocida 0818 fur from the plasmid pSS908 were shown. As shown in Figure 1A, the S. Choleraesuis fur mutant (S984) harboring pSS664 displayed a yellow-orange colony, whereas the S. Choleraesuis parent strain (S340) harboring pSS664 and S. Choleraesuis Δfur expressing P. multocida 0818 fur from the plasmid pSS908 displayed no changes in appearance. S. Choleraesuis Δfur with the P. multocida fur gene and the S. Choleraesuis parent strain expressing vector pSS664 exhibited a minor difference of protein bands above 70 kDa compared with S. Choleraesuis Δfur with vector pSS664 (Figure 1B).
Figure 1. Identification of the fur gene of P. multocida. (A) Detection of the secreted siderophores. Three S. Choleraesuis strains (S340 harboring the control plasmid pSS664, S984 (Δfur) harboring pSS664, and S984 harboring pSS908) were cultured on CAS plates, and the colony color was observed. (B) OMP profiles of S. Choleraesuis. OMPs were extracted from S340 (pSS664), S984 (Δfur) (pSS664), and S984 (pSS908) and subsequently subjected to SDS-PAGE. Coomassie blue staining was applied to visualize the OMPs. M refers to the protein marker.
Construction and Characterization of fur Mutants
Construction of the P. multocida Δfur mutant was performed using the suicide T-vector pYA4278 through suicide plasmid-mediated homologous recombination. The genotype was confirmed via PCR using three pairs of primers: 1&2, 3&4, and 5&6 (Figure 2A and Table 2 for primer sequences). Primers 3&6 were both designed based on the sequences of the P. multocida genome outside of the homology arms of fur. A fragment of fur (PM0352) was amplified using primers 1&2 (441 bp). Primers 4&5 were designed to target the middle of the KanR cassette. The two DNA segments amplified using primer pairs 3&4 and 5&6 (~650 bp) were present in the Δfur mutant but not in the parent strain (P. multocida 0818). In contrast, the complete fur sequence (1&2) was present in the parent strain but not in the Δfur mutant. The positive control 16S RNA sequence (501 bp) was amplified from both strains (Figure 2B). To further characterize the Δfur mutant, the Fur protein of P. multocida was synthesized and purified, and a polyclonal mouse anti-Fur antibody was prepared for western blot assays. Fur was detected in the parent strain and the complementary S985 (pSS911) strain but not in S985 (Figure 2C), indicating that the fur gene was successfully deleted from S985.
Figure 2. Construction of the Δfur mutant in P. multocida strain 0818. (A) Schematic strategy for deletion of the target fur gene (PM0352). The targeted fur gene of P. multocida 0818 was replaced with a KanR cassette through homologous recombination. Primers 1, 2, 3, 4, 5, and 6 were designed for the selection and characterization of the mutant clones. (B) Identification of the constructed Δfur mutant. The parent strain and S985 (Δfur) were identified using three primer pairs: 1&2, 3&4, and 5&6. M refers to the DNA marker; 16S RNA denotes amplification of the positive control gene in both strains. (C) Detection of Fur expression in the P. multocida strains. The parent strain, S985 (Δfur), and S985 (pSS911) were harvested at an OD600 of 1.0 following growth in BHI media, and the expression of Fur was subsequently detected in these strains using the mouse anti-Fur antibody via western blotting.
Phenotype Characterization of the P. multocida fur Mutant
To evaluate the influence of fur deletion in P. multocida, the growth rates, OMP profiles, and serum complement sensitivity of the parent strain (P. multocida 0818), S985 (P. multocida Δfur), and S985 (pSS911) were determined. The parent strain, S985 (P. multocida Δfur), and S985 (pSS911) exhibited similar and typical growth and proliferation patterns in a normal nutrient-rich (fresh BHI medium) environment (0-8 h), consisting of a (plateau) lag phase (0–2 h), a logarithmic phase (2–6 h), and a stationary phase (6–8 h; Figure 3A). However, when the iron ions in the medium were depleted, the S985 (P. multocida Δfur) curve displayed a decline of growth due to the lack of effective iron control (12–14 h; Figure 3B). Both the parent strain and S985 (pSS911) exhibited similar, typical growth rates in both untreated and heat-treated healthy duck serum; however, S985 (P. multocida Δfur) was rapidly killed in untreated serum but grew in heat-treated serum. We also observed that S985 did not grow to the extent of the parent and complemented strains in heat-treated serum (Table 3).
Figure 3. Phenotypic characterization of the P. multocida Δfur mutant. (A,B) Growth curves of the parent strain and the Δfur mutant. P. multocida was cultured in BHI broth supplemented with or without chloramphenicol, and the OD600 value of each strain was subsequently measured every 2 h over a period of 20 h, no significant differences for bacterial growth curves were observed. (C,D) Resistance to polymyxin B (C) or DOC (D). The parent strain, S985 (Δfur) and the complement strain S985 (pSS911) were cultured in BHI broth with or without 0.5 μg/ml polymyxin B and 250 μg/ml DOC for 1 h, and the survival ratio of each strain was calculated as the mean CFUs of the treated group divided by the mean CFUs of the untreated group. The data in (C,D) are expressed as the means ± SD analyzed at significance levels of 0.05 (*) and 0.001 (***).
Table 3. Serum bactericidal assay.
Alterations in the cell membrane protein profiles might reduce bacterial resistance to bile salts and cationic antimicrobial peptides; therefore, we examined the susceptibility of S985 (P. multocida Δfur) to polymyxin B and DOC. The Δfur mutant was highly sensitive to polymyxin B and DOC. After treatment with polymyxin B or DOC, the survival ratio of strain S985 (P. multocida Δfur) was significantly lower than that of both the parent strain and the complement strain, S985 (pSS911) (Figures 3C,D).
Wild-type P. multocida, grown under iron-limited conditions, and the P. multocida S985 fur mutant strain, grown under normal or iron-deficient conditions, both showed increased abundance of OMPs over 70 kDa (Figure 4), suggesting that these proteins are iron-Fur dependent and regulated by Fur or are associated with iron uptake.
Figure 4. Fur-regulated IROMPs in P. multocida. Outer membrane profiles of P. multocida. P. multocida strain grown under iron-rich (+Fe, BHI broth) and iron-limited (-Fe, BHI broth plus 100 μM 2,2′-dipyridyl, DIP) conditions.
Colonization of the P. multocida Parent Strain and the Δfur Mutant Strain
We evaluated the colonization of the liver, spleen, and lungs of ducks 2 days after the oral inoculation with P. multocida Δfur compared with the wild-type. We observed that the colonization levels of P. multocida Δfur in the liver, spleen, and lungs were all lower than those of the wild-type (Figure 5). Interestingly, for the parent strain, the level of colonization in the lungs was higher than those in the liver and spleen; however, for the fur mutant strain, the bacterial load in the lungs was lower compared with those in the liver and spleen.
Figure 5. Colonization of the P. multocida parent strain and the Δfur mutant strain. Colonization of Sheldrake duck tissues (n = 5) by the P. multocida parent strain and S985 (Δfur) after 2 days of oral inoculation with 108 CFUs of the respective strain. Each point represents one duck. The data are expressed as the means ± SD analyzed with significance levels set to 0.01 (**) and 0.001 (***).
Virulence of the P. multocida Δfur Mutant in the Duck Host
To determine the effects of the P. multocida Δfur mutant on virulence, the LD50 of the parent strain and fur mutant strain were evaluated in Sheldrake ducks, the natural host of P. multocida 0818. P. multocida Δfur administered orally was attenuated. The LD50 of the Δfur mutant was 1.65 × 108 (Table 4), which represented an approximate 146-fold increase in the LD50 compared to that of the wild-type strain, which had an LD50 of 1.13 × 106.
Table 4. Determination of the LD50 of P. multocida 0818 and the Δfur mutant.
Immune Protection of P. multocida Δfur Mutants in the Duck Host
Ducks orally immunized with P. multocida Δfur survived an oral challenge with 500 × the LD50 of wild-type P. multocida 0818 for 10 days after the second immunization, whereas all control ducks died within 2 days (Table 5). The whole-bacterial antigen- and OMP-specific responses of the serum IgY and bile IgA levels of ducks in samples prepared on day 0 prior to immunization and days 10 and 20 post-immunization were measured by indirect ELISA. As shown in Figure 6, no specific serum IgY or bile IgA responses to P. multocida antigens or OMPs were present 1 day prior to immunization. In contrast, the immunized group exhibited strong stimulation of serum IgY and bile IgA antibodies against whole-bacterial antigens and demonstrated significant IgY and IgA titres in the blood and bile, respectively. Similar to the results for whole-bacterial antigens, OMP antigen-specific serum IgY and bile IgA responses were stronger in the vaccinated group than in the control group.
Table 5. Protection rate conferred by the Δfur mutant.
Figure 6. Antibody responses induced by the Δfur mutant in ducks. Ducks (n = 5) were orally immunized with the Δfur mutant strain twice at an interval of 10 days. Subsequently, the responses of serum IgY specific to whole-bacterial antigens (A) and OMPs (C) and the responses of bile IgA specific to whole-bacterial antigen (B) and OMPs (D) were measured using indirect ELISA. The data are expressed as the means ± SD analyzed with significance levels set to 0.01 (**) and 0.001 (***).
Discussion
P. multocida causes acute infectious zoonotic diseases in humans and animals, including fowl cholera in birds, atrophic rhinitis in pigs, haemorrhagic septicaemia in cattle, snuffles in rabbits, and wound abscesses or meningitis in humans (33). These diseases pose significant threats to the food industry. Global regulators modulate the transcription and expression of hundreds of genes, including virulence factors. Therefore, studies of P. multocida global regulators provide a better understanding of immune escape mechanisms in hosts and facilitate the identification of toxins or released toxic by-products as well as other molecular mechanisms underlying host-pathogen interactions, providing significant aid to vaccine development. In-depth analyses performed in recent years have revealed several global regulators in P. multocida, such as PhoP/PhoQ (21), Crp (cAMP receptor protein) (28), Fis (nucleoid-associated protein) (34), and FnrP (35, 36). Fur is another global regulator that regulates the expression of iron acquisition genes and many other virulence genes.
Fur acts as either a repressor or an activator to regulate gene expression by binding as a Fe2+-Fur complex to promoter regions (Fur boxes) upstream of its target genes. Fe2+-Fur activation might occur directly, creating a competition for DNA binding sites with other inhibitors (37). For example, the bacterial ferritin gene bfrB in Pseudomonas aeruginosa and the hilD virulence factor in S. Typhimurium are both activated by Fe2+-Fur (38, 39). The indirect activation of Fur might also occur. Fe2+-Fur negatively regulates certain small regulatory RNAs. RyhB, a small non-coding RNA (sRNA), triggers the mRNA degradation through base pairing via RNase E and RNase III in E. coli (40). In P. multocida, a DNA microarray analysis showed that the transcriptional regulation of 174 genes was altered under iron-limited growth conditions, which includes virulence genes involved in iron uptake from the environment (41). In the present study, profiling of P. multocida outer membrane proteins revealed that the OMP bands of parent strains grown in iron-deficient environments were similar to the OMP bands of Δfur mutant strains grown under iron-limited or iron-sufficient conditions (Figure 3), suggesting that the expression of iron-regulated outer membrane proteins (IROMPs) is regulated by Fur and iron-restricted conditions.
As a global regulator, Fur plays essential roles in virulence. The absence or impairment of the fur gene leads to the attenuation of bacterial virulence. For example, Δfur mutation in Campylobacter jejuni has a significant effect on its colonization of the gastrointestinal tract in chickens (42). Fur proteins of N. meningitidis and N. gonorrhoeae are involved in disease pathogenesis(43). In the present study, we observed that the virulence of S985 (Δfur) decreased 146-fold after oral inoculation, indicating that the Fur protein plays a role in the regulation of virulence factors of P. multocida in the host iron-limiting environment. However, the decreased virulence observed with the P. multocida Δfur mutant strain was markedly less pronounced than in other common strains, such as Edwardsiella ictaluri (6) and S. Typhimurium (44). Thus, other genes of P. multocida might replace the fur iron regulation function to minimize adverse effects in fur-deleted strains. Moreover, the virulence of the P. multocida 0818 Δfur mutant strain differed from that of another P. multocida Δfur mutant strain constructed in 2001. The 2001 P. multocida Δfur mutant strain exhibited the same virulence as its P. multocida 0818 parent strain (LD50 = 5 CFU per animal) when the bacteria were intraperitoneally inoculated in Swiss mice (45). Similar results has been also observed for S. Typhimurium Δfur mutant when inoculated intraperitoneally and orally into the mouse model (44), indicating that the infective route, the environmental conditions, and the selected host affect the virulence of the bacteria.
Live oral vaccines encounter gastric acid in the stomach, the weak alkaline environment, bile salt (e.g., DOC), and antimicrobial peptides (e.g., polymyxin B) in the intestinal fluid, en route to the site of host cells targeted for colonization. We observed the increased susceptibility of S985 (Δfur) to polymyxin B and DOC (data not shown), consistent with other observations that the fur mutant strains of S. Typhimurium and the uropathogenic E. coli strain CFT073 show significantly greater sensitivity to oxidative stress compared to their wild-type parents (46, 47). S985 (Δfur) died more rapidly under insufficient nutrient conditions when grown in BHI over 12 h (Figure 3B) and showed a more obvious decrease in the colonization of these organs (Figure 4). Although fur mutant strains exhibited certain defects in host organs, the bacterial load and persistence of S985 (Δfur) in host organs was sufficient to induce humoral immune responses (Figure 4) and provided sufficient protection against a lethal dose challenge of wild-type P. multocida.
In addition to achieving adequate homologous protection, cross-protection has become a criterion for evaluating vaccine quality. Many researchers have suggested that the growth of P. multocida in vivo, particularly in iron-restricted environments, results in iron-regulated OMPs that are adequate cross-protective antigens (48–50). The P. multocida fur mutant continuously synthesized IROMPs, resulting in exposure to more epitopes that continuously stimulate the host's immune system to produce antibodies (Figures 6A–D). However, these high levels of IgG are not fully protective under challenges with high doses of virulent P. multocida. Although the P. multocida Δfur mutant strain is not a highly efficient and safe live-attenuated vaccine, fur is a global regulatory gene that controls the expression of virulence factors. Accordingly, the fur gene might be deleted together with other genes, such as the regulators phoP and crp or fis (essential for capsule production) (34) to increase attenuation and enhance cross-protection effects and immunogenicity of the vaccine.
Author Contributions
QK and QL designed the experiments and analyzed the data. YH and PL carried out the experiments. QL, YH and QK wrote the paper.
Funding
This work was supported by grants 31472179 and 31570928 from the Natural Science Foundation of China and the fundamental research funds for the Central Universities (SWU117061, SWU117062).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We thank Dr. Matthew Bellefleur for his help in language revision and correction of grammatical errors.
References
1. Ebanks RO, Goguen M, Knickle L, Dacanay A, Leslie A, Ross NW, et al. Analysis of a ferric uptake regulator (Fur) knockout mutant in Aeromonas salmonicida subsp salmonicida. Vet Microbiol. (2013) 162:831–41. doi: 10.1016/j.vetmic.2012.10.038
2. Becerra G, Merchan F, Blasco R, Igeno MI. Characterization of a ferric uptake regulator (Fur)-mutant of the cyanotrophic bacterium Pseudomonas pseudoalcaligenes CECT5344. J Biotechnol. (2014) 190:2–10. doi: 10.1016/j.jbiotec.2014.03.030
3. Porcheron G, Dozois CM. Interplay between iron homeostasis and virulence: Fur and RyhB as major regulators of bacterial pathogenicity. Vet Microbiol. (2015) 179:2–14. doi: 10.1016/j.vetmic.2015.03.024
4. Lill R. Function and biogenesis of iron-sulphur proteins. Nature (2009) 460:831–8. doi: 10.1038/nature08301
5. Fillat MF. The FUR (ferric uptake regulator) superfamily: diversity and versatility of key transcriptional regulators. Arch Biochem Biophys. (2014) 546:41–52. doi: 10.1016/j.abb.2014.01.029
6. Santander J, Golden G, Wanda SY, Curtiss R. Fur-regulated iron uptake system of Edwardsiella ictaluri and its influence on pathogenesis and immunogenicity in the catfish host. Infect Immun. (2012) 80:2689–703. doi: 10.1128/IAI.00013-12
7. Coy M, Neilands JB. Structural dynamics and functional domains of the fur protein. Biochemistry (1991) 30:8201–10. doi: 10.1021/bi00247a016
8. Johnson M, Cockayne A, Williams PH, Morrissey JA. Iron-responsive regulation of biofilm formation in Staphylococcus aureus involves fur-dependent and fur-independent mechanisms. J Bacteriol. (2005) 187:8211–5. doi: 10.1128/JB.187.23.8211-8215.2005
9. Sun FJ, Gao H, Zhang YQ, Wang L, Fang N, Tan YF, et al. Fur is a repressor of biofilm formation in Yersinia pestis. PLoS ONE (2012) 7:e52392. doi: 10.1371/journal.pone.0052392
10. Banin E, Vasil ML, Greenberg EP. Iron and Pseudomonas aeruginosa biofilm formation. Proc Natl Acad Sci USA. (2005) 102:11076–81. doi: 10.1073/pnas.0504266102
11. Horsburgh MJ, Ingham E, Foster SJ. in Staphylococcus aureus, fur is an interactive regulator with PerR, contributes to virulence, and is necessary for oxidative stress resistance through positive regulation of catalase and iron homeostasis. J Bacteriol. (2001) 183:468–75. doi: 10.1128/JB.183.2.468-475.2001
12. Mey AR, Wyckoff EE, Kanukurthy V, Fisher CR, Payne SM. Iron and fur regulation in Vibrio cholerae and the role of fur in virulence. Infect Immun. (2005) 73:8167–78. doi: 10.1128/IAI.73.12.8167-8178.2005
13. Rea RB, Gahan CGM, Hill C. Disruption of putative regulatory loci in Listeria monocytogenes demonstrates a significant role for fur and PerR in virulence. Infect Immun. (2004) 72:717–27. doi: 10.1128/IAI.72.2.717-727.2004
14. Torres VJ, Attia AS, Mason WJ, Hood MI, Corbin BD, Beasley FC, et al. Staphylococcus aureus fur regulates the expression of virulence factors that contribute to the pathogenesis of Pneumonia. Infect Immun. (2010) 78:1618–28. doi: 10.1128/IAI.01423-09
15. Harper M, Boyce JD, Adler B. Pasteurella multocida pathogenesis: 125 years after Pasteur. Fems Microbiol Lett. (2006) 265:1–10. doi: 10.1111/j.1574-6968.2006.00442.x
16. Gonzalez CT, Maheswaran SK. The role of induced virulence factors produced by Pasteurella haemolytica in the pathogenesis of bovine pneumonic pasteurellosis: review and hypotheses. Br Vet J. (1993) 149:183–93. doi: 10.1016/S0007-1935(05)80088-0
17. Heddleston KL, Gallagher JE, Rebers PA. Fowl cholera: gel diffusion precipitin test for serotyping Pasteruella multocida from avian species. Avian Dis. (1972) 16:925–36. doi: 10.2307/1588773
18. Sellyei B, Thuma A, Volokhov D, Varga Z. Comparative analysis of Pasteurella multocida isolates from acute and chronic fowl cholera cases in Hungary during the period 2005 through 2010. Avian Dis. (2017) 61:457–65. doi: 10.1637/11674-051817-Reg.1
19. Ahmad TA, Rammah SS, Sheweita SA, Haroun M, El-Sayed LH. Development of immunization trials against Pasteurella multocida. Vaccine (2014) 32:909–17. doi: 10.1016/j.vaccine.2013.11.068
20. Liu Q, Liu Q, Zhao X, Liu T, Yi J, Liang K, et al. Immunogenicity and cross-protective efficacy induced by outer membrane proteins from Salmonella typhimurium mutants with truncated LPS in mice. Int J Mol Sci. (2016) 17:416. doi: 10.3390/ijms17030416
21. Xiao KP, Liu Q, Liu XY, Hu YL, Zhao XX, Kong QK. Identification of the Avian Pasteurella multocida phoP gene and evaluation of the effects of phoP deletion on virulence and immunogenicity. Int J Mol Sci. (2016) 17:12. doi: 10.3390/ijms17010012
22. Roland K, Curtiss RIII, Sizemore D. Construction and evaluation of a delta cya delta crp Salmonella typhimurium strain expressing avian pathogenic Escherichia coli O78 LPS as a vaccine to prevent airsacculitis in chickens. Avian Dis. (1999) 43:429–41. doi: 10.2307/1592640
23. Kong Q, Yang J, Liu Q, Alamuri P, Roland KL, Curtiss RIII. Effect of deletion of genes involved in lipopolysaccharide core and O-antigen synthesis on virulence and immunogenicity of Salmonella enterica serovar typhimurium. Infect Immun. (2011) 79:4227–39. doi: 10.1128/IAI.05398-11
24. Bosse JT, Durham AL, Rycroft AN, Kroll JS, Langford PR. New plasmid tools for genetic analysis of Actinobacillus pleuropneumoniae and other Pasteurellaceae. Appl Environ Microb. (2009) 75:6124–31. doi: 10.1128/AEM.00809-09
25. May BJ, Zhang Q, Li LL, Paustian ML, Whittam TS, Kapur V. Complete genomic sequence of Pasteurella multocida, Pm70. Proc Natl Acad Sci USA. (2001) 98:3460–5. doi: 10.1073/pnas.051634598
26. Schwyn B, Neilands JB. Universal chemical assay for the detection and determination of siderophores. Anal Biochem. (1987) 160:47–56. doi: 10.1016/0003-2697(87)90612-9
27. Curtiss RIII, Wanda SY, Gunn BM, Zhang X, Tinge SA, Ananthnarayan V, et al. Salmonella enterica serovar typhimurium strains with regulated delayed attenuation in vivo. Infect Immun. (2009) 77:1071–82. doi: 10.1128/IAI.00693-08
28. Zhao XX, Liu Q, Xiao KP, Hu YL, Liu XY, Li YY, et al. Identification of the crp gene in avian Pasteurella multocida and evaluation of the effects of crp deletion on its phenotype, virulence and immunogenicity. BMC Microbiol. (2016) 16:125. doi: 10.1186/s12866-016-0739-y
29. Edwards RA, Keller LH, Schifferli DM. Improved allelic exchange vectors and their use to analyze 987P fimbria gene expression. Gene (1998) 207:149–57. doi: 10.1016/S0378-1119(97)00619-7
30. Huang C, Liu Q, Luo Y, Li P, Liu Q, Kong Q. Regulated delayed synthesis of lipopolysaccharide and enterobacterial common antigen of Salmonella Typhimurium enhances immunogenicity and cross-protective efficacy against heterologous Salmonella challenge. Vaccine (2016) 34:4285–92. doi: 10.1016/j.vaccine.2016.07.010
31. Chung JY, Wilkie I, Boyce JD, Townsend KM, Frost AJ, Ghoddusi M, et al. Role of capsule in the pathogenesis of fowl cholera caused by Pasteurella multocida serogroup A. Infect Immun. (2001) 69:2487–92. doi: 10.1128/IAI.69.4.2487-2492.2001
32. Wilkie IW, Grimes SE, O'Boyle D, Frost AJ. The virulence and protective efficacy for chickens of Pasteurella multocida administered by different routes. Vet Microbiol. (2000) 72:57–68. doi: 10.1016/S0378-1135(99)00187-X
33. Boyce JD, Adler B. How does Pasteurella multocida respond to the host environment? Curr Opin Microbiol. (2006) 9:117–22. doi: 10.1016/j.mib.2005.12.005
34. Steen JA, Steen JA, Harrison P, Seemann T, Wilkie I, Harper M, et al. Fis is essential for capsule production in Pasteurella multocida and regulates expression of other important virulence factors. PLoS Pathog. (2010) 6:e1000750. doi: 10.1371/journal.ppat.1000750
35. Uhlich GA, McNamara PJ, Iandolo JJ, Mosier DA. Cloning and characterization of the gene encoding Pasteurella haemolytica FnrP, a regulator of the Escherichia coli silent hemolysin sheA. J Bacteriol. (1999) 181:3845–8.
36. Uhlich GA, McNamara PJ, Iandolo JJ, Mosier DA. FnrP interactions with the Pasteurella haemolytica leukotoxin promoter. Fems Microbiol Lett. (2000) 186:73–7. doi: 10.1111/j.1574-6968.2000.tb09084.x
37. Isabella V, Wright LF, Barth K, Spence JM, Grogan S, Genco CA, et al. cis- and trans-acting elements involved in regulation of norB (norZ), the gene encoding nitric oxide reductase in Neisseria gonorrhoeae. Microbiology (2008) 154:226–39. doi: 10.1099/mic.0.2007/010470-0
38. Wilderman PJ, Sowa NA, FitzGerald DJ, FitzGerald PC, Gottesman S, Ochsner UA, et al. Identification of tandem duplicate regulatory small RNAs in Pseudomonas aeruginosa involved in iron homeostasis. Proc Natl Acad Sci USA. (2004) 101:9792–7. doi: 10.1073/pnas.0403423101
39. Teixido L, Carrasco B, Alonso JC, Barbe J, Campoy S. Fur activates the expression of Salmonella enterica pathogenicity island 1 by directly interacting with the hilD operator in vivo and in vitro. PLoS ONE (2011) 6:e19711. doi: 10.1371/journal.pone.0019711
40. Masse E, Escorcia FE, Gottesman S. Coupled degradation of a small regulatory RNA and its mRNA targets in Escherichia coli. Genes Dev. (2003) 17:2374–83. doi: 10.1101/gad.1127103
41. Paustian ML, May BJ, Kapur V. Pasteurella multocida gene expression in response to iron limitation. Infect Immun. (2001) 69:4109–15. doi: 10.1128/IAI.69.6.4109-4115.2001
42. Palyada K, Threadgill D, Stintzi A. Iron acquisition and regulation in Campylobacter jejuni. J Bacteriol. (2004) 186:4714–29. doi: 10.1128/JB.186.14.4714-4729.2004
43. Yu C, Genco CA. Fur-mediated global regulatory circuits in pathogenic Neisseria species. J Bacteriol. (2012) 194:6372–81. doi: 10.1128/JB.00262-12
44. Garcia-del Portillo F, Foster JW, Finlay BB. Role of acid tolerance response genes in Salmonella typhimurium virulence. Infect Immun. (1993) 61:4489–92.
45. Bosch M, Tarrago R, Garrido ME, Campoy S, de Henestrosa ARF, de Rozas AMP, et al. Expression of the Pasteurella multocida ompH gene is negatively regulated by the Fur protein. Fems Microbiol Lett. (2001) 203:35–40. doi: 10.1111/j.1574-6968.2001.tb10817.x
46. Kim JN, Kwon YM. Genetic and phenotypic characterization of the RyhB regulon in Salmonella Typhimurium. Microbiol Res. (2013) 168:41–9. doi: 10.1016/j.micres.2012.06.007
47. Porcheron G, Habib R, Houle S, Caza M, Lepine F, Daigle F, et al. The small RNA RyhB contributes to siderophore production and virulence of uropathogenic Escherichia coli. Infect Immun. (2014) 82:5056–68. doi: 10.1128/IAI.02287-14
48. Rimler RB. Cross-protection factor(s) of Pasteurella multocida: passive immunization of turkeys against fowl cholera caused by different serotypes. Avian Dis. (1987) 31:884–7. doi: 10.2307/1591046
49. Rimler RB, Rhoades KR. Solubilization of membrane-associated cross-protection factor(s) of Pasteurella multocida. Avian Dis. (1989) 33:258–63. doi: 10.2307/1590841
Keywords: Pasteurella multocida, fur, iron uptake, attenuated live vaccine, immune responses
Citation: Liu Q, Hu Y, Li P and Kong Q (2019) Identification of Fur in Pasteurella multocida and the Potential of Its Mutant as an Attenuated Live Vaccine. Front. Vet. Sci. 6:5. doi: 10.3389/fvets.2019.00005
Received: 14 September 2018; Accepted: 10 January 2019;
Published: 04 February 2019.
Edited by:
Subhash Verma, Chaudhary Sarwan Kumar Himachal Pradesh Krishi Vishvavidyalaya, IndiaReviewed by:
Filip Boyen, Ghent University, BelgiumLisa Bielke, The Ohio State University, United States
Copyright © 2019 Liu, Hu, Li and Kong. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Qingke Kong, a29uZ3Fpa2lAMTYzLmNvbQ==
†These authors have contributed equally to this work