Original Research ARTICLE
Using Transcriptome Analysis to Identify Genes Involved in Switchgrass Flower Reversion
- 1College of Agronomy, Northwest A&F University, Yangling, China
- 2State Key Laboratory of Crop Stress Biology for Arid Areas, Yangling, China
- 3Yangling Vocational & Technical College, Yangling, China
- 4Institute of Cotton Research, Chinese Academy of Agricultural Sciences, Anyang, China
Floral reversion is a process in which differentiated floral organs revert back to vegetative organs. Although this phenomenon has been described for decades, the underlying molecular mechanisms remain unclear. In this study, we found that immature switchgrass (Panicum virgatum) inflorescences can revert to neonatal shoots when incubated on a basal medium with benzylaminopurine. We used anatomical and histological methods to verify that these shoots were formed from floret primordia through flower reversion. To further explore the gene regulation of floral reversion in switchgrass, the transcriptome of reversed, unreversed, and uncultured immature inflorescences were analyzed and 517 genes were identified as participating in flower reversion. Annotation using non-redundant databases revealed that these genes are involved in plant hormone biosynthesis and signal transduction, starch and sucrose metabolism, DNA replication and modification, and other processes crucial for switchgrass flower reversion. When four of the genes were overexpressed in Arabidopsis thaliana, vegetative growth was facilitated and reproductive growth was inhibited in transgenic plants. This study provides a basic understanding of genes regulating the floral transition in switchgrass and will promote the research of floral reversion and flower maintenance.
Flowering is an essential process in the angiosperm life cycle and facilitates the transmission of hereditary components to offspring via sexual reproduction (Scutt and Vandenbussche, 2014). Mature plants generally produce determinate blossoms and seeds (Irish, 2010). However, in some plant species, the differentiated floral meristems can also be reversed into a vegetative phase when the environment becomes adverse for flowering (Battey and Lyndon, 1990; Tooke et al., 2005; Wang et al., 2010; Sen et al., 2013). This ‘reverse development’ is known as flower reversion, and has been described in both dicots and monocots (Battey and Lyndon, 1990; Tooke et al., 2005; Trevaskis et al., 2007). Studying flower reversion has greatly contributed to our understanding of the molecular mechanisms involved in flower formation and flower maintenance (Tooke et al., 2005; Irish, 2010; O’Maoileidigh et al., 2014).
The first detailed description of floral reversion was published in 1990 (Battey and Lyndon, 1990). Since then, researchers have studied the underlying reversion mechanisms using mutants with distorted flower development (Scutt and Vandenbussche, 2014). In Arabidopsis thaliana, the LEAFY (LFY) gene is required for flower determination, and studies have further indicated it plays a significant role in suppressing flower reversion. Knocking out LFY and Agamous (AG) function in A. thaliana leads to floral meristem reversion under a short photoperiod (Okamuro et al., 1996). Studies have also indicated that LFY and APETALA 1 (AP1) promote flower formation in A. thaliana by inhibiting the expression of AGAMOUS-LIKE 24 (AGL24), a flower formation suppressor (Yu et al., 2004). The MADS-box genes are another group of important regulators in flower development (Fornara et al., 2003). Downregulated expression of the MADS-box gene TM29 in tomato (Solanum lycopersicum L.) caused parthenocarpic fruit development and floral reversion (Ampomah-Dwamena et al., 2002). In maize (Zea mays L.), the indeterminate 1 (ID1) gene was found to be a floral reversion suppressor, since the non-functional mutant id1-m1 produced plantlets instead of florets on the tassel branches (Colasanti et al., 1998). These single mutant studies have contributed to our knowledge of flower reversion, yet a comprehensive understanding of gene expression patterns during flower reversion is still missing.
RNA sequencing (RNA-seq) is a powerful technique to examine transcriptome changes under specified spatial-temporal conditions and is especially useful for investigating insufficiently studied phenomena or species. In plants, RNA-seq has identified many important genes involved in metabolite biosynthesis (Wang et al., 2012; Yang et al., 2014), biotic and abiotic stress resistance (Liu et al., 2011), tillering (Palmer et al., 2012), flower development (Singh and Jain, 2014; Zhang et al., 2014), and fruit formation (Jiang et al., 2015).
Switchgrass (Panicum virgatum) is considered a model bioenergy crop due to its high biomass production and broad adaptability (Parrish et al., 2012). We found that switchgrass INFs underwent flower reversion when cultured on SIM (Chai et al., 2012), but produced flowers in vitro when cultured on Murashige and Skoog’s basal medium (MS; Murashige and Skoog, 1962). To identify the genes involved in switchgrass flower reversion, we compared the gene expression patterns of INFs incubated on SIM and MS media and uncultured INFs. Our primary aim was to gain an overall understanding of genes regulating flower reversion, which could benefit switchgrass breeding for high yields of biomass in the future.
Evidence of Flower Reversion
The SIM culture induced neonatal shoots in switchgrass INFs (Figure 1A). As shown in Figure 1A, the emerged shoots were connected with unreversed florets, and seem to be derived from the same rachilla, yet the base was clearly separated. There were shortened pedicel-like organs at the bottom of shoots (Figures 1B,C). Moreover, two small leaves were bent inward at the base of the shoot (Figure 1B) and were similar to glumes of the floret (Figure 1C).
FIGURE 1. Evidence of flower reversion in switchgrass. (A) Reversed shoots and unreversed florets. (B) Single reversed shoot separated. (C) A normal floret from switchgrass plant. (D) Longitudinal section of a reversed shoot. (E) Longitudinal section of a half-reversed floret. (F) Longitudinal section of an unreversed floret. Pd, pedicel; pdl, pedicel-like; gu, glume; gul, glume-like; pi, pistil; pil, pistil-like; sml, stamen-like; an, anther; st, style; ov, ovule. (A–C) Bar = 2 mm; (D,E) bar = 0.5 mm.
Next, we investigated the histological structure of the reversed florets (Figure 1D). In longitudinal sections of normal florets, the stamens and pistils were fully developed (Figure 1F) with separated styles and the ovule clearly visible (Figure 1F). In the half-reversed florets, the pistils were under-developed and spherical-shaped (Figure 1E); the stamens were normal (Figure 1E). In the reversed florets, the pistils developed into a deformed tissue without a style or stigma, and the stamens degenerated into pin-like tissues (Figure 1D). These results indicate that the emerged shoots on SIM were generated from floret primordia on the INF, and were most likely initiated from the pistil primordia.
Flower Reversion of Switchgrass in Different Stages
We harvested the 1 cm long INF at elongation stage 4 (E4; Figures 2A,B) (Moore et al., 1991; Burris et al., 2009). After 20 days of incubation on SIM medium, floret primordia were further developed, from which the compact, glume-like organs were formed (Figure 2C). In contrast, explants incubated on the MS medium developed stretched peduncles and young florets (Figure 2D). After 30 days, the glume-like organs on SIM further elongated and transformed into leaf-like tissues (Figure 2E). On MS medium, the peduncles elongated and the florets were initiated (Figure 2F). After 60 days of incubation, INFs were completely reversed into neonatal shoots on SIM medium (Figure 2G). Florets on MS medium developed into mature flowers with anthers and stigmas (Figure 2H).
FIGURE 2. Incubation of the experimental material. (A) Switchgrass plants used in this research. (B) Immature inflorescences used for incubation. (C) Explant cultured on SIM medium for 20 days. (D) Explant cultured on MS medium for 20 days. (E) Explant cultured on SIM medium for 30 days. (F) Explant cultured on MS medium for 30 days. (G) Explant cultured on SIM medium for 60 days. (H) Explant cultured on MS medium for 60 days. Bars = 1 cm.
RNA Sequencing and Identification of Novel Transcripts
To investigate the molecular mechanism of flower reversion, we performed RNA-seq on reversing flowers on SIM medium (simplified as REV), in vitro flowers on MS medium (simplified as FLO), and uncultured immature inflorescences (simplified as INF). Considering that gene regulation of flower reversion occurred prior to the phenotypic alteration, we used explants incubated for 20 days in the subsequent transcriptome analysis.
Sequencing the switchgrass samples (n = 12) generated a total of 62.24 Gb of sequencing data, comprising 308.14 million clean reads. To guarantee accurate data analysis, we filtered out all low-quality reads and obtained 225.29 million high-quality reads from 45.51 Gb of data (Table 1). The high-quality reads were mapped to the Switchgrass reference genome1 (version 4.1) and assembled separately. The transcripts were merged to a uniform set of transcripts for all 12 samples. As a result, 63,915 loci containing 93,800 consensus transcripts were detected. Compared to the reference annotation, 33,499 loci were missed while 13,491 were newly detected (Supplementary Table S1). Statistical analysis showed that the percentage of missed loci ranged from 32.7 to 36.0%, and the percentage of novel loci ranged from 9.4 to 12.3% (Supplementary Table S1).
Mining of Differentially Expressed Genes
Pairwise comparison between REV and FLO was expected to show differences in gene expression DEGs related to flower formation and reversion. Comparison of INF and REV was expected to generate DEGs related to floral reversion and the in vitro culture environment response, and comparison between INF and FLO was expected to generate DEGs related to in vitro flower development and in vitro culture environment response. To explore genes involved in switchgrass flower reversion, we combined results of the pairwise comparisons and used DEGs from REV vs. FLO and REV vs. INF, but not FLO vs. INF, as candidate genes (Figure 3).
Among all pairwise comparisons, REV vs. FLO produced the smallest number of DEGs (994), followed by FLO versus INF (7,100) and REV versus INF (10,388). The large difference in the number of DEGs of the three comparisons indicated that REV and FLO were the most similar materials, and in vitro culture resulted in a more dramatic change in expression patterns than flower reversion. Through the filtration approach described above, 517 genes were identified as significant in the induction of floral reversion. These genes are hereby referred to as ‘filtered-DEGs’ (Figure 3 and Supplementary Table S2).
Gene expression pattern analysis showed that 496 of the 517 filtered-DEGs were upregulated both in REV vs. FLO and REV vs. INF, suggesting the filtered-DEGs likely promoted the induction of flower reversion (detailed gene lists are provided in Supplementary Table S2).
Expression Pattern Verification Using qPCR
To verify the efficacy of filtered-DEGs, we performed qPCR with 12 of the high-expression filtered-DEGs (Supplementary Presentation S2). The results demonstrated that all examined genes showed expression trends similar to the RNA-seq results, with variation in the relative expression levels of some genes (Table 2). Upregulation of the genes 3NG053500, 5NG047700, and J517000 was more prominent in qPCR, while with the genes 1KG404800 and 1KG545400, it was less prominent. The difference might be caused by different qPCR and RNA-seq measurements of expression levels (Wagner et al., 2012; Biassoni and Raso, 2014) and differences in amplification efficiency of the primers.
Gene Ontology Classification of Filtered-DEGs
Gene Ontology (GO) enrichment of the filtered-DEGs provided insight into the function of target genes in plant developmental processes. Overall, 380 of 517 (67.7%) filtered-DEGs were annotated by at least one of the three categories of the GO function classification. Among these filtered-DEGs, 372 genes were related to cellular components, 356 genes were involved in molecular function, and 376 genes were involved in biological processes (Figure 4).
FIGURE 4. Gene Ontology classification of predicted-DEGs. A total of 380 DEGs were annotated by at least one of the three categories: biological process, cellular component, and molecular function. The x-axis indicates the level 2 GO terms, and the y-axis means the number of genes.
Within ‘cellular components,’ the most enriched filtered-DEGs were ‘cell,’ ‘cell part,’ ‘organelle,’ and ‘membrane.’ Within ‘molecular function,’ the most enriched items were ‘binding,’ ‘catalytic activity,’ ‘signal transducer activity,’ ‘structural molecule activity,’ and ‘transcription regulator activity’ terms. Within ‘biological process’, the most enriched items were ‘cellular process,’ ‘response to stimulus,’ ‘biological regulation,’ ‘metabolic process,’ and ‘regulation of biological process’ (Figure 4). Furthermore, the GO terms ‘reproduction,’ ‘reproductive process,’ ‘signaling,’ and ‘growth’ were also greatly enriched. These results suggest that signal transduction is essential for switchgrass flower reversion, and that most of the biological processes involved in switchgrass flower reversion are related to the cell membrane.
COG Classification of Predicted-DEGs
Using the Clusters of Orthologous Groups (COG) database of proteins, 214 genes of 517 (41.3%) filtered-DEGs were annotated and assigned to 21 functional clusters (Figure 5). The highest annotation frequency was ‘general function prediction only’ (97 genes, 45.3%), followed by ‘signal transduction mechanisms’ (94 genes, 43.9%), ‘posttranslational modification, protein turnover, and chaperones’ (65 genes, 30.4%), ‘secondary metabolites biosynthesis, transport and catabolism’ (65 genes, 30.4%), and ‘inorganic ion transport and metabolism’ (37 genes, 17.3%). We can deduce that signal transduction-related genes are essential for switchgrass flower reversion. Furthermore, there were 10 and 21 genes upregulated in functional clusters of ‘cell cycle control, cell division, and chromosome partitioning’ and ‘transcription,’ respectively.
FIGURE 5. Cluster of Orthologous Groups database of proteins (COG) functional classification of the predicted-DEGs. A total of 214 DEGs were annotated and assigned to 21 functional clusters. The x-axis indicates the COG categories as listed on the right, and y-axis means the number of genes.
Transcription Factors Implicated in Switchgrass Floral Reversion
Transcription factors (TFs) are a group of genes with different DNA-binding domains that influence the stress response and developmental processes in plants. We analyzed predicted-DEGs using a TF prediction tool2 (Jin et al., 2017), and 37 genes were annotated. Of these genes, 18 were assigned to the WRKY family, seven were assigned to the C2H2 family, seven were assigned to the NAC family, and each of the remaining five were assigned to bHLH, ERF, HSF, MYB, and TCP families. The WRKY TF family has the largest number of DEGs, which suggested that they play important roles in switchgrass flower reversion. Analysis showed that all 37 genes were upregulated (Supplementary Table S3).
KEGG Pathway Annotation of Predicted-DEGs
Pathway-based KEGG analysis of filtered-DEGs annotated 106 (20.5%) genes, of which 72 were enriched by 65 pathways. The most enriched pathway was ‘phenylpropanoid biosynthesis’ (12 genes, 2.3%), followed by ‘plant hormone signal transduction’ (9 genes, 1.7%), ‘plant–pathogen interaction’ (8 genes, 1.5%), and ‘MAPK signaling pathway - plant’ (6 genes, 1.1%) (Supplementary Table S4).
Of the 72 enriched filtered-DEGs, 15 genes were functionally related to plant hormone biosynthesis or signal transduction (Table 3). In the ‘zeatin biosynthesis’ and ‘plant hormone signal transduction’ pathways, six genes related to cytokinin biosynthesis or signal transduction were enriched (Table 3). We added 3 mg.L-1 benzylaminopurine (BAP) to the medium for switchgrass flower reversion induction. Genes involved in zeatin biosynthesis were annotated as cytokinin dehydrogenase (CKX, Pavir.5NG047700 and Pavir.J496300) and cis-zeatin O-glucosyltransferase (CISZOG, Pavir.7KG284900 and Pavir.7KG285000) respectively, which function in cytokinin degradation or cis-zeatin O-glucosylation, instead of as promoters of zeatin biosynthesis. Pavir.2KG175300 was annotated as a cytokinin receptor located at the membrane (Yamada et al., 2001) and Pavir.3NG010100 was annotated as a type-A Arabidopsis response regulator downstream of the cytokinin signal response (Hwang et al., 2012), and both were upregulated in reverse-developing explants. Three filtered-DEGs were annotated as auxin related genes (Table 3), of which, Pavir.2NG493200 was annotated as catalyzing indole-3-acetamide to indole-3-acetate (IAA) (Pollmann et al., 2003), and Pavir.3NG328700 and Pavir.4KG369800 were annotated as positive effectors of cell expansion downstream of the auxin signaling pathway (Spartz et al., 2012). In ethylene biosynthesis and signaling-related pathways, two genes were enriched (Table 3). Pavir.7NG301900 was annotated as catalyzing S-adenosyl-L-methionine to 1-aminocyclopropane-1-carboxylate, which is required for ethylene synthesis (Arteca and Arteca, 1999), and Pavir.9NG018300 was annotated as an ethylene response factor 1 (ERF1) TF involved in ethylene signal transduction. Additionally, two genes (Pavir.2NG586800 and Pavir.9KG267400) were annotated as Jasmonate ZIM-domain (JAZ) proteins and were involved in JA signaling (Yu et al., 2018), one gene was annotated as NPR1 and was involved in salicylic acid signaling, and one gene was annotated as SNRK2 and was involved in ABA signaling (Table 3, Fujii and Zhu, 2009).
TABLE 3. Kyoto Encyclopedia of Genes and Genomes (KEGG) enriched DEGs associated with hormone biosynthesis and signaling.
Gene Function Investigation of Filtered-DEGs
To verify the function of the filtered-DEGs in flower reversion, we selected four filtered-DEGs with high expression and large expression alteration refer to qPCR and RNA-seq results. The open reading frames were cloned into the pGreenII-based expression vector (Figure 6A and Supplementary Table S5; Hellens et al., 2000). Genes were transformed into the wild-type Colombia Arabidopsis by an Agrobacterium-mediated floral dip method (Clough and Bent, 1998). Gene transformation was verified by genotyping, gene expression was examined by semi-qPCR (Supplementary Presentation S1), and the phenotypes of screened offspring were further analyzed.
FIGURE 6. Investigation of the gene function of filtered-DEGs. (A) Expression vector employed. (B) Transgenic Pavir.8NG161900 plants compared with wild-type Arabidopsis. (C) Transgenic Pavir.7NG306800 and Pavir.1KG404800 plants compared with wild-type Arabidopsis. (D) Flowers of transgenic plants of Pavir.3NG053500 compared with wild-type Arabidopsis.
Ectopic expression of gene Pavir.8NG161900 (cysteine-rich RLK) caused delayed flowering time in Arabidopsis compared with wild-type plants (Figure 6B). Moreover, rosette and cauline leaves were enlarged, although the leaf number was unchanged (Figure 6B). Overexpression of Pavir.3NG053500 (dirigent protein) led to extremely shortened stamen filaments and undeveloped petals (Figure 6D). Furthermore, the transgenic plant was sterile and no seeds were produced when crossed and reciprocal crossed with wild-type plants. In transgenic Pavir.7NG306800 (F3H-2) and Pavir.1KG404800 (putative ripening-related protein 2) Arabidopsis plants, neither flowering structure nor timing was changed (Kim et al., 2008). However, the flower number was largely reduced (Figure 6C) and rosette leaves were enlarged (Figure 6C) compared to wild-type plants.
While flower reversion is an abnormal developmental process in plants, it allows the study of flower development from an interesting perspective. In this study, INFs from the switchgrass cultivar ‘Alamo’ incubated on SIM medium showed flower reversion, while those incubated on MS medium produced in vitro flowers. Anatomical and histological analyses showed that the emerged shoots originated from flower reversion and initiated from the pistil primordia in the floret primordium. Pairwise gene expression comparisons among REV, FLO, and INF generated 517 DEGs considered to regulate switchgrass flower reversion.
Cytokinin is an important plant hormone that regulates plant growth and development, and is considered to promote flower development. In Arabidopsis, exogenously applied BAP promoted flower formation in wild-type Columbia grown plants under 8 h short days (D’Aloia et al., 2011). In rice (Oryza sativa), LONELY GUY (LOG) encodes a cytokinin-activating enzyme and GRAIN NUMBER1 (GN1) encodes a cytokinin oxidase (CKX2), the log mutations caused small meristems and the gn1 mutations generated large inflorescence meristems and more grains (Yoshida and Nagato, 2011). In oilseed rape (Brassica napus), the content of all isoprenoid cytokinins increased significantly during vernalization, an essential process for inflorescence formation (Tarkowská et al., 2012). However, cytokinin can inhibit floral organ development. Overexpression of a cytokinin synthase isopentenyltransferase 4 (AtIPT4) in Arabidopsis increased cytokinin levels, and flowers of the transgene plants were abnormally developed (Li et al., 2010). AP1 is an important gene for the establishment of floral meristems, and functions as a suppressor of LOG and activator of the cytokinin degradation gene cytokinin oxidase 3 (CKX3). This suggests that plants will avoid cytokinin accumulation during flower meristem formation (Han et al., 2014). In the present study, exogenous application of BAP induced flower reversion in switchgrass. Gene expression analysis showed that two genes (Pavir.2KG175300 and Pavir.3NG010100) involved in cytokinin signal transduction were upregulated in REV. Four genes (Pavir.5NG047700, Pavir.J496300, Pavir.7KG284900, and Pavir.7KG285000) that degrade or inactivate cytokinin were also upregulated in REV, possibly due to feedback regulation from high concentrations of the exogenously applied cytokinin analog.
Additionally, the expression levels of genes related to biosynthesis and signal transduction of other hormones were also affected. Ethylene is an important hormone in plant reproductive development and fruit ripening, but may also play a role in inhibiting stamens and promoting carpel development in the early stage of melon (Cucumis melo) flower formation (Switzenberg et al., 2014). In switchgrass, genes associated with ethylene biosynthesis (Pavir.7NG301900) and signal transduction-related gene (Pavir.9NG018300) were upregulated in REV explants, suggesting that ethylene-related genes are involved in inhibiting flower development in switchgrass flower reversion. Besides, high levels of ethylene increased SAUR expression, a downstream auxin response gene, in Arabidopsis (Markakis et al., 2013). Two SAUR genes and an IAA biosynthesis-related gene were also upregulated in REV explants. JA is a widely studied hormone involved in the plant response to biotic and abiotic stressors, and reportedly contributes to rice spikelet morphogenesis by preventing OsMYC2 from activating OsMADS1, an E-class gene crucial for spikelet development (Cai et al., 2014). Two orthologs of JAZ in switchgrass were upregulated in REV explants, suggesting they function as flower development inhibitors in the reversion process. In addition, the expression of salicylic acid and ABA signaling-related genes was also altered. To conclude, plant hormone regulation networks were influenced by the addition of BAP in REV, and the reconstruction of plant hormone biosynthesis and signal transduction affected switchgrass flower reversion.
As mentioned above, our understanding of flower reversion is based on the analysis of relevant mutants (AGL24, LFY, and AG in A. thaliana; TM29 and SFT in S. lycopersicum; FBP2 in Petunia hybrid; IFA1, ID1, ZFL1, and ZFL2 in Z. mays; Laudencia-Chingcuanco and Hake, 2002; Molinero-Rosales et al., 2004). However, flower reversion is not simply a process of reversed flower development; genes that distinctively function in flower reversion were not detected in flower organ mutant studies. Comparison of flower reversion-related genes reported in previous research and the 517 DEGs obtained in this study revealed no overlap, further illustrating the difference between flower reversion and flower development.
Switchgrass flower reversion is a complex process. We examined functions of filtered-DEGs, four of which were overexpressed in Arabidopsis. Among them, Pavir.8NG161900 and Pavir.3NG053500 inhibited reproductive development, while Pavir.7NG306800 and Pavir.1KG404800 promoted vegetative development. No plant was observed with complete flower reversion in any transgene event; this might be attributed to the complexity of flower reversion in switchgrass, which is regulated by multiple genes (517 possible related genes) and the overexpression of a single gene is insufficient to induce the complete process of flower reversion. Nevertheless, flower reversion involves multiple biological processes, including the promotion of vegetable growth, suppression of flower formation, and dedifferentiation of flower organs. Although overexpression of the filtered-DEGs did not show complete flower reversion, vegetative and reproductive growth were affected, both of which are important in the flower reversion process.
Using the Illumina platform, we analyzed the gene expression of immature switchgrass inflorescences incubated on SIM and MS media, and uncultured INFs. Comparison of gene expression patterns identified 517 genes involved in switchgrass floral reversion. Based on the annotations of the NR, Swiss-Prot, GO, COG, and KEGG databases, we found that signal transduction and metabolism are essential for switchgrass flower reversion. Switchgrass flower reversion is a complex process that involves multiple biological processes, including the promotion of vegetative growth, suppression of flower formation, and dedifferentiation of flower organs; this was supported by the overexpression phenotypes of some of the filtered-DEGs in Arabidopsis.
Materials and Methods
Plant Materials and Explant Culture
The lowland switchgrass cultivar ‘Alamo’ was used in this study. Plants were grown in an outdoor field located at the Northwest A&F University in Yangling, Shaanxi, China (108.072°E, 34.295°N). Shoot apexes were harvested at the E4 stage in early July (Moore et al., 1991). The shoot apexes were surface sterilized with 70% ethanol and 8% sodium hypochlorite. After rinsing three times with sterile distilled water, shoot apexes were then cut by 0.5 cm at both ends and split longitudinally. Both halves of the shoot apex were placed, section downward, on solid MS (Murashige and Skoog, 1962) or SIM (Chai et al., 2012) medium and incubated under constant temperature (25°C) and a 20 h light/4 h dark photoperiod. The MS medium consisted of MS salts, organics, vitamins, 7.5 g/L-1 agar, and 30 g/L-1 sucrose. The SIM medium was MS medium with 3 mg/L-1 benzylaminopurine (BAP, Sanland, Fujian, China). More than 90% of the explants on the SIM medium induced neonatal shoots from flower reversion after 60 days of incubation.
Microscopic Structure Analysis of Switchgrass Floral Reversion
Shoot clumps produced on SIM medium were harvested and dismembered under a stereoscopic microscope (Nikon SMZ1500, Japan), and their morphological and anatomical characteristics were identified and described. Ten explants were analyzed per treatment.
Separated shoots were fixed with FAA (formalin: acetic acid: 70% alcohol = 5:5:90) stationary liquid for 6 h and stained with hematoxylin (Maikun, Shanghai, China) for 3 days, followed by color differentiation in tap water for 1 h (Avwioro, 2011). Stained materials were dehydrated with a gradient alcohol series (20 min each of 70, 80, and 90%, followed by 100% twice for 15 min each), cleared with xylene, and infiltrated with melted paraffin wax. Then, they were embedded in a paraffin block and cut into 10 μm slices using a paraffin microtome (JinHuaHuiYou HY-202A, China). Slices were dewaxed with xylene, post-mounted with Permount Mounting Medium (HuShi, Shanghai, China), and observed under an optical microscope (Chongqing UOP, UB203i, China).
RNA Isolation and mRNA-seq
The explants incubated for 20 days on MS and SIM media were flash frozen in liquid nitrogen together with uncultured INFs. Four replications were included in each treatment, and five explants were employed for each replicate.
Total RNA was isolated using TRIzol reagent (Invitrogen, Carlsbad, CA, United States) according to the manufacturer’s instructions. RNA quality and quantity were determined using electrophoresis on a 1.5% agarose gel using an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, United States) and a Nanodrop 2000 Spectrophotometer (Thermo Fisher Scientific, Waltham, MA, United States), respectively. mRNA was purified using a NEBNext Poly(A) mRNA Magnetic Isolation Module (E7490; New England Biolabs, Ipswich, MA, United States). The sequencing library for each sample was constructed using the NEBNext mRNA Library Prep Master Mix Set for Illumina (E6110; New England Biolabs) and the NEBNext Multiplex Oligos for Illumina (E7500; New England Biolabs). The insertion length was set to 180 bp and verified by electrophoresis on a 1.8% agarose gel. After quantifying the sequencing libraries with Library Quantification Kit-Illumina GA Universal (KK4824; Kapa Biosystems, Inc., Wilmington, MA, United States), the libraries were used for cluster station generation on Illumina cBot (Illumina, San Diego, CA, United States). Paired-end sequencing was performed on an Illumina HiSeq 2500 (Illumina).
Gene Expression Analysis and Mining of DEGs
Paired-end reads with a length of 202 bp were assessed with Q20, Q30, and GC content. Reads containing >20% low quality bases (quality score < 20) or >5% N were filtered out before assembly. All filtered transcriptome datasets were mapped versus the switchgrass genomic assembly3 (version 4.1; Casler et al., 2011) with Hisat2 (Kim et al., 2015) software and assembled by the Stringtie software without the ‘-e’ option to predict novel transcripts (Pertea et al., 2015, 2016). The assembled transcripts were merged together under the direction of GFF annotation (version 4.1) using Stringtie’s merge function to create a uniform set of transcripts for all 12 samples. The mapped reads were assembled again using Stringtie software with the uniform set of transcripts as a direction, and the ‘-e’ and ‘-B’ options were used to restrict novel transcript prediction and generate input file for DEG analysis (Pertea et al., 2016). The fold change was calculated using the R package ‘Ballgown’ (Fu et al., 2018), and DEGs were filtered with p-value ≤ 0.01 and fold change ≥ 2 as restrictions.
Annotation of DEGs
The software program BLAST (Altschul et al., 1997) was used to provide the most inclusive functional description of the assembled sequences with a threshold E-value ≤ 1e-5. Annotation information was based on four sources: (1) the NR protein database (Pruitt et al., 2007) from the NCBI, (2) Swiss-Prot (Apweiler et al., 2004) from the Universal Protein Resource (UniProt), (3) the high-level functions and biological system pathway database KEGG (Kanehisa et al., 2004), and (4) the Cluster of Orthologous Groups (COG) database of proteins (Tatusov et al., 2000). We also obtained the GO annotation of genes based on the NR annotation using the Blast2GO program with E-value ≤ 1e-6 (Conesa et al., 2005) and GO functional classification for all genes using WEGO software (Ye et al., 2006).
Real-Time qPCR Analysis
Total RNAs of reversed, unreversed, and uncultured INFs (three replications for each treatment) were isolated as described above. Using the PrimeScript RT reagent kit (Takara, Kyoto, Japan), the first strand of cDNA was synthesized according to the attached protocol. Real-time qPCR was performed using SYBR Premix EX Taq II Kit (Takara) on a QuantStudio 3 Flex Real Time System (Thermo Fisher Scientific). The total reaction volume was 20 μl, and the primers used are listed in Supplementary Table S6. The PCR reaction was performed as follows: 10 min at 95°C and 40 cycles of 94°C for 20 s and 60°C for 30 s. Data generated was analyzed using the ΔΔCt method (Biassoni and Raso, 2014). Switchgrass gene ef1-α was used as an internal control (Gimeno et al., 2014).
Gene Transformation of Arabidopsis
A binary vector based on pGreenII with a double 35S promoter for gene expression and a bar gene for selection was used (Hellens et al., 2000). A recombinant vector was introduced into the Agrobacterium tumefaciens strain GV3101 harboring the help plasmid pSoup and used for genetic transformation. A. tumefaciens was cultured at 28°C in liquid YEB medium with shaking at 250 rpm, and cells were harvested by centrifugation (5,000 rpm for 15 min at room temperature) and resuspended in infiltration medium to OD600 of approximately 0.8 (Clough and Bent, 1998). The gene transformation of Arabidopsis plants was performed according to the floral dip method (Clough and Bent, 1998) and seeds were harvested approximately 3–4 weeks later. After drying for 7 days at room temperature, seeds were sowed in well-watered soil and keep at 4°C for 2 days together with the plates, and the culture was transferred to normal conditions (Sanchez-Serrano and Salinas, 2014). After 7 days, seedlings with fully extended cotyledons were sprayed with 0.05% phosphinothricin (Sigma-Aldrich, St. Louis, MO, United States), and the surviving plants were further examined by PCR.
The raw sequence reads are available at the NIH Short Read Archive database (accession number SRP144197).
XY and WY designed the experiments. WY, ZA, SF, LM, and XK conducted the experiments. ZC helped to conduct experiments. WY and LS wrote the manuscript. All authors read and approved the manuscript.
This work was supported by the National Natural Science Foundation of China (31171607) and National Natural Science Foundation Committee.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2018.01805/full#supplementary-material
ABA, abscisic acid; COG, Cluster of Orthologous Groups database of proteins; DEGs, differentially expressed genes; FLO, immature inflorescence produced floret on MS medium; GO, Gene Ontology; INF, immature inflorescence; JA, jasmonic acid; KEGG, Kyoto Encyclopedia of Genes and Genomes; MS, Murashige-Skoog; NCBI, National Center for Biotechnology Information; NR, non-redundant; REV, immature inflorescence produced flower reversion on SIM; RNA-seq, RNA sequencing; SIM, shoot induction medium; TF, transcription factor.
- ^ https://genome.jgi.doe.gov/portal/pages/dynamicOrganismDownload.jsf?organism=Pvirgatum
- ^ http://planttfdb.cbi.pku.edu.cn/prediction.php
- ^ https://genome.jgi.doe.gov/portal/pages/dynamicOrganismDownload.jsf?organism=Pvirgatum
Altschul, S. F., Madden, T. L., Schaffer, A. A., Zhang, J. H., Zhang, Z., Miller, W., et al. (1997). Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 25, 3389–3402. doi: 10.1093/nar/25.17.3389
Ampomah-Dwamena, C., Morris, B. A., Sutherland, P., Veit, B., and Yao, J. L. (2002). Down-regulation of TM29, a tomato SEPALLATA homolog, causes parthenocarpic fruit development and floral reversion. Plant Physiol. 130, 605–617. doi: 10.1104/Pp.005223
Arteca, J. M., and Arteca, R. N. (1999). A multi-responsive gene encoding 1-aminocyclopropane-1-carboxylate synthase (ACS6) in mature Arabidopsis leaves. Plant Mol. Biol. 39, 209–219. doi: 10.1023/A:1006177902093
Burris, J. N., Mann, D. G., Joyce, B. L., and Stewart, C. N. (2009). An improved tissue culture system for embryogenic callus production and plant regeneration in switchgrass (Panicum virgatum L.). Bioenergy Res. 2, 267–274. doi: 10.1007/s12155-009-9048-8
Casler, M. D., Tobias, C. M., Kaeppler, S. M., Buell, C. R., Wang, Z. Y., Cao, P., et al. (2011). The switchgrass genome: tools and strategies. Plant Genome 4:273. doi: 10.3835/plantgenome2011.10.0026
Chai, G., Xu, K., Wang, Y., Sun, F., Song, L., Liu, S., et al. (2012). Establishment of a highly efficient regeneration system of switchgrass artificial spike buds. Acta Prataculturae Sinica 4, 98–104.
Colasanti, J., Yuan, Z., and Sundaresan, V. (1998). The indeterminate gene encodes a zinc finger protein and regulates a leaf-generated signal required for the transition to flowering in maize. Cell 93, 593–603. doi: 10.1016/S0092-8674(00)81188-5
Conesa, A., Gotz, S., Garcia-Gomez, J. M., Terol, J., Talon, M., and Robles, M. (2005). Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 21, 3674–3676. doi: 10.1093/bioinformatics/bti610
D’Aloia, M., Bonhomme, D., Bouché, F., Tamseddak, K., Ormenese, S., Torti, S., et al. (2011). Cytokinin promotes flowering of Arabidopsis via transcriptional activation of the FT paralogue TSF. Plant J. 65, 972–979. doi: 10.1111/j.1365-313X.2011.04482.x
Fujii, H., and Zhu, J. K. (2009). Arabidopsis mutant deficient in 3 abscisic acid-activated protein kinases reveals critical roles in growth, reproduction, and stress. Proc. Natl. Acad. Sci. U.S.A. 106, 8380–8385. doi: 10.1073/pnas.0903144106
Gimeno, J., Eattock, N., Deynze, A. V., and Blumwald, E. (2014). Selection and validation of reference genes for gene expression analysis in switchgrass (Panicum virgatum) using quantitative real-time RT-PCR. PLoS One 9:e91474. doi: 10.1371/journal.pone.0091474
Han, Y., Zhang, C., Yang, H., and Jiao, Y. (2014). Cytokinin pathway mediates APETALA1 function in the establishment of determinate floral meristems in Arabidopsis. Proc. Natl. Acad. Sci. U.S.A. 111, 6840–6845. doi: 10.1073/pnas.1318532111
Hellens, R. P., Edwards, E. A., Leyland, N. R., Bean, S., and Mullineaux, P. M. (2000). pGreen: a versatile and flexible binary Ti vector for Agrobacterium-mediated plant transformation. Plant Mol. Biol. 42, 819–832. doi: 10.1023/A:1006496308160
Jiang, L., Yan, S. S., Yang, W. C., Li, Y. Q., Xia, M. X., Chen, Z. J., et al. (2015). Transcriptomic analysis reveals the roles of microtubule-related genes and transcription factors in fruit length regulation in cucumber (Cucumis sativus L.). Sci. Rep. 5:8031. doi: 10.1038/Srep08031
Jin, J. P., Tian, F., Yang, D. C., Meng, Y. Q., Kong, L., Luo, J. C., et al. (2017). PlantTFDB 4.0: toward a central hub for transcription factors and regulatory interactions in plants. Nucleic Acids Res. 45, D1040–D1045. doi: 10.1093/nar/gkw982
Li, X. G., Su, Y. H., Zhao, X. Y., Li, W., Gao, X. Q., and Zhang, X. S. (2010). Cytokinin overproduction-caused alteration of flower development is partially mediated by CUC2 and CUC3 in Arabidopsis. Gene 450, 109–120. doi: 10.1016/j.gene.2009.11.003
Liu, B., Jiang, G. F., Zhang, Y. F., Li, J. L., Li, X. J., Yue, J. S., et al. (2011). Analysis of transcriptome differences between resistant and susceptible strains of the citrus red mite Panonychus citri (Acari: tetranychidae). PLoS One 6:e28516. doi: 10.1371/journal.pone.0028516
Markakis, M. N., Boron, A. K., Van, L. B., Saini, K., Cirera, S., Verbelen, J. P., et al. (2013). Characterization of a small auxin-up RNA (SAUR)-like gene involved in Arabidopsis thaliana development. PLoS One 8:e82596. doi: 10.1371/journal.pone.0082596
Molinero-Rosales, N., Latorre, A., Jamilena, M., and Lozano, R. (2004). Single flower truss regulates the transition and maintenance of flowering in tomato. Planta 218, 427–434. doi: 10.1007/s00425-003-1109-1
Moore, K. J., Moser, L. E., Vogel, K. P., Waller, S. S., Johnson, B. E., and Pedersen, J. F. (1991). Describing and quantifying growth stages of perennial forage grasses. Agron. J. 83, 1073–1077. doi: 10.2134/agronj1991.00021962008300060027x
Okamuro, J. K., denBoer, B. G. W., Lotys-Prass, C., Szeto, W., and Jofuku, K. D. (1996). Flowers into shoots: photo and hormonal control of a meristem identity switch in Arabidopsis. Proc. Natl. Acad. Sci. U.S.A. 93, 13831–13836. doi: 10.1073/pnas.93.24.13831
Palmer, N. A., Saathoff, A. J., Kim, J., Benson, A., Tobias, C. M., Twigg, P., et al. (2012). Next-generation sequencing of crown and rhizome transcriptome from an upland, tetraploid switchgrass. Bioenerg Res. 5, 649–661. doi: 10.1007/s12155-011-9171-1
Parrish, D. J., Casler, M. D., and Monti, A. (2012). “The evolution of switchgrass as an energy crop,” in Switchgrass, A Valuable Biomass Crop for Energy, ed. A. Monti (London: Springer-Verlag), 1–28.
Pertea, M., Kim, D., Pertea, G. M., Leek, J. T., and Salzberg, S. L. (2016). Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. Nat. Protoc. 11, 1650–1667. doi: 10.1038/nprot.2016.095
Pertea, M., Pertea, G. M., Antonescu, C. M., Chang, T. C., Mendell, J. T., and Salzberg, S. L. (2015). StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat. Biotechnol. 33, 290–295. doi: 10.1038/nbt.3122
Pollmann, S., Neu, D., and Weiler, E. W. (2003). Molecular cloning and characterization of an amidase from Arabidopsis thaliana capable of converting indole-3-acetamide into the plant growth hormone, indole-3-acetic acid. Phytochemistry 62, 293–300. doi: 10.1016/S0031-9422(02)00563-0
Pruitt, K. D., Tatusova, T., and Maglott, D. R. (2007). NCBI reference sequences (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res. 35, D61–D65. doi: 10.1093/nar/gkl842
Sen, S., Krishnasamy, V., Krishnaveni, S., Vijayalakshmi, C., Rajamani, K., and Gnanam, R. (2013). Evidence of floral reversion during in vitro shoot regeneration from young inflorescence of the medicinal geophyte Gloriosa superba L. Plant Tissue Cult. Biotechnol. 23, 115–120. doi: 10.3329/ptcb.v23i1.15568
Spartz, A. K., Lee, S. H., Wenger, J. P., Gonzalez, N., Itoh, H., Inzé, D., et al. (2012). The SAUR19 subfamily of small auxin up RNA genes promote cell expansion. Plant J. 70, 978–990. doi: 10.1111/j.1365-313X.2012.04946.x
Switzenberg, J. A., Little, H. A., Hammar, S. A., and Grumet, R. (2014). Floral primordia-targeted ACS (1-aminocyclopropane-1-carboxylate synthase) expression in transgenic Cucumis melo implicates fine tuning of ethylene production mediating unisexual flower development. Planta 240, 797–808. doi: 10.1007/s00425-014-2118-y
Tarkowská, D., Filek, M., Biesaga-Kościelniak, J., Marcińska, I., Macháèková, I., Krekule, J., et al. (2012). Cytokinins in shoot apices of Brassica napus plants during vernalization. Plant Sci. 187, 105–112. doi: 10.1016/j.plantsci.2012.02.003
Tatusov, R. L., Galperin, M. Y., Natale, D. A., and Koonin, E. V. (2000). The COG database: a tool for genome-scale analysis of protein functions and evolution. Nucleic Acids Res. 28, 33–36. doi: 10.1093/Nar/28.1.33
Trevaskis, B., Tadege, M., Hemming, M. N., Peacock, W. J., Dennis, E. S., and Sheldon, C. (2007). Short vegetative phase-like MADS-box genes inhibit floral meristem identity in barley. Plant Physiol. 143, 225–235. doi: 10.1104/pp.106.090860
Wagner, G. P., Kin, K., and Lynch, V. J. (2012). Measurement of mRNA abundance using RNA-seq data: RPKM measure is inconsistent among samples. Theory Biosci. 131, 281–285. doi: 10.1007/s12064-012-0162-3
Wang, X. J., Xu, R. H., Wang, R. L., and Liu, A. Z. (2012). Transcriptome analysis of Sacha inchi (Plukenetia volubilis L.) seeds at two developmental stages. BMC Genomics 13:716. doi: 10.1186/1471-2164-13-716
Yamada, H., Suzuki, T., Terada, K., Takei, K., Ishikawa, K., Miwa, K., et al. (2001). The Arabidopsis AHK4 histidine kinase is a cytokinin-binding receptor that transduces cytokinin signals across the membrane. Plant Cell Physiol. 42, 1017–1023. doi: 10.1093/pcp/pce127
Yang, Y., Guo, J., Yan, P., Li, Y., Liu, K., Gao, P., et al. (2014). Transcriptome profiling identified multiple jasmonate ZIM-domain proteins involved in the regulation of alkaloid biosynthesis in tobacco BY-2 Cells. Plant Mol. Biol. Rep. 33, 153–166. doi: 10.1007/s11105-014-0740-2
Yu, X., Chen, G., Tang, B., Zhang, J., Zhou, S., and Hu, Z. (2018). The Jasmonate ZIM-domain protein gene SlJAZ2 regulates plant morphology and accelerates flower initiation in Solanum lycopersicum plants. Plant Sci. 267, 65–73. doi: 10.1016/j.plantsci.2017.11.008
Zhang, L., Wang, L., Yang, Y., Cui, J., Chang, F., Wang, Y., et al. (2014). Analysis of Arabidopsis floral transcriptome: detection of new florally expressed genes and expansion of Brassicaceae-specific gene families. Front. Plant Sci. 5:802. doi: 10.3389/fpls.2014.00802
Keywords: switchgrass, floral reversion, cytokinin, flower maintenance, transcriptome
Citation: Yongfeng W, Aiquan Z, Fengli S, Mao L, Kaijie X, Chao Z, Shudong L and Yajun X (2018) Using Transcriptome Analysis to Identify Genes Involved in Switchgrass Flower Reversion. Front. Plant Sci. 9:1805. doi: 10.3389/fpls.2018.01805
Received: 02 August 2018; Accepted: 20 November 2018;
Published: 04 December 2018.
Edited by:Nunzio D’Agostino, The Council for Agricultural Research and Economics, Italy
Reviewed by:Hamed Bostan, North Carolina State University, United States
Valentino Ruggieri, Centre for Research in Agricultural Genomics (CRAG), Spain
Copyright © 2018 Yongfeng, Aiquan, Fengli, Mao, Kaijie, Chao, Shudong and Yajun. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xi Yajun, email@example.com