Horizontal gene transfer and the evolution of transcriptional regulation in Escherichia coli

Infant age inversely correlates with gut carriage of resistance genes, reflecting modifications in microbial carbohydrate metabolism during early life

The infant gut microbiome is increasingly recognized as a reservoir of antibiotic resistance genes, yet the assembly of gut resistome in infants and its influencing factors remain largely unknown. We characterized resistome in 4132 metagenomes from 963 infants in six countries and 4285 resistance genes were observed. The inherent resistome pattern of healthy infants (N = 272) could be distinguished by two stages: a multicompound resistance phase (Months 0-7) and a tetracycline-mupirocin-β-lactam-dominant phase (Months 8-14). Microbial taxonomy explained 40.7% of the gut resistome of healthy infants, with Escherichia (25.5%) harboring the most resistance genes. In a further analysis with all available infants (N = 963), we found age was the strongest influencer on the resistome and was negatively correlated with the overall resistance during the first 3 years (p < 0.001). Using a random-forest approach, a set of 34 resistance genes could be used to predict age (R 2 = 68.0%). Leveraging microbial host inference analyses, we inferred the age-dependent assembly of infant resistome was a result of shifts in the gut microbiome, primarily driven by changes in taxa that disproportionately harbor resistance genes across taxa (e.g., Escherichia coli more frequently harbored resistance genes than other taxa). We performed metagenomic functional profiling and metagenomic assembled genome analyses whose results indicate that the development of gut resistome was driven by changes in microbial carbohydrate metabolism, with an increasing need for carbohydrate-active enzymes from Bacteroidota and a decreasing need for Pseudomonadota during infancy. Importantly, we observed increased acquired resistance genes over time, which was related to increased horizontal gene transfer in the developing infant gut microbiome. In summary, infant age was negatively correlated with antimicrobial resistance gene levels, reflecting a composition shift in the gut microbiome, likely driven by the changing need for microbial carbohydrate metabolism during early life.

Unveiling distinct genetic features in multidrug-resistant Escherichia coli isolated from mammary tissue and gut of mastitis induced mice

Escherichia coli is one of the major pathogens causing mastitis in lactating mammals. We hypothesized that E. coli from the gut and mammary glands may have similar genomic characteristics in the causation of mastitis. To test this hypothesis, we used whole genome sequencing to analyze two multidrug resistant E. coli strains isolated from mammary tissue (G2M6U) and fecal sample (G6M1F) of experimentally induced mastitis mice. Both strains showed resistance to multiple (>7) antibiotics such as oxacillin, aztreonam, nalidixic acid, streptomycin, gentamicin, cefoxitin, ampicillin, tetracycline, azithromycin and nitrofurantoin. The genome of E. coli G2M6U had 59 antimicrobial resistance genes (ARGs) and 159 virulence factor genes (VFGs), while the E. coli G6M1F genome possessed 77 ARGs and 178 VFGs. Both strains were found to be genetically related to many E. coli strains causing mastitis and enteric diseases originating from different hosts and regions. The G6M1F had several unique ARGs (e.g., QnrS1, sul2, tetA, tetR, emrK, blaTEM-1/105, and aph(6)-Id, aph(3″)-Ib) conferring resistance to certain antibiotics, whereas G2M6U had a unique heat-stable enterotoxin gene (astA) and 7192 single nucleotide polymorphisms. Furthermore, there were 43 and 111 unique genes identified in G2M6U and G6M1F genomes, respectively. These results indicate distinct differences in the genomic characteristics of E. coli strain G2M6U and G6M1F that might have important implications in the pathophysiology of mammalian mastitis, and treatment strategies for mastitis in dairy animals.

Learning from the unknown: exploring the range of bacterial functionality

Determining the repertoire of a microbe's molecular functions is a central question in microbial biology. Modern techniques achieve this goal by comparing microbial genetic material against reference databases of functionally annotated genes/proteins or known taxonomic markers such as 16S rRNA. Here, we describe a novel approach to exploring bacterial functional repertoires without reference databases. Our Fusion scheme establishes functional relationships between bacteria and assigns organisms to Fusion-taxa that differ from otherwise defined taxonomic clades. Three key findings of our work stand out. First, bacterial functional comparisons outperform marker genes in assigning taxonomic clades. Fusion profiles are also better for this task than other functional annotation schemes. Second, Fusion-taxa are robust to addition of novel organisms and are, arguably, able to capture the environment-driven bacterial diversity. Finally, our alignment-free nucleic acid-based Siamese Neural Network model, created using Fusion functions, enables finding shared functionality of very distant, possibly structurally different, microbial homologs. Our work can thus help annotate functional repertoires of bacterial organisms and further guide our understanding of microbial communities.

The Architecture of Metabolic Networks Constrains the Evolution of Microbial Resource Hierarchies

Microbial strategies for resource use are an essential determinant of their fitness in complex habitats. When facing environments with multiple nutrients, microbes often use them sequentially according to a preference hierarchy, resulting in well-known patterns of diauxic growth. In theory, the evolutionary diversification of metabolic hierarchies could represent a mechanism supporting coexistence and biodiversity by enabling temporal segregation of niches. Despite this ecologically critical role, the extent to which substrate preference hierarchies can evolve and diversify remains largely unexplored. Here, we used genome-scale metabolic modeling to systematically explore the evolution of metabolic hierarchies across a vast space of metabolic network genotypes. We find that only a limited number of metabolic hierarchies can readily evolve, corresponding to the most commonly observed hierarchies in genome-derived models. We further show how the evolution of novel hierarchies is constrained by the architecture of central metabolism, which determines both the propensity to change ranks between pairs of substrates and the effect of specific reactions on hierarchy evolution. Our analysis sheds light on the genetic and mechanistic determinants of microbial metabolic hierarchies, opening new research avenues to understand their evolution, evolvability, and ecology.

Early effects of gene duplication on the robustness and phenotypic variability of gene regulatory networks

Background

Research on gene duplication is abundant and comes from a wide range of approaches, from high-throughput analyses and experimental evolution to bioinformatics and theoretical models. Notwithstanding, a consensus is still lacking regarding evolutionary mechanisms involved in evolution through gene duplication as well as the conditions that affect them. We argue that a better understanding of evolution through gene duplication requires considering explicitly that genes do not act in isolation. It demands studying how the perturbation that gene duplication implies percolates through the web of gene interactions. Due to evolution’s contingent nature, the paths that lead to the final fate of duplicates must depend strongly on the early stages of gene duplication, before gene copies have accumulated distinctive changes.

Methods

Here we use a widely-known model of gene regulatory networks to study how gene duplication affects network behavior in early stages. Such networks comprise sets of genes that cross-regulate. They organize gene activity creating the gene expression patterns that give cells their phenotypic properties. We focus on how duplication affects two evolutionarily relevant properties of gene regulatory networks: mitigation of the effect of new mutations and access to new phenotypic variants through mutation.

Results

Among other observations, we find that those networks that are better at maintaining the original phenotype after duplication are usually also better at buffering the effect of single interaction mutations and that duplication tends to enhance further this ability. Moreover, the effect of mutations after duplication depends on both the kind of mutation and genes involved in it. We also found that those phenotypes that had easier access through mutation before duplication had higher chances of remaining accessible through new mutations after duplication.

Conclusion

Our results support that gene duplication often mitigates the impact of new mutations and that this effect is not merely due to changes in the number of genes. The work that we put forward helps to identify conditions under which gene duplication may enhance evolvability and robustness to mutations.

Bacterial plasmid-associated and chromosomal proteins have fundamentally different properties in protein interaction networks

Plasmids facilitate horizontal gene transfer, which enables the diversification of pathogens into new anatomical and environmental niches, implying that plasmid-encoded genes can cooperate well with chromosomal genes. We hypothesise that such mobile genes are functionally different to chromosomal ones due to this ability to encode proteins performing non-essential functions like antimicrobial resistance and traverse distinct host cells. The effect of plasmid-driven gene gain on protein-protein interaction network topology is an important question in this area. Moreover, the extent to which these chromosomally- and plasmid-encoded proteins interact with proteins from their own groups compared to the levels with the other group remains unclear. Here, we examined the incidence and protein-protein interactions of all known plasmid-encoded proteins across representative specimens from most bacteria using all available plasmids. We found that plasmid-encoded genes constitute ~ 0.65% of the total number of genes per bacterial sample, and that plasmid genes are preferentially associated with different species but had limited taxonomical power beyond this. Surprisingly, plasmid-encoded proteins had both more protein-protein interactions compared to chromosomal proteins, countering the hypothesis that genes with higher mobility rates should have fewer protein-level interactions. Nonetheless, topological analysis and investigation of the protein-protein interaction networks' connectivity and change in the number of independent components demonstrated that the plasmid-encoded proteins had limited overall impact in > 96% of samples. This paper assembled extensive data on plasmid-encoded proteins, their interactions and associations with diverse bacterial specimens that is available for the community to investigate in more detail.

Overlapping and Distinct Functions of the Paralogous PagR Regulators of Bacillus anthracis

The Bacillus anthracis pagA gene, encoding the protective antigen component of anthrax toxin, is part of a bicistronic operon on pXO1 that codes for its own repressor, PagR1. In addition to the pagAR1 operon, PagR1 regulates sap and eag, two chromosome genes encoding components of the surface layer, a mounting structure for surface proteins involved in virulence. Genomic studies have revealed a PagR1 paralog, PagR2, encoded by a gene on pXO2. The amino acid sequences of the paralogues are 71% identical and show similarity to the ArsR family of transcription regulators. We determined that the expression of either rPagR1 or rPagR2 in a ΔpagR1 pXO1⁺/pXO2^- (PagR1-PagR2) background repressed the expression of pagA, sap, eag, and a newly discovered target, atxA, encoding virulence activator AtxA. Despite the redundancy in PagR1 and PagR2 function, we determined that purified rPagR1 bound DNA corresponding to the control regions of all four target genes and existed as a dimer in cell lysates, whereas rPagR2 exhibited weak binding to the DNA of the pagA and atxA promoters, did not bind sap or eag promoter DNA, and did not appear as a dimer in cell lysates. A single amino acid change in PagR2, S81Y, designed to match the native Y81 of PagR1, allowed for DNA-binding to the sap and eag promoters. Moreover, the S81Y mutation allowed for the detection of PagR2 homomultimers in coaffinity purification experiments. Our results expand our knowledge of the roles of the paralogues in B. anthracis gene expression and provide a potential mechanistic basis for differences in the functions of these repressors. IMPORTANCE The protective antigen component of the anthrax toxin is essential for the delivery of the enzymatic components of the toxin into host target cells. The toxin genes and other virulence genes of B. anthracis are regulated by multiple trans-acting regulators that respond to a variety of host-related signals. PagR1, one such trans-acting regulator, connects the regulation of plasmid-encoded and chromosome-encoded virulence genes by controlling both protective antigen and surface layer protein expression. Whether PagR2, a paralog of PagR1, also functions as a trans-acting regulator was unknown. This work advances our knowledge of the complex model of virulence regulation in B. anthracis and furthers our understanding of the intriguing evolution of this pathogen.

System Principles Governing the Organization, Architecture, Dynamics, and Evolution of Gene Regulatory Networks

Synthetic biology aims to apply engineering principles for the rational, systematical design and construction of biological systems displaying functions that do not exist in nature or even building a cell from scratch. Understanding how molecular entities interconnect, work, and evolve in an organism is pivotal to this aim. Here, we summarize and discuss some historical organizing principles identified in bacterial gene regulatory networks. We propose a new layer, the concilion, which is the group of structural genes and their local regulators responsible for a single function that, organized hierarchically, coordinate a response in a way reminiscent of the deliberation and negotiation that take place in a council. We then highlight the importance that the network structure has, and discuss that the natural decomposition approach has unveiled the system-level elements shaping a common functional architecture governing bacterial regulatory networks. We discuss the incompleteness of gene regulatory networks and the need for network inference and benchmarking standardization. We point out the importance that using the network structural properties showed to improve network inference. We discuss the advances and controversies regarding the consistency between reconstructions of regulatory networks and expression data. We then discuss some perspectives on the necessity of studying regulatory networks, considering the interactions’ strength distribution, the challenges to studying these interactions’ strength, and the corresponding effects on network structure and dynamics. Finally, we explore the ability of evolutionary systems biology studies to provide insights into how evolution shapes functional architecture despite the high evolutionary plasticity of regulatory networks.

Determination of the two-component systems regulatory network reveals core and accessory regulations across Pseudomonas aeruginosa lineages

Pseudomonas aeruginosa possesses one of the most complex bacterial regulatory networks, which largely contributes to its success as a pathogen. However, most of its transcription factors (TFs) are still uncharacterized and the potential intra-species variability in regulatory networks has been mostly ignored so far. Here, we used DAP-seq to map the genome-wide binding sites of all 55 DNA-binding two-component systems (TCSs) response regulators (RRs) across the three major P. aeruginosa lineages. The resulting networks encompass about 40% of all genes in each strain and contain numerous new regulatory interactions across most major physiological processes. Strikingly, about half of the detected targets are specific to only one or two strains, revealing a previously unknown large functional diversity of TFs within a single species. Three main mechanisms were found to drive this diversity, including differences in accessory genome content, as exemplified by the strain-specific plasmid in IHMA87 outlier strain which harbors numerous binding sites of conserved chromosomally-encoded RRs. Additionally, most RRs display potential auto-regulation or RR-RR cross-regulation, bringing to light the vast complexity of this network. Overall, we provide the first complete delineation of the TCSs regulatory network in P. aeruginosa that will represent an important resource for future studies on this pathogen.

Transcription Inhibitors with XRE DNA-Binding and Cupin Signal-Sensing Domains Drive Metabolic Diversification in Pseudomonas

Bacteria of the Pseudomonas genus, including the major human pathogen Pseudomonas aeruginosa, are known for their complex regulatory networks and high number of transcription factors, which contribute to their impressive adaptive ability. However, even in the most studied species, most of the regulators are still uncharacterized.

Transcription factors (TFs) are instrumental in the bacterial response to new environmental conditions. They can act as direct signal sensors and subsequently induce changes in gene expression leading to physiological adaptation. Here, by combining transcriptome sequencing (RNA-seq) and cistrome determination (DAP-seq), we studied a family of eight TFs in Pseudomonas aeruginosa. This family, encompassing TFs with XRE-like DNA-binding and cupin signal-sensing domains, includes the metabolic regulators ErfA, PsdR, and PauR and five so-far-unstudied TFs. The genome-wide delineation of their regulons identified 39 regulatory interactions with genes mostly involved in metabolism. We found that the XRE-cupin TFs are inhibitors of their neighboring genes, forming local, functional units encoding proteins with functions in condition-specific metabolic pathways. Growth phenotypes of isogenic mutants highlighted new roles for PauR and PA0535 in polyamines and arginine metabolism. The phylogenetic analysis of this family of regulators across the bacterial kingdom revealed a wide diversity of such metabolic regulatory modules and identified species with potentially higher metabolic versatility. Numerous genes encoding uncharacterized XRE-cupin TFs were found near metabolism-related genes, illustrating the need of further systematic characterization of transcriptional regulatory networks in order to better understand the mechanisms of bacterial adaptation to new environments.

IMPORTANCE Bacteria of the Pseudomonas genus, including the major human pathogen Pseudomonas aeruginosa, are known for their complex regulatory networks and high number of transcription factors, which contribute to their impressive adaptive ability. However, even in the most studied species, most of the regulators are still uncharacterized. With the recent advances in high-throughput sequencing methods, it is now possible to fill this knowledge gap and help the understanding of how bacteria adapt and thrive in new environments. By leveraging these methods, we provide an example of a comprehensive analysis of an entire family of transcription factors and bring new insights into metabolic and regulatory adaptation in the Pseudomonas genus.

CoryneRegNet 7, the reference database and analysis platform for corynebacterial gene regulatory networks

We present the newest version of CoryneRegNet, the reference database for corynebacterial regulatory interactions, available at www.exbio.wzw.tum.de/coryneregnet/. The exponential growth of next-generation sequencing data in recent years has allowed a better understanding of bacterial molecular mechanisms. Transcriptional regulation is one of the most important mechanisms for bacterial adaptation and survival. These mechanisms may be understood via an organism's network of regulatory interactions. Although the Corynebacterium genus is important in medical, veterinary and biotechnological research, little is known concerning the transcriptional regulation of these bacteria. Here, we unravel transcriptional regulatory networks (TRNs) for 224 corynebacterial strains by utilizing genome-scale transfer of TRNs from four model organisms and assigning statistical significance values to all predicted regulations. As a result, the number of corynebacterial strains with TRNs increased twenty times and the back-end and front-end were reimplemented to support new features as well as future database growth. CoryneRegNet 7 is the largest TRN database for the Corynebacterium genus and aids in elucidating transcriptional mechanisms enabling adaptation, survival and infection.

Dynamics of genetic variation in transcription factors and its implications for the evolution of regulatory networks in Bacteria

The evolution of regulatory networks in Bacteria has largely been explained at macroevolutionary scales through lateral gene transfer and gene duplication. Transcription factors (TF) have been found to be less conserved across species than their target genes (TG). This would be expected if TFs accumulate mutations faster than TGs. This hypothesis is supported by several lab evolution studies which found TFs, especially global regulators, to be frequently mutated. Despite these studies, the contribution of point mutations in TFs to the evolution of regulatory network is poorly understood. We tested if TFs show greater genetic variation than their TGs using whole-genome sequencing data from a large collection of Escherichia coli isolates. TFs were less diverse than their TGs across natural isolates, with TFs of large regulons being more conserved. In contrast, TFs showed higher mutation frequency in adaptive laboratory evolution experiments. However, over long-term laboratory evolution spanning 60 000 generations, mutation frequency in TFs gradually declined after a rapid initial burst. Extrapolating the dynamics of genetic variation from long-term laboratory evolution to natural populations, we propose that point mutations, conferring large-scale gene expression changes, may drive the early stages of adaptation but gene regulation is subjected to stronger purifying selection post adaptation.

Analyses of Livestock-Associated Staphylococcus aureus Pan-Genomes Suggest Virulence Is Not Primary Interest in Evolution of Its Genome

Staphylococcus aureus is not only part of normal flora but also an opportunistic pathogen relevant to microbial genomics, public health, and veterinary medicine. In addition to being a well-known human pathogen, S. aureus causes various infections in economically important livestock animals such as cows, sheep, goats, and pigs. There are very few studies that have examined the pan-genome of S. aureus or the host-specific strains' pan-genomes. We report on livestock-associated S. aureus' (LA-SA) pan-genome and suggest that virulence is not the primary interest in evolution of its genome. LA-SA' complete genomes were retrieved from the NCBI and pan-genome was constructed by high-speed Roary pipeline. The pan-genome size was 4637 clusters, whereas 42.46% of the pan-genome was associated with the core genome. We found 1268 genes were associated with the strain-unique genome, and the remaining 1432 cluster with the accessory genome. COG (clusters of orthologous group of proteins) analysis of the core genes revealed 34% of clusters related to metabolism responsible for amino acid and inorganic ion transport (COG categories E and P), followed by carbohydrate metabolism (category G). Virulent gene analysis revealed the core genes responsible for antiphagocytosis and iron uptake. The fluidity of pan-genome was calculated as 0.082 ± 0.025. Importantly, the positive selection analysis suggested a slower rate of evolution among the LA-SA genomes. We call for comparative microbial and pan-genome research between human and LA-SA that can help further understand the evolution of virulence and thus inform future microbial diagnostics and drug discovery.

Species-specific recruitment of transcription factors dictates toxin expression

Tight and coordinate regulation of virulence determinants is essential for bacterial biology and involves dynamic shaping of transcriptional regulatory networks during evolution. The horizontally transferred two-partner secretion system ExlB-ExlA is instrumental in the virulence of different Pseudomonas species, ranging from soil- and plant-dwelling biocontrol agents to the major human pathogen Pseudomonas aeruginosa. Here, we identify a Cro/CI-like repressor, named ErfA, which together with Vfr, a CRP-like activator, controls exlBA expression in P. aeruginosa. The characterization of ErfA regulon across P. aeruginosa subfamilies revealed a second conserved target, the ergAB operon, with functions unrelated to virulence. To gain insights into this functional dichotomy, we defined the pan-regulon of ErfA in several Pseudomonas species and found ergAB as the sole conserved target of ErfA. The analysis of 446 exlBA promoter sequences from all exlBA+ genomes revealed a wide variety of regulatory sequences, as ErfA- and Vfr-binding sites were found to have evolved specifically in P. aeruginosa and nearly each species carries different regulatory sequences for this operon. We propose that the emergence of different regulatory cis-elements in the promoters of horizontally transferred genes is an example of plasticity of regulatory networks evolving to provide an adapted response in each individual niche.

Excessive Promoters as Silencers of Genes Horizontally Acquired by Escherichia coli

Horizontally acquired genes are usually transcriptionally inactive, although most of them are associated with genomic loci enriched with promoter-like sequences forming “promoter islands.” We hypothesized that lateral DNA transfer induces local mutagenesis, accumulating AT base pairs and creating promoter-like sequences, whose occupancy with RNA polymerase and a specific silencer H-NS suppresses the transcription of foreign genes. Error-prone mutagenesis was implemented for the “promoter island” of a foreign gene appY and the promoter region of an inherent gene dps. Derivatives with changed transcriptional activity were selected using a reporter plasmid pET28_eGFP. Only one cycle of mutagenesis with negative selection suppressed the activity of the main dps promoter to the background level due to a single substitution in its -10 element, while positive selection gave a sequence with improved -35 element, thus testifying feasibility of the approach. The same suppression for appY was achieved by three cycles, while eightfold transcription activation required nine iterations of mutagenesis. In both cases, the number of potential start points decreased resulting in an ordinary regulatory region with only one dominant promoter in the case of positive selection. Efficiency of H-NS binding remained virtually unchanged in all mutant constructs. Based on these findings we conclude that excessive promoters can adversely affect transcription by providing a platform for interference between several RNA polymerase molecules, which can act as a silencer at promoter-dense regions.

Methodology of predicting novel key regulators in ovarian cancer network: a network theoretical approach

Background

Identification of key regulator/s in ovarian cancer (OC) network is important for potential drug target and prevention from this cancer. This study proposes a method to identify the key regulators of this network and their importance.

Methods

The protein-protein interaction (PPI) network of ovarian cancer (OC) is constructed from curated 6 hundred genes from standard six important ovarian cancer databases (some of the genes are experimentally verified). We proposed a method to identify key regulators (KRs) from the complex ovarian cancer network based on the tracing of backbone hubs, which participate at all levels of organization, characterized by Newmann-Grivan community finding method. Knockout experiment, constant Potts model and survival analysis are done to characterize the importance of the key regulators in regulating the network.

Results

The PPI network of ovarian cancer is found to obey hierarchical scale free features organized by topology of heterogeneous modules coordinated by diverse leading hubs. The network and modular structures are devised by fractal rules with the absence of centrality-lethality rule, to enhance the efficiency of signal processing in the network and constituting loosely connected modules. Within the framework of network theory, we device a method to identify few key regulators (KRs) from a huge number of leading hubs, that are deeply rooted in the network, serve as backbones of it and key regulators from grassroots level to complete network structure. Using this method we could able to identify five key regulators, namely, AKT1, KRAS, EPCAM, CD44 and MCAM, out of which AKT1 plays central role in two ways, first it serves as main regulator of ovarian cancer network and second serves as key cross-talk agent of other key regulators, but exhibits disassortive property. The regulating capability of AKT1 is found to be highest and that of MCAM is lowest.

Conclusions

The popularities of these key hubs change in an unpredictable way at different levels of organization and absence of these hubs cause massive amount of wiring energy/rewiring energy that propagate over all the network. The network compactness is found to increase as one goes from top level to bottom level of the network organization.

Evolutionary and functional patterns of shared gene neighbourhood in fungi

Gene clusters comprise genomically co-localized and potentially co-regulated genes that tend to be conserved across species. In eukaryotes, multiple examples of metabolic gene clusters are known, particularly among fungi and plants. However, little is known about how gene clustering patterns vary among taxa or with respect to functional roles. Furthermore, mechanisms of the formation, maintenance and evolution of gene clusters remain unknown. We surveyed 341 fungal genomes to discover gene clusters shared by different species, independently of their functions. We inferred 12,120 cluster families, which comprised roughly one third of the gene space and were enriched in genes associated with diverse cellular functions. Additionally, most clusters did not encode transcription factors, suggesting that they are regulated distally. We used phylogenomics to characterize the evolutionary history of these clusters. We found that most clusters originated once and were transmitted vertically, coupled to differential loss. However, convergent evolution-that is, independent appearance of the same cluster-was more prevalent than anticipated. Finally, horizontal gene transfer of entire clusters was somewhat restricted, with the exception of those associated with secondary metabolism. Altogether, our results provide insights on the evolution of gene clustering as well as a broad catalogue of evolutionarily conserved gene clusters whose function remains to be elucidated.

Draft Genome Sequence of Acidovorax sp. Strain NB1, Isolated from a Nitrite-Oxidizing Enrichment Culture

Here, we report the draft genome sequence of Acidovorax sp. strain NB1, isolated from an enrichment culture of nitrite-oxidizing bacteria (NOB). Genes involved in denitrification were found in the draft genome of NB1. The closest strain to NB1 based on genomic relatedness is Acidovorax sp. strain GW101-3H11, with 91.5% average nucleotide identity.

Here, we report the draft genome sequence of Acidovorax sp. strain NB1, isolated from an enrichment culture of nitrite-oxidizing bacteria (NOB). Genes involved in denitrification were found in the draft genome of NB1. The closest strain to NB1 based on genomic relatedness is Acidovorax sp. strain GW101-3H11, with 91.5% average nucleotide identity.

Life at Home and on the Roam: Genomic Adaptions Reflect the Dual Lifestyle of an Intracellular, Facultative Symbiont

"Candidatus Synechococcus feldmannii" is a facultative intracellular symbiont of the Atlanto-Mediterranean sponge Petrosia ficiformis. Genomic information of sponge-associated cyanobacteria derives thus far from the obligate and extracellular symbiont "Candidatus Synechococcus spongiarum." Here we utilized a differential methylation-based approach for bacterial DNA enrichment combined with metagenomics to obtain the first draft genomes of "Ca. Synechococcus feldmannii." By comparative genomics, we revealed that some genomic features (e.g., iron transport mediated by siderophores, eukaryotic-like proteins, and defense mechanisms, like CRISPR-Cas [clustered regularly interspaced short palindromic repeats-associated proteins]) are unique to both symbiont types and absent or rare in the genomes of taxonomically related free-living cyanobacteria. These genomic features likely enable life under the conditions found inside the sponge host. Interestingly, there are many genomic features that are shared by "Ca. Synechococcus feldmannii" and free-living cyanobacteria, while they are absent in the obligate symbiont "Ca. Synechococcus spongiarum." These include genes related to cell surface structures, genetic regulation, and responses to environmental stress, as well as the composition of photosynthetic genes and DNA metabolism. We speculate that the presence of these genes confers on "Ca. Synechococcus feldmannii" its facultative nature (i.e., the ability to respond to a less stable environment when free-living). Our comparative analysis revealed that distinct genomic features depend on the nature of the symbiotic interaction: facultative and intracellular versus obligate and extracellular. IMPORTANCE Given the evolutionary position of sponges as one of the earliest phyla to depart from the metazoan stem lineage, studies on their distinct and exceptionally diverse microbial communities should yield a better understanding of the origin of animal-bacterium interactions. While genomes of several extracellular sponge symbionts have been published, the intracellular symbionts have, so far, been elusive. Here we compare the genomes of two unicellular cyanobacterial sponge symbionts that share an ancestor but followed different evolutionary paths-one became intracellular and the other extracellular. Counterintuitively, the intracellular cyanobacteria are facultative, while the extracellular ones are obligate. By sequencing the genomes of the intracellular cyanobacteria and comparing them to the genomes of the extracellular symbionts and related free-living cyanobacteria, we show how three different cyanobacterial lifestyles are reflected by adaptive genomic features.

Comparative Genomics of Nitrogen Cycling Pathways in Bacteria and Archaea

Despite the explosion of metagenomic sequencing data, using -omics data to predict environmental biogeochemistry remains a challenge. One or a few genes (referred to as marker genes) in a metabolic pathway of interest in meta-omic data are typically used to represent the prevalence of a biogeochemical reaction. This approach often fails to demonstrate a consistent relationship between gene abundance and an ecosystem process rate. One reason this may occur is if a marker gene is not a good representative of a complete pathway. Here, we map the presence of 11 nitrogen (N)-cycling pathways in over 6000 complete bacterial and archaeal genomes using the Integrated Microbial Genomes database. Incomplete N-cycling pathways occurred in 39% of surveyed archaeal and bacterial species revealing a weakness in current marker-gene analyses. Furthermore, we found that most organisms have limited ability to utilize inorganic N in multiple oxidation states. This suggests that inter-organism exchange of inorganic N compounds is common, highlighting the importance of both community composition and spatial structure in determining the extent of recycling versus loss in an ecosystem.

Network Pharmacology Deciphering Mechanisms of Volatiles of Wendan Granule for the Treatment of Alzheimer's Disease

Objective

To explore the mechanisms of the volatiles of Wendan granule (WDG) for the treatment of Alzheimer's disease, network pharmacology method integrating absorption, distribution, metabolism, and excretion (ADME) screening, target fishing, network constructing, pathway analysing, and correlated diseases prediction was applied.

Methods

Twelve small molecular compounds of WDG were selected as the objects from 74 volatiles with the relative abundances above 2 %, and their ADME parameters were collected from Traditional Chinese Medicine Systems Pharmacology platform (TCMSP), and the corresponding targets, genes, pathways, and diseases were predicted according to the data provided by TCMSP, DrugBank, Uniport, and the Database for Annotation, Visualization, and Integrated Discovery (DAVID). Then the related pathways and correlation analysis were explored by the Kyoto Encyclopedia and Genomes (KEGG) database. Finally, the networks of compound target, target pathway, and pathway disease of WDG were constructed by Cytoscape software.

Results

Twelve compounds interacted with 49 targets, of which top three targets were gamma-aminobutyric acid receptor subunit alpha-1 (GABRA1), prostaglandin G/H synthase 2 (PGHS-2), and sodium-dependent noradrenaline transporter. Interestingly, these targets were highly associated with depression, insomnia, and Alzheimer's disease that mainly corresponded to mental and emotional illnesses.

Conclusion

The integrated network pharmacology method provides precise probe to illuminate the molecular mechanisms of the main volatiles of WDG for relieving senile dementia related syndromes, which will also facilitate the application of traditional Chinese medicine as an alternative or supplementary to conventional treatments of AD, as well as follow-up studies such as upgrading the quality standard of clinically applied herbal medicine and novel drug development.

Transcriptional Regulation in Archaea: From Individual Genes to Global Regulatory Networks

Archaea are major contributors to biogeochemical cycles, possess unique metabolic capabilities, and resist extreme stress. To regulate the expression of genes encoding these unique programs, archaeal cells use gene regulatory networks (GRNs) composed of transcription factor proteins and their target genes. Recent developments in genetics, genomics, and computational methods used with archaeal model organisms have enabled the mapping and prediction of global GRN structures. Experimental tests of these predictions have revealed the dynamical function of GRNs in response to environmental variation. Here, we review recent progress made in this area, from investigating the mechanisms of transcriptional regulation of individual genes to small-scale subnetworks and genome-wide global networks. At each level, archaeal GRNs consist of a hybrid of bacterial, eukaryotic, and uniquely archaeal mechanisms. We discuss this theme from the perspective of the role of individual transcription factors in genome-wide regulation, how these proteins interact to compile GRN topological structures, and how these topologies lead to emergent, high-level GRN functions. We conclude by discussing how systems biology approaches are a fruitful avenue for addressing remaining challenges, such as discovering gene function and the evolution of GRNs.

Inter-replicon Gene Flow Contributes to Transcriptional Integration in the Sinorhizobium meliloti Multipartite Genome

Integration of newly acquired genes into existing regulatory networks is necessary for successful horizontal gene transfer (HGT). Ten percent of bacterial species contain at least two DNA replicons over 300 kilobases in size, with the secondary replicons derived predominately through HGT. The Sinorhizobium meliloti genome is split between a 3.7 Mb chromosome, a 1.7 Mb chromid consisting largely of genes acquired through ancient HGT, and a 1.4 Mb megaplasmid consisting primarily of recently acquired genes. Here, RNA-sequencing is used to examine the transcriptional consequences of massive, synthetic genome reduction produced through the removal of the megaplasmid and/or the chromid. Removal of the pSymA megaplasmid influenced the transcription of only six genes. In contrast, removal of the chromid influenced expression of ∼8% of chromosomal genes and ∼4% of megaplasmid genes. This was mediated in part by the loss of the ETR DNA region whose presence on pSymB is due to a translocation from the chromosome. No obvious functional bias among the up-regulated genes was detected, although genes with putative homologs on the chromid were enriched. Down-regulated genes were enriched in motility and sensory transduction pathways. Four transcripts were examined further, and in each case the transcriptional change could be traced to loss of specific pSymB regions. In particularly, a chromosomal transporter was induced due to deletion of bdhA likely mediated through 3-hydroxybutyrate accumulation. These data provide new insights into the evolution of the multipartite bacterial genome, and more generally into the integration of horizontally acquired genes into the transcriptome.

Phylogenetic diversification and developmental implications of poly-(R)-3-hydroxyalkanoate gene cluster assembly in prokaryotes

Many polyhydroxyalkanoates (PHAs) system genes, such as phaC, phaA, phaB, phaR, phaP and phaZ, are often found to be organised in the form of operon-like clusters. In this study, a genome survey was performed to identify such clustered PHA systems among 256 prokaryotic organisms. These data were then used to generate a comprehensive 16S rRNA species tree depicting the phylogenetic distribution of the observed clusters with diverse gene arrangements. In addition, the gene occurrences and physical linkages between PHA system genes were quantitatively estimated. From this, we identified a centrally connected hub gene, i.e. the phaC gene of PHA. Furthermore, a comparative investigation was performed between the clusters of PHA and glycogen, which decoded the role of the hub gene in the cluster organisation of both systems. Together, these findings suggest that the highly connected hub gene might contribute substantively towards the organisation and maintenance of the gene network connectivity in the clusters, particularly in the storage reserve systems.

Kullback Leibler divergence in complete bacterial and phage genomes

The amino acid content of the proteins encoded by a genome may predict the coding potential of that genome and may reflect lifestyle restrictions of the organism. Here, we calculated the Kullback–Leibler divergence from the mean amino acid content as a metric to compare the amino acid composition for a large set of bacterial and phage genome sequences. Using these data, we demonstrate that (i) there is a significant difference between amino acid utilization in different phylogenetic groups of bacteria and phages; (ii) many of the bacteria with the most skewed amino acid utilization profiles, or the bacteria that host phages with the most skewed profiles, are endosymbionts or parasites; (iii) the skews in the distribution are not restricted to certain metabolic processes but are common across all bacterial genomic subsystems; (iv) amino acid utilization profiles strongly correlate with GC content in bacterial genomes but very weakly correlate with the G+C percent in phage genomes. These findings might be exploited to distinguish coding from non-coding sequences in large data sets, such as metagenomic sequence libraries, to help in prioritizing subsequent analyses.

Horizontally Acquired Biosynthesis Genes Boost Coxiella burnetii's Physiology

Coxiella burnetii, the etiologic agent of acute Q fever and chronic endocarditis, has a unique biphasic life cycle, which includes a metabolically active intracellular form that occupies a large lysosome-derived acidic vacuole. C. burnetii is the only bacterium known to thrive within such an hostile intracellular niche, and this ability is fundamental to its pathogenicity; however, very little is known about genes that facilitate Coxiella's intracellular growth. Recent studies indicate that C. burnetii evolved from a tick-associated ancestor and that the metabolic capabilities of C. burnetii are different from that of Coxiella-like bacteria found in ticks. Horizontally acquired genes that allow C. burnetii to infect and grow within mammalian cells likely facilitated the host shift; however, because of its obligate intracellular replication, C. burnetii would have lost most genes that have been rendered redundant due to the availability of metabolites within the host cell. Based on these observations, we reasoned that horizontally derived biosynthetic genes that have been retained in the reduced genome of C. burnetii are ideal candidates to begin to uncover its intracellular metabolic requirements. Our analyses identified a large number of putative foreign-origin genes in C. burnetii, including tRNA^Glu2 that is potentially required for heme biosynthesis, and genes involved in the production of lipopolysaccharide—a virulence factor, and of critical metabolites such as fatty acids and biotin. In comparison to wild-type C. burnetii, a strain that lacks tRNA^Glu2 exhibited reduced growth, indicating its importance to Coxiella's physiology. Additionally, by using chemical agents that block heme and biotin biosyntheses, we show that these pathways are promising targets for the development of new anti-Coxiella therapies.

Control of MarRAB Operon in Escherichia coli via Autoactivation and Autorepression

Choice of network topology for gene regulation has been a question of interest for a long time. How do simple and more complex topologies arise? In this work, we analyze the topology of the marRAB operon in Escherichia coli, which is associated with control of expression of genes associated with conferring resistance to low-level antibiotics to the bacterium. Among the 2102 promoters in E. coli, the marRAB promoter is the only one that encodes for an autoactivator and an autorepressor. What advantages does this topology confer to the bacterium? In this work, we demonstrate that, compared to control by a single regulator, the marRAB regulatory arrangement has the least control cost associated with modulating gene expression in response to environmental stimuli. In addition, the presence of dual regulators allows the regulon to exhibit a diverse range of dynamics, a feature that is not observed in genes controlled by a single regulator.

Comparative genome-scale modelling of Staphylococcus aureus strains identifies strain-specific metabolic capabilities linked to pathogenicity

Staphylococcus aureus is a preeminent bacterial pathogen capable of colonizing diverse ecological niches within its human host. We describe here the pangenome of S. aureus based on analysis of genome sequences from 64 strains of S. aureus spanning a range of ecological niches, host types, and antibiotic resistance profiles. Based on this set, S. aureus is expected to have an open pangenome composed of 7,411 genes and a core genome composed of 1,441 genes. Metabolism was highly conserved in this core genome; however, differences were identified in amino acid and nucleotide biosynthesis pathways between the strains. Genome-scale models (GEMs) of metabolism were constructed for the 64 strains of S. aureus These GEMs enabled a systems approach to characterizing the core metabolic and panmetabolic capabilities of the S. aureus species. All models were predicted to be auxotrophic for the vitamins niacin (vitamin B3) and thiamin (vitamin B1), whereas strain-specific auxotrophies were predicted for riboflavin (vitamin B2), guanosine, leucine, methionine, and cysteine, among others. GEMs were used to systematically analyze growth capabilities in more than 300 different growth-supporting environments. The results identified metabolic capabilities linked to pathogenic traits and virulence acquisitions. Such traits can be used to differentiate strains responsible for mild vs. severe infections and preference for hosts (e.g., animals vs. humans). Genome-scale analysis of multiple strains of a species can thus be used to identify metabolic determinants of virulence and increase our understanding of why certain strains of this deadly pathogen have spread rapidly throughout the world.

Evolutionary change and phylogenetic relationships in light of horizontal gene transfer

Horizontal gene transfer has, over the past 25 years, become a part of evolutionary thinking. In the present paper I discuss horizontal gene transfer (HGT) in relation to contingency, natural selection, evolutionary change speed and the Tree-of-Life endeavour, with the aim of contributing to the understanding of the role of HGT in evolutionary processes. In addition, the challenges that HGT imposes on the current view of evolution are emphasized.

Conserved Units of Co-Expression in Bacterial Genomes: An Evolutionary Insight into Transcriptional Regulation

Genome-wide measurements of transcriptional activity in bacteria indicate that the transcription of successive genes is strongly correlated beyond the scale of operons. Here, we analyze hundreds of bacterial genomes to identify supra-operonic segments of genes that are proximal in a large number of genomes. We show that these synteny segments correspond to genomic units of strong transcriptional co-expression. Structurally, the segments contain operons with specific relative orientations (co-directional or divergent) and nucleoid-associated proteins are found to bind at their boundaries. Functionally, operons inside a same segment are highly co-expressed even in the apparent absence of regulatory factors at their promoter regions. Remote operons along DNA can also be co-expressed if their corresponding segments share a transcriptional or sigma factor, without requiring these factors to bind directly to the promoters of the operons. As evidence that these results apply across the bacterial kingdom, we demonstrate them both in the Gram-negative bacterium Escherichia coli and in the Gram-positive bacterium Bacillus subtilis. The underlying process that we propose involves only RNA-polymerases and DNA: it implies that the transcription of an operon mechanically enhances the transcription of adjacent operons. In support of a primary role of this regulation by facilitated co-transcription, we show that the transcription en bloc of successive operons as a result of transcriptional read-through is strongly and specifically enhanced in synteny segments. Finally, our analysis indicates that facilitated co-transcription may be evolutionary primitive and may apply beyond bacteria.

HosA, a MarR Family Transcriptional Regulator, Represses Nonoxidative Hydroxyarylic Acid Decarboxylase Operon and Is Modulated by 4-Hydroxybenzoic Acid

Members of the Multiple antibiotic resistance Regulator (MarR) family of DNA binding proteins regulate transcription of a wide array of genes required for virulence and pathogenicity of bacteria. The present study reports the molecular characterization of HosA (Homologue of SlyA), a MarR protein, with respect to its target gene, DNA recognition motif, and nature of its ligand. Through a comparative genomics approach, we demonstrate that hosA is in synteny with nonoxidative hydroxyarylic acid decarboxylase (HAD) operon and is present exclusively within the mutS-rpoS polymorphic region in nine different genera of Enterobacteriaceae family. Using molecular biology and biochemical approach, we demonstrate that HosA binds to a palindromic sequence downstream to the transcription start site of divergently transcribed nonoxidative HAD operon and represses its expression. Furthermore, in silico analysis showed that the recognition motif for HosA is highly conserved in the upstream region of divergently transcribed operon in different genera of Enterobacteriaceae family. A systematic chemical search for the physiological ligand revealed that 4-hydroxybenzoic acid (4-HBA) interacts with HosA and derepresses HosA mediated repression of the nonoxidative HAD operon. Based on our study, we propose a model for molecular mechanism underlying the regulation of nonoxidative HAD operon by HosA in Enterobacteriaceae family.

Following the Footsteps of Chlamydial Gene Regulation

Regulation of gene expression ensures an organism responds to stimuli and undergoes proper development. Although the regulatory networks in bacteria have been investigated in model microorganisms, nearly nothing is known about the evolution and plasticity of these networks in obligate, intracellular bacteria. The phylum Chlamydiae contains a vast array of host-associated microbes, including several human pathogens. The Chlamydiae are unique among obligate, intracellular bacteria as they undergo a complex biphasic developmental cycle in which large swaths of genes are temporally regulated. Coupled with the low number of transcription factors, these organisms offer a model to study the evolution of regulatory networks in intracellular organisms. We provide the first comprehensive analysis exploring the diversity and evolution of regulatory networks across the phylum. We utilized a comparative genomics approach to construct predicted coregulatory networks, which unveiled genus- and family-specific regulatory motifs and architectures, most notably those of virulence-associated genes. Surprisingly, our analysis suggests that few regulatory components are conserved across the phylum, and those that are conserved are involved in the exploitation of the intracellular niche. Our study thus lends insight into a component of chlamydial evolution that has otherwise remained largely unexplored.

Reconstruction of novel transcription factor regulons through inference of their binding sites

Background

In most sequenced organisms the number of known regulatory genes (e.g., transcription factors (TFs)) vastly exceeds the number of experimentally-verified regulons that could be associated with them. At present, identification of TF regulons is mostly done through comparative genomics approaches. Such methods could miss organism-specific regulatory interactions and often require expensive and time-consuming experimental techniques to generate the underlying data.

Results

In this work, we present an efficient algorithm that aims to identify a given transcription factor’s regulon through inference of its unknown binding sites, based on the discovery of its binding motif. The proposed approach relies on computational methods that utilize gene expression data sets and knockout fitness data sets which are available or may be straightforwardly obtained for many organisms. We computationally constructed the profiles of putative regulons for the TFs LexA, PurR and Fur in E. coli K12 and identified their binding motifs. Comparisons with an experimentally-verified database showed high recovery rates of the known regulon members, and indicated good predictions for the newly found genes with high biological significance. The proposed approach is also applicable to novel organisms for predicting unknown regulons of the transcriptional regulators. Results for the hypothetical protein Dde0289 in D. alaskensis include the discovery of a Fis-type TF binding motif.

Conclusions

The proposed motif-based regulon inference approach can discover the organism-specific regulatory interactions on a single genome, which may be missed by current comparative genomics techniques due to their limitations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0685-y) contains supplementary material, which is available to authorized users.

Horizontal gene transfer in an acid mine drainage microbial community

Background

Horizontal gene transfer (HGT) has been widely identified in complete prokaryotic genomes. However, the roles of HGT among members of a microbial community and in evolution remain largely unknown. With the emergence of metagenomics, it is nontrivial to investigate such horizontal flow of genetic materials among members in a microbial community from the natural environment. Because of the lack of suitable methods for metagenomics gene transfer detection, microorganisms from a low-complexity community acid mine drainage (AMD) with near-complete genomes were used to detect possible gene transfer events and suggest the biological significance.

Results

Using the annotation of coding regions by the current tools, a phylogenetic approach, and an approximately unbiased test, we found that HGTs in AMD organisms are not rare, and we predicted 119 putative transferred genes. Among them, 14 HGT events were determined to be transfer events among the AMD members. Further analysis of the 14 transferred genes revealed that the HGT events affected the functional evolution of archaea or bacteria in AMD, and it probably shaped the community structure, such as the dominance of G-plasma in archaea in AMD through HGT.

Conclusions

Our study provides a novel insight into HGT events among microorganisms in natural communities. The interconnectedness between HGT and community evolution is essential to understand microbial community formation and development.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1720-0) contains supplementary material, which is available to authorized users.

An Arabidopsis Transcriptional Regulatory Map Reveals Distinct Functional and Evolutionary Features of Novel Transcription Factors

Transcription factors (TFs) play key roles in both development and stress responses. By integrating into and rewiring original systems, novel TFs contribute significantly to the evolution of transcriptional regulatory networks. Here, we report a high-confidence transcriptional regulatory map covering 388 TFs from 47 families in Arabidopsis. Systematic analysis of this map revealed the architectural heterogeneity of developmental and stress response subnetworks and identified three types of novel network motifs that are absent from unicellular organisms and essential for multicellular development. Moreover, TFs of novel families that emerged during plant landing present higher binding specificities and are preferentially wired into developmental processes and these novel network motifs. Further unveiled connection between the binding specificity and wiring preference of TFs explains the wiring preferences of novel-family TFs. These results reveal distinct functional and evolutionary features of novel TFs, suggesting a plausible mechanism for their contribution to the evolution of multicellular organisms.

The loose evolutionary relationships between transcription factors and other gene products across prokaryotes

Background

Tests for the evolutionary conservation of associations between genes coding for transcription factors (TFs) and other genes have been limited to a few model organisms due to the lack of experimental information of functional associations in other organisms. We aimed at surmounting this limitation by using the most co-occurring gene pairs as proxies for the most conserved functional interactions available for each gene in a genome. We then used genes predicted to code for TFs to compare their most conserved interactions against the most conserved interactions for the rest of the genes within each prokaryotic genome available.

Results

We plotted profiles of phylogenetic profiles, p-cubic, to compare the maximally scoring interactions of TFs against those of other genes. In most prokaryotes, genes coding for TFs showed lower co-occurrences when compared to other genes. We also show that genes coding for TFs tend to have lower Codon Adaptation Indexes compared to other genes.

Conclusions

The co-occurrence tests suggest that transcriptional regulation evolves quickly in most, if not all, prokaryotes. The Codon Adaptation Index analyses suggest quick gene exchange and rewiring of transcriptional regulation across prokaryotes.

Electronic supplementary material

The online version of this article (doi:10.1186/1756-0500-7-928) contains supplementary material, which is available to authorized users.

Longissimus dorsi transcriptome analysis of purebred and crossbred Iberian pigs differing in muscle characteristics

BACKGROUND

The two main genetic types in Iberian pig production show important phenotypic differences in growth, fattening and tissue composition since early developmental stages. The objective of this work was the evaluation of muscle transcriptome profile in piglets of both genetic types, in order to identify genes, pathways and regulatory factors responsible for their phenotypic differences. Contemporary families coming from pure Iberian pigs (IB) or from crossing with Duroc boars (DU×IB) were generated. Piglets (14 from each genetic type) were slaughtered at weaning (28 days) and longissimus dorsi was sampled for composition and gene expression studies. RNA was obtained and hybridized to Affymetrix Porcine Genechip expression arrays.

RESULTS

Loin muscle chemical composition showed significant differences between genetic types in intramuscular fat content (6.1% vs. 4.3% in IB and DUxIB animals, respectively, P = 0.009) and in saturated (P = 0.019) and monounsaturated fatty acid proportions (P = 0.044). The statistical analysis of gene expression data allowed the identification of 256 differentially expressed (DE) genes between genetic types (FDR < 0.10), 102 upregulated in IB and 154 upregulated in DU×IB. Transcript differences were validated for a subset of DE genes by qPCR. We observed alteration in biological functions related to extracellular matrix function and organization, cellular adhesion, muscle growth, lipid metabolism and proteolysis. Candidate genes with known effects on muscle growth were found among the DE genes upregulated in DU×IB. Genes related to lipid metabolism and proteolysis were found among those upregulated in IB. Regulatory factors (RF) potentially involved in the expression differences were identified by calculating the regulatory impact factors. Twenty-nine RF were found, some of them with known relationship with tissue development (MSTN, SIX4, IRX3), adipogenesis (CEBPD, PPARGC1B), or extracellular matrix processes (MAX, MXI1). Correlation among the expression of these RF and DE genes show relevant differences between genetic types.

CONCLUSION

These results provide valuable information about genetic mechanisms determining the phenotypic differences on growth and meat quality between the genetic types studied, mainly related to the development and function of the extracellular matrix and also to some metabolic processes as proteolysis and lipid metabolism. Transcription factors and regulatory mechanisms are proposed for these altered biological functions.

Cross-species gene-family fluctuations reveal the dynamics of horizontal transfers

Prokaryotes vary their protein repertoire mainly through horizontal transfer and gene loss. To elucidate the links between these processes and the cross-species gene-family statistics, we perform a large-scale data analysis of the cross-species variability of gene-family abundance (the number of members of the family found on a given genome). We find that abundance fluctuations are related to the rate of horizontal transfers. This is rationalized by a minimal theoretical model, which predicts this link. The families that are not captured by the model show abundance profiles that are markedly peaked around a mean value, possibly because of specific abundance selection. Based on these results, we define an abundance variability index that captures a family's evolutionary behavior (and thus some of its relevant functional properties) purely based on its cross-species abundance fluctuations. Analysis and model, combined, show a quantitative link between cross-species family abundance statistics and horizontal transfer dynamics, which can be used to analyze genome ‘flux’. Groups of families with different values of the abundance variability index correspond to genome sub-parts having different plasticity in terms of the level of horizontal exchange allowed by natural selection.

Building a better bacillus: the emergence of Mycobacterium tuberculosis

The genus Mycobacterium is comprised of more than 150 species that reside in a wide variety of habitats. Most mycobacteria are environmental organisms that are either not associated with disease or are opportunistic pathogens that cause non-transmissible disease in immunocompromised individuals. In contrast, a small number of species, such as the tubercle bacillus, Mycobacterium tuberculosis, are host-adapted pathogens for which there is no known environmental reservoir. In recent years, gene disruption studies using the host-adapted pathogen have uncovered a number of “virulence factors,” yet genomic data indicate that many of these elements are present in non-pathogenic mycobacteria. This suggests that much of the genetic make-up that enables virulence in the host-adapted pathogen is already present in environmental members of the genus. In addition to these generic factors, we hypothesize that molecules elaborated exclusively by professional pathogens may be particularly implicated in the ability of M. tuberculosis to infect, persist, and cause transmissible pathology in its host species, Homo sapiens. One approach to identify these molecules is to employ comparative analysis of mycobacterial genomes, to define evolutionary events such as horizontal gene transfer (HGT) that contributed M. tuberculosis-specific genetic elements. Independent studies have now revealed the presence of HGT genes in the M. tuberculosis genome and their role in the pathogenesis of disease is the subject of ongoing investigations. Here we review these studies, focusing on the hypothesized role played by HGT loci in the emergence of M. tuberculosis from a related environmental species into a highly specialized human-adapted pathogen.

Promoter islands as a platform for interaction with nucleoid proteins and transcription factors

Seventy-eight promoter islands with an extraordinarily high density of potential promoters have been recently found in the genome of Escherichia coli. It has been shown that RNA polymerase binds internal promoters of these islands and produces short oligonucleotides, while the synthesis of normal mRNAs is suppressed. This quenching may be biologically relevant, as most islands are associated with foreign genes, which expression may deplete cellular resources. However, a molecular mechanism of silencing with the participation of these promoter-rich regions remains obscure. It has been demonstrated that all islands interact with histone-like protein H-NS--a specific sentinel of foreign genes. In this study, we demonstrated the inhibitory effect of H-NS using Δhns mutant of Escherichia coli and showed that deletion of dps, encoding another protein of bacterial nucleoid, tended to decrease rather than increase the amount of island-specific transcripts. This observation precluded consideration of promoter islands as sites for targeted heterochromatization only and a computer search for the binding sites of 53 transcription factors (TFs) revealed six proteins, which may specifically regulate their transcriptional output.

Regulation of Salmonella enterica pathogenicity island 1 (SPI-1) by the LysR-type regulator LeuO

LeuO is a quiescent LysR-type regulator belonging to the H-NS regulon. Activation of leuO transcription represses expression of pathogenicity island 1 (SPI-1) in Salmonella enterica serovar Typhimurium and inhibits invasion of epithelial cells. Loss of HilE suppresses LeuO-mediated downregulation of SPI-1. Activation of leuO transcription reduces the level of HilD protein, and loss of HilE restores the wild type HilD level. Hence, LeuO-mediated downregulation of SPI-1 may involve inhibition of HilD activity by HilE, a view consistent with the fact that HilE is a HilD inhibitor. In vivo analyses using β-galactosidase fusions indicate that LeuO activates hilE transcription. In vitro analyses by slot blotting, electrophoretic mobility shift analysis and DNase I footprinting show that LeuO binds the hilE promoter region. Although residual SPI-1 repression by LeuO is observed in the absence of HilE, the LeuO-HilE-HilD 'pathway' appears to be the major mechanism. Because both leuO and SPI-1 are repressed by H-NS, activation of leuO transcription may provide a backup mechanism for SPI-1 repression under conditions that impair H-NS-mediated silencing.

Promoters of Escherichia coli versus promoter islands: function and structure comparison

Expression of bacterial genes takes place under the control of RNA polymerase with exchangeable σ-subunits and multiple transcription factors. A typical promoter region contains one or several overlapping promoters. In the latter case promoters have the same or different σ-specificity and are often subjected to different regulatory stimuli. Genes, transcribed from multiple promoters, have on average higher expression levels. However, recently in the genome of Escherichia coli we found 78 regions with an extremely large number of potential transcription start points (promoter islands, PIs). It was shown that all PIs interact with RNA polymerase in vivo and are able to form transcriptionally competent open complexes both in vitro and in vivo but their transcriptional activity measured by oligonucleotide microarrays was very low, if any. Here we confirmed transcriptional defectiveness of PIs by analyzing the 5'-end specific RNA-seq data, but showed their ability to produce short oligos (9-14 bases). This combination of functional properties indicated a deliberate suppression of transcriptional activity within PIs. According to our data this suppression may be due to a specific conformation of the DNA double helix, which provides an ideal platform for interaction with both RNA polymerase and the histone-like nucleoid protein H-NS. The genomic DNA of E.coli contains therefore several dozen sites optimized by evolution for staying in a heterochromatin-like state. Since almost all promoter islands are associated with horizontally acquired genes, we offer them as specific components of bacterial evolution involved in acquisition of foreign genetic material by turning off the expression of toxic or useless aliens or by providing optimal promoter for beneficial genes. The putative molecular mechanism underlying the appearance of promoter islands within recipient genomes is discussed.

Indirect and suboptimal control of gene expression is widespread in bacteria

Gene regulation in bacteria is usually described as an adaptive response to an environmental change so that genes are expressed when they are required. We instead propose that most genes are under indirect control: their expression responds to signal(s) that are not directly related to the genes' function. Indirect control should perform poorly in artificial conditions, and we show that gene regulation is often maladaptive in the laboratory. In Shewanella oneidensis MR-1, 24% of genes are detrimental to fitness in some conditions, and detrimental genes tend to be highly expressed instead of being repressed when not needed. In diverse bacteria, there is little correlation between when genes are important for optimal growth or fitness and when those genes are upregulated. Two common types of indirect control are constitutive expression and regulation by growth rate; these occur for genes with diverse functions and often seem to be suboptimal. Because genes that have closely related functions can have dissimilar expression patterns, regulation may be suboptimal in the wild as well as in the laboratory.

Signal correlations in ecological niches can shape the organization and evolution of bacterial gene regulatory networks

Transcriptional regulation plays a significant role in the biological response of bacteria to changing environmental conditions. Therefore, mapping transcriptional regulatory networks is an important step not only in understanding how bacteria sense and interpret their environment but also to identify the functions involved in biological responses to specific conditions. Recent experimental and computational developments have facilitated the characterization of regulatory networks on a genome-wide scale in model organisms. In addition, the multiplication of complete genome sequences has encouraged comparative analyses to detect conserved regulatory elements and infer regulatory networks in other less well-studied organisms. However, transcription regulation appears to evolve rapidly, thus, creating challenges for the transfer of knowledge to nonmodel organisms. Nevertheless, the mechanisms and constraints driving the evolution of regulatory networks have been the subjects of numerous analyses, and several models have been proposed. Overall, the contributions of mutations, recombination, and horizontal gene transfer are complex. Finally, the rapid evolution of regulatory networks plays a significant role in the remarkable capacity of bacteria to adapt to new or changing environments. Conversely, the characteristics of environmental niches determine the selective pressures and can shape the structure of regulatory network accordingly.

Expression divergence between Escherichia coli and Salmonella enterica serovar Typhimurium reflects their lifestyles

Escherichia coli K12 is a commensal bacteria and one of the best-studied model organisms. Salmonella enterica serovar Typhimurium, on the other hand, is a facultative intracellular pathogen. These two prokaryotic species can be considered related phylogenetically, and they share a large amount of their genetic material, which is commonly termed the "core genome." Despite their shared core genome, both species display very different lifestyles, and it is unclear to what extent the core genome, apart from the species-specific genes, plays a role in this lifestyle divergence. In this study, we focus on the differences in expression domains for the orthologous genes in E. coli and S. Typhimurium. The iterative comparison of coexpression methodology was used on large expression compendia of both species to uncover the conservation and divergence of gene expression. We found that gene expression conservation occurs mostly independently from amino acid similarity. According to our estimates, at least more than one quarter of the orthologous genes has a different expression domain in E. coli than in S. Typhimurium. Genes involved with key cellular processes are most likely to have conserved their expression domains, whereas genes showing diverged expression are associated with metabolic processes that, although present in both species, are regulated differently. The expression domains of the shared "core" genome of E. coli and S. Typhimurium, consisting of highly conserved orthologs, have been tuned to help accommodate the differences in lifestyle and the pathogenic potential of Salmonella.

Transcription factor family-based reconstruction of singleton regulons and study of the Crp/Fnr, ArsR, and GntR families in Desulfovibrionales genomes

Accurate detection of transcriptional regulatory elements is essential for high-quality genome annotation, metabolic reconstruction, and modeling of regulatory networks. We developed a computational approach for reconstruction of regulons operated by transcription factors (TFs) from large protein families and applied this novel approach to three TF families in 10 Desulfovibrionales genomes. Phylogenetic analyses of 125 regulators from the ArsR, Crp/Fnr, and GntR families revealed that 65% of these regulators (termed reference TFs) are well conserved in Desulfovibrionales, while the remaining 35% of regulators (termed singleton TFs) are species specific and show a mosaic distribution. For regulon reconstruction in the group of singleton TFs, the standard orthology-based approach was inefficient, and thus, we developed a novel approach based on the simultaneous study of all homologous TFs from the same family in a group of genomes. As a result, we identified binding for 21 singleton TFs and for all reference TFs in all three analyzed families. Within each TF family we observed structural similarities between DNA-binding motifs of different reference and singleton TFs. The collection of reconstructed regulons is available at the RegPrecise database (http://regprecise.lbl.gov/RegPrecise/Desulfovibrionales.jsp).