TREND: a platform for exploring protein function in prokaryotes based on phylogenetic, domain architecture and gene neighborhood analyses

Peroxisomal biogenesis factor 11 as a novel target to trigger lipid biosynthesis and salt stress resistance in oleaginous Tetradesmus obliquus

To overcome economic challenges in microalgal biofuel production, this study investigates the overexpression of peroxisome-localized peroxisomal biogenesis factor 11 (PEX11) to enhance lipid biosynthesis and improve salt stress resistance in Tetradesmus obliquus, aiming to advance microalgal biofuel production. Transgenic strains PEX11-2-1 and PEX11-2-2 exhibited a 2.13- and 2.51-fold increase in neutral lipid content and more cellular lipid droplets compared to WT, along with lipid yield and biomass escalating to 255.45 and 815.15 mg/L, respectively. This enhancement resulted from the redistribution of carbon precursors, increased intracellular reactive oxygen species, enhanced NADPH synthesis, and upregulation of lipid synthesis genes. Additionally, PEX11 improved salt stress tolerance by upregulating the expression of stress-responsive genes, including SnRK2 and PYRC. Fatty acid profile alterations, with increases in saturated fatty acids C16:0 and monounsaturated fatty acids C18:1, and decreases in polyunsaturated fatty acids, facilitated high-quality biofuel production. These findings highlight novel insights for advancing microalgae-based biorefinery.

The Psu protein of phage satellite P4 inhibits transcription termination factor ρ by forced hyper-oligomerization

Many bacteriophages modulate host transcription to favor expression of their own genomes. Phage satellite P4 polarity suppression protein, Psu, a building block of the viral capsid, inhibits hexameric transcription termination factor, ρ, by presently unknown mechanisms. Our cryogenic electron microscopy structures of ρ-Psu complexes show that Psu dimers clamp two inactive, open ρ rings and promote their expansion to higher-oligomeric states. ATPase, nucleotide binding and nucleic acid binding studies revealed that Psu hinders ρ ring closure and traps nucleotides in their binding pockets on ρ. Structure-guided mutagenesis in combination with growth, pull-down, and termination assays further delineated the functional ρ-Psu interfaces in vivo. Bioinformatic analyses revealed that Psu is associated with a wide variety of phage defense systems across Enterobacteriaceae, suggesting that Psu may regulate expression of anti-phage genes. Our findings show that modulation of the ρ oligomeric state via diverse strategies is a pervasive gene regulatory principle in bacteria.

FlhE functions as a chaperone to prevent formation of periplasmic flagella in Gram-negative bacteria

The bacterial flagellum, which facilitates motility, is composed of ~20 structural proteins organized into a long extracellular filament connected to a cytoplasmic rotor-stator complex via a periplasmic rod. Flagellum assembly is regulated by multiple checkpoints that ensure an ordered gene expression pattern coupled to the assembly of the various building blocks. Here, we use epifluorescence, super-resolution, and transmission electron microscopy to show that the absence of a periplasmic protein (FlhE) prevents proper flagellar morphogenesis and results in the formation of periplasmic flagella in Salmonella enterica. The periplasmic flagella disrupt cell wall synthesis, leading to a loss of normal cell morphology resulting in cell lysis. We propose that FlhE functions as a periplasmic chaperone to control assembly of the periplasmic rod, thus preventing formation of periplasmic flagella.

Ubiquitous purine sensor modulates diverse signal transduction pathways in bacteria

Purines and their derivatives control intracellular energy homeostasis and nucleotide synthesis, and act as signaling molecules. Here, we combine structural and sequence information to define a purine-binding motif that is present in sensor domains of thousands of bacterial receptors that modulate motility, gene expression, metabolism, and second-messenger turnover. Microcalorimetric titrations of selected sensor domains validate their ability to specifically bind purine derivatives, and evolutionary analyses indicate that purine sensors share a common ancestor with amino-acid receptors. Furthermore, we provide experimental evidence of physiological relevance of purine sensing in a second-messenger signaling system that modulates c-di-GMP levels.

GenoFig: a user-friendly application for the visualization and comparison of genomic regions

MOTIVATION

Understanding the molecular evolutionary history of organisms usually requires visual comparison of genomic regions from related species or strains. Although several applications already exist to achieve this task, they are either too old, too limited, or too complex for most user's needs.

RESULTS

GenoFig is a graphical application for the visualization of prokaryotic genomic regions, intended to be as easy to use as possible and flexible enough to adapt to a variety of needs. GenoFig allows the personalized representation of annotations extracted from GenBank files in a consistent way across sequences, using regular expressions. It also provides several unique options to optimize the display of homologous regions between sequences, as well as other more classical features such as sequence GC percent or GC-skew representations. In summary, GenoFig is a simple, free, and highly configurable tool to explore the evolution of specific genomic regions in prokaryotes and to produce publication-ready figures.

AVAILABILITY AND IMPLEMENTATION

Genofig is fully available at https://forgemia.inra.fr/public-pgba/genofig under a GPL 3.0 license.

AnnoView enables large-scale analysis, comparison, and visualization of microbial gene neighborhoods

The analysis and comparison of gene neighborhoods is a powerful approach for exploring microbial genome structure, function, and evolution. Although numerous tools exist for genome visualization and comparison, genome exploration across large genomic databases or user-generated datasets remains a challenge. Here, we introduce AnnoView, a web server designed for interactive exploration of gene neighborhoods across the bacterial and archaeal tree of life. Our server offers users the ability to identify, compare, and visualize gene neighborhoods of interest from 30 238 bacterial genomes and 1672 archaeal genomes, through integration with the comprehensive Genome Taxonomy Database and AnnoTree databases. Identified gene neighborhoods can be visualized using pre-computed functional annotations from different sources such as KEGG, Pfam and TIGRFAM, or clustered based on similarity. Alternatively, users can upload and explore their own custom genomic datasets in GBK, GFF or CSV format, or use AnnoView as a genome browser for relatively small genomes (e.g. viruses and plasmids). Ultimately, we anticipate that AnnoView will catalyze biological discovery by enabling user-friendly search, comparison, and visualization of genomic data. AnnoView is available at http://annoview.uwaterloo.ca.

FlhE functions as a chaperone to prevent formation of periplasmic flagella in Gram-negative bacteria

The bacterial flagellum is an organelle utilized by many Gram-negative bacteria to facilitate motility. The flagellum is composed of a several µm long, extracellular filament that is connected to a cytoplasmic rotor-stator complex via a periplasmic rod. Composed of ∼20 structural proteins, ranging from a few subunits to several thousand building blocks, the flagellum is a paradigm of a complex macromolecular structure that utilizes a highly regulated assembly process. This process is governed by multiple checkpoints that ensure an ordered gene expression pattern coupled to the assembly of the various flagellar building blocks in order to produce a functional flagellum. Using epifluorescence, super-resolution STED and transmission electron microscopy, we discovered that in Salmonella , the absence of one periplasmic protein, FlhE, prevents proper flagellar morphogenesis and results in the formation of periplasmic flagella. The periplasmic flagella disrupt cell wall synthesis, leading to a loss of the standard cell morphology resulting in cell lysis. We propose a model where FlhE functions as a periplasmic chaperone to control assembly of the periplasmic rod to prevent formation of periplasmic flagella. Our results highlight that bacteria evolved sophisticated regulatory mechanisms to control proper flagellar assembly and minor deviations from this highly regulated process can cause dramatic physiological consequences.

NADP+-dependent isocitrate dehydrogenase as a novel target for altering carbon flux to lipid accumulation and enhancing antioxidant capacity in Tetradesmus obliquus

Regulatory complexities in lipogenesis hinder the harmonization of metabolic carbon precursors towards lipid synthesis. Exploring regulatory complexities in lipogenesis, this study identifies NADP+-dependent isocitrate dehydrogenase (IDH) in Tetradesmus obliquus as a key factor. Overexpression IDH in strains ToIDH-1 and ToIDH-2 resulted in a 1.69 and 1.64-fold increase in neutral lipids, respectively, compared to the wild type, with lipid yield reaching 234.56 and 227.17 mg/L. Notably, despite slower growth, the cellular biomass augmented to 790.67 mg/L. Metabolite analysis indicated a shift in carbon precursors from protein to lipid and carbohydrate synthesis. Morphological observations revealed increases in the volume and number of lipid droplets, alongside a change in the fatty acid profile favoring monounsaturated and saturated fatty acids. Furthermore, IDH overexpression enhanced NADPH production and antioxidant activity, thereby further boosting lipid accumulation when combined with salt stress. This study suggests a pathway for improved lipogenesis and algal growth via metabolic engineering.

Diverse domain architectures of CheA histidine kinase, a central component of bacterial and archaeal chemosensory systems

IMPORTANCE

We found that in contrast to the best-studied model organisms, such as Escherichia coli and Bacillus subtilis, most bacterial and archaeal species have a CheA protein with a different domain composition. We report variations in CheA architecture, such as domain duplication and acquisition as well as class-specific domain composition. Our results will be of interest to those working on signal transduction in bacteria and archaea and lay the foundation for experimental studies.

Identification of nitric oxide mediated defense signaling and its microRNA mediated regulation during Phytophthora capsici infection in black pepper

Nitric oxide plays a significant role in the defense signaling during pathogen interaction in plants. Quick wilt disease is a devastating disease of black pepper, and leads to sudden mortality of pepper vines in plantations. In this study, the role of nitric oxide was studied during Phytophthora capsici infection in black pepper variety Panniyur-1. Nitric oxide was detected from the different histological sections of P. capsici infected leaves. Furthermore, the genome-wide transcriptome analysis characterized typical domain architect and structural features of nitrate reductase (NR) and nitric oxide associated 1 (NOA1) gene that are involved in nitric oxide biosynthesis in black pepper. Despite the upregulation of nitrate reductase (Pn1_NR), a reduced expression of Pn1_NOA1 was detected in the P. capsici infected black pepper leaf. Subsequent sRNAome-assisted in silico analysis revealed possible microRNA mediated regulation of Pn1_NOA mRNAs. Furthermore, sRNA/miRNA mediated cleavage on Pn1_NOA1 mRNA was validated through modified 5' RLM RACE experiments. Several hormone-responsive cis-regulatory elements involved in stress response was detected from the promoter regions of Pn_NOA1, Pn_NR1 and Pn_NR2 genes. Our results revealed the role of nitric oxide during stress response of P. capsici infection in black pepper, and key genes involved in nitric oxide biosynthesis and their post-transcriptional regulatory mechanisms.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-024-01414-z.

Ubiquitous purine sensor modulates diverse signal transduction pathways in bacteria

Purines and their derivatives are key molecules for controlling intracellular energy homeostasis and nucleotide synthesis. In eukaryotes, including humans, purines also act as signaling molecules that mediate extracellular communication and control key cellular processes, such as proliferation, migration, differentiation, and apoptosis. However, the signaling role of purines in bacteria is largely unknown. Here, by combining structural and sequence information, we define a purine-binding motif, which is present in sensor domains of thousands of bacterial receptors that modulate motility, gene expression, metabolism and second messenger turnover. The screening of compound libraries and microcalorimetric titrations of selected sensor domains validated their ability to specifically bind purine derivatives. The physiological relevance of purine sensing was demonstrated in a second messenger signaling system that modulates c-di-GMP levels.

Amine-recognizing domain in diverse receptors from bacteria and archaea evolved from the universal amino acid sensor

Bacteria possess various receptors that sense different signals and transmit information to enable an optimal adaptation to the environment. A major limitation in microbiology is the lack of information on the signal molecules that activate receptors. Signals recognized by sensor domains are poorly reflected in overall sequence identity, and therefore, the identification of signals from the amino acid sequence of the sensor alone presents a challenge. Biogenic amines are of great physiological importance for microorganisms and humans. They serve as substrates for aerobic and anaerobic growth and play a role of neurotransmitters and osmoprotectants. Here, we report the identification of a sequence motif that is specific for amine-sensing sensor domains that belong to the Cache superfamily of the most abundant extracellular sensors in prokaryotes. We identified approximately 13,000 sensor histidine kinases, chemoreceptors, receptors involved in second messenger homeostasis and Ser/Thr phosphatases from 8,000 bacterial and archaeal species that contain the amine-recognizing motif. The screening of compound libraries and microcalorimetric titrations of selected sensor domains confirmed their ability to specifically bind biogenic amines. Mutants in the amine-binding motif or domains that contain a single mismatch in the binding motif had either no or a largely reduced affinity for amines. We demonstrate that the amine-recognizing domain originated from the universal amino acid-sensing Cache domain, thus providing insight into receptor evolution. Our approach enables precise "wet"-lab experiments to define the function of regulatory systems and therefore holds a strong promise to enable the identification of signals stimulating numerous receptors.

Diverse domain architectures of CheA histidine kinase, a central component of bacterial and archaeal chemosensory systems

Chemosensory systems in bacteria and archaea are complex, multi-protein pathways that enable rapid cellular responses to environmental changes. The CheA histidine kinase is a central component of chemosensory systems. In contrast to other histidine kinases, it lacks a sensor (input) domain and utilizes dedicated chemoreceptors for sensing. CheA is a multi-domain protein; in model organisms as diverse as Escherichia coli and Bacillus subtilis, it contains five single-copy domains. Deviations from this canonical domain architecture have been reported, however, a broad genome-wide analysis of CheA diversity is lacking. Here, we present results of a genomic survey of CheA domain composition carried out using an unbiased set of thousands of CheA sequences from bacteria and archaea. We found that four out of five canonical CheA domains comprise a minimal functional unit (core domains), as they are present in all surveyed CheA homologs. The most common deviations from a classical five-domain CheA architecture are the lack of a P2/CheY-binding domain, which is missing from more than a half of CheA homologs and the acquisition of a response regulator receiver (CheY-like) domain, which is present in ~35% of CheA homologs. We also document other deviations from classical CheA architecture, including bipartite CheA proteins, domain duplications and fusions, and reveal that phylogenetically defined CheA classes have pre-dominant domain architectures. This study lays a foundation for a better classification of CheA homologs and identifies targets for experimental investigations.

Sensing preferences for prokaryotic solute binding protein families

Solute binding proteins (SBPs) are of central physiological relevance for prokaryotes. These proteins present substrates to transporters, but they also stimulate different signal transduction receptors. SBPs form a superfamily of at least 33 protein Pfam families. To assess possible links between SBP sequence and the ligand recognized, we have inspected manually all SBP three-dimensional structures deposited in the protein data bank and retrieved 748 prokaryotic structures that have been solved in complex with bound ligand. These structures were classified into 26 SBP Pfam families. The analysis of the ligands recognized revealed that most families possess a preference for a compound class. There were three families each that bind preferentially saccharides and amino acids. In addition, we identified families that bind preferentially purines, quaternary amines, iron and iron-chelating compounds, oxoanions, bivalent metal ions or phosphates. Phylogenetic analyses suggest convergent evolutionary events that lead to families that bind the same ligand. The functional link between chemotaxis and compound uptake is reflected in similarities in the ligands recognized by SBPs and chemoreceptors. Associating Pfam families with ligand profiles will be of help to design experimental strategies aimed at the identification of ligands for uncharacterized SBPs.

Origin and functional diversification of PAS domain, a ubiquitous intracellular sensor

Signal perception is a key function in regulating biological activities and adapting to changing environments. Per-Arnt-Sim (PAS) domains are ubiquitous sensors found in diverse receptors in bacteria, archaea, and eukaryotes, but their origins, distribution across the tree of life, and extent of their functional diversity are not fully characterized. Here, we show that using sequence conservation and structural information, it is possible to propose specific and potential functions for a large portion of nearly 3 million PAS domains. Our analysis suggests that PAS domains originated in bacteria and were horizontally transferred to archaea and eukaryotes. We reveal that gas sensing via a heme cofactor evolved independently in several lineages, whereas redox and light sensing via flavin adenine dinucleotide and flavin mononucleotide cofactors have the same origin. The close relatedness of human PAS domains to those in bacteria provides an opportunity for drug design by exploring potential natural ligands and cofactors for bacterial homologs.

Sm-like protein Rof inhibits transcription termination factor ρ by binding site obstruction and conformational insulation

Transcription termination factor ρ is a hexameric, RNA-dependent NTPase that can adopt active closed-ring and inactive open-ring conformations. The Sm-like protein Rof, a homolog of the RNA chaperone Hfq, inhibits ρ-dependent termination in vivo but recapitulation of this activity in vitro has proven difficult and the precise mode of Rof action is presently unknown. Our electron microscopic structures of ρ-Rof and ρ-RNA complexes show that Rof undergoes pronounced conformational changes to bind ρ at the protomer interfaces, undercutting ρ conformational dynamics associated with ring closure and occluding extended primary RNA-binding sites that are also part of interfaces between ρ and RNA polymerase. Consistently, Rof impedes ρ ring closure, ρ-RNA interactions, and ρ association with transcription elongation complexes. Structure-guided mutagenesis coupled with functional assays confirmed that the observed ρ-Rof interface is required for Rof-mediated inhibition of cell growth and ρ-termination in vitro. Bioinformatic analyses revealed that Rof is restricted to Pseudomonadota and that the ρ-Rof interface is conserved. Genomic contexts of rof differ between Enterobacteriaceae and Vibrionaceae, suggesting distinct modes of Rof regulation. We hypothesize that Rof and other cellular anti-terminators silence ρ under diverse, but yet to be identified, stress conditions when unrestrained transcription termination by ρ would be lethal.

Transmembrane Transcription Regulators Are Widespread in Bacteria and Archaea

To adapt and proliferate, bacteria must sense and respond to the ever-changing environment. Transmembrane transcription regulators (TTRs) are a family of one-component transcription regulators that respond to extracellular information and influence gene expression from the cytoplasmic membrane. How TTRs function to modulate expression of their target genes while localized to the cytoplasmic membrane remains poorly understood. In part, this is due to a lack of knowledge regarding the prevalence of TTRs among prokaryotes. Here, we show that TTRs are highly diverse and prevalent throughout bacteria and archaea. Our work demonstrates that TTRs are more common than previously appreciated and are enriched within specific bacterial and archaeal phyla and that many TTRs have unique transmembrane region properties that can facilitate association with detergent-resistant membranes. IMPORTANCE One-component signal transduction systems are the major class of signal transduction systems among bacteria and are commonly cytoplasmic. TTRs are a group of unique one-component signal transduction systems that influence transcription from the cytoplasmic membrane. TTRs have been implicated in a wide array of biological pathways critical for both pathogens and human commensal organisms but were considered to be rare. Here, we demonstrate that TTRs are in fact highly diverse and broadly distributed in bacteria and archaea. Our findings suggest that transcription factors can access the chromosome and influence transcription from the membrane in both archaea and bacteria. This study challenges thus the commonly held notion that signal transduction systems require a cytoplasmic transcription factor and highlights the importance of the cytoplasmic membrane in directly influencing signal transduction.

Amine recognizing domain in diverse receptors from bacteria and archaea evolved from the universal amino acid sensor

Bacteria contain many different receptor families that sense different signals permitting an optimal adaptation to the environment. A major limitation in microbiology is the lack of information on the signal molecules that activate receptors. Due to a significant sequence divergence, the signal recognized by sensor domains is only poorly reflected in overall sequence identity. Biogenic amines are of central physiological relevance for microorganisms and serve for example as substrates for aerobic and anaerobic growth, neurotransmitters or osmoprotectants. Based on protein structural information and sequence analysis, we report here the identification of a sequence motif that is specific for amine-sensing dCache sensor domains (dCache_1AM). These domains were identified in more than 13,000 proteins from 8,000 bacterial and archaeal species. dCache_1AM containing receptors were identified in all major receptor families including sensor kinases, chemoreceptors, receptors involved in second messenger homeostasis and Ser/Thr phosphatases. The screening of compound libraries and microcalorimetric titrations of selected dCache_1AM domains confirmed their capacity to specifically bind amines. Mutants in the amine binding motif or domains that contain a single mismatch in the binding motif, had either no or a largely reduced affinity for amines, illustrating the specificity of this motif. We demonstrate that the dCache_1AM domain has evolved from the universal amino acid sensing domain, providing novel insight into receptor evolution. Our approach enables precise "wet"-lab experiments to define the function of regulatory systems and thus holds a strong promise to address an important bottleneck in microbiology: the identification of signals that stimulate numerous receptors.

Five copper homeostasis gene clusters encode the Cu-efflux resistome of the highly copper-tolerant Methylorubrum extorquens AM1

Background

In the last decade, the use of copper has reemerged as a potential strategy to limit healthcare-associated infections and to control the spread of multidrug-resistant pathogens. Numerous environmental studies have proposed that most opportunistic pathogens have acquired antimicrobial resistance in their nonclinical primary habitat. Thus, it can be presumed that copper-resistant bacteria inhabiting a primary commensal niche might potentially colonize clinical environments and negatively affect the bactericidal efficacy of Cu-based treatments. The use of copper in agricultural fields is one of the most important sources of Cu pollution that may exert selection pressure for the increase of copper resistance in soil and plant-associated bacteria. To assess the emergence of copper-resistant bacteria in natural habitats, we surveyed a laboratory collection of bacterial strains belonging to the order Rhizobiales. This study proposes that Methylorubrum extorquens AM1 is an environmental isolate well adapted to thrive in copper-rich environments that could act as a reservoir of copper resistance genes.

Methods

The minimal inhibitory concentrations (MICs) of CuCl₂ were used to estimate the copper tolerance of eight plant-associated facultative diazotrophs (PAFD) and five pink-pigmented facultative methylotrophs (PPFM) belonging to the order Rhizobiales presumed to come from nonclinical and nonmetal-polluted natural habitats based on their reported source of isolation. Their sequenced genomes were used to infer the occurrence and diversity of Cu-ATPases and the copper efflux resistome of Mr. extorquens AM1.

Results

These bacteria exhibited minimal inhibitory concentrations (MICs) of CuCl₂ ranging between 0.020 and 1.9 mM. The presence of multiple and quite divergent Cu-ATPases per genome was a prevalent characteristic. The highest copper tolerance exhibited by Mr. extorquens AM1 (highest MIC of 1.9 mM) was similar to that found in the multimetal-resistant model bacterium Cupriavidus metallidurans CH34 and in clinical isolates of Acinetobacter baumannii. The genome-predicted copper efflux resistome of Mr. extorquens AM1 consists of five large (6.7 to 25.7 kb) Cu homeostasis gene clusters, three clusters share genes encoding Cu-ATPases, CusAB transporters, numerous CopZ chaperones, and enzymes involved in DNA transfer and persistence. The high copper tolerance and the presence of a complex Cu efflux resistome suggest the presence of relatively high copper tolerance in environmental isolates of Mr. extorquens.

Photoactive Yellow Protein Represents a Distinct, Evolutionarily Novel Family of PAS Domains

Photoactive yellow protein (PYP) is a model photoreceptor. It binds a p-coumaric acid as a chromophore, thus enabling blue light sensing. The first discovered single-domain PYP from Halorhodospira halophila has been studied thoroughly in terms of its structural dynamics and photochemical properties. However, the evolutionary origins and biological role of PYP homologs are not well understood. Here, we show that PYP is an evolutionarily novel domain family of the ubiquitous PAS (Per-Arnt-Sim) superfamily. It likely originated from the phylum Myxococcota and was then horizontally transferred to representatives of a few other bacterial phyla. We show that PYP is associated with signal transduction either by domain fusion or by genome context. Key cellular functions modulated by PYP-initiated signal transduction pathways likely involve gene expression, motility, and biofilm formation. We identified three clades of the PYP family, one of which is poorly understood and potentially has novel functional properties. The Tyr42, Glu46, and Cys69 residues that are involved in p-coumaric acid binding in the model PYP from H. halophila are well conserved in the PYP family. However, we also identified cases where substitutions in these residues might have led to neofunctionalization, such as the proposed transition from light to redox sensing. Overall, this study provides definition, a newly built hidden Markov model, and the current genomic landscape of the PYP family and sets the stage for the future exploration of its signaling mechanisms and functional diversity. IMPORTANCE Photoactive yellow protein is a model bacterial photoreceptor. For many years, it was considered a prototypical model of the ubiquitous PAS domain superfamily. Here, we show that, in fact, the PYP family is evolutionarily novel, restricted to a few bacterial phyla and distinct from other PAS domains. We also reveal the diversity of PYP-containing signal transduction proteins and their potential mechanisms.

The Repertoire of Solute-Binding Proteins of Model Bacteria Reveals Large Differences in Number, Type, and Ligand Range

Solute-binding proteins (SBPs) are of central physiological relevance for bacteria. They are located in the extracytosolic space, where they present substrates to transporters but also stimulate different types of transmembrane receptors coordinating compound uptake with signal transduction. SBPs are a superfamily composed of proteins recognized by 45 Pfam profiles. The definition of SBP profiles for bacteria is hampered by the fact that these Pfam profiles recognize sensor domains for different types of signaling proteins or cytosolic proteins with alternative functions. We report here the retrieval of the SBPs from 49 bacterial model strains with different lifestyles and phylogenetic distributions. Proteins were manually curated, and the ligands recognized were predicted bioinformatically. There were very large differences in the number and type of SBPs between strains, ranging from 7 SBPs in Helicobacter pylori 26695 to 189 SBPs in Sinorhizobium meliloti 1021. SBPs were found to represent 0.22 to 5.13% of the total protein-encoding genes. The abundance of SBPs was largely determined by strain phylogeny, and no obvious link with the bacterial lifestyle was noted. Most abundant (36%) were SBPs predicted to recognize amino acids or peptides, followed by those expected to bind different sugars (18%). To the best of our knowledge, this is the first comparative study of bacterial SBP repertoires. Given the importance of SBPs in nutrient uptake and signaling, this study enhances the knowledge of model bacteria and will permit the definition of SBP profiles of other strains. IMPORTANCE SBPs are essential components for many transporters, but multiple pieces of more recent evidence indicate that the SBP-mediated stimulation of different transmembrane receptors is a general and widespread signal transduction mechanism in bacteria. The double function of SBPs in coordinating transport with signal transduction remains to a large degree unexplored and represents a major research need. The definition of the SBP repertoire of the 49 bacterial model strains examined here, along with information on their cognate ligand profiles forms the basis to close this gap in knowledge. Furthermore, this study provides information on the forces that have driven the evolution of transporters with different ligand specificities in bacteria that differ in phylogenetics and lifestyle. This article is also a first step in setting up automatic algorithms that permit the large-scale identification of the SBP repertoire in proteomes.

Genome-guided prediction of acid resistance mechanisms in acidophilic methanotrophs of phylogenetically deep-rooted Verrucomicrobia isolated from geothermal environments

Verrucomicrobia are a group of microorganisms that have been proposed to be deeply rooted in the Tree of Life. Some are methanotrophs that oxidize the potent greenhouse gas methane and are thus important in decreasing atmospheric concentrations of the gas, potentially ameliorating climate change. They are widespread in various environments including soil and fresh or marine waters. Recently, a clade of extremely acidophilic Verrucomicrobia, flourishing at pH < 3, were described from high-temperature geothermal ecosystems. This novel group could be of interest for studies about the emergence of life on Earth and to astrobiologists as homologs for possible extraterrestrial life. In this paper, we describe predicted mechanisms for survival of this clade at low pH and suggest its possible evolutionary trajectory from an inferred neutrophilic ancestor. Extreme acidophiles are defined as organisms that thrive in extremely low pH environments (≤ pH 3). Many are polyextremophiles facing high temperatures and high salt as well as low pH. They are important to study for both providing fundamental insights into biological mechanisms of survival and evolution in such extreme environments and for understanding their roles in biotechnological applications such as industrial mineral recovery (bioleaching) and mitigation of acid mine drainage. They are also, potentially, a rich source of novel genes and pathways for the genetic engineering of microbial strains. Acidophiles of the Verrucomicrobia phylum are unique as they are the only known aerobic methanotrophs that can grow optimally under acidic (pH 2–3) and moderately thermophilic conditions (50–60°C). Three moderately thermophilic genera, namely Methylacidiphilum, Methylacidimicrobium, and Ca. Methylacidithermus, have been described in geothermal environments. Most of the investigations of these organisms have focused on their methane oxidizing capabilities (methanotrophy) and use of lanthanides as a protein cofactor, with no extensive study that sheds light on the mechanisms that they use to flourish at extremely low pH. In this paper, we extend the phylogenetic description of this group of acidophiles using whole genome information and we identify several mechanisms, potentially involved in acid resistance, including “first line of defense” mechanisms that impede the entry of protons into the cell. These include the presence of membrane-associated hopanoids, multiple copies of the outer membrane protein (Slp), and inner membrane potassium channels (kup, kdp) that generate a reversed membrane potential repelling the intrusion of protons. Acidophilic Verrucomicrobia also display a wide array of proteins potentially involved in the “second line of defense” where protons that evaded the first line of defense and entered the cell are expelled or neutralized, such as the glutamate decarboxylation (gadAB) and phosphate-uptake systems. An exclusive N-type ATPase F₀-F₁ was identified only in acidophiles of Verrucomicrobia and is predicted to be a specific adaptation in these organisms. Phylogenetic analyses suggest that many predicted mechanisms are evolutionarily conserved and most likely entered the acidophilic lineage of Verrucomicrobia by vertical descent from a common ancestor. However, it is likely that some defense mechanisms such as gadA and kup entered the acidophilic Verrucomicrobia lineage by horizontal gene transfer.

Methylation-Independent Chemotaxis Systems Are the Norm for Gastric-Colonizing Helicobacter Species

Many bacteria and archaea rely on chemotaxis signal transduction systems for optimal fitness. These complex, multiprotein signaling systems have core components found in all chemotactic microbes, as well as variable proteins found in only some species. We do not yet understand why these variations exist or whether there are specific niches that favor particular chemotaxis signaling organization. One variation is in the presence/absence of the chemotaxis methylation adaptation enzymes CheB and CheR. Genes for CheB and CheR are missing in the gastric pathogen Helicobacter pylori but present in related Helicobacter that colonize the liver or intestine. In this work, we asked whether there was a general pattern of CheB/CheR across multiple Helicobacter species. Helicobacter spp. all possess chemotactic behavior, based on the presence of genes for core signaling proteins CheA, CheW, and chemoreceptors. Genes for the CheB and CheR proteins, in contrast, were variably present. Niche mapping supported the idea that these genes were present in enterohepatic Helicobacter species and absent in gastric ones. We then analyzed whether there were differences between gastric and enterohepatic species in the CheB/CheR chemoreceptor target methylation sites. Indeed, these sites were less conserved in gastric species that lack CheB/CheR. Lastly, we determined that cheB and cheR could serve as markers to indicate whether an unknown Helicobacter species was of enterohepatic or gastric origin. Overall, these findings suggest the interesting idea that methylation-based adaptation is not required in specific environments, particularly the stomach. IMPORTANCE Chemotaxis signal transduction systems are common in the archaeal and bacterial world, but not all systems contain the same components. The rationale for this system variation remains unknown. In this report, comparative genomics analysis showed that the presence/absence of CheR and CheB is one main variation within the Helicobacter genus, and it is strongly associated with the niche of Helicobacter species: gastric Helicobacter species, which infect animal stomachs, have lost their CheB and CheR, while enterohepatic Helicobacter species, which infect the liver and intestine, retain them. This study not only provides an example that a chemotaxis system variant is associated with particular niches but also proposes that CheB and CheR are new markers distinguishing gastric from enterohepatic Helicobacter species.

The Natural Product Domain Seeker version 2 (NaPDoS2) webtool relates ketosynthase phylogeny to biosynthetic function

The Natural Product Domain Seeker (NaPDoS) webtool detects and classifies ketosynthase (KS) and condensation domains from genomic, metagenomic, and amplicon sequence data. Unlike other tools, a phylogeny-based classification scheme is used to make broader predictions about the polyketide synthase (PKS) and nonribosomal peptide synthetase (NRPS) genes in which these domains are found. NaPDoS is particularly useful for the analysis of incomplete biosynthetic genes or gene clusters, as are often observed in poorly assembled genomes and metagenomes, or when loci are not clustered, as in eukaryotic genomes. To help support the growing interest in sequence-based analyses of natural product biosynthetic diversity, here we introduce version 2 of the webtool, NaPDoS2, available at http://napdos.ucsd.edu/napdos2. This update includes the addition of 1417 KS sequences, representing a major expansion of the taxonomic and functional diversity represented in the webtool database. The phylogeny-based KS classification scheme now recognizes 41 class and subclass assignments, including new type II PKS subclasses. Workflow modifications accelerate run times, allowing larger datasets to be analyzed. In addition, default parameters were established using statistical validation tests to maximize KS detection and classification accuracy while minimizing false positives. We further demonstrate the applications of NaPDoS2 to assess PKS biosynthetic potential using genomic, metagenomic, and PCR amplicon datasets. These examples illustrate how NaPDoS2 can be used to predict biosynthetic potential and detect genes involved in the biosynthesis of specific structure classes or new biosynthetic mechanisms.

GeCoViz: genomic context visualisation of prokaryotic genes from a functional and evolutionary perspective

Synteny conservation analysis is a well-established methodology to investigate the potential functional role of unknown prokaryotic genes. However, bioinformatic tools to reconstruct and visualise genomic contexts usually depend on slow computations, are restricted to narrow taxonomic ranges, and/or do not allow for the functional and interactive exploration of neighbouring genes across different species. Here, we present GeCoViz, an online resource built upon 12 221 reference prokaryotic genomes that provides fast and interactive visualisation of custom genomic regions anchored by any target gene, which can be sought by either name, orthologous group (KEGGs, eggNOGs), protein domain (PFAM) or sequence. To facilitate functional and evolutionary interpretation, GeCoViz allows to customise the taxonomic scope of each analysis and provides comprehensive annotations of the neighbouring genes. Interactive visualisation options include, among others, the scaled representations of gene lengths and genomic distances, and on the fly calculation of synteny conservation of neighbouring genes, which can be highlighted based on custom thresholds. The resulting plots can be downloaded as high-quality images for publishing purposes. Overall, GeCoViz offers an easy-to-use, comprehensive, fast and interactive web-based tool for investigating the genomic context of prokaryotic genes, and is freely available at https://gecoviz.cgmlab.org.

Characterization of Opposing Responses to Phenol by Bacillus subtilis Chemoreceptors

Bacillus subtilis employs 10 chemoreceptors to move in response to chemicals in its environment. While the sensing mechanisms have been determined for many attractants, little is known about the sensing mechanisms for repellents. In this work, we investigated phenol chemotaxis in B. subtilis. Phenol is an attractant at low, micromolar concentrations and a repellent at high, millimolar concentrations. McpA was found to be the principal chemoreceptor governing the repellent response to phenol and other related aromatic compounds. In addition, the chemoreceptors McpC and HemAT were found to govern the attractant response to phenol and related compounds. Using chemoreceptor chimeras, McpA was found to sense phenol using its signaling domain rather than its sensing domain. These observations were substantiated in vitro, where direct binding of phenol to the signaling domain of McpA was observed using saturation transfer difference nuclear magnetic resonance. These results further advance our understanding of B. subtilis chemotaxis and further demonstrate that the signaling domain of B. subtilis chemoreceptors can directly sense chemoeffectors. IMPORTANCE Bacterial chemotaxis is commonly thought to employ a sensing mechanism involving the extracellular sensing domain of chemoreceptors. Some ligands, however, appear to be sensed by the signaling domain. Phenolic compounds, commonly found in soil and root exudates, provide environmental cues for soil microbes like Bacillus subtilis. We show that phenol is sensed as both an attractant and a repellent. While the mechanism for sensing phenol as an attractant is still unknown, we found that phenol is sensed as a repellent by the signaling domain of the chemoreceptor McpA. This study furthers our understanding of the unconventional sensing mechanisms employed by the B. subtilis chemotaxis pathway.

Multiple profile models extract features from protein sequence data and resolve functional diversity of very different protein families

Functional classification of proteins from sequences alone has become a critical bottleneck in understanding the myriad of protein sequences that accumulate in our databases. The great diversity of homologous sequences hides, in many cases, a variety of functional activities that cannot be anticipated. Their identification appears critical for a fundamental understanding of the evolution of living organisms and for biotechnological applications. ProfileView is a sequence-based computational method, designed to functionally classify sets of homologous sequences. It relies on two main ideas: the use of multiple profile models whose construction explores evolutionary information in available databases, and a novel definition of a representation space in which to analyse sequences with multiple profile models combined together. ProfileView classifies protein families by enriching known functional groups with new sequences and discovering new groups and subgroups. We validate ProfileView on seven classes of widespread proteins involved in the interaction with nucleic acids, amino acids and small molecules, and in a large variety of functions and enzymatic reactions. Profile-View agrees with the large set of functional data collected for these proteins from the literature regarding the organisation into functional subgroups and residues that characterise the functions. In addition, ProfileView resolves undefined functional classifications and extracts the molecular determinants underlying protein functional diversity, showing its potential to select sequences towards accurate experimental design and discovery of novel biological functions. On protein families with complex domain architecture, ProfileView functional classification reconciles domain combinations, unlike phylogenetic reconstruction. ProfileView proves to outperform the functional classification approach PANTHER, the two k-mer based methods CUPP and eCAMI and a neural network approach based on Restricted Boltzmann Machines. It overcomes time complexity limitations of the latter.

Amino acid sensor conserved from bacteria to humans

SignificanceAmino acids are the building blocks of life and important signaling molecules. Despite their common structure, no universal mechanism for amino acid recognition by cellular receptors is currently known. We discovered a simple motif, which binds amino acids in various receptor proteins from all major life-forms. In humans, this motif is found in subunits of calcium channels that are implicated in pain and neurodevelopmental disorders. Our findings suggest that γ-aminobutyric acid-derived drugs bind to the same motif in human proteins that binds natural ligands in bacterial receptors, thus enabling future improvement of important drugs.

Bacterial Sirtuins Overview: An Open Niche to Explore

Sirtuins are deacetylase enzymes widely distributed in all domains of life. Although for decades they have been related only to histones deacetylation in eukaryotic organisms, today they are considered global regulators in both prokaryotes and eukaryotes. Despite the important role of sirtuins in humans, the knowledge about bacterial sirtuins is still limited. Several proteomics studies have shown that bacterial sirtuins deacetylate a large number of lysines in vivo, although the effect that this deacetylation causes in most of them remains unknown. To date, only the regulation of a few bacterial sirtuin substrates has been characterized, being their metabolic roles widely distributed: carbon and nitrogen metabolism, DNA transcription, protein translation, or virulence. One of the most current topics on acetylation and deacetylation focuses on studying stoichiometry using quantitative LC-MS/MS. The results suggest that prokaryotic sirtuins deacetylate at low stoichiometry sites, although more studies are needed to know if it is a common characteristic of bacterial sirtuins and its biological significance. Unlike eukaryotic organisms, bacteria usually have one or few sirtuins, which have been reported to have closer phylogenetic similarity with the human Sirt5 than with any other human sirtuin. In this work, in addition to carrying out an in-depth review of the role of bacterial sirtuins in their physiology, a phylogenetic study has been performed that reveals the evolutionary differences between sirtuins of different bacterial species and even between homologous sirtuins.

The confluence of big data and evolutionary genome mining for the discovery of natural products

This review covers literature between 2003-2021The development and application of genome mining tools has given rise to ever-growing genetic and chemical databases and propelled natural products research into the modern age of Big Data. Likewise, an explosion of evolutionary studies has unveiled genetic patterns of natural products biosynthesis and function that support Darwin's theory of natural selection and other theories of adaptation and diversification. In this review, we aim to highlight how Big Data and evolutionary thinking converge in the study of natural products, and how this has led to an emerging sub-discipline of evolutionary genome mining of natural products. First, we outline general principles to best utilize Big Data in natural products research, addressing key considerations needed to provide evolutionary context. We then highlight successful examples where Big Data and evolutionary analyses have been combined to provide bioinformatic resources and tools for the discovery of novel natural products and their biosynthetic enzymes. Rather than an exhaustive list of evolution-driven discoveries, we highlight examples where Big Data and evolutionary thinking have been embraced for the evolutionary genome mining of natural products. After reviewing the nascent history of this sub-discipline, we discuss the challenges and opportunities of genomic and metabolomic tools with evolutionary foundations and/or implications and provide a future outlook for this emerging and exciting field of natural product research.

Synthaser: a CD-Search enabled Python toolkit for analysing domain architecture of fungal secondary metabolite megasynth(et)ases

BACKGROUND

Fungi are prolific producers of secondary metabolites (SMs), which are bioactive small molecules with important applications in medicine, agriculture and other industries. The backbones of a large proportion of fungal SMs are generated through the action of large, multi-domain megasynth(et)ases such as polyketide synthases (PKSs) and nonribosomal peptide synthetases (NRPSs). The structure of these backbones is determined by the domain architecture of the corresponding megasynth(et)ase, and thus accurate annotation and classification of these architectures is an important step in linking SMs to their biosynthetic origins in the genome.

RESULTS

Here we report synthaser, a Python package leveraging the NCBI's conserved domain search tool for remote prediction and classification of fungal megasynth(et)ase domain architectures. Synthaser is capable of batch sequence analysis, and produces rich textual output and interactive visualisations which allow for quick assessment of the megasynth(et)ase diversity of a fungal genome. Synthaser uses a hierarchical rule-based classification system, which can be extensively customised by the user through a web application ( http://gamcil.github.io/synthaser ). We show that synthaser provides more accurate domain architecture predictions than comparable tools which rely on curated profile hidden Markov model (pHMM)-based approaches; the utilisation of the NCBI conserved domain database also allows for significantly greater flexibility compared to pHMM approaches. In addition, we demonstrate how synthaser can be applied to large scale genome mining pipelines through the construction of an Aspergillus PKS similarity network.

CONCLUSIONS

Synthaser is an easy to use tool that represents a significant upgrade to previous domain architecture analysis tools. It is freely available under a MIT license from PyPI ( https://pypi.org/project/synthaser ) and GitHub ( https://github.com/gamcil/synthaser ).

Identification and characterization of a non-canonical menaquinone-linked formate dehydrogenase

The Molybdenum/Tungsten-bispyranopterin guanine dinucleotides (Mo/W-bisPGD) family of Formate Dehydrogenases (FDHs) plays roles in several metabolic pathways ranging from carbon fixation to energy harvesting owing to their reaction with a wide variety of redox partners. Indeed, this metabolic plasticity results from the diverse structures, cofactor content, and substrates employed by partner subunits interacting with the catalytic hub. Here, we unveiled two non-canonical FDHs in Bacillus subtilis which are organized into two-subunit complexes with unique features, ForCE1 and ForCE2. We show that the ForC catalytic subunit interacts with an unprecedented partner subunit, ForE, and that its amino acid sequence within the active site deviates from the consensus residues typically associated with FDH activity, as a histidine residue is naturally substituted with a glutamine. The ForE essential subunit mediates the utilization of menaquinone as an electron acceptor as shown by the formate:menadione oxidoreductase activity of both enzymes, their copurification with menaquinone, and the distinctive detection of a protein-bound neutral menasemiquinone radical by multifrequency electron paramagnetic resonance (EPR) experiments on the purified enzymes. Moreover, EPR characterization of both FDHs reveals the presence of several [Fe-S] clusters with distinct relaxation properties and a weakly anisotropic Mo(V) EPR signature, consistent with the characteristic Mo/bisPGD cofactor of this enzyme family. Altogether, this work enlarges our knowledge of the FDH family by identifying a non-canonical FDH, which differs in terms of architecture, amino acid conservation around the Mo cofactor, and reactivity.

Diversification of LytM Protein Functions in Polar Elongation and Cell Division of Agrobacterium tumefaciens

LytM-domain containing proteins are LAS peptidases (lysostaphin-type enzymes, D-Ala-D-Ala metallopeptidases, and sonic hedgehog) and are known to play diverse roles throughout the bacterial cell cycle through direct or indirect hydrolysis of the bacterial cell wall. A subset of the LytM factors are catalytically inactive but regulate the activity of other cell wall hydrolases and are classically described as cell separation factors NlpD and EnvC. Here, we explore the function of four LytM factors in the alphaproteobacterial plant pathogen Agrobacterium tumefaciens. An LmdC ortholog (Atu1832) and a MepM ortholog (Atu4178) are predicted to be catalytically active. While Atu1832 does not have an obvious function in cell growth or division, Atu4178 is essential for polar growth and likely functions as a space-making endopeptidase that cleaves amide bonds in the peptidoglycan cell wall during elongation. The remaining LytM factors are degenerate EnvC and NlpD orthologs. Absence of these proteins results in striking phenotypes indicative of misregulation of cell division and growth pole establishment. The deletion of an amidase, AmiC, closely phenocopies the deletion of envC suggesting that EnvC might regulate AmiC activity. The NlpD ortholog DipM is unprecedently essential for viability and depletion results in the misregulation of early stages of cell division, contrasting with the canonical view of DipM as a cell separation factor. Finally, we make the surprising observation that absence of AmiC relieves the toxicity induced by dipM overexpression. Together, these results suggest EnvC and DipM may function as regulatory hubs with multiple partners to promote proper cell division and establishment of polarity.

Diverse ATPase Proteins in Mobilomes Constitute a Large Potential Sink for Prokaryotic Host ATP

Prokaryote mobilome genomes rely on host machineries for survival and replication. Given that mobile genetic elements (MGEs) derive their energy from host cells, we investigated the diversity of ATP-utilizing proteins in MGE genomes to determine whether they might be associated with proteins that could suppress related host proteins that consume energy. A comprehensive search of 353 huge phage genomes revealed that up to 9% of the proteins have ATPase domains. For example, ATPase proteins constitute ∼3% of the genomes of Lak phages with ∼550 kbp genomes that occur in the microbiomes of humans and other animals. Statistical analysis shows the number of ATPase proteins increases linearly with genome length, consistent with a large sink for host ATP during replication of megaphages. Using metagenomic data from diverse environments, we found 505 mobilome proteins with ATPase domains fused to diverse functional domains. Among these composite ATPase proteins, 61.6% have known functional domains that could contribute to host energy diversion during the mobilome infection cycle. As many have domains that are known to interact with nucleic acids and proteins, we infer that numerous ATPase proteins are used during replication and for protection from host immune systems. We found a set of uncharacterized ATPase proteins with nuclease and protease activities, displaying unique domain architectures that are energy intensive based on the presence of multiple ATPase domains. In many cases, these composite ATPase proteins genomically co-localize with small proteins in genomic contexts that are reminiscent of toxin-antitoxin systems and phage helicase-antibacterial helicase systems. Small proteins that function as inhibitors may be a common strategy for control of cellular processes, thus could inspire future biochemical experiments for the development of new nucleic acid and protein manipulation tools, with diverse biotechnological applications.

The Epipeptide Biosynthesis Locus epeXEPAB Is Widely Distributed in Firmicutes and Triggers Intrinsic Cell Envelope Stress

The epeXEPAB (formerly yydFGHIJ) locus of Bacillus subtilis encodes a minimalistic biosynthetic pathway for a linear antimicrobial epipeptide, EpeX, which is ribosomally produced and post-translationally processed by the action of the radical-SAM epimerase, EpeE, and a membrane-anchored signal 2 peptide peptidase, EpeP. The ABC transporter EpeAB provides intrinsic immunity against self-produced EpeX, without conferring resistance against extrinsically added EpeX. EpeX specifically targets, and severely perturbs the integrity of the cytoplasmic membrane, which leads to the induction of the Lia-dependent envelope stress response. Here, we provide new insights into the distribution, expression, and regulation of the minimalistic epeXEPAB locus of B. subtilis, as well as the biosynthesis and biological efficiency of the produced epipeptide EpeX*. A comprehensive comparative genomics study demonstrates that the epe-locus is restricted to but widely distributed within the phylum Firmicutes. The gene products of epeXEP are necessary and sufficient for the production of the mature antimicrobial peptide EpeX*. In B. subtilis, the epeXEPAB locus is transcribed from three different promoters, one upstream of epeX (PepeX) and two within epeP (PepeA1 and PepeA2). While the latter two are mostly constitutive, PepeX shows a growth phase-dependent induction at the onset of stationary phase. We demonstrate that this regulation is the result of the antagonistic action of two global regulators: The transition state regulator AbrB keeps the epe locus shut off during exponential growth by direct binding. This tight repression is relieved by the master regulator of sporulation, Spo0A, which counteracts the AbrB-dependent repression of epeXEPAB expression during the transition to stationary phase. The net result of these three -promoters is an expression pattern that ensures EpeAB-dependent autoimmunity prior to EpeX* production. In the absence of EpeAB, the general envelope stress response proteins LiaIH can compensate for the loss of specific autoimmunity by providing sufficient protection against the membrane-perturbating action of EpeX*. Hence, the transcriptional regulation of epe expression and the resulting intrinsic induction of the two corresponding resistance functions, encoded by epeAB and liaIH, are well balanced to provide a need-based immunity against mature EpeX*.

FlaGs and webFlaGs: discovering novel biology through the analysis of gene neighbourhood conservation

SUMMARY

Analysis of conservation of gene neighbourhoods over different evolutionary levels is important for understanding operon and gene cluster evolution, and predicting functional associations. Our tool FlaGs (standing for Flanking Genes) takes a list of NCBI protein accessions as input, clusters neighbourhood-encoded proteins into homologous groups using sensitive sequence searching, and outputs a graphical visualization of the gene neighbourhood and its conservation, along with a phylogenetic tree annotated with flanking gene conservation. FlaGs has demonstrated utility for molecular evolutionary analysis, having uncovered a new toxin-antitoxin system in prokaryotes and bacteriophages. The web tool version of FlaGs (webFlaGs) can optionally include a BLASTP search against a reduced RefSeq database to generate an input accession list and analyse neighbourhood conservation within the same run.

AVAILABILITY AND IMPLEMENTATION

FlaGs can be downloaded from https://github.com/GCA-VH-lab/FlaGs or run online at http://www.webflags.se/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

GCsnap: Interactive Snapshots for the Comparison of Protein-Coding Genomic Contexts

The biological function and evolutionary history of protein-coding genes are not only written in their nucleotide sequences. The comparison of genomic contexts throughout different lineages may highlight genomic mechanisms in the generation of new protein families, while the conservation of gene clusters may unravel, for instance, metabolic pathways. Various tools and databases exist that allow for the analysis and comparison of genomic contexts, but each has its own limitations. Online databases allow for quick comparisons, but only for those genomes for which data were pre-calculated. More advanced tools may allow for the comparison of any genome, but are often limited to a given phylogenetic kingdom or provide only a snapshot of the genomic contexts without further information about the genes involved. Here, we introduce GCsnap, a flexible Python-based tool that allows for the interactive comparison of the genomic contexts of protein-coding genes from any genome at any taxonomic level, integrating them with functional and structural information. By connecting the output to different protein databases, users can navigate through the different genomic contexts from a simple interactive platform, facilitating the further analysis of the contexts found. GCsnap is not limited to a single input format, can perform batch jobs and accepts protein classification maps. Results are stored in detailed, human and machine-readable files, and customizable, publication-ready figures. GCsnap is freely available from https://github.com/JoanaMPereira/GCsnap and can be set up easily on any computer.

The Campylobacter jejuni chemoreceptor Tlp10 has a bimodal ligand-binding domain and specificity for multiple classes of chemoeffectors

Campylobacter jejuni is a bacterial pathogen that is a common cause of enteritis in humans. We identified a previously uncharacterized type of sensory domain in the periplasmic region of the C. jejuni chemoreceptor Tlp10, termed the DAHL domain, that is predicted to have a bimodular helical architecture. Through two independent ligand-binding sites in this domain, Tlp10 responded to molecular aspartate, isoleucine, fumarate, malate, fucose, and mannose as attractants and to arginine, galactose, and thiamine as repellents. Tlp10 also recognized glycan ligands when present as terminal and intermediate residues of complex structures, such as the fucosylated human ganglioside GM1 and Lewis^a antigen. A tlp10 mutant strain lacking the ligand-binding sites was attenuated in its ability to colonize avian caeca and to adhere to cultured human intestinal cells, indicating the potential involvement of the DAHL domain in host colonization and disease. The Tlp10 intracellular signaling domain interacted with the scaffolding proteins CheV and CheW, which couple chemoreceptors to intracellular signaling machinery, and with the signaling domains of other chemoreceptors, suggesting a key role for Tlp10 in signal transduction and incorporation into sensory arrays. We identified the DAHL domain in other bacterial signal transduction proteins, including the essential virulence induction protein VirA from the plant pathogen Agrobacterium tumefaciens Together, these results suggest a potential link between Tlp10 and C. jejuni virulence.

Origins and Molecular Evolution of the NusG Paralog RfaH

The only universally conserved family of transcription factors comprises housekeeping regulators and their specialized paralogs, represented by well-studied NusG and RfaH. Despite their ubiquity, little information is available on the evolutionary origins, functions, and gene targets of the NusG family members. We built a hidden Markov model profile of RfaH and identified its homologs in sequenced genomes. While NusG is widespread among bacterial phyla and coresides with genes encoding RNA polymerase and ribosome in all except extremely reduced genomes, RfaH is mostly limited to Proteobacteria and lacks common gene neighbors. RfaH activates only a few xenogeneic operons that are otherwise silenced by NusG and Rho. Phylogenetic reconstructions reveal extensive duplications and horizontal transfer of rfaH genes, including those borne by plasmids, and the molecular evolution pathway of RfaH, from "early" exclusion of the Rho terminator and tightened RNA polymerase binding to "late" interactions with the ops DNA element and autoinhibition, which together define the RfaH regulon. Remarkably, NusG is not only ubiquitous in Bacteria but also common in plants, where it likely modulates the transcription of plastid genes.IMPORTANCE In all domains of life, NusG-like proteins make contacts similar to those of RNA polymerase and promote pause-free transcription yet may play different roles, defined by their divergent interactions with nucleic acids and accessory proteins, in the same cell. This duality is illustrated by Escherichia coli NusG and RfaH, which silence and activate xenogenes, respectively. We combined sequence analysis and recent functional and structural insights to envision the evolutionary transformation of NusG, a core regulator that we show is present in all cells using bacterial RNA polymerase, into a virulence factor, RfaH. Our results suggest a stepwise conversion of a NusG duplicate copy into a sequence-specific regulator which excludes NusG from its targets but does not compromise the regulation of housekeeping genes. We find that gene duplication and lateral transfer give rise to a surprising diversity within the only ubiquitous family of transcription factors.

CryoEM structure of the type IVa pilus secretin required for natural competence in Vibrio cholerae

Natural transformation is the process by which bacteria take up genetic material from their environment and integrate it into their genome by homologous recombination. It represents one mode of horizontal gene transfer and contributes to the spread of traits like antibiotic resistance. In Vibrio cholerae, a type IVa pilus (T4aP) is thought to facilitate natural transformation by extending from the cell surface, binding to exogenous DNA, and retracting to thread this DNA through the outer membrane secretin, PilQ. Here, we use a functional tagged allele of VcPilQ purified from native V. cholerae cells to determine the cryoEM structure of the VcPilQ secretin in amphipol to ~2.7 Å. We use bioinformatics to examine the domain architecture and gene neighborhood of T4aP secretins in Proteobacteria in comparison with VcPilQ. This structure highlights differences in the architecture of the T4aP secretin from the type II and type III secretion system secretins. Based on our cryoEM structure, we design a series of mutants to reversibly regulate VcPilQ gate dynamics. These experiments support the idea of VcPilQ as a potential druggable target and provide insight into the channel that DNA likely traverses to promote the spread of antibiotic resistance via horizontal gene transfer by natural transformation.