Construction of Signature-tagged Mutant Library in Mesorhizobium loti as a Powerful Tool for Functional Genomics

Rhizobia are nitrogen-fixing soil bacteria that establish endosymbiosis with some leguminous plants. The completion of several rhizobial genome sequences provides opportunities for genome-wide functional studies of the physiological roles of many rhizobial genes. In order to carry out genome-wide phenotypic screenings, we have constructed a large mutant library of the nitrogen-fixing symbiotic bacterium, Mesorhizobium loti, by transposon mutagenesis. Transposon insertion mutants were generated using the signature-tagged mutagenesis (STM) technique and a total of 29 330 independent mutants were obtained. Along with the collection of transposon mutants, we have determined the transposon insertion sites for 7892 clones, and confirmed insertions in 3680 non-redundant M. loti genes (50.5% of the total number of M. loti genes). Transposon insertions were randomly distributed throughout the M. loti genome without any bias toward G+C contents of insertion target sites and transposon plasmids used for the mutagenesis. We also show the utility of STM mutants by examining the specificity of signature tags and test screenings for growth- and nodulation-deficient mutants. This defined mutant library allows for genome-wide forward- and reverse-genetic functional studies of M. loti and will serve as an invaluable resource for researchers to further our understanding of rhizobial biology.Key words: Mesorhizobium loti, signature-tagged mutagenesis, mutant library, reverse genetics     

ITIS,a bioinformatics tool for accurate identification of transposon insertion sites using next-generation sequencing data

Background

Transposable elements constitute an important part of the genome and are essential in adaptive mechanisms. Transposition events associated with phenotypic changes occur naturally or are induced in insertional mutant populations. Transposon mutagenesis results in multiple random insertions and recovery of most/all the insertions is critical for forward genetics study. Using genome next-generation sequencing data and appropriate bioinformatics tool, it is plausible to accurately identify transposon insertion sites, which could provide candidate causal mutations for desired phenotypes for further functional validation.

Results

We developed a novel bioinformatics tool, ITIS (Identification of Transposon Insertion Sites), for localizing transposon insertion sites within a genome. It takes next-generation genome re-sequencing data (NGS data), transposon sequence, and reference genome sequence as input, and generates a list of highly reliable candidate insertion sites as well as zygosity information of each insertion. Using a simulated dataset and a case study based on an insertional mutant line from Medicago truncatula, we showed that ITIS performed better in terms of sensitivity and specificity than other similar algorithms such as RelocaTE, RetroSeq, TEMP and TIF. With the case study data, we demonstrated the efficiency of ITIS by validating the presence and zygosity of predicted insertion sites of the Tnt1 transposon within a complex plant system, M. truncatula.

Conclusion

This study showed that ITIS is a robust and powerful tool for forward genetic studies in identifying transposable element insertions causing phenotypes. ITIS is suitable in various systems such as cell culture, bacteria, yeast, insect, mammal and plant.

Electronic supplementary material

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

In Vivo Himar1-Based Transposon Mutagenesis of Francisella tularensis

Francisella tularensis is the intracellular pathogen that causes human tularemia. It is recognized as a potential agent of bioterrorism due to its low infectious dose and multiple routes of entry. We report the development of a Himar1-based random mutagenesis system for F. tularensis (HimarFT). In vivo mutagenesis of F. tularensis live vaccine strain (LVS) with HimarFT occurs at high efficiency. Approximately 12 to 15% of cells transformed with the delivery plasmid result in transposon insertion into the genome. Results from Southern blot analysis of 33 random isolates suggest that single insertions occurred, accompanied by the loss of the plasmid vehicle in most cases. Nucleotide sequence analysis of rescued genomic DNA with HimarFT indicates that the orientation of integration was unbiased and that insertions occurred in open reading frames and intergenic and repetitive regions of the chromosome. To determine the utility of the system, transposon mutagenesis was performed, followed by a screen for growth on Chamberlain's chemically defined medium (CDM) to isolate auxotrophic mutants. Several mutants were isolated that grew on complex but not on the CDM. We genetically complemented two of the mutants for growth on CDM with a newly constructed plasmid containing a nourseothricin resistance marker. In addition, uracil or aromatic amino acid supplementation of CDM supported growth of isolates with insertions in pyrD, carA, or aroE1 supporting the functional assignment of genes within each biosynthetic pathway. A mutant containing an insertion in aroE1 demonstrated delayed replication in macrophages and was restored to the parental growth phenotype when provided with the appropriate plasmid in trans. Our results suggest that a comprehensive library of mutants can be generated in F. tularensis LVS, providing an additional genetic tool to identify virulence determinants required for survival within the host.     

Genome-Scale Metabolic Network Validation of Shewanella oneidensis Using Transposon Insertion Frequency Analysis

Transposon mutagenesis, in combination with parallel sequencing, is becoming a powerful tool for en-masse mutant analysis. A probability generating function was used to explain observed miniHimar transposon insertion patterns, and gene essentiality calls were made by transposon insertion frequency analysis (TIFA). TIFA incorporated the observed genome and sequence motif bias of the miniHimar transposon. The gene essentiality calls were compared to: 1) previous genome-wide direct gene-essentiality assignments; and, 2) flux balance analysis (FBA) predictions from an existing genome-scale metabolic model of Shewanella oneidensis MR-1. A three-way comparison between FBA, TIFA, and the direct essentiality calls was made to validate the TIFA approach. The refinement in the interpretation of observed transposon insertions demonstrated that genes without insertions are not necessarily essential, and that genes that contain insertions are not always nonessential. The TIFA calls were in reasonable agreement with direct essentiality calls for S. oneidensis, but agreed more closely with E. coli essentiality calls for orthologs. The TIFA gene essentiality calls were in good agreement with the MR-1 FBA essentiality predictions, and the agreement between TIFA and FBA predictions was substantially better than between the FBA and the direct gene essentiality predictions.     

Regulatory and DNA Repair Genes Contribute to the Desiccation Resistance of Sinorhizobium meliloti Rm1021

Sinorhizobium meliloti can form a nitrogen-fixing symbiotic relationship with alfalfa after bacteria in the soil infect emerging root hairs of the growing plant. To be successful at this, the bacteria must be able to survive in the soil between periods of active plant growth, including when conditions are dry. The ability of S. meliloti to withstand desiccation has been known for years, but genes that contribute to this phenotype have not been identified. Transposon mutagenesis was used in combination with novel screening techniques to identify four desiccation-sensitive mutants of S. meliloti Rm1021. DNA sequencing of the transposon insertion sites identified three genes with regulatory functions (relA, rpoE2, and hpr) and a DNA repair gene (uvrC). Various phenotypes of the mutants were determined, including their behavior on several indicator media and in symbiosis. All of the mutants formed an effective symbiosis with alfalfa. To test the hypothesis that UvrC-related excision repair was important in desiccation resistance, uvrA, uvrB, and uvrC deletion mutants were also constructed. These strains were sensitive to DNA damage induced by UV light and 4-NQO and were also desiccation sensitive. These data indicate that uvr gene-mediated DNA repair and the regulation of stress-induced pathways are important for desiccation resistance.     

mariner-Based Transposon Mutagenesis of Rickettsia prowazekii

Rickettsia prowazekii, the causative agent of epidemic typhus, is an obligate intracellular bacterium that grows directly within the cytoplasm of its host cell, unbounded by a vacuolar membrane. The obligate intracytoplasmic nature of rickettsial growth places severe restrictions on the genetic analysis of this distinctive human pathogen. In order to expand the repertoire of genetic tools available for the study of this pathogen, we have employed the versatile mariner-based, Himar1 transposon system to generate insertional mutants of R. prowazekii. A transposon containing the R. prowazekii arr-2 rifampin resistance gene and a gene coding for a green fluorescent protein (GFP_UV) was constructed and placed on a plasmid expressing the Himar1 transposase. Electroporation of this plasmid into R. prowazekii resulted in numerous transpositions into the rickettsial genome. Transposon insertion sites were identified by rescue cloning, followed by DNA sequencing. Random transpositions integrating at TA sites in both gene coding and intergenic regions were identified. Individual rickettsial clones were isolated by the limiting-dilution technique. Using both fixed and live-cell techniques, R. prowazekii transformants expressing GFP_UV were easily visible by fluorescence microscopy. Thus, a mariner-based system provides an additional mechanism for generating rickettsial mutants that can be screened using GFP_UV fluorescence.     

Establishment of an Arbitrary PCR for Rapid Identification of Tn917 Insertion Sites in Staphylococcus epidermidis: Characterization of Biofilm-Negative and Nonmucoid Mutants

Transposon mutagenesis with the Enterococcus faecalis transposon Tn917 is a genetic approach frequently used to identify genes related with specific phenotypes in gram-positive bacteria. We established an arbitrary PCR for the rapid and easy identification of Tn917 insertion sites in Staphylococcus epidermidis with six independent, well-characterized biofilm-negative Tn917 transposon mutants, which were clustered in the icaADBC gene locus or harbor Tn917 in the regulatory gene rsbU. For all six of these mutants, short chromosomal DNA fragments flanking both transposon ends could be amplified. All fragments were sufficient to correctly identify the Tn917 insertion sites in the published S. epidermidis genomes. By using this technique, the Tn917 insertion sites of three not-yet-characterized biofilm-negative or nonmucoid mutants were identified. In the biofilm-negative and nonmucoid mutant M12, Tn917 is inserted into a gene homologous to the regulatory gene purR of Bacillus subtilis and Staphylococcus aureus. The Tn917 insertions of the nonmucoid but biofilm-positive mutants M16 and M20 are located in genes homologous to components of the phosphoenolpyruvate-sugar phosphotransferase system (PTS) of B. subtilis, S. aureus, and Staphylococcus carnosus, indicating an influence of the PTS on the mucoid phenotype in S. epidermidis.     

Degradation of Substituted Phenylurea Herbicides by Arthrobacter globiformis Strain D47 and Characterization of a Plasmid-Associated Hydrolase Gene, puhA

Arthrobacter globiformis D47 was shown to degrade a range of substituted phenylurea herbicides in soil. This strain contained two plasmids of approximately 47 kb (pHRIM620) and 34 kb (pHRIM621). Plasmid-curing experiments produced plasmid-free strains as well as strains containing either the 47- or the 34-kb plasmid. The strains were tested for their ability to degrade diuron, which demonstrated that the degradative genes were located on the 47-kb plasmid. Studies on the growth of these strains indicated that the ability to degrade diuron did not offer a selective advantage to A. globiformis D47 on minimal medium designed to contain the herbicide as a sole carbon source. The location of the genes on a plasmid and a lack of selection would explain why the degradative phenotype, as with many other pesticide-degrading bacteria, can be lost on subculture. A 22-kb EcoRI fragment of plasmid pHRIM620 was expressed in Escherichia coli and enabled cells to degrade diuron. Transposon mutagenesis of this fragment identified one open reading frame that was essential for enzyme activity. A smaller subclone of this gene (2.5 kb) expressed in E. coli coded for the protein that degraded diuron. This gene and its predicted protein sequence showed only a low level of protein identity (25% over ca. 440 amino acids) to other database sequences and was named after the enzyme it encoded, phenylurea hydrolase (puhA gene).     

Genes for phaseolotoxin synthesis are located on the chromosome ofPseudomonas syringae pv.phaseolicola

Transposon mutagenesis was used to locate the genes for phaseolotoxin biosynthesis inPseudomonas syringae pathovar.phaseolicola. Mutants unable to produce toxin were obtained that carried Tn5 on different chromosomal restriction fragments. None of the Tn5-induced nontoxigenic mutants carried the transposon in plasmid DNA. The insertion of Tn5 intotox DNA was confirmed by site-directed mutagenesis. The results reported here suggest the involvement of at least five chromosomal genes in phaseolotoxin biosynthesis. All of the toxinminus mutants retained both full pathogenicity on beans and resistance to the toxin.     

The Hermes transposon of Musca domestica and its use as a mutagen of Schizosaccharomyces pombe

Transposon mutagenesis allows for the discovery and characterization of genes by creating mutations that can be easily mapped and sequenced. Moreover, this method allows for a relatively unbiased approach to isolating genes of interest. Recently, a system of transposon based mutagenesis for Schizosaccharomyces pombe became available. This mutagenesis relies on Hermes, a DNA transposon from the house fly that readily integrates into the chromosomes of S. pombe. The Hermes system is distinct from the retrotransposons of S. pombe because it efficiently integrates into open reading frames. To mutagenize S. pombe, cells are transformed with a plasmid that contains a drug resistance marker flanked by the terminal inverted repeats of Hermes. The Hermes transposase expressed from a second plasmid excises the resistance marker with the inverted repeats and inserts this DNA into chromosomal sites. After S. pombe with these two plasmids grow 25 generations, approximately 2% of the cells contain insertions. Of the cells with insertions, 68% contain single integration events. The protocols listed here provide the detailed information necessary to mutagenize a strain of interest, screen for specific phenotypes, and sequence the positions of insertion.     

Identification and analysis of the Shewanella oneidensis major oxygen-independent coproporphyrinogen III oxidase gene

Shewanella oneidenesis MR-1 is a facultative anaerobe that can use a large number of electron acceptors including metal oxides. During anaerobic respiration, S. oneidensis MR-1 synthesizes a large number of c cytochromes that give the organism its characteristic orange color. Using a modified mariner transposon, a number of S. oneidensis mutants deficient in anaerobic respiration were generated. One mutant, BG163, exhibited reduced pigmentation and was deficient in c cytochromes normally synthesized under anaerobic condition. The deficiencies in BG163 were due to insertional inactivation of hemN1, which exhibits a high degree of similarity to genes encoding anaerobic coproporphyrinogen III oxidases that are involved in heme biosynthesis. The ability of BG163 to synthesize c cytochromes under anaerobic conditions, and to grow anaerobically with different electron acceptors was restored by the introduction of hemN1 on a plasmid. Complementation of the mutant was also achieved by the addition of hemin to the growth medium. The genome sequence of S. oneidensis contains three putative anaerobic coproporphyrinogen III oxidase genes. The protein encoded by hemN1 appears to be the major enzyme that is involved in anaerobic heme synthesis of S. oneidensis. The other two putative anaerobic coproporphyrinogen III oxidase genes may play a minor role in this process.     

Sequencing and analysis of plasmids pAV1 and pAV2 ofAcinetobacter venetianus VE-C3 involved in diesel fuel degradation

Acinetobacter venetianus strain VE-C3 was isolated in the Venice lagoon (Italy) as a strain able to degrade diesel fuel oil. This strain possesses genes of the alkane monoxygenase complex responsible forn-alkane degradation and carries two plasmids, pAV1 (10820 bp) and pAV2 (15135 bp), which were supposed from the analysis of Alk mutant strains to harbour genetic determinants for hydrocarbon degradation. In this work we determined the nucleotide sequence of both plasmids and showed the presence of a putative aldehyde dehydrogenase gene, essential for hydrocarbon degradation, on plasmid pAV2, and of an ORF similar toalkL gene present on pAV1 plasmid. These data, combined with genetic reports indicating that strains lacking one of the two plasmids or carrying transposon insertion on pAV1, are defective inn-alkane degradation, suggest a complex genomic organisation of genes involved in alkane degradation inA. venetianus VE-C3. In this bacterium these genes are carried by both the chromosome and the plasmids, while inAcinetobacter sp. strain ADP1 and M1 all the genes for alkane monoxygenase complex are located only on the chromosome.     

Rapid Transposon Liquid Enrichment Sequencing (TnLE-seq) for Gene Fitness Evaluation in Underdeveloped Bacterial Systems

Whole-genome fitness analysis in microbes that uses saturating transposon mutagenesis combined with massively parallel sequencing (Tn-seq) is providing a measure of the contribution of each gene to a given growth condition. With this technique, gene fitness profiles and essential genes are discovered by simultaneous analyses of whether the absence of each gene product alters the growth kinetics of the bacterium. Here we modify the standard Tn-seq procedure to simplify and shorten the process by including delivery of the transposon through conjugation and liquid culture enrichment of the mutant pool, creating transposon liquid enrichment sequencing (TnLE-seq). To illustrate the success of these modifications and the robustness of the procedure, analyses of gene fitness of two cultures of the strictly anaerobic bacterium Desulfovibrio vulgaris Hildenborough were performed, with growth on lactate as the electron donor and sulfate as the electron acceptor. These data demonstrate reproducibility and provide a base condition for analysis of fitness changes in deletion mutants and in various growth conditions. The procedural modifications will facilitate the application of this powerful genetic analysis to microbes lacking a facile genetic system. Pilot studies produced 2.5 × 10⁵ and 3.4 × 10⁵ unique insertion mutants in the anaerobe Desulfovibrio vulgaris Hildenborough grown under typical laboratory conditions in rich medium. These analyses provided two similar high-resolution maps of gene fitness across the genome, and the method was also applied to growth in minimal medium. These results were also compared to the coverage obtained with a ca. 13,000-member cataloged transposon library constructed by sequencing transposon insertion sites in individual mutants.     

Efficient transposition of IS204-derived plasmids in Streptomyces coelicolor

In order to study functional gene expression in Streptomyces coelicolor, a mini-transposon encoding the apramycin resistance gene aac(3)IV within its inverted repeat (IR) boundaries was constructed based on IS204, which was previously identified in the genome of Nocardia asteroides YP21. The mini-transposon and IS204 transposase gene were then put on a kanamycin-resistant conjugative plasmid pDZY101 that can only replicate in Escherichia coli. After mating with S. coelicolor A3(2) M145, resistant colonies arose efficiently on both apramycin and kanamycin plates. Plasmid rescue indicated that entire plasmids were inserted into the M145 genome with cleavage at an inverted repeat junction formed by the right inverted repeat (IRR) and the last 18 bp of the transposase gene, while the left inverted repeat (IRL) was untouched. Southern blot analysis of the mutants using an aac(3)IV gene probe showed that transposition of plasmid pDZY101 was genetically stable, with a single-copy insertion within the S. coelicolor M145 genome. Several mutagenesis libraries of S. coelicolor M145 were constructed using plasmid pDZY101 derivatives and the transposon insertion site was determined. The correlation between novel mutant phenotypes and previously uncharacterized genes was established and these transposon locations were widely scattered around the genome.     

A Universal Mariner Transposon System for Forward Genetic Studies in the Genus Clostridium

DNA transposons represent an essential tool in the armoury of the molecular microbiologist. We previously developed a catP-based mini transposon system for Clostridium difficile in which the expression of the transposase gene was dependent on a sigma factor unique to C. difficile, TcdR. Here we have shown that the host range of the transposon is easily extended through the rapid chromosomal insertion of the tcdR gene at the pyrE locus of the intended clostridial target using Allele-Coupled Exchange (ACE). To increase the effectiveness of the system, a novel replicon conditional for plasmid maintenance was developed, which no longer supports the effective retention of the transposon delivery vehicle in the presence of the inducer isopropyl β-D-1-thiogalactopyranoside (IPTG). As a consequence, those thiamphenicol resistant colonies that arise in clostridial recipients, following plating on agar medium supplemented with IPTG, are almost exclusively due to insertion of the mini transposon into the genome. The system has been exemplified in both Clostridium acetobutylicum and Clostridium sporogenes, where transposon insertion has been shown to be entirely random. Moreover, appropriate screening of both libraries resulted in the isolation of auxotrophic mutants as well as cells deficient in spore formation/germination. This strategy is capable of being implemented in any Clostridium species.     

Novel Insight into the Genetic Context of the cadAB Genes from a 4-chloro-2-methylphenoxyacetic Acid-Degrading Sphingomonas

The 2-methyl-4-chlorophenoxyacetic (MCPA) acid-degrader Sphingomonas sp. ERG5 has recently been isolated from MCPA-degrading bacterial communities. Using Illumina-sequencing, the 5.7 Mb genome of this isolate was sequenced in this study, revealing the 138 kbp plasmid pCADAB1 harboring the 32.5 kbp composite transposon Tn6228 which contains genes encoding proteins for the removal of 2,4-dichlorophenoxyacetic acid (2,4-D) and MCPA, as well as the regulation of this pathway. Transposon Tn6228 was confirmed by PCR to be situated on the plasmid and also exist in a circular intermediate state - typical of IS3 elements. The canonical tfdAα-gene of group III 2,4-D degraders, encoding the first step in degradation of 2,4-D and related compounds, was not present in the chromosomal contigs. However, the alternative cadAB genes, also providing the initial degradation step, were found in Tn6228, along with the 2,4-D-degradation-associated genes tfdBCDEFKR and cadR. Putative reductase and ferredoxin genes cadCD of Rieske non-heme iron oxygenases were also present in close proximity to cadAB, suggesting that these might have an unknown role in the initial degradation reaction. Parts of the composite transposon contain sequence displaying high similarity to previously analyzed 2,4-D degradation genes, suggesting rapid dissemination and high conservation of the chlorinated-phenoxyacetic acid (PAA)-degradation genotype among the sphingomonads.     

Isolation of Predation-Deficient Mutants of Bdellovibrio bacteriovorus by Using Transposon Mutagenesis

The availability of the complete genome sequence of Bdellovibrio bacteriovorus provides an opportunity for investigating genes that play a significant role in predation. Using two independently derived facultatively predatory Bdellovibrio strains, we have designed a method to cultivate and screen transposon insertion mutants in 96-well microtiter dishes. Transposon insertion mutants were produced by introducing the plasposon pRL27, which carries a mini-Tn5. Mutants have been screened for predatory activity using 96-well plates. Seventeen independent nonpredatory mutants have been isolated, and DNA flanking the insertion has been sequenced. BLAST analysis revealed that most of these interrupted DNA sequences do not code for known proteins or functions. Two of the inactivated genes were analyzed further: one was found to code for a putative serine protease and the other a probable protein involved in secretion through the outer membrane. The methods described here are the first for the generation and isolation of predation-deficient mutants using random-transposon-insertion mutagenesis. As more mutants are isolated and their gene products analyzed, more light will be shed on how this predator carries out its exclusive life processes and perhaps how these products, or the organism itself, can be used for therapeutic, agricultural, and/or other purposes.     

Soluble fumarate reductase isoenzymes from Saccharomyces cerevisiae are required for anaerobic growth

In Saccharomyces cerevisiae, the cytosolic and promitochondrial isoenzymes of fumarate reductase are encoded by the FRDS and OSM1 genes, respectively. The product of the OSM1 gene is reported to be required for growth in hypertonic medium. Simultaneous disruption of the FRDS and OSM1 genes resulted in the inability of the yeasts to grow anaerobically on glucose as a carbon source, and disruption of the OSM1 gene caused poor growth under anaerobic conditions. However, the disruption of both the FRDS and/or OSM1 genes had no effect on aerobic growth or growth under hypertonic conditions. These results suggest that the fumarate reductase isoenzymes in Saccharomyces cerevisiae are essential for anaerobic growth but not for growth under hypertonic conditions.     

Discovery of photosynthesis genes through whole-genome sequencing of acetate-requiring mutants of Chlamydomonas reinhardtii

Large-scale mutant libraries have been indispensable for genetic studies, and the development of next-generation genome sequencing technologies has greatly advanced efforts to analyze mutants. In this work, we sequenced the genomes of 660 Chlamydomonas reinhardtii acetate-requiring mutants, part of a larger photosynthesis mutant collection previously generated by insertional mutagenesis with a linearized plasmid. We identified 554 insertion events from 509 mutants by mapping the plasmid insertion sites through paired-end sequences, in which one end aligned to the plasmid and the other to a chromosomal location. Nearly all (96%) of the events were associated with deletions, duplications, or more complex rearrangements of genomic DNA at the sites of plasmid insertion, and together with deletions that were unassociated with a plasmid insertion, 1470 genes were identified to be affected. Functional annotations of these genes were enriched in those related to photosynthesis, signaling, and tetrapyrrole synthesis as would be expected from a library enriched for photosynthesis mutants. Systematic manual analysis of the disrupted genes for each mutant generated a list of 253 higher-confidence candidate photosynthesis genes, and we experimentally validated two genes that are essential for photoautotrophic growth, CrLPA3 and CrPSBP4. The inventory of candidate genes includes 53 genes from a phylogenomically defined set of conserved genes in green algae and plants. Altogether, 70 candidate genes encode proteins with previously characterized functions in photosynthesis in Chlamydomonas, land plants, and/or cyanobacteria; 14 genes encode proteins previously shown to have functions unrelated to photosynthesis. Among the remaining 169 uncharacterized genes, 38 genes encode proteins without any functional annotation, signifying that our results connect a function related to photosynthesis to these previously unknown proteins. This mutant library, with genome sequences that reveal the molecular extent of the chromosomal lesions and resulting higher-confidence candidate genes, will aid in advancing gene discovery and protein functional analysis in photosynthesis.