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

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.

Functional characterization of diverse type I-F CRISPR-associated transposons

CRISPR-Cas systems generally provide adaptive immunity in prokaryotes through RNA-guided degradation of foreign genetic elements like bacteriophages and plasmids. Recently, however, transposon-encoded and nuclease-deficient CRISPR-Cas systems were characterized and shown to be co-opted by Tn7-like transposons for CRISPR RNA-guided DNA transposition. As a genome engineering tool, these CRISPR-Cas systems and their associated transposon proteins can be deployed for programmable, site-specific integration of sizable cargo DNA, circumventing the need for DNA cleavage and homology-directed repair involving endogenous repair machinery. Here, we selected a diverse set of type I-F3 CRISPR-associated transposon systems derived from Gammaproteobacteria, predicted all components essential for transposition activity, and deployed them for functionality testing within Escherichia coli. Our results demonstrate that these systems possess a significant range of integration efficiencies with regards to temperature, transposon size, and flexible PAM requirements. Additionally, our findings support the categorization of these systems into functional compatibility groups for efficient and orthogonal RNA-guided DNA integration. This work expands the CRISPR-based toolbox with new CRISPR RNA-guided DNA integrases that can be applied to complex and extensive genome engineering efforts.

The IMG/M data management and analysis system v.7: content updates and new features

The Integrated Microbial Genomes & Microbiomes system (IMG/M: https://img.jgi.doe.gov/m/) at the Department of Energy (DOE) Joint Genome Institute (JGI) continues to provide support for users to perform comparative analysis of isolate and single cell genomes, metagenomes, and metatranscriptomes. In addition to datasets produced by the JGI, IMG v.7 also includes datasets imported from public sources such as NCBI Genbank, SRA, and the DOE National Microbiome Data Collaborative (NMDC), or submitted by external users. In the past couple years, we have continued our effort to help the user community by improving the annotation pipeline, upgrading the contents with new reference database versions, and adding new analysis functionalities such as advanced scaffold search, Average Nucleotide Identity (ANI) for high-quality metagenome bins, new cassette search, improved gene neighborhood display, and improvements to metatranscriptome data display and analysis. We also extended the collaboration and integration efforts with other DOE-funded projects such as NMDC and DOE Biology Knowledgebase (KBase).

OrthoDB v11: annotation of orthologs in the widest sampling of organismal diversity

OrthoDB provides evolutionary and functional annotations of genes in a diverse sampling of eukaryotes, prokaryotes, and viruses. Genomics continues to accelerate our exploration of gene diversity and orthology is the most precise way of bridging gene functional knowledge with the rapidly expanding universe of genomic sequences. OrthoDB samples the most diverse organisms with the best quality genomics data to provide the leading coverage of species diversity. This update of the underlying data to over 18 000 prokaryotes and almost 2000 eukaryotes with over 100 million genes propels the coverage to another level. This achievement also demonstrates the scalability of the underlying OrthoLoger software for delineation of orthologs, freely available from https://orthologer.ezlab.org. In addition to the ab-initio computations of gene orthology used for the OrthoDB release, the OrthoLoger software allows mapping of novel gene sets to precomputed orthologs and thereby links to their annotations. The LEMMI-style benchmarking of OrthoLoger ensures its state-of-the-art performance and is available from https://lemortho.ezlab.org. The OrthoDB web interface has been further developed to include a pairwise orthology view from any gene to any other sampled species. OrthoDB-computed evolutionary annotations as well as extensively collated functional annotations can be accessed via REST API or SPARQL/RDF, downloaded or browsed online from https://www.orthodb.org.

HGTree v2.0: a comprehensive database update for horizontal gene transfer (HGT) events detected by the tree-reconciliation method

HGTree is a database that provides horizontal gene transfer (HGT) event information on 2472 prokaryote genomes using the tree-reconciliation method. HGTree was constructed in 2015, and a large number of prokaryotic genomes have been additionally published since then. To cope with the rapid rise of prokaryotic genome data, we present HGTree v2.0 (http://hgtree2.snu.ac.kr), a newly updated version of our HGT database with much more extensive data, including a total of 20 536 completely sequenced non-redundant prokaryotic genomes, and more reliable HGT information results curated with various steps. As a result, HGTree v2.0 has a set of expanded data results of 6 361 199 putative horizontally transferred genes integrated with additional functional information such as the KEGG pathway, virulence factors and antimicrobial resistance. Furthermore, various visualization tools in the HGTree v2.0 database website provide intuitive biological insights, allowing the users to investigate their genomes of interest.

Viral community-wide auxiliary metabolic genes differ by lifestyles, habitats, and hosts

BACKGROUND

Viral-encoded auxiliary metabolic genes (AMGs) are important toolkits for modulating their hosts' metabolisms and the microbial-driven biogeochemical cycles. Although the functions of AMGs have been extensively reported in numerous environments, we still know little about the drivers that shape the viral community-wide AMG compositions in natural ecosystems. Exploring the drivers of viral community-wide AMG compositions is critical for a deeper understanding of the complex interplays among viruses, hosts, and the environments.

RESULTS

Here, we investigated the impact of viral lifestyles (i.e., lytic and lysogenic), habitats (i.e., water, particle, and sediment), and prokaryotic hosts on viral AMG profiles by utilizing metagenomic and metatranscriptomic techniques. We found that viral lifestyles were the most important drivers, followed by habitats and host identities. Specifically, irrespective of what habitats viruses came from, lytic viruses exhibited greater AMG diversity and tended to encode AMGs for chaperone biosynthesis, signaling proteins, and lipid metabolism, which could boost progeny reproduction, whereas temperate viruses were apt to encode AMGs for host survivability. Moreover, the lytic and temperate viral communities tended to mediate the microbial-driven biogeochemical cycles, especially nitrogen metabolism, in different manners via AMGs. When focusing on each lifestyle, we further found clear dissimilarity in AMG compositions between water and sediment, as well the divergent AMGs encoded by viruses infecting different host orders.

CONCLUSIONS

Overall, our study provides a first systematic characterization of the drivers of viral community-wide AMG compositions and further expands our knowledge of the distinct interactions of lytic and temperate viruses with their prokaryotic hosts from an AMG perspective, which is critical for understanding virus-host-environment interactions in natural conditions. Video Abstract.

A novel and dual digestive symbiosis scales up the nutrition and immune system of the holobiont Rimicaris exoculata

BACKGROUND

In deep-sea hydrothermal vent areas, deprived of light, most animals rely on chemosynthetic symbionts for their nutrition. These symbionts may be located on their cuticle, inside modified organs, or in specialized cells. Nonetheless, many of these animals have an open and functional digestive tract. The vent shrimp Rimicaris exoculata is fueled mainly by its gill chamber symbionts, but also has a complete digestive system with symbionts. These are found in the shrimp foregut and midgut, but their roles remain unknown. We used genome-resolved metagenomics on separate foregut and midgut samples, taken from specimens living at three contrasted sites along the Mid-Atlantic Ridge (TAG, Rainbow, and Snake Pit) to reveal their genetic potential.

RESULTS

We reconstructed and studied 20 Metagenome-Assembled Genomes (MAGs), including novel lineages of Hepatoplasmataceae and Deferribacteres, abundant in the shrimp foregut and midgut, respectively. Although the former showed streamlined reduced genomes capable of using mostly broken-down complex molecules, Deferribacteres showed the ability to degrade complex polymers, synthesize vitamins, and encode numerous flagellar and chemotaxis genes for host-symbiont sensing. Both symbionts harbor a diverse set of immune system genes favoring holobiont defense. In addition, Deferribacteres were observed to particularly colonize the bacteria-free ectoperitrophic space, in direct contact with the host, elongating but not dividing despite possessing the complete genetic machinery necessary for this.

CONCLUSION

Overall, these data suggest that these digestive symbionts have key communication and defense roles, which contribute to the overall fitness of the Rimicaris holobiont. Video Abstract.

AsgeneDB: a curated orthology arsenic metabolism gene database and computational tool for metagenome annotation

Arsenic (As) is the most ubiquitous toxic metalloid in nature. Microbe-mediated As metabolism plays an important role in global As biogeochemical processes, greatly changing its toxicity and bioavailability. While metagenomic sequencing may advance our understanding of the As metabolism capacity of microbial communities in different environments, accurate metagenomic profiling of As metabolism remains challenging due to low coverage and inaccurate definitions of As metabolism gene families in public orthology databases. Here we developed a manually curated As metabolism gene database (AsgeneDB) comprising 400 242 representative sequences from 59 As metabolism gene families, which are affiliated with 1653 microbial genera from 46 phyla. AsgeneDB achieved 100% annotation sensitivity and 99.96% annotation accuracy for an artificial gene dataset. We then applied AsgeneDB for functional and taxonomic profiling of As metabolism in metagenomes from various habitats (freshwater, hot spring, marine sediment and soil). The results showed that AsgeneDB substantially improved the mapping ratio of short reads in metagenomes from various environments. Compared with other databases, AsgeneDB provides more accurate, more comprehensive and faster analysis of As metabolic genes. In addition, we developed an R package, Asgene, to facilitate the analysis of metagenome sequencing data. Therefore, AsgeneDB and the associated Asgene package will greatly promote the study of As metabolism in microbial communities in various environments.

HGD: an integrated homologous gene database across multiple species

Homology is fundamental to infer genes' evolutionary processes and relationships with shared ancestry. Existing homolog gene resources vary in terms of inferring methods, homologous relationship and identifiers, posing inevitable difficulties for choosing and mapping homology results from one to another. Here, we present HGD (Homologous Gene Database, https://ngdc.cncb.ac.cn/hgd), a comprehensive homologs resource integrating multi-species, multi-resources and multi-omics, as a complement to existing resources providing public and one-stop data service. Currently, HGD houses a total of 112 383 644 homologous pairs for 37 species, including 19 animals, 16 plants and 2 microorganisms. Meanwhile, HGD integrates various annotations from public resources, including 16 909 homologs with traits, 276 670 homologs with variants, 398 573 homologs with expression and 536 852 homologs with gene ontology (GO) annotations. HGD provides a wide range of omics gene function annotations to help users gain a deeper understanding of gene function.

Genome sequence of the entomopathogenic Serratia entomophila isolate 626 and characterisation of the species specific itaconate degradation pathway

BACKGROUND

Isolates of Serratia entomophila and S. proteamaculans (Yersiniaceae) cause disease specific to the endemic New Zealand pasture pest, Costelytra giveni (Coleoptera: Scarabaeidae). Previous genomic profiling has shown that S. entomophila isolates appear to have conserved genomes and, where present, conserved plasmids. In the absence of C. giveni larvae, S. entomophila prevalence reduces in the soil over time, suggesting that S. entomophila has formed a host-specific relationship with C. giveni. To help define potential genetic mechanisms driving retention of the chronic disease of S. entomophila, the genome of the isolate 626 was sequenced, enabling the identification of unique chromosomal properties, and defining the gain/loss of accessory virulence factors relevant to pathogenicity to C. giveni larvae.

RESULTS

We report the complete sequence of S. entomophila isolate 626, a causal agent of amber disease in C. giveni larvae. The genome of S. entomophila 626 is 5,046,461 bp, with 59.1% G + C content and encoding 4,695 predicted CDS. Comparative analysis with five previously sequenced Serratia species, S. proteamaculans 336X, S. marcescens Db11, S. nematodiphila DH-S01, S. grimesii BXF1, and S. ficaria NBRC 102596, revealed a core of 1,165 genes shared. Further comparisons between S. entomophila 626 and S. proteamaculans 336X revealed fewer predicted phage-like regions and genomic islands in 626, suggesting less horizontally acquired genetic material. Genomic analyses revealed the presence of a four-gene itaconate operon, sharing a similar gene order as the Yersinia pestis ripABC complex. Assessment of a constructed 626::RipC mutant revealed that the operon confer a possible metabolic advantage to S. entomophila in the initial stages of C. giveni infection.

CONCLUSIONS

Evidence is presented where, relative to S. proteamaculans 336X, S. entomophila 626 encodes fewer genomic islands and phages, alluding to limited horizontal gene transfer in S. entomophila. Bioassay assessments of a S. entomophila-mutant with a targeted mutation of the itaconate degradation region unique to this species, found the mutant to have a reduced capacity to replicate post challenge of the C. giveni larval host, implicating the itaconate operon in establishment within the host.

Distribution and Genomic Variation of Thermophilic Cyanobacteria in Diverse Microbial Mats at the Upper Temperature Limits of Photosynthesis

Thermophilic cyanobacteria have been extensively studied in Yellowstone National Park (YNP) hot springs, particularly during decades of work on the thick laminated mats of Octopus and Mushroom springs. However, focused studies of cyanobacteria outside these two hot springs have been lacking, especially regarding how physical and chemical parameters along with community morphology influence the genomic makeup of these organisms. Here, we used a metagenomic approach to examine cyanobacteria existing at the upper temperature limit of photosynthesis. We examined 15 alkaline hot spring samples across six geographic areas of YNP, all with various physical and chemical parameters and community morphology. We recovered 22 metagenome-assembled genomes (MAGs) belonging to thermophilic cyanobacteria, notably an uncultured Synechococcus-like taxon recovered from a setting at the upper temperature limit of photosynthesis, 73°C, in addition to thermophilic Gloeomargarita. Furthermore, we found that three distinct groups of Synechococcus-like MAGs recovered from different temperature ranges vary in their genomic makeup. MAGs from the uncultured very-high-temperature (up to 73°C) Synechococcus-like taxon lack key nitrogen metabolism genes and have genes implicated in cellular stress responses that diverge from other Synechococcus-like MAGs. Across all parameters measured, temperature was the primary determinant of taxonomic makeup of recovered cyanobacterial MAGs. However, total Fe, community morphology, and biogeography played an additional role in the distribution and abundance of upper-temperature-limit-adapted Synechococcus-like MAGs. These findings expand our understanding of cyanobacterial diversity in YNP and provide a basis for interrogation of understudied thermophilic cyanobacteria. IMPORTANCE Oxygenic photosynthesis arose early in microbial evolution-approximately 2.5 to 3.5 billion years ago-and entirely reshaped the biological makeup of Earth. However, despite the span of time in which photosynthesis has been refined, it is strictly limited to temperatures below 73°C, a barrier that many other biological processes have been able to overcome. Furthermore, photosynthesis at temperatures above 56°C is limited to circumneutral and alkaline pH. Hot springs in Yellowstone National Park (YNP), which have a large diversity in temperatures, pH, and geochemistry, provide a natural laboratory to study thermophilic microbial mats and the cyanobacteria within. While cyanobacteria in YNP microbial mats have been studied for decades, a vast majority of the work has focused on two springs within the same geyser basin, both containing similar community morphologies. Thus, the drivers of cyanobacterial adaptations to the upper limits of photosynthesis across a variety of environmental parameters have been understudied. Our findings provide new insights into the influence of these parameters on both taxonomic diversity and genomic content of cyanobacteria across a range of hot spring samples.

Insights into Evolutionary, Genomic, and Biogeographic Characterizations of Chryseobacterium nepalense Represented by a Polyvinyl Alcohol-Degrading Bacterium, AC3

Chryseobacterium spp. are Gram-negative rods found ubiquitously in the environment, with certain species being reported as having unusual degrading properties. Polyvinyl alcohol (PVA) is used widely in industry but causes serious global environmental pollution. Here, we report the complete genome sequence of a novel bacterium, AC3, that efficiently degrades PVA. As the representative genome of Chryseobacterium nepalense, key genomic characteristics (e.g., mobile genetic elements, horizontal genes, genome-scale metabolic network, secondary metabolite biosynthesis gene clusters, and carbohydrate-active enzymes) were comprehensively investigated to reveal the potential genetic features of this species. Core genome phylogenetic analysis in combination with average nucleotide identity, average amino acid identity, and in silico DNA-DNA hybridization values provided an accurate taxonomic position of C. nepalense in the genus Chryseobacterium. Comparative genomic analysis of AC3 with closely related species suggested evolutionary dynamics characterized by a species-specific genetic repertoire, dramatic rearrangements, and evolutionary constraints driven by selective pressure, which facilitated the speciation and adaptative evolution of C. nepalense. Biogeographic characterization indicated that this species is ubiquitously distributed not only in soil habitats but also in a variety of other source niches. Bioinformatic analysis revealed the potential genetic basis of PVA degradation in AC3, which included six putative genes associated with the synthesis of PVA dehydrogenase, cytochrome c, oxidized PVA hydrolase, and secondary alcohol dehydrogenase. Our study reports the first complete genome of C. nepalense with PVA-degrading properties, providing comprehensive insights into the genomic characteristics of this species and increasing our understanding of the microbial degradation of PVA. IMPORTANCE Although PVA is a biodegradable polymer, the widespread use of PVA in global industrialization has resulted in serious environmental problems. To date, knowledge of effective and applicable PVA-degrading bacteria is limited, and thus, the discovery of novel PVA biodegraders is pertinent. Here, we isolated a novel bacterial strain, AC3, which efficiently degraded PVA. The complete genome of AC3 was sequenced as the first genome sequence of the species C. nepalense. Comparative genomic analysis was performed to comprehensively investigate the phylogenetic relationships, genome-scale metabolic network, key genomic characteristics associated with genomic evolution, evolutionary dynamics between AC3 and its close relatives, and biogeographic characterization of C. nepalense, particularly regarding the potential genetic basis of PVA degradation. These findings could advance our understanding of the genomic characteristics of C. nepalense and PVA bioremediation.

PAT: a comprehensive database of prokaryotic antimicrobial toxins

Antimicrobial toxins help prokaryotes win competitive advantages in intraspecific or interspecific conflicts and are also a critical factor affecting the pathogenicity of many pathogens that threaten human health. Although many studies have revealed that antagonism based on antimicrobial toxins plays a central role in prokaryotic life, a database on antimicrobial toxins remains lacking. Here, we present the prokaryotic antimicrobial toxin database (PAT, http://bioinfo.qd.sdu.edu.cn/PAT/), a comprehensive data resource collection on experimentally validated antimicrobial toxins. PAT has organized information, derived from the reported literature, on antimicrobial toxins, as well as the corresponding immunity proteins, delivery mechanisms, toxin activities, structural characteristics, sequences, etc. Moreover, we also predict potential antimicrobial toxins in prokaryotic reference genomes and show the taxonomic information and environmental distribution of typical antimicrobial toxins. These details have been fully incorporated into the PAT database, where users can browse, search, download, analyse and view informative statistics and detailed information. PAT resources have already been used in our prediction and identification of prokaryotic antimicrobial toxins and may contribute to promoting the efficient investigation of antimicrobial toxin functions, the discovery of novel antimicrobial toxins, and an improved understanding of the biological roles and significance of these toxins.

Immunopeptidomics-based design of mRNA vaccine formulations against Listeria monocytogenes

Listeria monocytogenes is a foodborne intracellular bacterial pathogen leading to human listeriosis. Despite a high mortality rate and increasing antibiotic resistance no clinically approved vaccine against Listeria is available. Attenuated Listeria strains offer protection and are tested as antitumor vaccine vectors, but would benefit from a better knowledge on immunodominant vector antigens. To identify novel antigens, we screen for Listeria peptides presented on the surface of infected human cell lines by mass spectrometry-based immunopeptidomics. In between more than 15,000 human self-peptides, we detect 68 Listeria immunopeptides from 42 different bacterial proteins, including several known antigens. Peptides presented on different cell lines are often derived from the same bacterial surface proteins, classifying these antigens as potential vaccine candidates. Encoding these highly presented antigens in lipid nanoparticle mRNA vaccine formulations results in specific CD8⁺ T-cell responses and induces protection in vaccination challenge experiments in mice. Our results can serve as a starting point for the development of a clinical mRNA vaccine against Listeria and aid to improve attenuated Listeria vaccines and vectors, demonstrating the power of immunopeptidomics for next-generation bacterial vaccine development.

Expansion of the global RNA virome reveals diverse clades of bacteriophages

High-throughput RNA sequencing offers broad opportunities to explore the Earth RNA virome. Mining 5,150 diverse metatranscriptomes uncovered >2.5 million RNA virus contigs. Analysis of >330,000 RNA-dependent RNA polymerases (RdRPs) shows that this expansion corresponds to a 5-fold increase of the known RNA virus diversity. Gene content analysis revealed multiple protein domains previously not found in RNA viruses and implicated in virus-host interactions. Extended RdRP phylogeny supports the monophyly of the five established phyla and reveals two putative additional bacteriophage phyla and numerous putative additional classes and orders. The dramatically expanded phylum Lenarviricota, consisting of bacterial and related eukaryotic viruses, now accounts for a third of the RNA virome. Identification of CRISPR spacer matches and bacteriolytic proteins suggests that subsets of picobirnaviruses and partitiviruses, previously associated with eukaryotes, infect prokaryotic hosts.

Physiological and transcriptomic analysis uncovers salinity stress mechanisms in a facultative crassulacean acid metabolism plant Dendrobium officinale

Dendrobium officinale is a precious medicinal Chinese herb that employs facultative crassulacean acid metabolism (CAM) and has a high degree of abiotic stress tolerance, but the molecular mechanism underlying the response of this orchid to abiotic stresses is poorly understood. In this study, we analyzed the root microstructure of D. officinale plantlets and verified the presence of chloroplasts by transmission electron microscopy. To obtain a more comprehensive overview of the molecular mechanism underlying their tolerance to abiotic stress, we performed whole-transcriptome sequencing of the roots of 10-month-old plantlets exposed to salt (NaCl) treatment in a time-course experiment (0, 4 and 12 h). The total of 7376 differentially expressed genes that were identified were grouped into three clusters (P < 0.05). Metabolic pathway analysis revealed that the expression of genes related to hormone (such as auxins, cytokinins, abscisic acid, ethylene and jasmonic acid) biosynthesis and response, as well as the expression of genes related to photosynthesis, amino acid and flavonoid metabolism, and the SOS pathway, were either up- or down-regulated after salt treatment. Additionally, we identified an up-regulated WRKY transcription factor, DoWRKY69, whose ectopic expression in Arabidopsis promoted seed germination under salt tress. Collectively, our findings provide a greater understanding of the salt stress response mechanisms in the roots of a facultative CAM plant. A number of candidate genes that were discovered may help plants to cope with salt stress when introduced via genetic engineering.

Multi-Omics Approaches and Resources for Systems-Level Gene Function Prediction in the Plant Kingdom

In higher plants, the complexity of a system and the components within and among species are rapidly dissected by omics technologies. Multi-omics datasets are integrated to infer and enable a comprehensive understanding of the life processes of organisms of interest. Further, growing open-source datasets coupled with the emergence of high-performance computing and development of computational tools for biological sciences have assisted in silico functional prediction of unknown genes, proteins and metabolites, otherwise known as uncharacterized. The systems biology approach includes data collection and filtration, system modelling, experimentation and the establishment of new hypotheses for experimental validation. Informatics technologies add meaningful sense to the output generated by complex bioinformatics algorithms, which are now freely available in a user-friendly graphical user interface. These resources accentuate gene function prediction at a relatively minimal cost and effort. Herein, we present a comprehensive view of relevant approaches available for system-level gene function prediction in the plant kingdom. Together, the most recent applications and sought-after principles for gene mining are discussed to benefit the plant research community. A realistic tabulation of plant genomic resources is included for a less laborious and accurate candidate gene discovery in basic plant research and improvement strategies.

Phase-variable Type I methyltransferase M.NgoAV from Neisseria gonorrhoeae FA1090 regulates phasevarion expression and gonococcal phenotype

The restriction-modification (RM) systems are compared to a primitive, innate, prokaryotic immune system, controlling the invasion by foreign DNA, composed of methyltransferase (MTase) and restriction endonuclease. The biological significance of RM systems extends beyond their defensive function, but the data on the regulatory role of Type I MTases are limited. We have previously characterized molecularly a non-canonical Type I RM system, NgoAV, with phase-variable specificity, encoded by Neisseria gonorrhoeae FA1090. In the current work, we have investigated the impact of methyltransferase NgoAV (M.NgoAV) activity on gonococcal phenotype and on epigenetic control of gene expression. For this purpose, we have constructed and studied genetic variants (concerning activity and specificity) within M.NgoAV locus. Deletion of M.NgoAV or switch of its specificity had an impact on phenotype of N. gonorrhoeae. Biofilm formation and planktonic growth, the resistance to antibiotics, which target bacterial peptidoglycan or other antimicrobials, and invasion of human epithelial host cells were affected. The expression of genes was deregulated in gonococcal cells with knockout M.NgoAV gene and the variant with new specificity. For the first time, the existence of a phasevarion (phase-variable regulon), directed by phase-variable Type I MTase, is demonstrated.

An Ancient Respiratory System in the Widespread Sedimentary Archaea Thermoprofundales

Thermoprofundales, formerly Marine Benthic Group D (MBG-D), is a ubiquitous archaeal lineage found in sedimentary environments worldwide. However, its taxonomic classification, metabolic pathways, and evolutionary history are largely unexplored because of its uncultivability and limited number of sequenced genomes. In this study, phylogenomic analysis and average amino acid identity values of a collection of 146 Thermoprofundales genomes revealed five Thermoprofundales subgroups (A-E) with distinct habitat preferences. Most of the microorganisms from Subgroups B and D were thermophiles inhabiting hydrothermal vents and hot spring sediments, whereas those from Subgroup E were adapted to surface environments where sunlight is available. H2 production may be featured in Thermoprofundales as evidenced by a gene cluster encoding the ancient membrane-bound hydrogenase (MBH) complex. Interestingly, a unique structure separating the MBH gene cluster into two modular units was observed exclusively in the genomes of Subgroup E, which included a peripheral arm encoding the [NiFe] hydrogenase domain and a membrane arm encoding the Na+/H+ antiporter domain. These two modular structures were confirmed to function independently by detecting the H2-evolving activity in vitro and salt tolerance to 0.2 M NaCl in vivo, respectively. The peripheral arm of Subgroup E resembles the proposed common ancestral respiratory complex of modern respiratory systems, which plays a key role in the early evolution of life. In addition, molecular dating analysis revealed that Thermoprofundales is an early emerging archaeal lineage among the extant MBH-containing microorganisms, indicating new insights into the evolution of this ubiquitous archaea lineage.

Draft Genomes Sequences of 11 Geodermatophilaceae Strains Isolated from Building Stones from New England and Indian Stone Ruins found at historic sites in Tamil Nadu, India

Metagenomic analysis of stone microbiome from samples collected in New England, USA and Tamil Nadu, India identified numerous Actinobacteria including Geodermatphilaceae. A culture-dependent approach was performed as a companion study with this culture-independent metagenomic analysis of these stone samples and resulted in the isolation of eleven Geodermatphilaceae strains (2 Geodermatophilus and 9 Blastococcus strains). The genomes of the 11 Geodermatphilaceae strains were sequenced and analyzed. The genomes for the two Geodermatophilus isolates, DF1-2 and TF2-6, were 4.45 and 4.75 Mb, respectively, while the Blastococcus genomes ranged in size from 3.98 to 5.48 Mb. Phylogenetic analysis, digital DNA:DNA hybridization (dDDH), and comparisons of the average nucleotide identities (ANI) suggest the isolates represent novel Geodermatophilus and Blastococcus species. Functional analysis of the Geodermatphilaceae genomes provides insight on the stone microbiome niche.

Complete Genome Sequence of Flavobacterium sediminilitoris YSM-43 T , Isolated from Tidal Sediment in Yeosu

Here, we report the complete genome sequence of Flavobacterium sediminilitoris YSM-43^T, isolated from a tidal flat in Yeosu, Republic of Korea. The whole genome consists of one circular chromosome of 3,913,692 bp. A total of 3,599 genes were predicted, comprising 3,537 coding DNA sequences (CDSs), 50 tRNAs, 9 rRNAs, and 3 noncoding RNAs (ncRNAs).

Genomic sequence of the non-pathogen Neisseria sp. strain MA1-1 with antibiotic resistance and virulence factors isolated from a head and neck cancer patient

Recent research has claimed virulence factors or antimicrobial resistance in commensal or non-pathogenic Neisseria spp. This study aimed to isolate and analyze commensal microorganisms related to the genus Neisseria from the oral cavity of a patient with head and neck cancer. We successfully isolated strain MA1-1 and identified its functional gene contents. Although strain MA1-1 was related to Neisseria flava based on 16S rRNA gene sequence similarity, genomic relatedness analysis revealed that strain MA1-1 was closely related to Neisseria mucosa, reported as a commensal Neisseria species. The strain MA1-1 genome harbored genes for microaerobic respiration and the complete core metabolic pathway with few transporters for nutrients. A number of genes have been associated with virulence factors and resistance to various antibiotics. In addition, the comparative genomic analysis showed that most genes identified in the strain MA1-1 were shared with other Neisseria spp. including two well-known pathogens, Neisseria gonorrhoeae and Neisseria meningitidis. This indicates that the gene content of intra-members of the genus Neisseria has been evolutionarily conserved and is stable, with no gene recombination with other microbes in the host. Finally, this study provides more fundamental interpretations for the complete gene sequence of commensal Neisseria spp. and will contribute to advancing public health knowledge.

Infant Gut Microbial Metagenome Mining of α-l-Fucosidases with Activity on Fucosylated Human Milk Oligosaccharides and Glycoconjugates

The gastrointestinal microbiota members produce α-l-fucosidases that play key roles in mucosal, human milk, and dietary oligosaccharide assimilation. Here, 36 open reading frames (ORFs) coding for putative α-l-fucosidases belonging to glycosyl hydrolase family 29 (GH29) were identified through metagenome analysis of breast-fed infant fecal microbiome. Twenty-two of those ORFs showed a complete coding sequence with deduced amino acid sequences displaying the highest degree of identity with α-l-fucosidases from Bacteroides thetaiotaomicron, Bacteroides caccae, Phocaeicola vulgatus, Phocaeicola dorei, Ruminococcus gnavus, and Streptococcus parasanguinis. Based on sequence homology, 10 α-l-fucosidase genes were selected for substrate specificity characterization. The α-l-fucosidases Fuc18, Fuc19A, Fuc35B, Fuc39, and Fuc1584 showed hydrolytic activity on α1,3/4-linked fucose present in Lewis blood antigens and the human milk oligosaccharide (HMO) 3-fucosyllactose. In addition, Fuc1584 also hydrolyzed fucosyl-α-1,6-N-acetylglucosamine (6FN), a component of the core fucosylation of N-glycans. Fuc35A and Fuc193 showed activity on α1,2/3/4/6 linkages from H type-2, Lewis blood antigens, HMOs and 6FN. Fuc30 displayed activity only on α1,6-linked l-fucose, and Fuc5372 showed a preference for α1,2 linkages. Fuc2358 exhibited a broad substrate specificity releasing l-fucose from all the tested free histo-blood group antigens, HMOs, and 6FN. This latest enzyme also displayed activity in glycoconjugates carrying lacto-N-fucopentaose II (Lea) and lacto-N-fucopentaose III (Lex) and in the glycoprotein mucin. Fuc18, Fuc19A, and Fuc39 also removed l-fucose from neoglycoproteins and human α-1 acid glycoprotein. These results give insight into the great diversity of α-l-fucosidases from the infant gut microbiota, thus supporting the hypothesis that fucosylated glycans are crucial for shaping the newborn microbiota composition. IMPORTANCE α-l-Fucosyl residues are frequently present in many relevant glycans, such as human milk oligosaccharides (HMOs), histo-blood group antigens (HBGAs), and epitopes on cell surface glycoconjugate receptors. These fucosylated glycans are involved in a number of mammalian physiological processes, including adhesion of pathogens and immune responses. The modulation of l-fucose content in such processes may provide new insights and knowledge regarding molecular interactions and may help to devise new therapeutic strategies. Microbial α-l-fucosidases are exoglycosidases that remove α-l-fucosyl residues from free oligosaccharides and glycoconjugates and can be also used in transglycosylation reactions to synthesize oligosaccharides. In this work, α-l-fucosidases from the GH29 family were identified and characterized from the metagenome of fecal samples of breastfed infants. These enzymes showed different substrate specificities toward HMOs, HBGAs, naturally occurring glycoproteins, and neoglycoproteins. These novel glycosidase enzymes from the breast-fed infant gut microbiota, which resulted in a good source of α-l-fucosidases, have great biotechnological potential.

Analysis of lineage-specific protein family variability in prokaryotes combined with evolutionary reconstructions

Background

Evolutionary rate is a key characteristic of gene families that is linked to the functional importance of the respective genes as well as specific biological functions of the proteins they encode. Accurate estimation of evolutionary rates is a challenging task that requires precise phylogenetic analysis. Here we present an easy to estimate protein family level measure of sequence variability based on alignment column homogeneity in multiple alignments of protein sequences from Clade-Specific Clusters of Orthologous Genes (csCOGs).

Results

We report genome-wide estimates of variability for 8 diverse groups of bacteria and archaea and investigate the connection between variability and various genomic and biological features. The variability estimates are based on homogeneity distributions across amino acid sequence alignments and can be obtained for multiple groups of genomes at minimal computational expense. About half of the variance in variability values can be explained by the analyzed features, with the greatest contribution coming from the extent of gene paralogy in the given csCOG. The correlation between variability and paralogy appears to originate, primarily, not from gene duplication, but from acquisition of distant paralogs and xenologs, introducing sequence variants that are more divergent than those that could have evolved in situ during the lifetime of the given group of organisms. Both high-variability and low-variability csCOGs were identified in all functional categories, but as expected, proteins encoded by integrated mobile elements as well as proteins involved in defense functions and cell motility are, on average, more variable than proteins with housekeeping functions. Additionally, using linear discriminant analysis, we found that variability and fraction of genomes carrying a given gene are the two variables that provide the best prediction of gene essentiality as compared to the results of transposon mutagenesis in Sulfolobus islandicus.

Conclusions

Variability, a measure of sequence diversity within an alignment relative to the overall diversity within a group of organisms, offers a convenient proxy for evolutionary rate estimates and is informative with respect to prediction of functional properties of proteins. In particular, variability is a strong predictor of gene essentiality for the respective organisms and indicative of sub- or neofunctionalization of paralogs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13062-022-00337-7.

Acaricidal Activity and Field Efficacy Analysis of the Potential Biocontrol Agent Bacillus vallismortis NBIF-001 against Spider Mites

In recent years, spider mites have caused considerable economic losses to global agriculture. However, currently available management strategies are limited because of the rapid development of resistance. In this study, Bacillus vallismortis NBIF-001 was isolated and evaluated for its acaricidal activity. NBIF-001 exhibited a significant lethal effect on spider mites within 48 h. The median lethal concentration (LC₅₀) of the culture powders (3.2 × 10¹⁰ CFU/g) was 50.2 µg/mL for Tetranychus urticae (red form), 18.0 µg/mL for T. urticae (green form), and 15.7 µg/mL for Panonychus citri (McGregor). Cultivation optimisation experiments showed that when the number of spores increased, fermentation toxicity also increased. Moreover, field experiments demonstrated that NBIF-001 performed well in the biocontrol of P. citri, which showed a similar corrected field efficacy with the chemical control (67.1 ± 7.9% and 71.1 ± 6.4% after 14 days). Genomics analysis showed that NBIF-001 contains 231 factors and seven gene clusters of metabolites that may be involved in its acaricidal activity. Further bioassays of the fermentation supernatants showed that 50× dilution treatments killed 72.5 ± 5.4% of the mites in 48 h, which was similar with those of the broth. Bioassays of the supernatant proteins confirmed that various proteins exhibited acaricidal activity. Five candidate proteins were expressed and purified successfully. The bioassays showed that the small protein BVP8 exhibited significant acaricidal activity with an LC₅₀ of 12.4 μg/mL (T. urticae). Overall, these findings suggest that B. vallismortis NBIF-001 is a potential biocontrol agent for spider mite management.

Interaction between Plasma Metabolomics and Intestinal Microbiome in db/db Mouse, an Animal Model for Study of Type 2 Diabetes and Diabetic Kidney Disease

Evidence has demonstrated that either metabolites or intestinal microbiota are involved in the pathogenesis of type 2 diabetes (T2D) and diabetic kidney disease (DKD). To explore the interaction between plasma metabolomics and intestinal microbiome in the progress of T2D-DKD, in the current study, we analyzed metabolomics in the plasma of db/db mice with liquid chromatography–mass spectrometry and also examined intestinal prokaryotes and entire gut microbiome dysbiosis at the genus level with both 16S rDNA and metagenomic sequencing techniques. We found that Negativibacillus and Rikenella were upregulated, while Akkermansia, Candidatus, Erysipelatoclostridium and Ileibacterium were downregulated in the colon of db/db mice compared with non-diabetic controls. In parallel, a total of 91 metabolites were upregulated, while 23 were downregulated in the plasma of db/db mice. The top five upregulated metabolites included D-arabinose 5-phosphate, estrone 3-sulfate, L-theanine, 3′-aenylic acid and adenosine 5′-monophosphate, and the five most significantly downregulated metabolites were aurohyocholic acid sodium salt, calcium phosphorylcholine chloride, tauro-alpha-muricholic acid sodium salt, galactinol and phosphocholine. These plasma metabolites were interacted with intestinal microbiomes, which are mainly involved in the pathways related to the biosynthesis of unsaturated fatty acids, fatty acid elongation, steroid biosynthesis, and D-arginine and D-ornithine metabolism. In the differential metabolites, N-acetyl-L-ornithine, ornithine and L-kyn could be metabolized by the correspondingly differential ontology genes in the intestinal metagenome. The current study thereby provides evidence for a gut–metabolism–kidney axis in the metabolism of db/db mice, in which the gut microbiome and circulating metabolomics interact, and suggests that information from this axis may contribute to our understanding of T2D and DKD pathogenesis.

Genome-Wide Characterization Reveals Variation Potentially Involved in Pathogenicity and Mycotoxins Biosynthesis of Fusarium proliferatum Causing Spikelet Rot Disease in Rice

Fusarium proliferatum is the primary cause of spikelet rot disease in rice (Oryza sativa L.) in China. The pathogen not only infects a wide range of cereals, causing severe yield losses but also contaminates grains by producing various mycotoxins that are hazardous to humans and animals. Here, we firstly reported the whole-genome sequence of F. proliferatum strain Fp9 isolated from the rice spikelet. The genome was approximately 43.9 Mb with an average GC content of 48.28%, and it was assembled into 12 scaffolds with an N₅₀ length of 4,402,342 bp. There is a close phylogenetic relationship between F. proliferatum and Fusarium fujikuroi, the causal agent of the bakanae disease of rice. The expansion of genes encoding cell wall-degrading enzymes and major facilitator superfamily (MFS) transporters was observed in F. proliferatum relative to other fungi with different nutritional lifestyles. Species-specific genes responsible for mycotoxins biosynthesis were identified among F. proliferatum and other Fusarium species. The expanded and unique genes were supposed to promote F. proliferatum adaptation and the rapid response to the host's infection. The high-quality genome of F. proliferatum strain Fp9 provides a valuable resource for deciphering the mechanisms of pathogenicity and secondary metabolism, and therefore shed light on development of the disease management strategies and detoxification of mycotoxins contamination for spikelet rot disease in rice.

Identification and whole-genome sequencing analysis of Vibrio vulnificus strains causing pearl gentian grouper disease in China

Vibrio vulnificus is a pathogenic bacterium that causes disease in marine fish, affecting fish farming and human health worldwide. In May 2021, in the Bohai Bay region, a disease broke out in commercially farmed pearl gentian grouper (♀Epinephelus fuscoguttatus × ♂Epinephelus lanceolatus), causing huge economic losses. The diseased fish had skin lesions, water accumulation in their abdomens, and showed tissue and organ damage. V. vulnificus biotype 2 has been reported in eels and other marine fish, but it is less reported in pearl gentian grouper. In this study, the pathogenic strain isolated from diseased fish was identified as V. vulnificus EPL 0201 biotype 2 on the basis of physiological and biochemical characteristics and the results of 16S rRNA gene and gyrB sequencing, virulence gene detection, and recursive infection experiments. To gain a comprehensive understanding of the pathogenicity and drug resistance of this strain, whole-genome sequencing was performed. Whole-genome analysis showed that the gene map of this strain was complete. The Virulence Factor Database annotation results showed that this strain had the key virulence factor genes vvhA and rtxA, which cause host disease. In addition, this strain had genes conferring resistance against cephalosporins, aminoglycosides, tetracyclines, and sulfonamides. Antimicrobial susceptibility testing confirmed the presence of these resistance genes identified in the genome. The results of this study show that V. vulnificus EPL 0201 biotype 2 is a multi-drug resistant strain with high pathogenicity.

The survivor strain: isolation and characterization of Phormidium yuhuli AB48, a filamentous phototactic cyanobacterium with biotechnological potential

Despite their recognized potential, current applications of cyanobacteria as microbial cell factories remain in early stages of development. This is partly due to the fact that engineered strains are often difficult to grow at scale. This technical challenge contrasts with the dense and highly productive cyanobacteria populations thriving in many natural environments. It has been proposed that the selection of strains pre-adapted for growth in industrial photobioreactors could enable more productive cultivation outcomes. Here, we described the initial morphological, physiological, and genomic characterization of Phormidium yuhuli AB48 isolated from an industrial photobioreactor environment. P. yuhuli AB48 is a filamentous phototactic cyanobacterium with a growth rate comparable to Synechocystis sp. PCC 6803. The isolate forms dense biofilms under high salinity and alkaline conditions and manifests a similar nutrient profile to Arthrospira platensis (Spirulina). We sequenced, assembled, and analyzed the P. yuhuli AB48 genome, the first closed circular isolate reference genome for a member of the Phormidium genus. We then used cultivation experiments in combination with proteomics and metabolomics to investigate growth characteristics and phenotypes related to industrial scale cultivation, including nitrogen and carbon utilization, salinity, and pH acclimation, as well as antibiotic resistance. These analyses provide insight into the biological mechanisms behind the desirable growth properties manifested by P. yuhuli AB48 and position it as a promising microbial cell factory for industrial-scale bioproduction[221, 1631].

Functional characterization of prokaryotic dark matter: the road so far and what lies ahead

Eight-hundred thousand to one trillion prokaryotic species may inhabit our planet. Yet, fewer than two-hundred thousand prokaryotic species have been described. This uncharted fraction of microbial diversity, and its undisclosed coding potential, is known as the "microbial dark matter" (MDM). Next-generation sequencing has allowed to collect a massive amount of genome sequence data, leading to unprecedented advances in the field of genomics. Still, harnessing new functional information from the genomes of uncultured prokaryotes is often limited by standard classification methods. These methods often rely on sequence similarity searches against reference genomes from cultured species. This hinders the discovery of unique genetic elements that are missing from the cultivated realm. It also contributes to the accumulation of prokaryotic gene products of unknown function among public sequence data repositories, highlighting the need for new approaches for sequencing data analysis and classification. Increasing evidence indicates that these proteins of unknown function might be a treasure trove of biotechnological potential. Here, we outline the challenges, opportunities, and the potential hidden within the functional dark matter (FDM) of prokaryotes. We also discuss the pitfalls surrounding molecular and computational approaches currently used to probe these uncharted waters, and discuss future opportunities for research and applications.

Molecular adaptations during viral epidemics

In 1977, the world witnessed both the eradication of smallpox and the beginning of the modern age of genomics. Over the following half-century, 7 epidemic viruses of international concern galvanized virologists across the globe and led to increasingly extensive virus genome sequencing. These sequencing efforts exerted over periods of rapid adaptation of viruses to new hosts, in particular, humans provide insight into the molecular mechanisms underpinning virus evolution. Investment in virus genome sequencing was dramatically increased by the unprecedented support for phylogenomic analyses during the COVID-19 pandemic. In this review, we attempt to piece together comprehensive molecular histories of the adaptation of variola virus, HIV-1 M, SARS, H1N1-SIV, MERS, Ebola, Zika, and SARS-CoV-2 to the human host. Disruption of genes involved in virus-host interaction in animal hosts, recombination including genome segment reassortment, and adaptive mutations leading to amino acid replacements in virus proteins involved in host receptor binding and membrane fusion are identified as the key factors in the evolution of epidemic viruses.

Features of smaller ribosomes in candidate phyla radiation (CPR) bacteria revealed with a molecular evolutionary analysis

The candidate phyla radiation (CPR) is a large bacterial group consisting mainly of uncultured lineages. They have small cells and small genomes, and they often lack ribosomal proteins uL1, bL9, and/or uL30, which are basically ubiquitous in non-CPR bacteria. Here, we comprehensively analyzed the genomic information on CPR bacteria and identified their unique properties. The distribution of protein lengths in CPR bacteria peaks at around 100-150 amino acids, whereas the position of the peak varies in the range of 100-300 amino acids in free-living non-CPR bacteria, and at around 100-200 amino acids in most symbiotic non-CPR bacteria. These results show that the proteins of CPR bacteria are smaller, on average, than those of free-living non-CPR bacteria, like those of symbiotic non-CPR bacteria. We found that ribosomal proteins bL28, uL29, bL32, and bL33 have been lost in CPR bacteria in a taxonomic lineage-specific manner. Moreover, the sequences of approximately half of all ribosomal proteins of CPR differ, in part, from those of non-CPR bacteria, with missing regions or specifically added regions. We also found that several regions in the 16S, 23S, and 5S rRNAs of CPR bacteria are lacking, which presumably caused the total predicted lengths of the three rRNAs of CPR bacteria to be smaller than those of non-CPR bacteria. The regions missing in the CPR ribosomal proteins and rRNAs are located near the surface of the ribosome, and some are close to one another. These observations suggest that ribosomes are smaller in CPR bacteria than those in free-living non-CPR bacteria, with simplified surface structures.

Analysis of Effector and Immunity Proteins of the GA2 Type VI Secretion Systems of Gut Bacteroidales

Three distinct genetic architectures (GAs) of Type VI secretion systems (T6SSs) have been described in gut Bacteroidales species, each with unique genes and characteristics. Unlike the GA3 T6SSs, potent antagonism has not yet been demonstrated for the GA1 or GA2 T6SSs. We previously showed that the GA2 T6SS loci are contained on integrative and conjugative elements and that there are five subtypes. Collectively, GA2 are the most prevalent Bacteroidales T6SSs in the human populations analyzed. In this study, we provide a comprehensive bioinformatic analysis of the three variable regions of GA2 T6SS loci, which encode toxic effector and immunity proteins. In total, we identified 63 distinct effectors encoded within 31 nonredundant GA2 loci, 18 of which do not have described motifs or predicted functions. We provide experimental evidence for toxin activity for four different GA2 effectors, showing that each functions only when present in the periplasm, and experimentally confirm their cognate immunity proteins. Our data demonstrate that each GA2 locus encodes at least three distinct effectors with targets in both the cytoplasm and the periplasm. The data also suggest that the effectors of a given locus are loaded onto the tube by different mechanisms, which may allow all three effectors encoded within a single GA2 locus with distinct antibacterial activity to be loaded onto a single T6 tube, increasing the antagonistic effect. IMPORTANCE Humans are colonized with many gut Bacteroidales species at high density, allowing for extensive opportunities for contact-dependent antagonism. To begin to understand the antagonistic potential of the GA2 T6SSs of the gut Bacteroidales, we performed bioinformatic and experimental analyses of the three divergent regions containing the toxin effector and immunity genes. We show that each GA2 T6SS locus encodes at least three distinct toxic effectors including toxins linked to Rhs and Hcp with cytoplasmic targets, and unlinked effectors with targets in the periplasm. The diversity and modality of effectors exceeds that of the GA1 or GA3 T6SS loci (M. J. Coyne, K. G. Roelofs, and L. E. Comstock, BMC Genomics 17:58, 2016, https://doi.org/10.1186/s12864-016-2377-z) and suggests that these T6SSs have the potential to be potent antibacterial weapons in the human gut.

Identification of a Putative CodY Regulon in the Gram-Negative Phylum Synergistetes

CodY is a dominant regulator in low G + C, Gram-positive Firmicutes that governs the regulation of various metabolic pathways and cellular processes. By using various bioinformatics analyses and DNA affinity precipitation assay (DAPA), this study confirmed the presence of CodY orthologues and corresponding regulons in Gram-negative Synergistetes. A novel palindromic sequence consisting of AT-rich arms separated by a spacer region of variable length and sequence was identified in the promoters of the putative codY-containing operons in Synergistetes. The consensus sequence from genera Synergistes and Cloacibacillus (5′-AATTTTCTTAAAATTTCSCTTGATATTTACAATTTT) contained three AT-rich regions, resulting in two palindromic sequences; one of which is identical to Firmicutes CodY box (5′-AATTTTCWGAAAATT). The function of the consensus sequence was tested by using a recombinant CodY protein (His-CodY_DSM) of Cloacibacillus evryensis DSM19522 in DAPA. Mutations in the central AT-rich sequence reduced significantly the binding of His-CodY_DSM, whereas mutations in the 5′ or 3′ end AT-rich sequence slightly reduced the binding, indicating that CodY_DSM could recognize both palindromic sequences. The proposed binding sequences were found in the promoters of multiple genes involved in amino acids biosynthesis, metabolism, regulation, and stress responses in Synergistetes. Thus, a CodY-like protein from Synergistetes may function similarly to Firmicutes CodY.

Genomic Insights Into the Interspecific Diversity and Evolution of Mobiluncus, a Pathogen Associated With Bacterial Vaginosis

Bacterial vaginosis (BV) is a common vaginal infection and has been associated with increased risk for a wide array of health issues. BV is linked with a variety of heterogeneous pathogenic anaerobic bacteria, among which Mobiluncus is strongly associated with BV diagnosis. However, their genetic features, pathogenicity, interspecific diversity, and evolutionary characters have not been illustrated at genomic level. The current study performed phylogenomic and comparative genomic analyses of Mobiluncus. Phylogenomic analyses revealed remarkable phylogenetic distinctions among different species. Compared with M. curtisii, M. mulieris had a larger genome and pangenome size with more insertion sequences but less CRISPR-Cas systems. In addition, these two species were diverse in profile of virulence factors, but harbored similar antibiotic resistance genes. Statistically different functional genome profiles between strains from the two species were determined, as well as correlations of some functional genes/pathways with putative pathogenicity. We also showed that high levels of horizontal gene transfer might be an important strategy for species diversification and pathogenicity. Collectively, this study provides the first genome sequence level description of Mobiluncus, and may shed light on its virulence/pathogenicity, functional diversification, and evolutionary dynamics. Our study could facilitate the further investigations of this important pathogen, and might improve the future treatment of BV.

EspF of Enterohemorrhagic Escherichia coli Enhances Apoptosis via Endoplasmic Reticulum Stress in Intestinal Epithelial Cells: An Isobaric Tags for Relative and Absolute Quantitation-Based Comparative Proteomic Analysis

There have been large foodborne outbreaks related to Enterohemorrhagic Escherichia coli (EHEC) around the world. Among its virulence proteins, the EspF encoded by locus of enterocyte effacement is one of the most known functional effector proteins. In this research, we infected the HT-29 cells with the EHEC wild type strain and EspF-deficient EHEC strain. Via the emerging technique isobaric tags for relative and absolute quantitation (iTRAQ), we explored the pathogenic characteristics of EspF within host cells. Our data showed that the differences regarding cellular responses mainly contained immune regulation, protein synthesis, signal transduction, cellular assembly and organization, endoplasmic reticulum (ER) stress, and apoptosis. Notably, compared with the EspF-deficient strain, the protein processing in the ER and ribosome were upregulated during wild type (WT) infection. Our findings proved that the EspF of Enterohemorrhagic Escherichia coli induced ER stress in intestinal epithelial cells; the ER stress-dependent apoptosis pathway was also activated within the host cells. This study provides insight into the virulence mechanism of protein EspF, which will deepen our general understanding of A/E pathogens and their interaction with host proteins.

Genomic and Chemical Profiling of B9, a Unique Penicillium Fungus Derived from Sponge

This study presented the first insights into the genomic and chemical profiles of B9, a specific Penicillium strain derived from sponges of the South China Sea that demonstrated the closest morphological and phylogenetic affinity to P. paxillin. Via the Illumina MiSeq sequencing platform, the draft genome was sequenced, along with structural assembly and functional annotation. There were 34 biosynthetic gene clusters (BGCs) predicted against the antiSMASH database, but only 4 gene clusters could be allocated to known BGCs (≥50% identities). Meanwhile, the comparison between B9 and P. paxillin ATCC 10480 demonstrated clear distinctions in morphology, which might be ascribed to the unique environmental adaptability of marine endosymbionts. In addition, two novel pyridinones, penicidihydropyridone A (2) and penicidihydropyridone B (3), were isolated from cultures of B9, and structurally characterized by nuclear magnetic resonance (NMR) and mass spectrometry (MS). The absolute configurations were confirmed by comparison of experimental and calculated electronic circular dichroism (ECD) curves. In addition, structure-based molecular docking indicated that both neo-pyridinones might block the programmed cell death protein 1(PD-1) pathway by competitively binding a programmed cell death 1 ligand 1(PD-L1) dimer. This was verified by the significant inhibition rates of the PD-1/L1 interaction. These indicated that Penicillium sp. B9 possessed a potential source of active secondary metabolites.

Candidatus Nitrosopolaris, a genus of putative ammonia-oxidizing archaea with a polar/alpine distribution

Ammonia-oxidizing archaea (AOA) are key players in the nitrogen cycle of polar soils. Here, we analyzed metagenomic data from tundra soils in Rásttigáisá, Norway, and recovered four metagenome-assembled genomes (MAGs) assigned to the genus 'UBA10452', an uncultured lineage of putative AOA in the order Nitrososphaerales ('terrestrial group I.1b'), phylum Thaumarchaeota. Analysis of other eight previously reported MAGs and publicly available amplicon sequencing data revealed that the UBA10452 lineage is predominantly found in acidic polar and alpine soils. In particular, UBA10452 MAGs were more abundant in highly oligotrophic environments such as mineral permafrost than in more nutrient-rich, vegetated tundra soils. UBA10452 MAGs harbour multiple copies of genes related to cold tolerance, particularly genes involved in DNA replication and repair. Based on the phylogenetic, biogeographic, and ecological characteristics of 12 UBA10452 MAGs, which include a high-quality MAG (90.8% complete, 3.9% redundant) with a nearly complete 16S rRNA gene, we propose a novel Candidatus genus, Ca. Nitrosopolaris, with four species representing clear biogeographic/habitat clusters.

Transcriptome analysis reveals regulation mechanism of methyl jasmonate-induced terpenes biosynthesis in Curcuma wenyujin

Curcuma wenyujin is the source plant of three traditional Chinese medicines, which have been widely used in clinical treatment over 1000 years. The content of terpenes, the major medicinal active ingredients, is relatively low in this plant. Studies have shown that MeJA can promote terpenes biosynthesis in plants. However, the mechanism underlying the effect of MeJA in C. wenyujin remains unclear. In this work, the transcriptome of C. wenyujin leaves with MeJA treatment was analyzed to elucidate the regulation mechanism of MeJA-mediated terpene biosynthesis. Based on the RNA-seq data, 7,246 unigenes were differentially expressed with MeJA treatment. Expression pattern clustering of DEGs revealed that unigenes, related to JA biosynthesis and signal transduction, responded to exogenous MeJA stimulation on the early stage and maintained throughout the process. Subsequently, unigenes related to terpene biosynthesis pathway showed a significant up-regulation with 6 h treatment. The analysis results suggested that MeJA induced the expression of JA biosynthesis genes (such as LOXs, AOSs, AOCs, OPRs, and MFPs) and JA signal transduction core genes (JAZs and MYCs) to activate JA signaling pathway. Meanwhile, downstream JA-responsive genes presented up-regulated expression levels such as AACT, HMGSs, HMGRs, DXSs, DXRs, MCTs, HDSs, and HDRs, thus promoting terpenes biosynthesis. The transcriptional expressions of these genes were validated by qRT-PCR. In addition, six CwTPS genes in response to MeJA were identified. With MeJA treatment, the expression levels of CwTPSs were increased as well as those of the transcription factors MYB, NAC, bZIP, WRKY, AP2/ERF, and HLH. These TFs might potentially regulate terpenes biosynthesis. These results provide insights for regulation mechanism of terpenes biosynthesis.

To kill or to be killed: pangenome analysis of Escherichia coli strains reveals a tailocin specific for pandemic ST131

BACKGROUND

Escherichia coli (E. coli) has been one of the most studied model organisms in the history of life sciences. Initially thought just to be commensal bacteria, E. coli has shown wide phenotypic diversity including pathogenic isolates with great relevance to public health. Though pangenome analysis has been attempted several times, there is no systematic functional characterization of the E. coli subgroups according to the gene profile.

RESULTS

Systematically scanning for optimal parametrization, we have built the E. coli pangenome from 1324 complete genomes. The pangenome size is estimated to be ~25,000 gene families (GFs). Whereas the core genome diminishes as more genomes are added, the softcore genome (≥95% of strains) is stable with ~3000 GFs regardless of the total number of genomes. Apparently, the softcore genome (with a 92% or 95% generation threshold) can define the genome of a bacterial species listing the critically relevant, evolutionarily most conserved or important classes of GFs. Unsupervised clustering of common E. coli sequence types using the presence/absence GF matrix reveals distinct characteristics of E. coli phylogroups B1, B2, and E. We highlight the bi-lineage nature of B1, the variation of the secretion and of the iron acquisition systems in ST11 (E), and the incorporation of a highly conserved prophage into the genome of ST131 (B2). The tail structure of the prophage is evolutionarily related to R2-pyocin (a tailocin) from Pseudomonas aeruginosa PAO1. We hypothesize that this molecular machinery is highly likely to play an important role in protecting its own colonies; thus, contributing towards the rapid rise of pandemic E. coli ST131.

CONCLUSIONS

This study has explored the optimized pangenome development in E. coli. We provide complete GF lists and the pangenome matrix as supplementary data for further studies. We identified biological characteristics of different E. coli subtypes, specifically for phylogroups B1, B2, and E. We found an operon-like genome region coding for a tailocin specific for ST131 strains. The latter is a potential killer weapon providing pandemic E. coli ST131 with an advantage in inter-bacterial competition and, suggestively, explains their dominance as human pathogen among E. coli strains.

Staphylococcal saoABC Operon Codes for a DNA-Binding Protein SaoC Implicated in the Response to Nutrient Deficit

Whilst a large number of regulatory mechanisms for gene expression have been characterised to date, transcription regulation in bacteria still remains an open subject. In clinically relevant and opportunistic pathogens, such as Staphylococcus aureus, transcription regulation is of great importance for host-pathogen interactions. In our study we investigated an operon, exclusive to staphylococci, that we name saoABC. We showed that SaoC binds to a conserved sequence motif present upstream of the saoC gene, which likely provides a negative feedback loop. We have also demonstrated that S. aureus ΔsaoB and ΔsaoC mutants display altered growth dynamics in non-optimal media; ΔsaoC exhibits decreased intracellular survival in human dermal fibroblasts, whereas ΔsaoB produces an elevated number of persisters, which is also elicited by inducible production of SaoC in ΔsaoBΔsaoC double mutant. Moreover, we have observed changes in the expression of saoABC operon genes during either depletion of the preferential carbon or the amino acid source as well as during acidification. Comparative RNA-Seq of the wild type and ΔsaoC mutant demonstrated that SaoC influences transcription of genes involved in amino acid transport and metabolism, and notably of those coding for virulence factors. Our results suggest compellingly that saoABC operon codes for a DNA-binding protein SaoC, a novel staphylococcal transcription factor, and its antagonist SaoB. We linked SaoC to the response to nutrient deficiency, a stress that has a great impact on host-pathogen interactions. That impact manifests in SaoC influence on persister formation and survival during internalisation to host cells, as well as on the expression of genes of virulence factors that may potentially result in profound alternations in the pathogenic phenotype. Investigation of such novel regulatory mechanisms is crucial for our understanding of the dynamics of interactions between pathogenic bacteria and host cells, particularly in the case of clinically relevant, opportunistic pathogens such as Staphylococcus aureus.

Soil microbial diversity and functional profiling of a tropical rainforest of a highly dissected low hill from the upper Itaya river basin revealed by analysis of shotgun metagenomics sequencing data

The tropical rainforest of a highly dissected low hill from the upper Itaya river basin belongs to the western Amazonia region. Some investigations on the biodiversity of these rainforests were more focused on animals and plants diversity. The soils of this region are composed of moderately fertile sediments deposited recently from the initiation of the Andean orogenesis in the Miocene until now. However, scientific information about the soil microbial and functional diversity is still missing. This report presents shotgun metagenomics sequencing data from soils of this rainforest type. A composite loamy soil sample was collected from a primary forest, and metagenomic DNA was purified with standardized methods. Furthermore, libraries were prepared and paired-end sequenced on the Illumina NextSeq 550 platform. Raw Illumina paired-end reads have been uploaded and analysed in the Metagenomics RAST server (MG-RAST). The raw sequence data in fastq format is available at NCBI's Sequence Read Archive (SRA) with accession number SRX12846710.

Introduction to the principles and methods underlying the recovery of metagenome-assembled genomes from metagenomic data

The rise of metagenomics offers a leap forward for understanding the genetic diversity of microorganisms in many different complex environments by providing a platform that can identify potentially unlimited numbers of known and novel microorganisms. As such, it is impossible to imagine new major initiatives without metagenomics. Nevertheless, it represents a relatively new discipline with various levels of complexity and demands on bioinformatics. The underlying principles and methods used in metagenomics are often seen as common knowledge and often not detailed or fragmented. Therefore, we reviewed these to guide microbiologists in taking the first steps into metagenomics. We specifically focus on a workflow aimed at reconstructing individual genomes, that is, metagenome‐assembled genomes, integrating DNA sequencing, assembly, binning, identification and annotation.

UG/Abi: a highly diverse family of prokaryotic reverse transcriptases associated with defense functions

Reverse transcriptases (RTs) are enzymes capable of synthesizing DNA using RNA as a template. Within the last few years, a burst of research has led to the discovery of novel prokaryotic RTs with diverse antiviral properties, such as DRTs (Defense-associated RTs), which belong to the so-called group of unknown RTs (UG) and are closely related to the Abortive Infection system (Abi) RTs. In this work, we performed a systematic analysis of UG and Abi RTs, increasing the number of UG/Abi members up to 42 highly diverse groups, most of which are predicted to be functionally associated with other gene(s) or domain(s). Based on this information, we classified these systems into three major classes. In addition, we reveal that most of these groups are associated with defense functions and/or mobile genetic elements, and demonstrate the antiphage role of four novel groups. Besides, we highlight the presence of one of these systems in novel families of human gut viruses infecting members of the Bacteroidetes and Firmicutes phyla. This work lays the foundation for a comprehensive and unified understanding of these highly diverse RTs with enormous biotechnological potential.

The c-di-GMP Phosphodiesterase PipA (PA0285) Regulates Autoaggregation and Pf4 Bacteriophage Production in Pseudomonas aeruginosa PAO1

In Pseudomonas aeruginosa PAO1, 41 genes encode proteins predicted to be involved in the production or degradation of c-di-GMP, a ubiquitous secondary messenger that regulates a variety of physiological behaviors closely related to biofilm and aggregate formation. Despite extensive effort, the entire picture of this important signaling network is still unclear, with one-third of these proteins remaining uncharacterized. Here, we show that the deletion of pipA, which produces a protein containing two PAS domains upstream of a GGDEF-EAL tandem, significantly increased the intracellular c-di-GMP level and promoted the formation of aggregates both on surfaces and in planktonic cultures. However, this regulatory effect was not contributed by either of the two classic pathways modulating biofilm formation, exopolysaccharide (EPS) overproduction or motility inhibition. Transcriptome sequencing (RNA-Seq) data revealed that the expression levels of 361 genes were significantly altered in a ΔpipA mutant strain compared to the wild type (WT), indicating the critical role of PipA in PAO1. The most remarkably downregulated genes were located on the Pf4 bacteriophage gene cluster, which corresponded to a 2-log reduction in the Pf4 phage production in the ΔpipA mutant. The sizes of aggregates in ΔpipA cultures were affected by exogenously added Pf4 phage in a concentration-dependent manner, suggesting the quantity of phage plays a part in regulating the formation of aggregates. Further analysis demonstrated that PipA is highly conserved across 83 P. aeruginosa strains. Our work therefore for the first time showed that a c-di-GMP phosphodiesterase can regulate bacteriophage production and provided new insights into the relationship between bacteriophage and bacterial aggregation. IMPORTANCE The c-di-GMP signaling pathways in P. aeruginosa are highly organized and well coordinated, with different diguanylate cyclases and phosphodiesterases playing distinct roles in a complex network. Understanding the function of each enzyme and the underlying regulatory mechanisms not only is crucial for revealing how bacteria decide the transition between motile and sessile lifestyles, but also greatly facilitates the development of new antibiofilm strategies. This work identified bacteriophage production as a novel phenotypic output controlled transcriptionally by a phosphodiesterase, PipA. Further analysis suggested that the quantity of phage may be important in regulating autoaggregation, as either a lack of phage or overproduction was associated with higher levels of aggregation. Our study therefore extended the scope of c-di-GMP-controlled phenotypes and discovered a potential signaling circuit that can be target for biofilm treatment.

Conservation and Evolution of the Sporulation Gene Set in Diverse Members of the Firmicutes

The current classification of the phylum Firmicutes (new name, Bacillota) features eight distinct classes, six of which include known spore-forming bacteria. In Bacillus subtilis, sporulation involves up to 500 genes, many of which do not have orthologs in other bacilli and/or clostridia. Previous studies identified about 60 sporulation genes of B. subtilis that were shared by all spore-forming members of the Firmicutes. These genes are referred to as the sporulation core or signature, although many of these are also found in genomes of nonsporeformers. Using an expanded set of 180 firmicute genomes from 160 genera, including 76 spore-forming species, we investigated the conservation of the sporulation genes, in particular seeking to identify lineages that lack some of the genes from the conserved sporulation core. The results of this analysis confirmed that many small acid-soluble spore proteins (SASPs), spore coat proteins, and germination proteins, which were previously characterized in bacilli, are missing in spore-forming members of Clostridia and other classes of Firmicutes. A particularly dramatic loss of sporulation genes was observed in the spore-forming members of the families Planococcaceae and Erysipelotrichaceae. Fifteen species from diverse lineages were found to carry skin (sigK-interrupting) elements of different sizes that all encoded SpoIVCA-like recombinases but did not share any other genes. Phylogenetic trees built from concatenated alignments of sporulation proteins and ribosomal proteins showed similar topology, indicating an early origin and subsequent vertical inheritance of the sporulation genes. IMPORTANCE Many members of the phylum Firmicutes (Bacillota) are capable of producing endospores, which enhance the survival of important Gram-positive pathogens that cause such diseases as anthrax, botulism, colitis, gas gangrene, and tetanus. We show that the core set of sporulation genes, defined previously through genome comparisons of several bacilli and clostridia, is conserved in a wide variety of sporeformers from several distinct lineages of Firmicutes. We also detected widespread loss of sporulation genes in many organisms, particularly within the families Planococcaceae and Erysipelotrichaceae. Members of these families, such as Lysinibacillus sphaericus and Clostridium innocuum, could be excellent model organisms for studying sporulation mechanisms, such as engulfment, formation of the spore coat, and spore germination.

Global Phylogeny of Mycobacterium avium and Identification of Mutation Hotspots During Niche Adaptation

Mycobacterium avium is separated into four subspecies: M. avium subspecies avium (MAA), M. avium subspecies silvaticum (MAS), M. avium subspecies hominissuis (MAH), and M. avium subspecies paratuberculosis (MAP). Understanding the mechanisms of host and tissue adaptation leading to their clinical significance is vital to reduce the economic, welfare, and public health concerns associated with diseases they may cause in humans and animals. Despite substantial phenotypic diversity, the subspecies nomenclature is controversial due to high genetic similarity. Consequently, a set of 1,230 M. avium genomes was used to generate a phylogeny, investigate SNP hotspots, and identify subspecies-specific genes. Phylogeny reiterated the findings from previous work and established that Mycobacterium avium is a species made up of one highly diverse subspecies, known as MAH, and at least two clonal pathogens, named MAA and MAP. Pan-genomes identified coding sequences unique to each subspecies, and in conjunction with a mapping approach, mutation hotspot regions were revealed compared to the reference genomes for MAA, MAH, and MAP. These subspecies-specific genes may serve as valuable biomarkers, providing a deeper understanding of genetic differences between M. avium subspecies and the virulence mechanisms of mycobacteria. Furthermore, SNP analysis demonstrated common regions between subspecies that have undergone extensive mutations during niche adaptation. The findings provide insights into host and tissue specificity of this genetically conserved but phenotypically diverse species, with the potential to provide new diagnostic targets and epidemiological and therapeutic advances.