Compendium-Wide Analysis of Pseudomonas aeruginosa Core and Accessory Genes Reveals Transcriptional Patterns across Strains PAO1 and PA14

Integration of 168,000 samples reveals global patterns of the human gut microbiome

The factors shaping human microbiome variation are a major focus of biomedical research. While other fields have used large sequencing compendia to extract insights requiring otherwise impractical sample sizes, the microbiome field has lacked a comparably sized resource for the 16S rRNA gene amplicon sequencing commonly used to quantify microbiome composition. To address this gap, we processed 168,464 publicly available human gut microbiome samples with a uniform pipeline. We use this compendium to evaluate geographic and technical effects on microbiome variation. We find that regions such as Central and Southern Asia differ significantly from the more thoroughly characterized microbiomes of Europe and Northern America and that composition alone can be used to predict a sample's region of origin. We also find strong associations between microbiome variation and technical factors such as primers and DNA extraction. We anticipate this growing work, the Human Microbiome Compendium, will enable advanced applied and methodological research.

The diadenosine tetraphosphate hydrolase ApaH contributes to Pseudomonas aeruginosa pathogenicity

The opportunistic bacterial pathogen Pseudomonas aeruginosa causes a wide range of infections that are difficult to treat, largely because of the spread of antibiotic-resistant isolates. Antivirulence therapy, í.e. the use of drugs that inhibit the expression or activity of virulence factors, is currently considered an attractive strategy to reduce P. aeruginosa pathogenicity and complement antibiotic treatments. Because of the multifactorial nature of P. aeruginosa virulence and the broad arsenal of virulence factors this bacterium can produce, the regulatory networks that control the expression of multiple virulence traits have been extensively explored as potential targets for antivirulence drug development. The intracellular signaling molecule diadenosine tetraphosphate (Ap4A) has been reported to control stress resistance and virulence-related traits in some bacteria, but its role has not been investigated in P. aeruginosa so far. To fill this gap, we generated a mutant of the reference strain P. aeruginosa PAO1 that lacks the Ap4A-hydrolysing enzyme ApaH and, consequently, accumulates high intracellular levels of Ap4A. Phenotypic and transcriptomic analyses revealed that the lack of ApaH causes a drastic reduction in the expression of several virulence factors, including extracellular proteases, elastases, siderophores, and quorum sensing signal molecules. Accordingly, infection assays in plant and animal models demonstrated that ApaH-deficient cells are significantly impaired in infectivity and persistence in different hosts, including mice. Finally, deletion of apaH in P. aeruginosa clinical isolates demonstrated that the positive effect of ApaH on the production of virulence-related traits and on infectivity is conserved in P. aeruginosa. This study provides the first evidence that the Ap4A-hydrolysing enzyme ApaH is important for P. aeruginosa virulence, highlighting this protein as a novel potential target for antivirulence therapies against P. aeruginosa.

Widespread fungal-bacterial competition for magnesium lowers bacterial susceptibility to polymyxin antibiotics

Fungi and bacteria coexist in many polymicrobial communities, yet the molecular basis of their interactions remains poorly understood. Here, we show that the fungus Candida albicans sequesters essential magnesium ions from the bacterium Pseudomonas aeruginosa. To counteract fungal Mg2+ sequestration, P. aeruginosa expresses the Mg2+ transporter MgtA when Mg2+ levels are low. Thus, loss of MgtA specifically impairs P. aeruginosa in co-culture with C. albicans, but fitness can be restored by supplementing Mg2+. Using a panel of fungi and bacteria, we show that Mg2+ sequestration is a general mechanism of fungal antagonism against gram-negative bacteria. Mg2+ limitation enhances bacterial resistance to polymyxin antibiotics like colistin, which target gram-negative bacterial membranes. Indeed, experimental evolution reveals that P. aeruginosa evolves C. albicans-dependent colistin resistance via non-canonical means; antifungal treatment renders resistant bacteria colistin-sensitive. Our work suggests that fungal-bacterial competition could profoundly impact polymicrobial infection treatment with antibiotics of last resort.

Pseudomonas aeruginosa transcriptome analysis of metal restriction in ex vivo cystic fibrosis sputum

Chronic Pseudomonas aeruginosa lung infections are a feature of cystic fibrosis (CF) that many patients experience even with the advent of highly effective modulator therapies. Identifying factors that impact P. aeruginosa in the CF lung could yield novel strategies to eradicate infection or otherwise improve outcomes. To complement published P. aeruginosa studies using laboratory models or RNA isolated from sputum, we analyzed transcripts of strain PAO1 after incubation in sputum from different CF donors prior to RNA extraction. We compared PAO1 gene expression in this "spike-in" sputum model to that for P. aeruginosa grown in synthetic cystic fibrosis sputum medium to determine key genes, which are among the most differentially expressed or most highly expressed. Using the key genes, gene sets with correlated expression were determined using the gene expression analysis tool eADAGE. Gene sets were used to analyze the activity of specific pathways in P. aeruginosa grown in sputum from different individuals. Gene sets that we found to be more active in sputum showed similar activation in published data that included P. aeruginosa RNA isolated from sputum relative to corresponding in vitro reference cultures. In the ex vivo samples, P. aeruginosa had increased levels of genes related to zinc and iron acquisition which were suppressed by metal amendment of sputum. We also found a significant correlation between expression of the H1-type VI secretion system and CFTR corrector use by the sputum donor. An ex vivo sputum model or synthetic sputum medium formulation that imposes metal restriction may enhance future CF-related studies.IMPORTANCEIdentifying the gene expression programs used by Pseudomonas aeruginosa to colonize the lungs of people with cystic fibrosis (CF) will illuminate new therapeutic strategies. To capture these transcriptional programs, we cultured the common P. aeruginosa laboratory strain PAO1 in expectorated sputum from CF patient donors. Through bioinformatic analysis, we defined sets of genes that are more transcriptionally active in real CF sputum compared to a synthetic cystic fibrosis sputum medium. Many of the most differentially active gene sets contained genes related to metal acquisition, suggesting that these gene sets play an active role in scavenging for metals in the CF lung environment which may be inadequately represented in some models. Future studies of P. aeruginosa transcript abundance in CF may benefit from the use of an expectorated sputum model or media supplemented with factors that induce metal restriction.

Comparative whole-genome analysis of China and global epidemic Pseudomonas aeruginosa high-risk clones

OBJECTIVES

The various sequence types (STs) of Pseudomonas aeruginosa (P. aeruginosa) high-risk clones (HiRiCs) have been sporadically reported in China, but the systematic analysis of genomes for these STs remains limited. This study aimed to address the evolutionary pathways underlying the emergence of HiRiCs and their routes of dissemination from Chinese and global perspectives.

METHODS

The phylogenetic analysis was performed based on 416 newly sequenced clinical P. aeruginosa strains from Guangdong (GD), published genome sequences of 282 Chinese isolates, and 868 HiRiCs isolates from other countries. The genomic comparison study of global HiRiC ST244 was conducted to detect the model of global dissemination and local separation driven by association regional-specific antibiotic resistance genes. Furthermore, the evolutionary route of the emerging, China-specific HiRiC ST1971 was explored using Most Recent Common Ancestor (MRCA) analysis.

RESULTS

Based on comparative genomics analysis, we found a clear geographical separation of ST244 isolates, yet with an association between ST244 isolates from GD and America. We identified a set of 38 AMR genes that contribute to the geographical separation in ST244, and we identified genetic determinants either positively (MexB) and negatively (opmD) associated with GD ST244. For the China-unique HiRiC ST1971, its evolutionary history across different continents before emerging as ST1971 in China was also deduced.

CONCLUSION

This study provides insight into the specific genetics underlying regional differences among globally disseminated P. aeruginosa HiRiCs (ST244) as well as new understanding of the dissemination and evolution of a regional HiRiC (ST1971). Understanding the genetics of these and other HiRiCs may assist in controlling their emergence and further spread.

Integration of 168,000 samples reveals global patterns of the human gut microbiome

Understanding the factors that shape variation in the human microbiome is a major goal of research in biology. While other genomics fields have used large, pre-compiled compendia to extract systematic insights requiring otherwise impractical sample sizes, there has been no comparable resource for the 16S rRNA sequencing data commonly used to quantify microbiome composition. To help close this gap, we have assembled a set of 168,484 publicly available human gut microbiome samples, processed with a single pipeline and combined into the largest unified microbiome dataset to date. We use this resource, which is freely available at microbiomap.org, to shed light on global variation in the human gut microbiome. We find that Firmicutes, particularly Bacilli and Clostridia, are almost universally present in the human gut. At the same time, the relative abundance of the 65 most common microbial genera differ between at least two world regions. We also show that gut microbiomes in undersampled world regions, such as Central and Southern Asia, differ significantly from the more thoroughly characterized microbiomes of Europe and Northern America. Moreover, humans in these overlooked regions likely harbor hundreds of taxa that have not yet been discovered due to this undersampling, highlighting the need for diversity in microbiome studies. We anticipate that this new compendium can serve the community and enable advanced applied and methodological research.

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.