BAGEL2: mining for bacteriocins in genomic data

Mining bacterial genomes for bacteriocins is a challenging task due to the substantial structure and sequence diversity, and generally small sizes, of these antimicrobial peptides. Major progress in the research of antimicrobial peptides and the ever-increasing quantities of genomic data, varying from (un)finished genomes to meta-genomic data, led us to develop the significantly improved genome mining software BAGEL2, as a follow-up of our previous BAGEL software. BAGEL2 identifies putative bacteriocins on the basis of conserved domains, physical properties and the presence of biosynthesis, transport and immunity genes in their genomic context. The software supports parameter-free, class-specific mining and has high-throughput capabilities. Besides building an expert validated bacteriocin database, we describe the development of novel Hidden Markov Models (HMMs) and the interpretation of combinations of HMMs via simple decision rules for prediction of bacteriocin (sub-)classes. Furthermore, the genetic context is automatically annotated based on (combinations of) PFAM domains and databases of known context genes. The scoring system was fine-tuned using expert knowledge on data derived from screening all bacterial genomes currently available at the NCBI. BAGEL2 is freely accessible at http://bagel2.molgenrug.nl.

Penicillin-binding protein folding is dependent on the PrsA peptidyl-prolyl cis-trans isomerase in Bacillus subtilis

Summary The PrsA protein is a membrane-anchored peptidyl-prolyl cis-trans isomerase in Bacillus subtilis and most other Gram-positive bacteria. It catalyses the post-translocational folding of exported proteins and is essential for normal growth of B. subtilis. We studied the mechanism behind this indispensability. We could construct a viable prsA null mutant in the presence of a high concentration of magnesium. Various changes in cell morphology in the absence of PrsA suggested that PrsA is involved in the biosynthesis of the cylindrical lateral wall. Consistently, four penicillin-binding proteins (PBP2a, PBP2b, PBP3 and PBP4) were unstable in the absence of PrsA, while muropeptide analysis revealed a 2% decrease in the peptidoglycan cross-linkage index. Misfolded PBP2a was detected in PrsA-depleted cells, indicating that PrsA is required for the folding of this PBP either directly or indirectly. Furthermore, strongly increased uniform staining of cell wall with a fluorescent vancomycin was observed in the absence of PrsA. We also demonstrated that PrsA is a dimeric or oligomeric protein which is localized at distinct spots organized in a helical pattern along the cell membrane. These results suggest that PrsA is essential for normal growth most probably as PBP folding is dependent on this PPIase.

Gene position within a long transcript as a determinant for stochastic switching in bacteria

How cultures of genetically identical cells bifurcate into distinct phenotypic subpopulations under uniform growth conditions is an important question in developmental biology of relevance even to relatively simple developmental systems, such as spore formation in bacteria. A growing Bacillus subtilis culture consists of either cells that are motile and can swim or cells that are non-motile and are chained together. In this issue of Molecular Microbiology, Cozy and Kearns show that the probability of a cell to become motile depends on the position of the sigD gene within the long (27 kb) motility operon. sigD encodes the alternative sigma factor sigma(D) that, together with RNA polymerase, drives expression of genes required for cell separation and the assembly of flagella. sigD is the penultimate gene of the B. subtilis motility operon and, in the control strain approximately, 70% of the cells are motile. When sigD was moved upstream within the operon, a larger fraction of cells became motile (up to 100%). This study highlights that the position of a gene within an operon can have a large impact on the control of gene expression. Furthermore, it suggests that RNA polymerase processivity or mRNA turnover can play important roles as sources of noise in bacterial development, and that gene position might be an unrecognized and possibly widespread mechanism to regulate phenotypic variation.

Adaptation of Hansenula polymorpha to methanol: a transcriptome analysis

BACKGROUND

Methylotrophic yeast species (e.g. Hansenula polymorpha, Pichia pastoris) can grow on methanol as sole source of carbon and energy. These organisms are important cell factories for the production of recombinant proteins, but are also used in fundamental research as model organisms to study peroxisome biology. During exponential growth on glucose, cells of H. polymorpha typically contain a single, small peroxisome that is redundant for growth while on methanol multiple, enlarged peroxisomes are present. These organelles are crucial to support growth on methanol, as they contain key enzymes of methanol metabolism.In this study, changes in the transcriptional profiles during adaptation of H. polymorpha cells from glucose- to methanol-containing media were investigated using DNA-microarray analyses.

RESULTS

Two hours after the shift of cells from glucose to methanol nearly 20% (1184 genes) of the approximately 6000 annotated H. polymorpha genes were significantly upregulated with at least a two-fold differential expression. Highest upregulation (> 300-fold) was observed for the genes encoding the transcription factor Mpp1 and formate dehydrogenase, an enzyme of the methanol dissimilation pathway. Upregulated genes also included genes encoding other enzymes of methanol metabolism as well as of peroxisomal beta-oxidation.A moderate increase in transcriptional levels (up to 4-fold) was observed for several PEX genes, which are involved in peroxisome biogenesis. Only PEX11 and PEX32 were higher upregulated. In addition, an increase was observed in expression of the several ATG genes, which encode proteins involved in autophagy and autophagy processes. The strongest upregulation was observed for ATG8 and ATG11.Approximately 20% (1246 genes) of the genes were downregulated. These included glycolytic genes as well as genes involved in transcription and translation.

CONCLUSION

Transcriptional profiling of H. polymorpha cells shifted from glucose to methanol showed the expected downregulation of glycolytic genes together with upregulation of the methanol utilisation pathway. This serves as a confirmation and validation of the array data obtained. Consistent with this, also various PEX genes were upregulated. The strong upregulation of ATG genes is possibly due to induction of autophagy processes related to remodeling of the cell architecture required to support growth on methanol. These processes may also be responsible for the enhanced peroxisomal beta-oxidation, as autophagy leads to recycling of membrane lipids. The prominent downregulation of transcription and translation may be explained by the reduced growth rate on methanol (td glucose 1 h vs td methanol 4.5 h).

Response of Bacillus cereus ATCC 14579 to challenges with sublethal concentrations of enterocin AS-48

Background

Enterocin AS-48 is produced by Enterococcus faecalis S48 to compete with other bacteria in their environment. Due to its activity against various Gram positive and some Gram negative bacteria it has clear potential for use as a food preservative. Here, we studied the effect of enterocin AS-48 challenges on vegetative cells of Bacillus cereus ATCC 14579 by use of transcriptome analysis.

Results

Of the 5200 genes analysed, expression of 24 genes was found to change significantly after a 30 min treatment with a subinhibitory bacteriocin concentration of 0.5 μg/ml. Most of up-regulated genes encode membrane-associated or secreted proteins with putative transmembrane segments or signal sequences, respectively. One operon involved in arginine metabolism was significantly downregulated. The BC4206-BC4207 operon was found to be the most upregulated target in our experiments. BC4206 codes for a PadR type transcriptional regulator, while BC4207 codes for a hypothetical membrane protein. The operon structure and genes are conserved in B. cereus and B. thuringiensis species, but are not present in B. anthracis and B. subtilis. Using real-time qPCR, we show that these genes are upregulated when we treated the cells with AS-48, but not upon nisin treatment. Upon overexpression of BC4207 in B. cereus, we observed an increased resistance against AS-48. Expression of BC4207 in B. subtilis 168, which lacks this operon also showed increased resistance against AS-48.

Conclusion

BC4207 membrane protein is involved in the resistance mechanism of B. cereus cells against AS-48.

Mechanisms and evolution of control logic in prokaryotic transcriptional regulation

A major part of organismal complexity and versatility of prokaryotes resides in their ability to fine-tune gene expression to adequately respond to internal and external stimuli. Evolution has been very innovative in creating intricate mechanisms by which different regulatory signals operate and interact at promoters to drive gene expression. The regulation of target gene expression by transcription factors (TFs) is governed by control logic brought about by the interaction of regulators with TF binding sites (TFBSs) in cis-regulatory regions. A factor that in large part determines the strength of the response of a target to a given TF is motif stringency, the extent to which the TFBS fits the optimal TFBS sequence for a given TF. Advances in high-throughput technologies and computational genomics allow reconstruction of transcriptional regulatory networks in silico. To optimize the prediction of transcriptional regulatory networks, i.e., to separate direct regulation from indirect regulation, a thorough understanding of the control logic underlying the regulation of gene expression is required. This review summarizes the state of the art of the elements that determine the functionality of TFBSs by focusing on the molecular biological mechanisms and evolutionary origins of cis-regulatory regions.

Directionality and coordination of dehydration and ring formation during biosynthesis of the lantibiotic nisin

The lantibiotic nisin is a potent antimicrobial substance, which contains unusual lanthionine rings and dehydrated amino acid residues and is produced by Lactococcus lactis. Recently, the nisin biosynthetic machinery has been applied to introduce lanthionine rings in peptides other than nisin with potential therapeutic use. Due to difficulties in the isolation of the proposed synthetase complex (NisBTC), mechanistic information concerning the enzymatic biosynthesis of nisin is scarce. Here, we present the molecular characterization of a number of nisin mutants that affect ring formation. We have investigated in a systematic manner how these mutations influence dehydration events, which are performed enzymatically by the dehydratase NisB. Specific mutations that hampered ring formation allowed for the dehydration of serine residues that directly follow the rings and are normally unmodified. The combined information leads to the conclusion that 1) nisin biosynthesis is an organized directional process that starts at the N terminus of the molecule and continues toward the C terminus, and 2) NisB and NisC are alternating enzymes, whose activities follow one after another in a repetitive way. Thus, the dehydration and cyclization processes are not separated in time and space. On the basis of these results and previous knowledge, a working model for the sequence of events in the maturation of nisin is proposed.

Ubiquitous late competence genes in Bacillus species indicate the presence of functional DNA uptake machineries

Natural competence for genetic transformation, i.e. the ability to take up DNA and stably integrate it in the genome, has so far only been observed in the bacterial kingdom (both in gram-negative and gram-positive species) and may contribute to survival under adverse growth conditions. Bacillus subtilis, the model organism for the Bacillus genus, possesses a well-characterized competence machinery. Phylogenetic analysis of several genome sequences of different Bacillus species reveals the presence of many, but not all genes potentially involved in competence and its regulation. The recent demonstration of functional DNA uptake by B. cereus supports the significance of our genome analyses and shows that the ability for functional DNA uptake might be widespread among Bacilli.

Strain-specific impact of PsaR of Streptococcus pneumoniae on global gene expression and virulence

Previous studies have indicated that PsaR of Streptococcus pneumoniae is a manganese-dependent regulator, negatively affecting the expression of at least seven genes. Here, we extended these observations by transcriptome and proteome analysis of psaR mutants in strains D39 and TIGR4. The microarray analysis identified three shared PsaR targets: the psa operon, pcpA and prtA. In addition, we found 31 genes to be regulated by PsaR in D39 only, most strikingly a cellobiose-specific phosphotransferase system (PTS) and a putative bacteriocin operon (sp0142-sp0146). In TIGR4, 14 PsaR gene targets were detected, with the rlrA pathogenicity islet being the most pronounced. Proteomics confirmed most of the shared gene targets. To examine the contribution of PsaR to pneumococcal virulence, we compared D39 and TIGR4 wild-type (wt) and psaR mutants in three murine infection models. During colonization, no clear effect was observed of the psaR mutation in either D39 or TIGR4. In the pneumonia model, small but significant differences were observed in the lungs of mice infected with either D39wt or DeltapsaR: D39DeltapsaR had an initial advantage in survival in the lungs. Conversely, TIGR4DeltapsaR-infected mice had significantly lower bacterial loads at 24 h only. Finally, during experimental bacteraemia, D39DeltapsaR-infected mice had significantly lower bacterial loads in the bloodstream than wt-infected mice for the first 24 h of infection. TIGR4DeltapsaR showed attenuation at 36 h only. In conclusion, our results show that PsaR of D39 and TIGR4 has a strain-specific role in global gene expression and in the development of bacteraemia in mice.

Optimization of protein secretion by Bacillus subtilis

The gram-positive bacterium Bacillus subtilis is widely known for its capacity to produce and secrete large amounts of industrially relevant proteins, mostly endogenous enzymes like proteases and lipases. The use of B. subtilis has many advantages, such as its GRAS status and the easy and inexpensive culturing methods that can result in very high cell densities. Over the years many patents have been filed regarding the optimization of protein secretion by B. subtilis. For almost every step in the production and secretion process, from promoter optimization to deletion of extracellular proteases, patents have been claimed. An overview of the current literature and patents on these subjects is given. We will discuss recent patents regarding the optimization of protein overexpression and secretion in B. subtilis. A patent claiming modification of B. subtilis SecA will be discussed in more detail. Another recent patent claims a positive effect of heterologous protein secretion upon reduced expression of the yusZ and/or yusX genes, encoding putative oligopeptidases. Improvements are being made continuously, although many depend on the character of the specific protein under study.

MOTIFATOR: detection and characterization of regulatory motifs using prokaryote transcriptome data

SUMMARY

Unraveling regulatory mechanisms (e.g. identification of motifs in cis-regulatory regions) remains a major challenge in the analysis of transcriptome experiments. Existing applications identify putative motifs from gene lists obtained at rather arbitrary cutoff and require additional manual processing steps. Our standalone application MOTIFATOR identifies the most optimal parameters for motif discovery and creates an interactive visualization of the results. Discovered putative motifs are functionally characterized, thereby providing valuable insight in the biological processes that could be controlled by the motif.

AVAILABILITY

MOTIFATOR is freely available at http://www.motifator.nl.

DISCLOSE : DISsection of CLusters Obtained by SEries of transcriptome data using functional annotations and putative transcription factor binding sites

BACKGROUND

A typical step in the analysis of gene expression data is the determination of clusters of genes that exhibit similar expression patterns. Researchers are confronted with the seemingly arbitrary choice between numerous algorithms to perform cluster analysis.

RESULTS

We developed an exploratory application that benchmarks the results of clustering methods using functional annotations. In addition, a de novo DNA motif discovery algorithm is integrated in our program which identifies overrepresented DNA binding sites in the upstream DNA sequences of genes from the clusters that are indicative of sites of transcriptional control. The performance of our program was evaluated by comparing the original results of a time course experiment with the findings of our application.

CONCLUSION

DISCLOSE assists researchers in the prokaryotic research community in systematically evaluating results of the application of a range of clustering algorithms to transcriptome data. Different performance measures allow to quickly and comprehensively determine the best suited clustering approach for a given dataset.

Bistability, epigenetics, and bet-hedging in bacteria

Clonal populations of microbial cells often show a high degree of phenotypic variability under homogeneous conditions. Stochastic fluctuations in the cellular components that determine cellular states can cause two distinct subpopulations, a property called bistability. Phenotypic heterogeneity can be readily obtained by interlinking multiple gene regulatory pathways, effectively resulting in a genetic logic-AND gate. Although switching between states can occur within the cells' lifetime, cells can also pass their cellular state over to the next generation by a mechanism known as epigenetic inheritance and thus perpetuate the phenotypic state. Importantly, heterogeneous populations can demonstrate increased fitness compared with homogeneous populations. This suggests that microbial cells employ bet-hedging strategies to maximize survival. Here, we discuss the possible roles of interlinked bistable networks, epigenetic inheritance, and bet-hedging in bacteria.

Characterization of the individual glucose uptake systems of Lactococcus lactis: mannose-PTS, cellobiose-PTS and the novel GlcU permease

According to previous reports, Lactococcus lactis imports glucose via two distinct phosphoenolpyruvate:phosphotransferase systems (mannose-PTS and cellobiose-PTS) and one or more unknown non-PTS permease(s). GlcU was identified as the sole non-PTS permease involved in the transport of glucose. Additionally, the biochemical properties of PTS(Man), PTS(Cel) and GlcU were characterized in double knockout mutants with glucose uptake restricted to a single system. Transport susceptibility to protonophores indicated that glucose uptake via GlcU is proton-motive force dependent. Competition assays revealed a high specificity of GlcU for glucose. Furthermore, the permease has low affinity for glucose and displays strong preference for the beta-anomer as shown by the profiles of consumption of the two glucose anomers studied by (13)C-NMR. Similar kinetic properties were found for PTS(Cel), while PTS(Man) is a high-affinity system recognizing equally well the two anomeric forms of glucose. Transcripts of the genes encoding the three transporters are present simultaneously in the parent strain NZ9000 as shown by reverse transcription-PCR. Investigation of the distribution of GlcU homologues among bacteria showed that these proteins are restricted to the low-GC Gram-positive Firmicutes. This work completes the identification of the glucose transport systems in L. lactis MG1363.

The twin-arginine translocation (Tat) systems from Bacillus subtilis display a conserved mode of complex organization and similar substrate recognition requirements

The twin arginine translocation (Tat) system transports folded proteins across the bacterial plasma membrane. In Gram-negative bacteria, membrane-bound TatABC subunits are all essential for activity, whereas Gram-positive bacteria usually contain only TatAC subunits. In Bacillus subtilis, two TatAC-type systems, TatAdCd and TatAyCy, operate in parallel with different substrate specificities. Here, we show that they recognize similar signal peptide determinants. Both systems translocate green fluorescent protein fused to three distinct Escherichia coli Tat signal peptides, namely DmsA, AmiA and MdoD, and mutagenesis of the DmsA signal peptide confirmed that both Tat pathways recognize similar targeting determinants within Tat signals. Although another E. coli Tat substrate, trimethylamine N-oxide reductase, was translocated by TatAdCd but not by TatAyCy, we conclude that these systems are not predisposed to recognize only specific Tat signal peptides, as suggested by their narrow substrate specificities in B. subtilis. We also analysed complexes involved in the second Tat pathway in B. subtilis, TatAyCy. This revealed a discrete TatAyCy complex together with a separate, homogeneous, approximately 200 kDa TatAy complex. The latter complex differs significantly from the corresponding E. coli TatA complexes, pointing to major structural differences between Tat complexes from Gram-negative and Gram-positive organisms. Like TatAd, TatAy is also detectable in the form of massive cytosolic complexes.

MINOMICS: visualizing prokaryote transcriptomics and proteomics data in a genomic context

We have developed MINOMICS, a tool that allows facile and in-depth visualization of prokaryotic transcriptomic and proteomic data in conjunction with genomics data. MINOMICS generates interactive linear genome maps in which multiple experimental datasets are displayed together with operon, regulatory motif, transcriptional promoter and transcriptional terminator information.

AVAILABILITY

MINOMICS is freely accessible at http://www.minomics.nl

Prosecutor: parameter-free inference of gene function for prokaryotes using DNA microarray data, genomic context and multiple gene annotation sources

Background

Despite a plethora of functional genomic efforts, the function of many genes in sequenced genomes remains unknown. The increasing amount of microarray data for many species allows employing the guilt-by-association principle to predict function on a large scale: genes exhibiting similar expression patterns are more likely to participate in shared biological processes.

Results

We developed Prosecutor, an application that enables researchers to rapidly infer gene function based on available gene expression data and functional annotations. Our parameter-free functional prediction method uses a sensitive algorithm to achieve a high association rate of linking genes with unknown function to annotated genes. Furthermore, Prosecutor utilizes additional biological information such as genomic context and known regulatory mechanisms that are specific for prokaryotes. We analyzed publicly available transcriptome data sets and used literature sources to validate putative functions suggested by Prosecutor. We supply the complete results of our analysis for 11 prokaryotic organisms on a dedicated website.

Conclusion

The Prosecutor software and supplementary datasets available at allow researchers working on any of the analyzed organisms to quickly identify the putative functions of their genes of interest. A de novo analysis allows new organisms to be studied.

Increased D-alanylation of lipoteichoic acid and a thickened septum are main determinants in the nisin resistance mechanism of Lactococcus lactis

Nisin is a post-translationally modified antimicrobial peptide produced by Lactococcus lactis which binds to lipid II in the membrane to form pores and inhibit cell-wall synthesis. A nisin-resistant (Nis(R)) strain of L. lactis, which is able to grow at a 75-fold higher nisin concentration than its parent strain, was investigated with respect to changes in the cell wall. Direct binding studies demonstrated that less nisin was able to bind to lipid II in the membranes of L. lactis Nis(R) than in the parent strain. In contrast to vancomycin binding, which showed ring-like binding, nisin was observed to bind in patches close to cell-division sites in both the wild-type and the Nis(R) strains. Comparison of modifications in lipoteichoic acid of the L. lactis strains revealed an increase in d-alanyl esters and galactose as substituents in L. lactis Nis(R), resulting in a less negatively charged cell wall. Moreover, the cell wall displays significantly increased thickness at the septum. These results indicate that shielding the membrane and thus the lipid II molecule, thereby decreasing abduction of lipid II and subsequent pore-formation, is a major defence mechanism of L. lactis against nisin.

LysM, a widely distributed protein motif for binding to (peptido)glycans

Bacteria retain certain proteins at their cell envelopes by attaching them in a non-covalent manner to peptidoglycan, using specific protein domains, such as the prominent LysM (Lysin Motif) domain. More than 4000 (Pfam PF01476) proteins of both prokaryotes and eukaryotes have been found to contain one or more Lysin Motifs. Notably, this collection contains not only truly secreted proteins, but also (outer-)membrane proteins, lipoproteins or proteins bound to the cell wall in a (non-)covalent manner. The motif typically ranges in length from 44 to 65 amino acid residues and binds to various types of peptidoglycan and chitin, most likely recognizing the N-acetylglucosamine moiety. Most bacterial LysM-containing proteins are peptidoglycan hydrolases with various cleavage specificities. Binding of certain LysM proteins to cells of Gram-positive bacteria has been shown to occur at specific sites, as binding elsewhere is hindered by the presence of other cell wall components such as lipoteichoic acids. Interestingly, LysM domains of certain plant kinases enable the plant to recognize its symbiotic bacteria or sense and induce resistance against fungi. This interaction is triggered by chitin-like compounds that are secreted by the symbiotic bacteria or released from fungi, demonstrating an important sensing function of LysMs.

The relative value of operon predictions

For most organisms, computational operon predictions are the only source of genome-wide operon information. Operon prediction methods described in literature are based on (a combination of) the following five criteria: (i) intergenic distance, (ii) conserved gene clusters, (iii) functional relation, (iv) sequence elements and (v) experimental evidence. The performance estimates of operon predictions reported in literature cannot directly be compared due to differences in methods and data used in these studies. Here, we survey the current status of operon prediction methods. Based on a comparison of the performance of operon predictions on Escherichia coli and Bacillus subtilis we conclude that there is still room for improvement. We expect that existing and newly generated genomics and transcriptomics data will further improve accuracy of operon prediction methods.

Assessment of CcpA-mediated catabolite control of gene expression in Bacillus cereus ATCC 14579

Background

The catabolite control protein CcpA is a transcriptional regulator conserved in many Gram-positives, controlling the efficiency of glucose metabolism. Here we studied the role of Bacillus cereus ATCC 14579 CcpA in regulation of metabolic pathways and expression of enterotoxin genes by comparative transcriptome analysis of the wild-type and a ccpA-deletion strain.

Results

Comparative analysis revealed the growth performance and glucose consumption rates to be lower in the B. cereus ATCC 14579 ccpA deletion strain than in the wild-type. In exponentially grown cells, the expression of glycolytic genes, including a non-phosphorylating glyceraldehyde-3-phosphate dehydrogenase that mediates conversion of D-glyceraldehyde 3-phosphate to 3-phospho-D-glycerate in one single step, was down-regulated and expression of gluconeogenic genes and genes encoding the citric acid cycle was up-regulated in the B. cereus ccpA deletion strain. Furthermore, putative CRE-sites, that act as binding sites for CcpA, were identified to be present for these genes. These results indicate CcpA to be involved in the regulation of glucose metabolism, thereby optimizing the efficiency of glucose catabolism. Other genes of which the expression was affected by ccpA deletion and for which putative CRE-sites could be identified, included genes with an annotated function in the catabolism of ribose, histidine and possibly fucose/arabinose and aspartate. Notably, expression of the operons encoding non-hemolytic enterotoxin (Nhe) and hemolytic enterotoxin (Hbl) was affected by ccpA deletion, and putative CRE-sites were identified, which suggests catabolite repression of the enterotoxin operons to be CcpA-dependent.

Conclusion

The catabolite control protein CcpA in B. cereus ATCC 14579 is involved in optimizing the catabolism of glucose with concomitant repression of gluconeogenesis and alternative metabolic pathways. Furthermore, the results point to metabolic control of enterotoxin gene expression and suggest that CcpA-mediated glucose sensing provides an additional mode of control in moderating the expression of the nhe and hbl operons in B. cereus ATCC 14579.

Supervised Lowess normalization of comparative genome hybridization data--application to lactococcal strain comparisons

Background

Array-based comparative genome hybridization (aCGH) is commonly used to determine the genomic content of bacterial strains. Since prokaryotes in general have less conserved genome sequences than eukaryotes, sequence divergences between the genes in the genomes used for an aCGH experiment obstruct determination of genome variations (e.g. deletions). Current normalization methods do not take into consideration sequence divergence between target and microarray features and therefore cannot distinguish a difference in signal due to systematic errors in the data or due to sequence divergence.

Results

We present supervised Lowess, or S-Lowess, an application of the subset Lowess normalization method. By using a predicted subset of array features with minimal sequence divergence between the analyzed strains for the normalization procedure we remove systematic errors from dual-dye aCGH data in two steps: (1) determination of a subset of conserved genes (i.e. likely conserved genes, LCG); and (2) using the LCG for subset Lowess normalization. Subset Lowess determines the correction factors for systematic errors in the subset of array features and normalizes all array features using these correction factors. The performance of S-Lowess was assessed on aCGH experiments in which differentially labeled genomic DNA fragments of Lactococcus lactis IL1403 and L. lactis MG1363 strains were hybridized to IL1403 DNA microarrays. Since both genomes are sequenced and gene deletions identified, the success rate of different aCGH normalization methods in detecting these deletions in the MG1363 genome were determined. S-Lowess detects 97% of the deletions, whereas other aCGH normalization methods detect up to only 60% of the deletions.

Conclusion

S-Lowess is implemented in a user-friendly web-tool accessible from . We demonstrate that it outperforms existing normalization methods and maximizes detection of genomic variation (e.g. deletions) from microbial aCGH data.

A minimal Tat system from a gram-positive organism: a bifunctional TatA subunit participates in discrete TatAC and TatA complexes

The Tat system transports folded proteins across bacterial and thylakoid membranes. In Gram-negative organisms, a TatABC substrate-binding complex and separate TatA complex are believed to coalesce to form an active translocon, with all three subunits essential for translocation. Most Gram-positive organisms lack a tatB gene, indicating major differences in organization and possible differences in mode of action. Here, we have studied Tat complexes encoded by the tatAdCd genes of Bacillus subtilis. Expression of tatAdCd in an Escherichia coli tat null mutant results in efficient export of a large, cofactor-containing E. coli Tat substrate, TorA. We show that the tatAd gene complements E. coli mutants lacking either tatAE or tatB, indicating a bifunctional role for this subunit in B. subtilis. Second, we have identified and characterized two distinct Tat complexes that are novel in key respects: a TatAdCd complex of approximately 230 kDa that is significantly smaller than the analogous E. coli TatABC complex (approximately 370 kDa on BN gels) and a separate TatAd complex. The latter is a discrete entity of approximately 270 kDa as judged by gel filtration chromatography, very different from the highly heterogeneous E. coli TatA complex that ranges in size from approximately 50 kDa to over 600 kDa. TatA heterogeneity has been linked to the varying size of Tat substrates being translocated, but the singular nature of the B. subtilis TatAd complex suggests that discrete TatAC and TatA complexes may form a single form of translocon.