Menu

Auto Tag: ENCODE

September 22, 2019

Whole-genome sequencing and comparative analysis of two plant-associated strains of Rhodopseudomonas palustris (PS3 and YSC3).

Rhodopseudomonas palustris strains PS3 and YSC3 are purple non-sulfur phototrophic bacteria isolated from Taiwanese paddy soils. PS3 has beneficial effects on plant growth and enhances the uptake efficiency of applied fertilizer nutrients. In contrast, YSC3 has no significant effect on plant growth. The genomic structures of PS3 and YSC3 are similar; each contains one circular chromosome that is 5,269,926 or 5,371,816?bp in size, with 4,799 or 4,907 protein-coding genes, respectively. In this study, a large class of genes involved in chemotaxis and motility was identified in both strains, and genes associated with plant growth promotion, such as nitrogen fixation-, IAA synthesis- and ACC deamination-associated genes, were also identified. We noticed that the growth rate, the amount of biofilm formation, and the relative expression levels of several chemotaxis-associated genes were significantly higher for PS3 than for YSC3 upon treatment with root exudates. These results indicate that PS3 responds better to the presence of plant hosts, which may contribute to the successful interactions of PS3 with plant hosts. Moreover, these findings indicate that the existence of gene clusters associated with plant growth promotion is required but not sufficient for a bacterium to exhibit phenotypes associated with plant growth promotion.

Read more
September 22, 2019

Complete sequence of kenaf (Hibiscus cannabinus) mitochondrial genome and comparative analysis with the mitochondrial genomes of other plants.

Plant mitochondrial (mt) genomes are species specific due to the vast of foreign DNA migration and frequent recombination of repeated sequences. Sequencing of the mt genome of kenaf (Hibiscus cannabinus) is essential for elucidating its evolutionary characteristics. In the present study, single-molecule real-time sequencing technology (SMRT) was used to sequence the complete mt genome of kenaf. Results showed that the complete kenaf mt genome was 569,915?bp long and consisted of 62 genes, including 36 protein-coding, 3 rRNA and 23 tRNA genes. Twenty-five introns were found among nine of the 36 protein-coding genes, and five introns were trans-spliced. A comparative analysis with other plant mt genomes showed that four syntenic gene clusters were conserved in all plant mtDNAs. Fifteen chloroplast-derived fragments were strongly associated with mt genes, including the intact sequences of the chloroplast genes psaA, ndhB and rps7. According to the plant mt genome evolution analysis, some ribosomal protein genes and succinate dehydrogenase genes were frequently lost during the evolution of angiosperms. Our data suggest that the kenaf mt genome retained evolutionarily conserved characteristics. Overall, the complete sequencing of the kenaf mt genome provides additional information and enhances our better understanding of mt genomic evolution across angiosperms.

Read more
September 22, 2019

Whole-genome resequencing and pan-transcriptome reconstruction highlight the impact of genomic structural Variation on secondary metabolite gene clusters in the grapevine Esca pathogen Phaeoacremonium minimum.

The Ascomycete fungus Phaeoacremonium minimum is one of the primary causal agents of Esca, a widespread and damaging grapevine trunk disease. Variation in virulence among Pm. minimum isolates has been reported, but the underlying genetic basis of the phenotypic variability remains unknown. The goal of this study was to characterize intraspecific genetic diversity and explore its potential impact on virulence functions associated with secondary metabolism, cellular transport, and cell wall decomposition. We generated a chromosome-scale genome assembly, using single molecule real-time sequencing, and resequenced the genomes and transcriptomes of multiple isolates to identify sequence and structural polymorphisms. Numerous insertion and deletion events were found for a total of about 1 Mbp in each isolate. Structural variation in this extremely gene dense genome frequently caused presence/absence polymorphisms of multiple adjacent genes, mostly belonging to biosynthetic clusters associated with secondary metabolism. Because of the observed intraspecific diversity in gene content due to structural variation we concluded that a transcriptome reference developed from a single isolate is insufficient to represent the virulence factor repertoire of the species. We therefore compiled a pan-transcriptome reference of Pm. minimum comprising a non-redundant set of 15,245 protein-coding sequences. Using naturally infected field samples expressing Esca symptoms, we demonstrated that mapping of meta-transcriptomics data on a multi-species reference that included the Pm. minimum pan-transcriptome allows the profiling of an expanded set of virulence factors, including variable genes associated with secondary metabolism and cellular transport.

Read more
September 22, 2019

Characterization of a novel blaKLUC variant with reduced ß-lactam resistance from an IncA/C group plasmid in a clinical Klebsiella pneumoniae isolate.

Similar to other CTX-M family enzymes, KLUC is a recently identified and emerging determinant of cefotaxime resistance that has been recovered from at least three Enterobacteriaceae species, including Kluyvera cryocrescens, Escherichia coli, and Enterobacter cloacae. Whether this extended-spectrum ß-lactamase (ESBL) has been disseminated among commonly isolated Enterobacteriaceae is worthy of further investigation. In this study, we screened 739 nosocomial Enterobacteriaceae isolates (240 Klebsiella pneumoniae and 499 E. coli strains) and found that one K. pneumoniae and four E. coli isolates harbored the blaKLUC gene. Three blaKLUC determinants isolated from E. coli were entirely identical to a blaKLUC-3 gene previously recovered in the same hospital. PFGE of four blaKLUC-harboring E. coli strains showed that prevalence of these determinants was most likely mediated by horizontal gene transfer but not clonal dissemination. However, the variant isolated from K. pneumoniae belonged to a novel member of the KLUC enzyme group. This newly identified enzyme (KLUC-5) has an amino acid substitution compared with previously identified KLUC-1 (G18S) and KLUC-3 (G240D). Antimicrobial susceptibility tests showed that KLUC-5 significantly reduced resistance activity to almost all the selected antimicrobials compared to previously identified KLUC-3. Site-directed mutagenesis showed that blaKLUC-5-D240G and blaKLUC-5-S18G significantly enhanced the MIC against its best substrate. Conjugation and S1-PFGE indicated that blaKLUC-5 was located on a transferable plasmid, which was further decoded by single-molecule, real-time sequencing. Comparative genome analysis showed that its backbone exhibited genetic homology to the IncA/C incompatibility group plasmids. A transposable element, ISEcp1, was detected 256-bp upstream of the blaKLUC-5 gene; this location was inconsistent with the previously identified blaKLUC-1 but congruent with the variants recovered from E. coli in the same hospital. These data provide evidence of the increasingly emerging KLUC group of ESBLs in China.

Read more
September 22, 2019

Phenotypic and genomic properties of Brachybacterium vulturis sp. nov. and Brachybacterium avium sp. nov.

Two strains, VM2412T and VR2415T, were isolated from the feces of an Andean condor (Vultur gryphus) living in Seoul Grand Park, Gyeonggi-do, South Korea. Cells of both strains were observed to be Gram-stain positive, non-motile, aerobic, catalase positive and oxidase negative. Growth was found to occur at 10-30°C, showing optimum growth at 30°C. The strains could tolerate up to 15% (w/v) NaCl concentration and grow at pH 6-9. The strains shared 99.3% 16S rRNA gene sequence similarity to each other but were identified as two distinct species based on 89.0-89.2% ANIb, 90.3% ANIm, 89.7% OrthoANI and 38.0% dDDH values calculated using whole genome sequences. Among species with validly published names, Brachybacterium ginsengisoli DCY80T shared high 16S rRNA gene sequence similarities with strains VM2412T (98.7%) and VR2415T (98.4%) and close genetic relatedness with strains VM2412T (83.3-83.5% ANIb, 87.0% ANIm, 84.3% OrthoANI and 27.8% dDDH) and VR2415T (82.8-83.2% ANIb, 86.7% ANIm, 83.9% OrthoANI and 27.2% dDDH). The major fatty acid of the two strains was identified as anteiso-C15:0 and the polar lipids consisted of phosphatidylglycerol, diphosphatidylglycerol, presumptively phosphatidylethanolamine and three unidentified glycolipids. Strain VR2415T also produced an unidentified phospholipid. The cell walls of the two strains contained meso-diaminopimelic acid as diagnostic diamino acid and the whole cell sugars were ribose, glucose, and galactose. The strains contained MK-7 as their predominant menaquinone. The genomes of strains VM2412T, VR2415T, and B. ginsengisoli DCY80T were sequenced in this study. The genomic G+C contents of strains VM2412T and VR2415T were determined to be 70.8 and 70.4 mol%, respectively. A genome-based phylogenetic tree constructed using an up-to-date bacterial core gene set (UBCG) showed that the strains formed a clade with members of the genus Brachybacterium, supporting their taxonomic classification into the genus Brachybacterium. Based on phenotypic and genotypic analyses in this study, strains VM2412T and VR2415T are considered to represent two novel species of the genus Brachybacterium and the names Brachybacterium vulturis sp. nov. and Brachybacterium avium sp. nov. are proposed for strains VM2412T (=KCTC 39996T = JCM 32142T) and VR2415T (=KCTC 39997T = JCM 32143T), respectively.

Read more
September 22, 2019

Large plasmidome of dairy Lactococcus lactis subsp. lactis biovar diacetylactis FM03P encodes technological functions and appears highly unstable.

Important industrial traits have been linked to plasmids in Lactococcus lactis.The dairy isolate L. lactis subsp. lactis biovar diacetylactis FM03P was sequenced revealing the biggest plasmidome of all completely sequenced and published L. lactis strains up till now. The 12 plasmids that were identified are: pLd1 (8277 bp), pLd2 (15,218 bp), pLd3 (4242 bp), pLd4 (12,005 bp), pLd5 (7521 bp), pLd6 (3363 bp), pLd7 (30,274 bp), pLd8 (47,015 bp), pLd9 (15,313 bp), pLd10 (39,563 bp), pLd11 (9833 bp) and pLd12 (3321 bp). Structural analysis of the repB promoters and the RepB proteins showed that eleven of the plasmids replicate via the theta-type mechanism, while only plasmid pLd3 replicates via a rolling-circle replication mechanism. Plasmids pLd2, pLd7 and pLd10 contain a highly similar operon involved in mobilisation of the plasmids. Examination of the twelve plasmids of L. lactis FM03P showed that 10 of the plasmids carry putative genes known to be important for growth and survival in the dairy environment. These genes encode technological functions such as lactose utilisation (lacR-lacABCDFEGX), citrate uptake (citQRP), peptide degradation (pepO and pepE) and oligopeptide uptake (oppDFBCA), uptake of magnesium and manganese (2 mntH, corA), exopolysaccharides production (eps operon), bacteriophage resistance (1 hsdM, 1 hsdR and 7 different hsdS genes of a type I restriction-modification system, an operon of three genes encoding a putative type II restriction-modification system and an abortive infection gene) and stress resistance (2 uspA, cspC and cadCA). Acquisition of these plasmids most likely facilitated the adaptation of the recipient strain to the dairy environment. Some plasmids were already lost during a single propagation step signifying their instability in the absence of a selective pressure.Lactococcus lactis FM03P carries 12 plasmids important for its adaptation to the dairy environment. Some of the plasmids were easily lost demonstrating that propagation outside the dairy environment should be minimised when studying dairy isolates of L. lactis.

Read more
September 22, 2019

Complete genome of streamlined marine actinobacterium Pontimonas salivibrio strain CL-TW6T adapted to coastal planktonic lifestyle.

Pontimonas salivibrio strain CL-TW6T (=KCCM 90105?=?JCM18206) was characterized as the type strain of a new genus within the Actinobacterial family Microbacteriaceae. It was isolated from a coastal marine environment in which members of Microbactericeae have not been previously characterized.The genome of P. salivibrio CL-TW6T was a single chromosome of 1,760,810 bp. Genomes of this small size are typically found in bacteria growing slowly in oligotrophic zones and said to be streamlined. Phylogenetic analysis showed it to represent a lineage originating in the Microbacteriaceae radiation occurring before the snowball Earth glaciations, and to have a closer relationship with some streamlined bacteria known through metagenomic data. Several genomic characteristics typical of streamlined bacteria are found: %G?+?C is lower than non-streamlined members of the phylum; there are a minimal number of rRNA and tRNA genes, fewer paralogs in most gene families, and only two sigma factors; there is a noticeable absence of some nonessential metabolic pathways, including polyketide synthesis and catabolism of some amino acids. There was no indication of any phage genes or plasmids, however, a system of active insertion elements was present. P. salivibrio appears to be unusual in having polyrhamnose-based cell wall oligosaccharides instead of mycolic acid or teichoic acid-based oligosaccharides. Oddly, it conducts sulfate assimilation apparently for sulfating cell wall components, but not for synthesizing amino acids. One gene family it has more of, rather than fewer of, are toxin/antitoxin systems, which are thought to down-regulate growth during nutrient deprivation or other stressful conditions.Because of the relatively small number of paralogs and its relationship to the heavily characterized Mycobacterium tuberculosis, we were able to heavily annotate the genome of P. salivibrio CL-TW6T. Its streamlined status and relationship to streamlined metagenomic constructs makes it an important reference genome for study of the streamlining concept. The final evolutionary trajectory of CL-TW6 T was to adapt to growth in a non-oligotrophic coastal zone. To understand that adaptive process, we give a thorough accounting of gene content, contrasting with both oligotrophic streamlined bacteria and large genome bacteria, and distinguishing between genes derived by vertical and horizontal descent.

Read more
September 22, 2019

Reprogramming of the antimycin NRPS-PKS assembly lines inspired by gene evolution.

Reprogramming of the NRPS/PKS assembly line is an attractive method for the production of new bioactive molecules. However, it is usually hampered by the loss of intimate domain/module interactions required for the precise control of chain transfer and elongation reactions. In this study, we first establish heterologous expression systems of the unique antimycin-type cyclic depsipeptides: JBIR-06 (tri-lactone) and neoantimycin (tetra-lactone), and engineer their biosyntheses by taking advantage of bioinformatic analyses and evolutionary insights. As a result, we successfully accomplish three manipulations: (i) ring contraction of neoantimycin (from tetra-lactone to tri-lactone), (ii) ring expansion of JBIR-06 (from tri-lactone to tetra-lactone), and (iii) alkyl chain diversification of JBIR-06 by the incorporation of various alkylmalonyl-CoA extender units, to generate a set of unnatural derivatives in practical yields. This study presents a useful strategy for engineering NRPS-PKS module enzymes, based on nature’s diversification of the domain and module organizations.

Read more
September 22, 2019

Discovery of new genes involved in curli production by a uropathogenic Escherichia coli strain from the highly virulent O45:K1:H7 lineage.

Curli are bacterial surface-associated amyloid fibers that bind to the dye Congo red (CR) and facilitate uropathogenic Escherichia coli (UPEC) biofilm formation and protection against host innate defenses. Here we sequenced the genome of the curli-producing UPEC pyelonephritis strain MS7163 and showed it belongs to the highly virulent O45:K1:H7 neonatal meningitis-associated clone. MS7163 produced curli at human physiological temperature, and this correlated with biofilm growth, resistance of sessile cells to the human cationic peptide cathelicidin, and enhanced colonization of the mouse bladder. We devised a forward genetic screen using CR staining as a proxy for curli production and identified 41 genes that were required for optimal CR binding, of which 19 genes were essential for curli synthesis. Ten of these genes were novel or poorly characterized with respect to curli synthesis and included genes involved in purine de novo biosynthesis, a regulator that controls the Rcs phosphorelay system, and a novel repressor of curli production (referred to as rcpA). The involvement of these genes in curli production was confirmed by the construction of defined mutants and their complementation. The mutants did not express the curli major subunit CsgA and failed to produce curli based on CR binding. Mutation of purF (the first gene in the purine biosynthesis pathway) and rcpA also led to attenuated colonization of the mouse bladder. Overall, this work has provided new insight into the regulation of curli and the role of these amyloid fibers in UPEC biofilm formation and pathogenesis.IMPORTANCE Uropathogenic Escherichia coli (UPEC) strains are the most common cause of urinary tract infection, a disease increasingly associated with escalating antibiotic resistance. UPEC strains possess multiple surface-associated factors that enable their colonization of the urinary tract, including fimbriae, curli, and autotransporters. Curli are extracellular amyloid fibers that enhance UPEC virulence and promote biofilm formation. Here we examined the function and regulation of curli in a UPEC pyelonephritis strain belonging to the highly virulent O45:K1:H7 neonatal meningitis-associated clone. Curli expression at human physiological temperature led to increased biofilm formation, resistance of sessile cells to the human cationic peptide LL-37, and enhanced bladder colonization. Using a comprehensive genetic screen, we identified multiple genes involved in curli production, including several that were novel or poorly characterized with respect to curli synthesis. In total, this study demonstrates an important role for curli as a UPEC virulence factor that promotes biofilm formation, resistance, and pathogenesis. Copyright © 2018 Nhu et al.

Read more
September 22, 2019

Evolution of the U.S. biological select agent Rathayibacter toxicus.

Rathayibacter toxicus is a species of Gram-positive, corynetoxin-producing bacteria that causes annual ryegrass toxicity, a disease often fatal to grazing animals. A phylogenomic approach was employed to model the evolution of R. toxicus to explain the low genetic diversity observed among isolates collected during a 30-year period of sampling in three regions of Australia, gain insight into the taxonomy of Rathayibacter, and provide a framework for studying these bacteria. Analyses of a data set of more than 100 sequenced Rathayibacter genomes indicated that Rathayibacter forms nine species-level groups. R. toxicus is the most genetically distant, and evidence suggested that this species experienced a dramatic event in its evolution. Its genome is significantly reduced in size but is colinear to those of sister species. Moreover, R. toxicus has low intergroup genomic diversity and almost no intragroup genomic diversity between ecologically separated isolates. R. toxicus is the only species of the genus that encodes a clustered regularly interspaced short palindromic repeat (CRISPR) locus and that is known to host a bacteriophage parasite. The spacers, which represent a chronological history of infections, were characterized for information on past events. We propose a three-stage process that emphasizes the importance of the bacteriophage and CRISPR in the genome reduction and low genetic diversity of the R. toxicus species.IMPORTANCERathayibacter toxicus is a toxin-producing species found in Australia and is often fatal to grazing animals. The threat of introduction of the species into the United States led to its inclusion in the Federal Select Agent Program, which makes R. toxicus a highly regulated species. This work provides novel insights into the evolution of R. toxicusR. toxicus is the only species in the genus to have acquired a CRISPR adaptive immune system to protect against bacteriophages. Results suggest that coexistence with the bacteriophage NCPPB3778 led to the massive shrinkage of the R. toxicus genome, species divergence, and the maintenance of low genetic diversity in extant bacterial groups. This work contributes to an understanding of the evolution and ecology of an agriculturally important species of bacteria. Copyright © 2018 Davis et al.

Read more
September 22, 2019

Integration of genomic data with NMR analysis enables assignment of the full stereostructure of neaumycin B, a potent inhibitor of glioblastoma from a marine-derived Micromonospora.

The microbial metabolites known as the macrolides are some of the most successful natural products used to treat infectious and immune diseases. Describing the structures of these complex metabolites, however, is often extremely difficult due to the presence of multiple stereogenic centers inherent in this class of polyketide-derived metabolites. With the availability of genome sequence data and a better understanding of the molecular genetics of natural product biosynthesis, it is now possible to use bioinformatic approaches in tandem with spectroscopic tools to assign the full stereostructures of these complex metabolites. In our quest to discover and develop new agents for the treatment of cancer, we observed the production of a highly cytotoxic macrolide, neaumycin B, by a marine-derived actinomycete bacterium of the genus Micromonospora. Neaumycin B is a complex polycyclic macrolide possessing 19 asymmetric centers, usually requiring selective degradation, crystallization, derivatization, X-ray diffraction analysis, synthesis, or other time-consuming approaches to assign the complete stereostructure. As an alternative approach, we sequenced the genome of the producing strain and identified the neaumycin gene cluster ( neu). By integrating the known stereospecificities of biosynthetic enzymes with comprehensive NMR analysis, the full stereostructure of neaumycin B was confidently assigned. This approach exemplifies how mining gene cluster information while integrating NMR-based structure data can achieve rapid, efficient, and accurate stereostructural assignments for complex macrolides.

Read more
September 22, 2019

The hpRNA/RNAi pathway is essential to resolve intragenomic conflict in the Drosophila male germline.

Intragenomic conflicts are fueled by rapidly evolving selfish genetic elements, which induce selective pressures to innovate opposing repressive mechanisms. This is patently manifest in sex-ratio (SR) meiotic drive systems, in which distorter and suppressor factors bias and restore equal transmission of X and Y sperm. Here, we reveal that multiple SR suppressors in Drosophila simulans (Nmy and Tmy) encode related hairpin RNAs (hpRNAs), which generate endo-siRNAs that repress the paralogous distorters Dox and MDox. All components in this drive network are recently evolved and largely testis restricted. To connect SR hpRNA function to the RNAi pathway, we generated D. simulans null mutants of Dcr-2 and AGO2. Strikingly, these core RNAi knockouts massively derepress Dox and MDox and are in fact completely male sterile and exhibit highly defective spermatogenesis. Altogether, our data reveal how the adaptive capacity of hpRNAs is critically deployed to restrict selfish gonadal genetic systems that can exterminate a species. Copyright © 2018 Elsevier Inc. All rights reserved.

Read more
September 22, 2019

Opposite polarity monospore genome de novo sequencing and comparative analysis reveal the possible heterothallic life cycle of Morchella importuna.

Morchella is a popular edible fungus worldwide due to its rich nutrition and unique flavor. Many research efforts were made on the domestication and cultivation of Morchella all over the world. In recent years, the cultivation of Morchella was successfully commercialized in China. However, the biology is not well understood, which restricts the further development of the morel fungus cultivation industry. In this paper, we performed de novo sequencing and assembly of the genomes of two monospores with a different mating type (M04M24 and M04M26) isolated from the commercially cultivated strain M04. Gene annotation and comparative genome analysis were performed to study differences in CAZyme (Carbohydrate-active enzyme) enzyme content, transcription factors, duplicated sequences, structure of mating type sites, and differences at the gene and functional levels between the two monospore strains of M. importuna. Results showed that the de novo assembled haploid M04M24 and M04M26 genomes were 48.98 and 51.07 Mb, respectively. A complete fine physical map of M. importuna was obtained from genome coverage and gene completeness evaluation. A total of 10,852 and 10,902 common genes and 667 and 868 endemic genes were identified from the two monospore strains, respectively. The Gene Ontology (GO) and KAAS (KEGG Automatic Annotation Serve) enrichment analyses showed that the endemic genes performed different functions. The two monospore strains had 99.22% collinearity with each other, accompanied with certain position and rearrangement events. Analysis of complete mating-type loci revealed that the two monospore M. importuna strains contained an independent mating-type structure and remained conserved in sequence and location. The phylogenetic and divergence time of M. importuna was analyzed at the whole-genome level for the first time. The bifurcation time of morel and tuber was estimated to be 201.14 million years ago (Mya); the two monospore strains with a different mating type represented the evolution of different nuclei, and the single copy homologous genes between them were also different due to a genetic differentiation distance about 0.65 Mya. Compared with truffles, M. importuna had an extension of 28 clusters of orthologous genes (COGs) and a contraction of two COGs. The two different polar nuclei with different degrees of contraction and expansion suggested that they might have undergone different evolutionary processes. The different mating-type structures, together with the functional clustering and enrichment analysis results of the endemic genes of the two different polar nuclei, imply that M. importuna might be a heterothallic fungus and the interaction between the endemic genes may be necessary for its complete life history. Studies on the genome of M. importuna facilitate a better understanding of morel biology and evolution.

Read more
September 22, 2019

Comparative genomics of Salmonella enterica serovar Montevideo reveals lineage-specific gene differences that may influence ecological niche association.

Salmonella enterica serovar Montevideo has been linked to recent foodborne illness outbreaks resulting from contamination of products such as fruits, vegetables, seeds and spices. Studies have shown that Montevideo also is frequently associated with healthy cattle and can be isolated from ground beef, yet human salmonellosis outbreaks of Montevideo associated with ground beef contamination are rare. This disparity fuelled our interest in characterizing the genomic differences between Montevideo strains isolated from healthy cattle and beef products, and those isolated from human patients and outbreak sources. To that end, we sequenced 13 Montevideo strains to completion, producing high-quality genome assemblies of isolates from human patients (n=8) or from healthy cattle at slaughter (n=5). Comparative analysis of sequence data from this study and publicly available sequences (n=72) shows that Montevideo falls into four previously established clades, differentially occupied by cattle and human strains. The results of these analyses reveal differences in metabolic islands, environmental adhesion determinants and virulence factors within each clade, and suggest explanations for the infrequent association between bovine isolates and human illnesses.

Read more
September 22, 2019

Distinct genomic features characterize two clades of Corynebacterium diphtheriae: Proposal of Corynebacterium diphtheriae subsp. diphtheriae subsp. nov. and Corynebacterium diphtheriae subsp. lausannense subsp. nov.

Corynebacterium diphtheriae is the etiological agent of diphtheria, a disease caused by the presence of the diphtheria toxin. However, an increasing number of records report non-toxigenic C. diphtheriae infections. Here, a C. diphtheriae strain was recovered from a patient with a past history of bronchiectasis who developed a severe tracheo-bronchitis with multiple whitish lesions of the distal trachea and the mainstem bronchi. Whole-genome sequencing (WGS), performed in parallel with PCR targeting the toxin gene and the Elek test, provided clinically relevant results in a short turnaround time, showing that the isolate was non-toxigenic. A comparative genomic analysis of the new strain (CHUV2995) with 56 other publicly available genomes of C. diphtheriae revealed that the strains CHUV2995, CCUG 5865 and CMCNS703 share a lower average nucleotide identity (ANI) (95.24 to 95.39%) with the C. diphtheriae NCTC 11397T reference genome than all other C. diphtheriae genomes (>98.15%). Core genome phylogeny confirmed the presence of two monophyletic clades. Based on these findings, we propose here two new C. diphtheriae subspecies to replace the lineage denomination used in previous multilocus sequence typing studies: C. diphtheriae subsp. lausannense subsp. nov. (instead of lineage-2), regrouping strains CHUV2995, CCUG 5865, and CMCNS703, and C. diphtheriae subsp. diphtheriae subsp. nov, regrouping all other C. diphtheriae in the dataset (instead of lineage-1). Interestingly, members of subspecies lausannense displayed a larger genome size than subspecies diphtheriae and were enriched in COG categories related to transport and metabolism of lipids (I) and inorganic ion (P). Conversely, they lacked all genes involved in the synthesis of pili (SpaA-type, SpaD-type and SpaH-type), molybdenum cofactor and of the nitrate reductase. Finally, the CHUV2995 genome is particularly enriched in mobility genes and harbors several prophages. The genome encodes a type II-C CRISPR-Cas locus with 2 spacers that lacks csn2 or cas4, which could hamper the acquisition of new spacers and render strain CHUV2995 more susceptible to bacteriophage infections and gene acquisition through various mechanisms of horizontal gene transfer.

Read more

Subscribe for blog updates:

Talk with an expert

If you have a question, need to check the status of an order, or are interested in purchasing an instrument, we're here to help.