Links between Transcription, Environmental Adaptation and Gene Variability in Escherichia coli: Correlations between Gene Expression and Gene Variability Reflect Growth Efficiencies

Can Transcriptomics Elucidate the Role of Regulation in Invasion Success?

When a species invades a novel environment, it must bridge the environment-phenotype mismatch in its new range to persist. Contemporary invasion biology research has focused on the role that trait variation and adaptation, and their underlying genomic factors, play in a species' adaptive potential, and thus facilitating invasion. Empirical studies have provided valuable insights into phenotypes that persist and arise in novel environments, coupled with 'omics tools that further the understanding of the contributions of genomic architecture in species establishment. Particularly, the use of transcriptomics to explore the role of plasticity in the initial stages of an invasion is growing. Here, we assess the role of various mechanisms relating to regulation and functional adaptation (often measured via the transcriptome) that support trait-specific plasticity in invasive species, allowing phenotypic variability without directly altering genomic diversity. First, we present a comprehensive review of the studies utilising transcriptomics in invasion biology. Second, we collate the evidence for and against the role of a range of regulatory processes in contributing to invasive species plasticity. Finally, we pose open questions in invasion biology where the use of transcriptome data may be valuable, as well as discuss the methodological limitations.

Functional Optimization in Distinct Tissues and Conditions Constrains the Rate of Protein Evolution

Understanding the main determinants of protein evolution is a fundamental challenge in biology. Despite many decades of active research, the molecular and cellular mechanisms underlying the substantial variability of evolutionary rates across cellular proteins are not currently well understood. It also remains unclear how protein molecular function is optimized in the context of multicellular species and why many proteins, such as enzymes, are only moderately efficient on average. Our analysis of genomics and functional datasets reveals in multiple organisms a strong inverse relationship between the optimality of protein molecular function and the rate of protein evolution. Furthermore, we find that highly expressed proteins tend to be substantially more functionally optimized. These results suggest that cellular expression costs lead to more pronounced functional optimization of abundant proteins and that the purifying selection to maintain high levels of functional optimality significantly slows protein evolution. We observe that in multicellular species both the rate of protein evolution and the degree of protein functional efficiency are primarily affected by expression in several distinct cell types and tissues, specifically, in developed neurons with upregulated synaptic processes in animals and in young and fast-growing tissues in plants. Overall, our analysis reveals how various constraints from the molecular, cellular, and species' levels of biological organization jointly affect the rate of protein evolution and the level of protein functional adaptation.

Cyclical heat stress during lactation influences the microstructure of the bovine mammary gland

This study aimed to evaluate the impact of heat stress on mammary epithelial cell (MEC) losses into milk, secretory mammary tissue structure, and mammary epithelial cell activity. Sixteen multiparous Holstein cows (632 ± 12 kg BW) approximately 100 d in milk housed in climate-controlled rooms were paired by body weight and randomly allocated to one of 2 treatments, heat stress (HS) or pair feeding thermoneutral (PFTN) using 2 cohorts. Each cohort was subjected to 2 periods of 4 d each. In period 1, both treatments had ad libitum access to a common total mixed ration and were exposed to a controlled daily temperature-humidity index (THI) of 64. In period 2, HS cows were exposed to controlled cyclical heat stress (THI: 74 to 80), while PFTN cows remained at 64 THI and daily dry matter intake was matched to HS. Cows were milked twice daily, and milk yield was recorded at each milking. Individual milk samples on the last day of each period were used to quantify MEC losses by flow cytometry using butyrophilin as a cell surface marker. On the final day of period 2, individual bovine mammary tissue samples were obtained for histomorphology analysis, assessment of protein abundance, and evaluation of gene expression of targets associated with cellular capacity for milk and milk component synthesis, heat response, cellular proliferation, and autophagy. Statistical analysis was performed using the GLIMMIX procedure of SAS. Milk yield was reduced by 4.3 kg by HS (n = 7) compared with PFTN (n = 8). Independent of treatment, MEC in milk averaged 174 cells/mL (2.9% of total cells). There was no difference between HS vs. PFTN cows for MEC shed or concentration in milk. Alveolar area was reduced 25% by HS, and HS had 4.1 more alveoli than PFTN. Total number of nucleated MEC per area were greater in HS (389 ± 1.05) compared with PFTN (321 ± 1.05); however, cell number per alveolus was similar between groups (25 ± 1.5 vs. 26 ± 1.4). There were no differences in relative fold expression for GLUT1, GLUT8, CSN2, CSN3, LALBA, FASN, HSPA5, and HSPA8 in HS compared with PFTN. Immunoblotting analyses showed a decrease abundance for phosphorylated STAT5 and S6K1, and an increase in LC3 II in HS compared with PFTN. These results suggest that even if milk yield differences and histological changes occur in the bovine mammary gland after 4 d of heat exposure, MEC loss into milk, nucleated MEC number per alveolus, and gene expression of nutrient transport, milk component synthesis, and heat stress related targets are unaffected. In contrast, the abundance of proteins related to protein synthesis and cell survival decreased significantly, while an upregulation of proteins associated with autophagy in HS compared with PFTN.

Regulation of bacterial gene expression by non-coding RNA: It is all about time!

Commensal and pathogenic bacteria continuously evolve to survive in diverse ecological niches by efficiently coordinating gene expression levels in their ever-changing environments. Regulation through the RNA transcript itself offers a faster and more cost-effective way to adapt than protein-based mechanisms and can be leveraged for diagnostic or antimicrobial purposes. However, RNA can fold into numerous intricate, not always functional structures that both expand and obscure the plethora of roles that regulatory RNAs serve within the cell. Here, we review the current knowledge of bacterial non-coding RNAs in relation to their folding pathways and interactions. We posit that co-transcriptional folding of these transcripts ultimately dictates their downstream functions. Elucidating the spatiotemporal folding of non-coding RNAs during transcription therefore provides invaluable insights into bacterial pathogeneses and predictive disease diagnostics. Finally, we discuss the implications of co-transcriptional folding andapplications of RNAs for therapeutics and drug targets.

Machine learning-assisted medium optimization revealed the discriminated strategies for improved production of the foreign and native metabolites

The composition of medium components is crucial for achieving the best performance of synthetic construction in genetically engineered cells. Which and how medium components determine the performance, e.g., productivity, remain poorly investigated. To address the questions, a comparative survey with two genetically engineered Escherichia coli strains was performed. As a case study, the strains carried the synthetic pathways for producing the aromatic compounds of 4-aminophenylalanine (4APhe) or tyrosine (Tyr), common in the upstream but differentiated in the downstream metabolism. Bacterial growth and compound production were examined in hundreds of medium combinations that comprised 48 pure chemicals. The resultant data sets linking the medium composition to bacterial growth and production were subjected to machine learning for improved production. Intriguingly, the primary medium components determining the production of 4PheA and Tyr were differentiated, which were the initial resource (glucose) of the synthetic pathway and the inducer (IPTG) of the synthetic construction, respectively. Fine-tuning of the primary component significantly increased the yields of 4APhe and Tyr, indicating that a single component could be crucial for the performance of synthetic construction. Transcriptome analysis observed the local and global changes in gene expression for improved production of 4APhe and Tyr, respectively, revealing divergent metabolic strategies for producing the foreign and native metabolites. The study demonstrated that ML-assisted medium optimization could provide a novel point of view on how to make the synthetic construction meet the designed working principle and achieve the expected biological function.

Growth rate-associated transcriptome reorganization in response to genomic, environmental, and evolutionary interruptions

The genomic, environmental, and evolutionary interruptions caused the changes in bacterial growth, which were stringently associated with changes in gene expression. The growth and gene expression changes remained unclear in response to these interruptions that occurred combinative. As a pilot study, whether and how bacterial growth was affected by the individual and dual interruptions of genome reduction, environmental stress, and adaptive evolution were investigated. Growth assay showed that the presence of the environmental stressors, i.e., threonine and chloramphenicol, significantly decreased the growth rate of the wild-type Escherichia coli, whereas not that of the reduced genome. It indicated a canceling effect in bacterial growth due to the dual interruption of the genomic and environmental changes. Experimental evolution of the reduced genome released the canceling effect by improving growth fitness. Intriguingly, the transcriptome architecture maintained a homeostatic chromosomal periodicity regardless of the genomic, environmental, and evolutionary interruptions. Negative epistasis in transcriptome reorganization was commonly observed in response to the dual interruptions, which might contribute to the canceling effect. It was supported by the changes in the numbers of differentially expressed genes (DEGs) and the enriched regulons and functions. Gene network analysis newly constructed 11 gene modules, one out of which was correlated to the growth rate. Enrichment of DEGs in these modules successfully categorized them into three types, i.e., conserved, responsive, and epistatic. Taken together, homeostasis in transcriptome architecture was essential to being alive, and it might be attributed to the negative epistasis in transcriptome reorganization and the functional differentiation in gene modules. The present study directly connected bacterial growth fitness with transcriptome reorganization and provided a global view of how microorganisms responded to genomic, environmental, and evolutionary interruptions for survival from wild nature.

Variation in heat shock protein 40 kDa relates to divergence in thermotolerance among cryptic rotifer species

Genetic divergence and the frequency of hybridization are central for defining species delimitations, especially among cryptic species where morphological differences are merely absent. Rotifers are known for their high cryptic diversity and therefore are ideal model organisms to investigate such patterns. Here, we used the recently resolved Brachionus calyciflorus species complex to investigate whether previously observed between species differences in thermotolerance and gene expression are also reflected in their genomic footprint. We identified a Heat Shock Protein gene (HSP 40 kDa) which exhibits cross species pronounced sequence variation. This gene exhibits species-specific fixed sites, alleles, and sites putatively under positive selection. These sites are located in protein binding regions involved in chaperoning and may therefore reflect adaptive diversification. By comparing three genetic markers (ITS, COI, HSP 40 kDa), we revealed hybridization events between the cryptic species. The low frequency of introgressive haplotypes/alleles suggest a tight, but not fully impermeable boundary between the cryptic species.

Purifying selection enduringly acts on the sequence evolution of highly expressed proteins in Escherichia coli

The evolutionary speed of a protein sequence is constrained by its expression level, with highly expressed proteins evolving relatively slowly. This negative correlation between expression levels and evolutionary rates (known as the E–R anticorrelation) has already been widely observed in past macroevolution between species from bacteria to animals. However, it remains unclear whether this seemingly general law also governs recent evolution, including past and de novo, within a species. However, the advent of genomic sequencing and high-throughput phenotyping, particularly for bacteria, has revealed fundamental gaps between the 2 evolutionary processes and has provided empirical data opposing the possible underlying mechanisms which are widely believed. These conflicts raise questions about the generalization of the E–R anticorrelation and the relevance of plausible mechanisms. To explore the ubiquitous impact of expression levels on molecular evolution and test the relevance of the possible underlying mechanisms, we analyzed the genome sequences of 99 strains of Escherichia coli for evolution within species in nature. We also analyzed genomic mutations accumulated under laboratory conditions as a model of de novo evolution within species. Here, we show that E–R anticorrelation is significant in both past and de novo evolution within species in E. coli. Our data also confirmed ongoing purifying selection on highly expressed genes. Ongoing selection included codon-level purifying selection, supporting the relevance of the underlying mechanisms. However, the impact of codon-level purifying selection on the constraints in evolution within species might be smaller than previously expected from evolution between species.

Escherichia coli transcription factors of unknown function: sequence features and possible evolutionary relationships

Organisms need mechanisms to perceive the environment and respond accordingly to environmental changes or the presence of hazards. Transcription factors (TFs) are required for cells to respond to the environment by controlling the expression of genes needed. Escherichia coli has been the model bacterium for many decades, and still, there are features embedded in its genome that remain unstudied. To date, 58 TFs remain poorly characterized, although their binding sites have been experimentally determined. This study showed that these TFs have sequence variation at the third codon position G+C content but maintain the same Codon Adaptation Index (CAI) trend as annotated functional transcription factors. Most of these transcription factors are in areas of the genome where abundant repetitive and mobile elements are present. Sequence divergence points to groups with distinctive sequence signatures but maintaining the same type of DNA binding domain. Finally, the analysis of the promoter sequences of the 58 TFs showed A+T rich regions that agree with the features of horizontally transferred genes. The findings reported here pave the way for future research of these TFs that may uncover their role as spare factors in case of lose-of-function mutations in core TFs and trace back their evolutionary history.

Transcriptomic profiling of Escherichia coli K-12 in response to a compendium of stressors

Environmental perturbations impact multiple cellular traits, including gene expression. Bacteria respond to these stressful situations through complex gene interaction networks, thereby inducing stress tolerance and survival of cells. In this paper, we study the response mechanisms of E. coli when exposed to different environmental stressors via differential expression and co-expression analysis. Gene co-expression networks were generated and analyzed via Weighted Gene Co-expression Network Analysis (WGCNA). Based on the gene co-expression networks, genes with similar expression profiles were clustered into modules. The modules were analysed for identification of hub genes, enrichment of biological processes and transcription factors. In addition, we also studied the link between transcription factors and their differentially regulated targets to understand the regulatory mechanisms involved. These networks validate known gene interactions and provide new insights into genes mediating transcriptional regulation in specific stress environments, thus allowing for in silico hypothesis generation.

The gut environment regulates bacterial gene expression which modulates susceptibility to bacteriophage infection

Abundance and diversity of bacteria and their viral predators, bacteriophages (phages), in the digestive tract are associated with human health. Particularly intriguing is the long-term coexistence of these two antagonistic populations. We performed genome-wide RNA sequencing on a human enteroaggregative Escherichia coli isolate to identify genes differentially expressed between in vitro conditions and in murine intestines. We experimentally demonstrated that four of these differentially expressed genes modified the interactions between E. coli and three virulent phages by either increasing or decreasing its susceptibility/resistance pattern and also by interfering with biofilm formation. Therefore, the regulation of bacterial genes expression during the colonization of the digestive tract influences the coexistence of phages and bacteria, highlighting the intricacy of tripartite relationships between phages, bacteria, and the animal host in intestinal homeostasis.

Bridging Tumorigenesis and Therapy Resistance With a Non-Darwinian and Non-Lamarckian Mechanism of Adaptive Evolution

Although neo-Darwinian (and less often Lamarckian) dynamics are regularly invoked to interpret cancer's multifarious molecular profiles, they shine little light on how tumorigenesis unfolds and often fail to fully capture the frequency and breadth of resistance mechanisms. This uncertainty frames one of the most problematic gaps between science and practice in modern times. Here, we offer a theory of adaptive cancer evolution, which builds on a molecular mechanism that lies outside neo-Darwinian and Lamarckian schemes. This mechanism coherently integrates non-genetic and genetic changes, ecological and evolutionary time scales, and shifts the spotlight away from positive selection towards purifying selection, genetic drift, and the creative-disruptive power of environmental change. The surprisingly simple use-it or lose-it rationale of the proposed theory can help predict molecular dynamics during tumorigenesis. It also provides simple rules of thumb that should help improve therapeutic approaches in cancer.

Escherichia coli: Physiological Clues Which Turn On the Synthesis of Antimicrobial Molecules

Zoonotic pathogens, like Shiga toxin-producing Escherichia coli (STEC) are a food safety and health risk. To battle the increasing emergence of virulent microbes, novel mitigation strategies are needed. One strategy being considered to combat pathogens is antimicrobial compounds produced by microbes, coined microcins. However, effectors for microcin production are poorly understood, particularly in the context of complex physiological responses along the gastro-intestinal tract (GIT). Previously, we identified an E. coli competitor capable of producing a strong diffusible antimicrobial with microcin-associated characteristics. Our objective was to examine how molecule production of this competitor is affected by physiological properties associated with the GIT, namely the effects of carbon source, bile salt concentration and growth phase. Using previously described liquid- and agar-based assays determined that carbon sources do not affect antimicrobial production of E. coli O103F. However, bile salt concentrations affected production significantly, suggesting that E. coli O103F uses cues along the GIT to modulate the expression of antimicrobial production. Furthermore, E. coli O103F produces the molecule during the exponential phase, contrary to most microcins identified to date. The results underscored the importance of experimental design to identify producers of antimicrobials. To detect antimicrobials, conventional microbiological methods can be a starting point, but not the gold standard.

Correlated chromosomal periodicities according to the growth rate and gene expression

Linking genetic information to population fitness is crucial to understanding living organisms. Despite the abundant knowledge of the genetic contribution to growth, the overall patterns/features connecting genes, their expression, and growth remain unclear. To reveal the quantitative and direct connections, systematic growth assays of single-gene knockout Escherichia coli strains under both rich and poor nutritional conditions were performed; subsequently, the resultant growth rates were associated with the original expression levels of the knockout genes in the parental genome. Comparative analysis of growth and the transcriptome identified not only the nutritionally differentiated fitness cost genes but also a significant correlation between the growth rates of the single-gene knockout strains and the original expression levels of these knockout genes in the parental strain, regardless of the nutritional variation. In addition, the coordinated chromosomal periodicities of the wild-type transcriptome and the growth rates of the strains lacking the corresponding genes were observed. The common six-period periodicity was somehow attributed to the essential genes, although the underlying mechanism remains to be addressed. The correlated chromosomal periodicities associated with the gene expression-growth dataset were highly valuable for bacterial growth prediction and discovering the working principles governing minimal genetic information.

Evolution and Plasticity of the Transcriptome Under Temperature Fluctuations in the Fungal Plant Pathogen Zymoseptoria tritici

Most species live in a variable environment in nature. Yet understanding the evolutionary processes underlying molecular adaptation to fluctuations remains a challenge. In this study we investigate the transcriptome of the fungal wheat pathogen Zymoseptoria tritici after experimental evolution under stable or fluctuating temperature, by comparing ancestral and evolved populations simultaneously. We found that temperature regimes could have a large and pervasive effect on the transcriptome evolution, with as much as 38% of the genes being differentially expressed between selection regimes. Although evolved lineages showed different changes of gene expression based on ancestral genotypes, we identified a set of genes responding specifically to fluctuation. We found that transcriptome evolution in fluctuating conditions was repeatable between parallel lineages initiated from the same genotype for about 60% of the differentially expressed genes. Further, we detected several hotspots of significantly differentially expressed genes in the genome, in regions known to be enriched in repetitive elements, including accessory chromosomes. Our findings also evidenced gene expression evolution toward a gain of robustness (loss of phenotypic plasticity) associated with the fluctuating regime, suggesting robustness is adaptive in changing environment. This work provides valuable insight into the role of transcriptional rewiring for rapid adaptation to abiotic changes in filamentous plant pathogens.

The highly conserved chromosomal periodicity of transcriptomes and the correlation of its amplitude with the growth rate in Escherichia coli

The growth rate, representing the fitness of a bacterial population, is determined by the transcriptome. Chromosomal periodicity, which is known as the periodic spatial pattern of a preferred chromosomal distance in microbial genomes, is a representative overall feature of the transcriptome; however, whether and how it is associated with the bacterial growth rate are unknown. To address these questions, we analysed a total of 213 transcriptomes of multiple Escherichia coli strains growing in an assortment of culture conditions varying in terms of temperature, nutrition level and osmotic pressure. Intriguingly, Fourier transform analyses of the transcriptome identified a common chromosomal periodicity of transcriptomes, which was independent of the variation in genomes and environments. In addition, fitting of the data to a theoretical model, we found that the amplitudes of the periodic transcriptomes were significantly correlated with the growth rates. These results indicated that the amplitude of periodic transcriptomes is a parameter representing the global pattern of gene expression in correlation with the bacterial growth rate. Thus, our study provides a novel parameter for evaluating the adaptiveness of a growing bacterial population and quantitatively predicting the growth dynamics according to the global expression pattern.

Hypermutation as an Evolutionary Mechanism for Achromobacter xylosoxidans in Cystic Fibrosis Lung Infection

Achromobacter xylosoxidans can cause chronic infections in the lungs of patients with cystic fibrosis (CF) by adapting to the specific environment. The study of longitudinal isolates allows to investigate its within-host evolution to unravel the adaptive mechanisms contributing to successful colonization. In this study, four clinical isolates longitudinally collected from two chronically infected patients underwent whole genome sequencing, de novo assembly and sequence analysis. Phenotypic assays were also performed. The isolates coming from one of the patients (patient A) presented a greater number of genetic variants, diverse integrative and conjugative elements, and different protease secretion. In the first of these isolates (strain A1), we also found a large deletion in the mutS gene, involved in DNA mismatch repair (MMR). In contrast, isolates from patient B showed a lower number of variants, only one integrative and mobilizable element, no phenotypic changes, and no mutations in the MMR system. These results suggest that in the two patients the establishment of a chronic infection was mediated by different adaptive mechanisms. While the strains isolated from patient B showed a longitudinal microevolution, strain A1 can be clearly classified as a hypermutator, confirming the occurrence and importance of this adaptive mechanism in A. xylosoxidans infection.

Mutation Rate Heterogeneity Increases Odds of Survival in Unpredictable Environments

Mutation rates affect both a population's present fitness and its capacity to adapt to future environmental changes. When the available genetic variability limits adaptation to environmental change, natural selection favors high mutations rates. However, constitutively high mutation rates compromise the fitness of a population in stable environments. This problem may be resolved if an increase in mutation rates is limited to times of stress, restricted to some genomic regions, and occurs only in a subpopulation of cells. Such within-population heterogeneity of mutation rates can result from genetic, environmental, and stochastic effects. The presence of subpopulations of transient mutator cells does not jeopardize the overall fitness of a population under stable environmental conditions. However, they can increase the odds of survival in changing environments because they represent reservoirs of increased genetic variability. This article presents evidence that such heterogeneity of mutation rates is more the norm than the exception.

The Diversity of Bacterial Lifestyles Hampers Bacteriophage Tenacity

Phage therapy is based on a simple concept: the use of a virus (bacteriophage) that is capable of killing specific pathogenic bacteria to treat bacterial infections. Since the pioneering work of Félix d’Herelle, bacteriophages (phages) isolated in vitro have been shown to be of therapeutic value. Over decades of study, a large number of rather complex mechanisms that are used by phages to hijack bacterial resources and to produce their progeny have been deciphered. While these mechanisms have been identified and have been studied under optimal conditions in vitro, much less is known about the requirements for successful viral infections in relevant natural conditions. This is particularly true in the context of phage therapy. Here, we highlight the parameters affecting phage replication in both in vitro and in vivo environments, focusing, in particular, on the mammalian digestive tract. We propose avenues for increasing the knowledge-guided implementation of phages as therapeutic tools.

Quorum Sensing versus Quenching Bacterial Isolates Obtained from MBR Plants Treating Leachates from Municipal Solid Waste

Quorum sensing (QS) is a mechanism dependent on bacterial density. This coordinated process is mediated by the synthesis and the secretion of signal molecules, called autoinducers (AIs). N-acyl-homoserine lactones (AHLs) are the most common AIs that are used by Gram-negative bacteria and are involved in biofilm formation. Quorum Quenching (QQ) is the interference of QS by producing hydrolyzing enzymes, among other strategies. The main objective of the present study was to identify QS and QQ strains from MBR wastewater treatment plants. A total of 99 strains were isolated from two Spanish plants that were intended to treat leachate from municipal solid waste. Five AHL producers were detected using AHL biosensor strains (Chromobacterium violaceum CV026 and Agrobacterium tumefaciens NT1). Fifteen strains of seventy-one Gram-positive were capable of eliminating or reducing at least one AHL activity. The analysis of 16S rRNA gene sequence showed the importance of the Pseudomonas genus in the production of biofilms and the relevance of the genus Bacillus in the disruption of the QS mechanism, in which the potential activity of lactonase or acylase enzymes was investigated with the aim to contribute to solve biofouling problems and to increase the useful lifespan of membranes.

Coordinated regulation of core and accessory genes in the multipartite genome of Sinorhizobium fredii

Prokaryotes benefit from having accessory genes, but it is unclear how accessory genes can be linked with the core regulatory network when developing adaptations to new niches. Here we determined hierarchical core/accessory subsets in the multipartite pangenome (composed of genes from the chromosome, chromid and plasmids) of the soybean microsymbiont Sinorhizobium fredii by comparing twelve Sinorhizobium genomes. Transcriptomes of two S. fredii strains at mid-log and stationary growth phases and in symbiotic conditions were obtained. The average level of gene expression, variation of expression between different conditions, and gene connectivity within the co-expression network were positively correlated with the gene conservation level from strain-specific accessory genes to genus core. Condition-dependent transcriptomes exhibited adaptive transcriptional changes in pangenome subsets shared by the two strains, while strain-dependent transcriptomes were enriched with accessory genes on the chromid. Proportionally more chromid genes than plasmid genes were co-expressed with chromosomal genes, while plasmid genes had a higher within-replicon connectivity in expression than chromid ones. However, key nitrogen fixation genes on the symbiosis plasmid were characterized by high connectivity in both within- and between-replicon analyses. Among those genes with host-specific upregulation patterns, chromosomal znu and mdt operons, encoding a conserved high-affinity zinc transporter and an accessory multi-drug efflux system, respectively, were experimentally demonstrated to be involved in host-specific symbiotic adaptation. These findings highlight the importance of integrative regulation of hierarchical core/accessory components in the multipartite genome of bacteria during niche adaptation and in shaping the prokaryotic pangenome in the long run.

High Prevalence of Quorum-Sensing and Quorum-Quenching Activity among Cultivable Bacteria and Metagenomic Sequences in the Mediterranean Sea

There is increasing evidence being accumulated regarding the importance of N-acyl homoserine lactones (AHL)-mediated quorum-sensing (QS) and quorum-quenching (QQ) processes in the marine environment, but in most cases, data has been obtained from specific microhabitats, and subsequently little is known regarding these activities in free-living marine bacteria. The QS and QQ activities among 605 bacterial isolates obtained at 90 and 2000 m depths in the Mediterranean Sea were analyzed. Additionally, putative QS and QQ sequences were searched in metagenomic data obtained at different depths (15-2000 m) at the same sampling site. The number of AHL producers was higher in the 90 m sample (37.66%) than in the 2000 m sample (4.01%). However, the presence of QQ enzymatic activity was 1.63-fold higher in the 2000 m sample. The analysis of putative QQ enzymes in the metagenomes supports the relevance of QQ processes in the deepest samples, found in cultivable bacteria. Despite the unavoidable biases in the cultivation methods and biosensor assays and the possible promiscuous activity of the QQ enzymes retrieved in the metagenomic analysis, the results indicate that AHL-related QS and QQ processes could be common activity in the marine environment.

Mutator genomes decay, despite sustained fitness gains, in a long-term experiment with bacteria

Understanding the extreme variation among bacterial genomes remains an unsolved challenge in evolutionary biology, despite long-standing debate about the relative importance of natural selection, mutation, and random drift. A potentially important confounding factor is the variation in mutation rates between lineages and over evolutionary history, which has been documented in several species. Mutation accumulation experiments have shown that hypermutability can erode genomes over short timescales. These results, however, were obtained under conditions of extremely weak selection, casting doubt on their general relevance. Here, we circumvent this limitation by analyzing genomes from mutator populations that arose during a long-term experiment with Escherichia coli, in which populations have been adaptively evolving for >50,000 generations. We develop an analytical framework to quantify the relative contributions of mutation and selection in shaping genomic characteristics, and we validate it using genomes evolved under regimes of high mutation rates with weak selection (mutation accumulation experiments) and low mutation rates with strong selection (natural isolates). Our results show that, despite sustained adaptive evolution in the long-term experiment, the signature of selection is much weaker than that of mutational biases in mutator genomes. This finding suggests that relatively brief periods of hypermutability can play an outsized role in shaping extant bacterial genomes. Overall, these results highlight the importance of genomic draft, in which strong linkage limits the ability of selection to purge deleterious mutations. These insights are also relevant to other biological systems evolving under strong linkage and high mutation rates, including viruses and cancer cells.

Transcriptome Remodeling of Acinetobacter baumannii during Infection and Treatment

Acinetobacter baumannii is an increasingly common multidrug-resistant pathogen in health care settings. Although the genetic basis of antibiotic resistance mechanisms has been extensively studied, much less is known about how genetic variation contributes to other aspects of successful infections. Genetic changes that occur during host infection and treatment have the potential to remodel gene expression patterns related to resistance and pathogenesis. Longitudinal sets of multidrug-resistant A. baumannii isolates from eight patients were analyzed by RNA sequencing (RNA-seq) to identify differentially expressed genes and link them to genetic changes contributing to transcriptional variation at both within-patient and population levels. The number of differentially expressed genes among isolates from the same patient ranged from 26 (patient 588) to 145 (patient 475). Multiple patients had isolates with differential gene expression patterns related to mutations in the pmrAB and adeRS two-component regulatory system genes, as well as significant differences in genes related to antibiotic resistance, iron acquisition, amino acid metabolism, and surface-associated proteins. Population level analysis revealed 39 genetic regions with clade-specific differentially expressed genes, for which 19, 8, and 3 of these could be explained by insertion sequence mobilization, recombination-driven sequence variation, and intergenic mutations, respectively. Multiple types of mutations that arise during infection can significantly remodel the expression of genes that are known to be important in pathogenesis.

Health care-associated multidrug-resistant Acinetobacter baumannii can cause persistent infections in patients, but bacterial cells must overcome host defenses and antibiotic therapies to do so. Genetic variation arises during host infection, and new mutations are often enriched in genes encoding transcriptional regulators, iron acquisition systems, and surface-associated structures. In this study, genetic variation was shown to result in transcriptome remodeling at the level of individual patients and across phylogenetic groups. Differentially expressed genes include those related to capsule modification, iron acquisition, type I pili, and antibiotic resistance. Population level transcriptional variation reflects genome dynamics over longer evolutionary time periods, and convergent transcriptional changes support the adaptive significance of these regions. Transcriptional changes can be attributed to multiple types of genomic change, but insertion sequence mobilization had a predominant effect. The transcriptional effects of mutations that arise during infection highlight the rapid adaptation of A. baumannii during host exposure.

Elucidation of the sequential transcriptional activity in Escherichia coli using time-series RNA-seq data

Functional genomics and gene regulation inference has readily expanded our knowledge and understanding of gene interactions with regards to expression regulation. With the advancement of transcriptome sequencing in time-series comes the ability to study the sequential changes of the transcriptome. Here, we present a new method to augment regulation networks accumulated in literature with transcriptome data gathered from time-series experiments to construct a sequential representation of transcription factor activity. We apply our method on a time-series RNA-Seq data set of Escherichia coli as it transitions from growth to stationary phase over five hours and investigate the various activity in gene regulation process by taking advantage of the correlation between regulatory gene pairs to examine their activity on a dynamic network. We analyse the changes in metabolic activity of the pagP gene and associated transcription factors during phase transition, and visualize the sequential transcriptional activity to describe the change in metabolic pathway activity originating from the pagP transcription factor, phoP. We observe a shift from amino acid and nucleic acid metabolism, to energy metabolism during the transition to stationary phase in E. coli.