Genetic and Functional Analyses of Cutibacterium Acnes Isolates Reveal the Association of a Linear Plasmid with Skin Inflammation.

Cutibacterium acnes is a commensal bacterium on the skin that is generally well-tolerated, but different strain types have been hypothesized to contribute to the disease acne vulgaris. To understand how some strain types might contribute to skin inflammation, we generated a repository of C. acnes isolates from skin swabs of healthy subjects and subjects with acne and assessed their strain-level identity and capacity to stimulate cytokine release. Phylotype II K-type strains were more frequent on healthy and nonlesional skin of subjects with acne than those isolated from lesions. Phylotype IA-1 C-type strains were increased on lesional skin compared with those on healthy skin. The capacity to induce cytokines from cultured monocyte-derived dendritic cells was opposite to this action on sebocytes and keratinocytes and did not correlate with the strain types associated with the disease. Whole-genome sequencing revealed a linear plasmid in high-inflammatory isolates within similar strain types that had different proinflammatory responses. Single-cell RNA sequencing of mouse skin after intradermal injection showed that strains containing this plasmid induced a higher inflammatory response in dermal fibroblasts. These findings revealed that C. acnes strain type is insufficient to predict inflammation and that carriage of a plasmid could contribute to disease.

Guild and Niche Determination Enable Targeted Alteration of the Microbiome.

Microbiome science has greatly contributed to our understanding of microbial life and its essential roles for the environment and human health1-5. However, the nature of microbial interactions and how microbial communities respond to perturbations remains poorly understood, resulting in an often descriptive and correlation-based approach to microbiome research6-8. Achieving causal and predictive microbiome science would require direct functional measurements in complex communities to better understand the metabolic role of each member and its interactions with others. In this study we present a new approach that integrates transcription and translation measurements to predict competition and substrate preferences within microbial communities, consequently enabling the selective manipulation of the microbiome. By performing metatranscriptomic (metaRNA-Seq) and metatranslatomic (metaRibo-Seq) analysis in complex samples, we classified microbes into functional groups (i.e. guilds) and demonstrated that members of the same guild are competitors. Furthermore, we predicted preferred substrates based on importer proteins, which specifically benefited selected microbes in the community (i.e. their niche) and simultaneously impaired their competitors. We demonstrated the scalability of microbial guild and niche determination to natural samples and its ability to successfully manipulate microorganisms in complex microbiomes. Thus, the approach enhances the design of pre- and probiotic interventions to selectively alter members within microbial communities, advances our understanding of microbial interactions, and paves the way for establishing causality in microbiome science.

A Genome-Scale Atlas Reveals Complex Interplay of Transcription and Translation in an Archaeon.

The scale of post-transcriptional regulation and the implications of its interplay with other forms of regulation in environmental acclimation are underexplored for organisms of the domain Archaea. Here, we have investigated the scale of post-transcriptional regulation in the extremely halophilic archaeon Halobacterium salinarum NRC-1 by integrating the transcriptome-wide locations of transcript processing sites (TPSs) and SmAP1 binding, the genome-wide locations of antisense RNAs (asRNAs), and the consequences of RNase_2099C knockout on the differential expression of all genes. This integrated analysis has discovered that 54% of all protein-coding genes in the genome of this haloarchaeon are likely targeted by multiple mechanisms for putative post-transcriptional processing and regulation, with about 20% of genes likely being regulated by combinatorial schemes involving SmAP1, asRNAs, and RNase_2099C. Comparative analysis of mRNA levels (transcriptome sequencing [RNA-Seq]) and protein levels (sequential window acquisition of all theoretical fragment ion spectra mass spectrometry [SWATH-MS]) for 2,579 genes over four phases of batch culture growth in complex medium generated additional evidence for the conditional post-transcriptional regulation of 7% of all protein-coding genes. We demonstrate that post-transcriptional regulation may act to fine-tune specialized and rapid acclimation to stressful environments, e.g., as a switch to turn on gas vesicle biogenesis to promote vertical relocation under anoxic conditions and modulate the frequency of transposition by insertion sequence (IS) elements of the IS200/IS605, IS4, and ISH3 families. Findings from this study are provided as an atlas in a public Web resource (https://halodata.systemsbiology.net). IMPORTANCE While the transcriptional regulation landscape of archaea has been extensively investigated, we currently have limited knowledge about post-transcriptional regulation and its driving mechanisms in this domain of life. In this study, we collected and integrated omics data from multiple sources and technologies to infer post-transcriptionally regulated genes and the putative mechanisms modulating their expression at the protein level in Halobacterium salinarum NRC-1. The results suggest that post-transcriptional regulation may drive environmental acclimation by regulating hallmark biological processes. To foster discoveries by other research groups interested in the topic, we extended our integrated data to the public in the form of an interactive atlas (https://halodata.systemsbiology.net).

synDNA-a Synthetic DNA Spike-in Method for Absolute Quantification of Shotgun Metagenomic Sequencing.

Microbiome studies have the common goal of determining which microbial taxa are present, respond to specific conditions, or promote phenotypic changes in the host. Most of these studies rely on relative abundance measurements to drive conclusions. Inherent limitations of relative values are the inability to determine whether an individual taxon is more or less abundant and the magnitude of this change between the two samples. These limitations can be overcome by using absolute abundance quantifications, which can allow for a more complete understanding of community dynamics by measuring variations in total microbial loads. Obtaining absolute abundance measurements is still technically challenging. Here, we developed synthetic DNA (synDNA) spike-ins that enable precise and cost-effective absolute quantification of microbiome data by adding defined amounts of synDNAs to the samples. We designed 10 synDNAs with the following features: 2,000-bp length, variable GC content (26, 36, 46, 56, or 66% GC), and negligible identity to sequences found in the NCBI database. Dilution pools were generated by mixing the 10 synDNAs at different concentrations. Shotgun metagenomic sequencing showed that the pools of synDNAs with different percentages of GC efficiently reproduced the serial dilution, showing high correlation (r = 0.96; R2 ≥ 0.94) and significance (P < 0.01). Furthermore, we demonstrated that the synDNAs can be used as DNA spike-ins to generate linear models and predict with high accuracy the absolute number of bacterial cells in complex microbial communities. IMPORTANCE The synDNAs designed in this study enable accurate and reproducible measurements of absolute amount and fold changes of bacterial species in complex microbial communities. The method proposed here is versatile and promising as it can be applied to bacterial communities or genomic features like genes and operons, in addition to being easily adaptable by other research groups at a low cost. We also made the synDNAs' sequences and the plasmids available to encourage future application of the proposed method in the study of microbial communities.

Zebra: Static and Dynamic Genome Cover Thresholds with Overlapping References.

Assigning taxonomy remains a challenging topic in microbiome studies, due largely to ambiguity of reads which overlap multiple reference genomes. With the Web of Life (WoL) reference database hosting 10,575 reference genomes and growing, the percentage of ambiguous reads will only increase. The resulting artifacts create both the illusion of co-occurrence and a long tail end of extraneous reference hits that confound interpretation. We introduce genome cover, the fraction of reference genome overlapped by reads, to distinguish these artifacts. We show how to dynamically predict genome cover by read count and examine our model in Staphylococcus aureus monoculture. Our modeling cleanly separates both S. aureus and true contaminants from the false artifacts of reference overlap. We next introduce saturated genome cover, the true fraction of a reference genome overlapped by sample contents. Genome cover may not saturate for low abundance or low prevalence bacteria. We assuage this worry with examination of a large human fecal data set. By compositing the metric across like samples, genome cover saturates even for rare species. We note that it is a threshold on saturated genome cover, not genome cover itself, which indicates a spurious reference hit or distant relative. We present Zebra, a method to compute and threshold the genome cover metric across like samples, a recurrence to estimate genome cover and confirm saturation, and provide guidance for choosing cover thresholds in real world scenarios. Standalone genome cover and integration into Woltka are available: https://github.com/biocore/zebra_filter, https://github.com/qiyunzhu/woltka. IMPORTANCE Taxonomic assignment, assigning sequences to specific taxonomic units, is a crucial processing step in microbiome analyses. Issues in taxonomic assignment affect interpretation of what microbes are present in each sample and may be associated with specific environmental or clinical conditions. Assigning importance to a particular taxon relies strongly on independence of assigned counts. The false inclusion of thousands of correlated taxa makes interpretation ambiguous, leading to underconstrained results which cannot be reproduced. The importance sometimes attached to implausible artifacts such as anthrax or bubonic plague is especially problematic. We show that the Zebra filter retrieves only the nearest relatives of sample contents enabling more reproducible and biologically plausible interpretation of metagenomic data.

Acetate reprograms gut microbiota during alcohol consumption.

Liver damage due to chronic alcohol use is among the most prevalent liver diseases. Alcohol consumption frequency is a strong factor of microbiota variance. Here we use isotope labeled [1-13C] ethanol, metagenomics, and metatranscriptomics in ethanol-feeding and intragastric mouse models to investigate the metabolic impacts of alcohol consumption on the gut microbiota. First, we show that although stable isotope labeled [1-13C] ethanol contributes to fatty acid pools in the liver, plasma, and cecum contents of mice, there is no evidence of ethanol metabolism by gut microbiota ex vivo under anaerobic conditions. Next, we observe through metatranscriptomics that the gut microbiota responds to ethanol-feeding by activating acetate dissimilation, not by metabolizing ethanol directly. We demonstrate that blood acetate concentrations are elevated during ethanol consumption. Finally, by increasing systemic acetate levels with glyceryl triacetate supplementation, we do not observe any impact on liver disease, but do induce similar gut microbiota alterations as chronic ethanol-feeding in mice. Our results show that ethanol is not directly metabolized by the gut microbiota, and changes in the gut microbiota linked to ethanol are a side effect of elevated acetate levels. De-trending for these acetate effects may be critical for understanding gut microbiota changes that cause alcohol-related liver disease.

The Ubiquitous Human Skin Commensal Staphylococcus hominis Protects against Opportunistic Pathogens.

Staphylococcus hominis is frequently isolated from human skin, and we hypothesize that it may protect the cutaneous barrier from opportunistic pathogens. We determined that S. hominis makes six unique autoinducing peptide (AIP) signals that inhibit the major virulence factor accessory gene regulator (agr) quorum sensing system of Staphylococcus aureus. We solved and confirmed the structures of three novel AIP signals in conditioned medium by mass spectrometry and then validated synthetic AIP activity against all S. aureus agr classes. Synthetic AIPs also inhibited the conserved agr system in a related species, Staphylococcus epidermidis. We determined the distribution of S. hominis agr types on healthy human skin and found S. hominis agr-I and agr-II were highly represented across subjects. Further, synthetic AIP-II was protective in vivo against S. aureus-associated dermonecrotic or epicutaneous injury. Together, these findings demonstrate that a ubiquitous colonizer of human skin has a fundamentally protective role against opportunistic damage. IMPORTANCE Human skin is home to a variety of commensal bacteria, including many species of coagulase-negative staphylococci (CoNS). While it is well established that the microbiota as a whole maintains skin homeostasis and excludes pathogens (i.e., colonization resistance), relatively little is known about the unique contributions of individual CoNS species to these interactions. Staphylococcus hominis is the second most frequently isolated CoNS from healthy skin, and there is emerging evidence to suggest that it may play an important role in excluding pathogens, including Staphylococcus aureus, from colonizing or infecting the skin. Here, we identified that S. hominis makes 6 unique peptide inhibitors of the S. aureus global virulence factor regulation system (agr). Additionally, we found that one of these peptides can prevent topical or necrotic S. aureus skin injury in a mouse model. Our results demonstrate a specific and broadly protective role for this ubiquitous, yet underappreciated skin commensal.

Host DNA Depletion in Saliva Samples for Improved Shotgun Metagenomics.

Host DNA makes up the majority of DNA in a saliva sample. Therefore, shotgun metagenomics can be an inefficient way to evaluate the microbial populations of saliva since often <10% of the sequencing reads are microbial. In this chapter, we describe a method to deplete human DNA from fresh or frozen saliva samples, allowing for more efficient shotgun metagenomic sequencing of the salivary microbial community.

Combining Functional Genomics and Whole-Genome Sequencing to Detect Antibiotic Resistance Genes in Bacterial Strains Co-Occurring Simultaneously in a Brazilian Hospital.

(1) Background: The rise of multi-antibiotic resistant bacteria represents an emergent threat to human health. Here, we investigate antibiotic resistance mechanisms in bacteria of several species isolated from an intensive care unit in Brazil. (2) Methods: We used whole-genome analysis to identify antibiotic resistance genes (ARGs) and plasmids in 34 strains of Gram-negative and Gram-positive bacteria, providing the first genomic description of Morganella morganii and Ralstonia mannitolilytica clinical isolates from South America. (3) Results: We identified a high abundance of beta-lactamase genes in resistant organisms, including seven extended-spectrum beta-lactamases (OXA-1, OXA-10, CTX-M-1, KPC, TEM, HYDRO, BLP) shared between organisms from different species. Additionally, we identified several ARG-carrying plasmids indicating the potential for a fast transmission of resistance mechanism between bacterial strains. Furthermore, we uncovered two pairs of (near) identical plasmids exhibiting multi-drug resistance. Finally, since many highly resistant strains carry several different ARGs, we used functional genomics to investigate which of them were indeed functional. In this sense, for three bacterial strains (Escherichia coli, Klebsiella pneumoniae, and M. morganii), we identified six beta-lactamase genes out of 15 predicted in silico as those mainly responsible for the resistance mechanisms observed, corroborating the existence of redundant resistance mechanisms in these organisms. (4) Conclusions: Systematic studies similar to the one presented here should help to prevent outbreaks of novel multidrug-resistant bacteria in healthcare facilities.

Intestinal α1-2-Fucosylation Contributes to Obesity and Steatohepatitis in Mice.

BACKGROUND & AIMS: Fucosyltransferase 2 (Fut2)-mediated intestinal α1- 2-fucosylation is important for host-microbe interactions and has been associated with several diseases, but its role in obesity and hepatic steatohepatitis is not known. The aim of this study was to investigate the role of Fut2 in a Western-style diet-induced mouse model of obesity and steatohepatitis. METHODS: Wild-type (WT) and Fut2-deficient littermate mice were used and features of the metabolic syndrome and steatohepatitis were assessed after 20 weeks of Western diet feeding. RESULTS: Intestinal α1-2-fucosylation was suppressed in WT mice after Western diet feeding, and supplementation of α1-2-fucosylated glycans exacerbated obesity and steatohepatitis in these mice. Fut2-deficient mice were protected from Western diet-induced features of obesity and steatohepatitis despite an increased caloric intake. These mice have increased energy expenditure and thermogenesis, as evidenced by a higher core body temperature. Protection from obesity and steatohepatitis associated with Fut2 deficiency is transmissible to WT mice via microbiota exchange; phenotypic differences between Western diet-fed WT and Fut2-deficient mice were reduced with antibiotic treatment. Fut2 deficiency attenuated diet-induced bile acid accumulation by altered relative abundance of bacterial enzyme 7-α-hydroxysteroid dehydrogenases metabolizing bile acids and by increased fecal excretion of secondary bile acids. This also was associated with increased intestinal farnesoid X receptor/fibroblast growth factor 15 signaling, which inhibits hepatic synthesis of bile acids. Dietary supplementation of α1-2-fucosylated glycans abrogates the protective effects of Fut2 deficiency. CONCLUSIONS: α1-2-fucosylation is an important host-derived regulator of intestinal microbiota and plays an important role for the pathogenesis of obesity and steatohepatitis in mice.

The sum is greater than the parts: exploiting microbial communities to achieve complex functions.

Multi-species microbial communities are ubiquitous in nature. The widespread prevalence of these communities is due to highly elaborated interactions among their members thereby accomplishing metabolic functions that are unattainable by individual members alone. Harnessing these communal capabilities is an emerging field in biotechnology. The rational intervention of microbial communities for the purpose of improved function has been facilitated in part by developments in multi-omics approaches, synthetic biology, and computational methods. Recent studies have demonstrated the benefits of rational interventions to human and animal health as well as agricultural productivity. Emergent technologies, such as in situ modification of complex microbial community and community metabolic modeling, represent an avenue to engineer sustainable microbial communities. In this opinion, we review relevant computational and experimental approaches to study and engineer microbial communities and discuss their potential for biotechnological applications.

A comparison of DNA/RNA extraction protocols for high-throughput sequencing of microbial communities.

One goal of microbial ecology researchers is to capture the maximum amount of information from all organisms in a sample. The recent COVID-19 pandemic, caused by the RNA virus SARS-CoV-2, has highlighted a gap in traditional DNA-based protocols, including the high-throughput methods the authors previously established as field standards. To enable simultaneous SARS-CoV-2 and microbial community profiling, the authors compared the relative performance of two total nucleic acid extraction protocols with the authors' previously benchmarked protocol. The authors included a diverse panel of environmental and host-associated sample types, including body sites commonly swabbed for COVID-19 testing. Here the authors present results comparing the cost, processing time, DNA and RNA yield, microbial community composition, limit of detection and well-to-well contamination between these protocols.

A comparison of DNA/RNA extraction protocols for high-throughput sequencing of microbial communities.

One goal among microbial ecology researchers is to capture the maximum amount of information from all organisms in a sample. The recent COVID-19 pandemic, caused by the RNA virus SARS-CoV-2, has highlighted a gap in traditional DNA-based protocols, including the high-throughput methods we previously established as field standards. To enable simultaneous SARS-CoV-2 and microbial community profiling, we compare the relative performance of two total nucleic acid extraction protocols and our previously benchmarked protocol. We included a diverse panel of environmental and host-associated sample types, including body sites commonly swabbed for COVID-19 testing. Here we present results comparing the cost, processing time, DNA and RNA yield, microbial community composition, limit of detection, and well-to-well contamination, between these protocols. Accession numbers: Raw sequence data were deposited at the European Nucleotide Archive (accession#: ERP124610) and raw and processed data are available at Qiita (Study ID: 12201). All processing and analysis code is available on GitHub ( github.com/justinshaffer/Extraction_test_MagMAX ). Methods summary: To allow for downstream applications involving RNA-based organisms such as SARS-CoV-2, we compared the two extraction protocols designed to extract DNA and RNA against our previously established protocol for extracting only DNA for microbial community analyses. Across 10 diverse sample types, one of the two protocols was equivalent or better than our established DNA-based protocol. Our conclusion is based on per-sample comparisons of DNA and RNA yield, the number of quality sequences generated, microbial community alpha- and beta-diversity and taxonomic composition, the limit of detection, and extent of well-to-well contamination.

Bacterial modification of the host glycosaminoglycan heparan sulfate modulates SARS-CoV-2 infectivity.

The human microbiota has a close relationship with human disease and it remodels components of the glycocalyx including heparan sulfate (HS). Studies of the severe acute respiratory syndrome coronavirus (SARS-CoV-2) spike protein receptor binding domain suggest that infection requires binding to HS and angiotensin converting enzyme 2 (ACE2) in a codependent manner. Here, we show that commensal host bacterial communities can modify HS and thereby modulate SARS-CoV-2 spike protein binding and that these communities change with host age and sex. Common human-associated commensal bacteria whose genomes encode HS-modifying enzymes were identified. The prevalence of these bacteria and the expression of key microbial glycosidases in bronchoalveolar lavage fluid (BALF) was lower in adult COVID-19 patients than in healthy controls. The presence of HS-modifying bacteria decreased with age in two large survey datasets, FINRISK 2002 and American Gut, revealing one possible mechanism for the observed increase in COVID-19 susceptibility with age. In vitro , bacterial glycosidases from unpurified culture media supernatants fully blocked SARS-CoV-2 spike binding to human H1299 protein lung adenocarcinoma cells. HS-modifying bacteria in human microbial communities may regulate viral adhesion, and loss of these commensals could predispose individuals to infection. Understanding the impact of shifts in microbial community composition and bacterial lyases on SARS-CoV-2 infection may lead to new therapeutics and diagnosis of susceptibility.

Transcriptional profiling of lung macrophages identifies a predictive signature for inflammatory lung disease in preterm infants.

Lung macrophages mature after birth, placing newborn infants, particularly those born preterm, within a unique window of susceptibility to disease. We hypothesized that in preterm infants, lung macrophage immaturity contributes to the development of bronchopulmonary dysplasia (BPD), the most common serious complication of prematurity. By measuring changes in lung macrophage gene expression in preterm patients at risk of BPD, we show here that patients eventually developing BPD had higher inflammatory mediator expression even on the first day of life. Surprisingly, the ex vivo response to LPS was similar across all samples. Our analysis did however uncover macrophage signature genes whose expression increased in the first week of life specifically in patients resilient to disease. We propose that these changes describe the dynamics of human lung macrophage differentiation. Our study therefore provides new mechanistic insight into both neonatal lung disease and human developmental immunology.

Interplay of Staphylococcal and Host Proteases Promotes Skin Barrier Disruption in Netherton Syndrome.

Netherton syndrome (NS) is a monogenic skin disease resulting from loss of function of lymphoepithelial Kazal-type-related protease inhibitor (LEKTI-1). In this study we examine if bacteria residing on the skin are influenced by the loss of LEKTI-1 and if interaction between this human gene and resident bacteria contributes to skin disease. Shotgun sequencing of the skin microbiome demonstrates that lesional skin of NS subjects is dominated by Staphylococcus aureus (S. aureus) and Staphylococcus epidermidis (S. epidermidis). Isolates of either species from NS subjects are able to induce skin inflammation and barrier damage on mice. These microbes promote skin inflammation in the setting of LEKTI-1 deficiency due to excess proteolytic activity promoted by S. aureus phenol-soluble modulin α as well as increased bacterial proteases staphopain A and B from S. aureus or EcpA from S. epidermidis. These findings demonstrate the critical need for maintaining homeostasis of host and microbial proteases to prevent a human skin disease.

Genomic and Transcriptomic Evidence Supports Methane Metabolism in Archaeoglobi.

Euryarchaeal lineages have been believed to have a methanogenic last common ancestor. However, members of euryarchaeal Archaeoglobi have long been considered nonmethanogenic and their evolutionary history remains elusive. Here, three high-quality metagenomic-assembled genomes (MAGs) retrieved from high-temperature oil reservoir and hot springs, together with three newly assembled Archaeoglobi MAGs from previously reported hot spring metagenomes, are demonstrated to represent a novel genus of Archaeoglobaceae, "Candidatus Methanomixophus." All "Ca Methanomixophus" MAGs encode an M methyltransferase (MTR) complex and a traditional type of methyl-coenzyme M reductase (MCR) complex, which is different from the divergent MCR complexes found in "Ca Polytropus marinifundus." In addition, "Ca Methanomixophus dualitatem" MAGs preserve the genomic capacity for dissimilatory sulfate reduction. Comparative phylogenetic analysis supports a laterally transferred origin for an MCR complex and vertical heritage of the MTR complex in this lineage. Metatranscriptomic analysis revealed concomitant in situ activity of hydrogen-dependent methylotrophic methanogenesis and heterotrophic fermentation within populations of "Ca Methanomixophus hydrogenotrophicum" in a high-temperature oil reservoir.IMPORTANCE Current understanding of the diversity, biology, and ecology of Archaea is very limited, especially considering how few of the known phyla have been cultured or genomically explored. The reconstruction of "Ca Methanomixophus" MAGs not only expands the known range of metabolic versatility of the members of Archaeoglobi but also suggests that the phylogenetic distribution of MCR and MTR complexes is even wider than previously anticipated.

Author Correction: Environmental stimuli drive a transition from cooperation to competition in synthetic phototrophic communities.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.