Construction of a mariner-based transposon vector for use in insertion sequence mutagenesis in selected members of the Rhizobiaceae
© Perry and Yost; licensee BioMed Central Ltd. 2014
Received: 13 July 2014
Accepted: 18 November 2014
Published: 30 November 2014
The Rhizobiaceae family of Gram-negative bacteria often engage in symbiosis with plants of economic importance. Historically, genetic studies to identify the function of individual genes, and characterize the biology of these bacteria have relied on the use of classical transposon mutagenesis. To increase the rate of scientific discovery in the Rhizobiaceae there is a need to adapt high-throughput genetic screens like insertion sequencing for use in this family of bacteria. Here we describe a Rhizobiaceae compatible MmeI-adapted mariner transposon that can be used with insertion sequencing for high-throughput genetic screening.
The newly constructed mariner transposon pSAM_Rl mutagenized R. leguminosarum, S. meliloti, and A. tumefaciens at a high frequency. In R. leguminosarum, mutant pools were generated that saturated 88% of potential mariner insertions sites in the genome. Analysis of the R. leguminosarum transposon insertion sequencing data with a previously described hidden Markov model-based method resulted in assignment of the contribution of all annotated genes in the R. leguminosarum 3841 genome for growth on a complex medium. Good concordance was observed between genes observed to be required for growth on the complex medium, and previous studies.
The newly described Rhizobiaceaee compatible mariner transposon insertion sequencing vector pSAM_Rl has been shown to mutagenize at a high frequency and to be an effective tool for use in high-throughput genetic screening. The construction and validation of this transposon insertion sequencing tool for use in the Rhizobiziaceae will provide an opportunity for researchers in the Rhizobiaceae community to use high-throughput genetic screening, allowing for significant increase in the rate of genetic discovery, particularly given the recent release of genome sequences from many Rhizobiaceae strains.
Insertion sequencing (INSeq) is a technique for high throughput forward genetic screening that has recently become a favorable approach to studying gene function at the genome scale ,. INSeq relies on the use of next-generation DNA sequencing to audit the presence of hundreds of thousands of unique transposon insertions present in a pool of mutants that collectively saturate that organism’s genome with transposition events –. In general, INSeq based methods can use two different methods to analyze gene function. The first relies on sequencing the transposon insertions sites in an input pool and an output pool of transposon mutants, and using the differential representation of mutants in each pool to infer the functional role of each gene with sufficient representation of insertion sites . The second method relies on creating a mutant pool sufficiently large and complex that it saturates the genome and allows for analysis of regions with statistically fewer, or no, insertions than expected using a non-parametric , Bayesian model , or hidden Markov model (HMM) based analysis ,. Both approaches have been applied to several species of bacteria to investigate genes involved in colonization of hosts –, resistance to antibiotics , characterizing metabolic pathways ,, deducing core essential genomes –, and recently, examining genes involved in colonizing soil environments .
The Rhizobiaceae is a family of alpha-proteobacteria containing three agriculturally important genera of soil bacteria: Rhizobium, Sinorhizobium and Agrobacterium. Members in these genera share a unique relationship with plant hosts. Rhizobium, and Sinorhizobium are both able to enter into an endosymbiotic mutualism with certain species of leguminous plants, in which the Rhizobia fix atmospheric nitrogen into a biologically available form for the plant in return for fixed carbon and energy . This symbiosis is particularly important in the context of agriculturally produced pulse crops, where the Rhizobium legume symbiosis affords farmers the ability to reduce the rate of synthetic nitrogen fertilizers application . Conversely, the relationship of Agrobacterium with its plant host is parasitic. In this symbiosis, Agrobacterium infects the tissues of a plant host and transforms specific virulence genes into the host’s DNA, resulting in tumorgenic growth with altered cellular metabolism that the bacteria then colonize . The formation of several galls at the stem root interface results in a plant infection known as crown gall, that can have a significant impact on the crop yield of stone fruits, berries, and nuts .
Genetic research in Rhizobium, Sinorhizobium, and Agrobacterium has relied heavily on the use of transposon mutagenesis screens. Perhaps the most commonly used transposon in the Rhizobiaceae is the Tn5 transposon –. The use of Tn5 genetic screens is numerous and has helped to elucidate genes involved in metabolism –, desiccation tolerance ,, and cell envelope physiology  for example. Implementation of transposon mutagenesis with the high-throughput techniques of INSeq promise to accelerate the rate at which genetic research in the Rhizobiaceae is currently performed. Furthermore, it would allow for comprehensive genome screens for genes involved in host interactions, metabolism, survival, and possibly plasmid maintenance, under any testable condition.
The mariner class of transposon is a host independent transposon that unlike the random insertion transposons such as Tn5 is known to specifically insert into an organism's genome at thymine-adenine (‘TA’) motifs . Because of this defined insertion preference, transposition events can be modeled in silico in any sequenced genome to understand the defined number of insertion locations that exist. This type of analysis can be further refined to examine insertions per gene or within any defined region of interest in the genome. Furthermore, using a transposon with a defined number of insertion loci allows for robust statistical analysis when used in a transposon insertion sequencing methodology. One such analysis package uses a HMM to predict the essentiality of every gene under a particular growth condition in an organism’s genome given a sufficiently dense INSeq data set . The advantages of this type of analysis is that it does not solely rely on the comparison of an input and output pool of mutants as it uses statistical inference on a defined number of insertion sites to infer an over represented or under represented number of insertions within a given region.
Here we describe the modification of a previous INSeq mutagenesis vector pSAM_Bt, which uses a mariner class transposon with modified IR elements, to selectively capture 15–16 bp of genomic DNA adjacent to the transposon insertion . The DNA fragment is then processed for next-generation DNA sequencing. The new INSeq vector uses a Rhizobiaceae specific promoter to drive the expression of the plasmid-borne transposase, and has had the erythromycin resistance gene cassette replaced with a neomycin resistance gene cassette (ntpII) for ease of selection in the Rhizobiaceae. The host range of transposition with pSAM_Rl was examined and an INSeq experiment was performed on the model organism Rhizobium leguminosarum bv. viciae 3841 (RLV3841) to examine the insertion density that the pSAM_Rl transposon could achieve. Furthermore, the functionality of the transposon in RLV3841 for INSeq was also examined.
Construction of pSAM_Rl and transposition frequency within the Rhizobiaceae
Analysis of pSAM_Rl transposition in the RLV3841 genome
Summary of pSAM_Rl transposition in RLV3841
Mean read countB
Median read countC
Distribution of phenotypic classes throughout the RLV3841 genome
Each RLV3841 replicon was analyzed using the Tn-HMM python module  to classify genes based on the observed density of Tn insertions in each gene within the mixed mutant cell population. For example, a gene which had no detectable Tn insertion sites following DNA sequencing from the pooled mutant DNA suggests cells carrying a mutation in this gene were not maintained in the population. These genes are described as conferring an essential phenotype. The procedure is described in greater detail in the Materials and Methods.
Summary of phenotypic classes across the RLV3841 genome
Genes per specific phenotypic classA
Specific essential genes within the RLV3841 genome
Selected genes in RLV3841 cultured on TY agar with Tn insertions substantially below expected levels
Mean read countC
Selected chromosomal house keeping genes
Putative ATP synthase I
F0F1 ATP synthase subunit A
F0F1 ATP synthase subunit C
F0F1 ATP synthase subunit B'
F0F1 ATP synthase subunit B
F0F1 ATP synthase subunit Epsilon
F0F1 ATP synthase subunit Beta
F0F1 ATP synthase subunit Gamma
F0F1 ATP synthase subunit Alpha
F0F1 ATP synthase subunit Delta
Primosome assembly protein PriA
Cell division protein FtsZ
Putative cell division protein FtsA
Putative cell division protein FtsQ
Putative cell division protein FtsW
Putative cell division protein FtsH
DNA topoisomerase IV subunit B
DNA gyrase subunit A
DNA polymerase III subunit Alpha
Putative DNA polymerase III subunit Alpha
DNA polymerase III subunit Beta
Putative single-stranded-DNA-specific exonuclease
DNA-directed RNA polymerase subunit Beta
Putative DNA-directed RNA polymerase subunit Beta'
Putative DNA-directed RNA polymerase subunit Alpha
Putative ATP-dependant helicase
Putative ATP-dependant helicase
Putative replicative DNA helicase
Selected plasmid borne genes
Putative replication protein A
Putative replication protein B
Putative replication protein C
Putative replication protein A
Putative replication protein B
Putative replication protein C
Putative RepA replication protein
Putative RepB replication protein
Putative RepC replication protein
Putative replication partitioning protein
Putative replication partitioning protein
Putative replication initiation protein RepC
Putative replication protein RepA
Putative replication protein RepB
Putative replication initiation protein RepC
Putative replication protein
Putative replication protein B
Putative replication initiation protein RepC
Previously experimentally confirmed genes with a TY- defective phenotype
Putative carboxy-terminal processing protease precursor
Conserved hypothetical membrane protein
3-Oxoacyl acyl carrier protein synthase
Putative Na+/H+ antiporter subunit D
Examples of genes with higher than expected Tn insertions E
Putative lipid A oxidase
Putative transmembrane component of ABC transporter
Putative transmembrane protein
High throughput forward genetic screening is rapidly being adopted in a diverse range of organisms. INSeq, and similar high-throughput techniques, have been used in several bacterial species, including Salmonella, Mycobacterium, Haemophilus, Vibrio, Pseudomonas, Chronobacteria and Bacillus. The technique is very attractive as it allows for the high-throughput functional screening of an organism's genome under varied growth conditions. Here we present the design and validation of a transposon mutagenesis system that will allow for the application of high-throughput INSeq genetic screening for use within the Rhizobiaceae.
Transposon mutagenesis with pSAM_Rl in RLV3841
We have successfully demonstrated the implementation of a mariner class transposon to mutagenize selected species within the family Rhizobiaceae. We found that the MmeI-adapted mariner transposon harbored on pSAM_Rl could mutagenize R. leguminosarum, A. tumefaciens, and S. meliloti at a high frequency. In the case of RLV3841, it was observed that this high frequency of mutagenesis was dense enough to generate saturating libraries of transposon insertion mutants. We were able to generate transposon mutant libraries that saturated the RLV3841 genome with an insertion density of 0.88, which is higher then the insertion density of the data set used for validation of the HMM used for analysis . This suggests that the combination of mariner based transposon insertion sequencing with the Bayesian based HMM analysis would yield accurate and full genome level results. Furthermore, our analysis used approximately 2.3 M reads. And increased read depth could help increase confidence in the analysis, particularly in resolving the phenotypic classifications of essential and growth defective.
We also observed no bias of mariner transposon insertion across the RLV3841 chromosome and plasmids, suggesting that the transposon inserts randomly and therefore should allow for reliable whole-genome screening approaches. Analysis of read depth showed that there was also little bias in mean read-count per insertion, except in the case of pRL7, which was slightly higher than the other replicons. We suggest that this is a result of pRL7 being maintained at a higher copy number than the other megaplasmids and chromosome, and so sampling of the mutant pool DNA would results in a higher sampling of pRL7 transposon insertion tags than the other plasmids, or chromosome. When the increase in mean read counts is taken into consideration with the insertion density, it appears that although the higher copy number of pRL7 resulted in a higher mean read depth, it did not increase the insertions density, suggesting that the saturation of pRL7 had reached a plateau and an increased presence of pRL7 did not result in a corresponding increase in the number of unique insertions sites.
Analysis of RLV3841 INSeq using a hidden Markov model
Analysis of the TY INSeq data set with the Tn-HMM analysis package assigned accurate phenotypic classification to several genes thought to be essential housekeeping genes necessary for growth under normal conditions. The analysis showed that insertions in genes required for ATP synthesis were absent in the mutant pools. These genes are expected to be essential due to their central role in metabolism and their designation as essential in this assay supports the validity of using the MmeI-adapted mariner in conjunction with the HMM analysis as our INSeq methodology. Furthermore, visual investigation of the transposon insertion density around the region encoding the ATP synthase genes revealed a high insertion density leading up to and after the genes, further supporting the conclusion that the technique is robust and can discern regions of essentiality from those of other states at a high resolution.
In a few instances, a gene expected to be essential was observed to be neutral. RL4697 is annotated as a putative DNA polymerase III alpha sub-unit, and is therefore predicted to be required for proper DNA polymerase function; however, the gene was classified in the neutral category with 96.1% of all potential insertion sites in the locus observed to have insertions. We suggest that RL4697 may be misanotated, as RL1723, another DNA polymerase III alpha subunit, was observed to be essential (Table 3). This highlights another potential use of INSeq in the Rhizobiaceae for validation and quality improvement of genome annotations.
Five of the megaplasmids in the RLV3841 genome were observed to have a set of 3 plasmid replication genes that the INSeq approach identified as essential for plasmid replication and maintenance. The exception was pRL7 which has two sets of replication genes  and therefore functional redundancy may have complicated the the classification of the pRL7 rep genes into the phenotypic classes. The classification of the rep genes on each plasmid as essential provides validation of the INSeq method. Tn insertion within a replication locus would result in the loss of the plasmid from the mutant cell populations harvested for DNA extraction and INSeq DNA sequencing. This result highlights the value of the INSeq approach in the genetic characterization of novel plasmids, as the method is able to identify plasmid encoded genes that are required for plasmid replication and maintenance. Furthermore the method can identify plasmid-encoded genes that provide a fitness advantage to the host under specific growth conditions based on an observed low Tn insertion density.
The genes from previously described mutants with growth defects on TY medium were also examined. Previous work identified  several genes that are important for growth on TY medium. When we compared those results with the results produced by the INSeq analysis of RLV3841 grown on TY medium we observed good concordance with these previously published results. In our results, the four TY related genes were all observed to result in a growth-defective phenotype after ~16 generations of growth, which is in agreement that the interruption of these genes via mutagenesis will result in impaired growth on TY.
The construction and validation of the mariner pSAM_Rl transposon delivery vector as a transposon insertion sequencing tool for use in the Rhizobiziaceae will provide an opportunity for researchers in the Rhizobiaceae community to use a new high throughput genetic screening approach. There are many research opportunities within the Rhizobiaceae that could be examined using the INSeq methodology. For example, use of a INSeq approach in rhizobia to study genes required for rhizosphere colonization and plant infection will help to understand the competition problem observed in inoculant strains, by not only identifying new essential rhizosphere colonization genes, but also identifying mutations that provide a phenotypic growth advantage . Furthermore, the use of INSeq could be used to fully elucidate catabolic pathways, if saturating mutant pools were grown on minimal medium given a single carbon source, and contrasted with the results of similar experiments on rich media. In the near term, we will use the INSeq approach to increase our understanding of the gene networks involved in swarming physiology in RLV3841.
Materials and Methods
Bacterial strains, growth conditions and plasmids
Summary of bacterial strains and plasmids
R. leguminosarum bv. viciae 3841
S. meliloti RM1021
SU47 str-21 StrR
Plasmid-free derivative of A. tumefaciens strain C5 RifR
E. coli SM10 λpir
thi-1 thr leu tonA lacY supE recA::RP4-2-Tc::Mu KanR λpir
E. coli PIR1
F- ∆lac169 rpoS(Am) robA1 creC510 hsdR514 endA- recA1 uidA(∆MluI)::pir-116
AmpR EryR; RP4-oriT, oriR6K, mariner himar1C9 transposase with Bacteroides thetaiotamicron rpoD promoter, MmeI-adapted marinerIR elements
AmpR KanR; RP4-oriT, oriR6K
AmpR ; cloning vector
AmpR KanR; pGEM-T vector containing PCR amplified nptII gene from pSC189
AmpR; pGEM-T vector containing PCR amplified R. leguminosarum 3841 rpoD promoter region
AmpR KanR; EryR in pSAM_Bt replaced with nptII from pSC189
AmpR KanR; pSAM_Km with B. thetatiotamicron rpoD promoter replaced with R. leguminosarum 3841 rpoD promoter region
Construction of pSAM_Rl
Plasmid DNA was isolated using GenElute™ Plasmid Miniprep Kit (Sigma-Aldrich). The Neo/KanR cassette within pSC189  was PCR amplified using primers Tn189KmR_Fwd_XhoI and Tn189_Rev_XbaI (Additional file 1: Table S1) such that the XhoI and XbaI restriction enzyme sites were introduced on the 5’ and 3’ end, respectively. The 991 bp Neo/KanR PCR product was subsequently cloned into the pGEM®-T Easy Vector System creating plasmid pGEM::189KmR. The Neo/KanR cassette in pGEM::189KmR was digested with XhoI and XbaI restriction enzymes creating a 979 bp fragment that was directionally cloned into the pSAM_Bt mariner transposon. The resulting plasmid pSAM_Km was maintained in the E. coli strain PIR1(Invitrogen), which allows for high copy number maintenance of R6Kƴ-oriR plasmids.
Cloning of the RLV3841 rpoD promoter region was carried out by PCR amplifying a 366 bp region upstream of the rpoD gene start codon using primers Rlv_rpoD_Pro_Fwd and Rlv_rpoD_Pro_Rev (Additional file 1: Table S1). The rpoD promoter PCR product had a 5’ BamHI and 3’ NdeI restriction enzyme site introduced, and was subsequently cloned into the pGEM®-T Easy Vector System to create the plasmid pGEM::rpoD. The RLV3841 rpoD promoter region was then excised from pGEM::rpoD using NdeI and BamHI, and directionally cloned into pSAM_Km to create the vector pSAM_Rl (Figure 1). The new pSAM_Rl suicide vector carried an MmeI-adapted mariner transposon harbouring a NeoR/KanR cassette, and had the himar1C9 transposase transcriptionally fused to a RLV3841 rpoD promoter. For use in transposon mutagenesis, the pSAM_Rl construct was electroporated into E. coli SM10λpir (obtained from Dr. Peter Howard, University of Saskatchewan).
Testing pSAM_Rl transposition frequency
Transposition mutagenesis using pSAM_Rl was done in triplicate. Donor and recipient cells were grown in broth culture to an OD600 of approximately 0.8 and were pooled in a ratio of 1000 μL recipient to 500 μL of donor in a 1.5 mL microcentrifuge tube. The conjugation mixture was pelleted at 12,500 rpm for 3 min, washed once with 1000 μL 1X PBS, and resuspended in approximately 100 μL 1X PBS. The cell suspensions were then spotted onto pre-warmed TY agar plates and incubated at 30°C overnight. Conjugation spots were scraped and resuspended in 1000 μL of 1X PBS. Enumeration of transposon mutants was done using TY agar supplemented with 50 μg/mL Neo and the appropriate Rhizobiaceae counter selectable antibiotic (Table 4). Enumeration of total Rhizobiaceae was done on TY agar with the species specific selectable antibiotic.
Generating transposon mutant libraries for sequencing
Six independent conjugations of pSAM_Rl into RLV3841 were conducted on TY agar as described above. After 24 hours incubation at 30°C each of 6 mating spot was scraped and resuspended in 1 mL of 1X PBS and then pooled together in a final volume of 6 mL. For selection on TY agar, 1000 μL of resuspended cells were plated across 2 separate 245 × 245 mm2 (Corning) Neo + Str TY agar plates, in triplicate and incubated for ~48 hour at 30°C.
Following incubation, a faint film like growth was scraped off each plate and resuspended in 1 mL of 1×PBS, vortexed for 1 minutes, and then pelleted at 15, 000 RPM for 10 minutes. The supernatant of each resuspension was very viscous and still contained cells, it was equally aliquoted into 2 × 1.5 mL microcentrifuge tubes. The original pellet, and two tubes of supernatant, were then brought up to a final volume of 1000 μL with 1 M NaCl, vortexed thoroughly, and incubated on ice for 1 h to disrupt the viscous exopolysaccharide diffuse capsule to better collect the cells. The NaCl cell suspensions were then pelleted at 15, 000 RPM for 10 minutes, and the pellets from each replicate were pooled independently and resuspended in 1000 μL of TE buffer (pH 8.0). The resulting 3 mutant pools were used for independent DNA isolation and downstream library preparation.
Preparing sequencing libraries and DNA sequencing
Transposon insertion tags (Tn-tags) consisted of 53 bp of pSAM_Rl transposon sequence, including the 27 bp inverse repeat sequence, and 15–16 bp of adjacent genomic DNA. Library preparation was carried out independently for each of the 3 collected Tn-mutant pools. Tn-tags were prepared for DNA sequencing using a modified version of the INSeq method  to make the sequencing process amendable to the Ion Torrent PGM sequencing platform. Linear PCR products were amplified using the primer Ion Torrent BioSAM (Additional file 1: Table S1), with an annealing temperature of 58.6°C and 500 ng of template DNA. Linear PCR products were purified using a QIAquick PCR Purification Kit (Qiagen) according to the manufacturer's recommended protocol. The biotinylated linear PCR products were then bound to Pierce Streptavidin Magnetic Beads (Thermo Scientific) and enzymatic steps during library preparation were performed as described  with the substitution of Klenow (New England Biolabs), Random Primer 6 (New England Biolabs) and T4 DNA Ligase (New England Biolabs). Additionally, a custom library adapter, INSEQ_Adpt, was used in the adapter ligation step. The final PCR amplification of sequencing template was accomplished using fusion primers designed in accordance with Ion Amplicon Library Preparation (Fusion Method, Pub#: 4468326 Rev. C), using the PCR amplification conditions described in the INSeq methodology . The forward fusion primers IT_A_FP_1, IT_A_FP_2, and IT_A_FP_3 included IonXpress barcode sequences 1, 2, and 3 respectively, for downstream sequence separation. The reverse primer IT_trP1_FP was used in conjunction with a forward primer to introduce the trP1 sequencing adapter. The final sequencing template prepared from the Ion Amplicon library preparation was 187 bp in length. Sequencing template was gel purified using the Invitrogen E-Gel® SizeSelect™ system, and was analyzed using a Bioanalyzer High Sensitivity DNA Chip (Agilent Technologies) prior to sequencing for quality and molarity. The three technical replicates had a final concentration of 1.37, 1.82, and 1.23 μg/ul of sequencing library after size selection, respectively.
DNA sequencing was performed on the Ion Torrent PGM using 200 bp sequencing chemistry and a 316v2 sequencing chip. The total raw sequencing output of the Ion Torrent was 1.25, 1.73, and 1.74 million reads for each of the 3 replicates. The raw sequencing reads were then pooled for downstream data extraction an analysis.
Data extraction and transposon insertion analysis
The authors gratefully acknowledge the support for this work from Natural Sciences and Engineering Research Council of Canada (NSERC) Discovery Grant to CKY. We would also like to thank Dr. Dinah Tambalo for her helpful discussions throughout the project.
- Van Opijnen T, Camilli A: Transposon insertion sequencing: a new tool for systems-level analysis of microorganisms. Nat Rev Microbiol. 2013, 11: 1-8.View ArticleGoogle Scholar
- Barquist L, Boinett CJ, Cain AK: Approaches to querying bacterial genomes with transposon-insertion sequencing. RNA Biol. 2013, 10: 1161-1169. 10.4161/rna.24765.PubMed CentralView ArticlePubMedGoogle Scholar
- Opijnen T, Bodi KL, Camilli A: Tn-seq: high-throughput parallel sequencing for fitness and genetic interaction studies in microorganisms. Nat Methods. 2009, 6: 767-775. 10.1038/nmeth.1377.PubMed CentralView ArticlePubMedGoogle Scholar
- Goodman AL, McNulty NP, Zhao Y, Leip D, Mitra RD, Lozupone CA, Knight R, Gordon JI: Identifying Genetic Determinants Needed to Establish a Human Gut Symbiont in Its Habitat. Cell Host Microbe. 2009, 6: 279-289. 10.1016/j.chom.2009.08.003.PubMed CentralView ArticlePubMedGoogle Scholar
- Langridge GC, Phan M-D, Turner DJ, Perkins TT, Parts L, Haase J, Charles I, Maskell DJ, Peters SE, Dougan G, Wain J, Parkhill J, Turner AK: Simultaneous assay of every Salmonella Typhi gene using one million transposon mutants. Genome Res. 2009, 19: 2308-2316. 10.1101/gr.097097.109.PubMed CentralView ArticlePubMedGoogle Scholar
- Gawronski JD, Wong SMS, Giannoukos G, Ward DV, Akerley BJ: Tracking insertion mutants within libraries by deep sequencing and a genome-wide screen for Haemophilus genes required in the lung. Proc Natl Acad Sci. 2009, 106: 16422-16427. 10.1073/pnas.0906627106.PubMed CentralView ArticlePubMedGoogle Scholar
- Bishop AH, Rachwal PA, Vaid A: Identification of genes required by Bacillus thuringiensis for survival in soil by transposon-directed insertion site sequencing. Curr Microbiol. 2014, 68: 477-485. 10.1007/s00284-013-0502-7.View ArticlePubMedGoogle Scholar
- Zhang YJ, Ioerger TR, Huttenhower C, Long JE, Sassetti CM, Sacchettini JC, Rubin EJ: Global assessment of genomic regions required for growth in Mycobacterium tuberculosis. PLoS Pathog. 2012, 8: e1002946-10.1371/journal.ppat.1002946.PubMed CentralView ArticlePubMedGoogle Scholar
- DeJesus MA, Zhang YJ, Sassetti CM, Rubin EJ, Sacchettini JC, Ioerger TR: Bayesian analysis of gene essentiality based on sequencing of transposon insertion libraries. Bioinformatics. 2013, 29: 695-703. 10.1093/bioinformatics/btt043.PubMed CentralView ArticlePubMedGoogle Scholar
- Dejesus MA, Ioerger TR: A hidden Markov model for identifying essential and growth-defect regions in bacterial genomes from transposon insertion sequencing data. BMC Bioinformatics. 2013, 14: 1-12. 10.1186/1471-2105-14-303.View ArticleGoogle Scholar
- Chao MC, Pritchard JR, Zhang YJ, Rubin EJ, Livny J, Davis BM, Waldor MK: High-resolution definition of the Vibrio cholerae essential gene set with hidden Markov model-based analyses of transposon-insertion sequencing data. Nucleic Acids Res. 2013, 41: 9033-9048. 10.1093/nar/gkt654.PubMed CentralView ArticlePubMedGoogle Scholar
- Skurnik D, Roux D, Aschard H, Cattoir V, Yoder-Himes D, Lory S, Pier GB: A comprehensive analysis of in vitro and in vivo genetic fitness of Pseudomonas aeruginosa using high-throughput sequencing of transposon libraries. PLoS Pathog. 2013, 9: e1003582-10.1371/journal.ppat.1003582.PubMed CentralView ArticlePubMedGoogle Scholar
- Kamp HD, Patimalla-Dipali B, Lazinski DW, Wallace-Gadsden F, Camilli A: Gene fitness landscapes of Vibrio cholerae at important stages of its life cycle. PLoS Pathog. 2013, 9: e1003800-10.1371/journal.ppat.1003800.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang N, Ozer EA, Mandel MJ, Hauser R: Genome-wide identification of Acinetobacter baumannii genes necessary for persistence in the lung. MBio. 2014, 5: e01163-PubMed CentralPubMedGoogle Scholar
- Gallagher LA, Shendure J, Manoil C: Genome-scale identification of resistance functions in Pseudomonas aeruginosa using Tn-seq. MBio. 2011, 2: e00315-10.1128/mBio.00315-10.PubMed CentralView ArticlePubMedGoogle Scholar
- Brutinel ED, Gralnick JA: Anomalies of the anaerobic tricarboxylic acid cycle in Shewanella oneidensis revealed by Tn-seq. Mol Microbiol. 2012, 86: 273-83. 10.1111/j.1365-2958.2012.08196.x.View ArticlePubMedGoogle Scholar
- Griffin JE, Gawronski JD, Dejesus MA, Ioerger TR, Akerley BJ, Sassetti CM: High-resolution phenotypic profiling defines genes essential for mycobacterial growth and cholesterol catabolism. PLoS Pathog. 2011, 7: e1002251-10.1371/journal.ppat.1002251.PubMed CentralView ArticlePubMedGoogle Scholar
- Canals R, Xia X-Q, Fronick C, Clifton SW, Ahmer BMM, Andrews-Polymenis HL, Porwollik S, McClelland M: High-throughput comparison of gene fitness among related bacteria. BMC Genomics. 2012, 13: 212-10.1186/1471-2164-13-212.PubMed CentralView ArticlePubMedGoogle Scholar
- Klein BA, Tenorio EL, Lazinski DW, Camilli A, Duncan MJ, Hu LT: Identification of essential genes of the periodontal pathogen Porphyromonas gingivalis. BMC Genomics. 2012, 13: 578-10.1186/1471-2164-13-578.PubMed CentralView ArticlePubMedGoogle Scholar
- Phan M-D, Peters KM, Sarkar S, Lukowski SW, Allsopp LP, Moriel DG, Achard MES, Totsika M, Marshall VM, Upton M, Beatson SA, Schembri MA: The serum resistome of a globally disseminated multidrug resistant uropathogenic Escherichia coli clone. PLoS Genet. 2013, 9: e1003834-10.1371/journal.pgen.1003834.PubMed CentralView ArticlePubMedGoogle Scholar
- Baugh L, Gallagher LA, Patrapuvich R, Clifton MC, Gardberg AS, Edwards TE, Armour B, Begley DW, Dieterich SH, Dranow DM, Abendroth J, Fairman JW, Fox D, Staker BL, Phan I, Gillespie A, Choi R, Nakazawa-Hewitt S, Nguyen MT, Napuli A, Barrett L, Buchko GW, Stacy R, Myler PJ, Stewart LJ, Manoil C, Voorhis W, Van Voorhis WC: Combining functional and structural genomics to sample the essential Burkholderia structome. PLoS One. 2013, 8: e53851-10.1371/journal.pone.0053851.PubMed CentralView ArticlePubMedGoogle Scholar
- Fels SR, Zane GM, Blake SM, Wall JD: Rapid TnLE-seq for gene fitness evaluation in underdeveloped bacterial systems. Appl Environ Microbiol. 2013, 79: 7501-10.1128/AEM.02051-13.View ArticleGoogle Scholar
- Khatiwara A, Jiang T, Sung S-S, Dawoud T, Kim JN, Bhattacharya D, Kim H-B, Ricke SC, Kwon YM: Genome scanning for conditionally essential genes in Salmonella enterica serotype Typhimurium. Appl Environ Microbiol. 2012, 78: 3098-107. 10.1128/AEM.06865-11.PubMed CentralView ArticlePubMedGoogle Scholar
- Christen B, Abeliuk E, Collier JM, Kalogeraki VS, Passarelli B, Coller JA, Fero MJ, McAdams HH, Shapiro L: The essential genome of a bacterium. Mol Syst Biol. 2011, 7: 528-10.1038/msb.2011.58.PubMed CentralView ArticlePubMedGoogle Scholar
- Yanagi M, Yamasato K: Phylogenetic analysis of the family Rhizobiaceae and related bacteria by sequencing of 16S rRNA gene using PCR and DNA sequencer. FEMS Microbiol Lett. 1993, 107: 115-120. 10.1111/j.1574-6968.1993.tb06014.x.View ArticlePubMedGoogle Scholar
- Wielbo J: Rhizobial communities in symbiosis with legumes: genetic diversity, competition and interactions with host plants. Cent Eur J Biol. 2012, 7: 363-372. 10.2478/s11535-012-0032-5.Google Scholar
- Lupwayi NZ, Kennedy AC: Grain legumes in northern great plains. Agron J. 2007, 99: 1700-10.2134/agronj2006.0313s.View ArticleGoogle Scholar
- Pitzschke A: Agrobacterium infection and plant defense-transformation success hangs by a thread. Front Plant Sci. 2013, 4 (December): 519-PubMed CentralPubMedGoogle Scholar
- Pulawska J: Crown gall of stone fruits and nuts, economic significance and diversity of its causal agents: tumorigenic Agrobacterium spp. J Plant Pathol. 2010, 92: 87-98.Google Scholar
- Beringer JE, Beynon JL, Buchanan-Wollaston AV, Johnston AWB: Transfer of the drug-resistance transposon Tn5 to Rhizobium. Nature. 1978, 276: 633-634. 10.1038/276633a0.View ArticleGoogle Scholar
- Wilson KJ, Sessitsch A, Corbo JC, Giller KE, Akkermans AD, Jefferson RA: beta-Glucuronidase (GUS) transposons for ecological and genetic studies of rhizobia and other gram-negative bacteria. Microbiology. 1995, 141 (Pt 5): 1691-705. 10.1099/13500872-141-7-1691.View ArticlePubMedGoogle Scholar
- Reeve WG, Tiwari RP, Worsley PS, Dilworth MJ, Glenn AR, Howieson JG: Constructs for insertional mutagenesis, transcriptional signal localization and gene regulation studies in root nodule and other bacteria. Microbiology. 1999, 145 (Pt 6): 1307-1316. 10.1099/13500872-145-6-1307.View ArticlePubMedGoogle Scholar
- Ding H, Yip CB, Geddes BA, Oresnik IJ, Hynes MF: Glycerol utilization by Rhizobium leguminosarum requires an ABC transporter and affects competition for nodulation. Microbiology. 2012, 158 (Pt 5): 1369-78. 10.1099/mic.0.057281-0.View ArticlePubMedGoogle Scholar
- Oresnik IJ, Pacarynuk LA, O’Brien SAP, Yost CK, Hynes MF: Plasmid-encoded catabolic genes in Rhizobium leguminosarum bv. trifolii: evidence for a plant-inducible rhamnose locus Iinvolved in competition for nodulation. Mol Plant Microbe Interact. 1998, 11: 1175-1185. 10.1094/MPMI.1918.104.22.1685.View ArticleGoogle Scholar
- Geddes BA, Oresnik IJ: Inability to catabolize galactose leads to increased ability to compete for nodule occupancy in Sinorhizobium meliloti. J Bacteriol. 2012, 194: 5044-53. 10.1128/JB.00982-12.PubMed CentralView ArticlePubMedGoogle Scholar
- Vanderlinde EM, Harrison JJ, Muszyński A, Carlson RW, Turner RJ, Yost CK: Identification of a novel ABC transporter required for desiccation tolerance, and biofilm formation in Rhizobium leguminosarum bv. viciae 3841. FEMS Microbiol Ecol. 2010, 71: 327-40. 10.1111/j.1574-6941.2009.00824.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Vanderlinde EM, Muszynski A, Harrison JJ, Koval SF, Foreman DL, Ceri H, Kannenberg EL, Carlson RW, Yost CK: Rhizobium leguminosarum biovar viciae 3841, deficient in 27-hydroxyoctacosanoate-modified lipopolysaccharide, is impaired in desiccation tolerance, biofilm formation and motility. Microbiology. 2009, 155 (Pt 9): 3055-69. 10.1099/mic.0.025031-0.PubMed CentralView ArticlePubMedGoogle Scholar
- Vanderlinde EM, Magnus SA, Tambalo DD, Koval SF, Yost CK: Mutation of a broadly conserved operon (RL3499-RL3502) from Rhizobium leguminosarum biovar viciae causes defects in cell morphology and envelope integrity. J Bacteriol. 2011, 193: 2684-94. 10.1128/JB.01456-10.PubMed CentralView ArticlePubMedGoogle Scholar
- Plasterk RHA, Izsvák Z, Ivics Z: Resident aliens. Trends Genet. 1999, 15: 326-332. 10.1016/S0168-9525(99)01777-1.View ArticlePubMedGoogle Scholar
- Dong TG, Ho BT, Yoder-Himes DR: Identification of T6SS-dependent effector and immunity proteins by Tn-seq in Vibrio cholerae. Proc Natl Acad Sci. 2013, 110: 2623-2628. 10.1073/pnas.1222783110.PubMed CentralView ArticlePubMedGoogle Scholar
- Ding H, Hynes MF: Plasmid transfer systems in the rhizobia. Can J Microbiol. 2009, 55: 917-927. 10.1139/W09-056.View ArticlePubMedGoogle Scholar
- Beringer JE: R factor transfer in Rhizobium leguminosarum. J Gen Microbiol. 1974, 84: 188-98. 10.1099/00221287-84-1-188.PubMedGoogle Scholar
- Sambrook J, Fritsch EF, Maniatis T: Molecular Cloning: A Laboratory Manual. 1989, Cold Spring Harbour, New York, 2Google Scholar
- Poole PS, Blyth A, Reid CJ, Walters K: myo-Inositol catabolism and catabolite regulation in Rhizobium leguminosarum bv. viciae. Microbiology. 1994, 140: 2787-2795. 10.1099/00221287-140-10-2787.View ArticleGoogle Scholar
- Meade HM, Long SR, Ruvkun GB, Brown SE, Ausubel FM: Physical and genetic characterization of symbiotic and auxotrophic mutants of Rhizobium meliloti induced by transposon Tn5 mutagenesis. J Bacteriol. 1982, 149: 114-22.PubMed CentralPubMedGoogle Scholar
- Hynes MF, Simon R, Pühler A: The development of plasmid-free strains of Agrobacterium tumefaciens by using incompatibility with a Rhizobium meliloti plasmid to eliminate pAtC58. Plasmid. 1985, 13: 99-105. 10.1016/0147-619X(85)90062-9.View ArticlePubMedGoogle Scholar
- Miller VL, Mekalanos JJ: A novel suicide vector and its use in construction of insertion mutations: osmoregulation of outer membrane proteins and virulence determinants in Vibrio cholerae requires toxR. J Bacteriol. 1988, 170: 2575-2583.PubMed CentralPubMedGoogle Scholar
- Chiang SL, Rubin EJ: Construction of a mariner-based transposon for epitope-tagging and genomic targeting. Gene. 2002, 296: 179-185. 10.1016/S0378-1119(02)00856-9.View ArticlePubMedGoogle Scholar
- Goodman AL, Wu M, Gordon JI: Identifying microbial fitness determinants by insertion sequencing using genome-wide transposon mutant libraries. Nat Protoc. 2011, 6: 1969-80. 10.1038/nprot.2011.417.PubMed CentralView ArticlePubMedGoogle Scholar
- Martin M: Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet Journal. 2011, 17: 10-10.14806/ej.17.1.200.View ArticleGoogle Scholar
- Young JPW, Crossman LC, Johnston AWB: The genome of Rhizobium leguminosarum has recognizable core and accessory components. Genome. 2006, 7: R34-Google Scholar
- Langmead B, Trapnell C, Pop M, Salzberg SL: Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10: R25-10.1186/gb-2009-10-3-r25.PubMed CentralView ArticlePubMedGoogle Scholar
- Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R: The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009, 25: 2078-9. 10.1093/bioinformatics/btp352.PubMed CentralView ArticlePubMedGoogle Scholar
- Robinson JT, Thorvaldsdóttir H, Winckler W, Guttman M, Lander ES, Getz G, Mesirov JP: Integrative genomics viewer. Nat Biotechnol. 2011, 29: 24-26. 10.1038/nbt.1754.PubMed CentralView ArticlePubMedGoogle Scholar
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.