- Open Access
Field-grown miR156 transgenic switchgrass reproduction, yield, global gene expression analysis, and bioconfinement
© The Author(s) 2017
- Received: 10 July 2017
- Accepted: 19 October 2017
- Published: 30 November 2017
Genetic engineering has been effective in altering cell walls for biofuel production in the bioenergy crop, switchgrass (Panicum virgatum). However, regulatory issues arising from gene flow may prevent commercialization of engineered switchgrass in the eastern United States where the species is native. Depending on its expression level, microRNA156 (miR156) can reduce, delay, or eliminate flowering, which may serve to decrease transgene flow. In this unique field study of transgenic switchgrass that was permitted to flower, two low (T14 and T35) and two medium (T27 and T37) miR156-overexpressing ‘Alamo’ lines with the transgene under the control of the constitutive maize (Zea mays) ubiquitin 1 promoter, along with nontransgenic control plants, were grown in eastern Tennessee over two seasons.
miR156 expression was positively associated with decreased and delayed flowering in switchgrass. Line T27 did not flower during the 2-year study. Line T37 did flower, but not all plants produced panicles. Flowering was delayed in T37, resulting in 70.6% fewer flowers than controls during the second field year with commensurate decreased seed yield: 1205 seeds per plant vs. 18,539 produced by each control. These results are notable given that line T37 produced equivalent vegetative aboveground biomass to the controls. miR156 transcript abundance of field-grown plants was congruent with greenhouse results. The five miR156 SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) target genes had suppressed expression in one or more of the transgenic lines. Line T27, which had the highest miR156 overexpression, showed significant downregulation for all five SPL genes. On the contrary, line T35 had the lowest miR156 overexpression and had no significant change in any of the five SPL genes.
Because of the research field’s geographical features, this study was the first instance of any genetically engineered trait in switchgrass, in which experimental plants were allowed to flower in the field in the eastern U.S.; USDA-APHIS-BRS regulators allowed open flowering. We found that medium overexpression of miR156, e.g., line T37, resulted in delayed and reduced flowering accompanied by high biomass production. We propose that induced miR156 expression could be further developed as a transgenic switchgrass bioconfinement tool to enable eventual commercialization.
- Floral transition
- Gene flow
Switchgrass (Panicum virgatum L.) is a native North American perennial prairie grass mostly known for its use as a biofuel feedstock. The high biomass production, low input requirements, and its ability to be productive on marginal land are some features that make switchgrass an attractive cellulosic feedstock [1, 2]. However, the high degree of lignification of secondary cell walls (around 20% of switchgrass dry cell wall biomass) inhibits biomass conversion to fermentable sugars and biofuel in switchgrass, which, in turn, is an economic barrier to biofuel production [1–5]. Genetic engineering to reduce lignin levels in switchgrass cell walls appears to be essential for its optimal use as a biofuel crop [6–8]. Indeed, there are several success stories in producing transgenic switchgrass with altered lignification, which resulted in higher biofuel yield from field-grown biomass (e.g., [10, 11]), but the prospects of transgene flow from genetically engineered switchgrass is a regulatory concern. Transgene flow from switchgrass will likely need to be severely curtailed to facilitate the commercialization of transgenic varieties [6, 9]. This situation is especially pertinent in the eastern United States where switchgrass is endemic and common . Research has investigated several bioconfinement strategies, which include pollen ablation [13–15] and removal via site-specific recombinases [16, 17]. In addition, the delay or elimination of flowering itself could promote simultaneous improvements for a transgenic biomass crop such as switchgrass: it could decrease or eliminate pollen while simultaneously increase vegetative biomass [8, 18].
MicroRNAs (miRNAs) are an extensive class of small (20–24 nucleotides) regulatory RNAs that could be useful in genetic engineering to improve biofuel feedstocks by targeting stress responses, biomass production, and lignin content [19–31]. Specifically, miR156 targets the SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) transcription factor family, which is involved in the transition from vegetative to reproductive phases [32–35]. Overexpression of miR156 in switchgrass at low and moderate levels led to increased biomass and a non-flowering phenotype in the greenhouse . When two low and two moderate overexpressing lines were grown in the field, three of the lines flowered and one of these lines produced more biomass than the control . These results indicate that growth environment and gene expression play significant roles in the phenology of switchgrass. It should also be noted that the overexpression of miR156 at moderate to high levels led to an increase in saccharification efficiency and reduction in lignin content [36, 37].
Our research objectives in this study were to deploy a range of miR156-overexpressing switchgrass in a relevant field situation to closely examine flowering, reproduction, and biomass production. A field that would be considered a ‘marginal’ site (soils, fertility, and slope) on the Cumberland Plateau in Tennessee that is surrounded by forest enabled a 2-year study in which U.S. regulators allowed switchgrass plants to reproduce. In assessing a delayed/decreased flowering strategy for transgene bioconfinement of switchgrass, it was imperative to obtain two full flowering cycles in the field to gauge the degree of practical utility of this strategy. A transcriptomic study of the field-grown plants was performed to assess the influence of downstream genes impacted by miR156 expression, as well as any potential off-target effects, which are important for designing next-generation transgenic plants to further fine-tune the spatio-temporal expression of miR156 in switchgrass.
Field design and plant materials
Biomass and morphological characterization
Plants were checked weekly for the presence of panicles during both growing seasons, and first-date-to-flower was recorded. Aboveground biomass was harvested 10 cm above soil level after first frost (November) with plots pooled into a single harvest bag; the two ST2 plants from each plot were bagged separately from the surrounding plants per plot. All harvested biomass was oven-dried at 40 °C for 168 h, then dry biomass was tallied on a per-plot basis, and data were presented on a per-plant basis. Panicles were removed prior to harvest due to permit restrictions and bagged separately. Bags were stored in a greenhouse and allowed to air dry. Total panicle weights were recorded, averaged, and added to the average vegetative biomass weight to give total aboveground biomass production.
Panicles were counted during the removal process, and the lengths were measured for two randomly chosen panicles from each of five randomly selected plants per plot. A subsample of three panicles at the R4 stage of reproduction  was collected in September 2016 (year two) from each plot to tally flowers and spikelets per panicle.
The number of tillers per plant was tallied at each end-of-season harvest. Plant height (apex) was measured both before and after panicle removal. Leaf length, leaf width, stem diameter, and node number were taken at the end of the season on the two tallest tillers of each plant sampled. Leaf blade length and width were taken on the flag leaf or topmost mature leaf of each of the selected tillers. Tiller node number was counted from the soil line up, and representative internode diameter was taken using a Maxwell 150-mm digital caliper between the third and fourth nodes.
Seed collection and germination
After mature seeds were harvested from panicles, three subsamples per plant were tallied for 100-seed weight, then averaged. Seed number per plant was then derived by bulk seed weight and 100-seed weight. Seeds collected from transgenic lines or nontransgenic ‘Alamo’ controls were placed on solid MS basal medium , and germination percentage was calculated at 2 weeks after plating.
Cell wall characterization
End-of-season vegetative dry biomass was chipped to approximately 10-cm segments using a CS-4325 chipper shredder (Troy-Bilt, Valley City, Ohio) and then milled with a Wiley mill (Thomas Scientific, Model 4, Swedesboro, N.J.) through a 1-mm screen. Milled material was used to analyze the lignin content, syringyl-to-guaiacyl (S/G) monolignol ratio, and sugar release of the cell walls of each line by the National Renewable Energy Laboratory standard protocols. Lignin content and the S/G ratio were determined by pyrolysis molecular beam mass spectrometry as described in Sykes et al.  on an Extrel single-quadrupole molecular beam mass spectrometer. The peak intensities of lignin precursors were summed and used to estimate total lignin content. The S/G ratio was calculated by dividing the intensity of the syringyl peaks by the intensity of the guaiacyl peaks.
Sugar release was determined using the methods described in Selig et al. . Hydrolysis took place using the Ctec2 enzyme cocktail (Novozymes North America, Franklinton, NC). Released glucose levels were measured using the d-Glucose Assay Kit (glucose oxidase/peroxidase; GOPOD), and released xylose levels were determined by the d-Xylose Assay Kit (xylose dehydrogenase; XDH; Megasyme Intl., Bray, Ireland). Sugar release data were reported as grams of released sugar per gram of cell wall residue.
Microarray analysis was performed to determine downstream gene expression effects of miR156 overexpression. Three tillers were collected from each plot, resulting in four biological replicates for each of the four transgenic and ‘Alamo’ nontransgenic control lines. Total RNA was extracted from the combined tissues of randomly selected V3 stage tillers, as defined in Hardin et al. , from each line harvested on September 10, 2015 between 11:00 a.m. and 1:00 p.m. RNA was extracted using Tri-Reagent (Invitrogen, Carlsbad, Calif.) and subsequently cleaned and concentrated with the RNeasy® MinElute Cleanup Kit (Qiagen, Valencia, Calif.). Purified RNA (100 ng) was used for the expression analysis of each sample using a custom-designed switchgrass cDNA chip Pvi_cDNAa520831 (Affymetrix, Santa Clara, CA). Probe labeling, chip hybridization, and scanning were performed according to the manufacturer’s instructions for 3′ IVT PLUS Kit (Affymetrix). Data normalization among chips was conducted using the robust multichip average (RMA) . Gene selections based on Associative T test  were made using Matlab (MathWorks, Natick, MA). In this method, the background noise presented among replicates and technical noise during microarray experiments were measured by the residual presented among a group of genes whose residuals are homoscedastic. Genes whose residuals between the compared sample pairs that are significantly higher than the measured background noise level were considered to be differentially expressed. A selection threshold of 2 for transcript ratios and a Bonferroni-corrected P value threshold of 5.84201E−07 were used. The Bonferroni-corrected P value threshold was derived from 0.05/N in these analyses, where N is the number of probe sets on the chip. Microarray data are available in the ArrayExpress database accession number E-MTAB-5948 (http://www.ebi.ac.uk/arrayexpress).
Quantitative RT-PCR (qRT-PCR) analysis was used to assess transcript abundance of miR156 and its known target SPL genes. Total RNA was extracted using Tri-Reagent (Invitrogen) from V3 stage tillers collected mid-day on July 26, 2016. RNA samples were cleaned with the RNeasy® Mini Kit (Qiagen). The mature miR156 levels were determined using a highly sensitive stem-loop pulsed reverse transcription procedure  using a miR156-specific stem-loop primer. RT-PCR for SPL expression was performed using the High Capacity cDNA Reverse Transcription Kit (Applied Biosystems, Foster City, Calif.). SYBR Green (Applied Biosystems) was used as the reporter dye during qRT-PCR, and a QuantStudio™ 6 Flex Real-Time PCR System (Applied Biosystems) was used. The miR156 target gene transcript abundance qRT-PCR analysis included PvSPL1, PvSPL2, PvSPL3, and PvSPL6. miR156 expression was normalized using miR390 expression, and switchgrass PvUbq1 transcript abundance was used for normalization of data from each target gene with appropriate primers . Delta cycle threshold (ΔCt) was calculated by subtracting the target gene Ct from the Ct of the housekeeping gene (Housekeeping Ct—Target Ct = ΔCt).
SAS version 9.4 (SAS Institute Inc., Cary, NC) was used for all statistical analyses. A one-way ANOVA with Fisher’s least significant difference was used to compare means among lines within each year. Differences were considered significant when P values were less than or equal to 0.05.
miR156 overexpression levels affect flowering timing and reproductive effort
Flowering and reproduction of miR156-overexpressing switchgrass and the nontransgenic control in the field
Panicle number per plant
Panicle length (cm)
Spikelets per panicle
Percent seed germination
29.0 ± 1.6a
54.33 ± 1.69a
4.75 ± 3.47b
22.2 ± 1.1b
49.80 ± 1.29b
5.50 ± 1.89b
22.6 ± 1.7b
51.55 ± 1.67ab
22.75 ± 3.97a
0.0 ± 0.0d
9.9 ± 1.7c
26.77 ± 2.07c
0.25 ± 0.25b
103.5 ± 4.0a
73.34 ± 0.66a
27.5 ± 0.4a
34.75 ± 6.30a
60.6 ± 3.0c
61.46 ± 0.91c
24.6 ± 0.7c
15.25 ± 1.93b
98.8 ± 4.7a
68.01 ± 0.78b
25.8 ± 0.5bc
19.25 ± 3.33b
0.0 ± 0.0d
78.9 ± 7.5b
40.78 ± 1.26d
26.4 ± 1.0ab
18.25 ± 2.06b
All plants in the low overexpression lines flowered both years. T35 flowering phenology was delayed by 6 weeks relative to the control in year one, but was not delayed in year two (Fig. 2). T35 produced 22.1% fewer panicles, but were no different in length than the control (Table 1). The opposite was found in year two; T35 and the control produced the same number of panicles, but T35 panicles were shorter. However, T35 plants produced fewer flowers and seeds than the control for both years (Figs. 3, 4). Line T14 flowered at the same time as the control in year one and two weeks before the control in year two (Fig. 2). Although panicles emerged early in the season, they were fewer and smaller than control panicles (Table 1). T14 also produced fewer flowers and seeds than the control (Figs. 3, 4).
Seeds from the ‘Alamo’ non-transgenic control and transgenic lines were also collected and germinated. T35 was the only line to differ from the control in year one (18% higher germination), but there were no differences among transgenic lines in year 2, all of which had lower germination frequency than the control (Table 1).
Aboveground vegetative biomass production and plant morphology
Year one (2015) and year two (2016) end-of-season vegetative morphological data and cell wall characterization of miR156-overexpressing switchgrass and the wild-type control in the field
Leaf length (cm)
Leaf width (cm)
Internode diameter (mm)
Lignin (% CWR)
Sugar release (g/g CWR)
47.5 ± 3.7c
52.5 ± 1.1a
1.47 ± 0.02a
5.4 ± 0.1b
4.73 ± 0.08a
20.4 ± 0.4b
0.66 ± 0.01b
0.47 ± 0.00ab
30.9 ± 2.5c
36.6 ± 0.6c
1.15 ± 0.02c
4.9 ± 0.1c
4.25 ± 0.09b
21.3 ± 0.2a
0.69 ± 0.01a
0.44 ± 0.01b
41.2 ± 3.6c
48.2 ± 1.3b
1.34 ± 0.02b
5.1 ± 0.1bc
4.95 ± 0.11a
21.0 ± 0.1ab
0.69 ± 0.01a
0.49 ± 0.01a
193.2 ± 11.7a
15.3 ± 0.5e
0.34 ± 0.01e
5.5 ± 0.2b
1.27 ± 0.03d
20.5 ± 0.3b
0.59 ± 0.01c
0.44 ± 0.01b
106.0 ± 5.5b
27.3 ± 0.9d
0.73 ± 0.02d
7.8 ± 0.2a
3.04 ± 0.06c
20.5 ± 0.2b
0.59 ± 0.01c
0.49 ± 0.03a
112.2 ± 5.0b
52.6 ± 1.1a
1.17 ± 0.02b
8.0 ± 0.2c
5.36 ± 0.07b
23.2 ± 0.1b
0.66 ± 0.01ab
66.1 ± 3.7c
31.9 ± 1.1b
0.86 ± 0.03c
7.0 ± 0.2d
4.18 ± 0.13c
25.0 ± 0.2a
0.70 ± 0.02a
108.6 ± 6.8b
48.5 ± 0.8a
1.26 ± 0.03a
7.7 ± 0.1 cd
5.79 ± 0.07a
23.1 ± 0.1b
0.64 ± 0.03b
172.7 ± 16.0a
9.6 ± 0.4c
0.22 ± 0.01e
8.8 ± 0.3b
0.99 ± 0.03e
21.4 ± 0.4c
0.57 ± 0.00c
182.2 ± 6.6a
29.8 ± 1.1b
0.58 ± 0.02d
11.1 ± 0.2a
3.37 ± 0.06d
22.8 ± 0.2b
0.58 ± 0.01c
Cell wall composition (lignin content, digestibility, and sugar release) of the transgenic switchgrass lines had a few notable changes compared with the control. In both seasons, line T14 plant cells contained more lignin than the control (Table 2). T14, along with line T35 (both low overexpression lines), had higher S/G ratios than the control which suggests that they are more easily digestible (Table 2) . Both medium overexpression lines (T27 and T37) had lower S/G ratios than the control in both seasons. Transgenic lines did not differ from the control in sugar release (Table 2).
The expression levels of four SPL genes (PvSPL1, PvSPL2, PvSPL3, and PvSPL6) were also examined using quantitative RT-PCR to determine the effects of miR156 overexpression on its target genes in field-grown plants. All expression levels were examined on V3 stage vegetative tillers collected in year two (2016). The high variation among biological replicates resulted in no statistically significant differences for the expression levels of any SPL genes. Line T27 had the highest miR156 expression and showed the lowest PvSPL expression in general (Additional file 1: Figure S3).
The V3 stage tillers collected from the field in year one (2015) represent mid-season aboveground biomass for the global transcriptomic analysis (microarrays). Total RNA from all four transgenic lines and the ‘Alamo’ wild-type control was analyzed using Affymetrix microarray chips. Of the 85,587 probe sets examined, 14,507 were significantly up- or downregulated for one or more of the transgenic lines. Genes related to the miR156 pathway and flowering were chosen for further examination. Of 49 probe sets annotated as SPL according to known Arabidopsis thaliana and rice SPL sequences, eight SPL probes were found to be downregulated in open-flowering field conditions (Additional file 2: Table S1). SPL downregulation was negatively associated with mature miR156 overexpression (Fig. 6, Additional file 1: Figure S3). For the highest miR156 overexpression line T27, all eight SPL gene annotations were significantly downregulated (Additional file 2: Table S1). Six SPL gene annotations were downregulated in T37, which had the second highest miR156 overexpression. The expression of SPL genes appeared to have similar patterns to nontransgenics in the low overexpression lines; only two SPL genes were downregulated in T14, and none were downregulated in T35 (Additional file 2: Table S1). Probes corresponding to other important genes involved in flowering pathways, such as Arabidopsis AtFT (Flowering Locus T)/rice OsFTL (Flowering Locus T-Like) genes, were also significantly affected in miR156-overexpressing switchgrass (Additional file 2: Table S1).
Exploiting gene regulation by manipulating miRNAs could be useful in the sustained use of genetically engineered biofuel feedstocks to enhance desired traits such as abiotic and biotic stress responses, biomass yield, and lignin content [24–30]. miR156 targets the SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) transcription factor family, which is involved in many plant developmental processes including the vegetative to reproductive phase developmental transition [32–35, 47]. The overexpression of miR156 has been shown to delay flowering and increase biomass yield in multiple plant species [36, 48–51]. Arabidopsis thaliana plants engineered to overexpress miR156 had a moderate delay in flowering and an increase in total leaf number when grown under long days . A similar phenotype was seen in red clover (Trifolium pratense L.) engineered to overexpress miR156; transgenic red clover plants had an increased number of shoots and delayed flowering . Some transgenic events of switchgrass engineered to overexpress maize Corngrass1, a gene in the miR156 class of miRNAs, did not flower in a one-season California field trial, and weak overexpression levels did not affect biomass production . Transgenic switchgrass that overexpressed a rice miR156 precursor produced no flowering lines when grown in the greenhouse, and the low and medium overexpression lines produced more biomass than the control .
SPL downregulation causes delayed flowering in the field
Latitudinal origin and divergence of traits such as flowering time, growth and phenotype architecture, and disease susceptibility are used to classify switchgrass into either upland or lowland ecotypes [52–60]. Lowland switchgrass typically flowers later than varieties that originated in the north due to an elongated growth period . ‘Alamo,’ a lowland ecotype of switchgrass, typically flowers in mid–late June when grown in the southern United States . This study observed non-transgenic ‘Alamo’ switchgrass panicle production in mid- to late June for both growing seasons. Because the ‘Alamo’ nontransgenic control flowered in the same period as in the past studies [54, 60], a delayed flowering phenotype observed in transgenic lines can be contributed to miR156 overexpression rather than environmental effects. Transgenic lines T14, T35, and T37 flowered in the field. While this phenotype was different than the previous greenhouse study , the same was reported in a field study in Knoxville, Tenn. using the same miR156-overexpressing plants . Over the course of 3 years, T27 was the only line that did not produce panicles . SPL3 is an important upstream activator of floral meristem identity genes such as LEAFY, FRUITFULL, and APETALA1 , and the microarray revealed significant downregulation of SPL3 and APETALA1 in lines T27 and T37 (Additional file 2: Table S1). The medium overexpression lines were the only transgenic lines to have significant downregulation in SPL3, SPL4, and SPL5, which have overlapping functions to promote floral induction and transform the vegetative meristem to an inflorescence meristem [62, 63]. This downregulation of important SPL genes explains the delayed and non-flowering phenotypes of these two transgenic lines.
We observed that all transgenic lines produced shorter panicles than the control in year two, and lines T14 and T37 were also shorter in year one (Table 1). Overexpression of miR156 in rice resulted in short panicles with reduced spikelet and grain number . Line T37 was the only transgenic line to consistently produce fewer panicles and seeds than the control. While Xie et al.  found no difference in seed fertility, all miR156 switchgrass transgenic lines had lower seed germination than the control in year two (Table 1).
SPL downregulation results in altered plant phenotype
The trend in overexpression of miR156 in field-grown plants was consistent with that of previous greenhouse and field studies, as was the inverse relationship between miR156 and SPL gene target abundance (Fig. 6, Additional file 1: Figure S1) [36, 37]. Medium overexpression lines (T27 and T37) produced a high number of tillers, which is a common occurrence in plants overexpressing miR156 [36, 37, 48, 50, 51]. The high tiller number and short stature of T27 are most likely caused by a reduction in SPL1 and SPL2 expression (Additional file 1: Figure S4, Additional file 2: Table S1) which are important for tiller initiation and internode elongation . T27 and T37 tillers were thin compared to the control, and the leaves were smaller in both length and width for both lines (Table 2). When Arabidopsis thaliana was engineered to constitutively express miR156, plants produced leaves that resembled juvenile leaves in size, shape, and trichome production . miR156 promotes the expression of juvenile leaf traits by repressing SPL genes involved in plant maturation, such as SPL2/10 and SPL3/4/5, all of which were reduced in T27 and T37 (Additional file 2: Table S1) [34, 49, 62, 63, 65, 66]. The observed trends in vegetative biomass, height (without panicles), and tiller number were similar in ranking for year two data between this study and Baxter et al. , even though the latter study required panicle removal as a federal regulatory requirement in the field release permit. The high tiller number of line T37 without a reduction in height ‘rescued’ its biomass production. A miR156 overexpression level that falls between that of line T27 and line T37 (Fig. 6) would be ideal as it would most likely result in a non-flowering, high-yielding line. Such expression may be, ideally, triggered by environmental or developmental cues.
This two-year field study of miR156-overexpressing transgenic switchgrass is the first field experiment in the eastern U.S. in which USDA-APHIS-BRS regulators allowed open flowering. Thus, the present study was the first opportunity to closely examine the dynamics of switchgrass reproduction in the field using transgenic lines with a range of a miR156 expression. We found that medium overexpression levels of miR156 such as those in line T37 resulted in delayed and reduced flowering accompanied by high biomass production. Panicle size, seed production, and seed germination were also significantly reduced compared to the control. This outcome is the result of the downregulation of important miR156 SPL gene targets including SPL2/10 and SPL3/4/5. If miR156 overexpression was tied to developmental or environmental cues via conditional expression, then it could further optimize the use of miR156 overexpression as a bioconfinement tool.
CRJ drafted the manuscript, carried out the field study, collected phenotypic measurement data, and performed quantitative RT-PCR analysis. RJM assisted in experimental design, coordination of the study, and participated in drafting the manuscript. YT and JG performed the microarray analysis. RWS, GBT, and MFD completed the lignin and sugar release assays. YS contributed to the development of the study and writing of the grant proposal. ZYW provided the original plant clones and helped conceptualize the study. CNS conceived of the study and its design and coordination of the project, and participated in drafting the manuscript. All authors contributed to the text, data analysis, and interpretation of the results. All authors read and approved the final manuscript.
The authors would like to thank Kevin Hoyt, Martin Schubert, and the staff at the University of Tennessee Forest Resources and AgResearch and Education Center for their help in field maintenance. Thanks to Ben Wolfe, Marcus Laxton, Chance Losher, and Hayley Rideout for their assistance in data collection in the field. Thanks to Stacy Allen of Genomics Core at Noble Research Institute, LLC for microarray data generation. Thanks to Mitra Mazarei and Wusheng Liu for their indispensable advice and guidance.
The authors declare that they have no competing interests.
Availability of data and materials
The datasets used and analyzed during the current study are available from the corresponding author upon request. Microarray data will be deposited in the ArrayExpress database (http://www.ebi.ac.uk/arrayexpress).
Consent for publication
Ethics approval and consent to participate
Funding support was supplied by the Biotechnology Risk Assessment Grants Program (BRAG) Grant 2013-33522-20997 to CNS. The BRAG program is part of the United States Department of Agriculture—National Institute of Food and Agriculture. The transgenic plants were produced under the auspices of the BioEnergy Science Center (BESC), and some of the analyses were funded by BESC, which is a U.S. Department of Energy Bioenergy Research Center supported by the Office of Biological and Environmental Research in the DOE Office of Science. USDA Hatch funding also supported the research.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Sage RF, Peixoto M, Friesen P, Deen B. C4 bioenergy crops for cool climates, with special emphasis on perennial C4 grasses. J Exp Bot. 2015;66(14):4195–212.View ArticleGoogle Scholar
- McLaughlin SB, Kszos LA. Development of switchgrass (Panicum virgatum) as a bioenergy feedstock in the United States. Biomass Bioenergy. 2005;28(6):515–35.View ArticleGoogle Scholar
- Parrish DJ, Fike JH. The biology and agronomy of switchgrass for biofuels. Crit Rev Plant Sci. 2005;24(5–6):423–59.View ArticleGoogle Scholar
- Nigam PS, Singh A. Production of liquid biofuels from renewable resources. Prog Energy Combust. 2011;37(1):52–68.View ArticleGoogle Scholar
- Yang B, Wyman CE. Pretreatment: the key to unlocking low-cost cellulosic ethanol. Biofuel Bioprod Bioresour. 2008;2(1):26–40.View ArticleGoogle Scholar
- Stewart CN Jr. Biofuels and biocontainment. Nat Biotechnol. 2007;23(3):283–4.View ArticleGoogle Scholar
- Gressel J. Transgenics are imperative for biofuel crops. Plant Sci. 2008;174(1):246–63.View ArticleGoogle Scholar
- Jakob K, Zhou F, Paterson AH. Genetic improvement of C4 grasses as cellulosic biofuel feedstocks. In Vitro Cell Dev Plant. 2009;45:291–305.View ArticleGoogle Scholar
- Kausch AP, Hague J, Oliver M, Li Y, Daniell H, Mascia P, Watrud LS, Stewart CN Jr. Transgenic perennial biofuel feedstocks and strategies for bioconfinement. Biofuels. 2010;1(1):163–76.View ArticleGoogle Scholar
- Baxter HL, Mazarei M, Labbe N, Kline LM, Cheng Q, Windham MT, Mann DGJ, Fu C, Ziebell A, Sykes RW, Rodriguez M Jr, 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(7):914–24.View ArticleGoogle Scholar
- Baxter HL, Poovaiah CR, Yee KL, Mazarei M, Rodriguez M Jr, Thompson OA, Shen H, Turner GB, Decker SR, Sykes RW, Chen F. Field evaluation of transgenic switchgrass plants overexpressing PvMYB4 for reduced biomass recalcitrance. BioEnergy Res. 2015;8(3):910.View ArticleGoogle Scholar
- Kwit C, Stewart CN Jr. Gene flow matters in switchgrass (Panicum virgatum L.), a potential widespread biofuel feedstock. Ecol Appl. 2012;22(1):3–7.View ArticleGoogle Scholar
- Mariani C, Debeuckeleer M, Truettner J, Leemans J, Goldberg RB. Induction of male sterility in plants by a chimaeric ribonuclease gene. Nature. 1990;347:737–41.View ArticleGoogle Scholar
- Luo K, Duan H, Zhao D, Zheng X, Deng W, Chen Y, Stewart CN Jr, McAvoy R, Jian X, Wu Y, He A, Li Y. ‘GM-gene-deletor’: fused loxP-FRT recognition sequences dramatically improve the efficiency of FLP or CRE recombinase on transgene excision from pollen and seed of tobacco plants. Plant Biotechnol J. 2007;5(2):263–74.View ArticleGoogle Scholar
- 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. 2016;14(5):1281–90.View ArticleGoogle Scholar
- Moon HS, Abercrombie LL, Eda S, Blanvillain R, Thomson JG, Ow DW, Stewart CN Jr. Transgene excision in pollen using a codon optimized serine resolvase CinH-RS2 site-specific recombination system. Plant Mol Biol. 2011;75(6):621–31.View ArticleGoogle Scholar
- Somleva MN, Xuvv CA, Ryan KP, Thilmony R, Peoples O, Snell KD, Thomson J. Transgene autoexcision in switchgrass pollen mediated by the Bxb1 recombinase. BMC Biotechnol. 2014;14(1):79.View ArticleGoogle Scholar
- Kwit C, Moon HS, Warwick SI, Stewart CN Jr. Transgene introgression in crop relatives: molecular evidence and mitigation strategies. Trends Biotechnol. 2011;29(6):284–93.View ArticleGoogle Scholar
- Reinhart BJ, Weinstein EG, Rhoades MW, Bartel B, Bartel DP. MicroRNAs in plants. Gene Dev. 2000;16(13):1616–26.View ArticleGoogle Scholar
- Lagos-Quintana M, Raunhut R, Lendeckel W, Tuschl T. Identification of novel genes coding for small expressed RNAs. Science. 2001;294(5543):853–8.View ArticleGoogle Scholar
- Lau NC, Lim LP, Weinstein EG, Bartel DP. An abundant class of tiny RNAs with probable regulatory roles in Caenorhabditis elegans. Science. 2001;294(5543):858–62.View ArticleGoogle Scholar
- Lee RC, Ambros V. An extensive class of small RNAs in Caenorhabditis elegans. Science. 2001;294(5543):862–4.View ArticleGoogle Scholar
- Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116(2):281–97.View ArticleGoogle Scholar
- Zhang B, Pan X, Cobb GP, Anderson TA. Plant microRNA: a small regulatory molecule with big impact. Dev Biol. 2006;289(1):3–16.View ArticleGoogle Scholar
- Auer C, Frederick R. Crop improvement using small RNAs: applications and predictive ecological risk assessments. Trends Biotechnol. 2009;27(11):644–51.View ArticleGoogle Scholar
- Sun G. MicroRNAs and their diverse functions in plants. Plant Mol Biol. 2012;80(1):17–36.View ArticleGoogle Scholar
- Zhou M, Luo H. MicroRNA-mediated gene regulation: potential applications for plant genetic engineering. Plant Mol Biol. 2013;83(1):59–75.View ArticleGoogle Scholar
- Cui L-G, Shan J-X, Shi M, Gao J-P, Lin H-X. The miR156-SPL9-DFR pathway coordinates the relationship between development and abiotic stress tolerance in plants. Plant J. 2014;80(1):1108–17.View ArticleGoogle Scholar
- Trumbo JL, Zhang B, Stewart CN Jr. Manipulating microRNAs for improved biomass and biofuels from plant feedstocks. Plant Biotechnol J. 2015;13(3):337–54.View ArticleGoogle Scholar
- Zhang B, Wang Q. MicroRNA-based biotechnology for plant improvement. J Cell Physiol. 2015;230(1):1–5.View ArticleGoogle Scholar
- Djami-Tchatchou AT, Sanan-Mishra N, Ntushelo K, Dubery IA. Functional roles of microRNAs in agronomically important plants—potential as targets for crop improvement and protection. Front Plant Sci 2017;8:378.View ArticleGoogle Scholar
- Rhoades MW, Reinhart BJ, Lim LP, Burge CB, Bartel B, Bartel DP. Prediction of plant microRNA targets. Cell. 2002;110(4):513–20.View ArticleGoogle Scholar
- Poethig RS. Small RNAs and developmental timing in plants. Curr Opin Genet Dev. 2009;19(4):374–8.View ArticleGoogle Scholar
- Wu G, Park MY, Conway SR, Wang J-W, Weigel D, Poethig RS. The sequential action of miR156 and miR172 regulates developmental timing in Arabidopsis. Cell. 2009;138(4):750–9.View ArticleGoogle Scholar
- Matts J, Jagadeeswaran G, Roe BA, Sunkar R. Identification of microRNAs and their targets in switchgrass, a model biofuel plant species. J Plant Physiol. 2010;167(11):896–904.View ArticleGoogle Scholar
- Fu C, Sunkar R, Zhou C, Shen H, Zhang J-Y, Matts J, Wolf J, Mann DGJ, Stewart CN, Tang Y, Wang Z-Y. Overexpression of miR156 in switchgrass (Panicum virgatum L.) results in various morphological alterations and leads to improved biomass production. Plant Biotechnol J. 2012;10:443–52.View ArticleGoogle Scholar
- Baxter HL, Mazarei M, Dumitrache A, Rodriguez M, Natzke J, Gou J, Fu C, Sykes RW, Turner GB, Davis MF, Brown S, Davison B, Dixon RA, Wang Z-Y, Stewart CN Jr. Transgenic miR156 switchgrass in the field: growth, recalcitrance and rust susceptibility. Plant Biotechnol J. 2017. doi:10.1111/pbi.12747.Google Scholar
- Moore KJ, Moser LE, Vogel KP, Waller SS, Johnson BE, Pedersen JF. Describing and quantifying growth stages of perennial forage grasses. Agron J. 1991;83(6):1073–7.View ArticleGoogle Scholar
- Murashige T, Skoog F. A revised medium for rapid growth and bio assays with tobacco tissue cultures. Physiol Plant. 1962;15(3):473–97.View ArticleGoogle Scholar
- Sykes R, Yung M, Novaes E, Kirst M, Peter G, Davis M. High-throughput screening of plant cell-wall composition using pyrolysis molecular beam mass spectroscopy. In: Mielenz JR, editor. Biofuels: methods and protocols. New York City: Humana Press; 2009. p. 169–83.View ArticleGoogle Scholar
- Selig MJ, Tucker MP, Sykes RW, Reichel KL, Brunecky R, Himmel ME, Davis MF, Decker SR. Lignocellulose recalcitrance screening by integrated high-throughput hydrothermal pretreatment and enzymatic saccharification. Ind Biotechnol. 2010;6(2):104–11.View ArticleGoogle Scholar
- Hardin CF, Fu C, Hisano H, Xiao X, Shen H, Stewart CN Jr, Parrott W, Dixon RA, Wang Z-Y. Standardization of switchgrass sample collection for cell wall and biomass trait analysis. Bioenergy Res. 2013;6(2):755–62.View ArticleGoogle Scholar
- Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, Speed TP. Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003;4(2):249–64.View ArticleGoogle Scholar
- Dozmorov I, Centola M. An associative analysis of gene expression array data. Bioinformatics. 2003;19(2):204–11.View ArticleGoogle Scholar
- Varkonyi-Gasic E, Wu R, Wood M, Walton EF, Hellens RP. Protocol: a highly sensitive RT-PCR method for detection and quantification of microRNAs. Plant Methods. 2007;3(1):12.View ArticleGoogle Scholar
- Fu C, Mielenz JR, Xiao X, Ge Y, Hamilton CY, Rodriguez M Jr, Chen F, Foston M, Ragauskas A, Bouton J, Dixon RA, Wang Z-Y. Genetic manipulation of lignin reduces recalcitrance and improves ethanol production from switchgrass. Proc Natl Acad Sci USA. 2011;108(9):3803–8.View ArticleGoogle Scholar
- Preston JC, Hileman LC. Functional evolution in the plant SQUAMOSA-PROMOTER BINDING PROTEIN-LIKE (SPL) gene family. Front Plant Sci. 2013;4:80.Google Scholar
- Schwab R, Palatnik JF, Riester M, Schommer C, Schmid M, Weigel D. Specific effects of microRNAs on the plant transcriptome. Dev Cell. 2005;8(4):517–27.View ArticleGoogle Scholar
- Wu G, Poethig RS. Temporal regulation of shoot development in Arabidopsis thaliana by miR156 and its target SPL3. Development. 2006;133(18):3539–47.View ArticleGoogle Scholar
- Chuck GS, Tobias C, Sun L, Kraemer F, Li C, Dibble D, Arora R, Bragg JN, Vogel JP, Singh S, Simmons BA, Pauly M, Hake S. Overexpression of the maize Corngrass1 microRNA prevents flowering, improves digestibility, and increases starch content of switchgrass. Proc Natl Acad Sci USA. 2011;108(42):17550–5.View ArticleGoogle Scholar
- Zheng Q, Liu J, Goff BM, Dinkins RD, Zhu H. Genetic manipulation of miR156 for improvement of biomass production and forage quality in red clover. Crop Sci. 2016;56(3):1199–205.View ArticleGoogle Scholar
- McMillan C. Ecotypic differentiation within four North American prairie grasses. II. Behavioral variation within transplanted community fractions. Am J Bot. 1965;52(1):55–65.View ArticleGoogle Scholar
- Porter CL. An analysis of variation between upland and lowland switchgrass, Panicum virgatum L., in Central Oklahoma. Ecology. 1966;47(6):980–92.View ArticleGoogle Scholar
- Van Esbroeck GA, Hussey MA, Sanderson MA. Selection response and developmental bias for early and late panicle emergence in Alamo switchgrass. Crop Sci. 1998;38(2):342–6.View ArticleGoogle Scholar
- Casler MD. Ecotypic variation among switchgrass populations from the Northern USA. Crop Sci. 2005;45(1):388–98.View ArticleGoogle Scholar
- Casler MD, Vogel KP, Taliaferro CM, Wynia RL. Latitudinal adaptation of switchgrass populations. Crop Sci. 2004;44(1):293–303.View ArticleGoogle Scholar
- Casler MD, Vogel KP, Taliaferro CM, Ehlke NJ, Berdahl JD, Brummer EC, Kallenbach RL, West CP, Mitchell RB. Latitudinal and longitudinal adaptation of switchgrass populations. Crop Sci. 2007;47(1):2249–60.View ArticleGoogle Scholar
- Kiniry JR, Anderson LS, Johnson MV, Behrman KD, Brakie M, Burner D, Cordsiemon RL, Fay PA, Fritschi FB, Houx JH, Hawkes C, Juenger T, Kaiser J, Keitt TH, Lloyd-Reilley J, Maher S, Raper R, Scott A, Shadow A, West C, Wu Y, Zibilske L. Perennial biomass grasses and the Mason-Dixon Line: comparative productivity across latitudes in the Southern Great Plains. Bioenergy Res. 2013;6(1):276–91.View ArticleGoogle Scholar
- Milano ER, Lowry DB, Juenger TE. The genetic basis of upland/lowland ecotype divergence in switchgrass (Panicum virgatum). G3 Genes Genom Genet. 2016;6(11):3561–70.Google Scholar
- Grabowski PP, Evans J, Daum C, Deshpande S, Barry KW, Kennedy M, Ramstein G, Kaeppler SM, Buell CR, Jiang Y, Casler MD. Genome-wide associations with flowering time in switchgrass using exome-capture sequencing data. New Phytol. 2017;213(1):154–69.View ArticleGoogle Scholar
- Yamaguchi A, Wu M-F, Yang L, Wu G, Poethig RS, Wagner D. The microRNA-regulated SBP-box transcription factor SPL3 is a direct upstream activator of LEAFY, FRUITFULL, and APETALA1. Dev Cell. 2009;17(2):268–78.View ArticleGoogle Scholar
- Khan MRG, Ai X-Y, Zhang J-Z. Genetic regulation of flowering time in annual and perennial plants. WIREs RNA. 2014;5(3):347–59.View ArticleGoogle Scholar
- Xu M, Hu T, Zhao J, Park M-Y, Earley KW, Wu G, Yang L, Poethig RS. Developmental functions of miR156-regulated SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) genes in Arabidopsis thaliana. PLoS Genet. 2016;12(8):e1006263.View ArticleGoogle Scholar
- Xie K, Wu C, Xiong L. Genomic organization, differential expression, and interaction of SQUAMOSA promoter-binding-like transcription factors and microRNA156 in rice. Plant Physiol. 2006;142(1):280–93.View ArticleGoogle Scholar
- Wu Z, Cao Y, Yang R, Qi T, Hang Y, Lin H, Zhou G, Wang Z-Y, Fu C. Switchgrass SBP-box transcription factors PvSPL1 and 2 function redundantly to initiate side tillers and affect biomass yield of energy crop. Biotechnol Biofuels. 2016;9(1):101.View ArticleGoogle Scholar
- Shikata M, Koyama T, Mitsuda N, Ohme-Takagi M. Arabidopsis SBP-box genes SPL10, SPL11, and SPL2 control morphological change in association with shoot maturation in the reproductive phase. Plant Cell Physiol. 2009;50(12):2133–45.View ArticleGoogle Scholar