Characterization and analysis of the Burkholderia pseudomallei BsaN virulence regulon
© Chen et al.; licensee BioMed Central Ltd 2014
Received: 16 April 2014
Accepted: 18 July 2014
Published: 1 August 2014
Burkholderia pseudomallei is a facultative intracellular pathogen and the causative agent of melioidosis. A conserved type III secretion system (T3SS3) and type VI secretion system (T6SS1) are critical for intracellular survival and growth. The T3SS3 and T6SS1 genes are coordinately and hierarchically regulated by a TetR-type regulator, BspR. A central transcriptional regulator of the BspR regulatory cascade, BsaN, activates a subset of T3SS3 and T6SS1 loci.
To elucidate the scope of the BsaN regulon, we used RNAseq analysis to compare the transcriptomes of wild-type B. pseudomallei KHW and a bsaN deletion mutant. The 60 genes positively-regulated by BsaN include those that we had previously identified in addition to a polyketide biosynthesis locus and genes involved in amino acid biosynthesis. BsaN was also found to repress the transcription of 51 genes including flagellar motility loci and those encoding components of the T3SS3 apparatus. Using a promoter-lacZ fusion assay in E. coli, we show that BsaN together with the chaperone BicA directly control the expression of the T3SS3 translocon, effector and associated regulatory genes that are organized into at least five operons (BPSS1516-BPSS1552). Using a mutagenesis approach, a consensus regulatory motif in the promoter regions of BsaN-regulated genes was shown to be essential for transcriptional activation.
BsaN/BicA functions as a central regulator of key virulence clusters in B. pseudomallei within a more extensive network of genetic regulation. We propose that BsaN/BicA controls a gene expression program that facilitates the adaption and intracellular survival of the pathogen within eukaryotic hosts.
Melioidosis is a serious and often fatal infectious disease common to Southeast Asia and Northern Australia caused by the Gram-negative soil bacterium Burkholderia pseudomallei. B. pseudomallei is a highly versatile pathogen capable of surviving inside mammalian cells and in many environmental niches. The bacterium can infect numerous animal species, amoebae, nematodes, and tomato plants –, and has been previously found within the tissues of exotic grasses in Australia . The environmental origin of B. pseudomallei and its promiscuous host range have shaped the hypothesis that some of its genetic loci evolved in the rhizosphere as anti-predation determinants that subsequently promote “accidental” virulence in humans and animals. In recent years, important advances have been made in understanding the pathogenic mechanisms of B. pseudomallei including the roles of the Type III and Type VI Secretion Systems (T3SS, T6SS) –. B. pseudomallei contains three T3SSs and six T6SSs, but only T3SS3 (also referred to as the Burkholderia secretion apparatus, or T3SSBsa) and T6SS1 are critical for pathogenesis in mice and hamsters ,,.
Expression of the T3SS3 and T6SS1 gene clusters is tightly controlled, both temporally and spatially, during the B. pseudomallei intracellular lifecycle. We have identified a regulatory cascade that coordinately activates T3SS3 and T6SS1 gene expression in growth medium and in infected mammalian cells ,. At the top of the cascade is the TetR-type regulator BspR that stimulates the expression of bprP. The bspR gene is located on chromosome 1 of the B. pseudomallei genome whereas bprP is part of the T3SS3 gene cluster on chromosome 2 . The ToxR-like BprP in turn activates genes encoding the structural components of T3SS3, including the araC-type regulatory gene bsaN. BsaN is important for the activation of T3SS3 effector and translocon gene expression, and several regulatory genes including bprC and virAG, whose gene products control T6SS1 expression . The mechanisms through which these transcriptional regulators control the expression of their target genes are not understood. It is also unclear whether these regulators are acting directly on the identified target genes or through as yet undiscovered intermediary regulators, and whether additional host cell cofactors are involved that may serve as intracellular signals.
Compared to T3SSs in other pathogens such as Pseudomonas, Salmonella and Shigella, only a limited number of effectors have been identified for B. pseudomallei T3SS3. One of the effector proteins secreted by T3SS3 is BopE, which is annotated to exhibit guanine nucleotide exchange factor activity and has been reported to facilitate invasion of epithelial cells . bopA is generally assumed to encode a T3SS3 effector since it is located adjacent to bopE, although T3SS3-dependent secretion of BopA has never been verified. Functionally, BopA has been described to promote resistance to LC3-associated autophagy and a bopA mutation results in an intracellular replication defect ,. A third effector protein, BopC (BPSS1516), was recently shown to be secreted via T3SS3, and bopC mutants were reported to be less invasive in epithelial cells  and to exhibit delayed endosome escape and reduced intracellular growth in J774 murine macrophages .
To determine the full extent of the BsaN regulon and examine whether BsaN activates the expression of additional effector genes, we performed global transcriptome analysis of B. pseudomallei KHW wildtype (WT) and a ΔbsaN mutant strain using RNAseq. Our analysis shows that 111 genes are under the direct or indirect transcriptional control of BsaN. In addition to activating loci associated with T3SS3, we demonstrate that BsaN functions to repress transcription of other loci. Thus, BsaN functions as a central regulatory factor within a more extensive network to facilitate the intracellular lifecycle of B. pseudomallei.
Identification of the BsaN regulon through RNAseq analysis
BsaN (BPSS1546 in the reference B. pseudomallei K96243 genome) was previously shown to function as a central regulator of a hierarchical cascade that activates effector and translocon genes of T3SS3 as well as several associated regulatory genes ,. Furthermore, BsaN was shown to activate the expression of certain T6SS1-associated genes including the two-component regulatory system locus virAG (BPSS1494, 1495), and the bim actin motility genes (BPSS1490-1493). To gain further insight into the BsaN regulon, massive parallel sequencing was performed on cDNA prepared from RNA isolated from wild-type and ΔbsaN mutant strains. We had previously shown that complementation of our ΔbsaN mutant with a bsaN plasmid could restore the secretion of the BopE effector , showing that our complementation restored protein expression of the effectors and that the mutation was specific to bsaN and not due to off target effects.
Between 16 and 56 million reads (n = 2 from 3 combined cultures) were obtained that aligned to non-ribosomal genes in the KHW  genome (Additional file 1: Table S1). Reads of the technical replicates displayed high reproducibility (R-value) (Additional file 1: Table S1) demonstrating that variability was not introduced through sample preparation or sequencing errors. The K96243 reference genome was co-aligned for ease of gene annotation. The nucleotide sequences of chromosomes I and II are 99.3 and 99.1% identical, respectively.
List of 60 genes that are expressed 3-fold and higher in the wild-type versus Δ bsaN mutant strains (p < 0.01)
Gene locus ID
Actin based motility
Bim actin polymerization protein
NRPKS/PKS biosynthesis locus
Glyoxalase/bleomycin resistance protein/dioxygenase
Amino acid biosynthesis and sugar uptake
methionine biosynthesis protein MetW
2-amino-3-ketobutyrate CoA ligase
Periplasmic binding protein (ribose binding)
T6SS-1 response regulator
T6SS-1 His kinase
T3SS-3 AraC-type regulator
T3SS-3 response regulator
T3SS-3 HNS-type regulator
List of 51 genes that are expressed 3-fold and lower in the wild-type versus Δ bsaN mutant strains (p < 0.01)
Gene locus ID
Flagellar hook-associated protein
Non-ribosomal peptide synthase
PKS biosynthesis locus
−3.0 - (−6.1)
Fatty aldehyde dehydrogenase
Amino acid transporter
Deaminating oxidase subunit
Deaminating oxidase subunit
Deaminating oxidase subunit
Inner membrane fatty acid desaturase
Acyl carrier protein
Inner membrane fatty acid desaturase
Periplasmic ligand binding protein
Aromatic hydrocarbons catabolism-related reductase
Aromatic hydrocarbons catabolism-related dioxygenase
Cytochrome bo oxidase subunit
Cytochrome bo oxidase subunit
AraC-type regulator, adjacent to polyketide genes
Glycine betaine aldehyde dehydrogenase
Intriguingly, genes encoding the T3SS3 apparatus components were found to be repressed in the wildtype compared with the ΔbsaN mutant, suggesting a role for BsaN in limiting apparatus synthesis when translocon and effector genes are transcribed (Figure 1A, 1E, Table 2). Also repressed are polar flagellar motility loci on chromosome 1 including the flagellin genes fliC and fliD, as well as flagellar hook proteins flgL and flgK. Repression of these genes as well as motA (BPSL3309) and cheD (BPSS3302) were validated by qRT-PCR (Figure 1E). In Salmonella and other bacteria, motAB are key components of the flagellar motor complex . motAB in KHW are part of a chemotaxis (che) locus, which is repressed 2–2.9-fold (p < 0.01) as assessed by RNAseq. In addition, expression of a second polyketide biosynthesis locus (BPSS0303-BPSS0311) was reduced in a ΔbsaN mutant, possibly by repression of a co-localized araC-type regulatory gene, BPSS0336 (Table 2). However, down-regulation of this cluster could not be verified by qRT-PCR (data not shown). We were likewise unable to validate repression of BPSL2404-2405, which putatively encode transport and energy metabolism functions, respectively, in addition to BPSS1887-1888, which are postulated to encode oxidative enzymes for energy metabolism.
Additional loci implicated in lipid and energy metabolism are also repressed (Table 2). Catabolic genes encode a cytochrome o oxidase typically used by bacteria in an oxygen-rich environment , along with enzymes involved in the aerobic degradation of aromatic compounds and in the degradation of arginine. A gene involved in the synthesis of betaine, an osmoprotectant which serves to adapt Gram-negative bacteria to conditions of high osmolarity, is also repressed by BsaN.
BsaN together with chaperone BicA directly activate T3SS3 effector and T6SS1 regulatory genes
We have previously shown that expression of the two component regulatory system virAG and the genes from BPSS1520 (bprC) to BPSS1533 (bicA) in the T3SS3 cluster were regulated by BsaN in concert with the chaperone BicA . To determine whether BsaN/BicA activate these genes directly, bsaN and bicA open reading frames (orfs) from B. pseudomallei strain KHW were inserted into a plasmid downstream of an arabinose-inducible promoter on pMLBAD . These constructs were introduced into E. coli DH5α  along with an additional construct containing putative promoter regions of several BsaN target genes transcriptionally fused to lacZ on pRW50  or pRW50mob, which contains the oriT fragment for pOT182 . The effect of BsaN/BicA on promoter activity was then assessed by β-galactosidase activities.
Identification of transcriptional start sites and the sequence motif for BsaN/BicA activation
tssM is one of the highly activated genes in our RNAseq analysis (Table 1) confirming previous in vivo expression studies . However, despite the presence of the BsaN box upstream of the putative tssM operon (BPSS1512-1514), BsaN/BicA alone is not sufficient to activate tssM transcription in E. coli (Figure 3G). This suggests that tssM regulation is more complex and likely requires additional cis and/or trans-acting regulatory elements for activation.
Determining the sequence motif requirement for BsaN/BicA activation
BprP directly activates bsaN and bsaM
Analysis of BsaN/BicA-regulated virulence loci
BsaN/BicA directly induces the expression of three known T3SS3 effector loci, bopA, bopC and bopE. Recent studies suggest that the T3SS3 effectors BopC and BopE are involved in invasion of epithelial cells and endosome escape ,,, while BopA has been implicated in escape from autophagy . BopC was recently shown to be secreted via T3SS3 in B. pseudomallei K96243 , and our data confirm this (Additional file 1: Figure S1). In B. pseudomallei KHW, mutation of bopA, bopC or bopE individually resulted in no detectable difference in numbers of bacteria inside RAW264.7 mouse macrophages when measured 2 hr. after infection (Additional file 1: Figure S2A). Upon extended incubation times, however, the ΔbopA and the ΔbopACE strains exhibited an intracellular replication defect that was intermediate between levels observed for wildtype KHW and the ΔbsaM or ΔbsaN mutant derivatives. No differences in intracellular growth or host cell cytotoxicity were observed for the bopC or bopE mutant strains, although infection with the bopA or bopACE triple deletion mutants resulted in a decrease in cytotoxicity (Additional file 1: Figure S2B) that coincided with a reduction in the rate of intracellular replication (Additional file 1: Figure S2A), suggesting that intracellular replication results in host cell toxicity. This is in contrast to the T3SS3 ΔbsaM and the ΔbsaN regulatory mutants in strain KHW, which are limited in their ability to multiply intracellularly as previously reported (Additional file 1: Figure S2A).
Three BsaN/BicA-activated orfs are located between the T3SS3 and T6SS1 loci, and upstream of the T3SS3 effector gene bopC. We analyzed these orfs for potential roles in intracellular replication and cell-to-cell spread. BPSS1512 encodes TssM, was previously shown to be secreted independently of T3SS3 and T6SS1 and functions as a broad-base deubiquitinase, with activity on TNFR-associated factor-3, TNFR-associated factor-6, and IκBα . BPSS1513 is predicted to encode a short (97 aa) protein of unknown function and was not secreted under our assay conditions (Additional file 1: Figure S3A). folE (BPSS1514) encodes a putative GTP cyclohydrolase I, suggesting a role in tetrahydrofolate biosynthesis rather than in virulence. Consistent with this notion, Δ(BPSS1513-folE) mutant did not exhibit defects in cell-based virulence assays (Additional file 1: Figure S3B-E).
T3SSs and T6SSs play important roles in bacterial-host cell interactions ,. As each system is a complex structure encoded by 20 or more genes, it is expected that their expression and assembly would be tightly regulated. In B. pseudomallei, T3SS3 and T6SS1 gene clusters are highly induced following host cell infection , and their function is critical for virulence in animal models ,. T3SS3 has been shown to promote escape from endocytic vesicles, and T6SS1 plays a key role in promoting intercellular spread by fusion of adjacent cell membranes, leading to the formation of MNGCs that can be found in melioidosis patients . Upregulation of T3SS3 is mediated by a signalling cascade initiating from BspR through BprP, which in turn increases the expression of the AraC-type regulator BsaN. Using global transcriptome and promoter activation analysis, we have shown that the BsaN regulon occupies a central position in modulating the expression of T3SS3, T6SS1 and several additional loci that are likely involved in promoting virulence and intracellular survival.
Regulatory factors may act to control expression by acting directly on a given gene, or indirectly by modulating a regulatory intermediate. We found that BsaN in complex with the T3SS3 chaperone BicA directly controls the expression of 19 loci in a region of chromosome 2 containing T6SS1 and T3SS3 accessory genes (BPSS1494-BPSS1533). BsaN/BicA activated transcription of the operons encoding T3SS3 effector proteins, the BipBCD translocon complex, chaperones, and other transcriptional regulators, as well as two genes of unknown function (BPSS1513-1514). BsaN/BicA upregulates expression of T6SS1 by activating the transcription of the two component regulatory system loci virAG and bprC, which in turn induce the hcp and tssAB loci, encoding T6SS1 tube and sheath proteins ,. Interestingly, our RNAseq and qRT-PCR analyses revealed that BsaN also acts to repress transcription of T3SS3 apparatus genes in the bsaM and bsaN operons that are otherwise directly activated by the upstream regulator BprP. It is possible that BsaN mediates repression indirectly as the bsaM and bsaN intergenic region lacks a recognizable BsaN binding motif (see below). It is unlikely, however, that repression occurs due to decreased expression of bprP since its transcription is unchanged in a ΔbsaN mutant. Taken together, these findings demonstrate that BsaN plays a dual role in the regulation of T3SS3; one in coordinating translocon and effector transcription, and a second in preventing costly synthesis of T3SS3 apparatus components that are no longer required. Given the critical role of T3SS3 and T6SS1 in causing disease, BsaN/BicA could be considered a central regulator of B. pseudomallei mammalian virulence. Virulence studies in mice support this notion, since the ΔbsaN mutant was unable to cause disease  in contrast to the ΔbspR mutant, which produced a more chronic infection in mice compared to wildtype bacteria .
In addition to loci associated with T3SS3 and T6SS1, 41 other genes with potential roles in virulence were also found by RNAseq to be positively regulated by BsaN, most notably the bimBCAD intracellular motility operon and tssM. Regulation of bimA has been shown to be through virAG, explaining why no BsaN motif was identified for the operon. While bimA encodes an autotransporter protein that nucleates and polymerizes host cell actin to facilitate intracellular motility and cell-cell spread by the bacteria , the functions of the other loci in the bim operon are unknown. TssM has been shown to suppress host NFκB and Type I interferon pathways . TssM is expressed and secreted inside cells following infection with B. mallei, however, secretion occurs independently of T3SS3 and T6SS1 . BsaN was also found to activate expression of a putative non-ribosomal peptide synthase (NRPS)/polyketide synthase (PKS) biosynthesis locus. The diversity of polyketides, PKSs and NRPS/PKS hybrid systems was recently reviewed by Hertweck . The B. pseudomallei locus is similar in gene content to that of a recently described plasmid encoded NRPS/PKS system in the marine bacterium Alteromonas macleodii, which was suggested to produce a bleomycin-related antibiotic Unlike A. macleodii, the gene encoding the putative bleomycin-family resistance protein (BPSL2883) is not co-localized with the NRPS/PKS gene cluster, although they are similarly regulated by BsaN (Table 1).
BsaN is homologous to the Salmonella typhimurium InvF, Shigella flexneri MxiE and the Yersinia enterocholitica YsaB transcriptional regulators –. All belong to the AraC/XylS family of transcriptional regulators, which act in complex with a chaperone to activate their respective T3SS genes. The chaperones not only serve as cognate partners to the transcriptional activators but also pair with T3SS translocase proteins, which are secreted into the host membrane to facilitate the injection of effector proteins . We currently, have no understanding of the timed mechanism that frees BicA and allows it to partner with BsaN. The S. typhimurium chaperone SicA was shown to partition the translocase SipB and SipC, and it is sequestered by SipB . Once apparatus assembly is complete, translocases are secreted and SicA is free to complex and thus activate InvF. The InvF-SicA split feedback regulatory loop, which includes positive autoregulation of invF, is conserved in Y. enterocholitica. However, in S. flexneri MxiE-dependent activity is inhibited via sequestration by the T3SS substrate OspD1 when the apparatus is inactive . Only when OspD1 is secreted, can MxiE partner with its chaperone IpgC to activated transcription of effector genes.
Until recently, BopA and BopE were the only two known T3SS effector proteins in B. pseudomallei. The dearth of effectors is surprising when compared to other intracellular pathogens such as Shigella and Salmonella that are known to possess numerous effectors. We have independently identified BopC (BPSS1516) as a new T3SS3 effector based on its regulatory control by BsaN/BicA. bopC is transcribed in an operon encoding its chaperone (BPSS1517) and a transposase (BPSS1518) that are also activated by BsaN/BicA. Incidentally, we had previously predicted by a genome-wide screen that BPSS1516 would encode a T3SS effector based on genomic colocalization with T3SS chaperones . The BsaN regulatory motif we found in the promoters of the effectors was also recently reported to be associated with T3SS3 in a condition-dependent transcriptome study . Of the T3SS3-linked effector proteins; BopA, BopC and BopE, our results suggest that BopA is the most critical for promoting cellular infection, consistent with prior studies linking BopA to intracellular survival of B. pseudomallei and B. mallei,,. No cellular phenotype was evident following infection with ΔbopC or ΔbopE deletion mutants, and the ΔbopACE triple effector mutant was indistinguishable from the ΔbopA single deletion strain. As with bopE and bopC, no roles were observed for the BsaN-regulated effector candidate loci BPSS1513-1514 in cell-based virulence assays. BPSS1513 encodes a hypothetical protein and BPSS1514 is annotated as folE, a predicted GTP cyclohydrolase. Based on their genomic organization, the transcription of these loci is likely driven from the promoter upstream of BPSS1512 tssM. The secretion of HA-tagged BPSS1513 was not detected in in vitro secretion assays, although it is possible that the epitope tag could have interfered with secretion of BPSS1513, or that the assay was not performed at conditions permissive for secretion. It is intriguing why these three genes are placed under BsaN/BicA regulation by the bacterium. One possibility could be that they are important under specific stress conditions or during chronic infection.
Elucidating the scope of the BsaN regulon significantly enhances our understanding of B. pseudomallei pathogenic mechanisms. BsaN orchestrates the temporal and spatial expression of virulence determinants during progression through the intracellular lifecycle, promoting endosome escape and possibly evasion of autophagy through activation of T3SS3 effector loci, facilitating cell-cell spread by activation of T6SS1 and the bim intracellular motility loci, and suppressing cellular immunity via the action of the TssM ubiquitin hydrolase. BsaN also suppresses other loci that are potentially counterproductive following intracellular localization, such as the fla1 flagellar motility and chemotaxis locus, which could lead to activation of cellular immunity pathways through PAMP recognition. It is likely that the BsaN regulon and other virulence determinants that promote pathogenesis in higher mammals have been shaped primarily as a result of interactions with free-living protozoa, similar to what is believed to be the case for L. pneumophila. Indeed, many of the same BsaN-regulated systems, namely T3SS and T6SS, are thought to act as “anti-predation determinants” that facilitate endosome escape and promote survival within bacteriovorus amoebae by manipulating eukaryotic pathways that are conserved from protists to humans . The dual regulatory roles of BsaN – that of an activator and a suppressor – indicate that it is a key node in a regulatory program that successfully enables an environmental saprophyte to transition from the soil to surviving intracellularly.
Bacterial strains and culture conditions
List of strains used in this study
Relevant characteristic(s) a
Source or reference
Donor strain for conjugation
KHW, bsaN orf was deleted
KHW, bsaM orf was deleted
KHW, bopA orf was replaced by zeocin resistance gene, Zeor
KHW, bopE orf was deleted
KHW, bopC orf was replaced by zeocin resistance gene, Zeor
KHW, bopA and bopE orfs were deleted, bopC orf was replaced by zeocin resistance gene, Zeor
KHW, BPSS1513 and folE orfs were deleted
List of plasmids used in this study
Relevant characteristic(s) a
Source or reference
Broad host range vector containing inducible promoter Pbad, Tmr
Low copy-number vector containing a promoter-less lacZ gene, Tcr
Source of oriT sequence, Tcr
Conjugative, suicide vector containing sacB gene, Kmr
pRW50 containing oriT sequence from pOT182, Tcr
pMLBAD containing bspR orf from KHW, Tmr
pMLBAD containing bprP orf from KHW, Tmr
pMLBAD containing bprQ orf from KHW, Tmr
pMLBAD containing bprP and bprQ orfs from KHW, Tmr
pMLBAD containing shorter bsaN orf (GTG start) and bicA orf from KHW, Tmr
pMLBAD containing longer bsaN orf (ATG start) from KHW, Tmr
pMLBAD containing longer bsaN orf (ATG start) and bicA orf from KHW, Tmr
PbsaM-lacZ transcriptional fusion, pRW50mob containing bsaM upstream sequence from KHW, Tcr
PbsaN-lacZ transcriptional fusion, pRW50mob containing bsaN upstream sequence from KHW, Tcr
PbicA-lacZ transcriptional fusion, pRW50mob containing bicA upstream sequence from KHW, Tcr
PbicA1-lacZ transcriptional fusion, pRW50mob containing truncated bicA upstream sequence from KHW, type 1 deletion, Tcr
PbicA2-lacZ transcriptional fusion, pRW50mob containing truncated bicA upstream sequence from KHW, type 2 deletion, Tcr
PbopE-lacZ transcriptional fusion, pRW50mob containing bopE upstream sequence from KHW, Tcr
PbopA-lacZ transcriptional fusion, pRW50mob containing bopA upstream sequence from KHW, Tcr
PbopA1-lacZ transcriptional fusion, pRW50mob containing truncated bopA upstream sequence from KHW, type 1 deletion, Tcr
PbopA2-lacZ transcriptional fusion, pRW50mob containing truncated bopA upstream sequence from KHW, type 2 deletion, Tcr
PbprC-lacZ transcriptional fusion, pRW50mob containing bprB upstream sequence from KHW, Tcr
PbprD-lacZ transcriptional fusion, pRW50mob containing bprD upstream sequence from KHW, Tcr
PbprD1-lacZ transcriptional fusion, pRW50mob containing truncated bprD upstream sequence from KHW, type 1 deletion, Tcr
PbicA2-lacZ transcriptional fusion, pRW50mob containing truncated bicA upstream sequence from KHW, type 2 deletion, Tcr
PbprC-lacZ transcriptional fusion, pRW50mob containing bprC upstream sequence from KHW, Tcr
Ps1518-lacZ transcriptional fusion, pRW50 containing BPSS1518 upstream sequence from KHW, Tcr
Ps1518a-lacZ transcriptional fusion, pRW50 containing truncated BPSS1518 upstream sequence from KHW, type 1 deletion, Tcr
Ps1518b-lacZ transcriptional fusion, pRW50 containing truncated BPSS1518 upstream sequence from KHW, type 2 deletion, Tcr
PvirA-lacZ transcriptional fusion, pRW50 containing virA upstream sequence from KHW, Tcr
PvirA1-lacZ transcriptional fusion, pRW50 containing truncated virA upstream sequence from KHW, type 1 deletion, Tcr
PtssM-lacZ transcriptional fusion, pRW50 containing tssM upstream sequence from KHW, Tcr
Bacterial mutant construction
B. pseudomallei gene deletions were generated by allelic exchange. Approximately 1 kb fragments upstream and downstream of the target gene were amplified from genomic DNA and cloned into pK18mobsacB vector  simultaneously using In-Fusion PCR cloning kit (Clontech). The plasmids were introduced into B. pseudomallei strains by conjugation. Homologous recombination was then selected for by growing bacteria in LB + 15% sucrose to counter select the sacB gene in the pK18mobsacB plasmid backbone. Successful double cross-over clones were screened by colony PCR from kanamycin sensitive colonies.
Activation of potential promoters by regulators
The ability of regulators to directly activate the expression of promoters was examined in E. coli DH5α as described previously . Briefly, upstream regions of B. pseudomallei genes encompassing at least 100 bp of non-coding sequence upstream of the start codon were amplified from KHW genomic DNA and fused to the lacZ gene in pRW50 or pRW50mob to generate transcriptional fusions (Table 4). Coding sequences of regulators were amplified from KHW genomic DNA and cloned into the arabinose-inducible expression vector pMLBAD. The lacZ fusion plasmid and arabinose-inducible regulator plasmid were introduced into the E. coli DH5α. β-galactosidase activities arising from the expression of promoter-lacZ fusions were assessed. β-Galactosidase assays were performed and values were calculated as previously described .
Transcriptome analysis by RNAseq
Total RNA was extracted from three independently grown bacterial cultures that were combined at equal cell density in their exponential growth phase and quick frozen in dry ice-ethanol slurry. Approximately 2 × 109 ice cold cells were centrifuged at 3000 × g for 45 sec and 4°C and RNA was isolated from cell pellets using the RiboPure™-Bacteria Kit (Ambion). Stable RNAs were removed from 10 μg RNA using the MICROBExpress kit from Ambion. Absence of genomic DNA contamination was confirmed by PCR. Paired-end libraries for Illumina sequencing  were prepared using the TruSeq RNA sample preparation kit version 2.0 (Illumina) according to manufacturer's High Sample (HS) protocol albeit omitting the initial poly A selection step. Libraries were generated from 2 technical replicates using 350–500 ng enriched RNA from wildtype and ΔbsaN mutant strains as the starting material. Library preparation and sequencing was done by the UCLA Neuroscience Genomics Core (UNGC). Reads were aligned to chromosomes I and II of B. pseudomallei KHW (also called BP22) (RefSeq identification numbers NZ_CM001156.1 and NZ_CM001157.1) and B. pseudomallei K96243 (RefSeq identification numbers NC_006350.1 and NC_006351.1) as the annotated reference genome. The number of reads aligning to each genomic position on each strand was calculated and normalized using RPKM ([reads/kb of gene]/[million reads aligning to genome]). Differentially expressed genes identified by the log2 ratio of the differential between the wildtype and ΔbsaN RPKMs. Only, genes with a Δlog2 value of >1.5 and < −1.5 corresponding to 3-fold up or down regulated genes with an adjusted p value (padj) of <0.01 were considered for this study.
Measurement of B. pseudomallei gene expression by qRT- PCR
Expression of activated genes was confirmed by qRT-PCR of RNA prepared from bacteria grown in acidified RPMI. Gene repression was difficult to observe under these conditions; RNA for qRT-PCR analysis was therefore prepared from infected RAW264.7 cells using the following procedure: RAW264.7 cells (5 × 105 cells/well) were seeded and grown overnight in DMEM medium in 12 well plates. RAW264.7 cells were transferred to RPMI medium prior to infection and infected at MOI of 100:1. Bacterial RNA was isolated from infected RAW264.7 cells 4 hours post infection using TRIzol and PureLink RNA mini-kit (Invitrogen). cDNA was synthesized using 1 μg of RNA and the High Capacity Reverse Transcription Reagent Kit (Applied Biosystems). Transcripts were quantified using GoTaq qPCR Master Mix (Promega) in a BioRad iQ5 machine. Real-time PCR primers are listed in Additional file 1: Table S4. Relative RNA level of a particular gene in mutant strains was normalized to that of wild type using the 2−ΔΔCt method with 16S rRNA or recA as reference gene .
Mapping transcriptional start sites
The transcriptional start (+1) sites of the promoters were mapped by RNA ligase-mediated rapid amplification of cDNA ends (RLM-RACE) . The RLM-RACE was performed using GeneRacer Kit (Invitrogen) according to manufacturer’s instructions. The B. pseudomallei RNA was isolated as described previously .
Sequence motif predication and database search of motifs in the B. pseudomallei KHW genome
150 base pairs of nucleotide sequence upstream of the transcriptional starts of each gene was submitted to the bioinformatics tool – MEME (http://meme.nbcr.net/meme/cgi-bin/meme.cgi) for prediction of DNA motifs . The motif with the highest statistical significance (lowest E-value) was chosen and its data – in Position-Specific Probability Matrix format was submitted to MAST (http://meme.nbcr.net/meme/cgi-bin/mast.cgi) to search for the best matching positions in the upstream sequences of B. pseudomallei KHW genes .
Results were presented as mean ± standard deviation. Student’s t-test was used to find the significant differences between the means, defined as when p < 0.05 (*) and p < 0.01 (**).
We thank M. A. Valvano (University of Western Ontario) for pMLBAD plasmid, S.J. Busby (University of Birmingham) for pRW50 plasmid, S. Korbsrisate (Mahidol University) for BopC antibody, and M.P. Stevens (University of Edinburgh) for BopE antibody. This work is supported by grants T208A3105 from the Ministry of Education to YHG, NMRC/1221/2009 from the National Medical Research Council to YHG, an award from the Pacific Southwest Regional Center of Excellence in Biodefense and Emerging Infectious Diseases (NIH U54 A1065359) to JFM, and grant HDTRA1-11-1-0003 from the Defense Threat Reduction Agency to JFM. We would like to thank Isabelle Chen for her technical assistance.
- Galyov EE, Brett PJ, DeShazer D: Molecular insights into Burkholderia pseudomallei and Burkholderia mallei pathogenesis. Annu Rev Microbiol. 2010, 64: 495-517. 10.1146/annurev.micro.112408.134030.View ArticlePubMed
- Wiersinga WJ, van der Poll T, White NJ, Day NP, Peacock SJ: Melioidosis: insights into the pathogenicity of Burkholderia pseudomallei. Nat Rev Microbiol. 2006, 4 (4): 272-282. 10.1038/nrmicro1385.View ArticlePubMed
- Hasselbring BM, Patel MK, Schell MA: Dictyostelium discoideum as a model system for identification of Burkholderia pseudomallei virulence factors. Infect Immun. 2011, 79 (5): 2079-2088. 10.1128/IAI.01233-10.PubMed CentralView ArticlePubMed
- Inglis TJ, Rigby P, Robertson TA, Dutton NS, Henderson M, Chang BJ: Interaction between Burkholderia pseudomallei and Acanthamoeba species results in coiling phagocytosis, endamebic bacterial survival, and escape. Infect Immun. 2000, 68 (3): 1681-1686. 10.1128/IAI.68.3.1681-1686.2000.PubMed CentralView ArticlePubMed
- Lee YH, Chen Y, Ouyang X, Gan YH: Identification of tomato plant as a novel host model for Burkholderia pseudomallei. BMC Microbiol. 2010, 10: 28-10.1186/1471-2180-10-28.PubMed CentralView ArticlePubMed
- Kaestli M, Schmid M, Mayo M, Rothballer M, Harrington G, Richardson L, Hill A, Hill J, Tuanyok A, Keim P, Hartmann A, Currie BJ: Out of the ground: aerial and exotic habitats of the melioidosis bacterium Burkholderia pseudomallei in grasses in Australia. Environ Microbiol. 2012, 14 (8): 2058-2070. 10.1111/j.1462-2920.2011.02671.x.PubMed CentralView ArticlePubMed
- Burtnick MN, Brett PJ, Harding SV, Ngugi SA, Ribot WJ, Chantratita N, Scorpio A, Milne TS, Dean RE, Fritz DL, Peacock SJ, Prior JL, Atkins TP, Deshazer D: The Cluster 1 Type VI Secretion System is a Major Virulence Determinant in Burkholderia pseudomallei. Infect Immun. 2011, 79 (4): 1512-1525. 10.1128/IAI.01218-10.PubMed CentralView ArticlePubMed
- Chen Y, Wong J, Sun GW, Liu Y, Tan GY, Gan YH: Regulation of type VI secretion system during Burkholderia pseudomallei infection. Infect Immun. 2011, 79 (8): 3064-3073. 10.1128/IAI.05148-11.PubMed CentralView ArticlePubMed
- French CT, Toesca IJ, Wu TH, Teslaa T, Beaty SM, Wong W, Liu M, Schröder I, Chiou PY, Teitell MA, Miller JF: Dissection of the Burkholderia intracellular life cycle using a photothermal nanoblade. Proc Natl Acad Sci U S A. 2011, 108 (29): 12095-12100. 10.1073/pnas.1107183108.PubMed CentralView ArticlePubMed
- Stevens MP, Wood MW, Taylor LA, Monaghan P, Hawes P, Jones PW, Wallis TS, Galyov EE: An Inv/Mxi-Spa-like type III protein secretion system in Burkholderia pseudomallei modulates intracellular behaviour of the pathogen. Mol Microbiol. 2002, 46 (3): 649-659. 10.1046/j.1365-2958.2002.03190.x.View ArticlePubMed
- Sun GW, Lu J, Pervaiz S, Cao WP, Gan YH: Caspase-1 dependent macrophage death induced by Burkholderia pseudomallei. Cell Microbiol. 2005, 7 (10): 1447-1458. 10.1111/j.1462-5822.2005.00569.x.View ArticlePubMed
- Stevens MP, Haque A, Atkins T, Hill J, Wood MW, Easton A, Nelson M, Underwood-Fowler C, Titball RW, Bancroft GJ, Galyov EE: Attenuated virulence and protective efficacy of a Burkholderia pseudomallei bsa type III secretion mutant in murine models of melioidosis. Microbiology. 2004, 150 (Pt 8): 2669-2676. 10.1099/mic.0.27146-0.View ArticlePubMed
- Warawa J, Woods DE: Type III secretion system cluster 3 is required for maximal virulence of Burkholderia pseudomallei in a hamster infection model. FEMS Microbiol Lett. 2005, 242 (1): 101-108. 10.1016/j.femsle.2004.10.045.View ArticlePubMed
- Sun GW, Chen Y, Liu Y, Tan GY, Ong C, Tan P, Gan YH: Identification of a regulatory cascade controlling Type III Secretion System 3 gene expression in Burkholderia pseudomallei. Mol Microbiol. 2010, 76 (3): 677-689. 10.1111/j.1365-2958.2010.07124.x.View ArticlePubMed
- Stevens MP, Friebel A, Taylor LA, Wood MW, Brown PJ, Hardt WD, Galyov EE: A Burkholderia pseudomallei type III secreted protein, BopE, facilitates bacterial invasion of epithelial cells and exhibits guanine nucleotide exchange factor activity. J Bacteriol. 2003, 185 (16): 4992-4996. 10.1128/JB.185.16.4992-4996.2003.PubMed CentralView ArticlePubMed
- Cullinane M, Gong L, Li X, Lazar-Adler N, Tra T, Wolvetang E, Prescott M, Boyce JD, Devenish RJ, Adler B: Stimulation of autophagy suppresses the intracellular survival of Burkholderia pseudomallei in mammalian cell lines. Autophagy. 2008, 4 (6): 744-753. 10.4161/auto.6246.View ArticlePubMed
- Gong L, Cullinane M, Treerat P, Ramm G, Prescott M, Adler B, Boyce JD, Devenish RJ: The Burkholderia pseudomallei type III secretion system and BopA are required for evasion of LC3-associated phagocytosis. PLoS One. 2011, 6 (3): e17852-10.1371/journal.pone.0017852.PubMed CentralView ArticlePubMed
- Muangman S, Korbsrisate S, Muangsombut V, Srinon V, Adler NL, Schroeder GN, Frankel G, Galyov EE: BopC is a type III secreted effector protein of Burkholderia pseudomallei. FEMS Microbiol Lett. 2011, 323 (1): 75-82. 10.1111/j.1574-6968.2011.02359.x.View ArticlePubMed
- Srinon V, Muangman S, Imyaem N, Muangsombut V, Lazar Adler NR, Galyov EE, Korbsrisate S: Comparative assessment of the intracellular survival of the Burkholderia pseudomallei bopC mutant. J Microbiol. 2013, 51 (4): 522-526. 10.1007/s12275-013-2557-3.View ArticlePubMed
- Liu B, Koo GC, Yap EH, Chua KL, Gan YH: Model of differential susceptibility to mucosal Burkholderia pseudomallei infection. Infect Immun. 2002, 70 (2): 504-511. 10.1128/IAI.70.2.504-511.2002.PubMed CentralView ArticlePubMed
- Biggins JB, Liu X, Feng Z, Brady SF: Metabolites from the induced expression of cryptic single operons found in the genome of Burkholderia pseudomallei. J Am Chem Soc. 2011, 133 (6): 1638-1641. 10.1021/ja1087369.View ArticlePubMed
- Tang H, Braun TF, Blair DF: Motility protein complexes in the bacterial flagellar motor. J Mol Biol. 1996, 261 (2): 209-221. 10.1006/jmbi.1996.0453.View ArticlePubMed
- Cotter PA, Melville SB, Albrecht JA, Gunsalus RP: Aerobic regulation of cytochrome d oxidase (cydAB) operon expression in Escherichia coli: roles of Fnr and ArcA in repression and activation. Mol Microbiol. 1997, 25 (3): 605-615. 10.1046/j.1365-2958.1997.5031860.x.View ArticlePubMed
- Lefebre MD, Valvano MA: Construction and evaluation of plasmid vectors optimized for constitutive and regulated gene expression in Burkholderia cepacia complex isolates. Appl Environ Microbiol. 2002, 68 (12): 5956-5964. 10.1128/AEM.68.12.5956-5964.2002.PubMed CentralView ArticlePubMed
- Sambrook J, Maniatis T, Fritsch EF: Molecular cloning : a laboratory manual. 1989, Cold Spring Harbor Laboratory, Cold Spring Harbor, N.Y
- Lodge J, Fear J, Busby S, Gunasekaran P, Kamini NR: Broad host range plasmids carrying the Escherichia coli lactose and galactose operons. FEMS Microbiol Lett. 1992, 74 (2–3): 271-276. 10.1111/j.1574-6968.1992.tb05378.x.View ArticlePubMed
- DeShazer D, Brett PJ, Carlyon R, Woods DE: Mutagenesis of Burkholderia pseudomallei with Tn5-OT182: isolation of motility mutants and molecular characterization of the flagellin structural gene. J Bacteriol. 1997, 179 (7): 2116-2125.PubMed CentralPubMed
- Holden MT, Titball RW, Peacock SJ, Cerdeño-Tárraga AM, Atkins T, Crossman LC, Pitt T, Churcher C, Mungall K, Bentley SD, Sebaihia M, Thomson NR, Bason N, Beacham IR, Brooks K, Brown KA, Brown NF, Challis GL, Cherevach I, Chillingworth T, Cronin A, Crossett B, Davis P, DeShazer D, Feltwell T, Fraser A, Hance Z, Hauser H, Holroyd S, Jagels K: Genomic plasticity of the causative agent of melioidosis, Burkholderia pseudomallei. Proc Natl Acad Sci U S A. 2004, 101 (39): 14240-14245. 10.1073/pnas.0403302101.PubMed CentralView ArticlePubMed
- Shanks J, Burtnick MN, Brett PJ, Waag DM, Spurgers KB, Ribot WJ, Schell MA, Panchal RG, Gherardini FC, Wilkinson KD, Deshazer D: Burkholderia mallei tssM encodes a putative deubiquitinase that is secreted and expressed inside infected RAW 264.7 murine macrophages. Infect Immun. 2009, 77 (4): 1636-1648. 10.1128/IAI.01339-08.PubMed CentralView ArticlePubMed
- Teh BE, French CT, Chen Y, Chen IG, Wu TH, Sagullo E, Chiou PY, Teitell MA, Miller JF, Gan YH: Type three secretion system-mediated escape of Burkholderia pseudomallei into the host cytosol is critical for the activation of NFkappaB. BMC Microbiol. 2014, 14: 115-10.1186/1471-2180-14-115.PubMed CentralView ArticlePubMed
- Tan KS, Chen Y, Lim YC, Tan GY, Liu Y, Lim YT, Macary P, Gan YH: Suppression of host innate immune response by Burkholderia pseudomallei through the virulence factor TssM. J Immunol. 2010, 184 (9): 5160-5171. 10.4049/jimmunol.0902663.View ArticlePubMed
- Raymond B, Young JC, Pallett M, Endres RG, Clements A, Frankel G: Subversion of trafficking, apoptosis, and innate immunity by type III secretion system effectors. Trends Microbiol. 2013, 21 (8): 430-441. 10.1016/j.tim.2013.06.008.View ArticlePubMed
- Jani AJ, Cotter PA: Type VI secretion: not just for pathogenesis anymore. Cell Host Microbe. 2010, 8 (1): 2-6. 10.1016/j.chom.2010.06.012.PubMed CentralView ArticlePubMed
- Wong KT, Puthucheary SD, Vadivelu J: The histopathology of human melioidosis. Histopathology. 1995, 26 (1): 51-55. 10.1111/j.1365-2559.1995.tb00620.x.View ArticlePubMed
- Cascales E, Cambillau C: Structural biology of type VI secretion systems. Philos Trans R Soc Lond B Biol Sci. 2012, 367 (1592): 1102-1111. 10.1098/rstb.2011.0209.PubMed CentralView ArticlePubMed
- Stevens MP, Stevens JM, Jeng RL, Taylor LA, Wood MW, Hawes P, Monaghan P, Welch MD, Galyov EE: Identification of a bacterial factor required for actin-based motility of Burkholderia pseudomallei. Mol Microbiol. 2005, 56 (1): 40-53. 10.1111/j.1365-2958.2004.04528.x.View ArticlePubMed
- Hertweck C: The biosynthetic logic of polyketide diversity. Angew Chem Int Ed Engl. 2009, 48 (26): 4688-4716. 10.1002/anie.200806121.View ArticlePubMed
- Darwin KH, Miller VL: Type III secretion chaperone-dependent regulation: activation of virulence genes by SicA and InvF in Salmonella typhimurium. EMBO J. 2001, 20 (8): 1850-1862. 10.1093/emboj/20.8.1850.PubMed CentralView ArticlePubMed
- Kane CD, Schuch R, Day WA, Maurelli AT: MxiE regulates intracellular expression of factors secreted by the Shigella flexneri 2a type III secretion system. J Bacteriol. 2002, 184 (16): 4409-4419. 10.1128/JB.184.16.4409-4419.2002.PubMed CentralView ArticlePubMed
- Walker KA, Miller VL: Regulation of the Ysa type III secretion system of Yersinia enterocolitica by YsaE/SycB and YsrS/YsrR. J Bacteriol. 2004, 186 (13): 4056-4066. 10.1128/JB.186.13.4056-4066.2004.PubMed CentralView ArticlePubMed
- Deane JE, Abrusci P, Johnson S, Lea SM: Timing is everything: the regulation of type III secretion. Cell Mol Life Sci. 2010, 67 (7): 1065-1075. 10.1007/s00018-009-0230-0.PubMed CentralView ArticlePubMed
- Tucker SC, Galan JE: Complex function for SicA, a Salmonella enterica serovar typhimurium type III secretion-associated chaperone. J Bacteriol. 2000, 182 (8): 2262-2268. 10.1128/JB.182.8.2262-2268.2000.PubMed CentralView ArticlePubMed
- Parsot C, Ageron E, Penno C, Mavris M, Jamoussi K, d'Hauteville H, Sansonetti P, Demers B: A secreted anti-activator, OspD1, and its chaperone, Spa15, are involved in the control of transcription by the type III secretion apparatus activity in Shigella flexneri. Mol Microbiol. 2005, 56 (6): 1627-1635. 10.1111/j.1365-2958.2005.04645.x.View ArticlePubMed
- Tuanyok A, Auerbach RK, Brettin TS, Bruce DC, Munk AC, Detter JC, Pearson T, Hornstra H, Sermswan RW, Wuthiekanun V, Peacock SJ, Currie BJ, Keim P, Wagner DM: A horizontal gene transfer event defines two distinct groups within Burkholderia pseudomallei that have dissimilar geographic distributions. J Bacteriol. 2007, 189 (24): 9044-9049. 10.1128/JB.01264-07.PubMed CentralView ArticlePubMed
- Biggins JB, Ternei MA, Brady SF: Malleilactone, a polyketide synthase-derived virulence factor encoded by the cryptic secondary metabolome of Burkholderia pseudomallei group pathogens. J Am Chem Soc. 2012, 134 (32): 13192-13195. 10.1021/ja3052156.PubMed CentralView ArticlePubMed
- Lamont IL, Beare PA, Ochsner U, Vasil AI, Vasil ML: Siderophore-mediated signaling regulates virulence factor production in Pseudomonasaeruginosa. Proc Natl Acad Sci U S A. 2002, 99 (10): 7072-7077. 10.1073/pnas.092016999.PubMed CentralView ArticlePubMed
- Panina EM, Mattoo S, Griffith N, Kozak NA, Yuk MH, Miller JF: A genome-wide screen identifies a Bordetella type III secretion effector and candidate effectors in other species. Mol Microbiol. 2005, 58 (1): 267-279. 10.1111/j.1365-2958.2005.04823.x.View ArticlePubMed
- Ooi WF, Ong C, Nandi T, Kreisberg JF, Chua HH, Sun G, Chen Y, Mueller C, Conejero L, Eshaghi M, Ang RM, Liu J, Sobral BW, Korbsrisate S, Gan YH, Titball RW, Bancroft GJ, Valade E, Tan P: The condition-dependent transcriptional landscape of Burkholderia pseudomallei. PLoS Genet. 2013, 9 (9): e1003795-10.1371/journal.pgen.1003795.PubMed CentralView ArticlePubMed
- Whitlock GC, Estes DM, Young GM, Young B, Torres AG: Construction of a reporter system to study Burkholderia mallei type III secretion and identification of the BopA effector protein function in intracellular survival. Trans R Soc Trop Med Hyg. 2008, 102 (Suppl 1): S127-133. 10.1016/S0035-9203(08)70029-4.View ArticlePubMed
- Escoll P, Rolando M, Gomez-Valero L, Buchrieser C: From Amoeba to Macrophages: Exploring the Molecular Mechanisms of Legionella pneumophila Infection in Both Hosts. Curr Top Microbiol Immunol. 2013, 376: 1-34.PubMed
- Simon R, Priefer U, Pühler A: A broad range mobilization system for in vitro genetic engineering: Transposon mutagenesis in Gram-negative bacteria. Bio/Technology. 1983, 1: 784-791. 10.1038/nbt1183-784.View Article
- Schafer A, Tauch A, Jager W, Kalinowski J, Thierbach G, Puhler A: Small mobilizable multi-purpose cloning vectors derived from the Escherichia coli plasmids pK18 and pK19: selection of defined deletions in the chromosome of Corynebacterium glutamicum. Gene. 1994, 145 (1): 69-73. 10.1016/0378-1119(94)90324-7.View ArticlePubMed
- Zhang X, Bremer H: Control of the Escherichia coli rrnB P1 promoter strength by ppGpp. J Biol Chem. 1995, 270 (19): 11181-11189. 10.1074/jbc.270.19.11181.View ArticlePubMed
- Levin JZ, Yassour M, Adiconis X, Nusbaum C, Thompson DA, Friedman N, Gnirke A, Regev A: Comprehensive comparative analysis of strand-specific RNA sequencing methods. Nat Methods. 2010, 7 (9): 709-715. 10.1038/nmeth.1491.PubMed CentralView ArticlePubMed
- Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001, 25 (4): 402-408. 10.1006/meth.2001.1262.View ArticlePubMed
- Bensing BA, Meyer BJ, Dunny GM: Sensitive detection of bacterial transcription initiation sites and differentiation from RNA processing sites in the pheromone-induced plasmid transfer system of Enterococcus faecalis. Proc Natl Acad Sci U S A. 1996, 93 (15): 7794-7799. 10.1073/pnas.93.15.7794.PubMed CentralView ArticlePubMed
- Bailey TL, Elkan C: Fitting a mixture model by expectation maximization to discover motifs in biopolymers. Proc Int Conf Intell Syst Mol Biol. 1994, 2: 28-36.PubMed
- Bailey TL, Gribskov M: Combining evidence using p-values: application to sequence homology searches. Bioinformatics. 1998, 14 (1): 48-54. 10.1093/bioinformatics/14.1.48.View ArticlePubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.