Skip to main page content

• HOME
• Current Issue
• Archive
• Alerts
• About ASM
• Contact us
• Tech Support
• Journals.ASM.org

This Article Abstract Full Text (PDF) Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Reprints and Permissions Copyright Information Books from ASM Press MicrobeWorld Citing Articles Citing Articles via Google Scholar Google Scholar Articles by LeBlanc, P. M. Articles by Filion, M. Search for Related Content PubMed PubMed Citation Articles by LeBlanc, P. M. Articles by Filion, M. Agricola Articles by LeBlanc, P. M. Articles by Filion, M.

 Previous Article  |  Next Article 

Applied and Environmental Microbiology, July 2007, p. 4128-4134, Vol. 73, No. 13
0099-2240/07/$08.00+0     doi:10.1128/AEM.02590-06
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

Alteration of Soil Rhizosphere Communities following Genetic Transformation of White Spruce

Philippe M. LeBlanc,¹ Richard C. Hamelin,² and Martin Filion^1*

Université de Moncton, Department of Biology, Moncton, NB, Canada E1A 3E9,¹ Laurentian Forestry Centre, Natural Resources Canada, 1055 du PEPS Street, Sainte-Foy, QC, Canada G1V 4C7²

Received 7 November 2006/ Accepted 21 April 2007

Top ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION REFERENCES

 
The application of plant genetic manipulations to agriculture and forestry with the aim of alleviating insect damage through Bacillus thuringiensis transformation could lead to a significant reduction in the release of pesticides into the environment. However, many groups have come forward with very valid and important questions related to potentially adverse effects, and it is crucial to assess and better understand the impact that this technology might have on ecosystems. In this study, we analyzed rhizosphere soil samples collected from the first B. thuringiensis-transformed trees [with insertion of the CryIA(b) toxin-encoding gene] grown in Canada (Val-Cartier, QC, Canada) as part of an ecological impact assessment project. Using a robust amplified rRNA gene restriction analysis approach coupled with 16S rRNA gene sequencing, the rhizosphere-inhabiting microbial communities of white spruce (Picea glauca) genetically modified by biolistic insertion of the cryIA(b), uidA (beta-glucuronidase), and nptII genes were compared with the microbial communities associated with non-genetically modified counterparts and with trees in which only the genetic marker genes uidA and nptII have been inserted. Analysis of 1,728 rhizosphere bacterial clones (576 clones per treatment) using a Cramér-von Mises statistic analysis combined with a Monte Carlo comparison clearly indicated that there was a statistically significant difference (P < 0.05) between the microbial communities inhabiting the rhizospheres of trees carrying the cryIA(b), uidA, and nptII transgenes, trees carrying only the uidA and nptII transgenes, and control trees. Clear rhizosphere microbial community alterations due to B. thuringiensis tree genetic modification have to our knowledge never been described previously and open the door to interesting questions related to B. thuringiensis genetic transformation and also to the impact of commonly used uidA and nptII genetic marker genes.

Top ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION REFERENCES

 
The application of plant genetic manipulation to agriculture and forestry with the aim of increasing disease resistance or alleviating insect damage could lead to a significant reduction in the release of pesticides into the environment (27). However, considering the relatively recent development of these techniques and the lack of sufficient reliable data derived from environmental impact studies, the possible effects of transgenic DNA insertion into a plant's genome are yet to be fully evaluated. Very valid and important questions related to potentially adverse effects have been raised, and it is crucial to assess and better understand the impact that this technology might have on ecosystems.

Complex microbial communities living in close proximity to plant roots, better known as the rhizosphere, are sustained by the loss of organic matter through root exudation. It has been suggested that rhizospheres could be altered in response to plant genetic transformation (19, 28, 33). Comparative studies assessing whether there are differences between microbial communities living in the rhizospheres of genetically modified and non-genetically modified plants represent an important first step in determining if the presence of transgenic material can catalyze changes in the environment. Although a few studies have reported slight rhizosphere microbial diversity changes associated with transgenic plants compared with nontransgenic plants, thus far the differences have been either very minor, mainly attributed to environmental factors such as seasonal variation or field site (9-12, 14, 20, 21, 24, 32, 34, 42), or simply not statistically significant (13, 15, 23, 24, 29-31, 38, 40).

The forest industry represents an important economic sector in North America and in particular in Canada. Genetic improvement for wood quality and pest resistance has been used for decades, and genetically improved material has been widely deployed (44). More recently, genetic engineering has been applied to tree improvement as a shortcut to traditional breeding. The first documented genetic modification of a conifer was in 1991 (25). One of the particular interests in genetically modifying trees, among many other crops, has been Bacillus thuringiensis transformation to confer insect resistance through the insertion of cry genes. The expressed B. thuringiensis toxins demonstrate a remarkable specificity for their target host spectra, and nonconiferous plant systems expressing these toxins have had no significant effect on nontarget fauna, including plant pathogen populations (11), earthworms and nematodes (38), microbial activity (7), bacterial community structure, or total soil bacterial and fungal counts (7, 11, 38). But B. thuringiensis recombinant DNA and expressed Cry proteins were shown to be released into soils through root exudation and plant tissue decomposition, where they can remain intact and chemically active for extended periods of time (up to 240 days for the proteins) (36, 45, 48). This could be particularly important in conifers expressing a Cry toxin since release could occur over several decades.

The present study reports a survey of microbial biodiversity in an experimental plot where genetically modified white spruce (Picea glauca) and control white spruce were planted in 2000. Trees were genetically modified by biolistic insertion of a B. thuringiensis transgene [cryIA(b)] and the genetic marker genes uidA (encoding beta-glucuronidase [GUS]) and nptII [cryIA(b)/uidA/nptII trees] or only the marker genes uidA and nptII (uidA/nptII trees) (37). The purpose of insertion of the cryIA(b) transgene is expression of a delta-endotoxin which targets Lepidoptera insects, including the spruce budworm (Choristoneura fumiferana). Following ingestion and protease activation, the CryIA(b) toxin induces lesions in the epithelial walls of the insect's digestive tract, causing starvation and death (27). The trees are part of an environmental impact study to assess the impact of genetically modified trees. The objective of this study was to investigate whether changes in the microbial communities inhabiting the rhizosphere of white spruce are associated with the presence of transgenic constructs expressing the Cry protein or only markers. A robust culture-independent amplified rRNA gene restriction analysis (ARDRA) and a 16S rRNA gene sequencing-based approach were used to assess microbial composition and diversity.

Top ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION REFERENCES

 
Genetic transformation of white spruce.
Transgenic white spruce trees containing either cryIA(b), uidA, and nptII or only uidA and nptII were generated at Natural Resources Canada's Laurentian Forestry Centre in Quebec, QC, Canada (A. Seguin, unpublished data). Two embryogenic lines of white spruce were initiated from mature zygotic embryos using a procedure described by Cheliak and Klimaszewska (3). Maintenance and maturation of the tissues were performed as described previously (41). The synthetic construct that was described by Cheng et al. (4) was used in the present study. In the vectors the cry1A(b) gene is under the control of the CaMV 35S promoter (pKSB) or the maize ubiquitin promoter (pKUB) and is completed by the termination signal of the nopaline synthase gene. For selection and screening of the transgenic lines, the pBI426 vector, which contains a bifunctional fusion between the uidA and nptII genes, was used (6). All plasmid DNA were isolated and purified using CsCl gradient centrifugation as previously described (2). Microprojectile bombardment was carried out using a PDS1000 helium-driven system apparatus (Bio-Rad, Mississauga, ON, Canada) as described by Charest et al. (2), except that 0.6-µm gold particles were used. For cotransformation, plasmids were mixed together before adsorption onto the gold particles, resulting in 350 ng of DNA for each bombardment. In brief, 4- to 5-week-old cotyledonary stage somatic embryos (approximately 30 to 40) were picked and placed in the center of a petri dish containing induction/maintenance medium (41). The following day the embryos were bombarded using the following settings: 900-lb/in² rupture disk and 0.25-in gap. After bombardment, the somatic embryos were maintained on the induction medium for 10 days before they were transferred to selection medium (15 mg/liter kanamycin in maintenance medium). Kanamycin selection was maintained for up to 22 weeks with subculturing every 2 to 3 weeks. Reinduced embryonal masses (referred to as secondary somatic embryogenesis) were screened for GUS activity by histochemical staining. Emerging embryonal masses with GUS-positive sectors were cultivated independently (41). Regeneration of plants from stably transformed cultures was initiated approximately 4 to 6 months after DNA transformation.

Sampling.
Rhizosphere soil samples were collected in November 2003 from an experimental plot in Val Cartier (QC, Canada) where genetically modified and control white spruce trees were planted in 2000. The experimental design comprised three treatments: (i) trees transformed with cryIA(b), uidA, and nptII; (ii) trees transformed with uidA and nptII; and (iii) nontransformed controls. For the cryIA(b)/uidA/nptII and uidA/nptII treatments, replicate trees used in this study were derived from the same transformation event. The treatments were arranged in a completely randomized design. Soil samples (approximately 15 to 20 g [wet weight] each) were collected from the rhizospheres of six different trees (n = 6) to obtain a total of 18 soil samples. The samples were transported on ice to the laboratory and then lyophilized and cryopreserved at –80°C until total genomic soil DNA isolation.

Soil DNA isolation.
For each sample, DNA was extracted from 250-mg (dry weight) soil subsamples using a MoBio UltraClean soil DNA isolation kit by following the manufacturer's protocol (MoBio Laboratories, Solana Beach, CA). Isolated DNA was resuspended in 50 µl of S5 solution (MoBio Laboratories) and maintained at –20°C until 16S rRNA gene PCR amplification.

16S rRNA gene PCR amplification.
Soil DNA was PCR amplified using primers 968f and 1401r targeting a 433-bp portion of the 16S rRNA gene (35). Five microliters of 10x PCR buffer (QIAGEN, Mississauga, ON, Canada), 5 µl of a 5 µM solution of each primer, 1 µl of a 10 mM solution of each deoxyribonucleoside triphosphate, 2 µl of a 1/20 dilution of total soil DNA, 1.25 U of Taq DNA polymerase (QIAGEN), and 31.75 µl nanopure H₂O were combined to obtain a 50-µl (total volume) PCR mixture for each sample. PCR amplifications were performed using an MJ Research PTC-200 DNA Engine (MJ Research, Waltham, MA) with the following parameters: (i) initial denaturation for 3 min at 94°C; (ii) 35 cycles of denaturation for 30 s at 94°C, hybridization for 30 s at 60°C, and elongation for 45 s at 72°C; and (iii) final elongation for 5 min at 72°C. PCR products were then visualized using 1.5% agarose-1x Tris-acetate-EDTA gel electrophoresis, followed by ethidium bromide staining and UV transillumination using a Gel Logic 200 transilluminator (Eastman Kodak Company, Rochester, NY). Once amplicon sizes were confirmed, DNA cleanup was performed using a PCR cleanup kit (QIAGEN), and DNA was conserved at –20°C prior to cloning.

Shotgun cloning.
Cleaned amplicons were ligated into pCR4-TOPO vectors and inserted into chemically competent Escherichia coli cells as described in the One Shot TOPO-TA cloning kit protocol (Invitrogen, Carlsbad, CA). For each sample, 96 clones were randomly selected and cultured for 24 h at 37°C in 1.5 ml of Luria-Bertani broth supplemented with 100 µg/ml ampicillin.

Plasmid DNA extraction and ARDRA.
Plasmid DNA extraction was performed using a DirectPrep 96 miniprep kit (QIAGEN), and partial 16S rRNA gene inserts were PCR amplified using the PCR and protocols described above but with a 1/20 plasmid DNA dilution and PCR primers T3 and T7. After confirmation of amplicon sizes using gel electrophoresis, 13 µl of amplified DNA from each clone was separately digested for 2 h at 37°C using either 5 U of MspI or 5 U of HaeIII (New England Biolabs, Mississauga, ON, Canada) in 20-µl reaction mixtures using the appropriate restriction buffer and concentration, as suggested by the manufacturer (New England Biolabs). Restriction patterns were visualized by 3% agarose-1x Tris-boric acid-EDTA electrophoresis, followed by ethidium bromide staining and UV transillumination. Each different restriction pattern among the 1,728 patterns obtained was identified using the 1D Image Analysis software (v. 3.6; Kodak) and was defined as an operational taxonomical unit (OTU).

Sequencing.
Aliquots (1.5 µl) of extracted plasmidic DNA from at least one clone of each different OTU were sequenced with a CEQ 8000 genetic analysis system (Beckman Coulter, Fullerton, CA) using primer T7 according to the manufacturer's instructions.

Sequence analysis.
DNA sequences were edited, and consensus sequences were obtained using BIOEDIT v.7.0.4.1 (22). The presence of a chimera was verified using Bellerophon (http://foo.maths.uq.edu.au/huber/bellerophon.pl) (26) and CHIMERA CHECK (http://rdp.cme.msu.edu) (5), and suspected chimeric sequences were not included in further analysis. All the remaining sequences were screened against the sequences in GenBank (National Center for Biotechnology Information) using BLASTn v.2.2.12 (http://www.ncbi.nlm.nih.gov/BLAST/) (1) and RDP-II (5). Homologous sequences found in GenBank and RDP-II were used to construct a multiple-sequence alignment with ClustalX v.1.83 (46). A neighbor-joining analysis was performed using PAUP v.4.0b10 (43) with 1,000 bootstrap replicates. Phylogenetic trees were visualized using Treeview and were edited in Microsoft Word.

Statistical analysis.
Distance matrices were prepared using DNADIST and the Jukes-Cantor algorithm (PHYLIP v.3.64) (16). In order to test the null hypothesis that there are no difference in the three treatments, the Cramér-von Mises statistic was calculated, and a Monte Carlo test procedure was executed using the S-Libshuff program (v.1.0; http://www.plantpath.wisc.edu/fac/joh/S-libshuff.html) (39). The tests were conducted first for the three different treatments and then for all six individual rRNA gene libraries within each treatment. OTU frequency was determined by grouping all clones assigned to the same OTU. Richness, evenness, and dominance indices were calculated either manually or using EstimateS (v.7.5.05; http://purl.oclc.org/estimates). The indices retrieved included richness indicators such as the Shannon diversity index, the Chao-1 rarefaction estimator, the Margalef index, and the Alpha diversity and Shannon and Simpson evenness/dominance indices.

Nucleotide sequence accession numbers.
The nucleotide sequences determined in this study have been deposited in the NCBI database under accession no. DQ683995 to DQ684671.

Top ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION REFERENCES

 
High-molecular-weight DNA was recovered from every soil sample, and the 16S ribosomal cluster was successfully amplified by PCR. Restriction enzymes MspI and HaeIII were used to screen the 1,728 cloned amplicons for ARDRA. OTU grouping and sequencing based on restriction pattern were previously shown to be reproducible (17).

ARDRA was carried out on 1,728 shotgun-cloned partial 16S rRNA gene fragments (96 fragments per replicate). Restriction pattern data sets revealed 686 OTUs which were successfully sequenced, most of which had only few representative clones (Fig. 1). Of these OTUs, 314 were associated with the rhizospheres of cryIA(b)/uidA/nptII trees, 294 were associated with the rhizospheres of uidA/nptII trees, and 270 were associated with the rhizospheres of control trees. A total of 1,593 environmental clones were used for BLAST searches in the GenBank and RDP-II databases. This allowed us to classify unknowns into major bacterial groups. A high level of diversity was observed, and 16 major groups were represented (Table 1). The most frequently represented groups detected in all treatments were the Gammaproteobacteria (30.73 to 34.36% for the three treatments), Thermomicrobia (4.30 to 7.90%), Betaproteobacteria (4.67 to 6.13%), Chlamydiae (5.80 to 7.60%), and Alphaproteobacteria (5.00 to 7.20%). Within groups specifically associated with certain treatments, Bacteroidetes were found in only one replicate of the cryIA(b)/uidA/nptII treatment, whereas cyanobacteria were found in all but the uidA/nptII treatment. Genus incertae sedis WS3 and Nitrospira were also present in all treatments except the cryIA(b)/uidA/nptII treatment (Table 1). Certain OTUs were associated exclusively with certain treatments (Table 2). Twenty clones distributed among three OTUs were associated exclusively with the cryIA(b)/uidA/nptII treatment, 19 clones in three OTUs were associated exclusively with the uidA/nptII treatment, and 51 clones in seven OTUs were associated exclusively with the control treatment. A total of 27 clones in four OTUs shared association with both the cryIA(b)/uidA/nptII and uidA/nptII treatments, and 91 clones in 14 OTUs shared association with the uidA/nptII and control treatments.

View larger version (32K): [in this window] [in a new window]

FIG. 1. Observed OTU distribution for two transgenic P. glauca treatments and one nontransgenic P. glauca treatment.

View this table: [in this window] [in a new window]

TABLE 1. Relative clone abundance for the 18 different 16S rRNA gene libraries obtained from the rhizospheres of transgenic [cryIA(b)/uidA/nptII and uidA/nptII] and nontransgenic (control) trees with respect to different bacterial taxa

View this table: [in this window] [in a new window]

TABLE 2. OTUs systematically associated with certain treatments identified by BLAST nucleotide-nucleotide closest-neighbor comparisons in GenBank

To test the null hypothesis that the treatments had no effect on the microbial populations, a Cramér-von Mises statistic analysis combined with a Monte Carlo comparison was conducted. Combined sequence libraries for each of the three treatments were significantly different from each other (P < 0.05 adjusted for multiple comparisons at P < 0.0085) (Table 3). The null hypothesis was rejected for all treatment comparisons, including the cryIA(b)/uidA/nptII-control, cryIA(b)/uidA/nptII-uidA/nptII, and uidA/nptII-control comparisons. In each comparison, the test revealed significantly more differences among treatments than within treatments. Further analysis also revealed that most replicate sequence libraries for each treatment were significantly similar to each other (P < 0.05) (Table 4).

View this table: [in this window] [in a new window]

TABLE 3. Comparison of 16S rRNA gene sequence libraries from rhizospheres of transgenic P. glauca and nontransgenic P. glauca combined for treatmentsa

View this table: [in this window] [in a new window]

TABLE 4. Comparison of 16S rRNA gene sequence libraries from rhizospheres of transgenic P. glauca and nontransgenic P. glauca for replicatesa

Diversity indices were compiled for each of the three treatments (Table 5). All diversity indices increased from the control treatment to the uidA/nptII treatment to the cryIA(b)/uidA/nptII treatment.

View this table: [in this window] [in a new window]

TABLE 5. Diversity indices for bacterial OTUs associated with the rhizospheres of genetically modified and non-genetically modified P. glauca trees

Top ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION REFERENCES

 
The aims of this study were to characterize the rhizospheres of genetically modified and control P. glauca trees growing in an experimental plantation and to determine whether the treatments had an effect on microbial composition and diversity. Analysis of 1,728 rhizosphere bacterial clones (576 clones per treatment) using a highly discriminating culture-independent molecular approach indicated that there were statistically significant differences between the microbial communities inhabiting the rhizospheres of trees carrying the cryIA(b) and uidA/nptII transgenes, trees carrying only the uidA/nptII transgenes, and control trees. Inclusion of the Cramér-von Mises statistic allowed us to consider the species compositions of the rhizospheres while taking genetic diversity into account. In all cases, partial 16S rRNA gene cloned libraries were more similar within treatments than among treatments. However, before conclusions concerning the biological significance of these results are reached, the methodology and analyses should be discussed.

One can argue that the first potential weakness of this study is that the site was not replicated. We sampled the only experimental site with transgenic conifers in Canada, and therefore it was not possible to repeat our sampling at different sites. It is possible that our results represent an artifact of the site studied and cannot be generalized. However, the treatments were properly replicated within the experimental design at the study site, and the observed differences were significant for the treatments. This suggests that the observed differences were clearly associated with the treatments.

Another potential explanation for our observations is that the treatments were responsible for the observed differences but that factors associated with the treatments other than the actual transgenes were the real cause. For example, transformed P. glauca calli were screened using antibiotics to eliminate untransformed cells. Clearly, this could not be done with the untransformed callus cultures. This selection with antibiotics could have been responsible for the observed differences in microbial communities in our treatments. However, if this were the case, we would not expect to observe differences between the transformed trees containing the cryIA(b) and uidA/nptII constructs and the transformed trees containing only the uidA/nptII construct, since they were screened with the same antibiotics.

Another possibility is that the transformation events themselves, not the gene constructs, were responsible for the observed effects. Although this is entirely possible, it is difficult to imagine that this effect would be more important among treatments than within treatments. Therefore, our conclusion is that the transgenes themselves were responsible for the observed differences.

The analysis used in this study is relatively new. It is possible that the results obtained are in fact attributable to the method of analysis. Although the analysis is sound and intuitive, the differences observed could be biologically insignificant, even though they are statistically significant. We believe that this is not the case. The within-treatment comparisons were all much more similar (even though they were for different blocks) than the between-treatment comparisons. Also, the presence of multiple unique bacterial clones associated with specific treatments is indicative of a true treatment effect and could have biological importance. Unfortunately, the clones do not have close homology in the public databases for the gene sequenced. This makes it difficult to assess whether the observed trend could be extrapolated to bacterial groups.

In light of these findings, biological reasons rather than methodological reasons better support the results obtained and strongly suggest that trees within the treatments impacted the rhizosphere-inhabiting bacterial communities. Interestingly, the impact seems to be associated with both the B. thuringiensis-transformed trees, modified using the cryIA(b)/uidA/nptII gene set, and the uidA/nptII-modified trees, but the effects are different for these two treatments. The ecological diversity indices calculated also support this trend. Although this observation was surprising, there may be an explanation for it. A study conducted in 2003 established that a DNA construct composed of the nptII transgene remained detectable by PCR in potato field soil after 4 years of laboratory incubation (8). The two genetic markers used in the present study code for the production of GUS (uidA) and neomycin phosphotransferase II (nptII). They are now commonly used in plant genetic modification as means to assess the successful insertion of transgenic material, as well as to infer the level of expression. The uidA gene product catalyzes the cleavage of a wide variety of ß-glucuronides. As most plants assayed to date lack detectable glucuronidase activity, the incorporation of this gene in a plant genome allows detection of color formation in plant parts expressing this enzyme when chromogenic glucuronidase substrates are added to the growth media. The nptII gene codes for a small enzyme (25 kDa) catalyzing the phosphorylation of several antibiotics, including neomycin and kanamycin, which also allows bioassays measuring the success of a genetic modification to be conducted when suitable antibiotics are added to the growth media (18). Genetic transformation of plants with these markers results in exogenous DNA, RNA, and subsequent protein structures that are present within the plant tissue and could logically have the same potential to impact soil microbial communities as pest resistance transgenes would have.

The effects of B. thuringiensis-modified plants on soil microorganisms have been studied in other plant systems. Leaves of cotton genetically modified with B. thuringiensis [cryIA(b) and cryIA(c)] were placed in soil for decomposition, and microbial communities were monitored using community level physiological profiling and restriction fragment length polymorphism. Although some degree of variability was detected, the authors attributed the shift to the genetic modification of the plant itself rather than to the nature of the transgenic DNA (10). Plate counts revealed that minimal differences were observed when the microfloras of B. thuringiensis-producing genetically modified and B. thuringiensis (M-Trak)-treated potato plants were compared (11). The effects of B. thuringiensis [cryIA(b)]-transformed corn on microorganisms, earthworms, nematodes, protozoans, and fungi in soil were also studied. Using various techniques, no effect of the root exudate-released toxin or plant biomass on target organisms was measured (38). Studies using culture-independent approaches also yielded similar results; B. thuringiensis [cryIA(b)]-modified KeMingDao rice straw did not have significant persistent effects on the abundance of soil microorganisms when it was incubated in flooded paddy soil under laboratory conditions, although some significant differences were detected in the early stages of the experiment (47). With B. thuringiensis [cryIIIB(b)]-transformed corn there were also no differences between transgenic and nontransgenic microbial communities as determined using terminal restriction fragment length polymorphism analysis (7).

In our opinion, the main explanations supporting the clear microbial shift observed in this study compared to other studies with similar aims could be narrowed down to two possibilities. The first possibility is related to the high resolution of the methodological approach used. We used a robust culture-independent approach that accounted for the high microbial diversity usually found in rhizosphere soils, and an important number of 16S rRNA gene clones (1,593 clones in 686 OTUs) were identified and classified into major bacterial groups, which was usually not the case in most similar studies performed to date. The second possibility to explain why we found a significant difference is that rhizosphere microbial population shifts following B. thuringiensis transformation of plants have been studied before in different plant systems but not in the context of trees. However, it is difficult at this point to clearly define the potential implications for the trees, for the bacterial community outside the rhizosphere, and for soil processes in general.

Further studies are clearly needed to better understand the mechanisms involved in the microbial community structure shifts observed. Since the results obtained involved a single sampling date, it would also be interesting to track the shifts over a period of time to see if the differences are associated with certain growth stages of the plants or with seasonal variations. Microbial species or groups of interest could also be monitored quantitatively using recent advancements in real-time PCR. Finally, studies measuring functional genes to assess microbe functional diversity instead of taxonomic diversity could also add another perspective to this question.

 
This project was supported by the Canadian Biotechnology Strategy Fund.

We thank France LeBlanc and Josyanne Lamarche for technical assistance, as well as Josyanne Lamarche for input and review of the manuscript.

 
* Corresponding author. Mailing address: Université de Moncton, Department of Biology, Moncton, NB, Canada E1A 3E9. Phone: (506) 858-4329. Fax: (506) 858-4541. E-mail: martin.filion{at}umoncton.ca

Published ahead of print on 27 April 2007.

Top ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS DISCUSSION REFERENCES

 

1
1. Altschul, S. F., T. L. Madden, A. A. Schaffer, J. H. Zhang, Z. Zhang, W. Miller, and D. J. Lipman. 1997. Gapped Blast and PSI-Blast—a new generation of protein database search programs. Nucleic Acids Res. 25:3389-3402.[Abstract/Free Full Text]
2
2. Charest, P. J., Y. Devantier, and D. Lachance. 1996. Stable genetic transformation of Picea mariana (black spruce) via particle bombardment. In Vitro Cell. Dev. Biol. 32:91-99.
3
3. Cheliak, W. M., and K. Klimaszewska. 1991. Genetic variation in somatic embryogenic response in open-pollinated families of black spruce. Theor. Appl. Genet. 82:185-190.
4
4. Cheng, X. Y., R. Sardana, H. Kaplan, and I. Altosaar. 1998. Agrobacterium-transformed rice plants expressing synthetic cryIA(b) and cryIA(c) genes are highly toxic to striped stem borer and yellow stem borer. Proc. Natl. Acad. Sci. USA 95:2767-2772.[Abstract/Free Full Text]
5
5. Cole, J. R., B. Chai, R. J. Farris, Q. Wang, S. A. Kulam, D. M. McGarrell, G. M. Garrity, and J. M. Tiedje. 2005. The Ribosomal Database Project (RDP-II): sequences and tools for high-throughput rRNA analysis. Nucleic Acids Res. 33:D294-D296.[Abstract/Free Full Text]
6
6. Datla, R. S. S., J. K. Hammerlindl, L. E. Pelcher, W. L. Crosby, and G. Selvaraj. 1991. A bifunctional fusion between ß-glucuronidase and neomycin phosphotransferase: a broad-spectrum marker enzyme for plants. Gene 101:239-246.[CrossRef][Medline]
7
7. Devare, M. H., C. M. Jones, and J. E. Thies. 2004. Effect of Cry3Bb transgenic corn and tefluthrin on the soil microbial community: biomass, activity, and diversity. J. Environ. Qual. 33:837-843.[Abstract/Free Full Text]
8
8. de Vries, J., M. Heine, K. Harms, and W. Wackernagel. 2003. Spread of recombinant DNA by roots and pollen of transgenic potato plants, identified by highly specific biomonitoring using natural transformation of an Acinetobacter sp. Appl. Environ. Microbiol. 69:4455-4462.[Abstract/Free Full Text]
9
9. Di Giovanni, G. D., L. S. Watrud, R. J. Seidler, and F. Widmer. 1999. Comparison of parental and transgenic alfalfa rhizosphere bacterial communities using Biolog GN metabolic fingerprinting and enterobacterial repetitive intergenic consensus sequence-PCR (ERIC-PCR). Microb. Ecol. 37:129-139.[CrossRef][Medline]
10
10. Donegan, K. K., C. J. Palm, V. J. Fieland, L. A. Porteous, L. M. Ganio, D. L. Schaller, L. Q. Bucao, and R. J. Seidler. 1995. Changes in levels, species and DNA fingerprints of soil microorganisms associated with cotton expressing the Bacillus thuringiensis var. kurstaki endotoxin. Appl. Soil Ecol. 2:111-124.[CrossRef]
11
11. Donegan, K. K., D. L. Schaller, J. K. Stone, L. M. Ganio, G. Reed, P. B. Hamm, and R. J. Seidler. 1996. Microbial populations, fungal species diversity and plant pathogen levels in field plots of potato plants expressing the Bacillus thuringiensis var. tenebrionis endotoxin. Transgen. Res. 5:25-35.[CrossRef]
12
12. Donegan, K. K., R. J. Seidler, J. D. Doyle, L. A. Porteous, G. Digiovanni, F. Widmer, and L. S. Watrud. 1999. A field study with genetically engineered alfalfa inoculated with recombinant Sinorhizobium meliloti: effects on the soil ecosystem. J. Appl. Ecol. 36:920-936.[CrossRef]
13
13. Donegan, K. K., R. J. Seidler, V. J. Fieland, D. L. Schaller, C. J. Palm, L. M. Ganio, D. M. Cardwell, and Y. Steinberger. 1997. Decomposition of genetically engineered tobacco under field conditions—persistence of the proteinase inhibitor I product and effects on soil microbial respiration and protozoa, nematode and microarthropod populations. J. Appl. Ecol. 34:767-777.[CrossRef]
14
14. Dunfield, K. E., and J. J. Germida. 2001. Diversity of bacterial communities in the rhizosphere and root interior of field-grown genetically modified Brassica napus. FEMS Microb. Ecol. 38:1-9.[CrossRef]
15
15. Dunfield, K. E., and J. J. Germida. 2004. Impact of genetically modified crops on soil- and plant-associated microbial communities. J. Environ. Qual. 33:806-815.[Abstract/Free Full Text]
16
16. Felsenstein, J. 1989. PHYLIP—phylogeny inference package (version 3.2). Cladistics 5:164-166.
17
17. Filion, M., R. C. Hamelin, L. Bernier, and M. St-Arnaud. 2004. Molecular profiling of rhizosphere microbial communities associated with healthy and diseased black spruce (Picea mariana) seedlings grown in a nursery. Appl. Environ. Microbiol. 70:3541-3551.[Abstract/Free Full Text]
18
18. Flavell, R. B., E. Dart, R. L. Fuchs, and R. T. Fraley. 1992. Selectable marker genes: safe for plants? Bio/Technology 10:141-144.[CrossRef][Medline]
19
19. Gasson, M. J. 2000. Gene transfer from genetically modified food—commentary. Curr. Opin. Biotechnol. 11:505-508.[CrossRef][Medline]
20
20. Griffiths, B. S., I. E. Geoghegan, and W. M. Robertson. 2000. Testing genetically engineered potato, producing the lectins GNA and Con A, on non-target soil organisms and processes. J. Appl. Ecol. 37:159-170.[CrossRef]
21
21. Gyamfi, S., U. Pfeifer, M. Stierschneider, and A. Sessitsch. 2002. Effects of transgenic glufosinate-tolerant oilseed rape (Brassica napus) and the associated herbicide application on eubacterial and Pseudomonas communities in the rhizosphere. FEMS Microbiol. Ecol. 41:181-190.[Medline]
22
22. Hall, T. A. 1999. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp. Ser. 41:95-98.
23
23. Heuer, H., R. M. Kroppenstedt, J. Lottmann, G. Berg, and K. Smalla. 2002. Effects of T4 lysozyme release from transgenic potato roots on bacterial rhizosphere relative to communities are negligible natural factors. Appl. Environ. Microbiol. 68:1325-1335.[Abstract/Free Full Text]
24
24. Heuer, H., and K. Smalla. 1999. Bacterial phyllosphere communities of Solanum tuberosum L. and T4-lysozyme-producing transgenic variants. FEMS Microbiol. Ecol. 28:357-371.[CrossRef]
25
25. Huang, Y., A. M. Diner, and D. F. Karnosky. 1991. Agrobacterium rhizogenes-mediated genetic transformation and regeneration of a conifer: Larix decidua. In Vitro Cell. Dev. Biol. 27:201-207.
26
26. Huber, T., G. Faulkner, and P. Hugenholtz. 2004. Bellerophon; a program to detect chimeric sequences in multiple sequence alignments. Bioinformatics 20:2317-2319.[Abstract/Free Full Text]
27
27. Jouanin, L., M. Bonadebottino, C. Girard, G. Morrot, and M. Giband. 1998. Transgenic plants for insect resistance. Plant Sci. 131:1-11.
28
28. Kowalchuk, G. A., M. Bruinsma, and J. A. van Veen. 2003. Assessing responses of soil microorganisms to GM plants. Trends Ecol. Evol. 18:403-410.[CrossRef]
29
29. Lottmann, J., and G. Berg. 2001. Phenotypic and genotypic characterization of antagonistic bacteria associated with roots of transgenic and non-transgenic potato plants. Microb. Res. 156:75-82.[CrossRef]
30
30. Lottmann, J., H. Heuer, J. de Vries, A. Mahn, K. During, W. Wackernagel, K. Smalla, and G. Berg. 2000. Establishment of introduced antagonistic bacteria in the rhizosphere of transgenic potatoes and their effect on the bacterial community. FEMS Microb. Ecol. 33:41-49.[CrossRef][Medline]
31
31. Lottmann, J., H. Heuer, K. Smalla, and G. Berg. 1999. Influence of transgenic T4-lysozyme-producing potato plants on potentially beneficial plant-associated bacteria. FEMS Microb. Ecol. 29:365-377.[CrossRef]
32
32. Lukow, T., P. F. Dunfield, and W. Liesack. 2000. Use of the T-RFLP technique to assess spatial and temporal changes in the bacterial community structure within an agricultural soil planted with transgenic and non-transgenic potato plants. FEMS Microb. Ecol. 32:241-247.[CrossRef][Medline]
33
33. Lynch, J. M., A. Benedetti, H. Insam, M. P. Nuti, K. Smalla, V. Torsvik, and P. Nannipieri. 2004. Microbial diversity in soil: ecological theories, the contribution of molecular techniques and the impact of transgenic plants and transgenic microorganisms. Biol. Fertil. Soils 40:363-385.[CrossRef]
34
34. Maddaloni, M., F. Forlani, V. Balmas, G. Donini, L. Stasse, L. Corazza, and M. Motto. 1997. Tolerance to the fungal pathogen Rhizoctonia solani AG4 of transgenic tobacco expressing the maize ribosome-inactivating protein b-32. Transgen. Res. 6:393-402.[CrossRef]
35
35. Nubel, U., B. Engelen, A. Felske, J. Snaidr, A. Wieshuber, R. I. Amann, W. Ludwig, and H. Backhaus. 1996. Sequence heterogeneities of genes encoding 16S rRNAs in Paenibacillus polymyxa detected by temperature gradient gel electrophoresis. J. Bacteriol. 178:5636-5643.[Abstract/Free Full Text]
36
36. Palm, C. J., D. L. Schaller, K. K. Donegan, and R. J. Seidler. 1996. Persistence in soil of transgenic plant produced Bacillus thuringiensis var. kurstaki delta-endotoxin. Can. J. Microbiol. 42:1258-1262.
37
37. Pena, L., and A. Seguin. 2001. Recent advances in the genetic transformation of trees. Trends Biotechnol. 19:500-506.[CrossRef][Medline]
38
38. Saxena, D., and G. Stotzky. 2001. Bacillus thuringiensis (Bt) toxin released from root exudates and biomass of Bt corn has no apparent effect on earthworms, nematodes, protozoa, bacteria, and fungi in soil. Soil Biol. Biochem. 33:1225-1230.[CrossRef]
39
39. Schloss, P. D., B. R. Larget, and J. Handelsman. 2004. Integration of microbial ecology and statistics: a test to compare gene libraries. Appl. Environ. Microbiol. 70:5485-5492.[Abstract/Free Full Text]
40
40. Schmalenberger, A., and C. C. Tebbe. 2002. Bacterial community composition in the rhizosphere of a transgenic, herbicide-resistant maize (Zea mays) and comparison to its non-transgenic cultivar Bosphore. FEMS Microb. Ecol. 40:29-37.[Medline]
41
41. Séguin, A., D. Lachance, and P. J. Charest. 1996. Transient gene expression and stable genetic transformation into conifer tissues by microprojectile bombardment, p. 1-46. In K. Lindsey (ed.), Plant molecular biology manual, vol. B13. Kluwer Academic Publishers, Dordrecht, The Netherlands.
42
42. Siciliano, S. D., and J. J. Germida. 1999. Taxonomic diversity of bacteria associated with the roots of field-grown transgenic Brassica napus cv. Quest, compared to the non-transgenic B. napus cv. Excel and B. rapa cv. Parkland. FEMS Microbiol. Ecol. 29:263-272.[CrossRef]
43
43. Swofford, D. L. 2002. PAUP*: phylogenetic analysis using parsimony (*and other methods), 4th ed. Sinauer Associates, Sunderland, MA.
44
44. Tang, W., and R. J. Newton. 2003. Genetic transformation of conifers and its application in forest biotechnology. Plant Cell Rep. 22:1-15.[CrossRef][Medline]
45
45. Tapp, H., and G. Stotzky. 1998. Persistence of the insecticidal toxin from Bacillus thuringiensis subsp. kurstaki in soil. Soil Biol. Biochem. 30:471-476.[CrossRef]
46
46. Thompson, J. D., T. J. Gibson, F. Plewniak, F. Jeanmougin, and D. G. Higgins. 1997. The CLUSTAL-X Windows interface—flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 25:4876-4882.[Abstract/Free Full Text]
47
47. Wei-xiang, W., Y. Qing-fu, M. Hang, D. Xue-jun, and J. Wen-ming. 2004. Bt-transgenic rice straw affects the culturable microbiota and dehydrogenase and phosphatase activities in a flooded paddy soil. Soil Biol. Biochem. 36:289-295.[CrossRef]
48
48. Zwahlen, C., A. Hillbeck, P. Gugerli, and W. Nentwig. 2003. Degradation of the Cry1Ab protein within transgenic Bacillus thuringiensis corn tissue in the field. Mol. Ecol. 12:765-775.[CrossRef][Medline]

---

Applied and Environmental Microbiology, July 2007, p. 4128-4134, Vol. 73, No. 13
0099-2240/07/$08.00+0     doi:10.1128/AEM.02590-06
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

This Article Abstract Full Text (PDF) Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Reprints and Permissions Copyright Information Books from ASM Press MicrobeWorld Citing Articles Citing Articles via Google Scholar Google Scholar Articles by LeBlanc, P. M. Articles by Filion, M. Search for Related Content PubMed PubMed Citation Articles by LeBlanc, P. M. Articles by Filion, M. Agricola Articles by LeBlanc, P. M. Articles by Filion, M.