Small stress response proteins in Escherichia coli: proteins missed by classical proteomic studies

SearchDOGS bacteria, software that provides automated identification of potentially missed genes in annotated bacterial genomes

We report the development of SearchDOGS Bacteria, software to automatically detect missing genes in annotated bacterial genomes by combining BLAST searches with comparative genomics. Having successfully applied the approach to yeast genomes, we redeveloped SearchDOGS to function as a standalone, downloadable package, requiring only a set of GenBank annotation files as input. The software automatically generates a homology structure using reciprocal BLAST and a synteny-based method; this is followed by a scan of the entire genome of each species for unannotated genes. Results are provided in a HTML interface, providing coordinates, BLAST results, syntenic location, omega values (Ka/Ks, where Ks is the number of synonymous substitutions per synonymous site and Ka is the number of nonsynonymous substitutions per nonsynonymous site) for protein conservation estimates, and other information for each candidate gene. Using SearchDOGS Bacteria, we identified 155 gene candidates in the Shigella boydii sb227 genome, including 56 candidates of length < 60 codons. SearchDOGS Bacteria has two major advantages over currently available annotation software. First, it outperforms current methods in terms of sensitivity and is highly effective at identifying small or highly diverged genes. Second, as a freely downloadable package, it can be used with unpublished or confidential data.

The condition-dependent transcriptional landscape of Burkholderia pseudomallei

Burkholderia pseudomallei (Bp), the causative agent of the often-deadly infectious disease melioidosis, contains one of the largest prokaryotic genomes sequenced to date, at 7.2 Mb with two large circular chromosomes (1 and 2). To comprehensively delineate the Bp transcriptome, we integrated whole-genome tiling array expression data of Bp exposed to >80 diverse physical, chemical, and biological conditions. Our results provide direct experimental support for the strand-specific expression of 5,467 Sanger protein-coding genes, 1,041 operons, and 766 non-coding RNAs. A large proportion of these transcripts displayed condition-dependent expression, consistent with them playing functional roles. The two Bp chromosomes exhibited dramatically different transcriptional landscapes — Chr 1 genes were highly and constitutively expressed, while Chr 2 genes exhibited mosaic expression where distinct subsets were expressed in a strongly condition-dependent manner. We identified dozens of cis-regulatory motifs associated with specific condition-dependent expression programs, and used the condition compendium to elucidate key biological processes associated with two complex pathogen phenotypes — quorum sensing and in vivo infection. Our results demonstrate the utility of a Bp condition-compendium as a community resource for biological discovery. Moreover, the observation that significant portions of the Bp virulence machinery can be activated by specific in vitro cues provides insights into Bp's capacity as an “accidental pathogen”, where genetic pathways used by the bacterium to survive in environmental niches may have also facilitated its ability to colonize human hosts.

Bacterial transcriptomes are dynamic, context-specific and condition-dependent. Infection by the soil bacterium, Burkholderia pseudomallei, causes melioidosis, an often fatal infectious disease of humans and animals. Possessing a large multi-chromosomal genome, B. pseudomallei is able to persist and survive in a multitude of environments. To obtain a comprehensive overview of B. pseudomallei expressed transcripts, we initiated whole-genome transcriptome profiling covering a broad spectrum of conditions and exposures — a so-called “condition compendium”. Using the compendium, we confirmed many previously-annotated genes and operons, and also identified hundreds of novel transcripts including anti-sense transcripts and non-coding RNAs. By systematically examining genes exhibiting highly similar expression patterns, we ascribed putative functions to previously uncharacterized genes, and identified novel regulatory elements controlling these expression patterns. We also used the compendium to elucidate candidate virulence pathways associated with quorum-sensing and infection in mice. Our study showcases the power of a B. pseudomallei condition compendium as a valuable resource for understanding microbial physiology and the pathogenesis of melioidosis.

A component of gamma-radiation-induced cell death in E. coli is programmed and interlinked with activation of caspase-3 and SOS response

The current study deals with the molecular mechanism of radiation-induced cell death (RICD) in Escherichia coli. Irradiated E. coli cells displayed markers similar to those found in eukaryotic programmed cell death (PCD) such as caspase-3 activation and phosphatidylserine externalization. RICD was found to be suppressed upon pretreatment with sublethal concentrations of rifampicin or chloramphenicol, indicating the requirement of de novo gene expression. RICD was also found to be inhibited by cell permeable inhibitors of caspase-3 or poly (ADP-ribose) polymerase, indicating the involvement of PCD during RICD in E. coli. Radiation-induced SOS response was alleviated as observed with decrease in LexA level and also reduced cell filamentation frequency in the presence of caspase inhibitor. Further, the inhibitor-mediated rescue was not observed in single-gene knockouts of umuC, umuD, recB and ruvA, the genes which are associated with SOS response. This implies a linkage between SOS response and PCD in radiation-exposed E. coli cells.

The Escherichia coli CydX protein is a member of the CydAB cytochrome bd oxidase complex and is required for cytochrome bd oxidase activity

Cytochrome bd oxidase operons from more than 50 species of bacteria contain a short gene encoding a small protein that ranges from ∼30 to 50 amino acids and is predicted to localize to the cell membrane. Although cytochrome bd oxidases have been studied for more than 70 years, little is known about the role of this small protein, denoted CydX, in oxidase activity. Here we report that Escherichia coli mutants lacking CydX exhibit phenotypes associated with reduced oxidase activity. In addition, cell membrane extracts from ΔcydX mutant strains have reduced oxidase activity in vitro. Consistent with data showing that CydX is required for cytochrome bd oxidase activity, copurification experiments indicate that CydX interacts with the CydAB cytochrome bd oxidase complex. Together, these data support the hypothesis that CydX is a subunit of the CydAB cytochrome bd oxidase complex that is required for complex activity. The results of mutation analysis of CydX suggest that few individual amino acids in the small protein are essential for function, at least in the context of protein overexpression. In addition, the results of analysis of the paralogous small transmembrane protein AppX show that the two proteins could have some overlapping functionality in the cell and that both have the potential to interact with the CydAB complex.

The Escherichia coli Cpx envelope stress response regulates genes of diverse function that impact antibiotic resistance and membrane integrity

The Cpx envelope stress response mediates adaptation to stresses that cause envelope protein misfolding. Adaptation is partly conferred through increased expression of protein folding and degradation factors. The Cpx response also plays a conserved role in the regulation of virulence determinant expression and impacts antibiotic resistance. We sought to identify adaptive mechanisms that may be involved in these important functions by characterizing changes in the transcriptome of two different Escherichia coli strains when the Cpx response is induced. We show that, while there is considerable strain- and condition-specific variability in the Cpx response, the regulon is enriched for proteins and functions that are inner membrane associated under all conditions. Genes that were changed by Cpx pathway induction under all conditions were involved in a number of cellular functions and included several intergenic regions, suggesting that posttranscriptional regulation is important during Cpx-mediated adaptation. Some Cpx-regulated genes are centrally involved in energetics and play a role in antibiotic resistance. We show that a number of small, uncharacterized envelope proteins are Cpx regulated and at least two of these affect phenotypes associated with membrane integrity. Altogether, our work suggests new mechanisms of Cpx-mediated envelope stress adaptation and antibiotic resistance.

Overexpression of the KdpF membrane peptide in Mycobacterium bovis BCG results in reduced intramacrophage growth and altered cording morphology

Membrane peptides appear as an emerging class of regulatory molecules in bacteria, which can interact with membrane proteins, such as sensor kinases. To date, regulatory membrane peptides have been completely overlooked in mycobacteria. The 30 amino-acid-long KdpF peptide, which is co-transcribed with kdpABC genes and regulated by the KdpDE two-component system, is supposed to stabilize the KdpABC potassium transporter complex but may also exhibit unsuspected regulatory function(s) towards the KdpD sensor kinase. Herein, we showed by quantitative RT-PCR that the Mycobacterium bovis BCG kdpAB and kdpDE genes clusters are differentially induced in potassium-deprived broth medium or within infected macrophages. We have overexpressed the kdpF gene in M. bovis BCG to investigate its possible regulatory role and effect on mycobacterial virulence. Our results indicate that KdpF does not play a critical regulatory role on kdp genes expression despite the fact that KdpF interacts with the KdpD sensor kinase in a bacterial two-hybrid assay. However, overexpression of kdpF results in a significant reduction of M. bovis BCG growth in both murine and human primary macrophages, and is associated with a strong alteration of colonial morphology and impaired cording formation. To identify novel KdpF interactants, a mycobacterial library was screened using KdpF as bait in the bacterial two-hybrid system. This allowed us to identify members of the MmpL family of membrane proteins, known to participate in the biosynthesis/transport of various cell wall lipids, thus highlighting a possible link between KdpF and cell wall lipid metabolism. Taken together, these data suggest that KdpF overexpression reduces intramacrophage growth which may result from alteration of the mycobacterial cell wall.

Single transmembrane peptide DinQ modulates membrane-dependent activities

The functions of several SOS regulated genes in Escherichia coli are still unknown, including dinQ. In this work we characterize dinQ and two small RNAs, agrA and agrB, with antisense complementarity to dinQ. Northern analysis revealed five dinQ transcripts, but only one transcript (+44) is actively translated. The +44 dinQ transcript translates into a toxic single transmembrane peptide localized in the inner membrane. AgrB regulates dinQ RNA by RNA interference to counteract DinQ toxicity. Thus the dinQ-agr locus shows the classical features of a type I TA system and has many similarities to the tisB-istR locus. DinQ overexpression depolarizes the cell membrane and decreases the intracellular ATP concentration, demonstrating that DinQ can modulate membrane-dependent processes. Augmented DinQ strongly inhibits marker transfer by Hfr conjugation, indicating a role in recombination. Furthermore, DinQ affects transformation of nucleoid morphology in response to UV damage. We hypothesize that DinQ is a transmembrane peptide that modulates membrane-dependent activities such as nucleoid compaction and recombination.

Exposure of the bacterium Escherichia coli to DNA damaging agents induces the SOS response, which up-regulates gene functions involved in numerous cellular processes such as DNA repair, cell division, and replication. Most of the SOS regulated genes in E. coli have been characterized, but still there are several genes of unknown function. One of these uncharacterized genes is dinQ. In this work we characterize dinQ and two novel small RNAs, agrA and agrB, that regulate expression of dinQ. The DinQ peptide is localized in the inner membrane as a single transmembrane peptide of 27 amino acids. Small proteins of less than 50 amino acids are important in cellular processes such as regulation, signalling, and antibacterial action. Here we demonstrate that DinQ modulates recombination and transformation of nucleoid morphology in response to UV damage. Our results provide new insights into small hydrophobic peptides that could regulate important DNA metabolic processes dependent on the inner membrane of the cell.

EcoGene 3.0

EcoGene (http://ecogene.org) is a database and website devoted to continuously improving the structural and functional annotation of Escherichia coli K-12, one of the most well understood model organisms, represented by the MG1655(Seq) genome sequence and annotations. Major improvements to EcoGene in the past decade include (i) graphic presentations of genome map features; (ii) ability to design Boolean queries and Venn diagrams from EcoArray, EcoTopics or user-provided GeneSets; (iii) the genome-wide clone and deletion primer design tool, PrimerPairs; (iv) sequence searches using a customized EcoBLAST; (v) a Cross Reference table of synonymous gene and protein identifiers; (vi) proteome-wide indexing with GO terms; (vii) EcoTools access to >2000 complete bacterial genomes in EcoGene-RefSeq; (viii) establishment of a MySql relational database; and (ix) use of web content management systems. The biomedical literature is surveyed daily to provide citation and gene function updates. As of September 2012, the review of 37 397 abstracts and articles led to creation of 98 425 PubMed-Gene links and 5415 PubMed-Topic links. Annotation updates to Genbank U00096 are transmitted from EcoGene to NCBI. Experimental verifications include confirmation of a CTG start codon, pseudogene restoration and quality assurance of the Keio strain collection.

A connecter-like factor, CacA, links RssB/RpoS and the CpxR/CpxA two-component system in Salmonella

Background

Bacteria integrate numerous environmental stimuli when generating cellular responses. Increasing numbers of examples describe how one two-component system (TCS) responds to signals detected by the sensor of another TCS. However, the molecular mechanisms underlying this phenomenon remain poorly defined.

Results

Here, we report a connector-like factor that affects the activity of the CpxR/CpxA two-component system in Salmonella enterica serovar Typhimurium. We isolated a clone that induced the expression of a cpxP-lac gene fusion from a high-copy-number plasmid pool of random Salmonella genomic fragments. A 63-amino acid protein, CacA, was responsible for the CpxA/CpxR-dependent activation of the cpxP gene. The CpxR-activated genes cpxP and spy exhibited approximately 30% and 50% reductions in transcription, respectively, in a clean cacA deletion mutant strain in comparison to wild-type. From 33 response regulator (RR) deletion mutants, we identified that the RssB regulator represses cacA transcription. Substitution mutations in a conserved -10 region harboring the RNA polymerase recognition sequence, which is well conserved with a known RpoS -10 region consensus sequence, rendered the cacA promoter RpoS-independent. The CacA-mediated induction of cpxP transcription was affected in a trxA deletion mutant, which encodes thioredoxin 1, suggesting a role for cysteine thiol-disulfide exchange(s) in CacA-dependent Cpx activation.

Conclusions

We identified CacA as an activator of the CpxR/CpxA system in the plasmid clone. We propose that CacA may integrate the regulatory status of RssB/RpoS into the CpxR/CpxA system. Future investigations are necessary to thoroughly elucidate how CacA activates the CpxR/CpxA system.

Conserved small protein associates with the multidrug efflux pump AcrB and differentially affects antibiotic resistance

The AcrAB-TolC multidrug efflux pump confers resistance to a wide variety of antibiotics and other compounds in Escherichia coli. Here we show that AcrZ (formerly named YbhT), a 49-amino-acid inner membrane protein, associates with the AcrAB-TolC complex. Co-purification of AcrZ with AcrB, in the absence of both AcrA and TolC, two-hybrid assays and suppressor mutations indicate that this interaction occurs through the inner membrane protein AcrB. The highly conserved acrZ gene is coregulated with acrAB through induction by the MarA, Rob, and SoxS transcription regulators. In addition, mutants lacking AcrZ are sensitive to many, but not all, of the antibiotics transported by AcrAB-TolC. This differential antibiotic sensitivity suggests that AcrZ may enhance the ability of the AcrAB-TolC pump to export certain classes of substrates.

Predicting statistical properties of open reading frames in bacterial genomes

An analytical model based on the statistical properties of Open Reading Frames (ORFs) of eubacterial genomes such as codon composition and sequence length of all reading frames was developed. This new model predicts the average length, maximum length as well as the length distribution of the ORFs of 70 species with GC contents varying between 21% and 74%. Furthermore, the number of annotated genes is predicted with high accordance. However, the ORF length distribution in the five alternative reading frames shows interesting deviations from the predicted distribution. In particular, long ORFs appear more often than expected statistically. The unexpected depletion of stop codons in these alternative open reading frames cannot completely be explained by a biased codon usage in the +1 frame. While it is unknown if the stop codon depletion has a biological function, it could be due to a protein coding capacity of alternative ORFs exerting a selection pressure which prevents the fixation of stop codon mutations. The comparison of the analytical model with bacterial genomes, therefore, leads to a hypothesis suggesting novel gene candidates which can now be investigated in subsequent wet lab experiments.

The small protein CydX is required for function of cytochrome bd oxidase in Brucella abortus

A large number of hypothetical genes potentially encoding small proteins of unknown function are annotated in the Brucella abortus genome. Individual deletion of 30 of these genes identified four mutants, in BAB1_0355, BAB2_0726, BAB2_0470, and BAB2_0450 that were highly attenuated for infection. BAB2_0726, an YbgT-family protein located at the 3′ end of the cydAB genes encoding cytochrome bd ubiquinal oxidase, was designated cydX. A B. abortus cydX mutant lacked cytochrome bd oxidase activity, as shown by increased sensitivity to H₂O₂, decreased acid tolerance and increased resistance to killing by respiratory inhibitors. The C terminus, but not the N terminus, of CydX was located in the periplasm, suggesting that CydX is an integral cytoplasmic membrane protein. Phenotypic analysis of the cydX mutant, therefore, suggested that CydX is required for full function of cytochrome bd oxidase, possibly via regulation of its assembly or activity.

Deep sequencing uncovers numerous small RNAs on all four replicons of the plant pathogen Agrobacterium tumefaciens

Agrobacterium species are capable of interkingdom gene transfer between bacteria and plants. The genome of Agrobacterium tumefaciens consists of a circular and a linear chromosome, the At-plasmid and the Ti-plasmid, which harbors bacterial virulence genes required for tumor formation in plants. Little is known about promoter sequences and the small RNA (sRNA) repertoire of this and other α-proteobacteria. We used a differential RNA sequencing (dRNA-seq) approach to map transcriptional start sites of 388 annotated genes and operons. In addition, a total number of 228 sRNAs was revealed from all four Agrobacterium replicons. Twenty-two of these were confirmed by independent RNA gel blot analysis and several sRNAs were differentially expressed in response to growth media, growth phase, temperature or pH. One sRNA from the Ti-plasmid was massively induced under virulence conditions. The presence of 76 cis-antisense sRNAs, two of them on the reverse strand of virulence genes, suggests considerable antisense transcription in Agrobacterium. The information gained from this study provides a valuable reservoir for an in-depth understanding of sRNA-mediated regulation of the complex physiology and infection process of Agrobacterium.

Characterization of Zn(II)-responsive ribosomal proteins YkgM and L31 in E. coli

RT-PCR and DNA microarrays were used to probe for Zn(II)-responsive genes in E. coli cells that were made Zn(II) deficient. Microarray data revealed 114 genes were significantly up-regulated and 146 genes were significantly down-regulated in Zn(II) deficient conditions. The three most up-regulated genes were (1) znuA, which encodes for a periplasmic protein known to be involved with Zn(II) import, (2) yodA, which encodes for a periplasmic protein with unknown function, and (3) ykgM, which encodes for a ribosomal protein that is thought to be a paralog of ribosomal protein L31. YodA was over-expressed and purified as a maltose binding protein (MBP) fusion protein and shown to tightly bind 4 equivalents of Zn(II). Metal analyses showed that MBP-YkgM does not bind Zn(II). On the other hand, MBP-L31 tightly binds 1 equivalent of Zn(II). EXAFS studies on MBP-L31 suggest a ligand field of 1 histidine, 1 cysteine, and 2 additional N/O scatterers. Site-directed mutagenesis studies suggest that Cys16 coordinates Zn(II) in MBP-L31 and that the other three cysteines do not bind metal. These results are discussed in light of Zn(II) starvation model that has been postulated for B. subtilis.

Novel proteomic tools reveal essential roles of SRP and importance of proper membrane protein biogenesis

The signal recognition particle (SRP), which mediates cotranslational protein targeting to cellular membranes, is universally conserved and essential for bacterial and mammalian cells. However, the current understanding of the role of SRP in cell physiology and pathology is still poor, and the reasons behind its essential role in cell survival remain unclear. Here, we systematically analyzed the consequences of SRP loss in E. coli using time-resolved quantitative proteomic analyses. A series of snapshots of the steady-state and newly synthesized proteome unveiled three stages of cellular responses to SRP depletion, and demonstrated essential roles of SRP in metabolism, membrane potential, and protein and energy homeostasis in both the membrane and cytoplasm. We also identified a group of periplasmic proteins, including key molecular chaperones, whose localization was impaired by the loss of SRP; this and additional results showed that SRP is crucial for protein homeostasis in the bacterial envelope. These results reveal the extensive roles that SRP plays in bacterial physiology, emphasize the importance of proper membrane protein biogenesis, and demonstrate the ability of time-resolved quantitative proteomic analysis to provide new biological insights.

Using DNA microarrays to assay part function

In recent years, the capability of synthetic biology to design large genetic circuits has dramatically increased due to rapid advances in DNA synthesis technology and development of tools for large-scale assembly of DNA fragments. Large genetic circuits require more components (parts), especially regulators such as transcription factors, sigma factors, and viral RNA polymerases to provide increased regulatory capability, and also devices such as sensors, receivers, and signaling molecules. All these parts may have a potential impact upon the host that needs to be considered when designing and fabricating circuits. DNA microarrays are a well-established technique for global monitoring of gene expression and therefore are an ideal tool for systematically assessing the impact of expressing parts of genetic circuits in host cells. Knowledge of part impact on the host enables the user to design circuits from libraries of parts taking into account their potential impact and also to possibly modify the host to better tolerate stresses induced by the engineered circuit. In this chapter, we present the complete methodology of performing microarrays from choice of array platform, experimental design, preparing samples for array hybridization, and associated data analysis including preprocessing, normalization, clustering, identifying significantly differentially expressed genes, and interpreting the data based on known biology. With these methodologies, we also include lists of bioinformatic resources and tools for performing data analysis. The aim of this chapter is to provide the reader with the information necessary to be able to systematically catalog the impact of genetic parts on the host and also to optimize the operation of fully engineered genetic circuits.

Deep sequencing reveals as-yet-undiscovered small RNAs in Escherichia coli

BACKGROUND

In Escherichia coli, approximately 100 regulatory small RNAs (sRNAs) have been identified experimentally and many more have been predicted by various methods. To provide a comprehensive overview of sRNAs, we analysed the low-molecular-weight RNAs (< 200 nt) of E. coli with deep sequencing, because the regulatory RNAs in bacteria are usually 50-200 nt in length.

RESULTS

We discovered 229 novel candidate sRNAs (≥ 50 nt) with computational or experimental evidence of transcription initiation. Among them, the expression of seven intergenic sRNAs and three cis-antisense sRNAs was detected by northern blot analysis. Interestingly, five novel sRNAs are expressed from prophage regions and we note that these sRNAs have several specific characteristics. Furthermore, we conducted an evolutionary conservation analysis of the candidate sRNAs and summarised the data among closely related bacterial strains.

CONCLUSIONS

This comprehensive screen for E. coli sRNAs using a deep sequencing approach has shown that many as-yet-undiscovered sRNAs are potentially encoded in the E. coli genome. We constructed the Escherichia coli Small RNA Browser (ECSBrowser; http://rna.iab.keio.ac.jp/), which integrates the data for previously identified sRNAs and the novel sRNAs found in this study.

Precise manipulation of chromosomes in vivo enables genome-wide codon replacement

We present genome engineering technologies that are capable of fundamentally reengineering genomes from the nucleotide to the megabase scale. We used multiplex automated genome engineering (MAGE) to site-specifically replace all 314 TAG stop codons with synonymous TAA codons in parallel across 32 Escherichia coli strains. This approach allowed us to measure individual recombination frequencies, confirm viability for each modification, and identify associated phenotypes. We developed hierarchical conjugative assembly genome engineering (CAGE) to merge these sets of codon modifications into genomes with 80 precise changes, which demonstrate that these synonymous codon substitutions can be combined into higher-order strains without synthetic lethal effects. Our methods treat the chromosome as both an editable and an evolvable template, permitting the exploration of vast genetic landscapes.

A DNA damage checkpoint in Caulobacter crescentus inhibits cell division through a direct interaction with FtsW

Following DNA damage, cells typically delay cell cycle progression and inhibit cell division until their chromosomes have been repaired. The bacterial checkpoint systems responsible for these DNA damage responses are incompletely understood. Here, we show that Caulobacter crescentus responds to DNA damage by coordinately inducing an SOS regulon and inhibiting the master regulator CtrA. Included in the SOS regulon is sidA (SOS-induced inhibitor of cell division A), a membrane protein of only 29 amino acids that helps to delay cell division following DNA damage, but is dispensable in undamaged cells. SidA is sufficient, when overproduced, to block cell division. However, unlike many other regulators of bacterial cell division, SidA does not directly disrupt the assembly or stability of the cytokinetic ring protein FtsZ, nor does it affect the recruitment of other components of the cell division machinery. Instead, we provide evidence that SidA inhibits division by binding directly to FtsW to prevent the final constriction of the cytokinetic ring.

An integrated approach for finding overlooked genes in Shigella

BACKGROUND

The completion of numerous genome sequences introduced an era of whole-genome study. However, many genes are missed during genome annotation, including small RNAs (sRNAs) and small open reading frames (sORFs). In order to improve genome annotation, we aimed to identify novel sRNAs and sORFs in Shigella, the principal etiologic agents of bacillary dysentery.

METHODOLOGY/PRINCIPAL FINDINGS

We identified 64 sRNAs in Shigella, which were experimentally validated in other bacteria based on sequence conservation. We employed computer-based and tiling array-based methods to search for sRNAs, followed by RT-PCR and northern blots, to identify nine sRNAs in Shigella flexneri strain 301 (Sf301) and 256 regions containing possible sRNA genes. We found 29 candidate sORFs using bioinformatic prediction, array hybridization and RT-PCR verification. We experimentally validated 557 (57.9%) DOOR operon predictions in the chromosomes of Sf301 and 46 (76.7%) in virulence plasmid.We found 40 additional co-expressed gene pairs that were not predicted by DOOR.

CONCLUSIONS/SIGNIFICANCE

We provide an updated and comprehensive annotation of the Shigella genome. Our study increased the expected numbers of sORFs and sRNAs, which will impact on future functional genomics and proteomics studies. Our method can be used for large scale reannotation of sRNAs and sORFs in any microbe with a known genome sequence.

An expanding universe of small proteins

Historically, small proteins (sproteins) of less than 50 amino acids, in their final processed forms or genetically encoded as such, have been understudied. However, both serendipity and more recent focused efforts have led to the identification of a number of new sproteins in both Gram-negative and Gram-positive bacteria. Increasing evidence demonstrates that sproteins participate in a wide array of cellular processes and exhibit great diversity in their mechanisms of action, yet general principles of sprotein function are emerging. This review highlights examples of sproteins that participate in cell signaling, act as antibiotics and toxins, and serve as structural proteins. We also describe roles for sproteins in detecting and altering membrane features, acting as chaperones, and regulating the functions of larger proteins.

Identification of a chloroform-soluble membrane miniprotein in Escherichia coli and its homolog in Salmonella typhimurium

Two homologous 29 amino acid-long highly hydrophobic membrane miniproteins were identified in the Bligh-Dyer lipid extracts of Escherichia coli and Salmonella typhimurium using liquid chromatography/tandem mass spectrometry (LC/MS/MS). The amino acid sequences of the proteins were determined by collision-induced dissociation tandem mass spectrometry, in conjunction with a translating BLAST (tBLASTn) search, i.e., comparing the MS/MS-determined protein query sequence against the six-frame translations of the nucleotide sequences of the E. coli and S. typhimurium genomes. Further MS characterization revealed that both proteins retain the N-terminal initiating formyl-methionines. The methodologies described here may be amendable for detecting and characterizing small hydrophobic proteins in other organisms that are difficult to annotate or analyze by conventional methods.

Deep sequencing-based analysis of the anaerobic stimulon in Neisseria gonorrhoeae

BACKGROUND

Maintenance of an anaerobic denitrification system in the obligate human pathogen, Neisseria gonorrhoeae, suggests that an anaerobic lifestyle may be important during the course of infection. Furthermore, mounting evidence suggests that reduction of host-produced nitric oxide has several immunomodulary effects on the host. However, at this point there have been no studies analyzing the complete gonococcal transcriptome response to anaerobiosis. Here we performed deep sequencing to compare the gonococcal transcriptomes of aerobically and anaerobically grown cells. Using the information derived from this sequencing, we discuss the implications of the robust transcriptional response to anaerobic growth.

RESULTS

We determined that 198 chromosomal genes were differentially expressed (~10% of the genome) in response to anaerobic conditions. We also observed a large induction of genes encoded within the cryptic plasmid, pJD1. Validation of RNA-seq data using translational-lacZ fusions or RT-PCR demonstrated the RNA-seq results to be very reproducible. Surprisingly, many genes of prophage origin were induced anaerobically, as well as several transcriptional regulators previously unknown to be involved in anaerobic growth. We also confirmed expression and regulation of a small RNA, likely a functional equivalent of fnrS in the Enterobacteriaceae family. We also determined that many genes found to be responsive to anaerobiosis have also been shown to be responsive to iron and/or oxidative stress.

CONCLUSIONS

Gonococci will be subject to many forms of environmental stress, including oxygen-limitation, during the course of infection. Here we determined that the anaerobic stimulon in gonococci was larger than previous studies would suggest. Many new targets for future research have been uncovered, and the results derived from this study may have helped to elucidate factors or mechanisms of virulence that may have otherwise been overlooked.

Tiling array analysis of UV treated Escherichia coli predicts novel differentially expressed small peptides

Background

Despite comprehensive investigation, the Escherichia coli SOS response system is not yet fully understood. We have applied custom designed whole genome tiling arrays to measure UV invoked transcriptional changes in E. coli. This study provides a more complete insight into the transcriptome and the UV irradiation response of this microorganism.

Results

We detected a number of novel differentially expressed transcripts in addition to the expected SOS response genes (such as sulA, recN, uvrA, lexA, umuC and umuD) in the UV treated cells. Several of the differentially expressed transcripts might play important roles in regulation of the cellular response to UV damage. We have predicted 23 novel small peptides from our set of detected non-gene transcripts. Further, three of the predicted peptides were cloned into protein expression vectors to test the biological activity. All three constructs expressed the predicted peptides, in which two of them were highly toxic to the cell. Additionally, a remarkably high overlap with previously in-silico predicted non-coding RNAs (ncRNAs) was detected. Generally we detected a far higher transcriptional activity than the annotation suggests, and these findings correspond with previous transcription mappings from E. coli and other organisms.

Conclusions

Here we demonstrate that the E. coli transcriptome consists of far more transcripts than the present annotation suggests, of which many transcripts seem important to the bacterial stress response. Sequence alignment of promoter regions suggest novel regulatory consensus sequences for some of the upregulated genes. Finally, several of the novel transcripts identified in this study encode putative small peptides, which are biologically active.

Missing genes in the annotation of prokaryotic genomes

Background

Protein-coding gene detection in prokaryotic genomes is considered a much simpler problem than in intron-containing eukaryotic genomes. However there have been reports that prokaryotic gene finder programs have problems with small genes (either over-predicting or under-predicting). Therefore the question arises as to whether current genome annotations have systematically missing, small genes.

Results

We have developed a high-performance computing methodology to investigate this problem. In this methodology we compare all ORFs larger than or equal to 33 aa from all fully-sequenced prokaryotic replicons. Based on that comparison, and using conservative criteria requiring a minimum taxonomic diversity between conserved ORFs in different genomes, we have discovered 1,153 candidate genes that are missing from current genome annotations. These missing genes are similar only to each other and do not have any strong similarity to gene sequences in public databases, with the implication that these ORFs belong to missing gene families. We also uncovered 38,895 intergenic ORFs, readily identified as putative genes by similarity to currently annotated genes (we call these absent annotations). The vast majority of the missing genes found are small (less than 100 aa). A comparison of select examples with GeneMark, EasyGene and Glimmer predictions yields evidence that some of these genes are escaping detection by these programs.

Conclusions

Prokaryotic gene finders and prokaryotic genome annotations require improvement for accurate prediction of small genes. The number of missing gene families found is likely a lower bound on the actual number, due to the conservative criteria used to determine whether an ORF corresponds to a real gene.

Small RNAs and small proteins involved in resistance to cell envelope stress and acid shock in Escherichia coli: analysis of a bar-coded mutant collection

More than 80 small regulatory RNAs (sRNAs) and 60 proteins of 16 to 50 amino acids (small proteins) are encoded in the Escherichia coli genome. The vast majority of the corresponding genes have no known function. We screened 125 DNA bar-coded mutants to identify novel cell envelope stress and acute acid shock phenotypes associated with deletions of genes coding for sRNAs and small proteins. Nine deletion mutants (ssrA, micA, ybaM, ryeF, yqcG, sroH, ybhT, yobF, and glmY) were sensitive to cell envelope stress and two were resistant (rybB and blr). Deletion mutants of genes coding for four small proteins (yqgB, mgrB, yobF, and yceO) were sensitive to acute acid stress. We confirmed each of these phenotypes in one-on-one competition assays against otherwise-wild-type lacZ mutant cells. A more detailed investigation of the SsrA phenotype suggests that ribosome release is critical for resistance to cell envelope stress. The bar-coded deletion collection we generated can be screened for sensitivity or resistance to virtually any stress condition.