Jean-Christophe Lagier1, Fehmida Bibi2, Dhamodharan Ramasamy1, Esam I Azhar3, Catherine Robert1, Muhammad Yasir2, Asif A Jiman-Fatani4, Khalid Z Alshali5, Pierre-Edouard Fournier1, Didier Raoult6. 1. Unité de Recherche sur les Maladies Infectieuses et Tropicales Emergentes, UMR CNRS 7278 - IRD 198, Institut Hospitalo-Universitaire Méditerranée-Infection, Faculté de médecine, Aix-Marseille Université, France. 2. Special Infectious Agents Unit, King Fahd Medical Research Center, King Abdulaziz University, Jeddah, Saudi Arabia. 3. Special Infectious Agents Unit, King Fahd Medical Research Center, King Abdulaziz University, Jeddah, Saudi Arabia ; Department of Medical Laboratory Technology, Faculty of Applied Medical Sciences, King Abdulaziz University, Jeddah, Saudi Arabia. 4. Department of Medical Microbiology and Parasitology, Faculty of Medicine, King. 5. Department of Medicine, Faculty of Medicine, King Abdulaziz University, PO Box 80215, Jeddah, Postal code 21589, Saudi Arabia. 6. Unité de Recherche sur les Maladies Infectieuses et Tropicales Emergentes, UMR CNRS 7278 - IRD 198, Institut Hospitalo-Universitaire Méditerranée-Infection, Faculté de médecine, Aix-Marseille Université, France ; Special Infectious Agents Unit, King Fahd Medical Research Center, King Abdulaziz University, Jeddah, Saudi Arabia.
Abstract
Clostridium jeddahense strain JCD(T) (= CSUR P693 = DSM 27834) is the type strain of C. jeddahense sp. nov. This strain, whose genome is described here, was isolated from the fecal flora of an obese 24 year-old Saudian male (BMI=52 kg/m(2)). Clostridium jeddahense strain JCD(T) is an obligate Gram-positive bacillus. Here we describe the features of this organism, together with the complete genome sequence and annotation. The 3,613,503 bp long genome (1 chromosome, no plasmid) exhibits a G+C content of 51.95% and contains 3,462 protein-coding and 53 RNA genes, including 4 rRNA genes.
Clostridium jeddahense strain JCD(T) (= CSUR P693 = DSM 27834) is the type strain of C. jeddahense sp. nov. This strain, whose genome is described here, was isolated from the fecal flora of an obese 24 year-old Saudian male (BMI=52 kg/m(2)). Clostridium jeddahense strain JCD(T) is an obligate Gram-positive bacillus. Here we describe the features of this organism, together with the complete genome sequence and annotation. The 3,613,503 bp long genome (1 chromosome, no plasmid) exhibits a G+C content of 51.95% and contains 3,462 protein-coding and 53 RNA genes, including 4 rRNA genes.
strain JCDT (=CSUR P693 = DSM 27834), is the type strain of sp. nov. This bacterium is a Gram-positive, anaerobic, spore-forming indole, positive bacillus that was isolated from the stool of an obese 24 year-old Saudian individual, as a part of a culturomics study as previously reported.The usual parameters used to delineate a bacterial species include 16S rDNA sequence identity and phylogeny [1,2], genomic G + C content diversity, and DNA–DNA hybridization (DDH) [3,4]. Nevertheless, some limitations appeared notably because the cutoff values vary dramatically between species and genera [5]. The introduction of high-throughput sequencing techniques made genomic data for many bacterial species available [6]. We recently proposed a new method (taxono-genomics), which includes genomic data in a polyphasic approach to describe new bacterial species [6]. This strategy combines phenotypic characteristics, including MALDI-TOF MS spectrum, and genomic analysis [7-37].Here, we present a summary classification and a set of features for sp. nov. strain JCDT (=CSUR P693 = DSM 27834), together with the description of the complete genome sequencing and annotation. These characteristics support the circumscription of the speciesThe genus was created in 1880 [38] and consists of obligate anaerobic rod-shaped bacilli able to produce endospores [38]. More than 200 species have been described to date (http://www.bacterio.cict.fr/c/clostridium.html). Members of the genus are mostly environmental bacteria or associated with the commensal digestive flora of mammals. However, several are major human pathogens, including , and [38].
Classification and features
A stool sample was collected from an obese 24-year-old male Saudian volunteer patient living in Jeddah. The patient gave an informed and signed consent, and the agreement of the Ethical Committee of the King Abdulaziz University, King Fahd medical Research Centre, Saudi Arabia, and the local ethics committee of the IFR48 (Marseille, France) were obtained under agreement number 014-CEGMR-2-ETH-P and 09-022 respectively. The fecal specimen was preserved at -80°C after collection and sent to Marseille. Strain JCDT (Table 1) was isolated in July 2013 by anaerobic cultivation on 5% sheep blood-enriched Columbia agar (BioMerieux, Marcy l’Etoile, France) after a 5-day preincubation on blood culture bottle with rumen fluid. This strain exhibited a 97.3% nucleotide sequence similarity with strain DSM 1294 (Figure 1). This value was lower than the 98.7% 16S rRNA gene sequence similarity threshold recommended by Stackebrandt and Ebers to delineate a new species without carrying out DNA-DNA hybridization [2] and was in the 78. 4 to 98.9% range of 16S rRNA identity values observed among 41 species with validly published names [52].
Table 1
Classification and general features of strain JCDT according to the MIGS recommendations [39]
MIGS ID
Property
Term
Evidence codea
Current classification
Domain Bacteria
TAS [40]
Phylum Firmicutes
TAS [41-43]
Class Clostridia
TAS [44,45]
Order Clostridiales
TAS [46,47]
Family Clostridiaceae
TAS [46,48]
Genus Clostridium
IDA [46,49,50]
Species Clostridium jeddahense
IDA
Type strain JCDT
IDA
Gram stain
Positive
IDA
Cell shape
Rod
IDA
Motility
Motile
IDA
Sporulation
Sporulating
IDA
Temperature range
Mesophile
IDA
Optimum temperature
37°C
IDA
MIGS-6.3
Salinity
Unknown
IDA
MIGS-22
Oxygen requirement
Anaerobic
IDA
Carbon source
Unknown
IDA
Energy source
Unknown
IDA
MIGS-6
Habitat
Human gut
IDA
MIGS-15
Biotic relationship
Free living
IDA
MIGS-14
Pathogenicity Biosafety level Isolation
Unknown 2 Human feces
MIGS-4
Geographic location
Jeddah, Saudi Arabia
IDA
MIGS-5
Sample collection time
July 2013
IDA
MIGS-4.1
Latitude
21.422487
IDA
MIGS-4.1
Longitude
39.856184
IDA
MIGS-4.3
Depth
Surface
IDA
MIGS-4.4
Altitude
0 m above sea level
IDA
Evidence codes - IDA: Inferred from Direct Assay; TAS: Traceable Author Statement (i.e., a direct report exists in the literature); NAS: Non-traceable Author Statement (i.e., not directly observed for the living, isolated sample, but based on a generally accepted property for the species, or anecdotal evidence). These evidence codes are from the Gene Ontology project [51]. If the evidence is IDA, then the property was directly observed for a live isolate by one of the authors or an expert mentioned in the acknowledgements.
Figure 1
A consensus phylogenetic tree highlighting the position of strain JCDT relative to other type strains within the Clostridum genus. GenBank accession numbers are indicated in parentheses. Sequences were aligned using CLUSTALW, and phylogenetic inferences were obtained using the maximum-likelihood method in the MEGA software package. Numbers at the nodes are the percentages of bootstrap values from 500 replicates that support the node. was used as outgroup. The scale bar represents a 2% nucleotide sequence divergence.
Evidence codes - IDA: Inferred from Direct Assay; TAS: Traceable Author Statement (i.e., a direct report exists in the literature); NAS: Non-traceable Author Statement (i.e., not directly observed for the living, isolated sample, but based on a generally accepted property for the species, or anecdotal evidence). These evidence codes are from the Gene Ontology project [51]. If the evidence is IDA, then the property was directly observed for a live isolate by one of the authors or an expert mentioned in the acknowledgements.A consensus phylogenetic tree highlighting the position of strain JCDT relative to other type strains within the Clostridum genus. GenBank accession numbers are indicated in parentheses. Sequences were aligned using CLUSTALW, and phylogenetic inferences were obtained using the maximum-likelihood method in the MEGA software package. Numbers at the nodes are the percentages of bootstrap values from 500 replicates that support the node. was used as outgroup. The scale bar represents a 2% nucleotide sequence divergence.Four growth temperatures (25, 30, 37, 45°C) were tested; growth occurred between 25 and 37°C, but optimal growth was observed at 37°C, 24 hours after inoculation. No growth occurred at 45°C. Colonies were translucent and approximately 0.2 to 0.3 mm in diameter on 5% sheep blood-enriched Columbia agar (BioMerieux). Growth of the strain was tested on the same agar under anaerobic and microaerophilic conditions using GENbag anaer and GENbag microaer systems, respectively (BioMerieux), and in aerobic conditions, with or without 5% CO2. Growth was observed only anaerobically. No growth occurred in aerobic or microaerophilic conditions. Gram staining showed Gram-positive rods able to form spores (Figure 2). A motility test was positive. Cells grown on agar exhibit a mean diameter of 1 µm and a mean length of 1.22 µm in electron microscopy (Figure 3).
Figure 2
Gram stain of strain JCDT
Figure 3
Transmission electron micrograph of strain JCDT, taken using a Morgani 268D (Philips) at an operating voltage of 60kV.The scale bar represents 500 nm.
Gram stain of strain JCDTTransmission electron micrograph of strain JCDT, taken using a Morgani 268D (Philips) at an operating voltage of 60kV.The scale bar represents 500 nm.Strain JCDT exhibited neither catalase nor oxidase activity (Table 2). Using an API Rapid ID 32A strip (BioMerieux), positive reactions were obtained for indole production, alkaline phosphatase, arginine arylamidase, proline arylamidase, alanine arylamidase, glycine arylamidase, histidine arylamidase, glutamyl glutamic acid arylamidase and serine arylamidase. Negative reactions were obtained for arginine dihydrolase, α-galactosidase, β-galactosidase, α-glucosidase, β-glucosidase, α-arabinosidase, N-acetyl-β-glucosaminidase, glutamic acid decarboxylase, α-fucosidase, nitrate reduction, leucyl glycine arylamidase, fermentation of mannose and raffinose, urease, β-galactosidase-6-phosphatase, β-glucuronidase, phenylalanine arylamidase, leucine arylamidase, pyroglutamic acid arylamidase and tyrosine arylamidase. Using an API 50CH strip (Biomerieux), strain JCDT was asaccharolytic.
Table 2
Differential characteristics of JCDT, JC122 [11], FF1 [34], NCIMB 8052, B1, H10, DSM 753 and DSM 1294 [53]†.
Properties
C.jeddahense
C.sporosphaeroides
C.cellulolyticum
C.leptum
C.senegalense
C.dakarense
C.beijerinckii
C.difficile
Cell diameter (µm)
1.0
0.5-0.6
1.5
0.6-0.8
1.1
1.2
1.7
3.0
Oxygen requirement
Strictly anaerobic
Strictly anaerobic
Strictly anaerobic
Strictly anaerobic
Strictly anaerobic
Strictly anaerobic
Strictly anaerobic
Strictly anaerobic
Gram stain
Positive
Positive
Positive
Positive
Positive
Positive
Variable
Variable
Motility
Motile
Non Motile
Motile
Non Motile
Motile
Motile
Motile
Motile
Endospore formation
+
+
+
+
+
+
+
+
Indole
+
-
-
Na
-
+
Na
Na
Production of
Alkaline phosphatase
+
Na
Na
Na
-
+
Na
Na
Catalase
-
Na
-
Na
-
-
-
Na
Oxidase
-
Na
Na
Na
-
-
Na
Na
Nitrate reductase
-
-
-
Na
-
-
-
-
Urease
-
-
-
Na
-
-
-
Na
β-galactosidase
-
Na
Na
Na
-
-
Na
Na
N-acetyl-glucosamine
-
Na
Na
Na
+
Na
Na
Acid from
L-Arabinose
-
-
Na
Na
Na
-
+
-
Ribose
-
-
Na
W
Na
-
-
-
Mannose
-
-
Na
-
Na
-
+
+
Mannitol
-
-
Na
-
Na
-
+
+
Sucrose
-
-Na
Na
W
Na
-
+
+
D-glucose
-
Na
Na
Na
Na
+
+
Na
D-fructose
-
Na
Na
Na
Na
-
+
+
D-maltose
-
Na
Na
Na
Na
+
+
-
D-lactose
-
Na
Na
Na
Na
-
+
-
G+C content (%)
52
41
51
26.8
27.98
28
28
Habitat
Human gut
Environment
Compost
Human gut
Human gut
Human gut
Human gut
Human gut
†(Na = data not available; w = weak, v = variable reaction)
†(Na = data not available; w = weak, v = variable reaction)is susceptible to amoxicillin, amoxicillin-clavulanate, imipenem, metronidazole, doxycycline, rifampicin, vancomycin but resistant to ceftriaxone, ciprofloxacin and trimethoprim-sulfamethoxazole. The comparisons with other species are summarized in Table 2.Matrix-assisted laser-desorption/ionization time-of-flight (MALDI-TOF) MS protein analysis was carried out as previously described [54]. Briefly, a pipette tip was used to pick one isolated bacterial colony from a culture agar plate and spread it as a thin film on a MTP 384 MALDI-TOF target plate (Bruker Daltonics, Leipzig, Germany). Twelve distinct deposits from twelve isolated colonies were performed for strain JCDT. Each smear was overlaid with 2 µL of matrix solution (saturated solution of alpha-cyano-4-hydroxycinnamic acid) in 50% acetonitrile, 2.5% tri-fluoracetic acid, and allowed to dry for 5 minutes. Measurements were performed with a Microflex spectrometer (Bruker). Spectra were recorded in the positive linear mode for the mass range of 2,000 to 20,000 Da (parameter settings: ion source 1 (ISI), 20kV; IS2, 18.5 kV; lens, 7 kV). A spectrum was obtained after 675 shots with variable laser power. The time of acquisition was between 30 seconds and 1 minute per spot. The twelve JCDT spectra were imported into the MALDI BioTyper software (version 2.0, Bruker) and analyzed by standard pattern matching (with default parameter settings) against the main spectra of 3,769 bacteria, including 228 spectra from 96 species. The method of identification included the m/z from 3,000 to 15,000 Da. For every spectrum, a maximum of 100 peaks were compared with spectra in database. The resulting score enabled the identification of tested species, or not: a score ≥ 2 with a validly published species enabled identification at the species level, a score ≥ 1.7 but < 2 enabled identification at the genus level, and a score < 1.7 did not enable any identification. No significant MALDI-TOF score was obtained for strain JCDT against the Bruker database, suggesting that our isolate was not a member of a known species. We added the spectrum from strain JCDT to our database (Figure 4). Finally, the gel view showed the spectral differences with other members of the genus (Figure 5).
Figure 4
Reference mass spectrum from strain JCDT. Spectra from 12 individual colonies were compared and a reference spectrum was generated.
Figure 5
Gel view comparing strain JCDT to other species. The gel view displays the raw spectra of loaded spectrum files arranged as a pseudo-electrophoretic gel. The x-axis records the m/z value. The left y-axis displays the running spectrum number originating from subsequent spectra loading. The peak intensity is expressed by a grey scale scheme code. The grey scale bar on the right y-axis indicates the relation between the shade of grey a peak is displayed with and the peak intensity in arbitrary units. Species names are shown on the left.
Reference mass spectrum from strain JCDT. Spectra from 12 individual colonies were compared and a reference spectrum was generated.Gel view comparing strain JCDT to other species. The gel view displays the raw spectra of loaded spectrum files arranged as a pseudo-electrophoretic gel. The x-axis records the m/z value. The left y-axis displays the running spectrum number originating from subsequent spectra loading. The peak intensity is expressed by a grey scale scheme code. The grey scale bar on the right y-axis indicates the relation between the shade of grey a peak is displayed with and the peak intensity in arbitrary units. Species names are shown on the left.
Genome sequencing information
Genome project history
The organism was selected for sequencing on the basis of its phylogenetic position and 16S rDNA similarity to members of the genus , and is part of a study of the human digestive flora aiming at isolating all bacterial species in human feces [55]. It was the 101st genome of a species and the first genome of sp. nov. The GenBank accession number is CBYL00000000. The assembly consists of 104 contigs. Table 3 shows the project information and its association with MIGS version 2.0 compliance [39].
Table 3
Project information
MIGS ID
Property
Term
MIGS-31
Finishing quality
High-quality draft
MIGS-28
Libraries used
Paired end and Mate pair
MIGS-29
Sequencing platform
MySeq Illumina
MIGS-31.2
Fold coverage
94.91×
MIGS-30
Assemblers
Newbler
MIGS-32
Gene calling method
PRODIGAL
Genbank Date of Release
February 12, 2014
Genbank project ID
CBYL00000000
MIGS-13
Project relevance
Study of the human gut microbiome
Growth conditions and DNA isolation
sp. nov., strain JCDT (= CSUR P693 = DSM 27834) was grown on 5% sheep blood-enriched Columbia agar (BioMerieux) at 37°C in anaerobic atmosphere. Bacteria grown on three Petri dishes were harvested and resuspended in 4x100µL of TE buffer. Then, 200 µL of this suspension was diluted in 1ml TE buffer for lysis treatment that included a 30- minute incubation with 2.5 µg/µL lysozyme at 37°C, followed by an overnight incubation with 20 µg/µL proteinase K at 37°C. Extracted DNA was then purified using 3 successive phenol-chloroform extractions and ethanol precipitation at -20°C overnight. After centrifugation, the DNA was resuspended in 160 µL TE buffer.
Genome sequencing and assembly
Genomic DNA of was sequenced on a MiSeq sequencer (Illumina, Inc, San Diego CA 92121, USA) with 2 applications: paired end and mate pair. The paired end and the mate pair strategies were barcoded in order to be mixed respectively with 14 other genomic projects constructed according the Nextera XT library kit (Illumina) and 11 others projects with the nextera Mate pair kit (Illumina).The gDNA was quantified by a Qubit assay with the high sensitivity kit (Life technologies, Carlsbad, CA, USA) to 11.1 ng/µL and dilution was performed such that 1ng of each strain’s gDNA was used to construct the paired end library. The “tagmentation” step fragmented and tagged the DNA .Then limited cycle PCR amplification completed the tag adapters and introduced dual-index barcodes. After purification on Ampure beads (Life Technolgies, Carlsbad, CA, USA), the libraries were normalized on specific beads according to the Nextera XT protocol (Illumina). Normalized libraries are pooled into a single library for sequencing on the MiSeq. The pooled single strand library was loaded onto the reagent cartridge and then onto the instrument along with the flow cell. Automated cluster generation and paired-end sequencing with dual index reads was performed in a single 39-hour run at a 2x250 bp read length. Within this pooled run, the index representation was determined to be 7.3%. Total information of 5.3 Gbases was obtained from a 574 K/mm2 density with 95.4% (11,188,000 clusters) of the clusters passing quality control (QC) filters. From the genome sequencing process, the 753,292 produced Illumina reads for were filtered according to the read qualities.The mate pair library was constructed from 1 µg of genomic DNA using the Nextera Mate Pair Illumina guide. The genomic DNA sample is simultaneously fragmented and tagged with a mate pair junction adapter. The profile of the fragmentation was validated on an Agilent 2100 BioAnalyzer (Agilent Technologies, Inc., Santa Clara, CA, USA) with a DNA7500 labchip. The DNA fragments range in size from 1 kb up to 11 kb with a mean size of 7kb. No size selection was performed and 600 ng tagmented fragments were circularized. The larger circularized DNA molecules were physically sheared to smaller sized fragments with a mean size of 625 bp on the Covaris device S2 in microtubes (Woburn, MA, USA) .The library’s profile and the quantitation were visualized on a High Sensitivity Bioanalyzer LabChip. The libraries were normalized to 2 nM and pooled. After a denaturation step and dilution at 10 pM the pool of libraries was loaded onto the reagent cartridge and then onto the instrument along with the flow cell. Automated cluster generation and sequencing run was performed in a single 39-hour run at a 2x250 bp read length.Total information of 3.9 Gb was obtained from a 399 K/mm2 density with 97.9% (7,840,000 clusters) of the clusters passing quality control (QC) filters. Within this pooled run, the index representation for was determined to be 6.54%.From this genome sequencing process, the 501,426 produced Illumina reads for were filtered according to the read qualities.
Genome annotation
Open Reading Frames (ORFs) were predicted using Prodigal [56] with default parameters. However, the predicted ORFs were excluded if they spanned a sequencing gap region. The predicted bacterial protein sequences were searched against the GenBank [57] and Clusters of Orthologous Groups (COG) databases using BLASTP. The tRNAs and rRNAs were predicted using the tRNAScan-SE [58] and RNAmmer [59] tools, respectively. Signal peptides and numbers of transmembrane helices were predicted using SignalP [60] and TMHMM [61], respectively. Mobile genetic elements were predicted using PHAST [62] and RAST [63]. ORFans were identified if their BLASTP E-value was lower than 1e-03 for alignment length greater than 80 amino acids. If alignment lengths were smaller than 80 amino acids, we used an E-value of 1e-05. Such parameter thresholds have already been used in previous work to define ORFans. Artemis [64] and DNA Plotter [65] were used for data management and visualization of genomic features, respectively. The Mauve alignment tool (version 2.3.1) was used for multiple genomic sequence alignment [66].To estimate the mean level of nucleotide sequence similarity at the genome level between and 7 other members of the genus , we used the Average Genomic Identity Of gene Sequences (AGIOS) home-made software [6]. Briefly, this software combines the Proteinortho software [67] for detecting orthologous proteins between pairs of genomes, then retrieves the corresponding genes and determines the mean percentage of nucleotide sequence identity among orthologous ORFs using the Needleman-Wunsch global alignment algorithm. strain JCDT was compared to C. senegalense strain JC122, strain FF1, strain NCIMB 8052, strain B1, strain H10, strain DSM 753, and strain DSM 1294 (see Table 6B).
Table 6B
Genomic comparison of with 7 other species†
C.jeddahense
C.sporosphaeroides
C.cellulolyticum
C.dakarense
C.difficile
C.leptum
C.senegalense
C.beijerincki
C. jeddahense
3,462
1,573
876
816
847
1,030
770
1,044
C. sporosphaeroides
91.97
2,951
854
776
819
1,016
745
1,015
C. cellulolyticum
61.60
60.62
3,923
806
851
754
814
946
C. dakarense
57.30
56.34
65.70
5,020
1,271
665
1,110
1,142
C. difficile
57.56
56.80
65.65
77.74
3,390
714
1,098
1,171
C. leptum
67.98
68.07
61.94
58.55
58.84
3,591
651
780
C. senegalense
57.52
56.73
65.71
70.18
69.41
58.72
3,704
1,125
C. beijerincki
58.63
57.95
65.93
68.98
68.48
59.58
71.37
3,818
†Numbers of orthologous proteins shared between genomes (above diagonal), AGIOS values (below diagonal) and numbers of proteins per genome (bold numbers)
†Numbers of orthologous proteins shared between genomes (above diagonal), AGIOS values (below diagonal) and numbers of proteins per genome (bold numbers)
Genome properties
The genome is 3,613,503 bp long (1 chromosome, but no plasmid) with a 51.95% G+C content (Figure 6 and Table 4). Of the 3,515 predicted genes, 3,462 were protein-coding genes and 53 were RNAs, including 4 rRNAs. A total of 2,193 genes (62.38%) were assigned a putative function and 81 genes were identified as ORFans (2.3%). The properties and statistics of the genome are summarized in Tables 4 and 5. The distribution of genes into COG functional categories is presented in Table 5.
Figure 6
Graphical circular map of the chromosome. From the outside in: open reading frames oriented in the forward (colored by COG categories) direction, open reading frames oriented in the reverse (colored by COG categories) direction, RNA operon (red), and tRNAs (green), GC content plot, and GC skew (purple: negative values, olive: positive values).
Table 4
Nucleotide content and gene count levels of the genome
Attribute
Value
% of totala
Genome size (bp)
3,613,503
DNA G+C content (bp)
1,877,214
51.95
DNA Coding region (bp)
3,152,277
87.23
Number of replicons
1
Extra chromosomal element
0
Total genes
3,515
100
RNA genes
53
1.51
Protein-coding genes
3,462
98.49
Genes with function prediction
2,193
87.19
Genes assigned to COGs
2,515
71.55
Genes with peptide signals
135
3.84
Genes with transmembrane helices
887
25.23
a The total is based on either the size of the genome in base pairs or the total number of protein-coding genes in the annotated genome
Table 5
Number of genes associated with the 25 general COG functional categories
Code
Value
% agea
Description
J
154
4.45
Translation
A
0
0
RNA processing and modification
K
296
8.55
Transcription
L
138
3.98
Replication, recombination and repair
B
1
0.03
Chromatin structure and dynamics
D
24
0.69
Cell cycle control, mitosis and meiosis
Y
0
0
Nuclear structure
V
73
2.11
Defense mechanisms
T
156
4.5
Signal transduction mechanisms
M
116
3.35
Cell wall/membrane biogenesis
N
62
1.79
Cell motility
Z
0
0
Cytoskeleton
W
0
0
Extracellular structures
U
48
1.38
Intracellular trafficking and secretion
O
66
1.9
Posttranslational modification, protein turnover, chaperones
C
154
4.45
Energy production and conversion
G
237
6.84
Carbohydrate transport and metabolism
E
328
9.47
Amino acid transport and metabolism
F
56
1.61
Nucleotide transport and metabolism
H
92
2.66
Coenzyme transport and metabolism
I
85
2.45
Lipid transport and metabolism
P
164
4.74
Inorganic ion transport and metabolism
Q
53
1.53
Secondary metabolites biosynthesis, transport and catabolism
R
346
10
General function prediction only
S
195
5.63
Function unknown
-
947
27.35
Not in COGs
a The total is based on the total number of protein-coding genes in the annotated genome
Graphical circular map of the chromosome. From the outside in: open reading frames oriented in the forward (colored by COG categories) direction, open reading frames oriented in the reverse (colored by COG categories) direction, RNA operon (red), and tRNAs (green), GC content plot, and GC skew (purple: negative values, olive: positive values).a The total is based on either the size of the genome in base pairs or the total number of protein-coding genes in the annotated genomea The total is based on the total number of protein-coding genes in the annotated genome
Genome comparison with other genomes
We compared the genomes of JCDT, DSM 1294, DSM 753, C. beijerincki NCIMB 8052, H10, B1, DSM 25507, DSM 27086 (Table 6A).
Table 6A
Genomic comparison of with 7 other species†.
Species
Strain
Genome accessionnumber
Genome size (Mb)
G+C content
C. jeddahense
JCDT
CBYL00000000
3.61
51.95
C. sporosphaeroides
DSM 1294
ARTA01000000
3.17
53.5
C. cellulolyticum
H10
NC_011898
4.07
37.4
C. dakarense
DSM 27086
CBTZ010000000
3.73
27.98
C. difficile
B1
NC_017179
4.46
28.4
C. leptum
DSM 753
ABCB02000000
3.27
50.2
C. senegalense
DSM 25507
CAEV01000001
3.89
26.8
C. beijerincki
NCIMB 8052
NC_009617
6.0
29.0
† A: Species, Strain, GenBank accession number, genome size and G+C content of all compared genomes
† A: Species, Strain, GenBank accession number, genome size and G+C content of all compared genomesThe draft genome of (3.61 Mb) is larger than and (3.17 and 3.27 Mb respectively) but smaller than C. beijerincki, and (6.0, 4.07, 4.46, 3.89, 3.73 Mb respectively). It exhibits a higher G+C content than all other compared genome except (53.5%). has a higher gene content (3,462) than (2,951 and 3,390 respectively) but smaller than and (3,591, 3,818, 3,923, 3,704, 5,020 respectively). shared 1,573, 876, 816, 847, 1,030, 770 and 1,044 orthologous genes with and C. beijerincki respectively.When we compared with other species, AGIOS values ranged from 57.52 with to 91.97% with . Although the AGIOS value was elevated between and , we believe that the remarkable phenotypic differences, including motility, indole production (Table 2), and protein profile (Figure 7), enable the classification of as a new species.
Figure 7
Distribution of predicted genes of and 7 other species into COG categories. C.jdm= , C.spo= , C. lep= , C.bej = , C. cel = , C. diff = , C. sen = , C. dak = .
Distribution of predicted genes of and 7 other species into COG categories. C.jdm= , C.spo= , C. lep= , C.bej = , C. cel = , C. diff = , C. sen = , C. dak = .
Conclusion
On the basis of phenotypic, phylogenetic and genomic analyses (taxono-genomics), we formally propose the creation of sp. nov. that contains strain JCDT. This strain was isolated from the fecal flora of an obese 24 year-old Saudian individual living in Jeddah.
Description of sp. nov.
(jed.dah..en’.se L.gen. neutr. n. combination of Jeddah, the city in Saudi Arabia where the specimen was obtained from an obese Saudian patient sample.) Transparent colonies were 0.2 to 0.3 mm in diameter on blood-enriched agar. is a Gram-positive, obligate anaerobic, endospore-forming bacterium with a mean diameter of 1 µm. Optimal growth on axenic medium was observed at 37°C.is catalase negative and oxidase negative. Alkaline phosphatase, arginine arylamidase, proline arylamidase, alanine arylamidase, glycine arylamidase, histidine arylamidase, glutamyl glutamic acid arylamidase and serine arylamidase activities were positive. Arginine dihydrolase, α-galactosidase, β-galactosidase, α-glucosidase, β-glucosidase, α-arabinosidase, N-acetyl-β-glucosaminidase, glutamic acid decarboxylase, α-fucosidase, reduction of nitrate, leucyl glycine arylamidase, fermentation of mannose and raffinose, urease, β-galactosidase-6-phosphatase, β-glucuronidase, phenylalanine arylamidase, leucine arylamidase, pyroglutamic acid arylamidase and tyrosine arylamidase activities were negative. Asaccharolytic. Positive for indole. Cells are susceptible to amoxicillin, amoxicillin-clavulanate, imipenem, metronidazole, doxycycline, rifampicin, vancomycin but resistant to ceftriaxone, ciprofloxacin and trimethoprim-sulfamethoxazole.The G+C content of the genome is 51.95%. The 16S rDNA and genome sequences are deposited in GenBank under accession numbers HG726040 and CBYL00000000, respectively. The type strain is JCDT (= CSUR P693 = DSM 27834).
Authors: You Zhou; Yongjie Liang; Karlene H Lynch; Jonathan J Dennis; David S Wishart Journal: Nucleic Acids Res Date: 2011-06-14 Impact factor: 16.971
Authors: S Bakour; J Rathored; C I Lo; O Mediannikov; M Beye; C B Ehounoud; P Biagini; D Raoult; P-E Fournier; F Fenollar Journal: New Microbes New Infect Date: 2016-03-19
Authors: G A Durand; J-C Lagier; S Khelaifia; N Armstrong; C Robert; J Rathored; P-E Fournier; D Raoult Journal: New Microbes New Infect Date: 2016-02-10
Authors: E Guilhot; M T Alou; J-C Lagier; N Labas; C Couderc; J Delerce; A Diallo; C Sokhna; P-E Fournier; D Raoult; S Khelaifia Journal: New Microbes New Infect Date: 2017-02-02
Authors: M D Mbogning Fonkou; M Mailhe; S Ndongo; D Ricaboni; A Morand; F Cornu; M Tidjani Alou; M Bilen; C Andrieu; A Levasseur; F Cadoret; D Raoult Journal: New Microbes New Infect Date: 2018-06-18
Authors: Muhammad Yasir; Abeer M Ajlan; Shazi Shakil; Asif A Jiman-Fatani; Saad B Almasaudi; Muhammad Farman; Zainah M Baazeem; Rnda Baabdullah; Maha Alawi; Nabeela Al-Abdullah; Nashat A Ismaeel; Hani A Shukri; Esam I Azhar Journal: Sci Rep Date: 2018-10-04 Impact factor: 4.379
Authors: S Ndongo; M Beye; G Dubourg; T T Nguyen; C Couderc; D P Fabrizio; P-E Fournier; D Raoult; E Angelakis Journal: New Microbes New Infect Date: 2018-06-21