- Research article
- Open Access
The Riemerella anatipestifer M949_RS01035 gene is involved in bacterial lipopolysaccharide biosynthesis
Veterinary Research volume 49, Article number: 93 (2018)
In this study, the Riemerella anatipestifer mutant strain RA1062 was obtained by screening a random Tn4351 transposon mutant library. The mutant strain was unreactive with the anti-CH3 lipopolysaccharide monoclonal antibody, as demonstrated with an enzyme-linked immunosorbent assay, and its M949_RS01035 gene was inactivated. When cultured in trypticase soy broth, the late stage growth of the mutant RA1062 was significantly decreased. The mutant RA1062 was stained with crystal violet and presented a rough lipopolysaccharide phenotype, which differed from that of the wild-type strain CH3, suggesting that deletion of the M949_RS01035 gene resulted in defective lipopolysaccharide. Silver staining and Western blot analyses further confirmed that the RA1062 lipopolysaccharide had a deficiency in ladder-like binding pattern, as compared to lipopolysaccharide of the wild-type CH3 strain. In addition, the mutant RA1062 showed a higher susceptibility to complement-dependent killing, increased bacterial adhesion and invasion capacities to Vero cells, decreased blood bacterial loads, and attenuated virulence in infected ducks, when compared to the wild-type strain CH3. Moreover, RNA-Seq and real-time polymerase chain reaction analyses indicated that two genes were up-regulated and two were down-regulated in the mutant RA1062 genome. Furthermore, an animal protection experiment showed that immunization of ducks with inactivated RA1062 bacterin conferred effective cross-protection against challenge with the virulent R. anatipestifer serotypes 1, 2, and 10. This study presents evidence that the M949_RS01035 gene is involved in bacterial phenotype, virulence, and gene regulation in R. anatipestifer. The mutant strain RA1062 could be used as a cross-protective vaccine candidate.
Riemerella anatipestifer is a Gram-negative, non-motile, non-spore forming, and atrichous bacterium belonging to the family Flavobacteriaceae in the phylum Bacteroidetes . R. anatipestifer has received considerable attention because infection of this bacterium cause major economic losses to the duck farming industry through high mortality, weight loss, and high treatment costs . R. anatipestifer infection is a highly contagious disease that causes fibrinous pericarditis, airsacculitis, and perihepatitis in ducks . To date, 21 R. anatipestifer serotypes have been identified worldwide; however, there is currently no effective cross-protection between different serotypes [4, 5]. In China, serotypes 1, 2, and 10 are responsible for most outbreaks of R. anatipestifer infection .
A variety of vaccines to protect farm ducks from R. anatipestifer infection have been investigated, including the chaperonin GroEL and inactivated bacterin [7, 8]. On the other hand, some potential virulence factors of R. anatipestifer have been described, including OmpA, VapD, and CAMP cohemolysin [9, 10].
As with other Gram-negative bacteria, lipopolysaccharide (LPS) is probably one of the most important virulence factors of R. anatipestifer. LPS is the major outer membrane component of Gram-negative bacteria and typically comprise three structure domains: lipid A, a core oligosaccharide, and a polysaccharide O-antigen. The O-antigen is composed of repeated sequences of three to six sugar residues and is also the most structurally variable cell surface constituent. Pathogenic bacterial O-antigen plays important roles in avoiding phagocytosis and resisting the lytic action of the complement system [11,12,13,14]. The O-antigen has been also a basis for vaccine development against a variety of human pathogens, . In several species of bacteria, such as Escherichia coli and Salmonella enterica, the genes involved in LPS biogenesis have been well characterized [16, 17]. In R. anatipestifer, five genes (AS87_04050, M949_1556, M949_1603, M949_1360 and M949_RS01915) associated with LPS synthesis have been characterized in our previous studies [18,19,20,21,22].
In this report, we described the identification of the R. anatipestifer mutant strain RA1062, in which the M949_RS01035 gene was disrupted by insertion of the Tn4351 transposon. Furthermore, the LPS phenotype, bacterial virulence, gene regulation, and the cross-protection of the mutant strain RA1062 were characterized.
Materials and methods
The study protocol was approved by the Institutional Animal Care and Use Committee of Shanghai Veterinary Research Institute, the Chinese Academy of Agricultural Sciences (Approval No. Shvri-po-0072), and conducted in strict accordance with the recommendations outlined in the Guide for the Care and Use of Laboratory Animals. One-day old Cherry Valley ducks were obtained from Zhuang Hang Duck Farm (Shanghai, China) and housed in cages under a controlled temperature of 28–30 °C with water and food ad libitum under biosafety conditions.
Bacterial strains and growth conditions
The bacterial strains and plasmids used in this study are listed in Table 3. The mutant strain RA1062 was derived from the wild-type (WT) R. anatipestifer strain CH3 (GenBank accession no. CP006649), which is a serotype 1 strain, by insertion of the Tn4351 plasmid. R. anatipestifer strains were grown on tryptic soy agar (TSA; Difco Laboratories, Franklin Lakes, NJ, USA) or in tryptic soy broth (TSB, Difco) at 37 °C under an atmosphere of 5% CO2.
Escherichia coli strains were grown at 37 °C on Luria–Bertani (LB) plates or in LB broth. When necessary, antibiotics were added to the medium at the following concentrations: kanamycin at 50 μg/mL and erythromycin at 0.5 μg/mL for the mutant strain RA1062, and chloramphenicol at 30 μg/mL for E. coli strain BW19851 (pEP4351).
Indirect enzyme-linked immunosorbent assay (ELISA)
An indirect ELISA was used to screen the Tn4351 insertion mutant library for strains with defective reactivity to the anti-CH3 LPS monoclonal antibody (mAb), as described previously . Briefly, each well of a 96-well ELISA plate was coated with whole mutant cells at 109 CFU/well in 50 μL of carbonate-buffered saline (pH 9.6). The plates were heat-dried overnight at 55 °C. After washing three times with phosphate-buffered saline (PBS) containing 0.05% Tween 20 (PBST), the plates were blocked with PBS-5% skim milk at 37 °C for 2 h, then incubated with anti-CH3 LPS mAb, as the primary antibody, and horseradish peroxidase-conjugated anti-mouse immunoglobulin G (IgG; Tiangen Biotech (Beijing) Co., Ltd., Beijing, China), as the secondary antibody. The reaction was visualized by adding 100 μL of 3,30,5,50-tetramethyl benzidine (Tiangen) and stopped by the addition of 50 μL of 2 M H2SO4. The resulting optical density at 450 nm (OD450) was measured using a plate reader (Synergy 2; BioTek Instruments, Inc., Winooski, VT, USA). The WT strain CH3 was used as a positive coating control. The mutants with its coating wells presenting an OD450 reading more than 2.1 times lower than that of the positive coating control were collected for further characterization. All mutants were screened in triplicate.
Characterization of the mutant strain RA1062
The WT strain CH3 and the mutant strain RA1062 were identified by polymerase chain reaction (PCR) analysis with the primers 16S rRNA F/16S rRNA R, Erm-F/Erm-R, and RA1062-F/RA1062-R (Table 1). Bacterial mutants were constructed by transposon mutagenesis as described previously , with modifications. The E. coli BW19851 with the plasmid pEP4351 was used as the donor strain and R. anatipestifer CH3 as the recipient. For bacterial mating, both donor and recipient bacteria were grown to mid-logarithmic phase, mixed at a bacterial colony forming units (CFU) ratio of 1:2 and centrifuged at 5500 × g for 10 min. The bacterial pellet was washed and re-suspended with 10 mM MgSO4, and filtered through a Millipore membrane, which was then placed face up on TSA with 1 μg/mL erythromycin and 50 μg/mL kanamycin. Following overnight incubation at 30 °C, the bacteria were scraped off the filter, resuspended in 5 mL 10 mM MgSO4, and spread on TSA containing erythromycin and kanamycin to select for transconjugants. To confirm one insertion of the Tn4351 transposon in the mutant strain RA1062, Southern blot analysis was conducted as previously described . Briefly, genomic DNA of the mutant strain RA1062 was extracted using the TIANamp Bacteria DNA kit (Tiangen), digested with the endonuclease XbaΙ, then subjected to sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) and transferred to a nitrocellulose membrane. After washing with saline sodium citrate, the membrane was immobilized for 2 h at 80 °C. PCR was performed to amplify the transposon-specific probe representing the 410-bp IS4351 fragment from the plasmid pEP4351 using the primer pair Tn4351-F/Tn4351-R (Table 1). The probe was generated and hybridization was conducted using the DIG DNA labeling and detection kit (Roche Diagnostics USA, Indianapolis, IN, USA), according to the manufacturer’s protocol. The plasmid pEP4351 and genomic DNA of the WT strain CH3 were also subjected to hybridization analysis, for use as the positive and negative controls, respectively.
The nucleotide sequence surrounding the transposon insertion site was determined using inverse PCR, as described elsewhere . Briefly, genomic DNA was digested with the restriction enzyme HindIII and then treated with T4 ligase, which resulted in the formation of circular molecules. Primer pairs specific for Tn4351 (primers TN-1 and IS4351-F) were used to amplify the sequences adjacent to the insertion site using the LA PCR kit (TaKaRa Biotechnology (Dalian) Co., Ltd., Dalian, China). The nucleotide sequence was compared to sequence in the National Center for Biotechnology Information database using the BLASTX program .
Determination of the bacterial growth curves
The growth rates of the WT strain CH3 and the mutant strain RA1062 were determined and compared. Briefly, each bacterial strain was cultivated in TSB at 37 °C for 8 h with shaking. Equal amounts of each bacterial culture were then inoculated into 12-mL of fresh TSB medium at a ratio of 1:100 (v/v) and incubated at 37 °C, with shaking at 200 rpm, respectively. The bacterial growth rate was measured by counting the bacterial CFU at 2-h intervals for 16 h.
Characterization of the bacterial phenotype
The phenotypes of the WT strain CH3 and its mutant strain RA1062 were determined by crystal violet staining, as previously described  with slight modifications. The bacterial cells were grown in TSB at 37 °C for 8 h, then harvested by centrifugation at 5000 rpm for 5 min, suspended in sterile PBS to a density of 103 cells/mL, plated (0.1 mL) onto TSA, and incubated at 37 °C under an atmosphere of 5% CO2 for 36 h. The plates were then gently flooded with 5 mL of 0.5% crystal violet solution. After staining for 1 min, the excess stain was removed and the plates were examined immediately.
LPS extraction, silver staining, and Western blot analysis
Lipopolysaccharide was extracted from WT strain CH3 and mutant strain RA1062 cells using an LPS extraction Kit (iNtRON Biotechnology, Boca Raton, FL, USA), according to the manufacturer’s instruction. Purified LPS was analyzed by SDS-PAGE. Gels were stained with silver to visualize the presence of LPS  and then counter stained with Coomassie blue to exclude contaminating proteins.
For Western blot analysis, the purified LPS was separated by SDS-PAGE and then transferred onto nitrocellulose membranes (Millipore, Billerica, MA, USA), which were blocked overnight at 4 °C in PBS containing 5% nonfat milk and then washed with PBST. The blots were incubated with anti-CH3 LPS as the primary mAb and then with an IRDYE680CW-conjugated donkey anti-mouse IgG polyclonal antibody (LI-COR Biosciences, Lincoln, NE, USA) for 1 h. The blots were visualized with an Odyssey two-color infrared imaging system (LI-COR Biosciences).
Serum sensitivity assay
Bacterial susceptibility to normal duck sera was tested as described elsewhere , with slight modifications. Briefly, normal complement-sufficient duck sera (without anti-R. anatipestifer antibody) were collected from 18-old-day healthy Cherry Valley ducks, pooled, and filter-sterilized (0.22 μm). Pooled duck sera were diluted to 12.5%, 25%, and 50% (v/v) in pH 7.2 PBS. Then, 10-μL aliquots of bacterial suspension containing 106 CFU in PBS was added to 190 μL of serial diluted duck sera, 100% heat-inactivated serum, or PBS alone, respectively. The reaction mixtures were incubated at 37 °C under an atmosphere of 5% CO2 for 30 min. Afterward, the bacterial culture was serial diluted by tenfold and the bacterial count in each sample was calculated by plating onto TSA plates at 37 °C under an atmosphere of 5% CO2 for 36 h. The experiment was performed in triplicate.
Adhesion and invasion assays
Bacterial adhesion and invasion assays were performed using Vero cells (ATCC CCL-81) as described elsewhere . Briefly, Vero cells (105/well) were grown in 24-well tissue culture trays in Dulbecco’s modified Eagle medium (DMEM), containing 10% fetal bovine serum (Biowest, Nuaillé, France). Prior to infection, confluent monolayers of Vero cells were rinsed three times with sterile PBS and then infected with the WT strain CH3 or the mutant strain RA1062 at a multiplicity of infection of 100. Infected cells were incubated at 37 °C under an atmosphere of 5% CO2 for 1.5 h. Non-adherent bacteria were removed by washing three times with sterile PBS. The cell-adherent bacteria were enumerated following dispersion with PBS-0.1% trypsin, serial diluted tenfold, and spread onto TSA plates to estimate the number of bacteria (CFU) adhering to the cell monolayers. For the invasion assay, the cell culture, bacterial infection, and plating procedures were performed as described for the adherence assay. After bacterial infection, 1 mL of DMEM medium containing 100 μg/mL of gentamicin was added to each well and the plate was incubated at 37 °C under an atmosphere of 5% CO2 for additional 1 h to kill extracellular bacteria. All of samples were assayed in triplicate and the assay was conducted independently three times.
Bacterial virulence determination
The median lethal dose (LD50) of the WT strain CH3 and the mutant strain RA1062 was determined using 18-day-old Cherry Valley ducks as described elsewhere . The ducks were randomly divided into five groups (8 ducks/group), and inoculated intramuscularly with the appropriate bacterial strain at a dose of 106, 107, 108, 109, or 1010 CFU, respectively. Moribund ducks with clinical signs of diarrhoea, pyrexia, anorexia, stunted growth, respiratory signs, neurological abnormalities, or ocular signs were euthanized humanely with an intravenous injection of sodium pentobarbital at a dose of 120 mg/kg and counted as dead. On post-mortem, a yellow-white exudate and congestion can be seen throughout the body. Ducks were monitored daily for death rate for a period of 7 days post-infection to calculate LD50 value according to the improved Karber’s method . The diagnosis was finally confirmed by bacterial isolation from liver, spleen and brain samples of dead ducks.
To evaluate bacterial survival in vivo, 18-day-old Cherry Valley ducks were divided randomly into two groups (6 ducks/group) and injected intramuscularly with 5 × 108 CFU of the WT strain CH3 or the mutant strain RA1062. Blood samples were collected at 6, 12, 24, 48, and 72 h post-infection (hpi) (six ducks per group at each time point), diluted tenfold to appropriate concentrations, and plated in triplicate on TSA for bacterial counting.
Illumina sequencing for RNA-Seq and differential gene expression analyses
The quantity and quality of total RNA were assessed using the Agilent RNA 6000 Nano Kit and the Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). Total RNA was treated with the Ribo-Zero Magnetic Gold Kit (Epicentre® (an Illumina company), Madison, WI, USA) to remove ribosome RNA, then libraries were constructed using the TruSeq RNA Sample Prep Kit v2 (Illumina, Inc., San Diego, CA, USA) in accordance with the manufacturer’s instructions. The quality and quantity of the libraries were assessed by two methods: the distribution of the fragment size was checked using the Agilent 2100 bioanalyzer (Agilent DNA 1000 Reagents) and quantified using real-time qPCR (TaqMan Probe). The qualified libraries were amplified on cBot to generate the cluster on the flowcell (TruSeq PE Cluster Kit V3–cBot–HS, Illumina). Then, the amplified flowcell was sequenced via the paired-end method on the HiSeq2000 System (TruSeq SBS KIT-HS V3, Illumina), with a read length of 90 bp, which is the most common sequencing strategy.
The complete libraries were sequenced for 100 cycles on the Illumina HiSeq2000 system as described elsewhere . Image analysis and base calling were performed using Solexa pipeline version 1.8 (Off-Line Base Caller software, version 1.8) . Cleaned reads were aligned to the genome of R. anatipestifer strain CH3 using RNA Sequel software . Transcript levels were calculated as reads per kilobase of cDNA per million fragments mapped. Differentially expressed genes were analyzed using Cufflinks software (version 2.1.1) with significance at a fold change cut-off of 2.0  and considered statistically significant if the fold change was > 2.0 and false discovery rate of < 0.001.
Real-time quantitative PCR (qPCR) analysis
qPCR was performed to confirm the transcriptional levels of differentially expressed genes obtained by RNA-Seq analysis. Gene-specific primers (Table 1) were designed using primer3 online software version 0.4.0 . The expression levels of the l-lactate dehydrogenase encoding gene (ldh) were measured using the primer pair RA ldh-F/RA ldh-R (Table 1) and used as an internal control . Total RNA was isolated from the WT strain CH3 and the mutant strain RA1062 using TRIzol reagent (Invitrogen Corporation, Carlsbad, CA, USA), according to the manufacturer’s instructions. All RNA samples were treated with the TURBO DNA-free kit (Ambion/Life Sciences, Grand Island, NY, USA) to remove DNA contamination. cDNA was synthesized using PrimeScript RT Master Mix (Takara). qPCR was conducted using Go Taq qPCR Master Mix (Promega, Fitchburg, WI, USA) with the following parameters: 95 °C for 2 min, 40 cycles of 95 °C for 15 s, 55 °C for 15 s and 68 °C for 20 s, followed by one cycle of 95 °C for 15 s, 60 °C for 15 s and 95 °C for 15 s. Reactions were performed in triplicate and run on the Mastercycler ep realplex4 apparatus (Eppendorf AG, Hamburg, Germany). Quantification of transcriptional level was calculated according to the 2−ΔΔCt method.
Vaccination and challenge assays
To determine whether the mutant strain RA1062 conveyed effective cross-protection among different R. anatipestifer strains, the inactivated RA1062 and CH3 vaccines were prepared using the mutant strain RA1062 and the WT strain CH3, respectively. Briefly, strains RA1062 and CH3 were cultured separately in TSB at 37 °C for 10 h with shaking. The bacterial culture was then adjusted to 1010 CFU/mL and inactivated with 0.4% (vol/vol) formalin at 37 °C for 16 h, respectively. The vaccine was made by blending 3 volume inactivated strains (RA1062 or CH3) and 7 volumes of Montanide ISA 70 VG adjuvant (Seppic Shanghai Special Chemical Corporation, Shanghai, China) according to the manufacturer’s protocol. Each duck was subcutaneously injected in the neck with 0.3 mL of the vaccine containing 109 CFU bacterial cells of RA1062 or CH3 respectively. Cherry Valley ducks were divided randomly into three groups (24 ducks per group). The ducks in group 1 received two subcutaneous injections of the inactivated RA1062 vaccine, at a dose of 109 CFU. Each duck in group 2 received two subcutaneous injections of the inactivated CH3 vaccine, and those in group 3 received two subcutaneous injections of saline in adjuvant for use as controls. At 14 days after the second immunization, eight ducks from each group were challenged with R. anatipestifer strain WJ4 (serotype 1), Yb2 (serotype 2), or HXb2 (serotype 10) at 10 LD50, respectively. Ducks were monitored daily for clinical symptoms and death until 7 days post-infection. The protection rate was calculated as follows: [1 − (no. of dead ducks per group/total no. of ducks per group)] × 100. The experiment was repeated for three times.
All statistical analyses were performed using GraphPad Prism, version 5.0 for Windows software (GraphPad Software Inc., La Jolla, CA, USA). Adhesion and invasion assays, bacterial growth curves, bacterial loads in the blood of ducks, serum sensitivity assays, and RT-qPCR were two tailed, and a p-value of < 0.05 was considered statistically significant. Analysis of variance was used for comparisons of multiple groups.
Identification of the mutant strain RA1062
The mutant strain RA1062, which lacked reactivity with the anti-CH3 LPS mAb, was obtained by screening the transposon library using an indirect ELISA and identified by PCR amplification using the primer pairs 16S rRNA F/16S rRNA R, Erm-F/Erm-R, and RA1062-F/RA1062-R. As shown in Figure 1A, a 744-bp fragment of 16S rRNA gene was amplified from the WT strain CH3 (lane 1) and the mutant strain RA 1062 (lane 2), a 678-bp fragment of the M949_RS01035 gene was amplified from the WT strain CH3 (lane 4), but not from the mutant strain RA1062 due to the transposon insertion (lane 5), and a 644-bp fragment of the erm gene (contained in the Tn4351 transposon) was amplified from the mutant strain RA1062 (lane 8), but not from the WT strain CH3 (lane 7). The mutant strain RA1062 was confirmed to contain a single Tn4351 insertion in the chromosomal DNA by Southern blot analysis (Figure 1B, lane 2).
The transposon was inserted at nucleotide position 318 bp of the M949_RS01035 gene, which is 678 nucleotides in length and encodes intramembrane metalloprotease of the CPBP (CAAX proteases and bacteriocin-processing enzymes) family, which consists of 225 amino acids (Figure 1C). BLAST analysis showed that M949_RS01035 gene exists in R. anatipestifer serotype 1 strains CH3 and CH-1, as well as serotype 10 strain HXb2 (Additional file 1). qPCR analysis further confirmed that transcription of the M949_RS01035 gene was abolished in the mutant strain RA1062 (Figure 1D). Further investigation showed that transcription of chromosomally upstream M949_RS10475 gene, which encodes a hypothetical protein, and the downstream M949_RS01030 gene, which encodes a “prevent host death” protein, had no significant changes, as compared with the WT strain CH3 (Figure 1D).
Determination of bacterial growth curves
Bacterial growth curves of the WT strain CH3 and the mutant strain RA1062 in TSB medium were constructed. Compared to that of the WT strain CH3, the growth rate of the mutant strain RA1062 was similar at the early growth stage, but then significantly decreased at the late growth stage. After 12 h in TSB, growth of the WT strain CH3 reached a plateau of 5.12 × 109 CFU/mL, while the mutant strain RA1062 reached a plateau of 1.99 × 109 CFU/mL after 8 h in TSB. The bacterial numbers of the mutant strain RA1062 at the plateau was about 2.79-fold less than that of the WT strain CH3 (Figure 2).
Bacterial phenotype characterization by crystal violet staining
Crystal violet staining was used to characterize the bacterial phenotypes of the WT strain CH3 and the mutant strain RA1062. The result indicated that the WT strain CH3 was not stained by crystal violet (Figure 3A), while the mutant strain RA1062 was (Figure 3B), suggesting a change in its phenotype of the mutant strain RA1062.
Analysis of the bacterial LPS by silver staining and Western blot
LPS was purified from the WT strain CH3 and the mutant strain RA1062, then subjected to SDS-PAGE followed by silver staining and Western blot analysis. As shown in Figure 4A, LPS purified from the WT strain CH3 displayed a ladder-like pattern around 70 kDa (lane 2), while the ladder-like pattern was deficient in the mutant strain RA1062 LPS (lane 3) in the silver staining. Western blot analysis using the anti-CH3 LPS mAb detected a ladder-like pattern of O-antigen repeats in the CH3 strain (lane 2), but the pattern was defective in the mutant strain RA 1062 (Figure 4B). These results suggest that the R. anatipestifer M949_RS01035 gene is involved in LPS O-antigen biosynthesis.
Adhesion and invasion assays
To determine the role of the M949_RS01035 gene in bacterial adherence and invasion, we compared bacterial adhesion and invasion abilities of the WT strain CH3 and the mutant RA1062 on Vero cells. The results showed that the mutant strain RA1062 displayed approximately 21.92-fold enhanced adherence ability (Figure 5A) and 25.18-fold increased invasion capacity (Figure 5B) in comparison with those of the WT strain CH3 (p < 0.001).
The mutant strain RA1062 was more sensitive to normal duck sera
To determine whether the M949_RS01035 gene is involved in serum resistance of the WT strain CH3, we compared the abilities of the WT strain CH3 and the mutant strain RA1062 to resist the complement-mediated killing. The results showed that 25% diluted serum was effective in killing the mutant strain RA1062, but not the WT strain CH3, indicating that the mutant strain RA1062 was more sensitive to normal duck sera than the WT strain CH3 (Figure 6).
Determination of bacterial virulence
Bacterial virulence was evaluated based on LD50 determination using 18-day-old Cherry Valley ducks. The LD50 of the mutant strain RA1062 was 2.74 × 1010 CFU, which was more than 365-fold greater attenuated virulence than that of the WT strain CH3 (7.50 × 107 CFU). In addition, to confirm the role of the M949_RS01035 gene in systemic invasion and dissemination, bacterial loading in the blood of ducks infected with the WT strain and the mutant strain RA1062 was conducted. As shown in Figure 7, the bacterial loads in the blood of ducks infected with the WT strain CH3 continued to increase up to 72 hpi, while that in the blood of ducks infected with the mutant strain RA1062 decreased from 48 hpi, indicating a significant decrease in comparison with ducks infected with the WT strain CH3 at 48 and 72 hpi (p < 0.001). This experimental result further confirmed the attenuated virulence of the mutant strain RA1062.
Identification of differentially expressed genes
Strand-specific Illumina RNA-Seq analysis was used to identify differentially expressed genes between the WT strain CH3 and the mutant strain RA1062. In total, 12 genes were up-regulated and nine were down-regulated in the mutant strain RA1062 in comparison to the WT strain CH3 (Table 2), which were differentially expressed over fivefold based on RNA-Seq analysis. qPCR further confirmed that the M949_RS07300 and M949_RS04680 genes were up-regulated, and the M949_RS07580 and M949_RS03025 genes were down-regulated by more than fivefold at the transcriptional level. The proteins encoded by the M949_RS07300, M949_RS04680, and M949_RS03025 genes were annotated as TonB-dependent receptors. The coding product of the M949_RS07580 gene was annotated as the polysaccharide biosynthesis protein CapD, which is responsible for bacterial capsule polysaccharide biosynthesis. These results indicated that the M949_RS01035 gene regulated several genes responsible for bacterial polysaccharide biosynthesis.
To investigate the cross-protection against R. anatipestifer serotypes 1, 2, and 10, the RA1062-vaccinated ducks were challenged with virulent R. anatipestifer strains WJ4 (serotype 1), Yb2 (serotype 2), and HXb2 (serotype 10) at 10 LD50, respectively. Protection from challenge with strains WJ4, Yb2, and HXb2 was greater in RA1062-vaccinated ducks than those vaccinated with CH3. In contrast, all non-vaccinated ducks were dead within 7 days post-challenge. Three experimental replicates showed similar cross-protection results (Table 3). These results indicated that RA1062 provided cross-protection against challenge with R. anatipestifer serotypes 1, 2, and 10, which further confirmed its altered antigenicity.
This study was initially designed to discover genes of R. anatipestifer strain CH3 involved in LPS biosynthesis and to explore potential use of the mutant for development of a cross-protective vaccine. The mutant strain RA1062 was obtained by screening the random Tn4351 transposon mutant library with an indirect ELISA, which lacked reactivity with the anti-CH3 LPS mAb. Further investigations revealed that the phenotype of the mutant strain was altered with a decreased growth rate in TSB, increased sensitivity to normal duck serum, and increased bacterial adhesion and invasion capacities on Vero cells, as compared to the WT strain CH3. Sequence analysis revealed that the transposon was inserted into the M949_RS01035 gene of strain CH3 at 318 bp.
The protein encoded by the M949_RS01035 gene is annotated as an intramembrane metalloprotease of the CPBP (CAAX proteases and bacteriocin-processing enzymes) family, which exhibited a high degree of similarity with the Abi family proteins of Lactobacillus plantarum. The members of the Abi protein family are reportedly involved in bacteriocin self-immunity in streptococci . A BLAST search indicated that the M949_RS01035 gene exists in R. anatipestifer serotypes 1 strains CH3 and CH-1, as well as serotype 10 strain HXb2. Silver staining indicated a marked deficiency of the ladder-like pattern in RA1062 mutant LPS, as compared to CH3 LPS. Western-blot analysis further confirmed that the modified LPS was less reactive with the anti-CH3 LPS mAb, suggesting that the M949_RS01035 gene is associated with the LPS O-antigen biosynthesis.
It has been reported that the O-antigen side-chains are main components of LPS molecules of E. coli and Salmonella typhimurium involved in serum resistance by virtue of their anti-complement activities . In the present study, the mutant strain RA1062 exhibited increased sensitivity to normal duck serum. Therefore, we speculated that the higher sensitivity of the mutant strain RA1062 to normal duck serum may due to a defect to the O-antigen, leading to easier bacteria killing by serum. It is also well established that the O-antigen plays a significant role in bacterial virulence . Animal experiments further demonstrated that the virulence of mutant strain RA1062 was significantly decreased, as compared to the WT strain CH3. These experimental data demonstrate that bacterial O-antigen is not only involved in bacterial resistance to complement-mediated killing, but also an important virulence factor. It is established that bacterial adherence to the cell surface is the first step in the pathogenesis of microbial infections and functions as a critical determinant in colonization . In this study, we found that both adherence and invasion capacities of the mutant strain RA1062 to Vero cells were significantly increased, as compared to those of the WT strain CH3. This phenomenon is consistent with rough mutants of Brucella abortus, which are taken up in greater numbers by macrophages than the smooth parental strains . The enhanced ability of the mutant strain RA1062 to adhere to and invade Vero cells may also be dependent on the rough LPS phenotype, which employs different adherence and invasion mechanisms, as described for B. abortus .
Various antibiotics are currently used to prevent and control R. anatipestifer infection in ducks. However, the emergence of drug-resistant and multi-drug resistant strains due to the overuse of antibiotics has become a major barrier in R. anatipestifer treatment [2, 40]. Therefore, immunization of ducks with vaccine provides a valuable alternative to the use of antibiotics. Our study showed that ducks immunized with the mutant strain RA1062 were protected from challenge with R. anatipestifer strains WJ4 (serotype 1), Yb2 (serotype 2), and HXb2 (serotype 10) at ≥ 8/8, 6/8 and 8/8, respectively, suggesting that the RA1062 strain is an effective cross-protective vaccine candidate. A previous study reported that conserved LPS epitope(s) exist among different species and genera of non-enteric, Gram-negative, human pathogens . The ability of the mutant strain RA1062 to provide broad cross-protection may be due to the deficiency of LPS O-antigen, which leads to the presence of some conserved LPS epitopes. Another possibility is that some cross-immunogenic components of the R. anatipestifer outer membrane, such as outer membrane proteins, are exposed at the bacterial surface of the mutant strain RA1062. However, based on current experiments, the mechanism of the cross-protection induced by the mutant strain RA1062 vaccination remains unclear, thus further work is needed to clarify the nature of the cross-protection.
RNA-Seq has rapidly become the method of choice for the study of differential gene expression, as it enables the investigation and comparison of gene expression levels at unprecedented resolution . In our current RNA-Seq study, we found that 12 genes were up-regulated and nine were down-regulated by over fivefold in the mutant strain RA1062. qPCR verification further confirmed that two genes were up-regulated and two were down-regulated by over 5-fold. The proteins encoded by two up-regulated genes (M949_RS07300 and M949_RS04680) and one down–regulated gene (M949_RS03025) are TonB-dependent receptors. TonB-dependent receptors are members of a family of beta-barrel proteins located on the outer membrane of Gram-negative bacteria, which are known to be mainly involved in iron or vitamin B12 uptake . The crystallographic structures of two TonB-dependent receptors (FhuA and FepA) have recently been determined . A recent study has identified the TonB-dependent receptor TbdR1 in R. anatipestifer strain CH3, which is a cross-immunogenic protein among strains with different serotypes . Further investigation demonstrated that TbdR1 had a crucial role in the process of iron acquisition and complicated pathogenesis of R. anatipestifer strain CH3 . However, how these TonB-dependent receptors are regulated by the M949_RS01035 gene remains unknown, thus further studies are needed to elucidate the nature of this phenomenon. The M949_RS07580 gene, which encodes the polysaccharide biosynthesis protein CapD, was significantly down-regulated by over 30-fold, which was consistent with the findings of our previous study . It has been established that the polysaccharide biosynthesis protein CapD is a novel pathogenicity-associated determinant that is involved in the serum-resistance ability of H. parasuis and E. faecium [47, 48]. Protein sequence analysis indicated that there is a NAD-binding domain in the polysaccharide biosynthesis protein CapD, which is reportedly associated with LPS O-antigen biosynthesis . Additionally, the polysaccharide biosynthesis protein CapD is responsible for the biosynthesis of type 1 capsular polysaccharides in Staphylococcus aureus . Previous work has demonstrated that Group 1 and 4 capsules are related to LPS O-antigens . Based on our current knowledge, we speculated that the loss of the LPS O-antigen in the mutant strain RA1062 may be associated with the down-regulation of the M949_RS07580 gene, although further investigations are needed to elucidate the underlying mechanism.
In conclusion, we demonstrated that the M949_RS01035 gene is involved in bacterial phenotype, virulence, and gene regulation in R. anatipestifer. The mutant strain RA1062 could be used as a novel cross-protective vaccine candidate for further vaccine development.
enzyme-linked immunosorbent assay
polymerase chain reaction
colony forming unit
sodium dodecyl sulfate polyacrylamide gel electrophoresis
- LD50 :
median lethal doses
Segers P, Mannheim W, Vancanneyt M, De Brandt K, Hinz KH, Kersters K, Vandamme P (1993) Riemerella anatipestifer gen. nov., comb. nov., the causative agent of septicemia anserum exsudativa, and its phylogenetic affiliation within the Flavobacterium-Cytophaga rRNA homology group. Int J Syst Bacteriol 43:768–776
Sun N, Liu J-H, Yang F, Lin D-C, Li G-H, Chen Z-L, Zeng Z-L (2012) Molecular characterization of the antimicrobial resistance of Riemerella anatipestifer isolated from ducks. Vet Microbiol 158:376–383
Sandhu T (2008) Riemerella anatipestifer infection. In: Saif YM, Fadly AM, Glisson JR, McDougald LR, Nolan LK, Swane DE (eds) Diseases of poultry, 12th edn. Blackwell Publishing Ltd, Hoboken
Pathanasophon P, Phuektes P, Tanticharoenyos T, Narongsak W, Sawada T (2002) A potential new serotype of Riemerella anatipestifer isolated from ducks in Thailand. Avian Pathol 31:267–270
Pathanasophon P, Sawada T, Tanticharoenyos T (1995) New serotypes of Riemerella anatipestifer isolated from ducks in Thailand. Avian Pathol 24:195–199
Hu Q, Han X, Zhou X, Ding S, Ding C, Yu S (2010) Characterization of biofilm formation by Riemerella anatipestifer. Vet Microbiol 144:429–436
Han X, Hu Q, Ding S, Chen W, Ding C, He L, Wang X, Ding J, Yu S (2012) Identification and immunological characteristics of chaperonin GroEL in Riemerella anatipestifer. Appl Microbiol Biotechnol 93:1197–1205
Haiwen L (2013) Development and evaluation of a trivalent Riemerella anatipestifer-inactivated vaccine. Clin Vaccine Immunol 5:691–697
Chang CF, Hung PE, Chang YF (1998) Molecular characterization of a plasmid isolated from Riemerella anatipestifer. Avian Pathol 27:339–345
Hu Q, Han X, Zhou X, Ding C, Zhu Y, Yu S (2011) OmpA is a virulence factor of Riemerella anatipestifer. Vet Microbiol 150:278–283
Murray GL, Attridge SR, Morona R (2003) Regulation of Salmonella typhimurium lipopolysaccharide O antigen chain length is required for virulence; identification of FepE as a second Wzz. Mol Microbiol 47:1395–1406
Murray GL, Attridge SR, Morona R (2006) Altering the length of the lipopolysaccharide O antigen has an impact on the interaction of Salmonella enterica serovar typhimurium with macrophages and complement. J Bacteriol 188:2735–2739
Duerr CU, Zenk SF, Chassin C, Pott J, Gutle D, Hensel M, Hornef MW (2009) O-antigen delays lipopolysaccharide recognition and impairs antibacterial host defense in murine intestinal epithelial cells. PLoS Pathog 5:e1000567
Saldias MS, Ortega X, Valvano MA (2009) Burkholderia cenocepacia O antigen lipopolysaccharide prevents phagocytosis by macrophages and adhesion to epithelial cells. J Med Microbiol 58:1542–1548
Levine MM, Kotloff KL, Barry EM, Pasetti MF, Sztein MB (2007) Clinical trials of Shigella vaccines: two steps forward and one step back on a long, hard road. Nat Rev Microbiol 5:540–553
Austin EA, Graves JF, Hite LA, Parker CT, Schnaitman CA (1990) Genetic analysis of lipopolysaccharide core biosynthesis by Escherichia coli K-12: insertion mutagenesis of the rfa locus. J Bacteriol 172:5312–5325
Wang X, Quinn PJ (2010) Lipopolysaccharide: biosynthetic pathway and structure modification. Prog Lipid Res 49:97–107
Zou J, Wang X, Ding C, Tian M, Han X, Wang S, Yu S (2015) Characterization and cross-protection evaluation of M949_1603 gene deletion Riemerella anatipestifer mutant RA-M1. Appl Microbiol Biotechnol 99:10107–10116
Dou Y, Wang X, Yu G, Wang S, Tian M, Qi J, Li T, Ding C, Yu S (2017) Disruption of the M949_RS01915 gene changed the bacterial lipopolysaccharide pattern, pathogenicity and gene expression of Riemerella anatipestifer. Vet Res 48:6
Yu G, Wang X, Dou Y, Wang S, Tian M, Qi J, Li T, Ding C, Wu Y, Yu S (2016) Riemerella anatipestifer M949_1360 gene functions on the lipopolysaccharide biosynthesis and bacterial virulence. PLoS One 11:e0160708
Zou J, Wang X, Tian M, Cao S, Hou W, Wang S, Han X, Ding C, Yu S (2015) The M949_1556 gene plays a role on the bacterial antigenicity and pathogenicity of Riemerella anatipestifer. Vet Microbiol 177:193–200
Wang X, Ding C, Wang S, Han X, Hou W, Yue J, Zou J, Yu S (2014) The AS87_04050 gene is involved in bacterial lipopolysaccharide biosynthesis and pathogenicity of Riemerella anatipestifer. PLoS One 9:e109962
Hu Q, Zhu Y, Tu J, Yin Y, Wang X, Han X, Ding C, Zhang B, Yu S (2012) Identification of the genes involved in Riemerella anatipestifer biofilm formation by random transposon mutagenesis. PLoS One 7:e39805
Alvarez B, Secades P, Prieto M, McBride MJ, Guijarro JA (2006) A mutation in Flavobacterium psychrophilum tlpB inhibits gliding motility and induces biofilm formation. Appl Environ Microbiol 72:4044–4053
Basic Local Alignment Search Tool. https://blast.ncbi.nlm.nih.gov/Blast.cgi. Accessed 5 July 2016.
Tsai CM, Frasch CE (1982) A sensitive silver stain for detecting lipopolysaccharides in polyacrylamide gels. Anal Biochem 119:115–119
McQuillen DP, Gulati S, Rice PA (1994) Complement-mediated bacterial killing assays. Methods Enzymol 236:137–147
Reed LJ, Muench H (1938) A simple method of estimating fifty per cent endpoints. Am J Epidemiol 27:493–497
Wang X, Yue J, Ding C, Wang S, Liu B, Tian M, Yu S (2016) Deletion of AS87_03730 gene changed the bacterial virulence and gene expression of Riemerella anatipestifer. Sci Rep 6:22438
Whiteford N, Skelly T, Curtis C, Ritchie ME, Lohr A, Zaranek AW, Abnizova I, Brown C (2009) Swift: primary data analysis for the Illumina Solexa sequencing platform. Bioinformatics 25:2194–2199
Wilson GW, Stein LD (2015) RNASequel: accurate and repeat tolerant realignment of RNA-seq reads. Nucleic Acids Res 43:e122
Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L (2010) Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol 28:511–515
Primer 3 online software version 0.4.0. http://bioinfo.ut.ee/primer3-0.4.0/. Accessed 8 Sept 2016.
Huggett J, Dheda K, Bustin S, Zumla A (2005) Real-time RT-PCR normalisation; strategies and considerations. Genes Immun 6:279–284
Kjos M, Snipen L, Salehian Z, Nes IF, Diep DB (2010) The abi proteins and their involvement in bacteriocin self-immunity. J Bacteriol 192:2068–2076
Tomas JM, Ciurana B, Benedi VJ, Juarez A (1988) Role of lipopolysaccharide and complement in susceptibility of Escherichia coli and Salmonella typhimurium to non-immune serum. J Gen Microbiol 134:1009–1016
Sarkar S, Ulett GC, Totsika M, Phan M-D, Schembri MA (2014) Role of capsule and O antigen in the virulence of uropathogenic Escherichia coli. PLoS One 9:e94786
Pei J, Turse JE, Wu Q, Ficht TA (2006) Brucella abortus rough mutants induce macrophage oncosis that requires bacterial protein synthesis and direct interaction with the macrophage. Infect Immun 74:2667–2675
Zhang M, Han X, Liu H, Tian M, Ding C, Song J, Sun X, Liu Z, Yu S (2013) Inactivation of the ABC transporter ATPase gene in Brucella abortus strain 2308 attenuated the virulence of the bacteria. Vet Microbiol 164:322–329
Chen Y-P, Lee S-H, Chou C-H, Tsai H-J (2012) Detection of florfenicol resistance genes in Riemerella anatipestifer isolated from ducks and geese. Vet Microbiol 154:325–331
Campagnari AA, Spinola SM, Lesse AJ, Kwaik YA, Mandrell RE, Apicella MA (1990) Lipooligosaccharide epitopes shared among Gram-negative non-enteric mucosal pathogens. Microb Pathog 8:353–362
Finotello F, Di Camillo B (2015) Measuring differential gene expression with RNA-seq: challenges and strategies for data analysis. Brief Funct Genomics 14:130–142
Blanvillain S, Meyer D, Boulanger A, Lautier M, Guynet C, Denancé N, Vasse J, Lauber E, Arlat M (2007) Plant carbohydrate scavenging through TonB-dependent receptors: a feature shared by phytopathogenic and aquatic bacteria. PLoS One 2:e224
Ferguson AD, Deisenhofer J (2002) TonB-dependent receptors—structural perspectives. Biochim Biophys Acta 1565:318–332
Hu Q, Ding C, Tu J, Wang X, Han X, Duan Y, Yu S (2012) Immunoproteomics analysis of whole cell bacterial proteins of Riemerella anatipestifer. Vet Microbiol 157:428–438
Lu F, Miao S, Tu J, Ni X, Xing L, Yu H, Pan L, Hu Q (2013) The role of TonB-dependent receptor TbdR1 in Riemerella anatipestifer in iron acquisition and virulence. Vet Microbiol 167:713–718
Ali L, Spiess M, Wobser D, Rodriguez M, Blum HE, Sakιnç T (2016) Identification and functional characterization of the putative polysaccharide biosynthesis protein (CapD) of Enterococcus faecium U0317. Infect Genet Evol 37:215–224
Wang X, Xu X, Wu Y, Li L, Cao R, Cai X, Chen H (2013) Polysaccharide biosynthesis protein CapD is a novel pathogenicity-associated determinant of Haemophilus parasuis involved in serum-resistance ability. Vet Microbiol 164:184–189
Fallarino A, Mavrangelos C, Stroeher UH, Manning PA (1997) Identification of additional genes required for O-antigen biosynthesis in Vibrio cholerae O1. J Bacteriol 179:2147–2153
Lin WS, Cunneen T, Lee CY (1994) Sequence analysis and molecular characterization of genes required for the biosynthesis of type 1 capsular polysaccharide in Staphylococcus aureus. J Bacteriol 176:7005–7016
Whitfield C (2006) Biosynthesis and assembly of capsular polysaccharides in Escherichia coli. Annu Rev Biochem 75:39–68
The authors declare that they have no competing interests.
SY planned the experiments. YD, GY, and XW conducted the experiments. SW, TL, MT, JQ, CD, and SY analyzed and discussed the experimental results. YD and SY wrote the manuscript. All authors read and approved the final manuscript.
We thank Ms. Yingfang Lv for supporting the animal experiments.
This work was supported by the National Key Research and Development Program of China (No. 2016YFD0500805), the Shanghai Key Project on Agricultural Development through Science and Technology (2016HNG4-1) and Co-innovation of Science and Technology Innovation Project in Chinese Academy of Agricultural Sciences (CAAS-XTCX2016011-04-8).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.