Single-Nucleotide Polymorphism-Based Genetic Diversity Analysis of Clinical Pseudomonas aeruginosa Isolates

The core genome evolution of Lactobacillus crispatus as a driving force for niche competition in the human vaginal tract

The lower female reproductive tract is notoriously dominated by Lactobacillus species, among which Lactobacillus crispatus emerges for its protective and health-promoting activities. Although previous comparative genome analyses highlighted genetic and phenotypic diversity within the L. crispatus species, most studies have focused on the presence/absence of accessory genes. Here, we investigated the variation at the single nucleotide level within protein-encoding genes shared across a human-derived L. crispatus strain selection, which includes 200 currently available human-derived L. crispatus genomes as well as 41 chromosome sequences of such taxon that have been decoded in the framework of this study. Such data clearly pointed out the presence of intra-species micro-diversities that could have evolutionary significance contributing to phenotypical diversification by affecting protein domains. Specifically, two single nucleotide variations in the type II pullulanase gene sequence led to specific amino acid substitutions, possibly explaining the substantial differences in the growth performances and competition abilities observed in a multi-strain bioreactor culture simulating the vaginal environment. Accordingly, L. crispatus strains display different growth performances, suggesting that the colonisation and stable persistence in the female reproductive tract between the members of this taxon is highly variable.

Comparative Genomics of Pseudomonas aeruginosa Strains Isolated from Different Ecological Niches

The Pseudomonas aeruginosa genome can change to adapt to different ecological niches. We compared four genomes from a Mexican hospital and 59 genomes from GenBank from different niches, such as urine, sputum, and environmental. The ST analysis showed that high-risk STs (ST235, ST773, and ST27) were present in the genomes of the three niches from GenBank, and the STs of Mexican genomes (ST167, ST2731, and ST549) differed from the GenBank genomes. Phylogenetic analysis showed that the genomes were clustering according to their ST and not their niche. When analyzing the genomic content, we observed that environmental genomes had genes involved in adapting to the environment not found in the clinics and that their mechanisms of resistance were mutations in antibiotic resistance-related genes. In contrast, clinical genomes from GenBank had resistance genes, in mobile/mobilizable genetic elements in the chromosome, except for the Mexican genomes that carried them mostly in plasmids. This was related to the presence of CRISPR-Cas and anti-CRISPR; however, Mexican strains only had plasmids and CRISPR-Cas. bla_OXA-488 (a variant of bla_OXA50) with higher activity against carbapenems was more prevalent in sputum genomes. The virulome analysis showed that exoS was most prevalent in the genomes of urinary samples and exoU and pldA in sputum samples. This study provides evidence regarding the genetic variability among P. aeruginosa isolated from different niches.

Computationally Efficient Assembly of Pseudomonas aeruginosa Gene Expression Compendia

Thousands of Pseudomonas aeruginosa RNA sequencing (RNA-seq) gene expression profiles are publicly available via the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA). In this work, the transcriptional profiles from hundreds of studies performed by over 75 research groups were reanalyzed in aggregate to create a powerful tool for hypothesis generation and testing. Raw sequence data were uniformly processed using the Salmon pseudoaligner, and this read mapping method was validated by comparison to a direct alignment method. We developed filtering criteria to exclude samples with aberrant levels of housekeeping gene expression or an unexpected number of genes with no reported values and normalized the filtered compendia using the ratio-of-medians method. The filtering and normalization steps greatly improved gene expression correlations for genes within the same operon or regulon across the 2,333 samples. Since the RNA-seq data were generated using diverse strains, we report the effects of mapping samples to noncognate reference genomes by separately analyzing all samples mapped to cDNA reference genomes for strains PAO1 and PA14, two divergent strains that were used to generate most of the samples. Finally, we developed an algorithm to incorporate new data as they are deposited into the SRA. Our processing and quality control methods provide a scalable framework for taking advantage of the troves of biological information hibernating in the depths of microbial gene expression data and yield useful tools for P. aeruginosa RNA-seq data to be leveraged for diverse research goals. IMPORTANCE Pseudomonas aeruginosa is a causative agent of a wide range of infections, including chronic infections associated with cystic fibrosis. These P. aeruginosa infections are difficult to treat and often have negative outcomes. To aid in the study of this problematic pathogen, we mapped, filtered for quality, and normalized thousands of P. aeruginosa RNA-seq gene expression profiles that were publicly available via the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA). The resulting compendia facilitate analyses across experiments, strains, and conditions. Ultimately, the workflow that we present could be applied to analyses of other microbial species.

Comparative Genomics of Streptococcus oralis Identifies Large Scale Homologous Recombination and a Genetic Variant Associated with Infection

The viridans group streptococci (VGS) are a large consortium of commensal streptococci that colonize the human body. Many species within this group are opportunistic pathogens causing bacteremia and infective endocarditis (IE), yet little is known about why some strains cause invasive disease. Identification of virulence determinants is complicated by the difficulty of distinguishing between the closely related species of this group. Here, we analyzed genomic data from VGS that were isolated from blood cultures in patients with invasive infections and oral swabs of healthy volunteers and then determined the best-performing methods for species identification. Using whole-genome sequence data, we characterized the population structure of a diverse sample of Streptococcus oralis isolates and found evidence of frequent recombination. We used multiple genome-wide association study tools to identify candidate determinants of invasiveness. These tools gave consistent results, leading to the discovery of a single synonymous single nucleotide polymorphism (SNP) that was significantly associated with invasiveness. This SNP was within a previously undescribed gene that was conserved across the majority of VGS species. Using the growth in the presence of human serum and a simulated infective endocarditis vegetation model, we were unable to identify a phenotype for the enriched allele in laboratory assays, suggesting a phenotype may be specific to natural infection. These data highlighted the power of analyzing natural populations for gaining insight into pathogenicity, particularly for organisms with complex population structures like the VGS. IMPORTANCE The viridians group streptococci (VGS) are a large collection of closely related commensal streptococci, with many being opportunistic pathogens causing invasive diseases, such as bacteremia and infective endocarditis. Little is known about virulence determinants in these species, and there is a distinct lack of genomic information available for the VGS. In this study, we collected VGS isolates from invasive infections and healthy volunteers and performed whole-genome sequencing for a suite of downstream analyses. We focused on a diverse sample of Streptococcus oralis genomes and identified high rates of recombination in the population as well as a single genome variant highly enriched in invasive isolates. The variant lies within a previously uncharacterized gene, nrdM, which shared homology with the anaerobic ribonucleoside triphosphate reductase, nrdD, and was highly conserved among VGS. This work increased our knowledge of VGS genomics and indicated that differences in virulence potential among S. oralis isolates were, at least in part, genetically determined.

The Pseudomonas aeruginosa whole genome sequence: A 20th anniversary celebration

Toward the end of August 2000, the 6.3 Mbp whole genome sequence of Pseudomonas aeruginosa strain PAO1 was published. With 5570 open reading frames (ORFs), PAO1 had the largest microbial genome sequenced up to that point in time-including a large proportion of metabolic, transport and antimicrobial resistance genes supporting its ability to colonize diverse environments. A remarkable 9% of its ORFs were predicted to encode proteins with regulatory functions, providing new insight into bacterial network complexity as a function of network size. In this celebratory article, we fast forward 20 years, and examine how access to this resource has transformed our understanding of P. aeruginosa. What follows is more than a simple review or commentary; we have specifically asked some of the leaders in the field to provide personal reflections on how the PAO1 genome sequence, along with the Pseudomonas Community Annotation Project (PseudoCAP) and Pseudomonas Genome Database (pseudomonas.com), have contributed to the many exciting discoveries in this field. In addition to bringing us all up to date with the latest developments, we also ask our contributors to speculate on how the next 20 years of Pseudomonas research might pan out.

Pseudomonas aeruginosa: a clinical and genomics update

Antimicrobial resistance (AMR) has become a global medical priority that needs urgent resolution. Pseudomonas aeruginosa is a versatile, adaptable bacterial species with widespread environmental occurrence, strong medical relevance, a diverse set of virulence genes and a multitude of intrinsic and possibly acquired antibiotic resistance traits. P. aeruginosa causes a wide variety of infections and has an epidemic-clonal population structure. Several of its dominant global clones have collected a wide variety of resistance genes rendering them multi-drug resistant (MDR) and particularly threatening groups of vulnerable individuals including surgical patients, immunocompromised patients, Caucasians suffering from cystic fibrosis (CF) and more. AMR and MDR especially are particularly problematic in P. aeruginosa significantly complicating successful antibiotic treatment. In addition, antimicrobial susceptibility testing (AST) of P. aeruginosa can be cumbersome due to its slow growth or the massive production of exopolysaccharides and other extracellular compounds. For that reason, phenotypic AST is progressively challenged by genotypic methods using whole genome sequences (WGS) and large-scale phenotype databases as a framework of reference. We here summarize the state of affairs and the quality level of WGS-based AST for P. aeruginosa mostly from clinical origin.

Phenotypic and Genomic Comparison of the Two Most Common ExoU-Positive Pseudomonas aeruginosa Clones, PA14 and ST235

Genotyping of 2,882 Pseudomonas aeruginosa isolates that had been collected during the last 40 years identified the ExoU-positive lineages PA14 (ST253) and ST235 as the second and third most frequent clones in the P. aeruginosa population. Both clones were approximately 2-fold more frequently detected in animate habitats than in soil or aquatic habitats. While ST253 clone isolates were causing mainly acute and chronic infections in humans, ST235 isolates had been preferentially collected from hospitalized patients with severe acute infections, particularly, keratitis, urinary tract infections, burn wounds, and ventilator-associated pneumonia. The two major exoU clones differed substantially in the composition and flexibility of the accessory genome and by more than 8,000 amino acid sequences. Pronounced sequence variation between orthologs was noted in genes encoding elements of secretion systems and secreted effector molecules, including the type III secretion system, indicating the modes of action of the different clones. When comparing representatives of the two clones in batch culture, the PA14 strain orchestrated the quorum sensing circuitry for the expression of pathogenic traits and stopped growing in batch culture when it entered the stationary phase, but the quorum sensing-deficient ST235 strain expressed high type III secretion activity and continued to grow and to divide. In summary, unrestricted growth, high constitutive type III secretion activity, and facilitated uptake of foreign DNA could be major features that have made ST235 a global high-risk clone associated with poor outcomes of acute nosocomial infections.IMPORTANCE The ubiquitous and metabolically versatile environmental bacterium Pseudomonas aeruginosa can cause infections in a wide variety of hosts, including insects, plants, animals, and humans. P. aeruginosa is one of the ESKAPE (Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa, and Enterobacter species) pathogens that are the major cause of nosocomial infections in the United States and are a threat all over the world because of their capacity to become increasingly resistant to all available antibiotics. Most experimental work on P. aeruginosa has been performed with reference strains PAO1 and PA14, providing deep insight into key metabolic and regulatory pathways thought to be applicable to all P. aeruginosa strains. However, this comparative study on the two most common exoU-positive clones taught us that there are major lineages in the population such as the global high-risk clone ST235 that exhibit uncommon traits of lifestyle, genome mobility, and pathogenicity distinct from those in our knowledge gained from the studies with the reference strains.