COG database update: focus on microbial diversity, model organisms, and widespread pathogens

Epiphytic common core bacteria in the microbiomes of co-located green (Ulva), brown (Saccharina) and red (Grateloupia, Gelidium) macroalgae

BACKGROUND

Macroalgal epiphytic microbial communities constitute a rich resource for novel enzymes and compounds, but studies so far largely focused on tag-based microbial diversity analyses or limited metagenome sequencing of single macroalgal species.

RESULTS

We sampled epiphytic bacteria from specimens of Ulva sp. (green algae), Saccharina sp. (brown algae), Grateloupia sp. and Gelidium sp. (both red algae) together with seawater and sediment controls from a coastal reef in Weihai, China, during all seasons. Using 16S rRNA amplicon sequencing, we identified 14 core genera (consistently present on all macroalgae), and 14 dominant genera (consistently present on three of the macroalgae). Core genera represented ~ 0.7% of all genera, yet accounted for on average 51.1% of the bacterial abundances. Plate cultivation from all samples yielded 5,527 strains (macroalgae: 4,426) representing 1,235 species (685 potentially novel). Sequencing of selected strains yielded 820 non-redundant draft genomes (506 potentially novel), and sequencing of 23 sampled metagenomes yielded 1,619 metagenome-assembled genomes (MAGs), representing further 1,183 non-redundant genomes. 230 isolates and 153 genomes were obtained from the 28 core/dominant genera. We analyzed the genomic potential of phycosphere bacteria to degrade algal polysaccharides and to produce bioactive secondary metabolites. We predicted 4,451 polysaccharide utilization loci (PULs) and 8,810 biosynthetic gene clusters (BGCs). These were particularly prevalent in core/dominant genera.

CONCLUSIONS

Our metabolic annotations and analyses of MAGs and genomes provide new insights into novel species of phycosphere bacteria and their ecological niches for an improved understanding of the macroalgal phycosphere microbiome. Video Abstract.

Cataloging the Phylogenetic Diversity of Human Bladder Bacterial Isolates

Although the human bladder is reported to harbor unique microbiota, our understanding of how these microbial communities interact with their human hosts is limited, mostly owing to the lack of isolates to test mechanistic hypotheses. Niche-specific bacterial collections and associated reference genome databases have been instrumental in expanding knowledge of the microbiota of other anatomical sites, e.g., the gut and oral cavity. To facilitate genomic, functional, and experimental analyses of the human bladder microbiota, here we present a bladder-specific bacterial reference collection comprised of 1134 genomes. These genomes were culled from bacterial isolates obtained by a metaculturomic method from bladder urine collected by transurethral catheterization. This bladder-specific bacterial reference collection includes 196 different species, including representatives of major aerobes and facultative anaerobes, as well as some anaerobes. It captures 72.2 % of the genera found when we reexamined previously published 16S rRNA gene sequencing of 392 adult female bladder urine samples. Comparative genomic analysis found that the taxonomies and functions of the bladder microbiota shared more similarities with the vaginal microbiota than the gut microbiota. Whole-genome phylogenetic and functional analyses of 186 bladder E. coli isolates and 387 gut E. coli isolates supports the hypothesis that phylogroup distribution and functions of E. coli strains differ dramatically between these two very different niches. This bladder-specific bacterial reference collection is a unique resource that will enable hypothesis-driven bladder microbiota research and comparison to isolates from other anatomical sites.

Genome features and secondary metabolite potential of the marine symbiont Streptomyces sp. RS2

Streptomyces sp. RS2 was isolated from an unidentified sponge collected around Randayan Island, Indonesia. The genome of Streptomyces sp. RS2 consists of a linear chromosome of 9,391,717 base pairs with 71.9% of G + C content, 8270 protein-coding genes, as well as 18 rRNA and 85 tRNA loci. Twenty-eight putative secondary metabolites biosynthetic gene clusters (BGCs) were identified in the genome sequence. Nine of them have 100% similarity to BGCs for albaflavenone, α-lipomycin, coelibactin, coelichelin, ectoine, geosmin, germicidin, hopene, and lanthionine (SapB). The remaining 19 BGCs have low (< 50%) or moderate (50-80%) similarity to other known secondary metabolite BGCs. Biological activity assays of extracts from 21 different cultures of the RS2 strain showed that SCB ASW was the best medium for the production of antimicrobial and cytotoxic compounds. Streptomyces sp. RS2 has great potential to be a producer of novel secondary metabolites, particularly those with antimicrobial and antitumor activities.

Identification and genomic analysis of Pseudosulfitobacter koreense sp. nov. isolated from toxin-producing dinoflagellate Alexandrium pacificum

The bacterial strain AP-MA-4^T isolated from the marine dinoflagellate Alexandrium pacificum (KCTC AG60911), was subjected to a taxonomic analysis. Cells of strain AP-MA-4^T were Gram-stain-negative, aerobic, rod-shaped, optimum growth at 20 °C, pH 7.0, in the presence of 5% (w/v) NaCl. Strain AP-MA-4^T shared the highest 16S rRNA gene sequence similarity to Pseudosulfitobacter pseudonitzschiae DSM 26824^T (98.5%), followed by Ascidiaceihabitans donghaensis RSS1-M3^T (96.3%), Pseudoseohaeicola caenipelagi BS-W13^T (95.7%), and Sulfitobacter pontiacus CHLG 10^T (95.3%). Based on 16S rRNA phylogeny, strain AP-MA-4^T is phylogenetically closely related to Pseudosulfitobacter pseudonitzschiae (type species of Pseudosulfitobacter) and could be distinguished from the type species based on their phenotypic properties. The genome length of strain AP-MA-4^T was 3.48 Mbp with a 62.9% G + C content. The average nucleotide identity (ANI) and digital DNA-DNA hybridization (dDDH) values between strain AP-MA-4 ^T and its closely related type strains were 72.2-83.3 and 18.2-27.6%, respectively. Summed feature 8 (C_18:1ω7c and/or C_18:1ω6c) was identified the major fatty acids (> 10%). Phosphatidylglycerol (PG), phosphatidylethanolamine (PE), and phospholipid (PL) were demonstrated as the major polar lipids. The major respiratory quinone is ubiquinone-10 (Q-10). Based on genotypic and phenotypic features, strain AP-MA-4^T (= KCTC 92289^T = GDMCC 1.3585^T) represents a new Pseudosulfitobacter species, in which the name Pseudosulfitobacter koreense sp. nov. is proposed.

Characterization of the microbiome and volatile compounds in anal gland secretions from domestic cats (Felis catus) using metagenomics and metabolomics

Animals rely on volatile chemical compounds for their communication and behavior. Many of these compounds are sequestered in endocrine and exocrine glands and are synthesized by anaerobic microbes. While the volatile organic compound (VOC) or microbiome composition of glandular secretions has been investigated in several mammalian species, few have linked specific bacterial taxa to the production of volatiles or to specific microbial gene pathways. Here, we use metagenomic sequencing, mass-spectrometry based metabolomics, and culturing to profile the microbial and volatile chemical constituents of anal gland secretions in twenty-three domestic cats (Felis catus), in attempts to identify organisms potentially involved in host odor production. We found that the anal gland microbiome was dominated by bacteria in the genera Corynebacterium, Bacteroides, Proteus, Lactobacillus, and Streptococcus, and showed striking variation among individual cats. Microbiome profiles also varied with host age and obesity. Metabolites such as fatty-acids, ketones, aldehydes and alcohols were detected in glandular secretions. Overall, microbiome and metabolome profiles were modestly correlated (r=0.17), indicating that a relationship exists between the bacteria in the gland and the metabolites produced in the gland. Functional analyses revealed the presence of genes predicted to code for enzymes involved in VOC metabolism such as dehydrogenases, reductases, and decarboxylases. From metagenomic data, we generated 85 high-quality metagenome assembled genomes (MAGs). Of these, four were inferred to have high relative abundance in metagenome profiles and had close relatives that were recovered as cultured isolates. These four MAGs were classified as Corynebacterium frankenforstense, Proteus mirabilis, Lactobacillus johnsonii, and Bacteroides fragilis. They represent strong candidates for further investigation of the mechanisms of volatile synthesis and scent production in the mammalian anal gland.

Genomic characterization of coexisting anatoxin-producing and non-toxigenic Microcoleus subspecies in benthic mats from the Wolastoq, New Brunswick, Canada

The presence of toxigenic benthic cyanobacteria in riverine ecosystems is an increasing concern around the world. In 2018, the death of three dogs along the Wolastoq (also known as the Saint John River) in New Brunswick, Canada, was attributed to anatoxin exposure after they ingested benthic microbial mats found along the shore. Here, we shotgun sequenced the DNA of 15 non-axenic cyanobacterial isolates derived from four anatoxin-containing benthic mat samples associated with the dog deaths. Anatoxins were produced by some of the isolates, but not all. We retrieved near-complete Microcoleus metagenome-assembled genomes (MAGs) from the isolates that are closely related to anatoxin-producing Microcoleus from the Cardrona River (New Zealand), although the Microcoleus MAGs from the Wolastoq varied in the presence/absence of the anatoxin-a biosynthesis cluster. Sequence similarity at the genomic level suggests that toxigenic and non-toxigenic Microcoleus MAGs from the Wolastoq belong to the same species but are separate subspecies. The toxigenic and nontoxic Wolastoq Microcoleus subspecies coexisted in the mat samples in similar relative abundance. Overall genomic comparisons revealed that toxigenic Microcoleus MAGs are longer and code for more accessory genes than their non-toxigenic relatives, suggesting a differential responsiveness to changing environments, stress conditions and nutrient availability.

Comparative genomic analysis provides insights into taxonomy and temperature adaption of Aeromonas salmonicida

Aeromonas salmonicida has long been known as psychrophiles since it is mainly isolated from cold water fish, and recent reports have revealed the existence of mesophilic strains isolated from warm sources. However, the genetic differences between mesophilic and psychrophilic strains remain unclear due to few complete genomes of mesophilic strain are available. In this study, six A. salmonicida (2 mesophilic and 4 psychrophilic) were genome-sequenced, and comparative analyses of 25 A. salmonicida complete genomes were conducted. The ANI values and phylogenetic analysis revealed that 25 strains formed three independent clades, which were referred as typical psychrophilic, atypical psychrophilic and mesophilic groups. Comparative genomic analysis showed that two chromosomal gene clusters, related to lateral flagella and outer membrane proteins (A-layer and T2SS proteins), and insertion sequences (ISAs4, ISAs7 and ISAs29) were unique to the psychrophilic groups, while the complete MSH type IV pili were unique to the mesophilic group, all of which may be considered as lifestyle-related factors. The results of this study not only provide new insights into the classification, lifestyle adaption and pathogenic mechanism of different strains of A. salmonicida, but also contributes to the prevention and control of disease caused by psychrophilic and mesophilic A. salmonicida.

Metabolic cooperation between conspecific genotypic groups contributes to bacterial fitness

Microbial interactions are important for the survival of species and the stability of the microbial ecosystem. Although bacteria have diverse conspecific genotypes in the natural microbial ecosystem, little is known about whether wild-type strains within species would interact with each other and how the intraspecific interaction influences the growth of the species. In this work, using Lactobacillus acetotolerans, a dominant species with diverse conspecific genotypes in natural food fermentation ecosystems as a case, we studied the interactions between different genotypic groups of this species. In interspecific and intraspecific pairwise cocultures, the growth of L. acetotolerans decreased, but the increase of the phylogenetic similarity would reduce this negative effect, indicating a potential intraspecific interaction of this species. Meanwhile, the strain classification method affected the analysis of intraspecific interactions, which can be efficiently demonstrated using 99.5% average nucleotide identity (ANI) as the strain-level classification method. Using this ANI classification method, we revealed the population fitness significantly increased in cocultures of different genotypic groups. Facilitation involving 11 amino acids was identified between different ANI genotypic groups, which was beneficial for increasing population fitness. This work revealed that wild-type conspecific strains could interact with each other via cooperative metabolic changes and benefit each other to increase fitness. It shed new light on the survival and stability of species in natural microbial ecosystems.

G. de Lomana A, Gomes-Filho JV, Vêncio RZN, Moritz RL, Koide T, Baliga NS. 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).

The nickel-pincer coenzyme of lactate racemase: A case study of uncovering cofactor structure and biosynthesis

Cofactors are essential components of numerous enzymes, therefore their characterization by structural, biophysical, and biochemical approaches is crucial for understanding the resulting catalytic and regulatory mechanisms. In this chapter, we present a case study of a recently discovered cofactor, the nickel-pincer nucleotide (NPN), by demonstrating how we identified and thoroughly characterized this unprecedented nickel-containing coenzyme that is tethered to lactase racemase from Lactiplantibacillus plantarum. In addition, we describe how the NPN cofactor is biosynthesized by a panel of proteins encoded in the lar operon and describe the properties of these novel enzymes. Comprehensive protocols for conducting functional and mechanistic studies of NPN-containing lactate racemase (LarA) and the carboxylase/hydrolase (LarB), sulfur transferase (LarE), and metal insertase (LarC) used for NPN biosynthesis are provided for potential applications towards characterizing enzymes in the same or homologous families.

Whole-genome and comparative genome analysis of Mucor racemosus C isolated from Yongchuan Douchi

Mucor racemosus is the predominant fungal in the zhiqu stage of the fermentation of Yongchuan Douchi (Mucor-type), which plays an important role in the fermentation process of Yongchuan Douchi. However, there is a lack of information on the genetic analysis of M. racemosus. In this study, we isolated and identified M. racemosus C (accession no JAPEHQ000000000) from Yongchuan Douchi and analyzed the physiological indicators, then genomic information of the strain to perform a comprehensive analysis of its fermentation capacity and safety. M. racemosus C had neutral protease activity up to 68.051 U/mL at 30 °C and alkaline protease activity up to 57.367 U/mL at 25 °C. In addition, comparing the genomic data with the COGs database (NCBI), it was predicted that M. racemosus C undergoes extensive amino acid metabolism, making C suitable for the production of fermented foods (e.g., Douchi, Syoyu, and sufu). Finally, we performed virulence genes and resistance genes analysis, hemolysis experiment, aflatoxins assay, antibiotic resistance assay to evaluate the safety of M. racemosus C, and the results showed that M. racemosus C was safe, non-toxin-producing and non-hemolytic.

Transcriptome Analysis of Resistant and Susceptible Pecan (Carya illinoinensis) Reveals the Mechanism of Resistance to Black Spot Disease (Colletotrichum fioriniae)

Pecan, Carya illinoinensis (Wangenh.) K. Koch, is an important dried fruit and woody oil tree species grown worldwide. With continuous expansion of pecan cultivation, the frequency and scope of diseases, especially black spot disease, are increasing, damaging trees and reducing yields. In this study, the key factors in resistance to black spot disease (Colletotrichum fioriniae) were investigated between the high-resistance pecan variety "Kanza" and the low-resistance variety "Mahan". Leaf anatomy and antioxidase activities confirmed much stronger resistance to black spot disease in "Kanza" than in "Mahan". Transcriptome analysis indicated that the increased expression of genes associated with defense response, oxidation-reduction, and catalytic activity was involved in disease resistance. A connection network identified a highly expressed hub gene CiFSD2 (CIL1242S0042), which might participate in redox reactions to affect disease resistance. Overexpression of CiFSD2 in tobacco inhibited enlargement of necrotic spots and increased disease resistance. Overall, the expression of differentially expressed genes differed in pecan varieties with different levels of resistance to C. fioriniae infection. In addition, the hub genes associated with black spot resistance were identified and the functions clarified. The in-depth understanding of resistance to black spot disease provides new insights for early screening of resistant varieties and molecular-assisted breeding in pecan.

Dasania phycosphaerae sp. nov., isolated from phytoplankton sample from the south coast of the Republic of Korea

A novel bacterial strain, designated as GY-18T, was isolated from particulate material sampled from the Korean south coast. The bacterium was Gram-stain-negative, aerobic and rod-shaped. Oxidase and catalase were positive. Optimal growth was observed at 30 °C, pH 7.0, in the presence of 1 % (w/v) NaCl. Phylogenetic analyses based on 16S rRNA gene and whole-genome sequences suggested that strain GY-18T was a member of the genus Dasania , but represented a separate lineage. The 16S rRNA gene sequence of strain GY-18T showed high similarities to Dasania marina KCTC 12566T (98.4 %), Halioglobus japonicus KCTC 23429T (92.2 %), Marimicrobium arenosum KCTC 42300T (92.4 %) and Seongchinamella unica KCTC 62383T (92.2 %), respectively. The genome length of strain GY-18T was 3 808 450 bp with a DNA G+C content of 47.8 mol %. The major cellular fatty acids of strain GY-18T were summed feature 3 (comprising C16 : 1  ω7c and/or C16 : 1  ω6c), C16 : 0 and C10 : 0 3-OH (>10 %). Diphosphatidylglycerol, phosphatidylglycerol, phosphatidylethanolamine and phospholipids were shown to be the major polar lipids. The predominant isoprenoid quinone was Q-8. Based on phylogenetic, phenotypic, chemotaxonomic and genomic features, we propose that strain GY-18T represents a novel species in the genus Dasania of the family Spongiibacteraceae , with the proposed name Dasania phycosphaerae sp. nov. The type strain is GY-18T (=KCTC 92290T=GDMCC 1.3586T).

Using a whole genome co-expression network to inform the functional characterisation of predicted genomic elements from Mycobacterium tuberculosis transcriptomic data

A whole genome co-expression network was created using Mycobacterium tuberculosis transcriptomic data from publicly available RNA-sequencing experiments covering a wide variety of experimental conditions. The network includes expressed regions with no formal annotation, including putative short RNAs and untranslated regions of expressed transcripts, along with the protein-coding genes. These unannotated expressed transcripts were among the best-connected members of the module sub-networks, making up more than half of the 'hub' elements in modules that include protein-coding genes known to be part of regulatory systems involved in stress response and host adaptation. This data set provides a valuable resource for investigating the role of non-coding RNA, and conserved hypothetical proteins, in transcriptomic remodelling. Based on their connections to genes with known functional groupings and correlations with replicated host conditions, predicted expressed transcripts can be screened as suitable candidates for further experimental validation.

The evolution of morphological development is congruent with the species phylogeny in the genus Streptomyces

As the canonical model organism to dissect bacterial morphological development, Streptomyces species has attracted much attention from the microbiological society. However, the evolution of development-related genes in Streptomyces remains elusive. Here, we evaluated the distribution of development-related genes, thus indicating that the majority of these genes were ubiquitous in Streptomyces genomes. Furthermore, the phylogenetic topologies of related strict orthologous genes were compared to the species tree of Streptomyces from both concatenation and single-gene tree analyses. Meanwhile, the reconciled gene tree and normalization based on the number of parsimony-informative sites were also employed to reduce the impact of phylogenetic conflicts, which was induced by uncertainty in single-gene tree inference based merely on the sequence and the bias in the amount of phylogenetic information caused by variable numbers of parsimony-informative sites. We found that the development-related genes had higher congruence to the species tree than other strict orthologous genes. Considering that the development-related genes could also be tracked back to the common ancestor of Streptomyces, these results suggest that morphological development follows the same pattern as species divergence.

The phylogeny and metabolic potentials of an n-alkane-degrading Venatorbacter bacterium isolated from deep-sea sediment of the Mariana Trench

Recently, several reports showed that n-alkanes were abundant in the hadal zone, suggesting that n-alkanes could be an important source of nutrients for microorganisms in hadal ecosystems. To date, most of the published studies on the microbial capacity to degrade hydrocarbons were conducted only at atmospheric temperature and pressure (0.1 MPa), and little is known about whether and which microbes could utilize n-alkanes at in situ environmental conditions in the hadal zone, including low temperature and high hydrostatic pressure (especially >30 MPa). In this study, a piezotolerant bacterium, strain C2-1, was isolated from a Mariana Trench sediment at depth of 5,800 m. Strain C2-1 was able to grow at in situ temperature (4°C) and pressure (58 MPa) with n-alkanes as the sole carbon source. Phylogenetically, strain C2-1 and related strains (TMPB967, ST750PaO-4, IMCC1826, and TTBP476) should be classified into the genus Venatorbacter. Metagenomic analysis using ~5,000 publicly available datasets showed that Venatorbacter has a wide environmental distribution in seawater (38), marine sediments (3), hydrothermal vent plumes (2), Antarctic ice (1), groundwater (13), and marine sponge ecosystems (1). Most Venatorbacter species are non-obligate n-alkane degraders that could utilize, at a minimal, C16-C18 n-alkanes, as well as other different types of carbon substrates, including carbohydrates, amino acids, peptides, and phospholipids. The type II secretion system, extracellular proteases, phospholipase, and endonuclease of Venatorbacter species were robustly expressed in the metatranscriptomes of deep-sea hydrothermal vents, suggesting their important contribution to secondary productivity by degrading extracellular macromolecules. The identification of denitrifying genes suggested a genus-specific ecological potential that allowed Venatorbacter species to be active in anoxic environments, e.g., the oxygen-minimal zone (OMZ) and the deeply buried marine sediments. Our results show that Venatorbacter species are responsible for the degradation of hydrocarbon and extracellular macromolecules, suggesting that they may play an important role in the biogeochemistry process in the Trench ecosystems.

Defining the Infant Male Urobiome and Moving Towards Mechanisms in Urobiome Research

The urinary bladder harbors a community of microbes termed the urobiome, which remains understudied. In this study, we present the urobiome of healthy infant males from samples collected by transurethral catheterization. Using a combination of extended culture and amplicon sequencing, we identify several common bacterial genera that can be further investigated for their effects on urinary health across the lifespan. Many genera were shared between all samples suggesting a consistent urobiome composition among this cohort. We note that, for this cohort, early life exposures including mode of birth (vaginal vs. Caesarean section), or prior antibiotic exposure did not influence urobiome composition. In addition, we report the isolation of culturable bacteria from the bladders of these infant males, including Actinotignum schaalii, a bacterial species that has been associated with urinary tract infection in older male adults. Herein, we isolate and sequence 9 distinct strains of A. schaalii enhancing the genomic knowledge surrounding this species and opening avenues for delineating the microbiology of this urobiome constituent. Furthermore, we present a framework for using the combination of culture-dependent and sequencing methodologies for uncovering mechanisms in the urobiome.

Understudied proteins and understudied functions in the model bacterium Bacillus subtilis-A major challenge in current research

Model organisms such as the Gram-positive bacterium Bacillus subtilis have been studied intensively for decades. However, even for model organisms, no function has been identified for about one fourth of all proteins. It has recently been realized that such understudied proteins as well as poorly studied functions set a limitation to our understanding of the requirements for cellular life, and the Understudied Proteins Initiative has been launched. Of poorly studied proteins, those that are strongly expressed are likely to be important to the cell and should therefore be considered high priority in further studies. Since the functional analysis of unknown proteins can be extremely laborious, a minimal knowledge is required prior to targeted functional studies. In this review, we discuss strategies to obtain such a minimal annotation, for example, from global interaction, expression, or localization studies. We present a set of 41 highly expressed and poorly studied proteins of B. subtilis. Several of these proteins are thought or known to bind RNA and/or the ribosome, some may control the metabolism of B. subtilis, and another subset of particularly small proteins may act as regulatory elements to control the expression of downstream genes. Moreover, we discuss the challenges of poorly studied functions with a focus on RNA-binding proteins, amino acid transport, and the control of metabolic homeostasis. The identification of the functions of the selected proteins not only will strongly advance our knowledge on B. subtilis, but also on other organisms since many of the proteins are conserved in many groups of bacteria.

Exploring microbial functional biodiversity at the protein family level-From metagenomic sequence reads to annotated protein clusters

Metagenomics has enabled accessing the genetic repertoire of natural microbial communities. Metagenome shotgun sequencing has become the method of choice for studying and classifying microorganisms from various environments. To this end, several methods have been developed to process and analyze the sequence data from raw reads to end-products such as predicted protein sequences or families. In this article, we provide a thorough review to simplify such processes and discuss the alternative methodologies that can be followed in order to explore biodiversity at the protein family level. We provide details for analysis tools and we comment on their scalability as well as their advantages and disadvantages. Finally, we report the available data repositories and recommend various approaches for protein family annotation related to phylogenetic distribution, structure prediction and metadata enrichment.

Mixed heavy metal stress induces global iron starvation response

Multiple heavy metal contamination is an increasingly common global problem. Heavy metals have the potential to disrupt microbially mediated biogeochemical cycling. However, systems-level studies on the effects of combinations of heavy metals on bacteria are lacking. For this study, we focused on the Oak Ridge Reservation (ORR; Oak Ridge, TN, USA) subsurface which is contaminated with several heavy metals and high concentrations of nitrate. Using a native Bacillus cereus isolate that represents a dominant species at this site, we assessed the combined impact of eight metal contaminants, all at site-relevant concentrations, on cell processes through an integrated multi-omics approach that included discovery proteomics, targeted metabolomics, and targeted gene-expression profiling. The combination of eight metals impacted cell physiology in a manner that could not have been predicted from summing phenotypic responses to the individual metals. Exposure to the metal mixture elicited a global iron starvation response not observed during individual metal exposures. This disruption of iron homeostasis resulted in decreased activity of the iron-cofactor-containing nitrate and nitrite reductases, both of which are important in biological nitrate removal at the site. We propose that the combinatorial effects of simultaneous exposure to multiple heavy metals is an underappreciated yet significant form of cell stress in the environment with the potential to disrupt global nutrient cycles and to impede bioremediation efforts at mixed waste sites. Our work underscores the need to shift from single- to multi-metal studies for assessing and predicting the impacts of complex contaminants on microbial systems.

Botryosphaeriaceae gene machinery: Correlation between diversity and virulence

The Botryosphaeriaceae family comprises numerous fungal pathogens capable of causing economically meaningful diseases in a wide range of crops. Many of its members can live as endophytes and turn into aggressive pathogens following the onset of environmental stress events. Their ability to cause disease may rely on the production of a broad set of effectors, such as cell wall-degrading enzymes, secondary metabolites, and peptidases. Here, we conducted comparative analyses of 41 genomes representing six Botryosphaeriaceae genera to provide insights into the genetic features linked to pathogenicity and virulence. We show that these Botryosphaeriaceae genomes possess a large diversity of carbohydrate-active enzymes (CAZymes; 128 families) and peptidases (45 families). Botryosphaeria, Neofusicoccum, and Lasiodiplodia presented the highest number of genes encoding CAZymes involved in the degradation of the plant cell wall components. The genus Botryosphaeria also exhibited the highest abundance of secreted CAZymes and peptidases. Generally, the secondary metabolites gene cluster profile was consistent in the Botryosphaeriaceae family, except for Diplodia and Neoscytalidium. At the strain level, Neofusicoccum parvum NpBt67 stood out among all the Botryosphaeriaceae genomes, presenting a higher number of secretome constituents. In contrast, the Diplodia strains showed the lowest richness of the pathogenicity- and virulence-related genes, which may correlate with their low virulence reported in previous studies. Overall, these results contribute to a better understanding of the mechanisms underlying pathogenicity and virulence in remarkable Botryosphaeriaceae species. Our results also support that Botryosphaeriaceae species could be used as an interesting biotechnological tool for lignocellulose fractionation and bioeconomy.

Plasmid-Encoded Traits Vary across Environments

Plasmids are key mobile genetic elements in bacterial evolution and ecology as they allow the rapid adaptation of bacteria under selective environmental changes. However, the genetic information associated with plasmids is usually considered separately from information about their environmental origin. To broadly understand what kinds of traits may become mobilized by plasmids in different environments, we analyzed the properties and accessory traits of 9,725 unique plasmid sequences from a publicly available database with known bacterial hosts and isolation sources. Although most plasmid research focuses on resistance traits, such genes made up <1% of the total genetic information carried by plasmids. Similar to traits encoded on the bacterial chromosome, plasmid accessory trait compositions (including general Clusters of Orthologous Genes [COG] functions, resistance genes, and carbon and nitrogen genes) varied across seven broadly defined environment types (human, animal, wastewater, plant, soil, marine, and freshwater). Despite their potential for horizontal gene transfer, plasmid traits strongly varied with their host's taxonomic assignment. However, the trait differences across environments of broad COG categories could not be entirely explained by plasmid host taxonomy, suggesting that environmental selection acts on the plasmid traits themselves. Finally, some plasmid traits and environments (e.g., resistance genes in human-related environments) were more often associated with mobilizable plasmids (those having at least one detected relaxase) than others. Overall, these findings underscore the high level of diversity of traits encoded by plasmids and provide a baseline to investigate the potential of plasmids to serve as reservoirs of adaptive traits for microbial communities. IMPORTANCE Plasmids are well known for their role in the transmission of antibiotic resistance-conferring genes. Beyond human and clinical settings, however, they disseminate many other types of genes, including those that contribute to microbially driven ecosystem processes. In this study, we identified the distribution of traits genetically encoded by plasmids isolated from seven broadly categorized environments. We find that plasmid trait content varied with both bacterial host taxonomy and environment and that, on average, half of the plasmids were potentially mobilizable. As anthropogenic activities impact ecosystems and the climate, investigating and identifying the mechanisms of how microbial communities can adapt will be imperative for predicting the impacts on ecosystem functioning.

About the dark corners in the gene function space of Escherichia coli remaining without illumination by scientific literature

BACKGROUND

Although Escherichia coli (E. coli) is the most studied prokaryote organism in the history of life sciences, many molecular mechanisms and gene functions encoded in its genome remain to be discovered. This work aims at quantifying the illumination of the E. coli gene function space by the scientific literature and how close we are towards the goal of a complete list of E. coli gene functions.

RESULTS

The scientific literature about E. coli protein-coding genes has been mapped onto the genome via the mentioning of names for genomic regions in scientific articles both for the case of the strain K-12 MG1655 as well as for the 95%-threshold softcore genome of 1324 E. coli strains with known complete genome. The article match was quantified with the ratio of a given gene name's occurrence to the mentioning of any gene names in the paper. The various genome regions have an extremely uneven literature coverage. A group of elite genes with ≥ 100 full publication equivalents (FPEs, FPE = 1 is an idealized publication devoted to just a single gene) attracts the lion share of the papers. For K-12, ~ 65% of the literature covers just 342 elite genes; for the softcore genome, ~ 68% of the FPEs is about only 342 elite gene families (GFs). We also find that most genes/GFs have at least one mentioning in a dedicated scientific article (with the exception of at least 137 protein-coding transcripts for K-12 and 26 GFs from the softcore genome). Whereas the literature growth rates were highest for uncharacterized or understudied genes until 2005-2010 compared with other groups of genes, they became negative thereafter. At the same time, literature for anyhow well-studied genes started to grow explosively with threshold T10 (≥ 10 FPEs). Typically, a body of ~ 20 actual articles generated over ~ 15 years of research effort was necessary to reach T10. Lineage-specific co-occurrence analysis of genes belonging to the accessory genome of E. coli together with genomic co-localization and sequence-analytic exploration hints previously completely uncharacterized genes yahV and yddL being associated with osmotic stress response/motility mechanisms.

CONCLUSION

If the numbers of scientific articles about uncharacterized and understudied genes remain at least at present levels, full gene function lists for the strain K-12 MG1655 and the E. coli softcore genome are in reach within the next 25-30 years. Once the literature body for a gene crosses 10 FPEs, most of the critical fundamental research risk appears overcome and steady incremental research becomes possible.

Identification of a novel 5-aminomethyl-2-thiouridine methyltransferase in tRNA modification

The uridine at the 34th position of tRNA, which is able to base pair with the 3'-end codon on mRNA, is usually modified to influence many aspects of decoding properties during translation. Derivatives of 5-methyluridine (xm5U), which include methylaminomethyl (mnm-) or carboxymethylaminomethyl (cmnm-) groups at C5 of uracil base, are widely conserved at the 34th position of many prokaryotic tRNAs. In Gram-negative bacteria such as Escherichia coli, a bifunctional MnmC is involved in the last two reactions of the biosynthesis of mnm5(s2)U, in which the enzyme first converts cmnm5(s2)U to 5-aminomethyl-(2-thio)uridine (nm5(s2)U) and subsequently installs the methyl group to complete the formation of mnm5(s2)U. Although mnm5s2U has been identified in tRNAs of Gram-positive bacteria and plants as well, their genomes do not contain an mnmC ortholog and the gene(s) responsible for this modification is unknown. We discovered that MnmM, previously known as YtqB, is the methyltransferase that converts nm5s2U to mnm5s2U in Bacillus subtilis through comparative genomics, gene complementation experiments, and in vitro assays. Furthermore, we determined X-ray crystal structures of MnmM complexed with anticodon stem loop of tRNAGln. The structures provide the molecular basis underlying the importance of U33-nm5s2U34-U35 as the key determinant for the specificity of MnmM.

Differences in the regulatory strategies of marine oligotrophs and copiotrophs reflect differences in motility

Aquatic bacteria frequently are divided into lifestyle categories oligotroph or copiotroph. Oligotrophs have proportionately fewer transcriptional regulatory genes than copiotrophs and are generally non-motile/chemotactic. We hypothesized that the absence of chemotaxis/motility in oligotrophs prevents them from occupying nutrient patches long enough to benefit from transcriptional regulation. We first confirmed that marine oligotrophs are generally reduced in genes for transcriptional regulation and motility/chemotaxis. Next, using a non-motile oligotroph (Ca. Pelagibacter st. HTCC7211), a motile copiotroph (Alteromonas macleodii st. HOT1A3), and [¹⁴ C]l-alanine, we confirmed that l-alanine catabolism is not transcriptionally regulated in HTCC7211 but is in HOT1A3. We then found that HOT1A3 took 2.5-4 min to initiate l-alanine oxidation at patch l-alanine concentrations, compared to <30 s for HTCC7211. By modelling cell trajectories, we predicted that, in most scenarios, non-motile cells spend <2 min in patches, compared to >4 min for chemotactic/motile cells. Thus, the time necessary for transcriptional regulation to initiate prevents transcriptional regulation from being beneficial for non-motile oligotrophs. This is supported by a mechanistic model we developed, which predicted that HTCC7211 cells with transcriptional regulation of l-alanine metabolism would produce 12% of their standing ATP stock upon encountering an l-alanine patch, compared to 880% in HTCC7211 cells without transcriptional regulation.

Importance of mobile genetic element immunity in numerically abundant Trichodesmium clades

The colony-forming cyanobacteria Trichodesmium spp. are considered one of the most important nitrogen-fixing genera in the warm, low nutrient ocean. Despite this central biogeochemical role, many questions about their evolution, physiology, and trophic interactions remain unanswered. To address these questions, we describe Trichodesmium pangenomic potential via significantly improved genomic assemblies from two isolates and 15 new >50% complete Trichodesmium metagenome-assembled genomes from hand-picked, Trichodesmium colonies spanning the Atlantic Ocean. Phylogenomics identified ~four N2 fixing clades of Trichodesmium across the transect, with T. thiebautii dominating the colony-specific reads. Pangenomic analyses showed that all T. thiebautii MAGs are enriched in COG defense mechanisms and encode a vertically inherited Type III-B Clustered Regularly Interspaced Short Palindromic Repeats and associated protein-based immunity system (CRISPR-Cas). Surprisingly, this CRISPR-Cas system was absent in all T. erythraeum genomes, vertically inherited by T. thiebautii, and correlated with increased signatures of horizontal gene transfer. Additionally, the system was expressed in metaproteomic and transcriptomic datasets and CRISPR spacer sequences with 100% identical hits to field-assembled, putative phage genome fragments were identified. While the currently CO2-limited T. erythraeum is expected to be a 'winner' of anthropogenic climate change, their genomic dearth of known phage resistance mechanisms, compared to T. thiebautii, could put this outcome in question. Thus, the clear demarcation of T. thiebautii maintaining CRISPR-Cas systems, while T. erythraeum does not, identifies Trichodesmium as an ecologically important CRISPR-Cas model system, and highlights the need for more research on phage-Trichodesmium interactions.

Site Specialization of Human Oral Veillonella Species

Veillonella species are abundant members of the human oral microbiome with multiple interspecies commensal relationships. Examining the distribution patterns of Veillonella species across the oral cavity is fundamental to understanding their oral ecology. In this study, we used a combination of pangenomic analysis and oral metagenomic information to clarify Veillonella taxonomy and to test the site specialist hypothesis for the Veillonella genus, which contends that most oral bacterial species are adapted to live at specific oral sites. Using isolate genome sequences combined with shotgun metagenomic sequence data, we showed that Veillonella species have clear, differential site specificity: Veillonella parvula showed strong preference for supra- and subgingival plaque, while closely related V. dispar, as well as more distantly related V. atypica, preferred the tongue dorsum, tonsils, throat, and hard palate. In addition, the provisionally named Veillonella sp. Human Microbial Taxon 780 showed strong site specificity for keratinized gingiva. Using comparative genomic analysis, we identified genes associated with thiamine biosynthesis and the reductive pentose phosphate cycle that may enable Veillonella species to occupy their respective habitats. IMPORTANCE Understanding the microbial ecology of the mouth is fundamental for understanding human physiology. In this study, metapangenomics demonstrated that different Veillonella species have clear ecological preferences in the oral cavity of healthy humans, validating the site specialist hypothesis. Furthermore, the gene pool of different Veillonella species was found to be reflective of their ecology, illuminating the potential role of vitamins and carbohydrates in determining Veillonella distribution patterns and interspecies interactions.

Isolation and Genomics of Futiania mangrovii gen. nov., sp. nov., a Rare and Metabolically Versatile Member in the Class Alphaproteobacteria

Mangrove microorganisms are a major part of the coastal ecosystem and are directly associated with nutrient cycling. Despite their ecological significance, the collection of culturable mangrove microbes is limited due to difficulties in isolation and cultivation. Here, we report the isolation and genome sequence of strain FT118^T, the first cultured representative of a previously uncultivated order UBA8317 within Alphaproteobacteria, based on the combined results of 16S rRNA gene similarity, phylogenomic, and average amino acid identity analyses. We propose Futianiales ord. nov. and Futianiaceae fam. nov. with Futiania as the type genus, and FT118^T represents the type species with the name Futiania mangrovii gen. nov, sp. nov. The 16S rRNA gene sequence comparison reveals that this novel order is a rare member but has a ubiquitous distribution across various habitats worldwide, which is corroborated by the experimental confirmation that this isolate can physiologically adapt to a wide range of oxygen levels, temperatures, pH and salinity levels. Biochemical characterization, genomic annotation, and metatranscriptomic analysis of FT118^T demonstrate that it is metabolically versatile and active in situ. Genomic analysis reveals adaptive features of Futianiales to fluctuating mangrove environments, including the presence of high- and low-affinity terminal oxidases, N-type ATPase, and the genomic capability of producing various compatible solutes and polyhydroxybutyrate, which possibly allow for the persistence of this novel order across various habitats. Collectively, these results expand the current culture collection of mangrove microorganisms, providing genomic insights of how this novel taxon adapts to fluctuating environments and the culture reference to unravel possible microbe-environment interactions. IMPORTANCE The rare biosphere constitutes an essential part of the microbial community and may drive nutrient cycling and other geochemical processes. However, the difficulty in microbial isolation and cultivation has hampered our understanding of the physiology and ecology of uncultured rare lineages. In this study, we successfully isolated a novel alphaproteobacterium, designated as FT118^T, and performed a combination of phenotypic, phylogenetic, and phylogenomic analyses, confirming that this isolate represents the first cultured member of a previously uncultivated order UBA8317 within Alphaproteobacteria. It is a rare species with a ubiquitous distribution across different habitats. Genomic and metatranscriptomic analyses demonstrate that it is metabolically versatile and active in situ, suggesting its potential role in nutrient cycling despite being scarce. This work not only expands the current phylogeny of isolated Alphaproteobacteria but also provides genomic and culture reference to unravel microbial adaptation strategies in mangrove sediments and possible microbe-environment interactions.

Genome-scale model of Pseudomonas aeruginosa metabolism unveils virulence and drug potentiation

Pseudomonas aeruginosa is one of the leading causes of hospital-acquired infections. To decipher the metabolic mechanisms associated with virulence and antibiotic resistance, we have developed an updated genome-scale model (GEM) of P. aeruginosa. The model (iSD1509) is an extensively curated, three-compartment, and mass-and-charge balanced BiGG model containing 1509 genes, the largest gene content for any P. aeruginosa GEM to date. It is the most accurate with prediction accuracies as high as 92.4% (gene essentiality) and 93.5% (substrate utilization). In iSD1509, we newly added a recently discovered pathway for ubiquinone-9 biosynthesis which is required for anaerobic growth. We used a modified iSD1509 to demonstrate the role of virulence factor (phenazines) in the pathogen survival within biofilm/oxygen-limited condition. Further, the model can mechanistically explain the overproduction of a drug susceptibility biomarker in the P. aeruginosa mutants. Finally, we use iSD1509 to demonstrate the drug potentiation by metabolite supplementation, and elucidate the mechanisms behind the phenotype, which agree with experimental results.

Illumina RNA and SMRT Sequencing Reveals the Mechanism of Uptake and Transformation of Selenium Nanoparticles in Soybean Seedlings

Selenium (Se) is an essential element for mammals, and its deficiency in the diet is a global problem. Agronomic biofortification through exogenous Se provides a valuable strategy to enhance human Se intake. Selenium nanoparticles (SeNPs) have been regarded to be higher bioavailability and less toxicity in comparison with selenite and selenate. Still, little has been known about the mechanism of their metabolism in plants. Soybean (Glycine max L.) can enrich Se, providing an ideal carrier for Se biofortification. In this study, soybean sprouts were treated with SeNPs, and a combination of next-generation sequencing (NGS) and single-molecule real-time (SMRT) sequencing was applied to clarify the underlying molecular mechanism of SeNPs metabolism. A total of 74,662 nonredundant transcripts were obtained, and 2109 transcription factors, 9687 alternative splice events, and 3309 long non-coding RNAs (lncRNAs) were predicted, respectively. KEGG enrichment analysis of the DEGs revealed that metabolic pathways, biosynthesis of secondary metabolites, and peroxisome were most enriched both in roots and leaves after exposure to SeNPs. A total of 117 transcripts were identified to be putatively involved in SeNPs transport and biotransformation in soybean. The top six hub genes and their closely coexpressed Se metabolism-related genes, such as adenylylsulfate reductase (APR3), methionine-tRNA ligase (SYM), and chloroplastic Nifs-like cysteine desulfurases (CNIF1), were screened by WGCNA and identified to play crucial roles in SeNPs accumulation and tolerance in soybean. Finally, a putative metabolism pathway of SeNPs in soybean was proposed. These findings have provided a theoretical foundation for future elucidation of the mechanism of SeNPs metabolism in plants.

Gallionellaceae pangenomic analysis reveals insight into phylogeny, metabolic flexibility, and iron oxidation mechanisms

The iron-oxidizing Gallionellaceae drive a wide variety of biogeochemical cycles through their metabolisms and biominerals. To better understand the environmental impacts of Gallionellaceae, we need to improve our knowledge of their diversity and metabolisms, especially any novel iron oxidation mechanisms. Here, we used a pangenomic analysis of 103 genomes to resolve Gallionellaceae phylogeny and explore the range of genomic potential. Using a concatenated ribosomal protein tree and key gene patterns, we determined Gallionellaceae has four genera, divided into two groups-iron-oxidizing bacteria (FeOB) Gallionella, Sideroxydans, and Ferriphaselus with known iron oxidases (Cyc2, MtoA) and nitrite-oxidizing bacteria (NOB) Candidatus Nitrotoga with nitrite oxidase (Nxr). The FeOB and NOB have similar electron transport chains, including genes for reverse electron transport and carbon fixation. Auxiliary energy metabolisms including S oxidation, denitrification, and organotrophy were scattered throughout the Gallionellaceae FeOB. Within FeOB, we found genes that may represent adaptations for iron oxidation, including a variety of extracellular electron uptake (EEU) mechanisms. FeOB genomes encoded more predicted c-type cytochromes overall, notably more multiheme c-type cytochromes (MHCs) with >10 CXXCH motifs. These include homologs of several predicted outer membrane porin-MHC complexes, including MtoAB and Uet. MHCs are known to efficiently conduct electrons across longer distances and function across a wide range of redox potentials that overlap with mineral redox potentials, which can help expand the range of usable iron substrates. Overall, the results of pangenome analyses suggest that the Gallionellaceae genera Gallionella, Sideroxydans, and Ferriphaselus are primarily iron oxidizers, capable of oxidizing dissolved Fe2+ as well as a range of solid iron or other mineral substrates.

Complet+: a computationally scalable method to improve completeness of large-scale protein sequence clustering

A major challenge for clustering algorithms is to balance the trade-off between homogeneity, i.e., the degree to which an individual cluster includes only related sequences, and completeness, the degree to which related sequences are broken up into multiple clusters. Most algorithms are conservative in grouping sequences with other sequences. Remote homologs may fail to be clustered together and instead form unnecessarily distinct clusters. The resulting clusters have high homogeneity but completeness that is too low. We propose Complet+, a computationally scalable post-processing method to increase the completeness of clusters without an undue cost in homogeneity. Complet+ proves to effectively merge closely-related clusters of protein that have verified structural relationships in the SCOPe classification scheme, improving the completeness of clustering results at little cost to homogeneity. Applying Complet+ to clusters obtained using MMseqs2's clusterupdate achieves an increased V-measure of 0.09 and 0.05 at the SCOPe superfamily and family levels, respectively. Complet+ also creates more biologically representative clusters, as shown by a substantial increase in Adjusted Mutual Information (AMI) and Adjusted Rand Index (ARI) metrics when comparing predicted clusters to biological classifications. Complet+ similarly improves clustering metrics when applied to other methods, such as CD-HIT and linclust. Finally, we show that Complet+ runtime scales linearly with respect to the number of clusters being post-processed on a COG dataset of over 3 million sequences. Code and supplementary information is available on Github: https://github.com/EESI/Complet-Plus.

Insights into the genome of Methylobacterium sp. NMS14P, a novel bacterium for growth promotion of maize, chili, and sugarcane

A novel methylotrophic bacterium designated as NMS14P was isolated from the root of an organic coffee plant (Coffea arabica) in Thailand. The 16S rRNA sequence analysis revealed that this new isolate belongs to the genus Methylobacterium, and its novelty was clarified by genomic and comparative genomic analyses, in which NMS14P exhibited low levels of relatedness with other Methylobacterium-type strains. NMS14P genome consists of a 6,268,579 bp chromosome, accompanied by a 542,519 bp megaplasmid and a 66,590 bp plasmid, namely pNMS14P1 and pNMS14P2, respectively. Several genes conferring plant growth promotion are aggregated on both chromosome and plasmids, including phosphate solubilization, indole-3-acetic acid (IAA) biosynthesis, cytokinins (CKs) production, 1-aminocyclopropane-1-carboxylate (ACC) deaminase activity, sulfur-oxidizing activity, trehalose synthesis, and urea metabolism. Furthermore, pangenome analysis showed that NMS14P possessed the highest number of strain-specific genes accounting for 1408 genes, particularly those that are essential for colonization and survival in a wide array of host environments, such as ABC transporter, chemotaxis, quorum sensing, biofilm formation, and biosynthesis of secondary metabolites. In vivo tests have supported that NMS14P significantly promoted the growth and development of maize, chili, and sugarcane. Collectively, NMS14P is proposed as a novel plant growth-promoting Methylobacterium that could potentially be applied to a broad range of host plants as Methylobacterium-based biofertilizers to reduce and ultimately substitute the use of synthetic agrochemicals for sustainable agriculture.

Comparative Transcriptomics and Genomics from Continuous Axenic Media Growth Identifies Coxiella burnetii Intracellular Survival Strategies

Coxiella burnetii (Cb) is an obligate intracellular pathogen in nature and the causative agent of acute Q fever as well as chronic diseases. In an effort to identify genes and proteins crucial to their normal intracellular growth lifestyle, we applied a "Reverse evolution" approach where the avirulent Nine Mile Phase II strain of Cb was grown for 67 passages in chemically defined ACCM-D media and gene expression patterns and genome integrity from various passages was compared to passage number one following intracellular growth. Transcriptomic analysis identified a marked downregulation of the structural components of the type 4B secretion system (T4BSS), the general secretory (sec) pathway, as well as 14 out of 118 previously identified genes encoding effector proteins. Additional downregulated pathogenicity determinants genes included several chaperones, LPS, and peptidoglycan biosynthesis. A general marked downregulation of central metabolic pathways was also observed, which was balanced by a marked upregulation of genes encoding transporters. This pattern reflected the richness of the media and diminishing anabolic and ATP-generation needs. Finally, genomic sequencing and comparative genomic analysis demonstrated an extremely low level of mutation across passages, despite the observed Cb gene expression changes following acclimation to axenic media.

A single, peri-operative antibiotic can persistently alter the post-operative gut microbiome after Roux-en-Y gastric bypass

INTRODUCTION

Roux-en-Y gastric bypass (RYGB) significantly alters the gut microbiome and may be a mechanism for post-operative cardiovascular disease improvement. We have previously found an association between the class of peri-operative, intravenous antibiotic administered at the time of RYGB and the resolution rate of hypertension suggesting the gut microbiome as a mechanism. In this study, we performed a prospective study of RYGB to determine if a single intravenous antibiotic could alter the gastrointestinal microbial composition.

METHODS

Patients undergoing RYGB were randomized to a single, peri-operative antibiotic of intravenous cefazolin (n = 8) or clindamycin (n = 8). Stool samples were collected from four-time points: 2 weeks pre-op (- 2w), 2 days pre-op (- 2d), 2 weeks post-op (+ 2w) and 3 months post-op (+ 3m). Stool samples were processed for genomic DNA followed by Illumina 16S rRNA gene sequencing and shotgun metagenomic sequencing (MGS).

RESULTS

A total of 60 stool samples (- 2w, n = 16; - 2d, n = 15; + 2w, n = 16; + 3m, n = 13) from 16 patients were analyzed. 87.5% of patients were female with an average age of 48.6 ± 12.2 years and pre-operative BMI of 50.9 ± 23.3 kg/m². RYGB induced statistically significant differences in alpha and beta diversity. There were statistically significant differences in alpha diversity at + 2w and beta diversity at + 3m due to antibiotic treatment. MGS revealed significantly distinct gut microbiota with 11 discriminatory metagenomic assembled genomes driven by antibiotic treatment at 3 months post-op, including increased Bifidobacterium spp. with clindamycin.

CONCLUSION

RYGB induces significant changes in the gut microbiome at 2 weeks that are maintained 3 months after surgery. However, the single peri-operative dose of antibiotic administered at the time of RYGB induces unique and persisting changes to the gut microbiome that are antibiotic-specific. Increased Bifidobacterium spp. with clindamycin administration may improve the metabolic efficacy of RYGB when considering gut-microbiome driven mechanisms for blood pressure resolution.

Insight into sulfamethoxazole effects on aerobic denitrification by strain Pseudomonas aeruginosa PCN-2: From simultaneous degradation performance to transcriptome analysis

It is a well-established fact that aerobic denitrifying strains are profoundly affected by antibiotics, but bacterium performing simultaneous aerobic denitrification and antibiotic degradation is hardly reported. Here, a typical aerobic denitrifying bacterium Pseudomonas aeruginosa PCN-2 was discovered to be capable of sulfamethoxazole (SMX) degradation. The results showed that nitrate removal efficiency was decreased from 100% to 88.12%, but the resistance of strain PCN-2 to SMX stress was enhanced with the increment of SMX concentration from 0 to 100 mg/L. Transcriptome analysis revealed that the down-regulation of energy metabolism pathways rather than the denitrifying functional genes was responsible for the suppressed nitrogen removal, while the up-regulation of antibiotic resistance pathways (e.g., biofilm formation, multi-drug efflux system, and quorum sensing) ensured the survival of bacterium and the carrying out of aerobic denitrification. Intriguingly, strain PCN-2 could degrade SMX during aerobic denitrification. Seven metabolites were identified by the UHPLC-MS, and three degradation pathways (which includes a new pathway that has never been reported) was proposed combined with the expressions of drug metabolic genes (e.g., cytP450, FMN, ALDH and NAT). This work provides a mechanistic understanding of the metabolic adaption of strain PCN-2 under SMX stress, which provided a broader idea for the treatment of SMX-containing wastewater.

Pyridoxal 5'-phosphate synthesis and salvage in Bacteria and Archaea: predicting pathway variant distributions and holes

Pyridoxal 5’-phosphate or PLP is a cofactor derived from B6 vitamers and essential for growth in all known organisms. PLP synthesis and salvage pathways are well characterized in a few model species even though key components, such as the vitamin B6 transporters, are still to be identified in many organisms including the model bacteria Escherichia coli or Bacillus subtilis . Using a comparative genomic approach, PLP synthesis and salvage pathways were predicted in 5840 bacterial and archaeal species with complete genomes. The distribution of the two known de novo biosynthesis pathways and previously identified cases of non-orthologous displacements were surveyed in the process. This analysis revealed that several PLP de novo pathway genes remain to be identified in many organisms, either because sequence similarity alone cannot be used to discriminate among several homologous candidates or due to non-orthologous displacements. Candidates for some of these pathway holes were identified using published TnSeq data, but many remain. We find that ~10 % of the analysed organisms rely on salvage but further analyses will be required to identify potential transporters. This work is a starting point to model the exchanges of B6 vitamers in communities, predict the sensitivity of a given organism to drugs targeting PLP synthesis enzymes, and identify numerous gaps in knowledge that will need to be tackled in the years to come.

Systematic In Silico Assessment of Antimicrobial Resistance Dissemination across the Global Plasmidome

The emergence of pathogenic strains resistant to multiple antimicrobials is a pressing problem in modern healthcare. Antimicrobial resistance is mediated primarily by dissemination of resistance determinants via horizontal gene transfer. The dissemination of some resistance genes has been well documented, but few studies have analyzed the patterns underpinning the dissemination of antimicrobial resistance genes. Analyzing the %GC content of plasmid-borne antimicrobial resistance genes relative to their host genome %GC content provides a means to efficiently detect and quantify dissemination of antimicrobial resistance genes. In this work we automate %GC content analysis to perform a comprehensive analysis of known antimicrobial resistance genes in publicly available plasmid sequences. We find that the degree to which antimicrobial resistance genes are disseminated depends primarily on the resistance mechanism. Our analysis identifies conjugative plasmids as primary dissemination vectors and indicates that most broadly disseminated genes have spread from single genomic backgrounds. We show that resistance dissemination profiles vary greatly among antimicrobials, oftentimes reflecting stewardship measures. Our findings establish %GC content analysis as a powerful, intuitive and scalable method to monitor the dissemination of resistance determinants using publicly available sequence data.

Comparative transcriptomics and genomics from continuous axenic media growth identifies Coxiella burnetii intracellular survival strategies

Coxiella burnetii (Cb) is an obligate intracellular pathogen in nature and the causative agent of acute Q fever as well as chronic diseases. In an effort to identify genes and proteins crucial to their normal intracellular growth lifestyle, we applied a 'reverse evolution' approach where the avirulent Nine Mile Phase II strain of Cb was grown for 67 passages in chemically defined ACCM-D media and gene expression patterns and genome integrity from various passages was compared to passage number one following intracellular growth. Transcriptomic analysis identified a marked downregulation of the structural components of the type 4B secretion system (T4BSS), the general secretory (Sec) pathway, as well as 14 out of 118 previously identified genes encoding effector proteins. Additional downregulated pathogenicity determinants genes included several chaperones, LPS, and peptidoglycan biosynthesis. A general marked downregulation of central metabolic pathways was also observed, which was balanced by a marked upregulation of genes encoding transporters. This pattern reflected the richness of the media and diminishing anabolic, and ATP-generation needs. Finally, genomic sequencing and comparative genomic analysis demonstrated an extremely low level of mutation across passages, despite the observed Cb gene expression changes following acclimation to axenic media.

Tractidigestivibacter montrealensis sp. nov., a new member of human gut microbiota isolated from a healthy volunteer

Strain KD21T, isolated from the fecal sample of a healthy female volunteer, is a strictly anaerobic, non-motile, Gram-staining-positive, saccharolytic small rod that does not produce spores. Strain KD21T was able to grow in the range of temperature 28°C-37°C (optimum, 37 °C), pH 6.0-8.0 (optimum, pH 7.0), and with 0-5.0 g/l NaCl (optimum, 0 g/l NaCl). Bacteria cells reduced nitrates to nitrites. Its major fatty acids were C18:1ω9c, C16:0, C18:0, and summed in feature 8 (C18:1ω7c and/or C18:1ω6c). 16S rRNA gene phylogenetic analysis revealed that KD21T is a member of the genus Tractidigestivibacter and is distinct from any species with validly published names. The sequence showed 98.48% similarity with T. scatoligenes SK9K4T. The DNA G + C content of strain KD21T was 62.6 mol%. The DNA-DNA hybridization and OrthoANI values between strain KD21T and T. scatoligenes SK9K4T were 40.2% and 90.2%, respectively. Differences in phenotypic, phylogenetic, chemotaxonomic, and genomic characteristics indicated that strain KD21T represents a novel species within the genus Tractidigestivibacter. The name T. montrealensis sp. nov. is proposed and the type strain is KD21T (= CSUR Q8103T =  DSM 115111T).

Aging-associated augmentation of gut microbiome virulence capability drives sepsis severity

Prior research has focused on host factors as mediators of exaggerated sepsis-associated morbidity and mortality in older adults. This focus on the host, however, has failed to identify therapies that improve sepsis outcomes in the elderly. We hypothesized that the increased susceptibility of the aging population to sepsis is not only a function of the host, but also reflects longevity-associated changes in the virulence of gut pathobionts. We utilized two complementary models of gut microbiota-induced experimental sepsis to establish the aged gut microbiome as a key pathophysiologic driver of heightened disease severity. Further murine and human investigations into these polymicrobial bacterial communities demonstrated that age was associated with only subtle shifts in ecological composition, but an overabundance of genomic virulence factors that have functional consequence on host immune evasion.

One Sentence Summary

The severity of sepsis in the aged host is in part mediated by longevity-associated increases in gut microbial virulence.

Into the darkness of the microbial dark matter in situ activities through expression profiles of Patescibacteria populations

Patescibacteria form a highly diverse and widespread superphylum of uncultured microorganisms representing a third of the global microbial diversity. Most of our knowledge on Patescibacteria putative physiology relies on metagenomic mining and metagenome-assembled genomes, but the in situ activities and the ecophysiology of these microorganisms have been rarely explored, leaving the role of Patescibacteria in ecosystems elusive. Using a genome-centric metatranscriptomic approach, we analyzed the diel and seasonal gene transcription profiles of 18 Patescibacteria populations in brackish microbial mats to test whether our understanding of Patescibacteria metabolism allows the extrapolation of their in situ activities. Although our results revealed a circadian cycle in Patescibacteria activities, a strong streamlined genetic expression characterized the Patescibacteria populations. This result has a major consequence for the extrapolation of their physiology and environmental function since most transcribed genes were uncharacterized, indicating that the ecophysiology of Patescibacteria cannot be yet reliably predicted from genomic data.

eggNOG 6.0: enabling comparative genomics across 12 535 organisms

The eggNOG (evolutionary gene genealogy Non-supervised Orthologous Groups) database is a bioinformatics resource providing orthology data and comprehensive functional information for organisms from all domains of life. Here, we present a major update of the database and website (version 6.0), which increases the number of covered organisms to 12 535 reference species, expands functional annotations, and implements new functionality. In total, eggNOG 6.0 provides a hierarchy of over 17M orthologous groups (OGs) computed at 1601 taxonomic levels, spanning 10 756 bacterial, 457 archaeal and 1322 eukaryotic organisms. OGs have been thoroughly annotated using recent knowledge from functional databases, including KEGG, Gene Ontology, UniProtKB, BiGG, CAZy, CARD, PFAM and SMART. eggNOG also offers phylogenetic trees for all OGs, maximising utility and versatility for end users while allowing researchers to investigate the evolutionary history of speciation and duplication events as well as the phylogenetic distribution of functional terms within each OG. Furthermore, the eggNOG 6.0 website contains new functionality to mine orthology and functional data with ease, including the possibility of generating phylogenetic profiles for multiple OGs across species or identifying single-copy OGs at custom taxonomic levels. eggNOG 6.0 is available at http://eggnog6.embl.de.

The ESKAPE mobilome contributes to the spread of antimicrobial resistance and CRISPR-mediated conflict between mobile genetic elements

Mobile genetic elements (MGEs) mediate the shuffling of genes among organisms. They contribute to the spread of virulence and antibiotic resistance (AMR) genes in human pathogens, such as the particularly problematic group of ESKAPE pathogens. Here, we performed the first systematic analysis of MGEs, including plasmids, prophages, and integrative and conjugative/mobilizable elements (ICEs/IMEs), across all ESKAPE pathogens. We found that different MGE types are asymmetrically distributed across these pathogens, and that most horizontal gene transfer (HGT) events are restricted by phylum or genus. We show that the MGEs proteome is involved in diverse functional processes and distinguish widespread proteins within the ESKAPE context. Moreover, anti-CRISPRs and AMR genes are overrepresented in the ESKAPE mobilome. Our results also underscore species-specific trends shaping the number of MGEs, AMR, and virulence genes across pairs of conspecific ESKAPE genomes with and without CRISPR-Cas systems. Finally, we observed that CRISPR spacers found on prophages, ICEs/IMEs, and plasmids have different targeting biases: while plasmid and prophage CRISPRs almost exclusively target other plasmids and prophages, respectively, ICEs/IMEs CRISPRs preferentially target prophages. Overall, our study highlights the general importance of the ESKAPE mobilome in contributing to the spread of AMR and mediating conflict among MGEs.

Intestinal epithelial cell-derived components regulate transcriptome of Lactobacillus rhamnosus GG

Introduction

Intestinal epithelial cells (IECs) provide the frontline responses to the gut microbiota for maintaining intestinal homeostasis. Our previous work revealed that IEC-derived components promote the beneficial effects of a commensal and probiotic bacterium, Lactobacillus rhamnosus GG (LGG). This study aimed to elucidate the regulatory effects of IEC-derived components on LGG at the molecular level.

Methods

Differential gene expression in LGG cultured with IEC-derived components at the timepoint between the exponential and stationary phase was studied by RNA sequencing and functional analysis.

Results

The transcriptomic profile of LGG cultured with IEC-derived components was significantly different from that of control LGG, with 231 genes were significantly upregulated and 235 genes significantly down regulated (FDR <0.05). The Clusters of Orthologous Groups (COGs) and Gene Ontology (GO) analysis demonstrated that the predominant genes enriched by IEC-derived components are involved in nutrient acquisition, including transporters for amino acids, metals, and sugars, biosynthesis of amino acids, and in the biosynthesis of cell membrane and cell wall, including biosynthesis of fatty acid and lipoteichoic acid. In addition, genes associated with cell division and translation are upregulated by IEC-derived components. The outcome of the increased transcription of these genes is supported by the result that IEC-derived components significantly promoted LGG growth. The main repressed genes are associated with the metabolism of amino acids, purines, carbohydrates, glycerophospholipid, and transcription, which may reflect regulation of metabolic mechanisms in response to the availability of nutrients in bacteria.

Discussion

These results provide mechanistic insight into the interactions between the gut microbiota and the host.

Human small-intestinal gluten-degrading bacteria and its potential implication in celiac disease

Celiac disease (CeD) is an immune-mediated chronic disorder triggered by the ingestion of wheat gluten in genetically predisposed individuals. Gluten is a major food ingredient, infamously containing proline and glutamine-rich domains that are highly resistant to digestion by mammalian proteolytic enzymes. Thus, adhering to a gluten-free diet (GFD) is the only known treatment for CeD, albeit with many complications. Therefore, any therapy that eliminates the gluten immunogenic part before it reaches the small intestine is highly desirable. Probiotic therapy containing gluten-degrading bacteria (GDB) and their protease enzymes are possibly new approaches to treating CeD. Our study aimed to identify novel GDB from the duodenal biopsy of the first-degree relative (FDR) subjects (relatives of diseased individuals who are healthy but susceptible to celiac disease) with the potential to reduce gluten immunogenicity. Using the gluten agar plate technique, bacterial strains Brevibacterium casei NAB46 and Staphylococcus arlettae R2AA77 displaying glutenase activity were screened, identified, and characterized. Whole-genome sequencing found gluten-degrading prolyl endopeptidase (PEP) in the B. casei NAB46 genome and glutamyl endopeptidase (GEP) in the S. arlettae R2AA77 genome. Partially purified PEP has a specific activity of 1.15 U/mg, while GEP has a specific activity of 0.84 U/mg, which are, respectively, 6- and 9-fold times higher after concentrating the enzymes. Our results showed that these enzymes could hydrolyse immunotoxic gliadin peptides recognized in western blot using an anti-gliadin antibody. Additionally, a docking model was proposed for representative gliadin peptide PQPQLPYPQPQLP in the active site of the enzymes, where the residues of the N-terminal peptide extensively interact with the catalytic domain of the enzymes. These bacteria and their associated glutenase enzymes efficiently neutralize gliadin immunogenic epitopes, opening possibilities for their application as a dietary supplement in treating CeD patients.

A census of actin-associated proteins in humans

Actin filaments help in maintaining the cell structure and coordinating cellular movements and cargo transport within the cell. Actin participates in the interaction with several proteins and also with itself to form the helical filamentous actin (F-actin). Actin-binding proteins (ABPs) and actin-associated proteins (AAPs) coordinate the actin filament assembly and processing, regulate the flux between globular G-actin and F-actin in the cell, and help maintain the cellular structure and integrity. We have used protein-protein interaction data available through multiple sources (STRING, BioGRID, mentha, and a few others), functional annotation, and classical actin-binding domains to identify actin-binding and actin-associated proteins in the human proteome. Here, we report 2482 AAPs and present an analysis of their structural and sequential domains, functions, evolutionary conservation, cellular localization, abundance, and tissue-specific expression patterns. This analysis provides a base for the characterization of proteins involved in actin dynamics and turnover in the cell.

Diversity of the genus Cryobacterium and proposal of 19 novel species isolated from glaciers

The bacterial genus Cryobacterium includes at present 14 species that live in cryospheric environments. In this study, we analyzed 101 genomes of Cryobacterium with pure cultures obtained from GenBank. They could be classified into 44 species based on average nucleotide identity (ANI) analysis, showing the diversity of Cryobacterium. Among these, 19 strains in our laboratory were isolated from the glacier samples in China. The pairwise ANI values of these 19 strains and known species were <95%, indicating that they represented 19 novel species. The comparative genomic analysis showed significant differences in gene content between the two groups with a maximum growth temperature (T _max) of ≤ 20°C and a T _max of >20°C. A comprehensive and robust phylogenetic tree, including 14 known species and 19 novel species, was constructed and showed five phylogenetic branches based on 265 concatenated single-copy gene sequences. The T _max parameter had a strong phylogenetic signal, indicating that the temperature adaptation of Cryobacterium was largely through vertical transfer rather than horizontal gene transfer and was affected by selection. Furthermore, using polyphasic taxonomy combined with phylogenomic analysis, we proposed 19 novel species of the genus Cryobacterium by the following 19 names: Cryobacterium serini sp. nov., Cryobacterium lactosi sp. nov., Cryobacterium gelidum sp. nov., Cryobacterium suzukii sp. nov., Cryobacterium fucosi sp. nov., Cryobacterium frigoriphilum sp. nov., Cryobacterium cryoconiti sp. nov., Cryobacterium lyxosi sp. nov., Cryobacterium sinapicolor sp. nov., Cryobacterium sandaracinum sp. nov., Cryobacterium cheniae sp. nov., Cryobacterium shii sp. nov., Cryobacterium glucosi sp. nov., Cryobacterium algoritolerans sp. nov., Cryobacterium mannosilyticum sp. nov., Cryobacterium adonitolivorans sp. nov., Cryobacterium algoricola sp. nov., Cryobacterium tagatosivorans sp. nov., and Cryobacterium glaciale sp. nov. Overall, the taxonomy and genomic analysis can improve our knowledge of phenotypic diversity, genetic diversity, and evolutionary characteristics of Cryobacterium.

TransAAP: an automated annotation pipeline for membrane transporter prediction in bacterial genomes

Membrane transporters are a large group of proteins that span cell membranes and contribute to critical cell processes, including delivery of essential nutrients, ejection of waste products, and assisting the cell in sensing environmental conditions. Obtaining an accurate and specific annotation of the transporter proteins encoded by a micro-organism can provide details of its likely nutritional preferences and environmental niche(s), and identify novel transporters that could be utilized in small molecule production in industrial biotechnology. The Transporter Automated Annotation Pipeline (TransAAP) (http://www.membranetransport.org/transportDB2/TransAAP_login.html) is a fully automated web service for the prediction and annotation of membrane transport proteins in an organism from its genome sequence, by using comparisons with both curated databases such as the TCDB (Transporter Classification Database) and TDB, as well as selected Pfams and TIGRFAMs of transporter families and other methodologies. TransAAP was used to annotate transporter genes in the prokaryotic genomes in the National Center for Biotechnology Information (NCBI) RefSeq; these are presented in the transporter database TransportDB (http://www.membranetransport.org) website, which has a suite of data visualization and analysis tools. Creation and maintenance of a bioinformatic database specific for transporters in all genomic datasets is essential for microbiology research groups and the general research/biotechnology community to obtain a detailed picture of membrane transporter systems in various environments, as well as comprehensive information on specific membrane transport proteins.

Metaphor-A workflow for streamlined assembly and binning of metagenomes

Recent advances in bioinformatics and high-throughput sequencing have enabled the large-scale recovery of genomes from metagenomes. This has the potential to bring important insights as researchers can bypass cultivation and analyze genomes sourced directly from environmental samples. There are, however, technical challenges associated with this process, most notably the complexity of computational workflows required to process metagenomic data, which include dozens of bioinformatics software tools, each with their own set of customizable parameters that affect the final output of the workflow. At the core of these workflows are the processes of assembly-combining the short-input reads into longer, contiguous fragments (contigs)-and binning, clustering these contigs into individual genome bins. The limitations of assembly and binning algorithms also pose different challenges depending on the selected strategy to execute them. Both of these processes can be done for each sample separately or by pooling together multiple samples to leverage information from a combination of samples. Here we present Metaphor, a fully automated workflow for genome-resolved metagenomics (GRM). Metaphor differs from existing GRM workflows by offering flexible approaches for the assembly and binning of the input data and by combining multiple binning algorithms with a bin refinement step to achieve high-quality genome bins. Moreover, Metaphor generates reports to evaluate the performance of the workflow. We showcase the functionality of Metaphor on different synthetic datasets and the impact of available assembly and binning strategies on the final results.