Transcription Factors Exhibit Differential Conservation in Bacteria with Reduced Genomes

No Transcriptional Compensation for Extreme Gene Dosage Imbalance in Fragmented Bacterial Endosymbionts of Cicadas

Bacteria that form long-term intracellular associations with host cells lose many genes, a process that often results in tiny, gene-dense, and stable genomes. Paradoxically, the some of the same evolutionary processes that drive genome reduction and simplification may also cause genome expansion and complexification. A bacterial endosymbiont of cicadas, Hodgkinia cicadicola, exemplifies this paradox. In many cicada species, a single Hodgkinia lineage with a tiny, gene-dense genome has split into several interdependent cell and genome lineages. Each new Hodgkinia lineage encodes a unique subset of the ancestral unsplit genome in a complementary way, such that the collective gene contents of all lineages match the total found in the ancestral single genome. This splitting creates genetically distinct Hodgkinia cells that must function together to carry out basic cellular processes. It also creates a gene dosage problem where some genes are encoded by only a small fraction of cells while others are much more abundant. Here, by sequencing DNA and RNA of Hodgkinia from different cicada species with different amounts of splitting-along with its structurally stable, unsplit partner endosymbiont Sulcia muelleri-we show that Hodgkinia does not transcriptionally compensate to rescue the wildly unbalanced gene and genome ratios that result from lineage splitting. We also find that Hodgkinia has a reduced capacity for basic transcriptional control independent of the splitting process. Our findings reveal another layer of degeneration further pushing the limits of canonical molecular and cell biology in Hodgkinia and may partially explain its propensity to go extinct through symbiont replacement.

Genomic Insights of Alnus-Infective Frankia Strains Reveal Unique Genetic Features and New Evidence on Their Host-Restricted Lifestyle

The present study aimed to use comparative genomics to explore the relationships between Frankia and actinorhizal plants using a data set made of 33 Frankia genomes. The determinants of host specificity were first explored for "Alnus-infective strains" (i.e., Frankia strains belonging to Cluster Ia). Several genes were specifically found in these strains, including an agmatine deiminase which could possibly be involved in various functions as access to nitrogen sources, nodule organogenesis or plant defense. Within "Alnus-infective strains", Sp+ Frankia genomes were compared to Sp- genomes in order to elucidate the narrower host specificity of Sp+ strains (i.e., Sp+ strains being capable of in planta sporulation, unlike Sp- strains). A total of 88 protein families were lost in the Sp+ genomes. The lost genes were related to saprophytic life (transcriptional factors, transmembrane and secreted proteins), reinforcing the proposed status of Sp+ as obligatory symbiont. The Sp+ genomes were also characterized by a loss of genetic and functional paralogs, highlighting a reduction in functional redundancy (e.g., hup genes) or a possible loss of function related to a saprophytic lifestyle (e.g., genes involved in gas vesicle formation or recycling of nutrients).

Regulatory Interplay between RNase III and Antisense RNAs in E. coli: the Case of AsflhD and FlhD, Component of the Master Regulator of Motility

In order to respond to ever-changing environmental cues, bacteria display resilient regulatory mechanisms controlling gene expression. At the post-transcriptional level, this is achieved by a combination of RNA-binding proteins, such as ribonucleases (RNases), and regulatory RNAs, including antisense RNAs (asRNAs). Bound to their complementary mRNA, asRNAs are primary targets for the double-strand-specific endoribonuclease, RNase III. Taking advantage of our own and previously published transcriptomic data sets obtained in strains inactivated for RNase III, we selected several candidate asRNAs and confirmed the existence of RNase III-sensitive asRNAs for crp, ompR, phoP, and flhD genes, all encoding global regulators of gene expression in Escherichia coli. Using FlhD, a component of the master regulator of motility (FlhD₄C₂), as our model, we demonstrate that the asRNA AsflhD, transcribed from the coding sequence of flhD, is involved in the fine-tuning of flhD expression and thus participates in the control of motility.

A landscape of gene regulation in the parasitic amoebozoa Entamoeba spp

Entamoeba are amoeboid extracellular parasites that represent an important group of organisms for which the regulatory networks must be examined to better understand how genes and functional processes are interrelated. In this work, we inferred the gene regulatory networks (GRNs) in four Entamoeba species, E. histolytica, E. dispar, E. nuttalli, and E. invadens, and the GRN topological properties and the corresponding biological functions were evaluated. From these analyses, we determined that transcription factors (TFs) of E. histolytica, E. dispar, and E. nuttalli are associated mainly with the LIM family, while the TFs in E. invadens are associated with the RRM_1 family. In addition, we identified that EHI_044890 regulates 121 genes in E. histolytica, EDI_297980 regulates 284 genes in E. dispar, ENU1_120230 regulates 195 genes in E. nuttalli, and EIN_249270 regulates 257 genes in E. invadens. Finally, we identified that three types of processes, Macromolecule metabolic process, Cellular macromolecule metabolic process, and Cellular nitrogen compound metabolic process, are the main biological processes for each network. The results described in this work can be used as a basis for the study of gene regulation in these organisms.

Gene Regulatory Network Inference and Gene Module Regulating Virulence in Fusarium oxysporum

In this work, we inferred the gene regulatory network (GRN) of the fungus Fusarium oxysporum by using the regulatory networks of Aspergillus nidulans FGSC A4, Neurospora crassa OR74A, Saccharomyces cerevisiae S288c, and Fusarium graminearum PH-1 as templates for sequence comparisons. Topological properties to infer the role of transcription factors (TFs) and to identify functional modules were calculated in the GRN. From these analyzes, five TFs were identified as hubs, including FOXG_04688 and FOXG_05432, which regulate 2,404 and 1,864 target genes, respectively. In addition, 16 communities were identified in the GRN, where the largest contains 1,923 genes and the smallest contains 227 genes. Finally, the genes associated with virulence were extracted from the GRN and exhaustively analyzed, and we identified a giant module with ten TFs and 273 target genes, where the most highly connected node corresponds to the transcription factor FOXG_05265, homologous to the putative bZip transcription factor CPTF1 of Claviceps purpurea, which is involved in ergotism disease that affects cereal crops and grasses. The results described in this work can be used for the study of gene regulation in this organism and open the possibility to explore putative genes associated with virulence against their host.

Microbial adaptation to different environmental conditions: molecular perspective of evolved genetic and cellular systems

Microorganisms are ubiquitous on Earth and can inhabit almost every environment. In a complex heterogeneous environment or in face of ecological disturbance, the microbes adjust to fluctuating environmental conditions through a cascade of cellular and molecular systems. Their habitats differ from cold microcosms of Antarctica to the geothermal volcanic areas, terrestrial to marine, highly alkaline zones to the extremely acidic areas and freshwater to brackish water sources. The diverse ecological microbial niches are attributed to the versatile, adaptable nature under fluctuating temperature, nutrient availability and pH of the microorganisms. These organisms have developed a series of mechanisms to face the environmental changes and thereby keep their role in mediate important ecosystem functions. The underlying mechanisms of adaptable microbial nature are thoroughly investigated at the cellular, genetic and molecular levels. The adaptation is mediated by a spectrum of processes like natural selection, genetic recombination, horizontal gene transfer, DNA damage repair and pleiotropy-like events. This review paper provides the fundamentals insight into the microbial adaptability besides highlighting the molecular network of microbial adaptation under different environmental conditions.

RNA polymerases in strict endosymbiont bacteria with extreme genome reduction show distinct erosions that might result in limited and differential promoter recognition

Strict endosymbiont bacteria present high degree genome reduction, retain smaller proteins, and in some instances, lack complete functional domains compared to free-living counterparts. Until now, the mechanisms underlying these genetic reductions are not well understood. In this study, the conservation of RNA polymerases, the essential machinery for gene expression, is analyzed in endosymbiont bacteria with extreme genome reductions. We analyzed the RNA polymerase subunits to identify and define domains, subdomains, and specific amino acids involved in precise biological functions known in Escherichia coli. We also perform phylogenetic analysis and three-dimensional models over four lineages of endosymbiotic proteobacteria with the smallest genomes known to date: Candidatus Hodgkinia cicadicola, Candidatus Tremblaya phenacola, Candidatus Tremblaya Princeps, Candidatus Nasuia deltocephalinicola, and Candidatus Carsonella ruddii. We found that some Hodgkinia strains do not encode for the RNA polymerase α subunit. The rest encode genes for α, β, β', and σ subunits to form the RNA polymerase. However, 16% shorter, on average, respect their orthologous in E. coli. In the α subunit, the amino-terminal domain is the most conserved. Regarding the β and β' subunits, both the catalytic core and the assembly domains are the most conserved. However, they showed compensatory amino acid substitutions to adapt to changes in the σ subunit. Precisely, the most erosive diversity occurs within the σ subunit. We identified broad amino acid substitution even in those recognizing and binding to the -10-box promoter element. In an overall conceptual image, the RNA polymerase from Candidatus Nasuia conserved the highest similarity with Escherichia coli RNA polymerase and their σ70. It might be recognizing the two main promoter elements (-10 and -35) and the two promoter accessory elements (-10 extended and UP-element). In Candidatus Carsonella, the RNA polymerase could recognize all the promoter elements except the -10-box extended. In Candidatus Tremblaya and Hodgkinia, due to the α carboxyl-terminal domain absence, they might not recognize the UP-promoter element. We also identified the lack of the β flap-tip helix domain in most Hodgkinia's that suggests the inability to bind the -35-box promoter element.

Identifying Genes Devoted to the Cell Death Process in the Gene Regulatory Network of Ustilago maydis

Cell death is a process that can be divided into three morphological patterns: apoptosis, autophagy and necrosis. In fungi, cell death is induced in response to intracellular and extracellular perturbations, such as plant defense molecules, toxins and fungicides, among others. Ustilago maydis is a dimorphic fungus used as a model for pathogenic fungi of animals, including humans, and plants. Here, we reconstructed the transcriptional regulatory network of U. maydis, through homology inferences by using as templates the well-known gene regulatory networks (GRNs) of Saccharomyces cerevisiae, Aspergillus nidulans and Neurospora crassa. Based on this GRN, we identified transcription factors (TFs) as hubs and functional modules and calculated diverse topological metrics. In addition, we analyzed exhaustively the module related to cell death, with 60 TFs and 108 genes, where diverse cell proliferation, mating-type switching and meiosis, among other functions, were identified. To determine the role of some of these genes, we selected a set of 11 genes for expression analysis by qRT-PCR (sin3, rlm1, aif1, tdh3 [isoform A], tdh3 [isoform B], ald4, mca1, nuc1, tor1, ras1, and atg8) whose homologues in other fungi have been described as central in cell death. These genes were identified as downregulated at 72 h, in agreement with the beginning of the cell death process. Our results can serve as the basis for the study of transcriptional regulation, not only of the cell death process but also of all the cellular processes of U. maydis.

Determination and Dissection of DNA-Binding Specificity for the Thermus thermophilus HB8 Transcriptional Regulator TTHB099

Transcription factors (TFs) have been extensively researched in certain well-studied organisms, but far less so in others. Following the whole-genome sequencing of a new organism, TFs are typically identified through their homology with related proteins in other organisms. However, recent findings demonstrate that structurally similar TFs from distantly related bacteria are not usually evolutionary orthologs. Here we explore TTHB099, a cAMP receptor protein (CRP)-family TF from the extremophile Thermus thermophilus HB8. Using the in vitro iterative selection method Restriction Endonuclease Protection, Selection and Amplification (REPSA), we identified the preferred DNA-binding motif for TTHB099, 5′–TGT(A/g)NBSYRSVN(T/c)ACA–3′, and mapped potential binding sites and regulated genes within the T. thermophilus HB8 genome. Comparisons with expression profile data in TTHB099-deficient and wild type strains suggested that, unlike E. coli CRP (CRP_Ec), TTHB099 does not have a simple regulatory mechanism. However, we hypothesize that TTHB099 can be a dual-regulator similar to CRP_Ec.

Investigating the Role of FlhF Identifies Novel Interactions With Genes Involved in Flagellar Synthesis in Campylobacter jejuni

FlhF is a key protein required for complete flagellar synthesis, and its deletion results in the complete absence of a flagella and thus motility in Campylobacter jejuni. However, the specific mechanism still remains unknown. In this study, RNA-Seq, EMSAs, ChIP-qPCR and β-Galactosidase assays were performed to elucidate the novel interactions between FlhF and genes involved in flagellar synthesis. Results showed that FlhF has an overall influence on the transcription of flagellar genes with an flhF mutant displaying down-regulation of most flagellar related genes. FlhF can directly bind to the flgI promoter to regulate its expression, which has significant expression change in an flhF mutant. The possible binding site of FlhF to the flgI promoter was explored by continuously narrowing the flgI promoter region and performing further point mutations. Meanwhile, FlhF can directly bind to the promoters of rpoD, flgS, and fliA encoding early flagellin regulators, thereby directly or indirectly regulating the synthesis of class I, II, and III flagellar genes, respectively. Collectively, this study demonstrates that FlhF may directly regulate the transcription of flagellar genes by binding to their promoters as a transcriptional regulator, which will be helpful in understanding the mechanism of FlhF in flagellar biosynthetic and bacterial flagellation in general.

Trends in Symbiont-Induced Host Cellular Differentiation

Bacteria participate in a wide diversity of symbiotic associations with eukaryotic hosts that require precise interactions for bacterial recognition and persistence. Most commonly, host-associated bacteria interfere with host gene expression to modulate the immune response to the infection. However, many of these bacteria also interfere with host cellular differentiation pathways to create a hospitable niche, resulting in the formation of novel cell types, tissues, and organs. In both of these situations, bacterial symbionts must interact with eukaryotic regulatory pathways. Here, we detail what is known about how bacterial symbionts, from pathogens to mutualists, control host cellular differentiation across the central dogma, from epigenetic chromatin modifications, to transcription and mRNA processing, to translation and protein modifications. We identify four main trends from this survey. First, mechanisms for controlling host gene expression appear to evolve from symbionts co-opting cross-talk between host signaling pathways. Second, symbiont regulatory capacity is constrained by the processes that drive reductive genome evolution in host-associated bacteria. Third, the regulatory mechanisms symbionts exhibit correlate with the cost/benefit nature of the association. And, fourth, symbiont mechanisms for interacting with host genetic regulatory elements are not bound by native bacterial capabilities. Using this knowledge, we explore how the ubiquitous intracellular Wolbachia symbiont of arthropods and nematodes may modulate host cellular differentiation to manipulate host reproduction. Our survey of the literature on how infection alters gene expression in Wolbachia and its hosts revealed that, despite their intermediate-sized genomes, different strains appear capable of a wide diversity of regulatory manipulations. Given this and Wolbachia's diversity of phenotypes and eukaryotic-like proteins, we expect that many symbiont-induced host differentiation mechanisms will be discovered in this system.

Reverse Engineering of an Aspirin-Responsive Transcriptional Regulator in Escherichia coli

Bacterial transcription factors (TFs) are key devices for the engineering of complex circuits in many biotechnological applications, yet there are few well-characterized inducer-responsive TFs that could be used in the context of an animal or human host. We have deciphered the inducer recognition mechanism of two AraC/XylS regulators from Pseudomonas putida (BenR and XylS) for creating a novel expression system responsive to acetyl salicylate (i.e., aspirin). Using protein homology modeling and molecular docking with the cognate inducer benzoate and a suite of chemical analogues, we identified the conserved binding pocket of BenR and XylS. By means of site-directed mutagenesis, we identified a single amino acid position required for efficient inducer recognition and transcriptional activation. Whereas this modification in BenR abolishes protein activity, in XylS, it increases the response to several inducers, including acetyl salicylic acid, to levels close to those achieved by the canonical inducer. Moreover, by constructing chimeric proteins with swapped N-terminal domains, we created novel regulators with mixed promoter and inducer recognition profiles. As a result, a collection of engineered TFs was generated with an enhanced response to benzoate, 3-methylbenzoate, 2-methylbenzoate, 4-methylbenzoate, salicylic acid, aspirin, and acetylsalicylic acid molecules for eliciting gene expression in E. coli.

Dissecting the protein architecture of DNA-binding transcription factors in bacteria and archaea

Gene regulation at the transcriptional level is a central process in all organisms where DNA-binding transcription factors play a fundamental role. This class of proteins binds specifically at DNA sequences, activating or repressing gene expression as a function of the cell's metabolic status, operator context and ligand-binding status, among other factors, through the DNA-binding domain (DBD). In addition, TFs may contain partner domains (PaDos), which are involved in ligand binding and protein-protein interactions. In this work, we systematically evaluated the distribution, abundance and domain organization of DNA-binding TFs in 799 non-redundant bacterial and archaeal genomes. We found that the distributions of the DBDs and their corresponding PaDos correlated with the size of the genome. We also identified specific combinations between the DBDs and their corresponding PaDos. Within each class of DBDs there are differences in the actual angle formed at the dimerization interface, responding to the presence/absence of ligands and/or crystallization conditions, setting the orientation of the resulting helices and wings facing the DNA. Our results highlight the importance of PaDos as central elements that enhance the diversity of regulatory functions in all bacterial and archaeal organisms, and our results also demonstrate the role of PaDos in sensing diverse signal compounds. The highly specific interactions between DBDs and PaDos observed in this work, together with our structural analysis highlighting the difficulty in predicting both inter-domain geometry and quaternary structure, suggest that these systems appeared once and evolved with diverse duplication events in all the analysed organisms.

Reconstruction of Transcription Control Networks in Mollicutes by High-Throughput Identification of Promoters

Bacteria of the class Mollicutes have significantly reduced genomes and gene expression control systems. They are also efficient pathogens that can colonize a broad range of hosts including plants and animals. Despite their simplicity, Mollicutes demonstrate complex transcriptional responses to various conditions, which contradicts their reduction in gene expression regulation mechanisms. We analyzed the conservation and distribution of transcription regulators across the 50 Mollicutes species. The majority of the transcription factors regulate transport and metabolism, and there are four transcription factors that demonstrate significant conservation across the analyzed bacteria. These factors include repressors of chaperone HrcA, cell cycle regulator MraZ and two regulators with unclear function from the WhiA and YebC/PmpR families. We then used three representative species of the major clades of Mollicutes (Acholeplasma laidlawii, Spiroplasma melliferum, and Mycoplasma gallisepticum) to perform promoter mapping and activity quantitation. We revealed that Mollicutes evolved towards a promoter architecture simplification that correlates with a diminishing role of transcription regulation and an increase in transcriptional noise. Using the identified operons structure and a comparative genomics approach, we reconstructed the transcription control networks for these three species. The organization of the networks reflects the adaptation of bacteria to specific conditions and hosts.