The mosaic nature of intergenic 16S-23S rRNA spacer regions suggests rRNA operon copy number variation in Clostridium difficile strains

Clostridioides difficile infection: history, epidemiology, risk factors, prevention, clinical manifestations, treatment, and future options

SUMMARYClostridioides difficile infection (CDI) is one of the major issues in nosocomial infections. This bacterium is constantly evolving and poses complex challenges for clinicians, often encountered in real-life scenarios. In the face of CDI, we are increasingly equipped with new therapeutic strategies, such as monoclonal antibodies and live biotherapeutic products, which need to be thoroughly understood to fully harness their benefits. Moreover, interesting options are currently under study for the future, including bacteriophages, vaccines, and antibiotic inhibitors. Surveillance and prevention strategies continue to play a pivotal role in limiting the spread of the infection. In this review, we aim to provide the reader with a comprehensive overview of epidemiological aspects, predisposing factors, clinical manifestations, diagnostic tools, and current and future prophylactic and therapeutic options for C. difficile infection.

Accounting for 16S rRNA copy number prediction uncertainty and its implications in bacterial diversity analyses

16S rRNA gene copy number (16S GCN) varies among bacterial species and this variation introduces potential biases to microbial diversity analyses using 16S rRNA read counts. To correct the biases, methods have been developed to predict 16S GCN. A recent study suggests that the prediction uncertainty can be so great that copy number correction is not justified in practice. Here we develop RasperGade16S, a novel method and software to better model and capture the inherent uncertainty in 16S GCN prediction. RasperGade16S implements a maximum likelihood framework of pulsed evolution model and explicitly accounts for intraspecific GCN variation and heterogeneous GCN evolution rates among species. Using cross-validation, we show that our method provides robust confidence estimates for the GCN predictions and outperforms other methods in both precision and recall. We have predicted GCN for 592605 OTUs in the SILVA database and tested 113842 bacterial communities that represent an exhaustive and diverse list of engineered and natural environments. We found that the prediction uncertainty is small enough for 99% of the communities that 16S GCN correction should improve their compositional and functional profiles estimated using 16S rRNA reads. On the other hand, we found that GCN variation has limited impacts on beta-diversity analyses such as PCoA, NMDS, PERMANOVA and random-forest test.

Systems biology approach to functionally assess the Clostridioides difficile pangenome reveals genetic diversity with discriminatory power

SignificanceClostridioides difficile infections are the most common source of hospital-acquired infections and are responsible for an extensive burden on the health care system. Strains of the C. difficile species comprise diverse lineages and demonstrate genome variability, with advantageous trait acquisition driving the emergence of endemic lineages. Here, we present a systems biology analysis of C. difficile that evaluates strain-specific genotypes and phenotypes to investigate the overall diversity of the species. We develop a strain typing method based on similarity of accessory genomes to identify and contextualize genetic loci capable of discriminating between strain groups.

16S rDNA and ITS Sequence Diversity of Burkholderia mallei Isolated from Glanders-Affected Horses and Mules in India (2013-2019)

Glanders is a highly contagious and fatal infection of equids caused by the bacteria known as Burkholderia mallei. It is one of the notifiable equine diseases and is still present in Asia, South America and Africa. In India, glanders re-emerged in 2006, and thereafter, increasing numbers of cases were reported in different regions of the country. Between 2013 and 2019, 39 B. mallei were isolated from glanders-affected horses (n = 30) and mules (n = 9) from seven states of India such as Uttar Pradesh, Haryana, Delhi, Himachal Pradesh, Gujarat, Maharashtra and Tamil Nadu. In this study, the phylogenetic relationships of these isolates were assessed by sequence analysis of 16S rDNA gene and ITS region. Purified PCR-amplified products of 16S rDNA gene and ITS region were sequenced, aligned and phylogenetic trees were constructed using MEGA 11 software. Additionally, B. mallei 16S rDNA (n = 36) and ITS (n = 18) sequences available in the GenBank were also included for analysis to determine the diversity of older B. mallei isolates with recent Indian isolates. Both the phylogeny showed that the majority of the recent isolates from India are closely related to each other, but are genetically diverse from older isolates that originated from India. Nucleotide substitutions were also observed in a single and double position in 12 recent and two old Indian isolates. The study also indicates that similar B. mallei strains were responsible for glanders outbreaks in different states (Uttar Pradesh- Himachal Pradesh and Uttar Pradesh- Haryana) and this is due to the migration of infected animals from one state to another state. This study implies that 16S rDNA and ITS region may be used for molecular characterization of B. mallei associated with glanders in resource-limited settings.

Transition From PCR-Ribotyping to Whole Genome Sequencing Based Typing of Clostridioides difficile

Clostridioides difficile causes nosocomial outbreaks which can lead to severe and even life-threatening colitis. Rapid molecular diagnostic tests allow the identification of toxin-producing, potentially hypervirulent strains, which is critical for patient management and infection control. PCR-ribotyping has been used for decades as the reference standard to investigate transmission in suspected outbreaks. However, the introduction of whole genome sequencing (WGS) for molecular epidemiology provides a realistic alternative to PCR-ribotyping. In this transition phase it is crucial to understand the strengths and weaknesses of the two technologies, and to assess their correlation. We aimed to investigate ribotype prediction from WGS data, and options for analysis at different levels of analytical granularity. Ribotypes cannot be directly determined from short read Illumina sequence data as the rRNA operons including the ribotype-defining ISR fragments collapse in genome assemblies, and comparison with traditional PCR-ribotyping results becomes impossible. Ribotype extraction from long read Oxford nanopore data also requires optimization. We have compared WGS-based typing with PCR-ribotyping in nearly 300 clinical and environmental isolates from Switzerland, and in addition from the Enterobase database (n=1778). Our results show that while multi-locus sequence type (MLST) often correlates with a specific ribotype, the agreement is not complete, and for some ribotypes the resolution is insufficient. Using core genome MLST (cgMLST) analysis, there is an improved resolution and ribotypes can often be predicted within clusters, using cutoffs of 30-50 allele differences. The exceptions are ribotypes within known ribotype complexes such as RT078/RT106, where the genome differences in cgMLST do not reflect the ribotype segregation. We show that different ribotype clusters display different degrees of diversity, which could be important for the definition of ribotype cluster specific cutoffs. WGS-based analysis offers the ultimate resolution to the SNP level, enabling exploration of patient-to-patient transmission. PCR-ribotyping does not sufficiently discriminate to prove nosocomial transmission with certainty. We discuss the associated challenges and opportunities in a switch to WGS from conventional ribotyping for C. difficile.

Rapid Classification of Clostridioides difficile Strains Using MALDI-TOF MS Peak-Based Assay in Comparison with PCR-Ribotyping

Typing methods are needed for epidemiological tracking of new emerging and hypervirulent strains because of the growing incidence, severity and mortality of Clostridioides difficile infections (CDI). The aim of this study was the evaluation of a typing Matrix-Assisted Desorption/Ionization-Time of Flight Mass Spectrometry (MALDI-TOF MS (T-MALDI)) method for the rapid classification of the circulating C. difficile strains in comparison with polymerase chain reaction (PCR)-ribotyping results. Among 95 C. difficile strains, 10 ribotypes (PR1-PR10) were identified by PCR-ribotyping. In particular, 93.7% of the isolates (89/95) were grouped in five ribotypes (PR1-PR5). For T-MALDI, two classifying algorithm models (CAM) were tested: the first CAM involved all 10 ribotypes whereas the second one only the PR1-PR5 ribotypes. Better performance was obtained using the second CAM: recognition capability of 100%, cross-validation of 96.6% and agreement of 98.4% (60 correctly typed strains, limited to PR1-PR5 classification, out of 61 examined strains) with PCR-ribotyping results. T-MALDI seems to represent an alternative to PCR-ribotyping in terms of reproducibility, set up time and costs, as well as a useful tool in epidemiological investigation for the detection of C. difficile clusters (either among CAM included ribotypes or out-of-CAM ribotypes) involved in outbreaks.

Comparison of Clostridioides difficile strains from animals and humans: First results after introduction of C. difficile molecular typing and characterization at the Istituto Zooprofilattico Sperimentale of Piemonte, Liguria e Valle d'Aosta, Italy

PCR ribotypes (RTs027 and 078) are known causes of Clostridioides difficile infection (CDI) in humans. Molecular typing and characterization of 39 C. difficile strains isolated from samples from humas and animals in 2016-2018 indicated an overlap of RTs between community-acquired patients (CA-CDI) and domestic animals from the same geographical area; 14 RTs were identified: 12 RTs were positive for toxins A/B; RT078, RT080 and RT126 were also positive for binary toxin (CDT). Most of the RTs from the animals (RTs020, 078, 106, 126) were also detected in the samples from humans. Strains grouped into three clusters: cluster I included prevalently human strains, mainly RT 018; clusters II and III included strains from humans and animals, mainly RT078 and RT020. The CA-CDI strains suggested animals as a reservoir of C. difficile isolated together with other microorganisms from animals, highlighting the association of enteric pathogens as a cause of infection and death.

Does Intraspecific Variation in rDNA Copy Number Affect Analysis of Microbial Communities?

Amplicon sequencing of partial regions of the ribosomal RNA loci (rDNA) is widely used to profile microbial communities. However, the rDNA is dynamic and can exhibit substantial interspecific and intraspecific variation in copy number in prokaryotes and, especially, in microbial eukaryotes. As change in rDNA copy number is a common response to environmental change, rDNA copy number is not necessarily a property of a species. Variation in rDNA copy number, especially the capacity for large intraspecific changes driven by external cues, complicates analyses of rDNA amplicon sequence data. We highlight the need to (i) interpret amplicon sequence data in light of possible interspecific and intraspecific variation, and (ii) examine the potential plasticity in rDNA copy number as an important ecological factor to better understand how microbial communities are structured in heterogeneous environments.

Molecular epidemiology of predominant and emerging Clostridioides difficile ribotypes

There has been an increase in the incidence and severity of Clostridioides difficile infection (CDI) worldwide, and strategies to control, monitor, and diminish the associated morbidity and mortality have been developed. Several typing methods have been used for typing of isolates and studying the epidemiology of CDI; serotyping was the first typing method, but then was replaced by pulsed-field gel electrophoresis (PFGE). PCR ribotyping is now the gold standard method; however, multi locus sequence typing (MLST) schemes have been developed. New sequencing technologies have allowed comparing whole bacterial genomes to address genetic relatedness with a high level of resolution and discriminatory power to distinguish between closely related strains. Here, we review the most frequent C. difficile ribotypes reported worldwide, with a focus on their epidemiology and genetic characteristics.

Species-Specific Identification of Streptococcus based on DNA Marker in 16S-23S rDNA Internal Transcribed Spacer

Streptococcus is closely correspondent to human. The accurate species-specific identification method of Streptococcus is important for the bacteria clinical diagnosis, molecular epidemiological analysis, and microecological study. In the last decades, DNA markers are widely utilized for identification of prokaryotic species. However, 16S rDNA, the most popular bacterial DNA marker, cannot properly distinguish closely related Streptococcus species. In present study, we employed 16S-23S rRNA gene internal transcribed spacer (ITS) sequence to explore the species-specific DNA marker. We predicted the secondary structure of Streptococcus ITS sequence transcribed products. Then we identified that the specific and consensus sequences in the primary structure can be found occupying an individual subunit in the secondary structure, which explained the foundation of the mosaic-like structure of ITS. We evaluated the specificity of ITS in Streptococcus, and found that the specificity can be detected by a further analysis of a BLAST result. Then, we developed an identification procedure based on the ITS sequence. We verified the procedure by 500 ITS sequence. The accuracy rate of this procedure was 100% for Streptococcus at genus level, and 99.3% at species level. It suggested that ITS can be utilized to accurately identify Streptococcus at the species level. This work suggests that further exploration of ITS could be applied in other bacterial genera for identification and classification, which may be a useful topic for future microbiology studies.

Species-specific Identification of Vibrio sp. based on 16S-23S rRNA gene internal transcribed spacer

AIMS

To explore a prokaryotic species-specific DNA marker, 16S-23S rRNA gene internal transcribed spacer (ITS) sequence for identification and classification of Vibrio.

METHODS AND RESULTS

Five hundred and seventy four ITS sequences from 60 Vibrio strains were collected, then the primary and secondary structures of ITS sequence were analysed. The ITS was divided into several subunits, and the species-specificity of these subunits were evaluated by blast. The variable subunit of ITS showed high species-specificity. A protocol to identify a Vibrio species based on ITS analysis was developed and verified. Both the specificity and sensitivity were 100%. The phylogeny analysis of Vibrio based on ITS showed that ITS devised a better classification than 16S rDNA. Finally, an identification method of Vibrio based on ITS sequencing in food samples was developed and evaluated. The results of ITS sequencing were (100%) consistent with the results identified by ISO standard.

CONCLUSIONS

Vibrio could be accurately identified at the species level by using the ITS sequences.

SIGNIFICANCE AND IMPACT OF THE STUDY

The present study suggests that the ITS can be considered as a significant DNA marker for identification and classification of Vibrio species, and it posed a new path to screen the Vibrio in food sample.

Development and Implementation of Whole Genome Sequencing-Based Typing Schemes for Clostridioides difficile

Clostridioides difficile is an important nosocomial pathogen increasingly observed in the community and in different non-human reservoirs. The epidemiology and transmissibility of C. difficile has been studied using a variety of typing methods, including more recently developed whole-genome sequence (WGS) analysis that is becoming used routinely for bacterial typing worldwide. Here we review the schemes for WGS-based typing methods available for C. difficile and their applications in the field of human C. difficile infection (CDI). The two main approaches to discover genomic variations are single nucleotide variant (SNV) analysis and methods based on gene-by-gene comparisons (frequently called core genome or whole genome MLST, cgMLST, or wgMLST). SNV analysis currently provides the ultimate resolution, however, typing nomenclature and standardized methodology are missing. On the other hand, gene-by-gene approaches allow portability and standardized nomenclature, and are therefore becoming increasingly popular in bacterial epidemiology and outbreak investigation. Two commercial software packages (BioNumerics and Ridom SeqSphere+) and an open source database (EnteroBase) for allele and sequence type determination for C. difficile are currently available. Proof-of-concept WGS studies have already enabled advances in the investigation of the population structure of C. difficile species, microevolution within the epidemic strains, intercontinental transmission over time and in tracking of transmission events. WGS of clinical C. difficile isolates demonstrated a considerable genetic diversity suggesting diverse reservoirs for CDI. WGS was also shown to aid in resolving relapses and reinfections in recurrent CDI and has potential for use as a tool for assessing hospital infection prevention and control performance.

Relationship of diversity and the secondary structure in 16S-23S rDNA internal transcribed spacer: a case in Vibrio parahaemolyticus

The 16S-23S rDNA internal transcribed spacer (ITS) sequence, located in the rrn operon, has been analyzed and evaluated for use in phylogenetic analysis and the detection target of bacteria. The ITS region displays a high level of diversity, being present in multiple copies and displaying variability in both length and sequence, and it carries more phylogenetic information than 16S rDNA. However, appropriately identifying ITS regions to use in analyses is challenging. To solve this problem, we analyzed the ITS regions in Vibrio parahaemolyticus and predicted the secondary structure of each analogous rrn transcript. The genomic DNA of V. parahaemolyticus contains approximately 8-14 rrns, making it more complex than the sequences of most other bacterial species. We analyzed 216 ITSs, of which 206 ITSs come from 18 complete genomes, and 10 ITSs were identified in the present study. The subunits of each ITS were distinguished by their predicted secondary structures. We propose a refined backbone model of the V. parahaemolyticus ITS that can be applied to the sequences of other bacteria. The backbone includes C, V, tDNA and linker blocks. These blocks, which may represent true functional units, may be used as potential targets for phylogenetic analysis or molecular detection.

Development and validation of an internationally-standardized, high-resolution capillary gel-based electrophoresis PCR-ribotyping protocol for Clostridium difficile

PCR-ribotyping has been adopted in many laboratories as the method of choice for C. difficile typing and surveillance. However, issues with the conventional agarose gel-based technique, including inter-laboratory variation and interpretation of banding patterns have impeded progress. The method has recently been adapted to incorporate high-resolution capillary gel-based electrophoresis (CE-ribotyping), so improving discrimination, accuracy and reproducibility. However, reports to date have all represented single-centre studies and inter-laboratory variability has not been formally measured or assessed. Here, we achieved in a multi-centre setting a high level of reproducibility, accuracy and portability associated with a consensus CE-ribotyping protocol. Local databases were built at four participating laboratories using a distributed set of 70 known PCR-ribotypes. A panel of 50 isolates and 60 electronic profiles (blinded and randomized) were distributed to each testing centre for PCR-ribotype identification based on local databases generated using the standard set of 70 PCR-ribotypes, and the performance of the consensus protocol assessed. A maximum standard deviation of only ±3.8bp was recorded in individual fragment sizes, and PCR-ribotypes from 98.2% of anonymised strains were successfully discriminated across four ribotyping centres spanning Europe and North America (98.8% after analysing discrepancies). Consensus CE-ribotyping increases comparability of typing data between centres and thereby facilitates the rapid and accurate transfer of standardized typing data to support future national and international C. difficile surveillance programs.

Complete genome sequence of the Clostridium difficile laboratory strain 630Δerm reveals differences from strain 630, including translocation of the mobile element CTn5

BACKGROUND

Clostridium difficile strain 630Δerm is a spontaneous erythromycin sensitive derivative of the reference strain 630 obtained by serial passaging in antibiotic-free media. It is widely used as a defined and tractable C. difficile strain. Though largely similar to the ancestral strain, it demonstrates phenotypic differences that might be the result of underlying genetic changes. Here, we performed a de novo assembly based on single-molecule real-time sequencing and an analysis of major methylation patterns.

RESULTS

In addition to single nucleotide polymorphisms and various indels, we found that the mobile element CTn5 is present in the gene encoding the methyltransferase rumA rather than adhesin CD1844 where it is located in the reference strain.

CONCLUSIONS

Together, the genetic features identified in this study may help to explain at least part of the phenotypic differences. The annotated genome sequence of this lab strain, including the first analysis of major methylation patterns, will be a valuable resource for genetic research on C. difficile.

Recombination drives evolution of the Clostridium difficile 16S-23S rRNA intergenic spacer region

PCR-ribotyping, a typing method based on size variation in 16S-23S rRNA intergenic spacer region (ISR), has been used widely for molecular epidemiological investigations of C. difficile infections. In the present study, we describe the sequence diversity of ISRs from 43 C. difficile strains, representing different PCR-ribotypes and suggest homologous recombination as a possible mechanism driving the evolution of 16S-23S rRNA ISRs. ISRs of 45 different lengths (ranging from 185 bp to 564 bp) were found among 458 ISRs. All ISRs could be described with one of the 22 different structural groups defined by the presence or absence of different sequence modules; tRNAAla genes and different combinations of spacers of different lengths (33 bp, 53 bp or 20 bp) and 9 bp direct repeats separating the spacers. The ISR structural group, in most cases, coincided with the sequence length. ISRs that were of the same lengths had also very similar nucleotide sequence, suggesting that ISRs were not suitable for discriminating between different strains based only on the ISR sequence. Despite large variations in the length, the alignment of ISR sequences, based on the primary sequence and secondary structure information, revealed many conserved regions which were mainly involved in maturation of pre-rRNA. Phylogenetic analysis of the ISR alignment yielded strong evidence for intra- and inter-homologous recombination which could be one of the mechanisms driving the evolution of C. difficile 16S-23S ISRs. The modular structure of the ISR, the high sequence similarities of ISRs of the same sizes and the presence of homologous recombination also suggest that different copies of C. difficile 16S-23S rRNA ISR are evolving in concert.

Loop-mediated isothermal amplification assays for detecting Yersinia pseudotuberculosis in milk powders

Yersinia pseudotuberculosis is a Gram-negative foodborne pathogen that causes several diseases, such as enteritis, septicemia, and reactive arthritis. Loop-mediated isothermal amplification (LAMP) assay targeting the 16S-23S rDNA internal transcribed spacer (ITS) region was developed to detect Y. pseudotuberculosis in milk powder. The DNA amplification could be completed in 1 h, and detected by produced white precipitate visible to naked eyes. The detection limit of LAMP assay was 10(0) fg/reaction for genomic DNA, and 10(0) CFU/100 g milk powder coupled with 12 h enrichment. LAMP assay is 100 times more sensitive than conventional polymerase chain reaction method for detecting Y. pseudotuberculosis, and correctly identified 18 cases of Y. pseudotuberculosis contaminations from 236 commercial milk powder products. In conclusion, the developed LAMP assay may facilitate rapid detection of Y. pseudotuberculosis contaminations in agricultural and food products.

PRACTICAL APPLICATION

Rapid and accurate detection of Yersinia pseudotuberculosis in milk products.

Reprint of New opportunities for improved ribotyping of C. difficile clinical isolates by exploring their genomes

Clostridium difficile causes outbreaks of infectious diarrhoea, most commonly occurring in healthcare institutions. Recently, concern has been raised with reports of C. difficile disease in those traditionally thought to be at low risk i.e. community acquired rather than healthcare acquired. This has increased awareness for the need to track outbreaks and PCR-ribotyping has found widespread use to elucidate epidemiologically linked isolates. PCR-ribotyping uses conserved regions of the 16S rRNA gene and 23S rRNA gene as primer binding sites to produce varying PCR products due to the intergenic spacer (ITS1) regions of the multiple operons. With the explosion of whole genome sequence data it became possible to analyse the start of the 23S rRNA gene for a more accurate selection of regions closer to the end of the ITS1. However the following questions must still be asked: (i) Does the chromosomal organisation of the rrn operon vary between C. difficile strains? and (ii) just how conserved are the primer binding regions? Eight published C. difficile genomes have been aligned to produce a detailed database of indels of the ITS1's from the rrn operon sets. An iPad Filemaker Go App has been constructed and named RiboTyping (RT). It contains detail such as sequences, ribotypes, strain numbers, GenBank numbers and genome position numbers. Access to various levels of the database is provided so that details can be printed. There are three main regions of the rrn operon that have been analysed by the database and related to each other by strain, ribotype and operon: (1) 16S gene (2) ITS1 indels (3) 23S gene. This has enabled direct intra- and inter-genomic comparisons at the strain, ribotype and operon (allele) levels in each of the three genomic regions. This is the first time that such an analysis has been done. By using the RT App with search criteria it will be possible to select probe combinations for specific strains/ribotypes/rrn operons for experiments to do with diagnostics, typing and recombination of operons. Many more incomplete C. difficile whole genome sequencing projects are recorded in GenBank as underway and the rrn operon information from these can also be added to the RT App when available. The RT App will help simplify probe selection because of the complexity of the ITS1 in C. difficile even in a single genome and because other allele-specific regions (16S and 23S genes) of variability can be relationally compared to design extra probes to increase sensitivity.

Challenges for standardization of Clostridium difficile typing methods

Typing of Clostridium difficile facilitates understanding of the epidemiology of the infection. Some evaluations have shown that certain strain types (for example, ribotype 027) are more virulent than others and are associated with worse clinical outcomes. Although restriction endonuclease analysis (REA) and pulsed-field gel electrophoresis have been widely used in the past, PCR ribotyping is the current method of choice for typing of C. difficile. However, global standardization of ribotyping results is urgently needed. Whole-genome sequencing of C. difficile has the potential to provide even greater epidemiologic information than ribotyping.

Intra- and intergenomic variation of ribosomal RNA operons in concurrent Alteromonas macleodii strains

Biodiversity estimates based on ribosomal operon sequence diversity rely on the premise that a sequence is characteristic of a single specific taxon or operational taxonomic unit (OTU). Here, we have studied the sequence diversity of 14 ribosomal RNA operons (rrn) contained in the genomes of two isolates (five operons in each genome) and four metagenomic fosmids, all from the same seawater sample. Complete sequencing of the isolate genomes and the fosmids establish that they represent strains of the same species, Alteromonas macleodii, with average nucleotide identity (ANI) values >97 %. Nonetheless, we observed high levels of intragenomic heterogeneity (i.e., variability between operons of a single genome) affecting multiple regions of the 16S and 23S rRNA genes as well as the internally transcribed spacer 1 (ITS-1) region. Furthermore, the ribosomal operons exhibited intergenomic heterogeneity (i.e., variability between operons located in separate genomes) in each of these regions, compounding the variability. Our data reveal the extensive heterogeneity observed in natural populations of A. macleodii at a single point in time and support the idea that distinct lineages of A. macleodii exist in the deep Mediterranean. These findings highlight the potential of rRNA fingerprinting methods to misrepresent species diversity while simultaneously failing to recognize the ecological significance of individual strains.

New opportunities for improved ribotyping of C. difficile clinical isolates by exploring their genomes

Clostridium difficile causes outbreaks of infectious diarrhoea, most commonly occurring in healthcare institutions. Recently, concern has been raised with reports of C. difficile disease in those traditionally thought to be at low risk i.e. community acquired rather than healthcare acquired. This has increased awareness for the need to track outbreaks and PCR-ribotyping has found widespread use to elucidate epidemiologically linked isolates. PCR-ribotyping uses conserved regions of the 16S rRNA gene and 23S rRNA gene as primer binding sites to produce varying PCR products due to the intergenic spacer (ITS1) regions of the multiple operons. With the explosion of whole genome sequence data it became possible to analyse the start of the 23S rRNA gene for a more accurate selection of regions closer to the end of the ITS1. However the following questions must still be asked: (i) Does the chromosomal organisation of the rrn operon vary between C. difficile strains? and (ii) just how conserved are the primer binding regions? Eight published C. difficile genomes have been aligned to produce a detailed database of indels of the ITS1's from the rrn operon sets. An iPad Filemaker Go App has been constructed and named RiboTyping (RT). It contains detail such as sequences, ribotypes, strain numbers, GenBank numbers and genome position numbers. Access to various levels of the database is provided so that details can be printed. There are three main regions of the rrn operon that have been analysed by the database and related to each other by strain, ribotype and operon: (1) 16S gene (2) ITS1 indels (3) 23S gene. This has enabled direct intra- and inter-genomic comparisons at the strain, ribotype and operon (allele) levels in each of the three genomic regions. This is the first time that such an analysis has been done. By using the RT App with search criteria it will be possible to select probe combinations for specific strains/ribotypes/rrn operons for experiments to do with diagnostics, typing and recombination of operons. Many more incomplete C. difficile whole genome sequencing projects are recorded in GenBank as underway and the rrn operon information from these can also be added to the RT App when available. The RT App will help simplify probe selection because of the complexity of the ITS1 in C. difficile even in a single genome and because other allele-specific regions (16S and 23S genes) of variability can be relationally compared to design extra probes to increase sensitivity.

Structure of rrn operons in pathogenic non-cultivable treponemes: sequence but not genomic position of intergenic spacers correlates with classification of Treponema pallidum and Treponema paraluiscuniculi strains

This study examined the sequences of the two rRNA (rrn) operons of pathogenic non-cultivable treponemes, comprising 11 strains of T. pallidum ssp. pallidum (TPA), five strains of T. pallidum ssp. pertenue (TPE), two strains of T. pallidum ssp. endemicum (TEN), a simian Fribourg-Blanc strain and a rabbit T. paraluiscuniculi (TPc) strain. PCR was used to determine the type of 16S–23S ribosomal intergenic spacers in the rrn operons from 30 clinical samples belonging to five different genotypes. When compared with the TPA strains, TPc Cuniculi A strain had a 17 bp deletion, and the TPE, TEN and Fribourg-Blanc isolates had a deletion of 33 bp. Other than these deletions, only 17 heterogeneous sites were found within the entire region (excluding the 16S–23S intergenic spacer region encoding tRNA-Ile or tRNA-Ala). The pattern of nucleotide changes in the rrn operons corresponded to the classification of treponemal strains, whilst two different rrn spacer patterns (Ile/Ala and Ala/Ile) appeared to be distributed randomly across species/subspecies classification, time and geographical source of the treponemal strains. It is suggested that the random distribution of tRNA genes is caused by reciprocal translocation between repetitive sequences mediated by a recBCD-like system.

Maternal oral origin of Fusobacterium nucleatum in adverse pregnancy outcomes as determined using the 16S-23S rRNA gene intergenic transcribed spacer region

Fusobacterium nucleatum, a common Gram-negative anaerobe prevalent in the oral cavity, possesses the ability to colonize the amniotic cavity and the fetus. However, F. nucleatum may also be part of the vaginal microbiota from where it could reach the amniotic tissues. Due to the heterogeneity of F. nucleatum, consisting of five subspecies, analysis at the subspecies/strain level is desirable to determine its precise origin. The aims of this study were: (i) to evaluate the use of the 16S-23S rRNA gene intergenic transcribed spacer (ITS) region as a tool to differentiate subspecies of F. nucleatum, and (ii) to design a simplified technique based on the ITS to determine the origin of F. nucleatum strains associated with adverse pregnancy outcomes. Amplified fragments of the 16S-23S rRNA gene ITS region corresponding to the five subspecies of F. nucleatum were subjected to cloning and sequencing to characterize the different ribosomal operons of the subspecies. Distinctive length and sequence patterns with potential to be used for identification of the subspecies/strain were identified. These were used to evaluate the origin of F. nucleatum identified in neonatal gastric aspirates (swallowed amniotic fluid) by sequence comparisons with the respective oral and vaginal maternal samples. A simplified technique using a strain-specific primer in a more sensitive nested PCR was subsequently developed to analyse ten paired neonatal-maternal samples. Analysing the variable fragment of the ITS region allowed the identification of F. nucleatum subsp. polymorphum from an oral origin as potentially being involved in neonatal infections. Using a strain-specific primer, the F. nucleatum subsp. polymorphum strain was detected in both neonatal gastric aspirates and maternal oral samples in cases of preterm birth from mothers presenting with localized periodontal pockets. Interestingly, the same strain was not present in the vaginal sample of any case investigated. The 16S-23S rRNA gene ITS can be a useful tool to determine the origin of F. nucleatum. The results of this study strongly indicate that F. nucleatum subsp. polymorphum of oral origin could be involved with pregnancy complications.

Prophage carriage and diversity within clinically relevant strains of Clostridium difficile

Prophages are encoded in most genomes of sequenced Clostridium difficile strains. They are key components of the mobile genetic elements and, as such, are likely to influence the biology of their host strains. The majority of these phages are not amenable to propagation, and therefore the development of a molecular marker is a useful tool with which to establish the extent and diversity of C. difficile prophage carriage within clinical strains. To design markers, several candidate genes were analyzed including structural and holin genes. The holin gene is the only gene present in all sequenced phage genomes, conserved at both terminals, with a variable mid-section. This allowed us to design two sets of degenerate PCR primers specific to C. difficile myoviruses and siphoviruses. Subsequent PCR analysis of 16 clinical C. difficile ribotypes showed that 15 of them are myovirus positive, and 2 of them are also siphovirus positive. Antibiotic induction and transmission electron microscope analysis confirmed the molecular prediction of myoviruses and/or siphovirus presence. Phylogenetic analysis of the holin sequences identified three groups of C. difficile phages, two within the myoviruses and a divergent siphovirus group. The marker also produced tight groups within temperate phages that infect other taxa, including Clostridium perfringens, Clostridium botulinum, and Bacillus spp., which suggests the potential application of the holin gene to study prophage carriage in other bacteria. This study reveals the high incidence of prophage carriage in clinically relevant strains of C. difficile and correlates the molecular data to the morphological observation.

Comparison of two capillary gel electrophoresis systems for Clostridium difficile ribotyping, using a panel of ribotype 027 isolates and whole-genome sequences as a reference standard

PCR ribotyping is the most commonly used Clostridium difficile genotyping method, but its utility is limited by lack of standardization. In this study, we analyzed four published whole genomes and tested an international collection of 21 well-characterized C. difficile ribotype 027 isolates as the basis for comparison of two capillary gel electrophoresis (CGE)-based ribotyping methods. There were unexpected differences between the 16S-23S rRNA intergenic spacer region (ISR) allelic profiles of the four ribotype 027 genomes, but six bands were identified in all four and a seventh in three genomes. All seven bands and another, not identified in any of the whole genomes, were found in all 21 isolates. We compared sequencer-based CGE (SCGE) with three different primer pairs to the Qiagen QIAxcel CGE (QCGE) platform. Deviations from individual reference/consensus band sizes were smaller for SCGE (0 to 0.2 bp) than for QCGE (4.2 to 9.5 bp). Compared with QCGE, SCGE more readily distinguished bands of similar length (more discriminatory), detected bands of larger size and lower intensity (more sensitive), and assigned band sizes more accurately and reproducibly, making it more suitable for standardization. Specifically, QCGE failed to identify the largest ISR amplicon. Based on several criteria, we recommend the primer set 16S-USA/23S-USA for use in a proposed standard SCGE method. Similar differences between SCGE and QCGE were found on testing of 14 isolates of four other C. difficile ribotypes. Based on our results, ISR profiles based on accurate sequencer-based band lengths would be preferable to agarose gel-based banding patterns for the assignment of ribotypes.

Diversity of 16S-23S rDNA internal transcribed spacer (ITS) reveals phylogenetic relationships in Burkholderia pseudomallei and its near-neighbors

Length polymorphisms within the 16S-23S ribosomal DNA internal transcribed spacer (ITS) have been described as stable genetic markers for studying bacterial phylogenetics. In this study, we used these genetic markers to investigate phylogenetic relationships in Burkholderia pseudomallei and its near-relative species. B. pseudomallei is known as one of the most genetically recombined bacterial species. In silico analysis of multiple B. pseudomallei genomes revealed approximately four homologous rRNA operons and ITS length polymorphisms therein. We characterized ITS distribution using PCR and analyzed via a high-throughput capillary electrophoresis in 1,191 B. pseudomallei strains. Three major ITS types were identified, two of which were commonly found in most B. pseudomallei strains from the endemic areas, whereas the third one was significantly correlated with worldwide sporadic strains. Interestingly, mixtures of the two common ITS types were observed within the same strains, and at a greater incidence in Thailand than Australia suggesting that genetic recombination causes the ITS variation within species, with greater recombination frequency in Thailand. In addition, the B. mallei ITS type was common to B. pseudomallei, providing further support that B. mallei is a clone of B. pseudomallei. Other B. pseudomallei near-neighbors possessed unique and monomorphic ITS types. Our data shed light on evolutionary patterns of B. pseudomallei and its near relative species.

Toxins and Antibiotic Resistance in Staphylococcus aureus Isolated from a Major Hospital in Lebanon

Molecular characterization of Staphylococcus aureus is of both clinical and infection control importance. Virulence determinants using PCR and multiple drug resistance profiles were studied in 130 S. aureus isolates. PCR-RFLP analysis of the 16S-23S DNA spacer region was done to investigate the level of 16S-23S ITS (internal transcribed spacer) polymorphism. Methicillin-resistant S. aureus (MRSA), which represented 72% of the studied isolates, showed multiple drug resistance with 18% being resistant to 10-18 of the drugs used compared to a maximum resistance to 9 antibiotics with the methicillin sensitive S. aureus (MSSA) isolates. Exfoliative toxin A (ETA) was more prevalent than B (ETB) with virulent determinants being additionally detected in multiple drug-resistant isolates. 16S-23S ITS PCR-RFLP combined with sequencing of the primary product was successful in generating molecular fingerprints of S. aureus and could be used for preliminary typing. This is the first study to demonstrate the incidence of virulent genes, ACME, and genetic diversity of S. aureus isolates in Lebanon. The data presented here epitomize a starting point defining the major genetic populations of both MRSA and MSSA in Lebanon and provide a basis for clinical epidemiological studies.

Multidrug resistance in European Clostridium difficile clinical isolates

OBJECTIVES

Multidrug resistance and antibiotic resistance mechanisms were investigated in 316 Clostridium difficile clinical isolates collected during the first European surveillance on C. difficile in 2005.

METHODS

MICs of eight different antibiotics were determined using Etest. Reserpine- and carbonyl cyanide m-chlorophenylhydrazone-sensitive efflux was tested using the agar dilution method. Molecular analysis of the resistance mechanisms was performed using PCR assays, PCR mapping and sequencing.

RESULTS

One hundred and forty-eight C. difficile strains were resistant to at least one antibiotic and 82 (55%) were multidrug resistant. In particular, 48% of these isolates were resistant to erythromycin, clindamycin, moxifloxacin and rifampicin. New genetic elements or determinants conferring resistance to erythromycin/clindamycin or tetracycline were identified. Even if most multiresistant strains carried an erm(B) gene, quite a few were erm(B) negative. In-depth analysis of the underlying mechanism in these isolates was carried out, including analysis of 23S rDNA and the ribosomal proteins L4 and L22. Interestingly, resistance to rifampicin was observed in multidrug-resistant strains in association with resistance to fluoroquinolones. Mutations in the rpo(B) and gyrA genes were identified as the cause of resistance to these antibiotics, respectively.

CONCLUSIONS

Characterization of multidrug-resistant C. difficile clinical isolates shows that antibiotic resistance is changing, involving new determinants and mechanisms and providing this pathogen with potential advantages over the co-resident gut flora. The present paper provides, for the first time, a comprehensive picture of the different characteristics of multidrug-resistant C. difficile strains in Europe in 2005 and represents an important source of data for future comparative European studies.

Mechanisms behind variation in the Clostridium difficile 16S-23S rRNA intergenic spacer region

Clostridium difficile infection is an increasing problem in hospitals worldwide, mainly due to the recent emergence of a hypervirulent C. difficile strain. C. difficile PCR ribotyping, based on size variation of the 16S-23S rRNA intergenic spacer region (16S-23S ISR), is widely used in Europe for molecular epidemiological investigation. The mechanism underlying the 16S-23S ISR size variations in the genome of C. difficile is currently not completely understood. To elucidate this mechanism, isolates of six different PCR ribotypes were analysed by cloning and sequencing the 16S-23S ISR. A direct repeat, IB, of 9 bp was detected up to five times in the 16S-23S ISR in all 47 clones investigated. Thirty-five clones displayed differences either by ribotype or by nucleotide sequence. The sequences of the 16S-23S ISR of C. difficile showed a uniformly organized structure, composed of a tRNA(Ala) gene and spacers of 33 and 53 bp separated by the 9 bp direct repeat IB. The results of the study support the hypothesis that this composition is responsible for the length variations seen in the 16S-23S ISR, and indicate that these length variations result from slipped-strand mispairing and intra- and possibly interchromosomal homologous recombination.

Applicability of the 16S-23S rDNA internal spacer for PCR detection of the phytostimulatory PGPR inoculant Azospirillum lipoferum CRT1 in field soil

AIMS

To assess the applicability of the 16S-23S rDNA internal spacer regions (ISR) as targets for PCR detection of Azospirillum ssp. and the phytostimulatory plant growth-promoting rhizobacteria seed inoculant Azospirillum lipoferum CRT1 in soil.

METHODS AND RESULTS

Primer sets were designed after sequence analysis of the ISR of A. lipoferum CRT1 and Azospirillum brasilense Sp245. The primers fAZO/rAZO targeting the Azospirillum genus successfully yielded PCR amplicons (400-550 bp) from Azospirillum strains but also from certain non-Azospirillum strains in vitro, therefore they were not appropriate to monitor indigenous Azospirillum soil populations. The primers fCRT1/rCRT1 targeting A. lipoferum CRT1 generated a single 249-bp PCR product but could also amplify other strains from the same species. However, with DNA extracts from the rhizosphere of field-grown maize, both fAZO/rAZO and fCRT1/rCRT1 primer sets could be used to evidence strain CRT1 in inoculated plants by nested PCR, after a first ISR amplification with universal ribosomal primers. In soil, a 7-log dynamic range of detection (10(2)-10(8) CFU g(-1) soil) was obtained.

CONCLUSIONS

The PCR primers targeting 16S-23S rDNA ISR sequences enabled detection of the inoculant A. lipoferum CRT1 in field soil.

SIGNIFICANCE AND IMPACT OF THE STUDY

Convenient methods to monitor Azospirillum phytostimulators in the soil are lacking. The PCR protocols designed based on ISR sequences will be useful for detection of the crop inoculant A. lipoferum CRT1 under field conditions.

Reverse line blot hybridization and DNA sequencing studies of the 16S-23S rRNA gene intergenic spacer regions of five emerging pathogenic Nocardia species

The objective of this study was to examine DNA sequence polymorphisms in the 16S-23S rRNA gene intergenic spacer (ITS) regions of five emerging pathogenic Nocardia species: Nocardia beijingensis, Nocardia blacklockiae, Nocardia thailandica, Nocardia elegans and Nocardia vinacea. A set of six isolates belonging to the species of interest and 135 isolates belonging to other Nocardia species was studied. A PCR-based reverse line blot (RLB) hybridization assay incorporating species- or intraspecies ITS rRNA gene operon-specific probes was then developed for species identification. Substantial intraspecies sequence variation among different ITS operons was identified. Four sequence types of N. thailandica, eight sequence types of N. beijingensis (four types for each of two strains) and five sequence types of N. blacklockiae, N. elegans and N. vinacea were found. The results represent the first evidence of ITS sequence heterogeneity in emerging species of Nocardia. By incorporating species/operon-specific probes into a RLB assay, unique RLB patterns were identified for each of the species and every sequence type. The PCR/RLB assay demonstrated high specificity and showed promise in both the identification and genotyping of Nocardia species. More detailed studies of the polymorphism within the ITS locus may further advance our capacity to reliably identify and subtype medically important Nocardia species.

The vexed relationship between Clostridium difficile and inflammatory bowel disease: an assessment of carriage in an outpatient setting among patients in remission

OBJECTIVES

Comorbidity with Clostridium difficile may cause diagnostic delay in newly presenting inflammatory bowel disease (IBD) patients, trigger relapse in established disease, confound therapies, and serve as an indicator of an underlying defect in innate immunity. Retrospective analyses have suggested community acquisition; to address this we conducted a prospective analysis of C. difficile carriage in IBD patients using molecular methods specifically in an outpatient setting.

METHODS

Recruited participants had long-standing diagnoses of ulcerative colitis (n = 64) and Crohn's disease (n = 58), were in clinical remission, and had no recent exposure to antibiotics, corticosteroids, immunomodulatory drugs or recent hospitalization. Isolates were cultured from stools and confirmed by 16S sequencing. The antibiotic susceptibilities of the isolates were tested followed by further strain characterization by toxinotyping, ribotyping, and pulsed-field gel electrophoresis (PFGE).

RESULTS

The frequency of toxigenic C. difficile was higher in IBD patients than in healthy volunteers at 8.2 and 1.0%, respectively (P = 0.02 Fisher's exact test). All strains belonged to toxinotype 0 with rare subtypes of this group noted in five isolates and represented by an altered repressor genotype. Patients harbored a diverse range of toxigenic ribotype groups, including those previously associated with C. difficile-associated disease (CDAD) (R015, R005, and R020) and the rarer types R062, R050, and R003. Interestingly, common nosocomial groups were not identified. The considerable nonclonal distribution of distinct strains was further demonstrated by PFGE genomic fingerprinting. None of the study subjects experienced a clinical episode of CDAD during a 6-month period of follow-up.

CONCLUSIONS

Detection of C. difficile is increased in IBD outpatients in remission, and strain diversity is consistent with community acquisition from a multitude of sources.

16S-23S rDNA internal transcribed spacer regions in four Proteus species

Proteus is a Gram-negative, facultative anaerobic bacterium. In this study, 813 Proteus 16S-23S rDNA internal transcribed spacer (ITS) sequences were determined from 46 Proteus strains, including 388 ITS from 22 P. mirabilis strains, 211 ITS from 12 P. vulgaris strains, 169 ITS from 10 P. penneri strains, and 45 ITS from 2 P. myxofaciens strains. The Proteus strains carry mainly two types of ITS, ITS(Glu) (containing tRNA(Glu (UUC)) gene) and ITS(Ile+Ala) (containing tRNA(Ile (GAU)) and tRNA(Ala (UGC)) gene), and are in the forms of 28 variants with 25 genomic origins. The ITS sequences are a mosaic-like structure consisting of three conservative regions and two variable regions. The nucleotide identity of ITS subtypes in strains of the same species ranges from 96.2% to 100%. The divergence of Proteus ITS divergence was most likely due to intraspecies recombinations or horizontal transfers of sequence blocks. The phylogenetic relationship deduced from the second variable region of ITS sequences of the three facultative human pathogenic species P. mirabilis, P. vulgaris and P. penneri is similar with that based on 16S rDNA sequences, but has higher resolution to differentiate closely related P. vulgaris and P. penneri. This study is the first comprehensive study of ITS in four Proteus species and laid solid foundation for the development of high-throughput technology for quick and accurate identification of the important foodborne and nosocomial pathogens.

Unusual intragenomic and interspecific variability of Geobacillus 16S-23S rRNA internal transcribed spacers

The aim of this study was to evaluate the inter-and intraspecific as well as intragenomic variability of Geobacillus 16S–23S rRNA internal transcribed spacers without tRNA genes and to compare these sequences with sequences bearing tRNA genes. In this study the structural analysis was performed in a unique way because the length and the sequence of the structural blocks were adjusted to fit the structure of 16S–23S rRNA internal transcribed spacers of five different Geobacillus species. Our study demonstrated the mosaic-like structure of 16S–23S rRNA internal transcribed spacers in Geobacillus. Some characteristics of these spacers of geobacilli were not previously reported for other bacteria: unusually short conserved sequence in the 5′ end region, some identical conserved blocks in both 5′ and 3′ regions of 16S–23S rRNA internal transcribed spacers, the same sequence blocks in both 16S–23S and 23S–5S rRNA intergenic spacers. Our study demonstrated quite uniform arrangement of the sequence blocks in Geobacillus thermodenitrificans. This species diverged early in the phylogenetic tree of the genus Geobacillus. For the phylogenetically recent species Geobacillus kaustophilus and Geobacillus lituanicus the low inter-and intraspecific, but high intragenomic variability, as a consequence of recent phylogenetic events, was established.

Characterization of Clostridium difficile isolates using capillary gel electrophoresis-based PCR ribotyping

We have developed a Clostridium difficile PCR ribotyping method based on capillary gel electrophoresis and have compared it with conventional PCR ribotyping. A total of 146 C. difficile isolates were studied: five isolates were reference strains (PCR ribotypes 001, 014, 017, 027 and 053); 141 were clinical isolates comprising 39 Austrian PCR ribotypes collected in the period 2006–2007 at 25 Austrian healthcare facilities. Capillary gel electrophoresis yielded up to 11 fragments per isolate and 47 ribotype patterns. All but one of the five PCR ribotypes of reference strains were clearly reflected in the chromatograms of capillary-based typing. Capillary gel electrophoresis divided 24 isolates belonging to PCR ribotype type 014 into seven subgroups, whereas subtyping the same isolates using multiple-locus variable-number tandem-repeat analysis yielded three unrelated subgroups, without obvious correlation to sr subgroups. Using a web-based software program (http://webribo.ages.at), we were able to correctly identify these 014 isolates by simply allocating the seven subgroup patterns to one ribotype, i.e. to PCR ribotype 014. We consider capillary gel electrophoresis-based PCR ribotyping to be a way of overcoming the problems associated with inter-laboratory comparisons of typing results, while at the same time substantially diminishing the hands-on time for PCR ribotyping.