- Research article
- Open Access
Pollen-mediated gene flow from transgenic to non-transgenic switchgrass (Panicum virgatum L.) in the field
BMC Biotechnology volume 17, Article number: 40 (2017)
Switchgrass is C4 perennial grass species that is being developed as a cellulosic bioenergy feedstock. It is wind-pollinated and considered to be an obligate outcrosser. Genetic engineering has been used to alter cell walls for more facile bioprocessing and biofuel yield. Gene flow from transgenic cultivars would likely be of regulatory concern. In this study we investigated pollen-mediated gene flow from transgenic to nontransgenic switchgrass in a 3-year field experiment performed in Oliver Springs, Tennessee, U.S.A. using a modified Nelder wheel design. The planted area (0.6 ha) contained sexually compatible pollen source and pollen receptor switchgrass plants. One hundred clonal switchgrass ‘Alamo’ plants transgenic for an orange-fluorescent protein (OFP) and hygromycin resistance were used as the pollen source; whole plants, including pollen, were orange-fluorescent. To assess pollen movement, pollen traps were placed at 10 m intervals from the pollen-source plot in the four cardinal directions extending to 20 m, 30 m, 30 m, and 100 m to the north, south, west, and east, respectively. To assess pollination rates, nontransgenic ‘Alamo 2’ switchgrass clones were planted in pairs adjacent to pollen traps.
In the eastward direction there was a 98% decrease in OFP pollen grains from 10 to 100 m from the pollen-source plot (Poisson regression, F1,8 = 288.38, P < 0.0001). At the end of the second and third year, 1,820 F1 seeds were collected from pollen recipient-plots of which 962 (52.9%) germinated and analyzed for their transgenic status. Transgenic progeny production detected in each pollen-recipient plot decreased with increased distance from the edge of the transgenic plot (Poisson regression, F1,15 = 12.98, P < 0.003). The frequency of transgenic progeny detected in the eastward plots (the direction of the prevailing wind) ranged from 79.2% at 10 m to 9.3% at 100 m.
In these experiments we found transgenic pollen movement and hybridization rates to be inversely associated with distance. However, these data suggest pollen-mediated gene flow is likely to occur up to, at least, 100 m. This study gives baseline data useful to determine isolation distances and other management practices should transgenic switchgrass be grown commercially in relevant environments.
Switchgrass (Panicum virgatum L.) is a perennial C4 bunchgrass that is native to North America and is widely distributed from Canada to Mexico where it occasionally reaches “common” status in certain prairies [1,2,3,4], marshes , conservation reserve program settings [6, 7], and along roadsides and waste places [8, 9]. Elsewhere, it is typically not a large component of natural areas [10, 11] and is often found in much lower densities than those grown in agronomic settings. Traditionally switchgrass has been grown as a forage crop [12, 13], and because it is highly adaptable, it has served a variety of purposes including conservation, wildlife habitat, prairie restoration, erosion control, and in ornamental gardens [12, 14, 15]. Owing to its high biomass yield and suitability for marginal land, switchgrass has been the target of extensive development as a dedicated lignocellulosic biofuel feedstock . With the current drive to meet U.S. government-mandated benchmarks for energy production (16 billion gallons (1 gallon = 3.79 L) cellulosic biofuel by 2022; ), large-scale agronomic plantings of dedicated biofuel crops, including switchgrass, are projected to increase significantly through 2022 [14, 18, 19]. Although switchgrass research and development of conventional cultivars has significantly expanded over the last 25 years , extensive research efforts have been undertaken to produce transgenic switchgrass cultivars with increased biomass, enhanced saccharification efficiency, and modified lignin biosynthesis for improved biofuel properties [21,22,23,24,25,26,27,28,29,30]. Indeed, genetic engineering appears to be the only route to overcome biomass recalcitrance and sustainably improve biomass and biofuel production in tandem [31,32,33].
Although genetic engineering provides an opportunity to increase biomass yields and efficiencies in switchgrass production, transgenes from large-scale switchgrass plantings could escape to (1) non-transgenic switchgrass that might be used for livestock feed or other future uses, and (2) wild relatives in the landscape matrix [14, 31, 34]. Switchgrass is also being engineered for abiotic and biotic resistance traits (e.g. drought and salt tolerance, as well as pathogen resistance) , which could lead to increased distribution as well as shifted ecological niches if introgressed in conspecific or wild relative lineages. Switchgrass already exhibits a number of characteristics associated with invasiveness [36, 37]. These include C4 photosynthesis, a long canopy duration, few known pests and diseases, rapid growth in the spring (growing season), below-ground partitioning of nutrients in the fall (dormant season), and high water-use efficiency . The potential repercussions from habitat shifts include conservation concerns, with two different extremes being (1) the subsequent rapid spread of a transgene throughout the wild population and enhancing the invasiveness of the recipient population, or (2) extinction of the wild population [14,39,40,41,, 38–42]. The likely route of transgene escape from switchgrass is through pollen-mediated gene flow . It is well documented that wind-mediated pollen- and gene flow readily occurs in many grasses (reviewed by Kausch et al. ), and transgene escape and establishment has been documented with herbicide resistant creeping bentgrass (Agrostis stolonifera L.) [43,44,45]. Switchgrass shares similar characteristics (perennial, wind-pollinated, obligate outcrossing species) with creeping bentgrass, and considering field sites in the U.S. would be located near wild conspecifics and congeners, it is conceivable that pollen-mediated gene flow would occur . Yet to date relatively few studies have been published concerning switchgrass pollen biology and dispersal [46, 47]. The dearth of gene flow data, especially in the U.S. is largely to regulatory conditions: all regulated field trials to date require panicles to be removed at the R0 stage, obviating the ability to assess transgenic pollen movement and hybridization from transgenic pollen donors on non-transgenic pollen recipients.
It is imperative that we increase our existing knowledge on the drivers of transgene spread to satisfy imminent regulatory and conservation concerns regarding future transgenic bioenergy crops; indeed, stewardship of resources is a major reason. For switchgrass, we know little about pollen spread. To address the knowledge gap, we obtained release permits (11-007-106r and 13-046-104r-a1) through the United States Department of Agriculture’s Biotechnology Regulatory Services for the release of transgenic switchgrass expressing an orange fluorescent protein gene and a hygromycin resistance marker gene in order to perform a small multi-year pollen dispersal and cross-pollination field study. The study site was located in a secluded area surrounded by forest on all sides. Although the study was small in size, this experiment provided opportunity to gather much needed empirical data on transgenic switchgrass pollen dispersion and pollination distances under field conditions over a 3-year period in which permit conditions allowed plants to flower and set seed. These findings will aid in management practices for future cultivation of transgenic switchgrass, will be useful to both biosafety regulators and conservation biologists, and may also serve as a model for the incorporation of transgenic cultivars of other bioenergy feedstocks into the landscape.
Transgenic lines and strategy
We used ‘Alamo’ switchgrass, which is a viable lowland cultivar grown in the southeastern U.S. While switchgrass is known to have various ploidy levels, lowland varieties are diploidized ancient allotetraploids. Our strategy for estimating gene flow in switchgrass was based on using a relatively large number of pollen-source to create a windborne ‘pollen cloud’ for pollinate pollen recipient plants at various distances and directions. Given that switchgrass is putatively self-incompatible, our goal was for the vast majority of pollen produced by source plants to harbor the transgene. For that reason, we chose transgenic switchgrass lines containing multiple T-DNA insertions to ensure that the vast majority of pollen grains were transgenic. In addition, we also chose lines for which the orange fluorescent protein (OFP) marker gene was clearly detectable in pollen and vegetative tissues . Multiple transgenic events were screened for T-DNA copy number and OFP fluorescence  through Southern blot analysis and epifluorescent microscopy (Additional files 1 and 2). Based on the above selection criteria, transgenic lines ‘1–4’ and ‘2–9’ were chosen as pollen donors in the field experiment. Line 1–4’ possessed two transgene copies and line ‘2–9’ had six copies (Additional file 1), and exhibited OFP fluorescence when screened (Additional file 2). Moreover, we observed that practically every pollen grain was clearly OFP-positive, which ensured that the vast majority of pollen grains were transgenic. Assuming disomic inheritance and no association of different copies of the transgenes, 25% and 1.6% of non-transgenic pollen grains are expected in lines “1–4” and “2–9”, respectively (based on two and six transgene copies). Considering about 1/3 and 2/3 of the former and the latter plants in pollen donor population, the overall frequency of non-transgenic pollen should be about 9.4%. Since transgenic plants are expected to fertilize homozygous recessive (for the transgenic locus) pollen recipient plants, the same phenotypic frequencies are expected in hybrid seedlings. Of course, these are expectations on average, and it is possible that the apparent homogeneity of OFP pollen fluorescence observed in Additional file 2 is from sampling error, which is to say, we could have missed the non-transgenic pollen in the survey. Alternatively, the happenstance of transgene integration could have led to the observed pollen fluorescence distribution, which was apparently higher than the ~90% expected transgenic from the discussion above. We consider both scenarios in the Discussion section.
Plant growth and productivity
To fully assess the growth conditions, soil parameters such as water pH, buffer value and other micro nutrients were recorded, and the soil fertility was considered to be low-to-marginal (Additional file 3). The pollen-source plot mean tiller number (44 per plant) and height (107.5 cm) were comparable to plants in pollen-recipient plots (Table 1). The pollen-recipient plots had a wide range in means for both tiller number (7–80 tillers per plant) and plant height (74.0–208.6 cm). Above-ground biomass was measured for the 2014 season only (Table 1), and the mean weight of pollen-source plants (2.6 g per plant) was low compared to the pollen-recipient plots means (1.0–70.7 g per plant).
F1 progeny analysis
A total of 1,820 F1 seeds were collected from pollen-recipient plants and after cold-treatment, 962 germinated (52.9%). The total number of F1 seedlings analyzed for both the 2013 and 2014 seasons was 962, and a total of 204 F1 seedlings were deemed to be OFP positive then analyzed via OFP-gene specific PCR with congruent positive results. There were many seeds that did not germinate, which were not included in the final analysis as transgenic positive or negative. Nonetheless, most of these non-germinating seeds appeared to be small and pale with no detectable OFP-fluorescence in the seed; we deemed them to be non-viable. The number of F1 seedlings analyzed per pollen-receptor plant in the east direction was 717; 66 from the west, 91 from the south, and 88 from the north. Additional file 4, a sample gel from the PCR analysis, and Additional file 5, the summary Table S2, show that a portion of the progeny from pollen-recipient plants are not transgenic.
The east direction represented the long radius of the pollen-receptor plots as well as the direction of the prevailing wind. When normalized for pollen-receptor tiller number (Table 1), data indicated the number of F1 progeny collected from pollen-receptor plots to be inversely related to distance from the pollen-source plot (Fig. 1; Poisson regression, F1,15 = 12.98, P < 0.003). From the east direction, the majority of transgenic progeny were collected between 10 and 50 m plots (60.6%) compared to the remaining 40–100 m plots (39.4%). In addition, a greater than seven-fold decrease in transgenic progeny percentage was observed from the 10 m (79.2%) to 100 m (9.3%) plots. The percentages of transgenic F1 seedling observed in the west (10–30 m; 56.5–66.7%), south (10–30 m; 70.2–6.8%), and north (10–20 m; 32.0–34.2%) directions were similar to those observed in the east direction at similar distances (Additional file 5).
For all years, the Alamo 2 clonal type plants located in the pollen recipient plot began flowering in mid-July, while the ST1 clonal type plants in pollen-donor plots began booting the last week of July. Anthesis for both plant types overlapped starting in early August and continued until the last week of September of each growth season. During the pollen sampling period, the average wind speed was 5.2 km/h, the average dew point was 15.7 °C, the average temperature was 26.3 °C and the average humidity was 52.7% (Additional file 6) at the Oliver Springs study site (Additional file 7). OFP-tagged pollen counts collected from the eastward radius (direction of the prevailing wind) decreased as distance increased from the pollen-source plot (Fig. 2; Poisson regression, F1,8 = 288.38, P < 0.0001). These results indicate that the source pollen is capable of traveling at least 100 m, albeit at low frequency (ca. 0.01% of all switchgrass pollen from the source plants). More importantly, OFP-tagged pollen counts also declined with increased distance to the west (R 2 = 0.930), north (R 2 = 0.990), and south (R 2 = 0.993) directions from the source plants (Additional file 8). When OFP-tagged pollen counts (within 30 m) were compared across all the directions, the means of OFP-tagged pollen observed were highest in the eastward direction (104.3) as compared to the west (32.5), south (32.2) and north (29.4). The mean OFP-tagged pollen counts from 10 to 30 m away (104.3) from the source plants in the eastward direction were five-fold higher than mean pollen counts from 40 to 100 m (19.4).
The results presented here suggest that for switchgrass the majority of cross-pollination and transgenic pollen dispersal occur at short distances up to 50 m from the pollen source; however, a significant amount of pollen was detected at 50–100 m (Figs. 1 and 2). In addition, in the outermost pollen-receptor plot, 9.3% of the F1 seedlings recovered were transgenic. Our study is the first to record empirical data on switchgrass pollen dispersal and pollination distances in the field.
Transgenic cross-pollination percentages presented here in switchgrass as a function of distance from a transgenic pollen source are consistent with findings from other experimental studies concerned with transgene flow from a transgenic pollen source [49,50,51,52,53,54,55,56,57,58,59]. While the breeding systems and spatial scales involved in such studies vary widely, the primary conclusion is that the majority of transgenic cross-fertilization with nontransgenic counterparts occurs at short distances. In addition, studies addressing wind-pollinated grasses and having the largest spatial scales indicated transgenic cross-pollination percentages of 0.08% at 25 m from a transgenic source, and 0.006% at 250 m in maize , and levels declining to < 1% at 150 m in tall fescue . In our study, transgenic cross-pollination frequency declined with distance (beyond 50 m) from the transgenic source patch (Fig. 1), with significant reduction beginning at a distance 50 m from the transgenic source, falling to 9.3% at 100 m.
Our strategy was to purposefully use two clones that could not intraclonally self-pollinate, our experimental design was set up to gather worst-case scenario data: maximizing potential gene flow. A single clonal type of switchgrass in the pollen source plot and a different single clone as pollen receptor plants theoretically allow a singular mating choice given that switchgrass is putatively and obligate outcrosser. This situation is vastly different than that of a typical commercial field inhabited by a synthetic population of heterogeneous genotypes that would be expected to intermate with each other. Thus, our observed outcrossing rates by distance are likely many-fold higher than that expected in a commercial setting, even if all commercial plants are members of one cultivar. Nonetheless, declining transgenic pollen loads as a function of distance, and the possibility of selfing or apomixis, the former of which has been noted, rarely, in switchgrass [60, 61], may have led to nontransgenic seed production which could have influenced the frequency of observed transgenic progeny. Additionally, if the clone used for our receptor plants has any degree of self-compatibility, it is possible that, along with competitive interplays between the two pollen types and less-than-ideal overlap of floral phenologies, could have also influenced these results, which are similar to those exhibited by the wind-pollinated forage grass Festuca pratensis Huds. in a similarly designed, albeit nontransgenic, experiment . In that study, gene flow between the source and receptor plants of F. pratensis was 50% at a distance up to 15 m and showed a rapid decrease to 10% at 75 m, and 5% at 155 m .
Previously, we have not observed self-compatibility in our experimental clones in the greenhouse or field prior to these experiments. Nonetheless, we have also not isolated it from other genotypes the way they were isolated in these experiments. Thus, it is possible that selfing in these clones were practically below detection before. There have been prior reports of selfing in switchgrass [61, 63] dispelling the view that it is a completely obligate outcrosser. Those findings suggest there can be an occasional breakdown in the self-incompatibility system, which is the most likely explanation of the existence of some nontransgenic progeny. We estimated that at least ca. 90.6% of the pollen cloud was transgenic, with pollen microscopy evidence suggesting it may be even higher. Therefore, the 204 transgenic progeny (22%) from pollinations from pollen sink plants may be an underestimate of gene flow from the pollen-source population. Multiplying 204 by 0.094 (the expected non-transgenic pollen frequency from the transgenic plants) yields an estimated 19 additional progeny likely sired by transgenic plants. Therefore, of the 78% non-transgenic progeny from ‘pollen-sink’ plants, once we subtract the 19 progeny sired by transgenics, we are still left with 739 progeny (77%) coming from apparent self-pollination events. Thus, serendipitously, our results may have inadvertently led to increased evidence of selfing in switchgrass.
Our ability to empirically quantify pollen dispersal distances in switchgrass was arguably only feasible via the use of fluorescent protein tagged pollen source plants. In our case, the majority of pollen produced by our source plants expressed OFP, which distinguished switchgrass pollen from that of other heterospecifics at our field site and made fluorescent screening of pollen traps an efficient endeavor. Other methods utilized for the study of pollen movement have limitations, including destructive analysis or radioactive labels [44,65,, 62, 64–66]. Many researchers have utilized GFP-tagged pollen in transgenic plants to study gene flow [51, 67, 68]. In the present study, OFP was selected as the marker of choice since its expression can be visualized without sample destruction as has been shown with a different OFP variant in tobacco species . The OFP-tagged pollen method is more rapid and convenient than real-time PCR and other such DNA-based screening methods in which DNA extraction is a prerequisite. The OFP tracking technique has great potential as a tool for gene flow and risk assessment studies as well as for post-release monitoring of transgenic plants . The use of OFP and other fluorescent protein markers may also provide a viable option in estimating pollen flow of specific pollen donors in systems with concurrent heterospecific pollen flow and where palynological expertise to identify pollen grains to species is lacking.
Switchgrass pollen is diminutive with a mean diameter range of 36–48 μm (36.02 μm , 48.4 μm , and 43.7 μm ). Switchgrass pollen has short-lived viability after shedding, which also varies under environmental conditions. Greenhouse-grown ‘Alamo’ switchgrass produced pollen that was viable no longer than 20 min under sunny conditions and 150 min under cloudy conditions with a half-life of 4.3 min and 22.9 min; respectively (46). In another study, pollen from greenhouse grown switchgrass plants was exposed to ambient outdoor conditions and viability was completely lost after 60 min for most samples (half-life = 17.3 min) and 100 min for a single sample (47). Short viability periods require pollen to be dispersed quickly for pollination to occur over long distances, yet, depending on wind speeds and other environmental factors, a 60–100 min viability period may allow for pollen dispersal over a large area. Indeed, when switchgrass pollen dispersion was modeled based on 100 min viability, pollen was capable of moving up to 3.5 km in light wind and 6.5 km in stronger wind conditions; the majority of pollen is expected to travel < 200 m (47). In our study we observed a rapid decline in transgenic pollen deposition with increasing distance from the source plot; the majority of pollen deposition occurred within the first 50 m (Fig. 2, Additional file 8). The distance in which the majority of pollen was deposited was nearer to the pollen source than predicted by Ecker et al. ; however, the pattern of decline was similar. Our results may differ from the published model because of the secluded nature of our experimental site in which the forested border may have buffeted wind.
Nonetheless, these results are in agreement with pollen dispersal patterns observed in perennial ryegrass (Lolium perenne, ) where the bulk of pollen deposition occurred within 10 m from the source. Also, meadow fescue (Festuca pratensis ) pollen frequency declined rapidly within 75 m of the pollen source.
In a unique small field study in transgenic switchgrass, the first of its kind in which regulated switchgrass was allowed to flower in the field, we collected data on pollen movement and pollen-mediated transgenic hybridization. Over a 3-year period, the majority of pollen movement and cross-pollination occurred in close proximity to the pollen source, yet there was a significant amount of gene flow detected in our outermost sentinel plots (50–100 m) along with some apparent selfing. As U.S. regulatory agencies have set no allowable limits on the amount of adventitious presence for genetically engineered plants, the potentially acceptable limits of transgene flow is unknown in the U.S. for switchgrass and other crop species. If we used the E.U. upper limit of 0.09% [Regulation (EC) No 1830/2003; ], it would appear that in our worst-case scenario study, the gene flow observed here would have exceeded E.U. standards. Therefore, we could reasonably expect that at least a 100 m isolation distance would be required between transgenic edges and potential wild or cultivated nontransgenic plants that might receive transgenic pollen. For regulatory and risk assessment purposes, a spatially and genetically heterogeneous design could be useful to predict gene flow from commercial engineered switchgrass. Given the reproduction biology of this perennial species coupled with potentially fitness-enhancing traits (e.g. drought or salt tolerance, rapid growth rate, and disease resistance)  imbued via genetic engineering, transgene flow might require mitigation by an engineered bioconfinement system (e.g. selective male-sterility  or transgene excision .
Transgenic line selection and analysis
Switchgrass (‘Alamo’ ST1 clones) plants were genetically engineered through Agrobacterium-mediated transformation, as described in Burris et al. , with a vector containing the pporRFP (OFP) gene under the control of the maize ubiquitin (ZmUbi1) promoter and the hygromycin phosphotransferase (hph) selectable marker gene under the control of the rice actin (OsAct1) promoter . Transgenic T1 hemizygous plants were grown in a greenhouse (16/8 h day/night and 28/22 ° day/night). Transgenic lines were selectively chosen by evaluating T-DNA stable integration and transgene copy number using Southern blot analysis (Additional file 1) methodology as described in Moon et al. , with a few changes. Briefly, genomic DNA was extracted of T0 transgenic switchgrass events 1–4 and 2–9, along with other pporRFP transgenic switchgrass events. Ten micrograms of DNA from each sample was digested with BamHI, which cuts once in the T-DNA, to yield estimates of T-DNA insert number for each line. The blot was probed with 32P-labeled pporRFP gene fragment. To select the best OFP fluorescent events, visual analysis of OFP fluorescent levels was performed by epifluorescent microscopy for roots and sheaths as described by Burris et al.  and pollen (Additional file 2) epifluorescence microscopy [69, 75].
To examine transgenic pollen movement and pollination between transgenic and nontransgenic switchgrass plants, a 3-year field experiment was performed under USDA APHIS BRS release permits (11-007-106r and 13-046-104r). The field site in Oliver Springs, Tennessee U.S.A. was chosen because of its heavily forested borders that provided a natural barrier from wild or otherwise cultivated switchgrass plants. All experimental plants were transplanted by hand on September 20, 2011, and the total planted area was approximately 0.7 ha. A modified Nelder wheel design  consisted of a transgenic pollen-source plot and several pollen-recipient plots placed at 10 m intervals in the four cardinal directions from the pollen-source (Additional file 7). The pollen source and recipient plants had similar growth patterns and flowering times, which assured the highest possibility for transgene flow. The pollen-source plot contained 10 rows planted on 0.3 m centers, and each row contained 10 plants of the same type. Four rows were planted with clonally propagated plants from transgenic line ‘1–4’ (40 plants total) and the remaining six were planted clones from line ‘2–9’ (60 plants total). Post-transplant plant mortality led to the establishment of twenty-six ‘1–4’ plants and fifty ‘2–9’ plants, which remained throughout all experiments. Alamo 2 clones (cultivar ‘Alamo’) were chosen as the pollen-recipients to ensure cross-pollination would not occur between pollen-receptor plots. The long axis of the experiment was set up parallel with prevailing wind direction (eastward). Two receptor plants, 0.6 m apart from each other, were transplanted at distances in 10 m intervals. The northward plots were at 10 m and 20 m, southward and westward at 10 m to 30 m, and eastward from the source plants at10 m to 100 m (Additional file 7). Pollen-receptor plots suffered mortality as well. All plots had at least one plant and plots containing two receptor plants were: west 10 m, south 10 m, east 50 m, east 60 m, east 80 m, and east 90 m.
It is important to note that pollen-donor and pollen-recipient plants are both distinct clonal types originally derived from the cultivar ‘Alamo.’ These clonal types are sexually compatible with each other. However, clones of the same type are putatively not self-compatible. Specifically, pollen-donor plants are putatively not sexually compatible with other pollen-donor plants, and pollen-receptor plants are putatively not sexually compatible with pollen-receptor plants.
Plant growth and productivity
At the end of each growing season plant height, biomass production, and tiller number were recorded for both pollen-donor and –recipient plants (Table 1). Plant heights were represented by the height of the tallest tiller on each plant. Measurements were taken after panicles were removed for seed collection at the end of the growing season. Although seeds were not collected from pollen-source plants, for consistency panicles were removed prior to height measurement. Biomass measurements were taken as described in Baxter et al.  with a few modifications. For each plant, the senesced above-ground biomass was harvested 12.7 cm above the soil, oven-dried at 43 °C for approximately 96 h, and the dry biomass data were taken for each plant. After harvest, tiller numbers were counted for each plant in the field. To ensure plants were grown within typical fertility ranges, soil samples were collected from the source area as well several pollen-receptor plots (east at a distance of 25, 50, 70 and 100 m, and west, north and south at a distance of 15 m). A detailed soil analysis was carried out at the University of Tennessee Soil, Plant and Pest Center in Nashville, Tennessee. Soil properties such as water pH, buffer value, the amount of salts and trace elements were also analyzed (Additional file 3).
Pollen collection and screening
OFP pollen distribution measurements were taken by placing pollen traps adjacent to each pollen-receptor plot twice in September 2012 and three times in September 2013 for a total of five individual collections. Pollen traps were constructed by covering glass microscope slides (25 mm × 76 mm × 1 mm) with double-sided sticky adhesive tape (3 M136, 3 M Scotch Brand, St. Paul, Minn.) . Pollen traps were mounted horizontally onto vertical wooden stakes with a collection height of 1 m above the soil surface. For each sampling, the slides were exposed for six hour period to collect pollen. OFP pollen frequency was tallied using an epifluorescence microscope (Olympus Reflected Fluorescence system BX51, Olympus Corporation, Shinjuku, Tokyo, Japan) with a TxRed filter set (Texas Red®/Cy3.5, EX-D560/40,BS-595DCLP, EM-D630/60, Olympus Corporation, Shinjuku, Tokyo, Japan). OFP-tagged pollen was imaged using a color digital camera (Olympus QColor 5) with Qcapture imaging software (Q Imaging Corp., Burnaby, Canada). We used Poisson regression analysis (generalized linear mixed model using SAS 9.4; SAS Institute Inc., Cary, N.C.) to examine the relationship of OFP pollen counts and distance from the pollen-source plot. At the time of pollen collection, meteorological data such as temperature, wind speed, humidity, and dew point were recorded using an Ambient weather station located at the Cumberland Forest Unit Headquarters in Oliver Springs, Tenn. (Additional file 6).
Seed collection, germination and screening
At the end of the 2013 and 2014 growing seasons panicles were collected by hand from pollen-receptor plots and transported back to the laboratory; where seeds were cleaned, segregated by plot and stored at 4 °C for cold stratification. To estimate hybridization frequencies seeds were germinated in an environmental-controlled growth room (16/8 h day/night and 24/22 ° day/night). Seeds were placed in a Petri dishes lined with Whatman® 2 filter paper (Whatman International Ltd., Maidstone, England) and water was added as needed. Once germinated, the resultant seedlings were screened for the presence of OFP through epifluorescent microscopy as described by Burris et al. . One-to-two weeks post-germination seedlings were transferred to potting mix (Fafard® Peat-Lite Mix, Fafard, Inc., Anderson, S.C.) in an environmental controlled growth room (16/8 h day/night and 24/22 ° day/night). To further verify transgenic progeny seedlings were PCR-screened using primers (forward primer: TTTCAAAGCAAAGTGGGGTC; reverse primer: CACCATCTTCAAAGGTCATG) designed to amplify a 302 bp fragment of the pporRFP gene. Using a standard protocol genomic DNA was extracted from leaves  and PCR was performed using conditions described by Mann et al. . Poisson regression analysis was performed as previously described to evaluate the relationship between distance and transgenic progeny per pollen-recipient plant while taking tiller number per plant into account. For pollen-recipient plots containing two plants, the average number of transgenic progeny was calculated on a per-plant basis. The tallies of F1 germinated seedlings, transgenic F1 progeny, and tiller tallies were combined for both field seasons.
Weaver JE, Fitzpatrick TJ. The Prairie. Ecol Monogr. 1934;4:109–295.
Howe HF, Brown JS, Zorn‐Arnold B. A rodent plague on prairie diversity. Ecol Lett. 2002;5:30–6.
Baer SG, Collins SL, Blair JM, Knapp AK, Fiedler AK. Soil heterogeneity effects on tallgrass prairie community heterogeneity: an application of ecological theory to restoration ecology. Restor Ecol. 2005;13:413–24.
Haught JE, Myster RW. Effects of species, density, season and prairie-type on post-dispersal seed removal in Oklahoma. Am Midl Nat. 2008;159:482–8.
Ford MA, Grace JB. The interactive effects of fire and herbivory on a coastal marsh in Louisiana. Wetlands. 1998;18:1–8.
Mulkey VR, Owens VN, Lee DK. Management of switchgrass-dominated Conservation Reserve Program lands for biomass production in South Dakota. Crop Sci. 2006;46:712–20.
Adler PR, Sanderson MA, Weimer PJ, Vogel KP. Plant species composition and biofuel yields of conservation grasslands. Ecol Appl. 2009;19:2202–9.
Radford AE, Ahles HE, Bell CR. Manual of the vascular flora of the Carolinas. Chapel Hill: University of North Carolina Press; 1968.
Ahrens C, Ecker G, Auer C. The intersection of ecological risk assessment and plant communities: an analysis of Agrostis and Panicum species in the northeastern U.S. Plant Ecol. 2011;212:1629–42.
Grelen HE, Duvall VL. Common plants of longleaf pine-bluestem range. New Orleans: U.S. Department of Agriculture, Forest Service, Southern Forest Experiment Station; 1966. p. S0-23:96.
Byrd NA. Forestland Grazing: A Guide for Service Foresters in the South. Vol. 10. USDA, Forest Service, State and Private Forestry, Southeastern Area. 1980
Heaton E, Voigt T, Long SP. A quantitative review comparing the yields of two candidate C4 perennial biomass crops in relation to nitrogen, temperature and water. Biomass Bioenerg. 2004;27:21–30.
Parrish DJ, Fike JH. The biology and agronomy of switchgrass for biofuels. Crit Rev Plant Sci. 2005;24:423–59.
Kwit C, Stewart CN Jr. Gene flow matters in switchgrass (Panicum virgatum L.), a potential widespread biofuel feedstock. Ecol Appl. 2012;22:3–7.
Nageswara-Rao M, Soneji JR, Kwit C, Stewart CN Jr. Advances in biotechnology and genomics of switchgrass. Biotechnol Biofuels. 2013;6:77.
Sanderson MA, Reed RL, McLaughlin SB, Wullschleger SD, Conger BV, Parrish DJ, Wolf DD, Taliaferro C, Hopkins AA, Ocumpaugh WR, Hussey MA. Switchgrass as a sustainable bioenergy crop. Bioresource Technol. 1996;56:83–93.
US 110th Congress. Energy Independence and Security Act of 2007. Public Law 110–140; 2007. https://www.gpo.gov/fdsys/pkg/PLAW-110publ140/pdf/PLAW-110publ140.pdf. Accessed 28 Apr 2017.
Nageswara-Rao M, Hanson M, Agarwal S, Kwit C, Stewart CN Jr. Genetic diversity analysis of switchgrass (Panicum virgatum L.) populations using microsatellites and chloroplast sequences. Agroforest Syst. 2014;88:823–34.
Trainor AM, McDonald RI, Fargione J. Energy Sprawl is the largest driver of land use change in the United States. PLoS One. 2016;11, e0162269.
Casler MD, Vogel KP, Harrison M. Switchgrass germplasm resources. Crop Sci. 2015;55:2463–78.
Xu B, Escamilla-Treviño LL, Sathitsuksanoh N, Shen Z, Shen H, Zhang YHP, Dixon RA, Zhao B. Silencing of 4-coumarate:coenzyme A ligase in switchgrass leads to reduced lignin content and improved fermentable sugar yields for biofuel production. New Phytol. 2011;192:611–25.
Fu C, Xiao X, Xi Y, Ge Y, Chen F, Bouton J, Dixon RA, Wang ZY. Downregulation of cinnamyl alcohol dehydrogenase (CAD) leads to improved saccharification efficiency in switchgrass. BioEnerg Res. 2011;4:153–64.
Fu C, Mielenz JR, Xiao X, Ge Y, Hamilton CY, Rodriguez Jr M, Chen F, Foston M, Ragauskas A, Bouton J, Dixon RA, Wang ZY. Genetic manipulation of lignin reduces recalcitrance and improves ethanol production from switchgrass. Proc Natl Acad Sci USA. 2011;108:3803–8.
Saathoff AJ, Sarath G, Chow EK, Dien BS, Tobias CM. Down-regulation of cinnamyl-alcohol dehydrogenase in switchgrass by RNA silencing results in enhanced glucose release after cellulase treatment. PLoS One. 2011;6:e16416.
Shen H, Dixon RA, Parrott WA, Stewart CN Jr. Gateway-compatible vectors for high-throughput gene functional analysis in switchgrass (Panicum virgatum L.) and other monocot species. Plant Biotechnol J. 2012;10:226–36.
Baxter HL, Mazarei M, Labbe N, Kline LM, Cheng Q, Windham MT, Mann DGJ, Fu C, Ziebell A, Sykes RW, Rodriguez Jr M, Davis MF, Mielenz JR, Dixon RA, Wang Z-Y, Stewart CN Jr. Two-year field analysis of reduced recalcitrance transgenic switchgrass. Plant Biotechnol J. 2014;12:914–24.
Baxter HL, Poovaiah CR, Yee KL, Mazarei M, Rodriguez Jr M, Thompson OA, Shen H, Turner GB, Decker SR, Sykes RW, Chen F, Davis MF, Mielenz JR, Davison BH, Dixon RA, Stewart CN Jr. Field evaluation of transgenic switchgrass plants overexpressing PvMYB4 for reduced biomass recalcitrance. BioEnergy Res. 2015;8:910–21.
Poovaiah CR, Nageswara-Rao M, Soneji JR, Baxter H, Stewart CN Jr. Altered lignin biosynthesis using biotechnology to improve lignocellulosic biofuel feedstocks. Plant Biotechnol J. 2014;12:1163–73.
Wuddineh WA, Mazarei M, Zhang JY, Turner GB, Sykes RW, Decker SR, Davis MF, Udvardi MK, Stewart CN Jr. Identification and overexpression of a Knotted1-like transcription factor in switchgrass (Panicum virgatum L.) for lignocellulosic feedstock improvement. Front. Plant Sci. 2016;7:520.
Baxter HL, Alexander LW, Mazarei M, Haynes E, Turner GB, Sykes RW, Decker SR, Davis MF, Dixon RA, Wang ZY, Stewart CN Jr. Hybridization of downregulated-COMT transgenic switchgrass lines with field-selected switchgrass for improved biomass traits. Euphytica. 2016;209:341–55.
Stewart CN Jr. Biofuels and biocontainment. Nat Biotechnol. 2007;25:283–4.
Gressel J. Transgenics are imperative for biofuel crops. Plant Sci. 2008;174:246–63.
Jakob K, Zhou F, Paterson AH. Genetic improvement of C4 grasses as cellulosic biofuel feedstocks. In Vitro Cell Dev-Plant. 2009;45:291–305.
Kausch AP, Hague J, Oliver M, Li Y, Daniell H, Maschia P, Watrud LS, Stewart CN Jr. Transgenic biofuel feedstocks and strategies for bioconfinement. Biofuels. 2010;1:163–76.
Ridley CE, Mallory-Smith C. Gene flow and invasiveness in bioenergy systems. In: Quinn LD, Matlaga DP, Barney JN, editors. Bioenergy and Biological Invasions: Ecological, Agronomic and Policy Perspectives on Minimizing Risk. 2015;5:52-66.
Buddenhagen CE, Chimera C, Clifford P. Assessing biofuel crop invasiveness: a case study. PLoS ONE. 2009;4, e5261.
Chimera CG, Buddenhagen CE, Clifford PM. Biofuels: the risks and dangers of introducing invasive species. Biofuels. 2010;1:785–96.
Raghu S, Anderson RC, Daehler CC, Davis AS, Wiedenmann RN, Simberloff D, Mack RN. Adding biofuels to the invasive species fire? Science. 2006;313:1742.
Barney JN, DiTomaso JM. Nonnative species and bioenergy: are we cultivating the next invader? Bioscience. 2008;58:64–70.
Davis AS, Cousens RD, Hill J, Mack RN, Simberloff D, Raghu S. Screening bioenergy feedstock crops to mitigate invasion risk. Front Ecol Environ. 2010;8:533–9.
Barney JN. Bioenergy and invasive plants: quantifying and mitigating future risks. Invasive Plant Sci Manag. 2014;7:199–209.
Ecker G, Zalapa J, Auer C. Switchgrass (Panicum virgatum L.) genotypes differ between coastal sites and inland road corridors in the northeastern U.S. PLoS One. 2015;10:e0130414.
Belanger FC, Meagher TR, Day PR, Plumley K, Meyer WA. Interspecific hybridization between Agrostis stolonifera and related Agrostis species under field conditions. Crop Sci. 2003;43:240–6.
Watrud LS, Lee EH, Fairbrother A, Burdick C, Reichman JR, Bollman M, Storm M, King G, van De Water PK. Evidence for landscape-level, pollen-mediated gene flow from genetically modified creeping bentgrass with CP4 EPSPS as a marker. Proc Natl Acad Sci USA. 2004;101:14533–8.
Reichman JR, Watrud LS, Lee EH, Burdick CA, Bollman MA, Storm MJ, King GA, Mallory‐Smith CA. Establishment of transgenic herbicide‐resistant creeping bentgrass (Agrostis stolonifera L.) in nonagronomic habitats. Mol Ecol. 2006;15:4243–55.
Ge Y, Fu C, Bhandari H, Bouton J, Brummer EC, Wang Z. Pollen viability and longevity of switchgrass (Panicum virgatum L.). Crop Sci. 2011;51:2698–705.
Ecker G, Meyer T, Auer C. Pollen longevity and dispersion models for switchgrass. Crop Sci. 2013;53:1120–7.
Burris JN, Mann DGJ, Joyce BL, Stewart CN Jr. An improved tissue culture system for embryogenic callus production and plant regeneration in switchgrass. BioEnerg Res. 2009;2:267–74.
Ritala A, Nuutila AM, Aikasalo R, Kauppinen V, Tammisola J. Measuring gene flow in the cultivation of transgenic barley. Crop Sci. 2002;42:278–85.
Ma BL, Subedi KD, Reid LM. Extent of cross-fertilization in maize by pollen from neighboring transgenic hybrids. Crop Sci. 2004;44:1273–82.
Messeguer J, Marfà V, Català MM, Guiderdoni E, Melé E. A field study of pollen-mediated gene flow from Mediterranean GM rice to conventional rice and the red rice weed. Mol Breed. 2004;13:103–12.
Wang ZY, Lawrence R, Hopkins A, Bell J, Scott M. Pollen-mediated transgene flow in the wind-pollinated grass species tall fescue (Festuca arundinacea Schreb.). Mol Breed. 2004;14:47–60.
Darmency H, Klein EK, Gestat De Garanbé T, Gouyon P-H, Richard-Molard M, Muchembled C. Pollen dispersal in sugar beet production fields. Theor Appl Genet. 2009;118:1083–92.
Kim C-G, Park KW, Lee B, Kim DI, Park J-Y, Kim H-J, Park JE, An JH, Cho K-H, Jeong S-C, Choi KH, Harn CH, Kim HM. Gene flow from genetically modified to conventional chili pepper (Capsicum annuum L.). Plant Sci. 2009;176:406–12.
Kim C-G, Kim DI, Kim H-J, Park JE, Park KW, Jeong S-C, Choi KH, Kim HM. Assessment of gene flow from genetically modified anthracnose-resistant chili pepper (Capsicum annuum L.) to a conventional crop. J Plant Biol. 2009;52:251–8.
Van De Wiel CCM, Groeneveld RMW, Dolstra O, Kok EJ, Scholtens IMJ, Thissen JTNM, Smulders MJM, Lotz LAP. Pollen-mediated gene flow in maize tested for coexistence of GM and non-GM crops in the Netherlands: effect of isolation distances between fields. NJAS - Wageningen J Life Sci. 2009;56:405–23.
Rieben S, Kalinina O, Schmid B, Zeller SL. Gene flow in genetically modified wheat. PLoS One. 2011;6:e29730.
Scorza R, Kriss AB, Callahan AM, Webb K, Demuth M, Gottwald T. Spatial and temporal assessment of pollen- and seed-mediated gene flow from genetically engineered plum Prunus domestica. PLoS One. 2013;8, e75291.
Zhao XX, Tang T, Chen GM, Liu FX, Wang XL, Bu CP, Lu CM. Rationalizing the isolation distance needed for field trials involving genetically modified rapeseed (Brassica napus L.) in China. Chin Sci Bull. 2013;58:1558–67.
Martinez-Reyna JM, Vogel KP. Incompatibility systems in switchgrass. Crop Sci. 2002;42:1800–5.
Liu LL, Wu YQ. Identification of a selfing compatible genotype and mode of inheritance in switchgrass. Bioenerg Res. 2012;5:662–8.
Rognli OA, Nilsson N, Nurminiemi M. Effects of distance and pollen competition on gene flow in the wind-pollinated grass Festuca pratensis Huds. Heredity. 2000;85:550–60.
Casler MD, Tobias CM, Kaeppler SM, Buell CR, Wang Z-Y, Cao P, Schmutz J, Ronald P. The switchgrass genome: tools and strategies. Plant Genome. 2011;4:273–82.
Reinke DC, Bloom WL. Pollen dispersal in natural populations: A method for tracking individual pollen grains. Syst Bot. 1979;4:223–9.
Nilsson LA, Rabakonandrianina E, Pettersson B. Exact tracking of pollen transfer and mating in plants. Nature. 1992;360:666–8.
Song X, Munns K, Qiang S, Blackshaw R, Sharma R. Detection and quantification of 5 enolpyruvylshikimate-3-phosphate synthase (cp4 epsps) upon Brassica napus x Brassica juncea outcrossing using real-time PCR. Eur Food Res Technol. 2009;228:939–44.
Hudson LC, Chamberlain D, Stewart CN Jr. GFP-tagged pollen to monitor gene flow of transgenic plants. Mol Ecol Notes. 2001;1:321–4.
Moon HS, Halfhill MD, Hudson LC, Millwood RJ, Stewart CN Jr. Expression of green fluorescent protein in pollen of oilseed rape (Brassica napus L.) and its utility for assessing pollen movement in the field. Biotechnol J. 2006;1:1147–52.
Rice JH, Millwood RJ, Mundell RE, Chambers OD, Abercrombie LL, Davies HM, Stewart CN Jr. An orange fluorescent protein tagging system for real-time pollen tracking. BMC Res Notes. 2013;6:383–9.
Stewart CN Jr. Monitoring the presence and expression of transgenes in living plants. Trends Plant Sci. 2005;10:390–6.
Bragg LH. Pollen size variation in selected grass taxa. Ecology. 1969;50:124–7.
Giddings GD, Sackville Hamilton NR, Hayward MD. The release of genetically modified grasses. Part 1: pollen dispersal to traps in Lolium perenne. Theor Appl Genet. 1997;94:1000–6.
Moon HS, Abercrombie LL, Eda S, Blanvillain R, Thomson JG, Ow DW, Stewart CNJR. Transgene excision in pollen using a codon optimized serine resolvase CinH-RS2 site-specific recombinase system. Plant Mol Biol. 2011;75:621–31.
Millwood RJ, Moon HS, Poovaiah CR, Muthukumar B, Rice JH, Abercrombie JM, Abercrombie LL, Green WD, Stewart CN Jr. Engineered selective plant male sterility through pollen‐specific expression of the EcoRI restriction endonuclease. Plant Biotechnol J. 2015;14:1281–90.
Moon HS, Eda S, Saxton AM, Ow DW, Stewart CN Jr. An efficient and rapid transgenic pollen screening and detection method using flow cytometry. Biotechnol J. 2011;6:118–23.
Nelder JA. New kinds of systematic designs for spacing experiments. Biometrics. 1962;18:283–307.
Fonseca AE, Westgate ME, Doyle RT. Application of fluorescence microscopy and image analysis for quantifying dynamics of maize pollen shed. Crop Sci. 2002;42:2201–6.
Stewart CN Jr, Via LE. A rapid CTAB DNA isolation technique useful for RAPD fingerprinting and other PCR applications. Biotechniques. 1993;14:748–50.
Mann DGJ, LaFayette PR, Abercombie LL, King ZR, Mazarei M, Halter MC, Poovaiah CR, Baxter H, Shen H, Dixon RA, Parrott WA, Stewart CN Jr. Gateway-compatible vectors for high-throughput gene functional analysis in switchgrass (Panicum virgatum L.) and other monocot species. Plant Biotechnol J. 2012;10:226–36.
We thank numerous people who facilitated and assisted with field work and lab work, including Ben Wolfe, Marcus Laxton, Morgan Peck, John Alexander, and Jordan Skeen. We also thank Kevin Hoyt, Martin Schubert, and the staff at the University of Tennessee Institute of Agriculture, Forestry Resources AgResearch and Education Center for help with field site acquisition, preparation, and maintenance. Thanks also to Zeng-Yu Wang at the Samuel Roberts Noble Foundation, Ardmore, Oklahoma, USA for the gift of the ‘Alamo’ ST1 clone. This project was supported by Biotechnology Risk Assessment Grant Program competitive grant no. 2010-39211-21699 from the USDA National Institute of Food and Agriculture (NIFA). C. Neal Stewart Jr. also received support from the BioEnergy Science Center, a Bioenergy Research Center, supported by the Office of Biological and Environmental Research in the US Department of Energy Office of Science and the University of Tennessee. The experiment was performed under release-into-the-environment permits issued by United States Department of Agriculture Animal and Plant Health Inspection Service Biotechnology Regulatory Services and all regulations were followed.
This project was supported by Biotechnology Risk Assessment Grant (BRAG) Program competitive grant no. 2010-39211-21699 from the USDA National Institute of Food and Agriculture (USDA-NIFA) and the Agricultural Research Service (ARS), and USDA Hatch grants.
Availability of data and materials
The datasets supporting the conclusions of this article are included within the article and its additional files.
RM collected data for the 2014 field season and drafted the manuscript. RY collected data for the 2013 field season. M N-R planted the field and collected OFP pollen data and other data. E T-E, CRJ, and M H participated in seed collection, progeny analysis, and microscopy. JB produced the transgenic switchgrass and performed the Southern blot analysis. CNS, Jr. and CK conceptualized, designed, and supervised the study. All authors participated in writing and editing, and approved the final manuscript.
CNS is the plant biotechnology section editor for BMC Biotechnology. He played no editorial role in this paper.
Consent for publication
Ethics approval and consent to participate
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Southern blot analysis of T0 transgenic switchgrass events containing stable integration of the T-DNA. Blot was probed with 32P-labeled pporRFP gene fragment. Lanes labeled with ‘2BR’ contain transgenic line genomic DNA samples, lane ‘ST1 control” is the nontransgenic control genomic DNA sample, lane ‘2BR plasmid’ contains the plasmid positive control. Genomic DNA was digested BamHI genomic DNA (10 μg) from each plant sample. (PDF 46 kb)
Sample of orange fluorescent protein (OFP)-tagged switchgrass pollen from source plants. A), white light, B), TxRed filter, C), merged image. The scale bar is 400 micrometers. (PDF 61 kb)
Soil analysis from field site collected from the pollen-source plot and seven pollen-recipient plots. (PDF 6 kb)
A sample gel from PCR analysis for the pporRFP OFP gene that was performed on F1 switchgrass progeny (lanes 1–12) collected from field study pollen-recipient plots. A water only sample (W) and a nontransgenic pollen-recipient parental type (A-; Alamo II) were used a negative controls. A transgenic parental type sample (A+; line 2–9) and a plasmid with the pporRFP gene (P+) were used as positive control. A DNA size marker (M) was used to confirm expected band size (302 bp), and blank lanes (B) were used to separate sample types. (PDF 45 kb)
PCR progeny screen of F1 seedlings collected from pollen-recipient plots in the field. (PDF 6 kb)
Mean meteorological data collected for the growing seasons of 2012 and 2013 at Oliver Springs, Tennessee. (PDF 6 kb)
Experimental layout of field study located in Oliver Springs, Tennessee, U.S.A. Transgenic switchgrass (cultivar ‘Alamo’, clone ST1) served as a pollen-source (10 rows × 10 plants) while non-transgenic switchgrass (cultivar ‘Alamo’, clone Alamo 2) plants were placed in pollen-recipient plots planted at 10 m intervals. (PDF 68 kb)
Average number of orange fluorescent protein (OFP)-tagged switchgrass pollen grains detected as a function of distance in the north, south, and east directions (East, R 2 = 0.953; North, R 2 = 0.990; West, R 2 = 0.930; and South R 2 = 0.993) Pollen samples were collected in the field during the 2012 and 2013 growing seasons. (PDF 54 kb)