- Research article
- Open Access
Clonal analysis and virulent traits of pathogenic extraintestinal Escherichia coliisolates from swine in China
BMC Veterinary Researchvolume 8, Article number: 140 (2012)
Extraintestinal pathogenic Escherichia coli (ExPEC) can cause a variety of infections outside the gastrointestinal tract in humans and animals. Infections due to swine ExPECs have been occurring with increasing frequency in China. These ExPECs may now be considered a new food-borne pathogen that causes cross-infections between humans and pigs. Knowledge of the clonal structure and virulence genes is needed as a framework to improve the understanding of phylogenetic traits of porcine ExPECs.
Multilocus sequence typing (MLST) data showed that the isolates investigated in this study could be placed into four main clonal complexes, designated as CC10, CC1687, CC88 and CC58. Strains within CC10 were classified as phylogroup A, and these accounted for most of our porcine ExPEC isolates. Isolates in the CC1687 clonal complex, formed by new sequence types (STs), was classified as phylogroup D, with CC88 isolates considered as B2 and CC58 isolates as B1. Porcine ExPECs in these four clonal complexes demonstrated significantly different virulence gene patterns. A few porcine ExPECs were indentified in phylogroup B2, the phylogroup in which human ExPECs mainly exist. However some STs in the four clonal groups of porcine ExPECs were reported to cause extraintestinal infections in human, based on data in the MLST database.
Porcine ExPECs have different virulence gene patterns for different clonal complexes. However, these strains are mostly fell in phylogenentic phylogroup A, B1 and D, which is different from human ExPECs that concentrate in phylogroup B2. Our findings provide a better understanding relating to the clonal structure of ExPECs in diseased pigs and indicate a need to re-evaluate their contribution to human ExPEC diseases.
Escherichia coli (E. coli) are common bacteria that colonize the intestinal tract as normal flora and are known to cause a variety of diseases. Extraintestinal pathogenic E. coli (ExPEC) is epidemiologically and genetically different from intestinal pathogenic or commensal E. coli. A typical ExPEC does not cause enteric diseases; however, it can asymptomatically colonize the intestinal tract, or cause extraintestinal diseases in both normal and compromised hosts [1, 2]. The major virulence factors present in ExPECs are distinct from those found in other E. coli strains; ExPECs possess a variety of the combination of virulence-associated genes instead of a common virulence mechanism to upset the host defense system . ExPECs of human origin were defined by Johnson et al.  as E. coli isolates containing two or more virulence markers: papA/papC, sfa/foc, afa/dra, kpsMTII and iutA. The same criterion was also used to define ExPECs isolated from other animal species because fewer systematic studies were carried out on ExPEC from these animals with respect to virulence genes.
Within the Chinese pig industry, porcine ExPEC diseases have become more prevalent. This has corresponded to an increase in the number an increase in the number of resistant bacteria and immunocompromised diseases . The frequency of isolation of ExPECs from diseased pigs increased from 3.1% in 2004 to 14.6% in 2007. Moreover, strains possessing similar virulence profiles and serogroups to human ExPECs have been detected from our collection of porcine ExPECs, suggesting that these pathogens may have strong zoonotic potential . The virulence and zoonotic potential of porcine ExPECs may be linked to a specific clonal lineage within phylogenetic groups. It is important to understand the population structure of porcine ExPECs. However, the relevant information is rarely available.
Multilocus sequence typing (MLST) is a powerful method that can be applied to comparing genotypes and global surveillance of pathogens. MLST is based on the sequencing of 6–7 selected housekeeping genes and identification of polymorphic nucleotide sites. A combination of the alleles at the 6–7 loci provides a unique sequence types (ST) identifier that can be used to discriminate strains. A new implementation of algorithm, eBURST, has been shown to be effective in the analysis of MLST data. Related sample strains are grouped into a clonal complex (CC), and the STs which demonstrated the greatest number of single-locus variants (SLVs) were assigned as the founder. Compared with a dendrogram, eBURST provides more information about founder genotypes and patterns of the evolutionary descent.
In the present study, a collection of 81 ExPEC strains isolated from diseased pigs within China from 2006–2007 were subjected to MLST, analyzed via serogrouping and by virulence gene typing with the aim of improving our understanding of the clonal structure and virulence traits of swine ExPECs.
MLST and eBURST analysis of porcine extraintestinal E.coli
The CCs and the distributions of serogroups and virulence genes are shown in Table 1. Seven housekeeping genes of all 81 porcine ExPEC strains were successfully amplified and sequenced. A total of 43 STs was found, of which 29 were reported previously. Fourteen STs from 30 strains were identified as new STs and designated ST1687, 1688, 1690, 1691, 1692, 1694, 1696–1702, and 1731 (http://mlst.ucc.ie/mlst/dbs/Ecoli/GetTableInfo_html).
The two most frequent STs were ST1687 (19.8%) and ST10 (7.4%). eBURST analysis of 81 strains showed that there were five CCs, CC10, CC1687, CC88, CC46 and CC206, and two doublets. Nine isolates, including two of the newly identified STs (ST1698 and ST1701), were defined as a singletons. Further comparison of the porcine ExPECs with the entire MLST database indicated that CC46, CC206 and several singletons were derivatives of CC10, and that five singletons (ST224, ST58, ST1642, ST602 and ST683) were related to strains from CC58 (Figure 1). These related STs and CCs, which did not fall into a single CC but were grouped together according to the MLST database, were classified as CC10 and CC58 respectively in this study (Table 1). The strains within CC10 and their derivatives consisted of 26 isolates (32.1%), in which ST10 and ST46 occurred most frequently and accounted for 38.5% (n = 10) of STs.
E. coli can be classified as phylogroup A, B1, B2 or D according to the phylogenetic relationship of the sequences. Twenty-six isolates belonged to phylogroup A, 20 to phylogroup B1, 10 to phylogroup B2, and 25 to phylogroup D. The isolates belonging to phylogroups B2 and D accounted for 43.2% of all isolates. This was slightly lower than the proportion of the, 8isolates belonging to phylogroups A and B1 combined. The strains falling into phylogroup B2 accounted for 12.3% of all isolates. No significant correlation between serogroup and phylogenetic group was found, except that all isolates of serogroup O11 with ST1687, one of the newly identified STs, belonged to phylogenetic group D.
Distribution of serogroups in CCs
No difference in distribution of the five most frequent serogroups (O11, O8, O138, O161, and O101) among CCs was observed, except that serogroup O101 belonged to CC10, and serogroup O11 was predominant in CC1687. In CC58, doublets and singletons from several serogroups were dispersed among isolates of the same CC or ST, indicating that surface antigens of the strains in these CCs experienced robust recombination.
Association of CC with virulence genes
Association of the three major CCs with virulence genes was observed in our study. Both the kpsMT K1 and ompT genes were positively associated with CC1687 and CC58, respectively. The sfa and iss genes were negatively associated with CC1687 and CC10, respectively.
Geographical location and tissue origin of ExPECs
The eighty-one ExPECs were isolated from ten provinces in china, mainly from Hubei (n = 30) and Hunan (n = 29). The strains in CC1698 concentrate in Hunan province. Tissue origin of the samples varied, including lung, liver, spleen, pericardial sac, brain, blood, lymph, joint, and hydrops. The majority of the tissues from which ExPECs were collected were lung (46.9%), blood (12.3%) and pericardial sac (11.1%), The statistical analysis showed that no significant correlation was found between geographical location and tissue of origin.
Extraintestinal E. coli can lead to a variety of diseases outside the intestinal tract in a broad range of hosts. To date, most research has focused on the relationship between ExPECs in birds with those in human [7–9]. However, ExPECs of swine origin have not been well documented and investigated. The importance of swine ExPECs is not fully appreciated at present. This is mainly because the phylogenetic types present in swine differ considerably with the majority of ExPEC isolates in human, and because There have been much lower number of ExPECs isolated from pork compared with poultry meat [10, 11].
Studies by Dezfulian et al.  and Marynard et al.  showed that greater than 50% of ExPECs from diseased pigs belonged to phylogenetic groups A or B1, whereas most ExPECs from human clinical cases belonged to phylogroups B2 and/or D. A survey of ExPECs in food products carried out by Johnson et al.  in Minneapolis–St. Paul retail establishment in 1999–2000 revealed that no ExPECs were isolated from pork. However, ExPECs from pork have been detected in subsequent reports. A survey by Johnson et al.  in the same area reported that ExPEC was isolated from b 18.9% of samples of raw beef and pork in 2001–2003. Xia et al.  reported in 2006 that ExPECs were isolated from 8.3% of pork chops in Georgia, Maryland, Oregon, and Tennessee. Data from our clinical microbiology laboratory also showed that porcine ExPEC diseases increased five times from 2004 to 2007 in the 19 province of China . These investigations provide an important clue that the occurrence of ExPEC in pork has increased gradually over recent years. This increasing trend may be caused by the over- or misuse of antibiotics and/or bacteria diseases infecting immunocompromised pigs [15–18]. Multidrug resistance of porcine ExPEC strains and immunosuppressive diseases may promote the incidence of porcine ExPEC diseases in China. Therefore, swine ExPECs may comprise a significant proportion of ExPECs in the retail food sold throughout China, again highlighting the importance of why their molecular epidemiology must be better understood.
Most researchers have found that ExPECs of human and avian origin belong to phylogenetic group B2 and, to a lesser extent, to phylogroup D. Sporadically, these ExPECs have appeared in phylogroups A and B1. The majority of phylogroup A and B1 isolates are mainly the commensal and intestinal pathogenic strains. Moreno et al.  reported that ExPECs of phylogroups A and B1 presented low virulent potential and could cause sporadic diseases in compromised patients only. In our study, porcine ExPECs were mainly distributed among phylogroups A, B1 and D, but a few strains were classified as phylogroup B2. The discrepancy in phylogenetic types between swine ExPECs and human ExPECs could be explained by host specificity.
MLST analysis is a powerful phylogenetic tool with high reproducibility, thus facilitating global comparisons. By analyzing multiple housekeeping gene loci, strains with related phylogenetic background can be grouped and analyzed together . Comparison of porcine ExPECs with the MLST database showed that most of the clones observed in our study possessed the same allele profiles as human ExPEC clones reported elsewhere, especially those from South America, North America and Europe. Isolates of ExPECs in CC10 were first detected in Sweden and later seen in other parts of Europe. These were found to caused urinary infections in humans (http://mlst.ucc.ie/mlst/dbs/Ecoli/GetTableInfo_html). Moreno et al.  reported that ExPEC strains in this group caused extraintestinal infections in compromised patients. Some of the STs, such as ST520, ST218 and ST746 in CC10, were reported exclusively as diarrheagenic E. coli rather than ExPEC diseases. ExPEC strains in these clonal groups may be derived from intestinal pathogenic E. coli strains and have evolved into ExPEC through accumulation of virulence factors . Additionally, the new ST1687 and its related strains, with homogeneous phylogenetic background and belonging to the same serogroup, are described here for the first time. The large number of the novel STs suggests that they may prevail in Hunan and possess advantageous traits that have evolved independently. The molecular and epidemiological characteristics of the ST1687 strains require further surveillance and investigation.
The correlation between CC and virulence factors was also investigated in this study. Significant differences were detected between the four main clonal groups with respect to patterns of virulence genes. Strains in CC10 had significantly fewer iss genes, while strains in CC58 had more ompT genes. The iss gene is a plasmid-associated gene that has been reported to be widely distributed among avian extraintestinal pathogenic E. coli (APEC) . The ompT gene is a protease that cleaves colicin A E1 E2 E3. However, the role of iss and ompT in the pathogenesis of APECs or porcine ExPECs is unclear. Strains in CC1648 had few sfa and ibeA genes but more kptM K1 gene. Sfa is a Dr-binding adhesion confirmed to be associated with E.coli causing neonatal meningitis and urinary infection. The ibeA gene and K1 gene are reported to be important in the pathogenesis of NMECs . The presence or absence of these genes may be a trait inherited by clonal lineage and attributed to special pathogenic mechanisms. These results were consistent with those in a previous study on ExPECs from humans and other animals, in which ExPECs exhibited miscellaneous virulence gene patterns . Porcine ExPECs in different CCs probably also infected pigs via different mechanisms as human ExPECs. We also observed that the ExPEC isolated from diseased pigs possessed different virulence gene patterns as those isolated from retail pork. Xia et al.  reported that afa/draBC was not detected in 200 ExPEC strains isolated from pork, but this gene was detected in the four main clonal groups of porcine ExPECs. Therefore, the source of ExPECs detected in retail pork may be due to contamination during processing, storage, or transportation.
The presence of ExPECs in food products suggests that ExPECs may represent a new class of foodborne pathogens. Although pork is a common meat source that is consumed by a large number of people, the zoonotic potential of porcine ExPECs has not been thoroughly investigated. Studies on ExPEC strains in pork or diseased pigs have been limited to a few areas and sample pools [9, 10, 24]. There is a need to expand such studies to more countries to obtain a better picture of the profile of porcine ExPECs and its link to human diseases. Our study provides preliminary data to understand the association of clonality with virulence gene traits in porcine ExPECs. We also find that the predominant serogroups of porcine ExPECs are serogroup O11, O8, O138, O161, O101, and that clonal complex 10 and clonal complex 1687 are the main porcine ExPEC clonal lineages in China. The contribution of Porcine ExPECs to human ExPEC diseases could be more convincing and better explained if there is more information about the clonality, phylogenetic group, and virulent trait of commensal E. coli and human ExPECs in China.
Material and methods
Strains isolation and characterization
ExPEC (n = 81) strains were isolated from the samples of lung, brain, blood, spleen, and Pericardial sac of pigs with septicemia, meningitis and respiratory disease by the Clinic Microbiology Laboratory of Huazhong Agricultural University from 2006–2007. These samples came from the southern and central provinces of China. Each E. coli from the samples was collected from different pigs and confirmed with a rapid identification system (Biolog MicroStation). Mutiplex PCR described by Johnson et al.  was used to amplify five virulence genes: papA and/or papC (P fimbriae); sfa/foc (S and F1C fimbriae); afa/dra (Dr-antigen-binding adhesins); iutA (aerobactin receptor); and kpsMT II (group 2 capsular polysaccharide units). E. coli that possessed at least two of the five virulence genes were thought to be ExPECs. Serogrouping of porcine ExPEC isolates was performed with an E. coli serogrouping kit (China Institute of Veterinary Drug Control, Beijing, China). Phylogenetic grouping was conducted using the method of Clermont et al. . Triplex PCR was applied to detect chuA, yjaA and EspE4.C2 genes in swine ExPECs. Based on the presence or absence of these three genes, all ExPEC isolates were divided into four main phylogenetic groups according to the following criteria: isolates of chuA– and TspE4.C2– were assigned to group A; chuA–/TspE4.C2+ isolates were assigned to group B1. chuA+/yjaA + Isolates were assigned to group B2; chuA+/yjaA– isolates were assigned to group D. There are 35 different virulence genes that were thought to be associated with human and/or avian extraintestinal diseases (Table 1) that were assessed in this study based on reports from Johnson et al. , Chapman et al.  and Ewers et al. . The amplified products were analyzed using 0.8% agarose gel electrophoresis, and recorded with a gel documentation system. Results for each virulence factor that conflicted initially were investigated further and sequenced for verification. The details of these procedures are described in our previous paper .
Seven housekeeping genes (adk, fumC, gyrB, icd, mdh, purA and recA) were amplified uing primers and protocols outlined on the E. coli MLST website (http://mlst.ucc.ie/mlst/dbs/Ecoli). The PCRs were carried out in a 50-μl mixture containing 1 μl of DNA, 1 μl of each primer, 5 μl of primer star buffer, 4 μl of dNTP mixture, 38.5 μl of deionized distilled water, and 0.5 μl of primer star DNA polymerase. The amplified PCR products were confirmed with 0.8% (w/v) agarose gels through electrophoresis, and then sequenced by Sangon Biological Engineering Technology Inc (Shanghai, China). The sequences were compared with those of the seven housekeeping genes contained in the E. coli MLST database to obtain allelic profiles and sequence type (ST).
The proportion of distributions of the strains in clonal groups according to their serogroups and virulence genes were evaluated by Fisher’s exact test (two-tailed). Comparisons between different groups were tested, and if a significant difference was observed, comparisons between each group and all other groups combined were evaluated. A P-value less than 0.01 was considered significant. All statistical analyses were performed with SPSS version 17.0 (SPSS, Chicago, IL, USA).
Russo TA, Johnson JR: Medical and economic impact of extraintestinal infections due to Escherichia coli: focus on an increasingly important endemic problem. Microbes Infect. 2003, 5: 449-456. 10.1016/S1286-4579(03)00049-2.
Johnson JR, Russo TA: Molecular epidemiology of extraintestinal pathogenic (uropathogenic) Escherichia coli. Int J Med Microbiol. 2005, 295: 383-404. 10.1016/j.ijmm.2005.07.005.
Picard B, Garcia JS, Gouriou S, Duriez P, Brahimi N, Bingen E, Elion J, Denamur E: The link between phylogeny and virulence in Escherichia coli extraintestinal infection. Infect Immun. 1999, 67: 546-553.
Johnson JR, Murray AC, Gajewski A, Sullivan M, Snippes P, Kuskowski MA, Smith KE: Isolation and molecular characterization of nalidixic acid-resistant extraintestinal pathogenic Escherichia coli from retail chicken products. Antimicrob Agents Chemother. 2003, 47: 2161-2168. 10.1128/AAC.47.7.2161-2168.2003.
Tang XB, Tan C, Zhang X, Zhao ZQ, Xia X, Wu B, Guo AZ, Zhou R, Chen HC: Antimicrobial resistances of extraintestinal pathogenic Escherichia coli isolates from swine in China. Microb Pathog. 2011, 50: 207-212. 10.1016/j.micpath.2011.01.004.
Tan C, Tang XB, Zhang X, Ding Y, Zhao ZQ, Wu B, Cai XW, Liu ZF, He QG, Chen HC: Serogroups, Virulence genes, and phylogeny of extraintestinal pathogenic Escherichia coli isolates from diseased swine in China. Vet J2011, Epub ahead of print.
Moulin-Schouleur M, Schouler C, Tailliez P, Kao MR, Brée A, Germon P, Oswald E, Mainil J, Blanco M, Blanco J: Common virulence factors and genetic relationships between O18:K1:H7 Escherichia coli isolates of human and avian origin. J Clin Microbiol. 2006, 44: 3484-3492. 10.1128/JCM.00548-06.
Tivendale KA, Logue CM, Kariyawasam S, Jordan D, Hussein A, Li G, Wannemuehler Y, Nolan LK: Avian-pathogenic Escherichia coli strains are similar to neonatal meningitis E. coli strains and are able to cause meningitis in the rat model of human disease. Infect Immun. 2010, 78: 3412-3419. 10.1128/IAI.00347-10.
Kobayashi RK, Aquino I, Ferreira AL, Vidotto MC: EcoR. Phylogenetic Analysis and Virulence. Genotyping of Avian Pathogenic Escherichia coli Strains and Escherichia coli Isolates from Commercial Chicken Carcasses in Southern Brazil. Foodborne Pathog Dis. 2011, 8: 631-634. 10.1089/fpd.2010.0726.
Dezfulian H, Batisson I, Fairbrother M, Lau K, Nassar A, Szatmari G, Harel J: Presence and characterization of extraintestinal pathogenic Escherichia coli virulence genes in F165-positive E. coli strains isolated from diseased calves and pigs. J Clin Microbiol. 2003, 41: 1375-1385. 10.1128/JCM.41.4.1375-1385.2003.
Maynard C, Bekal S, Sanschagrin F, Levesque C, Brousseau R, Masson L, Larivière S, Harel J: Heterogeneity among virulence and antimicrobial resistance gene profiles of Extraintestinal Escherichia coli isolates of animal and human origin. J Clin Microbiol. 2004, 42: 5444-5452. 10.1128/JCM.42.12.5444-5452.2004.
Johnson JR, Delavari P, O’Bryan TT, Smith KE, Tatini S: Contamination of retail foods, particularly turkey, from community markets (Minnesota, 1999–2000) with antimi-crobial-resistant and extraintestinal pathogenic Escherichia coli. Foodborne Pathog Dis. 2005, 2: 38-49. 10.1089/fpd.2005.2.38.
Johnson JR, Kuskowski MA, Smith K, O’Bryan TT, Tatini S: Antimicrobial-resistant and extraintestinal pathogenic Escherichia coli in retail foods. J Infect Dis. 2005, 191: 1040-1049. 10.1086/428451.
Xia X, Meng J, Zhao S, Bodeis-Jones S, Gaines SA, Ayers SL, McDermott PF: Identification and antimicrobial resistance of extraintestinal pathogenic Escherichia coli from retail meats. J Food Prot. 2011, 74: 38-44. 10.4315/0362-028X.JFP-10-251.
Cai X, Chen H, Blackall J, Yin Z, Wang L, Liu Z, Jin M: Serological characterization of Haemophilus parasuis isolates from China. Vet Microbiol. 2005, 111: 231-236. 10.1016/j.vetmic.2005.07.007.
Normile D: Infectious diseases. WHO probes deadliness of China's pig-borne disease. Science. 2005, 309: 1308-1309.
Zhang YN, Peng J, Wang Q, Pei ZF, Zhang WJ, Niu ZX: Appearance of blaCMY-2) gene-positive Salmonella isolates of pig origin in China. Int J Antimicrob Agents. 2008, 31: 292-293. 10.1016/j.ijantimicag.2007.10.023.
Zhao Z, Wang C, Xue Y, Tang X, Wu B, Cheng X, He Q, Chen H: The occurrence of Bordetella bronchiseptica in pigs with clinical respiratory disease. Vet J. 2011, 188: 337-340. 10.1016/j.tvjl.2010.05.022.
Moreno E, Prats G, Planells I, Planes AM, Pérez T, Andreu A: Characterization of Escherichia coli isolates derived from phylogenetic groups A and B1 causing extraintestinal infection. Enferm Infecc Microbiol Clin. 2006, 24: 483-489. 10.1157/13092463.
Feil EJ, Li BC, Aanensen DM, Hanage WP, Spratt BG: eBURST: Inferring Patterns of Evolutionary Descent among Clusters of Related Bacterial Genotypes from Multilocus Sequence Typing Data. J Bacteriol. 2004, 186: 1518-1530. 10.1128/JB.186.5.1518-1530.2004.
Brzuszkiewicz E, Brüggemann H, Liesegang H, Emmerth M, Olschläger T, Nagy G, Albermann K, Wagner C, Buchrieser C, Emody L, Gottschalk G, Hacker J, Dobrindt U: How to become a uropathogen: comparative genomic analysis of extraintestinal pathogenic Escherichia coli strains. Proc Natl Acad Sci U S A. 2006, 103: 12879-12884. 10.1073/pnas.0603038103.
Rodriguez-siek E, Giddings W, Doetkott C, Johnson J, Nolan K: Characterizing the APEC pathotype. Vet Res. 2005, 36: 241-256. 10.1051/vetres:2004057.
Maruvada R, Kim KS: IbeA and OmpA of Escherichia coli K1 Exploit Rac1 Activation for Invasion of Human Brain Microvascular Endothelial Cells. Infect Immun. 2012, 80: 2035-2041. 10.1128/IAI.06320-11.
Girardeau JP, Lalioui L, Said AM, De Champs C, Le Bouguénec C: Extended virulence genotype of pathogenic Escherichia coli isolates carrying the afa-8 operon: evidence of similarities between isolates from humans and animals with extraintestinal Infections. J Clin Microbiol. 2003, 41: 218-226. 10.1128/JCM.41.1.218-226.2003.
Clermont O, Bonacorsi S, Bingen E: Rapid and simple determination of the Escherichia coli phylogenetic group. Appl Environ Microbiol. 2000, 66: 4555-4558. 10.1128/AEM.66.10.4555-4558.2000.
Johnson JR, Stell AL: Extended virulence genotypes of Escherichia coli strains from patients with urosepsis in relation to phylogeny and host compromise. J Infect Dis. 2000, 181: 261-72. 10.1086/315217.
Chapman TA, Wu XY, Barchia I, Bettelheim KA, Driesen S, Trott D, Wilson M, Chin JJ: Comparison of virulence gene profiles of Escherichia coli strains isolated from healthy and diarrheic swine. Appl Environ Microbiol. 2006, 72: 4782-4795. 10.1128/AEM.02885-05.
Ewers C, Li G, Wilking H, Kiessling S, Alt K, Antáo EM, Laturnus C, Diehl I, Glodde S, Homeier T, Böhnke U, Steinrück H, Philipp HC, Wieler LH: Avian pathogenic, uropathogenic, and newborn meningitis-causing Escherichia coli: How closely related are they?. Int J Med Microbiol. 2007, 297: 163-176. 10.1016/j.ijmm.2007.01.003.
This work was supported by grants from the National Natural Science Foundation of China (NO. 31030065) and the 863 Program (2012AA101601).
The authors declare that they have no competing interests.
Author Yi Ding, Xibiao Tang, Bin Wu, Zhuofei Xu, Lu Ping, Wugang Liu and Ruixuan Zhang performed the experiments. Author Weicheng Bei, Huanchun Chen and Chen Tan designed the experiments. Yi Ding, and Zhuofei Xu performed the analysis and interpretation of data. Yi Ding and Chen Tan wrote the manuscript. All the author approve submission of this manuscript to BMC Veterinary Research.
Electronic supplementary material
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.