Previous Article | Next Article ![]()
Applied and Environmental Microbiology, September 2006, p. 5895-5907, Vol. 72, No. 9
0099-2240/06/$08.00+0 doi:10.1128/AEM.00217-06
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Warwick HRI, University of Warwick, Warwick CV35 9EF, United Kingdom,1 Department of Biological Sciences, University of Warwick, Coventry CV4 7AL, United Kingdom,2 Division of Medical Microbiology, University of Liverpool, Duncan Building, Daulby Street, Liverpool L69 3GA, United Kingdom3
Received 27 January 2006/ Accepted 19 June 2006
|
|
|---|
|
|
|---|
Multilocus sequence typing (MLST) types bacteria based on the sequence of a number of housekeeping genes and is thought to reflect the phylogeny of strains more accurately than other methods (11). Also, because a number of genes at different locations on the chromosome are sequenced, the amount of recombination occurring in a population can be calculated. The method has been applied almost exclusively to important bacterial pathogens and not to other types of bacteria. Xenorhabdus strains have the potential to become commercially important genera, producing oral insecticidal toxins for use in biocontrol and novel antibiotics with the potential to control clinically important microorganisms. Although the type strains have been well studied, and both an Xenorhabdus nematophila strain and an Xenorhabdus bovienii strain have been sequenced (http://xenorhabdus.danforthcenter.org), the extent of diversity within the species is largely unknown. Therefore, whether populations exhibit enough diversity to warrant screening for novel toxins and/or antibiotics, not present in the type strains, is not yet known.
The xpt genes of X. nematophila PMFI296 are responsible for oral insecticidal activity against a range of Lepidopteran insects (24). The four xpt genes involved in this activity, xptA1, xptA2, xptB1, and xptC1, are tightly linked on the chromosome of X. nematophila PMFI296 (13). Analogous genes have been found in P. luminescens (3), Yersinia pestis (33) and Serratia entomophilia (8). To agree with current nomenclature for the toxin genes, two of the original names have been changed: the xptB1 gene has been renamed xptC1, and xptC1 has been named xptB1 (24). To date, most xpt-like genes have been found associated with lysogenic phage-like elements. In Photorhabdus strain W14, tcdA and tcdB are associated with the lysR phage gene and are expressed from phage-like promoters induced by mitomycin C (33). It has been postulated that tc genes in Photorhabdus may be part of a cryptic prophage (33). Indeed, tcdA and tcdB from P. luminescens W14 have been shown to produce phage-like particles when expressed in Escherichia coli (33). In Y. pestis, the xpt-like genes are interspersed with the lysogenic phage-like genes gp13 and gp19, which code for a porin and a lysine, respectively. In Serratia, although the sep genes are located on a plasmid (8), they are associated with a gp19-like gene. In X. nematophila PMFI296, the xpt genes are not associated with lysogenic phage genes and were perhaps horizontally acquired by other means. The aim of this work was to characterize United Kingdom Xenorhabdus isolates by sequencing key housekeeping genes and to determine whether typing schemes correlated to patterns of insecticidal activity. For selected strains the aim was to examine the context of xpt genes and recombination at housekeeping loci and fliC in order to give an indication of the amount of potential gene exchange in the populations.
|
|
|---|
Bacterial strain isolation.
Xenorhabdus strains were isolated from nematodes obtained from over 20 locations throughout the United Kingdom by an entrapment method using Galleria melonella larvae. Soil samples were collected and stored at 18°C for a maximum of 21 days. Leaf litter and large invertebrates were removed from the soil by hand, and approximately 30-g samples were placed in vented plastic sandwich boxes and baited with four G. melonella larvae. After 3, 7, and 14 days, the G. melonella larvae were observed for mortality and hemolymph from dead larvae was examined under a microscope. If rod-shaped bacteria were identified, the hemolymph sample was plated onto NBTA medium and incubated at 25°C for 72 h. To obtain pure cultures, dark blue and green colonies were selected and subcultured onto fresh NBTA plates, from which single colonies were picked and plated onto NA. In order to differentiate between the presumptive Xenorhabdus and Photorhabdus strains, they were grown for 16 h in 5 ml LB at 28°C. One milliliter of this culture was then placed in a luminometer (Wallac), and those that gave a positive signal (emitted light) were presumed to be Photorhabdus spp. and were not included in the study.
DNA manipulation techniques.
Chromosomal DNA was extracted from 1.5 ml of overnight X. nematophila culture by using a genomic DNA isolation kit (QIAGEN, Crawley, United Kingdom). PCR products were gel purified using a QIAquick PCR purification kit (QIAGEN) and products cloned using a pGEM Easy TA kit (Promega) according to the manufacturer's instructions. Cosmids and PCR products were sequenced using PCR, subcloning, and primer walking.
Sequencing of housekeeping genes and flagellin gene.
An internal fragment of the 16S rRNA, asd (accession no. AY317151.1), ompR (U07746.1), recA (U87924.1), serC (AY077462.1), fliC (X91047.1), and xptA3 genes was amplified by PCR. Each reaction was carried out in a 50-µl volume containing 1 µM of each primer (Table 1), 0.2 mM deoxynucleoside triphosphates, 200 ng chromosomal DNA, and 1 U Dynazyme (Flowgen). PCR cycles consisted of 94°C for 2 min and 25 cycles of 94oC for 15 s, 50 to 55oC (see Table 1 for value for each primer set) for 30 s, and 72°C for 45 s. The PCR products were purified and 200 ng of purified DNA used for sequencing with the same primers that had been used for the PCR amplification. Sequences of individual genes were edited using EditSeq (DNA*). Multiple-sequence alignments were generated using ClustalX and reedited to ensure correct alignment within Megalign. Sequences were compared using the packages Dnadist, Dnapars, Seqboot, Consense, and Neighbor within the Phylip suite of packages. Trees were viewed using Treeview.
|
View this table: [in a new window] |
TABLE 1. Primer sequences, expected product sizes, and annealing temperatures used in PCR amplification of the 16S rRNA gene, housekeeping genes, fliC, and xptA3
|
sawyer/geneconv) with silent polymorphisms only (r option). SplitsTrees were constructed using SplitsTree4 (http://www-ab.informatik.uni-tuebingen.de/software/jsplits) with the default settings. Bootscanning was implemented using Simplot 2.5 (http://sray.med.som.jhmi.edu/RaySoft/SimPlot) with the default settings, with the exception of a window size of 120 bp and the Dnapars algorithm.
Insect bioassays.
The activity of each strain against Pieris brassicae and Phaedon cochleariae was measured as previously described (24). The activity of each strain against Aedes aegypti was measured as described by Thomas et al. (32). All bioassays were carried out in duplicate, and positive insecticidal activity was determined as greater than 30% insect mortality within 3 days. Negative control samples had less than 5% mortality in repeat assays.
Southern blotting.
Chromosomal DNA (1 µg) was restricted and the fragments separated on a 0.7% agarose gel. The gel was blotted onto a positively charged nylon membrane (Roche) using 20x SSC (3 M NaCl, 300 mM Na-citrate, pH 7.0). Membranes were prehybridized in hybridization solution (5x SSC, 0.5% Roche blocking agent,1% N-lauryl sarcosine, 0.02% sodium dodecyl sulfate [SDS]) for 2 h at 68°C and hybridized in the same solution containing 20 to 50 ng µl1 of gel-purified labeled probe. Probes (approximately 500 bp in length) were generated by PCR using 3 µl digoxigenin DNA-labeling mix, 1 µM of each primer, 1 U Dynazyme (Flowgen), and 0.5 µg template. To create probes, a template consisting of the cloned PCR product containing the IS element was used. After hybridization, membranes were washed in 2x SSC containing 0.1% (wt/vol) SDS for 5 min at room temperature and then twice with 0.1x SSC containing 0.1% (wt/vol) SDS at 68°C for 15 min. Membranes were blocked in Tris (100 mM, pH 7.0) containing NaCl (150 mM) and 1% (wt/vol) blocking agent (Roche) for 1 h, before addition of antidigoxin-AP (Roche) as described by the manufacturers. Membranes were developed using nitroblue tetrazolium/5-bromo-4-chloro-3-indolylphosphate (Sigma, Gillingham, United Kingdom).
PCR mapping of xpt loci.
PCR was carried out using 1.7 U of the enzyme Expand HiFi (Roche) with 1 µM primer, 200 µM deoxynucleoside triphosphates, 0.5 µg DNA, and the recommended buffer. Cycling conditions consisted of 2 min at 94°C; 10 cycles of 94°C for 15 s, 60°C for 30 s, and 68°C for 8 min; 20 cycles of 94°C for 15 s, 60°C for 30 s, and 68°C for 8 min with a 5-s increment for each subsequent cycle; and 68°C for 10 min. PCR products were analyzed by gel electrophoresis on 0.7% agarose gels and characterized by restriction digestion and sequencing where appropriate.
|
|
|---|
![]() View larger version (19K): [in a new window] |
FIG. 1. Phylogenetic tree showing distances based on the 16S rRNA sequences obtained from isolates within insect-pathogenic nematodes that infected G. melonella. The analysis was carried out using the programs Dnadist (Jukes-Cantor maximum likelihood) and Neighbor (unweighted-pair group method using average linkages). The scale bar indicates 0.1% estimated sequence divergence. Sequences homologous to that of P. stuartii are indicated by asterisks. The strain shown is one representative of the group whose sequences differ by less than 2 bp.
|
|
View this table: [in a new window] |
TABLE 2. Genetic diversity at the four housekeeping loci in X. bovienii strains
|
|
View this table: [in a new window] |
TABLE 3. Allele type and insecticidal activity for Xenorhabdus strains
|
|
View this table: [in a new window] |
TABLE 4. Measure of recombination with homoplasy and Sawyer's test
|
![]() View larger version (13K): [in a new window] |
FIG. 2. Split decomposition analysis of housekeeping genes with SplitsTree. By this method, evolutionary data are canonically decomposed into a sum of weakly compatible splits and are represented by a splits graph. A, ompR; B, serC; C, asd; D, recA. Each scale bar represents the percent difference among the sequence types. The divergent asd and ompR alleles from strains X. bovienii GS and WSX97a were omitted from the analysis.
|
fliC analysis.
From the 193 strains sequenced, there were 41 unique fliC sequences (alleles), which were aligned using the ClustalW algorithm. The fliC gene showed a typical predicted flagellin protein structure, with fairly conserved NH2- and COOH-terminal sections and a highly variable internal section. The aligned sequences were analyzed by Sawyer's test and gave a maximum score of 75.8 (P = 0.0001), implying extensive recombination. Indeed, a total of 216 inner fragments (recombination events) were identified as being significant, 93 of which were still significant when Bonferroni corrected. The recombination events that were detected were throughout the fragment and were not localized to the internal highly variable region. Although it is unlikely that all putative recombination events are real, the actual number is still likely to be large, leading to a highly mosaic structure for this gene. This mosaic structure was further supported by analyzing the fliC alleles by using the bootscanning method (21) that is employed by the program Simplot. For such analysis, the 42 aligned fliC alleles were split into 12 groups, with alleles in each group showing no greater than 5% divergence. Using the bootscanning procedure, the nucleotide similarity of the consensus for each group was compared against the consensus of all the other groups (reference sequences) across the entire length of the gene. The resulting graph shows similarity between the sequence being studied and the reference sequences. Figure 3 shows examples of alleles that are composed of fragments from two (Fig. 3A) or three (Fig. 3B and C) other alleles. Despite this frequent recombination leading to complex mosaic genes, four alleles dominated, contributing to over 55% of the 193 strains. Presumably, despite frequent recombination, only the successful combinations survive and predominate. In order to study the origin of the highly variable central regions, the internal regions from position 377 to 722 of the aligned sequences were analyzed separately. There were 26 unique sequences, which could be split into 11 groups consisting of sequences showing less than 5% divergence. Only two groups showed any significant identity (e < 0.01) to any protein sequence in the database by BLASTX searching. Group 5 showed 32% identity to phase 1 flagellin from Salmonella, and group 9 showed 39% identity to a flagellin from Photorhabdus.
![]() View larger version (27K): [in a new window] |
FIG. 3. Bootscanning of nucleotide similarity of different fliC alleles, highlighting their mosaic structure. The graphs show examples of the similarity between the consensus sequence being studied (query sequence) and all other consensus sequences (reference sequences) across the entire length of the fliC gene PCR product. The query consensus sequences for each graph were composed of the following alleles: A, fliC27 and -40; (B) fliC4,-15, -16, and -29; and (C), fliC10 and -39. The reference sequences shown were composed of the following alleles: 1, fliC20, -21, -31, -33, -34, and -41; 2, fliC22 and -30; 3, fliC10 and -39; 4, fliC3, -8, -11, -12, -32, -35, and -36; 5, fliC1, -7, -13, -18, -19, and -23; 6, fliC4, -15, -16, and -29; 7, fliC3, -8, -11, -12, -32, -35, and -36; and 8, fliC5.
|
Characterization of the xpt PAI.
The xpt genes present in cHRIM1 (13) were detected by PCR in all seven X. nematophila strains. The area surrounding these genes was mapped in all seven strains by carrying out a combination of PCR and sequencing. The region to the right of the xpt loci in all strains shows the classical features of a pathogenicity island (PAI), with an integrase gene adjacent to the leucine (leu) tRNA gene and the core chromosomal genes uvrC and pgsA. This region is a common site for the insertion of PAIs and bacteriophages in many bacteria.
Upstream of the xptA to -C genes were the potential insect toxin genes xptD1 and xptE1 and a second exochitinase gene, chi2 (Fig. 4). Further upstream of this was a region of DNA that when translated showed homology to a peptide synthetase fragment and associated ATP binding cassette (ABC transporter) gene, homologues of which have been postulated to pump out peptide antibiotics from the cell (16). Beyond this were the core chromosomal genes fadR, nhaB, and dsbB, suggesting that this marked the right-hand junction (RHJ) of the PAI. The region containing the xpt genes and chitinase genes appeared, by PCR analysis, to be conserved in all seven X. nematophila strains, but small rearrangements or point mutations could not be discounted. Other regions, however, showed differences; at the RHJ, strain 0014 contained an approximately 500-bp deletion around the start of the helicase gene, and PMFI296 contained a region not present in the other strains, which was associated with a tnpA-like transposon. Deletions of parts of the peptide synthetase gene had occurred in both the 9965 and PMFI296/0014 strains (Fig. 4).
![]() View larger version (10K): [in a new window] |
FIG. 4. Structures of the xpt presumptive pathogenicity island in X. nematophila isolates, highlighting differences between the strains. Variations (a to c) in the left- and right-hand junctions of strains 9965, 0014, 0015, N38, L34, QQ2, and PMFI296 are highlighted. The positions of IS630-1 in strain PMFI296 and ISAS-1 in strains 0015, QQ2, N38, and L34 are indicated on the map line. The scale is in bp. , core chromosomal genes; , transposon-like genes; , other genes. Dotted lines show areas of deletion.
|
![]() View larger version (8K): [in a new window] |
FIG. 5. Comparison of the xpt loci of X. nematophila PMFI296 to those of X. bovienii GS. Areas of colinearity are marked with a dotted line. , core chromosomal genes; , other genes. Note that the homology between the strains extends beyond the core chromosomal genes into the phage-like genes at the RHJ of the putative xpt PAI.
|
![]() View larger version (78K): [in a new window] |
FIG. 6. Southern hybridization of restriction-digested bacterial chromosomal DNA with digoxigenin-labeled IS probes. (a) HindIII-digested DNA probed with ISAS-1; (b) EcoRI-digested DNA probed with IS630-1; (c) HindIII-digested DNA probed with IS1388; (d) HindIII-digested DNA probed with IS630-2; (e) HindIII-digested DNA probed with TnpA-1; (f) HindIII-digested DNA probed with TnpA-2; (g) HindIII-digested DNA probed with ISPSY-11; (h) HindIII-digested DNA probed with Xis. X. nematophila isolates: 1, L34; 2, N38; 3, QQ2; 4, 0015; 5, 0014; 6, PMFI296; 7, 9965. Markers are DNA cut with HindIII.
|
|
View this table: [in a new window] |
TABLE 5. Insertion sequences present in the xpt PAI of X. nematophilus strains
|
![]() View larger version (11K): [in a new window] |
FIG. 7. Structure of the xpt genes in X. bovienii isolate I73, highlighting classical features of a type 2 transposon. The same loci on the X. bovienii and P. luminescens chromosomes are also shown for comparison. The sequence of the predicted point of insertion (based on the X. bovienii GS sequence) is shown, with the duplicated target sequence highlighted in boldface. The scale is in bp. , core chromosomal genes; , transposon-like genes.
|
|
|
|---|
The four housekeeping loci examined in this study proved to be ideal candidates for use in MLST analysis, having low dN/dS ratios and only moderate variation (around 5% polymorphic sites). Even with analysis of only four loci, STs of the strains appeared to be more informative than analysis of 16S rRNA sequences alone. For example, the ST of X. bovienii strain GS, obtained from a novel nematode species in the American midwest (29), showed, as expected, that it was distinct from the strains isolated in the United Kingdom. However, analysis of the 16S rRNA of this strain could not differentiate it from the United Kingdom isolates. The United Kingdom isolate strain WSX97a also had a common 16S rRNA sequence but unique alleles at all four loci and therefore probably represents a rare United Kingdom strain, which would be interesting for further study. Surprisingly, the 16S rRNA sequences of strains NRF5e, BRF26b, and MER25b were distinct from those of all the other X. bovienii strains and showed closest homology to the 16S rRNA gene of P. stuartii. In contrast to this, housekeeping alleles of these three strains were common among the X. bovienii strains examined. Indeed, strain BRF26b shared an ST with nine other strains, and even the highly variable fliC gene in the three strains was identical to those in other X. bovienii strains. Therefore, it seems likely that the 16S rRNA gene had been subject to recombination with an external source, and this highlights the problems of typing bacteria based on a single gene.
Statistical evidence indicated that recombination had occurred at one of the four loci (serC), and the ISA of 0.256 suggested that only moderate interallele recombination had taken place. This is consistent with the fact that opportunity for genetic exchange would be limited due to the life cycle of the bacterium, which resides in monoculture in its host nematode and produces potent antimicrobial substances during insect infection, again limiting contact with other bacteria. Despite this, however, extensive recombination had taken place at the fliC locus, resulting in a highly variable central region, which corresponds to the exposed region of the flagella. This is similar to the situation with the fliC genes of bacterial pathogens which infect higher organisms and is thought to be a result of an "arms race" between the bacteria and the adaptive immune response of the host. However, for an insect pathogen such as X. bovienii, this is surprising, since insects possess only an innate immune system similar to the innate systems of higher organisms. The innate immune response is usually triggered by specific conserved patterns in bacteria. Indeed, Toll-like receptor 5, a protein responsible for triggering the innate response to bacterial infection, interacts with a conserved 13-amino-acid motif in the bacterial flagella (28). Therefore, the recognition of conserved domains in proteins is unlikely to lead to the diversity exhibited by the fliC genes of the X. bovienii isolates. Hence, other pressures must be at work in order to create this diversity. Similar diversity is seen in the fliC genes of Pseudomonas spp. (2). Like Xenorhabdus, most Pseudomonas species are not pathogens of higher organisms, and therefore it was hypothesized that the variation was a result of evading infection by bacteriophages. This could also be a possible explanation for the diversity seen in Xenorhabdus. All of the highly variable regions of the fliC genes examined showed little homology to any sequence in the database, implying that they were acquired from as-yet-undiscovered sequences or that frequent recombination had led to sequences that were no longer discernible. Acquisition of DNA from diverse sources in Xenorhabdus was also implied by the presence of IS elements that were similar to those from A. eutrophus and A. salmonicida in X. nematophila strains and a 16S rRNA gene similar to that of P. stuartii in three X. bovienii strains.
This work has also shown that the xpt genes are present on a PAI associated with a leu tRNA. The homologous tca genes in Photorhabdus have also recently been assigned to a putative PAI by virtue of the facts that it is unstable, contains phage-like genes, and is associated with aspV tRNA (25). However, the xpt PAI represents a more classical PAI with similarities to the high PAI of Yersinia spp. (20). Both the xpt PAI and the high PAI share a P4-like integrase at the right-hand junction, a core region containing pathogenicity genes, and an AT-rich left-hand junction containing many IS elements that has a tendency to delete (20). The left-hand junction of the PAI is likely to be between the ABC transporter gene and the dsbB gene. However, there is no direct repeat corresponding to the 3' end of the leu tRNA present in this region, which usually marks the end of a PAI. It is most likely that deletions have occurred in this area, removing the repeat, as seen with the Yersinia high PAI in E. coli (20). It has been postulated that PAIs are evolutionary relics, where, after acquisition and over time, functions such as mobility are gradually lost as the island becomes incorporated into the chromosome (6). Comparing the PAI locus with that of X. bovienii GS appears to support this, as phage-like genes in X. bovienii appear to have been replaced by xpt genes in X. nematophila as the locus is converted from a phage to a PAI. The presence of the peptide synthetase and associated ABC transporter in the PAI would give the bacterium and presumably its nematode host a competitive advantage by killing competing bacteria in the insect cadaver. It is interesting that novel colicin and immunity proteins are also found associated with insecticidal tc genes in Photorhabdus spp. (25).
In the PAI all genes are tightly grouped into two insect-pathogenic regions, a mobility region and a region putatively involved in peptide antibiotic synthesis and resistance, which are all separated by AT-rich junctions (Fig. 4). In Photorhabdus, the tca genes are associated with lysogenic phage genes, and therefore it is conceivable that these genes originated from an integrated phage. This work has shown that the xpt genes in X. nematophila are not associated with phage genes but show close linkage to IS elements and therefore may have originated on transposons. Thus, unlike the tc genes of Photorhabdus, the xpt genes are unlikely to be expressed from phage-like promoters or form phage-like particles. The xpt loci found in X. bovienii I73 would appear to support the hypothesis that the xpt genes originated from transposons. The IS elements that flank the I73 xpt genes are homologous to members of the IS6 family, which have been found to constitute composite transposons containing antibiotic resistance (14). Interestingly, the locus includes accessory regions outside of this putative transposon. A possible theory to explain the left-hand accessory region is that "local hopping" of an IS element (19) from the original transposon to a position further downstream would generate a larger transposon, capturing the intervening region. The original IS, now internal in the transposon, would no longer be required for mobility and could lead to degeneracy, as is the case with the two internal IS6 elements at this locus. The extra plasmid-like region to the right of the putative transposon could also have been acquired in a similar manner. Although no IS6-like element exists at the RHJ, a partial tnpA-like element at this location may contribute the necessary transpositional functions. To date IS elements in Xenorhabdus have not been studied. This work has shown that X. nematophila strains, as well as containing common Enterobacteriaceae IS elements such as IS630 and IS1388, also contain elements that are rarer in Enterobacteriaceae, such as ISAS1 which shows over 75% homology to IS elements in A. salmonicida yet only 33% to any IS currently identified in Photorhabdus. The many different types of IS elements present in X. nematophila strains and their high copy number would indicate a very fluid genome, perhaps playing a significant role in the evolution of strains. This is surprising for a bacterium that resides in the relatively sheltered environment of the nematode gut.
In the xpt PAI, most IS elements detected were present in the AT-rich junctions that flanked the groups of insect virulence and mobility genes. These areas may represent hot spots for insertion, thus allowing the PAI to evolve with the acquisition of new genes without the disruption of existing loci. Indeed, the presence of IS elements in these regions may enhance expression of the xpt loci. Strain PMFI296 had the highest insecticidal activity against P. brassicae, and expression of the xptA1 gene in this strain was shown to be approximately 10-fold greater than that of strain 9965 (data not shown). Since PMFI296 was the only strain to contain an IS630-like element downstream of the xptA1 gene, it seems likely that the higher insecticidal activity and increased expression of xptA1 is likely to be due to the presence of the IS630-like element. It is also possible that the presence of the ISAS1-like element in strains 0015, N38, QQ2, and L34 may increase expression of the upstream xptA2 operon. Since xptA1 and xptA2 have been shown to be involved in toxicity to different insects (24), the integration of IS elements in the xpt PAI may modulate pathogenicity towards different insect species and thus tailor the bacterium/nematode complex to inhabit slightly different niches.
The differences, even among closely related strains, in insecticidal activity of the X. bovienii strains observed in this study can be explained not only by the presence of these genes on mobile genetic elements but also by frequent loss of gene function, as indicated by the frameshift and IS insertion in the xptA3 genes detected in X. bovienii strains. Frameshift mutations were also observed in analogous genes in Y. pestis (17). Therefore, the mobility of these genes and the spectrum of activity of different xpt genes (24) would imply that, in a manner analogous to that for B. thuringiensis cry genes, a whole raft of xpt genes exhibiting different biological activity exist in different Xenorhabdus strains. Therefore, further analysis of novel Xenorhabdus strains is required in order to isolate xpt genes that target other important insect pests.
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»