Applications of genome-scale metabolic models to investigate microbial metabolic adaptations in response to genetic or environmental perturbations.

Environmental perturbations are encountered by microorganisms regularly and will require metabolic adaptations to ensure an organism can survive in the newly presenting conditions. In order to study the mechanisms of metabolic adaptation in such conditions, various experimental and computational approaches have been used. Genome-scale metabolic models (GEMs) are one of the most powerful approaches to study metabolism, providing a platform to study the systems level adaptations of an organism to different environments which could otherwise be infeasible experimentally. In this review, we are describing the application of GEMs in understanding how microbes reprogram their metabolic system as a result of environmental variation. In particular, we provide the details of metabolic model reconstruction approaches, various algorithms and tools for model simulation, consequences of genetic perturbations, integration of '-omics' datasets for creating context-specific models and their application in studying metabolic adaptation due to the change in environmental conditions.

Identification of an RNA sponge that controls the RoxS riboregulator of central metabolism in Bacillus subtilis.

sRNAs are a taxonomically-restricted but transcriptomically-abundant class of post-transcriptional regulators. While of major importance for adaption to the environment, we currently lack global-scale methodology enabling target identification, especially in species without known RNA hub proteins (e.g. Hfq). Using psoralen RNA cross-linking and Illumina-sequencing we identify RNA-RNA interacting pairs in vivo in Bacillus subtilis, resolving previously well-described interactants. Although sRNA-sRNA pairings are rare (compared with sRNA-mRNA), we identify a robust example involving the conserved sRNA RoxS and an unstudied sRNA RosA (Regulator of sRNA A). We show RosA to be the first confirmed RNA sponge described in a Gram-positive bacterium. RosA interacts with at least two sRNAs, RoxS and FsrA. The RosA/RoxS interaction not only affects the levels of RoxS but also its processing and regulatory activity. We also found that the transcription of RosA is repressed by CcpA, the key regulator of carbon-metabolism in B. subtilis. Since RoxS is already known to be transcriptionally controlled by malate via the transcriptional repressor Rex, its post-transcriptional regulation by CcpA via RosA places RoxS in a key position to control central metabolism in response to varying carbon sources.

A culture-independent sequence-based metagenomics approach to the investigation of an outbreak of Shiga-toxigenic Escherichia coli O104:H4.

IMPORTANCE: Identification of the bacterium responsible for an outbreak can aid in disease management. However, traditional culture-based diagnosis can be difficult, particularly if no specific diagnostic test is available for an outbreak strain. OBJECTIVE: To explore the potential of metagenomics, which is the direct sequencing of DNA extracted from microbiologically complex samples, as an open-ended clinical discovery platform capable of identifying and characterizing bacterial strains from an outbreak without laboratory culture. DESIGN, SETTING, AND PATIENTS: In a retrospective investigation, 45 samples were selected from fecal specimens obtained from patients with diarrhea during the 2011 outbreak of Shiga-toxigenic Escherichia coli (STEC) O104:H4 in Germany. Samples were subjected to high-throughput sequencing (August-September 2012), followed by a 3-phase analysis (November 2012-February 2013). In phase 1, a de novo assembly approach was developed to obtain a draft genome of the outbreak strain. In phase 2, the depth of coverage of the outbreak strain genome was determined in each sample. In phase 3, sequences from each sample were compared with sequences from known bacteria to identify pathogens other than the outbreak strain. MAIN OUTCOMES AND MEASURES: The recovery of genome sequence data for the purposes of identification and characterization of the outbreak strain and other pathogens from fecal samples. RESULTS: During phase 1, a draft genome of the STEC outbreak strain was obtained. During phase 2, the outbreak strain genome was recovered from 10 samples at greater than 10-fold coverage and from 26 samples at greater than 1-fold coverage. Sequences from the Shiga-toxin genes were detected in 27 of 40 STEC-positive samples (67%). In phase 3, sequences from Clostridium difficile, Campylobacter jejuni, Campylobacter concisus, and Salmonella enterica were recovered. CONCLUSIONS AND RELEVANCE: These results suggest the potential of metagenomics as a culture-independent approach for the identification of bacterial pathogens during an outbreak of diarrheal disease. Challenges include improving diagnostic sensitivity, speeding up and simplifying workflows, and reducing costs.

Lineage skewing and genome instability underlie marrow failure in a zebrafish model of GATA2 deficiency.

Inherited bone marrow failure associated with heterozygous mutations in GATA2 predisposes toward hematological malignancies, but the mechanisms remain poorly understood. Here, we investigate the mechanistic basis of marrow failure in a zebrafish model of GATA2 deficiency. Single-cell transcriptomics and chromatin accessibility assays reveal that loss of gata2a leads to skewing toward the erythroid lineage at the expense of myeloid cells, associated with loss of cebpa expression and decreased PU.1 and CEBPA transcription factor accessibility in hematopoietic stem and progenitor cells (HSPCs). Furthermore, gata2a mutants show impaired expression of npm1a, the zebrafish NPM1 ortholog. Progressive loss of npm1a in HSPCs is associated with elevated levels of DNA damage in gata2a mutants. Thus, Gata2a maintains myeloid lineage priming through cebpa and protects against genome instability and marrow failure by maintaining expression of npm1a. Our results establish a potential mechanism underlying bone marrow failure in GATA2 deficiency.

Deletion of TnAbaR23 results in both expected and unexpected antibiogram changes in a multidrug-resistant Acinetobacter baumannii strain.

Since the 2006 discovery of the Acinetobacter baumannii strain AYE AbaR1 resistance island, similar elements have been reported in numerous members of this species. As AbaR1 is distantly related to Tn7, we have renamed it TnAbaR1. TnAbaR transposons are known to carry multiple antibiotic resistance- and efflux-associated genes, although none have been experimentally studied en bloc. We deleted the TnAbaR transposon in A. baumannii A424, which we have designated TnAbaR23, and characterized independent deletion mutants DCO163 and DCO174. The NotI pulsed-field gel electrophoresis (PFGE) profile of strain DCO174 was consistent with targeted deletion of TnAbaR23 alone, but strain DCO163 apparently harbored a second large genomic deletion. Nevertheless, "subtractive amplification" targeting 52 TnAbaR and/or resistance-associated loci yielded identical results for both mutants and highlighted genes lost relative to strain A424. PCR mapping and genome sequencing revealed the entire 48.3-kb sequence of TnAbaR23. Consistent with TnAbaR23 carrying two copies of sul1, both mutants exhibited markedly increased susceptibility to sulfamethoxazole. In contrast, loss of tetAR(A) resulted in only a minor and variable increase in tetracycline susceptibility. Despite not exhibiting a growth handicap, strain DCO163 was more susceptible than strain DCO174 to 9 of 10 antibiotics associated with mutant-to-mutant variation in susceptibility, suggesting impairment of an undefined resistance-associated function. Remarkably, despite all three strains sharing identical gyrA and parC sequences, the ciprofloxacin MIC of DCO174 was >8-fold that of DCO163 and A424, suggesting a possible paradoxical role for TnAbaR23 in promoting sensitivity to ciprofloxacin. This study highlights the importance of experimental scrutiny and challenges the assumption that resistance phenotypes can reliably be predicted from genotypes alone.

Genomic analysis uncovers a phenotypically diverse but genetically homogeneous Escherichia coli ST131 clone circulating in unrelated urinary tract infections.

OBJECTIVES: To determine variation at the genome level in Escherichia coli ST131 clinical isolates previously shown to be phenotypically diverse. METHODS: The genomes of 10 ST131 isolates extensively characterized in previous studies were sequenced using combinations of Illumina and 454 sequencing technology. Whole-genome comparisons and phylogenetic comparisons were then performed across the strain set and with other closely related extraintestinal pathogenic E. coli (ExPEC) strain types. RESULTS: E. coli ST131 is overrepresented in a collection of clinical isolates, and there is large phenotypic variation amongst isolates. In contrast, genome sequencing of a selection of non-related clinical isolates shows almost no genetic variation between ST131 strains, and E. coli ST131 shows evidence of a genetically monomorphic pathogen showing a similar evolutionary trend to hypervirulent Clostridium difficile. CONCLUSIONS: A dominant circulating clone of E. coli ST131 has been identified in unrelated clinical urine samples in the UK. The clone splits into two distinct subgroups on the basis of antimicrobial resistance levels and carriage of extended-spectrum ß-lactamase plasmids. This provides the most comprehensive snapshot to date of the true molecular epidemiology of ST131 clinical isolates.

Defining bacterial species in the genomic era: insights from the genus Acinetobacter.

BACKGROUND: Microbial taxonomy remains a conservative discipline, relying on phenotypic information derived from growth in pure culture and techniques that are time-consuming and difficult to standardize, particularly when compared to the ease of modern high-throughput genome sequencing. Here, drawing on the genus Acinetobacter as a test case, we examine whether bacterial taxonomy could abandon phenotypic approaches and DNA-DNA hybridization and, instead, rely exclusively on analyses of genome sequence data. RESULTS: In pursuit of this goal, we generated a set of thirteen new draft genome sequences, representing ten species, combined them with other publically available genome sequences and analyzed these 38 strains belonging to the genus. We found that analyses based on 16S rRNA gene sequences were not capable of delineating accepted species. However, a core genome phylogenetic tree proved consistent with the currently accepted taxonomy of the genus, while also identifying three misclassifications of strains in collections or databases. Among rapid distance-based methods, we found average-nucleotide identity (ANI) analyses delivered results consistent with traditional and phylogenetic classifications, whereas gene content based approaches appear to be too strongly influenced by the effects of horizontal gene transfer to agree with previously accepted species. CONCLUSION: We believe a combination of core genome phylogenetic analysis and ANI provides an appropriate method for bacterial species delineation, whereby bacterial species are defined as monophyletic groups of isolates with genomes that exhibit at least 95% pair-wise ANI. The proposed method is backwards compatible; it provides a scalable and uniform approach that works for both culturable and non-culturable species; is faster and cheaper than traditional taxonomic methods; is easily replicable and transferable among research institutions; and lastly, falls in line with Darwin's vision of classification becoming, as far as is possible, genealogical.

Posttranscriptional Regulation in Response to Different Environmental Stresses in Campylobacter jejuni.

The survival strategies that Campylobacter jejuni (C. jejuni) employ throughout its transmission and infection life cycles remain largely elusive. Specifically, there is a lack of understanding about the posttranscriptional regulation of stress adaptations resulting from small noncoding RNAs (sRNAs). Published C. jejuni sRNAs have been discovered in specific conditions but with limited insights into their biological activities. Many more sRNAs are yet to be discovered as they may be condition-dependent. Here, we have generated transcriptomic data from 21 host- and transmission-relevant conditions. The data uncovered transcription start sites, expression patterns and posttranscriptional regulation during various stress conditions. This data set helped predict a list of putative sRNAs. We further explored the sRNAs' biological functions by integrating differential gene expression analysis, coexpression analysis, and genome-wide sRNA target prediction. The results showed that the C. jejuni gene expression was influenced primarily by nutrient deprivation and food storage conditions. Further exploration revealed a putative sRNA (CjSA21) that targeted tlp1 to 4 under food processing conditions. tlp1 to 4 are transcripts that encode methyl-accepting chemotaxis proteins (MCPs), which are responsible for chemosensing. These results suggested CjSA21 inhibits chemotaxis and promotes survival under food processing conditions. This study presents the broader research community with a comprehensive data set and highlights a novel sRNA as a potential chemotaxis inhibitor. IMPORTANCE The foodborne pathogen C. jejuni is a significant challenge for the global health care system. It is crucial to investigate C. jejuni posttranscriptional regulation by small RNAs (sRNAs) in order to understand how it adapts to different stress conditions. However, limited data are available for investigating sRNA activity under stress. In this study, we generate gene expression data of C. jejuni under 21 stress conditions. Our data analysis indicates that one of the novel sRNAs mediates the adaptation to food processing conditions. Results from our work shed light on the posttranscriptional regulation of C. jejuni and identify an sRNA associated with food safety.

Multifunctional Composite Hydrogels for Bacterial Capture, Growth/Elimination, and Sensing Applications.

Hydrogels are cross-linked networks of hydrophilic polymer chains with a three-dimensional structure. Owing to their unique features, the application of hydrogels for bacterial/antibacterial studies and bacterial infection management has grown in importance in recent years. This trend is likely to continue due to the rise in bacterial infections and antimicrobial resistance. By exploiting their physicochemical characteristics and inherent nature, hydrogels have been developed to achieve bacterial capture and detection, bacterial growth or elimination, antibiotic delivery, or bacterial sensing. Traditionally, the development of hydrogels for bacterial/antibacterial studies has focused on achieving a single function such as antibiotic delivery, antibacterial activity, bacterial growth, or bacterial detection. However, recent studies demonstrate the fabrication of multifunctional hydrogels, where a single hydrogel is capable of performing more than one bacterial/antibacterial function, or composite hydrogels consisting of a number of single functionalized hydrogels, which exhibit bacterial/antibacterial function synergistically. In this review, we first highlight the hydrogel features critical for bacterial studies and infection management. Then, we specifically address unique hydrogel properties, their surface/network functionalization, and their mode of action for bacterial capture, adhesion/growth, antibacterial activity, and bacterial sensing, respectively. Finally, we provide insights into different strategies for developing multifunctional hydrogels and how such systems can help tackle, manage, and understand bacterial infections and antimicrobial resistance. We also note that the strategies highlighted in this review can be adapted to other cell types and are therefore likely to find applications beyond the field of microbiology.

Bridging the N-terminal and middle domains in FliG of the flagellar rotor.

Flagella are necessary for bacterial movement and contribute to various aspects of virulence. They are complex cylindrical structures built of multiple molecular rings with self-assembly properties. The flagellar rotor is composed of the MS-ring and the C-ring. The FliG protein of the C-ring is central to flagellar assembly and function due to its roles in linking the C-ring with the MS-ring and in torque transmission from stator to rotor. No high-resolution structure of an assembled C-ring has been resolved to date, and the conformation adopted by FliG within the ring is unclear due to variations in available crystallographic data. Here, we use molecular dynamics (MD) simulations to study the conformation and dynamics of FliG in different states of assembly, including both in physiologically relevant and crystallographic lattice environments. We conclude that the linker between the FliG N-terminal and middle domain likely adopts an extended helical conformation in vivo, in contrast with the contracted conformation observed in some previous X-ray studies. We further support our findings with integrative model building of full-length FliG and a FliG ring model that is compatible with cryo-electron tomography (cryo-ET) and electron microscopy (EM) densities of the C-ring. Collectively, our study contributes to a better mechanistic understanding of the flagellar rotor assembly and its function.

Genomic Assembly of Clinical Candida glabrata (Nakaseomyces glabrata) Isolates Reveals within-Species Structural Plasticity and Association with In Vitro Antifungal Susceptibility.

The opportunistic human pathogen Candida glabrata has become an increasingly important threat to human health, with infections globally characterized by high mortality rates and multidrug resistance. To face this threat, more efficient diagnostic and therapeutic approaches are required, underpinning research to help define the intraspecies epidemiology, genetic variability, and therefore, diagnostic and therapeutic target stability. Previous comparative genetics studies conducted on limited numbers of strains only revealed partial resolution of chromosomal settings. In this study, by combining short- and long-read genome sequencing, phenotypic characterization, and comparative genomics over a large set of strains, we detected strict relationships between large chromosomal rearrangements and phylogenetic clades, genes subjected to different selective pressures, and new sets of genes associated with resistance to antifungals. Overall, these results not only provide a fundamental contribution to our knowledge of C. glabrata evolution and epidemiology but may also lay the foundations for the future development of tailored therapeutic approaches. IMPORTANCE The human pathogen Candida glabrata has become a global threat to human health, with infections characterized by high mortality and multidrug resistance. We have obtained nine fully assembled genomes from clinical isolates through a combination of short- and long-read sequencing approaches. The quality and completeness of such genomes and their subsequent comparison to the broadest set of genomes so far allowed us to pinpoint chromosomal rearrangements in several genomes and detect phylogenetic clades that were not associated with geographic location or isolation source. We identified a new set of genes associated with resistance to antifungals coding for adhesin or adhesin-like proteins, suggesting C. glabrata resists antifungals by forming aggregates or adhering to the host tissue. These results, which provide a fundamental contribution to our knowledge of C. glabrata evolution and epidemiology, may initiate the development of precision medicine interventions for patients with suspected or proven invasive fungal infections.

Genome sequences of three Acinetobacter baumannii strains assigned to the multilocus sequence typing genotypes ST2, ST25, and ST78.

Acinetobacter baumannii is an emerging opportunistic gram-negative pathogen responsible for hospital-acquired infections. A. baumannii epidemics described in Europe and worldwide were caused by a limited number of genotypic clusters of multidrug-resistant strains. Here, we report the availability of draft genome sequences for three multidrug-resistant A. baumannii strains assigned to multilocus sequence typing genotypes ST2, ST25, and ST78 that were more frequently isolated during outbreaks occurred in Greece, Italy, Lebanon, and Turkey.

Complete genome sequence of the Crohn's disease-associated adherent-invasive Escherichia coli strain HM605.

Adherent-invasive Escherichia coli strains are increasingly being associated with intestinal pathologies. Here we present the genome sequence of E. coli HM605, a strain isolated from colonic biopsy specimens of a patient with Crohn's disease.

Transcriptome analysis of avian pathogenic Escherichia coli O1 in chicken serum reveals adaptive responses to systemic infection.

Infections of avian pathogenic Escherichia coli (APEC) result in annual multimillion-dollar losses to the poultry industry. Despite this, little is known about the mechanisms by which APEC survives and grows in the bloodstream. Thus, the aim of this study was to identify molecular mechanisms enabling APEC to survive and grow in this critical host environment. To do so, we compared the transcriptome of APEC O1 during growth in Luria-Bertani broth and chicken serum. Several categories of genes, predicted to contribute to adaptation and growth in the avian host, were identified. These included several known virulence genes and genes involved in adaptive metabolism, protein transport, biosynthesis pathways, stress resistance, and virulence regulation. Several genes with unknown function, which were localized to pathogenicity islands or APEC O1's large virulence plasmid, pAPEC-O1-ColBM, were also identified, suggesting that they too contribute to survival in serum. The significantly upregulated genes dnaK, dnaJ, phoP, and ybtA were subsequently subjected to mutational analysis to confirm their role in conferring a competitive advantage during infection. This genome-wide analysis provides novel insight into processes that are important to the pathogenesis of APEC O1.

Modification of Bacteriophages to Increase Their Association with Lung Epithelium Cells In Vitro.

There is currently a renaissance in research on bacteriophages as alternatives to antibiotics. Phage specificity to their bacterial host, in addition to a plethora of other advantages, makes them ideal candidates for a broad range of applications, including bacterial detection, drug delivery, and phage therapy in particular. One issue obstructing phage efficiency in phage therapy settings is their poor localization to the site of infection in the human body. Here, we engineered phage T7 with lung tissue targeting homing peptides. We then used in vitro studies to demonstrate that the engineered T7 phages had a more significant association with the lung epithelium cells than wild-type T7. In addition, we showed that, in general, there was a trend of increased association of engineered phages with the lung epithelium cells but not mouse fibroblast cells, allowing for targeted tissue specificity. These results indicate that appending phages with homing peptides would potentially allow for greater phage concentrations and greater efficacy at the infection site.

YieJ (CbrC) mediates CreBC-dependent colicin E2 tolerance in Escherichia coli.

Colicin E2-tolerant (known as Cet2) Escherichia coli K-12 mutants overproduce an inner membrane protein, CreD, which is believed to cause the Cet2 phenotype. Here, we show that overproduction of CreD in a Cet2 strain results from hyperactivation of the CreBC two-component regulator, but CreD overproduction is not responsible for the Cet2 phenotype. Through microarray analysis and gene knockout and overexpression studies, we show that overexpression of another CreBC-regulated gene, yieJ (also known as cbrC), causes the Cet2 phenotype.

Comparison of CRISPR and Marker-Based Methods for the Engineering of Phage T7.

With the recent rise in interest in using lytic bacteriophages as therapeutic agents, there is an urgent requirement to understand their fundamental biology to enable the engineering of their genomes. Current methods of phage engineering rely on homologous recombination, followed by a system of selection to identify recombinant phages. For bacteriophage T7, the host genes cmk or trxA have been used as a selection mechanism along with both type I and II CRISPR systems to select against wild-type phage and enrich for the desired mutant. Here, we systematically compare all three systems; we show that the use of marker-based selection is the most efficient method and we use this to generate multiple T7 tail fibre mutants. Furthermore, we found the type II CRISPR-Cas system is easier to use and generally more efficient than a type I system in the engineering of phage T7. These results provide a foundation for the future, more efficient engineering of bacteriophage T7.

Exposure of Escherichia coli and Salmonella enterica serovar Typhimurium to triclosan induces a species-specific response, including drug detoxification.

OBJECTIVES: The use of triclosan within various environments has been linked to the development of multiple drug resistance (MDR) through the increased expression of efflux pumps such as AcrAB-TolC. In this work, we investigate the effect of triclosan exposure in order to ascertain the response of two species to the presence of this widely used biocide. METHODS: The transcriptomes of Salmonella enterica serovar Typhimurium SL1344 and Escherichia coli K-12 MG1655 after exposure to the MIC of triclosan (0.12 mg/L) were determined in microarray experiments. Phenotypic validation of the transcriptomic data included RT-PCR, ability to form a biofilm and motility assays. RESULTS: Despite important differences in the triclosan-dependent transcriptomes of the two species, increased expression of efflux pump component genes was seen in both. Increased expression of soxS was observed in Salmonella Typhimurium, however, within E. coli, decreased expression was seen. Expression of fabBAGI in Salmonella Typhimurium was decreased, whereas in E. coli expression of fabABFH was increased. Increased expression of ompR and genes within this regulon (e.g. ompC, csgD and ssrA) was seen in the transcriptome of Salmonella Typhimurium. An unexpected response of E. coli was the differential expression of genes within operons involved in iron homeostasis; these included fhu, fep and ent. CONCLUSIONS: These data indicate that whilst a core response to triclosan exposure exists, the differential transcriptome of each species was different. This suggests that E. coli K-12 should not be considered the paradigm for the Enterobacteriaceae when exploring the effects of antimicrobial agents.

Global responses of Escherichia coli to adverse conditions determined by microarrays and FT-IR spectroscopy.

The global gene expression and biomolecular composition in an Escherichia coli model strain exposed to 10 adverse conditions (sodium chloride, ethanol, glycerol, hydrochloric and acetic acid, sodium hydroxide, heat (46 degrees C), and cold (15 degrees C), as well as ethidium bromide and the disinfectant benzalkonium chloride) were determined using DNA microarrays and Fourier transform infrared (FT-IR) spectroscopy. In total, approximately 40% of all investigated genes (1682/4279 genes) significantly changed expression, compared with a nonstressed control. There were, however, only 3 genes (ygaW (unknown function), rmf (encoding a ribosomal modification factor), and ghrA (encoding a glyoxylate/hydroxypyruvate reductase)) that significantly changed expression under all conditions (not including benzalkonium chloride). The FT-IR analysis showed an increase in unsaturated fatty acids during ethanol and cold exposure, and a decrease during acid and heat exposure. Cold conditions induced changes in the carbohydrate composition of the cell, possibly related to the upregulation of outer membrane genes (glgAP and rcsA). Although some covariance was observed between the 2 data sets, principle component analysis and regression analyses revealed that the gene expression and the biomolecular responses are not well correlated in stressed populations of E. coli, underlining the importance of multiple strategies to begin to understand the effect on the whole cell.