BACKGROUND
clostridium difficile is an anaerobic, gram-positive, spore-forming bacterium responsible for infections ranging from mild diarrhea to pseudomembranous colitis, primarily in elderly patients exposed to antibiotics.
recently, this bacterium has become one of the most frequent microorganisms responsible for healthcare-associated infections in the united states  <cit> .
during  <dig>  the first outbreak of a strain with “hyperproducer ” of toxin a and b, referred to as the north american pulsed-field type one , restriction-endonuclease analysis type bi, or pcr ribotype  <dig> , was reported in north america .
since then, cases have been reported worldwide.
nevertheless, large outbreaks of nap1/bi/ <dig> are less reported in asia and latin america compared to north america and europe  <cit> .
however, cases of nap1/bi/ <dig> infection have been already described in hong kong and guangzhou, china  <cit> .
nap1/bi/ <dig> is characterized by an in vitro overproduction of toxins a and b as well as a binary toxin production, and an 18 bp deletion in tcdc gene  <cit> .
various epidemiological differences are noted in asia; for example, the most common circulating ribotypes are  <dig>   <dig>  and  <dig>  <cit> .
moreover, the toxin a-negative, b-positive  strain  <dig> has received wide attention  <cit> .
to the best of our knowledge, no report of c. difficile infection  outbreaks in china has been published, and even few of case reported so far in the english literature.
in the present study, we report the first nap1/bi/ <dig> outbreak in a hospital in mainland china, traced by a rare case of bloodstream infection  using whole genome.
this process of identifying outbreak was dramatic by a comprehensive retrospective investigation using temporal and spatial analysis of the cases and strain characterization including pcr-ribotyping, multi-locus sequence typing  and whole genome sequencing .
the identification of the outbreak immediately provides four key implications: 1) the threat of cdi in china is more serious than that previously believed because of the under diagnosis of cdi in hospitals, and more cases were ignored; 2) the epidemiology of cdi should be carefully investigated to check whether there have been more cases of nap1/bi/ <dig> infection  <cit> ; and 3) the infection control strategy should be strengthened.
4) wgs method may be helpful for investigating outbreaks and tracking strain transmission  <cit> .
methods
case definition and cdi diagnosis
this study was conducted in the department of traditional chinese medicine respiratory, a hospital of china from march  <dig> to april  <dig>  a cdi case-patient was defined as a patient with  <dig> or more unformed stools in  <dig> or fewer consecutive hours and a stool test resulted positive for c. difficile toxin  <cit> .
stool samples were tested for toxin by enzyme immunoassay , and were anaerobically cultured on cycloserine-cefoxitin-fructose agar.
suspected colonies were identified by agglutination with c. difficile latex reagent for the somatic antigen .
retrospectively epidemic investigation
a comprehensive retrospective epidemiological investigation was conducted.
first, clinical information from medical records, including hospital stays, admitting diagnoses, antibiotic use, and occurrence of diarrhea were collected.
using these retrospective case searches, in total,  <dig> isolates were collected.
among  <dig> patients with diarrhea,  <dig> were diagnosed as cdi  <cit> , including one patient with a c. difficile bloodstream infection.
symptom onset occurred after 48 h of hospitalization in  <dig> patients while  <dig> patients presented with diarrhea when admission.
second, the architectural structure of the ward and room arrangement were collected by field epidemiology investigation; at the same time, clinicians, nurses, and nurse-assistant were questioned about their behaviors during healthcare, paying particular attention to the overlap of hospital stays, use of the toilet and bathroom, and cleaning procedures.
validation and typing the strains in two outbreaks
in house-pcr was used to detect tcda and tcdb, using the primer sets according to lemee’s study  <cit> , whereas mlst was used to allele-typing seven housekeeping genes, including adk, atpa, dxr, glya, reca, soda and tpi  <cit> .
the isolates were characterized by a capillary gel electrophoresis-based ribotyping method using the nap1/bi/ <dig> strain uk <dig> as the reference  <cit> .
genetic diversity, whole genome sequencing and bioinformatics analysis
a 500 bp library was constructed and whole genome sequencing was performed using illumina hiseq  <dig> platform  to obtain a 100 bp paired-end read, with  <dig> , <dig>  bps in average for each sample   <cit> .
the multi-alignment sequences from read data and references  were obtained by the method that combines alignments from mappings to multiple reference sequences via realphy  <cit> , and then the phylogenetic tree was reconstructed via bayesian evolutionary analysis by beast  <cit> .
the genomic single nucleotide polymorphisms  between nap1/bi/ <dig> strains were detected by mapping the c. difficile nap1/bi/ <dig> reference genome r <dig> using soap <dig> software  <cit> .
the sites with quality of < <dig> were removed.
these snps were finally used to reconstruct the phylogenetic tree via bayesian evolutionary analysis by beast  <cit> .
the mutation pattern of nap1/bi/ <dig> was clustered by a hierarchical cluster algorithm using cluster/treeview software  <cit> .
RESULTS
identification the bsi case of c. difficile nap1/bi/ <dig> infection
the patient was admitted on august  <dig>   <dig> for community-acquired pneumonia.
he was treated by cefepime  and levofloxacin  intravenously for 7 days.
on august  <dig>  he displayed diarrhea, which progressively worsened.
on august  <dig>  his symptoms included fever, cramping, and abdominal discomfort.
his stool sample was tested positive for the c. difficile toxin and culture.
we simultaneously obtained his c. difficile isolate from an anaerobic blood culture.
he was successfully treated with oral metronidazole  for 13 days and then sequentially intravenous metronidazole  for 4 days, and he was discharged on september  <dig>  he was re-admitted three days later with persistent fever and diarrhea.
stool sample was positive again for toxin and culture but the blood culture was negative.
he was treated with oral vancomycin  for 12 days, and discharged on october  <dig>  clinical information and laboratory testing results of this patient  have been described in additional file 1: table s <dig>
retrospective study of clinical information of patients with cdi in the ward
this unusual bsi case due to c. difficile triggered a retrospective investigation of all cdi cases in the ward.
in our record, nineteen additional patients with cdi were identified from march  <dig> to april  <dig>  their hospitalization periods are described in fig.
<dig>  all the patients were elderly  and had received antimicrobial treatments before diarrhea .
it is worth noting that cdi patients were more likely to be treated with third generation cephalosporins and fluoroquinolones than patients negative for cdi in the ward at the same time .
two peaks of cdi occurred: one from march to july  <dig> and the other from july  <dig> to april  <dig> suggesting two potential outbreaks .fig.
<dig> occurrence of cdi in this ward from march  <dig> to may  <dig>  the x-axis represents the time line and the y-axis represents the patient id.
the green bars represent the period of hospital stay of each patient, the purple and yellow lines represent the dates of bouts of diarrhea and c. difficile isolation, respectively.
the isolate 13b from blood is marked by a red star .
the number of new cases in each month is displayed using vertical bars at the top, red for nap1/bi/ <dig> cases, and gray for remaining cases
all the positive strains of c. difficile isolated from stool samples were retained in our hospital.
the clinical information and laboratory testing results have been presented in additional file 1: table s <dig>  seventeen strains were a + b+ and five were a − b+.
mlst characterized these isolates into five sequence types : st <dig>  st <dig>  st <dig>  st <dig>  and st <dig>  eleven patients were infected by st <dig> strain  .
only two patients, p <dig> and p <dig> infected by st <dig> and st <dig>  respectively, already presented cdi-caused diarrhea by the time of admission, which may be community-acquired.
thus, all st <dig> infections were acquired within the hospital.
furthermore, all the st <dig> strains were confirmed as ribotype  <dig> by capillary gel electrophoresis fingerprint of ribotyping .fig.
<dig> capillary gel electrophoresis fingerprint of ribotyping of all c. difficile strains isolated from cdi patients in this study.
each lane represents as following: lane m, dna marker ; lane uk <dig>  reference strain uk <dig> ; other lanes, the  <dig> isolates from  <dig> cdi patients
genetic population analysis in outbreaks
all the isolated c. difficile strains were cultured and sequenced .
after mapping all reads to the reference genomes of r <dig>  and m <dig> , the multi-mapping results were combined and  <dig> ,190 bps aligned sequence was obtained from each genome, and then a bayesian tree was reconstructed .
the results revealed that all isolates were clearly characterized into five groups, as suggested by mlst , and the majority isolates from the outbreaks belonged to st <dig>  moreover, more than thousands of snps were detected between st <dig> and other types ; therefore, st <dig>  st <dig>  st <dig> and st <dig> presented large genome diversity compared with published st <dig> strains.
these results confirmed that the strains belonging to different sts are not genetically related as which has been demonstrated before  <cit> .
we also compared the snps within the same sts.
interestingly, when we scanned the genomes of st <dig>  most of them were close and tightly clustered together .
we discovered only one snp difference between isolates  <dig> and  <dig>  illustrating their close relationship.
these results suggest that st <dig> isolates might come from the same origin and were from patient to patient transmissions.
in addition to st <dig>  isolates within the other sts also clustered together.
for example, two isolates belonging to st <dig>  presented a genetic distance with seven snps, and three isolates of st <dig> presented one or two snps.
these results suggest that isolates from the same st may also come from the same origin.
these results are in accordance with a former study in which 0– <dig> snvs were identified between transmitted isolates obtained less than 124 days apart  <cit> .fig.
<dig>
a phylogenetic tree of  <dig> isolates from this study and two reference strains r <dig>  and m <dig>  reconstructed by beast based on the multi-mapping results of realphy.
the strain id is coordinate to patient id.
b snp pattern and relationship of nap1/bi/ <dig> isolates.
thirteen isolates from  <dig> cases  are clustered using nap1/bi/ <dig> strain r <dig> as the outgroup.
the light blue squares represent mutations compared to r20291
by combining the clinical and laboratory information, we confirmed two outbreaks of nap1/bi/ <dig> had occurred, one from march to july  <dig> and the other from december  <dig> to february  <dig> .fig.
<dig> the temporal and spatial transmission map with three potential risks for c. difficile nap1/bi/ <dig> spreading.
the transmission connections in two consecutive years were marked using blue and red lines, respectively.
the full lines represent the potential transfer route with both genomics and epidemic evidence, and dashed lines connected the patients among whom only have potential genetic links.
”?” represented the unknown nodes during the transmission and evolution.
the architectural structure, room arrangement, and distribution of nap1/bi/ <dig> cases in this ward were displayed at the bottom.
patients during the two outbreaks are distinguished by different colors of background: red for those in  <dig> and blue for those in  <dig> to 2014
whole genome analysis of nap1/bi/ <dig> isolates
after mapping the reads of nap1/bi/ <dig>  isolates to the reference genome r <dig>  a total of  <dig> snps were identified in our nap1/bi/ <dig> strains, including t1476g and g1514a located in rpob gene associated with fluoroquinolone-resistance .
all sequences from p <dig> were tightly clustered together, suggesting that these occurrences of cdi  were related as we assumed .
the route of genetic transfer is clearly from 13a to 13b then to 13c, indicated by the genetic pattern with one snp for each strain, although strain 13a presented an abnormal mutation.
the results of wgs analysis confirmed our subsection that strain 13b causing bsi might have originated from the gut .
the isolate 13c, responsible for cdi recurrence, was nearest to 13b from his blood.
however, this hypothesis ignores the presence of within host diversity.
the other possible alternation is that 13a, 13b and 13c may have existed simultaneously in the gut at both time points, and 13b caused bsi at the initial time.
since we cannot obtain multiple colony picks in this retrospective study, we empirically estimate that the gut infection lead to secondary bsi and gut persistence explains the relapse.
construction of the transmission map of nap1/bi/ <dig> for this outbreak
to reconstruct a putative transmission map, the clinical epidemiologic data were combined with the genetic data from results of wgs .
first, the temporal and spatial transmission links were constructed.
the patients with overlapping time spent in bed were considered as the potential transmission connections .
second, by conducting the field epidemiological investigation and consulting with heath care workers and patients, we assumed that the toilet and cleaning procedures were potential at high risk of transmission.
indeed, only one commode in the toilet and one mop were available in the ward, increasing the possible risk of the spore transmission among the patients.
last, we constructed the potential connections between the patients by genetic relationship among these strains.
by combining all the data, particularly the genetic relationships between isolates and possible contacts between patients, we reconstructed the transmission routes of all nap1/bi/ <dig> cases to explain the outbreaks.
furthermore, our study indicated a trend between decreasing number of cdi and increasing distance to the janitor’s equipment room and toilet .
genome sequences provide the detailed molecular evidence to connect part of the transmission routes.
first, phylogenetic analysis showed all nap1/bi/ <dig> isolates in this ward belonged to the same clone which existed in the ward in three years and caused two outbreaks, and we provided clear clues for explaining these outbreaks.
second, the isolates in these potential transfer routes shared most of the  <dig> snps, resulting in the direct connection among them .
for example, p <dig> only contains one snp difference with p <dig>  and this result indicates that these two patients may be infected by one strain or the pathogen may transfer from one to another.
similar results appeared in p <dig>  p <dig> and p <dig>  p <dig> and p <dig>  finally, wgs data also validate that the strains in two outbreaks in consecutive years were only different by  <dig> stable snps, indicating the infected strain may exist in the hospital and transfer again in the following year .
based on the sequence diversity of our data set, the estimated time to the most recent common ancestor of the c. difficile strain was december  <dig> .
the potential connection of isolates according to time, space distribution, and molecular clues were integrated to define the full map of transmission routes .
discussion
to our knowledge this is the first study investigating outbreaks of healthcare-associated c. difficile nap1/bi/ <dig> strain infection in china and reconstructing the transmission map by combining retrospective epidemiological analysis with whole-genome sequencing.
the outbreak caused by this strain has not been reported in china, especially via bloodstream infection, which has high mortality.
we provide evidence that wgs method could be used to tracing the source of infection and predicting the transmission routes.
c. difficile was the most commonly reported pathogen, accounting for  <dig>  % of healthcare-associated infections  <cit> .
the increased severity of cdi is primarily attributed to the emergence and rapid spread of the fluoroquinolone-resistant nap1/bi/ <dig> strain, which remains the most common clone worldwide  <cit> .
since the first large outbreak described in north america in  <dig>  the nap1/bi/ <dig> has become endemic in many hospital from north america, europe, and more recently in australia  <cit> .
while, for a long time, lacking of the reports of cdi caused by the strain, makes us confusing about the states of  <dig> in china.
three major reasons may have contributed to this situation: 1) the major diagnostic method for cdi available in chinese clinical laboratories is eia for toxins a/b, and very few hospitals can perform c. difficile culture, resulting in many undetected cdi cases; 2) the lack of anaerobic laboratory equipment and limited cdi knowledge by healthcare staff also contributes to the poor isolation compliance to infection control measures and reporting under diagnosis of cdi in china; and 3) there is no reference center for typing c. difficile strains in china and identification of the different pcr ribotypes requires the collection of reference strains.
therefore, distribution of c. difficile circulating clones including the nap/bi/ <dig> is not known in china.
pcr based sequencing method can detect the pathogen without reference strain by comparing with the public database.
thus, it becomes one of replacement to other method in lots of bacteria.
however, not like other hospital acquired infection  bacteria, such as mycobacterium tuberculosis and klebsiella pneumonia  <cit> , c. difficile is characterized by large genome diversity and instability increasing the difficulties to select proper primers for detection.
this characteristic will increase the difficulty to extend the existed molecular typing methods used in us or uk to develop country, and easy to cause clinic false-negative results.
so, the fact may be that: nap1/bi/ <dig> has existed in china before, yet it could not be identified immediately.
whole-genome sequencing can easily distinguish minor differences, providing valuable information for molecular typing, and pathogen or transmission route tracing  <cit> .
thus, the genome complication would boost the application of wgs in c. difficile research.
with the help of wgs, our research will enrich hai pathogen database, and hai outbreak tracing clues.
also, our research revealed the potential nap/bi/ <dig> threat for china will improve the completion of c. difficile nap/bi/ <dig> global epidemiology distribution  <cit> .
in this study, we investigated cdi outbreaks using medical personnel interview, medical record review and wgs for strain typing.
from the spatial analysis of the cases, we assumed that cdi patients are mainly concentrated near to the toilets, bathroom, and janitor’s equipment room , which may be suggested to improve infection control measures in the following health care.
however, we failed to obtain the isolates in environment.
we also conducted special consultation with nurse-assistant and determined that the chlorine-containing disinfectant concentration  used in the ward was lower than the concentration recommended for cdi infection control .
the epidemiological investigation revealed the absence of a dedicated area for cdi patient isolation, which was probably the most important contributor to the outbreak.
p <dig>  the first patient identified during the 2013– <dig> outbreak, was placed in different hospital beds in the days following the emergence of diarrhea.
in addition, this patient was neither isolated in a single room nor placed in contact precautions.
to the best of our knowledge, person-to-person transmission was the main route of cdi transmission.
however, the outbreak occurred during a long period spanning from  <dig> to  <dig>  suggesting the existence of a potential environmental reservoir.
for example, the hospitalized periods of patient four and seven had no overlap, and their rooms were quite apart away from each other , showing that there was no obvious epidemiological connection by conventional methods.
however, the wgs analysis revealed that there was only one snp difference between isolates  <dig> and  <dig> from these two patients , which provides an important epidemiology clue and help to form epidemiology hypothesis.
so, the method of wgs shows obvious potential superior to routine typing methods in epidemiology tracing.
for example, we may confirm the transfer route of p <dig> and the clear route of p <dig>  p <dig>  and p <dig>  which occurred during in the 2013– <dig> outbreak.
furthermore, isolated strains from most cases in the outbreaks were sequenced belonging to st <dig>  and there were also several strains belonging to other sts with sufficient genetic diversity to st <dig> in this study.
these genetically distinct isolates may indicate diverse alternative sources or more complicated transmission modes, since c. difficile is an ancient organism with diverse strains present in humans, animals, and food.
it is also noted that asymptomatic carriers might not be detected, which would have contributed to the transmission.
the literature review revealed that approximately 50 % of the molecular resulting clusters can be confirmed by contact tracing documenting transmission links, particularly for a longitudinal time, space, and population span, which are difficult to establish using classic epidemiological links  <cit> .
therefore, the molecular epidemiology based on wgs should be recommended for hai investigations, in addition to conventional field epidemiology investigations.
since we never found nap1/bi/ <dig> strain in the hospital before, this study is not carried out under a standard surveillance in the hospital, especially lack of pathogen information from environmental monitoring.
all clinical trials were drawn back to collect after we identified and validated the strain, the complete transfer route of all strains in hospital is still not very clear.
however, to our knowledge, this is the first case to identify c. difficile nap/bi/ <dig> outbreak by wgs method, which can be extended to other bacteria.
one of the most interesting features of the epidemic nap1/bi/ <dig> is high level of antibiotic resistance  <cit> .
in the present study, we found the outbreak strain in china is consistent but different to that reported in other countries .
for example, the rpob gene presents the special mutations in t1476g and g1514a compared with r <dig>  in addition to the previously fluoroquinolone-resistance site reported  <cit> .
in our cases, most patients were from the local community and were mainly admitted for chronic obstructive pulmonary disease and chronic bronchitis.
the misuse of antibiotics, including the overuse and inadequate use of antibiotics, actually increased the antibiotic pressure and risk for cdi.
in china, before  <dig>  patients suffering these diseases were usually prescribed antibiotics, particularly cephalosporins and fluoroquinolones.
additional file 1: figure s <dig> shows that cdi patients were more likely to be treated with third generation cephalosporins and fluoroquinolones than patients negative for cdi in the ward at the same time.
the results are in accordance with the previous publications  <cit> .
therefore, these observations highlight at least two noteworthy points: 1) primary medical clinicians should be aware of the potential cdi threat before the use of cephalosporins and fluoroquinolones to treat chronic obstructive pulmonary disease and chronic bronchitis of the elderly population and 2) the related family as well as healthcare workers caring for these elderly patients may also have potential c. difficile exposure risk, which should be particularly addressed.
as we mentioned above, our cases emerged after the strictest antibiotic regulation has been implemented by the chinese government in  <dig>  so, the regulations should be further enforced  <cit> .
in addition, the existence and outbreak of c. difficile in our research may underscore the importance of monitoring antibiotic proper usage and taking preventive action to reduce transmission of such multidrug-resistant organisms in the hospital setting, thereby preventing nosocomial outbreaks.
thus, we propose that c. difficile surveillance must be conducted in health cares in china, asia and other developing country, as an early warning of antibiotic resistance, allowing time for related detailed antibiotic interventions to be designed and implemented.
CONCLUSIONS
by genome sequencing, spatio-temporal analysis and field epidemiology investigation of the c. difficile nap1/bi/ <dig> outbreaks in  <dig> to  <dig> in this ward, we can estimate their complex transform network, revealed the possible modes of transmission.
based on their genetic diversity, we can assume that the toilets, bathroom, and janitor’s equipment room may be contaminated area, which may be suggested to improve infection control measures in the following health care.
additional file
additional file 1: supplementary figures and tables.
figure s <dig>  comparison of antibiotics use density  from c. difficile nap1/bi/ <dig> patients with that from all patients in this ward.
figure s <dig>  the phylogenetic relationship and divergence time of c. difficile nap1/bi/ <dig> isolates calculated by beast.
table s <dig>  epidemiology of cdi patients and molecular typing of c. difficile isolates.
table s <dig>  genome sequencing status of all isolates.
table s <dig>  the number of snps between different sts.
table s <dig>  snps of c. difficile nap1/bi/ <dig> isolates in this study compared to r <dig>
abbreviations
bsibloodstream infection
cdiclostridium difficile infection
haihospital acquired infection
mlstmulti-locus sequence typing
nap1the north american pulsed-field type one
snpsingle nucleotide polymorphism
stsequence type
wgswhole genome sequencing
funding
this work was supported by the national high technology research and development program , the national key technology r&d program , and the national natural science foundation of china .
availability of data and materials
the sequencing data described in this study have been deposited in the genbank database , and the sequencing methods and background information of isolates have been described in this article.
authors’ contributions
cc, yc, yc and hj designed the project.
cc, yc, hj, pd, hy, yz, fb and gg wrote the manuscript.
hj, hy, jw, and gh performed the strain isolation and molecular typing.
yc, zy, jz and qd performed the epidemiological study.
yz, yq performed the whole genome sequencing.
pd, wz, hw, nh, zw performed the bioinformatics analysis.
all authors have read and approved the final manuscript.
authors’ information
dr. hongbing jia is the chief of division of microbiology, department of clinical laboratory, china-japan friendship hospital.
his primary research interests are in clinical microbiology.
his group is devoting to special pathogens detection, bacterial resistant mechanisms, bacterial epidemiological research, especially focus on clostridium difficile.
competing interests
the authors declare that they have no competing interests.
consent for publication
not applicable.
ethics approval and consent to participate
this study was approved by the ethics committee of china-japan friendship hospital.

 clostridium difficile is an anaerobic, gram-positive, spore-forming bacterium responsible for infections ranging from mild diarrhea to pseudomembranous colitis, primarily in elderly patients exposed to antibiotics. recently, this bacterium has become one of the most frequent microorganisms responsible for healthcare-associated infections in the united states. during the first outbreak of a strain with “hyperproducer ” of toxin a and b, referred to as the north american pulsed-field type one, restriction-endonuclease analysis type bi, or pcr ribotype, was reported in north america. since then, cases have been reported worldwide. nevertheless, large outbreaks of nap1/bi/ are less reported in asia and latin america compared to north america and europe. however, cases of nap1/bi/ infection have been already described in hong kong and guangzhou, china. nap1/bi/ is characterized by an in vitro overproduction of toxins a and b as well as a binary toxin production, and an 18 bp deletion in tcdc gene. various epidemiological differences are noted in asia; for example, the most common circulating ribotypes are and. moreover, the toxin a-negative, b-positive strain has received wide attention. to the best of our knowledge, no report of c. difficile infection outbreaks in china has been published, and even few of case reported so far in the english literature. in the present study, we report the first nap1/bi/ outbreak in a hospital in mainland china, traced by a rare case of bloodstream infection using whole genome. this process of identifying outbreak was dramatic by a comprehensive retrospective investigation using temporal and spatial analysis of the cases and strain characterization including pcr-ribotyping, multi-locus sequence typing and whole genome sequencing. the identification of the outbreak immediately provides four key implications: 1) the threat of cdi in china is more serious than that previously believed because of the under diagnosis of cdi in hospitals, and more cases were ignored; 2) the epidemiology of cdi should be carefully investigated to check whether there have been more cases of nap1/bi/ infection; and 3) the infection control strategy should be strengthened. 4) wgs method may be helpful for investigating outbreaks and tracking strain transmission. methods case definition and cdi diagnosis this study was conducted in the department of traditional chinese medicine respiratory, a hospital of china from march to april a cdi case-patient was defined as a patient with or more unformed stools in or fewer consecutive hours and a stool test resulted positive for c. difficile toxin. stool samples were tested for toxin by enzyme immunoassay, and were anaerobically cultured on cycloserine-cefoxitin-fructose agar. suspected colonies were identified by agglutination with c. difficile latex reagent for the somatic antigen. retrospectively epidemic investigation a comprehensive retrospective epidemiological investigation was conducted. first, clinical information from medical records, including hospital stays, admitting diagnoses, antibiotic use, and occurrence of diarrhea were collected. using these retrospective case searches, in total, isolates were collected. among patients with diarrhea, were diagnosed as cdi, including one patient with a c. difficile bloodstream infection. symptom onset occurred after 48 h of hospitalization in patients while patients presented with diarrhea when admission. second, the architectural structure of the ward and room arrangement were collected by field epidemiology investigation; at the same time, clinicians, nurses, and nurse-assistant were questioned about their behaviors during healthcare, paying particular attention to the overlap of hospital stays, use of the toilet and bathroom, and cleaning procedures. validation and typing the strains in two outbreaks in house-pcr was used to detect tcda and tcdb, using the primer sets according to lemee’s study, whereas mlst was used to allele-typing seven housekeeping genes, including adk, atpa, dxr, glya, reca, soda and tpi. the isolates were characterized by a capillary gel electrophoresis-based ribotyping method using the nap1/bi/ strain uk as the reference. genetic diversity, whole genome sequencing and bioinformatics analysis a 500 bp library was constructed and whole genome sequencing was performed using illumina hiseq platform to obtain a 100 bp paired-end read, with,  bps in average for each sample. the multi-alignment sequences from read data and references were obtained by the method that combines alignments from mappings to multiple reference sequences via realphy, and then the phylogenetic tree was reconstructed via bayesian evolutionary analysis by beast. the genomic single nucleotide polymorphisms between nap1/bi/ strains were detected by mapping the c. difficile nap1/bi/ reference genome r using soap software. the sites with quality of < were removed. these snps were finally used to reconstruct the phylogenetic tree via bayesian evolutionary analysis by beast. the mutation pattern of nap1/bi/ was clustered by a hierarchical cluster algorithm using cluster/treeview software.  identification the bsi case of c. difficile nap1/bi/ infection the patient was admitted on august for community-acquired pneumonia. he was treated by cefepime and levofloxacin intravenously for 7 days. on august he displayed diarrhea, which progressively worsened. on august his symptoms included fever, cramping, and abdominal discomfort. his stool sample was tested positive for the c. difficile toxin and culture. we simultaneously obtained his c. difficile isolate from an anaerobic blood culture. he was successfully treated with oral metronidazole for 13 days and then sequentially intravenous metronidazole for 4 days, and he was discharged on september he was re-admitted three days later with persistent fever and diarrhea. stool sample was positive again for toxin and culture but the blood culture was negative. he was treated with oral vancomycin for 12 days, and discharged on october clinical information and laboratory testing results of this patient have been described in additional file 1: table s retrospective study of clinical information of patients with cdi in the ward this unusual bsi case due to c. difficile triggered a retrospective investigation of all cdi cases in the ward. in our record, nineteen additional patients with cdi were identified from march to april their hospitalization periods are described in fig.  all the patients were elderly and had received antimicrobial treatments before diarrhea. it is worth noting that cdi patients were more likely to be treated with third generation cephalosporins and fluoroquinolones than patients negative for cdi in the ward at the same time. two peaks of cdi occurred: one from march to july and the other from july to april suggesting two potential outbreaks.fig. occurrence of cdi in this ward from march to may the x-axis represents the time line and the y-axis represents the patient id. the green bars represent the period of hospital stay of each patient, the purple and yellow lines represent the dates of bouts of diarrhea and c. difficile isolation, respectively. the isolate 13b from blood is marked by a red star. the number of new cases in each month is displayed using vertical bars at the top, red for nap1/bi/ cases, and gray for remaining cases all the positive strains of c. difficile isolated from stool samples were retained in our hospital. the clinical information and laboratory testing results have been presented in additional file 1: table s seventeen strains were a + b+ and five were a − b+. mlst characterized these isolates into five sequence types: st st st st and st eleven patients were infected by st strain. only two patients, p and p infected by st and st respectively, already presented cdi-caused diarrhea by the time of admission, which may be community-acquired. thus, all st infections were acquired within the hospital. furthermore, all the st strains were confirmed as ribotype by capillary gel electrophoresis fingerprint of ribotyping.fig. capillary gel electrophoresis fingerprint of ribotyping of all c. difficile strains isolated from cdi patients in this study. each lane represents as following: lane m, dna marker; lane uk reference strain uk; other lanes, the isolates from cdi patients genetic population analysis in outbreaks all the isolated c. difficile strains were cultured and sequenced. after mapping all reads to the reference genomes of r and m, the multi-mapping results were combined and,190 bps aligned sequence was obtained from each genome, and then a bayesian tree was reconstructed. the results revealed that all isolates were clearly characterized into five groups, as suggested by mlst, and the majority isolates from the outbreaks belonged to st moreover, more than thousands of snps were detected between st and other types; therefore, st st st and st presented large genome diversity compared with published st strains. these results confirmed that the strains belonging to different sts are not genetically related as which has been demonstrated before. we also compared the snps within the same sts. interestingly, when we scanned the genomes of st most of them were close and tightly clustered together. we discovered only one snp difference between isolates and illustrating their close relationship. these results suggest that st isolates might come from the same origin and were from patient to patient transmissions. in addition to st isolates within the other sts also clustered together. for example, two isolates belonging to st presented a genetic distance with seven snps, and three isolates of st presented one or two snps. these results suggest that isolates from the same st may also come from the same origin. these results are in accordance with a former study in which 0– snvs were identified between transmitted isolates obtained less than 124 days apart.fig. a phylogenetic tree of isolates from this study and two reference strains r and m reconstructed by beast based on the multi-mapping results of realphy. the strain id is coordinate to patient id. b snp pattern and relationship of nap1/bi/ isolates. thirteen isolates from cases are clustered using nap1/bi/ strain r as the outgroup. the light blue squares represent mutations compared to r20291 by combining the clinical and laboratory information, we confirmed two outbreaks of nap1/bi/ had occurred, one from march to july and the other from december to february.fig. the temporal and spatial transmission map with three potential risks for c. difficile nap1/bi/ spreading. the transmission connections in two consecutive years were marked using blue and red lines, respectively. the full lines represent the potential transfer route with both genomics and epidemic evidence, and dashed lines connected the patients among whom only have potential genetic links. ”?” represented the unknown nodes during the transmission and evolution. the architectural structure, room arrangement, and distribution of nap1/bi/ cases in this ward were displayed at the bottom. patients during the two outbreaks are distinguished by different colors of background: red for those in and blue for those in to 2014 whole genome analysis of nap1/bi/ isolates after mapping the reads of nap1/bi/ isolates to the reference genome r a total of snps were identified in our nap1/bi/ strains, including t1476g and g1514a located in rpob gene associated with fluoroquinolone-resistance. all sequences from p were tightly clustered together, suggesting that these occurrences of cdi were related as we assumed. the route of genetic transfer is clearly from 13a to 13b then to 13c, indicated by the genetic pattern with one snp for each strain, although strain 13a presented an abnormal mutation. the results of wgs analysis confirmed our subsection that strain 13b causing bsi might have originated from the gut. the isolate 13c, responsible for cdi recurrence, was nearest to 13b from his blood. however, this hypothesis ignores the presence of within host diversity. the other possible alternation is that 13a, 13b and 13c may have existed simultaneously in the gut at both time points, and 13b caused bsi at the initial time. since we cannot obtain multiple colony picks in this retrospective study, we empirically estimate that the gut infection lead to secondary bsi and gut persistence explains the relapse. construction of the transmission map of nap1/bi/ for this outbreak to reconstruct a putative transmission map, the clinical epidemiologic data were combined with the genetic data from results of wgs. first, the temporal and spatial transmission links were constructed. the patients with overlapping time spent in bed were considered as the potential transmission connections. second, by conducting the field epidemiological investigation and consulting with heath care workers and patients, we assumed that the toilet and cleaning procedures were potential at high risk of transmission. indeed, only one commode in the toilet and one mop were available in the ward, increasing the possible risk of the spore transmission among the patients. last, we constructed the potential connections between the patients by genetic relationship among these strains. by combining all the data, particularly the genetic relationships between isolates and possible contacts between patients, we reconstructed the transmission routes of all nap1/bi/ cases to explain the outbreaks. furthermore, our study indicated a trend between decreasing number of cdi and increasing distance to the janitor’s equipment room and toilet. genome sequences provide the detailed molecular evidence to connect part of the transmission routes. first, phylogenetic analysis showed all nap1/bi/ isolates in this ward belonged to the same clone which existed in the ward in three years and caused two outbreaks, and we provided clear clues for explaining these outbreaks. second, the isolates in these potential transfer routes shared most of the snps, resulting in the direct connection among them. for example, p only contains one snp difference with p and this result indicates that these two patients may be infected by one strain or the pathogen may transfer from one to another. similar results appeared in p p and p p and p finally, wgs data also validate that the strains in two outbreaks in consecutive years were only different by stable snps, indicating the infected strain may exist in the hospital and transfer again in the following year. based on the sequence diversity of our data set, the estimated time to the most recent common ancestor of the c. difficile strain was december. the potential connection of isolates according to time, space distribution, and molecular clues were integrated to define the full map of transmission routes. discussion to our knowledge this is the first study investigating outbreaks of healthcare-associated c. difficile nap1/bi/ strain infection in china and reconstructing the transmission map by combining retrospective epidemiological analysis with whole-genome sequencing. the outbreak caused by this strain has not been reported in china, especially via bloodstream infection, which has high mortality. we provide evidence that wgs method could be used to tracing the source of infection and predicting the transmission routes. c. difficile was the most commonly reported pathogen, accounting for  % of healthcare-associated infections. the increased severity of cdi is primarily attributed to the emergence and rapid spread of the fluoroquinolone-resistant nap1/bi/ strain, which remains the most common clone worldwide. since the first large outbreak described in north america in the nap1/bi/ has become endemic in many hospital from north america, europe, and more recently in australia. while, for a long time, lacking of the reports of cdi caused by the strain, makes us confusing about the states of in china. three major reasons may have contributed to this situation: 1) the major diagnostic method for cdi available in chinese clinical laboratories is eia for toxins a/b, and very few hospitals can perform c. difficile culture, resulting in many undetected cdi cases; 2) the lack of anaerobic laboratory equipment and limited cdi knowledge by healthcare staff also contributes to the poor isolation compliance to infection control measures and reporting under diagnosis of cdi in china; and 3) there is no reference center for typing c. difficile strains in china and identification of the different pcr ribotypes requires the collection of reference strains. therefore, distribution of c. difficile circulating clones including the nap/bi/ is not known in china. pcr based sequencing method can detect the pathogen without reference strain by comparing with the public database. thus, it becomes one of replacement to other method in lots of bacteria. however, not like other hospital acquired infection bacteria, such as mycobacterium tuberculosis and klebsiella pneumonia, c. difficile is characterized by large genome diversity and instability increasing the difficulties to select proper primers for detection. this characteristic will increase the difficulty to extend the existed molecular typing methods used in us or uk to develop country, and easy to cause clinic false-negative results. so, the fact may be that: nap1/bi/ has existed in china before, yet it could not be identified immediately. whole-genome sequencing can easily distinguish minor differences, providing valuable information for molecular typing, and pathogen or transmission route tracing. thus, the genome complication would boost the application of wgs in c. difficile research. with the help of wgs, our research will enrich hai pathogen database, and hai outbreak tracing clues. also, our research revealed the potential nap/bi/ threat for china will improve the completion of c. difficile nap/bi/ global epidemiology distribution. in this study, we investigated cdi outbreaks using medical personnel interview, medical record review and wgs for strain typing. from the spatial analysis of the cases, we assumed that cdi patients are mainly concentrated near to the toilets, bathroom, and janitor’s equipment room, which may be suggested to improve infection control measures in the following health care. however, we failed to obtain the isolates in environment. we also conducted special consultation with nurse-assistant and determined that the chlorine-containing disinfectant concentration used in the ward was lower than the concentration recommended for cdi infection control. the epidemiological investigation revealed the absence of a dedicated area for cdi patient isolation, which was probably the most important contributor to the outbreak. p the first patient identified during the 2013– outbreak, was placed in different hospital beds in the days following the emergence of diarrhea. in addition, this patient was neither isolated in a single room nor placed in contact precautions. to the best of our knowledge, person-to-person transmission was the main route of cdi transmission. however, the outbreak occurred during a long period spanning from to suggesting the existence of a potential environmental reservoir. for example, the hospitalized periods of patient four and seven had no overlap, and their rooms were quite apart away from each other, showing that there was no obvious epidemiological connection by conventional methods. however, the wgs analysis revealed that there was only one snp difference between isolates and from these two patients, which provides an important epidemiology clue and help to form epidemiology hypothesis. so, the method of wgs shows obvious potential superior to routine typing methods in epidemiology tracing. for example, we may confirm the transfer route of p and the clear route of p p and p which occurred during in the 2013– outbreak. furthermore, isolated strains from most cases in the outbreaks were sequenced belonging to st and there were also several strains belonging to other sts with sufficient genetic diversity to st in this study. these genetically distinct isolates may indicate diverse alternative sources or more complicated transmission modes, since c. difficile is an ancient organism with diverse strains present in humans, animals, and food. it is also noted that asymptomatic carriers might not be detected, which would have contributed to the transmission. the literature review revealed that approximately 50 % of the molecular resulting clusters can be confirmed by contact tracing documenting transmission links, particularly for a longitudinal time, space, and population span, which are difficult to establish using classic epidemiological links. therefore, the molecular epidemiology based on wgs should be recommended for hai investigations, in addition to conventional field epidemiology investigations. since we never found nap1/bi/ strain in the hospital before, this study is not carried out under a standard surveillance in the hospital, especially lack of pathogen information from environmental monitoring. all clinical trials were drawn back to collect after we identified and validated the strain, the complete transfer route of all strains in hospital is still not very clear. however, to our knowledge, this is the first case to identify c. difficile nap/bi/ outbreak by wgs method, which can be extended to other bacteria. one of the most interesting features of the epidemic nap1/bi/ is high level of antibiotic resistance. in the present study, we found the outbreak strain in china is consistent but different to that reported in other countries. for example, the rpob gene presents the special mutations in t1476g and g1514a compared with r in addition to the previously fluoroquinolone-resistance site reported. in our cases, most patients were from the local community and were mainly admitted for chronic obstructive pulmonary disease and chronic bronchitis. the misuse of antibiotics, including the overuse and inadequate use of antibiotics, actually increased the antibiotic pressure and risk for cdi. in china, before patients suffering these diseases were usually prescribed antibiotics, particularly cephalosporins and fluoroquinolones. additional file 1: figure s shows that cdi patients were more likely to be treated with third generation cephalosporins and fluoroquinolones than patients negative for cdi in the ward at the same time. the results are in accordance with the previous publications. therefore, these observations highlight at least two noteworthy points: 1) primary medical clinicians should be aware of the potential cdi threat before the use of cephalosporins and fluoroquinolones to treat chronic obstructive pulmonary disease and chronic bronchitis of the elderly population and 2) the related family as well as healthcare workers caring for these elderly patients may also have potential c. difficile exposure risk, which should be particularly addressed. as we mentioned above, our cases emerged after the strictest antibiotic regulation has been implemented by the chinese government in so, the regulations should be further enforced. in addition, the existence and outbreak of c. difficile in our research may underscore the importance of monitoring antibiotic proper usage and taking preventive action to reduce transmission of such multidrug-resistant organisms in the hospital setting, thereby preventing nosocomial outbreaks. thus, we propose that c. difficile surveillance must be conducted in health cares in china, asia and other developing country, as an early warning of antibiotic resistance, allowing time for related detailed antibiotic interventions to be designed and implemented.  by genome sequencing, spatio-temporal analysis and field epidemiology investigation of the c. difficile nap1/bi/ outbreaks in to in this ward, we can estimate their complex transform network, revealed the possible modes of transmission. based on their genetic diversity, we can assume that the toilets, bathroom, and janitor’s equipment room may be contaminated area, which may be suggested to improve infection control measures in the following health care. additional file additional file 1: supplementary figures and tables. figure s comparison of antibiotics use density from c. difficile nap1/bi/ patients with that from all patients in this ward. figure s the phylogenetic relationship and divergence time of c. difficile nap1/bi/ isolates calculated by beast. table s epidemiology of cdi patients and molecular typing of c. difficile isolates. table s genome sequencing status of all isolates. table s the number of snps between different sts. table s snps of c. difficile nap1/bi/ isolates in this study compared to r abbreviations bsibloodstream infection cdiclostridium difficile infection haihospital acquired infection mlstmulti-locus sequence typing nap1the north american pulsed-field type one snpsingle nucleotide polymorphism stsequence type wgswhole genome sequencing funding this work was supported by the national high technology research and development program, the national key technology r&d program, and the national natural science foundation of china. availability of data and materials the sequencing data described in this study have been deposited in the genbank database, and the sequencing methods and background information of isolates have been described in this article. authors’ contributions cc, yc, yc and hj designed the project. cc, yc, hj, pd, hy, yz, fb and gg wrote the manuscript. hj, hy, jw, and gh performed the strain isolation and molecular typing. yc, zy, jz and qd performed the epidemiological study. yz, yq performed the whole genome sequencing. pd, wz, hw, nh, zw performed the bioinformatics analysis. all authors have read and approved the final manuscript. authors’ information dr. hongbing jia is the chief of division of microbiology, department of clinical laboratory, china-japan friendship hospital. his primary research interests are in clinical microbiology. his group is devoting to special pathogens detection, bacterial resistant mechanisms, bacterial epidemiological research, especially focus on clostridium difficile. competing interests the authors declare that they have no competing interests. consent for publication not applicable. ethics approval and consent to participate this study was approved by the ethics committee of china-japan friendship hospital.
by genome sequencing, spatio-temporal analysis and field epidemiology investigation, we can estimate their complex transform network and reveal the possible modes of transmission in this ward. the rapid spread of clostridium difficile nap1/bi/ has become one of the leading threats of healthcare-associated infections worldwide. integrated epidemiological investigation and whole genome sequencing of all strains, we constructed a nosocomial transmission map of these two c. difficile outbreaks and traced the origin of the infection. based on their genetic diversity, we can assume that the toilets, bathroom, and janitor’s equipment room may be contaminated area, which may be suggested to improve infection control measures in the following health care. surprisingly, we found the genome sequences were similar to c. difficile strain r indicating the occurrence of a rare c. difficile strain in china. through using whole genome sequencing, we succeeded to trace the origin of the bsi strain. clostridium difficilenap1/bi/027outbreaknosocomial infectionwhole genome sequencinghttp://dx.doi.org/ /501100002855ministry of science and technology of the people's republic of china2011aa02a1092012bah24f05jia hongbing http://dx.doi.org/ /501100001809national natural science foundation of china8120132281301402zhang wen chen chen issue-copyright-statement© the author 2016
 the rapid spread of clostridium difficile nap1/bi/ has become one of the leading threats of healthcare-associated infections worldwide. however, c. difficile infections have been rarely reported in asia, particularly in china.  in this study, we identified a rare c. difficile bloodstream infection from three isolates of a patient during repeated hospital admission. this finding triggered a retrospective epidemiological study to scan all cases and strains emerged from this ward during the past three years. using medical personnel interviews, medical record reviews and the genomic epidemiology, two outbreaks in and 2013– were identified. through using whole genome sequencing, we succeeded to trace the origin of the bsi strain. surprisingly, we found the genome sequences were similar to c. difficile strain r indicating the occurrence of a rare c. difficile strain in china. integrated epidemiological investigation and whole genome sequencing of all strains, we constructed a nosocomial transmission map of these two c. difficile outbreaks and traced the origin of the infection.  by genome sequencing, spatio-temporal analysis and field epidemiology investigation, we can estimate their complex transform network and reveal the possible modes of transmission in this ward. based on their genetic diversity, we can assume that the toilets, bathroom, and janitor’s equipment room may be contaminated area, which may be suggested to improve infection control measures in the following health care. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords clostridium difficilenap1/bi/027outbreaknosocomial infectionwhole genome sequencinghttp://dx.doi.org/ /501100002855ministry of science and technology of the people's republic of china2011aa02a1092012bah24f05jia hongbing http://dx.doi.org/ /501100001809national natural science foundation of china8120132281301402zhang wen chen chen issue-copyright-statement© the author 2016 