Prediction of transcription terminators in bacterial genomes

Detection of intrinsic transcription termination sites in bacteria: consensus from hairpin detection approaches

We compare the WebGeSTer and INtrinsic transcription TERmination hairPIN (INTERPIN) databases used for intrinsic transcription termination (ITT) site prediction in bacteria. The former deploys inverted nucleotide repeat detection for identification of RNA hairpin, while the latter a pair-potential function - the hairpin energy score evaluation being identical for both. We find INTERPIN more sensitive than WebGeSTer with about 6% and 51% additional predictions for ITTs in chromosomal and plasmid operons, respectively. INTERPIN hairpins are relatively shorter in length with ungapped stem, and even located in AT-rich segments, compared to GC-rich longer hairpins with a gapped stem in WebGeSTer. The GC%, length, and energy score from INTERPIN transcription units (TUs) are best inter-correlated while the lowest energy single hairpins from WebGeSTer, considered suitable for ITT, being the worst. Around 72% TUs from the two databases overlap, and ∼60% of all alternate ITT sites downstream of TUs overlap, of which 65% are cluster hairpins. This helps highlight hairpin features that can be used to identify termination sites in bacteria across different prediction methods. Overall, the pair-potential-function-based hairpins screened appear to be more consistent with the kinetic and thermodynamics processes of ITT known to date.Communicated by Ramaswamy H. Sarma.

Characterization and genome-informatic analysis of a novel lytic mendocina phage vB_PmeS_STP12 suitable for phage therapy pseudomonas or biocontrol

BACKGROUND

A novel lytic bacteriophage (phage) was isolated with Pseudomonas mendocina strain STP12 (P. mendocina) from the untreated site of Sewage Treatment Plant of Lovely Professional University, India. P. mendocina is a Gram-negative, rod-shaped, aerobic bacterium belonging to the family Pseudomonadaceae and has been reported in fifteen (15) cases of economically important diseases worldwide.

METHODS AND RESULTS

Here, a novel phage specifically infecting and killing P. mendocina strain STP12 was isolated from sewage sample using enrichment, spot test and double agar overlay (DAOL) method and was designated as vB_PmeS_STP12. The phage vB-PmeS-STP12 was viable at wide range of pH and temperature ranging from 4 to10 and - 20 to 70 °C respectively. Host range and efficiency of plating (EOP) analysis indicated that phage vB-PmeS-STP12 was capable of infecting and killing P. mendocina strain STP6 with EOP of 0.34. Phage vB_PmeS_STP12 was found to have a significant bacterial reduction (p < 0.005) at all the doses administered, particularly at optimal MOI of 1 PFU/CFU, compared to the control. Morphological analysis using high resolution transmission electron microscopy (HR-TEM) revealed an icosahedral capsid of ~ 55 nm in diameter on average with a short, non-contractile tail. The genome of vB_PmeS_STP12 is a linear, dsDNA containing 36,212 bp in size with a GC content of 58.87% harbouring 46 open reading frames (ORFs). The 46 predicted ORFs encode proteins with functional information categorized as lysis, replication, packaging, regulation, assembly, infection, immune, and hypothetical. However, the genome of vB_PmeS_STP12 appeared to be devoid of tRNAs, integrase gene, toxins genes, virulence factors, antimicrobial resistance genes (ARGs) and CRISPR arrays. The blast analysis with phylogeny revealed that vB_PmeS_STP12 is genetically similar to Pseudomonas phage PMBT14, Pseudomonas phage Almagne and Serratia phage Serbin with a highest identity of 74.00%, 74.93% and 59.48% respectively.

CONCLUSIONS

Taken together, characterization, morphological analysis and genome-informatics indicated that vB_PmeS_STP12 is podovirus morphotype belonging to the class Caudoviticetes, family Zobellviridae which appeared to be devoid of integrase gene, ARGs, CRISPR arrays, virulence factors and toxins genes, exhibiting stability and infectivity at wide range of pH (4 to10) and temperature (-20 to 70 °C), thereby making vB_PmeS_STP12 suitable for phage therapy or biocontrol. Based on the bibliometric analysis and data availability with respect to sequences deposited in GenBank, this is the first report of a phage infecting Pseudomonas mendocina.

Utilizing RNA-seq Data to Infer Bacterial Transcription Termination Sites and Validate Predictions

The transcription termination process is an important part of the gene expression process in the cell. It has been studied extensively, but many aspects of the mechanism are not well understood. The widespread availability of experimental RNA-seq data from high-throughput experiments provides a unique opportunity to infer the end of the transcription units genome wide. This data is available for both Rho-dependent and Rho-independent termination pathways that drive transcription termination in bacteria. Our book chapter gives an overview of the current knowledge of Rho-independent transcription termination mechanisms and the prediction approaches currently deployed to infer the termination sites. Thereafter, we describe our method that uses cluster hairpins to detect Rho-independent transcription termination sites. These clusters are a group of hairpins that lies at <15 bp from each other and are together capable of enforcing the termination process. The idea of a group of hairpins being extensively used for transcription termination is new, and results show that at least 52% of the total cases are of this type, while in the remaining cases, a single strong hairpin is capable of driving transcription termination. The reads derived from the RNA-seq data for corresponding bacteria have been used to validate the predicted sites. The predictions that match these RNA-seq derived sites have higher confidence, and we find almost 98% of the predicted sites, including alternate termination sites, to match the RNA-seq data. We discuss the features of predicted hairpins in detail for a better understanding of the Rho-independent transcription termination mechanism in bacteria. We also explain how users can use the tools developed by us to do transcription terminator predictions and design their experiments through genome-level visualization of the transcription termination sites from the precomputed INTERPIN database.

INTERPIN: A repository for intrinsic transcription termination hairpins in bacteria

The large-scale detection of putative intrinsic transcription terminators is limited to only a few bacteria currently. We discovered a group of hairpins, called cluster hairpins, present within 15 nucleotides from each other. These are expected to work in tandem to cause intrinsic transcription termination (ITT), while the single hairpin can do the same alone. Therefore, exploring these ITT sites and the hairpins across bacterial genomes becomes highly desirable. INTERPIN is the largest archived collection of in silico inferred ITT hairpins in bacteria, covering 12745 bacterial genomes and encompassing ten bacterial phyla for ∼25 million hairpins. Users can obtain details on operons, individual cluster, and single ITT hairpins that were screened therein. Integrated Genome Viewer (IGV) software interactively visualizes hairpin secondary and tertiary structures in the genomic context. We also discuss statistics for the occurrence of cluster or single hairpins and other termination alternatives while showing the validation of predicted hairpins against in vivo detected hairpins. The database is freely available at http://pallab.cds.iisc.ac.in/INTERPIN/. INTERPIN (database and software) can make predictions for both AT and GC-rich genomes, which has not been achieved by any other program so far. It can also be used to improve genome annotation as well as to get predictions to improve the understanding of the ITT pathway by further analysis.

Diversity of the type VI secretion systems in the Neisseria spp

Complete Type VI Secretion Systems were identified in the genome sequence data of Neisseria subflava isolates sourced from throat swabs of human volunteers. The previous report was the first to describe two complete Type VI Secretion Systems in these isolates, both of which were distinct in terms of their gene organization and sequence homology. Since publication of the first report, Type VI Secretion System subtypes have been identified in Neisseria spp. The characteristics of each type in N. subflava are further investigated here and in the context of the other Neisseria spp., including identification of the lineages containing the different types and subtypes. Type VI Secretion Systems use VgrG for delivery of toxin effector proteins; several copies of vgrG and associated effector / immunity pairs are present in Neisseria spp. Based on sequence similarity between strains and species, these core Type VI Secretion System genes, vgrG, and effector / immunity genes may diversify via horizontal gene transfer, an instrument for gene acquisition and repair in Neisseria spp.

Bioinformatic and experimental methods to identify and validate bacterial RNA-human RNA interactions

Ample evidence supports the importance of the microbiota on human health and disease. Recent studies suggest that extracellular vesicles are an important means of bacterial-host communication, in part via the transport of small RNAs (sRNAs). Bacterial sRNAs have been shown to co-precipitate with human and mouse RNA-induced silencing complex, hinting that some may regulate gene expression as eukaryotic microRNAs do. Bioinformatic tools, including those that can incorporate an sRNA's secondary structure, can be used to predict interactions between bacterial sRNAs and human messenger RNAs (mRNAs). Validation of these potential interactions using reproducible experimental methods is essential to move the field forward. This review will cover the evidence of interspecies communication via sRNAs, bioinformatic tools currently available to identify potential bacterial sRNA-host (specifically, human) mRNA interactions, and experimental methods to identify and validate those interactions.

A candidate multi-epitope vaccine against Lumpy skin disease

Lumpy skin disease (LSD) is a fulminant infectious disease that mostly affects cattle and causes considerable economic loss throughout the globe. This study was conducted to develop a new multi-epitope-based vaccine against LSD that can elicit immunological responses using an in silico reverse vaccinology approach. Initially, three antigenic proteins, protein E5, E3 ubiquitin-protein ligase LAP and 62 kDa protein, were manipulated to recognize potential T-cell and B-cell epitopes. To identify superior epitopes, a variety of bioinformatic techniques including antigenicity testing, transmembrane topology screening, allergenicity assessment, conservancy analysis, and toxicity evaluation were used. Finally, three new subunit vaccines (construct V1, V2 and V3) were developed employing the most effective epitopes, suitable adjuvants, pan HLA DR-binding epitope (PADRE) and linkers. Then, based on the antigenicity, solubility, and validation score of the 3D structures, construct V2 was chosen as one of the best candidate vaccines. The results of the molecular dynamic simulation and disulphide engineering indicated that the vaccine (construct V2) was stable. Additionally, the immunological simulation findings supported the vaccine candidate's ability to trigger humoral and cellular immune responses. Further validation of the proposed vaccine candidate may necessitate additional in vitro and in vivo investigations.

Inverse folding based pre-training for the reliable identification of intrinsic transcription terminators

It is well-established that neural networks can predict or identify structural motifs of non-coding RNAs (ncRNAs). Yet, the neural network based identification of RNA structural motifs is limited by the availability of training data that are often insufficient for learning features of specific ncRNA families or structural motifs. Aiming to reliably identify intrinsic transcription terminators in bacteria, we introduce a novel pre-training approach that uses inverse folding to generate training data for predicting or identifying a specific family or structural motif of ncRNA. We assess the ability of neural networks to identify secondary structure by systematic in silico mutagenesis experiments. In a study to identify intrinsic transcription terminators as functionally well-understood RNA structural motifs, our inverse folding based pre-training approach significantly boosts the performance of neural network topologies, which outperform previous approaches to identify intrinsic transcription terminators. Inverse-folding based pre-training provides a simple, yet highly effective way to integrate the well-established thermodynamic energy model into deep neural networks for identifying ncRNA families or motifs. The pre-training technique is broadly applicable to a range of network topologies as well as different types of ncRNA families and motifs.

Intrinsic transcriptional terminators are essential regulators in determining the 3’ end of transcripts in bacteria. The underlying mechanism involves RNA secondary structure, where nucleotides fold into a specific hairpin motif. Identifying terminator sequences in bacterial genomes has conventionally been approached with well-established energy models for structural motifs. However, the folding mechanism of transcription terminators is understood only partially, limiting the success of energy-model based identification.

Neural networks have been proposed to overcome these limitations. However, their adoption for predicting and identifying RNA secondary structure has been a double edged sword: Neural networks promise to learn features that are not represented by the energy models, while they are black boxes that lack explicit modeling assumptions and may fail to account for features that are well understandable based on decades-old energy models.

Here, we introduce a pre-training approach for neural networks that uses energy-model based inverse folding of structural motifs. As we demonstrate, this approach “brings back the energy model” to identify transcriptional terminators and overcomes the limitations of previous energy-model based predictions. Our approach works for diverse types of neural networks, and is suitable for the identification of structural motifs of many other RNA molecules beyond transcriptional terminators.

Reassembling a cannon in the DNA defense arsenal: Genetics of StySA, a BREX phage exclusion system in Salmonella lab strains

Understanding mechanisms that shape horizontal exchange in prokaryotes is a key problem in biology. A major limit on DNA entry is imposed by restriction-modification (RM) processes that depend on the pattern of DNA modification at host-specified sites. In classical RM, endonucleolytic DNA cleavage follows detection of unprotected sites on entering DNA. Recent investigation has uncovered BREX (BacteRiophage EXclusion) systems. These RM-like activities employ host protection by DNA modification, but immediate replication arrest occurs without evident of nuclease action on unmodified phage DNA. Here we show that the historical stySA RM locus of Salmonella enterica sv Typhimurium is a variant BREX system. A laboratory strain disabled for both the restriction and methylation activity of StySA nevertheless has wild type sequence in pglX, the modification gene homolog. Instead, flanking genes pglZ and brxC each carry multiple mutations (μ) in their C-terminal domains. We further investigate this system in situ, replacing the mutated pglZμ and brxCμ genes with the WT counterpart. PglZ-WT supports methylation in the presence of either BrxCμ or BrxC-WT but not in the presence of a deletion/insertion allele, ΔbrxC::cat. Restriction requires both BrxC-WT and PglZ-WT, implicating the BrxC C-terminus specifically in restriction activity. These results suggests that while BrxC, PglZ and PglX are principal components of the BREX modification activity, BrxL is required for restriction only. Furthermore, we show that a partial disruption of brxL disrupts transcription globally.

Horizontal gene transfer is a major driver of evolution and adaptation in bacteria. Genes from outside may be beneficial or dangerous to the receiving cell. Benefits include new food sources such as sugars, or new homes by adhesion, or new resistances, as to antibiotics. Dangers are posed by bacteriophages--viruses that take over the cell machinery, multiply, and release progeny to kill sister cells. Host-dependent restriction-modification systems enable defense that distinguishes relatives from strangers: using a modification pattern (M) carried by DNA bases added by the host cell to prevent restriction (R). Sisters and cousin cells will have the same protective pattern on DNA, while DNA of foreign origin will have the wrong M pattern and be restricted (R, rejected). Typically, restriction involves nuclease digestion. Here we address the enigmatic StySA RM system, one of the earliest to be genetically characterized. It is a variant of the newly recognized defense mechanism, BREX. BREX systems also track DNA history via modification pattern, but restrict by a novel, uncharacterized mechanism. Like other BREX family systems, StySA-BREX modification requires multiple components. When StySA-BREX transcription is unbalanced, we find global disruption of gene transcription. The disruption pattern does not suggest SOS-inducing damage to DNA.

Synthetic 3'-UTR valves for optimal metabolic flux control in Escherichia coli

As the design of genetic circuitry for synthetic biology becomes more sophisticated, diverse regulatory bioparts are required. Despite their importance, well-characterized 3′-untranslated region (3′-UTR) bioparts are limited. Thus, transcript 3′-ends require further investigation to understand the underlying regulatory role and applications of the 3′-UTR. Here, we revisited the use of Term-Seq in the Escherichia coli strain K-12 MG1655 to enhance our understanding of 3′-UTR regulatory functions and to provide a diverse collection of tunable 3′-UTR bioparts with a wide termination strength range. Comprehensive analysis of 1,629 transcript 3′-end positions revealed multiple 3′-termini classes generated through transcription termination and RNA processing. The examination of individual Rho-independent terminators revealed a reduction in downstream gene expression over a wide range, which led to the design of novel synthetic metabolic valves that control metabolic fluxes in branched pathways. These synthetic metabolic valves determine the optimal balance of heterologous pathways for maximum target biochemical productivity. The regulatory strategy using 3′-UTR bioparts is advantageous over promoter- or 5′-UTR-based transcriptional control as it modulates gene expression at transcription levels without trans-acting element requirements (e.g. transcription factors). Our results provide a foundational platform for 3′-UTR engineering in synthetic biology applications.

Clusters of hairpins induce intrinsic transcription termination in bacteria

Intrinsic transcription termination (ITT) sites are currently identified by locating single and double-adjacent RNA hairpins downstream of the stop codon. ITTs for a limited number of genes/operons in only a few bacterial genomes are currently known. This lack of coverage is a lacuna in the existing ITT inference methods. We have studied the inter-operon regions of 13 genomes covering all major phyla in bacteria, for which good quality public RNA-seq data exist. We identify ITT sites in 87% of cases by predicting hairpin(s) and validate against 81% of cases for which the RNA-seq derived sites could be calculated. We identify 72% of these sites correctly, with 98% of them located ≤ 80 bases downstream of the stop codon. The predicted hairpins form a cluster (when present < 15 bases) in two-thirds of the cases, the remaining being single hairpins. The largest number of clusters is formed by two hairpins, and the occurrence decreases exponentially with an increasing number of hairpins in the cluster. Our study reveals that hairpins form an effective ITT unit when they act in concert in a cluster. Their pervasiveness along with single hairpin terminators corroborates a wider utilization of ITT mechanisms for transcription control across bacteria.

iTerm-PseKNC: a sequence-based tool for predicting bacterial transcriptional terminators

MOTIVATION

Transcription termination is an important regulatory step of gene expression. If there is no terminator in gene, transcription could not stop, which will result in abnormal gene expression. Detecting such terminators can determine the operon structure in bacterial organisms and improve genome annotation. Thus, accurate identification of transcriptional terminators is essential and extremely important in the research of transcription regulations.

RESULTS

In this study, we developed a new predictor called 'iTerm-PseKNC' based on support vector machine to identify transcription terminators. The binomial distribution approach was used to pick out the optimal feature subset derived from pseudo k-tuple nucleotide composition (PseKNC). The 5-fold cross-validation test results showed that our proposed method achieved an accuracy of 95%. To further evaluate the generalization ability of 'iTerm-PseKNC', the model was examined on independent datasets which are experimentally confirmed Rho-independent terminators in Escherichia coli and Bacillus subtilis genomes. As a result, all the terminators in E. coli and 87.5% of the terminators in B. subtilis were correctly identified, suggesting that the proposed model could become a powerful tool for bacterial terminator recognition.

AVAILABILITY AND IMPLEMENTATION

For the convenience of most of wet-experimental researchers, the web-server for 'iTerm-PseKNC' was established at http://lin-group.cn/server/iTerm-PseKNC/, by which users can easily obtain their desired result without the need to go through the detailed mathematical equations involved.

Revisiting intrinsic transcription termination in mycobacteria: U-tract downstream of secondary structure is dispensable for termination

In Escherichia coli, G/C-rich hairpin structure followed by a U-tract in the 3' region of the nascent RNA are crucial determinants for intrinsic or factor independent transcription termination. In mycobacteria, there is a scarcity of such intrinsic terminators. However, secondary structures having G/C-rich stem devoid of any U's or with suboptimal U-tracts were identified earlier as terminators and found to be functional both in vitro and in vivo. Two different observations - that a mycobacterial RNA polymerase (RNAP) does not function at intrinsic terminators devoid of U-tracts and the identification of an altogether new motif for termination in mycobacteria necessitated re-examining a number of putative terminators for their function as terminators. When these in silico identified non-canonical terminators were subjected to experimental validation, they were found to dissociate RNA from the elongating RNAP. Termination is observed when the U-tracts were reduced, or totally absent both in vitro and in vivo. Our results, thus indicate that the presence of U-tract following the G/C-rich stem in an intrinsic terminator may not be an essential determinant for transcription termination in mycobacteria.

The Rho-Independent Transcription Terminator for the porA Gene Enhances Expression of the Major Outer Membrane Protein and Campylobacter jejuni Virulence in Abortion Induction

Campylobacter jejuni is a leading cause of foodborne illnesses worldwide. Its porA gene encodes the major outer membrane protein (MOMP) that is abundantly expressed and has important physiological functions, including a key role in systemic infection and abortion induction in pregnant animals. Despite the importance of porA in C. jejuni pathogenesis, mechanisms modulating its expression levels remain elusive. At the 3' end of the porA transcript, there is a Rho-independent transcription terminator (named T _porA in this study). Whether T _porA affects the expression and function of MOMP remains unknown and is investigated in this study. Green fluorescent protein (GFP) fusion constructs with the porA promoter at the 5' end and an intact T _porA or no T _porA at the 3' end of the gfp coding sequence revealed that both the transcript level of gfp and its fluorescence signals were more than 2-fold higher in the construct with T _porA than in the one without T _porA Real-time quantitative PCR (qRT-PCR) analysis of the porA mRNA and immunoblot detection of MOMP in C. jejuni showed that disruption of T _porA significantly reduced the porA transcript level and the expression of MOMP. An mRNA decay assay demonstrated that disruption of T _porA resulted in a shortened transcript half-life of the upstream gfp or porA gene, indicating that T _porA enhances mRNA stability. In the guinea pig model, the C. jejuni construct with an interrupted T _porA was significantly attenuated in abortion induction. Together, these results indicate that T _porA enhances the expression level of MOMP by stabilizing its mRNA and influences the virulence of C. jejuni.

Expanding the Toolbox of Broad Host-Range Transcriptional Terminators for Proteobacteria through Metagenomics

As the field of synthetic biology moves toward the utilization of novel bacterial chassis, there is a growing need for biological parts with enhanced performance in a wide number of hosts. Is not unusual that biological parts (such as promoters and terminators), initially characterized in the model bacterium Escherichia coli, do not perform well when implemented in alternative hosts, such as Pseudomonas, therefore limiting the construction of synthetic circuits in industrially relevant bacteria, for instance Pseudomonas putida. In order to address this limitation, we present here the mining of transcriptional terminators through functional metagenomics to identify novel parts with broad host-range activity. Using a GFP-based terminator trap strategy and a broad host-range plasmid, we identified 20 clones with potential terminator activity in P. putida. Further characterization allowed the identification of 4 unique sequences ranging from 58 to 181 bp long that efficiently terminate transcription in P. putida, E. coli, Burkholderia phymatum, and two Pseudomonas strains isolated from Antarctica. Therefore, this work presents a new set of biological parts useful for the engineering of synthetic circuits in Proteobacteria.

Bacterial Small RNAs in the Genus Herbaspirillum spp

The genus Herbaspirillum includes several strains isolated from different grasses. The identification of non-coding RNAs (ncRNAs) in the genus Herbaspirillum is an important stage studying the interaction of these molecules and the way they modulate physiological responses of different mechanisms, through RNA⁻RNA interaction or RNA⁻protein interaction. This interaction with their target occurs through the perfect pairing of short sequences (cis-encoded ncRNAs) or by the partial pairing of short sequences (trans-encoded ncRNAs). However, the companion Hfq can stabilize interactions in the trans-acting class. In addition, there are Riboswitches, located at the 5' end of mRNA and less often at the 3' end, which respond to environmental signals, high temperatures, or small binder molecules. Recently, CRISPR (clustered regularly interspaced palindromic repeats), in prokaryotes, have been described that consist of serial repeats of base sequences (spacer DNA) resulting from a previous exposure to exogenous plasmids or bacteriophages. We identified 285 ncRNAs in Herbaspirillum seropedicae (H. seropedicae) SmR1, expressed in different experimental conditions of RNA-seq material, classified as cis-encoded ncRNAs or trans-encoded ncRNAs and detected RNA riboswitch domains and CRISPR sequences. The results provide a better understanding of the participation of this type of RNA in the regulation of the metabolism of bacteria of the genus Herbaspirillum spp.

In silico Analysis Suggests Common Appearance of scaRNAs in Type II Systems and Their Association With Bacterial Virulence

In addition to its well-established defense function, CRISPR/Cas can also exhibit crucial non-canonical activity through endogenous gene expression regulation, which was found to mainly affect bacterial virulence. These non-canonical functions depend on scaRNA, which is a small RNA encoded outside of CRISPR array, that is typically flanked by a transcription start site (TSS) and a terminator, and is in part complementary to another small CRISPR/Cas-associated RNA (tracrRNAs). Identification of scaRNAs is however largely complicated by the scarcity of RNA-Seq data across different bacteria, so that they were identified only in a relatively rare CRISPR/Cas subtype (IIB), and the possibility of finding them in other Type II systems is currently unclear. This study presents the first effort toward systematic detection of small CRISPR/Cas-associated regulatory RNAs, where obtained predictions can guide future experiments. The core of our approach is ab initio detection of small RNAs from bacterial genome, which is based on jointly predicting transcription signals - TSS and terminators - and homology to CRISPR array repeat. Particularly, we employ our improved approach for detecting bacterial TSS, since accurate TSS detection is the main limiting factor for accurate small RNA prediction. We also explore how our predictions match to available RNA-Seq data and analyze their conservation across related bacterial species. In Type IIB systems, our predictions are consistent with experimental data, and we systematically identify scaRNAs throughout this subtype. Furthermore, we identify scaRNA:tracrRNA pairs in a number of IIA/IIC systems, where the appearance of scaRNAs co-occurs with the strains being pathogenic. RNA-Seq and conservation analysis show that our method is well suited for predicting CRISPR/Cas-associated small RNAs. We also find possible existence of a modified mechanism of CRISPR-associated small RNA action, which, interestingly, closely resembles the setup employed in biotechnological applications. Overall, our findings indicate that scaRNA:tracrRNA pairs are present in all subtypes of Type II systems, and point to an underlying connection with bacterial virulence. In addition to formulating these hypotheses, careful manual curation that we performed, makes an important first step toward fully automated predictor of CRISPR/Cas-associated small RNAs, which will allow their large scale analysis across diverse bacterial genomes.

A novel whole-cell biosensor of Pseudomonas aeruginosa to monitor the expression of quorum sensing genes

A novel whole-cell biosensor was developed to noninvasively and simultaneously monitor the in situ genetic activities of the four quorum sensing (QS) networks in Pseudomonas aeruginosa PAO1, including the las, rhl, pqs, and iqs systems. P. aeruginosa PAO1 is a model bacterium for studies of biofilm and pathogenesis while both processes are closely controlled by the QS systems. This biosensor worked well by selectively monitoring the expression of one representative gene from each network. In the biosensor, the promoter regions of lasI, rhlI, pqsA, and ambB (QS genes) controlled the fluorescent reporter genes of Turbo YFP, mTag BFP2, mNEON Green, and E2-Orange, respectively. The biosensor was successful in monitoring the impact of an important environmental factor, salt stress, on the genetic regulation of QS networks. High salt concentrations (≥ 20 g·L^-1) significantly downregulated rhlI, pqsA, and ambB after the biosensor was incubated for 17 h to 18 h at 37 °C, resulting in slow bacterial growth.

Increased leukotoxin production: Characterization of 100 base pairs within the 530 base pair leukotoxin promoter region of Aggregatibacter actinomycetemcomitans

Aggregatibacter actinomycetemcomitans leukotoxin (LtxA) is a major virulence factor that kills leukocytes permitting it’s escape from host immune surveillance. A. actinomycetemcomitans strains can produce high or low levels of toxin. Genetic differences reside in the “so called JP2” ltxA promoter region. These hyper-leukotoxin producing strains with the 530 bp deletion have been studied in detail. However, regions contained within the 530 bp deletion that could be responsible for modulation of leukotoxin production have not been defined. Here, we report, for the first time, on regions within the 530 bp that are responsible for high-levels of ltxA expression. We constructed a deletion of 530 bps in a primate isolate of A. actinomycetemcomitans, which produced leukotoxin equivalent to the JP2 strain. We then constructed sequential deletions in regions that span the 530 bps. Results indicated that expression of the ltxA transcript was reduced by a potential transcriptional terminator in promoter region 298 to 397 with a ΔG = −7.9 kcal/mol. We also confirmed previous findings that transcriptional fusion between the orfX region and ltxC increased ltxA expression. In conclusion, we constructed a hyper-leukotoxin producing A. actinomycetemcomitans strain and identified a terminator located in the promoter region extending from 298–397 that alters ltxA expression.

A Putative O-Linked β-N-Acetylglucosamine Transferase Is Essential for Hormogonium Development and Motility in the Filamentous Cyanobacterium Nostoc punctiforme

Most species of filamentous cyanobacteria are capable of gliding motility, likely via a conserved type IV pilus-like system that may also secrete a motility-associated polysaccharide. In a subset of these organisms, motility is achieved only after the transient differentiation of hormogonia, which are specialized filaments that enter a nongrowth state dedicated to motility. Despite the fundamental importance of hormogonia to the life cycles of many filamentous cyanobacteria, the molecular regulation of hormogonium development is largely undefined. To systematically identify genes essential for hormogonium development and motility in the model heterocyst-forming filamentous cyanobacterium Nostoc punctiforme, a forward genetic screen was employed. The first gene identified using this screen, designated ogtA, encodes a putative O-linked β-N-acetylglucosamine transferase (OGT). The deletion of ogtA abolished motility, while ectopic expression of ogtA induced hormogonium development even under hormogonium-repressing conditions. Transcription of ogtA is rapidly upregulated (1 h) following hormogonium induction, and an OgtA-GFPuv fusion protein localized to the cytoplasm. In developing hormogonia, accumulation of PilA but not HmpD is dependent on ogtA Reverse transcription-quantitative PCR (RT-qPCR) analysis indicated equivalent levels of pilA transcript in the wild-type and ΔogtA mutant strains, while a reporter construct consisting of the intergenic region in the 5' direction of pilA fused to gfp produced lower levels of fluorescence in the ΔogtA mutant strain than in the wild type. The production of hormogonium polysaccharide in the ΔogtA mutant strain is reduced compared to that in the wild type but comparable to that in a pilA deletion strain. Collectively, these results imply that O-GlcNAc protein modification regulates the accumulation of PilA via a posttranscriptional mechanism in developing hormogonia.IMPORTANCE Filamentous cyanobacteria are among the most developmentally complex prokaryotes. Species such as Nostoc punctiforme develop an array of cell types, including nitrogen-fixing heterocysts, spore-like akinetes, and motile hormogonia, that function in dispersal as well as the establishment of nitrogen-fixing symbioses with plants and fungi. These symbioses are major contributors to global nitrogen fixation. Despite the fundamental importance of hormogonia to the life cycle of filamentous cyanobacteria and the establishment of symbioses, the molecular regulation of hormogonium development is largely undefined. We employed a genetic screen to identify genes essential for hormogonium development and motility in Nostoc punctiforme The first gene identified using this screen encodes a eukaryotic-like O-linked β-N-acetylglucosamine transferase that is required for accumulation of PilA in hormogonia.

Synbiological systems for complex natural products biosynthesis

Natural products (NPs) continue to play a pivotal role in drug discovery programs. The rapid development of synthetic biology has conferred the strategies of NPs production. Synthetic biology is a new engineering discipline that aims to produce desirable products by rationally programming the biological parts and manipulating the pathways. However, there is still a challenge for integrating a heterologous pathway in chassis cells for overproduction purpose due to the limited characterized parts, modules incompatibility, and cell tolerance towards product. Enormous endeavors have been taken for mentioned issues. Herein, in this review, the progresses in naturally discovering novel biological parts and rational design of synthetic biological parts are reviewed, combining with the advanced assembly technologies, pathway engineering, and pathway optimization in global network guidance. The future perspectives are also presented.

DNA uptake sequences in Neisseria gonorrhoeae as intrinsic transcriptional terminators and markers of horizontal gene transfer

DNA uptake sequences are widespread throughout the Neisseria gonorrhoeae genome. These short, conserved sequences facilitate the exchange of endogenous DNA between members of the genus Neisseria. Often the DNA uptake sequences are present as inverted repeats that are able to form hairpin structures. It has been suggested previously that DNA uptake sequence inverted repeats present 3' of genes play a role in rho-independent termination and attenuation. However, there is conflicting experimental evidence to support this role. The aim of this study was to determine the role of DNA uptake sequences in transcriptional termination. Both bioinformatics predictions, conducted using TransTermHP, and experimental evidence, from RNA-seq data, were used to determine which inverted repeat DNA uptake sequences are transcriptional terminators and in which direction. Here we show that DNA uptake sequences in the inverted repeat configuration occur in N. gonorrhoeae both where the DNA uptake sequence precedes the inverted version of the sequence and also, albeit less frequently, in reverse order. Due to their symmetrical configuration, inverted repeat DNA uptake sequences can potentially act as bi-directional terminators, therefore affecting transcription on both DNA strands. This work also provides evidence that gaps in DNA uptake sequence density in the gonococcal genome coincide with areas of DNA that are foreign in origin, such as prophage. This study differentiates for the first time, to our knowledge, between DNA uptake sequences that form intrinsic transcriptional terminators and those that do not, providing characteristic features within the flanking inverted repeat that can be identified.

Genomics of Three New Bacteriophages Useful in the Biocontrol of Salmonella

Non-typhoid Salmonella is the principal pathogen related to food-borne diseases throughout the world. Widespread antibiotic resistance has adversely affected human health and has encouraged the search for alternative antimicrobial agents. The advances in bacteriophage therapy highlight their use in controlling a broad spectrum of food-borne pathogens. One requirement for the use of bacteriophages as antibacterials is the characterization of their genomes. In this work, complete genome sequencing and molecular analyses were carried out for three new virulent Salmonella-specific bacteriophages (UAB_Phi20, UAB_Phi78, and UAB_Phi87) able to infect a broad range of Salmonella strains. Sequence analysis of the genomes of UAB_Phi20, UAB_Phi78, and UAB_Phi87 bacteriophages did not evidence the presence of known virulence-associated and antibiotic resistance genes, and potential immunoreactive food allergens. The UAB_Phi20 genome comprised 41,809 base pairs with 80 open reading frames (ORFs); 24 of them with assigned function. Genome sequence showed a high homology of UAB_Phi20 with Salmonella bacteriophage P22 and other P22likeviruses genus of the Podoviridae family, including ST64T and ST104. The DNA of UAB_Phi78 contained 44,110 bp including direct terminal repeats (DTR) of 179 bp and 58 putative ORFs were predicted and 20 were assigned function. This bacteriophage was assigned to the SP6likeviruses genus of the Podoviridae family based on its high similarity not only with SP6 but also with the K1-5, K1E, and K1F bacteriophages, all of which infect Escherichia coli. The UAB_Phi87 genome sequence consisted of 87,669 bp with terminal direct repeats of 608 bp; although 148 ORFs were identified, putative functions could be assigned to only 29 of them. Sequence comparisons revealed the mosaic structure of UAB_Phi87 and its high similarity with bacteriophages Felix O1 and wV8 of E. coli with respect to genetic content and functional organization. Phylogenetic analysis of large terminase subunits confirms their packaging strategies and grouping to the different phage genus type. All these studies are necessary for the development and the use of an efficient cocktail with commercial applications in bacteriophage therapy against Salmonella.

Analysis of the First Temperate Broad Host Range Brucellaphage (BiPBO1) Isolated from B. inopinata

Brucella species are important human and animal pathogens. Though, only little is known about mobile genetic elements of these highly pathogenic bacteria. To date, neither plasmids nor temperate phages have been described in brucellae. We analyzed genomic sequences of various reference and type strains and identified a number of putative prophages residing within the Brucella chromosomes. By induction, phage BiPBO1 was isolated from Brucella inopinata. BiPBO1 is a siphovirus that infects several Brucella species including Brucella abortus and Brucella melitensis. Integration of the phage genome occurs adjacent to a tRNA gene in chromosome 1 (chr 1). The bacterial (attB) and phage (attP) attachment sites comprise an identical sequence of 46 bp. This sequence exists in many Brucella and Ochrobactrum species. The BiPBO1 genome is composed of a 46,877 bp double-stranded DNA. Eighty-seven putative gene products were determined, of which 32 could be functionally assigned. Strongest similarities were found to a temperate phage residing in the chromosome of Ochrobactrum anthropi ATCC 49188 and to prophages identified in several families belonging to the order rhizobiales. The data suggest that horizontal gene transfer may occur between Brucella and Ochrobactrum and underpin the close relationship of these environmental and pathogenic bacteria.

Effects of cooperation between translating ribosome and RNA polymerase on termination efficiency of the Rho-independent terminator

An experimental system was designed to measure in vivo termination efficiency (TE) of the Rho-independent terminator and position–function relations were quantified for the terminator tR2 in Escherichia coli. The terminator function was almost completely repressed when tR2 was located several base pairs downstream from the gene, and TE gradually increased to maximum values with the increasing distance between the gene and terminator. This TE–distance relation reflected a stochastic coupling of the ribosome and RNA polymerase (RNAP). Terminators located in the first 100 bp of the coding region can function efficiently. However, functional repression was observed when the terminator was located in the latter part of the coding region, and the degree of repression was determined by transcriptional and translational dynamics. These results may help to elucidate mechanisms of Rho-independent termination and reveal genomic locations of terminators and functions of the sequence that precedes terminators. These observations may have important applications in synthetic biology.

Revisiting operons: an analysis of the landscape of transcriptional units in E. coli

Background

Bacterial operons are considerably more complex than what were thought. At least their components are dynamically rather than statically defined as previously assumed. Here we present a computational study of the landscape of the transcriptional units (TUs) of E. coli K12, revealed by the available genomic and transcriptomic data, providing new understanding about the complexity of TUs as a whole encoded in the genome of E. coli K12.

Results and conclusion

Our main findings include that (i) different TUs may overlap with each other by sharing common genes, giving rise to clusters of overlapped TUs (TUCs) along the genomic sequence; (ii) the intergenic regions in front of the first gene of each TU tend to have more conserved sequence motifs than those of the other genes inside the TU, suggesting that TUs each have their own promoters; (iii) the terminators associated with the 3’ ends of TUCs tend to be Rho-independent terminators, substantially more often than terminators of TUs that end inside a TUC; and (iv) the functional relatedness of adjacent gene pairs in individual TUs is higher than those in TUCs, suggesting that individual TUs are more basic functional units than TUCs.

Electronic supplementary material

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

Genome-wide transcription start site profiling in biofilm-grown Burkholderia cenocepacia J2315

BACKGROUND

Burkholderia cenocepacia is a soil-dwelling Gram-negative Betaproteobacterium with an important role as opportunistic pathogen in humans. Infections with B. cenocepacia are very difficult to treat due to their high intrinsic resistance to most antibiotics. Biofilm formation further adds to their antibiotic resistance. B. cenocepacia harbours a large, multi-replicon genome with a high GC-content, the reference genome of strain J2315 includes 7374 annotated genes. This study aims to annotate transcription start sites and identify novel transcripts on a whole genome scale.

METHODS

RNA extracted from B. cenocepacia J2315 biofilms was analysed by differential RNA-sequencing and the resulting dataset compared to data derived from conventional, global RNA-sequencing. Transcription start sites were annotated and further analysed according to their position relative to annotated genes.

RESULTS

Four thousand ten transcription start sites were mapped over the whole B. cenocepacia genome and the primary transcription start site of 2089 genes expressed in B. cenocepacia biofilms were defined. For 64 genes a start codon alternative to the annotated one was proposed. Substantial antisense transcription for 105 genes and two novel protein coding sequences were identified. The distribution of internal transcription start sites can be used to identify genomic islands in B. cenocepacia. A potassium pump strongly induced only under biofilm conditions was found and 15 non-coding small RNAs highly expressed in biofilms were discovered.

CONCLUSIONS

Mapping transcription start sites across the B. cenocepacia genome added relevant information to the J2315 annotation. Genes and novel regulatory RNAs putatively involved in B. cenocepacia biofilm formation were identified. These findings will help in understanding regulation of B. cenocepacia biofilm formation.

Bioinformatics as a first-line approach for understanding bacteriophage transcription

Current approach to understanding bacteriophage transcription strategies during infection includes a combination of experimental and bioinformatics approaches, which is often time and resource consuming. Given the exponentially growing number of sequenced bacteriophage genomes, it becomes sensible asking to what extent one can understand bacteriophage transcription by using bioinformatics methods alone. We here argue that a suitable choice of computational methods may provide a highly efficient first-line approach for underst-anding bacteriophage transcription.

Campylobacter group II phage CP21 is the prototype of a new subgroup revealing a distinct modular genome organization and host specificity

BACKGROUND

The application of phages is a promising tool to reduce the number of Campylobacter along the food chain. Besides the efficacy against a broad range of strains, phages have to be safe in terms of their genomes. Thus far, no genes with pathogenic potential (e.g., genes encoding virulence factors) have been detected in Campylobacter phages. However, preliminary studies suggested that the genomes of group II phages may be diverse and prone to genomic rearrangements.

RESULTS

We determined and analysed the genomic sequence (182,761 bp) of group II phage CP21 that is closely related to the already characterized group II phages CP220 and CPt10. The genomes of these phages are comprised of four modules separated by very similar repeat regions, some of which harbouring open reading frames (ORFs). Though, the arrangement of the modules and the location of some ORFs on the genomes are different in CP21 and in CP220/CPt10. In this work, a PCR system was established to study the modular genome organization of other group II phages demonstrating that they belong to different subgroups of the CP220-like virus genus, the prototypes of which are CP21 and CP220. The subgroups revealed different restriction patterns and, interestingly enough, also distinct host specificities, tail fiber proteins and tRNA genes. We additionally analysed the genome of group II phage vB_CcoM-IBB_35 (IBB_35) for which to date only five individual contigs could be determined. We show that the contigs represent modules linked by long repeat regions enclosing some yet not identified ORFs (e.g., for a head completion protein). The data suggest that IBB_35 is a member of the CP220 subgroup.

CONCLUSION

Campylobacter group II phages are diverse regarding their genome organization. Since all hitherto characterized group II phages contain numerous genes for transposases and homing endonucleases as well as similar repeat regions, it cannot be excluded that these phages are genetically unstable. To answer this question, further experiments and sequencing of more group II phages should be performed.

A Search for Ribonucleic Antiterminator Sites in Bacterial Genomes: Not Only Antitermination?

BglG/LicT-like proteins are transcriptional antiterminators that prevent termination of transcription at intrinsic terminators by binding to ribonucleic antiterminator (RAT) sites and stabilizing an RNA conformation which is mutually exclusive with the terminator structure. The known RAT sites, which are located in intergenic regions of sugar utilization operons, show low sequence conservation but significant structural analogy. To assess the prevalence of RATs in bacterial genomes, we employed bioinformatic tools that describe RNA motifs based on both sequence and structural constraints. Using descriptors with different stringency, we searched the genomes of Escherichiacoli K12, uropathogenic E. coli and Bacillus subtilis for putative RATs. Our search identified all known RATs and additional putative RAT elements. Surprisingly, most putative RATs do not overlap an intrinsic terminator and many reside within open reading frames (ORFs). The ability of one of the putative RATs, which is located within an antiterminator-encoding ORF and does not overlap a terminator, to bind to its cognate antiterminator protein in vitro and in vivo was confirmed experimentally. Our results suggest that the capacity of RAT elements has been exploited during evolution to mediate activities other than antitermination, for example control of transcription elongation or of RNA stability.

Genome-Wide Transcriptional Profiling of Clostridium perfringens SM101 during Sporulation Extends the Core of Putative Sporulation Genes and Genes Determining Spore Properties and Germination Characteristics

The formation of bacterial spores is a highly regulated process and the ultimate properties of the spores are determined during sporulation and subsequent maturation. A wide variety of genes that are expressed during sporulation determine spore properties such as resistance to heat and other adverse environmental conditions, dormancy and germination responses. In this study we characterized the sporulation phases of C. perfringens enterotoxic strain SM101 based on morphological characteristics, biomass accumulation (OD600), the total viable counts of cells plus spores, the viable count of heat resistant spores alone, the pH of the supernatant, enterotoxin production and dipicolinic acid accumulation. Subsequently, whole-genome expression profiling during key phases of the sporulation process was performed using DNA microarrays, and genes were clustered based on their time-course expression profiles during sporulation. The majority of previously characterized C. perfringens germination genes showed upregulated expression profiles in time during sporulation and belonged to two main clusters of genes. These clusters with up-regulated genes contained a large number of C. perfringens genes which are homologs of Bacillus genes with roles in sporulation and germination; this study therefore suggests that those homologs are functional in C. perfringens. A comprehensive homology search revealed that approximately half of the upregulated genes in the two clusters are conserved within a broad range of sporeforming Firmicutes. Another 30% of upregulated genes in the two clusters were found only in Clostridium species, while the remaining 20% appeared to be specific for C. perfringens. These newly identified genes may add to the repertoire of genes with roles in sporulation and determining spore properties including germination behavior. Their exact roles remain to be elucidated in future studies.

Characterization of CPAF critical residues and secretion during Chlamydia trachomatis infection

CPAF (chlamydial protease-like activity factor), a Chlamydia serine protease, is activated via proximity-induced intermolecular dimerization that triggers processing and removal of an inhibitory peptide occupying the CPAF substrate-binding groove. An active CPAF is a homodimer of two identical intramolecular heterodimers, each consisting of 29-kDa N-terminal and 35-kDa C-terminal fragments. However, critical residues for CPAF intermolecular dimerization, catalytic activity, and processing were defined in cell-free systems. Complementation of a CPAF-deficient chlamydial organism with a plasmid-encoded CPAF has enabled us to characterize CPAF during infection. The transformants expressing CPAF mutated at intermolecular dimerization, catalytic, or cleavage residues still produced active CPAF, although at a lower efficiency, indicating that CPAF can tolerate more mutations inside Chlamydia-infected cells than in cell-free systems. Only by simultaneously mutating both intermolecular dimerization and catalytic residues was CPAF activation completely blocked during infection, both indicating the importance of the critical residues identified in the cell-free systems and exploring the limit of CPAF's tolerance for mutations in the intracellular environment. We further found that active CPAF was always detected in the host cell cytoplasm while nonactive CPAF was restricted to within the chlamydial inclusions, regardless of how the infected cell samples were treated. Thus, CPAF translocation into the host cell cytoplasm correlates with CPAF enzymatic activity and is not altered by sample treatment conditions. These observations have provided new evidence for CPAF activation and translocation, which should encourage continued investigation of CPAF in chlamydial pathogenesis.

Two Phages, phiIPLA-RODI and phiIPLA-C1C, Lyse Mono- and Dual-Species Staphylococcal Biofilms

Phage therapy is a promising option for fighting against staphylococcal infections. Two lytic phages, vB_SauM_phiIPLA-RODI (phiIPLA-RODI) and vB_SepM_phiIPLA-C1C (phiIPLA-C1C), belonging to the Myoviridae family and exhibiting wide host ranges, were characterized in this study. The complete genome sequences comprised 142,348 bp and 140,961 bp and contained 213 and 203 open reading frames, respectively. The gene organization was typical of Spounavirinae members, with long direct terminal repeats (LTRs), genes grouped into modules not clearly separated from each other, and several group I introns. In addition, four genes encoding tRNAs were identified in phiIPLA-RODI. Comparative DNA sequence analysis showed high similarities with two phages, GH15 and 676Z, belonging to the Twort-like virus genus (nucleotide identities of >84%); for phiIPLA-C1C, a high similarity with phage phiIBB-SEP1 was observed (identity of 80%). Challenge assays of phages phiIPLA-RODI and phiIPLA-C1C against planktonic staphylococcal cells confirmed their lytic ability, as they were able to remove 5 log units in 8 h. Exposure of biofilms to phages phiIPLA-RODI and phiIPLA-C1C reduced the amount of adhered bacteria to about 2 log units in both monospecies and dual-species biofilms, but phiIPLA-RODI turned out to be as effective as the mixture of both phages. Moreover, the frequencies of bacteriophage-insensitive mutants (BIMs) of Staphylococcus aureus and S. epidermidis with resistance to phiIPLA-RODI and phiIPLA-C1C were low, at 4.05 × 10(-7) ± 2.34 × 10(-9) and 1.1 × 10(-7) ± 2.08 × 10(-9), respectively. Overall, a generally reduced fitness in the absence of phages was observed for BIMs, which showed a restored phage-sensitive phenotype in a few generations. These results confirm that lytic bacteriophages can be efficient biofilm-disrupting agents, supporting their potential as antimicrobials against staphylococcal infections.

Vibrio vulnificus phage PV94 is closely related to temperate phages of V. cholerae and other Vibrio species

Background

Vibrio vulnificus is an important pathogen which can cause serious infections in humans. Yet, there is limited knowledge on its virulence factors and the question whether temperate phages might be involved in pathogenicity, as is the case with V. cholerae. Thus far, only two phages (SSP002 and VvAW1) infecting V. vulnificus have been genetically characterized. These phages were isolated from the environment and are not related to Vibrio cholerae phages. The lack of information on temperate V. vulnificus phages prompted us to isolate those phages from lysogenic strains and to compare them with phages of other Vibrio species.

Results

In this study the temperate phage PV94 was isolated from a V. vulnificus biotype 1 strain by mitomycin C induction. PV94 is a myovirus whose genome is a linear double-stranded DNA of 33,828 bp with 5′-protruding ends. Sequence analysis of PV94 revealed a modular organization of the genome. The left half of the genome comprising the immunity region and genes for the integrase, terminase and replication proteins shows similarites to V. cholerae kappa phages whereas the right half containing genes for structural proteins is closely related to a prophage residing in V. furnissii NCTC 11218.

Conclusion

We present the first genomic sequence of a temperate phage isolated from a human V. vulnificus isolate. The sequence analysis of the PV94 genome demonstrates the wide distribution of closely related prophages in various Vibrio species. Moreover, the mosaicism of the PV94 genome indicates a high degree of horizontal genetic exchange within the genus Vibrio, by which V. vulnificus might acquire virulence-associated genes from other species.

Mycobacterial RNA polymerase requires a U-tract at intrinsic terminators and is aided by NusG at suboptimal terminators

Intrinsic terminators, which encode GC-rich RNA hairpins followed immediately by a 7-to-9-nucleotide (nt) U-rich “U-tract,” play principal roles of punctuating and regulating transcription in most bacteria. However, canonical intrinsic terminators with strong U-tracts are underrepresented in some bacterial lineages, notably mycobacteria, leading to proposals that their RNA polymerases stop at noncanonical intrinsic terminators encoding various RNA structures lacking U-tracts. We generated recombinant forms of mycobacterial RNA polymerase and its major elongation factors NusA and NusG to characterize mycobacterial intrinsic termination. Using in vitro transcription assays devoid of possible mycobacterial contaminants, we established that mycobacterial RNA polymerase terminates more efficiently than Escherichia coli RNA polymerase at canonical terminators with imperfect U-tracts but does not terminate at putative terminators lacking U-tracts even in the presence of mycobacterial NusA and NusG. However, mycobacterial NusG exhibits a novel termination-stimulating activity that may allow intrinsic terminators with suboptimal U-tracts to function efficiently.

Bacteria rely on transcription termination to define and regulate units of gene expression. In most bacteria, precise termination and much regulation by attenuation are accomplished by intrinsic terminators that encode GC-rich hairpins and U-tracts necessary to disrupt stable transcription elongation complexes. Thus, the apparent dearth of canonical intrinsic terminators with recognizable U-tracts in mycobacteria is of significant interest both because noncanonical intrinsic terminators could reveal novel routes to destabilize transcription complexes and because accurate understanding of termination is crucial for strategies to combat mycobacterial diseases and for computational bioinformatics generally. Our finding that mycobacterial RNA polymerase requires U-tracts for intrinsic termination, which can be aided by NusG, will guide future study of mycobacterial transcription and aid improvement of predictive algorithms to annotate bacterial genome sequences.

Identification of a conserved branched RNA structure that functions as a factor-independent terminator

Anti-Q is a small RNA encoded on pCF10, an antibiotic resistance plasmid of Enterococcus faecalis, which negatively regulates conjugation of the plasmid. In this study we sought to understand how Anti-Q is generated relative to larger transcripts of the same operon. We found that Anti-Q folds into a branched structure that functions as a factor-independent terminator. In vitro and in vivo, termination is dependent on the integrity of this structure as well as the presence of a 3' polyuridine tract, but is not dependent on other downstream sequences. In vitro, terminated transcripts are released from RNA polymerase after synthesis. In vivo, a mutant with reduced termination efficiency demonstrated loss of tight control of conjugation function. A search of bacterial genomes revealed the presence of sequences that encode Anti-Q-like RNA structures. In vitro and in vivo experiments demonstrated that one of these functions as a terminator. This work reveals a previously unappreciated flexibility in the structure of factor-independent terminators and identifies a mechanism for generation of functional small RNAs; it should also inform annotation of bacterial sequence features, such as terminators, functional sRNAs, and operons.

The Bacillus subtilis TRAP protein can induce transcription termination in the leader region of the tryptophan biosynthetic (trp) operon independent of the trp attenuator RNA

In Bacillus subtilis, transcription of the tryptophan biosynthetic operon is regulated by an attenuation mechanism. When intracellular tryptophan levels are high, the TRAP protein binds to the 5′ leader region of the nascent trp mRNA and induces transcription termination prior to the structural genes. In limiting tryptophan, TRAP does not bind and the operon is transcribed. Two competing RNA secondary structures termed the antiterminator and terminator (attenuator) can form in the leader region RNA. In prior attenuation models, the only role of TRAP binding was to alter the RNA secondary structure to allow formation of the attenuator, which has been thought function as an intrinsic transcription terminator. However, recent studies have shown that the attenuator is not an effective intrinsic terminator. From these studies it was not clear whether TRAP functions independently or requires the presence of the attenuator RNA structure. Hence we have further examined the role of the attenuator RNA in TRAP-mediated transcription termination. TRAP was found to cause efficient transcription termination in the trp leader region in vivo when the attenuator was mutated or deleted. However, TRAP failed to induce transcription termination at these mutant attenuators in a minimal in vitro transcription system with B. subtilis RNA polymerase. Further studies using this system showed that NusA as well as the timing of TRAP binding to RNA play a role in the observed differences in vivo and in vitro.

Tiling array study of MNNG treated Escherichia coli reveals a widespread transcriptional response

The alkylating agent N-methyl-N'-nitro-N-nitrosoguanidine (MNNG) is known to trigger the adaptive response by inducing the ada-regulon – consisting of three DNA repair enzymes Ada, AlkB, AlkA and the enigmatic AidB. We have applied custom designed tiling arrays to study transcriptional changes in Escherichia coli following a MNNG challenge. Along with the expected upregulation of the adaptive response genes (ada, alkA and alkB), we identified a number of differentially expressed transcripts, both novel and annotated. This indicates a wider regulatory response than previously documented. There were 250 differentially-expressed and 2275 similarly-expressed unannotated transcripts. We found novel upregulation of several stress-induced transcripts, including the SOS inducible genes recN and tisAB, indicating a novel role for these genes in alkylation repair. Furthermore, the ada-regulon A and B boxes were found to be insufficient to explain the regulation of the adaptive response genes after MNNG exposure, suggesting that additional regulatory elements must be involved.

Proteins of novel lactic acid bacteria from Apis mellifera mellifera: an insight into the production of known extra-cellular proteins during microbial stress

BACKGROUND

Lactic acid bacteria (LAB) has been considered a beneficial bacterial group, found as part of the microbiota of diverse hosts, including humans and various animals. However, the mechanisms of how hosts and LAB interact are still poorly understood. Previous work demonstrates that 13 species of Lactobacillus and Bifidobacterium from the honey crop in bees function symbiotically with the honeybee. They protect each other, their hosts, and the surrounding environment against severe bee pathogens, bacteria, and yeasts. Therefore, we hypothesized that these LAB under stress, i.e. in their natural niche in the honey crop, are likely to produce bioactive substances with antimicrobial activity.

RESULTS

The genomic analysis of the LAB demonstrated varying genome sizes ranging from 1.5 to 2.2 mega-base pairs (Mbps) which points out a clear difference within the protein gene content, as well as specialized functions in the honeybee microbiota and their adaptation to their host. We demonstrate a clear variation between the secreted proteins of the symbiotic LAB when subjected to microbial stressors. We have identified that 10 of the 13 LAB produced extra-cellular proteins of known or unknown function in which some are arranged in interesting putative operons that may be involved in antimicrobial action, host interaction, or biofilm formation. The most common known extra-cellular proteins secreted were enzymes, DNA chaperones, S-layer proteins, bacteriocins, and lysozymes. A new bacteriocin may have been identified in one of the LAB symbionts while many proteins with unknown functions were produced which must be investigated further.

CONCLUSIONS

The 13 LAB symbionts likely play different roles in their natural environment defending their niche and their host and participating in the honeybee's food production. These roles are partly played through producing extracellular proteins on exposure to microbial stressors widely found in natural occurring flowers. Many of these secreted proteins may have a putative antimicrobial function. In the future, understanding these processes in this complicated environment may lead to novel applications of honey crop LAB proteins.

Isolation and characterization of phages with lytic activity against methicillin-resistant Staphylococcus aureus strains belonging to clonal complex 398

Some years ago, MRSA clonal complex (CC) 398 emerged, which spread extensively in livestock animals. People in contact with food production animals are at high risk of colonization. A reduction of MRSA CC398 in livestock might be achieved by application of virulent phages. However, there have not yet been any reports published on phages lysing MRSA CC398 strains. In this study, three virulent phages (PSa1, PSa2 and PSa3) with lytic activity against MRSA CC398 strains were isolated from German pig farms. Morphologically, the phages are members of the family Podoviridae, and they exhibited an identical host range. They lysed 52 (60 %) out of 86 tested MRSA CC398 strains representing 18 different spa types. While the PSa1 and PSa3 genomes have a similar size of approximately 17.5 kb, the PSa2 genome is somewhat larger (ca. 18.5 kb). Southern hybridization revealed strong DNA homologies between the phages, which was confirmed by sequence analysis of cloned restriction fragments and PCR products. Moreover, the whole PSa3 genomic sequence (17,602 bp) showed a close relationship to 44AHJD-like phages, which are not known to contain virulence-associated genes. To assess whether these phages might be candidates for applications, in vitro experiments were carried out in which the number of MRSA CC398 cells could be reduced by up to four log10 units. The phages were stable at a wide range of temperatures and pH values.

Complete Genome Sequence of Pseudomonas denitrificans ATCC 13867

Pseudomonas denitrificans ATCC 13867, a Gram-negative facultative anaerobic bacterium, is known to produce vitamin B12 under aerobic conditions. This paper reports the annotated whole-genome sequence of the circular chromosome of this organism.

Isolation, characterization and complete genome sequence of PhaxI: a phage of Escherichia coli O157 : H7

Bacteriophages are considered as promising biological agents for the control of infectious diseases. Sequencing of their genomes can ascertain the absence of antibiotic resistance, toxin or virulence genes. The anti-O157 : H7 coliphage, PhaxI, was isolated from a sewage sample in Iran. Morphological studies by transmission electron microscopy showed that it has an icosahedral capsid of 85-86 nm and a contractile tail of 115×15 nm. PhaxI contains dsDNA composed of 156 628 nt with a G+C content of 44.5 mol% that encodes 209 putative proteins. In MS analysis of phage particles, 92 structural proteins were identified. PhaxI lyses Escherichia coli O157 : H7 in Luria-Bertani medium and milk, has an eclipse period of 20 min and a latent period of 40 min, and has a burst size of about 420 particles per cell. PhaxI is a member of the genus 'Viunalikevirus' of the family Myoviridae and is specific for E. coli O157 : H7.

Enterococcal Rgg-like regulator ElrR activates expression of the elrA operon

The Enterococcus faecalis leucine-rich protein ElrA promotes virulence by stimulating bacterial persistence in macrophages and production of the interleukin-6 (IL-6) cytokine. The ElrA protein is encoded within an operon that is poorly expressed under laboratory conditions but induced in vivo. In this study, we identify ef2687 (renamed elrR), which encodes a member of the Rgg (regulator gene for glucosyltransferase) family of putative regulatory proteins. Using quantitative reverse transcription-PCR, translational lacZ fusions, and electrophoretic mobility shift assays, we demonstrate that ElrR positively regulates expression of elrA. These results correlate with the attenuated virulence of the ΔelrR strain in a mouse peritonitis model. Virulence of simple and double elrR and elrA deletion mutants also suggests a remaining ElrR-independent expression of elrA in vivo and additional virulence-related genes controlled by ElrR.

Role of MerH in mercury resistance in the archaeon Sulfolobus solfataricus

Crenarchaeota include extremely thermoacidophilic organisms that thrive in geothermal environments dominated by sulfidic ores and heavy metals such as mercury. Mercuric ion, Hg(II), inactivates transcription in the crenarchaeote Sulfolobus solfataricus and simultaneously derepresses transcription of a resistance operon, merHAI, through interaction with the MerR transcription factor. While mercuric reductase (MerA) is required for metal resistance, the role of MerH, an adjacent small and predicted product of an ORF, has not been explored. Inactivation of MerH either by nonsense mutation or by in-frame deletion diminished Hg(II) resistance of mutant cells. Promoter mapping studies indicated that Hg(II) sensitivity of the merH nonsense mutant arose through transcriptional polarity, and its metal resistance was restored partially by single copy merH complementation. Since MerH was not required in vitro for MerA-catalysed Hg(II) reduction, MerH may play an alternative role in metal resistance. Inductively coupled plasma-mass spectrometry analysis of the MerH deletion strain following metal challenge indicated that there was prolonged retention of intracellular Hg(II). Finally, a reduced rate of mer operon induction in the merH deletion mutant suggested that the requirement for MerH could result from metal trafficking to the MerR transcription factor.

Investigation of Mannheimia haemolytica bacteriophages relative to host diversity

AIMS

This study aimed to characterize the impact of lytic and temperate bacteriophages on the genetic and phenotypic diversity of Mannheimia haemolytica from feedlot cattle.

METHODS AND RESULTS

Strictly lytic phages were not detected from bovine nasopharyngeal (n = 689) or water trough (n = 30) samples, but Myoviridae- or Siphoviridae-like phages were induced from 54 of 72 M. haemolytica strains by mitomycin C, occasionally from the same strain. Phages with similar restriction fragment length polymorphism profiles (RFLP ≥70% relatedness) shared common host serotypes 1 or 2 (P < 0·0001). Likewise, phages with similar RFLP tended to occur in genetically related host bacteria (70-79% similarity). Host range assays showed that seven phages from host serotypes 1, 2 and 6 lysed representative strains of serotypes 1, 2 or 8. The genome of vB_MhM_1152AP from serotype 6 was found to be collinear with P2-like phage φMhaA1-PHL101.

CONCLUSIONS

Prophages are a significant component of the genome of M. haemolytica and contribute significantly to host diversity. Further characterization of the role of prophage in virulence and persistence of M. haemolytica in cattle could provide insight into approaches to control this potential respiratory pathogen.

SIGNIFICANCE AND IMPACT OF THE STUDY

This study demonstrated that prophages are widespread within the genome of M. haemolytica isolates and emphasized the challenge of isolating lytic phage as a therapeutic against this pathogen.

The genome of VP3, a T7-like phage used for the typing of Vibrio cholerae

The bacteriophage VP3 is used in a phage-biotyping scheme as one of the typing phages of Vibrio cholerae O1 biotype El Tor strains. Here, we have sequenced and analyzed its genome. The genome consists of 39,481 bp with an overall G + C content of 42.6 %. Fifty-two open reading frames (ORFs) were predicted. Within the genome, 17 highly conserved phage promoters and 6 rho-independent terminators were predicted. When assessed with Rluc as a reporter gene, 12 of 16 cloned VP3 promoters showed activity in the host strain V. cholerae biotype El Tor. Based on the temporal expression pattern detected using reverse transcription PCR (RT-PCR), VP3 ORFs can be classed into four groups, arranged according to their order in the VP3 genome. Terminators T1 and T6 are presumed to work efficiently. Sequencing of the typing phage VP3 of V. cholerae reveals its evolutionary subdivisions from the members of T7-like phages of Escherichia coli. Knowledge of VP3 expands the known host range of T7-like phages and will promote understanding the different infection mechanisms used by members of this genus.

Measurement and modeling of intrinsic transcription terminators

The reliable forward engineering of genetic systems remains limited by the ad hoc reuse of many types of basic genetic elements. Although a few intrinsic prokaryotic transcription terminators are used routinely, termination efficiencies have not been studied systematically. Here, we developed and validated a genetic architecture that enables reliable measurement of termination efficiencies. We then assembled a collection of 61 natural and synthetic terminators that collectively encode termination efficiencies across an ∼800-fold dynamic range within Escherichia coli. We simulated co-transcriptional RNA folding dynamics to identify competing secondary structures that might interfere with terminator folding kinetics or impact termination activity. We found that structures extending beyond the core terminator stem are likely to increase terminator activity. By excluding terminators encoding such context-confounding elements, we were able to develop a linear sequence-function model that can be used to estimate termination efficiencies (r = 0.9, n = 31) better than models trained on all terminators (r = 0.67, n = 54). The resulting systematically measured collection of terminators should improve the engineering of synthetic genetic systems and also advance quantitative modeling of transcription termination.