Whole genome analysis reveals new insights into the molecular characteristics of Clostridioides difficile NAP1/BI/027/ST1 clinical isolates in the People’s Republic of China

Abstract

Background

The epidemic new strain NAP1/BI/027/ST-1 of Clostridioides difficile (C. difficile) causes more severe coliti and a higher mortality rate than historical strains. However, C. difficile NAP1/BI/027/ST-1 (C. difficile RT027) infections have been rarely reported in Asia, particularly in China.

Purpose

The objective of this study was to strengthen the understanding of the molecular characterizations of C. difficile RT027 in China.

Patients and methods

Two C. difficile NAP1/BI/027/ST-1 were detected from two patients, and no additional isolates were found. Whole genome sequencing (WGS) was used to characterize two C. difficile RT027 isolates and control strain CD6 (from Hong Kong), and comparative genomic analysis was performed to compare genomic differences between seven isolates from Mainland China, CD6, and 10 isolates from North America and Europe.

Results

The comparative genomic analysis revealed that isolates obtained from Mainlan China were outside of the two epidemic lineages, FQR1 and FQR2, and might have decreased virulence and transmissibility for outbreak. Furthermore, unique SNP mutations were detected in isolates obtained from Mainland China, which may affect the biological function of C. difficile.

Conclusion

We speculate that C. difficile RT027 isolates in Mainland China may have different features, compared to those in North America and Europe.

Keywords:

• People's Republic of China
• multilocus sequence typing
• Clostridioides difficile infection
• whole genome sequencing

Introduction

Clostridioides difficile is widely recognized as an important diarrheal pathogen in North America and Europe. The mechanisms by which C. difficile causes severe disease in one individual, but silently colonizes another, are multifactorial but primarily associated with the use of antibiotics and other factors involving immune status, age, microbiota composition, and C. difficile strain types.¹ A hypervirulent strain of C. difficile, designated as NAP1/BI/027 or ST-1, caused more severe colitis and higher mortality than other types in North America and Europe since 2003.² Phylogeographic analysis of the whole genome sequencing data identified two major genetic lineages (FQR1 and FQR2) in C. difficile RT027, with independently acquired, identical mutations that convey high-level fluoroquinolone resistance (FQR), and traced their distinct patterns of global spread. The FQR1 lineage emerged in North America in 2001 and spread to South Korea and Switzerland, whereas the FQR2 lineage originated in 2003 in North America, but subsequently spread to the United Kingdom, continental Europe and Australia.³

The phylogenetic analysis has demonstrated that C. difficile is a genetically diverse lineage.⁴ Many different typing or fingerprinting methods have been applied to type and characterize C. difficile strains during an outbreak or infection, including PCR ribotyping and pulsed field gel electrophoresis (PFGE), among which PCR ribotyping has evolved as the standard typing method in Europe.⁵ Multilocus sequence typing (MLST) has gradually become popular in the molecular typing of multiple loci, which characterizes allelic variations of housekeeping genes, and assigns each particular allelic profile to a sequence type (ST).⁶

The use of whole-genome sequencing (WGS) has provided evidence for a higher degree of C. difficile strain diversity than previously acknowledged.⁷ WGS can be used to compare single nucleotide variants (SNVs) between isolates across the nonrepetitive core genome, which accounts for 80% of the 4.3 million-base pair C. difficile strain 630 reference genome.^8–10 WGS also offers additional benefits, including the reconstruction of long-term evolutionary histories and the in silico determination of virulence factors and antimicrobial resistance.¹¹

C. difficile RT027 first appeared in Mainland China in 2014, and subsequently caused sporadic cases, including a small healthcare-associated outbreak.^12–14 However, it has not replaced the other epidemic clones, such as ST54 (RT012), as the predominant clone across regions or the country.^15–17 This was quite different from North America and Europe, where hypervirulent C. difficile RT027 displaced endemic strains after it occurred. For example, RT027 has disseminated with a clear shift throughout Europe from UK to the Eastern Europe.^18,19 On the other hand, some countries have successfully controlled its spread and decreased its prevalence after taking measures.²⁰ The molecular epidemiology of C. difficile RT027 isolates from Mainland China may be different from other regions of the world. In the present study, we identified two clinical isolates of C. difficile, which were confirmed as C. difficile RT 027 by PCR ribotyping, NAP1 by PFGE and ST-1 by MLST and characterized by WGS. Then, comparative genomic analysis was performed to compare genomic differences between seven isolates from Mainland China and ten isolates from North America and Europe. This study will further our understanding of the molecular characterizations of C. difficile NAP1/BI/027/ST-1 in China.

Methods

Ethics

The present study is a retrospective study and ethical approval was given by the Medical Ethics Committee of The First Affiliated Hospital, School of Medicine, Zhejiang University (Reference Number: 2018–1020). None of the test results were used to alter individual care.

Study setting and general study design

This epidemiological study was conducted at the First Affiliated Hospital, School of Medicine, Zhejiang University, which is a 2,500-bed tertiary teaching hospital in Hangzhou, Zhejiang, China. Since September 2009, diagnosis of C. difficile infection (CDI) was carried out according to the clinician requirements using C. difficile culture and toxin gene tests.

Diarrhea was defined as three or more loose stools within 24 hours. Inpatients with diarrhea, whose stool samples were positive for both C. difficile culture and toxin gene tests and without evidence of another cause of diarrhea, were diagnosed as CDI. HA-CDI (Hospital-associated C. difficile infection) was defined as a patient with CDI symptoms onset more than 48 hours after admission to a healthcare-facility or with onset of symptoms in the community within 12 weeks following discharge from a health-care facility.²¹ CDI outbreak was defined as >2 isolates of the same genotype detected less than 7 days apart in one hospital either with onset of symptoms on the same ward, or accompanied by an increased CDI monthly incidence within the hospital.

Isolation and detection of C. difficile toxins

Stool samples were cultured on cycloserine-cefoxitin-fructose agar (CCFA) in an atmosphere composed of 80% N₂, 10% H₂ and 10% CO₂ at 37°C for 48 hours. The colonies were identified through the typical morphology and odor of C. difficile, and were further confirmed by Matrix-Assisted Laser Desorption-Ionization mass spectrometry (MALDI-MS, FlexControl 3.3-Microflex). Bacterial cultures from both stool specimens were negative for common diarrhea pathogens, including Salmonellaspp., Shigellaspp., Vibriospp., Staphylococcus aureus and Escherichia coli.

Genomic DNA was isolated from the two isolates using a QIAamp DNA Mini Kit (Qiagen NV, Hilden, Germany). The presence of tcdA and tcdB genes in two isolates was detected by PCR, as previously described.²² The presence of binary toxin genes cdtA, and cdtB were detected, as described by Stubbs et al.²³

A C. difficile NAP1/BI/027 strain, designated as CD6, which was kindly provided by Dr K. Y. Yuen (Department of Microbiology, Queen Mary Hospital, Hong Kong) was used as the positive control.²⁴

Capillary gel electrophoresis (CGE)-based PCR ribotyping and pulsed-field gel electrophoresis typing

The PCR-ribotyping of two isolates targeted the intergenic spacer region of C. difficile between the 16S and 23S rRNA genes.²⁵ The PCR products were analyzed by CGE. Strain CD6 was used as an internal control for PCR ribotyping. PFGE was performed, as previously described.²⁶

Multilocus sequence typing

MLST was performed by amplifying and sequencing seven housekeeping genes (adk, atpA, dxr, glyA, recA, sodA and tpi), according to previously described protocols.²⁷ MLST sequences were uploaded to the PubMLST database (http://pubmlst.org/cdifficile/) to determine the individual ST from the allelic combination.

Genome sequencing, assembly and annotation

The sequencing-quality genomic DNA of the three strains was extracted using a DNeasy Blood and Tissue Kit (Qiagen, Hilden, Germany), according to manufacturer’s instructions, and quantitated by Qubit 2.0. Whole-genome sequencing was performed on the Illumina HiSeq 2,500 System (Illumina, San Diego, CA, USA). Before performing the assembly, the sequencing reads were filtered to remove reads (<Q20) with low-quality base calls or reads that were similar to Illumina adapters. Subsequently, the raw reads of C. difficile RT027 strains were trimmed and mapped to the genome sequence of the C. difficile RT027 reference strain R20291 using BWA-SW to detect the genomic single nucleotide polymorphisms (SNPs) using a Genome Analysis Toolkit (GATK).^28,29 Finally, the genomes were submitted to the Rapid Annotation using Subsystem Technology (RAST) servers (http://rast.nmpdr.org/) for annotations.³⁰

Phylogenetic and comparative genomic analysis

In order to perform the comparative genomic analysis with previously published genome data, sequencing reads from representative C. difficile RT027 strains were downloaded from the National Center for Biotechnology Information (NCBI) genome database. As shown in Table 1, five strains reported by Jia et al¹⁴ were the first to cause a ward transmission in a hospital in Mainland China, whereas ten strains from the report of He et al³ were the node at the base of the star-like topology, and were associated with severe outbreaks in North America and Europe. SNPs were used to reconstruct the phylogenetic tree and maximum-likelihood tree by RAxML. Meanwhile, the heat-map reconstruction and analysis were performed using Roary Matrix for further comparative analysis.

Table 1 Information about the isolates in this study

Strain	Genome accession no.	Country or region	Year of strain isolated or reported	Reference
14,042,403		Mainland China	2014	This study
14,061,807		Mainland China	2014	This study
CD6		Hong Kong	2008	This study
P1	SRR1735363	Mainland China	2012	Jia et al^Citation14
P2	SRR1735364	Mainland China	2012	Jia et al^Citation14
P4	SRR1735366	Mainland China	2012	Jia et al^Citation14
P13	SRR1735375	Mainland China	2013	Jia et al^Citation14
P17	SRR1735381	Mainland China	2014	Jia et al^Citation14
BI-10	ERS017143	Pittsburgh, PA, USA	2001	He et al^Citation3
Can009	ERS017275	Montreal, QC, Canada	2003	He et al^Citation3
Can010	ERS017153	Montreal, QC, Canada	2003	He et al^Citation3
Can012	ERS017155	Montreal, QC, Canada	2003	He et al^Citation3
Lei013	ERS032938	Netherlands	2006	He et al^Citation3
Lei014	ERS032939	Netherlands	2006	He et al^Citation3
Lei015	ERS032940	Netherlands	2006	He et al^Citation3
Lei016	ERS032941	Netherlands	2006	He et al^Citation3
Lei018	ERS032943	Netherlands	2006	He et al^Citation3
Lei019	ERS032944	Netherlands	2006	He et al^Citation3

Antimicrobial susceptibility

Antimicrobial susceptibility testing was performed on Brucella agar plates containing 1 mg/L of vitamin K, 5 mg/L of hemin and 5% sheep red blood cells with the eight antimicrobial agents: metronidazole, vancomycin, clindamycin, erythromycin, linezolid, moxifloxacin, levofloxacin, and rifampicin. The minimum inhibitory concentration (MIC) breakpoints used were 8 mg/L for erythromycin, clindamycin, and the fluoroquinolones, and 32 mg/L for metronidazole, in accordance with Clinical Laboratory Stsandards Institute (CLSI) interpretative categories approved for anaerobic bacteria.³¹ For vancomycin, linezolid and rifampicin, in which the breakpoints were not available in CLSI documents, >2 mg/L, >4 mg/L and >32 mg/L were used as the breakpoints, respectively, according to the European Committee on Antimicrobial Susceptibility Testing (EUCAST) (http://www.eucast.org/clinical_breakpoints/). C. difficile ATCC 70,057 was used as a control.

Results

Epidemiological analysis and characterization of C. difficile isolates

In total, 421 (8.1%) non-duplicate toxigenic C. difficile isolates were identified from 5,171 patients suffering from diarrhea (parts data has been published) during the period of study (by the end of August 31, 2017).³² We have identified 33 STs among these isolates. The most prevalent STs were ST-54, ST-35, ST-37, and ST-3, followed by ST-2, ST-81, and ST-8. However, STs with binary toxin genes, such as ST-11, ST-5, and ST-201, were not common.³²

All positive strains of C. difficile isolated from stool samples were retained in our laboratory, from which two clinical strains were tested positive for both toxin A and toxin B, and binary toxin genes by PCR assay, which were further confirmed as ST-1 by MLST. Compared to CD6, these two isolates were assigned as C. difficile RT027 by PCR ribotyping and NAP1 type by PFGE (Figure 1A and B).

Figure 1 (A) PCR ribotyping of Clostridioides difficile (C. difficile) PCR ribotype 027 reference strains (A-01) and two clinical isolates (A-02 and B-03) identified in the present study. (B) Pulsed-field gel electrophoresis of C. difficile PCR ribotype 027 reference strains (A-01) and two clinical isolates (A-02 and B-03). Marker: Xba 1 was used for the reference marker, Salmonella H9812. Line A-01=CD06; Line A-02=14,042,403; Line B-03=14,061,807.

C. difficile strain 14,042,403 was isolated from the stool sample from a 71-year-old man with diarrhea during the hospitalization. The patient was treated with antibiotics for intra-abdominal infection before he was diagnosed as CDI. C. difficile strain 14,061,807 was isolated from an 83-year-old hospitalized man in the same hospital. This patient was also treated with different antibiotics for pulmonary infection and continued to suffer from persistent fever with persistent diarrhea.

General genomic features

After filtering to remove low quality and adapter contamination reads, 1,656,804, 5,048,333 and 5,248,180 paired reads were produced for CD6, 14,042,403 and 14,061,807 strains. Following the assembly, the average genome coverage was approximately 290× and the genome assembly produced 80,110,116 contigs, respectively, with an average GC content between 29.2% and 32.8%. The features of these three strains, including the assembly stats and gene content, are shown in Table S1.

Phylogenetic analyses

In order to further understand the phylogenetic relationship between the two isolates in the present study, paired-end DNA sequencing reads were mapped to the reference genomes of R20291. The two strains were highly similar with a difference of only four SNPs between them. Among these four different SNPs as shown in Figure 2A, two were in the non-coding regions, and two were in the coding regions (conserved hypothetical protein and ABC transporter ATP-binding protein). The phylogenetic tree revealed that these two strains were closely clustered into the same predominant evolutionary branch, which is similar to the other isolates obtained from Mainland China. However, the maximum-likelihood tree revealed different group patterns, in which two isolates were separated into the same clade, but were far away from the other strains (Figure 2B).

Figure 2 Patterns of molecular variation in 18 strains of C. difficile RT027, including the three isolates in this study, the other five strains from the People's Republic of China and the ten from North America and Europe. (A) Annotation by RAST revealed several SNPs, which were significantly different between the strains obtained from Mainland China and North America and Europe. The symbol arrow indicates the four different SNPs between two isolates in this study. (B) Maximum-likelihood phylogenetic analysis of based on core genome SNPs which were grouped into three different clades. The previously described strains of fluoroquinolone-resistant C. difficile RT027 were clustered into two clades (FQR1 and FQR 2), whereas the isolates from Mainland China were classified into a third clade.

In order to investigate the difference and relationship between C. difficile RT027 strains isolated from Mainland China and those obtained from North America and Europe, all WGS data of available C. difficile RT027 from Mainland China (by the end of October 2017) and ten isolates from North America and Europe were downloaded to align with the genome of the C. difficile RT027 strain R20291. A total of 628 SNPs were identified and grouped into three different clades (Figure 2B). The previously described strains of fluoroquinolone-resistant C. difficile RT027 were clustered into two clades (FQR1 and FQR 2), whereas the isolates from Mainland China were classified into a third clade.³ However, CD6, the strain isolated from Hong Kong, was grouped into the FQR2 lineage. Annotation by RAST revealed several SNPs, which were significantly different between the strains obtained from Mainland China and North America and Europe. These were located at significant protein-coding sites, including the conserved hypothetical protein, putative membrane protein, putative exported protein and two-component response regulator (Table S2). In addition, a total of 3,944 genes were identified (Figure 3), including 3,397 common core genes without soft core genes, 295 different shell genes, and 252 different cloud genes. The annotation for identification of the predicted function of proteins coded by different genes revealed that 12 kinds of proteins were only encoded by unique genes in the seven strains obtained from Mainland China, including excisionase and transposase from transposon Tn916 (Table S3). Bifunctional AAC/APH, which belongs to the aminoglycoside phosphotransferase family, and ribosomal-protein-alanine N-acetyltransferase, which acetylates the N-terminal alanine residues of specific ribosomal proteins, were the unique proteins in Chinese C. difficile RT027 isolates.

Figure 3 Heatmap reconstruction of 18 strains resulting from SNPs. A total of 3,944 genes were identified, including 3,397 common core genes without soft core genes, 295 different shell genes, and 252 different cloud genes. Two isolates in this study are labeled with point brackets.

Antimicrobial susceptibility

The antibiotic susceptibility patterns of C. difficile isolates are presented in Table 2. Strain 14,042,403 and 14,061,807 were resistant to most of the tested antibiotics, including fluoroquinolones with high MICs, but susceptible to metronidazole, vancomycin and linezolid.

Table 2 Antimicrobial susceptibility of the four C. difficile isolates

No.	Value (μg mL^−1)
	MZ	LZ	VA	CM	RI	LE	EM	MX
ATCC70057	0.023	0.5	0.38	1.5	<0.002	3	0.5	0.38
CD6	0.125	0.5	0.5	2	<0.002	>32	>256	>32
14,042,403	0.25	0.5	1.2	>256	>32	>32	>256	>32
14,061,807	0.75	0.75	1.5	>256	>32	>32	>256	>32

Discussion

In 2014, two strains of C. difficile NAP1/BI/027/ST-1 were successfully isolated from a tertiary hospital in East China. To the best of our knowledge, this is the first report of C. difficile NAP1/BI/027/ST-1 infection in this region.

C. difficile RT027 infection is not common in Mainland China, even in almost all regions of Asia. In the present study, we report two cases infected by C. difficile RT027, of which one patient recovered and another patient died. Although the two C. difficile RT027 strains were highly similar, with only differences in four SNPs, the genotypic characteristic demonstrates that the two isolates were not derived from the same clone, according to the previously estimated evolutionary rate of 0.74 SNP per year for C. difficile, indicating no direct transmission of C. difficile between patients hospitalized at the same time, although both cases had a history of sharing the same ward (Gastroenterology Department) for 2 days.⁷ It has been reported that exposure to fluoroquinolones was an independent risk factor for the spread of C. difficile RT027, while the consumption of fluoroquinolones in the People's Republic of China was similar to that in Europe and North America.^33,34 The two isolates were multi-drug resistant strains and also resistant to fluoroquinolones. All these evidences demonstrate the occurrence of C. difficile RT027 strains in this region, but there were no more isolates of C. difficile RT027 in the subsequent surveillance, which has been done since 2009.³² The review of relevant reports on C. difficile RT027 revealed no evidence to support the predominance of hypervirulent C. difficile RT027 strains in Mainland China, which was quite different from Europe or North America, where the prevalence of C. difficile RT027 emerged after it occurred.³⁵ This may be due to the fact that the lack of adequate laboratory diagnostic capacity and high-quality and multiple-center surveillance system for CDI, and the low submission rate of samples in the People's Republic of China. However, in this study, we tried to explain the difference by comparative genomic analysis.

Two distinct epidemic lineages of C. difficile RT027, FQR1 or FQR2, emerged in North America within a relatively short period after acquiring the same fluoroquinolone resistance-conferring mutations.³ The present results confirm that for all C. difficile isolates detected from Mainland China, all of them were outside of both FQR lineages, which might have decreased virulence and transmissibility for outbreak. In the study reported by He et al³ that the isolates outside of the both lineages represent the pre-epidemic C. difficile RT027 genetic background, from which FQR epidemic lineages emerge may support our hypothesis.³ Although the two isolates in this study were both resistant to levofloxacin and moxifloxacin while they had mutations in DNA gyrase subunit A gene gyrA (Asp87→Tyr), this mutation was different from an identical mutation (Thr82→Ile) in another C. difficile RT027.³⁶ A future study is needed to know if these mutations are hotspots within the so-called fluoroquinolone resistance determining region.

The present study revealed that several unique SNPs of C. difficile RT027 strains isolated in Mainland China were located at significant protein-coding sites that might affect the biological function of C. difficile. Hundreds of two-component response regulatory systems that allow organisms to sense and respond to environmental changes have now been found in bacterial genomes.³⁷ The sporulation program is one of two-component systems, and capsular polysaccharides play important roles in the mechanism for C. difficile to efficiently persist in the environment and transmit the disease.³⁸ SNP alterations in C. difficile strains from Mainland China change codons from ACT (threonine) to GCT (alanine) in amino acid sequences of the two-component response regulator involved in the sporulation program, and cause proline (CCG) to be replaced by leucine (CTG) in the capsular polysaccharide biosynthesis protein, respectively. However, although the mutation of particular protein-coding sites can affect protein function in two-component systems and capsular polysaccharides, it remains to be determined whether SNP mutations that result in amino acid substitutions are associated with the epidemiological features of C. difficile RT027 in the People's Republic of China.

Several special proteins have been detected in C. difficile isolates from Mainland China, including excisionase and transposase from transposon Tn916. Tn916 contains genes that encode a tyrosine recombinase and an excisionase, which are responsible for the integration and excision of an element.³⁹ Acting as one of the most widespread conjugative transposons, Tn916 is responsible for the dissemination of resistance genes, can be transferred from a donor cell to a recipient, and can regulate the resistance gene by transcriptional attenuation.^40,41 Further studies are needed to understand the molecular mechanism of SNP mutations in protein-coding sites that affect the features of C. difficile RT027 in the People's Republic of China.

The present study has several limitations. WGS results revealed significant genomic differences in C. difficile RT027 between Mainland China and North America and Europe, and the effects of SNP mutations on pathogenesis and transmission remain largely unclear. Secondly, only two C. difficile RT027 strains were identified and it is not enough data to draw conclusions about differences in bacterial strains on the ability of C. dfficile to cause epidemic disease in the People's Republic of China. However, we collected the whole genome data of RT027 reported in China, and combined with our data to analyze the difference between RT027 in the People's Republic of China and those in North America and Europe. Thirdly, we didn’t confirm the source of the C. difficile RT027 due to no environmental samples taking and no more C. difficile RT027 being identified.

In conclusion, two isolates of C. difficile NAP1/BI/027/ST-1 were identified from two patients, and no outbreak occurred. We speculate that the isolates of C. difficile RT027 isolated from Mainland China may have lower levels of transmissibility when compared to those in North America and Europe. The comparative genomic analysis revealed that C. difficile RT027 isolates from Mainland China were outside of the two epidemic lineages, FQR1 and FQR2, and might also have different epidemiological features. Unique SNP mutations were detected in isolates from Mainland China, leading to the probable alterations of the biological functions of C. difficile RT027. The results of the present study strengthen the understanding of the molecular characterizations of C. difficile RT027 in the People's Republic of China.

Availability of data and materials

The Whole Genome Sequencing of C. difficile isolates 14061807, CD6 and 14042403 have been deposited in GenBank under the following accession number: PNER00000000, PNES00000000, PNET0000000, respectively.

Author contributions

All authors made substantial contributions to conception and design, acquisition of data, or analysis and interpretation of data; took part in drafting the article or revising it critically for important intellectual content; gave final approval of the version to be published; and agree to be accountable for all aspects of the work.

Acknowledgments

The authors wish to thank Dr KY Yuen and WC Yam of the Department of Microbiology Queen Mary Hospital Faculty of Medicine, University of Hong Kong, who provided the PCR ribotype027 strain of C. difficile, CD6. This study was partially funded by grants from the National Key Research and Development Program of China (2017YFC1200203), the Mega-projects of Science research of China (2018ZX10712001-005) and the Zhejiang Provincial Natural Science Foundation of China under Grant No. LY18H040001.

Supplementary materials

Table S1 Genomic features of the three strains

	Paired reads	Total length (Mbp)	Read length (bp)	Average GC (%)	Contigs	Bases	Genes	CDS	tRNA	tmRNA
14,042,403	5,048,333	757.24,995*2	150	32.8	110	4,115,741	3719	3650	68	1
14,061,807	5,248,180	787.227*2	150	32.3675	116	4,118,129	3710	3653	56	1
CD6	1,656,804	248.5206*2	150	29.2257	80	4,129,319	3738	3689	48	1

Table S2 SNPs located within the protein-coding genes

Proteins	Code	Number	Position
2-isopropylmalate synthase	CBE02932	1
ATP-dependent RNA helicase	CBE05139	1
Capsular polysaccharide biosynthesis protein	CBE06153	1	3,128,507
Cell surface protein (S-layer precursor protein)	CBE06198	1	3,163,982
Conserved hypothetical protein	CBE03238	1	1,232,712
Conserved hypothetical protein	CBE06254	1	3,207,237
Cyclomaltodextrinase	CBE05153	1
DNA mismatch repair protein	CBE04805	1
DNA-directed RNA polymerase beta chain	CBE01567	1
Electron transfer flavoprotein beta-subunt	CBE02743	1
Hypothetical protein	CBE05871	1	2,931,889
Penicillin-binding protein	CBE02705	1	879,963
Peptide deformylase 2	CBE05821	1
Phosphoenolpyruvate-protein phosphotransferase	CBE06128	1
Putative 5-nitroimidazole reductase	CBE03766	2	1,547,479;1,547,553
Putative atp phosphoribosyltransferase regulatory subunit	CBE03922	1
Putative drug/sodium antiporters	CBE03848	2	1,600,434
Putative exported protein	CBE05745	1	2,848,439
Putative membrane protein	CBE02014	1	384,000
Putative penicillin-binding protein	CBE03184	2	1,202,866;1,203,554
Putative ribosomal protein	CBE01618	1
Putative signaling protein	CBE05464	1
Quinolinate synthetase A	CBE05440	1	2,649,551
Two-component response regulator	CBE04143	1	1,794,733

Table S3 Annotation proteins for the unique genes of C. difficile in mainland China

Annotation protein	Strains
	14,042,403	14,061,807	CD6	FQR1_BI_10	FQR2_Can009	FQR2_Can010	FQR2_Can012	Lei_013	Lei_014	Lei_015	Lei_016	Lei_018	Lei_019	P1	P13	P17	P2	P4
Helix-turn-helix	14042403_00781	14061807_01413												P1_01576	P13_02880	P17_02821	P2_02728	P4_01362
Hypothetical protein	14042403_00782	14061807_01412												P1_01575	P13_02879	P17_02820	P2_02727	P4_01361
Ribosomal-protein-alanine N-acetyltransferase	14042403_00783	14061807_01411												P1_01574	P13_02878	P17_02819	P2_02726	P4_01360
Bifunctional AAC/APH	14042403_00784	14061807_01410												P1_01573	P13_02877	P17_02818	P2_02725	P4_01359
Transposase, Mutator family	14042403_00785	14061807_01409												P1_01572	P13_02876	P17_02817	P2_02724	P4_01358
Internalin-J precursor	14042403_00786	14061807_01408												P1_01571	P13_02875	P17_02816	P2_02723	P4_01357
HTH-type transcriptional regulator SinR	14042403_00787	14061807_01407												P1_01570	P13_02874	P17_02815	P2_02722	P4_01356
Hypothetical protein	14042403_00788	14061807_01406												P1_01569	P13_02873	P17_02814	P2_02721	P4_01355
Excisionase from transposon Tn916	14042403_00789	14061807_01405												P1_01568	P13_02872	P17_02813	P2_02720	P4_01354
Transposase from transposon Tn916	14042403_00790	14061807_01404												P1_01567	P13_02871	P17_02812	P2_02719	P4_01353
Integrase core domain protein	14042403_01013	14061807_01055												P1_01579	P13_02883	P17_02824	P2_02731	P4_01365
Transposase	14042403_01014	14061807_01056												P1_01580	P13_02884	P17_02825	P2_02732	P4_01366

Disclosure

The authors report no conflicts of interest in this work.