BOG: R-package for Bacterium and virus analysis of Orthologous Groups

Viromes vs. mixed community metagenomes: choice of method dictates interpretation of viral community ecology

Background

Viruses, the majority of which are uncultivated, are among the most abundant biological entities on Earth. From altering microbial physiology to driving community dynamics, viruses are fundamental members of microbiomes. While the number of studies leveraging viral metagenomics (viromics) for studying uncultivated viruses is growing, standards for viromics research are lacking. Viromics can utilize computational discovery of viruses from total metagenomes of all community members (hereafter metagenomes) or use physical separation of virus-specific fractions (hereafter viromes). However, differences in the recovery and interpretation of viruses from metagenomes and viromes obtained from the same samples remain understudied.

Results

Here, we compare viral communities from paired viromes and metagenomes obtained from 60 diverse samples across human gut, soil, freshwater, and marine ecosystems. Overall, viral communities obtained from viromes were more abundant and species rich than those obtained from metagenomes, although there were some exceptions. Despite this, metagenomes still contained many viral genomes not detected in viromes. We also found notable differences in the predicted lytic state of viruses detected in viromes vs metagenomes at the time of sequencing. Other forms of variation observed include genome presence/absence, genome quality, and encoded protein content between viromes and metagenomes, but the magnitude of these differences varied by environment.

Conclusions

Overall, our results show that the choice of method can lead to differing interpretations of viral community ecology. We suggest that the choice of whether to target a metagenome or virome to study viral communities should be dependent on the environmental context and ecological questions being asked. However, our overall recommendation to researchers investigating viral ecology and evolution is to pair both approaches to maximize their respective benefits.

Large-scale comparative genomics to refine the organization of the global Salmonella enterica population structure

The White-Kauffmann-Le Minor (WKL) scheme is the most widely used Salmonella typing scheme for reporting the disease prevalence of the enteric pathogen. With the advent of whole-genome sequencing (WGS), in silico methods have increasingly replaced traditional serotyping due to reproducibility, speed and coverage. However, despite integrating genomic-based typing by in silico serotyping tools such as SISTR, in silico serotyping in certain contexts remains ambiguous and insufficiently informative. Specifically, in silico serotyping does not attempt to resolve polyphyly. Furthermore, in spite of the widespread acknowledgement of polyphyly from genomic studies, the prevalence of polyphyletic serovars is not well characterized. Here, we applied a genomics approach to acquire the necessary resolution to classify genetically discordant serovars and propose an alternative typing scheme that consistently reflect natural Salmonella populations. By accessing the unprecedented volume of bacterial genomic data publicly available in GenomeTrakr and PubMLST databases (>180 000 genomes representing 723 serovars), we characterized the global Salmonella population structure and systematically identified putative non-monophyletic serovars. The proportion of putative non-monophyletic serovars was estimated higher than previous reports, reinforcing the inability of antigenic determinants to depict the complexity of Salmonella evolutionary history. We explored the extent of genetic diversity masked by serotyping labels and found significant intra-serovar molecular differences across many clinically important serovars. To avoid false discovery due to incorrect in silico serotyping calls, we cross-referenced reported serovar labels and concluded a low error rate in in silico serotyping. The combined application of clustering statistics and genome-wide association methods demonstrated effective characterization of stable bacterial populations and explained functional differences. The collective methods adopted in our study have practical values in establishing genomic-based typing nomenclatures for an entire microbial species or closely related subpopulations. Ultimately, we foresee an improved typing scheme to be a hybrid that integrates both genomic and antigenic information such that the resolution from WGS is leveraged to improve the precision of subpopulation classification while preserving the common names defined by the WKL scheme.

Proteomic response of Turicibacter bilis MMM721 to chicken bile and its bile acids

Objective

Bile and its individual components, mainly bile acids, are important for digestion and drive bacterial community dynamics in the upper gastrointestinal tract of chickens. However, specific responses to bile acids have been characterized in only a few commensal bacteria, and it is unclear how other members of the microbiota respond to biliary stress. Here, we used label-free LC–MS/MS to assess the proteomic response of a common inhabitant of the chicken small intestine, Turicibacter bilis MMM721, to 24 h of growth in anaerobic growth media supplemented with 0.1% whole chicken bile, 0.1% taurochenodeoxycholic acid (TCDCA), or 0.1% taurocholic acid (TCA).

Results

Seventy, 46, and 10 differentially expressed proteins were identified in Turicibacter bilis MMM721 cultured with supplements of chicken bile, TCDCA, and TCA, respectively, when compared to unsupplemented controls. Many differentially expressed proteins were predicted to be involved in ribosomal processes, post-translational modifications and chaperones, and modifications to the cell surface. Ultimately, the T. bilis MMM721 response to whole bile and bile acids is complex and may relate to adaptations for small intestine colonization, with numerous proteins from a variety of functional categories being impacted.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13104-022-06127-8.

Osmotic stress induces long-term biofilm survival in Liberibacter crescens

Citrus greening, also known as Huanglongbing (HLB), is a devastating citrus plant disease caused predominantly by Liberibacter asiaticus. While nearly all Liberibacter species remain uncultured, here we used the culturable L. crescens BT-1 as a model to examine physiological changes in response to the variable osmotic conditions and nutrient availability encountered within the citrus host. Similarly, physiological responses to changes in growth temperature and dimethyl sulfoxide concentrations were also examined, due to their use in many of the currently employed therapies to control the spread of HLB. Sublethal heat stress was found to induce the expression of genes related to tryptophan biosynthesis, while repressing the expression of ribosomal proteins. Osmotic stress induces expression of transcriptional regulators involved in expression of extracellular structures, while repressing the biosynthesis of fatty acids and aromatic amino acids. The effects of osmotic stress were further evaluated by quantifying biofilm formation of L. crescens in presence of increasing sucrose concentrations at different stages of biofilm formation, where sucrose-induced osmotic stress delayed initial cell attachment while enhancing long-term biofilm viability. Our findings revealed that exposure to osmotic stress is a significant contributing factor to the long term survival of L. crescens and, possibly, to the pathogenicity of other Liberibacter species.

Genetic and metabolic signatures of Salmonella enterica subsp. enterica associated with animal sources at the pangenomic scale

BACKGROUND

Salmonella enterica subsp. enterica is a public health issue related to food safety, and its adaptation to animal sources remains poorly described at the pangenome scale. Firstly, serovars presenting potential mono- and multi-animal sources were selected from a curated and synthetized subset of Enterobase. The corresponding sequencing reads were downloaded from the European Nucleotide Archive (ENA) providing a balanced dataset of 440 Salmonella genomes in terms of serovars and sources (i). Secondly, the coregenome variants and accessory genes were detected (ii). Thirdly, single nucleotide polymorphisms and small insertions/deletions from the coregenome, as well as the accessory genes were associated to animal sources based on a microbial Genome Wide Association Study (GWAS) integrating an advanced correction of the population structure (iii). Lastly, a Gene Ontology Enrichment Analysis (GOEA) was applied to emphasize metabolic pathways mainly impacted by the pangenomic mutations associated to animal sources (iv).

RESULTS

Based on a genome dataset including Salmonella serovars from mono- and multi-animal sources (i), 19,130 accessory genes and 178,351 coregenome variants were identified (ii). Among these pangenomic mutations, 52 genomic signatures (iii) and 9 over-enriched metabolic signatures (iv) were associated to avian, bovine, swine and fish sources by GWAS and GOEA, respectively.

CONCLUSIONS

Our results suggest that the genetic and metabolic determinants of Salmonella adaptation to animal sources may have been driven by the natural feeding environment of the animal, distinct livestock diets modified by human, environmental stimuli, physiological properties of the animal itself, and work habits for health protection of livestock.

LdtR is a master regulator of gene expression in Liberibacter asiaticus

Huanglongbing or citrus greening disease is causing devastation to the citrus industry. Liberibacter asiaticus, an obligate intracellular pathogen of citrus, is one the causative agents of the disease. Most of the knowledge about this bacterium has been deduced from the in silico exploration of its genomic sequence. L. asiaticus differentially expresses genes during its transmission from the psyllid vector, Diaphorina citri, to the plant. However, the regulatory mechanisms for the adaptation of the bacterium into either hosts remain unknown. Here we show that LdtR, a MarR family transcriptional regulator, activates or represses transcription genome-wide. We performed a double approach to identify the components of the LdtR regulon: a transcriptome analysis in both the related bacterium Liberibacter crescens and citrus-infected leaves, strengthened with an in silico prediction of LdtR regulatory sites. Our results demonstrated that LdtR controls the expression of nearly 180 genes in L. asiaticus, distributed in processes such as cell motility, cell wall biogenesis, energy production, and transcription. These results provide new evidence about the regulatory network of L. asiaticus, where the differential expression of genes from these functional categories could be of great importance during the adaptation of the bacterium to either hosts.

Global Analysis and Comparison of the Transcriptomes and Proteomes of Group A Streptococcus Biofilms

Prokaryotes are thought to regulate their proteomes largely at the level of transcription. However, the results from this first set of global transcriptomic and proteomic analyses of paired microbial samples presented here show that this assumption is false for the majority of genes and their products in S. pyogenes. In addition, the tenuousness of the link between transcription and translation becomes even more pronounced when microbes exist in a biofilm or a stationary planktonic state. Since the transcriptome level does not usually equal the proteome level, the validity attributed to gene expression studies as well as proteomic studies in microbial analyses must be brought into question. Therefore, the results attained by either approach, whether RNA-seq or shotgun proteomics, must be taken in context and evaluated with particular care since they are by no means interchangeable.

To gain a better understanding of the genes and proteins involved in group A Streptococcus (GAS; Streptococcus pyogenes) biofilm growth, we analyzed the transcriptome, cellular proteome, and cell wall proteome from biofilms at different stages and compared them to those of plankton-stage GAS. Using high-throughput RNA sequencing (RNA-seq) and liquid chromatography-tandem mass spectrometry (LC-MS/MS) shotgun proteomics, we found distinct expression profiles in the transcriptome and proteome. A total of 46 genes and 41 proteins showed expression across the majority of biofilm time points that was consistently higher or consistently lower than that seen across the majority of planktonic time points. However, there was little overlap between the genes and proteins on these two lists. In line with other studies comparing transcriptomic and proteomic data, the overall correlation between the two data sets was modest. Furthermore, correlation was poorest for biofilm samples. This suggests a high degree of regulation of protein expression by nontranscriptional mechanisms. This report illustrates the benefits and weaknesses of two different approaches to global expression profiling, and it also demonstrates the advantage of using proteomics in conjunction with transcriptomics to gain a more complete picture of global expression within biofilms. In addition, this report provides the fullest characterization of expression patterns in GAS biofilms currently available.

IMPORTANCE Prokaryotes are thought to regulate their proteomes largely at the level of transcription. However, the results from this first set of global transcriptomic and proteomic analyses of paired microbial samples presented here show that this assumption is false for the majority of genes and their products in S. pyogenes. In addition, the tenuousness of the link between transcription and translation becomes even more pronounced when microbes exist in a biofilm or a stationary planktonic state. Since the transcriptome level does not usually equal the proteome level, the validity attributed to gene expression studies as well as proteomic studies in microbial analyses must be brought into question. Therefore, the results attained by either approach, whether RNA-seq or shotgun proteomics, must be taken in context and evaluated with particular care since they are by no means interchangeable.