IMG: the Integrated Microbial Genomes database and comparative analysis system

Evolutionary stasis of a deep subsurface microbial lineage

Sulfate-reducing bacteria Candidatus Desulforudis audaxviator (CDA) were originally discovered in deep fracture fluids accessed via South African gold mines and have since been found in geographically widespread deep subsurface locations. In order to constrain models for subsurface microbial evolution, we compared CDA genomes from Africa, North America and Eurasia using single cell genomics. Unexpectedly, 126 partial single amplified genomes from the three continents, a complete genome from of an isolate from Eurasia, and metagenome-assembled genomes from Africa and Eurasia shared >99.2% average nucleotide identity, low frequency of SNP's, and near-perfectly conserved prophages and CRISPRs. Our analyses reject sample cross-contamination, recent natural dispersal, and unusually strong purifying selection as likely explanations for these unexpected results. We therefore conclude that the analyzed CDA populations underwent only minimal evolution since their physical separation, potentially as far back as the breakup of Pangea between 165 and 55 Ma ago. High-fidelity DNA replication and repair mechanisms are the most plausible explanation for the highly conserved genome of CDA. CDA presents a stark contrast to the current model organisms in microbial evolutionary studies, which often develop adaptive traits over far shorter periods of time.

Genomic Variation Influences Methanothermococcus Fitness in Marine Hydrothermal Systems

Hydrogenotrophic methanogens are ubiquitous chemoautotrophic archaea inhabiting globally distributed deep-sea hydrothermal vent ecosystems and associated subseafloor niches within the rocky subseafloor, yet little is known about how they adapt and diversify in these habitats. To determine genomic variation and selection pressure within methanogenic populations at vents, we examined five Methanothermococcus single cell amplified genomes (SAGs) in conjunction with 15 metagenomes and 10 metatranscriptomes from venting fluids at two geochemically distinct hydrothermal vent fields on the Mid-Cayman Rise in the Caribbean Sea. We observed that some Methanothermococcus lineages and their transcripts were more abundant than others in individual vent sites, indicating differential fitness among lineages. The relative abundances of lineages represented by SAGs in each of the samples matched phylogenetic relationships based on single-copy universal genes, and genes related to nitrogen fixation and the CRISPR/Cas immune system were among those differentiating the clades. Lineages possessing these genes were less abundant than those missing that genomic region. Overall, patterns in nucleotide variation indicated that the population dynamics of Methanothermococcus were not governed by clonal expansions or selective sweeps, at least in the habitats and sampling times included in this study. Together, our results show that although specific lineages of Methanothermococcus co-exist in these habitats, some outcompete others, and possession of accessory metabolic functions does not necessarily provide a fitness advantage in these habitats in all conditions. This work highlights the power of combining single-cell, metagenomic, and metatranscriptomic datasets to determine how evolution shapes microbial abundance and diversity in hydrothermal vent ecosystems.

MicFunPred: A conserved approach to predict functional profiles from 16S rRNA gene sequence data

The 16S rRNA gene amplicon sequencing is a popular technique that provides accurate characterization of microbial taxonomic abundances but does not provide any functional information. Several tools are available to predict functional profiles based on 16S rRNA gene sequence data that use different genome databases and approaches. As variable regions of partially-sequenced 16S rRNA gene cannot resolve taxonomy accurately beyond the genus level, these tools may give inflated results. Here, we developed 'MicFunPred', which uses a novel approach to derive imputed metagenomes based on a set of core genes only, thereby minimizing false-positive predictions. On simulated datasets, MicFunPred showed the lowest False Positive Rate (FPR) with mean Spearman's correlation of 0.89 (SD = 0.03), while on seven real datasets the mean correlation was 0.75 (SD = 0.08). MicFunPred was found to be faster with low computational requirements and performed better or comparable when compared with other tools.

Transcriptomics Analysis Uncovers Transient Ceftazidime Tolerance in Burkholderia Biofilms

Burkholderia pseudomallei is an etiological agent of melioidosis, a severe community-acquired infectious disease. B. pseudomallei strain K96243 is sensitive to the drug ceftazidime (CAZ), but has been shown to exhibit transient CAZ tolerance when in a biofilm form. To investigate an observed shift in gene expression profile during CAZ tolerance condition and to better understand the mechanistic aspects of this transient tolerance, RNA-sequencing was performed on B. pseudomallei K96243 from the following three states: planktonic, biofilm, and planktonic shedding. Results indicated that the expression of 651 genes (10.97%) were significantly changed in both biofilm (resistant) and planktonic shedding (sensitive) cells in comparison to the planktonic state. The top four highly expressed genes identified in both states are associated with nitrosative stress response (BPSL2368), Fe-S homeostasis (BPSL2369), and nitrate respiration (BPSS1154 and BPSS1158). Additionally, five orthologous genes, BPSL2370-BPSL2374, implicated in Fe-S cluster biogenesis, and another gene, BPSL2863, involved in DNA-binding of the stress protein ferritin, were shown to increase expression by RT-qPCR. The shift in gene expression was especially prominent at the late stages of biofilm growth (72 and 96 h), specifically in the biofilm-challenged CAZ survivor cells. This suggested that in response to stress in a biofilm, differential expression of these genes may support development of the CAZ tolerance in Burkholderia. The application of iron chelator deferoxamine (DFO) to the biofilm caused a significant reduction in biofilm formation and associated CAZ tolerance. Therefore, the shift in Fe-S metabolism when B. pseudomallei is in a biofilm may help stabilize the levels of reactive oxygen species (ROS), thereby limiting tolerance to CAZ.

Using Functional Annotations to Study Pairwise Interactions in Urinary Tract Infection Communities

The behaviour of microbial communities depends on environmental factors and on the interactions of the community members. This is also the case for urinary tract infection (UTI) microbial communities. Here, we devise a computational approach that uses indices of complementarity and competition based on metabolic gene annotation to rapidly predict putative interactions between pair of organisms with the aim to explain pairwise growth effects. We apply our method to 66 genomes selected from online databases, which belong to 6 genera representing members of UTI communities. This resulted in a selection of metabolic pathways with high correlation for each pairwise combination between a complementarity index and the experimentally derived growth data. Our results indicated that Enteroccus spp. were most complemented in its metabolism by the other members of the UTI community. This suggests that the growth of Enteroccus spp. can potentially be enhanced by complementary metabolites produced by other community members. We tested a few putative predicted interactions by experimental supplementation of the relevant predicted metabolites. As predicted by our method, folic acid supplementation led to the increase in the population density of UTI Enterococcus isolates. Overall, we believe our method is a rapid initial in silico screening for the prediction of metabolic interactions in microbial communities.

Distribution of Denitrification among Haloarchaea: A Comprehensive Study

Microorganisms from the Halobacteria class, also known as haloarchaea, inhabit a wide range of ecosystems of which the main characteristic is the presence of high salt concentration. These environments together with their microbial communities are not well characterized, but some of the common features that they share are high sun radiation and low availability of oxygen. To overcome these stressful conditions, and more particularly to deal with oxygen limitation, some microorganisms drive alternative respiratory pathways such as denitrification. In this paper, denitrification in haloarchaea has been studied from a phylogenetic point of view. It has been demonstrated that the presence of denitrification enzymes is a quite common characteristic in Halobacteria class, being nitrite reductase and nitric oxide reductase the enzymes with higher co-occurrence, maybe due to their possible role not only in denitrification, but also in detoxification. Moreover, copper-nitrite reductase (NirK) is the only class of respiratory nitrite reductase detected in these microorganisms up to date. The distribution of this alternative respiratory pathway and their enzymes among the families of haloarchaea has also been discussed and related with the environment in which they constitute the major populations. Complete denitrification phenotype is more common in some families like Haloarculaceae and Haloferacaceae, whilst less common in families such as Natrialbaceae and Halorubraceae.

Genome Annotation of Poly(lactic acid) Degrading Pseudomonas aeruginosa, Sphingobacterium sp. and Geobacillus sp

Pseudomonas aeruginosa and Sphingobacterium sp. are well known for their ability to decontaminate many environmental pollutants while Geobacillus sp. have been exploited for their thermostable enzymes. This study reports the annotation of genomes of P. aeruginosa S3, Sphingobacterium S2 and Geobacillus EC-3 that were isolated from compost, based on their ability to degrade poly(lactic acid), PLA. Draft genomes of the strains were assembled from Illumina reads, annotated and viewed with the aim of gaining insight into the genetic elements involved in degradation of PLA. The draft genome of Sphinogobacterium strain S2 (435 contigs) was estimated at 5,604,691 bp and the draft genome of P. aeruginosa strain S3 (303 contigs) was estimated at 6,631,638 bp. The draft genome of the thermophile Geobacillus strain EC-3 (111 contigs) was estimated at 3,397,712 bp. A total of 5385 (60% with annotation), 6437 (80% with annotation) and 3790 (74% with annotation) protein-coding genes were predicted for strains S2, S3 and EC-3, respectively. Catabolic genes for the biodegradation of xenobiotics, aromatic compounds and lactic acid as well as the genes attributable to the establishment and regulation of biofilm were identified in all three draft genomes. Our results reveal essential genetic elements that facilitate PLA metabolism at mesophilic and thermophilic temperatures in these three isolates.

The ecological impacts of multiple environmental stressors on coastal biofilm bacteria

Ecological communities are increasingly exposed to multiple interacting stressors. For example, warming directly affects the physiology of organisms, eutrophication stimulates the base of the food web, and harvesting larger organisms for human consumption dampens top-down control. These stressors often combine in the natural environment with unpredictable results. Bacterial communities in coastal ecosystems underpin marine food webs and provide many important ecosystem services (e.g. nutrient cycling and carbon fixation). Yet, how microbial communities will respond to a changing climate remains uncertain. Thus, we used marine mesocosms to examine the impacts of warming, nutrient enrichment, and altered top-predator population size structure (common shore crab) on coastal microbial biofilm communities in a crossed experimental design. Warming increased bacterial α-diversity (18% increase in species richness and 67% increase in evenness), but this was countered by a decrease in α-diversity with nutrient enrichment (14% and 21% decrease for species richness and evenness, respectively). Thus, we show some effects of these stressors could cancel each other out under climate change scenarios. Warming and top-predator population size structure both affected bacterial biofilm community composition, with warming increasing the abundance of bacteria capable of increased mineralization of dissolved and particulate organic matter, such as Flavobacteriia, Sphingobacteriia, and Cytophagia. However, the community shifts observed with warming depended on top-predator population size structure, with Sphingobacteriia increasing with smaller crabs and Cytophagia increasing with larger crabs. These changes could alter the balance between mineralization and carbon sequestration in coastal ecosystems, leading to a positive feedback loop between warming and CO₂ production. Our results highlight the potential for warming to disrupt microbial communities and biogeochemical cycling in coastal ecosystems, and the importance of studying these effects in combination with other environmental stressors.

Dynamics of Baltic Sea phages driven by environmental changes

Phage predation constitutes a major mortality factor for bacteria in aquatic ecosystems, and thus, directly impacts nutrient cycling and microbial community dynamics. Yet, the population dynamics of specific phages across time scales from days to months remain largely unexplored, which limits our understanding of their influence on microbial succession. To investigate temporal changes in diversity and abundance of phages infecting particular host strains, we isolated 121 phage strains that infected three bacterial hosts during a Baltic Sea mesocosm experiment. Genome analysis revealed a novel Flavobacterium phage genus harboring gene sets putatively coding for synthesis of modified nucleotides and glycosylation of bacterial cell surface components. Another novel phage genus revealed a microdiversity of phage species that was largely maintained during the experiment and across mesocosms amended with different nutrients. In contrast to the newly described Flavobacterium phages, phages isolated from a Rheinheimera strain were highly similar to previously isolated genotypes, pointing to genomic consistency in this population. In the mesocosm experiment, the investigated phages were mainly detected after a phytoplankton bloom peak. This concurred with recurrent detection of the phages in the Baltic Proper during summer months, suggesting an influence on the succession of heterotrophic bacteria associated with phytoplankton blooms.

Oral microbial dysbiosis and its performance in predicting oral cancer

Dysbiosis of oral microbiome may dictate the progression of oral squamous cell carcinoma (OSCC). Yet, the composition of oral microbiome fluctuates by saliva and distinct sites of oral cavity and is affected by risky behaviors (smoking, drinking and betel quid chewing) and individuals' oral health condition. To characterize the disturbances in the oral microbial population mainly due to oral tumorigenicity, we profiled the bacteria within the surface of OSCC lesion and its contralateral normal tissue from discovery (n = 74) and validation (n = 42) cohorts of male patients with cancers of the buccal mucosa. Significant alterations in the bacterial diversity and relative abundance of specific oral microbiota (most profoundly, an enrichment for genus Fusobacterium and the loss of genus Streptococcus in the tumor sites) were identified. Functional prediction of oral microbiome shown that microbial genes related to the metabolism of terpenoids and polyketides were differentially enriched between the control and tumor groups, indicating a functional role of oral microbiome in formulating a tumor microenvironment via attenuated biosynthesis of secondary metabolites with anti-cancer effects. Furthermore, the vast majority of microbial signatures detected in the discovery cohort was generalized well to the independent validation cohort, and the clinical validity of these OSCC-associated microbes was observed and successfully replicated. Overall, our analyses reveal signatures (a profusion of Fusobacterium nucleatum CTI-2 and a decrease in Streptococcus pneumoniae) and functions (decreased production of tumor-suppressive metabolites) of oral microbiota related to oral cancer.

Alterations in bile acid metabolizing gut microbiota and specific bile acid genes as a precision medicine to subclassify NAFLD

Multiple mechanisms for the gut microbiome contributing to the pathogenesis of nonalcoholic fatty liver disease (NAFLD) have been implicated. Here, we aim to investigate the contribution and potential application for altered bile acids (BA) metabolizing microbes in NAFLD by post hoc analysis of whole metagenome sequencing (WMS) data. The discovery cohort consisted of 86 well-characterized patients with biopsy-proven NAFLD and 38 healthy controls. Assembly-based analysis was performed to identify BA-metabolizing microbes. Statistical tests, feature selection, and microbial coabundance analysis were integrated to identify microbial alterations and markers in NAFLD. An independent validation cohort was subjected to similar analyses. NAFLD microbiota exhibited decreased diversity and microbial associations. We established a classifier model with 53 differential species exhibiting a robust diagnostic accuracy [area under the receiver-operator curve (AUC) = 0.97] for detecting NAFLD. Next, eight important differential pathway markers including secondary BA biosynthesis were identified. Specifically, increased abundance of 7α-hydroxysteroid dehydrogenase (7α-HSDH), 3α-hydroxysteroid dehydrogenase (baiA), and bile acid-coenzyme A ligase (baiB) was detected in NAFLD. Furthermore, 10 of 50 BA-metabolizing metagenome-assembled genomes (MAGs) from Bacteroides ovatus and Eubacterium biforme were dominant in NAFLD and interplayed as a synergetic ecological guild. Importantly, two subtypes of patients with NAFLD were observed according to secondary BA metabolism potentials. Elevated capability for secondary BA biosynthesis was also observed in the validation cohort. These bacterial BA-metabolizing genes and microbes identified in this study may serve as disease markers. Microbial differences in BA-metabolism and strain-specific differences among patients highlight the potential for precision medicine in NAFLD treatment.

The extracellular contractile injection system is enriched in environmental microbes and associates with numerous toxins

The extracellular Contractile Injection System (eCIS) is a toxin-delivery particle that evolved from a bacteriophage tail. Four eCISs have previously been shown to mediate interactions between bacteria and their invertebrate hosts. Here, we identify eCIS loci in 1,249 bacterial and archaeal genomes and reveal an enrichment of these loci in environmental microbes and their apparent absence from mammalian pathogens. We show that 13 eCIS-associated toxin genes from diverse microbes can inhibit the growth of bacteria and/or yeast. We identify immunity genes that protect bacteria from self-intoxication, further supporting an antibacterial role for some eCISs. We also identify previously undescribed eCIS core genes, including a conserved eCIS transcriptional regulator. Finally, we present our data through an extensive eCIS repository, termed eCIStem. Our findings support eCIS as a toxin-delivery system that is widespread among environmental prokaryotes and likely mediates antagonistic interactions with eukaryotes and other prokaryotes.

Radical genome remodelling accompanied the emergence of a novel host-restricted bacterial pathogen

The emergence of new pathogens is a major threat to public and veterinary health. Changes in bacterial habitat such as a switch in host or disease tropism are typically accompanied by genetic diversification. Staphylococcus aureus is a multi-host bacterial species associated with human and livestock infections. A microaerophilic subspecies, Staphylococcus aureus subsp. anaerobius, is responsible for Morel’s disease, a lymphadenitis restricted to sheep and goats. However, the evolutionary history of S. aureus subsp. anaerobius and its relatedness to S. aureus are unknown. Population genomic analyses of clinical S. aureus subsp. anaerobius isolates revealed a highly conserved clone that descended from a S. aureus progenitor about 1000 years ago before differentiating into distinct lineages that contain African and European isolates. S. aureus subsp. anaerobius has undergone limited clonal expansion, with a restricted population size, and an evolutionary rate 10-fold slower than S. aureus. The transition to its current restricted ecological niche involved acquisition of a pathogenicity island encoding a ruminant host-specific effector of abscess formation, large chromosomal re-arrangements, and the accumulation of at least 205 pseudogenes, resulting in a highly fastidious metabolism. Importantly, expansion of ~87 insertion sequences (IS) located largely in intergenic regions provided distinct mechanisms for the control of expression of flanking genes, including a novel mechanism associated with IS-mediated anti-anti-sense decoupling of ancestral gene repression. Our findings reveal the remarkable evolutionary trajectory of a host-restricted bacterial pathogen that resulted from extensive remodelling of the S. aureus genome through an array of diverse mechanisms in parallel.

The emergence of new pathogens is a major threat to public and veterinary health. Some bacteria such as Staphylococcus aureus, have the capacity to infect many different host species including humans and livestock while others such as the closely-related S. aureus subsp. anaerobius, associated with a single type of pathology called Morel’s disease in small ruminants, are highly niche-restricted. However, our understanding of the genetic basis for such differences in bacterial host-tropism is very limited. Here, we discovered that S. aureus subsp. anaerobius evolved from an S. aureus ancestor and underwent an array of extensive changes to its genome that accompanied the transition to its current restricted lifestyle. We observed genome decay involving loss of function of hundreds of genes, large intra-chromosomal rearrangements affecting most of the genome, acquisition of a pathogenicity island, and expansion of large numbers of insertion sequences that are inserted at intergenic sites around the genome. Importantly, we found that IS elements affect the expression of neighbouring genes in different ways including a novel mechanism of IS-enabled disruption of ancestral gene repression. Taken together, we provide a remarkable example of radical genomic changes associated with evolutionary transition from a multi-host to highly restricted host ecology.

Functional Analysis of Phenazine Biosynthesis Genes in Burkholderia spp

Burkholderia encompasses a group of ubiquitous Gram-negative bacteria that includes numerous saprophytes as well as species that cause infections in animals, immunocompromised patients, and plants. Some species of Burkholderia produce colored, redox-active secondary metabolites called phenazines. Phenazines contribute to competitiveness, biofilm formation, and virulence in the opportunistic pathogen Pseudomonas aeruginosa, but knowledge of their diversity, biosynthesis, and biological functions in Burkholderia is lacking. In this study, we screened publicly accessible genome sequence databases and identified phenazine biosynthesis genes in multiple strains of the Burkholderia cepacia complex, some isolates of the B. pseudomallei clade, and the plant pathogen B. glumae We then focused on B. lata ATCC 17760 to reveal the organization and function of genes involved in the production of dimethyl 4,9-dihydroxy-1,6-phenazinedicarboxylate. Using a combination of isogenic mutants and plasmids carrying different segments of the phz locus, we characterized three novel genes involved in the modification of the phenazine tricycle. Our functional studies revealed a connection between the presence and amount of phenazines and the dynamics of biofilm growth in flow cell and static experimental systems but at the same time failed to link the production of phenazines with the capacity of Burkholderia to kill fruit flies and rot onions.IMPORTANCE Although the production of phenazines in Burkholderia was first reported almost 70 years ago, the role these metabolites play in the biology of these economically important microorganisms remains poorly understood. Our results revealed that the phenazine biosynthetic pathway in Burkholderia has a complex evolutionary history, which likely involved horizontal gene transfers among several distantly related groups of organisms. The contribution of phenazines to the formation of biofilms suggests that Burkholderia, like fluorescent pseudomonads, may benefit from the unique redox-cycling properties of these versatile secondary metabolites.

Programmable RNA editing with compact CRISPR-Cas13 systems from uncultivated microbes

Competitive coevolution between microbes and viruses has led to the diversification of CRISPR-Cas defense systems against infectious agents. By analyzing metagenomic terabase datasets, we identified two compact families (775 to 803 amino acids (aa)) of CRISPR-Cas ribonucleases from hypersaline samples, named Cas13X and Cas13Y. We engineered Cas13X.1 (775 aa) for RNA interference experiments in mammalian cell lines. We found Cas13X.1 could tolerate single-nucleotide mismatches in RNA recognition, facilitating prophylactic RNA virus inhibition. Moreover, a minimal RNA base editor, composed of engineered deaminase (385 aa) and truncated Cas13X.1 (445 aa), exhibited robust editing efficiency and high specificity to induce RNA base conversions. Our results suggest that there exist untapped bacterial defense systems in natural microbes that can function efficiently in mammalian cells, and thus potentially are useful for RNA-editing-based research.

Microbial Communities in Full-Scale Wastewater Treatment Systems Exhibit Deterministic Assembly Processes and Functional Dependency over Time

Microbial communities constitute the core component of biological wastewater treatment processes. We conducted a meta-analysis based on the 16S rRNA gene of temporal samples obtained from diverse full-scale activated sludge and anaerobic digestion systems treating municipal and industrial wastewater (collected in this study and published previously) to investigate their community assembly mechanism and functional traits over time, which are not currently well understood. The influent composition was found to be the main driver of the microbial community's composition, and relatively large proportions of specialist (26.1% and 18.6%) and transient taxa (67.2% and 68.1%) were estimated in both systems. Deterministic processes, especially homogeneous selection events (accounting for >53.8% of assembly events), were consistently identified as the dominant microbial community assembly mechanisms in both systems over time. Significant and strong correlations (Pearson's r = 0.51-0.92) were detected between the dynamics of the temporal community and the functional compositions in both systems, which suggests functional dependency. In contrast, the occurrence of sludge bulking and foaming in the activated sludge system led to an increase in stochastic assembly processes (i.e., limited dispersal and undominated events), a shift toward functional redundancy and less community diversity, a decreased community niche breadth index, and a more compact co-association network. This study illustrates that the mechanism of microbial community assembly and functional traits over time can be used to diagnose system performance and provide information on potential system malfunction.

Next generation sequencing approaches to evaluate water and wastewater quality

The emergence of next generation sequencing (NGS) is revolutionizing the potential to address complex microbiological challenges in the water industry. NGS technologies can provide holistic insight into microbial communities and their functional capacities in water and wastewater systems, thus eliminating the need to develop a new assay for each target organism or gene. However, several barriers have hampered wide-scale adoption of NGS by the water industry, including cost, need for specialized expertise and equipment, challenges with data analysis and interpretation, lack of standardized methods, and the rapid pace of development of new technologies. In this critical review, we provide an overview of the current state of the science of NGS technologies as they apply to water, wastewater, and recycled water. In addition, a systematic literature review was conducted in which we identified over 600 peer-reviewed journal articles on this topic and summarized their contributions to six key areas relevant to the water and wastewater fields: taxonomic classification and pathogen detection, functional and catabolic gene characterization, antimicrobial resistance (AMR) profiling, bacterial toxicity characterization, Cyanobacteria and harmful algal bloom identification, and virus characterization. For each application, we have presented key trends, noteworthy advancements, and proposed future directions. Finally, key needs to advance NGS technologies for broader application in water and wastewater fields are assessed.

Root Exudates Alter the Expression of Diverse Metabolic, Transport, Regulatory, and Stress Response Genes in Rhizosphere Pseudomonas

Plants live in association with microorganisms that positively influence plant development, vigor, and fitness in response to pathogens and abiotic stressors. The bulk of the plant microbiome is concentrated belowground at the plant root-soil interface. Plant roots secrete carbon-rich rhizodeposits containing primary and secondary low molecular weight metabolites, lysates, and mucilages. These exudates provide nutrients for soil microorganisms and modulate their affinity to host plants, but molecular details of this process are largely unresolved. We addressed this gap by focusing on the molecular dialog between eight well-characterized beneficial strains of the Pseudomonas fluorescens group and Brachypodium distachyon, a model for economically important food, feed, forage, and biomass crops of the grass family. We collected and analyzed root exudates of B. distachyon and demonstrated the presence of multiple carbohydrates, amino acids, organic acids, and phenolic compounds. The subsequent screening of bacteria by Biolog Phenotype MicroArrays revealed that many of these metabolites provide carbon and energy for the Pseudomonas strains. RNA-seq profiling of bacterial cultures amended with root exudates revealed changes in the expression of genes encoding numerous catabolic and anabolic enzymes, transporters, transcriptional regulators, stress response, and conserved hypothetical proteins. Almost half of the differentially expressed genes mapped to the variable part of the strains’ pangenome, reflecting the importance of the variable gene content in the adaptation of P. fluorescens to the rhizosphere lifestyle. Our results collectively reveal the diversity of cellular pathways and physiological responses underlying the establishment of mutualistic interactions between these beneficial rhizobacteria and their plant hosts.

Cutting the line: manipulation of plant immunity by bacterial type III effector proteases

Pathogens and their hosts are engaged in an evolutionary arms race. Pathogen-derived effectors promote virulence by targeting components of a host's innate immune system, while hosts have evolved proteins that sense effectors and trigger a pathogen-specific immune response. Many bacterial effectors are translocated into host cells using type III secretion systems. Type III effector proteases irreversibly modify host proteins by cleavage of peptide bonds and are prevalent among both plant and animal bacterial pathogens. In plants, the study of model effector proteases has yielded important insights into the virulence mechanisms employed by pathogens to overcome their host's immune response, as well as into the mechanisms deployed by their hosts to detect these effector proteases and counteract their effects. In recent years, the study of a larger number of effector proteases, across a wider range of pathogens, has yielded novel insights into their functions and recognition. One key limitation that remains is the lack of methods to detect protease cleavage at the proteome-wide level. We review known substrates and mechanisms of plant pathogen type III effector proteases and compare their functions with those of known type III effector proteases of mammalian pathogens. Finally, we discuss approaches to uncover their function on a system-wide level.

CAZymes in Maribacter dokdonensis 62-1 From the Patagonian Shelf: Genomics and Physiology Compared to Related Flavobacteria and a Co-occurring Alteromonas Strain

Carbohydrate-active enzymes (CAZymes) are an important feature of bacteria in productive marine systems such as continental shelves, where phytoplankton and macroalgae produce diverse polysaccharides. We herein describe Maribacter dokdonensis 62–1, a novel strain of this flavobacterial species, isolated from alginate-supplemented seawater collected at the Patagonian continental shelf. M. dokdonensis 62–1 harbors a diverse array of CAZymes in multiple polysaccharide utilization loci (PUL). Two PUL encoding polysaccharide lyases from families 6, 7, 12, and 17 allow substantial growth with alginate as sole carbon source, with simultaneous utilization of mannuronate and guluronate as demonstrated by HPLC. Furthermore, strain 62-1 harbors a mixed-feature PUL encoding both ulvan- and fucoidan-targeting CAZymes. Core-genome phylogeny and pangenome analysis revealed variable occurrence of these PUL in related Maribacter and Zobellia strains, indicating specialization to certain “polysaccharide niches.” Furthermore, lineage- and strain-specific genomic signatures for exopolysaccharide synthesis possibly mediate distinct strategies for surface attachment and host interaction. The wide detection of CAZyme homologs in algae-derived metagenomes suggests global occurrence in algal holobionts, supported by sharing multiple adaptive features with the hydrolytic model flavobacterium Zobellia galactanivorans. Comparison with Alteromonas sp. 76-1 isolated from the same seawater sample revealed that these co-occurring strains target similar polysaccharides but with different genomic repertoires, coincident with differing growth behavior on alginate that might mediate ecological specialization. Altogether, our study contributes to the perception of Maribacter as versatile flavobacterial polysaccharide degrader, with implications for biogeochemical cycles, niche specialization and bacteria-algae interactions in the oceans.

Metaproteogenomic Profiling of Chemosynthetic Microbial Biofilms Reveals Metabolic Flexibility During Colonization of a Shallow-Water Gas Vent

Tor Caldara is a shallow-water gas vent located in the Mediterranean Sea, with active venting of CO₂ and H₂S. At Tor Caldara, filamentous microbial biofilms, mainly composed of Epsilon- and Gammaproteobacteria, grow on substrates exposed to the gas venting. In this study, we took a metaproteogenomic approach to identify the metabolic potential and in situ expression of central metabolic pathways at two stages of biofilm maturation. Our findings indicate that inorganic reduced sulfur species are the main electron donors and CO₂ the main carbon source for the filamentous biofilms, which conserve energy by oxygen and nitrate respiration, fix dinitrogen gas and detoxify heavy metals. Three metagenome-assembled genomes (MAGs), representative of key members in the biofilm community, were also recovered. Metaproteomic data show that metabolically active chemoautotrophic sulfide-oxidizing members of the Epsilonproteobacteria dominated the young microbial biofilms, while Gammaproteobacteria become prevalent in the established community. The co-expression of different pathways for sulfide oxidation by these two classes of bacteria suggests exposure to different sulfide concentrations within the biofilms, as well as fine-tuned adaptations of the enzymatic complexes. Taken together, our findings demonstrate a shift in the taxonomic composition and associated metabolic activity of these biofilms in the course of the colonization process.

Metabolic versatility of the nitrite-oxidizing bacterium Nitrospira marina and its proteomic response to oxygen-limited conditions

The genus Nitrospira is the most widespread group of nitrite-oxidizing bacteria and thrives in diverse natural and engineered ecosystems. Nitrospira marina Nb-295^T was isolated from the ocean over 30 years ago; however, its genome has not yet been analyzed. Here, we investigated the metabolic potential of N. marina based on its complete genome sequence and performed physiological experiments to test genome-derived hypotheses. Our data confirm that N. marina benefits from additions of undefined organic carbon substrates, has adaptations to resist oxidative, osmotic, and UV light-induced stress and low dissolved pCO₂, and requires exogenous vitamin B₁₂. In addition, N. marina is able to grow chemoorganotrophically on formate, and is thus not an obligate chemolithoautotroph. We further investigated the proteomic response of N. marina to low (∼5.6 µM) O₂ concentrations. The abundance of a potentially more efficient CO₂-fixing pyruvate:ferredoxin oxidoreductase (POR) complex and a high-affinity cbb₃-type terminal oxidase increased under O₂ limitation, suggesting a role in sustaining nitrite oxidation-driven autotrophy. This putatively more O₂-sensitive POR complex might be protected from oxidative damage by Cu/Zn-binding superoxide dismutase, which also increased in abundance under low O₂ conditions. Furthermore, the upregulation of proteins involved in alternative energy metabolisms, including Group 3b [NiFe] hydrogenase and formate dehydrogenase, indicate a high metabolic versatility to survive conditions unfavorable for aerobic nitrite oxidation. In summary, the genome and proteome of the first marine Nitrospira isolate identifies adaptations to life in the oxic ocean and provides insights into the metabolic diversity and niche differentiation of NOB in marine environments.

Characterization of a novel roseophage and the morphological and transcriptional response of the sponge symbiont Ruegeria AU67 to infection

Sponges have recently been recognized to contain complex communities of bacteriophages; however, little is known about how they interact with their bacterial hosts. Here, we isolated a novel phage, called Ruegeria phage Tedan, and characterized its impact on the bacterial sponge symbiont Ruegeria AU67 on a morphological and molecular level. Phage Tedan was structurally, genomically and phylogenetically characterized to be affiliated with the genus Xiamenvirus of the family Siphoviridae. Through microscopic observations and transcriptomic analysis, we show that phage Tedan upon infection induces a process leading to metabolic and morphological changes in its host. These changes would render Ruegeria AU67 better adapted to inhabit the sponge holobiont due to an improved utilization of ecologically relevant energy and carbon sources as well as a potential impediment of phagocytosis by the sponge through cellular enlargement. An increased survival or better growth of the bacterium in the sponge environment will likely benefit the phage reproduction. Our results point towards the possibility that phages from host-associated environments require, and have thus evolved, different strategies to interact with their host when compared to those phages from free-living or planktonic environments.

Ultrafine particles altered gut microbial population and metabolic profiles in a sex-specific manner in an obese mouse model

Emerging evidence has highlighted the connection between exposure to air pollution and the increased risk of obesity, metabolic syndrome, and comorbidities. Given the recent interest in studying the effects of ultrafine particle (UFP) on the health of obese individuals, this study examined the effects of gastrointestinal UFP exposure on gut microbial composition and metabolic function using an in vivo murine model of obesity in both sexes. UFPs generated from light-duty diesel engine combustion of petrodiesel (B0) and a petrodiesel/biodiesel fuel blend (80:20 v/v, B20) were administered orally. Multi-omics approaches, including liquid chromatography-mass spectrometry (LC-MS) based targeted metabolomics and 16S rRNA gene sequence analysis, semi-quantitatively compared the effects of 10-day UFP exposures on obese C57B6 mouse gut microbial population, changes in diversity and community function compared to a phosphate buffer solution (PBS) control group. Our results show that sex-specific differences in the gut microbial population in response to UFP exposure can be observed, as UFPs appear to have a differential impact on several bacterial families in males and females. Meanwhile, the alteration of seventy-five metabolites from the gut microbial metabolome varied significantly (ANOVA p < 0.05) across the PBS control, B0, and B20 groups. Multivariate analyses revealed that the fuel-type specific disruption to the microbial metabolome was observed in both sexes, with stronger disruptive effects found in females in comparison to male obese mice. Metabolic signatures of bacterial cellular oxidative stress, such as the decreased concentration of nucleotides and lipids and increased concentrations of carbohydrate, energy, and vitamin metabolites were detected. Furthermore, blood metabolites from the obese mice were differentially affected by the fuel types used to generate the UFPs (B0 vs. B20).

The cyanobacterial taxis protein HmpF regulates type IV pilus activity in response to light

Motility is ubiquitous in prokaryotic organisms including the photosynthetic cyanobacteria where surface motility powered by type 4 pili (T4P) is common and facilitates phototaxis to seek out favorable light environments. In cyanobacteria, chemotaxis-like systems are known to regulate motility and phototaxis. The characterized phototaxis systems rely on methyl-accepting chemotaxis proteins containing bilin-binding GAF domains capable of directly sensing light, and the mechanism by which they regulate the T4P is largely undefined. In this study we demonstrate that cyanobacteria possess a second, GAF-independent, means of sensing light to regulate motility and provide insight into how a chemotaxis-like system regulates the T4P motors. A combination of genetic, cytological, and protein-protein interaction analyses, along with experiments using the proton ionophore carbonyl cyanide m-chlorophenyl hydrazine, indicate that the Hmp chemotaxis-like system of the model filamentous cyanobacterium Nostoc punctiforme is capable of sensing light indirectly, possibly via alterations in proton motive force, and modulates direct interaction between the cyanobacterial taxis protein HmpF, and Hfq, PilT1, and PilT2 to regulate the T4P motors. Given that the Hmp system is widely conserved in cyanobacteria, and the finding from this study that orthologs of HmpF and T4P proteins from the distantly related model unicellular cyanobacterium Synechocystis sp. strain PCC6803 interact in a similar manner to their N. punctiforme counterparts, it is likely that this represents a ubiquitous means of regulating motility in response to light in cyanobacteria.

Scarcity of scale-free topology is universal across biochemical networks

Biochemical reactions underlie the functioning of all life. Like many examples of biology or technology, the complex set of interactions among molecules within cells and ecosystems poses a challenge for quantification within simple mathematical objects. A large body of research has indicated many real-world biological and technological systems, including biochemistry, can be described by power-law relationships between the numbers of nodes and edges, often described as "scale-free". Recently, new statistical analyses have revealed true scale-free networks are rare. We provide a first application of these methods to data sampled from across two distinct levels of biological organization: individuals and ecosystems. We analyze a large ensemble of biochemical networks including networks generated from data of 785 metagenomes and 1082 genomes (sampled from the three domains of life). The results confirm no more than a few biochemical networks are any more than super-weakly scale-free. Additionally, we test the distinguishability of individual and ecosystem-level biochemical networks and show there is no sharp transition in the structure of biochemical networks across these levels of organization moving from individuals to ecosystems. This result holds across different network projections. Our results indicate that while biochemical networks are not scale-free, they nonetheless exhibit common structure across different levels of organization, independent of the projection chosen, suggestive of shared organizing principles across all biochemical networks.

Mucilaginibacter aquatilis sp. nov., Mucilaginibacter arboris sp. nov., and Mucilaginibacter ginkgonis sp. nov., novel bacteria isolated from freshwater and tree bark

Three Gram-stain-negative, strictly aerobic, non-motile and rod-shaped bacterial strains, designated as HME9299^T, HMF7410^T and HMF7856^T, were isolated from freshwater and tree bark collected in Yong-in, Republic of Korea. Strains HME9299^T, HMF7410^T and HMF7856^T exhibited the highest 16S rRNA gene sequence similarities of 97.2, 94.4 and 96.4 % to Mucilaginibacter daejeonensis Jip 10^T, Mucilaginibacter terrae CCM 8645^T and Mucilaginibacter phyllosphaerae PP-F2F-G21^T, respectively. Among themselves, the values were 94.1-95.7 %. Phylogenetic analysis of the 16S rRNA gene sequences of the three isolates revealed that they belonged to the genus Mucilaginibacter within the family Sphingobacteriaceae. The predominant fatty acids of three strains were summed feature 3 (comprising C_16 : 1  ω7c and/or C_16 : 1  ω6c) and iso-C_15 : 0. Strain HME9299^T contained a relatively large amount of C_16 : 1  ω5c. The predominant respiratory quinone was menaquinone-7. The genome sizes of strains HME9299^T, HMF7410^T and HMF7856^T were 4.33, 4.16 and 3.68 Mbp, respectively, and their DNA G+C contents were 41.6, 38.4 and 43.9 mol%, respectively. Based on the results of the phenotypic, genotypic, chemotaxonomic and phylogenetic investigation, three novel species, Mucilaginibacter aquatilis sp. nov, Mucilaginibacter arboris sp. nov. and Mucilaginibacter ginkgonis sp. nov., are proposed. The type strains are HME9299^T (=KCTC 42122^T=DSM 29146^T), HMF7410^T (=KCTC 62464^T=NBRC 113227^T) and HMF7856^T (=KCTC 72782^T=NBRC 114275^T), respectively.

Differences in the Microbial Composition of Hemodialysis Patients Treated with and without β-Blockers

β-blockers are commonly prescribed to treat cardiovascular disease in hemodialysis patients. Beyond the pharmacological effects, β-blockers have potential impacts on gut microbiota, but no study has investigated the effect in hemodialysis patients. Hence, we aim to investigate the gut microbiota composition difference between β-blocker users and nonusers in hemodialysis patients. Fecal samples collected from hemodialysis patients (83 β-blocker users and 110 nonusers) were determined by 16S ribosomal RNA amplification sequencing. Propensity score (PS) matching was performed to control confounders. The microbial composition differences were analyzed by the linear discriminant analysis effect size, random forest, and zero-inflated Gaussian fit model. The α-diversity (Simpson index) was greater in β-blocker users with a distinct β-diversity (Bray-Curtis Index) compared to nonusers in both full and PS-matched cohorts. There was a significant enrichment in the genus Flavonifractor in β-blocker users compared to nonusers in full and PS-matched cohorts. A similar finding was demonstrated in random forest analysis. In conclusion, hemodialysis patients using β-blockers had a different gut microbiota composition compared to nonusers. In particular, the Flavonifractor genus was increased with β-blocker treatment. Our findings highlight the impact of β-blockers on the gut microbiota in hemodialysis patients.

The effect of crop species on DNase-producing bacteria in two soils

Purpose

Extracellular deoxyribonucleases (exDNases) from microbial origin contribute substantially to the restriction of extracellular DNA (exDNA) in the soil. Hence, it is imperative to understand the diversity of bacterial species capable of performing this important soil function and how crop species influence their dynamics in the soil. The present study investigates the occurrence of DNase-producing bacteria (DPB) in leachate samples obtained from soils in which the crop species of alfalfa (Medicago sativa L.), canola (Brassica napus L.), soybean (Glycine max [L.] Merr.) and wheat (Triticum aestivum L.) were raised in a growth room.

Methods

Selective media containing methyl green indicator was used to screen for DPB from leachate samples, whereas the 16S rRNA sequence analysis was employed to identify the isolates.

Results

The proportion of culturable DPB ranged between 5.72 and 40.01%; however, we did observe specific crop effects that shifted throughout the growing period. In general, higher proportions of exDNase producers were observed when the soils had lower nutrient levels. On using the 16S rRNA to classify the DPB isolates, most isolates were found to be members of the Bacillus genera, while other groups included Chryseobacterium, Fictibacillus, Flavobacterium, Microbacterium, Nubsella, Pseudomonas, Psychrobacillus, Rheinheimera, Serratia and Stenotrophomonas. Five candidate exDNase/nuclease-encoding proteins were also identified from Bacillus mycoides genomes using online databases.

Conclusion

Results from this study showed that crop species, growth stage and soil properties were important factors shaping the populations of DPB in leachate samples; however, soil properties seemed to have a greater influence on the trends observed on these bacterial populations. It may be possible to target soil indigenous bacteria that produce exDNases through management to decrease potential unintended effects of transgenes originating from genetically modified organisms (GMOs) or other introduced nucleic acid sequences in the environment.

A Pseudoalteromonas Clade with Remarkable Biosynthetic Potential

Pseudoalteromonas species produce a diverse range of biologically active compounds, including those biosynthesized by nonribosomal peptide synthetases (NRPSs) and polyketide synthases (PKSs). Here, we report the biochemical and genomic analysis of Pseudoalteromonas sp. strain HM-SA03, isolated from the blue-ringed octopus, Hapalochlaena sp. Genome mining for secondary metabolite pathways revealed seven putative NRPS/PKS biosynthesis gene clusters, including those for the biosynthesis of alterochromides, pseudoalterobactins, alteramides, and four novel compounds. Among these was a novel siderophore biosynthesis gene cluster with unprecedented architecture (NRPS-PKS-NRPS-PKS-NRPS-PKS-NRPS). Alterochromide production in HM-SA03 was also confirmed by liquid chromatography-mass spectrometry. An investigation of the biosynthetic potential of 42 publicly available Pseudoalteromonas genomes indicated that some of these gene clusters are distributed throughout the genus. Through the phylogenetic analysis, a particular subset of strains formed a clade with extraordinary biosynthetic potential, with an average density of 10 biosynthesis gene clusters per genome. In contrast, the majority of Pseudoalteromonas strains outside this clade contained an average of three clusters encoding complex biosynthesis. These results highlight the underexplored potential of Pseudoalteromonas as a source of new natural products.IMPORTANCE This study demonstrates that the Pseudoalteromonas strain HM-SA03, isolated from the venomous blue-ringed octopus, Hapalochalaena sp., is a biosynthetically talented organism, capable of producing alterochromides and potentially six other specialized metabolites. We identified a pseudoalterobactin biosynthesis gene cluster and proposed a pathway for the production of the associated siderophore. A novel siderophore biosynthesis gene cluster with unprecedented architecture was also identified in the HM-SA03 genome. Finally, we demonstrated that HM-SA03 belongs to a phylogenetic clade of strains with extraordinary biosynthetic potential. While our results do not support a role of HM-SA03 in Hapalochalaena sp. venom (tetrodotoxin) production, they emphasize the untapped potential of Pseudoalteromonas as a source of novel natural products.

Draft genome sequences of three filamentous cyanobacteria isolated from brackish habitats

Brackish cyanobacterial genome sequences are relatively rare. Here, we report the 5.5 Mbp, 5.8 Mbp and 6.1 Mbp draft genomes of Spirulina sp. CCY15215, Leptolyngbya sp. CCY15150 and Halomicronema sp. CCY15110 isolated from coastal microbial mats on the North Sea beach of the island of Schiermonnikoog in the Netherlands. Large scale phylogenomic analyses reveal that Spirulina sp. CCY15215 is a large cell diameter cyanobacterium, whereas Leptolyngbya sp. CCY15150 and Halomicronema sp. CCY15110 are the first reported brackish genomes belonging to the LPP clade consisting primarily of Leptolyngbya, Plectonema and Phormidium spp. Further genome mining divulges that all new draft genomes contain, ggpS and ggpP, the genes responsible for synthesising glucosylglycerol (GG), a compatible solute found in moderately salt-tolerant cyanobacteria.

Starvation Alters Gut Microbiome in Black Soldier Fly (Diptera: Stratiomyidae) Larvae

Unlike for vertebrates, the impact of starvation on the gut microbiome of invertebrates is poorly studied. Deciphering shifts in metabolically active associated bacterial communities in vertebrates has led to determining the role of the associated microbiome in the sensation of hunger and discoveries of associated regulatory mechanisms. From an invertebrate perspective, such as the black soldier fly, such information could lead to enhanced processes for optimized biomass production and waste conversion. Bacteria associated with food substrates of black soldier fly are known to impact corresponding larval life-history traits (e.g., larval development); however, whether black soldier fly larval host state (i.e., starved) impacts the gut microbiome is not known. In this study, we measured microbial community structural and functional shifts due to black soldier fly larvae starvation. Data generated demonstrate such a physiological state (i.e., starvation) does in fact impact both aspects of the microbiome. At the phylum level, community diversity decreased significantly during black soldier fly larval starvation (p = 0.0025). Genus level DESeq2 analysis identified five genera with significantly different relative abundance (q < 0.05) across the 24 and 48 H post initiation of starvation: Actinomyces, Microbacterium, Enterococcus, Sphingobacterium, and Leucobacter. Finally, we inferred potential gene function and significantly predicted functional KEGG Orthology (KO) abundance. We demonstrated the metabolically active microbial community structure and function could be influenced by host-feeding status. Such perturbations, even when short in duration (e.g., 24 H) could stunt larval growth and waste conversion due to lacking a full complement of bacteria and associated functions.

Intestinal Dysbiosis in Young Cystic Fibrosis Rabbits

Individuals with cystic fibrosis (CF) often experience gastrointestinal (GI) abnormalities. In recent years, the intestinal microbiome has been postulated as a contributor to the development of CF-associated GI complications, hence representing a potential therapeutic target for treatment. We recently developed a rabbit model of CF, which is shown to manifest many human patient-like pathological changes, including intestinal obstruction. Here, we investigated the feces microbiome in young CF rabbits in the absence of antibiotics treatment. Stool samples were collected from seven- to nine-week-old CF rabbits (n = 7) and age-matched wild-type (WT) rabbits (n = 6). Microbiomes were investigated by iTag sequencing of 16S rRNA genes, and functional profiles were predicted using PICRUSt. Consistent with reports of those in pediatric CF patients, the fecal microbiomes of CF rabbits are of lower richness and diversity than that of WT rabbits, with a marked taxonomic and inferred functional dysbiosis. Our work identified a new CF animal model with the manifestation of intestinal dysbiosis phenotype. This model system may facilitate the research and development of novel treatments for CF-associated gastrointestinal diseases.

Highly decomposed organic carbon mediates the assembly of soil communities with traits for the biodegradation of chlorinated pollutants

To improve biodegradation strategies for chlorinated pollutants, the roles of soil organic matter and microbial function need to be clarified. It was hypothesised that microbial degradation of specific organic fractions in soils enhance community metabolic capability to degrade chlorinated pollutants. This field study used historic records of dieldrin concentrations since 1988 and established relationships between dieldrin dissipation and soil carbon fractions together with bacterial and fungal diversity in surface soils of Kurosol and Chromosol. Sparse partial least squares analysis linked dieldrin dissipation to metabolic activities associated with the highly decomposed carbon fraction. Dieldrin dissipation, after three decades of natural attenuation, was associated with increased bacterial species fitness for the decomposition of recalcitrant carbon substrates including synthetic chlorinated pollutants. These metabolic capabilities were linked to the decomposed carbon fraction, an important driver for the microbial community and function. Common bacterial traits among taxonomic groups enriched in samples with high dieldrin dissipation included their slow growth, large genome and complex metabolism which supported the notion that metabolic strategies for dieldrin degradation evolved in an energy-low soil environment. The findings provide new perspectives for bioremediation strategies and suggest that soil management should aim at stimulating metabolism at the decomposed, fine carbon fraction.

Analysis of the Yarrowia lipolytica proteome reveals subtle variations in expression levels between lipogenic and non-lipogenic conditions

Oleaginous yeasts have the ability to store greater than 20% of their mass as neutral lipids, in the form of triacylglycerides. The ATP citrate lyase is thought to play a key role in triacylglyceride synthesis, but the relationship between expression levels of this and other related enzymes is not well understood in the role of total lipid accumulation conferring the oleaginous phenotype. We conducted comparative proteomic analyses with the oleaginous yeast, Yarrowia lipolytica, grown in either nitrogen-sufficient rich media or nitrogen-limited minimal media. Total proteins extracted from cells collected during logarithmic and late stationary growth phases were analyzed by 1D liquid chromatography, followed by mass spectroscopy. The ATP citrate lyase enzyme was expressed at similar concentrations in both conditions, in both logarithmic and stationary phase, but many upstream and downstream enzymes showed drastically different expression levels. In non-lipogenic conditions, several pyruvate enzymes were expressed at higher concentration. These enzymes, especially the pyruvate decarboxylase and pyruvate dehydrogenase, may be regulating carbon flux away from central metabolism and reducing the amount of citrate being produced in the mitochondria. While crucial for the oleaginous phenotype, the constitutively expressed ATP citrate lyase appears to cleave citrate in response to carbon flux upstream from other enzymes creating the oleaginous phenotype.

Continental-Scale Gene Flow Prevents Allopatric Divergence of Pelagic Freshwater Bacteria

Allopatric divergence is one of the principal mechanisms for speciation of macro-organisms. Microbes by comparison are assumed to disperse more freely and to be less limited by dispersal barriers. However, thermophilic prokaryotes restricted to geothermal springs have shown clear signals of geographic isolation, but robust studies on this topic for microbes with less strict habitat requirements are scarce. Furthermore, it has only recently been recognized that homologous recombination among conspecific individuals provides species coherence in a wide range of prokaryotes. Recombination barriers thus may define prokaryotic species boundaries, yet, the extent to which geographic distance between populations gives rise to such barriers is an open question. Here, we investigated gene flow and population structure in a widespread species of pelagic freshwater bacteria, Polynucleobacter paneuropaeus. Through comparative genomics of 113 conspecific strains isolated from freshwater lakes and ponds located across a North–South range of more than 3,000 km, we were able to reconstruct past gene flow events. The species turned out to be highly recombinogenic as indicated by significant signs of gene transfer and extensive genome mosaicism. Although genomic differences increased with spatial distance on a regional scale (<170 km), such correlations were mostly absent on larger scales up to 3,400 km. We conclude that allopatric divergence in European P. paneuropaeus is minor, and that effective gene flow across the sampled geographic range in combination with a high recombination efficacy maintains species coherence.

Seeding Biochemistry on Other Worlds: Enceladus as a Case Study

The Solar System is becoming increasingly accessible to exploration by robotic missions to search for life. However, astrobiologists currently lack well-defined frameworks to quantitatively assess the chemical space accessible to life in these alien environments. Such frameworks will be critical for developing concrete predictions needed for future mission planning, both to determine the potential viability of life on other worlds and to anticipate the molecular biosignatures that life could produce. Here, we describe how uniting existing methods provides a framework to study the accessibility of biochemical space across diverse planetary environments. Our approach combines observational data from planetary missions with genomic data catalogued from across Earth and analyzed using computational methods from network theory. To demonstrate this, we use 307 biochemical networks generated from genomic data collected across Earth and "seed" these networks with molecules confirmed to be present on Saturn's moon Enceladus. By expanding through known biochemical reaction space starting from these seed compounds, we are able to determine which products of Earth's biochemistry are, in principle, reachable from compounds available in the environment on Enceladus, and how this varies across different examples of life from Earth (organisms, ecosystems, planetary-scale biochemistry). While we find that none of the 307 prokaryotes analyzed meet the threshold for viability, the reaction space covered by this process can provide a map of possible targets for detection of Earth-like life on Enceladus, as well as targets for synthetic biology approaches to seed life on Enceladus. In cases where biochemistry is not viable because key compounds are missing, we identify the environmental precursors required to make it viable, thus providing a set of compounds to prioritize for detection in future planetary exploration missions aimed at assessing the ability of Enceladus to sustain Earth-like life or directed panspermia.

Genomic Insights Into the Lifestyles of Thaumarchaeota Inside Sponges

Sponges are among the oldest metazoans and their success is partly due to their abundant and diverse microbial symbionts. They are one of the few animals that have Thaumarchaeota symbionts. Here we compare genomes of 11 Thaumarchaeota sponge symbionts, including three new genomes, to free-living ones. Like their free-living counterparts, sponge-associated Thaumarchaeota can oxidize ammonia, fix carbon, and produce several vitamins. Adaptions to life inside the sponge host include enrichment in transposases, toxin-antitoxin systems and restriction modifications systems, enrichments previously reported also from bacterial sponge symbionts. Most thaumarchaeal sponge symbionts lost the ability to synthesize rhamnose, which likely alters their cell surface and allows them to evade digestion by the host. All but one archaeal sponge symbiont encoded a high-affinity, branched-chain amino acid transporter system that was absent from the analyzed free-living thaumarchaeota suggesting a mixotrophic lifestyle for the sponge symbionts. Most of the other unique features found in sponge-associated Thaumarchaeota, were limited to only a few specific symbionts. These features included the presence of exopolyphosphatases and a glycine cleavage system found in the novel genomes. Thaumarchaeota have thus likely highly specific interactions with their sponge host, which is supported by the limited number of host sponge species to which each of these symbionts is restricted.

Trace gas oxidizers are widespread and active members of soil microbial communities

Soil microorganisms globally are thought to be sustained primarily by organic carbon sources. Certain bacteria also consume inorganic energy sources such as trace gases, but they are presumed to be rare community members, except within some oligotrophic soils. Here we combined metagenomic, biogeochemical and modelling approaches to determine how soil microbial communities meet energy and carbon needs. Analysis of 40 metagenomes and 757 derived genomes indicated that over 70% of soil bacterial taxa encode enzymes to consume inorganic energy sources. Bacteria from 19 phyla encoded enzymes to use the trace gases hydrogen and carbon monoxide as supplemental electron donors for aerobic respiration. In addition, we identified a fourth phylum (Gemmatimonadota) potentially capable of aerobic methanotrophy. Consistent with the metagenomic profiling, communities within soil profiles from diverse habitats rapidly oxidized hydrogen, carbon monoxide and to a lesser extent methane below atmospheric concentrations. Thermodynamic modelling indicated that the power generated by oxidation of these three gases is sufficient to meet the maintenance needs of the bacterial cells capable of consuming them. Diverse bacteria also encode enzymes to use trace gases as electron donors to support carbon fixation. Altogether, these findings indicate that trace gas oxidation confers a major selective advantage in soil ecosystems, where availability of preferred organic substrates limits microbial growth. The observation that inorganic energy sources may sustain most soil bacteria also has broad implications for understanding atmospheric chemistry and microbial biodiversity in a changing world.

Impacts of conductive materials on microbial community during syntrophic propionate oxidization for biomethane recovery

Propionate is one of the most important intermediates in anaerobic digestion, and its degradation requires a syntrophic partnership between propionate-oxidizing bacteria and hydrogenotrophic methanogens. Anaerobic digestion efficiency can be improved by direct interspecies electron transfer (DIET) through conductive materials. This study aimed to investigate the effects of DIET on syntrophic propionate oxidization under room temperature (20°C) and reveal the syntrophic partners. Firstly, conventional anaerobic consortium and conductive material-enriched consortium were tested for DIET under high H₂ partial pressure. The latter supplemented with granular activated carbon (GAC) can mitigate H₂ inhibition through DIET. Secondly, a DIET consortium was enriched for testing GAC and magnetite, both showed DIET facilitation. Microbial communities in GAC- and magnetite-supplemented reactors were similar. Syntrophic propionate-oxidizing bacteria, for example, Smithella (3.9%-9.9%) and a genus from the family Syntrophaceae (1.9%-3.6%) and methanogens Methanobacterium (30.3%-75.2%), Methanolinea (8.5%-25.2%), Methanosaeta (11.4%-36.7%), and Candidatus Methanofastidiosum (3.6%-6.6%), were predominant. Functional genes for cell mobility and membrane transport (3.3% and 9.5% in control reactor) increased with GAC (3.7% and 11.1%, respectively) and magnetite (3.7% and 10.9%, respectively) addition. Syntrophic propionate-oxidizing bacteria and methanogenesis partners were revealed by co-occurrence network, for example, Methanobacterium with Smithella, Syntrophobacter, Dechloromonas, and Trichococcus, signifying the importance of the syntrophic partnership in DIET environment. PRACTITIONER POINTS: DIET improved syntrophic propionate oxidization under room temperature condition (20°C). Microbial communities were similar for GAC- and magnetite-supplemented reactors, different with control reactor. Syntrophic propionate-oxidizing bacteria and methanogenesis partners were revealed by co-occurrence network. Methanobacterium and Smithella, Syntrophobacter, Dechloromonas, and Trichococcus were correlated.

Soil Acidobacteria Strain AB23 Resistance to Oxidative Stress Through Production of Carotenoids

Metagenomic studies revealed the prevalence of Acidobacteria in soils, but the physiological and ecological reasons for their success are not well understood. Many Acidobacteria exhibit carotenoid-related pigments, which may be involved in their tolerance of environmental stress. The aim of this work was to investigate the role of the orange pigments produced by Acidobacteria strain AB23 isolated from a savannah-like soil and to identify putative carotenoid genes in Acidobacteria genomes. Phylogenetic analysis revealed that strain AB23 belongs to the Occallatibacter genus from the class Acidobacteriia (subdivision 1). Strain AB23 produced carotenoids in the presence of light and vitamins; however, the growth rate and biomass decreased when cells were exposed to light. The presence of carotenoids resulted in tolerance to hydrogen peroxide. Comparative genomics revealed that all members of Acidobacteriia with available genomes possess the complete gene cluster for phytoene production. Some Acidobacteriia members have an additional gene cluster that may be involved in the production of colored carotenoids. Both colored and colorless carotenoids are involved in tolerance to oxidative stress. These results show that the presence of carotenoid genes is widespread among Acidobacteriia. Light and atmospheric oxygen stimulate carotenoid synthesis, but there are other natural sources of oxidative stress in soils. Tolerance to environmental oxidative stress provided by carotenoids may offer a competitive advantage for Acidobacteria in soils.

Metataxonomic profiling of bacterial communities and their predictive functional profiles in traditionally preserved meat products of Sikkim state in India

Traditionally preserved meat products are common food items in Sikkim state of India. We studied the high-throughput sequencing of four traditionally preserved meat products viz. beef kargyong, pork kargyong, yak satchu and khyopeh to profile the bacterial communities and also inferred their predictive functional profiles. Overall abundant OTUs in samples showed that Firmicutes was the abundant phylum followed by Proteobacteria and Bacteroidetes. Abundant species detected in each product were Psychrobacter pulmonis in beef kargyong, Lactobacillus sakei in pork kargyong, Bdellovibrio bacteriovorus and Ignatzschinera sp. in yak satchu and Lactobacillus sakei and Enterococcus sp. in khyopeh. Several genera unique to each product, based on analysis of shared OTUs contents, were observed among the samples except in khyopeh. Goods coverage recorded to 1.0 was observed, which reflected the maximum bacterial diversity in the samples. Alpha diversity metrics showed a maximum bacterial diversity in khyopeh and lowest in pork kargyong Community dissimilarities in the products were observed by PCoA plot. A total of 133 KEGG predictive functional pathways was observed in beef kargyong, 131 in pork kargyong, 125 in yak satchu and 101 in khyopeh. Metagenome contribution of the OTUs was computed using PICTRUSt2 and visualized by BURRITO software to predict the metabolic pathways. Several predictive functional profiles were contributed by abundant OTUs represented by Enterococcus, Acinetobacter, Agrobacterium, Bdellovibrio, Chryseobacterium, Lactococcus, Leuconostoc, Psychrobacter, and Staphylococcus.

Metagenomic insights into the metabolism and evolution of a new Thermoplasmata order (Candidatus Gimiplasmatales)

Thermoplasmata is a widely distributed and ecologically important archaeal class in the phylum Euryarchaeota. Because few cultures and genomes are available, uncharacterized Thermoplasmata metabolisms remain unexplored. In this study, we obtained four medium- to high-quality archaeal metagenome-assembled genomes (MAGs) from the filamentous fragments of black-odorous aquatic sediments (Foshan, Guangdong, China). Based on their 16S rRNA gene and ribosomal protein phylogenies, the four MAGs belong to the previously unnamed Thermoplasmata UBA10834 clade. We propose that this clade (five reference genomes from the Genome Taxonomy Database (GTDB) and four MAGs from this study) be considered a new order, Candidatus Gimiplasmatales. Metabolic pathway reconstructions indicated that the Ca. Gimiplasmatales MAGs can biosynthesize isoprenoids and nucleotides de novo. Additionally, some taxa have genes for formaldehyde and acetate assimilation, and the Wood-Ljungdahl CO₂ -fixation pathway, indicating a mixotrophic lifestyle. Sulfur reduction, hydrogen metabolism, and arsenic detoxification pathways were predicted, indicating sulfur-, hydrogen-, and arsenic-transformation potentials. Comparative genomics indicated that the H₄ F Wood-Ljungdahl pathway of both Ca. Gimiplasmatales and Methanomassiliicoccales was likely obtained by the interdomain lateral gene transfer from the Firmicutes. Collectively, this study elucidates the taxonomic and potential metabolic diversity of the new order Ca. Gimiplasmatales and the evolution of this subgroup and its sister lineage Methanomassiliicoccales.

Microbe-set enrichment analysis facilitates functional interpretation of microbiome profiling data

The commensal microbiome is known to influence a variety of host phenotypes. Microbiome profiling followed by differential abundance analysis has been established as an effective approach to study the mechanisms of host-microbiome interactions. However, it is challenging to interpret the collective functions of the resultant microbe-sets due to the lack of well-organized functional characterization of commensal microbiome. We developed microbe-set enrichment analysis (MSEA) to enable the functional interpretation of microbe-sets by examining the statistical significance of their overlaps with annotated groups of microbes that share common attributes such as biological function or phylogenetic similarity. We then constructed microbe-set libraries by query PubMed to find microbe-mammalian gene associations and disease associations by parsing the Disbiome database. To demonstrate the utility of our novel MSEA methodology, we carried out three case studies using publicly available curated knowledge resource and microbiome profiling datasets focusing on human diseases. We found MSEA not only yields consistent findings with the original studies, but also recovers insights about disease mechanisms that are supported by the literature. Overall, MSEA is a useful knowledge-based computational approach to interpret the functions of microbes, which can be integrated with microbiome profiling pipelines to help reveal the underlying mechanism of host-microbiome interactions.

Deciphering functional redundancy in the human microbiome

Although the taxonomic composition of the human microbiome varies tremendously across individuals, its gene composition or functional capacity is highly conserved - implying an ecological property known as functional redundancy. Such functional redundancy has been hypothesized to underlie the stability and resilience of the human microbiome, but this hypothesis has never been quantitatively tested. The origin of functional redundancy is still elusive. Here, we investigate the basis for functional redundancy in the human microbiome by analyzing its genomic content network - a bipartite graph that links microbes to the genes in their genomes. We find that this network exhibits several topological features that favor high functional redundancy. Furthermore, we develop a simple genome evolution model to generate genomic content network, finding that moderate selection pressure and high horizontal gene transfer rate are necessary to generate genomic content networks with key topological features that favor high functional redundancy. Finally, we analyze data from two published studies of fecal microbiota transplantation (FMT), finding that high functional redundancy of the recipient's pre-FMT microbiota raises barriers to donor microbiota engraftment. This work elucidates the potential ecological and evolutionary processes that create and maintain functional redundancy in the human microbiome and contribute to its resilience.

Comparisons of Transcriptome Profiles from Bacillus subtilis Cells Grown in Space versus High Aspect Ratio Vessel (HARV) Clinostats Reveal a Low Degree of Concordance

Although clinostats have long been used in space microbiology studies as ground-based analogs of spaceflight, few studies to date have systematically compared -omics data from clinostats versus spaceflight. This study compared the transcriptomic response of the Gram-positive bacterium Bacillus subtilis flown in space with corresponding transcriptomes derived from 2-D clinostat (High Aspect Ratio Vessel: HARV) experiments performed under the same conditions of bacterial strain, growth medium, temperature, and incubation time. High-quality total RNA (RNA Integrity Number >9.6) was isolated from multiple biological replicates from each treatment, transcripts were quantified by RNA-seq, and raw data was processed through a previously described standardized bioinformatics pipeline. Transcriptome data sets from spaceflight-grown and corresponding clinostat-grown cells were compared by using three different methods: (i) principal component analysis, (ii) analysis of differentially expressed genes, and (iii) gene set enrichment analysis of KEGG pathways. All three analyses found a low degree of concordance between the spaceflight and corresponding clinostat transcriptome data sets, ranging from 0.9% to 5.3% concordance. These results are in agreement with prior studies that also revealed low concordances between spaceflight and clinostat transcriptomes of the Gram-negative bacteria Rhodospirillum rubrum and Pseudomonas aeruginosa. The results are discussed from the perspective of several potential confounding factors, and suggestions are offered with the aim of achieving increased concordance between clinostat and spaceflight data.

Bacterial ecosystem functioning in organic matter biodegradation of different composting at the thermophilic phase

This study aimed to provide insights into prediction of composting ecological functioning through analyzing the critical bacterial populations and functions. The bacterial ecosystem functioning was essential, and cow dung, chicken manure, mushroom dreg and Chinese medicine residues were used as raw materials to quantify and predict the functioning of bacterial communities through synthetic spike-in standards accompanied Illumina sequencing and PICRUSt. Bacterial community of wheat straw and chicken manure compost (SCM) was similar to mushroom dreg and chicken manure compost (MCM), and Sinibacillus dominated in both treatments with the abundance of 20.73% and 41.36%, respectively. The correlation analysis between bacterial community and fluorescence EEM regional integration parameters showed that Lactobacillus (0.889), Enterococcus (0.888) and Erysipelothri (0.903) were positively correlated with P_V, n / P_III, n. The ontology analysis results showed that metabolism, genetic information processing, environmental information processing and cellular processes were the primary functions for bacterial community in all treatments.

Dysbiotic Lesional Microbiome With Filaggrin Missense Variants Associate With Atopic Dermatitis in India

Background: Atopic Dermatitis (AD) has been associated with the loss of function (LoF) mutations in Filaggrin (FLG) gene and increase in relative abundance of specific microbes in the lesional skin, predominantly in Caucasians. Our study aims to determine, in Indian AD patients, (a) the prevalence of FLG LoF and missense mutations, and (b) the nature and extent of dysbiosis and altered microbial pathways with and without mutations in FLG. AD patients (n = 34) and healthy controls (n = 54) were recruited from India in this study and shotgun sequencing was carried out in a subset of samples with adequate microbiome DNA concentration. Host DNA from the same subset of samples was subjected to FLG coding region sequencing and host-microbiome association was estimated. Results: The prevalence of FLG LoFs that are associated with AD globally were significantly lesser in our cases and controls (8.6%, 0%) than those reported in Europeans (27%, 2.6%). Staphylococcus aureus was present only on AD skin [abundance in Pediatric AD: 32.86%; Adult AD: 22.17%], but not on healthy skin on which Staphylococcus hominis (Adult controls: 16.43%, Adult AD: 0.20%; p = 0.002), Cutibacterium acnes (Adult controls:10.84%, Adult AD: 0.90%; p = 0.02), and Malassezia globosa (Adult controls: 8.89%, Adult AD: 0.005%; p = 0.001) were significantly more abundant. Microbial pathways mostly associated with skin barrier permeability, ammonia production and inflammation (Arginine and Proline metabolism, Histidine Metabolism and Staphylococcus aureus infection) were significantly enriched on AD skin metagenome. These pathways are also reported to impair antimicrobial peptide activity. Among AD patients with missense single nucleotide polymorphisms harboring "potentially damaging" alleles in FLG gene, damaging allele dosage was significantly (p < 0.02) positively correlated with relative abundance of phylum_Proteobacteria up to order_Pseudomonadales and negatively correlated with phylum_Firmicutes up to species_Staphylococcus aureus. Conclusion: Our study has provided evidence that host DNA profile is significantly associated with microbiome composition in the development of AD. Species and strain level analysis showed that the microbial pathways enriched in AD cases were mostly found in MRSA strains. These evidences can be harnessed to control AD by modulating the microbiome using a personalized strategy. Our findings on the association of FLG genotypes with the microbiome dysbiosis may pave the way for a personalized strategy to provide a more effective control of AD.

Blackwater biomethane recovery using a thermophilic upflow anaerobic sludge blanket reactor: Impacts of effluent recirculation on reactor performance

Thermophilic anaerobic digestion is a promising process for high-solid blackwater (BW) treatment due to improved hydrolysis rates, high methanogenesis efficiency, and pathogen removal, when compared with mesophilic treatment. In the present work, the effects of effluent recirculation (i.e., mixing) on thermophilic blackwater treatment were studied. A laboratory-scale thermophilic upflow anaerobic sludge blanket reactor was operated with and without effluent recirculation. The methanogenesis efficiency of the BW treatment increased from 45.0 ± 2.9% when effluent recirculation was applied to 56.7 ± 5.5% without effluent recirculation. Without effluent recirculation, the COD accumulation in the bioreactor was reduced from 17.2 to 3.8% and the effluent volatile fatty acids (VFA) concentration was reduced from 0.64 ± 0.18 to 0.15 ± 0.10 g/L. Further, both acetoclastic and hydrogenotrophic methanogenic activity increased from 101.3 ± 10.8 and 93.9 ± 6.1 to 120.4 ± 9.4 and 118.2 ± 13.2 mg CH₄-COD/(gVSS⋅d), respectively, after effluent recirculation was discontinued. The predominant methanogens changed from Methanothermobacter (67%) with effluent recirculation to Methanosarcina (62%) without effluent recirculation. As compared to the effluent recirculation conditions, the enhanced biomethane recovery and treatment performance without effluent recirculation can be attributed to the close proximity of bacteria and archaea groups and the reduced VFA accumulation. Predicted functional gene comparison showed higher prevalence of function for intermediate metabolite transportation (transporters, ATP-binding cassette (ABC) transporters, and two-component system) after discontinuing effluent circulation, which contributed to improved syntrophic propionate oxidation and syntrophic acetate oxidization and enhanced hydrogenotrophic methanogenesis.

METATRYP v 2.0: Metaproteomic Least Common Ancestor Analysis for Taxonomic Inference Using Specialized Sequence Assemblies-Standalone Software and Web Servers for Marine Microorganisms and Coronaviruses

We present METATRYP version 2 software that identifies shared peptides across the predicted proteomes of organisms within environmental metaproteomics studies to enable accurate taxonomic attribution of peptides during protein inference. Improvements include ingestion of complex sequence assembly data categories (metagenomic and metatranscriptomic assemblies, single cell amplified genomes, and metagenome assembled genomes), prediction of the least common ancestor (LCA) for a peptide shared across multiple organisms, increased performance through updates to the backend architecture, and development of a web portal (https://metatryp.whoi.edu). Major expansion of the marine METATRYP database with predicted proteomes from environmental sequencing confirms a low occurrence of shared tryptic peptides among disparate marine microorganisms, implying tractability for targeted metaproteomics. METATRYP was designed to facilitate ocean metaproteomics and has been integrated into the Ocean Protein Portal (https://oceanproteinportal.org); however, it can be readily applied to other domains. We describe the rapid deployment of a coronavirus-specific web portal (https://metatryp-coronavirus.whoi.edu/) to aid in use of proteomics on coronavirus research during the ongoing pandemic. A coronavirus-focused METATRYP database identified potential SARS-CoV-2 peptide biomarkers and indicated very few shared tryptic peptides between SARS-CoV-2 and other disparate taxa analyzed, sharing <1% peptides with taxa outside of the betacoronavirus group, establishing that taxonomic specificity is achievable using tryptic peptide-based proteomic diagnostic approaches.