The PATRIC Bioinformatics Resource Center: expanding data and analysis capabilities

Natural product biosynthetic potential reflects macroevolutionary diversification within a widely distributed bacterial taxon

IMPORTANCE

This is the most comprehensive study performed thus far on the biosynthetic potential within the Flavobacteriaceae family. Our findings reveal intertwined taxonomic and natural product biosynthesis diversification within the family. We posit that the carbohydrate, peptide, and secondary metabolism triad synergistically shaped the evolution of this keystone bacterial taxon, acting as major forces underpinning the broad host range and opportunistic-to-pathogenic behavior encompassed by species in the family. This study further breaks new ground for future research on select Flavobacteriaceae spp. as reservoirs of novel drug leads.

Complete genome sequences of Lactococcus lactis D1_2, a bacterium with antimicrobial properties isolated from peat soil

Lactococcus lactis is a beneficial lactic acid bacterium commonly studied for its probiotic properties and role in dairy production. Here, we present a complete genome of Lactococcus lactis D1_2, isolated from peat swamp forests. To discover the potential antimicrobial properties, the complete genome of the strain was sequenced and analyzed.

Decoding stoichiometric protein synthesis in E. coli through translation rate parameters

E. coli is one of the most widely used organisms for understanding the principles of cellular and molecular genetics. However, we are yet to understand the origin of several experimental observations related to the regulation of gene expression in E. coli. One of the prominent examples in this context is the proportional synthesis in multiprotein complexes where all of their obligate subunits are produced in proportion to their stoichiometry. In this work, by combining the next-generation sequencing data with the stochastic simulations of protein synthesis, we explain the origin of proportional protein synthesis in multicomponent complexes. We find that the estimated initiation rates for the translation of all subunits in those complexes are proportional to their stoichiometry. This constraint on protein synthesis kinetics enforces proportional protein synthesis without requiring any feedback mechanism. We also find that the translation initiation rates in E. coli are influenced by the coding sequence length and the enrichment of A and C nucleotides near the start codon. Thus, this study rationalizes the role of conserved and nonrandom features of genes in regulating the translation kinetics and unravels a key principle of the regulation of protein synthesis.

Long-term exposure to food-grade disinfectants causes cross-resistance to antibiotics in Salmonella enterica serovar Typhimurium strains with different antibiograms and sequence types

BACKGROUND

Disinfectants are important in the food industry to prevent the transmission of pathogens. Excessive use of disinfectants may increase the probability of bacteria experiencing long-term exposure and consequently resistance and cross-resistance to antibiotics. This study aims to investigate the cross-resistance of multidrug-resistant, drug-resistant, and drug-susceptible isolates of Salmonella enterica serovar Typhimurium (S. Typhimurium) with different sequence types (STs) to a group of antibiotics after exposure to different food-grade disinfectants.

METHODS

A panel of 27 S. Typhimurium strains with different antibiograms and STs were exposed to increasing concentrations of five food-grade disinfectants, including hydrogen peroxide (H2O2), benzalkonium chloride (BAC), chlorine dioxide (ClO2), sodium hypochlorite (NaClO), and ethanol. Recovered evolved strains were analyzed using genomic tools and phenotypic tests. Genetic mutations were screened using breseq pipeline and changes in resistance to antibiotics and to the same disinfectant were determined. The relative fitness of evolved strains was also determined.

RESULTS

Following exposure to disinfectants, 22 out of 135 evolved strains increased their resistance to antibiotics from a group of 14 clinically important antibiotics. The results also showed that 9 out of 135 evolved strains had decreased resistance to some antibiotics. Genetic mutations were found in evolved strains. A total of 77.78% of ST34, 58.33% of ST19, and 66.67% of the other STs strains exhibited changes in antibiotic resistance. BAC was the disinfectant that induced the highest number of strains to cross-resistance to antibiotics. Besides, H2O2 induced the highest number of strains with decreased resistance to antibiotics.

CONCLUSIONS

These findings provide a basis for understanding the effect of disinfectants on the antibiotic resistance of S. Typhimurium. This work highlights the link between long-term exposure to disinfectants and the evolution of resistance to antibiotics and provides evidence to promote the regulated use of disinfectants.

A metagenomics method for the quantitative detection of bacterial pathogens causing hospital-associated and ventilator-associated pneumonia

IMPORTANCE

The management of ventilator-associated pneumonia and hospital-acquired pneumonia requires rapid and accurate quantitative detection of the infecting pathogen. To this end, we propose a metagenomic sequencing assay that includes the use of an internal sample processing control for the quantitative detection of 20 relevant bacterial species from bronchoalveolar lavage samples.

Quantitative proteomics defines mechanisms of antiviral defence and cell death during modified vaccinia Ankara infection

Modified vaccinia Ankara (MVA) virus does not replicate in human cells and is the vaccine deployed to curb the current outbreak of mpox. Here, we conduct a multiplexed proteomic analysis to quantify >9000 cellular and ~80% of viral proteins throughout MVA infection of human fibroblasts and macrophages. >690 human proteins are down-regulated >2-fold by MVA, revealing a substantial remodelling of the host proteome. >25% of these MVA targets are not shared with replication-competent vaccinia. Viral intermediate/late gene expression is necessary for MVA antagonism of innate immunity, and suppression of interferon effectors such as ISG20 potentiates virus gene expression. Proteomic changes specific to infection of macrophages indicate modulation of the inflammatory response, including inflammasome activation. Our approach thus provides a global view of the impact of MVA on the human proteome and identifies mechanisms that may underpin its abortive infection. These discoveries will prove vital to design future generations of vaccines.

Sinomonas cellulolyticus sp. nov., isolated from Loktak lake

A novel cellulolytic strain JC656T was isolated from the rhizosphere soil of Alisma plantago-aquatica of floating island (Phumdis) of Loktak lake, Manipur, India. The 16S rRNA gene sequence similarities between strain JC656T and other Sinomonas type strains ranged between 98.5 and 97.3%, wherein strain JC656T exhibited the highest sequence similarity (98.5%) to Sinomonas notoginsengisoli KCTC 29237T. Colonies were yellow-colored and grew aerobically. Cells were gram-positive, rod-shaped and non-motile. The optimal growth of the strain JC656T occured at 28 °C and pH 7. Strain JC656T contained MK-9 as the predominant isoprenoid quinone and anteiso-C15:0, iso-C16:0 and anteiso-C17:0 as the major fatty acids. Diphosphatidylglycerol, phosphatidylinositol, phosphatidylglycerol, phosphatidylmonomethylethanolamine and a glycolipid were the polar lipids. Strain JC656T contained lysine, alanine, glutamine, diaminopimelic acid (DAP) and two unidentified amino acids as characteristic cell wall amino acids. The genome size of strain JC656T was 3.9 Mb with a DNA G + C content of 69.9 mol %. For the affirmation of the strain's taxonomic status, a detailed phylogenomic study was done. Based on its phylogenetic position and morphological, physiological, and genomic features, strain JC656T represents a new species of the genus Sinomonas, for which we propose the name Sinomonas cellulolyticus sp. nov. The type strain JC656T = (KCTC 49339T = NBRC 114142T).

Draft genome sequencing of halotolerant bacterium Salinicola sp. DM10 unravels plant growth-promoting potentials

In this study, strain DM10 was isolated from mangrove roots and characterized as a halotolerant plant growth-promoting bacterium. Strain DM10 exhibited the ability to solubilize phosphate, produce siderophore, show 1-aminocyclopropane-1-carboxylic acid deaminase activity, and hydrolyze starch. The rice plants subjected to a treatment of NaCl (200 mM) and inoculated with strain DM10 showed an improvement in the shoot length, root length, and dried weight, when compared to those exposed solely to saline treatment. The comprehensive genome sequencing of strain DM10 revealed a genome spanning of 4,171,745 bp, harboring 3626 protein coding sequences. Within its genome, strain DM10 possesses genes responsible for both salt-in and salt-out strategies, indicative of a robust genetic adaptation aimed at fostering salt tolerance. Additionally, the genome encodes genes involved in phosphate solubilization, such as the synthesis of gluconic acid, high-affinity phosphate transport systems, and alkaline phosphatase. In the genome of DM10, we identified the acdS gene, responsible for encoding 1-aminocyclopropane-1-carboxylate deaminase, as well as the amy1A gene, which encodes α-amylase. Furthermore, the genome of DM10 contains sequences associated with the iron (3+)-hydroxamate and iron uptake clusters, responsible for siderophore production. Such data provide a deep understanding of the mechanism employed by strain DM10 to combat osmotic and salinity stress, facilitate plant growth, and elucidate its molecular-level behaviors.

Insights from metagenome-assembled genomes on the genetic stability and safety of over-the-counter probiotic products

The demand for and acceptance of probiotics is determined by their quality and safety. Illumina NGS sequencing and analytics were used to examine eight marketed probiotics. Up to the species level, sequenced DNA was taxonomically identified, and relative abundances were determined using Kaiju. The genomes were constructed using GTDB and validated through PATRICK and TYGS. A FastTree 2 phylogenetic tree was constructed using several type strain sequences from relevant species. Bacteriocin and ribosomally synthesized polypeptide (RiPP) genes were discovered, and a safety check was performed to test for toxins, antibiotic resistance, and genetic drift genes. Except for two products with unclaimed species, the labeling was taxonomically correct. In three product formulations, Lactobacillus acidophilus, Limosilactobacillus reuteri, Lacticaseibacillus paracasei, and Bifidobacterium animalis exhibited two to three genomic alterations, while Streptococcus equinus was found in one. TYGS and GDTB discovered E. faecium and L. paracasei in distinctly different ways. All the bacteria tested had the genetic repertoire to tolerate GIT transit, although some exhibited antibiotic resistance, and one strain had two virulence genes. Except for Bifidobacterium strains, the others revealed a variety of bacteriocins and ribosomally synthesized polypeptides (RiPP), 92% of which were unique and non-homologous to known ones. Plasmids and mobile genetic elements are present in strains of L. reuteri (NPLps01.et_L.r and NPLps02.uf_L.r), Lactobacillus delbrueckii (NPLps01.et_L.d), Streptococcus thermophilus (NPLps06.ab_S.t), and E. faecium (NPLps07.nf_E.f). Our findings support the use of metagenomics to build better and efficient production and post-production practices for probiotic quality and safety assessment.

Proposal to Reclassify Companilactobacillus futsaii subsp. chongqingensis as a Later Heterotypic Synonym of Companilactobacillus futsaii subsp. futsaii

Previous studies have shown that Lactobacillus futsaii (now Companilactobacillus futsaii) can be subdivided at the subspecies level. The main purpose of this study is to explore whether this is correct by using a polyphasic taxonomic approach. Lactobacillus futsaii subsp. chongqingii was proposed and effectively published in 2019. The names L. futsaii subsp. chongqingensis corrig. and Lactobacillus futsaii subsp. futsaii were not validated until March 2023. However, in the reclassification of the genus Lactobacillus by Zheng et al. in April 2020, L. futsaii was transferred to Companilactobacillus as Companilactobacillus futsaii. So Lactobacillus futsaii subsp. chongqingensis and Lactobacillus futsaii subsp. futsaii should be transferred to Companilactobacillus futsaii now. In the present study, the relationship between L. futsaii subsp. chongqingensis and L. futsaii subsp. futsaii was re-evaluated. The type strains of L. futsaii subsp. chongqingensis and L. futsaii subsp. futsaii shared identical pheS and rpoA sequences, high dDDH value, similar phenotypic characteristics and fatty acid compositions, indicating that they belonged to the same subspecies. Here, we propose to reclassify Lactobacillus futsaii subsp. chongqingensis and Lactobacillus futsaii subsp. futsaii as Companilactobacillus futsaii subsp. chongqingensis comb. nov. and Companilactobacillus futsaii subsp. futsaii comb. nov., respectively, and Companilactobacillus futsaii subsp. chongqingensis as a later heterotypic synonym of Companilactobacillus futsaii subsp. futsaii.

Global pathogenomic analysis identifies known and candidate genetic antimicrobial resistance determinants in twelve species

Surveillance programs for managing antimicrobial resistance (AMR) have yielded thousands of genomes suited for data-driven mechanism discovery. We present a workflow integrating pangenomics, gene annotation, and machine learning to identify AMR genes at scale. When applied to 12 species, 27,155 genomes, and 69 drugs, we 1) find AMR gene transfer mostly confined within related species, with 925 genes in multiple species but just eight in multiple phylogenetic classes, 2) demonstrate that discovery-oriented support vector machines outperform contemporary methods at recovering known AMR genes, recovering 263 genes compared to 145 by Pyseer, and 3) identify 142 AMR gene candidates. Validation of two candidates in E. coli BW25113 reveals cases of conditional resistance: ΔcycA confers ciprofloxacin resistance in minimal media with D-serine, and frdD V111D confers ampicillin resistance in the presence of ampC by modifying the overlapping promoter. We expect this approach to be adaptable to other species and phenotypes.

Genome sequence of the virulent Aeromonas salmonicida atypical strain T30 isolated from sablefish with furunculosis

Aeromonas salmonicida, a Gram-negative bacterium, causes the disease furunculosis in multiple fish species. We present the complete genome sequence of the atypical A. salmonicida strain T30, which was isolated from furunculosis in sablefish in Manchester, WA, USA. Analyzing this genome will help to identify the bacterium's role in marine aquaculture.

Metagenomics for the microbiological diagnosis of hospital-acquired pneumonia and ventilator-associated pneumonia (HAP/VAP) in intensive care unit (ICU): a proof-of-concept study

BACKGROUND

Hospital-acquired and ventilator-associated-pneumonia (HAP/VAP) are one of the most prevalent health-care associated infections in the intensive care unit (ICU). Culture-independent methods were therefore developed to provide faster route to diagnosis and treatment. Among these, metagenomic next-generation sequencing (mNGS) has shown considerable promise.

METHODS

This proof-of-concept study describes the technical feasibility and evaluates the clinical validity of the mNGS for the detection and characterization of the etiologic agents causing hospital-acquired and ventilator-associated pneumonia. We performed a prospective study of all patients with HAP/VAP hospitalized in our intensive care unit for whom a bronchoalveolar lavage (BAL) was performed between July 2017 and November 2018. We compared BAL fluid culture and mNGS results of these patients.

RESULTS

A total of 32 BAL fluids were fully analyzed. Of these, 22 (69%) were positive by culture and all pathogens identified were also reported by mNGS. Among the culture-positive BAL samples, additional bacterial species were revealed by mNGS for 12 patients, raising the issue of their pathogenic role (colonization versus coinfection). Among BALF with culture-negative test, 5 were positive in mNGS test.

CONCLUSIONS

This study revealed concordant results for pneumonia panel pathogens between mNGS and culture-positive tests and identified additional pathogens potentially implicated in pneumonia without etiologic diagnosis by culture. mNGS has emerged as a promising methodology for infectious disease diagnoses to support conventional methods. Prospective studies with real-time mNGS are warranted to examine the impact on antimicrobial decision-making and clinical outcome.

Population genomic analysis of clinical ST15 Klebsiella pneumoniae strains in China

ST15 Klebsiella pneumoniae (Kpn) is a growing public health concern in China and worldwide, yet its genomic and evolutionary dynamics in this region remain poorly understood. This study comprehensively elucidates the population genomics of ST15 Kpn in China by analyzing 287 publicly available genomes. The proportion of the genomes increased sharply from 2012 to 2021, and 92.3% of them were collected from the Yangtze River Delta (YRD) region of eastern China. Carbapenemase genes, including OXA-232, KPC-2, and NDM, were detected in 91.6% of the studied genomes, and 69.2% of which were multidrug resistant (MDR) and hypervirulent (hv). Phylogenetic analysis revealed four clades, C1 (KL112, 59.2%), C2 (mainly KL19, 30.7%), C3 (KL48, 0.7%) and C4 (KL24, 9.4%). C1 appeared in 2007 and was OXA-232-producing and hv; C2 and C4 appeared between 2005 and 2007, and both were KPC-2-producing but with different levels of virulence. Transmission clustering detected 86.1% (n = 247) of the enrolled strains were grouped into 55 clusters (2-159 strains) and C1 was more transmissible than others. Plasmid profiling revealed 88 plasmid clusters (PCs) that were highly heterogeneous both between and within clades. 60.2% (n = 53) of the PCs carrying AMR genes and 7 of which also harbored VFs. KPC-2, NDM and OXA-232 were distributed across 14, 4 and 1 PCs, respectively. The MDR-hv strains all carried one of two homologous PCs encoding iucABCD and rmpA2 genes. Pangenome analysis revealed two major coinciding accessory components predominantly located on plasmids. One component, associated with KPC-2, encompassed 15 additional AMR genes, while the other, linked to OXA-232, involved seven more AMR genes. This study provides essential insights into the genomic evolution of the high-risk ST15 CP-Kpn strains in China and warrants rigorous monitoring.

CsrA coordinates the expression of ribosome hibernation and anti-σ factor proteins

Bacterial growth rate varies due to changing physiological signals and is fundamentally dependent on protein synthesis. Consequently, cells alter their transcription and translation machinery to optimize the capacity for protein production under varying conditions and growth rates. Our findings demonstrate that the post-transcriptional regulator CsrA in Escherichia coli controls the expression of genes that participate in these processes. During exponential growth, CsrA represses the expression of proteins that alter or inhibit RNA polymerase (RNAP) and ribosome activity, including the ribosome hibernation factors RMF, RaiA, YqjD, ElaB, YgaM, and SRA, as well as the anti-σ70 factor, Rsd. Upon entry into the stationary phase, RaiA, YqjD, ElaB, and SRA expression was derepressed and that of RMF, YgaM, and Rsd was activated in the presence of CsrA. This pattern of gene expression likely supports global protein expression during active growth and helps limit protein production to a basal level when nutrients are limited. In addition, we identified genes encoding the paralogous C-tail anchored inner membrane proteins YqjD and ElaB as robust, direct targets of CsrA-mediated translational repression. These proteins bind ribosomes and mediate their localization to the inner cell membrane, impacting a variety of processes including protein expression and membrane integrity. Previous studies found that YqjD overexpression inhibits cell growth, suggesting that appropriate regulation of YqjD expression might play a key role in cell viability. CsrA-mediated regulation of yqjD and ribosome hibernation factors reveals a new role for CsrA in appropriating cellular resources for optimum growth under varying conditions.IMPORTANCEThe Csr/Rsm system (carbon storage regulator or repressor of stationary phase metabolites) is a global post-transcriptional regulatory system that coordinates and responds to environmental cues and signals, facilitating the transition between active growth and stationary phase. Another key determinant of bacterial lifestyle decisions is the management of the cellular gene expression machinery. Here, we investigate the connection between these two processes in Escherichia coli. Disrupted regulation of the transcription and translation machinery impacts many cellular functions, including gene expression, growth, fitness, and stress resistance. Elucidating the role of the Csr system in controlling the activity of RNAP and ribosomes advances our understanding of mechanisms controlling bacterial growth. A more complete understanding of these processes could lead to the improvement of therapeutic strategies for recalcitrant infections.

Phylogenetic classification of natural product biosynthetic gene clusters based on regulatory mechanisms

The natural products (NPs) biosynthetic gene clusters (BGCs) represent the adapting biochemical toolkit for microorganisms to thrive different microenvironments. Despite their high diversity, particularly at the genomic level, detecting them in a shake-flask is challenging and remains the primary obstacle limiting our access to valuable chemicals. Studying the molecular mechanisms that regulate BGC expression is crucial to design of artificial conditions that derive on their expression. Here, we propose a phylogenetic analysis of regulatory elements linked to biosynthesis gene clusters, to classify BGCs to regulatory mechanisms based on protein domain information. We utilized Hidden Markov Models from the Pfam database to retrieve regulatory elements, such as histidine kinases and transcription factors, from BGCs in the MIBiG database, focusing on actinobacterial strains from three distinct environments: oligotrophic basins, rainforests, and marine environments. Despite the environmental variations, our isolated microorganisms share similar regulatory mechanisms, suggesting the potential to activate new BGCs using activators known to affect previously characterized BGCs.

Candidatus Kirkpatrickella diaphorinae gen. nov., sp. nov., an uncultured endosymbiont identified in a population of Diaphorina citri from Hawaii

Diaphorina citri is the hemipteran pest and vector of a devastating bacterial pathogen of citrus worldwide. In addition to the two core bacterial endosymbionts of D. citri, Candidatus Carsonella ruddii and Candidatus Profftella armatura, the genome of a novel endosymbiont and as of yet undescribed microbe was discovered in a Hawaiian D. citri population through deep sequencing of multiple D. citri populations. Found to be closely related to the genus Asaia in the family Acetobacteraceae by 16S rRNA gene sequence analysis, it forms a sister clade along with other insect-associated 16S rRNA gene sequences from uncultured bacterium found associated with Aedes koreicus and Sogatella furcifera. Multilocus sequence analysis confirmed the phylogenetic placement sister to the Asaia clade. Despite the culturable Asaia clade being the closest phylogenetic neighbour, attempts to culture this newly identified bacterial endosymbiont were unsuccessful. On the basis of these distinct genetic differences, the novel endosymbiont is proposed to be classified into a candidate genus and species 'Candidatus Kirkpatrickella diaphorinae'. The full genome was deposited in GenBank (accession number CP107052; prokaryotic 16S rRNA OP600170).

Elevation of Lentilactobacillus rapi subsp. dabitei Li et al. 2022 to the species level as Lentilactobacillus dabitei sp. nov

Lentilactobacillus rapi subsp. dabitei was proposed by Li et al. in 2022. The type strains of L. rapi subsp. dabitei and L. rapi subsp. rapi shared 93.1 % average nucleotide identity and 52.8 % digital DNA-DNA hybridization values. Strain IMAU80584T was proposed as a novel subspecies of L. rapi rather than a novel species of the genus Lentilactobacillus on the basis of similar phenotypic characteristics (including growth temperature and pH, tolerance to NaCl and features based on API 50CH and API ZYM). However, the phenotypic investigation performed by Li et al. was insufficient because some physiological and biochemical characteristics recommended by Mattarelli et al. were not included. In the present study, the taxonomic relationship between L. rapi subsp. dabitei and L. rapi subsp. rapi was re-evaluated. Based upon the data obtained in the present study, we propose to elevate L. rapi subsp. dabitei to the species level as Lentilactobacillus dabitei sp. nov. The type strain is IMAU80584T (=GDMCC 1.2566T=JCM 34647T).

An Eco-evolutionary Model on Surviving Lysogeny Through Grounding and Accumulation of Prophages

Temperate phages integrate into the bacterial genomes propagating along with the bacterial genomes. Multiple phage elements, representing diverse prophages, are present in most bacterial genomes. The evolutionary events and the ecological dynamics underlying the accumulation of prophage elements in bacterial genomes have yet to be understood. Here, we show that the local wastewater had 7% of lysogens (hosting mitomycin C-inducible prophages), and they showed resistance to superinfection by their corresponding lysates. Genomic analysis of four lysogens and four non-lysogens revealed the presence of multiple prophages (belonging to Myoviridae and Siphoviridae) in both lysogens and non-lysogens. For large-scale comparison, 2180 Escherichia coli genomes isolated from various sources across the globe and 523 genomes specifically isolated from diverse wastewaters were analyzed. A total of 15,279 prophages were predicted among 2180 E. coli genomes and 2802 prophages among 523 global wastewater isolates, with a mean of ~ 5 prophages per genome. These observations indicate that most putative prophages are relics of past bacteria-phage conflicts; they are "grounded" prophages that cannot excise from the bacterial genome. Prophage distribution analysis based on the sequence homology suggested the random distribution of E. coli prophages within and between E. coli clades. The independent occurrence pattern of these prophages indicates extensive horizontal transfers across the genomes. We modeled the eco-evolutionary dynamics to reconstruct the events that could have resulted in the prophage accumulation accounting for infection, superinfection immunity, and grounding. In bacteria-phage conflicts, the bacteria win by grounding the prophage, which could confer superinfection immunity.

The 2-methylpropene degradation pathway in Mycobacteriaceae family strains

Mycolicibacterium gadium IBE100 and Mycobacterium paragordonae IBE200 are aerobic, chemoorganoheterotrophic bacteria isolated from activated sludge from a wastewater treatment plant. They use 2-methylpropene (isobutene, 2-MP) as the sole source of carbon and energy. Here, we postulate a degradation pathway of 2-methylpropene derived from whole genome sequencing, differential expression analysis and peptide-mass fingerprinting. Key genes identified are coding for a 4-component soluble diiron monooxygenase with epoxidase activity, an epoxide hydrolase, and a 2-hydroxyisobutyryl-CoA mutase. In both strains, involved genes are arranged in clusters of 61.0 and 58.5 kbp, respectively, which also contain the genes coding for parts of the aerobic pathway of adenosylcobalamin synthesis. This vitamin is essential for the carbon rearrangement reaction catalysed by the mutase. These findings provide data for the identification of potential 2-methylpropene degraders.

Achieving quantitative reproducibility in label-free multisite DIA experiments through multirun alignment

DIA is a mainstream method for quantitative proteomics, but consistent quantification across multiple LC-MS/MS instruments remains a bottleneck in parallelizing data acquisition. One reason for this inconsistency and missing quantification is the retention time shift which current software does not adequately address for runs from multiple sites. We present multirun chromatogram alignment strategies to map peaks across columns, including the traditional reference-based Star method, and two novel approaches: MST and Progressive alignment. These reference-free strategies produce a quantitatively accurate data-matrix, even from heterogeneous multi-column studies. Progressive alignment also generates merged chromatograms from all runs which has not been previously achieved for LC-MS/MS data. First, we demonstrate the effectiveness of multirun alignment strategies on a gold-standard annotated dataset, resulting in a threefold reduction in quantitation error-rate compared to non-aligned DIA results. Subsequently, on a multi-species dataset that DIAlignR effectively controls the quantitative error rate, improves precision in protein measurements, and exhibits conservative peak alignment. We next show that the MST alignment reduces cross-site CV by 50% for highly abundant proteins when applied to a dataset from 11 different LC-MS/MS setups. Finally, the reanalysis of 949 plasma runs with multirun alignment revealed a more than 50% increase in insulin resistance (IR) and respiratory viral infection (RVI) proteins, identifying 11 and 13 proteins respectively, compared to prior analysis without it. The three strategies are implemented in our DIAlignR workflow (>2.3) and can be combined with linear, non-linear, or hybrid pairwise alignment.

Urine-mediated suppression of Klebsiella pneumoniae mucoidy is counteracted by spontaneous Wzc variants altering capsule chain length

Klebsiella pneumoniae is a hospital-associated pathogen primarily causing urinary tract infections (UTIs), pneumonia, and septicemia. Two challenging lineages include the hypervirulent strains, causing invasive community-acquired infections, and the carbapenem-resistant classical strains, most frequently isolated from UTIs. While hypervirulent strains are often characterized by a hypermucoid phenotype, classical strains usually present with low mucoidy. Since clinical UTI isolates tend to exhibit limited mucoidy, we hypothesized that environmental conditions may drive K. pneumoniae adaptation to the urinary tract and select against mucoid isolates. We found that both hypervirulent K. pneumoniae and classical Klebsiella UTI isolates significantly suppressed mucoidy when cultured in urine without reducing capsule abundance. A genetic screen identified secondary mutations in the wzc tyrosine kinase that overcome urine-suppressed mucoidy. Over-expressing Wzc variants in trans was sufficient to boost mucoidy in both hypervirulent and classical Klebsiella UTI isolates. Wzc is a bacterial tyrosine kinase that regulates capsule polymerization and extrusion. Although some Wzc variants reduced Wzc phospho-status, urine did not alter Wzc phospho-status. Urine does, however, increase K. pneumoniae capsule chain length diversity and enhance cell-surface attachment. The identified Wzc variants counteract urine-mediated effects on capsule chain length and cell attachment. Combined, these data indicate that capsule chain length correlates with K. pneumoniae mucoidy and that this extracellular feature can be fine-tuned by spontaneous Wzc mutations, which alter host interactions. Spontaneous Wzc mutation represents a global mechanism that could fine-tune K. pneumoniae niche-specific fitness in both classical and hypervirulent isolates. IMPORTANCE Klebsiella pneumoniae is high-priority pathogen causing both hospital-associated infections, such as urinary tract infections, and community-acquired infections. Clinical isolates from community-acquired infection are often characterized by a tacky, hypermucoid phenotype, while urinary tract isolates are usually not mucoid. Historically, mucoidy was attributed to capsule overproduction; however, recent reports have demonstrated that K. pneumoniae capsule abundance and mucoidy are not always correlated. Here, we report that human urine suppresses K. pneumoniae mucoidy, diversifies capsule polysaccharide chain length, and increases cell surface association. Moreover, specific mutations in the capsule biosynthesis gene, wzc, are sufficient to overcome urine-mediated suppression of mucoidy. These Wzc variants cause constitutive production of more uniform capsular polysaccharide chains and increased release of capsule from the cell surface, even in urine. These data demonstrate that K. pneumoniae regulates capsule chain length and cell surface attachment in response host cues, which can alter bacteria-host interactions.

Completed genomes for Rickettsia rickettsii isolated from ticks and quality controlled for motility phenotype

Complete genomes of Rickettsia rickettsii were sequenced with Illumina and PacBio technologies from low-passage isolates from ticks. These isolates were quality controlled for intact roaM, a regulator of actin-based motility that is negatively selected for in culture. The Sheila Smith strain was re-sequenced using the same methodology.

The ClpX protease is essential for inactivating the CI master repressor and completing prophage induction in Staphylococcus aureus

Bacteriophages (phages) are the most abundant biological entities on Earth, exerting a significant influence on the dissemination of bacterial virulence, pathogenicity, and antimicrobial resistance. Temperate phages integrate into the bacterial chromosome in a dormant state through intricate regulatory mechanisms. These mechanisms repress lytic genes while facilitating the expression of integrase and the CI master repressor. Upon bacterial SOS response activation, the CI repressor undergoes auto-cleavage, producing two fragments with the N-terminal domain (NTD) retaining significant DNA-binding ability. The process of relieving CI NTD repression, essential for prophage induction, remains unknown. Here we show a specific interaction between the ClpX protease and CI NTD repressor fragment of phages Ф11 and 80α in Staphylococcus aureus. This interaction is necessary and sufficient for prophage activation after SOS-mediated CI auto-cleavage, defining the final stage in the prophage induction cascade. Our findings unveil unexpected roles of bacterial protease ClpX in phage biology.

σ28-dependent small RNA regulation of flagella biosynthesis

Flagella are important for bacterial motility as well as for pathogenesis. Synthesis of these structures is energy intensive and, while extensive transcriptional regulation has been described, little is known about the posttranscriptional regulation. Small RNAs (sRNAs) are widespread posttranscriptional regulators, most base pairing with mRNAs to affect their stability and/or translation. Here, we describe four UTR-derived sRNAs (UhpU, MotR, FliX and FlgO) whose expression is controlled by the flagella sigma factor σ28 (fliA) in Escherichia coli. Interestingly, the four sRNAs have varied effects on flagellin protein levels, flagella number and cell motility. UhpU, corresponding to the 3´ UTR of a metabolic gene, likely has hundreds of targets including a transcriptional regulator at the top flagella regulatory cascade connecting metabolism and flagella synthesis. Unlike most sRNAs, MotR and FliX base pair within the coding sequences of target mRNAs and act on ribosomal protein mRNAs connecting ribosome production and flagella synthesis. The study shows how sRNA-mediated regulation can overlay a complex network enabling nuanced control of flagella synthesis.

Comparative genomics of Loigolactobacillus coryniformis with an emphasis on L. coryniformis strain FOL-19 isolated from cheese

Loigolactobacillus coryniformis is a member of lactic acid bacteria isolated from various ecological niches. We isolated a novel L. coryniformis strain FOL-19 from artisanal Tulum cheese and performed the whole-genome sequencing for FOL-19. Then, genomic characterization of FOL-19 against ten available whole genome sequences of the same species isolated from kimchi, silage, fermented meat, air of cowshed, dairy, and pheasant chyme was performed to uncover the genetic diversity and biotechnological potential of overall species. The average genome size of 2.93 ± 0.1 Mb, GC content of 42.96% ± 0.002, number of CDS of 2905 ± 165, number of tRNA of 56 ± 10, and number of CRISPR elements of 6.55 ± 1.83 was found. Both Type I and II Cas clusters were observed in L. coryniformis. No bacteriocin biosynthesis gene clusters were found. All strains harbored at least one plasmid except KCTC 3167. All strains were predicted to carry multiple IS elements. The most common origin of the IS elements was belong to Lactiplantibacillus plantarum. Comparative genomic analysis of L. coryniformis revealed hypervariability at the strain level and the presence of CRISPR/Cas suggests that L. coryniformis holds a promising potential for being a reservoir for new CRISPR-based tools. All L. coryniformis strains except PH-1 were predicted to harbor pdu and cbi-cob-hem gene clusters encoding industrially relevant traits of reuterin and cobalamin biosynthesis, respectively. These findings put a step forward for the genomic characterization of L. coryniformis strains for biotechnological applications via genome-guided strain selection to identify industrially relevant traits.

Chemoautotrophic production of gaseous hydrocarbons, bioplastics and osmolytes by a novel Halomonas species

BACKGROUND

Production of relatively low value, bulk commodity chemicals and fuels by microbial species requires a step-change in approach to decrease the capital and operational costs associated with scaled fermentation. The utilisation of the robust and halophilic industrial host organisms of the genus Halomonas could dramatically decrease biomanufacturing costs owing to their ability to grow in seawater, using waste biogenic feedstocks, under non-sterile conditions.

RESULTS

We describe the isolation of Halomonas rowanensis, a novel facultative chemoautotrophic species of Halomonas from a natural brine spring. We investigated the ability of this species to produce ectoine, a compound of considerable industrial interest, under heterotrophic conditions. Fixation of radiolabelled NaH14CO3 by H. rowanensis was confirmed in mineral medium supplied with thiosulfate as an energy source. Genome sequencing suggested carbon fixation proceeds via a reductive tricarboxylic acid cycle, and not the Calvin-Bensen-Bassham cycle. The mechanism of energy generation to support chemoautotrophy is unknown owing to the absence of an annotated SOX-based thiosulfate-mediated energy conversion system. We investigated further the biotechnological potential of the isolated H. rowanensis by demonstrating production of the gaseous hydrocarbon (bio-propane), bioplastics (poly-3-hydroxybutyrate) and osmolytes (ectoine) under heterotrophic and autotrophic CO2 fixation growth conditions.

CONCLUSIONS

This proof-of-concept study illustrates the value of recruiting environmental isolates as industrial hosts for chemicals biomanufacturing, where CO2 utilisation could replace, or augment, the use of biogenic feedstocks in non-sterile, industrialised bioreactors.

Mycobacteriaceae Phenome Atlas (MPA): A Standardized Atlas for the Mycobacteriaceae Phenome Based on Heterogeneous Sources

The bacterial family Mycobacteriaceae includes pathogenic and nonpathogenic bacteria, and systematic research on their genome and phenome can give comprehensive perspectives for exploring their disease mechanism. In this study, the phenotypes of Mycobacteriaceae were inferred from available phenomic data, and 82 microbial phenotypic traits were recruited as data elements of the microbial phenome. This Mycobacteriaceae phenome contains five categories and 20 subcategories of polyphasic phenotypes, and three categories and eight subcategories of functional phenotypes, all of which are complementary to the existing data standards of microbial phenotypes. The phenomic data of Mycobacteriaceae strains were compiled by literature mining, third-party database integration, and bioinformatics annotation. The phenotypes were searchable and comparable from the website of the Mycobacteriaceae Phenome Atlas (MPA, https://www.biosino.org/mpa/). A topological data analysis of MPA revealed the co-evolution between Mycobacterium tuberculosis and virulence factors, and uncovered potential pathogenicity-associated phenotypes. Two hundred and sixty potential pathogen-enriched pathways were found by Fisher's exact test. The application of MPA may provide novel insights into the pathogenicity mechanism and antimicrobial targets of Mycobacteriaceae.

Supplementary Information

The online version contains supplementary material available at 10.1007/s43657-023-00101-5.

RVFScan predicts virulence factor genes and hypervirulence of the clinical metagenome

Bacterial infections often involve virulence factors that play a crucial role in the pathogenicity of bacteria. Accurate detection of virulence factor genes (VFGs) is essential for precise treatment and prognostic management of hypervirulent bacterial infections. However, there is a lack of rapid and accurate methods for VFG identification from the metagenomic data of clinical samples. Here, we developed a Reads-based Virulence Factors Scanner (RVFScan), an innovative user-friendly online tool that integrates a comprehensive VFG database with similarity matrix-based criteria for VFG prediction and annotation using metagenomic data without the need for assembly. RVFScan demonstrated superior performance compared to previous assembly-based and read-based VFG predictors, achieving a sensitivity of 97%, specificity of 98% and accuracy of 98%. We also conducted a large-scale analysis of 2425 clinical metagenomic datasets to investigate the utility of RVFScan, the species-specific VFG profiles and associations between VFGs and virulence phenotypes for 24 important pathogens were analyzed. By combining genomic comparisons and network analysis, we identified 53 VFGs with significantly higher abundances in hypervirulent Klebsiella pneumoniae (hvKp) than in classical K. pneumoniae. Furthermore, a cohort of 1256 samples suspected of K. pneumoniae infection demonstrated that RVFScan could identify hvKp with a sensitivity of 90%, specificity of 100% and accuracy of 98.73%, with 90% of hvKp samples consistent with clinical diagnosis (Cohen's kappa, 0.94). RVFScan has the potential to detect VFGs in low-biomass and high-complexity clinical samples using metagenomic reads without assembly. This capability facilitates the rapid identification and targeted treatment of hvKp infections and holds promise for application to other hypervirulent pathogens.

Draft genome sequences of Pantoea sp. strains QMID1-QMID4 isolated from the midgut of Japanese honey bee (Apis cerana japonica)

We report the draft genome sequences of Pantoea sp. strains QMID1-QMID4 that were recovered from the midgut of Japanese honey bee (Apis cerana japonica). The strains possess the carotenoid biosynthetic gene cluster. The genome information expands our knowledge of their potential use as probiotics and/or prebiotics in honey bees.

Whole-genome sequencing and annotation of selected Enterobacteriaceae derived from commercial cucumber fermentation

We report five whole-genome sequences, along with annotations, representing Enterobacteriaceae from the genera Enterobacter (n = 1), Pantoea (n = 1), and Leclercia (n = 3) isolated from commercial cucumber fermentations performed in North Carolina and Minnesota, USA.

T9GPred: A Comprehensive Computational Tool for the Prediction of Type 9 Secretion System, Gliding Motility, and the Associated Secreted Proteins

Type 9 secretion system (T9SS) is one of the least characterized secretion systems exclusively found in the Bacteroidetes phylum, which comprises various environmental and economically relevant bacteria. While T9SS plays a central role in bacterial movement termed gliding motility, survival, and pathogenicity, there is an unmet need for a comprehensive tool that predicts T9SS, gliding motility, and proteins secreted via T9SS. In this study, we develop such a computational tool, Type 9 secretion system and Gliding motility Prediction (T9GPred). To build this tool, we manually curated published experimental evidence and identified mandatory components for T9SS and gliding motility prediction. We also compiled experimentally characterized proteins secreted via T9SS and determined the presence of three unique types of C-terminal domain signals, and these insights were leveraged to predict proteins secreted via T9SS. Notably, using recently published experimental evidence, we show that T9GPred has high predictive power. Thus, we used T9GPred to predict the presence of T9SS, gliding motility, and associated secreted proteins across 693 completely sequenced Bacteroidetes strains. T9GPred predicted 402 strains to have T9SS, of which 327 strains are also predicted to exhibit gliding motility. Further, T9GPred also predicted putative secreted proteins for the 402 strains. In a nutshell, T9GPred is a novel computational tool for systems-level prediction of T9SS and streamlining future experimentation. The source code of the computational tool is available in our GitHub repository: https://github.com/asamallab/T9GPred. The tool and its predicted results are compiled in a web server available at: https://cb.imsc.res.in/t9gpred/.

Single-strain behavior predicts responses to environmental pH and osmolality in the gut microbiota

Changes to gut environmental factors such as pH and osmolality due to disease or drugs correlate with major shifts in microbiome composition; however, we currently cannot predict which species can tolerate such changes or how the community will be affected. Here, we assessed the growth of 92 representative human gut bacterial strains spanning 28 families across multiple pH values and osmolalities in vitro. The ability to grow in extreme pH or osmolality conditions correlated with the availability of known stress response genes in many cases, but not all, indicating that novel pathways may participate in protecting against acid or osmotic stresses. Machine learning analysis uncovered genes or subsystems that are predictive of differential tolerance in either acid or osmotic stress. For osmotic stress, we corroborated the increased abundance of these genes in vivo during osmotic perturbation. The growth of specific taxa in limiting conditions in isolation in vitro correlated with survival in complex communities in vitro and in an in vivo mouse model of diet-induced intestinal acidification. Our data show that in vitro stress tolerance results are generalizable and that physical parameters may supersede interspecies interactions in determining the relative abundance of community members. This study provides insight into the ability of the microbiota to respond to common perturbations that may be encountered in the gut and provides a list of genes that correlate with increased ability to survive in these conditions. IMPORTANCE To achieve greater predictability in microbiota studies, it is crucial to consider physical environmental factors such as pH and particle concentration, as they play a pivotal role in influencing bacterial function and survival. For example, pH is significantly altered in various diseases, including cancers, inflammatory bowel disease, as well in the case of over-the-counter drug use. Additionally, conditions like malabsorption can affect particle concentration. In our study, we investigate how changes in environmental pH and osmolality can serve as predictive indicators of bacterial growth and abundance. Our research provides a comprehensive resource for anticipating shifts in microbial composition and gene abundance during complex perturbations. Moreover, our findings underscore the significance of the physical environment as a major driver of bacterial composition. Finally, this work emphasizes the necessity of incorporating physical measurements into animal and clinical studies to better understand the factors influencing shifts in microbiota abundance.

Predicting variable gene content in Escherichia coli using conserved genes

Having the ability to predict the protein-encoding gene content of an incomplete genome or metagenome-assembled genome is important for a variety of bioinformatic tasks. In this study, as a proof of concept, we built machine learning classifiers for predicting variable gene content in Escherichia coli genomes using only the nucleotide k-mers from a set of 100 conserved genes as features. Protein families were used to define orthologs, and a single classifier was built for predicting the presence or absence of each protein family occurring in 10%-90% of all E. coli genomes. The resulting set of 3,259 extreme gradient boosting classifiers had a per-genome average macro F1 score of 0.944 [0.943-0.945, 95% CI]. We show that the F1 scores are stable across multi-locus sequence types and that the trend can be recapitulated by sampling a smaller number of core genes or diverse input genomes. Surprisingly, the presence or absence of poorly annotated proteins, including "hypothetical proteins" was accurately predicted (F1 = 0.902 [0.898-0.906, 95% CI]). Models for proteins with horizontal gene transfer-related functions had slightly lower F1 scores but were still accurate (F1s = 0.895, 0.872, 0.824, and 0.841 for transposon, phage, plasmid, and antimicrobial resistance-related functions, respectively). Finally, using a holdout set of 419 diverse E. coli genomes that were isolated from freshwater environmental sources, we observed an average per-genome F1 score of 0.880 [0.876-0.883, 95% CI], demonstrating the extensibility of the models. Overall, this study provides a framework for predicting variable gene content using a limited amount of input sequence data. IMPORTANCE Having the ability to predict the protein-encoding gene content of a genome is important for assessing genome quality, binning genomes from shotgun metagenomic assemblies, and assessing risk due to the presence of antimicrobial resistance and other virulence genes. In this study, we built a set of binary classifiers for predicting the presence or absence of variable genes occurring in 10%-90% of all publicly available E. coli genomes. Overall, the results show that a large portion of the E. coli variable gene content can be predicted with high accuracy, including genes with functions relating to horizontal gene transfer. This study offers a strategy for predicting gene content using limited input sequence data.

Unveiling the microbial diversity and functional dynamics of Shiv Kund, Sohna hot spring, India through a shotgun metagenomics approach

In this research, we examined the microbial diversity in Sohna hot spring, Haryana, India using shotgun metagenome sequencing based on the Illumina Hiseq 4000 sequencing technology. The raw sequence data from metagenomic paired-end libraries were analysed for taxonomic classification, diversity, and functional annotation using MG-RAST online server. The results showed the presence of total of 57 phyla, 931 genera, and 2068 species, predominantly occupied by Moraxellaceae (Gammaproteobacteria). However, at the species level, we reported the presence of some representative pathogenic taxa, such as Acinetobacter baumannii and Moraxella osloensis. The functional annotation predicted at various levels based on SEED-based subsystem, KEGG ortholog identity (KO), Cluster of Orthologous Groups (COGs) database identified the predominance of genes associated with primary and secondary metabolism along with a crucial role in environmental and genetic signals, cellular communication, and cell signalling. Comparative Genome Analysis (CGA) using The Pathosystem Resource Integration Centre (PATRIC) tool based on genome annotation and assembly of the metagenomic libraries for representative taxon Acinetobacter baumannii (NCBI tax id:470) characterized the reads with a unique genome identifier of 470.20380 (A. baumannii DDLJ4) which is evolutionary closer to A. baumannii ATCC 470.17978 400667.7. In addition, the CARD database results about the presence of potential AMR pathotypes and the prevalence of adeABC, adeIJK, abeM gene-specific clusters that function as multidrug efflux pumps. Overall, the results provided a comprehensive insight into virulence and anti-microbial resistance mechanism and could be useful for developing potential drug targets against the possible AMR pathotypes.

A modular plasmid toolkit applied in marine bacteria reveals functional insights during bacteria-stimulated metamorphosis

A conspicuous roadblock to studying marine bacteria for fundamental research and biotechnology is a lack of modular synthetic biology tools for their genetic manipulation. Here, we applied, and generated new parts for, a modular plasmid toolkit to study marine bacteria in the context of symbioses and host-microbe interactions. To demonstrate the utility of this plasmid system, we genetically manipulated the marine bacterium Pseudoalteromonas luteoviolacea, which stimulates the metamorphosis of the model tubeworm, Hydroides elegans. Using these tools, we quantified constitutive and native promoter expression, developed reporter strains that enable the imaging of host-bacteria interactions, and used CRISPR interference (CRISPRi) to knock down a secondary metabolite and a host-associated gene. We demonstrate the broader utility of this modular system for testing the genetic tractability of marine bacteria that are known to be associated with diverse host-microbe symbioses. These efforts resulted in the successful conjugation of 12 marine strains from the Alphaproteobacteria and Gammaproteobacteria classes. Altogether, the present study demonstrates how synthetic biology strategies enable the investigation of marine microbes and marine host-microbe symbioses with potential implications for environmental restoration and biotechnology. IMPORTANCE Marine Proteobacteria are attractive targets for genetic engineering due to their ability to produce a diversity of bioactive metabolites and their involvement in host-microbe symbioses. Modular cloning toolkits have become a standard for engineering model microbes, such as Escherichia coli, because they enable innumerable mix-and-match DNA assembly and engineering options. However, such modular tools have not yet been applied to most marine bacterial species. In this work, we adapt a modular plasmid toolkit for use in a set of 12 marine bacteria from the Gammaproteobacteria and Alphaproteobacteria classes. We demonstrate the utility of this genetic toolkit by engineering a marine Pseudoalteromonas bacterium to study their association with its host animal Hydroides elegans. This work provides a proof of concept that modular genetic tools can be applied to diverse marine bacteria to address basic science questions and for biotechnology innovations.

Prophage identification and molecular analysis in the genomes of Pseudomonas aeruginosa strains isolated from critical care patients

Prophages are bacteriophages integrated into the bacterial host's chromosome. This research aims to analyze and characterize the existing prophages within a collection of 53 Pseudomonas aeruginosa strains from intensive care units (ICUs) in Portugal and Spain. A total of 113 prophages were localized in the collection, with 18 of them being present in more than one strain simultaneously. After annotation, five of them were discarded as incomplete, and the 13 remaining prophages were characterized. Of 13, 10 belonged to the siphovirus tail morphology group, 2 to the podovirus tail morphology group, and 1 to the myovirus tail morphology group. All prophages had a length ranging from 20,199 to 63,401 bp and a GC% between 56.2% and 63.6%. The number of open reading frames (ORFs) oscillated between 32 and 88, and in 3/13 prophages, more than 50% of the ORFs had an unknown function. With our findings, we show that prophages are present in the majority of the P. aeruginosa strains isolated from Portuguese and Spanish critically ill patients, many of them found in more than one circulating strain at the same time and following a similar clonal distribution pattern. Although a great sum of ORFs had an unknown function, number of proteins in relation to viral defense (anti-CRISPR proteins, toxin/antitoxin modules, proteins against restriction-modification systems) as well as to prophage interference into their host's quorum sensing system and regulatory cascades were found. This supports the idea that prophages have an influence in bacterial pathogenesis and anti-phage defense. IMPORTANCE Despite being known for decades, prophages remain understudied when compared to the lytic phages employed in phage therapy. This research aims to shed some light into the nature, composition, and role of prophages found within a set of circulating strains of Pseudomas aeruginosa, with special attention to high-risk clones. Given the fact that prophages can effectively influence bacterial pathogenesis, prophage basic research constitutes a topic of growing interest. Furthermore, the abundance of viral defense and regulatory proteins within prophage genomes detected in this study evidences the importance of characterizing the most frequent prophages in circulating clinical strains and in high-risk clones if phage therapy is to be used.

Distinct Origins and Transmission Pathways of blaKPC Enterobacterales across Three U.S. States

Carbapenem-resistant Enterobacterales (CRE) are among the most concerning antibiotic resistance threats due to high rates of multidrug resistance, transmissibility in health care settings, and high mortality rates. We evaluated the potential for regional genomic surveillance to track the spread of blaKPC-carrying CRE (KPC-CRE) by using isolate collections from health care facilities in three U.S. states. Clinical isolates were collected from Connecticut (2017 to 2018), Minnesota (2012 to 2018), and Tennessee (2016 to 2017) through the U.S. Centers for Disease Control and Prevention's Multi-site Gram-negative Surveillance Initiative (MuGSI) and additional surveillance. KPC-CRE isolates were whole-genome sequenced, yielding 255 isolates from 214 patients across 96 facilities. Case report data on patient comorbidities, facility exposures, and interfacility patient transfer were extracted. We observed that in Connecticut, most KPC-CRE isolates showed evidence of importation from outside the state, with limited local transmission. In Minnesota, cases were mainly from sporadic importation and transmission of blaKPC-carrying Klebsiella pneumoniae ST258, and clonal expansion of blaKPC-carrying Enterobacter hormaechei ST171, primarily at a single focal facility and its satellite facilities. In Tennessee, we observed transmission of diverse strains of blaKPC-carrying Enterobacter and Klesbiella, with evidence that most derived from the local acquisition of blaKPC plasmids circulating in an interconnected regional health care network. Thus, the underlying processes driving KPC-CRE burden can differ substantially across regions and can be discerned through regional genomic surveillance. This study provides proof of concept that integrating genomic data with information on interfacility patient transfers can provide insights into locations and drivers of regional KPC-CRE burden that can enable targeted interventions.

Description of Sporanaerobium hydrogeniformans gen. nov., sp. nov., an obligately anaerobic, hydrogen-producing bacterium isolated from Aravali hot spring in India

An obligately anaerobic bacterium XHS1971T, capable of degrading cellulose and xylan, was isolated from a sediment sample of Aravali hot spring, Ratnagiri, India. Cells of strain XHS1971T were Gram-stain-negative, spore-forming, motile, long-rods. Growth was observed at temperatures 30-50 °C (optimum 40-45 °C), pH 5.0-10.0 (optimum pH 8.0) and NaCl concentrations 0-0.5% (optimum 0%). Generation time of strain XHS1971T was 5 h under optimised growth conditions. Strain XHS1971T showed the ability to metabolise different complex and simple sugars constituting lignocellulosic biomass. Glucose was fermented majorly into hydrogen, formic acid, acetic acid, and ethanol, whereas carbon dioxide, butyric acid, lactic acid and succinic acid were produced in traces. 16S rRNA gene analysis of strain XHS1971T revealed < 94.5% homology with Cellulosilyticum lentocellum DSM5427T followed by Cellulosilyticum ruminicola JCM14822T, identifying strain as a distinct member of family Lachnospiraceae. The major cellular fatty acids (> 5%) were C14:0, C16:0, C18:0, and C16:1 ω7c. The genome size of the strain was 3.74 Mb with 35.3 mol% G + C content, and genes were annotated to carbohydrate metabolism, including genes involved in the degradation of cellulose and xylan and the production of hydrogen, ethanol and acetate. The uniqueness of strain was further validated by digital DNA-DNA hybridisation (dDDH), Average Nucleotide Identity (ANI), and Average Amino Acid Identity (AAI) values of 22%, 80%, and 63%, respectively, with nearest phylogenetic affiliates. Based on the detailed analyses, we propose a new genus and species, Sporanaerobium hydrogeniformans gen. nov., sp. nov., for strain XHS1971T (= MCC3498T = KCTC15729T = JCM32657T) within family Lachnospiraceae.

Emergent community architecture despite distinct diversity in the global whale shark (Rhincodon typus) epidermal microbiome

Microbiomes confer beneficial physiological traits to their host, but microbial diversity is inherently variable, challenging the relationship between microbes and their contribution to host health. Here, we compare the diversity and architectural complexity of the epidermal microbiome from 74 individual whale sharks (Rhincodon typus) across five aggregations globally to determine if network properties may be more indicative of the microbiome-host relationship. On the premise that microbes are expected to exhibit biogeographic patterns globally and that distantly related microbial groups can perform similar functions, we hypothesized that microbiome co-occurrence patterns would occur independently of diversity trends and that keystone microbes would vary across locations. We found that whale shark aggregation was the most important factor in discriminating taxonomic diversity patterns. Further, microbiome network architecture was similar across all aggregations, with degree distributions matching Erdos-Renyi-type networks. The microbiome-derived networks, however, display modularity indicating a definitive microbiome structure on the epidermis of whale sharks. In addition, whale sharks hosted 35 high-quality metagenome assembled genomes (MAGs) of which 25 were present from all sample locations, termed the abundant 'core'. Two main MAG groups formed, defined here as Ecogroup 1 and 2, based on the number of genes present in metabolic pathways, suggesting there are at least two important metabolic niches within the whale shark microbiome. Therefore, while variability in microbiome diversity is high, network structure and core taxa are inherent characteristics of the epidermal microbiome in whale sharks. We suggest the host-microbiome and microbe-microbe interactions that drive the self-assembly of the microbiome help support a functionally redundant abundant core and that network characteristics should be considered when linking microbiomes with host health.

Latent evolution of biofilm formation depends on life-history and genetic background

Adaptation to one environment can often generate phenotypic and genotypic changes which impact the future ability of an organism to thrive in other environmental conditions. In the context of host-microbe interactions, biofilm formation can increase survival rates in vivo upon exposure to stresses, like the host's immune system or antibiotic therapy. However, how the generic process of adaptation impacts the ability to form biofilm and how it may change through time has seldomly been studied. To do so, we used a previous evolution experiment with three strains of the Klebsiella pneumoniae species complex, in which we specifically did not select for biofilm formation. We observed that changes in the ability to form biofilm happened very fast at first and afterwards reverted to ancestral levels in many populations. Biofilm changes were associated to changes in population yield and surface polysaccharide production. Genotypically, mutations in the tip adhesin of type III fimbriae (mrkD) or the fim switch of type I fimbriae were shaped by nutrient availability during evolution, and their impact on biofilm formation was dependent on capsule production. Analyses of natural isolates revealed similar mutations in mrkD, suggesting that such mutations also play an important role in adaptation outside the laboratory. Our work reveals that the latent evolution of biofilm formation, and its temporal dynamics, depend on nutrient availability, the genetic background and other intertwined phenotypic and genotypic changes. Ultimately, it suggests that small differences in the environment can alter an organism's fate in more complex niches like the host.

Extensive screening reveals previously undiscovered aminoglycoside resistance genes in human pathogens

Antibiotic resistance is a growing threat to human health, caused in part by pathogens accumulating antibiotic resistance genes (ARGs) through horizontal gene transfer. New ARGs are typically not recognized until they have become widely disseminated, which limits our ability to reduce their spread. In this study, we use large-scale computational screening of bacterial genomes to identify previously undiscovered mobile ARGs in pathogens. From ~1 million genomes, we predict 1,071,815 genes encoding 34,053 unique aminoglycoside-modifying enzymes (AMEs). These cluster into 7,612 families (<70% amino acid identity) of which 88 are previously described. Fifty new AME families are associated with mobile genetic elements and pathogenic hosts. From these, 24 of 28 experimentally tested AMEs confer resistance to aminoglycoside(s) in Escherichia coli, with 17 providing resistance above clinical breakpoints. This study greatly expands the range of clinically relevant aminoglycoside resistance determinants and demonstrates that computational methods enable early discovery of potentially emerging ARGs.

Genome-wide and constrained ordination-based analyses of EC code data support reclassification of the species of Massilia La Scola et al. 2000 into Telluria Bowman et al. 1993, Mokoshia gen. nov. and Zemynaea gen. nov

Based on genome-wide data, Massilia species belonging to the clade including Telluria mixta LMG 11547T should be entirely transferred to the genus Telluria owing to the nomenclatural priority of the type species Telluria mixta. This results in the transfer of 35 Massilia species to the genus Telluria. The presented data also supports the creation of two new genera since peripherally branching Massilia species are distinct from Telluria and other related genera. It is proposed that 13 Massilia species are transferred to Mokoshia gen. nov. with the type species designated Mokoshia eurypsychrophila comb. nov. The species Massilia arenosa is proposed to belong to the genus Zemynaea gen. nov. as the type species Zemynaea arenosa comb. nov. The genome-wide analysis was well supported by canonical ordination analysis of Enzyme Commission (EC) codes annotated from genomes via pannzer2. This new approach was performed to assess the conclusions of the genome-based data and reduce possible ambiguity in the taxonomic decision making. Cross-validation of EC code data compared within canonical plots validated the reclassifications and correctly visualized the expected genus-level taxonomic relationships. The approach is complementary to genome-wide methodology and could be used for testing sequence alignment based data across genetically related genera. In addition to the proposed broader reclassifications, invalidly described species 'Massilia antibiotica', 'Massilia aromaticivorans', 'Massilia cellulosiltytica' and 'Massilia humi' are described as Telluria antibiotica sp. nov., Telluria aromaticivorans sp. nov., Telluria cellulosilytica sp. nov. and Pseudoduganella humi sp. nov., respectively. In addition, Telluria chitinolytica is reclassified as Pseudoduganella chitinolytica comb. nov. The use of combined genome-wide and annotation descriptors compared using canonical ordination clarifies the taxonomy of Telluria and its sibling genera and provides another way to evaluate complex taxonomic data.

Levilactobacillus brevis, autochthonous to cucumber fermentation, is unable to utilize citric acid and encodes for a putative 1,2-propanediol utilization microcompartment

The metabolic versatility of Levilactobacillus brevis, a heterofermentative lactic acid bacterium, could benefit environmentally compatible and low salt cucumber fermentation. The biodiversity of Lvb. brevis autochthonous to cucumber fermentation was studied using genotypic and phenotypic analyses to identify unique adjunct cultures. A group of 131 isolates autochthonous to industrial fermentations was screened using rep-PCR-(GTG)5 and a fermentation ability assay under varied combinations of salt (0 or 6%), initial pH (4.0 or 5.2), and temperature (15 or 30°C). No apparent similarities were observed among the seven and nine clusters in the genotypic and phenotypic dendrograms, respectively. A total of 14 isolates representing the observed biodiversity were subjected to comparative genome analysis. The autochthonous Lvb. brevis clustered apart from allochthonous isolates, as their genomes lack templates for citrate lyase, several putative hypothetical proteins, and some plasmid- and phage-associated proteins. Four and two representative autochthonous and allochthonous Lvb. brevis, respectively, were subjected to phenotype microarray analysis using an Omnilog. Growth of all Lvb. brevis strains was supported to various levels by glucose, fructose, gentiobiose, 1,2-propanediol, and propionic acid, whereas the allochthonous isolate ATCC14890 was unique in utilizing citric acid. All the Lvb. brevis genomes encode for 1,2-propanediol utilization microcompartments. This study identified a unique Lvb. brevis strain, autochthonous to cucumber, as a potential functional adjunct culture for commercial fermentation that is distinct in metabolic activities from allochthonous isolates of the same species.

Molecular characterization and in-depth genome analysis of Enterobacter sp. S-16

Enterobacter species are considered to be an opportunistic human pathogen owing to the existence of antibiotic-resistant strains and drug resides; however, the detailed analysis of the antibiotic resistance and virulence features in environmental isolates is poorly characterized. Here, in the study, we characterized the biochemical characteristics, and genome, pan-genome, and comparative genome analyses of an environmental isolate Enterobacter sp. S-16. The strain was identified as Enterobacter spp. by using 16S rRNA gene sequencing. To unravel genomic features, whole genome of Enterobacter sp. S-16 was sequenced using a hybrid assembly approach and genome assembly was performed using the Unicycler tool. The assembled genome contained the single conting size 5.3 Mbp, GC content 55.43%, and 4500 protein-coding genes. The genome analysis revealed the various gene clusters associated with virulence, antibiotic resistance, type VI secretion system (T6SS), and many stress tolerant genes, which may provide important insight for adapting to changing environment conditions. Moreover, different metabolic pathways were identified that potentially contribute to environmental survival. Various hydrolytic enzymes and motility functions equipped the strain S-16 as an active colonizer. The genome analysis confirms the presence of carbohydrate-active enzymes (CAZymes), and non-enzymatic carbohydrate-binding modules (CBMs) involved in the hydrolysis of complex carbohydrate polymers. Moreover, the pan-genome analysis provides detailed information about the core genes and shared genes with the closest related Enterobacter species. The present study is the first report showing the presence of YdhE/NorM in Enterobacter spp. Thus, the elucidation of genome sequencing will increase our understanding of the pathogenic nature of environmental isolate, supporting the One Health Concept.

Novel complete methanogenic pathways in longitudinal genomic study of monogastric age-associated archaea

BACKGROUND

Archaea perform critical roles in the microbiome system, including utilizing hydrogen to allow for enhanced microbiome member growth and influencing overall host health. With the majority of microbiome research focusing on bacteria, the functions of archaea are largely still under investigation. Understanding methanogenic functions during the host lifetime will add to the limited knowledge on archaeal influence on gut and host health. In our study, we determined lifelong archaea dynamics, including detection and methanogenic functions, while assessing global, temporal and host distribution of our novel archaeal metagenome-assembled genomes (MAGs). We followed 7 monogastric swine throughout their life, from birth to adult (1-156 days of age), and collected feces at 22 time points. The samples underwent gDNA extraction, Illumina sequencing, bioinformatic quality and assembly processes, MAG taxonomic assignment and functional annotation. MAGs were utilized in downstream phylogenetic analysis for global, temporal and host distribution in addition to methanogenic functional potential determination.

RESULTS

We generated 1130 non-redundant MAGs, representing 588 unique taxa at the species level, with 8 classified as methanogenic archaea. The taxonomic classifications were as follows: orders Methanomassiliicoccales (5) and Methanobacteriales (3); genera UBA71 (3), Methanomethylophilus (1), MX-02 (1), and Methanobrevibacter (3). We recovered the first US swine Methanobrevibacter UBA71 sp006954425 and Methanobrevibacter gottschalkii MAGs. The Methanobacteriales MAGs were identified primarily during the young, preweaned host whereas Methanomassiliicoccales primarily in the adult host. Moreover, we identified our methanogens in metagenomic sequences from Chinese swine, US adult humans, Mexican adult humans, Swedish adult humans, and paleontological humans, indicating that methanogens span different hosts, geography and time. We determined complete metabolic pathways for all three methanogenic pathways: hydrogenotrophic, methylotrophic, and acetoclastic. This study provided the first evidence of acetoclastic methanogenesis in archaea of monogastric hosts which indicated a previously unknown capability for acetate utilization in methanogenesis for monogastric methanogens. Overall, we hypothesized that the age-associated detection patterns were due to differential substrate availability via the host diet and microbial metabolism, and that these methanogenic functions are likely crucial to methanogens across hosts. This study provided a comprehensive, genome-centric investigation of monogastric-associated methanogens which will further improve our understanding of microbiome development and functions.

Complete genome sequence of probiotic Lactobacillus johnsonii 7409N31 isolated from a healthy Hanwoo calf

Lactobacillus johnsonii 7409N31 was isolated from the feces of a healthy 11-day-old Hanwoo calf from a farm in Geochang-gun, Gyeongsangnam-do, Korea. The genome of the strain was completely sequenced using the PacBio RSII sequencing system, and it was confirmed that it was composed of one circular chromosome. The size of the entire genome was 2,198,442 bp, and it had 35.01 mol% guanine + cytosine (G + C) content and 2,222 protein-coding sequences, 24 rRNA, 3 ncRNA, and 112 tRNA genes. Strain 7409N31 possessed genes encoding enzymes involved in the hydrolysis of both fibrous and non-fibrous carbohydrates. These data provide a comprehensive theoretical understanding for developing industrial probiotic feed additives that improve nutrient digestibility.

Metagenomics of culture isolates and insect tissue illuminate the evolution of Wolbachia, Rickettsia and Bartonella symbionts in Ctenocephalides spp. fleas

While fleas are often perceived simply as a biting nuisance and a cause of allergic dermatitis, they represent important disease vectors worldwide, especially for bacterial zoonoses such as plague (transmitted by rodent fleas) and some of the rickettsioses and bartonelloses. The cosmopolitan cat (Ctenocephalides felis ) and dog (Ctenocephalides canis ) fleas, as well as Ctenocephalides orientis (restricted to tropical and subtropical Asia), breed in human dwellings and are vectors of cat-scratch fever (caused by Bartonella spp.) and Rickettsia spp., including Rickettsia felis (agent of flea-borne spotted fever) and Rickettsia asembonensis , a suspected pathogen. These Rickettsia spp. are members of a phylogenetic clade known as the ‘transitional group’, which includes both human pathogens and arthropod-specific endosymbionts. The relatively depauperate flea microbiome can also contain other endosymbionts, including a diverse range of Wolbachia strains. Here, we present circularized genome assemblies for two C. orientis -derived pathogens (Bartonella clarridgeiae and R. asembonensis ) from Malaysia, a novel Wolbachia strain (w Cori), and the C. orientis mitochondrion; all were obtained by direct metagenomic sequencing of flea tissues. Moreover, we isolated two Wolbachia strains from Malaysian C. felis into tick cell culture and recovered circularized genome assemblies for both, one of which (w CfeF) is newly sequenced. We demonstrate that the three Wolbachia strains are representatives of different major clades (‘supergroups’), two of which appear to be flea-specific. These Wolbachia genomes exhibit unique combinations of features associated with reproductive parasitism or mutualism, including prophage WO, cytoplasmic incompatibility factors and the biotin operon of obligate intracellular microbes. The first circularized assembly for R. asembonensis includes a plasmid with a markedly different structure and gene content compared to the published plasmid; moreover, this novel plasmid was also detected in cat flea metagenomes from the USA. Analysis of loci under positive selection in the transitional group revealed genes involved in host–pathogen interactions that may facilitate host switching. Finally, the first B. clarridgeiae genome from Asia exhibited large-scale genome stability compared to isolates from other continents, except for SNPs in regions predicted to mediate interactions with the vertebrate host. These findings highlight the paucity of data on the genomic diversity of Ctenocephalides -associated bacteria and raise questions regarding how interactions between members of the flea microbiome might influence vector competence.

The genetic network underlying the evolution of pathogenicity in avian Escherichia coli

Introduction

Colibacillosis is a worldwide prevalent disease in poultry production linked to Escherichia coli strains that belong to the avian pathogenic E. coli (APEC) pathotype. While many virulence factors have been linked to APEC isolates, no single gene or set of genes has been found to be exclusively associated with the pathotype. Moreover, a comprehensive description of the biological processes linked to APEC pathogenicity is currently lacking.

Methods

In this study, we compiled a dataset of 2015 high-quality avian E. coli genomes from pathogenic and commensal isolates, based on publications from 2000 to 2021. We then conducted a genome-wide association study (GWAS) and integrated candidate gene identification with available protein-protein interaction data to decipher the genetic network underlying the biological processes connected to APEC pathogenicity.

Results

Our GWAS identified variations in gene content for 13 genes and SNPs in 3 different genes associated with APEC isolates, suggesting both gene-level and SNP-level variations contribute to APEC pathogenicity. Integrating protein-protein interaction data, we found that 15 of these genes clustered in the same genetic network, suggesting the pathogenicity of APEC might be due to the interplay of different regulated pathways. We also found novel candidate genes including an uncharacterized multi-pass membrane protein (yciC) and the outer membrane porin (ompD) as linked to APEC isolates.

Discussion

Our findings suggest that convergent pathways related to nutrient uptake from host cells and defense from host immune system play a major role in APEC pathogenicity. In addition, the dataset curated in this study represents a comprehensive historical genomic collection of avian E. coli isolates and constitutes a valuable resource for their comparative genomics investigations.

Stenotrophomonas maltophilia natural history and evolution in the airways of adults with cystic fibrosis

Introduction

Stenotrophomonas maltophilia is an opportunistic pathogen infecting persons with cystic fibrosis (pwCF) and portends a worse prognosis. Studies of S. maltophilia infection dynamics have been limited by cohort size and follow-up. We investigated the natural history, transmission potential, and evolution of S. maltophilia in a large Canadian cohort of 321 pwCF over a 37-year period.

Methods

One-hundred sixty-two isolates from 74 pwCF (23%) were typed by pulsed-field gel electrophoresis, and shared pulsotypes underwent whole-genome sequencing.

Results

S. maltophilia was recovered at least once in 82 pwCF (25.5%). Sixty-four pwCF were infected by unique pulsotypes, but shared pulsotypes were observed between 10 pwCF. In chronic carriage, longer time periods between positive sputum cultures increased the likelihood that subsequent isolates were unrelated. Isolates from individual pwCF were largely clonal, with differences in gene content being the primary source of genetic diversity objectified by gene content differences. Disproportionate progression of CF lung disease was not observed amongst those infected with multiple strains over time (versus a single) or amongst those with shared clones (versus strains only infecting one patient). We did not observe evidence of patient-to-patient transmission despite relatedness between isolates. Twenty-four genes with ≥ 2 mutations accumulated over time were identified across 42 sequenced isolates from all 11 pwCF with ≥ 2 sequenced isolates, suggesting a potential role for these genes in adaptation of S. maltophilia to the CF lung.

Discussion

Genomic analyses suggested common, indirect sources as the origins of S. maltophilia infections in the clinic population. The information derived from a genomics-based understanding of the natural history of S. maltophilia infection within CF provides unique insight into its potential for in-host evolution.